Optical Properties of Shallow Convective Clouds Diagnosed from a Bulk-Microphysics Large-Eddy Simulation

Joanna Slawinska Institute of Geophysics, University of Warsaw, Warsaw, Poland

Search for other papers by Joanna Slawinska in
Current site
Google Scholar
PubMed
Close
,
Wojciech W. Grabowski National Center for Atmospheric Research,* Boulder, Colorado

Search for other papers by Wojciech W. Grabowski in
Current site
Google Scholar
PubMed
Close
,
Hanna Pawlowska Institute of Geophysics, University of Warsaw, Warsaw, Poland

Search for other papers by Hanna Pawlowska in
Current site
Google Scholar
PubMed
Close
, and
Andrzej A. Wyszogrodzki National Center for Atmospheric Research,* Boulder, Colorado

Search for other papers by Andrzej A. Wyszogrodzki in
Current site
Google Scholar
PubMed
Close
Full access

* The National Center for Atmospheric Research is sponsored by the National Science Foundation

Corresponding author address: Wojciech W. Grabowski, NCAR/MMM, P.O. Box 3000, Boulder, CO 80307-3000. Email: grabow@ncar.ucar.edu

* The National Center for Atmospheric Research is sponsored by the National Science Foundation

Corresponding author address: Wojciech W. Grabowski, NCAR/MMM, P.O. Box 3000, Boulder, CO 80307-3000. Email: grabow@ncar.ucar.edu

1. Introduction

Large-eddy simulation (LES) models provide an indispensable tool to study processes within cloud- topped subtropical and trade wind boundary layers (e.g., Siebesma et al. 2003; Stevens et al. 2005, and references therein). Typically, LES models use bulk representation of cloud microphysics. Bulk approach assumes that warm (ice free) clouds are exactly at water saturation and that a cloud cannot exist in undersaturated conditions. It is well established from cloud observations (e.g., Stommel 1947; Warner 1955; Blyth 1993; Wang and Albrecht 1994) and cloud modeling (e.g., Brenguier and Grabowski 1993; Carpenter et al. 1998; Siebesma et al. 2003; Chosson et al. 2007) that shallow cumulus and stratocumulus clouds are strongly diluted by entrainment and that such dilution affects not only bulk thermodynamic properties [e.g., the liquid water content (LWC)] but also cloud microphysics (i.e., the spectrum of cloud droplets). For nonprecipitating clouds, conservation of total water (vapor plus liquid) and moist static energy determines bulk thermodynamic properties of cloud volumes diluted by entrainment of environmental air. Predicting changes of the cloud droplet spectra resulting from entrainment, on the other hand, requires additional constraints because situations where cloud water after homogenization is distributed over either the same number of smaller droplets (i.e., the homogeneous mixing scenario) or smaller number of droplets with the initial size (i.e., the extremely inhomogeneous mixing scenario; Baker and Latham 1979; Baker et al. 1980) are equally possible (see discussion in Andrejczuk et al. 2006).

The impact of entrainment and mixing on cloud droplet spectra has been shown to significantly affect mean radiative properties of a field of stratocumulus and cumulus clouds. Chosson et al. (2004) were first to show this for the case of stratocumulus. They examined the impact of microphysical transformation following entrainment and mixing by considering separately the homogeneous and extremely inhomogeneous mixing scenarios and showed that the mean cloud optical thickness derived by applying the homogeneous scheme was about 35% larger than for the extremely inhomogeneous mixing (see Chosson et al. 2007 for further discussion). This result prompted an investigation reported in Grabowski (2006, hereafter G06), where a similar issue was investigated in the context of convective–radiative quasi equilibrium mimicking the mean conditions on the earth and with the emphasis on the indirect aerosol effects. Two limits of the concentration of cloud droplets no were considered, either no = 100 cm−3, referred to as PRISTINE, or no = 1000 cm−3, referred to as POLLUTED. The 2D model used in G06’s study applied bulk microphysics and, similarly to Chosson et al. (2004, 2007), required additional assumptions concerning sizes of cloud droplets when coupled to the radiative transfer model. The relevant parameter is the effective radius, the ratio between the third and the second moment of the cloud droplets size distribution (e.g., Stephens 1978). Effective radius has been shown to be slightly larger than the mean volume radius (e.g., Martin et al. 1994; Pawlowska and Brenguier 2000). In G06’s study, the formulation of the effective radius in the diluted cloudy volumes turned out to be of critical importance, with the mean top-of-the-atmosphere (TOA) albedo (and thus the amount of solar energy reaching the surface) being the same in the PRISTINE case assuming the homogeneous mixing scenario and in the POLLUTED case with the extremely inhomogeneous mixing.

