Source author record

Richard E. Zeebe

Richard E. Zeebe appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

4works
3topics
2close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

4 published item(s)

preprint2022arXiv

Reduced variations in Earth's and Mars' orbital inclination and Earth's obliquity from 58 to 48 Myr ago due to solar system chaos

The dynamical evolution of the solar system is chaotic with a Lyapunov time of only $\sim$5 Myr for the inner planets. Due to the chaos it is fundamentally impossible to accurately predict the solar system's orbital evolution beyond $\sim$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 $\sim$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 $\sim$58 to $\sim$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 ($g$- and $s$-modes) 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 of $s$-modes, plus $g$-$s$-mode interactions in the inner solar system. The $g$-$s$-mode 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.

preprint2021arXiv

Stepsize errors in the $N$-body problem: discerning Mercury's true possible long-term orbits

Numerical integrations of the Solar System have been carried out for decades. Their results have been used, for example, to determine whether the Solar System is chaotic, whether Mercury's orbit is stable, or to help discern Earth's climate history. We argue that all of the past studies we consider in this work are affected by numerical chaos to different degrees, affecting the possible orbits and instability probability of Mercury, sometimes significantly. We show how to eliminate the effects of numerical chaos by resolving Mercury's pericentre passage. We also show that several higher order symplectic maps do not exhibit significant differences in resolving pericentre passage of Mercury (at fixed time step), making their advantages suspect for calculating long-term orbits. Resolving pericentre passage affects a wide array of orbital numerical studies, like exoplanet studies, studies of the galactic centre, and other $N$-body problems.

preprint2015arXiv

Dynamic stability of the Solar System: Statistically inconclusive results from ensemble integrations

Due to the chaotic nature of the Solar System, the question of its long-term stability can only be answered in a statistical sense, for instance, based on numerical ensemble integrations of nearby orbits. Destabilization of the inner planets, leading to close encounters and/or collisions can be initiated through a large increase in Mercury's eccentricity, with a currently assumed likelihood of ~1%. However, little is known at present about the robustness of this number. Here I report ensemble integrations of the full equations of motion of the eight planets and Pluto over 5 Gyr, including contributions from general relativity. The results show that different numerical algorithms lead to statistically different results for the evolution of Mercury's eccentricity (eM). For instance, starting at present initial conditions (eM ~= 0.21), Mercury's maximum eccentricity achieved over 5 Gyr is on average significantly higher in symplectic ensemble integrations using heliocentricthan Jacobi coordinates and stricter error control. In contrast, starting at a possible future configuration (eM ~= 0.53), Mercury's maximum eccentricity achieved over the subsequent 500 Myr is on average significantly lower using heliocentric than Jacobi coordinates. For example, the probability for eM to increase beyond 0.53 over 500 Myr is >90% (Jacobi) vs. only 40-55% (heliocentric). This poses a dilemma as the physical evolution of the real system - and its probabilistic behavior - cannot depend on the coordinate system or numerical algorithm chosen to describe it. Some tests of the numerical algorithms suggest that symplectic integrators using heliocentric coordinates underestimate the odds for destabilization of Mercury's orbit at high initial eM.

preprint2015arXiv

Highly stable evolution of Earth's future orbit despite chaotic behavior of the Solar System

Due to the chaotic nature of the Solar System, the question of its dynamic long-term stability can only be answered in a statistical sense, e.g. based on numerical ensemble integrations of nearby orbits. Destabilization, including catastrophic encounters and/or collisions involving the Earth, has been suggested to be initiated through a large increase in Mercury's eccentricity (eM), with an estimated probability of ~1%. However, it has recently been shown that the statistics of numerical Solar System integrations are sensitive to the accuracy and type of numerical algorithm. Here I report results from computationally demanding ensemble integrations (N=1,600 with slightly different initial conditions) at unprecedented accuracy based on the full equations of motion of the eight planets and Pluto over 5Gyr, including contributions from general relativity. The standard symplectic algorithm produced spurious results for highly eccentric orbits and during close encounters, which were hence integrated with a suitable Bulirsch-Stoer algorithm, specifically designed for these situations. The present study yields odds for a large increase in Mercury's eccentricity that are less than previous estimates. Strikingly, in two solutions Mercury continued on highly eccentric orbits (after reaching eM values >0.93) for 80-100Myr before colliding with Venus or the Sun. Most importantly, none of the 1,600 solutions led to a close encounter involving the Earth or a destabilization of Earth's orbit in the future. I conclude that Earth's orbit is dynamically highly stable for billions of years, despite the chaotic behavior of the Solar System.