Formation of Asteroid Families by Catastrophic Disruption: Simulations with Fragmentation and Gravitational Reaccumulation

This paper builds on preliminary work in which numerical simulations of the collisional disruption of large asteroids (represented by the Eunomia and Koronis family parent bodies) were performed and which accounted not only for the fragmentation of the solid body through crack propagation, but also for the mutual gravitational interaction of the resulting fragments. It was found that the parent body is first completely shattered at the end of the fragmentation phase, and then subsequent gravitational reaccumulations lead to the formation of an entire family of large and small objects with dynamical properties similar to those of the parent body. In this work, we present new and improved numerical simulations in detail. As before, we use the same numerical procedure, i.e., a 3D SPH hydrocode to compute the fragmentation phase and the parallel N-body code pkdgrav to compute the subsequent gravitational reaccumulation phase. However, this reaccumulation phase is now treated more realistically by using a merging criterion based on energy and angular momentum and by allowing dissipation to occur during fragment collisions. We also extend our previous studies to the as yet unexplored intermediate impact energy regime (represented by the Flora family formation) for which the largest fragment's mass is about half that of the parent body. Finally, we examine the robustness of the results by changing various assumptions, the numerical resolution, and different numerical parameters. We find that in the lowest impact energy regime the more realistic physical approach of reaccumulation leads to results that are statistically identical to those obtained with our previous simplistic approach. Some quantitative changes arise only as the impact energy increases such that higher relative velocities are reached during fragment collisions, but they do not modify the global outcome qualitatively. As a consequence, these new simulations confirm previous main results and still lead to the conclusion that: (1) all large family members must be made of gravitationally reaccumulated fragments; (2) the original fragment size distribution and their orbital dispersion are respectively steeper and smaller than currently observed for the real families, supporting recent studies on subsequent evolution and diffusion of family members; and (3) the formation of satellites around family members is a frequent and natural outcome of collisional processes.     

Numerical simulations of impacts involving porous bodies: II. Comparison with laboratory experiments

In this paper, we compare the outcome of high-velocity impact experiments on porous targets, composed of pumice, with the results of simulations by a 3D SPH hydrocode in which a porosity model has been implemented. The different populations of small bodies of our Solar System are believed to be composed, at least partially, of objects with a high degree of porosity. To describe the fragmentation of such porous objects, a different model is needed than that used for non-porous bodies. In the case of porous bodies, the impact process is not only driven by the presence of cracks which propagate when a stress threshold is reached, it is also influenced by the crushing of pores and compaction. Such processes can greatly affect the whole body's response to an impact. Therefore, another physical model is necessary to improve our understanding of the collisional process involving porous bodies. Such a model has been developed recently and introduced successfully in a 3D SPH hydrocode [Jutzi, M., Benz, W., Michel, P., 2008. Icarus 198, 242-255]. Basic tests have been performed which already showed that it is implemented in a consistent way and that theoretical solutions are well reproduced. However, its full validation requires that it is also capable of reproducing the results of real laboratory impact experiments. Here we present simulations of laboratory experiments on pumice targets for which several of the main material properties have been measured. We show that using the measured material properties and keeping the remaining free parameters fixed, our numerical model is able to reproduce the outcome of these experiments carried out under different impact conditions. This first complete validation of our model, which will be tested for other porous materials in the future, allows us to start addressing problems at larger scale related to small bodies of our Solar System, such as collisions in the Kuiper Belt or the formation of a family by the disruption of a porous parent body in the main asteroid belt.     

Fragment ejection velocities and the collisional evolution of asteroids

We have applied the algorithm developed by Petit and Farinella (Celest. Mech. 57, 1–28, 1993) to model the outcomes of impacts between asteroids of different sizes, to show that a crucial feature of these models is the assumed relationship between velocity and mass of fragments ejected after a shattering impact. Not only how the mean velocity depends upon mass is important to determine the extent of fragment reaccumulation, but also the distribution of velocities about the mean values. The available experimental evidence on this issue is still sparse, and does not constrain the collisional models well enough to allow us to make reliable predictions on the outcomes of impacts between bodies of size much larger than the laboratory targets. As a consequence, when the collisional outcome models are used as an input for simulations of the asteroid collisional history since the origin of the solar system, the results show a strong sensitivity to the assumed velocity vs mass relationship. This sensitivity is stronger in the diameter range (a few tens to a few hundreds of km) where the self-gravitational reaccumulation of fragments is most effective, but may also extend to much smaller sizes.     