As on the earth, the model cloudiness in G06 was dominated by shallow convection. However, since the model setup was designed with the deep convection in mind, the horizontal resolution was not adequate for shallow clouds. Moreover, the model applied 2D geometry, which also limits the generality of conclusions. It follows that it would be desirable to perform a similar investigation using a 3D model with a spatial resolution appropriate for shallow clouds, that is, using a LES model. This paper reports results of such an investigation. The next section briefly discusses the model and modeling setup. Section 3 provides details of the formulation of effective radius applying various assumptions concerning the local concentration of cloud droplets. Vertical profiles of the diagnosed effective radius are presented in section 4, together with results of radiation transfer model applied offline to the cloud data. Brief discussion of model results in section 5, including a comparison with observations reported in McFarlane and Grabowski (2007), concludes the paper.

2. The model and modeling setup

The model used in this study is the 3D anelastic semi-Lagrangian–Eulerian model EULAG documented in Smolarkiewicz and Margolin (1997, model dynamics), Grabowski and Smolarkiewicz (1996, model thermodynamics), and Margolin et al. (1999, subgrid-scale turbulent mixing). In this study, EULAG is set up to simulate steady-state trade wind shallow nonprecipitating convection observed during the Barbados Oceanographic and Meteorological Experiment (BOMEX; Holland and Rasmusson 1973) and recently used in the model intercomparison study described in Siebesma et al. (2003). In the BOMEX case, the 1.5-km-deep trade wind convection layer overlays 0.5-km-deep mixed layer near the ocean surface and is covered by 500-m-deep trade wind inversion layer. The cloud cover is about 10% and quasi-steady conditions are maintained by the prescribed large-scale subsidence, large-scale moisture advection, surface heat fluxes, and radiative cooling. The model setup is exactly as described in Siebesma et al. (2003). The horizontal/vertical grid length is 100 m/40 m and model time step is 3 s. As in Siebesma et al. (2003), the model is run for 6 h, and results from the last 4 h are used in the analysis.

In general, EULAG results fall within the spread of all models presented in Siebesma et al. (2003). It follows that EULAG simulation faithfully represents gross features of the trade wind convection observed in BOMEX, including mean thermodynamic and wind profiles, mean characteristics of the cloud field (e.g., the cloud cover and liquid water path), and mean profiles of the temperature and moisture fluxes between the surface and trade wind inversion. A single simulation (archived every 10 min of the model time) is used in the analysis presented in this paper. This is because, unlike in G06, effective radius and optical thickness are only diagnosed using various assumptions concerning the mixing between clouds and their environment. Furthermore, we apply a solar radiation transfer model (the same as in G06) to document the impact of the diagnosed effective radius on the mean cloud albedo and on the net solar flux at the surface.

3. Derivation of the effective radius from bulk model output

As suggested by Martin et al. [1994, Eq. (10) therein], the effective radius re is assumed proportional to the mean volume radius rυ:
i1520-0442-21-7-1639-e1
where k = 0.80 for the PRISTINE and k = 0.67 for the POLLUTED. The mean volume radius rυ, on the other hand, is locally defined as
i1520-0442-21-7-1639-e2
where ρw = 103 kg m−3 is the water density, and LWC and n are the local values of the liquid water content and concentration of cloud droplets, respectively. Note that EULAG predicts LWC, but not n, and additional assumptions are needed to predict local values of rυ and thus re. Four different assumptions concerning variability of n are reviewed below and are later used in the analysis of the model output.

a. Uniform cloud droplet concentration, r(u)υ, and r(u)e

In this case, it is assumed that the concentration of cloud droplets does not change within a cloud and it is taken in (2) as n = no (i.e., either 100 cm−3 for the PRISTINE or 1000 cm−3 for POLLUTED). This case was referred to as the homogeneous mixing scenario in G06. It was recognized, however, that this is an approximation because the homogeneous mixing involves also a change of droplet concentration, with the concentration in a diluted parcel smaller by a factor of χ, the proportion of the cloudy air in the mixture (see discussion below and in the appendix in G06).

b. Homogeneous mixing scenario, r(h)υ, and r(h)e

