In the future, data from the moderate resolution imaging spectroradiometer (MODIS) will be employed to map snow in an automated environment at a resolution of 250 m to 1 km. Using Landsat thematic mapper (TM) data, an algorithm, SNOMAP, has been developed to map snow-covered area. This algorithm will be used, with appropriate modification, with MODIS data following the launch of the first Earth Observing System (EOS) platform in 1998. SNOMAP has been shown to be successful in mapping snow in a variety of areas using TM data. However, significant errors may be present in mountainous areas due to effects of topography. To increase the accuracy of mapping global snow-covered area in the future using MODIS data, digital elevation model (DEM) data have been registered to TM data for parts of Glacier National Park, Montana, so that snow cover on mountain slopes can be mapped. This paper shows that the use of DEM data registered to TM data increases the accuracy of mapping snow-covered area. Using SNOMAP on a subscene within the 14 March 1991 TM scene of northwestern Montana, 215 km2 of snow is mapped when TM data are used alone to map the snow cover. We show that about 1062 km2 of snow are actually present as measured when the TM and DEM data are registered. Approximately five times more snow is present when the effects of topography are considered for this subscene, which is in a rugged area in Glacier National Park. A simple model has been developed to determine the relationship between terrain relief and the amount of correction that must be applied to map actual snow-covered area in Glacier National Park using satellite data alone.