Modeling studies investigating the causes of preferential depletion of silicic acid relative to nitrate during SERIES, a mesoscale iron enrichment in the NE subarctic Pacific

Numerical modeling experiments were conducted to examine the reasons for observed changes in the silicic acid ([Si(OH)4]) to nitrate ([NO3-]) drawdown ratio after the onset of algal iron stress during SERIES. During phytoplankton blooms and immediately after them, cells encounter a range of iron stress (between iron-replete and iron-deplete) and therefore show a range of growth rates. For these reasons, the potential influence of phytoplankton growth rate, under conditions of algal iron stress, on silicic acid and nitrate depletion were investigated in numerical experiments by altering the timing of a shift in the [Si(OH)4]:[NO3-] uptake ratio. These simulations suggested that the continued growth of iron-stressed phytoplankton at sub-maximum rates, with an elevated [Si(OH)4]:[NO3-] uptake ratio, induced depletion of silicic acid in the surface water and resulted in simultaneous limitation of growth by both iron and silicic-acid supply. Therefore, bottom-up control played an important role in terminating the phytoplankton bloom in SERIES. In the model simulations, the enhancement of diatom silicification due to increased rates of biomass-normalized silicic-acid uptake, led to increases in the export flux of opal after the onset of algal iron-stress and, consequently, it stimulated the silica pump. The regulation of both the [Si(OH)4]:[NO3-] uptake ratio and the growth rate of phytoplankton by iron supply are important factors that determine the relative consumption of silicic acid and nitrate upon iron stress, although the potential influence of a floristic shift in the diatom assemblage cannot be ruled out. These findings offer insights into the impact of iron fertilization, both artificial and natural, on the biogeochemical cycling of nutrients in high-nitrate, low-chlorophyll waters.