Anand Gnanadesikan

GFDL's CM2 global coupled climate models. Part I: Formulation and simulation characteristics

Abstract The formulation and simulation characteristics of two new global coupled climate models developed at NOAAs 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, wi...

The Role of Eddies in Determining the Structure and Response of the Wind-Driven Southern Hemisphere Overturning: Results from the Modeling Eddies in the Southern Ocean (MESO) Project

Abstract 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-t...

Evaluation of ocean carbon cycle models with data-based metrics

New radiocarbon and chlorofluorocarbon-11 data from the World Ocean Circulation Experiment are used to assess a suite of 19 ocean carbon cycle models. We use the distributions and inventories of these tracers as quantitative metrics of model skill and find that only about a quarter of the suite is consistent with the new data-based metrics. This should serve as a warning bell to the larger community that not all is well with current generation of ocean carbon cycle models. At the same time, this highlights the danger in simply using the available models to represent the state-of-the-art modeling without considering the credibility of each model.

GFDL's CM2 Global Coupled Climate Models. Part II: The Baseline Ocean Simulation

The current generation of coupled climate models run at the Geophysical Fluid Dynamics Laboratory (GFDL) as part of the Climate Change Science Program contains ocean components that differ in almost every respect from those contained in previous generations of GFDL climate models. This paper summarizes the new physical features of the models and examines the simulations that they produce. Of the two new coupled climate model versions 2.1 (CM2.1) and 2.0 (CM2.0), the CM2.1 model represents a major improvement over CM2.0 in most of the major oceanic features examined, with strikingly lower drifts in hydrographic fields such as temperature and salinity, more realistic ventilation of the deep ocean, and currents that are closer to their observed values. Regional analysis of the differences between the models highlights the importance of wind stress in determining the circulation, particularly in the Southern Ocean. At present, major errors in both models are associated with Northern Hemisphere Mode Waters and outflows from overflows, particularly the Mediterranean Sea and Red Sea.

The Southern Ocean biogeochemical divide

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 CO2 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 CO2 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 CO2 balance and the biological export production, by increasing surface nutrient drawdown in an ocean general circulation model. We demonstrate that atmospheric CO2 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.

Evaluating global ocean carbon models: The importance of realistic physics

A suite of standard ocean hydrographic and circulation metrics are applied to the equilibrium physical solutions from 13 global carbon models participating in phase 2 of the Ocean Carbon-cycle Model Intercomparison Project (OCMIP-2). Model-data comparisons are presented for sea surface temperature and salinity, seasonal mixed layer depth, meridional heat and freshwater transport, 3-D hydrographic fields, and meridional overturning. Considerable variation exists among the OCMIP-2 simulations, with some of the solutions falling noticeably outside available observational constraints. For some cases, model-model and model-data differences can be related to variations in surface forcing, subgrid-scale parameterizations, and model architecture. These errors in the physical metrics point to significant problems in the underlying model representations of ocean transport and dynamics, problems that directly affect the OCMIP predicted ocean tracer and carbon cycle variables (e.g., air-sea CO2 flux, chlorofluorocarbon and anthropogenic CO2 uptake, and export production). A substantial fraction of the large model-model ranges in OCMIP-2 biogeochemical fields (±25–40%) represents the propagation of known errors in model physics. Therefore the model-model spread likely overstates the uncertainty in our current understanding of the ocean carbon system, particularly for transport-dominated fields such as the historical uptake of anthropogenic CO2. A full error assessment, however, would need to account for additional sources of uncertainty such as more complex biological-chemical-physical interactions, biases arising from poorly resolved or neglected physical processes, and climate change.

Oceanic ventilation and biogeochemical cycling: Understanding the physical mechanisms that produce realistic distributions of tracers and productivity

