Numerical modeling of oblique hypervelocity impacts on strong ductile targets

Abstract– The majority of meteorite impacts occur at oblique incidence angles. However, many of the effects of obliquity on impact crater size and morphology are poorly understood. Laboratory experiments and numerical models have shown that crater size decreases with impact angle, the along‐range crater profile becomes asymmetric at low incidence angles, and below a certain threshold angle the crater planform becomes elliptical. Experimental results at approximately constant impact velocity suggest that the elliptical threshold angle depends on target material properties. Herein, we test the hypothesis that the threshold for oblique crater asymmetry depends on target material strength. Three‐dimensional numerical modeling offers a unique opportunity to study the individual effects of both impact angle and target strength; however, a systematic study of these two parameters has not previously been performed. In this work, the three‐dimensional shock physics code iSALE‐3D is validated against laboratory experiments of impacts into a strong, ductile target material. Digital elevation models of craters formed in laboratory experiments were created from stereo pairs of scanning electron microscope images, allowing the size and morphology to be directly compared with the iSALE‐3D craters. The simulated craters show excellent agreement with both the crater size and morphology of the laboratory experiments. iSALE‐3D is also used to investigate the effect of target strength on oblique incidence impact cratering. We find that the elliptical threshold angle decreases with decreasing target strength, and hence with increasing cratering efficiency. Our simulations of impacts on ductile targets also support the prediction from Chapman and McKinnon (1986) that cratering efficiency depends on only the vertical component of the velocity vector.     

Geophysical consequences of planetary-scale impacts into a Mars-like planet

All planetary bodies with old surfaces exhibit planetary-scale impact craters: vast scars caused by the large impacts at the end of Solar System accretion or the late heavy bombardment. Here we investigate the geophysical consequences of planetary-scale impacts into a Mars-like planet, by simulating the events using a smoothed particle hydrodynamics (SPH) model. Our simulations probe impact energies over two orders of magnitude (2 × 10²⁷-6 × 10²⁹ J), impact velocities from the planet’s escape velocity to twice Mars’ orbital velocity (6-50 km/s), and impact angles from head-on to highly oblique (0-75°). The simulation results confirm that for planetary-scale impacts, surface curvature, radial gravity, the large relative size of the impactor to the planet, and the greater penetration of the impactor, contribute to significant differences in the geophysical expression compared to small craters, which can effectively be treated as acting in a half-space. The results show that the excavated crustal cavity size and the total melt production scale similarly for both small and planetary-scale impacts as a function of impact energy. However, in planetary-scale impacts a significant fraction of the melt is sequestered at depth and thus does not contribute to resetting the planetary surface; complete surface resetting is likely only in the most energetic (6 × 10²⁹ J), slow, and head-on impacts simulated. A crater rim is not present for planetary-scale impacts with energies >10²⁹ J and angles ?45°, but rather the ejecta is more uniformly distributed over the planetary surface. Antipodal crustal removal and melting is present for energetic (>10²⁹ J), fast (>6 km/s), and low angle (?45°) impacts. The most massive impactors (with both high impact energy and low velocity) contribute sufficient angular momentum to increase the rotation period of the Mars-sized target to about a day. Impact velocities of >20 km/s result in net mass erosion from the target, for all simulated energies and angles. The hypothesized impact origin of planetary structures may be tested by the presence and distribution of the geochemically-distinct impactor material.     

Nonuniform cratering of the Moon and a revised crater chronology of the inner Solar System

