Search

With the aim to reconstruct the history of atmospheric nitrogen pollution in south-eastern Europe, the ice core records of ammonium and nitrate covering the years 1774-2009 were extracted from a 182 m long ice core drilled in 2009 at Mount Elbrus in the Caucasus, Russia.The Elbrus ammonium ice-core record indicates a 4-fold increase from pre 1780 to 1990-2010. In summer, when both the vertical atmospheric mixing and anthropogenic ammonia emissions dominated by agriculture sector are strengthened, a moderate an increase of ammonium already took place from 1774 and 1880. The summer ammonium increase accelerated a first time at the beginning of the 20th century, and a second time after 1950, levels reaching a plateau between 1980 and 2010. This ice-core trend is compared to estimated past anthropogenic NH3 emissions in Europe by using state-of-the-art atmospheric transport (FLEXPART model driven with 0.5°x0.5° ERA-5 reanalysis data) simulations of submicron aerosol. It is shown that the ammonium Elbrus ice core trend is in excellent agreement with estimates of past NH3 emissions from south-eastern Europe (South european Russia, Turkey, southern Ukraine). The ammonium level of ice deposited prior to agricultural activities is found to represent ~11-15% of the 1980-2010 level. This value is in the lower range of estimates at the global scale of the contribution of natural ammonia emission to total (natural plus anthropic) ammonia emissions.A potentially existing trend of nitrate in the Elbrus ice core remained insignificant prior to ~1945. The well-marked 3-fold increase seen in the Elbrus ice core between 1950 and 1980 followed by a stabilization of levels between 1980 and 2010 is consistent with estimated past emissions of NOx from south-eastern Europe. In agreement with the partitioning at the global scale of natural and anthropogenic emissions of NOx, nitrate levels in Elbrus ice deposited prior to 1940 represent 26% of those deposited at the end of the 20th century.Study supported by the Megagrant project (agreement № 075-15-2021-599, 8.06.2021)

Computed tomography (CT) is a nondestructive high-resolution way to investigate the three-dimensional structure of samples (ice, rock, etc.). The results of CT analysis of glacial cores consisting of firn and ice extracted on the Western plateau of the Elbrus Mountain (5100–5150 m a.s.l.) in the summer of 2017 are presented in the article. The core taken from the depth of 20.31–21.87 m and consisting of three sections (average length is 52 cm each) was analyzed. In order to maintain the natural negative temperature of the glacial core, a special cryothermos has been created. It conserved the temperature at the level of −25 °C. Data on the structural features of the samples and the three-dimensional pattern of the ice-firn density were obtained. Correlations between the density and some chemical elements had been established. The CT data made it possible also to determine sizes of ice crystals. Comparison of cross sections of cores with firn and ice thin sections (30 in total) has shown that the crystal structure is best displayed in the ice inter-layers since it is impossible to determine reliably sizes of the firn grains at the given survey resolution. Also, the use of the CT method made it possible to determine inclination of the firn layers within the ice core, which is caused by the inheritance of the slope of the surface microrelief and internal inhomogeneities of the firn thickness. Calculations showed that the angle of inclination of the layers varies from 6 to 9°.

In late March 2018, a large part of the Eastern Mediterranean experienced an extraordinary episode of African dust, one of the most intense in recent years, here referred to as the "Minoan Red" event. The episode mainly affected the Greek island of Crete, where the highest aerosol concentrations over the past 15 yeas were recorded, although impacts were also felt well beyond this core area. Our study fills a gap in dust research by assessing the multi-sectoral impacts of sand and dust storms and their socioeconomic implications. Specifically, we provide a multi-sectoral impact assessment of Crete during the occurrence of this exceptional African dust event. During the day of the occurrence of the maximum dust concentration in Crete, i.e. March 22nd, 2018, we identified impacts on meteorological conditions, agriculture, transport, energy, society (including closing of schools and cancellation of social events), and emergency response systems. As a result, the event led to a 3-fold increase in daily emergency responses compare to previous days associated with urban emergencies and wildfires, a 3.5-fold increase in hospital visits and admissions for Chronic Obstructive Pulmonary Disease (COPD) exacerbations and dyspnoea, a reduction of visibility causing aircraft traffic disruptions (eleven cancellations and seven delays), and a reduction of solar energy production. We estimate the cost of direct and indirect effects of the dust episode, considering the most affected socio-economic sectors (e.g. civil protection, aviation, health and solar energy production), to be between 3.4 and 3.8 million EUR for Crete. Since such desert dust transport episodes are natural, meteorology-driven and thus to a large extent unavoidable, we argue that the efficiency of actions to mitigate dust impacts depends on the accuracy of operational dust forecasting and the implementation of relevant early warning systems for social awareness., The authors gratefully acknowledge the COST Association for funding the COST Action inDust (CA16202) as well as the WMO Sand and Dust Storm Warning Advisory and Assessment System (SDS-WAS) and the ERA4CS DustClim and the AXA Research Fund for funding the AXA Chair on Sand and Dust Storms (hosted by the Barcelona Supercomputing Center). We thank the scientific team of the PROTEAS CSP facility for the feedback provided and T. Bojic and R. Burbidge for facilitating the access to the EUROCONTROL archive. We thank the scientific team of the PROTEAS CSP facility for providing DNI data for this case study. Thanks are due to FCT/MCTES for the financial support to CESAM (UIDP/50017/2020+UIDB/50017/2020) through national funds, and also to the Icelandic Research Fund for the grant no. 207057-051. Authors S. Kazadzis and P. Kosmopoulos would like to acknowledge the European Commission project EuroGEO e-shape (grant agreement No 820852). Also, International Cooperative for Aerosol Prediction (ICAP) and NASA mission researchers are gratefully for providing aerosol data for this study. Aurelio Tobias was supported by MCIN/AEI/10.13039/501100011033 (grant CEX2018-000794-S). S. Kutuzov acknowledges the Megagrant project (agreement No. 075-15-2021-599, 8.06.2021).

