Dynamical evolution of NEAs: Close encounters, secular perturbations and resonances

We discuss the main mechanisms affecting the dynamical evolution of Near-Earth Asteroids (NEAs) by analyzing the results of three numerical integrations over 1 Myr of the NEA (4179) Toutatis. In the first integration the only perturbing planet is the Earth. So the evolution is dominated by close encounters and looks like a random walk in semimajor axis and a correlated random walk in eccentricity, keeping almost constant the perihelion distance and the Tisserand invariant. In the second integration Jupiter and Saturn are present instead of the Earth, and the 3/1 (mean motion) and v ₆ (secular) resonances substantially change the eccentricity but not the semimajor axis. The third, most realistic, integration including all the three planets together shows a complex interplay of effects, with close encounters switching the orbit between different resonant states and no approximate conservation of the Tisserand invariant. This shows that simplified 3-body or 4-body models cannot be used to predict the typical evolution patterns and time scales of NEAs, and in particular that resonances provide some fast-track dynamical routes from low-eccentricity to very eccentric, planet-crossing orbits.On leave from the Department of Mathematics, University of Pisa, Via Buonarroti 2, 56127 Pisa, Italy, thanks to the G. Colombo fellowships of the European Space Agency.     

Physical investigation of the potentially hazardous Asteroid (144898) 2004 VD17

In this paper we present the observational campaign carried out at ESO NTT and VLT in April and May 2006 to investigate the nature and the structure of the near-Earth object (144898) 2004 VD17. In spite of a great quantity of dynamical information, according to which it will have a close approach with the Earth in the next century, the physical properties of this asteroid are largely unknown. We performed visible and near-infrared photometry and spectroscopy, as well as polarimetric observations. Polarimetric and spectroscopic data allowed us to classify 2004 VD17 as an E-type asteroid. A good agreement was also found with the spectrum of the aubrite meteorite Mayo Belwa. On the basis of the polarimetric albedo (pv=0.45) and of photometric data, we estimated a diameter of about 320 m and a rotational period of about 2 h. The analysis of the results obtained by our complete survey have shown that (144898) 2004 VD17 is a peculiar NEO, since it is close to the breakup limits for fast rotator asteroids, as defined by Pravec and Harris [Pravec, P., Harris, A.W., 2000. Icarus 148, 12-20]. These results suggest that a more robust structure must be expected, as a fractured monolith or a rubble pile in a “strength regime” [Holsapple, K.A., 2002. Speed limits of rubble pile asteroids: Even fast rotators can be rubble piles. In: Workshop on Scientific Requirements for Mitigation of Hazardous Comets and Asteroids, Washington, September, 2002].     

Orbital and thermal evolutions of four potential targets for a sample return space mission to a primitive near-Earth asteroid

In this paper, we present our study of the orbital and thermal evolutions, due to solar radiative heating, of four near-Earth asteroids (NEAs) considered as potential target candidates for sample return space missions to primitive asteroids. We used a dynamical model of the NEA population to estimate the most likely source region and orbital history of these objects. Then, for each asteroid, we integrated numerically over their entire lifetime a set of 14 initially indistinguishable orbit (clones), obtained by small variations of the nominal initial conditions. Using a thermal model, we then computed surface and sub-surface temperatures of these bodies during their dynamical history. Our aim is to determine whether these bodies are likely to have experienced high temperature level, and whether great temperature changes can be expected due to the orbital changes as well as their maximum and minimum values. Such information is important in the framework of sample return space missions whose goal is to bring back pristine materials. The knowledge of the temperature range of materials at different depth over the orbital evolution of potential targets can help defining sampling strategies that ensure the likelihood that unaltered material will be brought back. Our results suggest that for all the considered potential targets, the surface has experienced for some time temperatures greater than 400 K and at most 500 K with 50% probability. This probability drops rapidly with increasing temperature. Sub-surface materials at a depth of only 3 cm are much more protected from high temperature and generally do not reach temperatures exceeding 450 K (with 50% probability). They should thus be unaltered at this depth at least from a Sun-driven heating point of view. On the other hand, surface material for some of the considered objects can have a range of temperature which can make them less reliable as pristine materials. However, it is assumed here that the same material is constantly exposed to solar heat, while regolith turnover may occur. The latter can be caused by different processes such as seismic shaking and/or impact cratering. This would reduce the total time that materials are exposed to a certain temperature. Thus, it is very likely that a sample collected from any of the four considered targets, or any primitive NEA with similar dynamical properties, will have components that will be thermally unaltered as long as some of it comes from only 3 to 5 cm depth. Such a depth is not considered difficult to reach with some of the current designs of sampling devices.     

