Coulomb Failure Stress Change Research Articles

Surface Displacement and Source Model Separation of the Two Strongest Earthquakes During the 2019 Ridgecrest Sequence: Insights From InSAR, GPS, and Optical Data

AbstractSince the occurrence of the 2019 Mw 6.4 and Mw 7.1 Ridgecrest earthquake sequence in the Eastern California Shear Zone, coseismic deformation following the two earthquakes has been intensively studied, and source modeling with interferometric synthetic aperture radar (InSAR), Global Positioning System (GPS) and seismological datasets has been favored. However, we recently found that the coseismic modeling of the two earthquakes constrained by the dense near‐field Planet‐Lab optical measurements can be more detailed than the previously estimated, which requires an accurate Planet‐Lab displacement. In this study, we improve the long wavelength orbital ramps correction method for obtaining a more accurate Planet‐Lab horizontal displacement field. The corrected dense near‐field Planet‐Lab data are firstly used together with the intermediate‐field GPS data to invert and distinguish the fault slip distribution for the two earthquakes. The same scheme is used to simultaneously invert the InSAR, SAR, GPS, and optical datasets, to improve the constraints on seismic source parameters. Our inversion results show that the joint‐event slip model is rougher than the two single‐event slip models, but it has a more concentrated slip pattern and larger slip amplitude in some zones. We show that adding a near‐field constraint of the Planet‐Lab data in the combined‐data inversion can reduce the slip parameter uncertainty and enhance the model resolution. The Coulomb failure stress changes on the southeastern Blackwater, the southern Owens Valley, and the central Panamint Valley faults are enhanced by about 0.4–0.8 bar by the 2019 Ridgecrest earthquake sequence.

San Andreas Fault Stress Change Due To Groundwater Withdrawal in California's Central Valley, 1860‐2010

Changes in groundwater storage in California's Central Valley (CV) are considered partly responsible for vertical uplift surrounding the southern CV and for stress changes on nearby faults. Questions remain regarding the distribution of stress on the central San Andreas fault (CSAF) from recent and historical drawdown of the CV aquifer over the past 150 years (1860–2010). We combine groundwater storage change estimates for the 2006–2010 drought with a three-dimensional finite element model to estimate Coulomb failure stress change (ΔCFS) on the CSAF. We combine a simple parameterization of historical hydraulic head change for 1860–1960 with a CV hydrological model (1961–2003) to estimate visco-elastic effects on ΔCFS total and 2010 rate. We find that ΔCFS on the CSAF correlates positively with shallow seismicity and low frequency earthquakes for both short-term and long term stressing. Unloading uplift rates surrounding the southern CV only partially account for the observed vertical GPS rates.

Investigating Stress Transfer Between the Tuz Gölü Fault Zone and Hasan Dağ Volcano (Turkey)

Faulting, magmatism and volcanism are intrinsically linked by plate tectonics. Fault slip imparts stress changes to the surrounding crust and other faults and fractures. Volcano-tectonic hazard assessment in areas with long recurrence intervals of volcanic and tectonic activity requires an assessment of current stress levels. Here we investigate stress transfer between the Akhisar-Kiliç fault segment (AKFS) of the seismically active Tuz Gölü fault zone in the Central Anatolian Volcanic Province and the active Hasan Dağ volcanic complex. Current stress accumulation by protracted gradual slip on the AKFS since its most recent rupture (5.45 ka±0.16 BP) is quantified using the Coulomb Failure Stress change (ΔCFS). We calculate currently accrued ΔCFS magnitudes of between 2.5±0.2 and 15±0.5 MPa on the fault plane for published lower and upper-bound estimates of right-lateral slip rates, respectively. These changes are sufficient to promote failure of the segment. The M5.1 September 20, 2020 earthquake SW of Hasan Dağ occurred in a volume predicted by this study to have undergone fault unclamping by gradual slip along the AKFS. We also show that gradual slip of the AKFS contributes to the progressive unclamping of fractures and transtensional opening of potential magma pathways oriented perpendicularly to the AKFS both above and below Hasan Dağ’s magma reservoir, while pathways oriented parallel to the AKFS are being clamped. Earthquake moment magnitudes of between Mw5.94 and Mw6.76 due to hypothetical partial or complete rupture of the segment exacerbate these trends. The spatial pattern of Coulomb failure stress changes on Hasan Dağ’s magma reservoir is predominantly controlled by the location of rupture relative to the reservoir with the magnitude of the earthquake playing a subordinate role. We explore implications of our findings for the assessment of interconnected seismo-volcanic hazards and associated risks.