Supra-glacial material, including light-absorbing impurities (LAI) such as mineral dust of crustal and soil origin, black carbon, algae and cryoconite, reduce the reflectance of snow and glacier ice. The reduction depends on the amount of LAI and their physical and chemical properties, which vary spatially and temporally. Spectral reflectance data and snow and ice samples, containing LAI, were collected in the ablation zone of the Djankuat Glacier, Central Caucasus, Russia. The spectra of the samples containing mineral dust transported from deserts were characterized by negative visible near-infrared gradients and were different from the spectra of clean aged snow and exposed glacier ice and from the samples containing mineral dust produced locally. Geochemical and mineralogical analysis using X-ray diffraction and X-ray fluorescence spectrometry showed that samples containing desert dust were characterised by a high proportion of clay materials and such minerals as smectites, illite–smectites and palygorskite and by a smaller size of mineral particles. They were enriched in chromium, zinc and vanadium. The latter served as an indicator of dust transport over or origin from the oil-producing regions of the Middle East. There was a strong negative correlation between the amount of organic matter and mineral dust in the collected samples and the albedo of surfaces from which the samples were collected. The results suggested that organic matter reduced albedo more efficiently than mineral dust. The study highlighted the importance of supra-glacial material in changing the surface reflectivity of snow and glaciers in the Caucasus region.

Despite their high value and importance for various glaciological applications, detailed ice thickness measurements of alpine glaciers are still very limited. Knowledge of bedrock topography is essential for paleoglaciological studies. The Guliya ice cap located on the Tibetan Plateau is one of the highest and largest ice caps in mid-low latitude regions. A detailed ground-penetrating radar (GPR) survey was conducted on the Guliya ice cap in 2015 using 20 and 40 MHz frequency antennas. An empirical Bayesian kriging method was used for ice thickness interpolation and uncertainty assessment. GPR measurements revealed complex basal topography of the Guliya glacier with a maximum thickness of 371.12 ± 13 m. The internal reflections caused by changes in the dielectric properties were registered on the 40 MHz radargrams at the summit and were attributed to density variations. As a result of this fieldwork, one of the largest ice thickness datasets in High Mountain Asia was obtained. Guliya glacier elevation changes were assessed by differencing digital elevation models. The glacier gained mass from 2000 to 2015 with an average rate of 0.270 ± 0.11 m w.e. a−1at the summit and 0.279 ± 0.11 m w.e. a−1at the lower elevations.

Meteorological regime and glacier surface heat balance, GPR measurements of the ice thickness and seasonal snow cover were investigated in the crater of the Eastern Summit of Mt. Elbrus In the period from 18 to 30 August 2020 at 560 m a.s.l. On the base of preliminary data analysis, the predominance of fluctuations in the synoptic scale over the diurnal ones was revealed; high values of the average and maximum wind speed associated with the influence of jet currents and with the effects of leeward storms were identified; extremely high temporal variability of relative humidity and its very high deficit in cloudless conditions, which contributes to intensive evaporation and sublimation from the snow surface, were explored. The maximum ice thickness in the crater reaches 100 m, with an average of 45 m. A new 96.01 m ice core from glacier surface to bedrock has been recovered. The drilling speed varied from 11 to 1 m / h, decreasing on average with depth from 4.5 to 4.0 m / h. The thickness of the snow-firn strata is about 20 m, which is three times less than on the Western Plateau. The borehole temperature was measured. The temperature on the glacier bedrock was -0.6 °С. The calculated heat flux was 0.39 W/m2. Air sampling was carried out in the crater of the Eastern Summit of Elbrus and on the Garabashi glacier. The repeated measurement of the soil temperature in the fumarole field on the Elbrus Eastern Summit outer crater rim suggests that the temperature regime is stable.The research was carried out on the territory of the Elbrus National Park with the financial support of the Russian Science Foundation (project 17-17-01270).

The evolution of the Elbrus glacier complex, consisting of two dozen of glaciers, in the last two decades of the 20th century and at the beginning of the 21st century generally corresponded to the trend of a decrease in the glaciated area of ​​the whole Caucasus. Over the period 1960-2014, the area of ​​Elbrus glaciation decreased by approximately 15%, and over two decades 1997-2017 - by almost 11%. As of 2017, the area of ​​Elbrus glaciation was estimated to ca. 112 sq. km, its volume exceeded 5 cub. km. Elbrus glaciation contributes significantly to the formation of the hydrological regime in the region, and, therefore, may be considered as a major challenge ti the regional socio-economic development. The latter circumstance requires an accurate assessment of the glacial runoff, and, consequently, the calculation of the surface mass balance of the glacial complex. We use an energy balance model to calculate the current and future surface mass balance. The series of observations at the Terskol meteorological station, located fifteen kilometers from the southern spurs of Elbrus, and the Mestia meteorological station, located somewhat further, on the territory of Georgia on the southern slope of the Main Caucasian ridge, as well as data from automatic weather stations on Elbrus slopes and on Djankuat glacier a few tens of kilometers from Elbrus, were applied for model forcing to reproduce present surface mass balance. The modeling results were validated by comparison with the measured surfave mass balance components on Garabashi glacier, one of the glaciers on the southern slope of Elbrus. Climate projections until the end of the 21st century for the Elbrus region were composed on the basis of multi-model results of regional climate modeling within the CORDEX project for various scenarios.We demonstrate that simultaneous surface air temperature and insolation growth accompanied by decrease in precipitation, predicted by multi-model regional climate modeling and downscaled to the Central Caucasus area, will cause essential lifting of the equilibrium line altitude and shrinking of accumulation area. As a result, we must expect an accelerated degradation of Elbrus glaciation in forthcoming decades. The reported study was funded by RFBR and RS, project number 21-55-100003

The report highlights the results of first ice-core palynology studies from the Elbrus Western Plateau. The title of the highest point in Europe and the geographical location of Elbrus determine the diversity of natural conditions and, as a result, palynological spectra, which act as markers of seasonal vegetation, climate dynamics, fires and anthropogenic activities in the Mediterranean, southern European Russia, the Middle East, and North Africa.The 24-m ice core from the Elbrus Western Plateau collected in 2017 (5115 m a.s.l., 43о20′53,9′′ N, 42о25′36′′ E) covers the period 2012-2017. Pollen analysis revealed a significant number of biological markers contained in the ice core, including pollen and spores, fungi, algae, testate amoebae, feather barbules, microcharcoal, and black carbon.The obtained results show that taxonomic diversity and concentration of biomarkers in the ice core were determined by the seasons of the year and their inherent convective flows. Pollen assemblages are characterized by predominance of native Caucasian plant species. Among them pollen values of Picea forming the high-altitude forest belt in the Western Caucasus significantly exceed pollen frequency of Pinus growing near the upper timber line on Elbrus Mt in the Central Caucasus that suggests a westerlies of air masses and transfer of microparticles. A high abundance of non-pollen palynomorphs in pollen assemblages demonstrates a high potential for studying of human impact on mountain ecosystems. The first pollen data from the ice core evidences a promising resource of the high-altitude temperate glaciers as a flexible tool for atmospheric monitoring of microparticle transfer and fixing its seasonality and biotic relationships.This work was supported by the Russian Science Foundation, project № 17-17-01270.

