Representation of shallow grounding zones in an ice shelf-ocean model with terrain-following coordinates

Numerical modeling shows great potential as a method for investigating and predicting the future development of ice shelves in a warming climate. The quality of ice shelf-ocean models is continuously improving but some limitations remain. For models using a terrain-following vertical coordinate, one...

Full description

Bibliographic Details
Published in:Ocean Modelling
Main Authors: Schnaase, Frank, Timmermann, Ralph
Format: Article in Journal/Newspaper
Language:unknown
Published: Elsevier 2019
Subjects:
Online Access:https://oceanrep.geomar.de/id/eprint/49221/
https://doi.org/10.1016/j.ocemod.2019.101487
Description
Summary:Numerical modeling shows great potential as a method for investigating and predicting the future development of ice shelves in a warming climate. The quality of ice shelf-ocean models is continuously improving but some limitations remain. For models using a terrain-following vertical coordinate, one such limitation is the enforcement of a minimum water-column thickness beneath ice shelves by adjustment of bottom topography. How this local distortion of bathymetry from reality affects modeled melt rates and cavity circulation is unknown so far. To quantify this effect, simulations with the Finite Element Sea ice–ice shelf–Ocean Model (FESOM) were executed on four different grids with minimum water-column thicknesses of 20 m, 50 m, 100 m and 200 m. While we use a global model grid, modifications of bathymetry are applied only to Filchner–Ronne Ice Shelf. We show that the choice of minimum water-column thickness does not affect the total basal melt rate of this cold-water ice shelf but does in fact impact the distribution of melt rates with significant differences between experiments in the magnitude of melting near the grounding lines.