Dynamic Rupture Simulations Research Articles

Fault rock properties and conditions produce variance in slip during earthquake rupture propagation at the Nankai Trough

Although drilled samples of fault rocks have yielded information on frictional features of shallow subduction zones, the relationship of rupture propagation to the levels of friction and pore-fluid pressure remains uncertain. To investigate this topic, we performed dynamic rupture simulations along the megasplay fault that slipped during the 1944 Mw 8.0 Tonankai earthquake in the Nankai Trough. We used actual data from friction experiments on rocks from the fault segment and pre-existing pore pressures deduced from geophysical surveys for the shallow portion of 0–10 km depth along the fault. Simulations of low friction (friction coefficient ca. 0.04) produced large slip (about 30 m), whereas simulations using higher friction (friction coefficient ca. 0.2) suppressed the rupture. In simulations with low friction in which the pore-fluid pressure was nearly equal to the lithostatic stress, the slip decreased to about 25 m. However, when the simulations included slip-strengthening at shallow depth and higher friction, the slip still reached roughly 20 m. Such variability in slip during rupture propagation is caused by differences in the friction features and fluid pressure conditions of fault rocks, in which the friction features might be related to the mineral composition. Spatiotemporal heterogeneity in fault-rock type and their physical and hydraulic properties may fundamentally produce the complexity and variability of earthquake rupture propagation along the Nankai plate-subduction boundary.

Strong asymmetry in near-fault ground velocity during an oblique strike-slip earthquake revealed by waveform particle motions and dynamic rupture simulations

Strong asymmetry in near-fault ground velocity during an oblique strike-slip earthquake revealed by waveform particle motions and dynamic rupture simulations

Physics-based assessment of earthquake potential on the Anninghe-Zemuhe fault system in southwestern China

The seismic hazard of a fault system is controlled by the maximum possible earthquake magnitude it can host. However, existing methods to estimate maximum magnitudes can result in large uncertainties or ignore their temporal evolution. Here, we show how the maximum possible earthquake magnitude of a fault system can be assessed by combining high-resolution fault coupling maps with a physics-based model from three-dimensional dynamic fracture mechanics confirmed by dynamic rupture simulations. We demonstrate the method on the Anninghe-Zemuhe fault system in southwestern China, where dense near-fault geodetic data has been acquired. Our results show that this fault system currently has the potential to generate Mw7.0 earthquakes with maximum magnitudes increasing to Mw7.3 by 2200. These results are supported by the observed rupture extents and recurrence times of historical earthquakes and the b values of current seismicity. Our work provides a practical way to assess the earthquake potential of natural faults.

The 2022 Mw 6.6 Menyuan earthquake: An early-terminated runaway rupture by the complex fault geometry

The 2022 Mw 6.6 Menyuan earthquake, which caused unexpectedly large surface offsets (about 2.6–3.5 m) relative to the moment magnitude of this event, is one of the best-recorded strike-slip events in the northeastern Tibetan plateau. Although previous studies have revealed many kinematic features of this event, the mechanism and rupture dynamic characteristics underlying unexpectedly large surface dislocations need to be clarified. In this study, based on the grid search method, we conducted 61 dynamic rupture simulations considering terrain relief to screen out a preferred dynamic model that fits the geodetic, seismic, and geological observations well. The results show that the average stress drop of the nucleation part is relatively high (11 MPa), which could provide high initial rupture energy and give this event a chance to develop into a large earthquake. However, the rupture propagating eastward is slowed and terminated by two consecutive fault bends of about 10° each, and the westward propagating one might concentrate at a shallow depth (< 3 km) on the south secondary fault after it jumps over the stepover. These results imply that the complex fault geometry (the fault bend and stepover) might terminate the rupture early so that the observed magnitude of this event (Mw 6.6) is smaller than the estimated magnitude (Mw 7.1–7.2) based on the surface offsets. Our study again shows that, due to advances in computational power and observations (coverage and quality), physics-based dynamic rupture simulations have entered an era where it is possible to fit the observed data well, which is vital for further understanding the mechanisms of earthquakes.

Numerical simulation of 3-D seismic wave based on alternative flux finite-difference WENO scheme

