Secular Dynamics of Asteroids in the Inner Solar System

In the last three years we have carried out numerical and semi-analytical studies on the secular dynamical mechanisms in the region (semimajor axis a < 2 AU) where the NEA orbits evolve. Our numerical integrations (over a time span of a few Myr) have shown that: (i) the linear secular resonances with both the inner and the outer planets may play an important role in the dynamical evolution of NEAs; (ii) the apsidal secular resonance with Mars could provide an important dynamical transport mechanism by which asteroids in the Mars-crossing region eventually achieve Earth-crossing orbits; (iii) in this region, due to the interaction with the terrestrial planets, the Kozai resonance can occur at small inclinations, with the argument of perihelion ω librating around 0° or 180°, providing a temporary protection mechanism against close approaches to the planets. The location of the linear secular resonances in this zone has also been obtained by an automatic procedure using a semi-numerical method valid for all values of the inclinations and eccentricities of the small bodies, and also in the case of libration of the argument of perihelion. A map of the secular resonances in the (a, i) plane shows — in agreement with the numerical integrations — that all the resonances with the terrestrial and giant planets are present, and also that some of them overlap. Thus the way is now open to fully take into account secular resonances in modelling the dynamical evolution of NEAs. This revised version was published online in July 2006 with corrections to the Cover Date.     

Follow-up Ephemerides and the Accuracy of Preliminary Orbits

The accuracy of follow-up ephemerides and orbits of asteroids based on a short arc is investigated. The geocentric distance and its time derivative are the essential parameters determining the accuracy of predicted positions. The consequences for the optimal observing strategies are discussed. The node and inclination of the orbital plane are restricted enough to be useful for a preliminary sorting of candidates for linking.     

Dynamical Behavior of Asteroids Near Resonance: the 4:1 gap and the 7:2 Group

A comparative study of the evolution of the Sun–Jupiter–Asteroid system near the 4:1 and 7:2 resonances is performed by means of two techniques that proceed differently from the Hamiltonian corresponding to the planar restricted elliptic three-body problem. One technique is based on the classical Schubart averaging while the other is based on a mapping method in which the perturbing part of the Hamiltonian is expanded and the resulting terms are ordered according to a weight function that depends on the powers of eccentricities and the coefficients of the terms. For the mapping method the effect of Saturn on the asteroidal evolution is introduced and the degree of chaos is estimated by means of the Lyapunov time. Both methods are shown to lead to similar results and can be considered a suitable tool for describing the evolution of asteroids in the Kirkwood gap and the group corresponding to the 4:1 and 7:2 Jovian resonances, respectively. This revised version was published online in July 2006 with corrections to the Cover Date.     

Expansion of the disturbing function for high eccentricity and large amplitude of libration

Following some ideas, developed by Woltjer (1928), Message (1989), Yokoyama (1988, 1989) and Duriez (1990) an expansion of the disturbing function is given for high values of the eccentricity and large amplitude of libration. The classical expansion can be obtained as a particular case of the present model. Several asteroids with high eccentricity and large amplitude of libration are tested and the results are much better than those obtained from the classical theory.     

Resonant Periodic Orbits of Trans-Neptunian Objects

We study two and three-dimensional resonant periodic orbits, usingthe model of the restricted three-body problem with the Sun andNeptune as primaries. The position and the stability character ofthe periodic orbits determine the structure of the phase space andthis will provide useful information on the stability and longterm evolution of trans-Neptunian objects. The circular planarmodel is used as the starting point. Families of periodic orbitsare computed at the exterior resonances 1/2, 2/3 and 3/4 withNeptune and these are used as a guide to select the energy levelsfor the computation of the Poincaré maps, so that all basicresonances are included in the study. Using the circular planarmodel as the basic model, we extend our study to more realisticmodels by considering an elliptic orbit of Neptune and introducingthe inclination of the orbit. Families of symmetric periodicorbits of the planar elliptic restricted three-body problem andthe three-dimensional problem are found. All these orbitsbifurcate from the families of periodic orbits of the planarcircular problem. The stability of all orbits is studied. Althoughthe resonant structure in the circular problem is similar for allresonances, the situation changes if the eccentricity of Neptuneor the inclination of the orbit is taken into account. All theseresults are combined to explain why in some resonances there aremany bodies and other resonances are empty.     

An Explicit Symplectic Integrator for Cometary Orbits

A new symplectic algorithm is developed for cometary orbit integrations. The integrator can handle both high-eccentricity orbits and close encounters with planets. The method is based on time transformations for Hamiltonians separated into Keplerian and perturbation parts. The adaptive time-step of this algorithm depends on the distance from a centre and the magnitude of perturbations. The explicit leapfrog technique is simple and efficient.     

The Prograde C7 Resonance for Earth and Mars Satellite Orbits

