Using Probability Density Functions to Derive Consistent Closure Relationships among Higher-Order Moments

Vincent E. Larson Atmospheric Science Group, Department of Mathematical Sciences, University of Wisconsin—Milwaukee, Milwaukee, Wisconsin

Search for other papers by Vincent E. Larson in
Current site
Google Scholar
PubMed
Close
and
Jean-Christophe Golaz National Research Council, Naval Research Laboratory, Monterey, California

Search for other papers by Jean-Christophe Golaz in
Current site
Google Scholar
PubMed
Close
Restricted access

We are aware of a technical issue preventing figures and tables from showing in some newly published articles in the full-text HTML view.
While we are resolving the problem, please use the online PDF version of these articles to view figures and tables.

Abstract

Parameterizations of turbulence often predict several lower-order moments and make closure assumptions for higher-order moments. In principle, the low- and high-order moments share the same probability density function (PDF). One closure assumption, then, is the shape of this family of PDFs. When the higher-order moments involve both velocity and thermodynamic scalars, often the PDF shape has been assumed to be a double or triple delta function. This is equivalent to assuming a mass-flux model with no subplume variability. However, PDF families other than delta functions can be assumed. This is because the assumed PDF methodology is fairly general.

This paper proposes closures for several third- and fourth-order moments. To derive the closures, the moments are assumed to be consistent with a particular PDF family, namely, a mixture of two trivariate Gaussians. (This PDF is also called a double Gaussian or binormal PDF by some authors.) Separately from the PDF assumption, the paper also proposes a simplified relationship between scalar and velocity skewnesses. This PDF family and skewness relationship are simple enough to yield simple, analytic closure formulas relating the moments. If certain conditions hold, this set of moments is specifically realizable. By this it is meant that the set of moments corresponds to a real Gaussian-mixture PDF, one that is normalized and nonnegative everywhere.

This paper compares the new closure formulas with both large eddy simulations (LESs) and closures based on double and triple delta PDFs. This paper does not implement the closures in a single-column model and test them interactively. Rather, the comparisons are diagnostic; that is, low-order moments are extracted from the LES and treated as givens that are input into the closures. This isolates errors in the closures from errors in a single-column model. The test cases are three atmospheric boundary layers: a trade wind cumulus layer, a stratocumulus layer, and a clear convective case. The new closures have shortcomings, but nevertheless are superior to the double or triple delta closures in most of the cases tested.

Corresponding author address: Vincent E. Larson, Department of Mathematical Sciences, University of Wisconsin—Milwaukee, P.O. Box 413, Milwaukee, WI 53201. Email: vlarson@uwm.edu

Abstract

Parameterizations of turbulence often predict several lower-order moments and make closure assumptions for higher-order moments. In principle, the low- and high-order moments share the same probability density function (PDF). One closure assumption, then, is the shape of this family of PDFs. When the higher-order moments involve both velocity and thermodynamic scalars, often the PDF shape has been assumed to be a double or triple delta function. This is equivalent to assuming a mass-flux model with no subplume variability. However, PDF families other than delta functions can be assumed. This is because the assumed PDF methodology is fairly general.

This paper proposes closures for several third- and fourth-order moments. To derive the closures, the moments are assumed to be consistent with a particular PDF family, namely, a mixture of two trivariate Gaussians. (This PDF is also called a double Gaussian or binormal PDF by some authors.) Separately from the PDF assumption, the paper also proposes a simplified relationship between scalar and velocity skewnesses. This PDF family and skewness relationship are simple enough to yield simple, analytic closure formulas relating the moments. If certain conditions hold, this set of moments is specifically realizable. By this it is meant that the set of moments corresponds to a real Gaussian-mixture PDF, one that is normalized and nonnegative everywhere.

This paper compares the new closure formulas with both large eddy simulations (LESs) and closures based on double and triple delta PDFs. This paper does not implement the closures in a single-column model and test them interactively. Rather, the comparisons are diagnostic; that is, low-order moments are extracted from the LES and treated as givens that are input into the closures. This isolates errors in the closures from errors in a single-column model. The test cases are three atmospheric boundary layers: a trade wind cumulus layer, a stratocumulus layer, and a clear convective case. The new closures have shortcomings, but nevertheless are superior to the double or triple delta closures in most of the cases tested.

Corresponding author address: Vincent E. Larson, Department of Mathematical Sciences, University of Wisconsin—Milwaukee, P.O. Box 413, Milwaukee, WI 53201. Email: vlarson@uwm.edu

