Clouds and in the Tropical Transition Layer over Mesoscale Convective Systems

KATRINA S. VIRTS AND ROBERT A. HOUZE JR. Department of Atmospheric Sciences, University of Washington, Seattle, Washington

(Manuscript received 1 May 2015, in final form 5 August 2015)


Observations from A-Train satellites and other datasets show that mesoscale convective systems (MCSs) affect the and content of the tropical tropopause transition layer (TTL). The largest MCSs with reflectivity characteristics consistent with the presence of large stratiform and anvil regions have the greatest impact. Most MCSs are associated with in the TTL. Composites in MCS-relative co- ordinates indicate enhanced cloudiness and ice water content (IWC) extending toward the cold-point tro- popause (CPT), particularly in large and connected MCSs. Widespread anvils in the lower TTL are evident in the peak cloudiness diverging outward at those levels. Upper-tropospheric water vapor are enhanced near MCSs. Close to the centers of MCSs, water vapor is suppressed at TTL base, likely because of the combined effects of reduced moistening or dehydration at the higher TTL relative and sub- sidence above top. Weak moistening is observed near the CPT, consistent with sublimation of at the tops of the deepest MCSs. In the outflow region, moistening is observed in the lower TTL near the largest MCSs. Enhanced water vapor in the upper and lower TTL extends beyond the area of substantially enhanced cloudiness and IWC, in agreement with the observed radial outflow, indicating that MCSs are injecting water vapor into the environment and consistent with the possibility that MCS devel- opment may be favored by a premoistened environment.

1. Introduction is a positive maximum in the upper troposphere, with the strongest heating near anvil cloud top (Houze 1982). Deep convective clouds in the tropics sometimes ag- Observations from the CloudSat satellite, launched in gregate and develop extensive regions of stratiform pre- 2006, have permitted the analysis of the vertical structure cipitation and nonprecipitating anvil. Large, organized of clouds around the world. Cetrone and Houze (2009) cloud features of this type are called mesoscale convective and Yuan et al. (2011, hereafter YHH11) have used systems (MCSs; Houze 2014, chapter 9). Over half of CloudSat to examine characteristics of tropical MCS an- tropical and latent heating is produced by vils. The distribution of CloudSat reflectivities varies with MCSs (Yuan and Houze 2010, hereafter YH10). The la- anvil thickness and distance from the MCS center and is tent heating in MCSs is top heavy as a result of their sig- consistent with the detrainment of ice particles from the nificant stratiform precipitation areas (Houze 1982, 1989), MCS updrafts into the anvils and the sedimentation, or and this top-heavy heating is important in forcing the gravitational settling and eventual removal, of the larger large-scale tropical circulation (Hartmann et al. 1984; particles as the anvil age increases. Both studies noted Schumacher et al. 2004). The net radiative heating in MCSs contrasts between the anvil structures of and oceanic MCSs. Land MCSs exhibit broader and higher reflectivity distributions consistent with more intense and less extensive stratiform areas, while oceanic MCSs Denotes Open Access content. exhibit narrower distributions with modal intensity strongly decreasing with height, consistent with the de- velopment of broad stratiform areas. Corresponding author address: Katrina Virts, Department of Atmospheric Sciences, University of Washington, 408 ATG Bldg., MCS vertical structure in the troposphere has impli- Box 351640, Seattle, WA 98195-1640. cations for the systems’ impacts on the overlying tropical E-mail: [email protected] tropopause transition layer (TTL). The TTL extends

DOI: 10.1175/JAS-D-15-0122.1

Ó 2015 American Meteorological Society 4740 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 72 from around 14 km, just above the level of main convective 2006). These observations suggest that tropical convec- outflow, to the highly stable above 18 km tion affects the moisture budget of the TTL. (Fueglistaler et al. 2009; Randel and Jensen 2013). Air be- Analysis of the impact of deep convection on water low the TTL undergoes net radiative cooling as it subsides vapor profiles in the TTL has addressed the fundamental toward the ground, while air in the TTL tends to ascend question of whether the net contribution is to moisten or slowly while undergoing net radiative heating. The TTL dehydrate the environmental air. Relative with links the rising branch of the tropospheric Hadley circula- respect to ice increases with height in the upper tropo- tion, whose ascent is concentrated in areas of tropical moist sphere and TTL because of the decrease in saturation convection, and the broader rising branch of the strato- mixing ratio in the ascending cold air, and supersatura- spheric Brewer–Dobson Circulation. The TTL is charac- tion frequently occurs in the TTL (Luo et al. 2007; terized by extremely low (Anthes et al. 2008; Fueglistaler et al. 2009). Previous studies have empha- Kim and Son 2012), and air ascending through the TTL sized the importance of the ambient relative humidity undergoes freeze drying, producing the low stratospheric in determining the impact of convection: convection in water vapor concentrations noted by Brewer (1949). subsaturated regions moistens the environment, while in The mean level of neutral (LNB) in con- supersaturated regions net dehydration is observed, be- vective regions of the tropics lies between 12 and 14 km, cause excess water vapor is deposited onto the ice parti- just below TTL base (Takahashi and Luo 2014). Di- cles, which may be removed via sedimentation (Jensen vergent outflow from regions of tropical deep convection et al. 2007; Wright et al. 2009; Hassim and Lane 2010). occurs in the upper troposphere and at TTL base Chae et al. (2011), analyzing water vapor anomalies from (Highwood and Hoskins 1998; Folkins and Martin 2005). the Limb Sounder (MLS) below the cloud- Based on CloudSat observations of anvil characteristics, top heights identified in data from the Cloud– Takahashi and Luo (2012) concluded that maximum Lidar and Pathfinder Satellite Observations mass detrainment from convection is around 11 km, just (CALIPSO) satellite, observed strong moistening in the below the LNB determined from atmospheric soundings. upper troposphere and decreased moistening above The geographic distribution of cirrus in the upper tro- TTL base, shifting to dehydration above 16 km if the posphere strongly resembles the distribution of tropical tops of clouds were above this level, consistent with less precipitation, consistent with cirrus forming in and effective moistening with increasing background hu- spreading out from the upper portions of deep convec- midity. They furthermore emphasized that the mecha- tion (Virts et al. 2010). Convectively generated cirrus nism they posited is not restricted to convective towers also occur in the TTL (Massie et al. 2002; Luo and and cirrus outflow but could apply to any cirriform Rossow 2004), although in situ formation in ascending cloud; the mechanism could occur in situ or in the upper layers becomes increasingly important at higher alti- reaches of stratiform or anvil tops of MCSs. tudes (Riihimaki and McFarlane 2010; Virts et al. 2010). As mentioned above, a fraction of deep convective Trajectory analysis suggests that cirrus can persist up to clouds penetrate the CPT. The penetrating towers con- 2 days after convection weakens and be advected up to tain extremely cold, dry air and have been hypothesized 1000 km from the region of the parent convection (Luo as a mechanism for lower-stratospheric dehydration and Rossow 2004). As this occurs, sublimation (Danielsen 1993; Sherwood and Dessler 2000). How- of the smaller ice crystals adds vapor to the environment ever, modeling studies of such clouds indicate net hy- (Soden 2004). dration because they leave behind small ice particles Observations across the tropics indicate that some that undergo sublimation in the relatively warmer deep convection overshoots its LNB and penetrates into stratosphere (Chaboureau et al. 2007; Jensen et al. the TTL. Overshooting occurs over both and land, 2007). These modeling results have been confirmed by although the most vigorous convection is concentrated observational evidence (e.g., Corti et al. 2008; Khaykin over land (Alcala and Dessler 2002; Liu and Zipser 2005; et al. 2009). However, the horizontal scale of the pene- Zipser et al. 2006). Alcala and Dessler (2002) found that trating cores of active convection is very small, and it approximately 5% of convective elements observed by seems unlikely that they alone can account for the water the Tropical Rainfall Measuring Mission (TRMM) ra- vapor properties of the TTL. dar penetrated into the TTL. Takahashi and Luo (2012) Rossow and Pearl (2007) concluded that most con- reported that mean cloud-top height of the deepest vection penetrating near the CPT is associated with or- convective cores observed by CloudSat is 15 km, within ganized (i.e., large) MCSs; however, to date little the TTL. However, less than 1% of deep convective analysis has been attempted of MCS characteristics in clouds reach the cold-point tropopause (CPT; Alcala the TTL or how MCSs might affect the moisture content and Dessler 2002; Gettelman et al. 2002; Dessler et al. of the TTL. Yet MCSs dominate much of the upper-level DECEMBER 2015 V I R T S A N D HOUZE 4741

TABLE 1. Methodology used by YH10 to identify cloud and precipitation features and MCSs, based on MODIS Tb11 and AMSR-E AE_ Rain data. The cloud-top minimum Tb11RC1min is defined as the mean Tb11 of the coldest decile of the largest raining core (RC1). Adapted from YHH11.

Abbreviation Name Definition

HCC High cloud complex Region of MODIS Tb11 contained within a single 260-K isotherm HCS High cloud system Portion of HCC associated with a particular minimum value of Tb11 2 PF Precipitation feature Region of AMSR-E AE_Rain parameter surrounded by 1 mm h 1 contour RC Raining core Portion of any PF overlapping and/or located within an HCS 2 HRA Heavy rain area Portion of PF with rain rate greater than 6 mm h 1 MCS Mesoscale convective Any HCS whose largest RC satisfies the following criteria: system 1) Exceeds 2000 km2 in total area 2 2) Accounts for more than 70% of the total area with rain rate greater than 1 mm h 1 inside the HCS

