Greenland Ice Sheet Dynamics Research Articles

Holocene thinning in central Greenland controlled by the Northeast Greenland Ice Stream

Ice-core records from the interior of the Greenland ice sheet suggest widespread thinning during the Holocene. However, the recurring underestimation of this thinning in numerical models raises concerns about both the veracity of such reconstructions and the reliability of glaciological models. Recent work suggests the 8000-year-old Northeast Greenland Ice Stream (NEGIS), including a now-extinct northern tributary, may have been an early influence on Greenland ice-sheet dynamics. Yet, the inaccurate reproduction of NEGIS-like dynamics in most models hampers investigation of whether this feature played a role in Holocene ice-sheet thinning. Here we show that grounding-line retreat in northeast Greenland triggers elevation changes at the northern summit via ice-dynamic effects modulated by the paleo NEGIS system. In our simulations, fast ice-stream flow caused by transiently imposed reduced basal shear stress following the northeast retreat explains 55% ( ± 18%) of the estimated ice thinning, showing that ice-stream dynamics is one of the main drivers of the NGRIP Holocene surface elevation drop. Our findings show that the ice-flow in northeast Greenland plays a large role in ice-surface elevation changes in central Greenland.

Supraglacial lake drainage at a fast-flowing Greenlandic outlet glacier

Supraglacial lake drainage events influence Greenland Ice Sheet dynamics on hourly to interannual timescales. However, direct observations are rare, and, to date, no in situ studies exist from fast-flowing sectors of the ice sheet. Here, we present observations of a rapid lake drainage event at Store Glacier, west Greenland, in 2018. The drainage event transported 4.8 × 106 m3 of meltwater to the glacier bed in ∼5 h, reducing the lake to a third of its original volume. During drainage, the local ice surface rose by 0.55 m, and surface velocity increased from 2.0 m⋅d-1 to 5.3 m⋅d-1 Dynamic responses were greatest ∼4 km downstream from the lake, which we interpret as an area of transient water storage constrained by basal topography. Drainage initiated, without any precursory trigger, when the lake expanded and reactivated a preexisting fracture that had been responsible for a drainage event 1 y earlier. Since formation, this fracture had advected ∼500 m from the lake's deepest point, meaning the lake did not fully drain. Partial drainage events have previously been assumed to occur slowly via lake overtopping, with a comparatively small dynamic influence. In contrast, our findings show that partial drainage events can be caused by hydrofracture, producing new hydrological connections that continue to concentrate the supply of surface meltwater to the bed of the ice sheet throughout the melt season. Our findings therefore indicate that the quantity and resultant dynamic influence of rapid lake drainages are likely being underestimated.

A note on digital elevation model smoothing and driving stresses

Ice-flow fields, including the driving stress, provide important information on the current state and evolution of Antarctic and Greenland ice-sheet dynamics. However, computation of flow fields from continent-scale DEMs requires the use of smoothing functions and scales, the choice of which can be ad hoc. This study evaluates smoothing functions and scales for robust calculations of driving stress from Antarctic DEMs. Our approach compares a variety of filters and scales for their capacity to minimize the residual between predicted and observed flow direction fields. We find that a spatially varying triangular filter with a width of 8–10 ice thicknesses provides the closest match between the observed and predicted flow direction fields. We use the predicted flow direction fields to highlight artefacts in observed Antarctic velocities, demonstrating that comparison of multiple observational data sets has utility for quality control of continent-scale data sets.

Seasonal Evolution of the Subglacial Hydrologic System Modified by Supraglacial Lake Drainage in Western Greenland

AbstractThe impact of summer surface melt on Greenland Ice Sheet dynamics is modulated by the state of the subglacial hydrologic system. Studies of ice motion indicate that efficiency of the subglacial system increases over the melt season, decreasing the sensitivity of ice motion to surface melt inputs. However, the behavior of the subglacial hydrologic system is complex and some characteristics are still poorly constrained. Here we investigate the coevolution of subglacial hydrology and ice motion in the Pâkitsoq region of western Greenland during the 2011 melt season. We analyze measurements from 11 Global Positioning System stations, from which we derive ice velocity, longitudinal strain rates, and basal uplift, alongside observations of surface ablation and supraglacial lake drainages. We observe ice acceleration after the onset of local surface melting, followed by gradual ice deceleration, consistent with increasing subglacial efficiency. In the study area, supraglacial lake drainages cooccur with a change in regional strain rate patterns and ice deceleration, suggesting that lake drainages contribute to rapid subglacial reorganization. At lower ice surface elevations (below ~900 m above sea level), ice motion is correlated with both total basal uplift and its rate of change, while at higher elevations (~900–1,100 m above sea level), ice motion correlated only with the basal uplift rate. This pattern suggests that continued cavity growth or subglacial sediment dynamics may be important in the apparent increase in subglacial drainage efficiency at higher elevations in the ablation zone. Our results further suggest that transient subglacial behavior is important in the seasonal evolution of ice motion.

