Tidal effects in planetary systems are the main driver in the orbital migration of natural satellites. They result from physical processes occurring deep inside celestial bodies, whose effects are rarely observable from surface imaging. For giant planet systems, the tidal migration rate is determined by poorly understood dissipative processes in the planet, and standard theories suggest an orbital expansion rate inversely proportional to the power 11/2 in distance, implying little migration for outer moons such as Saturn's largest moon, Titan. Here, we use two independent measurements obtained with the Cassini spacecraft to measure Titan's orbital expansion rate. We find Titan migrates away from Saturn at 11.3 $\pm$ 2.0 cm/year, corresponding to a tidal quality factor of Saturn of Q $\simeq$ 100, and a migration timescale of roughly 10 Gyr. This rapid orbital expansion suggests Titan formed significantly closer to Saturn and has migrated outward to its current position. Our results for Titan and five other moons agree with the predictions of a resonance locking tidal theory, sustained by excitation of inertial waves inside the planet. The associated tidal expansion is only weakly sensitive to orbital distance, motivating a revision of the evolutionary history of Saturn's moon system. The resonance locking mechanism could operate in other systems such as stellar binaries and exoplanet systems, and it may allow for tidal dissipation to occur at larger orbital separations than previously believed.

Uranus and Neptune, and their diverse satellite and ring systems, represent the least explored environments of our Solar System, and yet may provide the archetype for the most common outcome of planetary formation throughout our galaxy. Ice Giants will be the last remaining class of Solar System planet to have a dedicated orbital explorer, and international efforts are under way to realise such an ambitious mission in the coming decades. In 2019, the European Space Agency released a call for scientific themes for its strategic science planning process for the 2030s and 2040s, known as Voyage 2050. We used this opportunity to review our present-day knowledge of the Uranus and Neptune systems, producing a revised and updated set of scientific questions and motivations for their exploration. This review article describes how such a mission could explore their origins, ice-rich interiors, dynamic atmospheres, unique magnetospheres, and myriad icy satellites, to address questions at the heart of modern planetary science. These two worlds are superb examples of how planets with shared origins can exhibit remarkably different evolutionary paths: Neptune as the archetype for Ice Giants, whereas Uranus may be atypical. Exploring Uranus' natural satellites and Neptune's captured moon Triton could reveal how Ocean Worlds form and remain active, redefining the extent of the habitable zone in our Solar System. For these reasons and more, we advocate that an Ice Giant System explorer should become a strategic cornerstone mission within ESA's Voyage 2050 programme, in partnership with international collaborators, and targeting launch opportunities in the early 2030s.

During its mission in the Saturn system, Cassini performed five close flybys of Dione. During three of them, radio tracking data were collected during the closest approach, allowing estimation of the full degree-2 gravity field by precise spacecraft orbit determination. The gravity field of Dione is dominated by $J_{2}$ and $C_{22}$, for which our best estimates are $J_{2} \times 10^6 = 1496 \pm 11$ and $C_{22} \times 10^6 = 364.8 \pm 1.8$ (unnormalized coefficients, 1-$\sigma$ uncertainty). Their ratio is $J_{2}/C_{22} = 4.102 \pm 0.044$, showing a significative departure (about 17-$\sigma$) from the theoretical value of $10/3$, predicted for a relaxed body in slow, synchronous rotation around a planet. Therefore, it is not possible to retrieve the moment of inertia directly from the measured gravitational field. The interior structure of Dione is investigated by a combined analysis of its gravity and topography, which exhibits an even larger deviation from hydrostatic equilibrium, suggesting some degree of compensation. The gravity of Dione is far from the expectation for an undifferentiated hydrostatic body, so we built a series of three-layer models, and considered both Airy and Pratt compensation mechanisms. The interpretation is non-unique, but Dione's excess topography may suggest some degree of Airy-type isostasy, meaning that the outer ice shell is underlain by a higher density, lower viscosity layer, such as a subsurface liquid water ocean. The data permit a broad range of possibilities, but the best fitting models tend towards large shell thicknesses and small ocean thicknesses.

