Effect of Intensity and Mode of Artificial Upwelling on Particle Flux and Carbon Export

Reduction of anthropogenic CO2 emissions alone will not sufficiently restrict global warming and enable the 1.5°C goal of the Paris agreement to be met. To effectively counteract climate change, measures to actively remove carbon dioxide from the atmosphere are required. Artificial upwelling has been proposed as one such carbon dioxide removal technique. By fueling primary productivity in the surface ocean with nutrient-rich deep water, it could potentially enhance downward fluxes of particulate organic carbon (POC) and carbon sequestration. In this study we investigated the effect of different intensities of artificial upwelling combined with two upwelling modes (recurring additions vs. one singular addition) on POC export, sinking matter stoichiometry and remineralization depth. We carried out a 39 day-long mesocosm experiment in the subtropical North Atlantic, where we fertilized oligotrophic surface waters with different amounts of deep water. The total nutrient inputs ranged from 1.6 to 11.0 μmol NO3– L–1. We found that on the one hand POC export under artificial upwelling more than doubled, and the molar C:N ratios of sinking organic matter increased from values around Redfield (6.6) to ∼8–13, which is beneficial for potential carbon dioxide removal. On the other hand, sinking matter was remineralized at faster rates and showed lower sinking velocities, which led to shallower remineralization depths. Particle properties were more favorable for deep carbon export in the recurring upwelling mode, while in the singular mode the C:N increase of sinking matter was more pronounced. In both upwelling modes roughly half of the produced organic carbon was retained in the water column until the end of the experiment. This suggests that the plankton communities were still in the process of adjustment, possibly due to the different response times of producers and consumers. There is thus a need for studies with longer experimental durations to quantify the responses of fully adjusted communities. Finally, our results revealed that artificial upwelling affects a variety of sinking particle properties, and that the intensity and mode with which it is applied control the strength of the effects.


INTRODUCTION
To limit global warming to between 1.5 and 2 • C as committed in the Paris Climate Agreement, reducing carbon dioxide (CO 2 ) emissions alone will most likely not suffice. Negative emission technologies, which actively remove CO 2 from the atmosphere, will be needed to achieve net zero CO 2 emissions (IPCC, 2018). Many such technologies focus on the ocean (GESAMP, 2019), which has the capacity to potentially store large amounts of extra carbon (Sabine et al., 2004). The world ocean has already absorbed roughly a third of cumulated anthropogenic CO 2 emissions (Khatiwala et al., 2013), thereby mitigating a substantial part of global warming. The ocean's storage capacity might be further enhanced by strengthening the carbon flux from the air-sea-interface to the deep ocean. One component of this flux is the biological carbon pump, which mediates the transport of organic carbon produced in the euphotic zone to the deep ocean. Once the carbon reaches the deep ocean, it remains out of touch with the atmosphere for decades to centuries (i.e., sequestered, Boyd et al., 2019). Some negative emission technologies focus on strengthening the biological carbon pump to enhance this natural carbon sink. In the oligotrophic ocean, increased vertical mixing and resulting nutrient-upwelling can enhance carbon fluxes to the deep ocean (Pedrosa-Pàmies et al., 2019). This principle is utilized by a negative emission technology called "artificial upwelling, " an approach in which nutrient-and CO 2 -rich deep ocean water is pumped to the surface layer with the goal of enhancing primary production and hence, carbon export.
The feasibility of artificial upwelling as a CO 2 -sequestration technique has been disputed by several modeling studies. Although it has the potential to sequester additional atmospheric carbon (Yool et al., 2009), Pan et al. (2015) found this potential to be rather small and connected to very high efforts; i.e., the upwelling of massive amounts of water (117 Sv) from 1,000 m depth would achieve an additional oceanic CO 2 uptake of 1 Gt C yr −1 . For comparison, in the year 2008 the oceans took up ∼2.5 Gt of anthropogenic carbon (Khatiwala et al., 2009). Furthermore, modeling studies have suggested that geophysical ocean-atmosphere feedbacks might lead to undesirable sideeffects, e.g., that stopping artificial upwelling after several decades of large-scale operation would increase global temperatures to levels higher than if it had not been applied at all (Oschlies et al., 2010;Keller et al., 2014;Kwiatkowski et al., 2015). Beside such geophysical feedbacks, other major unknowns are the biological and biogeochemical responses of pelagic communities to artificial upwelling. Basically all our knowledge about potential effects of artificial upwelling is based on models with a low degree of ecological complexity and simplified biogeochemistry (e.g., constant stoichiometry and particle properties). So far, there is little empirical research on the effects of artificial upwelling on the ecology and biogeochemistry of pelagic plankton communities.
The limited number of experimental studies so far have shown that artificial upwelling can significantly increase phytoplankton biomass in oligotrophic waters (McAndrew et al., 2007;Strohmeier et al., 2015;Giraud et al., 2016Giraud et al., , 2019Casareto et al., 2017). However, studies considering its effect on the export production and carbon sequestration processes are scarce, mostly due to the methodological difficulties of simulating artificial upwelling under close-to-natural conditions. Svensen et al. (2002) carried out a mesocosm nutrient enrichment experiment, in which they found that nutrient enrichment increases particulate organic carbon (POC) sedimentation and that the effect size is dependent on the frequency of nutrient additions. However, no study so far has examined in detail the properties of sinking particles, particularly sinking velocity and degradation rates, which together determine the remineralization depth and thus the potential for deep carbon export.
To create a net carbon sink using artificial upwelling, the export flux needs to overcompensate for the CO 2 upwelled with deep water . The most important export conditions determining carbon sequestration potential are (i) elemental stoichiometry of sinking particles and (ii) the remineralization depth of these particles. While artificially upwelled deep water usually contains more nitrogen and phosphorus than the surface water, it also contains more dissolved inorganic carbon (DIC), which could potentially outgas to the atmosphere. Hence, for artificial upwelling to increase the ocean's function as a net carbon sink, more carbon must be sequestered than is brought up as excess DIC (excess DIC = deep water DIC -surface water DIC). Carbon sequestration is thereby promoted by a high C:N ratio of sequestered particulate material (POC:PON). How much of the exported matter reaches the sequestration depth (i.e., the depth at which sinking material is regarded as sequestered) is determined by the remineralization length scale (RLS). It is a proxy for how deep particles can sink before being remineralized. Both parameters are dependent on the plankton community that develops as a response to artificial upwelling. Plankton community composition shapes not only the magnitude and C:N stoichiometry of the mass flux (Taucher et al., 2021), but also the characteristics of the sinking particles, such as their sinking velocities (SV) and carbon-specific remineralization rates (C remin ) (e.g., Legendre and Rivkin, 2002;Iversen and Ploug, 2010;Henson et al., 2012b;Turner, 2015;Bach et al., 2019). Community composition and succession are in turn controlled by the amount of added nutrients and the frequency with which they are supplied (Leibold et al., 1997;Jakobsen et al., 2015). This is why the intensity of artificial upwelling, as well as the mode of the application (i.e., the frequency of deep water addition) are important parameters that impact the export of organic matter.
Here, we present results from a mesocosm experiment where we studied the impact of different intensities and modes of artificial upwelling on arising food web dynamics and particulate matter export. We addressed the question whether increasing intensity of artificial upwelling leads to a higher mass flux and increased potential for carbon sequestration. Additionally, we compared the particle properties and their implications for deep carbon export of a singular upwelling event with a recurring mode of upwelling. The singular upwelling mode resembled an application of artificial upwelling that fertilizes each patch of water only once, e.g., by means of a moored wave pump (see e.g., Liu et al., 1999;Fan et al., 2016). The recurring mode on the other hand resembled an approach where the upwelling device drifts within a patch of water and fertilizes it over longer periods of time. Mesocosms are useful to study effects on whole pelagic communities, since they incorporate the responses of multiple trophic levels of a food web and allow for the assessment of a variety of ecological and biogeochemical parameters in an enclosed and well-characterized (eco-) system. Our study is the first to analyze the combined effects of upwelling intensity and upwelling mode on pelagic communities.

