National Center for Atmospheric Research: Model NCAR CCM2 (T42 L18) 1992
AMIP Representative(s)
Dr. David Williamson, National Center for Atmospheric Research, P.O. Box 3000, Boulder, Colorado 80307; Phone: +13034971372; Fax: +13034971324; email: wmson@ncar.ucar.EDU; World Wide Web URL: http://www.ncar.ucar.edu/.Model Designation
NCAR CCM2 (T42 L18) 1992Model Lineage
The NCAR Community Climate Model 2 (CCM2) is the historical descendant of the CCM1 model (cf. Williamson et al. 1987 [1] and Hack et al. 1989 [2]), but with most dynamical and physical parameterizations qualitatively changed. The NCAR CCM2 model is also substantially different from the NCAR GENESIS model.Model Documentation

Key documents are the NCAR CCM2 model description by Hack et al. (1993)
[3] and the user's guide by Bath et al. (1992)
[4]. Other papers that provide details on particular model features include Briegleb (1992)
[5], Briegleb et al. (1986)
[6], and Kiehl and Briegleb (1991)
[7] on the radiation parameterizations, Hack (1994)
[8] on the convection scheme, Holtslag and Boville (1993)
[9] on the simulation of boundarylayer diffusion, and Williamson and Rasch (1994)
[10] on the semiLagrangian transport scheme. Various aspects of the simulated climate with prescribed climatological sea surface temperatures are described by Kiehl et al. (1994)
[11] and Hack et al. (1994)
[12]. Model datasets available for analysis at NCAR (including those from the AMIP simulation) are summarized by Williamson (1993)
[13].
 The complete set of NCAR Technical Notes providing documentation on the CCM2 model and the CCM processor are available at the NCAR Technical Notes Catalog: http://www.ucar.edu/communications/technotes/technotes001100.shtml.WWW summaries of various CCM2 experiments are also available at URL http://www.cgd.ucar.edu/cms/ccm2/391/index.html. In addition, the model source (FORTRAN) code for an improved version of the CCM2 model (designated as CCM2.1, but with algorithms identical to those of the CCM2 model used in the AMIP experiment) can be downloaded by anonymous ftp at ftp.ucar.edu/pub/ccm.
