How precise is the orbit of asteroid (99942) Apophis and how probable is its collision with the Earth in 2036–2037?

The nearest in time close approach of potentially hazardous asteroid (99942) Apophis with the Earth will take place on April 13, 2029, when the minimum distance of the asteroid from the Earth’s center will be as small as 38 000 km. Such a close approach will result in substantial transformation of the asteroid’s orbit. The value of the perturbations depends on the minimum distance between the bodies during the approach. Among possible transformations of the orbit are those which result in new dangerous approaches and even in probable Apophis collisions with the Earth starting from 2036. At present, at least four solutions are known for the Apophis orbit which were obtained using all radar and most of available optical observations. The procedures of assigning weights to conditional equations and the models of the asteroid’s motion have differed to some extent when finding these solutions. Of considerable interest is the comparison of the found orbital parameters with the estimates of their accuracy, since small distinctions in their values result in considerable distinctions in the forecast of Apophis’ motion after 2029 and beyond. It is shown in the paper that the estimates of the probability of an Apophis collision with the Earth in 2036 differ by some orders of magnitude, according to various solutions. The influence of factors which were disregarded in the models of motion even more increases the uncertainty in forecasting the motion after 2029. More accurate forecasting can be achieved as a result of additional optical and, to a greater extent, a series of radar observations in 2013 and then in 2020–2021, and/or as a result of processing radio signals of the transmitter delivered to the Apophis surface or to the orbit of its artificial satellite, as it was proposed in a number of papers.     

Impact trajectories of the asteroid Apophis in the 21st century

The asteroid Apophis is one of the most hazardous near-Earth asteroids. As a result of the scattering of Apophis?? potential trajectories after its close approach in 2029, and its possible approach in 2036, there are many dangerous trajectories including impact trajectories after 2036. The purpose of this study is to identify and investigate these trajectories. We use the Everhart integrator; the DE405, DE423, and EPM2008 ephemerides; and two sets of initial data for Apophis (those collected by NASA in 2006 and by the IAA in 2010). More than 50 possible encounters in this century are presented, including 13 encounters between 2036 and 2050. The minimum geocentric distances obtained using a different ephemeris and initial conditions differ little between themselves. Analogous results in (Yeomans et al., 2009) are consistent with our results.     

On the Possibility of Deflecting an Asteroid from Collision with the Earth Using the Kinetic Method

The possibilities of deflecting an asteroid from its collision course with the Earth by changing its velocity with an impact are considered. Using the asteroid Apophis as an example, the time dependence of the positions and sizes of the keyholes leading to collision is studied. It has been found that the possibility of deflecting this asteroid usually exists, and the impact can be accomplished in principle, given the capabilities of modern space technology. A change in the velocity should be performed before the encounter of 2029 in order to use the gravitational maneuver effect. The possible accuracy of determining Apophis’ orbit and the keyholes that lead to collision and are associated with the resonance returns are considered.     

How selection and weighting of astrometric observations influence the impact probability. The case of asteroid (99942) Apophis

The aim of this paper is to show that in the case of a low probability of asteroid collision with the Earth, the appropriate selection and weighting of the data are crucial for the impact investigation and for analysing the impact possibilities using extensive numerical simulations. By means of the Monte Carlo special method, a large number of 'clone' orbits have been generated. A full range of orbital elements in the six-dimensional parameter space, that is, in the entire confidence region allowed by the observational material, has been examined. On the basis of 1000 astrometric observations of (99942) Apophis, the best solutions for the geocentric encounter distance of  6.065 ± 0.081 R_⊕  (without perturbations by asteroids) or  6.064 ± 0.095 R_⊕  (including perturbations by the four largest asteroids) were derived for the close encounter with the Earth on 2029 April 13. The present uncertainties allow for special configurations ('keyholes') during this encounter that may lead to very close encounters in future approaches of Apophis. Two groups of keyholes are connected with the close encounter with the Earth in 2036 (within the minimal distance of  5.7736−5.7763 R_⊕  on 2029 April 13) and 2037 (within the minimal distance of  6.3359–6.3488 R_⊕  ). The nominal orbits for our most accurate models run almost exactly in the middle of these two impact keyhole groups. A very small keyhole for the impact in 2076 has been found between these groups at the minimal distance of 5.97347   R_⊕  . This keyhole is close to the nominal orbit. The present observations are not sufficiently accurate to eliminate definitely the possibility of impact with the Earth in 2036 and for many years after.     

