Formation of condensable organic vapors from anthropogenic and biogenic volatile organic compounds (VOCs) is strongly perturbed by NO x in eastern China

. Oxygenated organic molecules (OOMs) are the crucial intermediates linking volatile organic compounds (VOCs) to secondary organic aerosols (SOAs) in the atmosphere, but comprehensive understanding of the characteristics of OOMs and their formation from VOCs is still missing. Ambient observations of OOMs using recently developed mass spectrometry techniques are still limited, especially in polluted urban atmospheres where VOCs and oxidants are extremely variable and complex. Here, we investigate OOMs, measured by a nitrate-ion-based chemical ionization mass spectrometer at Nanjing in eastern China, through performing positive matrix factorization on binned mass spectra (binPMF). The binPMF analysis reveals three factors about anthropogenic VOC (AVOC) daytime chemistry, three isoprene-related factors, three factors about biogenic VOC (BVOC) nighttime chemistry, and three factors about nitrated phenols. All factors are inﬂuenced by NO x in different ways and to different extents. Over 1000 non-nitro molecules have been identiﬁed and then reconstructed from the selected solution of binPMF, and about 72 % of the total signals are contributed by nitrogen-containing OOMs, mostly regarded as organic nitrates formed through peroxy radicals terminated by nitric oxide or nitrate-radical-initiated oxidations. Moreover, multi-nitrates account for about 24 % of the total signals, indicating the signiﬁcant presence of multiple generations, especially for isoprene (e.g., C 5 H 10 O 8 N 2 and C 5 H 9 O 10 N 3 ). Additionally, the distribution of OOM concentration on the carbon number conﬁrms their precursors are driven by AVOCs mixed with enhanced BVOCs during summer. Our results highlight the decisive role of NO x in OOM formation in densely populated areas, and we encourage more studies on the dramatic interactions between anthropogenic and biogenic emissions.

Abstract. Oxygenated organic molecules (OOMs) are the crucial intermediates linking volatile organic compounds (VOCs) to secondary organic aerosols (SOAs) in the atmosphere, but comprehensive understanding of the characteristics of OOMs and their formation from VOCs is still missing. Ambient observations of OOMs using recently developed mass spectrometry techniques are still limited, especially in polluted urban atmospheres where VOCs and oxidants are extremely variable and complex. Here, we investigate OOMs, measured by a nitrate-ion-based chemical ionization mass spectrometer at Nanjing in eastern China, through performing positive matrix factorization on binned mass spectra (binPMF). The binPMF analysis reveals three factors about anthropogenic VOC (AVOC) daytime chemistry, three isoprene-related factors, three factors about biogenic VOC (BVOC) nighttime chemistry, and three factors about nitrated phenols. All factors are influenced by NO x in different ways and to different extents. Over 1000 nonnitro molecules have been identified and then reconstructed from the selected solution of binPMF, and about 72 % of the total signals are contributed by nitrogen-containing OOMs, mostly regarded as organic nitrates formed through peroxy radicals terminated by nitric oxide or nitrate-radical-initiated oxidations. Moreover, multi-nitrates account for about 24 % of the total signals, indicating the significant presence of multiple generations, especially for isoprene (e.g., C 5 H 10 O 8 N 2 and C 5 H 9 O 10 N 3 ). Additionally, the distribution of OOM concentration on the carbon number confirms their precursors are driven by AVOCs mixed with enhanced BVOCs during summer. Our results highlight the decisive role of NO x in OOM formation in densely populated areas, and we encourage more studies on the dramatic interactions between anthropogenic and biogenic emissions.

