Multidecadal observations of the Antarctic ice sheet from restored analog radar records.

Airborne radar sounding can measure conditions within and beneath polar ice sheets. In Antarctica, most digital radar-sounding data have been collected in the last 2 decades, limiting our ability to understand processes that govern longer-term ice-sheet behavior. Here, we demonstrate how analog radar data collected over 40 y ago in Antarctica can be combined with modern records to quantify multidecadal changes. Specifically, we digitize over 400,000 line kilometers of exploratory Antarctic radar data originally recorded on 35-mm optical film between 1971 and 1979. We leverage the increased geometric and radiometric resolution of our digitization process to show how these data can be used to identify and investigate hydrologic, geologic, and topographic features beneath and within the ice sheet. To highlight their scientific potential, we compare the digitized data with contemporary radar measurements to reveal that the remnant eastern ice shelf of Thwaites Glacier in West Antarctica had thinned between 10 and 33% between 1978 and 2009. We also release the collection of scanned radargrams in their entirety in a persistent public archive along with updated geolocation data for a subset of the data that reduces the mean positioning error from 5 to 2.5 km. Together, these data represent a unique and renewed extensive, multidecadal historical baseline, critical for observing and modeling ice-sheet change on societally relevant timescales.

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.

A dynamic early East Antarctic Ice Sheet suggested by ice-covered fjord landscapes.

The first Cenozoic ice sheets initiated in Antarctica from the Gamburtsev Subglacial Mountains and other highlands as a result of rapid global cooling ∼34 million years ago. In the subsequent 20 million years, at a time of declining atmospheric carbon dioxide concentrations and an evolving Antarctic circumpolar current, sedimentary sequence interpretation and numerical modelling suggest that cyclical periods of ice-sheet expansion to the continental margin, followed by retreat to the subglacial highlands, occurred up to thirty times. These fluctuations were paced by orbital changes and were a major influence on global sea levels. Ice-sheet models show that the nature of such oscillations is critically dependent on the pattern and extent of Antarctic topographic lowlands. Here we show that the basal topography of the Aurora Subglacial Basin of East Antarctica, at present overlain by 2-4.5 km of ice, is characterized by a series of well-defined topographic channels within a mountain block landscape. The identification of this fjord landscape, based on new data from ice-penetrating radar, provides an improved understanding of the topography of the Aurora Subglacial Basin and its surroundings, and reveals a complex surface sculpted by a succession of ice-sheet configurations substantially different from today's. At different stages during its fluctuations, the edge of the East Antarctic Ice Sheet lay pinned along the margins of the Aurora Subglacial Basin, the upland boundaries of which are currently above sea level and the deepest parts of which are more than 1 km below sea level. Although the timing of the channel incision remains uncertain, our results suggest that the fjord landscape was carved by at least two iceflow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.

Evidence for elevated and spatially variable geothermal flux beneath the West Antarctic Ice Sheet.

Heterogeneous hydrologic, lithologic, and geologic basal boundary conditions can exert strong control on the evolution, stability, and sea level contribution of marine ice sheets. Geothermal flux is one of the most dynamically critical ice sheet boundary conditions but is extremely difficult to constrain at the scale required to understand and predict the behavior of rapidly changing glaciers. This lack of observational constraint on geothermal flux is particularly problematic for the glacier catchments of the West Antarctic Ice Sheet within the low topography of the West Antarctic Rift System where geothermal fluxes are expected to be high, heterogeneous, and possibly transient. We use airborne radar sounding data with a subglacial water routing model to estimate the distribution of basal melting and geothermal flux beneath Thwaites Glacier, West Antarctica. We show that the Thwaites Glacier catchment has a minimum average geothermal flux of ∼ 114 ± 10 mW/m(2) with areas of high flux exceeding 200 mW/m(2) consistent with hypothesized rift-associated magmatic migration and volcanism. These areas of highest geothermal flux include the westernmost tributary of Thwaites Glacier adjacent to the subaerial Mount Takahe volcano and the upper reaches of the central tributary near the West Antarctic Ice Sheet Divide ice core drilling site.

Evidence for a water system transition beneath Thwaites Glacier, West Antarctica.