We model the cratering of the Moon and terrestrial planets from the present knowledge of the orbital and size distribution of asteroids and comets in the inner Solar System, in order to refine the crater chronology method. Impact occurrences, locations, velocities and incidence angles are calculated semi-analytically, and scaling laws are used to convert impactor sizes into crater sizes. Our approach is generalizable to other moons or planets. The lunar cratering rate varies with both latitude and longitude: with respect to the global average, it is about 25% lower at (±65°N, 90°E) and larger by the same amount at the apex of motion (0°N, 90°W) for the present Earth-Moon separation. The measured size-frequency distributions of lunar craters are reconciled with the observed population of near-Earth objects under the assumption that craters smaller than a few kilometers in diameter form in a porous megaregolith. Varying depths of this megaregolith between the mare and highlands is a plausible partial explanation for differences in previously reported measured size-frequency distributions. We give a revised analytical relationship between the number of craters and the age of a lunar surface. For the inner planets, expected size-frequency crater distributions are calculated that account for differences in impact conditions, and the age of a few key geologic units is given. We estimate the Orientale and Caloris basins to be 3.73 Ga old, and the surface of Venus to be 240 Ma old. The terrestrial cratering record is consistent with the revised chronology and a constant impact rate over the last 400 Ma. Better knowledge of the orbital dynamics, crater scaling laws and megaregolith properties are needed to confidently assess the net uncertainty of the model ages that result from the combination of numerous steps, from the observation of asteroids to the formation of craters. Our model may be inaccurate for periods prior to 3.5 Ga because of a different impactor population, or for craters smaller than a few kilometers on Mars and Mercury, due to the presence of subsurface ice and to the abundance of large secondaries, respectively. Standard parameter values allow for the first time to naturally reproduce both the size distribution and absolute number of lunar craters up to 3.5 Ga ago, and give self-consistent estimates of the planetary cratering rates relative to the Moon.     

Non-intrusive measurements of crater growth

An experimental technique to measure crater growth is presented whereby a high speed video captures profiles of a crater forming after impact obtained using a vertical laser sheet centered on the impact point. Unlike previous so called “quarter-space experiments,” where projectiles were launched along a transparent Plexiglas sheet so that growth of half a crater could be viewed, the use of the laser sheet permits viewing changes in crater shape without any physical interference to the cratering process. This technique indicates that for low velocity impacts (<300 m/s) into 220 μm glass beads that are without cohesion and where the projectile is not disrupted, craters initially grow somewhat proportionally, but that later their depths remain essentially constant while their diameters continue to expand. In addition, these experiments indicate that as the impact velocity increases, the rate of growth and the transient depth to diameter ratio at the end of ejecta excavation decreases. These last two observations are probably due to the large time of penetration of the projectile, which becomes a significant fraction of the time of crater formation. This is contrary to the expectations for the scaling rules, which assumes a point source. Very high curtain angles (>45°) are also seen, and could be due to the low friction angle of the target. Significant crater modification, which is rarely seen in “quarter-space experiments,” is also observed and appears to be controlled by the dynamic angle of repose of the target. These latter observations indicate that differences in target friction angles may need to be considered when determining near rim ejecta-mass distributions and large-scale crater modification processes on the planets.     

High-resolution simulations of the impacts of asteroids into the venusian atmosphere III: further 3D models

We report on high-resolution three-dimensional calculations of oblique impacts into planetary atmospheres, specifically the atmosphere of Venus, extending the results of Korycansky et al. (2000, Icarus 146, 387-403; 2002, Icarus 157, 1-23). We have made calculations for impacts at 0°, 45°, and 60° from the vertical, different impactor velocities (10, 20, and 40 km s⁻¹), and different impactor masses and orientations. We present results for porous impactors using a simple model of porosity. We have investigated the sensitivity to initial conditions of the calculations [as a follow-up to the results found in Korycansky et al. (2002)] and resolution effects. For use in cratering calculations, we fit simple functions to the numerical results for mass and momentum that penetrate to a given altitude (column mass) and investigate the behavior of the fit coefficients as functions of impactor parameters such as mass, velocity, and impact angle. Generally speaking, the mass and momentum (and hence resulting crater diameters) depend primarily on impactor mass and mass of atmosphere encountered and weakly or not at all on other parameters such as impactor velocity, impact angle, or porosity. The column mass to which the last portion of the impactor penetrates is approximately equal to the mass of impactor at the top of the atmosphere before the impact takes place. Finally, we present the beginnings of a simplified but physically based model for the impactor and its fragments to reproduce the mass and momentum fluxes as a function of height during the impact.     

