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 investigate the dynamics of charged dust close to outer mean-motion resonances with planet Jupiter. The importance of the interplanetary magnetic field on the orbital evolution of dust is clearly demonstrated. New dynamical phenomena are found that do not exist in the classical problem of uncharged dust. We find changes in the orientation of the orbital planes of dust particles, an increased amount of chaotic orbital motions, sudden 'jumps' in the resonant argument, and a decrease in time of temporary capture due to the Lorentz force. Variations in the orbital planes of dust grain orbits are found to be related to the angle between the orbital angular momentum and magnetic axes of the heliospheric field and the rotation rate of the Sun. These variations are bound using a simplified model derived from the full dynamical problem using first order averaging theory. It is found that the interplanetary magnetic field does not affect the capture process, that is still dominated by the other non-gravitational forces. Our study is based on a dynamical model in the framework of the inclined circular restricted three-body problem. Additional forces include solar radiation pressure, solar wind drag, the Poynting-Robertson effect, and the influence of a Parker spiral type interplanetary magnetic field model. The analytical estimates are derived on the basis of Gauss' form of planetary equations of motion. Numerical results are obtained by simulations of dust grain orbits together with the system of variational equations. Chaotic regions in phase space are revealed by means of Fast Lyapunov Chaos Indicators.

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.

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.

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.