Subglacial Hydrological System Research Articles

Basal conditions of Denman Glacier from glacier hydrology and ice dynamics modeling

Abstract. Basal sliding in Antarctic glaciers is often modeled using a friction law that relates basal shear stresses to the effective pressure. As few ice sheet models are dynamically coupled to subglacial hydrology models, variability in subglacial hydrology associated with the effective pressure is often implicitly captured in the basal friction coefficient – an unknown parameter in the basal friction law. We investigate the impact of using effective pressures calculated from the Glacier Drainage System (GlaDS) model on basal friction coefficients calculated using inverse methods in the Ice-sheet and Sea-level System Model (ISSM) at Denman Glacier, East Antarctica, for the Schoof and Budd friction laws. For the Schoof friction law, a positive correlation emerges between the GlaDS effective pressure and basal friction coefficient in regions of fast ice flow. Using GlaDS effective pressures generally leads to smoother basal friction coefficients and basal shear stresses, and larger differences between the simulated and observed ice surface velocities compared with using an effective pressure equal to the ice overburden pressure plus the gravitational potential energy of the water. Compared with the Budd friction law, the Schoof friction law offers improved capabilities in capturing the spatial variations associated with known physics of the subglacial hydrology. Our results indicate that ice sheet model representation of basal sliding is more realistic when using direct outputs from a subglacial hydrology model, demonstrating the importance of coupling between ice sheet and subglacial hydrological systems. However, using our outputs we have also developed an empirical parameterization of effective pressure that improves the application of the Schoof friction law without requiring explicit hydrological modeling.

A parallel implementation of the confined–unconfined aquifer system model for subglacial hydrology: design, verification, and performance analysis (CUAS-MPI v0.1.0)

Abstract. The subglacial hydrological system affects (i) the motion of ice sheets through sliding, (ii) the location of lakes at the ice margin, and (iii) the ocean circulation by freshwater discharge directly at the grounding line or (iv) via rivers flowing over land. For modeling this hydrology system, a previously developed porous-media concept called the confined–unconfined aquifer system (CUAS) is used. To allow for realistic simulations at the ice sheet scale, we developed CUAS-MPI, an MPI-parallel C/C++ implementation of CUAS (MPI: Message Passing Interface), which employs the Portable, Extensible Toolkit for Scientific Computation (PETSc) infrastructure for handling grids and equation systems. We validate the accuracy of the numerical results by comparing them with a set of analytical solutions to the model equations, which involve two types of boundary conditions. We then investigate the scaling behavior of CUAS-MPI and show that CUAS-MPI scales up to 3840 MPI processes running a realistic Greenland setup on the Lichtenberg HPC system. Our measurements also show that CUAS-MPI reaches a throughput comparable to that of ice sheet simulations, e.g., the Ice-sheet and Sea-level System Model (ISSM). Lastly, we discuss opportunities for ice sheet modeling, explore future coupling possibilities of CUAS-MPI with other simulations, and consider throughput bottlenecks and limits of further scaling.

Constraints on the Timing and Extent of Deglacial Grounding Line Retreat in West Antarctica

AbstractProjections of Antarctica's contribution to future sea level rise are associated with significant uncertainty, in part because the observational record is too short to capture long‐term processes necessary to estimate ice mass changes over societally relevant timescales. Records of grounding line retreat from the geologic past offer an opportunity to extend our observations of these processes beyond the modern record and to gain a more comprehensive understanding of ice‐sheet change. Here, we present constraints on the timing and inland extent of deglacial grounding line retreat in the southern Ross Sea, Antarctica, obtained via direct sampling of a subglacial lake located 150 km inland from the modern grounding line and beneath >1 km of ice. Isotopic measurements of water and sediment from the lake enabled us to evaluate how the subglacial microbial community accessed radiocarbon‐bearing organic carbon for energy, as well as where it transferred carbon metabolically. Using radiocarbon as a natural tracer, we found that sedimentary organic carbon was microbially translocated to dissolved carbon pools in the subglacial hydrologic system during the 4.5‐year period of water accumulation prior to our sampling. This finding indicates that the grounding line along the Siple Coast of West Antarctica retreated more than 250 km inland during the mid‐Holocene (6.3 ± 1.0 ka), prior to re‐advancing to its modern position.

Outlet glacier flow response to surface melt: based on analysis of a high-resolution satellite data set