Coseismic and postseismic deformation of the 2016 Mw 6.0 Petermann ranges earthquake from satellite radar observations

Continental earthquakes, especially shallow earthquakes, can cause disasters in populated areas. On 20 May 2016, a moderate magnitude earthquake (Mw 6.0) ruptured on Northern Territory of Australia. With the aim of seismic hazard assessment, coseismic and postsesmic deformation of the Petermann earthquake has been determined using ALOS-2 ascending and Sentinel-1A descending images. The coseismic deformation field suggests that the earthquake ruptured in the NNW direction with a length of about 20 km, a rake angle of 57.1° and a moment magnitude of Mw 6.06. The maximum slip of 1.31 m can be observed at the depth of 1.25 km. The small baseline subset (SBAS) InSAR was utilized to obtain the postseismic deformation from 69 Sentinel-1A images in a descending orbit. The spatiotemporal evolution of the postseismic deformation presents a gradually decaying trend, with a maximum displacement of 5.4 cm within 858 days after the mainshock. Our afterslip model suggests that the maximum slip of 0.28 m occurred at the depth of 0.25 km with a slight shift to the northwest compared to the coseismic slip. The static Coulomb failure stress (CFS) changes at various depths were calculated, indicating the consistency between high loaded stress zones and aftershock distribution.

Dynamic Rupture Simulations of the 2008 7.9 Wenchuan Earthquake: Implication for Heterogeneous Initial Stress and Complex Multifault Geometry

AbstractThe 2008 Mw7.9 Wenchuan earthquake ruptured northwest‐dipping imbricate oblique reverse faults along the Longmenshan thrust belt at the eastern margin of the Tibetan Plateau, and developed one of the most complex coseismic rupture patterns for reverse faulting events ever reported in intraplate settings. To evaluate to what extent complex fault geometry and heterogeneities in initial stress can explain the characteristics of coseismic observations, we simulate spontaneous dynamic rupture propagation governed by slip‐weakening friction law on the geometrically complex multifault system loaded by heterogeneous stress. Our model reproduces many observed features including the multifault ruptures, locations of the maximum slip, Interferometric Synthetic Aperture Radar observations, source time functions and peak ground velocities. Our results suggest that the maximum horizontal principal stress in the northern part of Beichuan Fault (BFC) is rotated approximately counterclockwise compared to that in the southern part of BFC, yet this rotation is not sufficient to produce a supershear rupture. We infer that the fault core zone in the southwestern Beichuan Fault (BCF) may have been more severely damaged than the northeastern counterpart, which provides an explanation that the rupture in the northern BCF propagates faster than in the southern portion. Besides, we find that Coulomb failure stress changes on Wenchuan‐Maoxian Fault (WMF) may be counteracted by the aftershock slips of the Lixian Fault, which suggests Wenchuan earthquake could not significantly push the WMF closer to failure.

Spatiotemporal Evolution of Pore Pressure Changes and Coulomb Failure Stress in a Poroelastic Medium for Different Faulting Regimes

AbstractSpatial distribution of aftershocks has been explained by the co‐seismic static Coulomb Failure Stress changes (ΔCFS) to some extent. In practice, the earth's crust is porous medium saturated with fluid filling pores, cavities, cracks, and faults. The pore fluid flow plays an important role in the evolution of the stress field in porous medium of the earth. The ΔCFS can be influenced by the co‐seismic pore pressure changes and the post‐seismic pore fluid flow. In order to understand such impacts, we built 3D fault‐slip finite element models to study the evolution of the pore pressure changes and the ΔCFS based on the fully coupled poroelastic theory. The strike‐slip, reverse, and normal faulting models were investigated separately. The patterns of co‐seismic ΔCFS were similar to the elastic stress changes on a homogeneous half‐space, but there were considerable differences in details especially in the near‐fault area when considering co‐seismic pore pressure change. Due to post‐seismic pore fluid flow, the near‐field stress shadow narrowed gradually in the strike‐slip model. In the reverse faulting model, both the near‐field stress shadow and enhanced areas expanded. In the normal faulting model, the stress shadow near the epicenter reduced. The pore pressure changes decayed sharply at the beginning and then gradually approached to zero because the system evolves to drained conditions. These results can help us to understand the effect of pore fluid flow on Coulomb Failure Stress evolution and might improve the understanding of earthquake triggering and aftershock forecasting.

