Dynamic Rupture Process Research Articles

Detailed mapping of the surface rupture of the 12 February 2014 Yutian Ms7.3 earthquake, Altyn Tagh fault, Xinjiang, China

Large-scale detailed mapping plays a key role in revealing the rupture characteristics and mechanisms of strong earthquakes. Relatively few studies have been performed on the surface ruptures of large earthquakes in central and western Tibet due to its remote nature and high elevation. Based on high-resolution unmanned aerial vehicle (UAV) photography, we mapped the coseismic surface rupture of the 2014 Yutian Ms7.3 earthquake. Along the western Altyn Tagh fault system, the earthquake produced ≈37 km of surface rupture along the South Xor Kol fault (southern section S1), Xor Kol fault (central section S2) and Ashikule fault (northern section S3). Section S1 has a 16-km-long surface rupture with an average sinistral offset of 52±25 cm and a maximum offset of ~90 cm, while section S3 has a 14.2-km-long surface rupture with an average sinistral offset of 36±21 cm and a maximum offset of ~84 cm. A compilation of 5308 cracks yields an average crack width along the southern section of 85±71 cm and a maximum width of ~700 cm; the average width along the central section is 39±21 cm, and the maximum width is 243 cm; and the average width along the northern section is 61±44 cm with a maximum of ~340 cm. In addition, the average cumulative opening across rupture zone is 3.4±2.9 m along the southern section, with a maximum of ~17 m; 4.3±3.6 m along the central section, with a maximum of ~13 m; and 1.7±1.6 m along the northern section, with a maximum of ~6 m. Evidently, the average crack width and cumulative opening decrease towards bends and steps along the fault. A global synthesis of surface rupture distributions corresponding to strike-slip earthquakes indicates that the rupture zone is wider near the complex parts of fault geometries (such as bends, steps and fault bifurcations) than along straight sections, suggesting that the fault geometry has an obvious control on the surface rupture width. The widespread cracks at the intersection between the Xor Kol and South Xor Kol faults may indicate that an extensional regime is more likely to produce distributed offfault deformation, which provides an observational constraint for the numerical simulation of dynamic rupture on a fault. In addition to coseismic surface rupture, the Yutian earthquake also produced a large number of gravity-driven slides on alluvial fans with gentle slopes. The friction efficiency of the water-bearing salt layer beneath fans could decrease the sliding threshold and trigger instability under surface shaking. These distributed deformations and gravity-driven slides reflect the coupling between the rupture propagation and fault geometry and indicate that the rupture may have propagated in two directions along the Ashikule fault after passing through a step. Therefore, the investigation of coseismic surface rupture provides important observational constraints on the dynamic rupture process.

Resolving the Kinematics and Moment Release of Early Afterslip Within the First Hours Following the 2016 Mw 7.1 Kumamoto Earthquake: Implications for the Shallow Slip Deficit and Frictional Behavior of Aseismic Creep

AbstractAs stresses following rupture are dissipated continuous measurements of postseismic surface deformation provide insight into variations of the frictional strength of faults and the rheology of the lower crust and upper mantle. Due to the difficulty of capturing the earliest phase of afterslip, most analyses have focused on understanding postseismic processes over timescales of weeks to years. Here we investigate the kinematics, moment release, and frictional properties of the earliest phase of afterslip within the first hours following the 2016 Mw 7.1 Kumamoto earthquake using a network of 5‐minute sampled continuous Global Positioning System (GPS) stations. Using independent component analysis to filter the GPS data, we find that (1) early afterslip contributes only ~1% of total moment release within the first hour and 8% after 24 hr. This suggests that the lack of a coseismic slip deficit, which we estimate using standard geodetic data sets (e.g., InSAR, GPS, and pixel offsets) and which span the first 4 days of the postseismic period, is largely reflective of the dynamic rupture process and we can rule out contamination of moment release by early afterslip. (2) Early afterslip shows no evidence of a delayed nucleation or acceleration phase, where instead fault patches transition to immediate deceleration following rupture that is consistent with frictional relaxation under steady state conditions with dependence only on the sliding velocity. (3) There is a close correlation between the near‐field aftershocks and afterslip within the first hours following rupture, suggesting afterslip may still be an important possible triggering mechanism during the earliest postseismic period.

Space-time adaptive ADER discontinuous Galerkin schemes for nonlinear hyperelasticity with material failure

