Journal topic
Biogeosciences, 17, 35–45, 2020
https://doi.org/10.5194/bg-17-35-2020
Biogeosciences, 17, 35–45, 2020
https://doi.org/10.5194/bg-17-35-2020

Research article 03 Jan 2020

Research article | 03 Jan 2020

# Arctic (Svalbard islands) active and exported diatom stocks and cell health status

Arctic (Svalbard islands) active and exported diatom stocks and cell health status
Susana Agustí1, Jeffrey W. Krause2,3, Israel A. Marquez2,3, Paul Wassmann4, Svein Kristiansen4, and Carlos M. Duarte1,5 Susana Agustí et al.
• 1Red Sea Research Center, King Abdullah University of Science and Technology, Thuwal, 23955-6900, Kingdom of Saudi Arabia
• 2Dauphin Island Sea Lab, Dauphin Island, AL 36528-4603, USA
• 3Department of Marine Sciences, University of South Alabama, Mobile, AL 36688-0002, USA
• 4Department of Arctic and Marine Biology, UiT The Arctic University of Norway, 9037 Tromsø, Norway
• 5Arctic Research Centre, Department of Bioscience, Aarhus University, C.F. Mollers Alle 8, 8000 Aarhus C, Denmark

Correspondence: Susana Agustí (susana.agusti@kaust.edu.sa)

Abstract

Diatoms tend to dominate the Arctic spring phytoplankton bloom, a key event in the ecosystem including a rapid decline in surface-water pCO2. While a mass sedimentation event of diatoms at the bloom terminus is commonly observed, there are few reports on the status of diatoms' health during Arctic blooms and its possible role on sedimentary fluxes. Thus, we examine the idea that the major diatom-sinking event which occurs at the end of the regional bloom is driven by physiologically deteriorated cells. Here we quantify, using the Bottle-Net, Arctic diatom stocks below and above the photic zone and assess their cell health status. The communities were sampled around the Svalbard islands and encompassed pre- to post-bloom conditions. A mean of 24.2±6.7 % SE (standard error) of the total water column (max. 415 m) diatom standing stock was found below the photic zone, indicating significant diatom sedimentation. The fraction of living diatom cells in the photic zone averaged 59.4±6.3 % but showed the highest mean percentages (72.0 %) in stations supporting active blooms. In contrast, populations below the photic layer were dominated by dead cells (20.8±4.9 % living cells). The percentage of diatoms' standing stock found below the photic layer was negatively related to the percentage of living diatoms in the surface, indicating that healthy populations remained in the surface layer. Shipboard manipulation experiments demonstrated that (1) dead diatom cells sank faster than living cells, and (2) diatom cell mortality increased in darkness, showing an average half-life among diatom groups of 1.025±0.075 d. The results conform to a conceptual model where diatoms grow during the bloom until resources are depleted and supports a link between diatom cell health status (affected by multiple factors) and sedimentation fluxes in the Arctic. Healthy Arctic phytoplankton communities remained at the photic layer, whereas the physiologically compromised (e.g., dying) communities exported a large fraction of the biomass to the aphotic zone, fueling carbon sequestration to the mesopelagic and material to benthic ecosystems.

1 Introduction

Diatoms can support most of the Arctic primary production during the spring phytoplankton bloom (Krause et al., 2018), the key event setting the ecosystem and driving the intense carbon-uptake characteristic of the Arctic (Vaquer-Sunyer et al., 2013). However, silicic acid concentrations [Si(OH)4] are characteristically low in the European sector of the Arctic, due to the inflow of Si-depleted Atlantic water (Rey, 2012). In the Svalbard island region, Krause et al. (2018) showed diatoms to be limited by [Si(OH)4] at the spring bloom and suggested that silicon limitation could collapse a diatom bloom before nitrogen, when spring conditions favor diatoms instead of favoring the haptophyte Phaeocystis. A similar observation was made during the spring bloom in southern Greenland, whereby diatom depletion of [Si(OH)4] collapsed the bloom with ∼4µmol L−1 remaining nitrate (Krause et al., 2019).

The termination of the Arctic spring bloom is characterized by rapid sinking of diatom cells, leading to high sedimentary fluxes in the spring (Olli et al., 2002; Wassmann et al., 2006; Bauerfeind et al., 2009), precluding this production from being recycled in the upper ocean. The apparent rapid sinking of the senescent diatom bloom appears to sustain the depletion of CO2 in surface waters initiated by the bloom and drives strong atmospheric CO2 uptake (Bates and Mathis, 2009) as average pCO2 values post-bloom are typically below 300 ppm – with some values as low as 100 ppm (Takahashi et al., 2002; Holding et al., 2015).

Factors regulating diatom sedimentation have been explored for decades; however, there are few published reports on the status of diatoms' health in the Arctic during blooms and on the possible role deteriorated cell health status may play in driving sedimentary fluxes. Alou-Font et al. (2016) observed large variability in the health status of phytoplankton in the Canadian Arctic, influenced by the light and temperature conditions but not by nitrate concentration – typically thought to be the main yield-limiting nutrient. Silicon limitation has been shown to affect both autolysis (i.e., cell death) and the potential to form aggregates (which facilitate sinking) in Coscinodiscus wailesii cultures, whereas the latter was less pronounced under nitrogen limitation (Armbrecht et al., 2014).

