Scaling forces to asteroid surfaces: The role of cohesion

The scaling of physical forces to the extremely low ambient gravitational acceleration regimes found on the surfaces of small asteroids is performed. Resulting from this, it is found that van der Waals cohesive forces between regolith grains on asteroid surfaces should be a dominant force and compete with particle weights and be greater, in general, than electrostatic and solar radiation pressure forces. Based on this scaling, we interpret previous experiments performed on cohesive powders in the terrestrial environment as being relevant for the understanding of processes on asteroid surfaces. The implications of these terrestrial experiments for interpreting observations of asteroid surfaces and macro-porosity are considered, and yield interpretations that differ from previously assumed processes for these environments. Based on this understanding, we propose a new model for the end state of small, rapidly rotating asteroids which allows them to be comprised of relatively fine regolith grains held together by van der Waals cohesive forces.     

Motion of a Meteoroid Released from an Asteroid

Evidence of asteroid surface features as regolith grains and larger boulders implies resurfacing possibility due to external forces such as gravitational tidal force during close planet encounters. Motion of a meteoroid released from an asteroid in the gravitational fields of the asteroid and the Earth is modeled. We are interested mainly in a distance between the meteoroid and the asteroid as a function of the time. Applications to Itokawa and some close approaching NEAs are presented.     

Catastrophic fragmentation and formation of families: Preliminary results from a new numerical model

Preliminary results of an improved version of the semiempirical model for catastrophic break up processes developed by Paolicchi et al., (1989) are presented. Among the several changes with respect to the old version, the most important seem to be related to the new treatment of gravitational effects, including self-compression and reaccumulation of fragments. In particular, the new model is able to analyze processes involving both cm-sized objects, like those studied by means of laboratory experiments, as well as much larger bodies, for which self-gravitational effects are dominant; moreover, in this latter case the model seems in principle adequate to describe with the same physics very different phenomena, like the formation of plausible asteroid families and the creation of single, rapidly spinning, objects. This fact, if confirmed by refined analyses, may be of high importance for our general understanding of asteroid collisional evolution.     

Catastrophic disruption of pre-shattered parent bodies

In this paper, we analyze the effect of the internal structure of a parent body on its fragment properties following its disruption in different impact energy regimes. To simulate an asteroid breakup, we use the same numerical procedure as in our previous studies, i.e., a 3D SPH hydrocode to compute the fragmentation phase and the parallel N-body code pkdgrav to compute the subsequent gravitational re-accumulation phase. To explore the importance of the internal structure in determining the collisional outcome, we consider two different parent body models: (1) a purely monolithic one and (2) a pre-shattered one which consists of several fragments separated by damaged zones and small voids. We present here simulations spanning two different impact energy regimes—barely disruptive and highly catastrophic—corresponding to the formation of the Eunomia and Koronis families, respectively. As we already found for the intermediate energy regime represented by the Karin family, pre-shattered parent bodies always lead to outcome properties in better agreement with those of real families. In particular, the fragment size distribution obtained by disrupting a monolithic body always contains a large gap between the largest fragment and the next largest ones, whereas it is much more continuous in the case of a pre-shattered parent body. In the latter case, the ejection speeds of large fragments are also higher and a smaller impact energy is generally required to achieve a similar degree of disruption. Hence, unless the internal structure of bodies involved in a collision is known, predicting accurately the outcome is impossible. Interestingly, disrupting a pre-shattered parent body to reproduce the Koronis family yields a fragment size distribution characterized by four almost identical largest objects, as observed in the real family. This peculiar outcome has been found before in laboratory experiments but is obtained for the first time following gravitational re-accumulation. Finally, we show that material belonging to the largest fragments of a family originates from well-defined regions inside the parent body (the extent and location of which are dependent upon internal structure), despite the many gravitational interactions that occur during the re-accumulation process. Hence fragment formation does not proceed stochastically but results directly from the velocity field imparted during the impact.     

