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.

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.

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.

Over half of the recent mass loss from the Greenland ice sheet, and its associated contribution to global sea level rise, can be attributed to increased surface meltwater runoff, with the remainder a result of dynamical processes such as calving and ice discharge. It is therefore important to quantify the distribution of melting on the ice sheet if we are to adequately understand past ice sheet change and make predictions for the future. In this article, we present a novel semi-empirical approach for characterising ice sheet surface conditions using high-resolution synthetic aperture radar (SAR) backscatter data from the Sentinel-1 satellite. We apply a state-space model to nine sites within North-East Greenland to identify changes in SAR backscatter, and we attribute these to different surface types with reference to optical satellite imagery and meteorological data. A set of decision-making rules for labelling ice sheet melting states are determined based on this analysis and subsequently applied to previously unseen sites. We show that our method performs well in (1) recognising some of the ice sheet surface types such as snow and dark ice and (2) determining whether the surface is melting or not melting. Sentinel-1 SAR data are of high spatial resolution; thus, in developing a method to identify the state of the surface from these data, we improve our capability to understand the variation of ice sheet melting across time and space.

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.

Subglacial melt has important implications for ice-sheet dynamics. Locating and identifying subglacial lakes are expensive and time-consuming, requiring radar surveys or satellite methods. We explore three methods to identify source regions for lakes using seven continent-wide environmental characteristics that are sensitive to or influenced by ice-sheet temperature. A simple comparison of environmental properties at lake locations with their continent-wide distributions suggests a statistical relationship (high Kolmogorov-Smirnov statistic) between stable lake locations and ice thickness and surface temperatures, indicating melting under passive conditions. Active lakes, in contrast, show little correlation with direct thermally influenced parameters, instead exhibiting large statistical differences with horizontal velocity and bedrock elevation. More sophisticated techniques, including principal component analysis (PCA) and machine learning (ML) classification, provide better spatial identification of lake types. Positive PCA scores derived from the environmental characteristics correlate with stable lakes, whereas negative values correspond to active lakes. ML methods can also identify regions where subglacial lake melt sources are probable. While ML provides the most accurate classification maps, the combination of approaches adds deeper knowledge of the primary controls on lake formation and the environmental settings in which they are likely to be found.

Lake Vostok, East Antarctica, represents an extensive water surface at the base of the ice sheet. Snow, ice and atmospheric pressure loads applied anywhere within the lake area produce a hydrostatic response, involving deformations of the ice surface, ice–water interface and particle horizons. A modelling scheme is developed to derive height changes of these surfaces for a given load pattern. It is applied to a series of load scenarios, and predictions based on load fields derived from a regional climate model are compared to observational datasets. Our results show that surface height changes due to snow-buildup anomalies are damped over the lake area, reducing the spatial standard deviation by one-third. The response to air pressure variations, in turn, adds surface height variability. Atmospheric pressure loads may produce height changes of up to $\pm$4 cm at daily resolution, but decay rapidly with integration time. The hydrostatic load response has no significant impact neither on ICESat laser campaign biases determined over the lake area nor on vertical particle movements derived from GNSS observations.

The earliest airborne geophysical campaigns over Antarctica and Greenland in the 1960s and 1970s collected ice penetrating radar data on 35 mm optical film. Early subglacial topographic and englacial stratigraphic analyses of these data were foundational to the field of radioglaciology. Recent efforts to digitize and release these data have resulted in geometric and ice-thickness analysis that constrain subsurface change over multiple decades but stop short of radiometric interpretation. The primary challenge for radiometric analysis is the poorly-characterized compression applied to Z-scope records and the sparse sampling of A-scope records. Here, we demonstrate the information richness and radiometric interpretability of Z-scope records. Z-scope pixels have uncalibrated fast-time, slow-time, and intensity scales. We develop approaches for mapping each of these scales to physical units (microseconds, seconds, and signal to noise ratio). We then demonstrate the application of this calibration and analysis approach to a flight in the interior of East Antarctica with subglacial lakes and to a reflight of an East Antarctic ice shelf that was observed by both archival and modern radar. These results demonstrate the potential use of Z-scope signals to extend the baseline of radiometric observations of the subsurface by decades.

