Search

Only a few localised ice streams drain most of the ice from the Greenland Ice Sheet. Thus, understanding ice stream behaviour and its temporal variability is crucially important to predict future sea-level change. The interior trunk of the 700 km-long North-East Greenland Ice Stream (NEGIS) is remarkable due to the lack of any clear bedrock channel to explain its presence. Here, we present a 3-dimensional analysis of the folding and advection of its stratigraphic horizons, which shows that the localised flow and shear margins in the upper NEGIS were fully developed only ca 2000 years ago. Our results contradict the assumption that the ice stream has been stable throughout the Holocene in its current form and show that upper NEGIS-type development of ice streaming, with distinct shear margins and no bed topography relationship, can be established on time scales of hundreds of years, which is a major challenge for realistic mass-balance and sea-level rise projections.

The European Beyond EPICA project aims to extract a continuous ice core of up to 1.5 Ma, with a maximum age density of 20 kyr m−1 at Little Dome C (LDC). We present a 1D numerical model which calculates the age of the ice around Dome C. The model inverts for basal conditions and accounts either for melting or for a layer of stagnant ice above the bedrock. It is constrained by internal reflecting horizons traced in radargrams and dated using the EPICA Dome C (EDC) ice core age profile. We used three different radar datasets ranging from a 10 000 km2 airborne survey down to 5 km long ground-based radar transects over LDC. We find that stagnant ice exists in many places, including above the LDC relief where the new Beyond EPICA drill site (BELDC) is located. The modelled thickness of this layer of stagnant ice roughly corresponds to the thickness of the basal unit observed in one of the radar surveys and in the autonomous phase-sensitive radio-echo sounder (ApRES) dataset. At BELDC, the modelled stagnant ice thickness is 198±44 m and the modelled oldest age of ice is 1.45±0.16 Ma at a depth of 2494±30 m. This is very similar to all sites situated on the LDC relief, including that of the Million Year Ice Core project being conducted by the Australian Antarctic Division. The model was also applied to radar data in the area 10–15 km north of EDC (North Patch), where we find either a thin layer of stagnant ice (generally <60 m) or a negligible melt rate (<0.1 mm yr−1). The modelled maximum age at North Patch is over 2 Ma in most places, with ice at 1.5 Ma having a resolution of 9–12 kyr m−1, making it an exciting prospect for a future Oldest Ice drill site.

Future sea-level predictions require that the history of the Antarctic Ice Sheet is well understood and constrained by observations. Much of the ice sheets’ ice-dynamic properties are governed by processes at the ice-bed interface which can be imaged with radar sounding surveys. Here we use a combination of ultra-wideband radio-echo sounding data, satellite radar and laser altimetry data, as well as electromagnetic waveform modeling to characterize the properties of the ice base and the evolution of the subglacial morphology of the Jutulstraumen drainage basin (western Dronning Maud Land, Antarctica). Based on the classification of the bed topography, we reconstruct the step-by-step modifications the subglacial landscape has experienced since the beginning of the glaciation of Antarctica, 34 million years ago. Between 2017 and 2020, we find evidence of active episodic cascade-like subglacial water transport along the subglacial valley network. In addition, our high-resolution radio-echo sounding data reveal a cluster of anomalous basal ice units whose material properties we constrain by electromagnetic waveform modeling. Through this, we aim to derive the physical conditions at the ice base, and establish a link to the subglacial hydrology system. The combination of these observations will represent an important step towards a better understanding of large-scale ice-sheet dynamics in western Dronning Maud Land.

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 ∼24 000 km2 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).

Future sea-level predictions require that the history and physical state of the Antarctic ice sheet is well understood and constrained by observations. Much of the ice sheets’ ice-dynamic properties are governed by processes at the ice-bed interface which can be imaged with radar sounding surveys. Moreover, certain processes at the ice-sheet base can have an effect all the way to the ice surface, which in turn can be observed with satellites. Here we use a combination of ultra-wideband radio-echo sounding data, satellite radar and laser altimetry data to characterize the evolution of the subglacial morphology of the Jutulstraumen drainage basin (western Dronning Maud Land, Antarctica). Based on the classification of the bed topography, we reconstruct the step-by-step modifications the subglacial landscape has experienced since the beginning of the glaciation of Antarctica, 34 million years ago. In addition, between 2017 and 2020, we find evidence of active episodic cascade-like subglacial water transport along the subglacial valley network. The combination of these observations will represent an important step towards a better understanding of large-scale ice-sheet dynamics in western Dronning Maud Land.

