Search

Investigating the characteristics and transformation of water-soluble carbonaceous matter in the cryosphere regions is important for understanding biogeochemical process in the earth system. Water-soluble carbonaceous matter is a heterogeneous mixture of organic compounds that is soluble in aquatic environments. Despite its importance, we still lack systematic understanding for dissolved organic carbon (DOC) in several aspects including exact chemical composition and physical interactions with microorganisms, glacier meltwater. This review presents the chemical composition and physical properties of glacier DOC deposited through anthropogenic emission, terrestrial, and biogenic sources. We present the molecular composition of DOC and its effect over snow albedo and associated radiative forcings. Results indicate that DOC in snow/ice is made up of aromatic protein-like species, fulvic acid-like materials, and humic acid-like materials. Light-absorbing impurities in surface snow and glacier ice cause considerable albedo reduction and the associated radiative forcing is definitely positive. Water-soluble carbonaceous matter dominated the carbon transport in the high-altitude glacial area. Owing to prevailing global warming and projected increase in carbon emission, the glacial DOC is expected to release, which will have strong underlying impacts on cryosphere ecosystem. The results of this work have profound implications for better understanding the carbon cycle in high altitude cryosphere regions. A new compilation of globally distributed work is required, including large-scale measurements of glacial DOC over high-altitude cryosphere regions, to overcome and address the scientific challenges to constrain climate impacts of light-absorbing impurities related processes in Earth system and climate models. [ABSTRACT FROM AUTHOR]

Glacier surface albedo dominates glacier energy balance, thus strongly affecting the glacier mass balance. Glaciers in the Western Nyainqentanglha Mountains (WNM) experienced large mass losses in the past two decades, but long-term changes of glacier albedo and its drivers are less understood. In this study, we retrieved glacier albedo with MODIS reflectance data to characterize the spatiotemporal variability of albedo from 2001 to 2020. Air temperature, rainfall, snowfall and deposition of light-absorbing impurities (LAIs) were evaluated as potential drivers of the observed variability in glacier albedo. The results showed that: (1) the glacier albedo experienced large inter-annual fluctuations, with the mean albedo being 0.552 ± 0.002 and a clear decreasing trend of 0.0443 ± 2 × 10−4 dec−1 in the WNM. The fastest decline was observed in autumn and in the vicinity of the equilibrium line altitude, indicating an extended melt season and an expansion of the ablation region to higher elevation; (2) local meteorology and LAIs deposition are the main drivers of glacier albedo change, but their effects on seasonal albedos are different due to different glacier processes. Both air temperature and the balance between liquid and solid precipitation affect summer and autumn albedos due to glacier ablation. Air temperature is the main driver of spring and winter albedos due to sublimation and metamorphism of snow, while snowfall carried by westerlies has limited influence on these two seasonal albedos due to less snowfall. LAIs mainly affect spring albedo due to high concentration coupled with the southerly wind in spring. These findings highlight the significance of changes in glacier albedo and the key role of local meteorology and LAIs deposition in determining such changes, which play an important role in glaciological and cryosphere processes.

Using global Sentinel-1 radar backscatter data, we systematically map the locations of glaciers with surge-type activity during 2017–22. Patterns of pronounced increases or decreases in the strongest backscatter between two winter seasons often indicate large changes in glacier crevassing, which we treat here as a sign of surge-type activity. Validations against velocity time series, terminus advances and crevassing found in optical satellite images confirm the robustness of this approach. We find 115 surge-type events globally between 2017 and 2022, around 100 of which on glaciers already know as surge-type. Our data reveal a pronounced spatial clustering in three regions, (i) Karakoram, Pamirs and Western Kunlun Shan (~50 surges), (ii) Svalbard (~25) and (iii) Yukon/Alaska (~9), with only a few other scattered surges elsewhere. This spatial clustering is significantly more pronounced than the overall global clustering of known surge-type glaciers. The 2017–22 clustering may point to climatic forcing of surge initiation.

Measuring surface mass-balance in the accumulation areas of glaciers is challenging because of the high spatial variability of snow accumulation and the difficulty of conducting annual field glaciological measurements. Here, we propose a method that can solve both these problems for many locations. Ground-penetrating radar measurements and firn cores extracted from a site in the French Alps were first used to reconstruct the topography of a buried end-of-summer snow horizon from a past year. Using these data and surface elevation observations from LiDAR and Global Navigation Satellite System instruments, we calculated the submergence velocities over the period between the buried horizon and more recent surface elevation observations. The differences between the changes in surface elevation and the submergence velocities were then used to calculate the annual surface mass-balances with an accuracy of ±0.34 m w.e. Assuming that the submergence velocities remain stable over several years, the surface mass-balance can be reconstructed for subsequent years from the differences in surface elevation alone. As opposed to the glaciological method that requires substantial fieldwork year after year to provide only point observations, this method, once submergence velocities have been calculated, requires only remote-sensing data to provide spatially distributed annual mass-balances in accumulation areas.

Initialising model glaciers such that they match well with their real counterparts and are thus able to make more accurate predictions is an ongoing challenge in glacier modelling. We set out a data-assimilation approach using an ensemble Kalman filter in a 2D flowline example that provides one possible solution to this problem. We show that our approach is valid across a range of parameters and scenarios, including deliberately data-deficient or inaccurate ones, and leads to robust retrieval of the glacier bed. We also provide some suggestions for how best to use data assimilation within a mountain-glacier context.

The conducted studies allowed us to obtain a detailed picture of glaciation changes in the mountain regions of Russia, most provided with information about glaciers in the twentieth century. For the Caucasus glaciation, the data were obtained for the time periods 1911, 1952, 2000, 2014, 2018 and 2020; for Altai – for 1850, 1952, 2003 and 2018. For large glaciation nodes of the Central Altai – Katun, South and North Chuya Ridges additionally for 1968, 2008, and 2017. In both areas, a decrease in the area of glaciers since the beginning of the twentieth century and acceleration of the rate of reduction in the early twenty-first century have been revealed. The glaciers of the Caucasus and Altai reduced their size during this time by 46% and 48%, respectively. On average, in the twentieth century the glaciers of the Caucasus lost about 0.2% of the area per year, in the Altai – 0.15%, and in the early twentieth century 1.15 and 1.7% respectively. To study Kamchatka glaciers, we used data from the Glacier Inventory of the USSR (1950/1957) and images from different satellites in the period 2007–2019. The glaciers that were not previously registered in the Glacier Inventory of the USSR were identified. The greatest number of such glaciers is in the northern part of the Midnight Ridge, where out of 465 glaciers identified on modern satellite images, 216 were not included in the Glacier Inventory of the USSR. The area of glaciation in different regions of Kamchatka has changed extremely unevenly since the first cataloguing, which is associated with significant differences in glacier morphology. Glaciers of volcanic areas increased their size or remained stationary; here there is no tendency for glaciers to decrease due to the thick surface moraine consisting of volcanogenic material. Comparison of data from the Glacier Inventory of the USSR (as of the 1950s–70s) and the Glacier Inventory of Russia (2017–2019) shows a decrease in glaciation area from the mid-20th century to the end of the second decade of the 21st century in all mountainous regions of Russia except only the volcanic regions of Kamchatka. The area reduction ranges from 63% (Ural) to 13% (Kodar). The largest glacial systems of the Caucasus, Kamchatka and Altai have reduced their areas by 25, 22 and 39%, respectively.

