The large crater on the small Asteroid (2867) Steins

The maximum size of impact craters on finite bodies marks the largest impact that can occur short of impact induced disruption of the body. Recently attention has started to focus on large craters on small bodies such as asteroids and rocky and icy satellites. Here the large crater on the recently imaged Asteroid (2867) Steins (with crater diameter to mean asteroid radius ratio of 0.79) is shown to follow a limit set by other similar sized bodies with moderate macroporosity (i.e. fractured asteroids). Thus whilst large, the crater size is not novel, nor does it require Steins to possess an extremely large porosity. In one of the components of the binary Asteroid (90) Antiope there is the recently reported presence of an extremely large depression, possibly a crater, with depression diameter to mean asteroid radius ratio of ∼(1.4–1.62). This is consistent with the maximum size of a crater expected from previous observations of very porous rocky bodies (i.e. rubble-pile asteroids). Finally, a relationship between crater diameter (normalised to body radius) is proposed as a function of body porosity which suggests that the doubling of porosity between fractured asteroids and rubble-pile asteroids, nearly doubles the size (D/R value) of the largest crater sustainable on a rocky body.     

Radar observations of E-class Asteroids 44 Nysa and 434 Hungaria

We observed the E-class main-belt Asteroids (MBAs) 44 Nysa and 434 Hungaria with Arecibo Observatory's S-band (12.6 cm) radar. Both asteroids exhibit polarization ratios higher than those measured for any other MBA: Nysa, μc=0.50±0.02 and Hungaria, μc=0.8±0.1. This is consistent with the high polarization ratios measured for every E-class near-Earth asteroid (NEA) observed by Benner et al. [Benner, L.A.M., and 10 collegues, 2008. Icarus, submitted for publication] and suggests a common cause. Our estimates of radar albedo are 0.19±0.06 for Nysa and 0.22±0.06 for Hungaria. These values are higher than those of most MBAs and, when combined with their high polarization ratios, suggest that the surface bulk density of both asteroids is high. We model Nysa as an ellipsoid of dimension 113×67×65 km (±15%) giving an effective diameter Deff=79±10 km, consistent with previous estimates. The echo waveforms are not consistent with a contact binary as suggested by Kaasalainen et al. [Kaasalainen, M., Torppa, J., Piironen, J., 2002. Astron. Astrophys. 383, L19-L22]. We place a constraint on Hungaria's maximum diameter, Dmax?11 km consistent with previous size estimates.     

Extreme sensitivity of the YORP effect to small-scale topography

Radiation recoil (YORP) torques are shown to be extremely sensitive to small-scale surface topography, using numerical simulations. Starting from a set of “base objects” representative of the near-Earth object population, random realizations of three types of small-scale topography are added: Gaussian surface fluctuations, craters, and boulders. For each, the expected relative errors in the spin and obliquity components of the YORP torque caused by the observationally unresolved small-scale topography are computed. Gaussian power, at angular scales below an observational limit, produces expected errors of order 100% if observations constrain the surface to a spherical harmonic order l?10. For errors under 10%, the surface must be constrained to at least l=20. A single crater with diameter roughly half the object's mean radius, placed at random locations, results in expected errors of several tens of percent. The errors scale with crater diameter D as D² for D>0.3 and as D³ for D<0.3 mean radii. Objects that are identical except for the location of a single large crater can differ by factors of several in YORP torque, while being photometrically indistinguishable at the level of hundredths of a magnitude. Boulders placed randomly on identical base objects create torque errors roughly 3 times larger than do craters of the same diameter, with errors scaling as the square of the boulder diameter. A single boulder comparable to Yoshinodai on 25143 Itokawa, moved by as little as twice its own diameter, can alter the magnitude of the torque by factors of several, and change the sign of its spin component at all obliquities. Most of the total torque error produced by multiple unresolved craters is contributed by the handful of largest craters; but both large and small boulders contribute comparably to the total boulder-induced error. A YORP torque prediction derived from groundbased data can be expected to be in error by of order 100% due to unresolved topography. Small surface changes caused by slow spin-up or spin-down may have significant stochastic effects on the spin evolution of small bodies. For rotation periods between roughly 2 and 10 h, these unpredictable changes may reverse the sign of the YORP torque. Objects in this spin regime may random-walk up and down in spin rate before the rubble-pile limit is exceeded and fissioning or loss of surface objects occurs. Similar behavior may be expected at rotation rates approaching the limiting values for tensile-strength dominated objects.     

