We compute the inertial modes of a freely rotating two-layer planetary model with an ellipsoidal inviscid fluid core and a perfectly rigid mantle. We present a method to derive analytical formulae for the frequencies of the Free Core Nutation (FCN) and Chandler Wobble (CW) which are valid to all orders of the dynamical flattening of the core and mantle, and we show how the FCN and CW are the direct generalisation of the purely fluid Spin-Over mode (SO) and of the Eulerian Wobble (EW) to the case where the mantle can oscillate freely around a state of steady rotation. Through a numerical computation for an axisymmetric (oblate spheroidal) planet, we demonstrate that all other inertial modes of the steadily rotating fluid core are also free modes of the freely rotating two-layer planet.

Two years after InSight's arrival at Mars, the ExoMars 2020 mission will land on the opposite side of the red planet. Similarly to InSight, which carries the RISE (Rotation and Interior Structure Experiment) radio-science experiment, the ExoMars mission will have on board the Lander Radio-science (LaRa) experiment. The X-band transponders on RISE and LaRa, allowing for direct radio-link between the landers and stations on Earth, are dedicated to the investigation of Mars' deep interior through the precise measurement of the planet's rotation and orientation. The benefit of having LaRa after RISE for the determination of the Mars orientation and rotation parameters is demonstrated and the resulting improved constraints on the interior structure of Mars and, in particular, on its core are quantified via numerical simulations. In particular, we show that the amplitudes of the semi-annual prograde ($p_2)$ and the ter-annual retrograde ($r_3$) nutations will be determined with a precision of 6 and 4 milliarcseconds respectively by combining 700 days of RISE data with 700 days of LaRa data, about 35$\%$ more precise than what is expected from RISE alone. The impact of such an improvement on the determination of the core size of Mars is discussed and shown to be significant.

LaRa (Lander Radioscience) is an experiment on the ExoMars 2020 mission that uses the Doppler shift on the radio link due to the motion of the ExoMars platform tied to the surface of Mars with respect to the Earth ground stations (e.g. the deep space network stations of NASA), in order to precisely measure the relative velocity of the lander on Mars with respect to the Earth. The LaRa measurements shall improve the understanding of the structure and processes in the deep interior of Mars by obtaining the rotation and orientation of Mars with a better precision compared to the previous missions. In this paper, we provide the analysis done until now for the best realization of these objectives. We explain the geophysical observation that will be reached with LaRa (Length-of-day variations, precession, nutation, and possibly polar motion). We develop the experiment set up, which includes the ground stations on Earth (so-called ground segment). We describe the instrument, i.e. the transponder and its three antennas. We further detail the link budget and the expected noise level that will be reached. Finally, we detail the expected results, which encompasses the explanation of how we shall determine Mars' orientation parameters, and the way we shall deduce Mars' interior structure and Mars' atmosphere from them. Lastly, we explain briefly how we will be able to determine the Surface platform position.

This paper reviews habitability conditions for a terrestrial planet from the point of view of geosciences. It addresses how interactions between the interior of a planet or a moon and its atmosphere and surface (including hydrosphere and biosphere) can affect habitability of the celestial body. It does not consider in detail the role of the central star but focusses more on surface conditions capable of sustaining life. We deal with fundamental issues of planetary habitability, i.e. the environmental conditions capable of sustaining life, and the above-mentioned interactions can affect the habitability of the celestial body. We address some hotly debated questions including: - How do core and mantle affect the evolution and habitability of planets? - What are the consequences of mantle overturn on the evolution of the interior and atmosphere? - What is the role of the global carbon and water cycles? - What influence do comet and asteroid impacts exert on the evolution of the planet? - How does life interact with the evolution of the Earth's geosphere and atmosphere? - How can knowledge of the solar system geophysics and habitability be applied to exoplanets? In addition, we address the identification of preserved life tracers in the context of the interaction of life with planetary evolution.

The Earth is a rapidly rotating body. The centrifugal pull makes its shape resemble a flattened ellipsoid and Coriolis forces support waves in its fluid core, known as inertial waves. These waves can lead to global oscillations, or modes, of the fluid. Periodic variations of the Earth's rotation axis (nutations) can lead to an exchange of angular momentum between the mantle and the fluid core and excite these inertial modes. In addition to viscous torques that exist regardless of the shape of the boundaries, the small flattening of the core-mantle boundary (CMB) allows inertial modes to exert pressure torques on the mantle. These torques effectively couple the rigid-body dynamics of the Earth with the fluid dynamics of the fluid core. Here we present the first high resolution numerical model that solves simultaneously the rigid body dynamics of the mantle and the Navier-Stokes equation for the liquid core. This method takes naturally into account dissipative processes in the fluid that are ignored in current nutation models. We find that the Free Core Nutation (FCN) mode, mostly a toroidal fluid flow if the mantle has a large moment of inertia, enters into resonance with nearby modes if the mantle's moment of inertia is reduced. These mode interactions seem to be completely analogous to the ones discovered by Schmitt (2006) in a uniformly rotating ellipsoid with varying flattening.

