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.

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.

Several natural satellites of the giant planets have shown evidence of a global internal ocean, coated by a thin, icy crust. This crust is probably viscoelastic, which would alter its rotational response. This response would translate into several rotational quantities, i.e. the obliquity, and the librations at different frequencies, for which the crustal elasticity reacts differently. This study aims at modelling the global response of the viscoelastic crust. For that, I derive the time-dependency of the tensor of inertia, which I combine with the time evolution of the rotational quantities, thanks to an iterative algorithm. This algorithm combines numerical simulations of the rotation with a digital filtering of the resulting tensor of inertia. The algorithm works very well in the elastic case, provided the problem is not resonant. However, considering tidal dissipation adds different phase lags to the oscillating contributions, which challenge the convergence of the algorithm.

A Hamiltonian model is constructed for the spin axis of a planet perturbed by a nearby planet with both planets in orbit about star. We expand the planet-planet gravitational potential perturbation to first order in orbital inclinations and eccentricities, finding terms describing spin resonances involving the spin precession rate and the two planetary mean motions. Convergent planetary migration allows the spinning planet to be captured into spin resonance. With initial obliquity near zero, the spin resonance can lift the planet's obliquity to near 90 or 180 degrees depending upon whether the spin resonance is first or zero-th order in inclination. Past capture of Uranus into such a spin resonance could give an alternative non-collisional scenario accounting for Uranus's high obliquity. However we find that the time spent in spin resonance must be so long that this scenario cannot be responsible for Uranus's high obliquity. Our model can be used to study spin resonance in satellite systems. Our Hamiltonian model explains how Styx and Nix can be tilted to high obliquity via outward migration of Charon, a phenomenon previously seen in numerical simulations.

New Horizons mission observations show that the small satellites Styx, Nix, Kerberos and Hydra, of the Pluto-Charon system, have not tidally spun-down to near synchronous spin states and have high obliquities with respect to their orbit about the Pluto-Charon binary (Weaver et al. 2016). We use a damped mass-spring model within an N-body simulation to study spin and obliquity evolution for single spinning non-round bodies in circumbinary orbit. Simulations with tidal dissipation alone do not show strong obliquity variations from tidally induced spin-orbit resonance crossing and this we attribute to the high satellite spin rates and low orbital eccentricities. However, a tidally evolving Styx exhibits intermittent obliquity variations and episodes of tumbling. During a previous epoch where Charon migrated away from Pluto, the minor satellites could have been trapped in orbital mean motion inclination resonances. An outward migrating Charon induces large variations in Nix and Styx's obliquities. The cause is a commensurability between the mean motion resonance frequency and the spin precession rate of the spinning body. As the minor satellites are near mean motion resonances, this mechanism could have lifted the obliquities of all four minor satellites. The high obliquities need not be primordial if the minor satellites were at one time captured into mean motion resonances.

Most of the main planetary satellites of our Solar System are expected to be in synchronous rotation, the departures from the strict synchronicity being a signature of the interior. Librations have been measured for the Moon, Phobos, and some satellites of Saturn. I here revisit the theory of the longitudinal librations in considering that part of the interior is not hydrostatic, i.e. has not been shaped by the rotational and tidal deformations, but is fossil. This consideration affects the rotational behavior. For that, I derive the tensor of inertia of the satellite in splitting these two parts, before proposing an analytical solution that I validate with numerical simulations. I apply this new theory on Mimas and Epimetheus, for which librations have been measured from Cassini data. I show that the large measured libration amplitude of these bodies can be explained by an excess of triaxiality that would not result from the hydrostatic theory. This theory cannot explain the phase shift which has been measured in the diurnal librations of Mimas. This speaks against a solid structure for Mimas, i.e. Mimas could have a global internal ocean.

