Search

Climate-driven depletion of ocean oxygen strongly impacts the global cycles of carbon and nutrients as well as the survival of many animal species. One of the main uncertainties in predicting changes to marine oxygen levels is the regulation of the biological respiration demand associated with the biological pump. Derived from the Redfield ratio, the molar ratio of oxygen to organic carbon consumed during respiration (i.e., the respiration quotient, [Formula: see text]) is consistently assumed constant but rarely, if ever, measured. Using a prognostic Earth system model, we show that a 0.1 increase in the respiration quotient from 1.0 leads to a 2.3% decline in global oxygen, a large expansion of low-oxygen zones, additional water column denitrification of 38 Tg N/y, and the loss of fixed nitrogen and carbon production in the ocean. We then present direct chemical measurements of [Formula: see text] using a Pacific Ocean meridional transect crossing all major surface biome types. The observed [Formula: see text] has a positive correlation with temperature, and regional mean values differ significantly from Redfield proportions. Finally, an independent global inverse model analysis constrained with nutrients, oxygen, and carbon concentrations supports a positive temperature dependence of [Formula: see text] in exported organic matter. We provide evidence against the common assumption of a static biological link between the respiration of organic carbon and the consumption of oxygen. Furthermore, the model simulations suggest that a changing respiration quotient will impact multiple biogeochemical cycles and that future warming can lead to more intense deoxygenation than previously anticipated.

Surface ocean phosphate is commonly below the standard analytical detection limits, leading to an incomplete picture of the global variation and biogeochemical role of phosphate. A global compilation of phosphate measured using high-sensitivity methods revealed several previously unrecognized low-phosphate areas and clear regional differences. Both observational climatologies and Earth system models (ESMs) systematically overestimated surface phosphate. Furthermore, ESMs misrepresented the relationships between phosphate, phytoplankton biomass, and primary productivity. Atmospheric iron input and nitrogen fixation are known important controls on surface phosphate, but model simulations showed that differences in the iron-to-macronutrient ratio in the vertical nutrient supply and surface lateral transport are additional drivers of phosphate concentrations. Our study demonstrates the importance of accurately quantifying nutrients for understanding the regulation of ocean ecosystems and biogeochemistry now and under future climate conditions.

Metal dissolution from atmospheric aerosol deposition to the oceans is important in enhancing and inhibiting phytoplankton growth rates and modifying plankton community structure, thus impacting marine biogeochemistry. Here we review the current state of knowledge on the causes and effects of the leaching of multiple trace metals from natural and anthropogenic aerosols. Aerosol deposition is considered both on short timescales over which phytoplankton respond directly to aerosol metal inputs, as well as longer timescales over which biogeochemical cycles are affected by aerosols.

We examine climate change impacts on net primary production (NPP) and export production (sinking particulate flux; EP) with simulations from nine Earth system models (ESMs) performed in the framework of the fifth phase of the Coupled Model Intercomparison Project (CMIP5). Global NPP and EP are reduced by the end of the century for the intense warming scenario of Representative Concentration Pathway (RCP) 8.5. Relative to the 1990s, NPP in the 2090s is reduced by 2-16 % and EP by 7-18 %. The models with the largest increases in stratification (and largest relative declines in NPP and EP) also show the largest positive biases in stratification for the contemporary period, suggesting overestimation of climate change impacts on NPP and EP. All of the CMIP5 models show an increase in stratification in response to surface-ocean warming and freshening, which is accompanied by decreases in surface nutrients, NPP and EP. There is considerable variability across the models in the magnitudes of NPP, EP, surface nutrient concentrations and their perturbations by climate change. The negative response of NPP and EP to increasing stratification reflects primarily a bottom-up control, as upward nutrient flux declines at the global scale. Models with dynamic phytoplankton community structure show larger declines in EP than in NPP. This pattern is driven by phytoplankton community composition shifts, with reductions in productivity by large phytoplankton as smaller phytoplankton (which export less efficiently) are favored under the increasing nutrient stress. Thus, the projections of the NPP response to climate change are critically dependent on the simulated phytoplankton community structure, the efficiency of the biological pump and the resulting levels of regenerated production, which vary widely across the models. Community structure is represented simply in the CMIP5 models, and should be expanded to better capture the spatial patterns and climate-driven changes in export efficiency.

A large database of field estimates of phytoplankton community growth rates in natural populations was compiled and analyzed to determine the apparent temperature effect on phytoplankton community growth rate. We conducted an ordinary least squares regression to optimize the parameters in two commonly used growth-temperature relations (Arrhenius and Q10 models). Both equations fit the observational data equally with the optimized parameter values. The optimum apparent Q10 value was 1.47 ± 0.08 (95% confidence interval, CI). Microzooplankton grazing rates closely matched the temperature trends for phytoplankton growth. This likely reflects a dynamic adjustment of biomass and grazing rates by the microzooplankton to match their available food source, illustrating tight coupling of phytoplankton growth and microzooplankton grazing rates. The field-measured temperature effect and growth rates were compared with estimates from the satellite Carbon-based Productivity Model (CbPM) and three Earth System Models (ESMs), with model output extracted at the same month and sampling locations as the observations. The optimized, apparent Q10 value calculated for the CbPM was 1.51, with overestimation of growth rates. The apparent Q10 value in the Community Earth System Model (V1.0) was 1.65, with modest underestimation of growth rates. The GFDL-ESM2M and GFDL-ESM2G models produced apparent Q10 values of 1.52 and 1.39, respectively. Models with an apparent Q10 that is significantly greater than ~1.5 will overestimate the phytoplankton community growth response to the ongoing climate warming and will have spatial biases in estimated growth rates for the current era.

Numerical models of ocean biogeochemistry are relied upon to make projections about the impact of climate change on marine resources and test hypotheses regarding the drivers of past changes in climate and ecosystems. In large areas of the ocean, iron availability regulates the functioning of marine ecosystems and hence the ocean carbon cycle. Accordingly, our ability to quantify the drivers and impacts of fluctuations in ocean ecosystems and carbon cycling in space and time relies on first achieving an appropriate representation of the modern marine iron cycle in models. When the iron distributions from 13 global ocean biogeochemistry models are compared against the latest oceanic sections from the GEOTRACES program, we find that all models struggle to reproduce many aspects of the observed spatial patterns. Models that reflect the emerging evidence for multiple iron sources or subtleties of its internal cycling perform much better in capturing observed features than their simpler contemporaries, particularly in the ocean interior. We show that the substantial uncertainty in the input fluxes of iron results in a very wide range of residence times across models, which has implications for the response of ecosystems and global carbon cycling to perturbations. Given this large uncertainty, iron fertilization experiments based on any single current generation model should be interpreted with caution. Improvements to how such models represent iron scavenging and also biological cycling are needed to raise confidence in their projections of global biogeochemical change in the ocean.

