Coseismic Slip Distribution Research Articles

Synthetic aperture radar interferometry–based coseismic deformation and slip distribution of the 2022 Menyuan MS6.9 earthquake in Qinghai, China

On January 8, 2022, a 6.9 magnitude earthquake occurred in Menyuan County, Qinghai Province, with the epicenter located at the intersection of the Tuolaishan Fault and the Lenglongling Fault, which are part of the Qilian–Haiyuan fault zone. This study investigated the sliding characteristics and seismic mechanism of the earthquake to understand the activity and seismic risk of the fault on the northeastern margin of the Qinghai–Tibet Plateau. This paper analyzed Sentinel-1 synthetic aperture radar images to obtain the coseismic deformation field of the earthquake, which was then used to invert the slip distribution of the seismogenic fault and the coseismic Coulomb stress on the surrounding faults caused by the earthquake. It was found that the earthquake was primarily characterized by sinistral strike-slip movement. Along the satellite line of sight, the south wall of the fault had a maximum deformation of 0.62 m, and the north wall had a maximum deformation of 0.48 m. The coseismic slip distribution results indicated that the maximum slip of the earthquake was 4.51 m, and the moment magnitude was MW6.7. The Coulomb stress analysis showed that the 2016 Menyuan earthquake promoted the occurrence of the 2022 Menyuan earthquake.

Insight into the 1 December 2016 Mw 6.2 Juliaca Earthquake, Southern Peru, by InSAR Observations and Field Investigation

On 1 December 2016, an Mw 6.2 earthquake characterized by normal faulting occurred in the highlands of the central Andes in southern Peru, marking the region’s largest shallow event. The occurrence of the earthquake provides a significant chance to gain insight into the regional tectonic deformation and the seismogenic mechanism of the shallow normal-faulting earthquake, as well as the regional potential seismic risk. Here, we first utilize Sentinel-1A interferometric synthetic aperture radar (InSAR) data to extract the coseismic and postseismic deformation associated with this earthquake and then determine the detailed coseismic slip and postseismic afterslip distribution of this event. Coseismic modeling results indicate that the coseismic rupture is mainly characterized by normal faulting with some dextral strike-slip components. Most coseismic slip is confined to a depth range of 2–12 km, indicating an obvious slip deficit area in the shallow fault part. Further postseismic modeling reveals that the majority of afterslip is concentrated at depths of 0 to 5.4 km. The relatively shallow postseismic afterslip makes up for the coseismic slip deficit area to some extent. Through a joint analysis of the inversions, seismic data, and regional geology and geomorphology, we infer that the occurrence of this 2016 normal-faulting event is a result of regional gravitational collapse. In addition, we investigate the relationship between the 2016 earthquake and great historical earthquakes near the subduction zone of the central Andes and find that the 2016 event is likely promoted in advance by these events through our calculations of the coseismic and postseismic Coulomb stress changes. Finally, we should pay more attention to the nearby Falla Huaytacucho-Condoroma fault and the western segment of the Vilcanota Fault because of their relatively high stress loading.

High‐Resolution Interseismic Strain Mapping From InSAR Phase‐Gradient Stacking: Application to the North Anatolian Fault With Implications for the Non‐Uniform Strain Distribution Related to Coseismic Slip Distribution

AbstractHigh‐resolution interseismic strain mapping is important for studying faulting behavior and for assessing seismic hazards. Interferometric Synthetic Aperture Radar has been widely applied to measure interseismic deformation along active strike‐slip faults. However, phase unwrapping errors and over‐smoothing effects limit its ability to map the extremely‐high strain due to shallow creep. Here, without the involvement of ground‐based measurement, we perform phase‐gradient stacking on wrapped Sentinel‐1 interferograms to directly map the shear strain rates along the North Anatolian Fault (NAF) with unprecedented resolution. The derived high‐resolution strain‐rate map reveals five strain‐concentrated segments on the NAF, implying shallow creeps. We find that their spatial distribution coincides with the lower coseismic slip of previous earthquakes that occurred since 1939. The proposed method can be applied to other less‐studied strike‐slip faults to distinguish segments with shallow creep and strong coupling, and thus to better quantify the shallow strain budget and its associated seismic hazards.

