impact of diwali festival on aerosol optical properties ... · doi: 10.4209/aaqr.2017.04.0124...

Aerosol and Air Quality Research, 18: 522–532, 2018 Copyright © Taiwan Association for Aerosol Research ISSN: 1680-8584 print / 2071-1409 online doi: 10.4209/aaqr.2017.04.0124

Impact of Diwali Festival on Aerosol Optical Properties over an Urban City, Ahmedabad (India) Nisha Vaghmaria, Niyati Mevada, James Maliakal* Department of Physics of Physics, Electronics and Space Science, Gujarat University, Ahmedabad, 380009 Gujarat, India ABSTRACT

An attempt has been made to investigate changes in the characteristics of aerosol optical properties induced during Diwali/ New Year celebrations of 2012 over an urban city, Ahmedabad (India). Ground based measurements of average Aerosol Optical Depth (AOD) were carried out using a Microtops-II at Physics Department, Gujarat University, Ahmedabad during 10th to 16th Nov, 2012. AOD on the day just after Diwali is found to increase significantly; around 75% compared to pre-Diwali days for visible wavelength, remained at higher level for one more day (next day of New Year). Turbidity factor (β) also changed considerably, to a high value of 0.34 on the next day of Diwali, indicating a significant increase in aerosol loading associated with Diwali/ New Year festivities. Normal anti-correlation between Angstrom exponent (α) and turbidity parameter (β) is not seen on days following Diwali. The AOD and β values in the evening of the next day of Diwali increased to about 164% and 171% respectively compared to pre-Diwali days evenings; such an increase associated with the festivities has not been reported at any part of the country. Spectral variation of AOD also shows a significant change in the pattern during following two days of celebrations compared to pre-Diwali days. It is seen to follow power law more closely. Significant increase in AOD and change in spectral pattern suggest an increase of fine/accumulation mode particles due to Diwali/ New Year festivities. Keywords: Optical properties of aerosol; Aerosol Optical Depth (AOD); Angstrom parameters; Anthropogenic aerosol; Aerosol physics.

INTRODRCTION

Aerosols are minute particles in solid or liquid phase suspended in air. They occur over a wide range of sizes extending from 0.001 µm to about 100 µm. Physical and chemical properties of aerosols in the atmosphere depend on its source as well as mechanism of generation. Aerosols in the troposphere are produced directly by various natural processes such as winds, sea-sprays, volcanoes, forest fires etc. and by anthropogenic activities such as burning of fossil and biogenic fuels, forest clearing and activities such as construction, agricultural and transportation, etc. Another contribution to atmospheric aerosols comes from gas-to-particle conversion processes. Aerosols in the atmosphere larger than about 1 µm in size are mainly windblown dust and sea salt from sea-spray. Aerosols smaller than 1 µm are mostly formed by gas-to-particle conversion processes and also by incomplete burning processes. Bursting of fire crackers and burning of sparklers during festive celebrations *Corresponding author.

Tel.: +919727172824 E-mail address: [email protected]

inject lot of aerosols into the atmosphere locally and temporarily.

Apart from health effects and reducing visibility, aerosols play vital role in controlling radiation balance of earth-atmosphere system directly by scattering and absorbing solar radiation (Charlson et al., 1992; Andreae, 1995; Kaufman, 1998). They scatter incoming solar radiation back to space thereby enhance planetary albedo and exert a negative (cooling) radiative forcing (Twomey, 1977). Aerosols also cause negative radiative forcing indirectly through formation of clouds and thereby influencing cloud properties and their lifetime (Penner et al., 2003). On the other hand, certain aerosols such as black carbon and mineral dust cause positive radiative forcing (warming) by absorbing solar radiation (Kaufman et al., 2002; Lohmann and Feichter, 2005; IPCC, 2007). Direct and indirect effects of atmospheric aerosols on radiative forcing and cloud microphysics are strongly dependent on their size characteristics and chemical composition (Hegg et al., 1996; Haywood et al., 1997; Muller et al., 1999).

Atmospheric aerosols exhibit very large spatial as well as temporal variability. Distribution of aerosols in the atmosphere is decided basically by sources and sinks, but is also often influenced by prevailing meteorological conditions. Even though the natural aerosols play a very significant

Vaghmaria et al., Aerosol and Air Quality Research, 18: 522–532, 2018 523

role on a global scale, aerosols of anthropogenic origin are very important on a regional scale (Coakley et al., 1983; Kaufman and Fraser, 1983; Coakley and Cess, 1985; Kiehl and Breigleb, 1993; Andreae, 1995; Ramachandran et al., 2012). Recently Goel et al. (2017) found that incense smoke is a major source of particulate matter in the temple microenvironment of Kanpur city (India) and they observed PM10 mass concentrations of more than 2100 µg m–3 inside the temples. Impact of aerosols on climate is different from that of greenhouse gases, since lifetime of aerosols is much shorter and also due to the fact that the effects of aerosols are felt more on regional scale. Therefore, studies of aerosols at regional and local levels are very important.