On the prevention of a possible collision of asteroid Apophis with the Earth

Currently, there is some positive probability of a collision of the asteroid Apophis with the Earth in 2036. In this study, the problem of preventing the collision by correcting the asteroid’s orbit is examined. The characteristics of the impulsive correction are investigated, as well as the ways of its implementation by kinetic and nuclear impacts. Impulsive and weak effects are compared. Weak effects leading to slow changes in the asteroid’s orbit are considered to be more usable because of the potentially higher accuracy of this correction. The characteristics of the gravitational effect of the asteroid by a special spacecraft (SC) kept by its control jet engines at a certain point near the asteroid and gravitationally perturbing the motion of Apophis are analyzed. The change in the perigee radius of the Apophis orbit in 2036 and the SC mass consumption are examined as functions of the effect duration, the SC mass, its distance to the asteroid, the start time of the correction, and the velocity of the SC engine exhaust jet.     

Resonance orbits of near-Earth asteroids

Asteroid 99942 Apophis is one of the most hazardous NEAs (near-Earth asteroids) today. Some specific features of its travel are the possibility of repeated Earth approaches, loss of forecast precision due to trajectory dispersions, and nondeterministic motion. These specific features do not only characterize Apophis. Special methods are needed to find possible collision trajectories among these travels. These trajectories are located in the vicinity of resonance collision orbits.The present paper discusses methods of detecting hazardous trajectories in the event of nondeterministic motion and characterizing these trajectories as applied to asteroid Apophis, precision losses in the event of trajectory dispersions, conditions of determinacy losses, and hazardous trajectories in the vicinity of resonance orbits.     

On the impact of the Yarkovsky effect on Apophis’ orbit

In October 2009, a new set of optical observations of Apophis, a potentially hazardous asteroid, was published. These data have significantly expanded the interval of observations and their total number. In the article we compare the results of refinement of Apophis’ orbit made at the Jet Propulsion Laboratory (JPL, United States), the University of Pisa (Italy), and the Institute of Applied Astronomy (IAA) of the Russian Academy of Sciences with consideration for new observations. New orbits lead to a significant decrease in the probability of Apophis’ collision with the Earth in 2036. As a result of processing a large number of observations of asteroids approaching the Earth and main belt asteroids less than 40 km in size, with a large number of optical and, in many cases, radar observations in different oppositions, one of the authors revealed that additional acceleration affects their motion. This acceleration can be represented by the transversal component A ₂ in the orbital coordinate system. The presence of this acceleration can be interpreted as the Yarkovsky effect. The statistical properties of distribution of A ₂ for asteroids, for which it is determined quite reliably, evidence in favor of this interpretation. The value of additional acceleration for bodies the size of Apophis falls in the range ±10⁻¹³ AU/day². In this paper we have calculated the probability of Apophis colliding with the Earth in 2036 at different values of the transversal component of additional acceleration A ₂. For the resulting points, a plot of the probability of the collision against the A ₂ value has been constructed. At A ₂ = −8.748 × 10⁻¹⁴ AU/day² (and zero values of the radial A ₁ and normal A ₃ components) the nominal solution for Apophis’ orbit on April 13, 2029, is only 90 m from the middle of a “keyhole” 600 m in width, which leads to a collision of Apophis with the Earth in 2036. Since the scattering ellipse in the target plane in 2029 significantly overlaps the keyhole, the probability of collision at the given additional acceleration value is 0.0022. This result has been verified by the Monte Carlo method. Tests of 10000 random sets of orbital elements, which were found taking into account their correlation, have shown that 22 cases have resulted in virtual asteroids colliding with Earth in 2036. A plot of the probability of the collision against the value of A ₂ has been constructed.     

Orbital YORP and asteroid orbit evolution, with application to Apophis