Slow Slip Triggers the 2018 Mw 6.9 Zakynthos Earthquake Within the Weakly Locked Hellenic Subduction System, Greece

AbstractSlow slip events (SSEs) at subduction zones can precede large‐magnitude earthquakes and may serve as precursor indicators, but the triggering of earthquakes by slow slip remains insufficiently understood. Here, we combine geodetic, Coulomb wedge and Coulomb failure‐stress models with seismological data to explore the potential causal relationship between two SSEs and the 2018 Mw 6.9 Zakynthos Earthquake within the Hellenic Subduction System. We show that both SSEs released up to 10 mm of aseismic slip on the plate‐interface and were accompanied by an increase in upper‐plate seismicity rate. While the first SSE in late 2014 generated only mild Coulomb failure stress changes (≤3 kPa), that were nevertheless sufficient to destabilize faults of various kinematics in the overriding plate, the second SSE in 2018 caused stress changes up to 25 kPa prior to the mainshock. Collectively, these stress changes affected a highly overpressured and mechanically weak forearc, whose state of stress fluctuated between horizontal deviatoric compression and tension during the years preceding the Zakynthos Earthquake. We conclude that this configuration facilitated episodes of aseismic and seismic deformation that ultimately triggered the Zakynthos Earthquake.

Determining Stress State of Source Media with Identified Difference between Groundwater Level during Loading and Unloading Induced by Earth Tides

The groundwater level might be adopted as a useful tool to explore pre-seismic stress change in the earth crust, because it circulates in the deep crust and should be altered by the processes associated with the preparation of earthquakes. This work makes a new attempt that applies the load/unload response ratio (LURR) technique to study the stress state of the source media before the large earthquakes by calculating the ratio between the water levels during the loading and unloading phases. The change of Coulomb failure stress induced by earth tides in the tectonically preferred slip direction on the fault surface of the mainshock is adopted for differentiating the loading and unloading periods. Using this approach, we tested the groundwater level in the wells near the epicenters of some large earthquakes that occurred in the Sichuan-Yunnan region of southwest China. Results show that the LURR time series fluctuated narrowly around 1.0 for many years and reached anomalously high peaks 2~8 months prior to the mainshocks. For the earthquakes with multiple observation wells, the magnitude of the maximum values decreases with the distance from the epicenter. The underlying physics of these changes should be caused by the pre-seismic dilatancy. The corresponding volume variations in the crust could be observed in the geodetic time series in the same neighborhoods and during the same period.

Focal Mechanism Analysis of the September 25th, 2019 Mw 6.5 Ambon Earthquake and Its Implication for Seismotectonics

On September 25th, 2019, an Mw 6.5 earthquake occurred in Ambon, Maluku Province, Indonesia, and caused casualties and infrastructures damages. The epicenter located in a tectonically active region with the potential strike-slip and thrust faulting earthquake sources, yet the responsible fault is still not well understood. Based on focal mechanism solutions from available seismological agencies, i.e. USGS, GFZ, GCMT, and BMKG, the earthquake has a similar strike-slip focal mechanism, although there are discrepancies on detailed source parameters. To provide a better understanding of the earthquake mechanism and seismotectonic, we apply the Cut-and-Paste (CAP) focal mechanism inversion method to broadband seismic waveforms from regional and teleseismic distances. The CAP inversion results on the regional data grouped in different distance ranges show a robust strike-slip solution. We then refine the earthquake focal depth by performing the CAPtele inversion and resulted in a depth of 12 km with similar fault plane solution as the regionals. The ruptured fault plane is resolved by a directivity analysis using azimuthal pattern of the apparent source durations, which indicates an obvious unilateral rupture propagation toward SSE direction. Our result suggests the NNW-SSE orientated fault is the ruptured fault plane, which is also consistent with the near N-S distributed aftershocks. This fault is located in a narrow sea between Seram, Ambon and Haruku island and was not reported yet in previous studies. The Coulomb failure stress (CFS) changes analysis of the mainshock shows that the Ambon earthquake has promoted the off-fault aftershocks which occurred to the west of the ruptured fault.