Air over Ahmedabad is getting polluted by traffic jams, construction works, vehicle emissions, fumes from neighboring industrial areas, careless habits of burning waste etc. National Ambient Air Quality Standards of India prescribes annual average for PM10 as 60 µg m–3 (The Gazette of India, 2009); which is much less stringent compared to WHO’s air quality standards (WHO, 2006). Many parts of Ahmedabad city recorded annual average PM10 for the year 2012 much higher than the national average (GPCB, 2013); more than 80 µg m–3 over many parts of the city and more than 100 µg m–3 on certain pockets.

Various studies suggest exceptionally high level of pollution, much above the prescribed limit of National Ambient Air Quality standard levels in association with certain festivals and result in short term air pollution episodes. Impact of fireworks during festivals on short term perturbation in air pollution levels, especially in particulate matter, in relation to air quality and health problems have been studied by several researchers (Clark, 1997; Ravindra et al., 2003; Kulshrestha et al., 2004; Barman et al., 2008; Godri et al., 2010; Moreno et al., 2010; Thakur et al., 2010; Singh et al., 2010; Tiwari et al., 2012; Chatterjee et al., 2013). Recent investigation by Sahu and Kota (2017) found that PM2.5 is the dominant pollutant over New Delhi, and they suggested that mortality over New Delhi can be decreased by 6.2 and 6.5% by meeting the PM2.5 Indian standards and WHO set limits, respectively. Babu and Moorthy (2001) reported that fireworks are a source of soot in the atmosphere. Singh et al. (2003), Badarinath et al. (2009), Vyas and Saraswat (2012) and Singh et al. (2014) have shown an increase in AOD associated with Diwali celebrations over different parts of country. Increase in BC mass density following Diwali festival has been reported over different cities of the country by Babu and Moorthy (2001), Ramachandran and Rajesh (2007) and Nair et al. (2010). Singh et al. (2014) found a change in surface radiative forcing by –45% and atmospheric radiative forcing by +30% during Diwali days at Varanasi. They also observed an increase in both scattering coefficient as well as absorption coefficient. Such studies related to fireworks in other parts of the world have also been reported. An increase in particle number in accumulation mode (> 100 nm) has been witnessed during the Millennium Fireworks in Leipzig, Germany (Wehner et al., 2000). Wang et al. (2007) reported an increase of about 125% in PM10 levels over previous day due to fireworks display during Lantern Festival in Beijing.

Zhang et al. (2010) found that particle concentrations (PM10) during the peak hour of Chinese New Year’s firework event in Shanghai (China) were approximately three times higher than that on the preceding day, with a clear shift in the particles from nuclei mode (10–20 nm) and Aitken mode (20–100 nm) to accumulation mode (0.5–1.0 µm). Do et al. (2012) observed that the particles generated by the fireworks during annual Lantern Festival in Yanshui (Taiwan) possessed higher PM2.5-to-PM10 ratios than those present in the pre-fireworks display period. Investigation by Lin et al. (2016a) showed that Beehive firework display in YanShuei area of southern Taiwan during the Lantern Festival emits significant amounts of PM1 and PM2.5 and degrades short-term air quality and Lin et al. (2016b) found that hourly average concentrations of PM2.5, BC, total particle number, and ultrafine particle number during YanShui Beehive Firework Festival were 6.9, 2.3, 5.9, and 3.7 times greater than those during the same period on reference days.

Diwali is one of the important festivals celebrated enthusiastically every year during October/November all over India. Since New Year of Gujarat State falls on the next day of Diwali, the festival mood extends and intensifies all over the State. High concentration of anthropogenic aerosols along with gases such as SO2, NO2 etc. are injected into the atmosphere during festival period especially in urban areas by extensive burning of fire crackers. In this context, we investigated changes in aerosol optical properties over an urban city Ahmedabad (India) associated with Diwali/ New Year festivities of the year 2012. LOCATION AND CLIMATE OF STUDY AREA

Ahmedabad is a fast developing semi-arid urban city in Gujarat State (India) with very heavy traffic and lot of industries and thermal power plants in suburban areas. Gujarat State is bounded by Arabian Sea in the west and desert of Rajasthan in the northwest (Fig. 1). Climate of Ahmedabad is characterized by hot summer and dry winter. Average annual rainfall of Ahmedabad is 79 cm, in which, 95% occurs during monsoon season. Southwest monsoon brings moisture to the State during mid-June to mid-September and northerly winds bring mild chill during December to February.

Diwali celebrated throughout India on 13th Nov, 2012 and people of Gujarat celebrated the New Year on 14th Nov, 2012. On Diwali day (13th) evening, after lighting diya (lamps usually made from clay with cotton wick dipped in ghee or oil), fireworks started and continued to the next day early morning to welcome the New Year. INSTRMENTATION AND DATA UTILIZED

Basic measure of aerosol in atmosphere is Aerosol Optical Depth (AOD); which is a measure of depletion of solar radiation while passing through a column of atmosphere. Spectral dependence of depletion by particles in the atmosphere can be approximated by Angstrom power law (Angstrom, 1961) given below


Fig. 1. Location map of Gujarat State and Ahmedabad.

τ(λ) = βλ–α (1) where τ(λ) is the AOD at wavelength λ, β is the turbidity coefficient and α is the Angstrom exponent. On taking logarithm, the above equation can be written as lnτ(λ) = lnβ – αlnλ (2)

Angstrom exponent (α) is a measure of relative size