Melting and calving of glaciers and ice caps in Antarctica and Greenland could potentially contribute significantly to global sea level rise. Updates to existing outlines that provide critical glacier baseline information in both regions could help in the analysis of particular changes in glacier parameters such as area and volume from time-series inventories. Here we synthesize previously established techniques and apply new multi-source datasets to update glacier outlines in selected test areas of Antarctica and Greenland, as well as to reduce uncertainties and errors during the mapping process. The workflow includes mapping glacier boundaries, subdividing glaciers by watersheds and assigning glacier attributes. Complicated glacier scenarios and updating challenges in polar regions are discussed and demonstrated by representative case studies. For the first time in Antarctica, we analyze the effect of terminus types on mapped glacier areas, and in Greenland we compare the differences with glacier mapping results using Landsat OLI and ETM+. With new data sources, the methods described in this study might help to create glacier outlines on a larger scale in Antarctica and Greenland. Although data sources can be substituted, the enormous amount of manual labor required to update glacier inventories remains a significant challenge.

The uncertainty of glacier change projections is largely influenced by glacier models. In this study, we focus on temperature-index mass-balance (MB) models and their calibration. Using the Open Global Glacier Model (OGGM), we examine the influence of different surface-type dependent degree-day factors, temporal climate resolutions (daily, monthly) and downscaling options (temperature lapse rates, temperature and precipitation corrections) for 88 glaciers with in-situ observations. Our findings indicate that higher spatial and temporal resolution observations enhance MB gradient representation due to an improved calibration. The addition of surface-type distinction in the model also improves MB gradients, but the lack of independent observations limits our ability to demonstrate the added value of increased model complexity. Some model choices have systematic effects, for example weaker temperature lapse rates result in smaller projected glaciers. However, we often find counter balancing effects, such as the sensitivity to different degree-day factors for snow, firn and ice, which depends on how the glacier accumulation area ratio changes in the future. Similarly, using daily versus monthly climate data can affect glaciers differently depending on the shifting balance between melt and solid precipitation thresholds. Our study highlights the importance of considering minor model design differences to predict future glacier volumes and runoff accurately.

Decadal-scale, high-resolution geodetic measurements of glacier thinning have transformed our understanding of glacier response to climate change. Annual glacier mass balance can be estimated using remote-sensing proxies like snow-line altitude. These methods require field data for calibration, which are not available for most glaciers. Here we propose a method that combines multiple remotely-sensed proxies to obtain robust estimates of the annual glacier-wide balance using only remotely-sensed decadal-scale geodetic mass balance for calibration. The method is tested on Chhota Shigri, Argentière and Saint-Sorlin glaciers in the Himalaya and the Alps between 2001 and 2020, using four remotely-sensed proxies – the snow-line altitude, the minimum summer albedo over the glacier and two statistics of normalised difference snow index over the off-glacier area around the ablation zone. The reconstructed mass balance compares favourably with the corresponding glaciological field data (correlation coefficient 0.81 − 0.90, p < 0.001; root mean squared error 0.38 − 0.43 m w.e. a−1). The method presented may be useful to study interannual variability in mass balance on glaciers where no field data are available.

Shallow firn cores, in addition to a near-basal ice core, were recovered in 2018 from the Quelccaya ice cap (5470 m a.s.l) in the Cordillera Vilcanota, Peru, and in 2017 from the Nevado Illimani glacier (6350 m a.s.l) in the Cordillera Real, Bolivia. The two sites are ~450 km apart. Despite meltwater percolation resulting from warming, particle-based trace element records (e.g. Fe, Mg, K) in the Quelccaya and Illimani shallow cores retain well-preserved signals. The firn core chronologies, established independently by annual layer counting, show a convincing overlap indicating the two records contain comparable signals and therefore capture similar regional scale climatology. Trace element records at a ~1–4 cm resolution provide past records of anthropogenic emissions, dust sources, volcanic emissions, evaporite salts and marine-sourced air masses. Using novel ultra-high-resolution (120 μm) laser technology, we identify annual layer thicknesses ranging from 0.3 to 0.8 cm in a section of 2000-year-old radiocarbon-dated near-basal ice which compared to the previous annual layer estimates suggests that Quelccaya ice cores drilled to bedrock may be older than previously suggested by depth-age models. With the information collected from this study in combination with past studies, we emphasize the importance of collecting new surface-to-bedrock ice cores from at least the Quelccaya ice cap, in particular, due to its projected disappearance as soon as the 2050s.

Despite their importance for regional water resource planning and as indicators of climate change, records of in situ glacier mass balance remain short and spatially sparse in the Himalaya. Here, we present an updated series of in situ mass-balance measurements from Rikha Samba Glacier, Nepal, between 2011 and 2021. The updated in situ mass balance is −0.39 ± 0.32 m w.e. for this period. We use an energy-mass balance model to extend the annual mass-balance series back to 1974. The model is forced using daily meteorological variables from ERA5-Land reanalysis data that is linearly bias-corrected using observations from an automatic weather station situated near the glacier terminus. The modeled mass balance is consistent with the in situ mass-balance series measured 2011–2021 and with previous glaciological and geodetic estimates. The model results indicate a mass balance of −0.56 ± 0.27 m w.e. a−1 over the reconstruction period of 1974–2021, which is comparable to the mass losses experienced by other Himalayan glaciers during this time. An assessment of the sensitivity of the glacier mass balance to meteorological forcing suggests that a change in temperature of ±1 K has a stronger effect on the calculated mass balance compared to a ±20% change in either precipitation, or relative humidity, or solar radiation.