Numerical simulation study on the spontaneous rupture process and its influencing factors of the 2001 MS 8.1 Kunlun earthquake in China

On November 14, 2001, the MS 8.1 Kunlun earthquake occurred in the west of Kunlun Mountain pass in the East Kunlun fault zone of the Qinghai-Tibet Plateau, China. The field investigation and the results of kinematic inversion indicate that it is a strike-slip earthquake with a complex rupture process. The study of the spontaneous rupture process of this earthquake and its influencing factors is of great significance for understanding the occurrence mechanism of continental large earthquakes. First, based on the results of geological survey and kinematic inversion, a three-dimensional non-planar fault model is established. Then, using the curved-grid finite difference method (CG-FDM), the dynamic rupture process of the 2001 MS 8.1 Kunlun earthquake is numerically reproduced by adjusting the parameters for numerical simulation experiments. On this basis, the influence of background stress field, fault geometry and friction coefficient on the surface slip distribution, seismic moment rate function and rupture velocity are discussed respectively. The results show that the important reasons for the super-long rupture scale and super-shear rupture velocity of this earthquake may be the clockwise rotation of the maximum horizontal principal compressive stress along the NNE ∼ NE direction and the low friction coefficient distribution on the whole fault surface. The difference between the dynamic and static friction coefficients along the fault strike may be non-uniform, and it may be larger on the main fault than that on the secondary fault. This will result in a higher level of stress drop on the main fault, leading to the occurrence of super-shear rupture and the slip distribution with the Kusai Lake as the macroseismic epicenter. The large local geometric variation along the strike of the main fault may also have some influence on the distribution of co-seismic slip and the location of the super-shear rupture. In addition, the numerical simulation results in this paper confirm the conclusion that the East Kunlun fault zone has low friction coefficient.

High seismic velocity structures control moderate to strong induced earthquake behaviors by shale gas development

Moderate to strong earthquakes have been induced worldwide by shale gas development, however, it is still unclear what factors control their behaviors. Here we use local seismic networks to reliably determine the source attributes of dozens of M > 3 earthquakes and obtain a high-resolution shear-wave velocity model using ambient noise tomography. These earthquakes are found to occur close to the target shale formations in depth and along high seismic velocity boundaries. The magnitudes and co-seismic slip distributions of the 2018 Xingwen {M}_{{{{{{rm{L}}}}}}}5.7 and 2019 Gongxian {M}_{{{{{{rm{L}}}}}}}5.3 earthquakes are further determined jointly by seismic waveforms and InSAR data, and the co-seismic slips of these two earthquakes correlate with high seismic velocity zones along the fault planes. Thus, the distribution of high velocity zones near the target shale formations, together with the stress state modulated by hydraulic fracturing controls induced earthquake behaviors and is critical for understanding the seismic potentials of hydraulic fracturing.

Fault Geometry, Slip Distribution, and Potential Triggering of the 2022 Mw 6.2 Deadly Afghanistan Earthquake Revealed from Geodetic and Weather Data

Abstract The 22 June 2022 Mw 6.2 Khōst, Afghanistan, earthquake struck killing more than 1700 people and devastating the region. For studying this earthquake, we computed the coseismic deformation fields of the earthquake using the Sentinel-1 Terrain Observation with Progressive Scans Interferometric Synthetic Aperture Radar (InSAR). The InSAR results show that the maximum coseismic displacement in the satellite line of sight direction reaches up to 39 cm. We determined the geometric parameters of the fault and coseismic slip distribution from these InSAR measurements. The best-fitted fault model shows that the rupture occurred on a right-lateral strike-slip fault with a strike of 203.7° and a dip of 68°. The most slip is concentrated at a shallow depth within the upper 10 km with the maximum slip of ∼3 m at 2.5 km depth. The maximum slip produced by this earthquake is significantly larger than the slip produced by several other similar earthquakes with similar magnitudes, implying that the focused shallow slip is likely the reason for the significant damage in the earthquake. The heavy rainfall was recorded during the earthquake period, which resulted in complicated fringes in coseismic interferograms close to the earthquake in time. Because a positive spatial and temporal correlation with the earthquake occurrence can be seen, the rainfall may have potential contributions to the earthquake, which deserves additional analysis in future. Combined with the potential effects of the 2015 Mw 7.5 Hindu Kush deep-seated earthquake, the seismicity in Afghanistan is the result of the ongoing subduction of the Indian plate beneath the Eurasian plate along their west boundary.

