Modelling water flow under glaciers and ice sheets

Linear-viscous flow of temperate ice

Accurately modeling the deformation of temperate glacier ice, which is at its pressure-melting temperature and contains liquid water at grain boundaries, is essential for predicting ice sheet discharge to the ocean and associated sea-level rise. Central to such modeling is Glen's flow law, in which strain rate depends on stress raised to a power of n = 3 to 4. In sharp contrast to this nonlinearity, we found by conducting large-scale, shear-deformation experiments that temperate ice is linear-viscous (n ≈ 1.0) over common ranges of liquid water content and stress expected near glacier beds and in ice-stream margins. This linearity is likely caused by diffusive pressure melting and refreezing at grain boundaries and could help to stabilize modeled responses of ice sheets to shrinkage-induced stress increases.

A Consistent Framework for Coupling Basal Friction With Subglacial Hydrology on Hard-Bedded Glaciers

Below hard-bedded glaciers, both basal friction and distributed subglacial drainage are thought to be controlled by a network of cavities. Previous coupled hydro-mechanical models, however, describe cavity-driven friction and hydraulic transmissivity independently, resulting in a physically inconsistent cavity evolution between the two components of the models. Here, we overcome this issue by describing the hydro-mechanical system using a common cavity-evolution description, that governs both transient friction and hydraulic transmissivity. We show that our coupling approach is superior to previous formulations in explaining a unique observation record of glacier sliding speed from the French Alps. We find that, at multi-day to multi-decadal timescales, sliding speed can be expressed as a direct function of basal shear stress and water discharge, without accounting for water pressure, which simply adjusts to maintain the cavitation ratio needed to accommodate the water supply.

Radar-Sounding Characterization of the Subglacial Groundwater Table Beneath Hiawatha Glacier, Greenland

Radar-sounding surveys associated with the discovery of a large impact crater beneath Hiawatha Glacier, Greenland, revealed bright, flat subglacial reflections hypothesized to originate from a subglacial groundwater table. We test this hypothesis using radiometric and hydrologic analysis of those radar data. The dielectric loss between the reflection from the top of the basal layer and subglacial reflection and their reflectivity difference represent dual constraints upon the complex permittivity of the basal material. Either ice-cemented debris or fractured, well-drained bedrock explain the basal layer's radiometric properties. The subglacial reflector's geometry is parallel to isopotential hydraulic head contours, located 7.5-15.3 m below the interface, and 11 ± 7 dB brighter than the ice-basal layer reflection. We conclude that this subglacial reflection is a groundwater table and that its detection was enabled by the wide bandwidth of the radar system and unusual geologic setting, suggesting a path for future direct radar detection of subglacial groundwater elsewhere.

Spatial heterogeneity in subglacial drainage driven by till erosion

The distribution and drainage of meltwater at the base of glaciers sensitively affects fast ice flow. Previous studies suggest that thin meltwater films between the overlying ice and a hard-rock bed channelize into efficient drainage elements by melting the overlying ice. However, these studies do not account for the presence of soft deformable sediment observed underneath many West Antarctic ice streams, and the inextricable coupling that sediment exhibits with meltwater drainage. Our work presents an alternate mechanism for initiating drainage elements such as canals where meltwater films grow by eroding the sediment beneath. We conduct a linearized stability analysis on a meltwater film flowing over an erodible bed. We solve the Orr–Sommerfeld equation for the film flow, and we compute bed evolution with the Exner equation. We identify a regime where the coupled dynamics of hydrology and sediment transport drives a morphological instability that generates spatial heterogeneity at the bed. We show that this film instability operates at much faster time scales than the classical thermal instability proposed by Walder. We discuss the physics of the instability using the framework of ripple formation on erodible beds.

Evidence for a palaeo-subglacial lake on the Antarctic continental shelf

Subglacial lakes are widespread beneath the Antarctic Ice Sheet but their control on ice-sheet dynamics and their ability to harbour life remain poorly characterized. Here we present evidence for a palaeo-subglacial lake on the Antarctic continental shelf. A distinct sediment facies recovered from a bedrock basin in Pine Island Bay indicates deposition within a low-energy lake environment. Diffusive-advection modelling demonstrates that low chloride concentrations in the pore water of the corresponding sediments can only be explained by initial deposition of this facies in a freshwater setting. These observations indicate that an active subglacial meltwater network, similar to that observed beneath the extant ice sheet, was also active during the last glacial period. It also provides a new framework for refining the exploration of these unique environments.

While a number of subglacial lakes have been discovered beneath the modern Antarctic ice sheet, little is known regarding their prevalence during past glaciations. Here, Kuhn et al. present sedimentological evidence for an active palaeo-subglacial lake in Pine Island Bay during the last glaciation.

Greenland subglacial drainage evolution regulated by weakly connected regions of the bed

