A new digital elevation model of the Greenland ice sheet and surrounding rock outcrops has been produced at 1-km postings from a comprehensive suite of satellite remote sensing and cartographic data sets. Height data over the ice sheet were mainly from ERS-1 and Geosat radar altimetry. These data were corrected for a slope-dependent bias that had been identified in a previous study. The radar altimetry was supplemented with stereophotogrammetric data sets, synthetic aperture radar interferometry, and digitized cartographic maps over regions of bare rock and where gaps in the satellite altimeter coverage existed. The data were interpolated onto a regular grid with a spacing of ~1 km. The accuracy of the resultant digital elevation model over the ice sheet was assessed using independent and spatially extensive measurements from an airborne laser altimeter that had an accuracy of between 10 and 12 cm. In a comparison with the laser altimetry the digital elevation model was found to have a slope-dependent accuracy ranging from -1.04¿1.98 m to -0.06¿14.33 m over the ice sheet for a slope range of 0.0--1.0¿. The mean accuracy over the whole ice sheet was -0.33¿6.97 m. Over the bare rock areas the accuracy ranged from 20 to 200 m, dependent on the data source available. The new digital elevation model was used as an input data set for a positive degree day model of ablation. The new elevation model was found to reduce ablation by only 2% compared with using an older, 2.5-km resolution model, which suggests that resolution-induced errors in estimating ablation are no longer important. ¿ 2001 American Geophysical Union