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

The Geophysical Fluid Dynamics Laboratory (GFDL) has developed a coupled general circulation model (CM3) for the atmosphere, oceans, land, and sea ice. The goal of CM3 is to address emerging issues in climate change, including aerosol-cloud interactions, chemistry-climate interactions, and coupling between the troposphere and stratosphere. The model is also designed to serve as the physical system component of earth system models and models for decadal prediction in the near-term future-for example, through improved simulations in tropical land precipitation relative to earlier-generation GFDL models. This paper describes the dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component (AM3) of this model. Relative to GFDL AM2, AM3 includes new treatments of deep and shallow cumulus convection, cloud droplet activation by aerosols, subgrid variability of stratiform vertical velocities for droplet activation, and atmospheric chemistry driven by emissions with advective, convective, and turbulent transport. AM3 employs a cubed-sphere implementation of a finite-volume dynamical core and is coupled to LM3, a new land model with ecosystem dynamics and hydrology. Its horizontal resolution is approximately 200 km, and its vertical resolution ranges approximately from 70 m near the earth's surface to 1 to 1.5 km near the tropopause and 3 to 4 km in much of the stratosphere. Most basic circulation features in AM3 are simulated as realistically, or more so, as in AM2. In particular, dry biases have been reduced over South America. In coupled mode, the simulation of Arctic sea ice concentration has improved. AM3 aerosol optical depths, scattering properties, and surface clear-sky downward shortwave radiation are more realistic than in AM2. The simulation of marine stratocumulus decks remains problematic, as in AM2. The most intense 0.2% of precipitation rates occur less frequently in AM3 than observed. The last two decades of the twentieth century warm in CM3 by 0.328C relative to 1881-1920. The Climate Research Unit (CRU) and Goddard Institute for Space Studies analyses of observations show warming of 0.568 and 0.528C, respectively, over this period. CM3 includes anthropogenic cooling by aerosol-cloud interactions, and its warming by the late twentieth century is somewhat less realistic than in CM2.1, which warmed 0.668C but did not include aerosol-cloud interactions. The improved simulation of the direct aerosol effect (apparent in surface clear-sky downward radiation) in CM3 evidently acts in concert with its simulation of cloud-aerosol interactions to limit greenhouse gas warming.

[1] Observations show that strong equatorial volcanic eruptions have been followed by a pronounced positive phase of the Arctic Oscillation (AO) for one or two Northern Hemisphere winters. It has been previously assumed that this effect is forced by strengthening of the equator-to-pole temperature gradient in the lower stratosphere, caused by aerosol radiative heating in the tropics. To understand atmospheric processes that cause the AO response, we studied the impact of the 1991 Mount Pinatubo eruption, which produced the largest global volcanic aerosol cloud in the twentieth century. A series of control and perturbation experiments were conducted with the GFDL SKYHI general circulation model to examine the evolution of the circulation in the 2 years following the Pinatubo eruption. In one set of perturbation experiments, the full radiative effects of the observed Pinatubo aerosol cloud were included, while in another only the effects of the aerosols in reducing the solar flux in the troposphere were included, and the aerosol heating effects in the stratosphere were suppressed. A third set of perturbation experiments imposed the stratospheric ozone losses observed in the post-Pinatubo period. We conducted ensembles of four to eight realizations for each case. Forced by aerosols, SKYHI produces a statistically significant positive phase of the AO in winter, as observed. Ozone depletion causes a positive phase of the AO in late winter and early spring by cooling the lower stratosphere in high latitudes, strengthening the polar night jet, and delaying the final warming. A positive phase of the AO was also produced in the experiment with only the tropospheric effect of aerosols, showing that aerosol heating in the lower tropical stratosphere is not necessary to force positive AO response, as was previously assumed. Aerosol-induced tropospheric cooling in the subtropics decreases the meridional temperature gradient in the winter troposphere between 30°N and 60°N. The corresponding reduction of mean zonal energy and amplitudes of planetary waves in the troposphere decreases wave activity flux into the lower stratosphere. The resulting strengthening of the polar vortex forces a positive phase of the AO. We suggest that this mechanism can also contribute to the observed long-term AO trend being caused by greenhouse gas increases because they also weaken the tropospheric meridional temperature gradient due to polar amplification of warming.