This HDR-thesis is devoted to the study of the rotation of the natural satellites of the giant planets and of Mercury. These bodies have a resonant rotation. Most of the natural satellites rotate synchronously, showing the same hemisphere to their parent planet (1:1 spin-orbit resonance). The case of Mercury is unique since its spin rate is exactly 1.5 its mean motion (3:2 spin-orbit resonance). These two configurations are dynamical equilibria, reached after damping of the initial rotation of the relevant bodies. Thus, the rotation quantities are a signature of the interior, in particular of a putative global ocean. This manuscript divides into 3 parts. The first part is devoted to the synchronous resonance. It presents different models of rotation from a fully rigid body to a one with a global subsurfacic ocean. We always consider all the degrees of freedom simultaneously, using analytical and numerical resolutions. These models are applied on Titan, Callisto, Janus, Epimetheus, Mimas, Hyperion, and Io. The second part presents the resonant rotation of Mercury, target of the two space missions MESSENGER and BepiColombo. We reveal in particular how it got trapped into its 3:2 resonance. The final part presents an algorithm I have elaborated to tackle the rotational problems.

While it is accepted that the eccentricity of Mercury (0.206) favours entrapment into the 3:2 spin-orbit resonance, open is the question how and when the capture took place. A recent work by Makarov (2012) has demonstrated that trapping into this resonance is certain if the eccentricity is larger than 0.2, provided that we use a realistic tidal model, the one which is based on the Darwin-Kaula expansion of the tidal torque. The physics-based tidal model changes dramatically the statistics of the possible final spin states. First, we discover that after only one encounter with the spin-orbit 3:2 resonance this resonance becomes the most probable end-state. Second, if a capture into this (or any other) resonance takes place, the capture becomes final, several crossings of the same state being forbidden by our model. Third, within our model the trapping of Mercury happens much faster than previously believed: for most histories, 10 - 20 Myr are sufficient. Fourth, even a weak laminar friction between the solid mantle and a molten core would most likely result in a capture in the 2:1 or even higher resonance. So the principal novelty of our paper is that the 3:2 end-state is more ancient than the same end-state obtained when the constant time lag model is employed. The swift capture justifies our treatment of Mercury as a homogeneous, unstratified body whose liquid core had not yet formed by the time of trapping. We also provide a critical analysis of the hypothesis by Wieczorek et al. (2012) that the early Mercury might had been retrograde, whereafter it synchronised its spin and then accelerated it to the 3:2 resonance. Accurate processing of the available data on cratering does not support that hypothesis, while the employment of a realistic rheology invalidates a key element of the hypothesis, an intermediate pseudosynchronous state needed to spin-up to the 3:2 resonance.

We constructed a 6-degrees of freedom rotational model of Titan as a 3-layer body consisting of a rigid core, a fluid global ocean, and a floating ice shell. The ice shell exhibits partially-compensated lateral thickness variations in order to simultaneously match the observed degree-two gravity and shape coefficients. The rotational dynamics are affected by the gravitational torque of Saturn, the gravitational coupling between the inner core and the shell, and the pressure coupling at the fluid-solid boundaries. Between 10 and 13% of our model Titans have an obliquity (due to a resonance with the 29.5-year periodic annual forcing) that is consistent with the observed value. The shells of the successful models have a mean thickness of 130 to 140 km, and an ocean of ~250 km thickness. Our simulations of the obliquity evolution show that the Cassini obliquity measurement is an instantaneous one, and does not represent a mean value. Future measurements of the time derivative of the obliquity would help to refine the interior models. We expect in particular a variation of roughly 7 arcmin over the duration of the Cassini mission.

Miranda has a unusually high inclination ($I=4.338^\circ$), and its surface reveals signs of past endogenic activity. Investigations of the dynamical aspects of its orbital evolution suggest probable resonant processes, in particular with Umbriel, as an explanation for the present high inclination of Miranda. The tidal heating induced by gravitational interactions can lead to the rise of eccentricities and, consequently, to the increased dissipation of energy inside the satellite and higher internal temperatures. We study here the possible increase in eccentricities caused by orbital resonances and the resulting endogenic heating on Miranda taking into account its temperature dependent rheology. The coupled orbital-thermal evolution model was run with different rheological models and the thermal parameters starting form a cold thermal state, in radiative equilibrium with the environment. For the nominal parameters of the evolution scenarios studied, the resonances were not sufficient to rise neither the eccentricities nor the internal temperatures significantly. Lowest dissipation function $Q$ of around 100 and final eccentricity of $e\approx0.02$ were obtained during the resonance 3:1 with Umbriel.

