Bedrock overdeepenings exposed by continued glacial retreat can store precipitation and meltwater, potentially leading to the formation of new proglacial lakes. These lakes may pose threats of glacial lake outburst floods (GLOFs) in high mountain areas, particularly if new lakes form in geomorphological setups prone to triggering events such as landslides or moraine collapses. We present the first complete inventory for future glacial lakes in High Mountain Asia by computing the subglacial bedrock for ~100 000 glaciers and estimating overdeepening area, volume and impact hazard for the larger potential lakes. We detect 25 285 overdeepenings larger than 104 m2 with a volume of 99.1 ± 28.6 km3 covering an area of 2683 ± 773.8 km2. For the 2700 overdeepenings larger than 105 m2, we assess the lake predisposition for mass-movement impacts that could trigger a GLOF by estimating the hazard of material detaching from surrounding slopes. Our findings indicate a shift in lake area, volume and GLOF hazard from the southwestern Himalayan region toward the Karakoram. The results of this study can be used for anticipating emerging threats and potentials connected to glacial lakes and as a basis for further studies at suspected GLOF hazard hotspots.

Southern Ocean sea ice can exhibit widespread flooding and subsequent snow-ice formation, due to relatively thick snow covers compared to the total ice thickness. Considerable subkilometer scale variability in snow and ice thickness causes poorly constrained uncertainties in determining the amount of flooding that occurs. Using datasets of snow depth and ice thickness acquired in the Weddell Sea during austral winter 2013 (AWECS campaign) from three floes, we demonstrate large spatial variability of a factor 10 and 5 for snow and combined snow and ice thickness, respectively. The temporal evolution after the floe visit was recorded by automatic weather station and ice mass balance buoys. Using a physics-based, multi-layer snow/sea ice model in a one-dimensional and distributed mode to simulate the thermodynamic processes, we show that the distributed simulations, modeling flooding across the entire heterogeneous floe, produced vastly different amounts of flooding than one-dimensional single point simulations. Three times the flooding is produced in the one-dimensional simulation for the buoy location than distributed (floe-averaged) simulations. The latter is in close agreement with buoy observations. The results suggest that using point observations or one-dimensional simulations to extrapolate processes on the floe-scale can overestimate the amount of flooding and snow-ice formation.

Accurate knowledge of the ice thickness distribution and glacier bed topography is essential for predicting dynamic glacier changes and the future developments of downstream hydrology, which are impacting the energy sector, tourism industry and natural hazard management. Using AIR-ETH, a new helicopter-borne ground-penetrating radar (GPR) platform, we measured the ice thickness of all large and most medium-sized glaciers in the Swiss Alps during the years 2016–20. Most of these had either never or only partially been surveyed before. With this new dataset, 251 glaciers – making up 81% of the glacierized area – are now covered by GPR surveys. For obtaining a comprehensive estimate of the overall glacier ice volume, ice thickness distribution and glacier bed topography, we combined this large amount of data with two independent modeling algorithms. This resulted in new maps of the glacier bed topography with unprecedented accuracy. The total glacier volume in the Swiss Alps was determined to be 58.7 ± 2.5 km3 in the year 2016. By projecting these results based on mass-balance data, we estimated a total ice volume of 52.9 ± 2.7 km3 for the year 2020. Data and modeling results are accessible in the form of the SwissGlacierThickness-R2020 data package.

The validity of any glaciological paleo proxy used to interpret climate records is based on the level of understanding of their transfer from the atmosphere into the ice sheet and their recording in the snowpack. Large spatial noise in snow properties is observed, as the wind constantly redistributes the deposited snow at the surface routed by the local topography. To increase the signal-to-noise ratio and getting a representative estimate of snow properties with respect to the high spatial variability, a large number of snow profiles is needed. However, the classical way of obtaining profiles via snow-pits is time and energy-consuming, and thus unfavourable for large surface sampling programs. In response, we present a dual-tube technique to sample the upper metre of the snowpack at a variable depth resolution with high efficiency. The developed device is robust and avoids contact with the samples by exhibiting two tubes attached alongside each other in order to (1) contain the snow core sample and (2) to access the bottom of the sample, respectively. We demonstrate the performance of the technique through two case studies in East Antarctica where we analysed the variability of water isotopes at a 100 m and 5 km spatial scales.