Because of diatoms' obligate silicon requirement, its depletion in the water column would exclusively affect their physiology and, potentially, their biogeochemical fate. Lomas et al. (2019) recently demonstrated that polar diatoms have high elemental density (i.e., element content per unit biovolume) relative to low-latitude diatoms, and this is especially true for silicon content (consistent with results from Baines et al., 2010). Therefore, short-term changes in diatom physiology, e.g., responses to nutrient stress, may favor rapid sinking of polar diatoms much more than in temperate diatom species. While one could examine diatom cells from sediment traps, which is a standard approach used to explore diatoms' sinking fluxes, this methodology precludes accurate analysis of physiological health due to both the time required to collect cells (i.e., cell status can change) and trap fixatives (necessary to avoid “swimmers” from consuming sedimented material) lead to mortality of all cells. A new instrument, the Bottle-Net, has been applied to address this methodological gap. The Bottle-Net is a plankton net fitted inside a Rosette sampling system that can be used to collect plankton samples at depth without a prolonged deployment. This system was recently used to assess the stock and health status of microplankton in deep waters across the subtropical and tropical ocean (Agustí et al., 2015). Using the Bottle-Net at stations around the Svalbard islands, we examined diatom stocks within and below the photic layer and assessed their health status along contrasting stages of bloom development. We also conducted two exploratory experiments to test the hypotheses that dead diatom cells in the field sink faster than living ones, based on previous culture experiment results (Smayda, 1971), and that spring field diatoms can die rapidly upon falling below the photic layer.

2 Methods

## 2.1 Sampling and study area

The study was conducted between 17 and 29 May 2016 onboard the R/V Helmer Hanssen. The cruise started in the southwestern fjords of Svalbard islands transited northward toward Erik Eriksenstretet and then south towards stations near the polar front and the Barents Sea (Fig. 1).

Figure 1ARCEx cruise study area with the insert showing the station number and location (green dots) around the Svalbard islands.

Vertical profiles with a Seabird Electronics 911 plus CTD, provided with an oxygen sensor, fluorometer, turbidity meter and PAR sensor (Biospherical/LI-COR, SN 1060), were conducted at all stations sampled. Water samples were collected using 12 5 L Niskin bottles installed on a rosette sampler. Water samples were taken between the surface and the bottom (max. 500 m) for analysis of nutrients, diatom silica, productivity and other properties (Krause et al., 2018). We calculated the upper mixed layer (UPM) as the shallowest depth at which water density (σθ) differs from surface values by more than 0.05 kg m−3 (Mura et al., 1995).

At eight of the stations (Fig. 1), microphytoplankton samples were collected by using two Bottle-Net devices installed on the rosette sampler. The Bottle-Net is a new oceanographic device developed for the Malaspina 2010 circumnavigation expedition, which consists of a 20 µm conical plankton net housed in a cylindrical PVC pipe and is designed to be mounted in place of a Niskin bottle on the rosette sampler. The Bottle-Net cover (on top) hermetically opens and closes bottle remotely using the rosette's carousel bottle-firing mechanism, thereby initiating or terminating sample collection; the casing is open at the bottom to allow the water filtered through the internal plankton net to flow out (Agustí et al., 2015). The Bottle-Net is lowered with the top cover closed, opened at the desired bottom depth (Db, m) during the ascension of the rosette, and then the top cover is closed again at the upper depth (Du, m) of the water column to be sampled. This results in one integrated sample, from Db to Du, per deployment. Two Bottle-Nets were used mounted in the rosette sampling system, one to collect phytoplankton at the aphotic zone and the second to collect the community in the upper water column (photic layer). The two layers were selected by combining the information on light penetration (PAR sensor) and chlorophyll a fluorescence obtained during the downward CTD cast. The upper layer included the thickness of the photic layer to the depth where chlorophyll fluorescence faded away, which typically corresponded to very low levels of PAR (e.g.,  0.1 %–1 % of surface irradiance). For the aphotic zone, one Bottle-Net was remotely opened and started filtering water when the rosette reached the maximum depth at each station, and it collected until reaching the depth 10 m below the maximum depth of the photic zone. The second Bottle-Net was opened at the bottom of the photic layer and was kept open until reaching the water surface. Once on deck, the Bottle-Nets were gently rinsed with filtered seawater to retrieve the sample from the collector. Sampled volume was estimated as the product between the cross-sectional area of the mouth of the Bottle-Net and the vertical distance covered by the device from the start of the ascension to the closure of the top cover (Db to Du). The Bottle-Net presents an aspect ratio of 4, to avoid resuspension of materials filtered, displaying an efficiency of filtration of 96 % for deep tows (2000–4000 m) at towing velocities around to 30 m min−1, i.e., standard rosette retrieval velocities (Agustí et al., 2015).

## 2.2 Microplankton abundance and viability

Bottle-Net subsamples were processed to identify living and dead phytoplankton cells. The freshly collected samples were filtered onto 0.8 µm pore size black Nuclepore filters, stained with the BacLight Viability kit, placed on slides and frozen at −80C until examination under epifluorescence microscopy. Another fraction of the sample collected by the Bottle-Net was fixed with formalin for further examination at the laboratory. The observed diatoms were classified to genera. The percentage of living or dead cells relative to the total (i.e., dead plus living) was calculated for the total community and by genera.

The BacLight Viability kit (Molecular Probes, Invitrogen) is a double staining technique to test cell membrane permeability and is proven to be an effective method for determining phytoplankton viability (Llabrès and Agustí, 2008; Agustí et al., 2015). When excited with blue light under the epifluorescence microscope, living phytoplankton cells with intact membranes fluoresce green (SYTO 9, nucleic acid stain) and dead phytoplankton cells with compromised membranes fluoresce red (propidium iodine, nucleic acid stain). The samples were examined under blue light, most onboard the research vessel, using a Partec CyScope® high-power blue (470 nm) and green (528 nm) LED-illuminated epifluorescence microscope. In the laboratory at KAUST, all samples were examined using a Zeiss Axio Observer Z1 epifluorescence LED-illuminated microscope (Colibri 7 LED system). The fluorescence of the stained cells is well preserved at −80C for several months, and samples were transported frozen between the port of arrival (Tromsø, Norway) and KAUST.

## 2.3 Decay and sinking rates of living microphytoplankton cells

