Change Detection of Bare-Ice Albedo in the Swiss Alps
Total Page:16
File Type:pdf, Size:1020Kb
The Cryosphere, 13, 397–412, 2019 https://doi.org/10.5194/tc-13-397-2019 © Author(s) 2019. This work is distributed under the Creative Commons Attribution 4.0 License. Change detection of bare-ice albedo in the Swiss Alps Kathrin Naegeli1,2, Matthias Huss3,4, and Martin Hoelzle3 1Centre for Glaciology, Department of Geography and Earth Sciences, Aberystwyth University, Aberystwyth, Wales, UK 2Institute of Geography and Oeschger Centre for Climate Change Research, University of Bern, Bern, Switzerland 3Department of Geosciences, University of Fribourg, Fribourg, Switzerland 4Laboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich, Zurich, Switzerland Correspondence: Kathrin Naegeli ([email protected]) Received: 21 January 2018 – Discussion started: 29 January 2018 Revised: 14 November 2018 – Accepted: 7 January 2019 – Published: 1 February 2019 Abstract. Albedo feedback is an important driver of glacier 1 Introduction melt over bare-ice surfaces. Light-absorbing impurities strongly enhance glacier melt rates but their abundance, com- position and variations in space and time are subject to con- Glaciers are known to be excellent indicators of climate siderable uncertainties and ongoing scientific debates. In this change (IPCC, 2013). Increasing air temperatures and chang- study, we assess the temporal evolution of shortwave broad- ing precipitation patterns provoke snow lines to rise to higher band albedo derived from 15 end-of-summer Landsat scenes altitudes and thus a spatially greater exposure of bare-ice sur- for the bare-ice areas of 39 large glaciers in the western and faces. In connection with a general prolongation of the ab- southern Swiss Alps. Trends in bare-ice albedo crucially de- lation season, the increased climatic forcing causes an am- pend on the spatial scale considered. No significant negative plification of glacier melt (Kuhn, 1993; Oerlemans, 2001). temporal trend in bare-ice albedo was found on a regional However, these changing glacier characteristics trigger feed- to glacier-wide scale. However, at higher spatial scales, cer- back mechanisms, in particular positive albedo feedback, tain areas of bare ice, including the lowermost elevations and which enhances bare-ice melting (Tedesco et al., 2011; Box margins of the ablation zones, revealed significant darken- et al., 2012). Hence, the strongly negative mass balances of ing over the study period 1999 to 2016. A total glacier area many glaciers are not solely a direct signal of atmospheric of 13.5 km2 (equivalent to about 12 % of the average end- warming but result from a complex interplay of changes of-summer bare-ice area in the study area) exhibited albedo in climate forcing and related surface–atmosphere feedback trends significant at the 95 % confidence level or higher. Most mechanisms. of this area was affected by a negative albedo trend of about Currently, there is an ongoing debate about the occurrence −0.05 decade−1. Generally, bare-ice albedo exhibits a strong and rate of glacier and ice sheet darkening worldwide. While interannual variability, caused by a complex interplay of me- studies like those of Oerlemans et al.(2009), Wang et al. teorological conditions prior to the acquisition of the data, (2014), Takeuchi(2001) or Mernild et al.(2015) observed a local glacier characteristics and the date of the investigated darkening for one or several glaciers, in the European Alps, satellite imagery. Although a darkening of glacier ice was the Chinese and Nepalese Himalayas or in Greenland’s pe- found to be present over only a limited region, we emphasize ripheral glacierized areas over varying timescales, respec- that due to the recent and projected growth of bare-ice ar- tively, evidence of darkening from sectors of the Greenland eas and prolongation of the ablation season in the region, the Ice Sheet is less pronounced, leading to controversial discus- albedo feedback will considerably enhance the rate of glacier sions (cf. Box et al., 2012; Alexander et al., 2014; Dumont mass loss in the Swiss Alps in the near future. et al., 2014; Polashenski et al., 2015; Tedesco et al., 2016). The recalibration of the MODIS sensors leads to a reduction in spatial extent and statistical strength of albedo trends over the Greenland Ice Sheet (Casey et al., 2017). Moreover, the emergence of legacy contaminants, radionuclides or heavy Published by Copernicus Publications on behalf of the European Geosciences Union. 398 K. Naegeli et al.: Change detection of bare-ice albedo in the Swiss Alps metals contained in cryoconite holes at lower elevations on in the western and southern Swiss Alps from the local point Alpine glaciers (Bogdal et al., 2009; Pavlova et al., 2014; to the regional scale. Causes and external factors that might Steinlin et al., 2014, 2016; Baccolo et al., 2017) or outcrop- impact bare-ice albedo and explain its spatial and temporal ping ice containing high dust concentrations potentially as- evolution are discussed. sociated with paleo-climatic conditions in the ablation zone of the Greenland Ice Sheet (Wientjes and Oerlemans, 2010; Wientjes et al., 2012; Goelles and Boggild, 2017) may em- 2 Study sites and data phasize the impact of light-absorbing impurities on ice melt. In addition, the recent recognition of the potential role of bio- Our study focuses on 39 glaciers located in the western and logical impurities (Cook et al., 2017; Stibal et al., 2017; Ted- southern Swiss Alps (Fig.1). All of them are characterized stone et al., 2017) further indicates an amplification of this by a surface area of roughly 5 km2 or larger, and thus of- linkage between impurities and supra-glacial processes. fer a large-enough spatial extent to study the evolution of To date, most long-term studies either used point data bare-ice surfaces and related albedo changes. The investi- from automatic weather stations located in the ablation area gated glaciers vary considerably in size, ranging from about of a glacier (e.g. Oerlemans et al., 2009), coarsely spaced 5 km2 (Giétro, Schwarzberg) to almost 80 km2 (Aletsch), and satellite data from the Moderate Resolution Imaging Spec- span an elevation range from about 1850 m above sea level troradiometer (MODIS) (e.g. Stroeve et al., 2013; Mernild (m a.s.l.) to over 4500 m a.s.l. (Table 1). According to the et al., 2015), downscaled MODIS data (Dumont et al., 2012; most recent Swiss Glacier Inventory, the 39 glaciers cov- Sirguey et al., 2016; Davaze et al., 2018) or other remote ered a total area of 483 km2 in 2010 (Fischer et al., 2014). sensing data sets (e.g. Wang et al., 2014) to infer trends We used a Sentinel-2 scene (20 m spatial resolution) ac- in ice albedo. Mostly, studies validated the satellite-derived quired on the 23 August 2016 to manually adjust the glacier albedo values with in situ data measured at one to several outlines and to obtain up-to-date glacier extents totalling locations on the ground, which is not always ideal (Ryan 442 km2 (Paul et al., 2016). For our analysis, we excluded et al., 2017). However, for the limited size of Alpine glaciers heavily debris-covered parts, such as medial moraines or and the complex surrounding topography, the spatial reso- debris-covered glacier tongues, as we focus on the albedo lution of MODIS data (500 m) is not suitable, and no ap- of bare ice only. Hence, the area difference of 41 km2 be- propriate high-resolution albedo product is readily available. tween 2010 and 2016 does not solely stem from glacier re- Thus, studies focusing on alpine glaciers often base their treat, but is also due to our exclusion of all glacier areas with analysis on higher-resolution data sets to obtain informa- thick debris cover for this study, i.e. debris-covered glacier tion about glacier surface albedo and related changes and tongues (e.g. Zmutt, Unteraar, Zinal, Oberaletsch) or me- processes (Dumont et al., 2011; Fugazza et al., 2016; Di dial moraines (e.g. Aletsch). The obtained glacier outlines Mauro et al., 2017) or point-measurements from automatic for the year 2016 are used in all consecutive analysis; thus, weather stations, such as the long-term monitoring site on the glacier outline is kept constant over the study period and Vadret da Morteratsch, which revealed a point-based mean does not evolve with time. summer albedo decrease between 1996 and 2006 of 0.17 We used the Landsat Surface Reflectance Level-2 science (Oerlemans et al., 2009). However, when debating the dark- products of the USGS for Landsat 5 and 7 (TM/ETM+) and 8 ening of glaciers, a clear distinction between glacier-wide vs. (OLI) as a basis to obtain broadband shortwave albedo (see point-based investigations is necessary to be able to clearly Sect. 3.2). For Landsat TM and ETM+, the product is gener- separate a darkening effect due to a changing ratio of snow- ated from the specialized software Landsat Ecosystem Dis- covered to snow-free areas of a glacier from other processes turbance Adaptive Processing System (LEDAPS), whereas affecting the reflectivity of glacier surfaces. Moreover, sepa- the Landsat OLI product is based on the Landsat 8 Surface ration between albedo changes of bare ice or snow is required Reflectance Code (LaSRC). These data products consist of to correctly distinguish between differing processes and de- six (TM/ETM+) or seven (OLI) individual spectral bands in pendencies impacting snow and ice in particular ways. the wavelength range of around 440 to 2300 nm, with slight In this study, we use Landsat data to obtain spatially dis- deviations of the individual bandwidths for the specific sen- tributed bare-ice albedo for 39 glaciers, with a total area of sors. Detailed information about these products can be found 480 km2 (corresponding to about a quarter of the present in Masek et al.(2006) for Landsat TM/ETM+, and in Ver- glaciation of the European Alps), located in western and mote et al.(2016) for Landsat 8, as well as in the product southern Switzerland over the 17-year period 1999 to 2016.