Parameterizations for Cloud Overlapping and Shortwave Single-Scattering Properties for Use in General Circulation and Cloud Ensemble Models

Transcription

1 202 JOURNAL OF CLIMATE Parameterizations for Cloud Overlapping and Shortwave Single-Scattering Properties for Use in General Circulation and Cloud Ensemble Models MING-DAH CHOU AND MAX J. SUAREZ Laboratory for Atmospheres, NASA/Gpddard Sapce Flight Center, Greenbelt, Maryland CHANG-HOI HO Universities Space Research Association, Columbia, Maryland MICHAEL M.-H. YAN Science Systems and Applications, Inc., Seabrook, Maryland KYU-TAE LEE Department of Atmospheric and Environmental Sciences, Kangnung National University, Kangnung, South Korea (Manuscript received 7 December 996, in final form 6 June 997) ABSTRACT Parameterizations for cloud single-scattering properties and the scaling of optical thickness in a partial cloudiness condition have been developed for use in atmospheric models. Cloud optical properties are parameterized for four broad bands in the solar (or shortwave) spectrum; one in the ultraviolet and visible region and three in the infrared region. The extinction coefficient, single-scattering albedo, and asymmetry factor are parameterized separately for ice and water clouds. Based on high spectral-resolution calculations, the effective single-scattering coalbedo of a broad band is determined such that errors in the fluxes at the top of the atmosphere and at the surface are minimized. This parameterization introduces errors of a few percent in the absorption of shortwave radiation in the atmosphere and at the surface. Scaling of the optical thickness is based on the maximum-random cloud-overlapping approximation. The atmosphere is divided into three height groups separated approximately by the 400- and 700-mb levels. Clouds are assumed maximally overlapped within each height group and randomly overlapped among different groups. The scaling is applied only to the maximally overlapped cloud layers in individual height groups. The scaling as a function of the optical thickness, cloud amount, and the solar zenith angle is derived from detailed calculations and empirically adjusted to minimize errors in the fluxes at the top of the atmosphere and at the surface. Different scaling is used for direct and diffuse radiation. Except for a large solar zenith angle, the error in fluxes introduced by the scaling is only a few percent. In terms of absolute error, it is within a few watts per square meter.. Introduction Clouds have a large effect on the radiative heating cooling of the atmosphere and the earth s surface. Difficulties in modeling clouds involve the parameterizations of cloud generation, development, and dissipation, as well as cloud microphysical properties (water content, phase, particle shape, and size distribution) and optical properties (extinction coefficient, single-scattering albedo, scattering phase function). The microphysical properties are determined by the thermal and dynamical Corresponding author address: Dr. Ming-Dah Chou, Code 93, Laboratory for Atmospheres, NASA/Goddard Space Flight Center, Greenbelt, MD chou@climate.gsfc.nasa.gov properties of the environment, whereas the cloud optical properties are determined by cloud microphysical properties. In atmospheric models, cloud optical properties are commonly parameterized as functions of particle size and water content. Difficulties involved in the parameterization for cloud single-scattering properties are many. The single-scattering coalbedo of ice crystals is significantly larger than that of water droplets; thus, ice clouds have a larger absorption than water clouds for a given optical thickness. The optical properties of ice clouds are functions of the size, shape, and orientation of ice crystals, which vary over a large range and cannot be determined in current atmospheric models. The single-scattering coalbedo of ice and water particles vary rapidly with wavelength. Parameterization for the mean effective singlescattering coalbedo for a wide spectral interval depends 998 American Meteorological Society