Geodetic constraints of the 2015Mw6.5 Alor, East Indonesia earthquake: a strike-slip faulting in the convergent boundary

SUMMARYOn 2015 November 4, an Mw 6.5 earthquake struck the east region of Alor Island, eastern Indonesia. Here, we use Sentinel-1 and ALOS-2 data to explore the coseismic surface displacement of this earthquake. Based on the ALOS-2 coseismic interferograms, the first fault model and coseismic slip distribution of the 2015 Alor earthquake are presented in this study. The preferred slip model links to a blind south-southeast striking, south-southwest dipping strike-slip fault with amount of normal slip component and a peak slip of 2.09 m at 2.34 km of depth. Considering the fact that the ascending and descending InSAR prediction of the distributed slip model does not fit the InSAR observations at the northwest and southeast tips with the simple planar fault model. We tried to construct a strike-variable fault model to further improve data fit. The results of the calculated Coulomb stress change imply that the regions with positive CFS changes mainly located at the northwest and southeast extremities of the rupture of the Alor earthquake, and in the lobes north and south of the rupture. The most striking discovery from the InSAR observations of the Alor earthquake is that most of the displacements occurred on a fault whose existence was unknown before the earthquake.

Source Model of the 2023 Turkey Earthquake Sequence Imaged by Sentinel-1 and GPS Measurements: Implications for Heterogeneous Fault Behavior along the East Anatolian Fault Zone

On 6 February 2023, a devastating doublet of earthquakes with magnitudes of Mw 7.8 and Mw 7.6 successively struck southeastern Turkey near the border of Syria. The earthquake sequence represents the strongest earthquakes in Turkey during the past 80 years and caused an extensive loss of life and property. In this study, we processed Sentinel-1 and GPS data to derive the complete surface displacement caused by the earthquake sequence. The surface displacements were adopted to invert for the fault geometry and coseismic slip distribution on the seismogenic faults of the earthquake sequence. The results indicate that the coseismic rupture of the Turkey earthquake sequence was dominated by left-lateral strike slips with a maximum slip of ~10 m on the East Anatolian Fault Zone (EAFZ) and the Sürgü fault (SF). Significant surface ruptures are recognized based on the geodetic inversion, which is consistent with the analysis of post-earthquake satellite images. The cumulative released moment of the two earthquakes reached 9.62 × 1020 Nm, which corresponds to an event of Mw 7.95. Additionally, the interseismic fault slip rates and locking depths along the central and western segments of the EAFZ were estimated using the high-resolution long-term velocity field. The results reveal significant lateral variations of fault slip rates and locking depths along the central and western segments of the EAFZ. Generally, the estimated fault locking zone showed good spatial consistency with the coseismic fault rupture of the Mw 7.8 shock on the EAFZ. The static coulomb failure stress (CFS) change due to the Mw 7.8 earthquakes suggests that the subsequent Mw 7.6 event was certainly promoted by the Mw 7.8 shock. The stress transfers from the fault EAFZ to the fault SF were realized by unclamping the interface of the fault SF, which significantly reduces the effective normal stress on the fault plane. Large CFS increases in the western Puturge segment of the EAFZ, which was not ruptured in the 2020 Mw 6.8 and the 2023 Mw 7.8 earthquakes, highlight the future earthquake risk in this fault segment.

Slip distribution of the February 6, 2023 Mw 7.8 and Mw 7.6, Kahramanmaraş, Turkey earthquake sequence in the East Anatolian Fault Zone

