Articles | Volume 22, issue 7
Research article
12 Apr 2022
Research article |  | 12 Apr 2022

Influence of photochemical loss of volatile organic compounds on understanding ozone formation mechanism

Wei Ma, Zemin Feng, Junlei Zhan, Yongchun Liu, Pengfei Liu, Chengtang Liu, Qingxin Ma, Kang Yang, Yafei Wang, Hong He, Markku Kulmala, Yujing Mu, and Junfeng Liu

Volatile organic compounds (VOCs) tend to be consumed by atmospheric oxidants, resulting in substantial photochemical loss during transport. An observation-based model was used to evaluate the influence of photochemical loss of VOCs on the sensitivity regime and mechanisms of ozone formation. Our results showed that a VOC-limited regime based on observed VOC concentrations shifted to a transition regime with a photochemical initial concentration of VOCs (PIC-VOCs) in the morning. The net ozone formation rate was underestimated by 3 ppb h−1 (∼36 ppb d−1) based on the measured VOCs when compared with the PIC-VOCs. The relative contribution of the RO2 path to ozone production based on the PIC-VOCs accordingly increased by 13.4 %; in particular, the contribution of alkene-derived RO2 increased by approximately 10.2 %. In addition, the OH–HO2 radical cycle was obviously accelerated by highly reactive alkenes after accounting for photochemical loss of VOCs. The contribution of local photochemistry might be underestimated for both local and regional ozone pollution if consumed VOCs are not accounted for, and policymaking on ozone pollution prevention should focus on VOCs with a high reactivity.

1 Introduction

Ground surface ozone (O3) is an important atmospheric pollutant that is harmful to human health and is connected with respiratory, cardiovascular diseases, and premature mortality (Cohen et al., 2017). It is also harmful to vegetation growth. For example, it led to annual reductions in the yields of rice and wheat of 8 % and 6 %, respectively, and reduced forest biomass growth by 11 %–13 % in 2015 in China (Feng et al., 2019). Surface O3 concentrations have increased by 11.9 % over eastern China despite the air pollution control measures implemented in China from 2012 to 2017 (Dang and Liao, 2019). An economic loss of 0.09 % of Chinese gross domestic product (CNY 78 billion) is predicted for 2030 if policies against O3 pollution are not properly implemented (Xie et al., 2019). Therefore, urgent action to minimize O3 pollution in China is needed.

Tropospheric O3 is mainly produced from photochemical reactions between volatile organic compounds (VOCs) and nitrogen oxides (NOx: NO+NO2) (Seinfeld and Pandis, 2016; Liu et al., 2021). O3 is generated from a collision of O2 and O(3P) that is produced from photolysis of NO2 in the atmosphere. Peroxyl radicals (HO2 and RO2), which are produced from the oxidation of VOCs by the OH radical can efficiently convert NO (from the photolysis of NO2) to NO2, leading to a net O3 production by compensating for the titration of O3 by NO (Monks, 2005; Zhang et al., 2021a). Over the past 2 decades, a number of field observations focused on O3 pollution levels, and its precursors have been carried out in the Beijing–Tianjin–Hebei (BTH), Yangtze River Delta (YRD), and Pearl River Delta (PRD) regions (Wang et al., 2017; Li et al., 2019; Xue et al., 2014; Zhang et al., 2019). Due to the nonlinear relationship between O3 and its precursors and the variations in meteorological conditions, numerous studies have been performed to understand the sensitivity regime of O3 formation (Ling and Guo, 2014; Zhang et al., 2020), the photochemical process of O3 formation based on box models or observation-based models (OBMs) (He et al., 2019; Tan et al., 2019), and the sources of O3 using regional chemical transport models (Li et al., 2016b, c). Recently, the instantaneous production rate of the O3 formation process has attracted more attention; for example, studies examining radical recycling (OH–RO2–RO–HO2–OH) related to the production of O3 have been performed (Lu et al., 2017; Tan et al., 2017; Whalley et al., 2018). HCHO photolysis and alkene ozonolysis contributed approximately 85 % to the primary production of HO2 and HO radicals in Beijing, Shanghai and Guangzhou (Tan et al., 2019; Yang et al., 2017). The importance of HONO and HCHO photolysis for primary radical production has also been proposed in suburban and rural areas (Tan et al., 2017; Lu et al., 2012, 2013).

All of the OBM studies investigating the relationship between O3 and VOCs were based on measured datasets. However, VOCs are highly reactive to atmospheric oxidants, such as OH, NO3, and O3, among which OH is dominant. The lifetimes of some highly reactive VOCs, such as isoprene, are as short as only a few tens of minutes under typical daytime atmospheric conditions. The mixing ratios of VOCs observed at a sampling site are actually the residues of VOCs from emissions due to the photochemical loss during transport from the source site to the receptor site. If photochemically consumed VOCs are not considered, the O3 formation sensitivity and net O3 production may be misunderstood, and subsequent policymaking on O3 pollution prevention at regional or urban scales may be misguided. Thus, the photochemical age-based approach has been applied to evaluate the effect of photochemical processes on VOC measurements (Shao et al., 2011). This method was used to qualitatively or semi-quantitatively estimate the O3 formation process of the source-receptor (Gao et al., 2018) by calculating the O3 formation potential (OFP) (Han et al., 2017), identifying the critical species for O3 formation (Gao et al., 2021), or evaluating the VOC emissions ratio (Yuan et al., 2013). In evaluating the importance of initial VOCs to ozone production, Xie et al. (2008) found that the OFP at a Peking University site increased by 70 % after accounting for the photochemical loss of VOCs. Li et al. (2015) also showed that the OFPs of total NMHCs (excluding isoprene) increased by 16.1 % (from 59.6 to 69.2 ppb O3), 12.1 % (from 33.5 to 37.5 ppb O3), and 3.4 % (from 68.9 to 71.2 ppb O3) after correcting for photochemical loss in Gucheng, Quzhou, and Beijing, respectively. Gao et al. (2018) reported that the OFP could be underestimated by 23.4 % (62.4 ppb O3) in Beijing if the photochemical loss of VOCs is not considered. Zhan et al. (2021) found that based on measured VOCs, the OFP increased from 57.8 to 103.9 ppb using the initial VOCs. All the previous work was based on the maximum incremental reactivities (MIR) method. The application of such calculations using the MIR method is restricted to areas or episodes in which O3 formation is VOC sensitive (Carter, 1994). In the troposphere, the sensitivity of ozone formation to NOx and VOCs varies greatly, as evidenced by the wide range of OFP underestimations from ∼3 % to 70 % in previous work. In addition, the MIR values of VOC species for a specific region are calculated with the base scenario, in which NO concentration and other parameters are the values that correspond to the maximal incremental reactivity (IR). The fixed MIR values of different VOCs can neither reflect the nonlinear relationship between ozone and VOCs, involved in the complicated radical recycling (OH–RO2–RO–HO2–OH) related to the production of ozone, nor be used for analyzing the radical budget of the initial VOCs concentration. Thus, a quantitative analysis is necessary to explicitly understand the influence of photochemical loss of VOCs on ozone formation and its mechanism based on OBM studies, in which the dynamic atmospheric and meteorological conditions is accounted for.