The formation of asteroid satellites in large impacts: results from numerical simulations

We present results of 161 numerical simulations of impacts into 100-km diameter asteroids, examining debris trajectories to search for the formation of bound satellite systems. Our simulations utilize a 3-dimensional smooth-particle hydrodynamics (SPH) code to model the impact between the colliding asteroids. The outcomes of the SPH models are handed off as the initial conditions for N-body simulations, which follow the trajectories of the ejecta fragments to search for the formation of satellite systems. Our results show that catastrophic and large-scale cratering collisions create numerous fragments whose trajectories can be changed by particle-particle interactions and by the reaccretion of material onto the remaining target body. Some impact debris can enter into orbit around the remaining target body, which is a gravitationally reaccreted rubble pile, to form a SMAshed Target Satellite (SMATS). Numerous smaller fragments escaping the largest remnant may have similar trajectories such that many become bound to one another, forming Escaping Ejecta Binaries (EEBs). Our simulations so far seem to be able to produce satellite systems qualitatively similar to observed systems in the main asteroid belt. We find that impacts of 34-km diameter projectiles striking at 3 km s⁻¹ at impact angles of ∼30° appear to be particularly efficient at producing relatively large satellites around the largest remnant as well as large numbers of modest-size binaries among their escaping ejecta.     

The topographic limits of gravitationally bound, rotating sand piles

Rubble pile asteroids can attain shapes that are dramatically different from those of rotating, self-gravitating equilibrium fluids. A new numerical technique, called “seed growth,” is demonstrated for calculating three-dimensional bodies that are self-gravitating and rotating, and whose every surface is approximately at a constant angle, ?, with respect to the local horizontal. By altering the configuration of cusps, which are points along a constant longitude path where the surface angle changes sign but not magnitude, multiple solution shapes that satisfy the condition that all surface slopes are at a constant angle are possible. Five different cusp configurations are explored here, three of which yield solutions for 20°???30°. Rotational effects are explored, and it is found that for some solution shapes, the ratios of their shortest to longest dimensions, c/a, can fall outside the limits published in the literature for rotating, cohesionless, spheroidal bodies. Solution shapes show some similarities to observed small bodies, such as the saturnian satellite Atlas, the near-Earth Asteroid 1999 KW4, and some contact binary asteroids.     

The formation of asteroid satellites in large impacts: results from numerical simulations

We present results of 161 numerical simulations of impacts into 100-km diameter asteroids, examining debris trajectories to search for the formation of bound satellite systems. Our simulations utilize a 3-dimensional smooth-particle hydrodynamics (SPH) code to model the impact between the colliding asteroids. The outcomes of the SPH models are handed off as the initial conditions for N-body simulations, which follow the trajectories of the ejecta fragments to search for the formation of satellite systems. Our results show that catastrophic and large-scale cratering collisions create numerous fragments whose trajectories can be changed by particle-particle interactions and by the reaccretion of material onto the remaining target body. Some impact debris can enter into orbit around the remaining target body, which is a gravitationally reaccreted rubble pile, to form a SMAshed Target Satellite (SMATS). Numerous smaller fragments escaping the largest remnant may have similar trajectories such that many become bound to one another, forming Escaping Ejecta Binaries (EEBs). Our simulations so far seem to be able to produce satellite systems qualitatively similar to observed systems in the main asteroid belt. We find that impacts of 34-km diameter projectiles striking at 3 km s⁻¹ at impact angles of ∼30° appear to be particularly efficient at producing relatively large satellites around the largest remnant as well as large numbers of modest-size binaries among their escaping ejecta.     

Shape, size and multiplicity of main-belt asteroids: I. Keck Adaptive Optics survey

