# Impact of dust addition on Mediterranean plankton communities under present and future conditions of pH and temperature: an experimental overview

Frédéric Gazeau, Céline Ridame, France van Wambeke, Samir Alliouane, Christian Stolpe, Jean-Olivier Irisson, Sophie Marro, Jean-Michel Grisoni, Guillaume de Liège, Sandra Nunige, et al.

► To cite this version:

Frédéric Gazeau, Céline Ridame, France van Wambeke, Samir Alliouane, Christian Stolpe, et al. Impact of dust addition on Mediterranean plankton communities under present and future conditions of pH and temperature: an experimental overview. Biogeosciences, 2021, 18 (17), pp.5011-5034. 10.5194/bg-18-5011-2021. hal-03342476

HAL Id: hal-03342476

https://amu.hal.science/hal-03342476

Submitted on 16 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Impact of dust addition on Mediterranean plankton communities under present and future conditions of pH and temperature: an experimental overview

Frédéric Gazeau<fnref n="1" />, Céline Ridame<fnref n="2" />, France Van Wambeke<fnref n="3" />, Samir Alliouane<fnref n="1" />, Christian Stolpe<fnref n="1" />, Jean-Olivier Irisson<fnref n="1" />, Sophie Marro<fnref n="1" />, Jean-Michel Grisoni<fnref n="4" />, Guillaume De Liège<fnref n="4" />, Sandra Nunige<fnref n="3" />, Kahina Djaoudi<fnref n="3" />, Elvira Pulido-Villena<fnref n="3" />, Julie Dinasquet<fnref n="5" />,<fnref n="6" />, Ingrid Obernosterer<fnref n="6" />, Philippe Catala<fnref n="6" />, and Cécile Guieu<fnref n="1" />

**Correspondence:** Frédéric Gazeau (frederic.gazeau@imev-mer.fr)

Received: 4 June 2020 – Discussion started: 24 August 2020

Revised: 7 July 2021 – Accepted: 11 July 2021 – Published: 13 September 2021

**Abstract.** In low-nutrient low-chlorophyll areas, such as the Mediterranean Sea, atmospheric fluxes represent a considerable external source of nutrients likely supporting primary production, especially during periods of stratification. These areas are expected to expand in the future due to lower nutrient supply from sub-surface waters caused by climate-driven enhanced stratification, likely further increasing the role of atmospheric deposition as a source of new nutrients to surface waters. Whether plankton communities will react differently to dust deposition in a warmer and acidified environment remains; however, an open question. The potential impact of dust deposition both in present and future climate conditions was investigated in three perturbation experiments in the open Mediterranean Sea. Climate reactors (300 L) were filled with surface water collected in the Tyrrhenian Sea, Ionian Sea and in the Algerian basin during a cruise conducted in the frame of the PEACETIME project in May–June 2017. The experiments comprised two unmodified control tanks, two tanks enriched with a Saharan dust analogue and two tanks enriched with the dust analogue and maintained under warmer (+3 °C) and acidified (–0.3 pH unit) conditions. Samples for the analysis of an extensive number of biogeochemical parameters and processes were taken over

the duration (3–4 d) of the experiments. Dust addition led to a rapid release of nitrate and phosphate, however, nitrate inputs were much higher than phosphate. Our results showed that the impacts of Saharan dust deposition in three different basins of the open northwestern Mediterranean Sea are at least as strong as those observed previously, all performed in coastal waters. The effects of dust deposition on biological stocks were different for the three investigated stations and could not be attributed to differences in their degree of oligotrophy but rather to the initial metabolic state of the community. Ocean acidification and warming did not drastically modify the composition of the autotrophic assemblage, with all groups positively impacted by warming and acidification. Although autotrophic biomass was more positively impacted than heterotrophic biomass under future environmental conditions, a stronger impact of warming and acidification on mineralization processes suggests a decreased capacity of Mediterranean surface plankton communities to sequester atmospheric CO₂ following the deposition of atmospheric particles.

<footnote n="1">¹Sorbonne Université, CNRS, Laboratoire d'Océanographie de Villefranche, LOV, 06230 Villefranche-sur-Mer, France</footnote>

<footnote n="2">²CNRS-INSU/IRD/MNHN/UPMC, LOCEAN: Laboratoire d'Océanographie et du Climat: Expérimentation et Approches Numériques, UMR 7159, 75252 Paris CEDEX 05, France</footnote>

<footnote n="3">³Aix-Marseille Université, Université de Toulon, CNRS/INSU, IRD, MIO, UM 110, 13288 Marseille, France</footnote>

<footnote n="4">⁴Sorbonne Université, CNRS, Institut de la Mer de Villefranche, IMEV, 06230 Villefranche-sur-Mer, France</footnote>

<footnote n="5">⁵Scripps Institution of Oceanography, University of California San Diego, USA</footnote>

<footnote n="6">⁶CNRS, Sorbonne Université, Laboratoire d'Océanographie Microbienne, LOMIC, 66650 Banyuls-sur-Mer, France</footnote>
# 1 Introduction

Atmospheric deposition is well recognized as a significant source of micro- and macro-nutrients for surface waters of the global ocean (Duce et al., 1991; Jickells et al., 2005; Moore et al., 2013). The potential modulation of the biological carbon pump efficiency and the associated export of carbon by atmospheric deposition events are still poorly understood and quantified (Law et al., 2013). This is especially true for low-nutrient low-chlorophyll (LNLC) areas where atmospheric fluxes can play a considerable role in nutrient cycling and that represent 60% of the global ocean surface area (Longhurst et al., 1995) as well as 50% of global carbon export (Emerson et al., 1997). These regions are characterized by low availability of macronutrients (N, P) and/or micronutrients (trace metals, in particular Fe) that can severely limit or co-limit phytoplankton growth during long periods of the year.

The Mediterranean Sea is a typical example of these LNLC regions with overall surface chlorophyll $a$ concentrations below $0.2 \mu g L^{-1}$ all year round, except in the Ligurian Sea where relatively large blooms can be observed in late winter–early spring (Mayot et al., 2016). Recent estimates indicate that the atmospheric input of nutrients in the Mediterranean Sea is within the same order of magnitude as riverine inputs (Powley et al., 2017), and, therefore, a considerable external source of nutrients (Richon et al., 2018). Atmospheric deposition originates both from natural (mainly Saharan dust) and anthropogenic sources (e.g. Bergametti et al., 1989; Desboeufs et al., 2018). Dust deposition, mostly in the form of pulsed inputs, is mainly associated with wet deposition (Loÿe-Pilot and Martin, 1996). Ternon et al. (2010) reported an average annual dust flux over 4 years of $11.4 \ g m^{-2} yr^{-1}$ (average during the period 2003–2007) at the DYFAMED station in the northwestern Mediterranean Sea. In this region, the most important events reported in the 2010 decade amounted to ~$22 \ g m^{-2}$ (Bonnet and Guieu, 2006; Guieu et al., 2010b).

Atmospheric deposition provides new nutrients to surface waters (Guieu et al., 2010b; Kouvarakis et al., 2001; Markaki et al., 2003; Ridame and Guieu, 2002), Fe (Bonnet and Guieu, 2006) and other trace metals (Desboeufs et al., 2018; Guieu et al., 2010b; Theodosi et al., 2010), representing significant inputs likely supporting primary production in particular during the period of stratification in spring/summer (Bonnet et al., 2005; Ridame and Guieu, 2002), although no direct correlation between dust and ocean colour could be found from long series of satellite observation in that part of the Mediterranean Basin (Guieu and Ridame, 2020).

Previous micro- and mesocosm experiments have shown that wet dust deposition events in the northwestern Mediterranean Sea (the dominant deposition mode in that basin) are a stronger source of bioavailable nutrients compared to dry deposition. Wet deposition provides both new N and P while dry deposition supplies primarily P and, in contrast to wet

deposition, does not stimulate the growth of the autotrophic community with the exception of diazotrophs (Ridame et al., 2013), resulting in no significant increase in chlorophyll $a$ concentrations and primary production (Guieu et al., 2014a). In addition, wet dust deposition also modifies the bacterial assemblage leading to even stronger enhancements of heterotrophic production and respiration rates (Pulido-Villena et al., 2014). The carbon budget established from four artificial seeding experiments during the DUNE project (Guieu et al., 2014a) showed that by stimulating predominantly heterotrophic bacteria, atmospheric wet dust deposition can enhance the heterotrophic behaviour of these oligotrophic waters. This has the potential to reduce organic carbon export to deep waters during the winter mixing period (Pulido-Villena et al., 2008) and ultimately limit net atmospheric $CO_2$ draw-down.

Conversely, the deposition of lithogenic particle from Saharan dust can promote aggregation and ballast organic matter leading to enhanced vertical export of organic carbon (Bressac et al., 2014; Desboeufs et al., 2014; Louis et al., 2017a; Ternon et al., 2010). These lithogenic processes can represent a major part of the carbon export following a dust deposition event (up to 50% during the DUNE experiment; Bressac et al., 2014). Recently, Louis et al. (2017a) showed that Saharan dust deposition can also trigger the abiotic formation of transparent exopolymeric particles (TEPs), leading to the formation of organic-mineral aggregates, a formation process that is highly dependent on the quality and quantity of TEP precursors initially present in seawater.

In response to ocean warming and increased stratification, nutrient cycling in the open ocean is being and will continue to be perturbed in the next decades with regionally variable impacts (IPCC, 2019). Overall, LNLC areas are expected to expand in the future (Irwin and Oliver, 2009; Polovina et al., 2008) due to thermal-stratification-related reduction of the nutrient supply from sub-surface waters (Behrenfeld et al., 2006). As such, the role of atmospheric deposition as a source of new nutrients to surface waters might increase. Ongoing warming and acidification (IPCC, 2019) are also evidenced in the Mediterranean Sea (e.g. Kapsenberg et al., 2017; The Mermex group, 2011). Whether or not plankton communities will respond differently to dust deposition in future conditions is still largely unknown. Although dependent on resource availability, it is well known that remineralization by bacteria is subject to positive temperature control (López-Urrutia and Morán, 2007). Given that warming has no effect on primary productivity when plankton communities are nutrient limited (Marañón et al., 2018), temperature increase will most likely further push the balance towards net heterotrophy in oligotrophic areas.

In contrast, an in situ mesocosm experiment conducted during the summer stratification period in the northwestern Mediterranean Sea showed that the plankton community was not sensitive to ocean acidification under strong nutrient limitation (Maugendre et al., 2017, and references therein). A
batch experiment (Maugendre et al., 2015) showed that, under nutrient-depleted conditions in late winter, ocean acidification has a very limited impact on the plankton community and that small species (e.g. Cyanobacteria) might benefit from warming with a potential decrease in the export and energy transfer to higher trophic levels. In contrast, in more eutrophic (coastal) conditions, Sala et al. (2016) showed that ocean acidification had a positive effect on phytoplankton, especially on the pico and nano size classes. Similarly, Neale et al. (2014) showed that ocean acidification could lead to enhanced chlorophyll levels under low-light conditions with an opposite effect under high irradiance, in coastal communities of the Alboran Sea.