2 FEBRUARY 998 CHOU ET AL. 203 not only on the absorption by clouds but also on the water vapor within and above cloud layers (Ramaswamy and Freidenreich 992). All these problems make the parameterizations for cloud optical properties very difficult. Calculations of radiative heating are further complicated by the fact that clouds are not horizontally homogeneous and do not always cover the entire sky. The solar (or shortwave, SW) radiation algorithms used in atmospheric models are developed under the assumption that atmospheric layers are plane parallel and can only be applied to the case where an atmospheric layer is either totally cloud covered or cloud free. The simplest approach to dealing with the partially cloudy case is to assume that clouds in various layers overlap either randomly or maximally. The sky is then divided into sections. Within each section, clouds are homogeneous in a layer with cloud amount equal to either 0 or. Shortwave fluxes are calculated for each section and then weighted by the cloud amount to derive total fluxes. This approach requires a huge amount of computations. For the case of random cloud overlapping, it requires 2 n sets of calculations, where n is the number of cloud layers. For the case of maximum cloud overlapping, the computing time required is less. To reduce the computational burden, clouds in a partially covered layer are smeared to cover the entire layer, and the optical thickness is scaled by a function of the cloud amount (e.g., Kiehl et al. 994; Sud et al. 993). In addition to being a function of the cloud amount, scaling of the cloud optical thickness should also be a function of the optical thickness itself, the solar zenith angle, as well as the assumption made for cloud overlapping. In this study, we present parameterizations for cloud single-scattering properties, for the treatment of cloud overlapping, and for the scaling of cloud optical thickness. These parameterizations have been implemented in the GEOS (Goddard Earth Observing System) general circulation model (GCM) (Schubert et al. 993) and the Goddard cloud ensemble model (Tao et al. 996). The parameterizations for cloud single-scattering properties are extensions of those of Slingo (989) for water clouds and Fu (996) for ice clouds. The scaling of cloud optical thickness is based on the assumption of maximum cloud overlapping. Uncertainties in these cloud-radiation schemes are investigated by comparing the results with more detailed calculations. 2. Radiative transfer calculations In our radiation routine, the solar spectrum is divided into an ultraviolet (UV) and visible region ( 0.7 m) and an IR region (0.7 m 0 m). In the UV and visible region, we include absorption due to O 3 and aerosols and scattering due to gases, clouds, and aerosols. The UV and visible region is further divided into eight spectral intervals so that fluxes in the PAR (photosynthetically active radiation) ( TABLE. Spectral range and the weight in (6), h, for the four spectral bands. Note: In our radiation scheme, Band is divided into eight subbands. Ozone and Rayleigh extinction coefficients are given for each subband, and fluxes are computed for each subband. Band Spectral range (m) h (ice cloud) 2/3 /3 0 h (water cloud) m), UV-A ( m), UV-B ( m), and UV-C ( m) regions can be separately computed. An effective absorption coefficient for O 3, and an effective extinction coefficient for Rayleigh scattering are derived for each spectral interval. It is noted that, in this study, a single set of parameterizations for cloud optical properties is applied to the eight spectral intervals in the UV and visible region (Band shown in Table ). In the IR region, we include absorption due to water vapor, O 2,CO 2, clouds, and aerosols, and scattering due to clouds and aerosols. Calculations of the absorption due to water vapor and O 3 follow Chou and Lee (996), and those of the absorption due to O 2 and CO 2 follow Chou (990). This spectral region is divided into three bands (Table ). To avoid oscillation in the computed heating profile caused by the use of a small number of k intervals (Chou and Lee 996), each band is grouped into 0 k intervals, where k is the water vapor absorption coefficient. A one-parameter scaling approximation is applied to take into account the effect of temperature and pressure on k. To reduce computing time, we avoid applying a single multiple-scattering algorithm to the entire inhomogeneous atmosphere for flux calculations. Instead, we first compute the transmission and reflection functions for each layer and then use the two-stream adding method to compute fluxes in the atmosphere and at the surface (cf. Lacis and Hansen 974). In this approach, the transmission and reflection functions for each layer are separately computed for direct and diffuse incident radiation [cf. (4) and (5) of Chou 992]. In our radiation model, the reflection and transmission functions of a scattering layer are computed using the -Eddington approximation of Joseph et al. (976) for direct radiation and the two-stream approximation of Sagan and Pollack (967) for diffuse radiation. 3. Parameterizations for cloud single-scattering properties The cloud single-scattering properties that affect radiative transfer are the extinction coefficient, singlescattering albedo, and scattering phase function (or asymmetry factor). The extinction coefficient per unit mass of polydispersion cloud particles is given by 0 0

3 204 JOURNAL OF CLIMATE FIG.. Spectral distributions of the single-scattering coalbedo of ice and water clouds. The values of r i and r w are, respectively, the effective size of ice and water particles. (r)n(r) dr/c, () where (r) is the extinction cross section of a particle with a size r, the subscript is the wavelength, n(r)dr is the number of particles within the size range dr per unit volume, and C is the cloud water mass concentration. The single-scattering albedo is defined as the ratio of the scattering coefficient to the total extinction coefficient s s s a / /( ), (2) where the superscripts s and a denote scattering and absorption, respectively. The asymmetry factor is given by g P () d, (3) 2 where P is the scattering phase function of polydispersion cloud particles and is the cosine of the scattering angle. Radiation routines for atmospheric models usually use only a few broad spectral bands with effective optical properties parameterized for each band. The singlescattering coalbedo,, varies by several orders of magnitude in the solar spectrum. It is, therefore, difficult to derive an effective for a broad spectral band that can be applied to both strong and weak absorption situations. There are a number of approaches to deriving. Slingo and Schrecker (982) derived by weighting with the solar insolation and linearly averaging over a spectral band. Espinoza and Harshvardhan (996) used a square-root approximation to derive. Fu (996) used a mix of linear and logarithmic averaging depending on the strength of absorption in a given spectral band. The absorption of SW radiation by clouds depends not only on but also on water vapor and cloud ice/water amounts. It is clear that there is no unique method for deriving over a broad band. It can only be empirically determined based on the amounts of cloud particles and water vapor encountered in the earth s atmosphere. Let us define the linear and logarithmic averaging over a band as and ( ) ( ) s s (4) log( ) log( ) s s, (5) where s is the solar insolation at the top of the atmosphere (TOA) and is a narrow spectral interval, where optical properties can be treated as constants. The effective mean single-scattering coalbedo of a band is then computed from ( ) h( ) ( h)( ), (6) where h is the weight ranging from 0 to. The weight

