Climate change is causing Himalayan glaciers to shrink rapidly and natural hazards to increase, while downstream exposure is growing. Glacier shrinkage promotes the formation of glacial lakes, which can suddenly drain and produce glacier lake outburst floods (GLOFs). Bhutan is one of the most vulnerable countries globally to these hazards. Here we use remotely sensed imagery to quantify changes in supraglacial water storage on Tshojo Glacier, Bhutan, where previous supraglacial pond drainage events have necessitated downstream evacuation. Results showed a doubling of both total ponded area (104 529 m2 to 213 943 m2) and its std dev. (64 808 m2 to 158 550 m2) between the periods 1987–2003 and 2007–2020, which was predominantly driven by increases in the areas of the biggest ponds. These ponds drained regularly and have occupied the same location since at least 1967. Tshojo Glacier has remained in the first stage of proglacial lake development for 53 years, which we attribute to its moderate slopes and ice velocities. Numerical modelling shows that pond outbursts can reach between ~6 and 47 km downstream, impacting the remote settlement of Lunana. Our results highlight the need to better quantify variability in supraglacial water storage and its potential to generate GLOFs, as climate warms.

Ridge B is one of the least studied areas in Antarctica but has been considered to be a potential location for the oldest ice on Earth. Among important parameters for calculating where very old ice may exist, geothermal heat flux (GHF) is critical but poorly understood. Here, GHF is determined by quantifying the transitions between dry and wet basal conditions using a radioglaciological method applied to airborne radio-echo sounding data. GHF is then constrained by a thermodynamic model matched to the transitions. The results show that GHF in Ridge B varies locally and ranges from 48.5 to 65.1 mW m−2, with an average value of 58.0 mW m−2, which is consistent with the current known GHF constrained by subglacial lakes and derived from Vostok ice core temperature measurements. Our work highlights the value of considering local GHF when locating the oldest ice in this potential region or other regions.

At high elevations on the Greenland ice sheet meltwater percolates and refreezes in place, and hence does not contribute to mass loss. However, meltwater generation and associated surface runoff is occurring from increasingly higher altitudes, causing changes in firn stratigraphy that have led to the presence of near-surface ice slabs. These ice slabs force meltwater to flow laterally instead of percolating downwards. Here we present a simple, physics-based quasi-2-D model to simulate lateral meltwater runoff and superimposed ice (SI) formation on top of ice slabs. Using an Eulerian Darcy flow scheme, the model calculates how far meltwater can travel within a melt season and when it appears at the snow surface. Results show that lateral flow is a highly efficient runoff mechanism, as lateral outflow exceeds locally generated meltwater in all model gridcells, with total meltwater discharge sometimes reaching more than 30 times the average amount of in situ generated melt. SI formation, an important process in the formation and thickening of the ice slabs, can retain up to 40% of the available meltwater, and generally delays the appearance of visible runoff. Validating the model against field- or remote-sensing data remains challenging, but the results presented here are a first step towards a more comprehensive understanding and description of the hydrological system in the accumulation zone of the southwestern Greenland ice sheet.

Runoff from heavily glacierised catchments, its seasonality and the contribution from different storage units are highly relevant for assessing the seasonal and long-term water supply and flood prediction. Modelling studies on runoff from such basins often use simulated meteorological input (e.g. downscaling products) based on remote observations. We investigate the contribution of snow, firn and ice to runoff in the Vernagtferner basin (Ötztal Alps) from 2020 to 2022, using a physical modelling chain driven by a local observation network. We use the SNOWPACK model to simulate snow/ice development at observation sites and the Alpine3D model to calculate accumulation and melt at the catchment scale. Basin discharge is estimated using a gridded version of the HBV-ETH model. This approach largely reproduces observed glacier mass balances, while modelled and measured basin discharge are in good agreement. Snowmelt dominates discharge in the early melt season, while ice melt becomes increasingly important during summer. This is in strong contrast to the near-equilibrium mass balances in the 1980s, when ice melt played a minor role for annual discharge. The strong reduction of the accumulation area leads to a fundamental change from a snowmelt regime to an ice melt regime, which is especially pronounced in 2022.

