PHYSICAL MECHANISMS CONTROLLING SELF-AGGREGATION OF CONVECTION IN IDEALIZED NUMERICAL MODELING SIMULATIONS

Similar documents
PUBLICATIONS. Journal of Advances in Modeling Earth Systems

Radiative-Convective Instability

Radiative-Convective Instability

Radiative-Convective Instability. Kerry Emanuel Massachusetts Institute of Technology

Convective self-aggregation, cold pools, and domain size

Assessing the strength of self-aggregation feedbacks from in situ data

PUBLICATIONS. Journal of Advances in Modeling Earth Systems

A "New" Mechanism for the Diurnal Variation of Convection over the Tropical Western Pacific Ocean

Large-Eddy Simulations of Tropical Convective Systems, the Boundary Layer, and Upper Ocean Coupling

Why do GCMs have trouble with the MJO?

CHAPTER 8 NUMERICAL SIMULATIONS OF THE ITCZ OVER THE INDIAN OCEAN AND INDONESIA DURING A NORMAL YEAR AND DURING AN ENSO YEAR

Diurnal Timescale Feedbacks in the Tropical Cumulus Regime

Effect of Surface Fluxes versus Radiative Cooling on Tropical Deep Convection.

Presentation A simple model of multiple climate regimes

An Energy-Balance Analysis of Deep Convective Self-Aggregation above Uniform SST

Boundary layer equilibrium [2005] over tropical oceans

Large-Eddy Simulations of Tropical Convective Systems, the Boundary Layer, and Upper Ocean Coupling

Climate change with an iris-effect. Thorsten Mauritsen and Bjorn Stevens Max Planck Institute for Meteorology, Hamburg

Large-Eddy Simulations of Tropical Convective Systems, the Boundary Layer, and Upper Ocean Coupling

An energy-balance analysis of deep convective self-aggregation

Analysis of Cloud-Radiation Interactions Using ARM Observations and a Single-Column Model

2.1 Effects of a cumulus ensemble upon the large scale temperature and moisture fields by induced subsidence and detrainment

An Introduction to Physical Parameterization Techniques Used in Atmospheric Models

Moist static energy budget diagnostics for. monsoon research. H. Annamalai

Parameterizing large-scale dynamics using the weak temperature gradient approximation

An Introduction to Climate Modeling

Interactions among Cloud, Water Vapor, Radiation and. Large-scale Circulation in the Tropical Climate. Department of Atmospheric Sciences

ScienceDirect. Numerical modeling of atmospheric water content and probability evaluation. Part I

Hurricanes: Their physics and relationship to climate. Kerry Emanuel Massachusetts Institute of Technology

The Wavelength Dependence of the Gross Moist Stability and the Scale Selection in the Instability of Column-Integrated Moist Static Energy

Tropical cirrus and water vapor: an effective Earth infrared iris feedback?

Hurricanes are intense vortical (rotational) storms that develop over the tropical oceans in regions of very warm surface water.

WaVaCS summerschool Autumn 2009 Cargese, Corsica

Lecture 9: Climate Sensitivity and Feedback Mechanisms

Idealized Cloud-System Resolving Modeling for Tropical Convection Studies. Usama M. Anber

A new theory for moist convection in statistical equilibrium

UNIVERSITY OF READING School of Mathematics, Meteorology and Physics

8.2 Numerical Study of Relationships between Convective Vertical Velocity, Radar Reflectivity Profiles, and Passive Microwave Brightness Temperatures

How surface latent heat flux is related to lower-tropospheric stability in southern subtropical marine stratus and stratocumulus regions

2. Meridional atmospheric structure; heat and water transport. Recall that the most primitive equilibrium climate model can be written

Coupling with Ocean Mixed Layer leads to Intraseasonal variability in tropical deep convection evidence from Cloud-resolving Simulations

Journal of Advances in Modeling Earth Systems

Dynamics and Kinematics

KRISTIN LARSON AND DENNIS L. HARTMANN

The Ocean-Atmosphere System II: Oceanic Heat Budget

Cloud Radiative Feedbacks in GCMs : A Challenge for the Simulation of Tropical Climate Variability and Sensitivity

