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: Milanković Forcing in Deep Time
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
Award ID(s):
2034660
PAR ID:
10513535
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Paleoceanography and Paleoclimatology
Volume:
39
Issue:
5
ISSN:
2572-4517
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract The dynamical evolution of the solar system is chaotic with a Lyapunov time of only ∼5 Myr for the inner planets. Due to the chaos it is fundamentally impossible to accurately predict the solar system’s orbital evolution beyond ∼50 Myr based on present astronomical observations. We have recently developed a method to overcome the problem by using the geologic record to constrain astronomical solutions in the past. Our resulting optimal astronomical solution (called ZB18a) shows exceptional agreement with the geologic record to ∼58 Ma (Myr ago) and a characteristic resonance transition around 50 Ma. Here we show that ZB18a and integration of Earth’s and Mars’ spin vector based on ZB18a yield reduced variations in Earth’s and Mars’ orbital inclination and Earth’s obliquity (axial tilt) from ∼58 to ∼48 Ma—the latter being consistent with paleoclimate records. The changes in the obliquities have important implications for the climate histories of Earth and Mars. We provide a detailed analysis of solar system frequencies (gandsmodes) and show that the shifts in the variation in Earth’s and Mars’ orbital inclination and obliquity around 48 Ma are associated with the resonance transition and caused by changes in the contributions to the superposition ofsmodes, plusg–smode interactions in the inner solar system. Theg–smode interactions and the resonance transition (consistent with geologic data) are unequivocal manifestations of chaos. Dynamical chaos in the solar system hence not only affects its orbital properties but also the long-term evolution of planetary climate through eccentricity and the link between inclination and axial tilt. 
    more » « less
  3. Lacustrine strata are often among the highest-resolution terrestrial paleoclimate archives available. The manner in which climate signals are registered into lacustrine deposits varies, however, as a function of complex sedimentologic and diagenetic processes. The retrieval of reliable records of climatic forcing therefore requires a means of evaluating the potential influence of changing sedimentary transfer functions. Here, we use high-resolution X-ray fluorescence core scanning of the Wilkins Peak Member of the Green River Formation to characterize the long-term evolution of transfer functions in an ancient lacustrine record. Our analysis identifies a shift in the frequency distribution of Milankovitch-band variance between the lower and middle Wilkins Peak Member across a range of temporally calibrated elemental intensity records. Spectral analysis of the lower Wilkins Peak Member shows strong short eccentricity, obliquity, precession, and sub-Milankovitch−scale variability, while the middle Wilkins Peak Member shows strong eccentricity variability and reduced power at higher frequencies. This transition coincides with a dramatic decline in the number and volume of evaporite beds. We attribute this shift to a change in the Wilkins Peak Member depositional transfer function caused by evolving basin morphology, which directly influenced the preservation of bedded evaporite as the paleolake developed from a deeper, meromictic lake to a shallower, holomictic lake. The loss of bedded evaporite, combined with secondary evaporite growth, results in reduced obliquity- and precession-band power and enhanced eccentricity-band power in the stratigraphic record. These results underscore the need for careful integration of basin and depositional system history with cyclostratigraphic interpretation of the dominant astronomical signals preserved in the stratigraphic archive. 
    more » « less
  4. Marine δ18O data reveal astronomical forcing of the climate and cryosphere during the Miocene, when atmosphericPco2was on par with emissions scenarios over the next century. This inspired hypotheses for how Milankovitch cycles, ice-ocean interactions, and greenhouse gases influence ice volume. Mass balance controls for marine and terrestrial ice sheets differ, and proxy data collected far from Antarctica provide valuable but limited insight into regional processes. We evaluate clast abundance data from Antarctic marine sedimentary records, observing a strong signal of eccentricity and precession coincident with a terrestrial ice sheet and a clear obliquity signal at the margins of a marine ice sheet. These analyses are integrated with a synthesis of proxy data, and we argue that high variance in obliquity forcing (mediated and enhanced by the ocean and atmosphere) can inhibit ice sheet growth, even when insolation forcing is conducive to glaciation. This “obliquity disruption” explains cryosphere variability before the existence of large northern hemisphere ice sheets. 
    more » « less
  5. 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