Stress Transfer Along the Western Boundary of the Bayan Har Block on the Tibet Plateau From the 2008 to 2020 Yutian Earthquake Sequence in China

AbstractEight Ml ≥ 7.0 earthquakes have occurred around the Bayan Har block, NW Tibet, China, since 2000, resulting in a large number of casualties. Near the western boundary of the Bayan Har block, four Ml ≥ 6.0 Yutian earthquakes have occurred from 2008 to 2020. Stress interactions among them are comprehensively investigated by applying the ETAS (Epidemic‐Type Aftershock Sequence) model and calculating(Coulomb failure stress change) from a 3D linear viscoelastic mode. The combined (coseismic plus postseismic)induced by proceeding Yutian earthquakes on hypocenters of the 2012, 2014, and 2020 Yutian earthquakes are −1.5 × 10−4, 3.6 × 10−3, and 1.5 × 10−1 MPa, respectively. The background probabilities of the 2008, 2012, 2014, and 2020 Yutian earthquakes are 0.87, 0.97, 1.5 × 10−3, and 8.7 × 10−5, respectively. Combining those two independent approaches, we conclude that the 2008 and 2012 Yutian earthquakes are more like background earthquakes and that the 2014 and 2020 Yutian earthquakes were triggered by the proceeding Yutian earthquakes.

An integrated critical approach to off-fault strike-slip motion triggered by the 2011 Van mainshock (Mw 7.1), Eastern Anatolia (Turkey): New stress field constraints on subcrustal deformation

In this study, we retrieved the finite source characteristics of the October 23, 2011 Van earthquake (Mw 7.1) using the teleseismic waveforms to focus on the source location. The outstanding off-fault aftershock sequence of the Van mainshock was readily explained by calculating the Coulomb stress changes imparted to the surrounding crust. This may be accomplished through finite source modelling to examine the stress interaction between the fault, ruptured by the Van mainshock, and the surrounding fault(s) triggered by the same mainshock. In addition, to provide further support for the Coulomb failure stress changes in the off-fault area, centroid moment tensor (CMT) inversion of the off-fault aftershocks was performed and stress tensors were derived from their focal solutions. This identified the dominant fault slip, the constraints of the crustal stress fields and illuminated the crustal nature of the stress interaction. The off-fault aftershocks showed a strike-slip stress regime in rotational (to NW) and non-rotational (to N) stress fields of the upper and lower crusts, respectively. However, this was inconsistent with a horizontal compressional stress direction striking to the north. This suggests that a local source and/or rotation of lateral variation in stress magnitudes in crustal and sub-crustal structures strongly perturbed the regional stress field. It was also evident that these strike-slip aftershocks increased the intensity of stress in an off-fault area, NE of the source rupture. This revealed a uniquely triggered strike-slip motion, activated and rooted in the weak lower crust. We conclude that the Van mainshock rupture source area, associated with the stress changes imparted to the surrounding crust, had undergone anomalous modifications to generate distinctive off-fault aftershock responses in the entire crust, and also triggered and loaded the weak lower crust. We hypothesize that the strike-slip motion, the so called “transfer fault”, as a distinctly triggered slip event, was generated or selectively activated by subcrustal ductile processes in the absence of mantle lid beneath the study area. However, local slab fragmentation, tearing and cold mantle beneath the study area lead to paradigm changes in interpreting the strike-slip motion and subcrustal deformation. The presence of a small piece of oceanic lithosphere, consistent with fragmented, torn slab and cold mantle, may be an alternative hypothesis that remains to be tested. The Van earthquake, combined with careful examination of associated off-fault aftershocks, revealed new information about stress field constraints on subcrustal deformation. This investigation also provided insights into an important role of stress interaction, with a newly discovered transfer fault within the off-fault area, which extends through the entire crust beneath Lakes Van and Erçek areas.

The inferences on the Aegean geodynamic context from 30 October 2020 Samos earthquake (Mw:6.8)

