Chemical characterization of oxygenated organic compounds in the gas phase and particle phase using iodide CIMS with FIGAERO in urban air

The atmospheric processes under polluted environments involving interactions of anthropogenic pollutants and natural emissions lead to the formation of various and complex secondary products. Therefore, the characterization of oxygenated organic compounds in urban areas remains a pivotal issue in our understanding of the evolution of organic carbon. Here, we describe measurements of an iodide chemical ionization time-of-flight mass spectrometer installed with a Filter Inlet for Gases and AEROsols (FIGAERO-I-CIMS) in both the gas phase and the particle phase at an urban site in Guangzhou, a typical megacity in southern China, during the autumn of 2018. Abundant oxygenated organic compounds containing two to five oxygen atoms were observed, including organic acids, multi-functional organic compounds typically emitted from biomass burning, oxidation products of biogenic hydrocarbons and aromatics. Photochemistry played dominant roles in the formation of gaseous organic acids and isoprene-derived organic nitrates, while nighttime chemistry contributed significantly to the formation of monoterpene-derived organic nitrates and inorganics. Nitrogen-containing organic compounds occupied a significant fraction of the total signal in both the gas and particle phases, with elevated fractions at higher molecular weights. Measurements of organic compounds in the particle phase by FIGAERO-I-CIMS explained 24± 0.8 % of the total organic aerosol mass measured by aerosol mass spectrometer (AMS), and the fraction increased for more aged organic aerosol. The systematical interpretation of mass spectra of the FIGAERO-I-CIMS in the urban area of Guangzhou provides a holistic view of numerous oxygenated organic compounds in the urban atmosphere, which can serve as a reference for the future field measurements by FIGAERO-I-CIMS in polluted urban regions.


