Intercomparison of O3 formation and radical chemistry in the past decade at a suburban site in Hong Kong

Hong Kong, as one of the densely populated metropolises in East Asia, has been suffering from severe photochemical smog in the past decades, though the observed nitrogen oxides (NOx) and total volatile organic compounds (TVOCs) were significantly reduced. This study, based on the observation data in the autumns of 2007, 2013 and 2016, investigated the photochemical ozone (O3) formation and radical chemistry during the three sampling periods in Hong Kong with the aid of a photochemical box model incorporating the Master Chemical Mechanism (PBM–MCM). While the simulated locally produced O3 remained unchanged (p = 0.73) from 2007 to 2013, the observed O3 increased (p < 0.05) at a rate of 1.78 ppbv yr−1 driven by the rise in regionally transported O3 (1.77± 0.04 ppbv yr−1). Both the observed and locally produced O3 decreased (p < 0.05) from the VOC sampling days in 2013 to those in 2016 at a rate of −5.31± 0.07 and −5.52± 0.05 ppbv yr−1, respectively. However, a leveling-off (p = 0.32) was simulated for the regionally transported O3 during 2013–2016. The mitigation of autumn O3 pollution in this region was further confirmed by the continuous monitoring data, which have never been reported. Benefiting from the air pollution control measures taken in Hong Kong, the local O3 production rate decreased remarkably (p < 0.05) from 2007 to 2016, along with the lowering of the recycling rate of the hydroxyl radical (OH). Specifically, VOCs emitted from the source of liquefied petroleum gas (LPG) usage and gasoline evaporation decreased in this decade at a rate of −2.61± 0.03 ppbv yr−1, leading to a reduction of the O3 production rate from 0.51± 0.11 ppbv h−1 in 2007 to 0.10± 0.02 ppbv h−1 in 2016. In addition, solvent usage made decreasing contributions to both VOCs (rate=−2.29±0.03 ppbv yr−1) and local O3 production rate (1.22± 0.17 and 0.14± 0.05 ppbv h−1 in 2007 and 2016, respectively) in the same period. All the rates reported here were for the VOC sampling days in the three sampling campaigns. It is noteworthy that meteorological changes also play important roles in the inter-annual variations in the observed O3 and the simulated O3 production rates. Evaluations with more data in longer periods are therefore recommended. The analyses on the decadal changes of the local and regional photochemistry in Hong Kong in this study may be a reference for combating China’s nationwide O3 pollution in near future.

Abstract. Hong Kong, as one of the densely populated metropolises in East Asia, has been suffering from severe photochemical smog in the past decades, though the observed nitrogen oxides (NO x ) and total volatile organic compounds (TVOCs) were significantly reduced. This study, based on the observation data in the autumns of 2007, 2013 and 2016, investigated the photochemical ozone (O 3 ) formation and radical chemistry during the three sampling periods in Hong Kong with the aid of a photochemical box model incorporating the Master Chemical Mechanism (PBM-MCM). While the simulated locally produced O 3 remained unchanged (p = 0.73) from 2007 to 2013, the observed O 3 increased (p < 0.05) at a rate of 1.78 ppbv yr −1 driven by the rise in regionally transported O 3 (1.77 ± 0.04 ppbv yr −1 ). Both the observed and locally produced O 3 decreased (p < 0.05) from the VOC sampling days in 2013 to those in 2016 at a rate of −5.31 ± 0.07 and −5.52 ± 0.05 ppbv yr −1 , respectively. However, a leveling-off (p = 0.32) was simulated for the regionally transported O 3 during 2013-2016. The mitigation of autumn O 3 pollution in this region was further confirmed by the continuous monitoring data, which have never been reported. Benefiting from the air pollution control measures taken in Hong Kong, the local O 3 production rate decreased remarkably (p < 0.05) from 2007 to 2016, along with the lowering of the recycling rate of the hydroxyl radical (OH). Specifically, VOCs emitted from the source of liquefied petroleum gas (LPG) usage and gasoline evaporation decreased in this decade at a rate of −2.61 ± 0.03 ppbv yr −1 , leading to a reduction of the O 3 production rate from 0.51 ± 0.11 ppbv h −1 in 2007 to 0.10 ± 0.02 ppbv h −1 in 2016. In addition, solvent usage made decreasing contributions to both VOCs (rate = −2.29±0.03 ppbv yr −1 ) and local O 3 production rate (1.22 ± 0.17 and 0.14 ± 0.05 ppbv h −1 in 2007 and 2016, respectively) in the same period. All the rates reported here were for the VOC sampling days in the three sampling campaigns. It is noteworthy that meteorological changes also play important roles in the inter-annual variations in the observed O 3 and the simulated O 3 production rates. Evaluations with more data in longer periods are therefore recommended. The analyses on the decadal changes of the local and regional photochemistry in Hong Kong in this study may be a reference for combating China's nationwide O 3 pollution in near future.

