The Carancas meteorite impact crater,Peru: Geologic surveying and modeling of crater formation and atmospheric passage

Abstract— The recent Carancas meteorite impact event caused a worldwide sensation. An H4–5 chondrite struck the Earth south of Lake Titicaca in Peru on September 15, 2007, and formed a crater 14.2 m across. It is the smallest, youngest, and one of two eye‐witnessed impact crater events on Earth. The impact violated the hitherto existing view that stony meteorites below a size of 100 m undergo major disruption and deceleration during their passage through the atmosphere and are not capable of producing craters. Fragmentation occurs if the strength of the meteoroid is less than the aerodynamic stresses that occur in flight. The small fragments that result from a breakup rain down at terminal velocity and are not capable of producing impact craters. The Carancas cratering event, however, demonstrates that meter‐sized stony meteoroids indeed can survive the atmospheric passage under specific circumstances. We present results of a detailed geologic survey of the crater and its ejecta. To constrain the possible range of impact parameters we carried out numerical models of crater formation with the iSALE hydrocode in two and three dimensions. Depending on the strength properties of the target, the impact energies range between approximately 100–1000 MJ (0.024–0.24 t TNT). By modeling the atmospheric traverse we demonstrate that low cosmic velocities (12–14 kms^?1) and shallow entry angles (<20 °) are prerequisites to keep aerodynamic stresses low (<10 MPa) and thus to prevent fragmentation of stony meteoroids with standard strength properties. This scenario results in a strong meteoroid deceleration, a deflection of the trajectory to a steeper impact angle (40–60 °), and an impact velocity of 350–600 ms^?1, which is insufficient to produce a shock wave and significant shock effects in target minerals. Aerodynamic and crater modeling are consistent with field data and our microscopic inspection. However, these data are in conflict with trajectories inferred from the analysis of infrasound signals.     

GRAVITY RECONNAISSANCE AT THREE MAURITANIAN CRATERS OF EXPLOSIVE ORIGIN

We have carried out reconnaissance gravity surveys across three Mauritanian craters: Aouelloul, an undoubted meteorite crater; Tenoumer, a probable meteorite crater with a unique array of concentric dikes on its outer rim flanks containing xenoliths of country rock showing abundant shock artifacts; and Temimichat Ghallaman, a crater of possible meteorite impact origin. All three have residual negative gravity anomalies associated with their interiors. In all cases the gravity values return to “normal” immediately outside their rims. At Tenoumer the anomaly has the form and magnitude expected for a meteorite crater which has been subsequently in-filled with unconsolidated sediments to the level of the surrounding country. Maximum depth from the present crater floor to the bottom of the sedimentary fill (top of the original crater floor) is at least 750 feet. With a rim-rim diameter of 6,300 feet, the origin depth/diameter ratio of about 1:8 is virtually identical with that of Meteor Crater, Arizona. Temimichat, with a rim-rim diameter of 2,100 to 2,400 feet, is somewhat larger than has been previously reported. If it is meteoritic in origin the gravity data dictate a surprisingly shallow structure, with a depth from the present floor to the original crater floor of 150 feet maximum and an original depth/diameter ratio of 1:15. No positive evidence for an impact origin has yet been found for Temimichat. Aouelloul is also larger than generally reported, with a rim-rim diameter averaging 1,275 feet. As for Temimichat the gravity data dictate a remarkably shallow structure having a depth/diameter ratio of about 1:13. The combination of a shallow depth and a reasonably high rim apparently requires a scaled depth of burst for the impact event substantially in excess of 0.50, a value previously considered a maximum for explosive impacts. The morphological resemblance between Temimichat and Aouelloul is striking but, without additional evidence, this fact alone cannot be used to infer a similar origin.     

Evidence for a meteoritic origin of the September 15, 2007, Carancas crater

Abstract— On September 15th, 2007, around 11:45 local time in Peru, near the Bolivian border, the atmospheric entry of a meteoroid produced bright lights in the sky and intense detonations. Soon after, a crater was discovered south of Lake Titicaca. These events have been detected by the Bolivian seismic network and two infrasound arrays operating for the Comprehensive Nuclear‐Test‐Ban Treaty Organization, situated at about 80 and 1620 km from the crater. The localization and origin time computed with the seismic records are consistent with the reported impact. The entry elevation and azimuthal angles of the trajectory are estimated from the observed signal time sequences and back‐azimuths. From the crater diameter and the airwave amplitudes, the kinetic energy, mass and explosive energy are calculated. Using the estimated velocity of the meteoroid and similarity criteria between orbital elements, an association with possible parent asteroids is attempted. The favorable setting of this event provides a unique opportunity to evaluate physical and kinematic parameters of the object that generated the first actual terrestrial meteorite impact seismically recorded.     