Introduction
In urban air, atmospheric chemical processes are varied and complex, as the result of large emissions of both anthropogenic pollutants and biogenic volatile organic compounds, associated with strong interactions with each other (He et al., 2014;Karl et al., 2018;Shrivastava et al., 2019). Consequently, strong formation of secondary pollutants, e.g., ozone and secondary organic aerosol (SOA), are observed in urban and downwind regions . Oxygenated organic compounds are not fully accounted for in some earlier studies, which may explain some of the discrepancies between observations and models for many unaddressed issues in atmospheric chemistry. Oxygenated organic compounds are supposed to be the top candidates for missing OH reactivity observed in various environments including pristine rainforests and urbanized areas (Noelscher et al., 2016;Yang et al., 2016Yang et al., , 2017. The photolysis of carbonyls serves as a critical radical source driving ozone formation in highly polluted industrialized areas (Edwards et al., 2014;Liu et al., 2012;Xue et al., 2016). Although it was discovered a long time ago that oxygenated organic compounds make up a substantial fraction of submicron aerosol mass (Kroll and Seinfeld, 2008), enormous difficulty still exists in accurately predicting the formation and evolution of SOA (de Gouw et al., 2005;Hodzic et al., 2010;Volkamer et al., 2006).
One of the biggest obstacles to understanding the role of oxygenated organic compounds is the characterization of these extremely complicated and diverse chemicals which encompass tens of thousands of individual species spanning a wide range of volatility. Chemical ionization mass spectrometry (CIMS) is a powerful technique for the molecular-level characterization of oxygenated organic compounds because of the following advantages (Zhao, 2018): direct measurements and fast time response to capture the rapid temporal change of short-lifetime intermediates, soft ionization providing chemical information on a molecular level, and selective ionization ensuring measurements for specific classes of species. The iodide anion ionizes species mainly through adduction  and is used for the detection of oxygenated organic compounds, particularly organic compounds with two to five oxygen atoms (Lee et al., 2014;Lopez-Hilfiker et al., 2016;Riva et al., 2019). It has been shown that I-CIMS is an excellent technique to investigate the oxidation processes of volatile organic compounds (VOCs) and formation of SOA . Installed with a thermal desorption inlet that collects and heats aerosol to evaporate organic compounds, e.g., Filter Inlet for Gases and AEROsols (FIGAERO, Lopez-Hilfiker et al., 2014) and Micro-Orifice Volatilization Impactor (MOVI, Yatavelli et al., 2012), the CIMS instruments are capable of analyzing particle-phase species and gas-particle partitioning in a semicontinuous way Stolzenburg et al., 2018).
Although FIGAERO-CIMS has gained recent popularity in atmospheric chemistry research, much of the published work was done in chambers or in the laboratory (D'Ambro et al., 2017Hammes et al., 2019;Lopez-Hilfiker et al., 2015). As for the applications in field campaigns, most work has been mostly performed in forest or rural areas (Huang et al., 2019;Hunter et al., 2017;Lee et al., 2016Lee et al., , 2018b; measurements in the urban atmosphere by FIGAERO-CIMS are still limited (Le Breton et al., 2018b). Meanwhile, a sys-tematic analysis on mass spectra of FIGAERO-CIMS in the ambient air is imperative, for a more holistic view in investigating emissions and chemistry of oxygenated organic compounds. In this study, we present the measurement results using FIGAERO-I-CIMS during a coordinated campaign in Guangzhou, a megacity in the Pearl River region of China. We describe the experimental design, instrumentation setup, calibration and data processing for this instrument in the campaign. This work will provide a detailed interpretation of the mass spectra of oxygenated species in both the gas phase and the particle phase. The bulk chemical properties will also be discussed to provide an overview of organic compounds.

Measurement site and supporting data
Measurements were conducted during the coordinated campaign "Particles, Radicals and Intermediates from oxiDation of primary Emissions over the Great Bay Area" (PRIDE-GBA) in October and November 2018. The Great Bay Area (GBA) refers to a highly industrialized and urbanized area in southern China, including the two special administrative regions of Hong Kong and Macau, and nine cities surrounding the Pearl River estuary. Affected by the subtropical monsoon climate, the weather in the region was characterized by high temperatures and relative humidity (RH) as well as sufficient sunshine (total solar radiation of the Pearl River Delta region in the fall of 2016 was ∼ 1200 MJ/m 2 , Liu et al., 2018). The city of Guangzhou lies in the north of the GBA and south of the mountains. Therefore, the city is extensively influenced by both anthropogenic and biogenic emissions. The urban site was located at Guangzhou Institute of Geochemistry, Chinese Academy of Sciences (23.14 • N, 113.36 • E). Online instruments sampled from inlets set up in laboratories on the eighth floor or ninth floor (about 25 m above the ground).
In addition to FIGAERO-I-CIMS discussed later, measurement data from a suite of other instruments were also used in this work. A high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS, Aerodyne Research, Inc.) was deployed to provide chemical composition and many other parameters of ambient aerosol including f 60, liquid water content (LWC), particulate organic nitrate and elemental ratios (Hu et al., , 2018. The parameter f 60 is the ratio of the integrated signal at m/z 60 to the total signal of organic components and is used as a tracer for biomass burning emissions . LWC of aerosol was taken as the sum of water contributed by inorganic components predicted by the ISORROPIA II model and organic components calculated based on the organic hygroscopicity parameter (Fountoukis and Nenes, 2007;Guo et al., 2015). Based on AMS data, organic nitrate concentrations were determined by 2-3 times lower NO + 2 / NO + ratios 8457 for organic nitrate than inorganic nitrate (Fry et al., 2013). The calculation method of elemental ratios based on AMS data has been described elsewhere (Aiken et al., 2007;Canagaratna et al., 2015). Detailed information about AMS measurements from the PRIDE-GBA campaign is forthcoming in a separate paper. An online GC-MS/FID (Wuhan Tianhong Instrument Co., Ltd) and a proton transfer reaction time-of-flight mass spectrometer (PTR-ToF-MS, IONICON Analytik GmbH) (Yuan et al., 2017) served as the analytical techniques for measuring isoprene and other VOCs (e.g., monoterpenes, aromatics and a few oxygenated VOCs) , respectively. Trace gases (CO, O 3 , NO and NO 2 ) were measured by commercial gas monitors (Thermo Fisher Scientific Inc.) (Z. . Photolysis rates were measured by a PFS-100 photolysis spectrometer (Focused Photonics Inc.). Temperature and RH were measured by a Vantage Pro2 weather station (Davis Instruments Corp.). Time series and diurnal profiles of meteorological parameters, trace gases, the photolysis rate of NO 2 (j NO 2 ) and several important VOCs (isoprene, monoterpenes, toluene and benzene) are shown in Fig. S1. The temperature during the campaign was between 17 and 33 • C, with an average of 24 • C, and RH was between 27 % and 97 %, with an average of 70 %.

Experimental setup
Our instrument consists of a Filter Inlet for Gases and AEROsols (FIGAERO) and a time-of-flight chemical ionization mass spectrometer coupled with an iodide ionization source Lee et al., 2014;Lopez-Hilfiker et al., 2014). The FIGAERO is a multi-port inlet assembly following a two-step procedure alternating between gas mode, in which online measurements of gases and semicontinuous sampling of particle-phase species are conducted, and particle mode, in which particulate composition is investigated via thermal desorption Thornton et al., 2020). The iodide source is a "soft" ionization technique with little ionization-induced fragmentation and selective detection towards multi-functional organic compounds, providing elemental compositions for thousands of oxygenated compounds in the atmosphere (Hyttinen et al., 2018;Iyer et al., 2016;Lee et al., 2014;Riva et al., 2019). The sample air was drawn into the ion molecule reaction (IMR) chamber where it intersected and reacted with primary ions generated by flowing 2 mL/min 1000 ppm methyl iodide in 2.4 L/min N 2 through an X-ray source. The pressure in the IMR chamber was maintained at 370-390 mbar. Equipped with a long time-of-flight mass analyzer, our instrument was configured to measure singularly charged ions up to 603 Th with a mass resolving power of 10 000-11 000 (m/ m at 50 % height) during the campaign (Fig. S2).
Ambient air was continuously sampled through two inlets protruding about 1.5 m out of a window on the ninth floor of a building. One was a 3 m PFA tubing (1/4 in. o.d.) for gas-phase sampling, through which roughly 9 L/min air was drawn, and 2 L/min was directly taken into the instrument for gas measurements without removing particles, resulting in an inlet residence time of 0.24 s. The gas sampling line inside the room was covered by heat insulation associated with a heating cable to minimize condensation on the tubing surface. The other inlet for the particle phase was a 3.8 m metal tubing (3/8 in. o.d.) fitted with a PM 2.5 cyclone and a Nafion dryer (Perma Pure, model PD-07018T-12MSS) to reduce water content in the sampled air. The particle-phase inlet was drawn by a laminar flow at ∼ 8 L/min (Reynolds number of ∼ 1500), 3.8 L/min of which was collected on PTFE membrane filters (Zefluor®, Pall Inc., USA). The residence time was 1.3 s for the particle-phase sampling line. Semi-volatility and low-volatility compounds tend to interact with wall surfaces of both inlets and the IMR and thus extend response time . As accurate correction for wall losses remains impossible, no wall loss correction was performed in this study.
The FIGAERO worked in a cyclical 1 h pattern with two modes (Fig. S3): (i) measuring gas for the first 24 min while simultaneously collecting particles on the filter and (ii) then analyzing the particle-phase collection for another 36 min. In every 24 min gas mode, ambient air was measured for the first 21 min, followed by 3 min gas background by overflowing zero air at 5 L/min through a pinhole just in front of the IMR. The background measurements are inevitably influenced by wall interactions, especially for "sticky" species. Recently, Palm et al. (2019) proposed a new way to determine gas background ("fast background") by fast switching between the ambient air and background, which greatly improves accurate determination of the CIMS background. In the remaining 36 min, the components of the collected particles were thermally desorbed and introduced into the CIMS with 2 L/min N 2 carrier gas. The N 2 flow was ramped from ambient temperature to 175 • C in 12 min and held for another 20 min. A schematic diagram of working modes and temperature profile of FIGAERO heating in a single cycle are shown in Fig. S4. Particle background was determined every sixth 1 h running cycle in which ambient air passed over a filter (Parker Balston, model 9922-11-CQ) in front of the FIGAERO filter.

