Abstract. We present here a new solution for the astronomical computation of the insolation quantities on Earth spanning from −250 Myr to 250 Myr. This solution has been improved with respect to La93 (Laskar et al. 1993) by using a direct integration of the gravitational equations for the orbital motion, and by improving the dissipative contributions, in particular in the evolution of the Earth-Moon System. The orbital solution has been used for the calibration of the Neogene period (Lourens et al. 2004), and is expected to be used for age calibrations of paleoclimatic data over 40 to 50 Myr, eventually over the full Palaeogene period (65 Myr) with caution. Beyond this time span, the chaotic evolution of the orbits prevents a precise determination of the Earth's motion. However, the most regular components of the orbital solution could still be used over a much longer time span, which is why we provide here the solution over 250 Myr. Over this time interval, the most striking feature of the obliquity solution, apart from a secular global increase due to tidal dissipation, is a strong decrease of about 0.38 degree in the next few millions of years, due to the crossing of the s 6 + g 5 − g 6 resonance (Laskar et al. 1993). For the calibration of the Mesozoic time scale (about 65 to 250 Myr), we propose to use the term of largest amplitude in the eccentricity, related to g 2 − g 5 , with a fixed frequency of 3.200 /yr, corresponding to a period of 405 000 yr. The uncertainty of this time scale over 100 Myr should be about 0.1%, and 0.2% over the full Mesozoic era.

As the obliquity of Mars is strongly chaotic, it is not possible to give a solution for its evolution over more than a few million years. Using the most recent data for the rotational state of Mars, and a new numerical integration of the Solar System, we provide here a precise solution for the evolution of Mars' spin over 10 to 20 Myr. Over 250 Myr, we present a statistical study of its possible evolution, when considering the uncertainties in the present rotational state. Over much longer time span, reaching 5 Gyr, chaotic diffusion prevails, and we have performed an extensive statistical analysis of the orbital and rotational evolution of Mars, relying on Laskar's secular solution of the Solar System, based on more than 600 orbital and 200 000 obliquity solutions over 5 Gyr. The density functions of the eccentricity and obliquity are specified with simple analytical formulas. We found an averaged eccentricity of Mars over 5 Gyr of 0.0690 with standard deviation 0.0299, while the averaged value of the obliquity is 37.62 • with a standard deviation of 13.82 • , and a maximal value of 82.035 •. We find that the probability for Mars' obliquity to have reached more than 60 • in the past 1 Gyr is 63.0%, and 89.3% in 3 Gyr. Over 4 Gyr, the position of Mars' axis is given by a uniform distribution on a spherical cap limited by the obliquity 58.62 • , with the addition of a random noise allowing a slow diffusion beyond this limit. We can also define a standard model of Mars' insolation parameters over 4 Gyr with the most probable values 0.068 for the eccentricity and 41.80 • for the obliquity.

No abstract

No abstract

No abstract

No abstract

We study the global dynamics of the jovian Trojan asteroids by means of the frequency map analysis. We find and classify the main resonant structures that serve as skeleton of the phase space near the Lagrangian points. These resonances organize and control the long‐term dynamics of the Trojans. Besides the secondary and secular resonances, that have already been found in other asteroid sets in mean motion resonance (e.g. main belt, Kuiper belt), we identify a new type of resonance that involves secular frequencies and the frequency of the great inequality, but not the libration frequency. Moreover, this new family of resonances plays an important role in the slow transport mechanism that drives Trojans from the inner stable region to eventual ejections. Finally, we relate this global view of the dynamics with the observed Trojans, identify the asteroids that are close to these resonances and study their long‐term behaviour.

We develop an analytical Hamiltonian formalism adapted to the study of the motion of two planets in co-orbital resonance. The Hamiltonian, averaged over one of the planetary mean longitude, is expanded in power series of eccentricities and inclinations. The model, which is valid in the entire co-orbital region, possesses an integrable approximation modeling the planar and quasi-circular motions. First, focusing on the fixed points of this approximation, we highlight relations linking the eigenvectors of the associated linearized differential system and the existence of certain remarkable orbits like the elliptic Eulerian Lagrangian configurations, the Anti-Lagrange (Giuppone et al., 2010) orbits and some second sort orbits discovered by Poincar\'e. Then, the variational equation is studied in the vicinity of any quasi-circular periodic solution. The fundamental frequencies of the trajectory are deduced and possible occurrence of low order resonances are discussed. Finally, with the help of the construction of a Birkhoff normal form, we prove that the elliptic Lagrangian equilateral configurations and the Anti-Lagrange orbits bifurcate from the same fixed point L4.Comment: 25 pages. Accepted for publication in Celestial Mechanics and Dynamical Astronom