Interactive comment on “ Impacts of aerosols on weather and regional climate over the Pearl River Delta megacity area in China ” by Y .

(1) We have added a statement that “It should be pointed out that the use of visibility as a proxy for the aerosol content has certain uncertainties. In particular, the aerosol optical properties may also be dependent on the relative humidity (RH), since hygroscopic aerosols will increase their size as RH increases (Zhang et al., 2008)”. Nevertheless, several previous studies have suggested that the aerosol amount (PM) is well correlated to visibility with a correlation coefficient of above 0.8 over the PRD area (Wu et al., 2005; Deng et al., 2008). In the present study, we have already excluded the days with precipitation in producing the statistics of the daily mean visibility. We have C12083

Currently, the direct and indirect forcing of aerosols on climate is highly uncertain, representing the largest uncertainty in climate predictions (IPCC, 2007).Substantial effort has been made to examine the aerosol-cloud interaction under different atmospheric conditions using numerical model simulations.Several mechanisms have been proposed to elaborate the role of aerosols on cloud development and precipitation under different dynamic and thermodynamic scenarios.It has been suggested that the aerosol effects on cloud properties vary with convective available potential energy (CAPE) and wind shear (Lee et al., 2008) and depend on the morphology of clouds (Lee et al., 2010).Fan et al. (2007b) indicated that relative humidity plays an important role in the regulating aerosol effects.Another recent study by Fan et al. (2009) revealed a dominant role of wind shear on aerosol-cloud interaction.In addition, cold pool produced by evaporative cooling has been demonstrated to considerably modulate the influence of aerosols on the development of convective systems (van den Heever and Cotton, 2006;Tao et al., 2007).Khain (2009) further suggested that heat and condensate mass budgets need to be evaluated when considering the effect of aerosols on precipitation.In addition, it has been suggested that the cloud droplet effective radius may increase or decrease with aerosol loading, depending on cloud dynamic and thermodynamic conditions and aerosol properties (Yuan et al., 2008) The plausible effect of urban aerosols on lightning enhancement in thunderstorms was first suggested by Westcott (1995).Orville et al. (2001) examined cloud-to-ground lightning flashes from the National Lightning Detection Network (NLDN) in Houston, Texas for a twelve-year period and found a higher lightning flash density near the urban area.Over major urban areas of South Korea, enhancement of cloud-to-ground lightning under high aerosols loading was documented on the basis of the measured lightning flashes and particulate matter (PM) (Kar et al., 2009).However, quantitative assessment of the aerosol effect on lightning is lacking, because of simplistic treatments of aerosol-cloud interaction in most the microphysical schemes in atmospheric numerical models and the complexity in representing electrification and discharge processes in thunderstorms.Several mechanisms and parameterizations of cloud electrification have been proposed and developed on the basis of available laboratory results to simulate the structure and location of charge distributions and lightning occurrence in thunderstorms (Williams et al., 1991(Williams et al., , 1994;;Altaratz et al., 2005;Mansell et al., 2009).
According to the long-term ground lightning observation and results from the Lightning Imaging Sensor (LIS) database, the Pearl River Delta (PRD) area (113 • -114.5 • E, 21.5 • -23.5 • N) with a cluster of large cities including Guangzhou, Shenzhen, Hong Kong and Macau, exhibits frequent lightning strikes and has a lightning density of 31.4 flashes km −2 yr −1 , the largest value observed in China (Ma, 2005).Statistical studies revealed that lightning in the PRD area occurs mainly in the spring and summer with dominant negative strikes, consisting of 96.7 % of the total cloud-ground lightning (Yi et al., 2006).The thermodynamic and dynamical conditions in this area are favorable for cloud electrification and lightning formation, since the region is located in the subtropical and coastal area with abundant solar heating and moisture sources.In addition, urbanization in this area forms a cluster of large cities with the urban land surface characteristics such as increased roughness, decreased moisture availability, and decreased thermal inertia, leading to a prominent urban heat island effect.From the perspective of the large-scale dynamics, frequent intrusions of frontal systems from the northwest during the transition seasons between summer and winter and tropi-Introduction

Conclusions References
Tables Figures

Back Close
Full With fast economic and social development in the southern region of China since the 1970s, elevated anthropogenic pollutants from coal and biomass burning and industry emissions have caused severe episodes of air pollution in the PRD area.The occurrence of low visibility (less than 10 km) has remained highly frequent (150 days per year) since 1980, on the basis of the long-term visibility observation (Deng et al., 2008).Aerosol optical depth (AOD) retrieved from the Moderate-resolution Imaging Spectroradiometer (MODIS) satellite shows a typical value of larger than 0.6 in the PRD region, because of the existence of a persistent haze layer, which is also referred to as the Asian Brown Cloud (Wu et al., 2005).In addition to primary emissions of natural and anthropogenic origins, photochemical oxidation of anthropogenic and biogenic inorganic and organic compounds leads to nucleation and growth of secondary aerosols, contributing to PM pollution in this region (Suh et al., 2001;Zhang et al., 2002aZhang et al., , 2002bZhang et al., , 2004b;;Wang et al., 2010).Several major field campaigns of air quality studies, such as the Program of Regional Integrated Experiments of Air Quality in the Pearl River Delta (PRIDE-PRD) 2004 and 2008 (Zhang et al., 2008), have been conducted in this region to monitor the pollution situation and characterize the chemical and meteorological conditions responsible for accumulations of gaseous and PM pollutants of primary and secondary origins.
In this study, we perform an analysis of lightning, precipitation, and visibility and numerical modeling to elucidate the relationship between air pollution and thunderstorms and to quantify the aerosol indirect effects on cloud development, precipitation, and lightning over the PRD area.Seven-year measurements of precipitation, lightning flashes, and visibility from 2000 to 2006 are analyzed in this area.To assess the effects of aerosols on cloud processes, precipitation, lightning activity, a Weather Research and Forecasting model with a two-moment bulk microphysical scheme is employed to simulate a mesoscale convective system in this area.Sensitivity experiments Introduction

Conclusions References
Tables Figures

Back Close
Full have been performed to reflect aerosol conditions characteristic of both polluted and clean cases to further reveal the physical mechanism for the precipitation and lightning enhancement under the polluted aerosol condition.

Data sources and model description
Lightning flashes are taken from a local lightning monitoring system, which has been established since 2000.This represents a province-wide lighting detection (LD) network consisting of 16 lightning sensors to provide the coverage of lightning flashes for the entire Guangdong Province, with a focus on the Guangzhou megacity area.The LD network detects positive and negative cloud-to-ground lightning with a positioning error of less than 1 km and measures the electricity current associated with each flash.
The precipitation data are taken from 2000 rain gauges in a weather monitoring network of automatic weather (WP3103) stations, which monitor meteorological elements in real-time, including wind, pressure, rainfall, humidity, and etc.In addition, composite radar reflectivity data from five S-band 10 cm Doppler radars (Chinese CINRAD/SA) are used to validate the model simulation.
Although intensive measurements of aerosol properties (number concentration, size, and distribution, etc.) have been conducted in two major air quality field campaigns, i.e.PRIDE-PRD 2004 and 2008, long-term aerosol measurements have not continuously made in the region.Instead, we consider observations of satellite aerosol optical depth (AOD) and visibility in the region as a proxy for pollution conditions.High-resolution (1km) AOD at the 550 nm channel over the PRD area has been derived from MODIS by the Hong Kong University of Science and Technology (http://envf.ust.hk/); the derivative method has been previously described by Li et al. (2005).
Atmospheric visibility is affected by many factors, such as absorbing gaseous pollutants and concentrations, distributions and chemical composition of aerosols.In a homogeneous atmosphere, the observed visibility (R) is correlated with the extinction coefficient (β) through the Koschmieder formula: R = 3. 2006).Over the PRD area, visibility measurements are based on daily observations in 4 cities, including Guangzhou, Shenzhen, Zengcheng, and Huiyang.The daily visibility value is averaged from four measurements per day.The number of low visibility days caused by light fog (relative humidity >90 %), which is excluded from the present work, is less than 3 days per month over the PRD area since 2000, as previously reported by Wu et al. (2007).Hence, the visibility results used in this study directly correlate with the aerosol loading condition in the atmosphere.
As a fundamental component in the high-resolution cloud-resolving model, twomoment bulk microphysical parameterizations have been widely utilized (Morrison et al., 2005;van den Heever et al., 2006;Lee et al., 2008), considering their better representation of the hydrometeors spectrum, the advantage in investigation of the aerosol-cloud interaction compared with the single-moment bulk scheme, and more efficient computational performance than the spectral bin scheme.Li et al. (2008) implemented a two-moment bulk microphysics scheme into a three-dimensional cloud resolving Weather Research and Forecasting (CR-WRF) Model version 2.2.1, which is employed in the present work to investigate the effects of aerosols on cloud and precipitation processes.A detailed description of the CR-WRF has been provided elsewhere (Li et al., 2008), and only the feathers pertinent to this present work are discussed here.Briefly, the number concentration and mass mixing ratio of five hydrometeors, including cloud, rain, ice, snow, and graupel, are explicitly predicted in the scheme with the assumption of gamma size distribution (Mitchell et al., 1990).Thirty-four microphysical processes are considered in the scheme, including seven autoconversion parameterization options for the warm rain processes and three types of heterogeneous ice nucleation.In contrast to treating the number concentration of cloud droplets or cloud condensation nuclei (CCN) to be constant in many other two-moment bulk microphysical schemes, three new prognostic variables, i.e. number concentration, surface area, and mass mixing ratio of aerosols, are implemented in this scheme.The process of aerosol activation into cloud droplet is modeled using the classic K öhlar theory and the cloud number concentration is directly predicted from aerosol concentration and Introduction

Conclusions References
Tables Figures

Back Close
Full ambient relative humility.Under the assumption of a log-normal size distribution of aerosols, the model approach with forty-eight bins is used to depict the aerosol spectrum in the scheme.The performances of this microphysical scheme under different conditions have been evaluated and an inter-comparison with other available microphysical schemes in the WRF V2.2.1 has been reported previously (Li et al., 2008a(Li et al., , 2009)).

Analysis of observations of precipitation, lightning flashes, and visibility
Measurements of precipitation, lightning flashes, and visibility from 2000 to 2006 are analyzed to investigate the relationship between rainfall, lightning frequency, and aerosols over the PRD area.The high-resolution annual mean AOD from MODIS provides an overview of the aerosol distribution over the Guangdong Province.Figure 1a depicts the annual mean AOD in 2005, showing that the Guangzhou megacity area has a higher aerosol loading with the AOD value of larger than 0.6, compared with the vicinity of the PRD region, consistent with the previous studies (Wu et al., 2005).The correlation of daily heavy precipitation and lightning with visibility averaged over four months (March to June) from 2000 to 2006 are displayed in a scatter plot (Fig. 2).
In the present study, we focus on heavy precipitation, which is defined as a daily rainfall amount of greater than 25 mm.To exclude the seasonal factors contributing to the variations in precipitation and lightning, such as more frequent intrusions of tropical Introduction

Conclusions References
Tables Figures

Back Close
Full cyclones and larger intensity of solar heating induced convection in the summer, we focus on the period from March to June during the seven years. Figure 2a illustrates that daily heavy rainfall is inversely correlated with visibility; the linear Pearson correlation coefficient is −0.739.Similarly, Fig. 2b exhibits a negative correlation between daily averaged lightning flashes and visibility, with a correlation coefficient of −0.506.Hence, the analysis of the seven-year measurements of precipitation, lightning flashes, and visibility from 2000 to 2006 in the PRD region indicates that the large lightning density and heavy rainfall amount in the PRD area are closely linked with atmospheric aerosol loading and local anthropogenic pollution contributes to the occurrences of extreme weather events, including lightning and heavy rainfall.

Design of numerical simulations
To further elucidate the effects of aerosols on cloud processes, precipitation, and lightning activity, simulations using the CR-WRF model with the two-moment bulk microphysical scheme are conducted.A mesoscale convective system over the PRD area associated with a cold frontal passage on 28 March 2009 is simulated in this study.
Figure 3 shows infrared images of clouds from the Multi-functional Transport Satellite (MTSAT-1R), which is operated by the Japan Meteorological Agency to monitor the meteorological conditions in East Asia and West Pacific.The development of a frontal system is reflected by infrared images of clouds in term of the brightness temperature variance of the cloud top.At 05:00 UTC on 28 March, an embryonic convective cell appeared to the west of the Guangdong province.The convective system intensified and progressed eastwards, and at 09:00 UTC a mesoscale convective system swept through the PRD urban area.The extended planetary boundary layer due to the surface roughness and heat island effect in the urban region and high aerosol loading likely enforced the vertical motion of the system, responsible for the occurrence of the low brightness temperature (<65 • C).From the atmospheric sounding profile over Guangzhou city (23.2 • N, 113.3 • E), the CAPE was estimated to be 1055 J Kg −1 at 00:00 UTC, exhibiting an unstable thermodynamic condition for thunderstorm initialization along the frontal boundary (Fig. 4).Introduction

Conclusions References
Tables Figures

Back Close
Full The CR-WRF model is configured by a two-way interactive nested domain with two meteorology input files from 18:00 UTC, 27 March to 18:00 UTC, 28 March 2009.The outer domain has the size of 900 km × 600 km with a 3 km horizontal resolution and is centered at (23.5 • N, 114 • E).The nested domain provides a finer resolution of 1 km and covers a 300 km × 300 km area, centering at (22.8 • N, 113.5 • E), as illustrated in Fig. 6.The 1 • ×1 • NCEP Final Global Analyses data is used to establish the initial conditions and boundary conditions.To accurately represent city-scale surface physical processes over the urban groups in the PRD area, the Noah land surface model (LSM) coupled with the single-layer urban canopy model (Kusaka et al., 2001) is utilized in the CR-WRF simulation.The high-resolution land-use data from the Guangdong local Geographic Information System (GIS) database is adopted in the domain initialization, considering the large impact of the urban land-use on the convection development (van den Heever and Cotton, 2006).The Yonsei University (YSU) scheme is used to parameterize the boundary layer processes (Hong et al., 2006).No cumulus cloud parameterization is involved in the simulation.
To investigate the aerosol effect on thunderstorm and lightning, we conduct sensitive experiments with two aerosol scenarios: a polluted case (P-case) and clean case (Ccase), on the basis of atmospheric measurements conducted in this region (Zhang et al., 2008).Both cases contain two aerosol types, i.e. ammonium sulfate and sea salt.The identical dynamics and thermodynamics settings in the sensitivity studies preclude the impacts on the simulations from other factors, such as large-scale dynamics, solar heating, and humidity sources.According to the K öhler equation used for aerosol activation in the model, the difference in the chemical composition of aerosols is equivalent to a corresponding change in the CCN effective radius (Khain et al., 2005).The polluted continental aerosols are assumed to mainly contain ammonium sulfate and the clean aerosols are assumed to be mainly sea salt (Senfield and Pandis, 2006).Hence, ammonium sulfate and sea salt represent dominantly the fine mode and coarse mode in the aerosol spectrum, respectively.Background ammonium sulfate over the continent in the P-case is initiated as an exponentially decreasing profile with the maximal Introduction

Conclusions References
Tables Figures

Back Close
Full number concentration of 2200 cm −3 and a mass mixing ratio of 5.9 µg Kg −1 at the surface, while the initial aerosol profile of ammonium sulfate in the C-case is smaller, with the maximum value of 220 cm −3 , as shown in Fig. 6.The mass and number concentration of ammonium sulfate over the sea are assumed to be half of the values over the continent.To mimic urban aerosol pollution, a production rate of ammonium sulfate is assumed to be 0.5 µg Kg −1 h −1 over the PRD metropolitan area (Fan et al., 2005;Zhang et al., 2008).The background sea salt is initiated with a maximum number concentration of about 110 cm −3 and a mass mixing ratio of 0.36 µg Kg −1 in both P-and C-cases.A sea salt emission scheme considering wind velocity and relative humility over the sea surface is included in the model (Li et al., 2008a).

Radar reflectivity and precipitation
To derive the radar reflectivity, both the simulated mass and number concentration of raindrops, snow and graupel are considered to calculate the six-moment of hydrometeors, which is defined as the reflectivity factor.Comparison of radar reflectivity is made between model simulations and observations from five S-band (10 cm) Doppler radars (the locations of the radars are marked in Fig. 5). Figure 7  horizontally averaged over the PRD area from the gauge measurements and from simulations for the C-case and P-case.The simulated precipitation in the P-and C-cases occurs about one-hour earlier than the gauge measurement, but the precipitation rate reaches the peak value at 09:00 UTC in both simulations and the measurement.The maximal rainfall rate in the P-case is consistent with the gauge measurement, with a value of 7.6 mm h −1 at 09:00 UTC, while the peak rainfall rate of 6.5 mm h −1 is noticeably smaller in the C-case.It is interesting to note that, although the onset of precipitation between the P-case and C-case is similar, the P-case produces a larger maximal rainfall rate (by 14.5 %) than the C-case.
The enhancement in the surface rainfall is further evaluated through comparison of the spatial distribution of precipitation. Figure 9 depicts the accumulated precipitation from 09:00 UTC to 15:00 UTC over the PRD area.Both the P-case and C-case reproduce the precipitation belt to the north part of the PRD area, while the coverage of the simulated precipitation belts extend slightly to the south compared with the gauge measurement.A comparison between the P-case (Fig. 9b) and C-case (Fig. 9c) reveals that precipitation increases in most parts of the region under the high aerosol loading.Note that the enhancement in precipitation is non-uniform over the entire domain.For instance, in the C-case the region with accumulated precipitation greater than 100 mm inside the main rain belt is slightly larger than that in the P-case.The non-uniform responses of precipitation to aerosol loading in different locations may be attributable to local inhomogeneous thermodynamic conditions (Ntelekos et al., 2009;Khain, 2009).
In terms of the total accumulated rainfall integrated through the duration of the thunderstorm, the P-case and C-case produce 30.5 mm and 26.4 mm total precipitation, respectively, over the out domain.Both values are comparable with the corresponding gauge measurement of 27.9 mm, but the P-case yields a value of 16 % greater than that in the C-case.Hence, the model simulations reveal that elevated aerosol concentrations increase precipitation associated with the mesoscale convective system over the PRD area.Introduction

Conclusions References
Tables Figures

Back Close
Full We examine the detailed precipitation response to aerosols by considering several distinct precipitation categories.The precipitation is categorized into four levels by the daily precipitation amount, i.e. light (<10 mm d −1 ), moderate (10-25 mm d −1 ), heavy (25-50 mm d −1 ) and extreme heavy (>50 mm d −1 ) rain. Figure 10 shows the percentage of the grid areas under a certain precipitation category over the entire model domain (a) and the percentage of the precipitation amount under a certain category over the total precipitation amount (b).It is evident that light rain is distributed over lesser geographic grid areas and the percentage of light rain in the total precipitation amount also decreases in the P-case, revealing a reduction in light precipitation because of elevated aerosol loading.This result is accordance to a previous study by Qian et al. (2009), which suggested that heavy pollution in China suppresses light rain on the basis of observation and numerical modeling.On the contrary, the heavy rain and extremely heavy rain are both enhanced in the total precipitation amount and in the geographic distribution for the P-case.Hence, the model simulations indicate that while elevated aerosol loading suppresses light and moderate precipitation, but enhances heavy precipitation.

Lightning flashes and lightning potential index
Lightning occurrence and frequency are often considered as an important indicator of convective intensity, since charge separation and electrification of thunderstorms require the co-existence and interaction of both supercooled liquid and ice crystals (Williams et al., 1991;Williams and Zhang, 1996).Mesoscale convective systems produce ground flashes by an average rate of 42 flashes min −1 , which consists of about one-fourth of the annual lightning strikes globally (Goodman and MacGorman, 1986;Bond et al., 2001Bond et al., , 2002)).Considering the complicated microphysical processes involved in the electrical field buildup and lightning discharge, it is difficult to exactly represent the lightning occurrence in atmospheric models.In the present work, we evaluate the lightning activity associated with the mesoscale convective system on the basis of a lightning potential index (LPI), which has been developed by Yair et al. (2009).Introduction