Mineralogy,petrology, and distribution of meteorites at the Whitecourt crater,Alberta, Canada

The Whitecourt meteorite impact crater, Alberta, Canada is a rare example of a well‐preserved small impact structure, with which thousands of meteorite fragments are associated. As such, this crater represents a unique opportunity to investigate the effect of a low‐energy impact event on an impacting iron bolide. Excellent documentation of meteorite fragment locations and characteristics has generated a detailed distribution map of both shrapnel and regmaglypted meteorite types. The meteorites' distribution, and internal and external characteristics support a low‐altitude breakup of the impactor which caused atmospherically ablated (regmaglypted) meteorites to fall close to the crater and avoid impact‐related deformation. In contrast, shrapnel fragments sustained deformation at macro‐ and microscales resulting from the catastrophic disruption of the impactor. The impactor was significantly fragmented along pre‐existing planes of weakness, including kamacite lamellae and inclusions, resulting in a bias toward low‐mass (<100 g) fragments. Meteorite mineralogy was investigated and the accessory minerals were found to be dominated by sulfides and phosphides with rare carlsbergite, consistent with other low‐Ni IIIAB iron meteorites. Considerations of the total mass of meteoritic material recovered at the site relative to the probable fraction of the impactor that was preserved based on modeling suggests that the crater was formed by a higher velocity, lower mass impactor than previously inferred.     

Scaling of oblique impacts in frictional targets: Implications for crater size and formation mechanisms

Almost every meteorite impact occurs at an oblique angle of incidence, yet the effect of impact angle on crater size or formation mechanism is only poorly understood. This is, in large part, due to the difficulty of inferring impactor properties, such as size, velocity and trajectory, from observations of natural craters, and the expense and complexity of simulating oblique impacts using numerical models. Laboratory oblique impact experiments and previous numerical models have shown that the portion of the projectile’s kinetic energy that is involved in crater excavation decreases significantly with impact angle. However, a thorough quantification of planetary-scale oblique impact cratering does not exist and the effect of impact angle on crater size is not considered by current scaling laws. To address this gap in understanding, we developed iSALE-3D, a three-dimensional multi-rheology hydrocode, which is efficient enough to perform a large number of well-resolved oblique impact simulations within a reasonable time. Here we present the results of a comprehensive numerical study containing more than 200 three-dimensional hydrocode-simulations covering a broad range of projectile sizes, impact angles and friction coefficients. We show that existing scaling laws in principle describe oblique planetary-scale impact events at angles greater than 30° measured from horizontal. The displaced mass of a crater decreases with impact angle in a sinusoidal manner. However, our results indicate that the assumption that crater size scales with the vertical component of the impact velocity does not hold for materials with a friction coefficient significantly lower than 0.7 (sand). We found that increasing coefficients of friction result in smaller craters and a formation process more controlled by impactor momentum than by energy.     

Numerical modeling of oblique hypervelocity impacts on strong ductile targets

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

Constraints on the size of the Vredefort impact crater from numerical modeling

Abstract— The Vredefort structure in South Africa was created by a meteorite impact about two billion years ago. Since that time, the crater has been deeply eroded; so to estimate its original size, researchers have had to rely heavily upon comparison to other terrestrial impact structures. Recent estimates of the original crater diameter range from 160 km to as much as 400 km. In this study, we combined the capabilities of both hydrocode and finite-element modeling, using the former to predict where the pressure of an impact-generated shock wave would have been high enough to form planar deformation features (PDFs) and shatter cones and the latter to follow the subsequent displacement of these shock isobars during the collapse of the crater. We established constraints on the sizes of the projectile and the transient crater (and, thus, on the size of the final crater) by comparing the observed locations of PDFs around Vredefort to the results of our simulations of impacts by projectiles of various sizes. These simulations indicate that a rocky projectile with a diameter of ～10 km, impacting vertically at a velocity of 20 km/s generates shock pressures that are consistent with the distribution of PDFs around Vredefort. These projectile parameters correspond to a transient crater ～80 km in diameter or a final crater ～120–160 km in diameter. Allowing for uncertainties in our modeling procedures, we consider final craters 120 to 200 km in diameter to be consistent with the observed locations of PDFs at Vredefort. The shock pressure contour corresponding to the formation of shatter cones is almost horizontal near the surface, making the locations of these features less useful constraints on the crater size. However, they may provide a constraint on the amount of erosion that has occurred since the impact.     

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.     