Calibration experiments
Using various techniques, we calibrated dozens of chemical compounds in the laboratory. Table S1 summarizes the calibrated species and corresponding calibration methods. (1) Gas cylinders are commercially available for a few species (e.g., chlorine, hydrogen cyanide). The gaseous standards were diluted down to different concentrations and then introduced to the CIMS. (2) For those VOCs whose stan-dards are liquid or solid, solutions with known concentrations are made and then vaporized using the liquid calibration unit (LCU, IONICON Analytik GmbH) to provide gaseous standards. (3) Commercial permeation tubes are available for some species (e.g., nitric acid). (4) Some gaseous chemicals were generated in the laboratory. For example, isocyanic acid was generated from thermal decomposition of cyanuric acid in a diffusion cell (Li et al., 2021;, and dinitrogen pentoxide was generated via the reaction of ozone with excess nitrogen dioxide in a flow reactor (Bertram et al., 2009). (5) Compounds of low vapor pressure were calibrated through the FIGAERO . Briefly, certain amounts of target species dissolved in organic solvents (e.g., isopropanol or acetone) were deposited onto the PTFE filter of the FIGAERO using a syringe, and the droplet was then subjected to a temperature-programmed thermal desorption by N 2 gas. The sensitivity was determined as the integrated signals under thermogram profiles versus the amounts of deposited calibrant.
In addition to sensitivity calibration, the effect of humidity on the sensitivity for various species was investigated in the laboratory, some of which are shown in Fig. S5. Lowmolecular-weight acids, e.g., formic acid and nitric acid, tend to be more sensitive to the humidity changes than multi-functional compounds. A similar tendency of multifunctional compounds associated with less humidity dependence was also reported in previous work (Lee et al., 2014). Considering water vapor pressure in the IMR, our humiditydependent curves are generally consistent with those reported in Lee et al. (2014) (see detailed discussions in Sect. S3 in the Supplement).
In the later part of the campaign (after 22 October), an isotopically labeled formic acid (DCOOH, Cambridge Isotope Laboratories, Inc.) permeation tube held at constant temperature (65 • ) was mixed with 10 mL/min N 2 and continuously delivered into the entrance of the sampling inlet in order to derive a humidity dependence function from the field measurements. DCOOH signals during the campaign exhibited a humidity-dependent curve consistent with formic acid obtained in the laboratory (Fig. S5). We applied humidity correction to the species with the humidity-dependent curves determined in the laboratory (underlined species in Table S1). For other compounds, humidity correction was not applied, as there is no universal pattern of humidity dependence for all detected species, and multi-functional compounds that comprise the majority of the species measured by FIGAERO-I-CIMS are usually less influenced by humidity.
The measured concentration of DCOOH was steady after humidity correction was applied (Fig. S6g), indicating the stability of our instrument. In addition, we also performed field calibrations throughout the campaign to check the instrument status by spotting a solution mixture of levoglucosan, heptaethylene glycol and octaethylene glycol onto the FIGAERO filter every 2-3 d (Fig. S6). Multiple-point calibrations for these organic species were performed in the be-ginning and the end of the campaign. The concentration of the solution used in the first two calibration experiments was too high, so we prepared a new solution for calibrations in November. The relative changes of the determined calibration factors in November were within 50 % for the calibrated species.

Data processing
The TofWare software (version 3.0.3; Tofwerk AG, Switzerland) was used to conduct the high-resolution peak fitting for the mass spectra data of ToF-CIMS, including mass calibration, instrumental parameter optimization (peak shape and peak width) and bunch fitting of high-resolution peaks (Stark et al., 2015). In this study, the signals of ions were normalized to the sum signals of I − and H 2 OI − at 10 6 cps. Hourly particle-phase data were obtained by integrating the signals of various ions during each FIGAERO desorption period. Background-corrected signals were obtained by subtracting linearly interpolated background signals from ambient signals (and integrated signals) for ions in the gas (and particle) phase.
In order to determine the sensitivities of uncalibrated species, a voltage scanning procedure was performed from time to time throughout the campaign covering different times of the day Lopez-Hilfiker et al., 2016). Here, we selected four representative periods including morning, afternoon, evening and night on polluted days. By performing sigmoidal fitting on the remaining signals as a function of voltages, a dV 50 value of each ion from each period was determined, at which voltage half of one kind of ion dissociated . We observed a positive correlation between the sensitivities of the ions relative to maximum sensitivity and their average dV 50 values ( Fig. S7), consistent with previous studies Lopez-Hilfiker et al., 2016). This relationship was used to calculate response factors for uncalibrated species, after taking into account the relative transmission efficiency for the ions (see Sect. S1 in the Supplement for detailed analysis).