Conclusions References
Tables Figures

Back Close
Full The LPI parameterizes the potential of charge generation and separation in convective thunderstorm on the basis of the non-inductive graupel-ice mechanism, since there is a growing consensus that the occurrence of electrical charge transfer when graupel particles in the region of intense updraft collide with small ice crystal represents the dominant electrification mechanism in thunderstorms (Rakov and Uman, 2003).Under the assumption that the effective charge separation zone is between freezing level and −20 • C isotherms, determination of LPI mainly involves the simulated vertical wind component, temperature field, and the mass mixing ratios of hydrometeors explicitly predicted from the microphysical scheme, i.e.LPI = 1/V ∫∫∫εw 2 dxd yd z, where ε is decided by the partition of water content between liquid phase and ice phase.This method has been proved to be more accurate than other thermodynamic indices for lightning, such as the K-index and lifted index (Yair et al., 2009).
The temporal evolution of measured flashes from the local lightning detection network (the locations of the lightning sensors are marked in Fig. 5) is shown in Fig. 11a, including both positive and negative flashes. Figure 11a indicates that positive flashes dominate over the negative flashes and the occurrence of lightning flashes reaches the peak at 09:00 UTC.The time dependence LPI averaged over PRD calculated under the P-and C-cases is shown in Fig. 11b.Compared to the C-case, the evolution of LPI in the P-case exhibits a better consistency with the measured lightning flashes over PRD.The value of calculated LPI for the P-case and measured flashes are positive at 04:00 UTC, increase significantly after 07:00 UTC, and reach the peak around 09:00 UTC.In contrast, the predicted LPI for the C-case is noticeably delayed, compared to the P-case and measured lightning flashes.Over the entire thunderstorm duration LPI in the P-case is about 53 % higher than that in the C-case.
Figure 12  than 1.0 locates to the west of Guangzhou and migrates to the Guangzhou metropolitan and PRD area at 09:00 UTC.Such a feature, however, is absent from the calculated LPI for the C-case.It is also evident in Fig. 12 that there are discrepancies in the geographic distributions between the measured lightning activities and predicted LPI values to the northwest of PRD, which may be jointly explained by the limitation in the lightning detection system (which detects the cloud-to-ground flashes only), uncertainties in the model calculations, and the imperfect nature of the lightning potential parameterization.

