Search

Studies have demonstrated the potential of satellite thermal infrared observations to detect anomalous signals preceding large earthquakes. However, the lack of well-defined precursory characteristics and inherent complexity and stochasticity of the seismicity continue to impede robust earthquake forecasts. This study investigates the potential of pre-seismic thermal anomalies, derived from five satellite-based geophysical parameters, i.e., skin temperature, air temperature, total integrated column water vapor burden, outgoing longwave radiation (OLR), and clear-sky OLR, as valuable indicators for global earthquake forecasts. We employed a spatially self-adaptive multiparametric anomaly identification scheme to refine these anomalies, and then estimated the posterior probability of an earthquake occurrence given observed anomalies within a Bayesian framework. Our findings reveal a promising link between thermal signatures and global seismicity, with elevated forecast probabilities exceeding 0.1 and significant probability gains in some strong earthquake-prone regions. A time series analysis indicates probability stabilization after approximately six years. While no single parameter consistently dominates, each contributes precursory information, suggesting a promising avenue for a multi-parametric approach. Furthermore, novel anomaly indices incorporating probabilistic information significantly reduce false alarms and improve anomaly recognition. Despite remaining challenges in developing dynamic short-term probabilities, rigorously testing detection algorithms, and improving ensemble forecast strategies, this study provides compelling evidence for the potential of thermal anomalies to play a key role in global earthquake forecasts. The ability to reliably estimate earthquake forecast probabilities, given the ever-present threat of destructive earthquakes, holds considerable societal and ecological importance for mitigating earthquake risk and improving preparedness strategies. [ABSTRACT FROM AUTHOR]

[Display omitted] [ABSTRACT FROM AUTHOR]

Southeastern (SE) Tibet is one of the most seismically active regions in mainland China, but the spatial distribution of interseismic coupling that quantifies seismic hazard is unknown along most major faults. In this study, we constructed an elastic block model to invert Global Positioning System data for slip rates and locking coefficients along 20 major faults in SE Tibet. Our results identify 27 strongly coupled fault segments with locking coefficients >0.5, defined as potential seismogenic asperities, extending laterally for 36–330 km. Quantitative calculations of seismic moment budgets on these seismogenic asperities indicate that they are capable of generating Mw 6.4–7.7 earthquakes in the next few decades, of which the Anninghe, Daliangshan and Red River faults have the potential for Mw ≥ 7.5 earthquakes. The interseismic coupling model provides a component for probabilistic analysis of future seismic hazards in densely populated Southwest China. Plain Language Summary: In the period between earthquakes, crustal faults are either partially to fully stuck or slipping freely. In the former case, plate movement can cause a buildup of energy on the fault plane and eventually be released through earthquakes. Using Global Positioning System (GPS) to measure the movement of the Earth's surface, we can use models of crustal deformation to identify where and how quickly faults are accumulating energy and then make reasonable predictions about the locations and magnitudes of future earthquakes. Based on the above principles, in this study, we use GPS measurements to suggest the locations and magnitudes of potential earthquakes on 20 major faults in southeast Tibet, which includes most of the Sichuan and Yunnan provinces in southwestern China. We find that 27 fault segments with lengths of 36–330 km are stuck to some degree and accumulating energy. Currently, the accumulated energy can produce earthquakes of Mw 6.4–7.7, highlighting the severity of earthquake hazards in Southwest China. Key Points: We obtained a unified interseismic coupling model for 20 major fault systems in southeastern TibetWe identify 27 seismogenic asperities with locking coefficients >0.5 that could host Mw 6.4–7.7 earthquakesThe interseismic coupling model provides a component for probabilistic analysis of future seismic hazards in Southwest China [ABSTRACT FROM AUTHOR]

Strain accumulation in foreland zones is important for understanding mountain‐building processes and seismic hazards. Material heterogeneity in this border zone, particularly the contrast between sedimentary cover and basement, affects strain accommodation and fault behavior. One manifestation of this effect is depth separation of Interferometric Synthetic Aperture Radar (InSAR)‐derived slip models and corresponding aftershock clusters; hypotheses for this vertical separation remain controversial. In this study, we investigated strain accumulation in the Kalpin fold‐and‐thrust belt (KFTB) of southwestern Tian Shan, China, using the integration of InSAR measurements, teleseismic body‐waves, near‐field strong motion data (SM), and relocated aftershocks of the 2020 Mw 6.0 Kalpin earthquake. The SM modeling and analysis of post‐seismic InSAR time series are performed for the first time for this earthquake. Our results confirm a vertical separation of the mainshock and aftershock cluster, with the former on a weak décollement and the latter on faults within the basement. InSAR time series analysis shows that post‐seismic deformation was dominated by afterslip on a splay fault directly above the ruptured décollement. Static stress transfer of co‐seismic rupture cannot explain these observations. We speculate that fluid flow and high pore pressure along pre‐existing fault planes may have reduced the fault strength and been involved in the evolution of this aftershock cluster. We conclude that compressive strain in the KFTB is accommodated by a mixture of thin‐ and thick‐skin faulting and seismic deformation across the entire crustal thickness. Specifically, we suggest that the observed shortening is partly accommodated by infrequent large earthquakes on the weak décollement. Plain Language Summary: In foreland areas, strain can be accommodated by earthquakes, aseismic slip, or both. The contrasting material properties of the sediment and basement certainly affects this process, for example, the depth separation of different fault behavior. However, it is unclear how the fault behavior varies in depth from basement to weak sediment throughout different stages of the seismic cycle. 2020 Mw 6.0 Kalpin earthquake occurred in Kalpin fold‐and‐thrust belt, southwestern Tian Shan provides some insight on this matter. Here, we combined Interferometric Synthetic Aperture Radar, aftershock relocation, teleseismic body‐waves, and strong motion data to investigate the strain release behavior during and immediately after the earthquake. We found that strain release behavior differs between sediment and basement during and immediately after the earthquake. The mainshock ruptured the weak décollement and the aftershocks were primarily rooted in the basement, while the transient afterslip was mainly distributed in the sediments. Two findings are worth noting: (a) the afterslip occurred on a fault, whose stress was partially unloaded during the earthquake, suggesting the fault strength was reduced, and (b) the earthquake rupture of the décollement highlights how destructive earthquakes can initiate on the weak décollement. Key Points: The moderate 2020 Mw 6.0 Kalpin earthquake was initiated on a weak décollement in southwest Tian ShanSeismological data and Interferometric Synthetic Aperture Radar postseismic time series analysis confirm a vertical separation of the mainshock and aftershocksCombination of thin‐skin and thick‐skin deformation processes throughout the different stages of the seismic cycle [ABSTRACT FROM AUTHOR]

Pre-seismic anomalies have the potential to indicate imminent strong earthquakes in the short to medium terms. However, an improved understanding of the statistical significance between anomalies and earthquakes is required to develop operational forecasting systems. We developed a temporal integrated anomaly (TIA) method to obtain the temporal trends of multiparametric anomalies derived from the Atmospheric Infrared Sounder (AIRS) product before earthquakes. A total of 169 global earthquakes that occurred from 2006 to 2020 and had magnitudes of ≥7.0 and focal depths of ≤70 km were used to test this new method in a retrospective manner. In addition, 169 synthetic earthquakes were randomly generated to demonstrate the suppression capacity of the TIA method for false alarms. We identified four different TIA trends according to the temporal characteristics of positive and negative TIAs. Long-term correlation analyses show that the recognition ability was 12.4–28.4% higher for true earthquakes than for synthetic earthquakes (i.e., higher than that of a random guess). Incorporating 2–5 kinds of TIAs offered the best chance of recognizing imminent shocks, highlighting the importance of multiparameter anomalies. Although the TIA trend characteristics before the earthquakes were not unique, we identified certain unexplained pre-seismic phenomena within the remote sensing data. The results provide new insight into the relationships between pre-seismic anomalies and earthquakes; moreover, the recognition ability of the proposed approach exceeds that of random guessing. [ABSTRACT FROM AUTHOR]

The challenge of ruling out potential rupture nodal planes with opposite dip orientations during interferometric synthetic aperture radar (InSAR)-based kinematic inversions has been widely reported. Typically, slip on two or more different fault planes can match the surface deformation measurements equally well. The ambiguous choice of the nodal plane for the InSAR-based models was thought to be caused by InSAR's 1-D measurement and polar orbiting direction, leading to its poor sensitivity to north–south crustal motion. Through synthetic experiments and simulations, this paper quantitatively demonstrates the main reason of the ambiguous InSAR-based models, which confuse researchers in the small-to-moderate thrust earthquake cases investigation. We propose the inherent 1-D measurement is not the principle cause of the fault plane ambiguity, since models derived from the same InSAR data predict similar, but not identical, 3-D deformation patterns. They key to differentiating between these different models is to be able to resolve the small asymmetry in the surface deformation pattern, which may be smaller in amplitude than the typical noise levels in InSAR measurements. We investigate the fault geometry resolvability when using InSAR data with different noise levels through 'R' value. We find that the resolvability does not only rely on the InSAR noise, but also on the fault geometry itself (i.e. depth, dips angle and strike). Our result shows that it is impossible to uniquely determine the dip orientation of thrust earthquakes with M w < 6.0 and depth > 5.0 km with InSAR data at a noise level that is typical for mountain belts. This inference is independent from the specific data set (i.e. interferogram or time-series) and allows one to assess if one can expect to be able to resolve the correct fault plane at all. [ABSTRACT FROM AUTHOR]

