A two-component parameterization of marine ice-nucleating particles based on seawater biology and sea spray aerosol measurements in the Mediterranean Sea

<jats:p>Abstract. Ice-nucleating particles (INPs) have a large impact on the climate-relevant properties of clouds over the oceans. Studies have shown that sea spray aerosols (SSAs), produced upon bursting of bubbles at the ocean surface, can be an important source of marine INPs, particularly during periods of enhanced biological productivity. Recent mesocosm experiments using natural seawater spiked with nutrients have revealed that marine INPs are derived from two separate classes of organic matter in SSAs. Despite this finding, existing parameterizations for marine INP abundance are based solely on single variables such as SSA organic carbon (OC) or SSA surface area, which may mask specific trends in the separate classes of INP. The goal of this paper is to improve the understanding of the connection between ocean biology and marine INP abundance by reporting results from a field study and proposing a new parameterization of marine INPs that accounts for the two associated classes of organic matter. The PEACETIME cruise took place from 10 May to 10 June 2017 in the Mediterranean Sea. Throughout the cruise, INP concentrations in the surface microlayer (INPSML) and in SSAs (INPSSA) produced using a plunging aquarium apparatus were continuously monitored while surface seawater (SSW) and SML biological properties were measured in parallel. The organic content of artificially generated SSAs was also evaluated. INPSML concentrations were found to be lower than those reported in the literature, presumably due to the oligotrophic nature of the Mediterranean Sea. A dust wet deposition event that occurred during the cruise increased the INP concentrations measured in the SML by an order of magnitude, in line with increases in iron in the SML and bacterial abundances. Increases in INPSSA were not observed until after a delay of 3 days compared to increases in the SML and are likely a result of a strong influence of bulk SSW INPs for the temperatures investigated (T=-18 ∘C for SSAs, T=-15 ∘C for SSW). Results confirmed that INPSSA are divided into two classes depending on their associated organic matter. Here we find that warm (T≥-22 ∘C) INPSSA concentrations are correlated with water-soluble organic matter (WSOC) in the SSAs, but also with SSW parameters (particulate organic carbon, POCSSW and INPSSW,-16C) while cold INPSSA (T&lt;-22 ∘C) are correlated with SSA water-insoluble organic carbon (WIOC) and SML dissolved organic carbon (DOC) concentrations. A relationship was also found between cold INPSSA and SSW nano- and microphytoplankton cell abundances, indicating that these species might be a source of water-insoluble organic matter with surfactant properties and specific IN activities. Guided by these results, we formulated and tested multiple parameterizations for the abundance of INPs in marine SSAs, including a single-component model based on POCSSW and a two-component model based on SSA WIOC and OC. We also altered a previous model based on OCSSA content to account for oligotrophy of the Mediterranean Sea. We then compared this formulation with the previous models. This new parameterization should improve attempts to incorporate marine INP emissions into numerical models.                     </jats:p>

section_content: Abstract. Ice-nucleating particles (INPs) have a large impact on the climate-relevant properties of clouds over the oceans. Studies have shown that sea spray aerosols (SSAs), produced upon bursting of bubbles at the ocean surface, can be an important source of marine INPs, particularly during periods of enhanced biological productivity. Recent mesocosm experiments using natural seawater spiked with nutrients have revealed that marine INPs are derived from two separate classes of organic matter in SSAs. Despite this finding, existing parameterizations for marine INP abundance are based solely on single variables such as SSA organic carbon (OC) or SSA surface area, which may mask specific trends in the separate classes of INP. The goal of this paper is to improve the understanding of the connection between ocean biology and marine INP abundance by reporting results from a field study and proposing a new parameterization of marine INPs that accounts for the two associated classes of organic matter. The PEACETIME cruise took place from 10 May to 10 June 2017 in the Mediterranean Sea. Throughout the cruise, INP concentrations in the surface microlayer (INP SML ) and in SSAs (INP SSA ) produced using a plunging aquarium apparatus were continuously monitored while surface seawater (SSW) and SML biological properties were measured in parallel. The organic content of artificially generated SSAs was also evaluated. INP SML concentrations were found to be lower than those reported in the literature, presumably due to the oligotrophic nature of the Mediterranean Sea. A dust wet deposition event that occurred during the cruise increased the INP concentrations measured in the SML by an order of magnitude, in line with increases in iron in the SML and bacterial abundances. Increases in INP SSA were not observed until after a delay of 3 days compared to increases in the SML and are likely a result of a strong influence of bulk SSW INPs for the temperatures investigated (T = -18 • C for SSAs, T = -15 • C for SSW). Results confirmed that INP SSA are divided into two classes depending on their associated organic matter. Here we find that warm (T ≥ -22 • C) INP SSA concentrations are correlated with watersoluble organic matter (WSOC) in the SSAs, but also with SSW parameters (particulate organic carbon, POC SSW and INP SSW,-16C ) while cold INP SSA (T < -22 • C) are correlated with SSA water-insoluble organic carbon (WIOC) and SML dissolved organic carbon (DOC) concentrations. A relationship was also found between cold INP SSA and SSW nano-and microphytoplankton cell abundances, indicating that these species might be a source of water-insoluble or-Published by Copernicus Publications on behalf of the European Geosciences Union.

