Crater morphology in sandstone targets: The MEMIN impact parameter study

Abstract– Hypervelocity (2.5–7.8 km s^?1) impact experiments into sandstone were carried out to investigate the influence of projectile velocity and mass, target pore space saturation, target‐projectile density contrast, and target layer orientation on crater size and shape. Crater size increases with increasing projectile velocity and mass as well as with increasing target pore space saturation. Craters in water‐saturated porous targets are generally shallower and larger in volume and in diameter than craters from equivalent impacts into dry porous sandstone. Morphometric analyses of the resultant craters, 5–40 cm in diameter, reveal features that are characteristic of all of our experimental craters regardless of impact conditions (I) a large central depression within a fragile, light‐colored central part, and (II) an outer spallation zone with areas of incipient spallation. Two different mechanical processes, grain fragmentation and intergranular tensile fracturing, are recorded within these crater morphologies. Zone (I) approximates the shape of the transient crater formed by material compression, displacement, comminution, and excavation flow, whereas (II) is the result of intergranular tensile fracturing and spallation. The transient crater dimensions are reconstructed by fitting quadric parabolas to crater profiles from digital elevation models. The dimensions of this transient and of the final crater show the same trends: both increase in volume with increasing impact energy, and with increasing water saturation of the target pore space. The relative size of the transient crater (in percent of the final crater volume) decreases with increasing projectile mass and velocity, signifying a greater contribution of spallation on the final crater size when projectile mass and velocity are increased.     

Martian subsurface properties and crater formation processes inferred from fresh impact crater geometries

Abstract— The geometry of simple impact craters reflects the properties of the target materials, and the diverse range of fluidized morphologies observed in Martian ejecta blankets are controlled by the near‐surface composition and the climate at the time of impact. Using the Mars Orbiter Laser Altimeter (MOLA) data set, quantitative information about the strength of the upper crust and the dynamics of Martian ejecta blankets may be derived from crater geometry measurements. Here, we present the results from geometrical measurements of fresh craters 3–50 km in rim diameter in selected highland (Lunae and Solis Plana) and lowland (Acidalia, Isidis, and Utopia Planitiae) terrains. We find large, resolved differences between the geometrical properties of the freshest highland and lowland craters. Simple lowland craters are 1.5–2.0 times deeper (≥5s?_o difference) with >50% larger cavities (≥2s?_o) compared to highland craters of the same diameter. Rim heights and the volume of material above the preimpact surface are slightly greater in the lowlands over most of the size range studied. The different shapes of simple highland and lowland craters indicate that the upper ?6.5 km of the lowland study regions are significantly stronger than the upper crust of the highland plateaus. Lowland craters collapse to final volumes of 45–70% of their transient cavity volumes, while highland craters preserve only 25–50%. The effective yield strength of the upper crust in the lowland regions falls in the range of competent rock, approximately 9–12 MPa, and the highland plateaus may be weaker by a factor of 2 or more, consistent with heavily fractured Noachian layered deposits. The measured volumes of continuous ejecta blankets and uplifted surface materials exceed the predictions from standard crater scaling relationships and Maxwell's Z model of crater excavation by a factor of 3. The excess volume of fluidized ejecta blankets on Mars cannot be explained by concentration of ejecta through nonballistic emplacement processes and/or bulking. The observations require a modification of the scaling laws and are well fit using a scaling factor of ?1.4 between the transient crater surface diameter to the final crater rim diameter and excavation flow originating from one projectile diameter depth with Z = 2.7. The refined excavation model provides the first observationally constrained set of initial parameters for study of the formation of fluidized ejecta blankets on Mars.     

Mid‐sized complex crater formation in mixed crystalline‐sedimentary targets: Insight from modeling and observation

