## 1. Introduction

Marine stratus and stratocumulus clouds play an important role in global water and energy budget by increasing the local planetary albedo by 30%–50% while having relatively little impact on outgoing longwave radiation. These clouds are widespread and at any given time may cover much of the eastern subtropical Pacific and Atlantic, most of the Arctic Sea in summer, and large regions of the middle latitudes. Charlson et al. (1987) estimate that stratus and stratocumulus clouds are the only cloud type that cover over 25% of the world’s oceans. It is estimated (Ramanathan et al. 1989) that an increase of a few percent of cloud cover, or a comparable increase in stratocumulus cloud albedo, would counter the anticipated greenhouse warming of the next century, while similar decreases would double the warming.

Because of the large spatial coverage and persistence of marine stratocumulus, drizzle exerts a powerful influence on the structure and longevity of stratocumulus. The primary effects of drizzle are on vertical redistribution of heat and moisture (Paluch and Lenschow 1991; Feingold et al. 1996). This redistribution then influences microphysical properties and cloud macrophysical characteristics (cloud fraction, liquid water content), which can lead to cloud system breakup (Stevens et al. 1998) and significant changes in the surface radiation and moisture budgets. The latter effects were illustrated by vanZanten and Stevens (2005) who used in situ aircraft data from the Second Dynamics and Chemistry of Marine Stratocumulus (DYCOMS-II) field project to demonstrate that drizzle can be an important term in the water budget for marine stratocumulus. Austin et al. (1995) previously showed that local rain rates can be a significant fraction of the surface hydrologic balance—as much as 4–5 times larger than the local surface moisture flux. Drizzle is also an important component of the hydrologic cycle on climate scales, and an intercomparison of global climate models emphasized a systematic underestimate of the frequency of drizzle events over stratocumulus cloud regimes (Lau et al. 1996).

*Q*) and precipitation flux (

_{l}*R*) based on three different sets of parameters: (a) radar reflectivity,

*Z*, (b) radar reflectivity and Doppler velocity,

*V*, and (c) radar reflectivity and Doppler velocity spectrum width,

_{d}*σ*. As radar reflectivity represents the sixth moment of the drop size distribution (DSD), one can expect it to be correlated with other moments of the DSD, such as liquid water content

_{d}*Q*(third moment of DSD), or drizzle flux

_{l}*R*, which in stratocumulus clouds is proportional to the fourth DSD moment. Thus, a number of studies have been devoted to retrievals of

*Q*and

_{l}*R*in boundary layer stratocumulus based on radar reflectivity

*Z*alone. The success of the

*Q*retrievals depended on cloud type, but even more on the absence of drizzle, both in the cloud and below cloud base. The retrieval of

_{l}*Q*is rather straightforward in nondrizzling stratocumulus where cloud spectra are mostly unimodal and the contribution to reflectivity from the large droplet tail of the spectrum is minimal. A simple

_{l}*Z*–

*Q*relation in this case is justified (Sauvageot and Omar 1987; Frisch et al. 1995; Fox and Illingworth 1997):

*a*and

*b*depend on assumptions about the drop number concentration and the shape (mostly the width) of the drop spectrum. The task becomes more complicated once drizzle drops are present in significant numbers. Drizzle typically contributes little to

*Q*yet can profoundly influence reflectivity, which is proportional to the sixth moment of the droplet size distribution and thus sensitive to the large drop tail of the DSD. For this reason, radar reflectivity alone may not be sufficient for an accurate retrieval, especially in drizzling cases where a significant fraction of cloud liquid water is carried by small drops (

_{l}*r*< 25

*μ*m). To enhance the accuracy of

*Q*retrievals a number of studies have proposed using Doppler velocity measurements in addition to reflectivity (Frisch et al. 1995; Babb et al. 1999; Kollias et al. 2001a). Others have suggested combining radar observations with other cloud remote sensing instruments such as lidars (O’Connor et al. 2005) or passive microwave radiometers in an effort to constrain the retrieval (Liao and Sassen 1994; Frisch et al. 1998; Ovtchinnikov and Kogan 2000).

_{l}Using observable Doppler parameters, such as mean Doppler velocity and the Doppler velocity spectrum width, adds another dimension of complexity, since those parameters depend both on moments of the DSD (which have the direct bearing on microphysical parameters to be retrieved) and on the turbulent air velocity. From the standpoint of *microphysical* retrievals, the turbulent contribution represents noise that must be filtered in order to extract useful microphysical information. The “noise” signals in the mean Doppler velocity and the mean Doppler spectrum width have different magnitudes and will be investigated in more detail elsewhere. Here we will concentrate on the assessment of the accuracy of retrievals based on specified sets of Doppler radar parameters. Specifically we aim to determine 1) which Doppler parameters improve the accuracy of retrievals most significantly compared to retrievals based on *Z* alone, and 2) the maximum retrieval accuracy that can be achieved based on these parameters.

