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.     

Critical crater diameter and asteroid impact seismology

Abstract— If impact stress reverberation is the primary gradational process on an asteroid at global scales, then the largest undegraded crater records an asteroid's seismological response. The critical crater diameter D_crit is defined as the smallest crater whose formation disrupts all previous craters globally up to its size; it is solved for by combining relationships for crater growth and for stress attenuation. The computation for D_crit gives a simple explanation for the curious observation that small asteroids have only modest undegraded craters, in comparison to their size, whereas large asteroids have giant undegraded craters. D_crit can even exceed the asteroid diameter, in which case all craters are “local” and the asteroid becomes crowded with giant craters. D_crit is the most recent crater to have formed on a blank slate; when it is equated to the measured diameter of the largest undegraded crater on known asteroids, peak particle velocities are found to attenuate with the 1.2–1.3 power of distance—less attenuative than strong shocks, and more characteristic of powerful seismic disturbances. This is to be expected, since global degradation can result from seismic (cm s^?1) particle velocities on small asteroids. Attenuation, as modeled, appears to be higher on asteroids known to be porous, although these are also bodies for which different crater scaling rules might apply.     

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.     

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.     

A study of lunar impact crater size-distributions

Discrepancies in published crater frequency data prompted this study of lunar crater distributions. Effects modifying production size distributions of impact craters such as surface lava flows, blanketing by ejecta, superposition, infilling, and abrasion of craters, mass wasting, and the contribution of secondary and volcanic craters are discussed. The resulting criteria have been applied in the determination of the size distributions of unmodified impact crater populations in selected lunar regions of different ages. The measured cumulative crater frequencies are used to obtain a general calibration size distribution curve by a normalization procedure. It is found that the lunar impact crater size distribution is largely constant in the size range 0.3 km ?D ? 20 km for regions with formation ages between ≈ 3 × 10⁹ yr and ? 4 × 10⁹ yr. A polynomial of 4th degree, valid in the size range 0.8 km ?D ? 20 km, and a polynomial of 7th degree, valid in the size range 0.3 km ?D ? ? 20 km, have been approximated to the logarithm of the cumulative crater frequencyN as a function of the logarithm of crater diameterD. The resulting relationship can be expressed asN ～D ^α(D) where α is a function depending onD. This relationship allows the comparison of crater frequencies in different size ranges. Exponential relationships with constant α, commonly used in the literature, are shown to inadequately approximate the lunar impact crater size distribution. Deviations of measured size distributions from the calibration distribution are strongly suggestive of the existence of processes having modified the primary impact crater population.     

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.     

An analysis of differential impact melt-crater scaling and implications for the terrestrial impact record

Abstract— It has been known for some time that the volume of impact melt (V_m) relative to that of the transient cavity (V_tc) increases with the magnitude of the impact event. This paper investigates the influence that this phenomenon has on the nature of terrestrial impact craters. A model of impact melting is used to estimate the volume of melt produced during the impact of chondritic projectiles into granite targets at velocities of 15, 25, and 50 km S^?1. The dimensions of transient cavities formed under the same impact conditions are calculated from current crater-scaling relationships, which are derived from dimensional analysis of data from cratering experiments. Observed melt volumes at terrestrial craters are collated from the literature and are paired with the transient-cavity diameters (D_tc) of their respective craters; these diameters were determined through an established empirical relationship. The model and observed melt volumes have very similar trends with increasing transient-cavity diameter. This V_m-D_tc relationship is then used to make predictions regarding the nature of the terrestrial cratering record. In particular, with increasing size of the impact event, the depth of melting approaches the depth of the transient cavity. As a consequence, the base of the cavity, which ultimately would appear as an uplifted central structure in a complex crater, will record shock stresses that will increase up to a maximum of partial melting. Examination of the terrestrial record indicates a general trend for higher recorded shock levels in central structures at larger diameters; impact structures in the 100-km size range record partially melted and vesiculated parautochthonous target rocks in their centers. In addition, as the depth of melting approaches a depth equivalent to that attained by the base of the transient cavity, the floor of the transient cavity will have progressively less strength, with the result that cavity modification and uplift will not produce topographic central peaks. Again, the observed terrestrial record is not inconsistent with this prediction, and we offer differential melt scaling as a possible mechanism for the transition from central topographic peaks to rings with increasing crater diameter. Among other implications is the likelihood that impact basins in the 1000-km size range on the early Earth would not have the same multi-ring form as observed on the moon.     

