Earthquake energy dissipation in a fracture mechanics framework.

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.

The complex dynamics of the 2023 Kahramanmaras, Turkey, Mw 7.8-7.7 earthquake doublet.

The destructive 2023 moment magnitude (Mw) 7.8-7.7 earthquake doublet ruptured multiple segments of the East Anatolian Fault system in Turkey. We integrated multiscale seismic and space-geodetic observations with multifault kinematic inversions and dynamic rupture modeling to unravel the events' complex rupture history and stress-mediated fault interactions. Our analysis reveals three subshear slip episodes during the initial Mw 7.8 earthquake with a delayed rupture initiation to the southwest. The Mw 7.7 event occurred 9 hours later with a larger slip and supershear rupture on its western branch. Mechanically consistent dynamic models accounting for fault interactions can explain the unexpected rupture paths and require a heterogeneous background stress. Our results highlight the importance of combining near- and far-field observations with data-driven and physics-based models for seismic hazard assessment.

Dynamics, interactions and delays of the 2019 Ridgecrest rupture sequence.

The observational difficulties and the complexity of earthquake physics have rendered seismic hazard assessment largely empirical. Despite increasingly high-quality geodetic, seismic and field observations, data-driven earthquake imaging yields stark differences and physics-based models explaining all observed dynamic complexities are elusive. Here we present data-assimilated three-dimensional dynamic rupture models of California's biggest earthquakes in more than 20 years: the moment magnitude (Mw) 6.4 Searles Valley and Mw 7.1 Ridgecrest sequence, which ruptured multiple segments of a non-vertical quasi-orthogonal conjugate fault system1. Our models use supercomputing to find the link between the two earthquakes. We explain strong-motion, teleseismic, field mapping, high-rate global positioning system and space geodetic datasets with earthquake physics. We find that regional structure, ambient long- and short-term stress, and dynamic and static fault system interactions driven by overpressurized fluids and low dynamic friction are conjointly crucial to understand the dynamics and delays of the sequence. We demonstrate that a joint physics-based and data-driven approach can be used to determine the mechanics of complex fault systems and earthquake sequences when reconciling dense earthquake recordings, three-dimensional regional structure and stress models. We foresee that physics-based interpretation of big observational datasets will have a transformative impact on future geohazard mitigation.

Assessing Margin-Wide Rupture Behaviors Along the Cascadia Megathrust With 3-D Dynamic Rupture Simulations.

From California to British Columbia, the Pacific Northwest coast bears an omnipresent earthquake and tsunami hazard from the Cascadia subduction zone. Multiple lines of evidence suggests that magnitude eight and greater megathrust earthquakes have occurred - the most recent being 321 years ago (i.e., 1700 A.D.). Outstanding questions for the next great megathrust event include where it will initiate, what conditions are favorable for rupture to span the convergent margin, and how much slip may be expected. We develop the first 3-D fully dynamic rupture simulations for the Cascadia subduction zone that are driven by fault stress, strength and friction to address these questions. The initial dynamic stress drop distribution in our simulations is constrained by geodetic coupling models, with segment locations taken from geologic analyses. We document the sensitivity of nucleation location and stress drop to the final seismic moment and coseismic subsidence amplitudes. We find that the final earthquake size strongly depends on the amount of slip deficit in the central Cascadia region, which is inferred to be creeping interseismically, for a given initiation location in southern or northern Cascadia. Several simulations are also presented here that can closely approximate recorded coastal subsidence from the 1700 A.D. event without invoking localized high-stress asperities along the down-dip locked region of the megathrust. These results can be used to inform earthquake and tsunami hazards for not only Cascadia, but other subduction zones that have limited seismic observations but a wealth of geodetic inference.

Dynamic viability of the 2016 Mw 7.8 Kaikoura earthquake cascade on weak crustal faults.

We present a dynamic rupture model of the 2016 Mw 7.8 Kaikoura earthquake to unravel the event's riddles in a physics-based manner and provide insight on the mechanical viability of competing hypotheses proposed to explain them. Our model reproduces key characteristics of the event and constraints puzzling features inferred from high-quality observations including a large gap separating surface rupture traces, the possibility of significant slip on the subduction interface, the non-rupture of the Hope fault, and slow apparent rupture speed. We show that the observed rupture cascade is dynamically consistent with regional stress estimates and a crustal fault network geometry inferred from seismic and geodetic data. We propose that the complex fault system operates at low apparent friction thanks to the combined effects of overpressurized fluids, low dynamic friction and stress concentrations induced by deep fault creep.