The Potential Impact of Saharan Dust and Polluted Aerosols on Microbial Populations in the East Mediterranean Sea, an Overview of a Mesocosm Experimental Approach

Recent estimates of nutrient budgets for the Eastern Mediterranean Sea (EMS) indicate that atmospheric aerosols play a signiﬁcant role as suppliers of macro- and micro- nutrients to its Low Nutrient Low Chlorophyll water. Here we present the ﬁrst mesocosm experimental study that examines the overall response of the oligotrophic EMS surface mixed layer (Cretan Sea, May 2012) to two different types of natural aerosol additions, “pure” Saharan dust (SD, 1.6 mg l − 1 ) and mixed aerosols (A—polluted and desert origin, 1 mg l − 1 ). We describe the rationale, the experimental set-up, the chemical characteristics of the ambient water and aerosols and the relative maximal biological impacts that resulted from the added aerosols. The two treatments, run in triplicates (3 m 3 each), were compared to control-unamended runs. Leaching of ∼ 2.1–2.8 and 2.2–3.7 nmol PO 4 and 20–26 and change compared to the SD amendments (per a similar amount of mass addition), especially regarding BP and PP. We speculate that despite the co-limitation of P and N in the EMS, the additional N released by the A treatment may have triggered the relatively larger response in most of the rate and state variables as compared to SD. An implication of our study is that a warmer atmosphere in the future may increase dust emissions and inﬂuence the intensity and length of the already well stratiﬁed water column in the EMS and hence the impact of the aerosols as a signiﬁcant external source of new nutrients.

Recent estimates of nutrient budgets for the Eastern Mediterranean Sea (EMS) indicate that atmospheric aerosols play a significant role as suppliers of macro-and micro-nutrients to its Low Nutrient Low Chlorophyll water. Here we present the first mesocosm experimental study that examines the overall response of the oligotrophic EMS surface mixed layer (Cretan Sea, May 2012) to two different types of natural aerosol additions, "pure" Saharan dust (SD, 1.6 mg l −1 ) and mixed aerosols (A-polluted and desert origin, 1 mg l −1 ). We describe the rationale, the experimental set-up, the chemical characteristics of the ambient water and aerosols and the relative maximal biological impacts that resulted from the added aerosols. The two treatments, run in triplicates (3 m 3 each), were compared to control-unamended runs. Leaching of ∼2.1-2.8 and 2.2-3.7 nmol PO 4 and 20-26 and 53-55 nmol NO x was measured per each milligram of SD and A, respectively, representing an addition of ∼30% of the ambient phosphate concentrations. The nitrate/phosphate ratios added in the A treatment were twice than those added in the SD treatment. Both types of dry aerosols triggered a positive change (25-600% normalized per 1 mg l −1 addition) in most of the rate and state variables that were measured: bacterial abundance (BA), bacterial production (BP), Synechococcus (Syn) abundance, chlorophyll-a (chl-a), primary production (PP), and dinitrogen fixation (N 2 -fix), with relative changes among them following the sequence BP>PP≈N2-fix>chl-a≈BA≈Syn. Our results show that the "polluted" aerosols triggered a relatively larger biological

