Recent works demonstrated that the dynamics caused by the planetary oblateness coupled with the solar radiation pressure can be described through a model based on singly-averaged equations of motion. The coupled perturbations affect the evolution of the eccentricity, inclination and orientation of the orbit with respect to the Sun--Earth line. Resonant interactions lead to non-trivial orbital evolution that can be exploited in mission design. Moreover, the dynamics in the vicinity of each resonance can be analytically described by a resonant model that provides the location of the central and hyperbolic invariant manifolds which drive the phase space evolution. The classical tools of the dynamical systems theory can be applied to perform a preliminary mission analysis for practical applications. On this basis, in this work we provide a detailed derivation of the resonant dynamics, also in non-singular variables, and discuss its properties, by studying the main bifurcation phenomena associated to each resonance. Last, the analytical model will provide a simple analytical expression to obtain the area-to-mass ratio required for a satellite to deorbit from a given altitude in a feasible timescale.

Navigation satellites are known from numerical studies to reside in a dynamically sensitive environment, which may be of profound importance for their long-term sustainability. We derive the fundamental Hamiltonian of GNSS dynamics and show analytically that near-circular trajectories lie in the neighborhood of a Normally Hyperbolic Invariant Manifold (NHIM), which is the primary source of hyperbolicity. Quasi-circular orbits escape through chaotic transport, regulated by the NHIM's stable and unstable manifolds, following a power-law escape time distribution $P(t) \sim t^{-\alpha}$, with $\alpha \sim 0.8 - 1.5$. Our study is highly relevant for the design of satellite disposal trajectories, using manifold dynamics.

