A seafloor topography inversion iterative algorithm based on gravity forward modeling theory
The method of seafloor topography inversion from gravity data is an important way to obtain global seafloor topography.Aiming to solve the problem that traditional topography inversion algorithms need to indirectly introduce empirical parameters,this study,based on the principle of not adopting empirical formulas,established analytical observation equations between gravity anomalies and terrains by using the spatial domain method and the gravity orthogonal formulas,and solved the equations by the least squares method.In order to obtain the optimal solutions to the observation equations,Newtonian iterative relations were established and regularization parameters were introduced to enhance the convergence of the equations.Considering the influence of the edge effect on the regularization,a bilinear interpolation algorithm was used to encrypt the raw gravity anomaly data and the means of unconstrained grid expansion to weaken the error.Moreover,a 1'× 1'resolution seafloor topographic model was constructed in the Pacific Ocean(155°E～156°E,16°N～17°N),which was verified by the ship-surveyed bathymetry data,and the root mean squared error was reduced by about 12.7％compared with the traditional spatial-domain algorithm,gravity geology method,which verified the feasibility and accuracy of the proposed method.
万方数据
维普