Organic macromolecules constitute a high percentage of remote sea spray. They enter the atmosphere through adsorption onto bubbles followed by bursting at the ocean surface, and go on to influence the chemistry of the fine mode aerosol. We present a global estimate of mixed-layer macromolecular distributions, driven by offline marine systems model output. The approach permits estimation of oceanic concentrations and bubble film surface coverages for several classes of organic compound. Mixed layer levels are computed from the output of a global ocean ecodynamics model by relating the macromolecules to standard biogeochemical tracers. Steady state is assumed for labile forms, and for longer-lived components we rely on ratios to existing transported variables. Adsorption is then represented through conventional Langmuir isotherms, with equilibria deduced from laboratory analogs. Open water concentrations locally exceed one micromolar carbon for the total of proteins, polysaccharides and refractory heteropolycondensates. The shorter-lived lipids remain confined to regions of strong biological activity. Results are evaluated against available measurements for all compound types, and agreement is generally well within an order of magnitude. Global distributions are further estimated for both fractional coverage of bubble films at the air–water interface and the two-dimensional concentration excess. Overall, we show that macromolecular mapping provides a novel tool for the comprehension of oceanic surfactant patterns. These results may prove useful in planning field experiments and assessing the potential response of surface chemical behaviors to global change.

Selective removal of nitrogen (N) and phosphorus (P) from the marine dissolved organic matter (DOM) pool has been reported in several regional studies. Because DOM is an important advective/mixing pathway of carbon (C) export from the ocean surface layer and its non-Redfieldian stoichiometry would affect estimates of marine export production per unit N and P, we investigated the stoichiometry of marine DOM and its remineralization globally using a compiled DOM data set. Marine DOM is enriched in C and N compared to Redfield stoichiometry, averaging 317:39:1 and 810:48:1 for C:N:P within the degradable and total bulk pools, respectively. Dissolved organic phosphorus (DOP) is found to be preferentially remineralized about twice as rapidly with respect to the enriched C:N stoichiometry of marine DOM. Biogeochemical simulations with the Biogeochemical Elemental Cycling model using Redfield and variable DOM stoichiometry corroborate the need for non-Redfield dynamics to match the observed DOM stoichiometry. From our model simulations, preferential DOP remineralization is found to increase the strength of the biological pump by ∼9% versus the case of Redfield DOM cycling. Global net primary productivity increases ∼10% including an increase in marine nitrogen fixation of ∼26% when preferential DOP remineralization and direct utilization of DOP by phytoplankton are included. The largest increases in marine nitrogen fixation, net primary productivity, and carbon export are observed within the western subtropical gyres, suggesting the lateral transfer of P in the form of DOP from the productive eastern and poleward gyre margins may be important for sustaining these processes downstream in the subtropical gyres.

The ratio of carbon (C) to phosphorus (P) in marine phytoplankton is thought to be constant throughout the worlds'oceans. Known as the Redfield ratio, this relationship describes the links between carbon and phosphorus cycling and marine ecosystems. However, variations in the stoichiometry of phytoplankton have recently been identified, in particular strong latitudinal variability. Here we assess the impact of this variability in the C:P ratio of biomass on the C:P ratio of organic matter that is exported to the deep ocean using a biogeochemical inverse-model based on a data-constrained ocean circulation model and a global database of dissolved inorganic carbon and phosphate measurements. We identify global patterns of variability in the C:P ratios of exported organic matter, with higher values in the nutrient-depleted subtropical gyres, where organic matter export is relatively low, and lower ratios in nutrient-rich upwelling zones and high-latitude regions, where organic matter export is high. This suggests that total carbon export is relatively constant throughout the oceans, in agreement with recent estimates of carbon fluxes. We conclude that the latitudinal patterns of C:P in exported organic matter are consistent with the large-scale stoichiometric variations in phytoplankton C:P. We suggest that a future expansion of nutrient-depleted waters could result in a shift to more efficient C export that compensates for the expected decline in productivity.

The authors compare Community Earth System Model results to marine observations for the 1990s and examine climate change impacts on biogeochemistry at the end of the twenty-first century under two future scenarios (Representative Concentration Pathways RCP4.5 and RCP8.5). Late-twentieth-century seasonally varying mixed layer depths are generally within 10 m of observations, with a Southern Ocean shallow bias. Surface nutrient and chlorophyll concentrations exhibit positive biases at low latitudes and negative biases at high latitudes. The volume of the oxygen minimum zones is overestimated.The impacts of climate change on biogeochemistry have similar spatial patterns under RCP4.5 and RCP8.5, but perturbation magnitudes are larger under RCP8.5. Increasing stratification leads to weaker nutrient entrainment and decreased primary and export production (>30% over large areas). The global-scale decreases in primary and export production scale linearly with the increases in mean sea surface temperature. There are production increases in the high nitrate, low chlorophyll (HNLC) regions, driven by lateral iron inputs from adjacent areas. The increased HNLC export partially compensates for the reductions in non-HNLC waters (~25% offset). Stabilizing greenhouse gas emissions and climate by the end of this century (as in RCP4.5) will minimize the changes to nutrient cycling and primary production in the oceans. In contrast, continued increasing emission of CO2 (as in RCP8.5) will lead to reduced productivity and significant modifications to ocean circulation and biogeochemistry by the end of this century, with more drastic changes beyond the year 2100 as the climate continues to rapidly warm.

Ocean carbon uptake and storage simulated by the Community Earth System Model, version 1–Biogeochemistry [CESM1(BGC)], is described and compared to observations. Fully coupled and ocean-ice configurations are examined; both capture many aspects of the spatial structure and seasonality of surface carbon fields. Nearly ubiquitous negative biases in surface alkalinity result from the prescribed carbonate dissolution profile. The modeled sea–air CO2 fluxes match observationally based estimates over much of the ocean; significant deviations appear in the Southern Ocean. Surface ocean pCO2 is biased high in the subantarctic and low in the sea ice zone. Formation of the water masses dominating anthropogenic CO2 (Cant) uptake in the Southern Hemisphere is weak in the model, leading to significant negative biases in Cant and chlorofluorocarbon (CFC) storage at intermediate depths. Column inventories of Cant appear too high, by contrast, in the North Atlantic. In spite of the positive bias, this marks an improvement over prior versions of the model, which underestimated North Atlantic uptake. The change in behavior is attributable to a new parameterization of density-driven overflows. CESM1(BGC) provides a relatively robust representation of the ocean–carbon cycle response to climate variability. Statistical metrics of modeled interannual variability in sea–air CO2 fluxes compare reasonably well to observationally based estimates. The carbon cycle response to key modes of climate variability is basically similar in the coupled and forced ocean-ice models; however, the two differ in regional detail and in the strength of teleconnections.

This study used an ocean general circulation model to simulate the marine iron cycle in an investigation of how simulated distributions of weak iron-binding ligands would be expected to control dissolved iron concentrations in the ocean, with a particular focus on deep ocean waters. The distribution of apparent oxygen utilization was used as a proxy for humic substances that have recently been hypothesized to account for the bulk of weak iron-binding ligands in seawater. Compared to simulations using a conventional approach with homogeneous ligand distributions, the simulations that incorporated spatially variable ligand concentrations exhibited substantial improvement in the simulation of global dissolved iron distributions as revealed by comparisons with available field data. The improved skill of the simulations resulted largely because the spatially variable ligand distributions led to a more reasonable basin-scale variation of the residence time of iron when present at high concentrations. The model results, in conjunction with evidence from recent field studies, suggest that humic substances play an important role in the iron cycle in the ocean.

