Geophysical Fluid Dynamics Laboratory: Model GFDL CDG1 (R30 L14) 1992

Geophysical Fluid Dynamics Laboratory: Model GFDL CGG1 (R30 L14) 1992

AMIP Representative(s)

Dr. Richard Wetherald, Geophysical Fluid Dynamics Laboratory/NOAA, Princeton University, P.O. Box 308, Princeton, New Jersey 08540; Phone: +1-609-258-6515; Fax: +1-609-987-5063; e-mail: rw@GFDL.GOV; World Wide Web URL:

Model Designation

GFDL CDG1 (R30 L14) 1992

Model Lineage

One of several versions of global spectral models in use at the Geophysical Fluid Dynamics Laboratory, the GFDL model is applied primarily to climate studies. In its treatment of atmospheric dynamics, the model is similar to the GFDL Dynamical Extended-Range Forecasting (DERF) model (cf. Gordon and Stern 1982 [1] and Miyakoda and Sirutis 1986 [2]), but it displays a number of differences in numerical features (e.g., use of rhomboidal rather than triangular spectral truncation and a somewhat coarser vertical resolution), as well as in some model parameterizations (e.g., radiation and vertical diffusion).

Model Documentation

Key documentation of dynamical features is given by Gordon and Stern (1982) and Manabe and Hahn (1981) [3]. Physical parameterizations are described by Manabe (1969) [4], Holloway and Manabe (1971) [5], Manabe and Holloway (1975) [6], Manabe et al. (1965) [7], Broccoli and Manabe (1992) [8], Wetherald and Manabe (1988) [9], and Wetherald et al. (1991) [10].

Numerical/Computational Properties

Horizontal Representation

Spectral (spherical harmonic basis functions) with transformation to a Gaussian grid for calculation of nonlinear quantities and some physics.

Horizontal Resolution

Spectral rhomboidal 30 (R30), roughly equivalent to 2.25 x 3.75 degrees latitude-longitude.

Vertical Domain

Surface to about 15 hPa. For a surface pressure of 1000 hPa, the lowest atmospheric level is at about 997 hPa.

Vertical Representation

Finite-difference sigma coordinates.

Vertical Resolution

There are fourteen unevenly spaced sigma levels. For a surface pressure of 1000 hPa, four levels are below 800 hPa and four levels are above 200 hPa.

Computer/Operating System

The AMIP simulation was run on a Cray Y/MP computer using a single processor in a UNICOS operating environment.

Computational Performance

For the AMIP experiment, about 1 minute of Cray Y/MP computation time per simulated day.


For the AMIP simulation, the model atmosphere, soil moisture, and snow cover/depth are initialized for 1 January 1979 from a previous 10-year integration of the model with climatological sea surface temperatures (cf. Alexander and Mobley 1976) [11].

Time Integration Scheme(s)

A leapfrog semi-implicit scheme similar to that of Bourke (1974) [12] with an Asselin (1972) [13] frequency filter is used for time integration. The time step is 18 minutes for dynamics and physics, except for full calculations of all radiative fluxes, which are done once every 24 hours.


Orography is smoothed (see Orography). Negative moisture values (that arise because of spectral truncation) are filled by borrowing from nearest neighbors in the vertical, and horizontally in the east-west direction only.

Sampling Frequency

For the AMIP simulation, the model history is written once every 24 hours.

Dynamical/Physical Properties

Atmospheric Dynamics

Primitive-equation dynamics are expressed in terms of vorticity, divergence, surface pressure, specific humidity, and temperature (with a linearized correction for virtual temperature in diagnostic quantities, where applicable).


  • Linear fourth-order (del^4) horizontal diffusion is applied on constant pressure surfaces to vorticity, divergence, temperature, and specific humidity.

  • Second-order vertical diffusion with coefficients derived from mixing-length considerations is applied to momentum, heat, and moisture on the first 7 vertical levels above the surface (up to a height of about 5 km). The vertical diffusion is not stability-dependent.

Gravity-wave Drag

The parameterization of orographic gravity-wave drag follows linear theory, as formulated by Y. Hayashi (cf. Broccoli and Manabe 1992) [8]. The drag is given by the vertical divergence of the wave stress, which near the surface is equal to the product of the subgrid-scale orographic variance and a representative mountain wavenumber with a mass-weighted average of the atmospheric density, the Brunt-Vaisalla frequency, and the wind in the first three layers above the surface (centered on sigma = 0.9353, 0.97865, and 0.99665). The wave stress is assumed to be zero above a critical sigma level where the projection of the wind on the surface wind vector vanishes. Below this, the stress increases linearly to its surface value as the sigma-distance from this critical level increases.

Solar Constant/Cycles

The solar constant is the AMIP-prescribed value of 1365 W/(m^2). A seasonal, but not a diurnal cycle in solar forcing, is simulated.


