Using high resolution simulations of eight well-documented cloud cases in different climate regimes, this study investigated the statistical distributions of dynamic and thermodynamic variables in the cloud layer and examined various assumptions used by the current statistical cloud schemes. It is found that dynamic and thermodynamic variables skew differently in the cloud layer of shallow cumulus, stratocumulus, and deep convective clouds. Vertical velocity is positively skewed, but the skewed dynamic structure cannot account for the large skewness of positively skewed total mixing ratio q(t) and negatively skewed liquid water potential temperature theta(l). It is, thus, not physically sound to assume that the sub-grid variation of different variables follows the same skewed PDF. The simulations further show that the weighted standard deviations of q(t) and theta(l) have the same order of magnitude in all types of clouds, indicating that the variations of temperature and moisture are the equally important factors for sub-grid clouds. Thus, neglecting either one of them in a statistical cloud scheme may introduce significant bias in the parameterized clouds. Citation: Zhu, P., and P. Zuidema (2009), On the use of PDF schemes to parameterize sub-grid clouds, Geophys. Res. Lett., 36, L05807, doi: 10.1029/2008GL036817.