The expected mortality rates of living phytoplankton cells when transferred from the photic to the aphotic layer were examined at station no. 3 with vertical tows from the photic layer. An aliquot of the photic-layer microphytoplankton sample was resuspended in 2 L of 0.7 µm filtered surface water and incubated in the dark at 4 C for 7 d, simulating sedimentation from the photic layer into the aphotic layer. The community was sampled at the onset of the experiment and during set time intervals (i.e., 1, 3, 5 and 7 d). Immediately after sampling, cells were stained with the vital stain BacLight kit, then prepared and examined under an epifluorescence microscope (as described above) to quantify the proportion of living cells in the community. The half-life (i.e., the time required for the number of living cells to decline by 50 %) and the decay rate for each living-cell population were then calculated from the decline in living cells over time.

An experiment to test whether dead diatom cells sink faster than living cells was conducted shipboard using a sinking column (30 cm diameter, 1.35 m height, internal volume of 95 L). The chamber was placed on deck, filled with 20 µm filtered surface seawater, and left undisturbed for ∼1 h before starting the experiment. Microplankton collected in a vertical net tow (20 µm mesh) from the photic layer of Erik Eriksenstretet (close to the position of station no. 7) was resuspended in 1 L of 0.7 µm filtered seawater and gently added at the surface of the sinking column. A fresh subsample of the initial community, which was added to the surface of the chamber, was stained with the BacLight kit and the diatoms were examined for identification and the percentage of living or dead cells as described above. The samples at the bottom of the sinking column (sampling port located 1.35 m below the surface) were collected at intervals of time of 0 (time when the sample was added at the surface), 1, 4 and 12 h after the initial time, and were processed similarly to the initial community material.

3 Results

The stations sampled encompassed a spectrum of bloom conditions. Station 4 (off the western Svalbard shelf) waters were pre-bloom, as indicated by low diatom stocks, high dissolved inorganic nutrient concentrations (photic layer concentrations $\mathrm{Si}\left(\mathrm{OH}{\right)}_{\mathrm{4}}=\mathrm{4.15}±\mathrm{0.04}$µmol Si L−1, ${\mathrm{NO}}_{\mathrm{3}}+{\mathrm{NO}}_{\mathrm{2}}=\mathrm{9.43}±\mathrm{0.09}$µmol N L−1, Table 1) and relatively low stratification (Table 1). All other stations sampled were characterized by comparatively depleted nutrient concentrations (photic layer concentrations $\mathrm{Si}\left(\mathrm{OH}{\right)}_{\mathrm{4}}=\mathrm{0.99}±\mathrm{0.30}$µmol Si L−1, ${\mathrm{NO}}_{\mathrm{3}}+{\mathrm{NO}}_{\mathrm{2}}=\mathrm{1.93}±\mathrm{0.76}$µmol N L−1, Table 1), thereby representing communities that were either in advanced blooming stages or were senescent after blooming. Stations 6 (SW Svalbard shelf) and 8 (E Svalbard shelf) supported actively blooming diatom populations, with station 8 having the highest chlorophyll a concentration (10.5 µg Chl a L−1, as described in Krause et al., 2018), and a large fraction of living diatom cells (about 70 %, Table 1). Both locations had the highest stratification among the stations, as indicated by the low UPM values (Table 1). In contrast, station 9 (polar front) supported a senescent diatom population in post-bloom phase, as indicated by depleted nutrient pools and a low percentage of living diatom cells (46.0 %, Table 1). The highest mixing was observed at station sampled at the Barents Sea (Table 1). The percentage of living cells was not significantly correlated with the concentrations of NO3+NO2 (two-tailed test, $r=-\mathrm{0.54}$, P=0.17) or Si(OH)4 (two-tailed test, $r=-\mathrm{0.69}$, P=0.06).

Table 1Stations number and location, averaged (± SE) photic layer temperature, salinity, upper mixed layer (UPM) depth, nutrients, and measurements made with the Bottle Net (BN) in the photic and aphotic zones, indicating the depth of the tows, and the abundance and percentage of living diatoms found at the two layers.

Taxonomic classification under epifluorescence microscopy is not particularly accurate, but we were able to unambiguously identify different diatom genera and some species. The more abundant genera found in the samples were Thalassiosira spp., differentiated between large (L Thalassiosira) and small (Thalassiosira) morphotypes; Chaetoceros spp., with a large representation of Chaetoceros socialis; pennate diatoms including colonies of Fragilariopsis spp., Navicula pelagica and Pseudo-nitzschia spp.; less abundant but identifiable cells of Amphiprora hyperborean; and Coscinodiscus sp. among others.

The living (green fluorescence) and dead (red fluorescence) cells were clearly identified under the LED-illumination of the epifluorescence microscopes used (Fig. 2). The fraction of living diatom cells in the photic layer averaged 59.4±6.3 % but ranged broadly from 20.9 % in station 4, in pre-bloom state, to 72.0 % in station 5, which supported an active bloom. In contrast, the population sinking below the photic layer was comprised mostly of dead cells (20.8±4.9 % living cells, Fig. 2). Indeed, the fraction of living diatoms was consistently greater in the photic layer than in the diatom stock sinking below the photic layer (Wilcoxon ranked sign test, P=0.0078, Fig. 3), a pattern consistent across taxa found in at least four of the stations (large-celled Thalassiosira spp., P=0.02, N=4; Fragilariopsis spp., P=0.005, N=6; Chaetoceros spp., P=0.0054, N=6; Fig. 3), but the percent living cells in the photic layer and below this layer was not significantly different for the small-celled Thalassiosira spp. (P=0.09, N=6).