Anisotropic crystal fabrics in ice sheets develop as a consequence of deformation and hence record information of past ice flow. Simultaneously, the fabric affects the present-day bulk mechanical properties of glacier ice because the susceptibility of ice crystals to deformation is highly anisotropic. This is particularly relevant in dynamic areas such as fast-flowing glaciers and ice streams, where the formation of strong fabrics might play a critical role in facilitating ice flow. This fact is ignored in most state-of-the-art ice sheet models, and while their importance has been recognized years ago, accounting for fabrics evolution and their impact on the ice viscosity has only recently become feasible. Both, the application of such models in ice streams as well as their verification through in-situ observations are, however, still rare. We present an extensive dataset of fabric anisotropy derived from radar data recorded in the onset region of the Northeast Greenland Ice Stream by air-borne and ground-based systems. Our methods yield the horizontal anisotropy and are based on travel time anisotropy and splitting as well as birefringence-induced power modulation of radar signals. They complement each other and show good agreement. We compare these in-situ observations with the results obtained from a fabric-evolution model employed along flow tubes in the ice stream onset to discuss the fabric in light of past flow history and its significance for the current flow mechanics of the ice stream.

The landscape of Antarctica, hidden beneath kilometre-thick ice in most places, has been shaped by the interactions between tectonic and erosional processes. The flow dynamics of the thick ice cover deepened pre-formed topographic depressions by glacial erosion, but also preserved the subglacial landscapes in regions with moderate to slow ice flow. Mapping the spatial variability of these structures provides the basis for reconstruction of the evolution of subglacial morphology. This study focuses on the Jutulstraumen Glacier drainage system in Dronning Maud Land, East Antarctica. The Jutulstraumen Glacier reaches the ocean via the Jutulstraumen Graben, which is the only significant passage for draining the East Antarctic Ice Sheet through the western part of the Dronning Maud Land mountain chain. We acquired new bed topography data during an airborne radar campaign in the region upstream of the Jutulstraumen Graben to characterise the source area of the glacier. The new data show a deep relief to be generally under-represented in available bed topography compilations. Our analysis of the bed topography, valley characteristics and bed roughness leads to the conclusion that much more of the alpine landscape that would have formed prior to the Antarctic Ice Sheet is preserved than previously anticipated. We identify an active and deeply eroded U-shaped valley network next to largely preserved passive fluvial and glacial modified landscapes. Based on the landscape classification, we reconstruct the temporal sequence by which ice flow modified the topography since the beginning of the glaciation of Antarctica.

The North East Greenland Ice Stream clearly stands out in the surface velocity field of the ice flow of Greenland, with its sharp and narrow shear margins visible in the flow field almost up to the central divide. While the current extent and strength of the streaming can be determined from remotely sensed velocities of the ice surface, it is less known how the ice stream is affecting the deeper layers of ice in its catchment area, and how it may have evolved over time. The deformation of the ice due to streaming can be made visible by mapping the distortion of the isochronous stratigraphy of the ice. This has been done by an airborne radar survey centering on the location of the EGRIP drilling camp, carried out with the ultra wide band radar system (AWI UWB). The dense grid of profiles arranged mainly perpendicular to the ice flow reveals the imprint that the strong shearing leaves within the layering of the ice. Although the layers are tightly folded and distorted within the shear zones, it is possible to continuously trace reflections within the upper half of the ice column throughout the entire survey area. It can be shown that the intensity of the folding is linked to the strain rate field derived from the surface velocities, and that the deformation history of the ice is preserved in the folded layers, even after it is no longer affected by high strain rates. The advection patterns of the mapped stratigraphic features reveal how the streaming of the ice and the resulting local changes of surface topography may have affected the inflow into the stream and the position of the shear margins over time.

Crystal anisotropy of ice causes slight birefringence for electromagnetic waves. At the same time, the mechanical anisotropy amounts to several orders of magnitude, thus making fabric properties highly-relevant for internal deformation. To date, bulk anisotropy of glaciers and ice sheets can be determined by geophysical methods, such as polarimetric radar, or direct sampling from ice cores. A shortcoming has been so far that changes of bulk anisotropy could mainly be inferred at single point observations, but less so as continuous profiles. Here, we exploit the effect of birefringence caused by bulk anisotropy in co-polarized airborne radar data to determine the horizontal anisotropy across the North-East Greenland Ice Stream. We base our analysis on the fact that birefringence causes a second-order effect on radar amplitudes, which leads to a beat frequency in the low and medium frequency range (O(100 kHz)), which is proportional to the horizontal anisotropy. Complementing our radar analysis with direct fabric and dielectric property observations we can constrain the range of all three fabric eigenvalues as a function of space across and along the ice stream. Finally, we assess the effect of the inferred fabric distribution on the overall ice rheology in the context of ice stream dynamics. Our overall approach has the advantage that it can be applied to co-polarized radar systems, as commonly used in profiling surveys, and does not require dedicated cross-polarized radar set-up. This provides the opportunity to revisit older data, especially from Greenland and Antarctica, to map fabric anisotropy in ice-dynamically interesting regions.

