PDF Parameterization of Boundary Layer Clouds in Models with Horizontal Grid Spacings from 2 to 16 km

Many present-day numerical weather prediction (NWP) models are run at resolutions that permit deep convection. In these models, however, the boundary layer turbulence and boundary layer cloud features are still grossly underresolved. Underresolution is also present in climate models that use a multiscale modeling framework (MMF), in which a convection-permitting model is run in each grid column of a global general circulation model. To better represent boundary layer clouds and turbulence in convection-permitting models, a parameterization was developed that models the joint probability density function (PDF) of vertical velocity, heat, and moisture. Although PDF-based parameterizations are more complex and computationally expensive than many other parameterizations, in principle PDF parameterizations have several advantages. For instance, they ensure consistency of liquid (cloud) water and cloud fraction; they avoid using separate parameterizations for different cloud types such as cumulus and stratocumulus; and they have an appropriate formulation in the ‘‘terra incognita’’ in which updrafts are marginally resolved. In this paper, an implementation of a PDF parameterization is tested to see whether it improves the simulations of a state-of-the-art convection-permitting model. The PDF parameterization used is the Cloud LayersUnifiedBy Binormals (CLUBB)parameterization.The host cloud-resolving model usedis theSystem forAtmosphericModeling(SAM).SAMisrunbothwithandwithoutCLUBBimplementedinit.Simulations of two shallowcumulus(Cu)cases andtwo shallowstratocumulus(Sc)casesarerun ina 3D configurationat 2-, 4-, and 16-km horizontal grid spacings. Including CLUBB in the simulations improves some of the simulated fields—such as vertical velocity variance, horizontal wind fields, cloud water content, and drizzle water content—especially in the two Cu cases. Implementing CLUBB in SAM improves the simulations slightly at 2-km horizontal grid spacing, significantlyat4-kmgridspacing,andgreatlyat 16-kmgridspacing.Furthermore,thesimulationsthatinclude CLUBB exhibit a reduced sensitivity to horizontal grid spacing.

[1]  A. Arakawa The Cumulus Parameterization Problem: Past, Present, and Future , 2004 .

[2]  P. Bougeault,et al.  Modeling the Trade-Wind Cumulus Boundary Layer. Part II: A High-Order One-Dimensional Model , 1981 .

[3]  R. Neggers A Dual Mass Flux Framework for Boundary Layer Convection. Part II: Clouds , 2009 .

[4]  Kuan Xu,et al.  A PDF-Based Microphysics Parameterization for Simulation of Drizzling Boundary Layer Clouds , 2009 .

[5]  Vincent E. Larson,et al.  Small-Scale and Mesoscale Variability of Scalars in Cloudy Boundary Layers: One-Dimensional Probability Density Functions , 2001 .

[6]  A. P. Siebesma,et al.  A Large Eddy Simulation Intercomparison Study of Shallow Cumulus Convection , 2003 .

[7]  C. Bretherton,et al.  Subtropical Low Cloud Response to a Warmer Climate in a Superparameterized Climate Model. Part II: Column Modeling with a Cloud Resolving Model , 2009 .

[8]  Jean-Christophe Golaz,et al.  Large‐eddy simulation of the diurnal cycle of shallow cumulus convection over land , 2002 .

[9]  Evgueni I. Kassianov,et al.  The multi-scale aerosol-climate model PNNL-MMF: model description and evaluation , 2010 .

[10]  J. Wyngaard Toward Numerical Modeling in the “Terra Incognita” , 2004 .

[11]  Vincent E. Larson,et al.  A dynamic probability density function treatment of cloud mass and number concentrations for low level clouds in GFDL SCM/GCM , 2010 .

[12]  D. Randall,et al.  A Multiscale Modeling System: Developments, Applications, and Critical Issues , 2009 .

[13]  Kuan Xu,et al.  Simulation of Boundary-Layer Cumulus and Stratocumulus Clouds using a Cloud-Resolving Model With Low- and Third-Order Turbulence Closures , 2008 .

[14]  David A. Randall,et al.  Evaluation of the Simulated Interannual and Subseasonal Variability in an AMIP-Style Simulation Using the CSU Multiscale Modeling Framework , 2008 .

[15]  C. Bretherton,et al.  The University of Washington Shallow Convection and Moist Turbulence Schemes and Their Impact on Climate Simulations with the Community Atmosphere Model , 2009 .

[16]  Jon M. Reisner,et al.  A Study of Cloud Mixing and Evolution Using PDF Methods. Part I: Cloud Front Propagation and Evaporation , 2006 .

