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 debris-covered 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).

Modern geoinformatic techniques allow the automated creation of detailed glacier inventory data from glacier outlines and digital terrain models (DTMs). Once glacier entities are defined and an appropriate DTM is available, several methods exist to derive the inventory data (e.g. minimum, maximum and mean elevation; mean slope and aspect) for each glacier from digital intersection of both datasets. Compared to the former manual methods, the new grid-based statistical calculations are very fast and reproducible. The major aim of this contribution is to help in standardizing the related calculations to enhance the integrity of the Global Land Ice Monitoring from Space (GLIMS) database. The recommendations were prepared by a working group and also contribute to the European Space Agency project GlobGlacier. The document follows the former UNESCO manual for the production of the World Glacier Inventory published in 1970, identifies the potential pitfalls, and describes the differences from the former methods of compilation. The online background material for this paper (see http://www.glims.org) contains example scripts for calculation of each parameter and will be updated when required.

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.

The ongoing glacier shrinkage in the Alps requires frequent updates of glacier outlines to provide an accurate database for monitoring, modelling purposes (e.g. determination of run-off, mass balance, or future glacier extent), and other applications. With the launch of the first Sentinel-2 (S2) satellite in 2015, it became possible to create a consistent, Alpine-wide glacier inventory with an unprecedented spatial resolution of 10 m. The first S2 images from August 2015 already provided excellent mapping conditions for most glacierized regions in the Alps and were used as a base for the compilation of a new Alpine-wide glacier inventory in a collaborative team effort. In all countries, glacier outlines from the latest national inventories have been used as a guide to compile an update consistent with the respective previous interpretation. The automated mapping of clean glacier ice was straightforward using the band ratio method, but the numerous debris-covered glaciers required intense manual editing. Cloud cover over many glaciers in Italy required also including S2 scenes from 2016. The outline uncertainty was determined with digitizing of 14 glaciers several times by all participants. Topographic information for all glaciers was obtained from the ALOS AW3D30 digital elevation model (DEM). Overall, we derived a total glacier area of 1806±60 km² when considering 4395 glaciers >0.01 km². This is 14 % (−1.2 % a⁻¹) less than the 2100 km² derived from Landsat in 2003 and indicates an unabated continuation of glacier shrinkage in the Alps since the mid-1980s. It is a lower-bound estimate, as due to the higher spatial resolution of S2 many small glaciers were additionally mapped or increased in size compared to 2003. Median elevations peak around 3000 m a.s.l., with a high variability that depends on location and aspect. The uncertainty assessment revealed locally strong differences in interpretation of debris-covered glaciers, resulting in limitations for change assessment when using glacier extents digitized by different analysts. The inventory is available at https://doi.org/10.1594/PANGAEA.909133 (Paul et al., 2019).

Meltwater from glaciers in the European Alps plays an important role in hydropower production, and future glacier development is thus of economic interest. However, an up-to-date and alpine-wide inventory for accurate assessment of glacier changes or modelling of future glacier development has not hitherto been available. Here we present a new alpine-wide inventory (covering Austria, France, Italy and Switzerland) derived from ten Landsat Thematic Mapper (TM) scenes acquired within 7 weeks in 2003. Combined with the globally available digital elevation model from the Shuttle Radar Topography Mission, topographic inventory parameters were derived for each of the 3770 mapped glaciers, covering 2050 km². The area-class frequency distribution is very similar in all countries, and a mean northerly aspect (NW, N, NE) is clearly favoured (arithmetic counting). Mean glacier elevation is ~2900 m, with a small dependence on aspect. The total area loss since the previous glacier inventory (acquired around 1970±15 years) is roughly one-third, yielding a current area loss rate of ~2%a^–1. Digital overlay of the outlines from the latest Austrian glacier inventory revealed differences in the interpretation of glacier extents that prohibit change assessment. A comparison of TM-derived outlines with manually digitized extents on a high-resolution IKONOS image returned 1.5% smaller glaciers with TM.