skip to main content


Title: Twenty-first century ocean forcing of the Greenland ice sheet for modelling of sea level contribution
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
Award ID(s):
1916566
NSF-PAR ID:
10158504
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
The Cryosphere
Volume:
14
Issue:
3
ISSN:
1994-0424
Page Range / eLocation ID:
985 to 1008
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract

    Recent acceleration of Greenland's ocean‐terminating glaciers has substantially amplified the ice sheet's contribution to global sea level. Increased oceanic melting of these tidewater glaciers is widely cited as the likely trigger, and is thought to be highest within vigorous plumes driven by freshwater drainage from beneath glaciers. Yet melting of the larger part of calving fronts outside of plumes remains largely unstudied. Here we combine ocean observations collected within 100 m of a tidewater glacier with a numerical model to show that unlike previously assumed, plumes drive an energetic fjord‐wide circulation which enhances melting along the entire calving front. Compared to estimates of melting within plumes alone, this fjord‐wide circulation effectively doubles the glacier‐wide melt rate, and through shaping the calving front has a potential dynamic impact on calving. Our results suggest that melting driven by fjord‐scale circulation should be considered in process‐based projections of Greenland's sea level contribution.

     
    more » « less
  3. Abstract. Numerical simulations of the Greenland Ice Sheet (GrIS) over geologictimescales can greatly improve our knowledge of the critical factors drivingGrIS demise during climatically warm periods, which has clear relevance forbetter predicting GrIS behavior over the upcoming centuries. To assess thefidelity of these modeling efforts, however, observational constraints ofpast ice sheet change are needed. Across southwestern Greenland, geologicrecords detail Holocene ice retreat across both terrestrial-based and marine-terminating environments, providing an ideal opportunity to rigorouslybenchmark model simulations against geologic reconstructions of ice sheetchange. Here, we present regional ice sheet modeling results using theIce-sheet and Sea-level System Model (ISSM) of Holocene ice sheet historyacross an extensive fjord region in southwestern Greenland covering thelandscape around the Kangiata Nunaata Sermia (KNS) glacier and extendingoutward along the 200 km Nuup Kangerula (Godthåbsfjord). Oursimulations, forced by reconstructions of Holocene climate and recentlyimplemented calving laws, assess the sensitivity of ice retreat across theKNS region to atmospheric and oceanic forcing. Our simulations reveal thatthe geologically reconstructed ice retreat across the terrestrial landscapein the study area was likely driven by fluctuations in surface mass balancein response to Early Holocene warming – and was likely not influencedsignificantly by the response of adjacent outlet glaciers to calving andocean-induced melting. The impact of ice calving within fjords, however,plays a significant role by enhancing ice discharge at the terminus, leadingto interior thinning up to the ice divide that is consistent withreconstructed magnitudes of Early Holocene ice thinning. Our results,benchmarked against geologic constraints of past ice-margin change, suggestthat while calving did not strongly influence Holocene ice-margin migrationacross terrestrial portions of the KNS forefield, it strongly impactedregional mass loss. While these results imply that the implementation andresolution of ice calving in paleo-ice-flow models is important towardsmaking more robust estimations of past ice mass change, they also illustratethe importance these processes have on contemporary and future long-term icemass change across similar fjord-dominated regions of the GrIS. 
    more » « less
  4. Abstract

    The mass loss of the Greenland Ice Sheet is nearly equally partitioned between a decrease in surface mass balance from enhanced surface melt and an increase in ice dynamics from the acceleration and retreat of its marine-terminating glaciers. Much uncertainty remains in the future mass loss of the Greenland Ice Sheet due to the challenges of capturing the ice dynamic response to climate change in numerical models. Here, we estimate the sea level contribution of the Greenland Ice Sheet over the 21st century using an ice-sheet wide, high-resolution, ice-ocean numerical model that includes surface mass balance forcing, thermal forcing from the ocean, and iceberg calving dynamics. The model is calibrated with ice front observations from the past eleven years to capture the recent evolution of marine-terminating glaciers. Under a business as usual scenario, we find that northwest and central west Greenland glaciers will contribute more mass loss than other regions due to ice front retreat and ice flow acceleration. By the end of century, ice discharge from marine-terminating glaciers will contribute 50 ± 20% of the total mass loss, or twice as much as previously estimated although the contribution from the surface mass balance increases towards the end of the century.

     
    more » « less
  5. Abstract

    Thwaites Glacier (TG) plays an important role in future sea-level rise (SLR) contribution from the West Antarctic Ice Sheet. Recent observations show that TG is losing mass, and its grounding zone is retreating. Previous modeling has produced a wide range of results concerning whether, when, and how rapidly further retreat will occur under continued warming. These differences arise at least in part from ill-constrained processes, including friction from the bed, and future atmosphere and ocean forcing affecting ice-shelf and grounding-zone buttressing. Here, we apply the Ice Sheet and Sea-level System Model (ISSM) with a range of specifications of basal sliding behavior in response to varying ocean forcing. We find that basin-wide bed character strongly affects TG's response to sub-shelf melt by modulating how changes in driving stress are balanced by the bed as the glacier responds to external forcing. Resulting differences in dynamic thinning patterns alter modeled grounding-line retreat across Thwaites' catchment, affecting both modeled rates and magnitudes of SLR contribution from this critical sector of the ice sheet. Bed character introduces large uncertainties in projections of TG under equal external forcing, pointing to this as a crucial constraint needed in predictive models of West Antarctica.

     
    more » « less