Effect of volatiles and target lithology on the generation and emplacement of impact crater fill and ejecta deposits on Mars

Abstract— Impact cratering is an important geological process on Mars and the nature of Martian impact craters may provide important information as to the volatile content of the Martian crust. Terrestrial impact structures currently provide the only ground‐truth data as to the role of volatiles and an atmosphere on the impact‐cratering process. Recent advancements, based on studies of several well‐preserved terrestrial craters, have been made regarding the role and effect of volatiles on the impact‐cratering process. Combined field and laboratory studies reveal that impact melting is much more common in volatile‐rich targets than previously thought, so impact‐melt rocks, melt‐bearing breccias, and glasses should be common on Mars. Consideration of the terrestrial impact‐cratering record suggests that it is the presence or absence of subsurface volatiles and not the presence of an atmosphere that largely controls ejecta emplacement on Mars. Furthermore, recent studies at the Haughton and Ries impact structures reveal that there are two discrete episodes of ejecta deposition during the formation of complex impact craters that provide a mechanism for generating multiple layers of ejecta. It is apparent that the relative abundance of volatiles in the near‐surface region outside a transient cavity and in the target rocks within the transient cavity play a key role in controlling the amount of fluidization of Martian ejecta deposits. This study shows the value of using terrestrial analogues, in addition to observational data from robotic orbiters and landers, laboratory experiments, and numerical modeling to explore the Martian impact‐cratering record.     

Predictions for the LCROSS mission

Abstract— We describe the results of a variety of model calculations for predictions of observable results of the LCROSS mission to be launched in 2009. Several models covering different aspects of the event are described along with their results. Our aim is to bracket the range of expected results and produce a useful guide for mission planning. In this paper, we focus on several different questions, which are modeled by different methods. The questions include the size of impact crater, the mass, velocity, and visibility of impact ejecta, and the mass and temperature of the initial vapor plume. The mass and velocity profiles of the ejecta are of primary interest, as the ejecta will be the main target of the S‐S/C observations. In particular, we focus on such quantities as the amount of mass that reaches various heights. A height of 2 km above the target is of special interest, as we expect that the EDUS impact will take place on the floor of a moderate‐sized crater ?30 km in diameter, with a rim height of 1–2 km. The impact ejecta must rise above the crater rim at the target site in order to scatter sunlight and become visible to the detectors aboard the S‐S/C. We start with a brief discussion of crater scaling relationships as applied to the impact of the EDUS second stage and resulting estimated crater diameter and ejecta mass. Next we describe results from the RAGE hydrocode as applied to modeling the short time scale (t 0.1 s) thermal plume that is expected to occur immediately after the impact. We present results from several large‐scale smooth‐particle hydrodynamics (SPH) calculations, along with results from a ZEUS‐MP hydrocode model of the crater formation and ejecta mass‐velocity distribution. We finish with two semi‐analytic models, the first being a Monte Carlo model of the distribution of expected ejecta, based on scaling models using a plausible range of crater and ejecta parameters, and the second being a simple model of observational predictions for the shepherding spacecraft (S‐S/C) that will follow the impact for several minutes until its own impact into the lunar surface. For the initial thermal plume, we predict an initial expansion velocity of ?7 km s^?1, and a maximum temperature of ?1200 K. Scaling relations for crater formation and the SPH calculation predict a crater with a diameter of ?15 m, a total ejecta mass of ?10⁶kg, with ?10⁴kg reaching an altitude of 2 km above the target. Both the SPH and ZEUS‐MP calculations predict a maximum ejecta velocity of ?1 km s^?1. The semi‐analytic Monte Carlo calculations produce more conservative estimates (by a factor of ?5) for ejecta at 2 km, but with a large dispersion in possible results.     

Impact crater formation in icy layered terrains on Mars

Abstract— We present numerical simulations of crater formation under Martian conditions with a single near‐surface icy layer to investigate changes in crater morphology between glacial and interglacial periods. The ice fraction, thickness, and depth to the icy layer are varied to understand the systematic effects on observable crater features. To accurately model impact cratering into ice, a new equation of state table and strength model parameters for H₂O are fitted to laboratory data. The presence of an icy layer significantly modifies the cratering mechanics. Observable features demonstrated by the modeling include variations in crater morphometry (depth and rim height) and icy infill of the crater floor during the late stages of crater formation. In addition, an icy layer modifies the velocities, angles, and volumes of ejecta, leading to deviations of ejecta blanket thickness from the predicted power law. The dramatic changes in crater excavation are a result of both the shock impedance and the strength mismatch between layers of icy and rocky materials. Our simulations suggest that many of the unusual features of Martian craters may be explained by the presence of icy layers, including shallow craters with well‐preserved ejecta blankets, icy flow related features, some layered ejecta structures, and crater lakes. Therefore, the cratering record implies that near‐surface icy layers are widespread on Mars.     

