Tracer-based source apportioning of atmospheric organic carbon and 
the influence of anthropogenic emissions on secondary organic aerosol 
formation in Hong Kong

Abstract. Here we conducted comprehensive chemical characterization and source apportionment of 49 PM2.5 samples collected in Hong Kong. Besides the major aerosol constituents, 39 polar organic species, including 14 secondary organic aerosol (SOA) tracers of isoprene, monoterpenes, β-caryophyllene, and naphthalene, were quantified using gas chromatography-mass spectrometry (GC-MS). Six factors, i.e., SOA, secondary sulfate (SS), biomass burning (BB)/SOA, sea salt, marine vessels, and vehicle emissions, were apportioned by positive matrix factorization (PMF) as the major sources of ambient organic carbon (OC) in Hong Kong. The secondary formation, including OC from SOA, SS, and aging of BB plume, was the leading contributor to OC (51.4 %, 2.15 ± 1.37 μgC m−3) throughout the year. We then applied a tracer-based method (TBM) to estimate the SOA formation from the photo-oxidation of four selected precursors, and monoterpenes SOA was the most abundant. A Kintecus kinetic model was used to examine the formation channels of isoprene SOA, and the aerosol-phase ring-opening reaction of isoprene epoxydiols (IEPOX) was found to be the dominant formation pathway. Consistently, IEPOX tracers contributed 94 % of total GC-MS quantified isoprene SOA tracers. The TBM-estimated SOC (SOCTBM) and PMF-resolved SOC (SOCPMF) showed similar temporal trends; however, SOCTBM only accounted for 26.5 % of SOCPMF, indicating a large fraction of ambient SOA was from other reaction pathways/precursors. Results of Pearson's R and multivariate linear regression analysis showed that NOx processing played a key role in both daytime and nighttime SOA production in the region. SOA formation through nighttime NO3 oxidation of biogenic VOCs, especially monoterpenes, may have made a considerable contribution to the SOA loading in Hong Kong. Moreover, sulfate had a significant positive linear relationship with SOCPMF and SS-related SOC, and particle acidity was significantly correlated with SOC from the aging of BB.


Abstract. Here we conducted comprehensive chemical characterization and source apportionment of 49 PM 2.5 samples collected in Hong Kong. Besides the major aerosol constituents, 39 polar organic species, including 14 secondary organic aerosol (SOA) tracers of isoprene, monoterpenes, β-caryophyllene, and naphthalene, were quantified using gas chromatography-mass spectrometry (GC-MS). Six factors, i.e., SOA, secondary sulfate (SS), biomass burning (BB)/SOA, sea salt, marine vessels, and vehicle emissions, were apportioned by positive matrix factorization (PMF) as the major sources of ambient organic carbon (OC) in Hong Kong. The secondary formation, including OC from SOA, SS, and aging of BB plume, was the leading contributor to OC (51.4 %, 2.15 ± 1.37 µg C m −3 ) throughout the year. We then applied a tracer-based method (TBM) to estimate the SOA formation from the photo-oxidation of four selected precursors, and monoterpene SOA was the most abundant. A Kintecus kinetic model was used to examine the formation channels of isoprene SOA, and the aerosol-phase ring-opening reaction of isoprene epoxydiols (IEPOXs) was found to be the dominant formation pathway. Consistently, IEPOX tracers contributed 94 % to total GC-MS-quantified isoprene SOA tracers. The TBM-estimated secondary organic carbon (SOC TBM ) and PMF-apportioned SOC (SOC PMF ) showed similar temporal trends; however, SOC TBM only accounted for 26.5 % of SOC PMF , indicating a large fraction of ambient SOA was from other reaction pathways or precursors. Results of Pearson's R and multivari-ate linear regression analysis showed that NO x processing played a key role in both daytime and nighttime SOA production in the region. Moreover, sulfate had a significant positive linear relationship with SOC PMF and SS-related SOC, and particle acidity was significantly correlated with SOC from the aging of BB. more targeted control policies of PM 2.5 pollution (Hu et al., 2010;Huang et al., 2014;Schauer et al., 2007;Simoneit, 1999;Zheng et al., 2005). Huang et al. (2014) applied positive matrix factorization (PMF) to apportion the sources of OA at four urban locations in China, i.e., Beijing, Shanghai, Guangzhou, and Xi'an. They found that secondary formation accounted for a predominant fraction of OA (44 %-71 %) at all four sites. Hong Kong, a megacity located on the southern coast of China in the Pearl River Delta (PRD) region and a hub port for the southern Asia-Pacific region, has its unique organic carbon (OC) source characteristics. Hu et al. (2010) incorporated biogenic and anthropogenic SOA tracers and some primary organic aerosol (POA) markers into PMF and identified seven OA sources in Hong Kong. They found that 45 % of OC in Hong Kong during summertime was from secondary formation, and the number could reach up to 65 % on sampling days under regional pollution from the PRD area.
All these studies have illustrated the importance of secondary formation to OA in the ambient atmosphere. However, due to SOA's complex chemical composition and formation mechanisms, a precise prediction of SOA load from individual precursors at both regional and global scales is still challenging. An SOA tracer-based method (TBM) has been developed to partially solve this problem, which estimates the amount of SOA and secondary organic carbon (SOC) formed from the atmospheric oxidation of selected VOCs (i.e., isoprene, monoterpenes, β-caryophyllene, toluene, and naphthalene) using the mass ratios of tracer-to-SOA/SOC obtained from laboratory smog chamber experiments . However, TBM can only capture SOC formation from the above-listed VOC precursors, and it may underestimate the actual SOC levels in the ambient atmosphere due to the lack of SOA tracer-to-SOC ratio values of a broader range of OA precursors. Therefore, besides the SOA tracer-based method, we have also applied PMF to evaluate the contributions of SOC and primary emissions to OA in the region.
Many studies have reported the observational evidence of biogenic SOA enhancement induced by anthropogenic emissions, such as nitrogen oxides (NO x ) and sulfur dioxide (SO 2 ) (Huang et al., 2014;Xu et al., 2015;Rattanavaraha et al., 2016). NO x is one of the critical drivers of SOA formation through the photochemical oxidation of VOCs via peroxy radical pathways (Finlayson-Pitts and Pitts, 2000). Nitrogen dioxide (NO 2 ) reacts with ozone (O 3 ) to form an NO 3 radical, a critical nighttime gas oxidant. Several laboratory studies have reported high SOA yields from the oxidation of biogenic VOCs (BVOCs) by NO 3 radical (Fry et al., 2009;Ng et al., 2008). Some field studies also revealed that SOA formation from NO 3 oxidation of BVOCs occurs during both daytime and nighttime (Brown et al., 2013;Rollins et al., 2013). The effect of SO 2 on SOA formation was often explained in the context of particle acidity in laboratory studies, which promotes SOA production through acid-catalyzed heterogeneous reactions (Jang et al., 2002;Surratt et al., 2010). Sulfate was also suggested to enhance isoprene-SOA formation by acting as the nucleophiles, providing active aerosol surface area, and through the salting-in effect (Rattanavaraha et al., 2016;Xu et al., 2015). Recently,  proposed a new sulfate formation pathway in aqueous aerosols through NO 2 oxidation and ammonium neutralization, and synchronous enhancements of both nitrate and SOA production in aqueous aerosols were reported. These laboratory and field monitoring studies have shown that the abundance and chemical nature of ambient OA are significantly influenced by the complex interactions among source emissions, anthropogenic activities, atmospheric physical or chemical processes, and meteorological conditions .
In this study, we collected 49 PM 2.5 samples at an urban site in Hong Kong during a whole year period. Concentration levels of 39 polar organic species were quantified using gas chromatography-mass spectrometry (GC-MS), and their temporal and meteorological variations were evaluated. With the input of SOA tracers and primary source markers into PMF, we quantitatively assessed the contributions of various primary and secondary sources to OC. SOC formed from individual biogenic (i.e., isoprenes, monoterpenes, and βcaryophyllene) and anthropogenic VOCs (i.e., naphthalene) were estimated using the TBM. Finally, the impacts of anthropogenic pollutants (e.g., NO 2 , O 3 , NO 3 , SO 2 , and tropospheric odd oxygen (O x )) and PM 2.5 constituents (e.g., sulfate, acidity, and liquid water content) on total and individual SOCs estimated by both TBM and PMF were evaluated using Pearson's R analysis and multilinear regression model. This study provides comprehensive information on the sources of OA and SOA in Hong Kong as well as direct evidence of anthropogenic influences on the SOA formation in the region, which may serve as the scientific basis for the formulation of the PM 2.5 mitigation policy in the region.

