On the Automated Mapping of Snow Cover on Glaciers and Calculation of Snow Line Altitudes from Multi-Temporal Landsat Data

Mapping snow cover (SC) on glaciers at the end of the ablation period provides a possibility to rapidly obtain a proxy for their equilibrium line altitude (ELA) which in turn is a metric for the mass balance. Satellite determination of glacier snow cover, derived over large regions, can reveal its spatial variability and temporal trends. Accordingly, snow mapping on glaciers has been widely applied using several satellite sensors. However, as glacier ice originates from compressed snow, both have very similar spectral properties and standard methods to map snow struggle to distinguish snow on glaciers. Hence, most studies applied manual delineation of snow extent on glaciers. Here we present an automated tool, named ‘ASMAG’ (automated snow mapping on glaciers), to map SC on glaciers and derive the related snow line altitude (SLA) for individual glaciers using multi-temporal Landsat satellite imagery and a digital elevation model (DEM). The method has been developed using the example of the Ötztal Alps, where an evaluation of the method is possible using field-based observations of the annual equilibrium line altitude (ELA) and the accumulation area ratio (AAR) measured for three glaciers for more than 30 years. The tool automatically selects a threshold to map snow on glaciers and robustly calculates the SLA based on the frequency distribution of elevation bins with more than 50% SC. The accuracy of the SC mapping was about 90% and the SLA was determined successfully in 80% of all cases with a mean uncertainty of ±19 m. When cloud-free scenes close to the date of the highest snowline are available, a good to very good agreement of SC ratios (SCR)/SLA with field data of AAR/ELA are obtained, otherwise values are systematically higher/lower as useful images were often acquired too early in the summer season. However, glacier specific differences are still well captured. Snow mapping on glaciers is impeded by clouds and their shadows or when fresh snow is covering the glaciers, so that more frequent image acquisitions (as now provided by Sentinel-2) would improve results.

[1]  J. Dedieu,et al.  25 years (1981–2005) of equilibrium-line altitude and mass-balance reconstruction on Glacier Blanc, French Alps, using remote-sensing methods and meteorological data , 2008, Journal of Glaciology.

[2]  W. Haeberli,et al.  Observation-Based Estimates of Global Glacier Mass Change and Its Contribution to Sea-Level Change , 2016, Surveys in Geophysics.

[3]  Antoine Rabatel,et al.  Changes in glacier equilibrium-line altitude in the western Alps from 1984 to 2010: evaluation by remote sensing and modeling of the morpho-topographic and climate controls , 2013 .

[4]  Georg Kaser,et al.  Contribution potential of glaciers to water availability in different climate regimes , 2010, Proceedings of the National Academy of Sciences.

[5]  Y. Arnaud,et al.  Slight mass gain of Karakoram glaciers in the early twenty-first century , 2012 .

[6]  Antoine Rabatel,et al.  Can the snowline be used as an indicator of the equilibrium line and mass balance for glaciers in the outer tropics? , 2012, Journal of Glaciology.

[7]  Dorothy K. Hall,et al.  Characterization of Snow and Ice Reflectance Zones On Glaciers Using Landsat Thematic Mapper Data , 1987, Annals of Glaciology.

[8]  Y. Arnaud,et al.  Reconstructing the mass balance of Brewster Glacier, New Zealand, usingMODIS-derived glacier-wide albedo , 2016 .

[9]  Hiroyuki Fujisada,et al.  ASTER DEM performance , 2005, IEEE Transactions on Geoscience and Remote Sensing.

[10]  Andreas Kääb,et al.  Perspectives on the production of a glacier inventory from multispectral satellite data in Arctic Canada: Cumberland Peninsula, Baffin Island , 2005, Annals of Glaciology.

[11]  Bryan L. Isacks,et al.  Spectral mixture analysis of Landsat thematic mapper images applied to the detection of the transient snowline on tropical Andean glaciers , 1999 .

[12]  J. Dedieu,et al.  Using remote-sensing data to determine equilibrium-line altitude and mass-balance time series: validation on three French glaciers, 1994–2002 , 2005 .

