Search

Highlights • We study the collisional disruption of 100 km- diameter asteroids. • The resulting size-frequency (SFD) and velocity distributions are computed. • The effects of pore-crushing as well as friction are investigated. • The porous targets have a significantly higher impact strength than the rubble-pile parent bodies investigated previously. • SFDs resulting from a collision with a given specific energy are strongly dependent on the size scale. Abstract Understanding the collisional fragmentation and subsequent reaccumulation of fragments is crucial for studies of the formation and evolution of the small-body populations. Using an SPH / N-body approach, we investigate the size-frequency distributions (SFDs) resulting from the disruption of 100 km-diameter targets consisting of porous material, including the effects of pore-crushing as well as friction. Overall, the porous targets have a significantly higher impact strength (Q D *) than the rubble-pile parent bodies investigated previously (Benavidez et al., 2012) and show a behavior more similar to non-porous monolithic targets (Durda et al., 2007). Our results also confirm that for a given specific impact energy, the SFDs resulting from a parent body disruption are strongly dependent on the size scale. [ABSTRACT FROM AUTHOR]

We model large-scale ( ≈ 2000 km) impacts on a Mars-like planet using a Smoothed Particle Hydrodynamics code. The effects of material strength and of using different Equations of State on the post-impact material and temperature distributions are investigated. The properties of the ejected material in terms of escaping and disc mass are analysed as well. We also study potential numerical effects in the context of density discontinuities and rigid body rotation. We find that in the large-scale collision regime considered here (with impact velocities of 4 km/s), the effect of material strength is substantial for the post-impact distribution of the temperature and the impactor material, while the influence of the Equation of State is more subtle and present only at very high temperatures. [ABSTRACT FROM AUTHOR]

A giant impact is commonly thought to explain the dramatic contrast in elevation and crustal thickness between the two hemispheres of Mars known as the "Martian Dichotomy". Initially, this scenario referred to an impact in the northern hemisphere that would lead to a huge impact basin (dubbed the "Borealis Basin"), while more recent work has instead suggested a hybrid origin that produces the Dichotomy through impact-induced crust-production. The majority of these studies have relied upon impact scaling-laws inaccurate at such large-scales, however, and those that have included realistic impact models have utilised over-simplified geophysical models and neglected any material strength. Here we use a large suite of strength-including smoothed-particle hydrodynamics (SPH) impact simulations coupled with a more sophisticated geophysical scheme of crust production and primordial crust to simultaneously investigate the feasibility of a giant impact on either hemisphere of Mars to have produced its dichotomous crust distribution, and utilise spherical harmonic analysis to identify the best-fitting cases. We find that the canonical Borealis-forming impact is not possible without both excessive crust production and strong antipodal effects not seen on Mars' southern hemisphere today. Our results instead favour an impact and subsequent localised magma ocean in the southern hemisphere that results in a thicker crust than the north upon crystallisation. Specifically, our best-fitting cases suggest that the projectile responsible for the Dichotomy-forming event was of radius 500–750 km, and collided with Mars at an impact angle of 15–30 ° with a velocity of 1.2–1.4 times mutual escape speed (∼ 6–7 km/s). • Impact-induced crust production is significant in large-scale impacts. • If an impact caused the Martian Dichotomy, it must have occurred in the south. • A 1000–1500 km body impacting at an angle 15–30° and speed 6–7 km/s is most likely. • Material strength is important for planetary-scale impacts and cannot be neglected. • Spherical harmonics are an effective tool in analysing impact simulation results. [ABSTRACT FROM AUTHOR]

Highlights: [•] This work aims at supporting studies of the kinetic impactor concept. [•] We study numerically the momentum transferred by impacts on porous asteroids. [•] The effect of microporosity and/or macroporosity is investigated. [•] The momentum multiplication factor β is calculated for various impact conditions. [•] Our results confirm that β is small (β <2 for vimp<15km/s) for porous targets. [Copyright &y& Elsevier]