Shrubs shed light on 20th century Greenland Ice Sheet melting

The Greenland Ice Sheet (GrIS) is a key element of the global climate system and thus knowledge about its melting in the past is desirable. However, GrIS‐melt records only date back until 1979 and climate data based reconstructions covering the 20th century differ with respect to absolute values. To extend our knowledge about the Greenland Ice Sheet we examined the potential of shrub ring‐widths and wood‐anatomy as proxies for GrIS‐melt. We found significant correlations between shrub cell‐wall thickness and regional melt derived from passive microwave satellite brightness for 7% of the total GrIS area. A respective transfer function calibrated over the period 1979 to 2007 successfully passed model calibration‐verification tests and explained 42% of GrIS‐melt variability. Consequently, the first GrIS‐melt reconstruction based on shrub wood‐anatomy covering the period 1909 to 2012 is presented and compared against two temperature‐based reconstructions. For the period prior to 1930 the new reconstruction contrasts with existing literature but generally confirms that most recent record melt rates are amongst the highest since the early 20th century. We discuss the sensitivity of shrubs to several influencing factors besides summer temperature as possible reason for the observed differences and highlight the potential of using shrubs as multi‐parameter proxies within a network to increase our knowledge about 20th century Greenland Ice Sheet dynamics.

Laser altimetry reveals complex pattern of Greenland Ice Sheet dynamics

We present a new record of ice thickness change, reconstructed at nearly 100,000 sites on the Greenland Ice Sheet (GrIS) from laser altimetry measurements spanning the period 1993-2012, partitioned into changes due to surface mass balance (SMB) and ice dynamics. We estimate a mean annual GrIS mass loss of 243 ± 18 Gt ⋅ y(-1), equivalent to 0.68 mm ⋅ y(-1) sea level rise (SLR) for 2003-2009. Dynamic thinning contributed 48%, with the largest rates occurring in 2004-2006, followed by a gradual decrease balanced by accelerating SMB loss. The spatial pattern of dynamic mass loss changed over this time as dynamic thinning rapidly decreased in southeast Greenland but slowly increased in the southwest, north, and northeast regions. Most outlet glaciers have been thinning during the last two decades, interrupted by episodes of decreasing thinning or even thickening. Dynamics of the major outlet glaciers dominated the mass loss from larger drainage basins, and simultaneous changes over distances up to 500 km are detected, indicating climate control. However, the intricate spatiotemporal pattern of dynamic thickness change suggests that, regardless of the forcing responsible for initial glacier acceleration and thinning, the response of individual glaciers is modulated by local conditions. Recent projections of dynamic contributions from the entire GrIS to SLR have been based on the extrapolation of four major outlet glaciers. Considering the observed complexity, we question how well these four glaciers represent all of Greenland's outlet glaciers.

New constraints on Greenland ice sheet dynamics during the last glacial cycle: Evidence from the Uummannaq ice stream system

AbstractThis paper presents the first assessment of the Uummannaq ice stream system (UISS) in West Greenland. The UISS drained ~6% of the Greenland ice sheet (GrIS) at the Last Glacial Maximum (LGM). The onset of the UISS is a function of a convergent network of fjords which feed a geologically controlled trough system running offshore to the shelf break. Mapping, cosmogenic radiogenic nuclide (CRN) dating, and model output reveal that glacially scoured surfaces up to 1266 m above sea level (asl) in fjord‐head areas were produced by warm‐based ice moving offshore during the LGM, with the elevation of warm‐based ice dropping westwards to ~700 m asl as the ice stream trunk zone developed. Marginal plateaux with allochthonous blockfields suggest that warm‐based ice produced till and erratics up to ~1200 m asl, but CRN ages and weathering pits suggest this was pre‐LGM, with only cold‐based ice operating during the LGM. Deglaciation began on the outer shelf at ~14.8 cal. kyrs B.P., with Ubekendt Ejland becoming ice free at ~12.4 ka. The UISS then collapsed with over 100 km of retreat by ~11.4 ka–10.8 cal. kyrs B.P., a rapid and complex response to bathymetric deepening, trough widening, and sea‐level rise coinciding with rapidly increasing air temperatures and solar radiation, but which occurred prior to ocean warming at ~8.4 cal. kyrs B.P. Local fjord constriction temporarily stabilized the unzipped UISS margins at the start of the Holocene before ice retreat inland of the current margin at ~8.7 ka.

