Isotopic Composition of Glacier Ice and Meltwater in the Arid Parts of the Altai Mountains (Central Asia)

: The objective of this study is to reveal the isotopic composition of ice and meltwater in glaciated regions of South-Eastern Altai. The paper depicts differences between the isotopic composition of glacier ice from several types of glaciers and from various locations. Detected differences between the isotopic composition of glacier ice in diversiﬁed parts of the study region are related to local climate patterns. Isotopic composition of meltwater and isotopic separation for glacier rivers runoff showed that in the Tavan-Bogd massif, seasonal snow participates more in the formation of glacier runoff due to better conditions for snow accumulation on the surface of glaciers. In other research areas pure glacier meltwater prevails in runoff.


Introduction
The use of environmental isotope tracers (such as isotopic composition of natural waters) for the solution to hydrological and glaciological problems is one of the leading directions in the development of Earth sciences. Isotope studies in high mountain regions may be used for further research of glaciers features and mountain river runoff.
The South-Eastern Altai is characterized by arid climate and is a problem area in terms of water resources. Due to the small amount of precipitation runoff from glaciers in the South-East, Altai plays a special role in the economic activity of the local population. The study of high mountain glaciers and rivers for water balance research is a cornerstone for water management estimates and forecasting, especially within the framework of global climate change.
To achieve this aim, the following tasks were solved: evaluation and comparison of the isotopic characteristics of the glacier's ice in different parts of South-Eastern Altai and the evaluation of the isotopic composition and its time variability for the water of glacial rivers.
In Russia, stable isotopes research of glacier river runoff and glacier ice has been carried out for several years in the Central Caucasus, using examples of runoff from the Dzhankuat glacier (mainly by the Moscow State University) [1][2][3].
Isotopic research in the Altai mountains in general demands deep ice drilling. Earlier the ice core was investigated on Belukha Mountain [4,5]. The ice core was also sampled at the Tsambagarav massif (Mongolia) [6]. Apart from this ice-water, isotopic features and relations in this study area are not well presented.
The isotopic composition of water has been studied thoroughly in areas closest to the Altai Mountain areas, which are located in China. Earlier studies have determined the prevalence of glacial runoff in the total runoff from July to August and the prevalence of ground feeding in winter for the Urumqi River [7]. Similar studies of the isotopic composition of glacial-originating rivers were also performed in the north of the Tibetan mountains in China. The authors determined the isotopic characteristics of the components, which form the runoff, including snow, firn and glacial ice, and evaluated the role of glaciation systems in the feeding of the studied rivers [8][9][10][11]. Another isotopic separation study was carried out in a glacierized region of Southwest China [12]. In this part of China, the Lancang River has also been the object of complex stable isotopes research [13].
As stated above, Russian and Mongolian parts of the South-Eastern Altai are not well-provided by glacier ice and glacier runoff isotopic research. Additionally, this kind of research can be very useful especially within the framework of global climate change and further deglaciation.

Materials and Methods
Sample collection was carried out in the middle of the ablation season, from July to August. Over the course of nine years, more than 800 samples were drawn. Glacial ice was collected from the surface of glaciers into the plastic bags, then melted at an ambient temperature and poured into test tubes. Samples from watercourses were drawn directly into test tubes. Precipitation samples were collected immediately after precipitation to avoid evaporation The relative isotopic content is indicated by δ. SMOW (standard mean ocean water) is taken as the zero standard.
The analysis of isotopic characteristics for samples collected in 2012-2017 was carried out at the Laboratory for Climate and Environmental Change at the Arctic and Antarctic Research Institute with Cavity Ring-Down Spectroscopy (CRDS) using a a Picarro L2120-i gas analyser. Distilled tap water from St. Petersburg was used as a laboratory standard; it had the following characteristics: −9.79‰ in δ 18 O and −75.47‰ in δ 2 H relative to the IAEA "V-SMOW2" standard. The measurement precision was 0.05‰ for δ 18 O and 0.5‰ for δ 2 H, which is quite enough for this kind of research. Isotopic composition was rounded off to the nearest 0.1. In 2018, analysis was fulfilled at the Resource Center "X-ray Diffraction Research Methods" of the St. Petersburg State University Science Park using the same spectrometer. USGS45, USGS46 and GISP standards were used with the same measurement accuracy.
The relationship between the relative content of 18 O and D in precipitation was determined by the fractionation coefficients of these isotopes. The formula of linear relationship between 18 O and D was obtained by H. Craig, and it has the following form: This linear relation is designated as global meteoric water line (GMWL) [14]. In 2005, this equation was refined using long-term means from GNIP stations and the resulting GMWL was [15]: Studies of runoff in areas with developed glaciation are based on the fact that using stable isotopes as tracers, it is possible to divide the hydrograph of glacial rivers into components. The use of stable isotopes in glaciology and hydrology is based on natural differences in the isotopic composition of the components of the glacial rivers' runoff. Isotopic methods allow to carry out only two-component separation, i.e., to calculate the shares of the two most significant components in the total runoff can be calculated. This separation is carried out using the isotope balance formula, which generally has the form: R 18