Dynamical passage to approximate equilibrium shapes for spinning, gravitating rubble asteroids

Many asteroids are thought to be particle aggregates held together principally by self-gravity. Here we study — for static and dynamical situations — the equilibrium shapes of spinning asteroids that are permitted for rubble piles. As in the case of spinning fluid masses, not all shapes are compatible with a granular rheology. We take the asteroid to always be an ellipsoid with an interior modeled as a rigid-plastic, cohesion-less material with a Drucker-Prager yield criterion. Using an approximate volume-averaged procedure, based on the classical method of moments, we investigate the dynamical process by which such objects may achieve equilibrium. We first collapse our dynamical approach to its statical limit to derive regions in spin-shape parameter space that allow equilibrium solutions to exist. At present, only a graphical illustration of these solutions for a prolate ellipsoid following the Drucker-Prager failure law is available [Sharma, I., Jenkins, J.T., Burns, J.A., 2005a. Bull. Am. Astron. Soc. 37, 643; Sharma, I., Jenkins, J.T., Burns, J.A., 2005b. Equilibrium shapes of ellipsoidal soil asteroids. In: García-Rojo, R., Hermann, H.J., McNamara, S. (Eds.), Proceedings of the 5th International Conference on Micromechanics of Granular Media, vol. 1. A.A. Balkema, UK; Holsapple, K.A., 2007. Icarus 187, 500-509]. Here, we obtain the equilibrium landscapes for general triaxial ellipsoids, as well as provide the requisite governing formulae. In addition, we demonstrate that it may be possible to better interpret the results of Richardson et al. [Richardson, D.C., Elankumaran, P., Sanderson, R.E., 2005. Icarus 173, 349-361] within the context of a Drucker-Prager material. The graphical result for prolate ellipsoids in the static limit is the same as those of Holsapple [Holsapple, K.A., 2007. Icarus 187, 500-509] because, when worked out, his final equations will match ours. This is because, though the formalisms to reach these expressions differ, in statics, at the lowest level of approximation, volume-averaging and the approach of Holsapple [Holsapple, K.A., 2007. Icarus 187, 500-509] coincide. We note that the approach applied here was obtained independently [Sharma, I., Jenkins, J.T., Burns, J.A., 2003. Bull. Am. Astron. Soc. 35, 1034; Sharma, I., 2004. Rotational Dynamics of Deformable Ellipsoids with Applications to Asteroids. Ph.D. thesis, Cornell University] and it provides a general, though approximate, framework that is amenable to systematic improvements and is flexible enough to incorporate the dynamical effects of a changing shape, different rheologies and complex rotational histories. To demonstrate our technique, we investigate the non-equilibrium dynamics of rigid-plastic, spinning, prolate asteroids to examine the simultaneous histories of shape and spin rate for rubble piles. We have succeeded in recovering most results of Richardson et al. [Richardson, D.C., Elankumaran, P., Sanderson, R.E., 2005. Icarus 173, 349-361], who obtained equilibrium shapes by studying numerically the passage into equilibrium of aggregates containing discrete, interacting, frictionless, spherical particles. Our mainly analytical approach aids in understanding and quantifying previous numerical simulations.     

Mostly Dormant Comets and their Disintegration into Meteoroid Streams: A Review

The history of associating meteor showers with asteroidal-looking objects is long, dating to before the 1983 discovery that 3200 Phaethon moves among the Geminids. Only since the more recent recognition that 2003 EH1 moves among the Quadrantids are we certain that dormant comets are associated with meteoroid streams. Since that time, many orphan streams have found parent bodies among the newly discovered Near Earth Objects. The seven established associations pertain mostly to showers in eccentric or highly inclined orbits. At least 35 other objects are tentatively linked to streams in less inclined orbits that are more difficult to distinguish from those of asteroids. There is mounting evidence that the streams originated from discrete breakup events, rather than long episodes of gradual water vapor outgassing. If all these associations can be confirmed, they represent a significant fraction of all dormant comets that are in near-Earth orbits, suggesting that dormant comets break at least as frequently as the lifetime of the streams. I find that most pertain to NEOs that have not yet fully decoupled from Jupiter. The picture that is emerging is one of rapid disintegration of comets after being captured by Jupiter, and consequently, that objects such as 3200 Phaethon most likely originated from among the most primitive asteroids in the main belt, instead. They too decay mostly by disintegration into comet fragments and meteoroid streams. The disintegration of dormant comets is likely the main source of our meteor showers and the main supply of dust to the zodiacal cloud. Editorial handling: Frans Rietmeijer.     

