An analysis of balance velocities over the Greenland ice sheet and comparison with synthetic aperture radar interferometry

Balance velocities for the Greenland ice sheet have been calculated from a new digital elevation model (DEM), accumulation rates and an existing ice-thickness grid, using a fully two-dimensional finite-difference scheme. The pattern of velocities is compared with velocities derived from synthetic-aperture radar (SAR) interferometry for three different regions of the ice sheet. Differences between the two estimates of velocity highlight the respective strengths and weaknesses of the datasets and techniques used. A comparison with ten global positioning system-derived velocities indicates that the balance-velocity scheme and input datasets used here provide a remarkably good representation of the velocity distribution inland from the margins. These balance-velocity data, therefore, could help constrain numerical ice-sheet models. The balance velocities were found to be unreliable dose to the ice-sheet margins due to larger errors in ice thickness, surface slope and ablation rate in this region. Comparison of the balance velocities with SAR interferometry in the region of the "Northeast Greenland Ice Stream" indicates the importance of the smoothing distance that must be applied to the DEM before calculating balance velocities. A smoothing distance of 20 times the ice thickness gave good agreement between the two measures of velocity.