The rayed crater Zunil and interpretations of small impact craters on Mars

A 10-km diameter crater named Zunil in the Cerberus Plains of Mars created ∼¹⁰⁷ secondary craters 10 to 200 m in diameter. Many of these secondary craters are concentrated in radial streaks that extend up to 1600 km from the primary crater, identical to lunar rays. Most of the larger Zunil secondaries are distinctive in both visible and thermal infrared imaging. MOC images of the secondary craters show sharp rims and bright ejecta and rays, but the craters are shallow and often noncircular, as expected for relatively low-velocity impacts. About 80% of the impact craters superimposed over the youngest surfaces in the Cerberus Plains, such as Athabasca Valles, have the distinctive characteristics of Zunil secondaries. We have not identified any other large (?10 km diameter) impact crater on Mars with such distinctive rays of young secondary craters, so the age of the crater may be less than a few Ma. Zunil formed in the apparently youngest (least cratered) large-scale lava plains on Mars, and may be an excellent example of how spallation of a competent surface layer can produce high-velocity ejecta (Melosh, 1984, Impact ejection, spallation, and the origin of meteorites, Icarus 59, 234-260). It could be the source crater for some of the basaltic shergottites, consistent with their crystallization and ejection ages, composition, and the fact that Zunil produced abundant high-velocity ejecta fragments. A 3D hydrodynamic simulation of the impact event produced 10¹⁰ rock fragments ?10 cm diameter, leading to up to 10⁹ secondary craters ?10 m diameter. Nearly all of the simulated secondary craters larger than 50 m are within 800 km of the impact site but the more abundant smaller (10-50 m) craters extend out to 3500 km. If Zunil is representative of large impact events on Mars, then secondaries should be more abundant than primaries at diameters a factor of ∼1000 smaller than that of the largest primary crater that contributed secondaries. As a result, most small craters on Mars could be secondaries. Depth/diameter ratios of 1300 small craters (10-500 m diameter) in Isidis Planitia and Gusev crater have a mean value of 0.08; the freshest of these craters give a ratio of 0.11, identical to that of fresh secondary craters on the Moon (Pike and Wilhelms, 1978, Secondary-impact craters on the Moon: topographic form and geologic process, Lunar Planet. Sci. IX, 907-909) and significantly less than the value of ∼0.2 or more expected for fresh primary craters of this size range. Several observations suggest that the production functions of Hartmann and Neukum (2001, Cratering chronology and the evolution of Mars, Space Sci. Rev. 96, 165-194) predict too many primary craters smaller than a few hundred meters in diameter. Fewer small, high-velocity impacts may explain why there appears to be little impact regolith over Amazonian terrains. Martian terrains dated by small craters could be older than reported in recent publications.     

Morphologic classification of Martian craters and some implications

A computer data bank containing information on crater sizes, locations, and morphologies for all craters visible on Mariner 9 wide-angle mapping photography was used to construct a crater morphologic classification. Four general classes were constructed that can be interpreted to represent increasing degrees of crater degradation. Fresh class craters are nearly unmodified and consist of deep bowl-shaped craters and deep, flat-floored craters with terraced walls. The slightly modified class consists of deep flat-floored craters that usually have raised rims, but lack the terracing, central peaks, and hummocky floors indicative of unmodified impact crater morphology. Craters in the modified class are rimless and shallow and those in the ghost class are rimless and extremely shallow. Retention ages for fresh (i.e. unmodified) class craters on equatorial cratered terrain range from millions to billions of years, depending on the impact flux history used. If the trend is toward billions of years, then present degradation rates on Mars are low relative to earlier history and most craters in the degraded classes were probably modified in an early (>3.3 b.y.?) period.     

Detecting crater ejecta‐blanket boundaries and constraining source crater regions for boulder tracks and elongated secondary craters on Eros

Abstract– We present results of a numerical model of the dynamics of ejecta emplacement on asteroid 433 Eros. Ejecta blocks represent the coarsest fraction of Eros’ regolith and are important, readily visible, “tracer particles” for crater ejecta‐blanket units that may be linked back to specific source craters. Model results show that the combination of irregular shape and rapid rotation of an asteroid can result in markedly asymmetric ejecta blankets (and, it follows, ejecta block spatial distribution), with locally very sharp/distinct boundaries. We mapped boulder number densities in NEAR‐Shoemaker MSI images across a portion of a predicted sharp ejecta‐blanket boundary associated with the crater Valentine and confirm a distinct and real ejecta‐blanket boundary, significant at least at the 3‐sigma level. Using our dynamical model, we “back track” the landing trajectories of three ejecta blocks with associated landing tracks in an effort to constrain potential source regions where those blocks were ejected from Eros’ surface in impact events. The observed skip distances of the blocks upon landing on Eros’ surface and the landing speeds and elevation angles derived from our model allow us to estimate the coefficient of restitution, ε, of Eros’ surface for impacts of 10‐m‐scale blocks at approximately 5 m s^?1 impact speeds. We find mean values of ε of approximately 0.09–0.18.     

