Coupling framework (1.0) for the ice sheet model PISM (1.1.1) and the ocean model MOM5 (5.1.0) via the ice-shelf cavity module PICO

The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high resolutio...

Full description

Bibliographic Details
Main Authors: Kreuzer, Moritz, Reese, Ronja, Huiskamp, Willem Nicholas, Petri, Stefan, Albrecht, Torsten, Feulner, Georg, Winkelmann, Ricarda
Format: Text
Language:English
Published: 2020
Subjects:
Online Access:https://doi.org/10.5194/gmd-2020-230
https://gmd.copernicus.org/preprints/gmd-2020-230/
Description
Summary:The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high resolution configurations, limiting these studies to individual glaciers or regions over short time scales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM with the global ocean general circulation model MOM5 via the ice-shelf cavity module PICO. Since ice-shelf cavities are not resolved by MOM5, but parameterized with the box model PICO, the framework allows the ice sheet and ocean model to be run at resolution of 16 km and 3 degree, respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the entire Antarctic Ice Sheet and the Earth system over time spans on the order of centuries to millennia. In this study we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet model is calculated by PICO from modeled ocean temperatures and salinities at the depth of the continental shelf and, vice versa, the resulting mass and energy fluxes from the melting at the ice-ocean interface are transferred to the ocean model. Mass and energy fluxes are shown to be conserved to machine precision across the considered model domains. The implementation is computationally efficient as it introduces only minimal overhead. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions and time scales between the ice and ocean model in a generic way, and can thus be adopted to a wide range of model setups.