The dynamics of the Greenland Ice Sheet are affected by surface meltwater reaching the base of the ice, altering ice contact with the bedrock. Lack of understanding of this evolution hampers the ability to predict the effects of increasing temperatures on the Greenland Ice Sheet mass balance. Here we present a unique high-resolution study of ice velocity response to surface melting based on data from a COSMO-SkyMed satellite campaign over Upernavik Isstrøm (Northwest Greenland) for two months around the end of the 2014 melt season. We show that the velocity variations, due to both short-term (days) and seasonal variations in surface melt rates, are increasing in relative strength farther from the glacier terminus. Furthermore, we observe how ice dynamic response to frontal retreat, reaching several kilometres inland, can obscure the meltwater-induced velocity change close to the terminus. Future studies should consider the flow velocity dependence on the distance to the terminus, and local geometry, to distinguish subglacial hydrologic system changes from frontal processes and local basal conditions.

Tidewater-glacier response to supraglacial lake drainage

The flow speed of the Greenland Ice Sheet changes dramatically in inland regions when surface meltwater drains to the bed. But ice-sheet discharge to the ocean is dominated by fast-flowing outlet glaciers, where the effect of increasing surface melt on annual discharge is unknown. Observations of a supraglacial lake drainage at Helheim Glacier, and a consequent velocity pulse propagating down-glacier, provide a natural experiment for assessing the impact of changes in injected meltwater, and allow us to interrogate the subglacial hydrological system. We find a highly efficient subglacial drainage system, such that summertime lake drainage has little net effect on ice discharge. Our results question the validity of common remote-sensing approaches for inferring subglacial conditions, knowledge of which is needed for improved projections of sea-level rise.

Tunnel valley formation beneath deglaciating mid-latitude ice sheets: Observations and modelling

The geological record of landforms and sediments produced beneath deglaciating ice sheets offers insights into inaccessible glacial processes. Large subglacial valleys formed by meltwater erosion of sediments (tunnel valleys) are widespread in formerly glaciated regions such as the North Sea. Obtaining a better understanding of these features may help with the parameterisation of basal melt rates and the interplay between basal hydrology and ice dynamics in numerical models of past, present, and future ice-sheet configurations. However, the mechanisms and timescales over which tunnel valleys form remain poorly constrained. Here, we present a series of numerical modelling experiments, informed by new observations from high-resolution 3D seismic data (6.25 m bin size, ∼4 m vertical resolution), which test different hypotheses of tunnel valley formation and calculate subglacial water routing, seasonal water discharges, and the rates at which tunnel valleys are eroded beneath deglaciating ice sheets. Networks of smaller or abandoned channels, pervasive slump deposits, and subglacial landforms are imaged inside and at the base of larger tunnel valleys, indicating that these tunnel valleys were carved through the action of migrating smaller channels within tens of kilometres of the ice margin and were later widened by ice-contact erosion. Our model results imply that the drainage of extensive surface meltwater to the ice-sheet bed is the dominant mechanism responsible for tunnel valley formation; this process can drive rapid incision of networks of regularly spaced subglacial tunnel valleys beneath the fringes of retreating ice sheets within hundreds to thousands of years during deglaciation. Combined, our observations and modelling results identify how tunnel valleys form beneath deglaciating mid-latitude ice sheets and have implications for how the subglacial hydrological systems of contemporary ice sheets may respond to sustained climate warming.

Declining Basal Motion Dominates the Long‐Term Slowing of Athabasca Glacier, Canada

AbstractGlobally, glaciers are shrinking in response to climate change, with implications for global sea level rise as well as downstream ecosystems and water resources. Sliding at the ice‐bed interface (basal motion) provides a mechanism for glaciers to respond rapidly to climate change. While the short‐term dynamics of glacier basal motion (<10 years) have received substantial attention, little is known about how basal motion and its sensitivity to subglacial hydrology changes over long (>50 year) timescales—this knowledge is required for accurate prediction of future glacier change. We compare historical data with modern estimates from field and satellite data at Athabasca Glacier and show that the glacier thinned by 60 m (−21%) over 1961–2020. However, a concurrent increase in surface slope results in minimal change in the average driving stress (−6 kPa and −4%). These geometric changes coincide with relatively uniform slowing (−15 m a−1 and −45%). Simplified ice modeling suggests that declining basal motion accounts for most of this slow down (91% on average and 46% at minimum). A decline in basal motion can be explained by increasing basal friction resulting from geometric change in addition to increasing meltwater flux through a more efficient subglacial hydrologic system. These results highlight the need to include time‐varying dynamics of basal motion in glacier models and analyses. If these findings are generalizable, they suggest that declining basal motion reduces the flux of ice to lower elevations, helping to mitigate glacier mass loss in a warming climate.

