ALBERT

Description: Projection of the contribution of ice sheets to sea level change as part of the Coupled Model Intercomparison Project Phase 6 (CMIP6) takes the form of simulations from coupled ice sheet–climate models and stand-alone ice sheet models, overseen by the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). This paper describes the experimental setup for process-based sea level change projections to be performed with stand-alone Greenland and Antarctic ice sheet models in the context of ISMIP6. The ISMIP6 protocol relies on a suite of polar atmospheric and oceanic CMIP-based forcing for ice sheet models, in order to explore the uncertainty in projected sea level change due to future emissions scenarios, CMIP models, ice sheet models, and parameterizations for ice–ocean interactions. We describe here the approach taken for defining the suite of ISMIP6 stand-alone ice sheet simulations, document the experimental framework and implementation, and present an overview of the ISMIP6 forcing to be used by participating ice sheet modeling groups.

Description: Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between −7.8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between −6.1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

Description: Climate model projections have previously been used to compute ice shelf basal melt rates in ice sheet models, but the strategies employed – e.g., ocean input, parameterization, calibration technique, and corrections – have varied widely and are often ad hoc. Here, a methodology is proposed for the calculation of circum-Antarctic basal melt rates for floating ice, based on climate models, that is suitable for ISMIP6, the Ice Sheet Model Intercomparison Project for CMIP6 (6th Coupled Model Intercomparison Project). The past and future evolution of ocean temperature and salinity is derived from a climate model by estimating anomalies with respect to the modern day, which are added to a present-day climatology constructed from existing observational datasets. Temperature and salinity are extrapolated to any position potentially occupied by a simulated ice shelf. A simple formulation is proposed for a basal melt parameterization in ISMIP6, constrained by the observed temperature climatology, with a quadratic dependency on either the nonlocal or local thermal forcing. Two calibration methods are proposed: (1) based on the mean Antarctic melt rate (MeanAnt) and (2) based on melt rates near Pine Island's deep grounding line (PIGL). Future Antarctic mean melt rates are an order of magnitude greater in PIGL than in MeanAnt. The PIGL calibration and the local parameterization result in more realistic melt rates near grounding lines. PIGL is also more consistent with observations of interannual melt rate variability underneath Pine Island and Dotson ice shelves. This work stresses the need for more physics and less calibration in the parameterizations and for more observations of hydrographic properties and melt rates at interannual and decadal timescales.

Description: Fram Strait, the deepest gateway to the Arctic Ocean, is strongly influenced by eddy dynamics. Here we analyse the output from two eddy-resolving models (ROMS – Regional Ocean Modeling System; FESOM – Finite-Element Sea-ice Ocean Model) with around 1 km mesh resolution in Fram Strait, with a focus on their representation of eddy properties and dynamics. A comparison with mooring observations shows that both models reasonably simulate hydrography and eddy kinetic energy. Despite differences in model formulation, they show relatively similar eddy properties. The eddies have a mean radius of 4.9 and 5.6 km in ROMS and FESOM, respectively, with slightly more cyclones (ROMS: 54 %, FESOM: 55 %) than anticyclones. The mean lifetime of detected eddies is relatively short in both simulations (ROMS: 10 d, FESOM: 11 d), and the mean travel distance is 35 km in both models. More anticyclones are trapped in deep depressions or move toward deep locations. The two models show comparable spatial patterns of baroclinic and barotropic instability. ROMS has relatively stronger eddy intensity and baroclinic instability, possibly due to its smaller grid size, while FESOM has stronger eddy kinetic energy in the West Spitsbergen Current. Overall, the relatively good agreement between the two models strengthens our confidence in their ability to realistically represent the Fram Strait ocean dynamics and also highlights the need for very high mesh resolution.

