A novel method for inverting coseismic 3D surface deformation using InSAR considering the weight influence of the spatial distribution of GNSS points.

The combination of InSAR and GNSS data can provide more accurate measurements of coseismic 3D surface deformations, which are crucial in determining earthquake source parameters. Currently, the methods for inverting coseismic 3D surface deformations by integrating GNSS and InSAR data using the strain model (SM) (e.g., Simultaneous and integrated strain tensor estimation from geodetic measurements to obtain 3D displacement maps (SISTEM), SM-VCE) simply select GNSS points based on a certain number or distance, without considering the weight influence of the spatial distribution of G NSS points on the InSAR and GNSS data fusion results. To address this issue, this paper proposes a new method that uses Voronoi diagrams analyze the contribution of GNSS points to the deformation at target points. In this paper, GNSS points with high contribution can be used to construct observation value equations and optimize the internal weight of GNSS data. The proposed approach characterizes the contribution of GNSS data to target points by analyzing the changes in the Delaunay triangulation area formed before and after the target points are added to the GNSS data. Simulation experiments reveal that the proposed method outperforms conventional methods in accuracy and completeness. The proposed method is then applied to map the coseismic 3D surface deformations of the 2020 Mw6.5 Monte Cristo Range earthquake and find that the maximum deformation in the east–west, north–south, and vertical directions is about 13 cm, 11 cm, and 24 cm, respectively. [ABSTRACT FROM AUTHOR]