Experimental Setup
The mesocosm experiment was conducted from 5th of November 2018 and lasted for 39 days. Nine Kiel Off-Shore Mesocoms for Ocean Simulations (KOSMOS, see Riebesell et al., 2013 for technical information) were deployed in Gando Bay off the east coast of Gran Canaria (27 • 55.673 N,15 • 21.870 W). They enabled us to monitor the temporal development of the enclosed pelagic communities at in situ conditions. The mesocosms contained a 15 m long water column with a volume of ∼38 m 3 . We simulated artificial upwelling by replacing part of the oligotrophic mesocosm water with nutrient rich deep water. Our experimental design comprised one mesocosm as a control (i.e., no deep water addition), while the remaining eight received deep water additions in different modes and intensities. Four of these mesocosms were fertilized once with one big addition at the beginning of our study, i.e., the "singular treatments." The other four were fertilized eight times throughout the experiment with lower amounts of deep water per addition, but comparable total amounts of deep water summed up over the study period, i.e., the "recurring treatments." Both addition modes had four different levels of mixing ratio (low, medium, high, and extreme), with two respective mesocosms of different modes receiving similar amounts of deep water in total (Table 1).
We collected deep water for fertilization off the coast of Gran Canaria using a deep water collector, an opaque synthetic bag with 100 m 3 carrying capacity (see Taucher et al., 2017 for technical details), and moored it at the mesocosm site after each of two collections. The water was supposed to originate from ∼600 m depth, where NO 3 − concentrations are as high as 20-25 µmol L −1 (Llinás et al., 1994), in order to enhance primary and export production. However, due to loss of equipment under rough conditions we could not as planned collect water from that depth. Instead, we collected water from 330 m depth on the 26th of October 2018 (experimental day −10 = T-10) and from 280 m depth on the 28th of November (T23). Therefore, to reach the macronutrient concentrations necessary to achieve our planned experimental fertilization, we added nitrate (NO 3 − ), phosphate (PO 4 3− ) and silicate [Si(OH) 4 ] to the deep water prior to the first addition, resulting in concentrations of 25, 1.38, and 12.1 µmol L −1 , respectively. We calculated the amount of new N added to the mesocosms in the course of the whole experiment (Table 1) by calculating and adding up the net N inputs of each deep water replacement. For more details on the experimental setup and mesocosm activities (see Sswat et al., in prep.).

Sampling Procedure and Maintenance
Sampling for various parameters was generally carried out every second day (see Figure 1). A higher sampling frequency was pursued at the beginning of the experiment in order to cover the biological responses to the first two deep water additions in higher temporal resolution. The samples were either analyzed at our on-shore lab facilities at the Oceanic Platform of the Canary Islands (PLOCAN), the Marine Science and Technology Park (Parque Cientiìfico Tecnoloìgico Marino, PCTM) or the University of Las Palmas (ULPGC), or transported back to Kiel for analysis at the GEOMAR Helmholtz Centre for Ocean Research Kiel.
A CTD60M (Sea & Sun Technology GmbH, Trappenkamp, Germany) was cast to get depth profiles for temperature, salinity, density, pH, turbidity, oxygen (O 2 ) and photosynthetically active radiation. Bulk samples from the upper 13 m of the water column (WC) were taken for the analysis of primary productivity, chlorophyll a (Chl a), photosynthetic pigments, phyto-and microzooplankton analysis, prokaryotic heterotrophic production (PHP), and for the analysis of suspended particulate matter (PM wc ). The latter was composed of particulate organic and inorganic carbon (POC wc /PIC wc ), particulate organic nitrogen (PON wc ) and phosphorus (POP wc ), and biogenic silica (BSi wc ). Bulk samples were collected with integrated water samplers (HYDRO-BIOS Apparatebau GmbH, Kiel, Germany), which performed depth-integrated sampling from the upper 13 m of the water column. On each sampling day, around 40-60 L of bulk samples (i.e., 8-12 water samplers) were sampled per mesocosm. These were filled into 10 L carboys and stored dark and cool until arrival in the on-shore labs. Dissolved nutrients and carbonate chemistry (DIC and total alkalinity) were also sampled from the integrated water samplers.
Sedimented particulate matter was pumped out of the sediment trap (depth = 15 m) with a manual vacuum pump, not exceeding 0.3 bar during the process. Sediments were collected in 5 L glass bottles (Schott Scandinavia A/S, Kgs. Lyngby, Denmark) and stored in darkness until arrival in the lab, where they were subsampled for various parameters. The bottles were gently rotated to resuspend the material before homogeneous subsamples were taken. These subsamples were used for particle sinking velocity measurements and measurements of carbon-specific remineralization rates. For the latter, seven glass bottles (Schott, 310 mL volume) were additionally sampled with mesocosm water without headspace. They were used for incubating the sediment subsamples. The remaining sediment sample was weighed and analyzed for particulate carbon, nitrogen, phosphorus and biogenic silica content (see Boxhammer et al., 2016 for details of the method).
To prevent wall growth on the inside walls of the mesocosms, a ring-shaped wiper was pulled through each mesocosm roughly every 10 days (see Figure 1). This counteracted nutrient consumption by fouling organisms and their alteration of incoming light. For the latter reason, also the outside walls were cleaned twice during the study by divers equipped with brushes.
The sediment tubes of all mesocosms detached from the sediment traps on T30 due to strong winds and currents. They were reinforced and refitted on T31, however, no sediment Total new N added (µmol L −1 ) 0 1.6 3.1 5.6 9.8 1.6 3.1 6.2 11.0 FIGURE 1 | Timeline of the experiment depicting the sampling schedule, cleaning activities, and manipulations.
samples could be recovered that day and thus sedimented matter elemental composition, sinking velocity and remineralization rate could not be measured on T31. The tubes of the singular medium and high treatments and of the recurring extreme treatment disconnected again and had to be reattached a second time on T33. We were not able to sample the singular medium treatment that day, and the singular high and recurring extreme treatments had short incubation periods of only 18 h between the T33 and T35 samplings. We discarded all sinking velocity measurements from these days, since they showed systematically higher values than the measurements of the surrounding days. The mean seawater density inside all mesocosms between 0.5 and 15 m was slightly higher than in the surrounding water (1025.69 ± 0.11 and 1025.42 kg m −3 on T31, respectively), so that we assume no Atlantic water entered the small opening at the base of the sediment trap (ø ∼1 cm). An outflow of mesocosm water into the surrounding seems more likely under these circumstances. Nevertheless, we cannot rule out the possibility that e.g., due to wave and current action some Atlantic water entered the mesocosms. The surrounding water was, however, oligotrophic, and had substantially lower POC concentrations than any of the mesocosms on T33 according to our measurements (7.9 µmol POC L −1 in the Atlantic water compared to 35.7 ± 21.4 µmol POC L −1 in the mesocosms). This makes an influential contamination with biogenic material and/or species from the outside less probable, even if Atlantic water had entered the mesocosms. At the end of the experiment the water column in the mesocosms was mixed by pumping compressed air through the sediment hose. This might have altered the sediment flux and quality of the last sampling on T39.

