Profile image
By Ajay Shah (Reporter)
Contributor profile | More stories
Story Views

Last Hour:
Last 24 Hours:

The Diwali effect in Delhi air quality

Thursday, October 27, 2016 20:24
% of readers think this story is Fact. Add your two cents.

by Dhananjay Ghei, Arjun Gupta and Renuka Sane

As Diwali approaches, we have learned to worry about air quality. Over the last few years, several studies have noted the increase in pollution levels during the period of Diwali owing to increase in commercial activity and firework displays. However, as we show in our previous article, there is considerable variation in PM 2.5 levels in Delhi in terms of location/time/month:

  1. Time Effect: The effect of diwali is not uniform throughout the day and is more prevelant at particular time of the day than other times. We also need to adjust for the confounding effect of time: pollution levels are high during the night and low during the day.
  2. Location Effect: Several areas of Delhi are severly polluted throughout the time, whereas others see large variations in their pollution levels. All these reasons make it difficult to attribute the entire increase in PM2.5 on Diwali.
  3. Month Effect: The day of Diwali Festival varies in the Gregorian Calendar between the 17th October and 15th November every year. Existing pollution levels are already high when compared to the annual average. This is a confounding effect.

It is possible that the bad air that we see in Delhi at the time of Diwali is just the bad air quality in winter, and is not causally impacted upon by Diwali. In this article, we attempt to quantify the increase in the PM 2.5 levels during the Diwali period. Does Diwali have an impact upon air quality? If so, by how much?

