Martian cratering 8: Isochron refinement and the chronology of Mars

This paper reviews and refines the technique of dating martian surfaces by using impact-crater isochrons (defined as size distributions of impact craters on undisturbed martian surfaces of specified ages). In the 1970s, this system identified not only abundant ancient martian volcanic surfaces, but also extensive lava plains with ages of a few 10⁸ y-old; this dating was initially controversial but confirmed in the 1980s and 90s by martian meteorites. The present update utilizes updated estimates of the Mars/Moon cratering ratio (the most important calibration factor), improves treatment of gravity and impact velocity scaling effects, combines aspects of the crater size distribution data from earlier work by both Neukum and Hartmann, and for the first time applies a correction for loss of small meteoroids in the martian atmosphere from Popova et al. (2003, Meteorit. Planet. Sci. 38, 905-925). The updated isochrons are not radically different from the previous “2002 iteration” but fit observed data better and give somewhat older model ages for features dated from small craters (diameter D<100 m). Crater counts from young lava flows in various areas give good fits to the new isochrons over as much as 3 orders of magnitude in D, confirming the general isochron shape and giving crater retention ages in the range of some 10⁶ to some 10⁸ y, interpreted as lava flow ages. More complex, older units are also discussed. Uncertainties are greatest if only small craters (D?100 m) are used. Suggestions by other workers of gross uncertainties, due to local secondary craters and deposition/exhumation, are discussed; they do not refute our conclusions of significant volcanic, fluvial, and other geologic activity in the last few percent of martian geologic time or the importance of cratering as a tool for studying processes such as exhumation. Indeed, crater count data suggest certain very recent episodes of deposition, exhumation, and ice flow, possibly associated with obliquity cycles of ∼¹⁰⁷ y timescale. Evidence from ancient surfaces suggests higher rates of volcanism, fluvial activity, glaciation, and other processes in Noachian/Hesperian time than in Amazonian time.     

Crater densities on the satellites of Mars

The density of craters larger than 1 km in diameter has been determined for the entire surface of Phobos, and half that of Deimos. Densities of craters as small as 10 m on Phobos and 5 m on Deimos have been measured for small areas of the satellites. On both objects, crater densities are similar and yield plots which have slopes close to -1.9 on both incremental and cumulative log-log graphs. These densities are close to those expected to obtain under equilibrium conditions. They are also near the maximum predicted, based on the fragmentation lifetimes of the two objects: that is, the densities are near to the maximum possible before such objects are likely to suffer an impact severe enough to disrupt them. While the observed crater densities cannot be converted to absolute ages in any rigorous fashion, they can be understood if the flux at Mars has been similar to that at the Moon and if the surfaces that we see today generally date back to the end of the period of heavy bombardment some 4 billion years ago. It is extremely unlikely that the surfaces are younger than 1 billion years. There are no large areas on Phobos for which crater densities differ by more than a factor of 3 from the average.     

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.     

Small impact craters in the lunar regolith — Their morphologies,relative ages,and rates of formation

Apparently, there are two types of size-frequency distributions of small lunar craters (1–100 m across): (1) crater production distributions for which the cumulative frequency of craters is an inverse function of diameter to power near 2.8, and (2) steady-state distributions for which the cumulative frequency of craters is inversely proportional to the square of their diameters. According to theory, cumulative frequencies of craters in each morphologic category within the steady-state should also be an inverse function of the square of their diameters. Some data on frequency distribution of craters by morphologic types are approximately consistent with theory, whereas other data are inconsistent with theory.A flux of crater producing objects can be inferred from size-frequency distributions of small craters on the flanks and ejecta of craters of known age. Crater frequency distributions and data on the craters Tycho, North Ray, Cone, and South Ray, when compared with the flux of objects measured by the Apollo Passive Seismometer, suggest that the flux of objects has been relatively constant over the last 100 m.y. (within 1/3 to 3 times of the flux estimated for Tycho).Steady-state frequency distributions for craters in several morphologic categories formed the basis for estimating the relative ages of craters and surfaces in a system used during the Apollo landing site mapping program of the U.S. Geological Survey. The relative ages in this system are converted to model absolute ages that have a rather broad range of values. The range of values of the absolute ages are between about 1/3 to 3 times the assigned model absolute age.     

Distribution and relations of 4- to 10-km-diameter craters to global geologic units of Mars