distribution of aerosols. In general, α is found to reduce with increasing aerosol size; large value of α indicates relatively high ratio of smaller particles to larger particles, suggesting dominance of fine/accumulation mode particles over coarse mode particles. Typical values of α range from nearly zero for Sahelian/Saharan desert dust particles, which are dominated by coarse mode particles to > 2.0 for fresh smoke particles, which are dominated by accumulation mode aerosols (Holben et al., 1991; Kaufman et al., 1992; Eck et al., 1999; Reid et al., 1999). Turbidity parameter (β) represents atmospheric aerosol content in the vertical direction, β essentially is the AOD at λ = 1 µm. Generally, β varies between 0.0 and 0.5; β < 0.1 indicates a clean atmosphere, while β > 0.1 depicts a turbid atmosphere (Khoshsima et al., 2014). AOD, Angstrom exponent (α) and turbidity coefficient (β) are generally used to represent optical properties of aerosols.

A better empirical relationship between depletion of solar radiation due to aerosols and wavelength has been obtained by a second order polynomial fit (Eck et al., 1999; O‘Neill et al., 2001; Schuster et al., 2006) of the form lnτ(λ) = a2(lnλ)

2 + a1lnλ + a0 (3) where, the coefficient a2 accounts for the curvature of ln τ verses ln λ plot. This curvature indicates aerosol particle size distribution; negative curvature indicates significant contribution of fine mode particle and positive curvature indicates dominance of coarse mode particle. According to Schuster et al. (2006), when a2 – a1 is ≥ 2 implies a clear dominance of fine mode aerosol and a2 – a1 ≤ 1 implies

clear dominance of coarse mode aerosol. Ground based measurements of AOD were carried out

using a Microtops-II at Physics Department, Gujarat University, Ahmedabad (23.04°N, 72.54°E and 55 m amsl) during 10th to 16th Nov, 2012 at one hour interval from 0800 hr IST to 1700 hr IST. Observational site has comparatively better greenery and is away from industries. AOD measurements could not be continued beyond afternoon of 16th due to cloudy conditions. Low level winds during the period of study over Ahmedabad were nearly easterly with very low speed, less than 7 km hr–1. Dew point temperatures were about 15°C for almost all days of observation. Hysplit back trajectories have been used to see any change in source region of aerosol.

Microtops-II used for the present study is a 5 channel (380, 500, 675, 936 and 1020 nm) hand-held sunphotometer, with field view of 2.5°, designed by Solar Light Company, USA. Microtops calculates AOD values at each wavelength assuming validity of Beer-Lambert’s law, using channel’s signal, extra-terrestrial constants, atmospheric pressure, time and location. Solar distance correction is applied based on latitude, longitude and solar zenith angle. Optical depth due to Rayleigh scattering has been subtracted from the total optical depth to obtain AOD. Optical depth corrections from other processes such as O3 and NO2 absorption are ignored in Microtops. The 936 nm channel data has not been used for the present study, since it is intended to be used for the estimation of water vapour. More detailed description of the instrument and calibration procedures are available in Morys et al. (2001), Porter et al. (2001) and Ichoku et al. (2002). Typical errors in AOD measurements using Microtops-II sunphotometer are quite low (Ichoku et al., 2002). According to Kedia and Ramachandran (2011) cumulative error due to various sources in AOD measurements with Microtops-II is about 2–5% only, which is quite reasonable. RESULTS AND DISCUSSION Spectral Dependence of AOD

To see how long the impact of festival on AOD persists, the


spectral variations of AOD for forenoon (0900 hr), noon and afternoon (1600 hr) during pre-Diwali days (average of four days before Diwali night, 10th–13th) and three consecutive days after Diwali (14th–16th) are presented in Table 1. Most conspicuous feature of this data is the significant increase of AOD on following day of Diwali (14th) especially at visible wavelengths at all three times; 0.62 to 0.81, 0.55 to 0.82 and 0.58 to 1.17 respectively at forenoon, noon and afternoon. On this day, by 1600 hr, AOD380 reached more than double of that on the previous day, suggesting an accumulation of aerosols in the atmosphere. On 15th, during forenoon and noon hours AOD values were more than that of 14th and by the evening of 15th, AOD droped to about 0.77; suggesting an increased downsettling of aerosols and the reduced production. The AOD values dropped to pre-Diwali condition on 16th forenoon. Another feature to be noted from this table is that for pre-Diwali days, during both forenoon and afternoon hours, spectral variation follows power law to a good extent, but during noon, AOD1020 shows higher values, departing from power law. The higher AOD values seen in the longer wavelength are produced mainly by bigger sized particles. But, the spectral variation of AOD for noon plotted in Fig. 2 shows that for two days following Diwali (14th and 15th), spectral variation follows power law more closely; increase in AOD1020 seen during noon has disappeared. On 16th, atmospheric condition became as that of pre-Diwali days with with comparatively higher value of AOD1020.

To ascertain the impact of Diwali festival on wavelength dependence of AOD, hour to hour spectral variations of AOD for pre-Diwali days and following day of Diwali (14th) are shown in Fig. 3. A significant increase in AOD on 14th in comparison with pre-Diwali days especially for shorter wavelength is clearly visible in the figure. During morning and evening period of pre-Diwali days, spectral variation follows power law to a large extent, but during noon hours, AOD1020 increased compared to AOD675. On 14th, AOD increased with time especially in afternoon for all wavelengths and reached more than double compared to

