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.

North Pacific Carbon Cycle Response to Climate Variability on Seasonal to Decadal Timescales

Climate variability drives significant changes in the physical state of the North Pacific, and there may be important impacts of this variability on the upper ocean carbon balance across the basin. We address this issue by considering the response of seven biogeochemical ocean models to climate variability in the North Pacific. The models' upper ocean pCO(2) and air-sea CO(2) flux respond similarly to climate variability on seasonal to decadal timescales. Modeled seasonal cycles of pCO(2) and its temperature- and non-temperature-driven components at three contrasting oceanographic sites capture the basic features found in observations (Takahashi et al., 2002, 2006; Keeling et al., 2004; Brix et al., 2004). However, particularly in the Western Subarctic Gyre, the models have difficulty representing the temporal structure of the total pCO(2) seasonal cycle because it results from the difference of these two large and opposing components. In all but one model, the air-sea CO(2) flux interannual variability (1 sigma) in the North Pacific is smaller ( ranges across models from 0.03 to 0.11 PgC/yr) than in the Tropical Pacific ( ranges across models from 0.08 to 0.19 PgC/yr), and the time series of the first or second EOF of the air-sea CO(2) flux has a significant correlation with the Pacific Decadal Oscillation (PDO). Though air-sea CO(2) flux anomalies are correlated with the PDO, their magnitudes are small ( up to +/- 0.025 PgC/yr ( 1 sigma)). Flux anomalies are damped because anomalies in the key drivers of pCO(2) ( temperature, dissolved inorganic carbon (DIC), and alkalinity) are all of similar magnitude and have strongly opposing effects that damp total pCO(2) anomalies.

Physical and Biological Controls on the Latitudinal Asymmetry of Surface Nutrients and pCO(2) in the Central and Eastern Equatorial Pacific

Surface nutrients and dissolved inorganic carbon (DIC) in the central (CEP) and eastern equatorial Pacific (EEP) show much higher concentrations to the south than to the north of the equator. In this study, the physical and biological controls on this asymmetry are investigated using a coupled physical-biogeochemical model. Two numerical experiments are conducted to examine the effects of asymmetrical photosynthetic efficiency (a) due to asymmetrical iron supply about the equator. The experiment with asymmetrical photosynthesis produces improved results as compared with historical observations. A nitrate budget analysis suggests that in the EEP the divergence of upwelling waters controls the surface nitrate asymmetry with additional contribution from the South Equatorial Current (SEC) carrying nutrient-rich Peru upwelling water. The changes of a affect the surface nitrate distribution but not the overall asymmetry. The SEC further carries excess nitrate to the west and thus extends the asymmetry in the east to the CEP. In the CEP, however, stronger northward than southward transport tends to reduce the nitrate asymmetry, while the asymmetrical photosynthesis would help to maintain it. Similar processes also control the distributions of surface silicate and DIC in the equatorial Pacific, which is also affected by the air-sea CO(2) exchange. The asymmetrical photosynthesis influences the distribution of surface DIC, pCO(2), and the air-sea CO(2) flux, by redistributing about 20% CO(2) flux from the north to the south of the equator. Owing to the adjustment of air-sea CO(2) flux, however, the net surface DIC change is smaller than the direct change associated with primary production.