We are concerned with the numerical solution of a unified first order hyperbolic formulation of continuum mechanics that goes back to the work of Godunov, Peshkov and Romenski [65,68,97] (GPR model) and which is an extension of nonlinear hyperelasticity that is able to describe simultaneously nonlinear elasto-plastic solids at large strain, as well as viscous and ideal fluids. The proposed governing PDE system also contains the effect of heat conduction and can be shown to be symmetric and thermodynamically compatible, as it obeys the first and second law of thermodynamics. In this paper we extend the GPR model to the simulation of nonlinear dynamic rupture processes, which can be achieved by adding an additional scalar to the governing PDE system. This extra parameter describes the material damage and is governed by an advection-reaction equation, where the stiff and highly nonlinear reaction mechanisms depend on the ratio of the local equivalent stress to the yield stress of the material. The stiff reaction mechanisms are integrated in time via an efficient exponential time integrator. Due to the multiple spatial and temporal scales involved in the problem of crack generation and propagation, the model is solved on space–time adaptive Cartesian meshes using high order accurate discontinuous Galerkin finite element schemes endowed with an a posteriori subcell finite volume limiter. A key feature of our new model is the use of a twofold diffuse interface approach that allows the cracks to form anywhere and at any time, independently of the chosen computational grid, which is simply adaptive Cartesian (AMR). This is substantially different from many fracture modeling approaches that need to resolve discontinuities explicitly, such as for example dynamic shear rupture models used in computational seismology, where the geometry of the rupture fault needs to prescribed a priori. We furthermore make use of a scalar volume fraction function α that indicates whether a given spatial point is inside the solid (α=1) or outside (α=0), thus allowing the description of solids of arbitrarily complex shape. We show extensive numerical comparisons with experimental results for stress-strain diagrams of different real materials and for the generation and propagation of fracture in rocks and pyrex glass at low and high velocities. Overall, a very good agreement between numerical simulations and experiments is obtained. The proposed model is also naturally able to describe material fatigue.

The Near‐Fault Ground Motion Characteristics of Sustained and Unsustained Free Surface‐Induced Supershear Rupture on Strike‐Slip Faults

AbstractThe free surface is shown to be one of the key factors that may promote supershear rupture propagation on strike‐slip faults even if its initial shear stress is not larger enough as predicted by the Burridge‐Andrews mechanism. However, previous study has shown the free surface‐induced supershear rupture may be unsustained, which turns to sub‐Rayleigh rupture itself as the rupture propagates. We study the near‐field ground motion of sustained and unsustained supershear ruptures using the finite difference method based on the dynamic rupture processes of three vertical strike‐slip faults with the same initial stresses and different hypocenter depths. Both the unsustained supershear ruptures with shallower hypocenter depths show the sub‐Rayleigh characteristics in the peak ground velocity distribution, that is, strong amplitude is noticed beyond the end of the fault without any observed Mach cone. We observe that the arrival time differences between the maximum fault‐perpendicular and fault‐parallel velocity reveal the Mach cone clearly for the sustained supershear rupture. A distinct supershear phase in the high‐frequency seismograms is observed. We also compare the 70° dip strike‐slip case in which unsustained supershear rupture may still present sub‐Rayleigh characteristics in the near‐field ground motion, and the break of symmetry perturbs the normal stress as rupture propagates, which has a key influence in the rupture propagation and ground motion. Our work provides a new insight to understand the supershear rupture in nature earthquakes and the relationship between the ground motion and the rupture process on the fault.

Dynamic models of earthquake rupture along branch faults of the eastern San Gorgonio Pass region in California using complex fault structure