Introduction
Secondary organic aerosols (SOAs), as an important and complex component of submicron particles (Zhang et al., 2007;Jimenez et al., 2009;Huang et al., 2014), are fully involved in affecting climate (IPCC, 2013) and causing health risks (Nel, 2005;Lim et al., 2012). Volatile organic compounds (VOCs) are ubiquitous in the atmosphere and are recognized as the main precursors of SOAs (Hallquist et al., 2009;Ziemann and Atkinson, 2012). However, the missing intermediate processes from VOCs to SOAs are yet to be elucidated (Hallquist et al., 2009;Ehn et al., 2014).
Benefiting from state-of-the-art measurement technics (Bertram et al., 2011;Jokinen et al., 2012;Lee et al., 2014), many previously unreported oxygenated organic molecules (OOMs), as intermediates from VOCs to SOAs (Ziemann and Atkinson, 2012), have been discovered. Among OOMs, highly oxygenated organic molecules (HOMs), first observed in the gas phase at a boreal forest site  and reviewed by Bianchi et al. (2019), are so functionalized and have such low volatility that they can participate at the beginning of new particle formation (NPF) by stabilizing sulfuric acid Riccobono et al., 2014) or through clustering alone Bianchi et al., 2016), and they condense on existing particles and are thus responsible for a large fraction of SOAs . In addition to conventional VOC-to-OOM mechanisms summarized in the Master Chemical Mechanism (MCM) (http://mcm.york.ac.uk/, last access: 9 February 2021), recent studies have proposed new pathways, such as autoxidation (Crounse et al., 2013;Jokinen et al., 2014) and multigenerational oxidation (Rollins et al., 2012;Wang et al., 2020b), to form condensable vapors by adding oxygen atoms efficiently. The production of OOMs, especially HOMs, from precursors such as monoterpenes Jokinen et al., 2015;Berndt et al., 2016), sesquiterpenes , isoprene (Jokinen et al., 2015;Zhao et al., 2021), aromatics Molteni et al., 2018;Garmash et al., 2020), and alkanes  has been investigated in laboratories by using chemical ionization atmospheric pressure interface time-of-flight mass spectrometers with nitrate reagent ions (nitrate CI-APi-TOFs).
New insights and a general understanding about OOMs have been attained, yet many critical details about OOM formation and properties need to be addressed. First, the current kinetic descriptions of OOMs obtained from experiments are still limited, lacking for example individual H-shift rates for autoxidation and reaction rates of multigenerational products with oxidants. Furthermore, the complexity of the real atmosphere makes it more difficult to apply experimental results to ambient environments. The precursors compete for oxidants and vice versa, and their products will interact mechanistically in mixtures of atmospheric vapors (McFiggans et al., 2019;Heinritzi et al., 2020). However, in the labora-tory we usually study simple systems with a single precursor and a single oxidant. Moreover, most experiments are carried out for environments dominated by biogenic VOCs (BVOCs), while anthropogenic emissions receive less attention. In addition to classic anthropogenic VOCs (AVOCs), large amounts of primary emissions of oxygenated VOCs are also present in urban areas (Karl et al., 2018). The effect of NO x on OOMs is another key issue. NO x can terminate peroxy radicals (RO 2 ), outcompeting autoxidation propagation reactions and other bimolecular reactions (RO 2 + RO 2 , RO 2 + HO 2 ), and change the products' distribution and, consequently, size-dependently modulate the growth rates of organic aerosol particles . Additionally, NO x contributes non-linearly to atmospheric oxidants, which also influence the production of OOMs (Pye et al., 2019). It is anticipated that NO x plays a varied role in the formation of OOMs as well as of SOAs in different environments.
Therefore, more extensive OOMs observations are needed to validate the atmospheric implications of experiments, to be coupled with the global or regional model, and finally to comprehensively understand the fate of OOMs in the atmosphere. Until now, only a few ambient observations of OOMs using a nitrate CI-APi-TOF have been reported , and almost all of them focus on rural, forested, or remote atmospheres Massoli et al., 2018;Zhang et al., 2020;Beck et al., 2021).
The Yangtze River Delta (YRD) is one of the most developed regions in eastern China. Fine particulate matter, with an aerodynamic diameter smaller than 2.5 µm (PM 2.5 ), has been significantly reduced in eastern China since the implementation of the "Air Pollution Prevention and Control Action Plan" in 2013 , but (secondary) organic aerosols are still much more abundant than in clean areas Sun et al., 2020). Here we investigate condensable oxygenated organic vapors observed by a nitrate CI-APi-TOF in August-September 2019 at the Station for Observing Regional Processes of the Earth System (SORPES) in the western part of the YRD, an anthropogenicemissions-dominated environment Xu et al., 2017) mixed with enhanced biogenic emissions during summer (Wang et al., 2020a;Xu et al., 2021). A variety of oxidants Li et al., 2020; with numerous precursors (VOCs) suggest very complicated atmospheric oxidation processes and thousands of products (OOMs). Thereby, positive matrix factorization (PMF) (Paatero and Tapper, 1994) was applied to timeresolved mass spectra which had been pre-divided into small bins (binPMF; Zhang et al., 2019) to separate various sources or processes of OOMs. Combined with summarizing the ensemble chemical characteristics of OOMs, some interesting findings about the conversion of VOCs to OOMs were obtained.

Instrumentation
The nitrate CI-APi-TOF (Aerodyne Research Inc. and Tofwerk AG), combining a chemical ionization source (CI) and an atmospheric pressure interface time-of-flight mass spectrometer (APi-TOF) equipped with a long time-of-flight model (LTOF) with a mass resolution of 8000-12 000 Th Th −1 (Th denotes thomsons), was deployed to detect ambient sulfuric acid and OOMs. The ambient air was pulled into a laminar flow reactor, where the sample flow (10 L min −1 ) was surrounded by a purified airflow serving as the sheath flow (25 L min −1 ), through a stainless-steel tube (100 cm long, 3/4 in. diameter). Nitrate reagent ions were generated in the sheath flow by exposing air-containing nitric acid to a photoionizer X-ray (Model L9491, Hamamatsu, Japan). Detailed description of the instrument has been given elsewhere Jokinen et al., 2012). The data were acquired at a 1 Hz time resolution and analyzed with a tofTools package (version 6.11) based on MATLAB (Math-Works Inc.). Due to the diversity and unknown molecular structures of oxygenated organic compounds, standards for OOMs measurable by the nitrate CI-APi-TOF are still lacking. Like in other studies Trostl et al., 2016;Stolzenburg et al., 2018), an empirical method was used to quantify the concentrations of OOMs based on the ionization kinetics (pseudo-first-order reaction approximation) in the reaction tube of CI (Eq. 1) .
Here [OOM i ] is the concentration (molecules cm −3 ) of one OOM. On the right side of the equation, the numerator in the parentheses is the detected total signals (ions s −1 ) of one OOM charged by nitrate ions in adduct-forming or deprotonated ways, and the denominator is the sum of all reagent ion signals (ions s −1 ). First, a H 2 SO 4 -based calibration factor C, with a value of 4.2 × 10 9 molecules cm −3 , was obtained from a calibration using H 2 SO 4 (Kuerten et al., 2012) that proceeded taking into account the diffusion loss in the sampling line by assuming that all detected OOMs have the same ionization efficiency as H 2 SO 4 . The collision frequency of HOMs with nitrate clusters is comparable to that of sulfuric acid with nitrate clusters Hyttinen et al., 2015), yet the collision frequency of some moderately oxygenated molecules with nitrate clusters is relatively slower. Therefore, calibration by this method leads to a lower-limit estimate of OOM concentrations Trostl et al., 2016), but the accurate quantification of OOMs is not the main concern of this study, and the errors in the quantification of OOMs do not change our conclusions. Second, a mass-dependent transmission efficiency T i of the APi-TOF was inferred in a separate experiment by depleting the reagent ions with several perfluorinated acids . VOC precursors were measured by a proton transfer reaction time-of-flight mass spectrometer (PTR-TOF-MS; Ionicon Analytik, Innsbruck, Austria, TOF 1000 ultra). PM 2.5 was measured with a combined technique of light scattering photometry and beta radiation attenuation (Thermo Scientific SHARP Monitor Model 5030). The chemical compositions of PM 2.5 were determined on-line using a timeof-flight aerosol chemical speciation monitor (TOF-ACSM, Aerodyne Research Inc.). PMF analysis was further used to separate the organic aerosols (OAs) to primary organic aerosols and secondary organic aerosols (POAs and SOAs). The number concentrations of particles were measured by the scanning mobility particle sizer (SMPS) with a nano differential mobility analyzer (DMA) (4.0 to 63.8 nm) and long DMA (41.4 to 495.8 nm) and the aerodynamic particle sizer (APS) (0.5 to 18.0 µm). NO and NO 2 were measured using a chemiluminescence analyzer equipped with a blue-light converter (TEI, Model 42i-TL); O 3 , SO 2 , and CO were measured using the ultraviolet photometry, pulsed-UV fluorescence, and IR (infrared) photometry techniques (TEI, Model 49i, 43C, and 48C) respectively. Zero and span calibrations for trace gases were performed weekly during the campaign. Meteorological measurements including relative humidity (RH), wind speed, wind direction, and air temperature were recorded by an automatic weather station (Campbell Scientific, AG1000). J (O 1 D) was measured by an ultrafast charge-coupled device detector spectrometer, which is UVB enhanced (Meteorologie Consult GmbH, Germany).