Characterizing the coseismic slip behaviors of earthquakes could offer a better understanding of regional crustal deformation and future seismic potential assessments. On 18 December 2023, an Mw 6.0 earthquake occurred on the Lajishan–Jishishan fault system (LJFS) in the northeastern Tibetan Plateau, causing serious damage and casualties. The seismogenic fault hosting this earthquake is not well constrained, as no surface rupture was identified in the field. To address this issue, in this study, we use Interferometric Synthetic Aperture Radar (InSAR) data to investigate the coseismic surface deformation of this earthquake and invert both ascending and descending line-of-sight observations to probe the seismogenic fault and its slip characteristics. The InSAR observations show up to ~6 cm surface uplift caused by the Jishishan earthquake, which is consistent with the thrust-dominated focal mechanism. A Bayesian-based dislocation modeling indicates that two fault models, with eastern and western dip orientations, could reasonably fit the InSAR observations. By calculating the coseismic Coulomb failure stress changes (∆CFS) induced by both fault models, we find that the east-dipping fault scenario could reasonably explain the aftershock distributions under the framework of stress triggering, while the west-dipping fault scenario produced a negative ∆CFS in the region of dense aftershocks. Integrating regional geological structures, we suggest that the seismogenic fault of the Jishishan earthquake, which strikes NNE with a dip of 56° to the east, may be either the Jishishan western margin fault or a secondary buried branch. The optimal finite-fault slip modeling shows that the coseismic slip was dominated by reverse slip and confined to a depth range between ~5 and 15 km. The released seismic moment is 1.61 × 1018 N·m, which is equivalent to an Mw 6.07 earthquake. While the Jishishan earthquake ruptured a fault segment of approximately 20 km, it only released a small part of the seismic moment that was accumulated along the 220 km long Lajishan–Jishishan fault system. The remaining segments of the Lajishan–Jishishan fault system still have the capability to generate moderate-to-large earthquakes in the future. [ABSTRACT FROM AUTHOR]

Large block‐bounding faults on the Tibetan plateau are significant geological structures that accommodate tectonic movements and accumulate stress, leading to large earthquakes. Quantifying the interseismic slip deficit rate helps to better assess the earthquake potential. We combine available InSAR (2015–2020) and interseismic GPS data to determine fault coupling along 14 major block‐bounding faults. Spatially dense InSAR measurements remarkably improve the resolution of the coupling model. Combined with a GPS‐constrained block model, we examine the performance of the inversion approach with the stress constraint and the common Laplacian smoothing based on both synthetic tests and real data. We suggest that, for continental strike‐slip faults, adding the stress constraint can mitigate unphysical coupling distributions due to unreasonable assumptions or modeling artifacts, reducing the model uncertainty and improving the accuracy of the coupling model. This is particularly useful for segments featured by a highly heterogeneous distribution of coupling along the transition zone from locking to creeping region, partially‐coupling segment, and junction zone between main and subsidiary faults. We present a large‐scale fault coupling map along the major block‐bounding faults on the northeastern and eastern Tibetan plateau, highlighting the distinct degrees of fault coupling and lateral variations. The collage of coupling maps along different faults demonstrates the kinematic features over a broad time scale during earthquake cycles ranging from early to late interseismic phases, such as the segments ruptured during the 2001 Kokoxili earthquake and the 1920 Haiyuan earthquake. Plain Language Summary: Large faults at the boundary of blocks accommodate a significant portion of the relative movement between adjacent blocks and accumulate tectonic stress during the interseismic period. Consequently, these faults have a high potential to generate large earthquakes. To evaluate the distribution of strain accumulation, we analyze InSAR (2015–2020) and interseismic GPS measurements to construct a coupling model, which tells us the extent to which two blocks are locked together and accumulate stress before the next earthquake. We focus on 14 major faults on the eastern and northeastern Tibetan plateau. The addition of the spatially dense InSAR data helps to resolve the coupling distribution in greater detail, where GPS observations are limited. We evaluate different approaches to determine fault coupling, including stress constraints and smoothing techniques. Our findings indicate that adding a stress constraint performs better than the smoothing operator and helps mitigate physically unrealistic coupling distributions in the inversion using the spatial smoothing method. The stress constraint improves the accuracy of the coupling model. Based on the stress‐constrained inversion, we present a comprehensive fault coupling map in the northeastern and eastern Tibetan plateau, highlighting the spatial variations of fault coupling, which can be used in earthquake hazard assessment. Key Points: Stress constraints can mitigate unphysical coupling distributions and reduce the uncertainty of the coupling model for continental faultsWe highlight the improvement when using stress constraints in the case of sharp coupling contrasts along the strike‐slip faultWe provide a new coupling model along 14 major block‐bounding faults on the eastern and northeastern Tibetan plateau [ABSTRACT FROM AUTHOR]

Pre-seismic thermal anomalies have been widely reported, however, there is still no unified understanding of their causes or physical mechanisms. To clarify whether subsurface fluid convection and rock deformation directly related to seismic fault stress release can form thermal anomalies before an earthquake, we carried out thermo-hydro-mechanical (THM) coupled modeling of the Wenchuan earthquake, and considered the influence of maximum horizontal principal stress direction (MPSD) and fault dip. The results show that convective thermal flux (CT) anomalies of subsurface fluid and the deformation induced thermal flux (DT) of the rock mass caused by stress release of the seismic fault can explain the observed thermal anomalies along the fault. The calculated and observed thermal anomalies are comparable in spatial distribution, evolution, and magnitude. Specifically, for the Wenchuan earthquake, seismic fault stress release-induced CT and DT anomalies occur in the fault zone and in the adjacent hanging wall area. In the fault zone, CT and DT induced thermal anomalies were comparable, and superposition of the two was characterized by strong initial warming, followed by weaker warming and then short-term cooling. However, outside of the fault zone, CT induced thermal anomalies had much higher intensities than those induced by DT, and were characterized by band-like warming along the hanging wall of the fault. In addition, variations in the MPSD and fault dip cause the segmentation of thermal anomalies along the Wenchuan earthquake fault. From the southern to the northern segments of the fault, the intensity of thermal anomalies is significantly weakened with rotation of MPSD and an increase in fault dip. Our study provides evidence that seismic fault stress release is the cause of pre-seismic thermal anomalies; furthermore, the results indicate the CT and DT could be important physical mechanisms of pre-seismic thermal anomalies. [ABSTRACT FROM AUTHOR]

Knowledge of the slip behavior of the Altyn Tagh fault (ATF) has significant implications for our understanding of the tectonic deformation of the Tibetan Plateau. In this study, we process Global Positioning System (GPS) data spanning 2009–2017 across the western ATF, merge the solution with recently published GPS velocities, and obtain a dense GPS velocity field for northern Tibet. We introduce an elastic block model and estimate the fault slip rate, interseismic fault coupling (ISC), and seismic moment accumulation rate along the ATF. The estimated left‐lateral strike‐slip rate of the ATF decreases eastward from 12.8 ± 0.4 mm/a to 0.1 ± 0.2 mm/a. Results show a heterogeneous distribution of ISC along the fault, with the fault locking depth varying from 5 to 20 km. The seismic moment accumulation rate is 2.16–2.37 × 1018 N m/km along three segments between Sulamu Tagh and Aksay bends, where the accumulated seismic moment over the last five centuries could be balanced by earthquakes with magnitudes of Mw 7.7, Mw 7.6, and Mw 7.8, respectively. We calculate the ratio of the Indo‐Eurasia convergence accommodated by the ATF and find that 10.3% of the convergence is accommodated along the fault by lateral extrusion; quantitative analysis of the strain rate, however, shows that a large part of the northern Tibet is not well described by elastic block rotation. Consequently, we suggest that a hybrid kinematic model that includes both block‐like and continuous deformation is needed to better delineate the crustal deformation of the Tibetan Plateau. Key Points: Altyn Tagh fault slip rate and coupling inverted using new, dense interseismic GPS velocitiesThe three fault segments between Sulamu Tagh and Aksay bends have the potential to rupture in Mw 7.7, Mw 7.6, and Mw 7.8 earthquakesA hybrid kinematic model is needed to delineate the crustal deformation of the Tibetan Plateau [ABSTRACT FROM AUTHOR]

Previous studies using paleoseismic trenches and numerical models proposed two different fault rupture models for the Pingdingshan-Annanba segments of the Altyn Tagh Fault, which are both poorly validated because of a lack of measurements of surface rupture displacements. In this paper, we use high-resolution satellite imagery to measure offset geomorphic markers. Taking the fault geometry into consideration, the most likely paleoseismic rupture parameters are determined by comparing the results of the cumulative offset probability densities (COPDs) and the empirical relationships of seismic parameter. The results show that the most recent earthquake (1270A.D.-1775A.D.) along the Suoerkuli segment only ruptured the Suoerkuli Valley, with a surface rupture length of 140 km, an average horizontal displacement of 3.12 m and a moment magnitude (Mw) of 7.4–7.6. The penultimate earthquake (676A.D.-1347A.D.) ruptured both the Annanba segment and part of the Suoerkuli segment with an average horizontal displacement of greater than 5 m, a surface rupture length of 120–260 km and Mw 7.4–7.8. The rupture parameters of the above two earthquakes demonstrate that the Suoerkuli and Annanba segments are individual fault segments which can rupture independently with a surface rupture length more than 100 km. The Pingdingshan double bend to the west of the Suoerkuli segment is likely to be an endpoint that inhibits the rupture propagation towards further west. [ABSTRACT FROM AUTHOR]