INTRODUCTION
In low-nutrient low-chlorophyll (LNLC) marine environments, nutrient and trace metal inputs via atmospheric aerosols are considered important sources of macro and micro nutrients (Duce et al., 1991(Duce et al., , 2008Jickells et al., 2005;Kanakidou et al., 2012), fueling microbial production and influencing the bacterioplankton community structure (Moore et al., 2013;Guieu et al., 2014a;Chien et al., 2016;Rahav et al., 2016a). The Eastern Mediterranean Sea (EMS), located in the so-called dustbelt (Astitha et al., 2012), is considered extremely oligotrophic (reviewed in Siokou-Frangou et al., 2010) and is strongly influenced by natural desert sources that contain some of the highest atmospheric dust (mineral aerosol) concentrations near the Earth's surface (Klingmüller et al., 2016). The EMS region has been identified in the last decade as a hot-spot of climate change, showing a decrease in precipitation and an increase in the temperature and annual number of unusually hot days (Hoerling et al., 2012;Lelieveld et al., 2012;IPCC, 2014). It is suggested that the increasing temperatures together with the decreasing relative humidity of the last decade, have promoted soil drying, leading to increased dust emissions in the EMS, a process that is expected to continue in the future due to climate change (Klingmüller et al., 2016) and consequently supply more macro and micro nutrients into its surface oligotrophic water.
The EMS is also exposed to relatively high pollution levels of aerosols (Lelieveld et al., 2002), which interact with the already high background levels of natural mineral particles. The bioavailability of nutrients and trace metals from aerosols is related to both, the original aerosol chemical and mineralogical composition and the interactions and chemical transformations during atmospheric transport Baker and Jickells, 2006;Mackey et al., 2015). Recent studies have shown that the exposure to acid processes in the atmosphere can also change the phosphorus (P) bioavailability of dust/aerosols (Nenes et al., 2011;Bougiatioti et al., 2016).
Atmospheric deposition of nutrients to LNLC provinces is particularly important in regions where there is little input of new nutrients from other external sources (e.g., Jickells et al., 2005;Duce et al., 2008), as is the case in the EMS (Herut et al., 1999Krom et al., 2004). Such inputs are more significant during stratified periods, as they enable new production in a period of generally low autotrophic production (Guieu et al., 2010). In addition to new production, it has been shown that leachable nutrients from atmospheric inputs may enhance N 2 fixation, changes in phytoplankton species composition and carbon sequestration (e.g., Mills et al., 2004;Guo et al., 2012;Moore et al., 2013;Bressac et al., 2014;Guieu et al., 2014b;Rahav et al., 2016a). While a low availability of macronutrients (N, P) and metal micronutrients (e.g., Fe, Co) can limit or colimit phytoplankton growth in the ocean (e.g., Moore et al., 2013), high concentrations of some metals (e.g., Cu or Al) can be toxic to phytoplankton (e.g., Paytan et al., 2009;Jordi et al., 2012;Krom et al., 2016). Yet, the impact of these dust inputs on microbial populations has not been fully investigated in the EMS. The response of oligotrophic systems to dust inputs, whether as positive or as negative feedbacks to autotrophic and heterotrophic production (and thus to biogeochemical cycling), must therefore be further examined (Guieu et al., 2014a).
To date, these influences were assessed in the EMS mainly by on-board dust enrichment microcosm experiments Ternon et al., 2011) or in situ observations (Rahav et al., 2016a), which showed an increase in primary production and in heterotrophic bacterial activity, while no previous mesocosm experiment has been performed. The use of large-scale mesocosms may be particularly important to reduce the effects of bottle enclosure that may change the autotrophic-heterotrophic biomass ratio, especially in oligotrophic regions such as the EMS (Calvo-Díaz et al., 2011). Similarly, a mesocosm experiments in the Western Mediterranean Sea examined the addition of a dust analog (treated Saharan soil) to surface seawater (Guieu et al., 2014b). In that experiment, the Saharan dust amendment strongly stimulated primary production and algal biomass and altered the autotrophic phytoplankton communities (Guieu et al., 2014b). Furthermore, these additions affected the structure, diversity and functioning of the microbial food web (Pulido-Villena et al., 2014) while also increasing the N 2 fixation rates (Ridame et al., 2013).
Here we present an overview of a mesocosm experimental study that examines the microbial response of the oligotrophic EMS surface mixed layer (Cretan Sea, May 2012) to two different types of natural aerosol additions; "pure Saharan dust" and "mixed aerosols" (a natural mixture of desert dust and polluted particles). To the best of our knowledge, this is the first mesocosm study in which naturally collected aerosols were used in the Mediterranean Sea. The added particles were designed to mimic the impact of an intense dry atmospheric deposition to the upper mixed layer (1-1.5 mg of dust l −1 ) on the physiology and biomass of the ambient microbial populations of the EMS, using natural aerosols collected across the Levantine basin (Crete and Israel). The two aerosol types were compared to control runs in triplicates for a total duration of 8 days. This article provides a brief overview of the design, rationale, characteristics of the aerosols and the principal responses of the microbial variables in the mesocosm experiment, as part of a Research Topic entitled "Impact of aerosols (Saharan dust and mixed) on the East Mediterranean oligotrophic ecosystem, results from experimental studies." Other components of the mesocosm experiment, including microbial community structure and temporal dynamics, microbial biodiversity and function are discussed in other articles of this special issue (Guo et al., 2016;Rahav et al., 2016b; and others in this special issue).

Mesocosm Experimental Design and Sampling
An aerosol-enrichment mesocosm experiment was performed at the CRETACOSMOS facility of the Hellenic Centre for Marine Research (HCMR, www.cretacosmos.eu) in Crete, Greece, during the 10-18 of May 2012. The facility consists of a 350 m 3 land-based concrete pond, 5 m deep, supplied with continuous seawater flow-through in order to maintain ambient surface water temperature. The experiment was carried out using surface (∼10 m depth) seawater that was collected using a rotary submersible pump placed onboard the R/V Philia from a location 5 nautical miles north of Heraklion (Crete, Greece, 35 • 24.957 N, 25 • 14.441 E) at a bottom depth of 170 m during the 8-9 of May 2012. The collected seawater was equally distributed by gravity into nine food-grade polyethylene mesocosm bags to ensure the homogeneity of the collected seawater between bags. The mesocosms were mounted on aluminum frames (1.12 m diameter) attached to the pool's walls. Each mesocosm had a total volume of 3 m 3 (Figure 1). The mesocosms were gently mixed throughout the experiment, using an airlift pump to avoid stratification. The mesocosms were covered with a twolayer lid in order to protect them from natural atmospheric aerosol depositions during the experiment and mimic the light conditions at a 10 m water depth. HOBO data loggers (ONSET Corporation) were installed to measure the temperature and light (irradiation) in each bag as well as in the main pond.
On May 10th at 08:30, samples were collected from all the mesocosms prior to the dust/aerosol additions to serve as reference conditions. At 11:45, the dust/aerosols were added to the mesocosms and the first sampling was carried out 3 h later at 14:45. Each mesocosm was sampled daily during the morning (08:30) from May 11th to 15th, 2012, and then once every 2 days until May 18th, apart from parameters that required more frequent monitoring, which were sampled every day throughout the experiment. Acid-washed silicone tubes were used for transferring the samples collected at about 50 cm depth into similarly treated 20 L Nalgene containers.
Two different types of natural additions were performed in triplicates: "pure" Saharan dust, (labeled as SD) and mixed aerosol containing a natural mixture of desert dust and polluted European particles (labeled as A). The additions of SD and A resulted in a final concentration of 1.6 and 1 mg l −1 , respectively, in the mesocosm bags (Table 1). Triplicate control (labeled C) mesocosm treatments were run in parallel. Each replicate addition was performed by pouring a concentrated mixture of the aerosol/dust with filtered (0.2 µm) seawater into the bags. This mixture was prepared just prior to the addition in precleaned (10% hydrochloric acid washed) 100 ml polyethylene bottles. In order to mimic an intense dust storm effect, it was decided to add quantity equivalent to final concentration of 1.6 mg l −1 of aerosols in the mesocosm bags. Previous studies yielded a deposition of ∼1 mg of dust l −1 in the upper 5 m mixed layer during dust storm events in the EMS Rahav et al., 2016a). Yet, due to the limited mass availability of aerosols representing an European origin (Figure 2), we prepared a mixture of desert/mineral dust and polluted aerosols (treatment A), allowing for an addition equivalent to only 1 mg l −1 in the relevant mesocosms (Table 1, Figure 1). Aerosol chemical composition-Analyses of the chemical compositions of SD and A aerosols were carried out after total digestion with HF and aqua regia, following the procedure of ASTM (1983) and Herut et al. (2001). The digested samples were diluted in 25 ml with Milli-Q water and filtered through a Whatman 42 filter paper. Briefly, 1 ml of aqua regia solution and 4 ml of hydrofluoric acid were added to ca. 0.15 g of dry aerosol in 125 ml plastic bottles (that can withstand temperatures of up to 1300 • C on a sand bath). Prior to dilution, 5 ml of saturated boric acid was added to the resultant solution. The concentrations of trace elements were measured using an Agilent 280FS atomic absorption spectrometer, a graphite furnace Agilent 240Z AA and ICP-MS. Major elements were measured on an ICP-AES. The bottle blanks were usually less than 2% of the measured concentrations. The accuracy of the method was evaluated based on analyses of International Certified Reference Materials: MESS-4, SRM 2702, and IAEA 158. All elements gave results within 94-110% of the certified values, except for Cd, which gave an 80% recovery.