A Re-examination of the Crater near Crestone,Colorado*

Abstract The Crestone Crater is an elliptical bowl measuring 355 feet by 246 feet with a mean depth of 23 feet. It lies in unconsolidated sand on the surface of an alluvial fan at the base of the Sangre de Cristo Mountain Range in the San Luis Valley, Colorado (37° 54′ N, 105° 39′ W). Aerial photographs show the crater as a striking feature in its setting on a gently undulating terrain. We examined the site in August 1963 to appraise the possibility that it was formed by meteorite or comet impact. The crater and its vicinity were mapped at two-foot contour intervals, and two lines of auger-hole samples, eight feet deep, were collected across the crater. Sand from the rim and floor is similar in grain size and composition to that several miles to the north and south. It is barren of meteoritic debris, nickel-iron spherules, rock flour, and impact glass. The crater is less than half as deep relative to its diameter as known meteorite explosion craters. Furthermore, topographic profiles indicate that the crater does not form a depression in the land surface. The crater rim is a positive feature enclosing a basin that has a floor approximately level with the surface of the alluvial fan on which it lies. In the absence of any mineralogic or topographic evidence of impact or explosion, we conclude that this landform is not meteoritic or cometary in origin.     

Discovering research value in the Campo del Cielo,Argentina, meteorite craters

Abstract The Campo del Cielo meteorite crater field in Argentina contains at least 20 small meteorite craters, but a recent review of the field data and a remote sensing study suggest that there may be more. The fall occurred ～4000 years ago into a uniform loessy soil, and the craters are well enough preserved so that some of their parameters of impact can be determined after excavation. The craters were formed by multi-ton fragments of a type IA meteoroid with abundant silicate inclusions. Relative to the horizontal, the angle of infall was ～9°. Reflecting the low angle of infall, the crater field is elongated with apparent dimensions of 3 × 18.5 km. The largest craters are near the center of this ellipse. This suggests that when the parent meteoroid broke apart, the resulting fragments diverged from the original trajectory in inverse relation to their masses and did not undergo size sorting due to atmospheric deceleration. The major axis of the crater field as we know it extends along N63°E, but the azimuths of infall determined by excavation of Craters 9 and 10 are N83.5°E and N75.5°E, respectively. This suggests that the major axis of the crater field is not yet well determined. The three or four largest craters appear to have been formed by impacts that disrupted the projectiles, scattering fragments around the outsides of the craters and leaving no large masses within them; these are relatively symmetrical in shape. Other craters are elongated features with multi-ton masses preserved within them and no fragmentation products outside. There are two ways in which field research on the Campo del Cielo crater field is found to be useful. (1) Studies exist that have been used to interpret impact craters on planetary surfaces other than the Earth. This occurrence of a swarm of projectiles impacting at known angles and similar velocities into a uniform target material provides an excellent field site at which to test the applicability of those studies. (2) Individual craters at Campo del Cielo can yield the masses of the projectiles that formed them and their velocities, angles and azimuths of impact. From these data, there is a possibility to estimate parameters for the parent meteoroid at entry and, thus, learn enough about its orbit to judge whether or not it was compatible with an asteroidal origin. Preliminary indications are that it was. Campo del Cielo is a IA iron meteorite and Sikhote-Alin, an observed fall, is a IIB iron meteorite in Wasson's classification. The Sterlitamak iron, also an observed fall, is a medium octahedrite in the Prior-Hey classification. It would be interesting to compare their orbital parameters.     

The New Peruvian Meteorite Carancas: Mössbauer Spectroscopy and X-Ray Diffraction Studies

The Carancas meteorite fell on 15 September 2007 approximately 10 km south of Desaguadero, near Lake Titicaca, Peru, producing bright lights, clouds of dust in the sky and intense detonations. The Carancas meteorite is classified as a H4–5 ordinary chondrite with shock stage S3 and a degree of weathering W0. The Carancas meteorite is characterized by well defined chondrules composed either of olivine or pyroxene. The Mössbauer spectra show an overlapping of paramagnetic and magnetic phases. The spectra show two quadrupole doublets associated to olivine and pyroxene; and two magnetic sextets, associated with the primary phases kamacite/taenite and Troilite (Fe²⁺). Metal particles were extracted from the bulk powdered samples exhibit only kamacite and small amounts of the intergrowth tetrataenite/antitaenite. X-Ray diffractogram shows the primary phases olivine, pyroxene, troilite, kamacite, diopside and albite. Iron oxides has not been detected by Mössbauer spectroscopy or XRD as can be expected for a meteorite immediately recovered after its fall.     

