Coupling framework (1.0) for the PISM (1.1.4) ice sheet model and the MOM5 (5.1.0) ocean model via the PICO ice shelf cavity model in an Antarctic domain

Kreuzer, Moritz, Reese, Ronja, Huiskamp, Willem Nicholas, Petri, Stefan, Albrecht, Torsten, Feulner, Georg and Winkelmann, Ricarda (2021) Coupling framework (1.0) for the PISM (1.1.4) ice sheet model and the MOM5 (5.1.0) ocean model via the PICO ice shelf cavity model in an Antarctic domain. Geoscientific Model Development, 14 (6). pp. 3697-3714. ISSN 1991-9603

gmd-14-3697-2021.pdf - Published Version
Available under License Creative Commons Attribution 4.0.

Download (5MB) | Preview
Official URL:


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 timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Potsdam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3∘ respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of 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 component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice–ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups.

Item Type: Article
Additional Information: Funding information: This research has been supported by the Deutsche Forschungsgemeinschaft (DFG; grant nos. WI4556/4-1, WI4556/3-1, and WI4556/2-1), the Horizon 2020 programme (grant no. TiPACCs 820575), the German Federal Ministry of Education and Research (BMBF, FONA; grant nos. FKZ:01LP1925D, FKZ:01LP1504D, and FKZ:01LP1502C), NASA (grant no. NNX17AG65G), NSF (grant nos. PLR-1603799 and PLR-1644277), the European Regional Development Fund (ERDF), the German Federal Ministry of Education and Research (BMBF), and the Land Brandenburg. The publication of this article was funded by the Open Access Fund of the Leibniz Association.
Subjects: F700 Ocean Sciences
Department: Faculties > Engineering and Environment > Geography and Environmental Sciences
Depositing User: John Coen
Date Deposited: 18 Nov 2021 11:41
Last Modified: 18 Nov 2021 11:45

Actions (login required)

View Item View Item


Downloads per month over past year

View more statistics