In this study, The ENVISAT advanced synthetic aperture radar observations from 2003 to 2010 of a descending track covering an area of 100 km × 300 km were used to map the surface velocity field in northwestern Tibet. The derived line-of-sight (LOS) velocity map revealed that interseismic deformation was mainly located on the Altyn Tagh Fault (ATF) and other four immature subsidiary faults (i.e., Tashikule Fault, Muzitage-jingyuhe Fault, Heishibeihu Fault, and Woniuhu Fault). A 2D elastic screw dislocation model was used to interpret the interferometric synthetic aperture radar (InSAR) velocity profiles, which revealed the following results. (a) The oblique movement is partitioned between left-lateral slip at a rate of 6.3 ± 1.4 mm/y on the ATF and 5.9 ± 2.8 mm/y on the subsidiary faults. The low slip rate of the ATF indicates that the ATF does not drive the northeastward extrusion of material, with most of the extrusion occurring in the eastern interior of the plateau and the four subsidiary faults localizing the oblique convergence partitioned in the west. This can reasonably explain why catastrophic earthquakes and rapid slip do not occur all over along the ATF. (b) Based on the four subsidiary faults accommodating the oblique movement and the traces amalgamation with the EKLF (delineated Bayan Har plate boundary to the northeast), we concluded guardedly that the four subsidiary faults are the evoluting plate boundary of the Bayan Har block to the northwest. (c) The Tanan top-up structure had an uplift rate of ~0.6 mm/y at the south of the Tarim Basin. [ABSTRACT FROM AUTHOR]

Along the margins of orogenic plateaus, the viscous Earth structure and fault geometries play a primary role in controlling the tectonic evolution and earthquake generation. After the 2008 M w 7.9 Wenchuan earthquake, the long-standing debate regarding the tectonics producing and maintaining prominent topography across the Longmen Shan reignited. Post-seismic deformation, representing the surface strain history in response to lithospheric stress perturbations, provides important insights into the lithospheric rheology and active structures. Here, we construct a new 3-D post-seismic deformation model for the Wenchuan earthquake, invoking viscoelastic relaxation and afterslip. Our best-fitting model indicates that the steady-state viscosities of the lower crust in the region to the immediately west of the Songpan-Ganzi terrane and beneath the Songpan-Ganzi terrane are estimated to be 4.0 × 1018 and 1.0 × 1018 Pa s, respectively. Our results, combining geophysical and geodetic observations and model analyses, highlight the prevalent parallelism between the rheological and structural boundaries of the lower crust, which diverge northward away from the trend of the Longmen Shan fault at ∼20°. This diverging rheological structure and the partially coupled upper and lower crust have broad implications for the stress build-up, strain partitioning and deformation styles along the eastern Tibetan Plateau margin. [ABSTRACT FROM AUTHOR]

Based on the dense GPS velocity field in the northeastern margin of the Tibetan Plateau from 1999 to 2016, we have produced the deformation and strain characteristics of the Haiyuan fault and the Liupanshan fault. Estimated long-term slip rate along the Haiyuan-Liupanshan fault zones show a gradual decrease from 6.4 ± 1.6 mm/yr at the Tuolaishan fault to 2.9 ± 1.2 mm/yr at the Southern Liupanshan fault. Left-lateral thrusting movement was inverted for the Xiangshan-Tianjingshan fault (XS-TJS), which has an average slip rate of 2.1 ± 3.4 mm/yr during the study period. We also calculated the heterogeneous distribution of interseismic coupling along the fault zones. Our result also shows the locking depth of the Tianzhu seismic gap is ∼22 km. The slip rate deficit, the seismic moment accumulation rate, and the Coulomb stress accumulation rate are high on the fault planes, whereas the second invariant of the strain rate is low at the surface. The Liupanshan fault is locked to a depth of ∼23 km, and the corresponding seismic moment accumulation rate on the fault plane is high, while the strain rate at the surface is low. The accumulated strain along the Tianzhu seismic gap and the Liupanshan fault could be balanced by earthquakes with magnitudes of Mw7.9 and Mw7.4, considering the absence of large earthquakes over the last 1000 years and 1400 years respectively. The Haiyuan segments had ruptured during 1920 Haiyuan earthquake, and the estimated locking depth for period 1999–2016 is 5–10 km. Its seismic moment accumulation rate at depth is low and the strain rate at the surface is high. Our result indicates that 70% of the strike-slip along the Haiyuan segments transforms into thrusting along the Liupanshan fault, while the remaining 30% is related to the orogeny of the Liupanshan. For slip between the Haiyuan fault and the XS-TJS, about 27–34% of the slip is partitioned on the XS-TJS. [ABSTRACT FROM AUTHOR]

Earthquake forecasting aims to determine the likelihood of a damaging earthquake occurring in a particular area within a period of days to months. This provides ample preparation time for potential seismic hazards, resulting in significant socioeconomic benefits. Surface and atmospheric parameters derived from satellite thermal infrared observations have been utilized to identify pre-earthquake anomalies that may serve as potential precursors for earthquake forecasting. However, the correlation between these anomalies and impending earthquakes remains a significant challenge due to high false alarm and missed detection rates. To address this issue, we propose a spatially self-adaptive multiparametric anomaly identification scheme based on global strong earthquakes to establish the optimal recognition criteria. Each optimal parameter exhibits significant spatial variability within the seismically active region and indicates transient and subtle anomaly signals with a limited frequency of occurrences (<10 for most regions). In comparison to the fixed criterion for identifying anomalies, this new scheme significantly improves the positive Matthew's correlation coefficient (MCC) values from ~0.03 to 0.122–0.152. Additionally, we have developed a multi-parameter anomaly synthesis method based on the best MCC value of each parameter anomaly. On average, the MCC increased from 0.143 to 0.186, and there are now more earthquake-prone regions with MCC values > 0.5. Our research emphasizes the critical importance of a multiparametric system in earthquake forecasting, where each geophysical parameter can be assigned a distinct weight, and the findings specifically identify OLR, including all-sky and clear-sky ones, as the most influential parameter on a global scale, highlighting the potential significance of OLR anomalies for seismic forecasting. Encouraging results imply the effectiveness of utilizing multiparametric anomalies and provide some confidence in advancing our knowledge of operational earthquake forecasting with a more quantitative approach. [ABSTRACT FROM AUTHOR]

Based on the working principle of satellite radars, the earthquake deformation field measured by interferometric synthetic aperture (InSAR) is the projection of ground displacement associated with the seismogenic fault in the line of sight (LOS) of the satellite. However, LOS projections are complex, and are not only related to the ascending/descending modes and incidence angles of SAR data, but also related to the strike and motion senses of the fault. Even for the same earthquake, the LOS deformation derived from different ascending/descending data can be almost identical in one case, but quite different in another case, which makes the interpretation of InSAR seismic deformation and its comparison with field observations difficult. In this study, we undertook a quantitative analysis of the relationships between LOS observation sensitivity of InSAR and fault strike, fault motion sense, and ascending/descending modes, as well as 3D deformation fields. We studied the features and differences of the LOS deformation fields in different types of earthquakes using ascending/descending modes, with a particularly detailed analysis of the relations for a strike-slip type of earthquake. We also summarized the characteristics of LOS deformation fields of faults with different strikes and optimal observational data modes. Taking the strike-slip Yushu earthquake and the normal Gaize event as examples, we used SAR data of the ascending/descending modes to verify the results of quantitative calculations. These analyses will not only provide a more reasonable interpretation of InSAR seismic deformation fields and but also help understand the differences of seismic deformation fields revealed by data with different observational modes, therefore promoting the application of InSAR technology in seismology. [ABSTRACT FROM AUTHOR]

The July 3rd, 2015 Pishan, China Mw 6.4 earthquake occurred within the Hetian fold belt, a frontal thrust region between the Northwestern Tibetan Plateau and the Tarim basin. We investigate the fault geometry and the rupture process of the Pishan earthquake based on joint inversion of teleseismic body waves and InSAR measurements. Our results show that the top of the fault that ruptured in the Pishan earthquake is buried 5 ± 2 km beneath the surface and that the earthquake ruptured only a segment of a deeply-seated detachment lying under the Hetian fold belt. Our inferred coseismic slip model shows two slip asperities, located in the depth range of about 7–9 km and 10–11 km. Additionally, we resolve a slip deficit region between the slip asperities, which might indicate friction heterogeneity on the fault plane. The earthquake ruptured for ∼15 s, releasing a total seismic moment of 4.7 × 10 18 N m. We also show that the latest Pishan earthquake caused positive Coulomb stress changes on back thrust faults, potentially decreasing the time to the next earthquake on those faults. Based on our Pishan coseismic slip model, along with both the trend of the anticline system and the pattern of coseismic deformation, we envision an ongoing east-west extensional growth of the anticlines in the Western Kunlun frontal thrust region in addition to the northern underthrusting of the Northwestern Tibetan Plateau below the Tarim basin. [ABSTRACT FROM AUTHOR]

GPS-derived horizontal velocities, geologic fault slip rates and earthquake-derived fault slip vector azimuths are inverted simultaneously for fault coupling and slip rate deficit on the Haiyuan-Liupanshan fault in the northeastern margin of the Tibetan Plateau. Along the Haiyuan fault, the results show 3.2–6.2 mm/yr of left-lateral strike-slip, and the strike-slip transformed into thrusting deformation (2.8–3.5 mm/yr) along the Liupanshan fault. The results suggest full coupling down to ∼10 km along the Haiyuan fault. Significant portions from the Jinqianghe fault to the Maomaoshan fault, however, are locked to ∼23 km depth. The abrupt change in fault coupling coincides with the Tianzhu seismic gap. High slip rate deficit (3.0–4.5 mm/yr) and seismic moment accumulation rate are also interpreted along the seismic gap. The Liupanshan fault is locked to 15–20 km depth. A high seismic moment accumulation rate, low slip rate deficit (2.0–3.2 mm/yr) and scarcity of large seismic event over the last 1400 years may imply high strain accumulation on the fault. The results demonstrate that crustal deformation in the northeastern Tibetan Plateau is dominated by elastic block rotation. We also conclude that deformation due to fault coupling is limited along the near-field (less than 50 km with velocities less than 2 mm/yr) of the Haiyuan-Liupanshan fault. [ABSTRACT FROM AUTHOR]