Convection Trigger: A key to improving GCM MJO simulation? CRM Contribution to DYNAMO and AMIE

Lecture 10: Climate Sensitivity and Feedback

PUBLICATIONS. Journal of Advances in Modeling Earth Systems

Tropical Cyclone Genesis: What we know, and what we don t!

Interhemispheric climate connections: What can the atmosphere do?

October 1986 R. H. Johnson 721. Lower-Tropospheric Warming and Drying in Tropical Mesoscale Convective Systems:

Glaciology HEAT BUDGET AND RADIATION

More on Diabatic Processes

Response of convection to relative SST: Cloud-resolving simulations in 2D and 3D

Climate Feedbacks from ERBE Data

Moisture modes, cloud-radiative feedbacks and the MJO

Geophysics Fluid Dynamics (ESS228)

Single-Column Modeling, General Circulation Model Parameterizations, and Atmospheric Radiation Measurement Data

Parametrizing Cloud Cover in Large-scale Models

Tropospheric Moisture: The Crux of the MJO?

Climate Dynamics (PCC 587): Clouds and Feedbacks

- matter-energy interactions. - global radiation balance. Further Reading: Chapter 04 of the text book. Outline. - shortwave radiation balance

Do climate models over-estimate cloud feedbacks?

Understanding the Greenhouse Effect

Moisture modes, cloud-radiative feedbacks and the MJO

Tropical Convection and the Energy Balance at the Top of the Atmosphere

5. General Circulation Models

Correspondence between short and long timescale systematic errors in CAM4/CAM5 explored by YOTC data

Lecture 14. Marine and cloud-topped boundary layers Marine Boundary Layers (Garratt 6.3) Marine boundary layers typically differ from BLs over land

Chapter 7: Thermodynamics

Weak Temperature Gradient Simulations For Different Convective Environments

A Theory for Buoyancy and Velocity Scales in Deep Moist Convection

Quasi-equilibrium Theory of Small Perturbations to Radiative- Convective Equilibrium States

The Effect of Sea Spray on Tropical Cyclone Intensity

Impact of Explicit Atmosphere Ocean Coupling on MJO-Like Coherent Structures in Idealized Aquaplanet Simulations

Incorporation of 3D Shortwave Radiative Effects within the Weather Research and Forecasting Model

Moisture Modes and the Eastward Propagation of the MJO

Thermodynamics of Atmospheres and Oceans

Cloud-resolving simulation of TOGA-COARE using parameterized large-scale dynamics

A. Parodi 1, (1) CIMA Research Foundation, Italy. in cooperation with: K. A. Emanuel 2, and E. Foufoula-Georgiou 3 (2) EAPS, MIT, USA

MC-KPP: Efficient, flexible and accurate air-sea coupling

PUBLICATIONS. Journal of Advances in Modeling Earth Systems

The Amplification of East Pacific Madden Julian Oscillation Convection and Wind Anomalies during June November

NOTES AND CORRESPONDENCE. On the Radiative and Dynamical Feedbacks over the Equatorial Pacific Cold Tongue

Critical phenomena in atmospheric precipitation Supplementary Information

ZHIMING KUANG. (Manuscript received 16 July 2009, in final form 24 November 2009) ABSTRACT

Northern New England Climate: Past, Present, and Future. Basic Concepts

Numerical simulation of marine stratocumulus clouds Andreas Chlond

Project 3 Convection and Atmospheric Thermodynamics

Spectrum of Radiation. Importance of Radiation Transfer. Radiation Intensity and Wavelength. Lecture 3: Atmospheric Radiative Transfer and Climate

Testing the Fixed Anvil Temperature Hypothesis in a Cloud-Resolving Model

Development of a High-Resolution Coupled Atmosphere-Ocean-Land General Circulation Model for Climate System Studies

Interannual variability of top-ofatmosphere. CERES instruments

ENSO and April SAT in MSA. This link is critical for our regression analysis where ENSO and

Lecture 3: Atmospheric Radiative Transfer and Climate

Ocean-Atmosphere Interactions and El Niño Lisa Goddard

Lecture 7: The Monash Simple Climate

NSF 2005 CPT Report. Jeffrey T. Kiehl & Cecile Hannay

