Fast 3D inversion of gravity data using Lanczos bidiagonalization method

Reza Toushmalani, Hakim Saibi

Research output: Contribution to journalArticlepeer-review

22 Citations (Scopus)


This paper describes the application of a new inversion method to recover a three-dimensional density model from measured gravity anomalies. To attain this, the survey area is divided into a large number of rectangular prisms in a mesh with unknown densities. The results show that the application of the Lanczos bidiagonalization algorithm in the inversion helps to solve a Tikhonov cost function in a short time. The performance time of the inverse modeling greatly decreases by substituting the forward operator matrix with a matrix of lower dimension. A least-squares QR (LSQR) method is applied to select the best value of a regularization parameter. A Euler deconvolution method was used to avoid the natural trend of gravity structures to concentrate at shallow depth. Finally, the newly developed method was applied to synthetic data to demonstrate its suitability and then to real data from the Bandar Charak region (Hormozgan, south Iran). The 3D gravity inversion results were able to detect the location of the known salt dome (density contrast of −0.2 g/cm3) intrusion in the study area.

Original languageEnglish
Pages (from-to)4969-4981
Number of pages13
JournalArabian Journal of Geosciences
Issue number7
Publication statusPublished - Jul 2 2015
Externally publishedYes


  • 3D modeling
  • Bandar Charak
  • Gravity anomaly
  • Inverse problem
  • Lanczos bidiagonalization

ASJC Scopus subject areas

  • General Environmental Science
  • General Earth and Planetary Sciences


Dive into the research topics of 'Fast 3D inversion of gravity data using Lanczos bidiagonalization method'. Together they form a unique fingerprint.

Cite this