Using SAR interferometry on C band Envisat descending track and L band ALOS ascending track SAR images, respectively, we firstly obtain two coseismic deformation fields and one postseismic deformation of the 2010 Yushu earthquake, Tibet, China. In the meanwhile, we also obtain the azimuthal coseismic deformation of the Yushu event by Multi Aperture Interferometry (MAI) technique. With the 3 components of one-dimensional coseismic InSAR measurements, we resolve the complete 3-dimensional deformation of the 2010 Yushu event, which shows conformity and complexity to left lateral slip mechanism. The horizontal deformation is basically consistent with a sinistral slip event; whereas the vertical displacement does show certain level of complexity, which we argue is indicative of local fault geometry variation. Based on the InSAR data and elastic dislocation assumption, we invert for coseismic fault slip and early after-slip of the Yushu event. Our inversion results show major coseismic left lateral strike slip with only minor thrust component. The after-slip model fills most of the slip gaps left by the coseismic fault slip and finds a complementary slip distribution to the coseismic fault slip, which is a good indicator that future earthquake potential on the Yushu segment has been significantly reduced. [ABSTRACT FROM AUTHOR]

We invert InSAR from ascending and descending passes of ENVISAT/ASAR data to estimate the distribution of coseismic slip for the 2008 Mw 7.2 Yutian earthquake by separately using uniform finite element models (FEMs), non-uniform FEMs, and analytical elastic half-space (Okada) models. We use the Steepest Descent method and Laplacian smoothing to regularize and estimate the slip distribution for each of these different models. Fault surface ruptures interpolated from Quickbird optical images constrain the strike of the fault. The uniform FEM and Okada models assume a fault embedded in a homogeneous and elastic finite domain or half-space, respectively, while the non-uniform FEM describes a domain with spatially variable elastic properties according to the geological structure of Tibet plateau and Tarim. We find that the estimated slip distribution of the non-uniform FEM is similar to that of the uniform FEM but has a slightly larger moment (~2 %). However, more noticeable discrepancies occur between the slip distributions of the uniform FEM and Okada model, where the slip pattern estimated with the FEM looks more scattered and locates at a lower depth. We conclude that, for the case of the Yutian earthquake, the improvement brought by simulating the distribution of geologic material properties is insignificant. This suggests that the uniform geologic structure of Okada models may be sufficient for simulating large intra-continental earthquakes. [ABSTRACT FROM AUTHOR]

Abstract: We applied a joint inversion of teleseismic data and InSAR measurements, aiming at obtaining a robust rupture process and slip distribution of the 2010 Yushu earthquake, through reducing the trade-off between slip timing and location. With the condition of the final static deformation of InSAR measurements, the rupture history can be achieved with more confidence. Surface offset constraint is also applied to our inversion. Through synthetic data inversion we find that InSAR data can resolve better fault slip at near surface depth range of 0–15km than teleseismic data; yet it has much poorer resolution at near hypocenter area, where teleseismic data has a relatively better resolution. Our final joint inversion results show that the 2010 Yushu earthquake has an essential rupture time around 20s, during which 90% of the seismic moment has been released. Two peak energy releasing moments occur at 8 and 12s after the earthquake initiation, respectively. The fault slip breaks the surface at most segments of the total 76km fault model and the maximum slip reaches 2.2m. We also find two peak-slip asperities, one at near hypocentral area, mainly distributed at depth around 10–15km; the other distributed in a large area at the eastern segments at near surface. Our most robust estimation of the seismic moment reaches 2.30×1019 N.m, equivalent to a moment magnitude of M w 6.9. [Copyright &y& Elsevier]

The 5 September 2022 Mw 6.7 Luding earthquake occurred on the Moxi segment of the highly active Xianshuihe fault in eastern Tibet. Here, we constrain the coseismic slip by jointly inverting the coseismic displacements measured by Global Positioning System, seismometer and Interferometric Synthetic Aperture Radar. Along the Moxi fault, concentrated left‐lateral strike slip extends ∼30 km along the strike above 10 km depth, producing 0.7–1.0 m shallow slip. Clustered aftershocks and slip inversions suggest that the secondary conjugate Mozigou fault may also involve the rupture. Southward rupture propagation is likely arrested by the barrier‐like fault segments of the Xianshuihe fault, characterized by high interseismic coupling (>0.6) and reduced shear stress rate (<1 kPa/yr) due to interactions with surrounding large locked asperities. The distribution of aftershocks is highly correlated with the positive coseismic Coulomb failure stress changes , which bring the adjacent asperities on the Anninghe and Daliangshan faults ∼0.2 MPa closer to failure. Plain Language Summary: The Mw 6.7 Luding earthquake occurred in eastern Tibet on 5 September 2022, causing 93 casualties and 25 people missing, and substantial damage in the epicentral regions. In this study, we derive the static coseismic displacements from Global Positioning System, seismometer and Interferometric Synthetic Aperture Radar observations. We jointly invert these datasets for the coseismic slip distribution. Our results reveal that the coseismic slip is likely distributed on the conjugate Moxi and Mozigou faults, with the primary slip concentrated on the Moxi fault. We compare the coseismic slip with both the interseismic coupling distribution and shear stress rate on the Moxi fault. We find that the southward coseismic rupture termination is spatially coincident with the areas undergoing a low pre‐earthquake shear stress rate and a high interseismic coupling. The comparison suggests that the coseismic slip during the Luding earthquake is likely controlled by pre‐earthquake interactions with surrounding locked zones on the Anninghe and Daliangshan faults (DLS). We calculate the coseismic Coulomb failure stress changes and demonstrate that the aftershocks could be explained by the induced positive stress perturbations. We conclude that the Luding mainshock enhances earthquake hazards on the Anninghe and DLS. Key Points: Global Positioning System, seismometer, and Interferometric Synthetic Aperture Radar measurements reveal concentrated coseismic slip extending 30 km on the Moxi fault at 0–10 km depthInteractions with large locked asperities lower the shear stress rate to the south of earthquake and lead to the arrest of coseismic ruptureThe Luding event brought the Anninghe and Daliangshan faults ∼0.2 MPa closer to failure, representing enhanced earthquake hazard [ABSTRACT FROM AUTHOR]

The spatial distribution and characterization of a heavily damaged area can be determined by studying surface ruptures of seismogenic faults. If the distribution of surface ruptures can be obtained shortly after they occur, then areas heavily damaged by an earthquake can be readily identified. The information can then be used as a guide for earthquake relief programs. In this paper, an intensity offset-tracking method applied to an ALOS PALSAR image is used to map the Yushu earthquake rupture and to identify the faults activated by the earthquake. Azimuthal displacement analysis indicates that the surface rupture is about 55 km long, running from the epicenter to the southeast, trending N310°W, with a relative displacement of ∼1 m characterized by sinistral slip. The result of range displacement observations indicates that the north wall of the fault is dominated by decreases (i.e., uplift in line of sight observations), whereas in the south wall of the fault, the range displacement is dominated by increases (drops in line of sight observations). Given the position from which the images were recorded, this means that the north wall moves westward, and the south wall move eastward, i.e., left-lateral slip motion across the fault. Finally, an earthquake disaster assessment using computer-assisted image analysis software shows that buildings near the fault rupture have been destroyed most heavily; therefore, the shape of the heavily damage belt is controlled partially by the fault rupture's geometry and the damage degree relates to the magnitude of displacement field. [ABSTRACT FROM AUTHOR]

Abstract: On 21 March 2008, an Ms7.3 earthquake occurred at Yutian County, Xinjiang Uygur Autonomous Region, which is in the same year as 2008 Mw 7.9 Wenchuan earthquake. These two earthquakes both took place in the Bayar Har block, while Yutian earthquake is located in the west edge and Wenchuan earthquake is in the east. The research on source characteristics of Yutian earthquake can serve to better understand Wenchuan earthquake mechanism. We attempt to reveal the features of the causative fault of Yutian shock and its co-seismic deformation field by a sensitivity-based iterative fitting (SBIF) method. Our work is based on analysis and interpretation to high-resolution satellite (Quickbird) images as well as D-InSAR data from the satellite Envisat ASAR, in conjunction with the analysis of seismicity, focal mechanism solutions and active tectonics in this region. The result shows that the 22km long, nearly NS trending surface rupture zone by this event lies on a range-front alluvial platform in the Qira County. It is characterized by distinct linear traces and a simple structure with 1–3m-wide individual seams and maximum 6.5m width of a collapse fracture. Along the rupture zone are seen many secondary fractures and fault-bounded blocks by collapse, exhibiting remarkable extension. The co-seismic deformation affected a big range 100km×40km. D-InSAR analysis indicates that the interferometric deformation field is dominated by extensional faulting with a small strike-slip component. Along the causative fault, the western wall fell down and the eastern wall, that is the active unit, rose up, both with westerly vergence. The maximum subsidence displacement is ∼2.6m in the LOS, and the maximum uplift is 1.2m. The maximum relative vertical dislocation reaches 4.1m, which is 10km distant from the starting rupture point to south. The 42km-long seismogenic fault in the subsurface extends in NS direction as an arc, and it dipping angle changes from 70° near the surface to 52° at depth ∼10km. The slip on the fault plane is concentrated in the depth range 0–8km, forming a belt of length 30km along strike on the fault plane. There are three areas of concentrating slip, in which the largest slip is 10.5m located at the area 10km distant from the initial point of the rupture. [Copyright &y& Elsevier]

Focuses on the use of Synthetic Aperture Radar Interferometry (InSAR) in obtaining digital elevation data in different terrain and physiognomy regions. Relationship between elevation and phase difference; Selection of three different physiognomies and terrains as study regions; Feasibility of using InSAR in areas of with no vegetation and dry ground.