On February 6, 2023, two large earthquakes occurred near the Turkish town of Kahramanmaraş. The moment magnitude (Mw) 7.8 mainshock ruptured a 310 km-long segment of the left-lateral East Anatolian Fault, propagating through multiple releasing step-overs. The Mw 7.6 aftershock involved nearby left-lateral strike-slip faults of the East Anatolian Fault Zone, causing a 150 km-long rupture. We use remote-sensing observations to constrain the spatial distribution of coseismic slip for these two events and the February 20 Mw 6.4 aftershock near Antakya. Pixel tracking of optical and synthetic aperture radar data of the Sentinel-2 and Sentinel-1 satellites, respectively, provide near-field surface displacements. High-rate Global Navigation Satellite System data constrain each event separately. Coseismic slip extends from the surface to about 15 km depth with a shallow slip deficit. Most aftershocks cluster at major fault bends, surround the regions of high coseismic slip, or extend outward of the ruptured faults. For the mainshock, rupture propagation stopped southward at the diffuse termination of the East Anatolian fault and tapered off northward into the Pütürge segment, some 20 km south of the 2020 Mw 6.8 Elaziğ earthquake, highlighting a potential seismic gap. These events underscore the high seismic potential of immature fault systems.

Rupture Geometry and Slip Distribution of the Mw 7.2 Nippes Earthquake, Haiti, From Space Geodetic Data

AbstractOn 14 August 2021 the Mw 7.2 Nippes earthquake struck southern Haiti, rupturing a segment of the Enriquillo‐Plantain Garden Fault system (EPGF), a 300 km‐long strike‐slip fault system that accommodates half of the highly oblique convergence displacement between the Caribbean and the North American plates. We use coseismic surface displacements from Interferometric Synthetic Aperture Radar and Global Navigation Satellite System (GNSS) to estimate the geometry of the rupture through a systematic parametric exploration, determine its mechanism, and relate them to the regional tectonics derived from interseismic GNSS measurements. We show that the earthquake ruptured a north dipping fault (66 ± 4° dip) with a geodetically determined moment release that is 40% reverse and 60% strike‐slip. We cannot conclude whether this north‐dipping structure is the EPGF itself or a distinct fault running parallel to the EPGF. The rupture then evolved to the west on a vertical (86 ± 2° dip) fault parallel to the EPGF, the Ravine du Sud fault, with left‐lateral strike‐slip motion. The coseismic slip distribution of the 2010 Léogane and 2021 Nippes earthquakes, consistent with the transpressional interseismic strain rate field, show a segmentation of the Caribbean–North American plate boundary in southern Haiti and imply a revision in our understanding of the mode of earthquake rupture within the EPGF system.

Coseismic deformation and stress triggering of the 2021 MS6.4 Yangbi earthquake inverted from integrating GNSS and InSAR displacement fields

On 21 May 2021, the MS6.4 Yangbi earthquake occurred at the Tibetan Plateau’s southeastern margin along the Weixi-Qiaohou-Weishan fault. We obtained the fine coseismic deformation fields by integrating 37 near-field GNSS observations and Sentinel-1 InSAR observations of the ascending and descending orbits. Considering that the earthquake did not rupture to the surface, the existing earth fault data is difficult to effectively constrain the deep underground fault. We used aftershock precision location data to constrain fault strike, dip angle, and depth, effectively integrate GNSS and InSAR data, and obtained the fault model closest to the real fault. Based on the coseismic slip distribution model, we calculated the Coulomb stress generated from the earthquake and discussed the triggering relationship between the change of coseismic Coulomb stress and the aftershock distribution. The results show that the maximum rupture occurred at the southeast side of the epicenter and the overall unilateral rupture characteristic. We calculated the stress change when the Yangbi earthquake ruptured, significantly reducing the stress difference between the fault and the surrounding area, which has a significant unloading effect on the surrounding faults and reduces the seismic risk of the surrounding faults.

The 26 November 2019 Durrës earthquake, Albania: coseismic displacements and occurrence of slow slip events in the year following the earthquake

