Abstract. Ocean ecosystems are increasingly stressed by human-induced changes of their physical, chemical and biological environment. Among these changes, warming, acidification, deoxygenation and changes in primary productivity by marine phytoplankton can be considered as four of the major stressors of open ocean ecosystems. Due to rising atmospheric CO 2 in the coming decades, these changes will be amplified. Here, we use the most recent simulations performed in the framework of the Coupled Model Intercomparison Project 5 to assess how these stressors may evolve over the course of the 21st century. The 10 Earth system models used here project similar trends in ocean warming, acidification, deoxygenation and reduced primary productivity for each of the IPCC's representative concentration pathways (RCPs) over the 21st century. For the "businessas-usual" scenario RCP8.5, the model-mean changes in the 2090s (compared to the 1990s) for sea surface temperature, sea surface pH, global O 2 content and integrated primary productivity amount to +2.73 (±0.72) • C, −0.33 (±0.003) pH unit, −3.45 (±0.44) % and −8.6 (±7.9) %, respectively. For the high mitigation scenario RCP2.6, corresponding changes are +0.71 (±0.45) • C, −0.07 (±0.001) pH unit, −1.81 (±0.31) % and −2.0 (±4.1) %, respectively, illustrating the effectiveness of extreme mitigation strategies. Although these stressors operate globally, they display distinct regional patterns and thus do not change coincidentally. Large decreases in O 2 and in pH are simulated in global ocean intermediate and mode waters, whereas large reductions in primary production are simulated in the tropics and in the North Atlantic. Although temperature and pH projections are robust across models, the same does not hold for projections of subsurface O 2 concentrations in the tropics and global and regional changes in net primary productivity. These high uncertainties in projections of primary productivity and subsurface oxygen prompt us to continue intermodel comparisons to understand these model differences, while calling for caution when using the CMIP5 models to force regional impact models.

The ocean's biological pump strips nutrients out of the surface waters and exports them into the thermocline and deep waters. If there were no return path of nutrients from deep waters, the biological pump would eventually deplete the surface waters and thermocline of nutrients; surface biological productivity would plummet. Here we make use of the combined distributions of silicic acid and nitrate to trace the main nutrient return path from deep waters by upwelling in the Southern Ocean and subsequent entrainment into subantarctic mode water. We show that the subantarctic mode water, which spreads throughout the entire Southern Hemisphere and North Atlantic Ocean, is the main source of nutrients for the thermocline. We also find that an additional return path exists in the northwest corner of the Pacific Ocean, where enhanced vertical mixing, perhaps driven by tides, brings abyssal nutrients to the surface and supplies them to the thermocline of the North Pacific. Our analysis has important implications for our understanding of large-scale controls on the nature and magnitude of low-latitude biological productivity and its sensitivity to climate change.

