Articles | Volume 14, issue 11
https://doi.org/10.5194/gmd-14-6681-2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-14-6681-2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
02 Nov 2021
Development and technical paper |
| 02 Nov 2021
| 02 Nov 2021
Structural, petrophysical, and geological constraints in potential field inversion using the Tomofast-x v1.0 open-source code
Centre for Exploration Targeting (School of Earth Sciences), University of Western Australia, 35 Stirling Highway, 6009 Crawley, WA, Australia
Mineral Exploration Cooperative Research Centre, School of Earth Sciences, University of Western Australia, 35 Stirling Highway, 6009 Crawley, WA, Australia
now at: GeoRessources, Université de Lorraine, CNRS, 54000, Nancy, France
Vitaliy Ogarko
Centre for Exploration Targeting (School of Earth Sciences), University of Western Australia, 35 Stirling Highway, 6009 Crawley, WA, Australia
Mineral Exploration Cooperative Research Centre, School of Earth Sciences, University of Western Australia, 35 Stirling Highway, 6009 Crawley, WA, Australia
The International Centre for Radio Astronomy Research, University of Western Australia, 7 Fairway, 6009 Crawley, WA, Australia
Roland Martin
Laboratoire de Géosciences Environnement Toulouse GET, CNRS UMR 5563, Observatoire Midi-Pyrénées, Université Paul Sabatier, 14 avenue Edouard Belin, 31400, Toulouse, France
Mark Jessell
Centre for Exploration Targeting (School of Earth Sciences), University of Western Australia, 35 Stirling Highway, 6009 Crawley, WA, Australia
Mineral Exploration Cooperative Research Centre, School of Earth Sciences, University of Western Australia, 35 Stirling Highway, 6009 Crawley, WA, Australia
Mark Lindsay
Centre for Exploration Targeting (School of Earth Sciences), University of Western Australia, 35 Stirling Highway, 6009 Crawley, WA, Australia
Mineral Exploration Cooperative Research Centre, School of Earth Sciences, University of Western Australia, 35 Stirling Highway, 6009 Crawley, WA, Australia
Related authors
Lawrence A. Bird, Vitaliy Ogarko, Laurent Ailleres, Lachlan Grose, Jeremie Giraud, Felicity S. McCormack, David E. Gwyther, Jason L. Roberts, Richard S. Jones, and Andrew N. Mackintosh
EGUsphere, https://doi.org/10.5194/egusphere-2025-211, https://doi.org/10.5194/egusphere-2025-211, 2025
Short summary
Short summary
The terrain of the seafloor has important controls on the access of warm water to below floating ice shelves around Antarctica. Here, we present an open-source method to infer what the seafloor looks like around the Antarctic continent, and within these ice shelf cavities, using measurements of the Earth’s gravitational field. We present an improved seafloor map for the Vincennes Bay region in East Antarctica and assess its impact on ice melt rates.
Mark Douglas Lindsay, Vitaliy Ogarko, Jeremie Giraud, and Mosayeb Khademi
EGUsphere, https://doi.org/10.5194/egusphere-2024-3754, https://doi.org/10.5194/egusphere-2024-3754, 2025
Preprint archived
Short summary
Short summary
Geophysical data is used to understand the Earth's structure. Geophysical inversion is an optimisation technique that constructs a model of the Earth using geophysical data. Inversion is complex and requires constraints to produce a plausible model. Overfitting is a disadvantageous effect that can affect techniques that rely on optimisation. A new feature in an inversion platform reduces overfitting by informing inversion of spatial uncertainty contained within the geophysical data.
Vitaliy Ogarko, Kim Frankcombe, Taige Liu, Jeremie Giraud, Roland Martin, and Mark Jessell
Geosci. Model Dev., 17, 2325–2345, https://doi.org/10.5194/gmd-17-2325-2024, https://doi.org/10.5194/gmd-17-2325-2024, 2024
Short summary
Short summary
We present a major release of the Tomofast-x open-source gravity and magnetic inversion code that is enhancing its performance and applicability for both industrial and academic studies. We focus on real-world mineral exploration scenarios, while offering flexibility for applications at regional scale or for crustal studies. The optimisation work described in this paper is fundamental to allowing more complete descriptions of the controls on magnetisation, including remanence.
Jérémie Giraud, Guillaume Caumon, Lachlan Grose, Vitaliy Ogarko, and Paul Cupillard
Solid Earth, 15, 63–89, https://doi.org/10.5194/se-15-63-2024, https://doi.org/10.5194/se-15-63-2024, 2024
Short summary
Short summary
We present and test an algorithm that integrates geological modelling into deterministic geophysical inversion. This is motivated by the need to model the Earth using all available data and to reconcile the different types of measurements. We introduce the methodology and test our algorithm using two idealised scenarios. Results suggest that the method we propose is effectively capable of improving the models recovered by geophysical inversion and may be applied in real-world scenarios.
Jérémie Giraud, Hoël Seillé, Mark D. Lindsay, Gerhard Visser, Vitaliy Ogarko, and Mark W. Jessell
Solid Earth, 14, 43–68, https://doi.org/10.5194/se-14-43-2023, https://doi.org/10.5194/se-14-43-2023, 2023
Short summary
Short summary
We propose and apply a workflow to combine the modelling and interpretation of magnetic anomalies and resistivity anomalies to better image the basement. We test the method on a synthetic case study and apply it to real world data from the Cloncurry area (Queensland, Australia), which is prospective for economic minerals. Results suggest a new interpretation of the composition and structure towards to east of the profile that we modelled.
Guillaume Pirot, Ranee Joshi, Jérémie Giraud, Mark Douglas Lindsay, and Mark Walter Jessell
Geosci. Model Dev., 15, 4689–4708, https://doi.org/10.5194/gmd-15-4689-2022, https://doi.org/10.5194/gmd-15-4689-2022, 2022
Short summary
Short summary
Results of a survey launched among practitioners in the mineral industry show that despite recognising the importance of uncertainty quantification it is not very well performed due to lack of data, time requirements, poor tracking of interpretations and relative complexity of uncertainty quantification. To alleviate the latter, we provide an open-source set of local and global indicators to measure geological uncertainty among an ensemble of geological models.
Richard Scalzo, Mark Lindsay, Mark Jessell, Guillaume Pirot, Jeremie Giraud, Edward Cripps, and Sally Cripps
Geosci. Model Dev., 15, 3641–3662, https://doi.org/10.5194/gmd-15-3641-2022, https://doi.org/10.5194/gmd-15-3641-2022, 2022
Short summary
Short summary
This paper addresses numerical challenges in reasoning about geological models constrained by sensor data, especially models that describe the history of an area in terms of a sequence of events. Our method ensures that small changes in simulated geological features, such as the position of a boundary between two rock layers, do not result in unrealistically large changes to resulting sensor measurements, as occur presently using several popular modeling packages.
Mark Jessell, Jiateng Guo, Yunqiang Li, Mark Lindsay, Richard Scalzo, Jérémie Giraud, Guillaume Pirot, Ed Cripps, and Vitaliy Ogarko
Earth Syst. Sci. Data, 14, 381–392, https://doi.org/10.5194/essd-14-381-2022, https://doi.org/10.5194/essd-14-381-2022, 2022
Short summary
Short summary
To robustly train and test automated methods in the geosciences, we need to have access to large numbers of examples where we know
the answer. We present a suite of synthetic 3D geological models with their gravity and magnetic responses that allow researchers to test their methods on a whole range of geologically plausible models, thus overcoming one of the fundamental limitations of automation studies.
Mahtab Rashidifard, Jérémie Giraud, Mark Lindsay, Mark Jessell, and Vitaliy Ogarko
Solid Earth, 12, 2387–2406, https://doi.org/10.5194/se-12-2387-2021, https://doi.org/10.5194/se-12-2387-2021, 2021
Short summary
Short summary
One motivation for this study is to develop a workflow that enables the integration of geophysical datasets with different coverages that are quite common in exploration geophysics. We have utilized a level set approach to achieve this goal. The utilized technique parameterizes the subsurface in the same fashion as geological models. Our results indicate that the approach is capable of integrating information from seismic data in 2D to guide the 3D inversion results of the gravity data.
Jérémie Giraud, Mark Lindsay, Mark Jessell, and Vitaliy Ogarko
Solid Earth, 11, 419–436, https://doi.org/10.5194/se-11-419-2020, https://doi.org/10.5194/se-11-419-2020, 2020
Short summary
Short summary
We propose a methodology for the identification of rock types using geophysical and geological information. It relies on an algorithm used in machine learning called
self-organizing maps, to which we add plausibility filters to ensure that the results respect base geological rules and geophysical measurements. Application in the Yerrida Basin (Western Australia) reveals that the thinning of prospective greenstone belts at depth could be due to deep structures not seen from surface.
Evren Pakyuz-Charrier, Mark Jessell, Jérémie Giraud, Mark Lindsay, and Vitaliy Ogarko
Solid Earth, 10, 1663–1684, https://doi.org/10.5194/se-10-1663-2019, https://doi.org/10.5194/se-10-1663-2019, 2019
Short summary
Short summary
This paper improves the Monte Carlo simulation for uncertainty propagation (MCUP) method for 3-D geological modeling. Topological heterogeneity is observed in the model suite. The study demonstrates that such heterogeneity arises from piecewise nonlinearity inherent to 3-D geological models and contraindicates use of global uncertainty estimation methods. Topological-clustering-driven uncertainty estimation is proposed as a demonstrated alternative to address plausible model heterogeneity.
Jeremie Giraud, Mark Lindsay, Vitaliy Ogarko, Mark Jessell, Roland Martin, and Evren Pakyuz-Charrier
Solid Earth, 10, 193–210, https://doi.org/10.5194/se-10-193-2019, https://doi.org/10.5194/se-10-193-2019, 2019
Short summary
Short summary
We propose the quantitative integration of geology and geophysics in an algorithm integrating the probability of observation of rocks with gravity data to improve subsurface imaging. This allows geophysical modelling to adjust models preferentially in the least certain areas while honouring geological information and geophysical data. We validate our algorithm using an idealized case and apply it to the Yerrida Basin (Australia), where we can recover the geometry of buried greenstone belts.
Evren Pakyuz-Charrier, Mark Lindsay, Vitaliy Ogarko, Jeremie Giraud, and Mark Jessell
Solid Earth, 9, 385–402, https://doi.org/10.5194/se-9-385-2018, https://doi.org/10.5194/se-9-385-2018, 2018
Short summary
Short summary
MCUE is a method that produces probabilistic 3-D geological models by sampling from distributions that represent the uncertainty of the initial input dataset. This process generates numerous plausible datasets used to produce a range of statistically plausible 3-D models which are combined into a single probabilistic model. In this paper, improvements to distribution selection and parameterization for input uncertainty are proposed.
Vitaliy Ogarko and Mark Jessell
EGUsphere, https://doi.org/10.5194/egusphere-2025-1294, https://doi.org/10.5194/egusphere-2025-1294, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
We developed a new method to reconstruct underground rock layers from drillhole data, using an advanced algorithm to ensure geologically realistic results. By combining data from multiple drillholes, our approach reduces uncertainty and improves accuracy. Tested on South Australian data, it successfully predicted stratigraphy and highlighted ways to enhance data quality. This innovation makes geological analysis more reliable, aiding exploration and resource management.
Lawrence A. Bird, Vitaliy Ogarko, Laurent Ailleres, Lachlan Grose, Jeremie Giraud, Felicity S. McCormack, David E. Gwyther, Jason L. Roberts, Richard S. Jones, and Andrew N. Mackintosh
EGUsphere, https://doi.org/10.5194/egusphere-2025-211, https://doi.org/10.5194/egusphere-2025-211, 2025
Short summary
Short summary
The terrain of the seafloor has important controls on the access of warm water to below floating ice shelves around Antarctica. Here, we present an open-source method to infer what the seafloor looks like around the Antarctic continent, and within these ice shelf cavities, using measurements of the Earth’s gravitational field. We present an improved seafloor map for the Vincennes Bay region in East Antarctica and assess its impact on ice melt rates.
Mark Douglas Lindsay, Vitaliy Ogarko, Jeremie Giraud, and Mosayeb Khademi
EGUsphere, https://doi.org/10.5194/egusphere-2024-3754, https://doi.org/10.5194/egusphere-2024-3754, 2025
Preprint archived
Short summary
Short summary
Geophysical data is used to understand the Earth's structure. Geophysical inversion is an optimisation technique that constructs a model of the Earth using geophysical data. Inversion is complex and requires constraints to produce a plausible model. Overfitting is a disadvantageous effect that can affect techniques that rely on optimisation. A new feature in an inversion platform reduces overfitting by informing inversion of spatial uncertainty contained within the geophysical data.
Léonard Moracchini, Guillaume Pirot, Kerry Bardot, Mark W. Jessell, and James L. McCallum
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-154, https://doi.org/10.5194/gmd-2024-154, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
To facilitate the exploration of alternative hydrogeological scenarios, we propose to approximate costly physical simulations of contaminant transport by more affordable shortest distances computations. It enables to accept or reject scenarios within a predefined confidence interval. In particular, it can allow to estimate the probability of a fault acting as a preferential path or a barrier.
Vitaliy Ogarko, Kim Frankcombe, Taige Liu, Jeremie Giraud, Roland Martin, and Mark Jessell
Geosci. Model Dev., 17, 2325–2345, https://doi.org/10.5194/gmd-17-2325-2024, https://doi.org/10.5194/gmd-17-2325-2024, 2024
Short summary
Short summary
We present a major release of the Tomofast-x open-source gravity and magnetic inversion code that is enhancing its performance and applicability for both industrial and academic studies. We focus on real-world mineral exploration scenarios, while offering flexibility for applications at regional scale or for crustal studies. The optimisation work described in this paper is fundamental to allowing more complete descriptions of the controls on magnetisation, including remanence.
Jiateng Guo, Xuechuang Xu, Luyuan Wang, Xulei Wang, Lixin Wu, Mark Jessell, Vitaliy Ogarko, Zhibin Liu, and Yufei Zheng
Geosci. Model Dev., 17, 957–973, https://doi.org/10.5194/gmd-17-957-2024, https://doi.org/10.5194/gmd-17-957-2024, 2024
Short summary
Short summary
This study proposes a semi-supervised learning algorithm using pseudo-labels for 3D geological modelling. We establish a 3D geological model using borehole data from a complex real urban local survey area in Shenyang and make an uncertainty analysis of this model. The method effectively expands the sample space, which is suitable for geomodelling and uncertainty analysis from boreholes. The modelling results perform well in terms of spatial morphology and geological semantics.
Jérémie Giraud, Guillaume Caumon, Lachlan Grose, Vitaliy Ogarko, and Paul Cupillard
Solid Earth, 15, 63–89, https://doi.org/10.5194/se-15-63-2024, https://doi.org/10.5194/se-15-63-2024, 2024
Short summary
Short summary
We present and test an algorithm that integrates geological modelling into deterministic geophysical inversion. This is motivated by the need to model the Earth using all available data and to reconcile the different types of measurements. We introduce the methodology and test our algorithm using two idealised scenarios. Results suggest that the method we propose is effectively capable of improving the models recovered by geophysical inversion and may be applied in real-world scenarios.
Jérémie Giraud, Hoël Seillé, Mark D. Lindsay, Gerhard Visser, Vitaliy Ogarko, and Mark W. Jessell
Solid Earth, 14, 43–68, https://doi.org/10.5194/se-14-43-2023, https://doi.org/10.5194/se-14-43-2023, 2023
Short summary
Short summary
We propose and apply a workflow to combine the modelling and interpretation of magnetic anomalies and resistivity anomalies to better image the basement. We test the method on a synthetic case study and apply it to real world data from the Cloncurry area (Queensland, Australia), which is prospective for economic minerals. Results suggest a new interpretation of the composition and structure towards to east of the profile that we modelled.
Guillaume Pirot, Ranee Joshi, Jérémie Giraud, Mark Douglas Lindsay, and Mark Walter Jessell
Geosci. Model Dev., 15, 4689–4708, https://doi.org/10.5194/gmd-15-4689-2022, https://doi.org/10.5194/gmd-15-4689-2022, 2022
Short summary
Short summary
Results of a survey launched among practitioners in the mineral industry show that despite recognising the importance of uncertainty quantification it is not very well performed due to lack of data, time requirements, poor tracking of interpretations and relative complexity of uncertainty quantification. To alleviate the latter, we provide an open-source set of local and global indicators to measure geological uncertainty among an ensemble of geological models.
Richard Scalzo, Mark Lindsay, Mark Jessell, Guillaume Pirot, Jeremie Giraud, Edward Cripps, and Sally Cripps
Geosci. Model Dev., 15, 3641–3662, https://doi.org/10.5194/gmd-15-3641-2022, https://doi.org/10.5194/gmd-15-3641-2022, 2022
Short summary
Short summary
This paper addresses numerical challenges in reasoning about geological models constrained by sensor data, especially models that describe the history of an area in terms of a sequence of events. Our method ensures that small changes in simulated geological features, such as the position of a boundary between two rock layers, do not result in unrealistically large changes to resulting sensor measurements, as occur presently using several popular modeling packages.
Mark Jessell, Jiateng Guo, Yunqiang Li, Mark Lindsay, Richard Scalzo, Jérémie Giraud, Guillaume Pirot, Ed Cripps, and Vitaliy Ogarko
Earth Syst. Sci. Data, 14, 381–392, https://doi.org/10.5194/essd-14-381-2022, https://doi.org/10.5194/essd-14-381-2022, 2022
Short summary
Short summary
To robustly train and test automated methods in the geosciences, we need to have access to large numbers of examples where we know
the answer. We present a suite of synthetic 3D geological models with their gravity and magnetic responses that allow researchers to test their methods on a whole range of geologically plausible models, thus overcoming one of the fundamental limitations of automation studies.
Ranee Joshi, Kavitha Madaiah, Mark Jessell, Mark Lindsay, and Guillaume Pirot
Geosci. Model Dev., 14, 6711–6740, https://doi.org/10.5194/gmd-14-6711-2021, https://doi.org/10.5194/gmd-14-6711-2021, 2021
Short summary
Short summary
We have developed a software that allows the user to extract and standardize drill hole information from legacy datasets and/or different drilling campaigns. It also provides functionality to upscale the lithological information. These functionalities were possible by developing thesauri to identify and group geological terminologies together.
Mahtab Rashidifard, Jérémie Giraud, Mark Lindsay, Mark Jessell, and Vitaliy Ogarko
Solid Earth, 12, 2387–2406, https://doi.org/10.5194/se-12-2387-2021, https://doi.org/10.5194/se-12-2387-2021, 2021
Short summary
Short summary
One motivation for this study is to develop a workflow that enables the integration of geophysical datasets with different coverages that are quite common in exploration geophysics. We have utilized a level set approach to achieve this goal. The utilized technique parameterizes the subsurface in the same fashion as geological models. Our results indicate that the approach is capable of integrating information from seismic data in 2D to guide the 3D inversion results of the gravity data.
Lachlan Grose, Laurent Ailleres, Gautier Laurent, Guillaume Caumon, Mark Jessell, and Robin Armit
Geosci. Model Dev., 14, 6197–6213, https://doi.org/10.5194/gmd-14-6197-2021, https://doi.org/10.5194/gmd-14-6197-2021, 2021
Short summary
Short summary
Fault discontinuities in rock packages represent the plane where two blocks of rock have moved. They are challenging to incorporate into geological models because the geometry of the faulted rock units are defined by not only the location of the discontinuity but also the kinematics of the fault. In this paper, we outline a structural geology framework for incorporating faults into geological models by directly incorporating kinematics into the mathematical framework of the model.
Mark Jessell, Vitaliy Ogarko, Yohan de Rose, Mark Lindsay, Ranee Joshi, Agnieszka Piechocka, Lachlan Grose, Miguel de la Varga, Laurent Ailleres, and Guillaume Pirot
Geosci. Model Dev., 14, 5063–5092, https://doi.org/10.5194/gmd-14-5063-2021, https://doi.org/10.5194/gmd-14-5063-2021, 2021
Short summary
Short summary
We have developed software that allows the user to extract sufficient information from unmodified digital maps and associated datasets that we are able to use to automatically build 3D geological models. By automating the process we are able to remove human bias from the procedure, which makes the workflow reproducible.
Lachlan Grose, Laurent Ailleres, Gautier Laurent, and Mark Jessell
Geosci. Model Dev., 14, 3915–3937, https://doi.org/10.5194/gmd-14-3915-2021, https://doi.org/10.5194/gmd-14-3915-2021, 2021
Short summary
Short summary
LoopStructural is an open-source 3D geological modelling library with a model design allowing for multiple different algorithms to be used for comparison for the same geology. Geological structures are modelled using structural geology concepts and techniques, allowing for complex structures such as overprinted folds and faults to be modelled. In the paper, we demonstrate automatically generating a 3-D model from map2loop-processed geological survey data of the Flinders Ranges, South Australia.
Mark D. Lindsay, Sandra Occhipinti, Crystal Laflamme, Alan Aitken, and Lara Ramos
Solid Earth, 11, 1053–1077, https://doi.org/10.5194/se-11-1053-2020, https://doi.org/10.5194/se-11-1053-2020, 2020
Short summary
Short summary
Integrated interpretation of multiple datasets is a key skill required for better understanding the composition and configuration of the Earth's crust. Geophysical and 3D geological modelling are used here to aid the interpretation process in investigating anomalous and cryptic geophysical signatures which suggest a more complex structure and history of a Palaeoproterozoic basin in Western Australia.
Jérémie Giraud, Mark Lindsay, Mark Jessell, and Vitaliy Ogarko
Solid Earth, 11, 419–436, https://doi.org/10.5194/se-11-419-2020, https://doi.org/10.5194/se-11-419-2020, 2020
Short summary
Short summary
We propose a methodology for the identification of rock types using geophysical and geological information. It relies on an algorithm used in machine learning called
self-organizing maps, to which we add plausibility filters to ensure that the results respect base geological rules and geophysical measurements. Application in the Yerrida Basin (Western Australia) reveals that the thinning of prospective greenstone belts at depth could be due to deep structures not seen from surface.
Evren Pakyuz-Charrier, Mark Jessell, Jérémie Giraud, Mark Lindsay, and Vitaliy Ogarko
Solid Earth, 10, 1663–1684, https://doi.org/10.5194/se-10-1663-2019, https://doi.org/10.5194/se-10-1663-2019, 2019
Short summary
Short summary
This paper improves the Monte Carlo simulation for uncertainty propagation (MCUP) method for 3-D geological modeling. Topological heterogeneity is observed in the model suite. The study demonstrates that such heterogeneity arises from piecewise nonlinearity inherent to 3-D geological models and contraindicates use of global uncertainty estimation methods. Topological-clustering-driven uncertainty estimation is proposed as a demonstrated alternative to address plausible model heterogeneity.
Jeremie Giraud, Mark Lindsay, Vitaliy Ogarko, Mark Jessell, Roland Martin, and Evren Pakyuz-Charrier
Solid Earth, 10, 193–210, https://doi.org/10.5194/se-10-193-2019, https://doi.org/10.5194/se-10-193-2019, 2019
Short summary
Short summary
We propose the quantitative integration of geology and geophysics in an algorithm integrating the probability of observation of rocks with gravity data to improve subsurface imaging. This allows geophysical modelling to adjust models preferentially in the least certain areas while honouring geological information and geophysical data. We validate our algorithm using an idealized case and apply it to the Yerrida Basin (Australia), where we can recover the geometry of buried greenstone belts.
Evren Pakyuz-Charrier, Mark Lindsay, Vitaliy Ogarko, Jeremie Giraud, and Mark Jessell
Solid Earth, 9, 385–402, https://doi.org/10.5194/se-9-385-2018, https://doi.org/10.5194/se-9-385-2018, 2018
Short summary
Short summary
MCUE is a method that produces probabilistic 3-D geological models by sampling from distributions that represent the uncertainty of the initial input dataset. This process generates numerous plausible datasets used to produce a range of statistically plausible 3-D models which are combined into a single probabilistic model. In this paper, improvements to distribution selection and parameterization for input uncertainty are proposed.
Xiaojun Feng, Enyuan Wang, Jérôme Ganne, Roland Martin, and Mark W. Jessell
Solid Earth Discuss., https://doi.org/10.5194/se-2017-142, https://doi.org/10.5194/se-2017-142, 2018
Preprint withdrawn
J. Florian Wellmann, Sam T. Thiele, Mark D. Lindsay, and Mark W. Jessell
Geosci. Model Dev., 9, 1019–1035, https://doi.org/10.5194/gmd-9-1019-2016, https://doi.org/10.5194/gmd-9-1019-2016, 2016
Short summary
Short summary
We often obtain knowledge about the subsurface in the form of structural geological models, as a basis for subsurface usage or resource extraction. Here, we provide a modelling code to construct such models on the basis of significant deformational events in geological history, encapsulated in kinematic equations. Our methods simplify complex dynamic processes, but enable us to evaluate how events interact, and finally how certain we are about predictions of structures in the subsurface.
Related subject area
Solid Earth
Empirical modeling of tropospheric delays with uncertainty
CitcomSVE-3.0: a three-dimensional finite-element software package for modeling load-induced deformation and glacial isostatic adjustment for an Earth with a viscoelastic and compressible mantle
NSOAS24: a new global marine gravity model derived from multi-satellite sea surface slopes
GEOMAPLEARN 1.2: detecting structures from geological maps with machine learning – the case of geological folds
Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydromechanical problems with applications
Reconciling surface deflections from simulations of global mantle convection
A Bayesian framework for inferring regional and global change from stratigraphic proxy records (StratMC v1.0)
Three-dimensional analytical solution of self-potential from regularly polarized bodies in a layered seafloor model
A fast surrogate model for 3D Earth glacial isostatic adjustment using Tensorflow (v2.8.0) artificial neural networks
ShellSet v1.1.0 parallel dynamic neotectonic modelling: a case study using Earth5-049
FastIsostasy v1.0 – a regional, accelerated 2D glacial isostatic adjustment (GIA) model accounting for the lateral variability of the solid Earth
Broken Terrains v. 1.0: A supervised detection of fault-related lineaments on geological terrains
Automatic adjoint-based inversion schemes for geodynamics: reconstructing the evolution of Earth's mantle in space and time
Benchmarking the accuracy of higher-order particle methods in geodynamic models of transient flow
REHEATFUNQ (REgional HEAT-Flow Uncertainty and aNomaly Quantification) 2.0.1: a model for regional aggregate heat flow distributions and anomaly quantification
A new temperature–photoperiod coupled phenology module in LPJ-GUESS model v4.1: optimizing estimation of terrestrial carbon and water processes
A dynamic informed deep learning method for future estimation of laboratory stick-slip
High-precision 1′ × 1′ bathymetric model of Philippine Sea inversed from marine gravity anomalies
Deciphering past earthquakes from the probabilistic modeling of paleoseismic records – the Paleoseismic EArthquake CHronologies code (PEACH, version 1)
Modelling detrital cosmogenic nuclide concentrations during landscape evolution in Cidre v2.0
IMEX_SfloW2D v2: a depth-averaged numerical flow model for volcanic gas–particle flows over complex topographies and water
Simulation of a fully coupled 3D glacial isostatic adjustment – ice sheet model for the Antarctic ice sheet over a glacial cycle
AdaHRBF v1.0: gradient-adaptive Hermite–Birkhoff radial basis function interpolants for three-dimensional stratigraphic implicit modeling
PySubdiv 1.0: open-source geological modeling and reconstruction by non-manifold subdivision surfaces
Reconstructing tephra fall deposits via ensemble-based data assimilation techniques
ClinoformNet-1.0: stratigraphic forward modeling and deep learning for seismic clinoform delineation
Addressing challenges in uncertainty quantification: the case of geohazard assessments
DeepISMNet: three-dimensional implicit structural modeling with convolutional neural network
Towards automatic finite-element methods for geodynamics via Firedrake
MagmaFOAM-1.0: a modular framework for the simulation of magmatic systems
A global, spherical finite-element model for post-seismic deformation using Abaqus
SMAUG v1.0 – a user-friendly muon simulator for the imaging of geological objects in 3-D
CliffDelineaTool v1.2.0: an algorithm for identifying coastal cliff base and top positions
Capturing the interactions between ice sheets, sea level and the solid Earth on a range of timescales: a new “time window” algorithm
Spatial agents for geological surface modelling
RHEA v1.0: Enabling fully coupled simulations with hydro-geomechanical heterogeneity
Modelling of faults in LoopStructural 1.0
PALEOSTRIPv1.0 – a user-friendly 3D backtracking software to reconstruct paleo-bathymetries
LoopStructural 1.0: time-aware geological modelling
Sub3DNet1.0: a deep-learning model for regional-scale 3D subsurface structure mapping
Analytical solutions for mantle flow in cylindrical and spherical shells
Towards a model for structured mass movements: the OpenLISEM hazard model 2.0a
GO_3D_OBS: the multi-parameter benchmark geomodel for seismic imaging method assessment and next-generation 3D survey design (version 1.0)
PLUME-MoM-TSM 1.0.0: a volcanic column and umbrella cloud spreading model
HydrothermalFoam v1.0: a 3-D hydrothermal transport model for natural submarine hydrothermal systems
Synthetic seismicity distribution in Guerrero–Oaxaca subduction zone, Mexico, and its implications on the role of asperities in Gutenberg–Richter law
A new open-source viscoelastic solid earth deformation module implemented in Elmer (v8.4)
CobWeb 1.0: machine learning toolbox for tomographic imaging
pygeodyn 1.1.0: a Python package for geomagnetic data assimilation
IMEX_SfloW2D 1.0: a depth-averaged numerical flow model for pyroclastic avalanches
Jungang Wang, Junping Chen, and Yize Zhang
Geosci. Model Dev., 18, 1487–1504, https://doi.org/10.5194/gmd-18-1487-2025, https://doi.org/10.5194/gmd-18-1487-2025, 2025
Short summary
Short summary
The Global Navigation Satellite System (GNSS) is widely used for real-time monitoring and early warning of geohazards. Accurate modeling of tropospheric delays is critical to achieving high-precision GNSS solutions, and using external delays can improve real-time GNSS convergence times. Current empirical delay models only provide delay but not uncertainty. We propose a global empirical delay model with uncertainty and demonstrate its benefits in accelerating GNSS positioning convergence.
Tao Yuan, Shijie Zhong, and Geruo A
Geosci. Model Dev., 18, 1445–1461, https://doi.org/10.5194/gmd-18-1445-2025, https://doi.org/10.5194/gmd-18-1445-2025, 2025
Short summary
Short summary
Earth and other planets deform under various forces. Numerical modeling is critical in understanding the nature of various dynamic deformation processes. This article introduces a newly developed open-source package, CitcomSVE-3.0, which efficiently solves the viscoelastic deformation of planetary bodies. We present benchmark results against a semi-analytical code. With its accuracy and efficiency, CitcomSVE-3.0 could advance research in planetary and climatic sciences.
Shengjun Zhang, Xu Chen, Runsheng Zhou, and Yongjun Jia
Geosci. Model Dev., 18, 1221–1239, https://doi.org/10.5194/gmd-18-1221-2025, https://doi.org/10.5194/gmd-18-1221-2025, 2025
Short summary
Short summary
NSOAS24, a new global marine gravity model derived from multi-satellite altimetry missions, represents a significant advancement over its predecessor, NSOAS22. Through optimized processing procedures, NSOAS24 resolves previous issues and demonstrates improved accuracy. Compared to NSOAS22, NSOAS24 shows a reduction of approximately 0.7 mGal in standard deviation when validated against recent shipborne data. Notably, its accuracy now rivals internationally recognized models DTU21 and V32.1.
David Oakley, Christelle Loiselet, Thierry Coowar, Vincent Labbe, and Jean-Paul Callot
Geosci. Model Dev., 18, 939–960, https://doi.org/10.5194/gmd-18-939-2025, https://doi.org/10.5194/gmd-18-939-2025, 2025
Short summary
Short summary
In this work, we develop two automated workflows for identifying fold structures on geological maps using machine learning. In one method, we identify map patterns suggestive of folding based on pre-defined rules and apply a clustering algorithm to group those from the same fold together. In the other, we train a convolutional neural network to identify folds based on a set of training examples. We apply both methods to a set of synthetic maps and to real-world maps from two locations in France.
Yury Alkhimenkov and Yury Y. Podladchikov
Geosci. Model Dev., 18, 563–583, https://doi.org/10.5194/gmd-18-563-2025, https://doi.org/10.5194/gmd-18-563-2025, 2025
Short summary
Short summary
The accelerated pseudo-transient (APT) method is an efficient way to solve partial differential equations, particularly well-suited for parallel computing. This paper explores the APT method's effectiveness in solving elastic, viscoelastic, and hydromechanical problems, focusing on quasi-static conditions in 1D, 2D, and 3D. The study examines the best numerical settings for fast and accurate solutions. The paper shows how the APT method can handle complex problems in high-resolution models.
Conor P. B. O'Malley, Gareth G. Roberts, James Panton, Fred D. Richards, J. Huw Davies, Victoria M. Fernandes, and Sia Ghelichkhan
Geosci. Model Dev., 17, 9023–9049, https://doi.org/10.5194/gmd-17-9023-2024, https://doi.org/10.5194/gmd-17-9023-2024, 2024
Short summary
Short summary
We wish to understand how the history of flowing rock within Earth's interior impacts deflection of its surface. Observations exist to address this problem, and mathematics and different computing tools can be used to predict histories of flow. We explore how modeling choices impact calculated vertical deflections. The sensitivity of vertical motions at Earth's surface to deep flow is assessed, demonstrating how surface observations can enlighten flow histories.
Stacey Edmonsond and Blake Dyer
EGUsphere, https://doi.org/10.5194/egusphere-2024-2579, https://doi.org/10.5194/egusphere-2024-2579, 2024
Short summary
Short summary
The chemistry of sedimentary rocks is used to reconstruct past changes in Earth’s climate and biogeochemical cycles. Reconstructing global change requires merging stratigraphic proxy records from many locations, each of which may be incomplete, time-uncertain, and influenced by both global and local processes. StratMC uses Bayesian modeling to see through this complexity, building more accurate and testable reconstructions of global change from stratigraphic data.
Pengfei Zhang, Yi-an Cui, Jing Xie, Youjun Guo, Jianxin Liu, and Jieran Liu
Geosci. Model Dev., 17, 8521–8533, https://doi.org/10.5194/gmd-17-8521-2024, https://doi.org/10.5194/gmd-17-8521-2024, 2024
Short summary
Short summary
Our study proposes a three-dimensional analytical solution for the self-potential (SP) method to detect seafloor massive sulfide deposits. Using the mirror image method, we derived formulas for the SP generated by polarized bodies in layered media. We conducted experiments with a simulated environment and compared the results with our theoretical predictions. The findings confirmed a high degree of accuracy, demonstrating the reliability of our approach for faster and precise SP modeling.
Ryan Love, Glenn A. Milne, Parviz Ajourlou, Soran Parang, Lev Tarasov, and Konstantin Latychev
Geosci. Model Dev., 17, 8535–8551, https://doi.org/10.5194/gmd-17-8535-2024, https://doi.org/10.5194/gmd-17-8535-2024, 2024
Short summary
Short summary
A relatively recent advance in glacial isostatic adjustment modeling has been the development of models that include 3D Earth structure, as opposed to 1D structure. However, a major limitation is the computational expense. We have developed a method using artificial neural networks to emulate the influence of 3D Earth models to affordably constrain the viscosity parameter space. Our results indicate that the misfits are of a scale such that useful predictions of relative sea level can be made.
Jon B. May, Peter Bird, and Michele M. C. Carafa
Geosci. Model Dev., 17, 6153–6171, https://doi.org/10.5194/gmd-17-6153-2024, https://doi.org/10.5194/gmd-17-6153-2024, 2024
Short summary
Short summary
ShellSet is a combination of well-known geoscience software packages. It features a simple user interface and is optimised through the addition of a grid search input option (automatically searching for optimal models within a defined N-dimensional parameter space) and the ability to run multiple models in parallel. We show that for each number of models tested there is a performance benefit to parallel running, while two examples demonstrate a use case by improving an existing global model.
Jan Swierczek-Jereczek, Marisa Montoya, Konstantin Latychev, Alexander Robinson, Jorge Alvarez-Solas, and Jerry Mitrovica
Geosci. Model Dev., 17, 5263–5290, https://doi.org/10.5194/gmd-17-5263-2024, https://doi.org/10.5194/gmd-17-5263-2024, 2024
Short summary
Short summary
Ice sheets present a thickness of a few kilometres, leading to a vertical deformation of the crust of up to a kilometre. This process depends on properties of the solid Earth, which can be regionally very different. We propose a model that accounts for this often-ignored heterogeneity and run 100 000 simulation years in minutes. Thus, the evolution of ice sheets is modeled with better accuracy, which is critical for a good mitigation of climate change and, in particular, sea-level rise.
Michał Michalak, Christian Gerhards, and Peter Menzel
EGUsphere, https://doi.org/10.5194/egusphere-2024-2004, https://doi.org/10.5194/egusphere-2024-2004, 2024
Short summary
Short summary
This study presents a novel method for fault detection on geological terrains. Using synthetic models, we applied machine learning to classify terrain shape and nearby features. Testing on real borehole data validated its effectiveness across various fault orientations. The supervised approach represents a significant improvement over older methods that relied on simpler clustering techniques which were capable of identifying less orientations of faults.
Sia Ghelichkhan, Angus Gibson, D. Rhodri Davies, Stephan C. Kramer, and David A. Ham
Geosci. Model Dev., 17, 5057–5086, https://doi.org/10.5194/gmd-17-5057-2024, https://doi.org/10.5194/gmd-17-5057-2024, 2024
Short summary
Short summary
We introduce the Geoscientific ADjoint Optimisation PlaTform (G-ADOPT), designed for inverse modelling of Earth system processes, with an initial focus on mantle dynamics. G-ADOPT is built upon Firedrake, Dolfin-Adjoint and the Rapid Optimisation Library, which work together to optimise models using an adjoint method, aligning them with seismic and geologic datasets. We demonstrate G-ADOPT's ability to reconstruct mantle evolution and thus be a powerful tool in geosciences.
Rene Gassmöller, Juliane Dannberg, Wolfgang Bangerth, Elbridge Gerry Puckett, and Cedric Thieulot
Geosci. Model Dev., 17, 4115–4134, https://doi.org/10.5194/gmd-17-4115-2024, https://doi.org/10.5194/gmd-17-4115-2024, 2024
Short summary
Short summary
Numerical models that use simulated particles are a powerful tool for investigating flow in the interior of the Earth, but the accuracy of these models is not fully understood. Here we present two new benchmarks that allow measurement of model accuracy. We then document that better accuracy matters for applications like convection beneath an oceanic plate. Our benchmarks and methods are freely available to help the community develop better models.
Malte Jörn Ziebarth and Sebastian von Specht
Geosci. Model Dev., 17, 2783–2828, https://doi.org/10.5194/gmd-17-2783-2024, https://doi.org/10.5194/gmd-17-2783-2024, 2024
Short summary
Short summary
Thermal energy from Earth’s active interior constantly dissipates through Earth’s surface. This heat flow is not spatially uniform, and its exact pattern is hard to predict since it depends on crustal and mantle properties, both varying across scales. Our new model REHEATFUNQ addresses this difficulty by treating the fluctuations of heat flow within a region statistically. REHEATFUNQ estimates the regional distribution of heat flow and quantifies known structural signals therein.
Shouzhi Chen, Yongshuo H. Fu, Mingwei Li, Zitong Jia, Yishuo Cui, and Jing Tang
Geosci. Model Dev., 17, 2509–2523, https://doi.org/10.5194/gmd-17-2509-2024, https://doi.org/10.5194/gmd-17-2509-2024, 2024
Short summary
Short summary
It is still a challenge to achieve an accurate simulation of vegetation phenology in the dynamic global vegetation models (DGVMs). We implemented and coupled the spring and autumn phenology models into one of the DGVMs, LPJ-GUESS, and substantially improved the accuracy in capturing the start and end dates of growing seasons. Our study highlights the importance of getting accurate phenology estimations to reduce the uncertainties in plant distribution and terrestrial carbon and water cycling.
Enjiang Yue, Mengjiao Qin, Linshu Hu, Sensen Wu, and Zhenhong Du
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-46, https://doi.org/10.5194/gmd-2024-46, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
Laboratory earthquakes is important means to understand natural earthquakes while previous work focused on transient prediction, lacking future prediction capability. We propose a method and evaluate on data from lab experiments with different slip behaviors. It outperforms state-of-the-art methods in modeling slip moments, intervals and predictions beyond trained horizons especially for challenging slip scenarios, which is crucial for quasi-periodic geophysical process like seismicity.
Dechao An, Jinyun Guo, Xiaotao Chang, Zhenming Wang, Yongjun Jia, Xin Liu, Valery Bondur, and Heping Sun
Geosci. Model Dev., 17, 2039–2052, https://doi.org/10.5194/gmd-17-2039-2024, https://doi.org/10.5194/gmd-17-2039-2024, 2024
Short summary
Short summary
Seafloor topography, as fundamental geoinformation in marine surveying and mapping, plays a crucial role in numerous scientific studies. In this paper, we focus on constructing a high-precision seafloor topography and bathymetry model for the Philippine Sea (5° N–35° N, 120° E–150° E), based on shipborne bathymetric data and marine gravity anomalies, and evaluate the reliability of the model's accuracy.
Octavi Gómez-Novell, Bruno Pace, Francesco Visini, Joanna Faure Walker, and Oona Scotti
Geosci. Model Dev., 16, 7339–7355, https://doi.org/10.5194/gmd-16-7339-2023, https://doi.org/10.5194/gmd-16-7339-2023, 2023
Short summary
Short summary
Knowing the rate at which earthquakes happen along active faults is crucial to characterize the hazard that they pose. We present an approach (Paleoseismic EArthquake CHronologies, PEACH) to correlate and compute seismic histories using paleoseismic data, a type of data that characterizes past seismic activity from the geological record. Our approach reduces the uncertainties of the seismic histories and overall can improve the knowledge on fault rupture behavior for the seismic hazard.
Sébastien Carretier, Vincent Regard, Youssouf Abdelhafiz, and Bastien Plazolles
Geosci. Model Dev., 16, 6741–6755, https://doi.org/10.5194/gmd-16-6741-2023, https://doi.org/10.5194/gmd-16-6741-2023, 2023
Short summary
Short summary
We present the development of a code to simulate simultaneously the dynamics of landscapes over geological time and the evolution of the concentration of cosmogenic isotopes in grains throughout their transport from the slopes to the river outlets. This new model makes it possible to study the relationship between the detrital signal of cosmogenic isotope concentration measured in sediment and the erosion--deposition processes in watersheds.
Mattia de' Michieli Vitturi, Tomaso Esposti Ongaro, and Samantha Engwell
Geosci. Model Dev., 16, 6309–6336, https://doi.org/10.5194/gmd-16-6309-2023, https://doi.org/10.5194/gmd-16-6309-2023, 2023
Short summary
Short summary
We present version 2 of the numerical code IMEX-SfloW2D. With this version it is possible to simulate a wide range of volcanic mass flows (pyroclastic avalanches, lahars, pyroclastic surges), and here we present its application to transient dilute pyroclastic density currents (PDCs). A simulation of the 1883 Krakatau eruption demonstrates the capability of the numerical model to face a complex natural case involving the propagation of PDCs over the sea surface and across topographic obstacles.
Caroline J. van Calcar, Roderik S. W. van de Wal, Bas Blank, Bas de Boer, and Wouter van der Wal
Geosci. Model Dev., 16, 5473–5492, https://doi.org/10.5194/gmd-16-5473-2023, https://doi.org/10.5194/gmd-16-5473-2023, 2023
Short summary
Short summary
The waxing and waning of the Antarctic ice sheet caused the Earth’s surface to deform, which is stabilizing the ice sheet and mainly determined by the spatially variable viscosity of the mantle. Including this feedback in model simulations led to significant differences in ice sheet extent and ice thickness over the last glacial cycle. The results underline and quantify the importance of including this local feedback effect in ice sheet models when simulating the Antarctic ice sheet evolution.
Baoyi Zhang, Linze Du, Umair Khan, Yongqiang Tong, Lifang Wang, and Hao Deng
Geosci. Model Dev., 16, 3651–3674, https://doi.org/10.5194/gmd-16-3651-2023, https://doi.org/10.5194/gmd-16-3651-2023, 2023
Short summary
Short summary
We propose a Hermite–Birkhoff radial basis function (HRBF) formulation, AdaHRBF, with an adaptive gradient magnitude for continuous 3D stratigraphic potential field (SPF) modeling of multiple stratigraphic interfaces. In the linear system of HRBF interpolants constrained by the scattered on-contact attribute points and off-contact attitude points of a set of strata in 3D space, we add a novel optimization term to iteratively obtain the true gradient magnitude.
Mohammad Moulaeifard, Simon Bernard, and Florian Wellmann
Geosci. Model Dev., 16, 3565–3579, https://doi.org/10.5194/gmd-16-3565-2023, https://doi.org/10.5194/gmd-16-3565-2023, 2023
Short summary
Short summary
In this work, we propose a flexible framework to generate and interact with geological models using explicit surface representations. The essence of the work lies in the determination of the flexible control mesh, topologically similar to the main geological structure, watertight and controllable with few control points, to manage the geological structures. We exploited the subdivision surface method in our work, which is commonly used in the animation and gaming industry.
Leonardo Mingari, Antonio Costa, Giovanni Macedonio, and Arnau Folch
Geosci. Model Dev., 16, 3459–3478, https://doi.org/10.5194/gmd-16-3459-2023, https://doi.org/10.5194/gmd-16-3459-2023, 2023
Short summary
Short summary
Two novel techniques for ensemble-based data assimilation, suitable for semi-positive-definite variables with highly skewed uncertainty distributions such as tephra deposit mass loading, are applied to reconstruct the tephra fallout deposit resulting from the 2015 Calbuco eruption in Chile. The deposit spatial distribution and the ashfall volume according to the analyses are in good agreement with estimations based on field measurements and isopach maps reported in previous studies.
Hui Gao, Xinming Wu, Jinyu Zhang, Xiaoming Sun, and Zhengfa Bi
Geosci. Model Dev., 16, 2495–2513, https://doi.org/10.5194/gmd-16-2495-2023, https://doi.org/10.5194/gmd-16-2495-2023, 2023
Short summary
Short summary
We propose a workflow to automatically generate synthetic seismic data and corresponding stratigraphic labels (e.g., clinoform facies, relative geologic time, and synchronous horizons) by geological and geophysical forward modeling. Trained with only synthetic datasets, our network works well to accurately and efficiently predict clinoform facies in 2D and 3D field seismic data. Such a workflow can be easily extended for other geological and geophysical scenarios in the future.
Ibsen Chivata Cardenas, Terje Aven, and Roger Flage
Geosci. Model Dev., 16, 1601–1615, https://doi.org/10.5194/gmd-16-1601-2023, https://doi.org/10.5194/gmd-16-1601-2023, 2023
Short summary
Short summary
We discuss challenges in uncertainty quantification for geohazard assessments. The challenges arise from limited data and the one-off nature of geohazard features. The challenges include the credibility of predictions, input uncertainty, and assumptions’ impact. Considerations to increase credibility of the quantification are provided. Crucial tasks in the quantification are the exhaustive scrutiny of the background knowledge coupled with the assessment of deviations of assumptions made.
Zhengfa Bi, Xinming Wu, Zhaoliang Li, Dekuan Chang, and Xueshan Yong
Geosci. Model Dev., 15, 6841–6861, https://doi.org/10.5194/gmd-15-6841-2022, https://doi.org/10.5194/gmd-15-6841-2022, 2022
Short summary
Short summary
We present an implicit modeling method based on deep learning to produce a geologically valid and structurally compatible model from unevenly sampled structural data. Trained with automatically generated synthetic data with realistic features, our network can efficiently model geological structures without the need to solve large systems of mathematical equations, opening new opportunities for further leveraging deep learning to improve modeling capacity in many Earth science applications.
D. Rhodri Davies, Stephan C. Kramer, Sia Ghelichkhan, and Angus Gibson
Geosci. Model Dev., 15, 5127–5166, https://doi.org/10.5194/gmd-15-5127-2022, https://doi.org/10.5194/gmd-15-5127-2022, 2022
Short summary
Short summary
Firedrake is a state-of-the-art system that automatically generates highly optimised code for simulating finite-element (FE) problems in geophysical fluid dynamics. It creates a separation of concerns between employing the FE method and implementing it. Here, we demonstrate the applicability and benefits of Firedrake for simulating geodynamical flows, with a focus on the slow creeping motion of Earth's mantle over geological timescales, which is ultimately the engine driving our dynamic Earth.
Federico Brogi, Simone Colucci, Jacopo Matrone, Chiara Paola Montagna, Mattia De' Michieli Vitturi, and Paolo Papale
Geosci. Model Dev., 15, 3773–3796, https://doi.org/10.5194/gmd-15-3773-2022, https://doi.org/10.5194/gmd-15-3773-2022, 2022
Short summary
Short summary
Computer simulations play a fundamental role in understanding volcanic phenomena. The growing complexity of these simulations requires the development of flexible computational tools that can easily switch between sub-models and solution techniques as well as optimizations. MagmaFOAM is a newly developed library that allows for maximum flexibility for solving multiphase volcanic flows and promotes collaborative work for in-house and community model development, testing, and comparison.
Grace A. Nield, Matt A. King, Rebekka Steffen, and Bas Blank
Geosci. Model Dev., 15, 2489–2503, https://doi.org/10.5194/gmd-15-2489-2022, https://doi.org/10.5194/gmd-15-2489-2022, 2022
Short summary
Short summary
We present a finite-element model of post-seismic solid Earth deformation built in the software package Abaqus for the purpose of calculating post-seismic deformation in the far field of major earthquakes. The model is benchmarked against an existing open-source post-seismic model demonstrating good agreement. The advantage over existing models is the potential for simple modification to include 3-D Earth structure, non-linear rheologies and alternative or multiple sources of stress change.
Alessandro Lechmann, David Mair, Akitaka Ariga, Tomoko Ariga, Antonio Ereditato, Ryuichi Nishiyama, Ciro Pistillo, Paola Scampoli, Mykhailo Vladymyrov, and Fritz Schlunegger
Geosci. Model Dev., 15, 2441–2473, https://doi.org/10.5194/gmd-15-2441-2022, https://doi.org/10.5194/gmd-15-2441-2022, 2022
Short summary
Short summary
Muon tomography is a technology that is used often in geoscientific research. The know-how of data analysis is, however, still possessed by physicists who developed this technology. This article aims at providing geoscientists with the necessary tools to perform their own analyses. We hope that a lower threshold to enter the field of muon tomography will allow more geoscientists to engage with muon tomography. SMAUG is set up in a modular way to allow for its own modules to work in between.
Zuzanna M. Swirad and Adam P. Young
Geosci. Model Dev., 15, 1499–1512, https://doi.org/10.5194/gmd-15-1499-2022, https://doi.org/10.5194/gmd-15-1499-2022, 2022
Short summary
Short summary
Cliff base and top lines that delimit coastal cliff faces are usually manually digitized based on maps, aerial photographs, terrain models, etc. However, manual mapping is time consuming and depends on the mapper's decisions and skills. To increase the objectivity and efficiency of cliff mapping, we developed CliffDelineaTool, an algorithm that identifies cliff base and top positions along cross-shore transects using elevation and slope characteristics.
Holly Kyeore Han, Natalya Gomez, and Jeannette Xiu Wen Wan
Geosci. Model Dev., 15, 1355–1373, https://doi.org/10.5194/gmd-15-1355-2022, https://doi.org/10.5194/gmd-15-1355-2022, 2022
Short summary
Short summary
Interactions between ice sheets, sea level and the solid Earth occur over a range of timescales from years to tens of thousands of years. This requires coupled ice-sheet–sea-level models to exchange information frequently, leading to a quadratic increase in computation time with the number of model timesteps. We present a new sea-level model algorithm that allows coupled models to improve the computational feasibility and precisely capture short-term interactions within longer simulations.
Eric A. de Kemp
Geosci. Model Dev., 14, 6661–6680, https://doi.org/10.5194/gmd-14-6661-2021, https://doi.org/10.5194/gmd-14-6661-2021, 2021
Short summary
Short summary
This is a proof of concept and review paper of spatial agents, with initial research focusing on geomodelling. The results may be of interest to others working on complex regional geological modelling with sparse data. Structural agent-based swarming behaviour is key to advancing this field. The study provides groundwork for research in structural geology 3D modelling with spatial agents. This work was done with NetLogo, a free agent modelling platform used mostly for teaching complex systems.
José M. Bastías Espejo, Andy Wilkins, Gabriel C. Rau, and Philipp Blum
Geosci. Model Dev., 14, 6257–6272, https://doi.org/10.5194/gmd-14-6257-2021, https://doi.org/10.5194/gmd-14-6257-2021, 2021
Short summary
Short summary
The hydraulic and mechanical properties of the subsurface are inherently heterogeneous. RHEA is a simulator that can perform couple hydro-geomechanical processes in heterogeneous porous media with steep gradients. RHEA is able to fully integrate spatial heterogeneity, allowing allocation of distributed hydraulic and geomechanical properties at mesh element level. RHEA is a valuable tool that can simulate problems considering realistic heterogeneity inherent to geologic formations.
Lachlan Grose, Laurent Ailleres, Gautier Laurent, Guillaume Caumon, Mark Jessell, and Robin Armit
Geosci. Model Dev., 14, 6197–6213, https://doi.org/10.5194/gmd-14-6197-2021, https://doi.org/10.5194/gmd-14-6197-2021, 2021
Short summary
Short summary
Fault discontinuities in rock packages represent the plane where two blocks of rock have moved. They are challenging to incorporate into geological models because the geometry of the faulted rock units are defined by not only the location of the discontinuity but also the kinematics of the fault. In this paper, we outline a structural geology framework for incorporating faults into geological models by directly incorporating kinematics into the mathematical framework of the model.
Florence Colleoni, Laura De Santis, Enrico Pochini, Edy Forlin, Riccardo Geletti, Giuseppe Brancatelli, Magdala Tesauro, Martina Busetti, and Carla Braitenberg
Geosci. Model Dev., 14, 5285–5305, https://doi.org/10.5194/gmd-14-5285-2021, https://doi.org/10.5194/gmd-14-5285-2021, 2021
Short summary
Short summary
PALEOSTRIP has been developed in the framework of past Antarctic ice sheet reconstructions for periods when bathymetry around Antarctica differed substantially from today. It has been designed for users with no knowledge of numerical modelling and allows users to switch on and off the processes involved in backtracking and backstripping. Applications are broad, and it can be used to restore any continental margin bathymetry or sediment thickness and to perform basin analysis.
Lachlan Grose, Laurent Ailleres, Gautier Laurent, and Mark Jessell
Geosci. Model Dev., 14, 3915–3937, https://doi.org/10.5194/gmd-14-3915-2021, https://doi.org/10.5194/gmd-14-3915-2021, 2021
Short summary
Short summary
LoopStructural is an open-source 3D geological modelling library with a model design allowing for multiple different algorithms to be used for comparison for the same geology. Geological structures are modelled using structural geology concepts and techniques, allowing for complex structures such as overprinted folds and faults to be modelled. In the paper, we demonstrate automatically generating a 3-D model from map2loop-processed geological survey data of the Flinders Ranges, South Australia.
Zhenjiao Jiang, Dirk Mallants, Lei Gao, Tim Munday, Gregoire Mariethoz, and Luk Peeters
Geosci. Model Dev., 14, 3421–3435, https://doi.org/10.5194/gmd-14-3421-2021, https://doi.org/10.5194/gmd-14-3421-2021, 2021
Short summary
Short summary
Fast and reliable tools are required to extract hidden information from big geophysical and remote sensing data. A deep-learning model in 3D image construction from 2D image(s) is here developed for paleovalley mapping from globally available digital elevation data. The outstanding performance for 3D subsurface imaging gives confidence that this generic novel tool will make better use of existing geophysical and remote sensing data for improved management of limited earth resources.
Stephan C. Kramer, D. Rhodri Davies, and Cian R. Wilson
Geosci. Model Dev., 14, 1899–1919, https://doi.org/10.5194/gmd-14-1899-2021, https://doi.org/10.5194/gmd-14-1899-2021, 2021
Short summary
Short summary
Computational models of Earth's mantle require rigorous verification and validation. Analytical solutions of the underlying Stokes equations provide a method to verify that these equations are accurately solved for. However, their derivation in spherical and cylindrical shell domains with physically relevant boundary conditions is involved. This paper provides a number of solutions. They are provided in a Python package (Assess) and their use is demonstrated in a convergence study with Fluidity.
Bastian van den Bout, Theo van Asch, Wei Hu, Chenxiao X. Tang, Olga Mavrouli, Victor G. Jetten, and Cees J. van Westen
Geosci. Model Dev., 14, 1841–1864, https://doi.org/10.5194/gmd-14-1841-2021, https://doi.org/10.5194/gmd-14-1841-2021, 2021
Short summary
Short summary
Landslides, debris flows and other types of dense gravity-driven flows threaten livelihoods around the globe. Understanding the mechanics of these flows can be crucial for predicting their behaviour and reducing disaster risk. Numerical models assume that the solids and fluids of the flow are unstructured. The newly presented model captures the internal structure during movement. This important step can lead to more accurate predictions of landslide movement.
Andrzej Górszczyk and Stéphane Operto
Geosci. Model Dev., 14, 1773–1799, https://doi.org/10.5194/gmd-14-1773-2021, https://doi.org/10.5194/gmd-14-1773-2021, 2021
Short summary
Short summary
We present the 3D multi-parameter synthetic geomodel of the subduction zone, as well as the workflow designed to implement all of its components. The model contains different geological structures of various scales and complexities. It is intended to serve as a tool for the geophysical community to validate imaging approaches, design acquisition techniques, estimate uncertainties, benchmark computing approaches, etc.
Mattia de' Michieli Vitturi and Federica Pardini
Geosci. Model Dev., 14, 1345–1377, https://doi.org/10.5194/gmd-14-1345-2021, https://doi.org/10.5194/gmd-14-1345-2021, 2021
Short summary
Short summary
Here, we present PLUME-MoM-TSM, a volcanic plume model that allows us to quantify the formation of aggregates during the rise of the plume, model the phase change of water, and include the possibility to simulate the initial spreading of the tephra umbrella cloud intruding from the volcanic column into the atmosphere. The model is first applied to the 2015 Calbuco eruption (Chile) and provides an analytical relationship between the upwind spreading and some characteristic of the volcanic column.
Zhikui Guo, Lars Rüpke, and Chunhui Tao
Geosci. Model Dev., 13, 6547–6565, https://doi.org/10.5194/gmd-13-6547-2020, https://doi.org/10.5194/gmd-13-6547-2020, 2020
Short summary
Short summary
We present the 3-D hydro-thermo-transport model HydrothermalFoam v1.0, which we designed to provide the marine geosciences community with an easy-to-use and state-of-the-art tool for simulating mass and energy transport in submarine hydrothermal systems. HydrothermalFoam is based on the popular open-source platform OpenFOAM, comes with a number of tutorials, and is published under the GNU General Public License v3.0.
Marisol Monterrubio-Velasco, F. Ramón Zúñiga, Quetzalcoatl Rodríguez-Pérez, Otilio Rojas, Armando Aguilar-Meléndez, and Josep de la Puente
Geosci. Model Dev., 13, 6361–6381, https://doi.org/10.5194/gmd-13-6361-2020, https://doi.org/10.5194/gmd-13-6361-2020, 2020
Short summary
Short summary
The Mexican subduction zone along the Pacific coast is one of the most active seismic zones in the world, where every year larger-magnitude earthquakes shake huge inland cities such as Mexico City. In this work, we use TREMOL (sThochastic Rupture Earthquake ModeL) to simulate the seismicity observed in this zone. Our numerical results reinforce the hypothesis that in some subduction regions single asperities are responsible for producing the observed seismicity.
Thomas Zwinger, Grace A. Nield, Juha Ruokolainen, and Matt A. King
Geosci. Model Dev., 13, 1155–1164, https://doi.org/10.5194/gmd-13-1155-2020, https://doi.org/10.5194/gmd-13-1155-2020, 2020
Short summary
Short summary
We present a newly developed flat-earth model, Elmer/Earth, for viscoelastic treatment of solid earth deformation under ice loads. Unlike many previous approaches with proprietary software, this model is based on the open-source FEM code Elmer, with the advantage for scientists to apply and alter the model without license constraints. The new-generation full-stress ice-sheet model Elmer/Ice shares the same code base, enabling future coupled ice-sheet–glacial-isostatic-adjustment simulations.
Swarup Chauhan, Kathleen Sell, Wolfram Rühaak, Thorsten Wille, and Ingo Sass
Geosci. Model Dev., 13, 315–334, https://doi.org/10.5194/gmd-13-315-2020, https://doi.org/10.5194/gmd-13-315-2020, 2020
Short summary
Short summary
We present CobWeb 1.0, a graphical user interface for analysing tomographic images of geomaterials. CobWeb offers different machine learning techniques for accurate multiphase image segmentation and visualizing material specific parameters such as pore size distribution, relative porosity and volume fraction. We demonstrate a novel approach of dual filtration and dual segmentation to eliminate edge enhancement artefact in synchrotron-tomographic datasets and provide the computational code.
Loïc Huder, Nicolas Gillet, and Franck Thollard
Geosci. Model Dev., 12, 3795–3803, https://doi.org/10.5194/gmd-12-3795-2019, https://doi.org/10.5194/gmd-12-3795-2019, 2019
Short summary
Short summary
The pygeodyn package is a geomagnetic data assimilation tool written in Python. It gives access to the Earth's core flow dynamics, controlled by geomagnetic observations, by means of a reduced numerical model anchored to geodynamo simulation statistics. It aims to provide the community with a user-friendly and tunable data assimilation algorithm. It can be used for education, geomagnetic model production or tests in conjunction with webgeodyn, a set of visualization tools for geomagnetic models.
Mattia de' Michieli Vitturi, Tomaso Esposti Ongaro, Giacomo Lari, and Alvaro Aravena
Geosci. Model Dev., 12, 581–595, https://doi.org/10.5194/gmd-12-581-2019, https://doi.org/10.5194/gmd-12-581-2019, 2019
Short summary
Short summary
Pyroclastic avalanches are a type of granular flow generated at active volcanoes by different mechanisms, including the collapse of steep pyroclastic deposits (e.g., scoria and ash cones) and fountaining during moderately explosive eruptions. We present IMEX_SfloW2D, a depth-averaged flow model describing the granular mixture as a single-phase granular fluid. Benchmark cases and preliminary application to the simulation of the 11 February pyroclastic avalanche at Mt. Etna (Italy) are shown.
Cited articles
Ailleres, L., Jessell, M., de Kemp, E., Caumon, G., Wellmann, F., Grose, L., Armit, R., Lindsay, M., Giraud, J., Brodaric, B., Harrison, M., and Courrioux, G.: Loop – Enabling 3D stochastic geological modelling, ASEG Ext. Abstr., 2019, 1–3, https://doi.org/10.1080/22020586.2019.12072955, 2019.
Belge, M., Kilmer, M. E., and Miller, E. L.: Efficient determination of multiple regularization parameters in a generalized L-curve framework, Inverse Probl., 18, 314, https://doi.org/10.1088/0266-5611/18/4/314, 2002.
Bertsekas, D. P.: Nonlinear programming, 3rd edn., Athena Scientific, Belmont, Massachussetts, USA, ISBN: 978-1-886529-05-2, 2016.
Bhattacharyya, B. K.: Magnetic anomalies due to prism-shaped bodies with arbitrary polarization, Geophysics, 29, 517–531, 1964.
Bhattacharyya, B. K.: A generalized multibody model for inversion of magnetic anomalies, Geophysics, 45, 255–270, https://doi.org/10.1190/1.1441081, 1980.
Bijani, R., Lelièvre, P. G., Ponte-Neto, C. F., and Farquharson, C. G.: Physical-property-, lithology- and surface-geometry-based joint inversion using Pareto Multi-Objective Global Optimization, Geophys. J. Int., 209, 730–748, https://doi.org/10.1093/gji/ggx046, 2017.
Blakely, R. J.: Potential Theory in Gravity and Magnetic Applications, Cambridge University Press, Cambridge, 1995.
Boyd, S.: Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers, Found. Trends® Mach. Learn., 3, 1–122, https://doi.org/10.1561/2200000016, 2010.
Boyd, S., Parikh, N., Chu, E., Peleato, B., and Eckstein, J.: Distributedoptimization and statistical learning via the alternating direction methodof multipliers, Foundations and Trends® in Machine Learning, 3, 1–122, https://doi.org/10.1561/2200000016, 2011.
Brown, V., Key, K., and Singh, S.: Seismically regularized controlled-source electromagnetic inversion, Geophysics, 77, E57–E65, https://doi.org/10.1190/geo2011-0081.1, 2012.
Calcagno, P., Chilès, J. P., Courrioux, G., and Guillen, A.: Geological modelling from field data and geological knowledge, Part I. Modelling method coupling 3D potential-field interpolation and geological rules, Phys. Earth Planet. In., 171, 147–157, https://doi.org/10.1016/j.pepi.2008.06.013, 2008.
Carter-McAuslan, A., Lelièvre, P. G., and Farquharson, C. G.: A study of fuzzy c -means coupling for joint inversion, using seismic tomography and gravity data test scenarios, Geophysics, 80, W1–W15, https://doi.org/10.1190/geo2014-0056.1, 2015.
Craven, P. and Wahba, G.: Smoothing noisy data with spline functions, Numer. Math., 31, 377–403, https://doi.org/10.1007/BF01404567, 1978.
Čuma, M. and Zhdanov, M. S.: Massively parallel regularized 3D inversion of potential fields on CPUs and GPUs, Comput. Geosci., 62, 80–87, https://doi.org/10.1016/j.cageo.2013.10.004, 2014.
Čuma, M., Wilson, G. A., and Zhdanov, M. S.: Large-scale 3D inversion of potential field data, Geophys. Prospect., 60, 1186–1199, https://doi.org/10.1111/j.1365-2478.2011.01052.x, 2012.
de la Varga, M., Schaaf, A., and Wellmann, F.: GemPy 1.0: open-source stochastic geological modeling and inversion, Geosci. Model Dev., 12, 1–32, https://doi.org/10.5194/gmd-12-1-2019, 2019.
De Stefano, M., Golfré Andreasi, F., Re, S., Virgilio, M., and Snyder, F. F.: Multiple-domain, simultaneous joint inversion of geophysical data with application to subsalt imaging, Geophysics, 76, R69–R80, https://doi.org/10.1190/1.3554652, 2011.
Deal, M. M. and Nolet, G.: Nullspace shuttles, Geophys. J. Int., 124, 372–380, https://doi.org/10.1111/j.1365-246X.1996.tb07027.x, 1996.
Doetsch, J., Linde, N., Coscia, I., Greenhalgh, S. A., and Green, A. G.: Zonation for 3D aquifer characterization based on joint inversions of multimethod crosshole geophysical data, Geophysics, 75, G53–G64, https://doi.org/10.1190/1.3496476, 2010.
Dykstra, R. L.: An Algorithm for Restricted Least Squares Regression, J. Am. Stat. Assoc., 78, 837–842, https://doi.org/10.1080/01621459.1983.10477029, 1983.
Eckhardt, E. A.: Partnership between geology and geophysics in prospecting for oil, Geophysics, 5, 209–214, https://doi.org/10.1190/1.1441804, 1940.
Farquharson, C. G. and Oldenburg, D. W.: A comparison of automatic techniques for estimating the regularization parameter in non-linear inverse problems, Geophys. J. Int., 156, 411–425, https://doi.org/10.1111/j.1365-246X.2004.02190.x, 2004.
Fichtner, A. and Zunino, A.: Hamiltonian Nullspace Shuttles, Geophys. Res. Lett., 46, 644–651, https://doi.org/10.1029/2018GL080931, 2019.
Gallardo, L. A. and Meju, M. A.: Characterization of
heterogeneous near-surface materials by joint 2D inversion of dc
resistivity and seismic data, Geophys. Res. Lett., 30,
1–4, https://doi.org/10.1029/2003GL017370, 2003.
Gallardo, L. A. and Meju, M. A.: Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints, J. Geophys. Res.-Sol. Ea., 109, B03311, https://doi.org/10.1029/2003JB002716, 2004.
Gil, Y., David, C. H., Demir, I., Essawy, B. T., Fulweiler, R. W., Goodall, J. L., Karlstrom, L., Lee, H., Mills, H. J., Oh, J.-H., Pierce, S. A., Pope, A., Tzeng, M. W., Villamizar, S. R., and Yu, X.: Toward the Geoscience Paper of the Future: Best practices for documenting and sharing research from data to software to provenance, Earth Space Sci., 3, 388–415, https://doi.org/10.1002/2015EA000136, 2016.
Giraud, J., Jessell, M., Lindsay, M., Parkyuz-Charrier, E., and Martin, R.: Integrated geophysical joint inversion using petrophysical constraints and geological modelling, in: SEG Technical Program Expanded Abstracts 2016, pp. 1597–1601, Society of Exploration Geophysicists, Dallas, Texas, 2016.
Giraud, J., Pakyuz-Charrier, E., Jessell, M., Lindsay, M., Martin, R., and Ogarko, V.: Uncertainty reduction through geologically conditioned petrophysical constraints in joint inversion, Geophysics, 82, ID19–ID34, https://doi.org/10.1190/geo2016-0615.1, 2017.
Giraud, J., Pakyuz-Charrier, E., Ogarko, V., Jessell, M., Lindsay, M., and Martin, R.: Impact of uncertain geology in constrained geophysical inversion, ASEG Ext. Abstr., 2018, 1, https://doi.org/10.1071/ASEG2018abM1_2F, 2018.
Giraud, J., Lindsay, M., Ogarko, V., Jessell, M., Martin, R., and Pakyuz-Charrier, E.: Integration of geoscientific uncertainty into geophysical inversion by means of local gradient regularization, Solid Earth, 10, 193–210, https://doi.org/10.5194/se-10-193-2019, 2019a.
Giraud, J., Ogarko, V., Lindsay, M., Pakyuz-Charrier, E., Jessell, M., and Martin, R.: Sensitivity of constrained joint inversions to geological and petrophysical input data uncertainties with posterior geological analysis, Geophys. J. Int., 218, 666–688, https://doi.org/10.1093/gji/ggz152, 2019b.
Giraud, J., Ogarko, V., Lindsay, M., Pakyuz-Charrier, E., Jessell, M., and Martin, R.: Sensitivity of constrained joint inversions to geological and petrophysical input data uncertainties with posterior geological analysis, Geophys. J. Int., 218, 666–688, https://doi.org/10.1093/gji/ggz152, 2019c.
Giraud, J., Lindsay, M., Ogarko, V., Jessell, M., Martin, R., and Pakyuz-Charrier, E.: Integration of geoscientific uncertainty into geophysical inversion by means of local gradient regularization, Solid Earth, 10, 193–210, https://doi.org/10.5194/se-10-193-2019, 2019d.
Giraud, J., Lindsay, M., Jessell, M., and Ogarko, V.: Towards plausible lithological classification from geophysical inversion: honouring geological principles in subsurface imaging, Solid Earth, 11, 419–436, https://doi.org/10.5194/se-11-419-2020, 2020.
Guillen, A., Calcagno, P., Courrioux, G., Joly, A., and Ledru, P.: Geological modelling from field data and geological knowledge, Part II. Modelling validation using gravity and magnetic data inversion, Phys. Earth Planet. In., 171, 158–169, https://doi.org/10.1016/j.pepi.2008.06.014, 2008.
Hammond, G. E. and Lichtner, P. C.: Field-scale model for
the natural attenuation of uranium at the Hanford 300 Area using
high-performance computing, Water Resour. Res., 46, W09527, https://doi.org/10.1029/2009WR008819, 2010.
Hansen, P. C. and Johnston, P. R.: The L-Curve and its Use in the Numerical Treatment of Inverse Problems, in: Computational Inverse Problems in Electrocardiography, pp. 119–142, available at: https://www.sintef.no/globalassets/project/evitameeting/2005/lcurve.pdf (last access: 13 September 2021), 2001.
Hansen, P. C. and O'Leary, D. P.: The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems, SIAM J. Sci. Comput., 14, 1487–1503, https://doi.org/10.1137/0914086, 1993.
Heincke, B., Jegen, M., Moorkamp, M., Hobbs, R. W., and Chen, J.: An adaptive coupling strategy for joint inversions that use petrophysical information as constraints, J. Appl. Geophys., 136, 279–297, https://doi.org/10.1016/j.jappgeo.2016.10.028, 2017.
Hoerl, A. E. and Kennard, R. W.: Ridge Regression: Application to nonorthogonal problems, Technometrics, 12, 69–82, https://doi.org/10.1080/00401706.1970.10488634, 1970.
Jessell, M., Giraud, J., and Lindsay, M.: 3-D geological and petrophysical models with synthetic geophysics based on data from the Hamersley region (Western Australia), Zenodo [data set], https://doi.org/10.5281/zenodo.4431796, 2021a.
Jessell, M., Ogarko, V., de Rose, Y., Lindsay, M., Joshi, R., Piechocka, A., Grose, L., de la Varga, M., Ailleres, L., and Pirot, G.: Automated geological map deconstruction for 3D model construction using map2loop 1.0 and map2model 1.0, Geosci. Model Dev., 14, 5063–5092, https://doi.org/10.5194/gmd-14-5063-2021, 2021b.
Kamm, J., Lundin, I. A., Bastani, M., Sadeghi, M., and Pedersen, L. B.: Joint inversion of gravity, magnetic, and petrophysical data – A case study from a gabbro intrusion in Boden, Sweden, Geophysics, 80, B131–B152, https://doi.org/10.1190/geo2014-0122.1, 2015.
Komatitsch, D. and Martin, R.: An unsplit convolutional
perfectly matched layer improved at grazing incidence for the
seismic wave equation, Geophysics, 72, SM155-SM167, https://doi.org/10.1190/1.2757586, 2007.
Kostina, E., Saunders, M. A., and Schierle, I.:
Computation of Covariance Matrices for Constrainted Parameter
Estimation Problems using LSQR, Stanford Optim. Lab Rep.,
4026, 1–11, 2009.
Kumar, V., Grama, A., Gupta, A., and Karypis, G.: Introduction to Parallel Computing: Design and Analysis of Algorithms, Benjamin-Cummings Publishing Co., Inc., Redwood City, CA, 1994.
Lelièvre, P., Farquharson, C., and Hurich, C.: Joint inversion of seismic traveltimes and gravity data on unstructured grids with application to mineral exploration, Geophysics, 77, K1–K15, https://doi.org/10.1190/geo2011-0154.1, 2012.
Lelièvre, P. G. and Farquharson, C. G.: Integrated Imaging for Mineral Exploration, in: Integrated Imaging of the Earth: Theory and Applications, American Geophysical Union, 137–166, 2016.
Li, X. and Chouteau, M.: Three Dimensional Modeling in All Space, Surv. Geophys., 19, 339–368, https://doi.org/10.1023/A:1006554408567, 1998.
Li, Y. and Oldenburg, D. W.: 3-{D} inversion of magnetic data, Geophysics, 61, 394–408, 1996.
Li, Y. and Oldenburg, D. W.: 3-D inversion of gravity data, Geophysics, 63, 109–119, https://doi.org/10.1190/1.1444302, 1998.
Lindsay, M. D., Aillères, L., Jessell, M. W., de Kemp, E. A., and Betts, P. G.: Locating and quantifying geological uncertainty in three-dimensional models: Analysis of the Gippsland Basin, southeastern Australia, Tectonophysics, 546–547, 10–27, https://doi.org/10.1016/j.tecto.2012.04.007, 2012.
Lindsay, M. D., Perrouty, S., Jessell, M. W., and Aillères, L.: Making the link between geological and geophysical uncertainty: geodiversity in the Ashanti Greenstone Belt, Geophys. J. Int., 195, 903–922, https://doi.org/10.1093/gji/ggt311, 2013.
Martin, R., Komatitsch, D., Gedney, S. D., and Bruthiaux,
E.: A high-order time and space formulation of the unsplit perfectly
matched layer for the seismicwave equation using auxiliary
differential equations (ADE-PML), CMES-Comp. Model. Eng., 56, 17–42, https://doi.org/10.3970/cmes.2010.056.017, 2010.
Martin, R., Monteiller, V., Komatitsch, D., Perrouty, S., Jessell, M., Bonvalot, S., and Lindsay, M. D.: Gravity inversion using wavelet-based compression on parallel hybrid CPU/GPU systems: application to southwest Ghana, Geophys. J. Int., 195, 1594–1619, https://doi.org/10.1093/gji/ggt334, 2013.
Martin, R., Ogarko, V., Komatitsch, D., and Jessell, M.: Parallel three-dimensional electrical capacitance data imaging using a nonlinear inversion algorithm and Lp norm-based model, Measurement, 128, 428–445, https://doi.org/10.1016/j.measurement.2018.05.099, 2018.
Martin, R., Bodet, L., Tournat, V., and Rejiba, F.: Seismic wave propagation in nonlinear viscoelastic media using the auxiliary differential equation method, Geophys. J. Int., 216, 453–469, https://doi.org/10.1093/gji/ggy441, 2019.
Martin, R., Giraud, J., Ogarko, V., Chevrot, S., Beller, S., Gégout, P., and Jessell, M.: Three-dimensional gravity anomaly inversion in the Pyrenees using compressional seismic velocity model as structural similarity constraints, Geophys. J. Int., 225, 1063–1085, https://doi.org/10.1093/gji/ggaa414, 2020.
McLachlan, G. and Peel, D.: Finite Mixture Models, Wiley
Ser. Prob. Stat., 419,
https://doi.org/10.1198/tech.2002.s651, 2000.
Meju, M. A. and Gallardo, L. A.: Structural Coupling Approaches in Integrated Geophysical Imaging, pp. 49–67, 2016.
Moorkamp, M., Heincke, B., Jegen, M., Hobbs, R. W., and
Roberts, A. W.: Joint Inversion in Hydrocarbon Exploration, in:
Integrated Imaging of the Earth: Theory and Applications,
pp. 167–189, 2016.
Muñoz, G. and Rath, V.: Beyond smooth inversion: the use of nullspace projection for the exploration of non-uniqueness in MT, Geophys. J. Int., 164, 301–311, https://doi.org/10.1111/j.1365-246X.2005.02825.x, 2006.
Nettleton, L. L.: Geophysics, geology and oil finding, Geophysics, 14, 273–289, https://doi.org/10.1190/1.1437535, 1949.
Ogarko, V., Giraud, J., Martin, R., and Jessell, M.: Disjoint interval bound constraints using the alternating direction method of multipliers for geologically constrained inversion: Application to gravity data, Geophysics, 86, G1–G11, https://doi.org/10.1190/geo2019-0633.1, 2021a.
Ogarko, V., Giraud, J., and Martin, R.: Tomofast-x v1.0
source code, Zenodo [code], https://doi.org/10.5281/zenodo.4454220, 2021b.
Paasche, H. and Tronicke, J.: Cooperative inversion of 2D geophysical data sets: A zonal approach based on fuzzy c -means cluster analysis, Geophysics, 72, A35–A39, https://doi.org/10.1190/1.2670341, 2007.
Paige, C. C. and Saunders, M. A.: LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares, ACM T. Math. Software, 8, 43–71, https://doi.org/10.1145/355984.355989, 1982.
Pakyuz-Charrier, E., Lindsay, M., Ogarko, V., Giraud, J., and Jessell, M.: Monte Carlo simulation for uncertainty estimation on structural data in implicit 3-D geological modeling, a guide for disturbance distribution selection and parameterization, Solid Earth, 9, 385–402, https://doi.org/10.5194/se-9-385-2018, 2018.
Pakyuz-Charrier, E., Jessell, M., Giraud, J., Lindsay, M., and Ogarko, V.: Topological analysis in Monte Carlo simulation for uncertainty propagation, Solid Earth, 10, 1663–1684, https://doi.org/10.5194/se-10-1663-2019, 2019.
Parsekian, A. D., Singha, K., Minsley, B. J., Holbrook, W. S., and Slater, L.: Multiscale geophysical imaging of the critical zone, Rev. Geophys., 53, 1–26, https://doi.org/10.1002/2014RG000465, 2015.
Peng, R. D.: Reproducible Research in Computational Science, Science, 334, 1226–1227, https://doi.org/10.1126/science.1213847, 2011.
Portniaguine, O. and Zhdanov, M. S.: 3-D magnetic inversion with data compression and image focusing, Geophysics, 67, 1532–1541, https://doi.org/10.1190/1.1512749, 2002.
Rashidifard, M., Giraud, J., Ogarko, V., Jessell, M., and Lindsay, M.: Cooperative Inversion of Seismic and Gravity Data Using Weighted Structure-Based Constraints, in: NSG2020 3rd Conference on Geophysics for Mineral Exploration and Mining, European Association of Geoscientists & Engineers, online, pp. 1–5, 2020.
Ren, Z. and Kalscheuer, T.: Uncertainty and Resolution Analysis of 2D and 3D Inversion Models Computed from Geophysical Electromagnetic Data, Surv. Geophys., 41, 47–112, https://doi.org/10.1007/s10712-019-09567-3, 2019.
Rudin, L. I., Osher, S., and Fatemi, E.: Nonlinear total variation based noise removal algorithms, Physica D, 60, 259–268, https://doi.org/10.1016/0167-2789(92)90242-F, 1992.
Santos, E. T. F. and Bassrei, A.: L- and Θ-curve approaches for the selection of regularization parameter in geophysical diffraction tomography, Comput. Geosci., 33, 618–629, https://doi.org/10.1016/j.cageo.2006.08.013, 2007.
Scheidt, C., Li, L., and Caers, J.: Quantifying Uncertainty in Subsurface Systems, John Wiley & Sons, Inc., Hoboken, NJ, 2018.
Sun, J. and Li, Y.: Geophysical inversion using
petrophysical constraints with application to lithology
differentiation, SEG San Antonio 2011 Annual Meeting, 4–9 November 2012,
2644–2648, 2011.
Sun, J. and Li, Y.: Multidomain petrophysically constrained inversion and geology differentiation using guided fuzzy c-means clustering, Geophysics, 80, ID1–ID18, https://doi.org/10.1190/geo2014-0049.1, 2015.
Sun, J. and Li, Y.: Joint inversion of multiple geophysical data using guided fuzzy c -means clustering, Geophysics, 81, ID37–ID57, https://doi.org/10.1190/geo2015-0457.1, 2016.
Sun, J. and Li, Y.: Joint inversion of multiple geophysical and petrophysical data using generalized fuzzy clustering algorithms, Geophys. J. Int., 208, 1201–1216, https://doi.org/10.1093/gji/ggw442, 2017.
Sun, J., Li, Y., and Studies, M.: Joint inversion of multiple geophysical data: A petrophysical approach using guided fuzzy c-means clustering, SEG Las Vegas 2012 Annual Meeting, 18–23 September 2011, pp. 1–5, 2012.
Tarantola, A.: Inverse Problem Theory and Methods for Model Parameter Estimation, Society for Industrial and Applied Mathematics, Philadelphia, 2005.
Tarantola, A. and Valette, B.: Inverse Problems = Quest for Information, J. Geophys., 50, 159–170, https://doi.org/10.1038/nrn1011, 1982.
Tikhonov, A. N. and Arsenin, V. Y.: Solutions of Ill-Posed Problems, John Wiley, New York, 1977.
Tikhonov, A. N. and Arsenin, V. Y.: Solution of Ill-Posed Problems, Math. Comput., 32, 491, https://doi.org/10.2307/2006360, 1978.
Vasco, D. W.: Invariance, groups, and non-uniqueness: The discrete case, Geophys. J. Int., 168, 473–490, https://doi.org/10.1111/j.1365-246X.2006.03161.x, 2007.
Wegener, A.: Die Entstehung der Kontinente und Ozeane, Geogr. J., 61, 304, https://doi.org/10.2307/1781265, 1923.
Wellmann, F. and Caumon, G.: 3-D Structural geological
models: Concepts, methods, and uncertainties, C. Schmelzback (Ed.), pp. 1–121, Elsevier, Cambridge, Massachusetts, 2018.
Wiik, T., Nordskag, J. I., Dischler, E. Ø., and Nguyen, A. K.: Inversion of inline and broadside marine controlled-source electromagnetic data with constraints derived from seismic data, Geophys. Prospect., 63, 1371–1382, https://doi.org/10.1111/1365-2478.12294, 2015.
Yan, P., Kalscheuer, T., Hedin, P., and Garcia Juanatey, M. A.: Two-dimensional magnetotelluric inversion using reflection seismic data as constraints and application in the COSC project, Geophys. Res. Lett., 44, 3554–3563, https://doi.org/10.1002/2017GL072953, 2017.
Zhang, J. and Revil, A.: 2D joint inversion of geophysical data using petrophysical clustering and facies deformation, Geophysics, 80, M69–M88, 2015.
Zhdanov, M. S. and Lin, W.: Adaptive multinary inversion of gravity and gravity gradiometry data, Geophysics, 82, G101–G114, https://doi.org/10.1190/geo2016-0451.1, 2017.
Zhdanov, M. S., Gribenko, A., and Wilson, G.: Generalized joint inversion of multimodal geophysical data using Gramian constraints, Geophys. Res. Lett., 39, L09301, https://doi.org/10.1029/2012GL051233, 2012.
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 and provide theoretical details of the open-source algorithm we use.
We review different techniques to model the Earth's subsurface from geophysical data (gravity...
