Converging flow and anisotropy cause large-scale folding in Greenland's ice sheet

Shear margins in upper half of Northeast Greenland Ice Stream were established two millennia ago

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.

Three-dimensional topology dataset of folded radar stratigraphy in northern Greenland

We present a dataset of reconstructed three-dimensional (3D) englacial stratigraphic horizons in northern Greenland. The data cover four different regions representing key ice-dynamic settings in Greenland: (i) the onset of Petermann Glacier, (ii) a region upstream of the 79° North Glacier (Nioghalvfjerdsbræ), near the northern Greenland ice divide, (iii) the onset of the Northeast Greenland Ice Stream (NEGIS) and (iv) a 700 km wide region extending across the central ice divide over the entire northern part of central Greenland. In this paper, we promote the advantages of a 3D perspective of deformed englacial stratigraphy and explain how 3D horizons provide an improved basis for interpreting and reconstructing the ice-dynamic history. The 3D horizons are provided in various formats to allow a wide range of applications and reproducibility of results.

Noise Removal and Feature Extraction in Airborne Radar Sounding Data of Ice Sheets

The airborne ice-penetrating radar (IPR) is an effective method used for ice sheet exploration and is widely applied for detecting the internal structures of ice sheets and for understanding the mechanism of ice flow and the characteristics of the bottom of ice sheets. However, because of the ambient influence and the limitations of the instruments, IPR data are frequently overlaid with noise and interference, which further impedes the extraction of layer features and the interpretation of the physical characteristics of the ice sheet. In this paper, we first applied conventional filtering methods to remove the feature noise and interference in IPR data. Furthermore, machine learning methods were introduced in IPR data processing for noise removal and feature extraction. Inspired by a comparison of the filtering methods and machine learning methods, we propose a fusion method combining both filtering methods and machine-learning-based methods to optimize the feature extraction in IPR data. Field data tests indicated that, under different conditions of IPR data, the application of different methods and strategies can improve the layer feature extraction.

Thermodynamics of a fast-moving Greenlandic outlet glacier revealed by fiber-optic distributed temperature sensing

Measurements of ice temperature provide crucial constraints on ice viscosity and the thermodynamic processes occurring within a glacier. However, such measurements are presently limited by a small number of relatively coarse-spatial-resolution borehole records, especially for ice sheets. Here, we advance our understanding of glacier thermodynamics with an exceptionally high-vertical-resolution (~0.65 m), distributed-fiber-optic temperature-sensing profile from a 1043-m borehole drilled to the base of Sermeq Kujalleq (Store Glacier), Greenland. We report substantial but isolated strain heating within interglacial-phase ice at 208 to 242 m depth together with strongly heterogeneous ice deformation in glacial-phase ice below 889 m. We also observe a high-strain interface between glacial- and interglacial-phase ice and a 73-m-thick temperate basal layer, interpreted as locally formed and important for the glacier's fast motion. These findings demonstrate notable spatial heterogeneity, both vertically and at the catchment scale, in the conditions facilitating the fast motion of marine-terminating glaciers in Greenland.

Microstructures and Fabric Transitions of Natural Ice from the Styx Glacier, Northern Victoria Land, Antarctica

We investigated the microstructures of five ice core samples from the Styx Glacier, northern Victoria Land, Antarctica. Evidence of dynamic recrystallization was found in all samples: those at 50 m mainly by polygonization, and those at 170 m, largely by grain boundary migration. Crystallographic preferred orientations of all analyzed samples (view from the surface) typically showed a single cluster of c-axes normal to the surface. A girdle intersecting the single cluster occurs at 140–170 m with a tight cluster of a-axes normal to the girdle. We interpret the change of crystallographic preferred orientations (CPOs) at <140 m as relating to a combination of vertical compression, and shear on a horizontal plane, and the girdle CPOs at depths >140 m, as the result of horizontal extension. Based on the data obtained from the ground penetrating radar, the underlying bedrock topography of a nunatak could have generated the extensional stress regime in the study area. The results imply changeable stress regimes that may occur during burial as a result of external kinematic controls, such as an appearance of a small peak in the bedrock.

