Evaluation of the Transport Matrix Method for simulation of ocean biogeochemical tracers.

Kvale, Karin Frances , Khatiwala, Samar, Dietze, Heiner , Kriest, Iris and Oschlies, Andreas (2017) Evaluation of the Transport Matrix Method for simulation of ocean biogeochemical tracers. Open Access Geoscientific Model Development, 10 . pp. 2425-2445. DOI 10.5194/gmd-10-2425-2017.

[thumbnail of gmd-10-2425-2017.pdf]
gmd-10-2425-2017.pdf - Published Version
Available under License Creative Commons: Attribution 3.0.

Download (3MB) | Preview

Supplementary data:


Conventional integration of Earth system and ocean models can accrue considerable computational expenses, particularly for marine biogeochemical applications. "Offline" numerical schemes in which only the biogeochemical tracers are time stepped and transported using a pre-computed circulation field can substantially reduce the burden and are thus an attractive alternative. One such scheme is the "transport matrix method" (TMM), which represents tracer transport as a sequence of sparse matrix–vector products that can be performed efficiently on distributed-memory computers. While the TMM has been used for a variety of geochemical and biogeochemical studies, to date the resulting solutions have not been comprehensively assessed against their "online" counterparts. Here, we present a detailed comparison of the two. It is based on simulations of the state-of-the-art biogeochemical sub-model embedded within the widely used coarse-resolution University of Victoria Earth System Climate Model (UVic ESCM). The default, non-linear advection scheme was first replaced with a linear, third-order upwind-biased advection scheme to satisfy the linearity requirement of the TMM. Transport matrices were extracted from an equilibrium run of the physical model and subsequently used to integrate the biogeochemical model offline to equilibrium. The identical biogeochemical model was also run online. Our simulations show that offline integration introduces some bias to biogeochemical quantities through the omission of the polar filtering used in UVic ESCM and in the offline application of time-dependent forcing fields, with high latitudes showing the largest differences with respect to the online model. Differences in other regions and in the seasonality of nutrients and phytoplankton distributions are found to be relatively minor, giving confidence that the TMM is a reliable tool for offline integration of complex biogeochemical models. Moreover, while UVic ESCM is a serial code, the TMM can be run on a parallel machine with no change to the underlying biogeochemical code, thus providing orders of magnitude speed-up over the online model.

Document Type: Article
Research affiliation: OceanRep > GEOMAR > FB2 Marine Biogeochemistry > FB2-BM Biogeochemical Modeling
Kiel University > Kiel Marine Science
OceanRep > The Future Ocean - Cluster of Excellence
OceanRep > SFB 754
OceanRep > SFB 754 > B1
Refereed: Yes
Open Access Journal?: Yes
Publisher: Copernicus Publications (EGU)
Related URLs:
Projects: SFB754, Future Ocean, Opendap
Date Deposited: 08 Feb 2017 08:54
Last Modified: 06 Feb 2020 09:10
URI: https://oceanrep.geomar.de/id/eprint/36441

Actions (login required)

View Item View Item