Although most processes governing the surface mass balance on mountain glaciers are well understood, the causes and extent of spatial variability in accumulation remain poorly constrained. In the present study, we couple an energy balance–snow and firn mass-balance model to terrain-based modelling routines estimating mass redistribution by snowdrift, preferential deposition and avalanching. We find this newly coupled model improves the spatial accuracy of winter balance simulations on Storglaciären, Sweden, while retaining versatility and a low computational cost. Accumulation on Storglaciären is primarily driven by direct precipitation, which is locally increased due to small-scale orographic effects. Wind-driven snow transport leads to substantial deposition in the accumulation zone and slight erosion in the ablation zone. Avalanching is the smallest contributor to winter balance, but cannot be neglected. The role of mass transporting processes in maintaining the current mass equilibrium on Storglaciären highlights the necessity to understand the links between climatic predictors and accumulation in order to accurately assess climate sensitivity.

Characterising the structures within glaciers can give unique insight into ice motion processes. On debris-covered glaciers, traditional structural glaciological mapping is challenging because the lower glacier is hidden by the supraglacial debris layer. Here, we use high-resolution optical televiewer (OPTV) image logs from four boreholes drilled into Khumbu Glacier, Nepal, to overcome this limitation and investigate englacial structural features within a Himalayan debris-covered glacier. The OPTV logs show structural features that are up to an order of magnitude thinner than those observed at the glacier surface and reveal five structural units: (I) primary stratification of ice; (II) debris-rich planes that conform with the primary stratification; (III) water-healed crevasse traces; (IV) healed crevasse traces; and (V) steeply dipping planes of basally derived fine sediment near the glacier terminus. The OPTV logs also reveal that the primary stratification both decreases in dip with depth (by up to 56° over 20 m) and rotates with depth (by up to 100° over 20 m) towards parallelism with the proximal lateral moraine. This transformation and the presence of relict layers of basally derived sediment raised into an englacial position – possibly involving thrusting – near the glacier's now stagnant terminus reveal a previously more dynamic glacier regime.

Cryospheric change occurs in unequal spaces. Societies living near ice are divided by race, class, gender, geography, politics and other factors. Consequently, impacts of ice loss are not shared equally, and everyone experiences cryospheric changes differently. Responsibility for recent ice loss is also driven by a relatively small portion of humanity: those who emit the most greenhouse gases. Additionally, people who study the cryosphere come from institutions and societies where inequality is often systemic, making research on ice and snow a symptom of and contributor to social inequality. To better understand unequal effects of cryospheric change within and across diverse communities, including research communities, this paper focuses on three areas, drawing primarily from glacier-related work: (1) the social context of cryospheric changes; (2) attribution and responsibility for cryospheric changes and (3) imbalances in knowledge about the cryosphere. Addressing these dimensions of ice loss requires transdisciplinary approaches that connect research to societies and link glaciology and other cryospheric sciences with social sciences and humanities. These concepts, cases and suggestions to help address inequalities also reveal that no singular conceptualization of sustainability exists. Different societies, residents and researchers possess distinct understandings of and goals for ‘ice in a sustainable society’.

Global climatic changes in recent decades cause changes in the dynamics of mountain geosystems. Of interest is the response of nival-glacial inland geosystems, which is different from those in the pre-Oceanic regions. In 2013, the authors substantiated the expediency of separating the nival-glacial mountain range geosystems within the Baikal Rift Zone and Mongolian Altai into a meridional research transect, within which it is possible to establish latitudinal patterns of modern intracontinental glacier functioning. Over the past decade of ongoing studies of glacial forms within the transect, many new data have been collected and analyzed. The present article based on systematic generalizations of these materials. The analysis of the state of mountain glaciation in the south of Eastern Siberia and Mongolia in the last decades has revealed a number of peculiarities. The deglaciation that began in the 1970s–1980s continues in all the mountain systems under consideration. It is a reaction to global climatic changes. The greatest reduction is experienced by the near-slope glaciers and flattop glaciers and to a lesser extent by the cirque ice glaciers. They have a significant volume loss due to thinning, while changes in the area of such glaciers are less significant. This is the main difference between the intracontinental glaciers and the preoceanic glaciers. In general, many nival-glacial geosystems of transitional forms are formed. Quantitative characteristics are based on the calculation of the volume of the East Sayan glaciers measured by geophysical radar. Thus, over 120 years, the Peretolchina (northern) Glacier has decreased in length by a factor of 1.75, in area by a factor of 2.9, and in volume by a factor of 3.71. During the same period, the area of the neighboring Radde Glacier decreased from 0.43 to 0.09 km2; its thickness by 30 m; and, accordingly, its volume by 3 times. The warming of air temperature for Western Mongolia is 0.03–0.29°C/10 years and, for the Baikal region, 0.2–0.5°C/10 years. Warming of air temperatures in the ridges of the Mongolian Altai is noted up to 48° N, and southward the trend is unstable. [ABSTRACT FROM AUTHOR]

We present results from a 10-year study of the local features of mountain glacier dynamics in Central Asia—from the latitudes of the middle taiga of Cisbaikalia and Transbaikalia, through the mountains of the steppe zone of Mongolian Altai and the desert zone of China to the cold deserts of the Himalayas. Multitemporal satellite images were interpreted, making it possible to obtain quantitative information on changes in the area and length of glaciers in key areas of the studied mountain territories. These data are linked to climatic parameters, absolute altitude, permafrost, morphology and tectonics of the mountain ridges, and anthropogenic impact. The study revealed a general trend: the retreat of glaciers starting in the 1970s, with an increase in the average annual air temperature in the Northern Hemisphere by 1.4‒1.6°С. Accelerated glacier melting and a reduction in glacier length were observed from 1990 to 2020, when the most significant warming occurred. The following years with above-average temperatures are highlighted: 1990, 1995, 1997, 2002, 2005, 2005, 2007, 2008, 2011, 2012, 2016, 2019, and 2020. This period also witnessed an increase in the area and length of some glaciers studied in 2013‒2015 and in 2021 associated with abnormal winter precipitation in 2012, 2013, and 2020. It was found that in the mountains in the south of Eastern Siberia, the aforementioned years showed a slowdown in the reduction of the area of glaciers, an increase of the number of perennial snow patches, and an increase in the frequency of snow avalanches. It has been revealed that the glaciers in Siberia, Mongolian Altai, and the Altyn-Tag Range retreated from the terminal moraines of the Little Ice Age by an average of 500 m; on the Karlyktag Range, 800‒900 m; and in the Himalayas (Langtang area), from 1 to 4.5 km. [ABSTRACT FROM AUTHOR]

