skip to main content


Title: Estimating Greenland tidewater glacier retreat driven by submarine melting
Abstract. The effect of the North Atlantic Ocean on the Greenland Ice Sheet through submarine melting of Greenland's tidewater glacier calving fronts is thought to be a key driver of widespread glacier retreat, dynamic mass loss and sea level contribution from the ice sheet. Despite its critical importance, problems of process complexity and scale hinder efforts to represent the influence of submarine melting in ice-sheet-scale models. Here we propose parameterizing tidewater glacier terminus position as a simple linear function of submarine melting, with submarine melting in turn estimated as a function of subglacial discharge and ocean temperature. The relationship is tested, calibrated and validated using datasets of terminus position, subglacial discharge and ocean temperature covering the full ice sheet and surrounding ocean from the period 1960–2018. We demonstrate a statistically significant link between multi-decadal tidewater glacier terminus position change and submarine melting and show that the proposed parameterization has predictive power when considering a population of glaciers. An illustrative 21st century projection is considered, suggesting that tidewater glaciers in Greenland will undergo little further retreat in a low-emission RCP2.6 scenario. In contrast, a high-emission RCP8.5 scenario results in a median retreat of 4.2 km, with a quarter of tidewater glaciers experiencing retreat exceeding 10 km. Our study provides a long-term and ice-sheet-wide assessment of the sensitivity of tidewater glaciers to submarine melting and proposes a practical and empirically validated means of incorporating ocean forcing into models of the Greenland ice sheet.  more » « less
Award ID(s):
1916566
NSF-PAR ID:
10158508
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
The Cryosphere
Volume:
13
Issue:
9
ISSN:
1994-0424
Page Range / eLocation ID:
2489 to 2509
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Changes in ocean temperature and salinity are expected to be an important determinant of the Greenland ice sheet's future sea level contribution. Yet, simulating the impact of these changes in continental-scale ice sheet models remains challenging due to the small scale of key physics, such as fjord circulation and plume dynamics, and poor understanding of critical processes, such as calving and submarine melting. Here we present the ocean forcing strategy for Greenland ice sheet models taking part in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), the primary community effort to provide 21st century sea level projections for the Intergovernmental Panel on Climate Change Sixth Assessment Report. Beginning from global atmosphere–ocean general circulation models, we describe two complementary approaches to provide ocean boundary conditions for Greenland ice sheet models, termed the “retreat” and “submarine melt” implementations. The retreat implementation parameterises glacier retreat as a function of projected subglacial discharge and ocean thermal forcing, is designed to be implementable by all ice sheet models and results in retreat of around 1 and 15 km by 2100 in RCP2.6 and 8.5 scenarios, respectively. The submarine melt implementation provides estimated submarine melting only, leaving the ice sheet model to solve for the resulting calving and glacier retreat and suggests submarine melt rates will change little under RCP2.6 but will approximately triple by 2100 under RCP8.5. Both implementations have necessarily made use of simplifying assumptions and poorly constrained parameterisations and, as such, further research on submarine melting, calving and fjord–shelf exchange should remain a priority. Nevertheless, the presented framework will allow an ensemble of Greenland ice sheet models to be systematically and consistently forced by the ocean for the first time and should result in a significant improvement in projections of the Greenland ice sheet's contribution to future sea level change. 
    more » « less
  2. Abstract

    Submarine melting has been implicated in the accelerated retreat of marine‐terminating glaciers globally. Energetic ocean flows, such as subglacial discharge plumes, are known to enhance submarine melting in their immediate vicinity. Using observations and a large eddy simulation, we demonstrate that discharge plumes emit high‐frequency internal gravity waves that propagate along glacier termini and transfer energy to distant regions of the terminus. Our analysis of wave characteristics and their correlation with subglacial discharge forcing suggest that they derive their energy from turbulent motions within the discharge plume and its surface outflow. Accounting for the near‐terminus velocities associated with these waves increases predicted melt rates by up to 70%. This may help to explain known discrepancies between observed melt rates and theoretical predictions. Because the dynamical ingredients—a buoyant plume rising through a stratified ocean—are common to many tidewater glacier systems, such internal waves are likely to be widespread.

     
    more » « less
  3. Abstract

    Fjords are conduits for heat and mass exchange between tidewater glaciers and the coastal ocean, and thus regulate near‐glacier water properties and submarine melting of glaciers. Entrainment into subglacial discharge plumes is a primary driver of seasonal glacial fjord circulation; however, outflowing plumes may continue to influence circulation after reaching neutral buoyancy through the sill‐driven mixing and recycling, or reflux, of glacial freshwater. Despite its importance in non‐glacial fjords, no framework exists for how freshwater reflux may affect circulation in glacial fjords, where strong buoyancy forcing is also present. Here, we pair a suite of hydrographic observations measured throughout 2016–2017 in LeConte Bay, Alaska, with a three‐dimensional numerical model of the fjord to quantify sill‐driven reflux of glacial freshwater, and determine its influence on glacial fjord circulation. When paired with subglacial discharge plume‐driven buoyancy forcing, sill‐generated mixing drives distinct seasonal circulation regimes that differ greatly in their ability to transport heat to the glacier terminus. During the summer, 53%–72% of the surface outflow is refluxed at the fjord's shallow entrance sill and is subsequently re‐entrained into the subglacial discharge plume at the fjord head. As a result, near‐terminus water properties are heavily influenced by mixing at the entrance sill, and circulation is altered to draw warm, modified external surface water to the glacier grounding line at 200 m depth. This circulatory cell does not exist in the winter when freshwater reflux is minimal. Similar seasonal behavior may exist at other glacial fjords throughout Southeast Alaska, Patagonia, Greenland, and elsewhere.

     
    more » « less
  4. Abstract

    The uncertain response of marine terminating outlet glaciers to climate change at time scales beyond short-term observation limits models of future sea level rise. At temperate tidewater margins, abundant subglacial meltwater forms morainal banks (marine shoals) or ice-contact deltas that reduce water depth, stabilizing grounding lines and slowing or reversing glacial retreat. Here we present a radiocarbon-dated record from Integrated Ocean Drilling Program (IODP) Site U1421 that tracks the terminus of the largest Alaskan Cordilleran Ice Sheet outlet glacier during Last Glacial Maximum climate transitions. Sedimentation rates, ice-rafted debris, and microfossil and biogeochemical proxies, show repeated abrupt collapses and slow advances typical of the tidewater glacier cycle observed in modern systems. When global sea level rise exceeded the local rate of bank building, the cycle of readvances stopped leading to irreversible retreat. These results support theory that suggests sediment dynamics can control tidewater terminus position on an open shelf under temperate conditions delaying climate-driven retreat.

     
    more » « less
  5. Abstract

    At tidewater glacier termini, ocean‐glacier interactions hinge on two sources of freshwater—submarine melt and subglacial discharge—yet these freshwater fluxes are often unconstrained in their magnitude, seasonality, and relationship. With measurements of ocean velocity, temperature and salinity, fjord budgets can be evaluated to partition the freshwater flux into submarine melt and subglacial discharge. We apply these methods to calculate the freshwater fluxes at LeConte Glacier, Alaska, across a wide range of oceanic and atmospheric conditions during six surveys in 2016–2018. We compare these ocean‐derived fluxes with an estimate of subglacial discharge from a surface mass balance model and with estimates of submarine melt from multibeam sonar and autonomous kayaks, finding relatively good agreement between these independent estimates. Across spring, summer, and fall, the relationship between subglacial discharge and submarine melt follows a scaling law predicted by standard theory (melt ∼ discharge1/3), although the total magnitude of melt is an order of magnitude larger than theoretical estimates. Subglacial discharge is the dominant driver of variability in melt, while the dependence of melt on fjord properties is not discernible. A comparison of oceanic budgets with glacier records indicates that submarine melt removes 33%–49% of the ice flux into the terminus across spring, summer, and fall periods. Thus, melt is a significant component of the glacier's mass balance, and we find that melt correlates with seasonal retreat; however, melt does not appear to directly amplify calving.

     
    more » « less