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.     

Orbital evolution studies of planet-crossing asteroids

Numerical integrations of 26 orbits of planet-crossing astetoids of Apollo-Amor type have been performed, in a solar system model including the perturbations by the planets from Venus to Neptune. The 15:th order RADAU integrator (Everhart, 1985) has been used. Orbits for the asteroids 433 Eros, 887 Alinda, 1036 Ganymed, 1221 Amor, 1580 Betulia, 1627 Ivar, 1685 Toro, 1862 Apollo, 1863 Antinous, 1864 Daedalus, 1865 Cerberus, 1915 Quetzalcoatl and 1916 Boreas have been integrated over 100 000 years forward in time and for 1866 Sisyphus, 2102 Tantalus, 2201 Oljato, 2329 Orthos, 3360 1981 VA, 3552 1983 SA, 1981 EJ30, 1985 PA, 1985 WA, 1986 DA 1986 JK and 1986 RA a period of about 33 000 years has been covered. The orbital evolutions of these asteroids are discussed. This work is part of a larger study of the long-term orbital evolution of planet-crossing asteroids and will be continued within the project SPACEGUARD (Milani et al., 1987).     

Dynamics of Earth‐crossing asteroids: the protected Toro orbits

Some asteroids in Earth‐crossing orbits avoid close approaches by entering in a mean motion resonance whenever the distance between the two orbits is small. These orbits are ‘Toro class’ according to the classification of (Milani et al., 1989). This protection mechanism can be understood by a semi‐averaged model, in which the fast variables are removed and the dynamical variables are the critical argument and the semimajor axis, with dependence upon a slow parameter. The adiabatic invariant theory can be applied to this model and accounts for all the qualitative features of the orbits in this class, including the onset of the libration when the orbit distance is small. Because of the neglected perturbations by the other planets, this theory is approximate and the adiabatic invariant is conserved only with low accuracy moreover, the Toro state can be terminated by a close approach to another planet (typically Venus). “Would you tell me, please, which way I ought to go from here?” “That depends a good deal on where you want to get to,” said the Cat. Alice in Wonderland, L. Carroll     

A Simple Method for Numerical Calculations of the Evolution of Orbits of Near-Earth Asteroids

A simple method for numerical integration of the equations of motion of small bodies of the Solar System is proposed, which is especially efficient in studying the orbits with small perihelion distances. The evolution of orbits of 121 numbered asteroids with perihelion distances q < 1.2 AU is investigated over the time interval of years 2000–2100 with allowance made for the gravitational influence of nine planets and three largest asteroids. The circumstances of close encounters of asteroids with the Earth and other terrestrial planets are presented.     

Model orbits of asteroid 3040 Kozai

The orbital evolutions of the asteroid 3040 Kozai and model asteroids with similar orbits have been investigated. Their osculating orbits for an epoch 1991 December 10 were numerically integrated forward within the interval of 20,000 years, using a dynamical model of the solar system consisting of all inner planets, Jupiter, and Saturn.The orbit of the asteroid Kozai is stable. Its motion is affected only by long-period perturbations of planets. With change of the argument of perihelion of the asteroid Kozai, the evolution of the model asteroid orbits changes essentially, too. The model orbits with the argument of perihelion changed by the order of 10% show that asteroids with such orbital parameters may approach the Earth orbit, while asteroids with larger changes may even cross it, at least after 10,000 years. Long-term orbital evolution of asteroids with these orbital parameters is very sensitive on their angular elements.     

The large-scale structure of the asteroid belt

We examine the distributions of 2888 numbered minor planets over orbital inclination, eccentricity, and semimajor axis, and define 19 zones which we believe adequately to isolate the selection biases in survey programs of the physical properties of minor planets. Six numbered asteroids have exceptional orbits and fall into no zone. We also call attention to rather sharp upper limits, which become increasingly stringent at larger heliocentric distances, on orbital inclinations and eccentricity.     

The Trojan problem

There has been a renewed interest in the Trojan problem in recent years. Significant progress has been made in discovering and understanding dynamical features of motion of Jupiter's Trojan asteroids. The dynamics of hypothetical Trojan-type asteroids of other major planets has also been the subject of several recent investigations. This paper offers an overview on the current status of researches on real and hypothetical Trojan asteroids of the major planets. Results of analytical and numerical works are surveyed. Questions of dynamical properties, long-term evolution of orbits, stability regions around the triangular Lagrangian points are discussed among other problems of the Trojans.     

The asteroid ephemerides program at the U.S. Naval Observatory

