Radar and optical observations and physical modeling of near-Earth Asteroid 10115 (1992 SK)

We estimate Asteroid 1992 SK's physical properties from delay-Doppler images and Doppler-only echo spectra obtained during March 22-27, 1999, at Goldstone and from optical lightcurves obtained during February-March 1999 at Ond?ejov Observatory. The images span only about 15° of sky motion and are not strong, but they place up to twenty 40 m by 160 m pixels on the asteroid and have complete rotational phase coverage. Our analysis establishes that the radar observations are confined to subradar latitudes between −20° and −40°. The echo spectra and optical lightcurves span ∼80° of sky motion, which provides important geometric leverage on the pole direction. The lightcurves are essential for accurate estimation of the asteroid's shape and spin state. We estimate the asteroid's period to be 7.3182±0.0003 h and its pole direction to be at ecliptic longitude, latitude=(99°±5°,−3°±5°). The asteroid is about 1.4 km in maximum extent and mildly asymmetric, with an elongation of about 1.5 and relatively subdued topography. The OC radar albedo is 0.11±0.02 and the SC/OC ratio is 0.34±0.05. The current orbital solution permits accurate identification of planetary close approaches during 826-2690. We use our model to predict salient characteristics of radar images and optical lightcurves obtainable during the asteroid's March 2006 approach.     

A steady-state model of NEA binaries formed by tidal disruption of gravitational aggregates

We present results of a simulation of a steady-state binary near-Earth asteroid (NEA) population. This study combines previous work on tidal disruption of gravitational aggregates [Walsh, K.J., Richardson, D.C., 2006. Icarus 180, 201-216] with a Monte Carlo simulation of NEA planetary encounters. Evolutionary effects include tidal evolution and binary disruption from close planetary encounters. The results show that with the best known progenitor (small Main Belt asteroids) shape and spin distributions, and current estimates of NEA lifetime and encounter probabilities, that tidal disruption should account for approximately 1-2% of NEAs being binaries. Given the best observed estimate of a ∼15% binary NEA fraction, we conclude that there are other formation mechanisms that contribute significantly to this population. We also present the expected distribution of binary orbital and physical properties for the steady-state binary NEAs formed by tidal disruption. We discuss the effects on binary fraction and properties due to changes in the least constrained parameters, and other possible effects on our model that could account for differences between the presented results and the observed binary population. Finally, we model possible effects of a significant population of binaries migrating to the near-Earth population from the Main Belt.     

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.     

Optimal deflection of NEOs en route of collision with the Earth

