USGS Science⌕ Search

USGS · 70209450

History and dynamics of the Greater Yellowstone Glacial System during the last two glaciations

Abstract

The Greater Yellowstone Glacial System (GYGS) covered about 20,000 km2 at its maximum Pleistocene extent. The initiation, culmination, and ultimate decay of the GYGS involved complex interactions between several coalescent ice masses flowing from glacial source areas adjoining and including the Yellowstone Plateau. Here, we present an updated review of the history and dynamics of the GYGS during the penultimate (Bull Lake) and last (Pinedale) glaciations, drawing upon an integration of glacial-geologic mapping with >130 cosmogenic 10Be and 3He exposure ages. Bull Lake glacial deposits in greater Yellowstone are dated to ca. 150–140 ka and correlate with marine isotope stage 6. The Bull Lake glaciation extended well beyond the Pinedale along the southern and western GYGS margins, but Pinedale glaciers overrode Bull Lake ice limits on the north and east sides. The northeastward shift of the center of ice mass from Bull Lake to Pinedale time may be explained by uplift on the leading edge of the Yellowstone hotspot and subsidence on the trailing margin. In early Pinedale time (∼22–18 ka), ice buildup culminated in the high terrain of the Beartooth Uplift and High Absaroka Range. Glaciers from these source regions flowed onto the northeastern margin of the Yellowstone Plateau and advanced to terminal moraines beyond Clarks Fork Canyon and in Jackson Hole. By middle Pinedale time (∼18–16 ka), the Yellowstone Plateau ice cap surface had risen above the equilibrium-line altitude, stimulating orographic glacial buildup nourished by storms funneled eastward through the Snake River Plain. The plateau ice cap eventually thickened to >1000 m and joined glaciers from the Beartooth Uplift and Gallatin Range to form the northern Yellowstone outlet glacier. Terrain east and downwind of the ice cap crest was placed in a precipitation shadow, resulting in glacial recession in these regions. During the late Pinedale (∼16–13 ka), the plateau ice cap prograded southwestward toward the direction of moisture supply, leading to advances along the southern and western margins of the GYGS. Northern sectors of the plateau ice cap were nearly stagnant at this time. The Yellowstone region experienced widespread deglaciation ca. 15–14 ka in response to warming climate. Unloading of the ∼1 km-thick plateau ice cap and consequent release of pressure on the magmatic system beneath Yellowstone was not accompanied by volcanism, indicating that the magma chamber was not primed to erupt via decompression during the last deglaciation.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 43.36512572875844° to 45.72152152227954° latitude; -111.796875° to -108.8525390625° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Joseph M. Licciardi, Kenneth L. Pierce. 2018. History and dynamics of the Greater Yellowstone Glacial System during the last two glaciations. https://doi.org/10.1016/j.quascirev.2018.08.027

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

A 30,000-year continuous record of glacial history and water balance from Fish Lake, Utah documents the onset, development, and demise of LGM glaciation

Fish Lake provides a continuous, high-resolution freshwater paleoclimate archive of glacial and hydroclimate history from a catchment straddling the Colorado Plateau–Great Basin boundary in the western US. At >3600 m, the Fish Lake Plateau glaciates during stadials, making it an ideal site for tracking the region's glacial history. The high-resolution (>3.3 mm/decade) record contains multiple climate- and ecosystem-sensitive proxies. Geochronology relies on 210 Pb and AMS 14 C measurements on terrestrial macrofossils and charcoal, tephrochronology, and paleomagnetic measurements. Multiproxy analyses include indicators of glacial erosion (magnetic susceptibility, sediment density, % mineral matter, grainsize, scanning XRF) and biological productivity (biogenic silica, organic matter). These proxies show glacial inception during the local Last Glacial Maximum at 25.8 ka, coincident with Heinrich Stadial 2, with maximum down-valley ice expansion between 24.8 and 20.9 ka. Glaciers reached the shoreline but did not override the lake, and sedimentation rates did not exceed Holocene levels, as increased glacial erosion was offset by reduced productivity in the glaciated, turbid lake. Deglaciation began by 18.2 ka, was interrupted by a readvance under wetter conditions synchronous with Heinrich Stadial 1 and concluded with final catchment deglaciation ∼14.1 ka. After 13.7 ka, productivity proxies rise rapidly and remain high through the Holocene while erosion proxies drop abruptly. Comparison with the adjacent Bonneville pluvial lake-level record shows that glacial onset on the Fish Lake Plateau preceded the rise in Bonneville water levels at ∼24 ka by roughly a millennium, implying temperature decreases led P-E increases. Similarly, deglaciation onset preceded falling Bonneville water levels by several centuries.

