A Python multiprocessing approach for fast geostatistical simulations of subglacial topography

Realistically rough stochastic realizations of subglacial bed topography are crucial for improving our understanding of basal processes and quantifying uncertainty in sea-level rise projections with respect to topographic uncertainty. This can be achieved with Sequential Gaussian Simulation (SGS), which is used to generate multiple non-unique realizations of geological phenomena that sample the uncertainty space. However, SGS is very CPU intensive with a computational complexity of O( Nk ³ ), where N is the number of grid cells to simulate, and k is the number of neighboring points used for conditioning. This complexity makes SGS prohibitively time-consuming to implement at ice-sheet scales or fine resolutions. To reduce the time-cost, we implement and test a multiprocess version of SGS using Python’s multiprocessing module. By parallelizing the calculation of the weight parameters used in SGS, we achieve a speedup of 9.5 running on 16 processors for an N of 128,097. This speedup, as well as the speedup from using multiple processors, increases with N . This speed improvement makes SGS viable for large-scale topography mapping and ensemble ice-sheet modeling. Additionally, we have made our code repository and user tutorials publicly available (GitHub, Zenodo that others can use our multiprocess implementation of SGS on different datasets.