SUMMARY High-frequency non-physical oscillations may occur due to shock waves in seismic wavefield and dynamic rupture simulation. In this study, we introduced the alternative flux finite-difference weighted essentially non-oscillatory scheme to address potential shock wave issues in computational seismology effectively. The wavefield of the body-fitted curvilinear domain was accurately computed through conservative grid mapping, ensuring accurate implementation of free surface boundary conditions on irregular surfaces using characteristic boundary conditions and minimizing artificial boundary reflections with exponential decay absorbing layers. Finally, we compared our scheme with the GRTM for flat surfaces and the CGFDM3D-EQR for irregular surfaces to demonstrate its correctness and accuracy, and validated its non-oscillatory characteristics. The aforementioned scheme is anticipated to assume a significant function in simulating more intricate seismic wavefields or dynamic ruptures.

Linking 3D Long‐Term Slow‐Slip Cycle Models With Rupture Dynamics: The Nucleation of the 2014 Mw 7.3 Guerrero, Mexico Earthquake

AbstractSlow slip events (SSEs) have been observed in spatial and temporal proximity to megathrust earthquakes in various subduction zones, including the 2014 Mw 7.3 Guerrero, Mexico earthquake which was preceded by a Mw 7.6 SSE. However, the underlying physics connecting SSEs to earthquakes remains elusive. Here, we link 3D slow‐slip cycle models with dynamic rupture simulations across the geometrically complex flat‐slab Cocos plate boundary. Our physics‐based models reproduce key regional geodetic and teleseismic fault slip observations on timescales from decades to seconds. We find that accelerating SSE fronts transiently increase shear stress at the down‐dip end of the seismogenic zone, modulated by the complex geometry beneath the Guerrero segment. The shear stresses cast by the migrating fronts of the 2014 Mw 7.6 SSE are significantly larger than those during the three previous episodic SSEs that occurred along the same portion of the megathrust. We show that the SSE transient stresses are large enough to nucleate earthquake dynamic rupture and affect rupture dynamics. However, additional frictional asperities in the seismogenic part of the megathrust are required to explain the observed complexities in the coseismic energy release and static surface displacements of the Guerrero earthquake. We conclude that it is crucial to jointly analyze the long‐ and short‐term interactions and complexities of SSEs and megathrust earthquakes across several (a)seismic cycles accounting for megathrust geometry. Our study has important implications for identifying earthquake precursors and understanding the link between transient and sudden megathrust faulting processes.

Rupture Dynamics of Cascading Earthquakes in a Multiscale Fracture Network

AbstractFault‐damage zones comprise multiscale fracture networks that may slip dynamically and interact with the main fault during earthquake rupture. Using 3D dynamic rupture simulations and scale‐dependent fracture energy, we examine dynamic interactions of more than 800 intersecting multiscale fractures surrounding a listric fault, emulating a major listric fault and its damage zone. We investigate 10 distinct orientations of maximum horizontal stress, probing the conditions necessary for sustained slip within the fracture network or activating the main fault. Additionally, we assess the feasibility of nucleating dynamic rupture earthquake cascades from a distant fracture and investigate the sensitivity of fracture network cascading rupture to the effective normal stress level. We model either pure cascades or main fault rupture with limited off‐fault slip. We find that cascading ruptures within the fracture network are dynamically feasible under certain conditions, including: (a) the fracture energy scales with fracture and fault size, (b) favorable relative pre‐stress of fractures within the ambient stress field, and (c) close proximity of fractures. We find that cascading rupture within the fracture network discourages rupture on the main fault. Our simulations suggest that fractures with favorable relative pre‐stress, embedded within a fault damage zone, may lead to cascading earthquake rupture that shadows main fault slip. We find that such off‐fault events may reach moment magnitudes up to Mw ≈ 5.5, comparable to magnitudes that can be otherwise hosted by the main fault. Our findings offer insights into physical processes governing cascading earthquake dynamic rupture within multiscale fracture networks.

Fractal Spatial Distributions of Initial Shear Stress and Frictional Properties on Faults and Their Impact on Dynamic Earthquake Rupture