GT-57633 catalogue of Martian impact craters developed for evaluation of crater detection algorithms

Crater detection algorithms (CDAs) are an important subject of the recent scientific research. A ground truth (GT) catalogue, which contains the locations and sizes of known craters, is important for the evaluation of CDAs in a wide range of CDA applications. Unfortunately, previous catalogues of craters by other authors cannot be easily used as GT. In this paper, we propose a method for integration of several existing catalogues to obtain a new craters catalogue. The methods developed and used during this work on the GT catalogue are: (1) initial screening of used catalogues; (2) evaluation of self-consistency of used catalogues; (3) initial registration from three different catalogues; (4) cross-evaluation of used catalogues; (5) additional registrations and registrations from additional catalogues; and (6) fine-tuning and registration with additional data-sets. During this process, all craters from all major currently available manually assembled catalogues were processed, including catalogues by Barlow, Rodionova, Boyce, Kuzmin, and our previous work. Each crater from the GT catalogue contains references to crater(s) that are used for its registration. This provides direct access to all properties assigned to craters from the used catalogues, which can be of interest even to those scientists that are not directly interested in CDAs. Having all these craters in a single catalogue also provides a good starting point for searching for craters still not catalogued manually, which is also expected to be one of the challenges of CDAs. The resulting new GT catalogue contains 57,633 craters, significantly more than any previous catalogue. From this point of view, GT-57633 catalogue is currently the most complete catalogue of large Martian impact craters. Additionally, each crater from the resulting GT-57633 catalogue is aligned with MOLA topography and, during the final review phase, additionally registered/aligned with 1/256° THEMIS-DIR, 1/256° MDIM and 1/256° MOC data-sets. Accordingly, the resulting GT-57633 catalogue can successfully be used as a part of the framework for evaluation of CDAs.     

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.     

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.     

On the depths and shapes of the freshest kilometer-scale simple craters on the lunar maria: A new crater shape model

Recent work on the shapes of small, simple impact craters on the Moon has shown that the parabolic ideal does not well represent the vast majority of these craters. They are hyperbolic in shape and usually resemble a cone more than a parabola. A parabolic shape also does not fit the most commonly held archetype for simple craters in general (Linné), which is also hyperbolic. In addition, Linné itself may not be the best model for fresh simple craters, in terms of cross-sectional shape, although shape data to compare it to have heretofore been lacking. Here, the “free shadowfront method” for determining the shapes of simple craters is used to measure 64 fresh simple craters on five lunar maria to test both assumptions. Laser altimetry cross sections, available for many of the craters measured herein, are used to complement and spot-check the shadow measurement results, and thereby demonstrate the efficacy of the free shadowfront method. A new shape model is established, and two craters that better fit this model than Linné are identified. These are located at 24.45° N/328.12° E and 31.35° N/296.46° E and have diameters of 1.40 and 2.73 km, respectively. An apparent dichotomy between fresh simple craters smaller than 2.5 km and those larger than this is observed. Flat floors are found to be ubiquitous among the larger craters, but rare and small in extent in smaller ones. A slide in one crater which appears to be an incipient flat floor suggests a major mode of formation for these flat floors.     

Shock metamorphism and impact melting in small impact craters on Earth: Evidence from Kamil crater,Egypt

Kamil is a 45 m diameter impact crater identified in 2008 in southern Egypt. It was generated by the hypervelocity impact of the Gebel Kamil iron meteorite on a sedimentary target, namely layered sandstones with subhorizontal bedding. We have carried out a petrographic study of samples from the crater wall and ejecta deposits collected during our first geophysical campaign (February 2010) in order to investigate shock effects recorded in these rocks. Ejecta samples reveal a wide range of shock features common in quartz‐rich target rocks. They have been divided into two categories, as a function of their abundance at thin section scale: (1) pervasive shock features (the most abundant), including fracturing, planar deformation features, and impact melt lapilli and bombs, and (2) localized shock features (the least abundant) including high‐pressure phases and localized impact melting in the form of intergranular melt, melt veins, and melt films in shatter cones. In particular, Kamil crater is the smallest impact crater where shatter cones, coesite, stishovite, diamond, and melt veins have been reported. Based on experimental calibrations reported in the literature, pervasive shock features suggest that the maximum shock pressure was between 30 and 60 GPa. Using the planar impact approximation, we calculate a vertical component of the impact velocity of at least 3.5 km s^?1. The wide range of shock features and their freshness make Kamil a natural laboratory for studying impact cratering and shock deformation processes in small impact structures.     