4 FEBRUARY 998 CHOU ET AL. 205 FIG. 2. Same as Fig. except for the cloud asymmetry factor. h should be close to for the spectral bands with weak absorption and should decrease as absorption increases. It is to be determined empirically. Compared to, the extinction coefficient ( ) and the asymmetry factor (g ) vary rather smoothly with. Their effective mean values over a wide spectral band can be accurately approximated by s s and (7) g g s s. (8) Theoretical considerations and radiative transfer calculations have shown that cloud single-scattering properties are not significantly affected by details of the particle size distribution and can be adequately parameterized as functions of the effective particle size (Fu 996; Hu and Stamnes 993; Tsay et al. 989). Following Slingo and Schrecker (982), we parameterize the single-scattering properties for a broad spectral band by ao a /r e, (9) 2 bo br e br, 2 e and (0) 2 g co cr e cr, 2 e () where a, b, and c are regression coefficients, and r e is the effective particle size defined to be proportional to the ratio of the total volume of cloud particles to the total cross-sectional area, A c, of cloud particles. For spherical water droplets, the effective size is given by C rw rn(r) dr r n(r) dr, (2) 4w Ac where w is the density of water. For hexagonal ice crystals randomly oriented in space, the effective size is shown by Fu (996) to be 23 C ri, (3) 3 A where i is the density of ice. The cloud optical thickness,, is then given by Cz, (4) where z is the geometric thickness of a cloud layer. The spectral data,, and g calculated by Fu (996) for ice clouds and by Tsay et al. (989) for water clouds are used to derive,, and g from (4) (8). By assuming hexagonal ice crystals randomly oriented in space, Fu (996) computed the single-scattering parameters of ice clouds using the improved ray-tracing method of Yang and Liou (995). A total of 28 size distributions derived from in situ aircraft measurements were used in the calculations, which included samples from the First ISCCP Regional Experiment (FIRE) and Central Equatorial Pacific Experiment (CEPEX) field campaigns. The mean effective size of ice crystals, r i, ranges from 20 to 30 m. For water clouds, the single-scat- i c

5 206 JOURNAL OF CLIMATE TABLE 2. Coefficients a 0 and a of the parameterization, (9), for the cloud extinction coefficient. The units of a 0 and a are, respectively, m 2 g and m 2 g m. Spectral band Ice cloud Water cloud a 0 a a 0 a n n n n FIG. 3. The extinction coefficient of ice and water clouds as a function of the effective particle size. The values of r i and r w are, respectively, the effective size of ice and water particles. The data points are derived from (7), whereas the curves are the regression using (9). tering parameters of five water clouds were computed by Tsay et al. (989) from the Mie theory assuming spherical droplets and a log-normal size distribution. The effective radius of water droplets, r w, ranges from 4to20m. The size, shape, and refractive indices are different for ice crystals and water droplets. The extinction coefficient of ice clouds is smaller than that of water clouds, because ice crystals are much larger than water droplets. The single-scattering albedo and asymmetry factor of ice clouds and water clouds are also different, as shown in Figs. and 2. The ice cloud optical properties shown in the figures were computed by Fu (996) for a FIRE I (22 October 986) case and a CEPEX case, which have a mean effective size, r i,of45m and 97 m, respectively. For the water cloud, two log-normal size distributions with effective droplet radii, r w,of8 and 6 m are used to compute the cloud optical properties (Tsay et al. 989). It can be seen in the figures that the single-scattering coalbedo,, of ice crystals is larger than that of water droplets by nearly a factor of 0, while the asymmetry factor of ice crystals is smaller than that of water droplets for.4 m. Because of these differences, we parameterize the single-scattering properties separately for ice and water clouds. As in Slingo (989), we divide the solar spectrum into four wide bands and parameterize the single-scattering properties for these bands. The spectral ranges of these four bands are given in Table. There is one band in the ultraviolet and visible spectral regions and three in the IR region. Referring to Fig., the single-scattering coalbedo is small (0, weak absorption) in the first two bands but is large (3, strong absorption) in Band 4. Within individual bands, it varies by two orders of magnitude. Figure 2 shows that the asymmetry factor varies more strongly with wavelength for ice clouds than for water clouds. The extinction coefficient varies with wavelength weakly for both ice and water clouds, and the detailed spectral distributions are not shown in the figures. The single-scattering parameters,,, and g are computed, respectively, from (4), (5), (7), and (8) for the four spectral bands and all 33 cloud cases (28 ice clouds and 5 water clouds). By specifying various values of the weight h, the effective mean single-scattering albedo of a band are computed from (6), also for the four spectral bands and 33 cloud cases. These bandaveraged single-scattering parameters are then fit by (9) () separately for ice and water clouds. Fluxes are computed using these band-averaged single-scattering properties for a midlatitude summer atmosphere taken form McClatchey et al. (972) with clouds located at various heights in the atmosphere. Optimal values of h are then determined empirically (trial and error) such that the difference in the fluxes at the TOA and at the surface between the parameterization (four bands, Table ) and high-spectral resolution (00 cm ) calculations is minimized. For the high-spectral resolution calculations, the cloud single-scattering parameters are interpolated from that of Fu (996) and Tsay et al. (989), which have a spectral resolution of 000 cm in the UV spectral region and 50 cm in the near IR. Table shows the optimal values of h for the four spectral bands, separately for ice and water clouds. The weight