In this study, an OBM was used to evaluate the local O3 formation process in summer in Beijing based on concentrations of observed and photochemical initial concentrations of VOCs (PIC-VOCs). The O3NOx–VOC sensitivity, instantaneous O3 formation rate, and in situ O3 formation process were discussed. The aim of this study was to understand the possible influence of photochemical loss of VOCs on the formation sensitivity regime of O3 and how the photochemical loss of VOCs affects O3 formation. This study can provide new insight for better understanding atmospheric O3 pollution.

2 Methodology

2.1 Experimental section

Field observations were carried out on the Qingyuan campus of the Beijing Institute of Petrochemical Technology (BIPT, 39.73 N and 116.33 E) (Fig. S1 in the Supplement). Details on the observation site have been described in our previous work (Zhan et al., 2021). In short, the site is a typical suburban site in the Daxing District between 5th Ring Road and 6th Ring Road. The field campaign was carried out during 1–28 August 2019, when photochemistry was the most active and rainfall was rare in Beijing.

The concentrations of non-methane hydrocarbons (NMHCs) were detected by both a gas chromatography–flame ionization detector (GC/FID) and a single photon ionization mass spectrometer (SPI-MS 3000, Guangzhou Hexin Instrument Co., Ltd., China). A detailed description of the instrumentation can be found in previous publications (Zhan et al., 2021; Chen et al., 2020). The SPI-MS was also used to detect halohydrocarbons. More details on this instrument and its parameter settings have been described in previous studies (Zhang et al., 2019; Liu et al., 2020a). In short, a 0.002 in. thick polydimethylsiloxane (PDMS) membrane (Technical Products Inc., USA) was used to collect VOCs and diffuse them from the sample site to the detector under high vacuum conditions. Vacuum ultraviolet (VUV) light generated by a commercial D2 lamp (Hamamatsu, Japan) was utilized for ionization at 10.8 eV. For ion detection, two microchannel plates (MCPs, Hamamatsu, Japan) assembled with a chevron-type configuration were employed. This TOF-MS has a limit of detection (LOD) varying from 50 to 1 ppb with a 1 min time resolution for most trace gases without any preconcentration procedure. To verify the data compatibility of the SPI-MS and GC/FID, we compared the concentrations of toluene measured using these two different instruments (Fig. S2). The correlation coefficient was 0.9 (with a slope of 0.7), indicating that the concentrations of NMHCs were comparable using these two measurement techniques.

Oxygenated VOCs (OVOCs) were collected using 2,4-dinitrophenylhydrazine (DNPH)-coated silica gel cartridges (Sep-Pak, Waters) by an automatic sampling device with a sampling flow rate of 1.2 L min−1 and a duration of 2 h for each sample. Following this, the OVOCs were analyzed using high-performance liquid chromatography (HPLC, Inertsil ODS-P 5 µm 4.6×250 mm column, GL Sciences) with an acetonitrile–water binary mobile phase (Ma et al., 2019). To avoid possible contamination or desorption after sampling, cartridges were capped, placed into tightly closed plastic bags and kept in a refrigerator before analysis. The sampled cartridges were eluted with 5 mL acetonitrile and analyzed by HPLC as soon as possible after they were shipped back to the laboratory. This system was calibrated with eight-gradient standard solutions (TO11/IP-6A Aldehyde/Ketone-DNPH Mix, SUPELCO). The correlation coefficients were all greater than 0.999. The LOD for most OVOCs was approximately 10 ppt.

Trace gases, including NOx, SO2, CO, and O3, were measured using corresponding analyzers (Thermo Scientific, 42i, 43i, 48i, and 49i, respectively). The HONO concentration was measured using a homemade long-path absorption photometer (LOPAP) (Liu et al., 2020c). The meteorological parameters, including temperature (T), pressure (P), relative humidity (RH), wind speed, and wind direction, were measured by a weather station (AWS310, Vaisala). The photolysis rate (JNO2) was measured via continuous measurement of the actinic flux in the wavelength range of 285–375 nm using a JNO2 filter radiometer (JNO2 radiometer, Metcon).

2.2 Calculation of photochemical loss of VOCs

The photochemical loss of VOCs was calculated using the ratio method (Wiedinmyer et al., 2001; Yuan et al., 2013). The initial mixing ratio of a specific VOC was calculated using the following equations (Mckeen et al., 1996):


where [VOCi]t and [VOCi]t0 are the observed and initial concentrations of VOCi, respectively; ki is the second-order reaction rate between compound i and OH radical; and [OH] and Δt are the concentration of OH radical and the photochemical aging time, respectively. kE and kX are rate constants for the reaction between OH radicals and ethylbenzene (7.00×10-12 cm3 molec.−1 s−1) and xylene (1.87×10-11 cm3 molec.−1 s−1) (Atkinson and Arey, 2003), respectively. (X0/E0) is the initial mixing ratio between xylene and ethylbenzene, and (Xt/Et) is the mixing ratio between xylene and ethylbenzene at the observation time. In this study, we chose the mean concentrations of xylene and ethylbenzene at 05:00–06:00 LT as their initial concentrations before sunrise according to the ambient JNO2 (Fig. S3) to calculate the photochemical loss of OH exposure. In previous work (Shao et al., 2011; Zhan et al., 2021), the selection of ethylbenzene and xylene as tracers was justified for calculating ambient OH exposure under the following conditions: (1) the concentrations of xylene and ethylbenzene were well correlated (Fig. S4), which indicated that they were simultaneously emitted; (2) they had different degradation rates in the atmosphere; and (3) the calculated PICs were in good agreement with those calculated using other tracers, such as i-butene/propene (Fig. S5) (Zhan et al., 2021). To test the relative constant emission ratio from different sources, we chose benzene vs. acetylene and n-hexane vs. toluene as references, and the result is shown in Fig. S6. These ambient ratios could directly reflect their relative emission rates from sources (Goldan et al., 2000; Jobson et al., 2004). The linear correlation coefficients (R2) were generally higher than 0.7, which were equal to that reported by Shao et al. (2011). To further test the assumption that the emissions of xylene and ethylbenzene were constant throughout the day, their potential sources were calculated using a source–receptor model (the potential source contribution function, PSCF). As shown in Fig. S7, xylene and ethylbenzene showed similar distributions. In addition, the ratio of ethylbenzene/xylene at 05:00 and 06:00 LT was similar to that during the daytime. These results indicated that the emissions of xylene and ethylbenzene were constant throughout the day. The ratio of xylene to ethylbenzene and the OH exposure concentration are shown in Fig. S8. The results showed that the ratio of xylene to ethylbenzene increased gradually (07:00–12:00 LT), which is consistent with the trend of xylene and ethylbenzene. The OH exposure was from 0.82 to 8.1×106 molec. cm−3 h, with a mean daytime value of 4.3±1.9×106 molec. cm−3 h. Accordingly, the mean photochemical ages were 1.7±0.9 h using the mean daytime (8:00–17:00 LT) OH concentrations (4.3±3.1×106 molec. cm−3) calculated based on J(O1D) using the method reported in our previous work (Liu et al., 2020b, 2020c). This meant that VOCs would undergo obvious degradation even during a short range of transport in the atmosphere.

