Unprecedented insights into extents of biological responses to physical forcing in an Arctic sub-mesoscale filament by combining high-resolution measurement approaches

doi:10.21203/rs.3.rs-3376409/v1

Download PDF

Article

Unprecedented insights into extents of biological responses to physical forcing in an Arctic sub-mesoscale filament by combining high-resolution measurement approaches

https://doi.org/10.21203/rs.3.rs-3376409/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 08 Apr, 2024

Read the published version in Scientific Reports →

You are reading this latest preprint version

In Fram Strait, we combined underway-sampling using the remote-controlled Automated Filtration System for Marine Microbes (AUTOFIM) with CTD-sampling for eDNA analyses, and with high-resolution optical measurements in an unprecedented approach to determine variability in plankton composition in response to physical forcing in a sub-mesoscale filament. We determined plankton composition and biomass near the surface with a horizontal resolution of ~ 2 km, and addressed vertical variability at five selected sites. Inside and near the filament, plankton composition was tightly linked to the hydrological dynamics related to the presence of sea ice. The comprehensive data set indicates that sea-ice melt related stratification near the surface inside the sub-mesoscale filament resulted in increased sequence abundances of sea ice-associated diatoms and zooplankton near the surface. In analogy to the physical data set, the underway eDNA data, complemented with highly sampled phytoplankton pigment data suggest a corridor of 7 km along the filament with enhanced photosynthetic biomass and sequence abundances of sea-ice associated plankton. Thus, based on our data we extrapolated an area of 350 km² in Fram Strait with enhanced plankton abundances, possibly leading to enhanced POC export in an area that is around a magnitude larger than the visible streak of sea-ice.

Biological sciences/Ecology/Biooceanography

Biological sciences/Ecology/Ecosystem ecology

Earth and environmental sciences/Ocean sciences

Earth and environmental sciences/Ecology/Community ecology

Earth and environmental sciences/Ecology/Molecular ecology

The community structure in the marine realm is highly variable over temporal and geographical scales, related to variability in physicochemical drivers (e.g. temperature, salinity, nutrient availability, sea-ice coverage) and biological processes (phenological responses, advection and migration, population dynamics). Mesoscale hydrographic features (horizontal scale of 10–100 km) have been shown to impact biological and biogeochemical processes in various ways, including the modulation of nutrient supply to the euphotic zone leading to increased primary production [1], and community turnover due to vertical and horizontal transport [2, 3, 4]. Spatial and temporal turnover in species composition and biomass structure can affect food webs via losses and gains of trophic interactions, or by influencing the magnitude of existing interactions [3, 4, 5, 6], leading to high spatio-temporal variability in biogeochemical processes and carbon fluxes. Observations by satellite remote sensing, autonomous profiling floats, and towed instruments, in combination with improved modeling capacities highlighted the ubiquity and ecological relevance of oceanographic features at an even finer scale, the sub-mesoscale (0.1 -10km) scale [1, 7, 8, 9]. Such sub-mesoscale structures significantly differ in their physical properties from mesoscale structures and in terms of dynamics and impact. They often appear between mesoscale eddies [1, 10] or in the presence of strong horizontal gradients, for instance in the marginal ice zone (MIZ) [11].

There is growing interest in understanding the role of sub-mesoscale structures in shaping patchiness and variability of ecological networks and biogeochemical processes in order to better understand how ecosystems function, how their functionality varies with space and time, and how they respond to future environmental scenarios. Satellite observations revealed increasing chlorophyll a (Chl a) concentrations as a function of physical dynamics and intake of nutrients in sub-mesoscale fronts [12, 13]. Due to their small spatial and short temporal (hours to weeks) scales, in situ measurements, especially those addressing their impact on biological processes, are, however, difficult. Plankton biodiversity studies still largely depend on point sampling by CTD rosettes or net hauls from research vessels, with distances between stations usually too coarse to resolve the impact of sub-mesoscale features on small-scale plankton variability. Thus, our understanding of the impact of such small-scale dynamics on plankton composition and interactions is currently very limited. Sophisticated automated underway sampling-technology like AUTOFIM installed on board research vessels [14] can augment the current constraints with adequate horizontal resolution. Combining such technology with molecular plankton biodiversity studies based on 18S metabarcoding provides high-resolution information on plankton dynamics at various spatio-temporal scales [15].

Understanding linkages between sub-mesoscale features and ecosystem functionality is particularly important for marine areas, such as Fram Strait, that are frequently impacted by a highly dynamic and turbulent regime of mesoscale features [16, 11]. In eastern Fram Strait, warm saline Atlantic water is transported via the West Spitsbergen Current (WSC) northward, while cold and rather fresh Polar water is transported southward in western Fram Strait via the East Greenland current (EGC). A highly dynamic oceanographic regime in the zone between the two currents in combination with the semi-permanent sea-ice edge [11] and large horizontal density gradients observed in the marginal ice zone (MIZ) in the northern Fram Strait enhance the development of sub-mesoscale features in this area [11].

In July 2017 a sub-mesoscale filament was detected on satellite radar images due to the occurrence of a nearly straight 50 km long and 500 m wide streak of sea ice in the MIZ of Fram Strait. The filament was formed by a strong salinity gradient between Atlantic and Polar waters close to the ice edge. High-resolution physical measurements revealed a frontal system of two strong currents flowing in opposite directions along the filament, leading to a horizontal inward-flow from both sides and mixing mostly in the upper 100 m and subduction inside the filament [11]. In parallel to the physical measurements, we used a combination of automated underway-sampling and point sampling to collect samples for 18S meta-barcoding and optical data, addressing horizontal and vertical differences in the plankton composition within the Arctic sub-mesoscale filament. The underway-sampling was accomplished with a resolution of ~ 2 km, which is adequate to resolve sub-mesoscale linkages between the physical environment and plankton composition [11]. This high-resolution information on spatial variability in plankton composition was complemented by a high-resolution information on total Chl a-biomass and the contribution of major phytoplankton groups obtained from hyperspectral underway spectrophotometry [17]. Additionally, data from sampling with a CTD-rosette and an optical zooplankton recorder, the LOKI (Lightframe On-sight Key species Investigation system) [18] provided depth-resolved information on vertical plankton distribution in vicinity of the filament.