Aerosol Collection and Leaching Experiments
Dry deposited material was collected in Crete (Heraklion and Sambas) and Israel (Beit Yannay and Tel Shikmona, see Herut et al., 2002) during major Saharan dust storms and during periods that showed the transport of air masses was from Europe (Table 1, Figure 2). These two sites represent aerosols across the Levantine basin in the SE Mediterranean Sea. Aerosol origin and route was tracked by calculating 3-day back trajectories using the NOAA HYSPLIT model from the Air Resources Laboratory at 1000 and 3000 m altitude levels (Figure 2). We needed between 9 and 15 g of material for the additions to the different mesocosm replicates to obtain the final concentrations in the mesocosm bags (1 and 1.6 mg l −1 ), and therefore material collected on different sampling dates was pooled, representing the two types of treatments (Table 1). However, while for the SD treatment sufficient material was collected, only a limited amount could be collected to represent a clear dominated polluted European  derived aerosols. Therefore, treatment A represents a mixture of European as well as desert aerosols (Table 1).
Leaching experiments were performed according to two different methodologies using sterile surface seawater from different origins, in the nutrient laboratories of the Israel Oceanographic and Limnological Research (Haifa) and at the University of Leeds ( Figure S1). In Haifa, 250 ml plastic containers were pre-cleaned with 10% hydrochloric acid and washed afterwards with Milli-Q water. Each bottle was filled with filtered (0.2 µm), aged and poisoned (with 50 µL chloroform) SE Mediterranean surface seawater (collected with the R/V Shikmona, 40 km off the Israeli coast), to which 30 mg of SD or A powder was added. The bottles were covered with aluminum foil and shaken at room temperature for a total of 48 h. Subsamples (10 ml) collected at 0, 0.75, 2, 6, 24, and 48 h were pipetted using pre-cleaned (10% HCl) syringes through a 0.2 µm filter and immediately analyzed for nitrate+nitrite and phosphate concentrations. Nutrients were measured with a Seal Analytical AA-3 system (Kress et al., 2014;Ozer et al., in press). The precision level for nitrate+nitrite and phosphate was 0.02 and 0.003 µM, respectively. All analytical results were corrected against unamended blanks. The limit of detection (2 times the standard deviation of the blank) for the procedures was 0.075 µM for nitrate+nitrite and 0.008 µM for phosphate. For simplicity, we refer in the text to nitrate+nitrite as NO x and to ortho-phosphate as P. The quality assurance of the nutrient measurements was confirmed by the results of inter-comparison exercises (NOAA/NRC, JAPAN, QUASIMEME).
In Leeds, prior to the leaching of phosphate, 120 ml plastic containers were coated with iodine by adding a crystal or two of elemental iodine and placed in an oven at 40 • C for 10 min. The containers were then cooled to room temperature (22 • C) and swilled with Milli-Q water to remove all the excess I 2 . For the nitrate leaching, uncoated 250 ml plastic containers were used. For the phosphate leaching experiment, 50 ml of sterile seawater was added to each container in a biological safety cabinet. ∼6 mg of accurately weighed dust was then added and the containers were placed on a shaking table for 30 min, 2, 6, 24, and 49 h situated in a light proof box. The duration of the leaching experiment was based on the study of Mackey et al. (2012), which found an increasing amount of leached P after 24-48 h, and assumed stabilization. At each sampling, ∼7 ml were removed by syringe from the sample, filtered through 0.45 µm polycarbonate filters and stored at 4 • C for subsequent analysis. For the nitrate leaching experiment, 100 ml of sterile seawater were used and sampled in the same way as the other determinants, after 30 min of leaching. The analysis of all the nitrate and phosphate samples was carried out after the 48 h sample was collected. A series of 6 blank samples were run through the sampling procedure and the average values were subtracted from each sample analyzed (mean blank values were −3 ± 4 nmol P l −1 and 48 ± 50 nmol NO x l −1 ). Nutrient content was determined by standard SEAL AA-3 automated methods for phosphate (using the molybdate blue method), nitrate (as nitrite after Cd column reduction) and ammonium (using a fluorescence method). The precision of the replicate analysis was 1.8 ± 0.01 µM for phosphate, 6.0 ± 0.05 µM for nitrate and 5.75 ± 0.05 µM for ammonium.