By correlating the 1:25,000,000 geologic map of Mars of Scott and Carr (1977) with 4- to 10-km-diameter crater density data from Mariner 9 images, the average crater density for 23 of the equatorial geologic-geomorphic units on Mars was computed. The correlation of these two data sets was accomplished by digitizing both the crater density data and geologic map at the same scale and by comparing them in a computer. This technique assigns the crater density value found in the corresponding location on the geologic data set to a discrete computer file assigned each of the 23 geologic units. By averaging the crater density values accumulated in each file, an “average” crater density for each geologic unit was obtained. Condit believes these average crater density values are accurate indicators of the relative age of the geologic units considered. The statistical validity of these average values is strongest for the geologic units of the largest areal extent. The relative ages as obtained from the average crater density values for the seven largest geologic units, from youngest to oldest, are: Tharsis volcanic material, 21 ± 4 craters/10⁶km²; smooth plains material, 57 ± 14 craters/10⁶km²; rolling plains material, 66 ± 16 craters/10⁶km²; plains materials, 80 ± 17 craters/10⁶km²; ridged plains material, 128 ± 25 craters/10⁶km²; hilly and cratered material, 137 ± 38 craters/10⁶km²; and cratered plateau material, 138 ± 27 craters/10⁶km².     

Crater modification and geologic activity in Enceladus' heavily cratered plains: Evidence from the impact crater distribution

Although we can observe current activity on Saturn's satellite Enceladus with Cassini, insight into past activity is best achieved (for now) through studying the impact crater distributions. Furthermore, approximation of terrain ages can only be attained through calculations using crater densities and estimations of impact rates in the saturnian system. Here we focus on what the impact crater distribution in Enceladus' heavily cratered plains can tell us about Enceladus' geologic history. We use Cassini ISS images to count craters in the heavily cratered plains on Enceladus, along with Rhea, Dione, Tethys and Mimas as references, to develop and compare their size-frequency distributions. Comparisons of our counts show that Enceladus' cratered plains distribution is unique in that it appears to have a relative deficiency of craters for diameters ?2 km and ?6 km compared to the other satellites' heavily cratered plains. Our data also indicates that the impact crater density within the cratered plains changes with latitude. Specifically, both the north and south mid-latitude regions have approximately three times higher density than the equatorial region. We hypothesize that the “missing” small and large craters in Enceladus' cratered plains is due to a combination of viscous relaxation of the larger craters, and burial of the relaxed large craters and small craters by south polar plume and possibly E-ring material. We also conclude that the spatial density distribution is not consistent with recent polar wander.     

Morphology of lunar craters: A test of lunar erosional models

This paper is a test of published theoretical and experimental studies of crater erosion by micrometeorite bombardment which predict systematic variations in the morphology of lunar craters as a function of crater diameter and crater age. Numerical, ranking-type degradation classifications indicate that the craters on Mare Imbrium and Mare Tranquillitatus confirm these predictions by showing a systematic increase in degradation with decreasing diameter for craters smaller than a few kilometers in diameter but larger than the equilibrium diameter, and by showing fixed proportions of fresh, moderately degraded and very degraded craters under equilibrium conditions. Furthermore, the relative ages of the two mare surfaces may be determined using a diameter/mean-degradation-number curve. These determinations of relative age and process of crater erosion are both essentially independent of the traditionally studied crater diameter/frequency relationships. Morphologies of terra craters near Mare Humorum suggest a young, non-equilibrium crater population superposed on a perimordial population with about equilibrium proportions of fresh, moderately degraded and very degraded craters. The primordial population has been modified by pre-Imbrian or early Imbrian deposition of blanketing deposits. A comparative study of several crater degradation classifications indicates that all are essentially interchangeable.     

The solar system cratering record: Voyager 2 results at Uranus and implications for the origin of impacting objects

