TY - JOUR

T1 - Gravity Field Modeling Using Tesseroids with Variable Density in the Vertical Direction

AU - Lin, Miao

AU - Denker, Heiner

AU - Müller, Jürgen

N1 - Funding Information: We thank the editor and two anonymous reviewers for their constructive comments that helped to significantly improve the manuscript. This work was financially supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) − SFB 1128, projects C04 and C05. All figures were plotted by the Generic Mapping Tools (GMT; Wessel and Smith ).

PY - 2020/7

Y1 - 2020/7

N2 - We present an accurate method for the calculation of gravitational potential (GP), vector (GV), and gradient tensor (GGT) of a tesseroid, considering a density model in the form of a polynomial up to cubic order along the vertical direction. The method solves volume integral equations for the gravitational effects due to a tesseroid by the Gauss–Legendre quadrature rule. A two-dimensional adaptive subdivision technique, which automatically divides the tesseroids near the computation point into smaller elements, is applied to improve the computational accuracy. For those tesseroids having small vertical dimensions, an extension technique is additionally utilized to ensure acceptable accuracy, in particular for the evaluation of GV and GGT. Numerical experiments based on spherical shell models, for which analytical solutions exist, are implemented to test the accuracy of the method. The results demonstrate that the new method is capable of computing the gravitational effects of the tesseroids with various horizontal and vertical dimensions as well as density models, while the evaluation point can be on the surface of, near the surface of, outside the tesseroid, or even inside it (only suited for GP and GV). Thus, the method is attractive for many geodetic and geophysical applications on regional and global scales, including the computation of atmospheric effects for terrestrial and satellite usage. Finally, we apply this method for computing the topographic effects in the Himalaya region based on a given digital terrain model and the global atmospheric effects on the Earth’s surface by using three polynomial density models which are derived from the US Standard Atmosphere 1976.

AB - We present an accurate method for the calculation of gravitational potential (GP), vector (GV), and gradient tensor (GGT) of a tesseroid, considering a density model in the form of a polynomial up to cubic order along the vertical direction. The method solves volume integral equations for the gravitational effects due to a tesseroid by the Gauss–Legendre quadrature rule. A two-dimensional adaptive subdivision technique, which automatically divides the tesseroids near the computation point into smaller elements, is applied to improve the computational accuracy. For those tesseroids having small vertical dimensions, an extension technique is additionally utilized to ensure acceptable accuracy, in particular for the evaluation of GV and GGT. Numerical experiments based on spherical shell models, for which analytical solutions exist, are implemented to test the accuracy of the method. The results demonstrate that the new method is capable of computing the gravitational effects of the tesseroids with various horizontal and vertical dimensions as well as density models, while the evaluation point can be on the surface of, near the surface of, outside the tesseroid, or even inside it (only suited for GP and GV). Thus, the method is attractive for many geodetic and geophysical applications on regional and global scales, including the computation of atmospheric effects for terrestrial and satellite usage. Finally, we apply this method for computing the topographic effects in the Himalaya region based on a given digital terrain model and the global atmospheric effects on the Earth’s surface by using three polynomial density models which are derived from the US Standard Atmosphere 1976.

KW - Forward gravity modeling

KW - Gauss–Legendre quadrature

KW - Polynomial density model

KW - Tesseroids

KW - Topographic and atmospheric effects

UR - http://www.scopus.com/inward/record.url?scp=85083233149&partnerID=8YFLogxK

U2 - 10.1007/s10712-020-09585-6

DO - 10.1007/s10712-020-09585-6

M3 - Article

AN - SCOPUS:85083233149

VL - 41

SP - 723

EP - 765

JO - Surveys in geophysics

JF - Surveys in geophysics

SN - 0169-3298

IS - 4

ER -