Seasonal glacier dynamics are key to predicting hazards and glacier stability due to short-term events as well as improving glacier models. However, the short-term ice velocity variations remain poorly constrained for slow-moving glaciers in the Himalaya due to the scarcity of in situ observations and limitations of satellite data and methods. We present seasonal velocity variations of Drang Drung Glacier (western Himalaya) in 2021 using Sentinel-1 phase-based Interferometric Synthetic Aperture Radar and offset tracking. Smoothed velocity estimates reveal $\sim$400 % seasonal variability (3–13 m a$^{-1}$), with speedups in spring and autumn and slowdowns in summer and winter. We relate these patterns to changes in radar backscatter and seasonal widening of the proglacial stream observed in Planet imagery. To interpret the mechanisms, we simulate the evolution with the Subglacial Hydrology and Kinetic Transient Interactions model coupled to the Ice-sheet and Sea-level System Model. Results indicate that speedup–slowdown cycles and their upglacier migration are driven by meltwater-induced shifts in subglacial drainage efficiency. This study emphasizes the role of hydrology and basal sliding in Himalayan glacier dynamics, often oversimplified in existing models.

Quarrying is a significant, locally dominant glacial erosion process. For settings where glaciers cut into partially intact bedrock, prior work has hypothesized that it occurs when glaciers impose spatially concentrated loads to drive fracture growth in the underlying rock, linking pre-existing fractures to complete dislodgment. This prior work, however, has not rigorously explained how most of this process occurs or whether it can leave the bed with a form susceptible to subsequent quarrying. We use a numerical model that combines finite element and discrete element capabilities to calculate the co-evolution of stress, elastic deformation, and fracturing in a granite and a weak sandstone containing discontinuous prior fractures. We find that quarrying is achievable in situations with rapid glacier sliding, as expected from prior work, but only if additional factors contribute. These include, especially, transient episodes when loading increasingly concentrates on the lips of bedrock steps, imposition of shear traction by friction between entrained clasts and the bed, and exploitation of anisotropic structural weaknesses in the bedrock. Hydraulic fracturing can significantly reduce the loads needed for quarrying if low hydraulic transmissivity allows for large water pressure differences between saturated fractures and the adjacent subglacial water system.

The difference between the ice and water pressures, or the effective pressure, influences water flow and sliding at the ice-bed interface. Effective pressure is typically quantified with subglacial hydrology models because direct measurements of the subglacial environment are sparse. Active subglacial lakes provide an opportunity to constrain effective pressures with altimetry because subglacial water-volume changes manifest at the ice-sheet surface as elevation-change anomalies. Here, we develop a method for estimating effective pressures from altimetry data above active subglacial lakes. We synthesise a previous theory of subglacial lake effective pressure with an altimetry-based inverse method that relates elevation-change data to water-volume changes. We apply the method to elevation-change data from NASA’s ICESat-2 satellite altimetry mission over several active lakes in Antarctica. We find that deviations from flotation (zero effective pressure) are typically a negligible fraction of the overburden (e.g., $\sim$10 kPa), although larger deviations can arise when the ice viscosity is large. For example, effective pressures over subglacial lake Byrds10 in East Antarctica locally reached magnitudes on the order of the tensile strength of glacier ice (e.g., over 100 kPa). These effective pressure estimates can constrain subglacial hydrology models in regions with active subglacial lakes and provide new insights into glacier-bed dynamics.

Subglacial drainage models, often motivated by the relationship between hydrology and ice flow, sensitively depend on numerous unconstrained parameters. We explore using borehole water-pressure time series to calibrate the uncertain parameters of a popular subglacial drainage model, taking a Bayesian perspective to quantify the uncertainty in parameter estimates and in the calibrated model predictions. To reduce the computation time associated with Markov Chain Monte Carlo sampling, we construct a fast Gaussian process emulator to stand in for the subglacial drainage model. We first carry out a calibration experiment using synthetic observations consisting of model simulations with hidden parameter values as a demonstration of the method. Using real borehole water pressures measured in western Greenland, we find meaningful constraints on four of the eight model parameters and a factor-of-three reduction in uncertainty of the calibrated model predictions. These experiments illustrate Gaussian process-based Bayesian inference as a useful tool for calibration and uncertainty quantification of complex glaciological models using field data. However, significant differences between the calibrated model and the borehole data suggest that structural limitations of the model, rather than poorly constrained parameters or computational cost, remain the most important constraint on subglacial drainage modelling.

