In order to explore the constraint degree of GNSS and InSAR data on the Maduo earthquake under the layered earth model structure, and to understand the cause of the Maduo earthquake and the danger of the surrounding area, this paper uses D-InSAR technology to obtain the InSAR co-seismic deformation field of the Maduo earthquake on May 22,2021. Based on GNSS and InSAR data, the co-seismic slip model and fault plane stress distribution of the Maduo earthquake are jointly inverted. This paper calculates the Coulomb stress changes caused by 14 M ≥ 7 strong earthquakes co-seismic, post-seismic viscoelastic relaxation and inter-seismic tectonic stress loading of 19 fault segments in the Bayan Har block research area (96°E-106°E, 29°N-36°N) since 1900. The results show that: (1) The maximum Line Of Sight (LOS) deformation obtained based on the ascending and descending InSAR data was about 0.9 m. The maximum slip amount of GNSS + InSAR joint inversion is 4.87 m, which is located in the Dongcao Along Lake section (∼98.6°E). Combined with the stress drop distribution of the fault surface, it is found that the Maduo earthquake is dominated by strike-slip. (2) The historical strong earthquakes delayed the Maduo earthquake by nearly 103 years. The 1937 Huashixia earthquake, the 1947 Dari earthquake, the 1963 Dulan earthquake, and the 1973 Luhuo earthquake had the greatest stress unloading effect, which was equivalent to delaying the Maduo earthquake by 53 years, 21 years, 11 years, and 17 years, respectively. The cumulative Coulomb stress at the source location of the Maduo earthquake reached 0.0112 MPa. Most area of the fault sections such as the East Kunlun Fault, the Maduo-Gande Fault, the western and eastern sections of the Dari Fault, the western and eastern sections of the Qingshuihe Fault, the Ganzi-Yushu Fault, and the eastern section of the Kunlunshankou-Jiangcuo Fault are in the Coulomb stress loading area, which should be paid attention to.