InSAR coseismic deformation fields caused by the Maduo M_W7.3 earthquake occurring on May 22, 2021 were generated using the C-band Sentinel-1A/B SAR images with D-InSAR technology. The spatial characteristics, magnitude of coseismic deformation and segmentation of the seismogenic fault were analyzed. The surface rupture trace was depicted clearly by InSAR observations. In addition, the coseismic slip distribution inversion was carried out constrained by both ascending and descending InSAR deformation fields and relocated aftershocks to understand the characteristics of deep fault slip and geometry of the seismogenic fault. The regional stress disturbance was analyzed based on coseismic Coulomb stress change. The results show that the Maduo M_W7.3 earthquake occurred on a secondary fault within the Bayan Har block which is almost parallel to the main fault trace of the Kunlun Fault. According to field investigation, geological data and InSAR surface rupture traces, the seismogenic fault is confirmed to be the Kunlunshankou-Jiangcuo Fault. The rupture length of seismogenic fault is estimated to be~210km. The NWW direction is followed by the overall displacement field, which indicates a left-lateral strike-slip movement of seismogenic fault. The maximum displacement is about 0.9m in LOS direction observed by both ascending and descending InSAR data. The inversion result denotes that the strike of the seismogenic fault is 276°and the dip angle is 80°. The maximum slip is about 6m and the average rake is 4°. The predicted moment magnitude is M_W7.45, which is overall consistent with the result of GCMT. An obvious slip-concentrated area is located at the depth of 0~10km. The coseismic Coulomb stress change with the East Kunlun Fault as the receiver fault shows that the Maduo earthquake produced obvious stress increase near the eastern segment of the East Kunlun Fault. Thus the seismic risk increases based on the high interseismic strain rate along this segment, which should receive more attention. In addition, the coseismic Coulomb stress change with the Maduo-Gande Fault as the receiving fault indicates that the Maduo earthquake produced an obvious stress drop near the western part of the Maduo-Gande Fault, which indicates that the Maduo earthquake released the Coulomb stress of the Maduo-Gande Fault, and its seismic risk may be greatly reduced. However, there is a stress loading effect in the intersection area of the Maduo-Gande Fault and the Kunlunshankou-Jiangcuo Fault. Considering that aftershocks of Maduo earthquake will release excess energy, the greater earthquake risk may be reduced.

In the global scale, ten destructive earthquakes with magnitude larger than 7 happen on average each year. Yet the number of small earthquakes with limited or even no damage but recordable by seismographs(magnitude between 2.5 and 4.5)is over one million per year. In between, there are hundreds to thousands of earthquakes with moderate to strong magnitude(magnitude between 5.5 and 6.5)with notable destructiveness. The massive moderate to strong earthquakes are often less noticed or even overlooked, with only very few exceptions which caused human casualties and/or structure damages due to the very shallow focal depths. For medium earthquakes, the traditional seismology means can obtain the source mechanism solution of earthquake, but because of the inherent fuzziness of the source mechanism, it cannot distinguish the fault plane from the auxiliary nodal planes, because earthquakes of this magnitude usually do not produce surface rupture, and the result error is large, so it is not suitable for the study of medium and small earthquakes. It is of fundamental significance to further study the source fault of the moderate earthquakes, and more independent methods other than traditional seismology, such as satellite geodesy are needed. As one of the most applied satellite geodesy technique, interferometry of SAR(InSAR)satellite images are commonly used to obtain coseismic deformation related to earthquakes. InSAR has very high spatial sampling, though the temporal sampling is very low, which is several days to over a month depending on the satellite revisit span. The precision of coseismic deformation by InSAR can reach 2~3cm, which is good enough to obtain the surface deformation caused by a moderate earthquake. It is noted that InSAR coseismic measurements can detect 1-dimensional(1D)deformation along Line-of-Sight(LOS)direction. With multiple observing modes including descending and ascending, the InSAR deformation data is very useful for identifying surface ruptures, and for source fault plane discrimination. As a new geodetic observation technology, InSAR uses the elastic dislocation model to obtain source parameters, and the inversion results of fault parameters and slip distribution are more reliable. On September 24th, 2019, an M_W6.0 earthquake hit New Mirpur, Pakistan. The nearest known fault to the epicenter is the Main Frontal Thrust on its south side. We used the Sentinel-1A SAR imagery(TOPS-model)to reconstruct the InSAR coseismic deformation fields generated by the 2019 M_W6.0 Pakistan earthquake along the ascending and descending tracks. The ascending and descending deformation fields indicate that coseismic deformation is asymmetric by a trend of NW-SE in the south secondary fault of the Himalayan frontal thrust fault, with a maximum LOS displacement of~0.1m. The structures of ascending and descending deformation are highly consistent with each other, but the LOS displacement of southern side is obviously larger than the northern side. The continuous change of interference fringes between uplift and subsidence areas shows that there is no coherent phenomenon caused by excessively large deformation gradient or surface rupture, which indicates that the seismic fault rupture did not reach to the ground surface. Two initial fault models constrained by InSAR deformation, with a southwest-dipping and northeast-dipping fault, were utilized in the inversion. We finally determined the northeast-dipping fault as the seismogenic fault by joint inversion of ascending and descending observations, combined with tectonic setting. Our fault model suggests that an obvious slip concentrated area is located in the depth of 2~4km, with a peak slip of~0.64m and a mean rake angle of~125°. The north-dipping thrust motion with a small amount of strike-slip component dominated the faulting. The earthquake occurred in the low-dipping subduction zone between the Indian and Eurasian plates. The dip angle of the fault plane is relatively low. When the fault is ruptured, the upper wall thrust southwards and the north wall subducted northwards. Due to the compressional nappe structure, the front end of the upper wall was uplifted and the back end was stretched to become the subsidence area. Seismogenic fault is the south secondary fault of the Himalayan frontal thrust fault inferred from our coseismic fault model and rupture kinematic features. Active faults on the land have caused many large destructive earthquakes, resulting in surface faults and promoting the development of tectonic landforms. The detailed observation of coseismic surface rupture not only provides basic information for understanding the earthquake itself and estimating the earthquake recurrence period, but also helps to interpret the tectonic and geomorphic features in other areas. Since the M_W6.0 earthquake in Pakistan in 2019, no studies have been reported yet on this earthquake using InSAR technology, so the study of this earthquake provides a rare opportunity to assess the seismic risk of active thrust faults and to study the seismicity of northern Pakistan.

