A geomorphic analysis of Hale crater, Mars: The effects of impact into ice-rich crust

Hale crater, a 125 × 150 km impact crater located near the intersection of Uzboi Vallis and the northern rim of Argyre basin at 35.7°S, 323.6°E, is surrounded by channels that radiate from, incise, and transport material within Hale’s ejecta. The spatial and temporal relationship between the channels and Hale’s ejecta strongly suggests the impact event created or modified the channels and emplaced fluidized debris flow lobes over an extensive area (>200,000 km²). We estimate ∼10¹⁰ m³ of liquid water was required to form some of Hale’s smaller channels, a volume we propose was supplied by subsurface ice melted and mobilized by the Hale-forming impact. If 10% of the subsurface volume was ice, based on a conservative porosity estimate for the upper martian crust, 10¹² m³ of liquid water could have been present in the ejecta. We determine a crater-retention age of 1 Ga inside the primary cavity, providing a minimum age for Hale and a time at which we propose the subsurface was volatile-rich. Hale crater demonstrates the important role impacts may play in supplying liquid water to the martian surface: they are capable of producing fluvially-modified terrains that may be analogous to some landforms of Noachian Mars.     

In situ measurements of particle load and transport in dust devils

In situ (mobile) sampling of 33 natural dust devil vortices reveals very high total suspended particle (TSP) mean values of 296 mg m⁻³ and fine dust loadings (PM10) mean values ranging from 15.1 to 43.8 mg m⁻³ (milligrams per cubic meter). Concurrent three-dimensional wind profiles show mean tangential rotation of 12.3 m s⁻¹ and vertical uplift of 2.7 m s⁻¹ driving mean vertical TSP flux of 1689 mg m⁻³ s⁻¹ and fine particle flux of ∼1.0 to ∼50 mg m⁻³ s⁻¹. Peak PM10 dust loading and flux within the dust column are three times greater than mean values, suggesting previous estimates of dust devil flux might be too high. We find that deflation rates caused by dust devil erosion are ∼2.5-50 μm per year in dust devil active zones on Earth. Similar values are expected for Mars, and may be more significant there where competing erosional mechanisms are less likely.     

A lake in Uzboi Vallis and implications for Late Noachian-Early Hesperian climate on Mars

Uzboi Vallis (centered at ∼28°S, 323°E) is ∼400 km long and comprises the southernmost segment of the northward-draining Uzboi-Ladon-Morava (ULM) meso-scale outflow system that emerges from Argyre basin. Bond and Holden craters blocked the valley to the south and north, respectively, forming a Late Noachian-to-Hesperian paleolake basin that exceeded 4000 km³. Limited CRISM data suggest lake deposits in Uzboi and underlying basin floor incorporate relatively more Mg-clays and more Fe-clays, respectively. The short-lived lake overflowed and breached Holden crater’s rim at an elevation of −350 m and rapidly drained into the crater. Fan deltas in Holden extend 25 km from the breach and incorporate meter-sized blocks, and longitudinal grooves along the Uzboi basin floor are hundreds of meters long and average 60 m wide, suggesting high-discharge drainage of the lake. Precipitation-derived runoff rather than regional groundwater or overflow from Argyre dominated contributions to the Uzboi lake, although the failure of most tributaries to respond to a lowering of base level indicates their incision largely ended when the lake drained. The Uzboi lake may have coincided with alluvial and/or lacustrine activity in Holden, Eberswalde, and other craters in southern Margaritifer Terra, where fluvial/lacustrine activity may have required widespread, synoptic precipitation (rain or snow), perhaps associated with an ephemeral, global hydrologic system during the Late Noachian into the Hesperian on Mars.     

Dust devil sediment flux on Earth and Mars: Laboratory simulations