Fugitives from the Vesta family

Spectroscopic observations of Asteroid (4) Vesta and numerous members of the Vesta family located in the inner asteroid belt have determined that these objects have reflectance properties of basaltic material. A plausible hypothesis is that the surface of Vesta was punctured by large impacts in the past which dispersed fragments of its basaltic crust into space and produced one of the most prominent asteroid families ever created in the belt. Until recently, Vesta was the only known object in the asteroid belt which underwent differentiation and survived to the present epoch. Since 2000, many new small basaltic asteroids have been discovered in the inner and outer parts of the asteroid belt, possibly representing fragments from distinct differentiated bodies. These discoveries may help us to better understand the number and nature of objects in the inner Solar System that underwent geological differentiation. To investigate these issues we performed extensive numerical simulations whose aim was to reproduce, as precisely as possible, the dynamical evolution of Vesta's ejected fragments over timescales comparable to the family's age. Specifically, we numerically integrated the orbital evolution of 6600 test bodies with orbits that started within the Vesta family and dynamically evolved over 2 Gy. Our model included gravitational perturbation of all planets (except Mercury) and the Yarkovsky effect. The results show that a relatively large fraction of the original Vesta family members may have evolved out of the family borders defined by clustering algorithms and are now dispersed over the inner asteroid belt. We compared the orbital distribution of our model fragments with the orbital locations of known basaltic asteroids in various parts of the inner main belt to find that: (i) Most basaltic asteroids with semimajor axis located outside the Vesta family's borders in the inner main belt, including (809) Lundia and (956) Elisa, are most likely fugitives from the Vesta family that have evolved to their current orbits via various identified dynamical pathways. Our results also suggest that the Vesta family is at least ∼1 Gy old. (ii) Interestingly, orbits of many basaltic asteroids with , like those of (4796) Lewis and (5379) Abehiroshi, are displaced from the Vesta family to low inclinations and are not obtained in our simulations with sufficient efficiency. We propose that: (i) these small basaltic asteroids may be fragments of differentiated bodies other than (4) Vesta; or (ii) they were liberated from the Vesta's surface before (or during) the Late Heavy Bombardment epoch ∼3.8 Gy ago and their orbital inclinations separated from that of Vesta when secular resonances swept through the region.     

Catastrophic disruption of asteroids and family formation: a review of numerical simulations including both fragmentation and gravitational reaccumulations

In the last few years, thanks to the development of sophisticated numerical codes, a major breakthrough has been achieved in our understanding of the processes involved in small body collisions. In this review, we summarize the most recent results provided by numerical simulations, accounting for both the fragmentation of an asteroid and the gravitational interactions of the generated fragments. These studies have greatly improved our knowledge of the mechanisms that are at the origin of some observed features in the asteroid belt. In particular, the simulations have demonstrated that, for bodies larger than several kilometers, the collisional process not only involves the fragmentation of the asteroid but also the gravitational interactions between the ejected fragments. This latter mechanism can lead to the formation of large aggregates by gravitational reaccumulation of smaller fragments, and helps explain the presence of large members within asteroid families. Numerical simulations of the complete process have thus reproduced successfully for the first time the main properties of asteroid families, each formed by the disruption of a large parent body, and provided information on the possible internal structure of the parent bodies. A large amount of work remains necessary, however, to understand in deeper detail the physical process as a function of material properties and internal structures that are relevant to asteroids, and to determine in a more quantitative way the outcome properties such as fragment shapes and rotational states.     

Using Chebyshev polynomial interpolation to improve the computational efficiency of gravity models near an irregularly-shaped asteroid

