Received November 14, 2022; revised May 24, 2023; accepted June 27, 2023The characteristic feature of the Putorana Plateau is that the glaciological objects here are mostly represented by small glaciers and perennial snow patches. Their regime and morphology have so much common features that separation of these two categories of nival glacial formations from one another is extremely difficult problem. The distinctions between results of earlier studies carried out in the 1970s (the USSR Glacier Inventory estimated local resources at 22 glaciers with a total area of 2.5 km2) and at the beginning of the current century (in 2005, V.A. Sarana identified 61 glaciers with a total area of 7 km2) are too large to make any reliable conclusions about the current trends of the Putoran glaciation. In such conditions, the materials of rare field monitoring work performed on individual nival-glacial bodies become very valuable. Three expedition seasons of 2002–2004 included thorough geodetic and mass-balance measurements on 3 reference objects on the northern ledge of Lama Mts. – Prives (No. 30) Glacier, Marlborough (No. 31) Glacier and Strudoms snow patch. Similar field monitoring was repeated in summer’2019. The change in their configuration according to digital photography data, including the results of the analysis of satellite images, make possible to reveal that interannual fluctuations in the area of each glacier occur due to different vector displacements of its different sections along their entire perimeters. In 2019, the somewhat unexpected good budget state of all 3 monitored objects was detected. It manifested itself in growth of their areas and increased fraction of firn basin as well as in larger water equivalent of the firn residue as this followed from results of snow surveys. This fact contradicts the previous conclusion about the steady trend of deglaciation on the Putorana Plateau. The favourable condition for small glaciers occurred here in 2019 was mainly caused by reduced ablation owing to the weakened insolation that resulted from tremendous forest fires in Siberia, remarkable by their abnormal intensity in this year.

The main approaches to mountain glacier simulation were reviewed, and Oerlemans minimal model was chosen as a parameterization core of mountain glaciation in the Earth System models. The proposed model is based on a one-dimensional equation of glacier mass balance. The mass balance components are calculated with the use of a specially developed model of orographic precipitation, an algorithm for correction of incoming solar radiation onto an inclined ice surface, and schemes of calculation of turbulent heat-moisture exchange based on the Monin–Obukhov theory. The model was implemented for the Djankuat glacier (the Central Caucasus), for which a long measurement series is available. The model gives a good reproduction of the dynamics of glacier length over the period 1985–2020 based on measured mass balance values: –13 m/year, which is in practically perfect agreement with the field data. This means that the Oerlemans model can be used in Earth system models. The results of simulation based on the calculated mass balance showed a significant positive trend in ablation at a slight change in accumulation, which is also in agreement with the reality. However, in this case, the values of the annual thawing depth and the glacier contraction are twice as large and those observed in reality. The further development of the model (the inclusion of a snow cover block, the incorporation of debris mantle and mountain–valley circulation) will eliminate these shortcomings.

Mass balance bn is the most important indicator of a glacier evolution. However, after the decay of the USSR direct measurements of bn performed in 1985–1989 in the Inner Tien Shan, including the Sary-Tor Glacier in Ak-Shiyrak Massif, had been desisted. As a result the available series of the data were limited 1989. Measurements in this area were renewed only in 2015. This paper is devoted to restoring the continuity of the mass-balance series over the period of the gap in measurements and extending this series down to 1929, i.e. to the beginning of regular meteorological observations on the reference HMS Tien Shan (3660 m a.m.s.l.). Accumulation was reconstructed using a linear relationship of bn with the air temperature and precipitation sum. Reconstruction of ablation was based on its cubic relationship with the temperature (modified Krenke– Khodakov formula) or on two-parameter linear approximation using the air temperature and wind velocity. Thereby, the decade of direct instrumental measurements (1984/85–1988/89 and 2014/15–2018/19) resulted in deriving and analyzing continuous 90-year-long series of annual values of bn and its constituents, analytical type of referent glacio-meteorological equations being assumed unchanged in time. Reconstruction for the Sary-Tor Glacier reveals a dominant trend towards the mass loss with rare and short-time episodes of retarding the negative tendencies. The comparison made with the long series of mass balance of other glaciers in Asia indicates a certain degree of synchronicity, which is slightly disturbed in recent years: the degradation of Sary Tor Glacier tends to progress more intensively. Conclusions about its evolution are particularly relevant in connection with the assumption about the impact of the Kyrgyz-Canadian gold mining company «Kumtor Gold Company» on local ecosystems against the background of its interest in expanding the mining zone to the bowels of the Earth under the tongue of this glaciological object.