Bosumtwi impact structure,Ghana: Evidence for fluidized emplacement of the ejecta

The about 10.5 km diameter Bosumtwi impact crater is one of the youngest large impact structures on Earth. The crater rim is readily noticed on topographic maps or in satellite imagery. It defines a circular basin filled by water (Lake Bosumtwi) and lacustrine sediments. The morphology of this impact structure is also characterized by a circular plateau extending beyond the rim and up to 9–10 km from the center of the crater (about 2 crater radii). This feature comprises a shallow ring depression, also described as an annular moat, and a subdued circular ridge at its outer edge. The origin of this outermost feature could so far not be elucidated based on remote sensing data only. Our approach combines detailed topographic analysis, including roughness mapping, with airborne radiometric surveys (mapping near‐surface K, Th, U concentrations) and field observations. This provides evidence that the moat and outer ring are features inherited from the impact event and represent the partially eroded ejecta layer of the Bosumtwi impact structure. The characteristics of the outer ridge indicate that ejecta emplacement was not purely ballistic but requires ejecta fluidization and surface flow. The setting of Bosumtwi ejecta can therefore be considered as a terrestrial analog for rampart craters, which are common on Mars and Venus, and also found on icy bodies of the outer solar system (e.g., Ganymede, Europa, Dione, Tethys, and Charon). Future studies at Bosumtwi may therefore help to elucidate the mechanism of formation of rampart craters.     

Terrestrial impact: The record in the rocks*

Abstract— Approximately 130 terrestrial craters are currently known. They range up to 140 km, and perhaps as much as 200 km, in diameter and from Recent to ～2 billion years in age. The known sample, however, is highly biased to geologically young craters on the better known cratonic areas. The sample is also deficient in small (D < 20 km) craters compared to other planetary bodies. These biases are largely the result of active terrestrial geologic processes and their effects have to be considered when interpreting the record. The strength of the terrestrial cratering record lies in the availability of ground truth data, particularly on the structural and lithological nature of craters, which can be interpreted to understand and constrain large-scale impact processes. Some contributions include the definition of the concept of transient cavity formation and structural uplift during cratering events. Depths of excavation are poorly constrained, as very few terrestrial craters have preserved ejecta. Unlike their planetary counterparts, terrestrial impact craters are mostly recognized not by morphology but by the occurrence of characteristic shock metamorphic effects. Their study has led to models of shock wave attenuation and an understanding of the character and formation of various impact-lithologies, including impact melt rocks. They, in turn, aid in interpreting the nature of extraterrestrial samples, particularly samples from the lunar highlands. The recognition of diagnostic shock metamorphic effects and the signature of projectile contamination through geochemical anomalies in impact lithologies provide the basis for recognizing the impact signature in K/T boundary samples. The record also provides a basis for testing hypotheses of periodic cometary showers. Although inherently not suitable to define short wavelength periods in time due to relatively large uncertainties associated with crater ages, the current record shows no evidence of periodicity. Future directions in terrestrial impact studies will likely continue to focus on the K/T and related problems, including the recognition of other impact signatures in the stratigraphic record. Some emphasis will likely be given to the economic potential of craters and individual large structures, such as Sudbury, will provide an increasingly better understood context for interpreting planetary impact craters. To live up to the full potential of the record to constrain impact processes, however, more basic characterization studies are required, in addition to emphasis on topical areas of study.     

Hydrocode modeling of the Sierra Madera impact structure

Abstract— We present the first hydrocode simulations of the formation of the Sierra Madera structure (west Texas, USA), which was caused by an impact into a thick sedimentary target sequence. We modeled Sierra Madera using the iSALE hydrocode, and here we present two best‐fit models: 1) a crater with a rim (final crater) diameter of ?12 km, in agreement with previous authors' interpretations of the original structure, and 2) a crater ?16 km in diameter with increased postimpact erosion. Both models fit some of the geologic observational data, but discrepancies with estimates of peak shock pressure, extent of deformation, and stratigraphic displacement remain. This study suggests that Sierra Madera may be a larger crater than previously reported and illustrates some of the challenges in simulating impact deformation of sedimentary lithologies. As many terrestrial craters possess some amount of sedimentary rocks in the target sequence, numerical models of impacts into sedimentary targets are essential to our understanding of target rock deformation and the mechanics of crater formation.     