1.4 CONNECTIONS BETWEEN PV INTRUSIONS AND TROPICAL CONVECTION. Beatriz M. Funatsu* and Darryn Waugh The Johns Hopkins University, Baltimore, MD

Transcription:

3B.7 PHYSICAL MECHANISMS CONTROLLING SELF-AGGREGATION OF CONVECTION IN IDEALIZED NUMERICAL MODELING SIMULATIONS Allison A. Wing* and Kerry A. Emanuel Program in Atmospheres, Oceans, and Climate, Massachusetts Institute of Technology, Cambridge, MA. INTRODUCTION Moist convection in the tropical atmosphere is often organized into clusters containing many individual convective cells. The value of understanding how convection organizes cannot be overstated. Clusters of organized convection are ubiquitous in the tropics (e.g., Machado and Rossow, 993; Mapes and Houze Jr., 993) and have important impacts on weather and climate. They are responsible for much of the rainfall and cloudiness over the tropics, with approximately 5% of tropical rainfall due to mesoscale convective systems (Nesbitt et al., ). This allows tropical cloud clusters to modulate the radiative heating of the surface and atmosphere and influence the large-scale circulation and moisture distribution. In idealized modeling studies (e.g., Bretherton et al., 5) and observations (Tobin et al., 22, 23), there is a systematic dependence of mean humidity and radiative fluxes on the degree of convective aggregation. Therefore, understanding how and why tropical convection organizes is important for understanding both tropical and global climate variability, and climate sensitivity. Globally, 6.4% of tropical cloud clusters evolve into tropical cyclones each year (Hennon et al., 2), so tropical convective organization is also tightly linked to the problem of tropical cyclogenesis. Finally, understanding the mechanisms by which convection organizes may lead to insights on the Madden-Julian Oscillation, which can be considered convective organization on a large scale. The MJO has a direct impact on weather in the Indian and western Pacific Oceans and modulates tropical cyclone activity, yet, a complete theory for its existence and propagation characteristics remains elusive. In this study, the fundamental mechanism underlying the self-aggregation of convection is explored using a cloud-resolving model. The objective is to identify and quantify the interactions between the environment and the convection that allow the convection to spontaneously organize into a single cluster. This study addresses several questions: *Corresponding author address: Allison A. Wing, Massachusetts Institute of Technology 54-85, Cambridge, MA 239; email: awing@mit.edu How does self-aggregation evolve? Which feedbacks are important and what are their magnitudes? What physical mechanisms are behind the feedbacks controlling self-aggregation? 2. MODEL SIMULATIONS The model used is the System for Atmospheric Modeling, henceforth referred to as SAM (Khairoutdinov and Randall, 3). SAM is a three-dimensional cloud resolving model that employs the anelastic equations of motion. The prognostic thermodynamics variables are the total nonprecipitating water, total precipitating water, and the liquid water/ice static energy. The simulations discussed here were performed with a domain size of 768 x 768 km 2 with 64 vertical levels and rigid lid at 28 km. The horizontal resolution was 3 km, and doubly periodic lateral boundary conditions were employed. The model was initialized with a sounding from the domain average of a smaller domain run in radiative-convective equilibrium with white noise added in the boundary layer temperature field. There is no mean wind or other external forcing imposed. We used a fully interactive radiation scheme (RRTM), with solar insolation constant and equal to a value of 43.98 W/m 2 and no diurnal cycle. The surface sensible and latent heat fluxes were computed interactively. Finally, we performed simulations at fixed sea surface temperature (SST), with values between 295K and 32K, but most of the results shown here are for the simulation at 35 K. 3. EVOLUTION OF SELF-AGGREGATION Simulations of convection in radiative-convective equilibrium using three-dimensional cloud system resolving models often produce distributions of convection that are nearly random in space and in time. Figure a shows an example of this with a snapshot of

