Simplified geometry experiments

There have been three stages of ice sheet model intercomparisons based on simplified geometry experiments since the early 1990s [103].

EISMINT I (European Ice Sheet Modeling INiTiative) [27][1] was the first of these and involved only the isothermal shallow ice approximation (SIA). Both fixed margin and moving margin experiments were performed in EISMINT I, and various conclusions were drawn about the several numerical schemes used in the intercomparison. EISMINT I is superceded, however, by verification using the full variety of known exact solutions to the isothermal SIA [26]. The “rediscovery”, since EISMINT I, of the Halfar similarity solution with zero accumulation [104], and verification runs using that solution, already suffices to measure the isothermal SIA performance of PISM more precisely than would be allowed by comparison to EISMINT I results.

EISMINT II [19] pointed out interesting and surprising properties of the thermocoupled SIA. References [23], [105], [106], [40], [107], [22] each interpret the EISMINT II experiments and/or describe attempts to add more complete physical models to “fix” the (perceived and real) shortfalls of ice sheet model behavior on EISMINT II experiments. We believe that the discussion in [108], [40], [23] adequately explains the “spokes” in EISMINT II experiment F as a genuine fluid instability, while [41] and Appendix B of [22] adequately cautions against the continuum model that generates the “spokes” in EISMINT II experiment H. Thus we can move on from that era of controversy. In any case, PISM has built-in support for all of the published and unpublished EISMINT II experiments; these are described in the next subsection.

The ISMIP (Ice Sheet Model Intercomparison Project) [2] round of intercomparisons covers 2008–2013 (at least). There are four components of ISMIP substantially completed, namely HOM = Higher Order Models [109], [110], HEINO = Heinrich Event INtercOmparison [111], [84], MISMIP (below), and MISMIP3d (also below).

PISM participated in HEINO, but this ability is unmaintained. We believe the continuum problem described by HEINO, also used in EISMINT II experiment H (above), is not meaningfully approximate-able because of a required discontinuous jump in the basal velocity field. The continuum problem predicts infinite vertical velocity because of this jump ([22], Appendix B). Details of the numerical schemes and their results are irrelevant if the continuum model makes such a prediction. PISM offers the physical continuum model described in [22], an SIA+SSA hybrid, as an alternative to the continuum model used in ISMIP-HEINO and EISMINT II experiment H. Indeed the SIA+SSA hybrid is offered as a unified shallow model for real ice sheets (section Ice dynamics, the PISM view).

There is no current plan to support ISMIP-HOM [109], [110], but comparison of shallow PISM results to exact Stokes solutions is a goal for PISM evaluation.

A third and fourth ISMIP parts are the two parts of the Marine Ice Sheet Model Intercomparison Project, MISMIP [112] and MISMIP3D [113]. These experiments are supported in PISM, as described in subsections MISMIP and MISMIP3d below.

Intercomparison projects

Footnotes

[1]See http://homepages.vub.ac.be/~phuybrec/eismint.html
[2]See http://homepages.vub.ac.be/~phuybrec/ismip.html

Previous Up Next