3) Minimum cloud-top temperature above the RC (indicated by Tb11RC1min) is less than 220 K 4) More than 10% of RC is occupied by HRAs SMCS Separated MCS The largest RC of the MCS is part of a PF that contains less than three dominant RCs of any MCS CMCS Connected MCS The largest RC of the MCS is part of a PF that contains dominant RCs of at least three MCSs cloud coverage of the tropics. The stratiform portions of distribution of clouds, ice, and water vapor in MCSs in MCSs are orders of magnitude larger in scale than the the TTL is presented in section 5. Conclusions are in convective towers embedded in the MCSs, and the strat- section 6. iform regions are buoyant and have very cold tops rising and diverging over wide areas, suggesting that they might have a substantial effect on the moisture of the TTL. In 2. Data this paper, we investigate the water vapor and ice content a. Observations from A-Train satellites of MCSs in the upper troposphere and TTL to determine the impact of these larger convective entities on the water NASA’s A-Train constellation operates in a - content at upper levels. synchronous orbit, crossing the equator around 1:30 To identify MCSs, we use the dataset of YH10, who and 13:30 LT. The A-Train instruments used in this introduced a methodology for identifying MCSs in ob- study, including the YH10 MCS database derived from servations from NASA’s A-Train satellites (L’Ecuyer observations from two of those instruments, are de- and Jiang 2010): high, cold cloud tops observed by the scribed in this section. Moderate Resolution Imaging Spectroradiometer 1) MCS DATABASE (MODIS) and precipitating areas observed by the Ad- vanced Microwave Scanning Radiometer for Clouds and precipitating areas are two prominent Observing System (AMSR-E; see section 2 for details). MCS components readily observed by satellites. The Other instruments aboard A-Train satellites observe YH10 methodology begins with the identification of cloud vertical structure (CloudSat), including that of these components in data from two instruments on the optically thin cirrus layers (CALIPSO), and measure Aqua satellite: MODIS 10.8-mm brightness tempera- profiles of ice water content (IWC) and water vapor tures (TB11) and the AMSR-E AE-Rain product in the TTL (MLS). Collectively, these (Kummerow et al. 2001; Wilheit et al. 2003; Kummerow observations offer a unique to analyze the and Ferraro 2006). Their methodology is summarized in vertical structure of MCSs and their contribution to Table 1. clouds and moisture in the TTL. We analyze these ob- YH10 then divide MCSs into two categories. Sepa- servations of water vapor and ice in the vicinity of each rated MCSs (SMCSs) are those in which the pre- MCS in the YH10 database. cipitation feature (PF) containing the largest raining The remainder of the paper is organized as follows: the core (RC) contains at most one other dominant RC of YH10 MCS identification scheme is described in section any MCS. Long-lived MCSs in close vicinity sometimes 2, along with the other satellite datasets used in this study. merge to form larger precipitating systems (Williams The climatological distributions of MCSs and of clouds and Houze 1987; Mapes and Houze 1993). This type of and moisture in the TTL are shown in section 3.MCS system is designated by YH10 as a connected MCS vertical structure and microphysical characteristics in (CMCS), in which a rain area contains the largest RC of the troposphere are discussed in section 4, while the at least three MCSs. SMCSs occur more frequently and 4742 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 72 exhibit a greater range of sizes than CMCSs, so they are 4) MLS ICE WATER CONTENT AND WATER VAPOR further subdivided into small SMCSs (the smallest 25% The MLS aboard the satellite observes by area, with HCSs less than approximately 11 000 km2) microwave emissions in five spectral regions. We ana- and large SMCSs (the largest 25%, with HCSs greater lyze MLS version 3 vertical profiles of IWC and water than approximately 41 000 km2). In this study, we ana- vapor mixing ratio, whose validations are described by lyze all MCSs identified in the tropics (308N–308S) Read et al. (2007) and Wu et al. (2008), respectively. during 2007–10. Profiles are provided at 1.58 intervals along track and 2) CLOUDSAT REFLECTIVITY have vertical resolutions of 4 and 3 km, respectively, and horizontal resolutions of approximately 7 km CloudSat carries a 94-GHz cloud profiling radar across track and approximately 300 km along track. (CPR) that measures backscatter from clouds (Stephens The useful altitude ranges are from 215 to 83hPa for et al. 2002; Marchand et al. 2008). The backscatter IWC and above 316 hPa for water vapor. The un- profiles are available every 2.5 km along track and certainties for water vapor are 20% below 100 hPa and have a vertical resolution of 240 m at nadir. A cloud 10% at 100 hPa and above, while the uncertainty for mask is provided in the CloudSat geometric profile IWC is a factor of 2, largely owing to the uncertainty product (2B-GEOPROF). In this study, observations from particle size assumption (Jiang et al. 2012). The within clouds are identified by screening for cloud mask data are screened following the recommendations in values greater than or equal to 20. Livesey et al. (2011). The CloudSat profiles analyzed in this study are those MLS also observes , with identified by YHH11 as sampling some portion of an recorded biases of about 1.5 K in the TTL (Schwartz MCS. In that study and in YH10, it was noted that et al. 2008). For this paper, we apply an offset at each clouds with tops above 10 km identified by CloudSat level based on comparison with global posi- fall into two categories: precipitating clouds with bases tioning system (GPS) temperature observations, as near the surface and elevated anvils with bases pri- given by Chae et al. (2011). The water vapor and cor- marily above the freezing level. Accordingly, they de- rected temperature observations are then used to cal- fined anvil clouds as having base above 3 km and top culate relative humidity with respect to ice (referred to above 10 km. The same definition is used in this study. herein as relative humidity). Because it is dependent We will also present statistics on CloudSat profiles upon the resolution of the temperature observations, the sampling nonanvil regions of the MCS. This category is vertical resolution of the relative humidity is coarser dominated by profiles sampling the deep convective or than that of water vapor: approximately 5 km in the TTL stratiform precipitating areas of the system and is re- (Schwartz et al. 2008; Livesey et al. 2011). ferred to herein as the precipitating category.