Description of the sub-mesoscale filament

In late July 2017 a ~ 7 km wide sub-mesoscale filament occurred under a 50 km long and only 500 m wide nearly straight streak of sea ice at 2.5°E/79°N in the marginal ice zone of Fram Strait. It was characterized by a thin surface meltwater layer above a layer of denser water, which extended to > 250 m depth [11]. Outside the filament, meltwater occupied the upper 15 m of the water column, leading to a strong stratification of the surface ocean. Here, Polar water was located directly below the meltwater layer at 20–40 m depth, while Atlantic water occupied the deeper water layers below 50 m. Inside the filament, the meltwater layer was thicker and occupied the top ~ 25 m, below which Atlantic water was found. High-resolution physical measurements revealed a frontal system of two strong currents flowing in opposite directions along the filament, leading to a horizontal inward-flow from both sides and mixing mostly in the upper 100 m inside the filament. Furthermore, the authors hypothesized that denser waters at ca. 100 m in the filament were in the process of subduction [11].

Sampling

The samples were collected in July 2017 during the RV Polarstern cruise PS107 in Fram Strait (Fig. 1). To address horizontal variability in plankton biodiversity samples were collected on one hand underway on five transects with a spatial resolution of ~ 2.5 km (1.5 nm) at a depth of ~ 10 m with the Automated filtration device for marine microorganisms AUTOFIM [15], permanently installed on RV Polarstern. Additionally at five selected sites across the sub-mesoscale filament, samples from deeper water layers were taken with a rosette sampler equipped with 24 Niskin bottles (12 L per bottle) and sensors for Chl a fluorescence, temperature and salinity (CTD). Samples collected with the rosette were taken during the up-casts at 10, 20–30, 50, 100, 200 and 400 m depth. Subsamples of 2 L were transferred from the Niskin bottles into PVC bottles. Particulate organic matter for molecular analyses was collected by sequential filtration of each water sample through three mesh sizes (10 µm, 3 µm, 0.4 µm) on 45 mm diameter Isopore Membrane Filters at 200 mbar using a Millipore Sterifil filtration system (Millipore, USA). Here, two liters of seawater were collected and filtered on a filter with 0.4 µm pore size at 200 mbar. Subsequent to filtration, particulates are stored at -80°C until further processing in the laboratory.

Environmental parameters: Data for temperature and salinity are available via PANGAEA (doi.org/10.1594/PANGAEA.894189; doi.org/10.1594/PANGAEA.889535). Samples with a volume of ~ 50 ml were taken in parallel to AUTOFIM sampling from the ship's pump system and from the CTD. The samples were directly frozen at -80°C for subsequent nutrient analyses in the laboratory. Nutrients were analyzed with an Alliance Evolution III continuous flow autoanalyzer (Alliance Instruments GmbH, Freilassing, Germany). The water samples were measured unfiltered. Measurements were made simultaneously on five channels: Phosphate [19], silicate [20], nitrite, nitrate [21] and ammonium [22]. All measurements were calibrated with a five nutrient standard cocktail (All from Merck, traceable to SRM from NIST) diluted in artificial seawater (ASW), and ASW was used as wash-water between the samples. Data were all standardized by the same in-house reference material obtained from CTD water bottles. Each run we checked our standards with Reference Material for Nutrients in Seawater (CRM 7602-a + CRM 7603a) produced by NMIJ. Our standards and methods have been proven by inter-calibration exercises like ICES and Quasimeme. Nutrient concentrations are provided in the supplements (S1).

Total Chl a biomass and fractions of phytoplankton groups: Chl a concentrations and other phytoplankton pigments were determined using the high-pressure liquid chromatography (HPLC) on samples collected in parallel to the samples collected for eDNA analyses [17, 23]. These data were also used to develop an algorithm for obtaining a continuous dataset along the ship transect (at ~ 250 m resolution) of the concentration of Chl a and major phytoplankton pigments from underway hyperspectral spectrophotometry [17, 24]. For this study, we applied the diagnostic pigment analysis developed previously by Vidussi [25] and further refined by Uitz [26] to derive Chl a concentrations (conc.) of fucoxanthin, 19-hexanoyl fucoxanthin, 19-butanoyl fucoxanthin, peridinin, zeaxanthin, and chlorophyll b from HPLC samples. of the seven major phytoplankton groups that characterize the global ocean: Diatoms, prymnesiophytes, dinoflagellates, chlorophytes (which include the prasinophyte group), chrysophytes, cryptophytes, and prokaryotic microbes. We used the derived data set to calculate for each group its fraction on the total phytoplankton Chl a conc. (total Chl a). Looking at the phytoplankton group Chl a data considered within this study from the HPLC, diatoms contribute between 39%-70%, prymnesiophytes between 0–18%, dinoflagellates between 0–31%, chlorophytes between 5 to 46% to the total Chl a, while chrysophytes and prokaryotic phytoplankton have marginal contributions (0–4% and 1 to 2%). In order to derive a high resolution surface Chl a conc. data of the four major phytoplankton groups (diatoms, dinoflagellates, prymnesiophytes and chlorophytes) identified in the HPLC data set, we applied also the same diagnostic pigment analysis as for the HPLC data to the spectrophotometric pigment data set. For the three other phytoplankton groups no diagnostic pigments could be retrieved from the underway spectroscopy data, however, there contribution was found to be minimal for the Chl a in our study area (0–5%). We calculated the fraction of each phytoplankton group in respect to its Chl a compared to the total Chl a for both data sets.