In this case, it is assumed that the concentration of cloud droplets changes within a cloud only because of the dilution due to entrainment. It follows that n is taken as n = χno in (2), where, again, no = 100 cm−3 for the PRISTINE and 1000 cm−3 for POLLUTED. The proportion of the cloudy air in the mixture χ is approximated, assuming that the total water at a given location is a linear combination of the total water Q inside the adiabatic parcel raising from the cloud base and an unsaturated environmental air at this level:
i1520-0442-21-7-1639-e3
where Q = qυ + qc (qυ and qc are the water vapor and cloud water mixing ratios), and superscripts a and e refer to the adiabatic and environmental values, respectively. Note that (3) represents an extreme simplification of the entrainment problem in cumulus clouds (cf. Blyth 1993). Since there is no precipitation, Qa is equal to the total water at the cloud base. The environmental total water Qe = qeυ can be approximated as the mean over cloud-free grid boxes at this level. It follows from (3) that
i1520-0442-21-7-1639-e4
where additional limiting ensures only physically realizable values of χ. Adiabatic parcel analysis is performed using cloud-base conditions obtained by conditional sampling (i.e., considering only points with qc > 10−5 kg kg−1) at a height of 600 m, which is within the range of simulated cloud-base heights for various clouds during the 4-h analysis period. In summary, local values of χ are deduced from (4) and later used in (2) with n = χno to derive r(h)υ.

c. Intermediate mixing scenario, r(in)υ, and r(in)e

In this case, it is assumed that the concentration of cloud droplets changes within a cloud because of the dilution due to entrainment and that the change in droplet concentration is the same as the change of the mean volume radius cubed. Such a mixing scenario was suggested by direct numerical simulations of interfacial mixing at low turbulence levels discussed in Andrejczuk et al. (2006). In such a mixing scenario, the mean volume radius is prescribed as (see appendix in G06)
i1520-0442-21-7-1639-e5
where LWCa and raυ are the adiabatic liquid water content and the adiabatic mean volume radius [i.e., LWCa = 4/3 π (raυ)3 noρw]. The limiting in (5) ensures only physically realizable values of β.

d. Extremely inhomogeneous mixing scenario, r(ei)υ, and r(ei)e

In this case, it is assumed that entrainment changes only the concentration of cloud droplets without affecting the mean volume radius. In this case,
i1520-0442-21-7-1639-e6

The above four formulations of rυ result in various spatial distributions of re, separately for the PRISTINE and POLLUTED cases (i.e., choosing no as either 100 or 1000 cm−3). Selected statistics of these distributions are compared in the next section. Moreover, in each cloudy column, the cloud optical thickness can be derived as τ = 1.5/ρw∫LWC/re dz, with re defined in various ways as discussed above. Finally, a radiation transfer model can be applied to cloudy columns to assess the impact of various mixing scenarios on solar radiation.

4. Results

Figure 1 presents the contoured frequency by altitude diagrams (CFADs) of re derived according to the scenarios 1 to 4 and for PRISTINE and POLLUTED conditions. As expected, mean re increases with height in all scenarios. A noisy pattern near cloud tops in mixing scenarios 1–3 most likely reflects poor sampling of the deepest clouds. For a given formulation, POLLUTED cases show significantly smaller re when compared to the PRISTINE. Since the uniform scenario does not consider dilution of droplet concentration due to entrainment, it results in the smallest droplet sizes compared to other scenarios. The mean re at a given height gradually increases when one moves across mixing scenarios from uniform droplet concentration to the extremely inhomogeneous mixing. The effective radius in the extremely inhomogeneous mixing (i.e., the adiabatic re) is larger than the maximum re predicted by any other mixing scenario at the same height. This attests to the significant dilution of shallow convective clouds simulated in this study, which is best illustrated by a CFAD of the adiabatic fraction LWC/LWCa (not shown; cf. Fig. 7 in Siebesma et al. 2003). At any level, the uniform droplet concentration results in a histogram with the largest width. The difference between uniform droplet concentration and homogeneous mixing is significant, an effect neglected in G06. Of a particular relevance is the fact that the adiabatic re for the extremely inhomogeneous mixing in POLLUTED is close to the mean re for the uniform droplet concentration in PRISTINE; that is, re in both cases is in the range between 5 and 10 μm across most of the cloud depth.

Figure 2 shows probability density functions (PDFs) of the optical thickness τ for PRISTINE and POLLUTED and for uniform droplet concentration and extremely inhomogeneous mixing (i.e., the two scenarios that provide limiting sizes of re for either PRISTINE or POLLUTED case). Only model columns with liquid water path (LWP) larger than 5 × 10−3 kg m−2 are included in the analysis. Because of the limited vertical extent of BOMEX clouds, values of the optical thickness larger than 100 are rare, unlike in the case of convective–radiative equilibrium of G06, which featured deep as well as shallow clouds. The PDFs are similar to those in G06 (Fig. 5 therein), with maximum values for τ < 10 and long tails extending to values τ ∼ 100. As in Fig. 5 of G06, changing mixing scenario from the uniform droplet concentration to the extremely inhomogeneous mixing leads to a decrease of PDF for small τs and an increase for τ > 10. However, changes in the PDFs between PRISTINE and POLLUTED cases are smaller than in G06, consistent with only the first indirect effect considered here (i.e., changes of re alone), whereas both the first and the second indirect effects were considered in G06 (i.e., changes in both re and LWP).