Our evaluation is based on the concept of the Observing System Simulation Experiments (OSSEs) (Parsons and Dudhia 1997). Based on this concept, cloud radar parameters are obtained from data generated by the high-resolution Cooperative Institute for Mesoscale Meteorological Studies large-eddy simulation model with explicit microphysics (CIMMS LES EMP). Applying the OSSE framework for stratocumulus clouds, we quantitatively evaluate the errors of several cloud liquid water and drizzle flux retrievals. As both *V _{d}* and

*σ*are defined as intrinsic parameters of the DSD, and thus neglect the contribution from air turbulence in the sensed volume, our assessment should be considered as the lower limit on the retrieval errors.

_{d}## 2. Approach

### a. Model and data

The study is based on the CIMMS LES model, which combines 3D dynamics with an explicit (size resolving) formulation of liquid phase microphysical processes. The thermodynamic state is described in terms of virtual liquid water potential temperature and total water mixing ratio. Cloud physics processes are formulated based on prediction equations for cloud condensation nuclei and cloud/drizzle drops (19 and 25 bins, respectively). A detailed description of the model can be found in Kogan (1991), Kogan et al. (1995), and Khairoutdinov and Kogan (1999). Individual case studies and comparison of simulations with aircraft observations (Khairoutdinov and Kogan 1999; Liu et al. 2000) have demonstrated that the model can reasonably well reproduce major dynamical, radiative, and microphysical parameters. Indirect tests of a bulk drizzle parameterization derived from model DSDs (Khairoutdinov and Kogan 2000) showed good agreement with a large number of observational datasets (Wood et al. 2002; Wood 2005).

We simulated several cases of stratocumulus clouds observed during the Atlantic Stratocumulus Transition Experiment (ASTEX) field experiment in clean and polluted air masses. The simulated cloud layers represented cases with different intensities of drizzle *in the cloud* (drizzle is defined as drops in the 25–300-*μ*m radius range). Figures 1 and 2 show examples of cloud liquid water fields and cloud drop spectra superimposed with vertical wind velocities in model simulations for light and heavy drizzling conditions. Note the prevalence of bimodal cloud drop distributions in the heavy drizzle case.

From each simulation we extracted about 4000 to 6000 DSDs that were used to calculate cloud parameters, such as, for example, drop concentration, liquid water content, cloud and drizzle water content, radar reflectivity, and Doppler velocity. The set of DSDs, therefore, served as the source for deriving *Q _{l}* and

*R*retrievals using regression analysis and as a benchmark for evaluating them by comparing with the exact values of

*Q*and

_{l}*R*.

The range of cloud and drizzle parameters for all performed simulations is illustrated in Fig. 3 for datasets representing light (LD), moderate (MD), and heavy (HD) drizzle spectra. Since the cloud layer evolves significantly during the 3–6-h-long simulations, these datasets were further subdivided into subsets corresponding to a particular time of cloud evolution (e.g., LD5 refers to light drizzle case at 5 h into simulation). Table 1 shows cloud parameters and fraction of drizzle in cloud liquid water and reflectivity for subsets of data selected for analysis (LD ≡ LD1, MD ≡ MD1, and HD ≡ HD2).

### b. Definition of Doppler radar parameters

*Z*-weighted droplet terminal fall velocity

*averaged over the radar pulse volume*. The mean Doppler velocity of the pulse volume at time

*t*near point

**with coordinates (**

*x**x, y, z*) is

*f*(

_{υ}*r*) is the fall velocity of the drop with radius

*r*,

*Z*(

**) is radar reflectivity, and**

*x**F*(

_{Z}*x*)/

*Z*(

*x*) is the

*Z*-weighted drop terminal fall velocity at point

**. The bar denotes averaging over the radar pulse volume. For the K**

*x*_{α}-band vertically pointing Millimeter Cloud Radar (MMCR), which operates at 35 GHz with an 8.77-mm wavelength, the vertical gate size is 45 m.

^{1}The effective beamwidth for this radar is 0.2°–0.3°, which gives the radar a horizontal scanning dimension of about 10–50 m, depending on range. These radar pulse dimensions are comparable with those of a sampling volume required for a statistically robust determination of the drop size distribution

*n*(

**,**

*x**r*) in the full range of drop sizes, including drizzle drops. The latter, because of their low concentration, require an especially large sampling volume. The radar pulse and the DSD sampling volume dimensions are thus comparable with the grid dimensions in the LES model simulations (25 m in the vertical and 75 m in horizontal). For nondrizzling stratocumulus, the drop size distribution can be determined over a smaller sampling volume; in this case

*n*(

**,**

*x**r*) in (2) should represent DSDs averaged over the radar pulse volume.

*f*(

_{υ}*r*) for droplets in the drizzle size range is a linear function of

*r*, the Doppler velocity is essentially proportional to the ratio of the seventh to the sixth moment of the DSD (

*M*

_{7}/

*M*

_{6}). It is also worth noting that the drizzle flux is proportional to the fourth moment of the DSD (

*M*

_{4}).

