Subglacial Hydrology Research Articles

The feasibility of imaging subglacial hydrology beneath ice streams with ground-based electromagnetics

ABSTRACTSubglacial hydrologic systems in Antarctica and Greenland play a fundamental role in ice-sheet dynamics, yet critical aspects of these systems remain poorly understood due to a lack of observations. Ground-based electromagnetic (EM) geophysical methods are established for mapping groundwater in many environments, but have never been applied to imaging lakes beneath ice sheets. Here, we study the feasibility of passive- and active-source EM imaging for quantifying the nature of subglacial water systems beneath ice streams, with an emphasis on the interfaces between ice and basal meltwater, as well as deeper groundwater in the underlying sediments. We describe a suite of model studies that exam the data sensitivity as a function of ice thickness, water conductivity and hydrologic system geometry for models representative of a subglacial lake and a grounding zone estuary. We show that EM data are directly sensitive to groundwater and can image its lateral and depth extent. By combining the conductivity obtained from EM data with ice thickness and geological structure from conventional geophysical techniques, such as ground-penetrating radar and active seismic surveying, EM data have the potential to provide new insights on the interaction between ice, rock and water at critical ice-sheet boundaries.

Geothermal heat flux in the Amundsen Sea sector of West Antarctica: New insights from temperature measurements, depth to the bottom of the magnetic source estimation, and thermal modeling

AbstractFocused research on the Pine Island and Thwaites glaciers, which drain the West Antarctic Ice Shelf (WAIS) into the Amundsen Sea Embayment (ASE), revealed strong signs of instability in recent decades that result from variety of reasons, such as inflow of warmer ocean currents and reverse bedrock topography, and has been established as the Marine Ice Sheet Instability hypothesis. Geothermal heat flux (GHF) is a poorly constrained parameter in Antarctica and suspected to affect basal conditions of ice sheets, i.e., basal melting and subglacial hydrology. Thermomechanical models demonstrate the influential boundary condition of geothermal heat flux for (paleo) ice sheet stability. Due to a complex tectonic and magmatic history of West Antarctica, the region is suspected to exhibit strong heterogeneous geothermal heat flux variations. We present an approach to investigate ranges of realistic heat fluxes in the ASE by different methods, discuss direct observations, and 3‐D numerical models that incorporate boundary conditions derived from various geophysical studies, including our new Depth to the Bottom of the Magnetic Source (DBMS) estimates. Our in situ temperature measurements at 26 sites in the ASE more than triples the number of direct GHF observations in West Antarctica. We demonstrate by our numerical 3‐D models that GHF spatially varies from 68 up to 110 mW m−2.

Seasonal and Interannual Variability of Columbia Glacier, Alaska (2011–2016): Ice Velocity, Mass Flux, Surface Elevation and Front Position

Alaskan glaciers are among the largest contributors to sea-level rise outside the polar ice sheets. The contributions include dynamic discharge from marine-terminating glaciers which depends on the seasonally variable ice velocity. Columbia Glacier is a large marine-terminating glacier located in Southcentral Alaska that has been exhibiting pronounced retreat since the early 1980s. Since 2010, the glacier has split into two branches, the main branch and the west branch. We derived a 5-year record of surface velocity, mass flux (ice discharge), surface elevation and changes in front position using a dense time series of TanDEM-X synthetic aperture radar data (2011–2016). We observed distinct seasonal velocity patterns at both branches. At the main branch, the surface velocity peaked during late winter to midsummer but reached a minimum between late summer and fall. Its near-front velocity reached up to 14 m day−1 in May 2015 and was at its lowest speed of ~1 m day−1 in October 2012. Mass flux via the main branch was strongly controlled by the seasonal and interannual fluctuations of its velocity. The west branch also exhibited seasonal velocity variations with comparably lower magnitudes. The role of subglacial hydrology on the ice velocities of Columbia Glacier is already known from the published field measurements during summers of 1987. Our observed variability in its ice velocities on a seasonal basis also suggest that they are primarily controlled by the seasonal transition of the subglacial drainage system from an inefficient to an efficient and channelized drainage networks. However, abrupt velocity increase events for short periods (2014–2015 and 2015–2016 at the main branch, and 2013–2014 at the west branch) appear to be associated with strong near-front thinning and frontal retreat. This needs further investigation on the role of other potential controlling mechanisms. On the technological side, this study demonstrates the potential of high-resolution X-band SAR missions with a short revisit interval to examine glaciological variables and controlling processes.

