Upper ocean ecosystems dynamics and iron cycling in a global 3D model

A global 3-D marine ecosystem model with several key phytoplankton functional groups, multiple limiting nutrients, explicit iron cycling, and a mineral ballast/organic matter parameterization is run within a global ocean circulation model. The coupled biogeochemistry/ecosystem/circulation (BEC) model reproduces known basin-scale patterns of primary production, biogenic silica production, calcification, chlorophyll, macronutrient and dissolved iron concentrations. The model captures observed High Nitrate, Low Chlorophyll (HNLC) conditions in the Southern Ocean, subarctic and equatorial Pacific. Spatial distributions of nitrogen fixation are in general agreement with field data, with total N-fixation of 54 Tg N. Diazotrophs directly account for a small fraction of primary production (0.54%) but indirectly support 10% of primary production and nearly 8% of particulate organic carbon (POC) export. Diatoms disproportionately contribute to export of POC out of surface waters, but CaCO3 from the coccolithophores is the key driver of POC flux to the deep ocean in the model. An iron source from shallow ocean sediments is found critical in preventing iron limitation in shelf regions, most notably in the Arctic Ocean, but has a relatively localized impact. In contrast, global-scale primary production, export production, and nitrogen fixation are all sensitive to variations in atmospheric mineral dust inputs. The residence time for dissolved iron in surface waters is estimated to be a few years to a decade. Most of the iron utilized by phytoplankton is from subsurface sources supplied by mixing, entrainment, and ocean circulation. However, due to the short residence time of iron in the upper ocean, this subsurface iron pool is critically dependent on continual replenishment from atmospheric dust deposition and, to a lesser extent, lateral transport from shelf regions.