The core of every orbit determination process is the comparison between the measured observables and their predicted values, computed using the adopted mathematical models, and the minimization, in a least square sense, of their differences, known as residuals. In interplanetary orbit determination, Doppler observables, obtained by measuring the average frequency shift of the received carrier signal over a certain count time, are compared against their predicted values, usually computed by differencing two round-trip light-times. This formulation is known to be sensitive to round-off errors, caused by the use of finite arithmetic in the computation, giving rise to an additional noise in the residuals, called numerical noise, that degrades the accuracy of the orbit determination solution. This paper presents a mathematical model for the expected numerical errors in two-way and three-way Doppler observables, computed using the differenced light-time formulation. The model was validated by comparing its prediction to the actual noise in the computed observables, obtained by NASA/Jet Propulsion Laboratory's Orbit Determination Program. The model proved to be accurate within $3 \times 10^{-3} \,\text{mm/s}$ at $60 \,\text{s}$ integration time. Then it was applied to the case studies of Cassini's and Juno's nominal trajectories, proving that numerical errors can assume values up to $6 \times 10^{-2} \,\text{mm/s}$ at $60 \,\text{s}$ integration time, and consequently that they are an important noise source in the Doppler-based orbit determination processes. Three alternative strategies are proposed and discussed in the paper to mitigate the effects of numerical noise.

Ground-based astro-geodetic observations and atmospheric occultations, are two examples of observational techniques requiring a scrutiny analysis of atmospheric refraction. In both cases, the measured changes in observables are geometrically related to changes in the photon path and the light time of the received electromagnetic signal. In the context of geometrical optics, the change in the physical properties of the signal are related to the refractive profile of the crossed medium. Therefore, having a clear knowledge of how the refractivity governs the photon path and the light time evolution is of prime importance to clearly understand observational features. Analytical studies usually focused on spherically symmetric atmospheres and only few aimed at exploring the effect of the non-spherical symmetry on the observables. In this paper, we analytically perform the integration of the photon path and the light time of rays traveling across a planetary atmosphere. We do not restrict our attention to spherically symmetric atmospheres and introduce a comprehensive mathematical framework which allows to handle any kind of analytical studies in the context of geometrical optics. To highlight the capabilities of this new formalism, we carry out five realistic applications for which we derive analytical solutions. The accuracy of the method of integration is assessed by comparing our results to a numerical integration of the equations of geometrical optics in the presence of a quadrupolar moment $J_2$. This shows that the analytical solution leads to the determination of the light time and the refractive bending with relative errors at the level of one part in $10^8$ and one part in $10^5$, for typical values of the refractivity and the $J_2$ parameter at levels of $10^{-4}$ and $10^{-2}$, respectively.

Radio remote sensing of the heliosphere using spacecraft radio signals has been used to study the near-sun plasma in and out of the ecliptic, close to the sun, and on spatial and temporal scales not accessible with other techniques. Studies of space-time variations in the inner solar wind are particularly timely because of the desire to understand and predict space weather, which can disturb satellites and systems at 1AU and affect human space exploration. Here we demonstrate proof-of-concept of a new radio science application for spacecraft radio science links. The differing transfer functions of plasma irregularities to spacecraft radio up- and downlinks can be exploited to localize plasma scattering along the line of sight. We demonstrate the utility of this idea using Cassini radio data taken in 2001-2002. Under favorable circumstances we demonstrate how this technique, unlike other remote sensing methods, can determine center-of-scattering position to within a few thousandths of an AU and thickness of scattering region to less than about 0.02 AU. This method, applied to large data sets and used in conjunction with other solar remote sensing data such as white light data, has space weather application in studies of inhomogeneity and nonstationarity in the near-sun solar wind.