Nearly 75 years ago, Alfred C. Redfield observed a similarity between the elemental composition of marine plankton in the surface ocean and dissolved nutrients in the ocean interior. This stoichiometry, referred to as the Redfield ratio, continues to be a central tenet in ocean biogeochemistry, and is used to infer a variety of ecosystem processes, such as phytoplankton productivity and rates of nitrogen fixation and loss2-4. Model, field and laboratory studies have shown that different mechanisms can explain both constant and variable ratios of carbon to nitrogen and phosphorus among ocean plankton communities. The range of C/N/P ratios in the ocean, and their predictability, are the subject of much active research. Here we assess global patterns in the elemental composition of phytoplankton and particulate organic matter in the upper ocean, using published and unpublished observations of particulate phosphorus, nitrogen and carbon from a broad latitudinal range, supplemented with elemental data for surface plankton populations. We show that the elemental ratios of marine organic matter exhibit large spatial variations, with a global average that differs substantially from the canonical Redfield ratio. However, elemental ratios exhibit a clear latitudinal trend. Specifically, we observed a ratio of 195:28:1 in the warm nutrient-depleted low-latitude gyres, 137:18:1 in warm, nutrient-rich upwelling zones, and 78:13:1 in cold, nutrient-rich high-latitude regions. We suggest that the coupling between oceanic carbon, nitrogen and phosphorus cycles may vary systematically by ecosystem. © 2013 Macmillan Publishers Limited. All rights reserved.

Mineral aerosol deposition is recognized as the dominant source of iron to the open ocean and the solubility of iron in the dust aerosol is highly variable, with measurements ranging from 0.01-80%. Global models have difficulty capturing the observed variations in solubility, and have ignored the solubility dependence on aerosol size. We introduce two idealized physical models to estimate the size dependence of mineral aerosol solubility: a diffusion-controlled model and a surface-area-controlled model. These models produce differing time- and space-varying solubility maps for aerosol Fe and Al given the dust age at deposition, size-resolved dust entrainment fields, and the aerosol acidity. The resulting soluble iron deposition fluxes are substantially different, and more realistic, than a globally uniform solubility approximation. The surface-area-controlled solubility varies more than the diffusion-controlled solubility and better captures the spatial pattern of observed solubility in the Atlantic. However, neither of these two models explains the large solubility variation observed in the Pacific. We then examine the impacts of spatially variable, size-dependent solubility on marine biogeochemistry with the Biogeochemical Elemental Cycling (BEC) ocean model by comparing the modeled surface ocean dissolved Fe and Al with observations. The diffusion-based variable solubility does not significantly improve the simulation of dissolved Fe relative to a 5% globally uniform solubility, while the surface-area-based variable solubility improves the simulation in the North Atlantic but worsens it in the Pacific and Indian Oceans.

Biogeochemical cycling in the Southern Ocean (SO) plays a key role in the global sea-air CO2balance and in the ocean anthropogenic carbon inventory (Ito et al., 2010; Khatiwala et al., 2009; Sarmiento et al., 2004). Some previous studies suggest a decreasing trend in the Southern Ocean carbon sink (Le Quéré et al., 2007; Lovenduski et al., 2007; Wetzel et al., 2005). We investigate the interannual and decadal variations in sea-air CO2flux and phytoplankton production in the SO with hindcast simulations by an ocean biogeochemical model. Decreasing trends in sinking POC and primary production are found from 1979 to 2003, concurrent with a decreasing trend in carbon uptake from the atmosphere. Simulations show substantial interannual and decadal variability in productivity. The sea-air CO2flux is significantly correlated with sinking POC, especially in high productivity regions of the Southern Ocean. Both mixed layer depths and iron concentrations are important to the long-term trends in production and phytoplankton community structure. Sea ice cover also plays an important role at high latitudes. Variability in dust deposition in recent decades has little influence on total SO productivity and carbon uptake, however, there are regional impacts near dust source regions. Accurately representing mixed layer depths and their impacts on phytoplankton light stress are critical for understanding how climate change impacts SO ecosystems and biogeochemistry.

Phaeocystis antarctica is an important phytoplankton species in the Southern Ocean. We incorporated P. antarctica into the biogeochemical elemental cycling ocean model to study Southern Ocean ecosystem dynamics and biogeochemistry. The optimum values of ecological parameters for Phaeocystis were sought through synthesizing laboratory and field observations, and the model output was evaluated with observed chlorophyll a, carbon biomass, and nutrient distributions. Several factors have been proposed to control Southern Ocean ecosystem structure, including light adaptation, iron uptake capability, and loss processes. Optimum simulation results were obtained when P. antarctica had a relatively high α (P-I curve initial slope) value and a higher half-saturation constant for iron uptake than other phytoplankton. Simulation results suggested that P. antarctica had a competitive advantage under low irradiance levels, especially in the Ross Sea and Weddell Sea. However, the distributions of P. antarctica and diatoms were also strongly influenced by iron availability. Although grazing rates had an influence on total biomass, our simulations did not show a strong influence of grazing pressure in the competition between P. antarctica and diatoms. However, limited observations and the relative simplicity of zooplankton in our model suggest further research is needed. Overall, P. antarctica contributed ∼13% of annual primary production and ∼19% of sinking carbon export in the Southern Ocean (>40°S) in our best case simulation. At higher latitudes (>60°S) P. antarctica accounts for ∼23% of annual primary production and ∼30% of sinking carbon export.

The performance of 36 models (22 ocean color models and 14 biogeochemical ocean circulation models (BOGCMs)) that estimate depth-integrated marine net primary productivity (NPP) was assessed by comparing their output to in situ 14C data at the Bermuda Atlantic Time series Study (BATS) and the Hawaii Ocean Time series (HOT) over nearly two decades. Specifically, skill was assessed based on the models' ability to estimate the observed mean, variability, and trends of NPP. At both sites, more than 90% of the models underestimated mean NPP, with the average bias of the BOGCMs being nearly twice that of the ocean color models. However, the difference in overall skill between the best BOGCM and the best ocean color model at each site was not significant. Between 1989 and 2007, in situ NPP at BATS and HOT increased by an average of nearly 2% per year and was positively correlated to the North Pacific Gyre Oscillation index. The majority of ocean color models produced in situ NPP trends that were closer to the observed trends when chlorophyll-a was derived from high-performance liquid chromatography (HPLC), rather than fluorometric or SeaWiFS data. However, this was a function of time such that average trend magnitude was more accurately estimated over longer time periods. Among BOGCMs, only two individual models successfully produced an increasing NPP trend (one model at each site). We caution against the use of models to assess multiannual changes in NPP over short time periods. Ocean color model estimates of NPP trends could improve if more high quality HPLC chlorophyll-a time series were available.

The primary nutrients that limit marine phytoplankton growth rates include nitrogen (N), phosphorus (P), iron (Fe), and silicon (Si). Atmospheric transport and deposition provides a source for each of these nutrients to the oceans. We utilize an ocean biogeochemical model to examine the relative importance of these atmospheric inputs for ocean biogeochemistry and export production. In the current era, simulations with the biogeochemical elemental cycling ocean model suggest that globally, atmospheric Fe inputs could support roughly 50% of the Fe exported from the euphotic zone by sinking organic and inorganic particles. Variations in atmospheric iron inputs strongly impact spatial patterns of phytoplankton growth limitation and the areal extent of the high-nutrient, low-chlorophyll regions. Atmospheric inputs of N, Si, and P have smaller impacts, potentially accounting for 5.1%, 0.21%, and 0.12% of the biogenic export of these elements from the euphotic zone, respectively. Soluble Fe input from the atmosphere is sufficient to support most of the export production in many ocean regions, whether we use a spatially variable aerosol Fe solubility, or a globally constant 2% solubility. Regionally atmospheric N inputs can have significant impacts on marine biogeochemistry, potentially supporting >25% of the export production, an impact that is increasing due to human activities. Atmospheric Si and P inputs have only minimal impacts on marine ecosystem productivity and biogeochemistry, as these inputs are typically quite small relative to the flux of these nutrients from below the euphotic zone.