Sedimentary Signatures of Persistent Subglacial Meltwater Drainage From Thwaites Glacier, Antarctica

Subglacial meltwater drainage can enhance localized melting along grounding zones and beneath the ice shelves of marine-terminating glaciers. Efforts to constrain the evolution of subglacial hydrology and the resulting influence on ice stability in space and on decadal to millennial timescales are lacking. Here, we apply sedimentological, geochemical, and statistical methods to analyze sediment cores recovered offshore Thwaites Glacier, West Antarctica to reconstruct meltwater drainage activity through the pre-satellite era. We find evidence for a long-lived subglacial hydrologic system beneath Thwaites Glacier and indications that meltwater plumes are the primary mechanism of sedimentation seaward of the glacier today. Detailed core stratigraphy revealed through computed tomography scanning captures variability in drainage styles and suggests greater magnitudes of sediment-laden meltwater have been delivered to the ocean in recent centuries compared to the past several thousand years. Fundamental similarities between meltwater plume deposits offshore Thwaites Glacier and those described in association with other Antarctic glacial systems imply widespread and similar subglacial hydrologic processes that occur independently of subglacial geology. In the context of Holocene changes to the Thwaites Glacier margin, it is likely that subglacial drainage enhanced submarine melt along the grounding zone and amplified ice-shelf melt driven by oceanic processes, consistent with observations of other West Antarctic glaciers today. This study highlights the necessity of accounting for the influence of subglacial hydrology on grounding-zone and ice-shelf melt in projections of future behavior of the Thwaites Glacier ice margin and marine-based glaciers around the Antarctic continent.

Decadal Changes in Greenland Ice Sheet Firn Aquifers from Radar Scatterometer

Surface meltwater runoff is believed to be the main cause of the alarming mass loss in the Greenland Ice Sheet (GrIS); however, recent research has shown that a large amount of meltwater is not directly drained or refrozen but stored in the form of firn aquifers (FAs) in the interior of the GrIS. Monitoring the changes in FAs over the GrIS is of great importance to evaluate the stability and mass balance of the ice sheet. This is challenging because FAs are not visible on the surface and the direct measurements are lacking. A new method is proposed to map FAs during the 2010–2020 period by using the C-band Advanced Scatterometer (ASCAT) data based on the Random Forests classification algorithm with the aid of measurements from the NASA Operation IceBridge (OIB) program. Melt days (MD), melt intensity (MI), and winter mean backscatter (WM) parameters derived from the ASCAT data are used as the input vectors for the Random Forests classification algorithm. The accuracy of the classification model is assessed by ten-fold cross-validation, and the overall accuracy and Kappa coefficient are 97.49% and 0.72 respectively. The results show that FAs reached the maximum in 2015, and the accumulative area of FAs from 2010 to 2020 is 56,477 km2, which is 3.3% of the GrIS area. This study provides a way to investigate the long-term dynamics in FAs which have great significance for understanding the state of subsurface firn and subglacial hydrological systems.

Seismic Mapping of Subglacial Hydrology Reveals Previously Undetected Pressurization Event

AbstractUnderstanding the dynamic response of glaciers to climate change is vital for assessing water resources and hazards, and subglacial hydrology is a key player in glacier systems. Traditional observations of subglacial hydrology are spatially and temporally limited, but recent seismic deployments on and around glaciers show the potential for comprehensive observation of glacial hydrologic systems. We present results from a high‐density seismic deployment spanning the surface of Lemon Creek Glacier, Alaska. Our study coincided with a marginal lake drainage event, which served as a natural experiment for seismic detection of changes in subglacial hydrology. We observed glaciohydraulic tremor across the surface of the glacier that was generated by the subglacial hydrologic system. During the lake drainage, the relative changes in seismic tremor power and water flux are consistent with pressurization of the subglacial system of only the upper part of the glacier. This event was not accompanied by a significant increase in glacier velocity; either some threshold necessary for rapid basal motion was not attained, or, plausibly, the geometry of Lemon Creek Glacier inhibited speedup. This pressurization event would have likely gone undetected without seismic observations, demonstrating the power of cryoseismology in testing assumptions about and mapping the spatial extent of subglacial pressurization.