High‐resolution studies of double‐layered ejecta craters: Morphology,inherent structure,and a phenomenological formation model

The ejecta blankets of impact craters in volatile‐rich environments often possess characteristic layered ejecta morphologies. The so‐called double‐layered ejecta (DLE) craters are characterized by two ejecta layers with distinct morphologies. The analysis of high‐resolution image data, especially HiRISE and CTX, provides new insights into the formation of DLE craters. A new phenomenological excavation and ejecta emplacement model for DLE craters is proposed based on a detailed case study of the Martian crater Steinheim—a well‐preserved DLE crater—and studies of other DLE craters. The observations show that the outer ejecta layer is emplaced as medial and distal ejecta that propagate outwards in a debris avalanche or (if saturated with water) a debris flow mode after landing, overrunning previously formed secondary craters. In contrast, the inner ejecta layer is formed by a translational slide of the proximal ejecta deposits during the emplacement stage that overrun and superimpose parts of the outer ejecta layer. Based on our model, DLE craters on Mars are the result of an impact event into a rock/ice mixture that produces large amounts of shock‐induced vaporization and melting of ground ice, leading to high ejection angles, proximal landing positions, and an ejecta curtain with relatively wet (in terms of water in liquid form) composition in the distal part versus dryer composition in the proximal part. As a consequence, basal melting of ice components in the ejecta at the transient crater rim, which is induced by frictional heating and the enhanced pressure at depth, initiates an outwards directed collapse of crater rim material in a translational slide mode. Our results indicate that similar processes may also be applicable for other planetary bodies with volatile‐rich environments, such as Ganymede, Europa, and the Earth.     

A new method to determine the direction of impact: Asymmetry of concentric impact craters as observed in the field (Lockne), on Mars,in experiments,and simulations

Most impacts occur at an angle with respect to the horizontal plane. This is primarily reflected in the ejecta distribution, but at very low angle structural asymmetries such as elongation of the crater and nonradial development of the central peak become apparent. Unfortunately, impact craters with pristine ejecta layers are rare on Earth and also in areas with strong past or ongoing surface erosion on other planetary bodies, and the structural analysis of central peaks requires good exposures or even on‐site access to outcrop. However, target properties are known to greatly influence the shape of the crater, especially the relatively common target configuration of a weaker layer covering a more rigid basement. One such effect is the formation of concentric craters, i.e., a nested, deeper, inner crater surrounded by a shallow, outer crater. Here, we show that with decreasing impact angle there is a downrange shift of the outer crater with respect to the nested crater. We use a combination of (1) field observation and published 3‐D numerical simulation of one of the best examples of a terrestrial, concentric impact crater formed in a layered target with preserved ejecta layer: the Lockne crater, Sweden; (2) remote sensing data for three pristine, concentric impact craters on Mars with preserved ejecta layers further constraining the direction of impact; as well as (3) laboratory impact experiments, to develop the offset in crater concentricity into a complementary method to determine the direction of impact for layered‐target craters with poorly preserved ejecta layers.     

Distribution of ejecta from small impact craters

The present study focuses both on the influence of impact scale on ejecta expansion and on specific features of ejecta deposits around relatively small craters (i.e., those a few kilometers in width). The numerical model is based on the SOVA multimaterial multidimensional hydrocode, considering subaerial vertical impacts only, applying a 2‐D version of the code to projectiles of 100, 300, and 1000 m diameter. Ejecta can roughly be divided into two categories: “ballistic” ejecta and “convective” ejecta; the ballistic ejecta are the ejecta with which the air interacts only slightly, while the convective ejecta motion is entirely defined by the air flow. The degree of particle/air interaction can be defined by the time/length of particle travel before deceleration. Ejecta size‐distributions for the impacts modeled can be described by the same power law, but the size of maximum fragment increases with scale. There is no qualitative difference between the 100 m diameter projectile case and the 300 m diameter projectile impact. In both cases, fine ejecta decelerate in the air at a small distance from launching point and then rise to the stratosphere by air flows induced by the impacts. In the 1000 m‐scale impact, the mass of ejecta is so large that it moves the atmosphere itself to high altitudes. Thus, the atmosphere cannot decelerate even the fine ejecta and they consequently expand to the rarefied upper atmosphere. In the upper atmosphere, even fine ejecta move more or less ballistically and therefore may travel to high altitudes.     

