Glacier Balance Research Articles

The Role of the Snow Ratio in Mass Balance Change under a Warming Climate for the Dongkemadi Glacier, Tibetan Plateau

Abstract Understanding the effects of the snow ratio on glacier mass balance under variable climatic conditions is crucial for predicting how glaciers will respond to climate change, and for forecasting water supplies to surrounding lowland areas. Due to recent climate change, the historical annual snow ratio of the Dongkemadi (DKMD) Glacier showed a significant increasing trend (0.0538% a−1, p < 0.05), and an abrupt upward change in 1977 due to decreasing precipitation concentration. Snow ratios with fixed precipitation concentration and nonwarming climate scenarios were calculated to isolate the impact of the snow ratio on glacier mass balance. Under nonwarming conditions, the snow ratio showed little variability, ranging from 88.4% to 99.9%. Glacier modeling results comparing five snow ratio scenarios from 1961 to 2009 showed three main features as follows. (i) Glacier mass balance was low and more sensitive to a warming climate for lower snow ratio scenarios. (ii) The difference in mass balance between the scenarios fluctuated, but generally increased with time. Spatially, the ablation area change was larger (0.4 km2), and the equilibrium line altitude was higher (5.9 m) in scenarios with lower snow ratios. (iii) The change in net shortwave radiation was the main reason for changes in glacial melt, and the albedo played a key role in controlling the difference of glacier energy balance between snow ratio scenarios. Rain increment only accounted for about 20%–33% of meltwater increment. Overall, this study provides valuable information to evaluate how snow ratios impact the mass balance of glaciers with ongoing climate change.

Spatiotemporal heterogeneity and driving mechanisms of Himalayan glacier mass change in the early 21st century

Glaciers in the Himalayan mountain ranges contain huge amounts of water and are essential for the water security of their downstream basins. Due to the spatiotemporal heterogeneity and uncertainty of the glacier mass change (GMC) in the Himalayas, more refined and deeper investigations on the GMC and its driving forces are needed to complement the current research. We present the mass balance of Himalayan glaciers during 2000 to 2013 using the latest global digital elevation models; this method focuses on establishing a high-coverage result of the GMC and resolving the problems of timestamps and penetration uncertainties. The GMC of the entire Himalayas is −0.32 ± 0.15 m w.e.yr − 1; most of the lost mass drains into the Indus and Ganges river basins. Glacier mass loss gradually decreases from east to west, with the East and Central Himalayas exhibiting the highest heterogeneity in the vertical and horizontal directions, respectively. We propose two multivariate regression models with 12 potential driving factors that explain 63% to 89% of the GMC across different regions of the Himalayas. Precipitation significantly influences changes in glacier mass throughout the Himalayas, whereas debris coverage and glacier size have insignificant effects on changing the glacier mass. Our analysis reveals that seasonal climate-driven changes in the glacier mass of the Himalayas may threaten the downstream water supply.

Understanding monsoon controls on the energy and mass balance of glaciers in the Central and Eastern Himalaya

Abstract. The Indian and East Asian summer monsoons shape the melt and accumulation patterns of glaciers in High Mountain Asia in complex ways due to the interaction of persistent cloud cover, large temperature ranges, high atmospheric water content and high precipitation rates. Glacier energy- and mass-balance modelling using in situ measurements offers insights into the ways in which surface processes are shaped by climatic regimes. In this study, we use a full energy- and mass-balance model and seven on-glacier automatic weather station datasets from different parts of the Central and Eastern Himalaya to investigate how monsoon conditions influence the glacier surface energy and mass balance. In particular, we look at how debris-covered and debris-free glaciers respond differently to monsoonal conditions. The radiation budget primarily controls the melt of clean-ice glaciers, but turbulent fluxes play an important role in modulating the melt energy on debris-covered glaciers. The sensible heat flux decreases during core monsoon, but the latent heat flux cools the surface due to evaporation of liquid water. This interplay of radiative and turbulent fluxes causes debris-covered glacier melt rates to stay almost constant through the different phases of the monsoon. Ice melt under thin debris, on the other hand, is amplified by both the dark surface and the turbulent fluxes, which intensify melt during monsoon through surface heating and condensation. Pre-monsoon snow cover can considerably delay melt onset and have a strong impact on the seasonal mass balance. Intermittent monsoon snow cover lowers the melt rates at high elevation. This work is fundamental to the understanding of the present and future Himalayan cryosphere and water budget, while informing and motivating further glacier- and catchment-scale research using process-based models.