The resonance C7 is a 1:1 eccentricity (apsidal) resonance between the longitude of a satellite's pericentre and the mean longitude of the Sun. A previous paper by the author (Breiter, 1999) identified it as the strongest of the lunisolar apsidal resonances. After the reduction to a single degree of freedom, the problem is studied qualitatively for the prograde orbits around the Earth and Mars. Pitchfork, saddle-node, and saddle connection bifurcations give rise to a complicated phase flow, which may involve up to nine critical points. This revised version was published online in July 2006 with corrections to the Cover Date.     

A stability study of asteroid families near the 3:1 and 5:2 resonance with Jupiter

By using theD-criterion Lindblad (1992) has identified 14 asteroid families from a sample of 4100 numbered asteroids with proper elements from Milani and Kneevi (1990). Taxonomic types and other physical properties for a significant number of objects in five of the families show strong homogeneity within each family, further strengthening their internal relationship.To test the hypothesis of a common origin in, e.g., a catastrophic collision event, we have set out to integrate the orbits of the members of the Maria, Dora and Oppavia-Gefion families over some 10⁶ years. The mean distance for the Maria family is close to the 3:1 mean-motion resonance with Jupiter, while the other two families lie close to the 5:2 resonance.We used a simplified solar system model which included the perturbations by Jupiter and Saturn only and implemented Everhart's variable stepsize integrator RA15. All close encounters between the family members (within 0.1 AU) were recorded as well. Preliminary results from integrations over 4×10⁵ years are presented here.The statistics of close encounters show pronounced peaks for several members within each family, while for others no significant levels above the background of random encounters or even very low frequencies were found. This indicates a subclustering within the families. Quite a lot of very close (<0.005 AU) mutual encounters are found, which suggest that, at least for the larger members in a family, the mutual gravitational interactions could be of some importance for the real orbital evolutions.The encounter statistics between the Dora and Oppavia family members suggest a possible interrelationship between this two groups.     

Generalized Averaging Principle and the Secular Evolution of Planet Crossing Orbits

Planet crossing orbits give rise to mathematical singularities that make it not possible to apply the classical averaging principle to study the qualitative evolution of Near Earth Asteroids (NEAs). Recently this principle has been generalized to deal with crossings in a mathematical model with the planets on circular coplanar orbits. More accuracy is needed to compute the averaged evolution of planet crossing orbits for different purposes: computing reliable crossing times for the averaged motion, writing more precise proper elements and frequencies for NEAs, etc. In this paper we present the generalization of the averaging principle using a model where the eccentricity and the inclination of the planets are taken into account.     

A symplectic mapping approach to the study of the stochasticity in asteroidal resonances

In the case of the 2:1 and 3:2 resonances with Jupiter, it has not been yet possible to have a complete identification of all chaotic diffusion processes at work, mainly because the time scale of some of them are of an order still out of the reach of precise integrations. A planar Hadjidemetriou's mapping, using expansions valid for high eccentricities and scaled in order to accelerate the diffusion processes, was derived. The solutions obtained with the mapping show huge eccentricity variations in all orbits starting in the middle of the 2:I resonance, when the main short-period perturbations of Jupiter's orbit are considered. The solutions starting in the middle of the 3:2 resonance do not show any important diffusion.     

Angular momentum of binary asteroids: Implications for their possible origin

