Meteoritics & Planetary Science 41, Nr 10, 1509–1537 (2006) Abstract available online at http://meteoritics.org Martian subsurface properties and crater formation processes inferred from fresh impact crater geometries Sarah T. STEWART* and Gregory J. VALIANT Department of Earth and Planetary Sciences, Harvard University, 20 Oxford Street, Cambridge, Massachusetts 02138, USA *Corresponding author. E-mail: sstewart@eps.harvard.edu (Received 22 October 2005; revision accepted 30 June 2006) 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 (≥5σo difference) with >50% larger cavities (≥2σ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. INTRODUCTION Impact craters are a natural probe below planetary surfaces. The final shape and size of small, simple crater cavities are determined by the strength and physical properties (e.g., composition, layering, and porosity) of the subsurface materials. The final forms of larger craters are controlled by the gravitational collapse of the shock-weakened rock (Melosh and Ivanov 1999); the transition between the strength regime and gravity regime is proportional to the effective strength of the surface (O’Keefe and Ahrens 1993). Photogrammetry studies of Martian craters determined that the transition between crater formation regimes occurs at smaller crater diameters than expected based on comparisons to Earth, Mercury, and the Moon (Pike 1980, 1988). The smaller transition diameter and shallower crater depths have been interpreted as evidence for the presence of weak, layered, or volatile-rich materials on Mars (Cintala and MouginisMark 1980; Pike 1980). Observed regional variations in a wide range of crater morphologies suggest that Martian crater forms reflect subsurface heterogeneities (Barlow and Bradley 1990). Perhaps the most intriguing aspect of Martian craters is the range of morphologies observed in their ejecta blankets. © The Meteoritical Society, 2006. Printed in USA. 1509 1510 S. T. Stewart and G. J. Valiant All fresh and many older Martian craters larger than a few kilometers (Garvin et al. 2002) are surrounded by ejecta blankets that appear to have been fluidized and that then formed a ground-hugging flow. The observed continuous ejecta blankets have a layered appearance with a distal scarp, or rampart, and flow to a radial distance about twice as far as observed on Mercury and the Moon (Mouginis-Mark 1978). While certain rampart-type ejecta have been observed on icy satellites (Horner and Greeley 1982), the diversity of fluidized ejecta morphologies is unique to Mars. The ejecta features are thought to form by entrainment of liquid water (Barlow and Bradley 1990; Carr et al. 1977; Gault and Greeley 1978; Kuz’min et al. 1988; Mouginis-Mark 1981; Stewart et al. 2001, 2004) and/or interactions between the ejecta and the atmosphere (Barnouin-Jha and Schultz 1998; Schultz 1992; Schultz and Gault 1979). Therefore, excavated materials may also provide insight into the composition of the upper crust and perhaps the history of the Martian climate. However, the physical processes that control the formation of fluidized ejecta morphologies are still poorly understood (Barlow et al. 2005). Recent crater studies have used the Mars Orbiter Laser Altimeter (MOLA) global topography data to deduce more accurate size and shape scaling relationships (Garvin and Frawley 1998; Garvin et al. 2003), to study gradational forms (Garvin et al. 2000), to investigate the sedimentary history of the northern plains (Black and Stewart, Forthcoming; Boyce et al. 2005a, 2005b; Forsberg-Taylor et al. 2004; Meresse et al. 2005), and to infer the age of buried surfaces (Frey et al. 2002b). The topography data also allow for quantitative comparisons between final ejecta geometries and ejecta flow processes (Barnouin-Jha et al. 2005). Geometrical measurements of Martian craters may also be used to study the final stages of crater collapse. Crater scaling laws relate the impact conditions to the transient crater cavity and have been validated by impact and explosion experiments and numerical simulations (Holsapple 1993; O’Keefe and Ahrens 1993; Schmidt and Housen 1987). The observed crater cavities, however, are a result of the collapse of the transient cavity, which is still an area of intense study in impact-cratering mechanics. The collapse of simple craters is dominated by crater-wall slumping (Grieve and Garvin 1984), but the collapse of complex craters is much more complicated (Melosh and Ivanov 1999). The observed differences between Mars and other rocky bodies may be illuminated by close examination of Martian crater geometries and comparison to crater scaling laws. The MOLA topography data allow, for the first time, accurate measurement of the volume of ejecta around Martian craters. Initial studies discovered that the ejecta volume around some Martian craters is larger than can be explained by the size of the excavated cavity (Garvin and Frawley 1998). Excess ejecta volume could be explained by concentration of excavated materials into the continuous ejecta blanket, bulking, and/or incorporation of surface materials into the ejecta flow. The processes controlling the formation of fluidized ejecta depend critically on the total mass and source materials for the flow. Detailed comparisons between predicted excavation volumes and observed ejecta volumes will constrain the initial conditions for Martian ejecta flow dynamics. This study focuses on deriving quantitative information about the strength of the Martian crust, including regional variability and the formation of fluidized ejecta morphologies through analyses of fresh impact crater geometries. We measure several geometrical properties of impact craters using the MOLA data set and a new crater measurement interactive toolkit (Valiant and Stewart 2004). In the “Crater Measurement Methods” section, measurement methods are described and the accuracy and precision of the toolkit are determined from measurements of simulated craters on a range of Martian background terrains and at different MOLA track resolutions. In this work, we focus on fresh craters in highland and lowland plains, where crater geometry measurements are most accurate. In the “Crater Measurements Results” section, we describe large, resolved differences in fresh crater geometries between lowland plains regions (Utopia and Isidis basins and Acidalia Planitia) and highland plateau regions (Lunae and Solis Plana). In the “Implications for Martian Subsurface Properties” section, we infer the differences in material strength in the studied highland and lowland terrains using crater size scaling relationships and discuss the origin of near-surface layers of differing strength. In the “Ejecta Volumes and the Z-EDOZ Model” section, the volumes of ejected and uplifted materials are compared to crater excavation models and we discuss the implications for the formation of fluidized ejecta blankets. The major conclusions from this work are summarized in the “Conclusions” section, and a description of variables is included in Appendix A. CRATER MEASUREMENT METHODS Measurement Toolkit Crater measurements are executed using the HMars interactive toolkit, a graphical, cross-platform program written in the Interactive Data Language (IDL), shown in Fig. 1. The data used for crater measurements are the Precision Experiment Data Records (PEDR) altimetry profiles (version L in IAU2000 coordinates) from the Mars Orbiter Laser Altimeter (MOLA) experiment on the Mars Global Surveyor spacecraft (Smith et al. 2001). Digital elevation maps (DEMs) are derived from the altimetry profiles using the Delaunay triangulation function as implemented in IDL (Lee and Schachter 1980). DEMs are generated on the fly at the user-specified spatial resolution for the region of interest. Individual outlier altimetry tracks may Ejecta and subsurface properties from crater geometries be removed interactively, with recalculation of the DEM. To speed the generation of DEMs, the time sequence of altimetry profiles was reorganized into regional data sets on a 25 × 25 global grid with no change to the original PEDR data. Each MOLA altimetry point covers a spot on the surface ~168 m in diameter with a vertical accuracy of ~1 m (Smith et al. 2001). Measurements on individual craters are conducted at a minimum resolution of 40 pixels per crater rim diameter and usually at a resolution at the limit of the 300 m along-track spacing of the MOLA altimetry profiles. The crater features measured by HMars are illustrated in Fig. 2. Measurements of the crater rim diameter, DR, and rim height, HR, are derived from altimetry tracks. Measurements relative to the preimpact surface, i.e., the surface depth, dS, and the surface diameter, DS, are based on the DEMs. All volume measurements are calculated from the gridded data set. The crater depth from the rim, dR, is the sum of the surface depth and rim height. For each crater, the user interactively estimates the location of the crater center and rim diameter. The userestimated crater center is refined by convolving a ring with the topography gradient rather than fitting irregular rims with an ellipsoid. The crater radius, rim uplift height, and their variances are calculated by interpolating the rim location and maximum rim height along each altimetry track that passes within a specified fraction (e.g., 0.8) of the estimated crater radius to the estimated center. The highest track point near the rim, three points downrange, and three points uprange are fit with two quadratic functions to determine the intersection between the crater wall and the outer rim wall. The highest track point, which usually underestimates the true rim height, and the number of rim points used to calculate the rim diameter and rim height are also recorded in the crater data record. The background preimpact surface is derived from userspecified tie points beyond the edge of the ejecta blanket, shown as “+” in the example in Fig. 1. The tie points are interpolated by Delaunay triangulation across the crater cavity and ejecta blanket. We investigated different fitting procedures to define the preimpact surface (e.g., polynomial surfaces), but found that the Delaunay triangulation performed well on the simulated crater data sets (discussed below). In uneven terrain, particular care is required by the user to select background surface tie points that best represent the overall slope and roughness of the region. In addition, users may limit the area of integration surrounding the crater by tracing an outer edge to the region defined by the background tie points (outer solid white line in Fig. 1). When a background feature, e.g., a preexisting crater, lies just beyond the edge of the ejecta blanket, the best tie points for the preimpact surface may surround both craters. The edge of region definition then excludes the preexisting crater from the measurements. An example of a crater DEM with the background subtracted and edge of region definition is shown in Fig. 3a. 1511 The cavity volume from the preimpact surface, VCavity, the cavity volume from the rim, VRim, and the volume of ejecta and uplifted material, VAbove, are derived by integrating the DEM over the areas interior and exterior to the crater rim, respectively. Preimpact features or sparse data coverage may interfere with an accurate measurement of the volume of ejecta and uplifted material surrounding the crater. In these cases, one or more user-defined pie-shaped wedges may be used to exclude an area from integration (Figs. 1 and 3). The volume in the excluded region is set to zero, and the total volume measurements are adjusted to compensate for the missing area assuming axisymmetry. Measurement of the continuous ejecta blanket volume, VEjecta, requires knowledge of the profile of the deformed preimpact surface surrounding the crater. For simple bowlshaped craters, the surrounding topography is uplifted and forms a portion of the crater rim. In this case, the preimpact surface can be estimated by connecting a fraction of the rim height with the surface beyond the ejecta blanket. For larger complex craters, however, the process of crater collapse may result in little or no deflection of the preimpact surface outside the final crater rim (Melosh and Ivanov 1999). In this study focusing primarily on fresh, simple, and transition-sized craters with well-preserved rims, we estimate the maximum possible uplifted surface by a power law that lies just beneath the ejecta profile, as shown in Fig. 2. Hence, VEjecta provides a lower limit to the ejecta volume, and VAbove is the upper limit. If regions of the ejecta blanket are excluded with pie-shaped wedges, the value of VEjecta is also increased accordingly. Figure 3b displays the volume of ejecta after subtraction of the model uplifted surface. The uplifted surface height, hu, is defined by hu = 0.5HR(r/RR)−5.5 based on Martian cratering simulations (Stewart et al. 2004), where r is the radial distance from the crater center and RR is the crater rim radius. The uplift profile is discussed in more detail in the “Comparison to Observations and Improving the Scaling Relationships” section. In the example in Fig. 3b, a portion of the eastern rim has been oversubtracted, producing negative thickness in the ejecta blanket. The irregularities in the rim and gaps in the altimetry data result in volume errors near the crater rim, which are small compared to the overall uncertainty in the uplifted surface. As part of the data entry for each crater, the user traces the distal edge of the ejecta and, if a double or multiple layer rampart crater, the inner rampart (inner solid white line in Fig. 1). The user consults the Viking Mosaicked Digital Image Model (MDIM) and Thermal Emission Imaging System (THEMIS) imagery to aid the identification of features in the ejecta blanket. The HMars toolkit also includes several visualization aids such as generation of 2-D and 3-D shaded relief views of the DEM, viewing of track coverage, and viewing of along-track altimetry profiles and arbitrary profiles across the DEM. Following the guidelines of Barlow et al. (2000), the user classifies each crater ejecta type as either i) no rampart, ii) single rampart, iii) double rampart or 1512 S. T. Stewart and G. J. Valiant Fig. 1. A screen shot of the HMars interactive crater measurement program. Fig. 2. Measured crater attributes in HMars. DR = rim diameter; DS = surface diameter; dR = rim depth; dS = surface depth; HR = rim height; VCavity = volume of cavity below preimpact surface; VRim = volume of cavity below the rim level; VAbove = volume of ejecta and uplifted material; and VEjecta = volume of ejecta only. Profile of a DR = 8 km crater from Lunae Planum (11.6°N, 72.4°E). Vertical exaggeration (V.E.) is 24:1. multiple ramparts, iv) radial ejecta, or v) severely eroded ejecta. The gradation state of each crater is described qualitatively by the user as i) fresh, ii) moderately degraded (partial infill/partial ejecta erosion), iii) heavily degraded (significant infill and/or ejecta erosion), or iv) severely degraded (complete infill and rim erosion). The approximately 30 crater measurements and observations are stored in a data structure with a keyword and value format. The final crater database is a list of data structures in ASCII format. The user-defined traces of ejecta ramparts are stored in a binary IDL save file. Each crater measurement also generates a log file, which records all the user inputs and interactive selections. Each crater measurement may be recreated using the log file, e.g., to test different fitting algorithms for the same user inputs or to investigate an anomalous measurement. Ejecta and subsurface properties from crater geometries 1513 Fig. 3. a) A shaded relief perspective view of the crater DEM from Fig. 1 with the preimpact surface subtracted. Zeroed pixels on the DEM lie beyond the user-defined edge of the region. V.E. = 15:1. b) DEM of ejecta blanket with uplifted surface subtracted. V.E. = 37.1. DEM resolution is 350 m/pixel. Resolution and Measurement Tests The HMars toolkit was tested for systematic errors, resolution limits, and reproducibility. We conducted measurements of simulated craters on different background terrains and verified that no systematic offsets were present in the crater measurements. The resolution limits and standard deviation for each measurement were derived by varying the altimetry track density on simulated craters on different terrains. The reproducibility was tested by comparison to published crater geometry measurements and comparing measurements by different users. DEMs of simulated fresh craters were generated on a flat background surface at spatial resolutions comparable to the MOLA track density at high latitudes. The simulated crater sets are composed of 15 craters each in 8 size bins between 2 and 50 km rim diameter. The model craters were created using analytic functions for the height of the uplifted preimpact surface, described by (r/RR)−5.5, and crater cavity depth, modeled with (r/RR)4. The rim uplift model was derived from impact cratering calculations (Stewart et al. 2004). The rim uplift is sensitive to the target material properties, and the uncertainties in applying a rim uplift model to the Mars data are discussed in the “Comparison to Observations and Improving the Scaling Relationships” section. The simulated rim height was typically several percent of the rim diameter for small craters, dropping to about a percent for the largest craters. The simulated crater depth was a third of the crater rim diameter for the smallest craters, and decreased to several percent for the largest craters. The simulated crater set contains both single- and double-layer rampart ejecta. The exact ejecta thickness and range were randomly generated over a realistic range of parameters. The single-layer ejecta thickness was given by sin2(r2) for r = 0 to π0.5. Double-layer ejecta are described by a piecewise function with an exponential section and ending with a r0.2 profile. Examples of the simulated craters are shown in Fig. 4. The simulated craters’ DEMs were sampled to create synthetic altimetry track profiles. Three track densities were simulated: the low track density used actual altimetry tracks from the equatorial region on Mars, the medium track density had double the number of tracks, and high track density had four times the number of tracks. The three synthetic altimetry track densities are representative of the MOLA data coverage at the equator, mid-, and high latitudes. Four background surfaces were created to test the HMars toolkit. We considered an ideal flat background and three backgrounds generated by tiling patches from regions on Mars with few large craters that are representative of the study regions in this work. The backgrounds are centered at i) 33°N, 200°E (Arcadia Planitia); ii) 7°N, 290°E (Lunae Planum); and iii) 32°N, 98°E (Utopia Planitia). Each background surface was also sampled to create altimetry profiles at each track density. For each crater size and at each track density, the altimetry profiles of the craters were added to the altimetry profile of each background. The resulting synthetic altimetry tracks were used to generate synthetic DEMs using the same interpolation scheme as used on the real MOLA data. Example simulated crater fields 10 km in rim diameter are shown in Fig. 4, where the right half of each region also displays the track density. 1514 S. T. Stewart and G. J. Valiant Fig. 4. Simulated craters 10 km in rim diameter on Martian background terrains. a) Lunae Planum with low track density, centered at 7°N, 290°E. b) Utopia Planitia with high track density, centered at 32°N, 98°E. The footprint of each altimetry track overlies the right half of each image. Identical measurements were performed on each background and at each track density at a DEM resolution of 0.3 km/pixel for DR ≥ 6 km and 0.15 km/pixel for DR ≤ 4 km. When there were overlapping ejecta blankets, the pie-wedge tool was used to exclude these areas from the crater measurements. The results from error analyses of the simulated crater measurements are shown in Fig. 5. The statistics presented in Fig. 5 demonstrate that the MOLA data resolves populations of craters with DR ≥ 6 km at all latitudes and at all track densities. Measurement accuracy begins to degrade when the number of points around the rim falls below about ten, although reasonable measurements are possible with as few as four rim points on DR = 4 km craters. Large-scale topographic features in the background surface that cannot be interpolated accurately provide the greatest source of error on volume measurements of craters with DR > 20 km. At all sizes, the definition of the preimpact surface is the largest source of error in volume measurements. For craters with DR ≥ 6 km, the errors in measurement of the rim radius drop from 5% to less than 1% (Fig. 5b). For this measurement, the largest contribution to the rim radius error is the natural variability in the crater radius. On very smooth terrains, the error in the rim height is only a few percent (Figs. 5c and 5d), but on rough terrains, the error can be as large as 10–20%. The rim depth of the crater has an error of only a few percent, except for large craters on rough terrains, Ejecta and subsurface properties from crater geometries where the error rises to ~10% (Fig. 5e). The systematic errors in measurements of crater depth and rim height are dominated by the track coverage and the interpolation scheme to fit the rim crest along altimetry tracks. Because the systematic errors are relatively small, errors in real crater measurements are usually dominated by natural variability along the rim crest. Crater cavity measurements are also very robust, with ~10% error for DR ≥ 6 km across all terrains (Fig. 5f). If crater infill is less than or comparable to 10%, then cavity volume measurement errors are dominated by the interpolated DEM, and the effect of infill results in a negligible bias to the measurements. If crater infill is >10%, then the original cavity volume should be estimated by a model of the cavity profile immediately after the impact event. The errors in the volume of ejecta plus uplifted surface, VAbove, are larger than VCavity because of the larger area of interpolation for the preimpact surface (Fig. 5g). Nevertheless, on smooth surfaces, the systematic errors are only ~20%, rising to ~60% on rougher surfaces. For simulated craters with known uplifted surface profiles, measurements of VEjecta are only slightly worse than for VAbove (Fig. 5h). Since the profile of the uplifted surface is not known for real craters, the uplifted surface subtraction remains the largest source of error for ejecta volume measurements. We emphasize that the simulated crater analysis demonstrates that at all latitudes, populations of craters as small as DR = 2 km that are fresh and axisymmetric may be measured with typical errors <30%. There is almost no latitude dependence in the accuracy of crater population measurements using the MOLA data set for fresh craters with DR ≥ 2 km. Hence, even though MOLA tracks rarely cross through the center of small craters, the general properties of a group of small craters can be measured coarsely. Errors in individual crater measurements are larger and harder to quantify. Confidence in individual crater measurements can be gauged by comparing measurements between different research groups who have developed their own measurement toolkits and strategies. Hence, to further validate the HMars toolkit, we examined craters with published geometric measurements by Garvin et al. (2000) and Barnouin-Jha et al. (2005). Our measurements are in very good agreement with the published values, within the expected errors (Fig. 6). Note that rimdiameter measurements by Barnouin-Jha et al. (2005) are systematically larger than those of HMars, which has been confirmed by the authors (Barnouin-Jha, personal communication). The good agreement between individual crater geometry measurements by researchers using independent analytical tools suggests that high confidence may be placed in MOLA-based measurements of single fresh craters in the size range validated here (~5–12 km rim diameter). Finally, we investigated the reproducibility of measurements of fresh craters by different users. 1515 Comparisons between three different users yielded results within the expected errors. While individual crater measurements may vary by up to ~2 times the systematic 1σ errors presented in Fig. 5, all data trends presented in this work are reproducible. Therefore, we conclude that measurements of fresh impact craters using the HMars toolkit and the MOLA data set yields robust results. Furthermore, the systematic errors derived from analysis of simulated craters provide a means to discern the significance of observed differences in crater geometries between different terrains. We propose that impact crater measurement techniques be validated using sets of simulated craters. The simulated craters used in this study are available from the authors. The crater sets are distributed as IDL “save” files, composed of simulated tracks (latitude, longitude) at three different track densities; simulated altitudes (for each crater size, track density, and background); and simulated crater DEMs on a flat background. The geometric parameters for each simulated crater are tabulated for comparison to measurements. The set of simulated fresh craters used in this study does not encompass all of the challenges encountered in crater measurements, e.g., degraded and severely nonaxisymmetric craters. We find that tailored simulated crater sets should be generated to gauge the accuracy of measurements of specific types of craters on representative background terrains for the goals of a particular investigation. CRATER MEASUREMENTS RESULTS Study Regions and Definition of Fresh Craters We present the results from a survey of the freshest craters in Utopia Planitia-Elysium Mons region (U.P.), Isidis Planitia (I.P.), Acidalia Planitia (A.P.), Lunae Planum (L.P.), and Solis Planum (S.P.) (Fig. 7). The lowland surfaces and highland plateaus are relatively flat, so highly accurate geometrical measurements comparing highland and lowland surfaces are possible. A total of 501 potentially fresh craters were measured in these regions; heavily degraded craters were not measured. Craters at latitudes above 55°N were not considered because we observed significantly shallower crater depths in this population of craters. Craters were grouped by regions restricted to 86–150°E, 18–55°N for U.P., 81–100°E, 5–25°N for I.P., 308–355°E, 27–54°N for A.P., 286–300°E, 6–20°N for L.P., and 261–284°E, 15–30°S for S.P. In order to make more robust regional comparisons, craters in Isidis and Utopia were further restricted by the boundaries of the impact basins, shown in Fig. 7, centered at 87.4°E, 12.7°N and 112°E, 45°N with 1090 and 2100 km diameters, respectively (Frey et al. 1998; Thomson and Head 1999). Fresh craters were identified quantitatively by deep crater 1516 S. T. Stewart and G. J. Valiant Fig. 5. Error analyses for HMars measurements of simulated crater data sets. a) Mean number of crater rim points. b) Normalized standard deviation (1σm/actual value) in interpolated rim radius. c) Normalized standard deviation in mean interpolated rim height. d) Normalized standard deviation in maximum rim height. Ejecta and subsurface properties from crater geometries 1517 Fig. 5. Continued. Error analyses for HMars measurements of simulated crater data sets. e) Normalized standard deviation in maximum crater depth from rim. f) Normalized standard deviation in crater cavity volume. g) Normalized standard deviation in volume above background surface (VUplift + VEjecta). h) Normalized standard deviation in volume of ejecta. 1518 S. T. Stewart and G. J. Valiant Table 1. Power law fits to deepest crater depths, dS, and tallest crater rims, HR, in each region, with all variables in kilometers. Strength regime (DS < 10 km) Acidalia Planitia F d = 0.276 D S S R Gravity regime (DS > 7 km) F d = 0.579 D S S R Number of craters (Equation 1, c = 0.65) 29 0.782 0.849 0.276 0.560 F H = 0.070 D S Utopia Planitia F d = 0.333 D S S R F H = 0.120 D S F d = 0.423 D S S R 0.714 0.569 0.386 0.483 53 F H = 0.115 D S Isidis Planitia F d = 0.342 D S S R F H = 0.165 D S F d = 0.503 D S S R 0.744 0.874 0.305 0.409 24 F H = 0.063 D S Lunae Planum F d = 0.197 D S S R F H = 0.234 D S F d = 0.395 D S S R 0.667 0.915 0.347 0.674 48 F H = 0.058 D S Solis Planum F d = 0.404 D S S R F H = 0.082 D S F d = 0.304 D S S R 0.375 1.045 0.417 0.337 33 F H = 0.046 D S F H = 0.153 D S cavities and high rims, and qualitatively by imagery of the crater and ejecta blanket. We focused on craters with 4 < DR < 50 km; however, in a few cases, smaller craters were well resolved and those were included in our analyses. We did not include larger craters (DR > 50 km) because these craters are systematically older, and thus more degraded, than the fresher, smaller craters. Because all craters have experienced some level of degradation, we formulate criteria to identify the least degraded craters in each region. Since erosion softens the topography of impact craters, rim height and crater depth are strong indicators of preservation. To identify the least degraded craters in each region, we fit the maximum crater depth, F d , and maximum rim height, F H , as a function of s R diameter with two power laws, one for the strength regime and one for the gravity regime. The maximum value functions, which define the peak values for crater depth and rim height in each region, are given in Table 1. At each crater diameter, the best-preserved craters are defined by the following criteria: dS > c F d s [DS, region], and HR > c F HR [DS, region], and NR ≥ 4 (1) NR is the number of rim points. The number of rim points criterion, based on the results from the simulated crater analysis, establishes a minimum confidence level in the crater measurements. Since degradation and natural variability result in a continuum of values for rim height and crater depth, a cutoff point must be chosen to identify the freshest craters. Moderately fresh craters are defined by a cutoff value, c, of 0.4 and the freshest craters by a cutoff value of 0.65. The cutoff values were chosen to identify a sufficient number of the least-degraded craters in each region for statistical analyses. The power law fits to the strength and gravity regimes overlap between DS = 7 and 10 km. In this transition region, a particular crater is included in a set with value c when it satisfies the cutoff criteria in either regime. In Fig. 8, moderately fresh craters are plotted in black symbols and the 188 freshest craters are shown in colored symbols. Note that some moderately degraded craters with deep cavities do not have tall rims and vice-versa. The number of craters in each region identified by c = 0.65 is given in Table 1. The paraboloid crater shapes and comparisons to previously identified fresh crater populations on Earth, Mars, and the Moon (presented below) validate the approach taken here to identify fresh crater populations using the MOLA data. Results: Resolved Differences between Fresh Highland and Lowland Craters We find large, resolved differences between the freshest craters in the lowlands and highlands in their depth, rim height, cavity volumes, and volume above the preimpact surface. In Fig. 9a, the crater depth from the maximum rim point follows a continuous trend for the highland regions (triangle, diamond), but the lowland regions (circle, square, cross) have a break in the trend around the transition crater diameter, 7–10 km. In this and all subsequent plots, the fit line color corresponds to the symbol color for each region and the power law fits are given in Table 2. The maximum dR is significantly larger for strength-dominated craters in the lowlands compared to the highlands. A similar trend is Ejecta and subsurface properties from crater geometries 1519 Table 2. Crater geometry fits with DR and dS in kilometers. Region Strength regime Gravity regime Maximum dR (km) (Fig. 9) Acidalia Planitia Utopia Planitia Isidis Planitia Lunae Planum Solis Planum dR = 0.72 0.302D R 0.79 1.02 Strength regime Gravity regime VAbove (km3) (Fig. 14) 1.91 2.43 3.06 3.54 2.13 V Cavity = 2.56 0.126D S 0.29 0.43 0.57 0.52 V Above = 0.443D S V Above = 0.220D S V Above = 0.175 D S V Above = 0.207 D S V Above = 0.296 D S V Above = 0.092 D S V Above = 0.216D S 2.44 2.44 2.23 2.54 d R = 0.288D R d R = 0.213D R d R = 0.261D R d R = 0.532D R d R = 0.738D R d R = 0.494D R d R = 0.266D R d R = 0.302D R ds (km) (Fig. 10) V Above = 0.041 D S V Above = 0.015 D S V Above = 0.167 D S 0.59 0.26 2.07 VAbove/VCavity (Fig. 15) Acidalia Planitia Utopia Planitia Isidis Planitia Lunae Planum Solis Planum dS = dS = 0.38 0.425D S 0.39 0.492D S 0.27 0.47 0.16 dS = 0.38 0.384D S 0.41 0.41 0.56 0.62 V Cavity = 1.24 ± 0.53 V Cavity = 1.48 ± 0.88 V Cavity = 0.91 ± 0.39 V Cavity = 1.07 ± 0.43 V Cavity = 0.83 ± 0.34 DR/DS (Fig. 17) DR DS = – 0.070 1.28D R – 0.107 – 0.047 – 0.025 – 0.070 0.013 – 0.014 – 0.016 – 0.012 0.011 d S = 0.404D S d S = 0.351D S d S = 0.192D S d S = 0.175D S HR (km) (Fig. 12) d S = 0.691D S d S = 0.238D S d S = 0.469D S Acidalia Planitia Utopia Planitia Isidis Planitia Lunae Planum Solis Planum HR = 0.54 0.082D R 0.70 1.03 0.57 0.75 HR = 0.62 0.072D R 0.33 0.51 0.47 0.23 D R D S = 1.09D R D R D S = 1.18D R D R D S = 1.18D R D R D S = 1.14D R D R D S = 1.08D R H R = 0.076D R H R = 0.036D R H R = 0.081D R H R = 0.065D R H R = 0.175D R H R = 0.133D R H R = 0.091D R H R = 0.138D R D R D S = 1.47D R D R D S = 1.29D R D R D S = 1.18D R D R D S = 1.28D R VCavity (km3) (Fig. 13) Acidalia Planitia Utopia Planitia Isidis Planitia Lunae Planum Solis Planum VEjecta (km3) (Fig. 21) 2.43 0.152 D S 2.20 2.39 2.52 2.66 V Cavity = 2.56 0.126D S 2.59 2.24 2.64 2.20 V Cavity = V Ejecta = 0.593D R V Ejecta = 0.644D R 1.40 1.36 3.71 V Ejecta = 0.102 D R V Ejecta = 0.102 D R V Ejecta = 0.080 D R V Ejecta = 0.052 D R V Ejecta = 0.082 D R 2.37 2.46 2.41 2.49 2.05 V Cavity = 0.132D S V Cavity = 0.258D S V Cavity = 0.077D S V Cavity = 0.288D S V Cavity = 0.161D S V Cavity = 0.098D S V Cavity = 0.072D S V Ejecta = 0.004 D R V Ejecta = 0.001D R V Ejecta = 0.023D R 4.49 2.82 V Cavity = 0.178 D S observed for the average value of dR. The trends for crater rim depths reported by Garvin et al. (2003) and Pike (1980) are in good agreement with the highland craters (Fig. 9b). The power law fits in all regions converge as crater diameter increases. For DR ~ 8 km, the lowland crater depths are factors of 1.5 to 2.0 times greater than crater depths in Lunae Planum and Solis Planum (Fig. 9c). Lunae Planum is used as a reference in Fig. 9 and subsequent figures to illustrate the differences between fits to the data, because the population of craters in this region exhibits the least scatter of the regions studied here. In all figures of crater measurements versus diameter, the reader should refer to the simulated crater analyses presented in Fig. 5 to estimate the magnitude of the measurement error, which varies a function of crater diameter and terrain. The systematic measurement error, σm, for crater depth has a normalized standard deviation of <5% in this size range (Fig. 5e), so the observations are ≥10σm above the measurement error. However, the dominant error is the observed natural variability in each crater population. The standard deviation in the rim height measurement around a single crater is on average about 15%, which contributes to about a 4% error on the rim depth (which corresponds to about 1–2 symbol radii in Fig. 9). The observed measurement scatter of rim depth around the fit in each region, σο, varies from 4.5% in Isidis to 10% in Solis Planum, so the observed regional differences are ≥5σo detections. Most of the difference in crater rim depths, dR = dS + HR, 1520 S. T. Stewart and G. J. Valiant Fig. 6. Comparisons of crater geometry measurements from HMars (horizontal axis) to published values from Garvin et al. (2000) (X’s) and Barnouin-Jha et al. (2005) (diamonds) (vertical axis). Solid line is one-to-one correspondence. Fig. 7. Freshest well-resolved craters with diameters ≥4 km in Utopia-Elysium region = circles, Isidis impact basin = squares, Acidalia Planitia = crosses, Lunae Planum = triangles, and Solis Planum = diamonds. Filled black dots denote moderately fresh craters. The background is shaded MOLA topography in Mercator projection. is found in the depth from the preimpact surface level, dS (refer to Fig. 2). In Fig. 10, the surface depths follow the same general trend as the rim depths, but in this case, the difference between highlands and lowlands become more pronounced as the crater diameter decreases. As shown in Fig. 8, the maximum rim uplift heights are more tightly clustered than the surface depths for fresh craters. The difference in crater shape is apparent in typical MOLA track profiles for fresh craters ~9 km in rim diameter in each region (Fig. 11). The natural variation in the background surface elevation around the crater contributes 2% to the error in the surface depth, or about 1 symbol radius in Fig. 10. Although the rim height variations are smaller than crater depth variations, there are a couple of general trends. In the strength regime, the rim heights are similar in each region (Fig. 12a). For all crater sizes, the rim heights for the freshest craters in this study are similar to measurements of fresh lunar craters (Pike 1977) but systematically taller than globally averaged values for fresh Martian craters from Garvin et al. (2003) (Fig. 12b). Simple lunar craters display a nearly Ejecta and subsurface properties from crater geometries 1521 Fig. 8.The freshest craters (colored symbols) are defined by having both (a) deep crater cavities and (b) tall rim heights (Equation 1 with c = 0.65. Utopia basin = circles; Isidis basin = squares; Acidalia Planitia = crosses; Lunae Planum = triangles; Solis Planum = diamonds. Smaller black symbols denote moderately fresh craters in each region (c = 0.4). Fig. 9. The depth of the crater from maximum rim point. a) Power law fits to freshest craters in each region. Fit-line colors correspond to the color of symbols in each region. b) A comparison of rim depthto-rim diameter ratio to previously published crater depth functions. c) The ratio of power law fits in (a) to Lunae Planum fit. Utopia basin = circles; Isidis basin = squares; Acidalia Planitia = crosses; Lunae Planum = triangles; Solis Planum = diamonds. Fig. 10. The depth of the crater from average preimpact surface elevation. a) Power law fits to freshest craters in each region. b) Power law fits to ratio of surface depth to surface diameter. c) Ratio of power law fits in (a) to Lunae Planum fit. Utopia basin = circles; Isidis basin = squares; Acidalia Planitia = crosses; Lunae Planum = triangles; Solis Planum = diamonds. Fig. 11. Typical MOLA altimetry profiles passing near the center of fresh craters (DR ~ 9 km). Each dot represents an individual MOLA measurement. Crater locations: L.P. (293.0°E, 18.6°N); S.P. (277.0°E, 19.5°S); A.P. (325.8°E, 39.1°N); I.P. (88.0°E, 12.8°N); U.P. (109.6°E, 34.2°N). The vertical exaggeration is 10:1. The sharp dip in the profile near the Lunae Planum crater rim is a small crater which formed on the ejecta blanket. 1522 S. T. Stewart and G. J. Valiant is compared to ideal values for a paraboloid with the measured surface diameter and cavity volume (Fig. 16). All of the freshest craters have nearly paraboloid cavity shapes. Note that the highland craters are systematically slightly shallower than the lowland craters. Another useful shape parameter is the ratio of the rim diameter to the surface diameter, DR/DS, shown in Fig. 17, which will be used in the “Comparison to Observations and Improving the Scaling Relationships” section. Gravity regime craters all display similar values of DR/DS ~ 1.12. Strength-dominated craters in the Isidis and Utopia basins have slightly larger values of DR/ DS, indicating that the upper crater walls are less steep compared to the other regions. In summary, we find large, resolved differences between the freshest highland and lowland craters. The observed differences are larger than the natural variability in crater shapes and larger than measurement error limits using the MOLA data. There are significant differences in the crater depth (≥5σo) and cavity volumes (≥2σo) for strengthdominated craters. Rim heights and VAbove are slightly larger in the lowlands for all crater sizes. There is a population of craters in Utopia with excess VAbove compared to VCavity, which is discussed in more detail by Black and Stewart (Forthcoming). All fresh crater cavities have paraboloid shapes. The observed values for DR/DS are slightly lower than the typical value of 1.3 used in crater scaling (Melosh 1989, 1998). The larger dR/DR values for simple lowland craters compared to simple highland craters suggest that the surface materials have different effective strengths, which is discussed in the next section. IMPLICATIONS FOR MARTIAN SUBSURFACE PROPERTIES The definition of the freshest craters in this study identifies the least degraded and preferentially youngest craters in each region. Here we make the reasonable assumption that the freshest craters in each region were formed by the same impactor population. The freshest lowland craters are Amazonian and the freshest highland craters likely formed in the same epoch. Therefore, systematic differences in the observed crater population should reflect differences in the properties of the target surfaces. Crater Scaling Relationships The difference in surface properties may be quantified through the use of cratering scaling relationships, such as the commonly used π-group scaling laws (Melosh 1989). In this section, we derive the crater diameter and volume as a function of the impact parameters and the effective surface strength. These relationships are used in the next section to investigate the difference in surface strength between regions on Mars. constant rim height to diameter ratio, HR/DR, of 0.036. The lunar complex crater HR/DR power law is corrected for the difference in gravity between the two bodies by a factor of gMoon/gMars. Our measurement of rim heights interpolates the altimetry profiles to fit for the rim crest, and the excellent correlation to lunar crater profiles lends support to this approach. The rim height measurements by Garvin et al. (2003) may be biased by regions with shallower rim heights in the global average and/or a different approach to measurement of the rim crest. In the gravity regime, the mean interpolated rim heights of lowland craters are slightly larger than highland craters, where the power law fits differ by factors of 0.25 to 0.75 (Fig. 12c). The toolkit measurement errors of rim height are typically 10% or less (Fig. 5c). The observed mean normalized standard deviation of crater rim heights is 17%, averaging over all regions, so the measurement errors are dominated by the natural variability around the crater rim (Fig. 12d). While the differences are not dramatic, there is a robust trend for taller rim heights in the lowland, gravity-dominated, freshest craters. Significant differences are observed in crater cavity volume, as expected from the variations in crater depth. Gravity-dominated craters in each region have similar cavity volumes, but lowland strength-dominated craters have significantly larger cavity volumes compared to highland craters of the same crater diameter (Fig. 13). Lowland craters with DS ≤ 10 km have volumes >50% larger than highland craters. Typical measurement errors of VCavity are ~10% (Fig. 5f); hence, the volume measurements are well resolved (≥5σm). The observed scatter around the power law fits vary from 20 to 30%, so the volume differences are resolved by ≥2σo. The volume of the crater cavity measured from the peak rim height, VRim, shows the same trend as VCavity. The combined volume of ejecta and uplifted preimpact surface, VAbove, is systematically larger in the lowlands for all crater sizes, with the exception of small craters in Isidis. VAbove is 50–100% larger for lowland craters (Fig. 14). The measurement error for VAbove ranges between 20 and 60% (Fig. 5g). The volume difference is best resolved for crater diameters between 6 and 20 km. By mass conservation and neglecting the effects of bulking, the volume of the crater cavity below the preimpact surface is approximately equal to the combined volume of uplifted and ejected material (Fig. 15). As expected, the average values for the ratio, VAbove/ VCavity, are nearly 1. Notably, there is a population of craters with an excess of material above the preimpact surface. Ratios greater than 2 are more than 2σ away from the average for all craters excluding Utopia. There is a slight size dependence in the VAbove/VCavity ratio for craters in S.P. and I.P., but the other regions have a nearly constant average value for all crater sizes. The relative contributions to VAbove from VUplift and VEjecta are discussed in the “Ejecta Volumes and the Z-EDOZ Model” section. To verify that the data set of freshest craters does not have anomalously shaped craters, the depth from the surface Ejecta and subsurface properties from crater geometries 1523 [ Fig. 12. Mean interpolated rim heights. a) Power law fits to freshest craters in each region. b) Comparison to published Martian global rim height functions and lunar craters (Pike 1977). c) Ratio of power law fits in (a) to Lunae Planum fit. d) Normalized standard deviation of rim point measurements. Utopia basin = circles; Isidis basin = squares; Acidalia Planitia = crosses; Lunae Planum = triangles; Solis Planum = diamonds. Fig. 13. Cavity volume below preimpact surface. a) Power law fits to freshest craters in each region. b) Ratio of power law fits in (a) to Lunae Planum fit. Utopia basin = circles; Isidis basin = squares; Acidalia Planitia = crosses; Lunae Planum = triangles; Solis Planum = diamonds. Nondimensional scaling parameters, crater diameters, and volumes are defined by: 1.61gD P π 2 = -------------------2 vi Y π 3 = -------2 ρv i ρ 1/3 π D = D ST # ------$ ! m P" ρ π V = V T -----mP where π2 is the inverse Froude number, π3 is the inverse Cauchy number, DP is the projectile diameter, g = 3.72 m s−2 is the gravitational acceleration on Mars, and vi = 10 km s−1 is the average asteroid impact velocity on Mars (Bottke et al. 1994). The density of the projectile and surface are assumed to be equal with ρ = 2800 kg m−3, and mP is the mass of the projectile. VT is the volume of the transient cavity below the preimpact surface, DST is the diameter of the transient cavity at the preimpact surface, and Y is the effective strength of the surface. Fig. 14. The volume above the preimpact surface (combined ejecta and uplifted surface volume). a) Power law fits to freshest craters in each region. b) The ratio of power law fits in (a) to Lunae Planum fit. Utopia basin = circles; Isidis basin = squares; Acidalia Planitia = crosses; Lunae Planum = triangles; Solis Planum = diamonds. 1524 S. T. Stewart and G. J. Valiant Fig. 15. The ratio of volume above preimpact surface to cavity volume with mean values (solid lines) and log-linear fits to U.P. and S.P. (dashed lines). Utopia basin = circles; Isidis basin = squares; Acidalia Planitia = crosses; Lunae Planum = triangles; Solis Planum = diamonds. Fig. 16. The depth from the surface, dS, of the freshest crater cavities compared to expected depth for paraboloid cavity with observed volume, VCavity, and surface diameter, DS. Solid line is a one-to-one correlation. Utopia basin = circles; Isidis basin = squares; Acidalia Planitia = crosses; Lunae Planum = triangles; Solis Planum = diamonds. has an adjustment for the impact angle from the surface, φ, which is 45 degrees on average: – β ρ – 1/3 – 1/3 D ST = C 2 π 2 # --- $ ( sinφ ) ! m" (2) Fig. 17. The ratio of rim diameter-to-surface diameter for freshest craters with power law fits. Utopia basin = circles; Isidis basin = squares; Acidalia Planitia = crosses; Lunae Planum = triangles; Solis Planum = diamonds. Following the method of Melosh (1998), transient craters in the strength regime may be approximated by the π2 gravity scaling law. The final rim diameter separating the strength and gravity regimes may be scaled from the observed lunar transition diameter, D * Moon , by compensating for the differences in gravity and density of the surface: * D* Mars = D Moon [ ( ρg ) Moon ⁄ ( ρg ) Mars ] . Using ρMoon = * −3, g −2 2700 kg m Moon = 1.67 m s , and D Moon = 18.7 km, * D Mars = 8.1 km. During the final stage of crater formation, simple craters widen by crater wall slumping by a factor of CD = 1.2 (e.g., DS = CDDST), and the final rim diameter is larger than the surface diameter by a factor of CR = 1.3 (e.g., DR = CRDS) (Melosh 1989, 1998). Then, D R,Strength = C D C R D ST = 1.56D ST and D R,Gravity = ( D R,Strength ) 1.18 First we present the formulae for the transient and final crater diameter using the gravity scaling parameter π2. In the gravity regime, the surface diameter of the transient cavity scales as: πD = C2 π2 –β ⁄ ( D* ) 0.18 (3) where C2 = 1.6 and β = 0.22 for competent rock (Schmidt and Housen 1987). The scaled transient crater surface diameter where the scaling for the final rim diameter in the gravity regime is derived from the lunar crater population (Melosh 1998). The π2 scaling results for the transient and final crater diameters for a fixed impact velocity and variable projectile size are shown as the solid lines in Fig. 18. Alternatively, in the strength regime, the transient crater may be scaled by π3, where: Ejecta and subsurface properties from crater geometries 1525 πD = C3 π3 , # Y D ST = C 3 % --------& ! ρv 2" i $ –σ ρ # --- $ ! m" – 1/3 –σ 2 π V T = C 3 # -- $ ! 6" 2/3 π ( sinφ ) 1/3 . (4) – 2σ – µ/2 3# Y $ sin φ -- KD P % --------& 8 ! ρv 2" i (7) The scaling exponents are related to the coupling exponent, µ, of Holsapple and Schmidt (1987) by β = µ/(2 + µ) and σ = µ/ 2. The exponent σ is related to β by σ = β/(1 − β), so for competent rock σ = 0.28 for β = 0.22. The effective yield strength of the surface may be estimated using the transition crater rim diameter derived from cratering simulations by O’Keefe and Ahrens (1993), Y = 0.11D*ρg. D* (5) For = 8.1 km, Y = 9 MPa. The strength scaling curve intersects the gravity scaling curve at D* using a coefficient C3 = 0.88 for competent rock, which is shown as the dashed lines in Fig. 18. Note that, as expected, the diameters of strength-dominated craters are slightly smaller than those predicted from gravity scaling (solid lines). Because the projectile size is not an observable parameter, we present equations for the transient crater depth and volume and discuss collapse to the observed final crater shape in the next section. In the strength regime, the depth and volume of the transient crater cavity depend on the effective strength of the surface. The maximum depth from the preimpact surface of the transient cavity is related to the yield strength by (O’Keefe and Ahrens 1993): # Y $ – µ/2 1/3 ( sinφ ) . d ST = KD P % --------& 2 ! ρv i " (6) In Equation 7, dST was replaced by Equation 6. Note that according to coupling theory, σ = µ/2, so the exponent to the inverse Cauchy number is −3σ. The fitted values of µ/2 to simulations and σ to experiments are nearly equal. Here, we assume that the diameter and depth scaling constants (C3, σ, K, µ) are weakly dependent on the yield strength. The variation in C3 and σ is small over a wide range of materials (water, sand, and competent rock) (Schmidt and Housen 1987), and K and µ are constants fit to a 4 order of magnitude range of the inverse Cauchy number (O’Keefe and Ahrens 1993). Then, from Equation 7, VT ~ Y−3σ, DST ~ Y−σ, and dST ~ Y−σ. As the yield strength increases, the depth and diameter of the transient cavity decrease at a similar rate, and the overall cavity volume decreases by approximately V ~ 1/Y for σ ~ 1/3. In Fig. 18, the dotted line illustrates a 33% increase in the effective yield strength of the surface (Y = 12 MPa). The stronger surface results in smaller craters compared to the 9 MPa surface (dashed line) for the same impact parameters. Note that the final crater diameter curve is offset from the gravity regime curve although the transient crater diameter intersects the gravity regime curve. An effective yield strength of 12 MPa leads to a strength to gravity regime transition rim diameter of 10.5 km, which is around the upper end of the observed range of transition diameters on Mars. Regional Differences in Simple Crater Collapse To compare the scaling laws to the crater geometry measurements in the strength regime, the volume and depth of the final simple crater must be related to the transient crater cavity. We assume that the collapse process is similar for all simple craters in the size range considered here (about 4– 10 km surface diameter). Crater collapse decreases the cavity volume and increases the cavity diameter, such that D S = C D D ST V Cavity = C V V T (8) Here the dependence on impact angle is stated explicitly. K = 0.42 is derived from cratering simulations, and the fitted exponent µ/2 = 0.26 is reduced slightly from O’Keefe and Ahrens’ value of 0.28 to intersect more recently calculated 2 transient cavity depths of 6DP at Y ⁄ ( ρv i ) ~ 3 × 10−5, which is typical of Martian impacts (Stewart et al. 2004). Note that the value of µ = 0.52 fitted to calculations agrees very well with the range inferred from experiments (Holsapple and Schmidt 1987; Schmidt and Housen 1987). Using the volume of a paraboloid, V = πdD2/8, the scaling for the transient cavity volume is derived from Equation 4: #πd $ -! 8 ST" – 1/2 1/2 VT # Y $ – σ ρ – 1/3 – 1/3 = C 3 % --------& # --- $ ( sin φ ) ! m" 2 ! ρv i " Here VT is defined by Equation 7 and DST is defined by Equation 4. Since both the transient and final crater cavities are approximately paraboloid, CD and CV are related by: 2 2 π π -- d S D S = C V -- d ST D ST 8 8 2/3 # Y $ – 2σ 2 ρ – 2/3 # π -- d ST$ % --------& V T = C 3 # --- $ ( sin φ ) ! m" !8 " ! 2" ρv i dS 3 2 K π – 1/3 ------ D S = C V ----- # -- $ - D ST D ST C 3 ! 6" DS 1526 S. T. Stewart and G. J. Valiant Fig. 18. π-scaled transient surface diameter (gray lines) and final rim diameter (black lines) in the gravity regime (solid lines and inset) and strength regime (effective strength: Y = 9MPa is dashed; Y = 12 MPa is dotted) as a function of projectile diameter, DP, for 10 km s−1 impacts at 45 degrees. dS 3 2 ------ D S = C V ( 0.59D ST )D ST DS dS 3 3 3 ------ C D D ST = 0.59C V D ST DS 0.59C V C D = ----------------dS # ------$ ! D S" 1 -3 Fig. 19. Crater diameter and volume fits (Equations 8 and 9) for simple highland craters (black lines) and simple lowland craters (blue lines). The offset between the highland and lowland craters is a result of differences in both the crater cavity shape and crater collapse processes. The dashed red line indicates the shift to smaller volumes in the lowland fit (the solid blue line) for a change in effective yield strength from 9 to 12 MPa. (9) The derivation includes the substitution from Equations 6 and 7 of dST = 0.59DST for K = 0.42, C3 = 0.88, and σ = µ/2. Note that the coefficient relating dST and DST is weakly dependent on the effective yield strength. Scaling the final crater volume by 0 < CV < 1 includes collapse of the transient crater to the final crater shape and bulking of the slumped material. Reasonable values of CV will not be near 0 or 1. CD is constrained to be a modest increase of the transient surface diameter to the final surface diameter for simple craters based on experiments and simulations (e.g., 1.1 < CD < 1.3). Fits to the observed values of dS/DS are shown in Fig. 10. Because dS/DS is not constant, CD and CV cannot both be constant. Alternatively, the final diameter and cavity volume may have a weak power law dependence on the transient cavity diameter and volume in Equation 8. Using the dS/DS fit to Solis Planum (Table 2) in Equation 9, we fit CD and CV in Equation 8 to match the crater geometry measurements, shown in Fig. 19. The crater size scaling parameters are the same as used in the Y = 9 MPa case presented in the previous section (dashed line in Fig. 18). The highland data are fit with CV = 0.35 and 1.1 < CD < 1.3 for increasing DS (solid black line), CD = 1.2 and 0.5 > CV > 0.25 (dashed black line), or a continuous function of CD and CV over a similar range of values. Using the dS/DS fit to Utopia (Table 2), the lowland craters cannot be fit with the same CV = 0.35 derived for the highlands because the value of CD becomes nonphysical (0.95 < CD < 1.3) by Equation 9. CD cannot be less than 1 by definition. Hence, CV must be larger for lowland craters than highland craters. The lowland data are fit with CV = 0.6 and 1.1 < CD < 1.3 (solid blue line in Fig. 19) or CD = 1.2 and 0.7 > CV > 0.45 (dashed blue line) for increasing DS. Note that when the shape of the crater is constant (dS/DS = constant), changing the effective yield strength of the surface maintains the slope of VCavity versus DS in Fig. 19. For the same impactor population, decreasing the yield strength shifts the population of craters to larger volumes (up and to the right along the same line), while increasing the yield strength shifts the population to smaller volumes along the same line (Equation 7). Therefore, the offset between the highland and lowland craters in Fig. 19 must reflect the different dS/DS ratios for the highland and lowland populations. In Fig. 19, the dashed red line illustrates a 33% increase in the effective yield strength compared to the solid blue line. The increase in CV required to fit lowland data demonstrates that the freshest craters in the lowlands have experienced less transient cavity collapse compared to the freshest highland craters. Craters formed in stronger materials collapse less, suggesting that the lowland surfaces are Ejecta and subsurface properties from crater geometries stronger than the highland surfaces in this study. Craters in the lowlands preserve 45–70% of the transient crater volume compared to 25–50% in the highlands. The larger depth-todiameter (dS/DS) and rim height-to-diameter (HR/DR) ratios in the lowlands reflect stronger surface materials and less crater collapse compared to the highlands. Concurrently with our own work, Boyce et al. (2006) also inferred the presence of high target strength materials in Isidis and southwestern Utopia from crater depth-to-diameter relationships. Comparison to Simple Terrestrial and Lunar Craters The rim depth-to-rim diameter ratio (dR/DR), rim heightto-rim diameter ratio (HR/DR), and crater volume collapse coefficient (CV) for the freshest simple Martian lowland craters are similar to those for fresh simple craters on Earth. Estimates of the breccia lens and transient crater volumes at Meteor Crater (Arizona, USA) and Brent Crater (Ontario, Canada) indicate that CV ~ 0.5 (Grieve and Garvin 1984), in excellent agreement with the CV values for the Martian lowlands and at the upper limit for the range of CV values for the highland craters. The average value of dR/DR for Martian lowland craters with DR < 8 km is 0.20 ± 0.02(1σ) (Fig. 9b), which is similar to the average (dR/DR) = 0.15 ± 0.02 of six relatively uneroded terrestrial craters in the size range 0.39 ≤ DR ≤ 2.44 km (Fudali et al. 1980). Fresh simple lunar craters have dR/DR ~ 0.196, and the dR/DR value is independent of crater size in the strength regime (Pike 1977). The rim height ratio HR/DR for Martian lowland craters with DR < 8 km is 0.041 ± 0.005 (Fig. 12b). The six fresh terrestrial craters studied by Fudali et al. (1980) have very similar values, HR/DR = 0.046 ± 0.011. Measurements of fresh simple lunar craters yield HR/DR ~ 0.036, nearly independent of crater size (Fig. 12b) (Pike 1977). The nearly constant dR/DR and HR/DR ratios for simple craters is indicative of strength-dominated crater shape, where the gravity on the three different planetary bodies had a minor role during crater collapse. However, the fresh simple craters in the Martian highland plains do not have a constant dR/DR ratio. Instead, the dR/DR trend is similar to an extrapolation from the fits in the gravity regime. Three lines of evidence indicate that the surface materials are substantially stronger in the Martian lowlands compared to the highlands: 1) The observed fresh simple crater shapes indicate that craters in the Martian lowlands experienced less crater collapse compared to highland craters, and the lowland values of 0.45 < CV < 0.7 are similar to those of one of the best-preserved craters on Earth, Meteor Crater; 2) The HR/DR ratio for Martian simple craters is weakly independent of crater size. The ratio is similar to fresh lunar and terrestrial craters, and craters in Utopia and Isidis have higher crater rims in the 7 < DR < 11 km range compared to the other study regions; and 3) The dR/DR ratio is substantially larger for simple craters in the lowlands (by a 1527 factor of 1.5 to 2.0), nearly independent of crater size, and similar to fresh craters on Earth and the Moon. By contrast, the inconstant dR/DR ratios in the highlands follow extrapolation from the gravity regime. Regional Differences in Effective Surface Strength Quantifying the difference in the effective strength between the highland and lowland surfaces is difficult. One approach that provides a broad estimate for the strength of a planetary surface is to examine the simple-to-complex crater transition size. Using 1/g scaling from the Moon, the simple to complex transition is expected to be around 8 km rim diameter. However, Viking-based studies of 11 interior characteristics of 230 craters indicate that the global average is significantly smaller: 5.1 +3.1 km (Pike 1988). More recent – 1.9 measurements of ~6000 craters using the MOLA data suggest (primarily on the basis of accurate dR/DR measurements) that the global mean is centered at about 7 km with significant regional differences (Garvin et al. 2003; Garvin et al. 2002). Note that the simple and complex fits to dR/DR by Garvin et al. (2003) (Fig. 9) intersect at 5.4 km, closer to the 3.1 +0.7 km intersection between small and large crater depth – 0.5 to diameter fits derived by Pike (1988). The primary difficulties in determining the true depth around small craters using the MOLA data are the rare alignment of an altimetry profile through the centers of small craters and the ~300 m along track spacing of the altimetry points. Therefore, it is difficult to use the MOLA data alone to identify depth-todiameter transitions around small craters. Several research groups have noted that there are lowland regions with large simple craters up to DR ~ 10.2 km (Boyce et al. 2005a; Boyce et al. 2006; Garvin et al. 2000; Pike 1980), significantly beyond the average transition diameter. In this study, measurements of dR, HR, and VCavity indicate that the freshest lowland craters have a transition diameter as large as 10 km, with a transition range of 8–10 km. In the highland regions of this study, there is no break in crater parameters similar to the obvious transition in the lowland craters using the MOLA data. For these regions, we rely on previous work by Pike (1980) to estimate the transition diameter. First, note that the depth-to-diameter ratio for our highland regions follow the global trends established using the Viking imagery data (Fig. 9). Thus, the deep, simple craters identified in Isidis and Utopia are anomalous compared to the global crater population. Our measured highland crater data sets show no transition in the wellresolved size range ≥6 rim diameter and indicate that the transition size is less than 6 km. Pike’s global estimate of a 5.1 km strength to gravity transition rim diameter and 3.1 km intersection between power law fits in each region are presently the best estimate for the transition diameter of the highland regions in this study. Pike (1988) suggested that the global simple-to-complex 1528 S. T. Stewart and G. J. Valiant basement of Utopia basin has very similar properties to the Noachian infill, the cratering record supports a thick layer of infill rather than a thin (6.5 km, based on the depth of the transient cavity for final 10 km rim diameter craters at a mean impact velocity for Mars (Fig. 18). Boyce et al. (2006) suggest that an extensive, ~1–2 km thick, olivine-rich mafic-to-ultramafic rock unit identified by thermal emission spectroscopy in the Nili Fossae region (Hamilton and Christensen 2005; Hamilton et al. 2003; Hoefen et al. 2003) may be responsible for the relatively strong subsurface implied by the crater geometries in Isidis and Utopia. However, unless the olivine-rich unit is much thicker than currently observed, it alone cannot account for the crater observations. In addition, Hesperian volcanic flooding and sedimentary (Vastitas Borealis Formation) deposits and Amazonian volcanic deposits from Elysium form shallow layers within the upper 1 km of the Utopia basin (Thomson and Head 2001). The size range of observed fresh crater forms cannot be controlled by these young surface deposits. Instead, the stronger layer controlling the shape of 4 to 10 km surface diameter Utopian craters must be the primary sedimentary or volcanic fill that was deposited in Utopia near the end of the Noachian period. The Utopia basin has a positive free-air gravity anomaly similar to lunar mascons (Smith et al. 1999), and geophysical models of the primary infill in the basin are best fit with a layer of material about 20 km thick (Searls and Phillips 2004). Unless the Ejecta and subsurface properties from crater geometries study are not shallower than the lowland craters because of the presence of a buried stronger layer, an effect that has been seen on the Moon (Quaide and Oberbeck 1968) and in laboratory experiments (Piekutowski 1977). Because the depth-to-diameter ratios for highland craters follow a similar power law over the size range 3 to 50 km rim diameter, there is no observational evidence for a strong buried layer effecting the measurements of the freshest craters in this study. Although the estimates of surface strength provided here are coarse, analyses of crater geometries can contribute quantitative information about the regional geology on Mars. Using crater geometries to interpret surface strength relies upon the preservation of some fresh craters on Mars and the assumption of similar impactor populations for all fresh craters. Partial or differential degradation of the freshest craters at different locations on Mars would contribute to the error in determining the transition diameter between strengthand gravity-dominated craters and estimates of the magnitude of crater collapse. The inferred values for the collapse coefficient, cD, could be interpreted as lower limits, as gradational effects reduce the volume of the crater cavity. However, because the inferred collapse coefficient agrees well with fresh terrestrial craters, no systematic degradation effects are resolved by our analyses, and we infer that the preservation level of the freshest population of craters on Mars is suitable for geometrical analyses. Because the differences in geometries between craters in the highland and lowland regions are much larger than the observed natural variability in the measured crater parameters (which includes partial and heterogeneous degradation), we conclude that the inferred differences in crater collapse and target strength are resolved. As our understanding of crater collapse improves, the crater collapse volume may be a more useful parameter for inferring the surface strength in future studies. As an aside, we note that the inference of strong and weak surfaces from the distribution of Martian meteorites (Head et al. 2002) applies to surface layers and not subsurface layers as discussed above. Also, the shape of fresh crater cavities is not necessarily correlated with the youngest surface deposits, and care must be taken when associating crater geometric properties with the geologic terrain. EJECTA VOLUMES AND THE Z-EDOZ MODEL Almost all fresh craters in this study exhibit fresh rampart-type ejecta blankets. Ejecta volume measurements can be used to constrain theories for the formation of fluidized ejecta. The dynamics of Martian ejecta blankets have been compared to terrestrial debris flows and landslides (BarnouinJha et al. 2005; Ivanov 1996) and also to gas-supported flows similar to volcanic base surges (Mouginis-Mark 1981; Schultz 1992; Schultz and Gault 1979). In this section, we address the appropriate scaling for the amount of ballistic 1529 ejecta required to form the observed ejecta blankets around fresh Martian craters. We also assess the effects of incorporating surface materials into the ejecta flow and concentration of material in the continuous ejecta blanket by nonballistic processes. The Z-EDOZ Model Maxwell’s Z model is an empirical, analytical model used to describe the excavation flow field of the transient cavity (Maxwell 1973, 1977). This simple, predictive model has been used to describe the starting conditions for calculations of the ejecta blanket flow, assuming ballistic ejection (e.g., Ivanov 1996), and the collapse of the transient cavity to the final crater shape (Melosh and Ivanov 1999). Here we compare the observed ejecta and uplifted surface volumes to the predictions of the Z model. In the original Z model, target material emerging at a radius, rS, follow spiral trajectories from the impact point according to: r = r S ( 1 – cos θ ) 1 ⁄ (Z – 2) where r and θ are polar coordinates with the θ = 0 axis pointing vertically down. rS is the radial distance from the impact point along the preimpact surface and the point where the trajectory cuts the surface. The excavation flow is described by a collection of trajectories or stream tubes where rS varies from 0 to the radius of the transient cavity, RST. The value of Z defines the curvature in the flow field. Realistic limits on crater-scaling laws limit the range of Z to greater than 1.5 and less than 3.0 (Holsapple and Schmidt 1982). Comparisons to experiments and calculations show that 2.5 < Z < 3.0. For a constant value of Z, the ejection angle at the preimpact surface is the same throughout the transient cavity with a value of ϕ = arctan (Z − 2). Major refinements to the original Z model include centering the flow below the impact surface and varying the value of Z with time during the excavation of the transient cavity (Croft 1980; Thomsen et al. 1979; Thomsen et al. 1980). The depth of the effective center of Z model flow (EDOZ) is about one projectile diameter. Due to the momentum of the projectile, impact craters are similar to explosion craters with the charge buried at shallow depths (Melosh 1989). The depth of burial corresponds physically to the isobaric core in the pressure field generated by the shock wave from the impact (Ahrens and O’Keefe 1987; Pierazzo and Melosh 2000; Pierazzo et al. 1997). Moving the center of flow beneath the surface results in an increase in the total excavated volume and higher ejection angles. A schematic of the Z-EDOZ flow field is shown in Fig. 20a. Recent comparisons between the Z-EDOZ model and ejection angles from oblique impact experiment show that the early time flow field for 45° impacts deviates slightly from the axisymmetric 1530 S. T. Stewart and G. J. Valiant model (Anderson et al. 2004). At late times during crater excavation, the component of momentum-driven flow from the projectile trajectory is negligible and the assumption of axisymmetry provides a good comparison to azimuthally integrated crater measurements. In the constant Z-EDOZ model, the flow field originates from a constant depth, dZ, throughout the excavation process. Croft (1980) found that Z decreases with time and reaches a constant value between 2.5 and 3.0 at about half the transient cavity radius. Properties of the constant Z-EDOZ flow field are derived by Croft (1980; Appendix A.). The Z-EDOZ flow field intersects the preimpact surface with ejection angles that decrease with increasing radius: 2 tan∆ ----------- – tan ∆ + ( Z + 2 ) cos∆ ϕ = arctan ------------------------------------------------------1 ----------- – tan∆ ( Z – 1 ) cos∆ where ∆ is the angle between the horizontal plane at depth dZ and the point where the stream tube intersects the preimpact surface (Fig. 20a). When dZ = 0, ∆ = 0 and the ejection angle is reduced to a constant ϕ = arctan(Z − 2) as in the original Z-model. The angle ∆ needs to be solved numerically for the Z-EDOZ trajectory emerging at each rS, and the value of ∆ decreases with increasing rS. The excavated volume outlined by an Z-EDOZ trajectory emerging at rS is given by: π d Z 2π ( Z – 2 ) 1 + sin∆ 3 - V Z,r = -- ----- + ----------------------- ------------------- r S . S 3 r S 3 ( Z + 1 ) cos 3 ∆ Example Z-EDOZ trajectories are shown in Figs. 20c and 20e. Note that the Z-EDOZ stream tubes reach greater depths compared to the original Z model (Figs. 20b and 20d). The ejection velocity, ve, for each stream tube is usually described by the power law: R ST m v e = v ST # -------- $ ! rS " where 2.5 < m < 3.0 (Ivanov 1996), and vST is the velocity at the radius of the transient cavity, RST. By equating the kinetic energy in a stream tube with the potential energy, Ivanov derives v ST = kgR ST , where g is the gravitational acceleration and k = [(Z − 2)(Z + 1)]/[Z(Z + 2)]. For the size range of craters considered here, vST is typically several 10 s m s−1. The Z and Z-EDOZ model ejecta thickness around the transient crater is calculated assuming ballistic trajectories and shown as the dashed lines above the uplifted surface in Figs. 20b–e. The final crater rim radius, RR, is calculated from Equation 3 and shown as the vertical arrows in Figs. 20b–e for final crater radii of 2 and 12.5 km. Uplifted and ejected Fig. 20. a) A schematic of the Z-EDOZ model variables: dZ = depth to center of flow; θ, r = polar coordinates along flow trajectories; rS = horizontal distance where single trajectory meets surface; ∆ = angle from center of flow to surface at distance rS; Rh = distance from flow center to uplifted surface point (ru,hh); δ = angle from flow center to uplift point; RST = surface radius of transient crater. b–e) Z and Z-EDOZ model excavation zones (dashed lines), uplift profiles (solid line above surface) and ejecta thickness (dashed lines above uplift profiles) for final crater radii of 2 km and 12.5 km, normalized to RST. The thick line outlines the transient crater cavity, and the vertical arrow locates the final crater radii, RR (Equation 3). The inner and outer grey trajectories identify the source of ejecta reaching 3RR and 1RR, respectively. Ejecta and subsurface properties from crater geometries materials within RR collapse back into the final crater cavity. In Fig. 20, dimensions are scaled to the transient crater surface radii, RST. The thick solid line outlines the paraboloid transient crater cavity, and the dashed trajectory outlines the excavated region, where rS = RST. Comparison to Observations and Improving the Scaling Relationships The median radius of the distal edge of the continuous ejecta blankets measured around the freshest craters is about 3RR (mean 3.5RR with a 1σ scatter of 0.9RR). Here, we define the volume of ballistic ejecta in the continuous ejecta blanket, VBE, as the volume of ejecta that lands between 1RR and 3RR. In Fig. 20, the inner grey stream tube locates the source of ejecta reaching a distance of 3RR, and the outer grey stream tube reaches a distance of 1RR. The continuous ejecta volumes measured for the freshest craters in this study are presented in Fig. 21 with the ballistic ejecta volumes for Z = 2.7 for both the Z (black line) and Z-EDOZ (green lines) models. The volume of ejecta predicted by each model is significantly less than that inferred from the data. The original Z model is a very poor fit to the data. The Z-EDOZ model, with dz ~ 0.9DP using π2 scaling, is a better fit (solid green line). However, only with a very generous bulking factor of 2 (dashed green line) does the Z-EDOZ model approach the lowest ejecta volume measurements. The ballistic ejecta with a range between 1RR and 3RR is only a fraction of the total excavated volume, VEXC. Using the Z-EDOZ model with Z = 2.7 and dz ~ 0.9DP for final craters with diameters between 3 and 50 km, the fraction of ejecta emplaced ballistically near the transient crater rim at radii within the collapsed final crater ( RST. The range, ru, the uplift height, hu, are given by: hu = Rh sinδ − dZ ru = Rh cosδ Rh is the distance between the center of the flow, dZ, and the uplift location (ru, hu), and δ is the angle from the horizontal plane at depth dZ and point (ru, hu) (Fig. 20a). The angle δ is defined by: (Z + 1) ---------------(Z – 2) (Z – 2) ( Z + 1 ) ------------------------------- ( Z + 1 ) (Z – 2) sin δ = ( 1 – cos θ i ) + ( 1 – sin ∆ ) –1 where θi is the polar angle to the point where the stream tube intersects the transient cavity, e.g., at the intersection of the dot-dashed line and thick solid line in Fig. 20. Here, we assume the transient cavity is parabolic with a depth of 0.94RST (Equation 6). See Croft (1980) for a detailed derivation of the Z-EDOZ parameters. VUp is the volume of uplifted material between 1RR and 3RR. The total volume above the preimpact surface between 1RR and 3RR predicted by the Z-EDOZ model, VZ−Above, is the sum of VUp and VBE. The measured VAbove for the freshest Martian craters is compared to the Z and Z-EDOZ model predictions in Fig. 22a. Both models predict less volume above the preimpact surface than observed. Bulking the ballistic ejecta volume by a factor of two (VUp +2VBE) also falls short of the data. Around small simple craters, the near-rim ejecta will experience moderate bulking from fracturing, similar to the 5–10% bulking observed in the breccia lens within terrestrial crater cavities (Grieve and Garvin 1984). Distal ejecta around simple craters and ejecta surrounding large complex craters are heavily fractured and incorporate surface materials into the ejecta flow through ballistic sedimentation (Oberbeck 1975). The ejecta around the 1.8 km diameter Lonar crater, which formed within the Deccan Traps, India, displays features similar to the ground-hugging flow around Martian craters (Fudali et al. 1980; Stewart et al. 2005). At Lonar, the volume of the distal ejecta is about 30–40% larger than the bulk volume of the undisturbed country rock. Thus, increasing the total ejecta volume by a factor of 2 is a generous upper limit. The relative contribution of VBE to VZ−Above is shown in Fig. 22b. The volume of the continuous ejecta blanket around Martian craters was measured using an uplift profile of the form h = 0.5HR(r/RR)−5.5. The Z and Z-EDOZ models predict lower values of VEjecta/VAbove than observed, and the uplift 1532 S. T. Stewart and G. J. Valiant Fig. 21. a) Power law fits to the estimated ejecta volume. b) Ejecta volume predictions for original Z model (solid black line), Z-EDOZ model (solid green line), and Z-EDOZ model with volume increased by a factor of 2 (dashed green line). Ejecta volume (solid red line) and twice ejecta volume (dashed red line) for Z-EDOZ model with alternate crater size scaling factor CDCR = 1.4. Z = 2.7 in all cases. c) Ratio of power law fits in (a). Note fit to Utopia overlaps with fit to Acidalia. Utopia basin = circles; Isidis basin = squares; Acidalia Planitia = crosses; Lunae Planum = triangles; Solis Planum = diamonds. Fig. 22. Z model predictions compared to crater geometry measurements for (a) VAbove = VUplift + VEjecta, (b) VEjecta/VAbove, and (c) rim height-rim diameter ratio. Original Z model (solid black line), Z-EDOZ model (solid green line), and Z-EDOZ model with ejecta volume increased by a factor of 2 (dashed green line). Z-EDOZ model (solid red line) and Z-EDOZ model with twice ejecta volume (dashed red line) using alternate crater size scaling factor CDCR = 1.4. Z = 2.7 in all cases. Utopia basin = circles; Isidis basin = squares; Acidalia Planitia = crosses; Lunae Planum = triangles; Solis Planum = diamonds. profile requires some scrutiny. In the Z-EDOZ model, the height of the uplifted surface is about equal to the thickness of the ejecta at the final crater rim, in good agreement with the uplift profile coefficient of 0.5 used to analyze the data. Using a paraboloid transient cavity and Z = 2.7, the model slope of the uplifted surface is a −3 power law. If the uplifted surface occupies half the observed rim height, then a −3 power law lies above the observed crater topographic profiles near the crater rim. Note the steep exterior rim slopes in Fig. 11 and the “moat-like” depression around the crater rim in the example Utopian crater profile. Some of the uplifted surface may have slumped, moving outward from the crater rim, or been eroded as part of the ejecta emplacement process (e.g., Schultz and Gault 1979). Alternatively, the uplift profile predicted by the Z-EDOZ model is too high for Martian craters. The latter hypothesis is particularly problematic as the volume of the continuous ejecta blanket must increase as a result to match the measurements of VAbove. A −5.5 power law to analyze the data was used to maintain positive ejecta throughout the ejecta blanket, but this may have led to an overestimate of ejecta volume. None of the fresh crater profiles in Fig. 11 display simple ejecta profiles similar to the pure ballistic emplacement shown in Fig. 20. The steep nearrim topography may be a result of incorporation of the uplifted surface into the horizontal radial flow of the ejecta blanket, and the measured VEjecta includes a component of uplifted material. The Z-EDOZ model also predicts the height of the final crater rim, assuming ballistic ejecta and no postcrater formation relaxation. The sum of the model uplifted surface and ejecta thickness is shown with the crater rim height measurements, normalized to the final crater rim diameter in Fig. 22c. Again the original Z and Z-EDOZ models fall short of the data, and increasing the Z-EDOZ model rim ejecta thickness by a factor of two reaches the lower range of the observations. Therefore, we find that the original Z and Z-EDOZ models developed for terrestrial impact and explosion craters Ejecta and subsurface properties from crater geometries do not provide a good fit to fresh Martian impact craters. Note that the parameters controlling the distribution of ejecta are: Z, dZ ~ DP, the ejecta velocity scaling exponent m, and the scaling from the transient to final crater diameter CDCR. The values of Z and m are well constrained to lie between 2.5 and 3.0 and the model predictions are similar over this range. The diameter of the projectile DP is derived from the transient crater diameter using Equation 2 using the mean impact velocity on Mars. Varying the center of the flow field over 0.5DP < dZ < 1.5DP produces very similar results as presented in Figs. 21 and 22. On the other hand, the results of the Z and Z-EDOZ models are very sensitive to the scaling constants CDCR, which relate the transient crater surface diameter to the final crater rim diameter. As can be seen in Fig. 20, the amount of ballistic ejecta in the continuous ejecta blanket depends significantly on the final rim radius. In the standard scaling laws presented in the “Crater Scaling Relationships” section, CDCR = 1.2 × 1.3 = 1.56. Formation of transient crater cavities has been extensively studied in laboratory experiments and numerical simulations (e.g., O’Keefe and Ahrens 1993; Schmidt and Housen 1987). However, the final stage of crater formation, the collapse of the transient cavity to the final crater cavity, is not well understood (Melosh and Ivanov 1999). Using the crater geometry measurements, we can constrain the scaling from transient to final crater diameters on Mars. Recall that CD is derived for the craters in this study by fitting crater cavity volume and the transient cavity volume collapse coefficient (Equation 9). The fitted values of CD typically lie between 1.1 and 1.2. The ratio between the rim diameter and surface diameter for simple craters, CR, is measured to lie between 1.15 and 1.2 (Fig. 17). If the observed DR/DS ratio has not been significantly modified, the scaling coefficient CDCR is about 1.4 for the freshest Martian impact craters. The predictions of the Z-EDOZ model with CDCR = 1.4 are shown as red lines in Figs. 21 and 22. Using the new scaling factor, the Z-EDOZ model produces a much better fit to the observations. The calculated range for VAbove are in good agreement with the data. Note that the rim height predictions around simple craters are best fit by negligible bulking of the ejecta, which is expected around the rims of small craters. On the other hand, the rim heights of larger craters are better fit with significant ejecta bulking, which is a result of the increased fragmentation of ejecta landing near the final rim of complex craters. However, the predicted value of VEjecta is less than measured, even when bulking is taken into account. The steep near rim topography indicates that, even with the new crater size scaling factor, some of the uplifted material must have been incorporated in the near-rim ejecta flow. Based on comparison to the rim heights, ejecta volumes, and uplifted volumes on both highland and lowland terrains, a crater size scaling factor of CDCR = 1.4 and the Z-EDOZ 1533 model with dz ~ 1DP and Z = 2.7 may be used to predict the volumes of ejecta and uplifted material around fresh Martian craters. The refined model represents the first observationally constrained set of scaling laws for craters on Mars. The model provides a reliable set of initial parameters for studies of the dynamics of the ejecta blanket emplacement processes and the origin of fluidized ejecta morphologies on Mars. Alternatively, the poor fit of the Z-EDOZ model with CDCR = 1.56 could suggest that Martian ejecta blankets are not well characterized by ballistic emplacement. In order to match the observed geometric properties of Martian craters, the volume of material in the continuous ejecta blanket must draw on both the distal ejecta with ballistic range >3RR and proximal ejecta that land between the transient and final crater diameters. The observed continuous ejecta blanket has a volume almost equal to the entire predicted excavation volume. The high ejection velocities near the impact point result in large ballistic range. For example, the ballistic ejecta with a range between 3RR and 6RR have a volume of only 1/3VBE ~ 1/9VEXC. It is unlikely that the atmosphere can collect a significant fraction of the high velocity distal ejecta. Ejecta that lands near the transient crater rim may also contribute to the continuous ejecta blanket by horizontal motion away from the transient crater rim before the collapse to the final crater diameter. Again, it is unlikely that a significant fraction of ejecta landing within the final crater diameter contributes to the continuous ejecta blanket because the ejection and ballistic flight times are comparable to the crater collapse time (Melosh and Ivanov 1999), and there would be little time for significant horizontal flow before crater collapse. Finally, the resistance of the freshest ejecta blankets to eolian erosion suggests that the bulk ejecta cannot be too porous and friable. Therefore, volume of the continuous ejecta blanket cannot be explained solely by concentration of the excavated ejecta, and the standard size scaling laws require some refinement to match the observations. CONCLUSIONS Using the Mars Orbiter Laser Altimeter data set, we have measured the geometric properties of fresh craters in highland and lowland plains to investigate the emplacement processes of fluidized ejecta blankets and regional differences in the strength of the upper crust and crater collapse processes. Crater geometries provide unique insights into the properties of the Martian subsurface, and quantitative differences between terrains may be mapped using the geometry of fresh Martian craters. The major findings of this study are: 1. Using simulated craters on Martian background terrains, we demonstrate that the MOLA data has sufficient resolution to measure fresh crater geometries with high precision and accuracy on all craters with DR ≥ 6 km, and 1534 S. T. Stewart and G. J. Valiant some craters as small as DR = 2 km. Robust testing of the HMars crater measurement toolkit shows no systematic errors. Measurement errors are dominated by uneven background terrains. In this study of 3–50 km rim diameter fresh craters on highland and lowland plains, we resolve significant regional differences in crater depth and cavity volumes. Simple lowland craters are 1.5–2.0 times deeper (≥5σo difference) with >50% larger cavities (≥2σo) compared to highland craters of the same diameter. Rim heights and the volume of material above the preimpact surface are slightly larger in the lowlands over most of the size range studied. This work presents quantitative analyses of crater geometries to infer regional differences in subsurface properties. Using crater scaling laws for the volume of the transient crater cavity, we conclude that the freshest lowland craters in this study preserve a larger volume fraction (45–70%) of the transient cavity compared to highland craters (25–50%). Therefore, the effective strength of the upper crust in the lowland plains is stronger than in the highland plains. The presence of large simple craters up to 10.5 km in rim diameter in Utopia indicate that the effective yield strength of the upper several kilometers is similar to competent rock (approximately 9–12 MPa). In Isidis and Acidalia, similar effective strength is inferred from similar depthto-diameter ratios and cavity volumes. Unlike lowland craters, the geometries of fresh highland simple craters grade smoothly to large, gravity-dominated crater sizes. Based on previous studies using Viking data (Pike 1988), the transition from strength- to gravity-dominated craters in the highlands may be as small as about 3–5 km, implying effective strengths of about 3–6 MPa, in between the range of values for soils and competent rock. The large simple craters in Utopia reflect the strength of the Noachian infill in the Utopia basin. The weaker highland materials are consistent with weak, highly fractured, layered, and possibly volatile-rich Noachian deposits. Standard crater scaling laws and the Z-EDOZ crater excavation model underpredict the volume of the continuous ejecta blanket, volume of uplifted material, and the rim height by a factor of about 3. Bulking and concentration of ejected materials by nonballistic processes cannot explain the discrepancy between the model and observations. Using the geometry measurements of fresh Martian craters and analysis of crater collapse, we revise the scaling factor between the transient crater surface diameter and final crater rim diameter to CDCR = 1.4. Comparisons between the predicted volumes of ejecta and uplifted materials indicate significant deviations from ballistic emplacement, including evidence for postemplacement horizontal flow and incorporation of surface materials into the ejecta flow. 6. The crater geometry measurements are well-fit by the revised crater size scaling factor of CDCR = 1.4 and the Z-EDOZ model with Z = 2.7 and the center of flow at a depth of about 1 projectile diameter (dZ ~ 1DP). This crater scaling and excavation model may be used to predict the volumes of ejecta and uplifted material around fresh Martian craters. The refined model represents the first observationally constrained set of scaling laws for Mars. The model provides a reliable set of initial parameters for studies of the dynamics of the ejecta blanket emplacement processes and the origin of fluidized ejecta morphologies on Mars. Acknowledgments–This work was supported by NASA grant NAG5-13474. The authors thank Vaidya G. Rajagopalan, Zain Khalid, and Ben Black for their assistance with crater measurements. We appreciate the thoughtful reviews by J. Boyce and J. Garvin, which improved this work. Editorial Handling—Dr. Nadine Barlow REFERENCES Ahrens T. J. and O’Keefe J. D. 1987. Impact on the Earth, ocean and atmosphere. International Journal of Impact Engineering 5:13– 32. Anderson J. L. B., Schultz P. H., and Heineck J. T. 2004. Experimental ejection angles for oblique impacts: Implications for the subsurface flow-field. Meteoritics & Planetary Science 39:303–320. Barlow N., Stewart S., and Barnouin-Jha O. 2005. The role of volatiles and atmospheres on Martian impact craters (abstract). Eos 86:433. Barlow N. G., Boyce J. M., Costard F. M., Craddock R. A., Garvin J. B., Sakimoto S. E. H., Kuzmin R. O., Roddy D. J., and Soderblom L. A. 2000. Standardizing the nomenclature of Martian impact crater ejecta morphologies. Journal of Geophysical Research 105:26,733–26,738. Barlow N. G. and Bradley T. L. 1990. Martian impact craters: Correlations of the ejecta and interior morphologies with diameter, altitude, and terrain. Icarus 87:156–179. Barnouin-Jha O. S., Baloga S., and Glaze L. 2005. Comparing landslides to fluidized crater ejecta on Mars. Journal of Geophysical Research, doi:10.1029/2003JE002214. Barnouin-Jha O. S. and Schultz P. H. 1998. Lobateness of impact ejecta deposits from atmospheric interactions. Journal of Geophysical Research 103:25,739–25.756. Black B. A. and Stewart S. T. Forthcoming. Impact crater geometries provide evidence for ice-rich layers at low latitudes on Mars. Journal of Geophysical Research. Bottke W. F., Nolan M. C., Greenberg R., and Kolvoord R. A. 1994. Collisional lifetimes and impact statistics of near-Earth asteroids. In Hazards due to comets and asteroids, edited by Gehrels T. Tucson, Arizona: The University of Arizona Press. pp. 337–357. Boyce J. M., Mouginis-Mark P., and Garbeil H. 2005a. Ancient oceans in the northern lowlands of Mars: Evidence from impact crater depth/diameter relationships. Journal of Geophysical Research 110, doi:10.1029/2004JE002328. 2. 3. 4. 5. Ejecta and subsurface properties from crater geometries 1535 Boyce J. M., Mouginis-Mark P., Garbeil H., and Tornabene L. L. 2006. Deep impact crater in the Isidis and southwestern Utopia Planitia regions of Mars: High target material strength as a possible cause. Geophysical Research Letters 33, doi:10.1029/ 2005GL024462. Boyce J. M., Mouginis-Mark P. J., Garbeil H., and Soderblom L. A. 2005b. History of major degradational events in the highlands of Mars: Preliminary results from crater depth/diameter measurements (abstract #1055). 36th Lunar and Planetary Science Conference. CD-ROM. Buczkowski D. L., Frey H. V., Roark J. H., and McGill G. E. 2005. Buried impact craters: A topographic analysis of quasi-circular depressions, Utopia basin, Mars. Journal of Geophysical Research 110, doi:10.1029/2004JE002324. Carr M. H., Crumpler L. S., Cutts J. A., Greeley R., Guest J. E., and Masursky H. 1977. Martian impact craters and emplacement of the ejecta by surface flow. Journal of Geophysical Research 82: 4055–4065. Cintala M. J. and Mouginis-Mark P. J. 1980. Martian fresh crater depths: More evidence for subsurface volatiles? Geophysical Research Letters 7:329–332. Croft S. K. 1980. Cratering flow fields: Implications for the excavation and transient expansion stages of crater formation. Proceedings, 11th Lunar and Planetary Science Conference. pp. 2347–2378. Forsberg-Taylor N. K., Howard A. D., and Craddock R. A. 2004. Crater degradation in the Martian highlands: Morphometric analysis of the Sinus Sabaeus region and simulation modeling suggest fluvial processes. Journal of Geophysical Research 109, doi:10.1029/2004JE002242. Frey H., Hohner G., Wernecke A., Roark J., and Sakimoto S. 2002a. Constraining the thickness of ridged plains using relief from buried impact basins compared to exposed craters on adjacent cratered terrain (abstract #P32A-02). 2002 American Geophysical Union Spring Meeting. CD-ROM. Frey H. V., Doudnikoff C. E., and Mongeon A. M. 1991. Are Noachian-age ridged plains (Nplr) actually Early Hesperian in age. Proceedings, 21st Lunar and Planetary Science Conference. pp. 635–644. Frey H. V., Roark J. H., Shockey K. M., Frey E. L., and Sakimoto S. E. H. 2002b. Ancient lowlands on Mars. Geophysical Research Letters 29, doi:0.1029/2001GL013832. Frey H. V., Sakimoto S. E. H., and Roark J. H. 1998. MOLA topography and the Isidis basin: Constraints on basin center and ring diameters (abstract #1631). 29th Lunar and Planetary Science Conference. CD-ROM. Fudali R. F., Milton D. J., Fredriksson K., and Dube A. 1980. Morphology of Lonar crater, India: Comparisons and implications. The Moon and the Planets 23:493–515. Garvin J. B. and Frawley J. J. 1998. Geometric properties of Martian impact craters: Preliminary results from the Mars Orbiter Laser Altimeter. Geophysical Research Letters 25:4405–4408. Garvin J. B., Sakimoto S. E. H., and Frawley J. J. 2003. Craters on Mars: Global geometric properties from gridded MOLA topography (abstract #3277). The Sixth International Conference on Mars. CD-ROM. Garvin J. B., Sakimoto S. E. H., Frawley J. J., and Schnetzler C. 2000. North polar region craterforms on Mars: Geometric characteristics from the Mars Orbiter Laser Altimeter. Icarus 144:329–352. Garvin J. B., Sakimoto S. E. H., Frawley J. J., and Schnetzler C. 2002. Global geometric properties of Martian impact craters (abstract #1255). 33rd Lunar and Planetary Science Conference. CD-ROM. Gault D. E. and Greeley R. 1978. Exploratory experiments of impact craters formed in viscous-liquid targets: Analogs for Martian rampart craters? Icarus 34:486–495. Grieve R. A. F. and Garvin J. B. 1984. A geometric model for excavation and modification at terrestrial simple impact craters. Journal of Geophysical Research 89:11,561–11,572. Hamilton V. E. and Christensen P. R. 2005. Evidence for extensive, olivine-rich bedrock on Mars. Geology 33:433–436. Hamilton V. E., Christensen P. R., McSween H. Y., Jr., and Bandfield J. L. 2003. Searching for the source regions of Martian meteorites using MGS TES: Integrating Martian meteorites into the global distribution of igneous materials on Mars. Meteoritics & Planetary Science 38:871–885. Head J. N., Melosh H. J., and Ivanov B. A. 2002. Martian meteorite launch: High-speed ejecta from small craters. Science 298:1752– 1756. Hoefen T. M., Clark R. N., Bandfield J. L., Smith M. D., Pearl J. C., and Christensen P. R. 2003. Discovery of olivine in the Nili Fossae region of Mars. Science 302:627–630. Holsapple K. A. 1993. The scaling of impact processes in planetary sciences. Annual Review of Earth and Planetary Sciences 21: 333–373. Holsapple K. A. and Schmidt R. M. 1982. On the scaling of crater dimensions 2. Impact processes. Journal of Geophysical Research 87:1849–1870. Holsapple K. A. and Schmidt R. M. 1987. Point-source solutions and coupling parameters in cratering mechanics. Journal of Geophysical Research 92:6350–6376. Horner V. M. and Greeley R. 1982. Pedestal craters on Ganymede. Icarus 51:549–562. Ivanov B. A. 1996. Spread of ejecta from impact craters and the possibility of estimating the volatile content of the Martian crust. Solar System Research 30:43–58. Kuz’min R. O., Bobina N. N., Zabalueva E. V., and Shashkina V. P. 1988. Structural inhomogeneities of the Martian cryolithosphere. Solar System Research 22:121–133. Lee D. T. and Schachter B. J. 1980. Two algorithms for constructing a Delaunay triangulation. International Journal of Computer and Information Science 9:219–242. Malin M. C. and Edgett K. S. 2000. Sedimentary rocks of early Mars. Science 290:1927–1937. Maxwell D. E. 1973. Cratering flow and crater prediction methods. Technical Memorandum #TCAM 73–17. San Leandro, California: Physics International. 50 p. Maxwell D. E. 1977 Simple Z model of cratering, ejection, and the overturned flap. In Impact and explosion cratering, edited by Roddy D. J., Pepin R. O., and Merrill R. B. New York: Elmsford. pp. 983–1022. Melosh H. J. 1989. Impact cratering: A geologic process. New York: Oxford University Press. 253 p. Melosh H. J. 1998. Crater scaling program. http:// www.lpl.arizona.edu/tekton/crater.html. Melosh H. J. and Ivanov B. A. 1999. Impact crater collapse. Annual Review of Earth and Planetary Sciences 27:385–415. Meresse S., Costard F., and Mangold N. 2005. Evidence for removal episode on northern plains from the Martian fluidized ejecta volume (abstract #1748). 36th Lunar and Planetary Science Conference. CD-ROM. Mouginis-Mark P. 1981. Ejecta emplacement and modes of formation of Martian fluidized ejecta craters. Icarus 45:60–76. Mouginis-Mark P. J. 1978. Morphology of Martian rampart craters. Nature 272:691–694. O’Keefe J. D. and Ahrens T. J. 1993. Planetary cratering mechanics. Journal of Geophysical Research 98:17,011–17,028. 1536 S. T. Stewart and G. J. Valiant Oberbeck V. R. 1975. The role of ballistic erosion and sedimentation in lunar stratigraphy. Reviews of Geophysics and Space Physics 13:337–362. Piekutowski A. J. 1977. Cratering mechanisms observed in laboratory-scale high-explosive experiments. In Impact and explosion cratering, edited by Roddy D. J., Pepin R. O., and Merrill R. B. New York: Pergamon Press. pp. 67–102. Pierazzo E. and Melosh H. J. 2000. Melt production in oblique impacts. Icarus 145:252–261. Pierazzo E., Vickery A. M., and Melosh H. J. 1997. A reevaluation of impact melt production. Icarus 127:408–423. Pike R. J. 1977. Size-dependence in the shape of fresh impact craters on the Moon. In Impact and explosion cratering, edited by Roddy D. J., Pepin R. O., and Merrill R. B. New York: Pergamon Press. pp. 489–509. Pike R. J. 1980. Control of crater morphology by gravity and target type: Mars, Earth, Moon. Proceedings, 11th Lunar and Planetary Science Conference. pp. 2159–2189. Pike R. J. 1988. Geomorphology of impact craters on Mercury. In Mercury, edited by Vilas F., Chapman C. R., and Matthews M. S. Tucson, Arizona: The University of Arizona Press. pp. 165–273. Quaide W. L. and Oberbeck V. R. 1968. Thickness determinations of lunar surface layer from lunar impact craters. Journal of Geophysical Research 73:5247–5270. Rodriguez J. A. P., Sasaki S., Dohm J. M., Tanaka K. L., Strom B., Kargel J., Kuzmin R., Miyamoto H., Spray J. G., Fairen A. G., Komatsu G., Kurita K., and Baker V. 2005. Control of impact crater fracture systems on subsurface hydrology, ground subsidence, and collapse, Mars. Journal of Geophysical Research 110, doi:10.1029/2004JE002365 Schmidt R. M. and Housen K. R. 1987. Some recent advances in the scaling of impact and explosion cratering. International Journal of Impact Engineering 5:543–560. Schultz P. H. 1992. Atmospheric effects on ejecta emplacement. Journal of Geophysical Research 97:11,623–11,662. Schultz P. H. and Gault D. E. 1979. Atmospheric effects on Martian ejecta emplacement. Journal of Geophysical Research 84:7669– 7687. Scott D. H. and Tanaka K. L. 1986. Geologic map of the western equatorial region of Mars. U.S. Geological Survey Miscellaneous Investigations Series Map I-1802-A. Scale 1: 15,000,000. Searls M. L. and Phillips R. J. 2004. Mechanics of Utopia basin on Mars (abstract #1822). 35th Lunar and Planetary Science Conference. CD-ROM. Smith D. E., Zuber M. T., Frey H. V., Garvin J. B., Head J. W., Muhleman D. O., Pettengill G. H., Phillips R. J., Solomon S. C., Zwally H. J., Banerdt W. B., Duxbury T. C., Golombek M. P., Lemoine F. G., Neumann G. A., Rowlands D. D., Aharonson O., Ford P. G., Ivanov A. B., Johnson C. L., McGovern P. J., Abshire J. B., Afzal R. S. and Sun X. L. 2001. Mars Orbiter Laser Altimeter: Experiment summary after the first year of global mapping of Mars. Journal of Geophysical Research 106:23,689– 23,722. Smith D. E., Zuber M. T., Solomon S. C., Phillips R. J., Head J. W., Garvin J. B., Banerdt W. B., Muhleman D. O., Pettengill G. H., Neumann G. A., Lemoine F. G., Abshire J. B., Aharonson O., Brown C. D., Hauck S. A., Ivanov A. B., McGovern P. J., Zwally H. J., and Duxbury T. C. 1999. The global topography of Mars and implications for surface evolution. Science 284:1495–1503. Stewart S. T., Maloof A. C., Newsom H. E., Wright S. P., Louzada K. L., and Weiss B. P. 2005. Field observations of groundhugging ejecta flow at Lonar crater, India (abstract #3045). Workshop on the Role of Volatiles and Atmospheres on Martian Impact Craters. CD-ROM. Stewart S. T., O’Keefe J. D., and Ahrens T. J. 2001. The relationship between rampart crater morphologies and the amount of subsurface ice (abstract #2092). 32nd Lunar and Planetary Science Conference. Stewart S. T., O’Keefe J. D., and Ahrens T. J. 2004. Impact processing and redistribution of near-surface water on Mars. In Shock compression of condensed matter 2003, edited by Furnish M. D., Gupta Y. M., and Forbes J. W. Melville, New York: American Institute of Physics. pp. 1484–1487. Thomsen J. M., Austin M. G., Ruhl S. F., Schultz P. H., and Orphal D. L. 1979. Calculational investigation of impact cratering dynamics: Early time material motions. Proceedings, 10th Lunar and Planetary Science Conference. pp. 2741–2756. Thomsen J. M., Austin M. G., and Schultz P. H. 1980. The development of the ejecta plume in a laboratory-scale impact cratering event (abstract). Proceedings, 11th Lunar and Planetary Science Conference. pp. 1146–1148. Thomson B. J. and Head J. W. 1999. Utopia basin, Mars: A new assessment using Mars Orbiter Laser Altimeter (MOLA) data (abstract #1894). 30th Lunar and Planetary Science Conference. CD-ROM. Thomson B. J. and Head J. W. 2001. Utopia basin, Mars: Characterization of topography and morphology and assessment of the origin and evolution of basin internal structure. Journal of Geophysical Research 106:23,209–23,230. Valiant G. J. and Stewart S. T. 2004. Martian surface properties: Inferences from resolved differences in crater geometries (abstract #1293). 35th Lunar and Planetary Science Conference. CD-ROM. Appendix A. Variables. Variable c C2, C3 CD CR CV DP dR DR, DS dS dST DST Description Cutoff parameter for fresh crater definition π2, π3 scaling constants Ratio between DST and DS (simple craters) Ratio between DR and DS (simple craters) Crater volume collapse coefficient (simple craters) Projectile diameter Depth from rim Rim diameter, surface diameter Surface depth Transient crater surface depth Transient crater surface diameter Ejecta and subsurface properties from crater geometries Appendix A. Continued. Variables. Variable dZ D* Description Depth to flow center in Z-EDOZ model Transition rim diameter between strength and gravity regimes Observed maximum crater surface depth, as function of diameter and region Observed maximum rim height, as function of diameter and region Surface gravity Height of uplifted preimpact surface Rim height Ejection velocity scaling exponent Projectile mass Number of rim points from MOLA tracks Radial distance from crater center line; polar coordinate for Z-EDOZ model Distance between flow center, dZ, and uplift location (ru, hu) Rim radius, surface radius Surface radius of transient cavity Surface radius at uplift elevation hu Combined volume of ejecta and uplift Volume of ballistic ejecta with range between 1RR and 3RR Cavity volume below preimpact surface Ejection velocity for Z-EDOZ model stream tube Volume of ejecta Total excavated volume Projectile impact velocity Cavity volume below rim Ejection velocity for stream tube at transient crater radius Volume of transient crater cavity below preimpact surface Volume of uplifted material between 1RR and 3RR Combined uplift and ejecta volume between 1RR and 3RR predicted by Z-EDOZ model Excavated volume by a Z-EDOZ trajectory emerging at radius rS Effective yield strength of surface Z and Z-EDOZ model parameter π2 scaling exponent Angle from flow center to uplift point in Z-EDOZ model Angle between the horizontal at depth dZ and intersection of stream tube at surface Polar coordinate in Z-EDOZ model Scaling exponent coupling parameter Inverse Froude number Inverse Cauchy number Nondimensional crater diameter Nondimensional crater volume Density π3 scaling exponent Systematic measurement error Observed measurement scatter around a given fit Impact angle from surface Ejection angle from surface 1537 g hu HR m mP NR r Rh RR, RS RST ru VAbove VBE VCavity ve VEjecta VEXC vi VRim vST VT VUp VZ−Above Y Z β δ ∆ θ µ π2 π3 πD πV ρ σ σm σo φ ϕ Fd S FH R V Z,r S