Sediment Trap Material
At PLOCAN the sediment trap (ST) material was prepared for elemental analysis of POC ST , PON ST , and BSi ST by first of all separating the particles from the seawater. 3 mol L −1 ferric chloride (FeCl 3 ) were added to each 5 L bottle of sediment material to enhance flocculation and coagulation, followed by 3 mol L −1 NaOH addition to compensate for the decrease in pH (as described in detail in Boxhammer et al., 2016). After letting the material settle for 1 h, the supernatant was gently decanted. The remaining flocculated material was then centrifuged for 10 min at ∼5,200 g in a 6-16KS centrifuge (Sigma Laborzentrifugen GmbH, Osterode am Harz, Germany). An additional 10 min centrifugation step at ∼5,000 g in a 3K12 centrifuge (Sigma) resulted in compact sediment pellets, which were frozen at −20 • C and transported to Kiel for further processing. In Kiel, the pellets were freeze-dried to remove any leftover moisture and then ground in a cell mill (Edmund Bühler GmbH, Bodelshausen, Germany) to a fine homogeneous powder that was suitable for subsampling and further elemental analysis (Boxhammer et al., 2016). Subsamples for POC/N were weighed into tin capsules, acidified with 1 mol L −1 HCl, dried over night at 50 • C, and then measured in duplicate on a CN analyzer (Euro EA-CN, HEKAtech GmbH, Wegberg, Germany) according to Sharp (1974). To determine BSi concentrations, ∼2 mg subsamples of the sediment powder were measured spectrophotometrically following Hansen and Koroleff (1999). The sediment powder was stored dark and cool in glass vials. As some samples required multiple vials due to their high amounts, subsamples from all bottles were measured for each variable. If they differed among each other in elemental composition, the sample was pooled, homogenized, split up and remeasured until all bottles yielded the same results.

Water Column Samples
Water column samples were subsampled for elemental (POC WC /PON WC ) and pigment analysis (Chlorophyll a/Chl a) by collection on pre-combusted glass fiber filters (0.7 µm, Whatman) in our on-shore labs. POC and PON filters were acidified for ∼2 h to remove inorganic carbon, using 1 mol L −1 HCl, and dried over night at 60 • C in pre-combusted glass petri dishes. Filters for total particulate carbon (TPC) were dried without prior acidification. All filters were packed in tin cups (8 × 8 × 15 mm, LabNeed GmbH, Nidderau, Germany) and measured back in Kiel on a CN analyzer (Euro EA-CN, HEKAtech) as described for sediment C and N content above. Due to sample handling and/or measurement errors, measured POC wc concentrations were sometimes higher than those of TPC wc . Whenever the difference was greater than 10%, we report the TPC wc instead of the POC wc concentration (i.e., assuming that no PIC was present). Samples for Chl a were stored at −80 • C in cryovials until analysis in Kiel. They were extracted in acetone (100%) and homogenized with glass beads in a cell mill. After centrifugation (10 min, 5,200 rpm, 4 • C) they were filtered through 0.2 µm PTFE filters (VWR International GmbH, Darmstadt, Germany). Phytoplankton pigments, including Chl a in the supernatant were measured by an HPLC Ultimate 3,000 (Thermo Scientific GmbH, Schwerte, Germany).
Primary productivity was measured in the on-shore labs using a modified Nielsen (1952) 14 C uptake method described by Cermeño et al. (2012). Four water column subsamples (70 mL) per mesocosm were prefiltered and peaked with 2.96·10 5 Bq of 14 C-labeled sodium bicarbonate solution (NaH 14 CO 3 , PerkinElmer Inc., Waltham, United States). They were subsequently incubated in vitro for 24 h in a 12 h dark-light cycle. One of the subsamples was thereby covered with an opaque foil in order to measure the dark carbon uptake. Light intensity and temperature were set to in situ conditions based on CTD measurements. Thereafter, size fractions were filtered, the PIC and DIC fractions removed through acidification and the filters and filtrates were treated with a scintillation cocktail (Ultima Gold XR). Disintegrations per minute were counted on a scintillation counter (Beckmann LS-6500, Beckman Coulter Inc., Brea, United States). From these, knowing the concentration of added 14 C isotope and the measured in situ DIC concentration (see paragraph on DIC below), primary productivity rates (µmol C L −1 d −1 ) were calculated. For a more detailed description of the sampling and measurement procedure see Ortiz et al. (submitted).
PHP was estimated from rates of protein synthesis determined by the incorporation of tritiated leucine [( 3 H)leucine; Perkin Elmer] using the centrifugation method (Smith and Azam, 1992). Four subsamples (1 mL) and two trichloroacetic acid killed blanks were dispensed into screw-cap Eppendorf tubes. They were spiked with [ 3 H]leucine (final concentration: 20 nmol L −1 , specific activity 123 Ci mmol −1 ) and incubated at in situ temperature (21 • C) in the dark for 2-3 h. After the incubation, 100 µL of 50% trichloroacetic acid were added to the subsamples, which were kept with the blanks at −20 • C until centrifugation at 12 000 rpm for 20 min. The supernatant was carefully removed, and 1 mL of scintillation cocktail (Ultima Gold XR) was added to the Eppendorf tubes. They were stored in darkness for 24 h, after which the incorporated radioactivity was determined on a scintillation counter (Beckmann LS-6500). PHP was calculated using a conservative theoretical conversion factor of 1.55 kg C mol −1 Leu assuming no internal isotope dilution (Kirchman and Ducklow, 1993).
Samples for DIC measurements were filtered to remove particulate inorganic carbon (0.7 µm, Whatman) with an overflow of 1.5 times the volume of the filtrate. Inclusion of air in the filtration process was carefully avoided. The filtrate was fixed with mercuric chloride (HgCl 2 ). DIC concentrations were measured in triplicate on an AIRICA system (MARIANDA, Kiel, Germany) in Kiel, using a LI-COR LI-7000 Analyzer (LI-COR Biosciences GmbH, Bad Homburg, Germany). Certified reference materials (CRM batch 142, supplied by A. Dickson, Scripps Institution of Oceanography, United States) were used to determine the accuracy of DIC measurements.
Dissolved organic carbon (DOC) samples were filtered through pre-combusted GF/F filters (450 • C, 6 h) into polypropylene copolymer Nalgene TM bottles (Thermo Scientific), and stored at −20 • C. Prior to analysis, samples were acidified to a pH < 2 to remove inorganic carbon. Subsequently, DOC concentrations were measured by high temperature catalytic oxidation on a total organic carbon analyzer (TOC-V, Shimadzu Europa GmbH, Duisburg, Germany). The instrument was calibrated daily using potassium hydrogen phthalate (99.95-100.05%, p.a., Merck), which yielded an analytical precision of ± 1 µmol C L −1 . The accuracy was determined using certified reference materials (provided by D. A. Hansell, University of Miami, United States). Measured CRM concentrations were 43.39 ± 0.92 µmol C L −1 (n = 33), at a reference value of 42-45 µmol C L −1 .

