Collision Probability for Earth-Crossing Asteroids Using Orbital Ranging

We introduce new techniques for the computation of the collision probability for Earth-crossing asteroids in the case of short observational arcs and/or small numbers of observations. The techniques rely on the orbital element probability density computed using statistical orbital ranging. We apply the techniques to the Earth-crossing asteroid 1998 OX₄with non-vanishing collision probability in numerous close approaches after the year 2012 (inclusive). We study the invariance of the collision probability in transformations between different orbital element sets, and develop a Spearman rank correlation measure for the validity of the linear approximation. We introduce an optimized, fast version of the statistical ranging method.     

Asteroid orbits using phase-space volumes of variation

We present a statistical orbit computation technique for asteroids with transitional observational data, that is, a moderate number of data points spanning a moderate observational time interval. With the help of local least-squares solutions in the phase space of the orbital elements, we map the volume of variation as a function of one or more of the elements. We sample the resulting volume using a Monte Carlo technique and, with proper weights for the sample orbital elements, characterize the six-dimensional orbital-element probability density function. The volume-of-variation (VOV) technique complements the statistical ranging technique for asteroids with exiguous observational data (short time intervals and/or small numbers of observations) and the least-squares technique for extensive observational data. We show that, asymptotically, results using the new technique agree closely with those from ranging and least squares. We apply the technique to the near-Earth object 2004 HA₃₉, the main-belt object 2004 QR and the transneptunian object 2002 CX₂₂₄ recently observed at the Nordic Optical Telescope on La Palma, illustrating the potential of the technique in ephemeris prediction. The VOV technique helps us assess the phase transition in orbital-element probability densities, that is, the non-linear collapse of wide orbital-element distributions to narrow localized ones. For the three objects above, the transition takes place for observational time intervals of the order of 10 h, 5 d and 10 months, respectively, emphasizing the significance of the orbital-arc fraction covered by the observations.     

Asteroid identification over apparitions

We present a new method for the linking of scarce asteroid astrometry over apparitions, and apply it both to simulated and real data to prove its feasibility. Up to date, there has not been a robust method available to search for linkages between the approximately 50,000 provisionally designated sets of asteroid astrometry spanning less than two days. Unless such a scarce set of astrometry is linked to another set of astrometry, the underlying object can be considered lost as the ephemeris uncertainties are substantial. The new method, which can tackle the challenges, is based on Ranging, which is a fully nonlinear, statistical orbital inversion method. Ranging properly treats astrometric uncertainties and propagates the uncertainty to the resulting orbital-element probability density, which is sampled by a set of orbits. The new orbital-element-space multiple-address-comparison (oMAC) method uses dimensionality-reduction techniques and tree structures to efficiently search for overlapping probability densities in the orbital-element phase space. Overlapping probability densities indicate a candidate linkage between astrometric observation sets. To accept a candidate linkage, we have to find a many-body orbital solution which reproduces the observed positions within the observational uncertainties. To find the linking orbit, we use a multi-step approach starting from a Monte-Carlo generation of possible orbits in a reduced volume of the orbital-element phase space and ending with a least-squares orbital solution, which, in addition to the Sun's gravitation, also takes into account the gravitational influence of the relevant planets. The new multiple-address-comparison method has a loglinear computational complexity, that is, it scales as O(nlogn), where n is the number of included observation sets. It has recently also been implemented for the ephemeris-space multiple-address-comparison (eMAC) method, which is optimized for the short-term linking of scarce astrometry.     

OpenOrb: Open‐source asteroid orbit computation software including statistical ranging

