Eleven Winters
An open exploratory study · 2016–2026

Eleven Winters

Every winter, Delhi blames the weather, the farm fires in Punjab and Haryana, and Diwali. They arrive together, so nobody can say how much each one adds. We used the years when they didn't line up to pull them apart.

105 monitoring sites 39 cities 2,841,635 hours of PM2.5 Diwali 2026 forecast locked
The decade

Ten years of Delhi's air, drawn as filter tape

The monitors behind these numbers pull air through white glass-fibre tape and leave a grey spot for each sample. The dirtier the air, the darker the spot. Each strip below is one season, July to June, as recorded at the US Embassy in New Delhi. Every year the tape darkens from late October and stays dark into January. The marigold ticks are Diwali. White gaps are days the monitor did not report.

Delhi's decade, drawn as filter tape

Daily PM2.5 at the US Embassy in New Delhi from 2016 to 2026, one row per season. Every season darkens from late October to January.

Daily mean PM2.5 at the US Embassy, one strip per season from July to June, shaded by India's National AQI category. Marigold ticks mark Diwali.Source: US Embassy New Delhi monitor via OpenAQ; NAQI breakpoints from CPCB.
The design

Three things that change from year to year

Diwali follows the lunar calendar, so it moves across the burning season. In 2017 it fell on 19 October, before most farm fires. In 2020 it fell on 14 November, at their height. Delhi's firecracker rules changed four times in the same decade. And the timing of crop burning itself appears to have shifted. These changes let us compare like with like.

Diwali moves through the burning season, and the rules keep changing

Diwali dates from 2016 to 2026 placed on an October to November axis, with the firecracker rules for each year.

Diwali (Lakshmi Puja) dates, 2016–2026, and the firecracker rules in force in Delhi-NCR. The shaded band marks 1–15 November, commonly reported as the peak of paddy-residue burning.Source: Supreme Court, NGT and Delhi government orders as reported; panchang dates.
A · Weather

How much of the winter is weather?

For every monitor, a machine-learning model learns how PM2.5 responds to temperature, humidity, wind, rain and the height of the mixed layer of air near the ground. Re-running it hundreds of times with the weather shuffled gives the pollution each site would have had under typical weather.

Weather explains most of the rise from monsoon to winter in Kolkata, Lucknow or Mumbai: about 82% along the eastern Gangetic Plain and 81% outside it. In Punjab and Haryana, where paddy residue is burned, it explains only about 48%, and in Delhi-NCR, downwind, about 62%. Removing the weather narrows the winter gap between northern and other cities only from 45 to 33 µg/m³. Most of the north's winter problem is what is emitted, not the weather it is emitted into.

Weather explains part of the winter gap, not all of it

For each city, the winter minus monsoon PM2.5 gap before and after removing weather.

Winter (Nov–Jan) minus monsoon (Jul–Sep) PM2.5 by city, as observed (open circle) and after removing weather (filled). The number is the share of the gap that weather accounts for.Source: CPCB and US mission monitors; ERA5 via Open-Meteo; LightGBM weather normalisation.

Did the 2020 drop last?

Mostly. With the weather taken out, Delhi's PM2.5 was 14% below its 2016–2019 level during the lockdown. From mid-2020 on it stayed about 11% below. But since then it has barely moved: down about 0.5% a year.

Delhi stepped down in 2020 and has barely moved since

Delhi's monthly PM2.5 at the US Embassy, observed and weather-normalised, 2016 to 2026, with a step down in 2020.

Monthly PM2.5 at the US Embassy, observed and weather-normalised. The dashed line marks the change point found in the deseasonalised normalised series; the shaded band is the 2020 national lockdown.Source: US Embassy New Delhi; ERA5 via Open-Meteo; PELT change-point detection.
B · Diwali

Diwali night, with the weather and the farm smoke taken out

We built a baseline for every hour: what PM2.5 should have been given that night's weather, the time of year and how much farm smoke the wind was carrying from Punjab and Haryana. The model never sees Diwali, and every hour is predicted by a version of it that never saw that week. Diwali's excess is what the air held above that baseline from 6 pm to 3 am.

In the CPCB network median, Diwali night added between 127 and 392 µg/m³ in every year with data. Among 211 ordinary autumn nights scored the same way, the 99th percentile was 193.

Every Diwali night rises far above the weather-and-fire baseline

Hour by hour PM2.5 above the baseline around Diwali night, one line per year, coloured by the firecracker rules.

Hourly PM2.5 above the business-as-usual baseline (weather, calendar and wind-aware fire index) around Diwali night, CPCB Delhi network median, one line per year coloured by the firecracker rules. Shading marks 6 pm–3 am.Source: CPCB Delhi network via OpenAQ; NASA FIRMS; ERA5.

The rules made no visible difference

Mean Diwali-night excess was 192 µg/m³ in the two years without a bursting restriction, 303 with green crackers and 366 under a full ban. The two largest nights in the record came in 2021, under a full ban, and 2025, when green crackers were allowed again. The smallest came in 2017, when sales were suspended about ten days before an early Diwali. There are only a few Diwalis per regime, so this is not proof that bans cannot work, but it gives no sign that these did.