Large‐eddy simulations of the atmospheric boundary layer over an Alpine glacier: Impact of synoptic flow direction and governing processes

AbstractThe mass balance of mountain glaciers is of interest for several applications (e.g., local hydrology or climate projections), and turbulent fluxes can be an important contributor to glacier surface mass balance during strong melting events. The underlying complex terrain leads to spatial heterogeneity and non‐stationarity of turbulent fluxes. Owing to the contribution of thermally induced flows and gravity waves, exchange mechanisms are fully three‐dimensional, instead of only vertical. Additionally, glaciers have their own distinct microclimate, governed by a down‐glacier katabatic wind, which protects the glacier ice and interacts with the surrounding flows on multiple scales. In this study, we perform large‐eddy simulations with the Weather Research and Forecasting model at a horizontal grid spacing of 48 m to gain insight into the boundary‐layer processes over an Alpine valley glacier, the Hintereisferner. We choose two case studies from the Hintereisferner experiment measurement campaign with different synoptic wind directions (southwest and northwest). Model evaluation with an array of eddy‐covariance stations on the glacier tongue and surroundings reveals that the Weather Research and Forecasting model is able to simulate the general glacier boundary‐layer structure. Under a southwesterly airflow, the down‐glacier wind is supported by the synoptic wind parallel to the glacier axis, a stable boundary layer is present over the ice surface, and local processes govern the turbulence kinetic energy production. Under northwesterly airflow, a cross‐glacier valley flow and a breaking gravity wave lead to strong turbulent mixing and to the subsequent erosion of the glacier boundary layer. Stationarity analysis of the sensible heat flux suggests non‐stationary behaviour for both case study days, whereas non‐stationarity is highest on the northwesterly day during the gravity‐wave event. These results suggest that the synoptic wind direction has, in addition to upstream topography and the atmospheric stability, a strong impact on whether a local glacier boundary layer can form or not, influencing whether a glacier is able to maintain its own microclimate.

Influence of Topographic Shading on the Mass Balance of the High Mountain Asia Glaciers

Most studies attribute the glacier mass balance within High Mountain Asia (HMA) to climate change, ignoring the influence of its complex terrain. Knowledge of the influence of this complex terrain is crucial for understanding the spatial variability in its mass balance. However, there is a lack of any systematic assessment of this influence across HMA. Therefore, in this study, we used the glacier outlines and raster data (SRTM DEM, slope and aspect) to calculate the topographic shading of all 97,965 glaciers within HMA during the ablation period, which is regarded as a major index of the influence of complex terrain on the mass balance. The results showed that 27.19% of HMA glacier area was subjected to topographic shading, and regional differences were significant with respect to both their altitudinal and spatial distributions. The topographic shading contributed to the protection of the smallest glaciers from solar illumination. Furthermore, we found a significant correlation between the topographic shading and mass balance in these small north-facing glaciers. However, these small glaciers were most prevalent in the north-facing orientation, especially in West Kunlun, East Kunlun, Inner Tibet Plateau and Qilian Shan, where shading was found to increase with decreases in the glacier area. This indicates that complex terrain can affect the spatial distribution of the mass balance by altering the solar illumination pattern.

MODIS-observed spatiotemporal changes in surface albedo of Karakoram glaciers during 2000-2018

The role of albedo is very important in modulating the surface energy balance of glaciers. The main objective of this study is to assess the spatiotemporal variability in surface albedo of the Karakoram glaciers in Pakistan during the summer seasons (June, July and August) for the period from 2000-2018. We used Moderate Resolution Imaging Spectroradiometer (MODIS) data to estimate the amount of glacier surface albedo. We combined the MODIS Terra- and Aqua-derived albedo products to reduce the amount of cloud influence and to improve the estimation of glacier surface albedo. Our results indicate that the average annual decrease in albedo is ~0.041% during the summer. The decrease in albedo was relatively high during recent years, with an annual rate of decrease of ~0.45%. The decreasing trend in albedo is towards the north-western part of the Karakoram mountain range. Climate change is the potential cause of albedo variations in the study area. Albedo has a strong negative correlation with temperature (r = -0.811) and a strong positive correlation with precipitation (r = 0.809). The present study concludes that trend in decreasing albedo is higher during the recent years than the last decade and climate change is playing a vital role in it.

Processing of VENµS Images of High Mountains: A Case Study for Cryospheric and Hydro-Climatic Applications in the Everest Region (Nepal)