Stardust interstellar dust calibration: Hydrocode modeling of impacts on Al‐1100 foil at velocities up to 300 km s−1 and validation with experimental data

Abstract– We present initial results from hydrocode modeling of impacts on Al‐1100 foils, undertaken to aid the interstellar preliminary examination (ISPE) phase for the NASA Stardust mission interstellar dust collector tray. We used Ansys’ AUTODYN to model impacts of micrometer‐scale, and smaller projectiles onto Stardust foil (100 μm thick Al‐1100) at velocities up to 300 km s^?1. It is thought that impacts onto the interstellar dust collector foils may have been made by a combination of interstellar dust particles (ISP), interplanetary dust particles (IDP) on comet, and asteroid derived orbits, β micrometeoroids, nanometer dust in the solar wind, and spacecraft derived secondary ejecta. The characteristic velocity of the potential impactors thus ranges from <<1 to a few km s^?1 (secondary ejecta), approximately 4–25 km s^?1 for ISP and IDP, up to hundreds of km s^?1 for the nanoscale dust reported by Meyer‐Vernet et al. (2009) . There are currently no extensive experimental calibrations for the higher velocity conditions, and the main focus of this work was therefore to use hydrocode models to investigate the morphometry of impact craters, as a means to determine an approximate impactor speed, and thus origin. The model was validated against existing experimental data for impact speeds up to approximately 30 km s^?1 for particles ranging in density from 2.4 kg m^?3 (glass) to 7.8 kg m^?3 (iron). Interpolation equations are given to predict the crater depth and diameter for a solid impactor with any diameter between 100 nm and 4 μm and density between 2.4 and 7.8 kg m^?3.     

Impact craters in the northern hemisphere of Mars: Layered ejecta and central pit characteristics

Abstract— Mars Global Surveyor (MGS) and Mars Odyssey data are being used to revise the Catalog of Large Martian Impact Craters. Analysis of data in the revised catalog provides new details on the distribution and morphologic details of 6795 impact craters in the northern hemisphere of Mars. This report focuses on the ejecta morphologies and central pit characteristics of these craters. The results indicate that single‐layer ejecta (SLE) morphology is most consistent with impact into an ice‐rich target. Double‐layer ejecta (DLE) and multiple‐layer ejecta (MLE) craters also likely form in volatile‐rich materials, but the interaction of the ejecta curtain and target‐produced vapor with the thin Martian atmosphere may be responsible for the large runout distances of these ejecta. Pancake craters appear to be a modified form of double‐layer craters where the thin outer layer has been destroyed or is unobservable at present resolutions. Pedestal craters are proposed to form in an icerich mantle deposited during high obliquity periods from which the ice has subsequently sublimated. Central pits likely form by the release of vapor produced by impact into ice‐soil mixed targets. Therefore, results from the present study are consistent with target volatiles playing a dominant role in the formation of crater morphologies found in the Martian northern hemisphere.     

Numerical Modeling of the Largest Terrestrial Meteorite Craters

Multi-ring impact basins have been found on the surfaces of almost all planetary bodies in the Solar system with solid crusts. The details of their formation mechanism are still unclear. We present results of our numerical modeling of the formation of the largest known terrestrial impact craters. The geological and geophysical data on these structures accumulated over many decades are used to place constraints on the parameters of available numerical models with a dual purpose: (i) to choose parameters in available mechanical models for the crustal response of planetary bodies to a large impact and (ii) to use numerical modeling to refine the possible range of original diameters and the morphology of partially eroded terrestrial craters. We present numerical modeling results for the Vredefort, Sudbury, Chicxulub, and Popigai impact craters and compare these results with available geological and geophysical information.     

Modeling the formation of the K-Pg boundary layer