We propose a numerical method to compute the inertial modes of a container with near-spherical geometry based on the fully spectral discretisation of the angular and radial directions using spherical harmonics and Gegenbauer polynomial expansion respectively. This allows to solve simultaneously the Poincare equation and the no penetration condition as an algebraic polynomial eigenvalue problem. The inertial modes of an exact oblate spheroid are recovered to machine precision using an appropriate set of spheroidal coordinates. We show how other boundaries that deviate slightly from a sphere can be accommodated for with the technique of equivalent spherical boundary and we demonstrate the convergence properties of this approach for the triaxial ellipsoid.

Enceladus is characterised by a south polar hot spot associated with a large outflow of heat, the source of which remains unclear. We compute the viscous dissipation resulting from tidal and libration forcing in the moon's subsurface ocean using the linearised Navier-Stokes equation in a 3-dimensional spherical model. We conclude that libration is the dominant cause of dissipation at the linear order, providing up to about 0.001 GW of heat to the ocean, which remains insufficient to explain the (about) 10 GW observed by Cassini. We also illustrate how resonances with inertial modes can significantly augment the dissipation. Our work is an extension to Rovira-Navarro et al. [2019] to include the effects of libration. The model developed here is readily applicable to the study of other moons and planets.

The position of a Martian lander is affected by different aspects of Mars' rotational motions: the nutations, the precession, the length-of-day variations and the polar motion. These various motions have a different signature in a Doppler observable between the Earth and a lander on Mars' surface. Knowing the correlations between these signatures and the moments when these signatures are not null during one day or on a longer timescale is important to identify strategies that maximize the geophysical return of observations with a geodesy experiment, in particular for the ones on-board the future NASA InSight or ESA-Roscosmos ExoMars2020 missions. We provide first-order formulations of the signature of the rotation parameters in the Doppler and range observables. These expressions are functions of the diurnal rotation of Mars, the lander position, the planet radius and the rotation parameter. Additionally, the nutation signature in the Doppler observable is proportional to the Earth declination with respect to Mars. For a lander on Mars close to the equator, the motions with the largest signature in the Doppler observable are due to the length-of-day variations, the precession rate and the rigid nutations. The polar motion and the liquid core signatures have a much smaller amplitude. For a lander closer to the pole, the polar motion signature is enhanced while the other signatures decrease. We also numerically evaluate the amplitudes of the rotation parameters signature in the Doppler observable for landers on other planets or moons.

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.

Accurate positional measurements of planets and satellites are used to improve our knowledge of their dynamics and to infer the accuracy of planet and satellite ephemerides. In the framework of the FP7 ESPaCE project, we provide the positions of Mars, Phobos, and Deimos taken with the U.S. Naval Observatory 26-inch refractor during the 1971 opposition of the planet. These plates were measured with the digitizer of the Royal Observatory of Belgium and reduced through an optimal process that includes image, instrumental, and spherical corrections to provide the most accurate data. We compared the observed positions of the planet Mars and its satellites with the theoretical positions from INPOP10 and DE430 planetary ephemerides, and from NOE and MAR097 satellite ephemerides. The rms residuals in RA and Dec. of one position is less than 60 mas, or about 20 km at Mars. This accuracy is comparable to the most recent CCD observations. Moreover, it shows that astrometric data derived from photographic plates can compete with those of old spacecraft (Mariner 9, Viking 1 and 2).

The evolution and escape of the martian atmosphere and the planet's water inventory can be separated into an early and late evolutionary epoch. The first epoch started from the planet's origin and lasted $\sim$500 Myr. Because of the high EUV flux of the young Sun and Mars' low gravity it was accompanied by hydrodynamic blow-off of hydrogen and strong thermal escape rates of dragged heavier species such as O and C atoms. After the main part of the protoatmosphere was lost, impact-related volatiles and mantle outgassing may have resulted in accumulation of a secondary CO$_2$ atmosphere of a few tens to a few hundred mbar around $\sim$4--4.3 Gyr ago. The evolution of the atmospheric surface pressure and water inventory of such a secondary atmosphere during the second epoch which lasted from the end of the Noachian until today was most likely determined by a complex interplay of various nonthermal atmospheric escape processes, impacts, carbonate precipitation, and serpentinization during the Hesperian and Amazonian epochs which led to the present day surface pressure.

