Context. Upcoming science missions to Phobos will potentially provide unprecedented observations of Phobos’s orbit in the form of orbiter and/or lander tracking data. This will likely require an updating of the dynamical models currently used to invert this data, with the coupling between the satellite’s orbit and rotation being of particular importance. State-of-the-art ephemerides estimations for tidally locked satellites rely on a decoupled approach where translational models are combined with a simplified analytical representation of the moon’s rotation (typically a single-frequency periodic variation superimposed to a synchronous rotation). Aims. This paper investigates the coupled propagation of Phobos’s translational and rotational dynamics, and assesses the extent to which the most commonly used uncoupled model can emulate the results of the coupled integration, and what consequences the mismodeling has on the products of data inversion. Methods. We considered two models: a coupled model that propagates Phobos’s translational and rotational dynamics simultaneously, and an uncoupled model that assumes Phobos to be in a fully locked configuration with a once-per-orbit longitudinal libration. By simulating the dynamics for about ten years, first in a coupled and then in an uncoupled manner, we compared the results and used the coupled trajectory as simulated observations for an estimation of the different parameters using uncoupled translational dynamics. Results. For identical initial states, differences between the coupled and uncoupled trajectories were found to accumulate to 40 m, most predominantly in Phobos’s direction of motion. Longitudinal librations were misrepresented by the uncoupled model particularly around the frequencies of the normal mode, where forcings are amplified up to 3.6 × 10⁻³ degrees. Long-term latitudinal librations also arise from forcings due to coupling-induced changes in orbital inclination. The use of uncoupled models in data inversion results in true errors in the estimated parameters. In this case, we performed estimations of different lengths up to 1000 days to estimate Phobos’s initial state, once-per-orbit libration amplitude, and harmonic coefficients C_2,0 and C_2,2. Errors in dynamical parameters were found to be on the order of 10⁻³ degrees for the physical libration amplitude and of 10⁻⁵ for the harmonic coefficients (relative errors of around 0.1%). Conclusions. These true errors are one to three orders of magnitude above the formal errors expected from laser ranging measurements to a Phobos lander, which indicates that the typical single-frequency uncoupled model is not suitable for the proper inversion of such data. Refined rotation models will therefore be required, either by expanding the uncoupled model to multiple frequencies or by performing a fully coupled orbital-rotational propagation as proposed in this paper. We discuss the theoretical and practical limitations of an extended analytical parametrization in the specific case of tidally locked satellites, and advocate for the use of a fully coupled approach.

In 2031 the JUICE spacecraft will perform a multi-flyby tour of the Jovian system. Next to the radiometric tracking that is performed for navigation operations, the dedicated radio science instrument (3GM) collects high-accuracy radiometric measurements during the flybys. We investigate the capability of the radio science data to provide improved moon state knowledge during navigational operations. We introduce ephemeris updates from radio science data into our simulated navigation operations and examine the potential savings of statistical ΔV for corrective manoeuvres. A navigation orbit determination (OD) solution was simulated for the multi-flyby tour of JUICE, including the resulting state knowledge evolution of the Galilean moons. The OD was extended by an interface for external moon ephemeris updates, which was used to evaluate the impact of radio science generated external ephemerides on the statistical ΔV budgets for post-flyby cleanup manoeuvres. The moon state knowledge evolution during navigation operation showed a rapid reduction of the a-priori moon state uncertainty, for which the navigational tracking data coverage of the long, early tour arcs was identified as the driving factor. As a result of the longer tracking arcs, the moon state knowledge from navigation data results improves more quickly during the initial phase of the tour. Since the impact of moon state knowledge on the corrective manoeuvres is largest in this initial phase, the comparative analysis of the statistical ΔV cost shows that the adoption of radio science ephemeris products does not effectuate significant ΔV savings. Instead we showed that in order to achieve substantial ΔV savings improvements of Europa's and Ganymede's ephemerides are required ahead of JUICE's arrival. While the analysis concludes that data synergies are unlikely to benefit the navigational operations, it highlights other potential synergies between the navigation and radio science data. A comparatively strong signature of Io's dynamics was found in the simulated navigation data along the long early tour arcs, which could be leveraged for the benefit of the new global moon ephemeris solutions after JUICE.