Measurement of 33 PO 4 Uptake
Samples (10 ml) for PO 4 turnover times were collected every day and measured using 33 P-orthophosphate (Thingstad et al., 1993). Carrier-free 33 P-orthophosphate (PerkinElmer, specific activity: 370 MBq ml −1 ) was added to the samples at a final concentration of 20-79 pM. Samples used for the subtraction of the background and abiotic adsorption were fixed with 100% trichloroacetic acid (TCA) (final conc. 0.5%) before the isotope's addition. The samples were then incubated under subdued (laboratory) illumination. The incubation time varied between 1 and 20 min: short enough to assure a linear relationship between the fractions of the isotope adsorbed vs. the incubation time but long enough to reliably detect any isotope uptake above the background levels. Incubation was terminated by a cold chase of 100 mM KH 2 PO 4 (final conc. 1 mM). Subsamples (3.3 ml) were filtered in parallel onto 25 mm polycarbonate filters with 2, 0.6, and 0.2 µm pore sizes. All filters were placed on a Millipore 12 place manifold with Whatman (GF/C) glass fiber filters saturated with 100 mM KH 2 PO 4 as support. After filtration, the filters were placed in polyethylene vials with an Ultima Gold (Packard) scintillation cocktail and radio-assayed. After the radio-activities of the filter were corrected for those of the blank filter obtained from fixed samples, the phosphate turnover time (T [PO4] : h) was calculated as where f is the fraction (no dimension) of added isotope collected on the 0.2 µm filter after the incubation time (t: h).

Phosphate Concentrations in the Mesocosm Bags
Water samples were collected and analyzed immediately for their phosphate concentrations using the MAGIC method (Rimmelin and Moutin, 2005). The detection limit was 1.6 nM for phosphate.

Chlorophyll a
Seawater samples (500 ml) were passed through a Whatman GF/F filter (∼0.7 µm pore size) and extracted overnight (16 h) in 10 ml of 90% acetone solution in the dark (Holm-Hansen et al., 1965). Chlorophyll a concentrations were determined by the non-acidification method (Welschmeyer, 1994)

Picophytoplankton and Heterotrophic Bacterial Abundance
Samples for determining the picophytoplankton and heterotrophic bacterial abundance were collected every day throughout the experiment's duration. The samples were fixed with 0.2 µm filtered glutaraldehyde (a final concentration of 0.5%), kept at 4 • C for ∼45 min, flash-frozen in liquid nitrogen and then transferred to a −80 • C refrigerator until further processing. Frozen samples were thawed at room temperature and sub-samples were stained with SYBR Green I and incubated for 10 min in the dark, according to Vaulot and Marie (1999). Samples for picophytoplankton abundance were analyzed based on their auto-fluorescence signals, without pre-staining, using a FACSCalibur (Becton Dickinson) flow cytometer equipped with an air-cooled laser at 488 nm and a standard filter set-up. Flow cytometry data were acquired and processed with the Cell Quest Pro software (Becton Dickinson). An average estimated flow rate of 58 µL min −1 was used. The picophytoplankton carbon biomass was calculated from cell counts, assuming 175 fg C cell −1 for Synechococcus cells, 53 fg C cell −1 es (Campbell and Yentsch, 1989).

Primary Productivity (PP)
Photosynthetic carbon fixation rates were estimated using the 14 C incorporation method (Steemann-Nielsen, 1952). For each mesocosm, three light and one dark 320-mL polycarbonate bottles were filled with sample water during morning time (09:00-10:00 a.m.), inoculated with 5 µCi of NaH 14 CO 3 tracer (Perkin-Elmer) and incubated in a land-based tank for 3 h. At the end of the incubation, the spiked seawater samples were filtered through 0.2 µm polycarbonate filters (47 mm diameter) under low vacuum pressure (<150 mmHg) and the filters were collected for the determination of primary production rate. The filters were placed in 5-ml scintillation vials and were acidified with 1 ml of 0.1 N HCL in order to remove excess 14 C-bicarbonate overnight. After the addition of 4 ml scintillation cocktail (ULTIMA-GOLD), the radioactivity of the samples (disintegrations per minute, dpm) was measured in a Liquid Scintillation Counter (Packard Tri-Card 4000). Primary production rates were calculated by subtracting the dpm of the dark bottles from the respective light ones. We used a value of 26,400 mg C m −3 for the concentration of dissolved inorganic carbon (Triantaphyllou et al., 2010) and a value of 1.05 for the isotopic discrimination factor .