Laboratory simulations using the Arizona State University Vortex Generator (ASUVG) were run to simulate sediment flux in dust devils in terrestrial ambient and Mars-analog conditions. The objective of this study was to measure vortex sediment flux in the laboratory to yield estimations of natural dust devils on Earth and Mars, where all parameters may not be measured. These tests used particles ranging from 2 to 2000 μm in diameter and 1300 to 4800 kg m⁻³ in density, and the results were compared with data from natural dust devils on Earth and Mars. Typically, the cores of dust devils (regardless of planetary environment) have a pressure decrease of ∼0.1-1.5% of ambient atmospheric pressure, which enhances the lifting of particles from the surface. Core pressure decreases in our experiments ranged from ∼0.01% to 5.00% of ambient pressure (10 mbar Mars cases and 1000 mbar for Earth cases) corresponding to a few tenths of a millibar for Mars cases and a few millibars for Earth cases. Sediment flux experiments were run at vortex tangential wind velocities of 1-45 m s⁻¹, which typically correspond to ∼30-70% above vortex threshold values for the test particle sizes and densities. Sediment flux was determined by time-averaged measurements of mass loss for a given vortex size. Sediment fluxes of ∼10⁻⁶-10⁰ kg m⁻² s⁻¹ were obtained, similar to estimates and measurements for fluxes in dust devils on Earth and Mars. Sediment flux is closely related to the vortex intensity, which depends on the strength of the pressure decrease in the core (ΔP). This study found vortex size is less important for lifting materials because many different diameters can have the same ΔP. This finding is critical in scaling the laboratory results to natural dust devils that can be several orders of magnitude larger than the laboratory counterparts.     

Scale-dependent infrared radiative damping rates on Mars and their role in the deposition of gravity-wave momentum flux

Using a Curtis-matrix model of 15 μm CO₂ radiative cooling rates for the martian atmosphere, we have computed vertical scale-dependent IR radiative damping rates from 0 to 200 km altitude over a broad band of vertical wavenumbers ∣m∣ = 2π(1-500 km)⁻¹ for representative meteorological conditions at 40°N and average levels of solar activity and dust loading. In the middle atmosphere, infrared (IR) radiative damping rates increase with decreasing vertical scale and peak in excess of 30 days⁻¹ at ∼50-80 km altitude, before gradually transitioning to scale-independent rates above ∼100 km due to breakdown of local thermodynamic equilibrium. We incorporate these computed IR radiative damping rates into a linear anelastic gravity-wave model to assess the impact of IR radiative damping, relative to wave breaking and molecular viscosity, in the dissipation of gravity-wave momentum flux. The model results indicate that IR radiative damping is the dominant process in dissipating gravity-wave momentum fluxes at ∼0-50 km altitude, and is the dominant process at all altitudes for gravity waves with vertical wavelengths ?10-15 km. Wave breaking becomes dominant at higher altitudes only for “fast” waves of short horizontal and long vertical wavelengths. Molecular viscosity plays a negligible role in overall momentum flux deposition. Our results provide compelling evidence that IR radiative damping is a major, and often dominant physical process controlling the dissipation of gravity-wave momentum fluxes on Mars, and therefore should be incorporated into future parameterizations of gravity-wave drag within Mars GCMs. Lookup tables for doing so, based on the current computations, are provided.     

The SMART-1 lunar impact

The SMART-1 spacecraft impacted the Moon on 3rd September 2006 at a speed of 2 km s⁻¹ and at a very shallow angle of incidence (∼1°). The resulting impact crater is too small to be viewed from the Earth; accordingly, the general crater size and shape have been determined here by laboratory impact experiments at the same speed and angle of incidence combined with extrapolating to the correct size scale to match the SMART-1 impact. This predicts a highly asymmetric crater approximately 5.5-26 m long, 1.9-9 m wide, 0.23-1.5 m deep and 0.71-6.9 m³ volume. Some of the excavated mass will have gone into crater rim walls, but 0.64-6.3 m³ would have been ejecta on ballistic trajectories corresponding to a cloud of 2200-21,800 kg of lunar material moving away from the impact site. The shallow Messier crater on the Moon is similarly asymmetric and is usually taken as arising from a highly oblique impact. The light flash from the impact and the associated ejecta plume were observed from Earth, but the flash magnitude was not obtained, so it is not possible to obtain the luminous efficiency of the impact event.     

The rate of granule ripple movement on Earth and Mars

The rate of movement for 3- and 10-cm-high granule ripples was documented in September of 2006 at Great Sand Dunes National Park and Preserve during a particularly strong wind event. Impact creep induced by saltating sand caused ∼24 granules min⁻¹ to cross each cm of crest length during wind that averaged ∼9 m s⁻¹ (at a height well above 1 m), which is substantially larger than the threshold for saltation of sand. Extension of this documented granule movement rate to Mars suggests that a 25-cm-high granule ripple should require from hundreds to thousands of Earth-years to move 1 cm under present atmospheric conditions.     

