Applying InSAR And GNSS Data to Obtain 3-D Surface Deformations Based on Iterated Almost Unbiased Estimation And Laplacian Smoothness Constraint

Global navigation satellite system (GNSS) and interferometric synthetic aperture radar (InSAR) data are integrated to extract the 3-D surface deformations, which are of great significance for studying geological hazards. In this study, two major problems are focused on integration. For one thing, we propose an iterated almost unbiased estimation (IAUE) method to estimate the variance components of GNSS and InSAR for the case where the estimation of variance components of multisource data by traditional variance component estimation methods may be negative and inaccurate. For another, considering that heterogeneous data errors may lead to unstable 3-D solutions, we propose adding the Laplacian smoothness constraint (LSC) to the function model, which can smooth the solutions by minimizing the second derivative of the displacements. These two methods are abbreviated as IAUE-LSC. In the simulation experiment, the performance of traditional Helmert variance component estimation is first compared with IAUE. IAUE can not only converge more quickly, but also avoid negative variances. Furthermore, we find that the excessively large relative error ratio between GNSS and InSAR is an essential factor leading to the instability of the 3-D solutions. The IAUE-LSC method is immune to this instability and can obtain more stable results. In addition, the 2018 Hawaii case demonstrates that IAUE achieves improvements of 2.58, 2.77, and 7.69 cm in the east, north, and up directions relative to the traditional weighted least-squares method, while the combined IAUE-LSC achieves improvements of 2.29, 0.32, and 1.68 cm compared to the IAUE alone.