Near-Earth object velocity distributions and consequences for the Chicxulub impactor

An Öpik-based geometric algorithm is used to compute impact probabilities and velocity distributions for various near-Earth object (NEO) populations. The resulting crater size distributions for the Earth and Moon are calculated by combining these distributions with assumed NEO size distributions and a selection of crater scaling laws. This crater probability distribution indicates that the largest craters on both the Earth and the Moon are dominated by comets. However, from a calculation of the fractional probabilities of iridium deposition, and the velocity distributions at impact of each NEO population, the only realistic possibilities for the Chicxulub impactor are a short-period comet (possibly inactive) or a near-Earth asteroid. For these classes of object, sufficiently large impacts have mean intervals of 100 and 300 Myr respectively, slightly favouring the cometary hypothesis.     

The shape and appearance of craters formed by oblique impact on the Moon and Venus

Abstract— We surveyed the impact crater populations of Venus and the Moon, dry targets with and without an atmosphere, to characterize how the 3‐dimensional shape of a crater and the appearance of the ejecta blanket varies with impact angle. An empirical estimate of the impact angle below which particular phenomena occur was inferred from the cumulative percentage of impact craters exhibiting different traits. The results of the surveys were mostly consistent with predictions from experimental work. Assuming a sin²θ dependence for the cumulative fraction of craters forming below angle θ, on the Moon, the following transitions occur: >?45 degrees, the ejecta blanket becomes asymmetric; >?25 degrees, a forbidden zone develops in the uprange portion of the ejecta blanket, and the crater rim is depressed in that direction; >?15 degrees, the rim becomes saddle‐shaped; >?10 degrees, the rim becomes elongated in the direction of impact and the ejecta forms a “butterfly” pattern. On Venus, the atmosphere causes asymmetries in the ejecta blanket to occur at higher impact angles. The transitions on Venus are: >?55 degrees, the ejecta becomes heavily concentrated downrange; >?40 degrees, a notch in the ejecta that extends to the rim appears, and as impact angle decreases, the notch develops into a larger forbidden zone; >?10 degrees, a fly‐wing pattern develops, where material is ejected in the crossrange direction but gets swept downrange. No relationship between location or shape of the central structure and impact angle was observed on either planet. No uprange steepening and no variation in internal slope or crater depth could be associated with impact angle on the Moon. For both planets, as the impact angle decreases from vertical, first the uprange and then the downrange rim decreases in elevation, while the remainder of the rim stays at a constant elevation. For craters on Venus >?15 km in diameter, a variety of crater shapes are observed because meteoroid fragment dispersal is a significant fraction of crater diameter. The longer path length for oblique impacts causes a correlation of clustered impact effects with oblique impact effects. One consequence of this correlation is a shallowing of the crater with decreasing impact angle for small craters.     

Micron‐scale hypervelocity impact craters: Dependence of crater ellipticity and rim morphology on impact trajectory,projectile size,velocity, and shape

The interstellar collector on NASA's Stardust mission captured many particles from sources other than the interstellar dust stream. Impact trajectory may provide a means of discriminating between these different sources, and thus identifying/eliminating candidate interstellar particles. The collector's aerogel preserved a clear record of particle impact trajectory from the inclination and direction of the resultant tracks. However, the collector also contained aluminum foils and, although impact crater studies to date suggest only the most inclined impacts (>45° from normal) produce crater morphologies that indicate trajectory (i.e., distinctly elliptical), these studies have been restricted to much larger (mm and above) scales than are relevant for Stardust (μm). It is unknown how oblique impact crater morphology varies as a function of length scale, and therefore how well Stardust craters preserve details of impactor trajectory. Here, we present data from a series of impact experiments, together with complementary hydrocode modeling, that examine how crater morphology changes with impact angles for different‐sized projectiles. We find that, for our smallest spherical projectiles (2 μm diameter), the ellipticity and rim morphology provide evidence of their inclined trajectory from as little as 15° from normal incidence. This is most likely a result of strain rate hardening in the target metal. Further experiments and models find that variation in velocity and impactor shape complicate these trends, but that rim morphology remains useful in determining impact direction (where the angle of impact is >20° from normal) and may help identify candidate interstellar particle craters on the Stardust collector.     