ABSTRACT We investigate the influence of the heterogeneous slip-weakening distance (DC) in dynamic rupture simulations, in which DC is proportional to the fault irregularities. Specifically, we compare a heterogeneous fractal DC distribution to a uniform DC over the entire fault when the initial shear stress is also heterogeneous. We find that even small changes in the average value of DC (&amp;lt;1 mm) can lead to significant differences in the rupture evolution; that is, the average DC and the way DC is distributed determines if the rupture is a runaway, self-arrested, or nonpropagating. We find that the self-arrested ruptures differ from runaway ruptures in the amount of area characterized by large slips (asperities). Self-arrested ruptures match the Somerville et al. (1999) asperity criteria in which ∼25% of ruptured area radiate ∼45% of the total seismic moment. This criterion is not satisfied for runaway ruptures. For runaway ruptures, ∼50% of the ruptured area radiates about 70% of the seismic moment, indicating that the ruptured area is not linearly proportional to the seismic moment. Self-arrested ruptures are characterized by dynamic shear stress drops (SDs) in the range ∼2.9–5.5 MPa, whereas for runaway ruptures the dynamic SDs increase to values between ∼12 and 20 MPa. Self-arrested ruptures generated by fractal distributed DC resemble the rupture properties of observed earthquakes. In addition, results show that the conditions for self-arrested ruptures are connected to the decrease of residual energy at rupture boundaries.

Linked and fully coupled 3D earthquake dynamic rupture and tsunami modeling for the Húsavík–Flatey Fault Zone in North Iceland

Abstract. Tsunamigenic earthquakes pose considerable risks, both economically and socially, yet earthquake and tsunami hazard assessments are typically conducted separately. Earthquakes associated with unexpected tsunamis, such as the 2018 Mw 7.5 strike-slip Sulawesi earthquake, emphasize the need to study the tsunami potential of active submarine faults in different tectonic settings. Here, we investigate physics-based scenarios combining simulations of 3D earthquake dynamic rupture and seismic wave propagation with tsunami generation and propagation. We present time-dependent modeling of one-way linked and 3D fully coupled earthquakes and tsunamis for the ∼ 100 km long Húsavík–Flatey Fault Zone (HFFZ) in North Iceland. Our analysis shows that the HFFZ has the potential to generate sizable tsunamis. The six dynamic rupture models sourcing our tsunami scenarios vary regarding hypocenter location, spatiotemporal evolution, fault slip, and fault structure complexity but coincide with historical earthquake magnitudes. Earthquake dynamic rupture scenarios on a less segmented fault system, particularly with a hypocenter location in the eastern part of the fault system, have a larger potential for local tsunami generation. Here, dynamically evolving large shallow fault slip (∼ 8 m), near-surface rake rotation (± 20∘), and significant coseismic vertical displacements of the local bathymetry (± 1 m) facilitate strike-slip faulting tsunami generation. We model tsunami crest to trough differences (total wave heights) of up to ∼ 0.9 m near the town Ólafsfjörður. In contrast, none of our scenarios endanger the town of Akureyri, which is shielded by multiple reflections within the narrow Eyjafjörður bay and by Hrísey island. We compare the modeled one-way linked tsunami waveforms with simulation results using a 3D fully coupled approach. We find good agreement in the tsunami arrival times and location of maximum tsunami heights. While seismic waves result in transient motions of the sea surface and affect the ocean response, they do not appear to contribute to tsunami generation. However, complex source effects arise in the fully coupled simulations, such as tsunami dispersion effects and the complex superposition of seismic and acoustic waves within the shallow continental shelf of North Iceland. We find that the vertical velocity amplitudes of near-source acoustic waves are unexpectedly high – larger than those corresponding to the actual tsunami – which may serve as a rapid indicator of surface dynamic rupture. Our results have important implications for understanding the tsunamigenic potential of strike-slip fault systems worldwide and the coseismic acoustic wave excitation during tsunami generation and may help to inform future tsunami early warning systems.

Propagation of extended fractures by local nucleation and rapid transverse expansion of crack-front distortion.

Fractures are ubiquitous and can lead to the catastrophic material failure of materials. Although fracturing in a two-dimensional plane is well understood, all fractures are extended in and propagate through three-dimensional space. Moreover, their behaviour is complex. Here we show that the forward propagation of a fracture front occurs through an initial rupture, nucleated at some localized position, followed by a very rapid transverse expansion at velocities as high as the Rayleigh-wave speed. We study fracturing in a circular geometry that achieves an uninterrupted extended fracture front and use a fluid to control the loading conditions that determine the amplitude of the forward jump. We find that this amplitude correlates with the transverse velocity. Dynamic rupture simulations capture the observations for only a high transverse velocity. These results highlight the importance of transverse dynamics in the forward propagation of an extended fracture.

