Rupture Speed Research Articles

Effects of Dip Angle on Rupture Propagation Along Branch Fault Systems

ABSTRACT An important consideration in assessing seismic hazards is determining what is likely to happen when an earthquake rupture encounters a geometric complexity such as a branch fault. Previous studies showed parameters such as branch angle, stress orientation, and stress heterogeneity as key factors in the self-determined rupture path on branch faults. However, most of these studies were conducted in 2D or 3D with perfectly vertical faults. Therefore, in this study, we investigate the effects of dipping angle on rupture propagation along a branch fault system. We construct 3D finite-element meshes where we vary the dip angles (nine geometries in total) of the main and secondary faults, the stressing angle (Ψ=20°, 40°, and 65°), and the hypocenter location with nucleation on both the main and secondary segments. We find that for Ψ=40°, a rupture on the main fault is most likely to propagate across the branch intersection when the secondary fault is dipping. In addition, for Ψ=65°, a rupture on the secondary fault is most likely to propagate to the main fault when the secondary fault is shallowly dipping. This is caused by a fast rupture speed on the secondary fault and the dynamic stress effect that develops with the interaction of the free surface and the dipping secondary fault. These results indicate that dip angle is an important parameter in the determination of rupture path on branch fault systems, with potentially significant impact for seismic hazard, and should be considered in future dynamic rupture modeling studies.

Abnormal apparent supershear rupture with discontinuous jumping propagation during the 2023 Türkiye M7.8 earthquake

Earthquake ruptures along a single fault or along a connected system of faults are generally assumed to progress continuously. However, our analysis of the 2023 M7.8 Türkiye earthquake, using finite-fault joint source inversion, uncovered the occurrence of discontinuous rupture jumps. The main fault area adjacent to the splay fault where the earthquake started, and the deeper portion of the northeastern main fault segment exhibited triggered slip before the main rupture front arrived. Through seismic centroid analysis and finite-fault inversion, we estimated apparent rupture speeds within these slip patches reach approximately 6.0 km s-1, exceeding local S-wave velocity. The dynamic triggering mechanism induced the jumping rupture in these areas, resulting in an apparent rupture velocity surpassing the local shear wave velocity. These findings demonstrate the importance of dynamic triggering in adjacent fault systems during large earthquakes, influencing the extent and complexity of rupture propagation.

Dual effect of roughness during earthquake rupture sequences on faults with strongly rate-weakening friction

We numerically model earthquake rupture cycles on non-planar faults governed by rate and state friction with an enhanced dynamic weakening in the form of flash heating, investigating quantitatively the relationship between the fault geometry, stresses from the preceding earthquakes, and the rupture characteristics. On planar faults, earthquakes with similar magnitudes and stress drops periodically rupture the entire fault, propagating under a low background shear-to-normal stress ratio as self-healing slip pulses due to the large friction reduction, with a sub-Rayleigh rupture speed. Faults with low roughness levels generally show a similar pattern. Prior to some events, the stress ratio along the fault slightly increases, leading to ruptures with secondary slip pulses. As roughness increases, the stresses become more heterogeneous, leading to a more complex cycle with part of the ruptures arresting at restraining bends with a low stress ratio. Sequences of several partial ruptures lead to a significant stress accumulation on the locked fault segments. These are eventually released by crack-like earthquake ruptures with supershear propagation speed and significantly larger magnitudes, stress drops, and radiated energy than typical planar fault events. Therefore, while fault roughness provides a mechanism for rupture arrest, consistent with previous studies, it can also result in a substantial increase in earthquake magnitudes, which should be considered in hazard assessments.

Fault Orientation Trumps Fault Maturity in Controlling Coseismic Rupture Characteristics of the 2021 Maduo Earthquake

AbstractFault maturity has been proposed to exert a first order control on earthquake rupture, yet direct observations linking individual rupture to long‐term fault growth are rare. The 2021 Mw 7.4 Maduo earthquake ruptured the east‐growing end of the slow‐moving (∼1 mm/yr) Jiangcuo fault in north Tibet, providing an opportunity to examine the relation between rupture characteristics and fault structure. Here we combine field and multiple remote sensing techniques to map the surface rupture at cm‐resolution and document comprehensively on‐fault offsets and off‐fault deformation. The 158 km‐long surface rupture consists of misoriented structurally inherited N110°‐striking segments and younger optimally oriented N093°‐striking segments, relative to the regional stress field. Despite being comparatively newly formed, the ∼N093°‐striking fault segments accommodate more localized strain, with up to 3 m on‐fault left‐lateral slip and 25%–50% off‐fault deformation, and possibly faster rupture speed. These results are in contrast with previous findings showing more localized strain and faster rupture speed on more mature fault segments; instead, our observations suggest that fault orientation with respect to the regional stress can exert a more important control than fault maturity on coseismic rupture behavior when both factors are at play.