6 FEBRUARY 998 CHOU ET AL. 207 TABLE 3. Coefficients b 0, b, and b 2 of the parameterization, (0), for the cloud single-scattering coalbedo. The units of b and b 2 are, respectively, m and m 2. Spectral band Ice cloud Water cloud b 0 b b 2 b 0 b b h is for weak absorption bands and 0 for strong absorption bands. Figures 3 5 show the results of regression using (9) (). The coefficients a, b, and c are shown in Tables 2 4, respectively. It can be seen in Fig. 3 that the extinction coefficient varies weakly with spectral band but strongly with particle size. For large particles, it is independent of wavelength for both water and ice clouds. The results shown in Fig. 3 for ice clouds are indistinguishable among the four bands. Due to a large particle size, the extinction coefficient of ice clouds is significantly smaller than that of water clouds. It is interesting to note that in spite of a large difference in the particle size distribution, the ice particle extinction coefficients for the 28 clouds nearly fall onto a single curve of (9). Figure 4 shows the single-scattering coalbedo of Bands 2 4. The single-scattering coalbedo of Band is very small (0 5 ) and is not shown in the figure. For a given particle size, the single-scattering coalbedo varies by three orders of magnitude among the three IR bands. For a given spectral band, it is 8 times larger for ice clouds than for water clouds. The asymmetry factor, shown in Fig. 5, varies between 0.78 and 0.94 for different bands and particle sizes. Since the shape of size distribution has little effect on the single-scattering properties (Hu and Stamnes 993; Slingo and Schrecker 982), the results shown in Figs. 3 5 should be representative for a wide range of clouds with various particle-size distributions. Using the radiation model discussed in section 2, the SW heating of the cloud and the atmosphere computed with the parameterization is compared to high spectralresolution calculations. The temperature and humidity profiles used are typical of a midlatitude summer (McClatchey et al. 972). The surface albedo is set to 0.2. In addition to clouds, the absorption of SW radiation due to water vapor, O 3,CO 2, and O 2 are also included. The sky is set to be overcast (i.e., fractional cover ) with either a high ( mb), middle ( mb), or low ( mb) cloud. The high cloud is assumed to contain ice particles, and the others are water clouds. For the high cloud, two samples of the ice particle size distribution measured during CEPEX and FIRE field campaigns are used. The equivalent particle size, r i,is45m for the CEPEX case and 97 m for the FIRE case. For water clouds, two log-normal size distributions are used. The equivalent radii of the droplets, r w, are 8 and 6 m. Fluxes are calculated for a range of 0..0 in the cosine of the solar zenith angle, o, and in the cloud optical thickness in Band,. For other spectral bands, i, the optical thickness is scaled from i ( i / ), where i is the cloud extinction coefficient in band i. The percentage difference in the total (0.8 0 m) net downward flux at the TOA between the parameterization and the high spectral-resolution calculations is shown in Fig. 6 for the high, middle, and low clouds. It can be seen in the figure that the parameterization introduces a very small error (%) in all cases. Figure 7 shows the percentage error in the SW heating of the atmosphere. For the middle and low clouds, shown in the middle and bottom panels, the error is small (.5%). It increases to 4% for an optically thick cirrus cloud (upper panels). This relatively large error is related to the height of the clouds. Compared to water clouds, the absorption by water vapor above the cirrus cloud top is weak, which makes the heating of the atmosphere more sensitive to the parameterization for cirrus optical properties than the parameterization for water clouds. In addition to the heating of the entire atmospheric column, it is also important to understand the effect of the parameterization on the heating of the cloud layers. Figure 8 shows that the error in the heating of cirrus clouds reaches 4% (or 5.2 W m 2 ) for r i 45 m, TABLE 4. Coefficients c 0, c, and c 2 of the parameterization, (), for the cloud asymmetry factor. The units of c and c 2 are, respectively, m and m 2. Spectral band Ice cloud Water cloud c 0 c c 2 c 0 c c