A systematic scheme to develop dynamic earthquake rupture scenarios: a case study on the Wenchuan–Maoxian Fault in the Longmen Shan, China, thrust belt

The 2008 Wenchuan Mw 7.9 mainshock caused catastrophic destruction to cities along the northwestern margin of the Sichuan Basin. This earthquake did not activate the Wenchuan–Maoxian Fault (WMF) on the hinterland side and the conjugate buried Lixian Fault (LXF), but they could experience large earthquakes in the future. We propose a systematic scheme to develop scenario earthquakes for active fault systems with insufficient constrain of 3D fault geometries. We first performed stress tensor inversion to constrain the regional stress field. Then, we developed a new method to constrain fault geometries by inverting long-term slip rates under the given regional stress and applied it to the WMF. We conducted a set of 3D dynamic earthquake rupture simulations on the WMF and LXF to assess the scenarios of earthquake rupture processes. Several fault nucleation points, friction coefficients, and initial stress states are assessed, the general rupture patterns for these earthquake scenarios are evaluated, and finally, we find the scenarios that could fall into three groups. Depending on initial conditions, the dynamic rupture may start in the LXF, leading to magnitude-7.0 earthquakes, or start in the WMF, then cascade through the LXF, leading to magnitude-7.5 earthquakes, or both start and arrest in the WMF, leading to around magnitude-6.5 or -7.0 earthquakes. We find that the rupture starting on the reverse oblique-slip jumps to the strike-slip fault, but the reverse process is impeded.Graphical

The effect of crustal fluid for the spontaneous dynamic fault rupture simulation

The effect of crustal fluid for the spontaneous dynamic fault rupture simulation

Dynamic Rupture of the 2021 MW 7.4 Maduo Earthquake: An Intra‐Block Event Controlled by Fault Geometry

AbstractThe 2021 MW 7.4 Maduo event occurred within the Bayan Har block in eastern Tibet, which provides an opportunity to investigate the stress conditions and rheology of faults within the block. Here, we perform dynamic rupture simulations based on the finite element method to explore the physical conditions underlying this earthquake and the factors that controlled the rupture process. We construct the model with a nonplanar fault inferred from the Interferometric Synthetic Aperture Radar (InSAR) data and aftershocks sequence relocation. Our dynamic model is controlled by slip‐weakening friction law with initial stress on fault resolved from a uniform regional stress field. The preferred model produces an average slip of ∼2.2 m with a maximum slip of ∼4.0 m. There are three asperities distributed along the strike, which have captured the main features of the Maduo event. The simulation results are consistent with the static GPS coseismic surface displacements, InSAR data, and displacement waveforms recorded by high‐rate GNSS stations. By comparing the results with the planar fault model and rotated stress fields, we find that the fault geometry and regional stress field are the primary factors that control the rupture process of the event. Moreover, we infer that the unfavorable orientation and fault bend lead to minor slips on the branch fault. Furthermore, we investigate the potential mechanisms of supershear rupture on the eastern fault segment.

Ground-Motion Variability for Ruptures on Rough Faults

ABSTRACT Fault roughness influences earthquake rupture dynamics, seismic energy radiation, and, hence, resulting ground motion and its variability. Using 3D dynamic rupture simulations considering a range of rough-fault realizations, we investigate the effects of rupture complexity caused by fault roughness on ground-motion variability, that is, the variability of peak ground acceleration (PGA) and velocity (PGV) as a function of distance. In our analysis, we vary hypocenter locations (leading to unilateral and bilateral ruptures) and fault roughness amplitude to generate a set of magnitude M ≈ 7 strike-slip dynamic rupture simulations. Synthetic seismic waveforms computed on a dense set of surface sites (maximum resolved frequency 5.75 Hz) form our database for detailed statistical analyses. For unilateral ruptures, our simulations reveal that ground-shaking variability (in terms of PGA and PGV) remains nearly constant with increasing distance from the fault. In contrast, bilateral ruptures lead to slowly decreasing ground-motion variability with increasing distance in the near field (less than 20 km). The variability becomes almost constant at large fault distances. We also find that low-amplitude fault roughness leads to ruptures that are likely to generate higher PGA variability than events on faults with high-amplitude roughness. Increasing fault roughness distorts the radiation pattern, thereby reducing directivity effects and, hence, potentially lowering ground-motion variability. The average PGV variability from our rough-fault rupture models is consistent with estimates from empirical ground-motion models (GMMs). However, the average PGA variability exceeds the variability encoded in empirical GMMs by nearly 20%. Hence, our findings have implications for near-source ground-motion prediction in seismic hazard studies, because ground-motion variability depends on details of the earthquake rupture process and is larger than GMM estimates.

