Articles | Volume 21, issue 17
Atmos. Chem. Phys., 21, 13609–13630, 2021
Atmos. Chem. Phys., 21, 13609–13630, 2021

Research article 13 Sep 2021

Research article | 13 Sep 2021

In situ ozone production is highly sensitive to volatile organic compounds in Delhi, India

In situ ozone production is highly sensitive to volatile organic compounds in Delhi, India
Beth S. Nelson1, Gareth J. Stewart1, Will S. Drysdale1,2, Mike J. Newland1,d, Adam R. Vaughan1, Rachel E. Dunmore1, Pete M. Edwards1, Alastair C. Lewis1,2, Jacqueline F. Hamilton1, W. Joe Acton3,a, C. Nicholas Hewitt3, Leigh R. Crilley4,b, Mohammed S. Alam4, Ülkü A. Şahin5, David C. S. Beddows2,4, William J. Bloss4, Eloise Slater6, Lisa K. Whalley6,7, Dwayne E. Heard6, James M. Cash8, Ben Langford8, Eiko Nemitz8, Roberto Sommariva4, Sam Cox9, Shivani10,c, Ranu Gadi10, Bhola R. Gurjar11, James R. Hopkins1,2, Andrew R. Rickard1,2, and James D. Lee1,2 Beth S. Nelson et al.
  • 1Wolfson Atmospheric Chemistry Laboratories, Department of Chemistry, University of York, Heslington, York, YO10 5DD, UK
  • 2National Centre for Atmospheric Science, University of York, Heslington, York, YO10 5DD, UK
  • 3Lancaster Environment Centre, Lancaster University, Lancaster, LA1 4YW, UK
  • 4School of Geography, Earth and Environmental Sciences, University of Birmingham, Birmingham, B15 2TT, UK
  • 5İstanbul Üniversitesi-Cerrahpaşa Mühendislik Fakültesi Üniversite Mahallesi Bağlariçi Caddesi No:7, 34320 Avcılar, İstanbul, Turkey
  • 6School of Chemistry, University of Leeds, Leeds, LS2 9JT, UK
  • 7National Centre for Atmospheric Science, University of Leeds, Leeds, LS2 9JT, UK
  • 8UK Centre for Ecology and Hydrology, Penicuik, Midlothian, Edinburgh, EH26 0QB, UK
  • 9Research Software Engineering Team, University of Leicester, Leicester, LE1 7RH, UK
  • 10Department of Applied Sciences and Humanities, Indira Gandhi Delhi Technical University for Women, Kashmere Gate, New Delhi, Delhi, 110006, India
  • 11Indian Institute of Technology, Roorkee, Uttarakhand, 247667, India
  • anow at: School of Geography, Earth and Environmental Sciences, University of Birmingham, Birmingham, B15 2TT, UK
  • bnow at: Department of Chemistry, York University, Toronto, Ontario, M3J 1P3, Canada
  • cnow at: Department of Chemistry, Miranda House, University of Delhi, New Delhi, Delhi 110007, India
  • dnow at: ICARE-CNRS, 1 C Av. de la Recherche Scientifique, 45071 Orléans CEDEX 2, France

