Relative influence of meteorological conditions and aerosols on the lifetime of mesoscale convective systems Sudip Chakrabortya,1, Rong Fua, Steven T. Massieb, and Graeme Stephensc a

Jackson School of Geosciences, University of Texas at Austin, Austin, TX 78712; bLaboratory for Atmospheric and Space Physics, University of Colorado Boulder, Boulder, CO 80309; and cJet Propulsion Laboratory, Pasadena, CA 91109

Edited by John H. Seinfeld, California Institute of Technology, Pasadena, CA, and approved May 12, 2016 (received for review February 4, 2016)

Using collocated measurements from geostationary and polar-orbital satellites over tropical continents, we provide a large-scale statistical assessment of the relative influence of aerosols and meteorological conditions on the lifetime of mesoscale convective systems (MCSs). Our results show that MCSs’ lifetime increases by 3–24 h when vertical wind shear (VWS) and convective available potential energy (CAPE) are moderate to high and ambient aerosol optical depth (AOD) increases by 1 SD (1σ). However, this influence is not as strong as that of CAPE, relative humidity, and VWS, which increase MCSs’ lifetime by 3–30 h, 3–27 h, and 3–30 h per 1σ of these variables and explain up to 36%, 45%, and 34%, respectively, of the variance of the MCSs’ lifetime. AOD explains up to 24% of the total variance of MCSs’ lifetime during the decay phase. This result is physically consistent with that of the variation of the MCSs’ ice water content (IWC) with aerosols, which accounts for 35% and 27% of the total variance of the IWC in convective cores and anvil, respectively, during the decay phase. The effect of aerosols on MCSs’ lifetime varies between different continents. AOD appears to explain up to 20–22% of the total variance of MCSs’ lifetime over equatorial South America compared with 8% over equatorial Africa. Aerosols over the Indian Ocean can explain 20% of total variance of MCSs’ lifetime over South Asia because such MCSs form and develop over the ocean. These regional differences of aerosol impacts may be linked to different meteorological conditions. mesoscale convective systems

| aerosols | meteorological parameters

T

he hypothesis that aerosols may delay precipitation and increase cloud lifetime of shallow marine clouds (1) has motivated many researchers to study the aerosol indirect effect on convective clouds; however, the influence of aerosols on enhancing cloud lifetime has remained under debate. Mesoscale convective systems (MCSs) are deep convective clouds that cover several hundred kilometers. Previous studies have shown that aerosols affect deep convection, in particular that aerosols increase the number of smaller size cloud condensation nuclei (CCN) (2), which weaken coagulation and coalescence that form rain droplets, and consequently delay warm rainfall (3, 4). These processes allow more cloud droplets to rise above the freezing level and increase latent heat released due to glaciation (5), resulting in stronger updraft speed, enhanced cloud ice content (6), larger anvil size (5), and higher cloud top height (7). The top of the troposphere warms owing to the aerosol-induced changes in convective anvils (8). Although these aerosol effects have been seen in observations from field campaigns (9, 10), they have been undetectable on large spatial and multiyear scales. Rosenfeld et al. (11) have attributed this lack of detectability on the large scale of aerosol invigoration of convection to its variation with meteorological conditions and to the lack of knowledge of the relative humidity (RH) outside the clouds. The variation of aerosol effects on convection with meteorological parameters has been studied previously (11). For example, model simulation has shown that an increase in aerosol concentrations up to an optimal level can invigorate the MCSs under weak vertical wind shear (VWS) and higher RH but suppress the MCSs

7426–7431 | PNAS | July 5, 2016 | vol. 113 | no. 27

under strong VWS in a dry environment (12, 13). They found that, due to a significant enhancement in the convective available potential energy (CAPE), corresponding to an increase in RH from 50% to 70%, aerosol impact on ice crystal mass becomes pronounced, with a dramatic increase in the size of the anvils and the mass of ice crystals of the deep convection. However, such impacts are negligible when RH increases from 40% to 50%, due to little increase in CAPE (13). Moreover, consumption of CAPE for a given amount of rainfall is converted to an equal amount of kinetic energy that invigorates the convection (5, 14). These studies indicate that VWS, RH, and CAPE are important factors that can influence aerosol impacts on the MCSs. However, no quantitative assessment of the relative influence of aerosols versus these meteorological parameters on convective lifetime using satellite measurements has been established (11). The influence of aerosols on the MCSs is expected to vary in different phases of the convective life cycle. For example, Rosenfeld et al. (5) hypothesized that the impact of aerosol on deep convection is stronger and more prominent during the dissipating phase. Using cloud-resolving simulations, Fan et al. (15) found out that aerosol microphysical effects intensify the deep convection during the mature and decaying phases by forming a larger number of smaller and long-lasting particles, whereas additional latent heat released due to aerosols’ thermodynamic effect is responsible for invigorating the deep convections during the growing phase. Hence, the detection of aerosol impacts might not be visible until the mature phase, as no satellite can directly measure the thermodynamic properties. Significance Mesoscale convective systems (MCSs) are the primary source of precipitation over the tropics and midlatitudes, and their lifetime can have a large influence on the variability of rainfall, especially extreme rainfall that causes flooding. The hypothesis that aerosols can increase the lifetime of the MCSs by weakening or delaying precipitation has long been proposed, but we have not known whether that increase is significant on global and regional scales, and, if so, how it compares with the influence of meteorological conditions. We use multiyear collocated geostationary and polar-orbital satellite datasets to provide, to our knowledge, the first observational assessment of such an aerosol effect and its relative importance compared with other meteorological conditions in determining the variability of MCSs’ lifetime. Author contributions: S.C., R.F., and G.S. designed research; S.C. performed research; S.C. analyzed data; and S.C., R.F., and S.T.M. wrote the paper. The authors declare no conflict of interest. This article is a PNAS Direct Submission. 1

