The volcanic history of Mars: High-resolution crater-based studies of the calderas of 20 volcanoes

Determining absolute surface ages for bodies in the Solar System is, at present, only possible for Earth and Moon with radiometric dating for both bodies and biologic proxies such as fossils for Earth. Relative ages through cratering statistics are recognized as one of the most reliable proxies for relative ages, calibrated by lunar geologic mapping and Apollo program sample returns. In this work, we have utilized the Mars Reconnaissance Orbiter’s ConTeXt Camera’s images which provide the highest resolution wide-scale coverage of Mars to systematically crater-age-date the calderas of 20 of Mars’ largest volcanoes in order to constrain the length of time over which these volcanoes - and major volcanic activity on the planet, by extension - were active. This constitutes the largest uniform and comprehensive research on these features to date, eliminating unknown uncertainties by multiple researchers analyzing different volcanoes with varied data and methods. We confirm previous results that Mars has had active volcanism throughout most of its history although it varied spatially and temporally, with the latest large-scale caldera activity ending approximately 150 ma in the Tharsis region. We find a transition from explosive to effusive eruption style occurring in the Hesperian, at approximately 3.5 Ga ago, though different regions of the planet transitioned at different times. Since we were statistically complete in our crater counts to sizes as small as ∼60 m in most cases, we also used our results to study the importance of secondary cratering and its effects on crater size-frequency distributions within the small regions of volcanic calderas. We found that there is no “golden rule” for the diameters secondaries become important in crater counts of martian surfaces, with one volcano showing a classic field of secondaries ∼2 crater diameters from the center of its primary but not affecting the size-frequency distribution, and another clearly showing an influence but from no obvious primary.     

Evidence for Hesperian impact-induced hydrothermalism on Mars

Several hydrated silicate deposits on Mars are observed within craters and are interpreted as excavated Noachian material. Toro crater (71.8°E, 17.0°N), located on the northern edge of the Syrtis Major Volcanic Plains, shows spectral and morphologic evidence of impact-induced hydrothermal activity. Spectroscopic observations were used to identify extensive hydrated silicate deposits, including prehnite, chlorites, smectites, and opaline material, a suite of phases that frequently results from hydrothermal alteration in terrestrial craters and also expected on Mars from geochemical modeling of hydrothermal environments. When combined with altimetry and high-resolution imaging data, these deposits appear associated predominantly with the central uplift and with portions of the northern part of the crater floor. Detailed geologic mapping of these deposits reveals geomorphic features that are consistent with hydrothermal activity that followed the impact event, including vent-like and conical mound structures, and a complex network of tectonic structures caused by fluid interactions such as fractures and joints. The crater age has been calculated from the cumulative crater size-frequency distributions and is found to be Early Hesperian. The evidence presented here provides support for impact-induced hydrothermal activity in Toro crater, that extends phyllosilicate formation processes beyond the Noachian era.     

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.     

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.     

Millochau crater, Mars: Infilling and erosion of an ancient highland impact crater

The geology and stratigraphy of Millochau crater (21.4° S, 275° W), located in the highlands of Tyrrhena Terra, Mars, are documented through geomorphic analyses and geologic mapping. Crater size-frequency distributions and superposition relationships are used to constrain relative ages of geologic units and determine the timing and duration of the geologic processes that modified Millochau rim materials and emplaced deposits on Millochau's floor. Crater size-frequency distributions show a Middle Noachian age for rim materials and Middle Noachian to Early Hesperian ages for most of the interior deposits. Valley networks and gullies incised within Millochau's rim materials and interior wall, respectively, indicate fluvial activity was an important erosional process. Millochau contains an interior plateau, offset northeast of Millochau's center, which rises up to 400 m above the surrounding crater floor and slopes downward to the south and west. Layers exposed along the northern and eastern scarp boundaries of the plateau are tens to hundreds of meters thick and laterally continuous in MOC images. These layers suggest most materials within Millochau were emplaced by sedimentary processes (e.g., fluvial or eolian), with the potential for lacustrine deposition in shallow transient bodies of water and contributions of volcanic airfall. Mass wasting may have also contributed significant quantities of material to Millochau's interior, especially to the deposits surrounding the plateau. Superposition relationships combined with impact crater statistics indicate that most deposition and erosion of Millochau's interior deposits is ancient, which implies that fluvial activity in this part of Tyrrhena Terra is much older than in the eastern Hellas region. Eolian processes mobilized sediment to form complicated patterns of long- and short-wavelength dunes, whose emplacement is controlled by local topography. These deposits are some of the youngest within Millochau (Amazonian) and eolian modification may be ongoing.     