Equivalent Near-Field Corner Frequency Analysis of 3D Dynamic Rupture Simulations Reveals Dynamic Source Effects

Abstract Dynamic rupture simulations generate synthetic waveforms that account for nonlinear source and path complexity. Here, we analyze millions of spatially dense waveforms from 3D dynamic rupture simulations in a novel way to illuminate the spectral fingerprints of earthquake physics. We define a Brune-type equivalent near-field corner frequency (fc) to analyze the spatial variability of ground-motion spectra and unravel their link to source complexity. We first investigate a simple 3D strike-slip setup, including an asperity and a barrier, and illustrate basic relations between source properties and fc variations. Next, we analyze &amp;gt;13,000,000 synthetic near-field strong-motion waveforms generated in three high-resolution dynamic rupture simulations of real earthquakes, the 2019 Mw 7.1 Ridgecrest mainshock, the Mw 6.4 Searles Valley foreshock, and the 1992 Mw 7.3 Landers earthquake. All scenarios consider 3D fault geometries, topography, off-fault plasticity, viscoelastic attenuation, and 3D velocity structure and resolve frequencies up to 1–2 Hz. Our analysis reveals pronounced and localized patterns of elevated fc, specifically in the vertical components. We validate such fc variability with observed near-fault spectra. Using isochrone analysis, we identify the complex dynamic mechanisms that explain rays of elevated fc and cause unexpectedly impulsive, localized, vertical ground motions. Although the high vertical frequencies are also associated with path effects, rupture directivity, and coalescence of multiple rupture fronts, we show that they are dominantly caused by rake-rotated surface-breaking rupture fronts that decelerate due to fault heterogeneities or geometric complexity. Our findings highlight the potential of spatially dense ground-motion observations to further our understanding of earthquake physics directly from near-field data. Observed near-field fc variability may inform on directivity, surface rupture, and slip segmentation. Physics-based models can identify “what to look for,” for example, in the potentially vast amount of near-field large array or distributed acoustic sensing data.

A seismologically consistent surface rupture length model for unbounded and width-limited events

A new surface-rupture-length (SRL) relationship as a function of magnitude [Formula: see text], fault thickness, and fault dip angle is presented in this article. The objective of this study is to model the change in scaling between unbounded and width-limited ruptures. This is achieved through the use of seismological-theory-based relationships for the average displacement scaling and the aid of dynamic fault rupture simulations to constrain the rupture width scaling. The empirical data set used in the development of this relationship is composed of [Formula: see text] events ranging from [Formula: see text] 5 to [Formula: see text] and SRL 1.1 to 432 km. The dynamic rupture simulations data set includes [Formula: see text] events ranging from [Formula: see text] 4.9 to 8.2 and [Formula: see text] 1 to 655 km. For the average displacement [Formula: see text] scaling, three simple models and two composite models were evaluated. The simple average displacement models are a square root of the rupture area [Formula: see text], a down-dip width [Formula: see text], and a rupture length [Formula: see text] proportional model. The two composite models have a [Formula: see text] scaling for unbounded ruptures and transition to [Formula: see text] and [Formula: see text] scaling for width-limited events, respectively. The empirical data favors a [Formula: see text] scaling for the entire regime (unbounded and width-limited ruptures) followed by a [Formula: see text] scaling for unbounded that changes to [Formula: see text] scaling for width-limited ruptures. The selected models exhibit better predictive performance compared to linear [Formula: see text] type models, especially in the large magnitude range, which is dominated by width-limited events. A comparison with published SRL models shows consistent scaling for different fault types and tectonic environments.

The sharp turn: Backward rupture branching during the 2023 Mw 7.8 Kahramanmaraş (Türkiye) earthquake