Mass balance bn is the most important indicator of a glacier evolution. However, after the decay of the USSR direct measurements of bn performed in 1985-1989 in the Inner Tien Shan, including the Sary-Tor Glacier in Ak-Shiyrak Massif, had been desisted. As a result the available series of the data were limited 1989. Measurements in this area were renewed only in 2015. This paper is devoted to restoring the continuity of the mass-balance series over the period of the gap in measurements and extending this series down to 1929, i.e. to the beginning of regular meteorological observations on the reference HMS Tien Shan (3660 m a.m.s.l.). Accumulation was reconstructed using a linear relationship of bn with the air temperature and precipitation sum. Reconstruction of ablation was based on its cubic relationship with the temperature (modified Krenke-Khodakov formula) or on two-parameter linear approximation using the air temperature and wind velocity. Thereby, the decade of direct instrumental measurements (1984/85-1988/89 and 2014/15-2018/19) resulted in deriving and analyzing continuous 90-year-long series of annual values of bn and its constituents, analytical type of referent glacio-meteorological equations being assumed unchanged in time. Reconstruction for the Sary-Tor Glacier reveals a dominant trend towards the mass loss with rare and short-time episodes of retarding the negative tendencies. The comparison made with the long series of mass balance of other glaciers in Asia indicates a certain degree of synchronicity, which is slightly disturbed in recent years: the degradation of Sary Tor Glacier tends to progress more intensively. Conclusions about its evolution are particularly relevant in connection with the assumption about the impact of the Kyrgyz-Canadian gold mining company "Kumtor Gold Company" on local ecosystems against the background of its interest in expanding the mining zone to the bowels of the Earth under the tongue of this glaciological object.

Releases of avalanches are an important factor of snow accumulation for many mountain glaciers. However, the quantitative estimate of the share of avalanche snow in the incoming part of the mass balance has not yet been standardized due to the lack of a universal methodological scheme for the calculation and collection of actual material. In rare cases, e.g. for the Djankuat Glacier in the Caucasus, this problem is partly solved. Here, the conceptual basis of the required calculations is created, and the input information has been supplied for many years by the route snow-measuring profling performed in the spring close to the date of the maximum of seasonal snow reserves and accompanied by direct feld mapping. Annually, avalanche deposits attributed to additional snow income from outside the glacier are recorded in all high-altitude morphological zones, with the exception of the lowest (< 2850 m) hypsometric belt, although in the half-century history of monitoring there were cases when snow avalanches reached it. Te volume of avalanche deposits was calculated for 15 years afer 1991/92, and each time it correlated with the gross snow accumulation of the corresponding balance year. Te contribution of snow avalanches to the income of substances on the glacier varies from 1.8 to 10.0% and averages 4.7%. Paradoxical (but only at frst glance) conclusion based on the analysis of data on volumes of avalanche and total accumulation is the feedback of these indicators (r = -0.58). It should be noted that the share of avalanche contribution to alimentation of the glacier has a more convincing feedback with accumulation and background snow content of winter (r = -0.72). Probably, the revealed tendency of increasing role of the additional avalanche alimentation in low-snow years (and vice versa) is not limited by the only reference object, and it is true for any glacier prone to the avalanche influence.