The U.S. Naval Observatory has begun a program of ephemeris improvement and reference frame determination from the main belt asteroids. The program is, currently, starting out with a limited set of observations of the larger asteroids to determine the equator and equinox corrections for the USNO W1J00 transit circle observations catalog, and, if possible, improve the orbits of these asteroids based on this limited set of observations. For this project, transit circle observations of the Sun and the planets Mercury through Jupiter, are also being used to determine the equator, equinox, and ephemeris corrections, the next goal is to improve the orbits of the larger asteroids in the optical reference frame using observation series that cover a much longer period of time. This will allow the exploration of the differences between the dynamical reference frame based on radar observations of main belt asteroids and its relation with the optical reference frame. Other goals include the exploration of the mass distribution in the main asteroid belt from high precision observations, and the effect of this mass on the ephemerides of the major planets.     

Numerical Simulation of the Orbital Evolution of Near-Earth Asteroids Close to Mean Motion Resonances

The dynamics of near-Earth asteroids near mean motion resonances with the Earth or other planets is considered. The probability domains of the motion of some near-Earth asteroids close to low-order resonances are presented. The investigations have been carried out by means of a numerical integration of differential equations, taking into account the perturbations from the major planets and the Moon. For each investigated object an ensemble of 100 test particles with orbital elements nearby those of the nominal orbit has been constructed and its evolution has been retraced over the time interval (–3000, +3000 years). The initial set of orbits has been generated on the basis of probable variations of the initial orbital elements obtained from the least square analysis of observations.This revised version was published online in October 2005 with corrections to the Cover Date.     

Planetary close encounters: geometry of approach and post-encounter orbital parameters

Öpik's assumptions on the geometry of particle trajectories leading to and through planetary close encounters are used to compute the distribution of changes in heliocentric orbital elements that result from such encounters for a range of initial heliocentric orbits. Behaviour at encounter is assumed to follow two-body (particle—planet) gravitational scattering, while before and after encounter particle motion is only governed by the force of the Sun. Derivation of these distributions allows precise analysis of the probability of various outcomes in terms of the physical characteristics of the bodies involved. For example, they allow an explanation and prediction of the asymmetry of the extreme energy perturbations for different initial orbits. The formulae derived here may be applied to problems including the original accumulation of planets and satellites, and the continuing evolution of populations of small bodies, such as asteroids and comets.     

Surface Pyroxenes and Olivines of S-Type Asteroids as Studied by Modeling Their Reflectance Spectra

We model the reflectance spectra of SI–SVII-subtype asteroids. The spectra of minor planets contain little information as regards the abundance of free metal and the form of its existence on the asteroid surfaces. We investigate the properties of a set of probable spectra for the surface pyroxene and olivine of minor planets. Clinopyroxene with high calcium and iron contents and orthopyroxene with a ferrosilite molecular content of about 40% are typical of the S-type asteroids. The subtype number of S-type asteroids can be correlated with the bulk pyroxene composition. The forsterite molecular content in asteroid olivine lies within the probable range 40–74%. The prevalence of pyroxene over olivine and feldspar in the computed compositions of S-type asteroids suggests the probable presence of basalts on their surfaces.     

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.     

Asteroidal jet streams

Detailed studies of the asteroidal belt are of importance for clarifying whether the asteroids are fragments of an exploded planet or represent an intermediate state in the accretion of planets.A study of the Hirayama family Flora shows that it contains three groups of bodies travelling in almost identical orbits, thus constituting three jet streams.It is shown that the formation of jet streams is difficult to reconcile with the exploded planet view.In order to decide whether the formation of jet streams is a general phenomenon in the asteroidal belt, it is necessary to determine the orbits of a large number of small asteroids.     

Beobachtungen Kleiner Planeten. XIII

Es wird eine Übersicht gegeben über die im Jahr 1992 auf Tautenburger Schmidtplatten gefundenen Kleinen Planeten. Es wurden 518 Objekte beobachtet und für diese 2396 Positionen gerechnet. Es handelt sich um 408 Planetoiden mit provisorischer Bezeichnung, unter denen 274 neu vergebene Bezeichnungen sind, sowie um 110 numerierte Objekte. Für 172 Tautenburger Objekte wurden Bahnen aus einer Opposition gerechnet. Für 55 im Berichtszeitraum numerierte Planeten haben Tautenburger Positionen einen Beitrag geleistet. Acht Tautenburger Planetoiden wurden numeriert (Tab. 3). A summary is given about the Minor Planet survey performed in 1992 on Tautenburg Schmidt plates. 518 asteroids were observed and 2396 positions are calculated for them. These are 408 asteroids with provisional designations (274 of them have new designations) and 110 numbered objects. One-opposition orbits have been computed for 172 Tautenburg asteroids. Tautenburg observations could give a contribution to 55 planets numbered in the period of this report. Eight Tautenburg asteroids have been numbered (Tab. 3).     