Orbital and mission planning constraints for the deflection of NEOs impacting on Earth

This paper is the third in a series. Paper 1 presented the results of numerical modeling of deflections of NEOs in route of collision with the Earth. The model was applied to a variety of dynamical cases including both asteroidal and cometary NEOs. Paper 2 introduced the concept of “distributed deflection,” i.e., the possibility to provide the ΔV necessary to deflect an object with a succession of maneuvers each of which would have been insufficient per se to obtain the desired result. In both papers no assumptions were made on the physical composition and structure of the NEO, nor on the details of the possible deflection maneuvers from the point of view of mission analysis. Moreover, ΔV-plots were computed assuming only along-track impulses (both in the positive and negative directions), because it is easy to demonstrate that in general this is energetically the most favorable configuration. Also in the present paper no assumptions were made on the physical composition and structure of the NEO, even if order of magnitude considerations are made on the physical feasibility of a deflection, in terms of the internal strength of the NEO. We present here the results of an investigation on the mission requirements necessary to deflect an object (or contribute to a succession of deflecting maneuvers) in terms of accessibility of the spacecraft terminal orbit from Earth with the current launchers.     

Spin rate distribution of small asteroids

The spin rate distribution of main belt/Mars crossing (MB/MC) asteroids with diameters 3-15 km is uniform in the range from f=1 to 9.5 d⁻¹, and there is an excess of slow rotators with f<1 d⁻¹. The observed distribution appears to be controlled by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. The magnitude of the excess of slow rotators is related to the residence time of slowed down asteroids in the excess and the rate of spin rate change outside the excess. We estimated a median YORP spin rate change of ≈0.022 d⁻¹/Myr for asteroids in our sample (i.e., a median time in which the spin rate changes by 1 d⁻¹ is ≈45 Myr), thus the residence time of slowed down asteroids in the excess is ≈110 Myr. The spin rate distribution of near-Earth asteroids (NEAs) with sizes in the range 0.2-3 km (∼5 times smaller in median diameter than the MB/MC asteroids sample) shows a similar excess of slow rotators, but there is also a concentration of NEAs at fast spin rates with f=9-10 d⁻¹. The concentration at fast spin rates is correlated with a narrower distribution of spin rates of primaries of binary systems among NEAs; the difference may be due to the apparently more evolved population of binaries among MB/MC asteroids.     

近地小行星探测目标选择

从数以万计的近地小行星(NEAs)中选择探测目标,是一项非常重要的工作.选择小行星时需要考虑的主要问题是:NEAs的轨道参数、自旋周期、大小、科学价值、表面特征以及到达NEAs的时间.本文采用排除法从众多近地小行星中确定候选者,并根据科学目标的要求确定最终探测目标.     

Rotational fission of contact binary asteroids

The energetics and dynamics of contact binary asteroids as they approach and pass the rotational fission limit is studied. We presume that the asteroids are subject to an external torque, such as from the YORP effect, that increases their angular momentum. Furthermore, we assume the asteroids can be described by a fairly minimal model comprised of a sphere and ellipsoid resting on each other. The minimum energy configurations for contact binary asteroids at different levels of angular momentum are computed and discussed. We find distinct transitions between different configurations as the angular momentum of the system is increased. These indicate that rapidly rotating contact binary asteroids may seek out clearly different relative configurations than slowly rotating systems. We find a single end state of the systems prior to rotational fission, and distinct dynamical outcomes as a function of mass distribution and shape when the rotational fission limit is exceeded. Our theoretical results agree qualitatively with observed properties of near-Earth asteroids, and can be used to help explain the spin-rate barrier, contact binaries, and the observed morphology of most NEO binaries.