[13]  Rolf Weingartner,et al.  Toward mountains without permanent snow and ice: MOUNTAINS WITHOUT PERMANENT SNOW AND ICE , 2017 .

[14]  F. Paul,et al.  Comparison of mass balances for Vernagtferner, Oetzal Alps, as obtained from direct measurements and distributed modeling , 2009, Annals of Glaciology.

[15]  P. Teillet,et al.  On the Slope-Aspect Correction of Multispectral Scanner Data , 1982 .

[16]  T. Bolch,et al.  Multi-decadal mass balance series of three Kyrgyz glaciers inferred from transient snowline observations , 2017 .

[17]  R. D. Garg,et al.  Evaluation of vertical accuracy of open source Digital Elevation Model (DEM) , 2013, Int. J. Appl. Earth Obs. Geoinformation.

[18]  G. Kaser,et al.  A reanalysis of one decade of the mass balance series on Hintereisferner, Ötztal Alps, Austria: a detailed view into annual geodetic and glaciological observations , 2017 .

[19]  A. Lambrecht,et al.  Glacier changes in the Austrian Alps during the last three decades, derived from the new Austrian glacier inventory , 2007, Annals of Glaciology.

[20]  Dorothy K. Hall,et al.  Analysis of glacier facies using satellite techniques , 1991, Journal of Glaciology.

[21]  P. Holmlund,et al.  Historically unprecedented global glacier decline in the early 21st century , 2015 .

[22]  A. Lambrecht,et al.  Quantifying changes and trends in glacier area and volume in the Austrian Ötztal Alps (1969-1997-2006) , 2009 .

[23]  M. Minnaert.  The reciprocity principle in lunar photometry , 1941 .

[24]  P. Rastner,et al.  Historical analysis and visualization of the retreat of Findelengletscher, Switzerland, 1859–2010 , 2016 .

[25]  J. Dozier.  Spectral Signature of Alpine Snow Cover from the Landsat Thematic Mapper , 1989 .

[26]  E. Berthier,et al.  Seasonal changes in surface albedo of Himalayan glaciers from MODIS data and links with the annual mass balance , 2014 .

[27]  H. Rott.  Analyse der Schneeflächen auf Gletschern der tiroler Zentralalpen aus Landsat-Bildern , 1976 .

[28]  Georg Kaser,et al.  Measurements and Models of the Mass Balance of Hintereisferner , 1999 .

[29]  N. Eckert,et al.  Global glacier mass changes and their contributions to sea-level rise from 1961 to 2016 , 2019, Nature.

[30]  Georg Kaser,et al.  The Rofental: a high Alpine research basin (1890–3770 m a.s.l.) in the Ötztal Alps (Austria) with over 150 years of hydrometeorological and glaciological observations , 2017 .

[31]  Y. Arnaud,et al.  Linking glacier annual mass balance and glacier albedo retrieved from MODIS data , 2012 .

[32]  B. Markham,et al.  Summary of Current Radiometric Calibration Coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI Sensors , 2009 .

[33]  Edward Hanna,et al.  Identification of snow ablation rate, ELA, AAR and net mass balance using transient snowline variations on two Arctic glaciers , 2013, Journal of Glaciology.

[34]  Tobias Bolch,et al.  Climate change and glacier retreat in northern Tien Shan (Kazakhstan/Kyrgyzstan) using remote sensing data , 2007 .

[35]  Andreas Kääb,et al.  Glacier Remote Sensing Using Sentinel-2. Part II: Mapping Glacier Extents and Surface Facies, and Comparison to Landsat 8 , 2016, Remote. Sens..

[36]  N. Gorelick,et al.  Global Assessment of Supraglacial Debris‐Cover Extents , 2018, Geophysical Research Letters.

[37]  Y. Arnaud,et al.  Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas , 2012, Nature.

[38]  J. Graham Cogley,et al.  Geodetic and direct mass-balance measurements: comparison and joint analysis , 2009 .

