A PDF-Based Microphysics Parameterization for Shallow Cumulus Clouds

This is the published version. © Copyright 2014 American Meteorological Society (AMS). Permission to use figures, tables, and brief excerpts from this work in scientific and educational works is hereby granted provided that the source is acknowledged. Any use of material in this work that is determined to be “fair use” under Section 107 of the U.S. Copyright Act September 2010 Page 2 or that satisfies the conditions specified in Section 108 of the U.S. Copyright Act (17 USC §108, as revised by P.L. 94-553) does not require the AMS’s permission. Republication, systematic reproduction, posting in electronic form, such as on a web site or in a searchable database, or other uses of this material, except as exempted by the above statement, requires written permission or a license from the AMS. Additional details are provided in the AMS Copyright Policy, available on the AMS Web site located at (http://www.ametsoc.org/) or from the AMS at 617-227-2425 or copyrights@ametsoc.org.Unbiased calculations of microphysical process rates such as autoconversion and accretion in mesoscale numerical weather prediction models require that subgrid-scale (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 large-eddy 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 regional-model 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 grid-mean autoconversion rate. PDF assumptions have a much smaller impact on accretion, largely because of the near-linear 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 fixed-in-time and fixed-in-height PDFs give an acceptable level of accuracy, especially for the crucial autoconversion calculation