Eurasian ice-sheet dynamics and sensitivity to subglacial hydrology

ABSTRACTIce-stream dynamics are strongly controlled by processes taking place at the ice/bed interface where subglacial water both lubricates the base and saturates any existing, underlying sediment. Large parts of the former Eurasian ice sheet were underlain by thick sequences of soft, marine sediments and many areas are imprinted with geomorphological features indicative of fast flow and wet basal conditions. Here, we study the effect of subglacial water on past Eurasian ice-sheet dynamics by incorporating a thin-film model of basal water flow into the ice-sheet model SICOPOLIS and use it to better represent flow in temperate areas. The adjunction of subglacial hydrology results in a smaller ice-sheet building up over time and generally faster ice velocities, which consequently reduces the total area fraction of temperate basal ice and ice streaming areas. Minima in the hydraulic pressure potential, governing water flow, are used as indicators for potential locations of past subglacial lakes and a probability distribution of lake existence is presented based on estimated lake depth and longevity.

Spatially distributed runoff at the grounding line of a large Greenlandic tidewater glacier inferred from plume modelling

ABSTRACTUnderstanding the drivers of recent change at Greenlandic tidewater glaciers is of great importance if we are to predict how these glaciers will respond to climatic warming. A poorly constrained component of tidewater glacier processes is the near-terminus subglacial hydrology. Here we present a novel method for constraining near-terminus subglacial hydrology with application to marine-terminating Kangiata Nunata Sermia in South-west Greenland. By simulating proglacial plume dynamics using buoyant plume theory and a general circulation model, we assess the critical subglacial discharge, if delivered through a single compact channel, required to generate a plume that reaches the fjord surface. We then compare catchment runoff to a time series of plume visibility acquired from a time-lapse camera. We identify extended periods throughout the 2009 melt season where catchment runoff significantly exceeds the discharge required for a plume to reach the fjord surface, yet we observe no plume. We attribute these observations to spatial spreading of runoff across the grounding line. Persistent distributed drainage near the terminus would lead to more spatially homogeneous submarine melting and may promote more rapid basal sliding during warmer summers, potentially providing a mechanism independent of ocean forcing for increases in atmospheric temperature to drive tidewater glacier acceleration.

Extensive winter subglacial water storage beneath the Greenland Ice Sheet

AbstractSurface meltwater that reaches the base of the Greenland Ice Sheet exerts a fundamental impact on ice flow, but observations of catchment‐wide movement and distribution of subglacial water remain limited. Using radar‐sounding data from two seasons, we identify the seasonal distribution of subglacial water in western Greenland. Our analysis provides evidence of widespread subglacial water storage beneath Greenland in the wintertime. The winter storage is located primarily on bedrock ridges with higher bed elevations in excess of 200 m. During the melt season water moves to the subglacial troughs. This inverse relationship with topography indicates that the material properties of the glacier bed strongly influence subglacial drainage development. Both the spatial variations in bed properties and the initial state of the subglacial hydrology system at the start of the melt season lead to differing glacier dynamical responses to surface melting across the Greenland Ice Sheet.

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.

Glacial landscape evolution by subglacial quarrying: A multiscale computational approach