It should be noted that the kOH of isoprene is 9.98×10-11 cm3 molec.−1 s−1 at 298.15 K (Atkinson and Arey, 2003), almost 2 orders of magnitude greater than other VOCs. The ratio method assumes constant emissions for VOCs. However, the emission of isoprene greatly depends on temperature and solar irradiation intensity (Zhang et al., 2021b). In addition to accounting for photochemical loss, additional correction of daytime isoprene concentrations was performed using the average diurnal flux of isoprene emissions (Fig. S9; Zhang et al., 2021b). The emission of isoprene showed a clear unimodal curve, and the volume concentration of isoprene was calculated based on the daily emission curve using Eq. (S1) in the Supplement.

2.3 Observation-based model simulation

A box model based on the Master Chemical Mechanism (MCM3.3.1) and the Regional Atmospheric Chemical Mechanism (RACM2) was used in this study. The MCM3.3.1 was used to understand the instantaneous ozone formation process, and the RACM2 was used to depict the ozone isopleth due to its high computational efficiency (Sect. 2.4). Table S1 in the Supplement shows the model inputs. The model calculations were constrained with the measured meteorological parameters (RH, T, P, and JNO2) and the concentrations of trace gases, including inorganic species (NO, NO2, CO, SO2, and HONO) and 61 organic species (NMHCs (46), OVOCs (8), and halohydrocarbons (7)). The model was validated using the observed and simulated O3 concentrations, which showed good consistency, as shown in Fig. S10. The slope and correlation coefficients were 0.9 and 0.8, respectively (Fig. S11), indicating the validity of the model simulation. It is worth mentioning that the results of model simulation can sometimes be overestimated or underestimated to some extent, which has also been reported by previous studies (Zong et al., 2018; Zhang et al., 2020), but this did not affect our simulations of the ozone formation process and mechanisms because we constrained the ozone concentration during our simulations.

The ozone formation rate P(O3) can be quantified by the oxidation rate of NO to NO2 by peroxyl radicals (Tan et al., 2019), as expressed in Eq. (3). In this study, the modeled peroxyl radical concentrations were used to calculate the ozone production rate as follows

(3) P O 3 = k HO 2 + NO HO 2 NO + k RO 2 + NO RO 2 NO ,

where P(O3) is the ozone formation rate, [HO2] and [RO2] are the number concentrations of HO2 and RO2 radicals, kHO2+NO is the second reaction rate between HO2 and NO, and kRO2+NO is the second reaction rate for the reaction of RO2 and NO, which only produces RO and NO2. Once ozone forms, it will be consumed by OH, HO2, and alkenes. Additionally, some NO2 can react with OH, resulting in the formation of nitrate before photolysis. The chemical loss of both O3 and NO2 is considered in the calculation of the net ozone production rate (Tan et al., 2019),

(4) L O 3 = k O 3 + OH OH + k O 3 + HO 2 HO 2 + k O 3 + alkenes alkenes O 3 + k NO 2 + OH NO 2 OH ,

where L(O3) is the ozone chemical loss rate; [OH] is the number concentration of OH radical; kO3+OH, kO3+HO2, and kO3+alkenes are the second-order reaction rate constants between O3 and OH, HO2 and alkenes, respectively; and kNO2+OH is the second-order reaction rate constant between NO2 and OH. Finally, F(O3) is the net ozone formation rate calculated by the difference between P(O3) and L(O3), as expressed in Eq. (5),

(5) F O 3 = P O 3 - L O 3 .

2.4 Empirical kinetic modelling approach

The empirical kinetic modelling approach (EKMA) used in this work is a set of imaginary tests to reveal the dependence of photochemical oxidation products on the change in precursors. We set up a 30×30 matrix by reducing or increasing the measured VOCs and NOx concentrations in the model input. The resulting radical concentrations and ozone production rates were calculated correspondingly.

At this stage, the observed VOCs were grouped into different lumped species according to their RACM2 classification; more details can be found in a previous publication (Tan et al., 2017). The chemical model simulated photochemical reactions with input species for a time interval of 60 min, which was enough for NOx, OH, HO2, and RO2 to reach a steady state because the typical relaxation time of the chemical system is 5–10 min in summer (Tan et al., 2018). However, all the species and parameters were input at a 5 min interval by data interpolation to reduce simulation inconsistencies and large distortions of meteorological parameters at longer time intervals (Tan et al., 2018). The ozone production rate was calculated as described in Sect. 2.3. It is worth mentioning that the average survey data were selected as the baseline scenario in simulating the EKMA curve in this study.

3 Results and discussion

3.1 Overview of diurnal variation in O3, NOx, and TVOC

Figure 1 shows the average diurnal variation of concentrations in O3, NOx, and TVOC (including alkanes, alkenes, OVOCs, and halohydrocarbons) driven by emissions, photochemical reactions, and the evolution of the mixing layer height (MLH). The ozone concentration during the observation period was 44.8±27.2 ppb with a maximum of 119.1 ppb, as reported in our previous study (Zhan et al., 2021), which was generally comparable with the O3 concentrations during 2014–2018 (Ma et al., 2020). The O3 followed a unimodal curve with a minimum value (18.8±15.4 ppb) at 07:00 LT and then it increased to a maximum value (69.6 ppb) at 15:00 LT as photochemical ozone formed. In contrast, NOx reached its maximum concentration (39.7±14.2 ppb) at 07:00 LT and then decreased. After 07:00 LT, the mixing ratio of NO continuously dropped, while the concentration of NO2 decreased at first and then started to increase at 14:00 LT. The diurnal variations in the observed TVOCs were generally consistent with those of NO2. The observed TVOCs concentrations ranged from 2.2 to 23.2 ppb, with a mean value of 18.6±2.6 ppb. Compared to the concentrations (45.4±15.2 ppb) in the same period in August 2015 (Li et al., 2016a), the concentration of VOCs in Beijing was effectively reduced. However, the photochemical initial concentrations (PICs) of TVOCs, which varied from 2.2 to 27.8 ppb with a mean value of 24.5±2.1 ppb, showed a different diurnal curve compared with the observed concentrations. It slightly increased from 07:00 to 14:00 LT, which was similar to the diurnal variation of VOCs in previous work (Zhan et al., 2021). The average PIC-VOC was 6.9±0.5 ppb higher than the observed concentration of TVOCs, indicating an underestimated contribution of the local photochemistry of VOCs to O3 and organic aerosol formation.