Juice (JUpiter ICy moons Explorer) 3GM Radio Science Experiment will map the gravitational field of Ganymede with unprecedented accuracy and measure tidally-induced variations. These measurements will allow the characterization of its putative ocean and may resolve lateral variations in internal structure. Lateral variations cause an additional tidal signal that depends on their wavelength and amplitude. We show that shell thickness variations of (Formula presented.) the mean thickness produce an additional tidal signal (Formula presented.) times smaller than the main tidal signal, detectable given the accuracy of Juice. Using a Bayesian framework, we show that measuring differences between (Formula presented.) and (Formula presented.) constrain equator to pole shell thickness differences. Also measuring the degree-3 spherical harmonic signal due to degree 2 forcing constrain degree-1 and degree-3 structure. This demonstrates tidal tomography's potential to map three dimensional structure and supports its consideration for future missions.

We present an overview of the operations and engineering interface for Planetary Radio Interferometry and Doppler Experiment (PRIDE) radio astronomy observations as a scientific component of the ESA’s Jupiter Icy Moons Explorer (JUICE) mission, as well as other prospective planetary and space science missions. The article discusses advanced scheduling and planning methods that make it possible to create observing schedules for observations of specific spacecraft in concurrence with observations of natural radio sources. In order to put this into practice and find suitable natural background calibrator sources for PRIDE of JUICE mission, we developed planning and scheduling software. The conventional scheduling software for natural celestial radio sources is not set up to include spacecraft as observation targets in the necessary control files. Therefore, difficulties already arise during observation planning. We report on the development of new and the adaptation of existing routines used in astrophysical and geodetic VLBI for satellite scheduling and planning. The analysis of the PRIDE science observations led to improved observational planning, and the mission’s scheduling methodologies were studied using a systems engineering approach. In addition, we highlighted the new procedures, like finding charts for selecting calibrator radio sources over a range of frequency bands and the outcomes of those strategies for science operation planning. A simulation of the flyby of Venus during the cruise phase of the JUICE spacecraft, based on the Tudat software, is also presented, resulting in a promising opportunity to test PRIDE techniques and evaluate the improvements that PRIDE observables can make to natural bodies’ ephemerides. The first K _a-band (32 GHz) observations of the ESA’s BepiColombo by a radio telescope in the VLBI network, which employs a similar radio communications system as JUICE, were also demonstrated as a test case. The primary objective of these activities is to serve as a practice run for the upcoming operational PRIDE JUICE operations. We showcase the capabilities of the planning and scheduling software for other space missions.

In the coming decade, JUICE and Europa Clipper radio-science will yield the most accurate estimation to date of the Galilean moons’ physical parameters and ephemerides. JUICE's PRIDE (Planetary Radio Interferometry and Doppler Experiment) will help achieve such a solution by providing VLBI (Very Long Baseline Interferometry) observations of the spacecraft's lateral position, complementing nominal radio-science measurements. In this paper, we quantify how PRIDE VLBI can contribute to the moons’ ephemerides determination, in terms of attainable solution improvement and validation opportunities. To this end, we simulated VLBI data for JUICE, but also investigated the possibility to perform simultaneous tracking of JUICE and Europa Clipper, thus ultimately generating both single- and dual-spacecraft VLBI. We considered various tracking and data quality scenarios for both VLBI types, and compared the formal uncertainties provided by covariance analyses with and without VLBI. These analyses were performed for both global and local (i.e. per-flyby) estimations of the moons’ states, as eventually achieving a global solution first requires proceeding arc-per-arc. We showed that both single- and multi-spacecraft VLBI measurements only bring limited improvement to the global state estimation, but significantly contribute to the moons’ normal points (i.e. local states at flyby times), most notably in the out-of-plane direction. Additionally, we designed a validation plan exploiting PRIDE VLBI to progressively validate the classical radio-science solution, whose robustness and statistical realism is sensitive to modelling inconsistencies. By improving the local state estimations and offering various validation opportunities, PRIDE will be invaluable in overcoming possible dynamical challenges. It can therefore play a key role in reconstructing a global solution for the Galilean moons’ dynamics with the uncertainty levels promised by JUICE-Europa Clipper analyses. This, in turn, is critical to the accurate characterisation of tidal dissipation in the Jovian system, holding the key to the long-term evolution of the Galilean moons.