the outgoing longwave radiation, indicating the existence of high clouds. a) 2 Day OLR (W/m ) 32 28 26 Y (km) 24 22 8 6 4 b) X (km) 2 Day 8 OLR (W/m ) 32 28 26 Y (km) 24 22 We perform a sensitivity test at 35 K in which we prescribe the water vapor profile used in the longwave radiation calculation (but still use the modeled clouds) and find that this prevents aggregation. Prescribing the water vapor profile used in the shortwave radiation calculation does not prevent aggregation, indicating it is the longwave radiation - water vapor interaction that is key (in contrast to the results of Muller and Held (22), who found that the longwave-low cloud interaction was vital). A second sensitivity test in which we use the horizontal mean surface wind speed in the surface enthalpy flux calculation reveals that, at least for the simulation considered here, wind-dependent surface fluxes are necessary for aggregation. 8 6 moist regions increases. Therefore, when examining the initiation of self-aggregation, we focus on processes that can amplify the initial dry patch. 4. ANALYSIS FRAMEWORK 4 X (km) Figure : Snapshot of outgoing longwave radiation (OLR) at day (a) and day 8 (b) of a radiative-convective equilibrium simulation at 35 K. However, when certain conditions are met, the convection becomes organized into a single, intensely convecting moist clump surrounded by a broad region of dry subsiding air, a process termed selfaggregation (e.g., Bretherton et al., 5). In Figure b, a snapshot of the outgoing longwave radiation from a later time in the same simulation as Figure a, all the clouds are confined to a single cluster. Previous work indicated that cloud-water vaporradiation feedbacks that dry the drier air columns and moisten the moister air columns are essential to the self-aggregation process (Tompkins and Craig, 998; Bretherton et al., 5; Stephens et al., 8; Muller and Held, 22). Simulations at SST s of 3K, 33K, 35K and 37K self-aggregated, while the simulations at lower values of SST did not. We note in passing that the simulations at higher SSTs (i.e., 3 K) need a larger domain to aggregate. Aggregation is characterized by a dramatic increase in the domain averaged outgoing longwave radiation, which coincides with a dramatic decrease in the tropospheric humidity. An examination of the evolution of self-aggregation reveals that it begins as a dry patch that expands, eventually forcing all the convection into a single clump (Wing and Emanuel, 23). Over the first half of the simulation, the dry regions get progressively drier and larger; it is not until day 5 that the humidity in the We frame our analysis in terms of the budget for the spatial variance of vertically integrated frozen moist static energy (FMSE). Since column radiative flux conference and surface enthalpy fluxes are diabatic sources and sinks of vertically integrated FMSE, framing our analysis in this matter allows us to quantify these potential feedbacks. The frozen moist static energy, which is conserved during moist adiabatic processes in the model, including the freezing and melting of precipitation, is given by h = cp T + gz + Lv qv Lf qice () where Lf is the latent heat of fusion and qice represents all ice phase condensates. The budget equation for the spatial variance of vertically integrated frozen moist static energy is given by h 2 = h SEF + h NetSW 2 t +h NetLW h h ~uc h (2) where the primes indicate an anomaly from the horizontal mean. h indicates the density weighted vertical integral of FMSE. NetLW and Net SW are the column longwave and shortwave radiative flux convergences, respectively, given by NetLW = LWsfc LWtop, (3) NetSW = SWtop SWsfc. (4) and The surface enthalpy flux anomalies can be written