This paper presents results from a high spatial resolution survey of 33 main-belt asteroids with diameters >40 km using the Keck II Adaptive Optics (AO) facility. Five of these (45 Eugenia, 87 Sylvia, 107 Camilla, 121 Hermione, 130 Elektra) were confirmed to have satellite. Assuming the same albedo as the primary, these moonlets are relatively small (∼5% of the primary size) suggesting that they are fragments captured after a disruptive collision of a parent body or captured ejecta due to an impact. For each asteroid, we have estimated the minimum size of a moonlet that can positively detected within the Hill sphere of the system by estimating and modeling a 2-σ detection profile: in average on the data set, a moonlet located at 2/100×RHill (1/4×RHill) with a diameter larger than 6 km (4 km) would have been unambiguously seen. The apparent size and shape of each asteroid was estimated after deconvolution using a new algorithm called AIDA. The mean diameter for the majority of asteroids is in good agreement with IRAS radiometric measurements, though for asteroids with a D<200 km, it is underestimated on average by 6-8%. Most asteroids had a size ratio that was very close to those determined by lightcurve measurements. One observation of 104 Klymene suggests it has a bifurcated shape. The bi-lobed shape of 121 Hermione described in Marchis et al. [Marchis, F., Hestroffer, D., Descamps, P., Berthier, J., Laver, C., de Pater, I., 2005c. Icarus 178, 450-464] was confirmed after deconvolution. The ratio of contact binaries in our survey, which is limited to asteroids larger than 40 km, is surprisingly high (∼6%), suggesting that a non-single configuration is common in the main-belt. Several asteroids have been analyzed with lightcurve inversions. We compared lightcurve inversion models for plane-of-sky predictions with the observed images (9 Metis, 52 Europa, 87 Sylvia, 130 Elektra, 192 Nausikaa, and 423 Diotima, 511 Davida). The AO images allowed us to determine a unique photometric mirror pole solution, which is normally ambiguous for asteroids moving close to the plane of the ecliptic (e.g., 192 Nausikaa and 52 Europa). The photometric inversion models agree well with the AO images, thus confirming the validity of both the lightcurve inversion method and the AO image reduction technique.     

Impact cratering records of the mid-sized, icy saturnian satellites

Resolution of Voyager 1 and 2 images of the mid-sized, icy saturnian satellites was generally not much better than 1 km per line pair, except for a few, isolated higher resolution images. Therefore, analyses of impact crater distributions were generally limited to diameters (D) of tens of kilometers. Even with the limitation, however, these analyses demonstrated that studying impact crater distributions could expand understanding of the geology of the saturnian satellites and impact cratering in the outer Solar System. Thus to gain further insight into Saturn’s mid-sized satellites and impact cratering in the outer Solar System, we have compiled cratering records of these satellites using higher resolution CassiniISS images. Images from Cassini of the satellites range in resolution from tens m/pixel to hundreds m/pixel. These high-resolution images provide a look at the impact cratering records of these satellites never seen before, expanding the observable craters down to diameters of hundreds of meters. The diameters and locations of all observable craters are recorded for regions of Mimas, Tethys, Dione, Rhea, Iapetus, and Phoebe. These impact crater data are then analyzed and compared using cumulative, differential and relative (R) size-frequency distributions. Results indicate that the heavily cratered terrains on Rhea and Iapetus have similar distributions implying one common impactor population bombarded these two satellites. The distributions for Mimas and Dione, however, are different from Rhea and Iapetus, but are similar to one another, possibly implying another impactor population common to those two satellites. The difference between these two populations is a relative increase of craters with diameters between 10 and 30 km and a relative deficiency of craters with diameters between 30 and 80 km for Mimas and Dione compared with Rhea and Iapetus. This may support the result from Voyager images of two distinct impactor populations. One population was suggested to have a greater number of large impactors, most likely heliocentric comets (Saturn Population I in the Voyager literature), and the other a relative deficiency of large impactors and a greater number of small impactors, most likely planetocentric debris (Saturn Population II). Meanwhile, Tethys’ impact crater size-frequency distribution, which has some similarity to the distributions of Mimas, Dione, Rhea, and Iapetus, may be transitional between the two populations. Furthermore, when the impact crater distributions from these older cratered terrains are compared to younger ones like Dione’s smooth plains, the distributions have some similarities and differences. Therefore, it is uncertain whether the size-frequency distribution of the impactor population(s) changed over time. Finally, we find that Phoebe has a unique impact crater distribution. Phoebe appears to be lacking craters in a narrow diameter range around 1 km. The explanation for this confined “dip” at D = 1 km is not yet clear, but may have something to do with the interaction of Saturn’s irregular satellites or the capture of Phoebe.     