Glacier mass balance measurements, reconstructions and modeling are the precondition for assessing glacier sensitivity to regional climatic fluctuations. This paper presents new glaciological and geodetic mass balance data of Austre Grønfjordbreen located in the western part of Nordenskiöld Land in Central Spitsbergen. The average annual mass balance from 2014 to 2019 was −1.59 m w.e. The geodetic mass balance from 2008 to 2017 was −1.34 m w.e. The mass balance was also reconstructed by the temperature-index model from 2006 to 2020 and by spatially-distributed energy-balance models for 2011–2015 and 2019. We found a cumulative mass balance of −21.62 m w.e. over 2006–2020. The calculated mass-balance sensitivity to temperature was −1.04 m w.e. °C−1, which corresponds to the highest glacier mass balance sensitivity among Svalbard glaciers. Sensitivity to precipitation change was 0.10 m w.e. for a 10% increase in precipitation throughout the balance year. Comparing the results of the current study with other glacier mass balance assessments in Svalbard, we found that Austre Grønfjordbreen loses mass most rapidly due to its location, which is mostly influenced by the warm West Spitsbergen Current, small area and low elevation range.

Understanding the mechanisms controlling glacial retreat in the tropical Andes can strengthen future predictions of ice cover in the region. As glaciers are a dominant freshwater source in these regions, accurate ice cover predictions are necessary for developing effective strategies to protect future water resources. In this study, we investigated a 97-year dust record from two Nevado Illimani ice cores to determine the dominant factors controlling particle concentration and size distribution. In addition, we measured the area of a Nevado Illimani glacier (glacier n°8) using aerial photographs from 1956 and 2009. We identified two dustier periods during the 20th century (1930s–1940s and 1980s–2016), which were linked to reduced moisture transport from the Amazon basin. This promoted an unprecedented increase in the percentage of coarse dust particles (CPPn, ∅ > 10 μm) during the 1990s, as drier local conditions favored the emission and deposition of coarse particles on the glacier. Moisture advection from the Amazon basin to Nevado Illimani was influenced by tropical North Atlantic sea surface temperatures (TNA), which was supported by the correlation between TNA and CPPn (r = 0.52). Furthermore, glacial retreat has been accelerating since the 1980s, and a notable relationship between CPPn and the freezing level height (FLH, r = 0.41) was observed. This suggests that higher FLHs promote glacial retreat, which exposes fresh glacial sediments and facilitates the transport of coarse dust particles to the Nevado Illimani summit. Therefore, both the area of glacier n°8 and the ice core record of coarse dust particles were found to respond to climate variability—particularly to the warmer conditions across the southern tropical Andes and drier conditions over the Amazon basin.

Methane (CH4) is the third most powerful greenhouse gas. However, its warming potential is two orders of magnitude higher than of carbon dioxide and its residence time in the atmosphere is only 9.1 ± 0.9 years. It makes CH4 a good indicator of rapid climate variations both under natural conditions and due to the anthropogenic influence.The Elbrus ice core was drilled in 2009 on the Western Plato (43°20’53.9’’N, 42°25’36.0’’E) at elevation 5115 m a.s.l. It is 182 m long and is dated back to 280 ± 400 CE (Common Era). The CH4 mixing ratios were analyzed using a continuous flow analysis (CFA) system paired with optical-feedback cavity-enhanced absorption spectroscopy. The measurements campaign was organized at Institut des Géosciences de l'Environnement (IGE), Grenoble, France. This is a first high-resolution mid-latitude CH4 record. The record aims to better constrain the past evolution of mid-latitude methane sources.Here we present preliminary results of the methane concentration measurements of the Elbrus ice core in high-resolution (CFA CH4 record). We observe in situ production (max level 2900 ppb) and a baseline. We inspect a potential origin of the multiple spikes in the high-resolution record. Supposedly, either an in-situ production in the dust-rich layers occurred or a gas dissolution in the melt layers took place. However, the possibility of in-situ production during continuous gas extraction has to be further studied. The identified melt layers can serve as an indicator of interrupted stable water isotopic signal and may be supportive in the regional temperature reconstructions based on the Elbrus ice core record. A cleaned off the spikes record is inspected for the natural variability of the CH4 baseline concentration related to the short-term climate and methane emissions variability.

The ice-covered Europe's largest volcanic massif Elbrus (5,642 m) is a unique object for studying the reaction of mountain glaciers to climate changes. Elbrus glacial system contains more than 10% of the total ice volume in the Greater Caucasus. Elbrus glaciers influence on the recreation development. The rivers runoff from the Elbrus glaciers irrigates agricultural lands on steppe plains of the North Caucasus.The rate of glacier reduction in the late XX - early XXI centuries has increased significantly and in 1997-2017 Elbrus have lost 23% of its volume. Despite a number of glacier studies the mechanisms and quantitative characteristics surface mass exchange on Elbrus are still uncertain. Mass balance calculations were based on limited data. In particular, amount and distribution of snow accumulation, mass balance sensitivity to meteorological parameters under dramatic climate changes and other parameters remained unknown.Here we present the results of the detailed analysis of Garabashi glacier mass changes in 1982-2019 using glaciological and geodetic methods. Based on the new data of snow and ablation distribution the mass balance measurement system of Garabashi glacier was improved in 2018-2019. The mass balance over the studied period was also modelled using both temperature-index and distributed energy mass balance models calibrated by in situ measurements and albedo estimates from the remote sensing.The mass balance of the Garabashi glacier was close to zero or slightly positive in 1982-1997 and the cumulative mass balance was 1 m w.e. in this period. In 1997-2017 Garabashi glacier lost 12.58 m w.e. and 12.92 ± 0.95 m w.e. (−0.63 and −0.65 ± 0.05 m w.e. a−1) estimated by glaciological and geodetic method, respectively. Additional -1.7 m w.e. were lost in 2018-2019. This resulted in an area reduction by 14% and a loss of 27% of glacier volume. The observed glacier recession is driven by the pronounced increase in summer temperatures, especially since 1995, which is accompanied by nearly consistent precipitation rates The increase in incoming shortwave radiation, also played a significant role in the accelerated mass loss of glaciers in Caucasus. This study was supported by the RFBR grant 18-05-00838 a

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.

