Three-Dimensional Magnetotelluric Forward Modeling in the Mixed Space-Wavenumber Domain with Divergence Correction
- 1School of Geosciences and Info-Physics, Central South University, Changsha 410083, China
Abstract
Three-dimensional magnetotelluric (MT) forward modeling remains constrained by the difficulty of achieving both high numerical accuracy and computational efficiency. To address this limitation, we propose a mixed space-wavenumber domain method based on coupled equations for the magnetic vector potential and electric scalar potential under the Lorenz gauge. Our method incorporates a divergence correction scheme and employs the Gaussian Fourier transform technique to solve the three-dimensional MT forward problem. The vector-potential equation reduces the number of unknowns and computational cost, while simplifying the treatment of boundary conditions at conductivity interfaces. The divergence correction suppresses spurious solutions at anomalous-body boundaries in low-induction-number regimes. We use Gauss-FFT instead of the traditional FFT to reduce the errors due to aliasing and imposed periodicity as well as edge effects in the MT forward modeling. Numerical experiments demonstrate that both iterative solution time and memory consumption scale linearly with the number of computational nodes. On refined meshes, the proposed method requires approximately one-eighth of the iterative solution time and one-quarter of the memory required by the finite element method. Together, these improvements substantially enhance computational efficiency without compromising numerical accuracy. The method is therefore well suited for large-scale 3D MT forward modeling on laterally regular grids with layered or vertically heterogeneous media, as required in high-resolution mineral exploration and regional crustal-structure imaging, where both accuracy and efficiency are essential.
Related
- The remote reference as an instrumental variable in model space: time-series magnetotelluric inversion with noisy magnetic fields
- Latent Diffusion Modeling of Two- and Three-dimensional Fracture Networks with Sparse Observation Conditioning
- Three-Dimensional Numerical Modeling of Glacier Flow and U-Shaped Valley Evolution in Subtropical Taiwan
- A general lightweight global modeling framework for three-dimensional seismic exploration
- Insights into Atlantic-Pacific divergence in Arbacia (Echinoidea) using revised systematics of Neogene fossils from the Caribbean and southeastern United States
BibTeX
@article{tong2026threedimensional,
title = {{Three-Dimensional Magnetotelluric Forward Modeling in the Mixed Space-Wavenumber Domain with Divergence Correction}},
author = {Tiegang Tong and Yufei Yang and Diquan Li},
journal = {Mathematics},
year = {2026},
month = oct,
doi = {10.3390/math14213918},
url = {https://www.mdpi.com/2227-7390/14/21/3918},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Three-Dimensional Magnetotelluric Forward Modeling in the Mixed Space-Wavenumber Domain with Divergence Correction
AU - Tiegang Tong
AU - Yufei Yang
AU - Diquan Li
PY - 2026
DA - 2026/10/27
JO - Mathematics
DO - 10.3390/math14213918
UR - https://www.mdpi.com/2227-7390/14/21/3918
AB - Three-dimensional magnetotelluric (MT) forward modeling remains constrained by the difficulty of achieving both high numerical accuracy and computational efficiency. To address this limitation, we propose a mixed space-wavenumber domain method based on coupled equations for the magnetic vector potential and electric scalar potential under the Lorenz gauge. Our method incorporates a divergence correction scheme and employs the Gaussian Fourier transform technique to solve the three-dimensional MT forward problem. The vector-potential equation reduces the number of unknowns and computational cost, while simplifying the treatment of boundary conditions at conductivity interfaces. The divergence correction suppresses spurious solutions at anomalous-body boundaries in low-induction-number regimes. We use Gauss-FFT instead of the traditional FFT to reduce the errors due to aliasing and imposed periodicity as well as edge effects in the MT forward modeling. Numerical experiments demonstrate that both iterative solution time and memory consumption scale linearly with the number of computational nodes. On refined meshes, the proposed method requires approximately one-eighth of the iterative solution time and one-quarter of the memory required by the finite element method. Together, these improvements substantially enhance computational efficiency without compromising numerical accuracy. The method is therefore well suited for large-scale 3D MT forward modeling on laterally regular grids with layered or vertically heterogeneous media, as required in high-resolution mineral exploration and regional crustal-structure imaging, where both accuracy and efficiency are essential.
SN - 2227-7390
ER -