ABSTRACT. Deriving glacier outlines from satellite data has become increasingly popular in the past decade. In particular when glacier outlines are used as a base for change assessment, it is important to know how accurate they are. Calculating the accuracy correctly is challenging, as appropriate reference data (e.g. from higher-resolution sensors) are seldom available. Moreover, after the required manual correction of the raw outlines (e.g. for debris cover), such a comparison would only reveal the accuracy of the analyst rather than of the algorithm applied. Here we compare outlines for clean and debriscovered glaciers, as derived from single and multiple digitizing by different or the same analysts on very high-(1 m) and medium-resolution (30 m) remote-sensing data, against each other and to glacier outlines derived from automated classification of Landsat Thematic Mapper data. Results show a high variability in the interpretation of debris-covered glacier parts, largely independent of the spatial resolution (area differences were up to 30%), and an overall good agreement for clean ice with sufficient contrast to the surrounding terrain (differences $5%). The differences of the automatically derived outlines from a reference value are as small as the standard deviation of the manual digitizations from several analysts. Based on these results, we conclude that automated mapping of clean ice is preferable to manual digitization and recommend using the latter method only for required corrections of incorrectly mapped glacier parts (e.g. debris cover, shadow).

Abstract:With its temporal resolution of 10 days (five days with two satellites, and significantly more at high latitudes), its swath width of 290 km, and its 10 m and 20 m spatial resolution bands from the visible to the shortwave infrared, the European Sentinel-2 satellites have significant potential for glacier remote sensing, in particular mapping of glacier outlines and facies, and velocity measurements. Testing Level 1C commissioning and ramp-up phase data for initial sensor quality experiences, we find a high radiometric performance, but with slight striping effects under certain conditions. Through co-registration of repeat Sentinal-2 data we also find lateral offset patterns and noise on the order of a few metres. Neither of these issues will complicate most typical glaciological applications. Absolute geo-location of the data investigated was on the order of one pixel at the time of writing. The most severe geometric problem stems from vertical errors of the DEM used for ortho-rectifying Sentinel-2 data. These errors propagate into locally varying lateral offsets in the images, up to several pixels with respect to other georeferenced data, or between Sentinel-2 data from different orbits. Finally, we characterize the potential and limitations of tracking glacier flow from repeat Sentinel-2 data using a set of typical glaciers in different environments: Aletsch Glacier, Swiss Alps; Fox Glacier, New Zealand; Jakobshavn Isbree, Greenland; Antarctic Peninsula at the Larsen C ice shelf.

Abstract:Mapping of glacier extents from automated classification of optical satellite images has become a major application of the freely available images from Landsat. A widely applied method is based on segmented ratio images from a red and shortwave infrared band. With the now available data from Sentinel-2 (S2) and Landsat 8 (L8) there is high potential to further extend the existing time series (starting with Landsat 4/5 in 1982) and to considerably improve over previous capabilities, thanks to increased spatial resolution and dynamic range, a wider swath width and more frequent coverage. Here, we test and compare a variety of previously used methods to map glacier extents from S2 and L8, and investigate the mapping of snow facies with S2 using top of atmosphere reflectance. Our results confirm that the band ratio method works well with S2 and L8. The 15 m panchromatic band of L8 can be used instead of the red band, resulting in glacier extents similar to S2 (0.7% larger for 155 glaciers). On the other hand, extents derived from the 30 m bands are 4%-5% larger, indicating a more generous interpretation of mixed pixels. Mapping of snow cover with S2 provided accurate results, but the required topographic correction would benefit from a better orthorectification with a more precise DEM than currently used.

Abstract. In this study, we assess glacier area and length changes in mainland Norway from repeat Landsat TM/ETM+-derived inventories and digitized topographic maps. The multi-temporal glacier inventory consists of glacier outlines from three time ranges: 1947 to 1985 (GI n50 ), 1988 to 1997 (GI 1990), and 1999 to 2006 (GI 2000. For the northernmost regions, we include an additional inventory (GI 1900 ) based on historic maps surveyed between 1895 and 1907. Area and length changes are assessed per glacier unit, 36 subregions, and for three main parts of Norway: southern, central, and northern. The results show a decrease in the glacierized area from 2994 km 2 in GI n50 to 2668 km 2 in GI 2000 (total 2722 glacier units), corresponding to an area reduction of −326 km 2 , or −11 % of the initial GI n50 area. The average length change for the full epoch (within GI n50 and GI 2000 ) is −240 m. Overall, the comparison reveals both area and length reductions as general patterns, even though some glaciers have advanced. The three northernmost subregions show the highest retreat rates, whereas the central part of Norway shows the lowest change rates. Glacier area and length changes indicate that glaciers in maritime areas in southern Norway have retreated more than glaciers in the interior, and glaciers in the north have retreated more than southern glaciers. These observed spatial trends in glacier change are related to a combination of several factors such as glacier geometry, elevation, and continentality, especially in southern Norway.