To whom correspondence should be addressed. Email: [email protected].

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. 1073/pnas.1601935113/-/DCSupplemental.

www.pnas.org/cgi/doi/10.1073/pnas.1601935113

Results and Discussion The lifetime of the MCSs is influenced by CAPE (18–20), a useful measure for cloud buoyancy and vertical velocity (21, 22), moisture entrained from the warm atmospheric boundary layer

20 40 60 80 100 fAOD

D

0

20 40 60 80 100 RH850 [%]

-30

G

-24

-18

Lifetime_All_MCS fAOD 1% VWS 16% RH500 4%

Chakraborty et al.

Unexplained 19% CAPE 33% RH850 27%

VWS 200-900 [10-4 S-1]

20 40 60 80 100 fAOD

σfAOD*dlifetime/dfAOD

VWS 200-900 [10-4 S-1] 6 12 18 24 30 36

VWS 200-900 [10-4 S-1] 6 12 18 24 30 36

σfAOD*dlifetime/dfAOD

0

0

-12

6 12 18 24 30 36

0

B

σRH500*dlifetime/dRH500

E

0

6

12

20 40 60 80 100 fAOD

σVWS200-900*dlifetime/dVWS200-900

F

0

600 1800 3000 CAPE [J/kg]

-6

c

C

0

fAOD 0 20 40 60 80 100

A

σRH850*dlifetime/dRH850

VWS 200-900 [10-4 S-1] 6 12 18 24 30 36

VWS 200-900 [10-4 S-1] 6 12 18 24 30 36

σCAPE*dlifetime/dCAPE

(19, 20, 22), and lateral entrainment of free tropospheric moisture, which can dilute the buoyancy of the rising air and account for up to 33–50% of the rainfall (23). VWS contributes to organizing the storms, determining whether updraft and downdraft regions overlap, slantwise accent of the moist air, and precipitation (19, 24–27). We evaluate the influence of both the low-level VWS (26, 28) and the deep tropospheric VWS (see Methods). The former influences rainfall and total condensation (26), whereas the latter influences vertical velocity (26) and also MCS anvil formation (8, 26, 28–30). Past studies have shown that the response of deep convective clouds to aerosols is nonlinear. Over land, aerosol microphysical effects on deep convective clouds saturate and reverse at AOD > 0.3 (31). Moderate concentrations of aerosols (or CCNs) maximize the invigoration effect of aerosols, whereas higher aerosol concentrations can reduce the vigor of the convection (5). However, over humid land and oceans, more aerosol is required to suppress rainfall. Thus, a threshold of AOD ≈ 0.3 is likely to be optimal as a threshold of pollution to determine the influence of aerosols on the MCSs (32, 33). We use AOD > 0.15 over land and AOD > 0.3 over the ocean as the threshold for a polluted environment. To isolate the effect of aerosols on MCSs’ lifetime and compare it to the effect of these primary meteorological conditions, we evaluate how the lifetime of the MCSs changes with 1 SD (1σ) of CAPE, RH at 850 hPa (∼1.5 km above the sea level, RH850), RH at 500 hPa (∼5.5 km above the sea level, RH500), deep tropospheric VWS, and the fraction of the number of aerosol pixels with AOD (fAOD) greater than 0.15 (0.3) over land (the ocean) to the total number of AOD pixels within a range of 2° latitude/longitude from the boundary of the MCSs over all three regions (Methods and Supporting Information). Our analysis shows that using AOD > 0.3 as a threshold for polluted environment over both land and ocean does not significantly change the influence of aerosols on MCSs’ lifetime (Fig. S1). Also, the MCSs’ lifetime is not significantly correlated with the lower tropospheric VWS (Fig. S2). Moreover, response of MCSs’ lifetime to a 1σ increase in CAPE is weaker when we consider lower tropospheric VWS (Fig. S2 C and D), compared with that under the influence of deep VWS (Fig. 1A). Thus, we show our results in the main body of this work using AOD > 0.15 as the threshold for a polluted environment over land, and AOD > 0.3 as the threshold for a polluted environment over ocean, and also show results for deep tropospheric VWS.

18

20 40 60 80 100 RH850 [%]

24

30

