Fusion of Spatially Heterogeneous GNSS and InSAR Deformation Data Using a Multiresolution Segmentation Algorithm and Its Application in the Inversion of Slip Distribution

The fusion of global navigation satellite system (GNSS) and interferometric synthetic aperture radar (InSAR) deformation data can leverage the advantages of GNSS high temporal resolution and InSAR high spatial resolution, and obtain more abundant deformation data for constraints on geophysical structural and mechanical parameters. Existing studies seldom consider the spatial heterogeneity of largescale deformation data, which easily leads to obvious spatial aggregation of errors in the results of fusion. Here, we propose a novel multiresolution segmentation fusion (MRSF) method that uses a multiresolution segmentation algorithm to automatically classify the spatial heterogeneity of InSAR deformation data with similar deformation characteristics. We applied the MRSF method to the fusion of GNSS and InSAR deformation data covering the central valley aquifer system (CVAS) in southern California to verify its precision and robustness. Results show that the MRSF method can accurately reflect spatiotemporal evolution characteristics of displacement data and reliably estimate deformation for the times and locations of missing data. We then tested this method for geophysical parameter estimation by constructing three different sets of data, including dense GNSS sites, sparse GNSS sites, and sparse GNSS sites fused with InSAR data using MRSF, to invert the slip distribution of the Cascadia subduction zone. Results show that the inverted slip of the fused InSAR and GNSS data is comparable to that of the dense GNSS sites. Therefore, the MRSF method can obtain deformation results with high precision and high spatiotemporal resolution and effectively compensate for the lack of data caused by sparse GNSS sites during the geophysical inversion process.