Can Observations of Earthquake Scaling Constrain Slip Weakening?

Rachel E. Abercrombie, & James R. Rice

Published August 2005, SCEC Contribution #728

We use observations of earthquake source parameters over a wide magnitude range (MW~ 0–7) to place constraints on constitutive fault weakening. The data suggest a scale dependence of apparent stress and stress drop; both may increase slightly with earthquake size. We show that this scale dependence need not imply any difference in fault zone properties for different sized earthquakes. We select 30 earthquakes well-recorded at 2.5 km depth at Cajon Pass, California. We use individual and empirical Green's function spectral analysis to improve the resolution of source parameters, including static stress drop (Δσ) and total slip (S). We also measure radiated energy ES. We compare the Cajon Pass results with those from larger California earthquakes including aftershocks of the 1994 Northridge earthquake and confirm the results of Abercrombie (1995): μES/M0<<Δσ (where μ= rigidity) and both ES/M0 and Δσ increase as M0 (and S) increases. Uncertainties remain large due to model assumptions and variations between possible models, and earthquake scale independence is possible within the resolution. Assuming that the average trends are real, we define a quantity G'= (Δσ− 2μES/M0)S/2 which is the total energy dissipation in friction and fracture minus σ1S, where σ1 is the final static stress. If σ1=σd, the dynamic shear strength during the last increments of seismic slip, then G'=G, the fracture energy in a slip-weakening interpretation of dissipation. We find that G' increases with S, from ~103 J m−2 at S= 1 mm (M1 earthquakes) to 106–107 J m−2 at S= 1 m (M6). We tentatively interpret these results within slip-weakening theory, assuming G'≈G. We consider the common assumption of a linear decrease of strength from the yield stress (σp) with slip (s), up to a slip Dc. In this case, if either Dc, or more generally (σp−σd) Dc, increases with the final slip S we can match the observations, but this implies the unlikely result that the early weakening behaviour of the fault depends on the ultimate slip that the fault will sustain. We also find that a single slip-weakening function σF(s) is able to match the observations, requiring no such correlation. Fitting G' over S= 0.2 mm to 0.2 m with G'S1+n, we find n~ 0.3, implying a strength drop from peak σp−σF(S) Sn. This model also implies that slip weakening continues beyond the final slip S of typical earthquakes smaller than ~M6, and that the total strength drop σp−σd for large earthquakes is typically >20 MPa, larger than Δσ. The latter suggests that on average a fault is initially stressed below the peak strength, requiring stress concentration at the rupture front to propagate slipping.

Key Words
earthquake source parameters, fracture energy, seismic energy, slip-weakening, source scaling, stress drop

Citation
Abercrombie, R. E., & Rice, J. R. (2005). Can Observations of Earthquake Scaling Constrain Slip Weakening?. Geophysical Journal International, 162(2), 406-424. doi: 10.1111/j.1365-246X.2005.02579.x.