Isotopic Composition of Glacier Ice
Glaciers are one of the main water sources in arid South-Eastern Altai. To estimate their role in alimentation of mountain rivers, we required data on the isotopic composition of glacier ice. Sampling of glacier ice was carried out in three massifs: Tavan-Bogd, Mongun-Taiga in Russia and Tsambagarav in Mongolia ( Figure 1).
ponents. The use of stable isotopes in glaciology and hydrology is based on natural differences in the isotopic composition of the components of the glacial rivers' runoff. Isotopic methods allow to carry out only two-component separation, i.e., to calculate the shares of the two most significant components in the total runoff can be calculated.
This separation is carried out using the isotope balance formula, which generally has the form:

Isotopic Composition of Glacier Ice
Glaciers are one of the main water sources in arid South-Eastern Altai. To estimate their role in alimentation of mountain rivers, we required data on the isotopic composition of glacier ice. Sampling of glacier ice was carried out in three massifs: Tavan  The Tavan-Bogd massif is a large centre of modern glaciation. The northern slope of the massif is located in Russia. The south-eastern part belongs to Mongolia and the southwestern part belongs to China. The northern slope of the Tavan-Bogd massif is characterized by smaller areas of glaciation. In 2015 there were 16 glaciers on its territory with a total area of 23.46 km 2 . The average height of the firn border is 3335 m. One of the features of this massif is an increase in precipitation from east to west. The annual amount of precipitation in the highlands of the massif extends from 364 mm in the eastern part of the massif to 880 mm in the western part. The average height of the firn border decreases from 3415 m in the east to 3150 m in the west of the massif [16].
The Mongun-Taiga massif is located in the far east of Altai. The distribution of heat and moisture varies due to the significant height of the massif and the complexity of its The Tavan-Bogd massif is a large centre of modern glaciation. The northern slope of the massif is located in Russia. The south-eastern part belongs to Mongolia and the south-western part belongs to China. The northern slope of the Tavan-Bogd massif is characterized by smaller areas of glaciation. In 2015 there were 16 glaciers on its territory with a total area of 23.46 km 2 . The average height of the firn border is 3335 m. One of the features of this massif is an increase in precipitation from east to west. The annual amount of precipitation in the highlands of the massif extends from 364 mm in the eastern part of the massif to 880 mm in the western part. The average height of the firn border decreases from 3415 m in the east to 3150 m in the west of the massif [16].
The Mongun-Taiga massif is located in the far east of Altai. The distribution of heat and moisture varies due to the significant height of the massif and the complexity of its relief. Most precipitation falls on the windward slopes and in the high-altitude zone. In 2010, there were 32 glaciers in the massif with a total area of 20.3 km 2 (The Mongun-Taiga Mountain Massif, 2012). The average height of the firn border is 3390 m. The amount of precipitation on it was estimated at 269 mm [17].
Tsambagarav massif is located in the north-west of Mongolia. In 2015, there were 67 glaciers on the ridge with a total area of 68.4 km 2 ; the weighted average height of the firn border is 3748 m. There are large valley and cirque-valley glaciers on the north-eastern slope, the tongues of the largest of them descend to a height of 3000 m. The firn border on the glaciers of the north-eastern slope of Tsambagarav is located at an altitude 3500-3750 m.
Most of the samples were drawn in the ablation zones of valley glaciers. In the Mongun-Taiga massif, samples of ice from corrie glaciers and from the glacier dome on the top (3976 m.e.s.l.) were collected. Isotopic content of ice from different types of glaciers has noticeable differences. Box plots for δ 18 Table 1.
Mountain Massif, 2012). The average height of the firn border is 3390 m. The amount of precipitation on it was estimated at 269 mm [17].
Tsambagarav massif is located in the north-west of Mongolia. In 2015, there were 67 glaciers on the ridge with a total area of 68.4 km 2 ; the weighted average height of the firn border is 3748 m. There are large valley and cirque-valley glaciers on the north-eastern slope, the tongues of the largest of them descend to a height of 3000 m. The firn border on the glaciers of the north-eastern slope of Tsambagarav is located at an altitude 3500-3750 m.
Most of the samples were drawn in the ablation zones of valley glaciers. In the Mongun-Taiga massif, samples of ice from corrie glaciers and from the glacier dome on the top (3976 m.e.s.l.) were collected. Isotopic content of ice from different types of glaciers has noticeable differences. Box plots for δ 18      A maximum range of δ 18 O is observed in the ice samples from the ablation zones (ice tongue) of the valley glaciers. For corrie glaciers, mean δ 18 O values are lower, and the range is minimal. The samples of the glacier dome differ from other groups of samples significantly: the mean value is higher and the range of δ 18 O is not as large as the ranges from valley glacier samples. The reason of these differences will be discussed in part 4.
In Figure 3, box plots for δ 18 O of valley glaciers ice in three research areas are presented. Mean values are also presented in Table 2. The isotopic composition of glacier ice in the Tsambagarav massif was analysed more in a previous manuscript [18]. In Figure 3 Table 2, an extra 27 samples from Tsambagarav, which were collected in 2019, were added to the statistics.