Figure 2Photographs of the natural Arctic diatoms sampled with the Bottle-Net observed under epifluorescence microscopy and stained with the BacLight kit. (a) Colonies of Thalassiosira spp. showing green fluorescence corresponding to living cells. (b) Colonies of Fragilariopsis spp. showing dead cells (red fluorescence, vertical-left colonies) and living cells (green fluorescence, transversal-right colony). (c) Surface layer community, composed of multiple diatom genera (Chaetoceros spp., Fragilariopsis spp., Thalassiosira spp., pennates), showing a blend of living cells (green fluorescence) and dead cells (red fluorescence). (d) Aphotic zone sample showing dead colonies (red fluorescence) of Fragilariopsis spp. and Thalassiosira spp. (two-cell colony in the bottom right of the photo).

Figure 3Box plots showing the distribution of the percentage of living diatoms in the photic (blue) and aphotic (brown) layers. Percentage of living cells among (a) the total diatom community and (b) for the populations of the most abundant diatom taxa observed during the cruise. The asterisks indicate significant differences between photic and aphotic zones (P<0.05). Boxes encompass the central 50 % of the data, the horizontal line inside the box represents the median, and vertical bars encompass 90 % of the data.

Among stations, there was significant variability in the diatom assemblage structure. Earlier cruise stations were dominated by Fragilariopsis spp. and Chaetoceros spp. This changed from station 6 to 8, where communities were dominated by Fragilariopsis spp. and Thalassiosira spp. and were the areas with the highest diatom biomass observed (station no. 8, Fig. 4). Community composition changed at the polar front and Barents Sea stations (Fig. 4) with a larger contribution of other taxa, including Navicula pelagica (station no. 9, Fig. 4). The diversity of the diatoms found at the aphotic zone differed in several stations from that found at the photic layer (Fig. 4). The large-celled Thalassiosira sp. colonies dominated the aphotic community in several stations although they were not dominant at the photic community (Fig. 4). At station no. 4, the community sampled was more diverse at the aphotic layer than at the photic layer (Fig. 4). The stock of diatoms that had sunk below the photic layer comprised, on average, 24.2±6.7 % of the total water column stock, with the proportional contribution ranging considerably among groups (Fig. 5). The proportion of biomass of the large-celled Thalassiosira spp. colonies in the aphotic layer was the largest and Chaetoceros spp. the smallest (Fig. 5). Station no. 4 (pre-bloom status) had a larger proportion of diatom biomass in the aphotic layer, and station no. 8, diatom bloom station, had the lowest. At station no. 8, however, the photic-zone population of the dominant Thalassiosira species contained 54.8 % of living cells and was paralleled with a significant contribution of dead cells at the aphotic layer (Fig. 4), suggesting the collapse of the bloom had already been initiated despite the considerable photic-layer biomass. Similarly, Fragilariopsis senescence at the photic layer of station no. 3 (only 35.1 % of cells were alive at the photic layer) helps explain its larger contribution than in the aphotic layer (Fig. 4). There was a significant negative relationship between the percent of the diatom stock population that had sunk below the photic layer and the percent of living cells in the photic layer (R2=0.39, P<0.001, Fig. 5b), indicating that healthy, actively growing populations largely remain in the surface.

Figure 4Pie charts showing the diatom community within the photic and aphotic zones among stations. Colors correspond to different taxa, and the numbers indicate the percentage of cells relative to the total community.

Figure 5(a) The proportion (mean ± SE) of the water-column stock observed in the aphotic zone (relative to the sum of the aphotic and photic zones) among the different diatom taxa. (b) The relationship between the percentage of living diatoms cells among the different populations in the photic layer and the proportion of the water-column population stock observed in the aphotic zone. The line represents the fitted linear regression (R2=0.39, P<0.001).

The suggestion that dead diatom cells sink faster than living cells was tested experimentally. Initially, only 6.7 % of the cells of the Flagilariosis spp. and Thalassiosira spp. colonies dominating the community tested were dead. However, all cells settling to the bottom of the sedimentation chamber within 1 h of the experiment start were dead, including large Coscinodiscus sp. cells (Fig. 6). The population of cells settling to the chamber bottom 4 and 12 h following addition of the fresh, healthy community, was also largely dominated by dead Flagilariosis spp. and Thalassiosira spp. colonies, 82.2 % and 71.7 %, respectively. And the fraction of living cells which had settled the height of the chamber proportionally increased with time (Fig. 6). These experimental results indicated that dead diatom cells among the groups examined sink faster than living cells.

Figure 6Diatom cell viability among the sinking cells. The initial percentage of dead cells correspond to the fresh photic-zone Arctic microplankton sample (see Methods) and added to the surface of the sinking column (1.35 m height) at time 0. The percentages of dead cells at the bottom of the sinking column were collected at intervals times of 0, 1, 4 and 12 h after sample inoculum addition.

The experiment testing diatom survival in aphotic zone indicated that once diatom cells sink below the photic layer, they would die rapidly. The incubation was performed close to the temperature below the photic layer that averaged 2.978 C ± 0.003 at station no. 3, suggesting no thermal effects. The cell concentration at the onset of the experiment was 298 cells mL−1. The half-life (i.e., percent of living cells reduced to half) survival times were remarkably uniform across diatom taxa, ranging from 0.9 d, for Thalassiosira spp. to 1.3 d for Coscinodiscus sp., depending on species (Fig. 7). Once dead, the cells lysed; half-life periods for cell death and lysis after transfer into aphotic darkness increased from 1.6 d, for the smaller Flagilariosis spp. cells, to 5.3 d for the larger Thalassiosira spp. cells (Fig. 7).

Figure 7Decay in the cell abundance of living (blue diamonds) and total (orange squares) diatoms when exposed to aphotic zone light conditions. (a) Large-celled Thalassiosira spp. (b) Fragilariopsis spp. (c) Thalassiosira spp. (d) Pennate diatom. The solid black lines and equations show the fitted linear regressions for the percent of living cells (blue box, all fitted lines significant, p<0.05) and total population cells (orange box, none of the fitted lines were significant, p>0.05).

4 Discussion