Committed sea-level rise for the next century from Greenland ice sheet dynamics during the past decade

We use a three-dimensional, higher-order ice flow model and a realistic initial condition to simulate dynamic perturbations to the Greenland ice sheet during the last decade and to assess their contribution to sea level by 2100. Starting from our initial condition, we apply a time series of observationally constrained dynamic perturbations at the marine termini of Greenland's three largest outlet glaciers, Jakobshavn Isbræ, Helheim Glacier, and Kangerdlugssuaq Glacier. The initial and long-term diffusive thinning within each glacier catchment is then integrated spatially and temporally to calculate a minimum sea-level contribution of approximately 1 ± 0.4 mm from these three glaciers by 2100. Based on scaling arguments, we extend our modeling to all of Greenland and estimate a minimum dynamic sea-level contribution of approximately 6 ± 2 mm by 2100. This estimate of committed sea-level rise is a minimum because it ignores mass loss due to future changes in ice sheet dynamics or surface mass balance. Importantly, > 75% of this value is from the long-term, diffusive response of the ice sheet, suggesting that the majority of sea-level rise from Greenland dynamics during the past decade is yet to come. Assuming similar and recurring forcing in future decades and a self-similar ice dynamical response, we estimate an upper bound of 45 mm of sea-level rise from Greenland dynamics by 2100. These estimates are constrained by recent observations of dynamic mass loss in Greenland and by realistic model behavior that accounts for both the long-term cumulative mass loss and its decay following episodic boundary forcing.

Ice stream influence on West Greenland Ice Sheet dynamics during the Last Glacial Maximum

AbstractThis paper investigates the processes governing bedrock bedform evolution in ice sheet and ice stream areas in central West Greenland, and explores the evidence for a cross‐shelf ice stream at the Last Glacial Maximum (LGM). To the east of Sisimiut the formation of streamlined bedforms with high elongation ratios and high bedform density has been controlled by geological structure and topography in slow‐flowing ice sheet areas. At the coast, the effects of regional flow convergence, caused by coastal fjord orientation, routed ice into the Sisimiut/Itilleq area where it formed an ice stream onset zone. This funnelled ice into an offshore trough (Holsteinsborg Dyb), resulting in a southwesterly regional ice flow direction and the formation of a topographically routed ice stream (Holsteinsborg Isbrae). To the south of this, striae and bedform evidence show that local valley glaciers initially flowed east to west across the coast, but were later redirected by the Itilleq Fjord ice which turned southwestward due to diffluent flow and deflection by Holsteinsborg Isbrae. Roches moutonnées in this area have low elongation ratios and high bedform density, but do not provide unequivocal support for ice streaming, as they are a product of both bedrock structure and changes in ice flow direction, rather than enhanced flow velocities. Cosmogenic surface exposure ages limit maximum ice sheet surface elevation to ca. 755–810 m above sea level in this region. Such ice thickness enabled Holsteinsborg Isbrae to reach the mid/outer continental shelf during the LGM, and to contribute to the formation of a trough mouth fan and the Outer Hellefisk moraines. Initial deglaciation across this region was driven by rising sea level and increasing air temperatures prior to the Bølling Interstadial at ca. 14.5 cal. ka BP. Between 12 and 10 cal. ka BP both increased air and ocean temperatures post the Younger Dryas, and peak sea‐level rise up to the marine limit, caused accelerated thinning and marginal retreat through calving, although dating evidence suggests ice streams remained along the inner shelf/coast boundary until at least ca. 10 cal. ka BP, their longevity maintained by increased ice thickness and ice discharge. Copyright © 2009 John Wiley & Sons, Ltd.