Bans did not shrink the Diwali spike

Diwali night excess PM2.5 for each year from 2016 to 2025, against the range of ordinary autumn nights.

Diwali-night excess PM2.5 (6 pm–3 am, above the weather-and-fire baseline) by year. Filled points: CPCB network median, coloured by the firecracker rules. Open circles: US Embassy, with missing hours carried forward. The band covers 90% of ordinary October–December nights in the network median.Source: CPCB and US Embassy monitors via OpenAQ; NASA FIRMS; ERA5.

The monitor that goes dark

The US Embassy monitor is the most-used reference for Delhi. On Diwali it often repeats its last reading and then stops reporting during the worst hours. Late-night hours go missing 13 times as often as on an ordinary night, so any Diwali average built on this monitor alone leaves out the peak.

The reference monitor goes dark when the air is worst

Share of hours missing at the US Embassy monitor from 5 pm to 6 am, Diwali nights against ordinary nights.

Share of hours with no valid reading at the US Embassy monitor, Diwali nights against ordinary autumn nights.Source: US Embassy New Delhi via OpenAQ.

A fingerprint of fireworks

Compared with the same hours a week earlier, Diwali evening air is finer and more sulphurous. The PM2.5 to PM10 ratio rises from 0.41 to 0.69, and SO₂ rises by 37%. The traffic gases go the other way: NO₂ and CO both fall slightly, as fewer people drive. Fine particles and sulphur up, traffic down: that is fireworks.

The timing follows local custom. Delhi peaks around midnight. Mumbai and Hyderabad peak between 8 and 11 pm. Chennai peaks on the evening before, when Tamil Nadu celebrates Deepavali.

Fine particles and sulphur rise; traffic gases fall

Paired comparison of Diwali evening against reference evenings for PM2.5 to PM10 ratio, sulphur dioxide, carbon monoxide and nitrogen dioxide.

Diwali evening (6 pm–midnight) against the same hours on the five nights a week earlier, at CPCB Delhi sites. Thin lines are site-years; the thick line joins the medians.Source: CPCB Delhi monitors via OpenAQ.

The timing follows local custom

Median excess PM2.5 by hour around Diwali in Delhi, Kolkata, Mumbai, Hyderabad and Chennai.

Median hourly excess PM2.5 above the baseline around Diwali in five cities (Delhi: CPCB network median; others: US mission monitors). Dotted lines mark the midnights before and after Diwali day as dated in Delhi.Source: CPCB Delhi network and US mission monitors via OpenAQ; ERA5.
C · Farm fires

Fewer fires counted, but more of them at night

Official progress on stubble burning is measured in satellite fire counts. In Punjab and Haryana, VIIRS detections in October and November fell 94% between 2021 and 2025. But the same satellite also passes at 1:30 am, and the share of fires it catches then rose from 0.2% to 4.7%. More fires are still burning in the middle of the night, which is what you would see if farmers now light them after the afternoon pass.

Daytime detections collapsed while night detections rose

Afternoon-pass fire detections fall sharply after 2021 while the night-pass share rises from under one percent to several percent.

VIIRS S-NPP fire detections in Punjab and Haryana, October–November, split by overpass, after removing fixed heat sources. Diamonds mark official Punjab season counts.Source: NASA FIRMS VIIRS S-NPP; official counts as reported by PTI and the Government of India.

Each counted fire now carries more smoke

Delhi's weather-adjusted PM2.5 follows upwind fires over one to four days. From 2016 to 2021, 100 units of the fire index went with about 2.2% more PM2.5. In 2025 the same 100 units went with 25% more. The 2025 estimate is uncertain (95% interval 0 to 55%), but it is the pattern you would expect if each detected fire now stands for more burning.

Each detected fire now goes with more smoke in Delhi

Percent rise in Delhi PM2.5 per 100 fire index units by season, with confidence intervals; 2025 stands well above earlier years.

Percent change in Delhi's weather-adjusted PM2.5 per 100 units of the wind-aware fire index, by season, with 95% confidence intervals (HAC standard errors).Source: US Embassy New Delhi; NASA FIRMS; ERA5.

The ground monitors can't confirm it yet

If burning moved to the evening, monitors inside Punjab and Haryana should see more evening pollution during the harvest. They don't, so far. But the open archive is missing 2023 and 2024, and 2025 had so little burning that a change in its timing would barely move a city monitor. We report this as a test that could not yet be run, not as evidence against the shift.

Ground monitors in Punjab and Haryana do not yet show an evening shift

Harvest evening ratio relative to December for monitors in Punjab and Haryana by season.

Evening-to-midday PM2.5 ratio in the harvest window divided by the same ratio in December, per monitor and season. Ticks mark the median. Values above 1 mean the evening share rises at harvest.Source: CPCB monitors in Punjab, Haryana and Chandigarh via OpenAQ.
D · Forecast

A forecast we can be wrong about

Explaining the past is easy to overdo. So before Diwali 2026 we locked a forecast of the excess for each possible set of firecracker rules, published its checksum, and will score it after 8 November.

See the locked forecast

More

Explore any of the 105 monitors in the city explorer, read how every number is made in data and methods, or read the full paper.