Search

Reliable knowledge of ice discharge dynamics for the Greenland ice sheet via its ice streams is essential if we are to understand its stability under future climate scenarios as well as their dynamics in the past, especially when using numerical models for diagnosis and prediction. Currently active ice streams in Greenland have been well mapped using remote-sensing data while past ice-stream paths in what are now deglaciated regions can be reconstructed from the landforms they left behind. However, little is known about possible former and now defunct ice streams in areas still covered by ice. Here we use radio-echo sounding data to decipher the regional ice-flow history of the northeastern Greenland ice sheet on the basis of its internal stratigraphy. By creating a three-dimensional reconstruction of time-equivalent horizons, we map folds deep below the surface that we then attribute to the deformation caused by now-extinct ice streams. We propose that locally this ancient ice-!ow regime was much more focused and reached much farther inland than today’s and was deactivated when the main drainage system was reconfigured and relocated southwards. The insight that major ice streams in Greenland might start, shift or abruptly disappear will affect our approaches to understanding and modelling the past or future response of Earth’s ice sheets to global warming. Such behaviour has to be sufficiently reproduced by numerical models operating on the mid- to longer-term timescales to be considered adequate physical representations of the naturally occuring dynamic behaviour of ice streams.

We present a high-resolution airborne radar data set (EGRIP-NOR-2018) for the onset region of the Northeast Greenland Ice Stream (NEGIS). The radar data were acquired in May 2018 with the Alfred Wegener Institute's multichannel ultra-wideband (UWB) radar mounted on the Polar 6 aircraft. Radar profiles cover an area of similar to 24 000 km(2) and extend over the well-defined shear margins of the NEGIS. The survey area is centered at the location of the drill site of the East Greenland Ice-Core Project (EastGRIP), and several radar lines intersect at this location. The survey layout was designed to (i) map the stratigraphic signature of the shear margins with radar profiles aligned perpendicular to ice flow, (ii) trace the radar stratigraphy along several flow lines, and (iii) provide spatial coverage of ice thickness and basal properties. While we are able to resolve radar reflections in the deep stratigraphy, we cannot fully resolve the steeply inclined reflections at the tightly folded shear margins in the lower part of the ice column. The NEGIS is causing the most significant discrepancies between numerically modeled and observed ice surface velocities. Given the high likelihood of future climate and ocean warming, this extensive data set of new high-resolution radar data in combination with the EastGRIP ice core will be a key contribution to understand the past and future dynamics of the NEGIS. The EGRIP-NOR-2018 radar data products can be obtained from the PANGAEA data publisher (https://doi.pangaea.de/10.1594/PANGAEA. 928569; Franke et al., 2021a).

In the text, the citation that reads "Felix and King, 2011"should read "Ng and King, 2011". Likewise, the reference that reads "Felix N and King EC (2011) Kinematic waves in polar firn stratigraphy. Journal of Glaciology 57(206), 1119–1134. doi: 10.3189/002214311798843340.” Should read “Ng F and King EC (2011) Kinematic waves in polar firn stratigraphy. Journal of Glaciology 57(206), 1119–1134. doi: 10.3189/002214311798843340".

Radar sounding is a powerful geophysical approach for characterizing the subsurface conditions of terrestrial and planetary ice masses at local to global scales. As a result, a wide array of orbital, airborne, ground-based, and in situ instruments, platforms and data analysis approaches for radioglaciology have been developed, applied or proposed. Terrestrially, airborne radar sounding has been used in glaciology to observe ice thickness, basal topography and englacial layers for five decades. More recently, radar sounding data have also been exploited to estimate the extent and configuration of subglacial water, the geometry of subglacial bedforms and the subglacial and englacial thermal states of ice sheets. Planetary radar sounders have observed, or are planned to observe, the subsurfaces and near-surfaces of Mars, Earth's Moon, comets and the icy moons of Jupiter. In this review paper, and the thematic issue of the Annals of Glaciology on ‘Five decades of radioglaciology’ to which it belongs, we present recent advances in the fields of radar systems, missions, signal processing, data analysis, modeling and scientific interpretation. Our review presents progress in these fields since the last radio-glaciological Annals of Glaciology issue of 2014, the context of their history and future prospects.