AbstractGeologic data suggest that the Coachella Valley segment of the southern San Andreas fault (southern California, USA) is past its average recurrence time period. At its northern edge, this right-lateral fault segment branches into the Mission Creek and Banning strands of the San Andreas fault. Depending on how rupture propagates through this region, there is the possibility of a throughgoing rupture that could lead to the channeling of damaging seismic energy into the Los Angeles Basin. The fault structures and potential rupture scenarios on these two strands differ significantly, which highlights the need to determine which strand provides a more likely rupture path and the circumstances that control this rupture path. In this study, we examine the effect of different assumptions about fault geometry and initial stress pattern on the dynamic rupture process to test multiple rupture scenarios and thus investigate the most likely path(s) of a rupture that starts on the Coachella Valley segment. We consider three types of fault geometry based on the Southern California Earthquake Center Community Fault Model, and we create a three-dimensional finite-element mesh for each of them. These three meshes are then incorporated into the finite-element method code FaultMod to compute a physical model for the rupture dynamics. We use a slip-weakening friction law, and consider different assumptions of background stress, such as constant tractions and regional stress regimes with different orientations. Both the constant and regional stress distributions show that rupture from the Coachella Valley segment is more likely to branch to the Mission Creek than to the Banning fault strand. The fault connectivity at this branch system seems to have a significant impact on the likelihood of a throughgoing rupture, with potentially significant impacts for ground motion and seismic hazard both locally and in the greater Los Angeles metropolitan area.

Analysis of the role of branching angle in the dynamic rupture process on a 3-D branching fault system

The fault branching phenomenon, which may heavily influence the patterns of rupture propagation in fault systems, is one of the geometric complexities of fault systems that is widely observed in nature. In this study, we investigate the effect of the branching angle on the rupture inclination and the interaction between branch planes in two-fork branching fault systems by numerical simulation and theoretical analysis based on Mohr’s circle. A friction law dependent on normal stress is used, and special attention is paid to studying how ruptures on the upper and lower branch planes affect the stress and rupture on each other separately. The results show that the two branch planes affect each other in different patterns and that the intensity of the effect changes with the branching angle. The rupture of the lower branch plane has a negative effect on the rupture of the upper branch plane in the case of a small branching angle but has almost no negative effect in the case of a large branching angle. The rupture of the upper branch plane, however, suppresses the rupture of the lower branch plane regardless of whether the branching angle is large or small.

Analysis of rupture mode and acoustic emission characteristic of rock and coal samples with holes

Abstract In order to understand the rupture mechanism of rock and coal samples with holes, the acoustic emission (AE) tests of rock and coal samples with holes under uniaxial compression are done. Through the AE tests, the basic mechanical and AE variation rule of two samples in the total rupture process are obtained, the dynamic rupture process was observed and finally the spatial evolution and multi-fractal characteristic of AE are analyzed. The results show that the variation rule of AE events of two samples are coincident with the variation rule of stress. The uniaxial compressive strength of a rock sample with a hole is about five times than that of a coal sample with a hole, and the maximum AE pulsing counts (AEPC) of a rock sample with a hole is much larger than that of a coal sample with a hole. Due to the variation of lithology, the dynamic rupture process and the spatial evolution rule of AE events of rock and coal samples with holes are obviously different. But the distributions of the three-dimensional locations of the two samples are coincident with the macroscopic rupture morphology. Both of them have multi-fractal characteristic in the process of rupture, the multi-fractal spectrum width Δf(α) at the time of rupture is smaller than that before rupture and the Δf(α) before rupture is smaller than that after rupture. This indicates that the energy (E) before and after rupture is smaller than that at the time of rupture, the E after rupture is smaller than that before rupture and the E of a coal sample with a hole in each stage is less than that of the corresponding stages of a rock sample with a hole.

Analysis of progressive rupture process in surrounding rock for a deep homogeneous and circular opening

ABSTRACT The failure mode and scope of rock surrounding a rupture zone have important theoretical value for stability analysis and support design of excavated openings. This study uses the Mohr–Coulomb yield criteria to study the general pattern of rupture progress in rock surrounding circular openings by resolving the analytical solution of stress in elasto-plastic mechanics of the circular openings. The results indicate a dynamic rupture process that develops gradually in homogeneous and circular surrounding rock in the mode of circular annuli under conditions of hydrostatic pressure. The concept of a sub-rupture zone in the surrounding rock is also proposed, and the criterion for rupture of the sub-rupture zone is derived. The calculation of radial stress is undertaken for the outer boundary of the sub-rupture zone, and the calculation method of radii in the sub-rupture zone is obtained. The equality of radii of adjacent sub-rupture zones is taken as the criterion of cessation of rupturing in a rupture zone, and is the criterion upon which the analytical solutions of maximum radius and maximum radial stress are deduced based. The analytical solutions are highly consistent with the engineering example and numerical simulation results using FLAC3D™, which reveals the failure mechanism of surrounding rock and validates the accuracy of the theoretical calculations. The results provide a theoretical basis for optimizing support design and safety assessment in underground engineering projects.

