DECEMBER 2012
LEE ET AL.
3607
Effect of Aerosol on Cloud–Environment Interactions in Trade Cumulus SEOUNG-SOO LEE NOAA/Earth System Research Laboratory/Chemical Sciences Division, and CIRES, University of Colorado, Boulder, Colorado
GRAHAM FEINGOLD NOAA/Earth System Research Laboratory/Chemical Sciences Division, Boulder Colorado
PATRICK Y. CHUANG Department of Earth and Planetary Sciences, University of California, Santa Cruz, Santa Cruz, California (Manuscript received 25 January 2012, in final form 29 May 2012) ABSTRACT This study examines the role of aerosol in mediating interactions between a warm trade cumulus cloud system and the environment that spawns it. Numerical simulations of the observed and well-studied Rain in Cumulus over the Ocean (RICO) field experiment are performed. The results draw on simulations of 34-h duration so as to avoid conclusions based on transients. Simulations show that, on average, aerosol-perturbed clouds are initially deeper and more vigorous but that after about 14 h there is a reversal in this trend, and unperturbed clouds deepen relative to the perturbed clouds. Differences in cloud depth are about 100 m, and differences in vertical velocity variance are about 30%. After about 20 h, most cloud fields are statistically similar with the exception of rain rate and optical depth, which are lower and higher, respectively, in the highaerosol conditions. By sampling the model output at various points in the cloud system evolution, the mechanisms responsible for the initial differences and then convergence of most of the cloud field properties are addressed. Sensitivity tests indicate that responses are driven primarily by temperature profiles, rather than by humidity profiles, and that the general trend to homogenization of the bulk cloud field properties is robust for different forcings. Finally, the paper shows that even transient aerosol perturbations may endure beyond the duration of the perturbation itself, provided they persist long enough. Short-duration aerosol perturbations are unlikely to have much influence on the system.
1. Introduction Shallow, warm cumulus clouds play an important role in the transport of heat, moisture, and pollutants to the free troposphere. They also affect the planetary albedo and thus the Earth’s radiative budget (Sengupta et al. 1990). Furthermore, the manner in which they are represented in climate models has a strong bearing on climate sensitivity (Bony and Dufresne 2005). The prevalence of these clouds is well documented (Tiedtke et al. 1988; Slingo et al. 1994). In the trade wind regime, shallow cumulus clouds both influenced by and, in turn, affect environmental conditions,
Corresponding author address: Seoung-Soo Lee, NOAA/ Chemical Sciences Division, 325 Broadway, Boulder, CO 80305. E-mail:
[email protected] DOI: 10.1175/JAS-D-12-026.1 Ó 2012 American Meteorological Society
such as temperature, humidity, and wind shear (Malkus 1954; Squires 1958; Grabowski and Clark 1993; Zhao and Austin 2005). A distribution of buoyant plumes of varying magnitudes generates cumulus clouds that transport heat and moisture into the trade inversion to different extents. In this weakly stable transition region, characterized by strong mixing and downward transport of inversion air, evaporating cloud turrets humidify the warm and descending dry air, preconditioning it for subsequent convection. The upward flux of water associated with trade cumulus enhances surface evaporation, which supplies moisture to the system. Moisture is removed by precipitation, which is a strong function of cloud depth, but also controlled by cloud microphysical processes, themselves a function of the aerosol. This introduces the potential for the aerosol to influence the thermodynamic environment that supports the trade cumulus cloud system.
3608
JOURNAL OF THE ATMOSPHERIC SCIENCES
The role of the atmospheric aerosol in influencing cloud microphysical processes and perhaps modulating the role of shallow cumulus in the climate system has been the subject of much debate. While it is generally perceived that aerosol perturbations increase drop concentrations (Twomey 1977), suppress the formation of rainfall (Warner 1968; Albrecht 1989), and increase the lifetime of clouds, there is scant observational evidence to support these suppositions and even some to the contrary (Small et al. 2009). Finescale modeling studies have suggested that the aerosol may play multiple roles—increasing cloudiness in aerosol-poor conditions but reducing it in aerosol-rich conditions (e.g., Ackerman et al. 2004; Xue et al. 2008). Some modeling studies have shown that by delaying the onset of precipitation, higher aerosol concentrations transport water vapor deeper into the inversion layer and support deeper cumulus clouds (McFarquhar and Wang 2006; Stevens and Seifert 2008). These deeper clouds offset the effect of aerosol-induced suppression of collision–coalescence on rainfall. Selfregulation of this kind led Stevens and Feingold (2009) to the conclusion that in many cases, the multitude of internal feedbacks between microphysics and dynamics might reduce the potential for the aerosol to influence the cloud system. The system may therefore be perceived as ‘‘buffered’’ or robust to aerosol perturbations. Motivated by the importance of warm cumulus clouds in the climate system, this study focuses on understanding how aerosol perturbations affect the statistical properties of a cloud system for relatively long simulations (.24 h) and scales of 25 km, which are large enough for mesoscale circulations (i.e., mesobeta circulations) to develop. By integrating through transients, the goal is to ascertain just how robust the system is to aerosol perturbations. The paper deviates somewhat from the typical studies of aerosol–cloud interactions that tend to focus on microphysical responses to aerosol influences via traditional constructs (Twomey 1977; Warner 1968) and instead attempts to incorporate a more boundary layer, cloud-system-centric view. While these aforementioned processes drive the responses that we aim to quantify, it is becoming increasingly clear that microphysical changes carry dynamical consequences. Stated differently, by changing cloud microphysical processes, the aerosol changes the thermodynamic environment in which a cumulus cloud system lives and evolves. This view of the cloud system opens up a much richer spectrum of physics and nuance to the aerosol– cloud precipitation problem, and, as will be shown here, helps to clarify the role of the aerosol, rather than complicate it.
VOLUME 69
With this perspective we present a set of numerical experiments that will show that for the sounding considered, aerosol perturbations influence the thermodynamic environment and cloud field development in such a way that many perturbed and unperturbed cloud field properties tend to converge after sufficient time has elapsed. The hypothesis posed is stated as follows: the vertical development of cloud fields responds to aerosolinduced changes in atmospheric stability so that instability (i.e., potential energy) is consumed at different rates; after a transient period, however, the systems evolve to a similar state. Sensitivity tests will show that this response is robust to a variety of different model assumptions/simulated conditions. Some preliminary thoughts on how the time scale for this convergence might vary are also proffered.
2. Large-eddy simulation (LES) This study relies on numerical simulations to fulfill its aim. The Advanced Research Weather Research and Forecasting Model (ARW-WRF, version 3.1.1) is adopted for the simulations. In the ARW-WRF, a high-order monotonic advection scheme (Wang et al. 2009) and a double-moment bulk microphysical scheme (Feingold et al. 1998; Wang and Feingold 2009) are used. The representation of cloud condensation nuclei (CCN) and cloud and raindrop size distributions follows Wang and Feingold (2009). Briefly, prognostic equations are solved for CCN number concentration, number and mass mixing ratio of cloud droplets, and number and mass mixing ratio of raindrops. Droplet activation is calculated based on predicted supersaturation, an assumed aerosol size distribution (lognormal), and composition (ammonium sulfate). After activation, the CCN are tracked through the drop population and regenerated upon evaporation. The CCN concentration is reduced by drop collision–coalescence and by surface rain. No CCN sources are applied. Changes in the CCN are represented by changes in concentration alone; sensitivity to the aerosol size distribution and composition are not considered. For the trade cumulus regime, neglect of size distribution and composition variability is unlikely to have a significant impact on the results to be presented (Feingold 2003).
3. Case description a. Numerical experiments A number of simulations are performed for an observed case of trade cumuli during the Rain in Cumulus over the Ocean (RICO) field experiment (Rauber et al.
DECEMBER 2012
LEE ET AL.
2007). The horizontal domain length is set at 25 km for both the east–west (x) and north–south (y) directions, while the vertical domain length is set at 4 km to cover the planetary boundary layer (PBL). Periodic boundary conditions are imposed on the horizontal domain. The horizontal grid length (Dx and Dy) is 100 m, while the vertical grid length Dz is 40 m. A thermodynamic sounding composited by the Global Energy and Water Cycle Experiment Cloud System Study (GCSS) boundary layer working group (http:// www.knmi.nl/samenw/rico) based on soundings from the RICO field experiment is used. The u component of the wind is easterly, decreasing linearly from about 10 m s21 at the surface to about 2 m s21 at the model top (4 km), while the y component is northerly and constant at 3.8 m s21 throughout. Water vapor mixing ratio r and potential temperature u are nearly constant from the surface to 740 m. From 740 m up, r decreases from 13.8 to 1.8 g kg21, while u increases from 297.9 to 317 K at the model top. Unless otherwise stated, simulations use the same forcings as those prescribed by GCSS. These are briefly described in section 3b. The first simulation, referred to as the control (C) run, adopts an initial background aerosol number concentration of 100 mg21 (equivalent to 100 cm23 at an air density of 1 kg m23) and this number is constant over the PBL, following the GCSS specification. The C run lasts 34 h. The first 6 h 10 min of the simulation period is considered to be ‘‘spinup’’ and is excluded from analysis (the influence of the duration of the spinup is explored in section 5g). To examine aerosol effects on the cloud system and its environment, a new simulation is spawned at the end of the spinup period with the background aerosol number concentration enhanced by a factor of 2.5, yielding 250 mg21 at grid points with no clouds. This repeated run uses wind, pressure, temperature, and humidity at each grid point over the entire domain from the control run as initial conditions. This repeated run is referred to as ‘‘the high-aerosol’’ or H run. A summary of simulations is shown in Table 1. In addition to the control and high-aerosol runs, supplementary simulations are listed in Table 1. These additional simulations will be described in the following sections and the appendixes.
b. Surface and large-scale forcings and radiation All forcings are based on the GCSS prescription. Bulk formulas are used to represent the surface sensible, latent heat, and momentum fluxes. The sensible and latent heat fluxes are proportional to the wind velocity in the atmosphere immediately above the surface (;10 m s21) and to the difference in temperature or water vapor
3609
mixing ratio between the surface and the atmosphere immediately above it. The momentum fluxes are proportional to wind velocity in the atmosphere immediately above the surface. The large-scale subsidence and horizontal advection of temperature and moisture are also prescribed by the GCSS specifications. Radiation processes are not calculated explicitly but represented by a temporally and spatially constant advective tendency on temperature (cooling) with no consideration of the diurnal cycle. The cooling profile was derived by vanZanten et al. (2011) from an offline radiation calculation for the initial temperature and humidity sounding. The cooling rate is 2 K day21 close to the surface and decreases to about 1 K day21 in the free troposphere. The large-scale subsidence and advection of moisture are a function of height only.
4. Results A series of simulations are carried out with the GCSS sounding to test the hypothesis proposed in the introduction.
a. Time series and vertical profiles Figure 1 shows the time series of the domain-averaged liquid water path (LWP), vertical velocity variance w9w9, cloud fraction (CF), cloud-top height, convective available potential energy (CAPE), buoyancy flux w9u9, y the PBL top height, and precipitation rate for the C and H runs. The PBL top height is defined as the level where the increase in u is greatest over a 40-m vertical interval in each of the grid columns (similar to Rauber et al. 2007). Up to about 14 h, the H run exhibits larger LWP, w9w9, CF, and cloud-top height than the C run (Figs. 1a–d and 1g). However, between about 14 and 20 h, LWP, w9w9, CF, and cloud-top height are larger in the C run. After about 20 h, differences in these variables are insignificant. Relative differences between the C and H runs are not sensitive to the magnitude of the initial perturbations to the sounding (appendix A). To facilitate discussion, the approximate times at which there is a marked shift in C versus H behavior are delineated. The period up to 14 h, between 14 and 20 h, and from 20 h to the end of the simulations are referred to as the first, second, and third periods, respectively. CAPE represents the gravitational potential energy stored in the system that can be converted to turbulent kinetic energy, which explains the similar time evolution among CAPE, w9w9, w9u9, y and the PBL height in Figs. 1b and 1e–g. These variables also show larger values in the H run up to about 14 h, smaller values in the H run
CrH
H-at-24-h
Identical to that in the extended C run but starts at the end of the extended spinup period at 24 h Second period (14–20 h)
Identical to that in the C run but starts at the end of the spinup period at 6 h 10 min 48 h
H
Extended C
34 h
Period
C
Simulations
85 mg21
100 mg21 at the start of the 48-h simulation but 57 mg21 at the end of the extended spinup period at 24 h 143 mg21
100 mg21 at the start of the 34-h simulation but 96 mg21 at the end of the spinup period 240 mg21
Initial background aerosol
Identical to the H run r at the end of the first period
Identical to the extended C run r at the end of the extended spinup period at 24 h
Identical to the extended C run u at the end of the extended spinup period at 24 h
Identical to the C run u at the end of the first period
Follows GCSS specification
Identical to the C run r at the end of the spinup period
Follows GCSS specification
Initial r
Follows GCSS specification
Identical to the C run u at the end of the spinup period
Follows GCSS specification
Initial u
TABLE 1. Summary of simulations.
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Surface moisture fluxes
Included
Included
Included
Included
Included
Cooling from rain evaporation
Not included
Not included
Not included
Not included
Not included
Explicit radiation
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Magnitude of random perturbation imposed on the initial potential temperature
3610 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 69
Second period (14–20 h)
Second period (14–20 h)
Second period (14–20 h)
First and second periods (06 h 10 min–20 h)
First and second periods (06 h 10 min–20 h)
First and second periods (06 h 10 min–20 h)
First and second periods (06 h 10 min–20 h)
CuH1.0–1.5 run
CuH -other
C-rain-evap-off
H-rain-evap-off run
C-high-rain
H-high-rain
Period
CuH
Simulations
240 mg21
96 mg21
240 mg21
96 mg21
85 mg21
85 mg21
85 mg21
Initial background aerosol
Identical to the C run u at the end of the spinup period
Identical to the H run u at the end of the first period Identical to the H run u at the end of the first period only for a layer between 1.0 and 1.5 km Identical to the H run u at the end of the first period for layers except for a layer between 1.0 and 1.5 km Identical to the C run u at the end of the spinup period Identical to the C run u at the end of the spinup period Identical to the C run u at the end of the spinup period
Initial u
Identical to the C run r at the end of the spinup period
Identical to the C run r at the end of the spinup period
Identical to the C run r at the end of the spinup period
Identical to the C run r at the end of the spinup period
Identical to the C run r at the end of the first period
Identical to the C run r at the end of the first period
Identical to the C run r at the end of the first period
Initial r
TABLE 1. (Continued)
Increased by a factor of 5 as compared to the GCSS specification Increased by a factor of 5 as compared to the GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Surface moisture fluxes
Included
Included
Excluded
Excluded
Included
Included
Included
Cooling from rain evaporation
Not included
Not included
Not included
Not included
Not included
Not included
Not included
Explicit radiation
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Magnitude of random perturbation imposed on the initial potential temperature
DECEMBER 2012 LEE ET AL.
3611
First and second periods (06 h 10 min–20 h)
First and second periods (06 h 10 min–20 h)
From 8 to 20 h
From 8 to 20 h
From 10 to 20 h
From 10 to 20 h
From 7 to 10 h
From 7 to 10 h
From 7 to 10 h
34 h
H-radiation
C-at-hr-8
H-at-hr-8
C-at-hr-10
H-at-hr-10
HuC
HrC
H-C-second
Pertx2 C
Period
C-radiation
Simulations
94 mg21 at 7 h and 92 mg21 at 8 h 100 mg21 at the start of the 34-h simulation but 95 mg21 at the end of the spinup period
237 mg21
237 mg21
100 mg21
100 mg21
100 mg21
100 mg21
240 mg21
96 mg21
Initial background aerosol
Follows GCSS specification
Identical to the H run u at 7 h
Identical to the C run u at the end of the spinup period Identical to the C run u at the end of the spinup period Identical to the C run u at 8 h Identical to the H run u at 8 h Identical to the C run u at 10 h Identical to the H run u at 10 h u in the C run first at 7 h and then at 8 h Identical to the H run u at 7 h
Initial u
Follows GCSS specification
r in the C run first at 7 h and then at 8 h Identical to the H run r at 7 h
Identical to the C run r at 8 h Identical to the H run r at 8 h Identical to the C run r at 10 h Identical to the H run r at 10 h Identical to the H run r at 7 h
Identical to the C run r at the end of the spinup period
Identical to the C run r at the end of the spinup period
Initial r
TABLE 1. (Continued)
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification Follows GCSS specification Follows GCSS specification Follows GCSS specification Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Surface moisture fluxes
Included
Included
Included
Included
Included
Included
Included
Included
Included
Included
Cooling from rain evaporation
Not included
Not included
Not included
Not included
Not included
Not included
Not included
Not included
Included
Included
Explicit radiation
Magnitude increased by a factor of 2 compared to that in the C run
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification Follows GCSS specification Follows GCSS specification Follows GCSS specification Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Magnitude of random perturbation imposed on the initial potential temperature
3612 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 69
Identical to that in the Pert/2 C run but starts at the end of the spinup period at 6 h 10 min
Pert/2 H
Pert/2 C
Identical to that in the Pertx2 C run but starts at the end of the spinup period at 6 h 10 min 34 h
Period
Pertx2 H
Simulations
100 mg21 at the start of the 34-h simulation but 97 mg21 at the end of the spinup period 240 mg21
240 mg21
Initial background aerosol
Identical to the Pert/2 C run u at the end of the spinup period
Follows GCSS specification
Identical to the Pertx2 C run u at the end of the spinup period
Initial u
Identical to the Pert/2 C run r at the end of the spinup period
Follows GCSS specification
Identical to the Pertx2 C run r at the end of the spinup period
Initial r
TABLE 1. (Continued)
Follows GCSS specification
Follows GCSS specification
Follows GCSS specification
Surface moisture fluxes
Included
Included
Included
Cooling from rain evaporation
Not included
Not included
Not included
Explicit radiation
Magnitude decreased by a factor of 2 compared to that in the C run Magnitude decreased by a factor of 2 compared to that in the H run
Magnitude increased by a factor of 2 compared to that in the H run
Magnitude of random perturbation imposed on the initial potential temperature
DECEMBER 2012 LEE ET AL.
3613
3614
JOURNAL OF THE ATMOSPHERIC SCIENCES
FIG. 1. Time series of the domain-averaged (a) LWP, (b) w9w9, (c) CF, (d) cloud-top height, (e) CAPE, (f) w9u9, y (g) PBL-top height, and (h) precipitation rate. Vertical blue lines are added to mark start and end times of stages and periods.
VOLUME 69
DECEMBER 2012
LEE ET AL.
between about 14 and 20 h, and negligible differences after about 20 h. Figures 2 and 3 show the vertical profiles of liquid water content (LWC) and w9w9, respectively, averaged over the domain for each of the three periods and for the whole simulation period. The w9w9 and LWC profiles (Figs. 2 and 3) show similar vertical distributions. In the first period, both w9w9 and LWC are larger in the H run throughout the cloud layer, while in the second period both w9w9 and LWC show larger values in the C run, particularly in a layer between about 1.0 and 1.8 km. In the third period, differences in w9w9 and LWC are small compared to those in the first and second periods, consistent with similar time evolutions of LWP and w9w9 (Figs. 1a and 1b). Hence, one can assume with good confidence that the larger w9u9y induces larger w9w9, which in turn induces larger LWP, CF, and cloud-top height. The rain rate deserves special attention because of the influence of the aerosol on rain production. During the first period (6 h 10 min–14 h), the mean rain-rate differences are small as a result of higher LWP in the H run, offsetting aerosol suppression of collision–coalescence (Fig. 1h). During the second period (14–20 h), the largest precipitation difference is simulated, as enhancements in LWP in the C run work in unison with more efficient collision–coalescence. During the third period (after 20 h), precipitation differences are reduced but nevertheless are distinctly larger for the C run. Further analysis of the precipitation fields is deferred to the end of section 4.
b. First period To explain the reasons for the differences between the C and H runs (Fig. 1), u and r are examined. First, the CAPE difference (Fig. 1e) during the first period is examined. For this, the first period is subdivided into three stages, and u and r are considered in each of the stages. The first, second, and third stages correspond to periods up to 8 h, between hours 8 and 10, and between hours 10 and 14, respectively. The choice of the time intervals for the subdivisions is primarily based on the evolution of CAPE, but it is supported by a similar temporal evolution in other fields (Fig. 1). During the first and second stages, CAPE in the H run increases rapidly and then levels off at a magnitude significantly higher than that for the C run (Fig. 1e); during the third stage, it decreases. In contrast, CAPE in the C run increases much more slowly (Fig. 1e), and it is only at about 14 h, or around the beginning of the second period, that it exceeds the value of the (descending) CAPE in the H simulation. To understand the increasing and then steady CAPE differences between the C and H runs during the first
3615
and second stages, the vertical distributions of the domain-averaged r and u at 7 and 8 h are shown in Figs. 4a–d. Note that the distributions at 7 h are around the middle of the time period of the increasing CAPE in the H run. A smaller u develops in a layer between about 0.4 and 1.0 km in the H run and the difference in u between the H and C runs reaches about 0.2 K in the layer (Figs. 4b and 4d). This layer between about 0.4 and 1.0 km corresponds to the low-level cloud layer and is henceforth referred to as DZl. Sensitivity tests show that this smaller u at hours 7 and 8 increases CAPE in the H run during the first stage and promotes the larger CAPE during the second stage. These sensitivity tests are described in detail in appendix B. In addition, r is larger in the H run in DZl at both hours 7 and 8 and the difference in r between the H and C runs reaches about 0.2 g kg21 in the layer (Figs. 4a and 4c). Additional simulations demonstrate that it is the lower u, not the higher r, that accounts for most of the CAPE difference shown in Fig. 1e during the first and second periods (see appendix B). Does subcloud evaporation of rain play a role in explaining the evolution of C and H simulations? To address this, the C and H runs are repeated but with evaporative cooling associated with subcloud rain evaporation turned off. These are referred to as ‘‘the C-rain-evap-off run’’ and ‘‘the H-rain-evap-off run,’’ respectively. As in the C and H runs, the H-rain-evapoff run has larger LWC than the C-rain-evap-off run in the first period (Fig. 2e). Cloud-base height increases by about 100 m with subcloud evaporation of rain turned off. Further tests with rain evaporation turned off throughout the domain again produce similar trends to the C and H runs (not shown), reinforcing these results. Thus, the differences in latent heat distribution above cloud base associated with cloud liquid condensation and evaporation play a much more important role in the evolution of the H and C cases than does latent heat redistribution from rain evaporation below cloud base. Figure 5 shows the vertical distribution of the timeand area-averaged net condensation rate over each of the three stages during the first period for the C and H simulations. Note that negative net condensation represents net evaporation. We focus on the net condensation rate because calculations show that the net heating due to the net condensation is approximately one order of magnitude larger than that due to advection/ turbulent diffusion, and because in these simulations, radiative cooling is applied equally to both simulations. Figures 4 and 5 are useful for understanding the evolution of the cloud system through the three stages.
3616
JOURNAL OF THE ATMOSPHERIC SCIENCES
FIG. 2. Vertical profiles of the time- and domain-averaged LWC over (a),(e) first, (b),(f) second, (c) third, and (d) entire simulation periods.
VOLUME 69
DECEMBER 2012
3617
LEE ET AL.
FIG. 3. Vertical profiles of the time- and domain-averaged w9w9 over (a) first, (b) second, (c) third, and (d) entire simulation periods.
1) STAGE 1 (6
H
10
MIN
, T , 8 H)
Compared to the C run, the H run exhibits a lower net condensation rate in DZl (Fig. 5a), causing less heating in this layer (Figs. 4b and 4d) and thus smaller u, higher instability and CAPE, and lower convection inhibition (CIN). The time- and domain-averaged CIN over stage 1 is 5 and 7 J kg21 in the H and C runs, respectively. The lower net condensation rate is associated with larger r (Figs. 4a and 4c). The weaker heating and larger r in DZl in the H run are also simulated in the H-rain-evap-off run, confirming the negligible role of subcloud evaporation of precipitation (Figs. 4a and 4b). The relative cooling and moistening in the layer is due to two factors: (i) Smaller and more numerous droplets in the H run contribute to more efficient cooling and moistening of the air (Xue and Feingold 2006; Hill et al. 2009). The time-averaged cloud droplet number concentration and droplet radius over cloudy areas for C (H) are
25 (60) mg21 and 7.5 (5.4) mm. (ii) Similar amounts of condensed water but smaller and more numerous clouds in the H run have larger surface-to-volume ratios, leading to an increase in their susceptibility to entrainment of dry air, further enhancing cooling and moistening of the air (Xue and Feingold 2006). During the first stage, the time-averaged number and the time- and domain-averaged equivalent diameter of clouds are 338 (451) and 332 (297) m in the C (H) run, respectively.
2) STAGE 2 (8
H
, T , 10
H)
For the H run, the more efficient cooling generates larger instability, CAPE, and w9u9, y resulting in clouds with higher cloud tops. The elevated cloud-top heights can be seen in the elevated region of net evaporation in the H run (Fig. 5b) associated with convective divergence around cloud tops. However, there are
3618
JOURNAL OF THE ATMOSPHERIC SCIENCES
VOLUME 69
FIG. 4. Vertical profiles of the domain-averaged water vapor mixing ratio at (a) 7, (c) 8, (e) 10, and (g) 14 h. (b),(d), (f),(h) As in (a),(c),(e),(f), respectively, but for u . The top abscissa shows the differences (H minus C) in r and u fields.
a number of consequences of this efficient evaporation: first, a stronger humidification in a layer between about 1.4 and 1.8 km, which counters the accelerated evaporation; and second, the deeper convection expends more energy. In contrast, the slower rate of evaporation in the C run produces lower cloud-top heights and less CAPE consumption. While CAPE, w9w9, and w9u9y continue to rise steadily through this stage in the C run, they tend to level off in the H run, as seen in the Figs. 1b, 1e, and 1f. Moreover, the integrated net cooling in a layer
between about 1.0 and 1.5 km, which is similar in magnitude to that in the H run, generates more concentrated cooling (Fig. 5b) because of the lower cloudtop heights, resulting in destabilization of this layer (Fig. 4f). This layer between about 1.0 and 1.5 km corresponds to the midlevel cloud layer (maximum cloud-top heights are about 2.0 km for both the C and H runs) and is referred to as DZm. Again, the trends in net condensation rate in the C-rain-evap-off and H-rain-evap-off runs during stage 2 (Fig. 5d) are much
DECEMBER 2012
LEE ET AL.
3619
FIG. 4. (Continued)
the same as those for the C and H runs (Fig. 5b), demonstrating that more concentrated cooling in DZm in the C run is not influenced by the presence of subcloud rain evaporation.
3) STAGE 3 (10
H
, T , 14 H)
For the C run, the smaller u and thus larger instability in DZm sustains a steady rise in CAPE and w9u9y (Figs. 1e and 1f). Evaporative cooling (negative net condensation) for DZm continues to strengthen versus that in the H run (cf. Figs. 5b and 5c) and is expressed by increasing
(C vs H runs) u differences in the layer (cf. Figs. 4f and 4h). The general trend is for the destabilization of the C run and the stabilization of the H run. The H run experiences declining CAPE, and, by the end of stage 3, weakening w9w9 and w9u9. y Cloud fields such as LWP, CF, and cloud-top height have more inertia to these changes in stability and their decrease is delayed until the second period. As in stage 2, the trends in net condensation rate in the C-rain-evap-off and H-rain-evap-off runs during stage 3 (Fig. 5e) are similar to those for the C and H runs (Fig. 5c).
3620
JOURNAL OF THE ATMOSPHERIC SCIENCES
VOLUME 69
FIG. 5. Vertical profiles of the time- and domain-averaged net condensation rate over (a) first stage, (b),(d) second stage, (c),(e) third stage, and (f),(g) second and third periods. For effective visualization of the vertical distribution of net evaporation rate (i.e., net negative condensation) above ;1 km, different scales in the x axis are used for positive and negative net condensation rates. The top abscissa quantifies the net condensation rate differences (H minus C).
c. Second period INFLUENCE OF THERMODYNAMIC PROFILES The u and r differences between the C and H runs at the end of the first period (14 h) in Figs. 4g and 4h are largest in a layer between about 1.5 and 2 km. These are a result of moistening and cooling associated with deeper convection in the H run during the first period. This might be expected to favor subsequent clouds with higher cloud top-heights, larger depth, and thus larger LWC and LWP in the H run. However, the converse is true: during the second period, C-run clouds have, on average, larger LWC, LWP, and higher cloud-top heights (Figs. 1a, 1d, and 2b). This suggests that it is not differences in the layer between about 1.5 and 2 km, but
those in other layers that result in the LWC, LWP, and cloud-top height differences in the second period. As seen in Fig. 2b, larger LWP associated with larger CAPE and thus w9w9 (Figs. 1a, 1b, 1e, and 3b) for the C run during the second period (14–20 h) are mainly due to larger LWC in the layer between about 1.0 and 1.8 km. This layer corresponds to a mid- and upper-level cloud layer and is referred to as DZmu. To identify differences in u and r, which lead to the LWC differences in DZmu, the control run is repeated from the beginning to the end of the second period with different u and r conditions. Figure 6 depicts the vertical distribution of the time- and domain-averaged LWC for the repeated runs, as well as the C and H runs over the second period. The first of these repeated runs is identical to the control
DECEMBER 2012
LEE ET AL.
3621
FIG. 5. (Continued)
run from the beginning to the end of the second period, except that the initial r corresponds to that from the H run at each grid point at the last time step of the first period. This simulation is referred to as ‘‘CrH’’ (control with r associated with the H run). The larger LWC in DZmu is still simulated in the CrH run, similar to the C run (Fig. 6). In addition, as in the C run, the time- and domain-averaged CAPE is higher in the CrH run than in the H run (Table 2). It is also notable that the vertical distribution of differences in LWC between the CrH and H runs in DZmu is similar to that between the C and H runs. Thus, the CAPE and LWC differences between C and H in DZmu are not caused by differences in r fields.
Next, the C run is repeated from the beginning to the end of the second period, except with the u field from the H run at each grid point at the last time step of the first period, referred to as ‘‘CuH’’. The LWC increase in DZmu disappears and is similar to that in the H run. The time- and domain-averaged CAPE are similar in the CuH and H runs, in contrast to the larger difference between the C and H runs (Table 2). Thus, it is the u difference established during the first period that leads to the CAPE and LWC differences between the C and H runs in DZmu in the second period. The u profiles in Fig. 4h indicate that there is generally a larger instability in DZm at the end of the first period in
3622
JOURNAL OF THE ATMOSPHERIC SCIENCES
VOLUME 69
TABLE 2. Time- and domain-averaged CAPE for periods and stages.
Simulations C
H
CrH CuH CuH1.0–1.5 run CuH -other HuC HrC H-C-second
FIG. 6. Vertical profiles of the time- and domain-averaged LWC over the second period.
the C run. This is due to generally smaller u in the C run. In other layers, u tends to be lower in the H run. Hence, it is likely that the stability differences in other layers do not contribute to the larger LWC in the C run in the second period. To confirm this, the C run is repeated in the same manner as in the CuH but with the initial u from the H run at the last time step of the first period only for grid points in the layer DZm; the initial u values in layers other than DZm in this repeated run are the same as in the C run. This repeated run is referred to as ‘‘the CuH1.0–1.5 run.’’ A comparison among the CuH1.0–1.5 run, the CuH run, and the H run shows that CAPE and the LWC profiles in the CuH1.0–1.5 run are closer to those in the H run than in the CuH run (Table 2 and Fig. 6). Next, the C run is repeated in the same manner as in the CuH1.0–1.5 run but with the initial u from the H run at the last time step of the first period at all grid points except for those in DZm; that is, in this repeated run, the initial u is the same as in the C run in DZm, and in other layers it is the same as in the H run. This repeated run is referred to as ‘‘the CuH-other run.’’ The larger CAPE and LWC in DZmu are simulated in the CuH -other run as in the C run, and the differences in LWC between the CuH -other run and the H run are similar to those between the C and H runs (Table 2 and Fig. 6). These CuH1.0–1.5 and CuH -other runs demonstrate that it is the larger instability in DZm that causes the larger CAPE and LWC in DZmu in the C run during the second period.
Time period
Average CAPE (J kg21)
First stage (06 h 10 min–8 h) Second stage (8–10 h) Second period (14–20 h) First stage (06 h 10 min–8 h) Second stage (8–10 h) Second period (14–20 h) Second period (14–20 h) Second period (14–20 h) Second period (14–20 h) Second period (14–20 h) First stage (06 h 10 min–8 h) Second stage (8–10 h) Second stage (8–10 h) Second stage (8–10 h)
56 58 68 63 72 62 67 61 62 67 59 57 70 67
As in the C and H runs, the H-rain-evap-off run has smaller LWC than the C-rain-evap-off run in the second period (Fig. 2f). This demonstrates that during both the first and second periods, subcloud latent heat distribution from rain evaporation does not play an important role in the evolution of the H and C cases, and that the evolution is driven by above-cloud-base latent heat distribution through phase changes.
d. Third period and overview of the three periods It is notable that aerosol-initiated deviations in cloud and environmental properties decrease as time progresses; C versus H differences in mean CAPE, w9u9, y w9w9, LWP, and most other variables become small, even negligible, during the third period relative to their differences during the first and second periods (Fig. 1) (two exceptions to this convergence are the rain rate, Fig. 1h, and cloud optical depth, section 5d). As noted above, precipitation differences between C and H are closely related to LWP and drop concentration Nd differences. The mean precipitation rate is always higher for the C simulation, only marginally during period 1, and most significantly during period 2. To further examine the differences in precipitation, the frequency distributions of precipitation rates for periods 1, 2, and 3 are shown in Figs. 7a–c. During period 1, the H run is characterized by more frequent weak (,0.1 mm day21) rain, but the reverse is true in periods 2 and 3. In all three periods, the C run exhibits more frequent moderate to strong rain rates, particularly in period 2 when C rain is significantly larger. Note that the (very rare) strongest rain rates are always associated with a subset of the H-run clouds (see section 5e for further discussion).
DECEMBER 2012
LEE ET AL.
3623
FIG. 7. Distributions of precipitation rate frequency for the (a) first, (b) second, and (c) third periods.
5. Discussion Our contention, supported by the detailed analysis above, is that differences in the aerosol concentration create differences in environmental instability, which lead to cloud fields with different properties; these, in turn, remove instability at a rate commensurate with the buildup of instability. In the system perturbed by aerosol, low-level cooling allows CAPE to increase more rapidly, which leads to the earlier appearance of its peak in period 1 (Fig. 1e). But a consequence is that the cloud field consumes CAPE at a faster rate (Nober and Graf 2005). In the unperturbed system instability develops more slowly, but the attendant shallower cloud field eventually produces an instability at DZm that allows clouds to deepen (on average) more than their perturbed counterparts. The unperturbed and perturbed systems follow different temporal cycles, but eventually they converge to a similar state. Thus, aerosol perturbations generate dynamical responses that result in similar environmental and cloud field properties after a sufficiently long time period. It has previously been shown that increases in aerosol tend to speed up the cycle of shallow cumulus
convection by generating stronger instabilities earlier on and subsequently removing them at a faster rate. Jiang et al. (2009) calculated the birth and death rates of convective cells for the same RICO GCSS case (albeit using a different model) and showed that aerosol-perturbed systems have faster birth/death rates. Analysis of cloud size distributions in section 4b also supports earlier results from different models and different trade cumulus soundings that aerosolperturbed cloud systems comprise larger numbers of smaller clouds that have shorter lifetimes (Xue and Feingold 2006; Jiang et al. 2006). By perturbing thermodynamic profiles, the aerosol therefore appears to be able to change the rate at which instability is produced and consumed. The results raise a number of other interesting issues and questions; these are addressed below.
a. Role of precipitation Note that the results mentioned above are for cumulus clouds with low surface precipitation. However, it is likely that more strongly precipitating clouds would have a stronger influence on subcloud stability, and in this case precipitation may play a more important role in
3624
JOURNAL OF THE ATMOSPHERIC SCIENCES
VOLUME 69
FIG. 8. Time series of the domain-averaged (a) LWP, (b) w9w9, and (c) CAPE. Vertical blue lines are added to mark starting and ending time of stages and periods.
the cloud and LWC evolution and the effect of aerosol thereon. To investigate this, the C and H runs are repeated for the first and second periods, but with surface moisture fluxes increased by a factor of 5. These repeated runs are called C-high-rain and H-high-rain runs and have time- and domain-averaged rain rates that are 8 and 5 times larger than those in the C and H runs, respectively. Thus, the increased surface fluxes have a disproportionately stronger influence on the rain rate in the low aerosol simulation. This causes a much larger relative increase in the stabilization of the subcloud layer in the C-high-rain run compared to the C run than that between the H-high-rain run and the H run. Consequently, the differences in the midlevel cooling for these high-rain runs are much smaller than for the standard C and H runs during the first period, which in turn leads to significantly reduced differences (about 5 times smaller) in CAPE, w9w9, LWC, and LWP compared to those between the C and H simulations during the second period (Figs. 2b and 8). Stronger precipitation more effectively stabilizes the profile by heating the cloud layer and cooling the subcloud layer, and therefore it appears to hasten the convergence of the C and H simulations.
b. Role of radiation To examine the effect of the diurnal cycle of radiation, the C and H runs are repeated with radiation calculated explicitly using the National Center for Atmospheric Research (NCAR) Community Atmosphere Model
(CAM) radiation scheme for the first and second periods (recall that in the standard C and H runs, a constant cooling is applied to both C and H runs in lieu of radiation) These repeated runs are referred to as the C-radiation and H-radiation runs, respectively. In these repeated runs, the beginning and end of the simulations are assumed to correspond to 0610 and 2000 local solar time (LST) on 16 December 2004, since there are no designated LSTs for simulations in the GCSS specifications. With these assumed LSTs, the simulations capture the diurnal variation of radiation and its impact on clouds from sunrise through sunset to night. In the C and H runs, there are no twoway interactions between radiation and clouds, nor is there a diurnal cycle, whereas in the C-radiation and H-radiation runs, the diurnal cycle of incident shortwave radiation on clouds and two-way interactions among clouds, and shortwave and longwave radiation are taken into account (e.g., Wang and McFarquhar 2008a). Figures 2a, 2b, and 8 show that the qualitative nature of C versus H differences in CAPE, w9w9, LWP, and LWC evolution does not depend on the details of the radiation, at least not during the first and second periods. By 20 h, the mean CAPE fields have converged, although interestingly differences in LWP and w9w9 remain substantial, whereas they were significantly diminished in the standard C and H runs. This raises questions about the characteristic time scales for relaxation of the aerosol-perturbed system, which we will consider in section 5f.
DECEMBER 2012
LEE ET AL.
c. Simulations with a reset of aerosol Results have shown that aerosol-triggered differences in u or thermal instability (mediated by clouds) cause differences in the development of subsequent clouds. The aerosol triggers (i) low-level (about 0.4–1 km) cooling in the H run during stage 1 of period 1, which results in clouds with higher top heights; and (ii) concentrated cooling at midlevel (about 1–1.5 km) in the C run during stages 2 and 3 of period 1. This concentrated cooling develops clouds with higher LWP in the C run during the second period. Hence, it is possible that once aerosol-induced changes in environmental conditions are established, the system has ‘‘memory’’ of the aerosol perturbation. In particular, the increase in CAPE in the C run at the beginning of the second period (14 h) is likely to be caused by aerosol-related differences in environmental conditions that are set up long before the beginning of the second period. Two sets of the repeated C and H runs are performed to explore the role of the timing and duration of the aerosol perturbation in determining the convergence of the fields. The domain-averaged aerosol concentrations for the C and H runs are 91 and 235 mg21 at 8 h (i.e., the beginning of the second stage of period 1). To investigate the extent to which the system retains memory of the aerosol perturbation, aerosol number concentrations for the C and H runs are reset to 100 mg21 at 8 h at all grid points. This implies a brief duration of perturbation (from 6 h 10 min to 8 h), followed by a small increase in the CCN concentration of the C run (912100 mg21), and a more significant decrease from 235 to 100 mg21 for the H run. These repeated C and H runs are referred to as the C-at-hr-8 run and the H-at-hr-8 runs, respectively, and are extended to the end of the second period at 20 h. In these runs, from the second stage on, differences in cloud development between the C and H runs are caused only by differences in the environment induced by the aerosol perturbation during the first stage. As expected, CAPE, w9w9, LWP, CF, and cloudtop height in the C-at-hr-8 run (Fig. 9), generally follow those of C, although there are some differences during period 2. A larger difference can be seen in the H-at-hr-8 run (vs H), as a result of the more significant reduction in CCN concentration. At 8 h, which happens to be the time at which CAPE peaks in the H run, the switch to a lower aerosol concentration means that the faster convective turnover associated with high aerosol (Jiang et al. 2009) is no longer able to remove the instability as efficiently. As a result, H-at-hr-8 sustains CAPE, w9w9, LWP, CF, and cloud-top height at higher values than those in H. The convergence between C-at-hr-8 and H-at-hr-8 is still apparent, but it occurs along a different pathway.
3625
A second set of repeated runs, referred to as C-at-hr-10 and H-at-hr-10, is initiated at 10 h. In these runs the homogenization of the aerosol fields is from (domain averaged) aerosol concentrations of 89 (C) and 226 (H) to 100 mg21 (Table 1). Here too, the C-at-hr-10 simulation behaves much like C (Fig. 9); however, now the H-at-hr-10 fields are similar to the H fields. Similar to the H run, the longer exposure to high CCN concentrations establishes the low-level instability in DZl, which supports the more turbulent and thicker clouds. The depletion of CAPE, which starts soon after the switch to 100 mg21 aerosol concentrations and therefore still carries memory of the high aerosol, continues through the second period, resulting in converged cloud fields that are similar in the mean to their C and H counterparts. The repeated simulations in this section indicate that, if the aerosol perturbation is sufficiently long-lived, the accompanying changes in environmental conditions will influence cloud field evolution more so than the direct effect of the aerosol difference. Just how large a timeintegrated perturbation is required for this to hold true will likely be case dependent.
d. Aerosol influence on cloud optical depth Although not the focus of this study, it is instructive to consider the influence of aerosol perturbations on cloud optical depth t, which is closely related to albedo. In keeping with the relative importance of indirect effects in this paper, we consider simple calculations of domainaverage t based on the second moment of the drop size distributions (constrained by the bimodal lognormal function) and at a visible wavelength; we avoid higherorder calculations such as a posteriori three-dimensional radiative transfer, which would provide a rigorous assessment of aerosol indirect forcing. A time series of domain-average t for H and C simulations (Fig. 10) shows, not unexpectedly, that the t response is a modulation of the LWP response, because to first-order t ; LWP5/6 Nd1/3 . During the first period when increases in both LWP and Nd work in unison, t for the H run is significantly larger (55%) than for the C run. During period 2 the C run experiences significantly higher LWP, which almost makes up for the Nd differential, leading to a small (5%) increase in t for the H run. During period 3 when cloud fraction and LWP have converged, the average enhancement in t is about 30%; the ratio of Nd between the H and C simulations is 81 mg21/31 mg21 5 2.6, so that enhancement in t should be (2.6)1/3 5 1.38, or 38%. Therefore, only about 80% of the expected (Twomey 1977) increase is realized. Over the 34 h simulation, increases in LWP between the simulations are also approximately balanced by
3626
JOURNAL OF THE ATMOSPHERIC SCIENCES
VOLUME 69
FIG. 9. As in Figs. 1a–e, but for results from simulations with a reset of aerosol, as described in section 5c, together with those from the C and H runs.
decreases. The ratio of Nd between H and C is 75 mg21/ 32 mg21 5 2.34 and the expected increase in t is therefore 33%. Direct calculation indicates an increase in t ; 25%, so that again t increases are below those expected. Further analysis points to two possible sources of this discrepancy: (i) Horizontal variability in the local cloud optical properties within the 25 km domain. Idealized calculations that maintain the correct average LWP and Nd between the H and C simulations
but allow for different horizontal rearrangements of either LWP or Nd in the H and C cloud fields can be shown to result in both decreases (as in our case) or increases in the expected t response based on the average Nd differences as a result of the nonlinearity of t dependence on Nd. (ii) Inhomogeneities in the vertical distribution of cloud microphysical properties can also cause biases and nonlinear responses (Brenguier et al. 2000).
FIG. 10. Time series of domain-averaged cloud optical depth for the C and H simulations.
DECEMBER 2012
LEE ET AL.
3627
e. Aerosol influence on boundary layer depth and precipitation rate Previous work has suggested that the aerosol can influence the distribution of cumulus cloud-top heights (Xue and Feingold 2006; Wang and McFarquhar 2008b; Stevens and Seifert 2008; Jiang et al. 2009). The general trend is for an increase in boundary layer depth or cumulus cloud-top height with increasing aerosol as a result of elevation of the initial height at which precipitation forms and more efficient preconditioning of free-tropospheric air (e.g., Stevens 2007; Stevens and Seifert 2008; Nuijens et al. 2009). In this work, the increase in boundary layer depth with increasing aerosol occurs in the first period but is not a robust feature throughout the simulation (Fig. 1g). Note also that the changes in mean height do not necessarily reflect the changes in the maxima. For example, although Fig. 1d indicates a higher mean cloud-top height for the C run, during the second period, Fig. 2b clearly shows that during this same period, the deepest clouds are associated with the H simulation. Similar results can be seen in Xue and Feingold (2006, their Figs. 2d and 6f). Frequency distributions of cloud-top height for C and H runs (Fig. 11) show that in general the time series means in Fig. 1 are a faithful representation of the relative strength of convection from periods 1 through 3. However, as seen in Fig. 11, the H run always produces the highest cloud-top heights. This is supported by frequency distribution analysis of CAPE, LWP, and w9w9 (not shown), all of which exhibit a tail of more frequently occurring extrema associated with the H simulation, much like those in Fig. 11. Thus, the H simulation does produce a very small number of the deepest, most vigorous, and most strongly precipitating clouds (Figs. 7b, 7c, and 7d) in all periods, even when the mean cloud-top height is lower than in the C run.
f. Some thoughts on time scales for convergence There appears to be a general dearth of discussion in the literature on time scales for equilibration in response to an aerosol perturbation, particularly for shallow cumulus. Most early large-eddy simulations of boundary layer clouds were of short duration and limited to small domains, so that the studies were implicitly concerned with transients. Enhanced computational power has allowed simulations of much longer duration but only a handful of those have considered the influence of aerosol perturbations (e.g., Stevens and Seifert 2008; Wang and Feingold 2009; Bretherton et al. 2010; Kazil et al. 2011). Bretherton et al. (2010) discuss a 1–2-day thermodynamic adjustment period for stratocumulus
FIG. 11. Frequency distribution of cloud-top heights for C and H simulations for (a) first, (b) second, and (c) third periods. Note the persistent appearance of more frequently occurring high cloud-top heights in the H simulation, even when the mean values are similar (cf. Fig. 1d).
and focus on aerosol-induced (among others) evolution to different steady states. Whether similar aerosol-related bifurcations exist for the trade cumulus regime is unclear. While a rigorous study of factors influencing equilibration times must await further study, our assessment
3628
JOURNAL OF THE ATMOSPHERIC SCIENCES
VOLUME 69
FIG. 12. Time series of (a) LWP and (b) CAPE for extended C (solid) and H-at-24-h (dashed) runs from 24 to 48 h.
from the current study is that in the absence of largescale forcing, important controlling parameters would be (perhaps obviously) the timing and duration of the aerosol perturbation (Fig. 9) as well as the magnitude of precipitation (Fig. 8).
g. Simulations with a different spinup period The possibility that the difference between the C and H runs is particular to the choice of spinup period (6 h 10 min) is now raised. The end of the spinup period occurs long before the stabilization of the cloud system at 20 h (Fig. 1). How might results differ if the bifurcation into two different aerosol concentrations were to occur at a time when the system had already stabilized? To answer this question, the spinup period is extended to 24 h and the C run is extended to 48 h; this extended C run is referred to as ‘‘the extended C run.’’ As in the earlier H run, a new simulation is spawned that uses the wind, pressure, temperature, and humidity fields at each grid point from the extended C run as initial conditions at 24 h. The only difference is that this simulation has a background aerosol number concentration enhanced by a factor of 2.5, yielding 143 mg21 at grid points with no clouds. This simulation is referred to as ‘‘the H-at-24-hr run.’’ Figure 12 shows the LWP and CAPE evolution for the extended C and H-at-24-h runs from 24 h on. From 24 to 34 h LWP and CAPE are higher in the H-at-24-h run than in the extended C run. The situation reverses between 34 and 42 h, but after about 42 h, the fields stabilize, in qualitative agreement with Fig. 1. It is notable, however, that the magnitude of differences in LWP and CAPE from 24 to 42 h is smaller
than those in Fig. 1. It is also notable that the higher (lower) LWP and CAPE in the H-at-24 h run than in the extended C run is maintained for about 10 (8) from 24 (34) to 34 (42) h. These time periods of 10 and 8 h are longer than the duration of the first and second periods, respectively. Thus, although the magnitude of the LWP and CAPE differences is smaller, it takes longer to reach stabilization of LWP and CAPE compared to the simulations with the spinup time of 6 h 10 min. Nevertheless, the qualitative agreement with Fig. 1 provides confidence that the underlying mechanisms described herein are robust.
6. Summary and conclusions The prevalence of weakly precipitating trade cumulus and their importance for the climate system (e.g., Bony and Dufresne 2005) provides strong motivation for detailed examination of trade cumulus clouds and their response to changes in aerosol. We have examined here the influence of aerosol perturbations on long-duration (;30 h) trade cumulus simulations, specifically testing the hypothesis that the changes in thermodynamic profiles triggered by an aerosol perturbation cause changes in cloud vertical development that act to remove the differences. This concept is congruent with the buffered aerosol–cloud system discussed in Feingold and Siebert (2009), Stevens and Feingold (2009), and Lee and Feingold (2010). The aforementioned studies hypothesized that a system tends to reduce an aerosol impact by modifying environmental conditions or microphysical pathways.
DECEMBER 2012
LEE ET AL.
In the current paper, we have taken a detailed look at the sequence of events associated with the boundary layer response to an aerosol perturbation using the RICO GCSS case study as a basis for numerical experimentation. It is shown that increasing the aerosol concentration generates smaller droplets that evaporate more efficiently. This results in smaller and more numerous clouds with a larger surface-to-volume ratio, which through a feedback mechanism increases the susceptibility of clouds to entrainment of dry air (Small et al. 2009). Increases in aerosol therefore accelerate drop evaporation, resulting in stronger cooling and thus increased instability in the lower part of clouds near the beginning of simulations. The larger instability in the lower part of clouds and associated larger CAPE develop clouds with higher top heights, and higher LWC and LWP in the aerosol-perturbed simulation in the early period of simulations. These perturbed clouds with higher top heights result in elevated (in height) convective divergence and thus higher locations of peak evaporative cooling. Numerous sensitivity tests point to aerosol-induced perturbations in potential temperature dominating aerosol-induced perturbations in water vapor in terms of the response of the cloud system. In contrast, unperturbed simulations are characterized by a lower location of the peak in net cooling, which over time generates instability at altitudes that correspond to the middle part of perturbed clouds. The larger instability and associated larger buoyancy flux and CAPE at these altitudes invigorate the clouds and increase their LWP during the middle period of simulations. Thus, the system tends to adjust itself by modifying its environment in such a way that the effects of the initial aerosol perturbation are countered. This adjustment tends to occur after about 14 h, which cautions against short-term simulations when attempting to ascertain the effects of aerosol perturbations. After about 20 h, the aerosol-perturbed system stabilizes to a similar state to that in the unperturbed simulation. The result is that in the mean there are very small deviations in most cloud fields when considered over the entire simulation period. This result holds for different model initializations (Fig. A1) and different spinup periods (cf. Figs. 1 and 12), providing confidence in the underlying physical mechanisms described herein. One exception is the precipitation field. Although the time- and domain-averaged precipitation decreases with increasing aerosol, the precipitation response to aerosol is strongly modulated by cloud LWP or depth. The deepest clouds occur in the perturbed simulations, most noticeably in the first period (Fig. 11a), but also as rare
3629
events in periods 2 and 3 (Figs. 11b and 11c). This ‘‘tail’’ of deep clouds produces the largest (rare) precipitation rates (Figs. 7a–c). Thus, for a certain subset of the cloud population, aerosol-induced deepening of clouds can counter aerosol suppression of collision–coalescence. On average, however, the aerosol perturbation to microphysics tends to dominate for the 2.5-fold increase in background aerosol conditions between C and H. Stevens and Seifert (2008) demonstrated this compensating effect in their simulations and found that the deepening effect was able to overcome the microphysical effect under very similar drop concentrations but moister, more strongly precipitating conditions; an increase in drop concentration from 35 to 70 mg21 produced an increase in rain rate but further increases in drop concentration caused rain suppression, in agreement with the current simulations (in our simulations, cloudaveraged drop concentrations for C and H are 32 and 75 mg21, respectively). Prior work has indicated that the aerosol tends to increase cloudiness under more strongly precipitating conditions (as in Stevens and Seifert 2008), so the difference between the two studies is perhaps not unexpected. Another exception to convergence in cloud properties is the cloud optical depth t. Calculations of the domainaverage t show that when the LWP is constant (period 3), only about 75%–80% of the t (and to first-order cloud albedo) increase is realized because of horizontal and vertical spatial variability in cloud microphysical properties, and nonlinear dependence of t on drop concentration. In the weakly precipitating cloud simulations, instability changes are mostly a result of aerosol-induced changes in latent heat distribution associated with cloud liquid above cloud base. In contrast, in the more strongly precipitating cases, subcloud evaporation also plays a strong role in removing instability, with the result that the perturbed and unperturbed systems converge more rapidly than in the weakly precipitating cases. Thus, in general, both above-cloud-base changes in latent heat distribution induced by aerosol as well as changes in subcloud layers are important for improved understanding of interactions among aerosol, clouds, and environmental instability. What are the characteristic time scales for equilibration in response to an aerosol perturbation? In the current study, equilibration of fields occurs after about 20 h and is influenced to varying degrees by the timing/ duration of the perturbation and the magnitude of precipitation. For this case, an aerosol perturbation lasting 3 h 50 min causes changes in instability that leaves its imprint on the evolution of the cloud system regardless of the evolution of aerosol perturbation from 10 h on
3630
JOURNAL OF THE ATMOSPHERIC SCIENCES
VOLUME 69
FIG. A1. Time series of (a) LWP and (b) CAPE for C (solid) and H (dashed) simulations as in Fig. 1 (black lines) and for initial conditions with either half (green) or double (red) the initial temperature perturbations. Note the consistency between the three pairs of simulations.
(Fig. 9). Thus, the cloud system becomes insensitive to aerosol differences at some point in the evolution, and it is the earlier aerosol-induced changes in the thermodynamic environment and their feedbacks with cloud microphysics and dynamics that control the cloud field differences between aerosol-perturbed and unperturbed cases [see Lee (2011) for similar results pertaining to deep convective clouds]. A cautionary note: sensitivity to the duration of the aerosol perturbation could well differ if the perturbation were to be applied at a different stage of the cloud field evolution. These results demonstrate that even transient aerosol perturbations may linger longer than the duration of the perturbation itself, provided it persists for long enough. Short-duration aerosol perturbations are unlikely to have much influence on the system but further study is required to place these ideas on firmer footing. Acknowledgments. The authors thank NOAA’s Climate Goal Program and the NOAA/NSF Climate Process Team for supporting this work. P. Chuang thanks the Office of Naval Research (Grant N00014-08-1-0437) for its support of this work. NOAA’s HPCC is acknowledged for its computing support.
APPENDIX A Sensitivity to Initial Conditions In the spirit of ensemble simulations, we perform additional C and H simulations but modify the initial thermodynamic profile with random perturbations that are, respectively, half (Pert/2 C and Pert/2 H runs) and
double (Pertx2 C and Pertx2 H runs) the perturbations for the standard C and H simulations. Figure A1 presents a subset of the fields in Fig. 1, together with those from the modified initial conditions, and shows that the basic trends in C versus H time series are remarkably similar. While this ‘‘ensemble’’ of runs has a very limited number of members, Fig. A1 does suggest that the analysis presented herein is not associated with a very particular set of initial conditions.
APPENDIX B Further Sensitivity Tests Further sensitivity tests are conducted to investigate robustness of the results. The H run is repeated from 7 h (in the middle of the first stage) to the end of the second stage (at 10 h) in the first period to identify the cause of the persistent larger domain-averaged CAPE in the H run between 8 and 10 h (Fig. 1e). The first of the repeated runs, referred to as HuC, adopts u from the C run at each of grid points at 7 h. Then, when the simulation time in the HuC reaches 8 h, u from the C run at each of grid points at 8 h is adopted by HuC. As shown in Figs. 4b and 4d, this u from the C run has less instability below about 1 km at hours 7 and 8. Consequently, the difference in the time- and domain-averaged CAPE between the HuC and C runs averaged over the first stage is much smaller than that between the H and C runs (Table 2). During the second stage, the time- and domain-averaged CAPE is slightly smaller in the HuC run than in the C run (Table 2). To explore the effect of the r difference on the larger domain-averaged CAPE in the H run during the second
DECEMBER 2012
LEE ET AL.
stage, the H run is repeated again from 7 h to the end of the second stage by adopting r from the control run at each of the grid points, first at 7 h and then at 8 h. This repeated run is referred to as HrC. The time- and domainaveraged CAPE in this HrC run is larger than the C run over the second stage (Table 2). Finally, to explore the effect of the aerosol difference on the larger domain-averaged CAPE in the H run during the second stage of period 1, the H run is repeated again from 7 h to the end of the second stage by adopting the aerosol concentration from the C run at each of the grid points, first at 7 h and then at 8 h. This repeated run is referred to as the H-C-second run. The time- and domain-averaged CAPE in this H-C-second run is larger than the C run during the second stage. These repeated simulations demonstrate that the larger instability below about 1 km, established during the first stage of period 1 (between 06 h 10 min and 8 h), is the primary cause of the larger CAPE in the H case during the second stage (8–10 h). They also demonstrate that aerosol and r differences established during the first stage are not the cause of the CAPE differences between the H and C cases during the second stage. It is notable that the differences in cloud-top height and in the location of strong convective divergence around cloud tops between the HuC and the C runs are negligible compared to those between the H and C runs (Fig. 5b). This is because the CAPE difference between the HuC run and the C run is much smaller than that between the H and C runs (Table 2). This results in a similar net cooling in DZm between the HuC run and the C run (Fig. 5b), reinforcing the idea that the larger low-level instability in the H run established during the first stage is the main cause of differences in the net cooling in DZm between the H and C runs. REFERENCES Ackerman, A. S., M. P. Kirkpatrick, D. E. Stevens, and O. B. Toon, 2004: The impact of humidity above stratiform clouds on indirect aerosol climate forcing. Nature, 432, 1014–1017. Albrecht, B. A., 1989: Aerosols, cloud microphysics, and fractional cloudiness. Science, 245, 1227–1230. 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. Brenguier, J.-L., H. Pawlowska, L. Schu¨ ller, R. Preusker, J. Fischer, and Y. Fouquart, 2000: Radiative properties of boundary layer clouds: Droplet effective radius versus number concentration. J. Atmos. Sci., 57, 803–821. Bretherton, C. S., J. Uchida, and P. N. Blossey, 2010: Slow manifolds and multiple equilibria in stratocumulus-capped boundary layers. J. Adv. Model. Earth Syst., 2, doi:10.3894/ JAMES.2010.2.14.
3631
Feingold, G., 2003: Modeling of the first indirect effect: Analysis of measurement requirements. Geophys. Res. Lett., 30, 1997, doi:10.1029/2003GL017967. ——, and H. Siebert, 2009: Cloud–aerosol interactions from the micro to the cloud scale. Clouds in the Perturbed Climate System: Their Relationship to Energy Balance, Atmospheric Dynamics, and Precipitation, J. Heintzenberg and R. J. Charlson, Eds., The MIT Press, 319–338. ——, R. L. Walko, B. Stevens, and W. R. Cotton, 1998: Simulations of marine stratocumulus using a new microphysical parameterization scheme. Atmos. Res., 47, 505–528. Grabowski, W. W., and T. L. Clark, 1993: Cloud–environment interface instability. Part II: Extension to three spatial dimensions. J. Atmos. Sci., 50, 555–573. Hill, A. A., G. Feingold, and H. Jiang, 2009: The influence of entrainment and mixing assumption on aerosol–cloud interactions in marine stratocumulus. J. Atmos. Sci., 66, 1450– 1464. Jiang, H., H. Xue, A. Teller, G. Feingold, and Z. Levin, 2006: Aerosol effects on the lifetime of shallow cumulus. Geophys. Res. Lett., 33, L14806, doi:10.1029/2006GL026024. ——, G. Feingold, and I. Koren, 2009: Effect of aerosol on trade cumulus cloud morphology. J. Geophys. Res., 114, D11209, doi:10.1029/2009JD011750. Kazil, J., H. Wang, G. Feingold, A. D. Clarke, J. R. Snider, and A. R. Bandy, 2011: Chemical and aerosol processes in the transition from closed to open cells during VOCALS-Rex. Atmos. Chem. Phys. Discuss., 11, 4687–4748. Lee, S.-S., 2011: Dependence of aerosol-precipitation interactions on humidity in a multiple-cloud system. Atmos. Chem. Phys., 11, 2179–2196. ——, and G. Feingold, 2010: Precipitating cloud-system response to aerosol perturbations. Geophys. Res. Lett., 37, L23806, doi:10.1029/2010GL045596. Malkus, J. S., 1954: Some results of a trade-cumulus cloud investigation. J. Meteor., 11, 220–237. McFarquhar, G. M., and H. Wang, 2006: Effects of aerosols on trade wind cumuli over the Indian Ocean: Model simulations. Quart. J. Roy. Meteor. Soc., 132, 821–843. Nober, F. J., and H. F. Graf, 2005: A new convective cloud field model based on principles of self-organisation. Atmos. Chem. Phys., 5, 2749–2759. Nuijens, L., B. Stevens, and A. P. Siebesma, 2009: On the environment of shallow precipitating convection. J. Atmos. Sci., 66, 1962–1979. Rauber, R. M., and Coauthors, 2007: Rain in Shallow Cumulus over the Ocean: The RICO campaign. Bull. Amer. Meteor. Soc., 88, 1912–1928. Sengupta, S. K., R. M. Welch, M. S. Navar, T. A. Berendes, and D. W. Chen, 1990: Cumulus cloud field morphology and spatial patterns derived from high spatial resolution Landsat imagery. J. Appl. Meteor., 29, 1245–1267. Slingo, J., and Coauthors, 1994: Mean climate and transience in the tropics of the UGAMP GCM: Sensitivity to convective parametrization. Quart. J. Roy. Meteor. Soc., 120, 881–922. Small, J. D., P. Y. Chuang, G. Feingold, and H. Jiang, 2009: Can aerosol decrease cloud lifetime. Geophys. Res. Lett., 36, L16806, doi:10.1029/2009GL038888. Squires, P., 1958: Penetrative downdraughts in cumuli. Tellus, 10, 381–389. Stevens, B., 2007: On the growth of layers of nonprecipitating cumulus convection. J. Atmos. Sci., 64, 2916–2931.
3632
JOURNAL OF THE ATMOSPHERIC SCIENCES
——, and A. Seifert, 2008: Understanding macrophysical outcomes of microphysical choices in simulations of shallow cumulus convection. J. Meteor. Soc. Japan, 86, 143–162. ——, and G. Feingold, 2009: Untangling aerosol effects on clouds and precipitation in a buffered system. Nature, 461, 607–613, doi:10.1038/nature08281. Tiedtke, M., W. A. Heckley, and J. Slingo, 1988: Tropical forecasting at ECMWF: The influence of physical parametrization on the mean structure of forecasts and analyses. Quart. J. Roy. Meteor. Soc., 114, 639–664. Twomey, S., 1977: The influence of pollution on the shortwave albedo of clouds. J. Atmos. Sci., 34, 1149–1152. vanZanten, M. C., and Coauthors, 2011: Controls on precipitation and cloudiness in simulations of trade-wind cumulus as observed during RICO. J. Adv. Model. Earth Syst., 3, M06001, doi:10.1029/2011MS000056. Wang, H., and G. M. McFarquhar, 2008a: Large-eddy simulations of the diurnal cycle of shallow convection and cloudiness over the tropical Indian Ocean. Quart. J. Roy. Meteor. Soc., 134, 643–661. ——, and ——, 2008b: Modeling aerosol effects on shallow cumulus convection under various meteorological conditions observed over the Indian Ocean and implications for development of
VOLUME 69
mass-flux parameterizations for climate models. J. Geophys. Res., 113, D20201, doi:10.1029/2008JD009914. ——, and G. Feingold, 2009: Modeling mesoscale cellular structures and drizzle in marine stratocumulus. Part I: Impact of drizzle on the formation and evolution of open cells. J. Atmos. Sci., 66, 3237–3256. ——, W. C. Skamarock, and G. Feingold, 2009: Evaluation of scalar advection schemes in the Advanced Research WRF model using large-eddy simulations of aerosol–cloud interactions. Mon. Wea. Rev., 137, 2547–2558. Warner, J., 1968: A reduction in rainfall associated with smoke from sugar-cane fires: An inadvertent weather modification. J. Appl. Meteor., 7, 247–251. Xue, H., and G. Feingold, 2006: Large-eddy simulations of trade wind cumuli: Investigation of aerosol indirect effects. J. Atmos. Sci., 63, 1605–1622. ——, ——, and B. Stevens, 2008: Aerosol effects on clouds, precipitation, and the organization of shallow cumulus convection. J. Atmos. Sci., 65, 392–406. Zhao, M., and P. H. Austin, 2005: Life cycle of numerically simulated shallow cumulus clouds. Part I: Transport. J. Atmos. Sci., 62, 1269–1290.