Abstract— Large impact crater formation is an important geologic process that is not fully understood. The current paradigm for impact crater formation is based on models and observations of impacts in homogeneous targets. Real targets are rarely uniform; for example, the majority of Earth's surface is covered by sedimentary rocks and/or a water layer. The ubiquity of layering across solar system bodies makes it important to understand the effect target properties have on the cratering process. To advance understanding of the mechanics of crater collapse, and the effect of variations in target properties on crater formation, the first “Bridging the Gap” workshop recommended that geological observation and numerical modeling focussed on mid‐sized (15–30 km diameter) craters on Earth. These are large enough to be complex; small enough to be mapped, surveyed and modelled at high resolution; and numerous enough for the effects of target properties to be potentially disentangled from the effects of other variables. In this paper, we compare observations and numerical models of three 18–26 km diameter craters formed in different target lithology: Ries, Germany; Haughton, Canada; and El'gygytgyn, Russia. Based on the first‐order assumption that the impact energy was the same in all three impacts we performed numerical simulations of each crater to construct a simple quantitative model for mid‐sized complex crater formation in a subaerial, mixed crystalline‐sedimentary target. We compared our results with interpreted geological profiles of Ries and Haughton, based on detailed new and published geological mapping and published geophysical surveys. Our combined observational and numerical modeling work suggests that the major structural differences between each crater can be explained by the difference in thickness of the pre‐impact sedimentary cover in each case. We conclude that the presence of an inner ring at Ries, and not at Haughton, is because basement rocks that are stronger than the overlying sediments are sufficiently close to the surface that they are uplifted and overturned during excavation and remain as an uplifted ring after modification and post‐impact erosion. For constant impact energy, transient and final crater diameters increase with increasing sediment thickness.     

An empirical model for transient crater growth in granular targets based on direct observations

The present paper describes observations of crater growth up to the time of transient crater formation and presents a new empirical model for transient crater growth as a function of time. Polycarbonate projectiles were impacted vertically into soda-lime glass sphere targets using a single-stage light-gas gun. Using a new technique with a laser sheet illuminating the target [Barnouin-Jha, O.S., Yamamoto, S., Toriumi, T., Sugita, S., Matsui, T., 2007. Non-intrusive measurements of the crater growth. Icarus, 188, 506-521], we measured the temporal change in diameter of crater cavities (diameter growth). The rate of increase in diameter at early times follows a power law relation, but the data at later times (before the end of transient crater formation) deviates from the power law relation. In addition, the power law exponent at early times and the degree of deviation from a power law at later times depend on the target. In order to interpret these features, we proposed to modify Maxwell’s Z-model under the assumption that the strength of the excavation flow field decreases exponentially with time. We also derived a diameter growth model as: d(t)∝[1-exp(-βt)]^γ, where d(t) is the apparent diameter of the crater cavity at time t after impact, and β and γ are constants. We demonstrated that the diameter growth model could represent well the experimental data for various targets with different target material properties, such as porosity or angle of repose. We also investigated the diameter growth for a dry sand target, which has been used to formulate previous scaling relations. The obtained results showed that the dry sand target has larger degree of deviation from a power law, indicating that the target material properties of the dry sand target have a significant effect on diameter growth, especially at later times. This may suggest that the previously reported scaling relations should be reexamined in order to account for the late-stage behavior with the effect of target material properties.     

Cratering and modification of wet‐target craters: Projectile impact experiments and field observations of the Lockne marine‐target crater (Sweden)

Abstract— Marine impacts are one category of crater formation in volatile targets. At target water depths exceeding the diameter of the impactor, the zones of vaporization, melting, and excavation of the standard land‐target cratering model develop partially or entirely in the water column. The part of the crater that has a potential of being preserved (seafloor crater) may to a great extent be formed by material emplacement and excavation processes that are very different from land‐target craters. These processes include a high‐energy, water‐jet‐driven excavation flow. At greater water depths, the difference in strength of the target layers causes a concentric crater to evolve. The crater consists of a wide water cavity with a shallow excavation flow along the seabed surrounding a nested, deeper crater in the basement. The modification of the crater is likewise influenced by the water through its forceful resurge to fill the cavity in the water mass and the seafloor. The resurge flow is strongly erosive and incorporates both ejecta and rip‐up material from the seabed surrounding the excavated crater. A combination of field observations and impact experiments has helped us analyze the processes affecting the zone between the basement crater and the maximum extent of the water cavity. The resurge erosion is facilitated by fragmentation of the upper parts of the solid target caused by a) spallation and b) vibrations from the shallow excavation flow and, subsequently, c) the vertical collapse of the water cavity rim wall. In addition, poorly consolidated and saturated sediments may collapse extensively, possibly aided by a violent expansion of the pore water volume when it turns into a spray during passage of the rarefaction wave. This process may also occur at impacts into water‐saturated targets without an upper layer of seawater present. Our results have implications for impacts on both Earth and Mars, and possibly anywhere in the solar system where volatiles exist/have existed in the upper part of the target.     