Ice-sheet mass-balance estimates derived from repeat satellite-altimetry observations require accurate calculation of spatiotemporal variability in firn-air content (FAC). However, firn-compaction models remain a large source of uncertainty within mass-balance estimates. In this study, we investigate one process that is neglected in FAC estimates derived from firn-compaction models: enhanced layer thinning due to horizontal divergence. We incorporate a layer-thinning scheme into the Community Firn Model. At every time step, firn layers first densify according to a firn-compaction model and then thin further due to an imposed horizontal divergence rate without additional density changes. We find that horizontal divergence on Thwaites (THW) and Pine Island Glaciers can reduce local FAC by up to 41% and 18%, respectively. We also assess the impact of temporal variability of horizontal divergence on FAC. We find a 15% decrease in FAC between 2007 and 2016 due to horizontal divergence at a location that is characteristic of lower THW. This decrease accounts for 16% of the observed surface lowering, whereas climate variability alone causes negligible changes in FAC at this location. Omitting transient horizontal divergence in estimates of FAC leads to an overestimation of ice loss via satellite-altimetry methods in regions of dynamic ice flow.

As a result of climate warming, glacial meltwater discharge has been increasing in Greenland. During the summers of 2015 and 2016, there were rapid increases in discharge from Qaanaaq Glacier in northwestern Greenland. These discharges resulted in floods that destroyed the road linking the settlement of Qaanaaq to Qaanaaq Airport. Field measurements were performed and a numerical model of glacier runoff was developed to quantify these discharges. The high discharge associated with the 2015 flood, estimated at 9.1 m3 s−1 (hourly mean), resulted from intensive glacier melting due to warm air temperature and strong winds, while the high discharge associated with the 2016 flood resulted from heavy rainfall (90 mm d−1) that led to a peak discharge estimated at 19.9 m3 s−1. The developed model, when used to investigate future glacier runoff under warming conditions, revealed a nonlinear increase in glacial melt with increasing temperature. Additionally, the model forecasted a threefold increase in total summer discharge, owing to a 4 °C rise in temperature. Thus, this study quantified the impact of a changing climate on glacier runoff, which gives insight into future risks of flood hazards along the coast of Greenland.

Knowledge of glacier volume is crucial for ice flow modelling and predicting the impacts of climate change on glaciers. Rugged terrain, harsh weather conditions and logistic costs limit field-based ice thickness observations in the Himalaya. Remote-sensing applications, together with mathematical models, provide alternative techniques for glacier ice thickness and volume estimation. The objective of the present research is to assess the application of artificial neural network (ANN) modelling coupled with remote-sensing techniques to estimate ice thickness on individual glaciers with direct field measurements. We have developed two ANN models and estimated the ice thickness of Chhota Shigri Glacier (western Himalaya) on ten transverse cross sections and two longitudinal sections. The ANN model estimates agree well with ice thickness measurements from a ground-penetrating radar, available for five transverse cross sections on Chhota Shigri Glacier. The overall root mean square errors of the two ANN models are 24 and 13 m and the mean bias errors are ±13 and ±6 m, respectively, which are significantly lower than for other available models. The estimated mean ice thickness and volume for Chhota Shigri Glacier are 109 ± 17 m and 1.69 ± 0.26 km3, respectively.

In recent years, researchers have focused on the applications of uncrewed aerial vehicles (UAVs) in environmental remote sensing tasks. However, studies on glacier monitoring using UAV technology are relatively scarce, especially for high mountain glacier monitoring. To explore the feasibility of UAV technology for high mountain glaciers, four UAV surveys were deployed on two glaciers of the central Tibetan Plateau. Based on the images retrieved by UAV in 2017 and 2019, orthomosaics and digital elevation models were produced to quantify the length, area and elevation changes in the ablation zone of these two glaciers at different times. Additionally, we utilized several Landsat scenes to derive glacier changes over the last 30 years and combined these with the UAV data to assess the advantages and disadvantages of UAV technology in mountain glacier monitoring.