In the Central Himalayas, glaciers and snowmelt play an important hydrological role, as they ensure the availability of surface water outside the monsoon period. To compensate for the lack of field measurements in glaciology and hydrology, high temporal and spatial resolution optical remotely sensed data are necessary. The French–Israeli VENµS Earth observation mission has been able to complement field measurements since 2017. The aim of this paper is to evaluate the performance of different reflectance products over the Everest region for constraining the energy balance of glaciers and for cloud and snow cover mapping applied to hydrology. Firstly, the results indicate that a complete radiometric correction of slope effects such as the Gamma one (direct and diffuse illumination) provides better temporal and statistical metrics (R2 = 0.73 and RMSE = 0.11) versus ground albedo datasets than a single cosine correction, even processed under a fine-resolution digital elevation model (DEM). Secondly, a mixed spectral-textural approach on the VENµS images strongly improves the cloud mapping by 15% compared with a spectral mask thresholding process. These findings will improve the accuracy of snow cover mapping over the watershed areas downstream of the Everest region.

North Atlantic Cooling is Slowing Down Mass Loss of Icelandic Glaciers

AbstractIcelandic glaciers have been losing mass since the Little Ice Age in the mid‐to‐late 1800s, with higher mass loss rates in the early 21st century, followed by a slowdown since 2011. As of yet, it remains unclear whether this mass loss slowdown will persist in the future. By reconstructing the contemporary (1958–2019) surface mass balance of Icelandic glaciers, we show that the post‐2011 mass loss slowdown coincides with the development of the Blue Blob, an area of regional cooling in the North Atlantic Ocean to the south of Greenland. This regional cooling signal mitigates atmospheric warming in Iceland since 2011, in turn decreasing glacier mass loss through reduced meltwater runoff. In a future high‐end warming scenario, North Atlantic cooling is projected to mitigate mass loss of Icelandic glaciers until the mid‐2050s. High mass loss rates resume thereafter as the regional cooling signal weakens.

Density matters: ice compressibility and glacier mass estimation

Ice flow models typically assume that ice is incompressible, a reasonable assumption because ice density changes are indeed small and have a correspondingly small effect on the overall mass balance of glaciers and ice sheets. Given the immense volume of the ice sheets, however, even relatively small changes may influence global mean sea level to a degree that severely impacts humanity (Hauer and others, 2020). Here, we quantify the role of gravitational compression and thermal contraction in estimating ice sheet mass.

On the importance of subsurface heat flux for estimating the mass balance of alpine glaciers

Glaciers, as massive freshwater reservoirs, support the planet's living systems and have an impact on our daily lives, even for communities living far away. Ongoing and future climate change is predicted to have strong impacts on the mass balance of alpine glacier around the world. To understand the relationship between climate and glacier dynamics, a range of mass balance models are currently used. Most of these models however, ignore subsurface heat fluxes as a component of glacier mass balance. Here, we set out to investigate the importance of subsurface heat flux for the mass balance of an alpine glacier using a surface energy mass balance model (SEM) coupled with a multilayer subsurface heat conduction model (MSHCM) that resolves the subsurface glacier temperature. As a case study, we investigate the Urumqi Glacier No.1 in the Tianshan Mountains (NW China), which has a long and continuous time series of surface and subsurface glacier temperature measurements. We evaluate the results of both glacier temperature models (SEM and MSHCM) using these in situ observations and investigate the sensitivity of mass balance to five meteorological factors: air temperature, precipitation, incoming shortwave radiation, relative humidity, and wind speed. The mass balance of the glacier was simulated first by including the influence of subsurface heat flux, and second, the subsurface heat flux was neglected. Observed and simulated mass balance and the englacial temperature were found to be reasonably close in both cases. Furthermore, the mass balance was simulated with a zero surface temperature assumption, which resulted in a 6% overestimation of the summer ablation. We concluded that the mass balance of Urumqi Glacier No.1 was most sensitive to variations in temperature, followed by precipitation. Furthermore, our results show that subsurface heat flux in the ablation area can generally be neglected in estimating the mass balance of alpine glaciers during ablation season.

Reconstruction of the Historical (1750–2020) Mass Balance of Bordu, Kara-Batkak and Sary-Tor Glaciers in the Inner Tien Shan, Kyrgyzstan