Description: Recent modeling studies of ocean circulation in the southern Weddell Sea, Antarctica, project an increase over this century of ocean heat into the cavity beneath Filchner-Ronne Ice Shelf (FRIS). This increase in ocean heat would lead to more basal melting and a modification of the FRIS ice draft. The corresponding change in cavity shape will affect advective pathways and the spatial distribution of tidal currents, which play important roles in basal melting under FRIS. These feedbacks between heat flux, basal melting, and tides will affect the evolution of FRIS under the influence of a changing climate. We explore these feedbacks with a three-dimensional ocean model of the southern Weddell Sea that is forced by thermodynamic exchange beneath the ice shelf and tides along the open boundaries. Our results show regionally-dependent feedbacks that, in some areas, substantially modify the melt rates near the grounding lines of buttressed ice streams that flow into FRIS. These feedbacks are introduced by variations in meltwater production as well as the circulation of this meltwater within the FRIS cavity; they are influenced locally by sensitivity of tidal currents to water column thickness and non-locally by changes in circulation pathways that transport an integrated history of mixing and meltwater entrainment along flow paths. Our results highlight the importance of including explicit tidal forcing in models of future mass loss from FRIS and from the adjacent grounded ice sheet as individual ice stream grounding zones experience different responses to warming of the ocean inflow.

Description: Thinning rates of ice shelves vary widely around Antarctica, and basal melting is a major component of ice shelf mass loss. In this study, we present records of basal melting at a unique spatial and temporal resolution for East Antarctica, derived from autonomous phase-sensitive radars. These records show spatial and temporal variations of basal melting in 2017 and 2018 at Nivlisen, an ice shelf in central Dronning Maud Land. The annually averaged basal melt rates are in general moderate (∼0.8 m yr−1). Radar profiling of the ice shelf shows variable ice thickness from smooth beds to basal crevasses and channels. The highest basal melt rates (3.9 m yr−1) were observed close to a grounded feature near the ice shelf front. Daily time-varying measurements reveal a seasonal melt signal 4 km from the ice shelf front, at an ice draft of 130 m, where the highest daily basal melt rates occurred in summer (up to 5.6 m yr−1). In comparison with wind, air temperatures, and sea ice cover from reanalysis and satellite data, the seasonality in basal melt rates indicates that summer-warmed ocean surface water was pushed by wind beneath the ice shelf front. We observed a different melt regime 35 km into the ice shelf cavity, at an ice draft of 280 m, with considerably lower basal melt rates (annual average of 0.4 m yr−1) and no seasonality. We conclude that warm deep-ocean water at present has a limited effect on the basal melting of Nivlisen. On the other hand, a warming in surface waters, as a result of diminishing sea ice cover, has the potential to increase basal melting near the ice shelf front. Continuous in situ monitoring of Antarctic ice shelves is needed to understand the complex mechanisms involved in ice shelf–ocean interactions.

Description: Coupled ice sheet–ocean models capable of simulating moving grounding lines are just becoming available. Such models have a broad range of potential applications in studying the dynamics of marine ice sheets and tidewater glaciers, from process studies to future projections of ice mass loss and sea level rise. The Marine Ice Sheet–Ocean Model Intercomparison Project (MISOMIP) is a community effort aimed at designing and coordinating a series of model intercomparison projects (MIPs) for model evaluation in idealized setups, model verification based on observations, and future projections for key regions of the West Antarctic Ice Sheet (WAIS). Here we describe computational experiments constituting three interrelated MIPs for marine ice sheet models and regional ocean circulation models incorporating ice shelf cavities. These consist of ice sheet experiments under the Marine Ice Sheet MIP third phase (MISMIP+), ocean experiments under the Ice Shelf-Ocean MIP second phase (ISOMIP+) and coupled ice sheet–ocean experiments under the MISOMIP first phase (MISOMIP1). All three MIPs use a shared domain with idealized bedrock topography and forcing, allowing the coupled simulations (MISOMIP1) to be compared directly to the individual component simulations (MISMIP+ and ISOMIP+). The experiments, which have qualitative similarities to Pine Island Glacier Ice Shelf and the adjacent region of the Amundsen Sea, are designed to explore the effects of changes in ocean conditions, specifically the temperature at depth, on basal melting and ice dynamics. In future work, differences between model results will form the basis for the evaluation of the participating models.