[1] Differing models of the ocean circulation support different rates of ventilation, which in turn produce different distributions of radiocarbon, oxygen, and export production. We examine these fields within a suite of general circulation models run to examine the sensitivity of the circulation to the parameterization of subgridscale mixing and surface forcing. We find that different models can explain relatively high fractions of the spatial variance in some fields such as radiocarbon, and that newer estimates of the rate of biological cycling are in better agreement with the models than previously published estimates. We consider how different models achieve such agreement and show that they can accomplish this in different ways. For example, models with high vertical diffusion move young surface waters into the Southern Ocean, while models with high winds move more young North Atlantic water into this region. The dependence on parameter values is not simple. Changes in the vertical diffusion coefficient, for example, can produce major changes in advective fluxes. In the coarse-resolution models studied here, lateral diffusion plays a major role in the tracer budget of the deep ocean, a somewhat worrisome fact as it is poorly constrained both observationally and theoretically. INDEX TERMS: 4275 Oceanography: General: Remote sensing and electromagnetic processes (0689); 4532 Oceanography: Physical: General circulation; 4568 Oceanography: Physical: Turbulence, diffusion, and mixing processes; 4845 Oceanography: Biological and Chemical: Nutrients and nutrient cycling; KEYWORDS: biogeochemical cycles, particle export, vertical exchange

The Southern Hemisphere Westerlies in a Warming World: Propping Open the Door to the Deep Ocean

Abstract A coupled climate model with poleward-intensified westerly winds simulates significantly higher storage of heat and anthropogenic carbon dioxide by the Southern Ocean in the future when compared with the storage in a model with initially weaker, equatorward-biased westerlies. This difference results from the larger outcrop area of the dense waters around Antarctica and more vigorous divergence, which remains robust even as rising atmospheric greenhouse gas levels induce warming that reduces the density of surface waters in the Southern Ocean. These results imply that the impact of warming on the stratification of the global ocean may be reduced by the poleward intensification of the westerlies, allowing the ocean to remove additional heat and anthropogenic carbon dioxide from the atmosphere.

The GFDL CM3 Coupled Climate Model: Characteristics of the Ocean and Sea Ice Simulations

AbstractThis paper documents time mean simulation characteristics from the ocean and sea ice components in a new coupled climate model developed at the NOAA Geophysical Fluid Dynamics Laboratory (GFDL). The GFDL Climate Model version 3 (CM3) is formulated with effectively the same ocean and sea ice components as the earlier CM2.1 yet with extensive developments made to the atmosphere and land model components. Both CM2.1 and CM3 show stable mean climate indices, such as large-scale circulation and sea surface temperatures (SSTs). There are notable improvements in the CM3 climate simulation relative to CM2.1, including a modified SST bias pattern and reduced biases in the Arctic sea ice cover. The authors anticipate SST differences between CM2.1 and CM3 in lower latitudes through analysis of the atmospheric fluxes at the ocean surface in corresponding Atmospheric Model Intercomparison Project (AMIP) simulations. In contrast, SST changes in the high latitudes are dominated by ocean and sea ice effects absen...

Transient Response in a Z-Level Ocean Model That Resolves Topography with Partial Cells

Ocean simulations are in part determined by topographic waves with speeds and spatial scales dependent on bottom slope. By their very nature, discrete z-level ocean models have problems accurately representing bottom topography when slopes are less than the grid cell aspect ratio Dz/Dx. In such regions, the dispersion relation for topographic waves is inaccurate. However, bottom topography can be accurately represented in discrete zlevel models by allowing bottom-most grid cells to be partially filled with land. Consequently, gently sloping bottom topography is resolved on the scale of horizontal grid resolution and the dispersion relation for topographic waves is accurately approximated. In contrast to the standard approach using full cells, partial cells imply that all grid points within a vertical level are not necessarily at the same depth and problems arise with pressure gradient errors and the spurious diapycnal diffusion. However, both problems have been effectively dealt with. Differences in flow fields between simulations with full cells and partial cells can be significant, and simulations with partial cells are more robust than with full cells. Partial cells provide a superior representation of topographic waves when compared to the standard method employing full cells.