The climate history of the Greater Yellowstone Area shows changes on timescales ranging from seasons to millennia. Over thousands of years, the primary drivers of natural climate change are cyclical variations in solar radiation related to Earth’s orbit around the sun and associated changes in the amount of greenhouse gas in the atmosphere. Over years to centuries, the natural drivers of climate variability are volcanic activity, solar output, and coupled atmosphere-ocean circulation patterns. [high agreement, robust evidence]
The geologic record of the GYA indicates that the last glaciation (approximately 22,000-13,000 yr ago) was as much as 5-7°F (2.8-3.9°C) colder than the pre-industrial period (1850-1900). Two warm periods in the past are the early Holocene (11,500-7000 yr ago), which was about 1.8-3.6°F (1-2°C) warmer in summer than the pre-industrial period, and the Medieval Climate Anomaly (from 800 to 1300), characterized by prolonged droughts and slightly warmer summers than pre-industrial time. [high agreement, robust evidence]
The average temperature of the last two decades (2001-2020) is probably as high or higher than any period in the last 20,000 yr, and likely higher than previous glacial and interglacial periods in the last 800,000 yr. The current level of carbon dioxide in the atmosphere is the highest in the last 3.3 million years. [medium agreement, medium evidence]
What is Climate?
Climate differs from weather. Weather refers to atmospheric changes that occur over minutes to months and are reflected, for example, by the temperature, humidity, and precipitation at a location and particular time. Climate is the long-term average of weather over an extended time period, such as decades to centuries. In this report, we define the climate base period for comparison with future periods as the 1986 through 2005 average or mean. We chose this 20-year base period because 1) it captures observed global warming trends and, therefore, is a conservative (warm) baseline; and 2) climate model simulations of the historical period end in 2005 and projections of future climate in this Assessment begin in 2006.
The climate system describes all the interacting components that create Earth’s climate: the atmosphere (air), hydrosphere (water), the cryosphere (ice and permafrost), lithosphere (Earth’s upper rocky layer), and biosphere (living things). Climate change refers to shifts (e.g., decadal and longer) in the average or mean climate, which can be abrupt or gradual, as evidenced in historical and geological records discussed later in this chapter and in Chapter 3. A climate trend is a long-term trajectory of change in the mean climate. Climate variability refers to short-term departures from the mean state of the climate (note that climate variations are longer than individual weather events, spanning seasons or years). In the coming decades, climate change is projected to trend toward ever warmer conditions; however, as illustrated in Figure 2-1, climate variability may result in seasons and years that are warmer or colder than the 20-year means, just as occurs today.
Climate and Water Variables Discussed in the Assessment
The international climate science community uses over 50 essential physical, chemical, and biological variables to characterize the state of the Earth’s climate (WMOa undated). To qualify as an essential climate variable, the information about it must be 1) worldwide in coverage; 2) freely available; 3) quality controlled with appropriate documentation; and 4) considered relevant by an international panel of climate experts. Our report focuses on a small subset of the 50 essential climate variables that are relevant to the GYA:
Air temperature (referred to temperature in this report) is a measure of how hot or cold the air is with reference to some standard value. Seasonal variations in temperature result from latitudinal differences in the amount of solar radiation received at the Earth’s surface, contrasts in seasonal heating of land and oceans, and atmospheric circulation.
Precipitation is the quantity of water (liquid or solid) falling to the Earth’s surface at a specific place over a given period. Like temperature, precipitation varies from season to season and place to place and depends on coupled atmospheric-ocean circulation.
In addition to the climate variables, we also focus on other variables:
Snowfall and snowpack are measures of the amount and fate of solid winter precipitation. Snowfall is the amount of accumulated snow after a storm. It is measured in terms of the depth of solid water it contains. In mountainous and relatively dry areas like the GYA, 10 inches (25 cm) or more of snow is often needed to create 1 inch (2.5 cm) of liquid water when melted. Snowpack is the amount of snowfall that accumulates over the cold season. It also is measured by both depth (snow depth) and the amount of liquid water it stores (called snow water equivalent or SWE).
Streamflow (also called discharge) refers to water moving within a river measured by the volume of water passing a point in a given time. Streamflow is measured at gaging stations in units of cubic feet per second or cubic meters per second. In GYA, streamflow is strongly controlled by the seasonality of runoff from snowmelt.
Runoff is the depth of water uniformly distributed over an area, such as a watershed. It is the potential amount of water available for groundwater and streamflow.
Evapotranspiration is water lost through evaporation from bare soil and transpiration by plants. Potential evapotranspiration is the amount of evapotranspiration that would occur under unlimited water availability.
Drought is a prolonged period of dryness relative to long-term average conditions. The climatological community defines four types of drought: 1) meteorological drought occurs when unusually dry weather patterns persist over an area from days to months; 2) hydrological drought refers to low-water supply and usually occurs after many months of meteorological drought; 3) agricultural drought occurs when low soil moisture limits survival and production of crops and grazing lands; and 4) socioeconomic drought reflects the economic and social impact of a combination of hydrological and agricultural drought. In this report, we use the term drought, without distinguishing the type, but unless otherwise noted, we are referring to meteorological or hydrological drought.
Palmer Drought Severity Index (PDSI) is a standard measure of drought that combines temperature or potential evapotranspiration and precipitation data to quantify dryness or wetness relative to average or normal conditions. The PDSI describes soil moisture conditions (generally the top meter of soil).
Vapor pressure deficit is a measure of the drying capacity of the atmosphere based on air temperature and relative humidity. High vapor pressure deficits (i.e., high temperature combined with low humidity) can limit tree growth, increase their vulnerability to drought, and dry fuels, all potential contributors to wildfire.
The climate of the GYA is characterized by long, often bitterly cold winters. Summers are short and mild. May and June are generally the wettest months in the valleys; August is generally the driest. Snow is the primary form of winter precipitation.
The GYA’s climate is attributed primarily to its mid-latitude continental location, high average elevation, and distance from the Pacific and Gulf coasts. At approximately 44°N latitude, the region has long summer days and long winter nights. Even summer days are relatively cool, however, due to the high elevation of the GYA. GYA receives air masses not only from the Pacific Ocean to the west, but also from the Arctic Ocean to the north and Gulf of Mexico to the south. The relative contribution of these air masses and the moisture they entrain is reflected in seasonal temperature and precipitation patterns for any given year (Whitlock and Bartlein 1993).
Precipitation generally increases with elevation in GYA, as it does throughout the West. Cold, wet winters in the GYA reflect a combination of moisture carried by storms off the Pacific Ocean and frequent, cold Arctic air mass intrusions. Most of these storms are funneled northeastward along the Snake River Plain and the precipitation they carry is delivered as snow over the high mountains and plateaus of the GYA (Farnes 1997). Cold, dry weather in winter occurs when a sustained southward incursion of an Arctic air mass brings subzero temperatures (Fahrenheit) to the region. Winters are generally wetter in the Teton Range and western Yellowstone Plateau region than in the eastern GYA. Pacific storm systems, as well as moisture transported along the Rocky Mountain Front from the Gulf of Mexico, account for wet spring conditions in the region.
Summers in much of the GYA are typified by warm, dry conditions punctuated by thunderstorms. During summer, Pacific storm tracks shift well north of the GYA so summer rainfall is delivered by low-pressure centers and their related atmospheric disturbances (or fronts). Moisture in the northern and eastern GYA originates from the subtropical Gulf of Mexico, whereas that in the southwestern GYA comes from the subtropical Pacific Ocean.
Year-to-year and decadal climate variations that affect the GYA derive from recurring, global‑scale changes in atmosphere and ocean circulation patterns. The El Niño-Southern Oscillation (ENSO), for example, is a climate pattern set up by changes in sea-surface temperature and atmospheric pressure in the equatorial Pacific Ocean that can persist for several years. Warmer-than-normal equatorial ocean surface temperatures are associated with El Niño events, whereas colder surface temperatures are associated with La Niña events. ENSO influences storm tracks and pressure systems at mid-latitudes through atmospheric connections (called teleconnections) that, in turn, influence surface climate conditions across the West, including the GYA (Figure 2-2). The Pacific Decadal Oscillation (PDO) is a similar, multi-year pattern of climate variability forced by sea-surface temperature changes that occur on decadal scales. Phases of the PDO are identified by warm or cold ocean temperature patterns in the north Pacific Ocean. Even persistent decades of warmer and colder than normal sea-surface temperature in the North Atlantic Ocean known as the Atlantic Multi-decadal Oscillation (AMO) can interact with ENSO and PDO to affect long-term drought in the GYA (McCabe et al. 2004, 2008).
ENSO and PDO patterns alter the north-south position of Pacific storm tracks across western North America, which can result in large and contrasting variations in winter precipitation and air temperature that persist for short (~12-18 months) to long (decades) periods. The regional effects on precipitation from changes in the PDO are strongest along the Pacific Coast in the Pacific Northwest, whereas the greatest influence of ENSO is over the American Southwest and Southeast. On average, El Niño events bring warmer and drier conditions to the GYA whereas La Niña events bring cooler and wetter conditions; however, interaction with other atmosphere-ocean circulation processes often affect this generalized pattern. Thus, not all ENSO or PDO events have a similar effect on the climate of the GYA (Pederson et al. 2011a,b; Abatzoglou 2011; Pederson et al. 2013).
Past Climate Change
Natural climate change, ever ongoing, can be examined on many timescales. Over millions to 100s of millions of years, changes in the size and position of continents and ocean basins and related mountain uplift have shaped the Earth’s climate. Over tens to hundreds of thousands of years, repeated cycles of cold (called glacial periods or ice ages) and warmth (called interglacial periods) have been caused by seasonal and latitudinal variations in the amount of solar radiation received by the Earth.
Glacial-interglacial cycles result from continual changes in the tilt and wobble of Earth’s axis and in the elliptical orbit of the Earth around the sun. These recurring astronomical drivers of climate change, known as Milankovitch cycles, are the pacemaker of the ice ages (Ruddiman 2013). Levels of greenhouse gases in the atmosphere also varied with these cycles, amplifying the warmth of interglacial periods and the cold conditions of glacial intervals. The last ice age on the planet occurred between about 115,000 and 11,700 yr ago, with maximum glaciation between 27,000 and 19,000 yr ago in different regions (Clark et al. 2009). Global warming of 5-7°F (3-4°C) occurred between 19,000-11,000 yr ago, ushering in the current interglacial period, which is called the Holocene (the last 11,700 yr) (Clark et al. 2012; IPCC 2013).
Climate variations on timescales of centuries or less tend to be more regional in scale, but with different principal drivers (Ruddiman 2013):
Over decades to centuries, volcanic activity, changes in solar output, and global-scale changes in atmosphere-ocean circulation patterns have caused climate to vary. Notable examples of such variations have resulted in climate anomalies relative to a defined average, such as the persistent cold conditions and widespread glacial advances that occurred from about 1600-1850, known as the Little Ice Age, and the periods of warmth and drought that define the Medieval Climate Anomaly from about 800-1300.
Over interannual to decadal timescales, persistent atmosphere-ocean circulation patterns, such as ENSO and the PDO, are the important drivers of climate variations (discussed above).
The last 20,000 years
The climate of the GYA has varied widely over the last 20,000 yr—from the culmination of an ice age to periods that were warmer than the pre-industrial period. The climate history of the GYA, as interpreted from the geologic record and measured by observations, provides a useful context for perspective on the significance of current and projected climate changes.
GYA was extensively covered by ice during past glacial periods. Ice cover during the recent Pinedale glaciation (22,000-13,000 yr ago) and the previous Bull Lake glaciation (150,000-140,000 yr ago) are shown in Figure 2-3 (Licciardi and Pierce 2018). The Pinedale glaciation began when glaciers started to grow and expand in the Beartooth-Absaroka Mountains of northeastern GYA and Gallatin Range of northwestern GYA. By 15,000 yr ago, individual glaciers from the two regions had coalesced into a large Yellowstone ice cap centered over present-day Yellowstone Lake. Valley glaciers flowed from the ice cap down all the major river valleys. Geologists name the terminal ridges of gravel and boulders (moraines) deposited by these valley glaciers by their location (e.g., the Chico moraine, the Outer Jenny Lake moraine) and determine the age of the moraines using cosmogenic nuclide dating methods. From this information, geologists have determined that Yellowstone ice cap was asymmetrical; its maximum growth occurred to the southwest, indicative of the dominant source of precipitation from the direction of the Snake River Plain (Licciardi and Pierce 2018). Glaciers started to recede first in the northeast in the Clarks Fork drainage 19,800 yr ago, and last in the south at Jackson Lake 15,500 yr ago. Glacial ice was largely gone from the GYA by 12,000 yr ago.
The present-day landscape, river systems, and lakes of the GYA were formed largely during the Pinedale glaciation by the erosional and depositional processes associated with ice advance, melting, and recession (Good and Pierce 1996). The ruggedness of Teton and Wind River ranges exemplify glacial sculpting under former ice divides (Figure 2-4). The sagebrush-covered terraces within the major river valleys were created by high-volume braided rivers that flowed from melting glaciers and deposited coarse gravels beyond the ice margins. These porous gravels are the source of shallow, groundwater storage in most of the GYA’s river basins. Many lakes in the GYA (e.g., Jenny Lake, Jackson Lake, Fremont Lake) are dammed by moraines of gravel that were deposited at the terminus of valley glaciers. Other smaller lakes (e.g., Blacktail Pond, Swan Lake, Swamp Lake) were formed when blocks of ice buried under glacial debris melted with warming temperatures and created a depression on the land surface.
As the climate warmed and glaciers started to melt in the GYA, plants were able to colonize areas that had previously been covered by ice. Pollen buried in the sediment in Yellowstone’s lakes indicates that the first conifer to appear was juniper, probably common juniper (Juniperus communis), which established in a relatively open, tundra-like landscape. Next came Engelmann spruce (Picea engelmannii), followed by whitebark pine (Pinus albicaulis), limber pine (Pinus flexilis), and subalpine fir (Abies lasiocarpa) (Krause and Whitlock 2017). Lodgepole pine (Pinus contorta) was widespread after 11,000 yr ago, and Douglas-fir (Pseudotsuga menziesii) was the last conifer to arrive and expand its range after 9000 yr ago (Iglesias et al. 2018).
This sequence of forest development shows the capacity of the region’s conifers to respond to rising temperatures by adjusting their range and abundance over thousands of years. Similar responses will certainly take place in the future, but likely at a faster rate. Some native species (e.g., whitebark pine) may no longer find suitable climate in GYA and become regionally absent (Chang et al. 2014) and different species (e.g., Gambels oak [Quercus gambelii], western larch [Larix occidentalis], ponderosa pine [Pinus ponderosa]) may be better suited to future climate conditions. The rate of current climate change, however, is many times faster than what occurred in the past, and it is doubtful that species will be able to keep pace on a timescale relevant to forest management (Bartlein et al. 1997).
The current interglacial period, the Holocene, began as the latest of a series of interglacial periods. Two warm intervals in the Holocene serve as important benchmarks for evaluating future climate and ecological change in the GYA (Whitlock and Hostetler 2019). The first was a prolonged period from about 11,500 to 7000 yr ago (the early-Holocene period) when summers in the region were on average 1.8-3.6°F (1-2°C) warmer than the pre-industrial average (Kutzbach et al. 1998: Bartlein et al. 1998). The causes of this warming were increased solar radiation during the Northern Hemisphere summer resulting from slow Milankovitch variations in the tilt of the Earth’s axis and rising levels of greenhouse gases in the atmosphere.
Pollen records indicate that the early-Holocene period in the GYA was a time of expanded lodgepole pine forest and more Douglas-fir and aspen (Populus tremuloides) compared to present. The upper tree line lay at a higher elevation than at present in response to longer growing seasons, and lower tree line shifted upslope in response to drought (Whitlock 1993; Iglesias et al. 2018). Many of the small lakes and wetlands in northern Yellowstone National Park dried during the early Holocene, and fires were more frequent. Snow and ice fields at high elevations shrank in size and accumulated plant debris and artifacts that were preserved by ice during subsequent cold periods (see box). Longer summers at Yellowstone Lake likely resulted in earlier ice-off in spring and longer open-water conditions in fall (Thompson et al. 1998).
Patches of year-round ice are found at high elevations throughout the Greater Yellowstone Area (GYA), and scientists have discovered that some of these patches are thousands of years old and preserve valuable information about the past. Recent warming has resulted in substantial melting and shrinking of the ice bodies, exposing organic artifacts that have been frozen in the ice. GYA artifacts provide unique insights into the activities of ancient hunter-gatherers in the high mountains. For example, a 10,300-year-old atlatl dart, used in hunting big game, was recovered from a melting ice patch in northwestern Wyoming (photos).
These ice patches are also a valuable source of paleoenvironmental and paleoclimatic information (Chellman et al. 2021). In 2018, a 6-m-long ice core was taken from the same ice patch where the atlatl dart shaft was found. The core contained 29 layers of plant remains (e.g., seeds, pollen, needles, and organic matter), animal dung, and dust (photo). Radiocarbon dating revealed that these debris layers formed during periods of warm and/or dry conditions that occurred on average every 300 yr over the last 10,000 yr.
A nearby melting ice patch uncovered fossil logs of whitebark pine (Pinus albicaulis), indicating that during a warm period 5000 yr ago conifers grew at elevations at least 100 m (330 ft) above present-day tree line. Tree-ring analysis of the wood showed that this warm period persisted for about 800 yr. Scientists expect more discoveries as a warming climate continues to melt old ice patches. Uncovered debris and artifacts will help us better understand past high-elevation environments, as well as the people who lived there.
The last 1000 years
The second period of warmth since the last ice age, the Medieval Climate Anomaly (800-1300), occurred on most continents, although the underlying cause of warming is not fully understood. Tree-ring records and other studies from GYA offer regional information about past temperature, precipitation, summer drought, snowpack, and streamflow over the last 1000 yr. The Medieval Climate Anomaly was overall not as warm as the early Holocene and instead was characterized by multi-decadal periods of warm summer temperatures, low snowpack, and dry conditions, which are referred to as megadroughts (Pederson et al. 2011b; Martin et al. 2019; Heeter et al. 2021).
Megadroughts occurred in the GYA and across much of the western United States in the early 600s, late 800s, 1200s, and late 1500s (Williams et al. 2020). These dry periods led to more fires, desiccation of small lakes, reduced streamflow, an upslope shift in upper tree line, and reduced Old Faithful geyser activity (see box) (Meyer et al. 1995; Millspaugh et al. 2004; Pederson et al. 2011b; Hurwitz et al. 2020).
Old Faithful Geyser got its name in the 19th century because its eruptions were both regular and predictable. Recent years of low precipitation have resulted in less frequent eruptions of Old Faithful, and this slowdown has raised concerns from the public.
To investigate this change in eruption frequency, a team of scientists were given permission by Yellowstone National Park to collect 13 mineralized specimens of lodgepole pine (Pinus contorta) wood from the Old Faithful geyser mound (Hurwitz et al. 2020). The fact that trees at one time grew on the mound suggests that the geyser was not actively erupting at some point in the past. When eruptions at Old Faithful resumed, the trees were killed and preserved in mineral deposits. Radiocarbon dating of the wood samples show that tree establishment and associated eruption hiatus occurred in the early-13th through mid-14th centuries (1233-1362).
Independent climate studies based on tree-ring records indicate a severe and sustained drought across GYA in the mid-13th century at the time the trees grew on the Old Faithful mound. The scientists hypothesize that reduced precipitation limited the subsurface supply of water to the geyser basin causing a cessation in eruptions of Old Faithful for an extended period of time.
The Medieval Climate Anomaly was followed by a period of above-average snowpack, renewed glacial activity, and cool conditions called the Little Ice Age. The Little Ice Age occurred at different times around the world, and its beginning and end are variously defined (Mann 2003; Neukom et al. 2019). In this Assessment we use the period from about 1550-1850, but note that cooling events began as early as 1300 in the GYA (Heeter et al. 2021). Cooling during the Little Ice Age may have been triggered by heightened volcanic activity, decreased solar activity, a shift in atmosphere-ocean circulation patterns, or even increased forest cover (acting as a carbon sink) during times of human population decline (Mann 2003; Ruddiman 2013). Glaciers at high elevations in the Rocky Mountains were reactivated during this period (Carrara et al. 1987; Menounos et al. 2009), and annual snowpack was high in the GYA during the years of 1535-1550, 1600-1620, 1660-1790, and 1845-1895 (Pederson et al. 2011b). Following the Little Ice Age, the lowest snowpack of the last 1000 yr occurred from 1900 to 1949 and since the 1980s (see box).
The steady decline in snowpack since the 1980s (measured as the amount of liquid water [or snow water equivalent] on April 1) is a concern for natural resource managers and communities that depend on mountain snowpack for their water supply.
While the great ecological and societal importance of mountain snowpack is clear, the observational record of mountain snowpack variability is short. Thus, scientists used records of tree growth that are sensitive to changes in snowpack across the GYA to reconstruct April 1 snow water equivalent for over the past 800 yr (Pederson et al. 2011b).
The reconstruction (see figure) shows a significant decrease in snowpack during the 20th and early 21st centuries as compared to the previous 800 yr. During the Little Ice Age (circa 1550-1850; shown in blue shading), glaciers in GYA, like elsewhere in the northern Rocky Mountains, reached their greatest extent of the Holocene as a result of persistent above‑average snowpack and cool summers. Conversely, exceptionally low snowpack during the 1930s Dust Bowl drought (shown in red shading) and since the 1980s—both attributed in part to warm summer conditions—has not been observed since at least the Medieval Climate Anomaly (800-1300).
The tree-ring based reconstructions of snowpack in the GYA indicate that variations in summer temperature govern the overall amount of snowpack that persists over the long term (decades to centuries), whereas short-term differences (year-to-year to decadal) in snowpack are caused by variability in precipitation. The snowpack reconstruction implies that the recent decades of extremely low April 1 snow water equivalent relative to the last 800 yr are associated with regional warming; warming in the future will likely continue this trend (as discussed in Chapter 7).
The tree-ring based reconstructions of snowpack in the GYA indicate that variations in summer temperature govern the overall amount of snowpack that persists over the long term (decades to centuries), whereas short-term differences (year-to-year to decadal) in snowpack are caused by variability in precipitation.
The last 120 years
Observations over the last 120 yr (1900 to present) show long-term trends in temperature and precipitation with substantial year-to-year and decadal variability, including extreme dry and wet episodes relative to average conditions (Figure 2-5C). Some notable events in the GYA associated with trends and variability over the last 120 yr include:
Trends across the 120-year period.—Mean annual temperature and precipitation in the GYA have varied over the last 120 yr with a substantial range of year-to-year variability and extended periods that were drier or wetter and colder or warmer than average (Figures 2-5A and B). After an extended dry period from 1905-1945 that included the 1930s Dust Bowl drought, precipitation has been near or above the long-term average. GYA temperatures were below the long-term average before late 1920s and then increased during the Dust Bowl years. Temperatures then dropped to near average values until the late 1970s, when they started to increase substantially. The combination of changing temperature and precipitation resulted in variable drought conditions as characterized by the Palmer Drought Severity Index (PDSI). The PDSI shows a few extreme droughts in the past 120 yr, such as in the 1930s, 1988, and early 2000s (Figure 2-5C). Extreme cold and heavy snow events that were common in the late 19th century are now rare (see box).
Decadal-scale variability: the 1930s Dust Bowl drought.—Moisture variability across the GYA is evident as wet and dry conditions that lasted for decades (highlighted by 20-year smoothing average in Figure 2-5A and C). The tendency for moisture conditions to persist over extended periods presents unique challenges for resource managers and local communities. For example, sustained low precipitation, elevated temperatures, and drought conditions during the 1930s Dust Bowl event (orange highlighted boxes in Figure 2-5) resulted in years of elevated regional fire activity, severely reduced surface water resources and streamflow, and the foreclosure and sale of many farms and ranches around the GYA (Murphy 2003). In many USGS streamgage records in the GYA, the Dust Bowl drought still ranks as one of the most severe and sustained drought events on record.
Year-to-year variability: the 1988 Yellowstone National Park fires.—Unusually little precipitation fell in 1988 (red point, Figure 2-5A), when extensive forest fires swept through Yellowstone National Park. Average temperature was high and precipitation was low in 1988 (Figure 2-5B) resulting in severe drought, as indicated by the Palmer Drought Severity Index (PDSI; Figure 2-5C). PDSI is a measure of drought intensity that accounts for both the current weather and the cumulative effects of precipitation and temperature from previous months. (See the wildfire box in Chapter 3 for more information.)
A century before the Yellowstone fires of 1988, an extreme natural disturbance of a different type occurred: the Children’s Blizzard of 1888. The morning of January 12, 1888, was warm and calm across GYA and onto the Great Plains, but these conditions abruptly changed as an Arctic cold air mass enveloped the region, causing temperatures to plummet to subzero values. Children were making their normal trek to rural schools, but when the icy weather hit en route, some attempted to return home. Many didn’t make it back to their families and perished while stranded in the storm (Potter 2012). It is estimated that 250-500 individuals died in this event (Valle undated).
The Children’s Blizzard of 1888 was one of a series of severe, cold winter storms that swept the United States—from the Rocky Mountains to the East Coast—during the late 1800s and early 1900s. These winter storms were usually preceded by relatively warm weather and characterized by sudden drops in temperature and heavy snow.
Records from Fort Sheridan (now Mammoth) in Yellowstone National Park reported severe cold snaps and heavy snow starting January 3rd and extending to January 20, 1888, the year of the Children’s Blizzard. The lowest temperature recorded was -41°F (-40.5°C) on January 14, 1888. The high temperature on that same day was -25°F (-32°C). During that month, Mammoth recorded almost 30 inches (76 cm) of snow. Thirteen out of those 17 snow days experienced lows below -20°F (-29°C) (US National Archives and Records Administration undated).
The previous year, 1887, Montana ranchers experienced high cattle losses in what is now known as “The Great Die-Up.” Heavy snows, low temperatures, and strong winds created a thick crust of ice and snow that livestock could not break through to reach the sparse grasses beneath. Lack of food and exposure to the elements proved disastrous for Montana cattle (LeCain 2017).
Today, subzero cold weather is often associated with cold Arctic air and little snow. National Weather Service data from Bozeman, Montana, and Cody, Jackson, and Mammoth, Wyoming, show periods of extended subzero cold in the last 50 yr, but typically these periods received less than 5 inches (13 cm) of snow. At these four stations, 6 days in the last half century registered temperatures below -40°F (-40°C; the low recorded during the Children’s Blizzard of 1888). All 6 days occurred in Jackson WY (Climate Analyzer undated).
In the last decade (2010-2019), there have been only five times in the GYA when 8 inches (20 cm) or more of snow accumulated in a 48-hour period that also featured subzero drops in temperature. Four of these weather events occurred in Mammoth WY (in 2010, 2014, 2017, and 2019) with the lowest temperature of -29°F (-34°C) during the 2019 storm. The low temperature recorded in Bozeman MT during this same storm was -10°F (-23°C).
No GYA weather event in the last decade measures up to the 7 days of negative temperatures and 30 inches (76 cm) of snow that was recorded between the 3rd and 20th of January 1888. Our winters have gotten warmer and the absence of the extreme, extended sub-zero periods is an indication that the climate of GYA is changing.
Causes of Climate Change
The Earth’s energy balance is driven by solar radiation that is absorbed by land surface and oceans and radiated back to the atmosphere as heat (Figure 2-6). Greenhouse gas (GHG) molecules, like water vapor (H2O) and carbon dioxide (CO2), have chemical bond structures that trap some of the heat from the Earth’s surface that otherwise would escape back to space. In this way, GHGs promote the accumulation of heat in the lower atmosphere that is necessary to sustain life. (Without atmospheric water vapor and GHGs, the global temperature would be -0.4°F [-18°C], roughly 59°F [33°C] colder than present [WMOb undated].)
The heat-trapping capacity of GHGs has been known since 19th-century laboratory studies: increasing GHG concentrations increases temperature. The ability of a gas to trap heat is determined by the amount of the gas in the atmosphere, how long the gas lasts before breaking down, and the ability of the gas to absorb (or trap) energy. Water vapor is the most abundant GHG in the atmosphere but also one of the fastest to cycle. CO2 is the second most abundant GHG and has a lifetime of 300-1000 yr (NASAb undated); its concentration recently surpassed 415 parts per million (NOAAc undated). Concentrations of other GHGs in the atmosphere are lower than CO2, but they have far greater heat-trapping ability. For example, methane (CH4), which is measured in parts per billion, is 84 times more effective at trapping heat than CO2 but it only persists in the atmosphere for about a decade (NOAAc undated).
The past changes in climate discussed in this chapter are largely the result of natural drivers that affect the Earth’s energy and moisture balances in ways that result in cooling or warming. Additional human-caused or anthropogenic climate drivers—which can reinforce or attenuate the climate response to natural drivers—include changes in land cover, and increasing emissions of greenhouse gases, sulfate aerosols, and particulate matter like ash and soot. Since 1750, human-caused climate drivers have been rapidly increasing and, in the last century, their effect exceeds that of all natural climate drivers combined (see Chapter 4). The primary anthropogenic driver is the burning of fossil fuels, as described in detail in national and international climate assessments (IPCC 2013; USGCRP 2017; Blunden and Arndt 2019). Scientific agreement that humans are the cause of current climate change is overwhelming, as summarized by NASA (NASA undated):
The vast majority of actively publishing climate scientists—97%—agree that humans are causing global warming and climate change. Most of the leading science organizations around the world have issued public statements expressing this, including international and US science academies, the United Nations Intergovernmental Panel on Climate Change and a whole host of reputable scientific bodies around the world.
Concentrations of atmospheric CO2 have been directly measured since the 1950s at the Mauna Loa Observatory in Hawaii (Figure 2-7). The concentration exceeded 415 ppm in March 2021, by far the highest level in the past 800,000 yr when natural CO2 levels ranged between 180-290 ppm (EPICA Community Members 2004). The current level of CO2 also implies that today the Earth’s climate is warmer than the last 20,000 yr, and likely warmer than previous interglacial and glacial periods in the last 800,000 yr. Recent research based on analysis of Pliocene-age CO2 levels in deep ocean sediment cores suggests that there is more CO2 in the atmosphere than at any time in the past 3.3 million years (de la Vega et al. 2020). GHG levels in the atmosphere will continue to rise unless deliberate action is taken to reverse the trend through mitigation (IPCC 2018).
Recent research based on analysis of Pliocene-age CO2 levels in deep ocean sediment cores suggests that there is more CO2 in the atmosphere than at any time in the past 3.3 million years (de la Vega et al. 2020). GHG levels in the atmosphere will continue to rise unless deliberate action is taken to reverse the trend through mitigation.
Climate is the long-term average of weather and usually measured over a base period (e.g., 20 yr, from 1986 through 2005, in this report). Climate changes gradually or abruptly lead to different long-term trends and multi-decadal averages. Shorter (e.g., annual or decadal) variability is superimposed on long-term trends. Both trends and variability can change over time, and indeed they are related but should not be confused.
Seasonal temperature and precipitation in the GYA are governed by the relative contribution of air masses from the Pacific Ocean, Arctic Ocean, and Gulf of Mexico regions through the year. Winter and spring precipitation largely comes from Pacific storms, and summer (and sometimes spring) precipitation comes from subtropical sources in the Pacific and Gulf of Mexico. Year-to-year and decadal climate patterns, such as ENSO and PDO, are attributed to large-scale atmosphere‑ocean interactions and their influence on surface climate conditions in other regions. Given the inland location of the GYA, the relationship between ENSO and 20th-century climate variability in the GYA is relatively weak.
Climate change has occurred on all timescales in the Greater Yellowstone Area. Gradual changes over thousands of years are largely driven by cyclical variations in solar radiation related to Earth’s orbit around the sun and the natural variability in atmospheric greenhouse gases. Short-term variations occurring over years to centuries are related to changes in volcanic activity, solar output, and atmosphere-ocean circulation patterns.
The high elevations of the Greater Yellowstone area were covered by a large ice field cap from 22,000-13,000 yr ago, with glaciers flowing down all the major valleys to low elevations. The climate was 5-7°F (2.8-3.9°C) colder than the pre-industrial period. Past glaciations were responsible for shaping most of the landforms that we see in the region today.
A period of warming occurred from 11,500-7000 yr ago (the early-Holocene period), when summers were 1.8-3.6°F (1-2°C) warmer than the pre-industrial period. This was a time of vegetation change, drying wetlands, more fires, and shrinking snow fields.
The Medieval Climate Anomaly, from years 800 to 1300, was a time when summers were slightly warmer than the pre-industrial period. This period was characterized by decade-long droughts that brought more fires, lower streamflow, establishment of trees above present tree line, and even a near-century hiatus of geyser activity at Old Faithful. Notable droughts occurred in the early 600s, late 800s, 1200s, and late 1500s. The Medieval Climate Anomaly was followed by cold, snowy conditions in the Little Ice Age from about 1550-1850.
Warming globally and in the GYA over the 20th and 21st centuries is attributed to increased emission of anthropogenic greenhouse gases (e.g., CO2, CH4, and others) from the burning of fossil fuels. The average temperature of the last two decades (2001-2020) is probably higher than any period in the last 20,000 yr, and likely higher than previous interglacial or glacial periods in the last 800,000 yr. The current level of carbon dioxide in the atmosphere currently is the highest in the last 3.3 million years.
Abatzoglou JT. 2011. Influence of the PNA on declining mountain snowpack in the western United States. International Journal of Climatology 31:1135–42. https://doi.org/10.1002/joc.2137.
Bartlein PJ, Anderson PM, Anderson KH, Edwards ME, Thompson RS, Webb RS, Webb T III, Whitlock C. 1998. Paleoclimate simulations for North America over the past 21,000 years: features of simulated climate and comparisons with paleoenvironmental data. Quaternary Science Reviews 17:549-85.
Bartlein PJ, Whitlock C, Shafer S. 1997. Future climate in the Yellowstone National Park region and its potential impact on vegetation. Conservation Biology 11:782-92.
Blunden J, Arnd DS (editors). 2019. State of the climate in 2018. Bulletin of the American Meteorological Society 100(9):Si–S306. Available online https://doi.org/10.1175/2019BAMSStateoftheClimate.1. Accessed 13 May 2021.
Carrara, PE. 1987. Holocene and the latest Pleistocene glacial chronology, Glacier National Park, Montana. Canadian Journal of Earth Sciences 24:387–95.
Chang T, Hansen A, Piekielek N. 2014. Patterns and variability of suitable bioclimate habitat for Pinus albicaulis under multiple projected climate models. PLOS ONE 9(11):e111669.
Chellman NJ, Pederson GT, Lee CM, McWethy DB, Puseman K, Stone JR, Brown SR, McConnell JR. 2021. High elevation ice patch documents Holocene climate variability in the northern Rocky Mountains. Quaternary Science Advances 3:1000021. https://doi.org/10.1016/j.qsa.2020.100021.
Clark PU, Dyke AS, Shakun JD, Carlson AE, Clark J, Wohlfarth B, Mitrovic J, Hostetler SW, McCabe AM. 2009. The last glacial maximum. Science 325:710-14.
Clark PU, Shakun JD, Baker PA, Bartlein PJ, Brewer S, Brook EJ, Carlson AE, Cheng H, Kaufman DS, Liu Z, Marchitto TM, Mix AC, Morrill C, Otto-Bliesner B, Pahnke K, Russell JM, Whitlock C, Adkins JF, Blois J, Colman SC, Curry WN, Flower BP, He F, Johnson TC, Lynch-Stieglitz J, Markgraf V, McManus JF, Mitrovica JX, Moreno PI, Williams JW. 2012. Global climate evolution during the last deglaciation. Proceedings of the National Academy of Sciences USA 109 (19) E1134-E1142. doi:10.1073/pnas.1116619109.
Climate Analyzer. [undated]. The climate analyzer [website]. Available online http://www.climateanalyzer.org. Accessed 26 Feb 2021.
Davies B. 2020. Cosmic nuclide dating [website of ArcticGlaciers.org]. Available online www.antarcticglaciers.org/glacial-geology/dating-glacial-sediments-2/cosmogenic_nuclide_datin/. Accessed 16 Feb 2021.
de la Vega E, Chalk TB, Wilson PA, Bysani RP, Foster GL. 2020. Atmospheric CO2 during the Mid-Piacenzian Warm Period and the M2 glaciation. Scientific Reports 10:11002. https://doi.org/10.1038/s41598-020-67154-8.
EPICA Community Members. 2004. Eight glacial cycles from an Antarctic ice core. Nature 429:623-8.
Farnes P. 1997. The snows of Yellowstone. Yellowstone Science 5:8-11.
Good JD, Pierce KL. 1996. Interpreting the landscape: recent and ongoing geology of Grand Teton and Yellowstone national parks (2nd printing, revised and reprinted, 1998). Moose WY: Grand Teton Association. 58 p.
Heeter KJ, Rochner ML, Harley GL. 2021. Summer air temperature for the Greater Yellowstone Ecoregion (770-2019 CE) over 1250 years. Geophysical Research Letters 48:e2020GL092269. https//doi.org/10.1029/2020GL092269.
Hurwitz S, King JC, Pederson GT, Martin JT, Damby DE, Manga M, Hungerford JDG, Peek S. 2020. Yellowstone’s Old Faithful Geyser shut down by a severe thirteenth century drought. Geophysical Research Letters. doi: 10.1029/2020GL089871.
Iglesias V, Whitlock C, Krause TR, Baker RG. 2018. Reconstructing past ecosystem dynamics in the Greater Yellowstone Ecosystem region based on modern pollen-vegetation relationships. Journal of Biogeography 45:1768-80. doi:10.1111/jbl.13364.
[IPCC] International Panel on Climate Change. 2013. In: Stocker TF, Qin D, Plattner G-K, Tignor M, Allen SK, Boschung J, Nauels A, Xia Y, Bex V, Midgley PM, editors. Climate change 2013: the physical science basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge UK and New York NY: Cambridge University Press. 1535 p. Available online https://www.ipcc.ch/report/ar5/wg1/. Accessed 8 Mar 2021.
[IPCC] International Panel on Climate Change. 2018. In: Masson-Delmotte V, Zhai P, Pörtner H-O, Roberts D, Skea J, Shukla PR, Pirani A, Moufouma-Okia W, Péan C, Pidcock R, Connors S, Matthews JBR, Chen Y, Zhou X, Gomis MI, Lonnoy E, Maycock T, Tignor M, Waterfield T, editors. Global Warming of 1.5°C: an IPCC Special Report on the impacts of global warming of 1.5°C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty. 630 p. Available online https://www.ipcc.ch/site/assets/uploads/sites/2/2019/06/SR15_Full_Repor…. Accessed 8 Mar 2021.
Keeling RF, Keeling CD. 2017. Atmospheric monthly in situ CO2 data—Mauna Loa Observatory, Hawaii [webpage]. In Scripps CO2 program data; UC San Diego Library Digital Collections. Available online https://doi.org/10.6075/J08W3BHW. Accessed 9 Mar 2021.
Kennedy C. 2012 (Sep 28). NOAA Climate.Gov: El Niño and US winter weather [webpage]. Available online https://www.climate.gov/news-features/featured-images/el-ni%C3%B1o-and-us-winter-weather. Accessed 9 Mar 2021.
Krause TR, Whitlock C. 2017. Climatic and non-climatic controls shaping early postglacial conifer history in the northern Greater Yellowstone Ecosystem, USA. Journal of Quaternary Science 32:1022-36.
Kutzbach J, Gallimore R, Harrison S, Behling P, Selin R, Laarif F. 1998. Climate and biome simulations for the past 21,000 years. Quaternary Science Reviews 17:473-506.
LeCain TJ. 2017. The matter of history: how things create the past. Cambridge UK: Cambridge University Press. 366 p.
Le Treut H, Somerville R, Cubasch U, Ding Y, Mauritzen C, Mokssit A, Peterson T, Prather M. 2007. Historical Overview of Climate Change [chapter 1]. In: Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL, editors. Climate change 2007: the physical science basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. p 95-127. Cambridge UK and New York NY: Cambridge University Press.
Licciardi J, Pierce, KL. 2018. History and dynamics of the Greater Yellowstone glacial system during the last two glaciations. Quaternary Science Reviews 200:1-33. https://doi.org/10.1016/j.quascirev.2018.08.027.
Mann ME. 2003. Little Ice Age [section]. In: MacCracken MC, Perry JS, Munn T, editors. Encyclopedia of Global Environmental Change, volume 1: The Earth system—physical and chemical dimensions of global environmental change. p 504-9. Hoboken NJ: Wiley and Sons.
Martin JT, Pederson GT, Woodhouse CA, Cook ER, McCabe GJ, Wise EK, Erger P, Dolan L, McGuire M, Gangopadhyay S, Chase K, Littell JS, Gray ST, St. George S, Friedman J, Sauchy D, St Jacques J, King J. 2019. 1200 years of upper Missouri River streamflow reconstructed from tree rings. Quaternary Science Reviews 224:105971. https://doi.org/10.1016/j.quascirev.2019.105971.
McCabe GJ, Betancourt JL, Gray ST, Palecki MA, Hidalgo HG. 2008. Associations of multi-decadal sea-surface temperature variability with US drought. Quaternary International 188(1):31-40.
McCabe GJ, Palecki MA, Betancourt JL. 2004. Pacific and Atlantic Ocean influences on multidecadal drought frequency in the United States. Proceedings of the National Academy of Sciences USA 101:4136-41.
Menounos B, Osborn G, Clague JJ, Luckman BH. 2009. Latest Pleistocene and Holocene glacier fluctuations in western Canada. Quaternary Science Reviews 28:2049-74.
Meyer GA, Wells SG, Jull AJT. 1995. Fire and alluvial chronology in Yellowstone National Park: climatic and intrinsic controls on Holocene geomorphic processes. Geological Society of America Bulletin 107:1211-30.
Millspaugh SH, Whitlock C, Bartlein PJ. 2004. Postglacial fire, vegetation, and climate history of the Yellowstone-Lamar and Central Plateau provinces, Yellowstone National Park. In: Wallace L, editor. After the fires: the ecology of change in Yellowstone National Park. p 10-28. New Haven CT: Yale University Press.
Murphy M. 2003. Hope in hard times: new deal photographs of Montana, 1936-1942. Helena MT: Montana Historical Press Society. 256 p.
[NASA] National Aeronautics and Space Administration. [undated]. Do scientists agree on climate change [webpage]? Available online https://climate.nasa.gov/faq/17/do-scientists-agree-on-climate-change/. Accessed 9 Mar 2021.
Neukom R, Steiger N, Gomez-Navarro JJ, Want J, Werner JP. 2019. No evidence for globally coherent warm and cold periods over the preindustrial Common Era. Nature 571:550-4.
[NOAAa] National Oceanic and Atmospheric Administration. [undated]. National Centers for Environmental Information: Palmer Drought Severity Index data. Available online https://psl.noaa.gov/cgi-bin/data/timeseries/timeseries.pl?ntype=2&type…. Accessed 29 May 2020.
[NOAAa] National Oceanic and Atmospheric Administration. [undated]. National Centers for Environmental Information: Palmer Drought Severity Index data. Available online https: //psl.noaa.gov/cgi-bin/data/timeseries/timeseries.pl ?ntype=2&typediv=3&state=+48&averaged=11&division=1&year1=1895&year2=2019&anom=1&iseas=1&mon1=0&mon2=11&typeout=1&y1=&y2=&plotstyle=0&Submit=Create+Timeseries. Accessed 29 May 2020.
[NOAAb] National Oceanic and Atmospheric Administration. [undated]. The atmosphere: getting a handle on carbon dioxide [webpage]. Available online https://climate.nasa.gov/news/2915/the-atmosphere-getting-a-handle-on-c…. Accessed 9 Mar 2021.
[NOAAc] National Oceanic and Atmospheric Administration. [undated]. The NOAA annual greenhouse gas index [webpage]. Available online https://www.esrl.noaa.gov/gmd/aggi/aggi.html. Accessed 9 Mar 2021.
Pederson GT, Betancourt JL, McCabe GJ. 2013. Regional patterns and proximal causes of the recent snowpack decline in the Rocky Mountains. US Geophysical Research Letters 40:1811-6.
Pederson GT, Gray ST, Ault T, Marsh W, Fagre DB, Bunn AG, Woodhouse CA, Graumlich LJ. 2011a. Climatic controls on the snowmelt hydrology of the northern Rocky Mountains. Journal of Climate 24(6):1666-87.
Pederson GT, Gray ST, Woodhouse CA, Betancourt JL, Fagre DB, Littell JS, Watson E, Luckman BH, Graumlich LJ. 2011b. The unusual nature of recent snowpack declines in the North American Cordillera. Science 333:332-5.
Potter S. 2012. Retrospect: January 12, 1888: the Children’s Blizzard. Weatherwise 65(1):10-1. Available online https://doi.org/10.1080/00431672.2012.635992. Accessed 9 Mar 2021.
PRISM Climate Group. 2020. PRISM climate data [data source]. Available online http://prism.oregonstate.edu. Accessed 20 Dec 2020.
Rasmussen M, Anzick SL, Waters MR, Skoglund P, DeGiorgio M, Stafford Jr TW, Rasmussen S, Moltke I, Albrechtsen A, Doyle SM, and others. 2014. The genome of a late Pleistocene human from a Clovis burial site in western Montana. Nature 506:225–9.
Ruddiman WF. 2013. Earth’s climate: past and future. New York NY: WH Freeman and Company. 465 p.
Thompson RS, Hostetler SW, Bartlein PJ, Anderson KH. 1998. A strategy for assessing potential future changes in climate, hydrology, and vegetation in the western United States. US Geological Survey Circular 1153. Washington DC: US Government Printing Office. 20 p.
[USGCRP] US Global Change Research Program. 2017. Wuebbles DJ, Fahey DW, Hibbard KA, Dokken DJ, Steward BC, Maycock TK, editors. Climate science special report: fourth national climate assessment, vol 1. Washington DC: USGCRP. 470 p. doi:10.7930/J0J964J6.
US National Archives and Records Administration [undated]. Meteorological records, 1888. Record group 27; records of the Weather Bureau. Available by microfilm; see information at https://www.archives.gov/research/guide-fed-records/groups/027.html#27.2. Accessed 9 Mar 2021.
Valle D. [undated]. The Blizzards of 1888. National Weather Service Heritage. Available online https://vlab.ncep.noaa.gov/web/nws-heritage/-/the-children-s-blizzard. Accessed 9 Mar 2021.
Whitlock C. 1993. Postglacial vegetation and climate of Grand Teton and southern Yellowstone national parks. Ecological Monographs 63:173-98.
Whitlock C, Bartlein PJ. 1993. Spatial variations of Holocene climate change in the Yellowstone region. Quaternary Research 39(2):231-8.
Whitlock C, Hostetler S. 2019. Past warm periods provide vital benchmarks for understanding the future of the Greater Yellowstone Ecosystem. Yellowstone Science 27:72-6.
Williams AP, Cook ER, Smerdon JE, Cook BI, Abatzoglou JT, Bolles K, Baek SH, Badger AM, Livneh, B. 2020. Large contribution from anthropogenic warming to an emerging North American megadrought. Science 368:314-8.
[WMOa] World Meteorological Organization. [undated]. Essential climate variables [webpage]. Available online https://public.wmo.int/en/programmes/global-climate-observing-system/es…. Accessed 26 Dec 2020.
[WMOb] World Meteorological Organization. [undated]. Greenhouse gases [webpage]. Available online https://public.wmo.int/en/our-mandate/focus-areas/environment/greenhous…. Accessed 9 Mar 2021.
 Cosmogenic nuclide dating uses the interactions between cosmic rays and the atomic nuclides found in glacially transported boulders to provide age estimates for the rock’s exposure at the Earth’s surface (Davies 2020). In other words, cosmogenic nuclide dating determines how long the boulders in moraines have been at the surface, which in turn provides the age of glacier position.