Fig. 1. The rate of change of MCSs’ lifetime (hours, color shades) with (A) CAPE, (B) RH850, and (C) RH500 as a function of fAOD using AOD > 0.15 over land and AOD > 0.3 over the ocean as the threshold for a polluted environment and the deep tropospheric VWS, (D) fAOD as a function of RH850 and the deep tropospheric VWS, (E) fAOD as a function of CAPE and the deep tropospheric VWS, and (F) VWS as a function of the RH850 and fAOD. The rates represent a change of MCSs’ lifetime in hours associated with the variation of 1σ of that parameter. The bins with the number of samples less than 20 and insignificant rate of change at 95% confidence level are not shown. Note that fAOD = NPAOD > 0.15/(NPAOD > 0.15 + NPAOD < 0.15) over land and fAOD = NPAOD > 0.30/(NPAOD > 0.30+ NPAOD < 0.30) over ocean. NP is the number of pixels. The pie chart in G shows the fraction of variance of MCSs’ lifetime explained by the environmental variables using multiple linear regression for all of the MCSs.

PNAS | July 5, 2016 | vol. 113 | no. 27 | 7427

EARTH, ATMOSPHERIC, AND PLANETARY SCIENCES

Most of the studies discussed above were based on either model simulations or limited case studies from field campaigns. Observational assessments of how the influence of aerosols depends on meteorological conditions and the convective life cycle, and whether aerosols could have a significant influence on the lifetime of the MCSs at climate and large spatial scales, have not been available (11). Stevens and Feingold (16) suggested that climate models with parameterized cloud and turbulence schemes tend to overestimate the effect of aerosols on clouds compared with the models that resolve cloud processes and large eddies. The extrapolation of the results obtained from a limited number of field campaigns to global and climate scales is questionable. Hence, whether aerosols could increase the lifetime of the MCSs on large spatial and climate scales, and, if so, how such an effect would vary with meteorological conditions and the convective life cycle, has not been clear. Detection, measurement, and retrieval of the aerosols’ properties from the satellite data in the vicinity of deep convective clouds at various stages of a cloud’s lifecycle are challenging (17). Isolating the impact of aerosols from other meteorological conditions observationally requires large samples of the MCSs under similar meteorological conditions and at different phases of their life cycle (11). For this purpose, we use a suite of collocated geostationary and polar-orbital satellite measurements, along with reanalysis products, constrained by both physical principles in the atmospheric model and observations, to determine the phases of MCSs’ lifecycle, MCSs’ lifetime, and ice water content of the MCSs and associated ambient meteorological and aerosol conditions (see Methods). We examine 2,430 cases of the MCSs with collocated International Satellite Cloud Climatology Project (ISCCP) convective clouds and Moderate Resolution Imaging Spectroradiometer aerosol optical depth (MODIS AOD) measurements to estimate changes of MCSs’ lifetime with aerosol as a function of changing meteorological conditions. We estimate changes in MCSs’ lifetime with various meteorological parameters as a function of aerosol and other meteorological parameters. We can also estimate the relative influences of aerosols and various meteorological parameters on MCSs’ lifetime using multiple linear regression at different stages of the convective lifecycle and over different regions.

A

B

Lifetime_Growing

Unexplained 9% fAOD 0%

CAPE 36%

VWS 34% RH850 19%

RH500 2%

D

Unexplained 2% fAOD 4% VWS 12% RH500 0%

G

CAPE 39% RH850 43%

VWS 14%

CAPE 61% RH850 2%

H fAOD 26% VWS 0%

CAPE 9%

fAOD 24% RH850 46%

CAPE 57%

RH850 8%

CAPE 3%

RH850 64%

RH850 44% RH500 9%

VWS 9%

F

IWC_216_Mature

Unexplained 0%

Lifetime_Decaying Unexplained 5%

CAPE 25%

IZ_Mature

Unexplained 18% fAOD 3%

RH500 0%

IWC_216_Growing

Unexplained 18% fAOD 3% VWS 12%

E

IZ_Growing

C

Lifetime_Mature

Unexplained 0% fAOD 6% VWS 4% RH500 19%

IZ_Decaying

Unexplained 0% fAOD 35% VWS 12% RH500 1%

I

CAPE 51% RH850 1%

IWC_216_Decaying

Unexplained 41% fAOD 27%

CAPE 4% RH850 6% RH500 11%

Fig. 2. Fraction of the total variance of the MCS’s lifetime (A−C), IZ (D−F), and IWC_216 (G−I) explained by the environmental variables based on the multiple linear regression during the growing (A, D, and G), mature (B, E, and H), and decaying (C, F, and I) phases of the MCSs.