Bacterial Productivity (BP)
BP was estimated by the 3H-leucine method (Kirchmann et al., 1986), as modified by Smith and Azam (1992). For each mesocosm, duplicate SD, A and control samples were incubated with a mixture of L-[45 3H]-leucine (Perkin Elmer, 115 Ci mmol −1 ) and non-radioactive leucine to a final concentration of 20 nM. Samples were incubated for 2 h in the dark at the in-situ temperature, after which they were fixed and treated following the micro-centrifugation protocol (Smith and Azam, 1992), as described in detail by Van Wambeke et al. (2008). In brief, incubations were terminated after 2 h by the addition of trichloroacetic acid (TCA). Samples were then centrifuged at 16,000 × g and the resulting cell pellet was washed twice with 5% TCA and with 80% ethanol. Incorporation of 3H-leucine into the TCA-insoluble fraction was measured by liquid scintillation counting (Packard Tri-Carb 4000TR) after resuspension of the cell pellet in scintillation cocktail (Ultima-Gold). Bacterial production was calculated according to Kirchmann (1993), from 3H-leucine incorporation rates. Duplicate incubations had an analytical error <10%.
Concentration kinetics optimization was also performed to ensure that the bacterial growth was not limited by the concentration of leucine.
Dinitrogen (N 2 ) Fixation Rates 15 N 2 uptake measurements were performed using the 15 Nenriched seawater protocol described by Mohr et al. (2010), with minor modifications for the EMS (Rahav et al., 2013). The detailed methodology is presented in Rahav et al. (2016a) this Special Issue.

Dissolved Trace Metal Analysis
The samples used for the analysis of trace metals (each 1 L in volume) were collected in PTFE bottles (pre-treated overnight with 2N HNO 3 and rinsed afterwards with ultrapure Milli-Q water), and within 12 h from their collection the samples were filtered through pre-weighted nitrocellulose membrane filters (Millipore 0.45 µm pore size) under a clean laminar hood (class 100) in order to separate the particulate form from the dissolved form of the metals. The dissolved metals in the filtered samples were determined immediately after filtration through a process of pre-concentration accomplished by passing the sample through Chelex-100 resin columns for retaining the metals and by eluting them using 10 ml nitric acid 2N s.p. under the clean laminar hood (the preconcentration factor = 100). The pre-concentration procedure is a slight modification (Scoullos et al., 2007) of the procedure proposed by Riley and Taylor (1968) and Kingston et al. (1978). Trace metal concentrations (Cu, Pb, Zn, Mn, Al, Fe) in the eluates were determined by employing a Varian SpectrAA 200 Flame Atomic Absorption Spectrophotometer (FAAS) for Zn, Al and Fe and a Varian SpectrAA-640Z Graphite Furnace Atomic Absorption Spectrophotometer (GFAAS) with Zeeman background correction for Cu, Pb and Mn. The relative standard deviation (Sr = (S/χ) × 100) of the measurements that resulted from replicate (3-4) determinations and standard addition experiments was <5%.
to compare between the controls and the different treatments (P < 0.05) at the conclusion of the mesocosm experiment or between the results measured here and those compiled from the literature (Figures 5, 6). All tests were performed using the XLSTAT.

The Ambient Water and Initial Conditions
The characteristics of the ambient surface (10 m) water collected during the 8-9 of May north of the Island of Crete and the initial conditions of the experiment (1 day prior to the additions) were typical oligotrophic, representing the EMS offshore waters (Kress et al., 2014;Pitta et al., 2016) as well as other LNLC systems. All the measured state and rate parameters showed values typical of oligotrophic systems (Figure 3), which are consistent with previous observations of the spring conditions in the EMS (references in Figure 3). The chlorophyll a (chl-a) concentrations, primary production (PP) and bacterial production (BP) were at the lower range of the oligotrophic systems; 0.06 ± 0.00 µg l −1 , 0.42 ± 0.02 µg C l −1 h −1 , and 15.59 ± 7.57 ng C l −1 h −1 , respectively. Bacterial abundance (BA, 4.2 × 10 5 cells ml −1 ), Synechococcus abundance (Syn., 2.3 × 10 4 cells ml −1 ) and dinitrogen fixation (N 2 -fix, 0.21 ± 0.01 nmol N l −1 d −1 ) were within the oligotrophic range. The phosphate and phosphate turnover time were also characteristic to oligotrophic conditions, 0.012 µM and ∼1 h, respectively. The seawater temperature during the experiment was in the range of 19-20 • C, with no marked stratification of the water column inside the mesocosms.

