Geometric boundary conditions for modelling the velocity field of the Antarctic ice sheet

Abstract

This paper presents improved geometric boundary conditions (surface elevation and ice thickness) required as inputs to calculations of the surface-velocity field for the Antarctic ice sheet. A comparison of the two-dimensional horizontal velocity field obtained on the basis of conservation of mass (balance velocity) with the diagnostic velocity field calculated with an ice-sheet model (dynamic velocity) may yield information on shortcomings in the way the ice-sheet model describes the ice flow. Here, the surface-elevation grid is described in detail, as it has been generated specifically for such a study and represents a new standard in accuracy and resolution for calculating surface slopes. The digital-elevation model was generated on a 10 km grid size from over 20 000 000 height estimates obtained from eight 35 d repeat cycles of ERS-1 radar-altimeter data. For surface slopes less than 0.4°, the accuracy is better than 1.5 m. In areas of high surface slope (coastal and mountainous regions), the altimeter measurements have been supplemented with data taken from the Antarctic Digital Database. South of 81.5°, data from the SPRI folio map have been used. The ice-thickness grid was produced from a combination of a redigitization of the SPRI folio and the original radio-echo-sounding flight lines. For areas of grounded ice, the elevation of the bed was estimated from surface elevation and ice thickness. Significant differences (in excess of 25% of ice thickness) were obtained between an earlier digitization of the folio bed-elevation map and the data set derived here. Furthermore, a new value of 25.6 × 106 km3 was obtained for the total volume of the ice sheet and ice shelves, which is a reduction of 12% compared with the original estimate derived during the compilation of the SPRI folio. These differences will have an important influence on the results obtained by numerical ice-sheet models.