skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Eccentricity Paces Late Pleistocene Glaciations
Abstract Late Pleistocene glacial terminations are caused by rising atmospheric CO2occurring in response to atmospheric and ocean circulation changes induced by increased discharge from Northern Hemisphere ice sheets. While climate records place glacial terminations coincident with decreasing orbital precession, it remains unclear why a specific precession minimum causes a termination. We compare the orbital and ice volume configuration at each precession minima over the last million years to demonstrate that eccentricity, through its control on precession amplitude, period and coherence with obliquity, along with ice sheet size, determine whether a given precession minimum will cause a termination. We also demonstrate how eccentricity controls obliquity maxima and precession minima coherence, varying the duration of glaciations. Glaciations lasting ∼100 thousand years are controlled by Earth's eccentricity cycle of the same period, while the shortest (20–40 ka) and longest (155 ka) occupy the maxima and minimums of the 400 thousand year eccentricity cycle.  more » « less
Award ID(s):
2042495
PAR ID:
10614102
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
51
Issue:
12
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract. It is widely accepted that orbital variations areresponsible for the generation of glacial cycles during the latePleistocene. However, the relative contributions of the orbital forcingcompared to CO2 variations and other feedback mechanisms causing thewaxing and waning of ice sheets have not been fully understood. Testingtheories of ice ages beyond statistical inferences, requires numericalmodeling experiments that capture key features of glacial transitions. Here,we focus on the glacial buildup from Marine Isotope Stage (MIS) 7 to 6covering the period from 240 to 170 ka (ka: thousand years before present). Thistransition from interglacial to glacial conditions includes one of thefastest Pleistocene glaciation–deglaciation events, which occurred during MIS 7e–7d–7c (236–218 ka). Using a newly developed three-dimensional coupledatmosphere–ocean–vegetation–ice sheet model (LOVECLIP), we simulate thetransient evolution of Northern Hemisphere and Southern Hemisphere ice sheets duringthe MIS 7–6 period in response to orbital and greenhouse gas forcing. For arange of model parameters, the simulations capture the evolution of globalice volume well within the range of reconstructions. Over the MIS 7–6period, it is demonstrated that glacial inceptions are more sensitive toorbital variations, whereas terminations from deep glacial conditions needboth orbital and greenhouse gas forcings to work in unison. For someparameter values, the coupled model also exhibits a critical North Americanice sheet configuration, beyond which a stationary-wave–ice-sheettopography feedback can trigger an unabated and unrealistic ice sheetgrowth. The strong parameter sensitivity found in this study originates fromthe fact that delicate mass imbalances, as well as errors, are integratedduring a transient simulation for thousands of years. This poses a generalchallenge for transient coupled climate–ice sheet modeling, with suchcoupled paleo-simulations providing opportunities to constrain suchparameters. 
    more » « less
  2. Abstract During the last glacial period, the tropical hydrological cycle exhibited large variability across orbital and millennial timescales. However, the response of the Indian summer monsoon (ISM), its related impact on terrestrial ecosystems, and associated forcing mechanisms remain controversial. Here we present a marine record of pollen‐inferred vegetation changes suggesting that eastern India shifted from woody‐savanna mosaics during Marine Isotopic Stage 3 to grasslands during the Last Glacial Maximum resulting from large‐scale drying. Our data shows that ISM maximum is in phase with obliquity and precession maxima suggesting a dominant role of the Indian Ocean interhemispheric temperature gradient on glacial ISM variability. Persistent and abrupt dryland expansions of varying magnitude suggest rapid‐scale onset of aridity during Heinrich Stadial events and during the Toba eruption. We propose that the amplitude of ISM drought events are initiated by high latitude and volcanic forcings, although modulated by precession. 
    more » « less
  3. Abstract Astronomical cycles are strongly expressed in marine geological records, providing important insights into Earth system dynamics and an invaluable means of constructing age models. However, how various astronomical periods are filtered by the Earth system and the mechanisms by which carbon reservoirs and climate components respond, particularly in absence of dynamic ice sheets, is unclear. Using an Earth system model that includes feedbacks between climate, ocean circulation, and inorganic (carbonate) carbon cycling relevant to geological timescales, we systematically explore the impact of astronomically‐modulated insolation forcing and its expression in model variables most comparable to key paleoceanographic proxies (temperature, the δ13C of inorganic carbon, and sedimentary carbonate content). Temperature predominately responds to obliquity and is little influenced by the modeled carbon cycle feedbacks. In contrast, the cycling of nutrients and carbon in the ocean generates significant precession power in atmospheric CO2, benthic ocean δ13C, and sedimentary wt% CaCO3, while inclusion of marine sedimentary and weathering processes shifts power to the long eccentricity period. Our simulations produce reducedpCO2and dissolved inorganic carbon δ13C at long eccentricity maxima and, contrary to early Cenozoic marine records, CaCO3preservation in the model is enhanced during eccentricity modulated warmth. Additionally, the magnitude of δ13C variability simulated in our model underestimates marine proxy records. These model‐data discrepancies hint at the possibility that the Paleogene silicate weathering feedback was weaker than modeled here and that additional organic carbon cycle feedbacks are necessary to explain the full response of the Earth system to astronomical forcing. 
    more » « less
  4. null (Ed.)
    Abstract. Investigating the precise timing of regional-scale climate changes during glacial terminations and the interglacial periods that follow is key tounraveling the mechanisms behind these global climate shifts. Here, we present a high-precision time series of climate changes in the Austrian Alpsthat coincide with the later portion of Termination III (TIII), the entire penultimate interglacial (Marine Isotope Stage (MIS) 7), Termination IIIa(TIIIa), and the penultimate glacial inception (MIS 7–6 transition). Using state-of-the-art mass spectrometry techniques, we have constructed auranium-series chronology with relative age uncertainties averaging 1.7 ‰ (2σ) for our study period (247 to 191 thousand yearsbefore present, ka). Results reveal the onset of warming in the Austrian Alps associated with TIII at 242.5 ± 0.2 ka and theduration of MIS 7e warming between 241.8 and 236.7 (±0.6) ka. An abrupt shift towards higher δ18O values at216.8 ka marks the onset of regional warming associated with TIIIa. Two periods of high δ18O values (greater than−10 ‰ Vienna Pee Dee Belemnite (VPDB)) between 215.9–213.3 and 204.3–197.5 (±0.4) ka coincide with interglacial substages MIS 7c and 7a,respectively. Multiple fluorescent inclusions suggest a partial retreat of the local Alpine glacier during peak obliquity forcings at214.3 ± 0.4 ka. Two newly collected stalagmites from Spannagel Cave (SPA146 and 183) provide high-resolution replications of thelatter portion of the MIS 7a-to-6e transition. The resulting multi-stalagmite record reveals important chronological constraints on climate shifts inthe Austrian Alps associated with MIS 7 while offering new insight into the timing of millennial-scale changes in the North Atlantic realm leadingup to TIII and TIIIa. 
    more » « less
  5. Abstract Astronomical (or Milanković) forcing of the Earth system is key to understanding rhythmic climate change on time scales ≳104 y. Paleoceanographic and paleoclimatological applications concerned with past astronomical forcing rely on astronomical calculations (solutions), which represent the backbone of cyclostratigraphy and astrochronology. Here we present state‐of‐the‐art astronomical solutions over the past 3.5 Gyr. Our goal is to provide tuning targets and templates for interpreting deep‐time cyclostratigraphic records and designing external forcing functions in climate models. Our approach yields internally consistent orbital and precession‐tilt solutions, including fundamental solar system frequencies, orbital eccentricity and inclination, lunar distance, luni‐solar precession rate, Earth's obliquity, and climatic precession. Contrary to expectations, we find that the long eccentricity cycle (LEC) (previously assumed stable and labeled “metronome,” recent period ∼405 kyr), can become unstable on long time scales. Our results reveal episodes during which the LEC is very weak or absent and Earth's orbital eccentricity and climate‐forcing spectrum are unrecognizable compared to the recent past. For the ratio of eccentricity‐to‐inclination amplitude modulation (recent individual periods of ~2.4 and ~1.2 Myr, frequently observable in paleorecords) we find a wide distribution around the recent 2:1 ratio, that is, the system is not restricted to a 2:1 or 1:1 resonance state. Our computations show that Earth's obliquity was lower and its amplitude (variation around the mean) significantly reduced in the past. We therefore predict weaker climate forcing at obliquity frequencies in deep time and a trend toward reduced obliquity power with age in stratigraphic records. For deep‐time stratigraphic and modeling applications, the orbital parameters of our 3.5‐Gyr integrations are made available at 400‐year resolution. 
    more » « less