Generation and emplacement of fine-grained ejecta in planetary impacts

We report here on a survey of distal fine-grained ejecta deposits on the Moon, Mars, and Venus. On all three planets, fine-grained ejecta form circular haloes that extend beyond the continuous ejecta and other types of distal deposits such as run-out lobes or ramparts. Using Earth-based radar images, we find that lunar fine-grained ejecta haloes represent meters-thick deposits with abrupt margins, and are depleted in rocks ?1 cm in diameter. Martian haloes show low nighttime thermal IR temperatures and thermal inertia, indicating the presence of fine particles estimated to range from ∼10 μm to 10 mm. Using the large sample sizes afforded by global datasets for Venus and Mars, and a complete nearside radar map for the Moon, we establish statistically robust scaling relationships between crater radius R and fine-grained ejecta run-out r^* for all three planets. On the Moon, r^* ∼ R^−0.18 for craters 5-640 km in diameter. For Venus, radar-dark haloes are larger than those on the Moon, but scale as r^* ∼ R^−0.49, consistent with ejecta entrainment in Venus’ dense atmosphere. On Mars, fine-ejecta haloes are larger than lunar haloes for a given crater size, indicating entrainment of ejecta by the atmosphere or vaporized subsurface volatiles, but scale as R^−0.13, similar to the ballistic lunar scaling. Ejecta suspension in vortices generated by passage of the ejecta curtain is predicted to result in ejecta run-out that scales with crater size as R^1/2, and the wind speeds so generated may be insufficient to transport particles at the larger end of the calculated range. The observed scaling and morphology of the low-temperature haloes leads us rather to favor winds generated by early-stage vapor plume expansion as the emplacement mechanism for low-temperature halo materials.     

Impact megadomes and the origin of the martian crustal dichotomy

We show that a sufficiently energetic impact can generate a melt volume which, after isostatic adjustment and differentiation, forms a spherical cap of crust with underlying depleted mantle. Depending on impact energy and initial crustal thickness, a basin may be retained or impact induced crust may be topographically elevated. Retention of a martian lowland scale impact basin at impact energies ∼3 × 10²⁸-3 × 10²⁹ J requires an initial crustal thickness greater than 10 km. Formation of impact induced crust with size comparable to the martian highlands requires a larger impact energy, ∼1-3 × 10³⁰ J, and initial crustal thickness <20 km. Furthermore, we show that the boundary of impact induced crust can be elliptical due to a spatially asymmetric impact melt volume caused by an oblique impact. We suggest the term “impact megadome” for topographically elevated, impact induced crust and propose that processes involved in megadome formation may play an important role in the origin of the martian crustal dichotomy.     

Possible impact melt and debris flows at Tooting Crater, Mars

Candidate examples of impact melt flows and debris flows have been identified at Tooting crater, an extremely young (<2 Myr), 29 km diameter impact crater in Amazonis Planitia, Mars. Using HiRISE and CTX images, and stereo-derived digital elevation models derived from these images, we have studied the rim and interior wall of Tooting crater to document the morphology and topography of several flow features in order to constrain the potential flow formation mechanisms. Four flow types have been identified; including possible impact melt sheets and three types of debris flows. The flow features are all located within 2 km of the rim crest on the southern rim or lie on the southern interior wall of the crater ∼1500 m below the rim crest. Extensive structural failure has modified the northern half of the crater inner wall and we interpret this to have resulted in the destruction of any impact melt emplaced, as well as volatile-rich wall rock. The impact melt flows are fractured on the meter to decameter scale, have ridged, leveed lobes and flow fronts, and cover an area >6 km × 5 km on the southern rim. The debris flows are found on both the inner wall and rim of the crater, are ∼1-2 km in length, and vary from a few tens of meters to >300 m in width. These flows exhibit varying morphologies, from a channelized, leveed flow with arcuate ridges in the channel, to a rubbly flow with a central channel but no obvious levees. The flows indicate that water existed within the target rocks at the time of crater formation, and that both melt and fluidized sediment was generated during this event.     

Evaporation effects on the formation of martian gullies