The effect of the oceans on the terrestrial crater size‐frequency distribution: Insight from numerical modeling

Abstract— On Earth, oceanic impacts are twice as likely to occur as continental impacts, yet the effect of the oceans has not been previously considered when estimating the terrestrial crater size‐frequency distribution. Despite recent progress in understanding the qualitative and quantitative effect of a water layer on the impact process through novel laboratory experiments, detailed numerical modeling, and interpretation of geological and geophysical data, no definitive relationship between impactor properties, water depth, and final crater diameter exists. In this paper, we determine the relationship between final (and transient) crater diameter and the ratio of water depth to impactor diameter using the results of numerical impact models. This relationship applies for normal incidence impacts of stoney asteroids into water‐covered, crystalline oceanic crust at a velocity of 15 km s^?1. We use these relationships to construct the first estimates of terrestrial crater size‐frequency distributions (over the last 100 million years) that take into account the depth‐area distribution of oceans on Earth. We find that the oceans reduce the number of craters smaller than 1 km in diameter by about two‐thirds, the number of craters ?30 km in diameter by about one‐third, and that for craters larger than ?100 km in diameter, the oceans have little effect. Above a diameter of ?12 km, more craters occur on the ocean floor than on land; below this diameter more craters form on land than in the oceans. We also estimate that there have been in the region of 150 impact events in the last 100 million years that formed an impact‐related resurge feature, or disturbance on the seafloor, instead of a crater.     

Secondary craters and ejecta across the solar system: Populations and effects on impact-crater–based chronologies

We review the secondary-crater research over the past decade, and provide new analyses and simulations that are the first to model an accumulation of a combined primary-plus-secondary crater population as discrete cratering events. We develop the secondary populations by using scaling laws to generate ejecta fragments, integrating the trajectories of individual ejecta fragments, noting the location and velocity at impact, and using scaling laws to estimate secondary-crater diameters given the impact conditions. We also explore the relationship between the impactor size–frequency distribution (SFD) and the resulting secondary-crater SFD. Our results from these analyses indicate that the “secondary effect” varies from surface to surface and that no single conclusion applies across the solar system nor at any given moment in time—rather, there is a spectrum of outcomes both spatially and temporally, dependent upon target parameters and the impacting population. Surface gravity and escape speed define the spatial distribution of secondaries. A shallow-sloped impactor SFD will cause proportionally more secondaries than a steeper-sloped SFD. Accounting for the driving factors that define the magnitude and spatial distribution of secondaries is essential to determine the relative population of secondary craters, and their effect on derived surface ages.     

Investigating target versus impactor influences on Martian crater morphology at the simple‐complex transition

Comparing craters of identical diameter on a planet is an empirical method of studying the effects of different target and impactor properties while holding total impact energy nearly constant. We have analyzed the Martian crater population within a narrow diameter range (7 km < crater diameter < 9 km) at the simple‐complex crater transition using three approaches. We looked for correlations of morphology with surface geology using a global crater database and global geologic map. We examined selected regions in detail with high‐resolution images to further understand the relationship between crater morphology and bulk target properties. Finally, we examined craters in close proximity to each other in order to hold target properties constant, so that we could isolate impactor effects on crater morphology. We found a strong correlation between target properties and interior crater morphology, and we found little evidence that impactor properties (other than impact angle) affect crater appearance. Central uplift and wall slumping are enhanced for less consolidated targets. Layered targets affected both the excavation and modification stages of complex crater formation; the resulting craters have pseudoterraces, flat floors, and central pits.     

Initiation of Run-Out Flows on Venus by Oblique Impacts

