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.     

Martian cratering 12. Utilizing primary crater clusters to study crater populations and meteoroid properties

Images from Mars Global Surveyor and later images from Mars Reconnaissance Orbiter reveal that roughly half of the meteoroids striking Mars (at meter to few decameter crater diameters) fragment in the Martian atmosphere, producing small clusters of primary impact craters. Statistics of these “primary clusters” yield valuable information about important Martian phenomena and properties of interplanetary bodies, including meteoroid behavior in the Martian atmosphere, bulk strengths of bodies striking Mars, and the fraction of Martian “field secondary” craters, a datum that would improve crater count chronometry. Many Martian impactors fragment at altitudes significantly higher than 18 km above the mean surface of Mars, and we find that most bodies striking Mars and Earth have low bulk strengths, consistent with crumbly or highly fractured objects. Applying statistics of primary clusters at various elevations and independent diameter bins, we describe a technique to estimate the percentage of semirandomly scattered “field secondary” craters. Our provisional estimate of this percentage, in the diameter range ~250 m down to ~22 m, is ~40% to ~80% of the total impacts, with the higher percentages at smaller diameters. Our data argue against earlier suggestions of overwhelming dominance by either primaries or secondaries in this diameter range.     

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.     

Constraints on the Martian cratering rate based on the SNC meteorites and implications for Mars climatic history

Two constraints placed upon the cratering flux at Mars by the SNC meteorites are examined: crystallization ages as a constraint on surface ages and cosmic ray exposure ages and number of impacts as a constraint on absolute rates. The crystallization ages of the SNC meteorites appear to constrain the Martian cratering rate to be 4xLunar or more if the parent lavas are in the north of Mars and the number of SNC ejecting impacts are small. If the SNCs result from a single impact that formed the Lyot basin then the cratering rate must be at least 7xLunar or higher to produce a basin age less than the SNC crystallization age because the basin ages are themselves determined by crater counting. Assuming multiple uncorrelated impacts for SNC ejection from Mars over 10 million years a cratering rate of approximately 4xLunar is also found for ejecting impacts that form craters over 12km in diameter. Therefore, both crystallization ages and ejection ages and number of impacts appear consistent with a 4xLunar cratering rate at Mars. The effect on Martian chronologies of such a high cratering rate is to place the SNC crystallization ages partly within the epoch of channel formation on Mars and to extend this liquid water epoch over much of Mars history.     

K‐Ar dating of rocks on Mars: Requirements from Martian meteorite analyses and isochron modeling

Abstract— Radiometric age dating of Martian rocks and surfaces at known locations for which crater densities can be determined is highly desirable in order to fully understand Martian history. Performing K‐Ar age dating of igneous rocks on Mars by robots, however, presents technical challenges. Some of these challenges can be defined by examining Ar‐Ar data acquired on Martian meteorites, and others can be evaluated through numerical modeling of simulated K‐Ar isochrons like those that would be acquired robotically on Martian rocks. Excess ⁴⁰Ar is present in all shergottites. Thus for Martian rocks, the slopes of K‐Ar isochrons must be determined to reasonable precision in order to calculate reliable ages. Model simulations of possible isochrons give an indication of some requirements in order to define a precise rock age: Issues addressed here are: how many K‐Ar analyses should be made of rocks thought to have the same age; what range of K concentrations should these analyzed samples have; and what analytical uncertainty in K‐Ar measurements is desirable. Meteorite data also are used to determine the D/a² diffusion parameters for Ar in plagioclase and pyroxene separates of several shergottites and nakhlites. These data indicate the required temperatures and times for heating similar Martian rocks in order to extract Ar. Quantitatively extracting radiogenic ⁴⁰Ar could be difficult, and degassing cosmogenic Ar from mafic phases even more so. Considering all these factors, robotic K‐Ar dating of Martian rocks may be achievable, but will be challenging.     

Crater clusters on Mars: Shedding light on martian ejecta launch conditions

