In the framework of the planar and circular restricted three-body problem, we consider an asteroid that orbits the Sun in quasi-satellite motion with a planet. A quasi-satellite trajectory is a heliocentric orbit in co-orbital resonance with the planet, characterized by a non zero eccentricity and a resonant angle that librates around zero. Likewise, in the rotating frame with the planet it describes the same trajectory as the one of a retrograde satellite even though the planet acts as a perturbator. In the last few years, the discoveries of asteroids in this type of motion made the term "quasi-satellite" more and more present in the literature. However, some authors rather use the term "retrograde satellite" when referring to this kind of motion in the studies of the restricted problem in the rotating frame. In this paper we intend to clarify the terminology to use, in order to bridge the gap between the perturbative co-orbital point of view and the more general approach in the rotating frame. Through a numerical exploration of the co-orbital phase space, we describe the quasi-satellite domain and highlight that it is not reachable by low eccentricities by averaging process. We will show that the quasi-satellite domain is effectively included in the domain of the retrograde satellites and neatly defined in terms of frequencies. Eventually, we highlight a remarkable high eccentric quasi-satellite orbit corresponding to a frozen ellipse in the heliocentric frame. We extend this result to the eccentric case (planet on an eccentric motion) and show that two families of frozen ellipses originate from this remarkable orbit.

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.

We revisit the dynamics of Atlas. Using Cassini ISS astrometric observations spanning February 2004 to August 2013, Cooper et al. (2015) found evidence that Atlas is currently perturbed by both a 54:53 corotation eccentricity resonance (CER) and a 54:53 Lindblad eccentricity resonance (LER) with Prometheus. They demonstrated that the orbit of Atlas is chaotic, with a Lyapunov time of order 10 years, as a direct consequence of the coupled resonant interaction (CER/LER) with Prometheus. Here we investigate the interactions between the two resonances using the CoraLin analytical model (El Moutamid et al. 2014), showing that the chaotic zone fills almost all the corotation sites occupied by the satellite's orbit. Four 70:67 apse-type mean motion resonances with Pandora are also overlapping, but these resonances have a much weaker effect. Frequency analysis allows us to highlight the coupling between the 54:53 resonances, and confirms that a simplified system including the perturbations due to Prometheus and Saturn's oblateness only captures the essential features of the dynamics.

This paper describes several applications in astronomy and cosmology that are addressed using probabilistic modelling and statistical inference.

In the circular case of the coplanar Restricted Three-body Problem, we studied how the family of quasi-satellite (QS) periodic orbits allows to define an associated libration center. Using the averaged problem, we highlighted a validity limit of this one: for QS orbits with low eccentricities, the averaged problem does not correspond to the real problem. We do the same procedure to L 3 , L 4 and L 5 emerging periodic orbits families and remarked that for very high eccentricities F L4 and F L5 merge with F L3 which bifurcates to a stable family.

This paper tackles important aspects of comets dynamics from a statistical point of view. Existing methodology uses numerical integration for computing planetary perturbations for simulating such dynamics. This operation is highly computational. It is reasonable to wonder whenever statistical simulation of the perturbations can be much more easy to handle. The first step for answering such a question is to provide a statistical study of these perturbations in order to catch their main features. The statistical tools used are order statistics and heavy tail distributions. The study carried out indicated a general pattern exhibited by the perturbations around the orbits of the important planet. These characteristics were validated through statistical testing and a theoretical study based on Opik theory.

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.