The observed spectral variance is the sum of variances representing contributions from air turbulence *σ*^{2}_{w}, intrinsic microphysical variance due to the spread of drops terminal fall velocities *σ*^{2}_{d}, and the cross correlation between fluctuations of air and drop fall velocities *σ*^{2}* _{dw}*. The latter term is difficult to evaluate without information on the subgrid fluctuations in an LES model; however, our estimate based on resolvable scale fluctuations show that this term is on the same order of magnitude as

*σ*

^{2}

_{w}. In addition to the terms shown in (4), the expression for the observed spectral variance also includes other contributions, for example, due to drop oscillation/wobbling and finite width of the radar beam.

Estimates of these contributions have been made in many studies (see Doviak and Zrnic 1993 for review). Babb et al. (1999) derived an expression for the contributions of turbulence, modal diameter, and DSD shape to Doppler spectral density. They demonstrated that characteristic turbulent intensity, represented by the half-width of the vertical velocity distribution, can be recovered via minimization of a cost function. Kollias et al. (2001b) assumed a turbulence spectrum based on homogeneous energy dissipation, which can be integrated over relevant wavenumbers to obtain an expression for variance. This technique provided the best estimate of variance in the interior regions of updrafts and downdrafts, away from cloud boundaries where shear-generated turbulence tends to dominate. O’Connor et al. (2005) extended the method of Kollias et al. in order to separate the DSD and turbulent components of the variance. Their estimate of the turbulent contribution was a wind-speed-dependent fraction of the total variance calculated over a 30-s sampling time.

*V*and the Doppler spectral width

_{d}*σ*. For the intrinsic microphysical contribution, the latter is defined as

_{d}## 3. Results

*Q*retrieval based on

_{l}*Z*and

*V*we seek

_{d}*Q*in the form

_{l}*α, β, γ*are determined from a regression analysis of the LES data. By replacing the intrinsic microphysical parameter

*V*with the observable Doppler radar velocity

_{d}*V*

_{dv}=

*V*+

_{d}*w*, we can rewrite (6) in the form

*Z*and

_{d}*V*

_{dv}as an “effective” reflectivity

*Z*

_{eff}=

*Z*−

_{d}*γ*/

*β V*

_{dv}. Such a formulation is convenient because in the case of stratocumulus topped boundary layer the vertical velocity, horizontally averaged over the large time or spatial interval, is near zero (see, e.g., Frisch et al. 1995; Kollias and Albrecht 2000), thus permitting in principle a method for estimating the horizontal mean value of ln

*Q*:

_{l}Our analysis shows (Kogan et al. 2005) that the regression formulas based on power instead of exponential function yield similar approximation errors; however, when employing Doppler parameters the use of an exponential function is obviously advantageous for the reasons mentioned above.

### a. Errors in the retrieval of cloud liquid water

*Z*in Fig. 4 demonstrates that

*Q*can be reasonably well represented as a function of

_{l}*Z*(

_{m}*Z*is reflectivity in mm

_{m}^{6}m

^{−3}, while

*Z*is in dB

_{d}*Z*). The best fit in the form

*R*

^{2}= 0.941. Less than 10% of the data have errors outside the (−10%, +20%) interval for the whole range of

*Q*(see Fig. 5). The success of the retrieval in this case is primarily due to the relatively simple unimodal shape of the rather narrow drop spectra with relative drop spectrum dispersion

_{l}*σ*of about 0.25 (see Fig. 2). The mean drop radius,

*R*, for the LD case is rather small (7.5

_{m}*μ*m) and the mean precipitation flux is 0.3 mm day

^{−1}. Note that the

*Q*fraction in the liquid water content

_{r}*Q*(FQ

_{l}_{Qr}) is less than 0.1%, and the fraction of reflectivity, which comes from the drizzle part of the spectrum (FZ

_{Qr}), is <1% (Table 1). Obviously this is the main reason for the success of one-parameter (1P) retrieval in this case.

*Q*and

_{l}*Z*weakens when DSDs contain a larger fraction of drizzle drops that contribute increasingly to reflectivity (78% for HD; see Table 1) Analysis of the MD dataset reveals a significant scatter in the

*Q*−

_{l}*Z*scattergram indicating that retrievals of

*Q*based on

_{l}*Z*alone are rather inaccurate (

*R*

^{2}= 0.756). However, the accuracy of the

*Q*retrieval can be substantially increased when information on Doppler velocity is included. The top panel in Fig. 6 shows that a relationship in the following form results in a rather small degree of scatter and a quite accurate retrieval of

_{l}*Q*(

_{l}*R*

^{2}= 0.969):

The *Q _{l}* retrieval based on

*Z*alone in the heavily drizzling case HD is very poor (

*R*

^{2}= 0.181). Including

*V*in the HD case (bottom panel in Fig. 6) results in a significantly improved retrieval (

_{d}*R*= 0.618) relative to that based on

^{2}*Z*alone. However, the scatter in the HD case is larger than in MD case and

*R*

^{2}has decreased from 0.969 to 0.618. As evident from Table 1, the more numerous and larger drizzle drops in the HD case contribute appreciably both to

*Z*and

*V*, (mean fraction of drizzle contribution to

_{d}*Z*increased from 17% to 78%); however, the mean fraction of

*Q*in

_{r}*Q*increased only from 4% to 14%. The retrieval errors are not uniformly distributed over the range of

_{l}*Q*(Fig. 7). They can be as large as 100% for small values of

_{l}*Q*near cloud base; however, for larger values of

_{l}*Q*the standard deviation of the errors in the HD case is less than 20%–30%. For the moderate drizzle case MD the standard deviation of the errors is less than 10% for

_{l}*Q*> 0.2 g m

_{l}^{−3}and less than 30% for the whole range of

*Q*. The dependence of errors on drizzle is quite evident from histograms shown in Fig. 8. For heavy drizzle case about 35% of data points have errors larger than 25%, while for the medium drizzle case only 3% have errors this large. The use of Doppler spectrum width

_{l}*σ*instead of Doppler velocity affects the accuracy of the

_{d}*Q*retrieval rather insignificantly (Fig. 9), thus demonstrating that both Doppler parameters have approximately the same informational potential for microphysical retrieval. The decision of which to use should be based on such considerations as, for example, which parameter has smaller contribution from the air turbulence component, signal-to-noise ratio, etc.

_{l}*V*and Doppler spectrum width

_{d}*σ*on retrievals of

_{d}*Q*may at first seem surprising given the fact that these parameters are defined through higher moments of the DSD and thus should primarily characterize the tail of the spectrum. However, simple analysis shows that

_{l}*V*, for example, correlates well with the lower moments of the DSD. For a drop spectrum characterized by a lognormal distribution with modal radius

_{d}*r*and drop spectrum logarithmic width

_{0}*σ*, the

*k*th moment of the DSD is given by (see, e.g., Frisch et al. 1995)

*Q*as Eq. (6) and the results in Figs. 6 –8 demonstrate. The drop spectra in the LES simulations deviate from idealized lognormal distributions, and the relation between

_{l}*V*and

_{d}*M*

_{4}/

*M*

_{1}

^{3}are not as straightforward as in Eq. (12); nevertheless, our estimates based on LES-derived DSDs demonstrate that correlation between

*V*and

_{d}*M*

_{4}/

*M*

_{1}

^{3}is substantial. In the HD case, for example,

*R*

^{2}= 0.461.

### b. Errors in the retrieval of drizzle flux

The retrieval of drizzle flux *R* using *Z* and *V _{d}*, is more robust than retrieval of

*Q*, obviously because

_{l}*R*,

*Z*, and

*V*, all represent higher moments of the DSD (

_{d}*M*

_{4},

*M*

_{6}, and the ratio

*M*

_{7}/

*M*

_{6}, respectively). Thus, strong correlations between them are expected, and this is indeed the case for MD and HD datasets. In the moderate drizzle case MD the use of a 2P retrieval based on

*Z*and

*V*yields a nearly perfect correlation (

_{d}*R*

^{2}= 0.997; Fig. 10, top panel). The errors of the two-parameter (2P) retrieval for the MD case are shown in the bottom panel of Fig. 10. In this moderate drizzle case the errors are less than 5% in the whole drizzle flux range, except for drizzle rates less than 0.2 mm day

^{−1}.

For the heavy drizzle case HD *R ^{2}* increased from 0.794 for the one-parameter (1P) retrieval based on

*Z*only to 0.962 when the 2P retrieval based on

*Z*and

*V*is used (Fig. 11). The standard deviation of errors in this case is approximately in the 20%–40% range for the 1P retrieval but decreases to about 10% for the 2P retrieval (Fig. 12). As in the case of

_{d}*Q*retrieval, the errors of 2P retrievals based on

_{l}*Z*–

*V*and

_{d}*Z*–

*σ*(not shown) fall approximately into the same range.

_{d}## 4. Assessing retrievals based on drop spectra from LES, in situ observations, and analytical function approximations

The current study evaluates retrievals based on the OSSE approach using DSDs from an explicit LES model. The retrieval techniques may be also evaluated either by using in situ observations or approximations of cloud DSD by known analytical functions. Most of observational data are difficult to utilize due to the presence of drizzle drops within thin, nearly invisible stratocumulus clouds (D. Lilly 1999, personal communication). Even though the concentration of drizzle drops is very low, they dominate the reflectivity and can cause liquid water content retrievals based on reflectivity alone to perform poorly. This may explain why most of the *Q _{l}* retrievals were developed and tested for nondrizzling clouds.

*Q*and

_{l}*Z:*

Equation (13) is very close to the one obtained using LES model data [Eq. (9)], as the comparison in Fig. 13 demonstrates. The LES data provide somewhat higher (∼20%) values of *Q _{l}*, which may not be surprising as the simulation was based on one particular ASTEX case A209, which represented a rather deep, 350–400-m-thick cloud layer (see Fig. 1, top panel). The observations collected during 11 days of flights, on the other hand, included data from penetrations of thin clouds with samples of very low liquid water content.

Other retrieval methods can be obtained by approximating cloud DSD by known analytical functions. Frisch et al. (1995) developed a retrieval technique that assumes cloud and drizzle drop spectra are lognormally distributed. Using our LES-derived DSDs, we are able to evaluate the accuracy of the lognormal assumption and thus indirectly evaluate the accuracy potential of this class of retrievals. As the lognormal distribution function is defined by three parameters, we use the zeroth, first, and second moments of the LES-derived DSDs for definition of the lognormal distribution function parameters. We then compare higher moments (liquid water content, rain rate, and radar reflectivity) given by lognormal distribution function with the corresponding moments of the explicit LES-derived DSDs.

*M*

_{LN}and

*M*

_{LES}are moments given by the lognormal distribution function and from the explicit LES-derived DSDs, respectively.

Clearly, for lightly precipitating clouds (case 1), the LWC can be rather accurately approximated using the lognormal distribution function (mean ɛ ∼ 3%). The lognormal distribution function can similarly well represent the drizzle flux, although with somewhat larger error (mean ɛ ∼ 9%). The representation of radar reflectivity is very sensitive to the amount of drizzle; the error can increase sharply (mean ɛ ∼ 43%) even with a slight increase in drizzle rate (see Table 2, case 2). In the case of drizzling clouds, the radar reflectivity parameter is predominantly underestimated, obviously due to the fact that a single-mode lognormal distribution function has little or no information on drizzle drops far on the right from the mode.

In summary, the approximation of DSDs by a lognormal distribution function is justified when DSDs contain predominantly small cloud drops (*r* < 25 *μ*m) and only a fraction of drizzle drops. The retrievals of cloud liquid water similar to that of Frisch et al. (1995) should be quite accurate in this case. In cases when cloud drops span over a wide range of sizes, the approximation by a single-mode lognormal function may be inadequate and result in large errors.

## 5. Conclusions

We performed simulations of marine stratocumulus clouds observed during the Atlantic Stratocumulus Transition Experiment using the CIMMS large-eddy simulation model with size-resolving microphysics. Drop size distributions (DSDs) from these simulations represented a wide range of drizzling conditions and were used to evaluate the errors of retrievals of cloud microphysical parameters based on radar reflectivity, Doppler velocity, and Doppler spectrum width. For stratocumulus clouds with negligible amount of drizzle, the retrieval of cloud liquid water based on radar reflectivity alone is quite accurate and the parameters of the *Q _{l}–Z* relationship are in good agreement with the retrieval obtained from ASTEX observations by Fox and Illingworth (1997). When drizzle is present,

*Q*is poorly retrieved based on

_{l}*Z*alone; however, the retrieval is substantially improved when Doppler velocity or Doppler spectrum width is included. For

*Q*values larger than 0.2 g m

_{l}^{−3}, the standard deviation of errors is less than 10% in the moderate drizzle case; in the heavy drizzle case the errors are less than 20%–30%. The use of Doppler spectrum width

*σ*instead of Doppler velocity decreases the accuracy of the

_{d}*Q*retrieval only insignificantly, demonstrating that both Doppler parameters have approximately the same potential for improving microphysical retrievals.

_{l}The retrieval of precipitation flux *R* is generally more robust than *Q _{l}*, evidently because

*R*(proportional in stratocumulus clouds to the fourth moment of the DSD) is more closely correlated with the drizzle portion of the DSD than is

*Q*. In stratocumulus with heavy drizzle (

_{l}*R*> 2 mm day

^{−1})

*Z*–

*R*relationships can also be substantially improved by using the two-parameter retrievals. Errors of the two-parameter retrieval for the moderate drizzle case are less than 5%. For the heavy drizzle case, employing the two-parameter retrieval reduces the standard deviation of errors of the 1P retrieval from the 20%–40% range to about 10%. We emphasize that our error estimates represent the theoretical lower bound on retrieval errors, because the actual errors will inevitably increase, first and foremost, from uncertainties in estimation contributions from air turbulence. If the latter can be constrained and minimized (as in Babb et al. 1999; Kollias et al. 2001b; O’Connor et al. 2005), then the informational potential of radar reflectivity and Doppler parameters may be sufficient for substantial improvement in retrievals of cloud liquid water and precipitation flux under a wide range of drizzling conditions.

## Acknowledgments

This research was supported by the Office of Science (BER), U.S. Department of Energy, Grant DE-FG02-05ER64062, and the ONR Grants N00014-03-1-0304 and N00014-05-1-0550.

## REFERENCES

Austin, P., Wang Y. , Pincus R. , and Kujula V. , 1995: Precipitation in stratocumulus clouds: Observational and modeling results.

,*J. Atmos. Sci.***52****,**2329–2352.Babb, D. M., Verlinde J. , and Albrecht B. A. , 1999: Retrieval of cloud microphysical parameters from 94-GHz radar Doppler power spectra.

,*J. Atmos. Oceanic Technol.***16****,**489–503.Charlson, R. J., Lovelock J. E. , Andreae M. O. , and Warren S. G. , 1987: Ocean phytoplankton, atmospheric sulphur, cloud albedo and climate.

,*Nature***326****,**655–661.Clothiaux, E. E., and Coauthors, 1999: The Atmospheric Radiation Measurement Program cloud radars: Operational modes.

,*J. Atmos. Oceanic Technol.***16****,**819–827.Doviak, R. J., and Zrnic D. S. , 1993:

*Doppler Radar and Weather Observations*. Academic Press, 562 pp.Feingold, G., Stevens B. , Cotton W. R. , and Frisch A. S. , 1996: The relationship between drop in-cloud residence time and drizzle production in numerically simulated stratocumulus clouds.

,*J. Atmos. Sci.***53****,**1108–1122.Fox, N. I., and Illingworth A. J. , 1997: The retrieval of stratocumulus cloud properties by ground-based cloud radar.

,*J. Appl. Meteor.***36****,**485–492.Frisch, A. S., Fairall C. W. , and Snider J. B. , 1995: Measurements of stratus cloud and drizzle parameters in ASTEX with a Ka-band Doppler radar and a microwave radiometer.

,*J. Atmos. Sci.***52****,**2788–2799.Frisch, A. S., Feingold G. , Fairall C. W. , Uttal T. , and Snider J. B. , 1998: On cloud radar and microwave radiometer measurements of stratus cloud liquid water profiles.

,*J. Geophys. Res.***103****,**23195–23197.Khairoutdinov, M. F., and Kogan Y. L. , 1999: A large eddy simulation model with explicit microphysics: Validation against aircraft observations of a stratocumulus-topped boundary layer.

,*J. Atmos. Sci.***56****,**2115–2131.Khairoutdinov, M. F., and Kogan Y. L. , 2000: A new cloud physics parameterization in a large-eddy simulation model of marine stratocumulus.

,*Mon. Wea. Rev.***128****,**229–243.Kogan, Y. L., 1991: The simulation of a convective cloud in a 3-D model with explicit microphysics. Part I: Model description and sensitivity experiments.

,*J. Atmos. Sci.***48****,**1160–1189.Kogan, Y. L., Khairoutdinov M. P. , Lilly D. K. , Kogan Z. N. , and Liu Q. , 1995: Modeling of stratocumulus cloud layers in a large eddy simulation model with explicit microphysics.

,*J. Atmos. Sci.***52****,**2923–2940.Kogan, Y. L., Kogan Z. N. , and Mechem D. B. , 2005: The enhancement of radar retrievals by the use of higher moments of drop spectrum. Preprints,

*32d Conf. on Radar Meteorology*, Albuquerque, NM, Amer. Meteor. Soc., CD-ROM, P2R.6.Kollias, P., and Albrecht B. , 2000: The turbulence structure in a continental stratocumulus cloud from millimeter-wavelength radar observations.

,*J. Atmos. Sci.***57****,**2417–2434.Kollias, P., Albrecht B. A. , and Marks F. D. Jr., 2001a: Raindrop sorting induced by vertical drafts in convective clouds.

,*Geophys. Res. Lett.***28****,**2787–2790.Kollias, P., Albrecht B. A. , Lhermitte R. , and Savtchenko A. , 2001b: Radar observations of updrafts, downdrafts, and turbulence in fair-weather cumuli.

,*J. Atmos. Sci.***58****,**1750–1766.Lau, K-M., Kim J-H. , and Sud Y. , 1996: Intercomparison of hydrologic processes in AMIP GCMs.

,*Bull. Amer. Meteor. Soc.***77****,**2209–2227.Liao, L., and Sassen K. , 1994: Investigation of relationships between Ka-band radar reflectivity and ice and liquid water contents.

,*Atmos. Res.***34****,**231–248.Liu, Q., Kogan Y. L. , Lilly D. K. , Johnson D. W. , Innis G. E. , Durkee P. A. , and Nielson K. , 2000: LES modeling of ship track formation and its sensitivity to boundary layer structure.

,*J. Atmos. Sci.***57****,**2779–2791.Moran, K. P., Martner B. E. , Post M. J. , Kropfli R. A. , Welsh D. C. , and Wildener K. B. , 1998: An unattended cloud-profiling radar for use in climate research.

,*Bull. Amer. Meteor. Soc.***79****,**443–455.O’Connor, E. J., Hogan R. J. , and Illingworth A. J. , 2005: Retrieving stratocumulus drizzle parameters using Doppler radar and lidar.

,*J. Appl. Meteor.***44****,**14–27.Ovtchinnikov, M., and Kogan Y. L. , 2000: Retrieval of cloud liquid water profile in stratiform clouds from radar reflectivity measurements: Algorithm assessment using large-eddy simulations.

,*J. Geophys. Res.***105****,**17351–17359.Paluch, I. R., and Lenschow D. H. , 1991: Stratiform cloud formation in the marine boundary layer.

,*J. Atmos. Sci.***48****,**2141–2158.Parsons, D. B., and Dudhia J. , 1997: Observing System Simulation Experiments and objective analysis tests in support of the goals of the Atmospheric Radiation Measurement Program.

,*Mon. Wea. Rev.***125****,**2353–2381.Ramanathan, V., Cess R. D. , Harrison E. F. , Minnis P. , Barkstorm B. R. , Ahmad E. , and Hartman D. , 1989: Cloud-radiative forcing and climate: Results from the earth radiation budget experiment.

,*Science***243****,**57–63.Sauvageot, H., and Omar J. , 1987: Radar reflectivity of cumulus clouds.

,*J. Atmos. Oceanic Technol.***4****,**264–272.Stevens, B., Cotton W. R. , Feingold G. , and Moeng C-H. , 1998: Large-eddy simulations of strongly precipitating, shallow stratocumulus-topped boundary layers.

,*J. Atmos. Sci.***55****,**3616–3638.vanZanten, M. C., and Stevens B. , 2005: Observations of the structure of heavily precipitating marine stratocumulus.

,*J. Atmos. Sci.***62****,**4327–4342.Wood, R., 2005: Drizzle in stratiform boundary layer clouds. Part II: Microphysical aspects.

,*J. Atmos. Sci.***62****,**3034–3050.Wood, R., Field P. R. , and Cotton W. R. , 2002: Autoconversion rate bias in stratiform boundary layer cloud parameterizations.

,*Atmos. Res.***65****,**109–128.

Cloud droplet spectra in (top) light and (bottom) heavy drizzle simulations superimposed with isolines of the vertical velocity (m s^{−1}) in a vertical cross section (*y* = 1.425 km) through the cloud. At each spatial point the square box representing the drop mass distribution function *vs* log *r* is shown. The big box in the lower-left corner, shown for illustration, has two tick marks at the bottom corresponding to 10- and 100-*μ*m drop radius.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

Cloud droplet spectra in (top) light and (bottom) heavy drizzle simulations superimposed with isolines of the vertical velocity (m s^{−1}) in a vertical cross section (*y* = 1.425 km) through the cloud. At each spatial point the square box representing the drop mass distribution function *vs* log *r* is shown. The big box in the lower-left corner, shown for illustration, has two tick marks at the bottom corresponding to 10- and 100-*μ*m drop radius.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

Cloud droplet spectra in (top) light and (bottom) heavy drizzle simulations superimposed with isolines of the vertical velocity (m s^{−1}) in a vertical cross section (*y* = 1.425 km) through the cloud. At each spatial point the square box representing the drop mass distribution function *vs* log *r* is shown. The big box in the lower-left corner, shown for illustration, has two tick marks at the bottom corresponding to 10- and 100-*μ*m drop radius.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

Range of cloud parameters in the analyzed cases of stratocumulus cloud layers. The black square represents the mean and the error bars the std dev of a parameter.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

Range of cloud parameters in the analyzed cases of stratocumulus cloud layers. The black square represents the mean and the error bars the std dev of a parameter.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

Range of cloud parameters in the analyzed cases of stratocumulus cloud layers. The black square represents the mean and the error bars the std dev of a parameter.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

Scattergram of cloud liquid water as a function of reflectivity Z for the light drizzle case LD: *Q _{l}* in g m

^{−3},

*Z*in dB

_{d}*Z*, and

*Z*in mm

_{m}^{6}m

^{−3};

*R*

^{2}is the square of correlation coefficient, often referred to as coefficient of determination (COD).

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

Scattergram of cloud liquid water as a function of reflectivity Z for the light drizzle case LD: *Q _{l}* in g m

^{−3},

*Z*in dB

_{d}*Z*, and

*Z*in mm

_{m}^{6}m

^{−3};

*R*

^{2}is the square of correlation coefficient, often referred to as coefficient of determination (COD).

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

Scattergram of cloud liquid water as a function of reflectivity Z for the light drizzle case LD: *Q _{l}* in g m

^{−3},

*Z*in dB

_{d}*Z*, and

*Z*in mm

_{m}^{6}m

^{−3};

*R*

^{2}is the square of correlation coefficient, often referred to as coefficient of determination (COD).

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The cumulative distribution of *Q _{l}* retrieval errors for the LD case.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The cumulative distribution of *Q _{l}* retrieval errors for the LD case.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The cumulative distribution of *Q _{l}* retrieval errors for the LD case.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The retrieval of cloud liquid water as a function of reflectivity and Doppler velocity (*V _{d}*). (top) The moderate drizzle case MD and (bottom) the heavy drizzle case HD.

*Q*in g m

_{l}^{−3},

*Z*in dB

_{d}*Z*,

*V*in cm s

_{d}^{−1}.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The retrieval of cloud liquid water as a function of reflectivity and Doppler velocity (*V _{d}*). (top) The moderate drizzle case MD and (bottom) the heavy drizzle case HD.

*Q*in g m

_{l}^{−3},

*Z*in dB

_{d}*Z*,

*V*in cm s

_{d}^{−1}.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The retrieval of cloud liquid water as a function of reflectivity and Doppler velocity (*V _{d}*). (top) The moderate drizzle case MD and (bottom) the heavy drizzle case HD.

*Q*in g m

_{l}^{−3},

*Z*in dB

_{d}*Z*,

*V*in cm s

_{d}^{−1}.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The errors of retrieval of cloud liquid water as a function of *Q _{l}*. The solid and dashed black lines are the MD and HD mean errors; the shading areas represent the mean ±1 std dev. Light/dark gray shading corresponds to the HD/MD case, respectively.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The errors of retrieval of cloud liquid water as a function of *Q _{l}*. The solid and dashed black lines are the MD and HD mean errors; the shading areas represent the mean ±1 std dev. Light/dark gray shading corresponds to the HD/MD case, respectively.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The errors of retrieval of cloud liquid water as a function of *Q _{l}*. The solid and dashed black lines are the MD and HD mean errors; the shading areas represent the mean ±1 std dev. Light/dark gray shading corresponds to the HD/MD case, respectively.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The comparison of *Q _{l}* retrieval errors in the MD and HD cases based on two parameters,

*Z*and

_{d}*V*.

_{d}Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The comparison of *Q _{l}* retrieval errors in the MD and HD cases based on two parameters,

*Z*and

_{d}*V*.

_{d}Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The comparison of *Q _{l}* retrieval errors in the MD and HD cases based on two parameters,

*Z*and

_{d}*V*.

_{d}Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The comparison of *Q _{l}* retrieval errors in the HD case based on two parameters:

*Z*–

_{d}*V*(dashed) and

_{d}*Z*–

_{d}*σ*(solid).

_{d}Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The comparison of *Q _{l}* retrieval errors in the HD case based on two parameters:

*Z*–

_{d}*V*(dashed) and

_{d}*Z*–

_{d}*σ*(solid).

_{d}Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The comparison of *Q _{l}* retrieval errors in the HD case based on two parameters:

*Z*–

_{d}*V*(dashed) and

_{d}*Z*–

_{d}*σ*(solid).

_{d}Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

(top) The scatterplot of the retrieved vs exact drizzle flux for the MD case. (bottom) The errors of the 2P (*Z*–*V _{d}*) retrieval. The black line shows the mean error; the shading area represents the mean ±1 std dev.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

(top) The scatterplot of the retrieved vs exact drizzle flux for the MD case. (bottom) The errors of the 2P (*Z*–*V _{d}*) retrieval. The black line shows the mean error; the shading area represents the mean ±1 std dev.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

(top) The scatterplot of the retrieved vs exact drizzle flux for the MD case. (bottom) The errors of the 2P (*Z*–*V _{d}*) retrieval. The black line shows the mean error; the shading area represents the mean ±1 std dev.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The scatterplot of the retrieved vs exact drizzle flux for the HD case based on one and two (*Z*–*V _{d}*) parameters (to reduce clatter only a fraction of data points is shown).

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The scatterplot of the retrieved vs exact drizzle flux for the HD case based on one and two (*Z*–*V _{d}*) parameters (to reduce clatter only a fraction of data points is shown).

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The scatterplot of the retrieved vs exact drizzle flux for the HD case based on one and two (*Z*–*V _{d}*) parameters (to reduce clatter only a fraction of data points is shown).

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The errors of drizzle flux retrieval for the HD case. The black and white lines are the 1P and 2P mean errors; the shading areas represent the mean ±1 std dev. Light/dark gray shading corresponds to the 1P and 2P retrievals, respectively.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The errors of drizzle flux retrieval for the HD case. The black and white lines are the 1P and 2P mean errors; the shading areas represent the mean ±1 std dev. Light/dark gray shading corresponds to the 1P and 2P retrievals, respectively.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The errors of drizzle flux retrieval for the HD case. The black and white lines are the 1P and 2P mean errors; the shading areas represent the mean ±1 std dev. Light/dark gray shading corresponds to the 1P and 2P retrievals, respectively.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The comparison of retrievals based on LES data according to expression (6) and ASTEX observations by Fox and Illingworth (1997). The black line shows the ratio of LES to observational data.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The comparison of retrievals based on LES data according to expression (6) and ASTEX observations by Fox and Illingworth (1997). The black line shows the ratio of LES to observational data.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

The comparison of retrievals based on LES data according to expression (6) and ASTEX observations by Fox and Illingworth (1997). The black line shows the ratio of LES to observational data.

Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM603.1

Mean and std dev (in brackets) of drop spectra parameters for light (LD), moderate (MD), and heavy (HD) drizzling cases; *Q _{l}* and

*Q*: liquid and drizzle water content,

_{r}*N*and

_{c}*N*: total and drizzle concentration,

_{d}*R*and

_{m}*σ*: the mean radius and relative dispersion of drop spectrum,

*R*: drizzle flux,

*V*: Doppler velocity,

_{d}*Z*: reflectivity in dB

_{d}*Z*, FQ

_{Qr}and FZ

_{Qr}: fractions of

*Q*and

_{l}*Z*from

_{m}*Q*, respectively.

_{r}Mean and std dev of errors (%) in the liquid water content, drizzle flux, and radar reflectivity due to approximation of drop spectra by a lognormal distribution function. Case 2 contains more DSDs with larger drizzle rates compared to case 1 (mean and std dev of drizzle rates *R* are indicated in the first row).

^{1}

A more detailed description of the millimeter wavelength cloud radar parameters can be found in Moran et al. (1998) and Clothiaux et al. (1999).