Thwaites Glacier is one of the largest, most rapidly changing glaciers on Earth, and its landward-sloping bed reaches the interior of the marine West Antarctic Ice Sheet, which impounds enough ice to yield meters of sea-level rise. Marine ice sheets with landward-sloping beds have a potentially unstable configuration in which acceleration can initiate or modulate grounding-line retreat and ice loss. Subglacial water has been observed and theorized to accelerate the flow of overlying ice dependent on whether it is hydrologically distributed or concentrated. However, the subglacial water systems of Thwaites Glacier and their control on ice flow have not been characterized by geophysical analysis. The only practical means of observing these water systems is airborne ice-penetrating radar, but existing radar analysis approaches cannot discriminate between their dynamically critical states. We use the angular distribution of energy in radar bed echoes to characterize both the extent and hydrologic state of subglacial water systems across Thwaites Glacier. We validate this approach with radar imaging, showing that substantial water volumes are ponding in a system of distributed canals upstream of a bedrock ridge that is breached and bordered by a system of concentrated channels. The transition between these systems occurs with increasing surface slope, melt-water flux, and basal shear stress. This indicates a feedback between the subglacial water system and overlying ice dynamics, which raises the possibility that subglacial water could trigger or facilitate a grounding-line retreat in Thwaites Glacier capable of spreading into the interior of the West Antarctic Ice Sheet.

Radar for Europa Assessment and Sounding: Ocean to Near-Surface (REASON).

The Radar for Europa Assessment and Sounding: Ocean to Near-surface (REASON) is a dual-frequency ice-penetrating radar (9 and 60 MHz) onboard the Europa Clipper mission. REASON is designed to probe Europa from exosphere to subsurface ocean, contributing the third dimension to observations of this enigmatic world. The hypotheses REASON will test are that (1) the ice shell of Europa hosts liquid water, (2) the ice shell overlies an ocean and is subject to tidal flexing, and (3) the exosphere, near-surface, ice shell, and ocean participate in material exchange essential to the habitability of this moon. REASON will investigate processes governing this material exchange by characterizing the distribution of putative non-ice material (e.g., brines, salts) in the subsurface, searching for an ice-ocean interface, characterizing the ice shell's global structure, and constraining the amplitude of Europa's radial tidal deformations. REASON will accomplish these science objectives using a combination of radar measurement techniques including altimetry, reflectometry, sounding, interferometry, plasma characterization, and ranging. Building on a rich heritage from Earth, the moon, and Mars, REASON will be the first ice-penetrating radar to explore the outer solar system. Because these radars are untested for the icy worlds in the outer solar system, a novel approach to measurement quality assessment was developed to represent uncertainties in key properties of Europa that affect REASON performance and ensure robustness across a range of plausible parameters suggested for the icy moon. REASON will shed light on a never-before-seen dimension of Europa and - in concert with other instruments on Europa Clipper - help to investigate whether Europa is a habitable world.

Exploring the Interior of Europa with the Europa Clipper.

The Galileo mission to Jupiter revealed that Europa is an ocean world. The Galileo magnetometer experiment in particular provided strong evidence for a salty subsurface ocean beneath the ice shell, likely in contact with the rocky core. Within the ice shell and ocean, a number of tectonic and geodynamic processes may operate today or have operated at some point in the past, including solid ice convection, diapirism, subsumption, and interstitial lake formation. The science objectives of the Europa Clipper mission include the characterization of Europa's interior; confirmation of the presence of a subsurface ocean; identification of constraints on the depth to this ocean, and on its salinity and thickness; and determination of processes of material exchange between the surface, ice shell, and ocean. Three broad categories of investigation are planned to interrogate different aspects of the subsurface structure and properties of the ice shell and ocean: magnetic induction, subsurface radar sounding, and tidal deformation. These investigations are supplemented by several auxiliary measurements. Alone, each of these investigations will reveal unique information. Together, the synergy between these investigations will expose the secrets of the Europan interior in unprecedented detail, an essential step in evaluating the habitability of this ocean world.

Double ridge formation over shallow water sills on Jupiter's moon Europa.

Jupiter's moon Europa is a prime candidate for extraterrestrial habitability in our solar system. The surface landforms of its ice shell express the subsurface structure, dynamics, and exchange governing this potential. Double ridges are the most common surface feature on Europa and occur across every sector of the moon, but their formation is poorly understood, with current hypotheses providing competing and incomplete mechanisms for the development of their distinct morphology. Here we present the discovery and analysis of a double ridge in Northwest Greenland with the same gravity-scaled geometry as those found on Europa. Using surface elevation and radar sounding data, we show that this double ridge was formed by successive refreezing, pressurization, and fracture of a shallow water sill within the ice sheet. If the same process is responsible for Europa's double ridges, our results suggest that shallow liquid water is spatially and temporally ubiquitous across Europa's ice shell.

Ice mass loss sensitivity to the Antarctic ice sheet basal thermal state.

