Probing key organic substances driving new particle growth initiated by 1 iodine nucleation in coastal atmosphere

ABSTRACT


Introduction
Atmospheric new particle formation (NPF) contributes over half of global cloud condensation nuclei (CCN) (Merikanto et al., 2009) and thus influences cloud properties and Earth's radiation budget (Metzger et al., 2010). By deploying a highresolution chemical ionization mass spectrometer, recent laboratory and field studies have identified a group of highly oxygenated multifunctional organic molecules (HOMs) with high O/C ratios and low volatility from the reactions of volatile organic compounds (VOCs) such as monoterpenes (Ehn et al., 2014), sesquiterpenes  and alkenes (Mentel et al., 2015) with a hydroxyl radical (OH), ozone (O 3 ), nitrate radicals (NO 3 ) and chlorine atoms . These HOMs play an important role in particle nucleation and growth of continental NPF, as well as in the formation of secondary organic aerosols.
Unlike the deep understanding of continental HOMs, little is known about the role of organic compounds in the NPF in coastal or open-ocean atmosphere. The current state of knowledge is that the photolysis of molecular iodine (I 2 ) or iodomethane is the source of iodine oxides or oxoacids, the self-clustering of which could initiate NPF events with a particle number concentration sometimes exceeding 10 6 cm −3 Saiz-Lopez and Plane, 2004;Burkholder et al., 2004;Mahajan et al., 2010Mahajan et al., , 2012Sipilä et al., 2016;Stevanović et al., 2019;Kumar et al., 2018). But it is unknown if other species are needed to drive the growth of iodine clusters to CCN sizes in coastal or open-ocean atmosphere . Iodineinduced NPF (I-NPF) events were mostly reported in European coastlines Mahajan et al., 2010) and polar regions (Allan et al., 2015;Roscoe et al., 2015;Dall'Osto et al., 2018). In 2019 we provided evidence of I-NPF in the southeast coastline of China, based on particle number size distribution and iodine measurements . The focus of that paper  is, however, the speciation of organic iodine compounds in sizesegregated aerosol samples. Moreover, the use of a relatively low-resolution time-of-flight (TOF) mass analyzer and in vitro signal amplification approach in that paper did not allow the detection of the majority of non-aromatic organic iodine compounds. Organic iodine remains the most significant unknown in aerosol iodine chemistry at present .
Fourier transform ion cyclotron resonance (FT-ICR) coupled with soft ionization techniques such as electrospray ionization (ESI) and ambient pressure chemical ionization (APCI) allows characterization of complex organic mixtures at the molecular level due to its ultrahigh resolution and mass accuracy (Pratt and Prather, 2012). This technique has been used to examine molecular composition of organic aerosols (Schum et al., 2018;Zuth et al., 2018;Daellenbach et al., 2019;Xie et al., 2020) and cloud water (Zhao et al., 2013;Bianco et al., 2018). Studies investigating coastal organic aerosols have been rare. Virtually no study reported the characterization of organic compounds driving the growth of coastal or open-ocean new particles.
In this study, comprehensive chemical composition analyses were conducted on the size-segregated aerosol samples down to 10 nm, collected by 13-stage nano-MOUDI (nanomicro-orifice uniform deposit impactor) during the intense I-NPF days at a coastal site in China. Relative abundances of HSO − 4 , total iodine and total organic carbon (TOC) in 10-56 nm particles were compared between the I-NPF days and conventional continental NPF (C-NPF) days. In particular, using ultrahigh-resolution FT-ICR-MS, we conducted a nontarget analysis of particle-phase organic compounds to explore their molecular identity, formation mechanism and the role in new particle growth in the coastal atmosphere.

Sampling collection
The sampling site (29 • 29 N, 121 • 46 E) is near a small fishing village without permanent residents on the coastline of the East China Sea. It can be seen from the aerial photo ( Fig. S1a in the Supplement) that from the east to the west are the sea, intertidal zone, small paddy fields and mountain. The sampling site is about 40 and 200 m away at high tide and low tide, respectively. The classification of I-NPF event, C-NPF event or non-NPF was based on particle number size distributions (PNSDs) between 2 and 750 nm monitored from January to May 2018 by a scanning mobility particle spectrometer (SMPS; TSI DMA3081 and CPC3775; scanning range: 40-750 nm) and a neutral cluster air ion spectrometer (NAIS; scanning range: 2-42 nm). A nano-MOUDI sampling scheme was implemented according to the PNSD measurement. One set of nano-MOUDI samples was collected during the C-NPF days from 11 to 13 February; a second set was collected during overcast non-NPF days from 16 to 18 April; a third set was collected during the I-NPF days from 9 to 11 May. The PNSD during the three periods are shown in Fig. S2. Each set of nano-MOUDI samples was collected continuously for 72 h, during which NPF occurred on a daily basis, so that particle chemical composition of different event types can be obtained from offline analyses. Aluminum foil filters were used as sampling substrate to avoid the adsorption of gaseous species. For each set of nano-MOUDI samples, two nano-MOUDIs were placed side by side to collect 10-100 nm particles (on stages 10-13; other stages were silicon greased) and 100 nm-18 µm particles (on stages 1-9) separately, in order to reduce potential positive particlebounce artifacts. Three additional sets of blank samples were collected by placing a high-efficiency particulate air (HEPA) filter at the gas inlet of the nano-MOUDI. Detailed information on aerosol sample collection could be found in .

Sample preparation and analysis
Half of each filter was transferred into a 20 mL tapered plastic centrifuge tube, and 10 mL mixed solvent (1 : 1 v/v water and methanol; LCMS grade) is added. The mixture was sonicated for 40 min and filtered by a 0.2 µm PTFE membrane syringe filter. The filtrate was evaporated to almost dryness in a rotary evaporator below 40 • C and subsequently redissolved in 0.5 mL water. After being centrifuged for 30 min at 12 000 rpm, the supernatant was collected for total iodine (I) analysis by Agilent 7500a ICP-MS (Agilent Technologies, Santa Clara, CA, USA) and HSO − 4 analysis by UPLC-ESI-Q-TOF-MS. The measurements of HSO − 4 and total I were elaborated in our previous article Yu et al. (2019). Another half of each filter was extracted in the same way but used for TOC analysis by a TOC analyzer (Model TOC-5000A, Shimadzu, Japan) and nontarget MS analysis of organic compounds (OCs) by ESI-FT-ICR-MS (SolariX XR 9.4T instrument, Bruker Daltonics, Coventry, UK). Samples were infused by a syringe pump and analyzed in both positive (ESI+) and negative (ESI−) modes. ESI-FT-ICR MS operation conditions are included in the Supplement. Field blank sample extracts were analyzed following the same procedure.

FT-ICR MS data processing
A resolving power (m/ m 50 % ) 550 000 at m/z 300 of our FT-ICR-MS allows the determination of possible formulas for singly charged molecular ions. Only m/z values between 150-1000 that satisfy the signal-to-noise (S/N) ratio > 10 were considered. For each m/z value, several scientific rules were applied to calculate a reasonable elemental formula of natural organic molecule: the general formula is C 1-50 H 1-100 O 0-50 N 0-10 I 0-3 in the ESI+ mode; elemental ratios H/C, O/C and N/C are limited to 0.3-3, 0-3 and 0-1.3, respectively. The general formula is C 1-50 H 1-100 O 1-50 N 0-5 S 0-2 I 0-3 in the ESI− mode; elemental ratios H/C, O/C, N/C and S/C are limited to 0.3-3, 0-3, 0-0.5 and 0-0.2, respectively. Mass error must be smaller than 0.5 ppm. Formulas containing C, H, O, N, S and I isotopologues were removed from the formula lists. A formula with m/z > 500 was not reported if it did not belong to any CH 2 homologous series. For the formula C c H h O o N n S s I x , doublebond equivalents (DBEs) defined as DBE = 2c+2−h+n−x 2 were required to be nonnegative. Formula calculation was done following the same procedure for the three field blank samples. All formulas found in the field blank samples, regardless of peak intensity, were excluded from the formula lists of real samples. Aromaticity index (AI) is calculated If DBE AI ≤ 0 or C AI ≤ 0, then AI = 0. A threshold value of AI ≥ 0.5 provides an unambiguous minimum criterion for the presence of aromatic structure in a molecule (Yassine et al., 2014).

Organic compounds dominate the growth of new particles initiated by iodine nucleation
Although our offline technique did not allow us to probe nucleating cluster composition at ∼ 1.7 nm, four facts from our observation support that the NPF events from 9 to 11 May were initiated by iodine nucleation. First, strong I-NPF events were observed almost every sunny day in April and May, which was the growth and farming season of seaweed. HYSPLIT back-trajectories analysis (Draxler and Rolph, 2010) shows that air masses moved from the East China Sea to the sampling site during the I-NPF days from 8 to 10 May 2018 (Fig. S1b). Sea breeze was also expected to flow from the sea to the site in the daytime when the I-NPF events occurred. Second, the evolution of PNSD from 9 to 11 May was not like banana-shaped C-NPF observed on the winter days but rather was markedly similar to prior reports of iodine nucleation at European coastal sites (Mäkelä et al., 2002;Sipilä et al., 2016). Third, the production of 2-7 nm particles (N 2-7 ) during the C-NPFs followed a nearly identical variation with solar radiation (Fig. S2c), which is an indication that the C-NPFs was initiated by OH and H 2 SO 4 production dictated by solar radiation. However, this was not observed during the I-NPF events; instead, N 2-7 was anticorrelated with tidal height in the daytime (Fig. S2a). Fourth, probably the most important, mean total I in 10-56 nm particles during the I-NPF days (13.5 ng m −3 , Table 1) was 67 and 36 times higher than those during the C-NPF days (0.2 ng m −3 ) and non-event days (0.37 ng m −3 ), respectively. In the same size range, mean HSO − 4 concentration (0.2 µg m −3 ) during the I-NPF days was lower than that during the C-NPF days (0.5 µg m −3 ).
By assuming nanometer-sized particles are spherical with a density of 1.5 g cm −3 , we estimate from the PNSD data that aerosol mass in the 10-56 nm size range was enhanced by 3.0 and 1.3 µg m −3 at most by the selected I-NPF and C-NPF events ( Fig. S2b and d). The fraction of organic mass (OM) in the aerosol mass can be further calculated as (1.5×m TOC )/(m Total I +m − HSO 4 +1.5×m TOC )×100 % by assuming an OM/TOC ratio of 1.5. The result shows that mass fractions of OM are 95 %, 87 % and 68 %, respectively, in the size bins 10-18, 18-32 and 32-56 nm during the I-NPF days. Therefore, organic compounds dominate the aerosol mass in the 10-56 nm new particles during the I-NPF days and were critical for I-NPF to contribute to CCN. Our result is qualitatively consistent with previous measurements showing that nucleation mode particles initiated by iodine were composed of a remarkable fraction of organic compounds and sulfate (Mäkelä et al., 2002;Vaattovaara et al., 2006). The main purpose of this article is to identify these organic compounds during the I-NPF days. The OC composition during the C-NPF days is beyond the scope of this article.

Elemental composition of non-iodinated OC on the I-NPF days
Nontarget analysis of OC elemental composition was performed in detail on 10-18 nm, 32-56 nm, 180-560 nm and 3.2-5.6 µm particles during the I-NPF days. Elemental formulas in the four size bins can represent OC molecular composition of nucleation mode, Aitken mode, accumulation mode and coarse mode, respectively. OC formulas were divided into two categories: non-iodinated OC and iodinated OC. There are far more non-iodinated OC formulas than iodinated OC formulas in < 1 µm particles in terms of both formula number (Table 2) and relative intensity (Fig. 1). For example, 2831 non-iodinated OC formulas account for 96.6 % of OC total intensity in 10-18 nm particles, while 137 iodinated OC formulas account for the remaining 3.4 %. It means that non-iodinated OC dominates new particle growth during the I-NPF events. In this section, we first discuss chemical characteristics of non-iodinated OC, while the speciation of iodinated OC will be discussed in Sect. 3.4. The molecular formulas of non-iodinated OC were divided into seven elemental groups, CHO − , CHO + , CHON − , CHON + , CHOS − , CHONS − and CHN + . The number distribution of seven elemental groups for the four size bins is Table 1. Concentrations of total iodine (I), HSO − 4 and total organic carbon (TOC) in three size bins between 10-56 nm during the I-NPF, C-NPF and non-NPF days. For simplicity, only the sum of three size bins is are shown for the C-NPF and non-NPF days. BDL: below detection limit.

I-NPF
C-NPF Non-NPF 10-18 nm 18-32 nm 32-56 nm 10-56 nm 10-56 nm It should be noted that some formulas were repeatedly detected in ESI+ and ESI− modes. Some formulas detected in one size bin were also detected in another size bin. This is quantitatively shown in the first four rows of Table 2. For instance, 58 %, 25 % and 4 % of the formulas detected in 10-18 nm aerosols were also detected in 32-56 nm, 180-560 nm and 3.2-5.6 µm aerosols, respectively. In other words, the particles in neighboring size bins share more similarity in organic composition. An unexpected finding is that the number of non-iodinated OC formulas detected in 3.2-5.6 µm coarse particles (n = 266) is 1 order of magnitude lower than those of other bins. Reconstructed mass spectra of the seven elemental groups in ESI− and ESI+ modes are shown in Fig. S3 for the four size bins. CHON is the most commonly assigned elemental group in both ESI+ (2121 CHON + ) and ESI− (1268 CHON − ) modes, followed by the CHO group (857 CHO + formulas and 892 CHO − formulas). S-containing formulas are 357 CHOS − and 259 CHONS − . The formula number of the least common CHN + group is only 46. In terms of relative intensity, the CHON fraction in the ESI− mode decreases from 61 % of OC in the 10-18 nm bin to 20 % in the 180-560 nm bin (Fig. 1a), while the fractions of CHO and CHOS/CHONS increase with particle size. In the ESI+ mode, the fraction of CHON decreases from 88 % in 10-18 nm bin to 70 % in 180-560 nm bin, being always the dominant elemental group of non-iodinated OC (Fig. 1b). Low-molecular-weight (LMW) amines are important stabilizers in acid-base nucleation (Kurtén et al., 2008;Jen et al., 2014;Zheng et al., 2000;Yao et al., 2016), but their molecular ions are out of the mass range of our FT-ICR-MS. The CHN + formulas that we observed contained 9-50 C atoms and 1-7 N atoms, accounting for a negligible fraction, 1.7 %, of total intensity of all ESI+ formulas in the 10-18 nm particles.
Previous elemental composition studies using FT-ICR-MS were mostly conducted on PM 2.5 or PM 10 collected from marine (Schmitt-Kopplin et al., 2012;Bao et al., 2018;Ning et al., 2019), urban (Wu et al., 2019;Jiang et al., 2016), free troposphere (Schum et al., 2018;Mazzoleni et al., 2012) and forest sites (Kourtchev et al., 2013). In general, these studies showed that the numbers of CHO compounds were comparable with or more than those of CHON compounds. Our study shows clearly that elemental composition of aerosol OC is highly size dependent. New particle growth in the size range of 10-18 nm during the I-NPF event is dominated by the CHON elemental group, followed by CHO. The focus of this article narrows in on the identity and source of the CHON and CHO formulas in 10-18 nm particles, by comparing with those in the 180-560 nm size bin.

CHO formulas
There is a total of 531 CHO − formulas and 250 CHO + formulas in 10-18 nm particles. There are 54 CHO formulas commonly found in ESI+ and ESI− modes. In terms of relative intensity, CHO − compounds are more abundant than CHO + compounds (Fig. 3b, total intensity: 4.14 × 10 9 vs. 1.24 × 10 9 ). However, this is not indicative of absolute concentration of the two groups due to different ionization efficiency between ESI− and ESI+ modes. CHO − is characterized by a series of formulas with 20, 24, 28 and 33 C atoms, four or six O atoms, and one equivalent double bond (Fig. 2b). The total intensity of the top 10 formulas accounts for 30 % of all 531 formulas. Assuming CHO − formulas contain at least one carboxylic group (-COOH), the rest of their molecules should be saturated (DBE = 0) and contain two or four O atoms.
The above feature is not seen in either CHO + formulas in the 10-18 nm bin or CHO − formulas in the 180-560 nm bin.
There are more C 21 and C 24 formulas than other C subgroups in the CHO + formulas of the 10-18 nm bin (Fig. S4d), but none of them have exceptionally high intensity. The prominent formulas in the CHO − group in 180-560 nm particles have a relatively high unsaturation degree (DBE = 3-7; Fig. 2d). The relative intensities of subgroups according to C atom number in the CHO − formulas in the 180-560 nm bin are characterized by trimodal distribution with maximum intensity around C 9 , C 13 -C 16 and C 20 (Fig. 3d). The relative intensity of O atom subgroups is mono-modally distributed around O 7 (Fig. S5).

CHON formulas
As discussed earlier, CHON is the most abundant elemental group observed in the smallest size bin 10-18 nm. There is a total of 1005 CHON − formulas (total intensity 9.96 × 10 9 ) and 1139 CHON + formulas (6.45 × 10 9 ) in 10-18 nm bin. There are 355 CHON formulas commonly found in ESI+ and ESI− modes. A close examination of Figs. 2a and 3a reveals that CHON − is characterized by a series of C 18 and C 30 formulas with low DBE values (1-4). The 87 C 18 and 26 C 30 formulas account for 37 % of total intensity of CHON − . Such a feature is not seen for CHON + formulas that are rather uniformly distributed in the DBE vs. C diagram (Fig. S4a and  c). Generally speaking, CHON − compounds should contain a nitro-(-NO 2 ) or nitrooxy-(-ONO 2 ) group and are ionizable due to the presence of -COOH or hydroxy (-OH) (Lin et al., 2012). However, the presence of an amine group in CHON − formulas cannot be excluded. Take C 18 as an example: 51 out of 87 C 18 H h O o N − n formulas should contain at least one amine group, either because their O atom numbers are not large enough to allow the assignment of -NO 2 for all N atoms or because some formulas (25 out of 87) were also detected in ESI+ mode. In total, 51 C 18 H h O o N − n formulas   (Fig. S6b). All the collective evidence above reveals that nitrogen-containing organic compounds in 10-18 nm particles during the I-NPF days are partly composed of long-chain amino alcohols, amino acids and so on.
In summary, a series of very distinctive CHON − and CHO − formulas was observed in 10-18 nm new particles during the I-NPF days. These formulas are characterized by some specific numbers of C atoms (i.e., and equivalent double bonds (DBE = 1 for CHO − and 1-4 for CHON − ). The monomer-dimer distribution pattern that can arise from particle-phase oligomerization (Pospisilova et al., 2020) was not observed for these formulas in the mass spectra. We also assume that C 18,30 H h O o N n and C 20,24,28,33 H h O o are not labile intermediates like ester hydroperoxides that may undergo fast decomposition in the particles or during the sample preparation process (Zhao et al., 2018a, b). To the best of our knowledge, such CHON − and CHO − formulas have not been reported by previous aerosol studies. The chemical composition of new parti-cles is completely decoupled from the CHO − and CHON − formulas around C 10 , C 15 and C 20 in 180-560 nm submicron particles, which might be originated from continental terpene emissions. Fewer O atoms in C 18,30 H h O o N n and C 20,24,28,33 H h O o formulas than those in submicron aerosols indicate that these compounds should be more freshly emitted into the atmosphere. The discontinuous chemical composition and PNSD spectrum (Fig. S2b) below and above 50 nm particle size reflect the fact that the further growth of new particles beyond 50 nm in local I-NPF events cannot be monitored by our stationary sampling strategy.
On the other hand, we observed more complicated distributions of CHO + and CHON + formulas in 10-18 nm new particles that are of relatively small individual intensity and are rather uniformly distributed in DBE vs. C diagrams. Like CHON − and CHO − , those CHO + and CHON + formulas also possess a larger number of C atoms (C > 19) than their counterparts in 180-560 nm submicron aerosols (Fig. 3). In fact, 21 out of 30 of the most abundant CHON + formulas contain two or more N atoms; this ratio 21/30 is higher than those in CHON − formulas. Generally speaking, CHO + and CHON + formulas represent carbonyls / alcohols / epoxides and amino alcohols / amino acids, respectively. One interesting finding about CHO + and CHON + is that they tend to form [M + Na] + adducts in small aerosols and [M + H] + adducts in large aerosols (Fig. 1c). This indicates that the CHO + and CHON + compounds in new particles during the I-NPF days should possess different basic functional groups from those in submicron particles.

Possible precursors and formation mechanism of organic compounds detected in 10-18 nm new particles during the I-NPF days
It is unrealistic to simply propose one out of a large number of possible structures for a formula with large C atom number (e.g., ≥ 18). Our strategy is to first explore the possible precursors of the newly found C Sporadic spikes of 10-18 nm particles that can be an indication of cooking and traffic emissions were not seen in the PNSD spectrum, because such human activities were rare around the site during the sampling period. We thus also exclude the possibility of cooking and traffic emissions. Previous field measurements of marine NPF precursor focused on volatile species like iodine (Stevanović et al., 2019), iodomethanes , dimethyl sulfonic acid (Yvon et al., 1996;Barone et al., 1996;Barnes et al., 2006) and LMW amines (Ning et al., 2019;Ge et al., 2011). So far there is no report about aliphatic compounds with C number ≥ 18 in either gas-phase or new particles (Cochran et al., 2017;Bikkina et al., 2019). Therefore, we consulted the literature that reported chemical compounds isolated from bio-logical tissues of algae, plankton or other marine organisms. Potential precursors are listed in Table 3.

Fatty acids
Fatty acids (FAs) are widely found in animals, plants and microbes (Moss et al., 1995). Plants have a higher content of unsaturated FAs than animals. C 14 -C 24 fatty acids, including both saturated and unsaturated, have long been observed in seaweed (Dawczynski et al., 2007). Very long chain FAs (C 24 -C 36 ) have been isolated from green algae, chlorella kessleri, sponges and marine dinoflagellate (Litchfield et al., 1976;Řezanka and Podojil, 1984;Mansour et al., 1999). C 18 Oleic acid, linoleic acid and linolenic acid are the most commonly found unsaturated FAs in macro algae. FAs with a carbon chain shorter than C 20 were used by atmospheric chemists as organic tracers of atmospheric aerosols from microbe or kitchen emission (Simoneit and Mazurek, 1982;Zheng et al., 2000;Guo et al., 2003;Rogge et al., 1991;DeMott et al., 2018;Willoughby et al., 2016). In our study, no saturated FAs were detected in 10-18 nm particles. Only 1.5 % of CHO − formulas can be assigned to unsaturated FAs (that is, include two O atoms, 14-28 C atoms and DBE = 3-6). Other CHO compounds observed in 10-18 nm particles contain more than two O atoms and thus can be assigned as the oxidized derivatives of FAs.
Possible oxidation schemes of two typical C 18 (C 18 H 30 O 2 , α-linolenic acid, three C=C double bonds in the carbon chain) and C 28 unsaturated FAs (C 28 H 52 O 2 , two C=C double bonds), for instance, are proposed in Figs. S7 and S8. The reaction of an unsaturated FA after the emission into the atmosphere is initiated by OH addition to a C=C double bond and subsequent O 2 addition to form a peroxy radical (Atkinson et al., 1995;Calvert et al., 2000). Depending on the level of NO and reactivity, four competitive pathways are available for peroxy radicals to produce CHO or CHON formulas observed in our study: reaction with NO to form a -ONO 2 group (pathway 1) or an alkoxy radical that further reacts with O 2 to form a carbonyl (-C=O, pathway 2), reaction with RO 2 radicals to form a hydroxyl (-OH) or a -C=O group (pathway 3), and successive intermolecular H-shift/O 2 addition autoxidation (Crounse et al., 2013;Vereecken et al., 2015) (pathway 4).
Pathways 1 and 2 add -ONO 2 , -OH and -C=O groups to the carbon chain but do not reduce the DBE of the FA precursor. We propose that pathways 1 and 2 are preferred for those FAs with higher reactivity with NO (e.g., α-linolenic acid; Fig. S7). α-linolenic acid oxidation in the atmosphere via sequential occurrences of pathways 1 or 2 yields a series of oxygenated and nitrated organic compounds, among which C 18 H 31 NO 6 , C 18 H 31 NO 8 , C 18 H 31 NO 10 , C 18 H 32 N 2 O 10 and C 18 H 33 N 3 O 4 are found in 10-18 nm particles. These formulas explain the circles with DBE = 4 and C number = 18 shown in Fig. 2a (DBE vs. C atom number diagram). Table 3. Possible precursors and their presence in marine biological sources and our aerosol samples. ND: not detected.
ND Fatty alcohols C 30 -C 32 mono-and diunsaturated alcohols and diols Yellow-green algae (Volkman et al., 1992)  The net outcome of sequential pathway 3 and 4 reactions is to add -OH and -C=O groups and reduce the DBE of the FA precursor. We propose that pathways 3 and 4 are preferred for those FAs (e.g., C 28 FA C 28 H 52 O 2 ) with higher reactivity with RO 2 (Fig. S8). The end products are a series of C 28 H 52 O 6-8 , C 28 H 54 O 4-7 and C 28 H 56 O 6-8 compounds, which can explain the circles with C number = 28 and DBE = 1-3 in Fig. 2b.
In addition to fatty acids, fatty alcohols such as C 30 -C 32 mono-and diunsaturated alcohols and diols have been detected in yellow-green algae (Volkman et al., 1992). Although these unsaturated alcohols were not detected in our 10-18 nm particles, we suppose that they or their metabolites in algae may undergo similar reactions like unsaturated FA to generate condensable oxygenated and nitrated fatty alcohols in the atmosphere. Hydroxy fatty acids (HFAs) are important constituents of lipid in marine microalgae (Gelin et al., 1997), bacteria (Kim and Oh, 2013), seaweed (Kendel et al., 2013;Blokker et al., 1998) and leaf surface of higher plants (Pollard et al., 2008). Among them, two saturated HFAs, C 20 H 40 O 3 and C 32 H 64 O 4 , were found in our 10-18 nm aerosol sample with S/N ratios 15 and 28. However, because saturated hydroxy fatty acids are not oxidizable via the pathways proposed in our study, they are assumed un-likely to be precursors of other formulas observed in 10-18 nm particles.

Nonprotein amino acids and amino alcohols
Quantum chemical calculations have showed that amino acids like glycine, serine and threonine are potential participants in atmospheric nucleation via interaction with sulfuric acid (Elm et al., 2013;Ge et al., 2018;Li et al., 2020). However, we did not observe any of 20 essential amino acids in 10-18 nm in either ESI+ or ESI− modes. One reason may be that most essential amino acids have molecular weight less than 150, which is below mass scan range of our FT-ICR-MS.
There are a number of records in the literature about longchain nonprotein amino acids or amino alcohols isolated from marine organisms or plant biomembrane (Ishijima et al., 2016;VanMiddlesworth et al., 1992;Jares-Erijman et al., 1993;Tsukamoto et al., 2001). They include saturated dihydroxy amino acid (C 18 H 37 NO 4 , DBE = 1, simplifungin), monounsaturated polyhydroxy amino acids in sphigolipids  , DBE = 2-3), polyhydroxy amino alcohols (C 16-28 H 33-53 N 1-2 O 1-5 , DBE = 1-3, sphingosine and its natural metabolites) and polyunsaturated amino alcohols (C 18 H 31 NO and C 18 H 29 NO,. All of these for-mulas were detected in 10-18 nm aerosols with S/N in the range of 10-800. More importantly, all those compounds that contain at least one amine group and one C=C double bond can be precursors of observed CHON formulas containing an amine group via the pathways that we showed above. As an example, the oxidation scheme of an amino alcohol C 18 H 31 NO with 4 C=C double bonds in the carbon chain is illustrated in Fig. S9. Similar to C 28 FA, C 18 H 31 NO undergoes successive intermolecular H-shift/O 2 additions to produce a series of RO 2 radicals with a hydroperoxyl group (-OOH) in its carbon chain. The subsequent pathway 3 reactions, as well as the decomposition of -OOH groups, add -OH and -C=O groups in the carbon chain. Because C 18 H 31 NO possesses as many as four C=C double bonds, sequential pathways 3 and 4 reactions produce a large number of oxidation products, among which 57 are found in the formula list detected in 10-18 nm particles (Fig. S9). These products C 13 ,C 18 H 33 NO 4,[6][7][8][9][10] explain perfectly the presence of a series of formulas with C number = 18, DBE = 0-4 and a -NH 2 group shown in Fig. 2a.

Volatility estimation
Based on the reaction mechanisms proposed above, it is possible to estimate the volatility change from potential precursors to their oxidation products. First, from the list of elemental formulas detected in 10-18 nm particles, we select 49 formulas with high intensities, including 14 CHON − formulas with peak intensity > 1.00 × 10 8 , 23 CHON + formulas with peak intensity > 3.00 × 10 7 and 12 CHO − formulas (DBE = 1) with peak intensities > 3.00 × 10 7 . Possible combinations of -COOH, -ONO 2 , -C=O, C=C double bond, -NH 2 and -OH groups are searched for every formula obeying two simple rules: CHON − and CHO − formulas must possess a carboxyl or hydroxyl group; CHON + formulas must possess an amino group. Saturation concentration (C * ) values of the 49 formulas were then predicted for all combinations using a simple group contribution method developed by Pankow and Asher (2008). On the other hand, the C * of their possible precursors, including unsaturated FAs, fatty alcohols, nonprotein amino acids or amino alcohols, were predicted by the same method.
As we can see in Table S4 in the Supplement, C * values of the 49 formulas fall into the range of ELVOCs (extremely low volatility organic compounds; 3 × 10 −9 -3 × 10 −5 µg m −3 ) and even ULVOCs (ultra low volatility organic compounds; < 3 × 10 −9 µg m −3 ), while C * values of their precursors are in the range of SVOCs (semivolatile organic compounds; 0.3-300 µg m −3 ) or LVOCs (3 × 10 −5 -0.3 µg m −3 ). The addition of functional groups reduces the volatility of precursors by 2-7 orders of magnitude and thus makes their oxidation products condensable onto new particles during the I-NPF event days. According to the definition of Schervish and Donahue (2020) and Simon et al. (2020), ULVOCs can even drive pure biogenic nucleation. Therefore, the analysis of precursor-product volatility partly supports our hypothesis about the molecular identity and formation mechanism of the formulas detected in 10-18 nm particles. It should be noted that the volatility of VOC oxidation products can be assessed with numerous existing parameterizations, which require either exact functional groups or only the molecular formula (Peräkylä et al., 2019). Their estimation can vary by up to several orders of magnitude. But this will not change the conclusion drawn here.

Speciation of iodinated OC
Organic iodine compounds hold the key to understanding aerosol iodine chemistry and its role in regulating the recycling of halogens to the gas phase. We identified 440 iodinated OC formulas from the four size bins during the I-NPF days (Table 2). Eighty percent of the 440 formulas contain one I atom, and the rest of them contain two I atoms. In terms of relative intensity, iodinated OC accounts for 2.6 %-9.5 % of OC in fine particles, but its fraction in coarse particles increases to 23.3 % in ESI− mode and 45 % in ESI+ mode. The size distribution of seven iodinated OC groups (i.e., CHOI − , CHONI − , CHOSI − , CHONSI − , CHOI + , CHONI + and CHNI + ) resembles those of non-iodinated OC groups (Fig. 1). If we replace I atom(s) with H atom(s) in a formula, 107 out of 440 replaced formulas are also found in the noniodine OC formula list.
Iodinated OC with intensity > 1.00 × 10 7 in the four size bins were shown in Fig. 4. The DBE vs. C diagram for 10-18 nm particles is characterized by (1) nine C 18 H h O o N n I formulas with DBE = 1-4 and (2) C 9 H 16 NO 3 I and its C 10 -C 13 homologues. Because these formulas were detected in ESI+ mode, they are most likely iodinated amino acids. The 32-56 nm particles accommodate the most abundant iodinated OC formulas, in terms of both formula number and relative intensity. Prominent formulas include (1) diiodo acetic acid C 2 H 2 O 2 I 2 , diiodomethane CH 2 I 2 , (2) iodinated C 21 carbonyls C 21 H 39 OI and C 21 H 41 OI, (3) iodinated C 21,25,27,29 alcohols or ethers with DBE = 0, (4) iodinated C 10 and C 15 terpene and sesquiterpene oxidation products, and (5) iodinated organic sulfate C 8 H 17 N 2 SO 8 I and C 21 H 43 SO 4 I. In addition, C 9 H 10 NO 3 I detected in this size bin (S/N ratio: 35) can be tentatively assigned to an iodinated amino acid iodotyrosine that has been observed in seaweed (Yang et al., 2016), implying direct contribution from seaweed emission to new particles.
In 180-560 nm particles, the majority of iodinated OC are C 6-9 aromatic CHOI + compounds with AI > 0.5 and DBE = 5-7. Both C and O atom numbers of these iodinated OC are smaller than those of mono-modally distributed CHO + compounds around C 15 in the same particle size (Figs. 3d and S3b). This implies that iodine has a strong tendency to aromatic compounds in submicron aerosols due to electrophilic substitution on aromatic rings. In 3.2-5.6 µm particles, iodinated OC features C 4 -C 6 CHO + and CHON + compounds with DBE = 3-6, which again have fewer C atoms than non-iodinated OC. Supporting evidence from AI shows these compounds are not aromatic. Coastal 3.2-5.6 µm particles can be sea salt particles formed during bubble bursting of seawater (Russell et al., 2010;Schmitt-Kopplin et al., 2012;Quinn et al., 2014;Wilson et al., 2015). However, Hao et al. (2017) showed that iodinated OC products from the reaction between iodine and seawater are highly unsaturated carboxylic-rich polyphenols with DBE = 4-14 and C atoms = 10-30. It is thus apparent that iodinated OC in 3.2-5.6 µm particles was not directly from iodinated organic compounds in the seawater.

Atmospheric implications
Due to the 71 % ocean coverage of the Earth's surface, marine aerosol generation is important in determining the Earth's radiative budget and climate change. Of the limited number of studies reporting coastal NPF, most have focused on iodine emission, oxidation and nucleation in the early stage of NPF. In principle, abundant low-volatility condensing vapors other than iodine are required in coastal environments for the growth of iodine clusters to CCN. This article reveals a new group of important organic compounds involved in this process. It is most likely that their precursors are emitted mutually with iodine from direct exposure of coastal biota to either the atmosphere or biologically active sea surface. If this is true, we suggest the results in our location can be extrapolated to other iodine-rich coastal lo-cations, as long as iodine-NPF can be observed. More fundamental field, laboratory and modeling studies are needed to determine (1) exact emission sources and source rates of these organic precursors, (2) their gas-phase intermediates and oxidation mechanisms in the atmosphere, and (3) their quantitative contribution to global and regional CCN numbers.
Data availability. All of the datasets related to the chemical formulas detected in this work can be accessed via https://doi.org/10.3974/geodb.2020.03.26.V1 (Yu, 2020). Any other data used in this publication are available from the corresponding author Huan Yu (yuhuan@cug.edu.cn) upon request.
Author contributions. HY designed and conducted chemical analysis. YW and HY did data analysis and wrote the paper. XH conducted the field sampling. BJ and YL did the FT-ICR-MS analysis. BK, ML, DX, JC and JZY reviewed and revised the paper.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "Marine organic matter: from biological production in the ocean to organic aerosol particles and marine clouds (ACP/OS inter-journal SI)". It is not associated with a conference.
Financial support. This research has been supported by the National Science Foundation of China (grant nos. 41975831 and 41675124) and the National Key Research and Development Program of China (grant no. 2016YFC0203100).
Review statement. This paper was edited by Manuela van Pinxteren and reviewed by two anonymous referees.