Search

Over the past 60 years, scientists have strived to understand the past, present and future of the Antarctic Ice Sheet. One of the key components of this research has been the mapping of Antarctic bed topography and ice thickness parameters that are crucial for modelling ice flow and hence for predicting future ice loss and ensuing sea level rise. Supported by the Scientific Committee on Antarctic Research (SCAR), the Bedmap3 Action Group aims not only to produce new gridded maps of ice thickness and bed topography for the international scientific community, but also to standardize and make available all the geophysical survey data points used in producing the Bedmap gridded products. Here, we document the survey data used in the latest iteration, Bedmap3, incorporating and adding to all of the datasets previously used for Bedmap1 and Bedmap2, including ice-bed, surface and thickness point data from all Antarctic geophysical campaigns since the 1950s. More specifically, we describe the processes used to standardize and make these and future survey and gridded datasets accessible under the ‘Findable, Accessible, Interoperable and Reusable’ (FAIR) data principles. With the goals to make the gridding process reproducible and to allow scientists to re-use the data freely for their own analysis, we introduce the new SCAR Bedmap Data Portal (bedmap.scar.org, last access: 18 October 2022) created to provide unprecedented open access to these important datasets, through a user-friendly webmap interface. We believe that this data release will be a valuable asset to Antarctic research and will greatly extend the life cycle of the data held within it. Data are available from the UK Polar Data Centre: https://data.bas.ac.uk.

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. 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-flow 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 future approaches to understanding and modelling the response of Earth’s ice sheets to global warming.

Ice shelves surrounding the Antarctic perimeter moderate ice discharge towards the ocean through buttressing. Ice-shelf evolution and integrity depend on the local surface accumulation, basal melting and on the spatially variable ice-shelf viscosity. These components of ice-shelf mass balance are often poorly constrained by observations and introduce uncertainties in ice-sheet projections. Isochronal radar stratigraphy is an observational archive for the atmospheric, oceanographic and ice-flow history of ice shelves. Here, we predict the stratigraphy of locally accumulated ice on ice shelves with a kinematic forward model for a given atmospheric and oceanographic scenario. This delineates the boundary between local meteoric ice (LMI) and continental meteoric ice (CMI). A large LMI to CMI ratio hereby marks ice shelves whose buttressing strength is more sensitive to changes in atmospheric precipitation patterns. A mismatch between the steady-state predictions of the kinematic forward model and observations from radar can highlight inconsistencies in the atmospheric and oceanographic input data or be an indicator for a transient ice-shelf history not accounted for in the model. We discuss pitfalls in numerical diffusion when calculating the age field and validate the kinematic model with the full Stokes ice-flow model Elmer/Ice. The Roi Baudouin Ice Shelf (East Antarctica) serves as a test case for this approach. There, we find a significant east–west gradient in the LMI / CMI ratio. The steady-state predictions concur with observations on larger spatial scales (>10 km), but deviations on smaller scales are significant, e.g., because local surface accumulation patterns near the grounding zone are underestimated in Antarctic-wide estimates. Future studies can use these mismatches to optimize the input data or to pinpoint transient signatures in the ice-shelf history using the ever growing archive of radar observations of internal ice stratigraphy.

The dynamic mass loss of ice sheets due to ice flow constitutes one of the biggest uncertainties in projections of ice-sheet evolution and sea-level rise. One central, understudied aspect of ice flow is how the bulk orientation of the ice-crystal lattice (fabric) translates to the mechanical anisotropy of ice. Here, we present a comprehensive analysis of the depth-averaged fabric distribution and corresponding directional flow enhancement factors covering a large area of the Northeast Greenland Ice Stream onset region. Our results are based on a combination of methods applied to extensive airborne and ground-based radar surveys, ice-core observations, and numerical ice-flow modelling. They show a strong spatial variability of the horizontal fabric anisotropy and a rapid crystal reorganization on the order of hundreds of years coinciding with the ice-stream geometry. The ice in parts of the ice stream is found to be up to one order of magnitude harder for pure-shear deformation in flow direction compared to isotropic ice, while the shear margins are potentially softened for horizontal simple shear by a factor of two. The high spatial resolution and large-scale coverage of our results is unprecedented and significantly contributes towards more accurate ice-flow models and a better understanding of ice-stream dynamics. Moreover, our estimated time scale of crystal reorganization and methodological approaches can potentially also be transferred to the investigation of ice on extraterrestrial bodies and other geologic materials.

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.