Abstract— We are making an open‐source asteroid orbit computation software package called OpenOrb publicly available. OpenOrb is built on a well‐established Bayesian inversion theory, which means that it is to a large part complementary to orbit‐computation packages currently available. In particular, OpenOrb is the first package that contains tools for rigorously estimating the uncertainties resulting from the inverse problem of computing orbital elements using scarce astrometry. In addition to the well‐known least‐squares method, OpenOrb also contains both Monte‐Carlo (MC) and Markov‐Chain MC (MCMC; Oszkiewicz et al. [2009]) versions of the statistical ranging method. Ranging allows the user to obtain sampled, non‐Gaussian orbital‐element probability‐density functions and is therefore optimized for cases where the amount of astrometry is scarce or spans a relatively short time interval. Ranging‐based methods have successfully been applied to a variety of different problems such as rigorous ephemeris prediction, orbital element distribution studies for transneptunian objects, the computation of invariant collision probabilities between near‐Earth objects and the Earth, detection of linkages between astrometric asteroid observations within an apparition as well as between apparitions, and in the rigorous analysis of the impact of orbital arc length and/or astrometric uncertainty on the uncertainty of the resulting orbits. Tools for making ephemeris predictions and for classifying objects based on their orbits are also available in OpenOrb. As an example, we use OpenOrb in the search for candidate retrograde and/or high‐inclination objects similar to 2008 KV₄₂ in the known population of transneptunian objects that have an observational time span shorter than 30 days.     

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.     

Intrinsic nonlinearity and method of disturbed observations in inverse problems of celestial mechanics

Orbit determination from a small sample of observations over a very short observed orbital arc is a strongly nonlinear inverse problem. In such problems an evaluation of orbital uncertainty due to random observation errors is greatly complicated, since linear estimations conventionally used are no longer acceptable for describing the uncertainty even as a rough approximation. Nevertheless, if an inverse problem is weakly intrinsically nonlinear, then one can resort to the so-called method of disturbed observations (aka observational Monte Carlo). Previously, we showed that the weaker the intrinsic nonlinearity, the more efficient the method, i.e. the more accurate it enables one to simulate stochastically the orbital uncertainty, while it is strictly exact only when the problem is intrinsically linear. However, as we ascertained experimentally, its efficiency was found to be higher than that of other stochastic methods widely applied in practice. In the present paper we investigate the intrinsic nonlinearity in complicated inverse problems of Celestial Mechanics when orbits are determined from little informative samples of observations, which typically occurs for recently discovered asteroids. To inquire into the question, we introduce an index of intrinsic nonlinearity. In asteroid problems it evinces that the intrinsic nonlinearity can be strong enough to affect appreciably probabilistic estimates, especially at the very short observed orbital arcs that the asteroids travel on for about a hundredth of their orbital periods and less. As it is known from regression analysis, the source of intrinsic nonlinearity is the nonflatness of the estimation subspace specified by a dynamical model in the observation space. Our numerical results indicate that when determining asteroid orbits it is actually very slight. However, in the parametric space the effect of intrinsic nonlinearity is exaggerated mainly by the ill-conditioning of the inverse problem. Even so, as for the method of disturbed observations, we conclude that it practically should be still entirely acceptable to adequately describe the orbital uncertainty since, from a geometrical point of view, the efficiency of the method directly depends only on the nonflatness of the estimation subspace and it gets higher as the nonflatness decreases.     

Modeling orbital gamma‐ray spectroscopy experiments at carbonaceous asteroids

To evaluate the feasibility of measuring differences in bulk composition among carbonaceous meteorite parent bodies from an asteroid or comet orbiter, we present the results of a performance simulation of an orbital gamma‐ray spectroscopy (GRS) experiment in a Dawn‐like orbit around spherical model asteroids with a range of carbonaceous compositions. The orbital altitude was held equal to the asteroid radius for 4.5 months. Both the asteroid gamma‐ray spectrum and the spacecraft background flux were calculated using the MCNPX Monte‐Carlo code. GRS is sensitive to depths below the optical surface (to ≈20–50 cm depth depending on material density). This technique can therefore measure underlying compositions beneath a sulfur‐depleted (e.g., Nittler et al. 2001 ) or desiccated surface layer. We find that 3σ uncertainties of under 1 wt% are achievable for H, C, O, Si, S, Fe, and Cl for five carbonaceous meteorite compositions using the heritage Mars Odyssey GRS design in a spacecraft‐deck‐mounted configuration at the Odyssey end‐of‐mission energy resolution, FWHM = 5.7 keV at 1332 keV. The calculated compositional uncertainties are smaller than the compositional differences between carbonaceous chondrite subclasses.     

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.     