To document the impact of various mixing scenarios on mean optical properties of the cloud field, the two-stream radiation transfer model was applied column by column (i.e., in the independent column approximation mode) to the cloud model data combined with the diagnosed effective radius of cloud droplets (i.e., similarly to Chosson et al. 2004, 2007). The radiation transfer model, the same as in G06, comes from the National Center for Atmospheric Research (NCAR) Community Climate System Model (Kiehl et al. 1994). The solar constant is taken as 436 W m−2 (i.e., the nominal solar constant at an equinox divided by π to represent insolation averaged over the diurnal cycle) and a zero zenith angle is assumed. Figure 3 shows histograms of the TOA cloud albedo for PRISTINE and POLLUTED and for the uniform droplet concentration and the extremely inhomogeneous mixing (i.e., the two scenarios that provide limiting sizes of re for either PRISTINE or POLLUTED case). As in Fig. 2, only model columns with LWP larger than 5 × 10−3 kg m−2 are included in the analysis.1 As Fig. 3 shows, the range of TOA albedos for cloudy columns is between 0.1 and 0.7, with consistent changes for different local mixing scenarios. For instance, the change from uniform droplet concentration to extremely inhomogeneous mixing (i.e., with larger cloud droplets diagnosed in diluted cloud volumes) results in reduced frequencies for high albedos and corresponding increase of the peak for albedos around 0.1. This is true for both PRISTINE and POLLUTED cases.

Model results for various mixing scenarios are summarized in Table 1. The table shows averaged values of the mean optical thickness τ, the mean “effective” re, the mean cloud TOA albedo Acloudy, and the mean net solar flux at the surface SFcloudy, all calculated including only model columns with the LWP larger than 5 × 10−3 kg m−2, as in Figs. 2 and 3. The effective re is calculated in each column using the LWP and the optical thickness τ as re = 1.5 LWP/(ρwτ) [cf. Eq. (3) in G06]; its averaged value (for columns with the LWP larger than 5 × 10−3 kg m−2) is shown in the table. The conditionally averaged liquid water path is 8 × 10−2 kg m−2 (i.e., 3 to 4 times smaller than in the deep convection case of G06, Table 3 therein). In general, changes of the mean optical thickness and effective radius between PRISTINE and POLLUTED with various mixing scenarios in this study are similar to the changes reported in G06 (Table 3 therein). For instance, the mean optical thickness for PRISTINE with uniform droplet concentration is 11.5 and it is 15.8 for POLLUTED with extremely inhomogeneous mixing; that is, it changes by a factor of 1.4 between the two cases. Corresponding values are 32.1 and 57.4 in G06, which implies a factor of 1.8. For the mean effective re, corresponding values are 8.1 and 6.7 μm in this study and 8.3 and 7.4 μm in G06. Unlike in G06, the mean albedo for PRISTINE with uniform droplet concentration is smaller than for POLLUTED with extremely inhomogeneous mixing (0.33 versus 0.38), and the difference is comparable to the range observed for either PRISTINE or POLLUTED (from 0.27 to 0.33 or 0.38 to 0.45). For the SFcloudy, the differences between different mixing scenarios seem larger than in G06 (i.e., the change here is about 30 W m−2 from uniform droplet concentration to extremely inhomogeneous mixing compared to about 20 W m−2 in G06). One needs to keep in mind, however, that a different incoming solar flux was applied here than in G06 (436 versus 342 W m−2), and that results in Table 1 only include cloudy columns. Nevertheless, the overall impact of mixing scenarios in 3D LES simulations discussed here is consistent with results discussed in G06. In particular, POLLUTED clouds assuming extremely inhomogeneous mixing scenario have optical properties not far from PRISTINE clouds with uniform droplet concentration.

5. Discussion and conclusions

Results presented in this paper strongly support predictions of the low-resolution 2D study by G06 that the formulation of the effective radius in diluted volumes of shallow convective clouds plays an essential role in mean optical properties of a cloud field and thus in the mean TOA albedo of an ensemble of clouds. Such a conclusion is also consistent with results discussed in Chosson et al. (2004, 2007) for the case of stratocumulus. In the current investigation, similarly to G06, POLLUTED clouds assuming extremely inhomogeneous mixing scenario have optical properties similar to PRISTINE clouds with uniform droplet concentration.

