We study the presence of resonances in the region of space around the Earth. We consider a massless body (e.g, a dust particle or a small space debris) subject to different forces: the gravitational attraction of the geopotential, the effects of Sun and Moon. We distinguish different types of resonances: tesseral resonances are due to a commensurability involving the revolution of the particle and the rotation of the Earth, semi-secular resonances include the rates of variation of the mean anomalies of Moon and Sun, while secular resonances just depend on the rates of variation of the arguments of perigee and the longitudes of the ascending nodes of the perturbing bodies. We characterize such resonances, giving precise statements on the regions where the resonances can be found and provide examples of some specific commensurability relations.

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.

We consider the three-body mean motion resonance defined by the Jovian moons Io, Europa, and Ganymede, which is commonly known as the Laplace resonance. In particular, we construct approximate models for the evolution of the librating argument over the period of 100 years, focusing on its principal amplitude and frequency, and on the observed mean motion combinations associated with the quasi-resonant interactions. First, we numerically propagated the Cartesian equations of motion of the Jovian system for the period under examination, and by comparing the results with a suitable set of ephemerides, we derived the main dynamical effects on the target quantities. Using these effects, we built an alternative Hamiltonian formulation and used the normal forms theory to locate the resonance and to compute its main amplitude and frequency. From the Cartesian model we observe that on the timescale considered and with ephemerides as initial conditions, both the librating argument and the diagnostics are well approximated by considering the mutual gravitational interactions of Jupiter and the Galilean moons (including Callisto), and the effect of Jupiter's J2 harmonic. Under the same initial conditions, the Hamiltonian formulation in which Callisto and J2 are reduced to their secular contributions achieves larger errors for the quantities above, particularly for the librating argument. By introducing appropriate resonant variables, we show that these errors can be reduced by moving in a certain action-angle phase plane, which in turn implies the necessity of a tradeoff in the selection of the initial conditions. In addition to being a good starting point for a deeper understanding of the Laplace resonance, the models and methods described are easily generalizable to different types of multi-body mean motion resonances. They are also prime tools for studying the dynamics of extrasolar systems.

We study the dynamics of the de Sitter resonance, namely the stable equilibrium configuration of the first three Galilean satellites. We clarify the relation between this family of configurations and the more general Laplace resonant states. In order to describe the dynamics around the de Sitter stable equilibrium, a one-degree of freedom Hamiltonian normal form is constructed and exploited to identify initial conditions leading to the two families. The normal form Hamiltonian is used to check the accuracy in the location of the equilibrium positions. Besides, it gives a measure of how sensitive it is with respect to the different perturbations acting on the system. By looking at the phase-plane of the normal form, we can identify a \sl Laplace-like \rm configuration, which highlights many substantial aspects of the observed one.

The nearby space surrounding the Earth is densely populated by artificial satellites and instruments, whose orbits are distributed within the Low-Earth-Orbit region (LEO), ranging between 90 and 2 000 $km$ of altitude. As a consequence of collisions and fragmentations, many space debris of different sizes are left in the LEO region. Given the threat raised by the possible damages which a collision of debris can provoke with operational or manned satellites, the study of their dynamics is nowadays mandatory. This work is focused on the existence of equilibria and the dynamics of resonances in LEO. We base our results on a simplified model which includes the geopotential and the atmospheric drag. Using such model, we make a qualitative study of the resonances and the equilibrium positions, including their location and stability. The dissipative effect due to the atmosphere provokes a tidal decay, but we give examples of different behaviors, precisely a straightforward passage through the resonance or rather a temporary capture. We also investigate the effect of the solar cycle which is responsible of fluctuations of the atmospheric density and we analyze the influence of Sun and Moon on LEO objects.

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.

The increasing threat raised by space debris led to the development of different mathematical models and approaches to investigate the dynamics of small particles orbiting around the Earth. Such models and methods strongly depend on the altitude of the objects above Earth's surface, since the strength of the different forces acting on an Earth orbiting object (geopotential, atmospheric drag, lunar and solar attractions, solar radiation pressure, etc.) varies with the altitude of the debris. In this review, our focus is on presenting different analytical and numerical approaches employed in modern studies of the space debris problem. We start by considering a model including the geopotential, solar and lunar gravitational forces and the solar radiation pressure. We summarize the equations of motion using different formalisms: Cartesian coordinates, Hamiltonian formulation using Delaunay and epicyclic variables, Milankovitch elements. Some of these methods lead in a straightforward way to the analysis of resonant motions. In particular, we review results found recently about the dynamics near tesseral, secular and semi-secular resonances. As an application of the above methods, we proceed to analyze a timely subject namely the possible causes for the onset of chaos in space debris dynamics. Precisely, we discuss the phenomenon of overlapping of resonances, the effect of a large area-to-mass ratio, the influence of lunisolar secular resonances. We conclude with a short discussion about the effect of the dissipation due to the atmospheric drag and we provide a list of minor effects, which could influence the dynamics of space debris.