Abstract: The Veritas family is located in the outer main belt and is named after its apparent largest constituent, Asteroid (490) Veritas. The family age has been estimated by two independent studies to be quite young, around 8Myr. Therefore, current properties of the family may retain signatures of the catastrophic disruption event that formed the family. In this paper, we report on our investigation of the formation of the Veritas family via numerical simulations of catastrophic disruption of a 140-km-diameter parent body, which was considered to be made of either porous or non-porous material, and a projectile impacting at 3 or 5km/s with an impact angle of 0° or 45°. Not one of these simulations was able to produce satisfactorily the estimated size distribution of real family members. Based on previous studies devoted to either the dynamics or the spectral properties of the Veritas family, which already treated (490) Veritas as a special object that may be disconnected from the family, we simulated the formation of a family consisting of all members except that asteroid. For that case, the parent body was smaller (112km in diameter), and we found a remarkable match between the simulation outcome, using a porous parent body, and the real family. Both the size distribution and the velocity dispersion of the real reduced family are very well reproduced. On the other hand, the disruption of a non-porous parent body does not reproduce the observed properties very well. This is consistent with the spectral C-type of family members, which suggests that the parent body was porous and shows the importance of modeling the effect of this porosity in the fragmentation process, even if the largest members are produced by gravitational reaccumulation during the subsequent gravitational phase. As a result of our investigations, we conclude that it is very likely that the Asteroid (490) Veritas and probably several other small members do not belong to the family as originally defined, and that the definition of this family should be revised. Further investigations will be performed to better constrain the definitions and properties of other asteroid families of different types, using the appropriate model of fragmentation. The identification of very young families in turn will continue to serve as a tool to check the validity of numerical models. [Copyright &y& Elsevier]

Abstract: Numerical simulations of asteroid breakups, including both the fragmentation of the parent body and the gravitational interactions between the fragments, have allowed us to reproduce successfully the main properties of asteroid families formed in different regimes of impact energy, starting from a non-porous parent body. In this paper, using the same approach, we concentrate on a single regime of impact energy, the so-called catastrophic threshold usually designated by , which results in the escape of half of the target’s mass. Thanks to our recent implementation of a model of fragmentation of porous materials, we can characterize for both porous and non-porous targets with a wide range of diameters. We can then analyze the potential influence of porosity on the value of , and by computing the gravitational phase of the collision in the gravity regime, we can characterize the collisional outcome in terms of the fragment size and ejection speed distributions, which are the main outcome properties used by collisional models to study the evolutions of the different populations of small bodies. We also check the dependency of on the impact speed of the projectile. In the strength regime, which corresponds to target sizes below a few hundreds of meters, we find that porous targets are more difficult to disrupt than non-porous ones. In the gravity regime, the outcome is controlled purely by gravity and porosity in the case of porous targets. In the case of non-porous targets, the outcome also depends on strength. Indeed, decreasing the strength of non-porous targets make them easier to disrupt in this regime, while increasing the strength of porous targets has much less influence on the value of . Therefore, one cannot say that non-porous targets are systematically easier or more difficult to disrupt than porous ones, as the outcome highly depends on the assumed strength values. In the gravity regime, we also confirm that the process of gravitational reaccumulation is at the origin of the largest remnant’s mass in both cases. We then propose some power-law relationships between and both target’s size and impact speed that can be used in collisional evolution models. The resulting fragment size distributions can also be reasonably fitted by a power-law whose exponent ranges between −2.2 and −2.7 for all target diameters in both cases and independently on the impact velocity (at least in the small range investigated between 3 and 5km/s). Then, although ejection velocities in the gravity regime tend to be higher from porous targets, they remain on the same order as the ones from non-porous targets. [Copyright &y& Elsevier]

Abstract: 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. [Copyright &y& Elsevier]

Abstract: In this paper, we extend our Smooth Particle Hydrodynamics (SPH) impact code to include the effect of porosity at a sub-resolution scale by adapting the so-called model. Many small bodies in the different populations of asteroids and comets are believed to contain a high degree of porosity and the determination of both their collisional evolution and the outcome of their disruption requires that the effect of porosity is taken into account in the computation of those processes. Here, we present our model and show how porosity interfaces with the elastic-perfectly plastic material description and the brittle fracture model generally used to simulate the fragmentation of non-porous rocky bodies. We investigate various compaction models and discuss their suitability to simulate the compaction of (highly) porous material. Then, we perform simple test cases where we compare results of the simulations to the theoretical solutions. We also present a Deep Impact-like simulation to show the effect of porosity on the outcome of an impact. Detailed validation tests will be presented in a next paper by comparison with high-velocity laboratory experiments on porous materials [Jutzi et al., in preparation]. Once validated at small scales, our new impact code can then be used at larger scales to study impacts and collisions involving brittle solids including porosity, such as the parent bodies of C-type asteroid families or cometary materials, both in the strength- and in the gravity-dominated regime. [Copyright &y& Elsevier]