Tidal disruptions: A continuum theory for solid bodies

Although the theory of Roche 1847 for the tidal disruption limits of orbiting satellites assumes a fluid body, a length to diameter of exactly 2.07:1, and a particular body orientation, the theory is commonly applied to the satellites of the Solar System and to small asteroids and comets passing nearby a planet. Clearly these bodies are neither fluid nor generally are that elongated, so a more appropriate theory is needed. Here we present exact analytical results for the distortion and disruption limits of solid spinning ellipsoidal bodies subjected to tidal forces, using the Drucker-Prager strength model with zero cohesion. It is the appropriate model for dry granular materials such as sands and rocks, for rubble-pile asteroids and comets, and for larger satellites, asteroids and comets where the cohesion can be ignored. This study uses the same approach as the studies of spin limits for solid ellipsoidal bodies given in [Holsapple, K.A., 2001. Icarus 154, 432-448; Holsapple, K.A., 2004. Icarus 172, 272-303]. It is a static theory that predicts conditions for breakup and predicts the nature of the deformations at the limit state, but does not track the dynamics of the body as it comes apart. The strength is characterized by a single material parameter associated with an angle of friction, which can range from zero to 90°. The case with zero friction angle has no shear strength whatsoever, so it is then the model of a fluid or gas. The case of 90° represents a material that cannot fail in shear, but still has zero tensile strength. Typical dry soils have angles of friction of 30°-40°. Since the static fluid case is included in the theory as a special case, the classical results of Roche [Roche, E.A., 1847. Acad. Sci. Lett. Montpelier. Mem. Section Sci. 1, 243-262] and Jeans [Jeans, J.H., 1917. Mem. R. Astron. Soc. London 62, 1-48] are included and re-derived in their entirety; but the general solid case has much more variety and applicability. We consider both the spin-locked case, appropriate for most satellites of the Solar System; and the zero spin case, a possible case for a passing stray body. Detailed plots of many special cases are presented, in terms of shape, orientation and mass densities. A very typical result gives a closest approach d=1.5(ρ/ρP)^1/3R in terms of the planet radius R, and the satellite and planet mass densities ρ and ρP. We also use the theory to distinguish between conditions allowing global shape changes leading to new equilibrium states, or those leading to complete disruption. We apply the theory to the potentially hazardous Asteroid 99942 Apophis due to pass very near the Earth in 2029, and conclude it is extremely unlikely to experience any tidal readjustments during its passage. The states of many of the satellites of the Solar System are compared to the theory, and we find that all are well within their tidal disruption limits for expected values of the internal friction.     

Critical crater diameter and asteroid impact seismology

Abstract— If impact stress reverberation is the primary gradational process on an asteroid at global scales, then the largest undegraded crater records an asteroid's seismological response. The critical crater diameter D_crit is defined as the smallest crater whose formation disrupts all previous craters globally up to its size; it is solved for by combining relationships for crater growth and for stress attenuation. The computation for D_crit gives a simple explanation for the curious observation that small asteroids have only modest undegraded craters, in comparison to their size, whereas large asteroids have giant undegraded craters. D_crit can even exceed the asteroid diameter, in which case all craters are “local” and the asteroid becomes crowded with giant craters. D_crit is the most recent crater to have formed on a blank slate; when it is equated to the measured diameter of the largest undegraded crater on known asteroids, peak particle velocities are found to attenuate with the 1.2–1.3 power of distance—less attenuative than strong shocks, and more characteristic of powerful seismic disturbances. This is to be expected, since global degradation can result from seismic (cm s^?1) particle velocities on small asteroids. Attenuation, as modeled, appears to be higher on asteroids known to be porous, although these are also bodies for which different crater scaling rules might apply.     