The carbon dioxide concentration is the AMIP-prescribed value of 345 ppm. Zonally averaged seasonal mean ozone distributions are specified as a function of height after data of Hering and Borden (1965) [14]. These data are linearly interpolated for intermediate times. Radiative effects of water vapor, but not those of aerosols, also are included (see Radiation).


  • Shortwave Rayleigh scattering, and absorption in ultraviolet (wavelengths less than 0.35 micron) and visible (wavelengths between 0.5 to 0.7 micron) spectral bands by ozone, and in the near-infrared (wavelengths 0.7 to 4.0 microns) by water vapor are treated by a method similar to that of Lacis and Hansen (1974) [15]. Pressure corrections and multiple reflections between clouds and the surface are treated, but radiative effects of aerosols are not included.

  • Longwave radiation follows the method of Rodgers and Walshaw (1966) [16], as modified by Stone and Manabe (1968) [17]. Absorption by water vapor, carbon dioxide, and ozone is included. The 6.3-micron band, the rotation band, and the continuum of water vapor are subdivided into 19 subintervals, two of which contain the 15-micron carbon dioxide and 9.6-micron ozone absorption bands, with transmissions of the latter two gases multiplied together in overlapping bands. A random model is used to represent the absorptivity for each subinterval, and the Curtis-Godson approximation is used to estimate the effective pressure for absorption. The temperature dependence of line intensity is also incorporated. Carbon dioxide absorptivity is obtained from data of Burch et al. (1961) [18], with its temperature dependence following the model of Sasamori (1959) [19]. The ozone absorptivity is obtained from data of Walshaw (1957) [20].

  • Cloud-radiative interactions are treated as described by Wetherald and Manabe (1988) [9] and Wetherald et al. (1991) [10]. Cloud optical properties (absorptivity/ reflectivity in the ultraviolet-visible and near-infrared intervals, and emissivity in the longwave) depend on cloud height (high, middle, and low) and thickness. The values of shortwave properties are assigned following Rodgers (1967a) [21]; in the longwave, the emissivity of thin high clouds (above 10.5 km) is prescribed as 0.6 after Kondratiev (1972) [22], while all other clouds are treated as blackbodies (emissivity = 1.0). No partial cloudiness is accounted for in each grid box; clouds, therefore, are treated as fully overlapped in the vertical. See also Cloud Formation.


Penetrative convection is simulated via convective adjustment processes. If the lapse rate exceeds dry adiabatic, a convective adjustment restores the lapse rate to dry adiabatic, with conservation of dry static energy in the vertical. A moist convective adjustment scheme after Manabe et al. (1965) [7] also operates when the lapse rate exceeds moist adiabatic and the air is supersaturated. (For supersaturated stable layers, nonconvective large-scale condensation takes place--see Precipitation). In moist convective layers it is assumed that the intensity of convection is strong enough to eliminate the vertical gradient of potential temperature instantaneously, while conserving total moist static energy. It is further assumed that the relative humidity in the layer is maintained at 100 percent, owing to the vertical mixing of moisture, condensation, and evaporation from water droplets. Shallow convection is not explicitly simulated.

Cloud Formation

Cloud forms when the relative humidity of a vertical layer exceeds a height-dependent threshold, whether this results from large-scale condensation or moist convective adjustment (see Convection). (The threshold relative humidity varies linearly between 100 percent at the bottom of the model atmosphere to 90 percent at its top.) It is assumed that no condensed water is retained as cloud liquid water, but immediately precipitates (see Precipitation). Cloud type is defined according to height: high cloud forms above 10.5 km, middle cloud between 4.0 to 10.5 km, and low cloud between 0.0 to 4.0 km. Cloud may form in a single layer or in multiple contiguous layers; cloud occupying a single layer is treated as radiatively thin, and otherwise as radiatively thick. Clouds are assumed to fill the whole grid box (cloud fraction = 1), and therefore to be fully overlapped in the vertical. Cf. Wetherald and Manabe (1988) [9] and Wetherald et al. (1991) [10] for further details. See also Radiation.


Precipitation from large-scale condensation and from the moist convective adjustment process (see Convection) forms under supersaturated conditions. Subsequent evaporation of falling precipitation is not simulated.

Planetary Boundary Layer

The PBL is simulated only in a rudimentary way (e.g., above the surface layer vertical diffusion of momentum, heat, and moisture are not stability-dependent--see Diffusion). The PBL top is not computed, but its dynamics and physics are typically represented by the first 5 vertical levels above the surface (at about sigma = 0.997, 0.979, 0.935, 0.866, and 0.777). See also Surface Characteristics and Surface Fluxes.


Raw orography data obtained from the U.S. Navy dataset (cf. Joseph 1980 [23]) are smoothed by application of a Fourier filter in the east-west direction, and then are transformed to spectral space at the (R30) model resolution. Subgrid-scale orographic variances required for the gravity-wave drag parameterization are obtained from the same dataset (see Gravity-wave Drag).


