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 manuscript 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 one 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 surface ocean 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, sinking particles and vertical mixing roughly equivalent in productive regions of the CCE 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 biomasses are high, but that 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, with important implications from results that rely on non-ensemble model outputs. We also discuss the implications of our results for other data assimilation approaches.