These above-mentioned conditions are determined by the ISCCP convective tracking, MODIS, and the Modern-Era Retrospective analysis for Research and Applications (MERRA) datasets. Our calculations on the global tropical continental scale indicate that MCSs’ lifetime increases at the highest rates with an increase of CAPE under intermediate to high VWS values and polluted conditions. MCSs’ lifetime increases about 3–30 h for an increase of CAPE by 1σ (Fig. 1B). The multiple linear regression (Fig. 1F) also shows that CAPE explains 33% of the total variance of MCSs’ lifetime, larger than any other meteorological and aerosol conditions. RH850 is the second most influential parameter on MCSs’ lifetime, explaining 27% of its total variance (Fig. 1B). MCSs’ lifetime increases by 3–27 h per 1σ increase of RH850 under intermediate to high VWS values and polluted conditions. The lifetime of the MCSs also increases with RH500 by 3–9 h for 1σ increase of RH500 under polluted conditions (Fig. 1C), explaining only 4% of the total variance of MCSs’ lifetime. MCSs’ lifetime decreases by 6–9 h for a 1σ increase of VWS in dry environments with low fAOD values but increases by 3–30 h for a 1σ increase of VWS in the high RH850 environment. Overall, VWS explains 16% of the total variance of the lifetime of MCSs. It is the third most influential meteorological condition on MCCs’ lifetime. In addition, increases in MCSs’ lifetime with other meteorological and aerosol conditions (CAPE, RH850, RH500, and fAOD) are stronger at moderate to high levels of VWS (Fig. 1 A–E), possibly because precipitation efficiency decreases below 50% when VWS is above 20 × 10−4 s−1 (27), whereas weak to moderate VWS is associated with heavy convective rainfall (22, 34). The lifetime of MCSs appears to increase with aerosols by 3–15 h with 1σ increase of fAOD, but only under high RH850 and moderate to high VWS (Fig. 1D). MCSs’ lifetime decreases with fAOD under lower RH850 by 3–6 h per 1σ increase of fAOD. MCSs’ lifetime increases with fAOD by 3–24 h under high values of CAPE and moderate to high values of VWS (Fig. 1E), and also decreases at low values of CAPE and VWS at the rate of 3–6 h due to 1σ increase of fAOD. Overall, fAOD explains less than 1% of the total variance of MCSs’ lifetime. Thus, on the scale of global tropical continents, the lifetime of the MCSs is mainly linked to meteorological conditions and dominated by the CAPE 7428 | www.pnas.org/cgi/doi/10.1073/pnas.1601935113

and lower tropospheric humidity, which is consistent with previous studies (18, 19, 22, 23, 35). These results do not vary significantly for fAOD using AOD > 0.3 instead of AOD > 0.15 as the threshold for polluted environment over land (Fig. 1 and Fig. S1). However, the aerosol influence on MCSs’ lifetime, cloud ice water content, and convective anvils may not become detectable by satellites until the mature or decay phases, as the available satellite sensors cannot effectively detect changes of cloud particle size in the lower and middle troposphere inside of convection and cloud thermodynamic effects. Thus, we separately analyze the variations of MCSs’ lifetime with ambient conditions during the growing, mature, and decay phases. Fig. 2 A–C shows that CAPE dominates the variance of MCSs’ lifetime during the growing phase, explains 36% of its total variance, but becomes less influential during the mature and decay phases, explaining 25% and 9%, respectively, of the total variance. RH850 dominates the variance of MCSs’ lifetime during the mature and decay phases, explaining 46% and 44%, respectively, of the total variance. VWS explains about 4–34% of the total variance of the durations of these phases. The dominant influence of RH and CAPE over VWS reflects mainly the large geographic variations of MCSs’ lifetime with the ambient atmospheric humidity that is needed to sustain the MCSs. In comparison, aerosols explain 6% and 24% of the total variances of the mature and decay phase longevity, respectively, and do not appear to have detectable influence on the duration of the growing phase of the MCSs. Are the empirical relationships between MCSs’ lifetime and ambient meteorological and aerosol conditions shown above physically reasonable? Based on the hypothesized underlying mechanism of aerosol influence on convective lifetime, an increase in MCSs’ lifetime is due to increase of cloud ice, especially in the form of smaller cloud particles above the freezing level within convective cores and anvils. Such effects may be more apparent and detectable by satellite sensors, as slower sedimentation of a larger amount of smaller ice particles can sustain convective anvils for a longer time. Thus, we analyze the vertically integrated cloud ice content of the MCSs (IZ; see Methods, Supporting Information, and Fig. S4) derived from CloudSat (Fig. 2 D–F) Chakraborty et al.

Chakraborty et al.

A

VWS 13% RH500 12%

C

B

Equatorial Africa

Unexplained 0% fAOD 8%

CAPE 9%

RH850 58%

Eq. South America_Dry Season Unexplained 2% fAOD 20% VWS 29%

CAPE 31% RH850 0% RH500 18%

South Asia_Oceanic Parameters Unexplained 0% fAOD 20% VWS 12% RH500 15%

D

CAPE 8%

RH850 45%

Eq. South America_Wet Season Unexplained 15% fAOD 22% VWS 2% RH500 11%

CAPE 41% RH850 9%

Fig. 3. Fraction of the total variance of MCS’s lifetime over different geographic regions: (A) equatorial Africa, (B) South Asia, (C) equatorial South America during dry season, and (D) equatorial South America during wet season. For the Asian domain, we use the environmental parameters associated with the MCSs over the ocean and the AOD threshold of 0.3 for polluted conditions in the vicinity of the MCSs.