Sample collection
The PM 2.5 samples were collected on the 12th floor of Science Tower in the campus of Hong Kong Baptist University (114 • 15 E, 22 • 13 N; ∼ 40 m above the ground), which is a typical urban site. PM 2.5 samples were collected from 6 September 2011 to 16 August 2012, and a total of 49 samples were collected. A high-volume air sampler was used to collect PM 2.5 onto a quartz fiber filter (20 cm × 25 cm) at a flow rate of 1.13 m 3 min −1 for 24 h. The quartz fiber filters were prebaked at 550 • C for 24 h to remove organic contaminants. After sampling, the filters were immediately transferred to the laboratory and stored at −18 • C until analysis. 1 % trimethylchlorosilane, TMCS). For each aerosol sample, 20 cm 2 of the filter was cut into small pieces and sonicated for 10 min with 10 mL of distilled acetonitrile (HPLC grade); the extraction was repeated three times. The extracts were combined and filtered through a Millipore 0.45µm PTEE hydrophobic Teflon filter into a 50 mL round flask, concentrated to ∼ 0.5 mL by rotary evaporation, and transferred into a 5 mL reaction vial. The round flask was rinsed with 1 mL of acetonitrile three times, and the rinsing solvent was transferred into the reaction vial as well. The final extract was blown to dryness under a gentle stream of pure nitrogen gas at 40 • C and then derivatized with 100 µL of BSTFA and 50 µL of pyridine at 70 • C for 2 h. After the reaction vial cooled down to room temperature, 30 µL of tetracosane-d 50 (internal standard, 50 µg mL −1 in hexane) was added. The derivatives were analyzed by GC-MS. Two microliters of the derivatized sample or standard were injected and separated on an HP-5MS capillary column (30.0 m×250 µm×0.25 µm, Agilent J&W). The temperature program and instrument settings were adapted from the method used by Hu et al. (2008).
Saccharides, di-and tricarboxylic acids, 4-nitrocatechol, and cholesterol were identified and quantified using authentic standards. The SOA tracers were identified using surrogate compounds with similar structures and functional groups (Hu et al., 2008;Hu and Yu, 2013), and the detailed information is provided in Table 1. Recovery tests of these organic species were carried out by spiking the mixture of standards onto blank quartz filters, followed by the same sample extraction and analysis processes. Recoveries of the polar compounds were within the range of 80 % to 120 %. Analysis of hopanes has been reported in our previous study (Ma et al., 2019). Four hopanes, including 17α,21β-hopane, 17α,21β-22R-homhopane, 17α,21β-22S-homhopane, and 17α,21β-30-norhopane, were measured using an Agilent 6890N-5975 GC-MS with the thermal desorption (TD) method. Recoveries of four hopane standards ranged from 83 % to 98 %.