Overview of detected species in the mass spectra
We identify 1334 ions adducted with iodide from the mass spectra, among which 427 are charged closed-shell organic compounds containing only C, H and O elements (C x H y O z I − ) and 388 are charged closed-shell organic compounds containing C, H, O and N elements (C x H y N 1,2 O z I − ). For species with the formula of C x H y O z , x ranges from 1 to 20, y is an even number and no more than 2x + 2, and z is greater than or equal to 2. The range of carbon number x for the ions with C x H y N 1,2 O z is the same as the ions with C x H y O z . For species containing one nitrogen (C x H y NO z ), y is an odd number and less than 2x + 2; z is larger than or equal to 2. For species containing two nitrogen atoms (C x H y N 2 O z ), y is an even number and less than 2x + 1; z is larger than or equal to 4. Table 1 summarizes species discussed in the main text. Although iodide clusters with two nitrogen atoms and zero nitrogen atoms both lie on odd masses, they can be separated for certain ions with the current resolving power, as demonstrated by the peak fitting results of mass spectrum at m/z 311 (Fig. S8).
The campaign-averaged mass spectra of detected ions in the both gas and particle phases are shown in Fig. 1. In general, molecules in the particle phase have larger molecular weights compared to those in the gas phase. Signals in the mass range of 150-300 Th comprise a large fraction of gasphase compounds, and concentrations in the gas phase decrease quickly with m/z higher than 250 Th. In contrast, the detected signals in the particle phase are mainly distributed within the range of 200-320 Th.
Average nighttime (22:00-06:00) and daytime (10:00-18:00) mixing ratios for various species were shown in Fig. 2. Most species have higher concentrations during the daytime, especially for relatively volatile compounds in the gas phase, despite the fact that lower boundary layer height at night should increase nighttime concentration, as many primary gases behaved, e.g., CO ( Fig. S1) . The higher concentrations during the daytime for most species detected by FIGAERO-I-CIMS suggest the dominant role of photochemical-induced oxidation in forming these oxidized compounds. In addition to typical nocturnal species including nitryl chloride (ClNO 2 I − ), chlorine nitrate (ClONO 2 I − ) and dinitrogen pentoxide (N 2 O 5 I − ), higher concentrations for the ions of C 6 H 10 O 5 I − and C 6 H 12 O 5 I − were also observed, which will be discussed in the next section. A large number of particulate N-containing organic compounds increase at night as well, as shown by mass defect diagrams of C x H y O z and C x H y N 1,2 O z color-coded by the night-to-day ratios (Fig. S9).
Based on the mass spectra shown in Fig. 1, we identify a number of ions associated with high concentrations in both the gas phase and the particle phase. In the following Sect. 3.2-3.7, we will perform interpretation of the mass spectra by analyzing variability and correlation of these important ions, including monosaccharide-derived compounds (with brown tags in Fig. 1), oxygenated aromatics (with purple tags), organic acids (with pink tags), oxidation products of biogenic volatile organic compounds (BVOCs, with green tags), sulfur-containing compounds and inorganics (with blue tags). After going through detailed analysis at the species level, Sect. 3.8 will provide an overall picture about bulk chemical characteristics of detected organic compounds in terms of the distributions of average carbon oxidation states, carbon number and oxygen number. Lastly, Sect. 3.9 will compare our measurement of organic aerosol (OA) with AMS data.

Monosaccharide-derived compounds
C 6 H 10 O 5 and C 6 H 12 O 5 are highly correlated with each other in aerosol (r = 0.92), and they are two of the few C x H y O z compounds with higher concentrations at night. Previous work attributed them to monosaccharide-derived compounds emitted from biomass burning (Bhattarai et al., 2019;Qi et al., 2019;Reyes-Villegas et al., 2018;Simoneit et al., 1999).
In this campaign, C 6 H 10 O 5 was detected mostly in the particle phase (the fraction in the particle phase F p = 0.81 ± 0.09) with an average concentration of 0.073 ± 0.076 µg/m 3 . Its diurnal profile started increasing at dusk, reaching a peak at about midnight, and then fell off, as shown in Fig. 3. The mass fraction of C 6 H 10 O 5 in OA had a similar diurnal profile, and the ratios of C 6 H 10 O 5 to CO increased at night (from 0.17 ± 0.02 to 0.5 ± 0.03 µg m −3 /ppm, Fig. 3c), both suggesting enhanced emissions of this compound were related with combustion activities in the evening, e.g., residential biofuel burning for cooking as reported by some previous measurements in China Zhang et al., 2015). Furthermore, the time variations of particle-phase C 6 H 10 O 5 I − were very similar to those of the m/z 60 fragment in AMS mass spectra (Fig. 3a), which is an identified tracer of biomass burning OA produced from the thermal decomposition of levoglucosan and similar compounds on the vaporizer of AMS (Brege et al., 2018;Schneider et al., 2006). Therefore, C 6 H 10 O 5 was probably levoglucosan and its isomers (mannosan and galactosan), and C 6 H 12 O 5 was probably also a monosaccharide compound that had sources similar to C 6 H 10 O 5 .

Oxygenated aromatic compounds
Combustion activities emit a great deal of compounds besides saccharides that the I-CIMS instrument can detect, including nitro-aromatics and guaiacol derivatives (Gaston et al., 2016;Kong et al., 2021). Nitro-benzenediols (C 6 H 5 NO 4 I − ) as well as the highly correlated homologue methyl nitro-benzenediols (C 7 H 7 NO 4 I − ) (r = 0.88 in the particle phase) exhibited double peaks in their diurnal profiles (Fig. 4). One was in the evening, similar to levoglucosan (C 6 H 10 O 5 ). The other peak was at noon. The scatterplot of C 6 H 5 NO 4 as the function of C 6 H 10 O 5 exhibits two different slopes (Fig. 5): the lower slope at night (0.088 ± 0.005) indicates the contribution of biomass burning, while the higher slope during the daytime (0.26 ± 0.02) suggests there were other important sources for nitro-aromatics, potentially secondary formation from photooxidation of aromatics (Jenkin et al., 2003). Guaiacol derivatives may have similar sources with nitro-aromatics, as implied by the resemblance of the scatterplots of these two chemical classes versus levoglucosan (cf., Figs. S10 and 5).
Nitrophenols (C 6 H 5 NO 3 I − ), methyl nitrophenols (C 7 H 7 NO 3 I − ) and dinitrophenols (C 6 H 4 N 2 O 5 I − ) were the  Wang et al. (2016) most significant components of nitro-aromatics in the gas phase. Despite the fact that nitrated phenols could be formed by photochemical oxidation of their aromatic hydrocarbon precursors (Z. Yuan et al., 2016), none of them peaked in the daytime, consistent with a previous proposal on photolysis as the dominant loss pathway for these compounds Yuan et al., 2016). C 6 H 5 NO 3 and C 7 H 7 NO 3 peaked in the evening, suggesting important contributions of NO 3 -induced reactions and/or primary emissions. The peak time of C 6 H 4 N 2 O 5 was later than that of C 6 H 5 NO 3 , in agreement with dinitrophenols as  (Table 1) in the total ion signals for I-adduct ions measured in the gas phase (b) and the particle phase (c), respectively. the oxidation products from nitrophenols (Harrison et al., 2005). We also detected non-N-containing compounds that were identified as oxidation products of aromatics in the literature, including C 7 H 6 O 4 I − , C 7 H 8 O 4 I − and C 7 H 8 O 5 I − (Mehra et al., 2020;Schwantes et al., 2017). C 7 H 6 O 4 and C 7 H 8 O 4 correlated well with each other (r = 0.72 in the gas phase and 0.91 in the particle phase). High concentrations of C 7 H 6 O 4 and C 7 H 8 O 4 were mainly observed during the periods with lower NO x concentration, which was a contrast to the variations of nitrophenols (Fig. S10). In addition, the concentration ratios of C 7 H 8 O 4 I − and C 7 H 7 NO 3 I − are lower for higher NO x concentration (Fig. 5), consistent with the literature stating that the formation of C 7 H 6 O 4 and C 7 H 8 O 4 is suppressed at high NO x concentration . C 7 H 8 O 5 is reported as the ring-retaining oxidation product of C 7 H 8 O 4 , which is a typical oxidation product of toluene and cresol , as well as the ring-scission products of aromatic hydrocarbons with more carbon atoms, e.g., trimethyl benzenes (Mehra et al., 2020). Given that C 7 H 8 O 5 closely followed with C 7 H 8 O 4 (r = 0.93 in particles), toluene oxidation was probably the main contributor to this compound.