Recently, a method for the n-body analysis of the velocity change required to deflect a hazardous near-Earth object (NEO) was presented by Carusi et al. [Carusi, A., Valsecchi, G.B., D'Abramo, G., Boattini A., 2002. Icarus 159, 417-422]. We extent this method in order to optimize the velocity change vector instead of its along-track magnitude. From an application of both methods to a fictitious NEO we find Carusi's parallel approach to be reasonable for phases of unperturbed two-body motion. But, for orbit phases inhering third-body perturbations, i.e., for planetary close approaches or prior to a collision, the results obtained from the new method show the radial component of deflection impulse to play a major role. We show that a fivefold greater efficiency can be achieved by a deflection impulse being non-parallel to orbital velocity. The new method is applied to two possible 99942 Apophis impact trajectories in order to provide constraints for future Apophis deflection mission analysis.     

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.     

Physical properties of near-Earth Asteroid (33342) 1998 WT24

During its close Earth approach in 2001, the E-class near-Earth Asteroid (33342) 1998 WT24 was the focus of extensive radar, optical, and thermal infrared observations. We present a physical model of this object, estimated from Arecibo and Goldstone radar images that cover multiple rotations and span over 100° of sky motion. The asteroid has an equivalent diameter of 415±40 m and a diffuse radar scattering law that is identical in both senses of circular polarization, implying a surface that is extremely rough on centimeter-to-decimeter scales. The shape is dominated by three large basins, which may be impact craters or a relic of past dynamical disruption of the object. Analysis of YORP perturbations on WT24's spin state predicts that the asteroid's spin rate is decreasing at a rate of . Simply extrapolating this rate suggests that the asteroid will despin over the next 150 kyr and was spinning at its surface disruption rate 75 kyr ago, but the rotational evolution of WT24 depends on the surface's thermal properties and probably is more complex than a simple spin-down.     

The spin state of 433 Eros and its possible implications

In this paper, we show that Asteroid (433) Eros is currently residing in a spin-orbit resonance, with its spin axis undergoing a small-amplitude libration about the Cassini state 2 of the proper mode in the nonsingular orbital element sinI/2exp(?Ω), where I the orbital inclination and Ω the longitude of the node. The period of this libration is ?53.4 kyr. By excluding these libration wiggles, we find that Eros' pole precesses with the proper orbital plane in inertial space with a period of ?61.4 kyr. Eros' resonant state forces its obliquity to oscillate with a period of ?53.4 kyr between ?76° and ?89.5°. The observed value of ?89° places it near the latter extreme of this cycle. We have used these results to probe Eros' past orbit and spin evolution. Our computations suggest that Eros is unlikely to have achieved its current spin state by solar and planetary gravitational perturbations alone. We hypothesize that some dissipative process such as thermal torques (e.g., the so-called YORP effect) may be needed in our model to obtain a more satisfactory match with data. A detailed study of this problem is left for future work.     

The primordial excitation and clearing of the asteroid belt—Revisited

We have performed new simulations of two different scenarios for the excitation and depletion of the primordial asteroid belt, assuming Jupiter and Saturn on initially circular orbits as predicted by the Nice Model of the evolution of the outer Solar System [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465]. First, we study the effects of sweeping secular resonances driven by the depletion of the solar nebula. We find that these sweeping secular resonances are incapable of giving sufficient dynamical excitation to the asteroids for nebula depletion timescales consistent with estimates for solar-type stars, and in addition cannot cause significant mass depletion in the asteroid belt or produce the observed radial mixing of different asteroid taxonomic types. Second, we study the effects of planetary embryos embedded in the primordial asteroid belt. These embedded planetary embryos, combined with the action of jovian and saturnian resonances, can lead to dynamical excitation and radial mixing comparable to the current asteroid belt. The mass depletion driven by embedded planetary embryos alone, even in the case of an eccentric Jupiter and Saturn, is roughly 10-20× less than necessary to explain the current mass of the main belt, and thus a secondary depletion event, such as that which occurs naturally in the Nice Model, is required. We discuss the implications of our new simulations for the dynamical and collisional evolution of the main belt.     

Radar observations and a physical model of Asteroid 1580 Betulia

We report Arecibo (2380-MHz, 13-cm) observations of Asteroid 1580 Betulia in May-June 2002. We combine these continuous-wave Doppler spectra and delay-Doppler images with optical lightcurves from the 1976 and 1989 apparitions in order to estimate Betulia's shape and spin vector. We confirm the spin vector solution of Kaasalainen et al. [Kaasalainen, M., and 21 colleagues, 2004. Icarus 167, 178-196], with sidereal period P=6.13836 h and ecliptic pole direction (λ,β)=(136°,+22°), and obtain a model that resembles the Kaasalainen et al. convex-definite shape reconstruction but is dominated by a prominent concavity in the southern hemisphere. We find that Betulia has a maximum breadth of 6.59±0.66 km and an effective diameter of 5.39±0.54 km. These dimensions are in accord with reanalyzed polarimetric and radar data from the 1970s. Our effective diameter is 15% larger than the best radiometric estimate of Harris et al. [Harris, A.W., Mueller, M., Delbó, M., Bus, S.J., 2005. Icarus 179, 95-108], but this difference is much smaller than the size differences between past models. Considering orbits of test particles around Betulia, we find that this asteroid's unusual shape results in six equilibrium points close to its equatorial plane rather than the usual four points; two of these six points represent stable synchronous orbits while four are unstable. Betulia's close planetary encounters can be predicted for over four thousand years into the future.     

Fragmentation model dependence of collision cascades

Mass depletion of bodies through successive collisional disruptions (i.e., collision cascade) is one of the most important processes in the studies of the asteroids belt, the Edgeworth-Kuiper belt, debris disks, and planetary formation. The collisional disruption is divided into two types, i.e., catastrophic disruption and cratering. Although some studies of the collision cascades neglected the effect of cratering, it is unclear which type of disruption makes a dominant contribution to the collision cascades. In the present study, we construct a simple outcome model describing both catastrophic disruption and cratering, which has some parameters characterizing the total ejecta mass, the mass of the largest fragment, and the power-law exponent of the size distribution of fragments. Using this simple outcome model with parameters, we examine the model dependence of the mass depletion time in collision cascades for neglect of coalescence of colliding bodies due to high collisional velocities. We find the cratering collisions are much more effective in collision cascades than collisions with catastrophic disruption in a wide region of the model parameters. It is also found that the mass depletion time in collision cascades is mainly governed by the total ejecta mass and almost insensitive to the mass of the largest fragment and the power-law exponent of fragments for a realistic parameter region. The total ejecta mass is usually determined by the ratio of the impact energy divided by the target mass (i.e. Q-value) to its threshold value for catastrophic disruption, as well as in our simple model. We derive a mass depletion time in collision cascades, which is determined by of the high-mass end of collision cascades. The mass depletion time derived with our model would be applicable to debris disks and planetary formation.     

Linking the collisional history of the main asteroid belt to its dynamical excitation and depletion

The main belt is believed to have originally contained an Earth mass or more of material, enough to allow the asteroids to accrete on relatively short timescales. The present-day main belt, however, only contains ∼5×¹⁰⁻⁴ Earth masses. Numerical simulations suggest that this mass loss can be explained by the dynamical depletion of main belt material via gravitational perturbations from planetary embryos and a newly-formed Jupiter. To explore this scenario, we combined dynamical results from Petit et al. [Petit, J. Morbidelli, A., Chambers, J., 2001. The primordial excitation and clearing of the asteroid belt. Icarus 153, 338-347] with a collisional evolution code capable of tracking how the main belt undergoes comminution and dynamical depletion over 4.6 Gyr [Bottke, W.F., Durda, D., Nesvorny, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H., 2005. The fossilized size distribution of the main asteroid belt. Icarus 175, 111-140]. Our results were constrained by the main belt's size-frequency distribution, the number of asteroid families produced by disruption events from diameter D>100 km parent bodies over the last 3-4 Gyr, the presence of a single large impact crater on Vesta's intact basaltic crust, and the relatively constant lunar and terrestrial impactor flux over the last 3 Gyr. We used our model to set limits on the initial size of the main belt as well as Jupiter's formation time. We find the most likely formation time for Jupiter was 3.3±2.6 Myr after the onset of fragmentation in the main belt. These results are consistent with the estimated mean disk lifetime of 3 Myr predicted by Haisch et al. [Haisch, K.E., Lada, E.A., Lada, C.J., 2001. Disk frequencies and lifetimes in young clusters. Astrophys. J. 553, L153-L156]. The post-accretion main belt population, in the form of diameter D?1000 km planetesimals, was likely to have been 160±40 times the current main belt's mass. This corresponds to 0.06-0.1 Earth masses, only a small fraction of the total mass thought to have existed in the main belt zone during planet formation. The remaining mass was most likely taken up by planetary embryos formed in the same region. Our results suggest that numerous D>200 km planetesimals disrupted early in Solar System history, but only a small fraction of their fragments survived the dynamical depletion event described above. We believe this may explain the limited presence of iron-rich M-type, olivine-rich A-type, and non-Vesta V-type asteroids in the main belt today. The collisional lifetimes determined for main belt asteroids agree with the cosmic ray exposure ages of stony meteorites and are consistent with the limited collisional evolution detected among large Koronis family members. Using the same model, we investigated the near-Earth object (NEO) population. We show the shape of the NEO size distribution is a reflection of the main belt population, with main belt asteroids driven to resonances by Yarkovsky thermal forces. We used our model of the NEO population over the last 3 Gyr, which is consistent with the current population determined by telescopic and satellite data, to explore whether the majority of small craters (D<0.1-1 km) formed on Mercury, the Moon, and Mars were produced by primary impacts or by secondary impacts generated by ejecta from large craters. Our results suggest that most small craters formed on these worlds were a by-product of secondary rather than primary impacts.     

The fossilized size distribution of the main asteroid belt

Planet formation models suggest the primordial main belt experienced a short but intense period of collisional evolution shortly after the formation of planetary embryos. This period is believed to have lasted until Jupiter reached its full size, when dynamical processes (e.g., sweeping resonances, excitation via planetary embryos) ejected most planetesimals from the main belt zone. The few planetesimals left behind continued to undergo comminution at a reduced rate until the present day. We investigated how this scenario affects the main belt size distribution over Solar System history using a collisional evolution model (CoEM) that accounts for these events. CoEM does not explicitly include results from dynamical models, but instead treats the unknown size of the primordial main belt and the nature/timing of its dynamical depletion using innovative but approximate methods. Model constraints were provided by the observed size frequency distribution of the asteroid belt, the observed population of asteroid families, the cratered surface of differentiated Asteroid (4) Vesta, and the relatively constant crater production rate of the Earth and Moon over the last 3 Gyr. Using CoEM, we solved for both the shape of the initial main belt size distribution after accretion and the asteroid disruption scaling law . In contrast to previous efforts, we find our derived function is very similar to results produced by numerical hydrocode simulations of asteroid impacts. Our best fit results suggest the asteroid belt experienced as much comminution over its early history as it has since it reached its low-mass state approximately 3.9-4.5 Ga. These results suggest the main belt's wavy-shaped size-frequency distribution is a “fossil” from this violent early epoch. We find that most diameter D?120 km asteroids are primordial, with their physical properties likely determined during the accretion epoch. Conversely, most smaller asteroids are byproducts of fragmentation events. The observed changes in the asteroid spin rate and lightcurve distributions near D∼100-120 km are likely to be a byproduct of this difference. Estimates based on our results imply the primordial main belt population (in the form of D<1000 km bodies) was 150-250 times larger than it is today, in agreement with recent dynamical simulations.     

Computing the effects of YORP on the spin rate distribution of the NEO population

The overall change of NEO spin rate due to planetary encounters and YORP is evaluated by using a Monte Carlo model. A large sample of test objects mimicking a source population is evolved over a timescale comparable with the Solar System age until they reach a steady state spin distribution that should reproduce the current NEO distribution. The spin change due to YORP is computed for each body according to a simplified model based on Scheeres [Scheeres, D.J., 2007a. Icarus 188, 430-450].The steady state cumulative distribution of NEO spin rates obtained from our simulation nicely reproduces the observed one, once our results are biased to match the diameter distribution of the sample of objects included in the observational database. The excellent agreement strongly suggests that YORP is responsible for the concentration of spin at low rotation rates. In fact, in the absence of YORP the steady state population significantly deviates from the observed one. The spin evolution due to YORP is also so rapid for NEOs that the initial rotation rate distribution of any source population is quickly relaxed to that of the observed population. This has profound consequences for the study of NEO origin since we cannot trace the sources of NEOs from their rotation rate only.     

Minor element evidence that Asteroid 433 Eros is a space-weathered ordinary chondrite parent body

The NEAR mission to 433 Eros provided detailed data on the geology, mineralogy, and chemistry of this S-class asteroid [McCoy, T.J., Robinson, M.S., Nittler, L.R., Burbine, T.H., 2002. Chem. Erde 62, 89-121; Cheng, A.F., 1997. Space Sci. Rev. 82, 3-29] with a key science goal of understanding the relationship between asteroids and meteorites [Cheng, A.F., 1997. Space Sci. Rev. 82, 3-29; Gaffey, M.J., Burbine, T.H., Piatek, J.L., Reed, K.L., Chaky, D.A., Bell, J.F., Brown, R.H., 1993a. Icarus 106, 573-602]. Previously reported major element data revealed a bulk surface similar to that of ordinary chondrites, with the notable exception of sulfur, which was highly depleted [Trombka, J.I., and 23 colleagues, 2000. Science 289, 2101-2105; Nittler, L.R., and 14 colleagues, 2001. Meteorit. Planet. Sci. 36, 1673-1695]. The origin of this sulfur deficiency, and hence the fundamental nature of the asteroid's surface, has remained controversial. We report a new analysis of NEAR X-ray spectrometer data, indicating that Eros has Cr/Fe, Mn/Fe, and Ni/Fe ratios similar to ordinary chondrite meteorites of type LL or L. Chondritic levels of Cr, Mn, and Ni argue strongly against a partial melting explanation for the sulfur depletion. Instead, our results provide definitive evidence that Eros is a primitive body with composition and mineralogy similar to ordinary chondrites, but with a surface heavily modified by interactions with the solar wind and micrometeorites, processes collectively termed space weathering.     

The equilibrium of rubble-pile satellites: The Darwin and Roche ellipsoids for gravitationally held granular aggregates

Many new small moons of the giant planets have been discovered recently. In parallel, satellites of several asteroids, e.g., Ida, have been found. Strikingly, a majority of these new-found planetary moons are estimated to have very low densities, which, along with their hypothesized accretionary origins, suggests a rubble internal structure. This, coupled to the fact that many asteroids are also thought to be particle aggregates held together principally by self-gravity, motivates the present investigation into the possible ellipsoidal shapes that a rubble-pile satellite may achieve as it orbits an aspherical primary. Conversely, knowledge of the shape will constrain the granular aggregate's orbit—the closer it gets to a primary, both primary's tidal effect and the satellite's spin are greater. We will assume that the primary body is sufficiently massive so as not to be influenced by the satellite. However, we will incorporate the primary's possible ellipsoidal shape, e.g., flattening at its poles in the case of a planet, and the proloidal shape of asteroids. In this, the present investigation is an extension of the first classical Darwin problem to granular aggregates. General equations defining an ellipsoidal rubble pile's equilibrium about an ellipsoidal primary are developed. They are then utilized to scrutinize the possible granular nature of small inner moons of the giant planets. It is found that most satellites satisfy constraints necessary to exist as equilibrated granular aggregates. Objects like Naiad, Metis and Adrastea appear to violate these limits, but in doing so, provide clues to their internal density and/or structure. We also recover the Roche limit for a granular satellite of a spherical primary, and employ it to study the martian satellites, Phobos and Deimos, as well as to make contact with earlier work of Davidsson [Davidsson, B., 2001. Icarus 149, 375–383]. The satellite's interior will be modeled as a rigid-plastic, cohesion-less material with a Drucker–Prager yield criterion. This rheology is a reasonable first model for rubble piles. We will employ an approximate volume-averaging procedure that is based on the classical method of moments, and is an extension of the virial method [Chandrasekhar, S., 1969. Ellipsoidal Figures of Equilibrium. Yale Univ. Press, New Haven] to granular solid bodies.     

Simulation and analysis of the dynamics of binary near-Earth Asteroid (66391) 1999 KW4

Near-Earth Asteroid (66391) 1999 KW4 was the subject of the recently published first extensive radar imaging, shape and mutual orbit modeling, and physical and dynamical characterization of a binary asteroid. In this paper we present in detail our numerical simulation of KW4 behind that work. Our propagations of the system with some variation in estimated parameters cover the set of KW4's possible current dynamical states consistent with the body models and other information obtained directly from the observations. We also apply our implementation of this simulation capability to address some of the dynamical mechanisms by which KW4 may be moved into the more energetically excited of those possible current states, particularly solar gravity interaction. Through comparison of the results with certain features of the observation data, we conclude that the actual KW4 system is not in the most energetically relaxed configuration but must be moderately excited. The system occupies a generalized Cassini state 2 which is different from that considered in most previously published treatments of Cassini states in that it involves co-precession of the primary's spin frame and the mutual orbit rather than co-precession of a satellite's spin frame and that satellite's orbit about the primary. We present a simple analytical theory describing the system's dynamics, which should be applicable to any other binary systems, of which KW4 is representative, in which a massive, roughly oblate primary is spinning rapidly relative to the rate of its mutual orbit with an on-average synchronous, elongated secondary. We examine separately both the effect of the larger binary component's oblateness, and the effect of the smaller component's roughly triaxial ellipsoid shape. The simple analytical formulae obtained agree with full-detail numerical simulation results, and can be used for remote estimation of binary mass properties from observed system motion.     

2002 AA29: Earth's recurrent quasi-satellite?

We study the dynamical evolution of Asteroid 2002 AA₂₉. This object moves in the co-orbital region of the Earth and is the first known asteroid which experiences recurrent horseshoe-quasi-satellite transitions. The transitions between the HS and QS states are unique among other known Earth co-orbital asteroids and in the QS state 2002 AA₂₉ remains very close to Earth (within 0.2 AU for several decades [Connors, M., Chodas, P., Mikkola, S., Wiegert, P., Veillet, C., Innanen, K., 2002. Meteorit. Planet. Sci. 37, 1435-1441]). Based on results obtained analytically by Brasser et al. [Brasser, R., Heggie, D.C., Mikkola, S., 2004b. Celest. Mech. Dynam. Astron. 88, 123-152] we developed a simple analytical method to describe and analyze the motion of 2002 AA₂₉. We distinguish a few moments in time crucial for understanding its dynamics. Near 2400 and 2500 this object will be close to going through the maxima of the averaged disturbing function and it will either change its co-orbital regime by transition from the HS into QS state, or leave the librating mode. These approaches generate instability in the motion of 2002 AA₂₉. By means of 66 observations, covering a two-year interval, we extend the analysis of the long term evolution of this object presented by Connors et al. [Connors, M., Chodas, P., Mikkola, S., Wiegert, P., Veillet, C., Innanen, K., 2002. Meteorit. Planet. Sci. 37, 1435-1441] and Brasser et al. [Brasser, R., Innanen, K.A., Connors, M., Veillet, C., Wiegert, P., Mikkola, S., Chodas, P.W., 2004a. Icarus 171, 102-109]. Our analysis is based on a sample of 100 cloned orbits. We show that the motion of 2002 AA₂₉ is predictable in the time interval [−2600,7100] and outside of this interval the past and future orbital history can be studied using statistical methods.     

On the shapes and spins of “rubble pile” asteroids

We examine the shape of a “rubble pile” asteroid as it slowly gains angular momentum by YORP torque, to the point where “landsliding” occurs. We find that it evolves to a “top” shape with constant angle of repose from the equator up to mid-latitude, closely resembling the shapes of several nearly critically spinning asteroids imaged by radar, most notably (66391) 1999 KW4 [Ostro, S.J., Margot, J.-L., Benner, L.A.M., Giorgini, J.D., Scheeres, D.J., Fahnestock, E.G., Broschart, S.B., Bellerose, J., Nolan, M.C., Magri, C., Pravec, P., Scheirich, P., Rose, R., Jurgens, R.F., De Jong, E.M., Suzuki, S., 2006. Science 314, 1276-1280]. Similar calculations for non-spinning extremely prolate or oblate “rubble piles” show that even loose rubble can sustain shapes far from fluid equilibrium, thus inferences based on fluid equilibrium are generally useless for inferring bulk properties such as density of small bodies. We also investigate the tidal effects of a binary system with a “top shape” primary spinning at near the critical limit for stability. We find that very close to the stability limit, the tide from the secondary can actually levitate loose debris from the surface and re-deposit it, in a process we call “tidal saltation.” In the process, angular momentum is transferred from the primary spin to the satellite orbit, thus maintaining the equilibrium of near-critical spin as YORP continues to add angular momentum to the system. We note that this process is in fact dynamically related to the process of “shepherding” of narrow rings by neighboring satellites.     

The composition of M-type asteroids: Synthesis of spectroscopic and radar observations

We have conducted a radar-driven observational campaign of 22 main-belt asteroids (MBAs) focused on Bus–DeMeo Xc- and Xk-type objects (Tholen X and M class asteroids) using the Arecibo radar and NASA Infrared Telescope Facilities (IRTF). Sixteen of our targets were near-simultaneously observed with radar and those observations are described in a companion paper (Shepard, M.K., and 19 colleagues [2010]. Icarus, in press). We find that most of the highest metal-content asteroids, as suggested by radar, tend to exhibit silicate absorption features at both 0.9 and 1.9 μm, and the lowest metal-content asteroids tend to exhibit either no bands or only the 0.9 μm band. Eleven of the asteroids were observed at several rotational longitudes in the near-infrared and significant variations in continuum slope were found for nine in the spectral regions 1.1–1.45 μm and 1.6–2.3 μm. We utilized visible wavelength data (Bus, S.J., Binzel, R.P. [2002b]. Icarus 158, 146–177; Fornasier, S., Clark, B.E., Dotto, E., Migliorini, A., Ockert-Bell, M., Barucci, M.A. [2010]. Icarus 210, 655–673.) for a more complete compositional analysis of our targets. Compositional evidence is derived from our target asteroid spectra using two different methods: (1) a χ² search for spectral matches in the RELAB database, and (2) parametric comparisons with meteorites. This paper synthesizes the results of the RELAB search and the parametric comparisons with compositional suggestions based on radar observations. We find that for six of the seven asteroids with the highest iron abundances, our spectral results are consistent with the radar evidence (16 Psyche, 216 Kleopatra, 347 Pariana, 758 Mancunia, 779 Nina, and 785 Zwetana). Three of the seven asteroids with the lowest metal abundances, our spectral results are consistent with the radar evidence (21 Lutetia, 135 Hertha, 497 Iva). The remaining seven asteroids (22 Kalliope, 97 Klotho, 110 Lydia, 129 Antigone, 224 Oceana, 678 Fredegundis, and 771 Libera) have ambiguous compositional interpretations when comparing the spectral analogs to the radar analogs. The number of objects with ambiguous results from this multi-wavelength survey using visible, near-infrared, and radar wavelengths indicates that perhaps a third diagnostic wavelength region (such as the mid-infrared around 2–4 μm, the mid-infrared around 10–15 μm, and/or the ultraviolet around 0.2–0.4 μm) should be explored to resolve the discrepancies.     

V-type asteroids in the middle main belt

V-type asteroids are bodies whose surfaces are constituted of basalt. In the Main Asteroid Belt, most of these asteroids are assumed to come from the basaltic crust of Asteroid (4) Vesta. This idea is mainly supported by (i) the fact that almost all the known V-type asteroids are in the same region of the belt as (4) Vesta, i.e., the inner belt (semi-major axis 2.1<a<2.5 AU), (ii) the existence of a dynamical asteroid family associated to (4) Vesta, and (iii) the observational evidence of at least one large craterization event on Vesta's surface. One V-type asteroid that is difficult to fit in this scenario is (1459) Magnya, located in the outer asteroid belt, i.e., too far away from (4) Vesta as to have a real possibility of coming from it. The recent discovery of the first V-type asteroid in the middle belt (2.5<a<2.8 AU), (21238) 1995WV7 [Binzel, R.P., Masi, G., Foglia, S., 2006. Bull. Am. Astron. Soc. 38, 627; Hammergren, M., Gyuk, G., Puckett, A., 2006. ArXiv e-print, astro-ph/0609420], located at ∼2.54 AU, raises the question of whether it came from (4) Vesta or not. In this paper, we present spectroscopic observations indicating the existence of another V-type asteroid at ∼2.53 AU, (40521) 1999RL95, and we investigate the possibility that these two asteroids evolved from the Vesta family to their present orbits by a semi-major axis drift due to the Yarkovsky effect. The main problem with this scenario is that the asteroids need to cross the 3/1 mean motion resonance with Jupiter, which is highly unstable. Combining N-body numerical simulations of the orbital evolution, that include the Yarkovsky effect, with Monte Carlo models, we compute the probability that an asteroid of a given diameter D evolves from the Vesta family and crosses over the 3/1 resonance, reaching a stable orbit in the middle belt. Our results indicate that an asteroid like (21238) 1995WV7 has a low probability (∼1%) of having evolved through this mechanism due to its large size (D∼5 km), because the Yarkovsky effect is not sufficiently efficient for such large asteroids. However, the mechanism might explain the orbits of smaller bodies like (40521) 1999RL95 (D∼3 km) with ∼70-100% probability, provided that we assume that the Vesta family formed ?3.5 Gy ago. We estimate the debiased population of V-type asteroids that might exist in the same region as (21238) and (40521) (2.5<a?2.62 AU) and conclude that about 10 to 30% of the V-type bodies with D>1 km may come from the Vesta family by crossing over the 3/1 resonance. The remaining 70-90% must have a different origin.