We analyse seismic time series collected during experimental campaigns in the area of the David Glacier, Victoria Land, Antarctica, between 2003 and 2016. We observe hundreds of repeating seismic events, characterized by highly correlated waveforms (cross-correlation > 0.95), which mainly occur in the grounding zone, i.e. the region where the ice transitions from grounded ice sheet to freely floating ice shelf. The joint analysis of seismic events and observed local tidal measurements suggests that seismicity is not only triggered by a regular, periodic driver such as the ocean tides but also more likely by transient pulses. We consider potential environmental processes and their impact on the coupling between the glacier flow and the bedrock brittle failure. Among the environmental variables examined, our findings suggest that clustered and repeated seismic events may be related to transient episodes of ice-mass discharge correlated to a change in the subglacial hydrographic system that originates upstream of the glacier, lubricating the interface with the bedrock. This hypothesis is supported by the gravity variation observations provided by the GRACE satellite mission, which observed mass variations during periods characterized by seismic clustering.

Glacier collapse features, linked to subglacial cavities, are increasingly common on retreating Alpine glaciers. These features are hypothesized to result from glacier downwasting and subsurface ablation processes but the understanding regarding their distribution, formation and contribution to glacier mass loss remains limited. We present a Swiss-wide inventory of 223 collapse features observed over the past 50 years, revealing a sharp increase in their occurrence since the early 2000s. Using high-resolution digital elevation models, we derive a relationship between collapse feature area and ice ablation and estimate the Swiss-wide contribution of collapse features to glacier mass loss to be $19.8\times 10^6\,\text{m}^3$ of ice between 1971 and 2023. Based on extensive observations at Rhonegletscher, including surface displacement, ground-penetrating radar and drone-based elevation models, we quantify subsurface ablation rates of up to 27 cm d−1 and provide a detailed description of the collapse processes. We propose that glacier downwasting, enhanced energy supply through subglacial conduits and locally increased basal melt are key components to subglacial cavity growth. Our results highlight the importance of collapse features in the ongoing retreat of Alpine glaciers, stressing the need for further research to understand their formation and long-term implications for glacier dynamics under climate change.

Fast glacier motion is facilitated by slip at the ice-bed interface. For slip over rigid beds, areas of ice-bed separation (cavities) can exert significant control on slip dynamics. Analytic models of these systems assume that cavities instantaneously adjust to changes in slip and effective pressure forcings, but recent studies indicate transient forcings violate this—and other—underlying assumptions. To assess these incongruities, we conducted novel experiments emulating hard-bedded slip with ice-bed separation under periodic effective pressure transients. We slid an ice-ring over a sinusoidal bed while varying the applied overburden stress to emulate subglacial effective pressure cycles observed in nature and continuously recorded mechanical and geometric system responses. We observed characteristic lags and nonlinearities in system responses that were sensitive to forcing periodicity and trajectory. This gave rise to hysteresis not predicted in analytic theory, which we ascribed to a combination of geometric, thermal and rheologic processes. This framework corroborates other studies of transient glacier slip and we used it to place new constraints on transient phenomena observed in the field. Despite these divergences, average system responses converged toward model predictions, suggesting that analytic theory remains applicable for modeling longer-term behaviors of transiently forced slip with ice-bed separation.

Knowledge of how glaciations formed landscapes is particularly important as receding glaciers currently uncover subglacial landscapes that are prone to a series of natural hazards, but that also bring opportunities for hydropower or water resources. We present high-resolution (1 × 1 m) swath bathymetric data of a proglacial lake in front of the Rhonegletscher (Swiss Alps) that started to form in the early 2000s allowing a look into a freshly uncovered glacier bed and its characterized morphology in an overdeepened setting. The comparison of two surveys from 2015 and 2021 allows an unprecedented quantification of the accumulation and erosion processes in the central lake basin. This highly dynamic environment is characterized by iceberg calving, fluctuating outflow conduits, rapid sedimentation due to particle-laden meltwaters and dumped glacial debris. Assuming constant sediment yield, the Rhone Lake would persist for ~300 years. However, as intense glacier retreat continues in the coming decades, a chain of overdeepened lakes will be revealed that will act as long persisting sediment traps.

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.