Hydroxyl radical (OH) estimate
The OH concentration was calculated by applying Eq. (2), based on the assumption that gaseous sulfuric acid is mostly produced from the oxidation of SO 2 by OH and primarily lost by condensing onto particles: where k OH+SO 2 is a termolecular reaction constant for the rate-limiting step of the formation pathway of H 2 SO 4 in the atmosphere (Finlayson-Pitts and Pitts, 2000) and the condensation sink (CS) is the loss rate of H 2 SO 4 by condensation to the aerosol surface. The value of k OH+SO 2 is inferred from the IUPAC Task Group on Atmospheric Chemical Kinetic Data Evaluation (https://iupac-aeris.ipsl.fr/, last access: 9 August 2021). The value of CS was calculated following Eq.
(3) : where D is the diffusion coefficient of gaseous sulfuric acid, β m is a transition-regime correction factor dependent on the Knudsen number (Fuchs and Sutugin, 1971), and d p i and N i are the diameter and number concentration of particles in size bin i. It has been proved that (k OH+SO 2 · [SO 2 ] · [OH])/CS is a very reliable proxy for H 2 SO 4 during the day (Lu et al., 2019). The ozonolysis of alkenes can form stabilized Criegee intermediates (SCIs) in addition to OH, and SCIs can also oxidize SO 2 to form H 2 SO 4 . A previous study on H 2 SO 4 proxies in this site has revealed that the reactions of SO 2 with products from the ozonolysis of alkenes generate a moderate amount of nighttime sulfuric acid, with little effect on daytime sulfuric acid (Yang et al., 2021). Thus, OH may be overestimated during the nighttime. In this study, OH was used to calculate the production rates of RO 2 during the daytime (Fig. 4); the error in OH does not change the relative distribution of RO 2 from different precursors.

The binPMF approach
The binPMF approach has been used to analyze the measured high-resolution (HR) mass spectrometry data. Briefly, the raw spectra were divided into narrow bins with a width of 0.006 Th after mass calibration. The data matrix and error matrix were prepared according to the methods described by Zhang et al. (2019) for the PMF model inputs (Sect. S2 in the Supplement). Differently from the traditional PMF, which uses, for example, unit mass resolution (UMR) or HR data as input, binPMF still retains HR information as much as possible, avoids the uncertainty in HR peak fitting influencing the results of PMF, and separates the complex overlapping peaks for fitting. The PMF analysis in this work uses the Igor-based analyzing interface SoFi (solution finder, version 6.8) and ME-2 as described in Canonaco et al. (2013). After selecting the PMF solution, we fitted the HR peaks in each factor through tofTools.
3 Results and discussions Figure 1 shows temporal variation in OOMs and related parameters at the SORPES station in the northeastern sub-urb of Nanjing from 2 August to 6 September 2019. During the observation period, 22 of 35 d had maximum hourly temperatures above 30 • C and 29 d had maximum hourly J (O 1 D) values above 2 × 10 −5 s −1 . High temperatures and solar radiation indicate strong photochemistry, producing a large amount of ozone, with concentrations often exceeding 80 ppb. Even at night, the concentration of ozone is rarely lower than 10 ppb, resulting from the weak titration of low NO. At the same time, the reaction between ozone and high concentrations of NO 2 can provide sufficient NO 3 radicals, dominating nocturnal degradation of certain VOCs (Wayne et al., 1991). The elevated mixing level of total aromatic hydrocarbons is one of the main characteristics of the atmosphere in densely populated areas, in addition to which there should be many alkanes and alkenes which cannot be observed by the PTR-TOF-MS Xu et al., 2017). In the daytime with strong photochemical reactions (J (O 1 D) > 1 × 10 −5 s −1 ), we instead observed concentrations of isoprene higher than those of total aromatics (Fig. S1a in the Supplement). The complex mixtures of anthropogenic and biogenic VOCs can be oxidized through a variety of pathways to produce OOMs, of which some low-volatility components will condense into particles, forming organic aerosols. The concentrations of OOMs with a mass-to-charge ratio (m/z) below 360 Th are usually higher than 10 6 molecules cm −3 , and some can even reach up to 10 7 -10 8 molecules cm −3 . Clustered peaks on the spectra of OOMs and their clear daily variations imply a lot of chemical and physical dynamics information (Figs. 1d and S1b), which is the main aspect we want to explore in this work.
The binPMF analysis was performed to characterize the sources or processes of OOMs. A 14-factor solution was selected to interpret the data set, including 3 factors about AVOC daytime chemistry, 3 isoprene-related factors, 3 factors about BVOC nighttime chemistry, 3 factors about nitrated phenols (NPs), and 2 factors excluded from the following discussion. Of these 2 disregarded factors, one is mainly composed of fluorinated contaminations (F contaminations) and the other is mainly a mixture of nitrated phenols and fluorinated contaminations (mixed contaminations). When naming these factors, we prioritize the description of dominated species or their precursors, but if the precursors are complex mixtures, our naming highlights the characteristics of the chemical processes that drive certain factors. Although this may not be the optimal PMF solution, it still separates a lot of useful information. We also stress that the urban OOM mix is unlikely to be a perfect combination of independent, unchanging factors, which is an underlying assumption in the PMF algorithm. As such, there will be no solution which is complete and perfect, but we chose a solution which was able to provide us with interesting insights. Details of the PMF diagnostics are provided in Sect. S2 in the Supplement (Figs. S2-S6). For the convenience of discussions, we have grouped these factors based on shared characteristics of certain factors. ; (c) total aromatics (benzene + toluene + C 8 aromatics + C 9 aromatics + C 10 aromatics + styrene), isoprene, and OAs; and (d) mass spectra of the nitrate CI-APi-TOF with m/z in the range of 202-404 Th.