Table 1. AOD values during forenoon, noon and afternoon during pre-Diwali days and three consecutive days after Diwali (10–16, Nov. 2012).

WL in nm → 380 500 675 1020 Days↓ Forenoon 10..13 0.62 0.47 0.32 0.36 14 0.81 0.65 0.44 0.40 15 0.94 0.74 0.50 0.41 16 0.59 0.46 0.32 0.32

Noon 10..13 0.55 0.43 0.31 0.46 14 0.82 0.64 0.44 0.46 15 0.94 0.75 0.51 0.52 16 0.54 0.43 0.31 0.41

Afternoon 10..13 0.58 0.44 0.31 0.33 14 1.17 0.92 0.63 0.46 15 0.77 0.60 0.41 0.35 16 - - - -

Fig. 2. Spectral variation of AOD during noon for pre-Diwali and three consecutive days after Diwali.

Fig. 3. Hourly spectral variation of AOD for pre-Diwali and following day of Diwali.

pre-Diwali days by 1600 hrs especailly at visible wavelengths. Difference between AOD380 and AOD1020 in morning and evening hours for pre-Diwali days were nearly 0.27, but that has significantly changed on following day of Diwali as 0.5 in morning and 0.8 in evening. At the same time, increase in AOD1020 seen during noon period for pre-Diwali days are not observed on following day of Diwali. A similar pattern is observed on 15th also. This suggest a change in the size distribution pattern of aerosols in the atmosphere with an increase in smaller sized aerosols due to the fireworks for the following two days of Diwali. Zhao et al. (2014) reported that the firework displays have significant impact


on particle number concentrations in accumulation mode (100–1000 nm), especially in 200–500 nm. Vecchia et al. (2008) observed a significnat increase in particle number concentrations during the firework episode; up to 6.7 times in 1 hour for 0.5 < d < 1 µm size. Hourly Variation of AOD

Day-time hourly variation of AOD for all four wavelengths during 10th–16th Nov, 2012 is shown in Fig. 4. Hourly values of AOD for averaged 4 days before Diwali (10th–13th) presented in the figure shows that AOD values varied between 0.5 and 0.6 for 380 nm, 0.4 and 0.5 for 500 nm and 0.25 and 0.35 for 675 nm during pre-Diwali days. The AOD1020 shows a systematic variation with a peak during noon and minimum during morning and evening with a range of about 0.2.

On 14th morning, AOD values were higher compared to pre-Diwali days for all three visible wavelengths, especially for 380 nm, for which AOD increased from about 0.53 to about 0.83. This might be due to continuous firing of crackers and burning of sparklers during Diwali night (13th) and early morning of New Year day (14th). On 14th, AOD remained nearly same up to noon for visible wavelength, but thereafter AOD increased very significantly and systematically and reached to about 160% of pre-Diwali days value at 1700 hr for 380 nm. On the next day (15th) morning, AOD values dropped from previous day evening values, but values were more than that of previous day (14th) morning values. On this day also AOD remained nearly constant up to noon; thereafter it gradually decreased for all three visible wavelengths. On 16th, from morning onwards, AOD values

for all wavelengths were very close to pre-Diwali condition. Impact of fireworks on aerosols is clearly seen for two days after the festival with significant increase in AOD values especially for lower wavelengths. Changes in AOD1020 on the following day of Diwali (14th) with respect to pre-Diwali days were not very significant; evening AOD value has increased only about 0.2.

Hourly Variation of Angstrom Parameters

Hourly values of Angstrom parameters α and β during 10th to 16th Nov, 2012 have been estimated from plot of ln(AOD) versus ln(λ) for visible wavelengths. Such plots for pre-Diwali days (10th–13th) and for next day of Diwali (14th) for few typical hours are shown in Fig. 5. For pre-Diwali days, variation shows almost linear pattern at all times, but on the following day of Diwali (14th), variation indicates a non-linear nature with an obvious bulge. A similar pattern is observed on 15th also. This negative curvature seen on 14th and 15th at all times especially during morning and evening compared to pre-Diwali days indicates domination of fine/accumulation mode particles over coarse mode particles during the festive time (Eck et al., 1999; O’Neil et al., 2001; Kaskaoutis and Kambezidis, 2006).

Hourly variations of Angstrom parameters α and β for the period 10th to 16th Nov, 2012 are shown in Fig. 6. Angstrom exponent α for pre-Diwali days shows a systematic pattern, its value gradually decreased from 1.17 at 0800 hr to about 0.92 at 1300 hr, and then it increased to about 1.12 at 1700 hr. This indicates ratio of fine mode particle to coarse mode particle is high in morning and evening hours and gradually decreased to a minimum during noon period.

Fig. 4. Hourly variation of AOD for different wavelengths during 10th to 16th Nov. 2012.


Fig. 5. ln(AOD) versus ln(λ) plots before and following day of Diwali.

This dip observed in α during midday hours may be due to local convective activity, leading to a change in aerosol particle size distribution. During the festival period, Angstrom exponent α remained nearly constant about 1.1 throughout the day, suggesting that size distribution remained nearly same during the whole festival period.