Glaciers are important indictors of climate change as changes in glaciers physical features such as their area is in response to measurable evidence of fluctuating climate factors such as temperature, precipitation, and CO2. Although a general retreat of mountain glacier systems has been identified in relation to centennial trends toward warmer temperatures, there is the potential to extract a great deal more information regarding regional variations in climate from the mapping of the time history of the terminus position or surface area of the glaciers. The remote nature of glaciers renders direct measurement impractical on anything other than a local scale. Considering the sheer number of mountain glaciers around the globe, ground measurements of terminus position are only available for a small percentage of glaciers and ground measurements of glacier area are rare. In this project, changes in the terminal point and area of Franz Josef and Gorner glaciers were quantified in response to climate factors using satellite imagery taken by Landsat at regular intervals. Two supervised learning methods including a parametric method (multiple regression) and a nonparametric method (generalized additive model) were implemented to identify climate factors that impact glacier changes. Local temperature, CO2, and precipitation were identified as significant factors for predicting changes in both Franz Josef and Gorner glaciers. Spatiotemporal quantification of glacier change is an essential task to model glacier variations in response to global and local climate factors. This work provided valuable insights on quantification of surface area of glaciers using satellite imagery with potential implementation of a generic approach. [ABSTRACT FROM AUTHOR]

Using global Sentinel-1 radar backscatter data, we systematically map the locations of glaciers with surge-type activity during 2017–22. Patterns of pronounced increases or decreases in the strongest backscatter between two winter seasons often indicate large changes in glacier crevassing, which we treat here as a sign of surge-type activity. Validations against velocity time series, terminus advances and crevassing found in optical satellite images confirm the robustness of this approach. We find 115 surge-type events globally between 2017 and 2022, around 100 of which on glaciers already know as surge-type. Our data reveal a pronounced spatial clustering in three regions, (i) Karakoram, Pamirs and Western Kunlun Shan (~50 surges), (ii) Svalbard (~25) and (iii) Yukon/Alaska (~9), with only a few other scattered surges elsewhere. This spatial clustering is significantly more pronounced than the overall global clustering of known surge-type glaciers. The 2017–22 clustering may point to climatic forcing of surge initiation. [ABSTRACT FROM AUTHOR]

Glacier surface albedo dominates glacier energy balance, thus strongly affecting the glacier mass balance. Glaciers in the Western Nyainqentanglha Mountains (WNM) experienced large mass losses in the past two decades, but long-term changes of glacier albedo and its drivers are less understood. In this study, we retrieved glacier albedo with MODIS reflectance data to characterize the spatiotemporal variability of albedo from 2001 to 2020. Air temperature, rainfall, snowfall and deposition of light-absorbing impurities (LAIs) were evaluated as potential drivers of the observed variability in glacier albedo. The results showed that: (1) the glacier albedo experienced large inter-annual fluctuations, with the mean albedo being 0.552 ± 0.002 and a clear decreasing trend of 0.0443 ± 2 × 10−4 dec−1 in the WNM. The fastest decline was observed in autumn and in the vicinity of the equilibrium line altitude, indicating an extended melt season and an expansion of the ablation region to higher elevation; (2) local meteorology and LAIs deposition are the main drivers of glacier albedo change, but their effects on seasonal albedos are different due to different glacier processes. Both air temperature and the balance between liquid and solid precipitation affect summer and autumn albedos due to glacier ablation. Air temperature is the main driver of spring and winter albedos due to sublimation and metamorphism of snow, while snowfall carried by westerlies has limited influence on these two seasonal albedos due to less snowfall. LAIs mainly affect spring albedo due to high concentration coupled with the southerly wind in spring. These findings highlight the significance of changes in glacier albedo and the key role of local meteorology and LAIs deposition in determining such changes, which play an important role in glaciological and cryosphere processes. [ABSTRACT FROM AUTHOR]

Initialising model glaciers such that they match well with their real counterparts and are thus able to make more accurate predictions is an ongoing challenge in glacier modelling. We set out a data-assimilation approach using an ensemble Kalman filter in a 2D flowline example that provides one possible solution to this problem. We show that our approach is valid across a range of parameters and scenarios, including deliberately data-deficient or inaccurate ones, and leads to robust retrieval of the glacier bed. We also provide some suggestions for how best to use data assimilation within a mountain-glacier context. [ABSTRACT FROM AUTHOR]

Measuring surface mass-balance in the accumulation areas of glaciers is challenging because of the high spatial variability of snow accumulation and the difficulty of conducting annual field glaciological measurements. Here, we propose a method that can solve both these problems for many locations. Ground-penetrating radar measurements and firn cores extracted from a site in the French Alps were first used to reconstruct the topography of a buried end-of-summer snow horizon from a past year. Using these data and surface elevation observations from LiDAR and Global Navigation Satellite System instruments, we calculated the submergence velocities over the period between the buried horizon and more recent surface elevation observations. The differences between the changes in surface elevation and the submergence velocities were then used to calculate the annual surface mass-balances with an accuracy of ±0.34 m w.e. Assuming that the submergence velocities remain stable over several years, the surface mass-balance can be reconstructed for subsequent years from the differences in surface elevation alone. As opposed to the glaciological method that requires substantial fieldwork year after year to provide only point observations, this method, once submergence velocities have been calculated, requires only remote-sensing data to provide spatially distributed annual mass-balances in accumulation areas. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Theory and observation show that glacier-flow regimes characterized by high basal slip enhance the projection of topographic detail to the surface, motivating this investigation into the efficacy of using glacier surges to improve bed estimation. Here we adapt a Bayesian inversion scheme and apply it to real and synthetic data as a proof of concept. Synthetic tests show a reduction in mean RMSE between true and inferred beds by more than half, and an increase in the mean correlation coefficient of ~0.5, when data from slip- versus deformation-dominated regimes are used. Multi-epoch inversions, which partition slip- and deformation-dominated regimes, are shown to outperform inversions that average over these flow regimes thereby squandering information. Tests with real data from a surging glacier in Yukon, Canada, corroborate these results, while highlighting the challenges of limited or inconsistent data. With the growing torrent of satellite-based observations, fast-flow events such as glacier surges offer potential to improve bed estimation for some of the world's most dynamic glaciers.

The Tian Shan mountain range, known as the water towers of Central Asia, plays a key role in local water supply, yet large uncertainties remain about the amount of water that is stored in its glaciers. In this study, we assess the impact of the boundary conditions on ice thickness estimates using two inversion models: a mass conservation (MC) model and a basal shear stress (BS) model. We compare the widely used Randolph Glacier Inventory version 6 with the updated Glacier Area Mapping for Discharge from the Asian Mountains glacier inventory, as well as two digital elevation models (SRTM DEM and Copernicus DEM). The results show that the ice volume (in ~2000 CE) in the Tian Shan range is 661.0 ± 163.5 km3 for the MC model and 552.8 ± 85.3 km3 for the BS model. There are strong regional differences due to inventory, especially for glaciers in China (17–25%). However, the effect of different DEM sources on ice volume estimation is limited. By the end of the 21st century, the projected mass loss differences between inventories are higher than between adjacent emission scenarios, illustrating the vital importance of high-quality inventories. These differences should be carefully considered during water resource planning.

