A Novel Method for Reducing Uncertainty in Subglacial Topography: Implications for Greenland Ice Sheet Volume and Stability

Subglacial topography is a critical boundary condition for ice sheet models projecting past and future ice sheet–climate interactions. Contemporary ice-sheet-wide bed topography datasets are partially derived using mass conservation, but approximately 75% of the most widely used Greenland Ice Sheet (GrIS) dataset is based on simple interpolation of airborne radio-echo sounding (RES) measurements, such as kriging or streamline diffusion. Due to limited independent validation data, the errors and biases in this approach are poorly understood, creating largely unknown uncertainties in subglacial topography. Here, we interpolated synthetic RES observations of bed topography over ice-free areas with a known topography at a 5 m spatial resolution and quantify discrepancies. We found that the absolute error in kriged bed topography increases with distance from the input data, though at a reduced rate than previously estimated. The difference between an interpolated elevation estimate and the local mean elevation is a strong predictor of real bed errors (R2 = 0.72), with further improvement as input observation sparsity increases (R2 > 0.82). We propose a method to quantify and reduce uncertainty in kriged bed topography in sparsely surveyed regions, reducing uncertainty for at least 56% of the kriged interior at a 99% confidence interval. Our results suggest that subglacial depth is on average 5 m deeper than previous estimates, though individual areas may be shallower or deeper (σ = 41 m). Consequently, the area grounded below sea level is likely underestimated by 2%, increasing to 29% for regions deeper than 200 m. These findings have potential implications for the future stability of the GrIS under climate change.