A description of the crater Tsiolkovsky on the lunar far side

Although researchers in the last decade have been primarily concerned with the exotic findings of the more distant planets and moons in our solar system, as given by the Voyager series, there is still much work to be done on our nearer neighbours, including the Moon. This paper summarizes some basic age dating of a portion of the lunar surface, namely the mare in the crater Tsiolkovsky on the lunar far side.Using the Apollo 15 panoramic camera photographs, the cumulative crater frequency (N km^-2) relative to crater diameter (D) distribution has been obtained for the mare in the crater Tsiolkovsky. The diameter size range sampled was 0.07 km < D < 1 km. A total of 12 604 craters were counted and their average apparent diameters measured. There were 85 sample areas on the mare surface which were chosen at random, after exclusion of blanketed, volcanic or secondary cratered areas. It was found that a large proportion of the crater floor contains endogenic features, especially volcanic vents at approximately D = 0.3 km. An additional 7 areas of interest were also examined in detail for comparison with areas of purely primary impact craters. Evidence for up to 8 lava floodings can be detected from the size-frequency distributions although no visual data, e.g., flow lobes, can be seen on the mare surface.The total size-frequency distribution for all the areas is coincident with Neukum et al. (1975a and b) Calibration Distribution in the size range 0.25 km < D < 1 km which is at the smallest crater diameters that they obtained. Neukum et al. (1975a and b) give their distribution as a polynomial of 7th degree. However, in this present study a variation is indicated in the steepening of the curve for D < 0.1 km.The results also approximate (but only for D < 0.6 km) the distribution obtained by Shoemaker et al. (1970) in the range 100 m < D < 3 km where N ~ D ^-2.9. The best fit line reached for the data given here is N ~ D ^-2.682.Comparison of the distribution with plots for the maria at Apollo 11, 12, and 15 landing sites show that Tsiolkovsky mare is 3.51 ± 0.1 × 10⁹ yr old. This agrees with other workers (see Gornitz, 1973) who place it between Mare Tranquillitatis (Apollo 11 radiometric dating: 3.5 to 3.9 aeons) and Oceanus Procellarum (Apollo 12: 3.5 to 3.4 aeons). There are no rock samples from Tsiolkovsky to given an absolute age.This places Tsiolkovsky mare within the weighted mean of the age range (1.0 to 4.3 × 10⁹ yr old) of the maria on the Moon. From this it can be concluded that the processes producing the vast basalt outpourings seen on the Moon's face apply for the far side also and that there is a linking factor for the whole Moon.     

The geology of the Viking Lander 2 site revisited

Reevaluating the geologic history of the prior Mars landing sites provides important ground truth for recent and ongoing orbital missions. At the Viking 2 Lander (VL2) site, topographic measurements of relict landforms indicate that at least 100 m of sedimentary mantle material has been stripped away. The observed paucity of impact craters <100 m in diameter suggests that resurfacing processes (likely in the form of the recent deposition and removal of thin 1-10 m mantle layers) continue up to the present. A dearth of craters in the 100-500 m diameter range, however, also necessitates erosion of a thicker mantle layer. Partially inverted chains of secondary craters from nearby Mie Crater indicate that the mantle was already in place when the impact occurred. The density of craters superposed on Mie ejecta is consistent with a Late Hesperian age and provides a minimum age constraint for the mantle's emplacement. The thermophysical properties of the surface around VL2 as observed with Thermal Emission Imaging System (THEMIS) data indicate that the landing site occurs in an intracrater region that may typify mid to high northern latitude sites. Elevated thermal inertias of a pedestal crater superposed atop a larger pedestal crater suggest that rocky or indurated material can be created by impacts into sedimentary targets. Rock abundances at VL2 are consistent with the addition of impact-emplaced material from the missing small impact crater population documented in this study. Thus, the VL2 site may be a reasonable proxy for the landscape expected at the upcoming Phoenix Lander site.     

