Wavelet-Based Three-Dimensional Inversion for Geomagnetic Depth Sounding

The complexity of Earth’s structure poses a challenge to the multiscale detection capability of geophysics. In this paper, we present a new wavelet-based three-dimensional inversion method for geomagnetic depth sounding. This method is based on wavelet functions to transfer model parameters in the space domain into the wavelet domain. The model is represented by wavelet coefficients containing both large- and fine-scale information, enabling wavelet-based inversion to describe multiscale anomalies. L1-norm measurement is applied to measure the model roughness to accomplish the sparsity constraint in the wavelet domain. Meanwhile, a staggered-grid finite difference method in a spherical coordinate system is used to calculate the forward responses, and the limited-memory quasi-Newton method is applied to seek the solution of the inversion objective function. Inversion tests of synthetic data for multiscale models show that wavelet-based inversion is stable and has multiresolution. Although higher-order wavelets can lead to finer results, our tests present that a db6 wavelet is suitable for geomagnetic depth sounding inversion. The db6 inversion results of responses at 129 geomagnetic observatories around the world reveal a higher-resolution image of the mantle.