Sinking Velocity of Sediment Trap Particles
Sinking velocity and remineralization rates of sedimented particulate matter were measured in a temperature-controlled on-shore lab at PCTM. Sinking velocity was determined by video microscopy using the method described in Bach et al. (2012). Sediment subsamples were diluted with filtered seawater according to their particle density (1:25-1:100) and injected to a sinking chamber (a cuvette with the dimensions: 10 × 10 × 350 mm), which was mounted vertically on a FlowCam 8000 (Fluid Imaging Technologies Inc., Scarborough, United States). Gravitational settling of particles in the cuvette was subsequently monitored for 20 min. Measurements were carried out at in situ temperatures (∼21 • C) under ventilation to prevent a temperature gradient around the sinking chamber. Particles between 25 and 1,000 µm were identified at a frame rate of 15-20 fps. The FlowCam measures more than 60 visually assessed parameters, which were read into MATLAB (version R2018b) for data analysis. The MATLAB script described by Bach et al. (2012 see: "Evaluation of sinking velocities") was adjusted for the more recent FlowCam version and used for calculation of sinking velocities. By finding multiple captures of the same particle on a y-position-gradient and applying a linear regression model against time, sinking velocities were calculated. They were corrected for wall effects of the sinking chamber according to the equation given by Ristow (1997). Further data analysis was carried out with the programming software R (R Core Team, 2017) using RStudio (version 1.3.959) and the package "tidyverse" (Wickham et al., 2019). An optical proxy for particle porosity (P int ) was calculated according to Bach et al. (2019). It is essentially a measure of the brightness of a particle, scaled with its size.
Analysis of the distribution of sinking particle volume showed that most of the volume was contained in a small fraction of particles with high equivalent spherical diameter (ESD) in all measurements (largest 10% accounted for > 75% of total particle biovolume). The smaller size fractions, however, constituted the majority of particle counts (∼90% in the 25-200 µm fraction). To not give the more numerous small particles undue weight, we calculated a "weighted SV." It gives more weight to the high-volume, but underrepresented bigger particles, which consequently contain more biomass (see Supplementary Figure 2 for an illustration of weighted SV calculation). In order to do so, we determined the ESD for which at least 25% of the summed-up particle volume was contained in the smaller particles and 75% in the bigger ones. For this "weighted ESD" we calculated the corresponding "weighted SV" by fitting a linear model to the mean sinking velocities as a function of their ESD across different size classes (log-spaced). We used size class means to-againnot give the numerous small particles undue weight. This weighted sinking velocity (hereinafter referred to as sinking velocity or SV) is thus the sinking velocity that corresponds to the ESD which segments the measurement into two particle volume fractions containing 25 and 75% of the total particle volume. We calculated the "weighted particle porosity" in the same way (hereinafter referred to as porosity or P int ). We used the weighted SV for all following sinking velocity dependent calculations such as the remineralization length scale (see below).

Remineralization Rates of Sinking Particles
Remineralization rates were determined every 4-6 days. We therefore took water column samples, collected in four replicate and three control bottles per mesocosm [(4 + 3) bottles × 9 mesocosms = 63 bottles in total]. They were transported back to the temperature-controlled lab (set to in situ temperature based on CTD measurements), where they acclimatized in a water bath for 2 h. Then, 0.5-3 mL of sediment subsample of the respective mesocosm were added to the four replicate bottles. The control bottles were left untreated. To examine the rate at which sedimented POC was remineralized back to DIC, all bottles were incubated in the dark on a rotating plankton wheel (∼1 rpm) and Oxygen depletion over time was measured. O 2 measurements were carried out using a handheld optical measurement device (Fibox4 Trace, PreSens-Precision Sensing GmbH, Regensburg, Germany), measuring non-invasively on PSt3 optodes (PreSens) mounted inside the bottles. O 2 measurements were automatically corrected for temperature (measured in a dummy bottle) and atmospheric pressure by the Fibox4. We adjusted the optode salinity correction according to the daily observed salinity measurements of the CTD cast. The second O 2 measurement was carried out after the bottles had been mixed properly, about 2 h after the incubation start, and were repeated continuously in 2-6 h intervals. The incubations lasted between 19 and 43 h, during which O 2 measurements were done at least 7 and up to 16 times. Particles in the incubation bottles were then collected on pre-combusted glass fiber filters (0.7 µm, Whatman) and analyzed for their POC content the same way as the water column POC filters.
By dividing the O 2 consumption rate (r in µmol O 2 L −1 d −1 ) of the sedimented matter by its POC content at the end of the incubation (µmol C L −1 ), the carbon-specific remineralization rate of the sedimented particulate matter (C remin in d −1 ) was calculated according to: where RQ is the respiratory quotient (µmol C µmol O 2 −1 ), which is commonly used as 1-mol CO 2 produced:1-mol O 2 consumed (= 1) (Ploug and Grossart, 2000;Iversen and Ploug, 2013;Bach et al., 2019) and t (d) as the time interval from the start of the incubation until the start of the filtration. In order to distinguish the response of the sedimented matter from any background seawater oxygen consumption, the C remin rates of the sedimentcontaining bottles were corrected for the mean C remin rates in the blank bottles.