Nonuniform cratering of the Moon and a revised crater chronology of the inner Solar System

We model the cratering of the Moon and terrestrial planets from the present knowledge of the orbital and size distribution of asteroids and comets in the inner Solar System, in order to refine the crater chronology method. Impact occurrences, locations, velocities and incidence angles are calculated semi-analytically, and scaling laws are used to convert impactor sizes into crater sizes. Our approach is generalizable to other moons or planets. The lunar cratering rate varies with both latitude and longitude: with respect to the global average, it is about 25% lower at (±65°N, 90°E) and larger by the same amount at the apex of motion (0°N, 90°W) for the present Earth-Moon separation. The measured size-frequency distributions of lunar craters are reconciled with the observed population of near-Earth objects under the assumption that craters smaller than a few kilometers in diameter form in a porous megaregolith. Varying depths of this megaregolith between the mare and highlands is a plausible partial explanation for differences in previously reported measured size-frequency distributions. We give a revised analytical relationship between the number of craters and the age of a lunar surface. For the inner planets, expected size-frequency crater distributions are calculated that account for differences in impact conditions, and the age of a few key geologic units is given. We estimate the Orientale and Caloris basins to be 3.73 Ga old, and the surface of Venus to be 240 Ma old. The terrestrial cratering record is consistent with the revised chronology and a constant impact rate over the last 400 Ma. Better knowledge of the orbital dynamics, crater scaling laws and megaregolith properties are needed to confidently assess the net uncertainty of the model ages that result from the combination of numerous steps, from the observation of asteroids to the formation of craters. Our model may be inaccurate for periods prior to 3.5 Ga because of a different impactor population, or for craters smaller than a few kilometers on Mars and Mercury, due to the presence of subsurface ice and to the abundance of large secondaries, respectively. Standard parameter values allow for the first time to naturally reproduce both the size distribution and absolute number of lunar craters up to 3.5 Ga ago, and give self-consistent estimates of the planetary cratering rates relative to the Moon.     

Eminescu impact structure: Insight into the transition from complex crater to peak-ring basin on Mercury

Peak-ring basins represent an impact-crater morphology that is transitional between complex craters with central peaks and large multi-ring basins. Therefore, they can provide insight into the scale dependence of the impact process. Here the transition with increasing crater diameter from complex craters to peak-ring basins on Mercury is assessed through a detailed analysis of Eminescu, a geologically recent and well-preserved peak-ring basin. Eminescu has a diameter (∼125 km) close to the minimum for such crater forms and is thus representative of the transition. Impact crater size-frequency distributions and faint rays indicate that Eminescu is Kuiperian in age, geologically younger than most other basins on Mercury. Geologic mapping of basin interior units indicates a distinction between smooth plains and peak-ring units. Our mapping and crater retention ages favor plains formation by impact melt rather than post-impact volcanism, but a volcanic origin for the plains cannot be excluded if the time interval between basin formation and volcanic emplacement was less than the uncertainty in relative ages. The high-albedo peak ring of Eminescu is composed of bright crater-floor deposits (BCFDs, a distinct crustal unit seen elsewhere on Mercury) exposed by the impact. We use our observations to assess predictions of peak-ring formation models. We interpret the characteristics of Eminescu as consistent with basin formation models in which a melt cavity forms during the impact formation of craters at the transition to peak ring morphologies. We suggest that the smooth plains were emplaced via impact melt expulsion from the central melt cavity during uplift of a peak ring composed of BCFD-type material. In this scenario the ringed cluster of peaks resulted from the early development of the melt cavity, which modified the central uplift zone.     

Catalogue of large alluvial fans in martian impact craters

