Modeling

MISOMIP1

The first stage of MISOMIP activities (MISOMIP1, 2014-…) consisted of an idealized model intercomparison exercise with 3 components:

  1. MISMIP+: 3rd Marine Ice Sheet Model Intercomparison Project

Co-chairs: Stephen Cornford and Hilmar Gudmundsson

The following description has been adapted from Asay-Davis et al. (2016). It should be noted that the GMD manuscript remains the definitive description of the experiments and the abbreviated description that follows is meant only to present some motivation and a brief overview of the experiments.

The previous Marine Ice Sheet Model Intercomparison Projects, MISMIP and MISMIP3d, tested the capabilities of ice sheet models to simulate advance and retreat cycles under changes in ice softness and basal sliding, respectively. MISMIP tested flowline models in one horizontal dimension (1HD), while MISMIP3d required modeling in two horizontal dimensions (2HD). Each experiment taught the community a great deal about the numerical behavior of ice-sheet models of various types as well as the similarities and differences in the results they produced. In particular, the MISMIP results showed that steady-state grounding-line positions could differ markedly depending on the resolution, type of stress approximation, and discretization methods employed. Follow-up studies found that models with fixed grids (as opposed to those that track the grounding line in time) and without sub-grid-scale parameterizations of the grounding line require grounding-line resolution on the order of hundreds of meters to accurately reproduce grounding-line dynamics.

It was clear in discussions of a follow-up intercomparison exercise that the MISMIP3d experimental design had three shortcomings as a test of 2HD marine ice sheet models. First, it started from a steady state that was invariant in the cross-flow direction—that is, 1HD—and did not involve significant lateral stresses.  Second, the initial grounding lines of the shallow-shelf approximation (SSA) models were around 80km downstream from the Stokes models, but the grounding line only moved about 20 km in the perturbation experiment. That left an obvious question entirely unanswered: in a realistic simulation with the model parameters chosen to match geometry and velocity derived from observations, and thus with prescribed initial conditions, does the SSA provide a good approximation to the Stokes model? Third, grounding line migration was driven by changes to the basal traction field, rather than the ice shelf melting that is thought to be the dominant driver of present-day grounding-line retreat in West Antarctica.

MISMIP+ has been designed to address each of the shortcomings above. Regarding the first, the chosen geometry, based on Gudmundsson et al. (2012), results in strong lateral stresses that buttress the ice stream, and, given particular parameter choices, results in a stable grounding line crossing a retrograde slope, a configuration not possible in 1HD.  Regarding the second, modelers are free to choose certain model parameters so that their initial grounding line is close to that of a reference model, and in preliminary tests two models that bracketed the high resolution MISMIP3d results have been found to have grounding lines within a few kilometers of one another in steady state. Finally, extensive grounding line retreat is driven by sub-shelf melt rates.

The MISMIP+ experiments are initialized by running the model with no melting to steady state. The first MISMIP+ experiment (Ice0) continues the run with no melting for one hundred years to demonstrate that the initial state remains steady over this time period.  The second experiment (Ice1) explores the response of the ice sheet to a strong melt perturbation by prescribes a depth-dependent parameterization of basal melting with a peak melt rate of ~75 m/a, driving 100 (optionally 1000) years of grounding-line retreat.  After 100 years, the melt rate is set back to zero, allowing the grounding line to re-advance over 100 (optionally 900) years.  The third experiment (Ice2) emulates a large calving event by prescribing a strong melt rate (100 m/a) beyond a prescribed calving front for 100 (optionally 1000) years, followed by a removal of the perturbation for an additional 100 (optionally 900) years.

Evolution of the basal traction and ice shelf melt rate fields from a BISICLES MISMIP+ run.

2) ISOMIP+: 2nd Ice Shelf Ocean Model Intercomparison Project

Co-Chairs: Xylar Asay-Davis, Ben Galton-Fenzi

