skip to main content


Title: Solar System chaos and the Paleocene–Eocene boundary age constrained by geology and astronomy
Astronomical calculations reveal the Solar System’s dynamical evolution, including its chaoticity, and represent the backbone of cyclostratigraphy and astrochronology. An absolute, fully calibrated astronomical time scale has hitherto been hampered beyond ~50 million years before the present (Ma) because orbital calculations disagree before that age. Here, we present geologic data and a new astronomical solution (ZB18a) showing exceptional agreement from ~58 to 53 Ma. We provide a new absolute astrochronology up to 58 Ma and a new Paleocene–Eocene boundary age (56.01 ± 0.05 Ma). We show that the Paleocene–Eocene Thermal Maximum (PETM) onset occurred near a 405-thousand-year (kyr) eccentricity maximum, suggesting an orbital trigger. We also provide an independent PETM duration (170 ± 30 kyr) from onset to recovery inflection. Our astronomical solution requires a chaotic resonance transition at ~50 Ma in the Solar System’s fundamental frequencies.  more » « less
Award ID(s):
1658023
NSF-PAR ID:
10133086
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Science
Volume:
365
Issue:
6456
ISSN:
0036-8075
Page Range / eLocation ID:
926 to 929
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The chronology of the Paleocene-Eocene Thermal Maximum (PETM, ~56 Ma) remains disputed, hampering complete understanding of the possible trigger mechanisms of this event. Here we present an astrochronology for the PETM carbon isotope excursion from Howards Tract, Maryland a paleoshelf environment, on the mid-Atlantic Coastal Plain. Statistical evaluation of variations in calcium content and magnetic susceptibility indicates astronomical forcing was involved and the PETM onset lasted about 6 kyr. The astrochronology and Earth system modeling suggest that the PETM onset occurred at an extreme in precession during a maximum in eccentricity, thus favoring high temperatures, indicating that astronomical forcing could have played a role in triggering the event. Ca content data on the paleo-shelf, along with other marine records, support the notion that a carbonate saturation overshoot followed 
    more » « less
  2. Abstract

    Cyclostratigraphy and astrochronology are now at the forefront of geologic timekeeping. While this technique heavily relies on the accuracy of astronomical calculations, solar system chaos limits how far back astronomical calculations can be performed with confidence. High‐resolution paleoclimate records with Milankovitch imprints now allow reversing the traditional cyclostratigraphic approach: Middle Eocene drift sediments from Newfoundland Ridge are well‐suited for this purpose, due to high sedimentation rates and distinct lithological cycles. Per contra, the stratigraphies of Integrated Ocean Drilling Program Sites U1408–U1410 are highly complex with several hiatuses. Here, we built a two‐site composite and constructed a conservative age‐depth model to provide a reliable chronology for this rhythmic, highly resolved (<1 kyr) sedimentary archive. Astronomical components (g‐terms and precession constant) are extracted from proxy time‐series using two different techniques, producing consistent results. We find astronomical frequencies up to 4% lower than reported in astronomical solution La04. This solution, however, was smoothed over 20‐Myr intervals, and our results therefore provide constraints on g‐term variability on shorter, million‐year timescales. We also report first evidence that theg4g3“grand eccentricity cycle” may have had a 1.2‐Myr period around 41 Ma, contrary to its 2.4‐Myr periodicity today. Our median precession constant estimate (51.28 ± 0.56″/year) confirms earlier indicators of a relatively low rate of tidal dissipation in the Paleogene. Newfoundland Ridge drift sediments thus enable a reliable reconstruction of astronomical components at the limit of validity of current astronomical calculations, extracted from geologic data, providing a new target for the next generation of astronomical calculations.

     
    more » « less
  3. 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, plusgsmode interactions in the inner solar system. Thegsmode 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
  4. Abstract

    Reconstructions of the paleoclimate indicate that ancient climatic fluctuations on Earth are often correlated with variations in its orbital elements. However, the chaos inherent in the solar system’s orbital evolution prevents numerical simulations from confidently predicting Earth’s past orbital evolution beyond 50–100 Myr. Gravitational interactions among the Sun’s planets and asteroids are believed to set this limiting time horizon, but most prior works approximate the solar system as an isolated system and neglect our surrounding Galaxy. Here we present simulations that include the Sun’s nearby stellar population, and we find that close-passing field stars alter our entire planetary system’s orbital evolution via their gravitational perturbations on the giant planets. This shortens the timespan over which Earth’s orbital evolution can be definitively known by a further ∼10%. In particular, in simulations that include an exceptionally close passage of the Sun-like star HD 7977 2.8 Myr ago, new sequences of Earth’s orbital evolution become possible in epochs before ∼50 Myr ago, which includes the Paleocene–Eocene Thermal Maximum. Thus, simulations predicting Earth’s past orbital evolution before ∼50 Myr ago must consider the additional uncertainty from passing stars, which can open new regimes of past orbital evolution not seen in previous modeling efforts.

     
    more » « less
  5. Abstract

    The Paleocene‐Eocene Thermal Maximum (PETM, ∼56 Ma) presents a past analog for future global warming. Previous studies provided evidence for major loss of dissolved oxygen during the PETM, although understanding the degree and distribution of oxygen loss poses challenges. Magnetofossils produced by magnetotactic bacteria are sensitive to redox conditions in sediments and water columns, and have been used to reconstruct paleoredox conditions over a range of geological settings. Here, we present records of well‐preserved magnetofossils from cores along Walvis Ridge in the Southeast Atlantic that span the PETM across a depth transect (∼1,500–3,600 m paleodepth). Hysteresis, isothermal remanent magnetization curves, first‐order reversal curve diagrams, and low‐temperature magnetic measurements document large variations in magnetic properties of magnetofossils, which relate to time and water depth. Abundant magnetofossil grains are present within the studied sediments, and their morphologies change with paleodepth, as shown by transmission electron microscope observations. Magnetofossils from samples within the PETM onset at the deeper sites (∼2,600–3,600 m paleodepth) have lower coercivity values, a higher oxidation degree, and smaller grain sizes than those from shallower sites (∼1,500–1,800 m paleodepth), likely reflecting changes in paleoredox conditions at different paleodepths. We use the magnetofossil records to reconstruct relative changes in dissolved oxygen content at different water depths through the PETM, and suggest that ocean deoxygenation likely expanded downwards in the early stages of the PETM. We thus demonstrate the value of magnetofossil records for paleoenvironmental reconstructions over time and space, particularly for sediments that lack carbonate fossils.

     
    more » « less