Efficiency of Optimized Approaches for Gravity Operator Modeling

Abstract

Numerical tesseroid and radial-type approaches are presented and compared in terms of their efficiency for deriving the regional geoid height, vertical gravity, and gradiometric anomalies from sea floor topography grids. The vertical gradient function is particularly suitable for representing shorter wavelengths of gravity, typically less than 10 km. These two modeling methods were applied to the Great Meteor guyot in the Atlantic Ocean using its bathymetry. To optimize the computation of high-resolution gravity anomalies, the Armadillo, GSL, and OpenMP libraries were used to provide an environment for fast vector implementation, numerical integration for tesseroid calculation, and parallelization for loop iterations, resulting in a computation speed increase. The tesseroid and radial methods remain equivalent up to a resolution of about 1 min, with the radial method being faster when dealing with a large number of model points for the geoid. Aside from optimization enabling high-resolution gravity simulations, these fast modeling data can be used as the main operators in gravimetric inversion or to reduce the terrain effects in gravity observations, revealing gravity and sedimentary layers.