Search for: All records

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. 

Abstract Earthquakes are rupture-like processes that propagate along tectonic faults and cause seismic waves. The propagation speed and final area of the rupture, which determine an earthquake’s potential impact, are directly related to the nature and quantity of the energy dissipation involved in the rupture process. Here, we present the challenges associated with defining and measuring the energy dissipation in laboratory and natural earthquakes across many scales. We discuss the importance and implications of distinguishing between energy dissipation that occurs close to and far behind the rupture tip, and we identify open scientific questions related to a consistent modeling framework for earthquake physics that extends beyond classical Linear Elastic Fracture Mechanics. 

Abstract To better understand how normal stress heterogeneity affects earthquake rupture, we conducted laboratory experiments on a 760 mm poly (methyl‐mathacrylate) PMMA sample with a 25 mm “bump” of locally higher normal stress (∆σ_bt). We systematically varied the sample‐average normal stress () and bump prominence (). For bumps with lower prominence () the rupture simply propagated through the bump and produced regular sequences of periodic stick‐slip events. Bumps with higher prominence () produced complex rupture sequences with variable timing and ruptures sizes, and this complexity persisted for multiple stick‐slip supercycles. During some events, the bump remained locked and acted as a barrier that completely stopped rupture. In other events, a dynamic rupture front terminated at the locked bump, but rupture reinitiated on the other side of the bump after a brief pause of 0.3–1 ms. Only when stress on the bump was near critical did the bump slip and unload built up strain energy in one large event. Thus, a sufficiently prominent bump acted as a barrier (energy sink) when it was far from critically stressed and as an asperity (energy source) when it was near critically stressed. Similar to an earthquake gate, the bump never acted as a permanent barrier. In the experiments, we resolve the above rupture interactions with a bump as separate rupture phases; however, when observed through the lens of seismology, it may appear as one continuous rupture that speeds up and slows down. The complicated rupture‐bump interactions also produced enhanced high frequency seismic waves recorded with piezoelectric sensors. 

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. 

Abstract The interpretation of precursory seismicity can depend on a critical nucleation length scale h*, yet h* is largely unconstrained in the seismogenic crust. To estimate h* and associated earthquake nucleation processes at 2–7 km depths in Oklahoma, we studied seismic activity occurring prior to nine M 2.5–3.0 earthquakes that are aftershocks of the 3 September 2016 M 5.8 Pawnee, Oklahoma, earthquake. Four of the nine M 2.5–3.0 aftershocks studied did not have detectable seismicity within a 2 km radius of their hypocenters in the preceding 16 hr time windows. For the other five events, which did exhibit foreshock sequences, we estimated the static stress changes associated with each event of each sequence based on precise earthquake relocations and magnitude estimates. By carefully examining the spatiotemporal characteristics, we found all five of these M 2.5–3.0 aftershocks, and 70% of our studied events were plausibly triggered via static stress transfer from nearby earthquakes occurring hours to seconds earlier, consistent with the cascade nucleation model and a small h* in this region. The smallest earthquakes we could quantitatively study were M −1.5 events, which likely have 1–2 m rupture dimensions. The existence of these small events also supports a small nucleation length scale h*≤1  m, consistent with laboratory estimates. However, our observations cannot rule out more complicated earthquake initiation processes involving interactions between foreshocks and slow slip. Questions also remain as to whether aftershocks initiate differently from more isolated earthquakes. 

Abstract Earthquakes occur in clusters or sequences that arise from complex triggering mechanisms, but direct measurement of the slow subsurface slip responsible for delayed triggering is rarely possible. We investigate the origins of complexity and its relationship to heterogeneity using an experimental fault with two dominant seismic asperities. The fault is composed of quartz powder, a material common to natural faults, sandwiched between 760 mm long polymer blocks that deform the way 10 meters of rock would behave. We observe periodic repeating earthquakes that transition into aperiodic and complex sequences of fast and slow events. Neighboring earthquakes communicate via migrating slow slip, which resembles creep fronts observed in numerical simulations and on tectonic faults. Utilizing both local stress measurements and numerical simulations, we observe that the speed and strength of creep fronts are highly sensitive to fault stress levels left behind by previous earthquakes, and may serve as on-fault stress meters.