Microphysical properties
We analyze the modeled microphysical properties under the different aerosol conditions to gain insights into the effects of aerosols on cloud development, precipitation, and lightning activity.Table 1 summarizes the domain-averaged microphysical properties for the different hydrometeors.For the P-case, the number concentration of cloud droplets is much larger, but the effective radius of cloud droplets is smaller than those for the C-case.These differences in the properties of cloud droplets are a direct reflection of the aerosol condition, since elevated aerosol loading leads to more numerous cloud droplets, but a smaller effective radius in the P-case.In addition, the number concentration of rain drops is slightly larger in the C-case, but the P-case corresponds to a much larger size for rain drops.Since a smaller effective radius of cloud droplets suppresses the collection/coalescence process to form raindrops, this likely explains the lower concentration of rain drops in the P-case.The larger size of rain drops in the P-case is indicative of the contributions from melting graupel.It is also clear from Table 1 that the concentrations of ice crystals and snow are higher in the P-case, suggesting a hindered warm rain process, but an enhanced mixed process because of elevated aerosol loading.
The temporal variations of the vertical profiles of the four hydrometeors (i.e.cloud water, rain water, ice, and graupel) is displayed in Fig. 13.For each quantity, the mass mixing ratio of hydrometeors is integrated horizontally at a given altitude.Figure 13a 23289 Introduction