3) CALIPSO CLOUD OCCURRENCE b. ERA-Interim y The Cloud–Aerosol Lidar with Orthogonal Polar- The horizontal (u, ) and vertical velocities in v ization (CALIOP) is a two-wavelength polarization pressure coordinates ( ) in this study are from the lidar carried aboard CALIPSO. Cloud layer base and European Centre for Medium-Range Fore- top heights are identified from the lidar backscatter casts (ECMWF) interim reanalysis (ERA-Interim; and are provided at a resolution of 60 m in the vertical Dee et al. 2011). ERA-Interim fields are available four 8 in the TTL and 5 km along track. Following Fu et al. times daily at 1.5 resolution for five pressure levels in (2007), opaque features (those that completely at- the upper troposphere and TTL (200, 175, 150, 125, and y tenuate the lidar signal) are assumed to be deep 100 hPa). The u and fields are used to calculate the convective clouds and assigned a cloud base at radial —that is, the component of the horizontal Earth’s surface. While CloudSat’s estimated opera- wind directed toward or away from the center of tional sensitivity (approximately 232 to 230 dBZ)is each MCS. insufficient to see clouds with small IWC, CALIPSO c. Compositing technique can detect cloud layers with optical depths as low as 0.01 or less (Winker et al. 2007). In this study, we For the results presented in section 5, each MCS in the take advantage of the complementary capabilities of YH10 database is assigned a separate coordinate sys- these instruments, using CloudSat to examine the tem, where the origin corresponds to the center of the vertical structure of the deep convective clouds and largest raining core. Observations from CALIPSO, CALIPSO to analyze the extent of cirrus shields in MLS, and ERA-Interim are extracted and stored in the TTL. system-relative coordinates—that is, as a function of DECEMBER 2015 V I R T S A N D HOUZE 4743

FIG. 1. Density of (top) small SMCSs, (middle) large SMCSs, and (bottom) CMCSs for the years 2007–10, expressed as the number of systems observed in some portion of a 18318 grid box per year (note that the color scales differ). vertical height or pressure level and distance from the convection over the tropical continents and the warm MCS. For MLS and ERA-Interim, the dataset pressure pool. Clouds and IWC are more strongly focused in levels and horizontal sampling frequency (every 1.58) areas with high annual-mean precipitation and frequent are retained. As described earlier, CALIPSO observa- MCS occurrence (Fig. 1). Water vapor and relative hu- tions are available at higher spatial resolution; however, midity exhibit more diffuse distributions, indicating the in order to reduce noise, we calculate cloud occurrence advection of water vapor beyond the areas of frequent for 0.25 km vertical and 18 horizontal bins. No further deep convection as well as the effect of the large-scale spatial filtering is applied to the data. When indicated, temperature field. anomalies are calculated by subtracting the climatolog- The climatological, zonal-mean vertical structure of ical monthly mean at each location prior to compositing. tropical clouds and moisture is shown in Fig. 3. In the upper troposphere, clouds and ice are most frequent around 78N—the latitude of the rising branch of the 3. Climatology of MCSs and the TTL annual-mean Hadley circulation. A secondary maxi- The annual distribution in Fig. 1 indicates that all mum in IWC is observed south of the equator, associated three types of MCSs are observed throughout the deep with convection in the austral summer Hadley cell. IWC convective regions of the tropics—the tropical conti- concentrations decrease with height above ;170 hPa, nents, the intertropical convergence zones (ITCZs), and while cloud occurrence broadens and becomes more ubiquitously over the western Pacific warm pool. How- equatorially symmetric in the TTL, where cirrus can ever, as shown by YH10, both large and small SMCSs form either through the shearing off of convective anvils are more frequently observed over land or near-coastal or in situ in rising, adiabatically cooling layers. regions, while the merging of systems to produce Water vapor concentrations decrease by about three CMCSs occurs mostly over the , especially the orders of magnitude from the upper troposphere to the warm pool. lower stratosphere (Fig. 3c). This strong decrease with Annual-mean distributions of clouds and moisture at height reflects the temperature dependence of the 147 hPa, near TTL base, are shown in Fig. 2. Maxima in saturation mixing ratio, particularly at the low tem- each field are located over the regions of intense peratures observed in the TTL. The decrease in the 4744 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 72