Along-Strike Variation of Rupture Characteristics and Aftershock Patterns of the 2023 Mw 7.8 Türkiye Earthquake Controlled by Fault Structure

Abstract On 6 February 2023, an Mw 7.8 earthquake occurred along the East Anatolian fault zone (EAFZ) in southeastern Türkiye, representing the strongest earthquake in the region in nearly 80 yr. We investigate rupture characteristics and aftershock patterns of the earthquake through focal mechanism calculation, backprojection analysis, and finite-fault inversion. The results show bilateral rupture propagation of the mainshock with transient supershear speed in the southwest portion of the EAFZ, as well as shallower coseismic slip and abundant normal-faulting aftershocks in the same portion. We attribute these earthquake behaviors to the along-strike variation of fault structure of the EAFZ, which features a more complex fault geometry accompanied by numerous short normal faults in the southwest portion. These results shed light on fault segmentation, rupture speed variation, and slip partitioning along the EAFZ, advancing our understanding of fault structural control on earthquake behaviors in a complex multisegment fault system.

Supershear rupture with a two-scale damage model

An investigation of intersonic rupture is performed in a two-scale framework, with a damage model obtained by homogenization from microstructures with dynamically evolving mode II microcracks. The macroscopic fracture is analyzed in connection with the microcrack propagation. It is found that the supershear transition and the subsequent intersonic failure are achieved by collective effects of propagation and coalescence of sub-Rayleigh microcracks ahead the macroscopic rupture front. This is similar to the behavior experimentally observed for mode I dynamic brittle fracture (Guerra et al. 2012). Numerical simulations of a pre-cracked plate under shear loading are performed and characteristic features like the formation of shear shock waves, the mother-daughter or direct transition mechanisms, the specific evolutions of the rupture speed or the symmetry breaking near the crack tip are retrieved. Parametric studies allow to determine the influence of the size of the microstructure, the applied loading and the fracture properties of the solid. For the asymmetric impact of a Homalite plate problem (Rosakis et al., 1999), the numerical results obtained with the damage model are successfully confronted with the experimental observations.

Dynamics of episodic supershear in the 2023 M7.8 Kahramanmaraş/Pazarcik earthquake, revealed by near-field records and computational modeling

The 2023 M7.8 Kahramanmaraş/Pazarcik earthquake was larger and more destructive than what had been expected. Here we analyzed nearfield seismic records and developed a dynamic rupture model that reconciles different currently conflicting inversion results and reveals spatially non-uniform propagation speeds in this earthquake, with predominantly supershear speeds observed along the Narli fault and at the southwest (SW) end of the East Anatolian Fault (EAF). The model highlights the critical role of geometric complexity and heterogeneous frictional conditions in facilitating continued propagation and influencing rupture speed. We also constrained the conditions that allowed for the rupture to jump from the Narli fault to EAF and to generate the delayed backpropagating rupture towards the SW. Our findings have important implications for understanding earthquake hazards and guiding future response efforts and demonstrate the value of physics based dynamic modeling fused with near-field data in enhancing our understanding of earthquake mechanisms and improving risk assessment.

Tsunami potential in the Makran subduction zone: Amplification effects from earthquake rupture directivity and speed

Tsunami potential in the Makran subduction zone: Amplification effects from earthquake rupture directivity and speed

Kinematic Rupture Model of the 6 February 2023 Mw 7.8 Türkiye Earthquake from a Large Set of Near-Source Strong-Motion Records Combined with GNSS Offsets Reveals Intermittent Supershear Rupture

