Impact of bedrock description on modeling ice sheet dynamics

Recent glaciological surveys have revealed a significant increase of ice discharge from polar ice caps into the ocean. In parallel, ice flow models have been greatly improved to better reproduce current changes and forecast the future behavior of ice sheets. For these models, surface topography and bedrock elevation are crucial input parameters that largely control the dynamics and the ensuing overall mass balance of the ice sheet. For obvious reasons of inaccessibility, only sparse and uneven bedrock elevation data is available. This raw data is processed to produce Digital Elevation Models (DEMs) on a regular 5 km grid. These DEMs are used to constrain the basal boundary conditions of all ice sheet models. Here, by using a full-Stokes finite element code, we examine the sensitivity of an ice flow model to the accuracy of the bedrock description. In the context of short-term ice sheet forecast, we show that in coastal regions, the bedrock elevation should be known at a resolution of the order of one kilometer. Conversely, a crude description of the bedrock in the interior of the continent does not affect modeling of the ice outflow into the ocean. These findings clearly indicate that coastal regions should be prioritized during future geophysical surveys. They also indicate that a paradigm shift is required to change the current design of DEMs describing the bedrock below the ice sheets: they must give users the opportunity to incorporate high-resolution bedrock elevation data in regions of interest. Citation: Durand, G., O. Gagliardini, L. Favier, T. Zwinger, and E. le Meur (2011), Impact of bedrock description on modeling ice sheet dynamics, Geophys. Res. Lett., 38, L20501, doi: 10.1029/2011GL048892.