Hydrogeophysics is interdisciplinary by definition. As researchers strive to gain quantitative information on process-relevant subsurface parameters while integrating non-geophysical measurements, multi-physical geoscientific models are often developed that simulate the dynamic process and its geophysical response. Such endeavors are associated with considerable technical challenges due to coupling of different numerical models, which represents an initial hurdle for students and many practitioners. Even technically versatile users often end up with individually tailored solutions at the cost of scientific reproducibility.
We argue that the reproducibility of studies in computational hydrogeophysics, and therefore the advancement of the field itself, needs versatile open-source software. One example is pyGIMLi - a flexible and computationally efficient framework for modeling and inversion in geophysics. The library provides management for structured and unstructured 2D and 3D meshes, finite-element and finite-volume solvers, various geophysical forward operators, as well as a generalized Gauss-Newton based inversion framework.
In this contribution, we highlight some of the recent advances and use cases in research and education since its 1.0 release in 2017 (Rücker et al., 2017) including:
- generalized modeling and inversion frameworks for conventional, joint, time-lapse and process-based inversion
- geostatistical regularization operators for unstructured meshes (Jordi et al., 2018)
- improved constraints in the presence of petrophysical parameter transformations demonstrated by an estimation of water, ice, and air in partially frozen systems (Wagner et al., 2019)
- 3D visualization leveraging upon PyVista (Sullivan and Kaszynski, 2019)
- simulation of electrical streaming potentials
- complex-valued forward modeling and inversion of induced polarization
- forward modeling with anisotropic parameters
- availability for Mac OS
- improved API and documentation
Since the library is freely available and platform-compatible, it is also well suited for teaching. We demonstrate examples from Master level university courses and public outreach, where learners can interactively change model and acquisition parameters to study their influence on a hydrogeophysical process simulation. Finally, we would like to use this opportunity to discuss future developments with the community.
Jordi, C., Doetsch, J., Günther, T., Schmelzbach, C., & Robertsson, J. O. (2018). Geostatistical regularization operators for geophysical inverse problems on irregular meshes. Geophysical Journal International, 213(2), 1374–1386.
Rücker, C., Günther, T., Wagner, F.M., 2017. pyGIMLi: An open-source library for modelling and inversion in geophysics, Computers and Geosciences, 109, 106-123.
Sullivan, C., & Kaszynski, A. (2019). PyVista: 3D plotting and mesh analysis through a streamlined interface for the Visualization Toolkit (VTK). Journal of Open Source Software, 4(37), 1450.
Wagner, F. M., Mollaret, C., Günther, T., Kemna, A., & Hauck, C. (2019). Quantitative imaging of water, ice and air in permafrost systems through petrophysical joint inversion of seismic refraction and electrical resistivity data. Geophysical Journal International, 219(3), 1866–1875.