Stochastic Inversion of Transient Electromagnetic Data to Derive Aquifer Geometry and associated Uncertainties

Context

This repository contains all the data, scripts, and libraries needed to reproduce the results (except for the QGIS map) from the article "Stochastic Inversion of Transient Electromagnetic Data to Derive Aquifer Geometry and Associated Uncertainties" by Lukas Aigner, Hadrien Michel, Thomas Hermans, and Adrián Flores Orozco.

Reuse instructions

Functions and classes needed to run the scripts are located in the modules subfolder. Please find the scripts to reproduce the figures of the manuscript within the project folder. The scripts are configured to add a relative import path to os.sys

To build the Python environments:
There are several Conda environments (*.yml) provided; the correct environment to use depends on the use case (i.e. which parts of the provided code to execute).
For everything related to pyBEL1D please use the bel1d environment, for the sensitivity analysis please use the salib_tem.yml file, while all the other scripts can be run using the empypg.yml.

Special thanks

This work would not have been possible without many other open-source libraries, so please also have a look at the following repositories and consider citing the corresponding articles:

• https://github.com/hadrienmichel/pyBEL1D
• https://github.com/emsig/empymod
• https://github.com/gimli-org/gimli
• https://github.com/zperzan/pyDGSA

Licenses

The data in this repository is licensed under CC-BY, the majority of the code is licensed under the apache 2.0 license, with the exception of everything related to pyBEL1D which uses the BSD-2-clause.

Citation

If you find this work useful and consider publishing related work, please consider citing our article ("Stochastic Inversion of Transient Electromagnetic Data to Derive Aquifer Geometry and Associated Uncertainties") that is accepted with minor revisions in GJI.