The terrestrial cratering record: I. Current status of observations

The location, size, and principal characteristics of the currently known proven and probable terrestrial impact structures are tabulated. Of the 78 known probable structures, only 3 are Precambrian and the majority are <300 my in age. A survey of the variation in preservation with size and age indicates that, unless protected by sedimentary cover, a structure <20 km in diameter has a recognizable life of <600 my. The depth-diameter relationships of terrestrial structures are similar to lunar craters; however, it is believed that terrestrial craters were always shallower than their lunar counterparts. Complex structures formed in sedimentary targets are shallower than those in crystalline targets, and the transition from simple to complex crater morphology occurs in sedimentary strata at approximately one-half the diameter of the morphology transition in crystalline rocks. This is a reflection of target strength. Although observations indicate that crater size, target strength, and surface gravity are variables in the formation of complex craters, they do not permit an unequivocal choice between collapse and rebound processes for the formation of complex structures. It may be that both processes act together in the modification of crater morphology during the later stages of excavation. The major emphasis of recent shock metamorphic studies has been toward the development of models of cratering processes. An important contribution has been the identification, through meteoritic contamination in the melt rocks, of the type of bolide at a number of probable impact structures. This has served to strengthen the link between the occurrence of shock metamorphic effects and their origin by hypervelocity meteorite impact.     

Crater size-shape profiles for the moon and mercury: Terrain effects and interplanetary comparisons

New crater size-shape data were compiled for 221 fresh lunar craters and 152 youthful mercurian craters. Terraces and central peaks develop initially in fresh craters on the Moon in the 0–10 km diameter interval. Above a diameter of 65 km all craters are terraced and have central peaks. Swirl floor texture is most common in craters in the size range 20–30 km, but it occurs less frequently as terraces become a dominant feature of crater interiors. For the Moon there is a correlation between crater shape and geomorphic terrain type. For example, craters on the maria are more complex in terms of central peak and terrace detail at any given crater diameter than are craters in the highlands. These crater data suggest that there are significant differences in substrate and/or target properties between maria and highlands. Size-shape profiles for Mercury show that central peak and terrace onset is in the 10–20 km diameter interval; all craters are terraced at 65 km, and all have central peaks at 45 km. The crater data for Mercury show no clear cut terrain correlation. Comparison of lunar and mercurian data indicates that both central peaks and terraces are more abundant in craters in the diameter range 5–75 km on Mercury. Differences in crater shape between Mercury and the Moon may be due to differences in planetary gravitational acceleration (gMercury=2.3gMoon). Also differences between Mercury and the Moon in target and substrate and in modal impact velocity may contribute to affect crater shape.     

Geological overview and cratering model for the Haughton impact structure,Devon Island,Canadian High Arctic