AVOC daytime chemistry
The following daytime factors are characterized by C 6 -C 9 OOMs (Fig. 2a), considered to be derived from the oxidation of anthropogenic VOCs in this urban atmosphere, although we cannot completely exclude the presence of BVOC-derived OOMs, such as C 5 and C 10 OOMs.

Aro-OOM factor
The averaged double bond equivalent (DBE) of this factor is the largest among all factors (Table 1), with the main signals coming from compounds with DBE > 2 ( Fig. 2b) and consistent with the nature of the oxidation products of aromatics (Fig. 3a). Combined with the correlation with the production rates of OH-initiated primary RO 2 from aromatics calculated by Eq. (4) (Fig. 4), this factor is supposedly dominated by aromatic-derived OOMs (the Aro-OOM factor). This factor increases from 05:00 LT with a maximum at 10:00 LT and a sub-peak around 16:00 LT (Fig. 3e), following the diurnal variations in the P RO 2 of C 7 -C 10 aromatics ( Fig. 4b-d) but poorly correlated with the P RO 2 of benzene (Fig. 4a). Furthermore, OOMs with eight carbon atoms have the highest signal in this factor (Fig. 2a), derived from the most abundant C 8 aromatics + styrene RO 2 (Fig. 4f). Both of these findings can be explained by the fact that substituted aromatics have higher OH reactivity (Bloss et al., 2005) and higher HOM Y. Liu et al.: NO x -influenced formation of OOMs from AVOCs mixed with BVOCs yields Molteni et al., 2018) than their homologues with fewer carbon atoms. In terms of molecular formulas, the aromatic-derived OOMs have an overlap with monoterpene-derived OOMs (Mehra et al., 2020). Monoterpenes can contribute more C 10 OOMs than aromatics (P MT-RO 2 > P C 10 Aro-RO 2 ), but aromatics play a more important role in total in this factor since they provide more RO 2 in the urban atmosphere (Fig. 3f).
The main molecules of the Aro-OOM factor are summarized in Table S2 in the Supplement. The C x H 2x−5 O 6 N (x = [6, 12], of which C 8 H 11 O 6 N is the most intense) series can be produced by Reaction (R1a) of NO with the bicyclic peroxy radicals (HO−Ar−(O 2 ) 2 ), the key intermediates for aromatic oxidation proposed in the MCM (Bloss et al., 2005;Birdsall and Elrod, 2011). And here dihydroxy nitro-BTEX can be treated as an indicator of aromatic oxidation. In addition to the conventional products, C 9 H 13 O 7−9 N from the C x H 2x−5 O 7−9 N (x = [7, 13]) series is also significant in the OH-initiated and NO x -influenced oxidation experiments of 1,2,4-trimethylbenzene (Zaytsev et al., 2019) and of 1,3,5-trimethylbenzene (Tsiligiannis et al., 2019). More oxygenated compounds may come from autoxidation and multigenerational OH attacks. However, the effective OSc of this factor (Table 1) is lower than that of oxidation products of aromatics in recent laboratory studies (Zaytsev et al., 2019;Tsiligiannis et al., 2019;Wang et al., 2020c). We speculate that the abundances of NO x relative to oxidants and precursors in these experiments are not sufficient to reproduce the atmospheric conditions during our observation or that HOMs are more concentrated in aerosols due to the large condensation sink at this site . Although species with DBE < 3 ( Fig. 2b) in this factor are most likely produced from multiple OH attacks in aromatic oxidation, we cannot rule out the contribution of alkanes co-emitted with aromatics, such as the series C

Temp-related factor
This factor is named due to good correlation with temperature (Fig. 5) and shows a maximum intensity in the afternoon at around 15:00 LT (Fig. 3e). The Temp-related factor is the only one dominated by non-nitrogenous organics ( Fig. 3b and d) and has the highest effective OSc (Table 1) 5,11]), and C x H 2x−4 O 6 (x = [5, 10]) series are possibly products from RO 2 terminated by HO 2 (Reaction R2a) or closed-shell products from RO in Reactions (R3a) or (R3b). Temperature starts to rise at 06:00 LT (Fig. 12b), but this factor does not accumulate significantly until after about 10:00 LT (Fig. 3e), when the mixed level of NO is reduced to 1 ppb (Fig. 4f). This phenomenon suggests a probability of HO 2 -driven chemistry of this factor under low-NO conditions since NO can consume HO 2 and compete with HO 2 for RO 2 . Such low-NO atmospheric oxidation pathways have been suggested to be non-negligible in the afternoon in central Beijing (Newland et al., 2021).
A factor caused by similar chemical processes called the isoprene afternoon was discovered in the nitrate CI-APi-TOF data collected at a forest site in Centreville, Alabama, USA (Massoli et al., 2018), and correlated well with HO 2 , O 3 , and temperature. We also observed a number of isoprene oxidation products in the Temp-related factor (n C = 4 and (n C = 5 in Fig. 2a) 4,9]) series were also present in the light HOM factor which was supposed to be fragments from the oxidation of larger VOCs (e.g., monoterpenes) in Hyytiälä, Finland , while at the SORPES station, the C 6 -C 9 OOMs should mainly come from the oxidation of anthropogenic VOCs. At lower temperatures, the propensity of condensable organic molecules to condense into aerosols makes the concentration measured using the nitrate CI-APi-TOF lower. Thus, the total concentration of the Temp-related factor in the gas and aerosol phases was calculated based on a gas-particle equilibrium (Sect. S5 in the Supplement) and was still found to be temperature dependent (Fig. S7 in the Supplement), illustrating the temperature-influenced chemical process controlling the factor. For instance, unimolecular reaction rates like RO 2 H shifts increase qualitatively with temperature Frege et al., 2018).