Earth-based radar observations of the rotational dynamics of Mercury (Margot et al. 2012) combined with the determination of its gravity field by MESSENGER (Smith et al. 2012) give clues on the internal structure of Mercury, in particular its polar moment of inertia C, deduced from the obliquity (2.04 +/- 0.08) arcmin. The dynamics of the obliquity of Mercury is a very-long term motion (a few hundreds of kyrs), based on the regressional motion of Mercury's orbital ascending node. This paper, following the study of Noyelles & D'Hoedt (2012), aims at first giving initial conditions at any time and for any values of the internal structure parameters for numerical simulations, and at using them to estimate the influence of usually neglected parameters on the obliquity, like J3, the Love number k2 and the secular variations of the orbital elements. We use for that averaged representations of the orbital and rotational motions of Mercury, suitable for long-term studies. We find that J3 should alter the obliquity by 250 milli-arcsec, the tides by 100 milli-arcsec, and the secular variations of the orbital elements by 10 milli-arcsec over 20 years. The resulting value of C could be at the most changed from 0.346mR^2 to 0.345mR^2.

The Uranian satellite Miranda presents a high inclination (4.338{\deg}) and evidences of resurfacing. It is accepted since 20 years (e.g. Tittemore and Wisdom 1989, Malhotra and Dermott 1990) that this inclination is due to the past trapping into the 3:1 resonance with Umbriel. These last years there is a renewal of interest for the Uranian system since the Hubble Space Telescope permitted the detection of an inner system of rings and small embedded satellites, their dynamics being of course ruled by the main satellites. For this reason, we here propose to revisit the long-term dynamics of Miranda, using modern tools like intensive computing facilities and new chaos indicators (MEGNO and frequency map analysis). As in the previous studies, we find the resonance responsible for the inclination of Miranda and the secondary resonances associated, likely to have stopped the rise of Miranda's inclination at 4.5{\deg}. Moreover, we get other trajectories in which this inclination reaches 7{\deg}. We also propose an analytical study of the secondary resonances associated, based on the study by Moons and Henrard (1993).

This note tackles the problem of the rotation of Io with the 4-degrees of freedom Poincar\'e-Hough model. Io is modeled as a 2-layer body, i.e. a triaxial fluid core and a rigid outer layer. We show that the longitudinal librations should have an amplitude of about 30 arcseconds, independent of the composition of the core. We also estimate the tidal instability of the core, and show that should be slowly unstable.

This paper presents a study of the Poincar\'e-Hough model of rotation of the synchronous natural satellites, in which these bodies are assumed to be composed of a rigid mantle and a triaxial cavity filled with inviscid fluid of constant uniform density and vorticity. In considering an Io-like body on a low eccentricity orbit, we describe the different possible behaviors of the system, depending on the size, polar flattening and shape of the core. We use for that the numerical tool. We propagate numerically the Hamilton equations of the systems, before expressing the resulting variables under a quasi-periodic representation. This expression is obtained numerically by frequency analysis. This allows us to characterize the equilibria of the system, and to distinguish the causes of their time variations. We show that, even without orbital eccentricity, the system can have complex behaviors, in particular when the core is highly flattened. In such a case, the polar motion is forced by several degrees and longitudinal librations appear. This is due to splitting of the equilibrium position of the polar motion. We also get a shift of the obliquity when the polar flattening of the core is small.

The Cassini mission in the Saturnian system is an outstanding opportunity to improve our knowledge of the satellites of Saturn. The data obtained thanks to this mission must be confronted to theoretical models. This paper aims at modeling the rotation of Mimas, with respect to its possible internal structure. For that, we first build different interior models, in considering Mimas as composed of 2 rigid layers with different porosity. Then we simulate the rotational behavior of these models in a 3-degree of freedom numerical code, in considering complete ephemerides of a Mimas whose rotation is disturbed by Saturn. We also estimate the deviation of its longitudinal orientation due to tides. We expect a signature of the internal structure up to 0.53{\deg} in the longitudinal librations and an obliquity between 2 and 3 arcmin, the exact values depending on the interior. The longitudinal librations should be detectable, inverting them to get clues on the internal structure of Mimas is challenging.

