Abstract:Summary In this study, we developed a new method that can significantly accelerate the forward modeling of gravity fields generated by large-scale tesseroids while keeping the computational accuracy as high as possible. The cost of the high efficiency is that the method only works under the assumptions that (1) all tesseroids in the same latitude band have the same horizontal dimension, (2) the computation points are located at the same surface level and aligned with the horizontal centers of tesseroids, and (3) each tesseroid has a constant or linearly varying density. The new method firstly integrates the kernel function of the Newton’s volume integral analytically in the radial direction to eliminate its dependence on the vertical dimension of the tesseroid, and then expands the integrated kernel function into a Taylor series up to a certain order. Because the Taylor series expansion term of the integrated kernel function is an odd or even function of the difference between the longitudes of the tesseroid and computation point, there exist shifting or swapping symmetry relations among the gravity field of tesseroids. Consequently, the shifting or swapping symmetry is extended to the tesseroids with unequal vertical dimensions. Numerical experiments using the spherical shell model are conducted to verify the effectiveness of the new method. The results show that the computational speed of the new method is about 30 times faster than that of the traditional method which employs the Gauss-Legendre quadrature rule and a two-dimensional adaptive subdivision approach, while keeping almost the same computational accuracy. When applying the new method to an ice shell with unequal thicknesses, the results reveal that the relative errors of calculating V, Vz, and Vzz are smaller than 10−8, 10−6, and 10−4, respectively if the Taylor series expansion is truncated at order 4, while the computational time consumed by the new method is about 7 times less than that of the traditional method. Finally, the influence of the truncation order on the computational accuracy and the strategies for dividing the latitude band into several parts to further improve the accuracy are discussed.

Fast large-scale gravity modeling using a high-order compact cascadic multigrid strategy

Parallel Hierarchical Multiscale Modelling of Hydro-Mechanical Problems for Saturated Granular Soils

A Coupled SPFEM/DEM Approach for Multiscale Modeling of Large‐deformation Geomechanical Problems

Computational technology and strategy for large-scale high-accuracy gravity modeling

An efficient cascadic multigrid solver for 3D magnetotelluric forward modeling problems using potentials

Fast 3‐D large‐scale gravity and magnetic modeling using unstructured grids and an adaptive multilevel fast multipole method

Efficient Large-Scale 3D Gravity Modelling Using a Fast Evaluate Kernel Matrix Combined with Compressed Matrix Techniques

Gravity field forward modeling using tesseroids accelerated by Taylor series expansion and symmetry relations

An Efficient Multigrid Solver Based on a Four-Color Cell-Block Gauss-Seidel Smoother for 3D Magnetotelluric Forward Modeling

Three-dimensional forward modelling of gravity field vector and its gradient tensor using the compact difference schemes

3D Controlled‐source electromagnetic modelling in anisotropic media using secondary potentials and a cascadic multigrid solver

An Adaptive Conjugate Gradient Least-Squares Regularization (ACGLSR) Method for 3D Gravity Density Inversion

Fast 3D gravity and magnetic modelling using midpoint quadrature and 2D FFT

Scalable robust solvers for unstructured FE geodynamic modeling applications: Solving the Stokes equation for models with large localized viscosity contrasts

Rapid forward of gravity and tensor gravity data

Forward modeling of gravity anomalies based on cell mergence and parallel computing

Efficient three-dimensional magnetotelluric forward modeling based on a geometric multigrid preconditioner

A Fast Forward and Inversion Strategy for Three-Dimensional Gravity Field

Efficiency of Optimized Approaches for Gravity Operator Modeling

3D Monte Carlo geometry inversion using gravity data

Generalized Gauss-FFT 3D forward gravity modeling for irregular topographic mass having any 3D variable density contrast