Results presented here can be compared to the observational study of McFarlane and Grabowski (2007) where a CFAD of the effective radius and a PDF of the cloud optical thickness were derived from ground-based remote sensing of shallow tropical clouds over the Atmospheric Radiation Measurement Program (ARM) tropical western Pacific Nauru site using the technique described in McFarlane et al. (2002). As documented in McFarlane and Grabowski (2007, their Fig. 2), the CFAD of retrieved re has significant width, which increases with height, and tends to have a bimodal shape with wide peaks that limit the variability of re. The peak at small re represents cloud droplets with sizes not much different than those near the cloud base, suggesting a new activation of cloud droplets due to entrainment. The peak at large re, while representing much larger droplets, is still significantly smaller than the adiabatic size. CFADs derived from model results presented here have different shapes. The extremely inhomogeneous mixing scenario between the adiabatic cloud and its unsaturated environment predicts adiabatic re at a given height and is thus inconsistent with the observations. All three other mixing scenarios predict distributions of re with a significant width (as in observations) but with only one peak, somewhere between the adiabatic size and the cloud-base size. Despite these differences, the PDFs of the optical thickness are similar to the one derived from the observations (Fig. 4 in McFarlane and Grabowski 2007), with the maxima for low τs around 0.1 and PDFs decreasing to about 10−4 for τ = 100.