The results presented confirm that active and healthy diatom populations, as those actively growing during the spring bloom, are associated with relatively small stocks of fast-sinking diatoms. In contrast, unhealthy diatom populations, such as those present before blooming has initiated or in the senescent phase of the bloom, characterized by a large fraction of dead cells, support comparatively larger pools of sinking diatoms.

These observations are consistent with early reports, based mostly on laboratory cultures, indicating that dead diatom cells sink faster than living ones (Smayda, 1971). The experiment conducted, albeit at only one station, showed that dead cells sank much faster than living ones in a field assemblage with considerable diversity in species and in the physiological condition. Indeed, whereas the dominant populations tested were dominated by living healthy cells, only dead cells were collected at the bottom of the sedimentation chamber over the first few hours of the experiment, and the proportion of living cells collected increased over time. Moreover, our experimental assessment of diatom survival incubated at aphotic conditions suggested that once sinking below the photic layer, diatom cells could die at half-lives of 21.8 to 30.2 h across species. This result, although limited to one experiment, was consistent among the major genera and functional groups analyzed, and it reflected survival at in situ conditions. Smayda and Mitchell-Innes (1974) also reported the decrease in viable cells after darkness: “After 6 days of dark incubation, the number of viable cells of Chaetoceros curvisetus recognizable decreased from 760 to 240 cells per ml.”, representing a decay rate of 0.19 d−1 (i.e., 50 % loss of cells in 3.6 d) comparable to the rate reported here. Other studies also reported rates of living cell mortality in darkness close to those found here (Segovia and Berges, 2009; Agustí et al., 2015). The decay rates calculated for living or viable vegetative cells were more than 3 times faster than those observed for the total cells in the population. These included both viable and non-viable vegetative cells, which are, however, morphologically similar and could not be differentiated unless using specific methods to discriminate living from dead cells, such as the staining test used here. Phytoplankton vegetative cells do not survive under darkness (Smayda and Mitchell-Innes, 1974; Segovia and Berges, 2009; Segovia et al., 2003), and only resting spores and resting cells are able to survive in the dark (Ignatiades and Smayda, 1970; Smayda and Mitchell-Innes, 1974; Peters and Thomas, 1996). Our results reporting fast diatom cell death under aphotic conditions are contrasting with the expectation of high survival capacity of polar diatoms to darkness supported by existing evidence. Recent reports identified fast photosynthetic response to irradiance in diatoms sampled during the dark wintertime around the Svalbard islands (Kvernvik et al., 2018). Phenotypic selection for specific physiological properties allows polar diatoms to acclimate to low light and darkness (Lacour et al., 2019). Our experiment, however, was carried out in a spring community under a 24 h light: 0 h darkness photoperiod and with shallow mixing depths. The community tested was, therefore, likely photo-acclimated to a prolonged photoperiod and relatively high irradiance. Therefore, it is expected to respond differently under darkness than Arctic microalgae growing underneath the ice or under very short photoperiods and/or minimal irradiance levels (Lacour et al., 2019; Berge et al., 2015). Hence, Arctic phytoplankton are expected to show contrasting responses to prolonged darkness in the spring, when acclimated to longer photoperiods than in winter in which surviving cells are expected to be acclimated to short photoperiods.

Beyond the stressor of continuous darkness, the fast decay rates observed here under aphotic experimental conditions could also have been influenced, in addition to darkness, by the low nutrient availability and/or the possible presence of pathogens or parasites. This suggests that survival of natural populations below the photic layer may be lower than expected in the dark from axenic, high-nutrients in vitro studies with cultures, an aspect already pointed by Ignatiades and Smayda (1970). The [Si(OH)4] at station no. 3 was 0.75 µmol Si L−1 and [NO3+NO2] was 1.79 µmol N L−1; Lomas et al. (2019) demonstrated that polar diatoms have higher Si:N requirements (>1.5 mol : mol) than lower-latitude diatoms; thus nutrient data would predict the yield to be limited by Si relative to N. As the decline rates were derived from a single test, further experiments on cell decay rates of Arctic diatoms sampled in the spring under dark conditions will be required to confirm our results.

Cell abundance and health state observed were consistent with previous studies. In the Canadian Arctic, living cells in open water and ice-covered stations represented 57.3±5.8 % and 48.0±3.9 % (± SE), respectively (Alou-Font et al., 2016), which are similar to the proportions in our study. The percentage of living cells was higher during the bloom periods than the periods before and after (Alou-Font et al., 2016). This result appeared to be driven by light and low nutrient concentrations (Alou-Font et al., 2016). The quantification of the percent of living cells in our study helped to identify the different stages of the Arctic spring bloom among the stations sampled. A pre-bloom situation with low cell abundance and a small percentage of living cells was found at station no. 4 located further west of Svalbard islands, where silicic acid and nitrogen concentrations were higher and mixing was more significant than in other Arctic stations. The healthiest diatom community was observed at station no. 5, where the high stratification and Si(OH)4 concentration above the half saturation constant (KS) of 2 µM (from kinetic experiments in the same region by Krause et al., 2018) helped the diatoms to grow actively. The highest cell abundance was observed at station no. 8, but this persisted with a lower percentage of living diatoms, and the Si(OH)4 concentration well below the KS value indicated that the bloom was reaching the maximum capacity, although diatom sinking was still low. A post-bloom situation was identified at the polar front community, with similar percentages of living cells at the photic and aphotic zones as a result of high sinking induced by Si and nitrogen limitation, as suggested by the lower Si(OH)4 KS of 0.8 µM (Krause et al., 2018). The diatom community captured by the Bottle Net below the photic layer was consistent with the limited, but comparable, data obtained by sediment traps deployed in the area which also indicated Fragilariopsis and Thalassiosira species to be the dominant contributors to Si and biomass export (Krause et al., 2018).