Penetration of surface meltwater to the bed of the Greenland Ice Sheet each summer causes an initial increase in ice speed due to elevated basal water pressure, followed by slowdown in late summer that continues into fall and winter. While this seasonal pattern is commonly explained by an evolution of the subglacial drainage system from an inefficient distributed to efficient channelized configuration, mounting evidence indicates that subglacial channels are unable to explain important aspects of hydrodynamic coupling in late summer and fall. Here we use numerical models of subglacial drainage and ice flow to show that limited, gradual leakage of water and lowering of water pressure in weakly connected regions of the bed can explain the dominant features in late and post melt season ice dynamics. These results suggest that a third weakly connected drainage component should be included in the conceptual model of subglacial hydrology.

Surface meltwater draining to the bed of the Greenland Ice Sheet each summer causes ice flow changes inconsistent with the prevailing theory of channelizing subglacial drainage. Here, the authors show this is caused by limited, gradual leakage of water from previously ignored weakly connected regions of the bed.

Advances in modelling subglacial lakes and their interaction with the Antarctic ice sheet

Subglacial lakes have long been considered hydraulically isolated water bodies underneath ice sheets. This view changed radically with the advent of repeat-pass satellite altimetry and the discovery of multiple lake discharges and water infill, associated with water transfer over distances of more than 200 km. The presence of subglacial lakes also influences ice dynamics, leading to glacier acceleration. Furthermore, subglacial melting under the Antarctic ice sheet is more widespread than previously thought, and subglacial melt rates may explain the availability for water storage in subglacial lakes and water transport. Modelling of subglacial water discharge in subglacial lakes essentially follows hydraulics of subglacial channels on a hard bed, where ice sheet surface slope is a major control on triggering subglacial lake discharge. Recent evidence also points to the development of channels in deformable sediment in West Antarctica, with significant water exchanges between till and ice. Most active lakes drain over short time scales and respond rapidly to upstream variations. Several Antarctic subglacial lakes exhibit complex interactions with the ice sheet due to water circulation. Subglacial lakes can therefore-from a modelling point of view-be seen as confined small oceans underneath an imbedded ice shelf.

A decade of progress in observing and modelling Antarctic subglacial water systems

In the decade since the discovery of active Antarctic subglacial water systems by detection of subtle surface displacements, much progress has been made in our understanding of these dynamic systems. Here, we present some of the key results of observations derived from ICESat laser altimetry, CryoSat-2 radar altimetry, Operation IceBridge airborne laser altimetry, satellite image differencing and ground-based continuous Global Positioning System (GPS) experiments deployed in hydrologically active regions. These observations provide us with an increased understanding of various lake systems in Antarctica: Whillans/Mercer Ice Streams, Crane Glacier, Recovery Ice Stream, Byrd Glacier and eastern Wilkes Land. In several cases, subglacial water systems are shown to control ice flux through the glacier system. For some lake systems, we have been able to construct more than a decade of continuous lake activity, revealing internal variability on time scales ranging from days to years. This variability indicates that continuous, accurate time series of altimetry data are critical to understanding these systems. On Whillans Ice Stream, our results from a 5-year continuous GPS record demonstrate that subglacial lake flood events significantly change the regional ice dynamics. We also show how models for subglacial water flow have evolved since the availability of observations of lake volume change, from regional-scale models of water routeing to process models of channels carved into the subglacial sediment instead of the overlying ice. We show that progress in understanding the processes governing lake drainage now allows us to create simulated lake volume time series that reproduce time series from satellite observations. This transformational decade in Antarctic subglacial water research has moved us significantly closer to understanding the processes of water transfer sufficiently for inclusion in continental-scale ice-sheet models.

The distribution of basal water between Antarctic subglacial lakes from radar sounding

Antarctica's subglacial lakes have two end member geophysical expressions: as hydraulically flat, radar reflective regions highlighted in ice surface topography and radar sounding profiles ('definite lakes'), and as localized sites of elevation change identified from repeat elevation observations ('active lakes') that are often found in fast flowing ice streams or enhanced ice flow tributaries. While 'definite lakes' can be identified readily by high bed reflectivity in radar sounding, the identification and characterization of less distinct subglacial lakes and water systems with radar sounding are complicated by variable radio-wave attenuation in the overlying ice. When relying on repeat elevation observations, the relatively short times series and biased distribution of elevation observations, along with the episodic nature of 'active lake' outflow and replenishment, limit our understanding of how water flows under the ice sheet. Using recently developed methods for quantifying the radar scattering behaviour of the basal interface of the ice, we can avoid the problem of attenuation, and observe the plumbing of the subglacial landscape. In West Antarctica's Ross Sea Embayment, we confirm that extensive distributed water systems underlie these ice streams. Distributed water sheets are upstream in the onset regions of fast flow, while canal systems underly downstream regions of fast flow. In East Antarctica, we use specularity analysis to recover substantial hydraulic connectivity extending beyond previous knowledge, connecting the lakes already delineated by traditional radar sounding or surface elevation transients.