The long-term dynamics of the geostationary Earth orbits (GEO) is revisited through the application of canonical perturbation theory. We consider a Hamiltonian model accounting for all major perturbations: geopotential at order and degree two, lunisolar perturbations with a realistic model for the Sun and Moon orbits, and solar radiation pressure. The long-term dynamics of the GEO region has been studied both numerically and analytically, in view of the relevance of such studies to the issue of space debris or to the disposal of GEO satellites. Past studies focused on the orbital evolution of objects around a nominal solution, hereafter called the forced equilibrium solution, which shows a particularly strong dependence on the area-to-mass ratio. Here, we i) give theoretical estimates for the long-term behavior of such orbits, and ii) we examine the nature of the forced equilibrium itself. In the lowest approximation, the forced equilibrium implies motion with a constant non-zero average `forced eccentricity', as well as a constant non-zero average inclination, otherwise known in satellite dynamics as the inclination of the invariant `Laplace plane'. Using a higher order normal form, we demonstrate that this equilibrium actually represents not a point in phase space, but a trajectory taking place on a lower-dimensional torus. We give analytical expressions for this special trajectory, and we compare our results to those found by numerical orbit propagation. We finally discuss the use of proper elements, i.e., approximate integrals of motion for the GEO orbits.

We provide a detailed derivation of the analytical expansion of the lunar and solar disturbing functions. Although there exist several papers on this topic, many derivations contain mistakes in the final expansion or rather (just) in the proof, thereby necessitating a recasting and correction of the original derivation. In this work, we provide a self-consistent and definite form of the lunisolar expansion. We start with Kaula's expansion of the disturbing function in terms of the equatorial elements of both the perturbed and perturbing bodies. Then we give a detailed proof of Lane's expansion, in which the elements of the Moon are referred to the ecliptic plane. Using this approach the inclination of the Moon becomes nearly constant, while the argument of perihelion, the longitude of the ascending node, and the mean anomaly vary linearly with time. We make a comparison between the different expansions and we profit from such discussion to point out some mistakes in the existing literature, which might compromise the correctness of the results. As an application, we analyze the long--term motion of the highly elliptical and critically--inclined Molniya orbits subject to quadrupolar gravitational interactions. The analytical expansions presented herein are very powerful with respect to dynamical studies based on Cartesian equations, because they quickly allow for a more holistic and intuitively understandable picture of the dynamics.

We analyze the combined effect of Poynting-Robertson and solar wind drag on space debris. We derive a model within Cartesian, Gaussian and Hamiltonian frameworks. We focus on the geosynchronous resonance, although the results can be easily generalized to any resonance. By numerical and analytical techniques, we compute the drift in semi-major axis due to Poynting-Robertson and solar wind drag. After a linear stability analysis of the equilibria, we combine a careful investigation of the regular, resonant, chaotic behavior of the phase space with a long-term propagation of a sample of initial conditions. The results strongly depend on the value of the area-to-mass ratio of the debris, which might show different dynamical behaviors: temporary capture or escape from the geosynchronous resonance, as well as temporary capture or escape from secondary resonances involving the rate of variation of the longitude of the Sun. Such analysis shows that Poynting-Robertson and solar wind drag must be taken into account, when looking at the long-term behavior of space debris. Trapping or escape from the resonance can be used to place the debris in convenient regions of the phase space.

Using bifurcation theory, we study the secular resonances induced by the Sun and Moon on space debris orbits around the Earth. In particular, we concentrate on a special class of secular resonances, which depend only on the debris' orbital inclination. This class is typically subdivided into three distinct types of secular resonances: those occurring at the critical inclination, those corresponding to polar orbits, and a third type resulting from a linear combination of the rates of variation of the argument of perigee and the longitude of the ascending node. The model describing the dynamics of space debris includes the effects of the geopotential, as well as the Sun's and Moon's attractions, and it is defined in terms of suitable action-angle variables. We consider the system averaged over both the mean anomaly of the debris and those of the Sun and Moon. Such multiply-averaged Hamiltonian is used to study the lunisolar resonances which depend just on the inclination. Borrowing the technique from the theory of bifurcations of Hamiltonian normal forms, we study the birth of periodic orbits and we determine the energy thresholds at which the bifurcations of lunisolar secular resonances take place. This approach gives us physically relevant information on the existence and location of the equilibria, which help us to identify stable and unstable regions in the phase space. Besides their physical interest, the study of inclination dependent resonances offers interesting insights from the dynamical point of view, since it sheds light on different phenomena related to bifurcation 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 dynamics of small bodies around the Earth has gained a renewed interest, since the awareness of the problems that space debris can cause in the nearby future. A relevant role in space debris is played by lunisolar secular resonances, which might contribute to an increase of the orbital elements, typically of the eccentricity. We concentrate our attention on the lunisolar secular resonance described by the relation $2\dot{\omega}+\dot{\Omega}=0$, where $\omega$ and $\Omega$ denote the argument of perigee and the longitude of the ascending node of the space debris. We introduce three different models with increasing complexity. We show that the growth in eccentricity, as observed in space debris located in the MEO region at the inclination about equal to $56^\circ$, can be explained as a natural effect of the secular resonance $2\dot{\omega}+\dot{\Omega}=0$, while the chaotic variations of the orbital parameters are the result of interaction and overlapping of nearby resonances.