We present results from transient sensitivity studies with the Biogeochemical Elemental Cycling (BEC) ocean model to increasing anthropogenic atmospheric inorganic nitrogen (N) and soluble iron (Fe) deposition over the industrial era. Elevated N deposition results from fossil fuel combustion and agriculture, and elevated soluble Fe deposition results from increased atmospheric processing in the presence of anthropogenic pollutants and soluble Fe from combustion sources. Simulations with increasing Fe and increasing Fe and N inputs raised simulated marine nitrogen fixation, with the majority of the increase in the subtropical North and South Pacific, and raised primary production and export in the high-nutrient low-chlorophyll (HNLC) regions. Increasing N inputs alone elevated small phytoplankton and diatom production, resulting in increased phosphorus (P) and Fe limitation for diazotrophs, hence reducing nitrogen fixation (∼6%). Globally, the simulated primary production, sinking particulate organic carbon (POC) export. and atmospheric CO2 uptake were highest under combined increase in Fe and N inputs compared to preindustrial control. Our results suggest that increasing combustion iron sources and aerosol Fe solubility along with atmospheric anthropogenic nitrogen deposition are perturbing marine biogeochemical cycling and could partially explain the observed trend toward increased P limitation at station ALOHA in the subtropical North Pacific. Excess inorganic nitrogen ([NO3 −] + [NH4 +] − 16[PO4 3−]) distributions may offer useful insights for understanding changing ocean circulation and biogeochemistry.

We quantify the mechanisms governing interannual variability in the global, upper-ocean inorganic carbon system using a hindcast simulation (1979-2004) of an ecosystem-biogeochemistry model forced with time-evolving atmospheric physics and dust deposition. We analyze the variability of three key, interrelated metrics-air-sea CO2 flux, surface-water carbon dioxide partial pressure pCO2, and upper-ocean dissolved inorganic carbon (DIC) inventory-presenting for each metric global spatial maps of the root mean square (rms) of anomalies from a model monthly climatology. The contribution of specific driving factors is diagnosed using Taylor expansions and linear regression analysis. The major regions of variability occur in the Southern Ocean, tropical Indo-Pacific, and Northern Hemisphere temperate and subpolar latitudes. Ocean circulation is the dominant factor driving variability over most of the ocean, modulating surface dissolved inorganic carbon that in turn alters surface-water pCO2 and air-sea CO2 flux variability (global integrated anomaly rms of 0.34 Pg C yr-1). Biological export and thermal solubility effects partially damp circulation-driven pCO2 variability in the tropics, while in the subtropics, thermal solubility contributes positively to surface-water pCO2 and air-sea CO2 flux variability. Gas transfer and net freshwater inputs induce variability in the air-sea CO2 flux in some specific regions. A component of air-sea CO2 flux variability (global integrated anomaly rms of 0.14 Pg C yr-1) arises from variations in biological export production induced by variations in atmospheric iron deposition downwind of dust source regions. Beginning in the mid-1990s, reduced global dust deposition generates increased air-sea CO2 outgassing in the Southern Ocean, consistent with trends derived from atmospheric CO2 inversions. © 2008 Elsevier Ltd.

Depth-integrated primary productivity (PP) estimates obtained from satellite ocean color-based models (SatPPMs) and those generated from biogeochemical ocean general circulation models (BOGCMs) represent a key resource for biogeochemical and ecological studies at global as well as regional scales. Calibration and validation of these PP models are not straightforward, however, and comparative studies show large differences between model estimates. The goal of this paper is to compare PP estimates obtained from 30 different models (21 SatPPMs and 9 BOGCMs) to a tropical Pacific PP database consisting of ∼ 1000 14C measurements spanning more than a decade (1983-1996). Primary findings include: skill varied significantly between models, but performance was not a function of model complexity or type (i.e. SatPPM vs. BOGCM); nearly all models underestimated the observed variance of PP, specifically yielding too few low PP (< 0.2 g C m- 2 d- 1) values; more than half of the total root-mean-squared model-data differences associated with the satellite-based PP models might be accounted for by uncertainties in the input variables and/or the PP data; and the tropical Pacific database captures a broad scale shift from low biomass-normalized productivity in the 1980s to higher biomass-normalized productivity in the 1990s, which was not successfully captured by any of the models. This latter result suggests that interdecadal and global changes will be a significant challenge for both SatPPMs and BOGCMs. Finally, average root-mean-squared differences between in situ PP data on the equator at 140°W and PP estimates from the satellite-based productivity models were 58% lower than analogous values computed in a previous PP model comparison 6 years ago. The success of these types of comparison exercises is illustrated by the continual modification and improvement of the participating models and the resulting increase in model skill. © 2008 Elsevier B.V.

We present a generalized framework for assessing the skill of global upper ocean ecosystem-biogeochemical models against in-situ field data and satellite observations. We illustrate the approach utilizing a multi-decade (1979-2004) hindcast experiment conducted with the Community Climate System Model (CCSM-3) ocean carbon model. The CCSM-3 ocean carbon model incorporates a multi-nutrient, multi-phytoplankton functional group ecosystem module coupled with a carbon, oxygen, nitrogen, phosphorus, silicon, and iron biogeochemistry module embedded in a global, three-dimensional ocean general circulation model. The model is forced with physical climate forcing from atmospheric reanalysis and satellite data products and time-varying atmospheric dust deposition. Data-based skill metrics are used to evaluate the simulated time-mean spatial patterns, seasonal cycle amplitude and phase, and subannual to interannual variability. Evaluation data include: sea surface temperature and mixed layer depth; satellite-derived surface ocean chlorophyll, primary productivity, phytoplankton growth rate and carbon biomass; large-scale climatologies of surface nutrients, pCO2, and air-sea CO2 and O2 flux; and time-series data from the Joint Global Ocean Flux Study (JGOFS). Where the data is sufficient, we construct quantitative skill metrics using: model-data residuals, time-space correlation, root mean square error, and Taylor diagrams. © 2008 Elsevier B.V. All rights reserved.

We use measurements of ocean surface dissolved Al, a global Biogeochemical Elemental Cycling (BEC) ocean model, and the global Dust Entrainment and Deposition (DEAD) model to constrain dust deposition to the oceans. Our Al database contains all available measurements with best coverage in the Atlantic. Vertical profiles and seasonal data exist in limited regions. Observations show that surface dissolved Al is distributed similarly to the dust deposition predicted by DEAD and other models. There is an equatorial Atlantic Al maximum that decreases toward higher latitudes. There are high Al concentrations in the Mediterranean Sea and the Arabian Sea and low concentrations in the Pacific and the Southern Ocean. The ocean basins maintain more distinct Al profiles than Fe profiles in the upper ocean, consistent with a weaker biological influence on Al than Fe. The BEC-predicted surface dissolved Al compares relatively well with observations. The Al distribution reflects the combined effects of Al input from dust and Al removal by particle scavenging and biological uptake by diatoms. Model-observed biases suggest a southward shift of maximum dust deposition compared to current dust model predictions. DEAD appears to overestimate deposition north of 30°N in the Pacific and to underestimate deposition south of 30°N. Observed Al concentrations and the ocean model–predicted surface Al lifetime provide a semi-independent method to estimate oceanic dust deposition. This technique indicates that DEAD may overestimate dust deposition to the north equatorial Atlantic but underestimate in other Atlantic regions, the Southern Ocean, and the Arabian Sea. However, spatial variations in aerosol Al solubility may also contribute to the model-observation mismatch. Our results have implications for all dust-borne ocean nutrients including Fe and demonstrate the potential of marine geochemical data to constrain atmospheric aerosol deposition fields.