Aliph-OOM factor
This factor is dominated by organic nitrates (Fig. 3c and d), and contains the bulk of anthropogenic dinitrates and trinitrates.  Wang et al., 2020c), and nor have they been reported in the forest or rural environments Massoli et al., 2018). A reasonable assumption is that these saturated or nearly saturated compounds are the products of aliphatics (including alkanes, alkenes, aliphatic alcohol, etc.) during their oxidation affected intensively by NO x in the urban atmosphere. The Aliph-OOM factor has a broad afternoon peak lasting from 14:00 to 19:00 LT (Fig. 3e), suggesting that Note that MW is the molecular weight, OSc is the carbon oxidation state, O : C is the oxygen-to-carbon ratio, N : C is the nitrogen-to-carbon ratio, DBE is the double bond equivalent, C * is the saturation concentration, and log 10 (C * ) is the volatility. the formation of multi-nitrate requires enough OH exposure time.
Considering a simple scenario of alkane photo-oxidation under high-NO x conditions, the RO 2 generated from OH attack is completely terminated by NO (Fig. 6a). The chain-retaining products are C n H 2n O (one more carbonyl group than the precursor) and C n H 2n+1 O 3 N (one more nitrate group than the precursor), and the re-oxidation of these products is a repetition of the above process which is defined as the basic reaction scheme. The multiple-generation products Figure 4. Correlations of the Aro-OOM-dominated factor with the production rate of RO 2 from OH-initiated oxidation of (a) benzene (P C 6 Aro-RO 2 ), (b) toluene (P C 7 Aro-RO 2 ), (c) C 8 aromatics (P C 8 Aro-RO 2 ), (d) C 10 aromatics (P C 10 Aro-RO 2 ), and (e) monoterpenes (P MT-RO 2 ). All the dots are colored by J (O 1 D) to show the difference between day and night. The median diurnal patterns of this factor and related parameters are plotted in (f). of alkanes summarized in Fig. 6b are regarded as reference compounds, which we compare OOMs with to investigate other mechanisms that differ from those shown in Fig. 6a and b. Specifically, this comparison is performed between the reference molecule and OOMs with the same numbers of carbon, hydrogen, and nitrogen atoms but different numbers of oxygen atoms. The number of extra oxygen atoms (n O extra ) from each aliphatic OOM over its corresponding reference molecule was calculated by Eq. (5), that is, by subtracting carbonyl and nitrate oxygens from the molecule. Thus, n O extra can represent the additional oxygenated moieties such as the hydroxyl group (−OH), peroxy group (−OOH), and possibly ether group. These functional groups may come from RO isomerization (Orlando et al., 2003), the addition of OH to alkenes, pre-existing moieties in the precursor, RO 2 autoxidation, or specific RO 2 bimolecular termination reactions (RO 2 + RO 2 , RO 2 + HO 2 ).
As shown in Fig. 6c, aliphatic OOMs in this factor are mainly the third-generation products followed by the secondgeneration products, both of which have one or two oxygen-containing functional groups in addition to the carbonyls and nitrates. It should be noted that the first-generation (Fig. 6a) and basic products (Fig. 6b) here are underestimated due to the low sensitivity of the nitrate CI-APi-TOF to these compounds. The multifunctional products of aliphatics are condensable to form SOAs (correlation coefficients with SOAs shown in Fig. 5). Recent work has showed that autoxidation is more common than previously thought , and more studies are needed to explore the oxidation mechanisms of anthropogenic aliphatics and to evaluate their contribution to SOAs.

Photo-related factor
This factor is defined based on its correlation with J (O 1 D) (Fig. 5), having an apparent diurnal cycle with a peak at 12:00 LT (Fig. 7e). The major peak of the Photo-related factor is C 5 H 10 O 8 N 2 (Fig. 7a), most probably generated from double OH attack proceeding with double RO 2 + NO termination (Jenkin et al., 2015). C 5 H 10 O 8 N 2 can be also produced in an NO 3 + isoprene system (Ng et al., 2008;Zhao et al., 2021), whereas in this study, the nocturnal C 5 H 10 O 8 N 2 is principally from the isoprene-derived OOM (Isop-OOM) factor (Fig. 8b) which will be discussed later. Other peaks with n C ≤ 5, like C 5 H 7 O 7 N, C 4 H 7 O 6 N, and C 5 H 9 O 6 N, are also likely to be the isoprene products. The total signal of compounds with n C > 5 is not low, although their respective proportions are not as prominent as those of C 5 species (Fig. 7d), implying the contribution of other precursors together with isoprene. In addition, the relationship of this factor with isoprene and J (O 1 D) together (Fig. 5) reveals the effect of light-dependent emission of isoprene on it.

O x -and SOA-related factor
The atmospheric oxidation of VOCs produces low-volatility compounds, forming SOAs through gas-particle partitioning, and concurrently promotes ozone formation (Atkinson, 2000). Both SOAs and O x have long lifetimes (> 12 h), and their correlations have been extensively investigated (Herndon et al., 2008;Wood et al., 2010;Hu et al., 2016). The OOM factor related to ozone and SOAs together (Fig. 5), having slightly elevated concentrations during the daytime (Fig. 7e), is considered to be generated from this photochemical aging process. Apart from C 5 H 10 O 8 N 2 , other isoprene multi-nitrates are also present in this factor. C 5 H 9 O 10 N 3 , an isoprene hydroxyl trinitrate requiring at least two steps of oxidation found in the experimental study on isoprene oxidation by NO 3 (Zhao et al., 2021), of course does not appear in the Photo-related factor at all but is mostly apportioned into the O x -and SOA-related factor and the Isop-OOM factor ( Fig. 8c and d). Like the Photo-related factor, isoprene is a significant but not the only precursor of this factor (Figs. 2  and 7). The biggest peak of the O x -and SOA-related factor is an ion at m/z 264 with the formula C 6 H 5 O 3 N (HNO 3 NO 3 − ), identified as an adduct of nitrophenol (C 6 H 5 O 3 N) with nitrate dimmer (HNO 3 NO 3 − ). The time variation in C 6 H 5 O 3 N (HNO 3 NO 3 − ) is influenced by the reagent ions in addition to the atmospheric nitrophenol. So far, we do not know why this compound shares the same processes with others, but we performed a test that removed the bins with the unit m/z = 264 from the input matrix and still obtained this factor from the PMF model.