Chesapeake Bay impact structure: Morphology,crater fill,and relevance for impact structures on Mars

Abstract— The late Eocene Chesapeake Bay impact structure (CBIS) on the Atlantic margin of Virginia is one of the largest and best‐preserved “wet‐target” craters on Earth. It provides an accessible analog for studying impact processes in layered and wet targets on volatile‐rich planets. The CBIS formed in a layered target of water, weak clastic sediments, and hard crystalline rock. The buried structure consists of a deep, filled central crater, 38 km in width, surrounded by a shallower brim known as the annular trough. The annular trough formed partly by collapse of weak sediments, which expanded the structure to ?85 km in diameter. Such extensive collapse, in addition to excavation processes, can explain the “inverted sombrero” morphology observed at some craters in layered targets. The distribution of crater‐fill materials in the CBIS is related to the morphology. Suevitic breccia, including pre‐resurge fallback deposits, is found in the central crater. Impact‐modified sediments, formed by fluidization and collapse of water‐saturated sand and silt‐clay, occur in the annular trough. Allogenic sediment‐clast breccia, interpreted as ocean‐resurge deposits, overlies the other impactites and covers the entire crater beneath a blanket of postimpact sediments. The formation of chaotic terrains on Mars is attributed to collapse due to the release of volatiles from thick layered deposits. Some flat‐floored rimless depressions with chaotic infill in these terrains are impact craters that expanded by collapse farther than expected for similar‐sized complex craters in solid targets. Studies of crater materials in the CBIS provide insights into processes of crater expansion on Mars and their links to volatiles.     

Structural evidence from shock metamorphism in simple and complex impact craters: Linking observations to theory

Abstract— The structure of Canadian impact craters formed in crystalline rocks is analyzed using shock metamorphism and evidence for movement along shear zones. The analysis is based on an interpretation that, beyond the near field region, shock pressure attenuates down axis as P ? R^?2, in agreement with nuclear test and computed results, and as P ? R^?3 near the surface. In both simple and complex craters, the transient cavity is defined by the limit of fragmentation due to direct and reflected shock waves. The intersection of the transient cavity with hemispheric shock isobars indicates that the transient cavity has a parabolic form. Weakening by dilation during early uplift allows late stage slumping of the walls of simple craters. This is controlled by a spheroidal primary shear of radius r_s ～ 2d_t, where d_t is the depth of the transient crater due to excavation and initial compression. With increasing crater diameter, the size of the transient cavity decreases relative to the shock imprint, suggesting that fragmentation and excavation is limited by progressively earlier collapse of the margins under gravity. Central peak formation in complex craters may be initiated by relaxation of the shock‐compressed central parautochthone, so the primary shear, lubricated by friction melting, meets below the crater floor and drives the continuing upward motion.     

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.     

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.     

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.     

Transitional impact craters on the Moon: Insight into the effect of target lithology on the impact cratering process

We studied a data set of 28 well‐preserved lunar craters in the transitional (simple‐to‐complex) regime with the aim of investigating the underlying cause(s) for morphological differences of these craters in mare versus highland terrains. These transitional craters range from 15 to 42 km in diameter, demonstrating that the transition from simple to complex craters is not abrupt and occurs over a broad diameter range. We examined and measured the following crater attributes: depth (d), diameter (D), floor diameter (D_f), rim height (h), and wall width (w), as well as the number and onset of terraces and rock slides. The number of terraces increases with increasing crater size and, in general, mare craters possess more terraces than highland craters of the same diameter. There are also clear differences in the d/D ratio of mare versus highland craters, with transitional craters in mare targets being noticeably shallower than similarly sized highland craters. We propose that layering in mare targets is a major driver for these differences. Layering provides pre‐existing planes of weakness that facilitate crater collapse, thus explaining the overall shallower depths of mare craters and the onset of crater collapse (i.e., the transition from simple to complex crater morphology) at smaller diameters as compared to highland craters. This suggests that layering and its interplay with target strength and porosity may play a more significant role than previously considered.     

