<figure><img src="image_1.png" /></figure>

HAL
open science

# Analysis of Regional Changes in Geodetic Mass Balance for All CaucasusGlaciers over the Past Two Decades

Levan G. Tielidze, Vincent Jomelli, Gennady A. Nosenko

► To cite this version:

Levan G. Tielidze, Vincent Jomelli, Gennady A. Nosenko. Analysis of Regional Changes in Geodetic Mass Balance for All CaucasusGlaciers over the Past Two Decades. Atmosphere, 2022, 13 (2), 10.3390/atmos13020256. [https://doi.org/10.3390/atmos13020256]

HAL Id: hal-03662987

https://hal.science/hal-03662987

Submitted on 11 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Analysis of Regional Changes in Geodetic Mass Balance for All Caucasus Glaciers over the Past Two Decades

Levan G. Tielidze <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,* and Vincent Jomelli <fnref n="4" /> and Gennady A. Nosenko <fnref n="5" />

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Tielidze, L.G.; Jomelli, V.; Nosenko, G.A. Analysis of Regional Changes in Geodetic Mass Balance for All Caucasus Glaciers over the Past Two Decades. *Atmosphere* **2022**, *13*, 256. https://doi.org/10.3390/atmos13020256

Academic Editor: Lin Zhao

Received: 4 January 2022

Accepted: 30 January 2022

Published: 2 February 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** Glaciers and snow in the Caucasus are major sources of runoff for populated places in many parts of this mountain region. These glaciers have shown a continuous area decrease; however, the magnitude of mass balance changes at the regional scale need to be further investigated. Here, we analyzed regional changes in surface elevation (or thickness) and geodetic mass balance for 1861 glaciers (1186.1 ± 53.3 km<fnref n="2" />) between 2000 and 2019 from recently published dataset and outlines of the Caucasus glacier inventory. We used a debris-covered glacier dataset to compare the changes between debris-free and debris-covered glaciers. We also used 30 m resolution ASTER GDEM (2011) to determine topographic details, such as aspect, slope, and elevation distribution of glaciers. Results indicate that the mean rate of glacier mass loss has accelerated from 0.42 ± 0.61 m of water equivalent per year (m w.e. a⁻<fnref n="1" />) over 2000–2010, to 0.64 ± 0.66 m w.e. a⁻<fnref n="1" /> over 2010–2019. This was 0.53 ± 0.38 m w.e. a⁻<fnref n="1" /> in 2000–2019. Mass loss rates differ between the western, central, and eastern Greater Caucasus, indicating the highest mean annual mass loss in the western section (0.65 ± 0.43 m w.e. a⁻<fnref n="1" />) in 2000–2019 and much lower in the central (0.48 ± 0.35 m w.e. a⁻<fnref n="1" />) and eastern (0.38 ± 0.37 m w.e. a⁻<fnref n="1" />) sections. No difference was found between the northern and southern slopes over the last twenty years corresponding 0.53 ± 0.38 m w.e. a⁻<fnref n="1" />. The observed decrease in mean annual geodetic mass balance is higher on debris-covered glaciers (0.66 ± 0.17 m w.e. a⁻<fnref n="1" />) than those on debris-free glaciers (0.49 ± 0.15 m w.e. a⁻<fnref n="1" />) between 2000 and 2019. Thickness change values in 2010–2019 were 1.5 times more negative (0.75 ± 0.70 m a⁻<fnref n="1" />) than those in 2000–2010 (0.50 ± 0.67 m a⁻<fnref n="1" />) in the entire region, suggesting an acceleration of ice thinning starting in 2010. A significant positive trend of May–September air temperatures at two selected meteorological stations (Terskol and Mestia) along with a negative trend of October–April precipitation might be responsible for the negative mass balances and thinning for all Caucasus glaciers over the study period. These results provide insight into the change processes of regional glaciers, which is key information to improve glaciological and hydrological projections in the Caucasus region.

**Keywords:** geodetic mass balance; glacier thickness change; Greater Caucasus

## 1. Introduction

Glacier mass balance is an important variable for understanding the response of glaciers to climate change [1] and their contribution to regional water resources [2] as well as from the perspective of changes in global sea level [3–5]. Glacier mass balance can be derived by two different methods, such as the traditional glaciological method, and the geodetic method. Glaciological mass balance is based on in-situ measurements of accumulation and ablation using snow pits and stakes [1], while the geodetic mass balance can

<footnote n="1">¹ Antarctic Research Centre, Victoria University of Wellington, P.O. Box 600, Wellington 6140, New Zealand</footnote>

<footnote n="2">² School of Geography, Environment and Earth Sciences, Victoria University of Wellington, P.O. Box 600, Wellington 6140, New Zealand</footnote>

<footnote n="3">³ School of Natural Sciences and Medicine, Ilia State University, Cholokashvili Ave 3/5, Tbilisi 0162, Georgia</footnote>

<footnote n="4">⁴ Centre for Research and Teaching in Environmental Geoscience (CEREGE), Aix-Marseille University, CNRS, IRD, College De France, INRA, 13545 Aix-en-Provence, France; jomelli@cerege.fr</footnote>

<footnote n="5">5 Department of Glaciology, Institute of Geography, Russian Academy of Sciences, 29 Staromonetniy Pereulok, 119017 Moscow, Russia; gnosenko@gmail.com</footnote>

<footnote>* Correspondence: tielidzelevan@gmail.com</footnote>
be obtained from observations of elevation change derived from digital elevation models (DEMs) computed from stereo satellite images. Ground-based measurement gives glaciological mass balance relatively high accuracy at the point of measurement, which is then interpolated to the glacier area, while the main advantage of geodetic mass balance is the wide regional coverage, although, the uncertainties of geodetic mass balance measurements depend on the satellite source and is often relatively high [6]. However, close range sensing such as airborne or terrestrial laser scanning significantly reduces the uncertainty to an order of decimeters to meter [7,8]. Significant improvements have been made in the last decades to improve the accuracy of the geodetic approach [9,10]. Furthermore, this method has become widely used in recent years due to the possibility to monitor large glacierised areas [3,4].

Meltwater of glaciers and snow provides a water source for electricity generation, ecosystems, and communities in the Greater Caucasus. These glaciers also have a positive impact on the economy by being a major tourist attraction in this region. Several national parks on both sides of the Greater Caucasus (Georgia/Russia) welcome thousands of visitors each year. On the other hand, glacier hazards are relatively common in this region, leading to major casualties [11–13]. Thus, estimation of glacier mass change is of great significance in this region. Investigation of glacier mass balance will also provide important information for future research on various socio-economic impacts such as water resource management, risk assessments, and tourism in this mountain region.

Attempts to organize in-situ regular observations of glacier mass balance have been made since the early 1960s on both southern (Georgia) and northern (Russia) slopes of the Greater Caucasus [14,15]. The duration of those measurements was short for various reasons. Long-term glaciological observation campaigns (stake method) on the Russian side of the Greater Caucasus began in 1967 on Djankuat Glacier (43°11′57″ N; 42°45′24″ E) [16], which started in 1982 on Garabashi Glacier (43°18′19″ N; 42°28′13″ E) [17]. A similar monitoring program was started in 1968 on Buba (Tbilisa) Glacier (42°44′57″ N; 43°44′45″ E) on the Georgian side of the Greater Caucasus [18]. After the dissolution of the Soviet Union in 1990s, the monitoring program was interrupted in Georgia and no glaciological field measurements were carried out until 2007, when the mass balance monitoring programs were temporarily (2007–2010) initiated on Zopkhito Glacier (42°52′54″ N; 43°25′29″ E). One-year glaciological mass balance measurement was also initiated on Chalaati Glacier (43°7′49″ N; 42°42′30″ E) in 2011 [19]. The initiation of other continuous monitoring programs was difficult mainly due to hard political and social conditions resulting in more than 30 years of gaps in glaciological mass balance data. Nowadays, continuous series of long-term observations by traditional glaciological method continue only on Djankuat and Garabashi glaciers from the Russian side of the Greater Caucasus [20].

In the last few decades, research in the Greater Caucasus focused mainly on documenting changes in glacier area [21–24], indicating a significant shrinkage in regional ice coverage. Long-term ground-based measurements from Djankuat and Garabashi glaciers also suggest predominantly negative mass balances over the last few decades from this region [17,25]. Recently, geodetic mass balances have also been monitored within the local [26,27], and global [3,4] studies. The main limitations of local studies are incomplete area coverage, while the global studies always describe the Caucasus, with the Middle East making the regional analysis difficult. The vast majority of glaciers in the Middle East are very small, with mean area of <0.1 km² [28,29] and the climate is hot and arid (dry) [30], causing a much higher rate of glacier decline [28,29]. Thus, combining these two regions causes a bias with more negative values over the Caucasus region.

In this study, we use thickness and surface mass balance data from recently published global study by Hugonnet et al. (2021) [3]. This study offers a unique opportunity to analyze the behavior of individual glaciers in the Greater Caucasus as the data of mass balance and thickness changes are available for each specific glacier. Thus, the main goal of this study is to:
(i) analyze the behavior of glaciers in the Greater Caucasus separately from the Middle East;

(ii) estimate changes in surface elevation and geodetic mass balance for all Caucasus glaciers over the past two decades according to river basins and analyze these changes on sub-regional (e.g., western vs. eastern; northern vs. southern) and individual glaciers (e.g., debris-covered vs. debris-free);

(iii) compare this data with those from local existing ground-based and geodetic mass balance measurements along with local climate observations from high-mountain meteorological stations.

## 2. Study Area

Caucasus is a complex mountain system located between the Black and Caspian seas on the border of Europe and Asia consisting of two separate mountain systems such as the Greater Caucasus and the Lesser Caucasus (Transcaucasian Highlands). The Greater Caucasus is the higher and more extensive part of this region, extending for about 1100 km from northwest to southeast. The Lesser Caucasus is located approximately 100 km to the south and characterized by relatively lower elevations and warm climate, providing unfavourable conditions for the existence of modern glaciers. Thus, all contemporary glaciers of the Caucasian region belong to the Greater Caucasus.

According to geomorphological characteristics, altitudinal difference, and climatic conditions, the Greater Caucasus divided into three parts—the Western, Central, and Eastern. The central Greater Caucasus is situated between the mountain peaks of Elbrus (5642 m) and Kazbegi (5047 m) (Figure 1). Several summits in the central Greater Caucasus exceed 5000 m above sea level (a.s.l.): Dykh-Tau (5205 m), Shkhara (5203 m), Koshtan-Tau (5151 m), Jangha (Dzhangi-Tau) (5058 m), and Pushkin Peak (5034 m). The western section is the lowest part of the Greater Caucasus, with highest summit being Dombai-Ulgen (4046 m). The Eastern section of the Greater Caucasus is the most expanded part of this mountain region, with the highest summit being Mt. Bazardüzü (4466 m a.s.l.). The division of the Greater Caucasus into northern and southern macroslopes by location relative to the main crest of the range is also often used.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Glacier extent in 2014 in the Greater Caucasus according to the main river basins. Green circles correspond to changes in glacier surface mass balance (m w.e. a⁻¹) according to individual river basins. See also Table A1 to compare mass balance data relative to all river basins.</figcaption></figure>
The mountain range of the Greater Caucasus represents an orographic barrier to the flow of moisture brought from the west and southwest by Mediterranean and Atlantic cyclones. The greatest amount of precipitation falls on the southern slope of the western section. Here, the amount of annual precipitation at an altitude of 3000 m a.s.l. reaches 4000 mm. Precipitation gradually decreases from west to east, reaching about 1000 mm on the southern slopes of eastern Greater Caucasus. The northern slope is relatively more arid (dry) and cold than the southern one, resulting in ~500–1000 mm less annual precipitation and 1–2 °C lower summer temperatures [31,32]. The average regional lapse rate is minimum in winter (−2.3 °C per 1000 m) and maximum in summer (−5.2 °C per 1000 m) [33]. Overall, snow depth may reach 5–7 m in winter in several regions of the western part of the Greater Caucasus, such as the Upper Svaneti and northern Abkhazeti regions in Georgia [34].

Glacier decrease started from the last advances of the Little Ice Age (LIA) in the northern Greater Caucasus, dated approximatively to the late 1840s [35], though it may have been slightly earlier, in the 1810s, in the southern slopes [36]. The decrease rate was relatively small in centennial scale (~0.4% yr⁻¹ in 1911–2014) [22–24], while this has increased since the 1960s (0.5% yr⁻¹) [21] and reached its maximum rates over the past two decades (1.16% yr⁻¹) [37].

About 13.4% of the surface area of Caucasus glaciers was covered by debris in 2014 [38]. Nowadays, there are more than 2200 glaciers in the Greater Caucasus, covering an area of about 1060 km². The northern slopes of the mountain range contain more glaciers than the southern slopes. About 70% of the modern glaciated area is concentrated in the central section. Fourteen glaciers have an area > 10 km² resulting in ~220 km² or ~21% of the total glacier area in 2020. The Bezingi Glacier, with an area of ~39 km², was the largest glacier mapped in the 2020 database [37].

### 3. Data and Methods

In this study, we used geodetic mass balance and ice thickness data from Hugonnet et al. (2021), [3] obtained from 30 m resolution Terra Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model (GDEM) and ASTER L1A stereo images (each covering 60 km × 60 km) between 1 January 2000 and 31 December 2019 (downloaded from https://lpdaac.usgs.gov/products/ast_l1av003/) (Last access 13 December 2021). The specifically developed statistical methods are consistently described in Hugonnet et al. (2021) [3]. Outlines from the Caucasus glacier inventory from 2014 [21] were regarded as glacier boundaries in Hugonnet et al. (2021) [3]. These outlines for most glaciers in this area are identified according to the Landsat 8 Operational Land Imager (OLI) images and downloaded from the EarthExplorer website (http://earthexplorer.usgs.gov/) (Last access 13 December 2021).

The 5-year uncertainties by Hugonnet et al., 2021 [3], at the 95% confidence level, were lower than 0.5 m yr⁻¹ for glaciers larger than 1 km² and conservative for smaller glaciers. They, thus, validated the reliability of improved uncertainty approaches for volume change estimation down to the scale of individual glaciers. A detailed description of uncertainty calculation is given in Huggonnet et al., 2021 [3]. To check the accuracy of geodetic mass balance data by Hugonnet et al. (2021) [3] relative to the existing ground-based observation from the same period, we selected the Djankuat Glacier (Figure 2a) and extracted its glaciological mass balance data from the WGMS server [20]. For more comparison with other geodetic mass balance studies, we also selected the Kyukyurtyu, Ulluchiran, Karachaul, and Bolshoy Azau glaciers from the Elbrus Massif (Figure 2b), previously studied by Kutuzov et al. (2019) [27], based on two different DEMs from 1997 (aerial photography survey) and 2017 (Pléiades). We selected only these glaciers because they have almost similar internal boundaries in Hugonnet et al. (2021) [3] and Kutuzov et al. (2019) [27] and are comparable in terms of data analysis. The rest of the glaciers from Elbrus Massif were excluded for similar comparison due to the much different
internal boundaries, although the entire Elbrus Massif was also selected as a single outline with only outer/external boundaries.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Location of Djankuat (a) and Elbrus (b) glaciers. Glaciers with blue colour from the Elbrus Massif were used for further comparison (see Chapter 4.2). SPOT image 2019 and ALOS PALSAR DEM was used as a background. See also Figure 1 for relative location of these glaciers to the Greater Caucasus.</figcaption></figure>

We used a debris-covered glacier dataset from Tielidze et al. (2020a) [38] to compare geodetic mass balance and thickness changes between debris-free and debris-covered glaciers. This dataset was compiled based on the Glacier Classification Guidance from the Global Land Ice Measurements from Space (GLIMS) initiative for remote-sensing observations [39]. According to this guideline, all glaciers with a debris coverage of more than 10% were considered to be debris-covered and otherwise to be debris-free glaciers. We selected 34 glaciers from this dataset with an area of >2.0 km².

We also used the ASTER Global Digital Elevation Model (GDEM, 17 November 2011) version 3 to determine topographic details such as aspect, slope, and elevation distribution of glaciers. The DEM was downloaded from NASA LP DAAC Collections (http://earthexplorer.usgs.gov/) (Last access 13 December 2021).

Change in the mass (M) for all Caucasus glaciers over the last 20 years was calculated based on geodetic mass balance data by Hugonnet et al. (2021) [3] and the recently published Caucasus glacier inventory by Tielidze et al. (2021) [37] based on two different satellite sources from 2000 (Landsat) and 2020 (Sentinel).

This was calculated as follows:

$$M = dm \cdot \Delta t \cdot (S1 + S2)/2$$

where dm is the average surface balance (–0.53 m w.e. a⁻¹) (meter water equivalent per year, no mentioned in the rest part of this paper) for the entire Caucasus in 2000–2019 by Hugonnet et al. (2021) [3]; Δt is the duration of the observation period (20 years); S1 and S2 are the glacier areas at the beginning (1381.5 km² in 2000) and at the end (1069.1 km² in 2019/2020) of the observation period, according to Tielidze et al. (2021) [37].

We converted an ice volume (km³) to a mass of ice (in Gt) by using the following equation:

$$\text{Mass of ice (Gt)} = \text{Volume of ice (km}^3) \times \text{Density of ice (Gt/km}^3 \text{ or } 850 \pm 60 \text{ kg/m}^3)$$

Monthly climate data (temperature and precipitation) from two different meteorological stations were also used in this study. The Terskol station is located at an elevation of 2143 m a.s.l. on the Russian side of the Greater Caucasus, while the Mestia station is located at an elevation of 1441 m a.s.l. on the Georgian side (see location on Figure 1). Observations are available for the period of 2000–2018 for Terskol (temperatures and precipitation) and 2000–2015 for the Mestia station.
# 4. Results

## 4.1. *Mass Balance and Thickness Change for the Entire Region*

Changes in annual geodetic mass balance along with cumulative mass balance in the entire Greater Caucasus over the last twenty years are shown in Figure 3 (see also Table A1). We computed mean annual mass loss of 0.53 ± 0.38 m w.e. a⁻¹ for all glaciers in the Greater Caucasus between 2000 and 2019 with a cumulative mass change of -10.07 ± 7.22 m w.e. Calculated decrease in the mass of all glaciers in the Greater Caucasus over the past 20 years amounted to the ~12.99 km³ or 11.04 ± 0.78 Gt water (~0.58 Gt per year). The mean annual geodetic mass balance between 2000 and 2010 was -0.42 ± 0.61 m w.e. a⁻¹ for the entire mountain range, while it was -0.64 ± 0.66 m w.e. a⁻¹ between 2010 and 2019 (Figure 4). The rates of mass loss vary between the western, central, and eastern Greater Caucasus. This was relatively higher in the western section with 0.65 ± 0.43 m w.e. a⁻¹ in 2000–2019 (Figures 4 and 5a) and much lower in the eastern section (0.38 ± 0.37 m w.e. a⁻¹) (Figures 4 and 5c). The colour-coded map in Figure 5 shows spatial distribution of changes of surface mass for both slopes of western, central, and eastern Greater Caucasus glaciers in 2000–2019. No difference was found between the northern and southern slopes over the study period (0.53 ± 0.38 m w.e. a⁻¹ on both slopes). Mean annual mass balance on the glaciers of Russian Caucasus was slightly higher (-0.55 ± 0.39 m w.e. a⁻¹) than those on the Georgian counterparts (-0.51 ± 0.36 m w.e. a⁻¹). Glaciers in Azerbaijan experienced the lowest mass balance (-0.23 ± 0.29 m w.e. a⁻¹) over the last twenty years (Figure 4, Table A2). Mean annual surface mass balance of the north-and north-west-oriented glaciers was more negative than those of the south and south-western-oriented glaciers (Figure 6a). Figure 6b shows the mean annual geodetic mass balance rates for glaciers of seven size classes of the entire mountain range between 2000 and 2019 relative to the mean and minimum (terminus) elevations, indicating that glaciers with a size class of 0.5–1.0 km² experienced the highest mean annual mass loss (0.60 m w.e. a⁻¹). The large size class of glaciers with an area of 5.0–10.0 km² had the lowest annual mass loss (0.41 m w.e. a⁻¹) during the investigated period. Changes of geodetic mass balance and thickness for individual glaciers larger than 2.0 km² of the Greater Caucasus in 2000–2010, 2010–2019, and 2000–2019 are shown in Table A3.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Mean annual geodetic mass balance (m w.e. a⁻¹) and cumulative mass balance (m w.e.) for all Caucasus glaciers in 2000–2019.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Boxplots based on data of the individual glaciers showing changes in glacier annual geodetic mass balance (m w.e. a⁻¹) for all Caucasus glaciers according to the regions (western, central, and eastern), slopes (northern, southern, entire) and countries in 2000–2010 (blue), 2010–2019 (brown), and 2000–2019 (grey). Box gives lower and upper quartiles (vertical stems) along with the mean (black horizontal line inside of the box) and median (X symbol inside of the box) value of glacier mass change rate. The circles correspond to the outliers.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Changes in mean annual geodetic mass balance (m w.e. $a^{-1}$) (color coded) for western (a), central (b), and eastern (c) Greater Caucasus glaciers in 2000–2019.</figcaption></figure>

Correlation between the thinning rates and mean elevation for all Caucasus glaciers in observed periods suggests that lower elevation glaciers experienced higher annual surface lowering, i.e., the thinning rate varied significantly at an altitudinal range (Figure A1). The rate of thickness changes varies between investigated periods. We observed 1.5 times more negative values of thickness changes in 2010–2019 (0.75 ± 0.70 m $a^{-1}$) than in 2000–2010 (0.50 ± 0.67 m $a^{-1}$) for the entire mountain region, suggesting an acceleration of ice thinning rates starting in 2010’s (Table A1; Figure A2). The most negative values of thickness changes were observed in Russian Caucasus (0.65 ± 0.38 m $a^{-1}$) between 2000 and 2019, while the lowest thinning rate was observed in Azerbaijan (0.27 ± 0.31 m $a^{-1}$). Georgian glaciers experienced thickness changes of −0.60 ± 0.36 m $a^{-1}$ during the same time (Table A2;
Figure A2). In terms of individual river basins, the most negative elevation and mass balance change were observed between 2000 and 2010 in Bolshaya Laba basin (north-west Greater Caucasus) with $1.57 \pm 1.10 \text{ m a}^{-1}$ and $1.34 \pm 1.20 \text{ m w.e. a}^{-1}$ loss respectively (27 glaciers with a total area of $4.4 \pm 0.2 \text{ km}^2$) (Table A1, Figure 1).

<figure><img src="image_10.png" /><figcaption>**Figure 6.** Changes in mean annual mass balance (m w.e. a⁻¹) vs. mean and minimum elevation according to the aspect (a) and seven size classes (b) for all Caucasus glaciers in 2000–2019. Vertical stems show an uncertainty estimation of the mean values.</figcaption></figure>

A comparison of 34 selected glaciers (>2.0 km²) shows that 17 debris-covered glaciers experienced higher mass loss ($0.66 \pm 0.17 \text{ m w.e. a}^{-1}$) than other 17 debris-free glaciers ($0.49 \pm 0.15 \text{ m w.e. a}^{-1}$) (Figure 7). An average terminus elevation for all debris-covered glaciers was 2466 m a.s.l., while it was 2661 m a.s.l. for debris-free glaciers.

## 4.2. *Comparison with Previous Observations*

Comparison with glaciological measurements from Djankuat Glacier [20] shows that values by Hugonnet et al. (2021) [3] from the period of 2007–2015 match within the uncertainty, while the rates from the other periods (2001, 2003, 2006, and 2017) show different changes (Figure 8a). In a total of twenty years, the mean (relative) difference between both methods is $55 \pm 35\%$ with a *p*-value of 0.172 and a standard deviation of 0.14 (Figure 8b).

The geodetic mass balance data from selected glaciers (Kyukyurtyu, Ulluchiran, Karachaul, and Bolshoy Azau) from the Elbrus Massif in 2000–2019 by Hugonnet et al. (2021) [3] agree with those from geodetic measurements in 1997–2017 by Kutuzov et al. (2019) [27] (Figure 8c). In most cases, values by Hugonnet et al. (2021) [3] are slightly lower (except the Ulluchiran Glacier) than those from Kutuzov et al. (2019) [27]; however, values from selected glaciers match within the uncertainty. Overall, the annual geodetic mass balance between 2000 and 2019 was $-0.35 \pm 0.20 \text{ m w.e. a}^{-1}$ for the entire Elbrus Massif according to Hugonnet et al.
(2021) [3], while it was −0.55 ± 0.04 m w.e. a⁻¹ between 1997 and 2017 according to Kutuzov et al. (2019) [27]. Over the last two decades, the mean (relative) difference between both methods is 57 ± 56% with a standard deviation of 0.10.

4.3. *Temperature and Precipitation Change*

Mean annual May-September air temperatures at the Terskol weather station show a significant positive trend in the years 2000–2018, rising by 0.7 °C from 2000–2009 (10.1 °C) to 2010–2018 (10.8 °C) (Figure 9a). The warming trend is also positive for the Mestia weather station, rising by 0.9 °C from 2000–2009 (14.0 °C) to 2010–2015 (14.9 °C). A significant negative trend is seen in the October-April monthly precipitation data for the Terskol station in the years 2000–2017, decreasing by 132 mm from 2000–2009 (592 mm) to 2010–2017 (460 mm). This was also negative for the Mestia station, decreasing by 58 mm from 2000–2009 (622 mm) to 2010–2014 (564 mm) (Figure 9b).

<figure><img src="image_11.png" /></figure>

**Figure 7.** Changes in mean annual mass balance (m w.e. a⁻¹) for debris-free (a) and debris-covered (b) glaciers vs. terminus elevation in 2000–2019. Red dots indicate the terminus elevation of each glacier. Vertical stems show an uncertainty estimation of the mean values.
<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 8.** Comparison of mean annual mass balance (m w.e. $a^{-1}$) (a) and cumulative mass balance (m w.e.) (b) for Djankuat Glacier using glaciological (blue) [20] and geodetic (brown) [3] methods. The light gray color corresponds to the uncertainty of the geodetic method. Comparison of mean annual mass changes (m w.e. $a^{-1}$) for Elbrus glaciers using geodetic methods (c) by Kutuzov et al. (2019) [27] (blue) and Hugonnet et al. (2021) (brown) [3]. Vertical stems on panel c show an uncertainty estimation of the mean values.</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>**Figure 9.** Mean annual May-September air temperatures (a) and sum of the October-April precipitation (b) at the Terskol and Mestia meteorological stations over the last twenty years. The dotted lines are the trend showing a temperature/precipitation course with time.</figcaption></figure>

**5. Discussion**

Results obtained from Hugonnet et al. (2021) [3] show that negative values of geodetic mass balance dominated between 2000 and 2019, but with positive values in some river basins and individual glaciers in 2000–2010 (Table A1). There was a negative and accelerated trend in the second investigated period (2010–2019) for almost all individual glaciers, with very few exceptions (e.g., Adishi, Khalde, and Kolka glaciers) (Table A3).

Overall, observed mass loss and thickness lowering of all Caucasus glaciers with increased values since 2010 are in good agreement with recent findings of the latest glacier inventory from this region, resulting in accelerated glacier area loss over the past two decades (1.16% yr⁻¹), especially since 2014 (1.85% yr⁻¹) [37].

**5.1. Causes of Regional Change**

The decrease of surface mass balance of glaciers on different sections of the Greater Caucasus was not homogeneous. The highest mean annual mass loss between 2000 and 2019 was observed in the western Greater Caucasus, while this was relatively low in the central and eastern sections. The highest mass loss in the western section can be explained by low hypsometry corresponding the lowermost mean elevation of the glaciers. A comparison of decrease rates relative to the glacier mean elevation clearly shows this linkage (Figure 10). Furthermore, western Greater Caucasus experienced the highest increase of supra-glacial debris-covered area over the past several decades [38]. This part of the region is also
characterized by the highest glacier retreat after the eastern Greater Caucasus [37]. It is,
therefore, possible that glaciers are also rapidly thinning, favouring debris-covered area
over the ablation area [40,41].

<figure><img src="image_16.png" /><figcaption>**Figure 10.** Change in mean annual geodetic mass balance (m w.e. $a^{-1}$) as a function of mean elevation (m a.s.l.) in western (a), central (b), and eastern (c) Greater Caucasus between 2000 and 2019. Elevation bands in 200 m. Vertical stems show an uncertainty estimation of the mean values.</figcaption></figure>

Similar decrease rates between the northern and southern slopes over the study period ($0.53 \pm 0.38$ m w.e. $a^{-1}$ on both slopes) may be due to the influence of the eastern Greater Caucasus. Glaciers on the eastern Greater Caucasus only exist on the northern slopes (i.e., glaciers do not exist on the southern slopes), where we observed the lowest mass loss. This may offset (reduce) the overall regional mass balance decrease rates of the glaciers from the northern Greater Caucasus relative to their southern counterparts.

The fact that mean annual mass balance on the glaciers of the Russian Caucasus was slightly higher ($-0.55 \pm 0.39$ m w.e. $a^{-1}$) than those on the Georgian counterparts ($-0.51 \pm 0.36$ m w.e. $a^{-1}$) can be explained by more favorable conditions of Georgian glaciers due to the prevailing southwestern transport of precipitation [32]. Moreover, glaciers situated in the western-central section of the Georgian Caucasus exist in high cirques with a much steeper surface. Glaciers in the western-central section of the Russian side are, on average, less steep than on the Georgian side, mainly because most valley glacier tongues in Russia are longer and reach lower altitudes than the Georgian glaciers [38], causing more ice melt on large glacier tongues on the Russian side.

## 5.2. *Glacier Mass Change Relative to Aspect and Size Class*

Glacier mass loss was in a linear relationship between aspect and mean elevation (Figure 6a), but in a non-linear relationship between size class and mean elevation. For example, glaciers with a size class of 0.5–1.0 and 1.0–5.0 km² experienced the highest mean annual mass loss; however they did not experience the lowest mean elevation (Figure 6b). A more linear relationship was found between size class and terminus elevation. For example, glaciers with the same size class (0.5–1.0 and 1.0–5.0 km²) had lower terminus elevation than the rest of the glaciers with a size class of <0.5 km². The lowest terminus elevation, but not the highest mass loss, with all glaciers > 5 km² can be explained due to the large accumulation areas situated in high altitudes of large glaciers (e.g., Elbrus, Bezingi, Dykhsu, Karaugom, etc.). Hence, the mass loss for large glaciers occurred at the expense of glacier tongues; meanwhile, the majority of regional mass loss occurred at the expense of relatively small (<5.0 km²) glaciers.

## 5.3. *Debris-Covered and Debris-Free Glaciers*

A comparison of geodetic mass balance change for debris-covered and debris-free glaciers during the investigated period showed that mean annual mass balance decrease rates are higher for debris-covered glaciers (Figure 7, Table A3). Although the debris-
covered glaciers, in general, present negative rates, the supraglacial debris clearly reduces the thinning compared with debris-free surfaces at the same elevation [42]. However, it is likely that altitudinal difference between those two types of glaciers can change this assumption and play a significant role in glacier thinning. For example, an average terminus elevation for the investigated debris-covered glaciers was ~200 m lower (2466 m a.s.l.) than those for the debris-free glaciers (2661 m a.s.l.), suggesting relatively warm temperatures for debris-covered glaciers based on the average regional lapse rate in summer (−0.52 °C per 100 m) [33]. Supra-glacial debris thickness can also play a key role in changes of mass balance for individual glaciers. This is broadly confirmed that supra-glacial debris cover affects surface melt with increasing ablation in cases of thin debris cover (less than a few cm) or decreasing ablation under continuous thick debris cover [43–45]. It is also possible that ice cliffs and supra-glacial ponds play a key role by enhancing the total ablation of debris-covered glaciers in the Greater Caucasus. Ice cliffs have been identified as an important contribution to the ablation of debris-covered glaciers, even when their spatial extent is very small [46]. This phenomenon was studied by terrestrial photogrammetry and unmanned air vehicle (UAV) surveys confirming that ice cliffs have a net ablation rate about 3 times higher than the average glacier tongue surface [42]. We note that further studies with high-resolution imagery (e.g., UAV, LiDAR) are required to understand the influence of ice cliffs in the mass balance for individual debris-covered glaciers, or at the basin and region scale in the Greater Caucasus.

### 5.4. Difference between the Global and Local Studies

The difference between Hugonnet et al. (2021) [3] and Kutuzov et al. (2019) [27] from Elbrus glaciers may be due to the high-resolution DEMs used by Kutuzov (2019) [27] and the internal boundaries of the selected glaciers, which were very similar, but not the same. In addition, the study by Hugonnet et al. (2021) [3] covers 2000–2019, which is slightly different (1997–2017) than in Kutuzov et al. (2019) [27]. Overall, the difference between those measurements match within the uncertainties, while the difference in geodetic mass balance from Hugonnet et al. (2021) [3] and glaciological mass balance from ground-based measurement [20] of the Djankuat Glacier was relatively high. This can be due to the much higher accuracy of ground-based measurement [47] or lower accuracy of geodetic measurement [3]. However, the results from those two measurements match within the uncertainties and tightly match the first time-period (2000–2010), suggesting that they are comparable in terms of long-time measurement.

### 5.5. *Glacier Mass Balance Relative to Local Climate Change*

Glacier mass balance in the Greater Caucasus is very sensitive to summer temperatures and winter precipitation [47]. Accelerated mass balance decrease from the second investigated periods (2010–2019) may be explained by a sustained increase in summer temperatures and a decrease in winter precipitation associated with regional global warming. The average summer temperature in the mountainous region of the Greater Caucasus has increased by 0.5–0.7 °C over the past 30 years [48], which was also confirmed from our observation from Terskol and Mestia meteorological stations (Tables 1 and 2; Figure 11). In parallel of the increased summer temperatures, extension of melting period was observed in the northern Greater Caucasus. For example, melting on the Garabashi Glacier at the snowline altitude (3800–4000 m a.s.l.) ended up 2–3 weeks later (in early October) over the recent years [17]. Furthermore, the extension of ablation season over the last two decades was also confirmed by instrumental measurement from the southern Greater Caucasus (Zopkhito and Chalaati glaciers) (Figure 7 in [19]).
<table><caption>**Table 1.** Mean monthly air temperatures (°C) for the Terskol and Mestia meteorological stations.</caption><thead><tr><th></th><th>Jan</th><th>Feb</th><th>Mar</th><th>Apr</th><th>May</th><th>Jun</th><th>Jul</th><th>Aug</th><th>Sep</th><th>Oct</th><th>Nov</th><th>Dec</th></tr></thead><tbody><tr><td colspan="13" style="text-align: center; font-weight: bold;">Terskol</td></tr><tr><td>2000–2009</td><td>-8.6</td><td>-7.3</td><td>-3.3</td><td>1.2</td><td>6.4</td><td>10.2</td><td>12.9</td><td>12.4</td><td>8.5</td><td>4.0</td><td>-2.1</td><td>-6.4</td></tr><tr><td>2010–2018</td><td>-6.7</td><td>-5.4</td><td>-2.1</td><td>2.0</td><td>7.5</td><td>11.1</td><td>13.4</td><td>12.9</td><td>9.2</td><td>3.8</td><td>-1.8</td><td>-5.1</td></tr><tr><td>Change</td><td>2.0</td><td>1.9</td><td>1.1</td><td>0.8</td><td>1.1</td><td>0.9</td><td>0.6</td><td>0.5</td><td>0.6</td><td>-0.2</td><td>0.3</td><td>1.3</td></tr><tr><td colspan="13" style="text-align: center; font-weight: bold;">Mestia</td></tr><tr><td>2000–2009</td><td>-5.9</td><td>-3.3</td><td>0.1</td><td>5.4</td><td>10.4</td><td>13.7</td><td>16.7</td><td>16.7</td><td>12.6</td><td>8.1</td><td>1.3</td><td>-4.0</td></tr><tr><td>2010–2015</td><td>-3.7</td><td>-2.6</td><td>0.7</td><td>5.6</td><td>11.5</td><td>14.9</td><td>17.4</td><td>17.4</td><td>13.2</td><td>7.6</td><td>2.7</td><td>-2.9</td></tr><tr><td>Change</td><td>2.2</td><td>0.6</td><td>0.6</td><td>0.2</td><td>1.1</td><td>1.2</td><td>0.7</td><td>0.7</td><td>0.6</td><td>-0.5</td><td>1.3</td><td>1.2</td></tr></tbody></table>

<table><caption>**Table 2.** Sum of the monthly precipitation (mm) for the Terskol meteorological station in 2000–2009 and 2010–2017, and for the Mestia meteorological station in 2000–2009 and 2010–2014.</caption><thead><tr><th></th><th>Jan</th><th>Feb</th><th>Mar</th><th>Apr</th><th>May</th><th>Jun</th><th>Jul</th><th>Aug</th><th>Sep</th><th>Oct</th><th>Nov</th><th>Dec</th></tr></thead><tbody><tr><td colspan="13" style="text-align: center; font-weight: bold;">Terskol</td></tr><tr><td>2000–2009</td><td>530.1</td><td>692.7</td><td>1144.9</td><td>1005.2</td><td>1019.2</td><td>994.9</td><td>1089.0</td><td>929.7</td><td>1063.2</td><td>1167.4</td><td>807.0</td><td>568.0</td></tr><tr><td>2010–2017</td><td>584.4</td><td>211.1</td><td>531.2</td><td>815.6</td><td>834.1</td><td>735.4</td><td>738.8</td><td>810.9</td><td>798.2</td><td>794.4</td><td>331.1</td><td>413.0</td></tr><tr><td>Change</td><td>54.3</td><td>-481.6</td><td>-613.7</td><td>-189.6</td><td>-185.1</td><td>-259.5</td><td>-350.2</td><td>-118.8</td><td>-265.0</td><td>-373.0</td><td>-475.9</td><td>-155.0</td></tr><tr><td colspan="13" style="text-align: center; font-weight: bold;">Mestia</td></tr><tr><td>2000–2009</td><td>623</td><td>833.5</td><td>1150.8</td><td>1030.9</td><td>1058.2</td><td>974.7</td><td>927.8</td><td>945</td><td>930.2</td><td>1145.9</td><td>714.6</td><td>725.7</td></tr><tr><td>2010–2014</td><td>279.5</td><td>386.4</td><td>462.5</td><td>441.8</td><td>548.9</td><td>634.4</td><td>518.4</td><td>478.3</td><td>608.1</td><td>800.8</td><td>276.4</td><td>172.4</td></tr><tr><td>Change</td><td>-343.5</td><td>-447.1</td><td>-688.3</td><td>-589.1</td><td>-509.3</td><td>-340.3</td><td>-409.4</td><td>-466.7</td><td>-322.1</td><td>-345.1</td><td>-438.2</td><td>-553.3</td></tr></tbody></table>

<figure><img src="image_17.png" /><figcaption>**Figure 11.** Average May–September monthly air temperatures for the Terskol meteorological station in 2000–2009 and 2010–2018 (a) and for the Mestia meteorological station in 2000–2009 and 2010–2015 (b); Sum of the October–April precipitation for the Terskol meteorological station in 2000–2009 and 2010–2017 (c) and for the Mestia meteorological station in 2000–2009 and 2010–2014 (d).</figcaption></figure>
The increased incoming shortwave radiation, observed since the 1980s (10 W/m² over 10 years), may also play a triggering role in the accelerated mass loss of the Caucasus glaciers over the past several decades. It has been proposed that this trend is associated with a weakening of the processes of formation of high and low clouds, which is due to an increase in the frequency of anticyclones in the warm season [49].

Snow-ice-covered areas of the Greater Caucasus are subject to dust outbreak events originating from desert regions able to significantly decrease snow albedo. Several dust depositions events that originated from the Sahara between the 2009–2012 and 2018 period were observed in this region [50,51]. Dust events systematically shortened the snow cover duration, especially in the western Greater Caucasus. The shortening was higher for locations with higher accumulation and higher elevation [51] considering one of the significant indicators for glacier thinning in this region.

The decrease in albedo, along with increased temperatures over the ablation period causes up-glacier migration of the snowline. This has a strong impact on small glaciers that have practically lost their firn reserves over the past two decades. Consequently, the firn can no longer protect the ice, and the process of ice melting becomes even more intense due to the lower albedo (Figure 12).

<figure><img src="image_18.png" /><figcaption>**Figure 12.** Khasankoi glaciers (north-western Greater Caucasus) during the late summer without snowlines. Sentinel 2B image (04/09/2020) is used as a background.</figcaption></figure>

**6. Conclusions**

In this study, we used glacier mass balance and thickness data obtained from geodetic method by Hugonnet et al. (2021) [3] to clarify the response of region-wide glacier mass balance to climate change in the Greater Caucasus over the last two decades (2000–2019). We also analyzed the long-term spatio-temporal glacier surface elevation and mass change for the regional, river basin, and individual glaciers. Besides, we provided a key dataset (see Supplement) that can be useful for further studies, including, for example, hydrological modeling of runoff changes, as well as for the implementation of water management plans in this region. The main limitation of this study is related to the resolution of the data. This generally affects small glaciers more than larger complexes. Applying relatively high uncertainties was a conservative solution to this issue. This can be improved by using the high-resolution stereo imagery and continuous ground-based investigation in the future. However, the estimates of Djankuat and Elbrus glaciers from this dataset agree with glaciological and geodetic mass balance measurements from previous studies.

The main results indicated a mean regional glacier mass balance change of -0.53 ± 0.38 m w.e. a⁻¹ in 2000–2019 with an equivalent cumulative mass change of -10.07 ± 7.22 m w.e. Calculated loss in the mass of all glaciers in the Greater Caucasus over
the past 20 years amounted to the ~12.99 km³ or 11.04 ± 0.78 Gt water. Strong acceleration in mass loss and the thinning rates were observed since 2010, which likely occurred due to an increase in summer temperatures, accompanied by a reduction in winter precipitation.

The highest mass loss in the western Greater Caucasus was due to the lowest mean elevation of glaciers. We also hypothesized that the deposition events that originated from the Sahara between the 2009–2012 and 2018 periods may have increased glacier melting in the western section.

Lower terminus elevation in debris-covered glaciers indicated a strong decrease of surface mass balance relative to the high-terminated debris-free glaciers. Debris thickness, ice cliffs, and supraglacial ponds may also play a key role in enhancing the total ablation of debris-covered glaciers in the Greater Caucasus.

**Author Contributions:** Conceptualization, L.G.T. and V.J.; methodology, L.G.T. and G.A.N.; software, L.G.T.; validation, L.G.T., V.J. and G.A.N.; formal analysis, L.G.T., V.J. and G.A.N.; investigation, L.G.T., V.J. and G.A.N.; writing—original draft preparation, L.G.T.; review and editing, L.G.T., V.J. and G.A.N.; visualization, L.G.T., V.J. and G.A.N.; funding acquisition, V.J. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the French-Russian-Georgian collaborative project "DEGLaciation dans le grAnd Caucase—DEGLAC" (IRP 00008 Deglac).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are openly available in Zenodo https://doi.org/10.5281/zenodo.5816997 (Last access 2 February 2022).

**Acknowledgments:** We gratefully acknowledge the support of the two anonymous reviewers for useful suggestions and detailed comments which clearly enhanced the quality of the paper.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.
Appendix A

<table><caption>**Table A1.** Change in glacier thickness (m a⁻¹) and surface mass balance (m w.e. a⁻¹) for individual river basins of the Greater Caucasus in 2000–2010, 2010–2019, and 2000–2019. See also Figure 1 to check the location of all river basins relative to the Greater Caucasus.</caption><thead><tr><th rowspan="2">Main River Basin</th><th rowspan="2">Tributary River Basin</th><th colspan="2">Based on Landsat 2014<br>(Tielidze and Wheate, 2018) [21]</th><th colspan="3">Change of Glacier Thickness (m a<sup>-1</sup>)<br>Hugonnet et al. (2021) [3]</th><th colspan="3">Change in Glacier Surface Mass Balance (m w.e. a<sup>-1</sup>)<br>Hugonnet et al. (2021) [3]</th></tr><tr><th>Number</th><th>Area km<sup>2</sup></th><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th></tr></thead><tbody><tr><td rowspan="2">Mzimta</td><td>Belaya</td><td>41</td><td>6.5 ± 0.4</td><td>-1.05 ± 1.04</td><td>-1.01 ± 1.06</td><td>-1.03 ± 0.57</td><td>-0.89 ± 1.05</td><td>-0.86 ± 1.04</td><td>-0.88 ± 0.67</td></tr><tr><td>Malaya Laba</td><td>34</td><td>7.4 ± 0.5</td><td>-0.73 ± 0.88</td><td>-1.12 ± 1.02</td><td>-0.93 ± 0.52</td><td>-0.62 ± 0.83</td><td>-0.96 ± 0.97</td><td>-0.79 ± 0.57</td></tr><tr><td rowspan="7">Kuban</td><td>Bolshaya Laba</td><td>27</td><td>4.4 ± 0.2</td><td>-1.57 ± 1.10</td><td>-1.04 ± 1.34</td><td>-1.30 ± 0.67</td><td>-1.34 ± 1.20</td><td>-0.88 ± 1.31</td><td>-1.11 ± 0.84</td></tr><tr><td>Bolshoy Zelenchuk</td><td>59</td><td>24.6 ± 1.3</td><td>-1.12 ± 0.94</td><td>-0.94 ± 1.11</td><td>-1.03 ± 0.56</td><td>-0.95 ± 0.89</td><td>-0.80 ± 1.02</td><td>-0.88 ± 0.59</td></tr><tr><td>Maliy Zelenchuk</td><td>42</td><td>24.5 ± 1.1</td><td>-0.70 ± 0.67</td><td>-0.96 ± 0.73</td><td>-0.83 ± 0.38</td><td>-0.59 ± 0.62</td><td>-0.81 ± 0.69</td><td>-0.70 ± 0.40</td></tr><tr><td>Teberda</td><td>89</td><td>52.9 ± 2.5</td><td>-0.47 ± 0.59</td><td>-0.90 ± 0.64</td><td>-0.68 ± 0.34</td><td>-0.40 ± 0.53</td><td>-0.77 ± 0.60</td><td>-0.58 ± 0.34</td></tr><tr><td>Daut</td><td>23</td><td>4.8 ± 0.4</td><td>-0.49 ± 0.59</td><td>-0.77 ± 0.65</td><td>-0.63 ± 0.34</td><td>-0.42 ± 0.54</td><td>-0.65 ± 0.61</td><td>-0.54 ± 0.35</td></tr><tr><td>Uchkulan</td><td>68</td><td>19.1 ± 1.4</td><td>-0.38 ± 0.58</td><td>-0.99 ± 0.63</td><td>-0.68 ± 0.33</td><td>-0.32 ± 0.51</td><td>-0.84 ± 0.61</td><td>-0.58 ± 0.35</td></tr><tr><td>Ullukam</td><td>93</td><td>42.8 ± 2.2</td><td>-0.38 ± 0.57</td><td>-0.69 ± 0.59</td><td>-0.53 ± 0.32</td><td>-0.32 ± 0.50</td><td>-0.59 ± 0.55</td><td>-0.45 ± 0.31</td></tr><tr><td>Malka</td><td></td><td>12</td><td>52.2 ± 0.8</td><td>-0.38 ± 0.50</td><td>-0.45 ± 0.50</td><td>-0.42 ± 0.27</td><td>-0.33 ± 0.44</td><td>-0.39 ± 0.45</td><td>-0.36 ± 0.26</td></tr><tr><td>Baksan</td><td></td><td>147</td><td>146.2 ± 6.0</td><td>-0.51 ± 0.57</td><td>-0.74 ± 0.56</td><td>-0.62 ± 0.30</td><td>-0.44 ± 0.52</td><td>-0.63 ± 0.52</td><td>-0.53 ± 0.31</td></tr><tr><td rowspan="3">Chegem</td><td>Bashil-Auzusu</td><td>30</td><td>25.3 ± 1.0</td><td>-0.39 ± 0.64</td><td>-0.90 ± 0.61</td><td>-0.65 ± 0.33</td><td>-0.34 ± 0.56</td><td>-0.77 ± 0.57</td><td>-0.55 ± 0.33</td></tr><tr><td>Gara-Auzusu</td><td>15</td><td>24.6 ± 0.9</td><td>-0.30 ± 0.55</td><td>-0.65 ± 0.55</td><td>-0.47 ± 0.30</td><td>-0.26 ± 0.48</td><td>-0.55 ± 0.49</td><td>-0.40 ± 0.28</td></tr><tr><td>Bulungu</td><td>8</td><td>2.9 ± 0.2</td><td>0.04 ± 0.82</td><td>-0.90 ± 0.94</td><td>-0.43 ± 0.49</td><td>0.04 ± 0.74</td><td>-0.77 ± 0.90</td><td>-0.37 ± 0.48</td></tr><tr><td rowspan="5">Cherek</td><td>Cherek-Bezingskiy</td><td>38</td><td>72.9 ± 2.2</td><td>-0.24 ± 0.70</td><td>-0.51 ± 0.77</td><td>-0.37 ± 0.41</td><td>-0.21 ± 0.61</td><td>-0.43 ± 0.71</td><td>-0.32 ± 0.39</td></tr><tr><td>Cherek-Balkarskiy</td><td>80</td><td>96.5 ± 3.5</td><td>-0.51 ± 0.65</td><td>-0.47 ± 0.68</td><td>-0.49 ± 0.37</td><td>-0.43 ± 0.59</td><td>-0.40 ± 0.60</td><td>-0.42 ± 0.36</td></tr><tr><td>Psigansy</td><td>11</td><td>12.7 ± 0.6</td><td>-0.56 ± 0.73</td><td>-0.61 ± 0.75</td><td>-0.59 ± 0.35</td><td>-0.48 ± 0.66</td><td>-0.52 ± 0.67</td><td>-0.50 ± 0.36</td></tr><tr><td>Khanzidon</td><td>10</td><td>6.2 ± 0.4</td><td>-0.33 ± 0.70</td><td>-0.42 ± 0.74</td><td>-0.38 ± 0.38</td><td>-0.28 ± 0.66</td><td>-0.36 ± 0.73</td><td>-0.32 ± 0.40</td></tr><tr><td>Biliagikom</td><td>5</td><td>1.0 ± 0.1</td><td>-0.88 ± 0.74</td><td>-1.11 ± 0.76</td><td>-1.00 ± 0.41</td><td>-0.75 ± 0.67</td><td>-0.94 ± 0.70</td><td>-0.85 ± 0.42</td></tr><tr><td rowspan="4">Urukh</td><td>Urukh headwaters</td><td>31</td><td>25.0 ± 1.2</td><td>-0.54 ± 0.65</td><td>-0.59 ± 0.64</td><td>-0.57 ± 0.34</td><td>-0.46 ± 0.58</td><td>-0.50 ± 0.57</td><td>-0.48 ± 0.32</td></tr><tr><td>Karaugom</td><td>29</td><td>41.7 ± 1.4</td><td>-0.76 ± 0.72</td><td>-0.45 ± 0.71</td><td>-0.61 ± 0.39</td><td>-0.65 ± 0.67</td><td>-0.38 ± 0.62</td><td>-0.51 ± 0.38</td></tr><tr><td>Aigamuga</td><td>20</td><td>14.7 ± 0.8</td><td>-0.49 ± 0.79</td><td>-0.41 ± 0.79</td><td>-0.45 ± 0.42</td><td>-0.42 ± 0.71</td><td>-0.35 ± 0.70</td><td>-0.38 ± 0.41</td></tr><tr><td>Tseyadon</td><td>20</td><td>16.8 ± 0.6</td><td>-0.49 ± 0.92</td><td>-0.43 ± 0.88</td><td>-0.46 ± 0.48</td><td>-0.42 ± 0.81</td><td>-0.36 ± 0.77</td><td>-0.39 ± 0.44</td></tr><tr><td rowspan="3">Ardon</td><td>Sidan</td><td>1</td><td>0.03 ± 0.004</td><td>-0.07 ± 1.50</td><td>-0.74 ± 1.47</td><td>-0.41 ± 0.80</td><td>-0.06 ± 1.28</td><td>-0.63 ± 1.31</td><td>-0.35 ± 0.71</td></tr><tr><td>Vilsa</td><td>4</td><td>0.9 ± 0.08</td><td>-0.79 ± 0.76</td><td>-0.71 ± 0.72</td><td>-0.75 ± 0.40</td><td>-0.67 ± 0.71</td><td>-0.61 ± 0.68</td><td>-0.64 ± 0.44</td></tr><tr><td>Adaikom</td><td>7</td><td>4.6 ± 0.2</td><td>-0.33 ± 0.64</td><td>-0.79 ± 0.60</td><td>-0.56 ± 0.33</td><td>-0.28 ± 0.56</td><td>-0.67 ± 0.56</td><td>-0.47 ± 0.32</td></tr></tbody></table>
<table><caption>Table A1. *Cont.*</caption><thead><tr><th rowspan="2">Main River Basin</th><th rowspan="2">Tributary River Basin</th><th colspan="2">Based on Landsat 2014<br>(Tielidze and Wheate, 2018) [21]</th><th colspan="3">Change of Glacier Thickness (m a<sup>-1</sup>)<br>Hugonnet et al. (2021) [3]</th><th colspan="3">Change in Glacier Surface Mass Balance (m w.e. a<sup>-1</sup>)<br>Hugonnet et al. (2021) [3]</th></tr><tr><th>Number</th><th>Area km<sup>2</sup></th><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th></tr></thead><tbody><tr><td></td><td>Mamikhdon</td><td>15</td><td>3.7 ± 0.2</td><td>-0.82 ± 0.68</td><td>-0.79 ± 0.66</td><td>-0.81 ± 0.36</td><td>-0.70 ± 0.67</td><td>-0.67 ± 0.60</td><td>-0.69 ± 0.38</td></tr><tr><td></td><td>Nar</td><td>10</td><td>1.5 ± 0.07</td><td>-0.24 ± 0.66</td><td>-0.67 ± 0.62</td><td>-0.45 ± 0.34</td><td>-0.20 ± 0.61</td><td>-0.57 ± 0.59</td><td>-0.38 ± 0.36</td></tr><tr><td></td><td>Kasaidon</td><td>2</td><td>0.2 ± 0.02</td><td>-0.71 ± 0.91</td><td>-0.39 ± 0.87</td><td>-0.55 ± 0.47</td><td>-0.60 ± 0.81</td><td>-0.33 ± 0.76</td><td>-0.47 ± 0.44</td></tr><tr><td></td><td>Labagkomdon</td><td>1</td><td>0.06 ± 0.008</td><td>0.51 ± 1.01</td><td>-0.48 ± 0.93</td><td>0.02 ± 0.52</td><td>0.43 ± 0.87</td><td>-0.41 ± 0.81</td><td>0.01 ± 0.44</td></tr><tr><td></td><td>Baddon</td><td>8</td><td>2.6 ± 0.2</td><td>-0.66 ± 0.71</td><td>-0.45 ± 0.68</td><td>-0.55 ± 0.37</td><td>-0.56 ± 0.63</td><td>-0.38 ± 0.61</td><td>-0.47 ± 0.36</td></tr><tr><td></td><td>Arkhondon</td><td>4</td><td>3.5 ± 0.2</td><td>-0.95 ± 0.54</td><td>-0.63 ± 0.52</td><td>-0.79 ± 0.29</td><td>-0.81 ± 0.49</td><td>-0.54 ± 0.46</td><td>-0.67 ± 0.28</td></tr><tr><td></td><td>Fiagdon</td><td>25</td><td>6.2 ± 0.6</td><td>-0.23 ± 0.66</td><td>-0.50 ± 0.65</td><td>-0.36 ± 0.35</td><td>-0.20 ± 0.59</td><td>-0.42 ± 0.59</td><td>-0.31 ± 0.34</td></tr><tr><td></td><td>Gizeldon</td><td>24</td><td>25.2 ± 1.0</td><td>-0.46 ± 0.56</td><td>-0.43 ± 0.55</td><td>-0.44 ± 0.30</td><td>-0.39 ± 0.50</td><td>-0.36 ± 0.50</td><td>-0.38 ± 0.28</td></tr><tr><td></td><td>Tergi (Terek) headwaters</td><td>94</td><td>42.7 ± 3.0</td><td>-0.04 ± 0.68</td><td>-0.64 ± 0.70</td><td>-0.34 ± 0.37</td><td>-0.03 ± 0.60</td><td>-0.54 ± 0.65</td><td>-0.29 ± 0.36</td></tr><tr><td rowspan="2">Sunja right tributaries</td><td>Assa</td><td>10</td><td>2.1 ± 0.2</td><td>-0.13 ± 0.64</td><td>-0.35 ± 0.68</td><td>-0.24 ± 0.36</td><td>-0.11 ± 0.56</td><td>-0.30 ± 0.59</td><td>-0.20 ± 0.31</td></tr><tr><td>Arghuni</td><td>18</td><td>5.2 ± 0.4</td><td>-0.10 ± 0.66</td><td>-0.55 ± 0.67</td><td>-0.33 ± 0.36</td><td>-0.08 ± 0.58</td><td>-0.47 ± 0.60</td><td>-0.28 ± 0.33</td></tr><tr><td rowspan="2">Sulak</td><td>Sharo Argun</td><td>31</td><td>15.7 ± 0.8</td><td>-0.57 ± 0.83</td><td>-0.64 ± 0.82</td><td>-0.61 ± 0.43</td><td>-0.49 ± 0.75</td><td>-0.54 ± 0.73</td><td>-0.52 ± 0.41</td></tr><tr><td>Andiyskoye Koysu</td><td>35</td><td>11.8 ± 0.6</td><td>-0.70 ± 0.68</td><td>-0.68 ± 0.66</td><td>-0.69 ± 0.35</td><td>-0.60 ± 0.70</td><td>-0.58 ± 0.60</td><td>-0.59 ± 0.39</td></tr><tr><td></td><td>Avarskoye Koysu</td><td>34</td><td>5.9 ± 0.4</td><td>-0.53 ± 0.76</td><td>-0.54 ± 0.79</td><td>-0.53 ± 0.42</td><td>-0.45 ± 0.69</td><td>-0.46 ± 0.72</td><td>-0.45 ± 0.41</td></tr><tr><td></td><td>Samur</td><td>13</td><td>2.4 ± 0.2</td><td>-0.19 ± 0.54</td><td>-0.71 ± 0.55</td><td>-0.45 ± 0.30</td><td>-0.16 ± 0.47</td><td>-0.60 ± 0.51</td><td>-0.38 ± 0.28</td></tr><tr><td></td><td>Kusarchai</td><td>14</td><td>0.8 ± 0.1</td><td>-0.21 ± 0.62</td><td>-0.34 ± 0.59</td><td>-0.28 ± 0.32</td><td>-0.18 ± 0.54</td><td>-0.29 ± 0.55</td><td>-0.24 ± 0.30</td></tr><tr><td></td><td>Bzipi</td><td>16</td><td>3.1 ± 0.3</td><td>-1.07 ± 0.73</td><td>-0.22 ± 0.86</td><td>-0.65 ± 0.44</td><td>-0.91 ± 0.73</td><td>-0.19 ± 0.81</td><td>-0.55 ± 0.45</td></tr><tr><td></td><td>Kelasuri</td><td>2</td><td>0.9 ± 0.06</td><td>-1.04 ± 0.69</td><td>-1.11 ± 0.90</td><td>-1.07 ± 0.45</td><td>-0.89 ± 0.62</td><td>-0.94 ± 0.85</td><td>-0.91 ± 0.47</td></tr><tr><td></td><td>Kodori</td><td>154</td><td>42.0 ± 2.6</td><td>-0.46 ± 0.60</td><td>-0.96 ± 0.68</td><td>-0.71 ± 0.36</td><td>-0.39 ± 0.54</td><td>-0.81 ± 0.67</td><td>-0.60 ± 0.37</td></tr><tr><td></td><td>Enguri</td><td>289</td><td>223.6 ± 8.8</td><td>-0.31 ± 0.60</td><td>-0.84 ± 0.60</td><td>-0.57 ± 0.32</td><td>-0.26 ± 0.53</td><td>-0.71 ± 0.58</td><td>-0.49 ± 0.33</td></tr><tr><td></td><td>Khobistkali</td><td>11</td><td>0.6 ± 0.08</td><td>-1.08 ± 0.97</td><td>-0.71 ± 0.97</td><td>-0.89 ± 0.52</td><td>-0.92 ± 0.97</td><td>-0.60 ± 0.89</td><td>-0.76 ± 0.57</td></tr><tr><td></td><td>Rioni</td><td>117</td><td>50.6 ± 2.3</td><td>-0.64 ± 0.71</td><td>-0.61 ± 0.73</td><td>-0.62 ± 0.39</td><td>-0.54 ± 0.67</td><td>-0.52 ± 0.66</td><td>-0.53 ± 0.39</td></tr><tr><td></td><td>Liakhvi</td><td>12</td><td>2.0 ± 0.2</td><td>-0.22 ± 0.72</td><td>-0.76 ± 0.69</td><td>-0.49 ± 0.38</td><td>-0.18 ± 0.62</td><td>-0.65 ± 0.65</td><td>-0.42 ± 0.37</td></tr><tr><td></td><td>Aragvi</td><td>1</td><td>0.3 ± 0.04</td><td>0.06 ± 0.47</td><td>-0.56 ± 0.49</td><td>-0.25 ± 0.26</td><td>0.05 ± 0.40</td><td>-0.47 ± 0.43</td><td>-0.21 ± 0.23</td></tr><tr><td colspan="2">Total, Greater Caucasus</td><td>1861</td><td>1186.1 ± 53.3</td><td>-0.50 ± 0.67</td><td>-0.75 ± 0.70</td><td>-0.62 ± 0.37</td><td>-0.42 ± 0.61</td><td>-0.64 ± 0.66</td><td>-0.53 ± 0.38</td></tr></tbody></table>
<table><caption>**Table A2.** Change in glacier thickness (m a⁻¹) and surface mass balance (m w.e. a⁻¹) for Caucasus countries (Russia, Georgia, Azerbaijan) in 2000–2010, 2010–2019, and 2000–2019.</caption><thead><tr><th rowspan="2">Countries</th><th colspan="2">Based on Landsat 2014 (Tielidze and Wheate, 2018) [21]</th><th colspan="3">Change of Glacier Thickness (m a<sup>-1</sup>)<br>Hugonnet et al. (2021) [3]</th><th colspan="3">Change in Glacier Surface Mass Balance (m w.e. a<sup>-1</sup>)<br>Hugonnet et al. (2021) [3]</th></tr><tr><th>Number</th><th>Area km<sup>2</sup></th><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th></tr></thead><tbody><tr><td>Russia</td><td>1161</td><td>817.8 ± 35.1</td><td>-0.55 ± 0.70</td><td>-0.74 ± 0.73</td><td>-0.65 ± 0.38</td><td>-0.47 ± 0.63</td><td>-0.63 ± 0.67</td><td>-0.55 ± 0.39</td></tr><tr><td>Georgia</td><td>683</td><td>367.1 ± 16.8</td><td>-0.41 ± 0.64</td><td>-0.77 ± 0.67</td><td>-0.60 ± 0.36</td><td>-0.35 ± 0.58</td><td>-0.66 ± 0.64</td><td>-0.51 ± 0.36</td></tr><tr><td>Azerbaijan</td><td>17</td><td>1.2 ± 0.2</td><td>-0.17 ± 0.61</td><td>-0.37 ± 0.58</td><td>-0.27 ± 0.31</td><td>-0.14 ± 0.53</td><td>-0.31 ± 0.53</td><td>-0.23 ± 0.29</td></tr></tbody></table>

<table><caption>**Table A3.** Change in glacier thickness (m a⁻¹) and surface mass balance (m w.e. a⁻¹) for individual glaciers > 2.0 km² of the Greater Caucasus in 2000–2010, 2010–2019, and 2000–2019.</caption><thead><tr><th rowspan="2">Glacier Name</th><th rowspan="2">GLIMS_id</th><th rowspan="2">Area (km²)</th><th colspan="3">Change of Glacier Thickness (m a<sup>-1</sup>)<br>Hugonnet et al. (2021) [3]</th><th colspan="3">Change in Glacier Surface Mass Balance (m w.e. a<sup>-1</sup>)<br>Hugonnet et al. (2021) [3]</th></tr><tr><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th></tr></thead><tbody><tr><td colspan="9" style="text-align: left; font-weight: bold;">Glaciers with >10 km² Area</td></tr><tr><td>Bezingi</td><td>G043100E43030N</td><td>37.49</td><td>-0.22 ± 0.26</td><td>-0.78 ± 0.24</td><td>-0.50 ± 0.14</td><td>-0.19 ± 0.22</td><td>-0.67 ± 0.21</td><td>-0.43 ± 0.12</td></tr><tr><td>Dykhsu</td><td>G043183E42984N</td><td>27.53</td><td>-0.48 ± 0.27</td><td>-0.76 ± 0.26</td><td>-0.62 ± 0.15</td><td>-0.41 ± 0.23</td><td>-0.65 ± 0.23</td><td>-0.53 ± 0.14</td></tr><tr><td>Karaugom</td><td>G043748E42782N</td><td>23.99</td><td>-0.28 ± 0.28</td><td>-0.69 ± 0.26</td><td>-0.49 ± 0.15</td><td>-0.24 ± 0.24</td><td>-0.59 ± 0.23</td><td>-0.41 ± 0.14</td></tr><tr><td>Lekhziri</td><td>G042766E43151N</td><td>23.77</td><td>-0.46 ± 0.26</td><td>-1.44 ± 0.22</td><td>-0.95 ± 0.13</td><td>-0.39 ± 0.23</td><td>-1.22 ± 0.23</td><td>-0.81 ± 0.14</td></tr><tr><td>Dzhikiugankez</td><td>G042506E43356N</td><td>21.44</td><td>-1.01 ± 0.23</td><td>-1.45 ± 0.22</td><td>-1.23 ± 0.12</td><td>-0.86 ± 0.21</td><td>-1.23 ± 0.21</td><td>-1.04 ± 0.13</td></tr><tr><td>Agashtan</td><td>G043291E42928N</td><td>18.94</td><td>-1.04 ± 0.28</td><td>-1.02 ± 0.31</td><td>-1.03 ± 0.18</td><td>-0.89 ± 0.25</td><td>-0.86 ± 0.28</td><td>-0.88 ± 0.17</td></tr><tr><td>Azay Bolshoy</td><td>G042429E43293N</td><td>17.72</td><td>-0.39 ± 0.24</td><td>-0.69 ± 0.23</td><td>-0.54 ± 0.13</td><td>-0.33 ± 0.21</td><td>-0.59 ± 0.21</td><td>-0.46 ± 0.11</td></tr><tr><td>Ulluchiran</td><td>G042432E43372N</td><td>12.79</td><td>-0.08 ± 0.25</td><td>-0.29 ± 0.23</td><td>-0.19 ± 0.13</td><td>-0.07 ± 0.21</td><td>-0.25 ± 0.20</td><td>-0.16 ± 0.11</td></tr><tr><td>Midjirgi</td><td>G043165E43066N</td><td>12.74</td><td>-0.05 ± 0.30</td><td>-0.35 ± 0.33</td><td>-0.20 ± 0.18</td><td>-0.04 ± 0.26</td><td>-0.30 ± 0.28</td><td>-0.17 ± 0.15</td></tr><tr><td>Tsaneri Southern</td><td>G042999E43050N</td><td>12.32</td><td>-0.51 ± 0.29</td><td>-0.67 ± 0.29</td><td>-0.59 ± 0.16</td><td>-0.43 ± 0.25</td><td>-0.57 ± 0.26</td><td>-0.50 ± 0.15</td></tr><tr><td>Tseya</td><td>G043844E42752N</td><td>11.87</td><td>-0.21 ± 0.31</td><td>-1.04 ± 0.28</td><td>-0.63 ± 0.17</td><td>-0.18 ± 0.26</td><td>-0.89 ± 0.26</td><td>-0.53 ± 0.15</td></tr><tr><td>Tsaneri Northern</td><td>G043005E43089N</td><td>11.28</td><td>-0.46 ± 0.29</td><td>-0.45 ± 0.26</td><td>-0.46 ± 0.14</td><td>-0.39 ± 0.25</td><td>-0.38 ± 0.22</td><td>-0.39 ± 0.13</td></tr><tr><td>Irik</td><td>G042501E43334N</td><td>10.97</td><td>-0.34 ± 0.26</td><td>-0.75 ± 0.25</td><td>-0.54 ± 0.13</td><td>-0.29 ± 0.22</td><td>-0.63 ± 0.22</td><td>-0.46 ± 0.12</td></tr><tr><td colspan="9" style="text-align: left; font-weight: bold;">Glaciers with 5–10 km² area</td></tr><tr><td>Shaurty</td><td>G043051E43105N</td><td>9.89</td><td>-0.37 ± 0.29</td><td>-0.80 ± 0.26</td><td>-0.58 ± 0.15</td><td>-0.31 ± 0.25</td><td>-0.68 ± 0.24</td><td>-0.49 ± 0.14</td></tr><tr><td>Adishi</td><td>G043016E43017N</td><td>9.59</td><td>-0.07 ± 0.34</td><td>0.32 ± 0.33</td><td>0.12 ± 0.18</td><td>-0.06 ± 0.29</td><td>0.27 ± 0.28</td><td>0.10 ± 0.15</td></tr><tr><td>Kvitlodi</td><td>G042950E43105N</td><td>9.58</td><td>-0.10 ± 0.29</td><td>-0.88 ± 0.30</td><td>-0.49 ± 0.17</td><td>-0.08 ± 0.25</td><td>-0.75 ± 0.27</td><td>-0.42 ± 0.15</td></tr><tr><td>Terskol</td><td>G042488E43308N</td><td>9.48</td><td>-0.12 ± 0.26</td><td>-0.21 ± 0.25</td><td>-0.17 ± 0.14</td><td>-0.11 ± 0.23</td><td>-0.17 ± 0.21</td><td>-0.14 ± 0.12</td></tr></tbody></table>
<table><caption>Table A3. *Cont.*</caption><thead><tr><th rowspan="2">Glacier Name</th><th rowspan="2">GLIMS_id</th><th rowspan="2">Area (km²)</th><th colspan="3">Change of Glacier Thickness (m a⁻¹)<br>Hugonnet et al. (2021) [3]</th><th colspan="3">Change in Glacier Surface Mass Balance (m w.e. a⁻¹)<br>Hugonnet et al. (2021) [3]</th></tr><tr><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th></tr></thead><tbody><tr><td>Azay Maliy</td><td>G042449E43302N</td><td>9.45</td><td>0.00 ± 0.27</td><td>-0.22 ± 0.26</td><td>-0.11 ± 0.14</td><td>0.00 ± 0.23</td><td>-0.19 ± 0.22</td><td>-0.09 ± 0.12</td></tr><tr><td>Chalaati</td><td>G042713E43130N</td><td>9.24</td><td>-0.01 ± 0.28</td><td>-1.09 ± 0.27</td><td>-0.55 ± 0.15</td><td>-0.01 ± 0.24</td><td>-0.92 ± 0.25</td><td>-0.47 ± 0.13</td></tr><tr><td>Suatisi_Eastern</td><td>G044443E42706N</td><td>8.89</td><td>-0.11 ± 0.30</td><td>-0.86 ± 0.29</td><td>-0.48 ± 0.16</td><td>-0.09 ± 0.26</td><td>-0.73 ± 0.26</td><td>-0.41 ± 0.14</td></tr><tr><td>Khalde</td><td>G043043E42994N</td><td>8.59</td><td>-0.37 ± 0.31</td><td>0.11 ± 0.32</td><td>-0.13 ± 0.17</td><td>-0.31 ± 0.27</td><td>0.09 ± 0.27</td><td>-0.11 ± 0.15</td></tr><tr><td>Shkhelda</td><td>G042663E43151N</td><td>8.28</td><td>-0.44 ± 0.29</td><td>-1.13 ± 0.27</td><td>-0.79 ± 0.15</td><td>-0.38 ± 0.25</td><td>-0.96 ± 0.26</td><td>-0.67 ± 0.15</td></tr><tr><td>Midagrabin</td><td>G044390E42722N</td><td>8.16</td><td>-0.53 ± 0.31</td><td>-1.14 ± 0.28</td><td>-0.83 ± 0.16</td><td>-0.45 ± 0.27</td><td>-0.97 ± 0.26</td><td>-0.71 ± 0.15</td></tr><tr><td>Ullauzna</td><td>G043238E43064N</td><td>7.53</td><td>-0.09 ± 0.31</td><td>-0.27 ± 0.36</td><td>-0.18 ± 0.19</td><td>-0.08 ± 0.26</td><td>-0.23 ± 0.30</td><td>-0.16 ± 0.17</td></tr><tr><td>Bashil</td><td>G042901E43169N</td><td>7.35</td><td>-0.25 ± 0.31</td><td>-1.09 ± 0.28</td><td>-0.67 ± 0.16</td><td>-0.21 ± 0.27</td><td>-0.93 ± 0.25</td><td>-0.57 ± 0.14</td></tr><tr><td>Kvishi</td><td>G042497E43168N</td><td>7.33</td><td>0.04 ± 0.28</td><td>-1.03 ± 0.28</td><td>-0.50 ± 0.15</td><td>0.03 ± 0.24</td><td>-0.88 ± 0.26</td><td>-0.42 ± 0.13</td></tr><tr><td>Hrumkol</td><td>G043202E43023N</td><td>7.01</td><td>0.04 ± 0.30</td><td>-0.50 ± 0.32</td><td>-0.23 ± 0.18</td><td>0.03 ± 0.26</td><td>-0.42 ± 0.28</td><td>-0.20 ± 0.15</td></tr><tr><td>Mikelchiran</td><td>G042487E43376N</td><td>6.85</td><td>-0.15 ± 0.27</td><td>-0.44 ± 0.25</td><td>-0.30 ± 0.14</td><td>-0.13 ± 0.23</td><td>-0.38 ± 0.21</td><td>-0.25 ± 0.12</td></tr><tr><td>Psigansu</td><td>G043482E42972N</td><td>6.76</td><td>-0.26 ± 0.31</td><td>-0.90 ± 0.30</td><td>-0.58 ± 0.17</td><td>-0.22 ± 0.26</td><td>-0.77 ± 0.27</td><td>-0.49 ± 0.16</td></tr><tr><td>Maili</td><td>G044477E42719N</td><td>6.57</td><td>-0.19 ± 0.31</td><td>-0.44 ± 0.30</td><td>-0.31 ± 0.16</td><td>-0.16 ± 0.27</td><td>-0.38 ± 0.25</td><td>-0.27 ± 0.14</td></tr><tr><td>Djalovchatskiy</td><td>G041493E43296N</td><td>6.42</td><td>-1.07 ± 0.28</td><td>-2.01 ± 0.27</td><td>-1.54 ± 0.14</td><td>-0.91 ± 0.26</td><td>-1.71 ± 0.30</td><td>-1.31 ± 0.19</td></tr><tr><td>Bartuitsete</td><td>G043638E42825N</td><td>6.30</td><td>-0.65 ± 0.32</td><td>-0.81 ± 0.30</td><td>-0.73 ± 0.17</td><td>-0.56 ± 0.28</td><td>-0.68 ± 0.27</td><td>-0.62 ± 0.16</td></tr><tr><td>Lekhziri Central</td><td>G042751E43175N</td><td>6.27</td><td>0.09 ± 0.32</td><td>-0.88 ± 0.29</td><td>-0.39 ± 0.16</td><td>0.08 ± 0.27</td><td>-0.74 ± 0.26</td><td>-0.33 ± 0.14</td></tr><tr><td>Kulakchegemchiran</td><td>G042967E43132N</td><td>6.00</td><td>-0.41 ± 0.35</td><td>-1.12 ± 0.32</td><td>-0.77 ± 0.18</td><td>-0.35 ± 0.30</td><td>-0.95 ± 0.30</td><td>-0.65 ± 0.17</td></tr><tr><td>Karachaul</td><td>G042453E43381N</td><td>5.94</td><td>-0.07 ± 0.28</td><td>-0.30 ± 0.26</td><td>-0.19 ± 0.14</td><td>-0.06 ± 0.23</td><td>-0.26 ± 0.22</td><td>-0.16 ± 0.12</td></tr><tr><td>Gergeti</td><td>G044499E42692N</td><td>5.77</td><td>-0.03 ± 0.32</td><td>-0.65 ± 0.30</td><td>-0.34 ± 0.16</td><td>-0.03 ± 0.27</td><td>-0.55 ± 0.26</td><td>-0.29 ± 0.14</td></tr><tr><td>Alibekskiy</td><td>G041523E43281N</td><td>5.73</td><td>-1.12 ± 0.28</td><td>-1.30 ± 0.28</td><td>-1.21 ± 0.14</td><td>-0.95 ± 0.26</td><td>-1.10 ± 0.27</td><td>-1.03 ± 0.16</td></tr><tr><td>Kyukyurtlyu</td><td>G042404E43354N</td><td>5.69</td><td>0.07 ± 0.28</td><td>-0.24 ± 0.27</td><td>-0.08 ± 0.15</td><td>0.06 ± 0.24</td><td>-0.20 ± 0.23</td><td>-0.07 ± 0.13</td></tr><tr><td>Songutidontsete</td><td>G043809E42798N</td><td>5.66</td><td>-0.55 ± 0.34</td><td>-0.84 ± 0.32</td><td>-0.69 ± 0.18</td><td>-0.46 ± 0.30</td><td>-0.71 ± 0.29</td><td>-0.59 ± 0.17</td></tr><tr><td>Dolra</td><td>G042522E43171N</td><td>5.36</td><td>0.03 ± 0.29</td><td>-1.12 ± 0.30</td><td>-0.55 ± 0.16</td><td>0.02 ± 0.25</td><td>-0.95 ± 0.27</td><td>-0.47 ± 0.14</td></tr><tr><td>Ullukol</td><td>G042464E43376N</td><td>5.35</td><td>-0.14 ± 0.28</td><td>-0.37 ± 0.26</td><td>-0.26 ± 0.14</td><td>-0.12 ± 0.24</td><td>-0.32 ± 0.22</td><td>-0.22 ± 0.12</td></tr><tr><td>493v</td><td>G042811E43183N</td><td>5.32</td><td>-0.42 ± 0.34</td><td>-0.86 ± 0.31</td><td>-0.64 ± 0.17</td><td>-0.35 ± 0.29</td><td>-0.73 ± 0.28</td><td>-0.54 ± 0.16</td></tr><tr><td colspan="9">Glaciers with 2–5 km² area</td></tr><tr><td>Ushba Southern</td><td>G042633E43131N</td><td>4.72</td><td>-0.43 ± 0.33</td><td>-1.09 ± 0.32</td><td>-0.76 ± 0.17</td><td>-0.36 ± 0.28</td><td>-0.92 ± 0.29</td><td>-0.64 ± 0.16</td></tr><tr><td>Seri</td><td>G042838E43120N</td><td>4.65</td><td>-0.25 ± 0.34</td><td>-1.17 ± 0.32</td><td>-0.71 ± 0.18</td><td>-0.22 ± 0.29</td><td>-1.00 ± 0.30</td><td>-0.61 ± 0.17</td></tr><tr><td>Amanauz</td><td>G041626E43230N</td><td>4.63</td><td>-0.75 ± 0.29</td><td>-1.46 ± 0.29</td><td>-1.10 ± 0.15</td><td>-0.63 ± 0.26</td><td>-1.24 ± 0.29</td><td>-0.94 ± 0.17</td></tr></tbody></table>
<table><caption>Table A3. *Cont.*</caption><thead><tr><th rowspan="2">Glacier Name</th><th rowspan="2">GLIMS_id</th><th rowspan="2">Area (km²)</th><th colspan="3">Change of Glacier Thickness (m a⁻¹)<br>Hugonnet et al. (2021) [3]</th><th colspan="3">Change in Glacier Surface Mass Balance (m w.e. a⁻¹)<br>Hugonnet et al. (2021) [3]</th></tr><tr><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th></tr></thead><tbody><tr><td>Boko</td><td>G043719E42767N</td><td>4.62</td><td>-0.03 ± 0.34</td><td>-0.60 ± 0.33</td><td>-0.31 ± 0.19</td><td>-0.03 ± 0.29</td><td>-0.51 ± 0.28</td><td>-0.27 ± 0.16</td></tr><tr><td>Nageba</td><td>G042959E43036N</td><td>4.60</td><td>-0.16 ± 0.32</td><td>-0.24 ± 0.33</td><td>-0.20 ± 0.18</td><td>-0.14 ± 0.27</td><td>-0.21 ± 0.28</td><td>-0.17 ± 0.16</td></tr><tr><td>Kirtisho</td><td>G043596E42832N</td><td>4.51</td><td>-0.71 ± 0.33</td><td>-0.97 ± 0.31</td><td>-0.84 ± 0.18</td><td>-0.60 ± 0.29</td><td>-0.82 ± 0.28</td><td>-0.71 ± 0.17</td></tr><tr><td>Devdoraki</td><td>G044517E42715N</td><td>4.41</td><td>0.03 ± 0.32</td><td>-0.30 ± 0.30</td><td>-0.13 ± 0.16</td><td>0.02 ± 0.28</td><td>-0.25 ± 0.26</td><td>-0.11 ± 0.14</td></tr><tr><td>Tviberi</td><td>G042838E43099N</td><td>4.33</td><td>-0.30 ± 0.32</td><td>-1.46 ± 0.31</td><td>-0.88 ± 0.18</td><td>-0.26 ± 0.27</td><td>-1.24 ± 0.31</td><td>-0.75 ± 0.18</td></tr><tr><td>Bashkara</td><td>G042734E43194N</td><td>4.05</td><td>-0.48 ± 0.35</td><td>-1.02 ± 0.33</td><td>-0.75 ± 0.18</td><td>-0.41 ± 0.31</td><td>-0.87 ± 0.30</td><td>-0.64 ± 0.18</td></tr><tr><td>Sopruju Southern</td><td>G041572E43240N</td><td>3.96</td><td>-0.32 ± 0.29</td><td>-1.09 ± 0.29</td><td>-0.70 ± 0.14</td><td>-0.27 ± 0.25</td><td>-0.92 ± 0.27</td><td>-0.60 ± 0.14</td></tr><tr><td>Belalakaiskiy</td><td>G041561E4266N</td><td>3.88</td><td>-0.74 ± 0.30</td><td>-1.52 ± 0.29</td><td>-1.13 ± 0.15</td><td>-0.63 ± 0.26</td><td>-1.29 ± 0.30</td><td>-0.96 ± 0.18</td></tr><tr><td>Tanantsete</td><td>G043500E42888N</td><td>3.88</td><td>-0.53 ± 0.34</td><td>-0.23 ± 0.34</td><td>-0.38 ± 0.19</td><td>-0.45 ± 0.29</td><td>-0.20 ± 0.29</td><td>-0.32 ± 0.16</td></tr><tr><td>Djailik</td><td>G042922E43239N</td><td>3.86</td><td>-0.49 ± 0.38</td><td>-0.73 ± 0.35</td><td>-0.61 ± 0.19</td><td>-0.42 ± 0.32</td><td>-0.62 ± 0.30</td><td>-0.52 ± 0.17</td></tr><tr><td>Sugan</td><td>G043454E42991N</td><td>3.84</td><td>-0.42 ± 0.33</td><td>-0.76 ± 0.32</td><td>-0.59 ± 0.18</td><td>-0.36 ± 0.29</td><td>-0.65 ± 0.28</td><td>-0.50 ± 0.16</td></tr><tr><td>Tyutyun</td><td>G043249E43041N</td><td>3.79</td><td>-0.41 ± 0.35</td><td>-0.76 ± 0.41</td><td>-0.59 ± 0.22</td><td>-0.35 ± 0.30</td><td>-0.65 ± 0.36</td><td>-0.50 ± 0.20</td></tr><tr><td>Laila</td><td>G042534E42928N</td><td>3.73</td><td>-0.08 ± 0.30</td><td>-1.12 ± 0.29</td><td>-0.60 ± 0.16</td><td>-0.07 ± 0.26</td><td>-0.95 ± 0.27</td><td>-0.51 ± 0.14</td></tr><tr><td>Zaramag</td><td>G043845E42727N</td><td>3.67</td><td>-0.61 ± 0.36</td><td>-0.80 ± 0.33</td><td>-0.70 ± 0.19</td><td>-0.51 ± 0.31</td><td>-0.68 ± 0.29</td><td>-0.60 ± 0.18</td></tr><tr><td>Edena Bigger</td><td>G043380E42886N</td><td>3.66</td><td>-0.52 ± 0.33</td><td>-0.86 ± 0.31</td><td>-0.69 ± 0.18</td><td>-0.44 ± 0.29</td><td>-0.73 ± 0.28</td><td>-0.59 ± 0.17</td></tr><tr><td>Shkhara</td><td>G043106E42983N</td><td>3.61</td><td>-0.54 ± 0.37</td><td>-0.08 ± 0.37</td><td>-0.31 ± 0.20</td><td>-0.46 ± 0.32</td><td>-0.07 ± 0.32</td><td>-0.26 ± 0.17</td></tr><tr><td>Tanantsete</td><td>G043521E42874N</td><td>3.55</td><td>-0.80 ± 0.35</td><td>-0.83 ± 0.34</td><td>-0.81 ± 0.19</td><td>-0.68 ± 0.31</td><td>-0.70 ± 0.30</td><td>-0.69 ± 0.18</td></tr><tr><td>Gerti</td><td>G043276E43053N</td><td>3.50</td><td>-0.83 ± 0.34</td><td>-0.53 ± 0.40</td><td>-0.68 ± 0.22</td><td>-0.70 ± 0.30</td><td>-0.45 ± 0.35</td><td>-0.58 ± 0.19</td></tr><tr><td>Leadashti Northern</td><td>G042479E43178N</td><td>3.46</td><td>-0.17 ± 0.30</td><td>-1.08 ± 0.30</td><td>-0.63 ± 0.16</td><td>-0.15 ± 0.26</td><td>-0.92 ± 0.27</td><td>-0.53 ± 0.14</td></tr><tr><td>Khasautskiy</td><td>G041439E43313N</td><td>3.32</td><td>-0.95 ± 0.40</td><td>-1.47 ± 0.39</td><td>-1.21 ± 0.19</td><td>-0.81 ± 0.35</td><td>-1.25 ± 0.35</td><td>-1.03 ± 0.20</td></tr><tr><td>Shtulu Western</td><td>G043371E42899N</td><td>3.30</td><td>-0.79 ± 0.34</td><td>-0.81 ± 0.32</td><td>-0.80 ± 0.18</td><td>-0.67 ± 0.30</td><td>-0.69 ± 0.28</td><td>-0.68 ± 0.18</td></tr><tr><td>Buba</td><td>G043762E42752N</td><td>3.26</td><td>-0.72 ± 0.36</td><td>-0.59 ± 0.33</td><td>-0.65 ± 0.19</td><td>-0.61 ± 0.32</td><td>-0.50 ± 0.29</td><td>-0.55 ± 0.18</td></tr><tr><td>Mostotsete</td><td>G043421E42899N</td><td>3.24</td><td>-1.09 ± 0.34</td><td>-1.07 ± 0.33</td><td>-1.08 ± 0.19</td><td>-0.93 ± 0.31</td><td>-0.91 ± 0.30</td><td>-0.92 ± 0.20</td></tr><tr><td>Morde-auzh</td><td>G042112E43216N</td><td>3.07</td><td>-0.50 ± 0.29</td><td>-0.98 ± 0.31</td><td>-0.74 ± 0.16</td><td>-0.43 ± 0.25</td><td>-0.84 ± 0.28</td><td>-0.63 ± 0.16</td></tr><tr><td>Sopruju Northern</td><td>G041588E43250N</td><td>3.06</td><td>-0.44 ± 0.31</td><td>-0.95 ± 0.31</td><td>-0.70 ± 0.16</td><td>-0.37 ± 0.27</td><td>-0.81 ± 0.28</td><td>-0.59 ± 0.15</td></tr><tr><td>Kashkatash</td><td>G042693E43189N</td><td>3.05</td><td>-0.18 ± 0.34</td><td>-0.79 ± 0.32</td><td>-0.48 ± 0.17</td><td>-0.15 ± 0.29</td><td>-0.67 ± 0.29</td><td>-0.41 ± 0.16</td></tr><tr><td>Donguz Orun</td><td>G042509E43213N</td><td>3.01</td><td>-0.18 ± 0.38</td><td>-0.74 ± 0.37</td><td>-0.46 ± 0.19</td><td>-0.15 ± 0.32</td><td>-0.63 ± 0.32</td><td>-0.39 ± 0.17</td></tr><tr><td>Gvandra</td><td>G042078E43218N</td><td>2.95</td><td>-0.50 ± 0.29</td><td>-1.00 ± 0.34</td><td>-0.75 ± 0.18</td><td>-0.43 ± 0.26</td><td>-0.85 ± 0.31</td><td>-0.64 ± 0.17</td></tr><tr><td>Ushba Northern</td><td>G042613E43151N</td><td>2.89</td><td>0.06 ± 0.32</td><td>-0.80 ± 0.32</td><td>-0.37 ± 0.17</td><td>0.05 ± 0.27</td><td>-0.68 ± 0.28</td><td>-0.31 ± 0.15</td></tr><tr><td>Marukhi Northern</td><td>G041420E43362N</td><td>2.82</td><td>-0.52 ± 0.50</td><td>-1.60 ± 0.50</td><td>-1.06 ± 0.24</td><td>-0.44 ± 0.43</td><td>-1.36 ± 0.45</td><td>-0.90 ± 0.23</td></tr></tbody></table>
<table><caption>Table A3. *Cont.*</caption><thead><tr><th rowspan="2">Glacier Name</th><th rowspan="2">GLIMS_id</th><th rowspan="2">Area (km²)</th><th colspan="3">Change of Glacier Thickness (m a⁻¹)<br>Hugonnet et al. (2021) [3]</th><th colspan="3">Change in Glacier Surface Mass Balance (m w.e. a⁻¹)<br>Hugonnet et al. (2021) [3]</th></tr><tr><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th><th>2000–2010</th><th>2010–2019</th><th>2000–2019</th></tr></thead><tbody><tr><td>Lakra</td><td>G042468E43152N</td><td>2.69</td><td>-0.38 ± 0.32</td><td>-1.27 ± 0.32</td><td>-0.82 ± 0.16</td><td>-0.32 ± 0.27</td><td>-1.08 ± 0.30</td><td>-0.70 ± 0.17</td></tr><tr><td>Shtulu Eastern</td><td>G043395E42903N</td><td>2.67</td><td>-0.61 ± 0.34</td><td>-0.63 ± 0.33</td><td>-0.62 ± 0.19</td><td>-0.52 ± 0.29</td><td>-0.54 ± 0.29</td><td>-0.53 ± 0.17</td></tr><tr><td>Chalinginitsete</td><td>G043785E42817N</td><td>2.65</td><td>-1.01 ± 0.35</td><td>-1.03 ± 0.34</td><td>-1.02 ± 0.19</td><td>-0.86 ± 0.32</td><td>-0.87 ± 0.31</td><td>-0.87 ± 0.20</td></tr><tr><td>Tebulo</td><td>G045335E42577N</td><td>2.62</td><td>-0.25 ± 0.33</td><td>-0.25 ± 0.34</td><td>-0.25 ± 0.19</td><td>-0.21 ± 0.28</td><td>-0.21 ± 0.29</td><td>-0.21 ± 0.16</td></tr><tr><td>Tiutiurgy</td><td>G043081E43122N</td><td>2.62</td><td>-0.16 ± 0.34</td><td>-0.61 ± 0.34</td><td>-0.38 ± 0.19</td><td>-0.14 ± 0.29</td><td>-0.52 ± 0.30</td><td>-0.33 ± 0.16</td></tr><tr><td>Dvuyazichniy</td><td>G041519E43299N</td><td>2.62</td><td>-0.68 ± 0.31</td><td>-1.05 ± 0.31</td><td>-0.86 ± 0.16</td><td>-0.58 ± 0.27</td><td>-0.89 ± 0.28</td><td>-0.73 ± 0.15</td></tr><tr><td>Kichkinekol Boshoy</td><td>G042177E43228N</td><td>2.60</td><td>-0.69 ± 0.33</td><td>-0.99 ± 0.34</td><td>-0.84 ± 0.18</td><td>-0.59 ± 0.29</td><td>-0.84 ± 0.32</td><td>-0.71 ± 0.18</td></tr><tr><td>Ezengi</td><td>G042534E43202N</td><td>2.59</td><td>-0.02 ± 0.32</td><td>-1.13 ± 0.30</td><td>-0.58 ± 0.16</td><td>-0.02 ± 0.27</td><td>-0.96 ± 0.29</td><td>-0.49 ± 0.16</td></tr><tr><td>Tiutiu</td><td>G042879E43260N</td><td>2.56</td><td>-0.29 ± 0.71</td><td>-0.73 ± 0.70</td><td>-0.51 ± 0.36</td><td>-0.24 ± 0.60</td><td>-0.62 ± 0.61</td><td>-0.43 ± 0.31</td></tr><tr><td>Kolka</td><td>G044447E42733N</td><td>2.51</td><td>-0.83 ± 0.38</td><td>0.53 ± 0.35</td><td>-0.15 ± 0.19</td><td>-0.71 ± 0.34</td><td>0.45 ± 0.30</td><td>-0.13 ± 0.16</td></tr><tr><td>Nakhagibita</td><td>G043451E42962N</td><td>2.49</td><td>0.00 ± 0.35</td><td>-0.92 ± 0.34</td><td>-0.46 ± 0.19</td><td>0.00 ± 0.30</td><td>-0.78 ± 0.30</td><td>-0.39 ± 0.17</td></tr><tr><td>Karaugom</td><td>G043704E42797N</td><td>2.48</td><td>-0.78 ± 0.37</td><td>-0.41 ± 0.35</td><td>-0.59 ± 0.20</td><td>-0.66 ± 0.32</td><td>-0.35 ± 0.30</td><td>-0.50 ± 0.18</td></tr><tr><td>Zopkhito</td><td>G043422E42884N</td><td>2.45</td><td>-0.48 ± 0.35</td><td>-0.85 ± 0.34</td><td>-0.67 ± 0.20</td><td>-0.41 ± 0.30</td><td>-0.72 ± 0.30</td><td>-0.57 ± 0.18</td></tr><tr><td>Lailchala</td><td>G042581E42915N</td><td>2.38</td><td>-0.07 ± 0.33</td><td>-0.93 ± 0.32</td><td>-0.50 ± 0.17</td><td>-0.06 ± 0.28</td><td>-0.79 ± 0.29</td><td>-0.42 ± 0.15</td></tr><tr><td>Kharikhra</td><td>G042093E43100N</td><td>2.36</td><td>-0.42 ± 0.30</td><td>-0.47 ± 0.31</td><td>-0.45 ± 0.16</td><td>-0.36 ± 0.26</td><td>-0.40 ± 0.27</td><td>-0.38 ± 0.14</td></tr><tr><td>Garabashi</td><td>G042470E43302N</td><td>2.33</td><td>-0.45 ± 0.31</td><td>-0.61 ± 0.31</td><td>-0.53 ± 0.16</td><td>-0.38 ± 0.27</td><td>-0.52 ± 0.27</td><td>-0.45 ± 0.15</td></tr><tr><td>Mna</td><td>G044483E42698N</td><td>2.32</td><td>0.26 ± 0.34</td><td>-0.54 ± 0.34</td><td>-0.14 ± 0.18</td><td>0.22 ± 0.29</td><td>-0.46 ± 0.29</td><td>-0.12 ± 0.15</td></tr><tr><td>Namkvami</td><td>G043128E42966N</td><td>2.29</td><td>-0.53 ± 0.35</td><td>-0.15 ± 0.38</td><td>-0.34 ± 0.20</td><td>-0.45 ± 0.30</td><td>-0.12 ± 0.32</td><td>-0.29 ± 0.18</td></tr><tr><td>Koruldashi</td><td>G043157E42956N</td><td>2.29</td><td>-0.38 ± 0.36</td><td>-0.51 ± 0.39</td><td>-0.44 ± 0.21</td><td>-0.32 ± 0.31</td><td>-0.43 ± 0.34</td><td>-0.37 ± 0.19</td></tr><tr><td>Agbek</td><td>G041673E43238N</td><td>2.28</td><td>-0.72 ± 0.34</td><td>-0.82 ± 0.33</td><td>-0.77 ± 0.17</td><td>-0.61 ± 0.30</td><td>-0.70 ± 0.30</td><td>-0.65 ± 0.18</td></tr><tr><td>Shdavleri</td><td>G042337E43141N</td><td>2.28</td><td>-0.46 ± 0.32</td><td>-0.81 ± 0.32</td><td>-0.63 ± 0.17</td><td>-0.39 ± 0.28</td><td>-0.69 ± 0.28</td><td>-0.54 ± 0.16</td></tr><tr><td>Iriti</td><td>G042877E43141N</td><td>2.27</td><td>0.12 ± 0.42</td><td>-0.93 ± 0.40</td><td>-0.40 ± 0.21</td><td>0.10 ± 0.35</td><td>-0.79 ± 0.35</td><td>-0.34 ± 0.18</td></tr><tr><td>Djankuat</td><td>G042766E43192N</td><td>2.23</td><td>-0.47 ± 0.34</td><td>-1.31 ± 0.30</td><td>-0.89 ± 0.17</td><td>-0.40 ± 0.30</td><td>-1.11 ± 0.29</td><td>-0.76 ± 0.17</td></tr><tr><td>Khakel</td><td>G041843E43225N</td><td>2.22</td><td>-0.77 ± 0.30</td><td>-1.71 ± 0.31</td><td>-1.24 ± 0.16</td><td>-0.66 ± 0.28</td><td>-1.45 ± 0.34</td><td>-1.06 ± 0.20</td></tr><tr><td>Rtsivashki</td><td>G043434E42978N</td><td>2.13</td><td>-0.46 ± 0.36</td><td>-0.54 ± 0.35</td><td>-0.50 ± 0.19</td><td>-0.39 ± 0.31</td><td>-0.46 ± 0.30</td><td>-0.42 ± 0.17</td></tr><tr><td>Khasankoi-siuriulgen</td><td>G042386E43249N</td><td>2.08</td><td>-0.71 ± 0.32</td><td>-1.20 ± 0.31</td><td>-0.96 ± 0.16</td><td>-0.61 ± 0.28</td><td>-1.02 ± 0.29</td><td>-0.81 ± 0.17</td></tr><tr><td>Murkvami</td><td>G042813E43111N</td><td>2.08</td><td>-0.37 ± 0.35</td><td>-0.50 ± 0.32</td><td>-0.44 ± 0.18</td><td>-0.32 ± 0.30</td><td>-0.43 ± 0.28</td><td>-0.37 ± 0.16</td></tr><tr><td>Suatisi Central</td><td>G044416E42704N</td><td>2.07</td><td>0.41 ± 0.36</td><td>-0.75 ± 0.35</td><td>-0.17 ± 0.19</td><td>0.35 ± 0.31</td><td>-0.64 ± 0.31</td><td>-0.15 ± 0.16</td></tr><tr><td>Phastagtsete</td><td>G043655E42811N</td><td>2.06</td><td>-1.21 ± 0.37</td><td>-1.17 ± 0.36</td><td>-1.19 ± 0.20</td><td>-1.03 ± 0.35</td><td>-1.00 ± 0.34</td><td>-1.01 ± 0.23</td></tr><tr><td>Bulungy</td><td>G043105E43134N</td><td>2.02</td><td>-0.38 ± 0.34</td><td>-0.50 ± 0.39</td><td>-0.44 ± 0.21</td><td>-0.32 ± 0.29</td><td>-0.43 ± 0.34</td><td>-0.37 ± 0.19</td></tr></tbody></table>
<figure><img src="image_19.png" /></figure>

**Figure A1.** Change in in glacier thickness or surface elevation (m a⁻¹) as a function of glacier mean elevation (m a.s.l.) in the Greater Caucasus between 2000 and 2019. Elevation bands in 200 m. Vertical stems show an uncertainty estimation of the mean values.

<figure><img src="image_20.png" /></figure>

**Figure A2.** Changes in glacier thickness or surface elevation (m a⁻¹) for all Caucasus glaciers according to the regions (western, central, and eastern), slopes (northern, southern) (a), and countries (b) in 2000–2010, 2010–2019, and 2000–2019. Vertical stems showing an uncertainty estimation of the mean values.

## References

1. Cogley, J.G.; Hock, R.; Rasmussen, L.A.; Arendt, A.A.; Bauder, A.; Braithwaite, R.J.; Jansson, P.; Kaser, G.; Möller, M.; Nicholson, L. *Glossary of Glacier Mass Balance and Related Terms*; IHP-VII Technical Documents in Hydrology No. 86; IACS Contribution No. 2; UNESCO-IHP: Paris, France, 2011.

2. Huss, M.; Hock, R. A new model for global glacier change and sea-level rise. *Front. Earth Sci.* **2015**, *3*, 54. [CrossRef]

3. Hugonnet, R.; McNabb, R.; Berthier, E.; Menounos, B.; Nuth, C.; Girod, L.; Farinotti, D.; Huss, M.; Dussaillant, I.; Brun, F.; et al. Accelerated global glacier mass loss in the early twenty-first century. *Nature* **2021**, *592*, 726–731. [CrossRef] [PubMed]

4. Zemp, M.; Huss, M.; Thibert, E.; Eckert, N.; McNabb, R.; Huber, J.; Barandun, M.; Machguth, H.; Nussbaumer, S.U.; Gärtner-Roer, I.; et al. Global glacier mass changes and their contributions to sea-level rise from 1961 to 2016. *Nature* **2019**, *568*, 382–386. [CrossRef]

5. Hock, R.; Rasul, G.; Adler, C.; Cáceres, B.; Gruber, S.; Hirabayashi, Y.; Jackson, M.; Kääb, A.; Kang, S.; Kutuzov, S.; et al. High Mountain Areas. In *IPCC Special Report on Ocean and the Cryosphere in a Changing Climate (SROCC)*; IPCC: Geneva, Switzerland, 2019.
6. Nuth, C.; Kääb, A. Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change. *Cryosphere* **2011**, *5*, 271–290. [CrossRef]

7. Fischer, M.; Huss, M.; Kummert, M.; Hoelzle, M. Application and validation of long-range terrestrial laser scanning to monitor the mass balance of very small glaciers in the Swiss Alps. *Cryosphere* **2016**, *10*, 1279–1295. [CrossRef]

8. Klug, C.; Bollmann, E.; Galos, S.P.; Nicholson, L.; Prinz, R.; Rieg, L.; Sailer, R.; Stötter, J.; Kaser, G. Geodetic reanalysis of annual glaciological mass balances (2001–2011) of Hintereisferner, Austria. *Cryosphere* **2018**, *12*, 833–849. [CrossRef]

9. Basantes-Serrano, R.; Rabatel, A.; Vincent, C.; And Sirguey, P. An optimized method to calculate the geodetic mass balance of mountain glaciers. *J. Glaciol.* **2018**, *64*, 917–931. [CrossRef]

10. Vincent, C.; Cusicanqui, D.; Jourdain, B.; Laarman, O.; Six, D.; Gilbert, A.; Walpersdorf, A.; Rabatel, A.; Piard, L.; Gimbert, F.; et al. Geodetic point surface mass balances: A new approach to determine point surface mass balances on glaciers from remote sensing measurements. *Cryosphere* **2021**, *15*, 1259–1276. [CrossRef]

11. Evans, S.G.; Tutubalina, O.V.; Drobyshev, V.N.; Chernomorets, S.S.; McDougall, S.; Petrakov, D.A.; Hungr, O. Catastrophic detachment and high-velocity long-runout flow of Kolka Glacier, Caucasus Mountains, Russia in 2002. *Geomorphology* **2009**, *105*, 314–321. [CrossRef]

12. Chernomorets, S.S.; Petrakov, D.A.; Aleynikov, A.A.; Bekkiev, M.Y.; Viskhadzheva, K.S.; Dokukin, M.D.; Kalov, R.K.; Kidyaeva, V.M.; Krylenko, V.V.; Krylenko, I.V.; et al. The Outburst of Bashkara Glacier Lake (Central Caucasus, Russia) on 1 September 2017. *Kriosf. Zemli* **2018**, *XXII*, 61–70. [CrossRef]

13. Tielidze, L.G.; Kumladze, R.M.; Wheate, R.D.; Gamkrelidze, M. The Devdoraki Glacier Catastrophes, Georgian Caucasus. *Hung. Geogr. Bull.* **2019**, *68*, 21–35. [CrossRef]

14. Tsereteli, D.; Khazaradze, R.; Lomtatidze, G.; Inashvili, S.; Lashkhi, T.; Kurdghelaidze, G.; Kalandadze, G.; Chekurishvili, R. Glaciological observations on the Chalaati and Lechziri Glaciers (Upper Svaneti) in the spring of 1959. *Georgian Natl. Acad. Sci. —Work. Vakhushti Inst. Geogr.* **1962**, *XVIII*, 223–256. (In Georgian)

15. Menshutin, V.M. Mass balance of the Marukh Glacier (Western Caucasus) 1966–1968. *Data Glaciol. Stud.* **1970**, *17*, 292–296. (In Russian)

16. Boyarsky, I.Y. (Ed.) *Lednik Djankuat (Djankuat Glacier)*; Gidrometeoizdat: Saint Petersburg, Russia, 1978.

17. Rototaeva, O.V.; Nosenko, G.A.; Kerimov, A.M.; Kutuzov, S.S.; Lavrentiev, I.I.; Nikitin, S.A.; Kerimov, A.A.; Tarasova, L.N. Changes of the mass balance of the Garabashy Glacier, Mount Elbrus, at the turn of 20th and 21st centuries. *Ice Snow* **2019**, *59*, 5–22. [CrossRef]

18. Gigineishvili, G. (Ed.) *Tbilisa Glacier—Glacial-Hydro-Meteorological Investigation. Monograph*; Institute of Geography at Georgian National Academy of Sciences: Tbilisi, Georgia; Publishing house Metsniereba: Tbilisi, Georgia, 1986; p. 116. (In Russian)

19. Tielidze, L.G.; Svanadze, D.; Gadrani, L.; Asanidze, L.; Wheate, R.D.; Hamilton, G.S. A 54-year record of changes at Chalaati and Zopkhito glaciers, Georgian Caucasus, observed from archival maps, satellite imagery, drone survey and ground-based investigation. *Hung. Geogr. Bull.* **2020**, *69*, 175–189. [CrossRef]

20. WGMS; Zemp, M.; Nussbaumer, S.U.; GärtnerRoer, I.; Bannwart, J.; Paul, F.; Hoelzle, M. (Eds.) *Global Glacier Change Bulletin No. 4* (2018–2019); ISC (WDS)/IUGG (IACS)/UNEP/UNESCO/WMO; World Glacier Monitoring Service: Zurich, Switzerland, 2021; p. 278. [CrossRef]

21. Tielidze, L.G.; Wheate, R.D. The Greater Caucasus Glacier Inventory (Russia, Georgia and Azerbaijan). *Cryosphere* **2018**, *12*, 81–94. [CrossRef]

22. Lur'e, P.M.; Panov, V.D. Variations of Contemporary Glaciation of the Northern Slope of the Greater Caucasus in the 20th Century and the Forecast of Its Degradation in the 21st Century. *Russ. Meteorol. Hydrol.* **2014**, *39*, 254–259. [CrossRef]

23. Kotlyakov, V.M.; Khromova, T.E.; Nosenko, G.A.; Popova, V.V.; Chernova, L.P.; Murav'ev, A.Y. New Data on Current Changes in the Mountain Glaciers of Russia. In *Doklady Earth Sciences*; Pleiades Publishing: Warrensburg, MO, USA, 2015; Volume 464. [CrossRef]

24. Tielidze, L.G. Glacier change over the last century, Caucasus Mountains, Georgia, observed from old topographical maps, Landsat and ASTER satellite imagery. *Cryosphere* **2016**, *10*, 713–725. [CrossRef]

25. Rets, E.P.; Popovnin, V.V.; Toropov, P.A.; Smirnov, A.M.; Tokarev, I.V.; Chizhova, J.N.; Budantseva, N.A.; Vasil'chuk, Y.K.; Kireeva, M.B.; Ekaykin, A.A.; et al. Djankuat glacier station in the North Caucasus, Russia: A database of glaciological, hydrological, and meteorological observations and stable isotope sampling results during 2007–2017. *Earth Syst. Sci. Data* **2019**, *11*, 1463–1481. [CrossRef]

26. Zolotarev, E.A.; Kharkovets, E.G. Development of glaciers of Mount Elbrus after the Little Ice Age. *Ice Snow* **2012**, *52*, 15–22. [CrossRef]

27. Kutuzov, S.; Lavrentiev, I.; Smirnov, A.; Nosenko, G.; Petrakov, D. Volume Changes of Elbrus Glaciers From 1997 to 2017. *Front. Earth Sci.* **2019**, *7*, 153. [CrossRef]

28. Moussavi, M.; Zoej, M.; Vaziri, F.; Sahebi, M.; Rezaei, Y. A new glacier inventory of Iran. *Ann. Glaciol.* **2009**, *50*, 93–103. [CrossRef]

29. Azzoni, R.S.; Sarıkaya, M.A.; Fugazza, D. Turkish glacier inventory and classification from high-resolution satellite data. *Mediterr. Geosci. Rev.* **2020**, *2*, 153–162. [CrossRef]

30. Olawuyi, D. *Climate Change Law and Policy in the Middle East and North Africa Region*, 1st ed.; Routledge: London, UK, 2021; pp. 1–340. ISBN 9780367490324.
31. *World Atlas of Snow and Ice Resources*, V. 1.; (V. Kotlyakov editor) Russian Academy of Sciences: Moscow, Russian, 1997; 372p. (In Russian/English)

32. Volodicheva, N. The Caucasus. In *The Physical Geography of Northern Eurasia*; Shahgedanova, M., Ed.; Oxford University Press: Oxford, UK, 2002; pp. 350–376.

33. Kozachek, A.; Mikhalenko, V.; Masson-Delmotte, V.; Ekaykin, A.; Ginot, P.; Kutuzov, S.; Legrand, M.; Lipenkov, V.; Preunkert, S. Large-scale drivers of Caucasus climate variability in meteorological records and Mt El'brus ice cores. *Clim. Past* **2017**, *13*, 473–489. [CrossRef]

34. Sylvén, M.; Reinvang, R.; Andersone-Lilley, Ž. Climate Change in Southern Caucasus: Impacts on nature, people and society (WWF overview report). In *WWF Norway WWF Caucasus Programme*; WWF Norway: Washington, DC, USA, 2008.

35. Solomina, O.; Bushueva, I.; Dolgova, E.; Jomelli, V.; Alexandrin, M.; Mikhalenko, V.; Matskovsky, V. Glacier variations in the Northern Caucasus compared to climatic reconstructions over the past millennium. *Glob. Planet. Chang.* **2016**, *140*, 28–58. [CrossRef]

36. Tielidze, L.G.; Solomina, O.N.; Jomelli, V.; Dolgova, E.A.; Bushueva, I.S.; Mikhalenko, V.N.; Brauche, R.; ASTER, T. Change of Chalaati Glacier (Georgian Caucasus) since the Little Ice Age based on dendrochronological and Beryllium-10 data. *Ice Snow* **2020**, *60*, 453–470. [CrossRef]

37. Tielidze, L.G.; Nosenko, G.A.; Khromova, T.E.; Paul, F. Strong acceleration of glacier area loss in the Greater Caucasus over the past two decades. *Cryosphere Discuss* **2021**, *8*, 1–21. [CrossRef]

38. Tielidze, L.G.; Bolch, T.; Wheate, R.D.; Kutuzov, S.S.; Lavrentiev, I.I.; Zemp, M. Supra-glacial debris cover changes in the Greater Caucasus from 1986 to 2014. *Cryosphere* **2020**, *14*, 585–598. [CrossRef]

39. Rau, F.; Mauz, F.; Vogt, S.; Singh Khalsa, S.J.; Raup, B. *Illustrated GLIMS Glacier Classification Manual, Glacier Classification Guidance for the GLIMS Glacier Inventory*; 2005.

40. Benn, D.; Bolch, T.; Hands, K.; Gulley, J.; Luckman, A.; Nicholson, L.; Quincey, D.; Thompson, S.; Toumi, R.; Wiseman, S. Response of debris-covered glaciers in the Mount Everest region to recent warming, and implications for outburst flood hazards. *Earth-Sci. Rev.* **2012**, *114*, 156–174. [CrossRef]

41. Pratap, B.; Dobhal, D.; Mehta, M.; Bhambri, R. Influence of debris cover and altitude on glacier surface melting: A case study on Dokriani Glacier, central Himalaya, India. *Ann. Glaciol.* **2015**, *56*, 9–16. [CrossRef]

42. Brun, F.; Wagnon, P.; Berthier, E.; Shea, J.M.; Immerzeel, W.W.; Kraaijenbrink, P.D.A.; Vincent, C.; Reverchon, C.; Shrestha, D.; Arnaud, Y. Ice cliff contribution to the tongue-wide ablation of Changri Nup Glacier, Nepal, central Himalaya. *Cryosphere* **2018**, *12*, 3439–3457. [CrossRef]

43. Östrem, G. Ice Melting under a Thin Layer of Moraine, and the Existence of Ice Cores in Moraine Ridges. *Geogr. Ann.* **1959**, *41*, 228–230. [CrossRef]

44. Nicholson, L.I.; Benn, D.I. Calculating ice melt beneath a debris layer using meteorological data. *J. Glaciol.* **2006**, *52*, 463–470. [CrossRef]

45. Nicholson, L.; Benn, D.I. Properties of Natural Supraglacial Debris in Relation to Modelling Sub-debris Ice Ablation. *Earth Surf. Process. Landf.* **2013**, *38*, 490–501. [CrossRef]

46. Reid, T.; Brock, B. Assessing ice-cliff backwasting and its contribution to total ablation of debris-covered Miage glacier, Mont Blanc massif, Italy. *J. Glaciol.* **2014**, *60*, 3–13. [CrossRef]

47. Shahgedanova, M.; Stokes, C.R.; Gurney, S.D.; Popovnin, V.V. Interactions between mass balance, atmospheric circulation and recent climate change on the Djankuat Glacier, Caucasus Mountains, Russia. *J. Geophys. Res. Atmos.* **2005**, *110*, D04108. [CrossRef]

48. Toropov, P.A.; Aleshina, M.A.; Grachev, A.M. Large-scale climatic factors driving glacier recession in the Greater Caucasus, 20th–21st century. *Int. J. Climatol.* **2019**, *39*, 4703–4720. [CrossRef]

49. Toropov, P.A.; Mikhalenko, V.N.; Kutuzov, S.S.; Morozova, P.A.; Shestakova, A.A. Temperature and radiation regime of glaciers on slopes of the Mount Elbrus in the ablation period over last 65 years. *Ice Snow* **2016**, *1*, 5–19. [CrossRef]

50. Kutuzov, S.; Shahgedanova, M.; Mikhalenko, V.; Ginot, P.; Lavrentiev, I.; Kemp, S. High-resolution provenance of desert dust deposited on Mt. Elbrus, Caucasus in 2009–2012 using snow pit and firn core records. *Cryosphere* **2013**, *7*, 1481–1498. [CrossRef]

51. Dumont, M.; Tuzet, F.; Gascoin, S.; Picard, G.; Kutuzov, S.; Lafaysse, M.; Cluzet, B.; Nheili, R.; Painter, T.H. Accelerated snow melt in the Russian Caucasus mountains after the Saharan dust outbreak in March 2018. *J. Geophys. Res. Earth Surf.* **2020**, *125*, e2020JF005641. [CrossRef]