To date and to the best of our knowledge, there have been no attempts to evaluate the behaviour of plankton communities influenced by atmospheric deposition in the context of future temperature and pH changes. Such experiments were, therefore, conducted in the framework of the PEACETIME project (ProcEss studies at the Air-sEa Interface after dust deposition in the MEDITERRANEAN sea; http://peacetime-project.org/, last access: 12 August 2021) on board the R/V *Pourquoi Pas*? during May-June 2017. The project aimed at studying and parameterizing the chain of processes occurring in the Mediterranean Sea driven by atmospheric deposition events including under ongoing environmental changes (Guieu et al., 2020a). During the cruise, three perturbation experiments were conducted in climate reactors (300 L tanks) filled with surface water collected in the Tyrrhenian Sea (TYR), Ionian Sea (ION) and Algerian basin (FAST; Fig. 1). Six tanks were used to follow, simultaneously and with a high temporal resolution, the evolution of biological activity and stocks, nutrients, and dissolved organic matter as well as particle dynamics and export, following a dust deposition event simulated both under present environmental conditions and under a realistic climate change scenario for 2100 (ca. +3°C and -0.3 pH units; IPCC, 2013). In this paper, we will present the general set-up of the experiments and the evolution of nutrient and plankton communities (heterotrophic and autotrophic prokaryotes, photosynthetic eukaryotes, and micro-and meso-zooplankton). Other papers, related to these experiments in this special issue, focus on plankton metabolism (primary production, heterotrophic prokaryote production) and carbon export (Gazeau et al., 2021), microbial food web (Dinasquet et al., 2021), nitrogen fixation (Ridame et al., 2021), and on the release of insoluble elements (Fe, Al, REE, Th, Pa) from dust (Roy-Barman et al., 2021).

## 2 Material and methods

### 2.1 General set-up

Six experimental tanks (300 L; Fig. 2), in which the irradiance spectrum and intensity can be finely controlled and future ocean acidification and warming conditions can be fully

reproduced, were installed in a temperature-controlled container. The tanks are made of trace-metal free high-density polyethylene (HDPE) with a height of 1.09 m, a diameter of 0.68 m, a surface area of 0.36 m² and a volume of 0.28 m³. Each tank was equipped with a lid containing six rows of LEDs (Alpheus©). Each of these rows were composed of blue, green, cyan and white units in order to mimic the natural sun spectrum. At the conical base of each tank, a polyethylene (PE) bottle was screwed onto a polyvinyl chloride (PVC) valve that remained open during the duration of the whole experiment to collect the sinking material. Photosynthetically active radiation (PAR; 400–700 nm) and temperature were continuously monitored in each tank using respectively QSL-2100 Scalar PAR Irradiance Sensors (Biospherical Instruments©) and pt1000 temperature sensors (Metrohm©) connected to a D230 datalogger (Consort©).

Prior to the start of the experiments, tanks were cleaned following the protocol described by Bressac and Guieu (2013). Three sets of experiments were carried out at the long-duration stations ION, TYR and FAST, respectively, and comprised two unmodified control tanks (C1 and C2), two tanks enriched with Saharan dust (D1 and D2), and two tanks enriched with Saharan dust and maintained under warmer (+3 °C) and acidified (–0.3 pH unit) conditions (G1 and G2). The atmosphere above tanks C1, C2, D1 and D2 was flushed with ambient air (ca. 400 ppm, 6 L min⁻¹) and tanks G1 and G2 were flushed with air enriched with CO₂ (ca. 1000 ppm, 6 L min⁻¹) in order to prevent CO₂ degassing from the acidified tanks. CO₂ partial pressure (pCO₂) in both ambient air and CO₂-enriched air was monitored using two gas analysers (LI-820, LICOR©). The CO₂ concentration in the CO₂-enriched air was manually controlled through small injections of pure CO₂ (Air Liquide©) using a mass flow controller. Mixing in the tanks was ensured by a rotative PVC blade (9 rpm) mimicking natural turbulence.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Location of the sampling stations in the Mediterranean Sea on board the R/V *Pourquoi Pas*? during the PEACETIME cruise. Background shows satellite-derived surface chlorophyll *a* concentration averaged over the entire duration of the cruise (courtesy of Louise Rousselet).</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Diagram of an experimental tank (climate reactor).</figcaption></figure>

The tanks were filled by means of a peristaltic pump (Verder© VF40 with EPDM hose, flow of 1200 L h⁻¹) collecting seawater below the base of the boat at around 5 m, used to supply surface seawater continuously to a series of instruments during the entire campaign. In order to homogeneously fill the tanks, the flow was divided into six HDPE pipes distributing the water simultaneously into the different tanks. The procedure was started at the end of the day at all three stations and took approximately 2 h (including rinsing and initial sampling). While filling the tanks, samples were taken for the measurements of selected parameters (sampling time = t-12h before dust seeding; Table 1). After filling the tanks, seawater in tanks G1 and G2 was slowly warmed overnight using 500 W heaters, controlled by temperature-regulation units (COREMA©), to reach an offset of +3 °C. ¹³C-bicarbonate was added to all tanks at 04:00 LT (local time; Gazeau et al., 2021), and at 04:30 LT G1 and G2 were acidified by addition of CO₂-saturated filtered (0.2 µm) seawater (~ 1.5 L in 300 L; collected when filling the tanks at each station) to reach a pH offset of -0.3. Further samples for a range of parameters were taken (sampling time = t0; Table 1), followed by dust seeding carried out between 07:00 and 09:00 LT in tanks D1, D2, G1 and G2. The same dust analogue flux was applied as in the DUNE 2009 experiments described in Desboeufs et al. (2014). The dust was derived from the < 20 µm fraction of soil collected

in southern Tunisia (a major source for material transported and deposited in the northwestern Mediterranean) consisting of quartz (40%), calcite (30%) and clay (25%) with most particles (99%) smaller than 0.1 µm (Desboeufs et al., 2014). The collected material underwent an artificial chemical ageing process by addition of nitric and sulfuric acid (HNO₃ and H₂SO₄, respectively) to mimic cloud processes during atmospheric transport of aerosol with anthropogenic acid gases (Guieu et al., 2010a, and references therein). To mimic a realistic wet flux event for the Mediterranean of 10 g m⁻², 3.6 g of this analogue dust was quickly diluted in 2 L ultra-high-purity water (UHP water; 18.2 MΩ cm⁻¹ resistivity) and sprayed at the surface of the tanks using an all-plastic garden sprayer (duration = 30 min). The total N and P mass in the dust were 1.36±0.09% and 0.055±0.003%, respectively (see Desboeufs et al., 2014, for a full description of dust chemical composition). Biogeochemical parameters and processes measured during the experiments are listed in Table 1. The experiment lasted 3 d (72 h) at stations TYR and ION and 4 d (96 h) at station FAST, as constrained by the time available between stations. Seawater sampling was conducted 1 h (t1h), 6 h (t6h), 12 h (t12h), 24 h (t24h), 48 h (t48h) and 72 h (t72h) after dust additions in all three experiments with an additional sample after 96 h (t96h) at FAST. Acid-washed silicone tubes were used for transferring the water
collected by gravity from the tanks to the different vials or containers.

## 2.2 Analytical methods

### 2.2.1 Carbonate chemistry

Seawater samples for pH measurements were stored in 300 mL glass bottles with a glass stopper, pending analysis on board (within 2 h). Samples were transferred to 30 mL quartz cells, and absorbances at 434, 578 and 730 nm were measured at 25 °C on a Cary60 UV spectrophotometer (Agilent®) before and after addition of 50 µL of purified metacresol purple provided by Robert H. Byrne (University of South Florida, USA) following the method described by Dickson et al. (2007). pH on the total scale (pH~T~) was computed using the formula and constants of Liu et al. (2011). The accuracy of pH measurements (0.007 pH units) was estimated using a TRIS buffer solution (salinity 35, provided by Andrew Dickson, Scripps Institution of Oceanography, USA).

Seawater samples (500 mL) for total alkalinity ($A_T$) measurements were filtered on GF/F membranes and analysed onboard within one day. $A_T$ was determined potentiometrically using a Metrohm© titrator (Titrando 888) and a glass electrode (Metrohm©, Ecotrode Plus) calibrated using first NBS buffers (pH 4.0 and pH 7.0, to check that the slope was Nernstian) and then using a TRIS buffer solution (salinity 35, provided by Andrew Dickson, Scripps Institution of Oceanography, USA). Triplicate titrations were performed on 50 mL sub-samples at 25 °C, and $A_T$ was calculated as described by Dickson et al. (2007). Titrations of standard seawater provided by Andrew Dickson (Scripps Institution of Oceanography, USA; batch 151) yielded $A_T$ values within 5 µmol kg⁻¹ of the nominal value and a standard deviation of 1.5 µmol kg⁻¹ ($n = 40$).

All parameters of the carbonate chemistry were determined from pH~T~, $A_T$, temperature, salinity, and phosphate and silicate concentrations using the R package seacarb. Propagation of errors on computed parameters was performed using the new function "error" of this package, encompassing errors associated with the estimation of $A_T$ and pH~T~ as well as errors on the dissociation constants (Orr et al., 2018).

### 2.2.2 Nutrients

Seawater samples for dissolved nutrients were collected in polyethylene bottles after passage through sterile membrane filter capsules (Sartobran® 300; 0.2 µm) connected to the sampling tubes of each tank (Sartobran® 300; 0.2 µm) and analysed directly on board. Nitrate + nitrite ($NO_x$) and silicate ($Si(OH)_4$) measurements were conducted using a segmented flow analyser (AAIII HR Seal Analytical®) according to Aminot and Kérouel (2007) with a detection limit

of 0.05 µmol L⁻¹ for $NO_x$ and 0.08 µmol L⁻¹ for $Si(OH)_4$. In addition, at t-12h, $NO_x$ was also analysed by spectrometry at 540 nm, with a 1 m liquid waveguide capillary cell (LWCC), with a detection limit of ~10 nmol L⁻¹, and the reproducibility was ~6%. Ammonium concentrations in samples from t-12h were also measured on board using a Fluorimeter TD-700 (Turner Designs®) according to Holmes et al. (1999). This later method is based on the reaction of ammonia with orthophthaldehyde and sulfite and has a detection limit of 0.01 µmol L⁻¹. Dissolved inorganic phosphorus (DIP) concentrations were quantified using the LWCC method according to Pulido-Villena et al. (2010). The LWCC was 2.5 m long and the detection limit was 1 nmol L⁻¹.

### 2.2.3 Pigments

For pigment analysis, 2.5 L seawater from the tanks was filtered onto GF/F filters, immediately frozen in liquid nitrogen, and stored at -80 °C pending analysis at the SAPIGH analytical platform at the Institut de la Mer de Villefranche (IMEV, France). Filters were sonicated at -20 °C in 3 mL methanol (100 %) containing an internal standard (vitamin E acetate, Sigma®) and clarified one hour later by vacuum filtration through GF/F filters. The extracts were rapidly analysed (within 24 h) on a complete Agilent® Technologies 1200 series HPLC system. The pigments were separated and quantified as described in Ras et al. (2008).

### 2.2.4 Flow cytometry

For flow cytometry, samples (4.5 mL) were fixed with glutaraldehyde grade I (1 % final concentration) and incubated for 30 min at 4 °C, quick-frozen in liquid nitrogen and stored at -80 °C until analysis. Samples were thawed at room temperature. Counts were performed on a FACSCanto II flow cytometer (Becton Dickinson®) equipped with three air-cooled lasers: blue (argon 488 nm), red (633 nm) and violet (407 nm). Following Marie et al. (2010), *Synechococcus* spp. were distinguished by their strong orange fluorescence (585 ± 21 nm), and autotrophic pico- and nanoeukaryotes were distinguished by their scatter signals of red fluorescence (> 670 nm). For the enumeration of heterotrophic prokaryotes, cells were stained with SYBR Green I (Invitrogen – Molecular Probes) at 0.025 % ($v/v$) final concentration for 15 min at room temperature in the dark. Stained prokaryotic cells were distinguished and enumerated according to their right-angle light scatter (SSC) and green fluorescence at 530/30 nm. Heterotrophic prokaryotes were distinguished from autotrophic prokaryotes based on the green vs. red fluorescent signal. The same procedure was used for the enumeration of HNF, after staining with 0.05% ($v/v$) final SYBR Green I concentration for 15–30 min at room temperature in the dark (Christaki et al., 2011). Fluorescent beads (1.002 µm; Polysciences Europe®) were systematically added to all samples as internal standard. Cell
<figure><img src="image_5.png" /><figcaption>**Table 1.** List of parameters and processes investigated during the three experiments at stations TYR, ION and FAST. Corresponding papers are indicated. pH<sub>T</sub>: pH on the total scale, A<sub>T</sub>: total alkalinity, <sup>13</sup>C-C<sub>T</sub>: <sup>13</sup>C signature of dissolved inorganic carbon, NO<sub>x</sub>: nitrate + nitrite, DIP: dissolved inorganic phosphorus, Si(OH)<sub>4</sub>: silicate, DFe: dissolved iron, DAl: dissolved aluminium, Th-REE-Pa: thorium (<sup>230</sup>Th and <sup>232</sup>Th), rare earth elements and protactinium (<sup>231</sup>Pa), POC: particulate organic carbon, DOC: dissolved organic carbon, <sup>13</sup>C-DOC: <sup>13</sup>C signature of dissolved organic carbon, TEP: transparent exopolymer particles, NCP/CR: net community production and community respiration (oxygen based), <sup>14</sup>C-PP: primary production based on <sup>14</sup>C incorporation.</figcaption></figure>
concentrations were determined based on counts and flow rate, estimated with TruCount beads (BD biosciences©). The biomass of each group was estimated based on conversion equations and/or factors found in the literature (see Sect. 2.3.2).

### 2.2.5 Micro-phytoplankton and micro-heterotrophs

At t-12h, 500 mL samples were collected in glass vials and immediately preserved with 5% final concentration acidic Lugol's solution. Back at the Laboratoire d'Océanographie de Villefranche (LOV, France), 100 mL aliquots were transferred to sedimentation chambers (Utermohl) and counted under an inverted microscope at ×200 to ×400 magnification.

### 2.2.6 Meso-zooplankton

At the end of each experiment, the sedimentation bottles were removed, fixed with formaldehyde 4% (see Gazeau et al., 2021) and stored for analysis back at the Institut de la Mer de Villefranche, France. Subsequently, the valve at the base of each tank, which allowed retrieval of the sedimentation bottles without disturbance, was opened, and the remaining water inside the tanks (165–180 L at TYR; 172.5 L at ION and 150 L at FAST) was filtered through a 100 µm mesh size PVC sieve. The organisms retained were gently removed using a washing bottle filled with filtered seawater (0.2 µm) and transferred directly into a 250 mL bottle and fixed with 4% final concentration formaldehyde. These samples were processed using a ZooSCAN (Hydroptic©; Gorsky et al., 2010) at the PIQv platform of EMBRC France. Organisms were identified and counted using automatic classification with a reference dataset in EcoTaxa (https://ecotaxa OBS-vlfr.fr/, last access: 17 April 2020), followed by manual validation.

## 2.3 Data analyses

### 2.3.1 Nutrient inputs from dust

The maximum percentage of dust-born dissolved N and P was estimated based on initial N and P composition of the dust analogue (see Sect. 2.1; Desboeufs et al., 2014) and maximal concentrations observed in tanks D and G at t1h and t6h after seeding, as follows:

$$ \%_{\text{dissolution}} = \frac{\text{CONC}_{\text{max}} - \text{CONC}_{\text{init}}}{\text{CONC}_{\text{dust}}} \cdot 100, \quad (1) $$

where CONC is the concentration of the corresponding nutrient in each tank before seeding (t0), CONCmax corresponds to the concentration of the corresponding nutrient in each tank when nutrient concentration was at a maximum within the first 6 h after seeding, and CONCdust is the maximum potential concentration, assuming a 100% dissolution from dust analogue (based on dust content; Desboeufs et al., 2014; Sect. 2.1).

### 2.3.2 Autotrophic and heterotrophic biomass

Given that samples for micro-phytoplankton counts were taken only at t-12h, as a first approximation, autotrophic biomass was estimated as the sum of Synechococcus, autotrophic pico-eukaryote and nano-eukaryote biomass (based on flow cytometry). Conversion of abundances to carbon units was carried out assuming 250 fg C cell⁻¹ for Synechococcus (Kana and Glibert, 1987). The biovolume to carbon content relationship of Verity et al. (1992) was used for autotrophic pico- and nano-eukaryotes assuming a spherical shape and a diameter of 2 and 6 µm, respectively. Heterotrophic biomass was computed as the sum of heterotrophic prokaryote (HP) biomass and heterotrophic nanoflagellates (HNF) biomass. Conversion to carbon biomass was done assuming 20 fg C per cell (Lee and Fuhrman, 1987) for heterotrophic prokaryotes and 220 fg C µm⁻³ (Børsheim and Bratbak, 1987) with a spherical shape and 3 µm diameter for heterotrophic nanoflagellates. The ratio of autotrophic to heterotrophic biomass during the experiments was used to evaluate the trophic status of the investigated communities and its evolution. Finally, a proxy for micro-phytoplankton biomass (Bmicro) was estimated following Vidussi et al. (2001), as the sum of fucoxanthin and peridinin.

# 3 Results

## 3.1 Initial conditions

Initial conditions at the three sampling stations while filling the tanks (t-12h before seeding) are shown in Table 2. Total alkalinity concentrations, pHT, increased from west to east (Table 2). NOx and DIP concentrations followed different patterns with highest NOx values at station FAST and highest DIP concentrations at station TYR. Consequently, the lowest NOx: DIP ratio was measured at TYR (0.8), compared to ION and FAST (2.8 and 4.6, respectively). Ammonium concentrations ranged from 0.045 µmol L⁻¹ to below detection limit at FAST. Silicate concentrations were similar at stations TYR and ION (~1 µmol L⁻¹) and higher than at station FAST (0.64 µmol L⁻¹).

Very low chlorophyll a concentrations were measured at the three stations (0.063–0.072 µg L⁻¹). The proportion of the different major pigments (Fig. 3) indicates that phytoplankton communities were similar with a dominance of prymnesiophytes (i.e. 19'-hexanoyloxyfucoxanthin; Ras et al., 2008) followed by cyanobacteria (i.e. Zeaxanthin; Ras et al., 2008) at stations TYR and ION. In contrast, at station FAST, the plankton community was clearly dominated by photosynthetic prokaryotes (i.e. zeaxanthin and divinylchlorophyll a as proxies for cyanobacteria and prochlorophytes, respectively; Ras et al., 2008). At all three stations, the proportion of pigments representative of larger species
<table><caption>**Table 2.** Initial conditions (sampling time t-12h) at stations TYR, ION and FAST measured while filling the tanks. pH<sub>T</sub>: pH on the total scale, NO<sub>x</sub>: nitrate + nitrite, NH<sub>4</sub>: ammonium, DIP: dissolved inorganic phosphorus, Si(OH)<sub>4</sub>: silicate, TChla: total chlorophyll a, HNF: heterotrophic nanoflagellates. The three most important pigments in terms of concentration are also presented (19'-hexanoyloxyfucoxanthin, zeaxanthin and divinyl chlorophyll a). Biomasses of the different groups analysed through flow cytometry were estimated based on conversion equations and/or factors found in the literature (see Sect. 2.3). Autotrophic and heterotrophic biomass based on flow cytometry (fraction < 20 µm). Values below detection limits are indicated as < dl.</caption><thead><tr><th></th><th>Sampling station</th><th>TYR</th><th>ION</th><th>FAST</th></tr></thead><tbody><tr><td></td><td>Coordinates (decimal)</td><td>39.34° N, 12.60° E</td><td>35.49° N, 19.78° E</td><td>37.95° N, 2.90° N</td></tr><tr><td></td><td>Bottom depth (m)</td><td>3395</td><td>3054</td><td>2775</td></tr><tr><td></td><td>Day, month, year and time of sampling (local time)</td><td>17/05/2017 17:00</td><td>25/05/2017 17:00</td><td>02/06/2017 21:00</td></tr><tr><td></td><td>Temperature (°C)</td><td>20.6</td><td>21.2</td><td>21.5</td></tr><tr><td></td><td>Salinity</td><td>37.96</td><td>39.02</td><td>37.07</td></tr><tr><td rowspan="2">Carbonate chemistry</td><td>pH<sub>T</sub></td><td>8.04</td><td>8.07</td><td>8.03</td></tr><tr><td>Total alkalinity (µmol kg<sup>-1</sup>)</td><td>2529</td><td>2627</td><td>2443</td></tr><tr><td rowspan="5">Nutrients</td><td>NO<sub>x</sub> (nmol L<sup>-1</sup>)</td><td>14.0</td><td>18.0</td><td>59.0</td></tr><tr><td>NH<sub>4</sub><sup>+</sup> (µmol L<sup>-1</sup>)</td><td>0.045</td><td>0.022</td><td>&lt; dl</td></tr><tr><td>DIP (nmol L<sup>-1</sup>)</td><td>17.1</td><td>6.5</td><td>12.9</td></tr><tr><td>Si(OH)<sub>4</sub> (µmol L<sup>-1</sup>)</td><td>1.0</td><td>0.96</td><td>0.64</td></tr><tr><td>NO<sub>x</sub> / DIP (molar ratio)</td><td>0.8</td><td>2.5</td><td>4.6</td></tr><tr><td rowspan="4">Pigments</td><td>TChla (µg L<sup>-1</sup>)</td><td>0.063</td><td>0.066</td><td>0.072</td></tr><tr><td>19'-hexanoyloxyfucoxanthin (µg L<sup>-1</sup>)</td><td>0.017</td><td>0.021</td><td>0.016</td></tr><tr><td>Zeaxanthin (µg L<sup>-1</sup>)</td><td>0.009</td><td>0.006</td><td>0.036</td></tr><tr><td>Divinyl chlorophyll a (µg L<sup>-1</sup>)</td><td>~0</td><td>0</td><td>0.014</td></tr><tr><td rowspan="8">Flow cytometry</td><td>Autotrophic pico-eukaryotes (cell mL<sup>-1</sup>; biomass in µg C L<sup>-1</sup>)</td><td>347.8; 0.5</td><td>239.9; 0.4</td><td>701.0; 1.0</td></tr><tr><td>Autotrophic nano-eukaryotes (cell mL<sup>-1</sup>; biomass in µg C L<sup>-1</sup>)</td><td>150.5; 3.9</td><td>188.8; 4.8</td><td>196.6; 5.0</td></tr><tr><td>Synechococcus (cell mL<sup>-1</sup>; biomass in µg C L<sup>-1</sup>)</td><td>4972; 1.2</td><td>3037; 0.8</td><td>6406; 1.6</td></tr><tr><td>Autotrophic biomass (µg C L<sup>-1</sup>)</td><td>5.6</td><td>6.0</td><td>7.7</td></tr><tr><td>Heterotrophic prokaryotes abundance (× 10<sup>5</sup> cell mL<sup>-1</sup>)</td><td>4.79</td><td>2.14</td><td>6.15</td></tr><tr><td>HNF (abundance in cell mL<sup>-1</sup>)</td><td>110.1</td><td>53.6</td><td>126.2</td></tr><tr><td>Heterotrophic biomass (µg C L<sup>-1</sup>)</td><td>9.9</td><td>4.5</td><td>12.7</td></tr><tr><td>Microscopy</td><td>Pennate diatoms (abundance in cell L<sup>-1</sup>)</td><td>140</td><td>520</td><td>880</td></tr><tr><td></td><td>Centric diatoms (abundance in cell L<sup>-1</sup>)</td><td>200</td><td>380</td><td>580</td></tr><tr><td></td><td>Dinoflagellates (abundance in cell L<sup>-1</sup>)</td><td>2770</td><td>3000</td><td>3410</td></tr><tr><td></td><td>Autotrophic flagellates (abundance in cell L<sup>-1</sup>)</td><td>0</td><td>60</td><td>650</td></tr><tr><td></td><td>Ciliates (abundance in cell L<sup>-1</sup>)</td><td>270</td><td>380</td><td>770</td></tr></tbody></table>

(i.e. fucoxanthin and peridinin; diatoms and dinoflagellates respectively; Ras et al., 2008) were very small (< 5%).

At all stations, autotrophic nanoplankton contributed most to total biomass. Autotrophic and heterotrophic biomass and abundances were highest at station FAST, followed by ION for the autotrophs and TYR for heterotrophs (Table 2). Differences in standing stocks between stations where more pronounced for the heterotrophs. As a consequence, the ratio between autotrophic biomass and heterotrophic biomass ranged from ~0.6 at TYR and FAST to 1.3 at ION.

## 3.2 Conditions of irradiance, temperature and pH during the experiments

Irradiance levels during the experiments are shown in Fig. 4. Decrease in water transparency after dust addition was observed at all three stations with the lowest impact at station FAST where irradiance levels decreased by only 60 µmol photons m<sup>-2</sup> s<sup>-1</sup> after dust addition, reach-

ing similar levels to those observed for tanks D and G. At station TYR, a more pronounced decrease was observed in acidified and warmed tanks (G1 and G2) with a decrease in daily average maximum irradiance of ~60 and ~160 µmol photons m<sup>-2</sup> s<sup>-1</sup> compared to dust-amended tanks D and controls, respectively. Temperature control (Fig. 4) was not optimal showing deviations between replicates of treatment G of up to 1.0 °C (station FAST). Temperature in controls and D tanks displayed a daily cycle, increasing during the day and decreasing at night (Fig. 4). The differences between the warmed treatment (G) and the other tanks were +3, +3.2 and +3.6 °C at TYR, ION and FAST, respectively.

Addition of CO<sub>2</sub>-saturated filtered seawater led to a decrease in pH<sub>T</sub> from 8.05 ± 0.004 (average ± SD of C1, C2, D1 and D2 at t0) to 7.74 (average between G1 and G2) at station TYR, from 8.07 ± 0.002 to 7.78 at station ION and from 8.05 ± 0.001 to 7.72 at station FAST (Fig. 5). pH<sub>T</sub> levels remained more or less constant in the control and D tanks
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Proportion of the different pigments, as measured by high-performance liquid chromatography (HPLC) in pumped surface seawater for the three experiments (t-12h).</figcaption></figure>

during all three experiments with no clear impact of dust addition. In G tanks, pH levels gradually increased during the experiments with larger variability between duplicates. These increases remained moderate thanks to the flushing of CO₂-enriched air above the tanks (pCO₂ of 1017 ± 11, 983 ± 96 and 1023 ± 25 ppm at TYR, ION and FAST, respectively; data not shown). Partial pressure of CO₂ in ambient air was 410 ppm, similar for the three stations. In all experiments, the addition of ¹³C bicarbonate led to an increase in total alkalinity between 6 and 11 µmol kg⁻¹ at t0. Dust addition, right after t0 in tanks D and G, led to a A_T decrease between 8 and 16 µmol kg⁻¹ at t24h with no apparent effects of warming and acidification. Overall, no large changes in A_T were observed during the experiments (Fig. 5).

## 3.3 Changes in nutrient concentrations

Dust addition led to a rapid increase in NOₓ (~11 µmol L⁻¹ as observed during the first 6 h; Fig. 6; Table 3) at all three stations with no differences between treatments D and G. The corresponding percent dissolution of N from dust ranged between 94% and 99%. In contrast, maximum DIP release was much smaller, ranging between 20 and 37 nmol L⁻¹, with slightly higher values at FAST (31–37 nmol L⁻¹) as compared to the other stations. Percent dissolution for DIP corresponded to 9.2% to 17.3% of total phosphorus contained in dust. As a consequence, NOₓ: DIP ratios increased from initial values below 5 to above 300, within 6 h after dust seeding, in tanks D and G (Fig. 6).

After the rapid increase in N and P, both nutrients decreased with time. While nutrient variability was small in control tanks (NOₓ and DIP variations below 20 and 3 nmol L⁻¹, respectively), large decrease in both elements

occurred in dust-amended tanks (D and G; Table 4). Similar linear decrease in NOₓ were observed throughout the experiments at stations TYR and ION with no visible differences between tanks D and G. In contrast, at station FAST, a more pronounced decrease in NOₓ was observed in dust-amended (D and G) tanks, as well as in warmed and acidified tanks relative to the D treatment. Nevertheless, at all stations, NOₓ concentrations in D and G treatments remained far above ambient levels throughout the experiments (> 9 µmol L⁻¹). Abrupt decreases in DIP were observed during the three experiments after the initial increase. At station TYR, after 24 h, all DIP released from dust decreased to initial levels in tanks G while it took two more days to reach initial levels in tanks D. In contrast, at station ION, no clear difference in DIP dynamics was observed between treatments D and G, with concentrations that decreased rapidly during the first 24 h but remained above initial levels until the end of the experiment. At station FAST, similarly to station TYR, DIP decreased rapidly from t12h in treatment G, reaching levels close to initial conditions at the end of the experiment. DIP decrease was much lower in treatment D (Table 4) with concentrations maintained far above ambient levels throughout the experiment. As a consequence of the differences between NOₓ and DIP dynamics as well as differences among stations, the NOₓ: DIP ratio increased, with clear differences between stations (Fig. 6), and remained much higher than in the controls.

At all stations, silicate concentrations were higher in dust-amended tanks relative to the controls. At TYR, while concentrations remained stable in control tanks, they increased linearly with time in the other tanks (D and G) with no apparent effect of the imposed increase in temperature and decrease in pH (i.e. tanks G). The difference in Si(OH)₄ concentrations between dust-amended treatments (D and G) and controls was ~0.1 µmol L⁻¹ at the end of the experiment. At station ION, after an initial decrease in concentrations between t-12h and t0, concentrations increased in all tanks until the end of the experiment with higher values in dust-amended tanks (D and G) than in controls and no difference between D and G treatments. In contrast, at FAST, concentrations increased from t-12h to t48h (with higher values in dust-amended tanks) and decreased onward until the end of the experiment. At the end of the experiment (t96h), Si(OH)₄ concentrations were higher in the G treatment than in the D treatment, which were similar to the controls.

## 3.4 Changes in biological stocks

Temporal dynamics in biological parameters showed very different patterns at each station. At TYR, total chlorophyll a concentrations did not change in the dust-amended D tanks (Fig. 7) and even led to slightly decreased values 24 h after dust addition (e.g. -35% to -38% in D1 and D2, respectively, as compared to controls; Table 5). No clear effects of dust addition (tanks D vs. C) were detectable for
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Continuous measurements of temperature and irradiance level (PAR) in the six tanks during the experiments at TYR, ION and FAST. The dashed vertical line indicates the time of dust seeding (after t0).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** pH on the total scale (pH<sub>T</sub>) and total alkalinity (A<sub>T</sub>) measured in the six tanks during the experiments at TYR, ION and FAST. The dashed vertical line indicates the time of dust seeding (after t0). Error bars correspond to the standard deviation based on analytical triplicates.</figcaption></figure>

all groups based on pigment analyses (Fig. 7). Results obtained based on flow cytometry counts (Fig. 8) were coherent with these observations and showed stronger decreases in cell abundances for < 20 µm autotrophic groups in tanks D1 and D2 (–77% to –80%). In contrast, the abundance of heterotrophic prokaryotes (HPs) increased rapidly after dust addition under both ambient (+53%–68%) and future (+68%) environmental conditions, with no clear difference among treatments. In warmed and acidified tanks (G), strong discrepancies between the duplicates were observed for pigments and autotrophic cell abundances: tank G1 showed moderate increases for all variables with the exception of autotrophic pico-eukaryotes, while in G2 all variables responded strongly to dust addition with maximum relative changes of > 300%, with the exception of autotrophic nano-eukaryotes. While HNF abundances responded posi-

tively to the treatments in D1, D2 and G2, abundances increased sharply in tank G1 towards the end of the experiment.

At ION, clear differences between treatments were observed for almost all pigments and cell abundances (Figs. 7, 8). With the exception of autotrophic nano-eukaryotes and HNF, all variables (pigments and cell abundances) increased as a response to both dust addition and warmed and acidified conditions (Table 5). The maximum relative changes as compared to controls observed for total chlorophyll *a* were 109%–183% and 399%–426% in tanks D and G, respectively. The highest stimulation by dust addition was observed for *Synechococcus* with +317%–390% and +805%–1425% increases in abundances in D and G tanks respectively (Table 5). Autotrophic nano-eukaryotes and HNF abundances did not respond to dust addition under ambient conditions, but an increase in abundances occurred
<table><caption>**Table 3.** Maximum input of nitrate + nitrite (NOₓ) and dissolved inorganic phosphorus (DIP) released from Saharan dust in tanks D and G as observed from the discrete samples taken during the first 6 h after seeding. The estimated maximal percentage of dissolution is also presented (see Sect. 2.3.1 for details on the calculations).</caption><thead><tr><th rowspan="2"></th><th colspan="4">NO<sub>x</sub></th><th colspan="4">DIP</th></tr><tr><th>D1</th><th>D2</th><th>G1</th><th>G2</th><th>D1</th><th>D2</th><th>G1</th><th>G2</th></tr><tr><th></th><th colspan="8">µmol L<sup>-1</sup></th></tr><tr><th></th><th colspan="8">nmol L<sup>-1</sup></th></tr></thead><tbody><tr><th colspan="9">Maximum input</th></tr><tr><th>TYR</th><td>11.0</td><td>11.1</td><td>11.1</td><td>11.0</td><td>24.6</td><td>20.4</td><td>24.6</td><td>23.9</td></tr><tr><th>ION</th><td>11.2</td><td>11.6</td><td>11.2</td><td>11.3</td><td>23.3</td><td>22.0</td><td>19.6</td><td>22.9</td></tr><tr><th>FAST</th><td>11.3</td><td>11.1</td><td>11.1</td><td>11.2</td><td>30.8</td><td>31.3</td><td>36.9</td><td>29.8</td></tr><tr><th colspan="9">Maximum dissolution (%)</th></tr><tr><th>TYR</th><td>95</td><td>96</td><td>95</td><td>94</td><td>12</td><td>10</td><td>12</td><td>11</td></tr><tr><th>ION</th><td>96</td><td>99</td><td>96</td><td>97</td><td>11</td><td>10</td><td>9</td><td>11</td></tr><tr><th>FAST</th><td>97</td><td>97</td><td>95</td><td>97</td><td>15</td><td>15</td><td>17</td><td>14</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Nutrients (nitrate + nitrite): NOₓ, dissolved inorganic phosphorus: DIP, silicate: Si(OH)₄ and the molar ratio between NOₓ and DIP, measured in each tank during the experiments at TYR, ION and FAST. The dashed vertical line indicates the time of seeding (after t0).</figcaption></figure>

in treatment G. In contrast to observations at TYR, temperature and pH affected heterotrophic prokaryotes in all dust-amended tanks at station ION with a higher impact of dust addition under future environmental conditions.

At station FAST, all biological stocks increased strongly after dust addition (Figs. 7, 8 and Table 5). Total chlorophyll *a* increased exponentially until the end of the experiment with slightly lower values observed under ambient environmental conditions (+237%–318% in D tanks
<table><caption>**Table 4.** Removal rate of nitrate + nitrite (NOₓ) and dissolved inorganic phosphorus (DIP) in tanks D and G during the three experiments (TYR, ION and FAST). For NOₓ, rates were estimated based on linear regressions between maximum concentrations (i.e. after dust enrichment, at t1h or t6h) and final concentrations (t72h for TYR and ION and t96h for FAST). For DIP, rates were estimated based on linear regressions between maximum concentrations (i.e. after dust enrichment at t1h or t6h) and concentrations after stabilization were observed. This sampling time is shown in parentheses. All rates are expressed in nmol L⁻¹ h⁻¹.</caption><thead><tr><th rowspan="2"></th><th colspan="3">NOₓ</th><th colspan="3">DIP</th></tr><tr><th>TYR</th><th>ION</th><th>FAST</th><th>TYR</th><th>ION</th><th>FAST</th></tr></thead><tbody><tr><th>D1</th><td>-6.5</td><td>-8.6</td><td>-14.3</td><td>-0.4 (t72h)</td><td>-0.5 (t48h)</td><td>-0.2 (t96h)</td></tr><tr><th>D2</th><td>-1.0</td><td>-8.6</td><td>-13.5</td><td>-0.3 (t72h)</td><td>-0.8 (t24h)</td><td>-0.2 (t96h)</td></tr><tr><th>G1</th><td>-6.7</td><td>-13.1</td><td>-21.6</td><td>-1.3 (t24h)</td><td>-0.8 (t24h)</td><td>-1.5 (t24h)</td></tr><tr><th>G2</th><td>-0.8</td><td>-1.6</td><td>-25.2</td><td>-1.3 (t24h)</td><td>-1.6 (t24h)</td><td>-1.1 (t24h)</td></tr></tbody></table>

<table><caption>**Table 5.** Percent (%) maximum relative changes in tanks D and G as compared to controls (average between C1 and C2), for the experiments TYR, ION and FAST. The sampling time at which these maximum relative changes were observed is shown in brackets. Tchla refers to the concentration of total chlorophyll *a* and *B*<sub>micro</sub> to the biomass proxy of micro-phytoplankton (sum of Fucoxanthin and Peridinin, see Material and methods) based on high-performance liquid chromatography (HPLC). HP and HNF refer to heterotrophic prokaryote and heterotrophic nanoflagellate abundances, respectively, measured by flow cytometry.</caption><thead><tr><th rowspan="2">Experiment</th><th rowspan="2">Tank</th><th colspan="2">HPLC</th><th colspan="5">Flow cytometry</th></tr><tr><th>TChla</th><th>B<sub>micro</sub></th><th>Autotrophic pico-eukaryotes</th><th>Autotrophic nano-eukaryotes</th><th>Synechococcus</th><th>HP</th><th>HNF</th></tr></thead><tbody><tr><td>TYR</td><td>D1</td><td>-35 (t24h)</td><td>-33 (t12h)</td><td>-75 (t72h)</td><td>-80 (t1h)</td><td>-71 (t48h)</td><td>68 (t72h)</td><td>352 (t72h)</td></tr><tr><td>TYR</td><td>D2</td><td>-38 (t12h)</td><td>-39 (t24h)</td><td>-75 (t72h)</td><td>-80 (t1h)</td><td>-72 (t48h)</td><td>53 (t72h)</td><td>100 (t72h)</td></tr><tr><td>TYR</td><td>G1</td><td>60 (t72h)</td><td>52 (t72h)</td><td>-75 (t1h)</td><td>89 (t72h)</td><td>76 (t72h)</td><td>67 (t72h)</td><td>1095 (t72h)</td></tr><tr><td>TYR</td><td>G2</td><td>359 (t72h)</td><td>392 (t72h)</td><td>323 (t72h)</td><td>119 (t72h)</td><td>700 (t72h)</td><td>68 (t48h)</td><td>298 (t72h)</td></tr><tr><td>ION</td><td>D1</td><td>183 (t72h)</td><td>157 (t72h)</td><td>126 (t72h)</td><td>89 (t72h)</td><td>317 (t72h)</td><td>128 (t72h)</td><td>44 (t72h)</td></tr><tr><td>ION</td><td>D2</td><td>109 (t72h)</td><td>156 (t72h)</td><td>117 (t72h)</td><td>-59 (t1h)</td><td>390 (t72h)</td><td>133 (t72h)</td><td>27 (t72h)</td></tr><tr><td>ION</td><td>G1</td><td>399 (t72h)</td><td>454 (t72h)</td><td>458 (t72h)</td><td>256 (t72h)</td><td>805 (t72h)</td><td>176 (t72h)</td><td>175 (t72h)</td></tr><tr><td>ION</td><td>G2</td><td>426 (t72h)</td><td>612 (t72h)</td><td>510 (t72h)</td><td>292 (t72h)</td><td>1425 (t72h)</td><td>161 (t72h)</td><td>129 (t72h)</td></tr><tr><td>FAST</td><td>D1</td><td>318 (t96h)</td><td>356 (t96h)</td><td>113 (t96h)</td><td>208 (t72h)</td><td>348 (t96h)</td><td>27 (t96h)</td><td>-38 (t96h)</td></tr><tr><td>FAST</td><td>D2</td><td>237 (t96h)</td><td>322 (t96h)</td><td>91 (t96h)</td><td>219 (t72h)</td><td>197 (t96h)</td><td>40 (t48h)</td><td>-49 (t96h)</td></tr><tr><td>FAST</td><td>G1</td><td>399 (t96h)</td><td>415 (t96h)</td><td>198 (t72h)</td><td>274 (t72h)</td><td>357 (t48h)</td><td>61 (t48h)</td><td>243 (t24h)</td></tr><tr><td>FAST</td><td>G2</td><td>395 (t96h)</td><td>421 (t96h)</td><td>129 (t72h)</td><td>202 (t96h)</td><td>344 (t48h)</td><td>67 (t48h)</td><td>74 (t24h)</td></tr></tbody></table>

and ~+400% in G tanks). Prymnesiophytes (i.e. 19'-hexanoyloxyfucoxanthin) and diatoms (i.e. Fucoxanthin) appeared as the groups benefiting the most from dust addition with no large impacts of warming and acidification, while pelagophytes (i.e. 19'-butanoyloxyfucoxanthin) and green algae (i.e. total chlorophyll *b*) showed a stronger response in treatment G. Finally, although cyanobacteria (i.e. zeaxanthin) responded faster to dust addition under future environmental conditions (tank G), this effect attenuated towards the end of the experiment. In contrast to estimates based on pigments, increases in cell abundances did not generally last until the end of the experiments. While abundances of autotrophic pico-eukaryotes increased until t96h in treatment D, abundances sharply declined between t72h and t96h for this group in treatment G. The same trend was observed for *Synechococcus*, although discrepancies between duplicates in treatment D at t96h did not allow conclusions to be drawn on the behaviour of this group by the end of the experi-

ment. Abundances of autotrophic nano-eukaryotes declined sharply between t72h and t96h under present and future conditions. The decline in HP abundances occurred earlier during the experiment with moderate maximum relative differences as compared to controls at t48h. HP abundances declined very sharply between t48h and t96h in treatment G, reaching control levels, while this decline was less sharp under present environmental conditions. Finally, HNF dynamics during this experiment was hard to interpret given the large increase in abundances in only one duplicate of treatment G (t24h) followed by a gradual decline.

Abundances of meso-zooplankton at the end of the experiments showed relatively similar values at stations TYR and ION, while much higher levels were observed at station FAST (Fig. 9). As a consequence of large variability between duplicates at stations TYR and ION, no clear effects of treatments were detected. At station FAST, although the sample size was too low to statistically test for differ-
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Total chlorophyll *a* and major pigments, from high-performance liquid chromatography (HPLC) measurements, in each tank during the experiments at TYR, ION and FAST. The dashed vertical line indicates the time of seeding (after t0).</figcaption></figure>

ences, higher total abundances of meso-zooplankton species were observed in the dust-amended tanks with no differences between ambient and future conditions of temperature and pH. However, differences in abundance were visible between these two treatments for specific groups, with higher abundance of Harosa and lower abundance of Crustacea (other than copepods) and Mollusca in warmed and acidified tanks, respectively.

# 4 Discussion

## 4.1 Initial conditions

During this study, the mixed layer depth (MLD) was somewhat shallower at TYR (20 m) than at ION and FAST (~10 and ~15 m, respectively) at the time of the sampling (Van Wambeke et al., 2020). Such shallow MLDs are characteristic of the stratified and oligotrophic conditions encountered
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Abundance of autotrophic pico-eukaryotes, autotrophic nano-eukaryotes, *Synechococcus*, heterotrophic prokaryotes (HP), and heterotrophic nano-flagellates (HNF), measured by flow cytometry, in each tank during the experiments at TYR, ION and FAST. The evolution of autotrophic biomass (see Material and methods for details on the calculation) is also shown. The dashed vertical line indicates the time of seeding (after t0).</figcaption></figure>

in the western Mediterranean Basin in late spring/early summer (D'Ortenzio et al., 2005). Although direct measurements of NOₓ and DIP concentrations using nanomolar techniques (as performed in our study) are scarce in the Mediterranean Sea, the low levels measured during the cruise are in agree-

ment with DIP values reported for the three basins (Djaoudi et al., 2018) and with NOₓ and DIP concentrations measured in coastal waters of Corsica in late spring/early summer (Louis et al., 2017b; Pulido-Villena et al., 2014; Ri-dame et al., 2014). NOₓ: DIP molar ratios in surface wa-
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Abundances of meso-zooplankton species as measured in each tank at the end of the experiments at TYR, ION and FAST.</figcaption></figure>

ters were well below the Redfield ratio (16:1) and are also consistent with previous studies. The low NOₓ: DIP ratios and nutrient concentrations suggest that communities found at the three stations experienced N and P co-limitation at the start of the experiments, as previously shown by Tanaka et al. (2011). Nutrient enrichment experiments confirmed that, at the three sites, heterotrophic bacteria were mainly N–P co-limited (Van Wambeke et al., 2021). In contrast to N and P, dissolved Fe in surface seawater ranged from 1.5 nmol L⁻¹ at TYR to 2.5 nmol L⁻¹ at ION (Roy-Barman et al., 2021) and was unlikely limiting for biological activity as previously shown in the Mediterranean Sea under stratified conditions (Bonnet et al., 2005; Ridame et al., 2014).

The low total chlorophyll *a* concentrations in surface waters were typical for the western and central Mediterranean Sea in late spring/early summer, as estimated from remote sensing (Bosc et al., 2004) and from in situ measurements (Manca et al., 2004). While large species (i.e. dinatoms, dinoflagellates) represented only ~10% of the total chlorophyll *a* biomass, the composition of the smaller size phytoplankton communities differed substantially, with autotrophic nano-eukaryotes dominating at stations TYR and ION and a larger contribution from autotrophic picoeukaryotes and Cyanobacteria at station FAST. Due to their low competitiveness under nutrient limitation, the small contribution of large phytoplankton cells at the start of the experiment is a fingerprint of LNLC areas in general, and of surface Mediterranean waters in late spring and summer (Siokou-Frangou et al., 2010).

Biomass of both heterotrophic nanoflagellates and prokaryotes followed a west-to-east gradient (FAST > TYR > ION), with high relative contribution by heterotrophs at stations TYR and FAST (60% of biomass) while at ION autotrophs contributed 60% to plankton biomass. Accordingly, net community production (NCP) rates (Gazeau et al., 2021) showed an initial community close to metabolic balance (mean ± SE: -0.06 ± 0.09 µmol

O₂ L⁻¹ d⁻¹) at ION and the highest community respiration rates and consequently lowest NCP rates at station TYR (-1.9 µmol O₂ L⁻¹ d⁻¹) suggesting that the autotrophic plankton community was not very active and relied on regenerated nutrients, as shown by the high level of NH₄⁺ at the start of the experiment at TYR. In contrast, although slightly heterotrophic (Gazeau et al., 2021) and limited by the low amount of nutrients, the community at FAST showed the highest levels of ¹⁴C production and heterotrophic prokaryote production (Gazeau et al., 2021) as well as N₂ fixation (Ridame et al., 2021). Altogether, the heterotrophic signature of the three investigated stations, although closer to metabolic balance at ION, reflected typical biogeochemical conditions in the Mediterranean Sea during late spring to early summer (Regaudie-de-Gioux et al., 2009).

## 4.2 Critical assessment of the experimental system and methodology

The experimental tanks used in this study have been successfully validated in previous studies designed to investigate the inputs of macro- and micro-nutrients (e.g. NOₓ, DIP, DFe) and the export of organic matter, under close-to-abiotic conditions (natural seawater filtered onto 0.2 µm) following simulated wet dust events using the same analogue as used in our study (Bressac and Guieu, 2013; Louis et al., 2017a, 2018). Louis et al. (2017a, 2018) further investigated these impacts under lowered pH conditions resulting in a rapid increase in pH levels in the acidified filtered seawater due to CO₂ outgassing (from ~7.4 to ~7.7 in six days). In the present study, our experimental system further allowed atmospheric pCO₂ in addition to light and temperature (i.e. climate reactors) to be controlled. Thereby, this allowed CO₂ outgassing to be significantly reduced and pH levels to be maintained close to their targets. The regulation of atmospheric CO₂ was, however, consistently more efficient in tank G2 compared to G1 (Fig. 5), resulting in a small discrepancy in terms of pH (highest difference of 0.04 pH units between the two G tanks at FAST), possibly due to a potential leak or a longer flushing time above tank G1. Nevertheless, as no systematic differences in nutrient dynamics and biological response were observed between the two tanks, and these small differences in pH had no detectable effect on the obtained results.

The lids above tanks, equipped with LEDs in order to reproduce sunlight intensity and spectrum, were used for the first time during these experiments. While simulated intensities were close to estimates for the northwestern Mediterranean Sea at 5 m depth in June (~1100 µmol photons m⁻² s⁻¹; Bernard Gentili, personal communication, 2017) and fairly consistent between duplicates under control and dust-amended conditions, the largest differences were also observed between tanks G1 and G2. These discrepancies could result from small differences in PAR sensor calibration and/or of different turbidity related to the amount of particles remaining in the tanks. As for pH,
replication in terms of macronutrient dynamics and biological response appeared satisfactory (except at station TYR; see below).

Continuous measurements in the tanks showed that temperature was not spatially homogeneous, leading to significant differences among replicates. This was more pronounced for warmed tanks (treatment G) with a maximum average difference over the experimental period of 0.7 °C during the FAST experiment. As for pH and light, these discrepancies did not systematically lead to observable differences in the investigated stocks and processes between duplicates (except at TYR, see below).

The necessity to carry out the incubations in a clean container limited our possibility to set up additional replicates for the three treatments. As described above, differences between duplicates were, for the vast majority of studied variables and processes, lower than differences between treatments and appear robust considering the difficulty to incubate plankton communities for which slight differences in initial composition can translate into important differences in dynamics (Eggers et al., 2014). Nevertheless, important discrepancies were detected for autotrophic stocks (in particular Synechococcus) as well as HNF and processes (Gazeau et al., 2021) for the warmed and acidified treatment (tanks G1 and G2) at station TYR. The reason behind these differences is most likely the grazing impact of heterotrophic nanoflagellates on prokaryotic pico-plankton (Sherr and Sherr, 1994) in tank G1 where HNF abundance sharply increased during the experiment. Overall, while the methodology used in this study allowed the impacts of dust addition to be successfully evaluated under both present and future environmental conditions at two out of three tested waters, the discrepancies at station TYR prevent us from drawing any strong conclusion on the effect of dust addition on the dynamics of the community under future environmental conditions at that station.

## 4.3 Impact of dust addition under present environmental conditions

During all experiments, the observed increases in NOₓ and DIP a few hours after dust addition under present environmental conditions were similar to the enrichment obtained during the DUNE experiments at the surface of the mesocosms (~50 m³) after the simulation of a wet dust deposition using the same dust analogue and the same simulated flux (Pulido-Villena et al., 2014; Ridame et al., 2014). The intensity of the simulated wet deposition event (i.e. 10 g m⁻²) represents a high but realistic scenario, as several studies reported even higher short wet deposition events in this area of the Mediterranean Sea (Bonnet and Guieu, 2006; Loçe-Pilot and Martin, 1996; Ternon et al., 2010). Furthermore, based on previous studies reporting the mixing between dust and polluted air masses during the atmospheric transport of dust particles (e.g. Falkovich et al., 2001; Putaud et al., 2004),

we used an evapo-condensed dust analogue that mimics the processes taking place in the atmosphere prior to deposition, essentially the adsorption of inorganic and organic soluble species (e.g. sulfate and nitrate; see Guieu et al., 2010a, for further details). The imposed evapo-condensation processes are responsible for the large nitrate-releasing capacity of the dust particles used in our study. As a consequence, the addition of new nutrients from dust in our study and during the P and R DUNE experiments were much higher, especially for NOₓ, than those observed by Pitta et al. (2017, and references therein) and Ridame et al. (2014) following the simulation of a dry Saharan dust deposition event. This confirms that wet dust deposition is a more efficient source of bioavailable nutrients than dry dust deposition.

Although NOₓ and DIP increases after dust addition were similar in all experiments, the subsequent dynamics of these elements and the impacts on plankton community composition and functioning were drastically different. While NOₓ levels decreased moderately over the course of our experiments due to biological uptake, more abrupt decreases were observed for DIP released by dust, reaching values close to the ones observed in the controls, except at station FAST where concentrations were still above ambient levels at the end of the experiment.

Previous experiments on the effect of dust addition in the Mediterranean Sea showed significant increases in chlorophyll *a* concentrations (mean ~90% increase; Guieu and Ridame, 2020). Interestingly, no stimulation of autotrophic biomass and primary production rates (Gazeau et al., 2021) was observed in dust-amended tanks under present conditions at station TYR. To the best of our knowledge, this is the first experimental evidence of a complete absence of response from an autotrophic community following dust wet deposition. The absence of response from autotrophic stocks could be due to a tight top-down control by grazers hiding potential responses from the autotrophic community (Lekunberri et al., 2010; Marañón et al., 2010) and/or a competition for nutrients with heterotrophic prokaryotes (Marañón et al., 2010). Feliú et al. (2020) have shown that the mesozooplankton assemblage at TYR was clearly impacted by a dust event that took place nine days before sampling at that station as evidenced from particulate inventory of lithogenic proxies (Al, Fe) in the water column (Bressac et al., 2021), likely stimulating phytoplankton growth and consequently increasing the abundance of herbivorous grazers (copepods) and attracting carnivorous species well before the start of the experiment. Heterotrophic bacteria are also limited by inorganic nutrients, mainly DIP, in oligotrophic systems (Obernosterer et al., 2003; Van Wambeke et al., 2001). Recent studies have shown significant increases in heterotrophic bacterial abundance, respiration and/or production following dust deposition (and nutrient enrichment) in these areas (Lekunberri et al., 2010; Pitta et al., 2017; Pulido-Villena et al., 2008; Romero et al., 2011). Heterotrophs appear to be more stimulated by dust pulses than autotrophic plankton with an
increasing degree of oligotrophy, modulated by the com-
petition for nutrients between phytoplankton and bacteria
(Marañón et al., 2010). This response was reflected at station
TYR, with heterotrophic prokaryotes reacting quickly and
strongly to nutrient addition both in terms of abundances and
production rates (Gazeau et al., 2021). These two aforemen-
tioned hypotheses are not mutually exclusive, and the quick
response of heterotrophic prokaryotes to dust addition is co-
herent with the net heterotrophy at this station (see Sect. 4.1)
due to increases in community respiration and decreases in
net community production rates in dust-amended as com-
pared to control tanks (Gazeau et al., 2021). Hence, dust ad-
dition to surface waters strongly dominated by heterotrophs
leads to a reduction of the capacity of these communities to
export organic matter and sequester atmospheric CO2.

In contrast to the dynamics of the experiment at TYR, stimulation of primary producers was observed at stations ION and FAST under present conditions with overall higher impact than previous studies compiled by Guieu and Ridame (2020). The largest increase in chlorophyll *a* concentrations at station FAST is coherent with NOₓ decreases observed at this station. Interestingly, at FAST, DIP concentrations were still above ambient conditions at the end of the experiment. Maximum primary production rates (¹⁴C incorporation) at the end of the experiment suggest strong DIP recycling and the dominance of regenerated production towards the end of the experiment (Gazeau et al., 2021). Although, in some cases, *Synechococcus* appeared stimulated by dust addition (Herut et al., 2005; Lagaria et al., 2017; Paytan et al., 2009), Guieu et al. (2014b) showed that, based on the analysis of several aerosols addition studies, this group had generally weak responses to aerosol addition in contrast to nano- and micro-phytoplankton, suggesting that aerosol deposition may lead to an increase in larger phytoplankton. Yet, at stations ION and FAST, the increase in *Synechococcus* abundance in dust-amended tanks was the highest relative to those of pico- and nano-eukaryotes. In particular, at station ION, no clear response to nutrient enrichment was observed for nano-eukaryotes throughout the experiment. However, it must be stressed that our experiments were of a relatively short duration (3 to 4 d). The sharp increase in Fucoxanthin paralleled by a decrease in silica, at the end of the experiment at station FAST where DIP limitation was not yet apparent, suggests a delayed response of diatoms as compared to smaller taxa. The sharp decline in nano-eukaryote abundances in dust-amended tanks at the end of the FAST experiment further suggests that this group reacted quickly to nutrient enrichment and was progressively grazed and/or outcompeted by larger phytoplankton species.

While all groups of primary producers benefited from nu-
trient enrichment at FAST, the increases in heterotrophic
prokaryote abundances were moderate, leading to an increase
in net community production rates throughout the experi-
ment, reaching positive levels and a autotroph : heterotroph
ratio of 4, while control tanks remained below metabolic bal-

ance (Gazeau et al., 2021). At station ION, the situation was
intermediate with a similar enhancement of both autotrophic
and heterotrophic stocks and no clear changes in the ra-
tio between autotrophic and heterotrophic biomass (data not
shown), although the system evolved towards net autotrophy
at the end of the experiment in dust-amended tanks under
present environmental conditions (Gazeau et al., 2021).

Transfer of newly produced organic matter to higher
trophic levels in the different treatments was assessed
through the quantification of meso-zooplankton abundance
at the end of each experiment. Altogether it is not surprising
that an increase in meso-zooplankton abundances was only
detected at station FAST where the strongest enhancement of
primary production was observed. Such an increase in meso-
zooplankton abundance in the dust-amended as compared to
control treatment was observed during land-based mesocosm
experiments in the eastern Mediterranean Sea (Pitta et al.,
2017).

Finally, although no clear effects of dust deposition under present conditions were detectable on autotrophic prokaryotes at station TYR, the strongest increase in N₂ fixation rates was recorded at this station (Ridame et al., 2021). However, the potential impact of this process on NOₓ concentration is negligible compared to the very large stock of NOₓ present in the dust-amended tanks, as less than 1 nmol L⁻¹ d⁻¹ of NOₓ was produced through N₂ fixation (Ridame et al., 2021).

4.4 Impact of dust addition under future
environmental conditions

Few studies have investigated the release and fate of nutri-
ents from atmospheric deposition under climate conditions
as expected for the end of the century, and, to the best of
our knowledge, our study represents the first attempt to test
for the combined effect of ocean warming and acidifica-
tion on these processes. The study by Louis et al. (2018),
carried out with filtered (0.2 µm mesh size) natural seawa-
ter using the same dust analogue and flux as in the present
study, showed that even an extreme ocean acidification sce-
nario (~ -0.6 pH units) does not impact the bioavailability
of macro- and micro-nutrients (NOₓ, DIP and DFe) in the
oligotrophic northwestern Mediterranean Sea. Similar results
were found by Mélançon et al. (2016) in high-nutrient low-
chlorophyll (HNLC) waters of the northeastern Pacific, un-
der a moderate ocean acidification scenario (-0.2 pH units).
As no differences were observed for NOₓ and DIP con-
centrations within a few hours following dust addition un-
der present and future environmental conditions, our results
agree with these previous findings and further highlight the
absence of a direct effect of ocean warming (+3°C) on the
release of nutrients from atmospheric particles.
In contrast, different nutrient consumption dynamics were observed between ambient and warmed and acidified tanks. No impacts of warming and acidification could be observed for NOₓ at stations TYR and ION due to low net uptake rates compared to the large increase following dust addition. In contrast, at the most productive station FAST, as a consequence of strongly enhanced biological stocks (see thereafter) and metabolic rates (Gazeau et al., 2021), larger NOₓ consumption rates were shown under future environmental conditions.

The differences in DIP dynamics between the two dust-amended treatments were more complex to interpret. A clear feature of our experiments is that, in contrast to present-day pH and temperature conditions, all the stock of DIP released from dust was consumed at the end of the three experiments under future conditions. The rate of decrease differed depending on the station. While DIP dynamics were quite similar between tanks maintained under present and future environmental conditions at ION, warming and acidification induced a faster decrease in DIP at TYR and FAST, with a full consumption of the released DIP within 24 h. An interesting outcome at station TYR was that, despite the important discrepancies observed for autotrophic stocks and metabolic rates between the duplicates G1 and G2 (see Sect. 4.2), a similar dynamics was observed for DIP concentrations in these tanks. As heterotrophic prokaryote biomass and production rates (Gazeau et al., 2021) did not differ between these duplicate tanks, this further highlights the clear dominance of heterotrophic processes at this station, a dominance which was exacerbated by dust addition under future environmental conditions, leading to an even stronger heterotrophic state at the end of this experiment (Gazeau et al., 2021).

At station ION, large impacts of warming and acidification were found with twice the chlorophyll a concentrations than in the dust-amended D tanks. At this station, all autotrophic groups increased with ocean acidification and warming. *Synechococcus* and to a lesser extent picoplankaryotes showed the strongest response. Yet these differences in abundance did not lead to detectable changes in the composition of the autotrophic assemblage, with nanoplankaryotes largely dominating carbon biomass at the end of this experiment (62% in treatment G vs. 64% in treatment D). Although the ratio between autotrophic and heterotrophic biomass appeared to be positively impacted under future environmental conditions, reaching values of up to 2 at the end of the experiment, warming and acidification led to a decrease in net community production (Gazeau et al., 2021) suggesting that in the coming decades the capacity of surface seawater to sequester anthropogenic CO₂ will be lowered.

Similarly, at FAST, all phytoplankton groups were impacted positively by warming and acidification with the strongest changes detected for *Synechococcus* as compared to present environmental conditions. However, in contrast to station ION, all groups reached maximal abundances (and carbon biomass) after 3 d of incubations, thereafter drasti-

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Maximum relative change (%) of main biological stocks (TCHla: total chlorophyll a, HP: heterotrophic prokaryotes) and processes (BP: bacterial production; PP: ¹⁴C-based primary production; see Gazeau et al., 2021; BR: bacterial respiration (no data from this study); and N₂ fixation, Ridame et al., 2021) obtained during the present study at the three stations (TYR, ION and FAST) under ambient conditions of pH and temperature (open red squares) and future conditions (full green squares). Vertical extension of each squares are delimited by the range of responses observed among the duplicates for each treatment. The dotted green squares for station TYR highlight the large variability observed between duplicates for some parameters and processes that prevented drawing solid conclusions. Box plots (Med) represent the distribution of responses observed from studies conducted in the Mediterranean Sea, as compiled by Guieu and Ridame (2020).</figcaption></figure>

cally decreasing most likely as a consequence of DIP limitation (see above). It must be stressed that this pattern could not be observed from pigments as no samples were taken for these analyses after 3 d of incubation. Also, in contrast to station ION, the abundance of heterotrophic prokaryotes in the warmer and acidified treatment reached a maximum after 2 d of incubations and then decreased rapidly to reach levels observed in the control treatment. This suggests that heterotrophic prokaryotes were the first to suffer from DIP limitation and further highlights the dominance of autotrophs in terms of nutrient consumption at this station. Although the ratio between autotrophic and heterotrophic biomass increased under future environmental conditions at ION, Gazeau et al. (2021) reported on a decrease in net community production rates in this treatment as compared to ambient environmental conditions, suggesting that, in the future,
nutrient release from dust will lead to a lesser sequestration capacity of surface waters for atmospheric CO$_2$.

The positive effects of warming and acidification on the abundance of mostly small (< 20 µm) phytoplankton taxa, as observed at ION and FAST, are in line with previously published studies. Although the effect of ocean acidification on small autotrophic species shows a wide range (e.g. Dutkiewicz et al., 2015), there is increasing evidence that small phytoplankton species will be favoured in a warmer ocean (e.g. Chen et al., 2014; Daufresne et al., 2009; Morán et al., 2010). Our experimental protocol was not conceived to distinguish temperature from pH effects; however, results concur with those of Maugendre et al. (2015), which further suggested temperature over elevated CO$_2$ as the main driver of increased pico-phytoplankton abundance in the Mediterranean Sea.

These enhanced fertilizing effects on primary producers at ION and FAST, under future as compared to present environmental conditions, did not seem to reach higher trophic levels as no clear differences in meso-zooplankton abundances were observed between ambient and warmed and acidified tanks at the end of the experiments. The duration of our experiments was too short to carefully assess the proportion of newly formed organic matter consumed by meso-zooplankton species and its effect on their biomass, yet group-specific variations were observed. Finally, Gazeau et al. (2021) did not observe an additional impact of future environmental conditions on the export of organic matter after dust addition.

# 5 Conclusion

These experiments conducted during the PEACETIME cruise represent the first attempt to investigate the impacts of atmospheric deposition on surface plankton communities both under present and future environmental conditions. Despite a few experimental issues, the three experiments provided new insights on these potential impacts in the open Mediterranean Sea. Stark differences in the response to dust deposition were observed between the three investigated stations in the Tyrrhenian Sea, Ionian Sea and in the Alaskan basin. Given that the initial conditions at the three stations were very similar in terms of nutrient and chlorophyll concentrations, these differences seem to be rather a consequence of the initial metabolic states of the community (autotrophy vs. heterotrophy). In all three cases, nutrient addition from dust deposition did not strongly modify but rather exacerbated this initial state. Relative changes in the main parameters presented in this paper and processes presented in Gazeau et al. (2021) as a consequence of dust addition under present and future environmental conditions are shown in Fig. 10 and compared to the compilation of published data for the Mediterranean Sea from Guieu and Ridame (2020). At station TYR, under conditions of a clear dominance of

heterotrophs in the use of resources and potentially a higher top-down control from grazers, dust addition drove the community into an even more heterotrophic state with no detectable effect on primary producers. At station ION, where the community was initially closer to metabolic balance, both heterotrophic and autotrophic compartments benefited from dust-derived nutrients. At FAST, the station with the highest initial autotrophic production, addition of nutrients led to an increase in both compartments, but heterotrophic prokaryotes became quickly P-limited, and overall larger effects were observed for phytoplankton. Ocean acidification and warming did not have any detectable impact on the release of nutrients from atmospheric particles. Furthermore, these external drivers did not drastically modify the composition of the autotrophic assemblage with all groups benefiting from warmer and acidified conditions. However, although for two out of the three stations investigated, larger increases were observed for autotrophic as compared to heterotrophic stocks under future environmental conditions, a stronger impact of warming and acidification on mineralization processes (Gazeau et al., 2021) suggests that, in the future, the plankton communities of Mediterranean surface waters will have a decreased capacity to sequester atmospheric CO$_2$ following the deposition of atmospheric particles.

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

*Author contributions.* FG and CG designed and supervised the study. FG, CG, CR and KD sampled seawater from the experimental tanks during the experiments. JMG and GDL participated in the technical preparation of the experimental system, and SA, CS, JOI, SM, SN, KD and JD performed sample analyses. FG, CR and CG wrote the paper with contributions from all authors.

*Competing interests.* The authors declare that they have no conflict of interest.

*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

*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.** The authors thank the captain and the crew of the R/V *Pourquoi Pas*? for their professionalism and their work at sea. We thank Julia Uitz, Céline Dimier and the SAPIGH HPLC analytical service at Institut de la Mer de Villefranche (IMEV) for sampling and analysis of phytoplankton pigments; John Dolan for microscopic counting; and Lynne Macarez and the PIQv-platform of EMBRC-France, a national Research Infrastructure supported by ANR, under the reference ANR-10-INSB-02, for meso-zooplankton analyses.

**Financial support.** This study is a contribution to the PEACETIME project, a joint initiative of the MERMEX and ChArMEx components supported by CNRS-INSU, IFREMER, CEA, and Météo-France as part of the programme MISTRALS coordinated by INSU. PEACETIME is a contribution to SOLAS and IMBER international programmes. The project was endorsed as a process study by GEOTRACES. The project leading to this publication has also received funding from the European FEDER Fund under project 1166-39417. Julie Dinasquet was funded by a Marie Curie Actions-International Outgoing Fellowship (PIOF-GA-2013-629378).

**Review statement.** This paper was edited by Christine Klaas and reviewed by two anonymous referees.

## References

Aminot, A. and Kérouel, R.: Dosage automatique des nutriments dans les eaux marines : méthodes en flux continu, Editions Ifremer, méthodes d'analyse en milieu marin, 188 pp., 2007.

Behrenfeld, M. J., O'Malley, R. T., Siegel, D. A., McClain, C. R., Sarmiento, J. L., Feldman, G. C., Milligan, A. J., Falkowski, P. G., Letelier, R. M., and Boss, E. S.: Climate-driven trends in contemporary ocean productivity, Nature, 444, 752–755, 2006.

Bergametti, G., Dutot, A.-L., Buat-Ménard, P., Losno, R., and Remoudaki, E.: Seasonal variability of the elemental composition of atmospheric aerosol particles over the Northwestern Mediterranean, Tellus B, 41, 353–361, https://doi.org/10.3402/tellusb.v4i13.15092, 1989.

Bonnet, S., Guieu, C., Chiaverini, J., Ras, J., and Stock, A.: Effect of atmospheric nutrients on the autotrophic communities in a low nutrient, low chlorophyll system, Limnol. Oceanogr., 50, 1810–1819, https://doi.org/10.4319/lo.2005.50.6.1810, 2005.

Bonnet, S. and Guieu, C.: Atmospheric forcing on the annual iron cycle in the western Mediterranean Sea: A 1-year survey, J. Geophys. Res-Oceans, 111, C09010, https://doi.org/10.1029/2005JC003213, 2006.

Børsheim, K. Y. and Bratbak, G.: Cell volume to cell carbon conversion factors for a bacterivorous Monas sp. enriched from seawater, Mar. Ecol.-Prog. Ser., 36, 171–175, 1987.

Bosc, E., Bricaud, A., and Antoine, D.: Seasonal and inter-annual variability in algal biomass and primary production in the Mediterranean Sea, as derived from 4 years of SeaWiFS observations, Global Biogeochem. Cy., 18, GB1005, https://doi.org/10.1029/2003GB002034, 2004.

Bressac, M. and Guieu, C.: Post-depositional processes: What really happens to new atmospheric iron in the

ocean's surface?, Global Biogeochem. Cy., 27, 859–870, https://doi.org/10.1002/gbc.20076, 2013.

Bressac, M., Guieu, C., Doxaran, D., Bourrin, F., Desboeufs, K., Leblond, N., and Ridame, C.: Quantification of the lithogenic carbon pump following a simulated dust-deposition event in large mesocosms, Biogeosciences, 11, 1007–1020, https://doi.org/10.5194/bg-11-1007-2014, 2014.

Bressac, M., Wagener, T., Leblond, N., Tovar-Sánchez, A., Ridame, C., Albani, S., Guasco, S., Dufour, A., Jacquet, S., Dulac, F., Desboeufs, K., and Guieu, C.: Subsurface iron accumulation and rapid aluminium removal in the Mediterranean following African dust deposition, Biogeosciences Discuss. [preprint], https://doi.org/10.5194/bg-2021-87, in review, 2021.

Chen, B., Liu, H., Huang, B., and Wang, J.: Temperature effects on the growth rate of marine picoplankton, Mar. Ecol.-Prog. Ser., 505, 37–47, https://doi.org/10.3354/meps10773, 2014.

Christaki, U., Courties, C., Massana, R., Catala, P., Lebaron, P., Gasol, J. M., and Zubkov, M. V.: Optimized routine flow cytometric enumeration of heterotrophic flagellates using SYBR Green I, Limnol. Oceanogr.-Meth., 9, 329–339, https://doi.org/10.4319/lom.2011.9.329, 2011.

Daufresne, M., Lengfellner, K., and Sommer, U.: Global warming benefits the small in aquatic ecosystems, P. Natl. Acad. Sci. USA, 106, 12788–12793, https://doi.org/10.1073/pnas.0902080106, 2009.

Desboeufs, K., Leblond, N., Wagener, T., Bon Nguyen, E., and Guieu, C.: Chemical fate and settling of mineral dust in surface seawater after atmospheric deposition observed from dust seeding experiments in large mesocosms, Biogeosciences, 11, 5581–5594, https://doi.org/10.5194/bg-11-5581-2014, 2014.

Desboeufs, K., Bon Nguyen, E., Chevaillier, S., Triquet, S., and Dulac, F.: Fluxes and sources of nutrient and trace metal atmospheric deposition in the northwestern Mediterranean, Atmos. Chem. Phys., 18, 14477–14492, https://doi.org/10.5194/acp-18-14477-2018, 2018.

Dickson, A. G., Sabine, C. L., and Christian, J. R.: Guide to best practices for ocean CO₂ measurements, PICES, Sydney, 191 pp., 2007.

Dinasquet, J., Bigeard, E., Gazeau, F., Azam, F., Guieu, C., Marañón, E., Ridame, C., Van Wambeke, F., Obernosterer, I., and Baudoux, A.-C.: Impact of dust addition on the microbial food web under present and future conditions of pH and temperature, Biogeosciences Discuss. [preprint], https://doi.org/10.5194/bg-2021-143, in review, 2021.

Djaoudi, K., Van Wambeke, F., Coppola, L., D'Ortenzio, F., Helias-Nunige, S., Raimbault, P., Taillandier, V., Testor, P., Wagener, T., and Pulido-Villena, E.: Sensitive Determination of the Dissolved Phosphate Pool for an Improved Resolution of Its Vertical Variability in the Surface Layer: New Views in the P-Depleted Mediterranean Sea, Front. Mar. Sci., 5, 234, https://doi.org/10.3389/fmars.2018.00234, 2018.

D'Ortenzio, F., Iudicone, D., Montegut, C. de B., Testor, P., Antoine, D., Marullo, S., Santoleri, R., and Madec, G.: Seasonal variability of the mixed layer depth in the Mediterranean Sea as derived from in situ profiles, Geophys. Res. Lett., 32, L12605, https://doi.org/10.1029/2005GL022463, 2005.

Duce, R. A., Liss, P. S., Merrill, J. T., Atlas, E. L., Buat-Ménard, P., Hicks, B. B., Miller, J. M., Prospero, J. M., Arimoto, R., Church, T. M., Ellis, W., Galloway, J. N., Hansen, L., Jickells, T. D.,
Knap, A. H., Reinhardt, K. H., Schneider, B., Soudine, A., Tokos, J. J., Tsunogai, S., Wollast, R., and Zhou, M.: The atmospheric input of trace species to the world ocean, Global Biogeochem. Cy., 5, 193–259, https://doi.org/10.1029/91GB01778, 1991.

Dutkiewicz, S., Morris, J. J., Follows, M. J., Scott, J., Levitan, O., Dyhrman, S. T., and Berman-Frank, I.: Impact of ocean acidification on the structure of future phytoplankton communities, Nat. Clim. Change, 5, 1002–1006, https://doi.org/10.1038/nclimate2722, 2015.

Eggers, S. L., Lewandowska, A. M., e Ramos, J. B., Blanco-Ameijeiras, S., Gallo, F., and Matthiessen, B.: Community composition has greater impact on the functioning of marine phytoplankton communities than ocean acidification, Glob. Change Biol., 20, 713–723, https://doi.org/10.1111/gcb.12421, 2014.

Emerson, S., Quay, P., Karl, D., Winn, C., Tupas, L., and Landry, M.: Experimental determination of the organic carbon flux from open-ocean surface waters, Nature, 389, 951–954, https://doi.org/10.1038/40111, 1997.

Falkovich, A. H., Ganor, E., Levin, Z., Formenti, P., and Rudich, Y.: Chemical and mineralogical analysis of individual mineral dust particles, J. Geophys. Res.-Atmos., 106, 18029–18036, https://doi.org/10.1029/2000JD900430, 2001.

Feliú, G., Pagano, M., Hidalgo, P., and Carlotti, F.: Structure and function of epipelagic mesozooplankton and their response to dust deposition events during the spring PEACETIME cruise in the Mediterranean Sea, Biogeosciences, 17, 5417–5441, https://doi.org/10.5194/bg-17-5417-2020, 2020.

Gazeau, F., Van Wambeke, F., Marañón, E., Pérez-Lorenzo, M., Alliouane, S., Stolpe, C., Blasco, T., Leblond, N., Zäncker, B., Engel, A., Marie, B., Dinasquet, J., and Guieu, C.: Impact of dust addition on the metabolism of Mediterranean plankton communities and carbon export under present and future conditions of pH and temperature, Biogeosciences Discuss. [preprint], https://doi.org/10.5194/bg-2021-20, in review, 2021.

Gorsky, G., Ohman, M. D., Picheral, M., Gasparini, S., Stemmann, L., Romagnan, J.-B., Cawood, A., Pesant, S., García-Comas, C., and Prejger, F.: Digital zooplankton image analysis using the ZooScan integrated system, J. Plankton Res., 32, 285–303, https://doi.org/10.1093/plankt/fbp124, 2010.

Guieu, C. and Ridame, C.: Impact of atmospheric deposition on marine chemistry and biogeochemistry, in: Atmospheric Chemistry in the Mediterranean Region: Comprehensive Diagnosis and Impacts, edited by: Dulac, F., Sauvage, S., and Hamonou, E., Springer, Cham, Switzerland, 2020.

Guieu, C., Dulac, F., Desboeufs, K., Wagener, T., Pulido-Villena, E., Grisoni, J.-M., Louis, F., Ridame, C., Blain, S., Brunet, C., Bon Nguyen, E., Tran, S., Labiadh, M., and Dominici, J.-M.: Large clean mesocosms and simulated dust deposition: a new methodology to investigate responses of marine oligotrophic ecosystems to atmospheric inputs, Biogeosciences, 7, 2765–2784, https://doi.org/10.5194/bg-7-2765-2010, 2010a.

Guieu, C., Loye-Pilot, M. D., Benyahya, L., and Dufour, A.: Spatial variability of atmospheric fluxes of metals (Al, Fe, Cd, Zn and Pb) and phosphorus over the whole Mediterranean from a one-year monitoring experiment: Biogeochemical implications, Mar. Chem., 120, 164–178, https://doi.org/10.1016/j.marchem.2009.02.004, 2010b.

Guieu, C., Ridame, C., Pulido-Villena, E., Bressac, M., Desboeufs, K., and Dulac, F.: Impact of dust deposition on car-

bon budget: a tentative assessment from a mesocosm approach, Biogeosciences, 11, 5621–5635, https://doi.org/10.5194/bg-11-5621-2014, 2014a.

Guieu, C., Aumont, O., Paytan, A., Bopp, L., Law, C. S., Mahowald, N., Achterberg, E. P., Marañón, E., Salihoglu, B., Crise, A., Wagener, T., Herut, B., Desboeufs, K., Kanakidou, M., Olgun, N., Peters, F., Pulido-Villena, E., Tovar-Sanchez, A., and Völker, C.: The significance of the episodic nature of atmospheric deposition to Low Nutrient Low Chlorophyll regions, Global Biogeochem. Cy., 28, 1179–1198, https://doi.org/10.1002/2014GB004852, 2014b.

Guieu, C., D'Ortenzio, F., Dulac, F., Taillandier, V., Doglioli, A., Petrenko, A., Barrillon, S., Mallet, M., Nabat, P., and Desboeufs, K.: Introduction: Process studies at the air-sea interface after atmospheric deposition in the Mediterranean Sea – objectives and strategy of the PEACETIME oceanographic campaign (May–June 2017), Biogeosciences, 17, 5563–5585, https://doi.org/10.5194/bg-17-5563-2020, 2020a.

Guieu, C., Desboeufs, K., Albani, S., Alliouane, S., Aumont, O., Barbieux, M., Barrillon, S., Baudoux, A.-C., Berline, L., Bhairy, N., Bigeard, E., Bloss, M., Bressac, M., Brito, J., Carlotti, F., de Liege, G., Dinasquet, J., Djaoudi, K., Doglioli, A., D'Ortenzio, F., Doussin, J.-F., Duforet, L., Dulac, F., Dutay, J.-C., Engel, A., Feliu-Brito, G., Ferre, H., Formenti, P., Fu, F., Garcia, D., Garel, D., Gazeau, F., Giorio, C., Gregori, G., Grisoni, J.-M., Guasco, S., Guittonneau, J., Haëntjens, N., Heimburger, L.-E., Helias, S., Jacquet, S., Laurent, B., Leblond, N., Lefevre, D., Mallet, M., Marañón, E., Nabat, P., Nicosia, A., Obernosterer, I., Perez, L., M., Petrenko, A., Pulido-Villena, E., Raimbault, P., Ridame, C., Riffault, V., Rougier, G., Rousselet, L., Roy-Barman, M., Saiz-Lopez, A., Schmechtig, C., Sellegri, K., Siour, G., Taillandier, V., Tamburini, C., Thyssen, M., Tovar-Sanchez, A., Triquet, S., Uitz, J., Van Wambeke, F., Wagener, T., and Zaencker, B.: Biogeochemical dataset collected during the PEACETIME cruise, SEANOE [data set], https://doi.org/10.17882/75747, 2020b.

Herut, B., Zohary, T., Krom, M. D., Mantoura, R. F. C., Pitta, P., Psarra, S., Rassoulzadegan, F., Tanaka, T., and Frede Thingstad, T.: Response of East Mediterranean surface water to Saharan dust: On-board microcosm experiment and field observations, Deep-Sea Res. Pt. II, 52, 3024–3040, https://doi.org/10.1016/j.dsr2.2005.09.003, 2005.

Holmes, R. M., Aminot, A., Kérouel, R., Hooker, B. A., and Peterson, B. J.: A simple and precise method for measuring ammonium in marine and freshwater ecosystems, Can. J. Fish. Aquat. Sci., 56, 1801–1808, https://doi.org/10.1139/f99-128, 1999.

IPCC: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 1535 pp., 2013.

IPCC: IPCC Special Report on the Ocean and Cryosphere in a Changing Climate, edited by: Pörtner, H.-O., Roberts, D. C., Masson-Delmotte, V., Zhai, P., Tignor, M., Poloczanska, E., Mintenbeck, K., Alegría, A., Nicolai, M., Okem, A., Petzold, J., Rama, B., and Weyer, N. M., in press, 2019.
Irwin, A. J. and Oliver, M. J.: Are ocean deserts getting larger?, Geophys. Res. Lett., 36, L18609, https://doi.org/10.1029/2009gl039883, 2009.

Jickells, T. D., An, Z. S., Andersen, K. K., Baker, A. R., Berga-
metti, G., Brooks, N., Cao, J. J., Boyd, P. W., Duce, R. A.,
Hunter, K. A., Kawahata, H., Kubilay, N., laRoche, J., Liss,
P. S., Mahowald, N., Prospero, J. M., Ridgwell, A. J., Tegen,
I., and Torres, R.: Global Iron Connections Between Desert
Dust, Ocean Biogeochemistry, and Climate, Science, 308, 67–
71, https://doi.org/10.1126/science.1105959, 2005.

Kana, T. M. and Glibert, P. M.: Effect of irradiances up to 2000 µE m⁻² s⁻¹ on marine *Synechococcus* WH7803—I. Growth, pigmentation, and cell composition, Deep-Sea Res., 34, 479–495, https://doi.org/10.1016/0198-0149(87)90001-X, 1987.

Kapsenberg, L., Alliouane, S., Gazeau, F., Mousseau, L., and Gat-
tuso, J.-P.: Coastal ocean acidification and increasing total al-
kalinity in the northwestern Mediterranean Sea, Ocean Sci., 13,
411–426, https://doi.org/10.5194/os-13-411-2017, 2017.

Kouvarakis, G., Mihalopoulos, N., Tselepides, A., and Stavrakakis, S.: On the importance of atmospheric inputs of inorganic nitrogen species on the productivity of the Eastern Mediterranean Sea, Global Biogeochem. Cy., 15, 805–817, https://doi.org/10.1029/2001GB001399, 2001.

Lagaria, A., Mandalakis, M., Mara, P., Papageorgiou, N., Pitta, P., Tsiola, A., Kagiorgi, M., and Psarra, S.: Phytoplankton response to Saharan dust depositions in the Eastern Mediterranean Sea: A mesocosm study, Front. Mar. Sci., 3, 287, https://doi.org/10.3389/fmars.2016.00287, 2017.

Law, C. S., Brévière, E., de Leeuw, G., Garçon, V., Guieu, C., Kieber, D. J., Kontradowitz, S., Paulmier, A., Quinn, P. K., Saltzman, E. S., Stefels, J., and von Glasow, R.: Evolving research directions in Surface Ocean – Lower Atmosphere (SOLAS) science, Environ. Chem., 10, 1–16, https://doi.org/10.1071/EN12159, 2013.

Lee, S. and Fuhrman, J. A.: Relationships between Biovolume and Biomass of Naturally Derived Marine Bacterioplankton, Appl. Environ. Microb., 53, 1298–1303, 1987.

Lekunberri, I., Lefort, T., Romero, E., Vázquez-Domínguez, E., Romera-Castillo, C., Marrasé, C., Peters, F., Weinbauer, M., and Gasol, J. M.: Effects of a dust deposition event on coastal marine microbial abundance and activity, bacterial community structure and ecosystem function, J. Plankton Res., 32, 381–396, https://doi.org/10.1093/plankt/fbp137, 2010.

Liu, X., Patsavas, M. C., and Byrne, R. H.: Purification and Char-
acterization of meta-Cresol Purple for Spectrophotometric Sea-
water pH Measurements, Environ. Sci. Technol., 45, 4862–4868,
https://doi.org/10.1021/es200665d, 2011.

Longhurst, A., Sathyendranath, S., Platt, T., and Caverhill, C.: An estimate of global primary production in the ocean from satellite radiometer data, J. Plankton Res., 17, 1245–1271, https://doi.org/10.1093/plankt/17.6.1245, 1995.

López-Urrutia, A. and Morán, X. A. G.: Resource limitation of bacterial production distorts the temperature dependence of oceanic carbon cycling, Ecology, 88, 817–822, https://doi.org/10.1890/06-1641, 2007.

Louis, J., Pedrotti, M. L., Gazeau, F., and Guieu, C.: Experimental evidence of formation of transparent exopolymer particles (TEP) and POC export provoked by dust addition under cur-

rent and high pCO₂ conditions, PLOS ONE, 12, e0171980,
https://doi.org/10.1371/journal.pone.0171980, 2017a.

Louis, J., Guieu, C., and Gazeau, F.: Nutrient dynam-
ics under different ocean acidification scenarios in a
low nutrient low chlorophyll system: The Northwestern
Mediterranean Sea, Estuar. Coast. Shelf S., 186, 30–44,
https://doi.org/10.1016/j.ecss.2016.01.015, 2017b.

Louis, J., Gazeau, F., and Guieu, C.: Atmospheric nutri-
ents in seawater under current and high pCO₂ condi-
tions after Saharan dust deposition: Results from three
minicosm experiments, Prog. Oceanogr., 163, 40–49,
https://doi.org/10.1016/j-pocean.2017.10.011, 2018.

Lože-Pilot, M. D. and Martin, J. M.: Saharan Dust Input to the Western Mediterranean: An Eleven Years Record in Corsica, in: The Impact of Desert Dust Across the Mediterranean, edited by: Guerzoni, S. and Chester, R., Springer Netherlands, Dordrecht, 191–199, https://doi.org/10.1007/978-94-017-3354-0-18, 1996.

Manca, B., Burca, M., Giorgetti, A., Coatanoan, C., Garcia, M.-J., and Iona, A.: Physical and biochemical averaged vertical profiles in the Mediterranean regions: an important tool to trace the climatology of water masses and to validate incoming data from operational oceanography, J. Marine Syst., 48, 83–116, https://doi.org/10.1016/j.jmarsys.2003.11.025, 2004.

Marañón, E., Fernández, A., Mouriño-Carballido, B., Martínez-García, S., Teira, E., Cermeño, P., Chouciño, P., Huete-Ortega, M., Fernández, E., Calvo-Díaz, A., Morán, X. A. G., Bode, A., Moreno-Ostos, E., Varela, M. M., Patey, M. D., and Achterberg, E. P.: Degree of oligotrophy controls the response of microbial plankton to Saharan dust, Limnol. Oceanogr., 55, 2339–2352, https://doi.org/10.4319/lo.2010.55.6.2339, 2010.

Marañón, E., Lorenzo, M. P., Cermeño, P., and Mouriño-Carballido, B.: Nutrient limitation suppresses the temperature dependence of phytoplankton metabolic rates, ISME J., 12, 1836–1845, https://doi.org/10.1038/s41396-018-0105-1, 2018.

Marie, D., Simon, N., Guillou, L., Partensky, F., and Vaulot, D.: Flow cytometry analysis of marine picoplankton, in: living color: protocols in flow cytometry and cell sorting, edited by: Diamond, R. A. and DeMaggio, S., Springer, Berlin, 421–454, 2010.

Markaki, Z., Oikonomou, K., Kocak, M., Kouvarakis, G., Chanio-
taki, A., Kubilay, N., and Mihalopoulos, N.: Atmospheric de-
position of inorganic phosphorus in the Levantine Basin, east-
ern Mediterranean: Spatial and temporal variability and its role
in seawater productivity, Limnol. Oceanogr., 48, 1557–1568,
https://doi.org/10.4319/lo.2003.48.4.1557, 2003.

Maugendre, L., Gattuso, J.-P., Louis, J., de Kluijver, A., Marro, S., Soetaert, K., and Gazeau, F.: Effect of ocean warming and acidification on a plankton community in the NW Mediterranean Sea, ICES J. Mar. Sci., 72, 1744–1755, https://doi.org/10.1093/icesjms/fsu161, 2015.

Maugendre, L., Guieu, C., Gattuso, J.-P., and Gazeau, F.: Ocean acidification in the Mediterranean Sea: Pelagic mesocosm experiments. A synthesis, Estuar. Coast. Shelf S., 186, 1–10, https://doi.org/10.1016/j.ecss.2017.01.006, 2017.

Mayot, N., D'Ortenzio, F., Ribera d'Alcalà, M., Lavigne, H., and Claustre, H.: Interannual variability of the Mediterranean trophic regimes from ocean color satellites, Biogeosciences, 13, 1901–1917, https://doi.org/10.5194/bg-13-1901-2016, 2016.

Mélançon, J., Levasseur, M., Lizotte, M., Scarratt, M., Tremblay, J.-É., Tortell, P., Yang, G.-P., Shi, G.-Y., Gao, H., Semeniuk,
D., Robert, M., Arychuk, M., Johnson, K., Sutherland, N., Davelaar, M., Nemcek, N., Peña, A., and Richardson, W.: Impact of ocean acidification on phytoplankton assemblage, growth, and DMS production following Fe-dust additions in the NE Pacific high-nutrient, low-chlorophyll waters, Biogeosciences, 13, 1677–1692, https://doi.org/10.5194/bg-13-1677-2016, 2016.

Moore, C. M., Mills, M. M., Arrigo, K. R., Berman-Frank, I., Bopp, L., Boyd, P. W., Galbraith, E. D., Geider, R. J., Guieu, C., Jaccard, S. L., Jickells, T. D., La Roche, J., Lenton, T. M., Mahowald, N. M., Marañón, E., Marinov, I., Moore, J. K., Nakatsuka, T., Oschlies, A., Saito, M. A., Thingstad, T. F., Tsuda, A., and Ulloa, O.: Processes and patterns of oceanic nutrient limitation, Nat. Geosci., 6, 701–710, https://doi.org/10.1038/ngeo1765, 2013.

Morán, X. A. G., López-Urrutia, Á., Calvo-Díaz, A., and Li, W. K. W.: Increasing importance of small phytoplankton in a warmer ocean, Glob. Change Biol., 16, 1137–1144, https://doi.org/10.1111/j.1365-2486.2009.01960.x, 2010.

Neale, P. J., Sobrino, C., Segovia, M., Mercado, J. M., Leon, P., Cortés, M. D., Tuite, P., Picazo, A., Salles, S., Cabrerizo, M. J., Prasil, O., Montecino, V., Reul, A., and Fuentes-Lema, A.: Effect of CO2, nutrients and light on coastal plankton. I. Abiotic conditions and biological responses, Aquat. Biol., 22, 25–41, https://doi.org/10.3354/ab00587, 2014.

Obernosterer, I., Kawasaki, N., and Benner, R.: P-limitation of respiration in the Sargasso Sea and uncoupling of bacteria from P-regeneration in size-fractionation experiments, Aquat. Microb. Ecol., 32, 229–237, https://doi.org/10.3354/ame032229, 2003.

Orr, J. C., Epitalon, J.-M., Dickson, A. G., and Gatutuso, J.-P.: Routine uncertainty propagation for the marine carbon dioxide system, Mar. Chem., 207, 84–107, https://doi.org/10.1016/j.marchem.2018.10.006, 2018.

Paytan, A., Mackey, K. R. M., Chen, Y., Lima, I. D., Doney, S. C., Mahowald, N., Labiosa, R., and Post, A. F.: Toxicity of atmospheric aerosols on marine phytoplankton, P. Natl. Acad. Sci. USA, 106, 4601–4605, https://doi.org/10.1073/pnas.0811486106, 2009.

Pitta, P., Kanakidou, M., Mihalopoulos, N., Christodoulaki, S., Dimitriou, P. D., Frangoulis, C., Giannakourou, A., Kagiorgi, M., Lagaria, A., Nikolaou, P., Papageorgiou, N., Psarra, S., Santi, I., Tsapakis, M., Tsiola, A., Violaki, K., and Petihakis, G.: Saharan Dust Deposition Effects on the Microbial Food Web in the Eastern Mediterranean: A Study Based on a Mesocosm Experiment, Front. Mar. Sci., 4, 117, https://doi.org/10.3389/fmars.2017.00117, 2017.

Polovina, J. J., Howell, E. A., and Abecassis, M.: Ocean's least productive waters are expanding, Geophys. Res. Lett., 35, L03618, https://doi.org/10.1029/2007gl031745, 2008.

Powley, H. R., Krom, M. D., and Cappellen, P. V.: Understanding the unique biogeochemistry of the Mediterranean Sea: Insights from a coupled phosphorus and nitrogen model, Global Biogeochem. Cy., 31, 1010–1031, https://doi.org/10.1002/2017GB005648, 2017.

Pulido-Villena, E., Wagener, T., and Guieu, C.: Bacterial response to dust pulses in the western Mediterranean: Implications for carbon cycling in the oligotrophic ocean, Global Biogeochem. Cy., 22, GB1020, https://doi.org/10.1029/2007gb003091, 2008.

Pulido-Villena, E., Rerolle, V., and Guieu, C.: Transient fertilizing effect of dust in P-deficient LNLC surface ocean, Geophys. Res. Lett., 37, L01603, https://doi.org/10.1029/2009GL041415, 2010.

Pulido-Villena, E., Baudoux, A.-C., Obernosterer, I., Landa, M., Caparros, J., Catala, P., Georges, C., Harmand, J., and Guieu, C.: Microbial food web dynamics in response to a Saharan dust event: results from a mesocosm study in the oligotrophic Mediterranean Sea, Biogeosciences, 11, 5607–5619, https://doi.org/10.5194/bg-11-5607-2014, 2014.

Putaud, J.-P., Van Dingenen, R., Dell'Acqua, A., Raes, F., Matta, E., Decesari, S., Facchini, M. C., and Fuzzi, S.: Size-segregated aerosol mass closure and chemical composition in Monte Cimone (I) during MINATROC, Atmos. Chem. Phys., 4, 889–902, https://doi.org/10.5194/acp-4-889-2004, 2004.

Ras, J., Claustre, H., and Uitz, J.: Spatial variability of phytoplankton pigment distributions in the Subtropical South Pacific Ocean: comparison between in situ and predicted data, Biogeosciences, 5, 353–369, https://doi.org/10.5194/bg-5-353-2008, 2008.

Regaudie-de-Gioux, A., Vaquer-Sunyer, R., and Duarte, C. M.: Patterns in planktonic metabolism in the Mediterranean Sea, Biogeosciences, 6, 3081–3089, https://doi.org/10.5194/bg-6-3081-2009, 2009.

Richon, C., Dutay, J.-C., Dulac, F., Wang, R., Balkanski, Y., Nabat, P., Aumont, O., Desboeufs, K., Laurent, B., Guieu, C., Raimbault, P., and Beuvier, J.: Modeling the impacts of atmospheric deposition of nitrogen and desert dust-derived phosphorus on nutrients and biological budgets of the Mediterranean Sea, Prog. Oceanogr., 163, 21–39, https://doi.org/10.1016/j.pocean.2017.04.009, 2018.

Ridame, C. and Guieu, C.: Saharan input of phosphate to the oligotrophic water of the open western Mediterranean Sea, Limnol. Oceanogr., 47, 856–869, 2002.

Ridame, C., Guieu, C., and L'Helguen, S.: Strong stimulation of N₂ fixation in oligotrophic Mediterranean Sea: results from dust addition in large in situ mesocosms, Biogeosciences, 10, 7333–7346, https://doi.org/10.5194/bg-10-7333-2013, 2013.

Ridame, C., Dekaezemacker, J., Guieu, C., Bonnet, S., L'Helguen, S., and Malien, F.: Contrasted Saharan dust events in LNLC environments: impact on nutrient dynamics and primary production, Biogeosciences, 11, 4783–4800, https://doi.org/10.5194/bg-11-4783-2014, 2014.

Ridame, C., Dinasquet, J., Hallström, S., Bigeard, E., Riemann, L., Van Wambeke, F., Bressac, M., Pulido-Villena, E., Taillandier, V., Gazeau, F., Tover-Sanchez, A., Baudoux, A.-C., and Guieu, C.: N2 fixation in the Mediterranean Sea related to the composition of the diazotrophic community, and impact of dust under present and future environmental conditions, Biogeosciences Discuss. [preprint], https://doi.org/10.5194/bg-2021-190, in review, 2021.

Romero, E., Peters, F., Marrasé, C., Guadayol, Ò., Gasol, J. M., and Weinbauer, M. G.: Coastal Mediterranean plankton stimulation dynamics through a dust storm event: An experimental simulation, Estuar. Coast. Shelf S., 93, 27–39, https://doi.org/10.1016/j.ecss.2011.03.019, 2011.

Roy-Barman, M., Foliot, L., Douville, E., Leblond, N., Gazeau, F., Bressac, M., Wagener, T., Ridame, C., Desboeufs, K., and Guieu, C.: Contrasted release of insoluble elements (Fe, Al, rare earth elements, Th, Pa) after dust deposition in seawater:
a tank experiment approach, Biogeosciences, 18, 2663–2678, https://doi.org/10.5194/bg-18-2663-2021, 2021.

Sala, M. M., Aparicio, F. L., Balagué, V., Boras, J. A., Borrull, E., Cardelús, C., Cros, L., Gomes, A., López-Sanz, A., Malits, A., Martínez, R. A., Mestre, M., Movilla, J., Sarmento, H., Vázquez-Domínguez, E., Vaqué, D., Pinhassi, J., Calbet, A., Calvo, E., Gasol, J. M., Pelejero, C., and Marrasé, C.: Contrast-ing effects of ocean acidification on the microbial food web un-der different trophic conditions, ICES J. Mar. Sci., 73, 670–679, https://doi.org/10.1093/icesjms/fsv130, 2016.

Sherr, E. B. and Sherr, B. F.: Bacterivory and herbivory: Key roles of phagotrophic protists in pelagic food webs, Microb. Ecol., 28, 223–235, https://doi.org/10.1007/BF00166812, 1994.

Siokou-Frangou, I., Christaki, U., Mazzocchi, M. G., Montresor, M., Ribera d'Alcalá, M., Vaqué, D., and Zingone, A.: Plankton in the open Mediterranean Sea: a review, Biogeosciences, 7, 1543–1586, https://doi.org/10.5194/bg-7-1543-2010, 2010.

Tanaka, T., Thingstad, T. F., Christaki, U., Colombet, J., Cornet-Barthaux, V., Courties, C., Grattepanche, J.-D., Lagaria, A., Nedoma, J., Oriol, L., Psarra, S., Pujo-Pay, M., and Van Wambeke, F.: Lack of P-limitation of phytoplankton and heterotrophic prokaryotes in surface waters of three anticyclonic eddies in the stratified Mediterranean Sea, Biogeosciences, 8, 525–538, https://doi.org/10.5194/bg-8-525-2011, 2011.

Ternon, E., Guieu, C., Loÿe-Pilot, M.-D., Leblond, N., Bosc, E., Gasser, B., Miquel, J.-C., and Martín, J.: The impact of Saharan dust on the particulate export in the water column of the North Western Mediterranean Sea, Biogeosciences, 7, 809–826, https://doi.org/10.5194 bg-7-809-2010, 2010.

The Mermex group: Marine ecosystems' responses to climatic and anthropogenic forcings in the Mediterranean, Prog. Oceanogr., 91, 97–166, https://doi.org/10.1016/j-pocean.2011.02.003, 2011.

Theodosi, C., Markaki, Z., Tselepides, A., and Mihalopoulos, N.: The significance of atmospheric inputs of soluble and particulate major and trace metals to the eastern Mediterranean seawater, Mar. Chem., 120, 154–163, https://doi.org/10.1016/j.marchem.2010.02.003, 2010.

Van Wambeke, F., Goutx, M., Striby, L., Sempéré, R., and Vidussi, F.: Bacterial dynamics during the transition from spring bloom to oligotrophy in the northwestern Mediterranean Sea: relationships with particulate detritus and dissolved organic matter, Mar. Ecol.-Prog. Ser., 212, 89–105, 2001.

Van Wambeke, F., Taillandier, V., Deboeufs, K., Pulido-Villena, E., Dinasquet, J., Engel, A., Marañón, E., Ridame, C., and Guieu, C.: Influence of atmospheric deposition on biogeochemical cycles in an oligotrophic ocean system, Biogeosciences Discuss. [preprint], https://doi.org/10.5194/bg-2020-411, in review, 2020.

Van Wambeke, F., Pulido, E., Catala, P., Dinasquet, J., Djaoudi, K., Engel, A., Garel, M., Guasco, S., Marie, B., Nunige, S., Taillandier, V., Zäncker, B., and Tamburini, C.: Spatial patterns of ectoenzymatic kinetics in relation to biogeochemical properties in the Mediterranean Sea and the concentration of the fluorogenic substrate used, Biogeosciences, 18, 2301–2323, https://doi.org/10.5194/bg-18-2301-2021, 2021.

Verity, P. G., Robertson, C. Y., Tronzo, C. R., Andrews, M. G., Nelson, J. R., and Sieracki, M. E.: Relationships between cell volume and the carbon and nitrogen content of marine photosynthetic nanoplankton, Limnol. Oceanogr., 37, 1434–1446, https://doi.org/10.4319/lo.1992.37.7.1434, 1992.

Vidussi, F., Claustre, H., Manca, B. B., Luchetta, A., and Marty, J.-C.: Phytoplankton pigment distribution in relation to upper thermocline circulation in the eastern Mediterranean Sea during winter, J. Geophys. Res.-Oceans, 106, 19939–19956, https://doi.org/10.1029/1999JC000308, 2001.