Regulation of Phytoplankton Carbon to Chlorophyll Ratio By Light, Nutrients and Temperature in the Equatorial Pacific Ocean: A Basin-Scale Model

The complex effects of light, nutrients and temperature lead to a variable carbon to chlorophyll (C:Chl) ratio in phytoplankton cells. Using field data collected in the Equatorial Pacific, we derived a new dynamic model with a non-steady C:Chl ratio as a function of irradiance, nitrate, iron, and temperature. The dynamic model is implemented into a basin-scale ocean circulation-biogeochemistry model and tested in the Equatorial Pacific Ocean. The model reproduces well the general features of phytoplankton dynamics in this region. For instance, the simulated deep chlorophyll maximum (DCM) is much deeper in the western warm pool (similar to 100 m) than in the Eastern Equatorial Pacific (similar to 50 m). The model also shows the ability to reproduce chlorophyll, including not only the zonal, meridional and vertical variations, but also the interannual variability. This modeling study demonstrates that combination of nitrate and iron regulates the spatial and temporal variations in the phytoplankton C:Chl ratio in the Equatorial Pacific. Sensitivity simulations suggest that nitrate is mainly responsible for the high C:Chl ratio in the western warm pool while iron is responsible for the frontal features in the C:Chl ratio between the warm pool and the upwelling region. In addition, iron plays a dominant role in regulating the spatial and temporal variations of the C:Chl ratio in the Central and Eastern Equatorial Pacific. While temperature has a relatively small effect on the C:Chl ratio, light is primarily responsible for the vertical decrease of phytoplankton C:Chl ratio in the euphotic zone.

Significant Contribution of Large Particles to Optical Backscattering in the Open Ocean

The light scattering properties of oceanic particles have been suggested as an alternative index of phytoplankton biomass than chlorophyll-a concentration (chl-a), with the benefit of being less sensitive to physiological forcings (e.g., light and nutrients) that alter the intracellular pigment concentrations. The drawback of particulate scattering is that it is not unique to phytoplankton. Nevertheless, field studies have demonstrated that, to first order, the particulate beam-attenuation coefficient (c(p)) can track phytoplankton biomass. The relationship between c(p) and the particulate backscattering coefficient (b(bp)), a property retrievable from space, has not been fully evaluated, largely due to a lack of open-ocean field observations. Here, we present extensive data on inherent optical properties from the Equatorial Pacific surface waters and demonstrate a remarkable coherence in b(bp) and c(p). Coincident measurements of particle size distributions (PSDs) and optical properties of size-fractionated samples indicate that this covariance is due to both the conserved nature of the PSD and a greater contribution of phytoplankton-sized particles to b(bp) than theoretically predicted. These findings suggest that satellite-derived b(bp)could provide similar information on phytoplankton biomass in the open ocean as c(p).

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.

Phytoplankton Scales of Variability in the California Current System: 1. Interannual and Cross-Shelf Variability

In the California Current System, strong mesoscale variability associated with eddies and meanders of the coastal jet play an important role in the biological productivity of the area. To assess the dominant timescales of variability, a wavelet analysis is applied to almost nine years (October 1997 to July 2006) of 1-km-resolution, 5-day-averaged, Sea-viewing Wide Field-of-view Sensor (SeaWiFS) chlorophyll a (chl a) concentration data. The dominant periods of chlorophyll variance, and how these change in time, are quantified as a function of distance offshore. The maximum variance in chlorophyll occurs with a period of similar to 100-200 days. A seasonal cycle in the timing of peak variance is revealed, with maxima in spring/summer close to shore (20 km) and in autumn/winter 200 km offshore. Interannual variability in the magnitude of chlorophyll variance shows maxima in 1999, 2001, 2002, and 2005. There is a very strong out-of-phase correspondence between the time series of chlorophyll variance and the Pacific Decadal Oscillation (PDO) index. We hypothesize that positive PDO conditions, which reflect weak winds and poor upwelling conditions, result in reduced mesoscale variability in the coastal region, and a subsequent decrease in chlorophyll variance. Although the chlorophyll variance responds to basin-scale forcing, chlorophyll biomass does not necessarily correspond to the phase of the PDO, suggesting that it is influenced more by local-scale processes. The mesoscale variability in the system may be as important as the chl a biomass in determining the potential productivity of higher trophic levels.

Carbon-Based Ocean Productivity and Phytoplankton Physiology from Space

Ocean biogeochemical and ecosystem processes are linked by net primary production (NPP) in the ocean's surface layer, where inorganic carbon is fixed by photosynthetic processes. Determinations of NPP are necessarily a function of phytoplankton biomass and its physiological status, but the estimation of these two terms from space has remained an elusive target. Here we present new satellite ocean color observations of phytoplankton carbon (C) and chlorophyll (Chl) biomass and show that derived Chl:C ratios closely follow anticipated physiological dependencies on light, nutrients, and temperature. With this new information, global estimates of phytoplankton growth rates (mu) and carbon-based NPP are made for the first time. Compared to an earlier chlorophyll-based approach, our carbon-based values are considerably higher in tropical oceans, show greater seasonality at middle and high latitudes, and illustrate important differences in the formation and demise of regional algal blooms. This fusion of emerging concepts from the phycological and remote sensing disciplines has the potential to fundamentally change how we model and observe carbon cycling in the global oceans.

Interannual Variability in Timing of Bloom Initiation in the California Current System

