A Factorized Green-Operator Framework for Efficient Wavenumber-Domain Inversion of Magnetic Anomalies and Gradient Tensors
Clicks: 1
ID: 319000
2026
Article Quality & Performance Metrics
Overall Quality
Not rated
Combines reader engagement with the AI quality analysis. This
article has not been analysed, so there is no overall score —
reader engagement is measured and shown alongside.
Reader Engagement
0.0
/100
1 views
0 readers
AI Quality Assessment
Not analyzed
Readership in this journal
Ranked #210 of 220 articles by views in geophysical journal international
Most read
Least read
Bar heights use a square-root scale. Only the 120 most-read articles are drawn; the journal has 220 in total.
Mint this article as an NFT
Not yet mintedCreate a permanent, verifiable on-chain record of this article on the Scimatic Network. The NFT is held in your Journament account, and you can withdraw it to your own wallet at any time.
5
SUSD
one-off · no wallet required
Abstract
Summary Magnetic inversion is a key tool for imaging subsurface geological structures, but conventional 3-D magnetic inversion in the spatial domain is often limited by the computational and memory cost of large dense kernel matrices. Existing transformed-domain approaches improve efficiency, yet pseudo-3D implementations still rely on layer-by-layer accumulation and repeated Fourier transforms. In this study, we develop a unified wavenumber-domain framework for the forward modelling and inversion of total-field magnetic anomalies and magnetic gradient-tensor data. For regularly discretized rectangular prisms beneath a planar observation surface, the wavenumber-domain Green operator is reformulated into a factorized representation consisting of two explicitly stored diagonal/block-diagonal spectral factors and one implicitly applied separable horizontal operator. This implementation avoids repeated vertical layer superposition and reduces the forward evaluation to a single FFT/IFFT pair together with structured spectral multiplications. The factorized forward operator is then embedded in a Tikhonov-regularized inversion and solved through a Sherman-Morrison-Woodbury (SMW) reduced system. The transformed-domain data term is defined as an unweighted complex-valued least-squares residual, and its relation to the spatial-domain least-squares formulation is stated under the corresponding padding and truncation assumptions. Synthetic examples show that the method reproduces conventional spatial-domain responses and recovers the principal features of prescribed magnetization models under 5% Gaussian noise. For a 200×200×100 model, the forward modeling and core inversion times are 0.172 s and 31.73 s, respectively, on a standard laptop. Application to field data is used as a practical feasibility test and shows a data-consistent recovered magnetization distribution, but it should not be regarded as an independent geological validation of the recovered model. The current implementation assumes a planar observation surface, a regular FFT-compatible grid, and a spatially uniform magnetization direction. It does not yet address strong remanence, spatially variable magnetization, irregular topography, irregular acquisition geometries, depth weighting, focusing stabilizers, or geological constraints. Under these assumptions, the proposed framework provides an efficient, memory-economical, and scalable alternative for large-scale magnetic anomaly interpretation.
| Reference Key |
openalex_W7166298850
Use this key to autocite in the manuscript while using
SciMatic Manuscript Manager or Thesis Manager
|
|---|---|
| Authors | Sheng Liu, Shi Chen, Zhanfeng Huang, Yiju Tang, Fangchao Lu, Songbai Xuan, Shuanggen Jin |
| Journal | geophysical journal international |
| Year | 2026 |
| DOI |
10.1093/gji/ggag250
|
| URL | |
| Keywords | Keywords not found |
Citations
No citations found. To add a citation, contact the admin at info@scimatic.org
Comments
No comments yet. Be the first to comment on this article.