A robust estimation algorithm for the increasing breakdown point based on quasi-accurate detection and its application to parameter estimation of the GNSS crustal deformation model

Wei Qu,Hailu Chen,Qin Zhang,Yuan Gao,Qingliang Wang,Ming Hao

DOI: https://doi.org/10.1007/s00190-021-01574-w

IF: 4.809

2021-10-25

Journal of Geodesy

Abstract:Global navigation satellite system (GNSS) velocity fields generally contain outliers due to environmental interference, local crustal activities, or strong seismic activities, which significantly affect the parameter estimation of the crustal deformation model. The robust M estimation can be used to improve the accuracy of the parameter estimation while retaining as much GNSS observation information as possible. The current robust estimation still uses the least-squares residuals instead of real errors as the initial values of the equivalent weight function and the breakdown point of parameters estimation is less than 50%. First, a new automatic selection strategy is proposed for quasi-accurate observations (observations that are reliable and do not have outliers but require confirmations) using quasi-accurate detection, and the outliers are roughly identified, almost independent of the breakdown point. Second, the variance of the unit weight of the estimate of real error of quasi-accurate observations is used as the initial variance factor and the estimate of real error of quasi-accurate observations is used as the initial value for the iterative calculation of the robust M estimation to accurately identify and detect large, medium, and small outliers. We evaluate the rigorousness of the least-squares estimation, conventional robust estimation, median method, and new algorithm from four perspectives using simulated GNSS velocity fields, i.e., model parameter estimation, variance factor, distribution of the observation weights, and accuracy of forward deformation fields. The new algorithm effectively eliminates outliers and estimates the model parameters. When outliers are below 50%, the robustness of the new algorithm is better than that of the other three. When the proportion of outliers is higher than 50% (60%, 70%, and 80%), the other three algorithms break down, while the new algorithm yields better parameter estimates with an increasing breakdown point.

geochemistry & geophysics,remote sensing