23 FIG. 2. Annual-mean (a) CALIPSO cloud fraction, (b) MLS IWC (mg m ), (c) MLS water vapor mixing ratio (ppmv), and (d) MLS-relative humidity with respect to ice (%) at 147 hPa. Relative humidity values have been corrected for MLS temperature bias [see section 2a(4)]. infusion of water vapor into the environment by con- and Houze 1995; Masunaga et al. 2008). Cloud micro- vection also contributes to the water vapor distribution at physical processes can then be inferred from the re- higher altitudes. Despite the low water vapor concen- flectivity distributions. While YHH11 analyzed CFADs trations in the TTL, relative humidities are high (Fig. 3d), of anvil characteristics as functions of geographic re- with zonal-mean values in excess of 80% because of the gion, anvil thickness, and distance from the MCS center, extremely low temperatures. The high relative humidities we examine reflectivity characteristics for both the in the TTL promote cirrus occurrence. precipitating portions and the nonprecipitating anvils of Before considering the contribution of MCSs to each category of MCS (Fig. 4). The primary features of clouds and water vapor in the TTL, it is useful to analyze the CFADs are shared by each type of MCS: in the their characteristics in the troposphere, specifically their precipitating areas, peak reflectivities of 10–15 dBZ are vertical structure and microphysical attributes. In the observed in the melting layer (around 4–5 km). Re- following section, these characteristics are investigated flectivities decrease below the melting layer as a result using CloudSat observations. of signal attenuation by the precipitation. High re- flectivities dominate the distribution through most of the troposphere, indicating the presence of large drop- 4. Vertical structure of MCSs in the troposphere lets in the convective updraft below the melting layer The distribution of CloudSat reflectivities in MCSs and the formation of graupel and large ice particles as a function of height can be represented using con- above. In the upper troposphere, above 9–10 km, re- toured frequency by altitude diagrams (CFADs; Yuter flectivities generally decrease with height because of DECEMBER 2015 V I R T S A N D HOUZE 4745

anvils are lower, generally below 10 dBZ (Fig. 4). Peak reflectivities are observed in the lower portions of MCS anvils, around 6–10 km, and the median reflectivity de- creases with height to approximately 220 dBZ near 13 km. Vertical velocities in anvils are generally weak, such that ice particles tend to drift downward (Houze 2014, chapter 6). The presence of larger particles near anvil cloud base reflects this sedimentation, as well as ice particle growth by aggregation (McFarquhar and Heymsfield 1996). YHH11 further noted weaker re- flectivities and a narrowing of the reflectivity distribu- tion in anvils as distance from the MCS core increased, consistent with aging. While the CFADs for the three MCS categories share these basic similarities, key differences also exist. In the raining cores of small SMCSs, the reflectivity distribu- tion exhibits a more vertical orientation, with median reflectivities remaining high (greater than 0 dBZ)upto 13 km. In contrast, in large and connected MCSs, re- flectivities decrease more rapidly with height. To more clearly contrast the MCS types, difference plots are shown in Fig. 5. These plots were constructed by sub- tracting the CFADs for small SMCSs from those for large SMCSs, such that blue shading indicates the re- flectivities proportionately more likely to be observed in small systems. In Fig. 5, above the melting layer, high reflectivities are more commonly observed in small SMCSs, indicative of strong convection and less de- veloped stratiform regions. The reverse is observed above the melting layer in large SMCSs, where the narrower reflectivity distribution is consistent with the presence of stratiform rain areas. A narrower distribu- tion is also observed in the anvils of large SMCSs com- pared to those in small SMCSs, as shown by the pattern of blue shading in Fig. 5b. The greatest contrast is ob- served in the occurrence of high reflectivities above 9 km, which are much more common in the small SMCSs, as the larger particles formed in the intense convective updrafts find their way into the anvils of the small MCSs.

FIG. 3. Annual-mean, zonal-mean (a) CALIPSO cloud fraction, 23 (b) MLS ice water content (mg m ), (c) MLS water vapor mixing 5. MCS impacts on the TTL ratio (ppmv; note the log scale), and (d) MLS-relative humidity with respect to ice (%). Relative humidity values have been cor- Having examined MCS vertical structure in the tro- rected for MLS temperature bias [see section 2a(4)]. For (a), the posphere, we now investigate the impact of MCSs on the vertical coordinate is given in both pressure (hPa) and altitude (km, rounded to the nearest 0.5 km). overlying TTL. In this section, we present analyses of clouds, wind, IWC, and water vapor in the vicinity of the MCSs, with the variables composited in system-relative lower concentrations of particles of generally smaller coordinates as described in section 2c. dimension. a. Clouds and wind As air from the MCS raining cores is incorporated into the anvils, microphysical processes modify the Cross sections of composited cloud fraction in the characteristics of the ice particles. Reflectivities in MCS upper troposphere and TTL near MCSs are shown in 4746 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 72

FIG. 4. CFADs of CloudSat reflectivities in (a),(d) small SMCSs, (b),(e) large SMCSs, and (c),(f) CMCSs, using profiles sampling (top) the precipitating portions of the MCSs and (bottom) the anvil portions. Each CFAD is normalized such that its maximum value is 1.

the left column of Fig. 6. The basic cloud distribution is SMCS compared to a small SMCS. These differences similar for all MCS types: while the largest cloud frac- lessen with distance, and at 1500 km from the MCS, tions are observed lower in the troposphere, cloud oc- composite cloud fractions are ;0.2 in each case. At such currence is frequent in the upper troposphere near the large distances, CALIPSO is likely observing neigh- MCS center. About 74% of CALIPSO profiles within boring clouds. To reduce the impact of clouds that may 100 km of the center of an MCS observed a cloud top not be part of the MCSs, we subtract the climatological above 150 hPa, at TTL base, and about half of such monthly mean cloud fraction at each grid point prior to layers were opaque to the lidar (not shown). A cloud compositing. Anomalous cloud fractions for each MCS top above 120 hPa was observed in about 47% of the type are shown in the right column of Fig. 6. Here, the profiles. Above that level, cloud occurrence decreases contrast between small SMCSs and large and connected more rapidly with height. Beginning within about MCSs is even more striking. Enhanced cloudiness near 200 km of the MCS center and extending farther away, small SMCSs is observed through the depth of the TTL an outward bulge in cloud fraction indicates the but is confined to within about 200 km of the MCS spreading of cirrus anvils. This bulge is centered center, and the magnitude of the anomalies is small: less ;170 hPa near the MCS core and shifts upward to than 0.1. In large and connected MCSs, however, ;150 hPa at greater distances. cloudiness anomalies are larger and extend higher in the Large and connected MCSs are associated with sub- TTL. This result is in agreement with Rossow and Pearl stantially greater cloud fractions in the upper tropo- (2007), who reported an increasing likelihood of con- sphere and TTL than small SMCSs, particularly in the vective intrusions of the TTL with increasing cloud area within a few hundred kilometers of the MCS center, system size and that such intrusion was observed in al- where the probability of CALIPSO observing a cloud in most all systems with radii . 500 km. We also observe the upper troposphere is almost double near a large enhanced cloud occurrence to distances of 600 km and DECEMBER 2015 V I R T S A N D HOUZE 4747

MCSs are associated with both stronger and broader ascent than small SMCSs, and their anomalous di- vergence extends into the cloud-free environment. b. Ice Ice associated with MCSs can be derived from MLS observations, and IWC anomalies calculated as de- scribed in section 2c are shown in Fig. 7. IWC anomalies for the subset of levels above TTL base are shown in the right column. The largest ice concentrations and anomalies in these cross sections are in the upper tro- posphere near the MCS convective cores and extend upward into the TTL. Large and connected MCSs inject larger quantities of ice into the upper troposphere and

FIG. 5. Difference between CFADs for large and small SMCSs, TTL than small SMCSs. IWC decreases with distance using profiles sampling (a) the precipitating portions of the MCSs from the MCS core, indicating the removal of larger ice and (b) the anvil portions. Blue shading indicates where reflectivity particles via sedimentation and/or sublimation in the values are proportionately more likely to be observed in divergent airflow associated with the stratiform and a small SMCS. anvil regions. While peak enhanced cloudiness in the MCS outflow lies at or above TTL base, as illustrated by the overlaid cloud fraction anomalies, the largest IWC beyond in large and connected MCSs. The outward values are at the bottom of the cross section. This is bulge in the anomaly patterns is at a slightly higher al- partly a result of the relatively coarse vertical resolution titude than that in the composites in Fig. 6a, indicating of the IWC data. In addition, CALIPSO’s sensitivity that the MCS anvil cirrus and outflow are slightly higher enables it to detect cirrus layers with IWC too small to than what is typically observed in the tropics. be detected by MLS. As discussed in section 1, various studies have ana- lyzed how frequently deep convection reaches into the c. Water vapor tropical stratosphere. Our analysis indicates that 11% of CALIPSO profiles within 100 km of the center of an As shown in Fig. 2, water vapor near TTL base is more MCS indicate a cloud top above 95 hPa, near the CPT, broadly distributed than ice. The strong decrease of and approximately 1% indicate a cloud top above water vapor concentrations with altitude in the TTL 80 hPa, near the top of the TTL. Of the cloud layers dominates cross-sectional composites of water vapor with tops above 95 hPa within 100 km of the center of relative to MCS location (not shown). It is more in- an MCS, 65% have a sufficiently low optical depth that formative to present the composites as anomalies from they are not opaque to the lidar. These nonopaque the climatological monthly mean (Fig. 8). Each MCS layers have a median thickness of 1.8 km (not shown). type is associated with enhanced water vapor in the Thus, while some portion of convection in MCSs ex- upper troposphere. The anomalies decrease with dis- tends into the lower stratosphere, the majority of the tance from the MCS but extend radially outward beyond extremely high-topped clouds observed by CALIPSO the area of substantially enhanced cloudiness or IWC near MCSs are thin anvil clouds or cirrus layers over- (cf. Fig. 7). As noted above, the anomalous upper- lying the deep convective clouds, as discussed in tropospheric divergence also extends beyond the Garrett et al. (2004). cloudy region, particularly in large and connected The overlaid vectors in Fig. 6 show the radial winds MCSs. The true strength and extent of the MCS outflow and vertical velocities from ERA-Interim. Strong ascent can only be estimated from the reanalysis winds. Be- is observed in the convective core of the MCS, and its cause A-Train orbits repeat only every 16 days, we also magnitude is underestimated as a result of the relatively cannot determine whether premoistening of the upper coarse resolution of the reanalysis data. The vertical troposphere occurred prior to MCS formation. How- velocities weaken rapidly with height in the stably ever, the results in Fig. 8 indicate that water vapor stratified TTL. Upper-level divergence and outflow spreads laterally from the MCSs and into the sur- from the MCS are reflected in the outward direction of rounding cloud-free region prior to MCS dissipation. the composite radial winds. A similar pattern is ob- Previous studies have shown that diverging cirrus can served in the anomaly wind fields. Large and connected persist for up to 2 days after convective dissipation and 4748 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 72

FIG. 6. (left) Composite CALIPSO cloud fraction [colors and contours; contour interval (CI) 5 0.1] and ERA- Interim vertical velocity and horizontal radial wind (vectors) in the upper troposphere and TTL as a function of distance from MCS center, for (a),(b) small SMCSs, (c),(d) large SMCSs, and (e),(f) CMCSs. (right) Anomalies of cloud fraction (CI 5 0.05) and wind relative to the MCS center. be advected up to 1000 km and that elevated upper- occurs in response to strong divergence in the cloud tropospheric water vapor concentrations persist longer outflow. They point out that these water vapor patterns than the cirrus (Luo and Rossow 2004). are observed with any type of cloud in the TTL: deep Water vapor values in the TTL are too small to be seen convective, convectively generated cirriform, or in situ in the left column of Fig. 8, which is dominated by the cirriform. It is therefore reasonable that the processes moister upper troposphere. Anomalies for the subset of they describe should be associated with the MCSs dis- altitudes in the TTL (above ;150 hPa) are shown in the cussed herein, thus producing the negative water vapor right column of Fig. 8. Contrasting characteristics are anomalies in the lower TTL in Fig. 8. The cirriform upper- observed near the MCS centers (within a radius of about level clouds of MCSs and their anvils are produced by the 200 km) and at distances greater than about 400 km. Near combination of convective and mesoscale dynamics that the convective core, water vapor concentrations are be- characterize MCSs (Houze 2004, 2014, chapter 9). We low the climatological mean in the lowermost TTL, par- further observe weak moistening in the upper TTL above ticularly in large and connected MCSs. This altitude can the MCS centers, consistent with the observation of Corti be above, at, or below the cloud top, as indicated by the et al. (2008) that the deep clouds that extend upward overlaid cloud fraction anomalies and the discussion in beyond the CPT leave behind ice particles which sub- section 5a. Chae et al. (2011) noted rapidly decreasing limate and moisten the environment. moistening within clouds above about 14 km and weak Outside the convective core, at distances of about dehydration just below cloud-top level when the cloud 400 km from the MCS center and extending outward, extended above 16 km, as the environmental relative enhanced water vapor concentrations are observed in humidity became supersaturated. They also noted strong the lower TTL as well as the upper troposphere (Fig. 8), dry anomalies in the lower TTL when those altitudes lay indicating that outflow from the MCSs and evaporative just above the convective cloud tops, where subsidence processes are moistening the environment. The water DECEMBER 2015 V I R T S A N D HOUZE 4749

23 FIG. 7. (left) Composite MLS ice water content anomalies (colored shading; mg m ) as a function of distance from MCS center, for (a),(b) small SMCSs, (c),(d) large SMCSs, and (e),(f) CMCSs, overlaid with cloud fraction (contours; CI 5 0.05) and circulation anomalies (vectors) from Fig. 6. (right) As in (left), but for subset of altitudes in the TTL. vapor anomalies in this region decrease more rapidly moistening of a broad region by the smallest MCSs. with height than the environmental water vapor (cf. with Rather, it appears that small SMCSs tend to form in Fig. 3), consistent with decreased convective influence large-scale environments with a slightly moistened TTL. and less effective moistening in the TTL. The moisten- ing appears to extend upward to ;120 hPa; however, 6. Conclusions observations with finer vertical resolution would be needed to examine the precise vertical extent of the In this study, the vertical structure and microphysical moistening. As in the upper troposphere, enhanced composition of MCSs, as well as their role in detraining water vapor in the TTL extends laterally into cloud-free ice particles and water vapor in the TTL, have been regions. Weak dry anomalies are observed near the investigated using a suite of satellite and reanalysis ob- CPT, near the tops of the highest clouds (cf. with Fig. 6). servations. We have examined three types of MCSs, as The water vapor anomaly pattern near small SMCSs identified in observations from MODIS and AMSR-E exhibits different characteristics than that for large and (YH10): small and large separated MCSs and connected connected MCSs. Primarily, it is less robust, indicating systems mostly likely formed via the merging of multiple that it is the largest MCSs that are important to the MCSs. MCSs are ubiquitous in the convective regions of moistening of these upper levels. Positive anomalies the tropics (Fig. 1), including East Asia, the western associated with the small SMCSs are observed through Pacific warm pool, and Central and South America, most of the TTL, with peak values at 121 hPa. Laterally, which have been identified as major source regions the anomalies extend beyond the narrow core of en- for air parcels eventually entering the stratosphere hanced cloudiness and several hundred kilometers into (Schoeberl et al. 2013). SMCSs occur more frequently the region of suppressed cloudiness (i.e., beyond the red over the tropical landmasses, while the merging of line in Fig. 8b). As previously noted, the radial outflow longer-lived systems to form CMCSs is more common in associated with small SMCSs is weak; hence, it is un- the oceanic ITCZs and over the western Pacific warm realistic to interpret the anomaly pattern in Fig. 8b as pool (Fig. 1). 4750 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 72

FIG.8.AsinFig. 7, but colored shading is MLS water vapor mixing ratio anomalies (ppmv).

CloudSat reflectivity distributions in the precipitating but also a prominent outward bulge in cloud occur- regions of active MCSs indicate high reflectivities in the rence near TTL base associated with anvil cirrus. Most midtroposphere, indicative of large ice and graupel MCSs are associated with clouds in the TTL: 74% of particles forming in the convective updrafts and evolv- CALIPSO profiles within 100 km of the center of an ing into stratiform precipitation as the updrafts weaken MCS exhibited a cloud top above TTL base (about or being detrained into neighboring stratiform regions 150 hPa; Fig. 6). Large and connected MCSs are asso- of the MCSs (Fig. 4). Reflectivities in the MCS anvils are ciated with more frequent cloud occurrence in the TTL generally lower and decrease rapidly with altitude in the than small SMCSs, in agreement with the inference of upper troposphere, consistent with sedimentation and Rossow and Pearl (2007) that the likelihood of cloud aggregation within the anvils and the removal of larger penetrating upward into high levels increases with graupel prior to air entering the anvil, as previously system size. Rossow and Pearl (2007) and Takahashi shown by Cetrone and Houze (2009) and YHH11. Small and Luo (2014) both report that penetration of the TTL SMCSs exhibit characteristics of intense, young con- tends to occur in the growing stage of the cloud system, vection and are consequently more likely than large and although the latter noted significant penetration during connected MCSs to contain large ice and graupel par- the mature stage. Our results associate the most ex- ticles in the upper troposphere (Fig. 5). Larger MCSs tensive cloudiness in the TTL specifically with large exhibit narrower, more strongly tilted reflectivity dis- MCSs, in which anvils and overlying cirrus occur in tributions in the upper troposphere, consistent with the greater quantities. ‘‘Overshooting’’ convective ele- ‘‘size sorting’’ processes characteristic of mature anvils ments are actually less likely in these more widespread (McFarquhar and Heymsfield 1996). MCSs; for example, Virts and Houze (2015) showed Composites of cloud occurrence observed by that the larger MCSs have less vigorous convective CALIPSO near MCSs (Fig. 6) reveal that cloudiness elements, as indicated by occurrence. It is the anomalies in small SMCSs are primarily vertically ori- combination of convective and stratiform dynamics of ented and confined near the MCS cores, while large and MCSs that make them robust in affecting cloudiness in connected MCSs exhibit not only broader core regions the TTL. DECEMBER 2015 V I R T S A N D HOUZE 4751

MCSs inject ice into the upper troposphere and TTL, greater importance than isolated overshooting cells be- as shown by composites of anomalous MLS IWC with cause they are large enough in horizontal scale to be respect to MCS location (Fig. 7). The largest ice quan- highly impactful on the water vapor field. Quantification tities are near the centers of the MCSs, but elevated ice of the comparative contributions of convection of vari- concentrations also extend laterally into the stratiform ous sizes and organization to the total water vapor and anvil regions. Large and connected MCSs contrib- content of the TTL and lower stratosphere is needed but ute the greatest quantities of ice to the TTL. These is not straightforward given the water vapor data cur- systems contain predominately smaller ice crystals and rently available. less graupel than the small SMCSs (Fig. 5). The less MCS evolution cannot be examined using A-Train rapid sedimentation rates of the smaller particles allow satellite observations, which provide only snapshots in them to remain aloft in the stratiform and anvil regions. time. However, the results presented in this paper sug- We note that the IWC anomalies decrease nearly gest that small SMCSs are primarily young, developing monotonically with height in the upper troposphere and systems, while large and connected MCSs represent the TTL, while the cloud occurrence exhibits an outward mature stage of the MCS life cycle (Houze 2004, 2014, bulge associated with the spreading anvils. This differ- chapter 9). Small SMCSs contain vigorous convection ence is largely due to the coarser vertical resolution of that produces large ice particles, especially graupel. the IWC data (approximately 4 km in the TTL; Wu et al. These systems have limited stratiform and anvil regions, 2008) but is also consistent with lower ice concentrations but the anvils that do exist contain large particles de- in the higher cirrus layers observed by CALIPSO. trained from the convective core. Water vapor also de- Composites based on MLS water vapor mixing ratios trains from the convective region, moistening the upper (Fig. 8) demonstrate that MCSs moisten the upper tro- troposphere, while reduced moistening in convection posphere. Moreover, stronger moistening is associated penetrating the TTL combined with subsidence above with the large and connected MCSs. In the TTL, a cloud top produces negative water vapor anomalies near complex structure is observed that is consistent with TTL base. As MCSs mature (and in some cases merge), previous studies of dehydration and hydration at upper their convection becomes less vigorous, resulting in re- levels. Near the MCS cores, water vapor concentrations duced ice and graupel size. This ice is incorporated into are suppressed near TTL base. Higher relative humidity the extensive stratiform and anvil regions (either by in the TTL (Fig. 3) reduces the possible amount of detrainment or by cells weakening and becoming strat- moistening, and subsidence above cloud top has been iform), where it undergoes sedimentation and aggrega- observed to produce strong drying in the lower TTL, as tion, such that smaller ice particles dominate near cloud shown by Chae et al. (2011). Weak moistening is ob- top. At this stage, MCSs are characterized by extensive served in the upper TTL, near the CPT, as the deepest upper-level divergence. Anvils extend radially outward penetrating convection behind small ice crystals via this divergence, with peak cloud occurrence near that sublimate and moisten the environment. A reversed TTL base shifting to slightly higher altitudes with dis- pattern is observed from approximately 400 km out- tance from the MCS core. The MCSs thus moisten the ward, where the spreading anvils of large and connected upper troposphere via sublimation, and water vapor is MCSs moisten the lower TTL, and weak drying is ob- detrained laterally into the cloud-free air. At larger served at the CPT. The water vapor and radial wind distances from the convective core, where the above- anomalies in the upper troposphere and lower TTL ex- cloud subsidence is weaker, this moistening extends into tend into the cloud-free areas, indicating that the MCS the TTL. outflow is moistening the environment. We also em- Our discussion has focused on contrasts between phasize the possibility that a premoistened environment small SMCSs and large and connected MCSs. More favors the development of mature MCS stratiform and subtle differences exist between large SMCSs and anvil regions, but this cannot be investigated here owing CMCSs. As shown in Fig. 1, large SMCSs occur more to the limitations of satellite data. often over or near land, while most CMCSs are observed Our conclusion that large and connected MCSs more over the oceans. Thus, in a broad , differences strongly moisten the upper troposphere and lower TTL between these two MCS categories represent the dif- than small SMCSs contrasts with the notion that ferences between the largest MCSs over land and moistening of the TTL is primarily caused by isolated oceans. Examining Figs. 4 and 6–8, it can be seen that overshooting convective cells. The area covered by compared to large SMCSs, CMCSs are slightly larger overshooting convective cells is many orders of magni- (Yuan and Houze 2013), exhibit somewhat stronger tude smaller than the area occupied by MCSs. Our re- outflow in the upper troposphere, and are associated sults indicate that the largest, mature MCSs are of with larger ice and water vapor anomalies. However, 4752 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 72 these differences are slight, and our results indicate that Dessler, A. E., S. P. Palm, and J. D. Spinhirne, 2006: Tropical mature MCSs, whether they form over land or ocean, cloud-top height distributions revealed by the Ice, Cloud, and are major contributors to water vapor in the upper tro- Land Elevation Satellite (ICESat)/Geoscience Laser Altime- ter System (GLAS). J. Geophys. Res., 111, D12215, posphere and TTL. doi:10.1029/2005JD006705. Folkins, I., and R. V. Martin, 2005: The vertical structure of tropical Acknowledgments. The authors thank three reviewers convection and its impact on the budgets of water vapor and for their helpful suggestions and Beth Tully for processing ozone. J. Atmos. Sci., 62, 1560–1573, doi:10.1175/JAS3407.1. the graphics. This work was supported by the National Fu, Q., Y. Hu, and Q. Yang, 2007: Identifying the top of the tropical tropopause layer from vertical mass flux analysis and CALIPSO Aeronautics and Space Administration (Grants lidar cloud observations. Geophys.Res.Lett., 34, L14813, NNX10AM28G and NNX13AQ37G) and the U.S. De- doi:10.1029/2007GL030099. partment of Energy (Grant DE-SC008452). CALIPSO Fueglistaler, S., A. E. Dessler, T. J. Dunkerton, I. Folkins, Q. Fu, data are available from the NASA Langley Research and P. W. Mote, 2009: Tropical tropopause layer. Rev. Geo- Center Atmospheric Science Data Center (http://www- phys., 47, RG1004, doi:10.1029/2008RG000267. Garrett, T. J., A. J. Heymsfield, M. J. McGill, B. A. Ridley, D. G. CloudSat calipso.larc..gov/), data are available from Baumgardner, T. P. Bui, and C. R. Webster, 2004: Convective the CloudSat Data Processing Center (http://www. generation of cirrus near the tropopause. J. Geophys. Res., cloudsat.cira.colostate.edu/), MLS data are available 109, D21203, doi:10.1029/2004JD004952. from the Goddard Earth Sciences (GES) Data and In- Gettelman, A. E., M. L. Salby, and F. Sassi, 2002: Distribution and formation Services Center (DISC; http://mirador.gsfc. influence of convection in the tropical tropopause region. J. Geophys. Res., , 4080, doi:10.1029/2001JD001048. nasa.gov/), ERA-Interim fields are available from the 107 Hartmann, D. L., H. H. Hendon, and R. A. Houze Jr., 1984: European Centre for Medium-Range Weather Forecasts Some implications of the mesoscale circulations in tropi- (http://ecmwf.int), and MODIS and AMSR-E data are cal cloud clusters for large-scale dynamics and . available from the NASA data centers (http://ladsweb. J. Atmos. Sci., 41, 113–121, doi:10.1175/1520-0469(1984)041,0113: nascom.nasa.gov/ and http://nsidc.org/data/amsre). SIOTMC.2.0.CO;2. Hassim, M. E. E., and T. P. Lane, 2010: A model study on the influence of overshooting convection on TTL water vapour. REFERENCES Atmos. Chem. Phys., 10, 9833–9849, doi:10.5194/ Alcala, C. M., and A. E. Dessler, 2002: Observations of deep acp-10-9833-2010. convection in the tropics using the Tropical Rainfall Measur- Highwood, E. J., and B. J. Hoskins, 1998: The tropical tropopause. ing Mission (TRMM) precipitation radar. J. Geophys. Res., Quart. J. Roy. Meteor. Soc., 124, 1579–1602, doi:10.1002/ 107, 4792, doi:10.1029/2002JD002457. qj.49712454911. Anthes, R. A., and Coauthors, 2008: The COSMIC/FORMOSAT- Houze, R. A., Jr., 1982: Cloud clusters and large-scale vertical 3 mission: Early results. Bull. Amer. Meteor. Soc., 89, 313–333, motions in the tropics. J. Meteor. Soc. Japan, 60, 396–410. doi:10.1175/BAMS-89-3-313. ——, 1989: Observed structure of mesoscale convective systems Brewer, A. W., 1949: Evidence for a world circulation provided by and implications for large-scale heating. Quart. J. Roy. Meteor. the measurements of helium and water vapor distribution in Soc., 115, 425–461, doi:10.1002/qj.49711548702. the stratosphere. Quart. J. Roy. Meteor. Soc., 75, 351–363, ——, 2004: Mesoscale convective systems. Rev. Geophys., 42, doi:10.1002/qj.49707532603. RG4003, doi:10.1029/2004RG000150. Cetrone, J., and R. A. Houze Jr., 2009: Anvil clouds of tropical ——, 2014: Cloud Dynamics. 2nd ed. Elsevier, 432 pp. mesoscale convective systems in monsoon regions. Quart. Jensen, E. J., A. S. Ackerman, and J. A. Smith, 2007: Can over- J. Roy. Meteor. Soc., 135, 305–317, doi:10.1002/qj.389. shooting convection dehydrate the tropical tropopause layer? Chaboureau, J.-P., J.-P. Cammas, J. Duron, P. J. Mascart, N. M. J. Geophys. Res., 112, D11209, doi:10.1029/2006JD007943. Sitnikov, and H.-J. Voessing, 2007: A numerical study of tropical Jiang, J. H., and Coauthors, 2012: Evaluation of cloud and water cross-tropopause transport by convective overshoots. Atmos. vapor simulations in CMIP5 climate models using NASA Chem. Phys., 7, 1731–1740, doi:10.5194/acp-7-1731-2007. A-Train satellite observations. J. Geophys. Res., 117, D14105, Chae, J. H., D. L. Wu, W. G. Read, and S. C. Sherwood, 2011: The role doi:10.1029/2011JD017237. of tropical deep convective clouds on temperature, water vapor, Khaykin, S., and Coauthors, 2009: Hydration of the lower and dehydration in the tropical tropopause layer (TTL). Atmos. stratosphere by ice crystal over land convective sys- Chem. Phys., 11, 3811–3821, doi:10.5194/acp-11-3811-2011. tems. Atmos. Chem. Phys., 9, 2275–2287, doi:10.5194/ Corti, T., and Coauthors, 2008: Unprecedented evidence for deep acp-9-2275-2009. convection hydrating the tropical stratosphere. Geophys. Res. Kim, J., and S.-W. Son, 2012: Tropical cold-point tropopause: Lett., 35, L10810, doi:10.1029/2008GL033641. Climatology, seasonal cycle, and intraseasonal variability de- Danielsen, E. F., 1993: In situ evidence of rapid, vertical, irre- rived from COSMIC GPS radio occultation measurements. versible transport of lower tropospheric air into the lower J. Climate, 25, 5343–5360, doi:10.1175/JCLI-D-11-00554.1. tropical stratosphere by convective cloud turrets and by Kummerow, C., and R. Ferraro, 2006: EOS/AMSR-E level-2 larger-scale upwelling in tropical . J. Geophys. Res., rainfall. NASA Algorithm Theoretical Basis Doc., 10 pp. 98, 8665–8681, doi:10.1029/92JD02954. [Available online at http://www.ghcc.msfc.nasa.gov/AMSR/ Dee, D. P., and Coauthors, 2011: The ERA-Interim reanalysis: atbd/rainfall2atbd.doc.] Configuration and performance of the data assimilation system. ——, and Coauthors, 2001: The evolution of the Goddard profiling Quart. J. Roy. Meteor. Soc., 137, 553–597, doi:10.1002/qj.828. algorithm (GPROF) for rainfall estimation from passive DECEMBER 2015 V I R T S A N D HOUZE 4753

microwave sensors. J. Appl. Meteor., 40, 1801–1820, doi:10.1175/ measurements. J. Geophys. Res., 113, D15S11, doi:10.1029/ 1520-0450(2001)040,1801:TEOTGP.2.0.CO;2. 2007JD008783. L’Ecuyer, T. S., and J. H. Jiang, 2010: Touring the Sherwood, S. C., and A. E. Dessler, 2000: On the control of aboard the A-Train. Phys. Today, 63, 36–41, doi:10.1063/ stratospheric humidity. Geophys. Res. Lett., 27, 2513–2516, 1.3463626. doi:10.1029/2000GL011438. Liu, C., and E. J. Zipser, 2005: Global distribution of convection Soden, B. J., 2004: The impact of tropical convection and cirrus on penetrating the tropical tropopause. J. Geophys. Res., 110, upper tropospheric humidity: A Lagrangian analysis of satel- D23104, doi:10.1029/2005JD006063. lite measurements. Geophys. Res. Lett., 31, L20104, doi:10.1029/ Livesey, N. J., and Coauthors, 2011: Version 3.3 and 3.4 level 2 data 2004GL020980. quality and description document. Jet Propulsion Laboratory Stephens, G. L., and Coauthors, 2002: The CloudSat mission and JPL D-33509, 164 pp. [Available online at https://mls.jpl.nasa.gov/ the A-train. Bull. Amer. Meteor. Soc., 83, 1771–1790, doi:10.1175/ data/v3_data_quality_document.pdf.] BAMS-83-12-1771. Luo, Z., and W. B. Rossow, 2004: Characterizing tropical cirrus life Takahashi, H., and Z. Luo, 2012: Where is the level of neutral cycle, evolution, and interaction with upper-tropospheric water buoyancy for deep convection? Geophys. Res. Lett., 39, vapor using Lagrangian trajectory analysis of satellite ob- L15809, doi:10.1029/2011GL050692. servations. J. Climate, 17, 4541–4563, doi:10.1175/3222.1. ——, and ——, 2014: Characterizing tropical overshooting deep ——, D. Kley, and R. H. Johnson, 2007: Ten years of measurements convection from joint analysis of CloudSat and geostationary of tropical upper-tropospheric water vapor by MOZAIC. Part satellite observations. J. Geophys. Res. Atmos., 119, 112–121, I: Climatology, variability, transport, and relation to deep doi:10.1002/2013JD020972. convection. J. Climate, 20, 418–435, doi:10.1175/JCLI3997.1. Virts, K. S., and R. A. Houze Jr., 2015: Variation of lightning and Mapes, B. E., and R. A. Houze Jr., 1993: Cloud clusters and superclu- convective rain fraction in mesoscale convective systems of sters over the oceanic warm pool. Mon. Wea. Rev., 121, 1398–1415, the MJO. J. Atmos. Sci., 72, 1932–1944, doi:10.1175/ doi:10.1175/1520-0493(1993)121,1398:CCASOT.2.0.CO;2. JAS-D-14-0201.1. Marchand, R., G. G. Mace, T. Ackerman, and G. Stephens, 2008: ——, J. M. Wallace, Q. Fu, and T. P. Ackerman, 2010: Tropical Hydrometeor detection using Cloudsat—An earth-orbiting tropopause transition layer cirrus as represented by CALIPSO 94-GHz cloud radar. J. Atmos. Oceanic Technol., 25, 519–533, lidar observations. J. Atmos. Sci., 67, 3113–3129, doi:10.1175/ doi:10.1175/2007JTECHA1006.1. 2010JAS3412.1. Massie, S., A. Gettelman, W. Randel, and D. Baumgardner, 2002: Wilheit, T., C. D. Kummerow, and R. Ferraro, 2003: Rainfall al- Distribution of tropical cirrus in relation to convection. gorithms for AMSR-E. IEEE Trans. Geosci. Remote Sens., 41, J. Geophys. Res., 107, 4591, doi:10.1029/2001JD001293. 204–214, doi:10.1109/TGRS.2002.808312. Masunaga, H., M. Satoh, and H. Miura, 2008: A joint satellite and Williams, M., and R. A. Houze Jr., 1987: Satellite-observed global cloud-resolving model analysis of a Madden-Julian characteristics of winter monsoon cloud clusters. Mon. Wea. Oscillation event: Model diagnosis. J. Geophys. Res., 113, Rev., 115, 505–519, doi:10.1175/1520-0493(1987)115,0505: D17210, doi:10.1029/2008JD009986. SOCOWM.2.0.CO;2. McFarquhar, G. M., and A. J. Heymsfield, 1996: Microphysical Winker, D. M., W. H. Hunt, and M. J. McGill, 2007: Initial per- characteristics of three anvils sampled during the Central formance assessment of CALIOP. Geophys. Res. Lett., 34, Equatorial Pacific Experiment. J. Atmos. Sci., 53, 2401–2423, L19803, doi:10.1029/2007GL030135. doi:10.1175/1520-0469(1996)053,2401:MCOTAS.2.0.CO;2. Wright, J. S., R. Fu, and A. J. Heymsfield, 2009: A statistical Randel, W. J., and E. J. Jensen, 2013: Physical processes in the analysis of the influence of deep convection on water vapor tropical tropopause layer and their roles in a changing climate. variability in the tropical upper troposphere. Atmos. Chem. Nat. Geosci., 6, 169–176, doi:10.1038/ngeo1733. Phys., 9, 5847–5864, doi:10.5194/acp-9-5847-2009. Read, W. G., and Coauthors, 2007: Aura Microwave Limb Sounder Wu, D. L., and Coauthors, 2008: Validation of the Aura MLS cloud

upper tropospheric and lower stratospheric H2O and relative ice water content measurements. J. Geophys. Res., 113, humidity with respect to ice validation. J. Geophys. Res., 112, D15S10, doi:10.1029/2007JD008931. D24S35, doi:10.1029/2007JD008752. Yuan, J., and R. A. Houze Jr., 2010: Global variability of mesoscale Riihimaki, L. D., and S. A. McFarlane, 2010: Frequency and convective system anvil structure from A-Train satellite data. morphology of tropical tropopause layer cirrus from CALIPSO J. Climate, 23, 5864–5888, doi:10.1175/2010JCLI3671.1. observations: Are isolated cirrus different from those con- ——, and ——, 2013: Deep convective systems observed by nected to deep convection? J. Geophys. Res., 115, D18201, A-Train in the tropical Indo-Pacific region affected by the doi:10.1029/2009JD013133. MJO. J. Atmos. Sci., 70, 465–486, doi:10.1175/JAS-D-12-057.1. Rossow, W. B., and C. Pearl, 2007: 22-year survey of tropical ——, ——, and A. J. Heymsfield, 2011: Vertical structures of anvil convection penetrating into the lower stratosphere. Geophys. clouds of tropical mesoscale convective systems observed by Res. Lett., 34, L04803, doi:10.1029/2006GL028635. CloudSat. J. Atmos. Sci., 68, 1653–1674, doi:10.1175/ Schoeberl, M. R., A. E. Dessler, and T. Wang, 2013: Modeling 2011JAS3687.1. upper tropospheric and lower stratospheric water vapor Yuter, S. E., and R. A. Houze Jr., 1995: Three-dimensional kine- anomalies. Atmos. Chem. Phys., 13, 7783–7793, doi:10.5194/ matic and microphysical evolution of Florida cumulonimbus. acp-13-7783-2013. Part II: Frequency distributions of vertical velocity, reflectivity, Schumacher, C., R. A. Houze Jr., and I. Kraucunas, 2004: The trop- and differential reflectivity. Mon. Wea. Rev., 123, 1941–1963, ical dynamical response to latent heating estimates derived from doi:10.1175/1520-0493(1995)123,1941:TDKAME.2.0.CO;2. the TRMM Precipitation Radar. J. Atmos. Sci., 61, 1341–1358, Zipser, E. J., J. C. Daniel, C. Liu, S. W. Nesbitt, and D. P. Yorty, doi:10.1175/1520-0469(2004)061,1341:TTDRTL.2.0.CO;2. 2006: Where are the most intense on Earth? Schwartz, M. J., and Coauthors, 2008: Validation of the Aura Mi- Bull. Amer. Meteor. Soc., 87, 1057–1071, doi:10.1175/ crowave Limb Sounder temperature and geopotential height BAMS-87-8-1057.