Save
  • Abdella, K., and N. McFarlane, 1997: A new second-order turbulence closure scheme for the planetary boundary layer. J. Atmos. Sci., 54 , 18501867.

    • Search Google Scholar
    • Export Citation
  • Abdella, K., and A. C. Petersen, 2000: Third-order moment closure through a mass-flux approach. Bound.-Layer Meteor., 95 , 303318.

  • André, J. C., G. de Moor, P. Lacarrère, and R. du Vachat, 1978: Modeling the 24-hour evolution of the mean and turbulent structures of the planetary boundary layer. J. Atmos. Sci., 35 , 18611883.

    • Search Google Scholar
    • Export Citation
  • Bony, S., and K. A. Emanuel, 2001: A parameterization of the cloudiness associated with cumulus convection; evaluation using TOGA COARE data. J. Atmos. Sci., 58 , 31583183.

    • Search Google Scholar
    • Export Citation
  • Canuto, V. M., Y. Cheng, and A. Howard, 2001: New third-order moments for the convective boundary layer. J. Atmos. Sci., 58 , 11691172.

    • Search Google Scholar
    • Export Citation
  • Chung, T. J., 2002: Computational Fluid Dynamics. Cambridge University Press, 1012 pp.

  • Clarke, R. H., A. J. Dyer, R. R. Brook, D. G. Reid, and A. J. Troup, 1971: The Wangara experiment: Boundary layer data. Tech. Paper 19, Division of Meteorological Physics, CSIRO, Australia, 341 pp.

  • Cotton, W. R., and Coauthors, 2003: RAMS 2001: Current status and future directions. Meteor. Atmos. Phys., 82 , 529.

  • de Roode, S. R., P. G. Duynkerke, and A. P. Siebesma, 2000: Analogies between mass-flux and Reynolds-averaged equations. J. Atmos. Sci., 57 , 15851598.

    • Search Google Scholar
    • Export Citation
  • Deardorff, J. W., 1974: Three-dimensional numerical study of the height and mean structure of a heated planetary boundary layer. Bound.-Layer Meteor., 7 , 81106.

    • Search Google Scholar
    • Export Citation
  • Golaz, J-C., V. E. Larson, and W. R. Cotton, 2002a: A PDF-based model for boundary layer clouds. Part I: Method and model description. J. Atmos. Sci., 59 , 35403551.

    • Search Google Scholar
    • Export Citation
  • Golaz, J-C., V. E. Larson, and W. R. Cotton, 2002b: A PDF-based model for boundary layer clouds. Part II: Model results. J. Atmos. Sci., 59 , 35523571.

    • Search Google Scholar
    • Export Citation
  • Gryanik, V. M., and J. Hartmann, 2002: A turbulence closure for the convective boundary layer based on a two-scale mass-flux approach. J. Atmos. Sci., 59 , 27292744.

    • Search Google Scholar
    • Export Citation
  • Lappen, C-L., and D. A. Randall, 2001: Toward a unified parameterization of the boundary layer and moist convection. Part I: A new type of mass-flux model. J. Atmos. Sci., 58 , 20212036.

    • Search Google Scholar
    • Export Citation
  • Larson, V. E., R. Wood, P. R. Field, J-C. Golaz, T. H. Vonder Haar, and W. R. Cotton, 2001: Small-scale and mesoscale variability of scalars in cloudy boundary layers: One-dimensional probability density functions. J. Atmos. Sci., 58 , 19781994.

    • Search Google Scholar
    • Export Citation
  • Larson, V. E., J-C. Golaz, and W. R. Cotton, 2002: Small-scale and mesoscale variability in cloudy boundary layers: Joint probability density functions. J. Atmos. Sci., 59 , 35193539.

    • Search Google Scholar
    • Export Citation
  • Lewellen, W. S., and S. Yoh, 1993: Binormal model of ensemble partial cloudiness. J. Atmos. Sci., 50 , 12281237.

  • Mellor, G. L., 1977: The Gaussian cloud model relations. J. Atmos. Sci., 34 , 356358.

  • Mironov, D. V., V. M. Gryanik, V. N. Lykossov, and S. S. Zilitinkevich, 1999: Comments on “A new second-order turbulence closure scheme for the planetary boundary layer.”. J. Atmos. Sci., 56 , 34783481.

    • Search Google Scholar
    • Export Citation
  • Moeng, C-H., and Coauthors, 1996: Simulation of a stratocumulus-topped planetary boundary layer: Intercomparison among different numerical codes. Bull. Amer. Meteor. Soc., 77 , 261278.

    • Search Google Scholar
    • Export Citation
  • Pope, S. B., 2000: Turbulent Flows. Cambridge University Press, 771 pp.

  • Randall, D. A., Q. Shao, and C-H. Moeng, 1992: A second-order bulk boundary-layer model. J. Atmos. Sci., 49 , 19031923.

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

    • Search Google Scholar
    • Export Citation
  • Smith, R. N. B., 1990: A scheme for predicting layer clouds and their water content in a general circulation model. Quart. J. Roy. Meteor. Soc., 116 , 435460.

    • Search Google Scholar
    • Export Citation
  • Sommeria, G., and J. W. Deardorff, 1977: Subgrid-scale condensation in models of nonprecipitating clouds. J. Atmos. Sci., 34 , 344355.

    • Search Google Scholar
    • Export Citation
  • Tompkins, A. M., 2002: A prognostic parameterization for the subgrid-scale variability of water vapor and clouds in large-scale models and its use to diagnose cloud cover. J. Atmos. Sci., 59 , 19171942.

    • Search Google Scholar
    • Export Citation
  • Wang, S., and B. Stevens, 2000: Top-hat representation of turbulence statistics in cloud-topped boundary layers: A large eddy simulation study. J. Atmos. Sci., 57 , 423441.

    • Search Google Scholar
    • Export Citation
  • Zilitinkevich, S. S., V. M. Gryanik, V. N. Lykossov, and D. V. Mironov, 1999: Third-order transport and nonlocal turbulence closures for convective boundary layers. J. Atmos. Sci., 56 , 34633477.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1380 850 353
PDF Downloads 473 145 17