skip to main content


Title: Reduced Variations in Earth’s and Mars’ Orbital Inclination and Earth’s Obliquity from 58 to 48 Myr ago due to Solar System Chaos
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
Award ID(s):
2034660 2001022
NSF-PAR ID:
10370236
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astronomical Journal
Volume:
164
Issue:
3
ISSN:
0004-6256
Format(s):
Medium: X Size: Article No. 107
Size(s):
["Article No. 107"]
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract

    A giant-impact origin for the Moon is generally accepted, but many aspects of lunar formation remain poorly understood and debated. Ćuk et al. proposed that an impact that left the Earth–Moon system with high obliquity and angular momentum could explain the Moon’s orbital inclination and isotopic similarity to Earth. In this scenario, instability during the Laplace Plane transition, when the Moon’s orbit transitions from the gravitational influence of Earth’s figure to that of the Sun, would both lower the system’s angular momentum to its present-day value and generate the Moon’s orbital inclination. Recently, Tian & Wisdom discovered new dynamical constraints on the Laplace Plane transition and concluded that the Earth–Moon system could not have evolved from an initial state with high obliquity. Here we demonstrate that the Earth–Moon system with an initially high obliquity can evolve into the present state, and we identify a spin–orbit secular resonance as a key dynamical mechanism in the later stages of the Laplace Plane transition. Some of the simulations by Tian & Wisdom did not encounter this late secular resonance, as their model suppressed obliquity tides and the resulting inclination damping. Our results demonstrate that a giant impact that left Earth with high angular momentum and high obliquity (θ> 61°) is a promising scenario for explaining many properties of the Earth–Moon system, including its angular momentum and obliquity, the geochemistry of Earth and the Moon, and the lunar inclination.

     
    more » « less
  3. 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
  4. Abstract

    Reliable studies of the long-term dynamics of planetary systems require numerical integrators that are accurate and fast. The challenge is often formidable because the chaotic nature of many systems requires relative numerical error bounds at or close to machine precision (∼10−16, double-precision arithmetic); otherwise, numerical chaos may dominate over physical chaos. Currently, the speed/accuracy demands are usually only met by symplectic integrators. For example, the most up-to-date long-term astronomical solutions for the solar system in the past (widely used in, e.g., astrochronology and high-precision geological dating) have been obtained using symplectic integrators. However, the source codes of these integrators are unavailable. Here I present the symplectic integratororbitN(lean version 1.0) with the primary goal of generating accurate and reproducible long-term orbital solutions for near-Keplerian planetary systems (here the solar system) with a dominant massM0. Among other features,orbitN-1.0includesM0’s quadrupole moment, a lunar contribution, and post-Newtonian corrections (1PN) due toM0(fast symplectic implementation). To reduce numerical round-off errors, Kahan compensated summation was implemented. I useorbitNto provide insight into the effect of various processes on the long-term chaos in the solar system. Notably, 1PN corrections have the opposite effect on chaoticity/stability on a 100 Myr versus Gyr timescale. For the current application,orbitNis about as fast as or faster (factor 1.15–2.6) than comparable integrators, depending on hardware.1

    The orbitN source code (C) is available athttp://github.com/rezeebe/orbitN.

     
    more » « less
  5. Beaufort, Luc (Ed.)
    Abstract. The evolution of the Cenozoic cryosphere from unipolar to bipolar over the past 30 million years (Myr) is broadly known. Highly resolved records of carbonate (CaCO3) content provide insight into the evolution of regional and global climate, cryosphere, and carbon cycle dynamics. Here, we generate the first Southeast Atlantic CaCO3 content record spanning the last 30 Myr, derived from X-ray fluorescence (XRF) ln(Ca/Fe) data collected at Ocean Drilling Program Site 1264 (Walvis Ridge, SE Atlantic Ocean). We present a comprehensive and continuous depth and age model for the entirety of Site 1264 (~316 m; 30 Myr). This constitutes a key reference framework for future palaeoclimatic and palaeoceanographic studies at this location. We identify three phases with distinctly different orbital controls on Southeast Atlantic CaCO3 deposition, corresponding to major developments in climate, the cryosphere and the carbon cycle: (1) strong ~110 kyr eccentricity pacing prevails during Oligocene–Miocene global warmth (~30–13 Ma), (2) increased eccentricity-modulated precession pacing appears after the middle Miocene ClimateTransition (mMCT) (~14–8 Ma), and (3) pervasive obliquity pacing appears in the late Miocene (~7.7–3.3 Ma) following greater importance of high-latitude processes, such as increased glacial activity and high-latitude cooling. The lowest CaCO3 content (92 %–94 %) occurs between 18.5 and 14.5 Ma, potentially reflecting dissolution caused by widespread early Miocene warmth and preceding Antarctic deglaciation across the Miocene Climatic Optimum (~17–14.5 Ma) by 1.5 Myr. The emergence of precession pacing of CaCO3 deposition at Site 1264 after ~14 Ma could signal a reorganisation of surface and/or deep-water circulation in this region following Antarctic reglaciation at the mMCT. The increased sensitivity to precession at Site 1264 between 14 and 13 Ma is associated with an increase in mass accumulation rates (MARs) and reflects increased regional CaCO3 productivity and/or recurrent influxes of cooler, less corrosive deep waters. The highest carbonate content (%CaCO3) and MARs indicate that the late Miocene–early PlioceneBiogenic Bloom (LMBB) occurs between ~7.8 and 3.3Ma at Site 1264; broadly contemporaneous with the LMBB in the equatorial Pacific Ocean. At Site 1264, the onset of the LMBB roughly coincides with appearance of strong obliquity pacing of %CaCO3, reflecting increased high-latitude forcing. The global expression of the LMBB may reflect increased nutrient input into the global ocean resulting from enhanced aeolian dust and/or glacial/chemical weathering fluxes, due to enhanced glacial activity and increased meridional temperature gradients. Regional variability in the timing and amplitude of the LMBB may be driven by regional differences in cooling, continental aridification and/or changes in ocean circulation in the late Miocene. 
    more » « less