Given the range of bloom development represented among our stations, the results presented conform to a conceptual model in which nutrients, including Si (Rey, 2012; Krause et al., 2018), and a mixed layer drive the growth of diatoms during the Arctic spring bloom (Wassmann et al., 1997; Reigstad et al., 2002). For diatoms, Si depletion results in two potential physiological issues: yield limitation (i.e., diatom standing stock is too high to be supported by the available silicic acid) and intense kinetic and growth limitation (i.e., depleted silicic acid limits diatom Si uptake to such a degree that growth must slow, Krause et al., 2018). Thus, such a situation would stimulate mass sedimentation, suggested to be an evolutionary adaptation to help diatom communities persist when nutrients are limiting (Raven and Waite, 2004). A large fraction (30 %–50 %) of the diatom cells in the silicon-depleted photic layer were dead, pointing at acute silicic acid limitation as the likely factor triggering partial mortality, while the remaining cells are likely to be senescent. Unhealthy diatoms would then lose the capacity to actively regulate buoyancy that characterizes healthy diatom cells (Smayda, 1970), leading to rapid sinking of the bloom. Acute silicic acid limitation is identified, therefore, as the event leading to loss of the capacity to actively regulate buoyancy that characterizes diatom cells (Smayda, 1970) and rapid sinking of the bloom. The potential for rapid sedimentation is enhanced by higher silica quotas for polar diatoms (Lomas et al., 2019) compared to lower-latitude diatoms (Brzezinski, 1985).

Diatoms have been shown to have a remarkable metabolic capacity to regulate buoyancy, both maintaining zero and positive buoyancy (e.g., Villareal et al., 2014) involving regulation through the production of osmolytes (Gradmann and Boyd, 2002), which plays an important role in exploiting nutrient patchiness within the photic layer (Villareal et al., 2014). Diatom sinking rates are inversely related to growth rate, so that silicon depletion is expected to result in increased sinking rates, despite field diatoms reducing their silica per cell when kinetically limited by silicic acid (McNair et al., 2018). There is experimental demonstration that silicon depletion plays the most important role, compared to nitrogen or phosphorus, in triggering rapid sinking of diatom cells, indicating that biochemical aspects of silicon metabolism are particularly important to diatom buoyancy regulation (Bienfang et al., 1982). N : P ratios in this region do not suggest that phosphorus plays a limiting role in primary production, and when silicic acid is depleted enough nitrate remains to fuel growth of other phytoplankton groups (e.g., Phaeocystis, Krause et al., 2018). Once diatoms lose their capacity to regulate buoyancy and sink below the photic layer, they die rapidly and are unable to ascend back to the photic layer, resulting in the rapid sinking fluxes that drives high sedimentation rates characteristic of the termination of the Arctic spring bloom (Olli et al., 2002; Wassmann et al., 2006; Bauerfeind et al., 2009). Rapid sinking of the Arctic spring bloom, in turn, precludes carbon recycling in the photic layer, thereby leading to undersaturated pCO2 driving the large atmospheric CO2-uptake characteristic of the European sector of the Arctic during this season which does not equilibrate with the atmosphere until months later (Bates et al., 2009; Takahashi et al., 2002; Holding et al., 2015).

A large fraction of the total water column phytoplankton biomass was observed below the photic layer, representing on average 24 % ± 6.7 (± SE) of the surface diatom populations in the study area. This considerable proportion can be explained by high diatom export from the photic zone, as opposed to lateral advection. This is consistent with the high rates of biogenic silica (proxy for diatom biomass) export at stations 4, 7–8 and 10; rates were a factor of 4 higher than integrated diatom silica production in the upper water column and represented up to 40 % of the integrated diatom silica standing stock (Krause et al., 2018). These cruise trends are in agreement with the observation of large sinking events in the Arctic as reported for ice diatoms (Boetius et al., 2013; Aumack et al., 2014) associated with ice melting in the Arctic, which represent a large carbon supply to benthic communities in the Arctic shelves (Moran et al., 2005; Tamelander et al., 2006). While we do not report data for an ice-diatom assemblage, data for the same cruise showed that silicon-uptake rates of ice diatoms near stations nos. 7 and 8 were strongly limited by Si(OH)4 concentration in the surface waters, likely limiting their growth to a degree (Krause et al., 2018), and previous studies have noted their susceptibility to silicon limitation (Cota et al., 1990; Smith et al., 1990). Our results show that healthy phytoplankton communities remained at the photic layer, although dying communities exported a large fraction of the biomass (up to 65 %) to the aphotic zone.

In summary, the results presented here support a link between diatom cell health status, likely driven by progressive nutrient limitation, and sedimentation fluxes in the Arctic. Whereas the link between diatom health status and sinking rates has long been established (Smayda, 1971), the evidence corresponded to algal cultures in the laboratory. This is the first demonstration of these ideas using natural diatom communities in this region, partially due to the logistical challenges of assessing both viability and settling in the field. Our conceptual model suggests that deterioration of diatom health, such as that occurring when acute silicon limitation or other resource limitations along the spring bloom are reached, leads to loss of the capacity to regulate buoyancy and leads to rapid sinking, with cells exported below the photic layer. Understanding the role of cell health status, and the role of silicic acid depletion, in the regulation of diatom sinking rates is fundamental to mechanistically understand the biological pump in the Arctic and its response to future changes.

Data availability
Data availability.

Data are available upon request to the authors.

Author contributions
Author contributions.

SA and CMD conceived and designed the study. SA, CMD, IAM, JWK, PW and SK conducted analysis. All co-authors contributed to the writing of the paper, led by SA.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

This research was supported by King Abdullah University of Science and Technology through baseline funding BAS/1/1072-01-01 and BAS/1/1071-01-01 to Susana Agustí and Carlos M. Duarte, respectively. The ARCEx project funded by industry partners and the Research Council of Norway (project no. 228107) provided funding to Paul Wassmann, and the Dauphin Island Sea Lab funded Jeffrey W. Krause. We thank the science team and crew of the R/V Helmer Hanssen.

Financial support
Financial support.