[Display omitted] • AIRS data revealed pre-seismic surface temperature anomalies. • The high commission rate reduces the practicality of earthquake prediction. • The accuracy and omission rates of forecasting in Eastern Japan reached their peaks. • It's reasonable to apply this framework for assessing earthquake forecasting ability. There are growing observational evidences that various geophysical anomalies precede large earthquakes. However, the reliability of these anomalies for earthquake forecasting is controversial, and therefore more consistent assessment of forecasting ability is required. A framework for investigating pre-seismic anomaly detection using essential statistical indicators before global earthquakes is proposed. Surface temperature (ST) data from the Atmospheric Infrared Sounder (AIRS) sensor were used to realize this framework. First, seismic-related ST anomalies were identified, and then the statistical characteristics of forecasting ability for three indicators (accuracy, missed detection, and false alarm) were calculated in retrospective and prospective ways. The ST anomalies displayed some aggregation effects. Negative anomalies mainly concentrated on epicenters and to the north, while positive anomalies were found on the periphery; neither were strongly dependent on earthquake magnitude. The temporal evolution of forecasting metrics was relatively stable for the period 2010–2018. Mean accuracy, missed detection, and false alarm ratios were 6.01%, 1.60%, and 92.39%, respectively. Accuracy and missed detection ratios showed some spatial correlation and both peaked in the same area (e.g., eastern Japan); however, most areas showed very high false alarm ratios. Based on our findings, the combination of AIRS ST data and the Z-score anomaly detection algorithm to predict short-term earthquakes is currently not practical; the possibility of earthquake forecasting based on satellite thermal infrared measurements remains a huge challenge. However, our results confirmed the efficiency of this framework for quantitatively evaluating earthquake forecasting ability. This approach could be applied to various geophysical parameters and anomaly detection methods. [ABSTRACT FROM AUTHOR]

A strong motion seismometer is a kind of inertial sensor, and it can record middle- to high-frequency ground accelerations. The double-integration from acceleration to displacement amplifies errors caused by tilt, rotation, hysteresis, non-linear instrument response, and noise. This leads to long-period, non-physical baseline drifts in the integrated displacements. GNSS enables the direct observation of the ground displacements, with an accuracy of several millimeters to centimeters and a sample rate of 1 Hz to 50 Hz. Combining GNSS and a strong motion seismometer, one can obtain an accurate displacement series. Typically, a Kalman filter is adopted to integrate GNSS displacements and strong motion accelerations, using the empirical values of noise uncertainty. Considering that there are significantly different errors introduced by the above-mentioned tilt, rotation, hysteresis, and non-linear instrument response at different stations or at different times at the same station, it is inappropriate to employ a fixed noise uncertainty for strong motion accelerations. In this paper, we present a Sage–Husa Kalman filter, where the noise uncertainty of strong motion acceleration is adaptively estimated, to integrate GNSS and strong motion acceleration for obtaining the displacement series. The performance of the proposed method was validated by a shake table simulation experiment and the GNSS/strong motion co-located stations collected during the 2023 Mw 7.8 and Mw 7.6 earthquake doublet in southeast Turkey. The experimental results show that the proposed method enhances the adaptability to the variation of strong motion accelerometer noise level and improves the precision of integrated displacement series. The displacement derived from the proposed method was up to 28% more accurate than those from the Kalman filter in the shake table test, and the correlation coefficient with respect to the references arrived at 0.99. The application to the earthquake event shows that the proposed method can capture seismic waveforms at a promotion of 46% and 23% in the horizontal and vertical directions, respectively, compared with the results of the Kalman filter. [ABSTRACT FROM AUTHOR]

Shear strength of damage zone, representing the stress threshold for rupture initiation, is a critical parameter in faulting mechanics. Despite its significance, the damage‐zone's shear strength has not been estimated in natural earthquake ruptures. Here we employed coseismic deformation and strain, kinematic slip model, and finite element modeling to determine the elastic properties and peak shear stress of coseismic damage zones along the 2021 Mw 7.4 Maduo earthquake. Through the analysis of the lowest shear stress resulting in surface ruptures and the highest stress without surface rupture, we constrained the strength within a range of 7–17 MPa. Our result is consistent with strength (5–16 MPa) of sandstone samples from laboratory tests, demonstrating the validity of this estimation. Although factors such as fault maturity and confining pressure influence strength variation, the strength can directly reflect the stress threshold required for macroscopic surface rupture formation in fault damage zones dominated by sandstone. Plain Language Summary: Shear strength of a fault damage zone inform us its ability to withstand shear stress before surface rupture occurs. This information often provides insights into the earthquakes rupture hazards near the Earth's surface. Natural earthquakes provide a perfect opportunity for us to address the question of "what is the shear strength of a fault damage zone." To this end, we require a set of comparative references: the utmost shear stress that a damage zone can withstand without surface rupture and the minimal one necessary for the zone to generate surface rupture. In this investigation, we studied the 2021 Maduo earthquake because it had multiple distinct surface rupture segments in some places but not in others. By comparing shear stress in these segments, we could figure out how strong the damage zones were. We used real observations and numerical simulations to estimate how much stress these damage zones can withstand. Our results indicate that, for damage zones embedded in intact sandstone with a Young's modulus of 15 GPa, its strength ranges from 7 to 17 MPa. Significantly, our study is the first to reveal the shear strength of damage zone from a nutural earthquake using geodetic data. Key Points: Shear modulus and shear stress for coseismic damage zone were determined using coseismic deformation and finite element modelEstimated shear strength from overlap range of coseismic stress between surface rupture and non‐rupture parts ranges from 7 to 17 MPaShear strength exhibits positive and negative correlations with confining pressure and fault maturity, respectively [ABSTRACT FROM AUTHOR]

As the northern boundary of the Tibetan Plateau, the long Altyn Tagh fault (ATF) controls the regional tectonic environment, and the study of its long-term fault slip rate is key to understanding the tectonic evolution and deformation of the northern Tibetan Plateau. In this paper, we measure the fault slip rate of the western segment of the ATF using InSAR observations between 2015 to 2020. The Multi-Temporal Interferometric InSAR analysis is applied to obtain the two-dimensional fault-parallel and vertical displacement fields. The spatially dense InSAR observations clearly illustrate the asymmetrical pattern of displacement fields across the fault. Constrained by our InSAR observations, the fault slip rate and locking depth of the western segment of the ATF are inverted using four different models in a Bayesian framework. The two-layer viscoelastic model incorporating lateral heterogeneity of rheology in the lower crust indicates that the fault slip rate of the western ATF is estimated to be 9.8 ± 1.1 mm/yr (at 83.8°E across the ATF) and 8.6 ± 1.1 mm/yr (at 85.1°E), respectively, and the locking depth is 15.8 ± 4.3 km and 14.8 ± 4.9 km. Our new estimates generally agree with the previous estimates of fault slip rate constrained by GPS observations. We conclude that the contrast between the thickness of the elastic layer and the shear modulus of the Tibetan plateau and the Tarim basin jointly contribute to the asymmetric interseismic strain accumulation on the ATF. [ABSTRACT FROM AUTHOR]

We derived the joint slip models of the three major events in the 2016 Central Italy earthquake sequence by inverting strong-motion and InSAR datasets. b-values and the historic earthquake scarp offset were also investigated after processing the earthquake catalog and near-field digital elevation model data. The three major events gradually released seismic moments of 1.6 × 1018 Nm (Mw 6.1), 1.5 × 1018 Nm (Mw 6.1), and 1.1 × 1019 Nm (Mw 6.7), respectively. All the ruptures exhibit both updip and along-strike directivity, but differ in the along-strike propagation direction. The high b-value found beneath three mainshock hypocenters suggests possible fluid intrusions, explaining the cascading earthquake behavior. The cumulative surface scarp from past earthquakes shows rupturing features that are consistent with the 2016 earthquake sequence, suggesting a characteristic fault behavior. Under the assumption of the Gutenberg–Richter law, the slip budget closure test gives a maximum magnitude of Mw 6.7 and implies the seismic hazard from the largest event has been released in this sequence. [ABSTRACT FROM AUTHOR]

The article discusses the research on surface rupture zone of seismic faults in China using high-resolution remote sensing technology and mentions the quantitative research of earthquakes. Topics include the survey on geometric and kinematic characteristics of the rupture zone; the 7.2 earthquake at Yutian in Xinjiang, China; and the use of digital elevation model the surface rupture behavior.

Taking the 2017 Mw6.5 Jiuzhaigou earthquake as a case study, ionospheric disturbances (i.e., total electron content and TEC) and thermal infrared (TIR) anomalies were simultaneously investigated. The characteristics of the temperature of brightness blackbody (TBB), medium-wave infrared brightness (MIB), and outgoing longwave radiation (OLR) were extracted and compared with the characteristics of ionospheric TEC. We observed different relationships among the three types of TIR radiation according to seismic or aseismic conditions. A wide range of positive TEC anomalies occurred southern to the epicenter. The area to the south of the Huarong mountain fracture, which contained the maximum TEC anomaly amplitudes, overlapped one of the regions with notable TIR anomalies. We observed three stages of increasing TIR radiation, with ionospheric TEC anomalies appearing after each stage, for the first time. There was also high spatial correspondence between both TIR and TEC anomalies and the regional geological structure. Together with the time series data, these results suggest that TEC anomaly genesis might be related to increasing TIR. [ABSTRACT FROM AUTHOR]