Sea-level rise projections rely on accurate predictions of ice mass loss from Antarctica. Climate change promotes greater mass loss by destabilizing ice shelves and accelerating the discharge of upstream grounded ice. Mass loss is further exacerbated by mechanisms such as the Marine Ice Sheet Instability and the Marine Ice Cliff Instability. However, the effect of basal thermal state changes of grounded ice remains largely unexplored. Here, we use numerical ice sheet modeling to investigate how warmer basal temperatures could affect the Antarctic ice sheet mass balance. We find increased mass loss in response to idealized basal thawing experiments run over 100 years. Most notably, frozen-bed patches could be tenuously sustaining the current ice configuration in parts of George V, Adélie, Enderby, and Kemp Land regions of East Antarctica. With less than 5 degrees of basal warming, these frozen patches may begin to thaw, producing new loci of mass loss.

Extreme melt season ice layers reduce firn permeability across Greenland.

Surface meltwater runoff dominates present-day mass loss from the Greenland Ice Sheet. In Greenland's interior, porous firn can limit runoff by retaining meltwater unless perched low-permeability horizons, such as ice slabs, develop and restrict percolation. Recent observations suggest that such horizons might develop rapidly during extreme melt seasons. Here we present radar sounding evidence that an extensive near surface melt layer formed following the extreme melt season in 2012. This layer was still present in 2017 in regions up to 700 m higher in elevation and 160 km further inland than known ice slabs. We find that melt layer formation is driven by local, short-timescale thermal and hydrologic processes in addition to mean climate state. These melt layers reduce vertical percolation pathways, and, under appropriate firn temperature and surface melt conditions, encourage further ice aggregation at their horizon. Therefore, the frequency of extreme melt seasons relative to the rate at which pore space and cold content regenerates above the most recent melt layer may be a key determinant of the firn's multi-year response to surface melt.

Multisystem Synthesis of Radar Sounding Observations of the Amundsen Sea Sector From the 2004-2005 Field Season.

The Amundsen Sea Embayment of the West Antarctic Ice Sheet contains Thwaites and Pine Island Glaciers, two of the most rapidly changing glaciers in Antarctica. To date, Pine Island and Thwaites Glaciers have only been observed by independent airborne radar sounding surveys, but a combined cross-basin analysis that investigates the basal conditions across the Pine Island-Thwaites Glaciers boundary has not been performed. Here, we combine two radar surveys and correct for their differences in system parameters to produce unified englacial attenuation and basal relative reflectivity maps spanning both Pine Island and Thwaites Glaciers. Relative reflectivities range from -24.8 to +37.4 dB with the highest values beneath fast-flowing ice at the ice sheet margin. By comparing our reflectivity results with previously derived radar specularity and trailing bed echoes at Thwaites Glacier, we find a highly diverse subglacial landscape and hydrologic conditions that evolve along-flow. Together, these findings highlight the potential for joint airborne radar analysis with ground-based seismic and geomorphological observations to understand variations in the bed properties and cross-catchment interactions of ice streams and outlet glaciers.

Discovery of a hypersaline subglacial lake complex beneath Devon Ice Cap, Canadian Arctic.

Subglacial lakes are unique environments that, despite the extreme dark and cold conditions, have been shown to host microbial life. Many subglacial lakes have been discovered beneath the ice sheets of Antarctica and Greenland, but no spatially isolated water body has been documented as hypersaline. We use radio-echo sounding measurements to identify two subglacial lakes situated in bedrock troughs near the ice divide of Devon Ice Cap, Canadian Arctic. Modeled basal ice temperatures in the lake area are no higher than -10.5°C, suggesting that these lakes consist of hypersaline water. This implication of hypersalinity is in agreement with the surrounding geology, which indicates that the subglacial lakes are situated within an evaporite-rich sediment unit containing a bedded salt sequence, which likely act as the solute source for the brine. Our results reveal the first evidence for subglacial lakes in the Canadian Arctic and the first hypersaline subglacial lakes reported to date. We conclude that these previously unknown hypersaline subglacial lakes may represent significant and largely isolated microbial habitats, and are compelling analogs for potential ice-covered brine lakes and lenses on planetary bodies across the solar system.

Rapid submarine ice melting in the grounding zones of ice shelves in West Antarctica.

Enhanced submarine ice-shelf melting strongly controls ice loss in the Amundsen Sea embayment (ASE) of West Antarctica, but its magnitude is not well known in the critical grounding zones of the ASE's major glaciers. Here we directly quantify bottom ice losses along tens of kilometres with airborne radar sounding of the Dotson and Crosson ice shelves, which buttress the rapidly changing Smith, Pope and Kohler glaciers. Melting in the grounding zones is found to be much higher than steady-state levels, removing 300-490 m of solid ice between 2002 and 2009 beneath the retreating Smith Glacier. The vigorous, unbalanced melting supports the hypothesis that a significant increase in ocean heat influx into ASE sub-ice-shelf cavities took place in the mid-2000s. The synchronous but diverse evolutions of these glaciers illustrate how combinations of oceanography and topography modulate rapid submarine melting to hasten mass loss and glacier retreat from West Antarctica.