Abstract— The Haughton impact structure has been the focus of systematic, multi‐disciplinary field and laboratory research activities over the past several years. Regional geological mapping has refined the sedimentary target stratigraphy and constrained the thickness of the sedimentary sequence at the time of impact to ?1880 m. New ⁴⁰Ar–³⁹Ar dates place the impact event at ?39 Ma, in the late Eocene. Haughton has an apparent crater diameter of ?23 km, with an estimated rim (final crater) diameter of ?16 km. The structure lacks a central topographic peak or peak ring, which is unusual for craters of this size. Geological mapping and sampling reveals that a series of different impactites are present at Haughton. The volumetrically dominant crater‐fill impact melt breccias contain a calcite‐anhydrite‐silicate glass groundmass, all of which have been shown to represent impact‐generated melt phases. These impactites are, therefore, stratigraphically and genetically equivalent to coherent impact melt rocks present in craters developed in crystalline targets. The crater‐fill impactites provided a heat source that drove a post‐impact hydrothermal system. During this time, Haughton would have represented a transient, warm, wet microbial oasis. A subsequent episode of erosion, during which time substantial amounts of impactites were removed, was followed by the deposition of intra‐crater lacustrine sediments of the Haughton Formation during the Miocene. Present‐day intra‐crater lakes and ponds preserve a detailed paleoenvironmental record dating back to the last glaciation in the High Arctic. Modern modification of the landscape is dominated by seasonal regional glacial and niveal melting, and local periglacial processes. The impact processing of target materials improved the opportunities for colonization and has provided several present‐day habitats suitable for microbial life that otherwise do not exist in the surrounding terrain.     

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 MEMIN research unit: Scaling impact cratering experiments in porous sandstones

Abstract– The MEMIN research unit (Multidisciplinary Experimental and Modeling Impact research Network) is focused on analyzing experimental impact craters and experimental cratering processes in geological materials. MEMIN is interested in understanding how porosity and pore space saturation influence the cratering process. Here, we present results of a series of impact experiments into porous wet and dry sandstone targets. Steel, iron meteorite, and aluminum projectiles ranging in size from 2.5 to 12 mm were accelerated to velocities of 2.5–7.8 km s^?1, yielding craters with diameters between 3.9 and 40 cm. Results show that the target’s porosity reduces crater volumes and cratering efficiency relative to nonporous rocks. Saturation of pore space with water to 50% and 90% increasingly counteracts the effects of porosity, leading to larger but flatter craters. Spallation becomes more dominant in larger‐scale experiments and leads to an increase in cratering efficiency with increasing projectile size for constant impact velocities. The volume of spalled material is estimated using parabolic fits to the crater morphology, yielding approximations of the transient crater volume. For impacts at the same velocity these transient craters show a constant cratering efficiency that is not affected by projectile size.     

Comparison of central pit craters on Mars,Mercury, Ganymede,and the Saturnian satellites

We report on the first results of a large‐scale comparison study of central pit craters throughout the solar system, focused on Mars, Mercury, Ganymede, Rhea, Dione, and Tethys. We have identified 10 more central pit craters on Rhea, Dione, and Tethys than have previously been reported. We see a general trend that the median ratio of the pit to crater diameter (D_p/D_c) decreases with increasing gravity and decreasing volatile content of the crust. Floor pits are more common on volatile‐rich bodies while summit pits become more common as crustal volatile content decreases. Uplifted bedrock from below the crater floor occurs in the central peak upon which summit pits are found and in rims around floor pits, which may or may not break the surface. Peaks on which summit pits are found on Mars and Mercury share similar characteristics to those of nonpitted central peaks, indicating that some normal central peaks undergo an additional process to create summit pits. Martian floor pits do not appear to be the result of a central peak collapse as the median ratio of the peak to crater diameter (D_pk/D_c) is about twice as high for central peaks/summit pits than D_p/D_c values for floor pits. Median D_pk/D_c is twice as high for Mars as for Mercury, reflecting differing crustal strength between the two bodies. Results indicate that a complicated interplay of crustal volatiles, target strength, surface gravity, and impactor energy along with both uplift and collapse are involved in central pit formation. Multiple formation models may be required to explain the range of central pits seen throughout the solar system.     

Low crater frequencies and low model ages in lunar maria: Recent endogenic activity or degradation effects?