Basal freeze-on generates complex ice-sheet stratigraphy

Large, plume-like internal ice-layer-structures have been observed in radar images from both Antarctica and Greenland, rising from the ice-sheet base to up to half of the ice thickness. Their origins are not yet understood. Here, we simulate their genesis by basal freeze-on using numerical ice-flow modelling and analyse the transient evolution of the emerging ice-plume and the surrounding ice-layer structure as a function of both freeze-on rate and ice flux. We find good agreement between radar observations, modelled ice-plume geometry and internal layer structure, and further show that plume height relates primarily to ice-flux and only secondarily to freeze-on. An in-depth analysis, performed for Northern Greenland of observed spatial plume distribution related to ice flow, basal topography and water availability supports our findings regarding ice flux and suggests freeze-on is controlled by ascending subglacial water flow. Our results imply that widespread basal freeze-on strongly affects ice stratigraphy and consequently ice-core interpretations.

Subsurface ice-sheet radar images reveal large plume-shaped bodies rising from the base, with their origin not yet understood. Here, the authors show that freeze-on of water at the ice-sheet base combined with ice-flux explains the vertical extent, shape and structure of the observed plumes.

Recent Progress in Greenland Ice Sheet Modelling

PURPOSE OF REVIEW

This paper reviews the recent literature on numerical modelling of the dynamics of the Greenland ice sheet with the goal of providing an overview of advancements and to highlight important directions of future research. In particular, the review is focused on large-scale modelling of the ice sheet, including future projections, model parameterisations, paleo applications and coupling with models of other components of the Earth system.

RECENT FINDINGS

Data assimilation techniques have been used to improve the reliability of model simulations of the Greenland ice sheet dynamics, including more accurate initial states, more comprehensive use of remote sensing as well as paleo observations and inclusion of additional physical processes.

SUMMARY

Modellers now leverage the increasing number of high-resolution satellite and air-borne data products to initialise ice sheet models for centennial time-scale simulations, needed for policy relevant sea-level projections. Modelling long-term past and future ice sheet evolution, which requires simplified but adequate representations of the interactions with the other components of the Earth system, has seen a steady improvement. Important developments are underway to include ice sheets in climate models that may lead to routine simulation of the fully coupled Greenland ice sheet-climate system in the coming years.

Dynamic recrystallization during deformation of polycrystalline ice: insights from numerical simulations

The flow of glaciers and polar ice sheets is controlled by the highly anisotropic rheology of ice crystals that have hexagonal symmetry (ice lh). To improve our knowledge of ice sheet dynamics, it is necessary to understand how dynamic recrystallization (DRX) controls ice microstructures and rheology at different boundary conditions that range from pure shear flattening at the top to simple shear near the base of the sheets. We present a series of two-dimensional numerical simulations that couple ice deformation with DRX of various intensities, paying special attention to the effect of boundary conditions. The simulations show how similar orientations of c-axis maxima with respect to the finite deformation direction develop regardless of the amount of DRX and applied boundary conditions. In pure shear this direction is parallel to the maximum compressional stress, while it rotates towards the shear direction in simple shear. This leads to strain hardening and increased activity of non-basal slip systems in pure shear and to strain softening in simple shear. Therefore, it is expected that ice is effectively weaker in the lower parts of the ice sheets than in the upper parts. Strain-rate localization occurs in all simulations, especially in simple shear cases. Recrystallization suppresses localization, which necessitates the activation of hard, non-basal slip systems.This article is part of the themed issue 'Microdynamics of ice'.

Physical analysis of an Antarctic ice core-towards an integration of micro- and macrodynamics of polar ice

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. 450 m depth) with still small macroscopic strain is dominated by compression of bubbles and strong strain and recrystallization localization. Region 2 (approx. 450-1700 m 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. 1000 m 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-2030 m 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-2385 m depth) is an indicator of shear dominance. Region 5 (below approx. 2385 m 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.This article is part of the themed issue 'Microdynamics of ice'.