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.

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.

AI Radar Sensor: Creating Radar Depth Sounder Images Based on Generative Adversarial Network.

Significant resources have been spent in collecting and storing large and heterogeneous radar datasets during expensive Arctic and Antarctic fieldwork. The vast majority of data available is unlabeled, and the labeling process is both time-consuming and expensive. One possible alternative to the labeling process is the use of synthetically generated data with artificial intelligence. Instead of labeling real images, we can generate synthetic data based on arbitrary labels. In this way, training data can be quickly augmented with additional images. In this research, we evaluated the performance of synthetically generated radar images based on modified cycle-consistent adversarial networks. We conducted several experiments to test the quality of the generated radar imagery. We also tested the quality of a state-of-the-art contour detection algorithm on synthetic data and different combinations of real and synthetic data. Our experiments show that synthetic radar images generated by generative adversarial network (GAN) can be used in combination with real images for data augmentation and training of deep neural networks. However, the synthetic images generated by GANs cannot be used solely for training a neural network (training on synthetic and testing on real) as they cannot simulate all of the radar characteristics such as noise or Doppler effects. To the best of our knowledge, this is the first work in creating radar sounder imagery based on generative adversarial network.

A large impact crater beneath Hiawatha Glacier in northwest Greenland.

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

Holocene deceleration of the Greenland Ice Sheet.

Recent peripheral thinning of the Greenland Ice Sheet is partly offset by interior thickening and is overprinted on its poorly constrained Holocene evolution. On the basis of the ice sheet's radiostratigraphy, ice flow in its interior is slower now than the average speed over the past nine millennia. Generally higher Holocene accumulation rates relative to modern estimates can only partially explain this millennial-scale deceleration. The ice sheet's dynamic response to the decreasing proportion of softer ice from the last glacial period and the deglacial collapse of the ice bridge across Nares Strait also contributed to this pattern. Thus, recent interior thickening of the Greenland Ice Sheet is partly an ongoing dynamic response to the last deglaciation that is large enough to affect interpretation of its mass balance from altimetry.

A synthesis of the basal thermal state of the Greenland Ice Sheet.

The basal thermal state of an ice sheet (frozen or thawed) is an important control upon its evolution, dynamics and response to external forcings. However, this state can only be observed directly within sparse boreholes or inferred conclusively from the presence of subglacial lakes. Here we synthesize spatially extensive inferences of the basal thermal state of the Greenland Ice Sheet to better constrain this state. Existing inferences include outputs from the eight thermomechanical ice-flow models included in the SeaRISE effort. New remote-sensing inferences of the basal thermal state are derived from Holocene radiostratigraphy, modern surface velocity and MODIS imagery. Both thermomechanical modeling and remote inferences generally agree that the Northeast Greenland Ice Stream and large portions of the southwestern ice-drainage systems are thawed at the bed, whereas the bed beneath the central ice divides, particularly their west-facing slopes, is frozen. Elsewhere, there is poor agreement regarding the basal thermal state. Both models and remote inferences rarely represent the borehole-observed basal thermal state accurately near NorthGRIP and DYE-3. This synthesis identifies a large portion of the Greenland Ice Sheet (about one third by area) where additional observations would most improve knowledge of its overall basal thermal state.

Radiostratigraphy and age structure of the Greenland Ice Sheet.

Several decades of ice-penetrating radar surveys of the Greenland and Antarctic ice sheets have observed numerous widespread internal reflections. Analysis of this radiostratigraphy has produced valuable insights into ice sheet dynamics and motivates additional mapping of these reflections. Here we present a comprehensive deep radiostratigraphy of the Greenland Ice Sheet from airborne deep ice-penetrating radar data collected over Greenland by The University of Kansas between 1993 and 2013. To map this radiostratigraphy efficiently, we developed new techniques for predicting reflection slope from the phase recorded by coherent radars. When integrated along track, these slope fields predict the radiostratigraphy and simplify semiautomatic reflection tracing. Core-intersecting reflections were dated using synchronized depth-age relationships for six deep ice cores. Additional reflections were dated by matching reflections between transects and by extending reflection-inferred depth-age relationships using the local effective vertical strain rate. The oldest reflections, dating to the Eemian period, are found mostly in the northern part of the ice sheet. Within the onset regions of several fast-flowing outlet glaciers and ice streams, reflections typically do not conform to the bed topography. Disrupted radiostratigraphy is also observed in a region north of the Northeast Greenland Ice Stream that is not presently flowing rapidly. Dated reflections are used to generate a gridded age volume for most of the ice sheet and also to determine the depths of key climate transitions that were not observed directly. This radiostratigraphy provides a new constraint on the dynamics and history of the Greenland Ice Sheet. KEY POINTS: Phase information predicts reflection slope and simplifies reflection tracingReflections can be dated away from ice cores using a simple ice flow modelRadiostratigraphy is often disrupted near the onset of fast ice flow.

A reconciled estimate of ice-sheet mass balance.

We combined an ensemble of satellite altimetry, interferometry, and gravimetry data sets using common geographical regions, time intervals, and models of surface mass balance and glacial isostatic adjustment to estimate the mass balance of Earth's polar ice sheets. We find that there is good agreement between different satellite methods--especially in Greenland and West Antarctica--and that combining satellite data sets leads to greater certainty. Between 1992 and 2011, the ice sheets of Greenland, East Antarctica, West Antarctica, and the Antarctic Peninsula changed in mass by -142 ± 49, +14 ± 43, -65 ± 26, and -20 ± 14 gigatonnes year(-1), respectively. Since 1992, the polar ice sheets have contributed, on average, 0.59 ± 0.20 millimeter year(-1) to the rate of global sea-level rise.