GStatSim V1.0: a Python package for geostatistical interpolation and conditional simulation

By Emma MacKie

University of Florida

GStatSim is a Python package specifically designed for geostatistical interpolation and simulation. These tools are part of our ongoing effort to develop and adapt open-access geostatistical functions.

Launch Tool

You must login before you can run this tool.

Version 1.0 - published on 16 Oct 2023

Open source: license | download

View All Supporting Documents



Published on


The interpolation of geospatial phenomena is a common problem in Earth science applications that can be addressed with geostatistics, where spatial correlations are used to constrain interpolations. In certain applications, it can be particularly useful to a perform geostatistical simulation, which is used to generate multiple non-unique realizations that reproduce the variability in measurements and are constrained by observations. Despite the broad utility of this approach, there are few open-access geostatistical simulation software applications. To address this accessibility issue, we present GStatSim, a Python package for performing geostatistical interpolation and simulation. GStatSim is distinct from previous geostatistical tools in that it emphasizes accessibility for non-experts, geostatistical simulation, and applicability to remote sensing data sets. It includes tools for performing non-stationary simulations and interpolations with secondary constraints. This package is accompanied by a Jupyter Book with user tutorials and background information on different interpolation methods. These resources are intended to significantly lower the technological barrier to using geostatistics and encourage the use of geostatistics in a wider range of applications. We demonstrate the different functionalities of this tool for the interpolation of subglacial topography measurements in Greenland.

Code can also be found on GitHub:



This work was sponsored by the Earth Sciences Information Partners (ESIP) Lab, the National Science Foundation (NSF) GeoSMART program, the University of Florida Artificial Intelligence Initiative Center, and the National Science Foundation GEO OSE award #2324092.

Cite this work

Researchers should cite this work as follows:

  • MacKie, E. J., Field, M., Wang, L., Yin, Z., Schoedl, N., Hibbs, M., & Zhang, A. (2023). GStatSim V1. 0: a Python package for geostatistical interpolation and conditional simulation. Geoscientific Model Development, 16(13), 3765-3783. DOI -
  • Emma MacKie (2023), "GStatSim V1.0: a Python package for geostatistical interpolation and conditional simulation,"

    BibTex | EndNote