On October 30, 2020, an earthquake with a magnitude of 6.8 (Mw) struck the northern coast of Samos Island in the Kuşadası Gulf. The solution to the focal mechanism indicates that the earthquake of 30 October 2020 occurred on a normal fault with nodal planes of E-W strike; thus, indicating extension in N-S direction. The fault plane solutions show a N-S trending extension for normal faults, which are obtained by inverting the moment tensor waveforms of 23 earthquakes and the P-wave first motion polarities of 11 aftershocks. A normal fault stress regime of approximately N-S (N6°E) σ3 axis is given by the inversion of slip vectors measured at sites located on land in Kuşadası. The mean R value is 0.84, suggesting that the stress regime is triaxial extensional stress state. The inversion of the focal mechanism of earthquakes occurring on land and in the Gulf of Kuşadası describes an extensional stress regime active today, characterized by an approximately N-S (N9°E) σ3 axis. The calculated R value of 0.31 indicates a triaxial stress state. For the 30 October 2020 earthquake (Mw:6.8), the Coulomb failure stress change analysis shows a substantial reduction in stress in the N-S direction supporting the kinematic results. The N-S extension in Western Anatolia-Aegean is largely influenced by the relatively fast movement of the Hellenic trench southwards, related to the sinking of the African plate beneath Aegean.

Geodetic Constraints on Recent Subduction Earthquakes and Future Seismic Hazards in the Southwestern Coast of Mexico

AbstractThree major subduction earthquakes occurred on March 20, 2012 (Mw 7.4), February 16, 2018 (Mw 7.2), and June 23, 2020 (Mw 7.4) in the southwestern coast of Mexico, which caused fatalities, casualties, considerable damage, and raised safety concerns about future seismic hazards. We use satellite geodetic observations to invert for the slip distributions of the three events and then investigate their interactions. Coulomb Failure Stress (CFS) induced by their slip both on surrounding thrust and normal faults are calculated. The positive CFS changes, along with the spatial‐temporal seismicity evolution, approximate earthquake recurrence rate and interseismic coupling, collectively indicate an increased possibility of a near‐future rupture around the areas between the 2018 and 2020 events in Oaxaca. Furthermore, there is a lowered chance of shallow coastal or offshore normal earthquakes but an increased chance of inland normal ruptures.

Coulomb stress changes associated with the M7.3 Maduo earthquake and implications for seismic hazards

On May 21, 2021, the M7.3 Maduo earthquake occurred within the Bayan Har Block, one of the most seismically active regions in the Tibetan Plateau. This aroused attention to the seismic hazards in this block. Understanding the stress transfer to adjacent faults is important for assessing seismic hazards. In this study, based on the elastic dislocation theory and a multi-layered lithospheric model, we calculate the Coulomb failure stress changes caused by this event on the adjacent major faults and identify the segments that experienced significant increases in stress. Our preliminary results show that five segments, including the middle segment of the Dari Fault, the Tuosuo Lake segment, the Maqin–Maqu segment, the western segment of the Dari Fault, and the Yushu segment of the Garze–Yushu Fault, experienced significant increases in co-seismic stress, which exceed the proposed threshold (1 ​× ​104 ​Pa) for triggering an earthquake. The stress change pattern associated with the post-seismic viscoelastic relaxation effect is substantially consistent with that caused by the co-seismic effect. The segments with high combined (co-plus post-seismic) stress changes are identical to those indicated by the co-seismic stress increase. Based on the stress change in this study and those of previous studies on seismic activity and fault stress evolution, we emphasize the seismic potential of the Maqin–Maqu segment, the Tuosuo Lake segment, and the middle-west segments of the Dari Fault, which should receive increased attention.

Decomposition of elastic potential energy and a rational metric for aftershock generation

SUMMARY The occurrence of earthquakes can be regarded as shear fracture releasing the elastic potential energy stored in the Earth. The potential energy density of linear elastic forces is generally represented in the quadratic form of strain tensor components with the fourth-order coefficient tensor of elastic stiffness. When the material is isotropic, since the stiffness tensor is expressible as a linear combination of two independent symmetric tensors, we can decompose the elastic potential energy density into two independent parts, namely the volumetric part and the shearing part. By definition, the partial derivatives of the elastic potential energy density with respect to volumetric and shearing deformations give the corresponding generalized forces in the sense of Lagrangian mechanics: specifically, one-third of the first invariant of stress tensor (equivalent to the mean stress) to volumetric deformation and the square root of the second invariant of deviatoric stress tensor (equivalent to $\sqrt {3/2} $ times the octahedral shear stress). With these generalized forces instead of the normal and tangential stresses on a specific fault plane, we correctly represented the original concept of Coulomb's failure criterion (shear failure occurs when shearing stress is equal to shearing strength) and defined energetics-based failure stress (EFS). The change in EFS associated with the occurrence of a main fracture (ΔEFS) gives a rational metric for aftershock generation, which can be reduced to previously proposed various metrics in special cases. For example, when the level of background deviatoric stress is much higher than the magnitude of coseismic stress changes, the expression of ΔEFS is reduced to a similar form to the well-known Coulomb failure stress change (ΔCFS). Even in the energetics-based metric, the effects of pore-fluid pressure changes are essential. We theoretically examined the mechanical effects of induced and enforced pore-fluid pressure changes and elucidated that the difference between them is reflected in the focal mechanisms of aftershocks.