The cratering record at Uranus shows two different crater populations of different ages. The old crater population occurs on the heavily cratered surfaces of Oberon, Umbriel, and Miranda, while the younger one is found on Titania, Ariel and the resurfaced areas of Miranda. Since only the young population occurs on Titania, this satellite must have experienced a global resurfacing event which obliterated the older population prior to the impact of objects causing the younger one. The old crater population is characterized by an abundance of large craters and a relative paucity of small ones. The young crater population, however, has an abundance of small craters and a paucity of large ones relative to the old population. Furthermore, the abundance of small craters and the paucity of large craters increases with decreasing density. This change in the size distribution is consistent with a population of impactors that evolved with time by mutual collision, and therefore was probably in planetocentric orbits. In fact, both crater populations may be the result of accretional remnants in planetocentric orbits that evolved with time by mutual collisions. If so, then the higher crater density on Miranda compared to Oberon and Umbriel suggests that both Oberon and Umbriel were also resurfaced early in their histories.A comparison of the Solar System cratering record from Mercury to Uranus (19 AU) shows different crater populations at different locations in the Solar System. Computer simulations using a modified Holsapple-Schmidt crater scaling and short-period comet impact velocities to recover the projectile diameters from the cratering record produce different projectile populations in different parts of the Solar System. Furthermore, adjusting the Jovian crater curve to match that in the inner Solar System requires differences in the impact velocities that are unrealistic for objects in heliocentric orbits. These results suggest that the Solar System cratering record cannot be explained by a single family of objects in heliocentric orbits, e.g., comets. One possible explanation is that the cratering record is the result of different families of objects (possibly accretional remnants) indigenous to that region of the Solar System in which the different crater populations are found. Thus, in the inner Solar System, the impactors responsible for heavy bombardment were in heliocentric orbits with semimajor axes less than 3 AU. In the outer Solar System, they may have been in planetocentric orbits around each of the Jovian planets.     

Deducing the age of the dense Venus atmosphere

We show how crater size-density counts may be used to help constrain the history of the Venus atmosphere, based on the predictions of simple but reasonable models for crater production, surface erosion, and the effects of atmospheric drag and breakup on incident meteors in the Venus atmosphere. If the atmosphere is old, we may also be able to determine the importance of breakup as a mechanism for destroying incident meteors in a dense fluid. In particular, if the atmosphere is young, the old (uneroded) surfaces will have crater densities upward of 10^?4 km^?2 and a ratio of small (4 km) craters to large (128 km) craters near 10³. If the atmosphere is old and the breakup mechanism is dominant, absolute crater densities on Venus surfaces will be diminished by several orders of magnitude relative to the young atmosphere case. If atmospheric drag is dominant and the atmosphere is old, the absolute crater density will be lowered by perhaps an order of magnitude relative to the young atmosphere case, and the ratio of small to large craters will be reduced to a value near 10^1.5 according to the models. The comparison of crater populations on young, as well as old, surfaces on Venus can help in distinguishing the young and old atmosphere scenarios, especially since the situation may be complicated by currently undetermined erosional and tectonic processes. Once a large fraction of Venus surface has been imaged at kilometer resolution, as the VOIR project promises to do, it could be possible to make an early determination of the age of the Venus atmosphere.     

Automatic detection of sub-km craters in high resolution planetary images

Impact craters are among the most studied geomorphic planetary features because they yield information about the past geological processes and provide a tool for measuring relative ages of observed geologic formations. Surveying impact craters is an important task which traditionally has been achieved by means of visual inspection of images. The shear number of smaller craters present in high resolution images makes visual counting of such craters impractical. In this paper we present a method that brings together a novel, efficient crater identification algorithm with a data processing pipeline; together they enable a fully automatic detection of sub-km craters in large panchromatic images. The technical details of the method are described and its performance is evaluated using a large, 12.5 m/pixel image centered on the Nanedi Valles on Mars. The detection percentage of the method is ∼70%. The system detects over 35,000 craters in this image; average crater density is , but localized spots of much higher crater density are present. The method is designed to produce “million craters” global catalogs of sub-km craters on Mars and other planets wherever high resolution images are available. Such catalogs could be utilized for deriving high spatial resolution and high temporal precision stratigraphy on regional or even planetary scale.     

Viscous relaxation of impact craters on icy planetary surfaces: Determination of viscosity variation with depth

Spacecraft images show that the icy Galilean satellites have surfaces with very low topographic relief. Impact craters on Ganymede and Callisto are anomalously shallow and are characterized by sharp well-defined rims and domed floors. These morphological characteristics can be explained by viscous relaxation of topography on an icy crust in which the viscosity is uniform or decreases with depth. Under these conditions, large craters relax more rapidly than small craters, therefore explaining a possible underabundance of large craters. Viscous relaxation on an icy crust that is thin compared to the crater diameter or on a thick icy crust in which viscosity increases with depth could not produce this crater morphology and would result in the more rapid relaxation of small craters rather than large craters. The results of this study suggest that more detailed analysis of relaxing impact crater morphology may resolve the rate of viscosity decrease with depth and so provide evidence on the interior thermal evolution of icy planetary bodies.     