The internal ice stratigraphy as imaged by radar is an integrated archive of the atmospheric- oceanographic, and ice-dynamic history that the ice sheet has experienced. It provides an observational constraint for ice flow modeling that has been used for instance to predict age-depth relationships at prospective ice-coring sites in Antarctica’s interior. The stratigraphy is typically more disturbed and more difficult to image in coastal regions due to faster ice flow. Yet, knowledge of ice stratigraphy across ice shelf grounding lines and further seawards is important to help constrain ocean-induced melting and associated stability.Here, we present preliminary results of synthesizing information from radar stratigraphic characteristics from airborne and ground-based radar surveys that have been collected for specific projects starting from the 1990s onwards focusing on ice marginal zones of Antarctica. The key data is based on airborne surveys from the German Alfred Wegener Institute’s polar aircrafts equipped with a 150 MHz radar. In the meantime this system has been replaced by an ultra-wide band 150-520 MHz radar. The older data will provide a baseline with extensive coverage that can be used for model calibration and change detection over time. We aim to provide metrics of the radio stratigraphy (e.g. shape and slope of internal reflection horizons) as well as classified prevalent stratigraphy types that can be used to calibrate machine learning approaches such as simulation based inference. The data obtained will be integrated in coordination efforts within the SCAR AntArchitecture Action Group.

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 Alfred Wegener Institute’s multichannel ultra-wideband (UWB) radar mounted on the Polar6 aircraft. Radar profiles cover an area of ~24000 km2 and extend over the well-defined shear margins of the NEGIS. The survey area is centred 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 can not 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 modelled 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 at the PANGAEA Data Publisher (https://doi.pangaea.de/10.1594/PANGAEA.928569; Franke et al. 2021a).

Migration of subglacial water underneath thick Antarctic ice is difficult to observe directly but is known to influence ice flow dynamics. Here, we analyze a 6-year time series of displacement maps from differential Sentinel-1 SAR interferometry (DInSAR) in the upstream region of Jutulstraumen Glacier. Our results reveal short-term (between 12 days and 1 year) interconnected subsidence- and uplift events of the ice surface, which we interpret as a pressure response to the drainage and filling of subglacial lakes. This indicates an episodic cascade-like water transport with longer quiescent phases in a dynamically stable glacial setting. Abrupt events appear in the DInSAR time series and are confirmed by ICESat-2 altimetry. The events can be traced for a 1-year period along a urn:x-wiley:00948276:media:grl63164:grl63164-math-0001175 km flow path. We are able to observe the migration of subglacial water with unprecedented spatial and temporal resolution, providing a new observational baseline to further develop subglacial hydrological models.

Past temperature reconstructions from polar ice sheets are commonly based on stable water isotope records in ice-cores. However, despite major efforts in the understanding of the ice-core signal fo...

Ever since the first deep ice cores were drilled, it has been a challenge to determine their original, in-situ orientation. In general, the orientation of an ice core is lost as the drill is free to rotate during transport to the surface. For shallow ice cores, it is usually possible to match the adjacent core breaks, which preserves the orientation of the ice column. However, this method fails for deep ice cores, such as the EastGRIP ice core in Northeast Greenland. We provide a method to reconstruct ice core orientation using visual stratigraphy and borehole geometry. As the EastGRIP ice core is drilled through the Northeast Greenland Ice Stream, we use information about the directional structures to perform a full geographical re-orientation. We compared the core orientation with logging data from core break matching and the pattern of the stereographic projections of the crystals’c-axis orientations. Both comparisons agree very well with the proposed orientation method. The method works well for 441 out of 451 samples from a depth of 1375–2120 m in the EastGRIP ice core. It can also be applied to other ice cores, providing a better foundation for interpreting physical properties and understanding the flow of ice.

The Northeast Greenland Ice Stream (NEGIS) is the largest active ice stream on the Greenland Ice Sheet (GrIS) and a crucial contributor to the ice-sheet mass balance. To investigate the ice-stream dynamics and to gain information about the past climate, a deep ice core is drilled in the upstream part of the NEGIS, termed the East Greenland Ice-core Project (EastGRIP). Upstream flow can introduce climatic bias into ice cores through the advection of ice deposited under different conditions further upstream. This is particularly true for EastGRIP due to its location inside an ice stream on the eastern flank of the GrIS. Understanding and ultimately correcting for such effects requires information on the atmospheric conditions at the time and location of snow deposition. We use a two-dimensional Dansgaard–Johnsen model to simulate ice flow along three approximated flow lines between the summit of the ice sheet (GRIP) and EastGRIP. Isochrones are traced in radio-echo-sounding images along these flow lines and dated with the GRIP and EastGRIP ice-core chronologies. The observed depth–age relationship constrains the Monte Carlo method which is used to determine unknown model parameters. We calculate backward-in-time particle trajectories to determine the source location of ice found in the EastGRIP ice core and present estimates of surface elevation and past accumulation rates at the deposition site. Our results indicate that increased snow accumulation with increasing upstream distance is predominantly responsible for the constant annual layer thicknesses observed in the upper part of the ice column at EastGRIP, and the inverted model parameters suggest that basal melting and sliding are important factors determining ice flow in the NEGIS. The results of this study form a basis for applying upstream corrections to a variety of ice-core measurements, and the inverted model parameters are useful constraints for more sophisticated modelling approaches in the future.

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., Airborne ice‐penetrating radar data reveal the evolution of the subglacial morphology of the Jutulstraumen Glacier drainage system in western Dronning Maud Land. We identify various geomorphological patterns that are related to different stages of subglacial erosion and allow us to reconstruct the temporal sequence by which ice flow modified the topography since the beginning of glaciation of Antarctica.

Ice shelves are widely known to slow the transfer of Antarctic grounded ice to the ocean, especially if their flow is decelerated by local pinning points. Their longevity is influenced by variations in ice dynamics, surface accumulation and oceanic conditions in the ice shelf cavity. This is reflected in the ice shelf structure, which can be characterized by the shape of internal radar reflection horizons.We aim to map the internal ice shelf stratigraphy of ice shelves, starting with the narrow belt of ice-shelves in the Dronning Maud Land area. The final goal will be to evaluate these as a spatiotemporal archive of ice provenance and ice dynamics. The bulk of the data presented here were collected with AWI’s airborne multi frequency ultra-wideband radar and we combine these new observations with airborne and ground-based radar surveys from previous years. We present a consistent set of internal radar isochrones on a catchment scale for the Roi Baudoin area including the Ragnhild ice streams, the grounding-zone, the iceshelf and multiple ice rises. Using pattern matching technique we can link isochrones across different ice rises in the area, and hence provide first observational constraints on how ice rises jointly react to changes in atmospheric and oceanographic forcings. We also find a number of interesting features including dynamically induced dips in shear zones, truncating layers at the ice-shelf base, and the development of a meteoric ice layer distinguishing advected from newly accumulated ice in the iceshelf. The time series provided by radar observations over the last 10 years also offers the potential to map temporal changes. We use ice-flow modelling to provide age constraints for some internal layers and delineate portions within the shelf as a function of their advection history, hence marking areas of differing rheologies within the shelf. Taken together, this case study on a catchment scale is a primer to unravel the information stored in the isochronal stratigraphy of coastal Antarctica and contributes to international efforts (e.g., SCAR AntArchitecture) of mapping stratigraphy on ice sheet scales.

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.

The North East Greenland ice-stream (NEGIS) is the largest active ice-stream on the Greenland ice-sheet and is a crucial contributor to the ice-sheet mass balance. To investigate the ice-stream dynamics and to gain information about the past climate, a deep ice-core is drilled in the upstream part of the NEGIS, termed the East Greenland ice-core project (EastGRIP). Upstream flow effects introduce non-climatic bias in ice-cores and are particularly strong at EastGRIP due to high ice-flow velocities and the location in an ice-stream on the eastern flank of the Greenland ice-sheet. Understanding and ultimately correcting for such effects requires information on the source area and the local atmospheric conditions at the time of ice deposition. We use a two-dimensional Dansgaard-Johnsen model to simulate ice-flow along three approximated flow-lines between the summit of the ice-sheet and EastGRIP. Model parameters are determined using a Monte Carlo inversion by minimizing the misfit between modeled isochrones and isochrones observed in radio-echo-sounding images. We calculate backward-in-time particle trajectories to determine the source area of ice found in the EastGRIP core today and present estimates of surface elevation and past accumulation-rates at the deposition site. The thinning function and accumulated strain obtained from the modeled velocity field provide useful information on the deformation history in the EastGRIP ice. Our results indicate that increased accumulation in the upstream area is predominantly responsible for the constant annual layer thickness observed in the upper part of the ice column at EastGRIP. Inverted model parameters suggest that the imprint of basal melting and sliding is present in large parts along the flow profiles and that most internal ice deformation happens close to the bedrock. The results of this study can act as a basis for applying upstream corrections to a variety of ice-core measurements, and the model parameters can be useful constraints for more sophisticated modeling approaches in the future.

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 and compare it with numerical model results. 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 ice stream geometry and large ice surface velocities at the onset region of the Northeast Greenland Ice Stream (NEGIS) are not yet well reproduced by ice sheet models. The quantification of basal sliding and a parametrization of basal conditions remains a major gap. In this study, we assess the basal conditions of the onset region of the NEGIS in a systematic analysis of airborne ultra‐wideband radar data. We evaluate basal roughness and basal return echoes in the context of the current ice stream geometry and ice surface velocity. We observe a change from a smooth to a rougher bed where the ice stream widens, and a distinct roughness anisotropy, indicating a preferred orientation of subglacial structures. In the upstream region, the excess ice mass flux through the shear margins is evacuated by ice flow acceleration and along‐flow stretching of the ice. At the downstream part, the generally rougher bed topography correlates with a decrease in flow acceleration and lateral variations in ice surface velocity. Together with basal water routing pathways, this hints to two different zones in this part of the NEGIS: the upstream region collecting water, with a reduced basal traction, and downstream, where the ice stream is slowing down and is widening on a rougher bed, with a distribution of basal water toward the shear margins. Our findings support the hypothesis that the NEGIS is strongly interconnected to the subglacial water system in its onset region, but also to the subglacial substrate and morphology., Plain Language Summary: The Northeast Greenland Ice Stream (NEGIS) transports a large amount of ice mass from the interior of the Greenland Ice Sheet (GrIS) toward the ocean. The extent and geometry of the NEGIS are difficult to reproduce in current ice sheet models because many boundary conditions, such as the properties of the ice base, are not well known. In this study, we present new characteristics of the ice base from the onset region of the NEGIS derived by airborne radio‐echo sounding data. Our data yield a smooth and increasingly lubricated bed in the upstream part of our survey area, which enables the ice to accelerate. Our results confirm the hypothesis that the position of the ice stream boundaries are coupled to the subglacial hydrology system., Key Points: Basal roughness at the onset of the NEGIS hints to a geomorphic anisotropy and a change in the geomorphological regime. Basal water is funneled into the ice stream upstream and redistributed toward the shear margins further downstream. A smooth and progressively lubricated bed reduces basal traction and favors the acceleration of the NEGIS at its onset., A. P. Møller Foundation, US National Science Foundation, Alfred Wegener Institute, National Institute of Polar Research and Arctic Challenge for Sustainability, University of Bergen and Bergen Research Foundation, Swiss National Science Foundation, French Polar Institute Paul‐Emile Victor, Chinese Academy of Sciences and Beijing Normal University, NASA Operation IceBridge, NSF

The Northeast Greenland Ice Stream (NEGIS) is the largest active ice stream on the Greenland Ice Sheet (GrIS) and a crucial contributor to the ice-sheet mass balance. To investigate the ice-stream dynamics and to gain information about the past climate, a deep ice core is drilled in the upstream part of the NEGIS, termed the East Greenland Ice-Core Project (EastGRIP). Upstream flow effects introduce non-climatic bias in ice cores and are particularly strong at EastGRIP due to high ice-flow velocities and the location inside an ice stream on the eastern flank of the GrIS. Understanding and ultimately correcting for such effects requires information on the source area and the local atmospheric conditions at the time of ice deposition. We use a two-dimensional Dansgaard–Johnsen model to simulate ice flow along three approximated flow lines between the summit of the ice sheet (GRIP) and EastGRIP. Model parameters are determined using a Monte Carlo inversion by minimizing the misfit between modelled isochrones and isochrones observed in radio-echo-sounding (RES) images. We calculate backward-in-time particle trajectories to determine the source area of ice found in the EastGRIP ice core and present estimates of surface elevation and past accumulation rates at the deposition site. Our results indicate that increased accumulation in the upstream area is predominantly responsible for the constant annual layer thicknesses observed in the upper part of the ice column at EastGRIP. Inverted model parameters suggest that the imprint of basal melting and sliding is present in large segments along the flow profiles and that most internal ice deformation happens in the lower half of the ice column. The results of this study act as a basis for applying upstream corrections to a variety of ice-core measurements, and the model parameters are useful constraints for more sophisticated modelling approaches in the future.

The Northeast Greenland Ice Stream (NEGIS) is an important dynamic component for the total mass balance of the Greenland ice sheet, as it reaches up to the central divide and drains 12% of the ice sheet. The geometric boundary conditions and in particular the nature of the subglacial bed of the NEGIS are essential to understand its ice flow dynamics. We present a record of more than 8000 km of radar survey lines of multi-channel, ultra-wideband radio echo sounding data covering an area of 24 000 km2, centered on the drill site for the East Greenland Ice-core Project (EGRIP), in the upper part of the NEGIS catchment. Our data yield a new detailed model of ice-thickness distribution and basal topography in the region. The enhanced resolution of our bed topography model shows features which we interpret to be caused by erosional activity, potentially over several glacial–interglacial cycles. Off-nadir reflections from the ice–bed interface in the center of the ice stream indicate a streamlined bed with elongated subglacial landforms. Our new bed topography model will help to improve the basal boundary conditions of NEGIS prescribed for ice flow models and thus foster an improved understanding of the ice-dynamic setting.

Suture zones are abundant on Antarctic ice shelves and widely observed to impede fracture propagation, greatly enhancing ice-shelf stability. Using seismic and radar observations on the Larsen C Ice Shelf of the Antarctic Peninsula, we confirm that such zones are highly heterogeneous, consisting of multiple meteoric and marine ice bodies of diverse provenance fused together. Here we demonstrate that fracture detainment is predominantly controlled by enhanced seawater content in suture zones, rather than by enhanced temperature as previously thought. We show that interstitial seawater can reduce fracture-driving stress by orders of magnitude, promoting both viscous relaxation and the development of micro cracks, the incidence of which scales inversely with stress intensity. We show how simple analysis of viscous buckles in ice-penetrating radar data can quantify the seawater content of suture zones and their modification of the ice-shelf’s stress regime. By limiting fracture, enhancing stability and restraining continental ice discharge into the ocean, suture zones act as vital regulators of Antarctic mass balance., Suture zones are abundant on Antarctic ice shelves and widely observed to impede fracture propagation. Here we show that fracture detainment is principally controlled by the zones’ enhanced seawater contents, reducing fracture-driving stresses by orders of magnitude and therefore greatly enhancing stability.

In times of warming in polar regions, the prediction of ice sheet discharge is of utmost importance to society, because of its impact on sea level rise. In simulations the flow rate of ice is usually implemented as proportional to the differential stress to the power of the exponent n=3. This exponent influences the softness of the modeled ice, as higher values would produce faster flow under equal stress. We show that the stress exponent, which best fits the observed state of the Greenland Ice Sheet, equals n=4. Our results, which are not dependent on a possible basal sliding component of flow, indicate that most of the interior northern ice sheet is currently frozen to bedrock, except for the large ice streams and marginal ice. Ice in the polar ice sheets flows towards the oceans under its own weight. Knowing how fast the ice flows is of crucial importance to predict future sea level rise. The flow has two components: (1) internal shearing flow of ice and (2) basal motion, which is sliding along the base of ice sheets, especially when the ice melts at this base. To determine the first component we need to know how "soft" the ice is. By considering the flow velocities at the surface of the northern Greenland Ice Sheet and calculating the stresses that cause the flow, we determined that the ice is effectively softer than is usually assumed. Previous studies indicated that the base of the ice is thawed in large parts (up to about 50%) of the Greenland Ice Sheet. Our study shows that that is probably overestimated, because these studies assumed ice to be harder than it actually is. Our new assessment reduces the area with basal motion and thus melting to about 6-13% in the Greenland study area.

Impurities in polar ice cores have been studied so far mainly for the purpose of reconstructions of past atmospheric aerosol concentrations. However, impurities also critically influence physical properties of the ice matrix itself. To improve the data basis regarding the in-situ form of incorporation and spatial distribution of impurities in ice we used micro-cryo-Raman spectroscopy to identify the location, phase and composition of micrometer-sized inclusions in natural ice samples around the transition from marine isotope stage (MIS) 6 into 5e in the EDML ice core. The combination of Raman results with ice-microsctructure measurements and complementary impurity data provided by the standard analytical methods (IC, CFA, and DEP) allows for a more interdisciplinary approach interconnecting ice core chemistry and ice core physics. While the interglacial samples were dominated by sulfate salts - mainly gypsum, sodium sulfate (possibly thenardite) and iron-potassium sulfate (likely jarosite) - the glacial ice contained high numbers of mineral dust particles - in particular quartz, mica, feldspar, anatase, hematite and carbonaceous particles (black carbon). We cannot confirm cumulation of impurities in the grain boundary network as reported by other studies, neither micro-particles being dragged by migrating grain boundaries nor in form of liquid veins in triple junctions. We argue that mixing of impurities on millimeter scale and chemical reactions are facilitated by the deforming ice matrix. We review possible effects of impurities on physical properties of ice, however the ultimate identification of the deformation agent and the mechanism behind remains challenging.

The effect of grain size on strain rate of ice in the upper 2207 m in the North Greenland Eemian Ice Drilling (NEEM) deep ice core was investigated using a rheological model based on the composite flow law of Goldsby and Kohlstedt (1997, 2001). The grain size was described by both a mean grain size and a grain size distribution, which allowed the strain rate to be calculated using two different model end members: (i) the micro-scale constant stress model where each grain deforms by the same stress and (ii) the micro-scale constant strain rate model where each grain deforms by the same strain rate. The model results show that basal slip accommodated by grain boundary sliding produces almost all of the deformation in the upper 2207 m of the NEEM ice core, while dislocation creep (basal slip accommodated by non-basal slip) hardly contributes to deformation. The difference in calculated strain rate between the two model end members is relatively small. The calculated strain rate in the fine grained glacial ice (1419–2207 m) varies strongly with depth and is about 4–5 times higher than in the coarser grained Holocene ice (0–1419 m). Two peaks in strain rate are predicted at about 1980 and 2100 m of depth. The results from the rheological model and microstructures in the glacial ice indicate that fine grained layers in the glacial ice will act as internal preferential sliding zones in the Greenland ice sheet.

The Ninth Symposium on Polar Science/Ordinary sessions: [OM] Polar meteorology and glaciology, Wed. 5 Dec. / Entrance Hall (1st floor), National Institute of Polar Research

The Ninth Symposium on Polar Science/Ordinary sessions: [OM] Polar meteorology and glaciology, Wed. 5 Dec. / 3F Multipurpose Conference room, National Institute of Polar Research

Disturbances on the centimetre scale in the stratigraphy of the North Greenland Eemian Ice Drilling (NEEM) ice core (North Greenland) can be mapped by an optical line scanner as long as the ice has visual layering, such as, for example, cloudy bands. Different focal depths allow, to a certain extent, a three-dimensional view of the structures. In this study we present a detailed analysis of the visible folds, discuss their characteristics and frequency, and present examples of typical fold structures. We also analyse the structures with regard to the deformation boundary conditions under which they formed. The structures evolve from gentle waves at about 1500 m to overturned z folds with increasing depth. Occasionally, the folding causes significant thickening of layers. Their similar fold shape indicates that they are passive features and are probably not initiated by rheology differences between alternating layers. Layering is heavily disturbed and tracing of single layers is no longer possible below a depth of 2160 m. C axes orientation distributions for the corresponding core sections were analysed, where available, in addition to visual stratigraphy. The data show axial-plane parallel strings of grains with c axis orientations that deviate from that of the matrix, which shows a single maximum fabric at the depth where the folding occurs. Numerical modelling of crystal viscoplastic deformation and dynamic recrystallisation was used to improve the understanding of the formation of the observed structures during deformation. The modelling reproduces the development of bands of grains with a tilted-lattice orientation relative to the single maximum fabric of the matrix, and also the associated local deformation. We conclude from these results that the observed folding can be explained by formation of these tilted-lattice bands.

The effect of grain size on strain rate of ice in the upper 2207 m in the North Greenland Eemian Ice Drilling (NEEM) deep ice core was investigated using a rheological model based on the composite flow law of Goldsby and Kohlstedt (1997, 2001). The grain size was described by both a mean grain size and a grain size distribution, which allowed the strain rate to be calculated using two different model end-members: (i) the microscale constant stress model where each grain deforms by the same stress and (ii) the microscale constant strain rate model where each grain deforms by the same strain rate. The model results predict that grain-size-sensitive flow produces almost all of the deformation in the upper 2207 m of the NEEM ice core, while dislocation creep hardly contributes to deformation. The difference in calculated strain rate between the two model end-members is relatively small. The predicted strain rate in the fine-grained Glacial ice (that is, ice deposited during the last Glacial maximum at depths of 1419 to 2207 m) varies strongly within this depth range and, furthermore, is about 4–5 times higher than in the coarser-grained Holocene ice (0–1419 m). Two peaks in strain rate are predicted at about 1980 and 2100 m depth. The prediction that grain-size-sensitive creep is the fastest process is inconsistent with the microstructures in the Holocene age ice, indicating that the rate of dislocation creep is underestimated in the model. The occurrence of recrystallization processes in the polar ice that did not occur in the experiments may account for this discrepancy. The prediction of the composite flow law model is consistent with microstructures in the Glacial ice, suggesting that fine-grained layers in the Glacial ice may act as internal preferential sliding zones in the Greenland ice sheet.

Impurities control a variety of physical properties of polar ice. Their impact can be observed at all scales – from the microstructure (e.g., grain size and orientation) to the ice sheet flow behavior (e.g., borehole tilting and closure). Most impurities are likely to form micrometer-sized second phase inclusions. It has been suggested that these particles control the grain size of polycrystalline ice by the pinning of grain boundaries (Zener pinning), which should be reflected in the distribution of the inclusions in relation to the grain boundary network. We used an optical microscope to generate high-resolution large-scale maps (3 μm pix−1, 8 × 2 cm2) of the distribution of micro-inclusions in four samples from the EDML (Antarctica) and NEEM (Greenland) polar ice cores. The in-situ positions of more than 5000 μ-inclusions have been determined. A Raman microscope was used to confirm the extrinsic nature of a sample proportion of the mapped inclusions. A superposition of the 2D grain boundary network and μ-inclusion distributions show no significant correlations between grain boundaries and μ-inclusions. In particular, no signs could be found of grain boundaries harvesting μ-inclusions, no evidence of μ-inclusions inhibiting grain boundary migration by slow mode pinning could be detected. Consequences for our understanding of the impurity effect on ice microstructure and rheology are discussed.

An established rift in the Larsen C Ice Shelf, formerly constrained by a suture zone containing marine ice, grew rapidly during 2014 and is likely in the near future to generate the largest calving event since the 1980s and result in a new minimum area for the ice shelf. Here we investigate the recent development of the rift, quantify the projected calving event and, using a numerical model, assess its likely impact on ice shelf stability. We find that the ice front is at risk of becoming unstable when the anticipated calving event occurs.

Microstructures from deep ice cores reflect the dynamic conditions of the drill location as well as the thermodynamic history of the drill site and catchment area in great detail. Ice core parameters (crystal lattice-preferred orientation (LPO), grain size, grain shape), mesostructures (visual stratigraphy) as well as borehole deformation were measured in a deep ice core drilled at Kohnen Station, Dronning Maud Land (DML), Antarctica. These observations are used to characterize the local dynamic setting and its rheological as well as microstructural effects at the EDML ice core drilling site (European Project for Ice Coring in Antarctica in DML). The results suggest a division of the core into five distinct sections, interpreted as the effects of changing deformation boundary conditions from triaxial deformation with horizontal extension to bedrock-parallel shear. Region 1 (uppermost approx. 450m depth) with still small macroscopic strain is dominated by compression of bubbles and strong strain and recrystallization localization. Region 2 (approx. 450-1700m depth) shows a girdle-type LPO with the girdle plane being perpendicular to grain elongations, which indicates triaxial deformation with dominating horizontal extension. In this region (approx. 1000m depth), the first subtle traces of shear deformation are observed in the shape-preferred orientation (SPO) by inclination of the grain elongation. Region 3 (approx. 1700-2030m depth) represents a transitional regime between triaxial deformation and dominance of shear, which becomes apparent in the progression of the girdle to a single maximum LPO and increasing obliqueness of grain elongations. The fully developed single maximum LPO in region 4 (approx. 2030-2385m depth) is an indicator of shear dominance. Region 5 (below approx. 2385m depth) is marked by signs of strong shear, such as strong SPO values of grain elongation and strong kink folding of visual layers. The details of structural observations are compared with results from a numerical ice sheet model (PISM, isotropic) for comparison of strain rate trends predicted from the large-scale geometry of the ice sheet and borehole logging data. This comparison confirms the segmentation into these depth regions and in turn provides a wider view of the ice sheet. (c) 2016 The Author(s) Published by the Royal Society. All rights reserved. This research was funded by HGF grant no. VH-NG-802 to D.J. and I.W., as well as SPP 1158 DFG grant no. BO 1776/12-1 to I.W. S.H.F. acknowledges financial support from the Ramón y Cajal Research Grant RYC-2012-12167 of the Ministry of Economy and Competitiveness of Spain

Antarctica is surrounded by a variety of large, medium and small sized ice shelves, glacier tongues and coastal areas without offshore floating ice masses. We used the mosaic of the Radarsat-1 Antarctica Mapping Project (RAMP) Antarctic Mapping Mission 1 (AMM) to classify the coastline of Antarctica in terms of surface structure patterns close to the calving front. With the aid of an automated edge detection method, complemented by manual control, the surface structures of all ice shelves and glacier tongues around Antarctica were mapped. We found dense and less dense patterns of surface structures unevenly distributed over the ice shelves and ice tongues. Dense surface patterns are frequent on fast flowing ice masses (ice streams), whereas most ice shelves show a dense surface pattern only close to the grounding line. Flow line analyses on ten ice shelves reveal that the time of residence of the ice along a flow path and—associated with it—the healing of surface crevasses can explain the different surface structure distribution close to the grounding line and the calving front on many ice shelves. Based on the surface structures relative to the calving front within a 15 km-wide seaward strip, the ice shelf fronts can be separated into three classes. The resulting map of the classified calving fronts around Antarctica and their description provide a detailed picture of crevasse formation and the observed dominant iceberg shapes.