AMIP monthly sea surface temperature fields are prescribed, with daily values determined by linear interpolation.

Sea Ice

AMIP monthly sea ice extents are prescribed. Snow may accumulate on sea ice, but does not alter its thermodynamic properties. The surface temperature of sea ice is prognostically determined after Holloway and Manabe (1971) [5] from a surface energy balance (see Surface Fluxes) that includes a conduction heat flux from the ocean below. The conduction flux is proportional to the difference between the surface temperature of the ice and the subsurface ocean temperature (assumed to be at the melting temperature of sea ice, or 271.2 K), and the flux is inversely proportional to the constant ice thickness (2 m). The heat conductivity is assumed to be a constant equal to the value for pure ice, and there is no heat storage within the ice.

Snow Cover

Precipitation falls as snow if the air temperature at 350 m above the surface (sigma = 0.9575) is < 0 degrees C. Snow accumulates on both land and sea ice (see Sea Ice), and snow mass is determined prognostically from a budget equation that accounts for accumulation, melting, and sublimation. Sublimation is calculated as part of the surface evaporative flux. Snowmelt, determined from the excess surface energy available for a surface at temperature 273.2 K, contributes to soil moisture (cf. Holloway and Manabe 1971) [5]. Snow cover affects the albedo of the surface, but not its thermodynamic properties (see also Surface Characteristics, Surface Fluxes, and Land Surface Processes).

Surface Characteristics

  • Roughness lengths over ocean are determined from the surface wind stress after the method of Charnock (1955) [24]. Roughness lengths are prescribed constants for sea ice (1 x 10^-4 m) and land (0.1682 m).

  • Over oceans, the surface albedo depends on solar zenith angle (cf. Payne 1972 [25]). The albedo of snow-free sea ice is a constant 0.50. Albedos for snow-free land are obtained from the data of Posey and Clapp (1964) [26], and do not depend on solar zenith angle or spectral interval.

  • Snow cover modifies the local background albedo of the surface according to its depth, following Holloway and Manabe (1971) [5]. Poleward of 75 degrees latitude, permanent snow with albedo 0.80 is assumed. Equatorward of 75 degrees, snow albedo is set to 0.60 if the water-equivalent snow depth is at least a critical value of 0.01 m; otherwise, the albedo is a combination of the background and snow albedos weighted by the square root of the water-equivalent snow depth.

  • Longwave emissivity is prescribed as unity (blackbody emission) for all surfaces.

Surface Fluxes

  • Surface solar absorption is determined from the surface albedos, and longwave emission from the Planck equation, assuming blackbody emissivity (see Surface Characteristics).

  • Surface turbulent eddy fluxes follow Monin-Obukhov similarity theory, as formulated by Delsol et al. (1971) [27]. The momentum flux is proportional to the product of a drag coefficient, the wind speed, and the wind velocity vector at the lowest atmospheric level (at about sigma = 0.997). Surface sensible heat flux is proportional to the product of a transfer coefficient, the wind speed at the lowest atmospheric level, and the vertical difference between the temperature at the surface and that of the lowest level. The drag and transfer coefficients are functions of stability (bulk Richardson number) and surface roughness length (see Surface Characteristics).

  • The surface moisture flux is the product of potential evaporation and evapotranspiration efficiency beta. Potential evaporation is proportional to the product of the same transfer coefficient as for the sensible heat flux, the wind speed at the lowest atmospheric level, and the difference between the specific humidity at the lowest level and the saturated specific humidity for the local surface temperature and pressure. The evapotranspiration efficiency beta is prescribed to be unity over oceans, snow, and ice surfaces; over land, beta is a function of the ratio of soil moisture to the constant field capacity (see Land Surface Processes).

  • Above the constant-flux surface layer, the vertical diffusion of momentum, heat, and moisture is not stability-dependent (see Diffusion).

Land Surface Processes

  • Ground temperature is determined from a surface energy balance (see Surface Fluxes) without provision for soil heat storage.

  • Soil moisture is represented by the single-layer "bucket" model of Manabe (1969) [4], with field capacity everywhere 0.15 m. Soil moisture is increased by precipitation and snowmelt; it is depleted by surface evaporation, which is determined from a product of the evapotranspiration efficiency beta and the potential evaporation from a surface saturated at the local surface temperature and pressure (see Surface Fluxes). Over land, beta is given by the ratio of local soil moisture to a critical value that is 75 percent of field capacity, and is set to unity if soil moisture exceeds this value. Runoff occurs implicitly if soil moisture exceeds the field capacity.

Go to GFDL References

Return to GFDL Table of Contents

Return to Main Document Directory

Last update April 19, 1996. For further information, contact: Tom Phillips ( )

LLNL Disclaimers