Subglacial seismicity presents the opportunity to monitor inaccessible glacial beds at the epicentral location and time. Glaciers can be underlain by rock or till, a first order control on bed mechanics. Velocity-weakening, necessary for unstable slip, has been shown for each bed type, but is much stronger and evolves over more than an order of magnitude longer distances for till beds. Utilizing a de-stiffened double direct shear apparatus, we found conditions for instability at freezing temperatures and high slip rates for both bed types. During stick–slip stress-drops, we recorded acoustic emissions with piezoelectric transducers frozen into the ice. The two populations of event waveforms appear visually similar and overlap in their statistical features. We implemented a suite of supervised machine learning algorithms to classify the bed type of recorded waveforms and spectra, with prediction accuracy between 65–80%. The Random Forest Classifier is interpretable, showing the importance of initial oscillation peaks and higher frequency energy. Till beds have generally higher friction and resulting stress-drops, with more impulsive first arrivals and more high frequency content compared to rock emissions, but rock beds can produce many till-like events. Seismic signatures could enhance interpretation of bed conditions and mechanics from subglacial seismicity.

Ice sheet models use observations to infer basal shear stress, but the variety of methods and datasets available has resulted in a wide range of estimates. Radar-based metrics such as reflectivity and specularity content have been used to characterize subglacial hydrologic conditions that are linked to spatial variations in basal shear stress. We explore whether radar metrics can be used to inform models about basal shear stress. At Thwaites Glacier, West Antarctica, we sample basal shear stress inversions across a wide range of ice sheet models to see how the basal shear stress distribution changes in regions of varying relative reflectivity and specularity content. Our results reveal three key findings: (1) Regions of high specularity content exhibit lower mean basal shear stresses (2) Wet and bumpy regions, as characterized by high relative reflectivity and low specularity content, exhibit higher mean basal shear stresses (3) Models disagree about what basal shear stress should be at the onset of rapid ice flow and high basal melt where relative reflectivity and specularity content are low.

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.

Subglacial lake water-volume changes produce ice-elevation anomalies that provide clues about water flow beneath glaciers and ice sheets. Significant challenges remain in the quantitative interpretation of these elevation-change anomalies because the surface expression of subglacial lake activity depends on basal conditions, rate of water-volume change, and ice rheology. To address these challenges, we introduce an inverse method that reconstructs subglacial lake activity from altimetry data while accounting for the effects of viscous ice flow. We use a linearized approximation of a Stokes ice-flow model under the assumption that subglacial lake activity only induces small perturbations relative to a reference ice-flow state. We validate this assumption by accurately reconstructing lake activity from synthetic data that are produced with a fully nonlinear model. We then apply the method to estimate the water-volume changes of several active subglacial lakes in Antarctica by inverting data from NASA's Ice, Cloud, and land Elevation Satellite 2 (ICESat-2) laser altimetry mission. The results show that there can be substantial discrepancies (20% or more) between the inversion and traditional estimation methods due to the effects of viscous ice flow. The inverse method will help refine estimates of subglacial water transport and further constrain the role of subglacial hydrology in ice-sheet evolution.

Spatial variability in bed topography, characterized as bed roughness, impacts ice-sheet flow and organization and can be used to infer subglacial conditions and processes, yet is difficult to quantify due to sparse observations. Paleo-subglacial beds of formerly expanded glaciers found across the Antarctic continental shelf are well preserved, have relatively limited post-glacial sediment cover and contain glacial landforms that can be resolved at sub-meter vertical scales. We analyze high-resolution bathymetry offshore of Pine Island and Thwaites glaciers in the Amundsen Sea to explore spatial variability of bed roughness where streamlined subglacial landforms allow for the determination of ice-flow direction. We quantify bed roughness using std dev. and Fast Fourier Transform methods, each employed at local (100 km) and regional (101–2 km) scales and in along- and across-flow orientations to determine roughness expressions across spatial scales. We find that the magnitude of roughness is impacted by the parameters selected – which are often not sufficiently reported in studies – to quantify roughness. Important spatial patterns can be discerned from high-resolution bathymetry, highlighting both its usefulness in identifying patterns of streaming ice flow and underscores the need for a standardized way of characterizing topographic variability.

