The 2021 Mw7.3 Maduo (Qinghai, China) earthquake ruptured ∼160 km along a poorly known secondary fault inside the Bayanhar block on the northern Tibetan plateau, which is generally parallel to the Kunlun fault. Here we integrate the interseismic (2015–2020) and coseismic geodetic observations to quantify the interseismic strain rate, fault geometry and coseismic slip distribution. Our results reveal that the seismogenic fault is featured by the low (<20 nanostrain/yr) and nearly undetectable interseismic strain rate. Three‐dimensional displacement fields and coseismic strain maps demonstrate the spatial variations of rupture kinematics due to the change of fault geometry. Our study reveals the steeply north‐dipping fault geometry of the seismogenic fault. The majority of coseismic slip occurred between 0 and 15 km with slight shallow slip deficit, which not penetrates through the inferred elastic upper crust (∼20–25 km). Our study highlights the tectonic and geometric control on fault kinematics of the Maduo earthquake.

The 21 May 2021 Maduo earthquake was the largest event to occur on a secondary fault in the interior of the active Bayanhar block on the north-central Tibetan plateau in the last twenty years. A detailed kinematic study of the Maduo earthquake helps us to better understand the seismogenic environments of the secondary faults within the block, and its relationship with the block-bounding faults. In this study, firstly, SAR images are used to obtain the coseismic deformation fields. Secondly, we use a strain model-based method and steepest descent method (SDM) to resolve the three-dimensional displacement components and to invert the coseismic slip distribution constrained by coseismic displacement fields, respectively. The three-dimensional displacement fields reveal a dominant left-lateral strike-slip motion, local horizontal displacement variations and widely distributed near-fault subsidence/uplift deformation. We prefer a five-segment fault slip model, with well constrained fault geometry featuring different dip angles and striking, constrained by InSAR observations. The peak coseismic slip is estimated to be ~5 m near longitude 98.9°E at a depth of ~4–7 km. Overall, the distribution of the coseismic slip on the fault is highly correlated to the measured surface displacement offsets along the entire rupture. We observe the moderate shallow slip deficit and limited afterslip deformation following the Maduo earthquake, it may indicate the effects of off-fault deformation during the earthquake and stable interseismic creep on the fault. The occurrence of the Maduo earthquake on a subsidiary fault updates the importance and the traditional estimate of the seismic hazards for the Kunlun fault.

We model the time‐dependent postseismic displacements following the 2001 Mw7.8 Kokoxili earthquake, including both GPS (2001–2002 for near‐field and 2001–2010 for far‐field) and descending‐track InSAR line‐of‐sight time series (2003–2010) to study three postseismic deformation processes. The predicted deformation patterns and magnitude from poroelastic rebound are inconsistent with the geodetic observations. Far‐field postseismic deformation (>200 km from rupture) is primarily induced by upper mantle viscoelastic relaxation beneath Tibet and the Qaidam Basin and places a lower bound on transient and steady‐state viscosities on the order of 10¹⁹–10²⁰ Pas. Shallow stress‐driven afterslip (20 km) has an amplitude of ∼1 m during the first 3 years. Our results indicate that the consideration of deep afterslip increases our estimates of transient viscosity of the lower crust beneath Tibet and the Qaidam basin for this earthquake by as much as a factor of three. Our combined model incorporating viscoelastic relaxation and afterslip suggests that the effective transient and steady‐state viscosities in the Tibetan lower crust are 5 × 10¹⁸ Pas and 4 × 10¹⁹ Pas, respectively (transient viscosity = 2 × 10¹⁸ Pas without afterslip considered), while the effective transient and steady‐state viscosities below the Qaidam Basin area are 1 × 10¹⁹ Pas and 6 × 10¹⁹ Pas (transient viscosity = 4 × 10¹⁸ Pas without afterslip considered).

Interferometric synthetic aperture radar (InSAR) data from 6 Envisat ASAR descending tracks; spanning the 2003–2010 period; was used to measure interseismic strain accumulation across the Northeastern Tibetan Plateau. Mean line-of-sight (LOS) ratemaps are computed by stacking atmospheric-corrected and orbital-corrected interferograms. The ratemaps from one track with different atmospheric-corrected results or two parallel; partially overlapping tracks; show a consistent pattern of left-lateral motion across the fault; which demonstrates the MERIS and ECMWF atmospheric correction works satisfactorily for small stain measurement of this region; even with a limited number of interferograms. By combining the measurements of InSAR and GPS; a fine crustal deformation velocity and strain rate field was estimated on discrete points with irregular density depending on the fault location; which revealed that the present-day slip rate on the Haiyuan fault system varies little from west to east. A change (2–3 mm/year) in line-of-sight (LOS) deformation rate across the fault is observed from the Jinqianghe segment to its eastern end. Inversion from the cross-fault InSAR profiles gave a shallow locking depth of 3–6 km on the main rupture of the 1920 earthquake. We therefore infer that the middle-lower part of the seismogenic layer on the 1920 rupture is not yet fully locked since the 1920 large earthquake. Benefit from high spatial resolution InSAR data; a low strain accumulation zone with high strain rates on its two ends was detected; which corresponds to the creeping segment; i.e., the Laohushan fault segment. Contrary to the previous knowledge of squeezing structure; an abnormal tension zone is disclosed from the direction map of principal stress; which is consistent with the recent geological study. The distribution of principal stress also showed that the expanding frontier of the northeastern plateau has crossed the Liupan Shan fault zone; even arrived at the northeast area of the Xiaoguan Shan. This result agrees with the deep seismic reflection profile.

Taking the 2017 M_w6.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.