1. Introduction
It is now recognized that boundary layer clouds remain the largest source of uncertainty in current climate change projections (Cess et al. 1989; Bony and Dufresne 2005). Inadequate characterization of unresolved cloudiness severely affects climate predictions because of the vastly different radiative properties of stratocumulus versus shallow cumulus clouds. It is likely that the transition from stratocumulus to shallow cumulus regimes plays a key role in cloud–climate feedbacks (e.g., Teixeira et al. 2011) and it is thus crucial to understand the factors controlling this cloud transition. The stratocumulus to shallow cumulus transition has been investigated using a variety of approaches, ranging from large-scale statistical characterizations (often based on global observations) (Klein and Hartmann 1993; Wood and Hartmann 2006; Kawai and Teixeira 2010) to specific observational studies (Bretherton and Pincus 1995; Sandu et al. 2010; Mauger and Norris 2010) to mixed layer modeling studies (Bretherton and Wyant 1997) to high-resolution simulations (Krueger et al. 1995; Wyant et al. 1997).
For a better understanding of specific natural processes, it is often useful to have a hierarchy of models that allow insight into a particular problem at different levels (e.g., Held 2005). In the case of the stratocumulus to shallow cumulus transition, there seems to be a clear gap between the simplicity (and lack of causality) of observationally based statistical relations—for example, the linear relation between lower-tropospheric stability (LTS) and cloud fraction (CF) (e.g., Klein and Hartmann 1993)—and the complexity of high-resolution large-eddy simulation (LES) models. Even mixed layer models, which sit somewhere between LTS and LES models of the transition, suffer from particular problems: they can be quite complex—in a way that does not easily lead to a complete understanding of the model’s behavior—and yet they cannot even really be claimed as appropriate models of the transition, since the key variable that characterizes the transition, the cloud fraction, is not a variable in mixed layer models.
This is an unfortunate situation, and it is fair to say that simple models of the transition that do not involve the complexities of LES or mixed layer models but do include some measure of physical reasoning beyond empirical statistical correlations are clearly lacking. This is the main reason for us to propose the following simple stochastic model of the transition, where essentially we relate the decrease in cloud fraction during the transition to changes in the sea surface temperature (SST).
Our approach is to augment observed Lagrangian transition statistics with a simple model in an attempt to tease out some basic aspects of the transition. The simple model is constructed from simple physical (and dimensional) arguments that pertain to the changes experienced by the underlying physical processes across the transition. The Lagrangian data is taken from the analysis of Sandu et al. (2010), which details the statistics of CF along with the statistics of environmental factors that are often linked with the transition such as SST and large-scale divergence. They found that the transition, characterized by a decrease in CF during the first three days, is most correlated with an increase in SST; meanwhile, the large-scale divergence remains relatively constant during the same period.
2. A simple cloud transition model
As mentioned above, the development of this model is essentially driven by the following question: What is the simplest possible model that represents the key characteristics of the transition that still contains some basic physics and causality?
3. Comparison with observations
a. Description of Lagrangian data
The present model relies on observed Lagrangian statistics of SST and CF: SST as model input and CF for comparison with model output. These are taken from the Lagrangian analysis by Sandu et al. (2010) of stratocumulus to shallow cumulus cloud transitions in four eastern subtropical ocean regions, namely the northeast and southeast Pacific and Atlantic Oceans (NEP, SEP, NEA, and SEA). In their analysis, SST is taken from the interim European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-Interim), with spatial resolution of 1.5° × 1.5° and temporal resolution of 6 h, while CF is from collection 5 of the Moderate Resolution Imaging Spectroradiometer (MODIS) level-3 products, with spatial resolution of 1° × 1°, observed twice per day at 1030 and 1330 UTC.
The data of Sandu et al. (2010), reproduced in Fig. 1, shows that the median SST rises during the first 3 days and then levels off thereafter. However, other environmental factors, such as the large-scale divergence and the temperature above the boundary layer, remain essentially constant during the same period (Sandu et al. 2010). A decrease in lower-tropospheric stability is also observed, but this is associated with the increase in SST. Figure 1 shows that the CF decreases on average but exhibits considerable spread during the first 3 days.
b. Simple model with SST fluctuations
For simplicity, we assume that SST obeys the normal distribution, N(μSST, σSST), with mean μSST and standard deviation σSST taken from the Lagrangian analysis of Sandu et al. (2010); for a normal distribution, the median is equal to the mean and the interquartile range is equal to 1.35 times the standard deviation.
Figure 2 shows that the main features of the transition are reproduced by the model. The mean CF decreases systematically and the peak in the CF spread is a robust feature of the model, although these results are somewhat sensitive to the sharpness of the transition, β2 − β1. The values of β1 and β2 are fitted to the NEP data by eye and no additional tuning of the model is performed for the other three regions. The fits (β1, β2) ≈ (1.4, 2.1), (1.5, 2.0), and (1.6, 1.9) correspond to γ ≈ 0.78, 0.56, and 0.33, respectively.
c. Simple model with SST and wind fluctuations
To investigate the sensitivity of the results to the variability in wind speed, we extend the model to include the effects of wind fluctuations, V, appearing in the bulk aerodynamic transfer formula, (3). We introduce a stochastic component to the wind speed by assuming that V is independent of SST and V ~ N(7 m s−1, σV) for all time. Figure 3 shows that this simple stochastic model of the wind speed has a significant impact on the overall results: wind fluctuations increase the peak spread of CF without significantly changing the median CF. Except for the SEP region, the model shows remarkable good agreement with observations for both the CF median and interquartile range; the agreement occurs in three regions even though the key coefficients have been somewhat (loosely) optimized for just the NEP region. While we do not fully understand the discrepancy in the SEP region, we remark that the fit can be improved if F0 is lowered by 5%.
d. CF distribution
Using the NEP data and (β1, β2) ≈ (1.4, 2.1) as an indication of transition time scale, we find that the median β roughly coincides with β1 and β2 on day 0 and day 6, respectively. Although this result appears to imply a fairly smooth transition, PDFs of β and CF over the course of 6 days in Fig. 4 paint a more complex picture of how the model behaves. At any time, there are pretransition, transitioning, and posttransition air masses, corresponding to β < β1, β1 < β < β2, and β2 < β, respectively. While the PDF of β is Gaussian, the PDF of CF has two dominant peaks at C1 ≈ 1 and C2 ≈ 0.1 at any time.
Comparing the global (for all time) model CF PDF from the NEP region with the observed CF PDF from the NEA region (Mauger and Norris 2010) (see our Fig. 5), we note that the present model captures essential features of the cloudiness distribution in transition regions: the probabilities are concentrated in two peaks with small, but nonzero, probabilities in the intermediate states. However, these are not completely equal comparisons since the Lagrangian study of Mauger and Norris (2010) tracks air columns that revert to their initial CF after 6 days, but the transitions of Sandu et al. (2010) are characterized by a monotonic decrease of CF over 6 days.
e. Two-peak model
4. Concluding remarks
A lack of simple yet physically based (even if minimally so) models that are capable of representing the key characteristics of the stratocumulus to cumulus transition is the main motivation for the development of the simple model described in this paper. The model is essentially based on an assumed linear change of cloud fraction, during the transition, as a function of the ratio between two key parameters controlling the transition: the surface latent heat flux (which controls cumulus mixing) and the stratocumulus cloud-top longwave radiative flux. The surface latent heat flux (LHF) is calculated using the bulk formula with a normal PDF of stochastic sea surface temperature (SST) based on observed mean and variance from Sandu et al. (2010). The wind speed can also be based on a normal PDF. This forms the basis of a simple stochastic model, which reproduces essential features of the transition, including the peak in CF spread and the gradual decrease in CF median. The model shows good agreement with observations provided that the variability of SST and V, which influences LHF, is taken into account. According to the model, increasing the V fluctuations increases the spread in CF uniformly but leaves the median CF virtually unchanged.
Although the median CF decreases gradually over 6 days, each realization of the transition could occur more rapidly over narrower SST ranges. Such processes generate PDFs with two peaks (pre- and posttransition states) joined by infrequent transitioning states, consistent with observations. To decouple CF statistics from the precise role of controlling factors, a simple two-peak PDF model is shown to successfully predict the unknown interquartile range given the median.
Acknowledgments
The authors acknowledge the comments from three reviewers and the support provided by the Office of Naval Research, Marine Meteorology Program under Award N0001408IP20064; the NASA MAP Program; and the NOAA MAPP/CPO Program. This research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration.
REFERENCES
Bony, S., and J.-L. Dufresne, 2005: Marine boundary layer clouds at the heart of tropical cloud feedback uncertainties in climate models. Geophys. Res. Lett., 32, L20806, doi:10.1029/2005GL023851.
Bretherton, C. S., and R. Pincus, 1995: Cloudiness and marine boundary layer dynamics in the ASTEX Lagrangian experiments. Part I: Synoptic setting and vertical structure. J. Atmos. Sci., 52, 2707–2723.
Bretherton, C. S., and M. C. Wyant, 1997: Moisture transport, lower-tropospheric stability, and decoupling of cloud-topped boundary layers. J. Atmos. Sci., 54, 148–167.
Cess, R. D., and Coauthors, 1989: Interpretation of cloud–climate feedback as produced by 14 atmospheric general circulation models. Science, 245, 513–516.
Held, I. M., 2005: The gap between simulation and understanding in climate modeling. Bull. Amer. Meteor. Soc., 86, 1609–1614.
Kawai, H., and J. Teixeira, 2010: Probability density functions of liquid water path and cloud amount of marine boundary layer clouds: Geographical and seasonal variations and controlling meteorological factors. J. Climate, 23, 2079–2092.
Klein, S. A., and D. L. Hartmann, 1993: The seasonal cycle of low stratiform clouds. J. Climate, 6, 1587–1606.
Krueger, S. K., G. T. McLean, and Q. Fu, 1995: Numerical simulation of the stratus-to-cumulus transition in the subtropical marine boundary layer. Part I: Boundary-layer structure. J. Atmos. Sci., 52, 2839–2850.
Lilly, D. K., 1968: Models of cloud-topped mixed layers under a strong inversion. Quart. J. Roy. Meteor. Soc., 94, 292–309.
Mauger, G. S., and J. R. Norris, 2010: Assessing the impact of meteorological history on subtropical cloud fraction. J. Climate, 23, 2926–2940.
Sandu, I., B. Stevens, and R. Pincus, 2010: On the transitions in marine boundary layer cloudiness. Atmos. Chem. Phys., 10, 2377–2391.
Siebesma, P., and Coauthors, 2003: A large eddy simulation intercomparison study of shallow cumulus convection. J. Atmos. Sci., 60, 1201–1219.
Teixeira, J., and Coauthors, 2011: Tropical and subtropical cloud transitions in weather and climate prediction models: The GCSS/WGNE Pacific Cross-Section Intercomparison (GPCI). J. Climate, 24, 5223–5256.
Wood, R., and D. L. Hartmann, 2006: Spatial variability of liquid water path in marine low cloud: The importance of mesoscale cellular convection. J. Climate, 19, 1748–1764.
Wyant, M. C., C. S. Bretherton, H. A. Rand, and D. E. Stevens, 1997: Numerical simulations and a conceptual model of the stratocumulus to trade cumulus transition. J. Atmos. Sci., 54, 168–192.