In the California Current System the spring transition from poleward to equatorward alongshore wind stress heralds the beginning of upwelling-favorable conditions. The phytoplankton response to this transition is investigated using 8 years ( 1998-2005) of daily, 4-km resolution, Sea-viewing Wide Field of view Sensor ( SeaWiFS) chlorophyll a concentration data. Cluster analysis of the chlorophyll a time series at each location is used to separate the inshore upwelling region from offshore and oligotrophic areas. An objective method for estimating the timing of bloom initiation is used to construct a map of the mean bloom start date. Interannual variability in bloom timing and magnitude is investigated in four regions: 45 degrees N - 50 degrees N, 40 degrees N - 45 degrees N, 35 degrees N - 40 degrees N and 20 degrees N - 35 degrees N. Daily satellite derived wind data ( QuikSCAT) allow the timing of the first episode of persistently upwelling favorable winds to be estimated. Bloom initiation generally coincides with the onset of upwelling winds ( +/- 15 days). South of similar to 35 degrees N, where winds are southward year-round, the timing of increased chlorophyll concentration corresponds closely to timing of the seasonal increase in upwelling intensity. A 1-D model and satellite derived photosynthetically available radiation data are used to estimate time series of depth- averaged irradiance. In the far north of the region (> 46 degrees N) light is shown to limit phytoplankton growth in early spring. In 2005 the spring bloom in the northern regions (> 35 degrees N) had a "false start''. A sharp increase in chl a in February quickly receded, and a sustained increase in biomass was delayed until July. We hypothesize that this resulted in a mismatch in timing of food availability to higher trophic levels.

Seasonal Variability in the Bacteriolytic Capacity of the Deposit Feeder Arenicola Marina: Environmental Correlates

Although deposit-feeding macrofauna consume and digest sedimentary bacteria, it is unclear whether feeding rates and digestion efficiencies are high enough to significantly impact the composition and abundance of bacteria in marine sediments. It is likely that both feeding rates and efficiency of digestion vary markedly through space and time. We used a turbidimetric assay to compare the rate of bacteriolysis by digestive fluids collected seasonally from the deposit-feeding polychaete Arenicola marina. Under standardized, experimental conditions, bacteriolytic rates represent concentrations of lytic agents. This concentration was found to vary significantly throughout the year (p = 0.001), showing greater than a 2x range. Lytic agent concentration was positively correlated with bioavailable amino acid concentrations in the surface sediment (r = 0.85, p = 0.03) but showed no apparent relationship to other proxies for food resources (e.g, chl a), sediment temperature, or gut throughput time. In vitro, temperature has been shown to have a strong positive influence on bacteriolytic rate. Temperature has no influence, however, on the in situ concentration of lytic agent in gut fluids, thus it appears that compensation for this temperature dependence is unimportant. These findings, combined with previous kinetics studies with A. marina gut fluids, predict that the quantitative influence of deposit feeding on the microbial ecology of sediments will exhibit clear seasonal variation.

Physical-Biological Oceanographic Coupling Influencing Phytoplankton and Primary Production in the South China Sea

Two cruises were carried out in the summer and winter of 1998 to study coupled physical-chemical-biological processes in the South China Sea and their effects on phytoplankton stock and production. The results clearly show that the seasonal distributions of phytoplankton were closely related to the coupled processes driven by the East Asian Monsoon. Summer southwesterly monsoon induced upwelling along the China and Vietnam coasts. Several mesoscale cyclonic cold eddies and anticyclonic warm pools were identified in both seasons. In the summer, the upwelling and cold eddies, both associated with rich nutrients, low dissolved oxygen ( DO), high chlorophyll a (Chl a) and primary production ( PP), were found in the areas off the coast of central Vietnam, southeast of Hainan Island and north of the Sunda shelf, whereas in the winter they form a cold trough over the deep basin aligning from southwest to northeast. The warm pools with poor nutrients, high DO, low Chl a, and PP were found in the areas southeast of Vietnam, east of Hainan, and west of Luzon during the summer, and a northwestward warm jet from the Sulu Sea with properties similar to the warm pools was encountered during the winter. The phytoplankton stock and primary production were lower in summer due to nutrient depletion near the surface, particularly PO4. This phosphorus depletion resulted in phytoplankton species succession from diatoms to dinoflagellates and cyanophytes. A strong subsurface Chl a maximum, dominated by photosynthetic picoplankton, was found to contribute significantly to phytoplankton stocks and production.

Carbon-Based Primary Productivity Modeling With Vertically Resolved Photoacclimation

Net primary production (NPP) is commonly modeled as a function of chlorophyll concentration (Chl), even though it has been long recognized that variability in intracellular chlorophyll content from light acclimation and nutrient stress confounds the relationship between Chl and phytoplankton biomass. It was suggested previously that satellite estimates of backscattering can be related to phytoplankton carbon biomass (C) under conditions of a conserved particle size distribution or a relatively stable relationship between C and total particulate organic carbon. Together, C and Chl can be used to describe physiological state (through variations in Chl:C ratios) and NPP. Here, we fully develop the carbon-based productivity model (CbPM) to include information on the subsurface light field and nitracline depths to parameterize photoacclimation and nutrient stress throughout the water column. This depth-resolved approach produces profiles of biological properties (Chl, C, NPP) that are broadly consistent with observations. The CbPM is validated using regional in situ data sets of irradiance-derived products, phytoplankton chlorophyll: carbon ratios, and measured NPP rates. CbPM-based distributions of global NPP are significantly different in both space and time from previous Chl-based estimates because of the distinction between biomass and physiological influences on global Chl fields. The new model yields annual, areally integrated water column production of similar to 52 Pg C a(-1) for the global oceans.