We investigate the impact of viscoelastic tidal deformation of the Moon on the motion of a polar orbiter. The dissipative effects in the Moon’s interior, i.e., tidal phase lags, are modeled as Fourier series sampled at given frequencies associated with linear combinations of Delaunay arguments, the fundamental parameters describing the lunar motion around the Earth and the Sun. We implement the tidal model to evaluate the temporal lunar gravity field and the induced perturbation on the orbiter. We validate the numerical scheme via a frequency analysis of the perturbed orbital motion. We show that, in the case of the Lunar Reconnaissance Orbiter at a low altitude of less than 200 km, the main lunar tides and hence the potential Love numbers around the monthly and some multiple frequencies are dynamically separable. The omission of those effects in practice introduces a position error at the level of a few decimeters within 10 days.

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.

Context. Stellar occultations currently provide the most accurate ground-based measurements of the positions of natural satellites (down to a few kilometres for the Galilean moons). However, when using these observations in the calculation of satellite ephemerides, the uncertainty in the planetary ephemerides dominates the error budget of the occultation. Aims. We quantify the local refinement in the central planet's position achievable by performing Very Long Baseline Interferometry (VLBI) tracking of an in-system spacecraft temporally close to an occultation. We demonstrate the potential of using VLBI to enhance the science return of stellar occultations for satellite ephemerides. Methods. We identified the most promising observation and tracking opportunities offered by the Juno spacecraft around Jupiter as perfect test cases, for which we ran simulations of our VLBI experiment. Results. VLBI tracking at Juno's perijove close to a stellar occultation locally (in time) reduces the uncertainty in Jupiter's angular position in the sky to 250-400 m. This represents up to an order of magnitude improvement with respect to current solutions and is lower than the stellar occultation error, thus allowing the moon ephemeris solution to fully benefit from the observation. Conclusions. Our simulations showed that the proposed tracking and observation experiment can efficiently use synergies between ground- and space-based observations to enhance the science return on both ends. The reduced error budget for stellar occultations indeed helps to improve the moons'ephemerides, which in turn benefit planetary missions and their science products, such as the recently launched JUICE and upcoming Europa Clipper missions.

When reconstructing natural satellites' ephemerides from space missions' tracking data, the dynamics of the spacecraft and natural bodies are often solved for separately, in a decoupled manner. Alternatively, the ephemeris generation and spacecraft orbit determination can be performed concurrently. This method directly maps the available data set to the estimated parameters' covariances while fully accounting for all dynamical couplings. It thus provides a statistically consistent solution to the estimation problem, whereas this is not directly ensured with the decoupled strategy. For the Galilean moons in particular, the JUICE mission provides a unique, although challenging, opportunity for ephemerides improvement. For such a dynamically coupled problem, choosing between the two state estimation strategies will be influential. This paper quantifies the Galilean moons' state uncertainties attainable when applying a coupled estimation strategy to simulated JUICE data, and discusses the challenges that remain to be addressed to achieve such a coupled solution from real observations. We first provide a detailed, explicit formulation for the coupled approach, which was still missing in the literature although already used in past studies. We then assessed the relative performances of the two ephemerides generation techniques for the JUICE test case. To this end, we used both decoupled and coupled models on simulated JUICE radiometric data. We compared the resulting covariances for the Galilean moons' states, and showed that the decoupled approach yields slightly lower formal errors for the moons' tangential positions. However, the coupled model can reduce the state uncertainties by more than one order of magnitude in the radial direction (i.e. towards the central body). It also proved more sensitive to the dynamical coupling between Io, Europa and Ganymede, allowing the state solutions for the first two moons to fully benefit from JUICE orbital phase around Ganymede. On the other hand, we showed that the choice of state estimation methods does not strongly affect the moons' gravity field determination. Many issues still remain to be solved before a concurrent estimation strategy can be successfully applied, especially to reconstruct the moons’ dynamics over long timescales. Nonetheless, our analysis highlights promising ephemerides improvements and thus motivates future efforts to reach a coupled state solution for the Galilean moons.