We have identified two classes of crater clusters on Mars. One class is “small clusters” (crater diameter D∼ tens m, spread over few hundred m), fitting our earlier calculations for the breakup of weak stone meteoroids in the martian atmosphere [Popova, O.P., Nemtchinov, I.V., Hartmann, W.K., 2003. Meteorit. Planet. Sci. 38, 905-925]. The second class is “large clusters” (D∼ few hundred m, spread over 2 to 30 km), which do not fit any predictions for breakup of known meteoroid types. We consider a range of possible explanations. The best explanation relates to known, high-speed ejection of large, semi-coherent, fractured rock masses from the surface, as secondary debris from primary impacts. The clusters are probably due to breakup of partly fracture, few-hundred-meter scale weak blocks, especially during ascent (producing moderate lateral spreading velocities among the fragments during sub-orbital flight), and also during descent of the resulting swarm. These conclusions illuminate the launch conditions of martian meteorites, including fragmentation processes, although more work is needed on the lateral separation of fragments (during either atmosphere descent or ascent) due to the effects of volatiles in the projectiles. Martian meteorites probably come from smaller martian craters than the clusters' source craters. The latter probably have D?85 km, although we have not ruled out diameters as small as 15 km.     

Comparison of the crater morphology-size relationship for Mars,Moon, and Mercury

New central peak-crater size data for Mars shows that a higher percentage of relatively unmodified Martian craters have central peaks than do fresh lunar craters below a diameter of 30 km. For example, in the diameter range 10 to 20 km, 60% of studied Martian craters have central peaks compared to 26% for the Moon. Gault et al. (1975, J. Geophys. Res.80, 2444–2460) have demonstrated that central peaks occur in smaller craters on Mercury than on the Moon, and that this effect is due to the different gravity fields in which the craters formed. Similar differences when comparing Mars and the Moon show that gravity has affected the diameter at which central peaks form on Mars. Erosion on Mars, therefore, does not completely mask differences in crater interior structure that are caused by differences in gravity. Effects of Mars' higher surface gravity when compared to the Moon are not detected when comparing terrace and crater shape data. The morphology-crater size statistics also show that a full range of crater shapes occur on Mars, and craters tend to become more morphologically complex with increasing diameter. Comparisons of Martian and Mercurian crater data show differences which may be related to the greater efficacy of erosion on Mars.     

Comparison of central pit craters on Mars,Mercury, Ganymede,and the Saturnian satellites

We report on the first results of a large‐scale comparison study of central pit craters throughout the solar system, focused on Mars, Mercury, Ganymede, Rhea, Dione, and Tethys. We have identified 10 more central pit craters on Rhea, Dione, and Tethys than have previously been reported. We see a general trend that the median ratio of the pit to crater diameter (D_p/D_c) decreases with increasing gravity and decreasing volatile content of the crust. Floor pits are more common on volatile‐rich bodies while summit pits become more common as crustal volatile content decreases. Uplifted bedrock from below the crater floor occurs in the central peak upon which summit pits are found and in rims around floor pits, which may or may not break the surface. Peaks on which summit pits are found on Mars and Mercury share similar characteristics to those of nonpitted central peaks, indicating that some normal central peaks undergo an additional process to create summit pits. Martian floor pits do not appear to be the result of a central peak collapse as the median ratio of the peak to crater diameter (D_pk/D_c) is about twice as high for central peaks/summit pits than D_p/D_c values for floor pits. Median D_pk/D_c is twice as high for Mars as for Mercury, reflecting differing crustal strength between the two bodies. Results indicate that a complicated interplay of crustal volatiles, target strength, surface gravity, and impactor energy along with both uplift and collapse are involved in central pit formation. Multiple formation models may be required to explain the range of central pits seen throughout the solar system.     

Martian cratering revisited: Implications for early geologic evolution

Improved crater statistics from varied Martian terrains are compared to lunar crater populations. The distribution functions for the average Martian cratered terrain and the average lunar highlands over the diameter range 8–2000 km are quite similar. The Martian population is less dense by approximately 0.70 from 8 to 256 km diameter and diverges to proportionally lower densities at greater diameters. Crater densities on Martian “pure” terra give a lower limit to the Mars/Moon integrated crater flux of 0.75 since the last stabilization of the respective planetary crusts. The crater population >8 km diameter postdating the Martian northern plains is statistically indistinguishable from that population postdating the lunar maria. Monte Carlo simulations were performed to constrain plausible mechanisms of crater obliteration. The models demonstrate that if the crater density difference between the lunar and Martian terra has been due to resurfacing processes, random intercrater plains formation cannot be the sole process. If plains preferentially form in and obliterate larger craters, then the observed Martian distribution retains its “shape” as the crater density decreases. This result is consistent with the morphology of Martian intercrater plains.     