and in
In Figure 3, box plots for δ 18 O of valley glaciers ice in three research areas are presented. Mean values are also presented in Table 2. The isotopic composition of glacier ice in the Tsambagarav massif was analysed more in a previous manuscript [18]. In Figure 3 and in Table 2, an extra 27 samples from Tsambagarav, which were collected in 2019, were added to the statistics.  These isotopic characteristics of ice can be used in meltwater and river runoff research. In the next chapter, the isotopic characteristics of glacier-originated water will be presented.

Isotopic Composition of Glacier Meltwater
Most of the meltwater samples were taken in the Tavan-Bogd massif. Sampling was carried out in 2015 [19], 2018 and 2021 in the middle of ablation season. Meltwater samples were taken from the edges of most glaciers on the northern slope (Russian part of this massif). Snow and firn samples were also collected from the glacier surface for further isotopic separation. Figure 4 illustrates the scheme of the glaciers on the northern slope of Tavan-Bogd.  The range of isotopic composition decreases and the mean values of δ 18 O increase from Tavan-Bogd to Tsambagarav.
These isotopic characteristics of ice can be used in meltwater and river runoff research. In the next chapter, the isotopic characteristics of glacier-originated water will be presented.

Isotopic Composition of Glacier Meltwater
Most of the meltwater samples were taken in the Tavan-Bogd massif. Sampling was carried out in 2015 [19], 2018 and 2021 in the middle of ablation season. Meltwater samples were taken from the edges of most glaciers on the northern slope (Russian part of this massif). Snow and firn samples were also collected from the glacier surface for further isotopic separation. Figure 4 illustrates the scheme of the glaciers on the northern slope of Tavan-Bogd. The isotopic composition of the meltwater from the main glaciers is displayed in Table 3. Isotopic composition of the meltwater from the two largest glaciers (No. 11 and No. 14 at Figure 4  These isotopic differences are related with morphometric types of glaciers and f tures of accumulation, and these relations will be discussed in part 4 of this study. Meltwater samples have also been collected in other research areas. The isotopic co position of the runoff in the Tsambagarav massif was analysed in previous work [18] the Mongun-Taiga massif meltwater δ 18 O varies from −14.8‰ to −17.9‰, but most of m water samples have δ 18 O ‰ close to −16‰. Diagrams of δ 18 O-δ 2 H relations for ice and meltwater samples (Table S1) for all search areas are presented below ( Figure 5, Table 4). The Tavan-Bogd massif meltwa samples, which were collected in 2018, are not included due to abnormal conditions t year (ice tongues of most of the glaciers were covered by snow even in the middle of su mer).  These isotopic differences are related with morphometric types of glaciers and features of accumulation, and these relations will be discussed in part 4 of this study.
Meltwater samples have also been collected in other research areas. The isotopic composition of the runoff in the Tsambagarav massif was analysed in previous work [18]. In the Mongun-Taiga massif meltwater δ 18 O varies from −14.8‰ to −17.9‰, but most of meltwater samples have δ 18 O ‰ close to −16‰.
Diagrams of δ 18 O-δ 2 H relations for ice and meltwater samples (Table S1) for all research areas are presented below ( Figure 5, Table 4). The Tavan-Bogd massif meltwater samples, which were collected in 2018, are not included due to abnormal conditions that year (ice tongues of most of the glaciers were covered by snow even in the middle of summer).   As Figure 5 illustrates, meltwater from the Tavan-Bogd (Figure 5a) is rather far from the mean ice value in the comparison to other two massifs.
The slope of the regression line between δ 18 O and δ 2 H for Tavan-Bogd samples is close to GMWL (δ 2 H = 8.14δ 18 O + 10.9). The slope of the regression for Mongun-Taiga and Tsambagarav samples is close to meteoric water line in Arid Central Asia (δ 2 H = 7.67 δ 18 O + 3.32) [20] due to the higher aridity in these massifs.

Differences in the Isotopic Composition from Different Types of Glaciers in the Mongun-Taiga Massif
On the top of the glacier dome, summer snow with high δ 18 O melts and refreezes on the surface of glacier. This is why δ 18 O of samples from the dome, which were taken in the middle of summer, were maximal and the range of δ 18 O values was minimal.
Valley glaciers have large accumulation zones, which include corries and cirques where winter snow concentrates. That is why δ 18 O in valley glacier ice has the largest range. It is possible that the mean δ 18 O value of valley glacier ice reflects the mean isotopic characteristics of the accumulated precipitation. Valley glaciers play the main role in the feeding of the glacier rivers, so this value was used in further isotopic separation.
The mean isotopic content of the corrie glacier ice is lower, but close to the valley glaciers. The range is minimal. Minimal δ 18 O can be explained by better conditions for the accumulation of precipitation during cold seasons. Minimal range can be explained by water refreezing predominance according to small size of this kind of glacier. However, small number of samples (20) does not allow to make any strong conclusions.

Differences in Isotopic Composition of Valley Glaciers Ice from Different Research Areas
The maximum range of δ 18 O values and the minimum average δ 18 O is observed for ice samples on the territory of the Tavan-Bogd massif. Increase in average values and decrease in the δ 18 O range in glacial ice samples in the rest of the massifs is associated with an increase in aridity and, as a result, with a decrease in annual precipitation and an increase of summer precipitation share in accumulation (Table 5). Table 5. Precipitation amount and share of the summer precipitation in research areas.

Massif Precipitation Amount on the Firn Border, mm
Share of the Summer Precipitation, % Tavan-Bogd  364-660  80  Mongun-Taiga  269  81  Tsambagarav  270  89 Since the Tavan-Bogd massif is located to the west from the others, more precipitation falls on its territory due to the predominance of westerlies in the Altai mountains. Despite the predominance of summer precipitation typical for the whole territory of arid highlands, winter precipitation with low isotopic composition plays a greater part in the formation of ice here than on the other massifs. This is reflected by larger range of values and a minimum mean value of δ 18 O. Summer precipitation falls mainly in liquid form and weakly participates in the formation of ice. Winter precipitation occurs more often than in other studied massifs, and, therefore, is more involved in the accumulation, which is reflected in its isotopic composition.
The ice of the Tsambagarav massif glaciers, which is located to the southeast of the other sites and exists in the most arid conditions, has higher isotopic composition with a minimum range of δ 18 O. This indicates higher a contribution of precipitation during the warm season, since winter precipitation is practically absent.
The Comparison with OIPC data (model data about isotopic composition of precipitation, which is based on GNIP stations interpolation) shows that the ice of the glaciers of these massifs is formed primarily by the precipitation of the transitional seasons (autumn-spring), but for the Tsambagarav massif, as was noted above, precipitation of the warm season plays a slightly larger role. For the Tavan-Bogd massif, this statement is less relevant [21,22]. Thus, the isotopic composition of the glacial ice of the arid highlands glaciers reflects well the conditions for the accumulation of precipitation, from which these glaciers form.

Meltwater Isotopic Composition
The glaciation of the northern slope of Tavan-Bogd is represented primarily by small corrie glaciers. The isotopic composition of the meltwater shows that in the middle of the ablation season, snow from the surface of glaciers plays an important role in the glacial rivers feeding in the northern part of the massif, in contrast to the rivers of the other massifs. It is explained by the predominance of morphological types of glaciers, which are favourable for snow accumulation. The ratio of the snow and ice components in the glacial runoff varies depending on the size and morphological type of the glacier. For most small corrie glaciers with favourable snow accumulation conditions, the proportion of melt snow water in the glacial runoff is greater than in the valley. This is observed both in the values of the isotopic composition of meltwater and in its variability over the years: Table 3 shows that meltwater δ 18  However, for the non-valley glaciers share of snow meltwater in total runoff is more significant. Thus, the ability of glaciers to accumulate large masses of spring snow on their surface is their important hydrological role on the northern slope of Tavan-Bogd.
During detailed glaciological observations that have been conducted on the territory of the massif since 1999, a tendency toward a rapid retreat of the valley glaciers was observed in this area. So, it can be assumed that the role of seasonal snow in river nutrition will only increase with time. On the other hand, with continued retreat in glaciation, spring-summer snowfall will melt much faster since it remains longer on the surface of glaciers.
Isotopic composition of the glacier river water in the Mongun-Taiga and Tsambagarav massifs shows less impact of seasonal snow melting on total runoff, because most of the water samples in these massifs are closer to the mean glacier δ 18 O ( Figure 5).

Conclusions
Before isotopic separation, it is necessary to determine what exactly is meant by "glacial runoff". There are two main outlooks on its definition. According to the first, glacial runoff is all water formed as a result of melting on the surface of a glacier [23].
According to the second point of view glacial runoff is the runoff formed due to the melting of long-term reserves of firn and ice. This approach reflects the main hydrological role of glaciers, which consists of the accumulation of precipitation and their temporary exclusion from the water cycle with subsequent redistribution by years.
In this study, glacial runoff is understood only as melting of long-term reserves of firn and ice, which corresponds to the second approach. The choice of this approach is related to the fact that for studied glaciers it is important to understand, how ice and firn reserves influence the formation of runoff in high mountain regions. Seasonal precipitation will feed the mountain rivers even in the absence of glaciers.
Isotopic separation in the Tavan-Bogd massif was made for all streams near the glaciers edge. In the other massifs separation was carried out for main streams, which start in general at the edges of the valley glaciers, because corrie glaciers are not so widespread and do not produce large streams in these regions. As a clear glacier meltwater isotopic composition, the mean δ 18 O of valley glaciers ice was taken. As a seasonal snow meltwater isotopic composition, mean δ 18 O of snow from snow pits was used. The results are presented in Table 6. The results of isotopic separation prove suggestions from Section 4.3 of this article: seasonal snow from the surface of glaciers has more of a significant influence on glacier runoff in the Tavan-Bogd massif. The main reasons for this are the better conditions for accumulation of summer snow, according to the widespread occurrence of small corrie glaciers.
It can be concluded that the further degradation of glaciation will have a greater impact on the formation of glacial runoff in Mongun-Taiga and Tsambagarav, than in the Tavan-Bogd massif, where the contribution of seasonal snow is more significant.
In addition, the isotopic characteristics of different objects presented in this manuscript can be used in further hydrological and glaciological research.