The Effect of Tidal Interaction with a Gas Disk on Formation of Terrestrial Planets

We have performed N-body simulation on final accretion stage of terrestrial planets, including the effect of damping of eccentricity and inclination caused by tidal interaction with a remnant gas disk. As a result of runway and oligarchic accretion, about 20 Mars-sized protoplanets would be formed in nearly circular orbits with orbital separation of several to ten Hill radius. The orbits of the protoplanets would be eventually destabilized by long-term mutual gravity and/or secular resonance of giant gaseous planets. The protoplanets would coalesce with each other to form terrestrial planets through the orbital crossing. Previous N-body simulations, however, showed that the final eccentricities of planets are around 0.1, which are about 10 times higher than the present eccentricities of Earth and Venus. The obtained high eccentricities are the remnant of orbital crossing. We included the effect of eccentricity damping caused by gravitational interaction with disk gas as a drag force (“gravitational drag”) and carried out N-body simulation of accretion of protoplanets. We start with 15 protoplanets with 0.2M⊕ and integrate the orbits for 10⁷ years, which is consistent with the observationally inferred disk lifetime (in some runs, we start with 30 protoplanets with 0.1M⊕). In most runs, the damping time scale, which is equivalent to the strength of the drag force, is kept constant throughout each run in order to clarify the effects of the damping. We found that the planets' final mass, spatial distribution, and eccentricities depend on the damping time scale. If the damping time scale for a 0.2M⊕ mass planet at 1 AU is longer than 10⁸ years, planets grow to Earth's size, but the final eccentricities are too high as in gas-free cases. If it is shorter than 10⁶ years, the eccentricities of the protoplanets cannot be pumped up, resulting in not enough orbital crossing to make Earth-sized planets. Small planets with low eccentricities are formed with small orbital separation. On the other hand, if it is between 10⁶ and 10⁸ years, which may correspond to a mostly depleted disk (0.01-0.1% of surface density of the minimum mass model), some protoplanets can grow to about the size of Earth and Venus, and the eccentricities of such surviving planets can be diminished within the disk lifetime. Furthermore, in innermost and outermost regions in the same system, we often find planets with smaller size and larger eccentricities too, which could be analogous to Mars and Mercury. This is partly because the gravitational drag is less effective for smaller mass planets, and partly due to the “edge effect,” which means the innermost and outermost planets tend to remain without collision. We also carried out several runs with time-dependent drag force according to depletion of a gas disk. In these runs, we used exponential decay model with e-folding time of 3×10⁶ years. The orbits of protoplanets are stablized by the eccentricity damping in the early time. When disk surface density decays to ?1% of the minimum mass disk model, the damping force is no longer strong enough to inhibit the increase of the eccentricity by distant perturbations among protoplanets so that the orbital crossing starts. In this disk decay model, a gas disk with 10⁻⁴-10⁻³ times the minimum mass model still remains after the orbital crossing and accretional events, which is enough to damp the eccentricities of the Earth-sized planets to the order of 0.01. Using these results, we discuss a possible scenario for the last stage of terrestrial planet formation.     

From Astrometry to Celestial Mechanics: Orbit Determination with Very Short Arcs

Contemporary surveys provide a huge number of detections of small solar system bodies, mostly asteroids. Typically, the reported astrometry is not enough to compute an orbit and/or perform an identification with an already discovered object. The classical methods for preliminary orbit determination fail in such cases: a new approach is necessary. When the observations are not enough to compute an orbit we represent the data with an attributable (two angles and their time derivatives). The undetermined variables range and range rate span an admissible region of solar system orbits, which can be sampled by a set of Virtual Asteroids (VAs) selected by an optimal triangulation. The attributable results from a fit and has an uncertainty represented by a covariance matrix, thus the predictions of future observations can be described by a quasi-product structure (admissible region times confidence ellipsoid), which can be approximated by a triangulation with each node surrounded by a confidence ellipsoid. The problem of identifying two independent short arcs of observations has been solved. For each VA in the admissible region of the first arc we consider prediction at the time of the second arc and the corresponding covariance matrix, and we compare them with the attributable of the second arc with its own covariance. By using the penalty (increase in the sum of squares, as in the algorithms for identification) we select the VAs which can fit together both arcs and compute a preliminary orbit. Even two attributables may not be enough to compute an orbit with a convergent differential corrections algorithm. The preliminary orbits are used as first guess for constrained differential corrections, providing solutions along the Line Of Variations (LOV) which can be used as second generation VAs to further predict the observations at the time of a third arc. In general the identification with a third arc will ensure a least squares orbit, with uncertainty described by the covariance matrix.     