Monitoring the Hydrological Activities of Antarctic Subglacial Lakes Using CryoSat-2 and ICESat-2 Altimetry Data

Monitoring the hydrological activities of subglacial lakes is critical to understanding the subglacial hydrological system and evaluating the internal mass changes of the Antarctic ice sheet. Drainage or filling events of active lakes lead to elevation changes in the ice surface. These changes can be observed by satellite altimetry, but the monitoring must be conducted continuously since the water movements in active subglacial lakes may occur frequently. We used CryoSat-2 Baseline-D and ICESat-2 data from 2010 to 2020 to obtain the time series of the ice surface elevation changes for 17 active lakes. We also evaluated the uncertainty of the time series derived from the CryoSat-2 data by cross-validation. The mean and RMS of the biases between the CryoSat-2-based and ICESat-2-based time series are generally less than 0.3 m and 1.0 m, respectively. However, the mean and RMS are greater over the lakes with rough ice surfaces, such as Whillans6, KT1, Mac3, and Slessor23. The drainage and filling events continue exhibiting in the extended period and the hydrological activities of SLW, L12, Whillans6, L78, and Mac1 occurred periodically. Furthermore, we inferred the hydrological connections between the lakes combining simulated water pathways.

Subglacial Water Flow Over an Antarctic Palaeo‐Ice Stream Bed

AbstractThe subglacial hydrological system exerts a critical control on the dynamic behavior of the overlying ice because its configuration affects the degree of basal lubrication between the ice and the bed. Yet, this component of the glaciological system is notoriously hard to access and observe, particularly over timescales longer than the satellite era. In Antarctica, abundant evidence for past subglacial water flow over former ice‐sheet beds exists around the peripheries of the ice sheet including networks of huge channels carved into bedrock (now submarine) on the Pacific margin of West Antarctica. Here, we combine detailed bathymetric investigations of a channel system in Marguerite Trough, a major palaeo‐ice stream bed, with numerical hydrological modeling to explore subglacial water accumulation, routing and potential for erosion over decadal‐centennial timescales. Detailed channel morphologies from remotely operated vehicle surveys indicate multiple stages of localized incision, and the occurrence of potholes, some gigantic in scale, suggests incision by turbulent water carrying a significant bedload. Further, the modeling indicates that subglacial water is available during deglaciation and was likely released in episodic drainage events, from subglacial lakes, varying in magnitude over time. Our observations support previous assertions that these huge bedrock channel systems were incised over multiple glacial cycles through episodic subglacial lake drainage events; however, here we present a viable pattern for subglacial drainage at times when the ice sheet existed over the continental shelf and was capable of continuing to erode the bedrock substrate.

Long-period variability in ice-dammed glacier outburst floods due to evolving catchment geometry

Abstract. We combine a glacier outburst flood model with a glacier flow model to investigate decadal to centennial variations in outburst floods originating from ice-dammed marginal basins. Marginal basins can form due to the retreat and detachment of tributary glaciers, a process that often results in remnant ice being left behind. The remnant ice, which can act like an ice shelf or break apart into a pack of icebergs, limits a basin's water storage capacity but also exerts pressure on the underlying water and promotes drainage. We find that during glacier retreat there is a strong, nearly linear relationship between flood water volume and peak discharge for individual basins, despite large changes in glacier and remnant ice volumes that are expected to impact flood hydrographs. Consequently, peak discharge increases over time as long as there is remnant ice remaining in a basin, and peak discharge begins to decrease once a basin becomes ice-free. Thus, similar size outburst floods can occur at very different stages of glacier retreat. We also find that the temporal variability in outburst flood magnitude depends on how the floods initiate. Basins that connect to the subglacial hydrological system only after reaching flotation depth yield greater long-term variability in outburst floods than basins that are continuously connected to the subglacial hydrological system (and therefore release floods that initiate before reaching flotation depth). Our results highlight the importance of improving our understanding of both changes in basin geometry and outburst flood initiation mechanisms in order to better assess outburst flood hazards and their impacts on landscape and ecosystem evolution.

Proper orthogonal decomposition of ice velocity identifies drivers of flow variability at Sermeq Kujalleq (Jakobshavn Isbræ)