Can spatial statistics help decipher impact crater saturation?

Impact crater saturation equilibrium is a state where a surface is so densely cratered that a new crater cannot form without removing older craters and the observed crater density is in (quasi-)equilibrium. Whether densely cratered surfaces throughout the solar system are saturated for large, kilometer-sized craters has been debated for decades. This work explores if spatial statistics can provide insight if these crater distributions are in saturation equilibrium. The supposition is that crater distributions become more spatially uniform (more evenly spaced) as they reach saturation (Squyres et al. 1997 ). A numerical simulation of crater saturation is combined with observations of cratered terrains throughout the solar system to assess the utility of using spatial statistics. The numerical simulations examine spatial statistics and saturation equilibrium for crater distributions for various input population size-frequency distribution (SFD) slopes, along with a range in the effective crater erasure size, effectiveness of smaller craters erasing the rims of larger craters, and the amount of rim needed to recognize a crater. Simulations show that saturated terrains do become more spatially uniform, and that the degree of uniformity appears to be most dependent on the input SFD slope. When simulation results are compared to observed crater distributions, I find that large, kilometer-sized craters on densely cratered terrains throughout the solar system are likely in saturation equilibrium.     

The record of impact cratering on the great volcanic shields of the Tharsis region of Mars

Mariner 9 images of the four great volcanic shields of the Tharsis region of Mars show many circular craters ranging in diameter from 100mm to 20 km. Previous attempts to date the volcanoes from their apparent impact crater densities yielded a range of results. The principal difficulty is sorting volcanic from impact craters for diameters $\text{?1}\text{km}$. Many of the observed craters are aligned in prominent linear and concentric patterns suggestive of volcanic origin. In this paper an attempt is made to date areas of shield surface, covered with high resolution images using only scattered small (?1 km) craters of probable impact origin. Craters of apparent volcanic origin are systematically excluded from the dating counts.The common measure of age, deduced for all surfaces studied, is a calculated “crater age” F′ defined as the number of craters equal to or larger than 1 km in diameter per 10⁶km². The conclusions reached from comparing surface ages and their geological settings are: (1) Lava flow terrain surfaces with ages, F′, from 180 to 490 are seen on the four great volcanoes. Summit surfaces of similar ages, F′ = 360 to 420, occur on the rims of calderas of Arsia Mons, Pavonis Mons, and Olympus Mons. The summit of Ascraeus Mons is possibly younger; F′ is calculated to be 180 for the single area which could be dated. (2) One considerably younger surface, F′ < 110, is seen on the floor of Arsia Mon's summit caldera. (3) Nearly crater free lava flow terrain surfaces seen on Olympus Mons are estimated to be less than half the age of a summit surface. The summit caldera floor is similarly young. (4) The pattern of surface ages on the volcanoes suggests that their eruption patterns are similar to those of Hawaiian basaltic shields. The youngest surfaces seem concentrated on the mid-to-lower flanks and within the summit calderas. (5) The presently imaged sample of shield surface, though incomplete, clearly shows a broad range of ages on three volcanoes—Olympus, Arsia, and Pavonis Mons.Estimated absolute ages of impact dated surfaces are obtained from two previously published estimates of the history of flux of impacting bodies on Mars. The estimated ranges of age for the observed crater populations are 0.5 to 1.2b.y. and 0.07 to 0.2b.y. Areas which are almost certainly younger, less than 0.5 or 0.07b.y., are also seen. The spans of surface age derived for the great shields are minimum estimates of their active lifetimes, apparently very long compared to those of terrestrial volcanoes.     

A mechanism for the production of lunar crater rays

Observations of high resolution photographs of part of one of the prominent rays of the lunar crater Copernicus show that there is a concentration of small bright rayed and haloed craters within the ray. These craters contribute to the overall ray brightness; they have been measured and their surface distribution has been mapped. Sixty-two percent of the bright craters can be identified from study of high resolution photographs as concentric impact craters. These craters contain in their ejecta blankets, rocks from the lunar substrate that are brighter than the adjacent mare surface. It is concluded that the brightness of the large ray from the crater Copernicus is due to the composite effect of many small concentric impact craters with rocky ejecta blankets. If this is the dominant mechanism for the production of other rays from Copernicus and other large lunar craters, then rays may not contain significant amounts of ejecta from the central crater or from large secondary craters. They may in fact only reflect local excavation of mare substrate material by myriads of small secondary or tertiary impact craters.     

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.     

