Abstract. Knowledge of the ice thickness distribution of glaciers and ice caps is an important prerequisite for many glaciological and hydrological investigations. A wealth of approaches has recently been presented for inferring ice thickness from characteristics of the surface. With the Ice Thickness Models Intercomparison eXperiment (ITMIX) we performed the first coordinated assessment quantifying individual model performance. A set of 17 different models showed that individual ice thickness estimates can differ considerably – locally by a spread comparable to the observed thickness. Averaging the results of multiple models, however, significantly improved the results: on average over the 21 considered test cases, comparison against direct ice thickness measurements revealed deviations on the order of 10 ± 24 % of the mean ice thickness (1σ estimate). Models relying on multiple data sets – such as surface ice velocity fields, surface mass balance, or rates of ice thickness change – showed high sensitivity to input data quality. Together with the requirement of being able to handle large regions in an automated fashion, the capacity of better accounting for uncertainties in the input data will be a key for an improved next generation of ice thickness estimation approaches.

Abstract. Knowledge of supra-glacial debris cover and its changes remain incomplete in the Greater Caucasus, in spite of recent glacier studies. Here we present data of supra-glacial debris cover for 659 glaciers across the Greater Caucasus based on Landsat and SPOT images from the years 1986, 2000 and 2014. We combined semi-automated methods for mapping the clean ice with manual digitization of debris-covered glacier parts and calculated supra-glacial debris-covered area as the residual between these two maps. The accuracy of the results was assessed by using high-resolution Google Earth imagery and GPS data for selected glaciers. From 1986 to 2014, the total glacier area decreased from 691.5±29.0 to 590.0±25.8 km2 (15.8±4.1 %, or ∼0.52 % yr−1), while the clean-ice area reduced from 643.2±25.9 to 511.0±20.9 km2 (20.1±4.0 %, or ∼0.73 % yr−1). In contrast supra-glacial debris cover increased from 7.0±6.4 %, or 48.3±3.1 km2, in 1986 to 13.4±6.2 % (∼0.22 % yr−1), or 79.0±4.9 km2, in 2014. Debris-free glaciers exhibited higher area and length reductions than debris-covered glaciers. The distribution of the supra-glacial debris cover differs between the northern and southern and between the western, central and eastern Greater Caucasus. The observed increase in supra-glacial debris cover is significantly stronger on the northern slopes. Overall, we have observed up-glacier average migration of supra-glacial debris cover from about 3015 to 3130 m a.s.l. (metres above sea level) during the investigated period.

Abstract.A 182 m ice core was recovered from a borehole drilled into bedrock on the western plateau of Mt. Elbrus (43 • 20 53.9 N, 42 • 25 36.0 E; 5115 m a.s.l.) in the Caucasus, Russia, in 2009. This is the first ice core in the region that represents a paleoclimate record that is practically undisturbed by seasonal melting. Relatively high snow accumulation rates at the drilling site enabled the analysis of the intraseasonal variability in climate proxies. Borehole temperatures ranged from −17 • C at 10 m depth to −2.4 • C at 182 m. A detailed radio-echo sounding survey showed that the glacier thickness ranged from 45 m near the marginal zone of the plateau up to 255 m at the glacier center. The ice core has been analyzed for stable isotopes (δ 18 O and δD), major ions (K + , Na + , Ca 2+ , Mg 2+ , NH

Abstract. The first record of dust deposition events on Mt. Elbrus, Caucasus Mountains derived from a snow pit and a shallow firn core is presented for the 2009-2012 period. A combination of isotopic analysis, SEVIRI red-greenblue composite imagery, MODIS atmospheric optical depth fields derived using the Deep Blue algorithm, air mass trajectories derived using the HYSPLIT model and analyses of meteorological data enabled identification of dust source regions with high temporal (hours) and spatial (ca. 20-100 km) resolution. Seventeen dust deposition events were detected; fourteen occurred in March-June, one in February and two in October. Four events originated in the Sahara, predominantly in northeastern Libya and eastern Algeria. Thirteen events originated in the Middle East, in the Syrian Desert and northern Mesopotamia, from a mixture of natural and anthropogenic sources. Dust transportation from Sahara was associated with vigorous Saharan depressions, strong surface winds in the source region and mid-tropospheric southwesterly flow with daily winds speeds of 20-30 m s −1 at 700 hPa level. Although these events were less frequent than those originating in the Middle East, they resulted in higher dust concentrations in snow. Dust transportation from the Middle East was associated with weaker depressions forming over the source region, high pressure centred over or extending towards the Caspian Sea and a weaker southerly or southeasterly flow towards the Caucasus Mountains with daily wind speeds of 12-18 m s −1 at 700 hPa level. Higher concentrations of nitrates and ammonium characterised dust from the Middle East deposited on Mt. Elbrus in 2009 indicating contribution of anthropogenic sources. The modal values of particle size distributions ranged between 1.98 µm and 4.16 µm. Most samples were characterised by modal values of 2.0-2.8 µm with an average of 2.6 µm and there was no significant difference between dust from the Sahara and the Middle East.

We present a first version of the Svalbard ice‐free topography (SVIFT1.0) using a mass conserving approach for mapping glacier ice thickness. SVIFT1.0 is informed by more than 1 million point measurements, totalling more than 8,700 km of thickness profiles. SVIFT1.0 is publicly available and represents the geometric state around the year 2010. Our estimate for the total ice volume is 6,199 km3, equivalent to 1.5‐cm sea level rise. The thickness map suggests that 13% of the glacierized area is grounded below sea level. A complementary map of error estimates comprises uncertainties in the thickness surveys as well as in other input variables. Aggregated error estimates are used to define a likely ice‐volume range of 5,200–7,300 km3. The ice front thickness of marine‐terminating glaciers is a key quantity for ice loss attribution because it controls the potential ice discharge by iceberg calving into the ocean. We find a mean ice front thickness of 135 m for the archipelago (likely range 123–158 m).