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

Research output: Contribution to journalArticleResearchpeer review

Authors

Research Organisations

External Research Organisations

  • Shandong University of Science and Technology
View graph of relations

Details

Original languageEnglish
Pages (from-to)723-765
Number of pages43
JournalSurveys in geophysics
Volume41
Issue number4
Early online date26 Mar 2020
Publication statusPublished - 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

Cite this

Gravity Field Modeling Using Tesseroids with Variable Density in the Vertical Direction. / Lin, Miao; Denker, Heiner; Müller, Jürgen.
In: Surveys in geophysics, Vol. 41, No. 4, 07.2020, p. 723-765.

Research output: Contribution to journalArticleResearchpeer review

Lin M, Denker H, Müller J. Gravity Field Modeling Using Tesseroids with Variable Density in the Vertical Direction. Surveys in geophysics. 2020 Jul;41(4):723-765. Epub 2020 Mar 26. doi: 10.1007/s10712-020-09585-6, 10.1007/s10712-021-09654-4
Download
@article{063b1c8605ef4abda35592a449ea7ea3,
title = "Gravity Field Modeling Using Tesseroids with Variable Density in the Vertical Direction",
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{\textquoteright}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",
author = "Miao Lin and Heiner Denker and J{\"u}rgen M{\"u}ller",
note = "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 ).",
year = "2020",
month = jul,
doi = "10.1007/s10712-020-09585-6",
language = "English",
volume = "41",
pages = "723--765",
journal = "Surveys in geophysics",
issn = "0169-3298",
publisher = "Springer Netherlands",
number = "4",

}

Download

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 -

By the same author(s)