Numerical/Computational Properties
Horizontal Representation
Spectral (spherical harmonic basis functions) with transformation to a Gaussian grid for calculation of nonlinear quantities and most of the physics. Advection of water vapor is via shapepreserving semiLagrangian transport (SLT) on the Gaussian grid (cf. Williamson and Rasch 1994).Horizontal Resolution
Spectral triangular 42 (T42), roughly equivalent to 2.8 x 2.8 degrees latitudelongitude.Vertical Domain
Surface to 2.917 hPa. For a surface pressure of 1000 hPa, the lowest atmospheric level is at a pressure of about 992 hPa.Vertical Representation
Finite differences in hybrid sigmapressure coordinates after Simmons and Striifing (1981) [14], but modified to allow an upper boundary at nonzero (2.917 hPa) pressure. The verticaldifferencing formulation conserves global total energy in the absence of sources and sinks. See also Vertical Domain and Vertical Resolution.Vertical Resolution
There are 18 unevenly spaced hybrid sigmapressure levels. For a surface pressure of 1000 hPa, 4 levels are below 800 hPa and 7 levels are above 200 hPa.Computer/Operating System
The AMIP simulation was run on Cray 2 computers using multiple processors in the UNICOS environment.Computational Performance
For the AMIP experiment, about 7 minutes on a single processor of the Cray 2 computer per simulated day.Initialization
For the AMIP simulation, the model atmosphere is initialized from a previous model solution for 10 December 1978, and is then "spun up" to a simulated 1 January 1979 state. Snow cover/depth are prescribed from a mean January climatology only as a function of latitude (see Snow Cover). Annualaverage ground wetness is prescribed for 10 surface types in place of specifying soil moisture (see Surface Characteristics and Land Surface Processes).Time Integration Scheme(s)
A centered semiimplicit time integration scheme (cf. Simmons et al. 1978 [15] ) with an Asselin (1972)[16] frequency filter is used for many calculations, but horizontal and vertical diffusion (see Diffusion), the advection of water vapor by the SLT scheme (see Horizontal Representation), and adjustments associated with convection and largescale condensation (see Convection and Cloud Formation) are computed implicitly by a timesplitting procedure. The overall time step is 20 minutes for dynamics and physics, except for shortwave and longwave radiative fluxes and heating rates, which are calculated hourly (with longwave absorptivities and emissivities updated every 12 hourssee Radiation). Cf. Hack et al. (1993) [3] for further details.Smoothing/Filling
Orography is smoothed (see Orography). Because advection of moisture is treated by the SLT scheme (see Horizontal Representation) negative specific humidity values are avoided. In cases where negative mixing ratios would result from application of the countergradient term in the parameterization of nonlocal vertical diffusion of moisture in the planetary boundary layer (PBL) (see Diffusion, Planetary Boundary Layer, and Surface Fluxes), the countergradient term is not calculated. In addition, at each 20minute time step a "fixer" is applied to the surface pressure and water vapor so that the global average mass and moisture are conserved (cf. Williamson and Rasch 1994) [10].Sampling Frequency
For the AMIP simulation, the model history is written every 6 hours.Dynamical/Physical Properties
Atmospheric Dynamics
Primitiveequation dynamics are expressed in terms of vorticity, divergence, temperature, specific humidity, and the logarithm of surface pressure. Virtual temperature is used where applicable, and frictional/diffusive heating is included in the thermodynamic equation.Diffusion
 In the troposphere, linear biharmonic (del^4) horizontal diffusion (with coefficient 1 x 10^16 m^4/s) is applied to divergence and vorticity on hybrid sigmapressure surfaces, and to temperature on firstorder constant pressure surfaces (requiring that biharmonic diffusion of surface pressure also be calculated on the Gaussian grid). In the stratosphere linear secondorder (del^2) diffusion is applied to the same variables at the top three levels (with diffusivities increasing with height from 2.5 x 10^5 to 7.5 x 10^5 m^2 /s). In the top model layer, diffusion is enhanced by a factor of 10^3 on all spectral wave numbers that violate the CourantFriedrichsLewy (CFL) numerical stability criterion, based on the maximum wind speed.
 Above the PBL (see Planetary Boundary Layer) a secondorder, stabilitydependent local formulation of the vertical diffusion of momentum, heat, and moisture is adopted (cf. Smagorinsky et al. 1965
[17]). The mixing length is taken to be a constant 30 m, and the diffusivity is as given by Williamson et al. (1987)
[1] for unstable and neutral conditions and by Holtslag and Beljaars (1989)
[33] for stable conditions. Above the surface layer, but within the PBL under unstable conditions, mixing of heat and moisture (but not of momentum) is formulated as nonlocal diffusion, following Holtslag and Boville (1993)
[9]see
Surface Fluxes.
 Horizontal and vertical diffusion are calculated implicitly via time splitting apart from the solution of the semiimplicit dynamical equations (see Time Integration Scheme(s)).
Gravitywave Drag
Orographic gravitywave drag is parameterized after McFarlane (1987) [18]. The momentum drag is given by the vertical divergence of the wave stress, which is proportional to the product of the local squared amplitude of the gravity wave, the BruntVaisalla frequency, and the component of the local wind that is parallel to the flow at a nearsurface reference level. At this reference level, the wave amplitude is bound by the lesser of the subgridscale orographic variance (see Orography) or a wavesaturation value defined by the reference Froude number. Above this level, the gravitywave stress is assumed to be constant with height (zero vertical divergence), except in regions of wave saturation, where the amplitude is obtained from the local Froude number.Solar Constant/Cycles
The solar constant is the AMIPprescribed value of 1365 W/(m^2). Both seasonal and diurnal cycles in solar forcing are simulated over a repeatable solar year of exactly 365 days (i.e., leap years are not included).Chemistry
The carbon dioxide concentration is the AMIPprescribed value of 345 ppm. Monthly ozone volume mixing ratios derived by Chervin (1986) [19] from analyses of Dütsch (1978) [38] are linearly interpolated to obtain intermediate values every 12 hours. Radiative effects of oxygen and of water vapor, but not of aerosols, are also included (see Radiation).Radiation
 Shortwave scattering/absorption is parameterized by the deltaEddington approximation of Joseph et al. (1976)
