Element history of the Laplace resonance: A dynamical approach

Context. We consider the three-body mean motion resonance defined by the Jovian moons Io, Europa, and Ganymede, which is commonly known as the Laplace resonance. In terms of the moons' mean longitudes lambda(1) (Io), lambda(2) (Europa), and lambda(3) (Ganymede), this resonance is described by the librating argument phi(L) lambda(1) - 3 lambda(2) + 2 lambda(3)( )approximate to 180 degrees, which is the sum of phi(12) lambda(1) - 2 lambda(2) + pi(2) approximate to 180 degrees and phi(23) lambda(2) - 2 lambda(3) pi(2)( )approximate to 0 degrees, where pi(2) denotes Europa's longitude of perijove.Aims. In particular, we construct approximate models for the evolution of the librating argument phi(L) over the period of 100 yr, focusing on its principal amplitude and frequency, and on the observed mean motion combinations n(1) - 2n(2) and n(2) - 2n(3) associated with the quasi-resonant interactions above.Methods. First, we numerically propagated the Cartesian equations of motion of the Jovian system for the period under examination, and by comparing the results with a suitable set of ephemerides, we derived the main dynamical effects on the target quantities. Using these effects, we built an alternative Hamiltonian formulation and used the normal forms theory to precisely locate the resonance and to semi-analytically compute its main amplitude and frequency.Results. From the Cartesian model we observe that on the timescale considered and with ephemerides as initial conditions, both phi(L) and the diagnostics n(1) - 2n(2) and n(2) - 2n(3) are well approximated by considering the mutual gravitational interactions of Jupiter and the Galilean moons (including Callisto), and the effect of Jupiter's J(2) harmonic. Under the same initial conditions, the Hamiltonian formulation in which Callisto and J(2) are reduced to their secular contributions achieves larger errors for the quantities above, particularly for phi(L) . By introducing appropriate resonant variables, we show that these errors can be reduced by moving in a certain action-angle phase plane, which in turn implies the necessity of a tradeoff in the selection of the initial conditions.Conclusions. In addition to being a good starting point for a deeper understanding of the Laplace resonance, the models and methods described are easily generalizable to different types of multi-body mean motion resonances. Thus, they are also prime tools for studying the dynamics of extrasolar systems