The Chemical Composition and Nutrient Leachability of the Aerosols
The heavy metal concentrations were enriched in the SD compared to A, while the major element concentrations show similar values ( Table 2). The Si/Al ratios were at the upper range found in northern African mineral dusts, while the Fe/Ca ratios were relatively low (Formenti et al., 2011(Formenti et al., , 2014. The latter indicate a certain depletion in Al and Ca compared to typical African mineral dust. The Mn/Fe, Pb/Fe, Ni/Fe, and Cd/Fe ratios were similar to the equivalent ratios in the African aerosols sampled by Chien et al. (2016), while the Cu/Fe ratios in this study are somewhat higher. The Mn/Fe and P/Fe ratios (wt./wt.) in SD were similar to the Saharan soil used by Guieu et al. (2010) (0.015 vs. 0.015 for Mn/Fe and 0.014 vs. 0.017 for P/Fe) and to those measured by Chien et al. (2016) in African aerosols. The P/Fe, Pb/Fe, Cd/Fe, and Cu/Fe ratios in A were higher than in SD by 2.2, 1.7, 2.2, and 1.75, respectively, indicating enrichment in these elements probably due to a larger anthropogenic fraction. The amount of P leached was 2.1-2.7 and 2.2-3.7 nmol PO 4 per mg of SD and A, respectively (Table 3; Figure S1), representing an addition of ∼25-50% of the ambient concentrations. Similar P leached values were observed in other studies: 1.9 nmol phosphate per mg Saharan dust (Ridame and Guieu, 2002); 2.6-2.7 nmol phosphate per mg of Saharan dust (Chien et al., 2016); 4.2 nmol phosphate per mg of total suspended particles (TSP) sampled in the Red Sea (Mackey et al., 2012). The amount of nitrate leached was 20-26 and ∼54 nmol NO x per mg of SD and A, respectively (Table 3; Figure S1). The average leached N and P resulted in a distinctly lower N:P ratio in treatment SD (∼9:1) than in A (∼18:1) or 7:1 vs. 15:1, considering the larger leachability of P by Leeds laboratory ( Table 3). The higher leachable values of P obtained in Leeds compared to those retrieved in Haifa (∼30 and ∼75% for SD and A, respectively, Table 3; Figure S1) were probably attributed to the coating with iodine done in Leeds, which prevented adsorption into the walls of the containers. These leaching experiments were performed using sterile seawater, at different particle concentrations than those used in the mesocosm experiments. They also represent higher particle concentrations than those naturally found at the surface mixed layer after dust storms/aerosol depositions. Such differences, both in the seawater's biological reactivity and in the particle concentrations, may impact the amount of nutrient release (e.g., Ridame and Guieu, 2002;Mackey et al., 2012) and hence the leaching dynamics of N, P and other micronutrients. We therefore consider the experimental leaching amounts of N and P as an approximation of the total amount released in the mesocosms using the average values.
The turnover time of phosphorus (∼1 h), regardless of the treatment, represents extreme oligotrophic P starved conditions. The P released from the dust/aerosols was probably immediately consumed, as no increase in Tt was recorded even 3 h after addition.

Dissolved Trace Metal Concentrations
The presented dissolved trace metal concentrations were measured in the initial conditions (prior to additions) and also 3 and 24 h after addition. The amount of trace metals added to the experimental mesocosm after the SD and A additions is presented in Table 4. Fast enrichment (3 h after addition) was observed for Mn in both the SD and A treatments compared to the initial conditions and the control. A relatively fast enrichment of Ni, Fe, and Pb was observed in the SD treatment. Such a fast release of trace metals coincides with the observations in other studies Séguret et al., 2011). Generally, the dissolved concentrations of Mn, Ni, and Cu were similar to the values reported in Chien et al. (2016), while the Fe and Pb concentrations were enriched. The dissolved trace metal concentrations were similar to the range of values measured at the EMS offshore Israel (Figure 4).
Mn increased by ∼4 nM following the SD and ∼1 nM in the A treatments. This represents an increase by a factor of ∼2 and 1.3 in the ambient concentrations, respectively. Ni increased by ∼ 2 nM following the SD treatment. Pb and Fe increased by ∼1 and ∼6 nM in SD, respectively.

Biological Parameters
Both types of treatments triggered a positive change (relative to the unamended control mesocosms) in most of the performed rate and state measured parameters such as BA and BP, Syn abundance, chl-a, PP, and N 2 -fix. These changes are in agreement with other dust microcosm/mesocosm additions from the Mediterranean Sea and are presented as the maximal difference (treatment minus control, Figures 5A-E) or normalized to 1 mg l −1 of SD or A addition (Figures 6A,B). The added aerosol/dust concentrations (1 and 1.6 mg l −1 ) in our experiment fall within the linear dose-response range previously studied in the EMS (e.g., Herut et al., 2005; 0.2-4.9 mg l −1 ) and therefore supports a linear normalization. Overall, the maximal observed changes in the two treatments for the different parameters ranged from 25 Leaching experiments were carried in Haifa and Leeds as detailed in the Materials and Methods section and in Figure S1. to 660% per 1 mg l −1 of dust/aerosol addition ( Figure 6A). For SD treatment, the maximal change (increase) was observed after 48 h for BP (98%, 11.8 ± 1.2 ng C l −1 h −1 ), chl-a (27%, 33.6 ± 3.3 ng l −1 ) and N 2 -fix (223%, 0.21 ± 0.03 nmol N l −1 d −1 ) and 72 h for PP (61%, 0.28 ± 0.05 µg C l −1 h −1 ) and Syn. (23%, 2.75 ± 0.64 µg C l −1 ). BA showed a progressive increase (from 25%; 2.78 ± 0.43 µg C l −1 after 24 h to 48% at day 8) along the experiment. For A treatment, an immediate maximal change (increase) was observed after 3 h for BP (660%, 24.2 ± 8.5 ng C l −1 h −1 ), after 1 day for chl-a (47%, 33.2 ± 3.8 ng l −1 ), after 2 days for PP (102%, 0.29 ± 0.08 µg C l −1 h −1 ) and N 2 -fix (141%, 0.09 ± 0.03 nmol N l −1 d −1 ) and after 3 days for Syn.(41%, 2.89 ± 0.89 µg C l −1 ). BA showed a progressive increase (from 39%; 2.53 ± 0.46 ng C l −1 after 24 h to 63% at day 8) along the experiment.

DISCUSSION
Here we assess the impact of pure Saharan dust vs. mixed aerosols on the surface seawater autotrophic and heterotrophic microbial populations, mimicking the potential effects of an intense Saharan dust storm and a relatively intense mixed aerosol deposition. While the EMS is exposed to a relative high frequency of Saharan and other desert dust storms (Koçak et al., 2004;Ganor et al., 2010), which will probably increase in the future due to climate change , anthropogenic aerosols were found to be enriched in dust when it arrived at sampling sites after passing through populated and industrialized urban areas (Koçak et al., 2012). It has been observed that the anthropogenic component increased the trace metal content and changed their speciation in EMS aerosols (Kocak et al., 2010), and that such aerosols are enriched in nitrate compared to bioavailable phosphate (e.g., Chien et al., 2016). Therefore, any increase in anthropogenic aerosol deposition in the future and hence in the relative fraction of N deposition may influence the phytoplankton community structure in LNLC areas (Chien et al., 2016). Such a trend, of a relatively higher release of nitrate vs. phosphate, was observed here. While in the SD treatment more phosphate (by ∼34%) was added to the mesocosms compared to the A treatment, the opposite trend was observed for nitrate. The A treatment added more nitrate (by ∼40%) than the SD treatment. Consequently, the nitrate/phosphate ratios added in the A treatment were twice than those added in the SD treatment, 18 vs. 9 using the average or 15 vs. 7 considering a larger release of P ( Table 3). The phosphate turnover time prior to the additions of dust/aerosol was ∼1 h, a typical value of P deficiency in the EMS (Zohary and Robarts, 1998;Flaten et al., 2005;Tanaka et al., 2011). The SD and A additions did not trigger an increase in the phosphate turnover time caused by the release of phosphate from the dust/aerosol. By definition, the Tt equals the amount of bioavailable phosphate divided by its consumption (or uptake) rate. This suggests that the amount of phosphate added to the system from the dust/aerosol through leaching was relatively small and that it was rapidly removed by the microbial community. Bioassays in which a significant amount of phosphate was released showed an increase in phosphate Tt Tanaka et al., 2011). The amount of leachable phosphate and nitrate in the SD treatment, 3.8 and ∼23 nM, respectively, corresponds to a potential increase of ∼58 ng chl-a l −1 by using P:C ratios or ∼55 ng chl a l −1 by using N:C ratios of Redfield (1:16:106) and C:chl a ≈ 80 (wt., Behrenfeld et al., 2005). The observed chl-a increase is somewhat lower, ∼30 ng chl-a l −1 . However, the highest bacterial production rate was ∼3.3 nmol C l −1 h −1 , at 24 h after addition, which corresponds to ∼1.6 nmol P l −1 d −1 , using a C/P ratio of 50 for bacteria (Fagerbakke et al., 1996), although the latter is not constant (Zimmerman et al., 2014;Godwin and Cotner, 2015). Assuming that out of the released phosphate ∼ 1.6 nM was used by the bacteria and 2.2 nM was used by the autotrophs, the latter corresponds to the observed increase in chl-a (∼30 ng l −1 ). The lower observed chl-a values may also infer some grazing activity. A similar calculation for the A treatment shows a similar result, which corresponds to the observed increase in chl-a.
Additional trace metals, micronutrients (Fe, Zn, Mn, Co-not measured) or potentially toxic metals that have possibly been delivered by the dust were not assessed via leaching experiments in this study. However, these were measured 3 h and 1 day after the additions (Table 4), showing values similar to an open remote EMS station (100 km offshore Israel), to the Red Sea (Herut et al., unpublished) or, for part of the elements, to the mean ocean concentration, although the Mn, Fe and Zn values were significantly higher (Statham and Hart, 2005;Semeniuk et al., 2009;Moore et al., 2013).
The SD treatments were relatively enriched (p < 0.05) in dissolved Mn, Ni, and Fe concentrations compared to the A mesocosms, probably due to the larger mass addition (1.6 mg l −1 vs. 1 mg l −1 ). All the measured dissolved metal concentrations were well below the ecological ambient water quality criteria for chronic levels (Buchman, 2008). Metal bioavailability may determine phytoplankton productivity, especially by iron, copper, cobalt, zinc and nickel, via scarcity or toxic effects (Huertas et al., 2014). Their specific impact, including manganese, on biological processes is mainly related to the metalloproteins state (Cvetkovic et al., 2010). The enrichment of Mn in SD and A treatments as compared to C (p < 0.05) may serve as a micronutrient supporting photosynthesis (Huertas et al., 2014) and/or may reverse toxicity effects of Cu and Zn (Sunda, 1987). Although the bioavailability of the measured dissolved Fe is not certain (Rue and Bruland, 1995), the amount is high enough to exclude possible Fe limitation (e.g., Statham and Hart, 2005). Assuming that particulate Fe dissolution in seawater is ∼0.06% (0.03-0.17%, Blain et al., 2004;Chien et al., 2016), a maximum amount of ∼50 nmol Fe l −1 (900 nmol particulate Fe l −1 × 0.06%) was released from  Rahav et al., unpublished and this study (E). The red lines represents the observed values measured for SD and A during this study. The asterisks above the columns represent statistically significant differences (one-way ANOVA and a Fisher LSD means comparison test, P < 0.05) for mean values of P additions between stations. the SD addition, larger than in the A treatment (∼24 nmol Fe l −1 ), and probably supplying enough bioavailable Fe 2+ . The significantly higher (p = 0.02) N 2 -fix rates measured in the SD treatment (Rahav et al., 2016b) may be related to a relatively larger release of Fe in this treatment and its relatively lower nitrate/phosphate ratio (∼10:1). In addition, the significantly (p < 0.05) higher Ni concentrations in SD as compared to A or C, is known to play an important role in the cellular physiology of diazotrophic cyanobacteria such as Trichodesmium (Rodriguez and Ho, 2014) and supports the latter observation and the appearance of this N 2 -fixer in the SD treatment (Rahav et al., 2016b).  Ternon et al., 2011, C,D). The black line represents a 1:1 ratio between A and SD values. The asterisks above the columns represent statistically significant differences (one-way ANOVA and a Fisher LSD means comparison test, P < 0.05).
Both types of dry aerosol additions (SD and A) triggered a positive change (25-600% per 1 mg/L addition) in all rate and state measurements (Figures 5, 6), showing a general sequence (p < 0.05) of higher relative change in the rate parameters: N 2 -fix > BP > PP > Chla ≈ BA ≈ Syn for SD and BP > PP ≈ N 2 -fix > Chla ≈ BA ≈ Syn for A. Nevertheless, a larger normalized increase was observed in treatment A for all parameters, except for N 2 -fix ( Figure 6B). In addition, the BP dramatically increased within 3 h after addition (∼600%, yet reached 154% after 24 h), while no such response was observed for SD. The primary differences between the two treatments is the leachable molar ratio of nitrate/phosphate; almost twice in treatment A (∼15:1) than SD (∼9:1), and the additional amount of added leachable N by treatment A (∼30% more than SD) ( Table 3). No significant changes were observed after additions and between treatments in the dissolved organic nutrients, not likely to dominate the observed change in the heterotrophic bacterial activity. These differences in the relative impact (normalized to mass), between pure mineral/desert dust and mixed aerosols, were also observed by Ternon et al. (2011) for the addition of a dust analog (manipulated Saharan soil) and aerosols (total suspended material in air) collected by low volume sampler at sea (Figures 6C,D). While in the latter on-board experiment lower quantities of P were release by the aerosols (0.3-1.6 nmol P l −1 ) as compared to the Saharan analog (6 nmol P l −1 ), in our experiment similar amounts of P were added by A and SD, but larger quantities of N were added in A (Table 3). Yet, only the inorganic phase of nutrients was considered here, whereas the dissolved organic pools likely to be present in aerosols (e.g., Markaki et al., 2010) were not measured and should be further considered. Apart from the increased solubility in aerosols, the smaller size particles as compared to crustal mineral aerosols, increases the particle's surface/volume ratio and solubility (Baker and Jickells, 2006), implying a larger impact per mass of a mixed character as A treatment. Our data suggest that despite the co-limitation of P and N (Kress et al., 2005;Zohary et al., 2005), the additional N released by the A treatment may have triggered the relative larger response in most of the rate and state parameters. A similar observation was reported for the DUNE experiment in the Western Mediterranean Sea (Ridame et al., 2014) and was calculated for the data presented during the BOUM campaign across the Mediterranean Sea in which aerosol addition triggered stronger responses of the microbial community than Saharan dust (Tanaka et al., 2011;Ternon et al., 2011). Yet, the magnitude of the maximal net effect imposed by SD or A on the surface microbial populations (% change normalized to 1 mg l −1 ) is significantly lower than that reported from a microcosm experiment in the Ionian Sea during the BOUM campaign in summer 2008 (by ∼1 order of magnitude, Figure 6C, Ternon et al., 2011). When normalizing the responses from the BOUM study, Saharan analog or aerosols addition resulted in a general sequence of PP ≈ N 2 -fix ≈ Syn. > BA. Furthermore, the overall responses triggered by the aerosol addition in the BOUM campaign (microbial biomass and activity) were much higher than that reported for the Saharan analog amendments. This results in a high calculated aerosol/Saharan analog ratios in the BOUM study, ranging from ∼50 to 100% ( Figure 6D). These different responses by different aerosols should be further studied in detail considering all bioavailable macro and micro nutrients or other inhibiting elements.
The role atmospheric deposition plays in influencing bacterioplankton dynamics in LNLC regions is important for understanding (observationally and via ocean biogeochemical models) the current and future functioning of LNLC regions (Guieu et al., 2014a). Increasing temperatures leading to increased dust emissions in the EMS (Klingmüller et al., 2016) and anthropogenic activities contributing to the Mediterranean atmospheric chemical composition and aerosols solubility (exposure to lower pH levels during transport), are expected to continue in the future due to climate change, and thereby supply more macro and micro nutrients into its surface oligotrophic water. A warmer atmosphere may influence the intensity and length of the already well stratified water column in the EMS and hence the impact of the atmosphere as a significant external source of new macro and micro nutrients. We show here that the above trends may lead to the alteration of the heterotrophic/autotrophic relationships.

ACKNOWLEDGMENTS
This work was financed by the European Union Seventh Framework Program (FP7/2007(FP7/ -2013 under grant agreement no. 228224, "MESOAQUA: Network of leading MESOcosm facilities to advance the studies of future AQUAtic ecosystems from the Arctic to the Mediterranean" through grants to BH, IBF, ER, ZS, and MDK. The authors wish to thank G. Piperakis for his technical assistance, A. Konstantinopoulou for assistance with bacterial production analyses, D. Podaras and S. Diliberto for assistance during the experiment and N. Sekeris for his help with constructions and ideas on technical solutions. The captain and the crew of the R/V Philia are also thanked for their assistance during the transportation of water from the sea to the CRETACOSMOS facility. Funding was also provided by the Israel Science Foundation grants (996/08) to IBF and BH. Funding was also provided by Leverhulme Trust (grant RPG 406) to MDK.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmars. 2016.00226/full#supplementary-material Figure S1 | Leached nitrate (A) and phosphate (B) from the SD (circle) and A (triangle) particles as retrieved in the Leeds (white) and Haifa (red) laboratories.