Contribution of the Greenland Ice Sheet to sea level over the next millennium.

The Greenland Ice Sheet holds 7.2 m of sea level equivalent and in recent decades, rising temperatures have led to accelerated mass loss. Current ice margin recession is led by the retreat of outlet glaciers, large rivers of ice ending in narrow fjords that drain the interior. We pair an outlet glacier-resolving ice sheet model with a comprehensive uncertainty quantification to estimate Greenland's contribution to sea level over the next millennium. We find that Greenland could contribute 5 to 33 cm to sea level by 2100, with discharge from outlet glaciers contributing 8 to 45% of total mass loss. Our analysis shows that uncertainties in projecting mass loss are dominated by uncertainties in climate scenarios and surface processes, whereas uncertainties in calving and frontal melt play a minor role. We project that Greenland will very likely become ice free within a millennium without substantial reductions in greenhouse gas emissions.

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.

Complex Greenland outlet glacier flow captured.

The Greenland Ice Sheet is losing mass at an accelerating rate due to increased surface melt and flow acceleration in outlet glaciers. Quantifying future dynamic contributions to sea level requires accurate portrayal of outlet glaciers in ice sheet simulations, but to date poor knowledge of subglacial topography and limited model resolution have prevented reproduction of complex spatial patterns of outlet flow. Here we combine a high-resolution ice-sheet model coupled to uniformly applied models of subglacial hydrology and basal sliding, and a new subglacial topography data set to simulate the flow of the Greenland Ice Sheet. Flow patterns of many outlet glaciers are well captured, illustrating fundamental commonalities in outlet glacier flow and highlighting the importance of efforts to map subglacial topography. Success in reproducing present day flow patterns shows the potential for prognostic modelling of ice sheets without the need for spatially varying parameters with uncertain time evolution.

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.

Large subglacial lakes in East Antarctica at the onset of fast-flowing ice streams.

Water plays a crucial role in ice-sheet stability and the onset of ice streams. Subglacial lake water moves between lakes and rapidly drains, causing catastrophic floods. The exact mechanisms by which subglacial lakes influence ice-sheet dynamics are unknown, however, and large subglacial lakes have not been closely associated with rapidly flowing ice streams. Here we use satellite imagery and ice-surface elevations to identify a region of subglacial lakes, similar in total area to Lake Vostok, at the onset region of the Recovery Glacier ice stream in East Antarctica and predicted by ice-sheet models. We define four lakes through extensive, flat, featureless regions of ice surface bounded by upstream troughs and downstream ridges. Using ice velocities determined using interferometric synthetic aperture radar (InSAR), we find the onset of rapid flow (moving at 20 to 30 m yr(-1)) of the tributaries to the Recovery Glacier ice stream in a 280-km-wide segment at the downslope margins of these four subglacial lakes. We conclude that the subglacial lakes initiate and maintain rapid ice flow through either active modification of the basal thermal regime of the ice sheet by lake accretion or through scouring bedrock channels in periodic drainage events. We suggest that the role of subglacial lakes needs to be considered in ice-sheet mass balance assessments.