[39]  J. Oerlemans,et al.  Assessment of interannual variations in the surface mass balance of 18 Svalbard glaciers from the Moderate Resolution Imaging Spectroradiometer/Terra albedo product , 2007 .

[40]  Richard R. Irish,et al.  Landsat 7 automatic cloud cover assessment , 2000, SPIE Defense + Commercial Sensing.

[41]  E. Rodríguez,et al.  A Global Assessment of the SRTM Performance , 2006 .

[42]  Matthias Huss,et al.  Re-analysis of seasonal mass balance at Abramov glacier 1968-2014 , 2015 .

[43]  Omm.  Seasonal snow cover , 1970 .

[44]  C. Rixen,et al.  Shorter snow cover duration since 1970 in the Swiss Alps due to earlier snowmelt more than to later snow onset , 2016, Climatic Change.

[45]  M. Pelto.  Utility of late summer transient snowline migration rate on Taku Glacier, Alaska , 2011 .

[46]  N. Otsu.  A threshold selection method from gray level histograms , 1979 .

[47]  Pascal Peduzzi,et al.  Assessing high altitude glacier thickness, volume and area changes using field, GIS and remote sensing techniques: the case of Nevado Coropuna (Peru) , 2010 .

[48]  Andrea Fischer,et al.  Glaciers and climate change: Interpretation of 50 years of direct mass balance of Hintereisferner , 2010 .

[49]  Gerhard Krieger,et al.  Generation and performance assessment of the global TanDEM-X digital elevation model , 2017 .

[50]  W. Krabill,et al.  Penetration depth of interferometric synthetic‐aperture radar signals in snow and ice , 2001, Geophysical Research Letters.

[51]  A. Kääb,et al.  Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change , 2011 .

[52]  L. A. Rasmussen,et al.  Glossary of glacier mass balance and related terms , 2010 .

[53]  Richard Bamler,et al.  The SRTM Mission - A World-Wide 30 m Resolution DEM from SAR Interferometry in 11 Days , 1999 .

[54]  Xuehong Chen,et al.  Comparison of automatic thresholding methods for snow-cover mapping using Landsat TM imagery , 2013 .

[55]  S. Ekstrand,et al.  Landsat TM-based forest damage assessment : correction for topographic effects , 1996 .

[56]  J. Simões,et al.  Variations in annual snowline and area of an ice-covered stratovolcano in the Cordillera Ampato, Peru, using remote sensing data (1986–2014) , 2016 .

[57]  Masamu Aniya,et al.  The use of satellite and airborne imagery to inventory outlet glaciers of the Southern Patagonia Icefield, South America , 1996 .

[58]  G. Kaser,et al.  Seventy‐six years of mean mass balance rates derived from recent and re‐evaluated ice volume measurements on tropical Lewis Glacier, Mount Kenya , 2011 .

[59]  J. Oerlemans,et al.  Climate-model induced differences in the 21st century global and regional glacier contributions to sea-level rise , 2013, Climate Dynamics.

[60]  J. Dedieu,et al.  Spatio-temporal changes in glacier-wide mass balance quantified by optical remote sensing on 30 glaciers in the French Alps for the period 1983–2014 , 2016, Journal of Glaciology.

[61]  M. Hantel,et al.  Sensitivity of Alpine snow cover to European temperature , 2007 .

[62]  B. Menounos,et al.  An approach to derive regional snow lines and glacier mass change from MODIS imagery, western North America , 2013 .

[63]  K. Schulz,et al.  On the need for a time- and location-dependent estimation of the NDSI threshold value for reducing existing uncertainties in snow cover maps at different scales , 2017 .

[64]  Albert Rango,et al.  Temperature and emissivity separation from multispectral thermal infrared observations , 2002 .

[65]  D. Scherer,et al.  Modis derived equilibrium line altitude estimates for purogangri ice cap, tibetan plateau, and their relation to climatic predictors (2001–2012) , 2015 .

[66]  Renaud Mathieu,et al.  Detecting the equilibrium‐line altitudes of New Zealand glaciers using ASTER satellite images , 2009 .