Albedo is one of the variables controlling the mass balance of temperate glaciers. Multispectral imagers, such as MODerate Imaging Spectroradiometer (MODIS) on board the TERRA and AQUA satellites, provide a means to monitor glacier surface albedo. In this study, different methods to retrieve broadband glacier surface albedo from MODIS data are compared. The effect of multiple reflections due to the rugged topography and of the anisotropic reflection of snow and ice are particularly investigated. The methods are tested on the Saint Sorlin Glacier (Grandes Rousses area, French Alps). The accuracy of the retrieved albedo is estimated using both field measurements, at two automatic weather stations located on the glacier, and albedo values derived from terrestrial photographs. For summers 2008 and 2009, the root mean square deviation (RMSD) between field measurements and the broadband albedo retrieved from MODIS data at 250 m spatial resolution was found to be 0.052 or about 10% relative error. The RMSD estimated for the MOD10 daily albedo product is about three times higher. One decade (2000–2009) of MODIS data were then processed to create a time series of albedo maps of Saint Sorlin Glacier during the ablation season. The annual mass balance of Saint Sorlin Glacier was compared with the minimum albedo value (average over the whole glacier surface) observed with MODIS during the ablation season. A strong linear correlation exists between the two variables. Furthermore, the date when the average albedo of the whole glacier reaches a minimum closely corresponds to the period when the snow line is located at its highest elevation, thus when the snow line is a good indicator of the glacier equilibrium line. This indicates that this strong correlation results from the fact that the minimal average albedo values of the glacier contains considerable information regarding the relative share of areal surfaces between the ablation zone (i.e. ice with generally low albedo values) and the accumulation zone (i.e. snow with a relatively high albedo). As a consequence, the monitoring of the glacier surface albedo using MODIS data can provide a useful means to evaluate the interannual variability of the glacier mass balance. Finally, the albedo in the ablation area of Saint Sorlin Glacier does not exhibit any decreasing trend over the study period, contrasting with the results obtained on Morteratsch Glacier in the Swiss Alps.

Few glaciological field data are available on glaciers in the Hindu Kush–Karakoram–Himalayan (HKH) region, and remote sensing data are thus critical for glacier studies in this region. The main objectives of this study are to document, using satellite images, the seasonal changes of surface albedo for two Himalayan glaciers, Chhota Shigri Glacier (Himachal Pradesh, India) and Mera Glacier (Everest region, Nepal), and to reconstruct the annual mass balance of these glaciers based on the albedo data. Albedo is retrieved from Moderate Resolution Imaging Spectroradiometer (MODIS) images, and evaluated using ground based measurements. At both sites, we find high coefficients of determination between annual minimum albedo averaged over the glacier (AMAAG) and glacier-wide annual mass balance (B_a) measured with the glaciological method (R² = 0.75). At Chhota Shigri Glacier, the relation between AMAAG found at the end of the ablation season and B_a suggests that AMAAG can be used as a proxy for the maximum snow line altitude or equilibrium line altitude (ELA) on winter-accumulation-type glaciers in the Himalayas. However, for the summer-accumulation-type Mera Glacier, our approach relied on the hypothesis that ELA information is preserved during the monsoon. At Mera Glacier, cloud obscuration and snow accumulation limits the detection of albedo during the monsoon, but snow redistribution and sublimation in the post-monsoon period allows for the calculation of AMAAG. Reconstructed B_a at Chhota Shigri Glacier agrees with mass balances previously reconstructed using a positive degree-day method. Reconstructed B_a at Mera Glacier is affected by heavy cloud cover during the monsoon, which systematically limited our ability to observe AMAAG at the end of the melting period. In addition, the relation between AMAAG and B_a is constrained over a shorter time period for Mera Glacier (6 years) than for Chhota Shigri Glacier (11 years). Thus the mass balance reconstruction is less robust for Mera Glacier than for Chhota Shigri Glacier. However our method shows promising results and may be used to reconstruct the annual mass balance of glaciers with contrasted seasonal cycles in the western part of the HKH mountain range since the early 2000s when MODIS images became available.

The characterization of snow extent is critical for a wide range of applications. Since 1966, snow maps at different spatial resolutions have been produced using various satellite sensor images. Nowadays, the most widely used products are likely those derived from Moderate-Resolution Imaging Spectroradiometer (MODIS) data, which cover the whole Earth at a near-daily frequency. There are a variety of snow mapping methods for MODIS data, based on different methodologies and applied at different spatial resolutions. Up to now, all these products have been tested and evaluated separately. This study aims to compare the methods currently available for retrieving snow from MODIS data. The focus is on fractional snow cover, which represents the snow cover area at the subpixel level. We examine the two main approaches available for generating such products from MODIS data; namely, linear regression of the Normalized Difference Snow Index (NDSI) and spectral unmixing (SU). These two approaches have resulted in several methods, such as MOD10A1 (the NSIDC MODIS snow product) for NDSI regression, and MODImLAB for SU. The assessment of these approaches was carried out using higher resolution binary snow maps (i.e., showing the presence or absence of snow) at spatial resolutions of 10, 20, and 30 m, produced by SPOT 4, SPOT 5, and LANDSAT-8, respectively. Three areas were selected in order to provide landscape diversity: the French Alps (117 dates), the Pyrenees (30 dates), and the Moroccan Atlas (24 dates). This study investigates the impact of reference maps on accuracy assessments, and it is suggested that NDSI-based high spatial resolution reference maps advantage NDSI medium-resolution snow maps. For MODIS snow maps, the results show that applying an NDSI approach to accurate surface reflectance corrected for topographic and atmospheric effects generally outperforms other methods for the global retrieval of snow cover area. The improvements to the newer version of MOD10A1 (Collection 6) compared to the older version (Collection 5) are significant. Products based on SU provide a good alternative and more accurate retrieval of the snow fraction where wider ranges of land covers are concerned. The fusion process and its resulting 250 m spatial resolution product improve snow line retrieval. False detection in mixed pixels, probably due to the spectral variability associated with the various materials in the spectral mixture, has been identified as an area that will require improvement.