The mean specific mass balance of a glacier represents the direct link between a glacier and the local climate. Hence, it is intensively monitored throughout the world. In the Kyrgyz Tien Shan, glaciers are of crucial importance with regard to water supply for the surrounding areas. It is therefore essential to know how these glaciers behave due to climate change and how they will evolve in the future. In the Soviet era, multiple glaciological monitoring programs were initiated but these were abandoned in the nineties. Recently, they have been re-established on several glaciers. In this study, a reconstruction of the mean specific mass balance of Bordu, Kara-Batkak and Sary-Tor glaciers is obtained using a surface energy mass balance model. The model is driven by temperature and precipitation data acquired by combining multiple datasets from meteorological stations in the vicinity of the glaciers and tree rings in the Kyrgyz Tien Shan between 1750 and 2020. Multi-annual mass balance measurements integrated over elevation bands of 100 m between 2013 and 2020 are used for calibration. A comparison with WGMS data for the second half of the 20th century is performed for Kara-Batkak glacier. The cumulative mass balances are also compared with geodetic mass balances reconstructed for different time periods. Generally, we find a close agreement, indicating a high confidence in the created mass balance series. The last 20 years show a negative mean specific mass balance except for 2008–2009 when a slightly positive mass balance was found. This indicates that the glaciers are currently in imbalance with the present climatic conditions in the area. For the reconstruction back to 1750, this study specifically highlights that it is essential to adapt the glacier geometry since the end of the Little Ice Age in order to not over- or underestimate the mean specific mass balance. The datasets created can be used to get a better insight into how climate change affects glaciers in the Inner Tien Shan and to model the future evolution of these glaciers as well as other glaciers in the region.

Retrieval of Firn Thickness by Means of Polarisation Phase Differences in L-Band SAR Data

The knowledge of glacier zones’ extent and their temporal variations is fundamental for the retrieval of surface mass balance of glaciers and ice sheets. In this context, a key parameter is the firn line (FL), the lower boundary of the percolation zone, whose location is an indicator of time-integrated mass balance changes. Several approaches have been developed in the last decades to map the FL by means of Synthetic Aperture Radar (SAR) imagery, mainly exploiting backscatter intensities and their seasonal variation. In this paper, an alternative approach is proposed, based on co-polarisation phase differences (CPDs). In particular, CPDs are interpreted as the result of propagation through anisotropic firn layers and are, therefore, proposed as an indicator of the presence of firn. A model is employed to demonstrate the link between CPDs and firn depth, indicating the potential of polarimetric SAR to improve firn characterization beyond spatial extent and FL detection. The proposed approach is demonstrated on L-band airborne data, acquired on 21 May 2015 by the F-SAR sensor of DLR in West Greenland during the ARCTIC15 campaign, and validated with in-situ information available from other studies.

SBAS-InSAR-Based Analysis of Surface Deformation in the Eastern Tianshan Mountains, China

Due to the unique geographical characteristics of cold alpine and high-altitude regions, glaciers, permafrost, ground ice, rock glaciers, and other periglacial geomorphology have developed with fragile habitats, and these areas are often the birthplaces of many river basins and natural hazards. With global warming and the extensive cryogenesis and physical weathering, the thermal state of permafrost and the mass balance of glaciers have been changed, and thus it can be deduced that the surface deformation is of great concern. To obtain ground subsidence or uplift over a large area to understand local surface changes, the small baseline subset interferometric synthetic aperture radar (SBAS-InSAR) technique was applied to process 89-scene of Sentinel-1A images ranging from December 25, 2017 to January 2, 2021 to obtain surface deformation for these 3 years for the eastern Tianshan Mountains, China. The surface deformation characteristics of the area were analyzed to provide a basic dataset for environmental protection policies and mitigation or reduction of natural hazards in this region, and to verify the applicability of SBAS-InSAR technology in alpine and high-altitude areas. The results show that the SBAS-InSAR technique processing with sentinel-1A dataset cannot be effectively used to acquire ground deformation in areas covered by trees, scrub/shrub, glaciers, snow, and ground ice, where the decohered phenomenon is serious. In other regions, SBAS-InSAR can effectively measure surface subsidence or uplift. Surface deformation is significant throughout the study area, with rates ranging from −70.7 to 50.8 mm/a and with an average rate of 1.1 mm/a. There are obvious regions of uplift in the northwest, northeast, and central sections of the study area, with uplift greater than 155.73 mm in 3 years, and three obvious regions of subsidence in the northeast and west sections of the study area, with subsidence of at least −125.20 mm in 3 years. The remaining areas of deformation are scattered, with smaller amounts of settlement and uplift and with an isolated and sporadic distribution. Areas with elevations of 3,150 to 4,275 m.a.s.l., slopes of 15°–50°, and southwest, west, and northwest aspects are geologic disaster-prone regions and should receive more attention and more field monitoring. The results of this study have important implications for local environmental protection and hazard prevention.