A hypervelocity oblique impact results in a downrange-moving vapor cloud, a significant fraction of which is derived from the projectile. Since the vapor cloud expands to great extent and becomes very tenuous quickly on a planet with a thin or no atmosphere, it does not leave a well-defined geologic expression. The thick atmosphere of Venus, however, is sufficient to contain such a rapidly expanding vapor cloud. As a result of atmospheric interactions, impact vapor condenses and contributes to run-out flows around craters on Venus. Previous results of both laboratory experiments and simple semi-analytical calculations indicate that an impact-vapor origin can account for the morphology of run-out flows on Venus most consistently. However, the detailed dynamics and geologic record of downrange-moving impact vapor clouds in Venus's atmosphere are not understood quantitatively. To approach these problems, we carried out two-dimensional hydrocode calculations. Parametric studies of these hydrocode calculations yield simple scaling laws for both the total downrange travel distance and the final temperature of impact vapor clouds under conditions on Venus. Under typical impact conditions, impact vapor clouds travel downrange more than a crater radius prior to the completion of crater formation. Furthermore, the scaling law for the total travel distance is compared with observations for the downrange offset of the source regions of run-out flows around oblique craters. The results of this comparison suggest that energy/momentum-partitioning processes other than pure shock coupling may play important roles in hypervelocity impact at planetary scales. The results of hydrocode calculations also indicate that the terminal temperature of the impact vapor is close to the condensation temperatures of silicates, suggesting that two scenarios are possible for expected range of impact conditions: 1. Impact vapor condenses and forms run-out flows. 2. Impact vapor fails to condense and leaves no run-out flows. Consequently, natural variation in impact angle, velocity, and projectile composition may account for partial occurrence of run-out flows around impact craters on Venus.     

Atmospheric impacts, fragmentation, and small craters on Venus