We describe in this work a thorough study of the physical and orbital characteristics of extensively observed main-belt and trojan binaries, mainly taken from the LAOSA (Large Adaptive Optics Survey of Asteroids [Marchis, F., Baek, M., Berthier, J., Descamps, P., Hestroffer, D., Kaasalainen, M., Vachier, F., 2006c. In: Workshop on Spacecraft Reconnaissance of Asteroid and Comet Interiors. Abstract #3042]) database, along with a selection of bifurcated objects. Dimensionless quantities, such as the specific angular momentum and the scaled primary spin rate, are computed and discussed for each system. They suggest that these asteroidal systems might be the outcome of rotational fission or mass shedding of a parent body presumably subjected to an external torque. One of the most striking features of separated binaries composed of a large primary (Rp>100 km) with a much smaller secondary (Rs<20 km) is that they all have total angular momentum of ∼0.27. This value is quite close to the Maclaurin-Jacobi bifurcation (0.308) of a spinning fluid body. Alternatively, contact binaries and tidally locked double asteroids, made of components of similar size, have an angular momentum larger than 0.48. They compare successfully with the fission equilibrium sequence of a rotating fluid mass. In conclusion, we find that total angular momentum is a useful proxy to assess the internal structure of such systems.     

Novel Analysis of Tadpole and Horseshoe Orbits

A simple analysis of tadpole and horseshoe orbits is performed using a new perturbation technique. The first-order results are presented explicitely, extension to any perturbation order is delineated. This revised version was published online in July 2006 with corrections to the Cover Date.     

Earth impact probability of the Asteroid (25143) Itokawa to be sampled by the spacecraft Hayabusa

The Japanese spacecraft Hayabusa is planed to reach the Asteroid Itokawa in September 2005, and to bring back some samples of its surface to Earth in 2007. We have studied the future possible evolution of this asteroid by integrating numerically over 100 Myr a set of 39 initially indistinguishable orbits (clones), obtained either by small variations of the nominal initial conditions, or by using different computers (introducing different round-off errors). The results indicate that an Earth impact of this 500-m-size asteroid is likely within a million years, which is only a factor of four larger than the average impact frequency of asteroids of this size. The mission Hayabusa may thus sample a good candidate for being among the next 500-m-size Earth impactors.     

An Algebraic Method to Compute the Critical Points of the Distance Function Between Two Keplerian Orbits

We describe an efficient algorithm to compute all the critical points of the distance function between two Keplerian orbits (either bounded or unbounded) with a common focus. The critical values of this function are important for different purposes, for example to evaluate the risk of collisions of asteroids or comets with the Solar system planets. Our algorithm is based on the algebraic elimination theory: through the computation of the resultant of two bivariate polynomials, we find a 16th degree univariate polynomial whose real roots give us one component of the critical points. We discuss also some degenerate cases and show several examples, involving the orbits of the known asteroids and comets.      

Short-lived asteroids in the 7/3 Kirkwood gap and their relationship to the Koronis and Eos families

A population of 23 asteroids is currently observed in a very unstable region of the main belt, the 7/3 Kirkwood gap. The small size of these bodies—with the notable exception of (677) Aaltje (∼30 km)—as well as the computation of their dynamical lifetimes (3<T_D<172 Myr) shows that they cannot be on their primordial orbits, but were recently injected in the resonance. The distribution of inclinations appears to be bimodal, the two peaks being close to 2° and 10°. We argue that the resonant population is constantly being replenished by the slow leakage of asteroids from both the Koronis (I∼2°) and Eos (I∼10°) families, due to the drift of their semi-major axes, caused by the Yarkovsky effect. Assuming previously reported values for the Yarkovsky mean drift rate, we calculate the flux of family members needed to sustain the currently observed population in steady state. The number densities with respect to semi-major axis of the observed members of both families are in very good agreement with our calculations. The fact that (677) Aaltje is currently observed in the resonance is most likely an exceptional event. This asteroid should not be genetically related to any of the above families. Its size and the eccentricity of its orbit suggest that the Yarkovsky effect should have been less efficient in transporting this body to the resonance than close encounters with Ceres.     

Resonant dynamics of Medium Earth Orbits: space debris issues

The Medium Earth Orbit region is the home of the navigation constellations. It is shown how the orbits of these constellations of satellites are strongly affected by the so-called inclination dependent luni-solar resonances. The analytical theory of these resonances is recalled and a large set of numerical integrations is used to investigate the stability of the orbits of the constellations over very long time spans. The stability issue is important in the definition of possible disposal strategies for the constellation spacecraft, after their end-of-life. Two possible disposal strategies are envisaged involving either stable or unstable orbits (from the eccentricity growth point of view). In particular it is shown how, to have disposal orbits with moderately short lifetime (of the orders of 40–50 years), a very large disposal maneuver would be required to rise the inital eccentricity to values above ~ 0.3.     

Micrometeorites and Their Implications for Meteors

Micrometeorites (MMs) are extraterrestrial dust particles, in the size range 25–400 μm, recovered from the Earth’s surface. They have experienced a wide range of heating during atmospheric entry from completely molten spherules to particles heated to temperatures <300°C that have retained low temperature minerals. The majority of MMs have mineralogies, textures and compositions that strongly resemble components from chondritic meteorites suggesting these correspond to sporadic, low geocentric velocity meteors. Changes in MMs due to entry heating, however, have implications for meteoric processes in general that may allow the observed behaviour of meteors to be directly related to the material properties of their meteoroids.     

Numerical Evaluation of the General Yarkovsky Effect: Effects on Eccentricity and Longitude of Periapse

J. N. Spitale and R. Greenberg (2001, Icarus149, 222-234) developed a nonlinearized, finite-difference solution to the heat equation that yields orbital rates of change due to the Yarkovsky effect for small, spherical, bare-rock asteroids and used it to investigate changes in semimajor axis caused by the Yarkovsky effect. Here, we present results for changes in eccentricity and longitude of periapse. These results may be useful as benchmarks for simplified analytical solutions. Moreover, we explore a range of parameters, some of which are inaccessible to most other approaches. Instantaneous rates can be quite fast: For a 1-m scale body rotating with a 5-h period, de/dt can be as fast as 0.1 per million years (da/dt rates for similar test bodies were reported in J. N. Spitale and R. Greenberg (2001, Icarus149, 222-234)). For more typical rotation periods, these rates would be considerably slower. Output from our calculation method could be used in simulations of asteroid population evolution such as that by W. F. Bottke, D. P. Rubincam, and J. A. Burns (2000, Icarus145, 301-331). On long time scales, impacts would randomize the spin axis before significant orbital evolution could occur. Nevertheless, occasional favorable rotation states might persist long enough for substantial eccentricity changes to accumulate (1) if the body is decoupled from the main belt (e.g., many near-Earth asteroids), (2) if the population of very small (mm-scale) main-belt impactors is less than expected, or (3) if our numerical results are scaled up to km-size bodies.