Isop-OOM factor
The mass spectra of the Isop-OOM factor, as its name implies, are exclusively contributed by isoprene-derived compounds (Fig. 7c). C 5 H 10 O 8 N 2 contributes about 30 % of the intensity of this factor, and the dominance of C 5 H 10 O 8 N 2 was also found in the isoprene nitrate type-I factor in Centreville (Massoli et al., 2018). In addition to multi-nitrates (C 5 H 10 O 7−8 N 2 , C 5 H 8 O 6−9 N 2 , and C 5 H 9 O 10 N 3 summarized in Table S6 in the Supplement), several mononitrate series (C 4 H 7 O 5−7 N, C 5 H 9 O 4−9 N, C 5 H 7 O 5−8 N, and C 5 H 11 O 5−6 N) of this factor are also abundant in the isoprene nitrate type-II factor in Centreville (Massoli et al., 2018). Many of the isoprene nitrates here have been especially investigated in our previous observations in the YRD (Xu et al., 2021) and have been discovered in other field measurements (Lee et al., 2016;Massoli et al., 2018) and in many laboratories (Ng et al., 2008;Lambe et al., 2017). Generally, these compounds are second-and third-generation OH oxidation products of isoprene under high-NO x conditions (Wennberg et al., 2018).
The diurnal pattern of the Isop-OOM factor is relatively unclear (Fig. 7e), with obvious differences between mean and median values usually caused by plume events. This indicates that isoprene chemistry, usually varying evidently from day (OH-initiated) to night (NO 3 -initiated), is not the driver of this factor. This factor correlates positively with MVK-MACR and SOAs (r > 0.50, shown in Fig. 5) but not with isoprene and OH. It seems that these isoprene OOMs are produced elsewhere and then transported due to their longer lifetime determined by their relatively high volatility (Table 1). The Isop-OOM factors in the continental air masses are more intense than those in the coastal and YRD air masses (Fig. S8 in the Supplement), consistent with the spatial distribution of isoprene emissions (Sindelarova et al., 2014). An archetypal episode affected by continental air masses (13 August to 17 August 2019) is shown in Fig. 8. During this period, C 5 H 9 O 10 N 3 was almost entirely transported, while C 5 H 10 O 8 N 2 had strong in situ photochemical generation, in addition to the source of transport.

BVOC nighttime chemistry
The following nighttime factors are characterized by C 10 OOMs (Fig. 2a), which are identified as the oxidation products of monoterpenes. Except for the BVOC-OOM-I fac- tor (Fig. 9a), the contribution of isoprene-derived OOMs was much lower in these factors. Compared to the above isoprene-related factors, C 5 H 10 O 8 N 2 and C 5 H 9 O 10 N 3 were no longer significantly present in the following factors.

BVOC-OOM-I factor
The first nighttime factor has its maximum concentration at around 20:00 LT and decreases to very low values during the day. It is moderately correlated with the production rate of NO 3 radical (P NO 3 derived from Eq. 6) at night and reaches high intensity only under conditions of NO below 1 ppb (Fig. 10a), indicating a chemical process of NO 3 radical. The concentration of this factor is mainly from C 5 peaks, followed by C 6 -C 10 peaks (Fig. 9d), about 80 % of which are organic nitrates (Fig. 2c), designating the oxidations of isoprene and monoterpenes by NO 3 (BVOC-OOM-I). In the case of isoprene oxidation, the nitrate groups of C 5 H 9 O 4−8 N, C 5 H 7 O 5−8 N, and C 4 H 7 O 5−6 N series (summarized in Table S8 in the Supplement) are likely to come from the addition of NO 3 . Next, the C 5 H 10 O 8−9 N 2 and C 5 H 8 O 7−10 N 2 series are probably second-generation products. These compounds derived from an isoprene + NO 3 system have been discussed in previous work with laboratory (Kwan et al., 2012;Zhao et al., 2021) and ambient (Ayres et al., 2015;Xiong et al., 2015) data sets. Additionally, the C 6 -C 10 species are potentially the products of monoterpenes degraded by NO 3 .

BVOC-OOM-II factor
The second nighttime factor is intense at night and over 5 times lower during the day. Like the BVOC-OOM-I factor, this factor has high concentrations when NO is reduced to increase NO 3 availability (Fig. 10b), and about 80 % of compounds in this factor are organic nitrates (Fig. 2c). Accordingly, this may also be a factor strongly influenced by NO 3 . It is dominated by C 6 -C 10 OOMs, among which the highest intensity is at C 10 (Fig. 9d), coherent with the nature of monoterpene products (BVOC-OOM-II). This factor has weaker signals at C 15 which are plausibly the products of sesquiterpenes but could also be dimmers formed from Reaction (R4) (monoterpenes + isoprene or monoterpenes + C 5 monoterpenes fragments). Compared to the BVOC-OOM-I factor (Fig. 9d), this factor has more large mass molecules (C 10 ) and fewer small mass molecules (C 5 ), resulting in an effective volatility over 1 order of magnitude lower. A NO 3 -initiated factor, called the nighttime type-2 factor, has also been discovered in Hyytiälä, Finland , but the similar factor we found has a higher proportion of organic nitrates, due to the more abundant NO x here.
C 10 H 16 + NO 3 q O 2 − → C 10 H 16 NO 5 q → C 10 H 16 NO 4 q (Alkoxy) In terms of the fingerprint molecules of this factor (summarized in Table S9 in the Supplement), the C 10 H 15 O 5−12 N series are carbonyl products from precursor RO 2 or RO terminations, while the C 10 H 17 O 5−9 N series are alcohol or hydroperoxide products from precursor RO 2 terminations. The C 7 H 9 O 6−8 N, C 9 H 15 O 6−9 N, C 9 H 13 O 7−10 N, and C 8 H 13 O 7−8 N series are expected to be fragments. The closed-shell compounds mentioned above have been reported in the experiments of a monoterpene + NO 3 system (Nah et al., 2016;Faxon et al., 2018;Takeuchi and Ng, 2019).
It is noteworthy that a set of nitrogen-containing radicals, C 10 H 16 O 6−11 N (peak fittings are shown in Fig. S9 in the Supplement), is present in the BVOC-OOM-II factor. Starting from a generic monoterpene molecule with the formula C 10 H 16 , the NO 3 addition with fast O 2 addition results in a peroxy radical with the formula C 10 H 16 O 5 N. If the initial intermediate RO 2 is capable of proceeding via autoxidation by the formal addition of O 2 , we expect radicals, C 10 H 16 O 5+2x N (x denotes times of autoxidation performed) with an odd oxygen number, to be formed (Reaction R5a). In addition, peroxy radicals with an even oxygen number, C 10 H 16 O 6+2x N, are likely produced via a reaction chain (R5b): (1) RO 2 is propagated to RO through bimolecular reactions, and (2) RO isomerizes to an alcohol by internal H abstraction forming a carbon-centered radical (Orlando et al., 2003;Orlando and Tyndall, 2012); (3) the carbon-centered radical can again take up an oxygen molecule and follow the autoxidation route. The C 10 H 16 O 9 N radical is also moderately intense in the BVOC-OOM-I factor (Fig. 9a), testifying to the presence of NO 3 chemistry. These C 10 H 16 O 6−11 N rad- icals are also reported in the CLOUD chamber . In addition to C 10 radicals, a C 5 radical, C 5 H 8 O 5 N (peak fittings are shown in Fig. S9), is also found in the BVOC-OOM-II factor. C 5 H 8 O 5 N is possibly derived from the oxidation of isoprene initiated by NO 3 , as observed in the laboratory (Zhao et al., 2021). Another hypothesis is that C 5 H 8 O 5 N is formed from the fragmentation process of monoterpenes.

