Stochastic modeling of subglacial topography exposes uncertainty in water routing at Jakobshavn Glacier

Abstract Subglacial topography is an important feature in numerous ice-sheet analyses and can drive the routing of water at the bed. Bed topography is primarily measured with ice-penetrating radar. Significant gaps, however, remain in data coverage that require interpolation. Topographic interpolati...

Full description

Bibliographic Details
Published in:Journal of Glaciology
Main Authors: MacKie, Emma J., Schroeder, Dustin M., Zuo, Chen, Yin, Zhen, Caers, Jef
Format: Article in Journal/Newspaper
Language:English
Published: Cambridge University Press (CUP) 2020
Subjects:
Online Access:http://dx.doi.org/10.1017/jog.2020.84
https://www.cambridge.org/core/services/aop-cambridge-core/content/view/S0022143020000842
Description
Summary:Abstract Subglacial topography is an important feature in numerous ice-sheet analyses and can drive the routing of water at the bed. Bed topography is primarily measured with ice-penetrating radar. Significant gaps, however, remain in data coverage that require interpolation. Topographic interpolations are typically made with kriging, as well as with mass conservation, where ice flow dynamics are used to constrain bed geometry. However, these techniques generate bed topography that is unrealistically smooth at small scales, which biases subglacial water flowpath models and makes it difficult to rigorously quantify uncertainty in subglacial drainage patterns. To address this challenge, we adapt a geostatistical simulation method with probabilistic modeling to stochastically simulate bed topography such that the interpolated topography retains the spatial statistics of the ice-penetrating radar data. We use this method to simulate subglacial topography using mass conservation topography as a secondary constraint. We apply a water routing model to each of these realizations. Our results show that many of the flowpaths significantly change with each topographic realization, demonstrating that geostatistical simulation can be useful for assessing confidence in subglacial flowpaths.