At present, tracking data for planetary missions largely consists of radio observables: range-rate (Doppler), range and angular position (VLBI/\(\Delta \)DOR). Future planetary missions may use Interplanetary Laser Ranging (ILR) as a tracking observable. Two-way ILR will provide range data that are about 2 orders of magnitude more accurate than radio-based range data. ILR does not produce Doppler data, however. In this article, we compare the relative strength of radio Doppler and laser range data for the retrieval of parameters of interest in planetary missions, to clarify and quantify the science case of ILR, with a focus on geodetic observables. We first provide an overview of the near-term attainable quality of ILR, in terms of both the realization of the observable and the models used to process the measurements. Subsequently, we analyse the sensitivity of radio Doppler and laser range measurements in representative mission scenarios for parameters of interest. We use both an analytical approximation and numerical analyses of the relative sensitivity of ILR and radio Doppler observables for more general cases. We show that mm-precise range normal points are feasible for ILR, but mm-level accuracy and stability in the full analysis chain are unlikely to be attained, due to a combination of instrumental and model errors. We find that ILR has the potential for superior performance in observing signatures in the data with a characteristic period of greater than 0.33–1.65 hours (assuming 2–10 mm uncertainty for range and 10 \(\upmu \)m/s at 60 s for Doppler). This indicates that Doppler tracking will typically remain the method of choice for gravity field determination and spacecraft orbit determination in planetary missions. ILR data will be able to supplement the orbiter tracking data used for the estimation of parameters with a once-per-orbit signal. Laser ranging data, however, are shown to have a significant advantage for the retrieval of rotational and tidal characteristics from landers. Similarly, laser ranging data will be superior for the construction of planetary ephemerides and the improvement of solar system tests of gravitation, both for orbiter and for lander missions.

Next-generation planetary tracking methods, such as interplanetary laser ranging (ILR) and same-beam interferometry (SBI) promise an orders-of-magnitude increase in the accuracy of measurements of solar system dynamics. This requires a reconsideration of modelling strategies for the translational and rotational dynamics of natural bodies, to ensure that model errors are well below the measurement uncertainties. The influence of the gravitational interaction of the full mass distributions of celestial bodies, the so-called figure-figure effects, will need to be included for selected future missions. The mathematical formulation of this problem to arbitrary degree is often provided in an elegant and compact manner that is not trivially relatable to the formulation used in space geodesy and ephemeris generation. This complicates the robust implementation of such a model in operational software packages. We formulate the problem in a manner that is directly compatible with the implementation used in typical dynamical modelling codes: in terms of spherical harmonic coefficients and Legendre polynomials. An analytical formulation for the associated variational equations for both translational and rotational motion is derived. We apply our methodology to both Phobos and the KW4 binary asteroid system, to analyze the influence of figure-figure effects during estimation from next-generation tracking data. For the case of Phobos, omitting these effects during estimation results in relative errors of 0.42% and 0.065% for the \(\bar{C}_{20}\) and \(\bar{C}_{22}\) spherical harmonic gravity field coefficients, respectively. These values are below current uncertainties, but orders of magnitude larger than those obtained from past simulations for accurate tracking of a future Phobos lander, showing the need to apply the methodology outlined in this manuscript for selected future missions.

Context. The upcoming JUICE and Europa Clipper missions targeting Jupiter’s Galilean satellites will provide radio science tracking measurements of both spacecraft. Such data are expected to significantly help estimating the moons’ ephemerides and related dynamical parameters (e.g. tidal dissipation parameters). However, the two missions will yield an imbalanced dataset, with no flybys planned at Io, condensed over less than six years. Current ephemerides’ solutions for the Galilean moons, on the other hand, rely on ground-based astrometry collected over more than a century which, while being less accurate, bring very valuable constraints on the long-term dynamics of the system.Aims. An improved solution for the Galilean satellites’ complex dynamics could however be achieved by exploiting the existing synergies between these different observation sets.Methods. To quantify this, we merged simulated radio science data from both JUICE and Europa Clipper spacecraft with existing ground-based astrometric and radar observations, and performed the inversion in different configurations: either adding all available ground observations or individually assessing the contribution of different data subsets. Our discussion specifically focusses on the resulting formal uncertainties in the moons’ states, as well as Io’s and Jupiter’s tidal dissipation parameters.Results. Adding astrometry stabilises the moons’ state solution, especially beyond the missions’ timelines. It furthermore reduces the uncertainties in 1/Q(inverse of the tidal quality factor) by a factor two to four for Jupiter, and about 30–35% for Io. Among all data types, classical astrometry data prior to 1960 proved particularly beneficial. Overall, we also show that ground observations of Io add the most to the solution, confirming that ground observations can fill the lack of radio science data for this specific moon.Conclusions. We obtained a noticeable solution improvement when making use of the complementarity between all different observation sets. The promising results obtained with simulations thus motivate future efforts to achieve a global solution from actual JUICE and Clipper radio science measurements.