Current glaciation of the Central Caucasus is in the regressive stage of evolution. Observed long rising of the air temperature in the ablation seasons (almost 1 °C since 1968 on the Djankuat Glacier) together with insignificant growth of winter precipitation caused a long period of the glacier degradation: since 1871 the mass balance of the glacier, with rare exceptions, remained negative. During this time, the glacier has lost about 60 m w.e. mainly by areas of the smallest ice thickness, i.e. the steep slopes of the rocky framing of the fi n basin which become exposing from ice due to conditions of extreme instability. In parallel with the processes of freezing and thawing in cracks, this leads to a signifi ant increase in the frequency of rockfalls. In 2001–2003, huge rockfalls occurred on the glacier from the slopes, which had recently undergone deglaciation. The last one covered 4% of the glacier surface with about 70 000 m3 of clastic products. The marine cover differs significantly from the open ice surface in its thermal characteristics, which affect the structure of the thermal balance and the absolute values of ablation. The temperature inside the moraine is positive throughout the ablation season, even at a depth of 60 cm, but the absolute values of daily maximum melting are five times lower than those on open ice. The influence of the moraine material layer on the ablation is determined not only by its thickness, but also by the structure. The fi ld experiment conducted in 2011 demonstrated that melting under a layer of coarse-grained material (particle diameter 16–17 cm) is four times faster than under a layer of fi e-grained (7–8 cm) material equal in thickness. The study of the influence of the surface moraine and the account of its further growth in thickness and area made it possible to predict changes in the surface of the glacier by 2025. The climatic forecast was made using the HadCM3 model (A2 climate change scenario). Data on air temperature and precipitation were used to construct a predictive fi ld of the mass balance of the Djankuat Glacier for each year until 2025. Th s became possible by the use of a quantitative parameterized relationship between the mass balance of the entire glacier and each of its points, identifi d in the study of the total mass transfer of the glacier, as well as parameterized linear regression equations. The forecast of hypsometry of the Djankuat Glacier made it possible to determine that by 2025 the depression of the tongue areas covered by the moraine will be 7–15 m (or 42–45%) smaller than areas of the open ice. Th s means that at the present stage of the glaciation development the surface moraine plays a role in the evolution of the glacier, commensurable with the climate factor.

In this paper we present and discuss results of radio-echo sounding and modelling of ice thickness of the Djankuat Glacier. This glacier is the reference one for the Central Caucasus and was earlier studied comprehensively, but direct ice thickness measurements survey were not made so far. In 2012–2013, more than 20 km of ice-thickness measurement profiles were obtained using monopulse ground radar VIRL-6 with its central frequency of 20 MHz. Standard error of the ice-thickness measurement was 2.5%. Detailed maps of the ice-thickness and the bedrock topography based on these direct measurements were constructed for the first time. Its average ice-thickness is 31 m with the maximum of 105 m. Djankuat Glacier ice volume is 0.077±0.002 km3 when Djantugan plateau is not taken into account. Ice thickness map was supplemented by results of the ice thickness modelling using the GlabTop model. It was shown that the model reproduces the ice thickness distribution correctly, and a special approach was developed to correct model parameters using ice thickness measurements. Further on, we plan to use corrected GlabTop model for estimation glaciations of the Caucasus as well as to carry out ground radar sounding of the Djantugan plateau.

In order to evaluate the contribution of snow avalanches, coming from outside the glacier limits, into the seasonal accumulation, an original algorithm is proposed for glaciers with ongoing terrestrial mass balance programme. It was tested on the Djankuat Glacier, Caucasus. A case study was realized for 1998/99 balance year when accumulation value was the closest to its long-term mean over the 46-year-long monitoring period. Snow avalanche feeding turned out to be only 2.8 per cent of total accumulation that is about an order of magnitude smaller than assumed hitherto.

