Details
Original language | English |
---|---|
Pages (from-to) | 723-765 |
Number of pages | 43 |
Journal | Surveys in geophysics |
Volume | 41 |
Issue number | 4 |
Early online date | 26 Mar 2020 |
Publication status | Published - Jul 2020 |
Abstract
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.
Keywords
- Forward gravity modeling, Gauss–Legendre quadrature, Polynomial density model, Tesseroids, Topographic and atmospheric effects
ASJC Scopus subject areas
- Earth and Planetary Sciences(all)
- Geophysics
- Earth and Planetary Sciences(all)
- Geochemistry and Petrology
Cite this
- Standard
- Harvard
- Apa
- Vancouver
- BibTeX
- RIS
In: Surveys in geophysics, Vol. 41, No. 4, 07.2020, p. 723-765.
Research output: Contribution to journal › Article › Research › peer review
}
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 -