An improved double-Gaussian closure for the subgrid vertical velocity probability distribution function

Abstract

The vertical velocity probability distribution function (PDF) is analyzed throughout the depth of the lower atmosphere, including the subcloud and cloud layers, in four large-eddy simulation (LES) cases of shallow cumulus and stratocumulus. Double-Gaussian PDF closures are examined to test their ability to represent a wide range of turbulence statistics, from stratocumulus cloud layers characterized by Gaussian turbulence to shallow cumulus cloud layers displaying strongly non-Gaussian turbulence statistics. While the majority of the model closures are found to perform well in the former case, the latter presents a considerable challenge. A new model closure is suggested that accounts for high skewness and kurtosis seen in shallow cumulus cloud layers. The well-established parabolic relationship between skewness and kurtosis is examined, with results in agreement with previous studies for the subcloud layer. In cumulus cloud layers, however, a modified relationship is necessary to improve performance. The new closure significantly improves the estimation of the vertical velocity PDF for shallow cumulus cloud layers, in addition to performing well for stratocumulus. In particular, the long updraft tail representing the bulk of cloudy points is much better represented and higherorder moments diagnosed from the PDF are also greatly improved. However, some deficiencies remain owing to fundamental limitations of representing highly non-Gaussian turbulence statistics with a double- Gaussian PDF.