Pluto—The ninth planet's golden year: A meeting held at New Mexico State University,Las Cruces,New Mexico

The simple-to-complex transition for impact craters on Mars occurs at diameters between about 3 and 8 km. Ballistically emplaced ejecta surround primarily those craters that have a simple interior morphology, whereas ejecta displaying features attributable to fluid flow are mostly restricted to complex craters. Size-dependent characteristics of 73 relatively fresh Martian craters, emphasizing the new depth/diameter (d/D) data of D. W. G. Arthur (1980, to be submitted for publication), test two hypotheses for the mode of formation of central peaks in complex craters. In particular, five features appear sequentially with increasing crater size: first flat floors (3–4 km), then central peaks and shallower depths (4–5 km), next scalloped rims (? km), and lastly terraced walls (～8 km). This relative order indicates that a shallow depth of excavation and an unspecified rebound mechanism, not centripetal collapse and deep sliding, have produced central peaks and in turn have facilitated failure of the rim. The mechanism of formation of a shallow crater remains elusive, but probably operates only at the excavation stage of impact. This interpretation is consistent with two separate and complementary lines of evidence. First, field data have documented only shallow subsurface deformation and a shallow transient cavity in complex terrestrial meteorite craters and in certain surface-burst explosion craters; thus the shallow transient cavities of complex craters never were geometrically similar to the deep cavities of simple craters. Second, the average depths of complex craters and the diameters marking the transition from simple to complex craters on Mars and on three other terrestrial planets vary inversely with gravitational acceleration at the planetary surface, g, a variable more important in the excavation of a crater than in any subsequent modification of its geometry. The new interpretation is summarized diagrammatically for complex craters on all planets.     

The effect of target properties on crater morphology: Comparison of central peak craters on the Moon and Ganymede

Abstract— We examine the morphology of central peak craters on the Moon and Ganymede in order to investigate differences in the near‐surface properties of these bodies. We have extracted topographic profiles across craters on Ganymede using Galileo images, and use these data to compile scaling trends. Comparisons between lunar and Ganymede craters show that crater depth, wall slope and amount of central uplift are all affected by material properties. We observe no major differences between similar‐sized craters in the dark and bright terrain of Ganymede, suggesting that dark terrain does not contain enough silicate material to significantly increase the strength of the surface ice. Below crater diameters of ?12 km, central peak craters on Ganymede and simple craters on the Moon have similar rim heights, indicating comparable amounts of rim collapse. This suggests that the formation of central peaks at smaller crater diameters on Ganymede than the Moon is dominated by enhanced central floor uplift rather than rim collapse. Crater wall slope trends are similar on the Moon and Ganymede, indicating that there is a similar trend in material weakening with increasing crater size, and possibly that the mechanism of weakening during impact is analogous in icy and rocky targets. We have run a suite of numerical models to simulate the formation of central peak craters on Ganymede and the Moon. Our modeling shows that the same styles of strength model can be applied to ice and rock, and that the strength model parameters do not differ significantly between materials.     

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.     

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.     

Subsurface deformation of experimental hypervelocity impacts in quartzite and marble targets

Two impact cratering experiments on nonporous rock targets were carried out to determine the influence of target composition on the structural mechanisms of subsurface deformation. Projectiles of 2.5 mm diameter were accelerated to ~5 km s⁻¹ and impacted onto blocks of marble or quartzite. Subsurface deformation was mapped and analyzed on the microscale using thin sections of the bisected craters. Additionally, both experiments were modeled and the calculated strain zones underneath the craters were compared to experimental deformation features. Microanalysis shows that the formation of radial, tensile, and intragranular cracks is a common response of both nonporous materials to impact cratering. In the quartzite target, the subsurface damage is additionally characterized by highly localized deformation along shear bands with intense grain comminution, surrounded by damage zones. In contrast, the marble target shows closely spaced calcite twinning and cleavage activation. Crater diameter and depth as well as the damage lens underneath the crater are unexpectedly smaller in the marble target compared to the quartzite target, which is in contradiction to the marble's much weaker compressive and tensile strengths. However, numerical models result in craters that are similar in size as well as in strain accumulation at the end of transient crater formation, indicating that current models should still be viewed cautiously when compared to experimental details.     