Abstract. The increasing volume and spatio-temporal resolution of satellite-derived ice velocity data have created new exploratory opportunities for the quantitative analysis of glacier dynamics. One potential technique, proper orthogonal decomposition (POD), also known as empirical orthogonal functions, has proven to be a powerful and flexible technique for revealing coherent structures in a wide variety of environmental flows. In this study we investigate the applicability of POD to an openly available TanDEM-X/TerraSAR-X-derived ice velocity dataset from Sermeq Kujalleq (Jakobshavn Isbræ), Greenland. We find three dominant modes with annual periodicity that we argue are explained by glaciological processes. The primary dominant mode is interpreted as relating to the stress reconfiguration at the glacier terminus, known to be an important control on the glacier's dynamics. The second and third largest modes together relate to the development of the spatially heterogenous glacier hydrological system and are primarily driven by the pressurisation and efficiency of the subglacial hydrological system. During the melt season, variations in the velocity shown in these two subsidiary modes are explained by the drainage of nearby supraglacial melt ponds, as identified with a Google Earth Engine Moderate Resolution Imaging Spectroradiometer (MODIS) dynamic thresholding technique. By isolating statistical structures within velocity datasets and through their comparison to glaciological theory and complementary datasets, POD indicates which glaciological processes are responsible for the changing bulk velocity signal, as observed from space. With the proliferation of optical- and radar-derived velocity products (e.g. MEaSUREs, ESA CCI, PROMICE), we suggest POD, and potentially other modal decomposition techniques, will become increasingly useful in future studies of ice dynamics.

Seasonal dispersal of fjord meltwaters as an important source of iron and manganese to coastal Antarctic phytoplankton

Abstract. Glacial meltwater from the western Antarctic Ice Sheet is hypothesized to be an important source of cryospheric iron, fertilizing the Southern Ocean, yet its trace-metal composition and factors that control its dispersal remain poorly constrained. Here we characterize meltwater iron sources in a heavily glaciated western Antarctic Peninsula (WAP) fjord. Using dissolved and particulate ratios of manganese to iron in meltwaters, porewaters, and seawater, we show that surface glacial melt and subglacial plumes contribute to the seasonal cycle of iron and manganese within a fjord still relatively unaffected by climate-change-induced glacial retreat. Organic ligands derived from the phytoplankton bloom and the glaciers bind dissolved iron and facilitate the solubilization of particulate iron downstream. Using a numerical model, we show that buoyant plumes generated by outflow from the subglacial hydrologic system, enriched in labile particulate trace metals derived from a chemically modified crustal source, can supply iron to the fjord euphotic zone through vertical mixing. We also show that prolonged katabatic wind events enhance export of meltwater out of the fjord. Thus, we identify an important atmosphere–ice–ocean coupling intimately tied to coastal iron biogeochemistry and primary productivity along the WAP.

Surface Expression and Apparent Timing of Subglacial Lake Oscillations Controlled by Viscous Ice Flow

AbstractSubglacial lakes are bodies of water that form at the base of ice sheets and glaciers. Ice‐surface elevation above these lakes responds to water volume changes, providing one of few ways to monitor subglacial hydrological systems. Here, we use numerical models to compare surface elevation‐derived estimates of subglacial‐lake length, water‐volume change, and highstand or lowstand timing with their true values. Because ice flow influences the surface expression of lake‐volume change, the correspondence between these estimates and their true values depends strongly on ice thickness, volume‐change rate, and basal drag coefficient. For many realistic combinations of these factors, viscous relaxation of the ice‐sheet surface can render lake volume‐changes unobservable with altimetry. Our results highlight the need for new estimation methods that account for the effects of ice flow, as well as improvements to current resolution limitations that render some events unobservable with altimetry.

The distribution of glacial meltwater routes and associated murtoo fields in Finland

High-resolution LiDAR (Light detection and Ranging) -based digital elevation models (DEM) have greatly improved the mapping of glacial landforms and revealed new ones such as murtoos. Murtoos have extensive diversity in form, relief and size, and they often appear along meltwater routes. However, not all meltwater routes in the recently glaciated terrains include murtoos. We mapped different types of subglacial meltwater routes and the related distribution of murtoos in the Finnish part of Fennoscandian Ice Sheet (FIS). Subglacial meltwater routes represent a previously unknown extension of the subglacial hydrological system that supplements esker networks. Murtoo deposition along the routes is dictated by the marked concentration and routing of subglacial meltwater in a high-pressure environment outside the subglacial tunnel flow and channelized drainage zone. The main environments of murtoo route genesis include the margins of glacial lineation fields, lateral shear margins of ice streams or in between ice-flow sectors or corridors, confluence zones of ice stream onset areas, lee-sides of bedrock protrusions or thresholds, bedrock fracture valleys and potential subglacial lake inputs and outputs. This study adds to the knowledge of the development of subglacial drainage and emphasizes the role of murtoos as the potential missing link between the channelized and distributed subglacial drainage.