Changes in abundance and nature of microimpact craters on the surfaces of Australasian microtektites with distance from the proposed source crater location

Abstract– Over 4600 Australasian microtektites from 11 sediment cores along an N–S transect in the Central Indian Ocean have been investigated optically for microimpact features on their surfaces. Detailed scanning electron microscope examination of 68 microtektites along this transect shows 4091 such features. These samples are located between approximate distances of 3900–5000 km from the suggested impact site in Indochina and therefore constitute distal ejecta. The morphology of the microimpacts seems to show distinct variations with distance from the source crater. The total number of microcraters on each microtektite decreases drastically from North to South indicating systematic decrease in the spatial density of the ejecta, and decrease in collisional activity between microtektites with distance from the proposed source crater location. Closer to the proposed source crater location, the microcraters are predominantly small (few μm), pit bearing with radial and concentric cracks, suggestive of violent interparticle collisions. The scenario is reverse farther from the source crater with smaller numbers of impacted microtektites due to increased dispersion of the ejecta and the microcraters are large and shallow, implying gentle collisions with larger particles. These observations provide systematic ground truth for the processes that take place as the ejecta of a large oblique impact which generated the Australasian tektite strewn field is emplaced. The microimpacts appear to take place during the descent of the ejecta and their intensity and number density decrease as a function of the spatial density of the ejecta at any given place and with distance from the source region. These features could help understand processes that take place during ejecta emplacement on planets with substantial atmosphere such as Mars and Venus.     

Multivariate classification methods in planetary sciences

In this paper the problem of the classification of natural samples is discussed. An updated version of the G-mode multivariate statistical method for the classification of natural samples, applicable to a wide range of research fields, is discussed in this paper. This method allows an automatic classification in terms of homogeneous taxonomic units, without any a priori knowledge of the taxonomic structure of the natural observations; it provides informations on the different levels of classification present in the data set under study (classes and subclasses), on the level of information residing in each variable, on the level of similarity and/or difference among homogeneous classes.An earlier version of this method has been widely applied in planetary sciences, astrophysics and geological sciences. The authors give here a resumé of the most interesting results obtained in these different research fields, from the geochemistry of lunar samples to asteroids taxonomy to remote sensing of planetary surfaces. This method is extremely reliable and versatile, and it is suggested that its use be widespread whenever the problem of the classification of a large set of natural samples occurs.     

Polygonal impact craters in the Argyre region,Mars: Evidence for influence of target structure on the final crater morphology

Abstract— Impact craters that in plan view are distinctly polygonal rather than circular or elliptical are common on Mars and other planets (Öhman et al. 2005). Their actual formation mechanism, however, is somewhat debatable. We studied the polygonal craters of different degradational stages in the region of the Argyre impact basin, Mars. The results show that in the same areas, heavily degraded, moderately degraded, and fresh polygonal craters display statistically similar strike distributions of the straight rim segments. The fact that the strike distributions are not dependent on lighting conditions was verified by using two data sets (Viking and MOC‐WA) having different illumination geometries but similar resolutions. In addition, there are no significant differences in the amount of polygonality of craters in different degradational stages. These results clearly imply that large‐scale polygonality is not caused by degradation, but originates from the cratering process itself, concurring with the findings regarding lunar craters by Eppler et al. (1983). The straight rims of polygonal craters apparently reflect areal fracture patterns that prevail for a geologically long time.     

Floor-fractured crater models of the Sudbury Structure,Canada: Implications for initial crater size and crater modification