Quantifying Patterns of Supraglacial Debris Thickness and Their Glaciological Controls in High Mountain Asia

Mapping patterns of supraglacial debris thickness and understanding their controls are important for quantifying the energy balance and melt of debris-covered glaciers and building process understanding into predictive models. Here, we find empirical relationships between measured debris thickness and satellite-derived surface temperature in the form of a rational curve and a linear relationship consistently outperform two different exponential relationships, for five glaciers in High Mountain Asia (HMA). Across these five glaciers, we demonstrate the covariance of velocity and elevation, and of slope and aspect using principal component analysis, and we show that the former two variables provide stronger predictors of debris thickness distribution than the latter two. Although the relationship between debris thickness and slope/aspect varies between glaciers, thicker debris occurs at lower elevations, where ice flow is slower, in the majority of cases. We also find the first empirical evidence for a statistical correlation between curvature and debris thickness, with thicker debris on concave slopes in some settings and convex slopes in others. Finally, debris thickness and surface temperature data are collated for the five glaciers, and supplemented with data from one more, to produce an empirical relationship, which we apply to all glaciers across the entire HMA region. This rational curve: 1) for the six glaciers studied has a similar accuracy to but greater precision than that of an exponential relationship widely quoted in the literature; and 2) produces qualitatively similar debris thickness distributions to those that exist in the literature for three other glaciers. Despite the encouraging results, they should be treated with caution given our relationship is extrapolated using data from only six glaciers and validated only qualitatively. More (freely available) data on debris thickness distribution of HMA glaciers are required.

Shallow-Layers-Detection Ice Sounding Radar for Mapping of Polar Ice Sheets

The accumulation rate is a key parameter in computing the mass balance of glaciers and ice sheets to estimate sea level rise. A shallow-layers-detection ice sounding radar (SLDISR) is developed to measure the accumulation rate and shape of near-surface internal layers with high resolution. With a transmitting frequency from 500 to 2000 MHz, this frequency-modulated continuous wave (FMCW) radar provides a range resolution of about 16 cm in free space by using a Hanning window and a penetrating depth about 150 m under polar ice. The spectral analysis and coherent integration techniques are used to obtain a high processing gain and to improve the signal-to-noise ratio of the system. A phase-locked loop with wideband yttrium iron garnet (YIG) oscillator is applied to generate a sweeping chirp signal as an input source for the transmitter. A stable, low-frequency reference chirp signal is generated with a direct digital synthesizer (DDS) integrated in field-programmable gate array (FPGA). To reduce the high-speed requirement to the analog-to-digital converter (ADC), dechirp technology is adopted at the RF section of the receiver. The implementation of the digital unit is based on an FPGA chip. The designed radar has been successfully deployed in Antarctica during the 31st Chinese Antarctic Research Expedition (CHINARE 31) and CHINARE 33, mainly over the East Antarctic Ice Sheet (EAIS). The echograms indicate the effectiveness of the radar system on detecting clear internal reflecting horizons (IRHs) over ice sheets.

Detection of Spatial Shift in Flood Regime of the Kabul River Basin in Pakistan, Causes, Challenges, and Opportunities

Recent evidence of regional climate change impacts on hydrological cycle directed us to study the floods in a high elevated and rapidly urbanized river basin, the Kabul River basin (KRB), Pakistan, which is susceptible to frequent flooding. Therefore, we analyzed the changes in flood regime at various spatial and temporal scales and their possible causes, which is accomplished by using flood indicators, trend analysis, change point analysis, and hydrological modeling. The results showed that the northern and northwestern parts of the KRB were more exposed to flood hazard than the southern parts under long-term scenario (1961/64-2015). However, after the change points, the flood risk decreased in the northern and increased in the southern regions. This spatial shift increased the vulnerability of population to the flood hazard, because the majority of population resides in the southern region. The extreme precipitation has also increased, especially the maximum one-day rainfall and maximum five-day rainfall throughout the basin. Particularly, the major cause of the decrease in different flood indicators in the northern parts of the KRB is the corresponding decrease in the annual and monsoonal rainfall and corresponding positive mass balance of glaciers in the northern region after the occurrence of change point in flood regime. However, the major cause of the increase in flood hazard on the southern part of the KRB is associated with maximum five-day rainfall. A 68% variability of annual maximum flood for the Kabul River at Nowshera and an 84% variability of annual maximum flood for Bara River at Jhansi post are explained by maximum five-day rainfall. In addition, a considerable decrease in forests (–5.21%) and increase in the urban area (88.26%) from 1992–2015 also amplifies the risk of higher flood peaks. The results of hydrological modeling suggest that the six-hourly flood peak increased by 6.85% (1992–2010) and 4.81% (2010–2015) for the extreme flood of 2010 for the Kabul River at Nowshera. The flood peak per decade will increase by 8.6%, as compared to the flood peak under the land use scenario of 2010. Therefore, consideration of proper land use planning is crucial for sustainable flood management in the KRB.

