Spherical Planting Inversion of GRAIL Data

In large-scale potential field data inversion, constructing the kernel matrix is a time-consuming problem with large memory requirements. Therefore, a spherical planting inversion of Gravity Recovery and Interior Laboratory (GRAIL) data is proposed using the L1-norm in conjunction with tesseroids. S...

Full description

Bibliographic Details
Main Authors: Guangyin Lu, Dongxing Zhang, Shujin Cao, Yihuai Deng, Gang Xu, Yihu Liu, Ziqiang Zhu, Peng Chen
Format: Article
Language:English
Published: MDPI AG 2023-03-01
Series:Applied Sciences
Subjects:
Online Access:https://www.mdpi.com/2076-3417/13/5/3332
_version_ 1797615657665167360
author Guangyin Lu
Dongxing Zhang
Shujin Cao
Yihuai Deng
Gang Xu
Yihu Liu
Ziqiang Zhu
Peng Chen
author_facet Guangyin Lu
Dongxing Zhang
Shujin Cao
Yihuai Deng
Gang Xu
Yihu Liu
Ziqiang Zhu
Peng Chen
author_sort Guangyin Lu
collection DOAJ
description In large-scale potential field data inversion, constructing the kernel matrix is a time-consuming problem with large memory requirements. Therefore, a spherical planting inversion of Gravity Recovery and Interior Laboratory (GRAIL) data is proposed using the L1-norm in conjunction with tesseroids. Spherical planting inversion, however, is strongly dependent on the correct seeds’ density contrast, location, and number; otherwise, it can cause mutual intrusion of anomalous sources produced by different seeds. Hence, a weighting function was introduced to limit the influence area of the seeds for yielding robust solutions; moreover, it is challenging to set customized parameters for each seed, especially for the large number of seeds used or complex gravity anomalies data. Hence, we employed the “shape-of-anomaly” data-misfit function in conjunction with a new seed weighting function to improve the spherical planting inversion. The proposed seed weighting function is constructed based on the covariance matrix for given gravity data and can avoid manually setting customized parameters for each seed. The results of synthetic tests and field data show that spherical planting inversion requires less computer memory than traditional inversion. Furthermore, the proposed seed weighting function can effectively limit the seed influence area. The result of spherical planting inversion indicates that the crustal thickness of Mare Crisium is about 0 km because the Crisium impact may have removed all crust from parts of the basin.
first_indexed 2024-03-11T07:29:48Z
format Article
id doaj.art-76c91239a63c4c539a7c65ee59b2760e
institution Directory Open Access Journal
issn 2076-3417
language English
last_indexed 2024-03-11T07:29:48Z
publishDate 2023-03-01
publisher MDPI AG
record_format Article
series Applied Sciences
spelling doaj.art-76c91239a63c4c539a7c65ee59b2760e2023-11-17T07:22:25ZengMDPI AGApplied Sciences2076-34172023-03-01135333210.3390/app13053332Spherical Planting Inversion of GRAIL DataGuangyin Lu0Dongxing Zhang1Shujin Cao2Yihuai Deng3Gang Xu4Yihu Liu5Ziqiang Zhu6Peng Chen7School of Geosciences and Info-Physics, Central South University, Changsha 410083, ChinaSchool of Geosciences and Info-Physics, Central South University, Changsha 410083, ChinaSchool of Geosciences and Info-Physics, Central South University, Changsha 410083, ChinaSchool of Earth Sciences and Spatial Information Engineering, Hunan University of Science and Technology, Xiangtan 411201, ChinaSchool of Geosciences and Info-Physics, Central South University, Changsha 410083, ChinaSchool of Geosciences and Info-Physics, Central South University, Changsha 410083, ChinaSchool of Geosciences and Info-Physics, Central South University, Changsha 410083, ChinaSchool of Earth Sciences and Spatial Information Engineering, Hunan University of Science and Technology, Xiangtan 411201, ChinaIn large-scale potential field data inversion, constructing the kernel matrix is a time-consuming problem with large memory requirements. Therefore, a spherical planting inversion of Gravity Recovery and Interior Laboratory (GRAIL) data is proposed using the L1-norm in conjunction with tesseroids. Spherical planting inversion, however, is strongly dependent on the correct seeds’ density contrast, location, and number; otherwise, it can cause mutual intrusion of anomalous sources produced by different seeds. Hence, a weighting function was introduced to limit the influence area of the seeds for yielding robust solutions; moreover, it is challenging to set customized parameters for each seed, especially for the large number of seeds used or complex gravity anomalies data. Hence, we employed the “shape-of-anomaly” data-misfit function in conjunction with a new seed weighting function to improve the spherical planting inversion. The proposed seed weighting function is constructed based on the covariance matrix for given gravity data and can avoid manually setting customized parameters for each seed. The results of synthetic tests and field data show that spherical planting inversion requires less computer memory than traditional inversion. Furthermore, the proposed seed weighting function can effectively limit the seed influence area. The result of spherical planting inversion indicates that the crustal thickness of Mare Crisium is about 0 km because the Crisium impact may have removed all crust from parts of the basin.https://www.mdpi.com/2076-3417/13/5/3332spherical planting inversiontesseroidsGRAIL dataMare Crisiumseed weighting function
spellingShingle Guangyin Lu
Dongxing Zhang
Shujin Cao
Yihuai Deng
Gang Xu
Yihu Liu
Ziqiang Zhu
Peng Chen
Spherical Planting Inversion of GRAIL Data
Applied Sciences
spherical planting inversion
tesseroids
GRAIL data
Mare Crisium
seed weighting function
title Spherical Planting Inversion of GRAIL Data
title_full Spherical Planting Inversion of GRAIL Data
title_fullStr Spherical Planting Inversion of GRAIL Data
title_full_unstemmed Spherical Planting Inversion of GRAIL Data
title_short Spherical Planting Inversion of GRAIL Data
title_sort spherical planting inversion of grail data
topic spherical planting inversion
tesseroids
GRAIL data
Mare Crisium
seed weighting function
url https://www.mdpi.com/2076-3417/13/5/3332
work_keys_str_mv AT guangyinlu sphericalplantinginversionofgraildata
AT dongxingzhang sphericalplantinginversionofgraildata
AT shujincao sphericalplantinginversionofgraildata
AT yihuaideng sphericalplantinginversionofgraildata
AT gangxu sphericalplantinginversionofgraildata
AT yihuliu sphericalplantinginversionofgraildata
AT ziqiangzhu sphericalplantinginversionofgraildata
AT pengchen sphericalplantinginversionofgraildata