We conducted a systematic, global survey using Thermal Emission Imaging System Infrared (THEMIS IR) coverage (∼100 m/pixel) to search for large alluvial fans in impact craters on Mars. Our survey has focused on large fans (apron areas greater than ∼40 km², usually located in craters greater than 20 km in diameter) due to the resolution of the THEMIS images and Mars Orbiter Laser Altimeter (MOLA) coverage. We find that the host craters are found to have a distinctive diameter range from 30-150 km. The fans generally cluster in three geographic areas—southern Margaritifer Terra, southwestern Terra Sabaea, and southwestern Tyrrhena Terra, however several outliers do exist. The alluvial fans do not form in a particular orientation along the crater rim nor are they associated with the location of current high rim topography. Fan area magnitude and variability increase with crater diameter while fan concavity magnitude and variability increase with decreasing crater diameter. Smaller fan aprons in general have higher, more variable concavity. The source of the water forming these fans is uncertain given the challenges of accommodating the global distribution pattern and formation patterns within the craters.     

Machine cataloging of impact craters on Mars

This study presents an automated system for cataloging impact craters using the MOLA 128 pixels/degree digital elevation model of Mars. Craters are detected by a two-step algorithm that first identifies round and symmetric topographic depressions as crater candidates and then selects craters using a machine-learning technique. The system is robust with respect to surface types; craters are identified with similar accuracy from all different types of martian surfaces without adjusting input parameters. By using a large training set in its final selection step, the system produces virtually no false detections. Finally, the system provides a seamless integration of crater detection with its characterization. Of particular interest is the ability of our algorithm to calculate crater depths. The system is described and its application is demonstrated on eight large sites representing all major types of martian surfaces. An evaluation of its performance and prospects for its utilization for global surveys are given by means of detailed comparison of obtained results to the manually-derived Catalog of Large Martian Impact Craters. We use the results from the test sites to construct local depth-diameter relationships based on a large number of craters. In general, obtained relationships are in agreement with what was inferred on the basis of manual measurements. However, we have found that, in Terra Cimmeria, the depth/diameter ratio has an abrupt decrease at ∼38°S regardless of crater size. If shallowing of craters is attributed to presence of sub-surface ice, a sudden change in its spatial distribution is suggested by our findings.     

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.     

Mercury crater statistics from MESSENGER flybys: Implications for stratigraphy and resurfacing history

The primary crater population on Mercury has been modified by volcanism and secondary craters. Two phases of volcanism are recognized. One volcanic episode that produced widespread intercrater plains occurred during the period of the Late Heavy Bombardment and markedly altered the surface in many areas. The second episode is typified by the smooth plains interior and exterior to the Caloris basin, both of which have a different crater size-frequency distribution than the intercrater plains, consistent with a cratering record dominated by a younger population of impactors. These two phases may have overlapped as parts of a continuous period of volcanism during which the volcanic flux tended to decrease with time. The youngest age of smooth plains volcanism cannot yet be determined, but at least small expanses of plains are substantially younger than the plains associated with the Caloris basin. The spatial and temporal variations of volcanic resurfacing events can be used to reconstruct Mercury's geologic history from images and compositional and topographic data to be acquired during the orbital phase of the MESSENGER mission.     

Dating very young planetary surfaces from crater statistics: A review of issues and challenges

Determining the ages of young planetary surfaces relies on using populations of small, often sub-km diameter impact craters due to the higher frequency at which they form. Smaller craters however can be less reliable for estimating ages as their size-frequency distribution is more susceptible to alteration with debate as to whether they should be used at all. With the current plethora of meter-scale resolution images acquired of the lunar and Martian surfaces, small craters have been widely used to derive model ages to establish the temporal relation of recent geologic events. In this review paper, we discuss the many factors that make smaller craters particularly challenging to use and should be taken into consideration when crater counts are confined to small crater diameters. Establishing confidence in a model age ultimately requires an understanding of the geologic context of the surface being dated as reliability can vary considerably and limitations of the dating technique should be considered in applying ages to any geologic interpretation.     

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

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

Gully formation on Mars: Two recent phases of formation suggested by links between morphology, slope orientation and insolation history