Utah↗

Sedimentological reconstruction of hurricane storm surges from a carbonate reef coastline at Anegada, British Virgin Islands

Anegada is a small, low-lying carbonate island in the British Virgin Islands that is highly vulnerable to marine inundation from frequent hurricanes. A major challenge in forecasting the timing, intensity, and impacts of future hurricanes to islands such as Anegada is the limited short-term instrumental record that does not capture the full spectrum of inundation events from the frequent but low-magnitude to rare, but high-magnitude events. Washover deposits, which are transported onshore by hurricane-related storm surges and preserved within the onshore stratigraphy, can provide insight into the recurrence of these events. Here, we use grain size distributions, X-ray fluorescence (XRF) data, and an inverse sediment transport model to characterize a series of anomalous sand deposits from shallow hypersaline ponds to determine the variability in storm surges impacting Anegada over at least the last 70 years. The analyses were successfully applied to modern analogues of storms impacting Anegada's southern shoreline (Hurricanes Donna in 1960 at Category 3; Earl in 2010 at Category 4; and Irma in 2017 at Category 5). These deposits are identified by coarse-grained material and elevated marine geochemical signals (Ca, Sr, ln(Br/Ti)). Inverse modelled flow parameters (average flow depth = 0.89–1.46 m relative to the reef flat, average flow speed = 2.95–3.78 m/s) agree with observational data and eyewitness accounts, and the deposit associated with Hurricane Donna indicates the greatest impacts occurred along the southern shoreline (modelled flow depth = 1.24–1.46 m relative to the reef flat, flow speed = 3.49–3.78 m/s). Using the modern analogues as a basis for comparison, we infer the origin of an older sand deposit to be from Hurricane Dog in 1950, when it tracked by Anegada as a Category 3 storm. This methodology provides a transferable framework for interpreting storm recurrence in carbonate-dominated coastlines worldwide that could be used to inform hurricane hazard assessments and risk management efforts.

Anegada↗

Temperature and precipitation source variability and glacial dynamics in the southwestern United States at Fish Lake, Utah, since late MIS 4

An improved understanding of hydroclimate response to natural climate forcings is needed to evaluate impacts of changing seasonal circulation patterns on current and future drought in the southwestern United States (‘the Southwest’). However, few terrestrial records from the Southwest preserve changes in moisture availability continuously prior to the Last Glacial Maximum (LGM, ∼24.5 ka). Fish Lake, in the Upper Colorado River Basin, preserves a long (∼67 ka), continuous, high-resolution (∼3.3 cm/century) sedimentary record spanning from late Marine Isotope Stage (MIS) 4 to present. Here, we investigate the relative contributions of seasonal precipitation and moisture source using the hydrogen isotope value of the n -C 29 leaf wax alkane and changes in temperature using branched glycerol dialkyl glycerol tetraethers (brGDGTs). During the warm post glacial period (14 ka to present), a greater contribution of annual precipitation is sourced from the North American Monsoon (NAM) in response to declining winter insolation and rising greenhouse gases, weakening westerly storms and shifting them northward. MIS 3 (∼61 to 29 ka) is also influenced by NAM inputs, however, to a smaller magnitude than during the post glacial period when temperatures were warmer. However, stadial periods (MIS 2, 29 to 14.1 ka) show 2 H-enriched values, suggesting warm conditions and large contributions of NAM precipitation during the LGM, contradicting other proxies and regional records. The LGM is the only period in the record where glaciers are present in the catchment, thus we hypothesize glacial erosion reworked pre-LGM aged soil organic matter and biased the record during this time.

Utah↗