Massive identification of asteroids in three-body resonances

An essential role in the asteroidal dynamics is played by the mean motion resonances. Two-body planet–asteroid resonances are widely known, due to the Kirkwood gaps. Besides, so-called three-body mean motion resonances exist, in which an asteroid and two planets participate. Identification of asteroids in three-body (namely, Jupiter–Saturn–asteroid) resonances was initially accomplished by Nesvorný and Morbidelli (Nesvorný D., Morbidelli, A. [1998]. Astron. J. 116, 3029–3037), who, by means of visual analysis of the time behaviour of resonant arguments, found 255 asteroids to reside in such resonances. We develop specialized algorithms and software for massive automatic identification of asteroids in the three-body, as well as two-body, resonances of arbitrary order, by means of automatic analysis of the time behaviour of resonant arguments. In the computation of orbits, all essential perturbations are taken into account. We integrate the asteroidal orbits on the time interval of 100,000 yr and identify main-belt asteroids in the three-body Jupiter–Saturn–asteroid resonances up to the 6th order inclusive, and in the two-body Jupiter–asteroid resonances up to the 9th order inclusive, in the set of ～250,000 objects from the “Asteroids – Dynamic Site” (AstDyS) database. The percentages of resonant objects, including extrapolations for higher-order resonances, are determined. In particular, the observed fraction of pure-resonant asteroids (those exhibiting resonant libration on the whole interval of integration) in the three-body resonances up to the 6th order inclusive is ≈0.9% of the whole set; and, using a higher-order extrapolation, the actual total fraction of pure-resonant asteroids in the three-body resonances of all orders is estimated as ≈1.1% of the whole set.     

Impact solutions of Asteroid 2007 WD5 with Mars

A method for computing impact probabilities between asteroids and the planet Mars is presented that uses impact clones and validation analysis based on a normal distribution of computed errors. This method uses OrbFit software, and we present a calculation of the impact probabilities between Asteroid 2007 WD₅ and Mars, which passed within about 20,000 km of the martian surface on January 30, 2008. This method can be generalized for computing impact probabilities between asteroids and other planets including Earth. Presented method applies in principal the same technique already in use for years at the JPL NASA and by the group of researchers at the University of Pisa [Milani, A., Chesley, S.R., Sansaturio, M.E., Tommei, G., Valsecchi, G.B., 2005a. Icarus 173, 362].     

Identification of Asteroids in Two-Body Resonances

A catalog of asteroids in two-body orbital resonances with the planets of the Solar System has been created. The AstDyS database was a source of the input data, and all the numbered objects (467303 objects at the time of the analysis) were considered. The orbits were integrated in the framework of a pure gravitational problem considering all the planets of the Solar System and Pluto. From the analysis of the behavior of the resonant argument and the semimajor axis on the 100-kyr interval, the resonance membership and the libration type (pure or transient) were verified for each of the asteroids. A more accurate method to identify the resonant argument librations was developed on the basis of the analysis of mutual periodograms. We found 23251 resonant asteroids, ~36% of which (8397 objects) are in pure resonances.     

Quasi-satellite orbits in the general context of dynamics in the 1:1 mean motion resonance: perturbative treatment

Our investigation is motivated by the recent discovery of asteroids orbiting the Sun and simultaneously staying near one of the Solar System planets for a long time. This regime of motion is usually called the quasi-satellite regime, since even at the times of the closest approaches the distance between the asteroid and the planet is significantly larger than the region of space (the Hill’s sphere) in which the planet can hold its satellites. We explore the properties of the quasi-satellite regimes in the context of the spatial restricted circular three-body problem “Sun–planet–asteroid”. Via double numerical averaging, we construct evolutionary equations which describe the long-term behaviour of the orbital elements of an asteroid. Special attention is paid to possible transitions between the motion in a quasi-satellite orbit and the one in another type of orbits available in the 1:1 resonance. A rough classification of the corresponding evolutionary paths is given for an asteroid’s motion with a sufficiently small eccentricity and inclination.