as SEF = (i) ρl v (c E U) { q} + ρc p (c H U) { T } (ii) + ρl v {c E U} q + ρc p {c H U} T (iii) + ρl v (c E U) q + ρc p (c H U) T (iii) ρl v {(c E U) q } ρc p {(c H U) T }, (5) in which {A} represents a horizontal mean, A represents the anomaly from that mean, q = q T s q v, and T = T s T a. Term (i) in Equation 5 represents the part of the SEF anomaly due solely to variations in the surface wind speed. Term (ii) represents the part of the SEF anomaly due solely to variations in the air-sea enthalpy disequilibrium. Term (iii) represents the part of the SEF anomaly due to the product of variations in the surface wind speed and variations in the air-sea enthalpy disequilibrium, which we refer to as the eddy term. The surface exchange coefficients, c E and c H, vary strongly with the surface wind speed but only weakly with the air-sea disequilibrium over the range of values found in the simulations; therefore, we have combined the exchange coefficients with the surface wind speed when calculating horizontal means and anomalies with those means. Self-aggregation is associated with an increase in the variance of vertically integrated FMSE; therefore, processes that increase ĥ 2 favor self-aggregation. It is then clear from Equation 2 that if the correlation between the anomaly of a diabatic term and ĥ is positive, then there is an anomalous source of FMSE in a region of already high FMSE and thus a positive feedback on aggregation. For each of the terms in Equation 2, we take a daily average, and then horizontally average over 48 x 48 km 2 blocks to focus on the mesoscale organization. We then sort the blocks according to their column relative humidity, allowing us to examine how the terms evolve as a function of time and moisture-rank. Finally, we normalize each term by the horizontal mean of FMSE variance, {ĥ 2 }. Because ĥ 2 increases with time, normalizing in this manner makes it easier to interpret what is happening in the early stage of aggregation, when both the vertically integrated FMSE anomalies and forcing terms are small. 5. RESULTS Figure 2a is a hovmuller plot of the sum of the diabatic feedback terms (surface flux and column radiative flux convergence terms in Equation 2). The sum of the diabatic feedbacks is positive during the first twenty days of the simulation and strongest in the driest regions. With time, the positive values propagate toward moister regions, in association with the expanding dry patch. This evolution is consistent with our observation that self-aggregation begins as a dry patch that expands. In the moist regions, positive values persist throughout the simulation. We note that the diabatic terms are a negative feedback in the dry regions from day 3 to 7, but ĥ 2 is still increasing in those regions over that time period (not shown). This indicates that the horizontal convergence of the vertically integrated flux of FMSE, calculated as a residual from the rest of the budget, must be playing a role. We find that, in the intermediate stages of aggregation, the convergence of the flux of FMSE by the circulation is important in amplifying vertically integrated FMSE anomalies and promoting self-aggregation (Figure 2b). This is consistent with Bretherton et al. (5) and Muller and Held (22), who found that mesoscale circulations intensify the later stages of self-aggregation via an up gradient transfer of moist static energy, and will not be explored further here. Figure 3a-b reinforces the notion of competition between positive and negative feedbacks, as we see that the correlations of vertically integrated FMSE anomalies with the column radiative flux convergence anomalies are mostly positive during the first 6 days of the simulation (when the cluster is developing), while the correlations with the surface enthalpy flux anomalies are predominantly negative from day 2 to day 6. The total surface flux feedback (Figure 3b) is positive during the first twenty days of the simulation and is largest in the driest regions. In Figures 3c-f, the feedback terms are further decomposed. First, we examine the surface flux feedback and note that the correlation between vertically integrated FSME anomalies and surface enthalpy flux anomalies due to wind speed anomalies (Figure 3d) is mostly positive. This is the "WISHE" feedback; convective gustiness in the moist, intensely convective regions enhances the surface fluxes there. However, while WISHE is a positive feedback for aggregation in our simulations, it s primary role is to counteract a strongly negative surface flux feedback due to variations in the air-sea enthalpy disequilibrium. Because the simulations have a fixed, uniform sea surface temperature and the surface en-