Imaging of Seismogenic Asperities of the 2016 ML 6.0 Amatrice, Central Italy, Earthquake Through Dynamic Rupture Simulations

Numerical simulations are carried out for the dynamic rupture and wave propagation process of the 24th August 2016 ML 6.0 Amatrice, Italy, earthquake, using a boundary domain method (BDM), a hybrid method of boundary integral equation and finite difference methods. Dynamic rupture parameters of two seismogenic asperities are searched by iterative search through the comparison of the near-field ground motions. The preferred models indicate two asperities, aligned at around 4–5 km depth and separated from each other as well as from the initial rupture point. This requires a few supplementary patches connecting them, and that are less energetic than the asperities. The asperities are characterized by a radius of 2–3 km in the south (first to rupture) and of 2–4 km in the north (second to rupture), and the corresponding fracture energies of the asperities are (25.35 ± 0.63) × 1012 J and (38.05 ± 7.91) × 1012 J, respectively. These values are consistent with the scaling relation extrapolated from various analyses of large earthquakes. Although the parameter space of the search is limited due to the numerical performance of the dynamic rupture simulation, the proposed simple characterization of the earthquakes source confirms the scaling relation in fracture energy of the seismogenic asperities, which is essential for constructing mechanical earthquake source models.

Tsunami hazard assessment for atoll islands inside the South China Sea: A case study of the Xisha Archipelago

The potential tsunami hazard for atoll islands inside the South China Sea (SCS) has not been investigated previously. In this study, we assess the tsunami hazard for the Xisha Archipelago using a deterministic approach. The tsunami source considered in this study is hypothesized earthquakes in the Manila subduction zone. We employ COMCOT (COrnell Multi-grid Coupled Tsunami model) for tsunami simulations. Different with previous deterministic assessments, our study investigates the uncertainties caused by bathymetric resolution and rupture complexities (here refers to rupture speed and direction variations related to fault geometry complexity). Comparison between results of different resolution suggests that coarser grids may over/underestimate tsunami hazards around mid-ocean atoll islands and this depends on the surrounding bathymetry. But results of different datasets show that GEBCO may lead to a lower result than the multibeam data. We also investigate the uncertainties caused by different rupture speed and direction, which are highly affected by subducting seamounts/ridge. Different rupture direction will shift the hardest hit locations to areas, which are possibly underestimated in cases of instantaneous ruptures. Also, the distribution of slip and location of the rupture zone can affect the local tsunami wave height significantly. We suggest that near-island high-resolution bathymetric data and dynamic rupture process should be considered in future tsunami hazard assessment in locations with similar geological structures.

Theoretical and numerical investigations of floor dynamic rupture: A case study in Zhaolou Coal Mine, China

The floor dynamic rupture occurred in the 1306 longwall coal face (LCF) in Zhaolou Coal Mine was induced due to the superposition of static and dynamic stresses. In order to investigate the triggering conditions of the floor buckling, this study established a cusp catastrophe model based on elastic thin plate theory. During the occurrence of floor dynamic rupture process, the released energy was calculated, using elasticity modulus (E), Poisson ratio (μ), thickness (h), roadway width (a), maximum cohesion (σ0) and horizontal stress (σh) as influencing factors. Furthermore, by utilizing discrete element method (UDEC), and applying an artificial source in the model, a fish function was developed to harness the released kinetic energy (Ek) and its vibration curves. The results indicated that the main control factors of the floor dynamic rupture were the roadway width and horizontal stress, both of which had a positive correlation with the floor dynamic rupture risk. According to seismic computed tomography, the high horizontal stress in the 1306 LCF was verified through the evolution of P-wave velocity. The above mentioned findings, coupled with the catastrophe theory, UDEC and seismic computed tomography, could propose a certain reference for the pre-warning and prevention of the floor dynamic rupture.

Dynamic rupture of metal sheet subjected to laser irradiation and tangential subsonic airflow

ABSTRACT To reproduce the premature rupture process of metal sheet subjected to laser irradiation with subsonic airflow, which is an interesting phenomenon observed in the experiments given by Lawrence Livermore National Laboratory, a coupled numerical model considering the interaction and evolution of metal elastoplastic deformation and aerodynamic pressure profile is presented. With the thermal elastoplastic constitutive relationship and failure criterion, the simulated failure modes and dynamic rupture process are basically consistent with the experimental results, indicating plastic flow and multiple fracturing is the main failure mechanism. Compared with the case of non-airflow, subsonic airflow not only accelerates deformation, but also turns the bugle deformation, plastic strain and rupture mode into asymmetric.

