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.

The aim of this work is to analyze the orbital evolution of the mean eccentricity given by the Two-Line Elements (TLE) set of the Molniya satellites constellation. The approach is bottom-up, aiming at a synergy between the observed dynamics and the mathematical modeling. Being the focus the long-term evolution of the eccentricity, the dynamical model adopted is a doubly-averaged formulation of the third-body perturbation due to Sun and Moon, coupled with the oblateness effect on the orientation of the satellite. The numerical evolution of the eccentricity, obtained by a two-degree-of-freedom model assuming different orders in the series expansion of the third-body effect, is compared against the mean evolution given by the TLE. The results show that, despite being highly elliptical orbits, the second order expansion catches extremely well the behavior. Also, the lunisolar effect turns out to be non-negligible for the behavior of the longitude of the ascending node and the argument of pericenter. The role of chaos in the timespan considered is also addressed. Finally, a frequency series analysis is proposed to show the main contributions that can be detected from the observational data.

The icy moons are in the focus of the exploration plans of the leading space agencies because of the indications of water-based life and geological activity observed in a number of these objects. In particular, the presence of geyser-like jets of water near Enceladus' south pole has turned this moon of Saturn into a priority candidate to search for life and habitability features. This investigation proposes a set of trajectories between Halo orbits about Lagrangian points L1 and L2 in the Saturn-Enceladus Circular Restricted Three-Body Problem as science orbits for a future in situ mission at Enceladus. The design methodology is presented, followed by the analysis of the observational performance of the solutions. The conclusion is that the proposed orbits exhibit suitable features for their use in the scientific exploration of Enceladus, i.e., long transfer times, low altitudes, wide surface visibility windows and long times of overflight.

This contribution investigates the properties of a category of orbits around Enceladus. The motivation is the interest in the in situ exploration of this moon following the detection on behalf of Cassini of plumes of water and organic compounds close to its south pole. In a previous investigation, a set of heteroclinic transfers were designed between Halo orbits around the equilibrium points L1 and L2 of the circular restricted three-body problem with Saturn and Enceladus as primaries. The kinematical and geometrical characteristics of those trajectories makes them good candidates as science orbits for the extended observation of the surface of Enceladus: they are highly inclined, they approach the moon and they are maneuver free. However, the low heights above the surface and the strong perturbing effect of Saturn impose a more careful look at their dynamics, in particular regarding the influence of the polar flattening of the primaries. Therefore, those solutions are here reconsidered by employing a dynamical model that includes the effect of the oblateness of Saturn and Enceladus, individually and in combination. Substitutes of the Halo orbits around the equilibrium points L1 and L2 and their stable and unstable hyperbolic invariant manifolds are obtained in the perturbed models, and maneuver-free heteroclinic transfers are identified in the new framework. A systematic comparison with the corresponding solutions of the unperturbed problem shows that qualitative and quantitative features are not significantly altered when the oblateness of the primaries is taken into account, and that J2 of Saturn plays a larger role than the oblateness of Enceladus. From a mission perspective, the results confirm the scientific value of the solutions obtained in the classical circular restricted three-body problem and suggests that this simpler model can be used in a preliminary feasibility analysis.

The aim of this work is to provide an analytical model to characterize the equilibrium points and the phase space associated with the singly-averaged dynamics caused by the planetary oblateness coupled with the solar radiation pressure perturbations. A two-dimensional differential system is derived by considering the classical theory, supported by the existence of an integral of motion comprising semi-major axis, eccentricity and inclination. Under the single resonance hypothesis, the analytical expressions for the equilibrium points in the eccentricity-resonant angle space are provided, together with the corresponding linear stability. The Hamiltonian formulation is also given. The model is applied considering, as example, the Earth as major oblate body, and a simple tool to visualize the structure of the phase space is presented. Finally, some considerations on the possible use and development of the proposed model are drawn.

We present the main findings of a dynamical mapping performed in the Low Earth Orbit region. The results were obtained by propagating an extended grid of initial conditions, considering two different epochs and area-to-mass ratios, by means of a singly averaged numerical propagator. It turns out that dynamical resonances associated with high-degree geopotential harmonics, lunisolar perturbations and solar radiation pressure can open natural deorbiting highways. For area-to-mass ratios typical of the orbiting intact objects, these corridors can be exploited only in combination with the action exerted by the atmospheric drag. For satellites equipped with an area augmentation device, we show the boundary of application of the drag, and where the solar radiation pressure can be exploited.

The aim of this work is to highlight the crucial role that orbital resonances associated with solar radiation pressure can have in Low Earth Orbit. We review the corresponding literature, and provide an analytical tool to estimate the maximum eccentricity which can be achieved for well-defined initial conditions. We then compare the results obtained with the simplified model with the results obtained with a more comprehensive dynamical model. The analysis has important implications both from a theoretical point of view, because it shows that the role of some resonances was underestimated in the past, but also from a practical point of view in the perspective of passive deorbiting solutions for satellites at the end-of-life.

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.

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.

Numerical simulations carried out over the past decade suggest that the orbits of the Global Navigation Satellite Systems are unstable, resulting in an apparent chaotic growth of the eccentricity. Here we show that the irregular and haphazard character of these orbits reflects a similar irregularity in the orbits of many celestial bodies in our Solar System. We find that secular resonances, involving linear combinations of the frequencies of nodal and apsidal precession and the rate of regression of lunar nodes, occur in profusion so that the phase space is threaded by a devious stochastic web. As in all cases in the Solar System, chaos ensues where resonances overlap. These results may be significant for the analysis of disposal strategies for the four constellations in this precarious region of space.

It is known that most of the craters on the surface of the Moon were created by the collision of minor bodies of the Solar System. Main Belt Asteroids, which can approach the terrestrial planets as a consequence of different types of resonance, are actually the main responsible for this phenomenon. Our aim is to investigate the impact distributions on the lunar surface that low-energy dynamics can provide. As a first approximation, we exploit the hyberbolic invariant manifolds associated with the central invariant manifold around the equilibrium point L_2 of the Earth - Moon system within the framework of the Circular Restricted Three - Body Problem. Taking transit trajectories at several energy levels, we look for orbits intersecting the surface of the Moon and we attempt to define a relationship between longitude and latitude of arrival and lunar craters density. Then, we add the gravitational effect of the Sun by considering the Bicircular Restricted Four - Body Problem. As further exploration, we assume an uniform density of impact on the lunar surface, looking for the regions in the Earth - Moon neighbourhood these colliding trajectories have to come from. It turns out that low-energy ejecta originated from high-energy impacts are also responsible of the phenomenon we are considering.