Organic acids and related compounds
Organic acids were one of the most abundant species classes detected by I-CIMS (Fig. 1). Low-molecular-weight organic acids (e.g., formic, acetic, glycolic and pyruvic acid) constituted a significant fraction of signals in the mass spectra of the gas phase. As shown in Fig. 6 (and also Fig. S11), they had very similar temporal trends with diurnal maxima in the afternoon, indicating photochemical oxidation played a dominant role in their formation (de Gouw et al., 2018;Yuan et al., 2015).
In contrast to monocarboxylic acids, dicarboxylic acids partitioned mostly to the particle phase. As the dominant di- Figure 2. The ratios of concentrations at night (22:00-06:00) to concentrations during the day (10:00-18:00) for ions ranging from 150 to 400 Th in the gas phase (a) and the particle phase (b). The range of the y axis is set between 0 and 2 for clarity, although the ratios of some compounds are larger than 2. The numbers in boxes indicate the night / day ratios of tagged ions that exceed the y-axis ranges.
carboxylic acids in aerosol (Kawamura and Bikkina, 2016;Mellouki et al., 2015), 94 ± 5 % and 74 ± 13 % (mean ± 1 standard deviation of F p ) of C 2 H 2 O 4 and C 3 H 4 O 4 , attributed to oxalic and malonic acid, were found in the particle phase, respectively. The concentrations of C 4 H 6 O 4 were significantly lower compared to the C 2 and C 3 homologous series, but C 5 H 8 O 4 and C 6 H 10 O 4 had unexpected high abundance (Fig. 7). Additionally, C 5 H 8 O 4 and C 6 H 10 O 4 had considerable fractions in the gas phase (45 ± 13 % and 43± 11 %), which are significantly higher than their C 2 -C 3 homologous series. These two compounds were correlated well with each other in temporal variations (r = 0.97 and 0.91 in the gas phase and the particle phase, respectively), and their diurnal variations were different from those of oxalic and malonic acid (Fig. 6). Therefore, dicarboxylic acids may not be the dominant contributing species for the two compounds. C 5 H 8 O 4 and C 6 H 10 O 4 have been observed from previous studies on isoprene oxidation (Berndt et al., 2018(Berndt et al., , 2019, attributing them to epoxy hydroperoxyl carbonyl and accretion product, respectively. However, the relative contributions from these possibilities remain unclear. In  (Brege et al., 2018;Kawamura et al., 1996;Kawamura and Bikkina, 2016). Our measurements showed that temporal trends of C 4 H 4 O 4 and C 5 H 6 O 4 agreed well with those of aromatic hydrocarbons (Fig. S11b), and thus oxidation of aromatics could be an important contributor to C n H 2n−4 O 4 in the urban air.

Oxidation products of biogenic VOCs
In addition to high anthropogenic emissions of aromatics, terrestrial vegetation nearby also released significant amounts of biogenic VOCs (BVOCs) . During the campaign, the concentrations of isoprene at noon were between 0.1 and 1.5 ppb, whereas the range of daily maxima of monoterpenes was 0.05-2.5 ppb. Hence, a number of oxidation products of BVOCs were detected (Figs. 8 and S12).
The ion C 4 H 7 NO 5 I − was the most abundant N-containing C 4 organic compounds that was detected in the gas phase. Its daily maxima occurred in the afternoon and correlated moderately with methyl vinyl ketone (MVK) + methacrolein (MACR) measured by PTR-ToF-MS (Fig. 8b, r = 0.58). We consequently attributed C 4 H 7 NO 5 to MVK nitrates and MACR nitrates, which were reported as the second generation of organic nitrates formed from the oxidation of isoprene hydroxynitrates by OH in the presence of NO x (Fisher et al., 2016;Paulot et al., 2009). Strong correlations were observed between C 5 H 9 NO 4 I − , C 5 H 9 NO 5 I − and C 4 H 7 NO 5 I − (r = 0.93 and 0.80, respectively), which was in accordance with their similar formation pathways (Jacobs et al., 2014;Wennberg et al., 2018;Xiong et al., 2015). Hence, we expected these three compounds to be common oxidation products of isoprene in the polluted atmosphere. While in aerosol, 2-methylglyceric acid (C 4 H 8 O 4 ) is a commonly reported oxidation product of isoprene formed in high-NO x conditions (Surratt et al., 2010). We observed the corresponding ion C 4 H 8 O 4 I − contributing to OA especially in dry conditions with strong sunlight (Fig. S13). This evidence indicates that isoprene oxidation may contribute to C 4 H 8 O 4 , but potential contribution from other sources cannot be ruled out in urban areas.
In terms of monoterpenes, a reasonable correlation (Fig. S14a, r = 0.63) was found between the ions C 10 H 16 O 3 I − and C 10 H 16 O 2 H + measured by PTR-ToF-MS. C 10 H 16 O 2 H + was attributed to pinonaldehyde formed from the oxidation of monoterpenes (Glasius et al., 2000;Larsen et al., 2001;Mutzel et al., 2016). Therefore, we tentatively attribute C 10 H 16 O 3 I − to pinonic acid and its oxocarboxylic acid isomers, which are formed via the oxidation of pinonaldehyde (Fang et al., 2017). C 8 H 13 NO 6 also exhibited enhanced gas-phase formation during the day as pinonic acid did. The correlation coefficient of the two compounds was 0.71. In contrast to other monoterpene nitrates, particlephase C 8 H 11 NO 7 and C 10 H 15 NO 6 peaked at night and decreased during the daytime (Fig. S12), indicative of the role of NO 3 in producing organic nitrates as reported in the literature (Faxon et al., 2018). However, C 10 H 15 NO 6 in the gas phase showed a distinct diurnal profile with peak before the noon. Two possible types of compounds were pro- posed for C 10 H 15 NO 6 in previous studies: peroxyacetyl nitrate from pinonaldehyde (Faxon et al., 2018;Nah et al., 2016;Schwantes et al., 2020) or organic nitrates (Bean and Hildebrandt Ruiz, 2016;Boyd et al., 2015). Given the distinct diurnal profiles of C 10 H 15 NO 6 I − in the gas and particle phases and the fact that peroxyacetyl nitrate is supposed to dissociate during the FIGAERO heating (Slusher et al., 2004), we speculate that both compounds contributed to this ion. As shown in Fig. S15, C 8 H 12 O 4 and C 9 H 14 O 4 existed mostly in the particle phase (F p = 0.63 ± 0.11 and 0.67 ± 0.10, respectively). We interpret them as products of monoterpenes via photochemical processes, consistent with the interpretations presented in previous work (Mohr et al., 2013;Mutzel et al., 2015).

S-containing compounds
Organosulfates are concerned as important components of SOA (Hallquist et al., 2009;Surratt et al., 2007), and they can be detected by iodide anion via proton abstraction (Le Breton et al., 2018b;Lee et al., 2014). We detected the ion C 2 H 3 SO − 6 with a diurnal peak in the afternoon (Fig. 9). This ion was attributed to glycolic acid sulfate, as suggested by previous work (Galloway et al., 2009;Liao et al., 2015).
Abundant SO 3 I − was detected in particles, and it correlated well with the ion C 2 H 3 SO − 6 (Fig. 9b) and sulfates measured by AMS (Fig. S16). Previous work observed the sulfite ion radical ( q SO − 3 ) during the ionization of organosulfates in a liquid chromatography-electrospray ionization-tandem mass spectrometer (Huang et al., 2018). As a result, the ion SO 3 I − from FIGAERO-I-CIMS might be a potential indicator for the total organosulfates. However, more future work is needed for evaluating this possibility.
Other sulfate-related ions during gas-phase modes were also detected, including HSO − 4 (sulfuric acid) and CH 3 SO − 3 (methanesulfonic acid), which were enhanced in the gas phase during the daytime, in agreement with the notions of photochemistry-induced gas-phase oxidation (Brandt and van Eldik, 1995). However, these data were not available for quantification given that these low-volatility species would condense on our long gas sampling inlet. It should be noted that measuring sulfuric acid in the gas phase is difficult and generally requires a "wall-less" source design (Eisele and Tanner, 1993).

Inorganic compounds
There is a growing interest in N 2 O 5 and its product nitryl chloride (ClNO 2 ) because ClNO 2 is found to serve as a nocturnal reservoir of Cl radical and reactive nitrogen, and hence enhance the ozone formation the next day (Osthoff et al., 2008;Wang et al., 2016). Time series of N 2 O 5 and ClNO 2 exhibited two patterns. During most of the nights, N 2 O 5 started to increase quickly at sunset and lasted for only 2-3 h, and ClNO 2 increased in the meantime and ultimately reached its maximum at night, indicative of local formation of ClNO 2 . However, sometimes a high level of N 2 O 5 did not lead to an increase in ClNO 2 (tinted background in Fig. 10a), probably due to the lack of chloride salts on the aerosol. Other nocturnal species including ClONO 2 and Cl 2 were highly correlated with ClNO 2 as expected (r = 0.92 and 0.83, respectively), suggesting they had common formation mechanisms . HNO 3 I − was observed as one of the most abundant species in the mass spectra of FIGAERO-I-CIMS both in the gas phase and the particle phase. In the gas phase, the ion HNO 3 I − from I-CIMS has been used to quantify nitric acid (Lee et al., 2018a). The concentrations of gas-phase nitric acid peaked in the afternoon, suggesting photochemistry in the daytime as the dominant source for gas-phase nitric acid.
Previous studies suggested that HNO 3 I − from particlephase measurement by FIGAERO-I-CIMS can be indica- Figure 5. (a) Correlation between particle-phase C 6 H 5 NO 4 I − and C 6 H 10 O 5 I − . The data points are color-coded using the time of the day. Solid and dash lines represent the slopes during the nighttime and daytime, respectively. (b) Relative concentration of C 7 H 8 O 4 I − and C 7 H 7 NO 3 I − in the gas phase as a function of NO x concentration. The black line is the fitted curve using a double exponential function.
tive of nitrate in the particle phase . Here, the concentrations of HNO 3 I − in the particle phase were compared with particulate nitrate measured by AMS (Fig. 11c). A strong correlation was observed (r = 0.93), but the concentrations measured by FIGAERO-I-CIMS were higher (slope = 1.6), especially for higher concentrations of organic nitrates. Using a threshold of 1 µg/m 3 for organic nitrates, the slopes and correlations were higher for the data points with particulate organic nitrates larger than 1 µg/m 3 (slope = 1.8, r = 0.94) than those less than 1 µg/m 3 (slope = 1.1, r = 0.90). In short, our measurements suggest that HNO 3 I − in the particle phase from FIGAERO-I-CIMS is formed from thermal decomposition of both inorganic nitrate (e.g., NH 4 NO 3 ) and organic nitrates.

Bulk chemical properties of detected organic compounds
The above discussions on individual chemical groups provide insights into the identification of the mass spectra from FIGAERO-I-CIMS, along with sources and chemistry of Figure 6. Diurnal variations of organic acids in the gas phase (red) and particle phase (blue). The shaded area indicates 1 standard deviation. oxygenated organic compounds in the urban atmosphere. In this section and the following one, we will provide a bulk analysis of the detected organic compounds. Organic compounds detected by FIGAERO-I-CIMS were comprehensively characterized with detailed elementary composition in OS C − n C space (Fig. 12), which depicts the average oxidation states of carbon for closed-shell C x H y O z and C x H y N 1,2 O z compounds as a function of carbon number. The details in the calculation of OS C can be found in Sect. S2 in the Supplement. S-containing compounds were  omitted given their negligible variety and concentration compared to C x H y O z and C x H y N 1,2 O z . The average OS C in the particle phase was higher than that in the gas phase at the same carbon number, especially for carbon numbers between 2 and 10. This agrees with our expectation that more oxidized compounds would partition more strongly in aerosol, as indicated by larger fractions in particles (Fp) for higher OS C .
In addition, the average OS C generally increased for lower carbon number, as a result of functionalization and fragmentation during VOCs aging. However, there was a notable exception in C 5 which had a significantly reduced OS C , probably as the result of emissions of isoprene. The analysis of  The solid and dash lines show the fitted results for the dataset of pON less than 1 µg/m 3 and more than 1 µg/m 3 , respectively. The concentration of gaseous HNO 3 I − shown here only included the last 5 min of every gas-phase working mode, as high level of HNO 3 came out of aerosol, which then passed through the CIMS in a short time during particle analysis, and a substantial amount would subsequently accumulate on the inner surfaces, leading to a persistent carried over signal that was long enough to disturb the next gas measurement cycle (Palm et al., 2019). the OS C − n C space indicates that the large number of organic compounds measured by FIGAERO-I-CIMS are useful to characterize the evolution of organic compounds in the atmosphere. The distributions of carbon and oxygen numbers of organic compounds were also investigated, as shown in Fig. 13. Most abundant organic compounds measured by FIGAERO-I-CIMS were C 2 -C 3 compounds, which accounted for about 66 % in the gas phase and 56 % in the particle phase. It is unexpected that C 2 -C 3 compounds made up such a significant portion in the particle phase, indicating a non-negligible role of thermal decomposition from low volatility compounds such as accretion products or extremely low volatility organic compounds which were reported from FIGAERO measurements on SOA Stark et al., 2017). Organic compounds with carbon numbers over five constituted only 3 % in the gas phase, while they accounted for 30 % in the particle phase. The oxygen numbers of the majority of gaseous organic compounds were no more than 3. Organic compounds containing two to three oxygen atoms had the largest contribution in both the gas phase (96 %) and the particle phase (56 %). C x H y N 1,2 O z accounted for less than 10 % of the total oxygenated or-  ganic compounds. In the gas phase, compounds with five to six oxygen atoms accounted for 51 % of C x H y N 1,2 O z , indicative of the high levels of organic nitrates in the urban atmosphere. Nitrophenols also contributed significantly to C x H y N 1,2 O z compounds, as they accounted for 74 % of C x H y N 1,2 O z containing three oxygen atoms, which in turn contributed to 22 % of C x H y N 1,2 O z . In contrast, in the particle phase, the oxygen numbers of C x H y N 1,2 O z distributed relatively evenly, as the fractions of compounds with three to eight oxygen atoms were similar (between 12 % and 19 %). Compared to measurements in a forest in the southeastern United States (cf., Table S1 from Lee et al., 2016), the fractions of N-containing organic compounds with less than five oxygen atoms were significantly larger in our measurements as a result of higher concentrations of nitro-aromatics.
We further determined the fractions of N-containing organic compounds in total organic compounds as a function of m/z. It is clear that the observed fractions of N-containing organic compounds were higher for elevated m/z (Fig. 14) and N-containing ions commonly dominated at even nominal masses (Fig. S17). The gas-phase CHON ions within the m/z range of 250-350 Th accounted for about half of the organic compounds in this range. The fractions of CHON ions in the particle phase were somewhat smaller than those in the gas phase within the above m/z range but were comparable for higher m/z. A possible explanation for this is that functional groups of nitrate and nitro reduce less in vapor pressure for organic compounds than functional groups of carboxylic or oxygen-equivalent hydroxyl do (Capouet and Müller, 2006;Nannoolal et al., 2008;Pankow and Asher, 2008). Consequently, CHON compounds are generally more volatile than CHO compounds with similar molecular weights.
In the end, we determined the total concentration of Ncontaining organic compounds in the particle phase measured by FIGAERO-I-CIMS and compared it with the particulate organic nitrates derived from AMS (Fig. 15). Good agreement was achieved when the concentration of inorganic nitrate was relatively low, e.g., below 8 µg/m 3 . However, the discrepancies increased when inorganic nitrate concentration increased, which can affect the determination of organic nitrate from AMS. This encouraging result indicates that FIGAERO-I-CIMS is able to capture the variability of organic nitrates in the urban atmosphere, which can be helpful in understanding the sources and formation mechanism of these compounds.

Organic aerosol measurements
The total concentration of organic compounds in the particle phase measured by FIGAERO-I-CIMS was determined and compared with measurements of OA by AMS. The total organic compounds measured by FIGAERO-I-CIMS explained 24 ± 0.8 % (fitted slope ± 1 standard deviation) of the total OA on average (Fig. 16a), which is lower than the average fractions (∼ 50 %) reported previously in boreal and temperate forests Stark et al., 2017). The lower fractions determined here might be the result of larger contributions to OA from primary emissions in the urban air, which are composed of a large number of compounds with little signals in I-CIMS (Zhao et al., 2016). As shown in Fig. 16a, organic compounds measured by FIGAERO-I-CIMS accounted for higher fractions in OA concentrations by AMS for more aged OA, consistent with the fact that I-CIMS is more sensitive to oxygenated organic compounds with multiple functional groups (Lee et al., 2014;Lopez-Hilfiker et al., 2016). Furthermore, we expect this fraction to change with the relative contributions of primary emissions and secondary formation for organic compounds in the atmosphere. Similar trends were found in Le Breton et al. (2019), in which an acetate source was used. Acetate ions have been reported to selectively ionize highly oxygenated organic compounds as an iodide source does (Aljawhary et al., 2013).
Comparison of the van Krevelen diagrams between FIGAERO-I-CIMS and AMS also provides useful insights on the measurement of organic compounds in OA. The van Krevelen diagram has been used as a tool for analyzing functional groups and OA aging by plotting H / C ratios versus O / C ratios (Heald et al., 2010;Lambe et al., 2012). As shown in Fig. 16b, the data points of the bulk OA from FIGAERO-I-CIMS followed the same trend as the data points from AMS. However, the bulk OA measured by FIGAERO-I-CIMS only occupied a much smaller region with O / C ratios between 0.7 and 1.0. We further plotted all of the organic compounds in the H / C versus O / C space color-coded with their campaign-average concentrations (Fig. S18a). We observed that most data points from FIGAERO-I-CIMS distributed across the zone between the slope of 0 and −1.0. These observations provide additional evidence that FIGAERO-I-CIMS may only measure the more oxidized organic compounds in OA.
The correlation coefficients between the particle-phase concentrations at unit masses by FIGAERO-I-CIMS and OA mass concentration by AMS were calculated (Fig. S18b). The correlation coefficients were small for ions below m/z 200, as these ions contributed little to organic aerosol. Moderate and strong correlations (r > 0.7) were observed for the ions of m/z between 200 and 400 Th, implying that organic compounds with molecular weights of 100-300 g/mol may account for significant fractions in organic aerosol. The possible reason for the lower correlations for heavier compounds (m/z > 400 Th) with OA mass loadings is that these compounds might be related to specific sources or certain chemical processes, which might not contribute in large fractions to the total OA concentration.

