Spherical Gravity Forwarding of Global Discrete Grid Cells by Isoparametric Transformation

Shujin Cao, Peng Chen, Guangyin Lu (), Yihuai Deng, Dongxin Zhang and Xinyue Chen
Additional contact information
Shujin Cao: School of Earth Sciences and Spatial Information Engineering, Hunan University of Science and Technology, Xiangtan 411201, China
Peng Chen: School of Earth Sciences and Spatial Information Engineering, Hunan University of Science and Technology, Xiangtan 411201, China
Guangyin Lu: School of Geosciences and Info-Physics, Central South University, Changsha 410083, China
Yihuai Deng: School of Earth Sciences and Spatial Information Engineering, Hunan University of Science and Technology, Xiangtan 411201, China
Dongxin Zhang: School of Geosciences and Info-Physics, Central South University, Changsha 410083, China
Xinyue Chen: School of Earth Sciences and Spatial Information Engineering, Hunan University of Science and Technology, Xiangtan 411201, China

Mathematics, 2024, vol. 12, issue 6, 1-28

Abstract: For regional or even global geophysical problems, the curvature of the geophysical model cannot be approximated as a plane, and its curvature must be considered. Tesseroids can fit the curvature, but their shapes vary from almost rectangular at the equator to almost triangular at the poles, i.e., degradation phenomena. Unlike other spherical discrete grids (e.g., square, triangular, and rhombic grids) that can fit the curvature, the Discrete Global Grid System (DGGS) grid can not only fit the curvature but also effectively avoid degradation phenomena at the poles. In addition, since it has only edge-adjacent grids, DGGS grids have consistent adjacency and excellent angular resolution. Hence, DGGS grids are the best choice for discretizing the sphere into cells with an approximate shape and continuous scale. Compared with the tesseroid, which has no analytical solution but has a well-defined integral limit, the DGGS cell (prisms obtained from DGGS grids) has neither an analytical solution nor a fixed integral limit. Therefore, based on the isoparametric transformation, the non-regular DGGS cell in the system coordinate system is transformed into the regular hexagonal prism in the local coordinate system, and the DGGS-based forwarding algorithm of the gravitational field is realized in the spherical coordinate system. Different coordinate systems have differences in the integral kernels of gravity fields. In the current literature, the forward modeling research of polyhedrons (the DGGS cell, which is a polyhedral cell) is mostly concentrated in the Cartesian coordinate system. Therefore, the reliability of the DGGS-based forwarding algorithm is verified using the tetrahedron-based forwarding algorithm and the tesseroid-based forwarding algorithm with tiny tesseroids. From the numerical results, it can be concluded that if the distance from observations to sources is too small, the corresponding gravity field forwarding results may also have ambiguous values. Therefore, the minimum distance is not recommended for practical applications.

Keywords: gravity forward modelling; spherical hexagonal prisms; isoparametric transformation; discrete global grid system (search for similar items in EconPapers)
JEL-codes: C (search for similar items in EconPapers)
Date: 2024
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/2227-7390/12/6/885/pdf (application/pdf)
https://www.mdpi.com/2227-7390/12/6/885/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jmathe:v:12:y:2024:i:6:p:885-:d:1358598

Access Statistics for this article

Mathematics is currently edited by Ms. Emma He

More articles in Mathematics from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().