The Subglacial Antarctic Lakes Scientific Access (SALSA) Project accessed Mercer Subglacial Lake using environmentally clean hot-water drilling to examine interactions among ice, water, sediment, rock, microbes and carbon reservoirs within the lake water column and underlying sediments. A ~0.4 m diameter borehole was melted through 1087 m of ice and maintained over ~10 days, allowing observation of ice properties and collection of water and sediment with various tools. Over this period, SALSA collected: 60 L of lake water and 10 L of deep borehole water; microbes >0.2 μm in diameter from in situ filtration of ~100 L of lake water; 10 multicores 0.32–0.49 m long; 1.0 and 1.76 m long gravity cores; three conductivity–temperature–depth profiles of borehole and lake water; five discrete depth current meter measurements in the lake and images of ice, the lake water–ice interface and lake sediments. Temperature and conductivity data showed the hydrodynamic character of water mixing between the borehole and lake after entry. Models simulating melting of the ~6 m thick basal accreted ice layer imply that debris fall-out through the ~15 m water column to the lake sediments from borehole melting had little effect on the stratigraphy of surficial sediment cores.

Drilling and sampling are the most direct and effective methods available to study Antarctic subglacial lakes. Based on the Philberth probe, a Recoverable Autonomous Sonde (RECAS) allows for in situ lake water measurement and sampling, through the addition of an upper thermal tip and a cable recoiling mechanism. RECAS-200, a prototype of RECAS, has a drilling depth of 200 m, a surface supply voltage of 800 VAC and a downhole power of ~9.6 kW during drilling. In this study, a heating control system for RECAS-200 was designed. The system avoids the need for high-power step-down converters, by separating heating power from control power, thereby reducing the overall weight of the probe and avoiding the need to increase cable diameter. We also introduce a self-developed, small, solid-state, 800 VAC power regulator and a fuzzy PID temperature control algorithm. Their purpose was to manage the power adjustment of each heating element and to provide closed-loop temperature control of certain heating elements which can easily burn out due to overheating. Test results indicated that the proposed RECAS-200 heating control system met all our design specifications and could be easily assembled into the RECAS-200 probe.

Recent drilling successes on Rutford Ice Stream in West Antarctica demonstrate the viability of hot water drilling subglacial access holes to depths >2000 m. Having techniques to access deep subglacial environments reliably paves the way for subglacial lake exploration beneath the thick central West Antarctic Ice Sheet. An ideal candidate lake, overlain by ~2650 m of ice, identified by Centro de Estudios Científicos (CECs), Chile, has led to collaboration with British Antarctic Survey to access Subglacial Lake CECs (SLCECs). To conform with the Scientific Committee on Antarctic Research code of conduct, which provides a guide to responsible scientific exploration and stewardship of these pristine systems, any access drilling must minimise all aspects of contamination and disturbance of the subglacial environment. To meet these challenges, along with thicker ice and 2000 m elevation, pumping and water treatment systems developed for the Subglacial Lake Ellsworth project, together with new diesel generators, additional water heating and longer drill hose, are currently being integrated with the BEAMISH hot water drill. A dedicated test season near SLCECs will commission the new clean hot water drill, with testing and validation of all clean operating procedures. A subsequent season will then access SLCECs cleanly.