Assessment of the Seismic Hazards of the Marikina Valley Fault from 2019 Mw 6.1 Castillejos Earthquake and Historical Events

Abstract The 2019 Mw 6.1 Castillejos earthquake occurred in the Zambales range of the central Luzon Island in Philippines. No active fault was reported around the seismogenic zone according to previous investigations. This earthquake draws attention for assessment in seismic risk along the Marikina Valley fault system (MVFS) near the Manila dense metropolitan population. The Coulomb failure stress (CFS) change on the MVFS is estimated by the coseismic faulting model derived from the inversion of coseismic deformation field observed from the Differential Interferometric Synthetic Aperture Radar using both the Advanced Land Observing Satellite-2 and Sentinel-1 Synthetic Aperture Radar (SAR) images. The predicted CFS change is less than 0.5 kPa that implies insignificant Coulomb stress accumulation on the MVFS after the Mw 6.1 Castillejos event. However, the recorded 14 moderate and strong earthquakes in and around the Luzon islands caused significant CFS drop on the MVFS. This might delay the occurrence of the earthquake for 0.2–50 yr on the MVFS.

Investigation on Two Mw 3.6 and Mw 4.1 Earthquakes Triggered by Poroelastic Effects of Hydraulic Fracturing Operations Near Crooked Lake, Alberta

AbstractA coupled approach of fluid flow and geomechanics is proposed in this study to quantitatively understand the hydraulic fracturing‐induced poroelastic effects that activate pre‐existing faults and trigger earthquake swarms near Crooked Lake, Alberta. The 3D poroelastic simulation of the seismogenic fault zone is then conducted to characterize the pressure diffusion and stress perturbation that led to fault activation. Results show that the poroelastic effects on the high‐permeable damage zones of a conductive‐barrier fault triggered the sequential activation of the seismogenic fault in the basement and Winterburn Formation. In addition, the high‐permeable damage zones act as conduits for pore pressure diffusion along the fault, whereas the fault core functions as a barrier to prevent crossing flow. The poroelastic effects in terms of pressure diffusion and stress perturbation in response to fluid injection facilitated the fault slip and hence triggered the Mw 3.6 earthquake in the basement formation 40 days after the fracturing operations. Moreover, the mainshock created negative Coulomb failure stress changes, inhibiting further fault slip in the basement formation. Subsequently, the stimulation of another well facilitated poroelastic effects on the fault damage zone in top Winterburn Formation, reactivating the same fault with an Mw 4.1 earthquake 8 days after the initiation of treatments. It is essential to optimize the injection site selection near the existing faults to reduce risks of the induced earthquakes near Crooked Lake.

Normal Faulting Movement During the 2020 Mw 6.4 Yutian Earthquake: A Shallow Rupture in NW Tibet Revealed by Geodetic Measurements