Figure 1Overview of average diurnal variations in O3, NOx, and TVOC. The data represent measured results, except for those of the TVOC_PIC, which are calculated based on OH radical exposure. The data range is 1–28 August 2019.


3.2 Influence of photochemical loss of VOCs on the O3 formation sensitivity regime

The sensitivity of O3 formation is analyzed using the isopleth diagram generated from the EKMA model, which is widely used to qualitatively study O3-NOx-VOCs sensitivity. As described in Sect. 2.4, the concentrations of NO2 and VOCs were artificially scaled to ±75 % of the observed values to calculate the response of O3 concentration to an imaginary change in the concentrations of NO2 and VOC, with other constrained conditions remaining unchanged. Figure 2 shows the typical EKMA curves during our observations. The black stars and pentagons denote the observed concentrations of NOx and VOCs in the morning (09:00–10:00 LT) and at noon (14:00–15:00 LT), respectively, while the blue symbols are the corresponding values of PICs. Based on the measured data, O3 formation was in a VOC-limited regime in the morning and a NOx-limited regime in the afternoon. The black arrow indicates a linearly decreasing trend of NOx and VOCs from 09:00 to 15:00 LT in the chemical coordinate system, and ozone production shifted from VOC-limited to NOx-limited conditions from morning to afternoon, which was consistent with the mean diurnal profiles (Fig. 1). This was similar to the data reported in Wangdu (Tan et al., 2018). As expected, ozone production shifted from a VOC-limited regime (the observed VOCs) to a transition regime based on the PIC-VOCs in the morning. Ozone production clearly moved further to a NOx-limited regime in the afternoon after the photochemically consumed NOx and VOCs had been accounted for (Fig. 2). Because the average photochemical aging time was only 1.7±0.9 h, these results indicated that the O3 formation mechanism might typically be misdiagnosed, which misleads mitigation measures for O3 prevention if the consumed VOCs under real atmospheric conditions are not considered.

Figure 2Isopleth diagram of the ozone concentration as a function of the concentration of NOx and VOCs derived from an empirical kinetic modelling approach. The pentagons and stars indicate the status in the morning (09:00–10:00 LT) and at noon (14:00–15:00 LT), respectively. The black and blue colors represent the observed and corrected statuses, respectively.


3.3 Contribution of VOC species to O3 production

The time series of simulated OH, HO2, and RO2 concentrations were used to calculate the P(O3) and L(O3). The diurnally averaged P(O3) and L(O3) are shown in Fig. 3. Ozone formation can be divided into processes related to RO2+NO and HO2+NO (Sect. 2.3). According to their VOC precursors, peroxyl radical groups were divided into alkane-derived (ALKAP), alkene-derived (ALKEP), aromatic-derived (AROMP), isoprene-derived (ISOP), oxygenated-VOC-derived (OVOCP), and halohydrocarbon-derived (HALOP) RO2 and HO2P. The ozone destruction processes included the reaction between O3 and HOx (O3D1), the reaction between O1D and H2O (O3D2), the reaction between O3 and alkenes (O3D3), and the reaction between NO2 and OH (O3D4).

Based on the observed VOCs (or PIC-VOCs), a fast O3 production rate was observed at 14:00 LT (or 13:00 LT), with a diurnal maximum value of 16.1 (or 25.6) ppb h−1 (Fig. 3a and 3b), while the peak destruction rate was 6.4 (or 8.6) ppb h−1 at 15:00 LT (or 13:00 LT) (Fig. 3c and 3d). The average daytime P(O3) from 07:00 to 19:00 LT based on the initial concentrations of VOCs was 4.0±3.1 ppb h−1 higher than that based on the measured VOCs concentrations (Fig. 3b). At the same time, the F(O3) from 07:00 to 19:00 LT based on the initial concentrations of VOCs was also 3.0±2.1 ppb h−1 higher than the measured counterpart (Fig. S12). Thus, the net O3 production could be accumulatively underestimated by ∼36 ppb d−1 from 07:00 to 19:00 LT if the consumption of VOCs was not considered. This meant that the contribution of the local formation of O3 could be underestimated using the directly measured VOCs concentrations. It should be pointed out that it is better to compare O3 production with the true metric for O3 production. However, it is impossible to directly measure the true metric for O3 production in the atmosphere at the present time to know how well the method presented here corrects for that underestimation. In addition, the ozone concentrations must be constrained when simulating the ozone formation process (Lu et al., 2013; Tan et al., 2017). Thus, it is impossible to directly compare the ozone production based on PIC-VOCs with that using measured VOCs concentrations. Therefore, we alternatively compared the integrated net ozone production rates rather than ozone production or concentrations between the two scenarios. An upwind O3 and VOCs measurement combined with a trajectory analysis might provide an approach for checking the accuracy of our results. Alternatively, conducting a transient O3 production rate analysis after subtracting the transport of O3 with a regional model and/or satellite observation might be another option. Unfortunately, neither the upwind measurement nor the regional model simulation was available at the time of our study. To further check the accuracy of our results, we chose 4 August as a test case to explore the influence of the transport of ozone on a downwind site based on the trajectory analysis. As shown in Fig. S13, the mean ozone concentration of the downwind site (national monitoring station, NMS) was 27.6±21.9 ppb d−1 higher than that of the observation site (OS), which was slightly less than the difference (∼36 ppb d−1) between PIC-VOCs and observed VOCs and indirectly rationalized our results.

The HO2 path contributed 64.8 % to the total ozone formation on average, which was slightly higher than the reported value (57.0 %) in Wangdu (Tan et al., 2018), whereas the RO2 path, in which aromatics (9.4 %), alkenes (8.4 %), isoprene (7.8 %), alkanes (4.7 %), OVOCs (4.3 %), and halohydrocarbons (0.6 %) were the main contributors, contributed to the remaining part. For the PIC-VOCs, the dominant path of O3 production (51.7 %) was still the HO2 path, followed by the RO2 path related to alkenes (14.7 %), aromatics (12.8 %), and isoprene (11.7 %). The relative contribution of the RO2 path to P(O3) increased by 13.4 % compared with the measured VOCs, particularly alkene-derived RO2, which increased by 10.2 %. As shown in Fig. 3c and 3d, the destruction of total oxidants was dominated by the reaction between O3 and alkenes (O3D3) in the morning. It gradually shifted to the reaction between NO2 and OH (O3D4) from 11:00 to 16:00 LT and the photolysis of O3 followed by a reaction with water (O3D2) from 12:00 to 15:00 LT because O3 concentration increased while NO2 decreased (Fig. 3c). Figure S14 shows the percentages of the different paths of P(O3) and L(O3). The relative contributions of the reactions between O3 and alkenes (O3D3) and between NO2 and OH (O3D4) to the O3 sinks decreased when calculated based on PIC-VOCs compared with those of the measured VOCs, while they obviously increased for the other two paths, i.e., O3D1 and O3D2. The O3 destruction of the HOx and O3 reaction (O3D1) gradually increased with the continuous photochemical reaction. In addition, the maximum O3 formation rates of the RO2 derived from OVOCs and halohydrocarbons were 0.75 and 0.18 ppb h−1, respectively. These values could be underestimated due to the incomplete gas reaction mechanism of OVOCs and halohydrocarbons in MCM3.3.1. In general, the measured VOCs as model inputs could fail to truly reflect the oxidation capacity and underestimate the local formation of O3 and organic aerosols (Zhan et al., 2021).