In asteroid rendezvous missions, the dynamical environment near an asteroid's surface should be made clear prior to launch of the mission. However, most asteroids have irregular shapes,which lower the efficiency of calculating their gravitational field by adopting the traditional polyhedral method. In this work, we propose a method to partition the space near an asteroid adaptively along three spherical coordinates and use Chebyshev polynomial interpolation to represent the gravitational acceleration in each cell. Moreover, we compare four different interpolation schemes to obtain the best precision with identical initial parameters. An error-adaptive octree division is combined to improve the interpolation precision near the surface. As an example, we take the typical irregularly-shaped nearEarth asteroid 4179 Toutatis to demonstrate the advantage of this method; as a result, we show that the efficiency can be increased by hundreds to thousands of times with our method. Our results indicate that this method can be applicable to other irregularly-shaped asteroids and can greatly improve the evaluation efficiency.     

Asteroid families: Current situation

Being the products of energetic collisional events, asteroid families provide a fundamental body of evidence to test the predictions of theoretical and numerical models of catastrophic disruption phenomena. The goal is to obtain, from current physical and dynamical data, reliable inferences on the original disruption events that produced the observed families. The main problem in doing this is recognizing, and quantitatively assessing, the importance of evolutionary phenomena that have progressively changed the observable properties of families, due to physical processes unrelated to the original disruption events. Since the early 1990s, there has been a significant evolution in our interpretation of family properties. New ideas have been conceived, primarily as a consequence of the development of refined models of catastrophic disruption processes, and of the discovery of evolutionary processes that had not been accounted for in previous studies. The latter include primarily the Yarkovsky and Yarkovsky-O’Keefe-Radzvieski-Paddack (YORP) effects—radiation phenomena that can secularly change the semi-major axis and the rotation state. We present a brief review of the current state of the art in our understanding of asteroid families, point out some open problems, and discuss a few likely directions for future developments.     

The Yarkovsky effect is not responsible for small crater depletion on Eros and Itokawa

The near-Earth Asteroids Eros and Itokawa show a pronounced lack of small (?100 m) craters, the vast majority of which were formed during their time in the main belt, and this has been cited as possible evidence that small (?10 m) impactors are efficiently removed from the main belt by the Yarkovsky effect. Using well-tested models for the evolution of the main-belt size distribution and the evolution of crater populations on asteroid surfaces, I show that a pronounced lack of small impactors would require size-dependent removal far stronger than can result from the Yarkovsky effect (or any other known process). Furthermore, such strong removal would lead to wavelike perturbations in the main-belt and near-Earth asteroid size distributions that are inconsistent with their observed size distributions, as well as the cratering records on asteroid surfaces. A more likely explanation is that processes on asteroid surfaces, such as seismic shaking, are responsible for erasing small craters after they form.     

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.     

Dynamics of rotationally fissioned asteroids: Source of observed small asteroid systems

We present a model of near-Earth asteroid (NEA) rotational fission and ensuing dynamics that describes the creation of synchronous binaries and all other observed NEA systems including: doubly synchronous binaries, high-e binaries, ternary systems, and contact binaries. Our model only presupposes the Yarkovsky-O’Keefe-Radzievskii-Paddack (YORP) effect, “rubble pile” asteroid geophysics, and gravitational interactions. The YORP effect torques a “rubble pile” asteroid until the asteroid reaches its fission spin limit and the components enter orbit about each other (Scheeres, D.J. [2007]. Icarus 189, 370-385). Non-spherical gravitational potentials couple the spin states to the orbit state and chaotically drive the system towards the observed asteroid classes along two evolutionary tracks primarily distinguished by mass ratio. Related to this is a new binary process termed secondary fission - the secondary asteroid of the binary system is rotationally accelerated via gravitational torques until it fissions, thus creating a chaotic ternary system. The initially chaotic binary can be stabilized to create a synchronous binary by components of the fissioned secondary asteroid impacting the primary asteroid, solar gravitational perturbations, and mutual body tides. These results emphasize the importance of the initial component size distribution and configuration within the parent asteroid. NEAs may go through multiple binary cycles and many YORP-induced rotational fissions during their approximately 10 Myr lifetime in the inner Solar System. Rotational fission and the ensuing dynamics are responsible for all NEA systems including the most commonly observed synchronous binaries.     