Results
Hourly meteorological and air quality data (i.e., temperature, relative humidity (RH), O 3 , SO 2 , and NO 2 ) in the vicinity of the sampling site were collected by Hong Kong Environmental Protection Department (HKEPD) (http://envf.ust.hk/ dataview/gts/current/, last access: 16 June 2021). During the sampling period, the ambient temperature ranged from 14.52 to 31.01 • C, with an annual average of 24.17 ± 5.00 • C. The daily average of RH ranged from 52.94 % to 97.02 %, with a yearly average of 79.87 ± 10.54 %. Heavy rains are common in Hong Kong during summer, which effectively washes out the PM pollutants.
Hong Kong is located at the southeast edge of the Pearl River Delta (PRD) region. PRD is a rapidly developing area with intensive industrial activities. As described in our previous studies on the analysis of humic-like substances (HULISs) and water-soluble PM 2.5 -induced oxidative potential using the same set of PM 2.5 samples (Ma et al., 2019;Cheng et al., 2021), we carefully examined the air mass backward trajectories and categorized all sampling days into three groups, i.e., days mainly influenced by the regional pollution from the PRD region (regional days), days influenced by long-range regional transport of air mass from the northern and eastern China (LRT days), and days dominated by the locally generated pollutants (local days). The concentration levels of both PM 2.5 and O 3 and the spatial distribution patterns of SO 2 over the 18 Hong Kong air quality monitoring stations (http://envf.ust.hk/dataview/gts/current/, last access: 16 June 2021) on each sampling day were also checked to assist the classification. A summary of the classification of sampling days and the typical air mass backward trajectories under the three meteorological categories are presented in Table S2 and Fig. S2 in Ma et al. (2019), respectively.
The gas pollutants, i.e., O 3 , NO 2 , O x , and SO 2 , showed significantly higher average concentrations on regional days than those on LRT local days ( Table 2). The annual mean concentrations of O 3 , NO 2 , and SO 2 were 14.85 ± 8.69 ppb, 37.15 ± 9.76 ppb, and 4.45 ± 2.57 µg m −3 , respectively. Given that O 3 and NO 2 undergo a rapid photochemical conversion in the ambient atmosphere, the tropospheric odd oxygen O x (the sum of O 3 and NO 2 ) was calculated as an indicator of atmospheric oxidation capacity. As shown in Table 1, O x ranged from 49.72 to 145.90 µg m −3 , with a mean value of 99.31 ± 27.42 µg m −3 , indicating a high oxidation capacity of the Hong Kong atmosphere. The annual average concentrations of OC and EC were 4.18 ± 2.37 and 1.02 ± 0.54 µg C m −3 , respectively. Ambient OC levels observed on regional days (6.15 ± 2.51 µg C m −3 ) were about 2 times higher than those on LRT and local days; as for EC, it exhibited relative constant levels throughout the year (0.14-2.75 µg C m −3 ). This confirms that EC is mainly emitted locally in Hong Kong, and OC has some regional sources. Moreover, the OC/EC ratios of the collected samples ranged from 1.51 to 10.91, with an annual average value of 4.61, indicating that secondary formation could be a dominant source of OA in this region (Mancilla et al., 2015). Our previous study has observed that SOC contributed 45 % to OC in Hong Kong during the summer of 2006 (Hu et al., 2010).
Here the analysis was expanded to samples taken during the 1-year period to obtain a more comprehensive understanding of sources and their contributions to ambient OA in Hong Kong and the factors that impact ambient SOA formation.

Characterization of SOA tracers and other polar oxygenated organic compounds
The concentration levels of 39 organic species, including 14 SOA tracers, 12 saccharides, 11 di-and tricarboxylic acids, 4-nitrocatechol, and cholesterol, under different meteorological conditions, are listed in Table 1.
3.1.1 SOA tracers of isoprene, monoterpenes, β-caryophyllene, and naphthalene Seven isoprene SOA (Isop_ SOA ) tracers, i.e., 2methylglyceric acid, two methyltetrol isomers (2methylthreitol and 2-methylerythritol), three C 5 -alkene triol isomers (cis-2-methyl-1,3,4-trihydroxy-1-butene, 3-methyl-2,3,4-trihydroxy-1-butene, and trans-2-methyl-1,3,4-trihydroxy-1-butene), and 3-MeTHF-3,4-diols (including both cisand trans-3-methyltetrahydrofuran-3,4-diols) were identified and quantified. The sum of all Isop_ SOA tracers ranged from 1.67 to 117.17 ng m −3 , with the annual mean value of 22.78 ± 26.06 ng m −3 . Among the Isop_ SOA tracers, methyltetrols and C 5 -alkene triols were the most abundant, and they are suggested to be formed through the acid-catalyzed ring-opening reactions of isoprene epoxydiols (IEPOXs) under low-NO x condition Surratt et al., 2010). Higher concentrations of Isop_ SOA tracers were measured in summer and autumn than in winter and spring. This could be caused by the higher temperature, stronger solar radiation, and higher emission of isoprene in summer and autumn than in the other two seasons, which promoted the SOA formation from isoprene. This seasonal pattern is consistent with what was observed in other studies (Ding et al., 2012;Kleindienst et al., 2007;Lewandowski et al., 2008). However, if we compare the levels of isoprene tracers monitored at different sites during summer, the total amount of Isop_ SOA tracers measured in Hong Kong was about 5 times lower than those measured in several cities in the USA and a rural site (Wangqingsha (WQS) site) in the PRD area (Ding et al., 2012;Kleindienst et al., 2007;Lewandowski et al., 2008). This may be due to the different levels of isoprene, OH radical, and NO x at these sampling sites. The 3-MeTHF-3,4-diols, tracers formed through the intermolecular rearrangement of IEPOX under acidic conditions, was identified in Hong Kong PM 2.5 samples for the first time. It has an annual mean concentration of 0.23 ± 0.10 ng m −3 , which was about 70 times lower than that in Birmingham, USA (Rattanavaraha et al., 2016), but was comparable to what was observed at the WQS site in the PRD area (Q. F. . 2-Methylglyceric acid, an isoprene tracer formed from methacrylic acid epoxide (MAE) and hydroxymethylmethyl-α-lactone (HMML) under high-NO x conditions (Lin et al., 2013;Nguyen et al., 2015), presented a quite different temporal trend from those of the other six Isop_ SOA tracers, with the highest concentration in winter and then autumn, summer, and spring. Chamber studies suggested that MAE is an oxidation product resulting from the OH addition to methacryloyl peroxynitrate (MPAN), and its production is temperature dependent (Roberts and Bertman, 1992;Worton et al., 2013). Under higher temperatures, the loss of MPAN is dominated by thermal decomposition, which does not produce SOA tracers through the NO/NO 2 pathway. Under lower temperatures, thermal decomposition of MPAN is limited, and more MPAN reacts with OH to generate MAE. Therefore, the lower temperatures in winter would favor the production of MAE and the MAE Isop_ SOA tracers, such as 2-methylglyceric acid. Moreover, all Isop_ SOA tracers exhibited higher concentrations on regional days than LRT and local days. On regional days, air masses transported from the PRD area worsened the air quality in Hong Kong, and the higher levels of gaseous pollutants, e.g., O 3 , NO 2 , O x , and SO 2 (Table 2), promoted SOA formation.
Generally speaking, at an urban location with anthropogenic NO x emissions from automobiles and power plants, the generation of Isop_ SOA tracers from the MAE NO/NO 2 pathway should be more favored than the IEPOX HO 2 channel. However, in this study, 94 % of the total mass of the quantified Isop_ SOA tracers was produced through the IEPOX HO 2 pathway. A similar phenomenon was observed at the WQS site in the PRD region (Q. F. . Therefore, to better understand the influences of environmental factors on isoprene SOA formation in the region, we applied the kinetic models described by Eddingsaas et al. (2010), Worton et al. (2013), and Birdsall et al. (2014) to investigate the fate of both IEPOX and MAE in the atmosphere. Besides their degradation through acid-catalyzed ring-opening reactions on particles, IEPOX and MAE can also be oxidized in the gas phase or removed by dry deposition . We applied the Kintecus kinetic model (Ianni, 2015) to quantitatively evaluate the fractions of these two Isop_SOA intermediates that undergo gasphase oxidation, aerosol-phase acid-catalyzed ring-opening reaction, and dry deposition processes. Details of the model calculations are provided in the appendices. Figure 1 shows the comparison of the three elimination processes of IEPOX and MAE during the sampling period in Hong Kong. Given the high volatility of MAE (vapor pressure: 9.2 × 10 −5 atm) (Worton et al., 2013), it has a low tendency to partition into the particle phase and its uptake onto aqueous particles is mainly governed by Henry's law constant (k value of MAE to be 7.5 × 10 6 M atm −1 , which is 20 times lower than that of IEPOX (1.3 × 10 8 M atm −1 ). Moreover, Riedel et al. (2015) suggested that the heterogeneous reactive uptake coefficient of MAE (γ = 4.9 × 10 −4 ) through the ring-opening reaction was a factor of 30 lower than that of IEPOX. Therefore, as shown in Fig. 1, MAE was primarily eliminated by dry deposition (> 80 %) in the gas phase, and only a trivial fraction was degraded through the ringopening reactions (≤ 2 %). Our results on the fate of MAE were similar to those observed at the University of California Blodgett Forest Research Station (UC-BFRS) (Worton et al., 2013). However, our results on the relative contributions of these three degradation pathways to IEPOX loss were quite different from theirs, indicating a more sensitive response of IEPOX than MAE to the change of environmental oxidants and conditions. Given the high liquid water content (LWC; mean: 57.20 ± 37.15 µg m −3 ) and particle acidity (H + p ; mean: −0.28 ± 0.42) of PM 2.5 samples in this study (Table 2), the particle-phase ring-opening reaction (F rop ) was the dominant degradation pathway of IEPOX in the Hong Kong atmosphere (average: 97.6 %), and its loss through dry deposition and gas-phase oxidation is almost negligible. The F rop of IEPOX reported by Worton et al. (2013) was only 0.02 %, mainly due to the much lower LWC (mean: 0.4 µg m −3 ) and weaker H + p (pH mean: 4.4) of their PM 2.5 samples. These results demonstrated that particle-phase LWC and H + p played a more significant role in the atmospheric degradation of IEPOX than MAE. Results from the kinetic model simulation were strongly supported by the experimental finding of IEPOX tracers as the dominant Isop_ SOA tracers measured in Hong Kong. The average ratio of IEPOX tracers to MAE tracers was 16.54 (ranged from 3.00 to 71.58), and the average value of F rop−IEPOX /F rop−MAE was 191.92, confirming that the IEPOX HO 2 channel is the major formation pathway of isoprene SOA in the region.
Five SOA tracers of monoterpenes (Mono_ SOA ), i.e., 3-hydroxyglutaric acid, 3-hydroxy-4,4-dimethylglutaric acid, 3-methyl-1,2,3-butanetricarboxylic acid, 3isopropylpentanedioic acid, and 3-acetyl pentanedioic acid, were identified and quantified. Their summed concentrations ranged from 2.54 to 32.57 ng m −3 , with an annual average value of 10.76 ± 8.04 ng m −3 , which is comparable to that reported at the WQS site in the PRD region but lower than that measured in the USA (Ding et al., 2012;Kleindienst et al., 2007;Lewandowski et al., 2008). All Mono_ SOA tracers showed the highest level on regional days (mean: 18.00 ± 9.28 ng m −3 ), followed by LRT (mean: 10.31 ± 7.33 ng m −3 ) and local days (mean: 6.41 ± 2.75 ng m −3 ) (Table 1). Although a higher emission and faster photochemical degradation of monoterpenes are expected in summer due to the intense solar radiation and high temperature, higher levels of Mono_ SOA tracers were monitored in autumn and winter than in the other two seasons, similar to what was observed at the WQS site (Ding et al., 2014). This seasonal trend of monoterpene SOA tracers may be partly due to the lower mixing height and temperature during autumn and winter, which favored the partition of Mono_ SOA tracers into the aerosol phase. Moreover, most of the regional days were identified in autumn and winter. The higher levels of NO x , O 3 , O x , and SO 2 on regional days (Table 2) are also responsible for the enhanced monoterpene SOA production in autumn and winter. Among Mono_ SOA tracers, 3-hydroxyglutaric acid (3HGA) was the most abundant, contributing ∼ 60 % to the total mass of Mono SOA tracers. Smog chamber experiments showed that the production yield of 3-methyl-1,2,3-butanetricarboxylic acid (MBTCA) from α-pinene/NO x oxidation was significantly higher than those from the β-pinene/NO x and d-limonene/NO x experiments (Jaoui et al., 2005). Therefore, the ratio of 3HGA / MBTCA was used as a criterion to differentiate SOA from α-pinene and other monoterpenes Y. Cheng et al.: Tracer-based source apportioning of atmospheric organic carbon (Ding et al., 2014). The value of this ratio was obviously higher on regional days (8.58 ± 2.69) than those on LRT (6.64 ± 3.63) and local days (5.62 ± 3.14), indicating that monoterpenes other than α-pinene, such as β-pinene and d-limonene, might have a more significant contribution to SOA on regional days in the region.
β-Caryophyllinic acid is the SOA (Cary_ SOA ) tracer of β-caryophyllene, and it ranged from 0.49 to 5.82 ng m −3 , with an average annual mean value of 1.53 ± 1.07 ng m −3 . Similar to the other SOA tracers, β-caryophyllinic acid showed the highest concentrations on regional days (mean: 2.33 ± 1.21 ng m −3 ) than LRT (1.73 ± 1.16 ng m −3 ) and local days (0.94 ± 0.41 ng m −3 ) ( Table 1). For its seasonal trend, β-caryophyllinic acid also exhibited the highest concentration in autumn and winter than the other two seasons. The SOA tracer of toluene, 2,3-dihydroxy-4-oxopentanoid acid, was undetectable in this study, mainly due to its trace level in the Hong Kong atmosphere (Hu et al., 2008) and the limited sensitivity of GC with quadrapole MS. Even in our previous study on a batch of summer PM 2.5 samples using a more sensitive GC with ion trap MS, it was barely quantified with a concentration of less than 1 ng m −3 in most samples (Hu et al., 2008). Phthalic acid was suggested as the SOA tracer of naphthalene, given its abundance in both naphthalene SOA and ambient OA . With the awareness of the potential uncertainties, e.g., the primary origin of phthalic acid from biomass burning, we adopted phthalic acid as the SOA tracer of naphthalene representing the SOA formation from anthropogenic VOCs. The concentration levels of phthalic acid ranged from 0.80 to 16.42 ng m −3 , with an average of 4.31 ± 3.39 ng m −3 . Similar to the other SOA tracers, it also showed the highest concentrations on regional days (7.16 ± 3.61 ng m −3 ) than LRT (4.97 ± 3.30 ng m −3 ) and local days (2.26 ± 1.38 ng m −3 ) ( Table 1).

Saccharides and dicarboxylic acids
Twelve saccharides, i.e., levoglucosan, arabitol, fructose, meso-erythritol, sucrose, galactosan, mannitol, sorbitol, galactose, glucose, xylose, and xylitol, have been quantified. Of the 12 saccharides, levoglucosan, a tracer of BB, was by far the most abundant (range: 0.64-474.15 ng m −3 ; mean: 75.02 ± 111.43 ng m −3 ). It showed the highest levels on regional days (about 6 times higher than that on local days), especially during winter when BB activities in the PRD region were most frequent. Two primary saccharides, i.e., fructose and xylose, also exhibited the highest levels on regional days. They showed good correlations with levoglucosan (R 2 = 0.65 and 0.93), suggesting that they could be from BB as well.
Among the identified dicarboxylic acids, oxalic acid was the most abundant, followed by terephthalic acid, phthalic acid, malic acid, succinic acid, and others. Most dicarboxylic acids, including the five most abundant ones, showed higher levels on regional days; they were found with higher levels in winter and autumn as well. This temporal trend is similar to what we have observed for Mono_ SOA tracers and most saccharides, indicating that regional pollution had a dominant influence on the abundance of both primary and secondary aerosols in Hong Kong, far exceeding the influence of other environmental parameters, such as temperature and solar radiation. Atmospheric dicarboxylic acids have various sources. For example, oxalic acid was suggested to be secondarily formed from biogenic emissions and anthropogenic sources (e.g., BB and automobile exhaust) through both gas-phase reactions and in-cloud processing (Yu et al., 2005). Malic acid was suggested to be the photodegradation product of both succinic acid and biogenic SOA compounds (Hu and Yu, 2013). In this study, malic acid was found to be strongly correlated with 3HGA (R 2 = 0.96) and Mono_ SOA tracers (R 2 = 0.95) throughout the year, providing more evidence to the hypothesis that malic acid is a late-stage oxidation product of BVOCs, especially monoterpenes (Hu and Yu, 2013). Ambient terephthalic acid was mainly directly emitted from plastic waste incineration (Simoneit et al., 2005) and was used as a marker of waste incineration.
Besides dicarboxylic acids, two benzenetricarboxylic acids (i.e., 1,2,3-and 1,2,4-benzenetricarboxylic acids), 4nitrocatechol, and cholesterol were also quantified. The two benzenetricarboxylic acids were suggested to be the photo-degradation products of polycyclic aromatic hydrocarbons (PAHs) emitted from combustion activities (Kautzman et al., 2010). We have previously identified them in the water-soluble humic-like substance (HULIS) extracts of PM 2.5 samples collected in Beijing and Hong Kong (Ma et al., 2018(Ma et al., , 2019. The annual mean concentrations of 1,2,3-and 1,2,4-benzenetricarboxylic acids measured in this study were 2.27 ± 1.97 ng m −3 (range: 0.47-9.50 ng m −3 ) and 3.13 ± 2.68 ng m −3 (0.47-12.54 ng m −3 ), respectively, which were comparable to what was measured at the other four sites in the PRD region (X. . The 4nitrocatechol, which was secondarily generated from the photo-oxidation of naphthalene, was suggested as a tracer of atmospheric aging of a BB plume (Kitanovski et al., 2012). It strongly correlated with levoglucosan (R 2 = 0.88) and exhibited higher levels on regional days and during winter, which further confirmed its BB origin in the region. Therefore, the two benzenetricarboxylic acids and 4-nitrocatechol were included in PMF analysis as the SOA tracers of BB aging.

Source apportionment of organic aerosols
In this study, the US EPA PMF 5.0 was used to determine the major OA sources and to quantify their contributions to OC. Eighteen species were input into PMF, including EC, OC, Ni, V, major ions, and various primary and secondary organic tracers. Given their similar origins, some organic trac-ers were lumped together, and the lumped species were used as the fitting species in PMF. They were (1) C 5 -alkene triols, which is the sum of the three C 5 -alkene triol isomers; (2) IsopT, which is the sum of two methyltetrol isomers and 2-methyl glyceric acid; (3) MonoT, which is the sum of the five monoterpenes SOA tracers; and (4) Hopane, which is the sum of the four hopanes. Since C 5 -alkene triols were not in the SOA tracers list of the TBM , the lumped C 5 -alkene triols were used as a separated fitting species in PMF. PMF solutions were tested with four to eight factors. A hundred base runs were performed in each modeling run, and the run with the minimum Q value was selected.
The uncertainty values of each input species were calculated using the method described in our previous studies (Hu et al., 2010;Ma et al., 2016), which were set to be 20 % of the mean concentrations for OC and EC and 40 % of mean values for cations, anions, and all organic species. An extra modeling uncertainty of 10 % was used to account for possible temporal changes in the source profiles. The Q Robust /Q True ratio was 1.00, and scaled residuals were normally distributed between −0.2 and 0.2, indicating no influence of outliers on the solution. A hundred bootstrap runs were performed with a minimum correlation R value to examine the base run solution's stability and uncertainty. All bootstrapped factors were explicitly mapped to factors identified in base solution with no exception. In the displacement (DISP) assessment, no error was found, and the drop in Q value was less than 1 %, suggesting a stable solution. No swap factor appeared at dQ max = 4, indicating there was no considerable rotational ambiguity in the solution. Rotations were introduced to the solutions by adjusting the "FPEAK" value from −1 to +1, and the non-rotated solutions (FPEAK = 0.0) were considered to be the most interpretable ones. Moreover, a strong linear correlation between the measured and PMF-predicted OC (OC PMF ) (R 2 = 0.92) was observed, which also suggested a reliable PMF solution.
As shown in Fig. 2, the first factor was distinguished by high loadings of oxalate and biogenic SOA tracers, suggesting the secondary origin of this source. The second factor was dominated by large amounts of SO 2− 4 and NH + 4 , suggesting the process of secondary sulfate formation. In the third factor, about 90 % of levoglucosan was portioned into it, accompanied by 4-nitrocatechol, phthalic acid, and the two benzenetricarboxylic acids, indicating both the primary emission and aging of BB plume. Therefore, this factor was defined as BB and SOA (BB/SOA). The fourth factor was identified as vehicular emissions due to the large amounts of hopanes and EC portioned in this factor. The fifth factor has large amounts of Ni and V, which are signatures of residual oil combustion from the marine vessel class (Viana et al., 2009). It is well known that Hong Kong is one of the busiest container ports globally, which handles 50 % of the PRD's total cargo throughput. Therefore, the fifth factor was identified as marine vessels. The sixth factor has a high loading of Na + , Mg 2+ , and Ca 2+ , indicating the sea salt source.
The amount of OC apportioned to each factor in PMF analysis was considered the contribution of that source to ambient OC. Therefore, the two leading sources contributing to ambient OC in Hong Kong were BB (including both primary emission and aging process, OC BB : 27.9 %, 1.17 ± 1.99 µg C m −3 ) and SOA (SOC SOA : 27.5 %, 1.15 ± 0.82 µg C m −3 ), followed by marine vessels (OC marine : 15.6 %, 0.65 ± 0.58 µg C m −3 ), secondary sulfate (SS) (SOC SS : 14.5 %, 0.60 ± 0.46 µg C m −3 ), vehicle emissions (OC vehicle : 10.5 %, 0.44 ± 0.42 µg C m −3 ), and sea salt (OC sea : 4.0 %, 0.17 ± 0.19 µg C m −3 ) ( Table 2 and Fig. 3). Since a fraction of SOA from the aging of BB (SOC BB ) was apportioned into the BB/SOA factor, we calculated SOC BB using the following equation: where OC BB and [LEVO BB ] are the amounts of OC and levoglucosan portioned in the BB/SOA factor, respectively. Using levoglucosan as the tracer of primarily emitted BB OA, we calculated the amounts of primary organic carbon (POC) from BB (POC BB ) by dividing [LEVO] BB with 0.082, where 0.082 is the average ratio of levoglucosan to POC from the burning of major types of Chinese cereal straws (i.e., rice, wheat, and corn) obtained in combustion chamber experiments (Zhang et al., 2007). As cereal straws are one of the most common BB fuels in China, the above ratio (0.082) has been used to estimate BB contribution to POC in both Beijing (Zhang et al., 2008) and Hong Kong (Sang et al., 2011). Therefore, it was adopted to calculate POC BB in this study. Based on PMF results, the source-specific contributions to OC are presented in Table 2 and demonstrated in Fig. 3. The total SOC apportioned by PMF (SOC PMF ), i.e., the sum of SOC SOA , SOC SS , and SOC BB , accounted for 51.4 % (2.15 ± 1.37 µg C m −3 ) of OC in Hong Kong, with the secondary organic-rich sources (i.e., SOC SOA + SOC BB ) contributing 36.9 % (1.54 ± 1.13 µg C m −3 ) to total OC. Huang et al. (2014) also reported that secondary organic-rich sources accounted for 30 %-40 % of OC in Guangzhou, which is another PRD site. A higher level of SOC PMF and its contribution to OC were observed on regional days (3.27 ± 1.18 µg C m −3 , 57.4 %) compared to LRT (2.36 ± 1.54 µg C m −3 , 53.0 %) and local days (1.36 ± 0.81 µg C m −3 , 43.6 %). An even starker difference in the amounts of SOC BB between regional and local days was observed, which was 8 times higher on the regional days. This suggested that nonlocal sources were the dominant contributors to SOC BB . BB activities were intensive in the PRD region, especially during fall and winter. On regional days, freshly emitted and aged gaseous and aerosol-phase pollutants from the open burning of rice straws and other crops were transported from the northern PRD region into Hong Kong (Hu et al., 2010). Huang et al. (2014) examined the aging of BB plume at low temperatures. They found that the production of BB SOA was rapid at a typical OH Y. Cheng et al.: Tracer-based source apportioning of atmospheric organic carbon  radical concentration of wintertime China, and the amount of BB SOA may exceed BB POA in 4-14 h even at −10 • C. Given that the average temperatures in Hong Kong during autumn and winter were 26.15 and 17.76 • C, respectively; the formation of BB SOA should be achieved even quicker during the regional transport. As expected, SOC SOA also showed a higher average concentration on regional days (1.75 ± 0.75 µg C m −3 ) than on LRT (1.14 ± 0.82 µg C m −3 ) and local days (0.78 ± 0.65 µg C m −3 ), which is consistent with the trends of all SOA tracers. Although SOC from secondary inorganic-rich sources (SOC SS ) exhibited the highest levels (0.82 ± 0.38 µg C m −3 ) on regional days as well, its contribution to OC was relatively stable under the three synoptic conditions (Fig. 3). Several studies showed that SO 2 transported from the northern PRD region promoted secondary sulfate formation in Hong Kong through both gasphase and in-cloud oxidation pathways (Lu and Fung, 2016;Yu et al., 2005;Yuan et al., 2006). A recent study proposed that the sulfate formation in aqueous aerosols through NO 2 oxidation and ammonium neutralization can simultaneously enhance the production of both nitrate and SOA , which helps explain the considerable amount of SOC SS apportioned.
OC from the four primary sources, i.e., POC BB , OC marine , OC vehicle , and OC sea , accounted for 48.6 % of total OC throughout the year. Similar to SOC BB , POC BB showed a higher level (1.38 ± 1.75 µg C m −3 ) on regional days due to a large number of emissions from BB activities in the northern PRD area. OC vehicle remained a higher contribution on local days (15.6 %, 0.49 ± 0.46 µg C m −3 ), consistent with our previous finding that vehicle emission is a local pollution source (Hu et al., 2010). Similarly, marine vessels accounted for a greater amount and larger fraction of OC on local days (32.0 %, 1.00 ± 0.63 µg C m −3 ) than LRT (5.2 %, 0.23 ± 0.19 µg C m −3 ) and regional days (6.5 %, 0.37 ± 0.21 µg C m −3 ). On local days, the southeastern to southwestern wind brought pollutants from residual oil com-bustion from the ocean into Hong Kong, leading to a higher OC marine .
In summary, both secondary aerosol sources and air mass origins play important roles in atmospheric OC in Hong Kong. On regional days, air mass transported from the northern PRD area brought large amounts of air pollutants into Hong Kong, which promoted the SOA production from both anthropogenic emissions and BVOCs and resulted in a fraction of 57.4 % of OC being secondarily formed. On the other hand, local sources, including vehicle emissions and marine vessels, became more critical and significantly contributed to OC (56.4 %) on local days.

Estimation of SOC origin
To better understand the SOA precursors and their contributions to SOA/SOC in the region, we adopted a tracerbased method  to estimate the SOA/SOC formation from a group of selected biogenic and anthropogenic hydrocarbons, i.e., isoprene, monoterpenes, β-caryophyllene, and naphthalene. The mass ratio of tracer compounds to the total SOC (f SOC ) generated from individual VOC precursors was derived from smog chamber experiments Offenberg et al., 2007). By assuming the same f SOC value of the precursor under smog chamber conditions and in ambient air, one can use the quantified SOA tracer concentrations to estimate the amount of SOC from that precursor in the real atmosphere. It has been well noted that results obtained from this tracer-based method are subject to potential uncertainties from various aspects, e.g., the larger variation of precursor concentrations and more complicated environmental conditions in the real atmosphere than in smog chamber experiments; the decay of some tracer compounds during transport; mismatch of ambient and smog-chambergenerated SOA compositions, using surrogates other than ketopinic acid for the quantification of tracer compounds; and so on (Ding et al., 2014;Hu et al., 2008;Kleindienst et al., 2012Kleindienst et al., , 2007. However, using the tracer-based method, we can at least have a rough estimation of the key SOA precursors in the region, their contributions to ambient OC, and the amount of SOC from unknown precursors. Wang et al. (2013) noted that the SOA tracer-based method would significantly underestimate SOC_ Mono in the PRD region. Ding et al. (2014) gave a reasonable explanation that the mismatch of monoterpene tracers measured in ambient air and used to derive f SOC of monoterpenes in chamber studies may increase the uncertainty of SOC_ Mono . Thus, they picked the five Mono_ SOA tracers measured in their samples and derived the f SOC and f SOA values using the SOA tracers data and SOA/SOC concentrations reported by Offenberg et al. (2007). In this study, we only measured five out of nine monoterpene SOA tracers in Offenberg et al. (2007). Similar to Ding et al. (2014), to lower the uncertainty induced from the mismatch of SOA tracer compositions, we derived a f SOC_mono value of 0.047 based on the experimental data by Offenberg et al. (2007) and applied it to estimate SOC_ Mono . Many research groups have adopted this tracerbased method to assess SOC productions from the five studied VOCs at various locations in the world, and reasonable results have been obtained (Ding et al., 2012;Fu et al., 2014;Hu et al., 2008;Hu and Yu, 2013;Kleindienst et al., 2012Kleindienst et al., , 2007Lewandowski et al., 2008).
As shown in Table 2, SOC estimated by the tracer-based method (SOC TBM ) ranged from 0.11 to 1.53 µg C m −3 in Hong Kong, accounting for 3.8 % to 22.7 % of ambient OC levels. It exhibited the same trend as OC and SOC PMF (Fig. 4), i.e., with higher concentrations on regional days (0.81 ± 0.35 µg C m −3 ) than on LRT (0.50 ± 0.29 µg C m −3 ) and local days (0.28 ± 0.13 µg C m −3 ). Similar to our previous study, monoterpenes were found to be the most significant SOC contributor in the region, with SOC_ Mono ranging from 0.05 to 0.69 µg C m −3 and having an average concentration of 0.23 ± 0.17 µg C m −3 . SOC_ Iso and SOC_ Cary , on the other hand, were about 3 times smaller than SOC_ Mono and were 0.08 ± 0.09 and 0.07 ± 0.05 µg C m −3 , respectively.
Smog chamber experiments have been carried out to study the SOA yields from q OH oxidation, ozonolysis, and nitrate radical (NO 3 ) oxidation of monoterpenes and isoprene, and monoterpenes were found to be more effective in SOA production than isoprene (Lee et al., 2006a, b). Highly oxygenated organic molecules with low and extremely low volatilities were formed from the oxidation of monoterpenes and observed in both laboratory experiments and field measurements (Ehn et al., 2014;Jokinen et al., 2015;. Moreover, a synergistic O 3 + OH oxidation pathway of monoterpenes was recently proposed, which leads to the formation of extremely low-volatility oligomers and may result in even larger monoterpene SOA yields in the real atmosphere than what observed in the smog chamber experiments (Kenseth et al., 2018). Tsui et al. (2009) reported a total BVOC emission of 8.6 × 10 9 g C yr −1 in Hong Kong, with 40 % from monoterpenes and 30 % from isoprene. The remaining 30 % could be sesquiterpenes (e.g., β-caryophyllene) or other BVOCs. Therefore, the predominance of monoterpenes SOA in BVOCs-derived SOC is likely due to the combined effects of their high SOA yields and large emissions in the region. Like the SOA tracers, SOC from the four precursors all showed the highest level on regional days than those on LRT and local days (Table 2). On regional days, large amounts of VOC precursors and gaseous oxidants could be brought into Hong Kong through the regional transport of air masses from northern PRD and oxidized along the way. Conversely, on local days, the ocean breeze brings clean air masses from the South China Sea into Hong Kong, leading to a dilution effect of local air pollution. These results highlight that air mass origins play an important role in the SOC formation from both biogenic and anthropogenic VOCs. Given that SOC TBM is calculated based on the concentration levels of individual SOA tracers mea- sured in the ambient aerosols, it is reasonable that SOC attributed to each VOC precursor showed the same meteorological variations as their SOA tracers.
We observed similar temporal trends between SOC PMF and SOC TBM (R 2 = 0.71). However, SOC TBM only accounted for 26.5 % of SOC PMF , suggesting SOC must have been underestimated by the tracer-based method. A reasonable explanation is that secondary formation from nighttime reactions, multiphase reactions, and other SOA precursors are not considered in the SOA tracer-based method. Parameters used in the tracer-based method were derived from pure gas-phase photo-oxidation of VOC precursors in smog chambers (Kleindienst et al., , 2009. Therefore, it is better to be used as a complementary method with PMF in the source apportionment study of ambient OC, especially SOC.

Effects of anthropogenic influences on secondary aerosol formation
Increasing evidence from laboratory studies and ambient observations has shown that anthropogenic emissions can significantly affect SOA formation from terpenoids through multiple chemical processes in both daytime and nighttime (Xu et al., 2015;Zhang et al., 2018). We conducted the Pearson's R correlation analysis of all SOC terms (i.e., SOC_ Iso , SOC_ Mono , SOC_ Cary , SOC_ Nap , SOC TBM , SOC PMF , SOC BB , SOC SOA , and SOC SS ) with O 3 , NO 2 , SO 2 , O x , NO 3 , sulfate, particle acidity (H + p ), and particle liquid water content (LWC P ) ( Table 3). Details on the calculation of H + p and LWC P are presented in Appendix B. Since NO 3 was not directly monitored at HKEPD stations, its mixing ratio was estimated using the following equation: The numerator is the production of NO 3 (p[NO 3 ]) from O 3 and NO 2 , and the denominator is the reactivity of NO 3 for NO 3 -VOC reactions. From the IUPAC database, we obtained the temperature-dependent expression of k 1 (cm 3 molecules −1 s −1 , the production rate constant of NO 3 ) as 1.4 × 10 −13 exp(−2470/T ), where T is the ambient temperature in kelvin. Therefore, using k 1 and the measured concentration levels of [O 3 ] and [NO 2 ], we calculated p[NO 3 ] ( Table 2). Brown et al. (2016) reported a NO 3 -VOC reactivity of 6.5 ± 6.8 × 10 −3 s −1 in Hong Kong with a corresponding NO 3 lifetime of 2.5 min. NO 3 was then calculated as the ratio of p[NO 3 ] to this NO 3 reactivity value, and an annual mean level of 70 ± 47 ppt was estimated. As we mentioned earlier, O x is an indicator of atmospheric oxidation capacity. Five SOC terms, i.e., SOC_ Mono , SOC_ Nap , SOC TBM , SOC SOA , and SOC PMF , showed significant positive correlations with O x , especially SOC SOA and SOC PMF (R > 0.7, P < 0.01). However, only SOC SOA and SOC SS were found to be significantly correlated with O 3 (R > 0.50, P < 0.01). As for NO 2 , another critical component of O x , it exhibited statistically significant positive correlations with not only SOC SOA and SOC PMF but also several TBM estimated SOCs, including SOC TBM , SOC_ Mono , SOC_ Nap , and SOC_ Cary . This may be because SOA tracers used in TBM were produced from the photo-oxidation of these VOC precursors in the presence of NO x . The significant positive correlations between NO 2 and SOC SOA and SOC PMF also suggest that the daytime oxidation processes involving NO x are critical SOA formation pathways in the region. Significant correlations with R > 0.5 between NO 3 and SOC_ Mono , SOC SOA , SOC PMF , and SOC SS were also observed. BVOCs were found to account for > 80 % of the NO 3 reactivity in Hong Kong (Brown et al., (2018) and Xu et al. (2013) have reported an enhancement of nighttime monoterpenes SOA in the southeastern USA by NO 3 -monoterpene reactions. Therefore, our findings indicate that SOA formation through nighttime NO 3 oxidation of biogenic VOCs, especially monoterpenes, may have made a considerable contribution to the SOA loading in Hong Kong. However, more field measurement data, e.g., quantification of the particle-phase organic nitrates using real-time online mass spectrometry techniques, are needed to examine the impact of NO x processing on SOA formation in the region. Since NO 3 is a key precursor of nighttime production of HNO 3 , and nitrate is a significant component of secondary inorganic aerosols, it rationalized the correlations between NO 3 and SOC SS . Six SOC terms, i.e., SOC_ Mono , SOC_ Nap , SOC TBM , SOC SOA , SOC SS , and SOC PMF , showed significant positive correlations with sulfate, especially SOC SS and SOC PMF (R ≥ 0.8, P < 0.01). Given that sulfate is the key component of secondary inorganic aerosol, such a strong correlation between SOC SS and sulfate is expected. Moreover, several studies have suggested that sulfate also plays a dominant role in the production of aerosol-phase organosulfates through both nucleophilic addition reactions and the saltingin effect (Lin et al., 2012;Riva et al., 2015;Xu et al., 2015).
We then performed multivariate linear regression (MLR) analysis to obtain a quantitative and comprehensive understanding of the impacts of gaseous oxidants and aerosol characteristics on SOC TBM , SOC PMF , and the individual PMFapportioned SOCs (i.e., SOC SS , SOC SOA , and SOC_ BB ). Six parameters, namely, O 3 , NO 2 , NO 3 , sulfate, H + p , and LWC P , were included in the preliminary runs. However, the MLR results showed that O 3 was an insignificant factor for all SOC terms, even with negative regression coefficients. Pearson's R analysis also showed that SOC was more NO 2 dependent than O 3 . Therefore, it was excluded from the final MLR analysis, and the results are shown in Table 4.
We found that two parameters, i.e., sulfate and NO 3 , have statistically significant positive linear relationships (P ≤ 0.001) with SOC SS , and the regression coefficients were 0.913 and 0.234, respectively. The result is reasonable and consistent with what was observed from Pearson's R analysis, given that sulfate is the critical component in the PMF identified SS factor, and NO 3 is the precursor of nitrate through HNO 3 formation at nighttime. As for SOC_ BB , three parameters, i.e., NO 2 , NO 3 , and H + p , showed significant linear relationships with it (P < 0.01), with a regression coefficient of 0.639, −0.509, and 0.503, respectively. This indicates that a 1 mol L −1 increase in particle acidity was associated with a 0.503 µg C m −3 increase in SOC from BB aging. Phenols, which are produced from the combustion of lignin, are a typical class of gaseous compounds emitted in large amounts from BB (Bruns et al., 2016;Schauer et al., 2001). Recent laboratory studies indicate that phenols can undergo multiphase photochemical reactions in the atmosphere with the formation of nitrophenols and nitrocatechols (Finewax et al., 2018;Yu et al., 2014). Vione et al. (2001) observed the aqueous-phase photonitration of phenols, which was pH dependent with more nitro compounds generated at lower pH. Given the strong particle acidity (pH annual mean: −0.28) observed in the Hong Kong atmosphere, the formation of the 4-nitrocatechol and its analogs may be favored in the BB plume, which enhances SOC BB formation.
Both sulfate and NO 2 were found as the statistically significant factors that positively correlated with SOC PMF , with regression coefficients of 0.530 and 0.373, respectively (P < 0.001, Table 4). This suggests reducing the sulfate level by 1 µg m −3 and NO 2 level by 1 ppb could lower the total PMF-apportioned SOC by 0.530 and 0.373 µg C m −3 , respectively. NO 2 was also the most significant factor influencing SOC TBM , with a regression coefficient of 0.383 (P < 0.001), indicating that a decrease of NO 2 by 1 ppb can reduce SOC TBM by 0.383 µg C m −3 . As for SOC SOA , we found NO 3 as the most significant parameter (P < 0.01), and a decrease of 1 ppb NO 3 can lead to a reduction of SOC SOA by 0.384 µg C m −3 when holding other covariates unchanged. These results are consistent with what was observed from Pearson's R analysis, indicating the importance of NO x pro- Table 4. Results of multivariate linear analysis of PMF and TBM-resolved SOCs, NO 2 , NO 3 , sulfate, particle acidity (H + p ), and total particlephase liquate water content (LWC P ). **: P < 0.01; *: P < 0.05. Note that significant regressions are bold. cessing on both daytime and nighttime SOA production in the region.

Conclusions
In this study, we identified and quantitatively assessed the contributions of six primary and secondary sources to ambient OC in Hong Kong, and secondary formation was found to be the leading contributor. Anthropogenic emissions, including NO 2 , O x , NO 3 , and sulfate, significantly influenced SOA formation in the region. In particular, NO x processing in both daytime and nighttime has played a critical role. Although the ambient NO 2 level has dropped by 33.3 % from 1999 to 2019 (the government of Hong Kong Special Administrative Region (HKSAR), https://www.info.gov.hk/gia/general/ 202001/20/P2020012000874.htm, last access: 16 June 2021) and sulfate level in PM 2.5 was also lowered by about 30 % from 2000 to 2016 (Environmental Protection Department of Hong Kong, 2017), the roadside NO 2 level was still high. According to the 20-year air pollutants monitoring data released by HKSAR, the annual average concentrations of roadside NO 2 were much higher than the other gaseous pollutants, and they peaked during 2011-2013, which were 122 and 118 µg m −3 in 2011 and 2012, respectively. Although the annual ambient level of roadside NO 2 decreased to 80 µg m −3 in 2019, it is still 2 times higher than the annual objective level set by the HKSAR government, indicating a continuous significant impact of NO x on SOA formation in Hong Kong, especially in areas with heavy traffic load. Given that 90 % of the roadside NO 2 was from commercial vehicles, such as buses, trucks, minibuses, and so on, our results suggest that more stringent control of NO x emission from commercial vehicles is needed. This will benefit the community by reducing not only the background NO x levels but also the SOA pollution in Hong Kong.
Data availability. Raw data used in this study are archived at Hong Kong Baptist University and are available upon request by contacting the corresponding author.
Author contributions. YC and DH designed the study. YC did all the experiments and most of the data analysis. YM helped with regression analysis and data interpretation. YC drafted the article. DH helped with data analysis and interpretation and revised the article.
Competing interests. The authors declare that they have no conflict of interest. Review statement. This paper was edited by Nga Lee Ng and reviewed by three anonymous referees.