An official website of the United States government
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
The Effects of Characteristic Weakening Distance on Earthquake Nucleation Styles in Fully Dynamic Seismic Cycle Simulations
Abstract Earthquake nucleation is a crucial preparation process of the following coseismic rupture propagation. Under the framework of rate‐and‐state friction (RSF), it was found that the ratios of to parameters control whether earthquakes nucleate as an expanding crack or with a fixed length prior to the dynamic instability. However, the characteristic weakening distance controls the weakening efficiency of state variables in RSF and can influence the nucleation styles as well. Here we investigate the effects of 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. A larger (>0.75) is needed to produce expanding crack nucleation styles for relatively small , which suggests that fixed length nucleation styles may dominate on natural and laboratory faults. Furthermore, we find a more complex nucleation style when the nucleation site is not in the center of the asperity and identify a twin‐like nucleation style which includes two initial acceleration phases. We conclude that the earthquake nucleation style is strongly controlled by the value of . The possible dominance of fixed length nucleation styles suggests that the minimum size of earthquake rupture may be estimated at the early stage of the nucleation phase.
more »
« less
- Award ID(s):
- 1943742
- PAR ID:
- 10582169
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 129
- Issue:
- 12
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Determining conditions for earthquake slip on faults is a key goal of fault mechanics highly relevant to seismic hazard. Previous studies have demonstrated that enhanced dynamic weakening (EDW) can lead to dynamic rupture of faults with much lower shear stress than required for rupture nucleation. We study the stress conditions before earthquake ruptures of different sizes that spontaneously evolve in numerical simulations of earthquake sequences on rate‐and‐state faults with EDW due to thermal pressurization of pore fluids. We find that average shear stress right before dynamic rupture (aka shear prestress) systematically varies with the rupture size. The smallest ruptures have prestress comparable to the local shear stress required for nucleation. Larger ruptures weaken the fault more, propagate over increasingly under‐stressed areas due to dynamic stress concentration, and result in progressively lower average prestress over the entire rupture. The effect is more significant in fault models with more efficient EDW. We find that, as a result, fault models with more efficient weakening produce fewer small events and result in systematically lower b‐values of the frequency‐magnitude event distributions. The findings (a) illustrate that large earthquakes can occur on faults that appear not to be critically stressed compared to stresses required for slip nucleation; (b) highlight the importance of finite‐fault modeling in relating the local friction behavior determined in the lab to the field scale; and (c) suggest that paucity of small events or seismic quiescence may be the observational indication of mature faults that operate under low shear stress due to EDW.more » « less
-
Abstract This paper reviews laboratory observations of earthquake initiation and describes new experiments on a 3‐m rock sample where the nucleation process is imaged in detail. Many of the laboratory observations are consistent with previous work that showed a slow and smoothly accelerating earthquake nucleation process that expands to a critical nucleation length scaleLc, before it rapidly accelerates to dynamic fault rupture. The experiments also highlight complexities not currently considered by most theoretical and numerical models. This includes a loading rate dependency where a “kick” above steady state produces smaller and more abrupt initiation. Heterogeneity of fault strength also causes abrupt initiation when creep fronts coalesce on a stuck patch that is somewhat stronger than the surrounding fault. Taken together, these two mechanisms suggest a rate‐dependent “cascade up” model for earthquake initiation. This model simultaneously accounts for foreshocks that are a by‐product of a larger nucleation process and similarities between initialPwave signatures of small and large earthquakes. A diversity of nucleation conditions are expected in the Earth's crust, ranging from slip limited environments withLc< 1 m, to ignition‐limited environments withLc> 10 km. In the latter case,Lcfails to fully characterize the initiation process since earthquakes nucleate not because a slipping patch reaches a critical length but because fault slip rate exceeds a critical power density needed to ignite dynamic rupture.more » « less
-
Abstract Several regularly recurring moderate‐size earthquakes motivated dense instrumentation of the Parkfield section of the San Andreas fault (SAF), providing an invaluable near‐fault observatory. We present a seismo‐geodetic dynamic inversion of the 2004 Parkfield earthquake, which illuminates the interlinked complexity of faulting across time scales. Using fast‐velocity‐weakening rate‐and‐state friction, we jointly model coseismic dynamic rupture and the 90‐day evolution of postseismic slip in a 3D domain. We utilize a parallel tempering Markov chain Monte Carlo approach to solve this non‐linear high‐dimensional inverse problem, constraining spatially varying prestress and fault friction parameters by 30 strong motion and 12 GPS stations. From visiting 2 million models, we discern complex coseismic rupture dynamics that transition from a strongly radiating pulse‐like phase to a mildly radiating crack‐like phase. Both coseismic phases are separated by a shallow strength barrier that nearly arrests rupture and leads to a gap in the afterslip, reflecting the geologic heterogeneity along this segment of the SAF. Coseismic rupture termination involves distinct arrest mechanisms that imprint on afterslip kinematics. A backward propagating afterslip front may drive delayed aftershock activity above the hypocenter. Trade‐off analysis of the 10,500 best‐fitting models uncovers local correlations between prestress levels and the reference friction coefficient, alongside an anticorrelation between prestress and rate‐state parameters . We find that a complex, fault‐local interplay of dynamic parameters determines the nucleation, propagation, and arrest of both, co‐ and postseismic faulting. This study demonstrates the potential of inverse physics‐based modeling to reveal novel insights and detailed characterizations of well‐recorded earthquakes.more » « less
-
Several regularly recurring moderate-size earthquakes motivated dense instrumentation of the Parkfield section of the San Andreas fault, providing an invaluable near-fault observatory. We present a seismo-geodetic dynamic inversion of the 2004 Parkfield earthquake, which illuminates the interlinked complexity of faulting across time scales. Using fast-velocity-weakening rate-and-state friction, we jointly model 3D coseismic dynamic rupture and the 90-day evolution of postseismic slip. We utilize a parallel tempering Markov chain Monte Carlo approach to solve this non-linear high-dimensional inverse problem, constraining spatially varying prestress and fault friction parameters by 30 strong motion and 12 GPS stations. From visiting >2 million models, we discern complex coseismic rupture dynamics that transition from a strongly radiating pulse-like phase to a mildly radiating crack-like phase. Both coseismic phases are separated by a shallow strength barrier that nearly arrests rupture and leads to a gap in the afterslip. Coseismic rupture termination involves distinct arrest mechanisms that imprint on afterslip kinematics. A backward propagating afterslip front may drive delayed aftershock activity above the hypocenter. Analysis of the 10,500 best-fitting models uncovers local correlations between prestress levels and the reference friction coefficient, alongside an anticorrelation between prestress and rate-state parameters b−a. We find that a complex, fault-local interplay of dynamic parameters determines the nucleation, propagation, and arrest of both, co- and postseismic faulting. This study demonstrates the potential of inverse physics-based modeling to reveal novel insights and detailed characterizations of well-recorded earthquakes.more » « less