We provide an analytical approximation to the dynamics in each of the three most important low order secondary resonances (1:1, 2:1, and 3:1) bifurcating from the synchronous primary resonance in the gravitational spin-orbit problem. To this end we extend the perturbative approach introduced in Gkolias et. al. (2016), based on normal form series computations. This allows to recover analytically all non-trivial features of the phase space topology and bifurcations associated with these resonances. Applications include the characterization of spin states of irregular planetary satellites or double systems of minor bodies with irregular shapes. The key ingredients of our method are: i) the use of a detuning parameter measuring the distance from the exact resonance, and ii) an efficient scheme to `book-keep' the series terms, which allows to simultaneously treat all small parameters entering the problem. Explicit formulas are provided for each secondary resonance, yielding i) the time evolution of the spin state, ii) the form of phase portraits, iii) initial conditions and stability for periodic solutions, and iv) bifurcation diagrams associated with the periodic orbits. We give also error estimates of the method, based on analyzing the asymptotic behavior of the remainder of the normal form series.

In this work the orbital dynamics of Earth satellites about the geosynchronous altitude are explored, with primary goal to assess current mitigation guidelines as well as to discuss the future exploitation of the region. A thorough dynamical mapping was conducted in a high-definition grid of orbital elements, enabled by a fast and accurate semi-analytical propagator, which considers all the relevant perturbations. The results are presented in appropriately selected stability maps to highlight the underlying mechanisms and their interplay, that can lead to stable graveyard orbits or fast re-entry pathways. The natural separation of the long-term evolution between equatorial and inclined satellites is discussed in terms of post-mission disposal strategies. Moreover, we confirm the existence of an effective cleansing mechanism for inclined geosynchronous satellites and discuss its implications in terms of current guidelines as well as alternative mission designs that could lead to a sustainable use of the geosynchronous orbital region.

The slow deformation of terrestrial orbits in the medium range, subject to lunisolar resonances, is well approximated by a family of Hamiltonian flow with $2.5$ degree-of-freedom. The action variables of the system may experience chaotic variations and large drift that we may quantify. Using variational chaos indicators, we compute high-resolution portraits of the action space. Such refined meshes allow to reveal the existence of tori and structures filling chaotic regions. Our elaborate computations allow us to isolate precise initial conditions near specific zones of interest and study their asymptotic behaviour in time. Borrowing classical techniques of phase- space visualisation, we highlight how the drift is mediated by the complement of the numerically detected KAM tori.

In a realistic scenario, the evolution of the rotational dynamics of a celestial or artificial body is subject to dissipative effects. Time-varying non-conservative forces can be due to, for example, a variation of the moments of inertia or to tidal interactions. In this work, we consider a simplified model describing the rotational dynamics, known as the spin-orbit problem, where we assume that the orbital motion is provided by a fixed Keplerian ellipse. We consider different examples in which a non-conservative force acts on the model and we propose an analytical method, which reduces the system to a Hamiltonian framework. In particular, we compute a time parametrisation in a series form, which allows us to transform the original system into a Hamiltonian one. We also provide applications of our method to study the rotational motion of a body with time-varying moments of inertia, e.g. an artificial satellite with flexible components, as well as subject to a tidal torque depending linearly on the velocity.

We consider Earth satellite orbits in the range of semi-major axes where the perturbing effects of Earth's oblateness and lunisolar gravity are of comparable order. This range covers the medium-Earth orbits (MEO) of the Global Navigation Satellite Systems and the geosynchronous orbits (GEO) of the communication satellites. We recall a secular and quadrupolar model, based on the Milankovitch vector formulation of perturbation theory, which governs the long-term orbital evolution subject to the predominant gravitational interactions. We study the global dynamics of this two-and-a-half degrees-of-freedom Hamiltonian system by means of the fast Lyapunov indicator (FLI), used in a statistical sense. Specifically, we characterize the degree of chaoticity of the action space using angle-averaged normalized FLI maps, thereby overcoming the angle dependencies of the conventional stability maps. Emphasis is placed upon the phase-space structures near secular resonances, which are of first importance to the space debris community. We confirm and quantify the transition from order to chaos in MEO, stemming from the critical inclinations, and find that highly inclined GEO orbits are particularly unstable. Despite their reputed normality, Earth satellite orbits can possess an extraordinarily rich spectrum of dynamical behaviors, and, from a mathematical perspective, have all the complications that make them very interesting candidates for testing the modern tools of chaos theory.

We study the resonant dynamics in a simple one degree of freedom, time dependent Hamiltonian model describing spin-orbit interactions. The equations of motion admit periodic solutions associated with resonant motions, the most important being the synchronous one in which most evolved satellites of the Solar system, including the Moon, are observed. Such primary resonances can be surrounded by a chain of smaller islands which one refers to as secondary resonances. Here, we propose a novel canonical normalization procedure allowing to obtain a higher order normal form, by which we obtain analytical results on the stability of the primary resonances as well as on the bifurcation thresholds of the secondary resonances. The procedure makes use of the expansion in a parameter, called the detuning, measuring the shift from the exact secondary resonance. Also, we implement the so-called `book-keeping' method, i.e., the introduction of a suitable separation of the terms in orders of smallness in the normal form construction, which deals simultaneously with all the small parameters of the problem. Our analytical computation of the bifurcation curves is in excellent agreement with the results obtained by a numerical integration of the equations of motion, thus providing relevant information on the parameter regions where satellites can be found in a stable configuration.

The motion of a satellite around a planet can be studied by the Hill model, which is a modification of the restricted three body problem pertaining to motion of a satellite around a planet. Although the dynamics of the circular Hill model have been extensively studied in the literature, only few results about the dynamics of the elliptic model were known up to now, namely the equations of motion and few unstable families of periodic orbits. In the present study we extend these results by computing a large set of families of periodic orbits and their linear stability and classify them according to their resonance condition. Although most of them are unstable, we were able to find a considerable number of stable ones. By computing appropriate maps of dynamical stability, we study the effect of the planetary eccentricity on the stability of satellite orbits. We see that, even for large values of the planetary eccentricity, regular orbits can be found in the vicinity of stable periodic orbits. The majority of irregular orbits are escape orbits.