Search for: All records

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. 

Abstract. The ice sheet model intercomparison project for CMIP6 (ISMIP6) effort brings together the ice sheet and climate modeling communities to gain understanding of the ice sheet contribution to sea level rise. ISMIP6 conducts stand-alone ice sheet experiments that use space- and time-varying forcing derived from atmosphere–ocean coupled global climate models (AOGCMs) to reflect plausible trajectories for climate projections. The goal of this study is to recommend a subset of CMIP5 AOGCMs (three core and three targeted) to produce forcing for ISMIP6 stand-alone ice sheet simulations, based on (i) their representation of current climate near Antarctica and Greenland relative to observations and (ii) their ability to sample a diversity of projected atmosphere and ocean changes over the 21st century. The selection is performed separately for Greenland and Antarctica. Model evaluation over the historical period focuses on variables used to generate ice sheet forcing. For stage (i), we combine metrics of atmosphere and surface ocean state (annual- and seasonal-mean variables over large spatial domains) with metrics of time-mean subsurface ocean temperature biases averaged over sectors of the continental shelf. For stage (ii), we maximize the diversity of climate projections among the best-performing models. Model selection is also constrained by technical limitations, such as availability of required data from RCP2.6 and RCP8.5 projections. The selected top three CMIP5 climate models are CCSM4, MIROC-ESM-CHEM, and NorESM1-M for Antarctica and HadGEM2-ES, MIROC5, and NorESM1-M for Greenland. This model selection was designed specifically for ISMIP6 but can be adapted for other applications. 

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.