Averaging on Earth‐Crossing Orbits

The orbits of planetcrossing asteroids (and comets) can undergo close approaches and collisions with some major planet. This introduces a singularity in the Nbody Hamiltonian, and the averaging of the equations of motion, traditionally used to compute secular perturbations, is undefined. We show that it is possible to define in a rigorous way some generalised averaged equations of motion, in such a way that the generalised solutions are unique and piecewise smooth. This is obtained, both in the planar and in the threedimensional case, by means of the method of extraction of the singularities by Kantorovich. The modified distance used to approximate the singularity is the one used by Wetherill in his method to compute probability of collision. Some examples of averaged dynamics have been computed; a systematic exploration of the averaged phase space to locate the secular resonances should be the next step.'Alice sighed wearily. I think you might do something better with the time she said, than waste it asking riddles with no answers(Alice in Wonderland, L. Carroll)     

“Bumerang” Must Return to Earth

Solar System Research - The “Bumerang” space mission to the satellites of Mars is discussed. The purpose of the mission is to study the satellites and return of soil samples from Phobos...     

Earth Tides and Lense–Thirring Effect

The general relativistic Lense—Thirring effect can be measured by inspecting a suitable combination of the orbital residuals of the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II. The solid and ocean Earth tides affect the recovery of the parameter by means of which the gravitomagnetic signal is accounted for in the combined residuals. Thus an extensive analysis of the perturbations induced on these orbital elements by the solid and ocean Earth tides is carried out. It involves the l=2 terms for the solid tides and the l=2,3,4 terms for the ocean tides. The perigee of LAGEOS II turns out to be very sensitive to the l=3 part of the ocean tidal spectrum, contrary to the nodes of LAGEOS and LAGEOS II. The uncertainty in the solid tidal perturbations, mainly due to the Love number k ₂, ranges from 0.4% to 1.5%, while the ocean tides are uncertain at 5–15% level. The obtained results are used in order to check in a preliminary way which tidal constituents the Lense-Thirring shift is sensitive to. In particular it is tested if the semisecular 18.6-year zonal tide really does not affect the combined residuals. It turns out that, if modeled at the level of accuracy worked out in the paper, the l=2,4 m=0 and also, to a lesser extent, the l=3, m=0 tidal perturbations cancel out.This revised version was published online in October 2005 with corrections to the Cover Date.     

Asymmetric effects on Earth’s polar motion

Differential equations ruling the Earth’s polar motion are slightly asymmetric with respect to the pole coordinates. This is not only associated with the lack of axial symmetry around the Earth figure axis (triaxiality) but also with the longitude dependency of the pole tide (the main contribution). We propose a consistent handling of both asymmetric contributions, formulating a unique equation in the complex equatorial plane, of which we derive a general solution. Difference with respect to the usual symmetric solution is discussed and found significant in light of the present accuracy of the observed pole coordinates. For the same geophysical excitation, the prograde Chandler wobble is accompanied by a retrograde component up to 2 milliarcseconds (mas), transforming it in a slight elliptic motion. The asymmetric contribution is relatively larger in the geodetic excitation function, for Chandler wobble excitation mixes prograde and retrograde components of comparable level (1 mas).     

Fast Earth–Moon transfers with ballistic capture

This contribution deals with fast Earth–Moon transfers with ballistic capture in the patched three-body model. We compute ensembles of preliminary solutions using a model that takes into account the relative inclination of the orbital planes of the primaries. The ballistic capture orbits around the Moon are obtained relying on the hyperbolic invariant structures associated to the collinear Lagrangian points of the Earth–Moon system, and the Sun–Earth system portion of the transfers are quasi-periodic orbits obtained by a genetic algorithm. The trajectories are designed to be good initial guesses to search optimal cost-efficient short-time Earth–Moon transfers with ballistic capture in more realistic models.     

Transfers from the Earth to $$L_2$$ L 2 Halo orbits in the Earth–Moon bicircular problem

Celestial Mechanics and Dynamical Astronomy - The Cassini spacecraft discovered many close-in small satellites in Saturnian system, and several of them exhibit exotic orbital states due to...     

Extension of fast periodic transfer orbits from the Earth–Moon RTBP to the Sun–Earth–Moon Quasi-Bicircular Problem

