Complex Slip Distribution of the 2021 Mw 7.4 Maduo, China, Earthquake: An Event Occurring on the Slowly Slipping Fault

The 21 May 2021 Maduo earthquake occurred on the Kunlun Mountain Pass–Jiangcuo fault (KMPJF), a seismogenic fault with no documented large earthquakes. To probe its kinematics, we first estimate the slip rates of the KMPJF and Tuosuo Lake segment (TLS, ∼75 km north of the KMPJF) of the East Kunlun fault (EKLF) based on the secular Global Positioning System (GPS) data using the Markov chain Monte Carlo method. Our model reveals that the slip rates of the KMPJF and TLS are 1.7 ± 0.8 and 7.1 ± 0.3 mm/yr, respectively. Then, we invert high-resolution GPS and Interferometric Synthetic Aperture Radar observations to decipher the fault geometry and detailed coseismic slip distribution associated with the Maduo earthquake. The geometry of the KMPFJ significantly varies along strike, composed of five fault subsegments. The most slip is accommodated by two steeply dipping fault segments, with the patch of large sinistral slip concentrated in the shallow depth on a simple straight structure. The released seismic moment is ∼1.5×1020 N·m, equivalent to an Mw 7.39 event, with a peak slip of ∼9.3 m. Combining the average coseismic slip and slip rate of the main fault, an earthquake recurrence period of ∼1250−400+1120 yr is estimated. The Maduo earthquake reminds us to reevaluate the potential of seismic gaps where slip rates are low. Based on our calculated Coulomb failure stress, the Maduo earthquake imposes positive stress on the Maqin–Maqu segment of the EKLF, a long-recognized seismic gap, implying that it may accelerate the occurrence of the next major event in this region.

Surface Slip Distribution and Earthquake Rupture Model of the Fuyun Fault, China, Based on High-Resolution Topographic Data

The characteristics of earthquake surface ruptures, such as geometry, slip distribution, and coseismic deformation, contain important information about the earthquake rupture process, and so investigations and analyses of earthquake surface rupture have played a crucial role in modern earthquake hazard studies, especially with the increasing availability of high-resolution topographic and imagery data for tectono-geomorphic interpretation. In this study, we use Structure from Motion (SfM) photogrammetry to build a 1 m resolution digital elevation model (DEM) of the fault and combine this with filed observations to map the surface ruptures of the 1931 M8.0 Fuyun earthquake, China. These high-resolution topographic data enable to identify and measure the displaced gullies, and so the rupture locations and along-strike slip distribution are obtained in detail. Four paleoearthquake events are identified through the offset cluster characteristics. The coseismic offset of the 1931 Fuyun earthquake is extracted from the offset distribution, which shows four continuous undulations along the fault strike, corresponding to the four segments of surface rupture. Moreover, a high offset gradient is observed in the step area connected by the rupture segment. These findings, combined with the width and bending angle of the step area at the joint of the rupture segment, indicate that the 1931 Fuyun earthquake was a cascade rupture formed by four rupture segments. This study expands the available offset measurement data of Fuyun fault and confirms the applicability of high-resolution topographic data to active tectonic research.

Re-examination of Slip Distribution of the 2004 Sumatra–Andaman Earthquake (Mw 9.2) by the Inversion of Tsunami Data Using Green’s Functions Corrected for Compressible Seawater Over the Elastic Earth

We re-examined the slip distribution on faults of the 2004 Sumatra–Andaman (M 9.1 according to USGS) earthquake by the inversion of tsunami data with phase-corrected Green’s functions applied to linear long waves. The correction accounts for the effects of compressibility of seawater, elasticity of solid earth, and gravitational potential variation associated with the motion of mass to reproduce the delayed arrivals and the reversed phase of the first tsunami waves. We used sea surface height (SSH) data from satellite altimetry (SA) measurements along five tracks, and the tsunami waveforms recorded at tide gauges (TGs) and ocean bottom pressure gauges (OBPGs) in and around the Indian Ocean. The inversion results for both data sets for different rupture velocities (Vr) show that the reproducibility of the spatiotemporal SSHs and tsunami waveforms is improved by the phase corrections, although the effects are not so significant within the Indian Ocean. The best slip distribution model from joint inversion of SA, TG and OBPG data with Vr of 1.3 km/s shows the largest slips of 16–25 m off Sumatra Island, large slips of 2–11 m off the Nicobar Islands, and moderate slips of 2–6 m in the Andaman Islands. The inversion results reproduce the far-field tsunami waveforms well at distant stations even more than 13,000–25,000 km from the epicenter. The total source length is about 1400 km and the seismic moment is Mw 9.2, longer and larger than that of our previous estimates based on TG records.

Fault Geometry and Slip Distribution of the 2021 Mw 7.4 Maduo, China, Earthquake Inferred from InSAR Measurements and Relocated Aftershocks

A nearly 70 yr hiatus of major seismic activity in the central eastern Bayan Har block (BKB) ended on 22 May 2021, when a multislip-peak sinistral strike-slip earthquake struck western Maduo County, Qinghai. This earthquake, which ruptured the nearly 170 km long Kunlun Pass–Jiangcuo fault, is a rather unique event and offers a rare opportunity to probe the mechanical properties of the intraplate lithosphere of the central eastern BKB. Here, we inferred the fault geometry associated with the Maduo earthquake using Interferometric Synthetic Aperture Radar (InSAR), and relocated aftershocks and inverted the slip distribution through InSAR radar phases and range offsets. Our analysis revealed that the geometry of the fault varies along the strike: the southeastern end of the fault dips steeply to the northeast, whereas the northwestern end dips southwestward. Using the combined datasets to constrain a coseismic slip, we found that the 2021 Maduo event was dominated by sinistral strike-slip movement, with a slight normal-slip component at a shallow depth, rupturing the steep-dipping fault for nearly 170 km in length. Five asperities were detected along the fault strike in the shallow crust (0–12 km) with a peak slip of ∼4.2 m corresponding mostly to simple structures, namely, continuous and straight rupture segments, suggesting that the rupture propagated across geometrical barriers in a multiasperity way. Based on an analysis of the strain field and the focal mechanisms of both the 2021 Maduo earthquake and historical earthquakes that have occurred in the BKB, we propose that the fault zones within the BKB can also generate large earthquakes and have the ability to accommodate the ongoing eastward and northeastward penetration of the Indian plate into the Eurasian plate.

Slip distribution effect in spatial coulomb stress analysis (Case study: Palu earthquake on September 28, 2018)

On September 28, 2018, the Palu-Koro fault released the accumulated stress that caused the earthquake. An earthquake with magnitude 7.5 caused large and massive damage around Palu. There were many aftershocks along the Palu-Koro fault. This research aims to calculate a model of spatial Coulomb stress based on this event to find a correlation between mainshock and the aftershocks. The slip distribution was used as an input of the spatial stress Coulomb modeling to increase the accuracy. We use the Teleseismic Body-Wave Inversion method to calculate slip distribution along the fault plane. As a result, this earthquake was generated by the Palu-Koro fault movement with Mw 7.48, strike 350°, dip angle 67°, and rake -9°. There are three asperity zones along the fault plane located in the north and southern parts of the fault plane. The location of the most energy discharge is in the south asperity zone of the fault plane model with a maximum slip value of 1.65 meters. The spatial Coulomb stress change of this event shows that aftershocks concentration are in areas experiencing increased stress after the earthquake.