In order to investigate the formation of martian gullies and the stability of fluids on Mars, we examined about 120 gully images. Twelve HiRISE images contained a sufficient number of Transverse Aeolian Ridges (TARs) associated with the gullies to make the following measurements: overall gully length, length of the alcove, channel and apron, and we also measured the frequency of nearby TARs. Six of the 12 images examined showed a statistically significant negative correlation between overall gully length (alcove, channel and apron length) and TAR frequency. Previous experimental work from our group has shown that at temperatures below ∼200 K, evaporation rate increases by about an order of magnitude as wind speed increases from 0 to ∼15 m/s. Thus the negative correlations we observe between gully length and dune frequency can be explained by formation at temperatures below ∼200 K where wind speed/evaporation is a factor governing gully length. In these cases evaporation of the fluid carving the gully was a constraint on their dimensions. Cases where there is no correlation between gully length and TAR frequency, can be explained by formation at temperatures >200 K. The temperatures are consistent with Global Circulation Model and Thermal Emission Spectrometer (TES) data for these latitudes. The temperatures suggested by these trends are consistent with the fluid responsible for gully formation being a strong brine, such as Fe₂(SO₄)₃ which has a eutectic temperature of ∼200 K. We also find that formation timescales for gullies are 10⁵-10⁶ years.     

Adsorptive fractionation of HDO on JSC MARS-1 during sublimation with implications for the regolith of Mars

A chamber was constructed to simulate the boundary between the ice table, regolith and atmosphere of Mars and to examine fractionation between H₂O and HDO during sublimation under realistic martian conditions of temperature and pressure. Thirteen experimental runs were conducted with regolith overlying the ice. The thickness and characteristic grain size of the regolith layer as well as the temperature of the underlying ice was varied. From these runs, values for the effective diffusivity, taking into account the effects of adsorption, of the regolith were derived. These effective diffusivities ranged from 1.8 × 10⁻⁴ m² s⁻¹ to 2.2 × 10⁻³ m² s⁻¹ for bare ice and from 2.4 × 10⁻¹¹ m² s⁻¹ to 2.0 × 10⁻⁹ m² s⁻¹ with an adsorptive layer present. From these, latent heats of adsorption of 8.6 ± 2.6 kJ mol⁻¹ and 9.3 ± 2.8 kJ mol⁻¹ were derived at ice-surface temperatures above 223 ± 8 K and 96 ± 28 kJ mol⁻¹ and 104 ± 31 kJ mol⁻¹ respectively for H₂O and HDO were derived at colder temperatures. For temperatures below 223 K, the effective diffusivity of HDO was found to be lower than the diffusivity of H₂O by 40% on average, suggesting that the regolith was adsorptively fractionating the sublimating gas with a fractionation factor of 1.96 ± 0.74. Applying these values to Mars predicts that adsorbed water on the regolith is enriched in HDO compared to the atmosphere, particularly where the regolith is colder. Based on current observations, the D/H ratio of the regolith may be as high as 21 ± 8 times VSMOW at 12°S and L_S = 357° if the regolith is hydrated primarily by the atmosphere, neglecting any hydration from subsurface ice.     

Groundwater discharge and gully formation on martian slopes

Young gullies and gully deposits on walls of martian craters have been cited as evidence that liquid water flowed on the surface of Mars relatively recently. Effects of variable environmental conditions at the surface of Mars are modeled and applied to the case of groundwater emergence from shallow aquifers to investigate whether groundwater is a viable source to enable the erosion of these gullies. The model includes detailed treatment of ice growth in the aquifer. Model results indicate that groundwater discharge can be maintained under the current environmental conditions if the aquifer permeability is like that of terrestrial gravel or higher, if the aquifer is 350 K or warmer, or if the aquifer is a brine with a freezing point depressed to 250 K or below. Groundwater discharge cannot be maintained for the conservative case of a cold, pure water, semi-pervious aquifer. Cold (275 K) pure water pervious (gravel) aquifers, warm (350 K) pure water semi-pervious aquifers, and cold (275 K) CaCl₂ brine semi-pervious aquifers all exhibit a dependence of discharge on season, latitude and slope orientation in our modeling. Seasonal, latitudinal and azimuthal discharge variations are strongest for cold CaCl₂ brine semi-pervious aquifers, with discharges from this aquifer type favoring equator-facing slopes at mid and high southern latitudes. At all latitudes and slope azimuths under our nominal conditions, the cold pure water pervious aquifer, the cold pure water semi-pervious aquifer and the cold CaCl₂ brine semi-pervious aquifer all freeze completely shortly after the simulations are started. Discharge restarts in the summer for the cold pure water pervious aquifer and the cold brine aquifer, but discharge does not restart for the cold pure water semi-pervious aquifer. The warm pure water semi-pervious aquifer maintains daily seeps throughout the year at all but high latitudes. In the case of the cold pure water pervious aquifer, approximately 500,000 m³ of water could be discharged from a mid-latitude, 150-m thick aquifer with a 20-m wide seepage face orientated towards the equator or the pole after a single undermining-induced event before ice growth seals the seepage face. For a brine semi-pervious aquifer with the same dimensions, 200-300 m³ of water could be released from a mid-latitude 20-m wide equator-facing seepage face before the fresh exposure is sealed for the fall and winter seasons. Our results do not rule out groundwater emergence as a means of creating some recent gullies, but they indicate that rather special and perhaps unusual conditions would be required.     

