Modeling Responses of Diatom Productivity and Biogenic Silica Export to Iron Enrichment in the Equatorial Pacific Ocean

Using a three-dimensional physical-biogeochemical model, we have investigated the modeled responses of diatom productivity and biogenic silica export to iron enrichment in the equatorial Pacific, and compared the model simulation with in situ (IronEx II) iron fertilization results. In the eastern equatorial Pacific, an area of 540,000 km(2) was enhanced with iron by changing the photosynthetic efficiency and silicate and nitrogen uptake kinetics of phytoplankton in the model for a period of 20 days. The vertically integrated Chl a and primary production increased by about threefold 5 days after the start of the experiment, similar to that observed in the IronEx II experiment. Diatoms contribute to the initial increase of the total phytoplankton biomass, but decrease sharply after 10 days because of mesozooplankton grazing. The modeled surface nutrients (silicate and nitrate) and TCO(2) anomaly fields, obtained from the difference between the "iron addition'' and "ambient'' (without iron) concentrations, also agreed well with the IronEx II observations. The enriched patch is tracked with an inert tracer similar to the SF6 used in the IronEx II. The modeled depth-time distribution of sinking biogenic silica (BSi) indicates that it would take more than 30 days after iron injection to detect any significant BSi export out of the euphotic zone. Sensitivity studies were performed to establish the importance of fertilized patch size, duration of fertilization, and the role of mesozooplankton grazing. A larger size of the iron patch tends to produce a broader extent and longer-lasting phytoplankton blooms. Longer duration prolongs phytoplankton growth, but higher zooplankton grazing pressure prevents significant phytoplankton biomass accumulation. With the same treatment of iron fertilization in the model, lowering mesozooplankton grazing rate generates much stronger diatom bloom, but it is terminated by Si(OH)(4) limitation after the initial rapid increase. Increasing mesozooplankton grazing rate, the diatom increase due to iron addition stays at minimum level, but small phytoplankton tend to increase. The numerical model experiments demonstrate the value of ecosystem modeling for evaluating the detailed interaction between biogeochemical cycle and iron fertilization in the equatorial Pacific.

(Table 1) Water column biogenic silica concentrations and Si isotopic composition during Marion Dufresne cruise MD166 off South Africa

Southern Ocean biogeochemical processes have an impact on global marine primary production and global elemental cycling, e.g. by likely controlling glacial-interglacial pCO2 variation. In this context, the natural silicon isotopic composition (d30Si) of sedimentary biogenic silica has been used to reconstruct past Si-consumption:supply ratios in the surface waters. We present a new dataset in the Southern Ocean from a IPY-GEOTRACES transect (Bonus-GoodHope) which includes for the first time summer d30Si signatures of suspended biogenic silica (i) for the whole water column at three stations and (ii) in the mixed layer at seven stations from the subtropical zone up to the Weddell Gyre. In general, the isotopic composition of biogenic opal exported to depth was comparable to the opal leaving the mixed layer and did not seem to be affected by any diagenetic processes during settling, even if an effect of biogenic silica dissolution cannot be ruled out in the northern part of the Weddell Gyre. We develop a mechanistic understanding of the processes involved in the modern Si-isotopic balance, by implementing a mixed layer model. We observe that the accumulated biogenic silica (sensu Rayleigh distillation) should satisfactorily describe the d30Si composition of biogenic silica exported out of the mixed layer, within the limit of the current analytical precision on the d30Si. The failures of previous models (Rayleigh and steady state) become apparent especially at the end of the productive period in the mixed layer, when biogenic silica production and export are low. This results from (1) a higher biogenic silica dissolution:production ratio imposing a lower net fractionation factor and (2) a higher Si-supply:Si-uptake ratio supplying light Si-isotopes into the mixed layer. The latter effect is especially expressed when the summer mixed layer becomes strongly Si-depleted, together with a large vertical silicic acid gradient, e.g. in the Polar Front Zone and at the Polar Front.