Conclusions References
Tables Figures

Back Close
Full and b show shat that the amount of liquid cloud droplets is increasingly produced, when convection is initialized at 04:00 UTC in both P-and C-cases.This occurs since supersaturation occurs due to adiabatic cooling in the ascending air mass.The P-case yields considerably more liquid cloud water in the lower troposphere than the C-case; between 06:00 and 14:00 UTC the amount of liquid cloud water in the P-case sustains higher than that in the C-case.For the P-case, more aerosols are activated into cloud droplets and more water vapor condenses onto cloud droplets in the P-case partially explaining the enhanced cloud water amount.Furthermore, a smaller effective radius of cloud droplets and suppressed collection/coalescence processes lead to a less efficient sedimentation for smaller cloud droplets and a prolonged condensation process in the atmosphere.Clearly, the more abundant cloud water in the P-case resulting in a larger release of latent heat from droplet condensation, invigorating the convective development, to be discussed later.The enhanced cloud water above the freezing level (around 4 km) in the P-case is explained because of enhanced convection and smaller cloud droplets, which are easier to be lifted up.More supercooled water and enhanced convection promote the mixed-phase processes.Consequently, ice and graupel mass mixing ratios in the P-case are significantly increased compared to those in the C-case, as shown in Fig. 13e-f.Figure 13e shows that the region with ice crystals extends vertically to a higher altitude, indicating an elevated cloud top height.The maximum ice content is located above the −40 • C level, implying that homogeneous freezing of supercooled droplets represents a dominant mechanism for ice initiation.
It is interesting to note that there is relatively little difference in the rainwater amount between the C-case (Fig. 13c) and P-case (Fig. 13d).Such a behavior is rather surprising, since it is anticipated that there will be less rainwater due to suppressed collection/coalescence processes in the P-case.However, the less efficient conversion from cloud droplets to rain drops may be compensated by melting of graupel generated from the enhanced ice-phase process.As discussed above, the contribution of melting graupel to rainwater is supported by the predicted larger size of rain drops (Table 1) in the P-case.In addition, a larger size of hydrometeors implies a better chance to survive Introduction