Issues in research design
The opportunity to identify a Diwali effect comes from the fact that Diwali is a `moving holiday’ which takes place on a different day of each year. If this were not the case, it would be strongly correlated with changing climate.

Our ability to analyse these questions is greatly hampered by the lack of data. As of today, the data only runs from 1/2013 to 10/2016.

The air pollution caused by fireworks includes many contaminants. The data that we are studying covers only pm2.5.

Pollution levels on Diwali
The data used for the analysis comes from the US Consulate based in Chanakyapuri and the Central Pollution Control Board for 4 locations (R K Puram, Punjabi Bagh, Mandir Marg, Anand Vihar). The data consists of hourly PM 2.5 levels across the five locations from January 2013 to October 2016. We winsorise the data at 1% on both ends to remove the extreme tail values.

The effect of Diwali on pollution levels
We first estimate the effect of Diwali on daily data using an event study. We aggregate the hourly concentration of PM2.5, at each location, to arrive at the daily numbers. The day of the Lakshmi Puja is taken as the event day. Therefore, we get 3 events for each location. Next, we calculate the percentage change in PM2.5 concentration levels by differencing the logarithm of PM2.5 values. These are then re-indexed to show the cumulative change over a 20 day window.

Event study showing the change in PM2.5 around Diwali date (in days)

The solid line represents the average cumulative percentage change in PM2.5 values during the window, whereas the dashed line represents the confidence intervals calculated using the bootstrapped standard errors. We see that pollution levels start increasing one day before Diwali, and increase till two days after Diwali. It is also interesting to note that the increase in the pollution levels is significant during the two days after Diwali. This can be attributed to the fact that Diwali celebrations begin only on the night of Diwali, thereby leading to a significant increase the next day, as well as Diwali being celebrated over an extended period of time.

We now come at the same set of questions using a regression.

Contribution of Diwali on PM2.5: Regression analysis
Since Diwali is celebrated over a number of days we also define the following models:

  1. Diwali=t: Diwali
  2. Diwali={t-1:t+1}: 3 Days (day before Diwali, Diwali, day after Diwali)
  3. Diwali={t-1:t+2}: 4 Days (preceding day to two days after Diwali)

The model is as follows:

[ PM2.5_{it} = alpha + beta_1*Diwali_{t}+ beta_2*Diwali_{t}*l_{i} + m_t + h_t + l_i+epsilon_{it} ]

where, $i$ is location, and $t$ is time. Here, PM 2.5 is the hourly measured levels of the pollutant. The first model takes Diwali to be only the date of Diwali, second model defines the Diwali days from one day before to one day after and the third model considers Diwali from the preceding day to two days after Diwali. In addition, we have month ($m_t$), location ($l_i$), and hour ($h_t$) fixed effects. The base for the location interaction term is Anand Vihar. Robust standard errors are used for our analysis throughout.

Dependent variable:
Hourly PM2.5 Concentration
Diwali=t Diwali={t-1:t+1} Diwali={t-1:t+2}
(1) (2) (3)
Diwali -3.720 98.687 134.709
t = -0.177 t = 8.496*** t = 13.181***
Chanakyapuri*Diwali 17.270 -75.878 -87.035
t = 0.638 t = -5.100*** t = -6.692***
Mandir Marg*Diwali 73.078 -67.943 -66.844
t = 2.606*** t = -4.450*** t = -4.979***
Punjabi Bagh*Diwali 65.630 -49.033 -52.254
t = 2.374** t = -3.254*** t = -3.945***
R K Puram*Diwali 63.348 -54.228 -67.094
t = 2.291** t = -3.589*** t = -5.055***
Month FE Yes Yes Yes
Location FE Yes Yes Yes
Hour FE Yes Yes Yes
Observations 118,847 118,847 118,847
R2 0.264 0.264 0.266
Adjusted R2 0.264 0.264 0.266
F Statistic (df = 39; 118803) 1,091.020*** 1,094.274*** 1,103.673***

The first model (Column 1) shows that the baseline effect (i.e. at Anand Vihar) is not statistically different from non-Diwali days. For locations, other than Chanakyapuri, there is a differential effect on Diwali relative to Anand Vihar on Diwali. For instance, Diwali adds on an average 69.35 (73.07-3.72) µg/m3 PM2.5 particulate matter in air at Mandir Marg relative to Anand Vihar.

When we consider the second (Column 2) and third (Column 3) specifications, there is a statistically significant effect in Anand Vihar. The average particulate matter is 99 µg/m3 higher when we consider a two day Diwali, and 135 µg/m3 when we consider a three day Diwali period. While this may not seem much, given the already degraded air quality during these months, Diwali makes the pollution level reach alarming levels (>400, the monthly average in October November is around 340) which can have severe impacts on the health of people.

The Diwali effect is lower in other other locations relative to Anand Vihar. Thus, we see, that on the main day of Diwali, Anand Vihar is not too different from other days, while other locations have more pollutants relative to Anand Vihar. However, once we take into account 1-2 days after Diwali, we see that Anand Vihar is the most polluted location, and other locations have lower pollutants relative to Anand Vihar.

Very little is known, at present, about air quality and Diwali. Using the admittedly weak data resources, we have begun analysing this question here.

To the extent that these results are persuasive, they could help individuals plan strategies to avoid being in Delhi on these days. There is also a case for a Pigouvian tax on fireworks, in order to overcome the externality.

Previous work on Diwali, which helps us see other dimensions of Diwali, includes: Seasonal adjustment with Indian data: how big are the gains and how to do it by Rudrani Bhattacharya, Radhika Pandey, Ila Patnaik, Ajay Shah, and IEDs in Diwali and Toxic chemicals in Holi by Ajay Shah.


We encourage you to Share our Reports, Analyses, Breaking News and Videos. Simply Click your Favorite Social Media Button and Share.

Report abuse


Your Comments
Question   Razz  Sad   Evil  Exclaim  Smile  Redface  Biggrin  Surprised  Eek   Confused   Cool  LOL   Mad   Twisted  Rolleyes   Wink  Idea  Arrow  Neutral  Cry   Mr. Green

Top Stories
Recent Stories



Top Global

Top Alternative



Email this story
Email this story

If you really want to ban this commenter, please write down the reason:

If you really want to disable all recommended stories, click on OK button. After that, you will be redirect to your options page.