A programmable borehole measurement system was deployed in hot water drilled ice holes during the ‘Bed Access and Monitoring of Ice Sheet History’ (BEAMISH) project to drill to the bed of the Rutford Ice Stream in West Antarctica. This system operates autonomously (no live data) after deployment, and records borehole diameter (non-contact measurement), water column pressure, heading and inclination. Three cameras, two sideways looking and one vertical, are also included for visual inspection of hole integrity and sediments. The system is small, lightweight (~35.5 kg) and low power using only 6 ‘D’ cell sized lithium batteries, making it ideal for transport and use in remote field sites. The system is 2.81 m long and 165 mm in diameter, and can be deployed attached to the drill hose for measurements during drilling or on its own deployment line afterwards. The full system is discussed in detail, highlighting design strengths and weaknesses. Data from the BEAMISH project are also presented in the form of camera images showing hole integrity, and sensor data used to calculate borehole diameter through the full length of the hole. These data are used to show confidence in hole verticality and subsurface cavity development and connection.

Monitoring the tension in cables is significant in some ice drill and deepwater applications. Take our RECoverable Autonomous Sonde (RECAS) for example. It is able to melt a hole to the bottom of ice sheet and is able to move upwards. A winch is installed inside RECAS to release and recover the cable, whose tension needs to be monitored in real time in order to control the behavior of the winch. The high pressure of deep water and limited installation space pose great challenges in sensor development. In this paper, two editions of newly designed deepwater tension sensors are proposed. The first edition is based on a fresh hydraulic load module that operates in high pressure environment and the second edition tension, which aims to improve the accuracy, applies a newly designed watertight load module. Detailed force transmission and characteristic analysis of the sensors are carried out. The sensors have got through a series of experiments, including calibration experiments, pressure experiments and field experiments. The resultant accuracy of the second edition sensor, which has a better performance, is over 2% under the measuring range of 1000 kg and the dimension of the final sensor is as compact as 150 mm × 137 mm × 86 mm.

Deep-water ‘stable’ subglacial lakes likely contain microbial life adapted in isolation to extreme environmental conditions. How water is supplied into a subglacial lake, and how water outflows, is important for understanding these conditions. Isochronal radio-echo layers have been used to infer where melting occurs above Lake Vostok and Lake Concordia in East Antarctica but have not been used more widely. We examine englacial layers above and around Lake Ellsworth, West Antarctica, to establish where the ice sheet is ‘drawn down’ towards the bed and, thus, experiences melting. Layer drawdown is focused over and around the northwest parts of the lake as ice, flowing obliquely to the lake axis becomes afloat. Drawdown can be explained by a combination of basal melting and the Weertman effect, at the transition from grounded to floating ice. We evaluate the importance of these processes on englacial layering over Lake Ellsworth and discuss implications for water circulation and sediment deposition. We report evidence of a second subglacial lake near the head of the hydrological catchment and present a new high-resolution bed DEM and hydropotential model of the lake outlet zone. These observations provide insight into the connectivity between Lake Ellsworth and the wider subglacial hydrological system.

All radar power interpretations require a correction for attenuative losses. Moreover, radar attenuation is a proxy for ice-column properties, such as temperature and chemistry. Prior studies use either paired thermodynamic and conductivity models or the radar data themselves to calculate attenuation, but there is no standard method to do so; and, before now, there has been no robust methodological comparison. Here, we develop a framework meant to guide the implementation of empirical attenuation methods based on survey design and regional glaciological conditions. We divide the methods into the three main groups: (1) those that infer attenuation from a single reflector across many traces; (2) those that infer attenuation from multiple reflectors within one trace; and (3) those that infer attenuation by contrasting the measured power from primary and secondary reflections. To assess our framework, we introduce a new ground-based radar survey from South Pole Lake, comparing selected empirical methods to the expected attenuation from a temperature- and chemistry-dependent Arrhenius model. Based on the small surveyed area, lack of a sufficient calibration surface and low reflector relief, the attenuation methods that use multiple reflectors are most suitable at South Pole Lake.