Hourly variation of turbidity parameter β during pre-Diwali period shows nearly a constant value of 0.2. The β increased by about 40% in morning of next day of Diwali (14th) and from noon onwards β increased systematically with time and reached to more than 2.7 times at 1700 hr than that in pre-Diwali days. On 15th, β was even higher than that on 14th up to noon time, but in the afternoon hours it gradually decreased. On 16th morning, β reduced to that of pre-Diwali days. Present analysis of Angstrom parameters suggests that aerosol loading in the atmosphere has been significantly affected by the Diwali festival and particle size distribution has changed marginally. Slight increase in Angstrom exponent compared to pre-Diwali days especially

in noon hours indicates a slight increase in fine particle concentration in the atmosphere. This suggests the role of gas to particle conversion associated with enhanced release of gases during the fireworks, and possibility of direct generation of more fine mode soot particles. Variation of Daily Mean AOD and Angstrom Parameters

Average AOD for each day has been calculated using hourly AOD values from 0800 to 1700 hrs. Average AOD for individual days from 10th to 16th Nov, 2012 for all four wavelengths plotted in Fig. 7(a) shows a significant increase in AOD for visible wavelength from 13th to 14th and remained at high on 15th; then dropped to pre-Diwali condition on 16th. Vertical arrow in the figure indicates Diwali day (13th). On 14th, AOD increased by 76%, 79% and 72% respectively for 380 nm, 500 nm and 675 nm, with respect to mean of pre-Diwali days. This nature almost follows the pattern of fireworks. On evening of Diwali day


Fig. 6. Hourly variation of Angstrom parameters α and β during 10th to 16th Nov, 2012.

Fig. 7. (a) Day-to-day variation of average AOD (b) Spectral variation of average AOD during 10th to 16th Nov, 2012.

(13th), after lighting diya, fireworks started, and continued till evening of New Year day (14th). Increase in AOD1020 is not significant as that of visible wavelength; this suggests a significant increase in smaller sized particle concentration compared to bigger sized particle concentration. Singh et al. (2003) after analyzing aerosol characteristics during Diwali celebrations at Kanpur also observed more increase in AOD for 340 nm and 500 nm compared to 1020 nm. They reported that AOD increased from previous day of Diwali as against our present observation, and they attributed the practice of fireworks one or two days before Diwali, especially over urban areas for this increase. Singh et al. (2014) observed an enhancement of AOD500 over Varanasi during Diwali period compared to average for the period not affected by Diwali. Babu and Moorthy (2001) reported three times increase in mass concentration of black carbon (BC) associated with Deepavali (Diwali) festival at Trivandrum. Ramachandran and Rajesh (2007) found an increase in BC concentration during late night and early morning hours at Ahmedabad during Diwali time. Surface BC measurements can be compared with vertically integrated AOD values to a large extent during post-monsoon season because the boundary layer during this period is very shallow and winds are comparatively calm so that pollutants are trapped at lower levels.

Spectral variations of day-averaged AOD for the period 12th to 16th Nov, 2012 are shown in Fig. 7(b). Sharp increase in AOD at lower wavelength for the following two days of Diwali (14th and 15th) is clearly visible. For two consecutive days following Diwali, spectral variation of AOD follows the power law to good extent, which was not the case in pre-Diwali and post New Year days. AOD at lower wavelength has increased significantly compared to that at higher wavelength. This indicates that the Diwali/ New Year celebrations produced more fine/ accumulation mode particles in the atmosphere. This may be due to two factors as (i) the gases that are released in the form of SO2, NO2 etc. during fireworks might have converted into minute particles (ii) direct production of more fine mode soot particles during fireworks. Vecchi et al. (2008) and Kumar et al. (2016) also reported increased ambient fine particle concentrations in the atmosphere associated with firework displays.

The ln(day-averaged AOD) versus ln(λ) for a typical day before Diwali (12th) and next day of Diwali (14th) for visible wavelength are shown in Fig. 8. It is clear from the figures that on pre-Diwali day, variation follows almost a linear pattern, but next day of Diwali shows a negative curvature. The second degree polynomial fit (Eq. (3)) shows a2 changed from –0.216 on 12th to –0.669 on 14th and a2 – a1 increased from 1.144 on 12th to 1.333 on 14th. These results again indicate an increase in the dominance of fine/ accumulation mode particles over coarse mode particles in the atmosphere due to Diwali fireworks.

Temporal variations of Angstrom exponent α and turbidity parameter β calculated based on day-averaged AOD values for the period 10th to 16th Nov, 2012 are shown in Fig. 9. Vertical arrow in figure indicates the Diwali day (13th). Angstrom exponent changed only marginally, but turbidity


Fig. 8. ln(average AOD) versus ln(λ) plots for before Diwali and following day of Diwali.

factor increased significantly from 0.19 on Diwali day to 0.34 on next day (14th) of Diwali. On pre-Diwali days, αand β show a negative correlation as expected, but after Diwali, β also increased along with α.

To see whether wind pattern has any effect on changes

Fig. 9. Temporal variations of Angstrom parameters α and β during 10th to 16th Nov, 2012.

in aerosol characteristics, 24-hour back trajectories for 12th and 14th Nov, 2012 are presented in Fig. 10. It is clear from the figure that trajectories are similar, suggesting that changes observed on aerosol characteristics on 14th are only local origin. This confirms the impact of fireworks during Diwali festival on aerosol optical properties over Ahmedabad.