The area near Dome C, East Antarctica, is thought to be one of the most promising targets for recovering a continuous ice-core record spanning more than a million years. The European Beyond EPICA consortium has selected Little Dome C, an area ~35 km south-east of Concordia Station, to attempt to recover such a record. Here, we present the results of the final ice-penetrating radar survey used to refine the exact drill site. These data were acquired during the 2019–2020 Austral summer using a new, multi-channel high-resolution VHF radar operating in the frequency range of 170–230 MHz. This new instrument is able to detect reflections in the near-basal region, where previous surveys were unable to trace continuous horizons. The radar stratigraphy is used to transfer the timescale of the EPICA Dome C ice core (EDC) to the area of Little Dome C, using radar isochrones dating back past 600 ka. We use these data to derive the expected depth–age relationship through the ice column at the now-chosen drill site, termed BELDC. These new data indicate that the ice at BELDC is considerably older than that at EDC at the same depth, and that there is about 375 m of ice older than 600 ka at BELDC. Stratigraphy is well preserved to 2565 m, below which there is a basal unit with unknown properties. A simple ice flow model tuned to the isochrones suggests ages likely reach 1.5 Ma near 2525 m, ~40 m above the basal unit and ~240 m above the bed, with sufficient resolution (14±1 ka m−1) to resolve 41 ka glacial cycles.

The area near Dome C, East Antarctica, is thought to be one of the most promising targets for recovering a continuous ice-core record spanning more than a million years. The European Beyond EPICA consortium has selected Little Dome C (LDC), an area similar to 35 km southeast of Concordia Station, to attempt to recover such a record. Here, we present the results of the final ice-penetrating radar survey used to refine the exact drill site. These data were acquired during the 2019-2020 austral summer using a new, multi-channel high-resolution very high frequency (VHF) radar operating in the frequency range of 170-230 MHz. This new instrument is able to detect reflectors in the near-basal region, where previous surveys were largely unable to detect horizons. The radar stratigraphy is used to transfer the timescale of the EPICA Dome C ice core (EDC) to the area of Little Dome C, using radar isochrones dating back past 600 ka. We use these data to derive the expected depth-age relationship through the ice column at the now-chosen drill site, termed BELDC (Beyond EPICA LDC). These new data indicate that the ice at BELDC is considerably older than that at EDC at the same depth and that there is about 375m of ice older than 600 kyr at BELDC. Stratigraphy is well preserved to 2565 m, similar to 93% of the ice thickness, below which there is a basal unit with unknown properties. An ice-flow model tuned to the isochrones suggests ages likely reach 1.5 Myr near 2500 m, similar to 65m above the basal unit and similar to 265m above the bed, with sufficient resolution (19 +/- 2 kyrm(-1)) to resolve 41 kyr glacial cycles.

The North East Greenland Ice Stream (NEGIS) is delineated by well-defined shear margins, which are evident in the gradient of surface velocity field as well as in the surface topography, where they form troughs up to ten meters deep. In the upper part of the ice stream the margins appear not to be linked to bedrock topography. To understand this efficient system of mass transport towards the ocean it is essential to investigate the nature of the shear margins, as here very localized deformation decouples the inner ice stream from the slower flowing surrounding ice sheet. This process is influenced by several factors and feedback mechanisms, including the crystal fabric orientation, strain heating and localization of meltwater. In summary, the shear margins are area-wise a small part of the ice stream itself, but the processes leading to the localization of deformation are of similar importance for ice discharge as the processes enabling fast flow of the main trunk over the bed. We present results from an airborne radar survey with the AWI Ultra-Wide Band Radar system, covering an area 150 km upstream and 100 km downstream of the deep drilling site on the ice stream (EGRIP). Over the survey area the ice stream accelerates from 12 m/a to 75 m/a. We focus on the signatures of the shear margins in the radar data. In the regions of localized shear, the internal reflections in the radargrams show disturbances in the form of steep undulations, or chevron folds, which are intensified with ongoing shear. As the ice stream has been covered with 36 flow-perpendicular radar sections we are able to show the evolution of these characteristic signatures over the survey area, and thus, as an analog, over time. 3D-representations of the folded stratigraphic layers reveal how new folds are formed when the ice stream widens and how older structures are preserved in the outer part of the main trunk, where they are no longer subject to shear. Furthermore, we link the change of the shape