In this paper we investigate the formation of the Cretaceous-Paleogene (K-Pg) boundary layer through numerical modeling. The K-Pg layer is widely agreed to be composed of meteoritic material and target rock from the Chicxulub impact site, that has been ejected around the globe and mixed with local material during final deposition. The observed composition and thickness of the K-Pg boundary layer changes with azimuth and distance from the impact site. We have run a suite of numerical simulations to investigate whether we can replicate the observational data, with a focus on the distal K-Pg layer and the impact glasses at proximal sites such as Beloc, Haiti. Previous models of the K-Pg ejecta have assumed an initial velocity distribution and tracked the ejecta to its final destination. Here, we attempt to model the entire process, from impact to the arrival of the ejecta around the globe. Our models replicate the observed ejecta thickness at proximal sites, and the modeled ejecta is composed of sediments and silicate basement rocks, in agreement with observational data. Models that use a 45° impact angle are able to replicate the total ejecta and iridium volume at distal sites, and the majority of the ejecta is composed of meteorite and target sediments. Sub-vertical impacts generate too little iridium, and oblique impacts of ?30 degrees generate too much. However, in contrast to observations, models that involve ballistic transport of ejecta lead to ejecta thickness decreasing with increasing distance, and are unable to transport shocked minerals (quartz and zircon) from the Chicxulub basement rocks around the globe. We suggest that much of the K-Pg ejecta is transported non-ballistically, and that the most plausible mechanism is through re-distribution from a hot, expanding atmosphere. The results are important for future investigations of the environmental effects of the Chicxulub impact.     

An atmospheric blast/thermal model for the formation of high‐latitude pedestal craters

Abstract— Although tenuous, the atmosphere of Mars affects the evolution of impact‐generated vapor. Early‐time vapor from a vertical impact expands symmetrically, directly transferring a small percentage of the initial kinetic energy of impact to the atmosphere. This energy, in turn, induces a hemispherical shock wave that propagates outward as an intense airblast (due to high‐speed expansion of vapor) followed by a thermal pulse of extreme atmospheric temperatures (from thermal energy of expansion). This study models the atmospheric response to such early‐time energy coupling using the CTH hydrocode written at Sandia National Laboratories. Results show that the surface surrounding a 10 km diameter crater (6 km “apparent” diameter) on Mars will be subjected to intense winds (?200 m/s) and extreme atmospheric temperatures. These elevated temperatures are sufficient to melt subsurface volatiles at a depth of several centimeters for an ice‐rich substrate. Ensuing surface signatures extend to distal locations (?4 apparent crater diameters for a case of 0.1% energy coupling) and include striations, thermally armored surfaces, and/or ejecta pedestals—all of which are exhibited surrounding the freshest high‐latitude craters on Mars. The combined effects of the atmospheric blast and thermal pulse, resulting in the generation of a crater‐centered erosion‐resistant armored surface, thus provide a new, very plausible formation model for high‐latitude Martian pedestal craters.     

Reconstruction of the Morasko meteoroid impact—Insight from numerical modeling

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

Hydrothermal alteration at the Lonar Lake impact structure,India: Implications for impact cratering on Mars

Abstract— The 50,000 year old, 1.8 km diameter Lonar crater is one of only two known terrestrial craters to be emplaced in basaltic target rock (the 65 million year old Deccan Traps). The composition of the Lonar basalts is similar to martian basaltic meteorites, which establishes Lonar as an excellent analogue for similarly sized craters on the surface of Mars. Samples from cores drilled into the Lonar crater floor show that there are basaltic impact breccias that have been altered by post‐impact hydrothermal processes to produce an assemblage of secondary alteration minerals. Microprobe data and X‐ray diffraction analyses show that the alteration mineral assemblage consists primarily of saponite, with minor celadonite, and carbonate. Thermodynamic modeling and terrestrial volcanic analogues were used to demonstrate that these clay minerals formed at temperatures between 130°C and 200°C. By comparing the Lonar alteration assemblage with alteration at other terrestrial craters, we conclude that the Lonar crater represents a lower size limit for impact‐induced hydrothermal activity. Based on these results, we suggest that similarly sized craters on Mars have the potential to form hydrothermal systems, as long as liquid water was present on or near the martian surface. Furthermore, the Fe‐rich alteration minerals produced by post‐impact hydrothermal processes could contribute to the minor iron enrichment associated with the formation of the martian soil.     

The planforms of low‐angle impact craters in the northern hemisphere of Mars