7 208 JOURNAL OF CLIMATE FIG. 4. The single-scattering coalbedo of ice and water clouds as a function of the effective particle size. The values of r i and r w are, respectively, the effective size of ice and water particles. The data points are derived from (4) (6), whereas the curves are the regression using (0)..0, and o 0.5. This reduced cloud heating causes an increased heating of.6 W m 2 in the atmosphere, both above and below the cloud, leading to a reduced SW heating of 3.6 W m 2 for the total atmospheric column. This reduced heating is equivalent to 3% of the total atmospheric SW heating (including the cloud), as shown in the upper left panel of Fig. 7. The percentage error in the heating of the earth s surface is shown in Fig. 9. It is small even in the presence of high cirrus clouds. It is noted that the parameterization for g shown in the upper panel of Fig. 5 has a very small error at r i 45 and 97 m but has a relatively large error () at r i 80 m. When r i is set to 80 m, the results are similar to that shown in the upper panels of Figs Cloud overlapping Clouds could occur at various heights with fractional cover. Nearly all radiative transfer algorithms used in FIG. 5. The asymmetry factor of ice and water clouds as a function of the effective particle size. The values of r i and r w are, respectively, the effective size of ice and water particles. The data points are derived from (8), whereas the curves are the regression using (). atmospheric models apply only to a plane-parallel (i.e., horizontally homogeneous) atmosphere. Horizontal inhomogeneity is not allowed. A straightforward approach to dealing with a partial cloudiness situation is to divide the sky into sections. Within each section, an atmospheric layer is either free of clouds or filled totally with a homogeneous cloud. Radiative fluxes are then computed for each section, and the total SW heating is the sum of all sections weighted by the fractional cover of individual sections. Depending upon the number of cloud layers and the way these clouds overlap, computational costs of this approach could be huge. To simplify the computations, a cloud that partially fills a layer is usually smeared over the entire layer, and the optical thickness is adjusted by a factor dependent on the fractional cloud cover. As demonstrated by a number of studies (e.g., Harshvardhan and Randall 985; Cahalan et al. 994), the effect of on radiation is highly nonlinear, and it is not appropriate to scale (or cloud water

8 FEBRUARY 998 CHOU ET AL. 209 FIG. 6. Percentage difference in the net downward solar (0.8 0 m) flux at the TOA between the parameterization using (9) () and the high spectral-resolution (00 cm ) calculations for high, middle, and low clouds. The values of r i and r w are, respectively, the effective size of ice and water particles. The value of o is the cosine of the solar zenith angle, and is the optical thickness in Band. Positive values are shaded. path) linearly by cloud cover. In the National Center for Atmospheric Research (NCAR) community climate model CCM2, the optical thickness is scaled by f.5, where f is the fractional cloud cover (Kiehl et al. 994). In the GEOS GCM, the optical thickness was previously scaled by inverting a simple reflection function for clouds when illuminated by diffuse radiation and assuming no absorption by clouds (Sud et al. 993). As can be expected, the effective optical thickness corresponding to a cloud smeared to cover the entire sky is a function of, f,and the solar zenith angle, as well as whether scaling is based on reflection or absorption. Scaling of is further complicated by a large range of the way clouds in various layers overlap. FIG. 7. Same as Fig. 6 except for the absorption in the atmosphere (including clouds). a. Maximum-random overlapping Treatments of cloud overlapping in climate studies differ from model to model. In the GEOS GCM (Schubert et al. 993), a maximum-random overlapping scheme, shown in Fig. 0, is applied. Clouds are identified as high, middle, and low separated roughly by the 400- and 700-mb levels. Cloud layers close to each other are likely to be related. Therefore, cloud layers in each of the three height groups are assumed to be maximally overlapped, whereas those among different groups are assumed to be randomly overlapped (left panel of Fig. 0). Scaling of is applied only to the maximally overlapped clouds within each of the three height groups but not to the randomly overlapped clouds among different height groups (right panel of Fig. 0). Within each height group, clouds are smeared over the extent of the maximum cloud amount, f m, of that height group by scaling (right panel of Fig. 0). By assuming ran-