AbstractQuarrying of bedrock is a primary agent of subglacial erosion. Although the mechanical theory behind the process has been studied for decades, it has proven difficult to formulate the governing principles so that large‐scale landscape evolution models can be used to integrate erosion over time. The existing mechanical theory thus stands largely untested in its ability to explain postglacial topography. In this study we relate the physics of quarrying to long‐term landscape evolution with a multiscale approach that connects meter‐scale cavities to kilometer‐scale glacial landscapes. By averaging the quarrying rate across many small‐scale bedrock steps, we quantify how regional trends in basal sliding speed, effective pressure, and bed slope affect the rate of erosion. A sensitivity test indicates that a power law formulated in terms of these three variables provides an acceptable basis for quantifying regional‐scale rates of quarrying. Our results highlight the strong influence of effective pressure, which intensifies quarrying by increasing the volume of the bed that is stressed by the ice and thereby the probability of rock failure. The resulting pressure dependency points to subglacial hydrology as a primary factor for influencing rates of quarrying and hence for shaping the bedrock topography under warm‐based glaciers. When applied in a landscape evolution model, the erosion law for quarrying produces recognizable large‐scale glacial landforms: U‐shaped valleys, hanging valleys, and overdeepenings. The landforms produced are very similar to those predicted by more standard sliding‐based erosion laws, but overall quarrying is more focused in valleys, and less effective at higher elevations.

Bathymetric evolution of Tasman Glacier terminal lake, New Zealand, as determined by remote surveying techniques

Processes that drive iceberg calving at the margins of freshwater terminating glaciers are still poorly understood. This knowledge-gap is in part due to the challenge of obtaining good in situ data in a highly dynamic and dangerous environment. We are using emerging remote technologies, in the form of a remote controlled jet boat to survey bathymetry, and Structure from Motion (SfM) to characterize terminus morphology, to better understand relationships between lake growth and terminus evolution. Comparison of results between the jet boat mounted dual-frequency Garmin fish-finder with an Odom Echotrac DF3200 MKII with 200/38kHz dual-frequency transducer, showed that after a sound velocity adjustment, the remote survey obtained depth data within ±1m of the higher grade survey equipment. Water depths of up to 240m were recorded only 100m away from the terminus, and subaerial cliff height ranged from around 6 to 33m, with the central region of the terminus more likely to experience buoyancy. Subaqueous ice ramps are ephemeral features, and in 2015 multiple ice ramps extended out into the lake from the terminus by 100–200m. The consistent location of some of the subaqueous ramps between surveys may indicate that other processes, for example, subglacial hydrology, also influence evolving terminus morphology.

Ice‐margin and meltwater dynamics during the mid‐Holocene in the Kangerlussuaq area of west Greenland

Land‐terminating parts of the west Greenland ice sheet have exhibited highly dynamic meltwater regimes over the last few decades including episodes of extremely intense runoff driven by ice surface ablation, ponding of meltwater in an increasing number and size of lakes, and sudden outburst floods, or ‘jökulhlaups’, from these lakes. However, whether this meltwater runoff regime is unusual in a Holocene context has not been questioned. This study assembled high‐resolution topographical data, geological and landcover data, and produced a glacial geomorphological map covering ~1200 km2. Digital analysis of the landforms reveals a mid‐Holocene land‐terminating ice margin that was predominantly cold‐based. This ice margin underwent sustained active retreat but with multiple minor advances. Over c. 1000 years meltwater runoff became impounded within numerous and extensive proglacial lakes and there were temporary connections between some of these lakes via spillways. The ice‐dams of some of these lakes had several quasi‐stable thicknesses. Meltwater was apparently predominantly from supraglacial sources although some distributary palaeochannel networks and some larger bedrock palaeochannels most likely relate to mid‐Holocene subglacial hydrology. In comparison to the geomorphological record at other Northern Hemisphere ice‐sheet margins the depositional landforms in this study area are few in number and variety and small in scale, most likely due to a restricted sediment supply. They include perched fans and deltas and perched braidplain terraces. Overall, meltwater sourcing, routing and the proglacial runoff regime during the mid‐Holocene in this land‐terminating part of the ice sheet was spatiotemporally variable, but in a manner very similar to that of the present day.

Modeling Antarctic subglacial lake filling and drainage cycles