SUMMARY AND CONCLUSIONS

Impact of Diwali/New Year celebrations of the year 2012 on aerosol optical characteristics over an urban city Ahmedabad in India has been studied using Microtops measured AOD and the estimated Angstrom parameters. Measurements were conducted a few days before Diwali, days of Diwali and New Year and a few days after. Important results are summarized below. ● Increase in average AOD at visible wavelength on the

following day of Diwali compared to pre-Diwali days at

Fig. 10. 24-hour back trajectories for before Diwali and following day of Diwali.


Ahmedabad is much more compared to any other parts of the country that is reported. On the following day of Diwali, AOD reached to about 160% of pre-Diwali days at 1700 hr.

● Angstrom parameters show that aerosol loading in the atmosphere has been significantly affected by the Diwali/ New Year festival, but the particle size distribution has changed marginally, with a slight dominance of fine mode particle.

● Normal anti-correlation between Armstrong parameters α and β are not seen on two days after Diwali.

● Negative curvature of ln(AOD) versus ln(λ) plots on following two days of Diwali also suggests a dominance of fine particles over coarse mode particles; this may be due to enhanced gas to particle conversion or production of additional fine mode particles associated with fireworks.

ACKNOWLEDGEMENTS

Authors would like to thank Dr. S. Ramachandran, Professor, Physical Research Laboratory, Ahmedabad for fruitful discussions at various stages and Dr. V.M. Raval, Retired Professor, Gujarat University for going through the manuscript meticulously and suggesting suitable modifications. Authors also would like to thank all anonymous reviewers and editors for their valuable suggestions and comments. Financial support extended to Physics Department, Gujarat University by DST under FIST program and UGC-DRS under SAP program are duly acknowledged. REFERENCES Andreae, M.O. (1995). Climatic effects of changing

atmospheric aerosol levels. In World survey of climatology, Vol. 16, Henderson-Sellers, A. (Ed.), Elsevier, pp. 341–392.

Angstrom, A. (1961). Techniques of determining the turbidity of the atmosphere. Tellus 13: 214–223.

Babu, S.S. and Moorthy, K.K. (2001). Anthropogenic impact on aerosol black carbon mass concentration at a tropical coastal station: A case study. Curr. Sci. 81: 1208–1214.

Badarinath, K.V.S., Kharol, S.K., Sharma, A.R. and Prasad, V.K. (2009). Analysis of aerosol and carbon monoxide characteristics over Arabian Sea during crop residue burning period in the Indo-Gangetic Plains using multi- satellite remote sensing datasets. J. Atmos. Sol. Terr. Phys. 71: 1267–1276.

Barman, S.C., Singh, R., Negi, M.P.S. and Bhargava, S.K. (2008). Ambient air quality of Lucknow City (India) during use of fireworks on Diwali Festival. Environ. Monit. Assess. 137: 495–504.

Charlson, R.J., Schwartz, S.E., Hales, J.M., Cess, R.D., Coakley, Jr., J.A. Hansen J.E. and Hoffmann, D. (1992). Climate forcing by anthropogenic aerosols. Science 255: 423–440.

Chatterjee, A., Sarkar, C., Adak, A., Mukherjee, U., Ghosh,

S.K. and Raha, S. (2013). Ambient air quality during Diwali Festival over Kolkata – A mega-city in India. Aerosol Air Qual. Res. 13: 1133–1144.

Clark, H. (1997). Air pollution from fireworks. Atmos. Environ. 31: 2893–2894.

Coakley, J.A. and Cess, R.D. (1985). Response of the NCAR Community Climate Model to the radiative forcing by the naturally-occurring tropospheric aerosol. J. Atmos. Sci. 42: 1677–1692.

Coakley, J.A., Cess, R.D. and Yurevich, F.B. (1983). The effect of tropospheric aerosols on the Earth’s radiation budget; a parameterization for climate models. J. Atmos. Sci. 40: 116–138.

Do, T.M., Wang, C.F., Hsieh, Y.K. and Hsieh, H.F. (2012). Metals present in ambient air before and after a Firework Festival in Yanshui, Tainan, Taiwan. Aerosol Air Qual. Res. 12: 981–993.

Eck, T.F., Holben, B.N., Reid, J.S., Dubovik, O., Smirnov, A., O’Neill, N.T., Slutsker I. and Kinne, S. (1999). Wavelength dependence of the optical depth of biomass burning, urban, and desert dust aerosol. J. Geophys. Res. 104: 31333–31349.

Godri, K.J., Green, D.C., Fuller, G.W., Dall'osto, M., Beddows, D.C., Kelly, F.J., Harrison, R.M. and Mudway, I.S. (2010). Particulate oxidative burden associated with firework activity. Environ. Sci. Technol. 44: 8295–8301.

Goel, A., Wathore, R., Chakraborty, T. and Agrawal, M. (2017). Characteristics of exposure to particles due to incense burning inside temples in Kanpur, India. Aerosol Air Qual. Res. 17: 608–615.

GPCB (2013). Annual report 2012-13, Gujarat Pollution Control Board.