Starting from 80 families of low-energy fast periodic transfer orbits in the Earth–Moon planar circular Restricted Three Body Problem (RTBP), we obtain by analytical continuation 11 periodic orbits and 25 periodic arcs with similar properties in the Sun–Earth–Moon Quasi-Bicircular Problem (QBCP). A novel and very simple procedure is introduced giving the solar phases at which to attempt continuation. Detailed numerical results for each periodic orbit and arc found are given, including their stability parameters and minimal distances to the Earth and Moon. The periods of these orbits are between 2.5 and 5 synodic months, their energies are among the lowest possible to achieve an Earth–Moon transfer, and they show a diversity of circumlunar trajectories, making them good candidates for missions requiring repeated passages around the Earth and the Moon with close approaches to the last.     

On constructing the general Earth’s rotation theory

In the present paper the equations of the orbital motion of the major planets and the Moon and the equations of the three–axial rigid Earth’s rotation in Euler parameters are reduced to the secular system describing the evolution of the planetary and lunar orbits (independent of the Earth’s rotation) and the evolution of the Earth’s rotation (depending on the planetary and lunar evolution). Hence, the theory of the Earth’s rotation can be presented by means of the series in powers of the evolutionary variables with quasi-periodic coefficients with respect to the planetary–lunar mean longitudes. This form of the Earth’s rotation problem is compatible with the general planetary theory involving the separation of the short–period and long–period variables and avoiding the appearance of the non–physical secular terms.     

Earth–Mars transfers with ballistic escape and low-thrust capture

In this paper novel Earth–Mars transfers are presented. These transfers exploit the natural dynamics of n-body models as well as the high specific impulse typical of low-thrust systems. The Moon-perturbed version of the Sun–Earth problem is introduced to design ballistic escape orbits performing lunar gravity assists. The ballistic capture is designed in the Sun–Mars system where special attainable sets are defined and used to handle the low-thrust control. The complete trajectory is optimized in the full n-body problem which takes into account planets’ orbital inclinations and eccentricities. Accurate, efficient solutions with reasonable flight times are presented and compared with known results.     

Do planetary encounters reset surfaces of near Earth asteroids?

Processes such as the solar wind sputtering and micrometeorite impacts can modify optical properties of surfaces of airless bodies. This explains why spectra of the main belt asteroids, exposed to these ‘space weathering’ processes over eons, do not match the laboratory spectra of ordinary chondrite (OC) meteorites. In contrast, an important fraction of Near Earth Asteroids (NEAs), defined as Q-types in the asteroid taxonomy, display spectral attributes that are a good match to OCs. Here we study the possibility that the Q-type NEAs underwent recent encounters with the terrestrial planets and that the tidal gravity (or other effects) during these encounters exposed fresh OC material on the surface (thus giving it the Q-type spectral properties). We used numerical integrations to determine the statistics of encounters of NEAs to planets. The results were used to calculate the fraction and orbital distribution of Q-type asteroids expected in the model as a function of the space weathering timescale, t_sw (see main text for definition), and maximum distance, r^∗, at which planetary encounters can reset the surface. We found that t_sw ∼ 10⁶ yr (at 1 AU) and r^∗ ∼ 5R_pl, where R_pl is the planetary radius, best fit the data. Values t_sw < 10⁵ yr would require that r^∗ > 20R_pl, which is probably implausible because these very distant encounters should be irrelevant. Also, the fraction of Q-type NEAs would be probably much larger than the one observed if t_sw > 10⁷ yr. We found that t_sw ∝ q², where q is the perihelion distance, expected if the solar wind sputtering controls t_sw, provides a better match to the orbital distribution of Q-type NEAs than models with fixed t_sw. We also discuss how the Earth magnetosphere and radiation effects such as YORP can influence the spectral properties of NEAs.     

Earth–Mars halo to halo low thrust manifold transfers

Application of low thrust propulsion to interconnect ballistic trajectories on invariant manifolds associated with multiple circular restricted three body systems has been investigated. Sun-planet three body models have been coupled to compute the two ballistic trajectories, where electric propulsion is used to interconnect these trajectories as no direct intersection in the Poincarè sections exists. The ability of a low thrust to provide the energy change required to transit the spacecraft between two systems has been assessed for some Earth to Mars transfers. The approach followed consists in a planetary escape on the unstable manifold starting from a periodic orbit around one of the two collinear libration points near the secondary body. Following the planetary escape and the subsequent coasting phase, the electric thruster is activated and executes an ad-hoc thrusting phase. The complete transfer design, composed of the three discussed phases, and possible applications to Earth–Mars missions is developed where the results are outlined in this paper.     

