Quantifying biological carbon pump pathways with a data-constrained mechanistic model ensemble approach

Abstract

The ability to constrain the mechanisms that transport organic carbon into the deep ocean is complicated by the multiple physical, chemical, and ecological processes that intersect to create, transform, and transport particles in the ocean. In this paper we develop and parameterize a data-assimilative model of the multiple pathways of the biological carbon pump (NEMUROBCP). The mechanistic model is designed to represent sinking particle flux, active transport by vertically migrating zooplankton, and passive transport by subduction and vertical mixing, while also explicitly representing multiple biological and chemical properties measured directly in the field (including nutrients, phytoplankton and zooplankton taxa, carbon dioxide and oxygen, nitrogen isotopes, and 234Thorium). Using 30 different data types (including standing stock and rate measurements related to nutrients, phytoplankton, zooplankton, and non-living organic matter) from Lagrangian experiments conducted on 11 cruises from four ocean regions, we conduct an objective statistical parameterization of the model and generate 1 million different potential parameter sets that are used for ensemble model simulations. The model simulates in situ parameters that were assimilated (net primary production and gravitational particle flux) and parameters that were withheld (234Thorium and nitrogen isotopes) with reasonable accuracy. Model results show that gravitational flux of sinking particles and vertical mixing of organic matter from the euphotic zone are more important biological pump pathways than active transport by vertically migrating zooplankton. However, these processes are regionally variable, with sinking particles most important in oligotrophic areas of the Gulf of Mexico and California Current, sinking particles and vertical mixing roughly equivalent in productive coastal upwelling regions and the subtropical front in the Southern Ocean, and active transport an important contributor in the eastern tropical Pacific. We further find that mortality at depth is an important component of active transport when mesozooplankton biomass is high, but it is negligible in regions with low mesozooplankton biomass. Our results also highlight the high degree of uncertainty, particularly amongst mesozooplankton functional groups, that is derived from uncertainty in model parameters. Indeed, variability in BCP pathways between simulations for a specific location using different parameter sets (all with approximately equal misfit relative to observations) is comparable to variability in BCP pathways between regions. We discuss the implications of these results for other data-assimilation approaches and for studies that rely on non-ensemble model outputs.