An official website of the United States government
Abstract Many natural faults are believed to consist of velocity weakening (VW) patches surrounded by velocity strengthening (VS) sections. Numerical studies routinely employ this framework to study earthquake sequences including repeating earthquakes. In this laboratory study, we made a VW asperity, of lengthL, from a bare Poly(methyl methacrylate) PMMA frictional interface and coated the surrounding interface with Teflon to make VS fault sections. Behavior of this isolated asperity was studied as a function ofL(ranging from 100 to 400 mm) and the critical nucleation length, , which is inversely proportional to the applied normal stress (2–16 MPa). Consistent with recent numerical simulations, we observed aseismic slip for < 2, periodic slip for 2 < < 6, and non‐periodic slip for 10 < . Furthermore, we compared the experiments whereLwas contained by VS material to standard stick‐slip events whereLwas bounded by free surfaces (i.e.,L = the total sample length). The free surface case produced ∼10 times larger slip during stick‐slip events compared to the contained fault ruptures, even with identical . This disparity highlights how standard, complete‐rupture stick‐slip events differ from contained events expected in nature, due to both the free surface conditions and the heterogeneous normal stress along the fault near the free ends, as confirmed by Digital Image Correlation analysis. This study not only introduces the Teflon coating experimental technique for containing laboratory earthquake ruptures, but also highlights the utility of as a predictive parameter for earthquake behavior.
more »
« less
The effects of characteristic slip distance on earthquake nucleation styles in fully dynamic seismic cycles
Earthquake nucleation is a crucial preparation process of the following coseismic rupture propagation. Under the framework of rate-and-state friction, it was found that the ratios of a to b parameters control whether earthquakes nucleate as an expanding crack or a fixed length patch. However, as an essential parameter in earthquake physics, critical slip distance DRS controls the weakening efficiency of fault strength and can influence the nucleation styles. Here we investigate the effects of DRS on nucleation styles in the context of fully dynamic seismic cycles by evaluating the evolution of the nucleation zone quantitatively when it accelerates from the tectonic loading rate to seismic slip velocity. The inferred values of DRS from small-scale laboratory faults are 1-100 μm, several orders smaller than those obtained from geophysical observations on large natural faults. Considering the scale-dependence of widely observed DRS, the ratio of DRS to velocity weakening asperity size W is applied to substitute the absolute value of DRS in this study. We find when DRS/W is relatively large (~10-5), a/b=0.5 can separate two nucleation styles as found previously. For a relatively small DRS/W (~10-6), however, a/b larger than 0.7 is necessary to produce the typical expanding crack-like nucleation style. When DRS/W<4x10-7 and a/b<0.8, the fixed length nucleation style dominates. For some cases with a/b>0.75, the initial yielding phase accelerates to a considerable slip velocity just before the subsequent expanding fracture phase, which may explain the generation of foreshock activities. Specially, the first yielding phase is possible to trigger dynamic events without a secondary fracture phase. Furthermore, when the nucleation site is not in the middle of the asperity, large enough a/b (e.g., 0.8) could induce a complex nucleation style as well as abundant interseismic aseismic transients. We also recognize a special twin nucleation style that incorporates a failed acceleration phase. Our results reveal the critical role of DRS on earthquake nucleation styles and suggest that the fixed length nucleation style may be more common for the range of DRS/W (~10-4-~10-7) observed on natural and laboratory faults.
more »
« less
- Award ID(s):
- 1943742
- PAR ID:
- 10492003
- Publisher / Repository:
- 2023 AGU Annual Meeting
- Date Published:
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abundant heterogeneity has been documented on faults in nature across a wide range of length scales, including structural, mineralogical, and roughness variations. The role of complex heterogeneity on fault mechanics and frictional stability is not well established, and experiments investigating heterogeneity have typically incorporated a single source of heterogeneity. Here, we conduct rock friction experiments on rough, bimaterial faults that are creeping, or steadily sliding, to explore the role of lithological heterogeneity on fault mechanics and stability. When strong asperities juxtapose weak gouge, stable sliding occurs with a low friction coefficient, µ. Encounters of strong diabase asperities on talc gouge lined faults initiate dramatic increases in µ and transitions to unstable sliding characterized by frequent stick-slip events (StSE). Seismic moments and stress drops of StSE decrease with increasing asperity abundance. Stress is concentrated at asperities during encounters, increasing with decreasing asperity abundance and leading to extensive mechanical damage. Interactions between strong, velocity weakening asperities provide a model to explain the nucleation of seismic and aseismic slip events on nominally stable, creeping faults.more » « less
-
Abstract Understanding the dynamics of microearthquakes is a timely challenge with the potential to address current paradoxes in earthquake mechanics, and to better understand earthquake ruptures induced by fluid injection. We perform fully 3D dynamic rupture simulations caused by fluid injection on a target fault for Fault Activation and Earthquake Ruptures experiments generatingMw ≤ 1 earthquakes. We investigate the dynamics of rupture propagation with spatially variable stress drop caused by pore pressure changes and assuming different slip‐weakening constitutive parameters. We show that the spontaneous arrest of propagating ruptures is possible by assuming a high fault strength parameter S, that is, a high ratio between strength excess and dynamic stress drop. In faults with high S values (low rupturing potential), even minor variations inDc(from 0.45 to 0.6 mm) have a substantial effect on the rupture propagation and the ultimate earthquake size. Modest spatial variations of dynamic stress drop determine the rupture mode, distinguishing self‐arresting from run‐away ruptures. Our results suggest that several characteristics inferred for accelerating dynamic ruptures differ from those observed during rupture deceleration of a self‐arresting earthquake. During deceleration, a decrease of peak slip velocity is associated with a nearly constant cohesive zone size. Moreover, the residual slip velocity value (asymptotic value for a crack‐like rupture) decreases to nearly zero. This means that an initially crack‐like rupture becomes a pulse‐like rupture during spontaneous arrest. These findings highlight the complex dynamics of small induced earthquakes, which differ from solutions obtained from conventional crack‐like models of earthquake rupture.more » « less
-
Abstract In the quest to determine fault weakening processes that govern earthquake mechanics, it is common to infer the earthquake breakdown energy from seismological measurements. Breakdown energy is observed to scale with slip, which is often attributed to enhanced fault weakening with continued slip or at high slip rates, possibly caused by flash heating and thermal pressurization. However, seismologically inferred breakdown energy varies by more than six orders of magnitude and is frequently found to be negative-valued. This casts doubts about the common interpretation that breakdown energy is a proxy for the fracture energy, a material property which must be positive-valued and is generally observed to be relatively scale independent. Here, we present a dynamic model that demonstrates that breakdown energy scaling can occur despite constant fracture energy and does not require thermal pressurization or other enhanced weakening. Instead, earthquake breakdown energy scaling occurs simply due to scale-invariant stress drop overshoot, which may be affected more directly by the overall rupture mode – crack-like or pulse-like – rather than from a specific slip-weakening relationship.more » « less
-
Abstract Earthquake nucleation is a fundamental problem in earthquake science and has practical implications for forecasting seismic hazards. Laboratory experiments performed on large, meter‐scale fault systems offer unique insights into the nucleation process because the migration and expansion of the nucleation zone can be precisely detected, measured, and characterized using arrays of local strain and slip measurements. We report on a series of laboratory experiments conducted on a 1‐m direct shear machine. We sheared layers of quartz gouge between roughened acrylic forcing blocks over a range of normal stresses between 3 and 12 MPa, generating a spectrum of slip modes, ranging from aseismic creep to fast‐dynamic rupture. Co‐seismic slip, peak slip velocity, and high‐frequency acoustic energy content of laboratory earthquakes increases systematically with both cumulative fault slip and normal stress. Slower and smaller laboratory earthquake sequences have larger nucleation zones, creep more during their inter‐seismic period, and are deficient in high‐frequency energy compared to larger and faster rupture sequences. We find that the critical nucleation length scale,H*, scales inversely with cumulative fault slip and normal stress. A reduction inH*and an increase in event size can be explained by a decrease in the critical slip distance,Dc, or an increase in the frictional rate parameterb–aand is likely driven by shear localization. Together, our results indicate that homogeneous, mature fault zones that have undergone more cumulative fault slip are expected to have smallerH*and can more easily host dynamic instabilities, relative to immature faults.more » « less