Abstract. The growth and drainage of active subglacial lakes in Antarctica has previously been inferred from analysis of ice surface altimetry data. We use a subglacial hydrology model applied to a synthetic Antarctic ice stream to examine internal controls on the filling and drainage of subglacial lakes. Our model outputs suggest that the highly constricted subglacial environment of our idealized ice stream, combined with relatively high rates of water flow funneled from a large catchment, can combine to create a system exhibiting slow-moving pressure waves. Over a period of years, the accumulation of water in the ice stream onset region results in a buildup of pressure creating temporary channels, which then evacuate the excess water. This increased flux of water beneath the ice stream drives lake growth. As the water body builds up, it steepens the hydraulic gradient out of the overdeepened lake basin and allows greater flux. Eventually this flux is large enough to melt channels that cause the lake to drain. Lake drainage also depends on the internal hydrological development in the wider system and therefore does not directly correspond to a particular water volume or depth. This creates a highly temporally and spatially variable system, which is of interest for assessing the importance of subglacial lakes in ice stream hydrology and dynamics.

Greenland Ice Sheet seasonal and spatial mass variability from model simulations and GRACE (2003–2012)

Abstract. Improving the ability of regional climate models (RCMs) and ice sheet models (ISMs) to simulate spatiotemporal variations in the mass of the Greenland Ice Sheet (GrIS) is crucial for prediction of future sea level rise. While several studies have examined recent trends in GrIS mass loss, studies focusing on mass variations at sub-annual and sub-basin-wide scales are still lacking. At these scales, processes responsible for mass change are less well understood and modeled, and could potentially play an important role in future GrIS mass change. Here, we examine spatiotemporal variations in mass over the GrIS derived from the Gravity Recovery and Climate Experiment (GRACE) satellites for the January 2003–December 2012 period using a "mascon" approach, with a nominal spatial resolution of 100 km, and a temporal resolution of 10 days. We compare GRACE-estimated mass variations against those simulated by the Modèle Atmosphérique Régionale (MAR) RCM and the Ice Sheet System Model (ISSM). In order to properly compare spatial and temporal variations in GrIS mass from GRACE with model outputs, we find it necessary to spatially and temporally filter model results to reproduce leakage of mass inherent in the GRACE solution. Both modeled and satellite-derived results point to a decline (of −178.9 ± 4.4 and −239.4 ± 7.7 Gt yr−1 respectively) in GrIS mass over the period examined, but the models appear to underestimate the rate of mass loss, especially in areas below 2000 m in elevation, where the majority of recent GrIS mass loss is occurring. On an ice-sheet-wide scale, the timing of the modeled seasonal cycle of cumulative mass (driven by summer mass loss) agrees with the GRACE-derived seasonal cycle, within limits of uncertainty from the GRACE solution. However, on sub-ice-sheet-wide scales, some areas exhibit significant differences in the timing of peaks in the annual cycle of mass change. At these scales, model biases, or processes not accounted for by models related to ice dynamics or hydrology, may lead to the observed differences. This highlights the need for further evaluation of modeled processes at regional and seasonal scales, and further study of ice sheet processes not accounted for, such as the role of subglacial hydrology in variations in glacial flow.

Sliding velocity fluctuations and subglacial hydrology over the last two decades on Argentière glacier, Mont Blanc area

ABSTRACTThe subglacial observatory beneath the Argentière glacier provides a rare opportunity to study the interactions between glacier sliding velocity and subglacial runoff. The sliding velocity has been monitored in this cavity almost continuously since 1997 and the resulting data indicate a decrease in annual sliding velocities over the last two decades. We found close relationships between annual surface velocity, sliding velocity and ice thickness. These relationships indicate that the ice-flow velocity changes do not depend on subglacial water runoff changes at the annual timescale. The seasonal magnitudes of sliding also show a decrease over the last two decades. At the seasonal timescale, sliding velocity increases before or simultaneously with the large runoff increase in May, indicating a distributed drainage system. Conversely, at the end of the melt season, sliding velocity continues to decrease after the runoff returns to low winter values. The simultaneous increases of runoff and sliding velocity occur mainly before the spring transition. Later, sliding velocity generally appears not to be related to water inputs coming from the surface, except for some large accelerations after midAugust that are always associated with periods of rapidly increasing water inputs to the subglacial drainage system.

