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.