Stochastic Sub-ice-shelf Bathymetry for Thwaites, Crosson, and Dotson Ice Shelves

This repository contains the code to generate stochastic realizations of sub-ice-shelf bathymetry for Thwaites Eastern Ice Shelf, Dotson Ice Shelf, and Crosson Ice shelf, as well as 3 ensembles of realizations. Ensemble 1 includes geostatistical simulation of the Bouguer gravity disturbance and random background densities, ensemble 2 includes geostatistical simulation of the Bouguer disturbance but no background density sampling, and ensemble 3 uses a constant mean Bouguer disturbance and constant background density. The ensembles were generated at 2000 m resolution and were upscaled onto the original BedMachine 500 m grid. The results are in the following NetCDF files: Ensemble 1: dens_geoid_500.nc dens_geoid_2000.nc Ensemble 2: nodens_geoid_500.nc nodens_geoid_2000.nc Ensemble 3: krige_geoid_500.nc krige_geoid_2000.nc In order to reproduce the results see https://github.com/mjfield2/stochastic_bathymetry for usage instructions.