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

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 phytoplankto...

Full description

Bibliographic Details
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:
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
Description
Summary: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 ...