Impact crater formed on sintered snow surface simulating porous icy bodies

To improve the scaling parameter controlling the impact crater formation in the strength regime, we conducted impact experiments on sintered snow targets with the dynamic strength continuously changed from 20 to 200 kPa, and the largest crater size formed on small icy satellites was considered by using the revised scaling parameter. Ice and snow projectiles were impacted on a snow surface with 36% porosity at an impact velocity from 31 m s⁻¹ to 150 m s⁻¹. The snow target was sintered at the temperature from −5 °C to −18 °C, and the snow dynamic strength was changed with the sintering duration at each temperature. We found that the mass ejected from the crater normalized by the projectile mass, π_V, was related to the ratio of the dynamic strength to the impact pressure, , as follows: , where the impact pressure was indicated by P = ρ_tC_0tv_i/2 with the target density of ρ_t, when the impact velocity, v_i, was much smaller than the bulk sound velocity C_0t (typically 1.8 km s⁻¹ in our targets). The ratio of the largest crater diameter to the diameter of the target body, d_max/D, was estimated by calculating the crater diameter at the impact condition for catastrophic disruption and then compared to the observed d_max/D of jovian and saturnian small satellites, in order to discuss the formation condition of these large d_max/D in the strength regime.     

Sheath-limited unipolar induction in the solar wind

The collection of charged particles by electrodes in plasmas is controlled by the currentvoltage characteristics of the plasma sheath which forms at the electrode surface. This principle is applied to the steady-state electromagnetic interaction of the solar wind with moon-like bodies, or ‘solid-body’ interactions. In some cases the unipolar dynamo response of an electrically conducting body in the solar wind motional electric field can be controlled by sheath effects. This occurs for highly conducting bodies when the body radiusR is less than a critical valueR _c , with the result that no induced bow shock wave can form. For the Moon and MercuryR/R _c ?1, so that sheath effects do not limit their unipolar responses. The asteroids are found to be either too cold or too small to maintain steady-state induced magnetospheres. The Martian satellites, the irregular Jovian satellites, and the outer satellites of Saturn also haveR/R _c ?1. No bow shock waves should be generated by these bodies, unless they are highly magnetized or have large magnetic permeabilities. Unipolar induction heating of meteorite parent bodies in a primordial enhanced solar wind should not be inhibited by sheath effects, providedR?50 m.     

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.     

Orbital evolution of the Gefion and Adeona asteroid families: close encounters with massive asteroids and the Yarkovsky effect