Comets and small Kuiper belt objects are considered to be among the most primitive objects in the solar system as comets like C/1995 O1 Hale-Bopp are rich in highly volatile ices like CO. It has been suggested that early in the solar system evolution the precursors of both groups, the so-called icy planetesimals, were modified by both short-lived radiogenic heating and collisional heating. Here we employ 2D finite-difference numerical models to study the internal thermal evolution of these objects, where we vary formation time, radius and rock-to-ice mass fraction. Additionally we perform 3D SPH collision models with different impact parameters, thus considering both cratering and catastrophic disruption events. Combining the results of both numerical models we estimate under which conditions highly volatile ices like CO, CO 2 and NH 3 can be retained inside present-day comets and Kuiper belt objects. Our results indicate that for present-day objects derived from the largest post-collision remnant the internal thermal evolution controls the amount of remaining highly volatile ices, while for the objects formed from unbound post-collision material the impact heating is dominant. Finally we apply our results to present-day comets and Kuiper belt objects like 67P/Churyumov-Gerasimenko, C/1995 O1 Hale-Bopp and (486958) Arrokoth. [ABSTRACT FROM AUTHOR]

The current properties of small bodies provide important clues to their origin and history. However, how much small bodies were processed by past collisions and to what extent they retain a record of processes that took place during the formation and early evolution of the Solar System is still poorly understood. Here we study the degree of collisional heating and compaction by analysing the large set of previous simulations of small body break-ups by Jutzi et al. (2019), which used porous targets of 50–400 km in diameter and investigated a large range of impact velocities, angles as well as energies. We find that the degree of impact processing is generally larger than found in previous studies which considered smaller objects (e.g. Jutzi et al., 2017; Schwartz et al., 2018). However, there is a clear dichotomy in terms of impact processing: the escaping material always experiences stronger heating than the material bound to the largest remnant. Assuming they originate from the same parent body, some of the observed differences between the recently visited asteroids Ryugu and Bennu may be explained by a different location of the material eventually forming these asteroids in the original parent body. Our results also provide constraints on the initial size of cometary nuclei. • We study the degree of heating and compaction in collisional break-ups of small bodies (50 - 400 km) • We find that the degree of impact processing is generally larger than found in previous studies of smaller objects (1- 10 km) • There is a clear dichotomy in terms of impact processing of the escaping material vs. the gravitationally bound material • Our results provide constraints for the origin and initial properties of small bodies [ABSTRACT FROM AUTHOR]

Binary asteroid formation is a highly complex process, which has been highlighted with recent observations of satellites with unexpected shapes, such as the oblate Dimorphos by the NASA DART mission and the contact binary Selam by NASA's Lucy mission. There is no clear consensus on which dynamical mechanisms determine the final shape of these objects. In turn, we explore a formation pathway where spin-up and rotational failure of a rubble pile asteroid lead to mass-shedding and a wide circumasteroidal debris disk in which the satellite forms. Using a combination of smooth-particle hydrodynamical and N -body simulations, we study the dynamical evolution in detail. We find that a debris disk containing matter corresponding to a few percent of the primary asteroid mass extending beyond the fluid Roche limit can consistently form both oblate and bilobate satellites via a series of tidal encounters with the primary body and mergers with other gravitational aggregates. Principally, satellites end up prolate (elongated) and on synchronous orbits, accreting mainly in a radial direction while tides from the primary asteroid keep the shape intact. However, close encounters and mergers can break the orbital state, leading to orbital migration and deformation. Satellite–satellite impacts occurring in this regime have lower impact velocities than merger-driven moon formation in e.g. planetary rings, leading to soft impacts between differently sized, non-spherical bodies. The resulting post-merger shape of the satellite is highly dependent on the impact geometry. Only moons having experienced a prior mild or catastrophic tidal disruption during a close encounter with the primary asteroid can become oblate spheroids, which is consistent with the predominantly prolate observed population of binary asteroid satellites. • Wide debris disks facilitate formation of atypically shaped asteroid satellites. • In wide debris disks formation is driven via mergers. • Oblate shapes can only be reached via a combination of tidal disruptions and mergers. • Satellite mergers in orbit around an asteroid have low impact velocities. • Impact geometry and initial shapes highly affect the final shape. [ABSTRACT FROM AUTHOR]