Hydrocode Simulations of Chicxulub Crater Collapse and Peak-Ring Formation

We use hydrocode modeling to investigate dynamic models for the collapse of the Chicxulub impact crater. Our aim is to integrate the results from numerical simulations with kinematic models derived from seismic reflection and wide-angle velocity data to further our understanding of the formation of large impact craters. In our simulations, we model the collapse of a 100-km diameter, bowl-shaped cavity formed in comprehensively fractured crustal material. To facilitate wholesale collapse, we require that the strength of the target be significantly weakened. In the present model, we achieve this using acoustic fluidization, where strong vibrations produced by the expanding shock wave cause extreme pressure fluctuations in the target. At times and positions where the overburden pressure is sufficiently counteracted, the frictional resistance is reduced, enabling the rock debris to flow. Our simulations produce a collapsed crater that contains most of the features that we observe in the seismic data at Chicxulub. In particular, we observe a topographic peak ring, formed as material that is originally part of the central uplift collapses outward and is thrust over the inwardly collapsing transient crater rim. This model for peak-ring generation has not been previously demonstrated by numerical simulations and predicts that the peak ring is composed of deeply derived material and that the stratigraphy within the peak ring is overturned.     

Structural uplift and ejecta thickness of lunar mare craters: New insights into the formation of complex crater rims

We investigate the elevated crater rims of lunar craters. The two main contributors to this elevation are a structural uplift of the preimpact bedrock and the emplacement of ejecta on top of the crater rim. Here, we focus on five lunar complex mare craters with diameters ranging between 16 and 45 km: Bessel, Euler, Kepler, Harpalus, and Bürg. We performed 5281 measurements to calculate precise values for the structural rim uplift and the ejecta thickness at the elevated crater rim. The average structural rim uplift for these five craters amounts to S_RU = 70.6 ± 1.8%, whereas the ejecta thickness amounts to E_T = 29.4 ± 1.8% of the total crater rim elevation. Erosion is capable of modifying the ratio of ejecta thickness to structural rim uplift. However, to minimize the impact of erosion, the five investigated craters are young, pristine craters with mostly preserved ejecta blankets. To quantify how strongly craters were enlarged by crater modification processes, we reconstructed the dimensions of the transient crater. The difference between the transient crater diameter and the final crater diameter can extend up to 11 km. We propose reverse faulting and thrusting at the final crater rim to be one of the main contributing factors of forming the elevated crater rim.     

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.     

Modeling damage and deformation in impact simulations

Abstract— Numerical modeling is a powerful tool for investigating the formation of large impact craters but is one that must be validated with observational evidence. Quantitative analysis of damage and deformation in the target surrounding an impact event provides a promising means of validation for numerical models of terrestrial impact craters, particularly in cases where the final pristine crater morphology is ambiguous or unknown. In this paper, we discuss the aspects of the behavior of brittle materials important for the accurate simulation of damage and deformation surrounding an impact event and the care required to interpret the results. We demonstrate this with an example simulation of an impact into a terrestrial, granite target that produces a 10 km‐diameter transient crater. The results of the simulation are shown in terms of damage (a scalar quantity that reflects the totality of fragmentation) and plastic strain, both total plastic strain (the accumulated amount of permanent shear deformation, regardless of the sense of shear) and net plastic strain (the amount of permanent shear deformation where the sense of shear is accounted for). Damage and plastic strain are both greatest close to the impact site and decline with radial distance. However, the reversal in flow patterns from the downward and outward excavation flow to the inward and upward collapse flow implies that net plastic strains may be significantly lower than total plastic strains. Plastic strain in brittle rocks is very heterogeneous; however, continuum modeling requires that the deformation of the target during an impact event be described in terms of an average strain that applies over a large volume of rock (large compared to the spacing between individual zones of sliding). This paper demonstrates that model predictions of smooth average strain are entirely consistent with an actual strain concentrated along very narrow zones. Furthermore, we suggest that model predictions of total accumulated strain should correlate with observable variations in bulk density and seismic velocity.