This research has been supported by the KAUST (grant nos. BAS/1/1072-01-01 and BAS/1/1071-01-01), the Research Council of Norway (grant no. 228107) and the Dauphin Island Sea Lab (internal funding).

Review statement
Review statement.

This paper was edited by Koji Suzuki and reviewed by three anonymous referees.

References

Agustí, S., González-Gordillo, J. I., Vaqué, D., Estrada, M., Cerezo, M. I., Salazar, G., Gasol, J. M., and Duarte, C. M.: Ubiquitous healthy diatoms in the deep sea confirms deep carbon injection by the biological pump, Nat. Commun., 6, 7608, https://doi.org/10.1038/ncomms8608, 2015.

Alou-Font, E., Roy, S., Agustí, S., and Gosselin, M.: Cell viability, pigments and photosynthetic performance of Arctic phytoplankton in contrasting ice-covered and open-water conditions during the spring-summer transition, Mar. Ecol.-Progr. Ser., 3, 89–106, 2016.

Armbrecht, L. H., Smetacek, V., Assmy, P., and Klaas, C.: Cell death and aggregate formation in the giant diatom Coscinodiscus wailesii (Gran & Angst, 1931), J. Exp. Mar. Biol. Ecol., 452, 31–39, 2014.

Aumack, C. F., Juhl, A. R., and Krembs, C.: Diatom vertical migration within land-fast Arctic sea ice, J. Marine Syst., 139, 496–504, 2014.

Baines, S. B., Twining, B. S., Brzezinski, M. A., Nelson, D. M., and Fisher, N. S.: Causes and biogeochemical implications of regional differences in silicification of marine diatoms, Global Biogeochem. Cy., 24, GB4031, https://doi.org/10.1029/2010GB003856, 2010.

Bates, N. R. and Mathis, J. T.: The Arctic Ocean marine carbon cycle: evaluation of air–sea CO2 exchanges, ocean acidification impacts and potential feedbacks, Biogeosciences, 6, 2433–2459, https://doi.org/10.5194/bg-6-2433-2009, 2009.

Bauerfeind, E., Nöthig, E.-M., Beszczynska, A., Fahl, K., Kaleschke, L., Kreker, K., Klages, M., Soltwedel, T., Lorenzen, C., and Wegner, J.: Variations in vertical particle flux in the Eastern Fram Strait (79 N/4 E) during 2000–2005. Results from the Deep-Sea Long-Term observatory HAUSGARTEN, Deep-Sea Res. Pt. I, 56, 1471–1487, 2009.

Berge, J., Renaud, P. E., Darnis, G., Cottier, F., Last, K., Gabrielsen, T. M., Johnsen, G., Seuthe, L., Weslawski, J. M., Leu, E., and Moline, M.: In the dark: a review of ecosystem processes during the Arctic polar night, Progr. Oceanogr., 139, 258–271, 2015.

Bienfang, P. K., Harrison, P. J., and Quarmby, L. M.: Sinking rate response to depletion of nitrate, phosphate and silicate in four marine diatoms, Mar, Biol,, 67, 295–302, 1982.

Boetius, A., Albrecht, S., Bakker, K., Bienhold, C., Felden, J., Fernández-Méndez, M., Hendricks, S., Katlein, C., Lalande, C., Krumpen, T., and Nicolaus, M.: Export of algal biomass from the melting Arctic sea ice, Science, 339, 1430–1432, 2013.

Brzezinski, M. A.: The Si : C : N ratio of marine diatoms: Interspecific variability and the effect of some environmental variables, J. Phycol., 21, 347–357, 1985.

Cota, G. F., Anning, J. L., Harris, L. R., Harrison, W. G., and Smith, R. E.: Impact of ice algae on inorganic nutrients in seawater and sea ice in Barrow Strait, NWT, Canada, during spring, Can. J. Fish. Aquat. Sci., 47, 1402–1415, 1990.

Gradmann, D. and Boyd, C. M.: Impact of osmolytes on buoyancy of marine phytoplankton, Mar. Biol., 141, 605–618, 2002.

Holding, J., Duarte, C. M., Sanz-Martín, M., Mesa, E., Arrieta, J. M., Chierici, M., Hendriks, I., García-Corral, L. S., Regaudie-de-Gioux, A., Delgado, A., Reigstad, M., Wassmann, P., and Agustí, S.: Temperature dependence of CO2-enhanced primary production in the European Arctic Ocean, Nat. Clim. Change, 5, 1079–1082, 2015.

Ignatiades, L. and Smayda, T. J.: Autecological studies on the marine diatom Rhizosolenia fragilissima bergon. ii. enrichment and dark viability experiments 2, J. Phycol., 6, 357–364, 1970.

Krause, J. W., Duarte, C. M., Marquez, I. A., Assmy, P., Fernández-Méndez, M., Wiedmann, I., Wassmann, P., Kristiansen, S., and Agustí, S.: Biogenic silica production and diatom dynamics in the Svalbard region during spring, Biogeosciences, 15, 6503–6517, https://doi.org/10.5194/bg-15-6503-2018, 2018.

Krause, J. W., Schulz, I. K, Rowe, K. A., Dobbins, W., Winding, M., Sejr, M., Duarte, C. M., and Agustí, S.: Silicon limitation drives bloom termination and carbon sequestration in an Arctic bloom, Sci. Rep.-UK, 9, 8149, https://doi.org/10.1038/s41598-019-44587-4, 2019.

Kvernvik, A. C., Hoppe, C. J. M., Lawrenz, E., Prášil, O., Greenacre, M., Wiktor, J. M., and Leu, E.: Fast reactivation of photosynthesis in arctic phytoplankton during the polar night, J. Phycol., 54, 461–470, 2018.