A preliminary analysis of the orbit of the Mars Trojan asteroid (5261) Eureka

Observations and results of orbit determination of the first known Mars Trojan asteroid (5261) Eureka are presented. We have numerically calculated the evolution of the orbital elements, and have analyzed the behavior of the motion during the next 2 Myr. Strong perturbations by planets other than Mars seem to stabilize the eccentricity of the asteroid by stirring the high order resonances present in the elliptic restricted problem. As a result, the orbit appears stable at least on megayear timescales. The difference of the mean longitudes of Mars and Eureka and the semimajor axis of the asteroid form a pair of variables that essentially behave in an adiabatic manner, while the evolution of the other orbital elements is largely determined by the perturbations due to other planets.     

Spectral study of the Eunomia asteroid family : I. Eunomia

We present color ratio curves of the S-Asteroid 15 Eunomia, which have been extracted from high-precision photometric lightcurves obtained in three different VNIR wavelength bands at the Bochum Telescope, La Silla. The measured color ratio curves and near infrared spectra were used to derive a detailed surface composition model whose shape has been computed by V-lightcurve inversions. According to this analysis, the asteroid shows on one hemisphere a higher concentration of pyroxene, which causes an increased 440/700 nm and a reduced 940/700 nm reflectance ratio as well as a pronounced 2-μm absorption band. The remaining surface shows a higher concentration of olivine, leading to a reduced 440/700 nm and slightly increased 940/700 nm color ratio. In addition, we found that the maximum of the 440/700 nm color ratio curve coincide with the minimum of the 940/700 nm color ratio curve and vice versa. We demonstrate on the basis of USGS laboratory spectra that this anti-cyclical behavior can be explained by choosing Fe-rich olivine and a pyroxene with moderate Fe content as varying mineral phases. Furthermore, our observations confirm that 15 Eunomia is an irregular elongated and at least partially differentiated body. Previous spectral investigations of several smaller fragments of the Eunomia asteroid family revealed that the amount of fragments showing an increased pyroxene content exceeds the amount of pyroxene-poor fragments (Nathues, 2000, DLR Forschungsbericht, ISSN 1434-8454). This finding together with the observation that the major fraction of Eunomia's surface is enriched in olivine let us claim that a large fraction of the original pyroxene-enriched crust layer has been lost due to a major collision that created the Eunomia asteroid family. Significant spectral evidences, consistent with high concentrations of metals have been found neither in the rotational resolved spectra of 15 Eunomia nor in its fragments. This led to the conclusion that either no core consisting mainly of metals exists or that an eventual one has not been unearthed by the impact.     

Shapes and Pole Orientations of Asteroids (360) Carlova and (209) Dido

Several light-curves of asteroid (360) Carlova and (209) Dido in different epochs were analyzed to determine shapes and pole orientations by means of AM-method and least squares method. New values of a/b, b/c, λ _p and β _p for asteroid (360) Carlova were obtained, which are 1.52°, 1.5°, 120 ± 6° and 66 ± 7°, respectively. We report a first determination of the parameters of (209) Dido which are 1.3°, 1.1°, 221 ± 6° and 37 ± 3°, respectively.     

Secular perihelion advances of the inner planets and asteroid Icarus

A small effect expected from a recently proposed gravitational impact model (Wilhelm et al., 2013) is used to explain the remaining secular perihelion advance rates of the planets Mercury, Venus, Earth, Mars, and the asteroid (1566) Icarus—after taking into account the disturbances related to Newton’s Theory of Gravity. Such a rate was discovered by Le Verrier (1859) for Mercury and calculated by Einstein (1915, 1916) in the framework of his General Theory of Relativity (GTR). Accurate observations are now available for the inner Solar System objects with different orbital parameters. This is important, because it allowed us to demonstrate that the quantitative amount of the deviation from an 1/r potential is—under certain conditions—only dependent on the specific mass distribution of the Sun and not on the characteristics of the orbiting objects and their orbits. A displacement of the effective gravitational from the geometric centre of the Sun by about 4400 m towards each object is consistent with the observations and explains the secular perihelion advance rates.     