The following description has been adapted from Asay-Davis et al. (2015, under review). It should be noted that the GMD manuscript remains the definitive description of the experiments and the abbreviated description that follows is meant only to present some motivation and a brief overview of the experiments. The following description has been adapte

The ISOMIP+ experiments have been designed to make a number of improvements on the original ISOMIP experiments.  Whereas ISOMIP used highly idealized geometry (the ocean column at the grounding line was 200 m thick, the ice draft sloped linearly with latitude and was invariant with longitude, and the bedrock was perfectly flat), ISOMIP+ makes use of relatively complex geometry from MISMIP+ BISICLES simulations, including an ocean cavity that reaches zero thickness at the grounding line.  Where ISOMIP uses a velocity-independent, two-equation formulationof the melt boundary conditions, ISOMIP+ uses the velocity-dependent three-equation formulation more commonly used in realistic model configurations. ISOMIP specified ~10 km resolution, too coarse to resolve the ~9 km Rossby radius of deformation, and large values of the horizontal viscosity and diffusivities, leading to a laminar flow that evolved toward steady state without eddies or other fluctuations.   In contrast, ISOMIP+ runs will typically use smaller horizontal viscosity and diffusivities and higher resolution (~2 km), allowing for mesoscale eddies and unsteady flow. A smaller computational domain makes the experiments computationally feasible despite the higher resolution.  ISOMIP+ should provide more appropriate test cases than the original ISOMIP for realistic experiments, particularly for those focused on the Amundsen Sea region of WAIS.

ISOMIP+ prescribes five idealized experiments (Ocean0 through Ocean4) for ocean models that support ice-shelf cavities. Ocean0-2 have fixed topography while Ocean3 and Ocean4 have prescribed, evolving ice topography.  The experiments are designed to improve community understanding of the models themselves (both their capabilities and their limitation) as well as physical processes involved ice-ocean interactions.

ISOMIP+ has three main goals:
1. Provide comparable model results as a sanity check during model development,
2. Explore the ocean-model setup on its own before coupling to an ice-sheet model in the MISOMP1 experiments (see below),
3. Provide a basic setup from which a large variety of parameter and process studies could usefully be performed.

The experiments are:

  • Ocean0: 1-2 year run with warm forcing and static ice-shelf geometry,
  • Ocean1: 20-year run with cold-to-warm forcing and static ice-shelf geometry,
  • Ocean2: 20-year run with warm-to-cold forcing and static ice-shelf geometry,
  • Ocean3: 100-year run with warm forcing and retreating ice-shelf geometry,
  • Ocean4: 100-year run with cold forcing and advancing ice-shelf geometry.


A  transect of temperature  through  the  center of the  ISOMIP+  domain  from  example  simulations  from  the  POP2x ocean model, showing the evolution of the Ocean1 (top) and Ocean2 (bottom) simulations

Topography data sets for each experiment comes from a MISMIP+ simulation (see above) using BISICLES with SSA.

Example output from POP2x runs in the “common” (COM) configuration includes movies as well as NetCDF files demonstrating the format of the output grid. Participants may find it useful to use these files for their own output, replacing the POP2x results with their own.

3) MISOMIP1: 1st Marine Ice Sheet-Ocean Model Intercomparison Project

Co-Chairs: Xylar Asay-Davis, Helene Seroussi

The following description has been adapted from Asay-Davis et al. (2016).  It should be noted that the GMD manuscript remains the definitive description of the experiments and the abbreviated description that follows is meant only to present some motivation and a brief overview of the experiments.

The first Marine Ice Sheet-Ocean Model Intercomparison Project (MISOMIP1) combines MISMIP+ and ISOMIP+ (see above).  In some ways, the MISOMIP1 setup is similar to that of Goldberg et al. (2012a,b) in that it includes a narrow channel with strong ice-shelf buttressing and strong far-field restoring in the ocean.  MISOMIP1 differs from this previous work in having 1) steeper channel walls, meaning a stronger change in buttressing as the ice-shelf thickness changes, 2) a larger region of open ocean allowing for ocean dynamics both inside and outside the cavity, and 3) making use of  a bedrock topography with an upward-sloping region in the ice-flow direction, allowing us to investigate the possibility that thinning or other changes in the state of the ice sheet could trigger a marine ice-sheet instability (MISI).