SUMMARYOuter Albanides experienced a seismic sequence starting on 21 September 2019, with an Mw 5.6 earthquake, considered a foreshock, and culminated with the main shock on 26 November 2019, followed by a paramount aftershock activity. We propose a model for the coseismic slip distribution using InSAR, permanent, and campaign GNSS measurements. We tested two hypotheses: an earthquake on a thrust plane with the direction N160° and along with a back thrust. By varying the depth and dip angle for the first hypothesis and only the dip angle for the second hypothesis, we concluded the optimal solution is a blind thrust at a 15-km depth dipping eastward 40°, a maximum slip of 1.4 m and an Mw 6.38. A GNSS time-series obtained after 2020 shows two slow slip events (SSEs): the first one is 200 d after the main shock up to 26 d, and the second one is 300 d after the main shock up to 28 d. We tested three hypotheses: SSE along the basement thrust where the main shock has been localized, SSE along the flat formed by the detachment layer of the cover, and SSE along these two faults. We concluded that SSE occurred along the detachment layer or along the two faults.

Coseismic Deformation Field and Fault Slip Distribution Inversion of the 2020 Jiashi Ms 6.4 Earthquake: Considering the Atmospheric Effect with Sentinel-1 Data Interferometry

Due to some limitations associated with the atmospheric residual phase in Sentinel-1 data interferometry during the Jiashi earthquake, the detailed spatial distribution of the line-of-sight (LOS) surface deformation field is still not fully understood. This study, therefore, proposes an inversion method of coseismic deformation field and fault slip distribution, taking atmospheric effect into account to address this issue. First, an improved inverse distance weighted (IDW) interpolation tropospheric decomposition model is utilised to accurately estimate the turbulence component in tropospheric delay. Using the joint constraints of the corrected deformation fields, the geometric parameters of the seismogenic fault and the distribution of coseismic slip are then inverted. The findings show that the coseismic deformation field (long axis strike was nearly east-west) was distributed along the Kalpingtag fault and the Ozgertaou fault, and the earthquake was found to occur in the low dip thrust nappe structural belt at the subduction interface of the block. Correspondingly, the slip model further revealed that the slips were concentrated at depths between 10 and 20 km, with a maximum slip of 0.34 m. Accordingly, the seismic magnitude of the earthquake was estimated to be Ms 6.06. Considering the geological structure in the earthquake region and the fault source parameters, we infer that the Kepingtag reverse fault is responsible for the earthquake, and the improved IDW interpolation tropospheric decomposition model can perform atmospheric correction more effectively, which is also beneficial for the source parameter inversion of the Jiashi earthquake.

Fault locking behavior of the Longmenshan Thrust Belt preceding the 2008 Wenchuan earthquake

Fault locking inversion can clarify the occurrence mechanisms of large earthquakes and the seismic hazards of future earthquakes. Due to limited observations in mainland China, fault locking inversions have typically been unstable. To investigate a more suitable fault locking behavior of the Longmenshan Thrust Belt (LTB) preceding the 2008 Wenchuan earthquake, GNSS velocity field data during 1999–2007 and the least-squares collocation method (LSCM) are used in this study. Despite the limited stations, after interpolating the GNSS velocity with the LSCM, uncertainties are reduced to some extent during the inversion. Based on the interpolation data, the inversion shows two locked areas that are consistent with the coseismic slip distribution of the Wenchuan earthquake. The third locked area is located on the southern LTB. A checkerboard scheme is proposed, and historic earthquakes larger than M 3.0 since 1950 are analyzed to recheck our inversion. According to the inverted result, the Wenchuan and Lushan earthquakes are two independent events. Additionally, we propose that future destructive earthquakes may occur on the southern segment of the LTB. The maximum magnitude of a future possible earthquake is Mw 8.2.

Geodetic modelling of the 2022Mw 6.6 Menyuan earthquake: insight into the strain-partitioned northern Qilian Shan fault system and implications for regional tectonics and seismic hazards

SUMMARYStrain 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 Mw 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.

The 8 January 2022, Menyuan Earthquake in Qinghai, China: A Representative Event in the Qilian-Haiyuan Fault Zone Observed Using Sentinel-1 SAR Images