To better understand the role of iron in driving marine ecosystems, the Southern Ocean Iron Experiment (SOFeX) fertilized two surface water patches with iron north and south of the Antarctic Polar Front Zone (APFZ). Using 1-D coupled biological-physical simulations, we examine the biogeochemical dynamics that occurred both inside and outside of the fertilized patches during and shortly after the SOFeX field campaign. We focus, in particular, on three main issues governing the biological response to deliberate iron fertilization: the interaction among phytoplankton, light, macronutrient and iron limitation; dilution and lateral mixing between the fertilized patch and external, unfertilized waters; and the effect of varying mixed layer depth on the light field. At the patch south of the APFZ, sensitivity simulations with no dilution results in the maximum bloom magnitude, whereas dilution with external water extends the development of the north patch bloom by relieving silicon limitation. In model sensitivity studies for both sites, maximum chlorophyll concentration and dissolved inorganic carbon depletion inside the fertilized patches are inversely related to mixed layer depth, similar to the patterns observed across a number of iron fertilization field experiments. Our results suggest that Southern Ocean phytoplankton blooms resulting from natural or deliberate iron fertilization will tend to become iron-light co-limited unless the mixed layer depth is quite shallow. © 2007 Elsevier B.V. All rights reserved.

Analysis of satellite ocean color and wind speed data within the seasonal ice zone (SIZ) of the Southern Ocean sheds light on the physical processes that influence phytoplankton biomass distributions. A compilation of monthly averaged chlorophyll and percent sea ice cover data within the SIZ from 1997 to 2005 has been compared with monthly average wind speed data from 1999 to 2005. The size of the marginal ice zone (MIZ, areas of recent ice retreat) was fairly consistent from year to year, always peaking in December, with a mean area of 6.0 million km2. The mean area within the MIZ with phytoplankton blooms (chlorophyll exceeding 0.8 mg/m3) was 0.36 million km2. While the bloom areal extent seems small compared to the MIZ, in reality, because of gaps in the chlorophyll data, blooming regions comprise a much larger fraction of the MIZ. Considering only areas with valid chlorophyll data, the percentage of the MIZ with blooms was 17%, 21%, and 24% for the months of December, January, and February, respectively. December always has the largest MIZ area, but MIZ mean chlorophyll concentrations sometimes do not peak until February. Wind speed strongly impacts phytoplankton bloom dynamics within the MIZ. There is an inverse relation between wind speed and bloom occurrence, with blooms largely suppressed at high wind speeds. At low wind speeds (∼5 m/s), blooms are observed over about one third of the MIZ. Blooms are also much more frequent near the continent than in offshore waters, likely due to increased iron availability. Open ocean phytoplankton blooms in the Southern Ocean are likely to become iron-light co-limited except in regions where the mixed layer depth is relatively shallow.

Application of biogeochemical models to the study of marine ecosystems is pervasive, yet objective quantification of these models' performance is rare. Here, 12 lower trophic level models of varying complexity are objectively assessed in two distinct regions (equatorial Pacific and Arabian Sea). Each model was run within an identical one-dimensional physical framework. A consistent variational adjoint implementation assimilating chlorophyll-a, nitrate, export, and primary productivity was applied and the same metrics were used to assess model skill. Experiments were performed in which data were assimilated from each site individually and from both sites simultaneously. A cross-validation experiment was also conducted whereby data were assimilated from one site and the resulting optimal parameters were used to generate a simulation for the second site. When a single pelagic regime is considered, the simplest models fit the data as well as those with multiple phytoplankton functional groups. However, those with multiple phytoplankton functional groups produced lower misfits when the models are required to simulate both regimes using identical parameter values. The cross-validation experiments revealed that as long as only a few key biogeochemical parameters were optimized, the models with greater phytoplankton complexity were generally more portable. Furthermore, models with multiple zooplankton compartments did not necessarily outperform models with single zooplankton compartments, even when zooplankton biomass data are assimilated. Finally, even when different models produced similar least squares model-data misfits, they often did so via very different element flow pathways, highlighting the need for more comprehensive data sets that uniquely constrain these pathways.

Recent upward revisions in key sink/source terms for fixed nitrogen (N) in the oceans imply a short residence time and strong negative feedbacks involving denitrification and N fixation to prevent large swings in the ocean N inventory over timescales of a few centuries. We tested the strength of these feedbacks in a global biogeochemical elemental cycling (BEC) ocean model that includes water column denitrification and an explicit N fixing phytoplankton group. In the northern Indian Ocean and over longer timescales in the tropical Atlantic, we find strong stabilizing feedbacks that minimize changes in marine N inventory over timescales of ∼30–200 years. In these regions high atmospheric dust/iron inputs lead to phosphorus limitation of diazotrophs, and thus a tight link between N fixation and surface water N/P ratios. Maintenance of the oxygen minimum zones in these basins depends on N fixation driven export. The stabilizing feedbacks in other regions are significant but weaker owing to iron limitation of the diazotrophs. Thus Fe limitation appears to restrict the ability of N fixation to compensate for changes in denitrification in the current climate, perhaps leading the oceans to lose fixed N. We suggest that iron is the ultimate limiting nutrient leading to nitrogen being the proximate limiting nutrient over wide regions today. Iron stress was at least partially alleviated during more dusty, glacial times, leading to a higher marine N inventory, increased export production, and perhaps widespread phosphorus limitation of the phytoplankton community. The increased efficiency of the biological pump would have contributed to the glacial drawdown in atmospheric CO2.

We perform a sensitivity study with the Biogeochemical Elemental Cycling (BEC) ocean model to understand the impact of atmospheric inorganic nitrogen deposition on marine biogeochemistry and air-sea CO2 exchange. Simulations involved examining the response to three different atmospheric inorganic nitrogen deposition scenarios namely, Pre-industrial (22 Tg N/year), 1990s (39 Tg N/year), and an Intergovernmental Panel on Climate Change (IPCC) prediction for 2100, IPCC-A1FI (69 Tg N/year). Globally, the increasing N deposition had widespread, but modest effects on export production and air-sea CO2 exchange. The maximum increase in N deposition was 47 Tg N/year since Pre-industrial control for the IPCC-A1FI case, which had an increase in primary production (0.98 Gt C/year or 2%), export production (0.16 Gt C/year or 3%) and a decrease in atmospheric pCO2 of 1.66 ppm (0.6%) relative to the Pre-industrial control. In some regions, atmospheric N inputs supported >20% of the export production in the current era and >50% of the export production in the IPCC-A1FI case. As nitrogen deposition increased, N fixation decreased because the diazotrophs were outcompeted by diatoms and small phytoplankton under more N-replete conditions. This decrease in N fixation could partially counteract the ongoing increase in new nitrogen inputs via atmospheric N deposition.