In this paper, we focus on the possibility to test General Relativity in the Solar System with radioscience measurements. To this aim, we present a new software that simulates Range and Doppler signals directly from the space-time metric. This flexible approach allows one to perform simulations in General Relativity and in alternative metric theories of gravity. In a second step, a least-squares fit of the different initial conditions involved in the situation is performed in order to compare anomalous signals produced by a given alternative theory with the ones obtained in General Relativity. This software provides orders of magnitude and signatures stemming from hypothetical alternative theories of gravity on radioscience signals. As an application, we present some simulations done for the Cassini mission in Post-Einsteinian Gravity and in the context of MOND External Field Effect. We deduce constraints on the Post-Einsteinian parameters but find that the considered arc of the Cassini mission is not useful to constrain the MOND External Field Effect.

The laws of gravitation have been tested for a long time with steadily improving precision, leading at some moment of time to paradigmatic evolutions. Pursuing this continual effort is of great importance for science. In this communication, we focus on Solar System tests of gravity and more precisely on possible tests that can be performed with radio science observations (Range and Doppler). After briefly reviewing the current tests of gravitation at Solar System scales, we give motivations to continue such experiments. In order to obtain signature and estimate the amplitude of anomalous signals that could show up in radio science observables because of modified gravitational laws, we developed a new software that simulates Range/Doppler signals. We present this new tool that simulates radio science observables directly from the space-time metric. We apply this tool to the Cassini mission during its cruise from Jupiter to Saturn and derive constraints on the parameters entering alternative theories of gravity beyond the standard Parametrized Post Newtonian theory.

Context. Ceres is the most massive body of the asteroid belt and contains about 25 wt.% (weight percent) of water. Understanding its thermal evolution and assessing its current state are major goals of the Dawn Mission. Constraints on internal structure can be inferred from various observations. Especially, detailed knowledge of the rotational motion can help constrain the mass distribution inside the body, which in turn can lead to information on its geophysical history. Aims. We investigate the signature of the interior on the rotational motion of Ceres and discuss possible future measurements performed by the spacecraft Dawn that will help to constrain Ceres' internal structure. Methods. We compute the polar motion, precession-nutation, and length-of-day variations. We estimate the amplitudes of the rigid and non-rigid response for these various motions for models of Ceres interior constrained by recent shape data and surface properties. Results. As a general result, the amplitudes of oscillations in the rotation appear to be small, and their determination from spaceborne techniques will be challenging. For example, the amplitudes of the semi-annual and annual nutations are around ~364 and ~140 milli-arcseconds, and they show little variation within the parametric space of interior models envisioned for Ceres. This, combined with the very long-period of the precession motion, requires very precise measurements. We also estimate the timescale for Ceres' orientation to relax to a generalized Cassini State, and we find that the tidal dissipation within that object was probably too small to drive any significant damping of its obliquity since formation. However, combining the shape and gravity observations by Dawn offers the prospect to identify departures of non-hydrostaticity at the global and regional scale, which will be instrumental in constraining Ceres' past and current thermal state. We also discuss the existence of a possible Chandler mode in the rotational motion of Ceres, whose potential excitation by endogenic and/or exogenic processes may help detect the presence of liquid reservoirs within the asteroid.

A major goal of the BepiColombo mission to Mercury is the determination of the structure and state of Mercury's interior. Here the BepiColombo rotation experiment has been simulated in order to assess the ability to attain the mission goals and to help lay out a series of constraints on the experiment's possible progress. In the rotation experiment pairs of images of identical surface regions taken at different epochs are used to retrieve information on Mercury's rotation and orientation. The idea is that from observations of the same patch of Mercury's surface at two different solar longitudes of Mercury the orientation of Mercury can be determined, and therefore also the obliquity and rotation variations with respect to the uniform rotation. The estimation of the libration amplitude and obliquity through pattern matching of observed surface landmarks is challenging. The main problem arises from the difficulty to observe the same landmark on the planetary surface repeatedly over the MPO mission lifetime, due to the combination of Mercury's 3:2 spin-orbit resonance, the absence of a drift of the MPO polar orbital plane and the need to combine data from different instruments with their own measurement restrictions. By assuming that Mercury occupies a Cassini state and that the spacecraft operates nominally we show that under worst case assumptions the annual libration amplitude and obliquity can be measured with a precision of respectively 1.4 arcseconds (as) and 1.0 as over the nominal BepiColombo MPO lifetime with about 25 landmarks for rather stringent illumination restrictions. The outcome of the experiment cannot be easily improved by simply relaxing the observational constraints, or increasing the data volume.

Two papers recently published in Celestial Mechanics (Krasinsky 2006, and Krasinsky and Vasilyev 2006) have presented a model for Earth-rotation variations, called ERA-2005, based on numerical integration of a new set of equations for the rotation of a deformable Earth followed by a fit of the results of the integration to VLBI data. These papers claimed that this model was superior to any other existing model. The purpose of this Note is to bring to light fundamental errors in the derivation of the basic equations of the new theory, compounded by serious deficiencies in the process of fitting to the data; they make ERA-2005 unsuitable for consideration as a geophysics-based model of nutation and precession.