No abstract

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/).

[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,

[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,

A simple theory for the large-scale oceanic circulation is developed, relating pycnocline depth, Northern Hemisphere sinking, and low-latitude upwelling to pycnocline diffusivity and Southern Ocean winds and eddies. The results show that Southern Ocean processes help maintain the global ocean structure and that pycnocline diffusion controls low-latitude upwelling.

The Modeling Eddies in the Southern Ocean (MESO) project uses numerical sensitivity studies to examine the role played by Southern Ocean winds and eddies in determining the density structure of the global ocean and the magnitude and structure of the global overturning circulation. A hemispheric isopycnal-coordinate ocean model (which avoids numerical diapycnal diffusion) with realistic geometry is run with idealized forcing at a range of resolutions from coarse (2°) to eddy-permitting ( 1 ⁄6°). A comparison of coarse resolutions with fine resolutions indicates that explicit eddies affect both the structure of the overturning and the response of the overturning to wind stress changes. While the presence of resolved eddies does not greatly affect the prevailing qualitative picture of the ocean circulation, it alters the overturning cells involving the Southern Ocean transformation of dense deep waters and light waters of subtropical origin into intermediate waters. With resolved eddies, the surface-to-intermediate water cell extends farther southward by hundreds of kilometers and the deep-to-intermediate cell draws on comparatively lighter deep waters. The overturning response to changes in the winds is also sensitive to the presence of eddies. In noneddying simulations, changing the Ekman transport produces comparable changes in the overturning, much of it involving transformation of deep waters and resembling the mean circulation. In the eddypermitting simulations, a significant fraction of the Ekman transport changes are compensated by eddyinduced transport drawing from lighter waters than does the mean overturning. This significant difference calls into question the ability of coarse-resolution ocean models to accurately capture the impact of changes in the Southern Ocean on the global ocean circulation.

Modelling studies have demonstrated that the nutrient and carbon cycles in the Southern Ocean play a central role in setting the air-sea balance of CO(2) and global biological production. Box model studies first pointed out that an increase in nutrient utilization in the high latitudes results in a strong decrease in the atmospheric carbon dioxide partial pressure (pCO2). This early research led to two important ideas: high latitude regions are more important in determining atmospheric pCO2 than low latitudes, despite their much smaller area, and nutrient utilization and atmospheric pCO2 are tightly linked. Subsequent general circulation model simulations show that the Southern Ocean is the most important high latitude region in controlling pre-industrial atmospheric CO(2) because it serves as a lid to a larger volume of the deep ocean. Other studies point out the crucial role of the Southern Ocean in the uptake and storage of anthropogenic carbon dioxide and in controlling global biological production. Here we probe the system to determine whether certain regions of the Southern Ocean are more critical than others for air-sea CO(2) balance and the biological export production, by increasing surface nutrient drawdown in an ocean general circulation model. We demonstrate that atmospheric CO(2) and global biological export production are controlled by different regions of the Southern Ocean. The air-sea balance of carbon dioxide is controlled mainly by the biological pump and circulation in the Antarctic deep-water formation region, whereas global export production is controlled mainly by the biological pump and circulation in the Subantarctic intermediate and mode water formation region. The existence of this biogeochemical divide separating the Antarctic from the Subantarctic suggests that it may be possible for climate change or human intervention to modify one of these without greatly altering the other.

Abstract.Laboratory and field studies have revealed that iron has multiple roles in phytoplankton physiology, with particular importance for light-harvesting cellular machinery. However, although iron-limitation is explicitly included in numerous biogeochemical/ecosystem models, its implementation varies, and its effect on the efficiency of light harvesting is often ignored. Given the complexity of the ocean environment, it is difficult to predict the consequences of applying different iron limitation schemes. Here we explore the interaction of iron and nutrient cycles in an ocean general circulation model using a new, streamlined model of ocean biogeochemistry. Building on previously published parameterizations of photoadaptation and export production, the Biogeochemistry with Light Iron Nutrients and Gasses (BLING) model is constructed with only four explicit tracers but including macronutrient and micronutrient limitation, light limitation, and an implicit treatment of community structure. The structural simplicity of this computationallyinexpensive model allows us to clearly isolate the global effect that iron availability has on maximum light-saturated photosynthesis rates vs. the effect iron has on photosynthetic efficiency. We find that the effect on light-saturated photosynthesis rates is dominant, negating the importance of photosynthetic efficiency in most regions, especially the cold waters of the Southern Ocean. The primary exceptions to this occur in iron-rich regions of the Northern Hemisphere, where high light-saturated photosynthesis rates allow photosynthetic efficiency to play a more important role. In other words, the ability to efficiently harvest photons has little effect in regions where light-saturated growth rates are low. Additionally, we speculate that the phytoplankton cells domCorrespondence to: E. D. Galbraith (eric.galbraith@mcgill.ca) inating iron-limited regions tend to have relatively high photosynthetic efficiency, due to reduced packaging effects. If this speculation is correct, it would imply that natural communities of iron-stressed phytoplankton may tend to harvest photons more efficiently than would be inferred from ironlimitation experiments with other phytoplankton. We suggest that iron limitation of photosynthetic efficiency has a relatively small impact on global biogeochemistry, though it is expected to impact the seasonal cycle of plankton as well as the vertical structure of primary production.