Figure 3Mean diurnal profile of the instantaneous ozone production and destruction rate calculated from the MCM-OBM model (instantaneous ozone rate derived from observed VOCs in a and c and from PIC-VOCs in panels b and d). The upper panels present the speciation of the ozone formation rate. The lower panels present the speciation of the ozone destruction rate. The data range is 1–28 August 2019.


The budget of OH–HO2–RO2 radicals was further analyzed to understand the photochemical O3 formation process. The comparison of the radical budget derived from the observed and PIC-VOCs is shown in Fig. 4. The radical cycles are divided into radical sources (green boxes), radical sinks (black boxes), radical propagations (red circles), and equilibria between radical and reservoir species (yellow boxes). The numbers or percentages are the average formation rates (ppb h−1) or relative contributions of the corresponding reaction path based on the observed VOCs (outside the brackets) and the PIC-VOCs (inside the brackets) to a certain radical. The relative contributions of different radical paths based on the observed VOCs (outside the brackets) were comparable with those reported in Beijing, Shanghai, and Guangzhou (Tan et al., 2019), while variations were observed for some reaction paths based on the PIC-VOCs. For example, the reaction between ozone and alkenes based on initial VOC concentrations (percentages inside the brackets) contributed more to OH (from 7 % to 21 %) and HO2 radical production (from 6 % to 12 %), while photolysis of HONO and HCHO contributed less to the production of OH (from 76 % to 60 %) and HO2 radicals (from 44 % to 40 %), respectively. Other radical sources were consistent between the two scenarios. Interestingly, the average formation rates of OH, HO2, and RO2 radicals derived from the PIC-VOCs were obviously higher than those from the observed VOCs. In particular, the oxidation of NO by RO2 and HO2 increased by 1.6 and 1.3 ppb h−1, respectively. The enhanced oxidation rate of NO was equal to the increase in the average F(O3) in the analysis process above. This meant that the radical propagation of OH–RO2–HO2 sped up in the case of PIC-VOCs, subsequently accelerating the chemical loop of NO-NO2-O3. For the radical sinks and equilibria related to HNO4, RONO2, and PAN, the values were basically comparable between the two scenarios. In addition, the O3 formation from the RO2 path increased by 4.1 % (from 39.5 % to 43.6 %) in the simulation using the PIC-VOCs compared with the observed VOCs. The above budget analysis explained the observed increases in F(O3) (∼3 ppb h−1), which were mainly driven by the reaction of missed reactive VOCs, such as alkenes, with O3.

Figure 4Comparison of the OH–HO2–RO2 radical budget derived from the observed and PIC-VOCs under daytime conditions (07:00 to 19:00 LT). The green, black, red, and yellow boxes denote the sources of radicals, radical sinks, radical propagation, and racial equilibrium, respectively. The numbers or percentages outside and inside the brackets are the average formation rates (ppb h−1) or relative contributions to a specific radical of the corresponding reaction path based on observed VOCs and PIC-VOCs, respectively.


3.4 In situ O3 formation process

In addition to chemical processes, which can be simulated using the OBM-MCM model, transport processes, including horizontal, vertical transportation, and dry deposition processes (Tan et al., 2019), also have an important influence on the O3 concentration. Thus, the change in instantaneous ozone concentration can reflect the combined effect between photochemical and physical transport processes (Tan et al., 2019). This change can be expressed as follows:

(6) dO x d t = F O 3 + R O 3 ,

where dOx/dt is the O3 concentration change rate based on the measured data (ppb h−1), F(O3) is the net O3 formation rate (ppb h−1), and R(O3) indicates transportation (ppb h−1). A positive value of R(O3) indicates an inflow of O3 with air mass and vice versa. O3 was replaced with Ox (O3+ NO2) to correct the titration of O3 by NO (Pan et al., 2015).

Figure 5The variation in Ox concentration and formation rate during an O3 pollution episode (1 August). Panels (a) and (b) present the local ozone formation processes of the measured and PIC-VOCs, respectively.


The O3 budget analysis was performed during an O3 pollution episode (1 August). Figure 5 shows the simulated local ozone formation process based on the measured and PIC-VOCs. The hourly variation in O3 concentrations from 19:00 to 6:00 LT the next day was dominated by regional transportation without O3 formation, while local photochemical O3 formation could explain all or part of the O3 concentration change during the time window from 07:00 to 19:00 LT. The d(O3)/dt shows an increase from 07:00 to 15:00 LT. However, d(O3)/dt sharply changed to negative values at 16:00 LT, which was consistent with diurnal O3 (the O3 peaks at 15:00 LT) in Fig. 1.

The average daytime F(O3) based on the observed and photochemical initial concentrations was 6.4±4.0 and 8.9±6.7 ppb h−1, respectively. Photochemical O3 formation under both conditions started at 07:00 LT and reached maximum values of 12.6 and 19.6 ppb h−1 at 15:00 LT, respectively. The maximum daily value of P(O3) was higher than those in the urban areas of Japan, America, and England (Whalley et al., 2018; Ren et al., 2006; Griffith et al., 2016; Kanaya et al., 2009) and lower than those in the suburbs of Guangzhou (Lu et al., 2012) and the urban areas and suburbs of Beijing (Lu et al., 2013). Before 12:00 LT, the O3 formation rate based on the PIC-VOCs was slightly higher than that based on the measured VOCs, while both rates were within a range of 2.0–6.5 ppb h−1. From 12:00 to 17:00 LT, the O3 formation rate based on the PIC-VOCs and the observed concentration of VOCs greatly increased due to active photochemistry.

As shown in Fig. 5, the increased O3 concentration was larger than the local O3 photochemical production from 07:00 to 12:00 LT (R(O3) was positive). This was mainly because the nighttime residual layer (RL) is isolated from mixing with the nighttime surface layer under stable conditions (Tan et al., 2021). The RL usually contains an air mass with a higher ozone mixing ratio than in the surface layer. In the morning, surface heating causes mixing upward in the surface layer until the temperature inversion is eroded away and rapid mixing of pollutants throughout the surface and boundary layer occurs. However, R(O3) was negative in the afternoon, which indicated that the local O3 formation at the measurement site contributed to not only the changes in the in situ O3 concentration but also the O3 source of the downwind regions. This was more clearly shown in Fig. 4b under the PIC-VOC conditions. These results illustrated that local O3 photochemistry played a crucial role in both the local and regional O3 concentrations, which can be underestimated if consumed VOCs with high reactivities are ignored.