Ice shelves play a critical role in the long-term stability of ice sheets through their buttressing effect. The underlying bathymetry and cavity thickness are key inputs for modelling future ice sheet evolution. However, direct observation of sub-ice-shelf bathymetry is time-consuming, logistically risky, and in some areas simply not possible. Here we use new compilations of airborne and marine gravity, radar depth sounding, and swath bathymetry to provide new estimates of sub-ice-shelf bathymetry outboard of the rapidly changing West Antarctic Thwaites Glacier and beneath the adjacent Dotson and Crosson ice shelves. This region is of special interest, as the low-lying inland reverse slope of the Thwaites Glacier system makes it vulnerable to marine ice sheet instability, with rapid grounding line retreat observed since 1993 suggesting this process may be underway. Our results confirm a major marine channel >800 m deep extends tens of kilometres to the front of Thwaites Glacier, while the adjacent ice shelves are underlain by more complex bathymetry. Comparison of our new bathymetry with ice shelf draft reveals that ice shelves formed since 1993 comprise a distinct population where the draft conforms closely to the underlying bathymetry, unlike the older ice shelves, which show a more uniform depth of the ice base. This indicates that despite rapid basal melting in some areas, these recently floated parts of the ice shelf are not yet in dynamic equilibrium with their retreated grounding line positions and the underlying ocean system, a factor which must be included in future models of this region's evolution.

We report the discovery of a large impact crater beneath Hiawatha Glacier in northwest Greenland. From airborne radar surveys, we identify a 31-kilometer-wide, circular bedrock depression beneath up to a kilometer of ice. This depression has an elevated rim that cross-cuts tributary subglacial channels and a subdued central uplift that appears to be actively eroding. From ground investigations of the deglaciated foreland, we identify overprinted structures within Precambrian bedrock along the ice margin that strike tangent to the subglacial rim. Glaciofluvial sediment from the largest river draining the crater contains shocked quartz and other impact-related grains. Geochemical analysis of this sediment indicates that the impactor was a fractionated iron asteroid, which must have been more than a kilometer wide to produce the identified crater. Radiostratigraphy of the ice in the crater shows that the Holocene ice is continuous and conformable, but all deeper and older ice appears to be debris rich or heavily disturbed. The age of this impact crater is presently unknown, but from our geological and geophysical evidence, we conclude that it is unlikely to predate the Pleistocene inception of the Greenland Ice Sheet.

The catchments of Pine Island Glacier and Thwaites Glacier in the Amundsen Sea Embayment are two of the largest, most rapidly changing, and potentially unstable sectors of the West Antarctic Ice Sheet. They are also neighboring outlets, separated by the topographically unconfined eastern shear margin of Thwaites Glacier and the southwest tributary of Pine Island Glacier. This tributary begins just downstream of the eastern shear margin and flows into the Pine Island ice shelf. As a result, it is a potential locus of interaction between the two glaciers and could result in cross-catchment feedback during the retreat of either. Here, we analyze relative basal reflectivity profiles from three radar sounding survey lines collected using the UTIG HiCARS radar system in 2004 and CReSIS MCoRDS radar system in 2012 and 2014 to investigate the extent and character of ocean access beneath the southwest tributary. These profiles provide evidence of ocean access ~12 km inland of the 1992–2011 InSAR-derived grounding line by 2014, suggesting either retreat since 2011 or the intrusion of ocean water kilometers inland of the grounding line.