Summary
We deployed a FIGAERO-I-CIMS instrument to measure oxygenated organic compounds in both the gas phase and the particle phase at a representative urban site in China. The experimental design and instrumentation setup were described in detail, which goes above and beyond typical studies, including (1) performing sensitivity calibrations in the laboratory using multiple methods for multiple species; (2) performing voltage scanning for unknown compounds detected in the ambient air; and (3) performing humidity calibrations for multiple types of species, which we have not seen anyone do after Lee et al. (2014).
From the mass spectra, a number of important compounds in the urban atmosphere were identified. We detected high concentrations of several monosaccharide species (e.g., levoglucosan). They were potentially emitted from biomass burning, which also contributed to the enhancement of many nitro-aromatic species. Photochemistry was also found to be a strong source of nitro-aromatics. Low-molecular-weight organic acids were mainly observed in the gas phase, and observations supported daytime photochemistry as the dominant source. Different diurnal profiles for various BVOCderived organic nitrates were observed, reflecting their different formation pathways related to NO x chemistry (i.e., daytime photooxidation, nocturnal NO 3 reactions). Local formation of nitryl chloride was observed, highlighting the potential importance of nighttime chemistry in the urban region. Our measurements show that oxygenated organic compounds dominated the majority of detected species by FIGAERO-I-CIMS, in which CHO and CHON compounds both accounted for significant fractions. Nitrogen-containing organic compounds occupied a significant fraction of the total signals in both the gas and particle phases, with elevated fractions at higher molecular weights. The most abundant organic compounds were formic acid and multifunctional organic compounds containing three to five oxygen atoms. Organic compounds containing two to three carbon atoms accounted for over half of the total organic compounds in both the gas phase and the particle phase measured by FIGAERO-I-CIMS. During the campaign, the FIGAERO-I-CIMS measurements explained 24 ± 0.8 % of OA mass measured by AMS, but the fractions were higher for more aged organic aerosol. This evidence, along with the analysis of the van Krevelen plot, indicates that FIGAERO-I-CIMS is measuring the more oxidized fraction of OA in the urban air.
Our observations suggest that oxygenated organic compounds in urban environments are complicated in both sources and chemistry. Oxygenated organic compounds can be both emitted from various emission sources (e.g., vehicular emissions and biomass burning) and also secondarily produced in the atmosphere. The chemistry in forming and removing these oxygenated organic compounds can be associated with daytime and nocturnal reactions initiated by both anthropogenic and biogenic precursors with strong influences from NO x chemistry. This work demonstrates that the rich information in both gas and particle phases provided by FIGAERO-I-CIMS can greatly promote the understanding of emission and chemistry of organic carbon in the atmosphere of urban regions.
Data availability. The more detailed data can be provided by contacting the corresponding authors.
Author contributions. BY and MS designed the research. CSY, YL, ZLW, TGL, WWH, WC, CHW, CMW, SH, JPQ, BLW, CW, WS, XMW, ZYZ and XMW contributed to data collection. CSY performed the data analysis with contributions from WWH and WC. CSY and BY prepared the manuscript with contributions from JEK and other authors. All the authors reviewed the manuscript.
Competing interests. The authors declare that they have no conflict of interest.
Acknowledgements. The authors would like to thank Felipe Lopez-Hilfiker for helping to tune the CIMS instrument. The authors appreciate valuable comments from four anonymous reviewers. Review statement. This paper was edited by Arthur Chan and reviewed by four anonymous referees.