The unusual 80 km diameter Noachian-aged Asimov crater in Noachis Terra (46°S, 5°E) is characterized by extensive Noachian-Hesperian crater fill and a younger superposed annulus of valleys encircling the margins of the crater floor. These valleys provide an opportunity to study the relationships of gully geomorphology as a function of changing slope orientation relative to solar insolation. We found that the level of development of gullies was highly correlated with slope orientation and solar insolation. The largest and most complex gully systems, with the most well-developed fluvial landforms, are restricted to pole-facing slopes. In contrast, gullies on equator-facing slopes are smaller, more poorly developed and integrated, more highly degraded, and contain more impact craters. We used a 1D version of the Laboratoire de Météorologie Dynamique GCM, and slope geometries (orientation and angle), driven by predicted spin-axis/orbital parameter history, to assess the distribution and history of surface temperatures in these valleys during recent geological history. Surface temperatures on pole-facing slopes preferential for water ice accumulation and subsequent melting are predicted to occur as recently as 0.5-2.1 Ma, which is consistent with age estimates of gully activity elsewhere on Mars. In contrast, the 1D model predicts that water ice cannot accumulate on equator-facing slopes until obliquities exceed 45°, suggesting they are unlikely to have been active over the last 5 Ma. The correlation of the temperature predictions and the geological evidence for age differences suggests that there were two phases of gully formation in the last few million years: an older phase in which top-down melting occurred on equator-facing slopes and a younger more robust phase on pole-facing slopes. The similarities of small-scale fluvial erosion features seen in the gullies on Mars and those observed in gullies cut by seasonal and perennial snowmelt in the Antarctic Dry Valleys supports a top-down melting origin for these gullies on Mars.     

Titan’s global crater population: A new assessment

We report a revised crater population for Titan using Cassini RADAR data through January 2010 (flyby T65), and make a size-dependent correction for the incomplete coverage (～33%) using a Monte-Carlo model. Qualitatively, Titan’s landscape is more heavily cratered than Earth, but much less than Mars or Ganymede: the area fraction covered by craters is in fact comparable with that of Venus. Quantitative efforts to interpret crater densities for Titan as surface age have been confounded by widely divergent crater production rates proposed in the literature. We elucidate the specific model assumptions that lead to these differences (assumed projectile density, scaling function for simple crater diameter, and complex crater size exponent) and suggest these are reasonable bounding models, with the Korycansky and Zahnle (2005) model representing a crater retention age of ～1 Ga, and the Artemieva and Lunine (2005) model representing a crater retention age of ～200 Ma. These estimates are consistent with models of Titan’s evolution that predict a thickening of its crust 0.3–1.2 Gyr ago.     

Degradation of mid-latitude craters on Mars

Recent geomorphic, remote sensing, and atmospheric modeling studies have shown evidence for abundant ground ice deposits in the martian mid-latitudes. Numerous potential water/ice-rich flow features have been identified in craters in these regions, including arcuate ridges, gullies, and small flow lobes. Previous studies (such as in Newton Basin) have shown that arcuate ridges and gullies are mainly found in small craters (∼2-30 km in diameter). These features are located on both pole-facing and equator-facing crater walls, and their orientations have been found to be dependent on latitude. We have conducted surveys of craters >20 km in diameter in two mid-latitude regions, one in the northern hemisphere in Arabia Terra, and one in the southern hemisphere east of Hellas basin. In these regions, prominent lobes, potentially ice-rich, are commonly found on the walls of craters with diameters between ∼20-100 km. Additional water/ice-rich features such as channels, valleys, alcoves, and debris aprons have also been found in association with crater walls. In the eastern Hellas study region, channels were found to be located primarily on pole-facing walls, whereas valleys and alcoves were found primarily on equator-facing walls. In the Arabia Terra study region, these preferences are less distinct. In both study regions, lobate flows, gullies, and arcuate ridges were found to have pole-facing orientation preferences at latitudes below 45° and equator-facing orientation preferences above 45°, similar to preferences previously found for gullies and arcuate ridges in smaller craters. Interrelations between the features suggest they all formed from the mobilization of accumulated ice-rich materials. The dependencies of orientations on latitude suggest a relationship to differences in total solar insolation along the crater walls. Differences in slope of the crater wall, differences in total solar insolation with respect to wall orientation, and variations in topography along the crater rim can explain the variability in morphology of the features studied. The formation and evolution of these landforms may best be explained by multiple cycles of deposition of ice-rich material during periods of high obliquity and subsequent modification and transport of these materials down crater walls.     

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.