The ENE striking Longmu Co fault and the North Altyn Tagh left-lateral slip fault have led to the complex regional structure in the northwestern Tibetan Plateau, resulting in a series of normal faulting and strike slip faulting earthquakes. Using both the ascending and descending Sentinel-1A/B radar images, we depict the coseismic deformation caused by the 2020 Yutian Mw 6.4 earthquake with a peak subsidence of ~ 20 cm. We determine the seismogenic fault geometry by applying the Bayesian approach with a Markov Chain Monte Carlo sampling method, which can better characterize the posterior probability density functions of the source model parameters. The estimation results reveal that the earthquake is a normal faulting event with a moderate strike slip component. Based on the optimal fault geometry model, we extend the fault plane and invert for the distributed coseismic slip model. The optimal slip model shows that the coseismic slip is mainly concentrated at shallow depths of 3–10 km with a maximum slip of ~ 1.0 m. Our preferred geodetic coseismic model exhibits no surface rupture, which may likely be due to the shallow slip deficit in the uppermost crust. We calculate the combined loading effect of the Coulomb failure stress changes induced by the coseismic dislocations and postseismic viscoelastic relaxation of the 2008 Mw 7.1, 2012 Mw 6.4 and 2014 Mw 6.9 Yutian events. Our study demonstrates that the three preceding major Yutian shocks were insufficient to trigger the 2020 Yutian earthquake, which we consider perhaps reflects the natural release of elastic strain accumulated mainly through localized tectonic movement. We attribute the 2020 Yutian event to the release of extensional stress in a stepover zone controlled by the Longmu Co and the North Altyn Tagh sinistral strike slip fault systems. The seismic risk in the southwest end of the North Altyn Tagh fault has been elevated by the Yutian earthquake sequences, which require future attention.

Post-seismic crustal internal deformation in a layered earth model

SUMMARY This study introduces a novel method for computing post-seismic crustal internal deformation in a layered earth model. The surface dislocation Love number (DLN) calculated by the reciprocity theorem was implemented as the initial value. Furthermore, numerical integration of the value from the Earth's surface to the interior was undertaken to obtain the internal DLN. This method does not require a combination of the general solution and particular solution for the calculation of internal deformation above the seismic source, thus avoiding the loss of precision. When the post-seismic deformation within a certain period is calculated, the particular solutions at the beginning and end of the considered period cancel each other. This simplifies the calculation of post-seismic internal deformation. The numerical results depict that as the degrees increase, the post-seismic DLN reaches stability in a shorter interval of time. Thus, for improved efficiency of the post-seismic internal deformation calculation, the post-seismic DLNs should be calculated within 2000 degree and integrated with the coseismic results. As an application, the post-seismic Coulomb failure stress changes (∆CFS) induced by the 2011 Tohoku-Oki earthquake in the near field around the Japanese archipelagos and two major faults in Northeast China were simulated. The results exhibit that the ∆CFS values in the near field agree well with those simulated by the method in a half-space layered earth model, thus verifying the present method. The coseismic ∆CFS on the Mishan-Dunhua fault in Northeast China, as an example, is only 0.094–0.668 KPa. However, the ∆CFS caused by the viscoelastic relaxation of the mantle within 5 yr following the 2011 Tohoku-Oki event on the same fault exceeds the coseismic results. Therefore, the cumulative effect of the viscoelastic relaxation of the mantle is deserving of attention.

Sentinel-1 InSAR observations of co- and post-seismic deformation mechanisms of the 2016 Mw 5.9 Menyuan Earthquake, Northwestern China

On January 21, 2016, a Mw 5.9 earthquake occurred in Menyuan, Qinghai Province, in northwestern China, which is an earthquake-prone region located along the northeastern margin of the Tibetan Plateau. However, the source parameters, deformation mechanism, especially the influence of this earthquake on the surrounding regions and on post-seismic deformation, remains poorly understood. In this study, Sentinel-1 images from ascending and descending tracks were used to investigate the co- and post-seismic deformation fields. Uniform and distributed slip models were used to determine the source parameters. Finally, the co- and post-earthquake two-dimensional (2D) deformation decomposition, coulomb failure stress changes caused by the 2016 Menyuan earthquake, and the influence of this event on the surrounding region and active faults were analyzed. The results suggest that the 2016 Menyuan earthquake was triggered by a concomitant buried reverse fault in the northeastern part of the Northern Lenglongling Fault, which is a SW dipping sinistral thrust fault with a strike angle of 140°. The fault properties are consistent with the results of Interferometric Synthetic Aperture Radar (InSAR) 2D decompositions. The co-seismic deformation of the two satellite tracks have maximum line-of-sight uplifts of 55 mm and 65 mm, without notable subsidence. The post-seismic effect of the earthquake mainly affected the area between the Lenglongling Fault and the Northern Lenglongling Fault. The earthquake caused stress loading on the southwest of the western segment of the Lenglongling Fault and on the southeast of the eastern segment of the Northern Lenglongling Fault, as well as along these two faults.