Implementation and Evaluation of Physics‐Driven Dynamic Entrainment‐Mixing Parameterization in a Climate Model and Its Impact on Low‐Cloud Simulation
Explore the source record for details and available documents.
Engineering topics
Publications and source records attributed to Xu, Xiaoqi.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Entrainment and detrainment rates (ε and δ) constitute the most critical free parameters in mass flux schemes commonly employed for cumulus parameterizations. Recently, Zhu et al. (2021) introduced a new approach that utilizes aircraft observations to simultaneously estimate ε and δ for cumulus clouds, overcoming the limitation of other observation-based approaches that solely yield ε without offering insights into δ. This study aims to comprehensively evaluate the reliability of this new approach. First, evaluation using an Explicit Mixing Parcel Model demonstrates the capability of the new approach to back-calculate predetermined ε and δ based on the physical properties before and after the entrainment mixing. Second, evaluation using large-eddy simulations illustrates that the new approach yields consistent ε and δ profiles compared to the traditional approach. Sensitivity tests indicate a weak sensitivity of the estimated δ with the new approach to the entrained air source. A decrease in the proportion of cloudy air in the assumed detrained air leads to a reduction in the estimated δ, while ε remains unaffected. Finally, the most appropriate assumptions for entrained and detrained air are discussed. Estimating ε for cumulus parameterizations involves acquiring ambient air more than 500 m away from the cloud edge as entrained air. Due to implicit mean field approximations in the traditional approach, determining the optimal assumption for detrained air properties proves challenging. Finally, this study confirms the reliability of the new approach in estimating ε and δ, providing confidence in its application to extensive observational data and advancement in parameterization.
Explore the source record for details and available documents.
The entrainment–mixing process is key to the cloud’s dynamical and microphysical processes. However, to the best of our knowledge, no study has quantitatively linked entrainment rate and homogeneity of mixing. Additionally, their combined effects on cloud droplet relative dispersion and physical mechanisms are even more unclear. This study analyzed the aircraft measurements of cumulus clouds collected in the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems field campaign to identify these relationships. The relationship between the entrainment rate and homogeneous mixing degree was found to be moderately positive, and the underlying mechanisms related to mixing and evaporation time scales and environmental relative humidity were explored. The relationship between relative dispersion and entrainment rate changes from positive to negative with increasing entrainment rate, depending on the dominance of evaporation of large or small droplets. The relative dispersion first increases and then decreases with increasing homogeneous mixing degree, validating the results of previous model simulations. Further, the results shed new light on understanding the interactions between cloud microphysics and cloud dynamics through turbulent entrainment and mixing processes.
Parameterization of cloud microphysics is critical for accurate simulation of weather and climate, in which the characteristics of cloud particle spectrum inevitably further affect climate simulation by changing the cloud evolution and cloud radiation effects. The popular currently used double-moment cloud microphysics schemes in numerical models can predict the intercept (N 0 ) and slope (λ) parameters of cloud particle spectrum but cannot predict the spectral shape parameter (μ), which hinders accurate description of cloud physical processes in climate models. Therefore, in the present study, we built upon the ideas of previously developed triple-moment cloud microphysics scheme, considered radar reflectivity factor as the third predictor in addition to number and mass concentration, and introduced its related prediction equation to the Single Column Atmosphere Model Version 5.3 (SCAM5.3). Moreover, the relevant microphysical process formulas in the model were revised, and a triple-moment ice-cloud microphysics scheme was constructed to predict the μ of ice particle spectrum. Based on this model, a 29-day case of the Atmospheric Radiation Measurement Program in the summer of 1997 (ARM97) was simulated, and the differences in cloud fraction and radiation simulations between the double- and triple-moment schemes were analysed. The μ of ice particle spectrum predicted by the triple-moment scheme mainly ranged from 0 to 4; the peak value was around 2, and at least 85% of the μ values were greater than the default 0. Therefore, the μ=0 setting in the double-moment scheme is unreasonable. Moreover, the developed triple-moment ice-cloud microphysics scheme yielded a narrower ice particle spectrum, which was closer to the observation results of previous studies, than the double-moment cloud microphysics scheme (e.g., height=233 hPa). Furthermore, compared with the double-moment scheme, the triple-moment scheme achieved closer simulations of the cloud fraction observations, particularly for ice clouds in the upper levels. Considering the close association between cloud fraction and radiation, the error with the triple-moment scheme was smaller than that with the double-moment scheme regardless of the shortwave (surface downward shortwave radiation, net downward shortwave radiation at the top of the atmosphere, and net shortwave radiation at the surface) or longwave (surface downward longwave radiation and net longwave radiation at the surface) flux density. Finally, the improvement mechanism of the triple-moment scheme on cloud fraction and radiation simulations was explored. The major reason is that the triple-moment scheme weakens the autoconversion of ice crystal to form snow and enhances the growth process of ice crystal deposition, thereby increasing the mass concentration of ice crystals and ice cloud fraction. Overall, the developed triple-moment ice-cloud microphysics scheme can help improve the simulation ability of models for cloud-climate feedback and other processes.
Cloud microphysical properties are significantly affected by entrainment and mixing processes. However, it is unclear how the entrainment rate affects the relative dispersion of cloud droplet size distribution. Previously, the relationship between relative dispersion and entrainment rate was found to be positive or negative. To reconcile the contrasting relationships, the Explicit Mixing Parcel Model is used to determine the underlying mechanisms. When evaporation is dominated by small droplets, and the entrained environmental air is further saturated during mixing, the relationship is negative. However, when the evaporation of big droplets is dominant, the relationship is positive. Whether or not the cloud condensation nuclei are considered in the entrained environmental air is a key factor as condensation on the entrained condensation nuclei is the main source of small droplets. However, if cloud condensation nuclei are not entrained, the relationship is positive. If cloud condensation nuclei are entrained, the relationship is dependent on many other factors. High vertical velocity, relative humidity of environmental air, liquid water content, and low droplet number concentration are more likely to cause the negative relationship since new saturation is easier to achieve by evaporation of small droplets. Further, the signs of the relationship are not strongly affected by the turbulence dissipation rate, but the higher dissipation rate causes the positive relationship more significant for a larger entrainment rate. A conceptual model is proposed to reconcile the contrasting relationships. This work enhances the understanding of relative dispersion and lays a foundation for the quantification of entrainment-mixing mechanisms.
The parameterization of cloud entrainment rates has been problematic for many years, hindering the accurate representation of convective processes in large-scale models. Here, we extend our previous work on individual shallow convection to ensemble deep convection. Entrainment rates are estimated based on three-dimensional convective clouds from August 19-20 of 1999, during the Kwajalein Experiment, simulated using a high-resolution cloud-resolving model. They are found to be negatively correlated with both vertical velocity and buoyancy, and positively correlated with the vertical divergence of the vertical velocity and with the reciprocal of cloud radii. The physical mechanisms underlying these relationships are interpreted. It is found that the parameterizations with multiple properties perform better than those with a single property. Entrainment rates and relative humidity of entrained air are positively correlated at temperature higher than 0°C, but negatively correlated at temperature lower than 0°C. Relative humidity is also included in the parameterization of entrainment rates, which differs from our previous work on shallow cumulus clouds and other studies. Finally, two forms of parameterization for entrainment rate are recommended. The first treats the entrainment rate as a function of the vertical velocity and buoyancy for temperature higher than 0°C, but as a function of relative humidity and buoyancy for temperature lower than 0°C. The second involves an equation that relates entrainment rates to vertical velocity and buoyancy regardless of temperature.
Entrainment-mixing process significantly affects cloud micro/macrophysics and the evaluation of aerosol indirect effects. However, it remains an open question as to how to quantify entrainment-mixing mechanism for broad cloud droplet size distributions (CDSDs). Here, CDSDs with different spectral widths are used to initialize the Explicit Mixing Parcel Model to address this problem, and microphysical properties are compared for different CDSD widths. For relatively broad CDSDs, as the number concentration and liquid water content decrease, the volume-mean radius and mean radius increase, which is opposite to the scenario for relatively narrow CDSDs and the homogeneous/inhomogeneous conceptual model. For the relatively broad CDSDs, such relationships could be mistaken as inhomogeneous mixing with subsequent ascent in analysis of the conventional microphysical mixing diagram. This then causes traditional homogenous mixing degrees to be unrealistically negative and not applicable for relatively broad CDSDs. New measures are introduced to explicitly account for the effect of CDSD spectral shapes on quantifying entrainment-mixing mechanism. The new measures yield reasonable ranges of values that conform to the dynamical measures such as the Damköhler and transitional scale numbers, and their temporal variation can be explained physically. This study extends the measures of homogeneous mixing degrees from narrow to broad CDSDs, and sheds new light on the consideration of CDSD shapes in parameterizations of entrainment-mixing mechanism.
Explore the source record for details and available documents.
Entrainment-mixing mechanisms significantly affect cloud droplet number concentration, radius, and spectral shape. Quantitative examination of entrainment-mixing effects on cloud droplet spectral width is lacking. Here, we examine the effects of entrainment-mixing processes on cloud microphysics by 12,218 different setups, each simulated 10 times using the Explicit Mixing Parcel Model (EMPM) driven by the observational data from the Third Tibetan Plateau Atmospheric Scientific Experiment (TIPEX-III) campaign. Parameterizations of entrainment-mixing mechanisms are developed by relating homogeneous mixing degree to transition scale number that depends on the dissipation rate and droplet evaporation time scale. The correlation between relative dispersion of cloud droplet size distribution and homogeneous mixing degree changes from negative to positive with the decreasing homogeneous mixing degree. The different relationships are closely related to the competition between complete and partial droplet evaporation and the number concentration of small droplets, which are quantitatively described by two newly introduced dimensionless numbers. The competition is significantly affected by relative humidity and mixing fraction of entrained air as well as turbulence dissipation rate, but not much by cloud droplet number concentration. Especially, when relative humidity and dissipation rate are high, there is only a negative correlation. This study sheds new light on generalizing the homogeneous/inhomogeneous concept by considering relative dispersion, and also provides parameterizations of entrainment-mixing processes and relative dispersion for atmospheric models.
Numerical simulations often overpredict precipitation over the Tibetan Plateau (TP). To examine the factors causing precipitation overprediction, different parameterizations of liquid-phase microphysical processes (accretion, autoconversion, and entrainment-mixing) are implemented into the Morrison microphysics scheme to simulate a TP precipitation event in summer with the Weather Research and Forecasting (WRF) model. The general spatial distribution and temporal trend of precipitation are captured by all simulations, but the precipitation rate is overpredicted. The results from sensitivity experiments suggest that compared to other examined liquid-phase processes, the accretion process is more important in precipitation simulation over the TP region. Further investigation with the Heidke skill scores reveals that accretion parameterization that takes into account the raindrop size produces the most accurate results in terms of the total surface precipitation. This parameterization suppresses spurious accretion and does not produce liquid-phase precipitation until cloud droplets are big enough. Here, it is also confirmed that increasing the model resolution can reduce precipitation overprediction. Results from the case study are confirmed by the use of a one-month simulation.