Probing the Earth’s Interior with the LENA Detector

A future large-volume liquid scintillator detector such as the proposed 50 kton LENA (Low Energy Neutrino Astronomy) detector would provide a high-statistics measurement of terrestrial antineutrinos originating from β-decays of the uranium and thorium chains. Additionally, the neutron is scattered in the forward direction in the detection reaction . Henceforth, we investigate to what extent LENA can distinguish between certain geophysical models on the basis of the angular dependence of the geoneutrino flux. Our analysis is based on a Monte-Carlo simulation with different levels of light yield, considering an unloaded PXE scintillator. We find that LENA is able to detect deviations from isotropy of the geoneutrino flux with high significance. However, if only the directional information is used, the time required to distinguish between different geophysical models is of the order of severals decades. Nonetheless, a high-statistics measurement of the total geoneutrino flux and its spectrum still provides an extremely useful glance at the Earth’s interior.     

Obtaining the free frequencies of the non‐rigid Earth

In this paper, the mathematical algorithm elaborated by González, Getino and Farto (1998) is applied to four different nonrigid Earth models, in order to obtain the analytical expressions of the corresponding free frequencies. The solutions are studied, and the contributions of the different considered effects are evaluated. A numerical integration is also carried out, showing the validity of the obtained analytical solutions.     

Opposite polarity dusty plasma expansion in Earth’s mesosphere

The self-similar expansion method is applied on the fluid system of equations which describes a plasma system consisting of opposite polarity dust grains, positive ions and electrons. The resultant system of equations is solved numerically to study the properties of the plasma expansion of this system. It is found that the presence of the second species of the dust has a great effect on the properties of the expansion of the other species.     

The Translational–Rotational Motion of an Earth Spheroid Satellite

In this paper we consider the translational–rotational motion of a spheroid satellite in the gravitational field, taking into account the asphericity of the earth. The harmonic coefficients of the earth’s gravitational field are taken up to J ₄. The equations of motion are obtained in terms of the canonical elements of Delaunay-Andoyer. A first order solution is obtained using the perturbing technique of Lie series.     

Dynamics of the Earth’s core from VLBI observations

The Earth’s rotation is accompanied by free circadian oscillations of its liquid core in the inner cavity of the lower mantle, which perturb the angular momentum of the entire Earth and produce an additional free nutation of the celestial pole called free core nutation (FCN). Since this nutation causes resonances in the diurnal tides and in the expansions of luni—solar nutation, its study, especially an improvement of the FCN period, is of fundamental importance for the theory of the Earth’s rotation. We have determined the FCN parameters from a joint analysis of equidistant series of coordinates of the celestial pole obtained from the combined processing of VLBI observations on global networks of stations for the interval 1984.0–2008.4 by IERS (International Earth Rotation and Reference System Service, Paris, France) and NEOS (National Earth Orientation Service, Washington, USA). Applying a moving least-squares filter (MLSF) to these data has shown that the FCN period averaged over this time interval differs significantly from the theoretical one and its phase varies over a wide range. Using the mean quadratic collocation (MQC) method, we have obtained a new, more accurate stochastic FCN model. Its analysis by the envelope method has revealed long-term linear phase trends, calling into question not only the adopted FCN period but also its stability and, hence, the stability of the resonant effects in the Earth’s luni—solar nutation.     

Statistical Characteristics of Meteoroid Parameters in the Earth’s Atmosphere

Statistical characteristics of meteoroids with kinetic energy from 0.1 to 440 kt TNT are estimated based on NASA satellite observations made in 1994–2016. The distributions of the number of falling meteoroids are constructed and analyzed based on the values of their initial kinetic energy, initial velocity, initial mass, altitude, geographic coordinates of the maximum total radiated energy region, and the year of the fall. Correlation dependences “mass–initial kinetic energy,” “maximum total radiated energy region altitude–initial kinetic energy,” and “maximum total radiated energy region altitude–initial velocity (the square of the initial velocity)” are constructed.