Late Holocene sediment deposits in Pine Island Bay, West Antarctica, are hypothesized to be linked to intensive meltwater drainage during the retreat of the paleo-Pine Island Ice Stream after the Last Glacial Maximum. The uppermost sediment units show an abrupt transition from ice-proximal debris to a draped silt during the late Holocene, which is interpreted to coincide with rapid deglaciation. The small scale and fine sorting of the upper unit could be attributed to origins in subglacial meltwater; however the thickness and deposition rate for this unit imply punctuated- rather than continuous-deposition. This, combined with the deposit's location seaward of large, bedrock basins, has led to the interpretation of this unit as the result of subglacial lake outbursts in these basins. However, the fine-scale sorting of the silt unit is problematic for this energetic interpretation, which should mobilize and deposit a wider range of sediment sizes. To resolve this discrepancy, we present an alternative mechanism in which the silt was sorted by a distributed subglacial water system, stored in bedrock basins far inland of the grounding line, and subsequently eroded at higher flow speeds during retreat. We demonstrate that this mechanism is physically plausible given the subglacial conditions during the late Holocene. We hypothesize that similar silt units observed elsewhere in Antarctica downstream of bedrock basins could be the result of the same mechanism.

Mountain glaciers and ice sheets often host marginal and subglacial lakes that are hydraulically connected through subglacial drainage systems. These lakes exhibit complex dynamics that have been the subject of models for decades. Here we introduce and analyze a model for the evolution of glacial lakes connected by subglacial channels. Subglacial channel equations are supplied with effective pressure boundary conditions that are determined by a simple lake model. While the model can describe an arbitrary number of lakes, we solve it numerically with a finite element method for the case of two connected lakes. We examine the effect of relative lake size and spacing on the oscillations. Complex oscillations in the downstream lake are driven by discharge out of the upstream lake. These include multi-peaked and anti-phase filling–draining events. Similar filling–draining cycles have been observed on the Kennicott Glacier in Alaska and at the confluence of the Whillans and Mercer ice streams in West Antarctica. We further construct a simplified ordinary differential equation model that displays the same qualitative behavior as the full, spatially-dependent model. We analyze this model using dynamical systems theory to explain the appearance of filling–draining cycles as the meltwater supply varies.

Geophysicists that deploy seismic sensors in ablation zones of glaciers and ice sheets to record glaciogenic signatures can confront recording challenges caused by instrument melt-out or tilt. These challenges often require installing sensors in boreholes to delay melt-out, or securing sensors to structures that improve coupling. We show that some of these structures that were buried near a moulin at a snow-free site in the ablation zone of the Western Greenland Ice Sheet resonated as they became exposed, and caused their geophones to record temporally evolving, narrowband signals that mimic features of glaciogenic sources like moulin tremor. We quantify these artifacts with a mechanical model that shows instruments undergo structural resonance as they melt-out, at exposure rates that we predict from an ablation model (RACMO). These models reproduce general spectral features in our data, and enable us to estimate what instrument exposure reduces ice-to-sensor coupling enough to prevent icequake detection. Last, we use our resonance data to quantitatively measure how narrowband signals that originate from either artificial or glaciogenic sources will reduce the ability of certain waveform detectors (correlators) to capture transient seismic events, even if sensors remain coupled.

The ability to detect the surface expression of moving water beneath the Antarctic ice sheet by satellite has revealed a dynamic basal environment, with implications for regional ice dynamics, grounding-line stability, and fluxes of freshwater and nutrients to the Southern Ocean. Knowledge of subglacial activity on timescales important for near-term prediction of ice-sheet fluctuations (decadal to century) is limited by the short observational record of NASA's Ice, Cloud, and land Elevation Satellite (ICESat) laser altimetry mission used to generate the last continent-wide survey (2003–08). Here, we use synthetic aperture radar-interferometric-mode data from ESA's CryoSat-2 radar altimetry mission (2010–present), which samples 45 of the ICESat-derived subglacial lakes, to extend their time series to the end of 2016. The extended time series show that there have been surface-height changes at 20 of the 45 lakes since 2008, indicating that some of these features are persistent and potentially cyclic, while other features show negligible changes, suggesting these may be transient or nonhydrological features. Continued monitoring of active lakes for both height and velocity changes, as well as developing methods for identifying additional lakes, is critical to quantifying the full distribution of active subglacial lakes in Antarctica.

Subglacial 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.