We report the discovery of a large impact crater beneath Hiawatha Glacier in northwest Greenland. From airborne radar surveys, we identify a 31-kilometer-wide, circular bedrock depression beneath up to a kilometer of ice. This depression has an elevated rim that cross-cuts tributary subglacial channels and a subdued central uplift that appears to be actively eroding. From ground investigations of the deglaciated foreland, we identify overprinted structures within Precambrian bedrock along the ice margin that strike tangent to the subglacial rim. Glaciofluvial sediment from the largest river draining the crater contains shocked quartz and other impact-related grains. Geochemical analysis of this sediment indicates that the impactor was a fractionated iron asteroid, which must have been more than a kilometer wide to produce the identified crater. Radiostratigraphy of the ice in the crater shows that the Holocene ice is continuous and conformable, but all deeper and older ice appears to be debris rich or heavily disturbed. The age of this impact crater is presently unknown, but from our geological and geophysical evidence, we conclude that it is unlikely to predate the Pleistocene inception of the Greenland Ice Sheet.

We report the discovery of a large impact crater beneath Hiawatha Glacier in northwest Greenland. From airborne radar surveys, we identify a 31-kilometer-wide, circular bedrock depression beneath up to a kilometer of ice. This depression has an elevated rim that cross-cuts tributary subglacial channels and a subdued central uplift that appears to be actively eroding. From ground investigations of the deglaciated foreland, we identify overprinted structures within Precambrian bedrock along the ice margin that strike tangent to the subglacial rim. Glaciofluvial sediment from the largest river draining the crater contains shocked quartz and other impact-related grains. Geochemical analysis of this sediment indicates that the impactor was a fractionated iron asteroid, which must have been more than a kilometer wide to produce the identified crater. Radiostratigraphy of the ice in the crater shows that the Holocene ice is continuous and conformable, but all deeper and older ice appears to be debris rich or heavily disturbed. The age of this impact crater is presently unknown, but from our geological and geophysical evidence, we conclude that it is unlikely to predate the Pleistocene inception of the Greenland Ice Sheet.

We present a compilation of radio-echo sounding (RES) measurements of five radar systems (AWI, BAS, CReSIS, INGV and UTIG) around the EPICA Dome C (EDC) drill site, East Antarctica. The aim of our study is to investigate the differences of the various systems in their resolution of internal reflection horizons (IRHs) and bed topography, penetration depth, and capacity of imaging the basal layer.We address the questions of the compatibility of existing radar data for common interpretation, and the suitability of the individual systems for reconnaissance surveys. We find that the most distinct IRHs and IRH patterns can be identified and transferred between most data sets. Considerable differences between the RES systems exist in range resolution and depiction of the bottom-most region. Considering both aspects, which we judge as crucial factors in the search for old ice, the CReSIS and the UTIG systems are the most suitable ones. In addition to the RES data set comparison we calculate a synthetic radar trace from EDC density and conductivity profiles. We identify ten common IRHs in the measured RES data and the synthetic trace. Then we conduct a sensitivity study for which we remove certain peaks from the input conductivity profile. As a result the respective reflections disappear from the modeled radar trace. In this way, we establish a depth conversion of the measured travel-times of the IRHs. Furthermore, we use these sensitivity studies to investigate the cause of observed reflections. The identified IRHs are assigned ages from the EDC’s time scale. Due to the isochronous character of these conductivity-caused IRHs, they are a means to extend the Dome C age structure by tracing the IRHs along the RES profiles.

The evaluation of ice sheet models is one of the pressing problems in the study of ice sheets dynamics. Here we examine the question of how much isochronous information is contained within the publicly available CReSIS Greenland airborne radio-echo soundings dataset. We identify regions containing isochronous reflectors using ARESP algorithms [Sime et al., 2011]. We find that isochronous reflectors are present within 36% of the CReSIS RES englacial data by location, and 41% by total number of data. Between 1000 and 3000 m in depth, isochronous reflectors are present along more than 50% of the dataset flight path. Lower volumes of cold glacial period ice also correspond with more isochronous reflectors. We find good agreement between ARESP and continuity index [Karlsson et al., 2012] results, providing confidence in these findings. Ice structure datasets, based on data identified here, will be of use in evaluating ice sheet simulations and the assessment of past rates of snow accumulation.