Inversion of a glacier hydrology model

ABSTRACTThe subglacial hydrologic system exerts strong controls on the dynamics of the overlying ice, yet the parameters that govern the evolution of this system are not widely known or observable. To gain a better understanding of these parameters, we invert a spatially averaged model of subglacial hydrology from observations of ice surface velocity and outlet stream discharge at Kennicott Glacier, Wrangell Mountains, AK, USA. To identify independent parameters, we formally non-dimensionalize the forward model. After specifying suitable prior distributions, we use a Markov-chain Monte Carlo algorithm to sample from the distribution of parameter values conditioned on the available data. This procedure gives us not only the most probable parameter values, but also a rigorous estimate of their covariance structure. We find that the opening of cavities due to sliding over basal topography and turbulent melting are of a similar magnitude during periods of large input flux, though turbulent melting also exhibits the greatest uncertainty. We also find that both the storage of water in the englacial system and the exchange of water between englacial and subglacial systems are necessary in order to explain both surface velocity observations and the relative attenuation in the amplitude of diurnal signals between input and output flux observations.

Modeling the WorldView-derived seasonal velocity evolution of Kennicott Glacier, Alaska

ABSTRACTGlacier basal motion generates diurnal to multi-annual fluctuations in glacier velocity and mass flux. Understanding these fluctuations is important for prediction of future sea-level rise and for gaining insight into glacier physics and erosion. Here, we derive glacier velocity through cross-correlation of WorldView satellite imagery to document the evolution of ice surface velocity on Kennicott Glacier, Alaska, over the 2013 melt season. The summer speedup is spatially uniform over a ~12 km2area, over which the spring velocity varies significantly. Velocity increases by 1.4-fold to tenfold across the study domain, with larger values where spring velocities are low. To investigate the cross-glacier distribution of basal motion required to explain the observed surface speedup, we employ a two-dimensional cross-sectional glacier flow model. We find the model is insensitive to the spatial distribution of basal slip because stress gradient ice coupling diffuses the surface expression of the basal velocity field. While the temporal evolution of the subglacial hydrologic system is critical for predicting a glacier's response to meltwater inputs, our work suggests that glacier and ice-sheet models do not require a detailed representation of subglacial hydrology to accurately capture the spatial pattern of glacier speedup.

Effects of ice deformation on Röthlisberger channels and implications for transitions in subglacial hydrology

ABSTRACTAlong the base of glaciers and ice sheets, the sliding of ice over till depends critically on water drainage. In locations where drainage occurs through Röthlisberger channels, the effective pressure along the base of the ice increases and can lead to a strengthening of the bed, which reduces glacier sliding. The formation of Röthlisberger channels depends on two competing effects: (1) melting from turbulent dissipation opens the channel walls and (2) creep flow driven by the weight of the overlying ice closes the channels radially inward. Variation in downstream ice velocity along the channel axis, referred to as an antiplane shear strain rate, decreases the effective viscosity. The softening of the ice increases creep closure velocities. In this way, even a modest addition of antiplane shear can double the size of the Röthlisberger channels for a fixed water pressure or allow channels of a fixed radius to operate at lower effective pressure, potentially decreasing the strength of the surrounding bed. Furthermore, we show that Röthlisberger channels can be deformed away from a circular cross section under applied antiplane shear. These results can have broad impacts on sliding velocities and potentially affect the total ice flux out of glaciers and ice streams.

Subglacial hydrology at Rink Isbræ, West Greenland inferred from sediment plume appearance