Radio echo sounding of polar ice sheets provides important information on the ice bed topography and internal layers. These data have been used by scientists to create 3D maps of polar ice sheets for climate modeling as well as to reconstruct the climate history that dates back to hundreds of thousands of years. In this paper, we present the design, and development of three surface-based multi-channel radars in the VHF and UHF bands. We provide results from radar data multi-frequency and polarization radar data collected over the Greenland ice sheet. All the three radars shared the same digital waveform generator and digitizer, and were installed in and operated on a tracked vehicle. The radars are operated with 3 different antenna arrays designed for operation over 170-230 MHz, 180-340 MHz and 600-900 MHz. The results we sounded more than 2.7 km thick ice with radars operating at frequencies as high as 850 MHz with more than 40 dB signal-to-noise ratio.

In this letter, we report on the design, development, and field operation of a surface-based multi-channel ultraw- ideband (UWB) ultrahigh frequency (UHF) radar to measure ice thickness, basal conditions, and ice-shelf bottom melt rates. The radar concept is based on the recent success in sounding shallow low-loss ice (∼1 km) and measuring the ice-shelf melt rates with a 600–900-MHz low-power radar, referred to as the accumulation radar. Our proposed radar system operates over the same frequency band, from 600 to 900 MHz, with a peak transmit power of 800 W. We used a large and lightweight 16 m × 17 m antenna array arranged in a Mills cross-configuration to obtain the required radar sensitivity to sound more than 3-km- thick ice and image the internal layers at a fine vertical resolution of about 60 cm. We used the system at the East Greenland Ice-coring Project (EGRIP) site in Summer 2018 to collect data over ∼100 km of lines. We sounded about 2.8-km-thick ice with more than 40-dB signal-to-noise ratio and mapped the internal layers to a depth of almost 2.5 km. Our results show that an airborne or spaceborne radar operating at frequencies as high as 900 MHz with a large antenna array can be used to map large ice sheets in Greenland and Antarctica.

Radio echo sounding of polar ice sheets provides important information on the ice bed topography and internal layers. These data have been used by scientists to create 3-D maps of polar ice sheets for climate modeling as well as to reconstruct the climate history that dates back to hundreds of thousands of years. In this article, we present the design and development of three surface-based multichannel radars in the VHF and UHF bands. We provide results from radar data multifrequency and polarization radar data collected over the Greenland ice sheet. All the three radars shared the same digital waveform generator and digitizer, and were installed in and operated on a tracked vehicle. The radars are operated with three different antenna arrays designed for operation over 170 & x2013;230, 180 & x2013;340, and 600 & x2013;900 & x00A0;MHz. The results we obtained sounded more than 2.7 km thick ice with radars operating at frequencies as high as 850 MHz with more than 40 dB signal-to-noise ratio.

The North East Greenland Ice Stream (NEGIS) is an essential part of the Greenland mass balance, reaching up to the central divide of the ice sheet. Its shape is dominated by very distinct shear margins that appear not to be linked to bedrock topography in the upper part of the ice stream. To understand the effective mass transport within an ice stream it is necessary to investigate the nature of the shear margins, in which very localized deformation decouples the inner ice stream from the slower-flowing surrounding ice sheet. We present results from an airborne radar survey with the Alfred Wegener Institute Ultra Wide Band Radar system in the vicinity of the EGRIP ice core location, focusing on the signatures of the shear margins in the radar data. In regions of localized shear, the internal reflections show disturbances in the form of steep undulations, which are intensified with ongoing shear. We show the evolution of these characteristic signatures over the survey area. Our data covers the main trunk and the shear margins over 200 km along the ice stream, in which the flow velocity in the center increases from 12 m a–1 to 75 m a–1. We analyze the change in the shape of the internal reflections in the shear zones in combination with a strain rate field calculated from high-resolution flow velocities derived from TerraSAR-X SAR interferometry. In the vicinity of the EGRIP drilling site, 240 km downstream of the divide, the NEGIS widens and the shear margin in the northwest steps 6 km outward, which is visible in optical satellite imagery as well as in the satellite-derived velocity field. A special focus of our analysis lies on the step in the shear margin whose clearly visible radar-backscatter signatures indicate a newly forming vertical disturbance caused by a draw-down of the layers.