9 20 JOURNAL OF CLIMATE FIG. 8. Same as Fig. 6 except for the absorption in clouds. FIG. 9. Same as Fig. 6 except for the absorption at the surface. dom overlapping of clouds among different height groups, the atmosphere is then divided into 2 n sections, where n 3 is the number of height groups containing clouds (right panel of Fig. 0). Within each section, a layer is either totally cloud filled or cloud free. Fluxes are first computed for each section and then summed over all sections weighted by the fractional cover of individual sections. b. Scaling of optical thickness Scaling of cloud optical thickness in the GEOS GCM is based on the maximum-random cloud overlapping assumption shown in Fig. 0. To simplify the scaling of, we assume that reflection is more important to climate than absorption, since the former affects the heating of the earth atmosphere system, but the latter primarily redistributes the heating within the atmosphere as well as between the atmosphere and the surface. As stated in section 2, our radiation algorithm requires calculations of reflection/transmission of a cloud layer, separately for direct and diffuse radiation. By smearing cloud cover over the extent of the maximum cloud cover of respective height group, the optical thickness of a cloud layer is scaled by factors s and s, b s for direct radiation and (5) f s for diffuse radiation, (6) so that albedos of the layer remain the same before and after scaling, R( b, r e, o) f R(, r e, o) for direct radiation (7) and R ( f, r e) f R (, r e) for diffuse radiation, (8) where o is the cosine of the solar zenith angle, R is the albedo averaged over the entire solar spectrum for

10 FEBRUARY 998 CHOU ET AL. 2 FIG. 0. The maximum-random cloud-overlapping scheme. Clouds are identified as high, middle, and low separated roughly by the 400- and 700-mb levels. Clouds are assumed maximally overlapped within each height group and randomly overlapped among different height groups (left panel). Scaling of the optical thickness applies only to the maximally overlapped cloud layers (right panel). direct incident radiation, R is the albedo for diffuse incident radiation given by e 0 R (, r ) 2 R(, r, ) d, (9) e f is the normalized cloud cover given by f f/ f m, (20) and is the cosine of the zenith angle of an incident beam. Using the scaling of (5) and (6), a layer with a fractional cloud cover, f, and an optical thickness,, is reduced to a layer with a fractional cloud cover, f m, and equivalent optical thicknesses b for direct radiation and f for diffuse radiation. It is noted that the scaling is a function of the incident angle, and the optical thickness is scaled separately for the direct radiation and the diffuse radiation. Using values of r e representative of the water and ice cloud samples used in section 3, and g are computed from (0) and () for the individual bands given in Table. Cloud albedos R and R averaged over the solar spectrum are computed using the radiation model discussed in section 2. The scaling functions s and s are then derived from (5) (8) by giving f, R(, r e, o ), and R(, r e ). Sample calculations show that s and s depend only weakly on r e and the phase of cloud particles. Therefore, they are reduced to s ( f,, o ) and s ( f, ) by taking the mean values of the scaling functions for various particle sizes and phases. The above scaling ensures that the reflection of a cloud layer is the same before and after the scaling. However, errors are expected in the calculations of atmospheric and surface SW heating when there are other cloud layers present. In addition to the solar zenith angle, cloud cover, and optical thickness, an optimal cloud scaling should, in principle, also depend on the fractional cover and optical thickness of other cloud layers. This situation makes scaling of the optical thickness extremely difficult. To simplify the problem, we made simple adjustments to the scaling functions s and s, so that errors in the fluxes at the TOA and the surface are minimized for the simplest cloud situation, which has two cloud layers in only one of the three height groups. For the two-cloud situation, scaling is applied to the cloud with a smaller cloud cover (see right panel of Fig. 0). It is noted that scaling is not required when there is only one cloud layer in a given height group. The scaling functions are adjusted according to s ( f s) and (2) s ( f s), (22) where and are constants empirically (trial and error) chosen to minimize errors in the fluxes at the TOA and at the surface. The errors are defined as the differences between the flux computed with the scaling of and the flux computed with the sky divided into horizontally homogeneous sections (three sections in the two-cloud situation). The adjusted scaling functions preserve the