Conclusions References
Tables Figures

Back Close
Full evaporation during sedimentation below the cloud bottom, explaining the intense rainfall period between 09:00 and 11:00 UTC at the ground level (0 km in the figure).Also, there is an insignificant difference in the rainfall onset, despite of the suppressed warm rain process in the P-case.

Convective strength
Strong updraft (>10 m s −1 ), which is characteristic of thunderstorms in the developing and mature stages, represents a dominant factor in the formation of heavy precipitation and lightning (Rakov and Uman, 2005;Williams et al., 1991).We evaluate the vertical convection strength under different aerosol conditions to investigate the aerosols effect on cloud dynamics.Figure 14 shows the time series of the column maximum updraft and downdraft.The maximum upward and downward velocities are calculated at each column and averaged over the entire domain.It is evident from Fig. 14 that the variations in the updraft and downdraft are always in phase and both velocities reach the peak at 09:00 UTC.The updraft and downdraft velocities are larger in the P-case.The cross section of the vertical velocity at 10:00 UTC (Fig. 15) reveals more details on the differences in convection development between the two cases.The region with the updraft velocities greater than 10 m s −1 (between 112.8 • E and 113.0 • E) is significantly more prominent in the P-case than in the C-case.Figure 15 also reveals that the patterns and locations of the maximal updraft velocity exhibit noticeable differences between the two cases, indicating that an alteration of the structure and development of the convection system by aerosols.The release of latent heat and its vertical distribution are of critical importance in the feedback of microphysics on the cloud dynamics (Seifert and Beheng, 2005).Figure 16 shows the vertical profile of latent heat release.In the developing and mature stages of the mesoscale convective system, the latent heat release between 2 and 8 km is noticeably stronger in the P-case, because of a more efficient and prolonged condensation process, as discussed before.However, the cooling induced by evaporation and melting of hydrometeors below 2 km is not appreciably different between the two cases, 23291 Introduction