Calculation of Export Flux and Remineralization Depth
The mass flux to the sediment trap (POC ST , PON ST , BSi ST ) was calculated for each element from the amount measured in the sediment powder. The total content per sample was calculated and normalized to the volume of the mesocosm and the time between sample collection (48 h), which yielded the daily mass flux in µmol per liter mesocosm water (µmol L −1 d −1 ). Ten data points were missing due to the detachment of the sediment trap hoses on and around T30 (all mesocosms on T31 and the singular medium mesocosm on T33, see section "Sampling Procedure and Maintenance"). To calculate the cumulative POC ST flux ( POC ST ), daily POC ST fluxes were summed up, whereby the missing data points were interpolated using the two surrounding data points (T29 and T33, or T29 and T35 for the singular medium treatment). Cumulative mass fluxes are reported in µmol L −1 .
To find out how much of the produced organic carbon had been exported from the water column until a specific experimental day (Tx), total organic carbon in the water column (TOC WC ) was calculated, subtracted by the initial TOC WC concentration on T01 (= TOC WC ), and finally compared to the cumulative POC ST flux ( POC ST ).
Finally, the remineralization length scale (RLS, i.e., remineralization depth) was calculated, which is the quotient of sinking velocity and carbon specific remineralization rate.

RLS = SV C remin
It is the depth (m) by which 63% of the sinking organic particle flux have been remineralized (Cavan et al., 2017) and thus a proxy for the POC transfer efficiency to depth.