In this study, we use aerial photographs, satellite imagery and field observations to quantify changes in the area, terminus length, snowline elevation and surface elevation of eight glaciers in the Alexandra Fiord region, eastern Ellesmere Island, between 1959 and 2019. Comparisons to written and pictorial descriptions from the British Arctic Expedition extend the record of change in terminus position and surface elevation to 1875 for Twin Glacier. Glacier area at Alexandra Fiord decreased by a total of 15.77 ± 0.65 km2 (11.77 ± 0.49%) between 1959 and 2019, the mean end of summer snowline increased in elevation by 360 ± 84 m (8 ± 2 m a−1) between 1974 and 2019, and the glaciers thinned at an average rate of 0.60 ± 0.06 m a−1 between 2001 and 2018. Annual rates of terminus retreat were ~3–5 times higher over the period 1974–2019 compared to 1875–1974, and rates of thinning were ~2–3 times higher over 2001–18 compared to 1875–2001. Our results are consistent with rates of change determined for other glaciers of similar size on Ellesmere Island, and with accelerated rates of ice loss coincident with regional increases in air temperature of ~1.5°C since the early 1980s.

In this study, we explore small-scale (~1 to 20 km) thermal-refractive effects on basal geothermal heat flux (BGHF) at subglacial boundaries resulting from lateral thermal conductivity contrasts associated with subglacial topography and geologic contacts. We construct a series of two-dimensional, conductive, steady-state models that exclude many of the complexities of ice sheets in order to demonstrate the effect of thermal refraction. We show that heat can preferentially flow into or around a subglacial valley depending on the thermal conductivity contrast with underlying bedrock, with anomalies of local BGHF at the ice–bedrock interface between 80 and 120% of regional BGHF and temperature anomalies on the order of ±15% for the typical range of bedrock conductivities. In the absence of bed topography, subglacial contacts can produce significant heat flux and temperature anomalies that are locally extensive (>10 km). Thermal refraction can result in either an increase or decrease in the likelihood of melting and ice-sheet stability depending on the conductivity contrast and bed topography. While our models exclude many of the physical complexities of ice behavior, they illustrate the need to include refractive effects created by realistic geology into future glacial models to improve the prediction of subglacial melting and ice viscosity.

The Tibetan Plateau (TP) lakes are sensitive to climate change due to its seasonal ice cover, but few studies have paid attention to the freeze-thaw process of TP lakes and its key control parameters. By combining 216 simulation experiments using the LAKE2.0 model with the observations, we evaluated the effects of ice and snow albedo, ice (Kdi) and water (Kdw) extinction coefficients on the lake ice phenology, water temperature, sensible and latent heat fluxes. The reference experiment performs well in simulating the lake temperature, with a small positive bias increasing with depth, but it underestimates the ice thickness. The increase of ice albedo, snow albedo and Kdi induce a significant decrease in water temperature. Compared with the latent heat, the sensible heat flux is more sensitive to these three parameters. The ice thickness increases almost linearly with the increase of ice albedo but decreases with the increase of Kdi. The ice thickness and frozen days vary little with Kdw, but increasing Kdw can decrease the water temperature. Compared with the ice albedo, the Kdi and snow albedo have a large effect on the number of frozen days. This study brings to light the necessity to improve the parameterizations of the TP lakes freeze-thaw process.

Global warming is causing glaciers in the Caucasus Mountains and around the world to lose mass at an accelerated pace. As a result of this rapid retreat, significant parts of the glacierized surface area can be covered with debris deposits, often making them indistinguishable from the surrounding land surface by optical remote-sensing systems. Here, we present the DebCovG-carto toolbox to delineate debris-covered and debris-free glacier surfaces from non-glacierized regions. The algorithm uses synthetic aperture radar-derived coherence images and the normalized difference snow index applied to optical satellite data. Validating the remotely-sensed boundaries of Ushba and Chalaati glaciers using field GPS data demonstrates that the use of pairs of Sentinel-1 images (2019) from identical ascending and descending orbits can substantially improve debris-covered glacier surface detection. The DebCovG-carto toolbox leverages multiple orbits to automate the mapping of debris-covered glacier surfaces. This new automatic method offers the possibility of quickly correcting glacier mapping errors caused by the presence of debris and makes automatic mapping of glacierized surfaces considerably faster than the use of other subjective methods.