Photochemical aerosols were detected as high as 350 km of altitude in Pluto's atmosphere during the New Horizons fly-by. These aerosols are thought to affect Pluto's climate, by acting as cooling agents, and the colours of Pluto's surface, in particular in the dark regions named Cthulhu and Krun and at the North Pole. Pluto atmospheric and surface models have so far used the optical constants of Titan aerosol analogues (tholins), whereas their chemical composition is known to differ from that of Pluto aerosol analogues. In order to provide a new set of optical constants for Pluto tholins, we synthesized analogues of Pluto's aerosols and determined with spectroscopic ellipsometry their optical constants from 270 to 2100 nm. Three types of samples were produced from N₂:CH₄:CO gas mixtures differing in their CH₄:N₂ mixing ratio, representative of different altitudes in Pluto's current atmosphere or different seasons or epochs of Pluto. Our analysis shows a strong absorption by Pluto tholins in the UV and visible spectral ranges, with k index of a few 10⁻¹ at 270 nm, in agreement with N- and O-bearing organic molecules. Pluto tholins are less absorbent in the near-IR than in the UV–Vis wavelength range, with k of a few 10⁻³ between 600 and 2100 nm. Our comparative study highlights the dependency of n and k indices to the CH₄:N₂ mixing ratio. Aerosols formed at different altitudes in Pluto's atmosphere or during different seasons or epochs of Pluto will therefore affect the budget of Pluto radiative transfer differently. The optical constants presented in this study were tested with a Pluto surface model and with a model of light scattering. The surface modelling results highlight the suitability of these optical constants to reproduce Pluto compositional observations in the visible spectral range by MVIC and LEISA. The atmospheric modelling results conclude that Pluto tholins absorb 5 to 10 times less than Titan tholins at 500 nm, and this lower absorption is consistent with Alice observations of Pluto's haze.

Pluto's fly-by by the New Horizons spacecraft in July 2015 has revealed a dark reddish equatorial region, named Cthulhu, covered by a dark, non-icy material whose origin and composition have yet to be determined. It has been suggested that this material could form from the sedimentation of photochemical aerosols, originating from dissociation and ionisation processes in Pluto's high atmosphere (similarly to aerosols forming Titan's haze). This hypothesis is here further investigated by comparing New Horizons spectra collected both in the visible and the near-infrared to laboratory reflectance measurements of analogues of Pluto's aerosols (Pluto tholins). These aerosols were synthesised in conditions mimicking Pluto's atmosphere, and their optical and reflectance properties were determined, before being used in Hapke models. In particular, the single scattering albedo and phase function of Pluto tholins were retrieved through Hapke model inversion, performed from laboratory reflectance spectra collected under various geometries. From reconstructed reflectance spectra and direct comparison with New Horizons data, some of these tholins are shown to reproduce the photometric level (i.e. reflectance continuum) reasonably well in the near-infrared. Nevertheless, a misfit of the red visible slope still remains and tholins absorption bands present in the modelled spectra are absent in those collected by the New Horizons instruments. Several hypotheses are considered to explain the absence of these absorption features in LEISA data, namely high porosity effects or GCR irradiation. The formation of highly porous structures, which is currently our preferred scenario, could be promoted by either sublimation of ices initially mixed with the aerosols, or gentle deposition under Pluto's weak gravity.

Context. The apparent close encounters of two satellites in the plane of the sky, called mutual approximations, have been suggested as a different type of astrometric observation to refine the moons' ephemerides. The main observables are then the central instants of the close encounters, which have the advantage of being free of any scaling and orientation errors. However, no analytical formulation is available yet for the observation partials of these central instants, leaving numerical approaches or alternative observables (i.e. derivatives of the apparent distance instead of central instants) as options. Aims. Filling that gap, this paper develops an analytical method to include central instants as direct observables in the ephemerides estimation and assesses the quality of the resulting solution. Methods. To this end, the apparent relative position between the two satellites is approximated by a second-order polynomial near the close encounter. This eventually leads to an expression for mutual approximations' central instants as a function of the apparent relative position, velocity, and acceleration between the two satellites. Results. The resulting analytical expressions for the central instant partials were validated numerically. In addition, we ran a covariance analysis to compare the estimated solutions obtained with the two types of observables (central instants versus alternative observables), using the Galilean moons of Jupiter as a test case. Our analysis shows that alternative observables are almost equivalent to central instants in most cases. Accurate individual weighting of each alternative observable, accounting for the mutual approximation's characteristics (which are automatically included in the central instants' definition), is however crucial to obtain consistent solutions between the two observable types. Using central instants still yields a small improvement of 10-20% of the formal errors in the radial and normal directions (RSW frame), compared to the alternative observables' solution. This improvement increases when mutual approximations with low impact parameters and large impact velocities are included in the estimation. Conclusions. Choosing between the two observables thus requires careful assessment, taking into account the characteristics of the available observations. Using central instants over alternative observables ensures that the state estimation fully benefits from the information encoded in mutual approximations, which might be necessary depending on the application of the ephemeris solution.