Modelling the ecosystem response to iron fertilization in the subarctic NE Pacific: The influence of grazing, and Si and N cycling on CO2 drawdown

• Published in: Deep Sea Research Part II: Topical Studies in Oceanography
• Main Authors: Denman, K. L., Völker, Christoph, Pena, M. A., Rivkin, R. B.
• Format: Article in Journal/Newspaper
• Language: unknown
• Published: 2006
• Subjects: Subarctic
• Online Access: https://epic.awi.de/id/eprint/14950/; https://doi.org/10.1016/j.dsr2.2006.05.026; https://hdl.handle.net/10013/epic.25155

We have employed a coupled one-dimensional mixed layer /ecosystem /carbon cycle model to simulate both the normal annual cycle and the iron-fertilization experiment in the subarctic NE Pacific Ocean near Ocean Station P (50°N, 145°W) during summer 2002. We considered two size classes of phytoplankton, the larger representing diatoms, where each size class has a different degree of iron limitation, and compartments for nitrate, ammonium, microzooplankton, sinking detritus, and a prescribed annual cycle in mesozooplankton. The base ecosystem model is formulated in terms of nitrogen, but is coupled to sub-models of silicon and carbon. Diatoms formed aggregates during blooms that sink rapidly from the surface ocean, and diatoms also can be grazed by microzooplankton, consistent with observations. Using the same parameter set as for the base ecosystem model, we reproduce the basic responses to fertilization: an initial bloom of small phytoplankton (including calcifying coccolithophorids), followed rapidly by an increase of microzooplankton biomass; a continuing increase in diatoms that peak as silicate becomes limiting; and a later rapid sinking event of both carbon and silica particulates. Generally this sequence proceeds more rapidly in simulations than in situ. Simulations of the fertilization response show little sensitivity to the assumed fraction of small phytoplankton that are calcifiers, but a strong sensitivity to the assumed diatom uptake ratio Si:N. With an uptake ratio of 2.5, silicate is rapidly exhausted after fertilization, and 8 months later the regional pCO2 was 14 μatm higher than in the case with no fertilization (assuming no exchange with surrounding waters): for all other simulations the pCO2 anomaly is negative (indicating increased CO2 exchange from the atmosphere) but small, 25 μatm, suggesting a persistence for a single large-scale fertilization of less than 1 year. Simulated mixed layer Si:N drawdown ratios for different fixed diatom uptake ratios of Si:N illustrate the dangers of ...