Astrometric masses of 21 asteroids, and an integrated asteroid ephemeris

We apply the technique of astrometric mass determination to measure the masses of 21 main-belt asteroids; the masses of 9 Metis (1.03 ± 0.24 × 10^-11 M_⊙), 17 Thetis (6.17 ± 0.64 × 10^-13 M_⊙), 19 Fortuna (5.41 ± 0.76 × 10^-12 M_⊙), and 189 Phthia (1.87 ± 0.64 × 10^-14 M_⊙) appear to be new. The resulting bulk porosities of 11 Parthenope (12±4%) and 16 Psyche (46±16%) are smaller than previously-reported values. Empirical expressions modeling bulk density as a function of mean radius are presented for the C and S taxonomic classes. To accurately model the forces on these asteroids during the mass determination process, we created an integrated ephemeris of the 300 large asteroids used in preparing the DE-405 planetary ephemeris; this new BC-405 integrated asteroid ephemeris also appears useful in other high-accuracy applications.     

On the Possible Hazard on the Major Cities Caused by Asteroid Impact in the Pacific Ocean – II

Earlier, it was pointed out that asteroids with diameters close to 200 m that hit the Pacific Ocean once in 10⁴ y create large tsunami waves which could cause serious hazards to the cities and countries lying along the periphery. Here, a more detailed investigation is carried out in relation to tsunami hazards created by earthquakes. It is shown that the energy involved in the asteroid induced tsunami wave is some 300 times greater than large tsunami waves and 100 times greater than the Chile tsunami of 1960, the largest this century. Second, if the impact should take place in the Pacific south of Japan, so that the tsunami wave enters Tokyo Bay, the damage could cost 800 billion dollars, which is 6 times greater than the damage due to the Hanshin (Kōbe-Osaka) earthquake of 1995. This revised version was published online in July 2006 with corrections to the Cover Date.     

Stability of Surface Motion on a Rotating Ellipsoid

The dynamical environment on the surface of a rotating, massive ellipsoid is studied, with applications to surface motion on an asteroid. The analysis is performed using a combination of classical dynamics and geometrical analysis. Due to the small sizes of most asteroids, their shapes tend to differ from the classical spheroids found for the planets. The tri-axial ellipsoid model provides a non-trivial approximation of the gravitational potential of an asteroid and is amenable to analytical computation. Using this model, we study some properties of motion on the surface of an asteroid. We find all the equilibrium points on the surface of a rotating ellipsoid and we show that the stability of these points is intimately tied to the conditions for a Jacobi or MacLaurin ellipsoid of equilibria. Using geometrical analysis we can define global constraints on motion as a function of shape, rotation rate, and density, we find that some asteroids should have accumulation of material at their ends, while others should have accumulation of surface material at their poles. This study has implications for motion of a rover on an asteroid, and for the distribution of natural material on asteroids, and for a spacecraft hovering over an asteroid.     

Spatially resolved spectral observations of Asteroid 951 Gaspra

On October 29, 1991 the Galileo spacecraft encountered Asteroid 951 Gaspra with a telescopic CCD camera and a near infrared mapping spectrometer that provided the first optically resolved views of any asteroid. Data from these two sensors were combined to detect the spectral signature of iron bearing minerals on this S-type asteroid. A minimum of two spectral units were identified on 951 Gaspra, both containing a higher relative abundance of olivine than those found in ordinary chondrites. These data indicate that this S asteroid is an object that has undergone igneous differentiation processes. A 2.7 μm spectral feature was also detected on the surface of 951 Gaspra and may be due to the presence of structural OH.