Efficient 3-D Large-Scale Forward Modeling and Inversion of Gravitational Fields in Spherical Coordinates With Application to Lunar Mascons



Zhao, Guangdong, Chen, Bo, Uieda, Leonardo ORCID: 0000-0001-6123-9515, Liu, Jianxin, Kaban, Mikhail K, Chen, Longwei and Guo, Rongwen
(2019) Efficient 3-D Large-Scale Forward Modeling and Inversion of Gravitational Fields in Spherical Coordinates With Application to Lunar Mascons. JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 124 (4). pp. 4157-4173.

This is the latest version of this item.

Access the full-text of this item by clicking on the Open Access link.
[img] Text
zhao2019.pdf - Author Accepted Manuscript
Available under License : See the attached licence file.

Download (1MB)

Abstract

<jats:title>Abstract</jats:title><jats:p>An efficient forward modeling algorithm for calculation of gravitational fields in spherical coordinates is developed for 3‐D large‐scale gravity inversion problems. The 3‐D Gauss‐Legendre quadrature (GLQ) is used to calculate the gravitational fields of mass distributions discretized into tesseroids. Equivalence relations in the kernel matrix of the forward modeling are exploited to decrease storage and computation time. The numerical tests demonstrate that the computation time of the proposed algorithm is reduced by approximately 2 orders of magnitude, and the memory requirement is reduced by <jats:italic>N</jats:italic>′<jats:sub><jats:italic>λ</jats:italic></jats:sub> times compared with the traditional GLQ method, where <jats:italic>N</jats:italic>′<jats:sub><jats:italic>λ</jats:italic></jats:sub> is the number of the model elements in the longitudinal direction. These significant improvements in computational efficiency and storage make it possible to calculate and store the dense Jacobian matrix in 3‐D large‐scale gravity inversions. The equivalence relations can be applied to the Taylor series method or combined with the adaptive discretization to ensure high accuracy. To further illustrate the capability of the algorithm, we present a regional synthetic example. The inverted results show density distributions consistent with the actual model. The computation took about 6.3 hr and 0.88 GB of memory compared with about a dozen days and 245.86 GB for the traditional 3‐D GLQ method. Finally, the proposed algorithm is applied to the gravity field derived from the latest lunar gravity model GL1500E. Three‐dimensional density distributions of the Imbrium and Serenitatis basins are obtained, and high‐density bodies are found at the depths 10–60 km, likely indicating a significant uplift of the high‐density mantle beneath the two mascon basins.</jats:p>

Item Type: Article
Uncontrolled Keywords: gravitational field, spherical coordinate, forward modeling, inversion
Depositing User: Symplectic Admin
Date Deposited: 12 Oct 2020 09:09
Last Modified: 04 Sep 2023 18:17
DOI: 10.1029/2019JB017691
Open Access URL: https://doi.org/10.31223/osf.io/dzf9j
Related URLs:
URI: https://livrepository.liverpool.ac.uk/id/eprint/3071061

Available Versions of this Item