BVOC-OOM-III factor
The third nighttime factor (BVOC-OOM-III) is dominated by nitrogen-containing species with a contribution ratio of about 90 %, among which dinitrates account for more than 20 % (Fig. 2c). When grouped by carbon numbers, C 10 OOMs have the strongest signal. Unlike the above two NO 3related factors, this factor shows no correlation with P NO 3 but has positive correlation with NO, especially during the daytime (Fig. 10h). C 9 H 14 O, a typical product of NO-affected monoterpene oxidation (Calogirou et al., 1999), is found to be correlated with this factor (Fig. 10i). It is reasonable to infer that these organic nitrates may come from termina-tions of monoterpene RO 2 by NO. In addition to the elevated intensity during the nighttime, this factor still remains at a relatively high concentration in the morning, a concentration which is much higher than those of the two NO 3related factors (Fig. 9e). Owing to the suppression of RO 2 autoxidation by NO and the relatively insufficient oxidant in dark environments, the effective OSc of the BVOC-OOM-III factor is lower than those of other factors. Apart from the mononitrates summarized in Table S10 in the Supplement, the C 10 H 16 O 7−10 N 2 (dinitrates) and C 10 H 17 O 10 N 3 (a trinitrate charged by NO 3 − at m/z 401) are most likely the result of multiple-generation processes involving OH or NO 3 oxidation of monoterpenes and RO 2 + NO terminations. A similar factor, called terpene nitrates, has also been reported in Centreville, USA (Massoli et al., 2018), while in Hyytiälä, Finland , the daytime type-1 factor is related to NO.

Nitrated-phenol factors
Nitrated phenols are of concern because of their phytotoxicity (Rippen et al., 1987) and as important chromophores of brown carbon in aerosols (Desyaterik et al., 2013;Mohr et al., 2013). The sources of these highly volatile compounds are attributed to biomass burning, vehicle exhausts, and secondary gas-phase or aqueous-phase production (Harrison et al., 2005). Here we identified three factors about NPs, including the NP-I factor dominated by nitrophenol, the NP-II factor dominated by substituted nitrophenols, and the NP-III factor dominated by dinitrophenols. Although the mass spectrum of the NP-III factor is less pure than those of the NP-I and NP-II factors (Fig. 11), its time series follows well with C 6 H 4 O 5 N 2 (Fig. 11f), implying that this factor is driven by dinitrated-phenol chemistry. Since nitrated phenols have been broadly investigated and relatively clearly recognized (Harrison et al., 2005;Yuan et al., 2016;Wang et al., 2018b;Cheng et al., 2021), they are not discussed too much here. It seems that the chemistry of nitrated phenols is distinctive from that of other OOMs.

Ensemble chemical properties
After performing PMF analysis, over 1000 non-nitro molecules have been identified through HR peak fitting for each factor. The mean concentration of total non-nitro OOMs reconstructed from the selected PMF solution is about 2.1 × 10 8 molecules cm −3 . Ensemble chemical properties of these non-nitro OOMs are summarized in Fig. 12. The number of carbon atoms implies the precursor information of OOMs. C 5 OOMs, which principally consist of isoprene products benefiting from the high reactivity and intensive emissions of isoprene in summer, are the most abundant (Fig. 12c), while C 6 -C 9 OOMs are mostly likely formed from the oxidation of AVOCs such as aromatic and aliphatic series in the urban and suburban atmosphere, and as we expected, these AVOCderived OOMs account for about 50 % of the total signal (Fig. 12c). The intensity of OOMs decreases from C 7 to C 9 determined by the concentration distribution of precursors but becomes a plateau at C 10 ( Fig. 12c), indicating another source of C 10 OOMs, such as monoterpene oxidation. These results underscore the formation of SOA precursors from a mixture of anthropogenic and biogenic emissions, under ongoing forest cover increases (Wang et al., 2020a) in highly urbanized eastern China.
In addition to the anthropogenic VOCs, another humaninduced perturbation on the formation of OOMs is the NO xaffected chemistry of VOCs, i.e., RO 2 + NO terminations or NO 3 -initiated oxidations. As shown in Fig. 12c, about 72 % of OOMs are nitrogen-bearing compounds, regarded as organic nitrates within the allowed range of uncertainty. If isoprene nitrates are not included, organic nitrates peak at C 7 as do the nitrogen-free species, showing the significant production of organic nitrates through the AVOCs + NO x pathways. The NO x effect on AVOC-derived OOMs, typified by the Aro-OOM factor and the Aliph-OOM factor, is not shown in previous ambient measurements Lee et al., 2016;Massoli et al., 2018).
OOMs grouped by carbon numbers or nitrogen numbers consistently have high absolute concentrations in the daytime ( Fig. 12a and b), revealing the crucial role of photochemical progress, involving RO 2 + NO termination reactions, in OOMs. In addition, The C 5 -C 10 OOMs are enhanced again during 19:00-22:00 LT, and the nighttime peak of C 10 OOMs is even higher than their daytime peak (Fig. 12a). The nocturnal C 10 OOMs are more intense than C 9 OOMs (Fig. 12a), and there are more C 10 nitrates than C 9 nitrates (Fig. 12c). These results show the fate of VOCs degraded by NO 3 during the nighttime, and it is more important to monoterpenes. In contrast to nitrogen-free OOMs, organic nitrates are enriched through the reactions of BVOCs with NO 3 in the early evening (Fig. 12b), as indicated by three BVOC nighttime chemistry factors.
Apart from reflecting the influence of NO x , multi-nitrates also imply the multiple generations of VOC oxidation, and multiple generations are evident in the products of isoprene (e.g., C 5 H 10 O 8 N 2 and C 5 H 9 O 10 N 3 ) due to its two carbon-carbon double bonds. As products of mononitrates, multi-nitrates follow mononitrates in their diurnal variation, with double peaks initiated by OH and NO 3 respectively (Fig. 12b). Considering that the formation of organic nitrate is only a small branch of RO 2 +NO termination, the contribution of multi-step oxidation should be larger than that shown in Fig. 12c.