of the unexplainable variance of the MCSs over the Indian continent remains unclear. The dominant MCSs’ type and aerosol conditions over equatorial South America are sharply different between the wet (December to April) and dry (June to September) seasons. Thus, we separately evaluate the influences of meteorological and aerosol conditions on the MCSs for the dry or wet season. CAPE, VWS, and fAOD explain ∼31%, ∼29%, and ∼20%, respectively, of the total variance of MCSs’ lifetime during the dry season (Fig. 3C), suggesting a comparable role between these three variables. Ambient RH500 also explains a significant fraction of lifetime variance, likely because the drier middle troposphere is a stronger limiting factor for the occurrence of MCSs during the dry season than in the wet season over equatorial South America (36). During the wet season, CAPE and fAOD explain ∼41% and ∼22%, respectively, of the total variance. The link between RH and MCSs’ lifetime in South America is weaker compared with other tropical continents, presumably because humidity is not as strong a limiting factor as in the drier continents (Fig. 4B). Could different meteorological conditions explain variations of the relative environmental influences on the MCSs between different phases of the convective lifecycle and the three tropical regions? Fig. 4A shows that aerosol concentrations and meteorological conditions do not change significantly between the three phases of the MCSs’ lifecycle. Thus, the decrease of the total variance explained by meteorological conditions is unlikely to be caused by their changes during the MCSs’ lifetime. The increased fraction of the total variance of MCSs’ lifetime explained by aerosols during the decay phase is consistent with numerical model simulations (15, 37) and supports the cloud lifetime enhancement by the aerosols during that phase. This can be explained by the significant fraction of the variance of ice content inside the convective core and smaller cloud particles in convective anvils explained by fAOD. These effects could reduce sedimentation of ice particles and dissipation of convective anvils. We cannot rule out an influence of the reduced background aerosol loading on lifetime (5, 38). Fig. 4B shows that equatorial Africa is the driest region and has highly variable RH in the low and middle troposphere, as shown by the means and SDs of the RH850 and RH500. The region also has relatively low VWS, but is rich in aerosols (fAOD). Thus, moisture PNAS | July 5, 2016 | vol. 113 | no. 27 | 7429

EARTH, ATMOSPHERIC, AND PLANETARY SCIENCES

and ice water content at 216 hPa (∼12 km above the sea level, IWC_216) of convective anvils from the Aura Microwave Limb Sounder (MLS) datasets based on 966 cases of the MCSs (over the limited period of June 2006 for CloudSat to June 2008 for ISCCP-convective tracking data sets) with collocated ISCCP, CloudSat, and MLS measurements in Fig. 2 D–I. Notice that IZ represents the mass of larger ice particles above the freezing level (above 5 km) mainly within the convective cores and lower portion of the anvils, whereas IWC_216 represents the mass of smaller particles near the top of the convective anvil. Hence, their relationships with meteorological and aerosol conditions can be different. Fig. 2D shows that IZ’s variance is primarily explained by RH850 (43%) and CAPE (39%) during the growing phase. VWS explains 12% of the total variance of the IZ (Fig. 2D). The result is broadly expected based on the controlling factors of convection reported in the literature (18–20, 22). The variance of the cloud ice in the convective anvil (IWC_216, Fig. 2G) during the growing phase is explained mostly by CAPE, presumably because convective detrainment in the upper troposphere or anvils generally increases with buoyancy. During the mature phase, CAPE and VWS explain 57% and 14%, respectively, of the total variance of the IZ (Fig. 2E), presumably because of their influence on convective mass flux above the freezing level. In contrast, RH850 and fAOD appear to be most influential on the ice water content of the anvils, explaining 64% and 26%, respectively, of its total variance, presumably because of their strong influence on the amount and effective size of smaller ice particles (15). During the decay phase, CAPE and fAOD dominate the variance of IZ, explaining 51% and 35%, respectively, of its total variance. The fAOD dominates the IWC_216 variation, explaining 27% of its total variance. VWS and RH500 also have significant influence on IWC_216, each explaining 11% of the total variance in the decay phase. Such increasing influences of aerosols, middle tropospheric humidity, and deep tropospheric VWS are also physically plausible, as the MCSs at decay phase become increasingly detached from the lower troposphere. The analysis of these independently measured variables in Fig. 2 suggests a stronger influence of aerosols on the duration of the MCSs’ decay phase, probably due to the stronger influence of aerosols on cloud ice water content in both convective cores and anvils during this phase. CAPE, RH850, and VWS dominate the explanation of variance of MCSs’ lifetime and ice water content of the convective cores and anvils during the growing and mature phases. Generally, aerosols have a stronger influence on convective anvils than on ice water content of convective cores and MCSs’ lifetime, presumably because of the formation of a larger number of smaller ice particles (15), which are sensitive to aerosol loading. The smaller ice particles are more likely suspended in convective anvils, and thus can influence the anvil ice water content and lifetime. Fig. 3 shows the relative influences of ambient aerosols versus meteorological conditions on MCSs’ lifetime within each continent. Fig. 3A shows that meteorological parameters explain ∼92% and fAOD explains up to 8% of the total variance of the lifetime of the MCSs over equatorial Africa. Over South Asia, only 39% of the total variance of the MCSs’ lifetime can be explained by meteorological and aerosol conditions when we consider all of the MCSs (Fig. S3). However, we found that 45% of the MCSs formed and matured over the Indian Ocean, Bay of Bengal, and Arabian Sea (Fig. 4C). Evaluating the relationship between these MCSs when they were over the ocean and the associated ambient meteorological conditions and fAOD (using AOD > 0.3 as the threshold), we find that RH850 dominates the lifetime of these MCSs, explaining about 45% of the total variance; fAOD also has an important impact, explaining ∼20% of the total variance (Fig. 3B). RH500 and VWS explain 15% and 12%, respectively, of the total variance. What causes the large fraction