Automated Detection of Marine Glacier Calving Fronts Using the 2-D Wavelet Transform Modulus Maxima Segmentation Method

Changes in the calving front position of marine-terminating glaciers strongly influence the mass balance of glaciers, ice caps, and ice sheets. At present, quantification of frontal position change primarily relies on time-consuming and subjective manual mapping techniques, limiting our ability to understand changes to glacier calving fronts. Here we describe a newly developed automated method of mapping glacier calving fronts in satellite imagery using observations from a representative sample of Greenland’s peripheral marine-terminating glaciers. Our method is adapted from the 2-D wavelet transform modulus maxima (WTMM) segmentation method, which has been used previously for image segmentation in biomedical and other applied science fields. The gradient-based method places edge detection lines along regions with the greatest intensity gradient in the image, such as the contrast between glacier ice and water or glacier ice and sea ice. The lines corresponding to the calving front are identified using thresholds for length, average gradient value, and orientation that minimize the misfit with respect to a manual validation data set. We demonstrate that the method is capable of mapping glacier calving fronts over a wide range of image conditions (light to intermediate cloud cover, dim or bright, mélange presence, etc.). With these time series, we are able to resolve subseasonal to multiyear temporal patterns as well as regional patterns in glacier frontal position change.

Estimation of mass and energy balance of glaciers using a distributed energy balance model over the Chandra river basin (Western Himalaya)

AbstractThe ongoing glacier shrinking in the Himalayan region causes a significant threat to freshwater sustainability and associated future runoff. However, data on the spatial climatic contribution of glacier retreat is scanty in this region. To investigate the spatially distributed glacier surface energy and mass fluxes, a two‐dimensional mass balance model was developed and applied to the selected glaciers of the Chandra basin, in the Upper Indus Basin, Western Himalaya. This model is driven by the remote sensing data and meteorological variables measured in the vicinity of the Chandra basin for six hydrological years (October 2013 to September 2019). The modelled variables were calibrated/validated with the in‐situ observation from the Himansh station in the Chandra basin. We have derived air temperature (Ta) spatially using the multivariate statistical approach, which indicates a relative error of 0.02–0.05°C with the observed data. Additionally, the relative error between the modelled and observed radiation fluxes was <10.0 W m−2. Our study revealed that the Chandra basin glaciers have been losing its mass with a mean annual mass balance of −0.59 ± 0.12 m w.e.a−1for the six hydrological years. Results illustrated that the mean surface melt rate of the selected glaciers ranged from −5.1 to −2.5 m w.e.a−1that lies between 4500 and 5000 m a.s.l. The study revealed that the net radiation(RN)contributes ~75% in total energy (FM) during the melt season while sensible heat(HS), latent heat(Hl), and ground heat(HG)fluxes shared 15%, 8%, and 2%, respectively. Sensitivity analysis of the energy balance components suggested that the mass balance is highly sensitive to albedo and surface radiations in the study area. Overall, the proposed model performed well for glacier‐wide energy and mass balance estimation and confirms the utility of remote sensing data, which may help in reducing data scarcity in the upper reaches of the Himalayan region.

РЕКОНСТРУКЦИЯ БАЛАНСА МАССЫ ЛЕДНИКА САРЫ-ТОР ПО МЕТЕОРОЛОГИЧЕСКИМ ДАННЫМ

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 of Glaciers in Bhaga Basin, Western Himalaya: A Geospatial and Temperature-Weighted AAR based Model Approach

Mass Balance of Glaciers in Bhaga Basin, Western Himalaya: A Geospatial and Temperature-Weighted AAR based Model Approach