Conclusions
We have investigated the sources and characteristics of gasphase OOMs observed using a nitrate CI-APi-TOF at the SORPES station in the YRD of eastern China, an environment dominated by anthropogenic emissions with enhanced biogenic emissions during summer.
The binPMF analysis, which avoids the uncertainty introduced by high-resolution peak fitting to the input data matrix, was applied to deconvolve the complexity of the data set, and it resolved 14 factors, among which 12 factors have been discussed in detail. A morning factor (Aro-OOM), correlated with the production rates of RO 2 from aromatics, is characterized by unsaturated products of aromatics such as C x H 2x−5 O 6−9 N (x = [6,12]). An afternoon factor (Aliph-OOM), containing the bulk of C 6 −C 9 dinitrates and trinitrates such as C x H 2x−2 O 8 N 2 (x = [4, 13]) and C x H 2x O 8 N 2 (x = [4,8]), is assumed to be derived from aliphatic oxidation. A transported factor (Isop-OOM) correlates with MVK-MACR and SOAs and is exclusively dominated by isoprene nitrates (e.g., C 5 H 10 O 8 N 2 and C 5 H 9 O 10 N 3 ). A nighttime factor (BVOC-OOM-III), related to NO, is dominated by terpenes nitrates such as C 10 H 15 O 6 N, C 10 H 16 O 7−10 N 2 , and C 10 H 17 O 10 N 3 . In addition to the factors distinguished by precursors, several factors are driven by chemistry. A factor following J (O 1 D) (Photo-related), consisting of isoprene products mixed with others, is thought to be produced by in situ photochemistry. An afternoon factor (Temp-related), having the most abundant nitrogenfree OOMs such as C x H 2x−4 O 5−6 (x = [5, 10]), C x H 2x−2 O 5 (x = [5, 10]), and C x H 2x−6 O 4 (x = [5, 10]), is generated involving temperature-influenced chemistry. A daytime factor (O x -and SOA-related), correlated well with O x and SOAs, indicates the photochemical aging process. Two nighttime factors (BVOC-OOM-I and BVOC-OOM-II), benefiting from NO 3 and suppressed by NO, are considered to be produced from the NO 3 -initiated oxidation of BVOCs, and both of them have the fingerprint molecule, C 10 H 16 O 9 N. The remaining three factors are governed by nitrated phenols.
All of these factors from various precursors are influenced in different ways by NO x . Over 1000 non-nitro molecules have been identified and then reconstructed from the selected solution of binPMF, and about 72 % of the total signal is contributed by nitrogen-containing OOMs, almost regarded as organic nitrates formed through RO 2 + NO terminations or NO 3 -initiated oxidations. Moreover, multi-nitrates have a contribution ratio of about 24 % to total concentration, indicating the significant presence of multiple oxidation generations, especially for isoprene (e.g., C 5 H 10 O 8 N 2 and C 5 H 9 O 10 N 3 ). The nitrate CI-APi-TOF data set presented here highlights the decisive role of NO x chemistry on OOM formation in densely populated areas. The distributions of total observed OOMs at different n C values. OOMs on each carbon number are grouped by nitrogen number, and the total concentration fractions of each groups are reported in the pie chart. Since we selected peaks in the m/z range of 202-404 Th, OOMs with n C < 5 or n C > 10 detected by the nitrate CI-APi-TOF are underestimated.
The differences in OOMs observed in different environments are so clear, and the underlying causes for this are well worth considering. The precursors, oxidants, and formation pathways of OOMs change when moving from urbanized areas to pristine regions, as AVOC and NO x concentrations decrease and BVOC concentrations increase. This process can also occur under the trend of global warming and anthropogenic emissions mitigation, but we still know very little about it. Clarifying the variations in compositions, properties, and the formation efficiency of OOMs will help us to understand the evolution of SOA production during this process. In summary, our findings highlight the dramatic interactions between anthropogenic and biogenic emissions and encourage more investigations from a mechanistic point of view.
Appendix A: The selected solution for binPMF analysis of nitrate CI-APi-TOF data Figure   Data availability. Measurement data at the SORPES station, including OOM data and relevant trace gas and aerosol data as well as meteorological data, are available upon request from the corresponding author before the SORPES database is open to the public.
Author contributions. WN and AD designed this research. YL, YL, DG, CL, ZX, LC, TW, LW, PS, XQ, JW, and XC conducted the measurements. YL, WN, CY, YZ, DH, ZW, and DW analyzed the data. YL, WN, ND, ME, and AD wrote the paper.
Competing interests. Some authors are members of the editorial board of Atmospheric Chemistry and Physics. The peer-review process was guided by an independent editor, and the authors have also no other competing interests to declare.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.