Estimation and applicability of attenuation characteristics for source parameters and scaling relations in the Garhwal Kumaun Himalaya region, India

Source parameters of the small to moderate earthquakes are significant for understanding the dynamic rupture process, the scaling relations of the earthquakes and for assessment of seismic hazard potential of a region. In this study, the source parameters were determined for 58 small to moderate size earthquakes (3.0 ≤ Mw ≤ 5.0) occurred during 2007–2015 in the Garhwal-Kumaun region. The estimated shear wave quality factor (Qβ(f)) values for each station at different frequencies have been applied to eliminate any bias in the determination of source parameters. The Qβ(f) values have been estimated by using coda wave normalization method in the frequency range 1.5–16 Hz. A frequency-dependent S wave quality factor relation is obtained as Qβ(f) = (152.9 ± 7) f(0.82±0.005) by fitting a power-law frequency dependence model for the estimated values over the whole study region. The spectral (low-frequency spectral level and corner frequency) and source (static stress drop, seismic moment, apparent stress and radiated energy) parameters are obtained assuming ω−2 source model. The displacement spectra are corrected for estimated frequency-dependent attenuation, site effect using spectral decay parameter “Kappa”. The frequency resolution limit was resolved by quantifying the bias in corner frequencies, stress drop and radiated energy estimates due to finite-bandwidth effect. The data of the region shows shallow focused earthquakes with low stress drop. The estimation of Zúñiga parameter (ε) suggests the partial stress drop mechanism in the region. The observed low stress drop and apparent stress can be explained by partial stress drop and low effective stress model. Presence of subsurface fluid at seismogenic depth certainly manipulates the dynamics of the region. However, the limited event selection may strongly bias the scaling relation even after taking as much as possible precaution in considering effects of finite bandwidth, attenuation and site corrections. Although, the scaling can be improved further with the integration of large dataset of microearthquakes and use of a stable and robust approach.

​Demulsification performance, behavior and mechanism of different demulsifiers on the light crude oil emulsions

Chemical demulsifiers were widely applied for demulsification of crude oil emulsions which were generally formed with natural stabilizers during the extraction. Based on our previous study of the light crude oil emulsions stability, the demulsification performance of four classes of nonionic poly(ethylene oxide) poly(propylene oxide) copolymer demulsifiers and sodium dodecyl sulfate (SDS) as well as their dynamic demulsification process on the W/O light crude oil emulsion were analyzed by the multiple light scattering method. The demulsification mechanism of these demulsifiers were illuminated by the measurement of dynamic interfacial tension (IFT), dynamic interfacial viscoelasticity, rupture rate and microscopic changes of oil film by single droplet method. The results have shown that demulsification performance doesn’t depend on IFT, though effective demulsifiers could have a lower IFT around 1 mN/m. Multiple light scattering data and dynamic IFT clarified the demulsification process of efficient demulsifiers, quick diffusion and adsorption of demulsifiers, palpable change in interfacial property, droplet coalescence and sedimentation. The capacities of demulsifiers to penetrate into the interface and to decrease the interfacial viscoelasticity were the dominating demulsification mechanisms. The quantitative interpretation of multiple light scattering data associated with the dynamic IFT and interfacial viscoelasticity illuminated the demulsification mechanism of different demulsifiers. The rupture rate and dynamic rupture process of the oil film between the aqueous phases also vividly showed the differences in demulsification process. The nonionic demulsifiers AR902 and PR929 with aromatic groups has the highest speed and best performance due to the similar aromatic groups with those in asphaltenes, which produced the least stable oil film. This study provides a valuable guideline for the selection of demulsifier for light crude oil emulsions and enhanced demulsification mechanism studies.

Persistent rupture terminations at a restraining bend from slip rates on the eastern Altyn Tagh fault