Hydrocode modeling of oblique impacts: The fate of the projectile

Abstract— All impacts are oblique to some degree. Only rarely do projectiles strike a planetary surface (near) vertically. The effects of an oblique impact event on the target are well known, producing craters that appear circular even for low impact angles (>15° with respect to the surface). However, we still have much to learn about the fate of the projectile, especially in oblique impact events. This work investigates the effect of angle of impact on the projectile. Sandia National Laboratories' three‐dimensional hydrocode CTH was used for a series of high‐resolution simulations (50 cells per projectile radius) with varying angle of impact. Simulations were carried out for impacts at 90, 60, 45, 30, and 15° from the horizontal, while keeping projectile size (5 km in radius), type (dunite), and impact velocity (20 km/s) constant. The three‐dimensional hydrocode simulations presented here show that in oblique impacts the distribution of shock pressure inside the projectile (and in the target as well) is highly complex, possessing only bilateral symmetry, even for a spherical projectile. Available experimental data suggest that only the vertical component of the impact velocity plays a role in an impact. If this were correct, simple theoretical considerations indicate that shock pressure, temperature, and energy would depend on sin²θ, where θ is the angle of impact (measured from the horizontal). However, our numerical simulations show that the mean shock pressure in the projectile is better fit by a sin θ dependence, whereas shock temperature and energy depend on sin^3/2 θ. This demonstrates that in impact events the shock wave is the result of complex processes that cannot be described by simple empirical rules. The mass of shock melt or vapor in the projectile decreases drastically for low impact angles as a result of the weakening of the shock for decreasing impact angles. In particular, for asteroidal impacts the amount of projectile vaporized is always limited to a small fraction of the projectile mass. In cometary impacts, however, most of the projectile is vaporized even at low impact angles. In the oblique impact simulations a large fraction of the projectile material retains a net downrange motion. In agreement with experimental work, the simulations show that for low impact angles (30 and 15°), a downrange focusing of projectile material occurs, and a significant amount of it travels at velocities larger than the escape velocity of Earth.     

Scaling of oblique impacts in frictional targets: Implications for crater size and formation mechanisms

Almost every meteorite impact occurs at an oblique angle of incidence, yet the effect of impact angle on crater size or formation mechanism is only poorly understood. This is, in large part, due to the difficulty of inferring impactor properties, such as size, velocity and trajectory, from observations of natural craters, and the expense and complexity of simulating oblique impacts using numerical models. Laboratory oblique impact experiments and previous numerical models have shown that the portion of the projectile’s kinetic energy that is involved in crater excavation decreases significantly with impact angle. However, a thorough quantification of planetary-scale oblique impact cratering does not exist and the effect of impact angle on crater size is not considered by current scaling laws. To address this gap in understanding, we developed iSALE-3D, a three-dimensional multi-rheology hydrocode, which is efficient enough to perform a large number of well-resolved oblique impact simulations within a reasonable time. Here we present the results of a comprehensive numerical study containing more than 200 three-dimensional hydrocode-simulations covering a broad range of projectile sizes, impact angles and friction coefficients. We show that existing scaling laws in principle describe oblique planetary-scale impact events at angles greater than 30° measured from horizontal. The displaced mass of a crater decreases with impact angle in a sinusoidal manner. However, our results indicate that the assumption that crater size scales with the vertical component of the impact velocity does not hold for materials with a friction coefficient significantly lower than 0.7 (sand). We found that increasing coefficients of friction result in smaller craters and a formation process more controlled by impactor momentum than by energy.     

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.     

High-velocity dust collisions: forming planetesimals in a fragmentation cascade with final accretion

In laboratory experiments we determine the mass gain and loss in central collisions between centimetre- to decimetre-size SiO₂ dust targets and submillimetre- to centimetre-size SiO₂ dust projectiles of varying mass, size, shape and at different collision velocities up to ∼56.5 m s⁻¹. Dust projectiles much larger than 1 mm lead to a small amount of erosion of the target but decimetre targets do not break up. Collisions produce ejecta, which are smaller than the incoming projectile. Projectiles smaller than 1 mm are accreted by a target even at the highest collision velocities. This implies that net accretion of decimetre and larger bodies is possible. Independent of the original size of a considered projectile, after several collisions, all fragments will be of submillimetre size which might then be (re)accreted in the next collision with a larger body. The experimental data suggest that collisional growth through fragmentation and reaccretion is a viable mechanism to form planetesimals.     

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 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.     