Asteroid families are groupings of minor planets identified by clustering in their proper orbital elements; these objects have spectral signatures consistent with an origin in the break-up of a common parent body. From the current values of proper semimajor axes a of family members one might hope to estimate the ejection velocities with which the fragments left the putative break-up event (assuming that the pieces were ejected isotropically). However, the ejection velocities so inferred are consistently higher than N-body and hydro-code simulations, as well as laboratory experiments, suggest. To explain this discrepancy between today’s orbital distribution of asteroid family members and their supposed launch velocities, we study whether asteroid family members might have been ejected from the collision at low speeds and then slowly drifted to their current positions, via one or more dynamical processes. Studies show that the proper a of asteroid family members can be altered by two mechanisms: (i) close encounters with massive asteroids, and (ii) the Yarkovsky non-gravitational effect. Because the Yarkovsky effect for kilometer-sized bodies decreases with asteroid diameter D, it is unlikely to have appreciably moved large asteroids (say those with D > 15 km) over the typical family age (1-2 Gyr).For this reason, we numerically studied the mobility of family members produced by close encounters with main-belt, non-family asteroids that were thought massive enough to significantly change their orbits over long timescales. Our goal was to learn the degree to which perturbations might modify the proper a values of all family members, including those too large to be influenced by the Yarkovsky effect. Our initial simulations demonstrated immediately that very few asteroids were massive enough to significantly alter relative orbits among family members. Thus, to maximize gravitational perturbations in our 500-Myr integrations, we investigated the effect of close encounters on two families, Gefion and Adeona, that have high encounter probabilities with 1 Ceres, by far the largest asteroid in the main belt. Our results show that members of these families spreads in a of less than 5% since their formation. Thus gravitational interactions cannot account for the large inferred escape velocities.The effect of close encounters with massive asteroids is, however, not entirely negligible. For about 10% of the simulated bodies, close encounters increased the “inferred” ejection velocities from sub-100 m/s to values greater than 100 m/s, beyond what hydro-code and N-body simulations suggest are the maximum possible initial ejection velocity for members of Adeona and Gefion with D > 15 km. Thus this mechanism of mobility may be responsible for the unusually high inferred ejection speeds of a few of the largest members of these two families.To understand the orbital evolution of the entire family, including smaller members, we also performed simulations to account for the drift of smaller asteroids caused by the Yarkovsky effect. Our two sets of simulations suggest that the two families we investigated are relatively young compared to larger families like Koronis and Themis, which have estimated ages of about 2 Byr. The Adeona and Gefion families seems to be no more than 600 and 850 Myr old, respectively.     

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.     

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.     

A possible YORP effect on C and S Main Belt Asteroids

A rotating frequency analysis in a previous paper, showed that two samples of C and S-type asteroids belonging to the Main Belt, but not to any families, present two different values for the transition diameter to a Maxwellian distribution of the rotation frequency, respectively 48 and 33 km. In this paper, after a more detailed statistical analysis, aiming to verify that the result is physically relevant, we found a better estimate for the transition diameter, respectively D_C = 44 ± 2 km and D_S = 30 ± 1 km. The ratio between these estimated transition diameters, D_C/D_S = 1.5 ± 0.1, can be supported with the help of the YORP (Yarkovsky-O’Keefe-Radzievskii-Paddack) effect, although other physical causes cannot be completely ruled out.In this paper we have derived a simple scaling law for YORP which, taking into account the different average heliocentric distance, the bulk density, the albedo and the asteroid “asymmetry surface factor”, has enabled us to reasonably justify the ratio between the diameters transition of C-type and S-type asteroids. The same scaling law can be used to estimate a new ratio between the bulk densities of S and C asteroids samples (giving ρ_S/ρ_C ≈ 2.9 ± 0.3), and can explain why the asteroids near the transition diameter have about the same absolute magnitude. For C-type asteroids, using the found density ratio and other estimates of S-type density, it is also possible to estimate an average bulk density equal to 0.9 ± 0.1 g cm⁻³, a value compatible with icy composition. The suggested explanation for the difference of the transition diameters is a plausible hypothesis, consistent with the data, but it needs to be studied more in depth with further observations.     

Photometric survey of binary near-Earth asteroids

Photometric data on 17 binary near-Earth asteroids (15 of them are certain detections, two are probables) were analysed and characteristic properties of the near-Earth asteroid (NEA) binary population were inferred. We have found that binary systems with a secondary-to-primary mean diameter ratio Ds/Dp?0.18 concentrate among NEAs smaller than 2 km in diameter; the abundance of such binaries decreases significantly among larger NEAs. Secondaries show an upper size limit of Ds=0.5-1 km. Systems with Ds/Dp?0.5 are abundant but larger satellites are significantly less common. Primaries have spheroidal shapes and they rotate rapidly, with periods concentrating between 2.2 to 2.8 h and with a tail of the distribution up to ∼4 h. The fast rotators are close to the critical spin for rubble piles with bulk densities about 2 g/cm³. Orbital periods show an apparent cut-off at Porb∼11 h; closer systems with shorter orbital periods have not been discovered, which is consistent with the Roche limit for strengthless bodies. Secondaries are more elongated on average than primaries. Most, but not all, of their rotations appear to be synchronized with the orbital motion; nonsynchronous secondary rotations may occur especially among wider systems with Porb>20 h. The specific total angular momentum of most of the binary systems is similar to within ±20% and close to the angular momentum of a sphere with the same total mass and density, rotating at the disruption limit; this suggests that the binaries were created by mechanism(s) related to rotation near the critical limit and that they neither gained nor lost significant amounts of angular momentum during or since formation. A comparison with six small asynchronous binaries detected in the main belt of asteroids suggests that the population extends beyond the region of terrestrial planets, but with characteristics shifted to larger sizes and longer periods. The estimated mean proportion of binaries with Ds/Dp?0.18 among NEAs larger than 0.3 km is 15±4%. Among fastest rotating NEAs larger than 0.3 km with periods between 2.2 and 2.8 h, the mean proportion of such binaries is (66⁺¹⁰₋₁₂)%.     