ABSTRACT The 2023 Mw 7.8 southeast Türkiye earthquake was recorded by an unprecedentedly large set of strong-motion stations very close to its rupture, opening the opportunity to observe the rupture process of a large earthquake with fine resolution. Here, the kinematics of the earthquake source are inferred by finite-source inversion based on strong-motion records and coseismic offsets from permanent Global Navigation Satellite Systems stations. The strong-motion records at stations NAR and 4615, which are the closest to the splay fault (SPF) where the rupture initiated and which were previously interpreted to contain the signature of supershear rupture speeds, are successfully modeled here by a subshear rupture propagating unilaterally to the northeast. Once the rupture on the SPF reaches the east Anatolian fault (EAF), it propagates on the EAF bilaterally, extending about 120 km northeast and 180 km southwest. To the south, the depth extent of the rupture decreases, as it passes a bend of the EAF. Although the rupture velocity remains globally subshear along the EAF, we identify three portions of the fault where the rupture is transiently supershear. The transitions to supershear speed coincide with regions of reduced fault slip, which suggests supershear bursts generated by the failure of local rupture barriers. Toward the southwest termination, the rupture encircles an asperity before its failure, which is a feature that has been observed only on rare occasions. This unprecedented detail of the inversion was facilitated by the proximity to the fault and the exceptional density of the accelerometric network in the area.

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.

Complex multi-fault rupture and triggering during the 2023 earthquake doublet in southeastern Türkiye

Two major earthquakes (MW 7.8 and MW 7.7) ruptured left-lateral strike-slip faults of the East Anatolian Fault Zone (EAFZ) on February 6, 2023, causing >59,000 fatalities and ~$119B in damage in southeastern Türkiye and northwestern Syria. Here we derived kinematic rupture models for the two events by inverting extensive seismic and geodetic observations using complex 5-6 segment fault models constrained by satellite observations and relocated aftershocks. The larger event nucleated on a splay fault, and then propagated bilaterally ~350 km along the main EAFZ strand. The rupture speed varied from 2.5-4.5 km/s, and peak slip was ~8.1 m. 9-h later, the second event ruptured ~160 km along the curved northern EAFZ strand, with early bilateral supershear rupture velocity (>4 km/s) followed by a slower rupture speed (~3 km/s). Coulomb Failure stress increase imparted by the first event indicates plausible triggering of the doublet aftershock, along with loading of neighboring faults.

New physics of supersonic ruptures

AbstractUntil recently, it is believed that the rupture speed above the pressure wave is impossible since spontaneously propagating ruptures are driven by the energy released due to the rupture motion, which is transferred through the medium to the rupture tip region at the maximum speed equal to the pressure wave speed. However, the apparent violation of classic theories has been revealed by new experimental results demonstrating supersonic shear ruptures. This paper presents a detailed analysis of the recently discovered shear rupture mechanism (fan hinged), which suggests a new physics of energy supply to the tip of supersonic ruptures. The key element of this mechanism is the fan‐shaped structure of the head of extreme ruptures, which is formed as a result of an intense tensile cracking process with the creation of intercrack slabs that act as hinges between the shearing rupture faces. The fan structure is featured with the following extraordinary properties: extremely low friction approaching zero; amplification of shear stresses above the material strength at low applied shear stresses; creation of a self‐disbalancing stress state causing a spontaneous rupture growth; abnormally high energy release; generation of driving energy directly at the rupture tip which excludes the need to transfer energy through the medium. The fan mechanism operates in intact rocks at stress conditions corresponding to seismogenic depths and in pre‐existing extremely smooth interfaces due to identical tensile cracking processes at these conditions. This is Paper 1 (of two companion papers) which discusses the fan theory and extreme ruptures in experiments on extremely smooth interfaces. Paper 2 entitled “Fan‐hinged shear instead of frictional stick‐slip as the main and most dangerous mechanism of natural, induced and volcanic earthquakes in the earth's crust” considers extreme ruptures in intact rocks. Further study of this subject is a major challenge for deep underground science, earthquake and fracture mechanics, physics, and tribology.

Ground-Motion Variability from Kinematic Rupture Models and the Implications for Nonergodic Probabilistic Seismic Hazard Analysis

Abstract The variability of earthquake ground motions has a strong control on probabilistic seismic hazard analysis (PSHA), particularly for the low frequencies of exceedance used for critical facilities. We use a crossed mixed-effects model to partition the variance components from simulated ground motions of Mw 7 earthquakes on the Salt Lake City segment of the Wasatch fault zone. Total variability of simulated ground motions is approximately equivalent to empirical models. The high contribution from rupture speed suggests an avenue to reducing variability through research on the causes and predictions of rupture speed on specific faults. Simulations show a strong spatial heterogeneity in the variability that manifests from directivity effects. We illustrate the impact of this spatial heterogeneity on hazard using a partially nonergodic PSHA framework. The results highlight the benefit of accounting for directivity effects in nonergodic PSHA, in which models that account for additional processes controlling ground motions are paired with reductions in the modeled ground-motion variability.