Recently a number of studies have identified small lunar geologic structures to be <100 Ma in age using standard remote sensing techniques. Here we present new crater size frequency distributions (CSFDs) and model ages using craters D > 10 m for five small target units: one irregular mare patch (IMP) in Mare Nubium and four regions located on lunar wrinkle ridges in Mare Humorum. For comparison we also date another IMP found in a recent study in Mare Tranquillitatis (Braden et al. 2014 ). Absolute model age (AMA) derivation corresponds to 46 ± 5 Ma and 22 ± 1 Ma for Nubium and Sosigenes IMP, respectively. We show that for IMPs and in nearby control mare regions, similar production-like cumulative log–log SFD slopes of −3 are observed. In contrast, control mare regions in Mare Humorum exhibit shallower equilibrium slopes from −1.83 to −2. Three out of four wrinkle ridges appear to be in equilibrium but with crater lifetimes lower than on the corresponding maria. Low crater frequencies on one wrinkle ridge result in an age of 8.6 ± 1 Ma. This study region contains 80% fresh craters, which suggests that the crater population is still in production indicative of a recent resurfacing event.     

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.     

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.     

The large crater on the small Asteroid (2867) Steins

The maximum size of impact craters on finite bodies marks the largest impact that can occur short of impact induced disruption of the body. Recently attention has started to focus on large craters on small bodies such as asteroids and rocky and icy satellites. Here the large crater on the recently imaged Asteroid (2867) Steins (with crater diameter to mean asteroid radius ratio of 0.79) is shown to follow a limit set by other similar sized bodies with moderate macroporosity (i.e. fractured asteroids). Thus whilst large, the crater size is not novel, nor does it require Steins to possess an extremely large porosity. In one of the components of the binary Asteroid (90) Antiope there is the recently reported presence of an extremely large depression, possibly a crater, with depression diameter to mean asteroid radius ratio of ∼(1.4–1.62). This is consistent with the maximum size of a crater expected from previous observations of very porous rocky bodies (i.e. rubble-pile asteroids). Finally, a relationship between crater diameter (normalised to body radius) is proposed as a function of body porosity which suggests that the doubling of porosity between fractured asteroids and rubble-pile asteroids, nearly doubles the size (D/R value) of the largest crater sustainable on a rocky body.     

Geophysical signature of Målingen,the minor crater of the Lockne–Målingen doublet impact structure

Målingen is the 0.7 km wide minor crater associated to the 10 times larger Lockne crater in the unique Lockne–Målingen doublet. The craters formed at 458 Ma by the impact of a binary asteroid related to the well-known 470 Ma Main Belt breakup event responsible for a large number of Ordovician craters and fossil meteorites. The binary asteroid struck a target sequence including ~500 m of sea water, ~80 m of limestone, ~30 m of dark mud, and a peneplainized Precambrian crystalline basement. Although the Lockne crater has been extensively studied by core drillings and geophysics, little is known about the subsurface morphology of Målingen. We performed magnetic susceptibility and remanence, as well as density, measurements combined with gravity, and magnetic field surveys over the crater and its close vicinity as a base for forward magnetic and gravity modeling. The interior of the crater shows a general magnetic low of 90–100 nT broken by a clustered set of high-amplitude, short wavelength anomalies caused by bodies of mafic rock in the target below the crater and as allogenic blocks in the crater infill. The gravity shows a general −1.4 mgal anomaly over the crater caused by low-density breccia infill and fractured crystalline rocks below the crater floor. The modeling also revealed a slightly asymmetrical shape of the crater that together with the irregular ejecta distribution supports an oblique impact from the east, which is consistent with the direction of impact suggested for the Lockne crater.     

A general cratering-history model and its implications for the lunar highlands

Through analysis of a large number of Monte Carlo and Markov Chain simulations, a model for determining crater accumulation and crater obliteration histories has been derived. The model generally applies to populations of large craters. It predicts that the following relationships hold for subequilibrium-density crater populations: (1) the more negative the production function's exponent, α, (N～D^α) the lower the crater density at which the population size-frequency distribution will significantly depart from its production function; (2) the more negative the production function's exponent, the less obliteration a crater population will sustain after a set number of impacts. Application of the model to the lunar highlands implies (1) the production function for the large craters is highly structured, resembling the observed size-frequency distribution and not the function N～D^?2; (2) even the densely cratered highlands have not attained crater saturation or equilibrium. Direct simulations of the highlands' crater population supports the model's implications.