The above comparison suggests that the spatial variability of re in real clouds is more complicated than the variability predicted by simple models of microphysical transformations considered in this study. Variable in time and space transformations associated with the cloud–environment mixing (resulting in a superposition of CFADs shown in Fig. 1), combined with additional nucleation of cloud droplets above the cloud base either due to increasing updraft strength or due to entrainment (see discussion in McFarlane and Grabowski 2007), can all contribute to the observed variability. This subject will be further investigated using a newly developed two-moment warm rain bulk microphysics scheme (i.e., predicting both the number of cloud droplets and their mixing ratio), capable of predicting drizzle/rain formation as well as changes of the supersaturation field within shallow ice-free clouds (and thus the in-cloud nucleation), and including various scenarios of microphysical transformations due to turbulent entrainment and mixing (Morrison and Grabowski 2007, 2008). Such simulations will allow studying both the first and the second indirect effect using the BOMEX case considered here, as well as the precipitating shallow convection case based on observations during the Rain in Cumulus over the Ocean (RICO) experiment (see http://www.knmi.nl/samenw/rico/). Moreover, 3D radiative transfer should be applied to the cloud data as in Chosson et al. (2007). We will report results of these investigations in future publications.

Acknowledgments

This work was partially supported by the NOAA Grant NA05OAR4310107 (JS and WWG), by the 5.FP Project EVK2-CT2002-80010-CESSAR, KBN Grant SPUB-M, and PBS Grant ACCENT/295/2006 (JS and HP). Computer time was provided by NSF MRI Grants CNS-0421498, CNS-0420873, and CNS-0420985; NSF sponsorship of the National Center for Atmospheric Research and the University of Colorado; and a grant from the IBM Shared University Research (SUR) program. Comments on the manuscript by Hugh Morrison are acknowledged.

REFERENCES

  • Andrejczuk, M., W. W. Grabowski, S. P. Malinowski, and P. K. Smolarkiewicz, 2006: Numerical simulation of cloud–clear air interfacial mixing: Effects on cloud microphysics. J. Atmos. Sci., 63 , 32043225.

    • Search Google Scholar
    • Export Citation
  • Baker, M. B., and J. Latham, 1979: The evolution of droplet spectra and the rate of production of embryonic raindrops in small cumulus clouds. J. Atmos. Sci., 36 , 16121615.

    • Search Google Scholar
    • Export Citation
  • Baker, M. B., R. G. Corbin, and J. Latham, 1980: The influence of entrainment on the evolution of cloud droplet spectra: I. A model of inhomogeneous mixing. Quart. J. Roy. Meteor. Soc., 106 , 581598.

    • Search Google Scholar
    • Export Citation
  • Blyth, A. M., 1993: Entrainment in cumulus clouds. J. Appl. Meteor., 32 , 626641.

  • Brenguier, J-L., and W. W. Grabowski, 1993: Cumulus entrainment and cloud droplet spectra: A numerical model within a two-dimensional dynamical framework. J. Atmos. Sci., 50 , 120136.

    • Search Google Scholar
    • Export Citation
  • Carpenter Jr., R. L., K. K. Droegemeier, and A. M. Blyth, 1998: Entrainment and detrainment in numerically simulated cumulus congestus clouds. Part I: General results. J. Atmos. Sci., 55 , 34173432.

    • Search Google Scholar
    • Export Citation
  • Chosson, F., J-L. Brenguier, and M. Schröder, 2004: Radiative impact of mixing processes in boundary layer clouds. Proc. 14th Int. Conf. on Clouds and Precipitation, Bologna, Italy, International Association of Meteorology and Atmospheric Sciences, 371–374.

  • Chosson, F., J-L. Brenguier, and L. Schüller, 2007: Entrainment-mixing and radiative transfer simulation in boundary layer clouds. J. Atmos. Sci., 64 , 26702682.

    • Search Google Scholar
    • Export Citation
  • Grabowski, W. W., 2006: Indirect impact of atmospheric aerosols in idealized simulations of convective–radiative quasi equilibrium. J. Climate, 19 , 46644682.

    • Search Google Scholar
    • Export Citation
  • Grabowski, W. W., and P. K. Smolarkiewicz, 1996: On two-time-level semi-Lagrangian modeling of precipitating clouds. Mon. Wea. Rev., 124 , 487497.

    • Search Google Scholar
    • Export Citation
  • Holland, J. Z., and E. M. Rasmusson, 1973: Measurements of the atmospheric mass, energy, and momentum budgets over a 500-kilometer square of tropical ocean. Mon. Wea. Rev., 101 , 4455.

    • Search Google Scholar
    • Export Citation
  • Kiehl, J. T., J. J. Hack, and B. P. Briegleb, 1994: 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 , 2081520827.

    • Search Google Scholar
    • Export Citation
  • Margolin, L. G., P. K. Smolarkiewicz, and Z. Sorbjan, 1999: Large-eddy simulations of convective boundary layers using nonoscillatory differencing. Physica D, 133 , 390397.

    • Search Google Scholar
    • Export Citation
  • Martin, G. M., D. W. Johnson, and A. Spice, 1994: The measurement and parameterization of effective radius of droplets in warm stratocumulus clouds. J. Atmos. Sci., 51 , 18231842.

    • Search Google Scholar
    • Export Citation
  • McFarlane, S. A., and W. W. Grabowski, 2007: Optical properties of shallow tropical cumuli derived from ARM ground-based remote sensing. Geophys. Res. Lett., 34 .L06808, doi:10.1029/2006GL028767.

    • Search Google Scholar
    • Export Citation
  • McFarlane, S. A., K. F. Evans, and A. S. Ackerman, 2002: A Bayesian algorithm for the retrieval of liquid water cloud properties from microwave radiometer and millimeter radar data. J. Geophys. Res., 107 .4317, doi:10.1029/2001JD001011.

    • Search Google Scholar
    • Export Citation
  • Morrison, H., and W. W. Grabowski, 2007: Comparison of bulk and bin warm-rain microphysics models using a kinematic framework. J. Atmos. Sci., 64 , 28392861.

    • Search Google Scholar
    • Export Citation
  • Morrison, H., and W. W. Grabowski, 2008: Modeling supersaturation and subgrid-scale mixing with two-moment bulk warm microphysics. J. Atmos. Sci., 3 , 792812.

    • Search Google Scholar
    • Export Citation
  • Pawlowska, H., and J-L. Brenguier, 2000: Microphysical properties of stratocumulus clouds during ACE-2. Tellus, 52B , 868887.

  • Siebesma, A. P., and Coauthors, 2003: A large eddy simulation intercomparison study of shallow cumulus convection. J. Atmos. Sci., 60 , 12011219.

    • Search Google Scholar
    • Export Citation
  • Smolarkiewicz, P. K., and L. G. Margolin, 1997: On forward-in-time differencing for fluids: An Eulerian/semi-Lagrangian nonhydrostatic model for stratified flows. Atmos.–Ocean, 35 , 127152.

    • Search Google Scholar
    • Export Citation
  • Stephens, G. L., 1978: Radiation profiles in extended water clouds. II: Parameterization schemes. J. Atmos. Sci., 35 , 21232132.

  • Stevens, B., and Coauthors, 2005: Evaluation of large-eddy simulations via observations of nocturnal marine stratocumulus. Mon. Wea. Rev., 133 , 14431462.

    • Search Google Scholar
    • Export Citation
  • Stommel, H., 1947: Entrainment of air into a cumulus cloud. J. Meteor., 4 , 9194.

  • Wang, Q., and B. A. Albrecht, 1994: Observations of cloud-top entrainment in marine stratocumulus clouds. J. Atmos. Sci., 51 , 15301547.

    • Search Google Scholar
    • Export Citation
  • Warner, J., 1955: The water content of cumuliform cloud. Tellus, 7 , 449457.

Fig. 1.
Fig. 1.

CFADs of the diagnosed effective radius assuming (left) PRISTINE and (right) POLLUTED conditions. Panels represent (from top to bottom) uniform droplet concentration (u), homogeneous (h), intermediate (in), and extremely inhomogeneous (ei) mixing scenarios. Bin size is 0.5 μm. Frequency is contoured starting at 1% with contour interval of 2% and 4% for PRISTINE and POLLUTED, respectively. Note that for (ei) all contours collapse into a single profile.

Citation: Journal of Climate 21, 7; 10.1175/2007JCLI1820.1

Fig. 2.
Fig. 2.

PDFs of the optical thickness derived from model columns with LWP larger than 5 × 10−3 kg m−2 and for (left) PRISTINE and (right) POLLUTED conditions. Only data for uniform droplet concentration (u) and extremely inhomogeneous mixing (ei) are shown.

Citation: Journal of Climate 21, 7; 10.1175/2007JCLI1820.1

Fig. 3.
Fig. 3.

Histograms of cloud TOA albedo for model columns with LWP larger than 5 × 10−3 kg m−2 for (left) PRISTINE and (right) POLLUTED conditions. (top) Uniform droplet concentration and (bottom) extremely inhomogeneous mixing. The width of the albedo bin is 0.02. Mean albedo is shown in each panel.

Citation: Journal of Climate 21, 7; 10.1175/2007JCLI1820.1

Table 1.

Mean values of the optical thickness τ, re, TOA albedo Acloudy, and net solar flux at the surface SFcloudy for various mixing scenarios. Only model columns with LWP larger than 5 × 10−3 kg m−2 are included in the analysis. See text for details.

Table 1.

1

If all model columns are considered, the histogram is dominated by clear-sky values because the cloud cover in BOMEX is relatively low (around 10%; cf. Fig. 2a in Siebesma et al. 2003). The mean TOA albedo of the scene can be calculated as (1 − δ)Aclear + δAcloudy, where δ is the cloud cover and Aclear (≈0.05, mostly due to the assumed albedo of the ocean surface) and Acloudy are the mean albedos of clear-sky and cloudy columns, respectively.

Save
  • Andrejczuk, M., W. W. Grabowski, S. P. Malinowski, and P. K. Smolarkiewicz, 2006: Numerical simulation of cloud–clear air interfacial mixing: Effects on cloud microphysics. J. Atmos. Sci., 63 , 32043225.

    • Search Google Scholar
    • Export Citation
  • Baker, M. B., and J. Latham, 1979: The evolution of droplet spectra and the rate of production of embryonic raindrops in small cumulus clouds. J. Atmos. Sci., 36 , 16121615.

    • Search Google Scholar
    • Export Citation
  • Baker, M. B., R. G. Corbin, and J. Latham, 1980: The influence of entrainment on the evolution of cloud droplet spectra: I. A model of inhomogeneous mixing. Quart. J. Roy. Meteor. Soc., 106 , 581598.

    • Search Google Scholar
    • Export Citation
  • Blyth, A. M., 1993: Entrainment in cumulus clouds. J. Appl. Meteor., 32 , 626641.

  • Brenguier, J-L., and W. W. Grabowski, 1993: Cumulus entrainment and cloud droplet spectra: A numerical model within a two-dimensional dynamical framework. J. Atmos. Sci., 50 , 120136.

    • Search Google Scholar
    • Export Citation
  • Carpenter Jr., R. L., K. K. Droegemeier, and A. M. Blyth, 1998: Entrainment and detrainment in numerically simulated cumulus congestus clouds. Part I: General results. J. Atmos. Sci., 55 , 34173432.

    • Search Google Scholar
    • Export Citation
  • Chosson, F., J-L. Brenguier, and M. Schröder, 2004: Radiative impact of mixing processes in boundary layer clouds. Proc. 14th Int. Conf. on Clouds and Precipitation, Bologna, Italy, International Association of Meteorology and Atmospheric Sciences, 371–374.

  • Chosson, F., J-L. Brenguier, and L. Schüller, 2007: Entrainment-mixing and radiative transfer simulation in boundary layer clouds. J. Atmos. Sci., 64 , 26702682.

    • Search Google Scholar
    • Export Citation
  • Grabowski, W. W., 2006: Indirect impact of atmospheric aerosols in idealized simulations of convective–radiative quasi equilibrium. J. Climate, 19 , 46644682.

    • Search Google Scholar
    • Export Citation
  • Grabowski, W. W., and P. K. Smolarkiewicz, 1996: On two-time-level semi-Lagrangian modeling of precipitating clouds. Mon. Wea. Rev., 124 , 487497.

    • Search Google Scholar
    • Export Citation
  • Holland, J. Z., and E. M. Rasmusson, 1973: Measurements of the atmospheric mass, energy, and momentum budgets over a 500-kilometer square of tropical ocean. Mon. Wea. Rev., 101 , 4455.

    • Search Google Scholar
    • Export Citation
  • Kiehl, J. T., J. J. Hack, and B. P. Briegleb, 1994: 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 , 2081520827.

    • Search Google Scholar
    • Export Citation
  • Margolin, L. G., P. K. Smolarkiewicz, and Z. Sorbjan, 1999: Large-eddy simulations of convective boundary layers using nonoscillatory differencing. Physica D, 133 , 390397.

    • Search Google Scholar
    • Export Citation
  • Martin, G. M., D. W. Johnson, and A. Spice, 1994: The measurement and parameterization of effective radius of droplets in warm stratocumulus clouds. J. Atmos. Sci., 51 , 18231842.

    • Search Google Scholar
    • Export Citation
  • McFarlane, S. A., and W. W. Grabowski, 2007: Optical properties of shallow tropical cumuli derived from ARM ground-based remote sensing. Geophys. Res. Lett., 34 .L06808, doi:10.1029/2006GL028767.

    • Search Google Scholar
    • Export Citation
  • McFarlane, S. A., K. F. Evans, and A. S. Ackerman, 2002: A Bayesian algorithm for the retrieval of liquid water cloud properties from microwave radiometer and millimeter radar data. J. Geophys. Res., 107 .4317, doi:10.1029/2001JD001011.

    • Search Google Scholar
    • Export Citation
  • Morrison, H., and W. W. Grabowski, 2007: Comparison of bulk and bin warm-rain microphysics models using a kinematic framework. J. Atmos. Sci., 64 , 28392861.

    • Search Google Scholar
    • Export Citation
  • Morrison, H., and W. W. Grabowski, 2008: Modeling supersaturation and subgrid-scale mixing with two-moment bulk warm microphysics. J. Atmos. Sci., 3 , 792812.

    • Search Google Scholar
    • Export Citation
  • Pawlowska, H., and J-L. Brenguier, 2000: Microphysical properties of stratocumulus clouds during ACE-2. Tellus, 52B , 868887.

  • Siebesma, A. P., and Coauthors, 2003: A large eddy simulation intercomparison study of shallow cumulus convection. J. Atmos. Sci., 60 , 12011219.

    • Search Google Scholar
    • Export Citation
  • Smolarkiewicz, P. K., and L. G. Margolin, 1997: On forward-in-time differencing for fluids: An Eulerian/semi-Lagrangian nonhydrostatic model for stratified flows. Atmos.–Ocean, 35 , 127152.

    • Search Google Scholar
    • Export Citation
  • Stephens, G. L., 1978: Radiation profiles in extended water clouds. II: Parameterization schemes. J. Atmos. Sci., 35 , 21232132.

  • Stevens, B., and Coauthors, 2005: Evaluation of large-eddy simulations via observations of nocturnal marine stratocumulus. Mon. Wea. Rev., 133 , 14431462.

    • Search Google Scholar
    • Export Citation
  • Stommel, H., 1947: Entrainment of air into a cumulus cloud. J. Meteor., 4 , 9194.

  • Wang, Q., and B. A. Albrecht, 1994: Observations of cloud-top entrainment in marine stratocumulus clouds. J. Atmos. Sci., 51 , 15301547.

    • Search Google Scholar
    • Export Citation
  • Warner, J., 1955: The water content of cumuliform cloud. Tellus, 7 , 449457.

  • Fig. 1.

    CFADs of the diagnosed effective radius assuming (left) PRISTINE and (right) POLLUTED conditions. Panels represent (from top to bottom) uniform droplet concentration (u), homogeneous (h), intermediate (in), and extremely inhomogeneous (ei) mixing scenarios. Bin size is 0.5 μm. Frequency is contoured starting at 1% with contour interval of 2% and 4% for PRISTINE and POLLUTED, respectively. Note that for (ei) all contours collapse into a single profile.

  • Fig. 2.

    PDFs of the optical thickness derived from model columns with LWP larger than 5 × 10−3 kg m−2 and for (left) PRISTINE and (right) POLLUTED conditions. Only data for uniform droplet concentration (u) and extremely inhomogeneous mixing (ei) are shown.

  • Fig. 3.

    Histograms of cloud TOA albedo for model columns with LWP larger than 5 × 10−3 kg m−2 for (left) PRISTINE and (right) POLLUTED conditions. (top) Uniform droplet concentration and (bottom) extremely inhomogeneous mixing. The width of the albedo bin is 0.02. Mean albedo is shown in each panel.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 194 88 10
PDF Downloads 77 25 1