11 22 JOURNAL OF CLIMATE FIG.. Scaling functions for the cloud optical thickness, and, for direct (shown only for o 0.5) and diffuse radiation, respectively. The value of f is the normalized cloud cover defined by (20). shapes of s and s. The results for the scaling functions are saved as tables for later use. It is found that optimal values of and are 0.2 and 0.3, respectively. The results for the scaling functions and are shown in Fig.. It can be seen that ( f,, o ) f, and f as 0. For an optically thick cloud, and are not sensitive to. The error induced in flux calculations due to the scaling of is calculated for the case of maximum overlapping of two contiguous cloud layers within a given height group. The scaling uses (5) and (6), except s and s are replaced by and. The percentage errors in the net downward flux at the TOA and at the surface are shown in Figs. 2 4 for pairs of cloud layers in the high, middle, and low height groups, respectively. The midlatitude summer atmosphere of McClatchey et al. (972) is used. The high clouds are set to be in the mb and mb layers and contain ice particles with r i 97 m and 2 for each cloud layer, where is the optical thickness in Band. The FIG. 2. Percentage error in the net downward solar (0.8 0 m) flux at the TOA and at the surface induced by the scaling of the optical thickness for a pair of clouds ( mb and mb) in the high-cloud group. The value of f u is the fractional cover of the upper-cloud layer, and f i is the fractional cover of the lowercloud layer. The optical thickness in Band,, of each cloud layer is set to 2. The value of o is the solar zenith angle. Positive values are given by solid curves, and negative values by dotted curves. middle clouds are set to be in the mb and mb layers and contain water droplets with r w 0 m and 5 for each cloud layer. The low clouds are set to be in the mb and mb layers and contain water droplets with r w 0 m and 0 for each cloud layer. The single-scattering albedo and asymmetry factor are derived from (0) and () for the four bands shown in Table. The optical thickness in Bands 2 4 are scaled from that in Band ( ) by the ratio i / derived using (9), where the subscripts denote spectral bands, and i 2, 3, or 4. The results shown in the figures are for three solar zenith angles, 30, 60, and 75. The ordinate is the fractional cover of the upper cloud layer, f u, and the abscissa is the fractional cover of the lower cloud layer, f l. It can be seen in Figs. 2 4 that, generally, the per-

12 FEBRUARY 998 CHOU ET AL. 23 FIG. 3. Same as Fig. 2 except for a pair of clouds ( mb and mb) in the middle-cloud group and with 5 for each layer. FIG. 4. Same as Fig. 2 except for a pair of clouds ( mb and mb) in the low-cloud group and with 0 for each layer. centage error increases with increasing solar zenith angle. The errors are smaller when the upper-cloud amount is larger than the lower-cloud amount. In this situation, scaling of the optical thickness is applied to the lowercloud layer. When the cloud amounts of the two layers are equal, that is, f m f u f l, there is no need to scale the optical thickness, and there are no errors, as shown by the diagonal lines in the figures. Except for cases with large solar zenith angles, the error is less than a few percent. In all cases, the absolute errors are small both at the TOA and at the surface. In the GEOS GCM, the scaling functions derived for the two-cloud situation is applied universally to all situations even if there are more than two cloud layers in a given height group. It implicitly assumes that, in a given height group, accurate scaling of is not critical if there are two cloud layers above. 5. Conclusions We have developed parameterizations for computing cloud single-scattering properties and for scaling the optical thickness in a partial cloudiness situation. Due primarily to large differences in particle size and shape, differences in the single-scattering properties of ice and water clouds are large. Therefore, the extinction coefficient, single-scattering albedo, and asymmetry factor are parameterized separately for ice and water clouds as functions of the effective mean particle size. These parameterizations are variations from that of Slingo (989) for water clouds and that of Fu (996) for ice clouds. The shortwave spectrum is divided into four broad bands nearly identical to those used by Slingo (989), and the parameterizations for ice clouds are based on the high spectral-resolution single-scattering data from Fu (996). The single scattering coalbedo varies greatly with wavelength, and the approach to averaging the single-scattering coalbedo over a broad spectral band is empirically determined so that errors in flux calculations are minimized. Depending upon the strength of absorption, the averaging ranges between linear and logarithmic.