ˆ -h h S E F + h N etlw + h N ets W 8 8 7 7 6 5 4 3.5.5 uˆh 9 9.5 6 5 4 3.5 2 2 h 5 5 25 (a) Sum of diabatic terms 5 5 25 c (b) Convergence term: h h ~ u h c Figure 2: Sum of all diabatic correlation terms (a) and convergence term (b), h h ~ u h, (calculated as a residual) in ver2 tically integrated FMSE spatial variance budget, normalized at each time by {h }. Plotted as a function of time (y-axis) and moisture space (x-axis), where each term has been averaged over a day and a 48 48 km2 block, has units of days, and has been sorted according to block-averaged column relative humidity (CRH). On the x-axis, dry regions are on the left and moist regions are on the right. Results are from the simulation at 35 K. The black line is the h = contour. Note that the color bar saturates. thalpy flux is dominated by the latent heat component, the air-sea enthalpy disequilibrium essentially follows the boundary layer water vapor mixing ratio, which is diminished in the dry columns, enhancing the surface fluxes there (a negative feedback on self-aggregation). Finally, while not shown here, the "eddy term" involving the correlation between vertically integrated FMSE anomalies and the product of wind speed and air-sea disequilibrium anomalies reflects that the wind speed and disequilibrium anomalies are anti-correlated. Summing these components yields a total surface flux feedback that transitions from positive to negative (Figure 3b). Next, we examine the correlation between the vertically integrated FMSE anomalies and the column shortwave flux convergence anomalies (Figure 3e). This term is positive nearly everywhere, reflecting negative anomalies in NetSW in the dry regions (where h < ) and positive anomalies in NetSW in the moist regions (where h > ). This is due to increasing atmospheric absorption of shortwave radiation by water vapor as we move from dry to moist regions, although this is also modulated by cloud effects (not shown here). In the regions where there are low clouds, the clouds act to increase the column shortwave flux convergence, whereas high clouds in the very moistest regions block solar radiation from penetrating into the lower parts of the atmosphere, decreasing the column shortwave flux convergence. The feedback term involving column longwave flux convergence anomalies (Figure 3c) is more complicated, as it is can be either positive or negative. Ig- noring clouds for the moment, the column longwave flux convergence varies between dry and moist regions because variations in atmospheric water vapor determines variations in the longwave emissivity. Decreasing the water vapor causes competing responses of the longwave radiative fluxes. One effect of decreasing the water vapor in the upper troposphere (as occurs early in the simulation) is that the concentration of longwave emitters will decrease, which locally causes it to radiatively cool less (tending to decrease a sink of energy for the column). However, lower parts of the atmosphere will cool to space more effectively through a more transparent upper troposphere (tending to increase a sink of energy). In the first 2 days, the second effect "wins" and the longwave feedback is initially positive, amplifying the developing dry patch. Later, as the dry perturbation amplifies and the lower troposphere also becomes drier, the longwave feedback transitions from positive to negative because the emissivity has decreased to the extent that the dry regions are no longer able to cool to space effectively (decreasing a sink of energy in the low FMSE regions). A calculation of the longwave feedback term, h NetLW, using the clear sky longwave fluxes reveals that a large part of the positive longwave feedback at the beginning of the circulation (Figure 3c) is captured by the clear sky processes (Figure 4). Conversely, once the cluster is established (day 6 onward), the longwave feedback is strongly positive in the moistest regions primarily because the column longwave cooling is strongly reduced by the longwave opacity and low temperature

a) h N etlw + h N ets W 9.5 6 5 4 3 5 5 5 3 5 4 3 5 5 7 3.5 2 5 6 5 4 3.5 2 25 5.5 5 4 8.5 6 25 h ρl v { ce U } q + h ρcp { ch U } T f) 9 5 7 25 h N etsw e) 25 9 8 7.5 6 8 5 5 9.5 2.5 6 5 h ρl v (ce U ) { q} + h ρcp (ch U ) { T } 7 d) 8.5 2 25 9 4 h N etlw c).5 6 7.5 2 8 7 9 8 h S E F b) 5 4 3.5 2 5 5 25 Figure 3: Left column: correlation between vertically integrated FMSE anomalies and column radiative flux convergence anomalies (a: column radiative flux convergence, c: column longwave convergence, e: column shortwave convergence). Right column: correlation between vertically integrated FMSE anomalies and surface enthalpy flux anomalies (b: total surface enthalpy flux anomaly, d: anomaly due to surface wind speed anomalies, f: anomaly due to air-sea enthalpy disequilibrium anomalies). All terms have been averaged over each day and over 48 48 km2 blocks, normalized by {h 2 }, are from the simulation at 35 K, and have units of days. On the x-axis, dry regions are on the left and moist regions are on the right, sorted according to block-averaged column relative humidity (CRH). The black line is the h = contour, plotted as a reference. Note that the color bar saturates in a few places.