On the Kamchatka Peninsula in the Russian Far East, 405 glaciers with an estimated total mass of 49 Gt were reported in the 1970s. These have been retreating at an accelerated rate since the start of the 21st century. Because glacier studies in this region are scarce, ice loss and its influence on sea level rise and regional environments is poorly understood. In this study, we analyzed satellite data to quantify glacier mass change from 2000 to 2016 in six major glacier-covered regions on the peninsula. The mean rate of the glacier mass change over the study period was −0.46 ± 0.01 m w.e. a−1 (total mass change was −4.9 ± 0.1 Gt, −304.2 ± 9.1 Mt a−1), which is slightly lower than other regions in mid-latitude and subarctic zones. The mass loss accelerated from >−0.33 ± 0.02 m w.e. a−1 in the period 2000–2006/2010 to

Glacier change in Kenai Fjords National Park in southcentral Alaska affects local terrestrial, fresh water and marine ecosystems and will likely impact ecotourism. We used Landsat 4–8 imagery from 1984 through 2021 to manually map lower glacier ice margins for 19 maritime glaciers in Kenai Fjords National Park. Of these glaciers, six are tidewater, three are lake-terminating, six are land-terminating and four terminated in more than one environment throughout the study period. We used the mapped ice margins to quantify seasonal terminus position and areal change, including distinguishing between ice loss at glacier termini and along glacier margins. Overall, 13 glaciers substantially retreated (more than 2σ), 14 lost substantial area and only two underwent both net advance and area gain. The glaciers that had insubstantial length and area changes were predominantly tidewater. Cumulatively, the lower reaches of these 19 glaciers lost 42 km2 of ice, which was nearly evenly distributed between the terminus and the lateral margins. The rapid rate of glacier change and subsequent land cover changes are highly visible to visitors and locals at Kenai Fjords National Park, and this study quantifies those changes in terms of glacier length and area.

Sea-level projections depend sensitively on the parameterization used for basal slip in glacier flow models. During slip over rock-beds, ice-bed separation increases with slip velocity and basal water pressure. We present a method for using these variables and measured bed topography to estimate the average bed slope in contact with ice, ${\bar m}$. Three-dimensional numerical modeling of slip over small areas of former beds has shown that changes in ${\bar m}$ with increasing slip velocity and water pressure mimic changes in basal drag. Computed values of ${\bar m}$ can thus provide the form of the slip law that relates drag to velocity and water pressure, avoiding computationally expensive numerical modeling. The method is applied to 618 sections from four former glacier beds. Results generally show an increase in ${\bar m}$, and hence inferred basal drag, with slip velocity up to a limiting value, consistent with a regularized Coulomb slip law.

Glacial hydraulic tremor (GHT) can be monitored to observe changes in location and distribution of water flow beneath glacial ice, allowing the spatiotemporal evolution of subglacial hydrology to be studied continuously and remotely. We use frequency-dependent polarization analysis (FDPA) to classify types of GHT and assess its spatio-temporal extent beneath Rhonegletscher, Switzerland, in a continuous seismic record through 2018 and 2019 at three ice-proximal bedrock-based seismometers. We determine the frequency bands composing the GHT and calculate back azimuth angles pointing to a previously known subglacial channel. We also inspect the relationship between GHT seismic power and water discharge from the glacier to compare daily and seasonal shifts with the observed GHT. We observed the seasonal shift from a distributed system to a channelized system, and our dataset allowed comparison of channel locations within and across seasons, with implications for sediment evacuation and bed erosion. The successful use of this method to assess GHT previously on Taku glacier (the methods of which this project follows) and now Rhonegletscher shows that existing ice-proximal passive seismic installations can be used to easily and continuously monitor subglacial hydraulic activity.

We quantify the surface elevation changes along Yala Glacier in Langtang Valley, Nepal Himalaya, since 1981 using geodetic methods to understand the recent evolution and current state of small debris-free glaciers across the region. We analyse differential global positioning system measurements and aerial stereo imagery that were acquired along Yala Glacier in 2007, 2009, 2012 and 2015 to generate digital elevation models for each calculation period. Continuous surface lowering has mainly been observed across the down-glacier area during the calculation periods, although a large degree of variability exists, with this lowering trend propagating up-glacier in recent years. The area-weighted glacier mass balances range from −0.98 ± 0.27 to −0.26 ± 0.30 m w.e. a−1 for the five calculation periods (1981–2007, 2007–2009, 2009–2012, 2012–2015 and 2007–2015). These calculated mass-balance data reveal that Yala Glacier has undergone accelerated mass loss since the late 2000s, which is consistent with the results of previous in situ measurement and remote-sensing studies.

Melting and calving of glaciers and ice caps in Antarctica and Greenland could potentially contribute significantly to global sea level rise. Updates to existing outlines that provide critical glacier baseline information in both regions could help in the analysis of particular changes in glacier parameters such as area and volume from time-series inventories. Here we synthesize previously established techniques and apply new multi-source datasets to update glacier outlines in selected test areas of Antarctica and Greenland, as well as to reduce uncertainties and errors during the mapping process. The workflow includes mapping glacier boundaries, subdividing glaciers by watersheds and assigning glacier attributes. Complicated glacier scenarios and updating challenges in polar regions are discussed and demonstrated by representative case studies. For the first time in Antarctica, we analyze the effect of terminus types on mapped glacier areas, and in Greenland we compare the differences with glacier mapping results using Landsat OLI and ETM+. With new data sources, the methods described in this study might help to create glacier outlines on a larger scale in Antarctica and Greenland. Although data sources can be substituted, the enormous amount of manual labor required to update glacier inventories remains a significant challenge. [ABSTRACT FROM AUTHOR]

The uncertainty of glacier change projections is largely influenced by glacier models. In this study, we focus on temperature-index mass-balance (MB) models and their calibration. Using the Open Global Glacier Model (OGGM), we examine the influence of different surface-type dependent degree-day factors, temporal climate resolutions (daily, monthly) and downscaling options (temperature lapse rates, temperature and precipitation corrections) for 88 glaciers with in-situ observations. Our findings indicate that higher spatial and temporal resolution observations enhance MB gradient representation due to an improved calibration. The addition of surface-type distinction in the model also improves MB gradients, but the lack of independent observations limits our ability to demonstrate the added value of increased model complexity. Some model choices have systematic effects, for example weaker temperature lapse rates result in smaller projected glaciers. However, we often find counter balancing effects, such as the sensitivity to different degree-day factors for snow, firn and ice, which depends on how the glacier accumulation area ratio changes in the future. Similarly, using daily versus monthly climate data can affect glaciers differently depending on the shifting balance between melt and solid precipitation thresholds. Our study highlights the importance of considering minor model design differences to predict future glacier volumes and runoff accurately. [ABSTRACT FROM AUTHOR]

