Search

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. 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 both single- and dual-spacecraft VLBI, exploiting the potential simultaneous tracking of JUICE and Europa Clipper. We considered various tracking and data quality scenarios, and compared the formal uncertainties obtained with and without VLBI. This was performed for both global and local (i.e., per-flyby) estimations of the moons' states, as achieving a global solution first requires proceeding arc-per-arc. We showed that both single- and dual-spacecraft VLBI 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. By improving the local state estimations and offering various validation opportunities, PRIDE will be invaluable in overcoming possible dynamical modelling 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.

Planetary Radio Interferometry and Doppler Experiment (PRIDE) is a multi-purpose experimental technique aimed at enhancing the science return of planetary missions. The technique exploits the science payload and spacecraft service systems without requiring a dedicated onboard instrumentation or imposing on the existing instrumentation any special for PRIDE requirements. PRIDE is based on the near-field phase-referencing Very Long Baseline Interferometry (VLBI) and evaluation of the Doppler shift of the radio signal transmitted by spacecraft by observing it with multiple Earth-based radio telescopes. The methodology of PRIDE has been developed initially at the Joint Institute for VLBI ERIC (JIVE) for tracking the ESA's Huygens Probe during its descent in the atmosphere of Titan in 2005. From that point on, the technique has been demonstrated for various planetary and other space science missions. The estimates of lateral position of the target spacecraft are done using the phase-referencing VLBI technique. Together with radial Doppler estimates, these observables can be used for a variety of applications, including improving the knowledge of the spacecraft state vector. The PRIDE measurements can be applied to a broad scope of research fields including studies of atmospheres through the use of radio occultations, the improvement of planetary and satellite ephemerides, as well as gravity field parameters and other geodetic properties of interest, and estimations of interplanetary plasma properties. This paper presents the implementation of PRIDE as a component of the ESA's Jupiter Icy Moons Explorer (JUICE) mission., Comment: Accepted for publication in Space Science Reviews, v. 219 (special issue on the ESA's JUICE mission)

Context. The 98{\deg}-obliquity of Uranus is commonly attributed to giant impacts that occurred at the end of the planetary formation. This picture, however, is not devoid of weaknesses. Aims. On a billion-year timescale, the tidal migration of the satellites of Jupiter and Saturn has been shown to strongly affect their spin-axis dynamics. We aim to revisit the scenario of tilting Uranus in light of this mechanism. Methods. We analyse the precession spectrum of Uranus and identify the candidate secular spin-orbit resonances that could be responsible for the tilting. We determine the properties of the hypothetical ancient satellite required for a capture and explore the dynamics numerically. Results. If it migrates over 10 Uranus' radii, a single satellite with minimum mass 4e-4 Uranus' mass is able to tilt Uranus from a small obliquity and make it converge towards 90{\deg}. In order to achieve the tilting in less than the age of the Solar System, the mean drift rate of the satellite must be comparable to the Moon's current orbital expansion. Under these conditions, simulations show that Uranus is readily tilted over 80{\deg}. Beyond this point, the satellite is strongly destabilised and triggers a phase of chaotic motion for the planet's spin axis. The chaotic phase ends when the satellite collides into the planet, ultimately freezing the planet's obliquity in either a prograde, or plainly retrograde state (as Uranus today). Spin states resembling that of Uranus can be obtained with probabilities as large as 80%, but a bigger satellite is favoured, with mass 1.7e-3 Uranus' mass or more. Yet, a smaller ancient satellite is not categorically ruled out, and there is room for improving this basic scenario in future studies. Interactions among several pre-existing satellites is a promising possibility., Comment: Accepted for publication in Astronomy and Astrophysics

Optical navigation on a CubeSat must rely on the best extraction of the directions of some beacons from on-board images. We present an experiment on OPS-SAT, a CubeSat of the European Space Agency (ESA), that will characterize an onboard algorithm to this aim, named Angle-based Correlation (AbC). OPS-SAT is a 3-unit CubeSat with an Attitude Determination and Control System (ADCS) and an imager that have proved their reliability with typical performance at CubeSat scale. We selected a few star-fields that all present enough visible stars within a 10 ? field of view. When our experiment is run, OPS-SAT is pointed to the most convenient star-field at that time. There, the star-field is imaged and subwindows are extracted from the image around the expected location of each star, based on the attitude-quaternion reported by the ADCS. The AbC reconstructs the absolute direction of the central body, in principle unknown, which is the pointed known star in the experiment. The method intensively uses the quaternion algebra. The beacon location is first consolidated in the field of view with the AbC. Then, the field of view is finely positioned against the sky, again with the AbC. A covariance is associated with the found beacon direction. Our experiment with OPS-SAT manages the pointing and the imager, and processes the taken images. Then, it downloads the on-board computed absolute directions and their covariances, to be compared with the actual directions. After a campaign of intensive use of the experiment, the statistical performance of the algorithm will be established and compared to the on-board computed covariances. As a bonus, an assessment of OPS-SAT's inertial pointing stability will be available. The AbC can theoretically get rid of the Attitude Control Error (ACE) of the platform and of the Attitude Knowledge Error (AKE) estimated by the ADCS, and potentially converge to (...)

Orbital geophysical investigations of Enceladus are critical to understanding its energy balance. We identified key science questions for the geophysical exploration of Enceladus, answering which would support future assessment of Enceladus' astrobiological potential. Using a Bayesian framework, we explored how science requirements map to measurement requirements. We performed mission simulations to study the sensitivity of a single spacecraft and dual spacecraft configurations to static gravity and tidal Love numbers of Enceladus. We find that mapping Enceladus' gravity field, improving the accuracy of the physical libration amplitude, and measuring Enceladus' tidal response would provide critical constraints on the internal structure, and establish a framework for assessing Enceladus' long-term habitability. This kind of investigation could be carried out as part of a life search mission at little additional resource requirements., Comment: 21 pages, 6 figures. A paper submitted to PSJ, which is an extension of the white paper previously submitted to the decadal survey

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., Comment: Published in Nature Astronomy, SharedIt link: https://rdcu.be/b4I3Q

During the thirteen years in orbit around Saturn before its final plunge, the Cassini spacecraft provided more than ten thousand astrometric measurements. Such large amounts of accurate data enable the search for extremely faint signals in the orbital motion of the saturnian moons. Among these, the detection of the dynamical feedback of the rotation of the inner moons of Saturn on their respective orbits becomes possible. Using all the currently available astrometric data associated with Atlas, Prometheus, Pandora, Janus and Epimetheus, we first provide a detailed analysis of the Cassini Imaging Science Subsystem (ISS) data, with special emphasis on their statistical behavior and sources of bias. Then, we give updated estimates of the moons' averaged densities and try to infer more details about their interior properties by estimating the physical librations for Prometheus, Pandora, Epimetheus and Janus from anomalies in their apsidal precession. Our results are compatible with a homogeneous interior for Janus and Epimetheus, within the uncertainty of the measurements. On the other hand, we found some inconsistency for Pandora and Prometheus, which might result from a dynamical mismodeling of Saturn's gravity field. Last, we show how the synergistic introduction of libration measurements directly derived from imaging should allow the moons' moments of inertia to be better constrained., Comment: Accepted for publication in Icarus

Natural satellite systems present a large variety of orbital configurations in the solar system. While some are clearly the result of known processes, others still have largely unexplained eccentricity and inclination values. Iapetus has a still unexplained 3% orbital eccentricity and its orbital plane is tilted with respect to its local Laplace plane. On the other hand, astrometric measurements of saturnian moons have revealed high tidal migration rates, corresponding to a quality factor Q of Saturn of around 1600 for the mid-sized icy moons. We show how a past crossing of the 5:1 mean motion resonance between Titan and Iapetus may be a plausible scenario to explain Iapetus' orbit. We have carried out numerical simulations of the resonance crossing using an N-Body code as well as using averaged equations of motion. A large span of migration rates were explored for Titan and Iapetus was started on its local Laplace plane with a circular orbit. The resonance crossing can trigger a chaotic evolution of the eccentricity and the inclination of Iapetus. The outcome of the resonance is highly dependent on the migration rate. For a quality factor Q of over around 2000, the chaotic evolution of Iapetus in the resonance leads in most cases to its ejection, while simulations with a quality factor between 100 and 2000 show a departure from the resonance with post-resonant eccentricities spanning from 0 up to 15%, and free inclinations capable of reaching 11 degrees. Usually high inclinations come with high eccentricities but some simulations show elements compatible with Iapetus' current orbit. A quality factor between 100 and 2000 at the frequency of Titan would bring Titan and Iapetus into a 5:1 resonance, which would perturb Iapetus' eccentricity and inclination to values observed today. Such rapid tidal migration would have avoided Iapetus' ejection around 40 to 800 million years ago., Comment: 8 pages, 8 figures

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., Comment: 18 pages, 4 figures, accepted for publication in PSS

Accurate tracking of probes is one of the key points of space exploration. Range and Doppler techniques are the most commonly used. In this paper we analyze the impact of the transponder delay, $i.e.$ the processing time between reception and re-emission of a two-way tracking link at the satellite, on tracking observables and on spacecraft orbits. We show that this term, only partially accounted for in the standard formulation of computed space observables, can actually be relevant for future missions with high nominal tracking accuracies or for the re-processing of old missions. We present several applications of our formulation to Earth flybys, the NASA GRAIL and the ESA BepiColombo missions., Comment: 19 pages, 4 figures, 2 tables. Accepted for publication in Advances in Space Research

The closest ever fly-by of the Martian moon Phobos, performed by the European Space Agency's Mars Express spacecraft, gives a unique opportunity to sharpen and test the Planetary Radio Interferometry and Doppler Experiments (PRIDE) technique in the interest of studying planet - satellite systems. The aim of this work is to demonstrate a technique of providing high precision positional and Doppler measurements of planetary spacecraft using the Mars Express spacecraft. The technique will be used in the framework of Planetary Radio Interferometry and Doppler Experiments in various planetary missions, in particular in fly-by mode. We advanced a novel approach to spacecraft data processing using the techniques of Doppler and phase-referenced very long baseline interferometry spacecraft tracking. We achieved, on average, mHz precision (30 {\mu}m/s at a 10 seconds integration time) for radial three-way Doppler estimates and sub-nanoradian precision for lateral position measurements, which in a linear measure (at a distance of 1.4 AU) corresponds to ~50 m., Comment: 9 pages, 14 figures. Astronomy and Astrophysics, accepted on 2016/05/31

Tidal dissipation is the main driver of orbital evolution of natural satellites and a key point to understand the exoplanetary system configurations. Despite its importance, its quantification from observations still remains difficult for most objects of our own Solar system. In this work, we overview the method that has been used to determine, directly from observations, the tidal parameters, with emphasis on the Love number k2 and the tidal quality factor Q. Up-to-date values of these tidal parameters are summarized. Last, an assessment on the possible determination of the tidal ratio k2/Q of Uranus and Neptune is done. This may be particularly relevant for coming astrometric campaigns and future space missions focused on these systems., Comment: 11 pages, 3 figures, accepted in Celestial Mechanics and Dynamical Astronomy

Using astrometric observations spanning more than a century and including a large set of Cassini data, we determine Saturn's tidal parameters through their current effects on the orbits of the eight main and four coorbital moons. We have used the latter to make the first determination of Saturn's Love number, $k_2=0.390 \pm 0.024$, a value larger than the commonly used theoretical value of 0.341 (Gavrilov & Zharkov, 1977), but compatible with more recent models (Helled & Guillot, 2013) for which $k_2$ ranges from 0.355 to 0.382. Depending on the assumed spin for Saturn's interior, the new constraint can lead to a reduction of up to 80% in the number of potential models, offering great opportunities to probe the planet's interior. In addition, significant tidal dissipation within Saturn is confirmed (Lainey et al., 2012) corresponding to a high present-day tidal ratio $k_2/Q=(1.59 \pm 0.74) \times 10^{-4}$ and implying fast orbital expansions of the moons. This high dissipation, with no obvious variations for tidal frequencies corresponding to those of Enceladus and Dione, may be explained by viscous friction in a solid core, implying a core viscosity typically ranging between $10^{14}$ and $10^{16}$ Pa.s (Remus et al., 2012). However, a dissipation increase by one order of magnitude at Rhea's frequency could suggest the existence of an additional, frequency-dependent, dissipation process, possibly from turbulent friction acting on tidal waves in the fluid envelope of Saturn (Ogilvie & Li, 2004). Alternatively, a few of Saturn's moons might themselves experience large tidal dissipation.

During the last decade, the precision in the tracking of spacecraft has constantly improved. The discovery of few astrometric anomalies, such as the Pioneer and Earth flyby anomalies, stimulated further analysis of the operative modeling currently adopted in Deep Space Navigation (DSN). Our study shows that some traditional approximations lead to neglect tiny terms that could have consequences in the orbit determination of a probe in specific configurations such as during an Earth flyby. Therefore, we suggest here a way to improve the light time calculation used for probe tracking., Comment: 6 pages, 3 figures, submitted to proceedings of "International Symposium on Orbit Propagation and Determination, Lille 2011"

The origin of Saturn's inner mid-sized moons (Mimas, Enceladus, Tethys, Dione and Rhea) and Saturn's rings is debated. Charnoz et al. (2010) introduced the idea that the smallest inner moons could form from the spreading of the rings' edge while Salmon et al. (2010) showed that the rings could have been initially massive, and so was the ring's progenitor itself. One may wonder if the mid-sized moons may have formed also from the debris of a massive ring progenitor, as also suggested in Canup (2010). However, the process driving mid-sized moons accretion from the icy debris disks has not been investigated in details. In particular, this process does not seem able to explain the varying silicate contents of the mid-sized moons (from 6% to 57% in mass). Here, we explore the formation of large objects from a massive ice-rich ring (a few times Rhea's mass) and describe the fundamental properties and implications of this new process. Using a hybrid computer model, we show that accretion within massive icy rings can form all mid-sized moons from Mimas to Rhea. However in order to explain their current locations, intense dissipation within Saturn (with Qp<2000) would be required. Our results are consistent with a satellite origin tied to the rings formation at least 2.5 Gy ago, both compatible with either a formation concurrent to Saturn or during the Late Heavy Bombardment. Tidal heating related to high-eccentricity post-accretional episodes may induce early geological activity. If some massive irregular chunks of silicates were initially present within the rings, they would be present today inside the satellites' cores which would have accreted icy shells while being tidally expelled from the rings (via a heterogeneous accretion process)while those still present in the rings are interpreted as today Saturn's rings' propellers and ring-moons (like Pan or Daphnis)., Comment: 41 pages, 14 figures. Accepted in ICARUS

In this communication, we focus on the possibility to test GR with radioscience experiments. We present a new software that in a first step simulates the Range/Doppler signals directly from the space time metric (thus in GR and in alternative theories of gravity). In a second step, a least-squares fit of the involved parameters is performed in GR. This software allows one to get the order of magnitude and the signature of the modifications induced by an alternative theory of gravity on radioscience signals. As examples, we present some simulations for the Cassini mission in Post-Einsteinian gravity and with the MOND External Field Effect., Comment: 4 pages; Proceedings of "Les Rencontres de Moriond 2011 - Gravitation session"

Any model of tides is based on a specific hypothesis of how lagging depends on the tidal-flexure frequency. For example, Gerstenkorn (1955), MacDonald (1964), and Kaula (1964) assumed constancy of the geometric lag angle, while Singer (1968) and Mignard (1979, 1980) asserted constancy of the time lag. Thus, each of these two models was based on a certain law of scaling of the geometric lag. The actual dependence of the geometric lag on the frequency is more complicated and is determined by the rheology of the planet. Besides, each particular functional form of this dependence will unambiguously fix the appropriate form of the frequency dependence of the tidal quality factor, Q. Since at present we know the shape of the dependence of Q upon the frequency, we can reverse our line of reasoning and single out the appropriate actual frequency-dependence of the angular lag. This dependence turns out to be different from those employed hitherto, and it entails considerable alterations in the time scales of the tide-generated dynamical evolution. Phobos' fall on Mars is an example we consider., Comment: arXiv admin note: substantial text overlap with arXiv:astro-ph/0605521

Construction of a theory of orbits about a precessing oblate planet, in terms of osculating elements defined in a frame of the equator of date, was started in Efroimsky and Goldreich (2004) and Efroimsky (2005, 2006). We now combine that analytical machinery with numerics. The resulting semianalytical theory is then applied to Deimos over long time scales. In parallel, we carry out a purely numerical integration in an inertial Cartesian frame. The results agree to within a small margin, for over 10 Myr, demonstrating the applicability of our semianalytical model over long timescales. This will enable us to employ it at the further steps of the project, enriching the model with the tides, the pull of the Sun, and the planet's triaxiality. Another goal of our work was to check if the equinoctial precession predicted for a rigid Mars could have been sufficient to repel the orbits away from the equator. We show that for low initial inclinations, the orbit inclination reckoned from the precessing equator of date is subject only to small variations. This is an extension, to non-uniform precession given by the Colombo model, of an old result obtained by Goldreich (1965) for the case of uniform precession and a low initial inclination. However, near-polar initial inclinations may exhibit considerable variations for up to +/- 10 deg in magnitude. Nevertheless, the analysis confirms that an oblate planet can, indeed, afford large variations of the equinoctial precession over hundreds of millions of years, without repelling its near-equatorial satellites away from the equator of date: the satellite inclination oscillates but does not show a secular increase. Nor does it show secular decrease, a fact that is relevant to the discussion of the possibility of high-inclination capture of Phobos and Deimos.

In this paper we investigate the possibility of measuring the post-Newtonian general relativistic gravitomagnetic Lense-Thirring effect in the Jovian system of its Galilean satellites Io, Europa, Ganymede and Callisto in view of recent developments in processing and modelling their optical observations spanning a large time interval (125 years). The present day best observations have an accuracy between several kilometers to few tens of kilometers, which is just the order of magnitude of the Lense-Thirring shifts of the orbits of the Galilean satellites over almost a century. From a comparison between analytical development and numerical integration it turns out that, unfortunately, most of the secular component of the gravitomagnetic signature is removed in the process of fitting the initial conditions. Indeed, an estimation of the magnitude of the Lense-Thirring effect in the ephemerides residuals is given; the resulting residuals have a maximum magnitude of 20 meters only (over 125 years)., Comment: Latex, 10 pages, 4 tables, 3 figures, 27 references. Invited paper for a Special Issue of Int. J. Mod. Phys. D on the Lense-Thirring effect, D. Grumiller editor