[20] and Coakley et al. (1983)
[21] applied in 18 spectral intervals, as described by Briegleb (1992)
[5]. (These include 7 intervals between 0.20 and 0.35 micron to capture ozone HartleyHuggins band absorption and Rayleigh scattering; 1 interval between 0.35 to 0.70 micron to capture Rayleigh scattering and ozone Chappiusband and oxygen Bband absorption; 7 intervals between 0.70 and 5.0 microns to capture oxygen Aband and water vapor/liquid absorption; and 3 intervals between 2.7 and 4.3 microns to capture carbon dioxide absorption.) Following Slingo (1989)
[22], the shortwave optical properties of clouds for the deltaEddington approximation (optical depth, singlescattering albedo, asymmetry factor) are specified for 4 spectral ranges (with boundaries at 0.25, 0.69, 1.19, 2.38, and 4.0 microns). These properties depend on the specified effective droplet radius (10 microns) and the liquid water path (LWP), which is a prescribed nonlinear function of latitude and height (cf. Kiehl et al. 1994)
[11].
 Longwave absorption by ozone and carbon dioxide is treated by a broadband absorptance technique, following Ramanathan and Dickinson (1979) [23] and Kiehl and Briegleb (1991) [7]. A Voigt line profile (temperature) dependence is added to the pressure broadening of the absorption lines. Absorption by water vapor (and its overlap with that of ozone and carbon dioxide) are modeled as in Ramanathan and Downey (1986) [24]. Longwave broadband emissivity of clouds is a negative exponential function of LWP, with all clouds assumed to be randomly overlapped in the vertical. Cf. Hack et al. (1993) [3] for further details. See also Cloud Formation.
Convection
 If the atmosphere is moist adiabatically unstable, temperature/moisture column profiles are adjusted by a massflux convective parameterization (cf. Hack 1994)
[8]. The scheme utilizes a threelayer model that provides for convergence and entrainment in the lowest subcloud layer, cloud condensation and rainout in the middle layer, and limited detrainment in the top layer. This scheme is applied by working upward from the surface on three contiguous layers, and shifting up successively one layer at a time until the whole column is stabilized.
 The parameterization is based on simplified equations for the threelayer moist static energy that include (among other terms) the convective mass flux, a "penetration parameter" beta (ranging between 0 and 1) that regulates the detrainment of liquid water, and temperature and moisture perturbations furnished by the PBL parameterization (see Planetary Boundary Layer, Diffusion, and Surface Fluxes). Other free parameters in the scheme include minimum values for beta, for the vertical gradient of moist static energy, and for the depth of precipitating convection; a characteristic convective adjustment time scale; and a cloudwater to rainwater autoconversion coefficient. The parameter beta is determined by iteration, subject to constraints that it and the vertical gradient of moist static energy be at least their minimum values, that the convective mass flux be positive, and that the detrainment layer not be supersaturated. The profiles of convective mass flux, temperature, and moisture are then obtained, and the total convective precipitation rate is calculated by vertical integration of the convectivescale liquid water sink.
 If a layer in the stratosphere (i.e., at the top three vertical levels) is dry adiabatically unstable, the temperature is adjusted so that stability is restored under the constraint that sensible heat be conserved. Whenever two layers undergo this dry adjustment, the moisture is also mixed in a conserving manner. (In the model troposphere, vertical diffusion provides stabilizing mixing, and momentum is mixed as wellsee Diffusion). If a layer is supersaturated but stable, nonconvective condensation and precipitation result (see Precipitation).
Cloud Formation
 The cloud prediction formulation is as described by Kiehl et al. (1994)
[11]. Cloud amount is diagnostically determined from relative humidity, vertical velocity, atmospheric stability, and the convective precipitation rate, following a modified Slingo (1987)
[25] approach. Convective cloud, layer cloud, and lowlevel marine stratus/stratocumulus cloud associated with temperature inversions are treated. Clouds may form everywhere except in the surface layer (centered at sigma = 0.992). Nonconvective cloud below 700 hPa is classified as lowlevel cloud, that between 400 and 700 hPa as midlevel cloud, and that above 400 hPa as highlevel cloud.
 Convective cloud base and top are determined by the vertical extent of moist instability (see Convection). In each vertical column, the total fractional cloud amount is a logarithmic function of the convective precipitation rate, but is constrained to be between 0.2 and 0.8. The convective cloud fraction in each layer is determined assuming the cloud is distributed randomly in the vertical. For subsequent diagnosis of the fractional amount of nonconvective cloud (see below), the layer relative humidity is reduced proportional to the fraction of convective cloud present.
 In regions of upward vertical motion, the fraction of lowlevel layer cloud is a quadratic function of the difference between the reduced relative humidity (see above) and a constant threshold value (90 percent). The fraction of midlevel and highlevel layer cloud is a quadratic function of the difference between the reduced relative humidity and a threshold value that is a linear function of the squared BruntVaisalla frequency (i.e., it is proportional to the vertical stability).
 The fraction of marine stratus/stratocumulus is a function of the strength of the associated lowlevel inversion and the reduced relative humidity. Cf. Hack et al. (1993) [3] for further details. See also Radiation.
