ALBERT

Description: The most rapid loss of ice from the Antarctic Ice Sheet is observed where ice streams flow into the ocean and begin to float, forming the great Antarctic ice shelves that surround much of the continent. Because these ice shelves are floating, their thinning does not greatly influence sea level. However, they also buttress the ice streams draining the ice sheet, and so ice shelf changes do significantly influence sea level by altering the discharge of grounded ice. Currently, the most significant loss of mass from the ice shelves is from melting at the base (although iceberg calving is a close second). Accessing the ocean beneath ice shelves is extremely difficult, so numerical models are invaluable for understanding the processes governing basal melting. This paper describes the different ways in which ice shelf/ocean interactions are modeled and discusses emerging directions that will enhance understanding of how the ice shelves are melting now and how this might change in the future.

Description: Ice sheets and ice shelves are linked by the transition zone, the region where the grounded ice lifts off the bedrock and begins to float. Adequate resolution of the transition zone is necessary for numerically accurate ice sheet–ice shelf simulations. The required resolution depends on how the basal physics is parameterized. We propose a new, simple parameterization of the basal hydrology in a one-dimensional vertically integrated model. This parameterization accounts for connectivity between the basal hydrological system and the ocean in the transition zone. Our model produces a smooth transition between finite basal friction in the ice sheet and zero basal friction in the ice shelf. Through a set of experiments based on the Marine Ice Sheet Model Intercomparison Project (MISMIP), we show that a smoother basal shear stress, in addition to adding physical realism, significantly improves the numerical accuracy of our fixed-grid model, allowing for reliable grounding-line dynamics at resolutions ~1 km.

Description: Ice sheets and ice shelves are linked by the transition zone, the region where flow dominated by vertical shear stress makes a transition to flow dominated by extensional stress. Adequate resolution of the transition zone is necessary for numerically accurate ice sheet–ice shelf simulations. The required resolution depends on how the basal physics is parameterized. We propose a~new, simple parameterization of the effective pressure near the grounding line, combined with an existing friction law linking effective pressure to basal stress and sliding, in a one-dimensional, fixed-grid, vertically integrated model. This parameterization represents connectivity between the basal hydrological system and the ocean in the transition zone. Our model produces a smooth transition between finite basal friction in the ice sheet and zero basal friction in the ice shelf. In a set of experiments based on the Marine Ice Sheet Model Intercomparison Project (MISMIP), we show that with a smoother basal shear stress, the model yields accurate steady-state results at a fixed-grid resolution of ~1 km.

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 in 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 evaluation of the participating models.

Description: Jupiter's Great Red Spot (GRS) is one of its most distinct and enduring features, having been continuously observed since the 1800's. It currently spans the smallest latitude and longitude size ever recorded. Here we show analyses of 2014 Hubble spectral imaging data to study the color, structure and internal dynamics of this long-live storm.

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: 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: Ocean-induced melting below ice shelves is one of the dominant drivers for mass loss from the Antarctic Ice Sheet at present. An appropriate representation of subshelf melt rates is therefore essential for model simulations of marine-based ice sheet evolution. Continental-scale ice sheet models often rely on simple melt-parameterizations, in particular for long-term simulations, when fully coupled ice–ocean interaction becomes computationally too expensive. Such parameterizations can account for the influence of the local depth of the ice-shelf draft or its slope on melting. However, they do not capture the effect of ocean circulation underneath the ice shelf. Here we present the Potsdam Ice-shelf Cavity mOdel (PICO), which simulates the vertical overturning circulation in ice-shelf cavities and thus enables the computation of sub-shelf melt rates consistent with this circulation. PICO is based on an ocean box model that coarsely resolves ice shelf cavities and uses a boundary layer melt formulation. We implement it as a module of the Parallel Ice Sheet Model (PISM) and evaluate its performance under present-day conditions of the Southern Ocean. We identify a set of parameters that yield two-dimensional melt rate fields that qualitatively reproduce the typical pattern of comparably high melting near the grounding line and lower melting or refreezing towards the calving front. PICO captures the wide range of melt rates observed for Antarctic ice shelves, with an average of about 0.1 m a−1 for cold sub-shelf cavities, for example, underneath Ross or Ronne ice shelves, to 16 m a−1 for warm cavities such as in the Amundsen Sea region. This makes PICO a computationally feasible and more physical alternative to melt parameterizations purely based on ice draft geometry.

Description: Understanding ice shelf–ocean interaction is fundamental to projecting the Antarctic ice sheet response to a warming climate. Numerical ice shelf–ocean models are a powerful tool for simulating this interaction, yet are limited by inherent model weaknesses and scarce observations, leading to parameterisations that are unverified and unvalidated below ice shelves. We explore how different models simulate ice shelf–ocean interaction using the 2nd Ice Shelf–Ocean Model Intercomparison Project (ISOMIP+) framework. Vertical discretisation and resolution of the ocean model are shown to have a significant effect on ice shelf basal melt rate, through differences in the distribution of meltwater fluxes and the calculation of thermal driving. Z-coordinate models, which generally have coarser vertical resolution in ice shelf cavities, may simulate higher melt rates compared to terrain-following coordinate models. This is due to the typically higher resolution of the ice–ocean boundary layer region in terrain following models, which allows better representation of a thin meltwater layer, increased stratification, and as a result, better insulation of the ice from water below. We show that a terrain-following model, a z-level coordinate model and a hybrid approach give similar results when the effective vertical resolution adjacent to the ice shelf base is similar, despite each model employing different paradigms for distributing meltwater fluxes and sampling tracers for melting. We provide a benchmark for thermodynamic ice shelf–ocean interaction with different model vertical coordinates and vertical resolutions, and suggest a framework for any future ice shelf–ocean thermodynamic parameterisations.