Abstract：For the first time, we construct a three-dimensional finite element model for the 2021 Madoi MW7.3 earthquake to analyze the features of coseismic slip distribution. First, we process the high-precision GNSS coseismic deformation observations. Then, we construct a three-dimensional finite element model, which is validated by the elastic half-space Okada model. Finally, we invert the coseismic slip distribution model of the Madoi earthquake using the least squares method, by constraining 40 near-field and far-field GNSS coseismic displacement observations. The results show that the Madoi earthquake reptures with a maximum value of about 3.4 m are mainly distributed near the area of Yematan and Huanghe country. The results are consistent with the field investigation and can well explain the GNSS coseismic deformation observations.