1. Introduction
Stratocumulus clouds are the most common cloud type, covering approximately one-fifth of Earth’s surface in the annual mean, and as such have a large impact on Earth’s radiative budget (Klein and Hartmann 1993; Wood 2012). Therefore, it is important to understand the physical process that control stratocumulus cloud properties. The mixed-layer model (MLM) first proposed by Lilly (1968) has been a popular method to examine how specific physical processes impact stratocumulus cloud properties. These models are advantageous because they are computationally inexpensive and offer a quick and intuitive way to test hypotheses.
The underlying assumption of the MLM is that the stratocumulus-topped boundary layer (STBL) is well mixed; the assumption results in a zero-dimensional model that determines mixed-layer bulk properties (i.e., the STBL depth and the conserved variables: liquid water potential temperature and total water mixing ratio). Wood (2007, hereafter W07) developed an analytical equation that relates cloud thickness to the STBL depth and conserved variables of the MLM: liquid water potential temperature and total water mixing ratio. Utilizing the cloud thickness analytical equation coupled with an MLM, W07 analyzed the validity of the aerosol second indirect effect. Similarly, van der Dussen et al. (2014, hereafter VDD14) related liquid water path (LWP) to the STBL depth and the conserved variables of the MLM to determine an equilibrium value of the inversion stability parameter, beyond which a stratocumulus cloud will thin.
Both W07 and VDD14 assumed that the STBL remains well mixed as a result of the turbulence generated by longwave cooling at the top of the cloud. The liquid water lapse rate, defined as the rate at which liquid water mixing ratio changes with altitude, can be assumed to be adiabatic; the lapse rate is additionally assumed to be constant in height for the relatively thin stratocumulus clouds. Therefore, LWP is directly proportional to the square of the cloud thickness. However, the cloud thickness and LWP tendencies derived by W07 and VDD14, respectively, are not equivalent. The discrepancy stems from an error in W07’s derivation of the cloud-base-height response to changes in heat energy in the STBL. Thus, the goal of this paper is first to provide a derivation of the cloud thickness analytical equation that corrects W07’s derivation, thereby making the W07 cloud thickness tendency consistent with the VDD14 LWP tendency [provided as Eqs. (23) and (25) later]. Second, the paper uses large-eddy simulation to determine the accuracy of the LWP and cloud thickness tendency equations for a typical STBL.
2. Formulation of the cloud thickness and liquid water path tendencies
The tendencies of cloud thickness and liquid water path are formulated in terms of the two moist conserved variables:
W07 alternatively formulated the tendencies in terms of total water mixing ratio and liquid water static energy (
VDD14 arrived at the LWP tendency equation by formulating an equation for liquid water specific humidity at the top of the STBL. Here we choose to derive the equations in terms of cloud-base and inversion heights, similar to W07, as expressing the inversion height as function of the mass balance equation and having the cloud-base height respond to changes in heat and moisture content makes it easier to understand how the stratocumulus cloud layer would respond to the different physical factors such as entrainment or radiation.
a. Cloud thickness and liquid water path
b. Inversion height from mass balance
c. Cloud-base height from energy and moisture
1) Change in STBL moisture content at constant heat content
2) Change in STBL heat content at constant moisture content
d. Reconciling cloud thickness and liquid water path tendencies
3. Validation
To test how well the analytical equation of cloud thickness tendency performs, Eq. (23) was applied to two cases:
The first research flight (RF01) of DYCOMS-II (Stevens et al. 2003b) in a nocturnal STBL. The STBL was within the buoyancy reversal regime, which made the cloud deck particularly susceptible to dissolution due to run-away entrainment, which, in turn, made the simulations challenging (Stevens et al. 2003a).
The CGILS s12 case, which consists of a typical well-mixed stratocumulus over cool sea surface temperatures off the coast of California in June (Zhang et al. 2012).
In both cases, the mass, energy, and moisture tendencies of the STBL were obtained from output of the University of California, Los Angeles (UCLA) large-eddy simulation (LES; Stevens et al. 2005). The LES tendencies were then used as input to Eqs. (23) and (24) to compute the cloud thickness and LWP tendencies, respectively.
For the DYCOMS case, the vertical grid spacing is 5 m near the surface and the inversion with grid stretching in between and above the inversion, and the horizontal grid resolution is 50 m (Stevens et al. 2005). The CGILS case has vertical grid spacing that is 10 m near the surface and refined (10% per layer) to obtain a 5-m resolution near the inversion after which the grid is stretched again and the horizontal grid spacing was set at 25 m (Blossey et al. 2013). The cloud water content is a diagnostic variable based on the supersaturation, the cloud droplet number concentration is prescribed, and the droplets evolve into raindrops under the actions of the ambient flow and microphysical processes such as accretion and sedimentation (Seifert and Beheng 2006). An interactive radiation scheme was used for the CGILS case (Pincus and Stevens 2009). For the DYCOMS case, a parameterized radiation scheme was used (Stevens et al. 2003a). Very little impact on the BL evolution is therefore expected.
Note that unlike in W07 and VDD14, the analytical equation for cloud thickness tendency was not coupled to an MLM. Instead, the mass, energy, and moisture tendencies were obtained from the LES output as the goal of this paper is to provide a correct derivation of the cloud thickness tendency and not to test the validity of an MLM in simulating the STBL.
The cloud thickness derived from the analytical equation
4. Conclusions
A reconciliation of cloud thickness and LWP tendencies derived by W07 and VDD14 has been presented. W07’s derivation of the cloud-base-height response to changes in STBL heat content used the dry adiabatic lapse rate and neglected the fact that the temperature profile also shifts because of the addition or removal of heat to the STBL. Hence, when W07’s derivation is compared with the corrected response equation [Eq. (21)], W07’s derivation was found to overestimate the cloud-base-height response to changes in STBL heat content by about 22% for a typical cloud-base-height temperature of 286 K. Validation of the derived equations against LES results of the DYCOMS and CGILS s12 cases with constant and varying solar loading showed good agreement.
Following W07 and VDD14, the derived tendency equations [Eqs. (23)–(25)] can be coupled with the MLM formulation proposed by Lilly (1968). The MLM relates the heat and moisture tendencies to the different physical process occurring in the STBL, such as precipitation, entrainment, and radiation. Thus, Eqs. (23)–(25) coupled with the MLM can be utilized to study how different physical processes affect the cloud thickness.
The derived analytical Eq. (22) provides a direct relationship between cloud-base-height tendency and the heat and moisture tendencies; therefore, Eq. (22), for example, coupled with measurements of cloud-base height (e.g., with a ceilometer) could be used to validate observations of the heat and moisture budgets in the STBL. The stratocumulus cloud lifetime over land has seen renewed interest to enable accurate forecasting of solar power generation in coastal California, and the MLMs and cloud thickness tendency equations can provide insights into the importance of different terms in the moisture and heat budgets.
Acknowledgments
The authors thank the CPUC California Solar Initiative RD&D program for funding.
APPENDIX
Liquid Water Mixing Ratio Lapse Rate
REFERENCES
Albrecht, B. A., C. W. Fairall, D. W. Thomson, and A. B. White, 1990: Surface-based remote sensing of the observed and the adiabatic liquid water content of stratocumulus clouds. Geophys. Res. Lett.,17, 89–92, doi:10.1029/GL017i001p00089.
Blossey, P. N., and Coauthors, 2013: Marine low cloud sensitivity to an idealized climate change: The CGILS LES intercomparison. J. Adv. Model. Earth Syst.,5, 234–258, doi:10.1002/jame.20025.
Klein, S. A., and D. L. Hartmann, 1993: The seasonal cycle of low stratiform clouds. J. Climate, 6, 1587–1606, doi:10.1175/1520-0442(1993)006<1587:TSCOLS>2.0.CO;2.
Lilly, D. K., 1968: Models of cloud-topped mixed layers under a strong inversion. Quart. J. Roy. Meteor. Soc.,94, 292–309, doi:10.1002/qj.49709440106.
Nicholls, S., and J. Leighton, 1986: An observational study of the structure of stratiform cloud sheets: Part I. Structure. Quart. J. Roy. Meteor. Soc., 112, 431–460, doi:10.1002/qj.49711247209.
Pincus, R., and B. Stevens, 2009: Monte Carlo spectral integration: A consistent approximation for radiative transfer in large eddy simulations. J. Adv. Model. Earth Syst., 1, 1, doi:10.3894/JAMES.2009.1.1.
Seifert, A., and K. D. Beheng, 2006: A two-moment cloud microphysics parameterization for mixed-phase clouds. Part 1: Model description. Meteor. Atmos. Phys., 92 (1–2), 45–66, doi:10.1007/s00703-005-0112-4.
Stevens, B., and Coauthors, 2003a: On entrainment rates in nocturnal marine stratocumulus. Quart. J. Roy. Meteor. Soc., 129, 3469–3493, doi:10.1256/qj.02.202.
Stevens, B., and Coauthors, 2003b: Dynamics and Chemistry of Marine Stratocumulus—DYCOMS-II. Bull. Amer. Meteor. Soc., 84, 579–593, doi:10.1175/BAMS-84-5-579.
Stevens, B., and Coauthors, 2005: Evaluation of large-eddy simulations via observations of nocturnal marine stratocumulus. Mon. Wea. Rev., 133, 1443–1462, doi:10.1175/MWR2930.1.
van der Dussen, J. J., S. R. de Roode, and A. P. Siebesma, 2014: Factors controlling rapid stratocumulus cloud thinning. J. Atmos. Sci., 71, 655–664, doi:10.1175/JAS-D-13-0114.1.
Wood, R., 2007: Cancellation of aerosol indirect effects in marine stratocumulus through cloud thinning. J. Atmos. Sci., 64, 2657–2669, doi:10.1175/JAS3942.1.
Wood, R., 2012: Stratocumulus clouds. Mon. Wea. Rev., 140, 2373–2423, doi:10.1175/MWR-D-11-00121.1.
Zhang, M., C. S. Bretherton, P. N. Blossey, S. Bony, F. Brient, and J.-C. Golaz, 2012: The CGILS experimental design to investigate low cloud feedbacks in general circulation models by using single-column and large-eddy simulation models. J. Adv. Model. Earth Syst., 4, M12001, doi:10.1029/2012MS000182.