Conclusions References
Tables Figures

Back Close
Full indicating that the cold pool produced by evaporative cooling is less affected by aerosol leading.This vertical distribution of latent heat release contributes to destabilization of the atmosphere and further enhancement of the convective strength.

Conclusions
In this work, we present an analysis of long-term measurements of precipitation, visibility, and lightning in the Pearl River Delta (PRD) region, China, with a focus on the Guangzhou megacity area.Measurements of precipitation, lightning flashes, and visibility from 2000 to 2006 in this region are analyzed to assess the impacts of aerosols on cloud and precipitation processes.The statistical analysis shows that both heavy rainfall and lightning flashes over PRD are inversely correlated with visibility, with the correlation coefficients of −0.739 and −0.506, respectively.The results suggest that a large lightning flash density and a heavy rainfall amount in the PRD area are closely linked with atmospheric aerosol loading and local anthropogenic pollution contributes to the occurrences of extreme weather events, including lightning and heavy rainfall.
To further elucidate the effects of aerosols on cloud processes, precipitation, and lightning activity, the CR-WRF model with a two-moment bulk microphysical scheme is employed to simulate a the polluted aerosol condition.The results also reveal that elevated aerosol loading suppresses light and moderate precipitation (less than 25 mm per day), but enhances heavy precipitation (greater than 25 mm per day).
The LPI is considered to evaluate the lightning activity associated with the mesoscale convective system.The calculated LPI exhibits a temporal and spatial consistence with lightning flashes recorded by the local LD network.The simulations suggest that LPI is enhanced by about 50 % under the polluted aerosol condition.
An analysis of the microphysical properties of hydrometeors indicates that elevated aerosol loading yields more numerous cloud droplets, but a smaller effective radius, leading to a suppressed collection/coalescence process, but an enhanced mixed process.The model predicts considerably more cloud water in the lower troposphere for the polluted case, but indicates relatively little difference in the rainwater under different aerosol conditions.This later behavior may be explained by the fact that a less efficient conversion from cloud droplets to rain drops is compensated by melting of graupel generated from the enhanced ice-phase process under the polluted aerosol condition.The model predicts a larger size of raindrops under the polluted aerosol condition, further supporting the contribution to rainwater from melting graupel.The maximum ice content is located above the −40 • C level, implying that homogeneous freezing of supercooled droplets represents a dominant mechanism for ice initiation.
During the developing and mature stages of the mesoscale convective system, the latent heat release in the middle troposphere is enhanced under the polluted aerosol condition, because of a more efficient and prolonged condensation process, but the cold pool relevant to evaporative cooling is largely unaffected by the aerosol loading.Such a distribution of latent heat release may destabilize the atmosphere and enhance convection.
It should be pointed out that the complex aerosol chemistry, including the effects of water-soluble organics (Yu, 2000), was not treated in the CCN activation scheme, because of the computational limitation in the present CR-WRF framework.Also, the impacts of black carbon aerosols may be important on the regional climate change Introduction