We present a method for tracking radio-tagged pebbles and cobbles through subglacial meltwater channels under shallow temperate glaciers. Natural particles tagged with active radio transmitters were injected directly into a large subglacial channel 300 m up-glacier from the terminus of the Glacier d'Otemma, Switzerland. A roving antenna was developed to localise tagged particles planimetrically in subglacial and proglacial channel reaches (350 and 150 m long, respectively) using a probabilistic technique, delivering records of the change in particle location and transport distance over time with uncertainty. The roving antenna had a ±5−15 m planimetric precision, a 75% particle localisation rate and operated at a maximum ice depth of 47 m. Additionally, stationary supraglacial and proglacial antennas continuously monitored the passage of tagged particles through consecutive reaches of the channel, constraining the timing of particle transport events. The proglacial antenna system had a 98.1% detection rate and was operational to 0.89 m water depth during testing. Roving and stationary antenna records were combined to create a transport distance model for each particle, which may be used in conjunction with hydraulic data to investigate the kinematics of particle motion. When applied at scale in future studies, this method may be used to reveal the mechanisms and timescales of coarse sediment export from Alpine glaciers.

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.

Hydrologic pathways beneath ice sheets and glaciers play an important role in regulating ice flow. Antarctica has experienced, and will continue to experience, changes in ice dynamics and geometry, but the associated changes in subglacial hydrology have received less attention. Here, we use the GlaDS subglacial hydrology model to examine drainage evolution beneath an idealised Antarctic glacier in response to steepening ice surface slopes, accelerating ice velocities and subglacial lake drainages. Ice surface slope changes exerted a dominant influence, redirecting basal water to different outlet locations and substantially increasing channelised discharge crossing the grounding line. Faster ice velocities had comparatively negligible effects. Subglacial lake drainage results indicated that lake refilling times play a key role in drainage system evolution, with lake flux more readily accommodated following shorter refilling times. Our findings are significant for vulnerable Antarctic regions currently experiencing dynamic thinning since subglacial water re-routing could destabilise ice shelves through enhanced sub-shelf melting, potentially hastening irreversible retreat. These changes could also affect subglacial lake activity. We, therefore, emphasise that including a nuanced and complex representation of subglacial hydrology in ice-sheet models could provide critical information on the timing and magnitude of sea-level change contributions from Antarctica.

Glacier-bed characteristics that are poorly known and modeled are important in projected sea-level rise from ice-sheet changes under strong warming, especially in the Thwaites Glacier drainage of West Antarctica. Ocean warming may induce ice-shelf thinning or loss, or thinning of ice in estuarine zones, reducing backstress on grounded ice. Models indicate that, in response, more-nearly-plastic beds favor faster ice loss by causing larger flow acceleration, but more-nearly-viscous beds favor localized near-coastal thinning that could speed grounding-zone retreat into interior basins where marine-ice-sheet instability or cliff instability could develop and cause very rapid ice loss. Interpretation of available data indicates that the bed is spatially mosaicked, with both viscous and plastic regions. Flow against bedrock topography removes plastic lubricating tills, exposing bedrock that is eroded on up-glacier sides of obstacles to form moats with exposed bedrock tails extending downglacier adjacent to lee-side soft-till bedforms. Flow against topography also generates high-ice-pressure zones that prevent inflow of lubricating water over distances that scale with the obstacle size. Extending existing observations to sufficiently large regions, and developing models assimilating such data at the appropriate scale, present large, important research challenges that must be met to reliably project future forced sea-level rise.

Bed topography is a critical parameter for determining the modern-day and future dynamics of ice sheets and their outlet glaciers. This is because the topography controls the state of stress for glaciers. At glacier termini, topography can influence the timing of terminus retreat by controlling access to warm ocean waters and/or by influencing the ability of a glacier terminus to retreat over bed bumps (moraines). Inland from the terminus, the topography can also influence where glacier retreat and thinning can stabilize. In part, this is because of knickpoints in bed topography created through glacial erosion that may influence the extent to which thinning can diffuse inland for an individual glacier and thus, the timing and magnitude of long-term mass loss. Here we provide a review of the current literature on these topics. While much of the reviewed literature assumes that topography is stable on relevant timescales to humans, new research suggests that topography may change much faster than previously thought and this further complicates our ability to project future outlet glacier change.

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.