The optical properties of snow can be strongly modified by the presence of a variety of impurities including mineral dust and snow algae. We made use of measured concentration of snow algae and mineral dust to parameterize the BioSNICAR radiative transfer model. Surficial snow samples were gathered during a field campaign on 7th July 2020 at the Presena glacier (Rhaetian Alps). We collected 18 samples of surface snow containing different amount of snow algae and mineral dust. Through radiative transfer simulations we estimated an average broadband albedo reduction of 7.4 +/- 6.1 % and 35.3 +/- 7.4 % compared to clean snow, caused by snow algae and mineral dust presence, respectively. When we considered the combined effect of snow algae and dust, we estimated a broadband albedo reduction equal to 40.8 +/- 8.4 %. We estimated an average instantaneous radiative forcing induced by snow algae, mineral dust and both impurities equals to 42.3 (+/- 36.1) W/m(2), 203.7 (+/- 45.5) W/m(2), and 211.8 (+/- 45.9) W/m(2), respectively. Using BioSNICAR simulations, we also tested a series of narrowband spectral indices to determine the concentration of mineral dust and snow algae from multi- and hyper -spectral data. Results showed that most spectral indices used for snow algae mapping are correlated also with mineral dust concentration. We found that only an index correlates uniquely with snow algae: the scaled band integral at 680 nm. A new spectral index, namely the Green Blue Normalized Index, is therefore proposed to discriminate mineral dust from snow algae when both impurities are present. The high spectral resolution of current (e.g. PRISMA, EnMAP) and future (e.g. CHIME, SBG) hyperspectral satellite missions will be fundamental to decouple the effect of mineral dust and snow algae on the optical properties of snow. In fact, from those data it is possible to calculate all narrowband indices presented in this study.

For summer-accumulation-type glaciers, the glaciological literature is lacking studies on determining the snow line altitude (SLA) from optical images at the end of the summer as an indicator of the equilibrium line altitude (ELA). This paper presents a workflow for extracting the SLA from Landsat images based on the variation in the albedo with the altitude in the central line area of glaciers. The correlation of >0.8 at the 99% confidence level between the retrieved SLAs with ELAs derived from the interpolation of ground-based, mass balance measurements indicated that the workflow can be applied to derive the SLA from end-of-summer satellite data as an indicator of ELA. The ELA was under-estimated by the calculated SLA. The relationship between the end-of-summer SLA and the ELA depends on the intensity of glacier melting. Subsequently, the workflow was applied to the seven glaciers in the Eastern Tien Shan Mountains, and a time series of the SLA was obtained using 12 end-of-summer Landsat scenes from 1994 to 2016. Over the whole study period, a mean SLA of 4011.6 +/- 20.7 m above sea level (a.s.l.) was derived for the seven investigated glaciers, and an increasing SLA was demonstrated. The increase in SLAs was consistent for the seven glaciers from 1994 to 2016. Concerning the spatial variability, the east-west difference was prominent, and these differences exhibited a decreasing trend. The average SLA of each glacier is more influenced by its morpho-topographic variables. The interannual variations in the average SLA are mainly driven by the increasing summer air temperature, and the high correlation with the cumulative summer solid precipitation reflects the characteristics of the summer-accumulation-type glaciers.