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.

Crystal orientation fabric anisotropy causes directional hardening of the Northeast Greenland Ice Stream.

The dynamic mass loss of ice sheets constitutes one of the biggest uncertainties in projections of ice-sheet evolution. One central, understudied aspect of ice flow is how the bulk orientation of the crystal orientation fabric translates to the mechanical anisotropy of ice. Here we show the spatial distribution of the depth-averaged horizontal anisotropy and corresponding directional flow-enhancement factors covering a large area of the Northeast Greenland Ice Stream onset. Our results are based on airborne and ground-based radar surveys, ice-core observations, and numerical ice-flow modelling. They show a strong spatial variability of the horizontal anisotropy and a rapid crystal reorganisation on the order of hundreds of years coinciding with the ice-stream geometry. Compared to isotropic ice, parts of the ice stream are found to be more than one order of magnitude harder for along-flow extension/compression while the shear margins are potentially softened by a factor of two for horizontal-shear deformation.

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

The increasing catalogue of high-quality ice-penetrating radar data provides a unique insight in the internal layering architecture of the Greenland ice sheet. The stratigraphy, an indicator of past deformation, highlights irregularities in ice flow and reveals large perturbations without obvious links to bedrock shape. In this work, to establish a new conceptual model for the formation process, we analysed the radar data at the onset of the Petermann Glacier, North Greenland, and created a three-dimensional model of several distinct stratigraphic layers. We demonstrate that the dominant structures are cylindrical folds sub-parallel to the ice flow. By numerical modelling, we show that these folds can be formed by lateral compression of mechanically anisotropic ice, while a general viscosity contrast between layers would not lead to folding for the same boundary conditions. We conclude that the folds primarily form by converging flow as the mechanically anisotropic ice is channelled towards the glacier.