13 24 JOURNAL OF CLIMATE Compared to high spectral-resolution calculations, the parameterization for the cloud single-scattering properties introduces a small relative error of % at TOA. The error exceeds 0% in the SW heating by ice clouds. Due to a substantial compensation in the absorption within clouds and the regions above and below clouds, the relative errors in the absorption of SW radiation in the atmosphere and at the surface are less than a few percent. Multiple-scattering radiative transfer models used in atmospheric models are developed for plane-parallel atmospheres, and scaling of cloud optical thickness in an atmosphere with multiple cloud layers is necessary to make computations affordable. The scaling of optical thickness depends on the assumption applied to the overlapping of clouds at various heights. In the GEOS GCM, a maximum-random assumption for cloud overlapping is applied. Clouds are identified as high, middle, and low, separated approximately by the 400- and 700- mb levels. Clouds are assumed to be maximally overlapped within each height group but randomly overlapped among the three height groups. Scaling of the optical thickness is separately applied to each height group. Fluxes are then computed by dividing the atmosphere into a maximum of eight horizontally homogeneous sections. The fractional cloud amount of a layer in each section equals either 0 or. Scaling of the optical thickness is based on the case with two maximally overlapped cloud layers. Thus, it implies that flux calculations are not sensitive to the scaling of the optical thickness of a cloud layer with more than two cloud layers lying above. The scaling functions are calculated using a radiative transfer model and then empirically adjusted so that errors in the fluxes at the TOA and at the surface introduced by the scaling are minimized. They are separately applied to direct and diffuse radiation. Tables for these scaling functions were precomputed as functions of the cloud optical thickness, cloud amount, and, for the direct radiation case, the solar zenith angle. Flux errors introduced by the scaling are small, both at the TOA and at the surface. Except for cases with large solar zenith angles, the error is only a few percent. In terms of absolute error, it is within a few watts per meter squared. Acknowledgments. This work was supported by the Global Atmospheric Modeling and Analysis Program, Office of Mission to Planet Earth, NASA. The authors are grateful to Qiang Fu of Dalhousie University and Si-Chee Tsay of NASA/Goddard Space Flight Center for providing high spectral-resolution single-scattering data for ice clouds and water clouds. The clarity of the paper has been significantly improved due to suggestions made by the editor, James A. Coakley, Jr., and two reviewers. REFERENCES Cahalan, R. F., W. Ridgway, W. J. Wiscombe, and T. L. Bell, 994: The albedo of stratocumulus clouds. J. Atmos. Sci., 5, Chou, M.-D., 990: Parameterization for the absorption of solar radiation by O 2 and CO 2 with application to climate studies. J. Climate, 3, , 992: A solar radiation model for use in climate studies. J. Atmos. Sci., 49, , and K.-T. Lee, 996: Parameterizations for the absorption of solar radiation by water vapor and ozone. J. Atmos. Sci., 53, , and L. Kouvaris, 99: Calculations of transmission functions in the IR CO 2 and O 2 bands. J. Geophys. Res., 96, Espinoza, R. C., Jr., and Harshvardhan, 996: Parameterization of solar near-infrared radiative properties of cloudy layers. J. Atmos. Sci., 53, Fu, Q., 996: An accurate parameterization of the solar radiative properties of cirrus clouds for climate models. J. Climate, 9, Harshvardhan, and D. A. Randall, 985: Comments on The parameterization of radiation for numerical weather prediction and climate models. Mon. Wea. Rev., 3, Hu, Y. X., and K. Stamnes, 993: An accurate parameterization of the radiative properties of water clouds suitable for use in climate models. J. Climate, 6, Joseph, J. H., W. J. Wiscombe, and J. A. Weinman, 976: The delta- Eddington approximation for radiative flux transfer. J. Atmos. Sci., 33, Kiehl, J. T., J. J. Hack, and B. P. Briegleb, 994: The simulated earth radiation budget of the National Center for Atmospheric Research community climate model CCM2 and comparisons with the Earth Radiation Budget Experiment (ERBE). J. Geophys. Res., 99, Lacis, A. A., and J. E. Hansen, 974: A parameterization for the absorption of solar radiation in the earth s atmosphere. J. Atmos. Sci., 3, McClatchey, R. A., R. W. Fenn, J. E. A. Selby, F. E. Volz, and J. S. Garing, 972: Optical properties of the atmosphere. Tech. Note AFCRL , 08 pp. [NTIS N73842.] Ramaswamy, V., and S. M. Freidenreich, 992: A study of broadband parameterizations of the solar radiative interactions with water vapor and water drops. J. Geophys. Res., 97, Sagan, C., and J. B. Pollack, 967: Anisotropic nonconservative scattering and the clouds of Venus. J. Geophys. Res., 72, Schubert, S. D., R. B. Rood, and J. Pfaendtner, 993: An assimilated dataset for earth science applications. Bull. Amer. Meteor. Soc., 74, Slingo, A., 989: A GCM parameterization for the shortwave radiative properties of water clouds. J. Atmos. Sci., 46, , and H. M. Schrecker, 982: On the shortwave radiative properties of stratiform water clouds. Quart. J. Roy. Meteor. Soc., 08, Sud, Y. C., W. C. Chao, and G. K. Walker, 993: Dependence of rainfall on vegetation: Theoretical consideration, simulation experiments, observations, and inferences from simulated atmospheric soundings. J. Arid Environ., 25, 5 8. Tao, W.-K., S. Long, J. Simpson, C.-H. Sui, B. Ferrier, and M.-D. Chou, 996: Cloud-radiation mechanisms associated with a tropical and a midlatitude squall line. J. Atmos. Sci., 53, Tsay, S. C., K. Stamnes, and K. Jayaweera, 989: Radiative energy balance in the cloudy and hazy Arctic. J. Atmos. Sci., 46, Yang, P., and K.-N. Liou, 995: Light scattering by hexagonal ice crystals: Comparison of FDTD and geometric optics models. J. Opt. Soc. Amer., 2,