Introduction
Ground-level ozone (O 3 ) is one of the most representative air pollutants in photochemical smog, produced through photochemical reactions between volatile organic compounds (VOCs) and nitrogen oxides (NO x ) in the presence of sunlight (NRC;1992;Jacob et al., 1999;Guo et al., 2017). It is well documented that O 3 is harmful to human health (Bell et al., 2004), crops  and natural ecosystems (Ashmore, 2005). Through the last 30 years, extensive efforts have been made by the local and federal governments to alleviate the tropospheric O 3 pollution around the world (NRC;1992;NARSTO, 2000;. Effectiveness has gradually shown in some countries/regions, such as Switzerland, Germany, Ireland and eastern North America (Lefohn et al., 2010;Cui et al., 2011;Derwent et al., 2013;Parrish et al., 2014;Lin et al., 2017). In contrast, the O 3 levels in many places are still increasing or not decreasing at the expected rates, particularly in East Asia (Ding et al., 2008;Xu et al., 2008;Parrish et al., 2014;Xue et al., 2014a;. Hong Kong, as one of the densely populated metropolises in East Asia, has been suffering from severe photochemical smog in the past decades, though the locally emitted NO x and total VOCs (TVOCs) were significantly reduced (Xue et al., 2014a;Ou et al., 2015;Lyu et al., 2016a;. On the one hand, this indicates the nonlinear relationship between O 3 and its precursors. On the other hand, in addition to local O 3 formation, the observed O 3 in Hong Kong is also influenced by the regional transport due to the proximity of the highly industrialized Pearl River Delta (PRD) region. Earlier studies revealed that the local O 3 production is typically limited by VOCs in urban and some suburban areas in Hong Kong (Zhang et al., 2007;. Namely, cutting VOCs emissions will reduce O 3 production, while the reduction of NO x may cause an O 3 increment (Cheng et al., 2010Guo et al., 2011;. Previous studies also documented that photochemical O 3 formation is dependent upon the ratios between TVOCs and NO x (Sillman, 1999;Guo et al., 2013a;Ling et al., 2013), reactivity of VOC species (Zhang et al., 2007;Liu et al., 2008;Cheng et al., 2010) and the composition of NO x (i.e., relative abundances of NO 2 and NO) (Richter et al., 2005;Xu et al., 2008;. Moreover, located in the subtropical region, Hong Kong has relatively high temperature and strong solar radiation, which are favorable for local O 3 formation. For regional transport, studies (Wang et al., 2001;Ding et al., 2004; indicated that O 3 was generally built up in Hong Kong under the northerly winds, whereas it was often driven down by the sea breeze from the South China Sea (SCS) and by the southwest monsoon in warm seasons. The contribution of regional transport to O 3 in Hong Kong even reached 70 % under the dominance of the tropical cyclone (Huang et al., 2005), a typical synoptic condition conducive to severe O 3 pollution in the Northern Hemisphere (So and Wang, 2003;Huang et al., 2005;Lam et al., 2005). To improve the air quality in Hong Kong, a series of control measures aiming at restriction of VOC emissions have been implemented by Hong Kong government since 2007, which effectively reduced the concentrations of some VOCs, such as propane and i-/n-butanes emitted from taxis and public buses fueled by liquefied petroleum gas (LPG) (Lyu et al., 2016b), the aromatics mainly attributable to solvent usage, and the alkenes in association with diesel exhaust (Lyu et al., 2017a). As a result, Xue et al. (2014a) and Y.  found that the locally produced O 3 decreased. However, the regional and superregional transport of O 3 and its precursors from PRD and eastern China to Hong Kong had offset the decrease in the local O 3 production, resulting in an overall increase in the observed O 3 in Hong Kong from 2005 to 2013. Overall, the previous studies have greatly enhanced our understanding on O 3 pollution in Hong Kong, and details about the studies can be referred to in Table S1 in the Supplement. Despite many previous studies (Xue et al., 2014a(Xue et al., , 2016Ou et al., 2015;Lyu et al., 2016a;, the inter-annual variations in the O 3 formation regimes and radical chemistry have not yet been fully understood in Hong Kong. In addition, the online measurement data used in previous long-term O 3 study might hamper the exact understanding of the local O 3 formation mechanisms, due to the unavailability of many reactive VOCs, such as formaldehyde. Moreover, the trends of the local production and regional transport of O 3 were only updated to 2013 in previous studies (Xue et al., 2014a;. In fact, many measures were taken to reduce air pollutants' emissions in the latest years in Hong Kong and PRD. For examples, nearly 75 % of the old catalytic converters in LPG-fueled vehicles were renewed during September 2013-May 2014. A program to eliminate the pre-Euro IV diesel vehicles or to upgrade their emission standards to Euro IV was initiated in March 2014 and is still ongoing until 2019 at its third phase. In the PRD, the second stage of the clean air controlling program was implemented in (DGEPD, 2013. In 2014, the Guangdong provincial government launched an Action Plan for Air Pollution Prevention and Control (MEE PRC, 2014), putting the emphasis on the emission control of traffic, coalfired power plants and industrial sources. Investigations on the post-2013 variations in the local O 3 production in Hong Kong and the regional impacts provide a good opportunity for us to examine the effectiveness of these local and regional measures.
The objectives of this study were to re-examine the O 3 trend in the pre-2013 emissions and trace the O 3 evolution in the post-2013 emissions in Hong Kong and to explore the underlying mechanisms for the variations in O 3 formation and radical chemistry. With the aid of a photochemical box model, the locally produced and regionally transported O 3 , as well as its variation trends, were determined (see Sect. 2.5). Under the assumption that the local O 3 production in these years was changed due to a series of control measures in Hong Kong, we also aimed to evaluate the actual effectiveness of these control measures. China is suffering from severe O 3 pollution, almost the worst in the world. While O 3 began to decrease in most areas of North America and Europe, China's O 3 pollution was even aggravated in recent years. A series of air pollution control strategies have been implemented in China, though most of them were not specifically designed for O 3 abatement. Investigations on O 3 trends and the potential causes in Hong Kong would provide a good example of assessing the evolution of O 3 pollution and the effects of artificial interventions in China. In addition, the changes in the regional contribution to O 3 in Hong Kong determined in this study would throw light upon the variations in O 3 in China, particularly in south China. It is expected that this study would have some inspiration for the O 3 pollution control in other cities and regions in China.

Sampling site
Hong Kong is located on the southern coast of China with Guangdong Province to the north and the Pearl River estuary (PRE) to the northwest. The sampling site (22.29 • N, 113.94 • E), Tung Chung (TC), is in a newly developed suburban area in western Hong Kong, with a population of ∼ 77 400 in 2016 (CSD, 2011(CSD, , 2018. The urban center of Hong Kong is ∼ 20 km northeast of TC. Hong Kong is dominated by the subtropical oceanic monsoon climate. During warm seasons, the prevailing winds mainly come from the SCS at a relatively low speed (southwest winds). In cold seasons, the east and northeast winds are predominant. Generally, the sampling site receives relatively polluted air masses from mainland China, i.e., the PRD region, Yangtze River Delta region and even north China between October and March, when high O 3 levels are often observed . Therefore, the samplings were mainly conducted in October and November in this study, except for 4 out of 45 sampling days in September. The sampling site was close to a highway running to the Hong Kong International Airport (HKIA), and the HKIA was around 3 km to the north of the site. In addition, the local emissions from residential activities may modulate the air quality at this site. It was expected that O 3 at TC would be significantly influenced by NO emitted from the aforementioned sources. As shown in Fig. S1 in the Supplement, the O 3 titrated by NO ( O 3 −1) was equivalent to 29.94 ± 8.8 %, 26.7 ± 12.7 % and 32.5 ± 16.6 % of the photochemically formed O 3 ( O 3 −2) in the 2007, 2013 and 2016 sampling campaigns, respectively, confirming the importance of NO titration in modulating O 3 at the site. Figure 1 shows the locations of the sampling site (TC) and the 12 air quality monitoring stations in the PRD, which witnessed the evolution of air quality in the PRD over the last decade and are used to demonstrate the variations in regional O 3 in this study. More detailed description of the site can be found in our previous studies (Cheng et al., 2010;Jiang et al., 2010;Ling et al., 2013;Ou et al., 2015).

Continuous measurements of trace gases and collection of VOC and OVOC samples
Trace gases (SO 2 , CO, NO, NO 2 and O 3 ) and meteorological conditions were continuously measured at TC site for three autumn periods in 2007, 2013 and 2016 (see Table S2 Table S2 for the specific sampling dates). The three sampling periods were used as representatives of the autumns in the three years in this study, and the rationality will be discussed in Sect. 3.1. Trace gases were continuously measured at the TC air quality monitoring station operated by the Hong Kong Environmental Protection Department (HKEPD), ∼ 0.8 km from our sampling site. The instruments were the same as those used in the US air quality monitoring program (HKEPD, 2017a). Table S3 summarizes the instruments, analysis techniques, detection limits and the time resolutions for measurements of the trace gases. The high-resolution data were collected and averaged into the hourly averages. All the analyzers except the O 3 analyzer were zeroed daily by analyzing scrubbed ambient air and calibrated every 2 weeks by a span gas mixture with a NIST (National Institute of Standards and Technology) traceable standard, while the O 3 analyzer was calibrated using a transfer standard (Thermo Environmental Instruments (TEI) 49PS) every 2 weeks. Details about the quality assurance and control procedures can be found in Ling et al. (2016a). The meteorological parameters, including temperature, relative humidity, pressure, wind speed, wind direction, precipitation and solar radiation, were also continuously monitored by a mini weather station (Vantage Pro TM & Vantage Pro2 Plus TM weather stations, Davis Instruments) during the sampling periods. Data were integrated into 30 min averages by a built-in program in the weather station. The collection and analysis of VOCs and OVOCs were detailed in our previous studies (Guo et al., 2009;. Briefly, pre-cleaned and evacuated 2 L electropolished stainless-steel canisters were used to collect VOC samples. On O 3 episode days, a 1 h sample was collected each hour during the daytime (07:00-19:00 LT; all times throughout the paper are in local time), generating 13 samples per day, while five to seven 1 h samples were collected every other hour on non-O 3 episode days from 07:00 to 19:00 LT in the 2013 and 2016 sampling campaigns. However, 12 1 h samples were collected on each VOC sampling day between 07:00 and 18:00 in 2007, regardless of O 3 episodes or nonepisodes. The O 3 episode days were predicted prior to sampling based on weather forecast and numerical simulation of O 3 . Overall, the O 3 episodes were usually associated with high temperature, strong solar radiation, low humidity, and weak or northerly winds. A total of 414 canister samples, including 96 samples in 2007, 146 samples in 2013 and 172 samples in 2016, were collected and analyzed during the three sampling periods (Table S2).
In addition to VOC samples, OVOC samples were also collected on the same days as those for the collection of VOCs. Dinitrophenylhydrazine (DNPH)-silica cartridges (Waters Sep-Pak DNPH-Silica, Milford, MA) were used to collect the OVOC samples. An ozone scrubber (Sep-Pak; Waters Corporation, Milford, MA) was connected in front of the DNPH cartridge to prevent interference of ozone. The ozone scrubber was replaced every two OVOC samples. For each OVOC sample, air was drawn to pass the O 3 scrubber and the cartridge for 2 h (2.5 h in 2007 sampling campaign) at a flow rate of 0.5 L min −1 , which was controlled by a rotameter. During the sampling periods in 2013 and 2016, five to seven OVOC samples were collected every 2 h from 06:00 to 20:00 LT on both O 3 episode and non-episode days. In 2007, only two samples were collected on non-O 3 episode days at 10:30-13:00 and 13:00-15:30 and four samples between 08:00 and 18:00 on O 3 episode days. In total, 275 OVOC samples (28 in 2007, 124 in 2013 and 124 in 2016) were collected and analyzed in the three sampling campaigns (Table S2).

Chemical analysis 2.3.1 Analysis of VOCs
The concentrations of 48 speciated non-methane hydrocarbons (NMHCs) in the canisters were determined with an Entech model 7100 preconcentrator (Entech Instruments Inc., California, USA) coupled with a gas chromatography mass selective detector (model 5973N, Agilent Technologies, USA), a flame ionization detector, and an electron capture detector (GC-MSD-FID-ECD). The NMHCs were analyzed in the Guangzhou Institute of Geochemistry (GIG), Chinese Academy of Sciences for the samples collected in both 2007 and 2013, and in The Hong Kong Polytechnic University (HKPolyU) for the samples collected in 2016. It should be noted that the GC-MSD-FID-ECD system in the latter two institutes was the same as that at UCI, and intercomparisons were performed regularly among the three institutes, which showed reasonably good agreements Zeng et al., 2018). Detailed information about the analysis procedures and quality assurance and control can be found in Colman et al. (2001) and Simpson et al. (2010). Table S4 summarizes the limits of detection (LoDs), precisions and accuracies of the VOC analyses in the three institutes.
The OVOC samples were stored in a refrigerator at 4 • C after sampling. For analyses of OVOCs, the cartridges were eluted slowly with 2 mL of acetonitrile into a 2 mL volumetric flask. A high-performance liquid chromatography (HPLC) system (PerkinElmer series 2000, MA, USA) coupled with an ultraviolet (UV) detector operating at 360 nm was used for analysis. The instrument was calibrated using standards of five gradient concentrations covering the concentrations of interest for different OVOCs in ambient air. Good linear relationships (R 2 > 0.999) between the standard concentrations and responses of the instrument were obtained for the 16 analyzed OVOC species. The builtin computerized programs of quality control systems such as auto-linearization and autocalibration were used to guarantee the data quality. Detailed information about the analysis and quality control of OVOC samples was provided in Cheng et al. (2014), Cui et al. (2016) and Ling et al. (2016b). Due to the low detection rate of many OVOCs, this study only focused on formaldehyde, acetaldehyde, acetone and propionaldehyde, which had relatively high concentrations.

Positive matrix factorization (PMF)
PMF is a receptor model that has been extensively used for source apportionment of airborne particulate matters and VOCs (Lee et al., 1999;Brown et al., 2007). In this study, U.S. EPA PMF 5.0 model (U.S. EPA, 2017) was applied to identify the sources of O 3 precursors, according to Eq. (1) (Paatero, 1997;. where x ij is the measured concentration of the j th species in the ith sample, g ik represents the contribution of the kth source to the ith sample, f kj denotes the fraction of the j th species in the kth source, and e ij is the residual for the j th species in the ith sample. p stands for the total number of independent sources (Paatero, 2000a, b). The uncertainties of the concentrations applied to PMF were set in the same way as Polissar et al. (1998) and Reff et al. (2007). Values below or equal to the LoD were replaced by half of the LoDs and the uncertainties for these values were set as five-sixths of the corresponding LoDs. For the values greater than LoDs, the uncertainties were calculated as [(error fraction × concentration) 2 + (LoD) 2 ] 1/2 , where 10 % was assigned as the error fraction. Missing values (mainly due to maintenance or malfunction of the instruments) were replaced by the geometric mean of the measured values and their accompanying uncertainties were set as 4 times the geometric mean value. More details about the settings of the uncertainty were provided in Norris et al. (2008) and Zhang et al. (2012).
The model was run for 20 times with a random seed, and tests with different number of factors were conducted. The optimum solution was finally determined based on both a good fit to the observed data and the most reasonable and interpretable results according to the knowledge on the sources of O 3 precursors in Hong Kong Ou et al., 2015).

Observation-based model (OBM)
A photochemical box model coupled with the Master Chemical Mechanism (PBM-MCM) was used to simulate the photochemical O 3 formation on the VOC sampling days. In this study, MCM v3.2, a near-explicit chemical mechanism consisting of 5900 species and 16 500 reactions, which fully describes the homogeneous gas phase reactions in the atmosphere (Jenkin et al., 1997Saunders et al., 2003), was used. The observation data of temperature, relative humidity, O 3 , SO 2 , CO, NO, NO 2 and 52 C 2 -C 10 VOCs and OVOCs were input into the model. Specifically, the 52 VOCs and OVOCs included 19 alkanes, 16 alkenes, 13 aromatics and four OVOCs, as shown in Table S5, where the statistics of the mixing ratios of VOCs and OVOCs are also presented. Though previous studies (Guo et al., 2013b;Ling et al., 2016b) indicated that secondary formation dominated the sources of OVOCs in Hong Kong, the primary emissions could not be neglected. Therefore, formaldehyde, acetaldehyde, acetone and propionaldehyde with relatively high abundances were constrained to the observed concentrations in the model, while the other OVOCs with low concentrations and low detection rates were simulated by the model. Nitrous acid (HONO) was not monitored in this study. The average diurnal cycle of HONO mixing ratios measured at the same site in autumn in 2011  was input into the model to roughly represent its role in O 3 formation and atmospheric radical chemistry. Due to the data limitation, the trends of HONO at TC in the three sampling campaigns were not traceable. However, the measurements at a background site in Hong Kong indicated comparable levels of HONO (p > 0.1) between autumn in 2012 and in 2018 (unpublished data). Therefore, adopting the HONO measured in 2011 as the inputs of the simulations in the three sampling campaigns was likely a plausible assumption, despite some uncertainties. The model was also tailored to the real situations in Hong Kong. Specifically, the height of the planetary boundary layer was allowed to vary from 300 m at night to 1400 m at noon. The photolysis rates were calculated according to the measured solar radiations by the Tropospheric Ultraviolet and Visible radiation model (Madronich and Flocke, 1999;, with the detailed method described in Lyu et al. (2017b). In addition to the chemical processes, the exchange between the lower troposphere and free troposphere and dry deposition were also considered in the model. The concentrations of air pollutants in the free troposphere were set according to the observations at a mountainous site in Hong Kong (Lam et al., 2013). The dry deposition rates were adopted from the previous studies (Lam et al., 2013;Saunders et al., 2003). The other physical processes were not included in the model, which might lead to insufficient description of the transport. However, since the model was constrained to the observations, which included the transported air pollutants, the regional transport was partially considered. In addition, the observations at 07:00 on each day were used to initiate each day's modeling, through which the effect of regional transport before the daytime modeling was also considered. We admit that the PBM-MCM cannot perfectly reproduce the real atmospheric processes. However, it performed well in describing the in situ photochemistry in previous studies (Lam et al., 2013;Lyu et al., 2017b;. Actually, the deficiency of PBM-MCM in consideration of the atmospheric dynamics enabled us to assess the contributions of regional transport to O 3 in Hong Kong, based on the differences between the observed and simulated O 3 (Y. .

Simulation scenarios
The PBM-MCM simulates the in situ O 3 photochemistry based on the observed O 3 precursors. Figure S2 shows the average mixing ratios of some O 3 precursors in different wind sectors. The higher levels of CO, ethyne, ethane, propane and toluene under northwest winds indicated the transport of these species from the PRD to Hong Kong. Meanwhile, O 3 might also be transported to Hong Kong. Text S1 in the Supplement discusses the determination of the locally produced and regionally transported O 3 , as well as the uncertainties. Furthermore, to evaluate the contributions of VOC sources to the local O 3 production, two scenarios of model simulation were performed, i.e., scenario A and scenario B. The scenario A simulated the O 3 photochemistry in the whole air, which was constrained by the observed concentrations of all the O 3 precursors. The model simulations in scenario B (including six assumed sub-scenarios) were constrained by the concentrations of O 3 precursors with those contributed by individual sources being subtracted from the observed concentrations. Text S2 elaborates the setup of these scenarios. The simulated O 3 in scenario A was regarded as the locally produced O 3 , as the observed O 3 concentrations were not input to constrain the model. Bear in mind that the regional effects cannot be completely eliminated in this approach, due to the impacts of regional air on the observed concentrations of O 3 precursors. The differences between scenario A and scenario B reflected the contributions of the individual sources to the simulated O 3 production rate. It should be noted that due to the nonlinear relationships between O 3 and its precursors, the subtraction approach only qualitatively rather than quantitatively evaluated the contributions of VOC sources to O 3 production.
3 Results and discussion 3.1 Observation overview  (Cheng et al., 2010;. However, the increases in CO and the notable decrease in NO in 2013 could enhance the O 3 production. The higher O 3 in 2013 indicated that this effect overrode the reduction of TVOCs in influencing the O 3 production. In particular, the decrease in NO meant the reduced NO titration to O 3 , which has been recognized as a primary reason for O 3 increase in VOC-limited regime (Chou et al., 2006;. From 2013 to 2016, the decrease in O 3 was accompanied by the reductions of TVOCs and NO 2 , though CO remained increasing at the same time. NO 2 , as a direct source of O 3 through photolysis, plays an important role in modulating the O 3 variation. Though the causes of NO 2 reduction are unknown to us, it might be one of the critical factors contributing to the decline of O 3 in Hong Kong in recent years. Conversely, the increase in CO was also confirmed by the continuous monitoring data at TC, with a rate of 33.9 ± 0.7 ppbv yr −1 between 2013 and 2016. In fact, the consistent increasing trend (p < 0.05) was also observed at the roadside sites in Hong Kong (not shown here). While the causes of CO increase in Hong Kong may be complicated, the increased vehicle emission is a plausible explanation. Studies (Johnson, 2008;Yao et al., 2008) revealed that while the new engine technologies performed well in reducing NO x emission, they might lead to the increased emission of CO, with the application of a lower air-to-fuel ratio and engine temperature.
In addition, studies have confirmed that continental anticyclones and tropical cyclones are conducive to severe O 3 pollution in Hong Kong because these synoptic systems are often accompanied by northerly winds, high temperature, strong solar radiation, and relatively high pressure in Hong Kong (Ding et al., 2004;Huang et al., 2005;. Table S7 summarizes the number of O 3 episode days with tropical cyclone, continental anticyclone and lowpressure trough in the autumns of 2007, 2013 and 2016. In autumn 2007, eight, eight and one O 3 episode day(s) were found to be related to the tropical cyclone, continental anticyclone and low-pressure trough, respectively, with two O 3 episode days under the combined influence of tropical cy- clone and continental anticyclone. There were also 11 and five O 3 episode days in association with tropical cyclone and continental anticyclone in autumn 2013, respectively (H. . However, four out of the five episode days found in autumn 2016 were associated with tropical cyclone, with the other one related to low-pressure trough. Therefore, the lower O 3 and fewer O 3 episode days in 2016 were also benefited by the meteorological conditions. Figure 3 compares the simulated O 3 in scenario A and the observed O 3 on the VOC sampling days. Overall, both the magnitudes and the temporal patterns of the observed O 3 were reasonably reproduced, though the mean of the simulated O 3 (33.8 ± 1.9 ppbv) was slightly lower than the observed average (37.1 ± 2.0 ppbv). To quantitatively evaluate the model performance, the index of agreement (IOA) was used to examine the goodness of fit between simulated and observed O 3 . Within the range of 0-1, higher IOA represents better agreement between the simulated and observed values (Willmott, 1982). In this study, the overall IOA for the three sampling periods was 0.74, within the range of IOA (0.67-0.89) accepted by previous studies (Lyu et al., 2015. Good correlations (R 2 = 0.61) were also shown between the simulated and observed hourly O 3 . Bearing in mind the deficiencies of the box model in describing the atmospheric dynamics, we believed that the modeling results were acceptable, but special attention and explanation to the discrepancies between the simulated and observed O 3 were needed.

Model validation
It was found that the discrepancies were most likely caused by the transport processes, i.e., vertical and horizontal trans-port, which were not fully represented in the PBM-MCM model (George et al., 2013;Lakey et al., 2015;. For example, the simulated O 3 (maximum: 122.6 ppbv) was much higher than the observed O 3 (maximum: 44.3 ppbv) on 16 November 2007, when the strong southeast winds (wind direction: 90-180 • ) with the highest wind speed of 5.3 m s −1 prevailed in Hong Kong. The south sector winds from SCS might dilute the locally produced O 3 and the O 3 precursors/intermediates (such as the radicals), which were not constrained by the observations. The same circumstances were also observed on 27 October and 17 November 2007 and 11-12 September and 20 November 2013, with southeast winds dominating (74.4 %) during the daytime (Fig. 3). For those days with the simulated O 3 lower than the observed O 3 , i.e., 3 and 22-25 October 2013 and 6 November 2016, 69.3 % of the winds during the daytime came from the north (wind directions: 0-90 and 270-360 • ), which might transport the air masses laden with O 3 and/or O 3 precursors/intermediates not constrained to the observations from the inland PRD to the sampling site. The observed O 3 mixing ratios are plotted against the wind fields in Fig. S4. It is obvious that O 3 values were higher under the north winds, while they were lower in the south wind sectors, confirming the effects of dilution and regional transport of the south and north winds on O 3 pollution in Hong Kong, respectively.

Inter-annual variations in the locally produced and regional transported O 3
As discussed in Sect. 2.5, the simulated O 3 in scenario A could be regarded as the locally produced O 3 . Therefore, the differences between the observed O 3 and simulated O 3 in scenario A were treated as the regionally transported O 3 (Y. . It is noteworthy that some negative values were generated with this method, corresponding to  the dilution of the south winds to the locally produced O 3 as elaborated in Sect. 3.2.1. Figure 4 shows the hourly mixing ratios of the observed, local and regional O 3 at TC in daytime hours (07:00-19:00 LT) of the three sampling campaigns. Overall, the observed O 3 was mainly (88.7 ± 2.5 %) contributed by the local photochemical production, with regional transport only accounting for 11.3 ± 2.5 % of the observed daily maximum O 3 . However, regional transport was responsible for as high as 58.0 ± 5.4 % of the observed daily maximum O 3 in Hong Kong on the O 3 episode days when northerly winds prevailed, indicating the heavy O 3 burden superimposed by regional air masses from PRD. From 2007 to 2013, the simulated locally produced O 3 remained statistically unchanged (p > 0.1), in contrast to the increase in observed O 3 and regional O 3 at rates of 1.78 ± 0.05 ppbv yr −1 (p < 0.05) and 1.77±0.04 ppbv yr −1 (p < 0.05), respectively, the same trends as those reported by Y.    (Xue et al., 2014a;Lyu et al., 2017a;, and would be further evaluated in this study (Sect. 3.4). In addition, the emission controls in mainland China might contribute to the decrease in O 3 or at least lessen the regional O 3 burden in this period. For example, China's NO x emissions for the first time showed a decreasing trend from 2013, benefiting from the implementation of China's Clean Air Action Plan (Zheng et al., 2018). Furthermore, we looked into the monthly average O 3 observed at the 12 air quality monitoring stations across the inland PRD, including three regional monitoring stations, i.e., Tianhu, Wanqingsha and Jinguowan, and nine urban monitoring stations, i.e., Xiapu, Jinjuzui, Donghu, Tangjia, Liyuan, Huijingcheng, Zimaling, Luhu and Chengzhong (https://www.epd.gov.hk/epd/sc_chi/resources_ pub/publications/m_report.html, last access: 10 April 2019). As shown in Fig. 5, O 3 at these stations remained relatively stable (p = 0.68) during 2006-2013, which however showed a contrastively decreasing trend at a rate of −1.27 ± 0.25 ppbv yr −1 from 2013 to 2016. This corroborated our modeling results that the regional contribution to O 3 in Hong Kong ceased increasing or even began to decrease since 2013. While the substantial decrease in NO x was a plausible reason for the alleviated regional O 3 pollution, meteorological variations might also play roles in modulat-ing O 3 variations in these years (Li et al., 2019). However, analyses on the causes are out of the scope of this study. In addition to the reduced local formation and regional transport of O 3 , the more favorable meteorological conditions in 2016 might be another reason for the O 3 decrease, as discussed in Sect. 3.1.

Local production and destruction pathways of O 3 and the OH radical
3.3.1 In situ net O 3 production Figure 6 shows the average diurnal profiles of the simulated O 3 production and destruction pathways during the three sampling campaigns. Also shown are the average diurnal cycles of the simulated O 3 . The shift of the peaks between the net O 3 production rate and the simulated O 3 was due to the accumulation of the newly generated O 3 over time in the model, which was also true in the real situations. The re- actions between NO 2 and O 3 , leading to the formation of NO 3 and N 2 O 5 , in addition to dry deposition and aloft exchange, were the main depletions of the simulated O 3 in the late afternoon. Consistent with previous studies (Kanaya et al., 2009;Liu et al., 2012;Xue et al., 2014b), these pathways were not included in the calculation of the net O 3 production rate because we mainly focused on the photochemical processes in the hours when O 3 was accumulated. It was found that the reaction between HO 2 with NO dominated the O 3 production rates in all the cases, with an average rate of 3.7± 0.7 ppbv h −1 (56.5 ± 1.1 %, percentage of the total O 3 production rate, same below), 2.5 ± 0.3 ppbv h −1 (64.3 ± 0.8 %) and 1.4 ± 0.2 ppbv h −1 (67.7 ± 0.7 %) in the 2007, 2013 and 2016 sampling campaigns, respectively. In addition, the sum of the reaction rates between RO 2 radicals and NO contributed 3.0±0.6 ppbv h −1 (43.5±1.1 %), 1.5±0.2 ppbv h −1 (35.7 ± 0.8 %) and 0.7 ± 0.1 ppbv h −1 (32.3 ± 0.7 %) to the O 3 production rate in 2007, 2013 and 2016, respectively. The formation of HNO 3 though the reaction between OH and NO 2 served as the main scavenger pathway of O 3 , as NO 2 would be photolyzed and produce O 3 otherwise. On average, O 3 was consumed in this way at a rate of −1.3 ± 0.2 ppbv h −1 (80.7 ± 3.3 %, percentage of the total O 3 destruction rate, same below), −1.0 ± 0.1 ppbv h −1 (79.3 ± 1.8 %) and −0.6 ± 0.07 ppbv h −1 (81.6 ± 2.0 %) in 2007, 2013 and 2016, respectively. The photolysis of O 3 was the second contributor to O 3 destruction, with an average contribution of −0.11 ± 0.01 ppbv h −1 (8.5 ± 0.5 %) for the three sampling periods. In addition, the ozonolysis of unsaturated VOCs and the reactions between O 3 and radicals (OH and HO 2 ) were responsible for 3.5 ± 0.3 % and 1.7 ± 0.2 % of the total destruction rate of the locally produced O 3 , respectively.
Overall, the net local O 3 production rate decreased from 5.2 ± 1.

Recycling of OH radical
As one of the most important radicals in the atmosphere, OH initiates the oxidation of VOCs, leading to O 3 formation. Figure 7 presents the average diurnal profiles of the simulated OH and the formation and loss pathways dominating the recycling of OH during the three sampling periods, which roughly followed the typical pattern of the intensities of photochemical reactions, i.e., higher at noon and lower at the beginning and end of the day. On average, the simulated OH concentration was comparable (p = 0.4) between the 2007 sampling campaign (1.6±0.3×10 6 molecules cm −3 ) and the 2013 sampling campaign (1.5 ± 0.2 × 10 6 molecules cm −3 ), but it decreased (p < 0.05) to 1.0±0.2×10 6 molecules cm −3 in the 2016 sampling campaign.
As expected, the formation and loss rates of OH were basically balanced in all the cases. OH was mainly formed from the reaction of HO 2 + NO, which accounted for 69.8 ± 1.1 % of the total OH production rate over the three sampling campaigns. The photolysis of HONO ranked second in supplying OH with the contribution of 22.0 ± 1.4 %. As stated in Sect. 2.4.2, the average diurnal cycle of HONO measured at TC in 2011 was adopted in the simulations. To assess the uncertainties, we also calculated the HONO concentra- tions according to the measured HONO/NO x ratios and the NO x concentrations at TC in the three sampling campaigns (Fig. S5). The uncertainties in HONO concentrations and in the contributions of HONO to OH formation and loss rates are discussed in Text S3. The formation of OH from HONO photolysis was most efficient in the early morning, which was explained by the morning peak of HONO concentration, due to the nocturnal heterogeneous formation and the vehicle emissions in morning rush hours. Apart from the two dominant pathways, O 3 photolysis (6.3 ± 0.2 %), ozonolysis of unsaturated VOCs (1.5 ± 0.2 %) and H 2 O 2 photolysis (0.2±0.01 %) also made some contributions to the formation of OH, with the highest rates at noon or in the early afternoon when the productions of O 3 and H 2 O 2 were the most intensive. To sum up, the total formation rates of OH from the primary sources (photolysis of HONO, O 3 and H 2 O 2 , and ozonolysis of VOCs) were lower than the recycling rates of OH (HO 2 + NO) throughout the day at TC, consistent with the results in Xue et al. (2016) simulated at the same site. The dominant role of HO 2 + NO in OH formation at TC (average contribution of 69.8 ± 1.1 %) might be related to the abundant NO at this site. The same pathway was simulated and accounted for only 42.7 ± 0.2 % of the total OH formation rate at an island more than 40 km away from Hong Kong with very low NO concentrations, i.e., maximum of 0.56 ppbv (Y. . OH was mainly depleted by the reactions with VOCs (32.3 ± 1.2 %), NO 2 (31.9 ± 0.9 %), CO (19.3 ± 0.6 %) and NO (16.5 ± 1.1 %). The reaction rates of OH + NO (formation rates of HONO) had the highest values in the morning, approximately in line with the diurnal pattern of the HONO photolysis rates, which however were not completely balanced due to the constraint of HONO to observations in the model. The average net photolysis rates of HONO (differences between the HONO photolysis and formation rates) were 0.68 ± 0.21 × 10 6 , 0.70 ± 0.12 × 10 6 and 0.87 ± 0.12 × 10 6 molecules cm −3 s −1 in the 2007, 2013 and 2016 sampling campaigns, respectively. The losses of OH through the other pathways all exhibited the highest efficiencies at noon or in the early afternoon. It should be noted that the reaction between OH and NO 2 was not only the sink of OH but also a termination reaction in the photochemical system. In comparison, the termination reaction rates were lower than the OH formation rates from the primary sources (photolysis of HONO, O 3 and H 2 O 2 , and ozonolysis of VOCs) in the morning (07:00-10:00 LT), which were reversed in the following hours of the day due to the increases in OH concentrations.
Consistent with the variations in the local O 3 production, both the local formation and loss rates of OH decreased from 2007 to 2016 (p < 0.05), with much more obvious reductions in the later phase (2013)(2014)(2015)(2016). On the one hand, the continuous reduction of VOCs resulted in lower HO 2 and RO 2 concentrations (Fig. S6), hence the lower production rate of OH through the reaction of HO 2 + NO. At the same time, the destruction rates of OH also decreased due to the reductions of OH and the O 3 precursors, except for CO ( Fig. 7 and Table 1). The decreases in the OH production and destruction rates indicated that the propagation of the reaction cycles, namely the recycling of OH, became slower from 2007 to 2016. This also explained why the locally produced O 3 decreased in these 10 years, since O 3 is formed with the consumption and recycling of OH radical.

Source apportionment
To resolve the sources of O 3 precursors, 27 species, including CO, NO, NO 2 , 12 alkanes, four alkenes and eight aromatics, were applied to PMF for source apportionment. These species were of either high abundances or typical tracers of VOC sources in Hong Kong. Source apportionment was conducted for a total of 414 samples covering the three sampling periods, so that the uncertainty of the source apportionment results could be reduced, compared to separate source apportionments for each of the three sampling periods. Figure 8 shows the average profiles of the six sources resolved by PMF. The modeling errors were estimated with the bootstrap method integrated in PMF (Brown et al., 2015). Factor 1 was assigned as the combination of LPG usage and gasoline evaporation, in view of the high loadings of C 2 -C 5 hydrocarbons. Specifically, propane and i-/n-butanes are the main components of LPG in Hong Kong, and gasoline evaporation generally contains large quantities of i-/npentanes, in particularly i-pentane (Guo et al., 2013a;Lyu et al., 2017a). Factor 2 was characterized by moderate to high percentages of i-/n-pentanes and TEX (toluene, ethylbenzene and xylenes). These species are commonly seen in gasoline exhausts. Therefore, we defined this factor as gasoline exhausts. Both the third and fourth factors indicated solventrelated emissions. While Factor 3 likely represented household solvent usage, due to the dominance of hexane and hexane isomer (3-methylpentane) Ou et al., 2015), Factor 4 was more related to emissions from coatings and paints, in view of the dominance of the aromatics . Factor 5 was distinguished by the high concentrations of ethane, ethene, ethyne and benzene, together with the relatively heavy (C 7 -C 10 ) alkanes, which are typical species in diesel exhausts (Schauer et al., 1999;Kashdan et al., 2008;Sahoo et al., 2011). Therefore, this factor was designated as diesel exhausts. The last factor was denoted as biogenic emissions (BVOCs), due to the exclusive dominance of isoprene (Guenther, 2006). Figure S7 presents the total mixing ratio of VOCs emitted from each individual source extracted from PMF during the three sampling periods in Hong Kong. The VOC emissions from LPG usage and gasoline evaporation decreased significantly (p < 0.05) at a rate of −2.61 ± 0.03 ppbv yr −1 from 2007 to 2016. However, the VOCs in association with gasoline exhausts experienced an increase (rate = 1.32 ± 0.02 ppbv yr −1 , p < 0.05) in these years, indicating that the reduction of VOC emissions from LPG usage and gasoline evaporation was not attributable to the change in emissions of gasoline-fueled vehicles. Insight into the mixing ratios of propane and i-/n-butanes (LPG tracers) in this source revealed a significant decline from 3.51 ± 0.52 ppbv in the 2007 sampling campaign to 1.27 ± 0.11 ppbv in the 2016 sampling campaign. Therefore, the reduction of VOC emissions from LPG usage was most likely the reason of the decrease in VOCs allocated to the source of LPG usage and gasoline evaporation. In fact, it was confirmed by our previous studies (Lyu et al., 2016b;Yao et al., 2019) that the replacement of catalytic converters on LPG-fueled vehicles during September 2013-May 2014 effectively reduced the VOC emissions from LPG-fueled vehicles in Hong Kong. In addition, the variations in LPG usage in the inland PRD, where LPG was extensively used as vehicular and domestic fuels (Liu et al., 2008), might also contribute to the emission reduction of VOCs, in view of the decrease in LPG tracers in this source from 2007 (3.51 ± 0.52 ppbv) to 2013 (2.04 ± 0.27 ppbv), when no control was performed against LPG-fueled vehicle emissions in Hong Kong. The VOCs emitted from solvent usage (including the household solvent, coatings and paints) also decreased significantly (p < 0.05) from 2007 to 2016, likely benefiting from the actions taken to restrict the VOC contents in solvent products starting from 2007 (phase I) and 2010 (phase II) in Hong Kong (Lyu et al., 2017a). VOCs attributable to diesel exhausts decreased (p < 0.05) from the 2007 (2.6 ± 0.3 ppbv) to 2013 sampling campaign (2.0 ± 0.2 ppbv), which however were unchanged between 2013 and 2016 (2.2 ± 0.2 ppbv). In fact, a subsidy program has been implemented in Hong Kong since 2007 to progressively eliminate the pre-Euro IV diesel vehicles or to upgrade their emission standards to Euro IV (HKEPD, 2017b), and the effectiveness of this program in VOC reductions until 2013 was confirmed by Lyu et al. (2017a) with the online measurement data at the same site. However, while phase III of this program (2014-2019) is still ongoing, the VOCs emitted from diesel vehicles remained stable between the 2013 and 2016 sampling campaigns. This undesirable result might be due to the fact that the actions were mainly targeted at the pre-Euro, Euro I and Euro II diesel vehicles before 2013, whereas phase III of the program initiated in 2014 focused on the Euro III vehicles (HKEPD, 2017b(HKEPD, , 2018. Since the former were vehicles with higher emissions, it is not unreasonable that reduction of VOCs was more discernible between 2007 and 2013. Further, the effectiveness of the phase III program might be somewhat offset by the wearing-out of the pre-existing vehicles and the in-crease in diesel vehicle populations (Competition Commission, 2017). Further evaluation with more data in a longer period is recommended. At last, the increase in BVOCs from 2007 to 2013 but comparable levels between 2013 and 2016 seemed to be related to the lower (p < 0.05) temperature in the 2007 sampling campaign ( Fig. S8 and Table S6). In addition, the more frequent (62.8 %) southeast winds from SCS with higher wind speeds (2.3 ± 0.2 m s −1 ) might dilute BVOCs emitted from the terrestrial plants in the 2007 sampling campaign. Figure 9 presents the contributions of VOCs emitted from individual sources to the production and destruction rates of O 3 , as well as the simulated contributions to the O 3 mixing ratios. NO x was not included in these analyses because of its relatively high uncertainties in source apportionment results due to the short lifetimes. Consistent with the O 3 production and destruction in the whole air, the pathway of HO 2 + NO dominated over the reactions between RO 2 and NO in O 3 production for all the individual sources. The destruction of O 3 was mainly driven by NO 2 reacting with OH. For the net O 3 production rate, VOCs attributable to the coatings and paints made the largest contribution (0.38 ± 0.05 ppbv h −1 ), followed by gasoline exhausts (0.22 ± 0.03 ppbv h −1 ), LPG and gasoline evapora- tion (0.21 ± 0.03 ppbv h −1 ), BVOCs (0.19 ± 0.03 ppbv h −1 ), household solvent usage (0.15 ± 0.04 ppbv h −1 ) and diesel exhausts (0.13±0.01 ppbv h −1 ). Despite some peak shifts for the reasons illustrated in Sect. 3.3.1, the O 3 mixing ratios elevated by the individual sources followed the same pattern as the net O 3 production rates, with the highest O 3 enhancement (1.92 ± 0.21 ppbv) from the source of coatings and paints and the lowest increase from household solvent usage (0.86±0.06 ppbv) and diesel exhausts (0.83±0.06 ppbv). The contributions of source-specific VOCs to O 3 production, particularly the importance of solvent usage in O 3 formation in Hong Kong, were generally in line with previous studies Ou et al., 2015). This was actually expected according to the reactivity of major VOCs in each source. For example, the TEX in the source of coatings and paints (Fig. 8) has been identified to be of high O 3 formation potentials (Lau et al., 2010;Ling et al., 2011Ling et al., , 2013.

Source contributions to O 3 production
However, the PBM-MCM model simulations enabled us to quantitatively evaluate the contributions of VOC sources to O 3 production rates.
From a historical perspective, we found that the contribution of LPG usage and gasoline evaporation to O 3 production significantly decreased (p < 0.05) from the 2007 to 2016 sampling campaign (2007: 0.51±0.11 ppbv h −1 ; 2013: 0.20±0.03 ppbv h −1 ; 2016: 0.10±0.02 ppbv h −1 ), which coincided with the variations in VOCs emitted from LPGfueled vehicles as discussed above. Gasoline exhaust contributed much less (p < 0.05) to the net O 3 production rate in 2007 (0.02 ± 0.01 ppbv h −1 ) than those in 2013 (0.26 ± 0.05 ppbv h −1 ) and 2016 (0.27 ± 0.07 ppbv h −1 ), in line with the variations in VOCs emitted from this source. The reductions of VOC emissions from solvents also resulted in the consistent decrease in the net O 3 production rate from 1.22 ± 0.17 ppbv h −1 in the 2007 to 0.14 ± 0.05 ppbv h −1 in the 2016 sampling campaign. The O 3 production rates contributed by VOCs in diesel exhausts were reduced from 2007 (0.21±0.05 ppbv h −1 ) to 2013 (0.11±0.02 ppbv h −1 ) and remained unchanged thereafter (2016: 0.11 ± 0.02 ppbv h −1 ). The O 3 production rate traceable to BVOCs showed a significant increase from 2007 (0.04 ± 0.02 ppbv h −1 ) to 2016 (0.22±0.04 ppbv h −1 ) since the mixing ratios of BVOCs significantly increased (p < 0.05) in these years. It is noteworthy that the changes in meteorological conditions in these three sampling campaigns might also partially account for the variations in the source contributions to O 3 production. For example, the 2013 sampling campaign was characterized by the relatively higher temperature and lowest relative humidity among the three sampling periods, which favored O 3 formation in 2013 (Table S6). In addition, due to limited samples in this study, we recommend further assessments with more data in longer periods to be carried out in future study.

Conclusions
Photochemical pollution with high and increasing concentrations of O 3 has been an important environmental issue in south China. With the observation data of O 3 and its precursors at a suburban site in Hong Kong, downwind of south China, this study analyzed the inter-annual variations in O 3 and its photochemistry, as well as the contributions of VOC sources to the local O 3 production rates in 2007, 2013 and 2016. To our knowledge, this is the first time that a substantial alleviation of O 3 pollution in this region was identified between 2013 and 2016, in contrast to the repeatedly confirmed O 3 increase before 2013. In addition to the changes in meteorological conditions among the three sampling campaigns, the termination of the rise in regionally transported O 3 and the decrease in the local O 3 production rate contributed to the decline of O 3 in the later period. The emission reductions (particularly for NO x ) in mainland China starting from 2013, the year when China's Clean Air Action Plan was launched, might more or less play a role in ceasing the increase in regional O 3 . In Hong Kong, the replacement of catalytic converters and the constraints of VOC contents in solvent products led to the reductions of VOC emissions from LPG-fueled vehicles and solvent usage, respectively. As a result, the local O 3 production rate and the recycling rate of the OH radical decreased substantially from 2013 to 2016. Though the variations in meteorological conditions and the limited sample size might somewhat introduce uncertainties to the conclusions drawn from the present study, it is plausible that the local and regional interventions were effective on the control of O 3 pollution in Hong Kong. Nevertheless, studies with more data in longer periods should be conducted, not only in Hong Kong but also in mainland China where O 3 is still increasing in most of the territories. Data availability. The data are accessible at https: //drive.google.com/file/d/1rPH8TQMGvnGKjAwTYcNveqtCi_ zo019R/view?usp=sharing (Liu, 2019) Supplement. The supplement related to this article is available online at: https://doi.org/10.5194/acp-19-5127-2019-supplement.
Author contributions. HG and FJ initiated and designed the experiments, and XuL and XiL carried them out. XiL and YW developed the model code and performed the simulations. XuL and XiL prepared the paper and HG finalized the paper with contributions from all co-authors. XuL and XiL are contributed equally to this work.
Competing interests. The authors declare that they have no conflict of interest.