The origin of the Martian Crustal Dichotomy is a long-standing mystery since its discovery in the Mariner 9 era. Among various proposed hypotheses, a single giant impact origin (i.e. the Borealis impact) is the most well known, and the most studied. However, studies that include realistic impact models often adapt a simplified geological and geophysical model for predicting the final crustal distribution, while long-term mantle convection studies have mostly employed an over-simplified parametrization of the impact. Here we use a coupled SPH-thermochemical approach to first simulate an impact event, and then use the result of this realistic model as the initial condition for the long-term mantle convection model. We demonstrate that a giant impact collision results in a mantle-deep magma pond, which upon crystallization leads to more crust production on the impacted hemisphere. In other words, an impact-origin of Mars's southern highlands requires the giant impact to occur in the southern hemisphere. We find that both the impact scenario and the mantle properties affect the geometry of the impact-induced crust ("highlands") and the subsequent state of the interior, and that the formation of "highlands" extends beyond the initial magma pond. We show that a near head-on (15° from the normal) impact event with impactor radius intermediate between 500–750 km, together with a mantle viscosity of 1020 Pa s, can best reproduce the southern highlands of Mars with a geometry similar to that of present-day observations. • The formation mechanism of the Martian Crustal Dichotomy is investigated. • Coupled impact and interior evolution models are performed. • Crustal thickening of the impacted hemisphere forms the Martian southern highland. • An impactor radius between 500-750 km and a mantle viscosity of 1020 Pa s is favored. [ABSTRACT FROM AUTHOR]

We examine the relevance of tidal heating for large Trans-Neptunian Objects, with a focus on its potential to melt and maintain layers of subsurface liquid water. Depending on their past orbital evolution, tidal heating may be an important part of the heat budget for a number of discovered and hypothetical TNO systems and may enable formation of, and increased access to, subsurface liquid water. Tidal heating induced by the process of despinning is found to be particularly able to compete with heating due to radionuclide decay in a number of different scenarios. In cases where radiogenic heating alone may establish subsurface conditions for liquid water, we focus on the extent by which tidal activity lifts the depth of such conditions closer to the surface. While it is common for strong tidal heating and long lived tides to be mutually exclusive, we find this is not always the case, and highlight when these two traits occur together. We find cases where TNO systems experience tidal heating that is a significant proportion of, or greater than radiogenic heating for periods ranging from100′s of millions to a billion years. For subsurface oceans that contain a small antifreeze component, tidal heating due to very high initial spin states may enable liquid water to be preserved right up to the present day. Of particular interest is the Eris-Dysnomia system, which in those cases may exhibit extant cryovolcanism. [ABSTRACT FROM AUTHOR]

We have developed a method to investigate the whole mass–velocity distribution of impact fragments generated by catastrophic disruption of simulated planetesimals. Flash X-ray radiography including 12 iron particles for tracers was used to visualize the interior of the target, and the velocity distribution of the whole target was estimated by using the velocities of the tracers measured by X-ray images. High-velocity impact experiments in the laboratory and numerical simulations were conducted for four types of targets simulating rocky and frozen planetesimals at various specific energies, Q. These targets consisted of frozen clays with three different water contents ranging from 25 to 45 wt% and porous gypsum with a porosity of 51%. The shattering strength, Q S *, and the mass–velocity distribution (MVD) were studied for these targets. The Q S * of the frozen clays varied by a factor of 3–4 times, depending on the water content, and the Q S * for porous gypsum was almost the same as that for the frozen clays with lower water contents. The numerical impact simulations led to slightly different Q S * and MVD values for the frozen clay targets, possibly because of the partly ductile behavior of these samples. The MVDs resulting from the porous gypsum targets were well reproduced in the simulations. The cumulative mass of fragments with an ejection velocity slower than a specific velocity was examined to introduce a median velocity, v ⁎, charactering the mass–velocity distribution. The v ⁎ is defined as the velocity at which the cumulative mass corresponds to a half of the original target mass in the distribution. The v ⁎ values of the frozen clays were described by the empirical equation v ∗ = ε Q γ with almost the same ε and γ , irrespective of the water content, but the v ⁎ of porous gypsum was about 1/3 that of the frozen clays. These experimental results were well reproduced by the numerical simulations for both frozen clays and porous gypsum targets. The dispersion strength, Q D *, could be derived by comparing v ⁎ with the escape velocity, v esc , of a target body with an effective mass, M , and radius, R. From this, a semi-theoretical equation showing the dispersion strength was derived: Q D ∗ = 1 ε 2 GM R 1 / 2 1 / γ. Numerical simulations of catastrophic disruptions including self-gravity were conducted to directly determine Q D * at large scale. These calculations showed that the effective mass of the target body, which is used in the computation of v esc = v ⁎, should be a half of the original target mass, M = M target /2. Our results suggest that this approach for computing the semi-theoretical dispersion strength is suitable for bodies larger than ~10 km. • A flash X-ray radiography visualized fragments caused by impact disruption. • Mass-velocity distribution of impact fragments were experimentally obtained. • v ⁎ is defined as the velocity below which a half of fragments has, and determined. • Dispersion strength Q D ⁎ is estimated by v ⁎ and escape velocity of planetesimals. • Q D ⁎ estimated experimentally is consistent with that determined numerically. [ABSTRACT FROM AUTHOR]