Glacier surficial melt rates are commonly modelled using surface energy balance (SEB) models, with outputs applied to extend point-based mass-balance measurements to regional scales, assess water resource availability, examine supraglacial hydrology and to investigate the relationship between surface melt and ice dynamics. We present an improved SEB model that addresses the primary limitations of existing models by: (1) deriving high-resolution (30 m) surface albedo from Landsat 8 imagery, (2) calculating shadows cast onto the glacier surface by high-relief topography to model incident shortwave radiation, (3) developing an algorithm to map debris sufficiently thick to insulate the glacier surface and (4) presenting a formulation of the SEB model coupled to a subsurface heat conduction model. We drive the model with 6 years of in situ meteorological data from Kaskawulsh Glacier and Nàłùdäy (Lowell) Glacier in the St. Elias Mountains, Yukon, Canada, and validate outputs against in situ measurements. Modelled seasonal melt agrees with observations within 9% across a range of elevations on both glaciers in years with high-quality in situ observations. We recommend applying the model to investigate the impacts of surface melt for individual glaciers when sufficient input data are available.

To get a better overview of atmosphere-driven mass changes at the 79N Glacier (Nioghalvfjerdsfjorden Glacier), the largest outlet glacier of the northeast Greenland ice stream, the surface mass balance (SMB) is modeled by linking the COupled Snowpack and Ice surface energy and mass-balance model in PYthon (COSIPY) with the output of a regional atmospheric model (Polar WRF) for the years 2014–2018. After a manual model optimization, the model produces reliable results when compared to observations in the region and to values from the literature. High spatial resolution (1 km) simulations reveal strong interannual variability of the SMB. Stronger surface melting increased the ablation and runoff in years with high mass loss (2016 and 2017) whereas in other years (2015 and 2018) melting and refreezing inside the snowpack dominated the mass balance (MB). A cooler regional climate with higher snowfall-driven accumulation, higher albedo and reduced surface melt in the ablation period of 2018 resulted in a positive SMB in 2018, however, the annual total MB remained negative. The results suggest a promising new dataset for gaining more insights into mass-balance processes and their contribution to the acceleration of glacier retreat in northeast Greenland.

Stable isotope ratios (δ18O and δD) in Antarctic snow and ice are basic proxy indices of climate in ice core studies. The relation between the ratios has important indicative significance for moisture sources. In general, the fractionation characteristics of the two isotopes vary with different meteorological and topographical conditions. This paper presents the spatial and temporal distribution of meteoric water line (MWL) slopes along a traverse from the Zhongshan Station (ZSS) to Dome A in East Antarctica. It is found that the slopes decrease with the increasing distance inland from the coast and the lowest slope occurred at Dome A, where the long-range transported moisture dominates and clear sky snowing have an influence. The slopes in different layers of the snowpack showed a decreasing trend with depth and this is attributed to the fractionation during the interstitial sublimation and re-condensation processes of the water vapor. Frost flower development on the interior plateau surface can greatly alter the depth evolution of the MWL slope. The coastal snow pits also go through the post-depositional smoothing effect, but their influences are not so significant as the inland regions.

Photogrammetric reconstructions of the Aldegondabreen glacier on Svalbard from 17 archival terrestrial oblique photographs taken in 1910 and 1911 reveal a past volume of 1373.7 ± 78.2 · 106 m3; almost five times greater than its volume in 2016. Comparisons to elevation data obtained from aerial and satellite imagery indicate a relatively unchanging volume loss rate of − 10.1 ± 1.6 · 106 m3 a−1 over the entire study period, while the rate of elevation change is increasing. At this rate of volume loss, the glacier may be almost non-existent within 30 years. If the changes of Aldegondabreen are regionally representative, it suggests that there was considerable ice loss over the entire 1900s for the low elevation glaciers of western Svalbard. The 1910/11 reconstruction was made from a few of the tens of thousands of archival terrestrial photographs from the early 1900s that cover most of Svalbard. Further analysis of this material would give insight into the recent history and future prospects of the archipelago's glaciers. Photogrammetric reconstructions of this kind of material require extensive manual processing to produce good results; for more extensive use of these archival imagery, a better processing workflow would be required.