Less than 0.25 % of the 250 000 glaciers inventoried in the Randolph Glacier Inventory (RGI V.5) are currently monitored with in situ measurements of surface mass balance. Increasing this archive is very challenging, especially using time-consuming methods based on in situ measurements, and complementary methods are required to quantify the surface mass balance of unmonitored glaciers. The current study relies on the so-called albedo method, based on the analysis of albedo maps retrieved from optical satellite imagery acquired since 2000 by the MODIS sensor, on board the TERRA satellite. Recent studies revealed substantial relationships between summer minimum glacier-wide surface albedo and annual surface mass balance, because this minimum surface albedo is directly related to the accumulation–area ratio and the equilibrium-line altitude. On the basis of 30 glaciers located in the French Alps where annual surface mass balance data are available, our study conducted on the period 2000–2015 confirms the robustness and reliability of the relationship between the summer minimum surface albedo and the annual surface mass balance. For the ablation season, the integrated summer surface albedo is significantly correlated with the summer surface mass balance of the six glaciers seasonally monitored. These results are promising to monitor both annual and summer glacier-wide surface mass balances of individual glaciers at a regional scale using optical satellite images. A sensitivity study on the computed cloud masks revealed a high confidence in the retrieved albedo maps, restricting the number of omission errors. Albedo retrieval artifacts have been detected for topographically incised glaciers, highlighting limitations in the shadow correction algorithm, although inter-annual comparisons are not affected by systematic errors.

We investigated snow cover dynamics using time series of moderate (MODIS) to high (SPOT-4/5, Landsat-8) spatial resolution satellite imagery in a 3700 km² region of the southwestern French Alps. Our study was carried out in the context of the SPOT (Take 5) Experiment initiated by the Centre National d’Etudes Spatiales (CNES), with the aim of exploring the utility of high spatial and temporal resolution multispectral satellite imagery for snow cover mapping and applications in alpine ecology. Our three objectives were: (i) to validate remote sensing observations of first snow free day derived from the Normalized Difference Snow Index (NDSI) relative to ground-based measurements; (ii) to generate regional-scale maps of first snow free day and peak standing biomass derived from the Normalized Difference Vegetation Index (NDVI); and (iii) to examine the usefulness of these maps for habitat mapping of herbaceous vegetation communities above the tree line. Imagery showed strong agreement with ground-based measurements of snow melt-out date, although R² was higher for SPOT and Landsat time series (0.92) than for MODIS (0.79). Uncertainty surrounding estimates of first snow free day was lower in the case of MODIS, however (±3 days as compared to ±9 days for SPOT and Landsat), emphasizing the importance of high temporal as well as high spatial resolution for capturing local differences in snow cover duration. The main floristic differences between plant communities were clearly visible in a two-dimensional habitat template defined by the first snow free day and NDVI at peak standing biomass, and these differences were accentuated when axes were derived from high spatial resolution imagery. Our work demonstrates the enhanced potential of high spatial and temporal resolution multispectral imagery for quantifying snow cover duration and plant phenology in temperate mountain regions, and opens new avenues to examine to what extent plant community diversity and functioning are controlled by snow cover duration.

Snow depth has traditionally been estimated based on point measurements collected either manually or at automated weather stations. Point measurements, though, do not represent the high spatial variability in snow depths present in alpine terrain. Photogrammetric mapping techniques have progressed in recent years and are capable of accurately mapping snow depth in a spatially continuous manner, over larger areas and at various spatial resolutions. However, the strengths and weaknesses associated with specific platforms and photogrammetric techniques as well as the accuracy of the photogrammetric performance on snow surfaces have not yet been sufficiently investigated. Therefore, industry-standard photogrammetric platforms, including high-resolution satellite (Pléiades), airplane (Ultracam Eagle M3), unmanned aerial system (eBee+ RTK with SenseFly S.O.D.A. camera) and terrestrial (single lens reflex camera, Canon EOS 750D) platforms, were tested for snow depth mapping in the alpine Dischma valley (Switzerland) in spring 2018. Imagery was acquired with airborne and space-borne platforms over the entire valley, while unmanned aerial system (UAS) and terrestrial photogrammetric imagery was acquired over a subset of the valley. For independent validation of the photogrammetric products, snow depth was measured by probing as well as by using remote observations of fixed snow poles. When comparing snow depth maps with manual and snow pole measurements, the root mean square error (RMSE) values and the normalized median absolute deviation (NMAD) values were 0.52 and 0.47 m, respectively, for the satellite snow depth map, 0.17 and 0.17 m for the airplane snow depth map, and 0.16 and 0.11 m for the UAS snow depth map. The area covered by the terrestrial snow depth map only intersected with four manual measurements and did not generate statistically relevant measurements. When using the UAS snow depth map as a reference surface, the RMSE and NMAD values were 0.44 and 0.38 m for the satellite snow depth map, 0.12 and 0.11 m for the airplane snow depth map, and 0.21 and 0.19 m for the terrestrial snow depth map. When compared to the airplane dataset over a large part of the Dischma valley (40 km²), the snow depth map from the satellite yielded an RMSE value of 0.92 m and an NMAD value of 0.65 m. This study provides comparative measurements between photogrammetric platforms to evaluate their specific advantages and disadvantages for operational, spatially continuous snow depth mapping in alpine terrain over both small and large geographic areas.