Abstract
Abstract. 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 of their pre-industrial Southern Ocean habitat. Global opal export production thus drops to 75 % of its pre-industrial value by 2300. Model nutrients such as phosphate, silicate and nitrate build up along the Southern Ocean particle export pathway, but dissolved iron (for which ocean sources are held constant) increases in the upper ocean. This different behaviour of iron is attributed to a reduction of low-latitude NPP (and consequently, a reduction in both uptake and export and particle, including calcite scavenging), an increase in seawater temperatures (raising the solubility of particulate iron) and stratification that “traps” the iron near the surface. These results are meant to serve as a baseline for sensitivity assessments to be undertaken with this model in the future.
Highlights
It has become apparent in recent decades that the representation of elemental cycles of nutritive elements (N, P, Si, Fe) is important in order to simulate critical biogeochemical feedbacks (Heinze et al, 2019)
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
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
Summary
It has become apparent in recent decades that the representation of elemental cycles of nutritive elements (N, P, Si, Fe) is important in order to simulate critical biogeochemical feedbacks (Heinze et al, 2019). As far as we are aware, this model is unique among ecosystem models embedded within Earth system models in that it combines established parameterizations of a diatom functional type (i.e. fast growth rates and inefficient nutrient utilization, Tréguer et al, 2018) with a novel temperature-dependent opal dissolution parameterization and bottom-up and top-down competition for resources with explicit, fully independent calcifying phytoplankton and nitrogen fixing diazotrophs. These assumptions bias KMBM3 diatoms towards large species in the well-mixed high latitudes and early seasonal succession. A ballast model applied to the calcifiers (Kvale et al, 2015b) further differentiates KMBM3, as does the ability to operate in an “offline” parameter optimization framework (Kvale et al, 2017; Yao et al, 2019)
Sign-in/Register to view highlights & summary Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