Description: An increasing number of Southern Ocean models now include Antarctic ice shelf cavities, and simulate thermodynamics at the ice-shelf/ocean interface. This adds another level of complexity to Southern Ocean simulations, as ice shelves interact directly with the ocean and indirectly with sea ice. Here we present the first published model intercomparison and evaluation of present-day ocean/sea-ice/ice-shelf interactions, as simulated by two models: a circumpolar Antarctic configuration of MetROMS (ROMS: Regional Ocean Modelling System coupled to CICE: Community Ice CodE) and the global model FESOM (Finite Element Sea-ice/ice-shelf Ocean Model), where the latter is run at two different levels of horizontal resolution. From a circumpolar Antarctic perspective, we compare and evaluate simulated ice shelf basal melting and sub-ice shelf circulation, as well as sea ice properties and Southern Ocean water mass characteristics as they influence the sub-ice shelf processes. Despite their differing numerical methods, the two models produce broadly similar results, and share similar biases in many cases. Both models reproduce many key features of observations, but struggle to reproduce others, such as the high melt rates observed in the small warm-cavity ice shelves of the Amundsen and Bellingshausen Seas. Several differences in model design show a particular influence on the simulations. For example, FESOM's greater topographic smoothing can alter the geometry of some ice shelf cavities enough to affect their melt rates; this improves at higher resolution, since less smoothing is required. In the interior Southern Ocean, the vertical coordinate system affects the degree of water mass erosion due to spurious diapycnal mixing, with MetROMS' terrain-following coordinates leading to more erosion than FESOM's z-coordinates. Finally, increased horizontal resolution in FESOM leads to higher basal melt rates for small ice shelves, through a combination of stronger circulation and small-scale intrusions of warm water from offshore.

Description: Thinning rates of ice shelves vary widely around Antarctica and basal melting is a major component in ice shelf mass loss. In this study, we present records of basal melting, at unique spatial and temporal resolution for East Antarctica, derived from autonomous phase-sensitive radars. These records show spatial and temporal variations of ice shelf basal melting in 2017 and 2018 at Nivlisen, central Dronning Maud Land. The annually averaged melt rates are in general moderate (~ 0.8 m yr-1). Radar profiling of the ice-shelf shows variable ice thickness from smooth beds to basal crevasses and channels. The highest melt rates (3.9 m yr-1) were observed close to a grounded feature near the ice shelf front. Daily time-varying measurements reveal a seasonal melt signal 4 km from the ice shelf front, at an ice draft of 130 m, where the highest daily melt rates occurred in summer (up to 5.6 m yr-1). This seasonality indicates that summer-warmed ocean surface water was pushed by wind beneath the ice shelf front. We observed a different melt regime 35 km into the ice-shelf cavity, at an ice draft of 280 m, with considerably lower melt rates (annual average of 0.4 m yr-1) and no seasonality. We conclude that warm deep ocean water at present has limited effect on the basal melting of Nivlisen. On the other hand, a warming in surface waters, as a result of diminishing sea-ice cover has the potential to increase basal melting near the ice-shelf front. Many ice shelves like Nivlisen are stabilized by pinning points at their ice fronts and these areas may be vulnerable to future change.

Description: Recent modeling studies of ocean circulation in the southern Weddell Sea, Antarctica, project an increase over this century of ocean heat into the cavity beneath Filchner–Ronne Ice Shelf (FRIS). This increase in ocean heat would lead to more basal melting and a modification of the FRIS ice draft. The corresponding change in cavity shape will affect advective pathways and the spatial distribution of tidal currents, which play important roles in basal melting under FRIS. These feedbacks between heat flux, basal melting, and tides will affect the evolution of FRIS under the influence of a changing climate. We explore these feedbacks with a three-dimensional ocean model of the southern Weddell Sea that is forced by thermodynamic exchange beneath the ice shelf and tides along the open boundaries. Our results show regionally dependent feedbacks that, in some areas, substantially modify the melt rates near the grounding lines of buttressed ice streams that flow into FRIS. These feedbacks are introduced by variations in meltwater production as well as the circulation of this meltwater within the FRIS cavity; they are influenced locally by sensitivity of tidal currents to water column thickness (wct) and non-locally by changes in circulation pathways that transport an integrated history of mixing and meltwater entrainment along flow paths. Our results highlight the importance of including explicit tidal forcing in models of future mass loss from FRIS and from the adjacent grounded ice sheet as individual ice-stream grounding zones experience different responses to warming of the ocean inflow.