We study the dynamics of the space debris in regions corresponding to minor resonances; precisely, we consider the resonances 3:1, 3:2, 4:1, 4:3, 5:1, 5:2, 5:3, 5:4, where a j:l resonance (with j, l integers) means that the periods of revolution of the debris and of rotation of the Earth are in the ratio j/l. We consider a Hamiltonian function describing the effect of the geopotential and we use suitable finite expansions of the Hamiltonian for the description of the different resonances. In particular, we determine the leading terms which dominate in a specific orbital region, thus limiting our computation to very few harmonics. Taking advantage from the pendulum-like structure associated to each term of the expansion, we are able to determine the amplitude of the islands corresponding to the different harmonics. By means of simple mathematical formulae, we can predict the occurrence of splitting or overlapping of the resonant islands for different values of the parameters. We also find several cases which exhibit a transcritical bifurcation as the inclination is varied. These results, which are based on a careful mathematical analysis of the Hamiltonian expansion, are confirmed by a numerical study of the dynamical behavior obtained by computing the so-called Fast Laypunov Indicators. Since the Hamiltonian approach includes just the effect of the geopotential, we validate our results by performing a numerical integration in Cartesian variables of a more complete model including the gravitational attraction of Sun and Moon, as well as the solar radiation pressure.

We investigate the dynamics of satellites and space debris in external resonances, namely in the region outside the geostationary ring. Precisely, we focus on the 1:2, 1:3, 2:3 resonances, which are located at about 66 931.4 km, 87 705.0 km, 55 250.7 km, respectively. Some of these resonances have been already exploited in space missions, like XMM-Newton and Integral. Our study is mainly based on a Hamiltonian approach, which allows us to get fast and reliable information on the dynamics in the resonant regions. Significative results are obtained even by considering just the effect of the geopotential in the Hamiltonian formulation. For objects (typically space debris) with high area-to-mass ratio the Hamiltonian includes also the effect of the solar radiation pressure. In addition, we perform a comparison with the numerical integration in Cartesian variables, including the geopotential, the gravitational attraction of Sun and Moon, and the solar radiation pressure. We implement some simple mathematical tools that allows us to get information on the terms which are dominant in the Fourier series expansion of the Hamiltonian around a given resonance, on the amplitude of the resonant islands and on the location of the equilibrium points. We also compute the Fast Lyapunov Indicators, which provide a cartography of the resonant regions, yielding the main dynamical features associated to the external resonances. We apply these techniques to analyze the 1:2, 1:3, 2:3 resonances; we consider also the case of objects with large area-to-mass ratio and we provide an application to the case studies given by XMM-Newton and Integral.

We investigate the stability of motion close to the Lagrangian equilibrium points L4 and L5 in the framework of the spatial, elliptic, restricted three- body problem, subject to the radial component of Poynting-Robertson drag. For this reason we develop a simplified resonant model, that is based on averaging theory, i.e. averaged over the mean anomaly of the perturbing planet. We find temporary stability of particles displaying a tadpole motion in the 1:1 resonance. From the linear stability study of the averaged simplified resonant model, we find that the time of temporary stability is proportional to beta a1 n1 , where beta is the ratio of the solar radiation over the gravitational force, and a1, n1 are the semi-major axis and the mean motion of the perturbing planet, respectively. We extend previous results (Murray (1994)) on the asymmetry of the stability indices of L4 and L5 to a more realistic force model. Our analytical results are supported by means of numerical simulations. We implement our study to Jupiter-like perturbing planets, that are also found in extra-solar planetary systems.

In a dissipative system the time to reach an attractor is often influenced by the peculiarities of the model and in particular by the strength of the dissipation. In particular, as a dissipative model we consider the spin-orbit problem providing the dynamics of a triaxial satellite orbiting around a central planet and affected by tidal torques. The model is ruled by the oblateness parameter of the satellite, the orbital eccentricity, the dissipative parameter and the drift term. We devise a method which provides a reliable indication on the transient time which is needed to reach an attractor in the spin-orbit model; the method is based on an analytical result, precisely a suitable normal form construction. This method provides also information about the frequency of motion. A variant of such normal form used to parametrize invariant attractors provides a specific formula for the drift parameter, which in turn yields a constraint - which might be of interest in astronomical problems - between the oblateness of the satellite and its orbital eccentricity.