We combine Global Positioning System (GPS) velocity field with Interferometric Synthetic Aperture Radar (InSAR) results to study the interseismic deformation across the western Altyn Tagh fault (ATF) at longitude 86° E. GPS and InSAR data are consistent after correcting for the contribution from vertical deformation in the InSAR line-of-sight map. InSAR and GPS data identify an area of ∼2 mm a−1 sinistral shear and ∼6 mm a−1 of NS shortening located ∼150 km south of the ATF near the Manyi fault system. Excluding the data located in that area, Bayesian inversion of a 2-D profile across the ATF indicate a locking depth of 14.8 ± 3.5 km and a slip rate of 8.0 ± 0.4 mm a−1, lying at the lower range of previously published estimates. In addition, we find no significant offset between the fault at depth and the surface fault trace and no asymmetry of the interseismic profile that implicitly reveal lateral variations of the elastic strength across the ATF. Detailed analysis of InSAR profile across the fault show no distinguishable surface creep along the western ATF. Our study highlights how different data sets, data selection and model assumption might impact results on the ATF slip rate, locking depth and rheological contrast across the fault. [ABSTRACT FROM AUTHOR]

This study focuses on the crustal deformation and interseismic fault coupling along the strike-slip Kunlun fault, northern Tibet, whose western segment ruptured in the 2001 Mw 7.8 Kokoxili earthquake. We first integrated published Global Positioning System (GPS) velocity solutions and calculated strain rate fields covering the Kunlun fault. Our results show abnormally high post-earthquake strain rate values across the ruptures; furthermore, these exceed those in pre-earthquake data. Together with two tracks of interferometric synthetic aperture radar (InSAR) observations (2003–2010) and position time-series data from two continuous GPS sites, we show that the postseismic deformation of the Kokoxili earthquake may continue up to 2014; and that the postseismic transients of the earthquake affect the 2001–2014 GPS velocity solutions. We then processed the GPS data observed in 2014–2017 and obtained a dense interseismic velocity field for the northern Tibet. Using a fault dislocation model in a Bayesian framework, we estimated the slip rates and fault coupling on the Kunlun fault in 1991–2001 and 2014–2017. Results show an increase of slip rates and eastward migration of high fault coupling on the Kunlun fault after 2001. We propose the temporal variations are a result of the eastward accelerating movement, as a whole, of the Bayanhar block, whose boundaries were decoupled by several large earthquakes since 1997. Moreover, our results show the accumulated elastic strains along the Alake Lake-Tuosuo Lake segments could be balanced by an Mw 7.4–7.7 earthquake. [ABSTRACT FROM AUTHOR]

We use Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar (InSAR) data to quantify the interseismic coupling along the Xianshuihe‐Xiaojiang fault system (XXFS). Our results confirm 7–11 mm/a of left‐lateral strike‐slip motion along the XXFS. South of the Shimian County, high interseismic coupling is found down to 20 km‐depth along the Anninghe‐Zemuhe fault. A second 130 km‐long section, highly coupled down to 10–15 km‐depth, is identified along the northern Xiaojiang fault. North of the Kangding County, the Xianshuihe fault appears to be predominantly creeping at depth. Previously identified long‐lasting surface creep following the 1973 M 7.6 Luhuo earthquake has now ended, possibly marking the initiation of a new period of stress accumulation and progressive relocking of the fault. Farther south along the fault, we identify a ∼30 km‐long section (∼30.2°–30.4°N) creeping at ∼7 mm/yr. Furthermore, creep accelerated by ∼2 mm/yr during the 2008–2014 period, possibly as a result of the static Coulomb failure stress increment induced by the 2008 Mw 7.9 Wenchuan earthquake, that occurred ∼200 km northeast of it. The 2014 Mw 5.9 Kangding earthquake occurred along the creeping section at the end of the period of accelerated creep. Finally, taking advantage of the long historical records of past large earthquakes available since 1327, we quantitatively compare the moment budget along the XXFS. Results highlight that along the Anninghe, Zemuhe, and Xiaojiang faults, moment deficit for Mw > 7 earthquakes has now accumulated. Plain Language Summary: More probable location and size of future large earthquakes can be anticipated if we accurately know the present‐day distribution of interseismic fault coupling. Knowing date of the latest large earthquake, one can infer slip deficit and estimate the energy available for forthcoming seismic ruptures. Here, we focus on the Xianshuihe‐Anninghe‐Zemuhe‐Xiaojiang fault system in Tibet, which is the largest strike‐slip fault crossing southwestern China. Historical records attest numerous earthquakes during the last few centuries, but the fault has remained relatively silent in the past few decades. We use GPS data spanning nearly 20 years to investigate the present‐day interseismic fault coupling along the fault system. We find the Anninghe, Zemuhe, and Xiaojiang faults have accumulated enough energy for several Mw > 7 earthquakes. However, our results show that the northern segment called the Xianshuihe fault is not locked, suggesting a relatively low energy accumulation rate and seismic potential, despite a long history of frequent earthquakes. Finally, we identify a new ∼30 km‐long creeping section on the Xianshuihe fault, whose creep rate show a transient behavior, possibly induced by the effect a ∼200 km remote earthquake that occurred in 2008. Key Points: Highly heterogeneous interseismic coupling along the Xianshuihe‐Xiaojiang fault system, with the northern Xianshuihe fault predominantly creeping at depthSurface creep observed during 4 decades after the 1973 M 7.6 Luhuo earthquake has endedA new ∼30 km‐long transient creeping segment on the Xianshuihe fault, with creep rate increasing after the 2008 Mw 7.9 Wenchuan earthquake [ABSTRACT FROM AUTHOR]

The coseismic deformation field on both sides of the fault, especially the distribution and change characteristics of near-field deformation, not only provides important constraints for the fine inversion of the slip distribution model but also serves as an important basis for the anti-disruption defense of the cross-fault linear engineering facilities. In this paper, we used Sentinel-1 satellite data to obtain the coseismic deformation field of the Maduo earthquake by using InSAR and offset techniques. We quantitatively compared the coseismic displacement of the three types of data: InSAR, offset, and optical images. The results show that optical images and offset provided more robust near-fault (<2 km) deformation insights than InSAR, which exhibited irregular deformation patterns due to incoherence near the fault. The maximum relative displacements for InSAR and offset observations are ~2.8 m and 4 m, respectively. Then we tested various fault slip models with different data constraints, revealing that a combined inversion of GPS, InSAR, and offset data offers superior constraints on slip distribution. This integrative approach effectively captured both shallow and deep fault slip, particularly near the fault zone. The eastern branch fault model, jointly constrained by GPS, InSAR, and offset data, is the optimal coseismic slip distribution model for the Maduo earthquake, and the maximum slip is 5.55 m. [ABSTRACT FROM AUTHOR]

We investigated the spatial distribution of aseismic creep on the Laohushan‐Haiyuan fault using Global Positioning System (GPS) data (1999–2017) and Interferometric Synthetic Aperture Radar (InSAR) data (2003–2010). Comparisons among GPS, InSAR line‐of‐sight (LOS) rates, and leveling show that neither leveling nor GPS vertical velocities can fit the vertical signal mapped into the LOS, implying either complicated vertical crustal deformation in northeastern Tibet and/or complex error structures in the InSAR data. Thus, we combined horizontal GPS with high‐pass filtered InSAR data to produce a continuous LOS rate map crossing the fault. Our geodetic data reveal three creep sections along the fault. Both the restored LOS data and decomposed ascending and descending InSAR data highlight the fact that vertical motion can cause an overestimation of creep rate; we obtained a refined creep rate of 2.5 ± 0.4 mm/a on the Laohushan fault. We further identified a 10 km‐long, ∼3–5 mm/a creep section (∼104.2°E−104.3°E) and a 43 km‐long, ∼1–3 mm/a creep section (∼105.3°E−105.7°E) on the western and eastern Haiyuan fault respectively. Both are located on fault sections that ruptured during the 1920 M∼8 earthquake, suggesting that the 1920 earthquake was able to cross pre‐existing creep sections or that the fault shows heterogeneous relocking after large earthquakes, with creep lasting decades on some parts of the rupture. Fault coupling shows a highly variable rate of slip deficit accumulation along strike, suggesting that coupling might significantly evolve during the period between two large earthquakes. Plain Language Summary: Ample observations document slip on a fault plane is seismic and/or aseismic; the former situation refers to a fault slips suddenly and produce earthquakes, the latter means a fault slips gradually without generating seismic energy. Studying the fault‐creep activity contributes to better understand the seismic behavior of crustal faults. We focus on the Haiyuan fault system, which is located on the northeastern margin of Tibet. Several large earthquakes occurred along the fault system in the past, including the 1920 M∼8.0 Haiyuan and the 1927 M ∼8–8.3 Gulang earthquakes. Using Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS) data, we identified three creep sections along the fault. Our results show vertical deformation that has mapped into the InSAR line‐of‐sight (LOS) could leads to overestimation of the creep rate. Besides, our findings imply either large earthquakes can rupture creep regions or they can generate long‐lasting afterslip/creep. Key Points: Of three creep segments on the Laohushan‐Haiyuan fault, two are located along fault sections that ruptured during the 1920 M∼8 earthquakeVertical motion leads to overestimation of the creep rate on the Laohushan fault; our refined creep rate is 2.5 ± 0.4 mm/aOur findings support that large earthquakes can rupture pre‐existing creep regions or that they can generate long‐lasting afterslip or creep [ABSTRACT FROM AUTHOR]