Photon thrust from shape alone can produce quasi-secular changes in an asteroid's orbital elements. An asteroid in an elliptical orbit with a north–south shape asymmetry can steadily alter its elements over timescales longer than one orbital trip about the Sun. This thrust, called here orbital YORP (YORP = Yarkovsky–O'Keefe–Radzievskii–Paddack), operates even in the absence of thermal inertia, which the Yarkovsky effects require. However, unlike the Yarkovsky effects, which produce secular orbital changes over millions or billions of years, the change in an asteroid's orbital elements from orbital YORP operates only over the precession timescale of the orbit or of the asteroid's spin axis; this is generally only thousands or tens of thousands of years. Thus while the orbital YORP timescale is too short for an asteroid to secularly journey very far, it is long enough to warrant investigation with respect to 99942 Apophis, which might conceivably impact the Earth in 2036. A near-maximal orbital YORP effect is found by assuming Apophis is without thermal inertia and is shaped like a hemisphere, with its spin axis lying in the orbital plane. With these assumptions orbital YORP can change its along-track position by up to ±245 km, which is comparable to Yarkovsky effects. Though Apophis' shape, thermal properties, and spin axis orientation are currently unknown, the practical upper and lower limits are liable to be much less than the ±245 km extremes. Even so, the uncertainty in position is still likely to be much larger than the 0.5 km “keyhole” Apophis must pass through during its close approach in 2029 in order to strike the Earth in 2036.     

Spectral properties and composition of potentially hazardous Asteroid (99942) Apophis

The known close approach of Asteroid (99942) Apophis in April 2029 provides the opportunity for the case study of a potentially hazardous asteroid in advance of its encounter. The visible to near-infrared (0.55 to 2.45 μm) reflectance spectrum of Apophis is compared and modeled with respect to the spectral and mineralogical characteristics of likely meteorite analogs. Apophis is found to be an Sq-class asteroid that most closely resembles LL ordinary chondrite meteorites in terms of spectral characteristics and interpreted olivine and pyroxene abundances, although we cannot rule out some degree of partial melting. A meteorite analog allows some estimates and conjectures of Apophis' possible range of physical properties such as the grain density and micro-porosity of its constituent material. Composition and size similarities of Apophis with (25143) Itokawa suggest a total porosity of 40% as a “current best guess” for Apophis. Applying these parameters to Apophis yields a mass estimate of 2×¹⁰¹⁰ kg with a corresponding energy estimate of 375 Mt for its potential hazard. Substantial unknowns, most notably the total porosity, allow uncertainties in these mass and energy estimates to be as large as factors of two or three.     

Albedo and size determination of potentially hazardous asteroids: (99942) Apophis

The near-Earth object (99942) Apophis will make an extremely close approach to the Earth in 2029, and currently has approximately a one-in-45,000 chance of impacting our planet in 2036 (JPL Sentry, November 2006). Computation of the orbital evolution of this object is limited by insufficient knowledge of physical properties required to determine the role played by non-gravitational effects. Using polarimetric observations, we have obtained the first reliable determination of the albedo of Apophis, obtaining 0.33±0.08. We also derive an updated estimate of the asteroid's absolute magnitude: H=19.7±0.4. Using this albedo and H, we find that Apophis has a diameter of 270±60 m, slightly smaller than preliminary estimates based upon an assumed albedo. Our observations demonstrate the feasibility of polarimetric observations aimed at obtaining albedos and sizes of small, potentially hazardous asteroids.     

Treatment of star catalog biases in asteroid astrometric observations

In this paper, we discuss the detection of systematic biases in star positions of the USNO A1.0, A2.0, and B1.0 catalogs, as deduced from the residuals of numbered asteroid observations. We present a technique for the removal of these biases, and validate this technique by illustrating the resulting improvements in numbered asteroid residuals, and by establishing that debiased orbits predict omitted observations more accurately than do orbits derived from non-debiased observations. We also illustrate the benefits of debiasing to high-precision astrometric applications such as asteroid mass determination and collision analysis, including a refined prediction of the impact probability of 99942 Apophis. Specifically, we find the IP of Apophis to be lowered by nearly an order of magnitude to 4.5 × 10⁻⁶ for the 2036 close approach.     

The influence of the correlations between an asteroid’s orbital parameters on the estimation of the probability of planetary collision by the Monte Carlo method

