Reconciling surface ocean productivity, export fluxes and sediment composition in a global biogeochemical ocean model

This study focuses on an improved representation of the biological soft tissue pump in the global three-dimensional biogeochemical ocean model PISCES. We compare three parameterizations of particle dynamics: (1) the model standard version including two particle size classes, aggregation-disaggregation and prescribed sinking speed; (2) an aggregation-disaggregation model with a particle size spectrum and prognostic sinking speed; (3) a mineral ballast parameterization with no size classes, but prognostic sinking speed. In addition, the model includes a description of surface sediments and organic carbon early diagenesis. The integrated representation of material fluxes from the productive surface ocean down to the sediment-water interface allows taking advantage of surface ocean observations, sediment trap data and exchange fluxes at the sediment-water interface. The capability of the model to reproduce yearly averaged particulate organic carbon fluxes and benthic oxygen demand does at first order not dependent on the resolution of the particle size spectrum. Model results obtained with the standard version and with the one including a particle size spectrum and prognostic sinking speed are not significantly different. Both model versions overestimate particulate organic carbon between 1000 and 2000 m, while deep fluxes are of the correct order of magnitude. Predicted benthic oxygen fluxes correspond with respect to their large scale distribution and magnitude to data based estimates. Modeled particulate organic C fluxes across the mesopelagos are most sensitive to the intensity of zooplankton flux feeding. An increase of the intensity of flux feeding in the standard version results in lower mid- and deep-water particulate organic carbon fluxes, shifting model results to an underestimation of particulate organic carbon fluxes in the deep. The corresponding benthic oxygen fluxes are too low. The model version including the mineral ballast parameterization yields an improved fit between modeled and observed particulate organic carbon fluxes below 2000 m and down to the sediment-water interface. Our results suggest that aggregate formation alone might not be sufficient to drive an intense biological pump. The later is most likely driven by the combined effect of aggregate formation and mineral ballasting.