Sub-glacial canyon features up to 580 m deep between flat terraces were identified beneath Devon Ice Cap during a 2023 radar echo sounding (RES) survey. The largest canyon connects a hypothesized brine network near the Devon Ice Cap summit with the marine-terminating Sverdrup outlet glacier. This canyon represents a probable drainage route for the hypothesized water system. Radar bed reflectivity is consistently 30 dB lower along the canyon floor than on the terraces, contradicting the signature expected for sub-glacial water. We compare these data with backscattering simulations to demonstrate that the reflectivity pattern may be topographically induced. Our simulated results indicated a 10 m wide canal-like water feature is unlikely along the canyon floor, but smaller features may be difficult to detect via RES. We calculated basal temperature profiles using a 2D finite difference method and found the floor may be up to 18°C warmer than the terraces. However, temperatures remain below the pressure melting point, and there is limited evidence that the canyon floor supports a connected drainage system between the DIC summit and Sverdrup Glacier. The terrain beneath Devon Ice Cap demonstrates limitations for RES. Future studies should evaluate additional correction methods near complex terrain, such as RES simulation as we demonstrate here.

Observations of glacier melt and runoff are of fundamental interest in the study of glaciers and their interactions with their environment. Considerable recent interest has developed around distributed acoustic sensing (DAS), a sensing technique which utilizes Rayleigh backscatter in fiber optic cables to measure the seismo-acoustic wavefield in high spatial and temporal resolution. Here, we present data from a month-long, 9 km DAS deployment extending through the ablation and accumulation zones on Rhonegletscher, Switzerland, during the 2020 melt season. While testing several types of machine learning (ML) models, we establish a regression problem, using the DAS data as the dependent variable, to infer the glacier discharge observed at a proglacial stream gauge. We also compare two predictive models that only depend on meteorological station data. We find that the seismo-acoustic wavefield recorded by DAS can be utilized to infer proglacial discharge. Models using DAS data outperform the two models trained on meteorological data with mean absolute errors of 0.64, 2.25 and 2.72 m3 s−1, respectively. This study demonstrates the ability of in situ glacier DAS to be used for quantifying proglacial discharge and points the way to a new approach to measuring glacier runoff.

On the Greenland Ice Sheet, hydrofracture connects the supraglacial and subglacial hydrologic systems, coupling surface runoff dynamics and ice velocity. In recent decades, the growth of low-permeability ice slabs in the wet snow zone has expanded Greenland's runoff zone, but observations suggest that surface-to-bed connections are rare, because meltwater drains through crevasses into the porous firn beneath ice slabs. However, there is little quantitative evidence confirming the absence of surface-to-bed fracture propagation. Here, we use poromechanics to investigate whether water-filled crevasses in ice slabs can propagate vertically through an underlying porous firn layer. Based on numerical simulations, we develop an analytical estimate of the water injection-induced effective stress in the firn given the water level in the crevasse, ice slab thickness, and firn properties. We find that the firn layer substantially reduces the system's vulnerability to hydrofracture because much of the hydrostatic stress is accommodated by a change in pore pressure, rather than being transmitted to the solid skeleton. This result suggests that surface-to-bed hydrofracture will not occur in ice slab regions until all pore space proximal to the initial flaw has been filled with solid ice.

