Submitted as: development and technical paper 30 Mar 2021

Submitted as: development and technical paper | 30 Mar 2021

Review status: a revised version of this preprint was accepted for the journal GMD.

Structural, petrophysical and geological constraints in potential field inversion using the Tomofast-x v1.0 open-source code

Jeremie Giraud1,2, Vitaliy Ogarko3,4,*, Roland Martin5, Mark Jessell1,2, and Mark Lindsay1,2 Jeremie Giraud et al.
  • 1Centre for Exploration Targeting (School of Earth Sciences), University of Western Australia, 35 Stirling Highway, 6009 Crawley, Australia
  • 2Mineral Exploration Cooperative Research Centre, School of Earth Sciences, University of Western Australia, 35 Stirling Highway, WA Crawley 6009, Australia
  • 3The International Centre for Radio Astronomy Research, University of Western Australia, 7 Fairway, WA Crawley 6009, Australia
  • 4ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D) Australia
  • 5Laboratoire de Géosciences Environnement Toulouse GET, CNRS UMR 5563, Observatoire Midi-Pyrénées, Université Paul Sabatier, 14 avenue Edouard Belin, 31400 Toulouse, France
  • *Formerly at Centre for Exploration Targeting (School of Earth Sciences), University of Western Australia, 35 Stirling Highway, 6009 Crawley, Australia

Abstract. The quantitative integration of geophysical measurements with data and information from other disciplines is becoming increasingly important in answering the challenges of undercover imaging and of the modelling of complex areas. We propose a review of the different techniques for the utilisation of structural, petrophysical and geological information in single physics and joint inversion as implemented in the Tomofast-x open-source inversion platform. We detail the range of constraints that can be applied to the inversion of potential field data. The inversion examples we show illustrate a selection of scenarios using a realistic synthetic dataset inspired by real-world geological measurements and petrophysical data from the Hamersley region (Western Australia). Using Tomofast-x’s flexibility, we investigate inversions combining the utilisation of petrophysical, structural and/or geological constraints while illustrating the utilisation of the L-curve principle to determine regularisation weights. Our results suggest that the utilisation of geological information to derive disjoint interval bound constraints is the most effective method to recover the true model. It is followed by model smoothness and smallness conditioned by geological uncertainty, and cross-gradient minimisation.

Jeremie Giraud et al.

Status: final response (author comments only)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on gmd-2021-14', Anonymous Referee #1, 07 May 2021
    • AC1: 'Reply on RC1', Jeremie Giraud, 16 Jul 2021
  • RC2: 'Comment on gmd-2021-14', Mehrdad Bastani, 11 Jun 2021
    • AC2: 'Reply on RC2', Jeremie Giraud, 16 Jul 2021

Jeremie Giraud et al.

Jeremie Giraud et al.


Total article views: 499 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
334 155 10 499 5 5
  • HTML: 334
  • PDF: 155
  • XML: 10
  • Total: 499
  • BibTeX: 5
  • EndNote: 5
Views and downloads (calculated since 30 Mar 2021)
Cumulative views and downloads (calculated since 30 Mar 2021)

Viewed (geographical distribution)

Total article views: 444 (including HTML, PDF, and XML) Thereof 444 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 23 Jul 2021
Short summary
We review different techniques to model the Earth's subsurface from geophysical data (gravity field anomaly, magnetic field anomaly) using geological models and measurements of the rocks' properties. We show examples of application using idealised examples reproducing realistic features, while providing theoreitical details of the open-source algorithm we use.