4 Conclusions

In this study, we presented the local O3 formation process in August 2019 in Beijing based on the concentrations of observed VOCs and PIC-VOCs. The mean diurnal profile of O3 was unimodal with a peak at 15:00 LT, while NOx and observed TVOCs showed an opposite diurnal curve, and the PICs of TVOCs showed a different diurnal curve compared with that of the observed VOCs, with a slight increase from 07:00 to 14:00 LT. The EKMA curve indicated that instantaneous O3 production was dependent on the real-time concentrations of NOx and VOCs, i.e., the VOC-limited regime in the morning (09:00–10:00 LT) and the NOx-limited regime at noon (14:00–15:00 LT). The sensitivity regime of O3 formation could be misdiagnosed if the consumed VOCs are not considered; for example, the VOC-limited regime (observed) shift to a transition regime (PIC-VOCs) in the morning is ignored. The mean F(O3) based on PIC-VOCs was 3.0 ppb h−1 higher than that based on the measured VOCs, indicating that the underestimation of local photochemistry in the local O3 concentration could reach ∼36 ppb d−1 if the consumed VOCs are not accounted for. The mean ozone concentration of downwind site was 27.6 ppb d−1 higher than the observation site, slightly lower than the difference (∼36 ppb d−1) between PIC-VOCs and observed VOCs, which indirectly supported the accuracy of the above results. The radical budget analysis explained the observed increases in F(O3) (3 ppb h−1), which were mainly driven by the reaction of missed reactive VOCs, such as alkenes, with O3. In addition, the OH–HO2 radical cycle was obviously accelerated by highly reactive alkenes after the photochemical loss of VOCs was accounted for. Finally, the results of the in situ O3 formation process indicated that local O3 photochemical formation played a key role in both local and regional O3 concentrations. In conclusion, our results suggested that PIC-VOCs were more suitable than the observed VOC concentrations for diagnosing O3 formation sensitivity.

Data availability