Primary Production, Export Flux, and Stoichiometry
Primary productivity as well as Chl a and POC WC concentrations increased in all treatment mesocosms following the first deep water addition on T4. Phytoplankton blooms dominated by diatoms (Ortiz et al., submitted) developed in the singular treatments with intensity and duration depending on their upwelling intensity. The highest phytoplankton biomass (measured as Chl a concentration) was observed in the extreme singular treatment (11.2 µg L −1 on T9). Recurring upwelling sustained more stable phytoplankton biomass and productivity compared to singular upwelling. In the high and extreme recurring treatments primary productivity rates increased until well into the second half of the experiment (Figures 2A,B). Note that in most mesocosms a large fraction of produced organic matter was retained in the water column as suspended POC wc at the end of the experiment ( Figure 2C). This was particularly evident in the recurring treatments, which received nutrients until T32. In the extreme recurring treatment, the amount of POC WC at the end of the experiment was higher than its total cumulative POC mass flux (Figures 2C,E: 102 µmol L −1 POC WC vs. 82 µmol L −1 of cumulative POC ST , on T39).
The treatment differences observed in Chl a and primary productivity were mirrored in the POC mass flux. A postbloom export event occurred in the singular treatments 4-10 days after their bloom peak in the water column, with the time lag being longer at higher intensities of deep water addition. This is based on visual inspection (comparison of Figures 2B,D, e.g., extreme singular treatment bloom peak on T9, highest POC ST flux between T17-T21), since the temporal variation in the POC ST data did not allow a reliable calculation of time lags as e.g., in Stange et al. (2017). In contrast, the recurring mode led to POC ST maxima 20-30 days after initial fertilization ( Figure 2D). In both modes, the overall cumulated POC ST at the end of the experiment correlated positively with the cumulated primary productivity ( Figure 2F). The control mesocosm was excluded from this regression analysis, because it had unexplainably high export rates, although it did not receive any new N via deep water addition. When comparing the total amount of dissolved and particulate N at the beginning and end of the study period in the control, we found that approx. 4.2 µmol L −1 N sank out to the sediment trap, which were not accounted for in the initial N pools, according to our measurements. Apparently, there was an N source (or multiple) in the control, which we did not capture, e.g., fouling residues on the inside walls of the sediment trap or large swimmers (zooplankton, small fish), which were initially present but not accounted for. Other causes for the discrepancy might have been methodology-related, e.g., inaccurate sampling of patchily distributed PON WC . Regarding the test statistics, removing the control from the quantitative regression analysis in Figure 2F decreased the p-values (from 0.014 to 0.044 in the singular and from 0.015 to 0.060 in the recurring upwelling mode) and did not considerably affect the R-values (from 0.90 to The peaks in POC ST on T7 and T27 ( Figure 2D) in all mesocosms likely occurred due to the preceding cleaning of the inside mesocosm walls on T6 and T25, respectively. It appears that a carbon-rich biofilm was growing on the walls, which was removed during cleaning and subsequently sank into the sediment traps. This is also indicated by the high C:N ratios of POC flux on these days (compare Figures 2D, 3A). On T33 and T39, POC ST was lower than expected for most mesocosms, contrasting the time points before and after these days. For T33, the reason for this was probably a lower than usual accumulation period of material due to the preceding sediment tube detachments (see section "Sampling Procedure and Maintenance"). Many large and heavy particles might have sunk out before the reattachment of the tubes, leaving a higher proportion of slower sinking particles with lower POC content for the accumulation period. Regarding T39, the final fish net haul and the mixing of the water column through the sediment hose on T38 likely caused lower sedimentation rates of organic matter. The fish net haul by removing large and sticky aggregates, the mixing by resuspending already settled material.
In all treatment mesocosms except for the recurring low treatment more organic carbon was retained in the water column as POC and DOC than exported to the sediment trap as POC (ratio between POC ST and TOC WC < 1). The ratios did not substantially differ among the singular treatments, whereas they decreased with increasing upwelling intensity in the recurring mode ( Table 2).
In the singular upwelling treatments, the C:N ratios in the suspended particulate matter pool generally increased to higher than Redfield ratios on T11, 1 week after the deep water addition (Figure 3B). This occurred precisely when the blooming phytoplankton communities became nutrient limited (see Supplementary Figure 1 for inorganic nutrient concentrations), and was particularly prominent under high upwelling intensities. In the recurring treatments, the C:N ratios of POM WC showed a high temporal variation, they, however, tended to increase throughout the study period.
The C:N ratio of the mass flux was higher than the canonical Redfield ratio of 6.6 in all mesocosms throughout the experiment (Figure 3A). Increasing amounts of deep water addition, however, enhanced the mass flux C:N ratios further. We found a positive correlation between N addition and the C:N of the mass flux in the singular treatment between T11 and T21, when most of the post-bloom flux occurred ( Figure 3C). Furthermore, the Si:C ratio during this time was higher in the singular treatments compared to the recurring ones and the control (Figure 3D).

Particle Properties and Remineralization Depth
Remineralization rates of sinking particles ranged from 0.09 to 0.12 d −1 initially, and peaked in most mesocosms a week after the first deep water addition (T11). Thereafter, they decreased in all treatments ( Figure 4A). Sinking velocities were generally lower in the singular treatments than in the recurring ones ( Figure 4B). This difference between the two modes was confirmed by a t-test for mean sinking velocities calculated for days T7-T15 [two-sample t(6) = −4.019, p = 0.007]. Low sinking velocities in the singular mode were accompanied by high porosities in  particles > 90 µm ( Figure 4D). The difference in porosity between particles sinking under different upwelling modes is highlighted by the visual assessment of FlowCam pictures from T11 ( Figure 4E). Altogether, the singular treatments exported sedimented particulate matter which was carbon-and biogenic silica-rich and sank slowly during the export event. Particles > 100 µm were more porous compared to those in the recurring treatments. In contrast, the recurring treatments generally produced sinking particles which had lower C:N ratios, sank faster and were less porous than the ones in the singular treatments throughout the experiment.
We found that upwelling intensity had significant effects on sinking particle properties and on the remineralization depth during the time we observed increased mass flux (T11-T39, Figures 5A-C). Particles sampled from the sediment trap sank slowest and were remineralized most rapidly in mesocosms with high and extreme upwelling intensities. The fast remineralization rates correlated with high PHP in the water column during the same time period (Figure 5D). The RLS, as the quotient of SV and C remin , also decreased with increasing deep water fertilization.

Temporal Decoupling of Biomass Production and Particulate Organic Carbon Export
We found that increased intensities of artificial upwelling resulted in higher primary productivity and POC export from the surface. However, upwelling also resulted in a temporal decoupling between biomass production ( TOC WC ) and POC export (cumulative POC ST ). We found that less than half of the produced organic carbon was exported in most of our upwelling treatments during our experimental duration, while the rest was retained in the water column (see Table 2). This is consistent with findings by Taucher et al. (2017), who also fertilized Gran Canarian oligotrophic communities with a single pulse of nutrient-rich deep water, and found that roughly half of the produced particulate carbon was retained in the water FIGURE 5 | Correlations between mean carbon-specific remineralization rate (A), mean sinking velocity (B), mean remineralization length scale (C) and mean prokaryotic heterotrophic production (D) against cumulative N addition. It considers the time frame starting with the singular bloom export event until the end of the experiment (T11-T39), for which means of each of the parameters shown on y-axes were calculated. Since the recurring treatments received multiple deep water additions over time, mean cumulative N additions were calculated for the considered time period. Regression lines are plotted in the color of the corresponding upwelling mode (singular or recurring) and include the control treatment. Regression equations are provided in Supplementary Table 1. column for at least 25 days after the fertilization. Seemingly, artificial upwelling in oligotrophic surface waters does thus not lead to rapid export of most of the produced organic matter in the time scale of weeks.
Such decoupling between production and export occurs in a variety of oceanic regions (Cavan et al., 2015;Henson et al., 2019;Laws and Maiti, 2019), including eastern boundary upwelling systems (e.g., Kelly et al., 2018;Bach et al., 2020). It is thus no peculiarity of artificial upwelling. In natural upwelling systems reasons for a temporal decoupling can be e.g., persistent phytoplankton species in the water column (Bach et al., 2020) or weak trophic links in the planktonic food web (Stukel et al., 2011). A tight coupling of primary production and zooplankton consumers can generally lead to the export of high amounts of fast sinking fecal pellets, which can increase the POC export efficiency (Stukel et al., 2011;Le Moigne et al., 2016) and its transfer efficiency to depth (Steinberg and Landry, 2017). In our study, the response of the mesozooplankton community to the deep water fertilization was not very pronounced. Although a treatment effect on mesozooplankton biomass established toward the end of the experiment, values mostly stayed on a similar or lower level compared to before the first fertilization (∼20 µg C L −1 or lower, Spisla et al., in prep.). We also rarely observed fecal pellets in the FlowCam measurements of the sediment subsamples. The low mesozooplankton biomass and fecal pellet abundances indicate a limited impact of the mesozooplankton community on the POC flux. This hints at a potential caveat of short-term artificial upwelling as a tool for enhancing C sequestration: The fauna adapted to oligotrophic systems may not be able to readily consume great amounts of primary produced biomass and efficiently channel it into the export pathway. The dilution effect of the deep water addition amplifies the zooplankton handicap, making it even more difficult for them to catch up with the rapid phytoplankton reproduction.
Another factor that delayed the sedimentation of primary produced matter was surface microbial recycling.
Le Moigne et al. (2016) and Henson et al. (2019) found high bacterial abundances to be associated with high primary production and low export efficiency (i.e., only a small fraction of primary produced biomass is exported). Microbial recycling can channel particulate organic matter into the dissolved organic matter pool, thereby making it unavailable for metazoan consumers. Hence the breakdown of POC into DOC reduces the chances for either direct sinking of particles or incorporation and active transportation by zooplankton (Legendre and Le Fèvre, 1995). We found similar patterns of microbial activity as the above mentioned studies associated to the decoupling of primary production and export. Microbial production increased significantly with upwelling intensity in both upwelling modes (Figure 5D), thereby retaining organic matter in the surface microbial loop and making it unavailable for immediate export.
Due to the retention of organic matter in the water column, we were not able to capture the whole export response of our fertilized mesocosms, especially of the recurring highly fertilized ones. We argue that the high amounts of suspended POC (see Figure 2C) in the high and extreme recurring treatments at the end of the experiment represented a high potential for POC export, e.g., by particle aggregation and gravitational settling. Nevertheless, the retention of organic matter in the water column makes a quantitative mass flux comparison between the two upwelling modes difficult. Laws and Maiti (2019) found that the decoupling between primary and export production at the ALOHA time-series station disappears when taking into account data of monthly time intervals. We argue that the same is needed in future artificial upwelling studies. The experimental duration must be long enough to allow for the quantitative assessment of export fluxes as well as the longterm response of higher trophic levels to the enhanced surface primary production.

Upwelling Mode Shapes Sinking Particle Properties
The export events in the singular treatments were strongly linked to the preceding diatom blooms, as indicated by the elevated BSi:POC ratios in the sedimented matter ( Figure 3D and Supplementary Figure 1). They featured relatively porous particles, which sank slower than the ones in the recurring treatments (Figures 4B,D), and hence resulted in a low RLS ( Figure 4C). Multiple studies have shown that diatom blooms typically result in inefficient particle transfer to depth (e.g., Guidi et al., 2009;Henson et al., 2012a;Maiti et al., 2013), potentially due to the associated high particle porosities (Lam et al., 2011;Puigcorbé et al., 2015;Bach et al., 2019). Our data supports these observations, stressing the importance of a particle's porosity for its sinking velocity and thus for its potential to be transported to depth.
The influence of transparent exopolymer particles (TEP) may explain the low sinking velocities and high porosities in the singular treatments. TEP are exuded by microbial cells and facilitate the formation of aggregates (Passow, 2002;Engel et al., 2004). They can lower particle density and increase porosity, which in turn reduces particle sinking velocity (Azetsu-Scott and Passow, 2004;Mari et al., 2017). TEP production is favored by high primary productivity and nutrient limitation (Obernosterer and Herndl, 1995), both of which occurred in the higher singular treatments on T9 ( Figure 2B and Supplementary Figure 1). Additionally, the POC:PON ratio in the water column increased by > 60% from T9-T11 in the three highest singular treatments (mean values: 7.3 on T9-12.2 on T11, see Figure 4C), which indicates potentially enhanced TEP exudation. The same scenario has already been observed in a previous mesocosm experiment off Gran Canaria , in which TEP concentrations were measured. During that experiment, TEP increased when the blooming phytoplankton community became nutrient limited (compare Taucher et al., 2017(compare Taucher et al., , 2018. We thus argue that also in our experiment increased TEP exudation was likely the cause for the reduction of particle sinking velocity in the singular upwelling mode, rendering it less efficient regarding the potential for deep POC export than the recurring mode. Besides sinking velocity, respiration of sinking matter was another equally important factor determining the remineralization depth of sinking particles in our experiment. The temporal changes in C remin were equally high (about fourfold) as the changes in sinking velocity over time. We measured the highest remineralization rates during the initial export event on T11, shortly after phytoplankton communities had shifted to diatom dominance (T7-T9, Ortiz et al., submitted). This is in accordance with earlier studies, which speculated that the structure of the phytoplankton community affects remineralization rates (e.g., Guidi et al., 2015), and have shown that diatom-dominated systems can lead to the export of sinking particles which are quickly respired (Guidi et al., 2015;Bach et al., 2019). In our case, the mechanistic reason for this could be that the particles contained a large fraction of easily degradable organic matter. We believe that the increasing porosities from T11 onward ( Figure 4D) support this hypothesis. Although there are, to the best of our knowledge, no studies which found a direct positive effect of aggregate porosity on degradability, it stands to reason that higher porosities increase a particle's surfaceto-volume-ratio, thus possibly enhancing its susceptibility to microbial attachment and respiration. This is in line with a hypothesis by Francois et al. (2002), who postulate that the settling of loosely packed (i.e., highly porous) aggregates leads to higher respiration rates in the mesopelagic zone. Respiration of sinking particles during this time was higher in the singular compared to the recurring treatments, likely owing to higher porosities and possibly also higher TEP concentrations. Beside the phytoplankton community structure, the zooplankton community is another factor that determines particle properties (Cavan et al., 2019). Through the packaging of phytoplankton cells into dense fecal pellets, zooplankton can make organic matter less susceptible to respiration (Steinberg and Landry, 2017). As discussed above, mesozooplankton was not abundant in our experiment, and therefore potentially not capable of repackaging the high amounts of particulate matter in the high upwelling treatments. This left the sinking biomass loosely packed for the bacterial community to feed on. We suggest that the missing repackaging of particles and the resulting high lability of sinking organic matter was the reason for increasing remineralization with increasing upwelling intensity.
The differences in particle properties between our singular and recurring form of artificial upwelling indicate variability in the efficiency of their POC transfer to depth. Diatom spring blooms (Martin et al., 2011) and diatom communities in coastal upwelling systems (Abrantes et al., 2016) have been reported to promote efficient carbon deep export. However, the diatom blooms in our singular upwelling mode led to the export of porous, slow sinking particles, which were quickly remineralized by microbial processes. A single pulse of artificial upwelling thus resulted in a system likely to recycle the vast majority of the freshly produced organic matter in the surface ocean. Although particles in the recurring upwelling mode were respired similarly quickly, they were less porous and sank faster, which led to a slightly higher mean RLS in the recurring upwelling treatment ( Figure 5C). This suggests a more efficient transfer of produced organic matter to depth under recurring upwelling conditions compared to a singular upwelling pulse.

Potential of Artificial Upwelling for Carbon Sequestration
As discussed above, one of the major factors that influence the potential of artificial upwelling for net carbon sequestration is how deep particles sink before they are remineralized (RLS). Another factor is the C:N ratio of the sinking organic matter (POC:PON), or, more generally speaking, the ratio of C to the limiting nutrient. The influence of both of these factors on carbon sequestration efficiency will be discussed in the following. We emphasize that the discussion is based upon a highly simplified 1-dimensional view of the water column, and neglects 3 dimensional movements of water masses, carbon, and nutrients through the ocean.
We start this discussion with the influence of the ratio of carbon to the limiting nutrient. Since the latter was nitrogen in our case, as is typical for large parts of the Atlantic (Moore et al., 2013), we will focus on the C:N ratio in the following. Let us assume that over long enough time scales all upwelled nutrients will be sequestered (i.e., exported to the sequestration depth), and that the upwelled deep water contains excess DIC (excess DIC = deep water DIC -surface water DIC) and excess DIN (excess DIN = deep water DIN -surface water DIN) compared to the surface ocean. In this case, a system will act as a net carbon sink if it sequesters organic matter with a C:N ratio higher than that of the upwelled excess DIC and nitrogen. In the Canary Island region, excess DIC and excess DIN between surface and 1,000 m depth (here assumed as sequestration depth) are approximately ∼150 and ∼20 µmol L −1 , respectively (Llinás et al., 1994;González-Dávila et al., 2010). This corresponds to a C:N ratio of 7.5 of upwelled excess DIC and DIN. Accordingly, for the artificial upwelling approach to generate a net carbon sink, the C:N of sequestered organic matter would need to exceed 7.5 on average, until all N added via deep water has been sequestered.
We found that artificial upwelling led to an increase in C:N ratios of sedimented matter to well above 7.5. There was a trend of increasing POC ST :PON ST ratios over time in all mesocosms ( Figure 3A). The ratios in the singular treatment mesocosms correlated positively with upwelling intensity in the period of highest export (T11-T21, Figure 3C). It is likely that these high ratios would have either prevailed with depth or increased even further, owing to preferential N remineralization by heterotrophic microorganisms (Boyd and Trull, 2007).
It should be noted that the control mesocosm also displayed increasing C:N ratios during the experiment. This suggests that increasing C:N might have at least partly been an enclosure effect, which also occurred in absence of an artificial upwelling treatment. By enclosing a planktonic community inside a mesocosm, the vertical mixing is reduced from the depth of the mixed layer (20-120 m in the Canary Island region, see Troupin et al., 2010) to the depth of the mesocosm (15 m in our setup). This leads e.g., to an increase of depth-and time-integrated light intensity experienced by the phytoplankton community, which may affect phytoplankton photophysiology. Moreover, the community gets truncated as large grazers (>3 mm) are excluded from the food web, thus relieving lower trophic levels like microzooplankton of grazing pressure. Both effects could have caused changes in lower trophic level productivity in our mesocosm experiment, which was fueled by leftover nutrients enclosed at experiment start (see Supplementary Figure 1), thereby leading to an increase in POC WC and POC ST :PON ST ratios (Figures 2C, 3A). Nevertheless, the increase in sedimented matter C:N ratios was larger in those mesocosms that received deep water fertilizations, and scaled almost linearly with high upwelling rates.
Artificial upwelling can consequently enhance C:N ratios of sinking organic matter in the surface, which would favor potential carbon sequestration. This is especially the case when fertilization is carried out in single pulses. We cannot assess, however, if these ratios would have prevailed over longer time scales, i.e., until all added N had been exported with the remaining suspended particulate organic matter.
A mechanism that might counteract the positive effect of elevated C:N ratios is the observed decrease in RLS with increasing upwelling intensity in both upwelling modes. Particle sinking velocities decreased while at the same time respiration of sinking organic material accelerated with increasing nutrient input via deep water addition. Both effects resulted in a lower RLS and thus higher flux attenuation with increasing fertilization (Figures 5A-C). Earlier studies have shown that a higher fraction of export production makes it through the mesopelagic zone in less productive ecosystems and ends up in the deep ocean (e.g., Henson et al., 2012b). Our results match this finding, with less productive mesocosms showing a higher RLS than the very productive ones. These differences in RLS and thus particle transfer efficiency to depth, as well as the differences in C:N stoichiometries of sinking matter are illustrated in Figure 6.
It should be noted that our study could not assess several factors that might play an important role in determining the potential of artificial upwelling for carbon export and sequestration. First of all, our study did not account for any active biological or physical processes occurring further down the water column. Processes such as repackaging and FIGURE 6 | Scheme depicting the effects that upwelling intensity had on carbon sequestration-relevant parameters in our high and low upwelling scenarios. It visualizes the strength and C:N ratio of the export flux (POC ST :PON ST ), as well as the arrow size indicating strength of the recycling pathway and the efficiency of particle transport to depth. Note that we put more emphasis on the difference between POC:PON ratios than on the quantitative PON difference between high and low upwelling intensities. It should also be noted that the effects depicted here only reflect the initial responses to artificial upwelling in oligotrophic waters. They do not cover the ultimate fate of the upwelled inorganic nutrients and hence do not represent the long-term potential for carbon sequestration. consumption of settling particles by zooplankton (Turner, 2015;Stukel et al., 2019) or the active diel migration of mesopelagic biota  and references therein) play an important role in regulating particle transfer to depth and thus carbon sequestration (Buesseler et al., 2007;Robinson et al., 2010;Giering et al., 2014;Sanders et al., 2016;Cavan et al., 2019). Due to the limited vertical dimensions of our mesocosms, we were not able to resolve these complex processes occurring throughout the water column, and they are thus not incorporated in Figure 6. Secondly, due to the limited experimental duration and the decoupling between production and export of organic matter (especially in the high recurring treatments, see section "Temporal Decoupling of Biomass Production and POC Export"), we could not cover the whole export response period. We were therefore not able to draw quantitative conclusions for the carbon sequestration potential of artificial upwelling. It is likely that this potential would have increased, had the plankton communities had more time to adapt. For example, increased repackaging of suspended biomass into dense, fast-sinking fecal pellets could have favored POC deep export. Hence, we stress the need for monthly to seasonal experimental durations to validate the long-term potential of artificial upwelling for enhanced export production and carbon sequestration.

CONCLUSION AND OUTLOOK
Artificial upwelling had opposing effects on carbon export and potential transfer to the deep ocean. On the one hand, it enhanced the mass flux and C:N ratios of sinking matter. The latter is a key prerequisite for enhancing ocean carbon sequestration.
Biogeochemical modeling studies examining the feasibility of artificial upwelling in terms of carbon sequestration usually assume a static stoichiometry of C and N according to the Redfield-ratio (= 6.6). An important next step would be to reevaluate the findings of these studies with a higher than Redfield C:N ratio as we found in our experiment. On the other hand, artificial upwelling resulted in a shallower remineralization depth and a stronger temporal decoupling between POC production and its export. The grazer community was not able to capitalize on the enhanced primary production, and hence could not form the link between biomass production and vertical flux. Thus, artificial upwelling resulted in the export of relatively porous, slowly sinking particles, which were susceptible to remineralization and therefore of low potential for deep export and sequestration. Future work should thus focus on responses of the zooplankton community and allow it more time to react to phytoplankton growth. We furthermore found that the mode, with which upwelling is applied, is an important factor for particle properties and carbon export. We found that the postbloom export event in the singular mode produced carbon-rich particles, which were, however, respired faster and sank slower than under recurring upwelling conditions. In contrast, particle properties and the resulting RLS in the recurring mode were more favorable for POC deep export. Additionally, the continuously elevated primary and export production in the recurring mode might be more advantageous for linking the primary and export production of a slow reacting oligotrophic food web. Our assessment thus suggests that a recurring fertilization mode might be more suitable for net CO 2 removal from the atmosphere than a singular upwelling pulse. However, further research will be needed to verify if either the higher C:N ratios of sinking matter in the singular mode or the more suitable particle properties in the recurring mode are more critical for carbon dioxide removal.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in an online repository. The name of the repository and accession numbers are: PANGAEA, https://doi.pangaea.de/10.1594/PANGAEA.