Time-resolved studies of hypervelocity vertical impacts into porous particulate targets: Effects of projectile density on early-time coupling and crater growth

The depth and duration of energy and momentum coupling in an impact shapes the formation of the crater. The earliest stages of crater growth (when the projectile transfers its energy and momentum to the target) are unrecoverable when the event is described by late stage parameters, which collapse the initial conditions of the impact into a singular point in time and space. During the coupling phase, the details of the impact are mapped into the ejecta flow field. In this experimental study, we present new experimental and computational measurements of the ejecta distribution and crater growth extending from early times into main-stage ballistic flow for hypervelocity impacts over a range of projectile densities. Specifically, we assess the effect of projectile density on coupling depth and location in porous particulate (sand) targets. A non-invasive high-speed imaging technique is employed to capture the velocity of individual ejecta particles very early in the cratering event as a function of both time and launch position. These data reveal that the effects of early-stage coupling, such as non-constant ejection angles, manifest not only in early-time behavior but also extend to main-stage crater growth. Time-resolved comparisons with hydrocode calculations provide both benchmarking and insight into the parameters controlling the ejection process. Measurements of the launch position and metrics for the transient diameter to depth ratio as a function of time demonstrate non-proportional crater growth throughout much of excavation. Low-density projectiles couple closer to the surface, thereby leading to lower ejection angles and larger effective diameter to depth ratios. These results have implications for the ballistic emplacement of ejecta on planetary surfaces, and are essential to interpreting temporally resolved data from impact missions.     

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.     

Stony meteorite characterization by non‐destructive measurement of magnetic properties

Abstract— Four parameters of low‐field magnetic susceptibility (bulk value, frequency dependence, degree of anisotropy, and ellipsoid shape) have been determined for 321 stony meteorites from the National Collection of Canada. These parameters provide a basis for rapid, non‐destructive, and accurate meteorite classification as each meteorite class tends to have a distinct range of values. Chondrites show a clear trend of increasing bulk susceptibility from LL to L to H to E within the 3.6 to 5.6 logχ (in 10⁻⁹ m³/kg) range, reflecting increasing Fe‐Ni metal and Fe‐Ni sulfide content. Achondrite values range in logχ from 2.4 to 4.7 and primitive achondrites from 4.2 to 5.7. Frequency dependence is observed, using 19,000 Hz and 825 Hz, with variations in strength among meteorite classes and individual specimen dependence ranging from 1–25.6%. Degrees of anisotropy range from 1 to 53% with both oblate and prolate ellipsoids present. The aubrite class is marked by high degrees of anisotropy, low bulk magnetic susceptibility, and prolate fabric. Camel Donga is set apart from other eucrites, marked by higher bulk susceptibility, degree of anisotropy, and magnitude of oblate ellipsoid shape. The Shergotty, Nakhla, and Chassigny (SNC) meteorites show subclass distinction using frequency dependence and Chassigny is set apart with a relatively strong oblate fabric. The presence of both strong oblate and prolate fabrics among and within meteorite classes of chondritic and achondritic material points to a complex, multi‐mechanism origin for anisotropy, more so than previously thought, and likely dominated by impact processes in the later stages of stony parent body formation.     

Propagation of impact‐induced shock waves in porous sandstone using mesoscale modeling

Generation and propagation of shock waves by meteorite impact is significantly affected by material properties such as porosity, water content, and strength. The objective of this work was to quantify processes related to the shock‐induced compaction of pore space by numerical modeling, and compare the results with data obtained in the framework of the Multidisciplinary Experimental and Modeling Impact Research Network (MEMIN) impact experiments. We use mesoscale models resolving the collapse of individual pores to validate macroscopic (homogenized) approaches describing the bulk behavior of porous and water‐saturated materials in large‐scale models of crater formation, and to quantify localized shock amplification as a result of pore space crushing. We carried out a suite of numerical models of planar shock wave propagation through a well‐defined area (the “sample”) of porous and/or water‐saturated material. The porous sample is either represented by a homogeneous unit where porosity is treated as a state variable (macroscale model) and water content by an equation of state for mixed material (ANEOS) or by a defined number of individually resolved pores (mesoscale model). We varied porosity and water content and measured thermodynamic parameters such as shock wave velocity and particle velocity on meso‐ and macroscales in separate simulations. The mesoscale models provide additional data on the heterogeneous distribution of peak shock pressures as a consequence of the complex superposition of reflecting rarefaction waves and shock waves originating from the crushing of pores. We quantify the bulk effect of porosity, the reduction in shock pressure, in terms of Hugoniot data as a function of porosity, water content, and strength of a quartzite matrix. We find a good agreement between meso‐, macroscale models and Hugoniot data from shock experiments. We also propose a combination of a porosity compaction model (ε–α model) that was previously only used for porous materials and the ANEOS for water‐saturated quartzite (all pore space is filled with water) to describe the behavior of partially water‐saturated material during shock compression. Localized amplification of shock pressures results from pore collapse and can reach as much as four times the average shock pressure in the porous sample. This may explain the often observed localized high shock pressure phases next to more or less unshocked grains in impactites and meteorites.     

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.     