Dynamic Rupture Process of the 2023 Mw 7.8 Kahramanmaraş Earthquake (SE Türkiye): Variable Rupture Speed and Implications for Seismic Hazard

AbstractWe considered various non‐uniformities such as branch faults, rotation of stress field directions, and changes in tectonic environments to simulate the dynamic rupture process of the 6 February 2023 Mw 7.8 Kahramanmaraş earthquake in SE Türkiye. We utilized near‐fault waveform data, GNSS static displacements, and surface rupture to constrain the dynamic model. The results indicate that the high initial stress accumulated in the Kahramanmaraş‐Çelikhan seismic gap leads to the successful triggering of the East Anatolian Fault (EAF) and the supershear rupture in the northeast segment. Due to the complexity of fault geometry, the rupture speed along the southeastern segment of the EAF varied repeatedly between supershear and subshear, which contributed to the unexpectedly strong ground motion. Furthermore, the triggering of the EAF reminds us to be aware of the risk of seismic gaps on major faults being triggered by secondary faults, which is crucial to prevent significant disasters.

Study of Finite Element Simulation on the Mechano-Bactericidal Mechanism of Hierarchical Nanostructure Arrays.

Biomimetic nanostructures with bactericidal performance have become the research focus in constructing sterilization surfaces, but the mechano-bactericidal mechanism is still not fully understood, especially for the hierarchical nanostructure arrays with different heights. Herein, the interaction between Escherichia coli cells and nanostructure arrays was simulated by finite element, and the initial rupture points, i.e., critical action sites, of bacterial cells and the effects of nanostructure geometries on the cell rupture speed were analyzed based on the mechano-response of Escherichia coli cells on flat (identical heights) and hierarchical nanostructure arrays. The critical action sites of bacterial cells on nanostructure arrays are all at the three-phase junction zone of cell-liquid-nanostructure, but they are slightly shifted by the height difference ΔH of nanostructures on hierarchical nanopillar (NP)/nanosheet (NS) arrays, where the NP is higher than the NS. When ΔH < 20 nm, the site nears the NS corners, and when ΔH ≥ 20 nm, the site is consistent with that of the NP/NP array, i.e., the site locates at the three-phase junction zone of cell-liquid-high NP. In addition, except for decreasing the NP diameter, the NS thickness/width, or properly increasing the nanostructure spacing, the cell rupture can be accelerated via increasing the ΔH of nanostructures. ΔH = 40 nm is distinguished as the boundary for the effect of nanostructure ΔH on the cell rupture speed. When ΔH < 40 nm, the cell rupture speed rapidly increases as the ΔH increases; when ΔH ≥ 40 nm, the cell rupture speed reaches the maximum value and remains stable. This study provides a new strategy on how to design high-efficiency bactericidal surfaces.

How Does Thermal Pressurization of Pore Fluids Affect 3D Strike-Slip Earthquake Dynamics and Ground Motions?

ABSTRACT Frictional heating during earthquake rupture raises the fault-zone fluid pressure, which affects dynamic rupture and seismic radiation. Here, we investigate two key parameters governing thermal pressurization of pore fluids – hydraulic diffusivity and shear-zone half-width – and their effects on earthquake rupture dynamics, kinematic source properties, and ground motions. We conduct 3D strike-slip dynamic rupture simulations assuming a rate-and-state dependent friction law with strong velocity weakening coupled to thermal-pressurization of pore fluids. Dynamic rupture evolution and ground shaking are densely evaluated across the fault and Earth’s surface to analyze the variations of rupture parameters (slip, peak slip rate, rupture speed, and rise time), correlations among rupture parameters, and variability of peak ground velocity. Our simulations reveal how variations in thermal-pressurization affect earthquake rupture properties. We find that the mean slip and rise time decrease with increasing hydraulic diffusivity, whereas mean rupture speed and peak slip-rate remain almost constant. Mean slip, peak slip-rate, and rupture speed decrease with increasing shear-zone half-width, whereas mean rise time increases. Shear-zone half-width distinctly affects the correlation between rupture parameters, especially for parameter pairs (slip, rupture speed), (peak slip-rate, rupture speed), and (rupture speed, rise time). Hydraulic diffusivity has negligible effects on these correlations. Variations in shear-zone half-width primarily impact rupture speed, which then may affect other rupture parameters. We find a negative correlation between slip and peak slip-rate, unlike simpler dynamic rupture models. Mean peak ground velocities decrease faster with increasing shear-zone half-width than with increasing hydraulic diffusivity, whereas ground-motion variability is similarly affected by both the parameters. Our results show that shear-zone half-width affects rupture dynamics, kinematic rupture properties, and ground shaking more strongly than hydraulic diffusivity. We interpret the importance of shear-zone half-width based on the characteristic time of diffusion. Our findings may inform pseudodynamic rupture generators and guide future studies on how to account for thermal-pressurization effects.