Global Navigation Satellite System (GNSS) observations and ground-based timelapse photography obtained over the record-high 2019/2020 melt season are combined to characterise the flexure and fracture behaviour of a previously formed doline on George VI Ice Shelf, Antarctica. The GNSS timeseries shows a downward vertical displacement of the doline centre with respect to the doline rim of ~60 cm in response to loading from a central meltwater lake. The GNSS data also show a tens-of-days episode of rapid-onset, exponentially decaying horizontal displacement, where the horizontal distance between the doline rim and its centre increases by ~70 cm. We interpret this event as the initiation and/or widening of a fracture, aided by stress perturbations associated with meltwater loading in the doline basin. Viscous flexure modelling indicates that the meltwater loading generates tensile surface stresses exceeding 75 kPa. This, together with our timelapse photos of circular fractures around the doline, suggests the first such documentation of meltwater-loading-induced ‘ring fracture’ formation on an ice shelf, equivalent to the fracture type proposed as part of the chain-reaction lake drainage process involved in the 2002 breakup of the Larsen B Ice Shelf.

Glacier surges are opportunities to study large amplitude changes in ice velocities and accompanying links to subglacial hydrology. Although the surge phase is generally explained as a disruption in the glacier's ability to drain water from the bed, the extent and duration of this disruption remain difficult to observe. Here we present a combination of in situ and remotely sensed observations of subglacial water discharge and evacuation during the latter half of an active surge and subsequent quiescent period. Our data reveal intermittently efficient subglacial drainage prior to surge termination, showing that glacier surges can persist in the presence of channel-like subglacial drainage and that successive changes in subglacial drainage efficiency can modulate active phase ice dynamics at timescales shorter than the surge cycle. Our observations favor an explanation of fast ice flow sustained through an out-of-equilibrium drainage system and a basal water surplus rather than binary switching between states in drainage efficiency.

Surface height changes above three previously undetected subglacial lakes in northeastern Greenland are documented using CryoSat, DEMs and ICESat-2. Between 7 February and 6 March 2012, the central ice region (22.6 km2) above the largest lake dropped by ~37 m followed by a further drop of 12 m in the following 29 days. This implies a subglacial water outflow, or jökulhlaup, of at least 1 km3 at rates of hundreds of cubic meters per second. A comparable outflow occurred again between 23 July and 15 September 2019, with smaller outflows in the fall of 2014 and 2016. In contrast, a second smaller subglacial lake at a higher elevation had two subglacial outbursts of ~0.3 km3 in 2012 and 2019 but the lake filling was gradual and not strongly seasonal or episodic. Water remained in both lakes after the outflows but this may not be the case for the third smallest and lowest subglacial lake. While there appears to be some hydrological link between the three lakes, the flux of water moving under the ice in this area appears to be larger than would be expected from local summer melt. However, the source of the excess water remains uncertain.

The number of glacial lakes has grown globally concurrently with the retreat of glaciers in the last few decades, increasing the risk of potentially hazardous glacial lake outburst floods (GLOFs) and posing a threat to downstream communities. Norway has several known ice-dammed lakes that produce repeated GLOFs but as we show here, the existing GLOF database is incomplete and needs to be improved through continuous monitoring of glaciers and glacial lakes. This study examines the case of an ice cap in central Norway hosting at least four drainage-prone lakes. We reconstruct the sequence of lake drainage events through a combination of remote sensing, ground-truthing and citizen science while scrutinising the applicability of the PlanetScope imagery vs Sentinel-2 and Landsat-8 OLI products. As opposed to the Landsat imagery that often fails to resolve even the largest glacial lakes of Tystigbreen, both PlanetScope and Sentinel-2 are helpful in identifying previously unrecognised glacial lakes and undocumented drainage events. Our analysis suggests that a fusion of the two satellite products may be beneficial for automated tracking of glacial lake changes. We also demonstrate that local knowledge and systematic involvement of citizens in data collection have a potential to enrich GLOF databases.