Characteristics of oceanic impact‐induced large water waves—Re‐evaluation of the tsunami hazard

Abstract— The potential hazard of a meteorite impact in the ocean is controversial with respect to the destructive power of generated large ocean waves (tsunamis). We used numerical modeling of hypervelocity impact to investigate the generation mechanism and the characteristics of the resulting waves up to a distance of 100–150 projectile radii. The wave signal is primarily controlled by the ratio between projectile diameter and water depth, and can be roughly classified into deep‐water and shallow‐water impacts. In the latter, the collapse of the crater rim results in a wave signal similar to solitary waves, which propagate and decay in agreement with shallow‐water wave theory. The much more likely scenario for an asteroid impact on Earth is a relatively small body (much smaller than the water depth) striking the deep sea. In this case, the collapse of the transient crater results in a significantly different and much more complex wave signal that is characterized by strong nonlinear behavior. We found that such waves decay much more rapidly than previously assumed and cannot be treated as long waves. For this reason, the shallow‐water theory is not applicable for the computation of wave propagation, and more complex models (full solution of the Boussinesq equations) are required.     

Ernst Julius Öpik's (1916) note on the theory of explosion cratering on the Moon's surface—The complex case of a long‐overlooked benchmark paper

High‐velocity impact as a common phenomenon in planetary evolution was ignored until well into the twentieth century, mostly because of inadequate understanding of cratering processes. An eight‐page note, published in Russian by the young Ernst Julius Öpik, a great Estonian astronomer, was among the key selenological papers, but due to the language barrier, it was barely known and mostly incorrectly cited. This particular paper is here intended to serve as an explanatory supplement to an English translation of Öpik's article, but also to document an early stage in our understanding of cratering. First, we outline the historical–biographical background of this benchmark paper, and second, a comprehensive discussion of its merits is presented, from past and present perspectives alike. In his theoretical research, Öpik analyzed the explosive formation of craters numerically, albeit in a very simple way. For the first time, he approximated relationships among minimal meteorite size, impact energy, and crater diameter; this scaling focused solely on the gravitational energy of excavating the crater (a “useful” working approach). This initial physical model, with a rational mechanical basis, was developed in a series of papers up to 1961. Öpik should certainly be viewed as the founder of the numerical simulation approach in planetary sciences. In addition, the present note also briefly describes Nikolai A. Morozov as a remarkable man, a forgotten Russian scientist and, surprisingly, the true initiator of Öpik's explosive impact theory. In fact, already between 1909 and 1911, Morozov probably was the first to consider conclusively that explosion craters would be circular, bowl‐shaped depressions even when formed under different impact angles.     

The age of Wolfe Creek meteorite crater (Kandimalal), Western Australia

Wolfe Creek crater lies in northwestern Australia at the edge of the Great Sandy Desert. Together with Meteor Crater, it is one of the two largest craters on Earth from which meteorite fragments have been recovered. The age of the impact is poorly constrained and unpublished data places the event at about 300,000 years ago. In comparison, Meteor Crater is well constrained by exposure dating. In this paper, we present new ages for Wolfe Creek Crater from exposure dating using the cosmogenic nuclides ¹⁰Be and ²⁶Al, together with optically stimulated luminescence ages (OSL) on sand from a site created by the impact. We also present a new topographic survey of the crater using photogrammetry. The exposure ages range from ~86 to 128 ka. The OSL ages indicate that the age of the impact is most likely to be ~120 ka with a maximum age of 137 ka. Considering the geomorphic setting, the most likely age of the crater is 120 ± 9 ka. Last, we review the age of Meteor Crater in Arizona. Changes in production rates and scaling factors since the original dating work revise the impact age to 61.1 ± 4.8 ka, or ~20% older than previously reported.