Mechanism study on structural and mechanical properties of triblock copolymer hydrogels by coarse-grained molecular dynamics simulations

In this paper, the fracture mechanism of a hydrophobic-hydrophilic-hydrophobic triblock copolymer physical hydrogel system is investigated through customized coarse-grained molecular dynamics simulations. The stress-strain behavior of the system subject to uniaxial tensile stretch is examined at various strain rates. The stretching of hydrogel network can be described as follows: (1) for hydrophilic chains: entangled-Tightened-Alignment-Coiled; (2) for hydrophobic chains: separated sphere-Coalescence-Ellipsoid-Pulled out-Coiled. The coalescence speed and rupture speed are found to be two competing factors that dominate the rupture process of the physical hydrogel, and both are correlated with the strain rate parameter. When the strain rate is quite high, the hydrophobic chains will be pulled out quickly without a chance to coalesce or get entangled, resulting in inferior mechanical properties. The optimum ratio of hydrophobic segment and hydrophilic segment for the best gel mechanics is then studied, resulting that a segment ratio of 1:6:1 for the hydrophobic-hydrophilic-hydrophobic hydrogel that may lead to the strongest mechanical properties. In summary, we revealed that the architecture of triblock copolymer can be optimized, and thus the mechanical properties of these advanced polymer networks can be tuned for specific applications.

Laboratory earthquakes decipher control and stability of rupture speeds

Earthquakes are destructive natural hazards with damage capacity dictated by rupture speeds. Traditional dynamic rupture models predict that earthquake ruptures gradually accelerate to the Rayleigh wave speed with some of them further jumping to stable supershear speeds above the Eshelby speed (~sqrt{2} times S wave speed). However, the 2018 Mw 7.5 Palu earthquake, among several others, significantly challenges such a viewpoint. Here we generate spontaneous shear ruptures on laboratory faults to confirm that ruptures can indeed attain steady subRayleigh or supershear propagation speeds immediately following nucleation. A self-similar analysis of dynamic rupture confirms our observation, leading to a simple model where the rupture speed is uniquely dependent on a driving load. Our results reproduce and explain a number of enigmatic field observations on earthquake speeds, including the existence of stable subEshelby supershear ruptures, early onset of supershear ruptures, and the correlation between the rupture speed and the driving load.

Understanding the Rupture Kinematics and Slip Model of the 2021 Mw 7.4 Maduo Earthquake: A Bilateral Event on Bifurcating Faults

AbstractWe image the rupture process of the 2021 Mw 7.4 Maduo, Tibet earthquake using slowness‐enhanced back‐projection (BP) and joint finite fault inversion, which combines teleseismic broadband body waves, long‐period (166–333 s) seismic waves, and 3D ground displacements from radar satellites. The results reveal a left‐lateral strike‐slip rupture, propagating bilaterally on a 160 km long north‐dipping sub‐vertical fault system that bifurcates near its east end. About 80% of the total seismic moment occurs on the asperities shallower than 10 km, with a peak slip of 5.7 m. To simultaneously match the observed long‐period seismic waves and static displacements, potential deep slip is required, despite a tradeoff with the rigidity of the shallow crust. The deep slip existence, local crustal rigidity, and synthetic long‐period Earth response for Tibet earthquakes thus deserve further investigation. The WNW branch ruptures ∼75 km at ∼2.7 km/s, while the ESE branch ruptures ∼85 km at ∼3 km/s, though super‐shear rupture propagation possibly occurs during the ESE propagation from 12 to 20 s. Synthetic BP tests confirm overall sub‐shear rupture speeds and reveal a previously undocumented limitation caused by the signal interference between two bilateral branches. The stress analysis on the forks of the fault demonstrates that the pre‐compression inclination, rupture speed, and branching angle could explain the branching behavior on the eastern fork.

A High‐Order Accurate Summation‐By‐Parts Finite Difference Method for Fully‐Dynamic Earthquake Sequence Simulations Within Sedimentary Basins

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