Haywood, J.M, Roberts, D.L., Slingo, A., Edwards, J.M. and Shine, K. (1997). General circulation model calculations of the direct radiative forcing by anthropogenic sulphate and fossil fuel soot aerosol. J. Clim. 10: 1562–1577.

Hegg, D.A, Hobbs, P.V., Gasso, S., Nance, J.D. and Rango, A.L. (1996). Aerosol measurements in the Artic relevant to direct and indirect radiative forcing. J. Geophys. Res. 101: 23349–23363.

Holben, B.N., Eck T.F. and Fraser, R.S. (1991). Temporal and spatial variability of aerosol optical depth in the Sahel Region in relation to vegetation remote sensing, Int. J. Remote Sens. 12: 1147–1163.

Ichoku, I., Levy, R., Kaufman, Y.J., Remer, L.A., Li, R.R., Martins, V.J., Holben, B.N., Abuhassan, N., Slutsker, I., Eck T.F. and Pietras C. (2002). Analysis of the performance characteristics of the five-channel Microtops II Sun Photometer for measuring aerosol optical thickness and precipitable water vapor. J. Geophys. Res. 107: 4179.

IPCC (2007). Climate change 2007: The physical science basis, Chapter 2. Cambridge University Press, Cambridge, UK, Contributions of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change.

Kaskaoutis, D.G. and Kambezidis, H.D. (2006). Investigation on the wavelength dependence of the aerosol optical depth in Athens area. Q. J. R. Meteorol. Soc. 132: 2217–


2234. Kaufman, Y. (1998). Smoke cloud and radiation experiment

in Brazil (SCAR-B). J. Geophys. Res. 103: 31783–31808. Kaufman, Y.J. and Fraser, R.S. (1983). Light extinction by

aerosols during summer air pollution. J. Clim. Appl. Meteorol. 22: 1694–1725.

Kaufman, Y.J., Setzer, A., Ward, D., Tanre, D., Holben, B.N., Menzel, P., Pereira, M.C. and Rasmussen, R. (1992). Biomass burning airborne and spaceborne experiment in the Amazonas (BASE-A). J. Geophys. Res. 97: 14581–14599.

Kaufman, Y.J., Tanre, D. and Boucher, O. (2002). A satellite view of aerosols in the climate system. Nature 419: 215–223.

Kedia, S. and Ramachandran, S. (2011). Seasonal variations in aerosol characteristics over an urban location and a remote site in western India. Atmos. Environ. 45: 2120–2128.

Khoshsima, M., Bidokhti, A.A. and Ahmadi-Givi, F. (2014). Variation of aerosol optical depth and Angstrom parameters at a suburban location in Iran during 2009-2010. J. Earth Syst. Sci. 123: 187–199.

Kiehl, J.T. and Breigleb, B.P. (1993). The radiative roles of sulfate aerosols and greenhouse gases in climate forcing. Science 260: 311–314.

Kulshrestha, U.C., Nageswara Rao, T., Azhaguvel, S., and Kulshrestha, M.J. (2004). Emissions and accumulation of metals in the atmosphere due to crackers and sparkles during Diwali festival in India. Atmos. Environ. 38: 4421–4425.

Kumar, M., Singh, R.K., Murari, V., Singh, A.K., Singh, R.S. and Banerjee, T. (2016). Fireworks induced particle pollution: A spatio-temporal analysis. Atmos. Res. 180: 78–91.

Lin, C.C., Tsai, J.H., Huang, K.L., Yeh, C.K.J., Chen, H.L., Chen, S.J., Lee, J.T. and Hsieh, Y.C. (2016a). Characteristics of respirable particulate metals emitted by a Beehive Firework Display in YanShuei area of Southern Taiwan. Aerosol Air Qual. Res. 16: 2227–2236.

Lin, C.C., Yang, L.S. and Cheng, Y.H. (2016b). Ambient PM2.5, black carbon, and particle size-resolved number concentrations and the Angstrom exponent value of aerosols during the firework display at the Lantern Festival in Southern Taiwan. Aerosol Air Qual. Res. 16: 373–387

Lohmann, U. and Feichter, J. (2005). Global indirect aerosol effects: A review. Atmos. Chem. Phys. 5: 715–737.

Moreno, T., Querol, X., Alastuey, A., Minguillon, M.C., Pey, J., Rodriguez, S., Miro, J.V., Felis, C. and Gibbons, W. (2007). Recreational atmospheric pollution episodes: Inhalable metalliferous particles from firework displays. Atmos. Environ. 41: 913–922.

Morys, M., Mims III, F.M., Hagerup, S., Anderson, S.E., Baker, A., Kia, J. and Walkup, T. (2001). Design, calibration, and performance of Microtops II handheld ozone monitor and sun photometer. J. Geophys. Res. 106: 14573–14582.

Muller, D., Wandinger, U. and Ansmann, A. (1999). Micro-physical particle parameters from extinction and

backscatter Lidar data by inversion with regularization. Appl. Opt. 38: 1981–1999.

Nair, A.V., Mohan Kumar K. and Satheesh, S.K. (2010). Measurements of aerosols black carbon at an urban site in Southern India. Aerosol and Clouds: Climate Change Perspectives, IASTA conference, pp. 463–465.