A Late Pleistocene channelized subglacial meltwater system on the Atlantic continental shelf south of Ireland

The study of palaeo‐glacial landforms and sediments can give insights into the nature and dynamics of ice sheets. This is particularly the case with regards to the subglacial record, which is challenging to observe in contemporary glaciated settings and hence remains only partially understood. The subglacial hydrological system is an essential component of ice dynamics, where increased water pressure enhances ice motion and sediment deformation, thus reducing ice‐bed contact. Tunnel valleys are large, sinuous, steep‐sided incisions that, together with smaller scale meltwater channels, indicate subglacial meltwater discharge beneath large ice sheets. Through the use of high‐resolution marine geophysical data, a system of buried and exposed tunnel valleys, possible subglacial or proglacial meltwater channels and palaeo‐fluvial valleys have been identified across the shelf of the Celtic Sea between Ireland and Britain. The presence of steep‐sided and overdeepened tunnel valleys is indicative of a large channelized meltwater drainage system beneath the former Irish Sea Ice Stream, the most extensive ice stream to drain the last British–Irish Ice Sheet. After the rapid ice expansion across the Celtic Sea shelf around 28–26 ka, the tunnel valleys were carved into both bedrock and glacigenic sediments and are associated with rapid ice stream retreat northwards into the Irish Sea Basin between 25.6 and 24.3 ka. The presence of a major subglacial meltwater system on the relatively shallow shelf suggests that significant erosive meltwater discharge occurred during the last deglaciation and highlights the important contribution of meltwater to the retreat of the British–Irish Ice Sheet on the continental shelf.

Variations in esker morphology and internal architecture record time-transgressive deposition during ice margin retreat in Northern Ireland

The architecture and evolution of the subglacial hydrological system plays a key role in modulating ice flow. Eskers provide an opportunity to understand subglacial hydrology at a broader perspective than contemporary studies. Recent research has established a morphogenetic classification for eskers, but these studies have been limited to topographically simple regions of a single ice sheet. We present an updated map of esker distribution in Northern Ireland based on 5-m resolution elevation data. We also present a high-resolution map of the glacial geomorphology of SW Northern Ireland, based on ~0.4-m resolution elevation data. Ground Penetrating Radar data from four sites along the >20-km long Evishanoran Esker system in central Northern Ireland are combined with geomorphological observations to provide insight into depositional processes and controls on esker formation. Esker architecture indicates two styles of deposition, including an initial high energy flow event in a subglacial conduit and delta foreset deposition close to the ice sheet margin during ice margin retreat. These delta foreset deposits can be used to reconstruct former ice margins. We identify that local topographic complexity and geological structures (e.g., faults) are important controls on esker formation. The broad-scale esker architecture remains the same despite variable esker planform morphology, suggesting hydrological conditions alone cannot explain esker morphology. This study provides further evidence that morphogenetic relationships cannot be based solely on remote sensing data and must be supported by robust field observations, especially where post-glacial processes may distort esker morphology (e.g., peat infilling).

Cryoegg: development and field trials of a wireless subglacial probe for deep, fast-moving ice

AbstractSubglacial hydrological systems require innovative technological solutions to access and observe. Wireless sensor platforms can be used to collect and return data, but their performance in deep and fast-moving ice requires quantification. We report experimental results from Cryoegg: a spherical probe that can be deployed into a borehole or moulin and transit through the subglacial hydrological system. The probe measures temperature, pressure and electrical conductivity in situ and returns all data wirelessly via a radio link. We demonstrate Cryoegg's utility in studying englacial channels and moulins, including in situ salt dilution gauging. Cryoegg uses VHF radio to transmit data to a surface receiving array. We demonstrate transmission through up to 1.3 km of cold ice – a significant improvement on the previous design. The wireless transmission uses Wireless M-Bus on 169 MHz; we present a simple radio link budget model for its performance in cold ice and experimentally confirm its validity. Cryoegg has also been tested successfully in temperate ice. The battery capacity should allow measurements to be made every 2 h for more than a year. Future iterations of the radio system will enable Cryoegg to transmit data through up to 2.5 km of ice.