Mercury is the target of two space missions: MESSENGER (NASA) which orbit insertion is planned for March 2011, and ESA/JAXA BepiColombo, that should be launched in 2014. Their instruments will observe the surface of the planet with a high accuracy (about 1 arcsec for BepiColombo), what motivates studying its rotation. Mercury is assumed to be composed of a rigid mantle and an at least partially molten core. We here study the influence of the core-mantle interactions on the rotation perturbed by the solar gravitational interaction, by modeling the core as an ellipsoidal cavity filled with inviscid fluid of constant uniform density and vorticity. We use both analytical (Lie transforms) and numerical tools to study this rotation, with different shapes of the core. We express in particular the proper frequencies of the system, because they characterize the response of Mercury to the different solicitations, due to the orbital motion of Mercury around the Sun. We show that, contrary to its size, the shape of the core cannot be determined from observations of either longitudinal or polar motions. However, we highlight the strong influence of a resonance between the proper frequency of the core and the spin of Mercury that raises the velocity field inside the core. We show that the key parameter is the polar flattening of the core. This effect cannot be directly derived from observations of the surface of Mercury, but we cannot exclude the possibility of an indirect detection by measuring the magnetic field.

The Saturnian coorbital satellites Janus and Epimetheus present a unique dynamical configuration in the Solar System, because of high-amplitude horseshoe orbits, due to a mass ratio of order unity. As a consequence, they swap their orbits every 4 years, while their orbital periods is about 0.695 days. Recently, Tiscareno et al.(2009) got observational informations on the shapes and the rotational states of these satellites. In particular, they detected an offset in the expected equilibrium position of Janus, and a large libration of Epimetheus. We here propose to give a 3-dimensional theory of the rotation of these satellites in using these observed data, and to compare it to the observed rotations. We consider the two satellites as triaxial rigid bodies, and we perform numerical integrations of the system in assuming the free librations as damped. The periods of the three free librations we get, associated with the 3 dimensions, are respectively 1.267, 2.179 and 2.098 days for Janus, and 0.747, 1.804 and 5.542 days for Epimetheus. The proximity of 0.747 days to the orbital period causes a high sensitivity of the librations of Epimetheus to the moments of inertia. Our theory explains the amplitude of the librations of Janus and the error bars of the librations of Epimetheus, but not an observed offset in the orientation of Janus.

In Noyelles et al. (2008, Astron. Astrophys., 478, 959-970), a resonance involving the wobble of Titan is hinted. This paper studies this scenario and its consequences. The first step is to build an accurate analytical model that would help to find the likely resonances in the rotation of every synchronous body. In this model, I take the orbital eccentricity of the body into account, as well as its variable inclination with respect to Saturn's equator. Then an analytical study using the second fundamental model of the resonance is performed to study the resonance of interest. Finally, I study the dissipative consequences of this resonance. I find that this resonance may have increased the wobble of Titan by several degrees. For instance, if Titan's polar momentum C is equal to 0.355MR_T^2 (M and R_T being respectively Titan's mass and radius), the wobble might be forced to 41 degrees. Thanks to an original formula, I find that the dissipation associated with the forced wobble might not be negligible compared to the contribution of the eccentricity. I also suspect that, due to the forced wobble, Titan's period of rotation may be somewhat underestimated by observers. Finally, I use the analytical model presented in this paper to compute the periods of the free librations of the four Galilean satellites as well as the Saturnian satellite Rhea. For Io and Europa, the results are consistent with previous studies. For the other satellites, the periods of the free librations are respectively 186.37 d, 23.38 y and 30.08 y for Ganymede, 2.44 y, 209.32 y and 356.54 y for Callisto, and 51.84 d, 2.60 y and 3.59 y for Rhea.

We study the forced rotation of Titan seen as a rigid body at the equilibrium Cassini state, involving the spin-orbit synchronization. We used both the analytical and the numerical ways. We analytically determined the equilibrium positions and the frequencies of the 3 free librations around it, while a numerical integration associated to frequency analysis gave us a more synthetic, complete theory, where the free solution split from the forced one. We find a mean obliquity of 2.2 arcmin and the fundamental frequencies of the free librations of about 2.0977, 167.4883, and 306.3360 years. Moreover, we bring out the main role played by Titan's inclination on its rotation, and we suspect a likely resonance involving Titan's wobble.