Spin limits of Solar System bodies: From the small fast-rotators to 2003 EL61

Holsapple [Holsapple, K.A., 2001. Icarus 154, 432-448; Holsapple, K.A., 2004. Icarus 172, 272-303] determined the spin limits of bodies using a model for solid bodies without tensile or cohesive strength, but with the pressure-induced shear strengths characteristic of dry sands and gravels. That theory included the classical analyses for fluid bodies given by Maclaurin, Jacobi and others as a special case. For the general solid bodies, it was shown that there exists a very wide range of permissible shapes and spin limits; and explicit algebraic results for those limits were given. This paper gives an extension of those analyses to include geological-like materials that also have tensile and cohesive strength. Those strengths are necessary to explain the smaller, fast-rotating asteroids discovered in the last few years. I find that the spin limits for these more general solids have two limiting regimes: a strength regime for bodies with a diameter <3 km, and a gravity regime for the larger bodies with a diameter >10 km (which is the case covered by the earlier papers). I derive explicit algebraic forms for the dependence of the spin limits on shape, mass density and material strength properties. The comparison of the theory to the database for the spins of asteroids and trans-neptunian objects (TNO's) objects shows excellent agreement. For large bodies (diameter D>10 km), the presence of cohesive and/or tensile strength does not permit higher spin rates than would be allowed for rubble pile bodies. Thus, the fact that the spin rates of all large bodies is limited to periods greater than about 2 h does not imply that they are rubble piles. In contrast, for small bodies (D<10 km) the presence of even a very small amount of strength allows much more rapid spins. Small bodies might then be rubble piles but require a small amount of bonding. Finally, I make some remarks about the application of the theory to the TNO's and large asteroids, and question whether a common assumption by researchers that those bodies must take on relaxed fluid shapes is warranted. If not, then the densities and shapes required by that assumption are not valid. I use 2003 EL61 as a prime example.     

Multi-wavelength observations of Asteroid 2100 Ra-Shalom

We observed near-Earth asteroid (NEA) 2100 Ra-Shalom over a six-year period, obtaining rotationally resolved spectra in the visible, near-infrared, thermal-infrared, and radar wavelengths. We find that Ra-Shalom has an effective diameter of Deff=2.3±0.2 km, rotation period P=19.793±0.001 h, visual albedo pv=0.13±0.03, radar albedo , and polarization ratio μc=0.25±0.04. We used our radar observations to generate a three-dimensional shape model which shows several structural features of interest. Based on our thermal observations, Ra-Shalom has a high thermal inertia of ∼10³ J m⁻² s^−0.5 K⁻¹, consistent with a coarse or rocky surface and the inferences of others [Harris, A.W., Davies, J.K., Green, S.F., 1998. Icarus 135, 441-450; Delbo, M., Harris, A.W., Binzel, R.P., Pravec, P., Davies, J.K., 2003. Icarus 166, 116-130]. Our spectral data indicate that Ra-Shalom is a K-class asteroid and we find excellent agreement between our spectra and laboratory spectra of the CV3 meteorite Grosnaja. Our spectra show rotation-dependent variations consistent with global variations in grain size. Our radar observations show rotation-dependent variations in radar albedo consistent with global variations in the thickness of a relatively thin regolith.