Restraining double-bends along strike-slip faults inhibit or permit throughgoing ruptures depending on bend angle, length, and prior rupture history. Modeling predicts that for mature strike-slip faults in a regional stress regime characterized by simple shear, a restraining bend of >18° and >4 km length impedes propagating rupture. Indeed, natural evidence shows that the most recent rupture(s) of the Xorkoli section (90°–93°E) of the eastern Altyn Tagh fault (ATF) ended at large restraining bends. However, when multiple seismic cycles are considered in numerical dynamic rupture modeling, heterogeneous residual stresses enable some ruptures to propagate further, modulating whether the bends persistently serve as barriers. These models remain to be tested using observations of the cumulative effects of multiple earthquake ruptures. Here we investigate whether a large restraining double-bend on the ATF serves consistently as a barrier to rupture by measuring long-term slip rates around the terminus of its most recent surface rupture at the Aksay bend. Our results show a W—E decline in slip as the SATF enters the bend, as would be predicted from repeated rupture terminations there. Prior work demonstrated no Holocene slip on the central, most misoriented portion of the bend, while 19–79 m offsets suggest that multiple ruptures have occurred on the west side of the bend during the Holocene. Thus we conclude the gradient in the SATF's slip rate results from the repeated termination of earthquake ruptures there. However, a finite slip rate east of the bend represents the transmission of some slip, suggesting that a small fraction of ruptures may fully traverse or jump the double-bend. This agreement between natural observations of slip accumulation and multi-cycle models of fault rupture enables us to translate observed slip rates into insight about the dynamic rupture process of individual earthquakes as they encounter geometric complexities along faults.

Investigation into Mechanism of Floor Dynamic Rupture by Evolution Characteristics of Stress and Mine Tremors: A Case Study in Guojiahe Coal Mine, China

In order to explore the mechanism of floor dynamic rupture, the current study adopts a thin plate model to further investigate the condition of floor failure. One of the possible explanations could be floor buckling due to high horizontal stress and dynamic disturbance ultimately leading to rapid and massive release of elastic energy thus inducing dynamic rupture. Seismic computed tomography and 3D location were employed to explore the evolution characteristics of floor stress distribution and positions of mine tremors. In the regions of floor dynamic rupture, higher P‐wave velocity was recorded prior to the dynamic rupture. On the contrary, relatively lower reading was observed after the dynamic rupture thus depicting a high stress concentration condition. Meanwhile, evolution of mine tremors revealed the accumulation and subsequent release of energy during the dynamic rupture process. It was further revealed that dynamic rupture was induced due to the superposition of static and dynamic stresses: (i) the high static stress concentration due to frontal and lateral abutment stress from coal pillar and (ii) dynamic stress from the fracture and caving of coal pillar, hard roof, and key stratum. In the later part of this study, the floor dynamic rupture occurrence process would be reproduced through numerical simulations within a 0.6 sec time frame. The above‐mentioned findings would be used to propose a feasible mechanism for prewarning and prevention of floor dynamic rupture using seismic computed tomography and mine tremors 3D location.

Effect of Seismogenic Depth and Background Stress on Physical Limits of Earthquake Rupture Across Fault Step Overs

AbstractEarthquakes can rupture geometrically complex fault systems by breaching fault step overs. Quantifying the likelihood of rupture jump across step overs is important to evaluate earthquake hazard and to understand the interactions between dynamic rupture and fault growth processes. Here we investigate the role of seismogenic depth and background stress on physical limits of earthquake rupture across fault step overs. Our computational and theoretical study is focused on the canonical case of two parallel strike‐slip faults with large aspect ratio, uniform prestress and friction properties. We conduct a systematic set of 3‐D dynamic rupture simulations with different seismogenic depth, step over distance, and initial stresses. We find that the maximum step over distance Hc that a rupture can jump depends on seismogenic depth W and strength excess to stress drop ratio S, commonly used to evaluate probable rupture velocity, as , where n = 2 when Hc/W < 0.2 (or S > 1.5) and n = 1 otherwise. The critical nucleation size, largely controlled by frictional properties, has a second‐order effect on Hc. Rupture on the secondary fault is mainly triggered by the stopping phase emanated from the rupture end on the primary fault. Asymptotic analysis of the peak amplitude of stopping phases sheds light on the mechanical origin of the relations between Hc, W, and S, and leads to the scaling regime with n = 1 in far field and n = 2 in near field. The results suggest that strike‐slip earthquakes on faults with large seismogenic depth or operating at high shear stresses can jump wider step overs than observed so far in continental interplate earthquakes.