International audience

The Great Acceleration of the anthropogenic impact on the Earth system is marked by the ubiquitous distribution of anthropogenic materials throughout the global environment, including technofossils, radionuclides and the exponential increases of methane and carbon dioxide concentrations. However, personal care products as direct tracers of human domestic habits are often overlooked. Here, we present the first research combining fragrances, as novel personal care products, and polycyclic aromatic hydrocarbons (PAHs) as combustion and industrial markers, across the onset of the Great Acceleration in the Elbrus, Caucasus, ice core. This archive extends from the 1930s to 2005, spanning the profound changes in the relationship between humans and the environment during the twentieth century. Concentrations of both fragrances and PAHs rose throughout the considered period, reflecting the development of the Anthropocene. However, within this rising trend, remarkable decreases of the tracers track the major socioeconomic crises that occurred in Eastern Europe during the second half of the twentieth century.

International audience; Light absorbing particles, such as mineral dust, are a potent climate forcing agent. Many snow-covered areas are subject to dust outbreak events originating from desert regions able to significantly decrease snow albedo. Here, using a combination of Sentinel-2 imagery, in situ measurements and ensemble detailed snowpack simulations, we study the impact on snow cover duration of a major dust deposition event that occurred in the Caucasus in March 2018. This is, to the best our knowledge, the first study using ensemble approach and Sentinel-2 imagery to quantify the impact of a dust event on the snow cover evolution. We demonstrate that the calculation of the impact is strongly affected by the snow model uncertainties but that the March 2018 dust event systematically shortened the snow cover duration in Western Caucasus. The shortening is higher for location with higher accumulation and higher elevation (median values of 23+-7 days) than for location at lower elevation (median values of 15+-3 days). This is because for sites with higher location and higher accumulation, melt occurs later in the season when more incoming solar energy is available. This highlights the huge impact of a single one-day event on snow cover duration, and consequently, on the hydrology of a large region.

Rock glaciers and large ice-debris complexes are common in many mountain ranges and are especially prominent in semi-arid mountains such as the Andes or the Tien Shan. These features contain a significant amount of ice but their occurrence and evolution are not well known. Here, we present an inventory of the ice-debris complexes for the Ak-Shiirak, Tien Shan’s second largest glacierised massif, and a holistic methodology to investigate two characteristic and large ice-debris complexes in detail based on field investigations and remote sensing analysis using Sentinel-1 SAR data, 1964 Corona and recent high resolution stereo images. Overall, we found 74 rock glaciers and ice-debris complexes covering an area of 11.2 km2 (3.2% of the glacier coverage) with a mean elevation of about 3950m asl. Most of the complexes are located south-east of the main ridge of Ak-Shiirak. Ground penetrating radar (GPR) measurements reveal high ice content with the occurrence of massif debris-covered dead-ice bodies in the parts within the Little Ice Age glacier extent. These parts showed significant surface lowering, in some places exceeding 20m between 1964 and 2015. The periglacial parts are characterised by complex rock glaciers of different ages. These rock glaciers could be remnants of debris-covered ice located in permafrost conditions. They show stable surface elevations with no or only very low surface movement. However, the characteristics of the fronts of most rock glacier parts indicate slight activity and elevation gains at the fronts slight advances. GPR data indicated less ice content and slanting layers which coincide with the ridges and furrows and could mainly be formed by glacier advances under permafrost conditions. Overall, the ice content is decreasing from the upper to the lower part of the ice-debris complexes. Hence, these complexes, and especially the glacier-affected parts, should be considered when assessing the hydrological impacts of climate change.

Summary Comparison of two methods of measurements of snow cover thickness on the glacier Austre Grønfjordbreen, Svalbard was performed in the spring of 2014. These methods were the radar (500 MHz) observations and standard snow surveys. Measurements were conducted in 77 different points on the surface of the glacier. A good correlation (R2 = 0.98) was revealed. In comparison with the data of snow surveys, the radar measurements show a similar but more detailed pattern of the distribution of the snow cover depth. The discrepancy between the depths of snow cover on maps plotted from data of both methods did not exceed 30 cm in most parts of the glacier. The standard error of interpolation of the radar data onto the entire glacier surface amounts, on average, to 18 cm. This corresponds to the error of radar measurements of 18.8% when an average snow depth is about 160 cm and 9.4% at its maximum thickness of 320 cm. The distance between the measurement points at which the spatial covariance of the snow depth disappears falls between 236 and 283 m along the glacier, and between 117 and 165 m across its position. We compared the results of radar measurements of the pulse-delay time of reflections from the base of the snow cover with the data of manual probe measurements at 10 points and direct measurements of snow depth and average density in 12 snow pits. The average speed of radio waves propagation in the snow was determined as Vcr = 23.4±0.2 cm ns−1. This magnitude and the Looyenga and Kovacs formulas allowed estimating the average density of snow cover ρL = 353.1±13.1 kg m−3 and ρK = 337.4±12.9 kg m−3. The difference from average density measured in 12 pits ρav.meas = 387.4±12.9 kg m−3 amounts to −10.8% and −14.8%. In 2014, according to snow and radar measurements, altitudinal gradient of snow accumulation on the glacier Austre Grønfjordbreen was equal to 0.21 m/100 m, which is smaller than the average values (0.35 m/100 m). According to the results of snow measurements of 2011–2014, the average thickness of the snow cover on the glacier Austre Grønfjordbreen was by 17 cm greater than in 1979. In the very snowy year 2012, it was higher by 21.5 cm in comparison with the year 1979, and its spatial variability (standard deviation σН) had increased by 25.6 cm. Estimates of spatial and temporal variability of snow cover depth will be used to analyze the hydrothermal state of the glacier and its changes with regard to revealed features and climatic trends.