O’Neill, N.T., Dubovic, O. and Eck, T.F. (2001). Modified Angstrom exponent for the characterization of sub-micrometer aerosols. Appl. Opt. 40: 2368–2375.

Penner, J.E., Zhang, S.Y. and Chung, C.C. (2003). Soot and smoke may not warm climate. J. Geophys. Res. 108: 4657–4665.

Porter, J.N., Miller, M., Pietras, C., and Mottel, G. (2001). Ship-based sun photometer measurements using Microtops sun photometers. J. Atmos. Oceanic Technol. 18: 765–774.

Ramachandran, S. and Rajesh, T.A. (2007). Black carbon aerosols mass concentration over Ahmedabad, an urban location in western India: Comparison with urban site in Asia, Europe, Canada, and the United States. J. Geophys. Res. 112: D06211.

Ramachandran, S., Srivastava, R., Kedia, S. and Rajesh, T.A. (2012). Contribution of natural and anthropogenic aerosols to optical properties and radiative effects over an urban location. Environ. Res. Lett. 7: D034028.

Ravindra, K., Mor, S. and Kaushik, C.P. (2003). Short-term variation in air quality associated with firework events: A case study. J. Environ. Monit. 5: 260–264.

Reid, J.S., Eck., T.F., Christopher., S.A, Hobbs, P.V. and Holben, B. (1999). Use of the Ångstrom exponent to estimate the variability of optical and physical properties of aging smoke particles in Brazil. J. Geophys. Res. 104: 27473–27489.

Sahu, S.K. and Kota, S.H. (2017). Significance of PM2.5 air quality at the Indian capital. Aerosol Air Qual. Res. 17: 588–597.

Schuster G.L., Dubovik O. and Holben B.N. (2006). Angstrom exponent and bimodal aerosol size distribution. J. Geophys. Res. 111: D07207.

Singh, B.P., Srivastava, A.K., Tiwari, S., Singh, S., Singh, R.K., Bisht, D.S., Lal, D.M., Singh, A.K., Mall, R.K. and Srivastava, M.K. (2014). Radiative impact of fireworks at a tropical Indian location: A case study. Adv. Meteorol. 2014: 197072.

Singh, D.P, Gadi, R., Mandal, T.K., Dixit, C.K., Singh, K., Saud, T., Singh, N. and Gupta, P.K. (2010). Study of temporal variation in ambient air quality during Diwali festival in India, Environ. Monit. Assess. 169: 1–13.

Singh, R.P., Sagnik, D. and Holben, B. (2003). Aerosol behavior in Kanpur during Diwali festival. Curr. Sci. 84: 1302–1304.

Thakur, B., Chakraborty, S., Debsarkar, A., Chakraborty, S. and Srivastava, R.C. (2010). Air pollution from fireworks during festival of lights (Deepawali) in Howrah, India - A case study. Atmosfera 23: 347–365.

The Gazette of India (2009). No. 217, National Ambient Air Quality Standards, Central Pollution Control Board Notification, N. Delhi, 18th Nov. 2009.

Tiwari, S., Chate, D.M. and Srivastava, M.K. (2012).


Statistical evaluation of PM10 and distribution of PM1, PM2.5, and PM10 in ambient air due to extreme fireworks episodes (Deepawali festivals) in megacity Delhi. Nat. Hazard. 61: 521–531.

Twomey, S.A. (1977). The influence of pollution on the shortwave albedo of clouds. J. Atmos. Sci. 34: 1149–1152.

Vecchi, R., Bernardoni, V., Cricchio, D., D'Alessandro, A., Fermo, P., Lucarelli, F., Navad, S., Piazzalungab, A. and Vallia, G. (2008). The impact of fireworks on airborne particles. Atmos. Environ. 42: 1121–1132.

Vyas, B.M. and Saraswat, V. (2012). Studies of Atmospheric Aerosol’s Parameters during Pre-Diwali to Post–Diwali festival period over Indian Semi Arid Station i.e., Udaipur. Appl. Phys. Res. 4: 40–55.

Wang, Y., Zhuang, G., Xu, C. and An, Z. (2007). The air pollution caused by the burning of fireworks during the lantern festival in Beijing. Atmos. Environ. 41: 417–431.

Wehner, B, Wiedensohler, A., and Heintzenberg, J. (2000). Submicrometer aerosol size distributions and mass

concentration of the millennium fireworks 2000 in Leipzig, Germany. J. Aerosol Sci. 31: 1489–1493.

WHO (2006). Air quality guidelines for particulate matter, Ozone, Nitrogen Dioxide and sulpher dioxide, WHO/SDE/PHE/OEH/06.02.

Zhang, M., Wang, X., Chen, J., Cheng, T., Wang, T., Yang, X., Gong, Y., Geng, F. and Chen, C. (2010). Physical characterization of aerosol particles during the Chinese New Year’s firework events. Atmos. Environ. 44: 5191–5198.

Zhao, S., Yu, Y., Yin, D., Liu, N. and He, J. (2014). Ambient particulate pollution during Chinese Spring Festival in urban Lanzhou, Northwestern China. Atmos. Pollut. Res. 5: 335–343.

Received for review, April 24, 2017 Revised, October 5, 2017

Accepted, October 6, 2017

impact of diwali festival on aerosol optical properties ... · doi: 10.4209/aaqr.2017.04.0124...

Documents