Dynamic rupture propagation on geometrically complex fault with along-strike variation of fault maturity: insights from the 2014 Northern Nagano earthquake

Understanding the effect of the complex fault geometry on the dynamic rupture process and discriminating it from the complexity originating from the rheological properties of faults, is an essential subject in earthquake science. The 2014 Northern Nagano earthquake, which occurred near the end zone of a major active fault system, provided unique observations that enabled us to investigate both the detailed geometrical fault structure and surface deformation patterns as well as the temporal sequence led up from a prominent foreshock activity. We first develop a geometrical fault model with a substantial variation along strike, and a model for the regional stress field, which is well constrained by the observations. This significant along-strike variation in fault geometry probably reflects the difference of fault maturity at the end zone of the complex fault system. We used this model in order to conduct a set of dynamic rupture simulations using the highly efficient spatiotemporal boundary integral equation method. Based on our simulations, we show that the observed surface deformation can be reasonably explained as the effect of the non-planar fault geometry with a number of branch faults and bends.Graphical abstract.

Accurate estimation of seismic source parameters of induced seismicity by a combined approach of generalized inversion and genetic algorithm: Application to The Geysers geothermal area, California

AbstractThe accurate determination of stress drop, seismic efficiency, and how source parameters scale with earthquake size is an important issue for seismic hazard assessment of induced seismicity. We propose an improved nonparametric, data‐driven strategy suitable for monitoring induced seismicity, which combines the generalized inversion technique together with genetic algorithms. In the first step of the analysis the generalized inversion technique allows for an effective correction of waveforms for attenuation and site contributions. Then, the retrieved source spectra are inverted by a nonlinear sensitivity‐driven inversion scheme that allows accurate estimation of source parameters. We therefore investigate the earthquake source characteristics of 633 induced earthquakes (Mw 2–3.8) recorded at The Geysers geothermal field (California) by a dense seismic network (i.e., 32 stations, more than 17.000 velocity records). We find a nonself‐similar behavior, empirical source spectra that require an ωγ source model with γ > 2 to be well fit and small radiation efficiency ηSW. All these findings suggest different dynamic rupture processes for smaller and larger earthquakes and that the proportion of high‐frequency energy radiation and the amount of energy required to overcome the friction or for the creation of new fractures surface changes with earthquake size. Furthermore, we observe also two distinct families of events with peculiar source parameters that in one case suggests the reactivation of deep structures linked to the regional tectonics, while in the other supports the idea of an important role of steeply dipping faults in the fluid pressure diffusion.

Dynamic Rupture Modelling of the 1999 Düzce, Turkey Earthquake

The dynamic rupture process and near-source ground motion of the 1999 Mw 7.1 Du¨zce Earthquake are simulated. The fault rupture is governed by the slip-weakening friction model coupled to a three-dimensional viscoelastic wave equation. The problem is solved numerically by a 3-D dynamic rupture code that uses a generalized finite difference method. Initial parameterization of stress drop (∆τ) and strength excess (Se) for dynamic rupture calculations is obtained from the slip velocity distribution of a kinematic waveform inversion (KI) model by solving the elastodynamic equation with the kinematic slip as a boundary condition. Using the kinematic slip distribution and observed ground motion as constraints, a trial and error procedure was followed to define the stress parameterization. Preferred model describes the source in terms of stress with three asperities (located, respectively, at the deep, middle and shallow) and strong barriers between asperities. Se is as high as 19 Mpa at barriers between the three asperities and ∆τ is maximum about 40 Mpa at the deepest asperity. This heterogeneity in stress distribution produces abrupt jumps in rupture velocity, exhibiting locally apparent rupture speed exceeding the P wave velocity at the borders between barriers and asperities, due to sharp changes of fault strength and stress drop at those areas. Overall, consistent with other studies, the rupture propagation is dominated by supershear speed toward the eastern asperities and at shallow surface. Simulated surface rupture at the eastern fault is consistent with other studies; nevertheless, the western shallower parts did not rupture during the simulation, suggesting that those regions may have already broken during the 1999 Kocaeli event, which occurred three months earlier. Ground motion simulation catches the major characteristics of the observed waveforms. Distribution of simulated peak ground velocity (PGV) in low frequency (0.1–0.5 Hz.) inside the study area reveals the propagation pattern on the field, with PGV reaching to 1.2 and 2.2 m/s in the NS and EW components, respectively.