Subglacial hydrology models struggle to reproduce seasonal drainage patterns that are consistent with observed subglacial water pressures and surface velocities. We modify the standard sheet-flow parameterization within a coupled sheet–channel subglacial drainage model to smoothly transition between laminar and turbulent flow based on the locally computed Reynolds number in a physically consistent way (the ‘transition’ model). We compare the transition model to standard laminar and turbulent models to assess the role of the sheet-flow parameterization in reconciling observed and modelled water pressures under idealized and realistic forcing. Relative to the turbulent model, the laminar and transition models improve seasonal simulations by increasing winter water pressure and producing a more prominent late-summer water pressure minimum. In contrast to the laminar model, the transition model remains consistent with its own internal assumptions across all flow regimes. Based on the internal consistency of the transition model and its improved performance relative to the standard turbulent model, we recommend its use for transient simulations of subglacial drainage.

To better constrain meltwater transport and ice viscosity in temperate glaciers, particularly in ice stream shear margins, we use a custom permeameter to study the untested model relationship between the permeability of temperate ice and its liquid water content. The permeability of lab-made ice of two mean grain diameters (1.8 and 4.2 mm) is measured, and water content is controlled with the ice salinity and measured calorimetrically. Fluorescein dye is added to through-flowing, chilled water to highlight flow pathways through the ice after experiments. As predicted by a simple model, permeability increases with approximately the square of the water content and by about three orders of magnitude across water contents of 0.1–4.4%. However, permeability values are less than those of the model by average factors of 2.6 and 4.1 for the finer and coarser ice, respectively. This discrepancy is likely due to tortuous, truncated or air-clogged veins. The order-of-magnitude agreement between measured and modeled values may indicate that reduced permeability from these factors is nearly compensated by preferential flow in oversized veins that are isolated or arborescent. Both kinds of preferred flow pathways are observed but the latter only in fine-grained ice at water contents > 2%.

Active subglacial lakes provide a rare glimpse of the subglacial environment and hydrological processes at play. Several studies contributed to establishing active subglacial lake inventories and document lake drainage and connection, but few focused on the period between lake drainage when the melt production and transport contribute to the refilling of these lakes. In this study, we employ high-resolution CryoSat-2 altimetry data from 2010 to 2021 to compile an inventory of recharging lakes across Antarctica. We extract recharge rates from these lakes, which serve as a lower limit on subglacial melt production. These recharge rates are compared against predictions obtained by routing modelled subglacial meltwater at the ice-sheet's base. Our findings indicate that modelled recharge rates are consistent with observations in all but one of the investigated lakes, providing a lower bound on geothermal heat fluxes. Lake Cook E2 displays recharge rates far exceeding predictions, indicating that processes are taking place that are currently unaccounted for. Considering recharge in hydrologically connected lake networks instead of individually provides a stricter constraint on melt production. Recharge rates extracted from the Thwaites Lake system suggest that subglacial melt production has been underestimated.

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.

A comparison of major ion chemistry of subglacial boreholes and discharging subglacial waters reveals three fundamentally different glacier hydrochemical regimes. Subglacial waters from alpine glaciers have chemistry distinct from the subglacial waters of Greenland or Antarctica. Greenland and Antarctica also differ fundamentally from each other, with Greenland Ice Sheet waters, at least during the summer melt season, remaining dilute and unaffected by saturation reactions and Antarctic Ice Sheet waters controlled by a range of saturation states. Some Antarctic waters form concentrated brines, capable of depressing the freezing point by >10°C. While these waters have only been directly sampled where they rarely emerge, geophysical observations from Devon Ice Cap and Greenland show liquid water at the glacier bed in locations where ice is thin and slowly moving and a cold bed is otherwise expected. This raises the possibility that lithogenic subglacial brines could be widespread and that our existing subglacial hydrochemical measurements might be biased by seasonal sampling of freely discharging water. The potential for diverse ranges of subglacial environments under ice sheets suggests the need for new and ambitious sampling programs to characterize difficult to access subglacial waters and quantify their impact on glacier dynamics, geobiology and global geochemical cycling.

