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.

Almost all the planets of our solar system have moons. Each planetary system has however unique characteristics. The Martian system has not one single big moon like the Earth, not tens of moons of various sizes like for the giant planets, but two small moons: Phobos and Deimos. How did form such a system? This question is still being investigated on the basis of the Earth-based and space-borne observations of the Martian moons and of the more modern theories proposed to account for the formation of other moon systems. The most recent scenario of formation of the Martian moons relies on a giant impact occurring at early Mars history and having also formed the so-called hemispheric crustal dichotomy. This scenario accounts for the current orbits of both moons unlike the scenario of capture of small size asteroids. It also predicts a composition of disk material as a mixture of Mars and impactor materials that is in agreement with remote sensing observations of both moon surfaces, which suggests a composition different from Mars. The composition of the Martian moons is however unclear, given the ambiguity on the interpretation of the remote sensing observations. The study of the formation of the Martian moon system has improved our understanding of moon formation of terrestrial planets: The giant collision scenario can have various outcomes and not only a big moon as for the Earth. This scenario finds a natural place in our current vision of the early solar system when conditions were favorable for giant collisions to occur. The next step in exploration of Martian moon is a sample return mission to test the giant collision scenario for their origin, and to provide tests of models of early solar system dynamics since Mars may retain material exchanged between the inner and outer solar system.

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.

Enceladus's gravity and shape have been explained in terms of a thick isostatic ice shell floating on a global ocean, in contradiction of the thin shell implied by librations. Here we propose a new isostatic model minimizing crustal deviatoric stress, and demonstrate that gravity and shape data predict a $\rm{38\pm4\,km}$-thick ocean beneath a $\rm{23\pm4\,km}$-thick shell agreeing with -- but independent of -- libration data. Isostatic and tidal stresses are comparable in magnitude. South polar crust is only $7\pm4\rm\,km$ thick, facilitating the opening of water conduits and enhancing tidal dissipation through stress concentration. Enceladus's resonant companion, Dione, is in a similar state of minimum stress isostasy. Its gravity and shape can be explained in terms of a $\rm{99\pm23\,km}$-thick isostatic shell overlying a $\rm{65\pm30\,km}$-thick global ocean, thus providing the first clear evidence for a present-day ocean within Dione.