Correspondence: Beth S. Nelson ( and James D. Lee (


The Indian megacity of Delhi suffers from some of the poorest air quality in the world. While ambient NO2 and particulate matter (PM) concentrations have received considerable attention in the city, high ground-level ozone (O3) concentrations are an often overlooked component of pollution. O3 can lead to significant ecosystem damage and agricultural crop losses, and adversely affect human health. During October 2018, concentrations of speciated non-methane hydrocarbon volatile organic compounds (C2–C13), oxygenated volatile organic compounds (o-VOCs), NO, NO2, HONO, CO, SO2, O3, and photolysis rates, were continuously measured at an urban site in Old Delhi. These observations were used to constrain a detailed chemical box model utilising the Master Chemical Mechanism v3.3.1. VOCs and NOx (NO + NO2) were varied in the model to test their impact on local O3 production rates, P(O3), which revealed a VOC-limited chemical regime. When only NOx concentrations were reduced, a significant increase in P(O3) was observed; thus, VOC co-reduction approaches must also be considered in pollution abatement strategies. Of the VOCs examined in this work, mean morning P(O3) rates were most sensitive to monoaromatic compounds, followed by monoterpenes and alkenes, where halving their concentrations in the model led to a 15.6 %, 13.1 %, and 12.9 % reduction in P(O3), respectively. P(O3) was not sensitive to direct changes in aerosol surface area but was very sensitive to changes in photolysis rates, which may be influenced by future changes in PM concentrations. VOC and NOx concentrations were divided into emission source sectors, as described by the Emissions Database for Global Atmospheric Research (EDGAR) v5.0 Global Air Pollutant Emissions and EDGAR v4.3.2_VOC_spec inventories, allowing for the impact of individual emission sources on P(O3) to be investigated. Reducing road transport emissions only, a common strategy in air pollution abatement strategies worldwide, was found to increase P(O3), even when the source was removed in its entirety. Effective reduction in P(O3) was achieved by reducing road transport along with emissions from combustion for manufacturing and process emissions. Modelled P(O3) reduced by  20 ppb h−1 when these combined sources were halved. This study highlights the importance of reducing VOCs in parallel with NOx and PM in future pollution abatement strategies in Delhi.

1 Introduction

The majority of the world's population now lives in urban areas. This is projected to increase from 55 % in 2018 to 68 % of the global population by 2030, with 90 % of this growth occurring in Asia and Africa (Molina, 2021; United Nations, Department of Economic and Social Affairs, Population Division, 2019, 2018). Rapidly increasing industrialisation and urbanisation, coinciding with fast population growth, have led to worsening air quality in many of these densely populated regions. This is driven by the increasing emissions of nitrogen oxides (NOx= NO + NO2), largely associated with transport, and volatile organic compounds (VOCs), released from a diverse range of sources. Photochemical reactions in the atmosphere then lead to the formation of a wide range of important secondary pollutants, including ozone (O3) and secondary inorganic and organic aerosol.

Tropospheric O3 is both an air pollutant and an important greenhouse gas throughout the troposphere (Skeie et al., 2020). High levels of O3 can adversely affect vegetation, global crop yields (Avnery et al., 2011), and human health, with long-term exposure increasing the risk of death from cardiovascular and respiratory illnesses (Jerrett et al., 2009), and short-term exposure leading to the exacerbation of asthma in children (Thurston et al., 1997). O3 exposure has been linked to both acute and chronic pulmonary and cardiovascular health outcomes through both animal toxicological and human clinical studies, with one study showing statistically significant decreases in the lung function of adults on an average exposure of 70 ppbV of O3 across five 6.6 h windows (WHO, 2005, 2013; Schelegle et al., 2009; EPA, 2020; Fleming et al., 2018).

As a result of increased anthropogenic emissions, tropospheric O3 increased globally during the 20th century and has continued to rise regionally in Asia during the 21st century (Fleming et al., 2018; Royal Society, 2008; Lu et al., 2020). Background tropospheric O3 has also continued to increase (Parrish et al., 2014; Tarasick et al., 2019). Since the 1990s, surface O3 trends have varied by region (Gaudel et al., 2018; Cooper et al., 2020; Lu et al., 2020), but trends in the free troposphere have been overwhelmingly positive (Gaudel et al., 2020; Liao et al., 2021). Both satellite data and global chemical transport models have identified India and East Asia as the region with the greatest O3 increases between 1980–2016 (Ziemke et al., 2019), with the rate of change per decade between 2005–2016 more than double that of the rate between 1979–2005 (Ziemke et al., 2019).

Unlike other pollutants such as NOx and SO2, ground-level O3 is not directly emitted but is formed in the atmosphere from the photochemical processing of a range of reactive precursor species (Calvert et al., 2015). O3 reduction strategies are complicated by its non-linear relationship with its precursor species (NOx, CO, and VOCs); their reduction does not necessarily lead to a reduction in O3, and O3 production is also influenced by longer-lived gases such as CH4 and CO. The urban atmosphere can be classified as being either VOC-limited or NOx-limited. In a VOC-limited environment, O3 can be effectively reduced by reducing VOCs, whereas decreasing NOx may have a limited effect or even increase the local O3 production rate, P(O3) (Li et al., 2021). In a NOx-limited regime, reducing NOx emissions is the most efficient approach for reducing P(O3). Before implementing emission abatement strategies, it is therefore important to consider which regime prevails. Highly populated urban areas are commonly VOC-limited, meaning reducing NOx without also reducing VOC sources may potentially lead to an increase in local P(O3) (Sillman, 1999; Sillman et al., 1990) and hence to an increase in ground-level O3 concentrations.

In general, O3 formation is mediated by the reactions of peroxy radicals, RO2 and HO2, formed in the OH-initiated oxidation of VOCs (Reaction R1), with NO to produce NO2 (Reactions R2 and R4). NO2 is then rapidly photolysed back to NO, forming O(3P) (Reaction R5), which can rapidly react with O2, leading to O3 (Reaction R6). This recycling of NO to NO2 leads to net production of O3 (Calvert et al., 2015). A description of how net O3 production, P(O3), is calculated from observed and modelled concentrations can be found in Sect. 3.4.


The propensity of a particular VOC species to enhance O3 production is determined by its reactivity, structure, and ambient concentration (Jenkin et al., 2017). The sources of emissions, and their relative magnitudes, differ between cities, leading to differences in the role played by individual VOCs with respect to O3 production. Previous studies in megacities have found that a range of different VOC classes lead to ground-level O3 production. A comprehensive study of chemical processing in London, during the 2012 summer ClearfLo campaign, found biogenic and longer-chain VOCs from diesel sources to be of greatest importance to VOC–hydroxyl radical reactivity and thus in situ O3 formation (Dunmore et al., 2015; Whalley et al., 2016). Studies investigating the sensitivity of O3 production to different VOC classes in Shanghai (2006–2007) found that monoaromatic species dominated, accounting for 45 % of the total O3 production (Geng et al., 2008). Another study in Shanghai (summer 2009) identified significant contributions to O3 production from but-2-enes (Ran et al., 2011). Biogenic species, such as isoprene, were found to have the greatest impact on OH reactivity, k(OH), an indicator for O3 production, in Seoul in 2015 (Kim et al., 2016) and the Pearl River Delta in 2006 (Lou et al., 2010), whereas another study pointed to the importance of monoaromatics, followed by isoprene and anthropogenic alkenes to O3 production in Seoul in spring 2016 (Simpson et al., 2020). These studies incorporated a variety of chemical detail into their models, dependent on the breadth of their VOC measurement suite, with many studies not accounting for oxygenated or biogenic species other than isoprene (Lou et al., 2010). Zavala et al. (2020) show that the key contributors to VOC–OH reactivity in the Mexico City Metropolitan Area (MCMA) between the 1990s and 2019 have changed, with aromatic and alkene contributions decreasing with reduced VOC emissions from mobile source. Oxygenated VOCs from solvent consumption and personal care products dominate the VOC–OH reactivity in recent years, leading to sustained high O3 concentrations in the MCMA (McDonald et al., 2018; Zavala et al., 2020). Understanding which precursor species are key to O3 production in any given city allows governments to introduce measures to combat air quality problems (Molina, 2021).

Over the past two decades, advances in vehicle emission technology, along with improvements in residential heating, have led to decreased NOx emissions in industrialised and highly populated regions of the western world (Georgoulias et al., 2019). In many European cities, additional measures banning vehicle types in busy areas during weekdays, and upgrading the heavy goods vehicle (HGV) fleet, have led to further reductions in ambient NO2 (Font et al., 2019). In 2013, China introduced the Air Pollution Prevention and Control Action Plan (Zhang et al., 2019). New measures included the improvement of industrial emissions standards, the promotion of cleaner fuels to replace coal in the residential sector, and the removal of older vehicles from the roads. As a result of these new controls, NO2 emissions in Beijing have decreased by 32 % since 2012 (Krotkov et al., 2016; Liu et al., 2016; Miyazaki et al., 2017). Despite these successes, surface O3 pollution across China has continued to increase (M. Li et al., 2019; Lu et al., 2018a). The overall O3 formation potential (OFP) of VOCs has increased alongside this, despite reduction strategies leading to reduced emissions of alkanes and alkenes. This is explained by an increasing emission of VOC species with higher OFPs, such as toluene and xylenes, driven by solvent use and industrial processes (M. Li et al., 2019).

The megacity of Delhi, with an estimated population of over 28 million in 2018 (United Nations, Department of Economic and Social Affairs, Population Division, 2019), suffers from some of the world's poorest air quality (WHO, 2014). Significant efforts have been made over the past two decades to reduce the air pollution burden in Delhi, including the introduction of fuel quality standards, a new metro system to improve public transport, and an odd–even traffic number plate system (Bansal and Bandivadekar, 2013; Kumar et al., 2017). Since the late 1980s, steps to mitigate the impacts of vehicle and fuel emissions in India have been implemented. Initial interventions included switching to compressed natural gas (CNG) for auto rickshaws and buses in Delhi and other major cities. Since 2010, Bharat IV fuel quality standards have been implemented across cities in India, based on the 2003 Auto Fuel Policy. These changes have resulted in reduced annual NOx emissions nationally, compared with the projected emissions without policy implementation (Bansal and Bandivadekar, 2013).

A recent study of the implications of the SARS-CoV-2 pandemic on air quality compared pollutant levels of PM2.5, PM10, NO2, CO, and O3 in Delhi, and across other Indian cities, before and after a national lockdown (Jain and Sharma, 2020; Mahato et al., 2020; Sharma et al., 2020). The results showed significant improvements in post-lockdown air quality, with reductions in ambient NO2, PM2.5, and PM10 exceeding 50 % compared to business as usual. However, increased concentrations of ground-level O3 (>10 %) were also observed and were attributed to reductions of NO leading to reduced consumption of O3. Another study found, after removing meteorological biases, concentrations of NO2 and PM2.5 at urban background sites in Delhi to have reduced by  51 % and  5 % respectively, with O3 concentrations increasing by  8 % (Shi et al., 2021). Prior to the SARS-CoV-2 pandemic, studies have observed high concentrations of ambient O3 in New Delhi, with high levels associated with anti-cyclonic conditions (sunny and warm, stagnant winds and lower humidity) typical during October, the month in which the observational dataset used in this study was obtained (Jain et al., 2005; Sharma et al., 2016). Poor air quality is exacerbated in late October–November, when additional O3 precursor pollutants are emitted from regional agricultural biomass burning of crop residues within the wider area, and from firecrackers and the burning of effigies as part of seasonal festivities (Jain et al., 2014; Sawlani et al., 2019).

Several recent studies have examined the sources of VOCs in Delhi. Top-down and bottom-up approaches have shown gas-phase organic air pollution to be predominantly from petrol (gasoline) and diesel fuel sources. In two 2018 studies at two different urban sites in Delhi, Stewart et al. (2021a) and Wang et al. (2020) showed that 52 % and 57 % of the measured mixing ratio were from combined petrol and diesel sources. Smaller contributions to the overall VOC burden were found from solid fuel combustion, 16 % (Stewart et al., 2021a) and 27 % (Wang et al., 2020). These data were in line with Guttikunda and Calori (2013), who produced an inventory which showed that petrol and diesel sources were responsible for 65 % by mass of hydrocarbons in Delhi. VOC emissions have also been shown to be dominated by petrol and diesel sources ( 50 %) in a study which conducted positive matrix factorisation (PMF) analysis on proton transfer reaction mass spectrometer flux measurements taken in a follow-on campaign in early November 2018 at the Indira Gandhi Delhi Technical University for Women (IGDTUW) (Cash et al., 2021).

There have been several recent studies focused on understanding VOC emissions from sources specific to Delhi. Stewart et al. (2021b) studied the types of intermediate-volatility and semi-volatile gases released from solid fuels in Delhi to better understand their potential impact on air quality, and Stewart et al. (2021c) produced highly speciated non-methane VOC (NMVOC) emission factors from a range of solid fuel combustion sources characteristic to Delhi, which were used by Stewart et al. (2021d) in regional policy models and global chemical transport models. Stewart et al. (2021e) found that fuel wood, crop residue, cow dung cake, and municipal solid waste burning were shown to be 30, 90, 120, and 230 times more reactive with the OH radical, which can lead to O3 formation, than liquified petroleum gas (LPG), and may be one the factors for the high O3 levels, and overall poor air quality, observed.

To implement successful ground-level O3 reduction strategies, a good understanding of the non-linear, chemically complex processing of its precursor species is imperative. This paper presents a comprehensive in situ O3 production sensitivity study, using an extensive speciated VOC and oxygenated VOC (o-VOC) measurement suite obtained in Delhi during the Atmospheric Pollution and Human Health programme in an Indian Megacity (APHH-India) DelhiFlux project post-monsoon field campaign in October 2018. Measurements of NOx, CO, SO2, O3, HONO, 34 photolysis rates, pressure, temperature, and relative humidity complete the dataset. The observations were used to constrain a detailed zero-dimensional chemical box model, incorporating the Master Chemical Mechanism v3.1.1 (MCM:, last access: March 2020). Significant VOC classes and sources contributing to the OH reactivity (k(OH)) are identified by constraining the model to the full VOC suite. By modifying the constraints of these species, their contribution to in situ O3 production was investigated by comparing relative changes in the modelled rate of O3 production. The sensitivity of in situ O3 production to changes in NO, photolysis, and particulate matter (PM) was also investigated.

2 Experimental

2.1 Site description

The APHH-India DelhiFlux post-monsoon measurement campaign took place between 4 October and 5 November 2018. The primary measurement site was located on the campus of the IGDTUW, near Kashmere Gate in Old Delhi (28.67 N, 77.23E). The campus is an open area, with some green spaces, and is close to major roads and highways: 1.5 km north of the busy Chandni Chowk market, 0.6 km north of Old Delhi (Delhi Junction) Railway Station, and 0.3 km west of National Highway 44 (Fig. 1). The campus is mostly pedestrianised, with occasional traffic activity from delivery cars, taxis, and auto rickshaws. Inter State Bus Terminal (ISBT) is located <100 m away from the measurement site. IGDTUW facilitated the sampling of ambient air from a height of  5 m and measurements were made of a large range of VOCs, o-VOCs, NOx, CO, SO2, HONO, photolysis rates, and PM.

Figure 1The measurement site (IGDTUW) in Old Delhi, north of New Delhi (left). In closer detail (right), the site was just north of the Delhi Junction Railway Station, and west of National Highway 44. Map data © Google Maps 2020.

2.2 VOC and o-VOC measurements

Two gas chromatography (GC) instruments and a proton transfer reaction time-of-flight mass spectrometer, with quadrupole ion guide (PTR-QiTOF, Ionicon Analytik, Innsbruck) were deployed to measure an extensive range of VOCs and o-VOCs. GC instrumentation included a dual-channel gas chromatography flame ionisation detector, (DC)-GC-FID, measuring C2–C8 hydrocarbons and some o-VOCs, and a two-dimensional GC flame ionisation detector (GC × GC-FID) to measure larger hydrocarbon species (C6–C13). PTR-QiTOF measurements of a variety of o-VOCs completed the VOC measurement suite of 15 alkanes, 10 alkenes, 2 alkynes, 29 aromatics, 11 carbonyls, 2 alcohols, and 15 monoterpenes (see Table S1 in the Supplement). All three instruments operated between 11 and 27 October 2018. The two GC instruments shared an inlet, located approximately 5 m above ground level. The sample line from the inlet to the laboratory was made from 1/2 in. O.D. (9.53 mm I.D.) perfluoroalkoxy (PFA) and was heated to limit adsorption of compounds to surfaces. VOCs and o-VOCs were calibrated using 4 ppbV and a 4 ppmV (diluted before calibration) gas standard cylinders (National Physical Laboratory, UK) respectively, containing a variety of VOCs and o-VOCs. The linearity of the detector response to higher mixing ratios was tested on both instruments prior to the campaigns, as described by Stewart et al. (2021a).

The (DC)-GC-FID measurement period was 5 October–27 October 2018. A 500 mL sample pre-purge (flow rate of 100 mL min−1 for 15 min) preceded a 500 mL sample collection (flow rate of 25 mL min−1 for 20 min) using a Markes International CIA Advantage. The sample was passed through a cold glass finger (−30C) to remove water, before being adsorbed onto a Markes International Ozone Precursor dual-bed sorbent cold trap using a Markes International Unity 2 for pre-concentration. After sampling, the trap was heated (250 C for 3 min) to allow for thermal desorption of the sample and passed to the GC oven in helium carrier gas. The sample was split 50 : 50 and injected into two columns in the oven, allowing for the detection of both non-oxygenated (10 m × 0.53 mm LOWOX column) and oxygenated (50 m × 0.53 mm Al2O3 PLOT column) VOCs. The oven was held at 40 C for 5 min, then ramped up to 110 C at a rate of 13 C min−1, before a final ramp to 200 C at 7 C min−1, where it was held for 30 min (Hopkins et al., 2003).

The GC × GC-FID measurement period was 11 October–4 November 2018. A 2.1 L sample (flow rate 70 mL min−1 for 30 min) was collected and passed through a cold glass finger (−30C) to remove water. The sample was absorbed onto a TO-15/TO-17 air toxic cold sorbent trap in a Markes International Unity 2 for pre-concentration. The trap was heated (250 C for 5 min) to allow for thermal desorption and the sample injected down a transfer line. It was then refocused with liquid CO2 at the head of a non-polar BPX5 column (SGE Analytical, 15 m × 0.15 µm × 0.25 mm) held at 50 psi for 60 s. This was connected to a polar BPX50 column (SGE Analytical 2 m × 0.25 µm × 0.25 mm) held at 23 psi via a modulator held at 180 C (5 s modulation, Analytical Flow Products MDVG-HT). The oven was held for 2 min at 35 C, then ramped up to 130 C at a rate of 2.5 C min−1 and held for 1 min before a final ramp to 180 C at 10 C min−1 where it was held for 8 min (Dunmore et al., 2015; Stewart et al., 2021a).

The PTR-QiTOF-MS operated from 4 October–4 November 2018 at a flow rate of 20 L min−1, subsampling from the same inlet line as the GCs. The drift tube pressure and temperature were 3.5 mbar and 60 C, respectively, giving an E/N (the ratio between electric field strength, E, and buffer gas density, N, in the drift tube) of 120 Td. The PTR-QiTOF-MS subsampled from a 1/2 in. PFA inlet line positioned 5 m above ground next to the inlet used by the GC instruments. The PTR-QiTOF-MS was calibrated daily using a calibration gas (Apel-Riemer Environmental Inc., Miami) containing 18 compounds: methanol; acetonitrile; acetaldehyde; acetone; dimethyl sulfide; isoprene; methacrolein; methyl vinyl ketone; 2-butanol; benzene; toluene; 2-hexanone; m-xylene; heptanal; α-pinene; 3-octanone and 3-octanol at 1000 ppbV (±5 %); and β-caryophyllene at 500 ppbV (±5 %), diluted dynamically into zero air. More details on the PTR-QiTOF-MS can be found in Jordan et al. (2009).

2.3 Measurements of NOx, CO, SO2, O3, and HONO

Measurements of nitrogen oxides (NOx) were made using a dual-channel chemiluminescence instrument (Air Quality Designs Inc., Colorado, USA), and carbon monoxide (CO) was measured with a resonance fluorescent instrument (Model Al5002, Aerolaser GmbH, Germany). Both instruments were calibrated every 2–3 d throughout the campaign using standard gas cylinders from the National Physical Laboratories, UK. O3 was measured using an ozone analyser (49i, Thermo Scientific) with a limit of detection of 1 ppbV. The instrument setup and calibration methodology are as described by Squires et al. (2020). An SO2 analyser (43i, Thermo Scientific) with a limit of detection of 2 ppbV was used to provide measurements of SO2. HONO was measured with a commercially available long-path absorption photometer (LOPAP) described in Heland et al. (2001) according to the standard procedures outlined in Kleffmann and Wiesen (2008), with baseline measurements taken at regular intervals (8 h). In this study, HONO was also recorded using the PTR-QiTOF-MS with the protonated HONO ion observed at m/z 48.007. This signal was humidity corrected and calibrated against the LOPAP HONO measurements. HONO measurements using PTR-MS have been reported previously in Koss et al. (2018) and are described in more detail in Crilley et al. (2021).

2.4 Measurement of photolysis rates

The model was constrained with the measured photolysis frequencies j(O1D), j(NO2), and j(HONO), which were calculated from the measured wavelength-resolved actinic flux and published absorption cross sections and photodissociation quantum yields (Whalley et al., 2021). All other photolysis rates are parameterised as a function of solar zenith angle using a two-stream isotropic scattering model as described by Jenkin et al. (1997) and Saunders et al. (2003). In each case, clear-sky variation of a specific photolysis rate (j) with solar zenith angle (χ) can be described well by Eq. (1):

(1) j = l cos ( χ ) m × e - n sec ( χ ) ,

where the parameters l, m, and n optimised for each photolysis rate (see Table 2 in Saunders et al., 2003, and, last access: March 2020).

For species which photolyse at near-UV wavelengths (< 360 nm), such as HCHO and CH3CHO, the photolysis rates were calculated by scaling to the ratio of clear-sky j(O1D) to observed j(O1D) to account for attenuation by clouds and aerosol. For species which photolyse further into the visible, the ratio of clear-sky j(NO2) to observed j(NO2) was used.

2.5 Model description

A tailored zero-dimensional chemical box model of the lower atmosphere, incorporating a subset of the Master Chemical Mechanism (MCM v3.3.1; Jenkin et al., 2015) into the AtChem2 modelling toolkit (Sommariva et al., 2020), was used to identify the main drivers of in situ O3 production, P(O3), in Delhi. The MCM describes the detailed atmospheric chemical degradation of 143 VOCs, though 17 500 reactions of 6900 species. More details can be found on the MCM website (, last access: March 2020). Observations of 86 unique VOCs, NOx, CO, SO2, and total aerosol surface area (ASA), along with 34 observationally derived photolysis rates, temperature, pressure, and relatively humidity, were averaged or interpolated to 15 min data and used to constrain the model. Some measured VOCs are not described in the MCM. For these species, an existing mechanism in the MCM was used as a surrogate mechanism. Surrogate species were selected based on their structural similarly to the species of interest, and reaction rates with OH, HO2 and NO3 were amended to values found in the IUPAC atmospheric chemical kinetics database (Atkinson et al., 2006) and Atkinson and Arey (2003) (see Table S1 in the Supplement). A fixed deposition rate of 1.2×10-5 s−1 was used, giving model-generated species a lifetime of approximately 24 h.

The model was constrained to an adjusted value of the observed HONO to account for high surface concentrations and an expected decline in concentration with height. The rate of vertical transport of chemical species through turbulence, or Deardorff velocity (w), was calculated during a 2-week follow-on chemical flux campaign in early November 2018 at IGDTUW, after the ground-level measurement period ended. Measurements made from a 30 m tower allowed for the calculation of w, from which an hourly averaged diel has been used to determine an approximate rate of vertical HONO transport. Observed HONO concentrations ([HONO]meas) were converted to an adjusted HONO concentration ([HONO]input) with Eq. (2):

(2) [ HONO ] input = [ HONO ] meas e - j ( HONO ) t ,

where t is calculated as the time for [HONO]meas to travel to half the boundary layer height at the measured Deardorff velocity, w. The first-order loss of HO2 to aerosol surface area (k) was calculated with Eq. (3).

(3) k = ω A γ 4 ,

where ω is the mean molecular speed of HO2 (equal to 43 725 cm s−1 at 298 K), γ is the aerosol uptake coefficient (0.2 is used here as recommended by Jacob, 2000), and A is the measured aerosol surface area in cm2 cm−3 calculated (as NSD πd2) from hourly number size distributions (NSDs). Aerodynamic diameter NSDs, Na (>352 nm), were collected using a GRIMM 1.108 (portable laser aerosol spectrometer and dust monitor) and merged into concurrently measured electric mobility diameter NSD, Nd (14–640 nm – measured using a TSI scanning mobility particle sizer (SMPS), consisting of a 3080 Electrostatic Classifier, 3081 differential mobility analyser (DMA), and 3775 condensation particle counter (CPC)) using the algorithm developed by Beddows et al. (2010).

The base reference model was run for 5 d, with each day being constrained to the diel campaign-averaged observations, to allow for the spin-up of model-generated species. Only the fifth day was taken for analysis to ensure steady state was reached. A range of additional modelling scenarios, based on this base model, were used to investigate the sensitivity of O3 production, P(O3), to changes in VOCs, NOx, aerosol surface area, and/or photolysis rates:

  • Scenario 1: vary both total VOC and NOx by a factor. VOC and NOx observations were multiplied by a factor of 0.01, 0.1, 0.25, 0.5, 0.75, 0.9, 1.1, 1.25, 1.5, 1.75, and 2.

  • Scenario 2: vary individual primary VOCs and NOx by a factor. A VOC of interest and NOx observations were multiplied by a factor of 0.01, 0.1, 0.25, 0.5, 0.75, 0.9, 1.1, 1.25, 1.5, 1.75, and 2. Observed carbonyls (all oxygenated compounds excluding alcohols) were not constrained in this scenario to allow for secondary compounds to be varied with changes in concentrations of their precursor VOC.

  • Scenario 3: vary by VOC class and vary NOx. Individual VOC classes were multiplied by a factor of 0.2, 0.4, 0.6, 0.8, 1.2, 1.4, 1.6, 1.8, and 2. NOx was multiplied by a factor of 0.25, 0.5, and 0.75. All possible combinations of altered VOC class and NOx observations were constrained in the model. The remaining VOCs not grouped into the class of interest are constrained at their observed values. Observed carbonyls are not constrained in this scenario to allow for secondary compounds to be varied with changes in concentrations of their precursor VOC.

  • Scenario 4: increase individual VOC of interest by 5 %. The constrained concentration of a VOC species of interest was increased by 5 % by molar mass following the procedure carried out in Elshorbany et al. (2009). Comparison of the change in P(O3) against the base model allows for the determination of ΔP(O3)increm (Sect. 3.4). Observed carbonyls are not constrained in this scenario to allow for secondary compounds to be varied with changes in concentrations of their precursor VOC.

  • Scenario 5: vary ASA. ASA was multiplied by a factor of 0.1, 0.3, 0.5, 0.7, 0.9, 1.1, 1.3, 1.5, 1.7, 1.9, and 2. The difference in P(O3) between each model and the base model was examined.

  • Scenario 6: vary photolysis rates. All 34 photolysis rates were multiplied by a factor of 0.1, 0.3, 0.5, 0.7, 0.9, 1.1, 1.3, 1.5, 1.7, 1.9, and 2. The difference in P(O3) between each model and the base model was examined.

  • Scenario 7: vary both individual primary VOCs and NOx by Emissions Database for Global Atmospheric Research (EDGAR) source sector. VOCs were grouped and observed concentrations were divided into each source sector as described by the EDGAR v5.0 Global Air Pollutant Emissions and EDGAR v4.3.2_VOC_spec inventories (, last access: November 2020). The proportions of VOC and NOx concentrations attributed to each individual source were multiplied by a factor of 0.01, 0.1, 0.25, 0.5, 0.75, and 0.9. Observed carbonyls were not constrained in this scenario to allow for secondary production of these compounds to be varied with changes in concentrations of their precursor VOC. For each of these model runs, there were six variations where monoterpenes were assumed to be 0 %, 10 %, 25 %, 50 %, 75 %, and 100 % anthropogenic.

3 Results

3.1 Observed NOx, CO, and O3

The observed mixing ratio (ppbV) time series of NO, NO2, CO, and O3 during the campaign are presented in Fig. 2. High concentrations of NO and CO were observed towards the latter half of the month, with high daytime O3 concentrations observed throughout the entire measurement period.

Figure 2Observed mixing ratio time series of (from top to bottom): NO, NO2, CO, and O3 during October 2018.


The diel profile of NO, and to a lesser extent CO, has a U-shaped profile, with much higher concentrations observed at night compared to the day (Fig. 3). This profile results from a shallow and stagnant nocturnal boundary layer, described in more detail in Stewart et al. (2021a). Average daytime concentrations (06:00–18:00 IST, Indian standard time, roughly in line with sunrise and sunset times) of NO and CO were 58.8 ppbV and 1.2 ppmV, respectively. These campaign-averaged diel profiles have large standard deviations, representing the day-to-day variation in concentration throughout October (Fig. 2). Average nighttime concentrations (18:00–06:00) of NO and CO were 247.0 ppbV and 2.9 ppmV, respectively, with NO and CO concentrations occasionally exceeding 800 ppbV and 10 ppmV, respectively, in the latter half of the month (Fig. 2). The NO2 profile shows two peaks, at approximately 09:00 and 18:00, perhaps due to increased commuting traffic during these times. The O3 diel profile peaks at approximately 13:00, with a mean peak concentration of 78.3 ppbV. A rapid increase in O3 concentration is observed first thing in the morning from around 05:00, which levels off at around 06:00 before again rapidly increasing from 08:00 to a peak concentration  13:00. As the precursor species of O3 depend on light to undergo the chemical processing that leads to O3 production, this profile suggests that as the Sun rises rapid photochemical formation of O3 is initiated by the photolysis of precursors that have accumulated in the lower atmosphere during the nighttime. The early morning increase in O3 concentrations may also be influenced by the reduced NO titration as the boundary layer height increases (Fig. 3).

Figure 3October campaign-averaged diel profiles. The shaded ribbon represents the standard deviation of this average.


Average mixing ratios of NOx and VOCs observed during the post-monsoon Delhi campaign are comparable to those observed in 1970s Los Angeles. High concentrations of NOx ( 200 ppbV) and VOCs (e.g. benzene  10 ppbV; toluene  30 ppbV) led to large concentrations of O3 in the US city, averaging at around 400 ppbV (Pollack et al., 2013). Despite the high O3 concentrations observed in Delhi, O3 concentrations in post-monsoon Delhi did not exceed 150 ppbV and are of similar magnitude to observed O3 in Beijing, Shanghai, and Guangzhou despite much higher observed VOC and NOx concentrations (Tan et al., 2019). The lower O3 observed in Delhi compared to 1970s Los Angeles can be attributed to differences in both topography and meteorology. The isolated coastal city of Los Angeles is surrounded by mountains to the north and east, with the prevailing wind dominating from the coast. Due to its basin-like topography, and with a cool onshore sea breeze often creating a temperature inversion, the air mass circulates within the city and the transport of emissions out of the basin is impeded. Although landlocked Delhi lies to the southwest of the Himalaya, the city is very flat and resides far enough away from the mountain range to allow for the efficient transportation of air masses from the city. It is also important to consider that the very high concentrations of NOx and VOCs observed peak to comparable concentrations to Los Angeles during the evening and at night, where they are trapped due to a shallow, stagnant boundary layer, and there is little to no photochemical activity. O3 production rates in Delhi peak in the morning, when concentrations of pollutants, though still high, are much lower than at night due to rapid boundary layer expansion (see Sect. 3.4).

Although 12 key air pollutants, including O3, have prescribed national ambient air quality standards (Dube, 2009), only four have been identified for regular and continuous monitoring: sulfur dioxide (SO2), nitrogen dioxide (NO2), respirable suspended particulate matter (RSPM, or PM10), and fine particulate matter (PM2.5). The prescribed national standards for O3 are an 8-hourly limit of 100 µg m−3 ( 50 ppbV) and 1-hourly limit of 180 µg m−3 ( 90 ppbV). The observed O3 8 h averages between 09:00 and 17:00 throughout the measurement period are shown in Fig. 4, along with the hourly averaged maximum O3 concentrations.

Figure 4The 8 h mean O3 concentration between 09:00 and 17:00 IST (a) and hourly averaged maximum daily O3 concentration (b) during October 2018. The dashed black line shows the prescribed national standards for O3.


Our observations show that the 8 h O3 prescribed national standard is exceeded on 16 d during our 24 d measurement period (67 % of days), and the 1 h max is exceeded on 8 d (33 % of days). The published national standards state that any pollutant which exceeds the prescribed values for two consecutive days qualifies for regular and continuous monitoring. As there are up to four consecutive days in which the 8 h daytime mean O3 limit was breached during our campaign, this implies that O3 should be continuously and regularly monitored in this part of Old Delhi.

3.2 Observed VOCs

Alkanes were the predominant VOC class contributing to the total measured VOC mass concentration (42 %), consistent with previous observations at the site (Shivani et al., 2018), followed by alcohols (18 %), aromatics (17.8 %), and carbonyls (13.9 %). The percentage contributions of each VOC class to total measured VOCs, along with mean mass concentrations of the top 10 contributors, are presented in Fig. 5. The top individual species contributing to the VOC mass concentrations were ethanol (10.7 %), n-butane (9.9 %) and n-propane (7.9 %).

Figure 5Percentage contribution of different VOC classes to the total mean measured NMVOCs (a) and the mean concentrations of the top 10 contributors to total measured VOC concentrations (b) during the campaign in µg m−3. Colours correspond to the different VOC classes.


The general diel profile for all VOCs, excluding isoprene and some oxygenated species, is U-shaped (Fig. 6). This U-shaped profile is also observed for NO and CO. Concentrations are much higher during the night, and lower in the day, as they are concentrated in a stagnant, shallow boundary layer that forms over the city at night (Stewart et al., 2021a), and are subject to photochemical losses during the day. The U shape is less apparent for acetone, methanol, and ethanol. This may be a result of very high emissions of these compounds during the day and/or formation through secondary chemistry (Stewart et al., 2021a). Isoprene has a typical biogenic diurnal profile, in contrast to monoterpenes which have a similar trend to other anthropogenic species (see α-terpinene, Fig. 6). A separate study which conducted PMF analysis on the co-located PTR-QiTOF flux measurements, resolved two traffic-related factors (Cash et al., 2021). A large proportion of the total measured monoterpenes were resolved within traffic factors ( 60 %) and traffic emissions dominated throughout the campaign ( 50 % of total VOC emissions). Gkatzelis et al. (2021) have recently shown that monoterpene measurements in New York City can be predominantly apportioned to fragranced volatile consumer products (VCPs) and other anthropogenic sources (e.g. cooking and building materials), and that although biogenic sources are the dominant source of monoterpenes globally, in urban environments monoterpene fluxes from fragranced VCPs can compete with the emissions from local vegetation (see Sect. 3.8 for further discussion of anthropogenic sources of monoterpenes).

The standard deviation of all the aggregate VOC diel profiles is large due to large variations in concentrations from day to day throughout the month. Generally, high concentrations of all species are observed at the end of October and lower concentrations at the start (see Figs. S1 and S2 in the Supplement).

Figure 6Mean campaign-averaged diel mixing ratio profiles of selected VOCs during the campaign: ethane, ethene, n-butane, isoprene, benzene, α-terpinene, methanol, ethanol, acetone. The shaded ribbon represents 1 standard deviation from the mean. Colours correspond to the different VOC classes.


3.3 OH reactivity, k(OH)

Speciated VOC concentrations alone do not indicate which individual compounds are important for O3 formation. It is crucial to account for their reactivities with atmospheric oxidants, and their structure, to gain insight into each individual species and compound class contribution to in situ P(O3). All VOCs react with OH, leading to peroxy radical formation. These peroxy radical species (HO2 and RO2) mediate the conversion of NO to NO2, leading to P(O3) (Reactions R1–R4). The rate at which VOCs react with OH is thus the rate-determining step in the amount of O3 formed.

The chemical box model described in Sect. 2.5 was used to investigate the total OH reactivity, expressed as k(OH) – a first-order loss rate in units of s−1, of observed precursors to O3 (base reference model). NOx, CO, and individual VOC class contributions to k(OH) are presented in Fig. 7, along with the k(OH) of unmeasured species generated by the model (referred to as “model-generated species” from now on). VOCs and model-generated species represented 67.4 % of the total k(OH), with alkenes (9.6 %) and aromatics (8.8 %) being the largest VOC class contributors. The k(OH) value was higher during the night, showing an inverse relationship with boundary layer height. This is typical for urban environments where nighttime emissions are typically released into a shallow boundary layer (discussed further in Stewart et al., 2021a).

Nighttime boundary layer heights in Delhi were very low, leading to a clearly defined k(OH) profile with a maximum of  250 s−1 at around 21:00 and a minimum of  57 s−1 at 14:00. The campaign average boundary layer height range was approximately 39–1550 m. Generally, k(OH) in megacities is found to peak at around 06:00, consistent with morning emissions into a shallow boundary layer. A small peak is also observed around this time in Delhi, but the largest peak is calculated at around 21:00. NOx represents  35 %–40 % of the total k(OH) at these times, perhaps due to high volumes of traffic. This is supported by VOC traffic PMF factors described in Cash et al. (2021), which peaks during the evening ( 19:00–21:00), and account for  87 % of the total emissions at this time.

The values of k(OH) determined in this work are significantly higher than those observed in other megacities, with a daytime minimum k(OH) more than double that observed in Beijing during the summer of 2017 (Whalley et al., 2021). Previous studies in New York City, Mexico City, and Tokyo have observed k(OH) > 100 s−1 (Ren et al., 2006; Shirley et al., 2006; Yoshino et al., 2006). A maximum summertime k(OH) of 116 s−1 was observed in London during rush hour, but lower OH reactivities of 22–37 s−1 were more typical during the campaign (Whalley et al., 2016). In a study during summertime in Seoul, average k(OH) was approximately 15 s−1 in the afternoon, increasing to 20 s−1 at night (Kim et al., 2016). High OH reactivities were observed in the Pearl River Delta, with mean maximum reactivities of 50 s−1 at daybreak, and mean minimums of 20 s−1 observed at noon (Lou et al., 2010). Average summertime k(OH) observations from Beijing, China, observed a k(OH) maximum of  37 s−1 at around 06:00 and daily minimum of k(OH)  22 s−1 at 15:00 (Whalley et al., 2021), within a boundary layer range of approximately 150–1500 m.

Figure 7Diel profile of campaign-averaged VOCs, CO, and NOx contributions to OH reactivity, k(OH).


In Delhi, O3 concentrations rapidly increase from  08:00, peaking around  13:00 before declining in the afternoon (Fig. 3). A breakdown of the percentage contribution of each species class to morning k(OH) (08:00–12:00) is presented in Fig. 8. The average morning k(OH) is dominated by model-generated species (31.9 %), followed by NOx (21.5 %) and carbonyls (9.2 %).

Figure 8Relative contribution of species classes to average morning OH reactivity, k(OH), (08:00–12:00). The contribution of each group to the total k(OH) including NOx and model-generated species (a) and the contribution of each class to just the VOC proportion of k(OH) (b).


3.4 O3 production potentials

Model-generated species are the products of reactions of VOCs of all classes, meaning the overall contribution of individual classes of VOCs to k(OH) is underestimated in Figs. 7 and 8 (Sect. 3.3). To assess the true contribution of VOCs to k(OH) and production, several model runs where each constrained VOC class was increased by 5 % were compared to the base reference model (scenario 4, Sect. 2.5).

One way to assess the contributions of different VOCs to O3 production is to determine the change in O3 production, ΔP(O3), when the base reference model is compared to an adapted model where the VOC of interest is increased by 5 % (Elshorbany et al., 2009). The additional P(O3) resulting from the 5 % perturbed system was averaged over 24 h. For both the base and 5 % scenarios, the net rate of O3 production, P(O3), is calculated by subtracting the instantaneous rate of O3 loss, L(O3)inst, by the instantaneous rate of O3 production, P(O3)inst, via Eqs. (4)–(6). P(O3)inst was calculated by determining the NO2 production rate through reactions of NO with HO2 and RO2 (Whalley et al., 2018), assuming the production of a molecule of NO2 equates to the production of a molecule of O3 via Reaction (R5). NO2 loss processes which do not yield O3, such as removal by reaction with OH and reaction with acyl peroxy radicals to form peroxy acetyl nitrates (PANs), must also be accounted for. In these calculations, OH, HO2, and RO2 concentrations are generated from the observationally constrained model.

P(O3)inst is the rate of NO oxidation by HO2 and RO2 radicals:

(4) P ( O 3 ) inst = ( k HO 2 + NO HO 2 NO + i k RO 2 i + NO RO 2 NO ) .

L(O3)inst is the rate of loss of O3 though reactions with OH, HO2, and photolysis followed by a reaction with H2O vapour, along with the loss of NO2 through reactions with OH:

(5) L ( O 3 ) inst = j O 1 D O 3 f + k OH + O 3 OH O 3 + k HO 2 + O 3 O 3 HO 2 + k OH + NO 2 + M NO 2 OH M + i k RO 2 i + NO 2 + M [ RO 2 ] NO 2 M ,

where f is the fraction of O(1D) atoms (formed in the photolysis of O3) that react with H2O vapour to form OH, rather than undergo collisional stabilisation. Net O3 production, P(O3), can then be calculated as

(6) P ( O 3 ) = P ( O 3 ) inst - L ( O 3 ) inst .

The O3 production resulting from an incremental increase of each VOC between 08:00 and 12:00, P(O3)increm, was then calculated with Eq. (7):

(7) Δ P ( O 3 ) increm = Mean ( P ( O 3 ) i - P ( O 3 ) base ) ,

where P(O3)i is the mean O3 production between 08:00 and 12:00, calculated from the model run where the VOC of interest is increased by 5 %. Using this approach, the 10 VOCs contributing to the greatest change in P(O3) on an incremental increase, P(O3)increm, are isoprene, α-terpinene, ethene, toluene, propene, α-phellandrene, ethanol, 1,2,4-trimethylbenzene, CO, and β-ocimene (Table 1).

Table 1The top 10 highest ΔP(O3)increm VOCs, and their respective classes, averaged between 08:00 and 12:00.

Download Print Version | Download XLSX

3.5P(O3) sensitivities to VOC / NOx ratios

To better understand the complex non-linear chemistry at play, a box model is used to probe the chemical sensitivities of observed precursors to O3 formation. The model was run 144 times, each with adjusted concentrations of VOCs and NOx (scenario 1). The observed concentrations were multiplied by a factor to generate unique model runs from which P(O3) was calculated using Eqs. (4)–(6). The resultant P(O3) isopleth is presented in Fig. 9.

Figure 9Mean modelled O3 production, P(O3), isopleth between 08:00 and 12:00 based on varying VOC and NOx concentrations in the model. The red diamond at point (1,1) represents modelled P(O3) at observed VOC and NOx concentrations.


The modelled VOC–NOx P(O3) isopleth supports the assignment of Delhi being, on average, in a VOC-sensitive photochemical regime (Sillman et al., 1990), with the diel profile of O3 production peaking at 09:00 (see the Supplement, Fig. S3). Reducing NOx alone and maintaining VOC concentrations would result in an increase in P(O3). Therefore reducing, or even maintaining, O3 levels in the future will require a reduction in VOC emissions, if future emission control measures continue to target NOx emissions in Delhi. This is consistent with a study of observational data in Delhi, whereby a SARS-CoV-2 lockdown led to a reduction in NOx emissions and concentrations and an increase in the concentration of O3 (Jain and Sharma, 2020). To implement an efficient and realistic VOC reduction plan, the key VOCs contributing to P(O3) need to be identified, along with their sources.

As has been previously discussed, O3 concentrations limits of 50 ppbV were regularly exceeding during the campaign, with the maximum daily 8 h averages peaking at 88 ppbV (Fig. 4). To successfully reduce O3 to the limit of 50 ppbV, O3 production must be reduced by 56 %. This can be achieved by reducing NOx by 25 %, 50 %, and 75 % along with a concurrent reduction in VOCs of 48 %, 61 %, and 78 %, respectively. Alternatively, if VOCs were halved, NOx reductions could not exceed 29 % to peak O3 below 50 ppbV. To obtain a reduction in O3 production without reducing VOCs would require a NOx reduction of at least 92 %. However, it is important to consider that reducing emissions of NOx and VOCs by these values will only impact the in situ formation of O3. Regional O3 production leading to the transportation of O3 from outside of Delhi is beyond the scope of this analysis and must also be taken into consideration by policy makers.

3.6P(O3) sensitivity to VOCs: by class

The impact of changing VOC concentrations (by class) on mean P(O3) was investigated using scenario 3 (Sect. 2.5). For each model run, the constrained concentrations of all species in the class of interest were multiplied by a “VOC factor”. Modelled P(O3) was calculated upon changing “VOC factor” for CO and six different VOC classes: alkanes, alkenes, aromatics, monoterpenes, isoprene, and alcohols. Changes in P(O3) were found to be most sensitive to changes in aromatics, followed by the monoterpenes and alkenes (Fig. 10). Halving aromatic, monoterpene, and alkene concentrations independently (reducing a VOC factor of 1 to a VOC factor of 0.5, Fig. 10) reduced modelled P(O3) by 15.6 %, 13.1 %, and 12.9 %, respectively. However, future air quality control strategies are also likely to include a co-reduction in NOx. As we observed in Fig. 9, a reduction in NOx coupled with an insufficient reduction in VOCs is likely to increase P(O3) in Delhi, under VOC-limited conditions. On concurrently reducing VOC class and NOx by half (factor of 0.5, Fig. 10), aromatics, alkenes, and monoterpenes lead to the smallest increase in P(O3) (24.9 %, 29.8 %, and 35.5 % respectively). However, this still represents a significant increase in P(O3). This suggests targeting one VOC class alone in future pollutant reduction strategies is insufficient, and that reducing a source/sources emitting multiple VOC classes is important to avoid an increase in P(O3) on simultaneously reducing NOx.

Figure 10Modelled mean morning (08:00–12:00) change in P(O3) rate for CO and six VOC classes upon varying their concentrations by multiplying by the “VOC factor” (a) and by varying just NOx (b).


3.7 The impact of aerosol surface area on P(O3)

Aerosol particles in the troposphere can impact gas phase chemistry and photochemical activity, affecting P(O3), in a number of ways. Aerosols can interact with incoming solar radiation, either scattering or absorbing sunlight. The precise impact of aerosol is dependent on a range of factors including chemical composition, particle size distribution, and phase state. It has previously been shown that in highly polluted urban areas, attenuation of the actinic flux due to aerosol absorption can significantly reduce photolysis rates, and hence P(O3), by up to 30 % (Castro et al., 2001; Hollaway et al., 2019; Real and Sartelet, 2011; Wang et al., 2019). Conversely, aerosol scattering can potentially increase P(O3) by increasing the photolysis rates (Dickerson et al., 1997; He and Carmichael, 1999).

Aerosol can also participate in heterogeneous chemistry; i.e. uptake of radicals to the aerosol or reactions at the aerosol surface can affect gas phase radical budgets (George et al., 2015). A recent regional modelling study (K. Li et al., 2019) linked increasing O3 concentration trends between 2013 and 2017 in Beijing to decreasing PM2.5 concentrations. The study attributed this to the decreased uptake of HO2 radicals to aerosol particles, leading to increased HO2 available to participate in P(O3) (see Reaction R4). However, an experimental study (Tan et al., 2020) in the North China Plain did not observe the effect. The relationship between aerosol optical depth (AOD) and photolysis rates is strongly non-linear and has a larger impact at low solar zenith angles (i.e. in the morning/evening/high latitudes) (Wang et al., 2019).

To investigate the potential impact of aerosol related processes in Delhi, ASA and photolysis rates were varied independently in the model using a scaling factor (“Factor”, Fig. 11), and the resulting impact on P(O3) calculated. The first-order loss of HO2 to aerosol surface area (k) was calculated using Eq. (3) (scenario 5, Sect. 2.5).

Figure 11Impact of varying photolysis rates and ASA on mean modelled morning O3 production rates, ΔP(O3) (08:00–12:00).


Changes in ASA, with respect to HO2 uptake, were found to have minimal impact on the modelled mean morning P(O3) rate (08:00–12:00), as shown in Fig. 11, in line with the observations of Tan et al. (2020). The lifetime of HO2 is dominated by its reaction with NO, within the campaign-averaged NO concentrations observed. Changes to the photolysis rates have a large impact on P(O3) (Fig. 11), with P(O3) increasing roughly linearly with increasing photolysis rates. The impact of changes to AOD on photolysis rates are non-linear and dependent on solar zenith angle (SZA) (i.e. time of day). This work shows that the majority of daily P(O3) occurs between 08:00 and 12:00, representing an SZA range of 35–80. Aerosol loading is expected to reduce in Delhi going into the future if India implements air pollution controls in line with other countries. Based on the work of Wang et al. (2019), reductions in PM, leading to reductions in AOD, could potentially increase j(NO2) by up to 30 % at midday and up to 100 % at 08:00. This will lead to increased photochemical activity and higher P(O3).

From an air quality perspective, it is important to consider reducing not only VOCs but also NOx concentrations in Delhi. Along with this, reductions in PM should also be considered. As we have seen in Fig. 11, reducing ASA does not significantly impact P(O3) through heterogenous chemistry. However, an ASA reduction will result in a decrease in AOD, which will in turn increase photolysis rates. A recent study by Chen et al. (2021) modelled the impact of AOD on photolysis in Delhi during November 2018, and found that halving the AOD results in a 14 % increase in j(NO2). In Delhi, total VOCs would need to be reduced by  50 % to maintain current P(O3) rates if NOx were halved and photolysis rates increased by 14 %. However, the impact of reducing aerosol on photolysis rates is complex. Reductions in aerosol at ground level may lead to either increased or decreased photolysis rates near the surface, depending on the scattering properties of the aerosol. Changes in photolysis rates from increased or decreased aerosol loading may also vary throughout the depth of the boundary layer. The box model assumes photolysis rates are uniform throughout the boundary layer and that aerosols are well mixed. A more detailed study on the temporal and spatial patterns of aerosol in Delhi and its impact on photolysis rates is required to fully assess the aerosol impact on in situ O3 production (Castro et al., 2001).

It is important to note that this study focuses on the sensitivity of P(O3) to ASA through HO2 uptake only. Additional chemical consequences and feedbacks of decreasing aerosol, such as changes to HONO concentrations, have not been accounted for here. With the high levels of NO2 observed in Delhi, the heterogenous conversion of NO2 to HONO on particle surfaces may be an important mechanism (Liu et al., 2014; Lee et al., 2016; Tong et al., 2016; Lu et al., 2018b). HONO reductions from decreased ASA may lead to reduced OH radical formation in Delhi; thus, the impact of ASA reduction on P(O3) may be underestimated in this work.

3.8 Discussion

If future air pollution controls in Delhi follow the air quality strategies prevalent or planned in many other countries, including the EU, there is a danger that urban O3 concentrations could significantly increase, unless careful consideration is given to the specific atmospheric chemistry occurring in the city. Many other major cities across the globe have focused their air pollution abatement strategies on reducing NOx and particulate emissions from traffic sources. Urban NOx emissions (at tailpipe) are likely to decline over time, as a result of improved exhaust gas treatments, the turnover of the fleet to newer, less polluting vehicles, and the increasing prevalence of electric vehicles (Molina, 2021). While this is also likely to reduce ambient VOC concentrations through reduction in tailpipe and evaporative emissions, this reduction may be smaller than for NOx. The magnitude of resulting changes in P(O3) from reduced road transport emissions may depend on the proportion of VOCs emitted from road transport relative to other, non-vehicular sources. As demonstrated here, reductions in NOx without sufficiently reducing VOC emissions may lead to large increases in P(O3) under a continued VOC-limited regime. Reducing traffic emissions will also likely lead to reduced aerosol loading in Delhi. Our study suggests a reduction in aerosol surface area will have very little direct effect on P(O3) via heterogenous chemistry, as HO2 reactivity is dominated by its reaction with NO. However, an ASA reduction is likely to increase the amount of sunlight reaching the boundary layer and hence photolysis rates. This will lead to a subsequent increased P(O3), as discussed in more detail by Chen et al. (2021).

P(O3) in Delhi was found to be most sensitive to reductions in aromatics and alkenes, and so monitoring the abundance and knowing the sources of these compounds in Delhi is essential for implementing effective pollutant reduction strategies to avoid a future rise in P(O3). As these classes are thought to come mainly from traffic sources, it is possible that reducing road transport emissions may reduce traffic-derived VOCs sufficiently, along with reducing NOx, to perturb a large increase in P(O3). However, the proportion of VOCs in Delhi, and at IGDTUW, attributed to traffic sources is uncertain, and thus the extent to which road transport reductions will impact P(O3) is unknown.

More information of the effect of reductions in NOx and VOCs by source can be obtained by varying sources described by emissions inventories. Anthropogenic emissions of NOx, CO, and VOCs by source in Delhi are available from the EDGAR v5.0 Global Air Pollutant Emissions and EDGAR v4.3.2_VOC_spec inventories ( According to the inventory, all pollutants investigated in this analysis can be almost entirely described by five source sectors (Table 2): road transport (RT); railways, pipelines, and off-road transport (RPORT); energy for buildings (EB); combustion for manufacturing (CM); and process emissions (PE). However, it should be noted here that the EDGAR emissions inventory describes a coarse, low-spatial-resolution city-wide representation of emissions from Delhi. For this analysis, we assume the EDGAR sector split ratios for Delhi are representative of those at the IGDTUW measurement site.

Model-constrained concentrations were varied by their contributions to sources in the EDGAR emissions inventory (scenario 7, Sect. 2.5). Data from the inventory represent average annual emissions from source sectors, last updated for VOCs in 2012 and 2008 for monoterpenes, with 0.1× 0.1 spatial resolution. VOCs were assumed to be described in full by all sources available from the inventory, with no biogenic influence. Isoprene was excluded from the analysis as it is assumed to have an entirely biogenic source. As a significant anthropogenic signature is seen from the monoterpenes, and because of their significance for P(O3) in this study, 50 % are assumed to be from anthropogenic sources found in the EDGAR emission inventory for the purpose of this analysis. A full list and description of sources used from the EDGAR inventory can be found in the Supplement. Figure 12 shows the change in P(O3) rate when reducing VOCs and NOx contributions from individual source sectors.

Table 2Relative proportion of pollutants emitted from five EDGAR emission inventory source sectors: road transport (RT); railways, pipelines, and off-road transport (RPORT); energy for buildings (EB); combustion for manufacturing (CM); and process emissions (PE).

* The percentage of monoterpenes attributed to these sources in varied in the model and is assumed to be 50 % in the base case.

Download Print Version | Download XLSX

Figure 12Change in P(O3) on reducing RT, CM, PE, and EB (a); and (b) reducing road transport simultaneously along with either (i) CM, (ii) PE, (iii) EB, (iii) CM, or PE. The shaded region, bounded by a dashed line, represents the variability in ΔP(O3) when 0 %–100 % of observed monoterpenes contribute to anthropogenic EDGAR source sectors that include PE. The solid lines represent the base case whereby 50 % of monoterpenes are assumed to be from anthropogenic sources described in the EDGAR inventory (the other 50 % are assumed to be from biogenic sources).


Based on this analysis, reducing the RT source in isolation results in increased P(O3), whilst reducing CM, PE, and EB independently leads to decreased P(O3). This is explained by CM, PE, and EB being major sources for VOCs, and RT also being a major source of NOx. Although reducing RT also reduces some VOCs, particularly aromatic and higher alkane species, there is not sufficient VOC reduction to compensate for the large co-reduction in NOx, leading to increased P(O3) (Fig. 10). Reducing RT by 100 % would still lead to a modelled P(O3) increase of  40 ppb h−1 (Fig. 12a). It is important to acknowledge the uncertainties in this analysis. The EDGAR emissions inventory attributes a low percentage split of VOCs to the RT source sector at the city level, perhaps underestimating the proportion of VOCs that would be reduced with RT reductions at the measurement site.

Air quality mitigation strategies should therefore focus on reducing RT, alongside one or more additional major VOC sources. For example, Fig. 12b shows P(O3) begins to reduce when RT emissions are reduced along with PE. These reductions are even greater when RT is reduced with both PE and CM sources. When RT is reduced by 50 % alongside PE and PE + CM, modelled P(O3) is reduced by  10 ppb h−1 and  20 ppb h−1, respectively. The VOCs that contribute the most to PE and CM emissions are n-butane, propane, alcohols, toluene, and xylenes. Although alcohols and alkanes were not identified as the key VOC classes contributing to P(O3) (Sect. 3.6), they are the two largest classes of VOCs observed in Delhi by mass, making up 42 % and 18 % of the total measured VOCs, respectively (Fig. 5). Reducing the PE source thus leads to large reductions in VOCs by mass, helping to drive down P(O3). However, it is worth noting that the effectiveness of reducing the RT + PE source on modelled P(O3) is dependent on the proportion of anthropogenic monoterpene emissions in Delhi. According to the EDGAR emission inventory, 98.1 % of anthropogenic monoterpenes in Delhi are attributed to process emissions (PE) (Table 2). These emissions include sources such as emissions from chemical industry, and other industrial processes, and include solvent emissions and emissions from product use (see Table S3 in the Supplement). As emissions from these sources are grouped together in the EDGAR inventory, the exact sources that monoterpenes are attributed to cannot be identified. The sensitivity of ΔP(O3) from reducing process emission sources (PE, RT + PE, and RT + PE + CM) is shown by the shaded regions in Fig. 12, where the dashed lines represent the sensitivity limits where the observed monoterpenes are between 0 % and 100 % anthropogenic (as opposed to biogenic). There is relatively little impact on P(O3) on reducing RT + PE when monoterpenes are assumed to have an entirely biogenic source. However, it is clear that although the degree to which reducing process emissions along with road transport in this study impacts P(O3) cannot be accurately determined, even if no monoterpenes are reduced within this source, reducing it does not negatively impact P(O3). It is also important to consider possible underestimations for the monoterpene contribution to RT. Monoterpene observations in Delhi were strongly correlated with CO emissions, suggesting anthropogenic sources (Stewart et al., 2021a). The EDGAR emissions inventory assigns 0 % of the anthropogenic monoterpenes in the inventory to the RT source sector (Table 2). An analysis of the PTR-QiTOF flux data, obtained at the IGDTUW measurement site directly after the concentration measurement period ended, suggests  60 % of the monoterpenes observed could be attributed to traffic factors (Cash et al., 2021). A study by Wang et al, 2020 suggested vehicular and burning sources may contribute to the anthropogenic emissions of biogenic molecules. Other possible sources of anthropogenic monoterpenes in Delhi are emissions from cooking herbs and spices, and from fragrances and personal care products (Klein et al., 2016; McDonald et al., 2018).

It should be noted that the EDGAR inventory data used in this study are cropped to the Delhi area, and thus some anthropogenic sources may be missing. One example of this is that of agricultural burning, which is frequent in areas surrounding Delhi and across the Indo-Gangetic Plain (IGP) (Jat and Gurjar, 2021; Kulkarni et al., 2020). It is likely that regional sources of O3 and its precursors will have a significant effect on Delhi's air quality, and vice versa. There are likely to be many missing sources from the EDGAR emissions inventory, along with discrepancies between the top-down and bottom-up data. In addition, we assume the ratios of VOCs and NOx in each source sector for the Delhi region are representative of the ratios emitted near IGDTUW. Whilst the OH reactivity analysis provides an instantaneous assessment of a VOC's reactivity, based on local observations with a small spatial footprint, the EDGAR Delhi-wide sectoral split represents an aggregate of a wider region and may underrepresent the RT VOC contribution observed at the measurement site. As a result, the magnitude of the modelled increased P(O3) when road transport is reduced may represent a worst-case scenario.

4 Conclusions

A detailed chemical box model constrained to an extensive observational dataset of 86 VOCs, 34 photolysis rates, NO, NO2, CO, SO2, HONO, temperature, pressure, and relative humidity was used to explore the sensitivity of photochemical O3 production, P(O3), to VOCs and NOx in the Indian megacity of Delhi. The urban measurement site at the Indira Gandhi Delhi Technical University for Women is determined to be in a VOC-limited chemical regime. Our analysis examined the sensitivity of VOC classes to mean morning P(O3), and the aromatic VOC class was identified as being the most important, with a 50 % reduction in ambient concentrations leading to a reduction in modelled morning P(O3) of 15.6 %, followed by monoterpenes and alkenes (13.1 % and 12.9 %, respectively). The direct impact of the aerosol burden on P(O3) was found to be negligible with regards to heterogeneous radical uptake. However, P(O3) was sensitive to increasing photolysis, which may result from decreasing particulate matter that is likely to arise from future emission reduction strategies. Though it is important to reduce NOx and particulate matter in all abatement strategies, reducing NOx without reducing VOCs was found to significantly increase P(O3). VOCs and NOx were also evaluated by their respective contributions to different source sectors as defined in the EDGAR emissions inventory. Reducing emissions from road transport sources alone led to increased P(O3), even when the source was removed in its entirety. This suggests a balanced approach to pollution reduction strategies is required, and multiple sources should be targeted simultaneously for effective reduction, incorporating major monoaromatic, alkene, and anthropogenic monoterpene VOC sources. Reducing road transport emissions along with combustion from manufacturing and process emissions was found to be an effective way of reducing P(O3). Monoterpenes are found to have a significant impact on P(O3) and showed a diurnal profile consistent with other anthropogenic VOCs. Future work should be carried out to determine the sources and fraction of anthropogenic emissions to the observed monoterpene concentrations measured in Delhi, with a recent study suggesting 60 % of the monoterpenes observed at the site are coming from traffic related sources (Cash et al., 2021). To further understand the complex urban atmospheric chemistry in Delhi, more model constraints are required, such as the measurement of radical species and k(OH). The potential effects of chlorine chemistry were also not investigated in this study and may have an important role in the local chemistry due to the prevalence of waste burning in the city (Gunthe et al., 2021). Satellite observations and aircraft measurements may also help develop our understanding of the regional impacts of regional agricultural biomass burning on Delhi from neighbouring states. This work highlights that a careful approach, considering the complexities of chemical processing in the urban atmosphere, is required for effective air quality improvement strategies.

Data availability

Data used in this study can be accessed from the CEDA archive: (Nemitz et al., 2020).


The supplement related to this article is available online at:

Author contributions

BSN prepared the manuscript with contributions from all authors. BSN, GJS, WSD, ARV, RED, WJA, LRC, MSA, BL, EN, and JRH provided measurements and data processing of the comprehensive suite of atmospheric pollutants used in this study. MJN, PME, ES, LKW, DEH, RS, and SC provided modelling assistance and model code development used in the chemical box model. ÜAS, DCSB, LKW, MSA, and WJB contributed to the data processing of measured chemical species. JMC contributed to scientific discussion, critical to the findings of this work. S, RG, BRG, and EN assisted with the logistics of lab setup at the measurement site, and data analysis. ACL, JFH, CNH, WJB, JRH, ARR, and JDL provided overall guidance to experimental setup and design, along with assistance in operating instrumentation, and data analysis and interpretation.

Competing interests

The authors declare that they have no conflict of interest.


The paper does not discuss policy issues, and the conclusions drawn in the paper are based on interpretation of results by the authors and in no way reflect the viewpoint of the funding agencies or institutions authors are affiliated with.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


The authors acknowledge Pawel Misztal and Brian Davison for their involvement in operating the PTR-QiTOF, and Tuhin Mandal at CSIR-National Physical Laboratory for his support in facilitating the measurement sites used in this project. This work was supported by the Newton Bhabha fund administered by the UK Natural Environment Research Council through the DelhiFlux and ASAP projects of the Atmospheric Pollution and Human Health in an Indian Megacity (APHH-India) programme. The authors gratefully acknowledge the financial support provided by the UK Natural Environment Research Council and the Earth System Science Organization, Ministry of Earth Sciences, Government of India, under the Indo-UK Joint Collaboration (DelhiFlux). Beth S. Nelson and Gareth J. Stewart acknowledge the NERC SPHERES doctoral training programme for studentships. James M. Cash is supported by a NERC E3 DTP studentship.

Financial support

This research has been supported by the Natural Environment Research Council (grant nos. NE/P016502/1 and NE/P01643X/1) and the Govt. of India, Ministry of Earth Sciences (grant no. MoES/16/19/2017/APHH, DelhiFlux).

Review statement

This paper was edited by Joshua Fu and reviewed by Keding Lu and one anonymous referee.


Atkinson, R. and Arey, J.: Atmospheric Degradation of Volatile Organic Compounds Atmospheric Degradation of Volatile Organic Compounds, Chem. Rev., 103, 4605–4638,, 2003. 

Atkinson, R., Baulch, D. L., Cox, R. A., Crowley, J. N., Hampson, R. F., Hynes, R. G., Jenkin, M. E., Rossi, M. J., Troe, J., and IUPAC Subcommittee: Evaluated kinetic and photochemical data for atmospheric chemistry: Volume II – gas phase reactions of organic species, Atmos. Chem. Phys., 6, 3625–4055,, 2006 (data available at:, last access: August 2018). 

Avnery, S., Mauzerall, D. L., Liu, J., and Horowitz, L. W.: Global crop yield reductions due to surface ozone exposure: 2. Year 2030 potential crop production losses and economic damage under two scenarios of O3 pollution, Atmos. Environ., 45, 2297–2309,, 2011. 

Bansal, G. and Bandivadekar, A.: Overview of India's Vehicle Emissions Control Program Past Successes and Future Prospects, Int. Counc. Clean Transp., 1–180, available at: (last access: 1 September 2021), 2013. 

Beddows, D. C. S., Dall'osto, M., and Harrison, R. M.: An Enhanced Procedure for the Merging of Atmospheric Particle Size Distribution Data Measured Using Electrical Mobility and Time-of-Flifht Analysers, Aerosol Sci. Tech., 44, 930–938,, 2010. 

Calvert, J. G., Orlando, J. J., Stockwell, W. R., and Wallington, T. J.: The Mechanisms of Reactions Influencing Atmospheric Ozone, Oxford University Press, New York, 2015. 

Cash, J. M., Nemitz, E., Di Marco, C., Mullinger, N. J., Heal, M. R., Acton, W. J. F., Hewitt, C. N., Misztal, P. K., Mandal, T. K., Shivani, Gadi, R., Gurjar, B. R., and Langford, B.: Source apportionment of VOC fluxes in Old Delhi using high resolution proton transfer reaction mass spectrometry, in preparation, 2021. 

Castro, T., Madronich, S., Rivale, S., Muhlia, A., and Mar, B.: The influence of aerosols on photochemical smog in Mexico City, Atmos. Environ., 35, 1765–1772,, 2001. 

Chen, Y., Beig, G., Archer-Nicholls, S., Drysdale, W., Acton, J., Lowe, D., Nelson, B. S., Lee, J. D., Ran, L., Wang, Y., Wu, Z., Sahu, S. K., Sokhi, R. S., Singh, V., Gadi, R., Hewitt, C. N., Nemitz, E., Archibald, A., McFiggins, G., and Wild, O.: Avoiding high ozone pollution in Delhi, India, Faraday Discuss., 226, 502–514,, 2021. 

Crilley, L. R, Acton, J. W., Alam, M. S., Sommariva, R., Kramer, L. J., Drysdale, W., Lee, J., Nemitz E., Langford, B., Hewitt, C. N., and Bloss, W. J.: Measurements of nitrous acid in Delhi: insights into sources, in preparation, 2021. 

Cooper, O. R., Schultz, M. G., Schröder, S., Chang, K-L., Gaudel, A., Carbajal Benítez, G., Cuevas, E., Fröhlich, M., Galbally, I. E., Molloy, S., Kubistin, D., Lu, X., McClure-Begley, A., Nédélec, P., O'Brien, J., Oltmans, S. J., Petropavlovskikh, I., Ries, L., Senik, I., Sjöberg, K., Solberg, S., Spain, G. T., Spangl, W., Steinbacher, M., Tarasick, D., Thouret, V., and Xu, X.: Multi-decadal surface ozone trends at globally distributed remote locations, Elem. Sci. Anth., 8, 23,, 2020. 

Dickerson, R. R., Kondragunta, S., Stenchikov, G., Civerolo, K. L., Doddridge, B. G., and Holben, B. N.: The impact of aerosols on solar ultraviolet radiation and photochemical smog, Science, 278, 827–830,, 1997. 

Dube, R.: Revised national air quality standards, The Gazette of India, II.III, 1–4, 2009. 

Dunmore, R. E., Hopkins, J. R., Lidster, R. T., Lee, J. D., Evans, M. J., Rickard, A. R., Lewis, A. C., and Hamilton, J. F.: Diesel-related hydrocarbons can dominate gas phase reactive carbon in megacities, Atmos. Chem. Phys., 15, 9983–9996,, 2015. 

Elshorbany, Y. F., Kleffmann, J., Kurtenbach, R., Rubio, M., Lissi, E., Villena, G., Gramsch, E., Rickard, A. R., Pilling, M. J., and Wiesen, P.: Summertime photochemical ozone formation in Santiago, Chile, Atmos. Environ., 43, 6398–6407,, 2009. 

EPA (Center for Public Health and Environmental Assessment, Office of Research and Development, Environmental Protection Agency): Integrated Science Assessment (ISA) for Ozone and Related Photochemical Oxidants, U.S. Environmental Protection Agency, Research Triangle Park, North Carolina, USA, available at: (last access: September 2021), 2020 

Fleming, Z. L., Doherty, R. M., Von Schneidemesser, E., Malley, C. S., Cooper, O. R., Pinto, J. P., Colette, A., Xu, X., Simpson, D., Schultz, M. G., Lefohn, A. S., Hamad, S., Moolla, R., Solberg, S., and Feng, Z.: Tropospheric ozone assessment report: Present-day ozone distribution and trends relevant to human health, Elem. Sci. Anth., 6, 12,, 2018. 

Font, A., Guiseppin, L., Blangiardo, M., Ghersi, V., and Ghersi Fuller, G. W.: A tale of two cities: is air pollution improving in Paris and London?, Environ. Pollut., 249, 1–12,, 2019. 

Gaudel, A., Cooper, O. R., Ancellet, G., Barret, B., Boynard, A., Burrows, J. P., Clerbaux, C., Coheur, P. F., Cuesta, J., Cuevas, E., Doniki, S., Dufour, G., Ebojie, F., Foret, G., Garcia, O., Granados-Muñoz, M. J., Hannigan, J. W., Hase, F., Hassler, B., Huang, G., Hurtmans, D., Jaffe, D., Jones, N., Kalabokas, P., Kerridge, B., Kulawik, S., Latter, B., Leblanc, T., Le Flochmoën, E., Lin, W., Liu, J., Liu, X., Mahieu, E., McClure-Begley, A., Neu, J. L., Osman, M., Palm, M., Petetin, H., Petropavlovskikh, I., Querel, R., Rahpoe, N., Rozanov, A., Schultz, M. G., Schwab, J., Siddans, R., Smale, D., Steinbacher, M., Tanimoto, H., Tarasick, D. W., Thouret, V., Thompson, A. M., Trickl, T., Weatherhead, E., Wespes, C., Worden, H. M., Vigouroux, C., Xu, X., Zeng, G., and Ziemke, J.: Tropospheric Ozone Assessment Report: Present-day distribution and trends of tropospheric ozone relevant to climate and global atmospheric chemistry model evaluation, Elem. Sci. Anth., 6, 39,, 2018. 

Gaudel, A., Cooper, O. R., Chang, K-L., Bourgeois, I., Ziemke, J. R., Stode, S. A., Oman, L. D., Sellitto, P., Nédélec, P., Blot, R., Thouret, V., and Granier, C.: Aircraft observations since the 1990s reveal increases of tropospheric ozone at multiple locations across the Northern Hemisphere, Sci. Adv., 6, eaba8272,, 2020. 

Geng, F., Tie, X., Xu, J., Zhou, G., Peng, L., Gao, W., Tang, X., and Zhao, C.: Characterizations of ozone, NOx, and VOCs measured in Shanghai, China, Atmos. Environ., 42, 6873–6883,, 2008. 

George, C., Ammann, M., D'Anna, B., Donaldson, D. J., and Nizkorodov, S. A.: Heterogeneous Photochemistry in the Atmosphere, Chem. Rev., 115, 4218–4258,, 2015. 

Georgoulias, A. K., van der A, R. J., Stammes, P., Boersma, K. F., and Eskes, H. J.: Trends and trend reversal detection in 2 decades of tropospheric NO2 satellite observations, Atmos. Chem. Phys., 19, 6269–6294,, 2019. 

Gkatzelis, G. I., Coggon, M. M., Mcdonald, B. C., Peischl, J., Gilman, J. B., Aikin, K. C., Robinson, M. A., Canonaco, F., Prevot, A. S. H., and Trainer, M.: Observations Confirm that Volatile Chemical Products Are a Major Source of Petrochemical Emissions in U.S. Cities, Environ. Sci. Technol., 55, 4332–4343,, 2021. 

Gunthe, S. S., Liu, P., Panda, U., Raj, S. S., Sharma, A., Darbyshire, E., Reyes-Villegas, E., Allan, J., Chen, Y., Wang, X., Song, S., Pöhlker, M. L., Shi, L., Wang, Y., Kommula, S. M., Liu, T., Ravikrishna, R., McFiggans, G., Mickley, L. J., Martin, S. T., Pöschl, U., Andreae, M. O., and Coe, H.: Enhanced aerosol particle growth sustained by high continental chlorine emission in India, Nat. Geosci., 14, 77–84,, 2021. 

Guttikunda, S. K. and Calori, G.: A GIS based emissions inventory at 1 km × 1 km spatial resolution for air pollution analysis in Delhi, India, Atmos. Environ., 67, 101–111,, 2013. 

He, S. and Carmichael, G. R.: Sensitivity of photolysis rates and ozone production in the troposphere to aerosol properties, J. Geophys. Res.-Atmos., 104, 26307–26324,, 1999. 

Heland, J., Kleffmann, J., Kurtenbach, R., and Wiesen, P.: A new instrument to measure gaseous nitrous acid (HONO) in the atmosphere, Environ. Sci. Technol., 35, 3207–3212,, 2001. 

Hollaway, M., Wild, O., Yang, T., Sun, Y., Xu, W., Xie, C., Whalley, L., Slater, E., Heard, D., and Liu, D.: Photochemical impacts of haze pollution in an urban environment, Atmos. Chem. Phys., 19, 9699–9714,, 2019. 

Hopkins, J. R., Lewis, A. C., and Read, K. A.: A two-column method for long-term monitoring of non-methane hydrocarbons (NMHCs) and oxygenated volatile organic compounds (o-VOCs), J. Environ. Monit., 5, 8–13,, 2003. 

Jacob, D. J.: Heterogeneous chemistry and tropospheric ozone, Atmos. Environ., 34, 2131–2159,, 2000. 

Jain, N., Bhatia, A., and Pathak, H.: Emission of air pollutants from crop residue burning in India, Aerosol Air Qual. Res., 14, 422–430,, 2014. 

Jain, S. and Sharma, T.: Social and travel lockdown impact considering coronavirus disease (Covid-19) on air quality in megacities of india: Present benefits, future challenges and way forward, Aerosol Air Qual. Res., 20, 1222–1236,, 2020. 

Jain, S. L., Arya, B. C., Kumar, A., Ghude, S. D., and Kulkarni, P. S.: Observational study of surface ozone at New Delhi, India, Int. J. Remote Sens., 26, 3515–3524,, 2005. 

Jat, R. and Gurjar, B. R.: Contribution of different source sectors and source regions of Indo-Gangetic Plain in India to PM2.5 pollution and its short-term health impacts during peak polluted winter, Atmos. Pollut. Res., 12, 89–100,, 2021. 

Jenkin, M. E., Saunders, S. M., and Pilling, M. J.: The tropospheric degradation of volatile organic compounds: A protocol for mechanism development, Atmos. Environ., 31, 81–104,, 1997. 

Jenkin, M. E., Young, J. C., and Rickard, A. R.: The MCM v3.3.1 degradation scheme for isoprene, Atmos. Chem. Phys., 15, 11433–11459,, 2015. 

Jenkin, M. E., Derwent, R. G., and Wallington, T. J.: Photochemical ozone creation potentials for volatile organic compounds: Rationalization and estimation, Atmos. Environ., 163, 128–137,, 2017. 

Jerrett, M., Burnett, R. T., Pope, C. A., Ito, K., Thurston, G., Krewski, D., Shi, Y., Calle, E., and Thun, M.: Long-Term Ozone Exposure and Mortality, N. Engl. J. Med., 360, 1085–1095,, 2009. 

Jordan, A., Haidacher, S., Hanel, G., Hartungen, E., Märk, L., Seehauser, H., Schottkowsky, R., Sulzer, P., and Märk, T. D.: A high resolution and high sensitivity proton-transfer-reaction time-of-flight mass spectrometer (PTR-TOF-MS), Int. J. Mass Spectrom, 286, 122–128,, 2009. 

Kim, S., Sanchez, D., Wang, M., Seco, R., Jeong, D., Hughes, S., Barletta, B., Blake, D. R., Jung, J., Kim, D., Lee, G., Lee, M., Ahn, J., Lee, S. D., Cho, G., Sung, M. Y., Lee, Y. H., Kim, D. B., Kim, Y., Woo, J. H., Jo, D., Park, R., Park, J. H., Hong, Y. D., and Hong, J. H.: OH reactivity in urban and suburban regions in Seoul, South Korea-an East Asian megacity in a rapid transition, Faraday Discuss., 189, 231–251,, 2016. 

Kleffmann, J. and Wiesen, P.: Technical Note: Quantification of interferences of wet chemical HONO LOPAP measurements under simulated polar conditions, Atmos. Chem. Phys., 8, 6813–6822,, 2008. 

Klein, F., Farren, N. J., Bozzetti, C., Daellenbach, K. R., Kilic, D., Kumar, N. K., Pieber, S. M., Solwik, J. G., Tuthil, R. N., Hamilton, J. F., Baltensperger U., Prévôt, A. S. H., and El Haddad, I.: Indoor terpene emissions from cooking with herbs and pepper and their secondary organic aeosol production potential, Sci. Rep.-UK, 6, 36623,, 2016. 

Koss, A. R., Sekimoto, K., Gilman, J. B., Selimovic, V., Coggon, M. M., Zarzana, K. J., Yuan, B., Lerner, B. M., Brown, S. S., Jimenez, J. L., Krechmer, J., Roberts, J. M., Warneke, C., Yokelson, R. J., and de Gouw, J.: Non-methane organic gas emissions from biomass burning: identification, quantification, and emission factors from PTR-ToF during the FIREX 2016 laboratory experiment, Atmos. Chem. Phys., 18, 3299–3319,, 2018. 

Krotkov, N. A., McLinden, C. A., Li, C., Lamsal, L. N., Celarier, E. A., Marchenko, S. V., Swartz, W. H., Bucsela, E. J., Joiner, J., Duncan, B. N., Boersma, K. F., Veefkind, J. P., Levelt, P. F., Fioletov, V. E., Dickerson, R. R., He, H., Lu, Z., and Streets, D. G.: Aura OMI observations of regional SO2 and NO2 pollution changes from 2005 to 2015, Atmos. Chem. Phys., 16, 4605–4629,, 2016. 

Kulkarni, S. H., Ghude, S. D., Jena, C., Karumuri, R. K., Sinha, B., Sinha, V., Kumar, R., Soni, V. K., and Khare, M.: How Much Does Large-Scale Crop Residue Burning Affect the Air Quality in Delhi?, Environ. Sci. Technol., 54, 4790–4799,, 2020. 

Kumar, P., Gulia, S., Harrison, R. M., and Khare, M.: The influence of odd–even car trial on fine and coarse particles in Delhi, Environ. Pollut., 225, 20–30,, 2017. 

Lee, J. D., Whalley, L. K., Heard, D. E., Stone, D., Dunmore, R. E., Hamilton, J. F., Young, D. E., Allan, J. D., Laufs, S., and Kleffmann, J.: Detailed budget analysis of HONO in central London reveals a missing daytime source, Atmos. Chem. Phys., 16, 2747–2764,, 2016. 

Li, K., Jacob, D. J., Liao, H., Shen, L., Zhang, Q., and Bates, K. H.: Anthropogenic drivers of 2013–2017 trends in summer surface ozone in China, P. Natl. Acad. Sci. USA, 116, 422–427,, 2019. 

Li, K., Jacob, D. J., Liao, H., Qiu, Y., Shen, L., Zhai, S., Bates, K. H., Sulprizio, M. P., Song, S., Lu, X., Zhang, Q., Zheng, B., Zhang, Y., Zhang, J., Lee, H. C., and Kuk, S. K.: Ozone pollution in the North China Plain spreading into the late-winter haze season, P. Natl. Acad. Sci. USA, 118, 1–7,, 2021. 

Li, M., Zhang, Q., Zheng, B., Tong, D., Lei, Y., Liu, F., Hong, C., Kang, S., Yan, L., Zhang, Y., Bo, Y., Su, H., Cheng, Y., and He, K.: Persistent growth of anthropogenic non-methane volatile organic compound (NMVOC) emissions in China during 1990–2017: drivers, speciation and ozone formation potential, Atmos. Chem. Phys., 19, 8897–8913,, 2019. 

Liao, Z., Ling, Z., Gao, M., Sun, J., Zhao, W., Ma, P., Quan, J., and Fan, S.: Tropospheric Ozone Variability Over Hong Kong Based on Recent 20 years (2000–2019) Ozonesonde Observation, J. Geophys. Res.-Atmos., 126, e2020JD033054,, 2021. 

Liu, F., Zhang, Q., Van Der A, R. J., Zheng, B., Tong, D., Yan, L., Zheng, Y., and He, K.: Recent reduction in NOx emissions over China: Synthesis of satellite observations and emission inventories, Environ. Res. Lett., 11, 114002,, 2016. 

Liu, Z., Yang, Y., Costabile, F., Amoroso, A., Zhao, C., Huey, L. G., Stickel, R., Liao, J., and Zhu, T.: Evidence of aerosols as a media for rapid daytime HONO production over China, Environ. Sci. Technol., 48, 14386–14391, 2014. 

Lou, S., Holland, F., Rohrer, F., Lu, K., Bohn, B., Brauers, T., Chang, C. C., Fuchs, H., Häseler, R., Kita, K., Kondo, Y., Li, X., Shao, M., Zeng, L., Wahner, A., Zhang, Y., Wang, W., and Hofzumahaus, A.: Atmospheric OH reactivities in the Pearl River Delta – China in summer 2006: measurement and model results, Atmos. Chem. Phys., 10, 11243–11260,, 2010. 

Lu, X., Hong, J., Zhang, L., Cooper, O. R., Schultz, M. G., Xu, X., Wang, T., Gao, M., Zhao, Y., and Zhang, Y.: Severe Surface Ozone Pollution in China: A Global Perspective, Environ. Sci. Technol. Lett., 5, 487–494,, 2018a. 

Lu, X., Wang, Y., Li, J., Shen, L., and Fung, J. C. H.: Evidence of herterogenous HONO formation from aerosols and the regional photochemical impacts of this HONO source, Environ. Res. Lett., 13, 114002,, 2018b. 

Lu, X., Zhang, L., Wang, X., Gao, M., Li, K., Zhang, Y., Yue, X., and Zhang, Y.: Rapid increases in warm-season surface ozone and resulting health impact in China since 2013, Environ. Sci. Technol. Lett., 7, 240–247, 2020. 

Mahato, S., Pal, S., and Ghosh, K. G.: Effect of lockdown amid COVID-19 pandemic on air quality of the megacity Delhi, India, Sci. Total Environ., 730, 139086,, 2020. 

McDonald, B. C., De Gouw, J. A., Gilman, J. B., Jathar, S. H., Akherati, A., Cappa, C. D., Jimenez, J. L., Lee-Taylor, J., Hayes, P. L., McKeen, S. A., and Cui, Y. Y.: Volatile chemical products emerging as largest petrochemical source of urban organic emissions, Science, 359, 760–764, 2018. 

Miyazaki, K., Eskes, H., Sudo, K., Boersma, K. F., Bowman, K., and Kanaya, Y.: Decadal changes in global surface NOx emissions from multi-constituent satellite data assimilation, Atmos. Chem. Phys., 17, 807–837,, 2017. 

Molina, L. T.: Introductory Lecture: Air Quality in Megacities, Faraday Discuss., 226, 9–52,, 2021. 

Nemitz, E., Acton, W. J., Alam, M. S., Drysdale, W. S., Dunmore, R. E., Hamilton, J. F., Hopkins, J. R., Langford, B., Nelson, B. S., Stewart, G. S., Vaughan, A. R., and Whalley, L. K.: (APHH India) Megacity Delhi atmospheric emission quantification, assessment and impacts (DelhiFlux), CEDA Archive [data set],, last access: March 2020. 

Parrish, D. D., Lamarque, J. F., Naik, V., Horowitz, L., Shindell, D. T., Staehelin, J., Derwent, R., Cooper, O. R., Tanimoto, H., Volz-Thomas, A., Gilge, S., Scheel, H. E., Steinbacher, M., and Fröhlich, M.: Long-term changes in lower tropospheric baseline ozone concentrations: Comparing chemistry-climate models and observations at northern midlatitudes, J. Geophys. Res., 119, 5719–5736,, 2014. 

Pollack, I. B., Ryerson, T. B., Trainer, M., Neuman, J. A., Roberts, J. M., and Parrish, D. D.: Trends in ozone, its precursors, and related secondary oxidation products in Los Angeles, California: A synthesis of measurements from 1960 to 2010, J. Geophys. Res.-Atmos., 118, 5893-5911,, 2013. 

Ran, L., Zhao, C. S., Xu, W. Y., Lu, X. Q., Han, M., Lin, W. L., Yan, P., Xu, X. B., Deng, Z. Z., Ma, N., Liu, P. F., Yu, J., Liang, W. D., and Chen, L. L.: VOC reactivity and its effect on ozone production during the HaChi summer campaign, Atmos. Chem. Phys., 11, 4657–4667,, 2011. 

Real, E. and Sartelet, K.: Modeling of photolysis rates over Europe: impact on chemical gaseous species and aerosols, Atmos. Chem. Phys., 11, 1711–1727,, 2011. 

Ren, X., Brune, W. H., Mao, J., Mitchell, M. J., Lesher, R. L., Simpas, J. B., Metcalf, A. R., Schwab, J. J., Cai, C., Li, Y., Demerjian, K. L., Felton, H. D., Boynton, G., Adams, A., Perry, J., He, Y., Zhou, X., and Hou, J.: Behavior of OH and HO2 in the winter atmosphere in New York City, Atmos. Environ., 40, 252–263,, 2006. 

Royal Society: Ground-level ozone in the 21st century: future trends, impacts and policy implications, available at: (last access: July 2021), 2008. 

Saunders, S. M., Jenkin, M. E., Derwent, R. G., and Pilling, M. J.: Protocol for the development of the Master Chemical Mechanism, MCM v3 (Part A): tropospheric degradation of non-aromatic volatile organic compounds, Atmos. Chem. Phys., 3, 161–180,, 2003. 

Sawlani, R., Agnihotri, R., Sharma, C., Patra, P. K., Dimri, A. P., Ram, K., and Verma, R. L.: The severe Delhi SMOG of 2016: A case of delayed crop residue burning, coincident firecracker emissions, and atypical meteorology, Atmos. Pollut. Res., 10, 868–879,, 2019. 

Schelegle, E. S., Morales, C. A., Walby, W. F., Marion, S., and Allen, R. P.: 6.6-hour inhalation of ozone concentration from 60 to 87 parts per million in healthy humans, Am. J. Resp. Crit. Care, 180,, 2009. 

Sharma, A., Sharma, S. K., Rohtash, and Mandal, T. K.: Influence of ozone precursors and particulate matter on the variation of surface ozone at an urban site of Delhi, India, Sustain. Environ. Res., 26, 76–83,, 2016. 

Sharma, S., Zhang, M., Anshika, Gao, J., Zhang, H., and Kota, S. H.: Effect of restricted emissions during COVID-19 on air quality in India, Sci. Total Environ., 728, 138878,, 2020. 

Shi, Z., Song, C., Liu, B., Lu, G., Xu, J., Van Vu, T., Elliott, R. J. R., Li, W., Bloss, W. J., and Harrison, R. M.: Abrupt but smaller than expected changes in surface air quality attributable to COVID-19 lockdowns, Sci. Adv., 7, eabd6696,, 2021. 

Shirley, T. R., Brune, W. H., Ren, X., Mao, J., Lesher, R., Cardenas, B., Volkamer, R., Molina, L. T., Molina, M. J., Lamb, B., Velasco, E., Jobson, T., and Alexander, M.: Atmospheric oxidation in the Mexico City Metropolitan Area (MCMA) during April 2003, Atmos. Chem. Phys., 6, 2753–2765,, 2006. 

Shivani, Gadi, R., Sharma, S. K., Mandal, T. K., Kumar, R., Mona, S., Kumar, S., and Kumar, S.: Levels and sources of organic compounds in fine ambient aerosols over National Capital Region of India, Environ. Sci. Pollut. Res., 25, 31071–31090, 2018. 

Sillman, S.: The relation between ozone, NOx and hydrocarbons in urban and polluted rural environments, Atmos. Environ., 33, 1821–1845,, 1999. 

Sillman, S., Logan, A., and Wofsy, C.: The Sensitivity of Ozone to Nitrogen Oxides and Hydrocarbons in Regional Ozone Episodes, J. Geophys. Res., 95, 1837–1851, 1990. 

Simpson, I. J., Blake, D. R., Blake, N. J., Meinardi, S., Barletta, B., Hughes, S. C., Fleming, L. T., Crawford, J. H., Diskin, G. S., Emmons, L. K., Fried, A., Guo, H., Peterson, D. A., Wisthaler, A., Woo, J.-H., Barré, J., Gaubert, B., Kim, J., Kim, M. J., Kim, Y., Knote, C., Mikoviny, T., Pusede, S. E., Schroeder, J. R., Wang, Y., Wennberg, P. O., and Zeng, L.: Characterization, sources and reactivity of volatile organic compounds (VOCs) in Seoul and surrounding regions during KORUS-AQ, Elem. Sci. Anth., 8, 37,, 2020. 

Skeie, R. B., Myhre, G., Hodnebrog, Ø., Cameron-Smith, P. J., Deushi, M., Hegglin, M. I., Horowitz, L. W., Kramer, R. J., Michou, M., Mills, M. J., Olivié, D. J. K., O'Connor, F. M., Payter, D., Samset, B. H., Sellar, A., Shindell, D., Takemura, T., Tilmes, S., and Wu, T.: Historical total ozone radiative forcing dervied from CMIP6 simulations, npj Clim. Atmos. Sci., 3, 32,, 2020. 

Sommariva, R., Cox, S., Martin, C., Borońska, K., Young, J., Jimack, P. K., Pilling, M. J., Matthaios, V. N., Nelson, B. S., Newland, M. J., Panagi, M., Bloss, W. J., Monks, P. S., and Rickard, A. R.: AtChem (version 1), an open-source box model for the Master Chemical Mechanism, Geosci. Model Dev., 13, 169–183,, 2020. 

Squires, F. A., Nemitz, E., Langford, B., Wild, O., Drysdale, W. S., Acton, W. J. F., Fu, P., Grimmond, C. S. B., Hamilton, J. F., Hewitt, C. N., Hollaway, M., Kotthaus, S., Lee, J., Metzger, S., Pingintha-Durden, N., Shaw, M., Vaughan, A. R., Wang, X., Wu, R., Zhang, Q., and Zhang, Y.: Measurements of traffic-dominated pollutant emissions in a Chinese megacity, Atmos. Chem. Phys., 20, 8737–8761,, 2020. 

Stewart, G. J., Nelson, B. S., Drysdale, W. S., Acton, W. J. F., Vaughan, A. R., Hopkins, J. R., Dunmore, R. E., Hewitt, C. N., Nemitz, E., Mullinger, N., Langford, B., Shivani, Reyes-Villegas, E., Gadi, R., Rickard, A. R., Lee, J. D., and Hamilton, J. F.: Sources of non-methane hydrocarbons in surface air in Delhi, India, Faraday Discuss., 409–431,, 2021a. 

Stewart, G. J., Nelson, B. S., Acton, W. J. F., Vaughan, A. R., Farren, N. J., Hopkins, J. R., Ward, M. W., Swift, S. J., Arya, R., Mondal, A., Jangirh, R., Ahlawat, S., Yadav, L., Sharma, S. K., Yunus, S. S. M., Hewitt, C. N., Nemitz, E., Mullinger, N., Gadi, R., Sahu, L. K., Tripathi, N., Rickard, A. R., Lee, J. D., Mandal, T. K., and Hamilton, J. F.: Emissions of intermediate-volatility and semi-volatile organic compounds from domestic fuels used in Delhi, India, Atmos. Chem. Phys., 21, 2407–2426,, 2021b. 

Stewart, G. J., Acton, W. J. F., Nelson, B. S., Vaughan, A. R., Hopkins, J. R., Arya, R., Mondal, A., Jangirh, R., Ahlawat, S., Yadav, L., Sharma, S. K., Dunmore, R. E., Yunus, S. S. M., Hewitt, C. N., Nemitz, E., Mullinger, N., Gadi, R., Sahu, L. K., Tripathi, N., Rickard, A. R., Lee, J. D., Mandal, T. K., and Hamilton, J. F.: Emissions of non-methane volatile organic compounds from combustion of domestic fuels in Delhi, India, Atmos. Chem. Phys., 21, 2383–2406,, 2021c. 

Stewart, G. J., Nelson, B. S., Acton, W. J. F., Vaughan, A. R., Hopkins, J. R., Yunus, S. S. M., Hewitt, C. N., Wild, O., Nemitz, E. G., Gadi, R., Sahu, L. K., Mandal, T. K., Gurjar, B. R., Rickard, A. R., Lee, J. D., and Hamilton, J. F.: Emission estimates and inventories of non-methane organic compounds from anthropogenic burning sources in India, Atmos. Environ.-X, 11, 100115,, 2021d. 

Stewart, G. J., Nelson, B. S., Acton, W. J. F., Vaughan, A. R., Hopkins, J. R., Yunus, S. S. M., Hewitt, C. N., Nemitz, E., Mullinger, N., Gadi, R., Rickard, A. R., Lee, J. D., Mandal, T. K., and Hamilton, J. F.: Comprehensive organic emission profiles, secondary organic aerosol production potential, and OH reactivity of domestic fuel combustion in Delhi, India, Environ. Sci. Atmos., 1, 104–117,, 2021e. 

Tan, Z., Lu, K., Jiang, M., Su, R., Wang, H., Lou, S., Fu, Q., Zhai, C., Tan, Q., Yue, D., Chen, D., Wang, Z., Xie, S., Zeng, L., and Zhang, Y.: Daytime atmospheric oxidation capacity in four Chinese megacities during the photochemically polluted season: a case study based on box model simulation, Atmos. Chem. Phys., 19, 3493–3513,, 2019. 

Tan, Z., Hofzumahaus, A., Lu, K., Brown, S. S., Holland, F., Huey, L. G., Kiendler-Scharr, A., Li, X., Liu, X., Ma, N., Min, K. E., Rohrer, F., Shao, M., Wahner, A., Wang, Y., Wiedensohler, A., Wu, Y., Wu, Z., Zeng, L., Zhang, Y., and Fuchs, H.: No Evidence for a Significant Impact of Heterogeneous Chemistry on Radical Concentrations in the North China Plain in Summer 2014, Environ. Sci. Technol., 54, 5973–5979,, 2020. 

Tarasick, D., Galbally, I. E., Cooper, O. R., Schultz, M. G., Ancellet, G., Leblanc, T., Wallington, T. J., Ziemke, J., Liu, X., Steinbacher, M., Staehelin, J., Vigouroux, C., Hannigan, J. W., García, O., Foret, G., Zanis, P., Weatherhead, E., Petropavlovskikh, I., Worden, H., Osman, M., Liu, J., Chang, K.-L., Gaudel, A., Lin, M., Granados-Muñoz, M., Thompson, A. M., Oltmans, S. J., Cuesta, J., Dufour, G., Thouret, V., Hassler, B., Trickl, T., and Neu, J. L.: Tropospheric Ozone Assessment Report: Tropospheric ozone from 1877 to 2016, observed levels, trends and uncertainties, Elem. Sci. Anth., 7, 39,, 2019. 

Thurston, G. D., Lippmann, M., Scott, M. B., and Fine, J. M.: Summertime haze air pollution and children with asthma, Am. J. Resp. Crit. Care, 155, 654–660,, 1997. 

Tong, S., Hou, S., Zhang, Y., Chu, B., Liu, Y., He, H., Zhao, P., and Ge, M.: Exploring the nitrous acid (HONO) formation mechanims in winter Beijing: direct emissions and heterogenous production in urban and suburban areas, Faraday Discuss., 189, 213–230, 2016. 

United Nations, Department of Economic and Social Affairs, Population Division: The World's Cities in 2018 – Data booklet (ST/ESA/ SER.A/417), 34 pp., available at: (last access: September 2021), 2018. 

United Nations, Department of Economic and Social Affairs, Population Division: World Urbanization Prospects: The 2018 Revision (ST/ESA/SER.A/420), available at: (last access: September 2021), 2019. 

Wang, L., Slowik, J. G., Tripathi, N., Bhattu, D., Rai, P., Kumar, V., Vats, P., Satish, R., Baltensperger, U., Ganguly, D., Rastogi, N., Sahu, L. K., Tripathi, S. N., and Prévôt, A. S. H.: Source characterization of volatile organic compounds measured by proton-transfer-reaction time-of-flight mass spectrometers in Delhi, India, Atmos. Chem. Phys., 20, 9753–9770,, 2020. 

Wang, W., Li, X., Shao, M., Hu, M., Zeng, L., Wu, Y., and Tan, T.: The impact of aerosols on photolysis frequencies and ozone production in Beijing during the 4-year period 2012–2015, Atmos. Chem. Phys., 19, 9413–9429,, 2019. 

Whalley, L. K., Stone, D., Bandy, B., Dunmore, R., Hamilton, J. F., Hopkins, J., Lee, J. D., Lewis, A. C., and Heard, D. E.: Atmospheric OH reactivity in central London: observations, model predictions and estimates of in situ ozone production, Atmos. Chem. Phys., 16, 2109–2122,, 2016. 

Whalley, L. K., Stone, D., Dunmore, R., Hamilton, J., Hopkins, J. R., Lee, J. D., Lewis, A. C., Williams, P., Kleffmann, J., Laufs, S., Woodward-Massey, R., and Heard, D. E.: Understanding in situ ozone production in the summertime through radical observations and modelling studies during the Clean air for London project (ClearfLo), Atmos. Chem. Phys., 18, 2547–2571,, 2018. 

Whalley, L. K., Slater, E. J., Woodward-Massey, R., Ye, C., Lee, J. D., Squires, F., Hopkins, J. R., Dunmore, R. E., Shaw, M., Hamilton, J. F., Lewis, A. C., Mehra, A., Worrall, S. D., Bacak, A., Bannan, T. J., Coe, H., Percival, C. J., Ouyang, B., Jones, R. L., Crilley, L. R., Kramer, L. J., Bloss, W. J., Vu, T., Kotthaus, S., Grimmond, S., Sun, Y., Xu, W., Yue, S., Ren, L., Acton, W. J. F., Hewitt, C. N., Wang, X., Fu, P., and Heard, D. E.: Evaluating the sensitivity of radical chemistry and ozone formation to ambient VOCs and NOx in Beijing, Atmos. Chem. Phys., 21, 2125–2147,, 2021. 

WHO (World Health Organization Regional Office for Europe): Air quality guidelines global update 2005: Particulate matter, ozone, nitrogen dioxide and sulfur dioxide, available at: (last access: March 2020), 2005. 

WHO (World Health Organization): Review of evidence on health aspects of air pollution – REVIHAAP final technical report, available at: (last access: March 2020), 2013. 

WHO (World Health Organization): Global Urban Ambient Air Pollution Database, available at: (last access: 31 July 2019), 2014. 

Yoshino, A., Sadanaga, Y., Watanabe, K., Kato, S., Miyakawa, Y., Matsumoto, J., and Kajii, Y.: Measurement of total OH reactivity by laser-induced pump and probe technique-comprehensive observations in the urban atmosphere of Tokyo, Atmos. Environ., 40, 7869–7881,, 2006. 

Zavala, M., Brune, W. H., Velasco, E., Retama, A., Cruz-Alavez, L. A., and Molina, L. T.: Changes in ozone production and VOC reactivity in the atmosphere of the Mexico City Metropolitan Area, Atmos. Environ., 238, 117747,, 2020.  

Zhang, Q., Zheng, Y., Tong, D., Shao, M., Wang, S., Zhang, Y., Xu, X., Wang, J., He, H., Liu, W., Ding, Y., Lei, Y., Li, J., Wang, Z., Zhang, X., Wang, Y., Cheng, J., Liu, Y., Shi, Q., Yan, L., Geng, G., Hong, C., Li, M., Liu, F., Zheng, B., Cao, J., Ding, A., Gao, J., Fu, Q., Huo, J., Liu, B., Liu, Z., Yang, F., He, K., and Hao, J.: Drivers of improved PM2.5 air quality in China from 2013 to 2017, P. Natl. Acad. Sci. USA, 116, 24463–24469,, 2019. 

Ziemke, J. R., Oman, L. D., Strode, S. A., Douglass, A. R., Olsen, M. A., McPeters, R. D., Bhartia, P. K., Froidevaux, L., Labow, G. J., Witte, J. C., Thompson, A. M., Haffner, D. P., Kramarova, N. A., Frith, S. M., Huang, L.-K., Jaross, G. R., Seftor, C. J., Deland, M. T., and Taylor, S. L.: Trends in global tropospheric ozone inferred from a composite record of TOMS/OMI/MLS/OMPS satellite measurements and the MERRA-2 GMI simulation , Atmos. Chem. Phys., 19, 3257–3269,, 2019. 

Short summary
Ozone production at an urban site in Delhi is sensitive to volatile organic compound (VOC) concentrations, particularly those of the aromatic, monoterpene, and alkene VOC classes. The change in ozone production by varying atmospheric pollutants according to their sources, as defined in an emissions inventory, is investigated. The study suggests that reducing road transport emissions alone does not reduce reactive VOCs in the atmosphere enough to perturb an increase in ozone production.
Final-revised paper