Vertical deformation analysis based on combined adjustment for GNSS and leveling data

A method is proposed to fuse the velocity data of the global navigation satellite system (GNSS) and leveling height via combined adjustment with constraints. First, stable GNSS–leveling points are uniformly selected, and the constraints of the geodetic height change velocity and normal height change velocity are given. Then, the GNSS vertical velocities and leveling height difference are used as observations of combined adjustment, and robust least-squares estimation are used to estimate the velocities of the unknown points. Finally, a vertical movement model is established with the GNSS vertical velocities and leveling vertical velocities obtained via combined adjustment. Data from the second-order leveling network and GNSS control points in Shandong Province are taken as test data, and eight calculation schemes are used for discussion. One of the schemes, the bifactor robust combined adjustment method based on variance component estimation with two kinds of vertical velocity constraints achieves the optimal results. The method applied in the scheme can be recommended for data fusion of GNSS and leveling, further improving the reliability of vertical crustal movement in Shandong Province.