Black carbon (BC), emitted by fossil fuel combustion and biomass burning, is the second largest man-made contributor to global warming after carbon dioxide (Bond et al., 2013). However, limited information exists on its past emissions and atmospheric variability. In this study, we present the first high-resolution record of refractory BC (rBC, including mass concentration and size) reconstructed from ice cores drilled at a high-altitude eastern European site in Mt. Elbrus (ELB), Caucasus (5115 m a.s.l.). The ELB ice core record, covering the period 1825–2013, reflects the atmospheric load of rBC particles at the ELB site transported from the European continent with a larger rBC input from sources located in the eastern part of Europe. In the first half of the 20th century, European anthropogenic emissions resulted in a 1.5-fold increase in the ice core rBC mass concentrations with respect to its level in the preindustrial era (before 1850). The summer (winter) rBC mass concentrations increased 5-fold (3.3-fold) in 1960–1980, followed by a decrease until ∼ 2000. Over the last decade, the rBC signal for summertime slightly increased. We have compared the signal with the atmospheric BC load simulated using past BC emissions (ACCMIP and MACCity inventories) and taken into account the contribution of different geographical regions to rBC distribution and deposition at the ELB site. Interestingly, the observed rBC variability in the ELB ice core record since the 1960s is not in perfect agreement with the simulated atmospheric BC load. Similar features between the ice core rBC record and the best scenarios for the atmospheric BC load support anthropogenic BC increase in the 20th century being reflected in the ELB ice core record. However, the peak in BC mass concentration observed in ∼ 1970 in the ice core is estimated to occur a decade later from past inventories. BC emission inventories for the period 1960s–1970s may be underestimating European anthropogenic emissions. Furthermore, for summertime snow layers of the 2000s, the slightly increasing trend of rBC deposition likely reflects recent changes in anthropogenic and biomass burning BC emissions in the eastern part of Europe. Our study highlights that the past changes in BC emissions of eastern Europe need to be considered in assessing ongoing air quality regulation.

This study describes the analysis of changes in area and volume of the Mt.Elbrus glacier system, Central Caucasus from 1997 to 2017. It is based on helicopter-borne ice thickness measurements, comparison of high-resolution imagery and two digital elevation models (DEMs) with 10 m resolution. More than 250 km of ground-penetrating radar (GPR) profiles of ice thickness with reliable reflections were obtained. The total volume of Mt. Elbrus glaciers was 5.03 ± 0.85 km3 of ice in 2017. Our results show that 68% of the total ice volume is concentrated below 4,000 m a.s.l. where the average ice thickness was 44.6 ± 7.3 m, 18% of the volume lies within 4,000–4,500 m a.s.l. (thickness of 41.2 ± 7.3 m), and just 14% lies above 4,500 m a.s.l. (thickness of 29.7 ± 6.7 m). The glacier-covered area of Mt. Elbrus decreased from 125.76 ± 0.65 km2 in 1997 to 112.20 ± 0.58 km2 in 2017, a reduction of 10.8%. Over the same period the volume decreased by 22.8%. The mass balance of the Elbrus glaciers decreased by −0.55 ± 0.04 m w.e. a−1 from 1997 to 2017. Mass balance on west-oriented glaciers is less negative than on east-and south-oriented glaciers where mass balance is most negative. The mass balance of the east-oriented Djikiugankez glacier decreased at the fastest average rate (−0.97 ± 0.07 m w.e. a−1). This glacier contains 28% of the total Elbrus glacier system ice volume, most of which is concentrated below 4,000 m a.s.l. Only one small glacier on the western slope demonstrated mass gain. Our results match well with the long term direct mass balance measurements on the Garabashi glacier on Elbrus which lost 12.58 m w.e. and 12.92 ± 0.95 m w.e. between 1997 and 2017 estimated by glaciological and geodetic method, respectively. The rate of Elbrus glacier mass loss tripled in 1997-2017 compared with the 1957-1997 period.

This study reports on the glaciochemistry of a deep ice core (182 m long) drilled in 2009 at Mount Elbrus (43°21′ N, 42°26′ E; 5115 m above sea level) in the Caucasus, Russia. Radiocarbon dating of the particulate organic carbon fraction in the ice suggests a basal ice age of ~ 1670 ± 400 cal yr BP. Based on chemical stratigraphy, the upper 168.6 m of the core were dated by counting annual layers. The seasonally resolved chemical records cover the years 1774–2009 (Common Era), thus, being useful to reconstruct many aspects of atmospheric pollution in central Europe from pre-industrial times to present-day. After having examined the extent to which the arrival of large dust plumes originating from Sahara and Middle East modifies the chemical composition of the Elbrus (ELB) snow and ice layers, we focus on the sulfur pollution. The ELB sulfate levels indicate a four- and six-fold increase from 1774–1900 to 1980–1995 in winter and summer, respectively. Remaining close to 116 ± 28 ppb during the nineteen century, the summer sulfate levels started to rise at a mean rate of ~ 6 ppb per year from 1920 to 1950. The summer sulfate increase accelerated between 1950 and 1975 (11 ppb per year), levels reaching a maximum between 1980 and 1990 (730 ± 152 ppb) and subsequently decreasing to 630 ± 130 ppb at the beginning of the twenty first century. Long-term sulfate trends observed in the ELB ice cores are compared with those previously obtained in Alpine ice, the most important difference consists in a more pronounced decrease of the sulfur pollution over the three last decades in western than central Europe.

