Explicit silicate cycling in the Kiel Marine Biogeochemistry Model version 3 (KMBM3) embedded in the UVic ESCM version 2.9

• Published in: Geoscientific Model Development
• Main Authors: Kvale, Karin, Keller, David P., Koeve, Wolfgang, Meissner, Katrin J., Somes, Christopher J., Yao, Wanxuan, Oschlies, Andreas
• Format: Article in Journal/Newspaper
• Language: English
• Published: Copernicus Publications 2021
• Subjects: article; Verlagsveröffentlichung; Antarc*; Antarctic; Southern Ocean
• Online Access: https://doi.org/10.5194/gmd-14-7255-2021; https://noa.gwlb.de/receive/cop_mods_00058953; https://noa.gwlb.de/servlets/MCRFileNodeServlet/cop_derivate_00058563/gmd-14-7255-2021.pdf; https://gmd.copernicus.org/articles/14/7255/2021/gmd-14-7255-2021.pdf

We describe and test a new model of biological marine silicate cycling, implemented in the Kiel Marine Biogeochemical Model version 3 (KMBM3), embedded in the University of Victoria Earth System Climate Model (UVic ESCM) version 2.9. This new model adds diatoms, which are a key component of the biological carbon pump, to an existing ecosystem model. This new model combines previously published parameterizations of a diatom functional type, opal production and export with a novel, temperature-dependent dissolution scheme. Modelled steady-state biogeochemical rates, carbon and nutrient distributions are similar to those found in previous model versions. The new model performs well against independent ocean biogeochemical indicators and captures the large-scale features of the marine silica cycle to a degree comparable to similar Earth system models. Furthermore, it is computationally efficient, allowing both fully coupled, long-timescale transient simulations and “offline” transport matrix spinups. We assess the fully coupled model against modern ocean observations, the historical record starting from 1960 and a business-as-usual atmospheric CO2 forcing to the year 2300. The model simulates a global decline in net primary production (NPP) of 1.4 % having occurred since the 1960s, with the strongest declines in the tropics, northern midlatitudes and Southern Ocean. The simulated global decline in NPP reverses after the year 2100 (forced by the extended RCP8.5 CO2 concentration scenario), and NPP returns to 98 % of the pre-industrial rate by 2300. This recovery is dominated by increasing primary production in the Southern Ocean, mostly by calcifying phytoplankton. Large increases in calcifying phytoplankton in the Southern Ocean offset a decline in the low latitudes, producing a global net calcite export in 2300 that varies only slightly from pre-industrial rates. Diatom distribution moves southward in our simulations, following the receding Antarctic ice front, but diatoms are outcompeted by calcifiers across most ...