We compared elastic moduli in polar firn derived from diving wave refraction seismic velocity analysis, firn-core density measurements and microstructure modelling based on firn-core data. The seismic data were obtained with a small electrodynamic vibrator source near Kohnen Station, East Antarctica. The analysis of diving waves resulted in velocity–depth profiles for different wave types (P-, SH- and SV-waves). Dynamic elastic moduli of firn were derived by combining P- and S-wave velocities and densities obtained from firn-core measurements. The structural finite-element method (FEM) was used to calculate the components of the elastic tensor from firn microstructure derived from X-ray tomography of firn-core samples at depths of 10, 42, 71 and 99 m, providing static elastic moduli. Shear and bulk moduli range from 0.39 to 2.42 GPa and 0.68 to 2.42 GPa, respectively. The elastic moduli from seismic observations and the structural FEM agree within 8.5% for the deepest achieved values at a depth of 71 m, and are within the uncertainty range. Our observations demonstrate that the elastic moduli of the firn can be consistently obtained from two independent methods which are based on dynamic (seismic) and static (tomography and FEM) observations, respectively, for deeper layers in the firn below ∼10 m depth.

The large Ice sheets, Greenland and Antartica, are two key players for understanding the future effects of climate change when it comes to sea level rise. A considerable uncertainty is their dynamic response to changing boundary conditions, e.g. increased melting at the surface percolating to their base or warmer ocean water temperatures underneath ice shelves. To decipher ongoing processes and reveal their development, glaciologists often turn to the observation of internal structures in and basal conditions of ice sheets using radio-echo sounding (RES) techniques. Over the last twenty years the radar techniques employed for these observation have seen a considerable improvement, such that today’s ice-penetrating radar has little in common with the former RES techniques, where data were still stored on photographic films. The newest of such systems, in operation since 2016, is AWI’s ultrawideband radar , operating in the radio (150-600 MHz) and microwave frequency range, an advancement from the original MCoRDS/I system developed by the Center for Remote Sensing of Ice Sheets (CReSIS). With 24 elements in use for the AWI UWB, the lateral and vertical resolution for imaging the interior of ice sheets at kilometers depth has been brought to the range of the sub-meter scale. This talk introduces the scientific objectives motivating the development of this system, presents its technical aspects, including data recording and processing, and finally shows first results from the last Greenland and Antarctic field seasons.

The European Project for Ice Coring in Antarctica (EPICA) ice core was drilled between 2001 and 2006 at the Kohnen Station, Antarctica. During the drilling process the borehole was logged repeatedly. Repeated logging of the borehole shape is a means of directly measuring the deformation of the ice sheet not only on the surface but also with depth, and to derive shear strain rates for the lower part, which control the volume of ice transported from the inner continent towards the ocean. The logging system continuously recorded the tilt of the borehole with respect to the vertical (inclination) as well as the heading of the borehole with respect to magnetic north (azimuth) by means of a compass. This dataset provides the basis for a 3-D reconstruction of the borehole shape, which is changing over time according to the predominant deformation modes with depth. The information gained from this analysis can then be evaluated in combination with lattice preferred orientation, grain size and grain shape derived by microstructural analysis of samples from the deep ice core. Additionally, the diameter of the borehole, which was originally circular with a diameter of 10 cm, was measured. As the ice flow velocity at the position of the EDML core is relatively slow (about 0.75 m/a), the changes of borehole shape between the logs during the drilling period were very small and thus difficult to interpret. Thus, the site has been revisited in the Antarctic summer season 2016 and logged again using the same measurement system. The change of the borehole inclination during the time period of 10 years clearly reveals the transition from a pure shear dominated deformation in the upper part of the ice sheet to shear deformation at the base. We will present a detailed analysis of the borehole parameters and the deduced shear strain rates in the lower part of the ice sheet. The results are discussed with respect to ice microstructural data derived from the EDML ice core. Microstructural

Ice cores at a fist glance are almost ideal paleoclimate archives. They not only record temperatures of the past in a relatively direct way, they also keep the atmospheric composition intact. In the end it is only a question how smart we are in using appropriate analytic techniques to unravel the complete atmospheric record – greenhouse gases, aerosols and dust. Another advantage of ice cores is that we can count on very high resolution in time and although we have no direct dating method we can fairly confidently ascribe an absolute time to samples retrieved from a certain depth. From the first ice cores drilled in the sixties until today results from ice cores have helped us to better understand the climate system in particular its inherent global teleconnections. Ice cores are also our climate conscience – they show us without doubt how humankind is altering climate boundary conditions i.e. the greenhouse gas concentrations at an unprecedented rate. In my presentation I will briefly touch on techniques from drilling to analytics and show some of the most important results.