Asteroidal catastrophic collisions simulated by hypervelocity impact experiments

We report the results of six impact fragmentation experiments carried out with free-falling macroscopic targets of different compositions and shapes, and with projectile velocities close to 9 km/sec, i.e., significantly higher than the sound velocity in the target materials. The data have been examined by deriving the mass and shape distributions of the fragments, by reconstructing two of the shattered targets in order to study the geometry of the fracture surfaces, and by analyzing the properties of the fine-grained high-velocity ejecta. The fragment mass distributions show clearly that the degree of target fragmentation depends strongly on the impact parameter. Apart from the few largest fragments, these distributions are well represented by two power laws with different exponents, connected at a size of about 1 cm. The fragment shapes are generally in good agreement with those observed in previous experiments, and no significant shape vs size dependence has been found down to sizes of the order of 0.1 mm. The fragments tend to become larger and possibly more irregular in shape when they are generated farther from the impact point. The fracture surfaces are oriented roughly along meridians and parallels (with the pole at the impact point) when the target is spherical, but are clustered around the symmetry planes when the target is ellipsoidal. Fine-grained particles, with typical sizes and velocities of 0.01 cm and 1 km/sec, respectively, are ejected at low-elevation angles and in a rather collimated way, starting both from the neighborhood of the impact point and from regions of incipient cracking. Particular attention has been paid to a comparison between these results and the observed properties of the outcomes of asteroidal catastrophic collisions, like the dynamical families and the small inner planet crossing objects. While the collisional theory for the origin of families is fully consistent with the experimental results (with some indication for a significant role of the parent asteroid's self-gravitation), the elongated shapes of several Apollo-Amor objects are much rarer among the laboratory fragments, and thus appear to require a different explanation.     

Impact ejection,spallation, and the origin of meteorites

Recent discoveries suggest that some meteorites have originated from major planets or satellites. Although it has been suggested that a large primary impact event might eject rock fragments as secondaries, it was previously supposed that material ejected at several kilometers per second would be highly shocked or perhaps melted. It is shown that a small amount of material (0.01 to 0.05 projectile mass) may be ejected at high velocity shock pressures. The approach utilizes observations of stress-wave propagation from large underground explosions to predict stresses and particle velocities in the near-surface environment. The largest fragments ejected at any velocity are spalls that originate from the target planet's surface. The spall size is proportional to the radius of the primary impactor and the target tensile strength and inversely proportional to ejection velocity. The shock level in the spalls is low, typically half of the dynamic crushing strength of the rock. The model also predicts the aspect ratio of the spalled fragments, the angle of ejection, and the sizes and shock level of other fragments originating deeper in the target. Comparison with data from laboratory experiments, the Ries Crater, and secondary crater sizes shows generally good agreement, although the observed fragment size at ejection velocities greater than 1 km/sec is considerably smaller than the simple version of the theory predicts. The theory indicates that although significant masses of solid material could be ejected from the Moon or Mars by large meteorite impacts, the fragments ejected from ca. 30-km-diameter craters are at most a few tens of meters in diameter if the most optimistic assumptions are made. The maximum fragment diameter is more likely to be about a meter. This theory, however, applies rigorously only up to ejection velocities of ca 1 km/sec. Further numerical extensions are necessary before film conclusions can be drawn, especially for Martian ejecta.     

Reconstruction of the Morasko meteoroid impact—Insight from numerical modeling

The Morasko strewn field located near Poznań, Poland comprises seven impact craters with diameters ranging from 20 to 90 m, all of which were formed in glacial sediments around 5000 yr ago. Numerous iron meteorites have been recovered in the area and their distribution suggests a projectile with the trajectory from NE to SW. Similar impact events producing crater strewn fields on average happen every 500 yr and pose a serious risk for modern civilization, which is why it is of utmost importance to study terrestrial strewn fields in detail. In this work, we investigate the Morasko meteoroid passage through the atmosphere, the distribution of its fragments on the ground, and the process of forming individual craters by means of numerical modeling. By combining atmospheric entry modeling, Pi‐group scaling of transient crater size and hydrocode simulations of impact processes, we constructed a comprehensive model of the Morasko strewn field formation. We determined the preatmospheric parameters of the Morasko meteoroid. The entry mass is between 600 and 1100 tons, the velocity range is between 16 and 18 km s^?1, and the trajectory angle is 30–40^°. Such entry velocities and trajectory angles do not deviate from typical values for near‐Earth asteroids, although the initial mass we determined can be considered as small. Our studies on velocities and masses of crater‐forming fragments showed that the biggest Morasko crater was formed by a projectile about 1.5 m in diameter with the impact velocity ~10 km s^?1. Environmental consequences of the Morasko impact event are very localized.     

Velocity of finer fragments from impact

We describe the method and the result of a new experiment to obtain velocity distribution of fine ejecta fragments, from a few to a hundred microns in size, produced from basalt targets by impacts of nylon projectiles at a velocity of 3.7 km s⁻¹. The size distribution of holes perforated by the ejecta fragments on thin films and foils placed around the targets was investigated, and the size-velocity relation was determined with the aid of an empirical formula for threshold penetration (McDonnell and Sullivan, Hypervelocity Impacts in Space, Unit for Space Sciences, University of Kent, 1992). The velocity of the fastest fragments, at a given size, is from the extrapolation of the size-velocity relation for 1–100 mm fragments (Nakamura and Fujiwara, Icarus92, 132–146, 1991; Nakamura et al, Icarus100, 127–135, 1992). The laboratory results are also compared with those obtained from the study of secondary craters around large lunar craters (Vickery, Icarus67, 224–236, 1986, Geophys. Res. Lett. 14, 726–729, 1987). All these data provide a smooth size-velocity relationship in the normalized fragment size range of four orders of magnitude.     

Cumulative damage in strength-dominated collisions of rocky asteroids: Rubble piles and brick piles

Laboratory impact experiments were performed to investigate the conditions that produce large-scale damage in rock targets. Aluminum cylinders (6.3 mm diameter) impacted basalt cylinders (69 mm diameter) at speeds ranging from 0.7 to 2.0 km/s. Diagnostics included measurements of the largest fragment mass, velocities of the largest remnant and large fragments ejected from the periphery of the target, and X-ray computed tomography imaging to inspect some of the impacted targets for internal damage. Significant damage to the target occurred when the kinetic energy per unit target mass exceeded roughly of the energy required for catastrophic shattering (where the target is reduced to one-half its original mass). Scaling laws based on a rate-dependent strength were developed that provide a basis for extrapolating the results to larger strength-dominated collisions. The threshold specific energy for widespread damage was found to scale with event size in the same manner as that for catastrophic shattering. Therefore, the factor of four difference between the two thresholds observed in the lab also applies to larger collisions. The scaling laws showed that for a sequence of collisions that are similar in that they produce the same ratio of largest fragment mass to original target mass, the fragment velocities decrease with increasing event size. As a result, rocky asteroids a couple hundred meters in diameter should retain their large ejecta fragments in a jumbled rubble-pile state. For somewhat larger bodies, the ejection velocities are sufficiently low that large fragments are essentially retained in place, possibly forming ordered “brick-pile” structures.     

Impact and explosion crater ejecta,fragment size,and velocity

A model was developed for the mass distribution of fragments that are ejected at a given velocity for impact and explosion craters. The model is semiempirical in nature and is derived from (1) numerical calculations of cratering and the resultant mass versus ejection velocity, (2) observed ejecta blanket particle size distributions, (3) an empirical relationships between maximum ejecta fragment size and crater diameter, (4) measurements of maximum ejecta size versus ejecta velocity, and (5) an assumption on the functional form for the distribution of fragments ejected at a given velocity. This model implies that for planetary impacts into competent rock, the distribution of fragments ejected at a given velocity is broad; e.g., 68% of the mass of the ejecta at a given velocity contains fragments having a mass less than 0.1 times a mass of the largest fragment moving at that velocity. Using this model, we have calculated the largest fragment that can be ejected from asteroids, the Moon, Mars, and Earth as a function of crater diameter. The model is unfortunately dependent on the size-dependent ejection velocity limit for which only limited data are presently available from photography of high explosive-induced rock ejecta. Upon formation of a 50-km-diameter crater on an atmosphereless planet having the planetary gravity and radius of the Moon, Mars, and Earth, fragments having a maximum mean diameter of ≈30, 22, and 17 m could be launched to escape velocity in the ejecta cloud. In addition, we have calculated the internal energy of ejecta versus ejecta velocity. The internal energy of fragments having velocities exceeding the escape velocity of the moon (～2.4 km/sec) will exceed the energy required for incipient melting for solid silicates and thus, the fragments ejected from Mars and the Earth would be melted.     

High-resolution simulations of clump–clump collisions using SPH with particle splitting

We investigate, by means of numerical simulations, the phenomenology of star formation triggered by low-velocity collisions between low-mass molecular clumps. The simulations are performed using a smoothed particle hydrodynamics code which satisfies the Jeans condition by invoking on-the-fly particle splitting.
Clumps are modelled as stable truncated (non-singular) isothermal, i.e. Bonnor–Ebert, spheres. Collisions are characterized by M ₀ (clump mass), b (offset parameter, i.e. ratio of impact parameter to clump radius) and     (Mach number, i.e. ratio of collision velocity to effective post-shock sound speed). The gas subscribes to a barotropic equation of state, which is intended to capture (i) the scaling of pre-collision internal velocity dispersion with clump mass, (ii) post-shock radiative cooling and (iii) adiabatic heating in optically thick protostellar fragments.
The efficiency of star formation is found to vary between 10 and 30 per cent in the different collisions studied and it appears to increase with decreasing M ₀, and/or decreasing b , and/or increasing     . For   b < 0.5  collisions produce shock-compressed layers which fragment into filaments. Protostellar objects then condense out of the filaments and accrete from them. The resulting accretion rates are high,     , for the first     . The densities in the filaments,     , are sufficient that they could be mapped in NH₃ or CS line radiation, in nearby star formation regions.     

Karin cluster formation by asteroid impact

Insights into collisional physics may be obtained by studying the asteroid belt, where large-scale collisions produced groups of asteroid fragments with similar orbits and spectra known as the asteroid families. Here we describe our initial study of the Karin cluster, a small asteroid family that formed 5.8±0.2 Myr ago in the outer main belt. The Karin cluster is an ideal ‘natural laboratory’ for testing the codes used to simulate large-scale collisions because the observed fragments produced by the 5.8-Ma collision suffered apparently only limited dynamical and collisional erosion. To date, we have performed more than 100 hydrocode simulations of impacts with non-rotating monolithic parent bodies. We found good fits to the size-frequency distribution of the observed fragments in the Karin cluster and to the ejection speeds inferred from their orbits. These results suggest that the Karin cluster was formed by a disruption of an ≈33-km-diameter asteroid, which represents a much larger parent body mass than previously estimated. The mass ratio between the parent body and the largest surviving fragment, (832) Karin, is ≈0.15-0.2, corresponding to a highly catastrophic event. Most of the parent body material was ejected as fragments ranging in size from yet-to-be-discovered sub-km members of the Karin cluster to dust grains. The impactor was ≈5.8 km across. We found that the ejections speeds of smaller fragments produced by the collision were larger than those of the larger fragments. The mean ejection speeds of >3-km-diameter fragments were . The model and observed ejection velocity fields have different morphologies perhaps pointing to a problem with our modeling and/or assumptions. We estimate that ∼5% of the large asteroid fragments created by the collision should have satellites detectable by direct imaging (separations larger than 0.1 arcsec). We also predict a large number of ejecta binary systems with tight orbits. These binaries, located in the outer main belt, could potentially be detected by lightcurve observations. Hydrocode modeling provides important constraints on the interior structure of asteroids. Our current work suggests that the parent asteroid of the Karin cluster may have been an unfractured (or perhaps only lightly fractured) monolithic object. Simulations of impacts into fractured/rubble pile targets were so far unable to produce the observed large gap between the first and second largest fragment in the Karin cluster, and the steep slope at small sizes (≈6.3 differential index). On the other hand, the parent asteroid of the Karin cluster was produced by an earlier disruptive collision that created the much larger, Koronis family some 2-3 Gyr ago. Standard interpretation of hydrocode modeling then suggests that the parent asteroid of the Karin cluster should have been formed as a rubble pile from Koronis family debris. We discuss several solutions to this apparent paradox.     

Impacts into weak dust targets under microgravity and the formation of planetesimals

We carried out 16 collision experiments in the drop tower in Bremen, Germany. Dust projectiles and solid projectiles of several mm in size impacted a dust target 5 cm in depth and width at velocities between 3.5 and 21.5 m/s. For solid impactors we found significant mass loss on the front (impact) side of the target. Mass loss depended on the impact velocity and projectile type (solid sphere or dust) and was up to 35 times the projectile mass for targets of the lowest tensile strength. Typical fragment velocities on the front side of the target ranged from 3 to 12 cm/s. The ejecta velocity was independent of the impact velocity but it increased with projectile mass. On the back side of the target (opposite to the impact side) mass was ejected from the target above a certain threshold impact velocity. Ejection velocity on the back side increased with impact velocity and is larger for solid projectiles than for dust projectiles. In one case a slightly stronger target gained mass in a slow dust-dust collision. We verified that collisions of dust projectiles with compact, very strong dust targets lead to a more massive target accreting part of the projectile. Applied to planetesimal formation, the experiments suggest that the maximum possible ejecta velocity from a body of several cm in size after a collision is small. Ejecta were slow enough that they were reaccreted by means of gas flow if large pores were part of the body's morphology. While very weak bodies cannot grow in the primary collision at the given velocities, this can lead to growth by secondary collisions. Slight compression, which could result from preceding collisions, might lead to immediate growth of a body in slow collisions by adding projectile mass.     

Cratering erosion of planetary embryos

I have performed 3D numerical hydrodynamic simulations of impacts of stony projectiles on stony planar targets in a range of impact velocities from 1.25 to 60 km/s. The projectile and target masses ejected at speeds greater than some given values have been calculated. This provided a possibility to determine impact erosion of a target which undergoes bombardment with comparatively small bodies. The relative losses of target masses and masses of retained projectile material have been averaged over impact angles and approximated by analytical formulas as functions of impact and escape velocities. The balance between escaped material of a target and retained material of a projectile determines growth or reduction of a target mass. The target cratering erosion predominates over the projectile retention when the impacts have velocities of more than 3-5 times the escape velocity of a target. The results can be applied to collisions of planetary embryos with planetesimals, which have higher velocities than embryo-embryo impacts. Estimates for impact velocities 1-10 km/s show that while large embryos accrete planetesimals smaller embryos erode and can completely vanish or partly lose their silicate shells if they are differentiated. Application of calculated erosion efficiency to Mercury made it possible to test a hypothesis (Vityazev, A.V., Pechernikova, G.V., Safronov, V.S. [1988]. Formation of Mercury and removal of its silicate shell. In: Vilas, F., Chapman, C.R., Matthews, M.S. (Eds.), Mercury. Univ. Arizona Press., Tucson, pp. 667−669) that differentiated massive proto-Mercury has lost its mantle due to collisions with objects of moderate sizes. It turned out that in order for this to happen, relative collision velocities must exceed 25 km/s. As alternatives to the widely-known hypothesis of a giant impact on a massive proto-Mercury, other possibilities are considered, which do not require such high speeds. The first one is formation of a number of small-sized metal-rich embryos which lose their silicate shells due to cratering erosion. The second is that a small proto-Mercury was metallic and gained its mantle at the latest stage of its accumulation when it grew so large that the erosion became ineffective.     

Modeling global impact effects on middle-sized icy bodies: applications to Saturn's moons

Disrupted terrains that form as a consequence of giant impacts may help constrain the internal structures of planets, asteroids, comets and satellites. As shock waves and powerful seismic stress waves propagate through a body, they interact with the internal structure in ways that may leave a characteristic impression upon the surface. Variations in peak surface velocity and tensile stress, related to landform degradation and surface rupture, may be controlled by variations in core size, shape and density. Caloris Basin on Mercury and Imbrium Basin on the Moon have disturbed terrain at their antipodes, where focusing is most intense for an approximately symmetric spheroid. Although, the icy saturnian satellites Tethys, Mimas, and Rhea possess giant impact structures, it is not clear whether these structures have correlated disrupted terrains, antipodal or elsewhere. In anticipation of high-resolution imagery from Cassini, we investigate antipodal focusing during giant impacts using a 3D SPH impact model. We first investigate giant impacts into a fiducial 1000 km diameter icy satellite with a variety of core radii and compositions. We find that antipodal disruption depends more on core radius than on core density, suggesting that core geometry may express a surface signature in global impacts on partially differentiated targets. We model Tethys, Mimas, and Rhea according to their image-derived shapes (triaxial for Tethys and Mimas and spherical for Rhea), varying core radii and densities to give the proper bulk densities. Tethys shows greater antipodal values of peak surface velocity and peak surface tensile stress, indicating more surface damage, than either Mimas or Rhea. Results for antipodal and global fragmentation and terrain rupture are inconclusive, with the hydrocode itself producing global disruption at the limits of model resolution but with peak fracture stresses never exceeding the strength of laboratory ice.