Estimates of sources/sinks of carbon dioxide (CO2) at the Earth's surface are commonly made using atmospheric CO2 inverse modeling, terrestrial and oceanic biogeochemical modeling, and inventory-based studies. In this study, we compare sea-air CO2 fluxes from the Time-Dependent Inverse (TDI) atmosphere model and the marine Biogeochemical Elemental Cycling (BEC) model to study the processes involved in ocean carbon cycling at subbasin scales. A dust generation and transport model, based on analyzed meteorology and terrestrial vegetation cover, is also used to estimate the interannual variability in dust and iron deposition to different ocean basins. Overall, a fairly good agreement is established between the TDI and BEC model results for the net annual patterns and seasonal cycle of sea-air CO2 exchange. Sensitivity studies with the ocean biogeochemical model using increased or reduced atmospheric iron inputs indicate the relative sensitivity of air-sea CO2 exchange. The simulated responses to changes in iron inputs are not instantaneous (peak response after ∼2−3 years). The TDI model derived seasonal cycles for the Southern Ocean (South Atlantic) are better matched by the BEC model by increasing (decreasing) iron inputs through atmospheric aerosols. Our results suggest that some of the interannual variability in TDI model air-sea CO2 fluxes during the past decade may be explainable by dust variability that relaxes/increases iron limitation in high-nitrate, low-chlorophyll (HNLC) ocean regions.

A global ocean biogeochemical model is used to quantify the sensitivity of marine biogeochemistry and air–sea CO2 exchange to variations in dust deposition over decadal timescales. Estimates of dust deposition generated under four climate states provide a large range in total deposition with spatially realistic patterns; transient ocean model experiments are conducted by applying a step-function change in deposition from a current climate control. Relative to current conditions, higher dust deposition increases diatom and export production, nitrogen fixation and oceanic net CO2 uptake from the atmosphere, while reduced dust deposition has the opposite effects. Over timescales less than a decade, dust modulation of marine productivity and export is dominated by direct effects in high-nutrient, low-chlorophyll regions, where iron is the primary limiting nutrient. On longer timescales, an indirect nitrogen fixation pathway has increased importance, significantly amplifying the ocean biogeochemical response. Because dust iron input decouples carbon cycling from subsurface macronutrient supply, the ratio of the change in net ocean CO2 uptake to change in export flux is large, 0.45–0.6. Decreasing dust deposition and reduced oceanic CO2 uptake over the next century could provide a positive feedback to global warming, distinct from feedbacks associated with changes in stratification and circulation.

In this paper, we review the state of the art and major challenges in current efforts to incorporate biogeochemical functional groups into models that can be applied on basin-wide and global scales, with an emphasis on models that might ultimately be used to predict how biogeochemical cycles in the ocean will respond to global warming. We define the term "biogeochemical functional group" to refer to groups of organisms that mediate specific chemical reactions in the ocean. Thus, according to this definition, "functional groups" have no phylogenetic meaning-these are composed of many different species with common biogeochemical functions. Substantial progress has been made in the last decade toward quantifying the rates of these various functions and understanding the factors that control them. For some of these groups, we have developed fairly sophisticated models that incorporate this understanding, e.g. for diazotrophs (e.g. Trichodesmium), silica producers (diatoms) and calcifiers (e.g. coccolithophorids and specifically Emiliania huxleyi). However, current representations of nitrogen fixation and calcification are incomplete, i.e., based primarily upon models of Trichodesmium and E. huxleyi, respectively, and many important functional groups have not yet been considered in open-ocean biogeochemical models. Progress has been made over the last decade in efforts to simulate dimethylsulfide (DMS) production and cycling (i.e., by dinoflagellates and prymnesiophytes) and denitrification, but these efforts are still in their infancy, and many significant problems remain. One obvious gap is that virtually all functional group modeling efforts have focused on autotrophic microbes, while higher trophic levels have been completely ignored. It appears that in some cases (e.g., calcification), incorporating higher trophic levels may be essential not only for representing a particular biogeochemical reaction, but also for modeling export. Another serious problem is our tendency to model the organisms for which we have the most validation data (e.g., E. huxleyi and Trichodesmium) even when they may represent only a fraction of the biogeochemical functional group we are trying to represent. When we step back and look at the paleo-oceanographic record, it suggests that oxygen concentrations have played a central role in the evolution and emergence of many of the key functional groups that influence biogeochemical cycles in the present-day ocean. However, more subtle effects are likely to be important over the next century like changes in silicate supply or turbulence that can influence the relative success of diatoms versus dinoflagellates, coccolithophorids and diazotrophs. In general, inferences drawn from the paleo-oceanographic record and theoretical work suggest that global warming will tend to favor the latter because it will give rise to increased stratification. However, decreases in pH and Fe supply could adversely impact coccolithophorids and diazotrophs in the future. It may be necessary to include explicit dynamic representations of nitrogen fixation, denitrification, silicification and calcification in our models if our goal is predicting the oceanic carbon cycle in the future, because these processes appear to play a very significant role in the carbon cycle of the present-day ocean and they are sensitive to climate change. Observations and models suggest that it may also be necessary to include the DMS cycle to predict future climate, though the effects are still highly uncertain. We have learned a tremendous amount about the distributions and biogeochemical impact of bacteria in the ocean in recent years, yet this improved understanding has not yet been incorporated into many of our models. All of these considerations lead us toward the development of increasingly complex models. However, recent quantitative model intercomparison studies suggest that continuing to add complexity and more functional groups to our ecosystem models may lead to decreases in predictive ability if the models are not properly constrained with available data. We also caution that capturing the present-day variability tells us little about how well a particular model can predict the future. If our goal is to develop models that can be used to predict how the oceans will respond to global warming, then we need to make more rigorous assessments of predictive skill using the available data. © 2006 Elsevier Ltd. All rights reserved.

Satellite remote sensing estimates of surface chlorophyll, temperature, wind speed, and sea ice cover are examined in the region of the Southern Ocean Iron Experiment (SOFeX). Our objectives are to place SOFeX into a regional context and highlight regional mesoscale spatial and monthly temporal variability. SOFeX fertilized two patches with iron, one south of the Antarctic Polar Front (PF) and one north of the PF but south of the Subantarctic Front (SAF). Satellite-observable phytoplankton blooms developed in both patches. The spring sea ice retreat near the south patch site was delayed in the 2001–2002 season, in turn delaying the naturally occurring, modest spring bloom in this region. Ambient surface chlorophyll concentrations for the area surrounding the southern patch during January 2002 are low (mean 0.26 mg/m3) compared with climatological January values (0.42 mg/m3). Regions east and west at similar latitudes exhibited higher mean chlorophyll concentrations (0.79 and 0.74 mg/m3, respectively). These modest phytoplankton blooms were likely stimulated by melting sea ice via changes in the light-mixing regime and release of iron and were smaller in magnitude than the iron-induced bloom within the SOFeX southern patch (>3 mg/m3). Iron inputs from melting ice may drive much of the natural spatial and temporal variability within the seasonal ice zone. Mean chlorophyll concentrations surrounding the SOFeX northern patch site were similar to climatological values during the SOFeX season. The northern patch was stretched into a long, thin filament along the southern boundary of the SAF, likely increasing the mixing/dilution rate with surrounding waters.