Precipitation
Subgridscale precipitation is generated in unstable conditions by the moist convective scheme (see Convection). Gridscale precipitation forms as a result of supersaturation under stable conditions. In this case, the moisture is adjusted so that the layer is just saturated, with the excess condensing as precipitation; the layer temperature is adjusted according to the associated latent heat release. (Moisture and temperature are mutually adjusted in two iterations.) Subsequent evaporation of falling precipitation is not simulated. Cf. Hack et al. (1993) [3] for details.Planetary Boundary Layer
The PBL height is determined by iteration at each 20minute time step following the formulation of Troen and Mahrt (1986) [26]; the height is a function of the critical bulk Richardson number for the PBL, uv winds and virtual temperature at the PBL top, and the 10meter virtual temperature, which is calculated from the temperature and moisture of the surface and of the lowest atmospheric level (at sigma = 0.992) following Geleyn (1988) [27]. Within the PBL, there is nonlocal diffusion of heat and moisture after Holtslag and Boville (1993) [9]; otherwise (and under all conditions for momentum), properties are mixed by the stabilitydependent local diffusion that applies in the model's free atmosphere. See also Diffusion and Surface Fluxes.Orography
 Raw orography is obtained from the U.S. Navy dataset with resolution of 10 minutes arc on a latitude/longitude grid (cf. Joseph 1980
[28]). These data are areaaveraged to a 1 x 1degree grid, interpolated to a T119 Gaussian grid, spectrally truncated to the model's T42 Gaussian grid, and then spectrally filtered to reduce the amplitude of the smallest scales.
 The subgridscale orographic variances required for the gravitywave drag parameterization (see Gravitywave Drag) are also obtained from the U.S. Navy dataset. For the spectral T42 model resolution, the variances are first evaluated on a 2 x 2degree grid, assuming they are isotropic. Then the variances are binned to the T42 Gaussian grid (i.e., all values whose latitude and longitude centers fall within each Gaussian grid box are averaged together), and are smoothed twice with a 121 spatial filter. Values over ocean are set to zero.
Ocean
AMIP monthly sea surface temperature fields are prescribed, with intermediate values determined at every 20minute time step by linear interpolation.Sea Ice
AMIP monthly sea ice extents are prescribed, with intermediate values determined at every 20minute time step by linear interpolation. The temperature of the ice is predicted by the same fourlayer scheme as used for soil temperature (see Land Surface Processes), but with a fixed temperature (2 degrees C) of the underlying ocean rather than a zeroflux condition, as the lower boundary condition. The four layer thicknesses are all 0.5 m, and the ice density, heat capacity, and conductivity are specified uniform constants; however, daily snow cover that is prescribed from climatology (see Snow Cover) alters the thermodynamic properties and thickness of the top layer in proportion to the relative mass of snow and ice. Cf. Hack et al. (1993) [3] for further details.Snow Cover
Snow cover (expressed as an equivalent depth of water) is prescribed as a function of latitude and longitude from the mean January and July data of Forderhase et al. (1980) [29] that are bilinearly interpolated to the T42 Gaussian grid. Intermediate daily values are obtained by assuming a singleharmonic annual variation. Snow cover is prescribed on sea ice as well as land, and affects the albedo, the roughness and wetness, and the thermodynamics of the surface (see Surface Characteristics). In addition, sublimation of snow contributes to the surface evaporative flux (see Surface Fluxes); however, because ground wetness is prescribed, snowmelt does not affect soil hydrology. See also Sea Ice and Land Surface Processes.Surface Characteristics
 Over land, surface wetness fractions, roughness lengths, and albedos are derived from the Matthews (1983)