Crater geometry and ejecta thickness of the Martian impact crater Tooting

Abstract— We use Mars Orbiter Laser Altimeter (MOLA) topographic data and Thermal Emission Imaging System (THEMIS) visible (VIS) images to study the cavity and the ejecta blanket of a very fresh Martian impact crater ?29 km in diameter, with the provisional International Astronomical Union (IAU) name Tooting crater. This crater is very young, as demonstrated by the large depth/diameter ratio (0.065), impact melt preserved on the walls and floor, an extensive secondary crater field, and only 13 superposed impact craters (all 54 to 234 meters in diameter) on the ?8120 km² ejecta blanket. Because the pre‐impact terrain was essentially flat, we can measure the volume of the crater cavity and ejecta deposits. Tooting crater has a rim height that has >500 m variation around the rim crest and a very large central peak (1052 m high and >9 km wide). Crater cavity volume (i.e., volume below the pre‐impact terrain) is ?380 km³ the volume of materials above the pre‐impact terrain is ?425 km³. The ejecta thickness is often very thin (<20 m) throughout much of the ejecta blanket. There is a pronounced asymmetry in the ejecta blanket, suggestive of an oblique impact, which has resulted in up to ?100 m of additional ejecta thickness being deposited down‐range compared to the up‐range value at the same radial distance from the rim crest. Distal ramparts are 60 to 125 m high, comparable to the heights of ramparts measured at other multi‐layered ejecta craters. Tooting crater serves as a fresh end‐member for the large impact craters on Mars formed in volcanic materials, and as such may be useful for comparison to fresh craters in other target materials.     

Martian subsurface properties and crater formation processes inferred from fresh impact crater geometries

Abstract— The geometry of simple impact craters reflects the properties of the target materials, and the diverse range of fluidized morphologies observed in Martian ejecta blankets are controlled by the near‐surface composition and the climate at the time of impact. Using the Mars Orbiter Laser Altimeter (MOLA) data set, quantitative information about the strength of the upper crust and the dynamics of Martian ejecta blankets may be derived from crater geometry measurements. Here, we present the results from geometrical measurements of fresh craters 3–50 km in rim diameter in selected highland (Lunae and Solis Plana) and lowland (Acidalia, Isidis, and Utopia Planitiae) terrains. We find large, resolved differences between the geometrical properties of the freshest highland and lowland craters. Simple lowland craters are 1.5–2.0 times deeper (≥5s?_o difference) with >50% larger cavities (≥2s?_o) compared to highland craters of the same diameter. Rim heights and the volume of material above the preimpact surface are slightly greater in the lowlands over most of the size range studied. The different shapes of simple highland and lowland craters indicate that the upper ?6.5 km of the lowland study regions are significantly stronger than the upper crust of the highland plateaus. Lowland craters collapse to final volumes of 45–70% of their transient cavity volumes, while highland craters preserve only 25–50%. The effective yield strength of the upper crust in the lowland regions falls in the range of competent rock, approximately 9–12 MPa, and the highland plateaus may be weaker by a factor of 2 or more, consistent with heavily fractured Noachian layered deposits. The measured volumes of continuous ejecta blankets and uplifted surface materials exceed the predictions from standard crater scaling relationships and Maxwell's Z model of crater excavation by a factor of 3. The excess volume of fluidized ejecta blankets on Mars cannot be explained by concentration of ejecta through nonballistic emplacement processes and/or bulking. The observations require a modification of the scaling laws and are well fit using a scaling factor of ?1.4 between the transient crater surface diameter to the final crater rim diameter and excavation flow originating from one projectile diameter depth with Z = 2.7. The refined excavation model provides the first observationally constrained set of initial parameters for study of the formation of fluidized ejecta blankets on Mars.     

Determination of Mars crater geometric data: Insights from high-resolution digital elevation models

