When dealing with satellites orbiting a central body on a highly elliptical orbit, it is necessary to consider the effect of gravitational perturbations due to external bodies. Indeed, these perturbations can become very important as soon as the altitude of the satellite becomes high, which is the case around the apocentre of this type of orbit. For several reasons, the traditional tools of celestial mechanics are not well adapted to the particular dynamic of highly elliptical orbits. On the one hand, analytical solutions are quite generally expanded into power series of the eccentricity and therefore limited to quasi-circular orbits [17, 25]. On the other hand, the time-dependency due to the motion of the third-body is often neglected. We propose several tools to overcome these limitations. Firstly, we have expanded the disturbing function into a finite polynomial using Fourier expansions of elliptic motion functions in multiple of the satellite's eccentric anomaly (instead of the mean anomaly) and involving Hansen-like coefficients. Next, we show how to perform a normalization of the expanded Hamiltonian by means of a time-dependent Lie transformation which aims to eliminate periodic terms. The difficulty lies in the fact that the generator of the transformation must be computed by solving a partial differential equation involving variables which are linear with time and the eccentric anomaly which is not time linear. We propose to solve this equation by means of an iterative process.

Recent studies have shown that the medium-Earth orbit (MEO) region of the Global Navigation Satellite Systems is permeated by a devious network of lunisolar secular resonances, which can interact to produce chaotic and diffusive motions. The precarious state of the four navigation constellations, perched on the threshold of instability, makes it understandable why all past efforts to define stable graveyard orbits, especially in the case of Galileo, were bound to fail; the region is far too complex to allow of an adoption of the simple geosynchronous disposal strategy. We retrace one such recent attempt, funded by ESA's General Studies Programme in the frame of the GreenOPS initiative, that uses a systematic parametric approach and the straightforward maximum-eccentricity method to identify long-term stable regions, suitable for graveyards, as well as large-scale excursions in eccentricity, which can be used for post-mission deorbiting of constellation satellites. We then apply our new results on the stunningly rich dynamical structure of the MEO region toward the analysis of these disposal strategies for Galileo, and discuss the practical implications of resonances and chaos in this regime. We outline how the identification of the hyperbolic and elliptic fixed points of the resonances near Galileo can lead to explicit criteria for defining optimal disposal strategies.

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

It has long been suspected that the Global Navigation Satellite Systems exist in a background of complex resonances and chaotic motion; yet, the precise dynamical character of these phenomena remains elusive. Recent studies have shown that the occurrence and nature of the resonances driving these dynamics depend chiefly on the frequencies of nodal and apsidal precession and the rate of regression of the Moon's nodes. Woven throughout the inclination and eccentricity phase space is an exceedingly complicated web-like structure of lunisolar secular resonances, which become particularly dense near the inclinations of the navigation satellite orbits. A clear picture of the physical significance of these resonances is of considerable practical interest for the design of disposal strategies for the four constellations. Here we present analytical and semi-analytical models that accurately reflect the true nature of the resonant interactions, and trace the topological organization of the manifolds on which the chaotic motions take place. We present an atlas of FLI stability maps, showing the extent of the chaotic regions of the phase space, computed through a hierarchy of more realistic, and more complicated, models, and compare the chaotic zones in these charts with the analytical estimation of the width of the chaotic layers from the heuristic Chirikov resonance-overlap criterion. As the semi-major axis of the satellite is receding, we observe a transition from stable Nekhoroshev-like structures at three Earth radii, where regular orbits dominate, to a Chirikov regime where resonances overlap at five Earth radii. From a numerical estimation of the Lyapunov times, we find that many of the inclined, nearly circular orbits of the navigation satellites are strongly chaotic and that their dynamics are unpredictable on decadal timescales.