We investigate the coseismic and post-seismic deformation due to the 6 February 2018 M w 6.4 Hualien earthquake to gain improved insights into the fault geometries and complex regional tectonics in this structural transition zone. We generate coseismic deformation fields using ascending and descending Sentinel-1A/B InSAR data and GPS data. Analysis of the aftershocks and InSAR measurements reveal complex multifault rupture during this event. We compare two fault model joint inversions of SAR, GPS and teleseismic body waves data to illuminate the involved seismogenic faults, coseismic slip distributions and rupture processes. Our preferred fault model suggests that both well-known active faults, the dominantly left-lateral Milun and Lingding faults, and previously unrecognized oblique-reverse west-dipping and north-dipping detachment faults, ruptured during this event. The maximum slip of ∼1.6 m occurred on the Milun fault at a depth of ∼2–5 km. We compute post-seismic displacement time series using the persistent scatterer method. The post-seismic range-change fields reveal large surface displacements mainly in the near-field of the Milun fault. Kinematic inversions constrained by cumulative InSAR displacements along two tracks indicate that the afterslip occurred on the Milun and Lingding faults and the west-dipping fault just to the east. The maximum cumulative afterslip of 0.4–0.6 m occurred along the Milun fault within ∼7 months of the main shock. The main shock-induced static Coulomb stress changes may have played an important role in driving the afterslip adjacent to coseismic high-slip zones on the Milun, Lingding and west-dipping faults. [ABSTRACT FROM AUTHOR]

The Türkiye-Syria earthquake doublet of 6 February 2023 (Mw 7.8 at 01:17 UTC and Mw 7.6 at 10:24 UTC) resulted in extensive damage and tens of thousands of casualties. We present the surface displacements of the two earthquakes from synthetic aperture radar (SAR) offset tracking measurements. We extracted the geometric parameters of the rupture faults from the surface displacements and early aftershock distribution, based on which we inverted the coseismic slip distributions. We then calculated Coulomb stress to investigate the triggering relationship between the earthquakes and stress transfer to neighbouring faults and regions. The coseismic ruptures of the earthquake doublet were predominantly left-lateral strike-slip motions distributed between 0 and 15 km depth. The maximum fault slip reached > 8 m (Mw 7.8) and almost 10 m (Mw 7.6). The coseismic deformation and fault slip motion are consistent with the overall westward extrusion of the Anatolian Plate relative to the Eurasian and Arabian plates. The Mw 7.8 earthquake increased Coulomb failure stress at the hypocenter of the Mw 7.6 earthquake, implying that the Mw 7.8 event had a strong positive causative effect. Moreover, coseismic stress perturbations revealed a positive Coulomb stress effect on the middle Puturge Fault, northern Dead Sea Fault Zone (DSFZ), Yesemek Fault, Antakya Fault, and Turkoglu Fault, indicating an increasing seismic hazard in these regions. [ABSTRACT FROM AUTHOR]

The 2021 Mw 7.4 Maduo earthquake (eastern Tibetan Plateau), with a surface deformation zone that is characterized by multiple distinct surface rupture segments separated by segments with no surface ruptures, provides an excellent opportunity to constrain the strain threshold (minimum strain) for the formation of surface ruptures. Through sub‐pixel correlation of pre‐ and post‐event SPOT‐6/7 satellite images (pixel = 1.5 m), we derive surface displacement and shear‐strain fields of the Maduo earthquake. By quantifying the minimum shear strain along the surface rupture segments and the maximum shear strain along the segments with no surface ruptures, we estimate a rupture strain threshold, which ranges from 0.8% to 1.8% and appears to decrease with structural maturity of the causative fault. This threshold is generally consistent with a rupture strain limit (0.5%–1.5%) of intact rocks from laboratory measurements, and is higher than a commonly‐assumed strain threshold (0.5%) for inelastic deformation. Plain Language Summary: During large earthquakes, two sides of the fault dislocate relative to each other and create a narrow surface deformation zone in between. Whether surface rupture can be produced depends on the strain, a unitless measurement of deformation, within the surface deformation zone. For natural earthquakes, the minimum strain that is needed to create surface rupture (the rupture strain threshold) has not yet been determined. Here we investigate such strain for the surface rupture of 2021 Mw 7.4 Maduo earthquake at the eastern Tibetan Plateau. By analyzing the Maduo coseismic surface strain field and its correlation with spatial extent of the coseismic surface rupture, we determine that the rupture strain threshold ranges from 0.8% to 1.8%, and it appears to decrease with structural maturity of the causative fault. This threshold is generally consistent with a rupture strain limit (0.5%–1.5%) of intact rocks from laboratory measurements, and is higher than a commonly‐assumed strain threshold (0.5%) for inelastic deformation. Our study, for the first time, provides a direct measurement of rupture strain threshold based on natural earthquake studies. Key Points: We derive the 2021 Maduo coseismic surface displacement and strain fields through sub‐pixel correlation of SPOT‐6/7 optical imagesThe strain threshold for the formation of Maduo surface rupture is 0.8%–1.8%, and it appears to decrease with the causative fault maturityThis value is consistent with a laboratory‐derived rupture threshold 0.5%–1.5% and is higher than a commonly‐assumed inelastic threshold 0.5% [ABSTRACT FROM AUTHOR]

The ascending and descending interferometric synthetic aperture radar data are used to investigate the fault rupture and slip model of the 2018 Mw 7.5 Sulawesi, Indonesia, earthquake. The best fitting slip model indicates that this earthquake ruptured not only a segment extending the Palu fault to the north but also a northwestern segment offshore. The slip on the onshore fault is predominant left‐lateral strike slip. The slip on the offshore fault is dominated by normal faulting with a maximum slip of ~6.3 m. The newly discovered offshore normal faulting is likely to be the cause of the tsunami after the shock. Combined with previous geomorphic, tectonic, geodetic, and modeling studies, we suggest that the kinematics of the Palu fault maintained the same style of faulting from north to south, which resulted from an oblique extension occurred on an east dipping fault at depths. The deformation pattern of NW Sulawesi is dominated by this slip mechanism. Plain Language Summary: On 28 September 2018, a large earthquake (Mw 7.5) struck the Palu city, Sulawesi Island, eastern Indonesia, where three tectonic plates, the Indo‐Australian, Philippine Sea and Sunda Plates, interact. Using a joint analysis of ALOS‐2 and Landsat 8 satellite images, we found coseismic slip partitioning during this earthquake. We propose a conceptual structural model to interpret the slip mechanism, where oblique extension occurred on an east dipping fault at depth and the slip is partitioned into left‐lateral strike‐slip motion at the ground surface and normal faulting at the seabed surface. A similar small‐scale geological system is observed in the Palu basin along the same fault system, which provides good evidence for this mechanism. Our study has important implications for understanding the geohazards and deformation mechanism of faults in eastern Indonesia. Key Points: This earthquake ruptured one onshore segment with predominant strike‐slip motion and an offshore fault with a significant normal slipOffshore normal faulting is responsible for the tsunami after the earthquakeThis event revealed the kinematics of the Palu fault, which can be used to interpret the geomorphological expression in the Palu basin [ABSTRACT FROM AUTHOR]

Abstract The Mw6.9 Yushu earthquake occurred on 14 April 2010, in Qinghai, China; the largest aftershock, a Mw5.8 event, occurred west of the mainshock on 29 May 2010 (˜40 days later). The aftershock had a different focal mechanism from the mainshock. Furthermore, seismicity after 29 May showed different spatial characteristics in terms of focal depth and distribution direction. To better understand the faulting and the relationship between these two events, we derived the whole displacement field caused by the Yushu mainshock and the Mw5.8 aftershock based on multi-perspective Interferometric Synthetic Aperture Radar (InSAR) data. We then conducted a robust inversion of the slip distribution jointly constrained by InSAR and GPS data. The results indicate that the Mw5.8 aftershock produced a separated deformation field with significant displacement changes of up to ˜4–6 cm, which indicates another intersecting ruptured fault at the west end of the Yushu fault. The slip distribution shows a 75-km NW rupture with a maximum slip of ˜2.1 m at a depth of ˜0–10 km on the main Yushu fault, and a 20 km NE rupture with peak slip of ˜0.4 m at a depth of ˜5–15 km on a vertical hidden fault. Both events showed a dominant left-lateral component. The total rupture length associated with the 2010 Yushu earthquake sequence reached ˜95 km. By calculating Coulomb stress changes, we confirmed that the mainshock triggered the Mw5.8 aftershock. Our results imply that the increased stress at the western end of the Yushu fault caused by the mainshock rupture may have played an important role in transferring the rupture plane from the Yushu fault to the NE hidden fault. [ABSTRACT FROM AUTHOR]

Strain partitioning between strike-slip faults in mountains and thrust faults in the foreland is a typical mountain building process to accommodate oblique plate convergence. Studying the geometry and movement of such strain-partitioned fault systems is key to understanding the mountain building process and related seismic hazards. The 2022 M w 6.6 Menyuan earthquake is the largest strike-slip earthquake to have ruptured the northern Qilian Shan fault system in the modern geodetic era. We combined satellite and field observations to determine the fault geometry and coseismic slip distribution in the Menyuan earthquake, and link the distribution of coseismic slip with the pattern of interseismic strain accumulation within the northern Qilian Shan from our geodetic slip model. We find that the Menyuan earthquake ruptured a 25 km-long section of the left-lateral Longlongling Fault between the surface and 7 km depth. The maximum slip was 4 m at 3–4 km depth. Damage to a high-speed railway tunnel recorded a fault offset of 2.7 m at a depth of 200 m compared to 2.5–3.0 m on the surface, suggesting that dispersion of the rupture through unconsolidated shallow sediments was limited, at least at the tunnel site. We also determined the pattern of interseismic deformation prior to the earthquake using Interferometric Synthetic Aperture Radar and Global Navigation Satellite System data. We found the interseismic geodetic data can be explained by the oblique movement of a low-angle décollement beneath the Qilian Shan rather than a strain partitioning fault system. We suggested that the strike-slip faults and foreland thrusts are separated by a creeping décollement, which would act as a barrier to stop the cascading rupture of the strike-slip and thrust fault. [ABSTRACT FROM AUTHOR]