Citizen science and related engagement programmes have proliferated in recent years throughout the sciences but have been reasonably limited in the cryospheric sciences. In the semiarid Andes we at the Centro de Estudios Avanzados en Zonas Áridas have developed a range of initiatives together with the wider community and stakeholder institutions to improve our understanding of the role snow and ice play in headwater catchments. In this paper we reflect on ongoing engagement with communities living and working in and near study sites of cryospheric science in northern Chile as a strategy that can both strengthen the research being done and empower local communities. [ABSTRACT FROM AUTHOR]

Mountain glaciers are crucial sources of fresh water, contributing directly and indirectly to water, energy and food supplies for hundreds of millions of people. Assessing the impact of diminishing glacial meltwater contributions to the security of this resource is critical as we seek to manage and adapt to changing freshwater dynamics in a warming world. Both water quantity and quality influence water (in)security, so understanding the fluxes of water, sediment and contaminants through glacial and proglacial systems is required for holistic assessment of meltwater contribution to downstream resource security. In this paper we consider the socio-environmental role of and pressures on glacier-fed waters, discuss key research priorities for the assessment of both the quantity and quality of meltwater and reflect on the importance of situating our understanding within a transdisciplinary and inclusive research landscape. [ABSTRACT FROM AUTHOR]

Basal sliding of glaciers and ice sheets remains a source of uncertainty in simulating the long‐term evolution of ice masses. In particular, the response of ice flow to changes in driving stress depends strongly on the value of the exponent m in nonlinear friction laws (e.g., Weertman's law), which is poorly constrained by observations. Here we constrain the friction law at a natural scale on Argentière Glacier (French Alps, hard‐bed), taking advantage of well‐resolved observations of glacier mass balance, geometry and basal sliding over time spans that include large changes in driving stress. By combining three different independent methods based on (a) surface velocity inversion, (b) transient length change modeling, and (c) direct local sliding measurements, we consistently find a value of m = 3.1 ± 0.3. We suggest that Weertman's law is suitable for modeling the long‐term evolution of hard‐bedded glaciers and ice sheets. Plain Language Summary: Modeling the evolution of glaciers and ice sheets under climate variability requires accurate estimation of ice flow velocities by numerical models. One of the most challenging components of these models is the representation of basal sliding at the rock‐ice interface, generally described by relationships between stress and sliding speed. These relationships are mainly derived from laboratory experiments and theoretical studies, and their ability to represent nature still needs to be evaluated. Because direct observation beneath glaciers and ice sheets is difficult, few studies have attempted to validate sliding models from natural scale observations. In this context, our study provides a rare constraint, based on observations on an alpine glacier, on the law that should be used to model glacier sliding on clean bedrock (so‐called hard‐bed glaciers). We show that a simple power law performs well in explaining long‐term glacier behavior for a power exponent of ∼3.1. We suggest that this exponent should be adopted in non‐linear power laws incorporated into ice flow models that perform future projections of hard‐bedded glaciers and ice sheet evolution. Key Points: We suggest that a Weertman friction law is suitable for multi‐decadal simulations of the hard‐bedded glaciersCombining three independent methods and data sets we infer a Weertman's friction exponent of 3.1 ± 0.3The temporal evolution of glacier fronts is found to be strongly sensitive to the friction law exponent [ABSTRACT FROM AUTHOR]

Characterising the structures within glaciers can give unique insight into ice motion processes. On debris-covered glaciers, traditional structural glaciological mapping is challenging because the lower glacier is hidden by the supraglacial debris layer. Here, we use high-resolution optical televiewer (OPTV) image logs from four boreholes drilled into Khumbu Glacier, Nepal, to overcome this limitation and investigate englacial structural features within a Himalayan debris-covered glacier. The OPTV logs show structural features that are up to an order of magnitude thinner than those observed at the glacier surface and reveal five structural units: (I) primary stratification of ice; (II) debris-rich planes that conform with the primary stratification; (III) water-healed crevasse traces; (IV) healed crevasse traces; and (V) steeply dipping planes of basally derived fine sediment near the glacier terminus. The OPTV logs also reveal that the primary stratification both decreases in dip with depth (by up to 56° over 20 m) and rotates with depth (by up to 100° over 20 m) towards parallelism with the proximal lateral moraine. This transformation and the presence of relict layers of basally derived sediment raised into an englacial position – possibly involving thrusting – near the glacier's now stagnant terminus reveal a previously more dynamic glacier regime. [ABSTRACT FROM AUTHOR]

Despite their importance for regional water resource planning and as indicators of climate change, records of in situ glacier mass balance remain short and spatially sparse in the Himalaya. Here, we present an updated series of in situ mass-balance measurements from Rikha Samba Glacier, Nepal, between 2011 and 2021. The updated in situ mass balance is −0.39 ± 0.32 m w.e. for this period. We use an energy-mass balance model to extend the annual mass-balance series back to 1974. The model is forced using daily meteorological variables from ERA5-Land reanalysis data that is linearly bias-corrected using observations from an automatic weather station situated near the glacier terminus. The modeled mass balance is consistent with the in situ mass-balance series measured 2011–2021 and with previous glaciological and geodetic estimates. The model results indicate a mass balance of −0.56 ± 0.27 m w.e. a−1 over the reconstruction period of 1974–2021, which is comparable to the mass losses experienced by other Himalayan glaciers during this time. An assessment of the sensitivity of the glacier mass balance to meteorological forcing suggests that a change in temperature of ±1 K has a stronger effect on the calculated mass balance compared to a ±20% change in either precipitation, or relative humidity, or solar radiation. [ABSTRACT FROM AUTHOR]

Although most processes governing the surface mass balance on mountain glaciers are well understood, the causes and extent of spatial variability in accumulation remain poorly constrained. In the present study, we couple an energy balance–snow and firn mass-balance model to terrain-based modelling routines estimating mass redistribution by snowdrift, preferential deposition and avalanching. We find this newly coupled model improves the spatial accuracy of winter balance simulations on Storglaciären, Sweden, while retaining versatility and a low computational cost. Accumulation on Storglaciären is primarily driven by direct precipitation, which is locally increased due to small-scale orographic effects. Wind-driven snow transport leads to substantial deposition in the accumulation zone and slight erosion in the ablation zone. Avalanching is the smallest contributor to winter balance, but cannot be neglected. The role of mass transporting processes in maintaining the current mass equilibrium on Storglaciären highlights the necessity to understand the links between climatic predictors and accumulation in order to accurately assess climate sensitivity. [ABSTRACT FROM AUTHOR]

Shallow firn cores, in addition to a near-basal ice core, were recovered in 2018 from the Quelccaya ice cap (5470 m a.s.l) in the Cordillera Vilcanota, Peru, and in 2017 from the Nevado Illimani glacier (6350 m a.s.l) in the Cordillera Real, Bolivia. The two sites are ~450 km apart. Despite meltwater percolation resulting from warming, particle-based trace element records (e.g. Fe, Mg, K) in the Quelccaya and Illimani shallow cores retain well-preserved signals. The firn core chronologies, established independently by annual layer counting, show a convincing overlap indicating the two records contain comparable signals and therefore capture similar regional scale climatology. Trace element records at a ~1–4 cm resolution provide past records of anthropogenic emissions, dust sources, volcanic emissions, evaporite salts and marine-sourced air masses. Using novel ultra-high-resolution (120 μm) laser technology, we identify annual layer thicknesses ranging from 0.3 to 0.8 cm in a section of 2000-year-old radiocarbon-dated near-basal ice which compared to the previous annual layer estimates suggests that Quelccaya ice cores drilled to bedrock may be older than previously suggested by depth-age models. With the information collected from this study in combination with past studies, we emphasize the importance of collecting new surface-to-bedrock ice cores from at least the Quelccaya ice cap, in particular, due to its projected disappearance as soon as the 2050s. [ABSTRACT FROM AUTHOR]

Decadal-scale, high-resolution geodetic measurements of glacier thinning have transformed our understanding of glacier response to climate change. Annual glacier mass balance can be estimated using remote-sensing proxies like snow-line altitude. These methods require field data for calibration, which are not available for most glaciers. Here we propose a method that combines multiple remotely-sensed proxies to obtain robust estimates of the annual glacier-wide balance using only remotely-sensed decadal-scale geodetic mass balance for calibration. The method is tested on Chhota Shigri, Argentière and Saint-Sorlin glaciers in the Himalaya and the Alps between 2001 and 2020, using four remotely-sensed proxies – the snow-line altitude, the minimum summer albedo over the glacier and two statistics of normalised difference snow index over the off-glacier area around the ablation zone. The reconstructed mass balance compares favourably with the corresponding glaciological field data (correlation coefficient 0.81 − 0.90, p < 0.001; root mean squared error 0.38 − 0.43 m w.e. a−1). The method presented may be useful to study interannual variability in mass balance on glaciers where no field data are available. [ABSTRACT FROM AUTHOR]

Cryospheric change occurs in unequal spaces. Societies living near ice are divided by race, class, gender, geography, politics and other factors. Consequently, impacts of ice loss are not shared equally, and everyone experiences cryospheric changes differently. Responsibility for recent ice loss is also driven by a relatively small portion of humanity: those who emit the most greenhouse gases. Additionally, people who study the cryosphere come from institutions and societies where inequality is often systemic, making research on ice and snow a symptom of and contributor to social inequality. To better understand unequal effects of cryospheric change within and across diverse communities, including research communities, this paper focuses on three areas, drawing primarily from glacier-related work: (1) the social context of cryospheric changes; (2) attribution and responsibility for cryospheric changes and (3) imbalances in knowledge about the cryosphere. Addressing these dimensions of ice loss requires transdisciplinary approaches that connect research to societies and link glaciology and other cryospheric sciences with social sciences and humanities. These concepts, cases and suggestions to help address inequalities also reveal that no singular conceptualization of sustainability exists. Different societies, residents and researchers possess distinct understandings of and goals for 'ice in a sustainable society'. [ABSTRACT FROM AUTHOR]

We present geodetic mass-balance estimates for ten glaciers (22.6 ± 1.1 km2) around Volcán Domuyo between 1962 and 2020 (and 46 glaciers covering 29 ± 1.5 km2 between 1984 and 2020), derived from airborne, ASTER and Pléiades imagery. Overall, we find a slightly negative mass balance (−0.15 ± 0.09 m w.e. a–1) for the entire 1962–2020 time span. A closer inspection of sub-periods reveals, however, an increasingly negative mass balance with time. The Domuyo glaciers shifted from a moderately positive mass balance of +0.28 ± 0.13 m w.e. a–1 between 1962 and 1984, to a strongly negative mass balance of −0.99 ± 0.19 m w.e. a–1 between 2012 and 2020. An increase in summer temperatures and a decrease in winter precipitation during the last four decades are likely drivers of the observed glacier changes. We support this finding by implementing a minimal glacier model, relying solely on monthly precipitation and air temperatures. The mass-balance evolution detected in Volcán Domuyo is consistent with other sites spread across the Central Andes, suggesting rapidly increasing glacier wastage impacts are occurring at a geographically wider scale.

We present spatially distributed seasonal and annual surface mass balances of Wolverine Glacier, Alaska, from 2016 to 2020. Our approach accounts for the effects of ice emergence and firn compaction on surface elevation changes to resolve the spatial patterns in mass balance at 10 m scale. We present and compare three methods for estimating emergence velocities. Firn compaction was constrained by optimizing a firn model to fit three firn cores. Distributed mass balances showed good agreement with mass-balance stakes (RMSE = 0.67 m w.e., r = 0.99, n = 41) and ground-penetrating radar surveys (RMSE = 0.36 m w.e., r = 0.85, n = 9024). Fundamental differences in the distributions of seasonal balances highlight the importance of disparate physical processes, with anomalously high ablation rates observed in icefalls. Winter balances were found to be positively skewed when controlling for elevation, while summer and annual balances were negatively skewed. We show that only a small percent of the glacier surface represents ideal locations for mass-balance stake placement. Importantly, no suitable areas are found near the terminus or in elevation bands dominated by icefalls. These findings offer explanations for the often-needed geodetic calibrations of glaciological time series.

Glacier motion responds dynamically to changing meltwater inputs, but the multi-decadal response of basal sliding to climate remains poorly constrained due to its sensitivity across multiple timescales. Observational records of glacier motion provide critical benchmarks to decode processes influencing glacier dynamics, but multi-decadal records that precede satellite observation and modern warming are rare. Here we present a record of motion in the ablation zone of Saskatchewan Glacier that spans seven decades. We combine in situ and remote-sensing observations to inform a first-order glacier flow model used to estimate the relative contributions of sliding and internal deformation on dynamics. We find a significant increase in basal sliding rates between melt-seasons in the 1950s and those in the 1990s and 2010s and explore three process-based explanations for this anomalous behavior: (i) the glacier surface steepened over seven decades, maintaining flow-driving stresses despite sustained thinning; (ii) the formation of a proglacial lake after 1955 may support elevated basal water pressures; and (iii) subglacial topography may cause dynamic responses specific to Saskatchewan Glacier. Although further constraints are necessary to ascertain which processes are of greatest importance for Saskatchewan Glacier's dynamic evolution, this record provides a benchmark for studies of multi-decadal glacier dynamics.

The Desert Andes contain >4500 ice masses, but only a handful are currently being monitored. We present the mass changes of the small mountain glacier Agua Negra (1 km2) and of the rest of glaciers in the Jáchal river basin. Remote-sensing data show Agua Negra glacier lost 23% of its area during 1959–2019. Glaciological measurements during 2014–2021 indicate an average annual mass balance of −0.52 m w.e. a−1, with mean winter and summer balances of 0.80 and −1.33 m w.e. a−1, respectively. The Equilibrium Line Altitude (ELA) is estimated to be 5100 ± 100 m a.s.l., which corresponds to an Accumulation Area Ratio (AAR) of 0.28 ± 0.21. Geodetic data from SRTM X and Pléiades show a doubling of the loss rate from −0.32 ± 0.03 m w.e. a−1 in 2000–2013, to −0.66 ± 0.06 m w.e. a−1 in 2013–2019. Comparatively, the ice losses for the entire Jáchal river basin (25 500 km2) derived from ASTER show less negative values, −0.11 ± 16 m w.e. a−1 for 2000–2012 and −0.23 ± 14 m w.e. a−1 for 2012–2018. The regional warming trend since 1979 and a recent decline in snow accumulation are probably driving the observed glacier mass balance.

An end of summer snowline (EOSS) photographic dataset for Aotearoa New Zealand contains over four decades of equilibrium line altitude (ELA) observations for more than 50 index glaciers. This dataset provides an opportunity to create a climatological ELA reference series that has several applications. Our work screened out EOSS sites that had low temporal coverage and also removed limited observations when the official survey did not take place. Snowline data from 41 of 50 glaciers in the EOSS dataset were retained and included in a normalised master snowline series that spans 1977–2020. Application of the regionally representative normalised master snowline series in monthly and seasonally resolved climate response function analyses showed consistently strong relationships with austral warm-season temperatures for land-based stations west of the Southern Alps and the central Tasman Sea. There is a trend towards higher regional snowlines since the 1990s that has been steepening in recent decades. If contemporary decadal normalised master snowline series trends are maintained, the average Southern Alps snowline elevation will be displaced at least 200 m higher than normal by the 2025–2034 decade. More frequent extremely high snowlines are expected to drive more extreme cumulative mass-balance losses that will reduce the glacierised area of Aotearoa New Zealand.

We investigate relationships between synoptic-scale atmospheric variability and the mass-balance of 13 Andean glaciers (located 16–55° S) using Pearson correlation coefficients (PCCs) and multiple regressions. We then train empirical glacier mass-balance models (EGMs) in a cross-validated multiple regression procedure for each glacier. We find four distinct glaciological zones with regard to their climatic controls: (1) The mass-balance of the Outer Tropics glaciers is linked to temperature and the El Niño-Southern Oscillation (PCC ⩽ 0.6), (2) glaciers of the Desert Andes are mainly controlled by zonal wind intensity (PCC ⩽ 0.9) and the Antarctic Oscillation (PCC ⩽0.6), (3) the mass-balance of the Central Andes glaciers is primarily correlated with precipitation anomalies (PCC ⩽ 0.8), and (4) the glacier of the Fuegian Andes is controlled by winter precipitation (PCC ≈ 0.7) and summer temperature (PCC ≈ −0.9). Mass-balance data in the Lakes District and Patagonian Andes zones, where most glaciers are located, are too sparse for a robust detection of synoptic-scale climatic controls. The EGMs yield R2 values of ~ 0.45 on average and ⩽ 0.74 for the glaciers of the Desert Andes. The EGMs presented here do not consider glacier dynamics or geometry and are therefore only suitable for short-term predictions.

Mountain glaciers, including the poorly examined glaciers of Lake Karakul Basin, are considered the most vulnerable part of the cryosphere directly reacting to the changing climatic conditions. The most recent aerial photographing and single field surveys of the glaciation in Lake Karakul Basin wre carried out as late as in 1953. The investigation of the basin glaciers therefore represents great scientific importance, including in terms of glaciology and climatology. The article describes the calculations of the surface mass balance(S MB) of Glacier No.139 in Lake Karakul Basin allowing, based on the findings,t o assess the degree of climate change impact on the glaciation in the target zone,as well as the possible changes in the state of glaciers in the future. Durin g2018-2019, the mass balance of Glacier No.139 in Lake Karakul Basin amountedt o -0.26 m WE. Glacier No.139 is located in the basin’ ssouthwestern part. T he data obtained from the Karakul Weather Station allowed establishing that the climate in the target watershed is harshly cold with little snow in winter an dmild summer, and the bulk of precipitation occurring in the warm season. Due tothe fact that in the past and until recently no glacier mass observations wereco nducted in the target zone, the research findings can serve as initial data t oexpand the scientific and applied knowledge, as well as contribute to increasing thea ccuracy of glacier dynamics modeling in the area.findings can serve as initial data to expand the scientific and applied knowledge, as well as contribute to increasing the accuracy of glacier dynamics modeling in the area. the dynamics of glaciers in this region.

Regional-scale surface energy balance (SEB) models of glacier melt require forcing by coarse-gridded data from reanalysis and/or global climate models that need to be downscaled to glacier scale. As on-glacier meteorological observations are rare, it generally remains unknown how exact the reanalysis and downscaled data are for the local-scale SEB modeling. We address this question by evaluating the performance of reanalysis from the European Centre for Medium-Range Weather Forecasts (ERA5 and ERA5-Land reanalysis), with and without downscaling, at four glaciers in Western Canada with available onglacier meteorological measurements collected over different summer seasons. We dynamically downscale ERA5 with the Weather Research and Forecasting (WRF) model at 3.3 km and 1.1 km grid spacing. We find that the SEB model, forced separately with the observations and the two reanalyses, yields <10%difference in simulated total melt energy and shows strong correlations (>0.79) in simulated timeseries of daily melt energy at each site. The good performance of the reanalysis-derived melt energy is partly due to cancellation of biases between overestimated incoming shortwave radiation and substantially underestimated wind speed and subsequently turbulent heat fluxes. Downscaling with WRF improves the simulation of wind speed, while other meteorological variables show similar performance to ERA5 without downscaling. The choice of WRF physics parameterization schemes is shown to have a relatively large impact on the simulations of SEB components, but a smaller impact on the modeled total melt energy. The results increase our confidence in dynamical downscaling with WRF for long-term glacier melt modeling in this region. [ABSTRACT FROM AUTHOR]