Context. The strong tidal dissipation in the couple Jupiter-Io is spread to all the moons involved in the Laplace resonance (Io, Europa, and Ganymede), leading to a migration of their orbits. Aims. We aim to characterize the future behavior of the Galilean satellites over the Solar System lifetime and to quantify the stability of the Laplace resonance. Since tidal dissipation makes possible the exit from the current resonances or capture into new ones, we investigate the capture of Callisto into resonance. Methods. We perform hundreds of propagations using an improved version of a recent semi-analytical model. As Ganymede moves outwards, it approaches the 2:1 resonance with Callisto, inducing a temporary chaotic motion in the system. For this reason, we draw a statistical picture of the outcome of the resonant encounter. Results. The system can settle into two distinct outcomes: A) a chain of three 2:1 two-body resonances (Io-Europa, Europa-Ganymede and Ganymede-Callisto), or B) a resonant chain involving the 2:1 two-body resonance Io-Europa plus at least one pure 4:2:1 three-body resonance, most frequently between Europa, Ganymede and Callisto. In case A (56\% of the simulations), the Laplace resonance is always preserved and the eccentricities remain confined to small values below 0.01. In case B (44\% of the simulations), the Laplace resonance is generally disrupted and the eccentricities of Ganymede and Callisto can increase up to about 0.1, making this configuration unstable and driving the system into new resonances. Conclusion. From our results, the capture of Callisto into resonance appears to be extremely likely (100\% of our simulations). Assuming the most recent estimate of the dissipation between Io and Jupiter, the resonant encounter happens at about 1.5 Gyrs from now. Therefore, the stability of the Laplace resonance is guaranteed at least up to about 1.5 Gyrs.

The Galilean satellites' dynamics has been studied extensively during the last century. In the past it was common to use analytical expansions in order to get simple models to integrate, but with the new generation computers it became prevalent the numerical integration of very sophisticated and almost complete equations of motion. In this article we aim to describe the resonant and secular motion of the Galilean satellites through a Hamiltonian, depending on the slow angles only, obtained with an analytical expansion of the perturbing functions and an averaging operation. In order to have a model as near as possible to the actual dynamics, we added perturbations and we considered terms that in similar studies of the past were neglected, such as the terms involving the inclinations and the Sun's perturbation. Moreover, we added the tidal dissipation into the equations, in order to investigate how well the model captures the evolution of the system.

The JUICE mission will perform detailed measurements of the properties of the Galilean moons, with a nominal in-system science-mission duration of about 3.5 years. Using both the radio tracking data, and (Earth- and JUICE-based) optical astrometry, the dynamics of the Galilean moons will be measured to unprecedented accuracy. This will provide crucial input to the determination of the ephemerides and physical properties of the system, most notably the dissipation in Io and Jupiter. The data from Planetary Radio Interferometry and Doppler Experiment (PRIDE) will provide the lateral position of the spacecraft in the ICRF. In this article, we analyze the relative quantitative influence of the JUICE-PRIDE observables to the determination of the ephemerides of the Jovian system and the associated physical parameters. We perform a covariance analysis for a broad range of mission and system characteristics. We analyze the influence of VLBI data quality, observation planning, as well as the influence of JUICE orbit determination quality. This provides key input for the further development of the PRIDE observational planning and ground segment development. Our analysis indicates that the VLBI data are especially important for constraining the dynamics of Ganymede and Callisto perpendicular to their orbital planes. Also, the use of the VLBI data makes the uncertainty in the ephemerides less dependent on the error in the orbit determination of the JUICE spacecraft itself. Furthermore, we find that optical astrometry data of especially Io using the JANUS instrument will be crucial for stabilizing the solution of the normal equations. Knowledge of the dissipation in the Jupiter system cannot be improved using satellite dynamics obtained from JUICE data alone, the uncertainty in Io's dissipation obtained from our simulations is similar to the present level of uncertainty.

Aims. Numerous trans-Neptunian objects are known to be in mean-motion resonance with Neptune. We aim to describe their long-term orbital evolution (both past and future) by means of a one-degree-of-freedom secular model. In this paper, we focus only on objects with a semi-major axis larger than 50 astronomical units (au). Methods. For each resonant object considered, a 500 000-year numerical integration is performed. The output is digitally filtered to get the parameters of the resonant secular model. Their long-term (Giga-year) orbital evolution is then represented by the level curves of the secular Hamiltonian. Results. For the majority of objects considered, the mean-motion resonance has little impact on the long-term trajectories (the secular dynamics is similar to a non-resonant one). However, a subset of objects is strongly affected by the resonance, producing moderately-high-amplitude oscillations of the perihelion distance and/or libration of the argument of perihelion around a fixed centre. Moreover, the high perihelion distance of the object 2015 FJ345 is plainly explained by long-term resonant dynamics, allowing us to also deduce its orbital elements at the time of capture in resonance (at least 15 million years ago). The same type of past evolution is expected for 2014 FZ71.