Multiple lines of evidence indicate that the 2023 Mw 7.8 Kahramanmaraş (Türkiye) earthquake started on a splay fault, then branched bilaterally onto the nearby East Anatolian Fault (EAF). This rupture pattern includes one feature previously deemed implausible, called backward rupture branching: rupture propagating from the splay fault onto the SW EAF segment through a sharp corner (with an acute angle between the two faults). To understand this feature, we perform 2.5-D dynamic rupture simulations considering a large set of possible scenarios. We find that both subshear and supershear ruptures on the splay fault can trigger bilateral ruptures on the EAF, which themselves can be either subshear, supershear, or a mixture of the two. In most cases, rupture on the SW segment of the EAF starts after rupture onset on its NE segment: the SW rupture is triggered by the NE rupture. Only when the EAF has initial stresses very close to failure can its SW segment be directly triggered by the initial splay-fault rupture, earlier than the activation of the NE segment. These results advance our understanding of the mechanisms of multi-segment rupture and the complexity of rupture processes, paving the way for a more accurate assessment of earthquake hazards.

3D Dynamic Rupture Modeling of the 6 February 2023, Kahramanmaraş, Turkey Mw 7.8 and 7.7 Earthquake Doublet Using Early Observations

Abstract The 2023 Turkey earthquake sequence involved unexpected ruptures across numerous fault segments. We present 3D dynamic rupture simulations to illuminate the complex dynamics of the earthquake doublet. Our models are constrained by observations available within days of the sequence and deliver timely, mechanically consistent explanations of the unforeseen rupture paths, diverse rupture speeds, multiple slip episodes, heterogeneous fault offsets, locally strong shaking, and fault system interactions. Our simulations link both earthquakes, matching geodetic and seismic observations and reconciling regional seismotectonics, rupture dynamics, and ground motions of a fault system represented by 10 curved dipping segments and embedded in a heterogeneous stress field. The Mw 7.8 earthquake features delayed backward branching from a steeply branching splay fault, not requiring supershear speeds. The asymmetrical dynamics of the distinct, bilateral Mw 7.7 earthquake are explained by heterogeneous fault strength, prestress orientation, fracture energy, and static stress changes from the previous earthquake. Our models explain the northward deviation of its eastern rupture and the minimal slip observed on the Sürgü fault. 3D dynamic rupture scenarios can elucidate unexpected observations shortly after major earthquakes, providing timely insights for data-driven analysis and hazard assessment toward a comprehensive, physically consistent understanding of the mechanics of multifault systems.

The break of earthquake asperities imaged by distributed acoustic sensing.

Rupture imaging of megathrust earthquakes with global seismic arrays revealed frequency-dependent rupture signatures1-4, but the role of high-frequency radiators remains unclear3-5. Similar observations of the more abundant crustal earthquakes could provide critical constraints but are rare without ultradense local arrays6,7. Here we use distributed acousticsensing technology8,9 to image the high-frequency earthquake rupture radiators. By converting a 100-kilometre dark-fibre cable into a 10,000-channel seismic array, we image four high-frequency subevents for the 2021 Antelope Valley, California, moment-magnitude 6.0 earthquake. After comparing our results with long-period moment-release10,11 and dynamic rupture simulations, we suggest that the imaged subevents are due to the breaking of fault asperities-stronger spots or pins on the fault-that substantially modulate the overall rupture behaviour. An otherwise fading rupture propagation could be promoted by the breaking of fault asperities in a cascading sequence. This study highlights how we can use the extensive pre-existing fibre networks12 as high-frequency seismic antennas to systematically investigate the rupture process of regional moderate-sized earthquakes. Coupled with dynamic rupture modelling, it could improve our understanding of earthquake rupture dynamics.

The 2016 Menyuan Earthquake: The Largest Self‐Arrested Crustal Earthquake Ever Observed

AbstractThe 2016 Mw 5.9 Menyuan earthquake occurred near the Haiyuan fault system in the northeastern Tibetan plateau. Although some kinematic properties of this event have been investigated, more dynamic properties still need to be clarified. First, we use the empirical Green's function method to obtain this event's apparent source time functions (ASTFs), which present a bell shape with a duration of 3s and corresponding smooth spectra. We then conduct dynamic rupture simulations, which are well constrained by the seismic and geodetic observations, to explore more in‐depth source properties. The results show that the characteristics of ASTFs can be explained by the self‐arresting rupture model, which implies that the rupture propagation of this event may spontaneously terminate before reaching the barriers. This result suggests that the initial conditions when this earthquake nucleated may have controlled the rupture process. Our work contributes to understanding the stopping mechanism of some moderate earthquake ruptures.