A global three-dimensional marine ecosystem model with several key phytoplankton functional groups, multiple limiting nutrients, explicit iron cycling, and a mineral ballast/organic matter parameterization is run within a global ocean circulation model. The coupled biogeochemistry/ecosystem/circulation (BEC) model reproduces known basin-scale patterns of primary and export production, biogenic silica production, calcification, chlorophyll, macronutrient and dissolved iron concentrations. The model captures observed high nitrate, low chlorophyll (HNLC) conditions in the Southern Ocean, subarctic and equatorial Pacific. Spatial distributions of nitrogen fixation are in general agreement with field data, with total N-fixation of 55 Tg N. Diazotrophs directly account for a small fraction of primary production (0.5%) but indirectly support 10% of primary production and 8% of sinking particulate organic carbon (POC) export. Diatoms disproportionately contribute to export of POC out of surface waters, but CaCO3 from the coccolithophores is the key driver of POC flux to the deep ocean in the model. An iron source from shallow ocean sediments is found critical in preventing iron limitation in shelf regions, most notably in the Arctic Ocean, but has a relatively localized impact. In contrast, global-scale primary production, export production, and nitrogen fixation are all sensitive to variations in atmospheric mineral dust inputs. The residence time for dissolved iron in the upper ocean is estimated to be a few years to a decade. Most of the iron utilized by phytoplankton is from subsurface sources supplied by mixing, entrainment, and ocean circulation. However, owing to the short residence time of iron in the upper ocean, this subsurface iron pool is critically dependent on continual replenishment from atmospheric dust deposition and, to a lesser extent, lateral transport from shelf regions.

Satellite ocean color data from the Sea Viewing Wide Field of View Sensor (SeaWiFS) are used to investigate phytoplankton bloom dynamics at the Antarctic Polar Front (PF). Satellite sea surface temperature (SST) data are used to map the location of the PF at weekly timescales. Elevated chlorophyll within the PF often appears as a narrow band that occupies only a portion of the SST gradient across the PF. Phytoplankton blooms within the PF occur most frequently during the month of December and are unevenly distributed within the Southern Ocean. Elevated chlorophyll concentrations at the PF are most frequently observed where the current is interacting with large topographic features. Mesoscale physical processes, including meander-induced upwelling and increased eddy mixing, where the PF encounters large topographic features likely leads to increased nutrient flux to surface waters in these regions. The highest mean chlorophyll values associated with the PF occur where the front comes into contact with relatively shallow waters along the North Scotia Ridge and at Kerguelen Plateau. Iron input from sedimentary sources likely plays an important role in these regions. Over seasonal timescales it appears likely that light-limitation prevents phytoplankton blooms at the PF during winter and spring months. PF blooms are observed most commonly during December when surface radiation peaks and mixed layer depths are rapidly shoaling. Even during December, when the light regime would seem to be favorable, PF blooms are largely restricted to regions where enhanced nutrient fluxes to surface waters due to frontal dynamics are likely. During late summer, nutrient limitation due to depletion of iron and possibly silicate largely prevent blooms at the PF. In the fall, deepening mixed layers would provide some relief from nutrient limitation but likely lead again to light-limitation of growth rates and the prevention of blooms. © 2002 Elsevier Science B.V. All rights reserved.

The role of iron in regulating the flux of carbon through the surface layer of the ocean has become increasingly apparent during the past 15 years. Before that time, the analytical challenges of measuring trace (parts per trillion) iron concentrations from iron ships using gear suspended on an iron wire precluded oceanographers from making accurate measurements. Laboratory experiments were invariably conducted with samples that were seriously contaminated with elevated iron concentrations. We now recognize, through greatly improved methodologies, that iron is a key regulator of phytoplankton primary production throughout the ocean. Small changes in iron concentration may produce large variations in the export of particulate organic carbon from the ocean's sunlit surface layer into deep-sea sediments. These variations in carbon export may occur over glacial/interglacial cycles at a scale sufficient to influence the flux of carbon dioxide from the atmosphere to the ocean. Such processes have been hypothesized to be an important driver of the changes in atmospheric carbon dioxide concentration that have been recorded in ice cores over the past 400,000 years.

A global marine ecosystem mixed-layer model is used to study iron cycling and nutrient-limitation patterns in surface waters of the world ocean. The ecosystem model has a small phytoplankton size class whose growth can be limited by N, P, Fe, and/or light, a diatom class which can also be Si-limited, and a diazotroph phytoplankton class whose growth rates can be limited by P, Fe, and/or light levels. The model also includes a parameterization of calcification by phytoplankton and is described in detail by Moore et al. (Deep-Sea Res. II, 2002). The model reproduces the observed high nitrate, low chlorophyll (HNLC) conditions in the Southern Ocean, subarctic Northeast Pacific, and equatorial Pacific, and realistic global patterns of primary production, biogenic silica production, nitrogen fixation, particulate organic carbon export, calcium carbonate export, and surface chlorophyll concentrations. Phytoplankton cellular Fe/C ratios and surface layer dissolved iron concentrations are also in general agreement with the limited field data. Primary production, community structure, and the sinking carbon flux are quite sensitive to large variations in the atmospheric iron source, particularly in the HNLC regions, supporting the Iron Hypothesis of Martin (Paleoceanography 5 (1990) 1-13). Nitrogen fixation is also strongly influenced by atmospheric iron deposition. Nitrogen limits phyuoplankton growth rates over less than half of the world ocean during summer months. Export of biogenic carbon is dominated by the sinking particulate flux, but detrainment and turbulent mixing account for 30% of global carbon export. Our results, in conjunction with other recent studies, suggest the familiar paradigm that nitrate inputs to the surface layer can be equated with particulate carbon export needs to be expanded to include multiple limiting nutrients and modes of export. © 2001 Elsevier Science Ltd. All rights reserved.

A new marine ecosystem model designed for the global domain is presented, and model output is compared with field data from nine different locations. Field data were collected as part of the international Joint Global Ocean Flux Study (JGOFS) program, and from historical time series stations. The field data include a wide variety of marine ecosystem types, including nitrogen- and iron-limited systems, and different physical environments from high latitudes to the mid-ocean gyres. Model output is generally in good agreement with field data from these diverse ecosystems. These results imply that the ecosystem model presented here can be reliably applied over the global domain. The model includes multiple potentially limiting nutrients that regulate phytoplankton growth rates. There are three phytoplankton classes, diatoms, diazotrophs, and a generic small phytoplankton class. Growth rates can be limited by available nitrogen, phosphorus, iron, and/or light levels. The diatoms can also be limited by silicon. The diazotrophs are capable of nitrogen fixation of N2 gas and cannot be nitrogen-limited. Calcification by phytoplankton is parameterized as a variable fraction of primary production by the small phytoplankton group. There is one zooplankton class that grazes the three phytoplankton groups and a large detrital pool. The large detrital pool sinks out of the mixed layer, while a smaller detrital pool, representing dissolved organic matter and very small particulates, does not sink. Remineralization of the detrital pools is parameterized with a temperature-dependent function. We explicitly model the dissolved iron cycle in marine surface waters including inputs of iron from subsurface sources and from atmospheric dust deposition. © 2001 Elsevier Science Ltd. All rights reserved.