The probability of an asteroid colliding with a planet can be estimated by the Monte Carlo method, in particular, through the statistical simulation of the possible initial conditions for the motion of an asteroid based on the probability density distribution set by the respective covariance matrix to be further projected with the orbital model onto the supposed time point of the collision. Hence, the collision probability is calculated as the ratio between the number of projected (virtual) asteroids striking the planet and their total number. The main problem is that different elements of the initial conditions (orbit or state vector) are correlated and, therefore, cannot be simulated independently. These correlations are reflected in the nondiagonal covariance matrix of the solution. The matrix is diagonalized by an orthogonal transformation. In the uncertainty domain constructed from the diagonal matrix elements, the initial values for each of the six orbital elements are simulated independently from the other elements, but with the accounting for their normal distribution. The program for calculating the normal distribution is based on the central limit theorem. Each sample of the initial values for the six orbital elements is transferred to the initial reference frame using an inverse transformation. Then, numerical integration is used to track the asteroid’s motion along the respective orbit to predict a possible impact event. Asteroids 99942 Apophis and 2007 WD5 are used as examples to show that disregarding the correlations when diagonalizing the covariance matrix to set the initial conditions may seriously distort the collision probability estimates. The paper gives the probabilities of the collisions of Apophis with the Earth and asteroid 2007 WD5 with Mars calculated by the author from observation sets showing nonzero collision probabilities. The author’s estimates are compared to those calculated by NASA.     

Apophis: the Story Behind the Scenes

On December 20, 2004 the Minor Planet Center issued the Minor Planet Electronic Circular (MPEC) 2004-Y25 announcing the discovery of a new Near Earth Asteroid (NEA) with designation 2004 MN₄. Only two days later, when the Christmas holidays were about to begin, it was already apparent that this asteroid, currently known as Apophis, would be notorious: our close-approach monitoring system, CLOMON2, was already showing a Virtual Impactor (VI) in 2029 reaching the level 2 in the Torino Scale, the first asteroid to reach that level since our monitoring system had been operational. However, this was just the beginning of what it was to come in the subsequent days. In this paper we will give an overview of the NEODyS-CLOMON2 system and provide the details on how Apophis’ collision scenario evolved, the way NEODyS’ team handled it and the crazy 2004’ Christmas holidays we had due to this unexpected guest.     

Asteroid close encounters characterization using differential algebra: the case of Apophis

A method for the nonlinear propagation of uncertainties in Celestial Mechanics based on differential algebra is presented. The arbitrary order Taylor expansion of the flow of ordinary differential equations with respect to the initial condition delivered by differential algebra is exploited to implement an accurate and computationally efficient Monte Carlo algorithm, in which thousands of pointwise integrations are substituted by polynomial evaluations. The algorithm is applied to study the close encounter of asteroid Apophis with our planet in 2029. To this aim, we first compute the high order Taylor expansion of Apophis’ close encounter distance from the Earth by means of map inversion and composition; then we run the proposed Monte Carlo algorithm to perform the statistical analysis.     

The orbit of asteroid (99942) Apophis as determined from optical and radar observations

The results of improving the orbit accuracy for the asteroid Apophis and the circumstances of its approach to Earth in 2029 are described. Gravitational perturbations from all of the major planets and Pluto, Ceres, Pallas, and Vesta are taken into account in the equations of motion of the asteroid. Relativistic perturbations from the Sun and perturbations due to the oblateness of the Sun and Earth and due to the light pressure are also included in the model. Perturbations from the Earth and Moon are considered separately. The coordinates of the perturbing bodies are calculated using DE405. The phase correction and the gravitational deflection of light are taken into account. The numerical integration of the equations of motion and equations in variations is performed by the 15th-order Everhart method. The error of the numerical integration over the 2005–2029 interval, estimated using forward and backward computations, is not more than 3 × 10^?11 AU. Improved coordinates and velocities at epoch JD2454200.5 (April 10, 2007) were obtained applying the weighted leastsquares fit. For the period from March 15, 2004, to August 16, 2006, 989 optical and 7 radar observations were used. The resulting system represents the optical observations with an error of 0.37 (66 conditional equations were rejected). The residuals of the radar observations are an order, or more, smaller than their errors. The system of Apophis’ elements and the estimates of their precision obtained in this study are in perfect agreement with the results published by other authors. The minimum Apophis-Earth distance is about 38 200 km on April 13, 2029. This estimate agrees to within 20 km with those calculated based on other published systems of elements. The effect of some model components on the minimum distance is estimated.     

Asteroid Apophis: Evaluating the impact hazards of such bodies