We review the methods and data sets used to determine morphometric parameters related to the depth (e.g., rim height and cavity depth) and diameter of Martian craters over the past ~45 yr, and discuss the limitations of shadow length measurements, photoclinometry, Earth-based radar, and laser altimetry. We demonstrate that substantial errors are introduced into crater depth and diameter measurements that are inherent in the use of 128th-degree gridded Mars Orbiter Laser Altimeter (MOLA) topography. We also show that even the use of the raw MOLA Precision Engineering Data Record (PEDR) data can introduce errors in the measurement of craters a few kilometers in diameter. These errors are related to the longitudinal spacing of the MOLA profiles, the along-track spacing of the individual laser shots, and the MOLA spot size. Stereophotogrammetry provides an intrinsically more accurate method for measuring depth and diameter of craters on Mars when applied to high-resolution image pairs. Here, we use 20 stereo Context Camera (CTX) image pairs to create digital elevation models (DEMs) for 25 craters in the diameter range 1.5–25.6 km and cover the latitude range of 25^° S to 42^° N. These DEMs have a spatial scale of ~24 m per pixel. Six additional craters, 1.5–3.1 km in diameter, were studied using publically available DEMs produced from High-Resolution Imaging Science Experiment (HiRISE) image pairs. Depth/diameter and rim height were determined for each crater, as well as the azimuthal variation of crater rim height in 1-degree increments. These data indicate that morphologically fresh Martian craters at these diameters are significantly deeper for a given size than previously reported using Viking and MOLA data, most likely due to the improvement in spatial resolution provided by the CTX and HiRISE data.     

Distribution of small channels on the Martian surface

The distribution of small channels on Mars has been mapped from Mariner 9 images, at the 1:5 000 000 scale, by the author. The small channels referred to here are small valleys ranging in width from the resolution limit of the Mariner 9 wide-angle images (～1 km) to about 10 km. The greatest density of small band occurs in dark cratered terrain. This dark zone forms a broad subequatorial band around the planet. The observed distribution may be the result of decreased small-channel visibility in bright areas due to obscuration by a high albedo dust or sediment mantle. Crater densities within two small-channel segments show crater size-frequency distributions consistent with those of the oldest of the heavily cratered plains units. Such crater densities coupled with the almost exclusive occurrence of small channels in old cratered terrain and the generally degraded appearance of small channels in the high-resolution images (～100 m) imply a major episode of small-channel formation early in Martian geologic history.     

Relative crater production rates on planets

Dynamical histories of planetesimals in specified orbits, calculated by Wetherill (1975) and others, have estimates of relative numbers of impacts on different planets. These impact rates, $\text{F}$, are converted to crater production rates, F, by means of tables developed in this paper. Conversions are dependent on impact velocity and surface gravity. Crater retention ages can then be derived from (crater density)/(crater production rate). Such calculations of impact rates and their histories give the only basis, independent of sample dating, for establishing absolute geologic histories of the planets, contrary to published implications that this can be done by comparison of photos alone. A survey of the results, from orbits of interplanetary objects studied to date, indicates that the terrestial planets have crater production rates within a factor ten of each other, and that planet's crater retention ages can probably be determined with a factor of ±3. Further calculations of orbital histories of additional interplanetary bodies are suggested to put photogeologic analyses from spacecraft imagery on a firmer basis.Applications to Mars, as an example, using least-squares fits to crater-count data, suggest an average age of 0.3 to 3 b.y. for two types of channels. The Tharsis volcanics are found to be slightly younger than the channels (strongly confirmed by photomorphology since they are not cut by channels) and Olympus Mons is about 0.06 to 0.6 b.y. old, contrary to recent assertions that Olympus Mons is 2.5 b.y. old and most Martian volcanic provinces older than 3 b.y. Data strongly support the hypothesis that Martian channels formed in a fluvial climate that persisted on Mars until the Tharsis volcanism caused a change in the Martian obliquity state, as outlined by Toon, Ward, and Burns (1977).     

Characterization of fluvial activity in Parana Valles using different age-dating techniques

Martian valley networks provide the best evidence that the climate on Mars was different in the past. Although these features are located primarily in heavily cratered terrain of Noachian age (>3.7 Ga), the ages of the features and the time when they were active is not well understood. From superposed craters several recent global studies determined that most valley networks formed during the Late Noachian to Early Hesperian; however, there were some disparities between the techniques. In this study, our principal objective was to test the reliability of the different age-dating techniques to better understand their accuracy and limitations. We applied these techniques to Parana Valles using a variety of high-resolution images taken from different instruments that allow us to identify smaller craters (D > 125 m) while providing sufficient coverage to support a statistically reliable sampling of crater populations, which is necessary to reduce the uncertainties in age determination. Our results indicate that Parana Valles formed during the Early Hesperian Period but that the crater density (D > 353 m) is heterogeneous inside the Parana Valles basin. The crater population decreases from the headwaters downstream recording a resurfacing event that is most likely related to the erosion of downstream sub-basins. The terrain near the source area is Late Noachian to Early Hesperian in age while terrains closer to the outlet are Early to Late Hesperian in age. Crater densities (D > 125 m) inside the valley are also heterogeneous and record several resurfacing events on the valley floor. Where the width of the valley network narrows to <2 km we found evidence of an Amazonian age eolian deposit that is a relatively thin layer of only few meters that was probably deposited as a result of topographic influences. Our results validate the reliability of several proposed age-dating techniques, but we also determined the accuracy and applicability of these techniques. Our results also demonstrate that crater populations can be used to not only determine the relative ages of valley networks, but also to map the distribution of sedimentary materials and the extent of resurfacing events that occurred after valley network formation.     