MISOMIP1 prescribes two coupled ice sheet-ocean experiments (IceOcean1 and IceOcean2), each with two parts. We expect the MISOMIP1 experiment to play an analogous role in evaluating coupled ice sheet-ocean systems to that of the ISOMIP projects for standalone ocean models with ice-shelf cavities and the MISMIP projects for ice-sheet models.

IceOcean1 begins with the ice-sheet steady state that also served as the initial conditions for the Ice0, Ice1 and Ice2 experiments. Unlike in ISOMIP+, IceOcean1 does not include a dynamic calving criterion. Ice is allowed to become as thin as the ice sheet and ocean components permit (potentially zero thickness) without calving.

For the first one hundred years, warm restoring (maximum temperature of 1 C) is applied near the ocean’s northern boundary, which is expected to induce strong melting and subsequent rapid ice retreat.  After one hundred years, the restoring abruptly switches to cold conditions (constant temperature of -1.9 C) at the ocean’s northern boundary, which should shut reduce melting within a decade and allow ice to re-advance for the remaining 100 years of simulation.

Specifying calving is a major problem in the design of MISOMIP1. There was general agreement in the community that ice-sheet models have not been shown to behave reliably with dynamic calving, nor is there any consensus about which calving parameterizations are appropriate or physically realistic.  Nevertheless, we felt that it was important for testing the robustness of the ice-sheet and ocean components in MISOMIP1 that there be an experiment with a dynamic, sheer cliff at the calving front. We include an optional coupled experiment, IceOcean2, that is identical to IceOcean1 except that it includes dynamic calving in the ice-sheet component.


Example results from a POPSICLES simulation of the IceOcean1 (top) and IceOcean2 (bottom) experiment showing transects through the center of the domain of ocean temperature and ice topography.

Example output from POPSICLES runs includes movies and example NetCDF results.  Participants may find it useful to use these files for their own output, replacing the POPSICLES results with their own.

  • IceOcean1_COM_ocean_POPSICLES.nc (1.8G)
  • IceOcean1_COM_ice_POPSICLES.nc (coming soon)
  • IceOcean2_COM_ocean_POPSICLES.nc (1.8G)
  • IceOcean2_COM_ice_POPSICLES.nc (coming soon)
For further details about these experiments, we refer to the accompanying experimental design:
Asay-Davis, X. S., Cornford, S. L., Durand, G., Galton-Fenzi, B. K., Gladstone, R. M., Gudmundsson, G. H., Hattermann, T., Holland, D. M., Holland, D., Holland, P. R., Martin, D. F., Mathiot, P., Pattyn, F., and Seroussi, H. (2016). 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), Geoscientific Model Development, 9, 2471–2497, https://doi.org/10.5194/gmd-9-2471-2016

MISOMIP2

The second stage of MISOMIP activities (MISOMIP2, 2019-…) proposes standalone ocean and ice experiments, as well as coupled ice-ocean experiments. Previous acronyms MISMIP+ and ISOMIP+ will be discontinued, and all experiments are collectively referred to as MISOMIP2. Experiments will use realistic geometries and forcing, with particular focus on the Amundsen and Weddell seas. We invite contributions from global and regional ocean model configurations, if ice shelves are included. We invite contributions from circum-Antarctic or regional ice model setups, with a basal melt parameterization of your choice. We adopt a ‘come as you are’ approach, with an option for common experiments. The project will assemble and provide a ‘MIPkit’ alongside the publication of the experimental protocol. The MIPkit will contain pre-gridded data of ocean and ice observations for model-data comparison, as well as Python and Matlab packages for basic pre-processing and plotting of MISOMIP2 output. Further details will be made available at a later stage. A publication that outlines the experimental design is currently in preparation, with the table below summarising the basic framework.