@article{ReeseWinkelmannGudmundsson2018,
  author    = {Reese, Ronja and Winkelmann, Ricarda and Gudmundsson, Gudmundur Hilmar},
  title     = {Grounding-line flux formula applied as a flux condition in numerical simulations fails for buttressed Antarctic ice streams},
  series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union},
  volume    = {12},
  journal   = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union},
  number    = {10},
  publisher = {Copernicus},
  address   = {G{\"o}ttingen},
  issn      = {1994-0416},
  doi       = {10.5194/tc-12-3229-2018},
  pages     = {3229 -- 3242},
  year      = {2018},
  abstract  = {Currently, several large-scale ice-flow models impose a condition on ice flux across grounding lines using an analytically motivated parameterisation of grounding-line flux. It has been suggested that employing this analytical expression alleviates the need for highly resolved computational domains around grounding lines of marine ice sheets. While the analytical flux formula is expected to be accurate in an unbuttressed flow-line setting, its validity has hitherto not been assessed for complex and realistic geometries such as those of the Antarctic Ice Sheet. Here the accuracy of this analytical flux formula is tested against an optimised ice flow model that uses a highly resolved computational mesh around the Antarctic grounding lines. We find that when applied to the Antarctic Ice Sheet the analytical expression provides inaccurate estimates of ice fluxes for almost all grounding lines. Furthermore, in many instances direct application of the analytical formula gives rise to unphysical complex-valued ice fluxes. We conclude that grounding lines of the Antarctic Ice Sheet are, in general, too highly buttressed for the analytical parameterisation to be of practical value for the calculation of grounding-line fluxes.},
  language  = {en}
}
@article{ReeseAlbrechtMengeletal.2018,
  author    = {Reese, Ronja and Albrecht, Torsten and Mengel, Matthias and Asay-Davis, Xylar and Winkelmann, Ricarda},
  title     = {Antarctic sub-shelf melt rates via PICO},
  series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union},
  volume    = {12},
  journal   = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union},
  number    = {6},
  publisher = {Copernicus},
  address   = {G{\"o}ttingen},
  issn      = {1994-0416},
  doi       = {10.5194/tc-12-1969-2018},
  pages     = {1969 -- 1985},
  year      = {2018},
  abstract  = {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 sub-shelf 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 ma(-1) for cold sub-shelf cavities, for example, underneath Ross or Ronne ice shelves, to 16 ma(-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.},
  language  = {en}
}