[17]  Stephen B. Pope,et al.  Calculations of Swirl Combustors Using Joint Velocity-Scalar Probability Density Function Method , 1997 .

[18]  John S. Kain,et al.  The Kain–Fritsch Convective Parameterization: An Update , 2004 .

[19]  D. Lilly,et al.  On entrainment rates in nocturnal marine stratocumulus , 2003 .

[20]  S. J. Weiss,et al.  Some Practical Considerations Regarding Horizontal Resolution in the First Generation of Operational Convection-Allowing NWP , 2008 .

[21]  C. Bretherton,et al.  Large-Eddy Simulations of a Drizzling, Stratocumulus-Topped Marine Boundary Layer , 2009 .

[22]  Brian M. Griffin,et al.  Multi-variate probability density functions with dynamics for cloud droplet activation in large-scale models: single column tests , 2010 .

[23]  A. Tompkins A Prognostic Parameterization for the Subgrid-Scale Variability of Water Vapor and Clouds in Large-Scale Models and Its Use to Diagnose Cloud Cover , 2002 .

[24]  M. Khairoutdinov,et al.  A New Cloud Physics Parameterization in a Large-Eddy Simulation Model of Marine Stratocumulus , 2000 .

[25]  A. Arakawa,et al.  Toward unification of the multiscale modeling of the atmosphere , 2011 .

[26]  A. P. Siebesma,et al.  A Combined Eddy-Diffusivity Mass-Flux Approach for the Convective Boundary Layer , 2007 .

[27]  David A. Randall,et al.  Toward a Unified Parameterization of the Boundary Layer and Moist Convection. Part I: A New Type of Mass-Flux Model , 2001 .

[28]  D. Randall,et al.  Simulations of the Atmospheric General Circulation Using a Cloud-Resolving Model as a Superparameterization of Physical Processes , 2005 .

[29]  Kevin W. Manning,et al.  Experiences with 0–36-h Explicit Convective Forecasts with the WRF-ARW Model , 2008 .

[30]  D. Randall,et al.  A cloud resolving model as a cloud parameterization in the NCAR Community Climate System Model: Preliminary results , 2001 .

[31]  M. Köhler,et al.  A Dual Mass Flux Framework for Boundary Layer Convection. Part I: Transport , 2009 .

[32]  Christopher S. Bretherton,et al.  Climate sensitivity and cloud response of a GCM with a superparameterization , 2006 .

[33]  L. Gérard,et al.  An integrated package for subgrid convection, clouds and precipitation compatible with meso‐gamma scales , 2007 .

[34]  D. Randall,et al.  Cloud resolving modeling of the ARM summer 1997 IOP: Model formulation, results, uncertainties, and sensitivities , 2003 .

[35]  Vincent E. Larson,et al.  A PDF-Based Model for Boundary Layer Clouds. Part I: Method and Model Description , 2002 .

[36]  A. Sobel,et al.  The Global Circulation of the Atmosphere , 2021 .

[37]  João Paulo Teixeira,et al.  An eddy‐diffusivity/mass‐flux parametrization for dry and shallow cumulus convection , 2004 .

[38]  Sylvie Malardel,et al.  A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction , 2009 .

[39]  G. Thompson,et al.  Impact of Cloud Microphysics on the Development of Trailing Stratiform Precipitation in a Simulated Squall Line: Comparison of One- and Two-Moment Schemes , 2009 .

[40]  Todd D. Ringler,et al.  Exploring a Multiresolution Modeling Approach within the Shallow-Water Equations , 2011 .

[41]  C. Bretherton,et al.  Evaluation of Large-Eddy Simulations via Observations of Nocturnal Marine Stratocumulus , 2005 .

[42]  James O. Pinto,et al.  Mesoscale modeling of springtime Arctic mixed-phase stratiform clouds using a new two-moment bulk microphysics scheme , 2005 .

[43]  H. Jonker,et al.  Large-Eddy Simulation: How Large is Large Enough? , 2004 .

[44]  Song‐You Hong,et al.  Improvement of the K-profile Model for the Planetary Boundary Layer based on Large Eddy Simulation Data , 2003 .

[45]  M. Dubey,et al.  The potential impacts of pollution on a nondrizzling stratus deck : Does aerosol number matter more than type? , 2008 .

[46]  V. Larson,et al.  Using Probability Density Functions to Derive Consistent Closure Relationships among Higher-Order Moments , 2005 .

[47]  David A. Randall,et al.  A second-order bulk boundary-layer model , 1992 .

[48]  J. Deardorff,et al.  Subgrid-Scale Condensation in Models of Nonprecipitating Clouds , 1977 .