The recession of mountain glaciers around the world has been linked to anthropogenic climate change and small glaciers (e.g. <2km2) are thought to be particularly vulnerable, with reports of their disappearance from several regions. However, the response of small glaciers to climate change can be modulated by non-climatic factors such as topography and debris cover and there remain a number of regions where their recent change has evaded scrutiny. This paper presents results of the first multi-year remote sensing survey of glaciers in the Kodar Mountains, the only glaciers in SE Siberia, which we compare to previous glacier inventories from this continental setting that reported total glacier areas of 18.8km2 in ca. 1963 (12.6km2 of exposed ice) and 15.5km2 in 1974 (12km2 of exposed ice). Mapping their debris-covered termini is difficult but delineation of debris-free ice on Landsat imagery reveals 34 glaciers with a total area of 11.72±0.72km2 in 1995, followed by a reduction to 9.53±0.29km2 in 2001 and 7.01±0.23km2 in 2010. This represents a ~44% decrease in exposed glacier ice between ca. 1963 and 2010, but with 40% lost since 1995 and with individual glaciers losing as much as 93% of their exposed ice. Thus, although continental glaciers are generally thought to be less sensitive than their maritime counterparts, a recent acceleration in shrinkage of exposed ice has taken place and we note its coincidence with a strong summer warming trend in the region initiated at the start of the 1980s. Whilst smaller and shorter glaciers have, proportionally, tended to shrink more rapidly, we find no statistically significant relationship between shrinkage and elevation characteristics, aspect or solar radiation. This is probably due to the small sample size, limited elevation range, and topographic setting of the glaciers in deep valleys-heads. Furthermore, many of the glaciers possess debris-covered termini and it is likely that the ablation of buried ice is lagging the shrinkage of exposed ice, such that a growth in the proportion of debris cover is occurring, as observed elsewhere. If recent trends continue, we hypothesise that glaciers could evolve into a type of rock glacier within the next few decades, introducing additional complexity in their response and delaying their potential demise.

The Buordakh Massif, in the Cherskiy Range of northeast Siberia, contains mountains over 3000 in and, despite its and climate, numerous glaciers. This paper presents a glacier inventory for the region and documents some 80 glaciers, which range in size from 0.1 to 10.4 km(2) (total glacierized area is ca. 70 km(2)). The inventory is based on mapping derived from Landsat 7 ETM+ satellite imagery from August 2001, augmented with data from field investigations obtained at that time. The glaciers in this region are of the 'firn-less,' cold, continental type, and their mass balance relies heavily on the formation of superimposed ice. The most recent glacier maximum extents have also been delineated, and these are believed to date from the Little Ice Age (ca. A.D. 1550-1850). Glacier areal extent has reduced by some 14.8 km(2) (ca. 17%) since this most. recent maximum. Of the 80 glaciers catalogued, 49 have undergone a measurable retreat from their most recent maximum extent.

Superficial moraines grew in size during the entire 32-year-long period of direct monitoring of water and ice balance of the Djankuat Glacier in the Caucasus. The total area of debris cover on the glacier increased from 0.104 km(2) (3% of the entire glacier surface) in 1968 to 0.266 km(2) (8% of the glacier) in 1996. Such rapid dynamics of moraine formation greatly influences the ablation rate and distorts fields of mass-balance components. Sub-debris thawing can be calculated by means of a model, which describes the role of debris cover for the thermal properties of a glacier. Its meltwater equivalent depends mainly on debris thickness. In 1983 and 1994 the debris cover was repeatedly mapped over the whole glacier portion that was covered with morainic material. Sub-moraine ablation increases (vs. pure ice surface) under the thin, less than ca. 7-8 cm, debris layer, whereas the thicker debris cover reduces the liquid runoff due to its shielding effect. Zones differing due to their hydrological effect are depicted on the glacier map and the degree of debris influence on ablation is estimated quantitatively. As a whole runoff from debris-covered parts of the Djankuat Glacier has diminished due to the dominant shielding effect. Variation of the terminus is also shown to be dependent on the evolution of superficial moraine.

The relationship between the thickness of a glacier stream and its surface slope is discussed. Morphometrical features of the glacier surface serve as a basis for estimating the maximum thickness the glacier can attain. A method for indirect calculation of ice thickness is hereby suggested. A detailed large-scale map permits determination of the total glacier volume. The Djankuat Glacier, a representative valley glacier in the Central Caucasus, is used as an example to illustrate the applicability of the proposed calculation scheme. Its probable thickness and volume are estimated on the grounds of photogrammetrical surveys, a topographic map and spatial pattern of ice flow velocity on the surface. Several values of directly measured ice thickness (by various methods such as thermal drilling, radar sounding, gravimetry etc .) verify the calculation results. As a result, ice thickness over the entire glacier area has been mapped. Methodical limitations and accuracy of indirect estimates are also discussed.