Data are available upon request to Yongchun Liu (


The supplement related to this article is available online at:

Author contributions

WM contributed to the methodology, data curation, and writing of the original draft. ZF contributed to the methodology, investigation, data curation, and writing of the original draft. JZ contributed to methodology, investigation, and data curation. YL contributed to the conceptualization, investigation, data curation, reviewing and editing the text, supervision, and funding acquisition. PFL contributed to the methodology, investigation, data curation, and reviewing and editing the text. CL contributed to the methodology, investigation, data curation, and reviewing and editing the text. QM contributed to methodology, investigation, and data curation. KY contributed to the methodology, investigation, and data curation. YW contributed to the methodology, investigation, resources, and data curation. HH contributed to the acquisition of resources and reviewing and editing the text. MK contributed to the methodology and reviewing the text. YM contributed to the conceptualization, methodology, data curation, and reviewing and editing the text. JL contributed to the conceptualization, methodology, data curation, reviewing and editing the text, and supervision.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


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


This research was financially supported by the Strategic Priority Research Program of the Chinese Academy of Sciences and the Beijing University of Chemical Technology.

Financial support

This research has been supported by the National Natural Science Foundation of China (grant nos. 41877306, 92044301, and 21976190), the Ministry of Science and Technology of the People's Republic of China (grant no. 2019YFC0214701) and the specific research fund of the Innovation Platform for Academicians of Hainan Province (grant no. YSPTZX202205).

Review statement

This paper was edited by Kelley Barsanti and reviewed by two anonymous referees.


Atkinson, R. and Arey, J.: Atmospheric degradation of volatile organic compounds, Chem. Rev., 103, 4605–4638,, 2003. 

Carter, W. P. L.: Development of Ozone Reactivity Scales for Volatile Organic Compounds, J. Air Waste Manage. Assoc., 44, 881–899,, 1994. 

Chen, T., Liu, J., Liu, Y., Ma, Q., G. Y., Zhong, C., Jiang, H., Chu, B., Zhang, P., Ma, J., Liu, P., Wang, Y., Mu, Y., and He, H.: Chemical characterization of submicron aerosol in summertime Beijing: A case study in southern suburbs in 2018, Chemosphere, 247, 125918,, 2020. 

Cohen, A. J., Brauer, M., Burnett, R., Anderson, H., Frostad., J., Estep, K., Balakrishnan, K., Brunekreef, B., Dandona, L., Dandona, R., Feigin, V., Freedman, G., Hubbell, B., Jobling, A., Kan, H., Knibbs,L., Liu, Y., Martin, R., Morawska, L., Arden Pope III, C., Shin, H., Straif, K., Shaddick, G., Thomas, M., van Dingenen, R., van Donkelaar, A., Vos, T., Murray, C., and Forouzanfar, M.: Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the Global Burden of Diseases Study 2015, Lancet, 389, 1907–1918,, 2017. 

Dang, R. and Liao, H.: Radiative Forcing and Health Impact of Aerosols and Ozone in China as the Consequence of Clean Air Actions over 2012–2017, J. Geophys. Res. Lett., 46, 12511–12519,, 2019. 

Feng, Z., De Marco, A., Anav, A., Gualtieri, M., Sicard, P., Tian, H., Fornasier, F., Tao, F., Guo, A., and Paoletti, E.: Economic losses due to ozone impacts on human health, forest productivity and crop yield across China, Environ. Int., 131, 104966,, 2019. 

Gao, J., Zhang, J., Li, H., Li, L., Xu, L., Zhang, Y., Wang,Z., Wang, X., Zhang, W., Chen, Y., Cheng, X,. Zhang, H., Peng, L., Chai, F., and Wei, Y.: Comparative study of volatile organic compounds in ambient air using observed mixing ratios and initial mixing ratios taking chemical loss into account – A case study in a typical urban area in Beijing, Sci. Total Environ., 628–629, 791–804,, 2018. 

Gao, Y., Li, M., Wan, X., Zhao, X., Wu, Y., Liu, X., and Li, X.: Important contributions of alkenes and aromatics to VOCs emissions, chemistry and secondary pollutants formation at an industrial site of central eastern China, Atmos. Environ., 244, 117927,, 2021. 

Goldan, P. D., Parrish, D. D., Kuster, W. C., Trainer, M., Mckeen, S. A., Holloway, J., Jobson, B. T., Sueper, D. T., and Fehsenfeld, F. C.: Airborne measurements of isoprene, CO, and anthropogenic hydrocarbons and their implications, J. Geophys. Res.-Atmos., 105, 9091–9105,, 2000. 

Griffith, S. M., Hansen, R. F., Dusanter, S., Michoud, V., Gilman, J. B., Kuster, W. C., Veres, P. R., Graus, M., de Gouw, J. A., Roberts, J., Young, C., Washenfelder, R., Brown, S. S., Thalman, R., Waxman, E., Volkamer, R., Tsai, C., Stutz, J., Flynn, J. H., Grossberg, N., Lefer, B., Alvarez, S. L., Rappenglueck, B., Mielke, L. H., Osthoff, H. D., and Stevens, P. S.: Measurements of hydroxyl and hydroperoxy radicals during CalNex-LA: Model comparisons and radical budgets, J. Geophys. Res.-Atmos., 121, 4211–4232,, 2016. 

Han, D., Wang, Z., Cheng, J., Wang, Q., Chen, X., and Wang, H.: Volatile organic compounds (VOCs) during non-haze and haze days in Shanghai: characterization and secondary organic aerosol (SOA) formation, Environ. Sci. Pollut. R., 24, 18619–18629,, 2017. 

He, Z., Wang, X., Ling, Z., Zhao, J., Guo, H., Shao, M., and Wang, Z.: Contributions of different anthropogenic volatile organic compound sources to ozone formation at a receptor site in the Pearl River Delta region and its policy implications, Atmos. Chem. Phys., 19, 8801–8816,, 2019. 

Jobson, B. T., Berkowitz, C. M., Kuster, W. C., Goldan, P. D., Williams, E. J., Fesenfeld, F. C., Apel, E. C., Karl, T., Lonneman, W. A., and Riemer, D.: Hydrocarbon source signatures in Houston, Texas: Influence of the petrochemical industry, J. Geophys. Res.-Atmos., 109, D23405,, 2004. 

Kanaya, Y., Pochanart, P., Liu, Y., Li, J., Tanimoto, H., Kato, S., Suthawaree, J., Inomata, S., Taketani, F., Okuzawa, K., Kawamura, K., Akimoto, H., and Wang, Z. F.: Rates and regimes of photochemical ozone production over Central East China in June 2006: a box model analysis using comprehensive measurements of ozone precursors, Atmos. Chem. Phys., 9, 7711–7723,, 2009. 

Li, J., Wu, R., Li, Y., Hao, Y., and Xie, S. Z., Liming;: Effects of rigorous emission controls on reducing ambient volatile organic compounds in Beijing, China, Sci. Total Environ., 557, 531–541,, 2016a. 

Li, J., Yang, W., Wang, Z., Chen, H., Hu, B., Li, J., Sun, Y., Fu, P., and Zhang, Y.: Modeling study of surface ozone source-receptor relationships in East Asia, Atmos. Res., 167, 77–88,, 2016b. 

Li, L., Xie, S., Zeng, L., Wu, R., and Li, J.: Characteristics of volatile organic compounds and their role in ground-level ozone formation in the Beijing-Tianjin-Hebei region, China, Atmos. Environ., 113, 247–254,, 2015. 

Li, L., An, J., Shi, Y., Zhou, M., Yan, R., Huang, C., Wang, H., Lou, S., Wang, Q., Lu, Q., and Wu, J.: Source apportionment of surface ozone in the Yangtze River Delta, China in the summer of 2013, Atmos. Environ., 144, 194–207,, 2016c. 

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. 

Ling, Z. H. and Guo, H.: Contribution of VOC sources to photochemical ozone formation and its control policy implication in Hong Kong, Environ. Sci. Policy, 38, 180–191,, 2014. 

Liu, J., Li, X., Tan, Z., Wang, W., and Zhang, Y.: Assessing the Ratios of Formaldehyde and Glyoxal to NO2 as Indicators of O3–NOx–VOC Sensitivity, Environ. Sci. Technol., 55, 10935–10945,, 2021. 

Liu, Y., Yan, C., Feng, Z., Zheng, F., and Kulmala, M.: Continuous and comprehensive atmospheric observations in Beijing: a station to understand the complex urban atmospheric environment, Big Earth Data, 4, 295–321,, 2020a. 

Liu, Y., Ni, S., Jiang, T., Xing, S., Zhang, Y., Bao, X., Feng, Z., Fan, X., Zhang, L., and Feng, H.: Influence of Chinese New Year overlapping COVID-19 lockdown on HONO sources in Shijiazhuang, Sci. Total Environ., 745, 141025,, 2020b. 

Liu, Y., Zhang, Y., Lian, C., Yan, C., Feng, Z., Zheng, F., Fan, X., Chen, Y., Wang, W., Chu, B., Wang, Y., Cai, J., Du, W., Daellenbach, K. R., Kangasluoma, J., Bianchi, F., Kujansuu, J., Petäjä, T., Wang, X., Hu, B., Wang, Y., Ge, M., He, H., and Kulmala, M.: The promotion effect of nitrous acid on aerosol formation in wintertime in Beijing: the possible contribution of traffic-related emissions, Atmos. Chem. Phys., 20, 13023–13040,, 2020c. 

Lu, K. D., Rohrer, F., Holland, F., Fuchs, H., Bohn, B., Brauers, T., Chang, C. C., Häseler, R., Hu, M., Kita, K., Kondo, Y., Li, X., Lou, S. R., Nehr, S., Shao, M., Zeng, L. M., Wahner, A., Zhang, Y. H., and Hofzumahaus, A.: Observation and modelling of OH and HO2 concentrations in the Pearl River Delta 2006: a missing OH source in a VOC rich atmosphere, Atmos. Chem. Phys., 12, 1541–1569,, 2012. 

Lu, K. D., Hofzumahaus, A., Holland, F., Bohn, B., Brauers, T., Fuchs, H., Hu, M., Häseler, R., Kita, K., Kondo, Y., Li, X., Lou, S. R., Oebel, A., Shao, M., Zeng, L. M., Wahner, A., Zhu, T., Zhang, Y. H., and Rohrer, F.: Missing OH source in a suburban environment near Beijing: observed and modelled OH and HO2 concentrations in summer 2006, Atmos. Chem. Phys., 13, 1057–1080,, 2013. 

Lu, X., Chen, N., Wang, Y., Cao, W., Zhu, B., Yao, T., Fung, J. C. H., and Lau, A. K. H.: Radical budget and ozone chemistry during autumn in the atmosphere of an urban site in central China, J. Geophys. Res.-Atmos., 122, 3672–3685,, 2017. 

Ma, J., Kamal, J., Li, V., and Lam, J.: Effects of China's current Air Pollution Prevention and Control Action Plan on air pollution patterns, health risks and mortalities in Beijing 2014–2018, Chemosphere, 260, 127572,, 2020. 

Ma, Z., Liu, C., Zhang, C., Liu, P., Ye, C., Xue, C., Zhao, D., Sun, J., Du, Y., Chai, F., and Mu, Y.: The levels, sources and reactivity of volatile organic compounds in a typical urban area of Northeast China, J. Environ. Sci., 79, 121–134,, 2019. 

Mckeen, S. A., Liu, S. C., Hsie, E. Y., Lin, X., Bradshaw, J. D., Smyth, S., Gregory, G. L., and Blake, D. R.: Hydrocarbon ratios during PEM-WEST A: A model perspective, J. Geophys Res.-Atmos., 101, 2087–2109,, 1996. 

Monks, P. S.: Gas-Phase Radical Chemistry in the Troposphere, Chem. Soc. Rev., 34, 376–395,, 2005. 

Pan, X., Kanaya, Y., Tanimoto, H., Inomata, S., Wang, Z., Kudo, S., and Uno, I.: Examining the major contributors of ozone pollution in a rural area of the Yangtze River Delta region during harvest season, Atmos. Chem. Phys., 15, 6101–6111,, 2015. 

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, S252–S263,, 2006. 

Seinfeld, J. H. and Pandis, S. N.: Atmospheric chemistry and physics: from air pollution to climate change, 3rd edn., John Wiley & Sons Inc., Hoboken, New Jersey, ISBN 9781119221166, 2016. 

Shao, M., Wang, B., Lu, S., Yuan, B., and Wang, M.: Effects of Beijing Olympics Control Measures on Reducing Reactive Hydrocarbon Species, Environ. Sci. Technol., 45, 514–519,, 2011. 

Tan, Z., Fuchs, H., Lu, K., Hofzumahaus, A., Bohn, B., Broch, S., Dong, H., Gomm, S., Häseler, R., He, L., Holland, F., Li, X., Liu, Y., Lu, S., Rohrer, F., Shao, M., Wang, B., Wang, M., Wu, Y., Zeng, L., Zhang, Y., Wahner, A., and Zhang, Y.: Radical chemistry at a rural site (Wangdu) in the North China Plain: observation and model calculations of OH, HO2 and RO2 radicals, Atmos. Chem. Phys., 17, 663–690,, 2017. 

Tan, Z., Lu, K., Jiang, M., Su, R., Dong, H., Zeng, L., Xie, S., Tan, Q., and Zhang, Y.: Exploring ozone pollution in Chengdu, southwestern China: A case study from radical chemistry to O3-VOC-NOx sensitivity, Sci. Total Environ., 636, 775–786,, 2018. 

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. F., Ma, X. F., Lu, K. D., Jiang, M. Q., Zou, Q., Wang, H. C., Zeng, L. M., and Zhang, Y. H.: Direct evidence of local photochemical production driven ozone episode in Beijing: A case study, Sci. Total Environ., 800, 148868,, 2021. 

Wang, T., Xue, L., Brimblecombe, P., Lam, Y. F., Li, L., and Zhang, L.: Ozone pollution in China: A review of concentrations, meteorological influences, chemical precursors, and effects, Sci. Total Environ., 575, 1582–1596,, 2017. 

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. 

Wiedinmyer, C., Friedfeld, S., Baugh, W., Greenberg, J., Guenther, A., Fraser, M., and Allen, D.: Measurement and analysis of atmospheric concentrations of isoprene and its reaction products in central Texas, Atmos. Environ., 35, 1001–1013,, 2001. 

Xie, Y., Dai, H., Zhang, Y., Wu, Y., Hanaoka, T., and Masui, T.: Comparison of health and economic impacts of PM2.5 and ozone pollution in China, Environ. Int., 130, 104881,, 2019. 

Xue, L. K., Wang, T., Gao, J., Ding, A. J., Zhou, X. H., Blake, D. R., Wang, X. F., Saunders, S. M., Fan, S. J., Zuo, H. C., Zhang, Q. Z., and Wang, W. X.: Ground-level ozone in four Chinese cities: precursors, regional transport and heterogeneous processes, Atmos. Chem. Phys., 14, 13175–13188,, 2014. 

Yang, Y., Shao, M., Keßel, S., Li, Y., Lu, K., Lu, S., Williams, J., Zhang, Y., Zeng, L., Nölscher, A. C., Wu, Y., Wang, X., and Zheng, J.: How the OH reactivity affects the ozone production efficiency: case studies in Beijing and Heshan, China, Atmos. Chem. Phys., 17, 7127–7142,, 2017. 

Yuan, B., Hu, W. W., Shao, M., Wang, M., Chen, W. T., Lu, S. H., Zeng, L. M., and Hu, M.: VOC emissions, evolutions and contributions to SOA formation at a receptor site in eastern China, Atmos. Chem. Phys., 13, 8815–8832,, 2013. 

Zhan, J., Feng, Z., Liu, P., He, X., He, Z., Chen, T., Wang, Y., He, H., Mu, Y., and Liu, Y.: Ozone and SOA formation potential based on photochemical loss of VOCs during the Beijing summer, Environ. Pollut., 285, 117444,, 2021. 

Zhang, K., Li, L., Huang, L., Wang, Y., Huo, J., Duan, Y., Wang, Y., and Fu, Q.: The impact of volatile organic compounds on ozone formation in the suburban area of Shanghai, Atmos. Environ., 232, 117511,, 2020. 

Zhang, K., Huang, L., Li, Q., Huo, J., Duan, Y., Wang, Y., Yaluk, E., Wang, Y., Fu, Q., and Li, L.: Explicit modeling of isoprene chemical processing in polluted air masses in suburban areas of the Yangtze River Delta region: radical cycling and formation of ozone and formaldehyde, Atmos. Chem. Phys., 21, 5905–5917,, 2021a. 

Zhang, M., Gao, W., Yan, J., Wu, Y., Marandino, C. A., Park, K., Chen, L., Lin, Q., Tan, G., and Pan, M.: An integrated sampler for shipboard underway measurement of dimethyl sulfide in surface seawater and air, Atmos. Environ., 209, 86–91,, 2019. 

Zhang, Q., Li, L., Zhao, W., Wang, X., Jiang, L., Liu, B., Li, X., and Lu, H.: Emission characteristics of VOCs from forests and its impact on regional air quality in Beijing, China Environ. Sci., 41, 622–632,, 2021b. 

Zong, R. H., Xue, L. K., Wang, T. B., and Wang, W. X.: Inter-comparison of the Regional Atmospheric Chemistry Mechanism (RACM2) and Master Chemical Mechanism (MCM) on the simulation of acetaldehyde, Atmos. Environ., 186, 144–149,, 2018. 

Short summary
The influence of photochemical loss of volatile organic compounds (VOCS) on O3 formation is investigated using an observation-based model. The sensitivity regime of ozone formation might be misdiagnosed due to the photochemical loss of VOCs in the atmosphere. The contribution of local photochemistry is underestimated regarding O3 pollution when one does not consider the photochemical loss of VOCs.
Final-revised paper