Chl a conc. data for all phytoplankton (and for major phytoplankton groups) were obtained from Sentinel-3A sensor OLCI (downloaded in full resolution (300m) from EUMETSAT at https://www.eumetsat.it/new-release-l2-reprocessed-dataset/; EUMETSAT 2022) within the area and timing of the study and used to upscale Chl a point measurements from HPLC water sample analysis and underway spectrophotometry.

LOKI-casts and image analyses

LOKI casts were conducted from 400m to the surface at each of the five stations, that were also sampled via CTD for eDNA analyses. The system continuously takes images (max. 19 frames sec-1) during the up-cast from mesozooplankton organisms that are concentrated by a 150µm plankton net, leading to a flow-through chamber with a 6.1 mp camera (for a detailed description of LOKI see [18]). All images were loaded into the LOKI browser, a software that assigns optical parameters (hue factors, gray scale, skewness etc.) to each image and links the images to the respective metadata. Then, the image quality was enhanced and double takes of objects were removed using the software ZooMi. Subsequently, the images were uploaded to the EcoTaxa website (ecotaxa.obs-vlfr.fr), an application that facilitates the annotation, i.e. assignment of taxonomic categories, to the organisms presented on the images. Due to the high resolution of the LOKI images, it was often possible to identify families and sometimes-even species, thus tackling their fine scale distribution in the water column.

DNA-isolation

Isolation of genomic DNA from the field samples was carried out using the NucleoSpin Plant Kit (Machery-Nagel, Germany) following the manufacturer’s protocol. Samples were processed around one to two weeks subsequent to sampling. The resulting DNA-extracts were stored at -20°C.

Illumina-Sequencing 18S rDNA

For Illumina-Sequencing, a fragment of the 18S rDNA containing the hypervariable V4 region was amplified with the primer set 528iF(GCGGTAATTCCAGCTCC) and 964iR(ACTTTCGTTCTTGATYRR) [15]. All PCRs (polymerase chain reaction) had a final volume of 50 µL and contained 0.02 U Phusion Polymerase (Thermo Fisher, Germany), the 10-fold polymerase buffer according to manufacturer’s specification, 0.8 mM (each) dNTP (Eppendorf, Germany), 0.2 µmol L^− 1 of each Primer and 1µL of template DNA. PCR amplification was performed in a thermal cycler (Eppendorf, Germany) with an initial denaturation (94°C, 2 min) followed by 35 cycles of denaturation (94°C, 20 sec), annealing (58°C, 30 sec), and extension (68°C, 30 sec) with a single final extension (68°C, 10 min). The PCR products were purified from an agarose gel 1% [w/v] with the NucleoSpin Gel Kit (Machery-Nagel, Germany) and Minelute PCR Purification kit (Qiagen, Germany). Subsequent to purification of the 18S rDNA fragment the DNA concentration of the samples was determined using the Quantus Fluorometer (Promega, USA). Prior to the library preparation, the DNA fragments were diluted with TE-buffer to a concentration of 0.2 ng/µL. The library preparation was based on the 16S metagenomic protocol of Illumina (Illumina, USA). Finally, sequencing of the DNA-fragments was carried out using a MiSeq-Sequencer (Illumina, USA). Raw sequences had an approximate length of ~ 200 bp. Sequences generated in this study have been deposited via GfBio [27] in the European Nucleotide Archive (ENA) with the accession number PRJEB66268.

Sequence analyses and annotation

Further bioinformatic processing of the raw sequences was done with the dada2 package v.1.18 in R v.4.0 [28]. Realignment of the sequences, as well as the removal of the forward and reverse primer, was done by the bioinformatic tool Cutadapt v.3.4 [29]. To improve the evaluation of sequences, low-quality-3’ ends were trimmed based on a visual review of the quality plots [30]. By filtering the sequences based on the expected error, a minimum quality was guaranteed for sequence pairs. The remaining probable sequencing errors were identified using the error profile and rectified by dada2 [28]. Remaining sequence pairs were merged to amplicon sequence variants (ASVs) by dada2, when having an overlap of at least 25 base pairs without a mismatch (Callahan et al., 2016). The Protist Ribosomal Reference database (PR2) was used to taxonomically classify each ASV [31]. Potential chimeres, caused by the polymerase chain reaction, were discarded [28]. Further statistical analyses were performed using R v.4.0.3 (R Core Team, 2021) in R Studio v.1.2.5001 (R Studio Team, 2019). The rrarefy function from the vegan package v.2.5-7 was used for rarefaction and normalization [32]. Visualizations and further calculations were done with ampvis2 v.2.6.7, ggplot2 v.3.3.5, FactoMineR v.2.4, factoextra v.1.0.7.999, cluster v.2.1.4 and vegan [33, 34, 35, 36, 37]. Further editing (e.g.increasing the size of text and legends in graphs) were performed with Inkscape (version 1.1.1). The mapping of the stations as well as the nutrient distribution in the surface samples was done using QGIS (QGIS.org, 2022). The significance level was set at p < 0.05 for all calculations.

Environmental conditions

Based on the high-resolution physical characterization of the sub-mesoscale filament [11], we grouped the samples collected in an area inside the filament (IF), outside East and outside West of the filament (OE and OW). The area inside the filament includes a corridor of ~ 7 km along the 500 m wide streak of sea-ice (Fig. 1).

Environmental parameters, including temperature, salinity, and concentrations of nutrient- and Chl a differed in surface waters at a depth of ~ 10 m between samples collected inside and outside the filament (Fig. 2). Temperatures outside west of the filament (OW) were lower than inside the filament, while temperatures outside east of the filament (OE) were higher than inside the filament. Salinity values west of the filament were lower than inside the filament, while the values east of the filament were comparable to the values measured inside the filament. In general, nitrate concentrations were low with maximum levels of 0.17 µmol/l in the entire study area pointing towards depletion of nutrients during a past bloom period. The lowest nitrate levels occurred inside the filament. Photosynthetic biomass reflected by Chl a concentrations was slightly lower in the surface layer inside the filament (Fig. 2). Supplement S1 provides an overview of the environmental parameter measured for this study and the assignment of samples to the different spatial areas identified across the filament.

Vertical and horizontal plankton patterns in the sub-mesoscale frontal system

After quality filtering Illumina sequencing generated 5.058.797 high quality reads of the 18S rDNA V4-region from 74 samples collected near the surface on five different transects, and at six depths (10–400 m) at five CTD-stations across the sub-mesoscale filament. In a first step, community profiles were grouped based on Jaccard’s distances, into a non-metrical MDS ordination plot. There was a clear separation of 18S community profiles according to sampling depth and sampling location (Fig. 3). They grouped into four significantly distinct clusters (Fig. 3), supported by an ANOSIM (R = 0.81; p = 0.0001). Except for station 16, the 18S sequence composition in samples collected in the upper 30 m of the water column (cluster I and II) was significantly different from samples collected below (cluster III and IV), and significantly different between inside (cluster I) and outside (cluster II) of the filament. Cluster I contains samples collected mainly inside the filament, while cluster II contains samples outside the filament.

This clustering was the basis for determining the horizontal expansion on the biological response on physical forcing across the sub-mesoscale filament. Similar to the physical characterization of the study area, we identified a corridor of ~ 7 km width along the sea ice streak in which the plankton community composition was different from that of the adjacent community around ~ 2 km away (Fig. 4). The samples in cluster III originated from water-depths between 50–100 m, and samples in cluster IV were collected at most stations between 100–400 m.

Taxonomic composition of Eukaryotic plankton communities

Metazoa, Dinoflagellata and Ochrophyta constituted more than 50% of the eukaryotic microbial sequences in all samples of this study (Fig. 5A). Differences in the relative sequence abundance of these three taxonomic groups best explain segregation of the plankton community profiles in the MDSplot. The relative sequence abundance of Metazoa was highest in the upper 30 m of the water column inside the filament (p = 0.001), while sequence abundances of Dinoflagellata and Ochrophyta (mainly Bacillariophyta in this data set) were lower in the upper water column inside the filament compared to that outside the filament (p = 0.001). These data are supported by characteristic pigment data providing information on the contribution of different taxonomic groups to Chl a biomass. In cluster I diatoms and dinoflagellates were the main contributors to Chl a biomass, where diatoms contributed 40%-60% and dinoflagellates 25–35% to total Chl a biomass. In cluster II diatoms contributed similarly as in cluster I to Chl a, but here, chlorophytes were more abundant, making up between 20 (east of the front) to 50% (west of the front). The rest of the Chl a (~ 10% each) outside the filament was split up between dinoflagellates (mostly only east of the front) and prymnesiophytes (supplement S2).

Metazoan eDNA read abundances consisted mainly of the crustacean class Maxillopoda and the genera Calanus, Oithona, Metridia and Pseudocalanus. The use of LOKI allowed us for the first time to directly correlate sequence based information on zooplankton composition and distribution with optical counts. The metazoan read abundances were highly correlated to optical surveys of zooplankton abundances via the LOKI (Fig. 6). Moreover, the zooplankton composition was similar as derived from molecular and optical data. Both approaches thus confirm microscopic counts from multi net samples, which show higher zooplankton abundances in the upper 30 m inside the filament compared to outside the filament [4].

Among the Bacillariophyta, the relative sequence abundance of the sea ice related pelagic diatoms Actinocyclus curvatulus and Fragilariopsis cylindrus was higher in the upper water column inside the filament than outside the filament (p = 0.001). Inside the filament, at stations Stn 10 & Stn 12, the relative sequence abundances of A. curvatulus in the upper 10 m and in the Chl a max were in a similar range, while it was lower in the Chl a max than in the upper 10 m at Stn 14 outside the filament. In contrast, Chaetoceros socialis appeared with low relative sequence abundances in the upper water column inside the filament as compared to both, outside the filament and deep water layers at stations Stn 10, Stn 12, Stn 14 and Stn 18 (p = 0.017). This species is known to be a major contributor to ice-free Arctic phytoplankton communities [38]. It contributed more than 80% of all Bacillariophyta at a depth of 400 m at Stn 10 & 12, inside the filament, corresponding to around ~ 8% of all sequences identified at this depth (Fig. 5B).

Among the Dinophyceae, the relative sequence abundance of Gyrodinium fusiforme, Woloszynskia halophila and Pelagodinium beii was the highest in the upper 30 m of the water column (p = 0.030) (Fig. 5C). The sequences of Woloszynskia halophila were significantly more abundant in the 10 m samples from inside the filament (p = 0.007), while Gyrodinium heleveticum sequences were significantly more abundant in the 10 m samples from outside the filament (p = 0.042) (Fig. 5C). Gyrodinium fusiforme sequence abundances were not significantly different between samples from inside and outside the filament.

Ciliophora, Cercozoa, Haptophyta (mainly Phaeocystis sp.), Picozoa, Choanoflagellida, Chlorophyta contributed less than Dinoflagellata and Ochrophyta, but their relative abundances also showed a distinct pattern, being lower in samples from the upper water column inside the filament, compared to samples from outside the filament. The relative sequence abundances of Haptophyta, Picozoa and Chlorophyta were alleviated in some of the samples collected below 100 m, suggesting enhanced physical export of these taxa that are known to have low sinking rates [39, 40]. However, standard deviation in the deep community profiles was high (Fig. 5A).

There is growing evidence that marine physical features such as sub-mesoscale eddies or frontal systems impact ambient ecosystems, their biomes, and related biogeochemical processes by physical forcing that leads to changes in the availability of nutrients and/or distribution of organisms in a frontal system [3, 4]. In this study, we combined for the first time high-resolution measurements of the physical environment with horizontally high-resolution surveys of molecular plankton diversity and Chl a biomass composition and vertical profiles of physical and biological parameters for characterizing the effects of a marine physical feature on the surrounding ecosystem in the marginal ice zone of Fram Strait (Arctic Ocean). This unprecedented approach allowed us to estimate the spatial extent in which the planktonic community is influenced by a hydrographic sub-mesoscale structure. We observed significant horizontal and vertical differences in the 18S plankton community profiles and photosynthetic biomass that aligned with water mass distribution and physical forcing in the filament. Here, sea-ice melt related stratification might have been a key factor driving the community composition in the filament. Fram Strait is the main gate for sea ice export from the Arctic Ocean, and large parts of the western Fram Strait are ice-covered throughout the year. The eastern Fram Strait is mostly ice-free and the southern location of the ice edge in this area is highly variable. In consequence, environmental conditions in the eastern Fram Strait are severely impacted by the presence or absence of the marginal ice zone (MIZ) and its related processes, such as sea ice melt and meltwater stratification [41]. The latter is supposed to be of central importance for setting patterns of Arctic phytoplankton distribution [42]. This is supported by the findings of this study, that describe enhanced sequence abundances of sea-ice associated diatoms, zooplankton and enhanced POC- export within the sub-mesoscale feature that was strongly stratified by meltwater occupying the upper 25 m inside the filament.

In general, the ASVs identified in this data-set represented all major taxonomic eukaryotic microbial groups expected to be found in plankton samples in Fram Strait such as stramenopiles, dinoflagellates, Ciliophora, Haptophyta or Chlorophyta [43, 44]. The physical-chemical settings provided an environment for the development of plankton communities, which were on one hand horizontally distinct inside and outside the filament. On the hand, they were vertically distinct, with surface communities being different from those in the underlying water layers. Phytoplankton communities have a regularly recurring annual succession pattern [45, 46]. In Fram Strait, blooms of the diatom Chaetoceros spp. together with other centric diatoms are characteristic for pelagic phytoplankton blooms during early summer, followed by increasing abundances of small picoplankton such as Picozoa or Mammiellophyceae [46, 47]. The higher presence of 18S sequences associated with Chaetoceros sociales outside the filament in conjunction with small picoplankton taxa and elevated biomass of Chlorophyta suggested the presence of a late bloom-community of pelagic origin at 10 m depth during the observation period in vicinity of the filament. Moreover, the presence of significant shares of Chaetoceros socialis sequences at a depth of 400 m only at stations inside the filament furthermore suggests enhanced export in this area of the filament via a combination of physical forcing amplifying the biological carbon pump. This assumption is in line with enhanced carbon flux inside the filament [3], which was hypothesized to be a consequence of a subduction process ongoing in the area [11]. In contrast to the distribution pattern of Chaetoceros socialis, sequences of Actinocyclus curvatulus, had higher relative sequence abundance inside the filament compared to outside the filament. The occurrence of this species is tightly linked to the presence of sea ice [48], suggesting it to be an indicator for a sea ice related community having elevated growth in the meltwater layer inside the filament, which went down to 25 m. Thus, Chl a biomass inside the filament might be associated with growth of sea-ice associated diatoms near the streak of sea-ice in the center of the filament.

Among the dinoflagellates, we found besides marine also fresh-water taxa, such as Gyrodinium helveticum [49, 50]. One reason for the occurrence of freshwater species in the surface waters could be the melting of sea ice, which transfers fresh-water taxa from the melt ponds to the upper water column. This would also explain the higher relative abundance of Woloszynskia halophila inside the filament, This species is a brackish Dinophyceae species usually found in the Baltic Sea [51], but also in sea-ice associated systems [52]. Although we find these fresh-water and brackish taxa, it is unclear if they stay alive under marine conditions or if we only detect the remains of the algae. The heterotrophic marine species Gyrodinium spirale found with higher sequence abundance inside the filament is often associated with diatom blooms in pelagic systems [53]. At high latitudes, sea-ice melt might serve as an iron source to promote phytoplankton blooms in its vicinity [54], and it is known that heterotrophic Gyrodinium species, like G. spirale and G. fusiforme, prefer iron-enriched blooms [55]. Thus, we postulate that alleviated sequence abundance of G. spirale as well as higher sequence abundances of sea-ice related diatom species inside the filament might be related to higher sea-ice melt related iron concentrations in a small geographical area inside the filament. Moreover, iron-enriched blooms are preferential food sources for zooplankton [56], which had alleviated sequence and species abundances inside the filament. In summary, the community composition inside the front system might indirectly suggest iron-intake from sea-ice near melting sea ice, though this study is missing information on iron concentrations in the study area. Future studies in the marginal ice-zone should assess the iron-flux from sea-ice to improve our mechanistic understanding of linkages between sea-ice melt and plankton productivity in the Arctic Ocean.

In this study the eDNA data for zooplankton correlated surprisingly well with the optical in situ surveys via LOKI. Here, it was certainly an advantage that LOKI captured the fine scale distribution of zooplankton which allowed to relate the zooplankton density within an interval from 5 m above to 5 m below each water sampling depth. Net samples integrate over much larger, pre-fixed intervals, and thus correlations between vertical point measurements in the water column and integrated net samples likely would have resulted in a much larger variation. The zooplankton taxa found in this study via eDNA analyses, such as Calanus., Oithona., Maxillopoda, Metridia sp and Pseudocanalus. are all well known as contributors to plankton communities in Fram Strait [4], and corresponded well to the image data. Since we are aware that the 18S rRNA gene is of limited value for a detailed and reliable taxonomic identification of zooplankton, we did not carry out correlation analyses at species level. Our results nevertheless suggest that using the 18S rRNA gene has the potential to serve as a proxy for zooplankton abundances, and should thus be included in a smart holistic long-term observation strategy. Further evaluations of eDNA based zooplankton surveys with optical counts, however, are needed to test this conclusion.

In summary, our data suggest that sea ice melt and meltwater stratification is promoting elevated abundance and growth of ice-related taxa inside the filament. The presence and higher relative abundances of taxa, which are known to be tightly linked to high sea-ice concentrations, coincided with higher relative sequence abundance of metazoa. The distribution pattern of 18S eDNA sequences associated with metazoa are in line with the findings of the assessment of zooplankton composition via optical in-situ observations described in this study, but also net-tows that found higher abundances, of zooplankton inside the filament [4]. The concentration of zooplankton inside the filament in conjunction with the enhanced growth of sea-ice associated phytoplankton is likely the result of the convergence of surface water associated with the sub-mesoscale filament, which also led to the accumulation of sea ice [11]. Despite differences in sequence abundances of phyto- and zooplankton inside and outside the filament, Chl a concentrations at 10 m were only slightly different. Assuming zooplankton grazing, nearly similar Chl a concentrations inside and outside the filament might point towards higher primary productivity inside the filament, as suggested previously based on differences in nutrient concentrations [3]. The reproductive activity of the zooplankton inside the filament might be a consequence of the availability of sea-ice associated phytoplankton, which might have a better food-quality for zooplankton compared to phytoplankton dominated by Chaetoceros sp. [57]. Thus, alleviated carbon export inside the filament might be a combination of enhanced carbon export due to zooplankton grazing on sea-ice associated phytoplankton and enhanced export of biomass due to subduction processes inside the filament.

In conclusion, only by combining high-resolution eDNA analysis near the water surface, targeted vertical profiling at CTD stations and optical observations via the LOKI, we were able to assess the spatial extent of ecosystem change at both horizontal and vertical scales near the sub-mesoscale structure. We were able to show that changes in biology within the filament were closely linked to physical structuring and occupied a similar space. The high-resolution data from the upper water column showed a clearer picture of plankton community boundaries near the filament than the five CTD stations and the vertical optical surveys [3, 4]. The latter, in turn, provided important information on the vertical distribution of plankton communities and carbon flux in the study area. By combining the information from the horizontal and vertical surveys, we assume that sea-ice melt and associated physical and biological processes in a sub-mesoscale filament may have led to elevated POC-export in an overall area of ~ 350 km² of Fram Strait during summer of 2017, illustrating spatially far reaching consequences of sea-ice melt for marine ecosystems. This is consistent with previous observations from long-term sediment trap observations demonstrating higher export of POC in regions with seasonal sea ice cover or near the ice edge compared to ice-free regions, although further high-resolution studies in the marginal ice zone are needed to make more precise estimates of the extent of spatial impacts of near-ice-edge processes on marine ecosystems.

Acknowledgements

We thank the captain and crew of RV Polarstern for their excellent support during expedition PS107. This work was accomplished as part of the grants AWI_PS107_4 and AWI_PS107_9. We are very grateful to the chief scientist Ingo Schewe and all other cruise participants with whom we developed this interdisciplinary “Front-Study'' in a communal effort. Moreover, we are very grateful to Kerstin Korte and Swantje Ziemann for excellent technical support. This study was conducted in the framework of the HGF Infrastructure Program FRAM of the Alfred-Wegener-Institute Helmholtz Centre for Polar and Marine Research. NH was supported by the Russian–German Research Cooperation QUARCCS funded by the German Ministry for Education and Research (BMBF) under grant 03F0777A. AB was partly funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Projektnummer 268020496 – TRR 172, within the Transregional Collaborative Research Center “ArctiC Amplification: Climate Relevant Atmospheric and SurfaCe Processes, and Feedback Mechanisms (AC)3 in subproject C03. Ship time was provided under grant AWI_PS107_07.

Author contributions

KM, WJvA and NH conducted the fieldwork. JW conducted the eDNA data analyses. JW and KM drafted the manuscript. MG carried out the nutrient analyses. BN and NH conducted the analyses of the LOKI-data. SN conducted the quality check and annotation of 18S sequencing data. AB conducted calculations of Chl a for phytoplankton groups from HPLC and ACS data and supervised the related fieldwork. All authors contributed to improving the final manuscript, by contributions to the scientific interpretation of the data and the discussion of results.

Data availability statement (mandatory)

Raw sequencing data are deposited in ENA under accession number PRJEB66268. All other data are available via PANGAEA or provided as part of this manuscript.

Lévy, M., Klein, P. & Treguier, A. M. Impact of sub-mesoscale physics on production and subduction of phytoplankton in an oligotrophic regime. Journal of marine research 59, 535–565 (2001).
Acha, E. M., Piola, A., Iribarne, O. & Mianzan, H. Ecological processes at marine fronts: oases in the ocean. Springer (2015).
Fadeev, E., et al. Sub-mesoscale physicochemical dynamics directly shape bacterioplankton community structure in space and time. Limnology and Oceanography 66, 2901–2913 (2021).
Kaiser, P., et al. Effects of a Sub-mesoscale Oceanographic Filament on Zooplankton Dynamics in the Arctic Marginal Ice Zone. Frontiers in Marine Science 8, 853 (2021).
Bartley, T. J., et al. Food web rewiring in a changing world. Nature ecology & evolution 3, 345–354 (2019).
Pecuchet, L., et al. Spatio-temporal dynamics of multi‐trophic communities reveal ecosystem‐wide functional reorganization. Ecography 43, 197–208 (2020).
Klein, P., et al. Upper ocean turbulence from high-resolution 3D simulations. Journal of Physical Oceanography 38, 1748–1763 (2008).
Klein, P. & Lapeyre, G. The oceanic vertical pump induced by mesoscale and sub-mesoscale turbulence. Annu. Rev. Mar. Sci 1, 351–375 (2009).
Bracher, A., et al. High-resolution water column phytoplankton composition across the Atlantic Ocean from ship-towed vertical undulating radiometry. Frontiers in Marine Science 7, 235 (2020).
Lévy, M., Franks, P. J. & Smith, K. S. The role of sub-mesoscale currents in structuring marine ecosystems. Nature communications 9, 1–16 (2018).
von Appen, W. J., et al. Observations of a sub-mesoscale cyclonic filament in the marginal ice zone. Geophysical Research Letters 45, 6141–6149 (2018).
Liu, X., & Levine, N. M. Enhancement of phytoplankton chlorophyll by sub-mesoscale frontal dynamics in the North Pacific Subtropical Gyre. Geophysical research letters 43, 1651–1659 (2016).
Guo, M., Xiu, P., Chai, F. & Xue, H. Mesoscale and sub-mesoscale contributions to high sea surface chlorophyll in subtropical gyres. Geophysical Research Letters 46, 13217–13226 (2019).
Metfies, K., et al. High-resolution monitoring of marine protists based on an observation strategy integrating automated on-board filtration and molecular analyses. Ocean Sci. 12, 1237–1247 (2016).
Metfies, K., et al. Uncovering the intricacies of microbial community dynamics at Helgoland Roads at the end of a spring bloom using automated sampling and 18S meta-barcoding. Plos One 15, e0233921 (2020).
Kawasaki, T. & Hasumi, H.. The inflow of Atlantic water at the Fram Strait and its interannual variability. Journal of Geophysical Research: Oceans 121, 502–519 (2016).
Liu, Y., et al. Retrieval of phytoplankton pigments from underway spectrophotometry in the Fram Strait. Remote Sensing 11, 318 (2019).
Schulz, J., et al. Imaging of plankton specimens with the lightframe on-sight keyspecies investigation (LOKI) system. Journal of the European Optical Society - Rapid publications 5, 10017s (2010).
Murphy, J. & Riley, J. P. A modified single solution method for the determination of phosphate in natural waters. Analytica chimica acta 27, 31–36 (1962).
Strickland, J. D. H. & Parsons, T. R. A practical handbook of seawater analysis. Research Board of Canada, Bulletin. 167, 65 (1972).
Grasshoff, K., Kremling, K. & Ehrhardt, M. (Eds.). Methods of seawater analysis (John Wiley & Sons, 2009).
Kérouel, R. & Aminot, A. Fluorometric determination of ammonia in sea and estuarine waters by direct segmented flow analysis. Marine Chemistry 57, 265–275 (1997).
Liu, Y., et al. Phytoplankton pigment concentration measured by HPLC during POLARSTERN cruise PS107. Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Bremerhaven, PANGAEA, https://doi.org/10.1594/PANGAEA.894860 (2018).
Liu, Y., et al. Phytoplankton pigment concentration estimated from underway AC-S particulate absorption data during POLARSTERN cruise PS107. PANGAEA, https://doi.org/10.1594/PANGAEA.898100 (2019).
Vidussi, F., Claustre, H., Manca, B.B., Luchetta, A. & Marty, J.C., Phytoplankton pigment distribution in relation to upper thermocline circulation in the eastern Mediterranean Sea during winter. J. Geophys. Res. 106, 19939–19956 (2001).
Uitz, J., Claustre, H., Morel, A., Hooker S. B. Vertical distribution of phytoplankton communities in open ocean: An assessment based on surface chlorophyll. Journal of Geophysical Research. Oceans 111, DOI: 10.1029/2005JC0032072006 (2006).
Diepenbroek, M., et al. Informatik 2014 – Big Data Komplexität meistern. GI-Edition: Lecture Notes in Informatics (LNI) – Proceedings. (ed. GI) 232, 1711–1724. (Köllen Verlag; 2014).
Callahan, B. J., et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nature methods 13, 581–583 (2016).
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. Journal 17, 10–12 (2011).
Pauli, N. C., et al. (2021). Selective feeding in Southern Ocean key grazers—diet composition of krill and salps. Communications biology 4, 1–12 (2021).
Guillou, L., et al. The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote small sub-unit rRNA sequences with curated taxonomy. Nucleic acids research 41, D597-D604 (2012).
Oksanen, J., et al. Package ‘vegan’. Community ecology package, version, 2(9), 1–295 (2013).
Andersen, K. S., Kirkegaard, R. H., Karst, S. M., & Albertsen, M. ampvis2: an R package to analyse and visualise 16S rRNA amplicon data. BioRxiv, 299537 (2018).
Wickham, H. ggplot2. Wiley interdisciplinary reviews: computational statistics 3, 180–185 (2011).
Lê, S., Josse, J. & Husson, F. FactoMineR: an R package for multivariate analysis. Journal of statistical software 25, 1–18 (2008).
Kassambara, A. & Mundt, F. Factoextra: Extract and Visualize the Results of Multivariate Data Analyses. R Package Version 1.0.7. https://CRAN.R-project.org/package=factoextra (2020).
Maechler, M., Rousseeuw, P., Struyt, A. & Hubert, M. cluster: Cluster Analysis Basics and Extensions. R Package Version 2.0.1. (2015)
Gradinger, R.R. & Baumann, M.E.M. Distribution of phytoplankton communities in relation to the large-scale hydrographical regime in the Fram Strait. Marine Biology 111, 311–321 (1991).
Wolf, C., Iversen, M., Klaas, C. & Metfies, K. Limited sinking of Phaeocystis during a 12 days sediment trap study. Molecular Ecology 25, 3428–3435 (2016).
Bachy, C., et al. (2022). Phytoplankton Surveys in the Arctic Fram Strait Demonstrate the Tiny Eukaryotic Alga Micromonas and Other Picoprasinophytes Contribute to Deep Sea Export. Microorganisms 10, 961 (2022).
Nöthig, E. M., et al. Summertime chlorophyll a and particulate organic carbon standing stocks in surface waters of the Fram Strait and the Arctic Ocean (1991–2015). Frontiers in Marine Science 7, 350 (2020).
Lester, C. W., Wagner, T. J. W., McNamara, D. E. & Cape, M. R. The Influence of Meltwater on Phytoplankton Blooms Near the Sea-Ice Edge. Geophysical Research Letters 48, e2020GL091758 (2021).
Kilias, E., Wolf, C., Nöthig, E-M., Peeken, I. & Metfies, K. Protist distribution in the Western Fram Strait in Summer 2010 based on 454-Pyrosequencing of 18S rDNA. J. Phycol. 49, 996–1010 (2013).
Metfies, K., von Appen W-J., Kilias, E., Nicolaus, A. & Nöthig, E-M. (2016) Biogeography and photosynthetic biomass of Arctic marine pico-eukaroytes during summer of the record sea ice minimum 2012. PLOS ONE 11, e0148512 (2016).
Kaese, L., et al. Metabarcoding analysis suggests that flexible food web interactions in the eukaryotic plankton community are more common than specific predator-prey relationships at Helgoland Roads. ICES Journal of Marine Science 78, 3372–3383 (2021).
Wietz, M., et al. (2021). The polar night shift: seasonal dynamics and drivers of Arctic Ocean microbiomes revealed by autonomous sampling. ISME Communications 1, 76 (2021).
Hodal, H. & Kristiansen, S. (2008). The importance of small-celled phytoplankton in spring blooms at the marginal ice zone in the northern Barents Sea. Deep Sea Research Part II: Topical Studies in Oceanography, 55(20–21), 2176–2185.
Oksman, M., Juggins, S., Miettinen, A., Witkowski, A. & Weckström, K. The biogeography and ecology of common diatom species in the northern North Atlantic, and their implications for paleoceanographic reconstructions. Marine Micropaleontology 148, 1–28 (2019).
Takano, Y. & Horiguchi, T. Surface ultrastructure and molecular phylogenetics of four unarmored heterotrophic dinoflagellates, including the type species of the genus Gyrodinium (Dinophyceae). Phycological Research 52, 107–116 (2004).
Logares, R., Shalchian-Tabrizi, K., Boltovskoy, A. & Rengefors, K. Extensive dinoflagellate phylogenies indicate infrequent marine–freshwater transitions. Molecular Phylogenetics and Evolution 45, 887–903 (2007).
Jaanus, A., et al. Distribution patterns of isomorphic cold-water dinoflagellates (Scrippsiella/Woloszynskia complex) causing ‘red tides’ in the Baltic Sea. Hydrobiologia 554, 137–146 (2006).
Bachy, C., López-García, P., Vereshchaka, A. & Moreira, D. Diversity and vertical distribution of microbial eukaryotes in the snow, sea ice and seawater near the North Pole at the end of the polar night. Frontiers in microbiology 2, 106 (2011).
Sherr, E. B. & Sherr, B. F. Heterotrophic dinoflagellates: a significant component of microzooplankton biomass and major grazers of diatoms in the sea. Marine Ecology Progress Series 352, 187–197 (2007).
Wang, S., Bailey, D., Lindsay, K., Moore, J. K. & Holland, M. Impact of sea ice on the marine iron cycle and phytoplankton productivity. Biogeosciences 11, 4713–4731 (2014).
Saito, H., Ota, T., Suzuki, K., Nishioka, J. & Tsuda, A. (2006). Role of heterotrophic dinoflagellate Gyrodinium sp. in the fate of an iron induced diatom bloom. Geophysical Research Letters 33, https://doi.org/10.1029/2005GL025366 (2006).
Olsen, Y., et al. A comparative study of responses in plankton food web structure and function in contrasting European coastal waters exposed to experimental nutrient addition. Limnology and oceanography 51, 488–503 (2006).
Hargave, B.T. & Geen G.H. (1970). Effect of copepod grazing on two natural phytoplankton populations. Journal of the Fisheries Board of Canada 27, (1970).

No competing interests reported.

Download PDF

Journal Publication

published 08 Apr, 2024

Read the published version in Scientific Reports →

Editorial decision: Revision requested
13 Nov, 2023
Reviews received at journal
25 Oct, 2023
Reviewers agreed at journal
15 Oct, 2023
Reviewers invited by journal
14 Oct, 2023
Editor assigned by journal
12 Oct, 2023
Editor invited by journal
10 Oct, 2023
Submission checks completed at journal
10 Oct, 2023
First submitted to journal
21 Sep, 2023

You are reading this latest preprint version

Unprecedented insights into extents of biological responses to physical forcing in an Arctic sub-mesoscale filament by combining high-resolution measurement approaches

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Material and Methods

Results

Discussion

Declarations

References

Additional Declarations

Supplementary Files

Status:

Journal Publication

Version 1