[30] 1 x 1degree, 32type vegetation data set. These values are aggregated to the 10 surface types (including land ice) distinguished by the model and are averaged over the T42 Gaussian grid boxes.
 The prescribed annualaverage wetness fractions for the 10 surface types range between 0.01 for deserts to 1.0 for open water and ice and snowcovered surfaces. On land, the surface wetness is weighted by the local fractional area of snow, which depends both on snow depth and the surface roughness length (to account for uneven coverage of vegetation). The heat capacity and conductivity of six distinguished landsurface thermal types also depend on surface wetness. See also Sea Ice and Land Surface Processes
 Over land, the surface roughness length ranges from 0.04 m for tundra to 1.0 m for evergreen forest. The roughness is a uniform 1 x 10^4 m over ocean, and 0.04 m over ice surfaces (cf. Hack et al. 1993
[3] for further details).
 The snowfree land albedos are constants (independent of time or moisture conditions) for the 10 distinguished surfaces. These are composed of five quantities: the fraction of strong zenithangle dependence and four surface albedos (for two zenith angles and spectral intervals 0.2 to 0.7 micron and 0.7 to 4.0 microns). The land albedo is altered by snow cover: it is an average of the background albedo and the snow albedo (which depends on surface temperature for the diffuse beam and on solar zenith angle for the direct beam) that is weighted by the fractional snow cover (see above). Over the oceans, surface albedos are prescribed to be 0.025 for the directbeam (with sun overhead) and 0.06 for the diffusebeam component of radiation; the directbeam albedo varies with solar zenith angle. The albedo of ice surfaces is a function of surface temperature. Cf. Briegleb et al. (1986)
[6], Briegleb (1992)
[5] and Dickinson et al. (1986)
[31] for further details.
 The longwave emissivity is set to unity (blackbody emission) for all surfaces.
Surface Fluxes
 Surface solar absorption is determined from surface albedo, and longwave emission from the Planck equation with prescribed surface emissivity of 1.0 (see Surface Characteristics).
 Turbulent vertical eddy fluxes of momentum, heat, and moisture are expressed as bulk formulae, following MoninObukhov similarity theory. The values of wind, temperature, and humidity required for the bulk formulae are taken to be those at the lowest atmospheric level (sigma = 0.992), which is assumed to be within a constantflux surface layer. The drag and transfer coefficients in the bulk formulae are functions of roughness length (see Surface Characteristics) and stability (bulk Richardson number), following the method of Louis et al. (1981)
[32] for neutral and unstable conditions, and Holtslag and Beljaars (1989)
[33] for stable conditions. The bulk formula for the surface moisture flux also includes a prescribed surface wetness fraction (see Surface Characteristics and Land Surface Processes) that determines the evaporation realized as a fraction of potential evaporation from a saturated surface.
 Above the surface layer, but within the PBL (see Planetary Boundary Layer) under unstable conditions, mixing of heat and moisture (but not of momentum) is formulated as nonlocal vertical diffusion by eddies with length scales of the order of the PBL depth (cf. Deardorff 1972 [34]). Under these conditions, a countergradient term that depends on the surface flux, a convective vertical velocity scale, and the PBL height is added to the eddy diffusivity coefficient of heat and moisture. Within the stable and neutral PBL (and under all conditions for momentum), the same stabilitydependent local vertical diffusion as is utilized in the model's free atmosphere applies (see Diffusion). Cf. Holtslag and Boville (1993) [9] for further details.
Land Surface Processes
 Soil temperature is determined from heat conduction in a fourlayer model. Layer thicknesses vary spatially, depending on the penetration depth of solar forcing, which is a function of forcing period and of the heat capacity and conductivity specified for each of the 10 distinguished surface types (see Surface Characteristics). The thicknesses of the bottom three soil layers are specified according to the local penetration depth of solar forcing with periods of 1 day, 2 weeks, and 1 year respectively. The thickness of the top soil layer is specified so that the diurnal range and phase of the surface temperature compares well with observations cited by Bhumralkar (1975)
[35]. The top layer's thickness and heat capacity/conductivity are altered by snow (see Snow Cover) in proportion to the relative masses of snow and soil; these thermodynamic properties are also a linear function of ground wetness fraction (see below). The heat conduction equation is solved by a backward implicit CrankNicholson numerical scheme (cf. Smith 1965
[36] and Washington and Verplank 1986
[37]), with the net surface energy balance being the upper boundary condition (see Surface Fluxes), and with zero heat flux specified at the lower boundary.
 Soil moisture is not predicted, and precipitation and snowmelt therefore do not affect the land surface hydrology. Instead, constant wetness fractions are specified for the 10 distinguished landsurface types (see Surface Characteristics). These wetness fractions affect the heat capacity and conductivity of the soil (see above), and they constrain the magnitude of the surface evaporative flux (see Surface Fluxes). Cf. Hack et al. (1993) [3] for further details.
Return to NCAR Table of Contents
Return to Main Document Directory
Last update September 18, 1996. For further information, contact: Tom Phillips ( phillips@tworks.llnl.gov )
UCRLID116384