ABSTRACTMarine-terminating outlet glaciers discharge most of the Greenland ice sheet's mass through frontal ablation and meltwater runoff. While calving can be estimated by in situ and remote sensing observations, submarine melting and subglacial meltwater transport are more challenging to quantify. Here we investigate the subglacial hydrology of Rink Isbræ, a fast-flowing West Greenland tidewater glacier, using time-lapse photography, modeled runoff estimates and daily satellite imagery from 2007 to 2011. We find that sediment plumes appear episodically at four distinct locations across the terminus, and last between 2 h and 17 d. This suggests short-term variability in discharge and the existence of persistent pathways. The seasonal onset of sediment plumes occurs before supraglacial lake drainages, shortly after the onset of runoff, and only after the wintertime ice mélange has begun disintegrating. Plumes were also visible after the cessation of runoff (23 ± 5 d), which is indicative of subglacial storage. The lack of either a seasonal velocity change or a correspondence between meltwater availability and plume occurrence suggests that the subglacial system persists in a state of inefficient drainage. Subglacial hydrology at tidewater glaciers is of critical importance in understanding dynamics at the ice front.

Linking glacially modified waters to catchment-scale subglacial discharge using autonomous underwater vehicle observations

Abstract. Measurements of near-ice (< 200 m) hydrography and near-terminus subglacial hydrology are lacking, due in large part to the difficulty in working at the margin of calving glaciers. Here we pair detailed hydrographic and bathymetric measurements collected with an autonomous underwater vehicle as close as 150 m from the ice–ocean interface of the Saqqarliup sermia–Sarqardleq Fjord system, West Greenland, with modeled and observed subglacial discharge locations and magnitudes. We find evidence of two main types of subsurface glacially modified water (GMW) with distinct properties and locations. The two GMW locations also align with modeled runoff discharged at separate locations along the grounded margin corresponding with two prominent subcatchments beneath Saqqarliup sermia. Thus, near-ice observations and subglacial discharge routing indicate that runoff from this glacier occurs primarily at two discrete locations and gives rise to two distinct glacially modified waters. Furthermore, we show that the location with the largest subglacial discharge is associated with the lighter, fresher glacially modified water mass. This is qualitatively consistent with results from an idealized plume model.

Complex Greenland outlet glacier flow captured.

The Greenland Ice Sheet is losing mass at an accelerating rate due to increased surface melt and flow acceleration in outlet glaciers. Quantifying future dynamic contributions to sea level requires accurate portrayal of outlet glaciers in ice sheet simulations, but to date poor knowledge of subglacial topography and limited model resolution have prevented reproduction of complex spatial patterns of outlet flow. Here we combine a high-resolution ice-sheet model coupled to uniformly applied models of subglacial hydrology and basal sliding, and a new subglacial topography data set to simulate the flow of the Greenland Ice Sheet. Flow patterns of many outlet glaciers are well captured, illustrating fundamental commonalities in outlet glacier flow and highlighting the importance of efforts to map subglacial topography. Success in reproducing present day flow patterns shows the potential for prognostic modelling of ice sheets without the need for spatially varying parameters with uncertain time evolution.

High basal melting forming a channel at the grounding line of Ross Ice Shelf, Antarctica

AbstractAntarctica's ice shelves are thinning at an increasing rate, affecting their buttressing ability. Channels in the ice shelf base unevenly distribute melting, and their evolution provides insight into changing subglacial and oceanic conditions. Here we used phase‐sensitive radar measurements to estimate basal melt rates in a channel beneath the currently stable Ross Ice Shelf. Melt rates of 22.2 ± 0.2 m a−1 (>2500% the overall background rate) were observed 1.7 km seaward of Mercer/Whillans Ice Stream grounding line, close to where subglacial water discharge is expected. Laser altimetry shows a corresponding, steadily deepening surface channel. Two relict channels to the north suggest recent subglacial drainage reorganization beneath Whillans Ice Stream approximately coincident with the shutdown of Kamb Ice Stream. This rapid channel formation implies that shifts in subglacial hydrology may impact ice shelf stability.