Orbit computation for transneptunian objects

Using statistical orbital ranging, we systematically study the orbit computation problem for transneptunian objects (TNOs). We have automated orbit computation for large numbers of objects, and, more importantly, we are able to obtain orbits even for the most sparsely observed objects (observational arcs of a few days). For such objects, the resulting orbit distributions include a large number of high-eccentricity orbits, in which TNOs can be perturbed by close encounters with Neptune. The stability of bodies on the computed orbits has therefore been ascertained by performing a study of close encounters with the major planets. We classify TNO orbit distributions statistically, and we study the evolution of their ephemeris uncertainties. We find that the orbital element distributions for the most numerous single-apparition TNOs do not support the existence of a postulated sharp edge to the belt beyond 50 AU. The technique of statistical ranging provides ephemeris predictions more generally than previously possible also for poorly observed TNOs.     

Impact orbits of the asteroid 2009 FJ with the Earth

We show how to calculate the impact orbits of dangerous asteroids using the freely available the OrbFit software, and compare our results with impact orbits calculated using Sitarski??s independent software (Sitarski, 1999; 2000; 2006). The new method is tested on asteroid 2009 FJ. Using the OrbFit package to integrate alternate orbits along the line of variation (Milani et al., 2002; 2005a; 2005b), we identify impact orbits and can plot paths of risk for the Earth or any other body in the Solar System. We present the orbital elements of asteroid 2009 FJ and its ephemerides, along with uncertainties, for the next 100 years. This paper continues a long-term research program on impact solutions for asteroids (Wlodarczyk, 2007; 2008; 2009).     

Cosmic Roulette: Comets In The Main Belt Asteroid Region

We have produced top ten ranked lists of impact velocity, mainbelt asteroid region dwell times and impact probabilities for a selection of short period comets. The comet with the combined highest ranking with respect to impact probability and impact velocity is Comet C/1766 G1 Helfenzrieder. Since it is not clear that this comet still exists, the highest ranked, presently active, comet with respect to the likelihood of suffering impacts from meter-sized objects while in the main belt asteroid region is Comet 28P/Neujmin 1. We find no evidence to support the existence of a distinctive sub-set of the short period comets liable to show repeated outburst or splitting behavioursdue to small body, meter-sized, asteroid impacts. This revised version was published online in July 2006 with corrections to the Cover Date.     

Neptune Trojans: a Revisit and a New Membertwo

Up to now, 17 Neptune Trojan asteroids have been detected with their orbits being well determined by continuous observations. This paper analyzes systematically their orbital dynamics. Our results show that except for two temporary members with relatively short lifespans on Trojan orbits, the vast majority of Neptune Trojans located within their orbital uncertainties may survive in the solar system age. The escaping probability of Neptune Trojans, through slow diffusion in the orbital element space in 4.5 billion years, is estimated to be ～50%. The asteroid 2012 UW177 classified as a Centaur asteroid by the IAU Minor Planet Center currently is in fact a Neptune Trojan. Numerical simulations indicate that it is librating on the tadpole-shaped orbit around the Neptune's L₄ point. It was captured into the current orbit approximately 0.23 million years ago, and will stay there for at least another 1.3 million years in the future. Its high inclination of i ≈ 54^° not only makes it the most inclined Neptune Trojan, but also makes it exhibit the complicated and interesting co-orbital transitions between the leading and trailing Trojans via the quasi-satellite orbit phase.     

Development of an observational error model

In calculating the orbit of a minor planet with a least-squares algorithm, current practice is to assume that all observations of a given era have the same uncertainty, and that the errors in these observations are uncorrelated. These assumptions are unrealistic; and they lead to sub-optimal orbits.Our objective is to develop and validate an observational error model that provides realistic estimates of the uncertainties and correlations in asteroid observations. When used to populate the covariance matrix of the least-squares algorithm, the resulting orbits are shown to more accurately and precisely represent asteroid trajectories.     