The continuous collision between Tian Shan and Tarim Basin causes not only the uplift of mountains, but also the earthquakes across the entire Tian Shan, particularly in the transient zone from mountains to the adjacent basins, where the critical infrastructures and residents are seriously under threat from these earthquake hazards. On 19^th January, 2020, an earthquake occurred in the Kalpintage fold thrust belt in the southwest Tian Shan foreland. We call this event the 2020 M_W6.0 Kalpintage earthquake, which is the first moderate earthquake captured by modern geodetic measurement techniques. This event therefore provides a rare opportunity to look into the local tectonics and seismic risk in southwest Tian Shan. In this study, we obtained the coseismic deformation of 2020 M_W6.0 Kalpintage earthquake from Sentinel-1A SAR and strong motion data, and then inverted its kinematic slip model. We derived the InSAR interferograms from both ascending and descending tracks. Both of them present similar deformation patterns, two deformation peaks over the Kalpintage anticline. That means: 1)The surface deformation is dominated by vertical displacement, and 2)the coseismic rupture plane is highly suspected to be the shallowly dipping decollement at the base of the sediment cover. We got the 3-D displacements of 6 strong motion stations by double integrating the strong motion acceleration signals. The result shows tiny displacement on the strong motion stations, except for the Xikeer station, which locates at the front of the Kalpintage anticline, where the InSAR interferograms are seriously incoherent. Two slip models can equally fit to the ascending and descending InSAR interferograms: One is a strike slip model with strike of N-S, the other is a thrust model with strike of E-W. This ambiguity in the slip models for the M_W6.0 Kalpintage earthquake is caused by 1)the extremely small dip angles of the causative fault, 2)the inherent shortcomings of the InSAR measurements i.e. the 1-D measurements along the line of sight, the polar orbiting direction of the SAR satellite, and 3)the serious atmospheric delay due to contrasting topography in southwest Tian Shan. We did not distinguish the two ambiguous models with InSAR data due to the weak constraints of InSAR for this event. However, the two quite different slip models show the same spatial dimension and position beneath the Kalpintage anticline, also the same seismic slip vector moving toward the Tarim Basin. We then presumed the two slip models refer to the same fault plane, the weak decollement at the base of the sediment cover, and its rupture released the compressive strain in this fold and thrust belt in the southwest Tian Shan front. The confusing problem is neither the strike slip model nor the thrust model can explain the displacement derived from strong motion. The simple error estimates show small uncertainty in the strong-motion-derived displacement, but we cannot really know the real errors without the comparison to the collocated continuous GNSS observation. Because of the discrepancy between the strong motion displacement and InSAR-derived slip model, we speculate the inelastic deformation occurred in front of the Kalpintage anticline where thick weak sediments exist. We think this earthquake ruptured the decollements in the lower sediments bounded by the adjacent anticlines, which are uplifted in this event. The M_W6.0 Kalpintage earthquake balanced the stress accommodated during the convergence of the Tian Shan and Tarim Basin. We managed to explain all of the ruptures in the southwest Tian Shan by combining the regional tectonic, geophysical data and the available earthquake catalogues with good quality and then estimated the earthquake hazards. The earthquakes, including 1902 M_W7.7 karshigar, 1996 M_W6.3 Jiashi, 1997—2003 Jiashi sequence and 2020 M_W6.0 Kalpintage earthquake, can be explained in one frame, the underthrusting of the Tarim Basin toward the southwest Tian Shan. Our calculation suggests that a M_W7.0+ event could be generated around Kalpintage anticline belt if without barriers on the decollements.