A meteorite impact crater field in eastern Bavaria? A preliminary report

Abstract— Numerous circular depressions north of Burghausen in eastern Bavaria, with diameters ranging from meters to tens of meters in size and dispersed over an area of at least 11 times 7 km, are suspected to have an extraterrestrial origin since they resemble other small meteorite impact craters. The depressions are bowl‐shaped, have high circularity and a characteristic rim. Most of them were formed in unconsolidated glacial gravels and pebbles intermixed with fine‐grained sand and clay. Magnetic investigations reveal weak anomalies with amplitudes of less than ±10 nano Tesla (nT). In some cases, the origins of the anomalies are suspected to be due to human activity within the structures. So far, no traces of meteoritic material have been detected. An evident archaeological or local geological explanation for the origin of the craters does not exist. A World War I and II explosive origin can be excluded since trees with ages exceeding 100 years can be found in some craters. One crater was described in 1909. Carbon‐14 dating of charcoal found in one crater yielded an age of 1790 ± 60 years. Hence, a formation by meteorite impacts that occurred in Celtic or early medieval times should be considered. A systematic archaeological excavation of some structures and an intensified search for traces of meteoritic material are planned.     

Rays and secondary craters of Tycho

The large, fresh crater Tycho in the nearside lunar highlands has an extensive system of bright rays covering approximately 560,000 km², containing dense clusters of secondary craters. Examination of crater densities in several clusters shows that Tycho produced almost 10⁶ secondary craters larger than 63 m diameter. This is a lower limit, because small crater densities are reduced, most likely by mass wasting. We estimate a crater erasure rate of 2-6 cm/Myr, varying with crater size, and consistent with previous results. This process has removed many small craters, and it is probable that the original number of secondary craters formed by Tycho was higher. Also, we can only identify distant secondaries of Tycho where they occur in bright rays. Craters on Mars and Europa also formed large numbers of secondaries, but under possibly ideal conditions for spallation as a mechanism to produce high-velocity ejecta fragments. The results from Tycho show that large numbers of such fragments can be produced even from impact into a heavily fragmented target on which spallation is expected to be less important.     

Transport and emplacement of crater and basin deposits

Material is ejected from impact craters in ballastic trajectories; it impacts first near the crater rim and then at progressively greater ranges. Ejecta from craters smaller than approximately 1 km is laid predominantly on top of the surrounding surface. With increasing crater size, however, more and more surrounding surface will be penetrated by secondary cratering action and these preexisting materials will be mixed with primary crater ejecta. Ejecta from large craters and especially basin forming events not only excavate preexisting, local materials, but also are capable of moving large amounts of material away from the crater. Thus mixing and lateral transport give rise to continuous deposits that contain materials from within and outside the primary crater. As a consequence ejecta of basins and large highland craters have eroded and mixed highland materials throughout geologic time and deposited them in depressions inside and between older crater structures.Because lunar mare surfaces contain few large craters, the mare regolith is built up by successive layers of predominantly primary ejecta. In contrast, the lunar highlands are dominated by the effects of large scale craters formed early in lunar history. These effects lead to thick fragmental deposits which are a mixture of primary crater material and local components. These deposits may also properly be named regolith though the term has been traditionally applied only to the relatively thin fine grained surficial deposit on mare and highland terranes generated during the past few billion year. We believe that the surficial highland regolith - generated over long periods of time - rests on massive fragmental units that have been produced during the early lunar history.     

Monte Carlo calculations of lunar Regolith thickness distributions

Regolith thickness distributions associated with crater populations observed on selected maria surfaces have been calculated using a Monte Carlo computer technique. The calculations assume that the crater type produced and the volume of debris ejected and added to the growing regolith depends on the ratio of crater diameter and regolith thickness present at the time and place of formation of each crater. Calculated thickness distributions obtained are in agreement with those estimated using a previously described statistical method based on the morphology of small lunar craters. Additionally, the Monte Carlo calculations accurately predict the size frequency distributions of the same types of small, fresh lunar craters used in the statistical method. The model employed is therefore realistic. Furthermore, the model calculations presented are shown to have value (a) in predicting the thickness of the regolith from crater populations at various lunar sites, (b) relative dating applications in which crater populations are compared, and (c) in interpreting the origin and history of regolith deposits at specific locations.