The Whole Antarctic Ocean Model (WAOM v1.0): Development and evaluation

The Regional Ocean Modeling System (ROMS), including an ice shelf component, has been applied on a circum-Antarctic domain to derive estimates of ice shelf basal melting. Significant improvements made compared to previous models of this scale are the inclusion of tides and a horizontal spatial resolution of 2 km, which is sufficient to resolve on-shelf heat transport by bathymetric troughs and eddy-scale circulation. We run the model with ocean–atmosphere–sea ice conditions from the year 2007 to represent nominal present-day climate. We force the ocean surface with buoyancy fluxes derived from sea ice concentration observations and wind stress from ERA-Interim atmospheric reanalysis. Boundary conditions are derived from the ECCO2 ocean state estimate; tides are incorporated as sea surface height and barotropic currents at the open boundary. We evaluate model results using satellite-derived estimates of ice shelf melting and established compilations of ocean hydrography. The Whole Antarctic Ocean Model (WAOM v1.0) qualitatively captures the broad scale difference between warm and cold regimes as well as many of the known characteristics of regional ice–ocean interaction. We identify a cold bias for some warm-water ice shelves and a lack of high-salinity shelf water (HSSW) formation. We conclude that further calibration and development of our approach are justified. At its current state, the model is ideal for addressing specific, process-oriented questions, e.g. related to tide-driven ice shelf melting at large scales

Mapping Antarctic Suspension Feeder Abundances and Seafloor Food-Availability, and Modeling Their Change After a Major Glacier Calving

Seafloor communities are a critical part of the unique and diverse Antarctic marine life. Processes at the ocean-surface can strongly influence the diversity and abundance of these communities, even when they live at hundreds of meters water depth. However, even though we understand the importance of this link, there are so far no quantitative spatial predictions on how seafloor communities will respond to changing conditions at the ocean surface. Here, we map patterns in abundance of important habitat-forming suspension feeders on the seafloor in East Antarctica, and predict how these patterns change after a major disturbance in the icescape, caused by the calving of the Mertz Glacier Tongue. We use a purpose-built ocean model for the time-period before and after the calving of the Mertz-Glacier Tongue in 2010, data from satellites and a validated food-availability model to estimate changes in horizontal flux of food since the glacier calving. We then predict the post-calving distribution of suspension feeder abundances using the established relationships with the environmental variables, and changes in horizontal flux of food. Our resulting maps indicate strong increases in suspension feeder abundances close to the glacier calving site, fueled by increased food supply, while the remainder of the region maintains similar suspension feeder abundances despite a slight decrease in total food supply. The oceanographic setting of the entire region changes, with a shorter ice-free season, altered seafloor currents and changes in food-availability. Our study provides important insight into the flow-on effects of a changing icescape on seafloor habitat and fauna in polar environments. Understanding these connections is important in the context of current and future effects of climate change, and the mapped predictions of the seafloor fauna as presented for the study region can be used as a decision-tool for planning potential marine protected areas, and for focusing future sampling and monitoring initiatives

Vertical Processes and Resolution Impact Ice Shelf Basal Melting: A Multi-Model Study

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

Vertical processes and resolution impact ice shelf basal melting: A multi-model study

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. © 2020 The Author

Modelling the Response of Ice Shelf Basal Melting to Different Ocean Cavity Environmental Regimes

We present simulation results from a version of the Regional Ocean Modeling System modified for ice shelf/ocean interaction, including the parameterisation of basal melting by molecular diffusion alone. Simulations investigate the differences in melting for an idealised ice shelf experiencing a range of cold to hot ocean cavity conditions. Both the pattern of melt and the location of maximum melt shift due to changes in the buoyancy-driven circulation, in a different way to previous studies. Tidal forcing increases both the circulation strength and melting, with the strongest impact on the cold cavity case. Our results highlight the importance of including a complete melt parameterisation and tidal forcing. In response to the 2.4 degrees C ocean warming initially applied to a cold cavity ice shelf, we find that melting will increase by about an order of magnitude (24 x with tides and 41 x without tides)

Intercomparison of Antarctic ice-shelf, ocean, and sea-ice interactions simulated by MetROMS-iceshelf and FESOM 1.4

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 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 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 coordinate leading to more erosion than FESOM's z coordinate. 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

Experimental design for three interrelated marine ice sheet and ocean model intercomparison projects: MISMIP v. 3 (MISMIP +), ISOMIP v. 2 (ISOMIP +) and MISOMIP v. 1 (MISOMIP1)

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

Future Sea Level Change Under Coupled Model Intercomparison Project Phase 5 and Phase 6 Scenarios From the Greenland and Antarctic Ice Sheets

Projections of the sea level contribution from the Greenland and Antarctic ice sheets (GrIS and AIS) rely on atmospheric and oceanic drivers obtained from climate models. The Earth System Models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6) generally project greater future warming compared with the previous Coupled Model Intercomparison Project phase 5 (CMIP5) effort. Here we use four CMIP6 models and a selection of CMIP5 models to force multiple ice sheet models as part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We find that the projected sea level contribution at 2100 from the ice sheet model ensemble under the CMIP6 scenarios falls within the CMIP5 range for the Antarctic ice sheet but is significantly increased for Greenland. Warmer atmosphere in CMIP6 models results in higher Greenland mass loss due to surface melt. For Antarctica, CMIP6 forcing is similar to CMIP5 and mass gain from increased snowfall counteracts increased loss due to ocean warming

Insights on the vulnerability of Antarctic glaciers from the ISMIP6 ice sheet model ensemble and associated uncertainty

The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from −5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty