A PDFBased Microphysics Parameterization for Shallow Cumulus Clouds [Journal of the Atmospheric Sciences]
(Journal of the Atmospheric Sciences Via Acquire Media NewsEdge) ABSTRACT
Unbiased calculations of microphysical process rates such as autoconversion and accretion in mesoscale numerical weather prediction models require that subgridscale (SGS) variability over the model grid volume be taken into account. This variability can be expressed as probability distribution functions (PDFs) of microphysical variables.Using dynamically balanced largeeddy simulation (LES)model results froma case ofmarine trade cumulus, the authors develop PDFs of the cloud water, droplet concentration, and rainwater variables (qc, Nc, andqr). Both 1Dand 2Djoint PDFs (JPDFs) are presented. The authors demonstrate that accounting for the JPDFs results inmore accurate process rates for a regionalmodel grid size. Bias in autoconversion and accretion rates are presented, assuming different formulations of the JPDFs. Approximating the 2D PDF using a product of individual 1D PDFs overestimates the autoconversion rates by an order of magnitude, whereas neglecting the SGS variability altogether results in a drastic underestimate of the gridmean autoconversion rate. PDF assumptions have a much smaller impact on accretion, largely because of the nearlinear dependence of the variables in the accretion rate formula and the relatively weak correlation between qc and qr over the small LES grid volumes. The latter is attributed to the spatial decorrelation in the vertical between the two fields. Although the full PDFs are both height and time dependent, results suggest that fixedintime and fixedinheight PDFs give an acceptable level of accuracy, especially for the crucial autoconversion calculation.
(ProQuest: ... denotes formulae omitted.)
1. Introduction
Shallow marine trade cumulus are one of the most prevalent cloud types on Earth and play an important role in establishing the thermodynamic structure of the lower atmosphere in the trade latitudes. To first order, the trade inversion structure is established through a balance between the cooling and moistening effects accompanying the evaporation of trade cumulus updrafts and warming and drying from subsidence (Riehl et al. 1951; Stevens 2007). By reorganizingmoisture through the boundary layer, shallow cumulus (both marine and continental) may play a role in preconditioning the atmosphere for episodes of deeper convection (Rickenbach and Rutledge 1998;Wu et al. 2009). Although precipitation from deep convection has received the most attention, studies have shown that shallow cumulus clouds are the most common type of precipitating clouds (by sheer number) over the western tropical Pacific (e.g., Rickenbach and Rutledge 1998; Johnson et al. 1999). Satellite observations fromthe National Aeronautics and Space Administration (NASA) Tropical Rainfall Measuring Mission (TRMM) project support this conclusion about the ubiquity of shallow marine cumulus. Short and Nakamura (2008) found modes of cloudtop height at 23 km (corresponding to trade cumulus) and near 5 km (the cumulus congestus mode), with the shallow mode nearly always present. The shallow and congestus modes together contribute about 20% of the total tropical oceanic precipitation, although the bulk of the rainfall tends to come from the deeper congestus clouds.
Because shallow cumulus cloud tops typically lie below 4 km in altitude (and below the 08C isotherm), their microphysical processes are governed by the wellestablished theory of warmrainmicrophysics, which forms the basis of cloud physics parameterizations used in numerical models. The accuracy of the representation of microphysical and dynamical processes in numerical models depends on the grid scale. Shortrange (limited area) forecast models give better results when horizontal grid spacing approaches about 1 km. For example, the Met Office model has shown that finer grid spacings lead to substantial improvements in forecasting rainfall events, especially when the synoptic forcing is strong (Lean et al. 2008). Grid spacings of ;1 km may not be sufficient to fully resolve individual convective clouds; nevertheless, higher resolution aids in representing structures such as cold pools (SavicJovcic and Stevens 2008; Wang and Feingold 2009) that accompany mesoscale organization. More realistic kinetic energy spectra are also captured with horizontal grid spacing less than 10 km (Terasaki et al. 2009).
The use of highresolution simulations may be currently practical in the research mode; however, present computer capacity precludes cloudresolving representations of moist convection in longterm, realtime forecasts over large areas. In addition, higher resolution alone does not guarantee improved forecasts, as low forecasting skill is commonly attributed to inadequacies in parameterizations of moist physical processes such as microphysical and convective parameterizations. It remains vitally important to improve the accuracy of traditional microphysical and convective parameterizations.
Previous microphysical parameterizations have been developed for application at a point in the cloud, generally taken to be represented by the model grid point. The larger grid spacings used by mesoscale models [;(220) km] may contain substantial unresolved subgridscale (SGS) variability. For nonlinear process rates like autoconversion, using gridmean values of the microphysical variables can lead to substantial bias in process rates (Pincus and Klein 2000; Larson et al. 2001a; Wood et al. 2002; Larson et al. 2012) and radiative (Cahalan et al. 1994; Jeffery and Austin 2003) quantities. This bias is a consequence of Jensen's inequality [see the comprehensive discussion in Larson et al. (2001a)], which, for either convex or concave functions, is expressed as f (x) 6¼ f (x), where the overbar represents the gridvolume average. In other words, evaluating the process rate using the gridvolume average variables is not equal to themean of the process rates calculated everywhere in the grid (assuming one knows the SGS distribution of the variables). By extension, these biases in microphysical process rates can lead to biases in cloud systemproperties like precipitation rate, cloud liquid water path, and boundary layer depth (Larson et al. 2012). These differences presumably extend to the mesoscale structure of the cloud system as well via the influence of precipitation on mesoscale variability (Comstock et al. 2005; SavicJovcic and Stevens 2008;Wang and Feingold 2009). The bias potentially includes feedbacks onto the aerosol field, with variations in coalescence processing (Mechem et al. 2006) accompanying differences in precipitation rate. As the grid dimensions decrease, subgridscale variability decreases and the gridmean process rates (e.g., autoconversion) will approach the values calculated simply from the gridpoint variables (see, e.g., Khairoutdinov and Kogan 2000).
Unbiased process rates can be evaluated by integrating over the SGS probability distribution functions (PDFs) for the parameterization variables. Previous studies have established PDFs of atmospheric variables from observations (Grossman 1984; Eck and Mahrt 1991; Price and Wood 2002; Lee et al. 2010; Larson et al. 2001a,b) and numerical models (Lewellen and Yoh 1993; Xu and Randall 1996; Wang and Stevens 2000). For process rates that are a function of multiple variables, the knowledge of the joint probability distribution function (JPDF) of the different variables becomes an essential component of the parameterization.
The magnitude of the SGS variability bias is strongly dependent upon the nonlinearity of the process rates and the character of the SGS cloud variability reflected in the shape of the PDFs. Thus, an accurate cloud parameterization designed for use in a mesoscale model requires that cloud variability be taken into account. SGS cloud variability is most commonly represented using a PDFbased approach, where SGS distributions of different quantities are represented as probability density functions (Tompkins 2002; Golaz et al. 2002; Cheng and Xu 2009). At this point in time, the most sophisticated of these PDFbased approaches is the Cloud LayersUnified by Binormals (CLUBB) parameterization (Golaz et al. 2002; Larson and Golaz 2005). CLUBB assumes a joint PDF of vertical velocity, liquid water potential temperature, and total water.
Any PDFbased approach makes assumptions about the shape of the distribution. Most PDF parameterizations have assumed a prescribed analytical form for the PDF: delta functions (Randall et al. 1992; Lappen and Randall 2001),Gaussian (Sommeria andDeardorff1977), and double Gaussian (Lewellen and Yoh 1993; Golaz et al. 2002). CLUBB employs joint doubleGaussian distributions in utqtw space (Golaz et al. 2002), which are constrained from in situ aircraftobservations and largeeddy simulation (LES) output (Larson et al. 2002). Cloud observations gathered over a wide variety of conditions show substantial variations in cloud water distributions, which may be unimodal or multimodal, and either symmetrical or skewed (Tompkins 2002). Although a wide variety of PDFs have been tested, the data sources used to derive the PDFs have significant limitations because of the difficulties in obtaining comprehensive and integrated datasets over the wide range of scales. Most previous observational studies have included either in situ measurements by aircraft(Davis et al. 1996; Wood and Field 2000; Larson et al. 2001a), or retrievals from satellite instrumentation (Wielicki and Parker 1994; Barker 1996; Barker et al. 1996). Each of these data sources has its limitations. For example, aircraftprovides a onedimensional path through a cloud field, possibly leading to problems of undersampling and unrepresentative sampling. The currently available satellite data also have limitations with respect to retrieving vertical profiles of water vapor and cloud microphysical variables. Previous work has concentrated on constraining PDFs of thermodynamic quantities, for example, utqtw space, with condensate associated with only a portion of that joint parameter space.Until recently (e.g., Larson andGriffin 2013; Lebsock et al. 2013; Boutle et al. 2014), relatively little research has been performed on characterizing distributions of the cloud microphysical variables themselves [e.g., mixing ratios of cloud and rainwater (qc, qr) and cloud and rain drop concentration (Nc, Nr) in order to constrain the joint SGS variability of the microphysical quantities].
Cloudresolving models (CRMs) can provide more complete datasets that can be used to evaluate the assumptions underlying PDF approaches. Compared to observations, 3D CRM datasets are dynamically balanced1 and respond to the specified thermodynamic conditions (Lewellen and Yoh 1993; Xu and Krueger 1991; Xu and Randall 1996). However, CRMs still resolve cloudscale motions on a rather coarse grid (;1 km), and they neglect cloud variability at scales smaller than this because of the simple ''binary'' assumption of either zero or overcast fractional cloud cover in each grid cell. Another limitation of typical CRMs is the use of generic bulk microphysics parameterizations, which may not be specifically tuned to peculiarities of particular cloud types of interest (e.g., trade cumulus).
Weavoid these drawbacks by employing an LESmodel, configured with grid spacings an order of magnitude finer than those typically used by CRMs. Furthermore, the LES employed in this study uses a bulkmicrophysical parameterization2 developed from sizeresolving (bin) microphysics simulations especially formulated for shallow cumulus clouds (Kogan 2013). Output datasets obtained from the bulk microphysics LES model are used to develop PDFs of microphysical variables and scaling (closure) relationships with gridmean variables. The closure relationship specifies how gridmean variables in the mesoscale model determine the particular PDFs used to represent SGS variability. We demonstrate that using JPDFs of cloud microphysical variables obtained from LES results in greatly more accurate process rates for typical mesoscale model grid sizes.We begin by first quantifying 1DPDFs. Although substantial previous work exists along these lines, we nevertheless first explore 1D PDFs because of our finding that some JPDFs can be well expressed by the product of two 1D PDFs.
2. Model and dataset
All simulations used to develop and evaluate the PDFs were based on a version of the Cooperative Institute for Mesoscale Meteorological Studies (CIMMS) LES (Kogan et al. 1995; Khairoutdinov and Kogan 1999) called the System for Atmospheric Modeling Bulk Microphysics (SAMBM;Kogan 2013). The dynamical core of SAMBM consists of the System for Atmospheric Modeling (SAM), developed byM.Khairoutdinov (Khairoutdinov and Randall 2003). The bulk formulation of microphysics inSAMBMuses the microphysical parameterization for shallow cumulus described in Kogan (2013). Similar to the bulk parameterization of Khairoutdinov and Kogan (2000), this new parameterization is based onmultiple nonlinear regression of output from LES run with sizeresolving (bin)microphysics. The bulk autoconversion, accretion, and selfcollection rates, along with expressions for terminal velocity, are all formulated from the bin microphysics results. The new parameterization has been tested against simulations using the CIMMS LES with explicit microphysics (SAMEX; Kogan et al. 2012) in case studies of northeast Atlantic marine stratocumulus [the Atlantic Stratocumulus Experiment (ASTEX); Albrecht et al. 1995] and marine trade cumulus based on the Rain in Cumulus over the Ocean (RICO) field campaign (vanZanten et al. 2011).
Initial profiles, largescale vertical velocity, and tendencies of temperature and moisture corresponding to average conditions over the 3week period from 16 December 2004 to 8 January 2005 of the RICO field campaign were used for initialization and forcing of the LES [described in detail in the LES model intercomparison study by vanZanten et al. (2011)]. Initial cloud condensation nuclei (CCN) concentration was assumed to be 104.4cm23, a number similar to the concentrations measured by the Passive Cavity Aerosol Spectrometer Probe (PCASP) during the RICO flight RF11 on 7 January 2005 (J. Hudson 2010, personal communication). CCN concentration is a prognostic variable in the SAMBMand thus is allowed to evolve via advection, turbulence, activation, and drop evaporation; it also evolves because of coalescence processing over the course of the simulation [see Eq. (1) in Kogan (2013)]. The model horizontal and vertical grid spacings for the simulation were 100 and 40 m, respectively, with a total of 512 3 512 3 100 grid points (51.2 3 51.2 3 4km3 domain). The dynamical time step was 2 s. The simulation was run for 24 h, and model output from the last 12 h was employed in constructing the PDFs. Over the 24h simulation, this single case produced thousands of clouds having a wide range of depth, updraftintensity, and precipitation production.
3. Onedimensional PDFs
We first obtain the 1D (single variable) PDFs. Analysis of the 1D PDFs helps in understanding the variation of PDFs in the vertical, as well as the dependence of the PDFs on cloud thickness, which is one of the major factors controlling precipitation and the distribution of cloud microphysical parameters. As Fig. 1a shows, precipitation rapidly increases in the second half (1224 h) of the simulation, a period which contributes about 86% of the total accumulated surface precipitation. Since our focus is on precipitation parameterization, the analysis dataset is selected from clouds evolving during the second half of the simulation. The analysis included a total of 610 clouds, each with a horizontal cross section3 larger than 100 grid points.
Figure 1b shows the evolution of the mean and standard deviation of cloud top heights. The bulk effect of the cumulus clouds is to deepen the boundary layer, which is reflected in the increase of cloud top with time. Trade cumulus act to warm and moisten the subcloud layer and to cool and moisten the cloud layer (Siebesma et al. 2003; Stevens 2007). As the environment destabilizes with time, clouds grow higher, produce greater amounts of liquid water content, and lead to larger precipitation rates and total accumulated precipitation.
To evaluate contribution of clouds of different thickness to precipitation, we divide the dataset into four groups, each containing approximately equal numbers of clouds (an exact division is not possible because of discretization of cloudtop heights using the 40m vertical grid).Distribution of clouds by cloudtop heights is shown in Fig. 2. The first quartile (group 1) includes clouds with tops lower than 2.1 km, the second quartile (group 2) is between 2.14 and 2.26 km, group 3 is between 2.30 and 2.42km, and group 4 has tops higher than 2.46 km (the four groups are also referred to as G1G4).
The cloudthickness discretization introduces up to 20% variation in the number of clouds in each group (Fig. 3a). Figures 3b and 3c show that the mean horizontal cross section and volume of a cloud in the G4 category are more than double that of any of the other cloud categories. The mean precipitation flux at cloud base of a cloud in G4 is about twice that of a cloud in G3, and three times that of a G2 cloud (Fig. 3d). The effects of both larger precipitation rate and cross section of clouds in the G4 category combine to result in a precipitation contribution from the G4 clouds (Figs. 3e,f) substantially larger than the contribution from the other cloud categories. The total precipitation flux produced by all clouds in G1 is only 3.2% of the total amount from G2, G3, and G4 (Fig. 3f). For this reason, the development of our PDF parameterization will be based on the analysis of clouds from G2 to G4. Even though the precipitation contribution from the smaller G1 clouds is insignificant, nevertheless, their contribution to cloud cover is substantial (about 16%; Fig. 3b), and these smaller clouds cannot be neglected in studies addressing cloud albedo. The latter topic, however, is beyond our current investigation.
For PDF calculations, cloud variables in the 4km model vertical domain were subdivided into ten 400mthick layers. Using the 400mthick layers substantially increases the statistical confidence of the calculations compared to calculating PDFs for every model vertical level. The clouds in category G2 span four of these analysis layers, while G3 and G4 cloud categories span five and six layers, respectively. Layers 1 and 2 aremostly below cloud base, layers 36 cover cloud base and the middle of the cloud layer, and layers 7 and 8 represent regions near cloud top (for G3 and G4 categories).
We calculated PDFs of three cloud variables: cloud drop number concentration Nc, cloud water mixing ratio qc, and rainwater mixing ratio qr. These quantities enter the expressions for microphysical autoconversion and accretion rates, and accounting for their variability (as described by PDFs) is an important aspect of calculating unbiased microphysical process rates over a mesoscale NWPgrid. PDFs of the different microphysical variables were found by first normalizing by the mean value over each 400m layer. The PDFs of normalized parameters were binned into 15 regularly spaced categories over the normalized range from 0 to 3.
a. PDF of clouddrop concentration
The vertical structure of the Nc distributions over the 10 vertical layers is shown by contouring values of equal probability (Fig. 4). This approach is conceptually similar to the contoured frequency by altitude diagram (CFAD) plots of Yuter and Houze (1995), though here we apply it to microphysical variables rather than reflectivity. The PDF at each level was calculated if the number of cloudy points in the dataset4 exceeded a threshold value of 100. The bottom layer and the upper two layers did not contain sufficient cloud points above this threshold and therefore are not included in the PDF calculations. Qualitatively, the PDFs for G2G4 exhibit similar structure with a pronounced maximum near cloud base and increasing PDF width toward cloud top, especially pronounced for the deepest clouds (G4). This similarity may explain why the use of a PDF derived for one particular group (G4) provides a good approximation of gridaveraged conversion rates for all clouds in G2G4 (see section 5).
More details about PDF shape and form can be deduced from the plots of individual PDFs shown in Fig. 5 for G2G4. The lowerlayer PDFs (solid black lines in Fig. 5) are slightly negatively skewed for all categories (G2G4), while the distributions in the middle and upper layers of G3 and G4 display a slight positive skewness. The distributions broaden with height, with the increased frequency of large drop concentrations being possibly due to enhanced activation of small cloud drops as the updraftintensifies with altitude. The noted skewness of the distributions appears to be inconsequential, and a Gaussian function of the following form provides an accurate approximation (R2 $ 0.94 for all layers, except for the topmost layer, where R250.81; see Table 1) of the PDFs for any given layer:
...(1)
The parameter y0 here is used to scale the function f(x) to 1; the other parameters of the distribution shown in Table 1 exhibit consistent height dependence. In particular, the width (w) of the distribution increases with height, and the mode (xc) and amplitude (A) slightly decrease with height.
b. PDF of cloud water mixing ratio The vertical structure of the PDFs for cloud water mixing ratio is shown in Fig. 6. The distributions are asymmetrical, reflecting substantial positive skewness. The vertical dependence of the PDFs is weaker than for Nc and shows little systematic variation, except perhaps a slightly narrower distribution in themiddle of the cloud. Individual PDFs are shown in Fig. 7. Weapproximate the cloud water PDFs for all layers using a lognormal function:
...(2)
The parameters here have the same definitions as for the Gaussian function [Eq. (1)]. As an example, Table 2 shows the bestfit parameters for the lognormal approximation to the cloud water distribution in layer 5 for G2G4.
c. PDF of rainwater mixing ratio
The PDFs of rainwater mixing ratio shown in Figs. 8 and 9 are similar to those for the cloud watermixing ratio. Specifically, the rainwater distributions are asymmetrical with positive skewness. However, the range of rainwater distributions is narrower than for cloud water, and the distributions display a somewhat more pronounced broadening near cloud top, especially for G3 and G4. As for cloud water, the lognormal function provides a good fit (R2$0.95) for the simulationderived rainwaterPDFs. Row 2 of Table 2 shows parameters of the lognormal distribution in layer 5 for G2G4.
4. Twodimensional PDFs
Formulas for microphysical process rates typically involve at least two different microphysical variables. For example, an expression for autoconversion rate (Khairoutdinov and Kogan 2000) has the form
...(3)
suggesting the importance of using a JPDF to capture the mutual variability of these quantities. Whereas in section 3 we considered only onedimensional PDFs, here we look at the full twodimensional PDFs for qc and Nc (input variables for autoconversion) and for qc and qr (input variables for accretion). We also assess the extent to which the full twodimensional PDF can be approximated by the product of two onedimensional PDFs associated with the individual variables.
Figure 10 shows the JPDF of cloud mixing ratio qc and cloud drop number concentration Nc, the microphysical variables employed in some autoconversion parameterizations. For comparison we also show the 2D PDF, which is a product of the two onedimensional PDFs of qc and Nc. The importance of using the fully twodimensional JPDFs is obvious. The PDFs in the right panels exhibit distinct ridges aligned parallel to the x and y axes, which are artifacts of approximating the 2D PDF as the products of the two 1D PDFs [i.e., the functional form f(x)g(y)]. The 2D JPDFs on the leftpanels, on the other hand, are fully twodimensional PDFs of the form h(x, y) and do not contain the artificial ridge structures. The shape of the PDFs shows an elongated distribution on the x axis (qc), compared to more narrow distribution of the Nc parameter on the y axis. The difference between f(x)g(y) and h(x, y) indicates that the two quantities covary and thus are not independent. For example, the modal value of qc increases with increasing Nc, whereas the modal value of qc is fixed in the products of the two 1D PDFs.We note that the 2D JPDFs are rather consistent across all four layers, suggesting the possibility that cloud variability in qc andNc,might be approximated by a single, heightindependent JPDF.
The fully twodimensional JPDFs of qc and qr (the quantities used to calculate accretion in bulk models) in Fig. 11 show a somewhat different behavior than for qc and Nc. The JPDFs exhibit greater variation with height, a result of the differences in the qr distribution at different layers. In spite of the more complicated vertical variability, representing the 2D PDF by a product of the individual 1D PDFs is a better approximation for qc and qr than for qc and Nc. (Fig. 11, right panels). This result reflects the weaker correlation between qc and qr compared to that between qc and Nc. From the physical point of view, the increase in qc due to condensation frequently involves droplet activation and a corresponding increase in Nc at the same point. On the other hand, growth of qr proceeds through the coalescence of small droplets (autoconversion), sedimentation of the embryonic precipitation drops, and collection of cloud water while falling through the cloud. As a result of these processes, the maximum in the qc profile occurs in the upper third of the cloud, whereas the qr profilemaximumoccurs lower in the cloud. The correlation between cloud and rainwatermay further weaken in subsaturated downdraftregions where cloud drops evaporate more easily, but rain drops continue to grow by coalescence (Kogan et al. 2012). This spatial decorrelation between qc and qr over small grid volumes (the current LES uses 100m horizontal and 40mvertical grid) does not contradict the fact that deeper clouds with more cloud water produce, as a whole, more rain (see Fig. 3). CloudSat data, on the other hand, exhibit a robust correlation between cloud and rainwater (Lebsock et al. 2013). This apparent dilemma as to the correlation or noncorrelation of qc and qr is ultimately an issue of scale, specifically, the sampling volume over which the correlation is evaluated. Small sampling volumes such as our LES grid volume or the sampling volumes associated with aircraftobservations (Boutle et al. 2014) tend to exhibit weak correlation between qc and qr. The larger analysis volume used by Lebsock et al. is based on CloudSat observations; the latter are characterized by a vertical sampling interval of 480m and a horizontal pixel size at the surface of 1.4 3 1.7 km2. This CloudSat sampling volumeis three orders of magnitude larger than our LES grid volume and has an order of magnitude greater grid spacing in the vertical. The finer grids resolve the vertical decorrelation between qc and qr; the coarser sampling volumes do not.
In the next section, we apply these PDFs to evaluate gridmean microphysical process rates. The qualitative difference between the two PDF formulations (the ''exact'' JPDF and the 2DPDF approximated by the product of two 1D PDFs) is especially pronounced for the variables qc and Nc, resulting in a substantial impact on the gridaveraged autoconversion rates. Compared to autoconversion, the accretion rate, which employs qc and qr, is less sensitive to details about how the PDF is formulated.
5. Testing of the PDF parameterization a. General formulation
Let us introduce dimensionless parameters f and c, defined as
...(4)
where qc, Nc, and qr are horizontal averages of qc, qr, and Nc over the cloudy areas at each 400m layer. The gridmean conversion rates averaged over the cloudy areas are then given by double integrals; for example, the autoconversion rate is expressed as
...(5)
Here V is the JPDF of variables qc and Nc, and V0 is its nondimensional analog in the u and c space defined from relationship
...(6)
Introducing the PDF enhancement factor Dau, which is a direct measure of subgrid variability,
...(7)
and we can rewrite Eq. (5) as
...(8)
The local autoconversion rate is a highly nonlinear function of its two parameters. In trade wind shallow cumulus and marine boundary layer stratocumulus, it may be expressed as (Kogan 2013)
...(9)
Here qc is given in kilograms per kilogram and Nc is given per cubic centimeters. With a ' 4.2 and b ' 23, we can expect substantial difference between the exact expression given by Eq. (8) and the approximation based on gridaveraged parameters, which neglects contribution from the PDF enhancement factor Dau:
...(10)
The expression for the accretion rate has a functional form similar to the one for the autoconversion rate but involves qc and qr (Kogan 2013):
...(11)
The derivation for accretion rate is similar to the autoconversion rate:
...(12)
where C is the JPDF of variables qc and qr .
There are two substantial differences between the gridmean accretion rate in Eq. (12) and the expression for the gridmean autoconversion inEq. (5). The first is that correlation between variables qc and qr is weaker than correlation between qc and Nc. This is evidenced in Fig. 11 by the greater similarity (especially over the lower layers) between the JPDF C(qc, qr) and the distributions expressed as products of 1D PDFs of independent variables qc and qr. As a result, the general function C(qc, qr) can be approximated as
...(13)
and the expression for the gridmean accretion rate Eq. (12) can be rewritten as
...(14)
Here the overbars represent gridmean averages obtained using 1D PDFs of corresponding variables.
The second difference is that the exponents g and d in the formulation of the accretion rate are close to 1, which makes the dependence Eq. (11) nearly linear.
In this case, approximation (14) can be further simplified:
...(15)
We evaluate the formulation (14) of using the 1D individual PDFs to approximate JPDF of qc and qr and the formulation (15) where SGS variability is completely neglected in the section below using data from our simulation output.
b. Testing of different PDF approximations
The impact of using different PDF formulations in calculations of the gridmean microphysical process rates was evaluated by first calculating the conversion rates at each grid point and then averaging them for each layer. These calculations (one for autoconversion and one for accretion) establish the benchmark process rates. Each analysis layer consists of 5123 5123 10 points; however, for the RICO cloud system, the cloud fraction is about 0.1, so the number of cloudy points in each layer is only about 2 3 105. Based on the values of qc, Nc, and qr at each grid point, we calculate the exact gridpoint autoconversion and accretion rates, which are then averaged for each 400m layer to provide exact conversion rates in each layer. Applying this procedure for each hourly simulation output resulted in about 200 data points. These exact values of gridmean autoconversion and accretion are compared with values calculated based on PDFs and layeraveraged parameters according to Eqs. (4)(15).
The parameterized conversion rates are calculated for each of the three following approaches: 1) based on the JPDF; 2) based on the PDF approximated by a product of two 1D PDFs; and 3) based on the layermean variables qc,Nc, and qr. While calculations 1 and 2 use one or another formulation of PDFs for obtaining gridaveraged conversion rates, calculation 3 completely neglects SGS variability and, therefore, does not use PDFs at all.
Figure 12 shows the comparison of autoconversion and accretion process rates calculated using these three approaches. Accounting for the full variability via using the JPDF results in the most accurate values for both autoconversion and accretion rates (au2D and ac2D). Using a 2D PDF formed as a product of two 1D PDFs (au1Dand ac1D) overestimates the autoconversion rates by an order of magnitude. Neglecting variability completely (auMean and acMean), on the contrary, underestimates the autoconversion rate (Fig. 12a).
The differences between accretion rates using the three different approaches are smaller (Fig. 12b), even in approach 3, where SGS variability is neglected. As mentioned before, this weaker sensitivity to SGS variability relative to the sensitivity found in the autoconversion calculation stems from the fact that 1) the correlation between variables qc and qr is weaker than between qc and Nc and 2) the accretion process rate in Eq. (12) is given by a nearly linear function of qc and qr.
The errors of various approximations can be more clearly assessed in a plot showing ratios of approximated averages to the exact values (Fig. 13). The full JPDF provides the most accurate conversion rates, with only 7.5% and 10.5% mean error for autoconversion and accretion rates, respectively.5 Approximating the JPDF as a product of the individual 1D PDFs produces a substantial deterioration in autoconversion rates but shows a smaller effect on accretion (Fig. 13b). The range of bias increases for larger autoconversion rates, where it may lead to substantial errors in precipitation rate. Neglecting variability altogether (Fig. 13c) results in about 60% 80% underestimation of autoconversion, while large values of accretion are underestimated by up to 40%and small values are overestimated by up to 100%.
The comparisons shown in Figs. 12 and 13 were made using 1D and 2D joint formulations of PDFs, which demonstrated the superior performance of using a JPDF. However, both the JPDF and the less accurate formulation based on 1D PDFs were, nevertheless, calculated using the exact finescale distributions of cloud parameters at each time step and at each layer. In short, the PDFs employed in the above calculations were both time and height dependent, that is, they were ideally tuned to the clouds they were attempting to represent. In largescale models, accurate, timedependent information about PDFs is not available, necessitating the use of a simplified, a priori prescribed form of a PDF. We made sensitivity tests to evaluate the usefulness and accuracy of formulating such simplified JPDFs. In one case, we assumed the JPDFs were constant in time but varied in height. To obtain these fixedintime JPDFs, we used all 13 hourly data outputs from the 1224h simulation and calculated 10 PDFs for each 400m layer. In the second case, the JPDF was constant in both time and height; specifically, JPDFs V(qc,Nc) and C(qc, qr) from a fixed cloud layer were used to approximate JPDFs at all layers. Sensitivity experiments showed that the best accuracy is achieved when the JPDFs from cloud layer 6 (corresponding to the 20002400m cloud layer) were used as a proxy for all cloudy layers. Layer 6 lies in the upper third of the cloud layer and corresponds to the region of maximum updraftand precipitation growth.
Both simplified formulations of JPDFs were used in calculations of conversion rates that were compared with calculations based on the most accurate JPDFs, which depend both on time and height (referred to as exact JPDF). The results shown in Fig. 14a demonstrate that using a JPDF that is fixed in time but variable in height results in accuracy comparable to using exact JPDFs. The errors in estimating autoconversion rates (mean of 10.1% with standard deviation of 8.5%) are about the same as using the exact time and heightdependent JPDF (7.5% 67.8%). The errors in estimation of accretion rates are also comparable (13.9% 610.7% versus 10.5% 68.2% in the exact JPDF case). The use of an even simpler JPDFfixed in both time and altitude (Fig. 14b)produces rather accurate results for autoconversion rates (mean of 10.0% with standard deviation of 7.5%). For accretion rates, the ''fixed height'' JPDF formulation results in greater error (mean of 34.8%with standard deviation of 35.8%), which is unevenly distributed. For the smallest 50% accretion rates, the errors can be as large as 100%, while for larger accretion rates, the errors are in the 20%30% range.
6. Conclusions
Unbiased calculations of microphysical process rates such as autoconversion and accretion in mesoscaleNWP models require that SGS variability over the model grid volume be taken into account. This variability can be expressed as PDFs of microphysical variables, and the process rates can be integrated over these PDFs in order to obtain the unbiased gridaveraged process rates. Using dynamically balanced LES results from a case of shallow marine trade cumulus, we develop PDFs of the individual variables qc, Nc, and qr, as well as joint PDFs (JPDFs) of pairs of variables (qc, Nc) and (qc, qr). The latter are used to calculate autoconversion and accretion, respectively, in bulk models. Both 1D PDFs and 2D JPDFs are best formulated as functions of their nondimensional parameters, that is, normalized using the layermean parameter values. Nondimensionalizing the variables reduces the height dependence and, hence, allows for more robust parameterization of PDFs.
Onedimensional PDFs are useful for microphysical parameterizations based on single moments of the drop size distribution. We demonstrate that these 1D PDFs can be approximated using analytical functions and analyze the variation of their parameters in the vertical. The 1D PDF results anchor our study to the previous observational and numerical modeling studies (discussed in section 1) characterizing PDFs. The 1D PDFs of droplet concentration Nc can be well approximated by aGaussian function, which broadens somewhat with height. Cloud water (qc) and rain (qr) water mixing ratio distributions are highly skewed and can be represented by lognormal distributions, with relatively weak height dependence.
We find that JPDFs of qc and Nc (inputs to the autoconversion) cannot be represented by a simple product of their individual 1D PDF; only the ''exact'' JPDF captures correctly the covariation between qc and Nc. The product of the individual 1D PDFs of qc and qr (inputs to the accretion rate formula), on the other hand, is a better representation of the full joint JPDF of qc and qr. This lack of strong correlation between qc and qr arises from the vertical decorrelation of the two fields, since cloudwater tends to be at a maximumover the upper part of the cloud and rainwater is greatest in the lower part of the cloud. For shallow boundary layer clouds, these details in the vertical distribution of qc and qr are not resolved by large sampling volumes such as those from CloudSat.
The process rates calculated using different formulations of the PDFs are compared to the benchmark process rate calculated using the pointbypoint variability over the LES grid. Specifically, we calculate the bias in autoconversion and accretion rate. Using the JPDFs that are both height and time dependent produces the most accurate values for both autoconversion and accretion. Approximating the 2D JPDF by using a product of individual 1D PDFs overestimates the autoconversion rates by an order of magnitude, whereas neglecting the SGS variability altogether results in a drastic underestimate of the autoconversion rate. Differences in accretion rate among the different PDF assumptions are small compared to autoconversion rate differences, a result largely arising from the weaker correlation between qc and qr, and the nearlinearity of the accretion rate formula.
Because of the impracticality of using PDFs that are evolving in time during the mesoscale model simulation, we evaluated two additional PDF approximations: 1) fixed in time but height dependent and 2) fixed in time and fixed in height. The latter option would be the easier to implement in mesoscale NWP because it would necessitate only a single JPDF form for each of the two pairs of microphysical variables. Results suggest that use of a single (fixed in time and height) JPDF form for each pair of variables gives an acceptable level of accuracy, especially for autoconversion. We thus recommend, and supply in the appendix, two 15 3 15 tables corresponding to JPDFs of the variable inputs for the autoconversion and accretion rate calculations. This recommendation is justified based on Fig. 14, which shows the errors are essentially the same as using the exact JPDFs for the autoconversion rates. Similarly, small error values are associated with large values of accretion, but small accretion rates are accompanied by substantially greater error. The influence of this accretion rate error on cloud evolution cannot be determined simply by comparing instantaneous error values; rather, the influence of the errors needs to be evaluated over the course of the full timedependent mesoscale numerical simulation.
More accurate results can be obtained by using heightdependent JPDFs. In a numerical model one would need to use six tables, one for each 400mthick layer ranging fromcloud base to the highest cloud top (layers 38). This approach would be straightforward in current NWP models. For brevity, we do not show these tables in the appendix but can supply them on request.
As a final comment, we note that the PDFs developed in this paper have a characteristic scale associated with the LES domain from which they are derived. In this case, the domain size corresponds to a horizontal domain size of 51.2 3 51.2km2. Our results (not shown here) demonstrated that major cloud system properties (cloud cover, time series of surface accumulated precipitation, and mean turbulent kinetic energy) are very similar in simulations with domains of 51.2 3 51.2km2 and 25.6325.6km2 (256 points versus 512 points in each horizontal direction). Future efforts will characterize the scaledependent variability of microphysical variables in more detail using several representative mesoscale model grid sizes.
Acknowledgments. The authors are grateful to three anonymous reviewers for many constructive comments.
This investigation was supported by NOAA/Office of Oceanic and Atmospheric Research under NOAA University of Oklahoma Cooperative Agreements NA17RJ1227 and NA08OAR4320904, the U.S. Department of Commerce, ONR Grants N000141110439 and N000141110518, the National Oceanic and AtmosphericAdministration (NOAA) Climate ProgramOffice (CPO) Climate Prediction Program for the Americas/ Earth System Science Program (CPPA/ESS) Grant NA10OAR4310160, and the Department of Energy Office of Science Grant DESC0006736. The computing for this project was performed at the OU Supercomputing Center for Education and Research (OSCER) at the University of Oklahoma.
1 By ''dynamically balanced'' we mean a microphysical dataset that is produced across a realistic distribution of thermodynamical parameters: the latter can be generated in a 3D CRM or LES simulation.
2 The ''bulk'' microphysics parameterization is based on the prediction of moments of the drop size distribution function, which represent bulk cloud characteristics such as cloud/rain liquid water content and cloud/rain drop concentration.
3 Horizontal cross section is defined from liquid water path field using a threshold of 40 gm22
4 Note that the dataset included output from 13 time levels from the 1224h simulation.
5 These remaining errors are due to the use of rather coarse resolution (only 15 3 15 bins) and limiting the PDFnormalized variables' range to 3 times the mean value. Discretizing the PDF in greater detail (higher resolution and a wider range) would reduce the errors, but at the expense of a greater number of PDF bins to be integrated over. The selected number of bins (15) is a good compromise between accuracy and complexity of the discretized PDF.
REFERENCES
Albrecht, B. A., C. S. Bretherton, D. W. Johnson,W. H. Schubert, and A. S. Frisch, 1995: The Atlantic stratocumulus transition experimentASTEX. Bull. Amer. Meteor. Soc., 76, 889904.
Barker, H. W., 1996: A parameterization for computing gridaveraged solar fluxes for inhomogeneous marine boundary layer clouds. Part I: Methodology and homogeneous biases. J. Atmos. Sci., 53, 22892303.
_____, B. A. Wielicki, and L. Parker, 1996: A parameterization for computing gridaveraged solar fluxes for inhomogeneous marine boundary layer clouds. Part II: Validation using satellite data. J. Atmos. Sci., 53, 23042316.
Boutle, I. A., S. J. Abel, P. G. Hill, and C. J. Morcrette, 2014: Spatial variability of liquid cloud and rain: Observations and microphysical effects. Quart. J.Roy. Meteor. Soc., doi:10.1002/qj.2140, in press.
Cahalan, R. F., W. Ridgway, W. J. Wiscombe, and T. L. Bell, 1994: The albedo of fractal stratocumulus clouds. J. Atmos. Sci., 51, 24342455.
Cheng, A., and K.M. Xu, 2009: A PDFbased microphysics parameterization for simulation of drizzling boundary layer clouds. J. Atmos. Sci., 66, 23172334.
Comstock, K. M., C. S. Bretherton, and S. E. Yuter, 2005: Mesoscale variability and drizzle in southeast Pacific stratocumulus. J. Atmos. Sci., 62, 37923807.
Davis, A., A. Marshak, W. Wiscombe, and R. Cahalan, 1996: Scale invariance of liquid water distributions in marine stratocumulus. Part I: Spectral properties and stationarity issues. J. Atmos. Sci., 53, 15381558.
Eck, M., and L. Mahrt, 1991: A formulation for boundarylayer cloud cover. Ann. Geophys., 9, 716724.
Golaz, J.C., V. E. Larson, and W. R. Cotton, 2002: A PDFbased model for boundary layer clouds. Part I: Method and model description. J. Atmos. Sci., 59, 35403551.
Grossman, R. L., 1984: Bivariate conditional sampling of moisture flux over a tropical ocean. J. Atmos. Sci., 41, 32383253.
Jeffery, C. A., and P. H. Austin, 2003: Unified treatment of thermodynamic and optical variability in a simple model of unresolved low clouds. J. Atmos. Sci., 60, 16211631.
Johnson, R. H., T. M. Rickenbach, S. A. Rutledge, P. E. Ciesielski, and W. H. Schubert, 1999: Trimodal characteristics of tropical convection. J. Climate, 12, 23972418.
Khairoutdinov, M. F., and Y. L. Kogan, 1999: A large eddy simulation model with explicit microphysics: Validation against aircraftobservations of a stratocumulustopped boundary layer. J. Atmos. Sci., 56, 21152131.
_____, and _____, 2000: A new cloud physics parameterization for largeeddy simulation models of marine stratocumulus. Mon. Wea. Rev., 128, 229243.
_____, and D. A. Randall, 2003: Cloud resolving modeling of the ARM summer 1997 IOP: Model formulation, results, uncertainties, and sensitivities. J. Atmos. Sci., 60, 607625.
Kogan, Y. L., 2013: A cumulus cloud microphysics parameterization for cloudresolving models. J. Atmos. Sci., 70, 1423 1436.
_____, M. P. Khairoutdinov, D. K. Lilly, Z. N. Kogan, and Q. Liu, 1995: Modeling of stratocumulus cloud layers in a large eddy simulation model with explicit microphysics. J. Atmos. Sci., 52, 29232940.
_____,D. B. Mechem, and K. Choi, 2012: Effects of seasalt aerosols on precipitation in simulations of shallow cumulus. J. Atmos. Sci., 69, 463483.
Lappen, C.L., and D. A. Randall, 2001: Towards a unified parameterization of the boundary layer and moist convection. Part I: A new type of massflux model. J. Atmos. Sci., 58, 2021 2036.
Larson, V. E., and J.C. Golaz, 2005: Using probability density functions to derive consistent closure relationships among higherorder moments. Mon. Wea. Rev., 133, 10231042.
_____, and B. M. Griffin, 2013: Analytic upscaling of a local microphysics scheme. Part I: Derivation. Quart. J. Roy. Meteor. Soc., 139, 4657.
_____, R. Wood, P. R. Field, J.C. Golaz, T. H. Vonder Harr, and W. R. Cotton, 2001a: Systematic biases in the microphysics and thermodynamics of numerical models that ignore subgridscale variability. J. Atmos. Sci., 58, 11171128.
_____, _____, _____, _____, _____, and _____, 2001b: Smallscale and mesoscale variability of scalars in cloud boundary layers: Onedimensional probability density functions. J. Atmos. Sci., 58, 19781994.
_____, J.C. Golaz, and W. R. Cotton, 2002: Smallscale and mesoscale variability in cloudy boundary layers: Joint probability density functions. J. Atmos. Sci., 59, 35193539.
_____,D. P. Schanen,M.Wang, M. Ovchinnikov, and S. Ghan, 2012: PDF parameterization of boundary layer clouds in models with horizontal grid spacings from 2 to 16 km. Mon. Wea. Rev., 140, 285306.
Lean, H. W., P. A. Clark, M. Dixon, N. M. Roberts, A. Fitch, R. Forbes, and C. Halliwell, 2008: Characteristics of highresolution versions of the Met Office Unified Model for forecasting convection over the United Kingdom. Mon. Wea. Rev., 136, 34083424.
Lebsock,M., H. Morrison, and A. Gettelman, 2013: Microphysical implications of cloudprecipitation covariance derived from satellite remote sensing. J. Geophys. Res. Atmos., 118, 6521 6533, doi:10.1002/jgrd.50347.
Lee, S., B. H. Kahn, and J. Teixeira, 2010: Characterization of cloud liquid water content distributions from CloudSat. J. Geophys. Res., 115, D20203, doi:10.1029/2009JD013272.
Lewellen, W. S., and S. Yoh, 1993: Binormal model of ensemble partial cloudiness. J. Atmos. Sci., 50, 12281237.
Mechem, D. B., P. C. Robinson, and Y. L. Kogan, 2006: Processing of cloud condensation nuclei by collisioncoalescence in a mesoscale model. J. Geophys. Res., 111, D18204, doi:10.1029/ 2006JD007183.
Pincus, R., and S. A. Klein, 2000: Unresolved spatial variability and microphysical process rates in largescale models. J. Geophys. Res., 105, 27 05927 065.
Price, J. D., and R. Wood, 2002: Comparison of probability density functions for total specific humidity and saturation deficit humidity, and consequences for cloud parameterization. Quart. J. Roy. Meteor. Soc., 128, 20592072.
Randall, D. A., Q. Shao, and C.H. Moeng, 1992: A secondorder bulk boundarylayer model. J. Atmos. Sci., 49, 19031923.
Rickenbach, T. M., and S. A. Rutledge, 1998: Convection inTOGA COARE: Horizontal scale, morphology, and rainfall production. J. Atmos. Sci., 55, 27152729.
Riehl, H., T. C. Yeh, J. S. Malkus, and N. E. la Seur, 1951: The northeast trade of the Pacific Ocean. Quart. J. Roy. Meteor. Soc., 77, 598626.
SavicJovcic, V., and B. Stevens, 2008: The structure and mesoscale organization of precipitating stratocumulus. J. Atmos. Sci., 65, 15871605.
Short, D. A., and K. Nakamura, 2008: TRMM radar observations of shallow precipitation over the tropical oceans. J. Climate, 13, 41074124.
Siebesma, A. P., and Coauthors, 2003: A large eddy simulation intercomparison study of shallow cumulus convection. J. Atmos. Sci., 60, 12011219.
Sommeria, G., and J. W. Deardorff, 1977: Subgridscale condensation in models of nonprecipitating clouds. J. Atmos. Sci., 34, 344355. Stevens, B., 2007: On the growth of layers of nonprecipitating cumulus convection. J. Atmos. Sci., 64, 29162931.
Terasaki, K., H. L. Tanaka, and M. Satoh, 2009: Characteristics of the kinetic energy spectrum of NICAM model atmosphere. SOLA, 5, 180183.
Tompkins, A. M., 2002: A prognostic parameterization for the subgridscale variability of water vapor and clouds in largescale models and its use to diagnose cloud cover. J. Atmos. Sci., 59, 19171942.
vanZanten, M. C., and Coauthors, 2011: Controls on precipitation and cloudiness in simulations of tradewind cumulus as observed during RICO. J. Adv. Model. Earth Syst., 3, M06001, doi:10.1029/2011MS000056.
Wang, H., and G. Feingold, 2009: Modeling mesoscale cellular structures and drizzle in marine stratocumulus. Part I: Impact of drizzle on the formation and evolution of open cells. J. Atmos. Sci., 66, 32373256.
Wang, S., and B. Stevens, 2000: Tophat representation of turbulence statistics in cloudtopped boundary layers: A large eddy simulation study. J. Atmos. Sci., 57, 423441.
Wielicki, B. A., and L. Parker, 1994: Frequency distributions of cloud liquid water path in oceanic boundary layer cloud as a function of regional cloud fraction. Preprints, Eighth Conf. on Atmospheric Radiation, Nashville, TN, Amer.Meteor. Soc., 415417.
Wood, R., and P. R. Field, 2000: Relationships between total water, condensed water and cloud fraction examined using aircraftdata. J. Atmos. Sci., 57, 18881905.
_____, _____, and W. R. Cotton, 2002: Autoconversion rate bias in stratiform boundary layer cloud parameterizations. Atmos. Res., 65, 109128.
Wu, C.M., B. Stevens, and A. Arakawa, 2009: What controls the transition from shallow to deep convection? J. Atmos. Sci., 66, 17931806.
Xu, K.M., and S. K. Krueger, 1991: Evaluation of cloudiness parameterizations using a cumulus ensemble model. Mon. Wea. Rev., 119, 342367.
_____, and D. A. Randall, 1996: Evaluation of statistically based cloudiness parameterizations used in climate models. J. Atmos. Sci., 53, 31033119.
Yuter, S. E., and R. A. Houze Jr., 1995: Threedimensional kinematic and microphysical evolution of Florida cumulonimbus. Part II: Frequency distribution of vertical velocity, reflectivity, and differential reflectivity. Mon.Wea. Rev., 123, 19411963.
YEFIM L. KOGAN
Cooperative Institute for Mesoscale Meteorological Studies, University of Oklahoma, Norman, Oklahoma
DAVID B. MECHEM
Atmospheric Science Program, Department of Geography, University of Kansas, Lawrence, Kansas
(Manuscript received 1 July 2013, in final form 1 November 2013)
Corresponding author address: Yefim Kogan, Cooperative Institute for Mesoscale Meteorological Studies, University of Oklahoma, 120 David L. Boren Blvd., Suite 2100, Norman, OK 730727304.
Email: ykogan@ou.edu
DOI: 10.1175/JASD130193.1
APPENDIX
Tabulated Values of the Joint Probability Distribution Functions
As results shown in Fig. 14b demonstrated, using JPDFs that are fixed in time and height provides accuracy comparable to using ''exact'' JPDFs for autoconversion rates (errors about 10%) and about 30% errors for the whole range of accretion rates. Over the most important range of middle and large values of accretion rates, the accuracy is higher and comparable to that of the exact JPDF formulation. The most simple and practical way to implement PDF approach in mesoscale models, therefore, is to use values of JPDFs V(qc,Nc) and C(qc, qr) provided in Tables A1 and A2. These JPDFs were calculated for layer 6 (corresponding to the 2000 2400m cloud layer) using all 13 hourly data outputs from the 1224 h simulation. An even more accurate but still practical approach is to use six tables corresponding to layers 38 for each JPDF and thus capture the PDFs' dependence on height (see Fig. 14a for results of the errors calculations). These tables are not provided here because of space limitations, but they are available from the corresponding author by request.
(c) 2014 American Meteorological Society
[ Back To TMCnet.com's Homepage ]