Analysis of satellite ocean color, sea surface temperature, and sea ice cover data reveals consistent patterns between biological production, iron availability, and physical forcings in the Southern Ocean. The consistency of these patterns, in conjunction with information on physical conditions during the last glacial maximum (LGM), enables estimates of export production at the LGM. The LGM Southern Ocean experienced increased wind speeds, colder sea surface and atmospheric temperatures, increased deposition of atmospheric dust, and a greatly expanded winter sea ice cover. These variations had strong effects on Southern Ocean ecology and on air-sea fluxes of CO2. The seasonal ice zone (SIZ) was much larger at the LGM (30 million km2) than at present (19 million km2). The Antarctic Polar Front (PF) likely marked the northern boundary of this expanded SIZ throughout the Southern Ocean, as it does today in the Drake Passage region. A large northward shift in the position of the PF during glacial times is unlikely due to topographic constraints. North of the PF, the increased flux of aeolian dust during glacial times altered phytoplankton species composition and increased export production, and as a result this region was a stronger sink for atmospheric CO2 than in the modern ocean. South of the PF, interactions between the biota and sea ice strongly influence air-sea gas exchange over seasonal timescales. The combined influence of melting sea ice and increased aeolian dust flux (with its associated iron) increased both primary and export production by phytoplankton over daily-monthly timescales during austral spring/summer, resulting in a strong flux of CO2 into the ocean. Heavy ice cover would have minimized air-sea gas exchange over much of the rest of the year. Thus, an increased net flux of CO2 into the ocean is likely during glacial times, even in areas where annual primary production declined. We estimate that export production in the Southern Ocean as a whole was increased by 2.9-3.6 Gt C yr−1 at the LGM, relative to the modern era. Altered seasonal sea ice dynamics would further increase the net flux of CO2 into the ocean. Thus the Southern Ocean was a strong sink for atmospheric CO2 and contributed substantially to the lowering of atmospheric CO2 levels during the last ice age.

The United States Southern Ocean Joint Global Ocean Flux Study (JGOFS), also known as AESOPS (Antarctic Environment and Southern Ocean Process Study), focused on two distinct regions. The first was the Ross-Sea continental shelf, where a series of six cruises collected a variety of data from October 1996 through February 1998. The second area was the southwest Pacific sector of the Southern Ocean, spanning the Antarctic Circumpolar Current (ACC) at ∼ 170°W. Data were collected within this region during five cruises from September 1996 through March 1998, as well as during selected transits between New Zealand and the Ross Sea. The first results of these cruises are described in this issue. The Ross-Sea investigation extensively sampled the area along 76°30′S to elucidate the temporal patterns and processes that contribute to making this one of the Antarctic's most productive seas. Hydrographic distributions confirm that stratification is initiated early in October within the polynya, generating an environment that is favorable for phytoplankton growth. Significant spatial variations in mixed-layer depths, the timing of the onset of stratification, and the strength of the stratification existed throughout the growing season. Nutrient concentrations reflected phytoplankton uptake, and reached their seasonal minimal in early February. Chlorophyll concentrations were maximal in early January, whereas productivity was maximal in late November, which reflects the temporal uncoupling between growth and biomass accumulation in the region. Independent estimates of biogenic export suggest that majority of the flux occurred in late summer and was strongly uncoupled from phytoplankton growth. The ACC region exhibited seasonal changes that in some cases were greater than those observed in the Ross Sea.

Satellite ocean color data from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) were used to examine distributions of chlorophyll concentration within the Southern Ocean for the period October 1997 through September 1998. Over most of the Southern Ocean, mean chlorophyll concentrations remained quite low (< 0.3–0.4 mg m−3). Phytoplankton blooms where chlorophyll concentration exceeded 1.0 mg m−3 were observed in three general areas, which included coastal/shelf waters, areas associated with the seasonal sea ice retreat, and the vicinity of the major Southern Ocean fronts. These chlorophyll distribution patterns are consistent with an iron-limited system. Mean chlorophyll concentrations from SeaWiFS are compared with values from the coastal zone color scanner (CZCS). The SeaWiFS global chlorophyll algorithm works better than the CZCS in Southern Ocean waters. Primary production in the Southern Ocean was estimated with the vertically generalized production model of Behrenfeld and Falkowski [1997]. Annual primary production in the Southern Ocean (>30°S) was estimated to be 14.2 Gt C yr−1, with most production (∼80%) taking place at midlatitudes from 30° to 50°S. Primary production at latitudes >50°S was estimated to be 2.9 Gt C yr−1. This is considerably higher than previous estimates based on in situ data but less than some recent estimates based on CZCS data. Our estimated primary production is sufficient to account for the observed Southern Hemisphere seasonal cycle in atmospheric O2 concentrations.

SeaWiFS estimates of surface chlorophyll concentrations are reported for the region of the U.S. JGOFS study in the Southern Ocean (∼ 170°W, 60°S). Elevated chlorophyll was observed at the Southern Ocean fronts, near the edge of the seasonal ice sheet, and above the Pacific-Antarctic Ridge. The elevated chlorophyll levels associated with the Pacific-Antarctic Ridge are surprising since even the crest of the ridge is at depths > 2000 m. This elevated phytoplankton biomass is likely the result of mesoscale physical-biological interactions where the Antarctic Circumpolar Current (ACC) encounters the ridge. Four cruises surveyed this region between October 1997 and March 1998, as part of the U.S. JGOFS. Satellite-derived chlorophyll concentrations were compared with in situ extracted chlorophyll measurements from these cruises. There was good agreement (r² of 0.72, from a linear regression of shipboard vs. satellite chlorophyll), although SeaWiFS underestimated chlorophyll concentrations relative to the ship data.

The location of the Antarctic Polar Front (PF) was mapped over a 7-year period (1987–1993) within images of satellite-derived sea surface temperature. The mean path of the PF is strongly steered by the topographic features of the Southern Ocean. The topography places vorticity constraints on the dynamics of the PF that strongly affect spatial and temporal variability. Over the deep ocean basins the surface expression of the PF is weakened, and the PF meanders over a wide latitudinal range. Near large topographic features, width and temperature change across the front increase, and large-scale meandering is inhibited. Elevated mesoscale variability is seen within and downstream of these areas and may be the result of baroclinic instabilities initiated where the PF encounters large topographic features. The strong correlations between topography and PF dynamics can be understood in the context of the planetary potential vorticity (PPV orf/H) field. Mean PPV at the PF varies by more than a factor of 2 along its circumpolar path. However, at the mesoscale the PF remains within a relatively narrow range of PPV values around the local mean. Away from large topographic features, the PF returns to a preferred PPV value of ∼25 × 10−9 m−1 s−1 despite large latitudinal shifts. The mean paths of the surface and subsurface expressions of the PF are closely coupled over much of the Southern Ocean.

The path of the Antarctic Polar Front (PF) is mapped using satellite sea surface temperature data from the NOAA/NASA Pathfinder program. The mean path and variability of the PF are strongly influenced by bathymetry. Meandering intensity is weaker where the bathymetry is steeply sloped and increases in areas where the bottom is relatively flat. There is an inverse relationship between meandering intensity and both the width of the front and the change in temperature across it There is a persistent, large separation between the surface and subsurface expressions of the PF at Ewing Bank on the Falkland Plateau.

Rhizosolenia real abundance, distribution and chemical composition were studied on two cruises in the central North Pacific eyre in order to determine large-scale distribution patterns and contribution to upward nitrogen (N) flux. These macroscopic diatom mats are composed of multiple species of Rhizosolenia that exploit subsurface nitrate pools by vertically migrating below the euphotic zone. Although numerically dominated by the small-diameter species. R.fallax (73-95% of total numbers), mat biovolume was dominated by large-diameter (>50 μm diameter) Rhizostdenia spp. (85- 99% of total volume). Integrated mat abundance was substantially higher when mats accumulated at the surface during calm weather (