of high clouds (compare Figure 4 and Figure 3c). ĥ NetLW Clear Sky 9 8 7 6 5 4 3 2.5.5 by NCAR s Computational and Information Systems Laboratory, sponsored by the National Science Foundation. This work was supported by NSF grants 32244, 3648, 85639. The first author received additional support from the MIT Joint Program on the Science and Policy of Global Change. Much of this work was published as part of Wing and Emanuel (23) and is adapted here in accordance with the Creative Commons Attribution License. 5 5 25 9. REFERENCES Figure 4: Same as Figure 3c, but for clear sky longwave radiative fluxes. 7. SUMMARY AND IMPLICATIONS FOR TEMPER- ATURE DEPENDENCE We introduced a new analysis framework based on the budget for the spatial variance of vertically integrated FMSE, which allowed us to meaningfully quantify each feedback throughout the evolution of self-aggregation. Each of the feedbacks discussed have comparable magnitudes, indicating that each physical mechanism is important for some stage of self-aggregation. Notably, the mechanisms that amplify the initial dry patch and control the evolution to an aggregated state are different from those that maintain the cluster once it is established. The fact that the longwave - water vapor feedback can be positive or negative depending on the balance of opposing responses to a moisture perturbation suggests it as a candidate for explaining the temperature dependence of aggregation. The results of the sensitivity test in which a prescribed water vapor profile is used in the radiation calculation also point to longwave radiation, rather than shortwave. A companion study, Emanuel et al. (23), included the longwave - water vapor feedback discussed here in a simple two-layer model (which also included a representation of convection and large-scale vertical velocity), and showed that the model becomes unstable at high temperatures. The instability results because, at high temperature, the emissivity of the lower troposphere is large enough (due to high water vapor concentration) that variations of its radiative cooling depend primarily on variations in upper tropospheric water vapor. 8. ACKNOWLEDGEMENTS We thank Marat Khairoutdinov for providing SAM, the cloud resolving model. We also acknowledge high-performance computing support from Bluefire and Yellowstone (ark:/8565/d7wd3xhc) provided Bretherton, C. S., P. N. Blossey, and M. Khairoutdinov, 5: An energy-balance analysis of deep convective self-aggregation above uniform SST. J. Atmos. Sci., 62, 4237 4292. Emanuel, K., A. A. Wing, and E. M. Vincent, 23: Radiative-convective instability. J. Adv. Model. Earth Syst., 5, doi:.2/23ms 27. Hennon, C. C., C. N. Helms, K. R. Knapp, and A. R. Bowen, 2: An objective algorithm for detecting and tracking tropical cloud clusters: Implications for tropical cycloge prediction. J. Atmos. Oceanic Technol., 28, 7 8. Khairoutdinov, M. F. and D. A. Randall, 3: Cloud resolving modeling of the ARM summer 997 IOP: Model formulation, results, uncertainties, and sensitivities. J. Atmos. Sci., 6, 67 625. Machado, L. A. T. and W. B. Rossow, 993: Structural characteristics and radiative properties of tropical cloud clusters. Monthly Weather Review, 2, 3234 326. Mapes, B. E. and R. A. Houze Jr., 993: Cloud clusters and superclusters over the oceanic warm pool. Monthly Weather Review, 2, 398 45. Muller, C. J. and I. M. Held, 22: Detailed investigation of the self-aggregation of convection in cloud resovling simulations. J. Atmos. Sci., 69, 255 2565. Nesbitt, S. W., E. J. Zipser, and D. J. Cecil, : A census of precipitation features in the tropics using TRMM: Radar, ice scattering, and lightning observations. J. Climate, 3, 487 46. Stephens, G. L., S. van den Heever, and L. Pakula, 8: Radiative-convective feedbacks in idealized states of radiativeconvective equilibrium. J. Atmos. Sci., 65, 3899 396. Tobin, I., S. Bony, C. E. Holloway, J. Y. Grandpeix, G. Seze, D. Coppin, S. J. Woolnough, and R. Roca, 23: Does convective aggregation need to be represented in cumulus parameterizations? J. Adv. Model. Earth Syst., 5, doi:/.2/jame.2 47. Tobin, I., S. Bony, and R. Roca, 22: Observational evidence for relationships between the degree of aggregation of deep convection, water vapor, surface fluxes, and radiation. J. Climate, 25, 6885 694. Tompkins, A. and G. Craig, 998: Radiative-convective equilibrium in a three-dimensional cloud-ensemble model. Q. J. R. Meteorol. Soc., 24, 273 297. Wing, A. A. and K. A. Emanuel, 23: Physical mechanisms controlling self-aggregation of convection in idealized numerical modeling simulations. J. Adv. Model. Earth Syst., 5, doi:.2/23ms 269.