The Kaskawulsh Glacier is an iconic outlet draining the icefields of the St. Elias Mountains in Yukon, Canada. We determine and attempt to interpret its catchment-wide mass budget since 2007. Using SPOT5/6/7 data we estimate a 2007–18 geodetic balance of −0.46 ± 0.17 m w.e. a−1. We then compute balance fluxes and observed ice fluxes at nine flux gates to examine the discrepancy between the climatic mass balance and internal mass redistribution by glacier flow. Balance fluxes are computed using a fully distributed mass-balance model driven by downscaled and bias-corrected climate-reanalysis data. Observed fluxes are calculated using NASA ITS_LIVE surface velocities and glacier cross-sectional areas derived from ice-penetrating radar data. We find the glacier is still in the early stages of dynamic adjustment to its mass imbalance. We estimate a committed terminus retreat of ~23 km under the 2007–18 climate and a lower bound of 46 km3 of committed ice loss, equivalent to ~15% of the total glacier volume.

The thickness of a supraglacial layer is critical to the magnitude and time frame of glacier melt. Field-based, short pulse, ground-penetrating radar (GPR) has successfully measured debris thickness during a glacier's melt season, when there is a strong return from the ice–debris interface, but profiling with GPR in the absence of a highly reflective ice interface has not been explored. We investigated the performance of 960 MHz signals over 2 km of transects on Changri Nup Glacier, Nepal, during the post-monsoon. We also performed laboratory experiments to interpret the field data and investigate electromagnetic wave propagation into dry rocky debris. Laboratory tests confirmed wave penetration into the glacier ice and suggest that the ice–debris interface return was missing in field data because of a weak dielectric contrast between solid ice and porous dry debris. We developed a new method to estimate debris thicknesses by applying a statistical approach to volumetric backscatter, and our backscatter-based calculated thickness retrievals gave reasonable agreement with debris depths measured manually in the field (10–40 cm). We conclude that, when melt season profiling is not an option, a remote system near 1 GHz could allow dry debris thickness to be estimated based on volumetric backscatter.

Glacier mass-balance observations at seasonal resolution have been performed since 1914 at two sites on Claridenfirn, Switzerland. The measurements are the longest uninterrupted records of glacier mass balance worldwide. Here, we provide a complete re-analysis of the 106-year series (1914–2020), focusing on both point and glacier-wide mass balance. The approaches to evaluate and homogenize the direct observations are described in detail. Based on conservative assumptions, average uncertainties of $\pm$0.25 m w.e. are estimated for glacier-wide mass balances at the annual scale. It is demonstrated that long-term variations in mass balance are clearly driven by melting, whereas decadal changes in accumulation are uncorrelated with mass balance and can only be relevant in short periods. Mass change of Claridenfirn is impacted by dry calving at a frontal ice cliff. Considerations of ice volume flux at a cross-profile reveal long-term variations in frontal ice loss accounting for $\sim$9% of total annual ablation on average. The effect of changes in frontal ablation mostly explains $\lt$10% of the mass-balance difference relative to the period 1960–1990, but accounts for $\sim$20% in 2010–2020. Glacier mass changes are discussed in the context of observations throughout the European Alps indicating that Claridenfirn is regionally representative.

We develop a two-dimensional, plan-view formulation of ice-shelf flow and viscoelastic ice-shelf flexure. This formulation combines, for the first time, the shallow-shelf approximation for horizontal ice-shelf flow (and shallow-stream approximation for flow on lubricated beds such as where ice rises and rumples form), with the treatment of a thin-plate flexure. We demonstrate the treatment by performing two finite-element simulations: one of the relict pedestalled lake features that exist on some debris-covered ice shelves due to strong heterogeneity in surface ablation, and the other of ice rumpling in the grounding zone of an ice rise. The proposed treatment opens new venues to investigate physical processes that require coupling between the longitudinal deformation and vertical flexure, for instance, the effects of surface melting and supraglacial lakes on ice shelves, interactions with the sea swell, and many others.