Conclusions References
Tables Figures

Back Close
Full over PRD (Yu et al., 2001).Future work is required to develop computationally effective parameterizations to account for the effects of aerosol compositions on CCN activation and the direct radiative forcing on cloud dynamics and regional climate by absorbing aerosols.In addition, further modeling experiments are needed to examine the impacts of aerosols on different types of convective system in the PRD area, such as typhoons which also contribute to heavy rainfall and lightning formation.Introduction

Conclusions References
Tables Figures

Back Close
Full  Full , have been closely linked with modification of cloud systems of diverse scales, ranging from isolated convective storms (van den Heever Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | cal cyclones (typhoons) developed over the western Pacific in the summertime foster strong vertical motion and development of intensive convective systems for heavy precipitation and lightning formation.

Figure
Figure 1b displays the annual cloud-ground (CG) flash density distribution in 2005 over the Guangdong Province.The lightning flashes exhibit a high flash density in the Guangzhou megacity area, with the highest flash density of 40 flashes km −1 yr −1 .A comparison between the lightning flashes and AOD show a similarity in the geographic distributions of the two quantities over the PRD area, with the largest values coinciding with the Guangzhou megacity area.The correlation of daily heavy precipitation and lightning with visibility averaged over four months (March to June) from 2000 to 2006 are displayed in a scatter plot (Fig.2).In the present study, we focus on heavy precipitation, which is defined as a daily rainfall amount of greater than 25 mm.To exclude the seasonal factors contributing to the Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | displays the measured and modeled distributions of the maximum radar reflectivity at two time periods.Comparison between the observation and simulation exhibits a general agreement in the distributions and developments of the maximum radar reflectivity, showing a northeastsouthwest orientated convective boundary along the frontal system.Most of the measured radar reflectivity along the frontal boundary is reproduced by the simulation.In particular, the simulation is consistent with the observation on the location of the reflectivity value of greater than 30 dBz, corresponding to the location of active convection development of the thunderstorm with heavy precipitation in the north of the PRD area between latitudes 23• N and 23.5• N.Precipitation measurements from 2000 gauge stations in the Guangdong province are compared with the WRF simulation.Figure 8 displays hourly rainfall rate Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | displays measured one-hour lightning flashes from 08:00 UTC to 09:00 UTC, along with the calculated LPI at 08:00 UTC in the two cases.The observation shows the highest frequency of lightning occurs to the west of Guangzhou City at 08:00 UTC and progressed slightly southeastward at 09:00 UTC.The P-case reproduces this measured characteristics, showing that the area of the LPI values larger Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | mesoscale convective system occurring on 28 Match 2009 in the Guangzhou megacity area.Comparison between observations and model simulations exhibits a general agreement in the distributions and developments of the maximum radar reflectivity.Most of the measured radar reflectivity along the frontal boundary is reproduced by the simulation, in particular, the location of the most active convection development of the thunderstorm.The modeled maximal rainfall rate and spatial distribution of precipitation are also in agreement with measurements from rain gauges over this region.Model sensitivity experiments reveal that elevated aerosol concentrations increase precipitation associated with the mesoscale convective system over the PRD area.The simulations suggest that the total precipitation is enhanced by about 16 % under Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Fan, J., Yuan, T., Comstock, J. M., Ghan, S., Khain, A., Leung, L. R., Li, Z., Martins, V. J., and Ovchinnikov, M.: Dominant role by vertical wind shear in regulating aerosol effects on deep convective clouds, J. Geophys.Res., 114, D22206, doi:10.1029/2009jd012352,2009.Farias, W. R. G., Pinto Jr, O., Naccarato, K. P., and Pinto, I. R. C. A.: Anomalous lightning activity over the Metropolitan Region of S ão Paulo due to urban effects, Atmos.Res., 91Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Figure 1.

Fig. 1 .
Fig. 1.(a) Annual mean AOD from the MODIS satellite in 2005 over southern China.(b) Annual mean CG flash density distribution in 2005 from the local lightning detection network in the Guangdong Province.

Fig. 2 .
Fig. 2. (a) Correlation between daily visibility and the total amount of daily heavy rainfall over PRD from 2000 to 2006.Both daily visibility and the total amount of daily heavy rainfall are averaged over the four-month period (from March to June).The total amount of daily heavy rainfall in each month is a summation of the daily rainfall amount from each rain gauge station with a daily rainfall greater than 25 mm.The line represents a linear regression through all data.(b) Correlation between daily visibility and daily lightning flash density (flashes km −2 ) from 2000 to 2006.Both daily visibility and the daily lightning flash density are averaged over the four-month period (from March to June).The line represents a linear regression through all data.

Fig. 3 .
Fig. 3. Infrared images of clouds from the MTSAT satellite, showing the development of a frontal system in term of the brightness temperature variance of the cloud top.

Fig. 4 .
Fig. 4. Atmospheric sounding in Guangzhou (113.33 • E, 23.17 • N) at 00:00 UTC on 28 March 2009.The black solid line corresponds to the temperature, and the black dash line represents the dew point temperature.

Fig. 6 .
Fig. 6.Vertical profiles of the number concentrations of ammonium sulfate in the polluted (P-) and clean (C-) cases for model initialization.

Fig. 8 .Fig. 9 .
Fig. 8. Temporal evolution of the domain-averaged rainfall rate.The red dashed line represents the gauge measurement, the blue solid line corresponds to the P-case, and the green solid line corresponds to the C-case.

Fig. 10 .
Fig. 10.Probability distribution functions of the four different rainfall categories, i.e. light (<10 mm d −1 ), moderate (10-25 mm d −1 ), heavy (25-50 mm d −1 ) and extreme heavy (>50 mm d −1 ) rain, simulated for the P-and C-cases.(a) Percentage of grid areas under a certain precipitation category over the entire model domain.(b) Percentage of the precipitation amount under a certain category over the total precipitation amount.Dark blue corresponds to the P-case, and red corresponds to the C-case.