Soon after the discovery of asteroid 99942 Apophis, it was classified as a potentially hazardous object with a high probability of an impact on the Earth in 2029. Although subsequent observations have substantially reduced the probability of a collision, it has not been ruled out; moreover, similar-sized asteroids in orbits intersecting the Earth’s orbit may well be discovered in the near future. We conduct a numerical simulation of an atmospheric passage and an impact on the Earth’s surface of a stony cosmic body with a diameter of 300 m and kinetic energy of about 1000 Mt, which roughly corresponds to the parameters of the asteroid Apophis, at atmospheric entry angles of 90° (vertical stroke), 45°, and 30°. The simulation is performed by solving three-dimensional equations of hydrodynamics and radiative transfer equations in the approximations of radiative heat conduction and volume emission. The following hazards are considered: an air shock wave, ejecta from the crater, thermal radiation, and ionospheric disturbances. Our calculations of the overpressure and wind speed on the Earth’s surface show that the zone of destruction of the weakest structures can be as large as 700–1000 km in diameter; a decrease in the flight path angle to the surface leads to a marked increase in the area affected by the shock wave. The ionospheric disturbances are global in nature and continue for hours: at distances of several thousand kilometers at altitudes of more than 100 km, air density disturbances are tens of percent and the vertical and horizontal velocity components reach hundreds of meters per second. The impact of radiation on objects on the Earth’s surface is estimated by solving the equation of radiative transfer along rays passing through a luminous area. In clear weather, the size of the zone where thermal heating may ignite wood can be as large as 200 km, and the zone of individual fire outbreaks associated with the ignition of flammable materials can be twice as large. In the 100-km central area, which is characterized by very strong thermal damage, there is ignition of structures, roofs, clothes, etc. The human hazardous area increases with the decrease in the trajectory angle, and people may experience thermal effects at distances of up to 250–400 km from the crater.     

Properties of collision trajectories of asteroids with the earth

Close encounters of asteroids with the Earth can lead to a transition to resonant orbits, close encounters, and collisions in the future. The structure of the sets of possible collisions of asteroids with the Earth is similar to a fractal structure thanks to resonant returns of asteroids. Areas of semi-major axes and positions and sizes of gaps leading to collisions with the Earth for Apophis asteroids 2007 VK184 and 2011 AG5 are given.     

On the possibility of the guidance of small asteroids to dangerous celestial bodies using the gravity-assist maneuver

In this paper, the method of changing the trajectories of hazardous asteroids with orbits known for some years to be on a possible collision course with the Earth is considered. The method relies on the use of small asteroids (asteroid-projectiles) directed at hazardous celestial bodies by giving the projectile a sufficiently small velocity impulse ensuring the Earth gravity assist. As a result, the asteroid-projectile vector can be controllably changed over a wide range. Apophis is considered as an example of the target asteroid. The technical feasibility of this method is discussed. It is noted that despite the potential use of this elegant method, its practical implementation requires further research and development.     

Solar radiation pressure on (99942) Apophis

Near-Earth asteroid (99942) Apophis currently resides among the top positions on the list of objects with small, yet non-zero impact probability with the Earth. For that reason an unusual observational and theoretical effort has been dedicated to precisely characterize its future orbit. Here we discuss orbital perturbation of Apophis due to incident and reflected solar radiation pressure (SRP). We both revisit recent analytical estimate of the SRP effects for this body and also formulate a numerical approach allowing us to compute the SRP orbital perturbation under general assumptions. Contrary to some previous results, we show that SRP has a much smaller effect on the Apophis trajectory than does the thermal re-radiation force which produces the Yarkovsky effect. When the Yarkovsky effect becomes constrained enough in the future, our approach may be used to improve the orbit determination for this asteroid.     

On the accuracy of the orbit of asteroid (99942) APOPHIS at the time of its encounter with the Earth in 2029

We estimate the effect of trajectory measurement errors on the orbital parameters of asteroid Apophis determined from improvements. For this purpose, based on all of the optical and radar observations available to date, we have computed a nominal orbit of the asteroid. The scatter ellipsoid of the initial conditions of motion has been obtained by two methods. In the first, universally accepted method, the scatter ellipsoid is calculated by assuming a linear dependence of the errors in the parameters being determined on observational errors. In the second method, the scatter region of the orbital parameters around the nominal-orbit parameters is determined by the Monte Carlo method. We show that the region determined by the latter method at the initial epoch differs only slightly from the scatter ellipsoid for the linear approximation. We estimate the sizes of the projections of the corresponding regions onto the target plane at the time of the closest encounter of the asteroid with the Earth in 2029. The projections are approximated by ellipses. Our computations have shown that the ellipse has the following sizes: 389.6 km for the semimajor axis and 16.4 km for the semiminor axis in the linear case and 330.0 and 11.1 km, respectively, in the nonlinear case.