Preprints
https://doi.org/10.5194/egusphere-2022-1224
https://doi.org/10.5194/egusphere-2022-1224
 
18 Nov 2022
18 Nov 2022
Status: this preprint is open for discussion.

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

Emma Johanne MacKie1, Michael Field1, Lijing Wang2, Zhen Yin2, Nathan Schoedl3,4,5, Matthew Hibbs3, and Allan Zhang5 Emma Johanne MacKie et al.
  • 1Department of Geological Sciences, University of Florida, Gainesville, 32601, USA
  • 2Department of Geological Sciences, Stanford University, Stanford, 94301, USA
  • 3Department of Computer and Information Science and Engineering, University of Florida, Gainesville, 32601, USA
  • 4Department of Mathematics, University of Florida, Gainesville, 32601, USA
  • 5Department of Statistics, University of Florida, Gainesville, 32601, USA

Abstract. The interpolation of geospatial phenomena is a common problem in Earth sciences applications that can be addressed with geostatistics, where spatial correlations are used to constrain interpolations. In certain applications, it can be particularly useful to perform geostatistical simulation, which is used to generate multiple non-unique realizations that reproduce the variability of measurements and are constrained by observations. Despite the broad utility of this approach, there are few open-access geostatistical simulation software. To address this accessibility issue, we present GStatSim, a Python package for performing geostatistical interpolation and simulation. GStatSim is distinct from previous geostatistics 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.

Emma Johanne MacKie et al.

Status: open (until 13 Jan 2023)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse

Emma Johanne MacKie et al.

Emma Johanne MacKie et al.

Viewed

Total article views: 204 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
136 63 5 204 5 4
  • HTML: 136
  • PDF: 63
  • XML: 5
  • Total: 204
  • BibTeX: 5
  • EndNote: 4
Views and downloads (calculated since 18 Nov 2022)
Cumulative views and downloads (calculated since 18 Nov 2022)

Viewed (geographical distribution)

Total article views: 194 (including HTML, PDF, and XML) Thereof 194 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
1
 
 
 
 
Latest update: 06 Dec 2022
Download
Short summary
Earth scientists often have to fill in spatial gaps in measurements. This gap-filling or interpolation can be accomplished with geostatistical methods, where the statistical relationships between measurements are used to inform how these gaps should be filled. Despite the broad utility of these methods, there are few freely available geostatistics software. We present GStatSim, a Python package for performing different geostatistical interpolation methods.