Abstract— We have surveyed Martian impact craters greater than 5 km in diameter using Viking and thermal emission imaging system (THEMIS) imagery to evaluate how the planform of the rim and ejecta changes with decreasing impact angle. We infer the impact angles at which the changes occur by assuming a sin²θ dependence for the cumulative fraction of craters forming below angle θ. At impact angles less than ?40° from horizontal, the ejecta become offset downrange relative to the crater rim. As the impact angle decreases to less than ?20°, the ejecta begin to concentrate in the cross‐range direction and a “forbidden zone” that is void of ejecta develops in the uprange direction. At angles less than ?10°, a “butterfly” ejecta pattern is generated by the presence of downrange and uprange forbidden zones, and the rim planform becomes elliptical with the major axis oriented along the projectile's direction of travel. The uprange forbidden zone appears as a “V” curving outward from the rim, but the downrange forbidden zone is a straight‐edged wedge. Although fresh Martian craters greater than 5 km in diameter have ramparts indicative of surface ejecta flow, the ejecta planforms and the angles at which they occur are very similar to those for lunar craters and laboratory impacts conducted in a dry vacuum. The planforms are different from those for Venusian craters and experimental impacts in a dense atmosphere. We interpret our results to indicate that Martian ejecta are first emplaced predominantly ballistically and then experience modest surface flow.     

Fates of satellite ejecta in the Saturn system

We use conventional numerical integrations to assess the fates of impact ejecta in the Saturn system. For specificity we consider impact ejecta launched from four giant craters on three satellites: Herschel on Mimas, Odysseus and Penelope on Tethys, and Tirawa on Rhea. Speeds, trajectories, and size of the ejecta are consistent with impact on a competent surface (“spalls”) and into unconsolidated regolith. We do not include near-field effects, jetting, or effects peculiar to highly oblique impact. Ejecta are launched at velocities comparable to or exceeding the satellite's escape speed. Most ejecta are swept up by the source moon on time-scales of a few to several decades, and produce craters no larger than 19 km in diameter, with typical craters in the range of a few km. As much as 17% of ejecta reach satellites other than the source moon. Our models generate cratering patterns consistent with a planetocentric origin of most small impact craters on the saturnian icy moons, but the predicted craters tend to be smaller than putative Population II craters. We conclude that ejecta from the known giant craters in the saturnian system do not fully account for Population II craters.     

Testing landslide and atmospheric-effects models for the formation of double-layered ejecta craters on Mars

Double-layered ejecta (DLE) craters are distinctive among the variety of crater morphologies observed on Mars, but the mechanism by which they form remains under debate. We assess two ejecta emplacement mechanisms: (1) atmospheric effects from ejecta curtain-induced vortices or a base surge and (2) ballistic emplacement followed by a landslide of ejecta assisted by either surface- or pore-ice. We conduct a morphological analysis of the ejecta facies for three DLE craters which impacted into irregular pre-existing topography. We find that the unique topographic environments affected the formation of grooves and the inner facies, and thus appear to be inconsistent with an atmospheric-effects origin but are supportive of the landslide hypothesis. We distinguish between the two landslide models (lubrication by either surface- or pore-ice) by assessing relationships between DLE crater ejecta and morphologic features indicative of buried ice deposits, including sublimation pits, ring-mold craters, expanded secondary craters, and excess ejecta craters. The association of DLE craters with these features suggests that surface ice was present at the time of the impacts that formed the DLE craters. We also compare the Froude numbers of DLE crater ejecta to landslides, and find that the ejecta of DLE craters are kinematically and frictionally similar to terrestrial landslides that overran glaciers. This suggests that the grooves on DLE craters may plausibly form through the same shear/splitting mechanism as the landslides. In summary, our analysis supports the hypothesis that DLE craters form through meteoroid impacts into decameters-thick surface ice deposits (emplaced during periods of higher obliquity) followed by ejecta sliding on the ice.     

Early fracturing and impact residue emplacement: Can modelling help to predict their location in major craters?

Abstract— Understanding the nature and composition of larger extraterrestrial bodies that may collide with the Earth is important. One source of information lies in the record of ancient impact craters, some of which have yielded chemical information as to the impacting body. Many deeply eroded craters have no remaining melt sheet or ejecta yet may contain impactor residue within basement fractures. The emplacement mechanism for fractionated siderophile residues is likely to be gaseous, although, melt droplets and some solid materials may survive. For breccia‐ and melt‐filled fractures to contain extraterrestrial material, they must form very early in the impact process. Most current numerical models do not dwell on the formation and location of early major fractures, although, fractures in and around small craters on brittle glass exposed to hypervelocity impact in low Earth orbit have been successfully simulated. Modelling of fracture development associated with larger craters may help locate impact residues and test the models themselves.