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.

Spherical radial basis functions model: approximating an integral functional of an isotropic Gaussian random field

Simulation of Gaussian random field in a ball

Regional Gravity Field Modeling Using Band-Limited SRBFs: A Case Study in Colorado

Gravity field modeling in mountainous areas based on band-limited SRBFs

Intrinsic Random Functions on the sphere

SRBF_Soft: a Python-based open-source software for regional gravity field modeling using spherical radial basis functions based on the data-adaptive network design methodology

Isotropic Variogram Matrix Functions on Spheres

Modeling random isotropic vector fields on the sphere: theory and application to the noise in GNSS station position time series

Adaptive Mesh-Free Approach for Gravity Inversion Using Modified Radial Basis Function

Adaptive selection strategy of shape parameters for LRBF for solving partial differential equations

Spherical Gravity Forwarding of Global Discrete Grid Cells by Isoparametric Transformation

Gaussian Variational Schemes on Bounded and Unbounded Domains

Stable and Accurate Least Squares Radial Basis Function Approximations on Bounded Domains

Non-Parametric Approximations for Anisotropy Estimation in Two-dimensional Differentiable Gaussian Random Fields

Non-Homogeneity Estimation and Universal Kriging on the Sphere

On a generalized Gaussian radial basis function: Analysis and applications

Estimating covariance functions of multivariate skew-Gaussian random fields on the sphere

Flexible Basis Representations for Modeling Large Non-Gaussian Spatial Data

A two-level approach to implicit surface modeling with compactly supported radial basis functions

Exploring substitution random functions composed of stationary multi-Gaussian processes

A stochastic model based on Gaussian random fields to characterize the morphology of granular objects