Although over 600 Antarctic subglacial lakes have been identified using radar and satellite observations, the bathymetry and bed properties, which are key to understanding conditions within the lake, have been determined in very few localities. We present measurements of water column thickness and lakebed properties from Lago Subglacial CECs (SLC), located beneath 2653 m of ice at the Rutford-Institute-Minnesota divide in Antarctica. Seismic profiles indicate a maximum water column thickness of 301.3 ± 1.5 m, at the widest part of the lake, with an estimated lake volume of 2.5 ± 0.3 km3. Seismic imaging and measurements of the reflection strength at the ice base and lakebed indicate >15 m of high-porosity fine-grained sediment in the central section of the lakebed, consistent with a depositional sequence with an age of up to 0.5 Ma. These observations, along with previous radar measurements and modelling, indicate a low-energy sedimentary environment with a long water-residence time. As such, SLC is a suitable target for exploration via direct access to recover sediment records of ice sheet and climate history and investigate microbial life with long periods of isolation.

Water pressure beneath glaciers influences ice velocity. Subglacial hydrology models are helpful for gaining insight into basal conditions, but models depend on unconstrained parameters, and a current challenge is reproducing elevated water pressures in winter. We eliminate terms related to englacial storage, opening by sliding, and melt due to changes in the pressure-melting-point temperature, to create a minimalist version of the Subglacial Hydrology And Kinetic, Transient Interactions (SHAKTI) model, and apply this model to Helheim Glacier in east Greenland to explore the winter base state of the subglacial drainage system. Our results suggest that meltwater produced at the bed alone supports active winter drainage with large areas of elevated water pressure and preferential drainage pathways, using a continuum approach that allows for transitions between flow regimes. Transmissivity varies spatially over several orders of magnitude from 10−4 to 103 m2s−1, with regions of weak transmissivity representing poorly connected regions of the system. Bed topography controls the location of primary drainage pathways, and high basal melt rates occur along the steep valley walls. Frictional heat from sliding is a dominant source of basal melt; different approaches for calculating basal shear stress produce significantly different basal melt rates and subglacial discharge.

This study presents the first light detection and ranging (LiDAR)-based morphometric description of a water blister from a past ice-sheet bed caused by rapid supraglacial drainage. The blister formed during the rapid early Holocene deglaciation of the Fennoscandian Ice Sheet (FIS). It is located in southwest Finland within a subglacial meltwater route interpreted to represent the transition from a distributed to a channelized drainage system. A LiDAR digital elevation model was supplemented with sedimentological and ground-penetrating radar data on blister outflow channels and sedimentology of downflow polymorphous mounds and ridges (PMRs). Unlike the water blisters recorded from the rapid drainage of supraglacial lakes on the Greenland Ice Sheet, the smaller blister size here was either due to crevasse or moulin drainage, or was a supraglacial lake drainage that tapped into a pre-existing, relatively efficient drainage system and related semi-sorted sediments, promoting rapid drainage and reworking of PMRs along the meltwater route. The preservation potential or exposure probability of blister marks is presumably low but they can provide important information about evolution of subglacial drainage systems that is of value to modern interpretations of glacial hydrology.

The rapid formation of glacial lakes is one of the most conspicuous landscape changes caused by atmospheric warming in glacierised regions. However, relatively little is known about the history and current state of glacial lakes in the High Arctic. This study aims to address this issue by providing the first inventory of glacial lakes in Svalbard, focusing in particular on the post-Little Ice Age evolution of glacial lakes and their typology. To do so, we used aerial photographs and satellite imagery together with archival topographic data from 1936 to 2020. The inventory comprises the development of 566 glacial lakes (146 km2) that were still in direct contact with glaciers during the period 2008–2012. The results show a consistent increase in the total area of glacial lakes from the 1930s to 2020 and suggest an apparent link between climatic and geological factors, and the formation of specific lake dam types: moraine, ice, or bedrock. We also detected 134 glacial lake drainage events that have occurred since the 1930s. This study shows that Svalbard has one of the highest rates of glacial lake development in the world, which is an indicator of the overall dynamics of landscape change in the archipelago in response to climate change.