Dynamic modeling of sequences of earthquakes and aseismic slip (SEAS) provides a self‐consistent, physics‐based framework to connect, interpret, and predict diverse geophysical observations across spatial and temporal scales. Amid growing applications of SEAS models, numerical code verification is essential to ensure reliable simulation results but is often infeasible due to the lack of analytical solutions. Here, we develop two benchmarks for three‐dimensional (3D) SEAS problems to compare and verify numerical codes based on boundary‐element, finite‐element, and finite‐difference methods, in a community initiative. Our benchmarks consider a planar vertical strike‐slip fault obeying a rate‐ and state‐dependent friction law, in a 3D homogeneous, linear elastic whole‐space or half‐space, where spontaneous earthquakes and slow slip arise due to tectonic‐like loading. We use a suite of quasi‐dynamic simulations from 10 modeling groups to assess the agreement during all phases of multiple seismic cycles. We find excellent quantitative agreement among simulated outputs for sufficiently large model domains and small grid spacings. However, discrepancies in rupture fronts of the initial event are influenced by the free surface and various computational factors. The recurrence intervals and nucleation phase of later earthquakes are particularly sensitive to numerical resolution and domain‐size‐dependent loading. Despite such variability, key properties of individual earthquakes, including rupture style, duration, total slip, peak slip rate, and stress drop, are comparable among even marginally resolved simulations. Our benchmark efforts offer a community‐based example to improve numerical simulations and reveal sensitivities of model observables, which are important for advancing SEAS models to better understand earthquake system dynamics.
- NSF-PAR ID:
- 10417587
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Date Published:
- Journal Name:
- Bulletin of the Seismological Society of America
- Volume:
- 113
- Issue:
- 2
- ISSN:
- 0037-1106
- Page Range / eLocation ID:
- 499 to 523
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract 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.
-
more » « less
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.
-
Abstract On 19 September 2022, a major earthquake struck the northwestern Michoacán segment along the Mexican subduction zone. A slip model is obtained that satisfactorily explains geodetic, teleseismic, and tsunami observations of the 2022 event. The preferred model has a compact large-slip patch that extends up-dip and northwestward from the hypocenter and directly overlaps a 1973 Mw 7.6 rupture. Slip is concentrated offshore and below the coast at depths from 10 to 30 km with a peak value of ∼2.9 m, and there is no detected coseismic slip near the trench. The total seismic moment is 3.1×1020 N·m (Mw 7.6), 72% of which is concentrated in the first 30 s. Most aftershocks are distributed in an up-dip area of the mainshock that has small coseismic slip, suggesting near-complete strain release in the large-slip patch. Teleseismic P waveforms of the 2022 and 1973 earthquakes are similar in duration and complexity with high cross-correlation coefficients of 0.68–0.98 for long P to PP signal time windows, indicating that the 2022 earthquake is a quasi-repeat of the 1973 earthquake, possibly indicating persistent frictional properties. Both the events produced more complex P waveforms than comparable size events along Guerrero and Oaxaca, reflecting differences in patchy locking of the Mexican megathrust.more » « less
-
more » « less
Abstract We present an efficient numerical method for earthquake sequences in 2D antiplane shear that incorporates wave propagation. A vertical strike‐slip fault governed by rate‐and‐state friction is embedded in a heterogeneous elastic half‐space discretized using a high‐order accurate Summation‐by‐Parts finite difference method. Adaptive time‐stepping is applied during the interseismic periods; during coseismic rupture we apply a non‐stiff method, enabling a variety of explicit time stepping methods. We consider a shallow sedimentary basin and explore sensitivity to spatial resolution and the switching criteria used to transition between solvers. For sufficient grid resolution and switching thresholds, simulations results remain robust over long time scales. We explore the effects of full dynamics and basin depth and stiffness, making comparisons with quasi‐dynamic counterparts. Fully‐dynamic ruptures generate higher stresses, faster slip rates and rupture speeds, producing seismic scattering in the bulk. Because single‐event dynamic simulations penetrate further into sediments compared to the quasi‐dynamic simulations, we hypothesize that the incorporation of inertial effects would produce sequences of only surface‐rupturing events. However, we find that subbasin ruptures can still emerge with elastodynamics, for sufficiently compliant basins. We also find that full dynamics can increase the frequency of surface‐rupturing events, depending on basin depth and stiffness. These results suggest that an earthquake's potential to penetrate into shallow sediments should be viewed through the lens of the earthquake sequence, as it depends on basin properties and wave‐mediated effects, but also on self‐consistent initial conditions obtained from seismogenic cycling.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1785/0120220066