Insights into the Martian mantle: The age and isotopics of the meteorite fall Tissint

The recent witnessed fall of the meteorite Tissint represents the delivery of a pristine new sample from the surface of Mars. This meteorite provides an unprecedented opportunity to study a variety of aspects about the planet's evolution. Using the Rb–Sr and Sm–Nd isotopic systems, we determined that Tissint, a depleted shergottite, has a crystallization age of 574 ± 20 Ma, an initial ε¹⁴³Nd = +42.2 ± 0.5, and an initial ⁸⁷Sr/⁸⁶Sr = 0.700760 ± 11. These initial Nd and Sr isotopic compositions suggest that Tissint originated from a mantle source on Mars that is distinct from the source reservoirs of the other Martian meteorites. The known crystallization ages, geochemical characteristics, ejection ages, and ejection dynamics of Tissint and other similarly grouped Martian meteorites suggest that they are likely derived from a source crater up to approximately 90 km in diameter with an age of approximately 1 Ma that is located on terrain that is approximately 600 million years old.     

Layered ejecta craters and the early water/ice aquifer on Mars

Abstract— A model for emplacement of deposits of impact craters is presented that explains the size range of Martian layered ejecta craters between 5 km and 60 km in diameter in the low and middle latitudes. The impact model provides estimates of the water content of crater deposits relative to volatile content in the aquifer of Mars. These estimates together with the amount of water required to initiate fluid flow in terrestrial debris flows provide an estimate of 21% by volume (7.6 × 10⁷km³) of water/ice that was stored between 0.27 and 2.5 km depth in the crust of Mars during Hesperian and Amazonian time. This would have been sufficient to supply the water for an ocean in the northern lowlands of Mars. The existence of fluidized craters smaller than 5 km diameter in some places on Mars suggests that volatiles were present locally at depths less than 0.27 km. Deposits of Martian craters may be ideal sites for searches for fossils of early organisms that may have existed in the water table if life originated on Mars.     

Highly siderophile element (HSE) abundances in the mantle of Mars are due to core formation at high pressure and temperature

Highly siderophile elements (HSEs) can be used to understand accretion and core formation in differentiated bodies, due to their strong affinity for FeNi metal and sulfides. Coupling experimental studies of metal–silicate partitioning with analyses of HSE contents of Martian meteorites can thus offer important constraints on the early history of Mars. Here, we report new metal–silicate partitioning data for the PGEs and Au and Re across a wide range of pressure and temperature space, with three series designed to complement existing experimental data sets for HSE. The first series examines temperature effects for D(HSE) in two metallic liquid compositions—C‐bearing and C‐free. The second series examines temperature effects for D(Re) in FeO‐bearing silicate melts and FeNi‐rich alloys. The third series presents the first systematic study of high pressure and temperature effects for D(Au). We then combine our data with previously published partitioning data to derive predictive expressions for metal–silicate partitioning of the HSE, which are subsequently used to calculate HSE concentrations of the Martian mantle during continuous accretion of Mars. Our results show that at midmantle depths in an early magma ocean (equivalent to approximately 14 GPa, 2100 °C), the HSE contents of the silicate fraction are similar to those observed in the Martian meteorite suite. This is in concert with previous studies on moderately siderophile elements. We then consider model calculations that examine the role of melting, fractional crystallization, and sulfide saturation/undersaturation in establishing the range of HSE contents in Martian meteorites derived from melting of the postcore formation mantle. The core formation modeling indicates that the HSE contents can be established by metal–silicate equilibrium early in the history of Mars, thus obviating the need for a late veneer for HSE, and by extension volatile siderophile elements, or volatiles in general.