HiRISE views enigmatic deposits in the Sirenum Fossae region of Mars

HiRISE images together with other recent orbital data from Mars define new characteristics of enigmatic Hesperian-aged deposits in Sirenum Fossae that are mostly 100-200 m thick, drape kilometers of relief, and often display generally low relief surfaces. New characteristics of the deposits, previously mapped as the “Electris deposits,” include local detection of meter-scale beds that show truncating relationships, a generally light-toned nature, and a variably blocky, weakly indurated appearance. Boulders shed by erosion of the deposits are readily broken down and contribute little to talus. Thermal inertia values for the deposits are ∼200 J m⁻² K⁻¹ s^−1/2 and they may incorporate hydrated minerals derived from weathering of basalt. The deposits do not contain anomalous amounts of water or water ice. Deflation may dominate degradation of the deposits over time and points to an inventory of fine-grained sediment. Together with constraints imposed by the regional setting on formation processes, these newly resolved characteristics are most consistent with an eolian origin as a loess-like deposit comprised of redistributed and somewhat altered volcanic ash. Constituent sediments may be derived from airfall ash deposits in the Tharsis region. An origin directly related to airfall ash or similar volcanic materials is less probable and emplacement by alluvial/fluvial, impact, lacustrine, or relict polar processes is even less likely.     

Minimum estimates of the amount and timing of gases released into the martian atmosphere from volcanic eruptions

Volcanism has been a major process during most of the geologic history of Mars. Based on data collected from terrestrial basaltic eruptions, we assume that the volatile content of martian lavas was typically ∼0.5 wt.% water, ∼0.7 wt.% carbon dioxide, ∼0.14 wt.% sulfur dioxide, and contained several other important volatile constituents. From the geologic record of volcanism on Mars we find that during the late Noachian and through the Amazonian volcanic degassing contributed ∼0.8 bar to the martian atmosphere. Because most of the outgassing consisted of greenhouse gases (i.e., CO₂ and SO₂) warmer surface temperatures resulting from volcanic eruptions may have been possible. Our estimates suggest that ∼1.1 × 10²¹ g (∼8 ± 1 m m⁻²) of juvenile water were released by volcanism; slightly more than half the amount contained in the north polar cap and atmosphere. Estimates for released CO₂ (1.6 × 10²¹ g) suggests that a large reservoir of carbon dioxide is adsorbed in the martian regolith or alternatively ∼300 cm cm⁻² of carbonates may have formed, although these materials would not occur readily in the presence of excess SO₂. Up to ∼120 cm cm⁻² (2.2 × 10²⁰ g) of acid rain (H₂SO₄) may have precipitated onto the martian surface as the result of SO₂ degassing. The hydrogen flux resulting from volcanic outgassing may help explain the martian atmospheric D/H ratio. The amount of outgassed nitrogen (∼1.3 mbar) may also be capable of explaining the martian atmospheric ¹⁵N/¹⁴N ratio. Minor gas constituents (HF, HCl, and H₂S) could have formed hydroxyl salts on the surface resulting in the physical weathering of geologic materials. The amount of hydrogen fluoride emitted (1.82 × 10¹⁸ g) could be capable of dissolving a global layer of quartz sand ∼5 mm thick, possibly explaining why this mineral has not been positively identified in spectral observations. The estimates of volcanic outgassing presented here will be useful in understanding how the martian atmosphere evolved over time.     