On 8 January 2022, a Ms 6.9 earthquake occurred in Menyuan, Qinghai, China. This event provided important geodetic data before and after the earthquake, facilitating the investigation of the slip balance along the seismogenic faults to understand seismogenic behavior and assess seismic risk. In this study, we obtained the interseismic (2016–2021) and coseismic deformation fields of the 2022 earthquake using Sentinel-1 synthetic aperture radar (SAR) images and estimated the slip rate, fault locking, and coseismic slip of the seismogenic faults. The results indicated that the seismogenic fault of the 2022 Menyuan earthquake, i.e., the Tuolaishan–Lenglongling Fault, had shallow locked areas before the earthquake; its long-term slip rate could reach 6 ± 1.2 mm/yr. The earthquake ruptured a sinistral strike-slip fault with a high dip angle; the maximum slip magnitude reached 3.47 m, with a moment magnitude of 6.6. The area of coseismic slip > 1.5 m was equivalent to the range of the isoline, with a locking value of 0.6. The interseismic locking region can limit the approximate scope of the coseismic slip distribution. The 2022 Menyuan earthquake released energy that had accumulated over 482 years in the stepover region between the Lenglongling and Tuolaishan faults. The accumulated elastic strain power of the Tuolaishan Fault was equivalent to an Mw 6.79 earthquake. These circumstances in terms of the strain energy balance demonstrate that interseismic locking, as constrained from the geodetic data, and the elapsed time from the previous paleoseismic event are useful for earthquake location and energy predictions.

A New Method Applied for the Determination of Relative Weight Ratios Under the TensorFlow Platform When Estimating Coseismic Slip Distribution

AbstractWhen estimating a coseismic slip distribution with multiple observation types, the relative weight ratios (the contribution ratios of all kinds of actual observations in the joint inversion, which should be positive by default) and regularization parameter critically impact the inversion accuracy. In this paper, we propose a new method for determining the relative weights of multiple observations of jointing inversion slip distributions. This method regards the observations and the relative weight ratios as the training data sets and training parameters, respectively; in addition, the constructed function was optimized by the gradient descent method and the exponential decay learning rate method under the TensorFlow platform (GDED). We compared the GDED method with the Helmert variance component estimation method (HVCE) and Akaike's Bayesian information criterion method (ABIC) by designing eight simulation experiments. These simulations included changing the errors and resolutions of the observations, the grid size, the strike angle and the dip angle of the fault, the moment magnitudes of earthquakes, the number of observation types, and the complexity of subduction‐zone earthquakes. We analyzed and discussed the characteristics of the above three methods in determining the relative weight ratios during the joint inversion process. The results of these simulation experiments showed that the GDED method and the ABIC method perform much better than the HVCE method when negative variance occurs in the HVCE method, while the accuracies of the slip distributions produced by the three methods were similar when no negative variance occurred; the GDED method can be used to prevent the occurrence of negative variance compared to the HVCE method and is more time efficient than the ABIC method. Furthermore, we applied these methods to invert the source parameters of the 2015 Mw 8.3 Illapel earthquake (in Chile), and the inversion results showed that the slip of the Illapel earthquake was mainly distributed at depths of approximately 10–50 km. The maximum slip of the Illapel earthquake was approximately 7.9 m, distributed 50–70 km northwest of the epicenter at depths of approximately 10–20 km.

Earthquake Cycle Deformation Associated With the 2021 MW 7.4 Maduo (Eastern Tibet) Earthquake: An Intrablock Rupture Event on a Slow‐Slipping Fault From Sentinel‐1 InSAR and Teleseismic Data

AbstractIn the continents, the importance of earthquakes that occur away from major block‐bounding faults is still debated. The 21 May 2021 MW ∼ 7.4 Maduo earthquake occurred on a secondary fault away from previously‐identified major block boundaries. Here we use 7 years of Sentinel‐1 Interferometric Synthetic Aperture Radar (InSAR) time series (between October 2014 and November 2021) to determine the distribution of coseismic slip and early postseismic afterslip following the Maduo earthquake, and the preceding interseismic strain accumulation. We devised a 13‐segment 3‐D fault geometry constrained by the SAR range offsets and the distribution of relocated aftershocks and used a Bayesian method incorporating von Karman regularization to solve for coseismic slip and afterslip models. We also used teleseismic waveforms as a standalone inversion to show the rupture evolution in space and time during the earthquake, finding that it propagates bilaterally with three notable rupture episodes. Our preferred coseismic self‐similar slip model shows a moderate shallow slip deficit, with the majority of moment release occurring in the depth interval of 1–10 km. The coseismic slip deficit is taken up in part by afterslip at shallow (<4 km) depths that grows linearly with time during the first ∼6 months, and at >10 km depths where afterslip grows logarithmically with time. We suggest that this heterogeneity is likely controlled by spatial variations in fault friction related to lithology. We discuss the implications for seismic hazard away from major tectonic block boundaries in light of our observations of the earthquake cycle on this intrablock fault.

