Search for: All records

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 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 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 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. 

SUMMARY Earthquake ruptures are generally considered to be cracks that propagate as fracture or frictional slip on pre-existing faults. Crack models have been used to describe the spatial distribution of fault offset and the associated static stress changes along a fault, and have implications for friction evolution and the underlying physics of rupture processes. However, field measurements that could help refine idealized crack models are rare. Here, we describe large-scale laboratory earthquake experiments, where all rupture processes were contained within a 3-m long saw-cut granite fault, and we propose an analytical crack model that fits our measurements. Similar to natural earthquakes, laboratory measurements show coseismic slip that gradually tapers near the rupture tips. Measured stress changes show roughly constant stress drop in the centre of the ruptured region, a maximum stress increase near the rupture tips and a smooth transition in between, in a region we describe as the earthquake arrest zone. The proposed model generalizes the widely used elliptical crack model by adding gradually tapered slip at the ends of the rupture. Different from the cohesive zone described by fracture mechanics, we propose that the transition in stress changes and the corresponding linear taper observed in the earthquake arrest zone are the result of rupture termination conditions primarily controlled by the initial stress distribution. It is the heterogeneous initial stress distribution that controls the arrest of laboratory earthquakes, and the features of static stress changes. We also performed dynamic rupture simulations that confirm how arrest conditions can affect slip taper and static stress changes. If applicable to larger natural earthquakes, this distinction between an earthquake arrest zone (that depends on stress conditions) and a cohesive zone (that depends primarily on strength evolution) has important implications for how seismic observations of earthquake fracture energy should be interpreted. 

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 scaleL_c, 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 withL_c< 1 m, to ignition‐limited environments withL_c> 10 km. In the latter case,L_cfails 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.