La2010: a new orbital solution for the long-term motion of the Earth

TLDR

The earlier La2004 orbital solution had limitations that the new La2010 aims to overcome. This study introduces a new orbital solution for Earth's motion covering the period from 0 to –250 Myr. The solution is obtained by fitting parameters and initial conditions to a 1‑Myr segment of the high‑accuracy INPOP08 ephemeris and by revising the precession equations so they are no longer averaged over the Earth–Moon orbital motion. La2010 remains accurate for more than 50 Myr into the past or future, but its precision rapidly.

Abstract

We present here a new solution for the astronomical computation of the orbital motion of the Earth spanning from 0 to -250 Myr. The main improvement with respect to the previous numerical solution La2004 (Laskar et al. 2004) is an improved adjustment of the parameters and initial conditions through a fit over 1 Myr to a special version of the high accurate numerical ephemeris INPOP08 (Fienga et al. 2009). The precession equations have also been entirely revised and are no longer averaged over the orbital motion of the Earth and Moon. This new orbital solution is now valid over more than 50 Myr in the past or in the future with proper phases of the eccentricity variations. Due to chaotic behavior, the precision of the solution decreases rapidly beyond this time span, and we discuss the behavior of various solutions beyond 50 Myr. For paleoclimate calibrations, we provide several different solutions that are all compatible with the most precise planetary ephemeris. We have thus reached the time where geological data are now required to discriminate among planetary orbital solutions beyond 50 Myr.