A Method for the Determination of an Intermediate Orbit from Three Positions of the Small Body on the Celestial Sphere

A new method is suggested for finding the preliminary orbit from three complete measurements of the angular coordinates of a celestial body developed by analogy with the classic Lagrange–Gauss method. The proposed method uses the intermediate orbit that we had constructed in an earlier paper based on two position vectors and the corresponding time interval. This intermediate orbit allows for most of the perturbations in the motion of the body. Using the orbital motion of asteroid 1566 Icarus as an example, we compare the results obtained by applying the classic and the new method. The comparison shows the new method to be highly efficient for studying perturbed motion. It is especially efficient if applied to high-precision observational data covering short orbital arcs.     

Radar observations of asteroid 1998 ML14

Abstract— Goldstone and Arecibo delay‐Doppler radar imaging of asteroid 1998 ML 14 shortly after its discovery reveals a 1 km diameter spheroid with prominent topography on one side and subdued topography on the other. The object's radar and optical properties are typical for S‐class near‐Earth asteroids. The gravitational slopes of a shape model derived from the images and assumed to have a uniform density are shallow, exceeding 30° over only 4% of the surface. If 1998 ML14's density distribution is uniform, then its orbital environment is similar to a planetary body with a spheroidal gravitational field and is relatively stable. Integration of a radar‐refined orbit reveals that the 1998 apparition was the asteroid's closest approach to Earth since at least 1100 and until 2283, when it approaches to within 2.4 lunar distances. Outside of that time interval, orbit uncertainties based on the present set of observations preclude reliable prediction.     

Transneptunian Object Ephemeris Service (TNOEPH)

We present a web service called TNOEPH ( http://asteroid.lowell.edu/)for ephemeris uncertainty prediction and dynamical classification ofshort-arc trans\-neptunian objects (TNOs). User-supplied observations are transformed to a rigorous sky-plane uncertainty map using the technique of statistical orbital ranging. We show examples of thegrowth of ephemeris uncertainty with time, and give the probabilities ofdifferent dynamical classifications for a few short-arc TNOs.     

Peculiarities of the motion of asteroid 99942 Apophis

The Apophis asteroid attracted the attention of scientists immediately after its discovery in 2004, because the initially determined orbit of this asteroid assumes a possible collision with Earth in April 2029. The size of Apophis is about several hundred meters, and its collision with Earth might result in a large regional or even global catastrophe. At present, the trajectory of Apophis has been calculated more accurately, and a collision in 2029 has ruled out; the asteroid will pass Earth at a distance of about 37 000 km from its center. However, close approaches or collisions are possible after 2029, including the most probable in 2036. The risk of collision in 2036 is well known and actively examined by the scientists. In this study, we consider the peculiarities of the asteroid motion associated with its approach in 2029 and with a possible close approach in 2036. The trajectories scatter during the approaches and the loss of accuracy is associated with these scatterings. As a result, the trajectory of Apophis may become nondeterministic after 2036; that is, it cannot now be determined unambiguously. Although such events are very unlikely, it is interesting to examine a variety of alternative variants of Apophis’ close approaches and collisions with Earth immediately after 2036. The effects of small variations in the asteroid velocity at different moments in time after its impact with a certain mass are discussed.     

Application of the transformation of variables technique for uncertainty mapping in nonlinear filtering

This paper addresses the impact nonlinear observations of state variables have on uncertainty accuracy associated with state estimation algorithms. The transformation of variables technique is applied to exactly map probability density functions (PDFs) between domains completely spanned by different combinations of basis vectors. The technique allows for proper generation of the likelihood density when converting from measurement to state variable space and for association of a present state distribution with prior observation data. The exact mapping of probability distribution functions between domains and proper characterization of prior knowledge allows for Bayesian estimation to be appropriately carried out. A Bayes filter utilizing the technique is developed which uses the technique to map the uncertainty in time for generation of the prior density and in space for generation of the likelihood density. The filter is compared with conventional nonlinear filtering techniques in multiple scenarios to demonstrate the utility and insight offered for object tracking and parameter estimation applications.     

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.