We use high-resolution three-dimensional numerical models of aerodynamically disrupted asteroids to predict the characteristic properties of small impact craters on Venus. We map the mass and kinetic energy of the impactor passing though a plane near the surface for each simulation, and find that the typical result is that mass and energy sort themselves into one to several strongly peaked regions, which we interpret as more-or-less discrete fragments. The fragments are sufficiently well separated as to imply the formation of irregular or multiple craters that are quite similar to those found on Venus. We estimate the diameters of the resulting craters using a scaling law derived from the experiments of Schultz and Gault (1985, J. Geophys. Res. 90 (B5), 3701-3732) of dispersed impactors into targets. We compare the spacings and sizes of our estimated craters with measured diameters tabulated in a Venus crater database (Herrick and Phillips, 1994a, Icarus 111, 387-416; Herrick et al., 1997, in: Venus II, Univ. of Arizona Press, Tucson, AZ, pp. 1015-1046; Herrick, 2003, http://www.lpi.usra.edu/research/vc/vchome.html) and find quite satisfactory agreement, despite the uncertainty in our crater diameter estimates. The comparison of the observed crater characteristics with the numerical results is an after-the-fact test of our model, namely the fluid-dynamical treatment of large impacts, which the model appears to pass successfully.     

Laboratory simulation of the herringbone pattern associated with lunar secondary crater chains

V-shaped ridge components of the herringbone pattern associated with lunar secondary crater chains have been simulated by simultaneous and nearly simultaneous impact of two projectiles near one another. The impact velocities and angles of the projectiles were similar to those of the fragments that produced secondary craters found at various ranges from large lunar craters.Variables found to affect the included angles of the V-shaped ridges are: relative time of impact of the projectiles, impact angle, relative projectile mass, and azimuth angle of the crater chain relative to the projection of the flight line onto the target surface. The functional relationships between the forms of the ridges and many of these variables are similar to those observed for lunar V-shaped ridges.Comparison of the magnitudes of the ridge angles of both laboratory crater pairs and secondary crater chains of the crater Copernicus implies that material was ejected from Copernicus at angles in excess of 60°, measured from the normal, to form many of Copernicus' satellitic craters. Moreover, other independent calculations presented indicate that many of the fragments that produced secondary craters also ricocheted to produce tertiary craters.Application of the study to identification of isolated secondary craters and to the determination of the origin of large lunar craters is discussed.     

Impacts onto H2O ice: Scaling laws for melting, vaporization, excavation, and final crater size

Shock-induced melting and vaporization of H₂O ice during planetary impact events are widespread phenomena. Here, we investigate the mass of shock-produced liquid water remaining within impact craters for the wide range of impact conditions and target properties encountered in the Solar System. Using the CTH shock physics code and the new 5-phase model equation of state for H₂O, we calculate the shock pressure field generated by an impact and fit scaling laws for melting and vaporization as a function of projectile mass, impact velocity, impact angle, initial temperature, and porosity. Melt production nearly scales with impact energy, and natural variations in impact parameters result in only a factor of two change in the predicted mass of melt. A fit to the π-scaling law for the transient cavity and transient-to-final crater diameter scaling are determined from recent simulations of the entire cratering process in ice. Combining melt production with π-scaling and the modified Maxwell Z-model for excavation, less than half of the melt is ejected during formation of the transient crater. For impact energies less than about 2 × 10²⁰ J and impact velocities less than about 5 km s⁻¹, the remaining melt lines the final crater floor. However, for larger impact energies and higher impact velocities, the phenomenon of discontinuous excavation in H₂O ice concentrates the impact melt into a small plug in the center of the crater floor.     

Modeling the Ries‐Steinheim impact event and the formation of the moldavite strewn field

Abstract— Using detailed geological, petrographic, geochemical, and geographical constraints we have performed numerical modeling studies that relate the Steinheim crater (apparent diameter D_a = 3.8 km), the Ries crater (D_a = 24 km) in southern Germany, and the moldavite (tektite) strewn field in Bohemia and Moravia (Czech Republic), Lusatia (East Germany), and Lower Austria. The moldavite strewn field extends from ～200 to 450 km from the center of the Ries to the east‐northeast forming a fan with an angle of ～57°. An oblique impact of a binary asteroid from a west‐southwest direction appears to explain the locations of the craters and the formation and distribution of the moldavites. The impactor must have been a binary asteroid with two widely separated components (some 1.5 and 0.15 km in diameter, respectively). We carried out a series of three‐dimensional hydrocode simulations of a Ries‐type impact. The results confirm previous results suggesting that impacts around 30–50° (from the horizontal) are the most favorable angles for near‐surface melting, and, consequently for the formation of tektites. Finally, modeling of the motion of impact‐produced tektite particles through the atmosphere produces, in the downrange direction, a narrow‐angle distribution of the moldavites tektites in a fan like field with an angle of ～75°. An additional result of modeling the motion of melt inside and outside the crater is the preferred flow of melt from the main melt zone of the crystalline basement downrange towards the east‐northeast rim. This explains perfectly the occurrence of coherent impact melt bodies (some tens of meters in size) in a restricted zone of the downrange rim of the Ries crater. The origin of these melt bodies, which represent chemically a mixture of crystalline basement rocks similar to the main melt mass contained (as melt particles <0.5 m in size) in the suevite, do not occur at any other portion of the Ries crater rim and remained enigmatic until now. Although the calculated distribution of moldavites still deviates to some degree from the known distribution, our results represent an important step toward a better understanding of the origin and distribution of the high‐velocity surface melts and the low‐velocity, deep‐seated melt resulting from an oblique impact on a stratified target.     

On the Shape of an Impact Crater at High-Speed Impacts

This paper presents the results of an experimental study on the geometric parameters of craters that originated at the impact interaction of polyethylene projectiles with a massive organic-glass target. The impactor speed ranged from 2.7 to 6.28 km/s. The shapes of the resulting craters are determined. Using statistical analysis and scaling theory, we obtained the dependence of the diameter and depth of the crater on defining parameters. The critical impact energy resulting in the catastrophic breakup of the target is estimated.     

Modeling crater populations on Venus and Titan

We describe a model for crater populations on planets and satellites with dense atmospheres, like those of Venus and Titan. The model takes into account ablation (or mass shedding), pancaking, and fragmentation. Fragmentation is assumed to occur due to the hydrodynamic instabilities promoted by the impactors’ deceleration in the atmosphere. Fragments that survive to hit the ground make craters or groups thereof. Crater sizes are estimated using standard laws in the gravity regime, modified to take into account impactor disruption. We use Monte Carlo methods to pick parameters from appropriate distributions of impactor mass, zenith angle, and velocity. Good fits to the Venus crater populations (including multiple crater fields) can be found with reasonable values of model parameters. An important aspect of the model is that it reproduces the dearth of small craters on Venus: this is due to a cutoff on crater formation we impose, when the expected crater would be smaller than the (dispersed) object that would make it. Hydrodynamic effects alone (ablation, pancaking, fragmentation) due to the passage of impactors through the atmosphere are insufficient to explain the lack of small craters. In our favored model, the observed number of craters (940) is produced by ∼5500 impactors with masses , yielding an age of (1-σ uncertainty) for the venusian surface. This figure does not take into account any uncertainties in crater scaling and impactor population characteristics, which probably increase the uncertainty to a factor of two in age.We apply the model with the same parameter values to Titan to predict crater populations under differing assumptions of impactor populations that reflect present conditions. We assume that the impactors (comets) are made of 50% porous ice. Predicted crater production rates are ≈190 craters . The smallest craters on Titan are predicted to be in diameter, and ≈5 crater fields are expected. If the impactors are composed of solid ice (density ), crater production rates increase by ≈70% and the smallest crater is predicted to be in diameter. We give cratering rates for denser comets and atmospheres 0.1 and 10 times as thick as Titan's current atmosphere. We also explicitly address leading-trailing hemisphere asymmetries that might be seen if Titan's rotation rate were strictly synchronous over astronomical timescales: if that is the case, the ratio of crater production on the leading hemisphere to that on the trailing hemisphere is ≈4:1.     

Numerical modeling of asteroid survivability and possible scenarios for the Morokweng crater‐forming impact

The fate of the impactor is an important aspect of the impact‐cratering process. Defining impactor material as surviving if it remains solid (i.e., does not melt or vaporize) during crater formation, previous numerical modeling and experiments have shown that survivability decreases with increasing impact velocity, impact angle (with respect to the horizontal), and target density. Here, we show that in addition to these, impactor survivability depends on the porosity and shape of the impactor. Increasing impactor porosity decreases impactor survivability, while prolate‐shaped (polar axis > equatorial axis) impactors survive impact more so than spherical and oblate‐shaped (polar axis < equatorial axis) impactors. These results are used to produce a relatively simple equation, which can be used to estimate the impactor fraction shocked to a given pressure as a function of these parameters. By applying our findings to the Morokweng crater‐forming impact, we suggest impact scenarios that explain the high meteoritic content and presence of unmolten fossil meteorites within the Morokweng crater. In addition to previous suggestions of a low‐velocity and/or high‐angled impact, this work suggests that an elongated and/or low porosity impactor may also help explain the anomalously high survivability of the Morokweng impactor.     

Azimuthal impact directions from oblique impact crater morphology

Planetary impact craters have a high degree of radial symmetry. This hampers efforts to identify the azimuthal impact direction for most craters – the radially symmetric component of an impact crater swamps any asymmetries that may be present. We demonstrate how the asymmetric component can be isolated and the direction of the asymmetries quantified using a two-dimensional eigenfunction expansion over a circular domain. The complex coefficients of expansion describe the magnitude and phase (angular alignment) of each term. From the analysis of hypervelocity impact craters formed in the laboratory, with impact angles ranging from 0° to 50° from the surface normal, we show that asymmetries which reveal the impact direction are still present at just 10° from the surface normal, and that the phase of one complex coefficient of expansion, c ₁₁, indicates the impact direction. Analysis of the lunar crater Hadley shows bilateral symmetry in the radially asymmetric component, which may be due to oblique impact. The 31-km lunar ray crater Kepler has morphological features that indicate the azimuthal impact direction. Coefficient c ₁₁ gives an azimuthal impact direction similar to that expected from the morphology, although post-impact gravitational collapse and slumping obscure the result to some degree. Ray craters may provide a means of testing the method for smaller 'simple' craters when data are available.