We have generated a more than 500 km long postseismic deformation rate map and cumulative displacement time series in the central Bayankala Block of the Tibetan Plateau using ENVISAT/ASAR data from 2003 to 2010 by the π-RATE stacking algorithm. This rate map spans a period of ∼7.2 years and reveals that postseismic motion of 2001 Kokoxili earthquake exhibits a striking signal, dominating crustal deformation of the central Bayankala Block with a cross-fault magnitude ∼9–11 mm/yr in line of sight (LOS) (∼93.1°E). The southern and northern parts of the postseismic deformation field exhibit different patterns and variable magnitudes, reflecting asymmetry of the displacement distribution. Postseismic motion affects eastward extrusion of the central Bayankala Block, which reaches ∼15 km north of the Ganzi-Yushu fault. To further investigate viscoelastic relaxation of the lower crust in this area after approximately two years, E-M, E-M-M and E-S models are constructed. The result shows that the best fit viscosity for the lower crust is about 1 × 10 19  Pa·s. Comparison between cumulative displacements resolved by these three models shows that viscoelastic relaxation of the lower crust makes the most significant contribution to postseismic stress relaxation after 2001 event. [ABSTRACT FROM AUTHOR]

On 8 August 2017, a Ms7.0 earthquake stroke the city of Jiuzhaigou, Sichuan, China. The Jiuzhaigou earthquake occurred on a buried fault in the vicinity of three well-known active faults and this event offers a unique opportunity to study tectonic structures in the epicentral region and stress transferring. Here we present coseismic displacement field maps for this earthquake using descending and ascending Sentinel-1A Interferometric Synthetic Aperture Radar (InSAR) data. Deformation covered an area of approximately 50 × 50 km, with a maximum line-of-sight (LOS) displacement of ~22 cm in ascending and ~14 cm in descending observations on the west side of the source fault. Based on InSAR and Global Positioning System (GPS) measurements, both separately and jointly, we constructed a one-segment model to invert the coseismic slip distribution and dip angle of this event. Our final fault slip model suggests that slip was concentrated at an upper depth of 15 km; there was a maximum slip of ~1.3 m and the rupture was dominated by a left-lateral strike-slip motion. The inverted geodetic moment was approximately 6.75 × 10 18  Nm, corresponding to a moment magnitude of Mw6.5, consistent with seismological results. The calculated static Coulomb stress changes indicate that most aftershocks occurred in stress increasing zones caused by the mainshock rupture; the Jiuzhaigou earthquake has brought the western part of the Tazang fault 0.1–0.4 MPa closer to failure, indicating an increasing seismic hazard in this region. The Coulomb stress changes caused by the 2008 Mw7.8 Wenchuan earthquake suggest that stress loading from this event acted as a trigger for the Jiuzhaigou earthquake. [ABSTRACT FROM AUTHOR]

An earthquake of seismic moment magnitude (Mw)6.3 occurred in southern Taiwan on February 6th, 2016. In this paper, we use differential interferometric synthetic aperture radar (D-InSAR) technology and data from the Sentinel-1A(S-1A)/IW radar satellite to estimate the coseismic deformation of this event. On the basis of ascending data, InSAR analysis reveals a circle-like uplift area, 45 km in diameter, 20 km northwest of the epicenter, that has a maximum line of sight (LOS) displacement of 12 cm. South and east of this uplift area, minor subsidence occurred, while descending data indicated uplift in the west and sink in the east, with maximum values of 8.0 cm and 6.0 cm, respectively. Displacements are continuously distributed across the entire deformed region, implying that the seismogenic fault did not pierce the Earth’s surface. Based on these InSAR displacements coupled with GPS observations, we inverted the distribution of source slip using an elastic half-space fault model both separately and jointly. Four kinds of inversion results consistently show a slip concentration to the northwest of the epicenter, with maximum slip in the range 0.35–0.55 m confined to depths of 6–15 km in the crust. The source fault is dominated by thrust with left-lateral strike slip; joint inversion using InSAR ascending and descending data, as well as GPS data, resulted in maximum slip of 0.44 m, within the range determined by separate inversions. Inversion results show that this was a Mw6.25 earthquake, consistent with measurements collected by the United States Geological Survey (USGS). [ABSTRACT FROM AUTHOR]

The density of GPS measurements is usually one of the key issues in resolving 3-D coseismic deformation field from integrating GPS and Interferometric Synthetic Aperture Radar (InSAR) measurements with pure mathematic interpolation methods An approach that combines the elastic dislocation model with the Best Quadratic Unbiased Estimator (BQUE) or a robust estimation method named IGG (Institute of Geodesy and Geophysics) is proposed to reconstruct 3-D coseismic deformation field, in which only a small amount of GPS data is needed to produce a reasonable initial 3-D coseismic deformation. Then the BQUE and IGG are used to weight the InSAR and GPS measurements to avoid computational issues caused by the negative variance problem and to decrease the impact from gross errors. The Wenchuan earthquake is used to test the proposed method. We find that the developed method makes it possible to use only a few GPSs and InSAR data to recover the 3-D coseismic deformation field, which offers extensive future usage for measuring earthquake deformation, particularly in some tectonically active regions with sparse GPS measurements. [ABSTRACT FROM AUTHOR]

Coseismic deformation fields caused by the moment magnitude (Mw)7.8 mainshock and Mw7.3 aftershock of the 2015 Gorkha-Nepal earthquake are obtained by analyzing Sentinel-1A/IW ascending and descending interferometry data. Results show that the deformation field associated with the Mw7.8 mainshock roughly resembles a prolate ellipse, extending from the epicenter about 20° east by south. The main region of deformation is about 160 km by 110 km, comprising a large southern area of uplift, and a small northern area of subsidence. Assuming that rupture occurred in a homogeneous elastic half-space, the coseismic fault slip models of the mainshock and aftershock are inverted based on a shallow dip fault constrained by the three data sets, Sentinel-1A/IW descending data, ascending data, and ALOS-2 descending data, separately or in combination. Mainshock slip distributions generated from all three data sets are similar, and inversion constrained by all three in combination reveal a comprehensive fault slip model. Indeed, coseismic slip is mainly distributed within a narrow 40 km zone to the north of the Main Frontal Trust (MFT), and at 6–15 km subsurface depth. In addition, the maximum slip in this event was about 5.1 m, the Mw7.8 mainshock ruptured the deep part of the seismogenic zone, while the region between the southern boundary of the rupture area and the MFT remained locked. Therefore, a considerable earthquake risk remains to the south of Kathmandu. The inverted coseismic slip of the Mw7.3 aftershock was concentrated in a small area, close to, and southeast of the epicenter, with maximum displacement of about 3 m. Finally, because there is no overlap between the two slip areas of the mainshock and aftershock, the gap between them, about 15 km in length, has additional potential to generate future earthquakes. [ABSTRACT FROM AUTHOR]

We obtained the rupture process and slip distribution of the 2022 Mw6.6 Menyuan earthquake by jointly inverting accelerogram data and InSAR measurements. The near-field InSAR measurements provide good constraints on the shallow slip distributions (<6 km). The accelerogram data enable us to better resolve the deeper coseismic slip (>6 km). The combination of two types of data provided improved constrains on slip distribution of the 2022 Menyuan earthquake. The results from joint inversion of InSAR and accelerogram data reveal a 26-km-long rupture length, which roughly agrees with the mapped length from the optically identified surface rupture trace and the InSAR deformation field. We imaged a major asperity with a dimension of 14 × 6 km at 4 km depth updip of the hypocenter. The maximum slip is estimated to be 3.8 m at 4 km depth. The duration of the 2022 Menyuan earthquake is ~14 s, and 90% of the seismic moment is released in the first 10 s. The total seismic moment is estimated to be 1.31 × 1 × 1019 N·m, equivalent to a moment magnitude of Mw6.7. Our results highlight that the moderate but shallow rupture during the 2022 Menyuan earthquake could intensify the seismic damage on the surface, confirmed by field investigations. [ABSTRACT FROM AUTHOR]

Normal faulting in orogens is usually associated with gravitational potential energy (GPE) due to elevation differences and contrasting material rheology between mountains and foreland areas. The 2008–2020 Yutian normal earthquake sequence in the West Kunlun Mountains has been closely studied over the last decade, but the mechanisms of normal faulting in the region are still an open question as both gravitational force and releasing step-overs have been proposed to explain the extension in West Kunlun. We investigated the normal faulting mechanism through (1) co- and postseismic kinematic models constrained by interferometric synthetic aperture radar (InSAR) and (2) analytical force balance modeling with a thin viscous approximation. We used the InSAR time series technique to map the postseismic deformation within 3 years following the 2008 Yutian earthquake. The InSAR-derived afterslip model revealed a clear slip polarity transition from strike-slip to normal faulting, suggesting a tight connection between the normal faulting and bounding strike-slip faults. The global navigation satellite system (GNSS)-constrained thin viscous sheet model gave an estimate of ~1022 Pa∙s for the effective viscosity of the West Kunlun lithosphere, which is generally in line with previous estimates that used a similar methodology. The stress field predicted by the thin viscous model primarily features north–south compression rather than east–west extension, suggesting that normal faulting is caused not by the GPE contrast across the margin of West Kunlun, but instead by the releasing step-overs bounded by strike-slip faults. The main contribution of this paper is therefore identifying the primary causes of normal faulting in West Kunlun and improving our understanding of the active tectonics of the NW Tibetan Plateau. • The postseismic InSAR time series following the 2008 Mw 7.1 Yutian earthquake is presented. • The slip transition from strike-slip to normal faulting was revealed by slip models. • Extension in West Kunlun is due to releasing step-overs rather than gravitational force. [ABSTRACT FROM AUTHOR]