Abstract The pattern of radial and concentric offset dikes at Sudbury strongly resembles fracture patterns in certain volcanically modified craters on the Moon. Since the Sudbury dikes apparently formed shortly after the impact event, this resemblance suggests that early endogenic modification at Sudbury was comparable to deformation in lunar floor-fractured craters. Although regional deformation has obscured many details of the Sudbury Structure, such a comparison of Sudbury with lunar floor-fractured craters provides two alternative models for the original size and surface structures of the Sudbury basin. First, the Sudbury date pattern can be correlated with fractures in the central peak crater Haldane (36 km in diameter). This comparison indicates an initial Sudbury diameter of between 100 and 140 km but requires loss of a central peak complex for which there is little evidence. Alternatively, comparison of the Sudbury dikes with fractures in the two-ring basin Schrödinger indicates an initial Sudbury diameter of at least ～ 180 km, which is in agreement with other recent estimates for the size of the Sudbury Structure. In addition to constraining the size and structure of the original Sudbury crater, these comparisons also suggest that crater modification may reflect different deformation mechanisms at different sizes. Most lunar floor-fractured craters are attributed to deformation over a shallow, crater-centered intrusion; however, there is no evidence for such an intrusion at Sudbury. Instead, melts from the evolving impact melt sheet probably entered fractures formed by isostatically-induced flexure of the crater floor. Since most of the lunar floor-fractured craters are too small (<100-km diameter) to induce significant isostatic adjustment, crater modification by isostatic uplift apparently is limited to only the largest of craters, whereas deformation over igneous intrusions dominates the modification of smaller craters.     

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.     

Malvinas (Falkland) Plateau structure versus Mjølnir crater: Geophysical workflow template for proposed marine impact craters

A diagnostic geophysical‐based template, supported by modelling, is suggested to be used prior to, or in combination with geological/drilling data, when proposing a marine impact crater. The latter refers to impacts occurring in a marine setting and resulting in structures that are currently partially or totally underwater. The methodology is based on the well‐documented Mjølnir crater in the Barents Sea. The template has been developed in conjunction with the recently proposed and debated impact crater on the Malvinas (Falkland) Plateau in the South Atlantic. Despite their different sizes, their comparison adds to the ambiguous nature of the Malvinas structure and shows that the integrated analysis of seismic and potential field data and modelling is crucial for any interpretation of a marine impact crater without relevant geological information. The proposed workflow template utilizes all available geophysical data and is composed of a series of iterative steps, including a range of alternative nonimpact interpretations that must be discussed and accounted for. Subsequently, further iterative geophysical modelling is required to support and decipher the impact related processes. A more complex impact crater model and additional impact crater features can be resolved by physical property modelling. In all cases, a close spatial correspondence of the defined impact structure with potential field anomalies is a necessity to establish a causal relationship. We suggest that the diagnostic workflow template provides a methodology to be applied to future studies of the Malvinas structure, as well as to proposed marine (and, with minor adaptions, to nonmarine) impact craters in general.     

Review of the Barringer crater studies and views on the crater’s origin

The first scientific studies of Barringer crater, Arizona, USA (also known as the Coon Butte crater), began more than a century ago; however, views on the crater’s origin have been contradictory. At the beginning of the 20th century, D.M. Barringer, a mining engineer, became interested in the possibility of finding large useable iron masses in this crater and searched for these masses for more than 25 years, standing up for the idea of the crater’s meteoric origin, contrary to the objections of opponents who tried to indicate that the crater was caused by terrestrial geological processes. Mining, accompanied by different scientific works, made it possible to obtain reliable data on the structure and impact origin of the crater; however, attempts to find meteorite iron deposits in this crater were unsuccessful. Barringer crater was the first object on the Earth where purposeful studies were performed for many decades and made it possible to develop many criteria of the impact origin of circular geological structures and mechanisms of formation of these structures, as well as to compare this crater with similar morphostructures on the surfaces of other planets. These studies have played an important role in the formation and development of the theory of impact cratering, which has been generally acknowledged in present-day science.     

Ejecta size-velocity relation derived from the distribution of the secondary craters of kilometer-sized craters on Mars

The relation between the size and velocity of impact crater ejecta has been studied by both laboratory experiments and numerical modeling. An alternative method, used here, is to analyze the record of past impact events, such as the distribution of secondary craters on planetary surfaces, as described by Vickery (Icarus 67 (1986) 224; Geophys. Res. Lett. 14 (1987) 726). We first applied the method to lunar images taken by the CLEMENTINE mission, which revealed that the size-velocity relations of ejecta from craters 32 and 40 km in diameter were similar to those derived by Vickery for a crater 39 km in diameter. Next, we studied the distribution of small craters in the vicinity of kilometer-sized craters on three images from the Mars Orbiter Camera (MOC) on board the Mars Global Surveyor (MGS). If these small craters are assumed to be secondaries ejected from the kilometer-sized crater in each image, the ejection velocities are of hundreds of meters per second. These data fill a gap between the previous results of Vickery and those of laboratory studies.