Academia.eduAcademia.edu

Displacement-Length Scaling of Faults on Earth, Mars, and Beyond

2004

Faults on Mars exhibit ~1/5 of the maximum displacement D max of terrestrial faults for any given length L. This systematic shift to smaller D max /L ratios results from the gravity dependence of shear driving stress, yield strength of rock at the fault tipline, and stiffness (modulus) of the surroundings. To first order, the D max /L ratio scales with gravity, so faults of any type on smaller planets and satellites will accumulate smaller displacements per unit length than faults on Earth. This reduction is modulated by the lithology and pore-pressure state of the lithosphere.

Lunar and Planetary Science XXXV (2004) 1157.pdf DISPLACEMENT-LENGTH SCALING OF FAULTS ON EARTH, MARS, AND BEYOND. Richard A. Schultz1, Chris H. Okubo1, and Scott J. Wilkins1,2 1Geomechanics–Rock Fracture Group, Department of Geological Sciences/ 172, Mackay School of Earth Sciences and Engineering, University of Nevada, Reno, NV 89557–0138 (www.mines.unr.edu/geo-eng/geomech; [email protected]); 2Now with Shell International Exploration and Production, P.O. Box 481, Houston, Texas 77001. Background: For an individual fault (having a central well-slipped portion bounded by frictionally stronger end zones: a “Dugdale-Barenblatt” model) [8], Dmax 2 1±ν 2 σ = C σ d ± σ y 1 ± cos π σ d L E 2 y (1) in which ν is Poisson’s ratio, E is Young’s modulus (both of the surrounding rock), L is fault length, σd is the cumulative shear driving stress on the fault, and σy is the yield strength of rock bounding the fault. Plausible values of σy /σd are 2–3 [e.g., 7,8] so that C = 0.4–0.6. 10 4 10 2 10 -1 10 γ basalt (wet) basalt (dry) SS (wet) SS (dry) -3 Mars Mars Mars Mars 10 10 3 -2 Earth basalt (wet) Earth SS (wet) = Introduction: Precision measurements of the maximum displacement (“offset,” Dmax) and map lengths L of surface-breaking faults on Mars [1,2] and Mercury [3,4] demonstrate that less displacement per unit length is accumulated along faults on these planets than along terrestrial ones. For example, normal faults from Tempe Terra (Mars) and thrust faults from Arabia (Mars) show Dmax/L ratios of 6.7 x 10–3 [2] and 6 x 10–3 [3], respectively. Thrust faults from Mercury also show Dmax/L ratios of 6.5 x 10–3 [4]. Typical values for terrestrial faults (normal, strike-slip, or thrust) are ~1–5 x 10–2 [5,6]. Data from the literature for normal faults from Earth and Mars are shown in Fig. 1. In this abstract we demonstrate the key role played by gravity in displacement-length (D-L) scaling of faults and, by implication, other types of fractures such as joints and dikes, on other planetary bodies. We incorporate gravity explicitly into the standard D-L scaling relations for faults [7,8] and find that the systematic shift toward smaller maximum displacements for normal and thrust faults on Mars and Mercury is related to the reduced gravity on these planets relative to the Earth. Equation (1) shows that the Dmax/L ratio for faults depends on three key factors [7,8]: ♦ Driving Stress. This is the shear stress leading to slip and displacement along the fault, represented as σd = (∆µ σn), where ∆µ is the difference between static and residual (or maximum and dynamic) friction coefficients for single-slip events [7,9], and σn is the magnitude of effective (compressive) normal stress resolved from far-field tectonic stresses onto the fault plane. For faults, the relevant quantity is the cumulative shear driving stress [7,10] associated with the total (cumulative) geologic offset along the fault. ♦ Rock Strength. The shear yield strength of unfaulted rock at the fault’s tipline modulates fault propagation and hence, the Dmax/L ratio. Stronger rock requires greater near-tip stresses to break, leading to larger values of displacement along the fault [7,8,11]. ♦ Rock Stiffness. The stiffness of surrounding unfaulted rock is given by S ∝ (1–ν 2)/E. The various moduli (Young’s, shear, and deformation [12,13]) are interchangable here. As modulus increases, stiffness decreases and displacement along a fault decreases. Dmax, m Summary: Faults on Mars exhibit ~1/5 of the maximum displacement Dmax of terrestrial faults for any given length L. This systematic shift to smaller Dmax/L ratios results from the gravity dependence of shear driving stress, yield strength of rock at the fault tipline, and stiffness (modulus) of the surroundings. To first order, the Dmax/L ratio scales with gravity, so faults of any type on smaller planets and satellites will accumulate smaller displacements per unit length than faults on Earth. This reduction is modulated by the lithology and pore-pressure state of the lithosphere. 10 1 10 – 1 10 10 2 10 3 10 4 10 5 10 6 L, m Fig. 1. Displacement-length data from Earth (L < 10 km, small circles, normal faults, after [16]) and Mars (Tempe Terra normal faults, L > 10 km). Upper two curves calculated for terrestrial data; lower four curves calculated for Mars. Curve order (planet, rock type, pore-pressure state) given by caption. Lunar and Planetary Science XXXV (2004) 1157.pdf Results: We find that linear scaling, with n = 1 (giving a constant ratio of Dmax/L as on Fig. 1 [6,7]) requires constant values of shear driving stress, rock (yield) strength, and stiffness (modulus) across the length scale of the fault population. Our calculations assume a depth of 100 m over which these parameters were averaged; varying the depth interval does not alter the scaling conclusions. Planetary gravity enters into each of these three factors, but more subtly than simply as a ratio of planetary gravities (e.g., gMars/gEarth) because the reduction in Dmax/L exceeds the gravity ratio (Fig. 1). Driving (shear) stress is reduced on lower gravity planets because the resolved effective normal stress σn, and the horizontal and vertical far-field stresses (which promote fault slip), are all dependent on σv = ρgz, in which ρ is average (wet or dry) rock density and z is depth. Rock strength is calculated from the Hoek-Brown (Mohr) envelope for rock masses [12,13] assuming a Rock Mass Rating (RMR) of 50, appropriate for wet conditions and typical near-surface fracture densities; dry conditions are considered by increasing RMR to 65, following standard practice [e.g., 12]. For the basalt, ρ = 2900 kg m–3, mi = 22 and the unconfined compressive strength of the intact rock material σc = 250 MPa; for a weaker rock mass (such as tuff or poorly indurated sedimentary material), the parameters for sandstone are used (ρ = 2200 kg m–3, mi = 19, σc = 100 MPa). The yield strength is taken to be the peak shear strength of the rock mass averaged over the 100-m depth interval used. Rock strength increases with planetary gravity for the same Mars Earth SS Basalt D/L 0.01 Wet Dry Wet Sandstone Dry Basalt 0.001 1000 Σσd σy or σd 100 depth range below the surface and is analogous to the well-known dependence of peak strength and confining pressure observed in experiments. Young’s modulus is calculated from the known increase with depth from Earth [14], starting with a value at the surface corresponding to the deformation modulus calculated from RMR [12,13]. Because density and the pore-pressure state of the crust (that modulate E) depend on gravity, Young’s modulus decreases with gravity for equivalent conditions and depth ranges below the surface. Values of the main parameters for Mars, normalized by the values for the terrestrial wet basaltic rock mass, are shown in Fig. 2. For the same conditions of rock type (e.g., basaltic rock mass) and fluid-saturated crustal rocks (i.e., “wet” conditions), g reduces Dmax by gMars/gEarth = 0.38 (via the driving stress term). Yield strength in shear scales with gravity, with the Martian basaltic rock mass ~0.5 of the corresponding terrestrial one. Modulus decreases with decreasing g, to a normalized value of ~0.84 for the Martian case. The combined effect of g on all three key factors discussed above (0.38*0.5*0.84 = 0.16) is a reduction in Dmax/L of about a factor of 5–6, consistent with the data (Figs. 1 and 2). Scatter in D-L datasets arises from several sources [e.g., 5,6,8] including mechanical interaction and linkage, which are both found on Earth and Mars [1,15].Variations in rock type or pore-pressure state appear small in relation to these other factors. Conclusions and Implications: The systematically smaller values of displacement for faults on Mars and Mercury are real—not an artifact of measurement technique or data resolution. We infer that faults on Venus should accumulate somewhat smaller displacements than their terrestrial counterparts given a ~10% reduction in gravity (g = 8.8 m s–2) relative to the Earth. Faults on the icy satellites of Jupiter, Saturn, and beyond should also scale with gravity, with particular values of the Dmax/L ratio depending on appropriate values of near-tip ice strength and ice stiffness. 10 1 σy 0.1 σd σd* σy / σd 0.01 10 8 6 4 E, GPa 0 20 References: [1] Schultz, JGR 102, 12,009–12,015, 1997. [2] Wilkins et al., GRL 29, 1884, doi:10.1029/2002GL015391, 2002. [3] Watters et al., Geology 26, 991–994, 1998. [4] Watters et al., GRL 27, 3659–3662, 2000. [5] Cowie and Scholz, JSG 14, 1149–1156, 1992. [6] Clark and Cox, JSG 18, 147–152, 1996. [7] Cowie and Scholz, JSG 14, 1133– 1148, 1992. [8] Schultz and Fossen, JSG 24, 1389–1411, 2002. [9] Scholz, Nature 391, 37–42, 1998. [10] Schultz, GRL 30, 1593, doi:10.1029/2002GL016681, 2003. [11] Wibberley et al., Earth Planet. Sci. 331, 419–425, 2000. [12] Bieniawski, Engineering Rock Mass Classifications, Wiley, 1989. [13] Schultz, JSG 18, 1139–1149, 1996. [14] Rubin, Bull. Volcanol. 52, 302–319, 1990. [15] Schultz, JSG 21, 985– 993, 1999. [16] Schlische et al., Geology 24, 683–686, 1996. 40 Fig. 2. Parameters for D-L scaling from Mars normalized to Earth. [Abstract published in: Lunar and Planetary Science, XXXV, CD-ROM, Lunar and Planetary Institute, Houston (2004)].