A

μ

B

σ Growing Mature Decaying

100 90 80 70 60 50 40 30 20 10 0

μ

1500

σ Eq. Africa S. Asia Eq. S. America

2000 100 80

2000 1500

60 1000

1000 40 500 0 RH850[%]

Fraction of number of MCSs originating over land [%]; MCS lifetime [hr].

C

0

VWS[10-4S- 1] CAPE (J/kg) fAOD[%] RH500[%]

140

500

20

0 RH850[%]

VWS[10-4S- 1] CAPE (J/kg) fAOD[%] RH500[%]

Fraction of number of MCSs originated over land Mean lifetime over ocean Mean lifetime over land

120 100 80 60 40 20 0 Eq. Africa

S. Asia

Eq. S. America

availability could be the primary limiting factor for lifetime of the MCSs (19), and the RH850 and RH500 together explain about 70% of their total variance. The MCSs over equatorial South America are surrounded by the lowest fAOD (when we consider fAOD over the land sectors of all three regions) and VWS (Fig. 4B) and the highest amount of RH850. Such favorable meteorological (Fig. 1B) and aerosol conditions may explain the stronger connection between fAOD and MCSs’ lifetime over South America (Fig. 3 C and D) than those over equatorial Africa and South Asia. In addition, MCSs over equatorial South America are smaller in size (mean radius ∼168 km) and have lesser numbers of convective cores [mean number of convective cores (NCC) ∼3] and shorter lifetime ( 0.3 as the threshold for a polluted environment.

other words, 45% of the systems that end up decaying over the land spend much of their lifetime over the ocean. Throughout our analysis, fAOD explains high fractions of total variance of the MCSs’ lifetime when fAOD values are low, for example, over equatorial South America (with the lowest continental fAOD over all of the regions), over the Indian Ocean, Bay of Bengal, and Arabian Sea, and during the decaying stage. Many environmental parameters, such as the CAPE (18–21), RH (13, 20, 22, 23), and VWS (12, 22, 26, 28, 34), have been known to influence MCSs’ lifetime. However, to our knowledge, this is the first satellite-based global tropical continental-scale assessment of the relative roles of meteorological conditions versus aerosols in determining the variations of cloud ice and lifetime of the MCSs using collocated geostationary satellites, closely synchronized A-train polar orbital satellite, and MERRA reanalysis datasets on MCSs’ lifetime. These results show that aerosol can either increase MCSs’ lifetime under higher relative humidity and CAPE and moderate-to-high deep tropospheric VWS or decrease MCSs’ lifetime under lower relative humidity and CAPE and low VWS. On the global tropical continental scale, an increase in MCSs’ lifetime with aerosol is most apparent during the decay phase, presumably by increasing the ice water content of convective anvils and convective cores. Such a dependence of the aerosol effect on specific meteorological conditions and phases of the MCSs’ lifecycle presumably explains the difficulty of establishing aerosol effects on a global scale with variations in meteorological conditions. It also appears to broadly explain a stronger influence of aerosols on MCSs’ lifetime in humid equatorial South America during the wet season and over the Indian Ocean/Bay Bengal with lower aerosol loading, versus the relatively dry conditions with higher aerosol loading over equatorial Africa and dry season South America. Thus, the approach developed in this study can potentially provide a key to reconcile the differences of the aerosol effects on the MCSs reported in various literature. MCSs produce heavy rainfall, and thus their lifetime should be important for determining rainfall/flood and diabatic heating for atmospheric circulations. However, our understanding of aerosol impacts on MCSs’ lifetime has been very limited, in large part, due to lack of adequate long-term observations on large scales. The large samples of the MCSs observed by a suite of satellite Chakraborty et al.

study. For South Asia, only monsoon months (June−September) of every year are considered. Domains of analysis, details of the collocation technique, and related error estimation processes are given in Chakraborty et al. (41). Once a collocated MCS is observed, we compute fAOD using the MODIS AOD pixels (5, 31–33) within an area of 2° latitude/longitude in the vicinity of an MCS. We determine MCSs’ lifetime, phase of convective life cycle, and different convective properties from the ISCCP DX dataset (42), and calculate the CAPE, RH850, RH500, and VWS (28) from the MERRA data (41). We integrate the mean vertical reflectivity profile (Z) above the freezing level to represent the columnar IWC (43) or IZ. For details of the calculations and equations, please see Supporting Information. We have used multiple linear regression statistics (41) based on the linear combination of the predictors (CAPE, RH850, RH500, VWS, and fAOD) to form a set of predictive equations to estimate the variance of MCS explained lifetime (see Supporting Information).

We collocate MODIS and ISCCP datasets to identify the convections that are observed by different satellites between January 2003 and June 2008, to have sufficient numbers of the MCSs needed for the statistical analysis used in this

ACKNOWLEDGMENTS. We acknowledge the providers of the ISCCP, CloudSat, TRMM, Aura MLS, Aqua MODIS, and MERRA datasets. S.C. and R.F. were supported by NASA Aura Science Team Grant (NNX1172G) and the Office of Biological & Environmental Research within the Department of Energy, Office of Science Grant (DE-SC0011117). S.T.M. and S.C. are supported by NASA CALIPSO/CLOUDSAT Grant NNX14AO85G. The supercomputer at the University of Texas has been used to store and analyze the data.

1. Albrecht BA (1989) Aerosols, cloud microphysics, and fractional cloudiness. Science 245(4923):1227–1230. 2. Rosenfeld D, Woodley WL (2000) Deep convective clouds with sustained supercooled liquid water down to −37.5 °C. Nature 405(6785):440–442. 3. Rosenfeld D, Lensky IM (1998) Satellite-based insights into precipitation formation processes in continental and maritime convective clouds. Bull Am Meteorol Soc 79(11):2457–2476. 4. Rosenfeld D, et al. (2007) The Chisholm firestorm: Observed microstructure, precipitation and lightning activity of a pyro-cumulonimbus. Atmos Chem Phys 7(3): 645–659. 5. Rosenfeld D, et al. (2008) Flood or drought: How do aerosols affect precipitation? Science 321(5894):1309–1313. 6. Lin JC, Matsui T, Pielke RA, Kummerow C (2006) Effects of biomass-burning-derived aerosols on precipitation and clouds in the Amazon Basin: A satellite-based empirical study. J Geophys Res Atmos 111(D19):D19204. 7. Koren I, Kaufman YJ, Rosenfeld D, Remer LA, Rudich Y (2005) Aerosol invigoration and restructuring of Atlantic convective clouds. Geophys Res Lett 32(14):L14828. 8. Koren I, Remer LA, Altaratz O, Martins JV, Davidi A (2010) Aerosol-induced changes of convective cloud anvils produce strong climate warming. Atmos Chem Phys 10(10): 5001–5010. 9. Bister M, Kulmala M (2011) Anthropogenic aerosols may have increased upper tropospheric humidity in the 20th century. Atmos Chem Phys 11(9):4577–4586. 10. Lindsey DT, Fromm M (2008) Evidence of the cloud lifetime effect from wildfireinduced thunderstorms. Geophys Res Lett 35(22):L22809. 11. Rosenfeld D, et al. (2014) Global observations of aerosol-cloud-precipitation-climate interactions. Rev Geophys 52(4):750–808. 12. Fan JW, et al. (2009) Dominant role by vertical wind shear in regulating aerosol effects on deep convective clouds. J Geophys Res Atmos 114(D22):D22206. 13. Fan JW, Zhang RY, Li GH, Tao WK (2007) Effects of aerosols and relative humidity on cumulus clouds. J Geophys Res Atmos 112(D14):D14204. 14. Rosenfeld D (2006) Aerosol-cloud interactions control of Earth radiation and latent heat release budgets. Space Sci Rev 125(1-4):149–157. 15. Fan J, et al. (2013) Microphysical effects determine macrophysical response for aerosol impacts on deep convective clouds. Proc Natl Acad Sci USA 110(48):E4581–E4590. 16. Stevens B, Feingold G (2009) Untangling aerosol effects on clouds and precipitation in a buffered system. Nature 461(7264):607–613. 17. Koren I, Feingold G, Remer LA (2010) The invigoration of deep convective clouds over the Atlantic: Aerosol effect, meteorology or retrieval artifact? Atmos Chem Phys 10(18):8855–8872. 18. Mapes BE (1993) Gregarious tropical convection. J Atmos Sci 50(13):2026–2037. 19. Houze RA (2004) Mesoscale convective systems. Rev Geophys 42(4):RG4003. 20. Houze RA (1993) Cloud Dynamics (Academic, San Diego, CA). 21. Zipser EJ (2003) Some views on “hot towers” after 50 years of tropical field programs and two years of TRMM data. Cloud Systems, Hurricanes, and the Tropical Rainfall Measuring Mission (TRMM): A Tribute to Dr. Joanne Simpson, eds Tao W-K, Adler R (Am Meteorol Soc, Boston), pp 49–58. 22. Cotton WR, Anthes RA (1989) Storm and Cloud Dynamics (Academic, San Diego). 23. Langhans W, Yeo K, Romps DM (2015) Lagrangian investigation of the precipitation efficiency of convective clouds. J Atmos Sci 72(3):1045–1062.

24. Moncrieff MW (1978) Dynamical structure of 2-dimensional steady convection in constant vertical shear. Q J R Meteorol Soc 104(441):543–567. 25. Kingsmill DE, Houze RA (1999) Kinematic characteristics of air flowing into and out of precipitating convection over the west Pacific warm pool: An airborne Doppler radar survey. Q J R Meteorol Soc 125(556):1165–1207. 26. Weisman ML, Rotunno R (2004) “A theory for strong long-lived squall lines” revisited. J Atmos Sci 61(4):361–382. 27. Cotton WR, Anthes RA (1989) Storm and Cloud Dynamics (Academic, San Diego), p 560. 28. Petersen WA, Fu R, Chen MX, Blakeslee R (2006) Intraseasonal forcing of convection and lightning activity in the southern Amazon as a function of cross-equatorial flow. J Clim 19(13):3180–3196. 29. Kilroy G, Smith RK, Wissmeier U (2014) Tropical convection: The effects of ambient vertical and horizontal vorticity. Q J R Meteorol Soc 140(682):1756–1770. 30. Harrison SJ (1992) Book reviews: Fundamentals of Weather and Climate by Robin Mcilveen. Scott Geogr Mag 108(2):133. 31. Koren I, Martins JV, Remer LA, Afargan H (2008) Smoke invigoration versus inhibition of clouds over the Amazon. Science 321(5891):946–949. 32. Lee J, Kim J, Yang P, Hsu NC (2012) Improvement of aerosol optical depth retrieval from MODIS spectral reflectance over the global ocean using new aerosol models archived from RONET inversion data and tri-axial ellipsoidal dust database. Atmos Chem Phys 12(15):7087–7102. 33. Livingston JM, et al. (2014) Comparison of MODIS 3 km and 10 km resolution aerosol optical depth retrievals over land with airborne sunphotometer measurements during ARCTAS summer 2008. Atmos Chem Phys 14(4):2015–2038. 34. Maddox RA, Chappell CF, Hoxit LR (1979) Synoptic and meso-alpha scale aspects of flash flood events. Bull Am Meteorol Soc 60(2):115–123. 35. Del Genio AD, Wu JB (2010) The role of entrainment in the diurnal cycle of continental convection. J Clim 23(10):2722–2738. 36. Li WH, Fu R (2006) Influence of cold air intrusions on the wet season onset over Amazonia. J Clim 19(2):257–275. 37. Lohmann U, Feichter J (2005) Global indirect aerosol effects: A review. Atmos Chem Phys 5(3):715–737. 38. Van Den Heever SC, Cotton WR (2007) Urban aerosol impacts on downwind convective storms. J Appl Meteorol Climatol 46(6):828–850. 39. Riemann-Campe K, Fraedrich K, Lunkeit F (2009) Global climatology of Convective Available Potential Energy (CAPE) and Convective Inhibition (CIN) in ERA-40 reanalysis. Atmos Res 93(1-3):534–545. 40. Bhowmik SKR, Sen Roy S, Kundu PK (2008) Analysis of large-scale conditions associated with convection over the Indian monsoon region. Int J Climatol 28(6):797–821. 41. Chakraborty S, Fu R, Wright JS, Massie ST (2015) Relationships between convective structure and transport of aerosols to the upper troposphere deduced from satellite observations. J Geophys Res Atmos 120(13):6515–6536. 42. Machado LAT, Rossow WB, Guedes RL, Walker AW (1998) Life cycle variations of mesoscale convective systems over the Americas. Mon Weather Rev 126(6):1630–1654. 43. Matrosov SY (2009) A method to estimate vertically integrated amounts of cloud ice and liquid and mean rain rate in stratiform precipitation from radar and auxiliary data. J Appl Meteorol Climatol 48(7):1398–1410.

Methods

Chakraborty et al.

PNAS | July 5, 2016 | vol. 113 | no. 27 | 7431

EARTH, ATMOSPHERIC, AND PLANETARY SCIENCES

sensors over global tropical continents have enabled us to evaluate aerosol effects during different phases of MCSs’ lifecycle under similar meteorological conditions. In doing so, this work has advanced our capability to evaluate whether or not aerosols can increase convective lifetimes on the climate scale and to identify the favorable meteorological conditions for aerosols to affect the lifetime of the MCSs. These are fundamental questions that have motivated many studies for decades. However, we did not address the influence of different types of aerosols (such as absorbing or scattering) on MCSs’ lifetime, nor directly observe the effect of aerosols on cloud microphysics. Use of aerosol type data from in situ measurements, Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation as well as aerosol index data from the Ozone Monitoring Instrument along with ground-based measurements could be the next step to address these limitations.

Relative influence of meteorological conditions and aerosols on the lifetime of mesoscale convective systems.

Using collocated measurements from geostationary and polar-orbital satellites over tropical continents, we provide a large-scale statistical assessmen...
837KB Sizes 0 Downloads 5 Views