The surface exchange of heat, moisture and momentum between the atmosphere and land, ocean or ice surfaces are treated with a bulk exchange formulation. We present a description of each surface exchange separately. Although the functional forms of the exchange relations are identical, we present the descriptions of these components as developed and represented in the various subroutines in CAM 3.0. The differences in the exchange expressions are predominantly in the definition of roughness lengths and exchange coefficients. The description of surface exchange over ocean follows from Bryan et al. , and the surface exchange over sea ice is discussed in chapter 6. Over lakes, exchanges are computed by a lake model embedded in the land surface model described in the following section.
In CAM 3.0, the NCAR Land Surface Model (LSM)  has been
replaced by the Community Land Model CLM2 . This new
model includes components treating hydrological and biogeochemical
processes, dynamic vegetation, and biogeophysics. Because of the
increased complexity of this new model and since a complete
description is available online, users of CAM 3.0 interested in CLM
should consult this documentation at
http://www.cgd.ucar.edu/tss/clm/. A discussion is provided here only of the component of CLM which controls surface exchange processes.
Land surface fluxes of momentum, sensible heat, and latent heat are calculated from Monin-Obukhov similarity theory applied to the surface (i.e. constant flux) layer. The zonal and meridional momentum fluxes (kg ms), sensible heat (W m) and water vapor (kg ms) fluxes between the surface and the lowest model level are:
For the vegetated fraction of the grid, and , where and are the air temperature and specific humidity within canopy space. For the non-vegetated fraction, and , where and are the air temperature and specific humidity at ground surface. These terms are described by Dai et al. .
The aerodynamic roughness is used for wind, while the thermal roughness is used for heat and water vapor. In general, is different from , because the transfer of momentum is affected by pressure fluctuations in the turbulent waves behind the roughness elements, while for heat and water vapor transfer no such dynamical mechanism exists. Rather, heat and water vapor must ultimately be transferred by molecular diffusion across the interfacial sublayer. Over bare soil and snow cover, the simple relation from Zilitinkevich  can be used :
The roughness is proportional to canopy height, and is also affected by fractional vegetation cover, leaf area index, and leaf shapes. The roughness is derived from the simple relationship , where is the canopy height. Similarly, the zero-plane displacement height is proportional to canopy height, and is also affected by fractional vegetation cover, leaf area index, and leaf shapes. The simple relationship is used to obtain the height.
(1) Turbulence scaling parameters
A length scale (the Monin-Obukhov length) is defined by
The stability parameter is defined as
(2) Flux-gradient relations 
The flux-gradient relations are given by:
Under very unstable conditions, the flux-gradient relations are taken from Kader and Yaglom :
To ensure the functions and are continuous, the simplest approach (i.e., without considering any transition regions) is to match the above equations at for and for .
Under very stable conditions (i.e., ), the relations are taken from Holtslag et al. :
(3) Integral forms of the flux-gradient relations
Integration of the wind profile yields:
Integration of the potential temperature profile yields:
The expressions for the specific humidity profiles are the same as those for potential temperature except that ( ), and are replaced by ( ), and respectively. The stability functions for are
Note that the CLM code contains extra terms involving , , and for completeness. These terms are very small most of the time and hence are omitted in Eqs. 4.434 and 4.435.
In addition to the momentum, sensible heat, and latent heat fluxes, land surface albedos and upward longwave radiation are needed for the atmospheric radiation calculations. Surface albedos depend on the solar zenith angle, the amount of leaf and stem material present, their optical properties, and the optical properties of snow and soil. The upward longwave radiation is the difference between the incident and absorbed fluxes. These and other aspects of the land surface fluxes have been described by Dai et al. .
The bulk formulas used to determine the turbulent fluxes of momentum (stress), water (evaporation, or latent heat), and sensible heat into the atmosphere over ocean surfaces are
In (4.439), the transfer coefficients between the ocean surface and the atmosphere are computed at a height and are functions of the stability, :
Over oceans, m under all conditions and m for , m for , which are given in Large and Pond . The momentum roughness length depends on the wind speed evaluated at 10 m as
The transfer coefficients in (4.439) and (4.440) depend on the stability following (4.441)-(4.444), which itself depends on the surface fluxes (4.445) and (4.446). The transfer coefficients also depend on the momentum roughness, which itself varies with the surface fluxes over oceans (4.447). The above system of equations is solved by iteration.