section_content: Ice-nucleating particles (INPs) are a subset of aerosol particles that are required for the heterogeneous nucleation of ice particles in the atmosphere. While extremely rare (Rogers et al., 1998), INPs greatly control the ice content of clouds, which is crucial to a range of climate-relevant characteristics including precipitation onset, lifetime, and radiative forcing (Verheggen et al., 2007). Despite their importance, the knowledge of INP sources and concentrations, particularly in marine regions, remains low as evidenced by the large uncertainties in modelled radiative properties of clouds (McCoy et al., 2015 (McCoy et al.,, 2016;; Franklin et al., 2013). 

While the ice-nucleating (IN) ability of marine SSA particles is less efficient than their terrestrial counterparts (De-Mott et al., 2016), modelling studies have shown that marine INPs are of particular importance in part due to the lack of other INP sources in such remote regions (Burrows et al., 2013; Vergara-Temprado et al., 2017). For this reason, recent studies have been conducted to better understand which SSA particles contribute to the marine INP population as well as the relationship between SSA emission and ecosystem productivity. Results from these studies suggest that the IN ability of SSAs is linked to the biological productivity of source waters, with higher productivity leading to greater IN activity (DeMott et al., 2016; Bigg, 1973; Schnell and Vali, 1976). For example, it has been shown that both the cell surface and organic exudate of the marine diatom Thalassiosira pseudonana can promote freezing at conditions relevant to mixed-phase clouds (Knopf et al., 2011; Wilson et al., 2015). More recently, mesocosm studies on phytoplankton blooms using two separate in-lab SSA-generation techniques have furthered the understanding of the connection between ocean biology and the IN activity of SSAs (McCluskey et al., 2017). In-depth chemical analysis of the artificially generated SSAs during this set of experiments has revealed marine INPs may be related to two classes of organic matter: a regularly occurring surface-active molecule type related to dissolved organic carbon (DOC) and long-chain fatty acids, and an episodic heat-labile microbially derived type (McCluskey et al., 2018a). 

As the understanding of the connection between ocean biology and marine INPs has improved, parameterizations for predicting marine INP abundance using readily available ocean parameters have been proposed. Wilson and co-authors (Wilson et al., 2015) identified a temperaturedependent relationship between TOC and ice-nucleating entity (INE) number concentrations in the surface microlayer (SML) from samples collected in the North Atlantic and Arctic ocean basins. They then extended this relationship from the ocean to the atmosphere to predict the abundance of INPs in SSAs based on model estimates of marine organic carbon aerosol concentrations. The parameterization was tested for the first time on field measurements of marine aerosol over the North Atlantic at Mace Head and was found to overestimate INP abundance in pristine marine aerosols by a factor of 4 to 100 at -15 and -20 (McCluskey et al., 2019). 

Despite the recent progress made in the understanding of marine INPs, there remains much room for improvement. To date, previous parameterizations have only been tested in the two field studies mentioned in the previous paragraph, underscoring the need for more real-world observations. Furthermore, the field studies conducted so far have taken place in regions of the ocean where biological productivity is high (i. e. North Atlantic and Southern Ocean). As modelling work has shown that the link between ocean biology and SSA organic content properties in oligotrophic waters differs from those in highly productive regions (Burrows et al., 2014) there is a need for more measurements in waters with low primary productivity. Finally, despite the finding that marine INPs may exist as two separate populations, no model has yet been proposed to account for this. 

This paper addresses the current gaps in the knowledge of marine INPs by ( 1
section_name: Introduction

section_content: In the frame of the PEACETIME project (http: //peacetime-project. org/, last access: 11 March 2021), an oceanographic campaign took place aboard the French research vessel (R/V) Pourquoi Pas? between 10 May-10 June 2017 with the purpose of investigating the processes that occur at the air-sea interface in the Mediterranean Sea. The cruise started in La Seyne, France, and travelled in a clockwise fashion between 35 to 42 • latitude and 0 to 20 • longitude. The observations and process studies performed on board in both the whole water column and the atmosphere are described elsewhere (Freney et al., 2020). Here, we focus on the measurements conducted to describe the SML, surface seawater (SSW), and aerosol properties.
section_name: Methods

section_content: SSW properties presented here were obtained from sampling at depths of 20 cm and 5 m. First, 21 parameters including various chemical properties, microbial assemblages, hydrological properties, and optical properties were monitored using the ship's underway system that continuously collected seawater at 5 m under the ship using a large peristaltic pump (Verder VF40 with EPDM hose). These measurements included counts of specific microbial classes (e. g. Synechococcus, Prochlorococcus, picoeukaryotes, nanoeukaryotes, microphytoplankton, high-phycoerythin-containing cells, coccolithophores, cryptophytes), as well as seawater biovolume, chlorophyll a (chl a), and POC concentrations. Chl a was determined from the particulate absorption spectrum line height at 676 nm after adjusting to PEACETIME chl a from high-performance liquid chromatography (HPLC). POC was estimated from the particulate attenuation at 660 nm using an empirical relationship specific to PEACETIME (POC = 1405. 1 × c p (660) -52. 4). For enumeration of phytoplankton cells, an automated Cytosense flow cytometer (Cytobuoy, NL) operating at a time resolution of 1 h was connected to the continuous underway seawater system. Particles were carried in a laminar-flow-filtered seawater sheath fluid and subsequently detected with forward scatter and sideward scatter as well as fluorescence in the red (FLR > 652 nm) and orange (FLO 552-652 nm) ranges. Distinction between highly concentrated picophytoplankton and cyanobacteria groups and lower concentrated nano-and microphytoplankton was accomplished using two trigger levels (trigger level FLR 7. 34 mV, sampling speed of 4 mm 3 s -1 analysing 0. 65 ± 0. 18 cm 3 and trigger level FLR 14. 87 mV at a speed of 8 mm 3 s -1 analysing 3. 57 ± 0. 97 cm 3 ). 

The second set of SSW measurements were made on seawater collected at ∼ 20 cm depth from a pneumatic boat that was periodically deployed at a distance of 2 km from the R/V to avoid contamination. The SSW was manually collected using acid-cleaned borosilicate bottles. From these discrete samples, microbial composition and cell abundance of the SSW was monitored as described in a companion paper (Tovar-Sanchez et al., 2019). Measurements included heterotrophic bacteria counts, high nucleic acid and low nucleic acid bacteria (HNA and LNA bacteria, respectively), total non-cyanobacteria-like cells (NCBL), cyanobacteria like cells (CBL), and total phytoplankton concentration (NCBL+ CBL). These were further segregated into size classes of small, medium, and large which roughly correspond to the pico-, nano-, and microsize classifications from the underway measurements. Trace metals (i. e. Cd, Co, Cu, Fe, Ni, Mo, V, Zn, Pb) were analysed by inductively coupled plasma mass spectrometry, although here we only report on Fe. Finally the abundance of DOC and marine gel-like particles, including transparent exopolymer particles (TEPs) and Coomassie stainable particles (CSPs), were also measured as described in the literature (Engel, 2009).
section_name: Surface seawater (SSW)

section_content: At the same time SSW samples were manually collected on the pneumatic boat, SML samples were also collected using a glass plate sampling method which has been previously described in the literature (Tovar-Sanchez et al., 2019). The glass plate was cleaned overnight with acid and rinsed with ultrapure Milli-Q water. Roughly 100 dips of the glass were conducted to collect 500 mL of SML water into 0. 5 L acidcleaned low-density polyethylene plastic bottles. The samples were then acidified on board to pH < 2 with ultrapuregrade hydrochloric acid in a class-100 HEPA laminar-flow hood. The same measurements done for the SSW samples (see above, Sect. 2. 1) were then made on the SML samples. Enrichment factor was calculated for relevant properties as the ratio of SML to SSW:

In addition to biological measurements, concentrations of immersion freezing mode INPs in SML samples (and a small number of SSW samples, n = 4) were measured between 22 May-7 June using an offline method described previously (Stopelli et al., 2014). Briefly, prior to acidification of the SML samples, additional aliquots were separated and stored in Corning Falcon 15 mL conical tubes and frozen at -20 • C until analysis. Before INP measurement, each aliquot was gradually defrosted and distributed into an array of 26 Eppendorf tubes filled up to 200 µL. The array was then immersed inside an LED-based Ice Nuclei Detection Apparatus (LINDA), and the number of ice-nucleating particles per litre (INP L -1 ) of SML water followed the method described in Stopelli et al. (2014), which was originally formulated by Vali (1971) :

where N total is the total number of tubes, N unfrozen is the total number of unfrozen tubes, and V tube is the volume of sample in each tube. The number of unfrozen tubes is calculated by first blank correcting the number of frozen tubes and then subtracting that value from the total number of tubes. We calculated uncertainty as the binomial proportion confidence interval (95 %) using the Wilson score interval. Samples were not corrected for salinity in this study.
section_name: Surface microlayer

section_content: Sea spray aerosols were generated using a sea spray generation apparatus which has been described previously (Schwier et al., 2015; Schwier et al., 2017). The characteristics of the setup were selected to mimic Fuentes et al. (2010). These parameters (water flow rates, plunging water depth, etc. ) have been shown to mimic nascent SSAs well. The apparatus consists of a 10 L glass tank with a plunging jet system. A continuous flow of seawater collected at 5 m depth using the ship's underway seawater circulating system (described above) was supplied to the apparatus. Particle-free air was passed perpendicular to the water surface at a height of 1 cm to send a constant airflow across the surface of the water. Aerosols were then either dried with a 1 m long silica dryer for online instrumentation (see Sect. 2. 3. 3), with a 30 cm silica gel dryer cascade impactor sampling with subsequent chemical analysis, or were sampled directly from the sea spray generator onto filters for INP analysis.
section_name: Artificially generated sea spray aerosol

section_content: Aerosol particles were also sampled onto PM1 quartz fibre filters mounted on a four-stage cascade impactor (10 L min -1 ) on a daily basis (24 h duration). Samples were then extracted in Milli-Q water by sonication for 30 min for the analysis of water-soluble components. The main inorganic ion abundance (i. e. SO 2- 4, NO - 3, NH + 4, Na +, Cl -, K +, Mg 2+, Ca 2+ ) was analysed via ion chromatography. An Ion-Pac CS16 3 × 250 mm Dionex separation column with gradient methanesulfonic acid elution was used for cations, while an IonPac AS11 2 × 250 mm Dionex column with gradient potassium hydroxide elution was used for anions. Watersoluble organic carbon (WSOC) and water-insoluble organic carbon (WIOC) were also determined. WSOC was measured after water extraction using a high-temperature catalytic oxidation instrument (Shimadzu; TOC 5000 A). Total organic carbon (which we now refer to as OC) was measured using a Multi N/C 2100 elemental analyser (Analytik Jena, Germany) with a furnace solids module. The analysis was performed on an 8 mm diameter filter punch, pre-treated with 40 µL of H 3 PO 4 (20 % v v -1 ) to remove contributions from inorganic carbon. WIOC was determined as the difference between OC and WSOC. Finally, we calculated organic mass fraction of SSA (OMSS) by taking the ratio of OM / (OM + SeaSalt), where OM is the sum of WSOM and WIOM, calculated as WSOM = WSOC × 1. 8 and WIOM = WIOC × 1. 4, and SeaSalt is the sum of inorganic ion abundance as determined above.
section_name: Offline PM1 filter analysis

section_content: INP concentrations were determined from filter-based samples of total suspended particles over a 24 h duration daily or from the average of two filters (day and night). The volume sampled on each filter averaged 8. 95 × 10 3 ± 2. 26 m 3 of air. The concentration of INPs in the SSAs was determined for the condensation freezing mode using a Dynamic Filter Processing Chamber (DFPC), which has been used in multiple previous studies and found to agree well with other INP monitoring instruments (DeMott et al., 2018; Hiranuma et al., 2019; McCluskey et al., 2018b). A full description of the instrument can be found in the literature (DeMott et al., 2018). Briefly, bulk SSAs formed using the plunging aquarium apparatus were impacted onto 47 mm nitrocellulose filters which were then placed on a metal plate coated with a smooth surface of Vaseline. Air entered the chamber and was sent through a cooling coil allowing it to become saturated with respect to water. Different supersaturations with respect to ice and liquid water can be obtained by controlling the temperatures of the filter and the air flowing across the filter. Filter air temperature combinations were set three different ways, all resulting in a supersaturation with respect to liquid water of 1. 02. The filter temperatures were -18, -22, and -25 • C (-15. 9, -19. 6, and -22. 3 • C for air temperature). Under these conditions, condensation freezing is expected to be the dominant freezing mode for INPs. It has been reported (Vali et al., 2015) that condensation freezing and immersion freezing are not distinguishable from one another. Filters were processed inside the DFPC for 15 min and monitored for formation of ice crystals upon activation of INPs. Based on sampling time and flow rate, the number of INPs per volume were calculated. We report an uncertainty of ±30 % based on previous reports of the DFPC (DeMott et al., 2018).
section_name: INPs

section_content: Particle size distribution and number concentrations of aerosols generated with the plunging apparatus were monitored using a custom-made differential mobility particle sizer (DMPS) preceded by a 1 µm size-cut impactor and X-ray neutralizer (TSI Inc. ). Total counts from the DMPS system were checked using a condensation particle counter (CPC, TSI3010). Using the DMPS, a total of 25 size bins ranging between 10-500 nm (dry particle electrical mobility diameter) were scanned over a 10 min time period. We then aver-aged the size distributions across each DFPC sampling period. Comparison of the total CPC-based SSA number concentration to the SSA number concentration derived from the DMPS revealed near unity, indicating nearly all of the particle number concentrations were captured by the DMPS. While studies typically present INP concentrations normalized by total SSA surface area, this was not possible in our experiment as the size distribution of supermicrometre particles was not monitored. However, in the supporting information, we do present a theoretical surface area normalized INP SSA calculation for comparison with other studies. The theoretical distribution was based on in situ particle number concentration measurements at Mace Head and open-ocean eddy correlation flux measurements from the Eastern Atlantic (Table S1 in the Supplement; Ovadnevaite et al., 2014), with the resulting surface area distribution shown in Fig. S1 in the Supplement.
section_name: Size distribution measurements

section_name: Results

section_content: Ice-nucleating particle characteristics were determined for the SSW, SML, and SSAs. Figure 1a at T = -15 • C). The enhancement on 4 June occurred on the same day as a dust deposition event which led to an enrichment of iron in the SML relative to the underlying water (see Sect. 3. 2). While only four SSW samples were analysed for INP concentrations, they exhibited similar concentrations and trends to those seen in the SML, with an observed maximum on 4 June (2. 4 × 10 3 INP L -1 at T = -15. 0 • C). Based on these four samples, no significant enrichment of INPs was observed in the SML compared to SSW, except during the dust deposition event when the SML concentration was enriched by a factor of 4. 5. 

Figure 1b shows the concentration of ice-nucleating particles in SSAs (INP SSA ) normalized by SSA particle concentration for particles with diameters between 0. 1 and 0. 5 µm at three different temperatures as observed by the DFPC. It should be noted that INP SSA measurements were conducted at colder temperatures than for the INP SML measurements due to differences between the LINDA and DFPC instruments. In general, the highest concentrations of INP SSA were observed at the beginning of the voyage, with an initial value of 4. 0 × 10 -9 INP SSA,-25 • C SSA -1 observed on 13 May, increasing to a maximum observed value of 1. 3 × 10 -8 INP SSA,-25 • C SSA -1 on 20 May. After May 20, a considerable drop in INP S SA, -25 • C con-centrations was observed. Concentrations remained low, albeit with slight fluctuations, before increasing again to 5. 2 × 10 -9 INP SSA,-25 • C SSA -1 on 7 June. It is also worth noting that the highest concentrations of INPs active at -18 • C (INP SSA,-18C SSA -1 ) were observed on this day. The increase in INP concentrations around the time of the dust deposition event in early June is similar to the trend observed for seawater INPs, albeit with a lag of at least 1 d (no observations of INP SSA were made on 6 June). 

Figure 2 shows the comparison of observed INP concentrations at various temperatures in the SML and SSW with those reported in previous studies. The concentrations we report here are lower than those from Arctic seawater samples reported by Irish et al. (2017 Irish et al. (, 2019) ) and from Arctic and North Atlantic seawater samples reported by Wilson et al. (2015). The difference can likely be attributed to the fact that eutrophic Arctic and North Atlantic seawater is more biologically active than the oligotrophic Mediterranean Sea. Our values agree well with those reported by Gong et al. (2020) who calculated INP concentrations in mid-latitude seawater off the coast of Cabo Verde. The authors of that study also posited that the low INP concentrations relative to Irish et al. (2017 Irish et al. (, 2019) ) and Wilson et al. (2015) were due to the lower biological activity of the oligotrophic seawater near Cabo Verde. As we did not measure the size distribution of particles larger than 500 nm, we cannot directly compare our INP SSA abundance to values cited in previous studies, where concentrations are typically normalized by SSA surface area which is dominated by supermicrometre particles. However, we were able to calculate a theoretical surface area distribution for particles between 0. 5-10 µm based on previous studies. The resulting surface area normalized INP concentrations and comparison with literature values is shown in the supporting information (Fig. S2 in the Supplement).
section_name: INPs in the seawater and SSAs

section_content: As described in the methods section, various seawater biogeochemical properties were monitored throughout the voyage for the SSW and SML. Plots of selected continuous measurements from the R/V's underway sampling system and discrete measurements from the pneumatic boat of relevant biogeochemical values are found in the supporting information (Figs. S3 and S4 in the Supplement, respectively). Biogeochemical properties are described in more detail in our companion papers (Guieu et al., 2020a; Tovar-Sanchez et al., 2019), and seawater gel properties will be discussed in an upcoming paper. Here, we present a broad summary of observed conditions. In general, surface waters were characterized by oligotrophic conditions as expected for the season. Bacteria concentrations ranged between 2 × 10 5 and 7 × 10 5 cells mL -1 in the SSW and were greatest at the start and end periods of the voyage. NCBL abundance followed a similar trend and ranged between 4. 0 × 10 2 -4. 0 × 10 3 cells mL -1. Observed DOC values ranged between 700-900 µgC L -1 and POC values between 42-80 µgC L -1, and these were within the range of expected values for the oligotrophic Mediterranean (540-860 µgC L -1 for DOC and 9. 6-104 µgC L -1 for POC; Pujo-Pay et al., 2011). SSW TEP concentrations ranged between 1. 2 × 10 6 and 1. 1 × 10 7 particles L -1, with CSP between 5. 6 × 10 6 and 9. 3 × 10 6 particles L -1, and will be discussed in a future paper. 

Enrichment factors (EFs) in the SML relative to the SSW remained low with an average of 1. 10 for DOC, 1. 07 for bacteria, and 1. 17 for NCBL. As POC was not measured in the SML, we cannot report its EF. TEP was typically enriched relative to the SSW, with an average EF of 4. 5, while CSP EF was on average 2. 7. Of importance, the dust deposition event that occurred on 4 June lead to a drastic increase in SML dissolved iron relative to the SSW (EF ∼ 800). This deposition event had important impacts on the biology of the surface seawaters, which is the focus of another paper (Freney et al., 2020). As a result, TEP EF increased to 17, bacteria EF increased to 1. 5, and NBCL EF to 2. 4. We next discuss the correlations between INP abundances and biogeochemical properties in the following sections.
section_name: Correlations between INP and biogeochemical conditions

section_content: Table 1 shows the correlation statistics between INP SML,-15 dust event (for micro-NCBL SSW and TEP SSW ) or due to the dust event (for TOC SSW, DOC SSW, and nanoeukaryote cell abundances). This points to a non-proportional transfer of each species from the bulk seawater to the SML relative to one another. Given the high p values and weak correlation coefficients, it is likely that INP SML is not strongly driven by the properties of the underlying SSW. Rather, we posit that INPs in the SML are dictated by SML properties, as shown in the following paragraph. 

Figure 4 shows scatterplots of statistically significant relationships between INP SML,-15 • C concentrations and various SML properties. INP SML,-15 • C is most strongly positively correlated with dissolved iron (r = 0. 99), TEP EF (r = 0. 95), and bacteria EF (r = 0. 93). However, these relationships are skewed by the outlier due to the drastic increase in iron observed on 4 June (Fig. S2a ) from the dust deposition event, as described previously. It is difficult to discriminate between the dust and biological impact on the INP SML,-15 • C, as dust is known to have good INP properties while also being capable of fertilizing the surface ocean with dissolved iron, leading to concomitant increases in biological activity. It is also possible that the dust deposition led to increased abundance of terrestrial OC, which would exhibit different INP activity. When considering days before the dust event, INP SML,-15 • C is only significantly correlated with dissolved iron (r = 0. 91) and TOC in the SML (r = -0. 93). We note that while no longer statistically signif-icant for pre-dust days, moderate correlations were still observed between INP SML,-15 • C and total NCBL (r = 0. 48), HNA bacteria (r = 0. 78), and total bacteria (r = 0. 64). Previous reports examining the correlation between INP and microbial abundance have yielded mixed results. For example, a report of INPs in Arctic SML and SSW found no statistically significant relationship between the temperature at which 10 % of droplets had frozen and bacteria or phytoplankton abundances in bulk SSW and SML samples (Irish et al., 2017). However, recent mesocosm studies using nutrient-enriched seawater found that INP abundances between -15 and -25 • C in the aerosol phase were positively correlated with aerosolized bacterial abundance (McCluskey et al., 2017). 

A previous study by Wilson and co-authors presented an INP parameterization (hereafter termed W15) based on a positive relationship between seawater TOC and INP abundance in Arctic, North Pacific, and Atlantic SML and SSW (Wilson et al., 2015). Total organic carbon in the SML (TOC SML µgC L -1 ), derived here as the sum of POC in the SSW (POC SSW ) and DOC in the SML (DOC SML ), was poorly correlated with INP SML,-15 • C (r = 0. 31, p = 0. 50). Figure 5 indicating the TOC SML in Mediterranean waters is less IN active at these temperatures than predicted by the W15 parameterization. 

In agreement with our findings, a recent study found that the W15 model overpredicted observed INP concentrations in the aerosol phase during two separate mesocosm experiments (McCluskey et al., 2017) by assuming the INP / TOC ratio in the SML was preserved in the aerosol phase. The authors of that study speculated that the overprediction by the W15 model was due to the fact that it does not account for the complex transfer mechanism of organic matter from the SML to the aerosol phase. Our results here show that the overprediction by W15 persists even when calculating INPs in the SML and therefore the overprediction may be due to other factors beyond the transfer of organic matter from the SML to the atmosphere. We stress, however, that the TOC value used in this study was derived using DOC SML and POC SSW values as POC measurements in the SML were not conducted. As there typically exists an enrichment of organic matter in the SML relative to the bulk seawater, it is possi- ble that the POC SSW we used to calculate TOC SML was below the actual POC content in the SML, thus underestimating TOC SML. However, if this was the case, a higher abundance of TOC SML would only further increase the overprediction of W15 relative to our observations. Finally, it is possible that the oligotrophic nature of Mediterranean waters results in a pool of TOC with a different chemical composition than what is observed in more biologically productive waters such as the Arctic and Atlantic. For example, the pool of TOC during this study was dominated by DOC and featured low POC content, presumably due to low biological productivity. 

In
section_name: Correlations between INP SML abundance and seawater properties

section_content: In the following section, we compare INP SSA at various temperatures with seawater and SSA properties. Submicrometre particle concentrations ranged between 1000-3000 particles cm -3 (Fig. S5 in the Supplement), and its dependence on seawater biology is further explored in a separate paper (Sellegri et al., under revision). For comparison with seawater properties, INP SSA was first normalized by SSA particle concentration (0. 1 < D p < 0. 5 µm). 

Table 2 shows the correlation statistics between INP SSA normalized by SSA particle number concentration and select conditions in the SML for relationships that were statistically significant. Figure 6 shows the corresponding scatter plots for these relationships. We also tested for correlations on days not affected by the dust event (i. e. days before 4 June), and their statistics are shown in parentheses in Table 2 and Fig. 6. Surprisingly, there were no significant correlations between INP SSA,-18C and conditions in the SML, including TEP and CSP abundance and enrichment factors, bacteria abundance and enrichment factors and nor with INP SML as measured by the LINDA instrument. This is somewhat unexpected considering INPs in the SML at -15 • C were correlated with SML phytoplankton and bacteria counts, which are all expected to transfer efficiently from the SML to the aerosol phase, an assumption widely used in the modelling community. Similarly, -22 • C INP SSA had no significant correlations with SML variables, except for TEP EF which was positively correlated (r = 0. 69) when only considering days before the dust deposition event. At -25 • C, INP SSA was found to be significantly correlated with DOC SML and TOC SML on all days (r = 0. 82 and r = 0. 80 for DOC and TOC, respectively). When examining only pre-dust-event days, the significant correlations included DOC enrichment as well as nano-and micro-CBL. 

Table 3 and the corresponding scatter plots in Fig. 7 show that a weak correlation exists between INP SSA active at -18 • C and POC SSW for all days but becomes significant and strong for days not including the dust event. This points to the possible interference of a different class of organic carbon (e. g. terrestrial OC) or the introduction of some other IN-active material (e. g. dissolved iron) which masks the impact of the original pool of POC SSW on INP concentrations. INP SSA,-18C are also significantly correlated with INP SSW,-16C (results not shown), but with a sample size of n = 4 this finding requires further validation. Nonetheless, this result could indicate that INP SSA at this tempera- 2. R 2 values for all days are shown in each plot, with R 2 values for days not including the dust deposition event (i. e. days before 4 June) in parentheses. R 2 for statistically significant relationships are shown in bold. 

Table 2. Correlation statistics between INP SSA and properties in the SML, where p is the p-value test for significance and r is the Pearson correlation coefficient. Values in parentheses are calculated for days before the dust deposition event (i. e. days before 4 June). Values that are not statistically significant (p > 0. 05) are italicized.
section_name: Correlations between INP SSA abundance and observed SSAs and seawater conditions

section_content: p all days (p pre-dust ) r all days (r pre-dust )

-18 Table 4 and Fig. 8 show the significant correlations between INP SSA and SSA properties. A time series of SSA chemical properties is shown in Fig. S6 in the Supplement. A positive correlation was observed between INP SSA,-18C and SSA OC as well as in the ratio of SSA water-soluble organic carbon to organic carbon (WSOC / OC). The correlation between WSOC / OC and INP SSA,-18C makes sense given the finding that INP SSA,-18C was correlated with POC SSW. A higher WSOC / OC value would suggest a higher fraction of soluble organics which would be expected to transfer to the atmosphere from the bulk SSW rather than the SML due to their high solubility. 

Figure 8 and Table 4 also show that INP SSA,-25C had a significant correlation with WIOC and organic mass fraction of sea spray (OMSS) (r = 0. 58 and r = 0. 65, respectively). As mentioned above, INP SSA,-25C was found to be correlated with various microbes in the SSW, specifically Prochlorococcus, coccolithophores, and nano-and micro-NCBL. Phytoplankton are known for their ability to produce extracellular polymeric substances (Thornton, 2014), and a previous mesocosm experiment showed microbially derived longchain fatty acids were efficiently ejected from the seawater as SSAs, increasing the fraction of highly aliphatic WIOC (Cochran et al., 2017). A separate paper discusses the trend and controls on SSA chemical composition, linking the different classes of organic carbon in submicrometre SSAs to seawater chemical and biological properties (Freney et al., 2020). In this work, OMSS was linked to POC SSW and the coccolithophore cell abundance. In light of this and given the correlation of INP SSA,-25C with seawater microbial abundance and with SSA OMSS and WIOC, it seems likely that INP SSA at this temperature is related to the exudates of phy- (Wilson et al., 2015, and McCluskey et al., 2018a, respectively). To improve upon existing models, we formulated various parameterizations consisting of different time periods, features, and number of components for temperature ranges. Predictor features were chosen based upon their correlation with INP concentrations as described in the previous section. Single-component parameterizations in which INPs across all three temperatures were linked with the same features were compared with two-component parameterizations in which INPs were split into warm and cold categories, each having their own predictor features. Finally, we developed and compared an altered version of the W15 model to account for the oligotrophic seawater of the Mediterranean Sea, as the existing model was formulated from observations of eutrophic waters. An altered version of the MC18 model for oligotrophy is presented in the Supplement (Fig. S7 in the Supplement), based on calculations of INP concentrations normalized by theoretical total SSA surface area. Each parameterization was recalculated using data across all days of the cruise as well as for only days before the dust deposition event in order to determine the impact of the dust event on the ability to predict INPs. The complete set of parameterizations and their associated fit metrics (R 2 and R 2 adj. ) are given in Table S2 in the Supplement. 

Figure 9a shows observed vs. predicted INP SSA for the W15 model. Similar to our results for seawater INPs (Fig. 5 ), a large overprediction is found relative to our observations when using W15. We also present re-calculated best-fit lines to data using the same features as in W15 (i. e. OC) in order to account for possible changes due to the oligotrophic nature of the Mediterranean Sea. We term this parameterization the altered Wilson fit for oligotrophy, which is given by the following:

The results for this fit are shown in Fig. 9a

Figure 9b and c show the results of our single-component model using POC SSW and the two-part model which uses SSA WIOC and OC and considers the separate temperature classes of INPs. The adjusted R 2 for each model on the log-transformed INP abundance were R 2 adj = 0. 402 for the single-component model using POC SSW and R 2 adj = 0. 60 for the two-component model using OC and WIOC. This result reveals that the two-component method performs as well as the altered Wilson parameterization. Each parameterization's fit to the data is improved when considering pre-dust ±3. 5×10 6 INP (gC) -1. This points to a depletion in the abundance of INP-active material by a factor of 16 as it transfers from the seawater to the SSAs, which is typically assumed to be negligible in modelling studies. However, when available, using a ratio of INP SSW / TOC SSW to predict sea-spray-originating INPs in the atmosphere seems a better approach than using the ratio INP SSW /NaCl SSW. Finally, we remind the readers that the two-component parameterization uses results of SSA chemistry for submicrometre particles only. As previous studies have shown that the overwhelming majority of SSA OC is found in the submicrometre phase (Gantt and Meskhidze, 2013), we argue that our analysis of WIOC, WSOC, and OC concentrations in submicrometre SSAs is representative of the whole size range of SSAs. 

In this paper we have presented results from the month-long PEACETIME cruise which took place in the Mediterranean Sea during the spring of 2017, which was characterized with a dust wet deposition event that occurred towards the end of the cruise. We find that the INP concentrations measured in the seawater are in agreement with previous studies on oligotrophic waters (Gong et al., 2020). We next investigated the As the cold INPs are also correlated to the SSW nano-and micro-NCBL cell abundance as well, we hypothesize that these classes of phytoplankton produce surface-active waterinsoluble organic matter that is IN active at these temperatures and is transferred to the atmosphere via the SML. Unfortunately, we do not have measurements of the "colder" temperatures of INPs in the SML to check this hypothesis. 

We finally proposed a single-component model linking INP / SSA to POC SSW and a two-component model linking warm INPs to SSA OC and cold INPs to SSA WIOC. Both models utilize features that are readily approximated either from satellite data, biogeochemical models, or from existing parameterizations and observations (Aumont et al., 2015; Rasse et al., 2017; Albert et al., 2010). We then showed these parameterizations fit the data much better than previous single-component model based solely on OC content (W15) developed from studies of more biologically active waters. We also re-calculated parameterizations based on SSA OC content but for the oligotrophic Mediterranean Sea. The parameterization using SSA OC content fits almost as well as the two-component model using SSA OC and WIOC. However, given the results of correlation analysis with SSA properties as well as results from previous studies indicating a dual composition of INPs, we believe the two-component model should help improve attempts to incorporate marine INP emissions into numerical models.
section_name: Variable

section_content: Data availability. Underlying research data are being used by researcher participants of the PEACETIME campaign to prepare other papers, and therefore data are not publicly accessible at the time of publication. The Guieu et al. (2020b) study will be accessible at https://www. seanoe. org/data/00645/75747/, once the special issue is completed (all papers should be published by June 2021). 

Supplement. The supplement related to this article is available online at: https://doi. org/10. 5194/acp-21-4659-2021-supplement.

section_content: Author contributions. KS designed the emission experiments, and KS and AN performed them. AN performed the seawater, SML, and sea spray aerosol sample INP analysis, with the help of FB and GS. JVT performed the formal INP data analysis and merged data set investigation, with help from KS and EF. BZ and AT-S performed SML sampling and analysed the sample together with ARR. AE, BZ, and ATS contributed to formal analysis of SML data. MR and EF analysed the chemical composition of sea spray aerosol samples. MT, IO, and JD contributed to the analysis of the chemical and biological composition of the SSW. CG contributed to funding acquisition and management of the PEACETIME project. JVT and KS wrote the original draft. All authors contributed to writing and revision of the work. 

Competing interests. The authors declare that they have no conflict of interest. Special issue statement. This article is part of the special issue "Atmospheric deposition in the low-nutrient-low-chlorophyll (LNLC) ocean: effects on marine life today and in the future (ACP/BG interjournal SI)". It is not associated with a conference. 

Acknowledgements. This study is a contribution to the PEACE-TIME project (https://doi. org/10. 17600/17000300 and http:// peacetime-project. org, last access: 19 March 2021), a joint initiative of the MERMEX and ChArMEx projects and to the Sea2Cloud project. The PEACETIME cruise was endorsed as a process study by GEOTRACES and is also a contribution to IMBER and SOLAS. Sea2Cloud is endorsed by SOLAS. We thank the captain and the crew of the R/V Pourquoi Pas? for their professionalism and their work at sea. The underway optical instrumentation was provided by Emmanuel Boss's group. Review statement. This paper was edited by Anne Perring and reviewed by three anonymous referees.