The formulation and simulation characteristics of two new global coupled climate models developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) are described. The models were designed to simulate atmospheric and oceanic climate and variability from the diurnal time scale through multicentury climate change, given our computational constraints. In particular, an important goal was to use the same model for both experimental seasonal to interannual forecasting and the study of multicentury global climate change, and this goal has been achieved.Two versions of the coupled model are described, called CM2.0 and CM2.1. The versions differ primarily in the dynamical core used in the atmospheric component, along with the cloud tuning and some details of the land and ocean components. For both coupled models, the resolution of the land and atmospheric components is 2°latitude ϫ 2.5°longitude; the atmospheric model has 24 vertical levels. The ocean resolution is 1°in latitude and longitude, with meridional resolution equatorward of 30°becoming progressively finer, such that the meridional resolution is 1/3°at the equator. There are 50 vertical levels in the ocean, with 22 evenly spaced levels within the top 220 m. The ocean component has poles over North America and Eurasia to avoid polar filtering. Neither coupled model employs flux adjustments.The control simulations have stable, realistic climates when integrated over multiple centuries. Both models have simulations of ENSO that are substantially improved relative to previous GFDL coupled models. The CM2.0 model has been further evaluated as an ENSO forecast model and has good skill (CM2.1 has not been evaluated as an ENSO forecast model). Generally reduced temperature and salinity biases exist in CM2.1 relative to CM2.0. These reductions are associated with 1) improved simulations of surface wind stress in CM2.1 and associated changes in oceanic gyre circulations; 2) changes in cloud tuning and the land model, both of which act to increase the net surface shortwave radiation in CM2.1, thereby reducing an overall cold bias present in CM2.0; and 3) a reduction of ocean lateral viscosity in the extratropics in CM2.1, which reduces sea ice biases in the North Atlantic. Both models have been used to conduct a suite of climate change simulations for the 2007 Intergovernmental Panel on Climate Change (IPCC) assessment report and are able to simulate the main features of the observed warming of the twentieth century. The climate sensitivities of the CM2.0 and CM2.1 models are 2.9 and 3.4 K, respectively. These sensitivities are defined by coupling the atmospheric components of CM2.0 and CM2.1 to a slab ocean model and allowing the model to come into equilibrium with a doubling of atmospheric CO 2 . The output from a suite of integrations conducted with these models is freely available online (see http://nomads.gfdl.noaa.gov/).

Ocean ecosystems are increasingly stressed by human-induced changes of their physical, chemical and biological environment. Among these changes, warming, acidification, deoxygenation and changes in primary productivity by marine phytoplankton can be considered as four of the major stressors of open ocean ecosystems. Due to rising atmospheric CO₂ in the coming decades, these changes will be amplified. Here, we use the most recent simulations performed in the framework of the Coupled Model Intercomparison Project 5 to assess how these stressors may evolve over the course of the 21st century. The 10 Earth System Models used here project similar trends in ocean warming, acidification, deoxygenation and reduced primary productivity for each of the IPCC's representative concentration parthways (RCP) over the 21st century. For the "business-as-usual" scenario RCP8.5, the model-mean changes in 2090s (compared to 1990s) for sea surface temperature, sea surface pH, global O₂ content and integrated primary productivity amount to +2.73 °C, −0.33 pH unit, −3.45% and −8.6%, respectively. For the high mitigation scenario RCP2.6, corresponding changes are +0.71 °C, −0.07 pH unit, −1.81% and −2.0% respectively, illustrating the effectiveness of extreme mitigation strategies. Although these stressors operate globally, they display distinct regional patterns. Large decreases in O₂ and in pH are simulated in global ocean intermediate and mode waters, whereas large reductions in primary production are simulated in the tropics and in the North Atlantic. Although temperature and pH projections are robust across models, the same does not hold for projections of sub-surface O₂ concentrations in the tropics and global and regional changes in net primary productivity

[1] We present a new synthesis of the oceanic cycles of organic carbon, silicon, and calcium carbonate. Our calculations are based on a series of algorithms starting with satellite-based primary production and continuing with conversion of primary production to sinking particle flux, penetration of particle flux to the deep sea, and accumulation in sediments. Regional and global budgets from this synthesis highlight the potential importance of shelves and near-shelf regions for carbon burial. While a high degree of uncertainty remains, this analysis suggests that shelves, less than 50 m water depths accounting for 2% of the total ocean area, may account for 48% of the global flux of organic carbon to the seafloor. Our estimates of organic carbon and nitrogen flux are in generally good agreement with previous work while our estimates for CaCO 3 and SiO 2 fluxes are lower than recent work. Interannual variability in particle export fluxes is found to be relatively small compared to intra-annual variability over large domains with the single exception of the dominating role of El Niño-Southern Oscillation variability in the central tropical Pacific. Comparison with available sediment-based syntheses of benthic remineralization and burial support the recent theory of mineral protection of organic carbon flux through the deep ocean, pointing to lithogenic material as an important carrier phase of organic carbon to the deep seafloor. This work suggests that models which exclude the role of lithogenic material would underestimate the penetration of POC to the deep seafloor by approximately 16-51% globally, and by a much larger fraction in areas with low productivity. Interestingly, atmospheric dust can only account for 31% of the total lithogenic flux and 42% of the lithogenically associated POC flux, implying that a majority of this material is riverine or directly erosional in origin.Citation: Dunne, J. P., J. L. Sarmiento, and A. Gnanadesikan (2007), A synthesis of global particle export from the surface ocean and cycling through the ocean interior and on the seafloor, Global Biogeochem. Cycles, 21, GB4006,

Nitrogen fixation is crucial for maintaining biological productivity in the oceans, because it replaces the biologically available nitrogen that is lost through denitrification. But, owing to its temporal and spatial variability, the global distribution of marine nitrogen fixation is difficult to determine from direct shipboard measurements. This uncertainty limits our understanding of the factors that influence nitrogen fixation, which may include iron, nitrogen-to-phosphorus ratios, and physical conditions such as temperature. Here we determine nitrogen fixation rates in the world's oceans through their impact on nitrate and phosphate concentrations in surface waters, using an ocean circulation model. Our results indicate that nitrogen fixation rates are highest in the Pacific Ocean, where water column denitrification rates are high but the rate of atmospheric iron deposition is low. We conclude that oceanic nitrogen fixation is closely tied to the generation of nitrogen-deficient waters in denitrification zones, supporting the view that nitrogen fixation stabilizes the oceanic inventory of fixed nitrogen over time.

The physical climate formulation and simulation characteristics of two new global coupled carbon–climate Earth System Models, ESM2M and ESM2G, are described. These models demonstrate similar climate fidelity as the Geophysical Fluid Dynamics Laboratory’s previous Climate Model version 2.1 (CM2.1) while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4p1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in El Niño–Southern Oscillation being overly strong in ESM2M and overly weak in ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to total heat content variability given its lack of long-term drift, gyre circulation, and ventilation in the North Pacific, tropical Atlantic, and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to surface circulation given its superior surface temperature, salinity, and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. The overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon–climate models.

[1] We present new empirical and mechanistic models for predicting the export of organic carbon out of the surface ocean by sinking particles. To calibrate these models, we have compiled a synthesis of field observations related to ecosystem size structure, primary production and particle export from around the globe. The empirical model captures 61% of the observed variance in the ratio of particle export to primary production (the pe ratio) using sea-surface temperature and chlorophyll concentrations (or primary productivity) as predictor variables. To describe the mechanisms responsible for pe-ratio variability, we present size-based formulations of phytoplankton grazing and sinking particle export, combining them into an alternative, mechanistic model. The formulation of grazing dynamics, using simple power laws as closure terms for small and large phytoplankton, reproduces 74% of the observed variability in phytoplankton community composition wherein large phytoplankton augment small ones as production increases. The formulation for sinking particle export partitions a temperature-dependent fraction of small and large phytoplankton grazing into sinking detritus. The mechanistic model also captures 61% of the observed variance in pe ratio, with large phytoplankton in high biomass and relatively cold regions leading to more efficient export. In this model, variability in primary productivity results in a biomass-modulated switch between small and large phytoplankton pathways.Citation: Dunne, J. P., R. A. Armstrong, A. Gnanadesikan, and J. L. Sarmiento (2005), Empirical and mechanistic models for the particle export ratio, Global Biogeochem. Cycles, 19, GB4026,