EGU22-8167, updated on 28 Mar 2022
EGU General Assembly 2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.

The GeoStat-Framework: Create your geostatistical model in Python!

Sebastian Müller1,2, Lennart Schüler1,3, Alraune Zech1,4, and Falk Heße1,2
Sebastian Müller et al.
  • 1Helmholtz Centre for Environmental Research - UFZ, Computational Hydrosystems, Leipzig, Germany
  • 2Institute of Earth and Environmental Sciences, University Potsdam, Potsdam, Germany
  • 3Center for Advanced Systems Understanding (CASUS), Görlitz, Germany
  • 4Department of Earth Sciences, University Utrecht, Utrecht, Netherlands

The GeoStat Framework is a coherent ecosystem of Python packages for geostatistical applications and subsurface simulations. 
It provides an easily usable open source collection of software packages. They are well documented, including exhaustive hands-on guides and examples for helping non-programmers and non-domain-experts getting started. The main applications of the packages are:

  • GSTools & PyKrige - spatial random field generation, kriging, data normalization & transformation, and geostatistical analyses based on variogram methods
  • AnaFlow - (semi-)analytical solutions for specific groundwater-flow scenarios and the extended Generalized Radial Flow model
  • welltestpy - store, manipulate, and analyze well-based field testing campaigns with a focus on estimating parameters of subsurface heterogeneity from pumping test data.
  • ogs5py - pre-processing, operating, and post-processing of subsurface flow and transport simulations by providing a Python API for the FEM solver OpenGeoSys 5

With this collection of flexible toolboxes we aim to close the gap of missing software for real-world applications in the field of geostatistics. 
Especially GSTools is the first comprehensive Python-toolbox for covariance models, field generation, kriging, variogram estimation, data normalization and transformation. The formerly independed project PyKrige, developed by Benjamin Murphy, has been migrated to the GeoStat-Framework und we started to build a common Rust backend for the numerical heavy lifting called GSTools-Core.

We will show a set of complex workflow examples, like temperature trend analysis and pumping test ensemble simulations, and we will give an overview of the provided functionality and an outlook for the future. All workflows are made accessible as GeoStat-Examples repositories.


  • GSTools: Müller, S., et. al., 2021 (in review)
  • AnaFlow: Müller, S., et. al., 2021
  • welltestpy: Müller, S., et. al., 2021
  • ogs5py: Müller, S., et. al., 2020


  • Website:
  • GitHub:
  • GeoStat-Examples:

How to cite: Müller, S., Schüler, L., Zech, A., and Heße, F.: The GeoStat-Framework: Create your geostatistical model in Python!, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8167,, 2022.