The volcanic history of central Elysium Planitia: Implications for martian magmatism

Central Elysium Planitia (CEP) is located south of Elysium Mons. Back to the era of the Viking orbiters, clues accumulated in favor of recent volcanism in relation with ground water release and the formation of long sub-parallel fissures. Four aqueous flood channel systems emanate from linear fissures. Recent eruptions of low viscosity lavas originate from these fissures and from low shield volcanoes. The objective of this paper is to constrain the volcanic history of this region, and to determine the chronological relationships with fluvial/erosional processes. New observations (e.g., new shield volcanoes and one new fluvial event) are summarized on a context map. Thirty-five surfaces have been dated from the count of about 15,000 impact craters. Ages have been cross-checked with relative stratigraphy when possible. A probabilistic approach has been introduced to compare similar ages and define periods of volcanic activity. Our results confirm that some volcanic features are extremely recent (∼2 My). Active periods are found at 2.5-3 My, 4.3 My, 13.5-16.2 My, 19 My, 21-32 My, 58 My, 71 My, 85-95 My, 134 My, 173 My and 234 My, not excluding the possibility that some of the gaps would be filled with additional crater counts. The volcanic activity thus extended for at least the last 250 My. The lava volumes have been estimated from the topographic modeling of the floor of depressions filled up by volcanic products, including the volumes of several large crater cavities buried under lavas (>20% of the total volume). Our new estimation of the total lava volume is 1.5 ± 0.2 × 10⁵ km³. This value corresponds to an average thickness of one hundred meters of lavas for the young volcanic plain. As a consequence, the total eruption rate at CEP, defined as the total volume of lava divided by the time of emplacement 1.4 × 10⁻²-1.8 × 10⁻² m³/s is lower than values typically estimated for terrestrial hot spots or large igneous provinces, suggesting longer inactive periods. The concept of mantle plumes responsible for terrestrial flood volcanism may not be applicable to the case of CEP and the mechanism proposed in Schumacher and Dreuer (2007) offers a plausible alternative to explain our observations.     

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.     

Stratigraphy, mineralogy, and origin of layered deposits inside Terby crater, Mars

The 174 km diameter Terby impact crater (28.0°S-74.1°E) located on the northern rim of the Hellas basin displays anomalous inner morphology, including a flat floor and light-toned layered deposits. An analysis of these deposits was performed using multiple datasets from Mars Global Surveyor, Mars Odyssey, Mars Express and Mars Reconnaissance Orbiter missions, with visible images for interpretation, near-infrared data for mineralogical mapping, and topography for geometry. The geometry of layered deposits was consistent with that of sediments that settled mainly in a sub-aqueous environment, during the Noachian period as determined by crater counts. To the north, the thickest sediments displayed sequences for fan deltas, as identified by 100 m to 1 km long clinoforms, as defined by horizontal beds passing to foreset beds dipping by 6-10° toward the center of the Terby crater. The identification of distinct sub-aqueous fan sequences, separated by unconformities and local wedges, showed the accumulation of sediments from prograding/onlapping depositional sequences, due to lake level and sediment supply variations. The mineralogy of several layers with hydrated minerals, including Fe/Mg phyllosilicates, supports this type of sedimentary environment. The volume of fan sediments was estimated as >5000 km³ (a large amount considering classical martian fan deltas such as Eberswalde (6 km³)) and requires sustained liquid water activity. Such a large sedimentary deposition in Terby crater is characteristic of the Noachian/Phyllosian period during which the environment favored the formation of phyllosilicates. The latter were detected by spectral data in the layered deposits of Terby crater in three distinct layer sequences. During the Hesperian period, the sediments experienced strong erosion, possibly enhanced by more acidic conditions, forming the current morphology with three mesas and closed depressions. Small fluvial valleys and alluvial fans formed subsequently, attesting to late fluvial processes dated as late Early to early Late Hesperian. After this late fluvial episode, the Terby impact crater was submitted to aeolian processes and permanent cold conditions with viscous flow features. Therefore, the Terby crater displays, in a single location, geologic features that characterize the three main periods of time on Mars, with the presence of one of the thickest sub-aqueous fan deposits reported on Mars. The filling of Terby impact crater is thus one potential “reference geologic cross-section” for Mars stratigraphy.     

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.