Structural Heterogeneity of the Alaska‐Aleutian Forearc: Implications for Interplate Coupling and Seismogenic Behaviors

AbstractWe determine high‐resolution 3‐D seismic images (Vp, Vs, and Vp/Vs) of the Alaska‐Aleutian forearc region, which are used to further estimate the serpentinization volume fraction of the subducting Pacific slab. Our results show that the hydration degree of the lower plate is high in the Shumagin Gap but relatively low in the Semidi and Kodiak segments, which roughly shows a trend of eastward decrease that is correlated with an increase of the interseismic slip deficit rate. The variation of slab serpentinization degree derived from our Vp model shows an uneven decreasing trend, well consistent with the along‐strike segmentation of interplate coupling and coseismic slip distributions of the 2020 Mw 7.8 Simeonof and 2021 Mw 8.2 Chignik megathrust earthquakes. The 2020 Mw 7.6 Sand Point earthquake took place within the strongly hydrated slab, whose occurrence was probably related to the variation of slip deficit rate that is probably affected by the changes in plate hydration. In addition, slow earthquakes occurred between the hydrated slab with pre‐subduction structures (e.g., bending faults, seamounts, and fracture zones) and the fluid‐saturated overriding crust and mantle wedge, where the intraslab seismicity and dehydration are intensive. These features suggest that lower‐plate hydration and dehydration may affect the interplate coupling and seismogenic behaviors in subduction zones.

Airborne LiDAR-Based Mapping of Surface Ruptures and Coseismic Slip of the 1955 Zheduotang Earthquake on the Xianshuihe Fault, East Tibet

ABSTRACTSurface ruptures and coseismic slip distributions of large earthquakes are the keys to understanding earthquake rupture processes, analyzing rupture history of associated faults, and assessing earthquake hazards. Detailed mapping of surface ruptures of large historical earthquakes is needed but is difficult in remote regions. The 1955 Ms 7.5 Zheduotang earthquake occurred in a prominent restraining bend of the central sinistral Xianshuihe fault and within a high-relief and densely vegetated mountain range. This study characterizes the 1955 earthquake surface rupture via 50 cm resolution airborne light detection and ranging (LiDAR) data, combined with field investigations. Our mapping of geomorphic features showed that relatively fresh mole tracks and fault scarps are still preserved beneath the dense vegetation cover. The 1955 earthquake ruptured at least 30 km in length, which consisted of four sections separated by stepovers or changes in the strike. Overall, the multistranded rupture was more complex than that of other historical earthquakes to the northwest along the Xianshuihe fault but consistent with its local structural setting. We collected 47 LiDAR-derived and 48 field-based left-lateral (with minor vertical) measurements. The clustering of the smallest offsets suggests that the average sinistral coseismic displacement of the 1955 earthquake was ∼2.1 m based on LiDAR-derived data and ∼1.5 m based on field measurements. This difference highlights the ambiguity and difficulties associated with surface rupture investigations of historical earthquakes decades after the event. The rupture length (∼30 km) and average sinistral displacement (1.5–2.1 m) suggest a moment magnitude of 6.6–7.4 based on empirical relationships of strike-slip earthquakes worldwide. The magnitude is smaller than the widely accepted Ms 7.5 in the catalog, suggesting that the previously reported magnitude was possibly overestimated. Our data have implications for the seismic hazard evaluation of the Sichuan–Tibet railway, which passes through the northern part of the ruptured fault.