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.

Animal Guts as Nonideal Chemical Reactors: Partial Mixing and Axial Variation in Absorption Kinetics

Animal guts have been idealized as axially uniform plug-flow reactors (PFRs) without significant axial mixing or as combinations in series of such PFRs with other reactor types. To relax these often unrealistic assumptions and to provide a means for relaxing others, I approximated an animal gut as a series of n continuously stirred tank reactors (CSTRs) and examined its performance as a Function of n. For the digestion problem of hydrolysis and absorption in series, I suggest as a first approximation that a tubular gut of length L and diameter D comprises n=L/D tanks in series. For n greater than or equal to 10, there is little difference between performance of the nCSTR model and an ideal PFR in the coupled tasks of hydrolysis and absorption. Relatively thinner and longer guts, characteristic of animals feeding on poorer forage, prove more efficient in both conversion and absorption by restricting axial mixing, in the same total volume, they also give a higher rate of absorption. I then asked how a fixed number of absorptive sites should be distributed among the n compartments. Absorption rate generally is maximized when absorbers are concentrated in the hindmost few compartments, but high food quality or suboptimal ingestion rates decrease the advantage of highly concentrated absorbers. This modeling approach connects gut function and structure at multiple scales and can be extended to include other nonideal reactor behaviors observed in real animals.