Spectral-infinite-element simulations of magnetic anomalies

Hom Nath Gharti, Jeroen Tromp

Research output: Contribution to journalArticlepeer-review

7 Scopus citations


SUMMARY We implement a spectral-infinite-element method (SIEM) to compute magnetic anomalies by solving a discretized form of the Poisson/Laplace equation. The SIEM combines the highly accurate spectral-element method with the mapped-infinite element method, which reproduces an unbounded domain accurately and efficiently. This combination is made possible by coupling Gauss-Legendre-Lobatto quadrature in spectral elements with Gauss-Radau quadrature in infinite elements along the infinite directions. Our method has two distinct advantages over traditional methods. First, the higher-order discretization accurately renders complex magnetized heterogeneities. Second, since the computation time is independent of the number of observation points, the method is efficient for very large models. We illustrate the accuracy and efficiency of our method by comparing calculated magnetic anomalies for various magnetized heterogeneities with corresponding analytical and commonly used computational solutions. We conclude with a practical example involving a complex 3-D model of an ore mine.

Original languageEnglish (US)
Pages (from-to)1656-1667
Number of pages12
JournalGeophysical Journal International
Issue number3
StatePublished - Feb 27 2019

All Science Journal Classification (ASJC) codes

  • Geophysics
  • Geochemistry and Petrology


  • Magnetic anomalies: modelling and interpretation
  • Magnetic properties
  • Numerical modelling
  • Rock and mineral magnetism


Dive into the research topics of 'Spectral-infinite-element simulations of magnetic anomalies'. Together they form a unique fingerprint.

Cite this