Lacour, T., Morin, P. I., Sciandra, T., Donaher, N., Campbell, D. A., Ferland, J., and Babin, M.: Decoupling light harvesting, electron transport and carbon fixation during prolonged darkness supports rapid recovery upon re-illumination in the Arctic diatom Chaetoceros neogracilis, Polar Biol., 1–13, 2019.

Llabrés, M. and Agustí, S.: Extending the cell digestion assay to quantify dead phytoplankton cells in cold and polar Waters, Limnol. Oceanogr.-Meth., 6, 659–666, 2008.

Lomas, M. W., Baer, S. E., Acton, S., and Krause, J. W.: Pumped up by the Cold: Elemental Quotas and Stoichiometry of Polar Diatoms, Front. Mar. Sci., 6, 286, https://doi.org/10.3389/fmars.2019.00286, 2019.

McNair, H. M., Brzezinski, M. A., and Krause, J. W.: Diatom populations in an upwelling environment decrease silica content to avoid growth limitation, Environ. Microbiol., 20, 4184–4193, 2018.

Moran, S. B., Kelly, R. P., Hagstrom, K., Smith, J. N., Grebmeier, J. M., Cooper, L. W., Cota, G. F., Walsh, J. J., Bates, N. R., Hansell, D. A., and Maslowski, W.: Seasonal changes in POC export flux in the Chukchi Sea and implications for water column-benthic coupling in Arctic shelves, Deep-Sea Res. Pt. II, 52, 3427–3451, 2005.

Mura, M. P., Satta, M. P., and Agustí, S.: Water-mass influences on summer Antarctic phytoplankton biomass and community structure, Polar Biol., 15, 15–20, 1995.

Olli, K., Wassmann, P., Ratkova, T., Arashkevich, E., and Pasternak, A.: Seasonal variation in vertical flux of biogenic matter in the marginal ice zone and the central Barents Sea, J. Marine Sys., 38, 189–204, 2002.

Peters, E. and Thomas, D. N.: Prolonged darkness and diatom mortality I: Marine Antarctic species, J. Exp. Mar. Biol. Ecol., 207, 25–41, 1996.

Raven, J. A. and Waite, A. M.: The evolution of silicification in diatoms: inescapable sinking and sinking as escape?, New Phytol., 162, 45–61, 2004.

Reigstad, M., Wassmann, P., Riser, C. W., Øygarden, S., and Rey, F.: Variations in hydrography, nutrients and chlorophyll a in the marginal ice-zone and the central Barents Sea, J. Marine Syst., 38, 9–29, 2002.

Rey, F.: Declining silicate concentrations in the Norwegian and Barents Seas, ICES J. Mar. Sci., 69, 208–212, 2012.

Segovia, M. and Berges, J. A.: Inhibition of caspase-like activities prevents the appearance of reactive oxygen species and dark-induced apoptosis in the unicellular chlorophyte Dunaliella tertiolecta, J. Phycol., 45, 1116–1126, https://doi.org/10.1111/j.1529-8817.2009.00733.x, 2009.

Segovia, M., Haramaty, L., Berges, J. A., and Falkowski, P. G.: Cell Death in the Unicellular Chlorophyte Dunaliella tertiolecta. A Hypothesis on the Evolution of Apoptosis in Higher Plants and Metazoans, Plant Physiol., 132, 9–105, https://doi.org/10.1104/pp.102.017129, 2003.

Smayda, T. J.: The suspension and sinking of phytoplankton in the sea, Oceanogr. Mar. Biol., 8, 353–414, 1970.

Smayda, T. J.: Normal and accelerated sinking of phytoplankton in the sea, Mar. Geol., 11, 105–122, 1971.

Smayda, T. J. and Mitchell-Innes, B.: Dark survival of autotrophic, planktonic marine diatoms, Mar. Biol., 25, 195–202, 1974.

Smith, R. E., Harrison, W. G., Harris, L. R., and Herman, A. W.: Vertical fine structure of particulate matter and nutrients in sea ice of the high Arctic, Can. J. Fish. Aquat. Sci., 47, 1348–1355, 1990.

Takahashi, T., Sutherland, S. C., Sweeney, C., Poisson, A., Metzl, N., Tilbrook, B., Bates, N., Wanninkhof, R., Feely, R. A., Sabine, C., and Olafsson, J.: Global sea–air CO2 flux based on climatological surface ocean pCO2, and seasonal biological and temperature effects, Deep-Sea Res. Pt. II, 49, 1601–1622, 2002.

Tamelander, T., Renaud, P. E., Hop, H., Carroll, M. L., Ambrose Jr., W. G., and Hobson, K. A.: Trophic relationships and pelagic–benthic coupling during summer in the Barents Sea Marginal Ice Zone, revealed by stable carbon and nitrogen isotope measurements, Mar. Ecol.-Prog. Ser., 310, 33–46, 2006.

Vaquer-Sunyer, R., Duarte, C. M., Holding, J., Regaudie-de-Gioux, A., García-Corral, L. S., Reigstad, M., and Wassmann, P.: Seasonal patterns in Arctic planktonic metabolism (Fram Strait – Svalbard region), Biogeosciences, 10, 1451–1469, https://doi.org/10.5194/bg-10-1451-2013, 2013.

Villareal, T. A., Pilskaln, C. H., Montoya, J. P., and Dennett, M.: Upward nitrate transport by phytoplankton in oceanic waters: balancing nutrient budgets in oligotrophic seas, PeerJ, 2, e302, https://doi.org/10.7717/peerj.302, 2014.

Wassmann, P., Egge, J. K., Reigstad, M., and Aksnes, D. L.: Influence of dissolved silicate on vertical flux of particulate biogenic matter, Mar. Pollut. Bull., 33, 10–21, 1997.

Wassmann, P., Slagstad, D., Rise, C. W., and Reigstad, M.: Modelling the ecosystem dynamics of the Barents Sea including the marginal ice zone II. Carbon flux and interannual variability, J. Marine Syst., 59, 1–24, 2006.