Ice cores are one of the most valuable paleo-archives. Records from the ice cores can provide information not only about the amount of dust in the atmosphere but also about dust sources and its changes in the past. A 182 m long ice core has been recovered at the western plateau of Mt. Elbrus (5115 m elevation) in 2009. This record was extended with the shallow ice core drilling in 2013. Here we present analysis of the concentrations of Ca2+, a commonly used proxy of dust, recorded in Elbrus ice core over the period 1774–2013. The calcium record reveals a quasi decadal variability with a general increasing trend. Using multiple regression analysis we found a statistically significant spatial correlation of the Elbrus Ca2+ summer concentrations and precipitation and soil moisture content in Levant region (specifically Syria and Iraq). The Ca2+ record also correlates with drought index in North Africa (r = 0.69 p r = 0.71 p 2+ and Pacific circulation indexes (PDO, SOI and Niño 4) which may indicate that due to increase in frequency of extreme El Niño and La Niña events with climate warming the influence of circulation over the Pacific Ocean extended to the Middle East.

Glaciers distinct from the Greenland and Antarctic ice sheets cover an area of approximately 706,000 square kilometres globally1, with an estimated total volume of 170,000 cubic kilometres, or 0.4 metres of potential sea-level-rise equivalent2. Retreating and thinning glaciers are icons of climate change3 and affect regional runoff4 as well as global sea level5,6. In past reports from the Intergovernmental Panel on Climate Change, estimates of changes in glacier mass were based on the multiplication of averaged or interpolated results from available observations of a few hundred glaciers by defined regional glacier areas7–10. For data-scarce regions, these results had to be complemented with estimates based on satellite altimetry and gravimetry11. These past approaches were challenged by the small number and heterogeneous spatiotemporal distribution of in situ measurement series and their often unknown ability to represent their respective mountain ranges, as well as by the spatial limitations of satellite altimetry (for which only point data are available) and gravimetry (with its coarse resolution). Here we use an extrapolation of glaciological and geodetic observations to show that glaciers contributed 27 ± 22 millimetres to global mean sea-level rise from 1961 to 2016. Regional specific-mass-change rates for 2006–2016 range from −0.1 metres to −1.2 metres of water equivalent per year, resulting in a global sea-level contribution of 335 ± 144 gigatonnes, or 0.92 ± 0.39 millimetres, per year. Although statistical uncertainty ranges overlap, our conclusions suggest that glacier mass loss may be larger than previously reported11. The present glacier mass loss is equivalent to the sea-level contribution of the Greenland Ice Sheet12, clearly exceeds the loss from the Antarctic Ice Sheet13, and accounts for 25 to 30 per cent of the total observed sea-level rise14. Present mass-loss rates indicate that glaciers could almost disappear in some mountain ranges in this century, while heavily glacierized regions will continue to contribute to sea-level rise beyond 2100. The largest collection so far of glaciological and geodetic observations suggests that glaciers contributed about 27 millimetres to sea-level rise from 1961 to 2016, at rates of ice loss that could see the disappearance of many glaciers this century.

Debris cover on glaciers can significantly alter melt, and hence, glacier mass balance and runoff. Debris coverage typically increases with shrinking glaciers. Here, we present data on debris cover and its changes for 559 glaciers located in different regions of the Greater Caucasus mountains based on 1986, 2000 and 2014 Landsat and SPOT images. Over this time period, the total glacier area decreased from 691.5 km2 to 590.0 km2 (0.52 % yr−1). Thereby, the debris covered area increased from ~ 11 to ~ 24 % on the northern, and from ~ 4 to 10 % on the southern macro-slope between 1986 and 2014. Overall, we found 18 % debris cover for the year 2014. With the glacier shrinkage, debris-covered area and the number of debris-covered glaciers increased as a function of elevation, slope, aspect, glacier morphological type, Little Ice Age moraines, and lithology.

A 182 m ice core has been recovered from a borehole drilled through the glacier to the bedrock at the Western Plateau of Mt Elbrus (43°20'53.9'' N, 42°25'36.0'' E; 5115 m a.s.l.), the Caucasus, Russia, in 2009. This is the first ice core in the region which represents a paleoclimate record practically undisturbed by seasonal melting. Relatively high snow accumulation rate at the drilling site enabled analysis of the intra-seasonal climate proxies' variability. Borehole temperatures ranged from −17 °C at 10 m depth and −2.4 °C at 182 m. A detailed radio-echo sounding survey showed that the glacier thickness ranged from 45 m near marginal zone of the plateau up to 255 m at the central part. The ice core has been analyzed for stable isotopes (δ18O and δ D), major ions (K+, Na+, Ca2+, Mg2+, NH4+, SO42-, NO3-, Cl-, F-), succinic acid (HOOCCH2COOH), and tritium content. The mean annual net accumulation rate was estimated from distinct annual oscillations of δ18O, δ D, succinic acid, and NH4+ and is 1455 mm w.e. for the last 140 years. Using annual layer counting also for the dating of the ice core, a good agreement with the absolute markers of the tritium 1963 bomb test time horizon located at the core depth of 50.7 m w.e. and the sulfate peak of the Katmai eruption (1912) at 87.7 m w.e. was obtained. According to mathematical modeling results, the bottom ice age at the maximal glacier depth is predicted to be about 660 years BP. As the 2009 borehole was situated downstream of this point, the estimated bottom ice age of the drilling site does not exceed 350–400 years BP. Taking into account the information that we have acquired on the Western Plateau Elbrus glacier and first results of the ice core analysis, these data can be used to reconstruct the atmospheric history of the European region.

Data collected between 1974 and 2016 from snow pits and core samples from two Peruvian ice fields demonstrate the effect of the recent warming over the tropical Andes, augmented by El Nino, on the preservation of the climate record. As the 0°C isotherm is approaching the summit of the Quelccaya ice cap in the Andes of southern Peru (5,670 meters above sea level (masl)), the distinctive seasonal δ18O oscillations in the fresh snow deposited within each thermal year are attenuated at depth due to melting and percolation through the firn. This has become increasingly pronounced over 43 years. In the Andes of northern Peru, the ice field on the col of Nevado Huascaran (6050 masl) has retained its seasonal δ18O variations at depth due to its higher elevation. During the 2015/2016 El Nino, snow on Quelccaya and Huascaran was isotopically (δ18O) enriched and the net sum of accumulation over the previous year (NSA) was below the mean for non–El Nino years, particularly on Quelccaya (up to 64% below the mean) which was more pronounced than the NSA decrease during the comparable 1982/1983 El Nino. Interannual large-scale oceanic and middle to upper-level atmospheric temperatures influence δ18O in precipitation on both ice fields, although the influences are variably affected by strong El Nino–Southern Oscillation events, especially on Quelccaya. The rate of ice wastage along Quelccaya's margin was dramatically higher during 2015/2016 compared with that of the previous 15 years, suggesting that warming from future El Ninos may accelerate mass loss on Peruvian glaciers.

Changes in the map area of 498 glaciers located on the Main Caucasus ridge (MCR) and on Mt. Elbrus in the Greater Caucasus Mountains (Russia and Georgia) were assessed using multispectral ASTER and panchromatic Landsat imagery with 15 m spatial resolution in 1999/2001 and 2010/2012. Changes in recession rates of glacier snouts between 1987–2001 and 2001–2010 were investigated using aerial photography and ASTER imagery for a sub-sample of 44 glaciers. In total, glacier area decreased by 4.7 ± 2.1% or 19.2 ± 8.7 km2 from 407.3 ± 5.4 km2 to 388.1 ± 5.2 km2. Glaciers located in the central and western MCR lost 13.4 ± 7.3 km2 (4.7 ± 2.5%) in total or 8.5 km2 (5.0 ± 2.4%) and 4.9 km2 (4.1 ± 2.7%) respectively. Glaciers on Mt. Elbrus, although located at higher elevations, lost 5.8 ± 1.4 km2 (4.9 ± 1.2%) of their total area. The recession rates of valley glacier termini increased between 1987–2000/01 and 2000/01–2010 (2000 for the western MCR and 2001 for the central MCR and Mt.~Elbrus) from 3.8 ± 0.8, 3.2 ± 0.9 and 8.3 ± 0.8 m yr−1 to 11.9 ± 1.1, 8.7 ± 1.1 and 14.1 ± 1.1 m yr−1 in the central and western MCR and on Mt. Elbrus respectively. The highest rate of increase in glacier termini retreat was registered on the southern slope of the central MCR where it has tripled. A positive trend in summer temperatures forced glacier recession, and strong positive temperature anomalies in 1998, 2006, and 2010 contributed to the enhanced loss of ice. An increase in accumulation season precipitation observed in the northern MCR since the mid-1980s has not compensated for the effects of summer warming while the negative precipitation anomalies, observed on the southern slope of the central MCR in the 1990s, resulted in stronger glacier wastage.

A 181.8 m ice core was recovered from a borehole drilled into bedrock on the western plateau of Mt El'brus (43°20′53.9′′ N, 42°25′36.0′′ E; 5115 m a.s.l.) in the Caucasus, Russia, in 2009 (Mikhalenko et al., 2015). Here, we report on the results of the water stable isotope composition from this ice core with additional data from the shallow cores. The distinct seasonal cycle of the isotopic composition allows dating by annual layer counting. Dating has been performed for the upper 126 m of the deep core combined with 20 m from the shallow cores. The whole record covers 100 years, from 2013 back to 1914. Due to the high accumulation rate (1380 mm w.e. year−1) and limited melting, we obtained isotopic composition and accumulation rate records with seasonal resolution. These values were compared with available meteorological data from 13 weather stations in the region and also with atmosphere circulation indices, back-trajectory calculations, and Global Network of Isotopes in Precipitation (GNIP) data in order to decipher the drivers of accumulation and ice core isotopic composition in the Caucasus region. In the warm season (May–October) the isotopic composition depends on local temperatures, but the correlation is not persistent over time, while in the cold season (November–April), atmospheric circulation is the predominant driver of the ice core's isotopic composition. The snow accumulation rate correlates well with the precipitation rate in the region all year round, which made it possible to reconstruct and expand the precipitation record at the Caucasus highlands from 1914 until 1966, when reliable meteorological observations of precipitation at high elevation began.

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 in the order of 10 ± 24% of the mean ice thickness (1-sigma estimate). For models relying on multiple data sets -- such as surface ice velocity fields, surface mass balance, or rates of ice thickness change – the results highlighted the sensitivity to input data consistency. 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.

Surge-Type glaciers have been observed in several mountain ranges of the world. Though Karakoram and Pamir are the hot spots for the occurrence of surge-Type glaciers in High Mountain Asia, few surge-Type glaciers also exist in Tien Shan. These have not been studied or reported in detail in the recent literature. We have identified 39 surge-Type glaciers and five tributary surges in Tien Shan either from available literature or by visual interpretation using available images from the period 1960 until 2014. Out of the 39 glaciers, 9 are confirmed as surge-Type, 13 are very probably surge-Type, and the remaining are possibly of surge-Type. Most of the surge-Type glaciers are located in Ak-Shiirak and Central Tien Shan. Compared with the normal glaciers of Tien Shan, the surge-Type glaciers are larger, cover higher ranges of elevations, and have shallower slopes. There is no significant difference in aspect. The largest surge events were observed in Central Tien Shan: North Inylchek Glacier (years 1996/1997) and Samoilowich Glacier (years 1992 until 2006) advanced several kilometers. The surge cycle was around 50 years for both of these glaciers. The advance was less pronounced for all other surge-Type glaciers during the period ca. 1960-2014. Some of the tributary glaciers behaved differently than the main glaciers in the sense that they continuously advanced during the entire period of our study, whereas the main glaciers have remained stable or retreated. Publisher PDF

First-ever ice core drilling at Mt. Kazbek (Caucasus Mountains) took place in the summer of 2014. A shallow ice core (18 m) was extracted from a plateau at ~4500 m a.s.l. in the vicinity of the Mt. Kazbek summit (5033 m a.s.l.). A detailed radar survey showed that the maximum ice thickness at this location is ~250 m. Borehole temperature of −7 °C was measured at 10 m depth. The ice core was analyzed for oxygen and deuterium isotopes and dust concentration. From the observed seasonal cycle, it was determined that the ice core covers the time interval of 2009–2014, with a mean annual snow accumulation rate of 1800 mm w. eq. Multiple melt layers have been detected. δ18O values vary from −25 to −5‰. The dust content was determined using a particle sizing and counting analyzer. The dust layers were investigated using scanning electron microscopy and X-ray diffraction analysis. Dust can be separated into two categories by its origin: local and distant. Samples reflecting predominantly local origin consisted mainly of magmatic rocks, while clay minerals were a characteristic of dust carried over large distances, from the deserts of the Middle East and Sahara. The calculated average dust flux over three years at Kazbek was of 1.3 mg/cm2 a−1. Neither δ18O nor dust records appear to have been affected by summer melting. Overall, the conditions on Kazbek plateau and the available data suggest that the area offers good prospects of future deep drilling in order to obtain a unique environmental record.

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-green-blue 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.

A significant desert dust deposition event occurred on Mt. Elbrus, Caucasus Mountains, Russia on 5 May 2009, where the deposited dust later appeared as a brown layer in the snow pack. An examination of dust transportation history and analysis of chemical and physical properties of the deposited dust were used to develop a new approach for high-resolution "provenancing" of dust deposition events recorded in snow pack using multiple independent techniques. A combination of SEVIRI red-green-blue composite imagery, MODIS atmospheric optical depth fields derived using the Deep Blue algorithm, air mass trajectories derived with HYSPLIT model and analysis of meteorological data enabled identification of dust source regions with high temporal (hours) and spatial (ca. 100 km) resolution. Dust, deposited on 5 May 2009, originated in the foothills of the Djebel Akhdar in eastern Libya where dust sources were activated by the intrusion of cold air from the Mediterranean Sea and Saharan low pressure system and transported to the Caucasus along the eastern Mediterranean coast, Syria and Turkey. Particles with an average diameter below 8 μm accounted for 90% of the measured particles in the sample with a mean of 3.58 μm, median 2.48 μm. The chemical signature of this long-travelled dust was significantly different from the locally-produced dust and close to that of soils collected in a palaeolake in the source region, in concentrations of hematite. Potential addition of dust from a secondary source in northern Mesopotamia introduced uncertainty in the "provenancing" of dust from this event. Nevertheless, the approach adopted here enables other dust horizons in the snowpack to be linked to specific dust transport events recorded in remote sensing and meteorological data archives.

Black carbon (BC), emitted by fossil fuel combustion and biomass burning, is the second largest man-made contributor to global warming after carbon dioxide (Bond et al., 2013). However, limited information exists on its past emissions and atmospheric variability. In this study, we present the first high-resolution record of refractory BC (rBC, including mass concentration and size) reconstructed from ice cores drilled at a high-altitude Eastern European site in Mt. Elbrus (ELB), Caucasus (5115 m a.s.l.). The ELB ice core record, covering the period 1825–2013, reflects the atmospheric load of rBC particles at the ELB site transported from the European continent with a larger rBC input from sources located in the Eastern part of Europe. In the first half of the 20th century, European anthropogenic emissions resulted in a 1.5-fold increase in the ice core rBC mass concentrations as respect to its level in the preindustrial era (before 1850). The rBC mass concentrations increased by a 5-fold in 1960–1980, followed by a decrease until ~ 2000. Over the last decade, the rBC signal for summer time slightly increased. We have compared the signal with the atmospheric BC load simulated using past BC emissions (ACCMIP and MACCity inventories) and taken into account the contribution of different geographical region to rBC distribution and deposition at the ELB site. Interestingly, the observed rBC variability in the ELB ice core record since the 1960s is not in perfect agreement with the simulated atmospheric BC load. Similar features between the ice core rBC record and the best scenarios for the atmospheric BC load support that anthropogenic BC increase in the 20th century is reflected in the ELB ice core record. However, the peak in BC mass concentration observed in ~ 1970 in the ice core is estimated to occur a decade later from past inventories. BC emission inventories for the period 1960s–1970s may be underestimating European anthropogenic emissions. Furthermore, for summer time snow layers of the last 2000s, the slightly increasing trend of rBC deposition likely reflects recent changes in anthropogenic and biomass burning BC emissions in the Eastern part of Europe. Our study highlights that the past changes in BC emissions of Eastern Europe need to be considered in assessing on-going air quality regulation.

The observed increase in summer temperatures and the related glacier downwasting has led to a noticeable decrease of frozen water resources in Central Asia with possible future impacts on the economy of all downstream countries in the region. Glaciers in the Ak Shyirak massif located in the Inner Tien Shan are not only affected by climate change but also impacted by the open pit gold mining of the Kumtor Gold Company. In this study glacier inventories referring to the years 2003 and 2013 were created for the Ak Shyirak massif based on satellite imagery. The 193 glaciers had a total area of 351.2 ± 5.6 km2 in 2013. Compared to 2003 the total glacier area decreased by 5.9 ± 3.4. During 2003–2013 the shrinkage rate of Ak Shyirak glaciers was twice than that in 1977–2003 and similar to shrinkage rates in Tien Shan frontier ranges. We assessed glacier volume in 2013 using volume–area (VA) scaling and GlabTop modelling approaches. Resulting values for the whole massif differ strongly the \{VA\} scaling derived volume is 30.0–26.4 km3 whereas the GlabTop derived volume accounts for 18.8–13.2 km3. Ice losses obtained from both approaches were compared to geodetically derived volume change. \{VA\} scaling underestimates ice losses between 1943 and 2003 whereas GlabTop reveals a good match for eight glaciers for the period 2003–2012. In comparison to radio echo soundings from three glaciers the GlabTop model reveals a systematic underestimation of glacier thickness with a mean deviation of 16. GlabTop tends to significantly underestimate ice thickness in accumulation areas but tends to overestimate ice thickness in the lowermost parts of glacier snouts. Direct technogenic impact is responsible for about 7 of area and 5 of mass loss for glaciers in the Ak Shyirak massif during 2003–2013. Therefore the increase of summer temperature seems to be the main driver of accelerated glacier shrinkage in the area.

Changes in the extent of glaciers and rates of glacier termini retreat in the eastern Terskey-Alatoo Range, the Tien Shan Mountains, Central Asia have been evaluated using the remote sensing techniques. Changes in the extent of 335 glaciers between the end of the Little Ice Age (LIA; mid-19th century), 1990 and 2003 have been estimated through the delineation of glacier outlines and the LIA moraine positions on the Landsat TM and ASTER imagery for 1990 and 2003 respectively. By 2003, the glacier surface area had decreased by 19% of the LIA value, which constitutes a 76 km(2) reduction in glacier surface area. Mapping of 109 glaciers using the 1965 1:25,000 maps revealed that glacier surface area decreased by 12.6% of the 1965 value between 1965 and 2003. Detailed mapping of 10 glaciers using historical maps and aerial photographs from the 1943-1977 period, has enabled glacier extent variations over the 20th century to be identified with a higher temporal resolution. Glacial retreat was slow in the early 20th century but increased considerably between 1943 and 1956 and then again after 1977. The post-1990 period has been marked by the most rapid glacier retreat since the end of the LIA. The observed changes in the extent of glaciers are in line with the observed climatic warming. The regional weather stations have revealed a strong climatic warming during the ablation season since the 1950s at a rate of 0.02-0.03 degrees Ca-1. At the higher elevations in the study area represented by the Tien Shan meteorological station, the summer warming was accompanied by negative anomalies in annual precipitation in the 1990s enhancing glacier retreat. However, trends in precipitation in the post-1997 period cannot be evaluated due to the change in observational practices at this station. Neither station in the study area exhibits significant long-term trends in precipitation. Crown Copyright (C) 2009 Published by Elsevier B.V. All rights reserved.