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

Research article 23 Mar 2020

Research article | 23 Mar 2020

# Population dynamics of modern planktonic foraminifera in the western Barents Sea

Population dynamics of modern planktonic foraminifera in the western Barents Sea
Julie Meilland1, Hélène Howa2, Vivien Hulot3,4, Isaline Demangel5,6, Joëlle Salaün7, and Thierry Garlan7 Julie Meilland et al.
• 1MARUM – Center for Marine Environmental Sciences, Leobener Str. 8, 28359 Bremen, Germany
• 2LPG-BIAF, UMR-CNRS 6112, University of Angers, Angers, France
• 3University of French Polynesia, UMR-241 EIO, Labex Corail, FAA'A, Tahiti, French Polynesia
• 4Département Ressources Biologiques et Environnement, Ifremer, UMR-241 EIO, Labex Corail, Vairao, Tahiti, French Polynesia
• 5Institute of Earth Sciences, University of Graz, NAWI Graz Geocenter, Graz, Austria
• 6Department of Geology, University of Lund, Lund, Sweden
• 7SHOM – Sciences et Techniques Marines/Géologie Marine, Brest, France

Correspondence: Julie Meilland (jmeilland@marum.de)

Abstract

This study reports on diversity and distribution of planktonic foraminifera (PF) in the Barents Sea Opening (BSO). Populations of PF living in late summer (collected by means of stratified plankton tows) and recently deposited individuals (sampled by interface corer) were compared. High abundances reaching up to 400 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$ in tow samples and 8000 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$ in surface sediments were recorded in the centre of the studied area while low abundances were observed in coastal areas, likely due to continental influences. The living and core-top assemblages are mainly composed of the same four species Neogloboquadrina pachyderma, Neogloboquadrina incompta, Turborotalita quinqueloba and Globigerinita uvula. The two species G. uvula and T. quinqueloba dominate the upper water column, whereas surface sediment assemblages display particularly high concentrations of N. pachyderma. The unusual dominance of G. uvula in the water sample assemblages compared to its low proportion in surface sediments might be the signature of (1) a seasonal signal due to summer phytoplankton composition changes at the BSO, linked to the increase in summer temperature at the study site, and/or (2) a signal of a larger timescale and wider geographical reach phenomenon reflecting poleward temperate/subpolar species migration and consecutive foraminiferal assemblage diversification at high latitudes due to global change. Protein concentrations were measured on single specimens and used as a proxy of individual carbon biomass. Specimens of all species show the same trend, a northward decrease in their size-normalized-protein concentration. This suggests that foraminiferal biomass is potentially controlled by different constituents of their organelles (e.g. lipids). The coupling of data from plankton tows, protein measurements and surface sediments allows us to hypothesize that PF dynamics (seasonality and distribution) are decoupled from their metabolism.

1 Introduction

Polar areas are sensitive to global temperature changes, particularly in the Arctic where warming occurs faster than in the rest of the world and has accelerated over the past 50 years (Shepherd, 2016). This Arctic amplification appears to be mainly caused by sea-ice loss under increasing CO2 (Dai et al., 2019). Recently increased advection of Atlantic Water into the Barents Sea modifies its physico-chemical properties (Smedsrud et al., 2013), which directly affect the entire ecosystem in the region. Higher temperatures lead to increased rates of planktonic primary production (Vaquer-Sunyer et al., 2013), and increased CO2 concentrations are expected to have a fertilization effect on marine autotrophs (Holding et al., 2015). Not only is productivity increasing, but spring and summer blooms are also occurring earlier in the European Arctic Ocean (Oziel et al., 2017). As a response, some taxa of calcifying groups (i.e. foraminifera, coccolithophores, molluscs and echinoderms; Beaugrand et al., 2013) exhibit a poleward movement in agreement with expected biogeographical changes under sea temperature warming. Both satellite images (Smyth et al., 2004; Burenkov et al., 2011) and in situ measurements (Dylmer et al., 2013; Giraudeau et al., 2016; AMAP, 2018) have shown rapid expansion of temperate species of coccolithophores in the Arctic. For example, Emiliania huxleyi shows a striking poleward shift ($>\mathrm{5}{}^{\circ }$) in the distribution of its blooms (Neukermans et al., 2018). This phenomenon, called “atlantification” (Årthun et al., 2012), is expected to impact every trophic level of the food web, from small phytoplankton species (Neukermans et al., 2018) to larger organisms (Dalpadado et al., 2012). Recent studies have investigated the ecology and biodiversity of planktonic foraminifera from the high-latitude North Atlantic (e.g. Schiebel et al., 2017). The species N. pachyderma comprises more than 90 % of recent assemblages (i.e. found in surface sediments) from the Polar region, north of Iceland (Kucera et al., 2005). Rather few studies on living planktonic foraminifera (PF) communities have focussed on (sub-)Arctic regions. Plankton tows also show that N. pachyderma is the dominant species in the Arctic (∼90 %) followed by T. quinqueloba (5 %) (Volkmann, 2000; Pados and Spielhagen, 2014). Through the compilation of population density profiles from 104 stratified plankton tow hauls collected in the Arctic and the North Atlantic oceans, Greco et al. (2019) investigated the ecology of N. pachyderma, in particular the variability of its habitat depth, and underlined the knowledge gap on its ecological preferences. In the western subpolar North Atlantic (Irminger Sea), the production of N. pachyderma shows two peaks, in spring and late summer, while winter shows a low production (Jonkers et al., 2010, 2013). The diversity of PF has increased in polar waters over the past decades, even though it remains low in comparison to lower latitudes, due to the poleward migration of warm-water species (Schiebel et al., 2017). A similar process appears to occur in the Southern Hemisphere (Meilland et al., 2016). Some species from lower latitudes are described as new components of poleward assemblages. The shift of PF assemblages to warmer conditions, since the pre-industrial stage, has been very recently highlighted more globally in the Northern Hemisphere (Jonkers et al., 2019). These major modifications in PF distribution patterns display changes more related to primary production than to water temperature itself (e.g. Jonkers et al., 2010; Schiebel et al., 2017). Planktonic foraminifera, being sensitive to ambient water geochemistry, are considered good indicators of the polar changing environments (Schiebel et al., 2017). More studies on living PF communities in the Arctic regions are needed to assess the spatial and temporal variability in their population dynamics and to better constrain today's polar and subpolar species ecological preferences.

Taking the opportunity of a cruise dedicated to the exploration of the physical oceanography of the western Barents Sea (MOCOSED 2014 cruise), we investigated the connections between the spatial variability of living PF, phytoplankton communities (Giraudeau et al., 2016) and the hydrological system through a south-to-north transect, between northern Norway and Spitsbergen [68–76] N. Along this transect, we compared PF living faunas (from plankton tow) to the assemblages found on the sea floor (from core-top sediments) in order to investigate eventual recent changes in their population dynamics. This latitudinal transect also gave us the opportunity to quantify protein concentrations of individual living PF in this area for the first time and along a physico-chemical gradient to see if and how it varies and explore how planktonic foraminifera from the same species may adjust to different environments.

2 Oceanographic setting

The studied area covers the western Barents Sea margin, i.e. Barents Sea Opening (BSO), where the surface and intermediate ocean circulation are characterized by the confrontation of the North Atlantic and the Arctic waters (Fig. 1). The seasonal and interannual dynamics of these two water masses, interacting with the complex topography of the western margin (Storfjordrenna and Bjørnøyrenna glacial troughs separated by shallow Spitsbergenbanken), determine the position of the Polar Front (Loeng, 1991). The Norwegian Atlantic Current (NwAC) carries Atlantic Water into the Barents Sea. Along the western Barents Sea margin, Atlantic Water is transported to the Fram Strait by the West Spitsbergen current (Skagseth et al., 2008; Oziel et al., 2017; Fig. 1).

Figure 1Sampling map of the MOCOSED cruise in the western Barents Sea with schematic surface circulation (red arrows: Atlantic Water; blue arrows: Arctic Water; green arrow: Norwegian Coastal Current; adapted from Oziel et al., 2017). Little green dots display the 32 CTD–Niskin stations along the south-to-north transect; large blue circles show the location of the nine MultiNet stations; medium red circles underline the five multicore stations.

In the southern part of the transect a strong thermohalocline clearly underlined a surface mixed layer of about 30–35 m depth. This cline slightly deepened northwards and blurred out north of 74.5 N, where no stratification was observed in the water column. From south to north of the transect (i) high temperatures and low salinities reflected the Norwegian Coastal Water (NwCW), which is also enriched in Chl a. The relatively warm NwCW (8.5 to 11 C) extended northwards up to 74.5 N, overlying the colder Norwegian Atlantic Water (NwAW). Less saline (33.5) to the south, NwCW became saltier (34.9) in the vicinity of the Spitsbergenbanken; (ii) at the northern end of the transect, the NwAW penetrated the Barents Sea through the Storfjordrenna trough with temperatures from 6 to 8 C and an open marine salinity of 35.1 (Giraudeau et al., 2016).

The Chl a content followed the hydrological pattern described above (Fig. 2c). Relatively high concentrations (mean $\approx \mathrm{0.8}\phantom{\rule{0.125em}{0ex}}\mathrm{mg}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$) were located in the surface mixed layer composed of NwCW. The highest values around 1.25 mg m−3 were recorded off the Norwegian coast. Chl a content decreased northwards (north of 74.5 N) and reached ≈0.4mg m−3 in the upper layer (0–60 m) of the well-mixed NwAW. The composition of the phytoplankton community observed in surface water at seven stations along the studied transect was essentially dominated by three algal groups: fuco-flagellates (25 % to 43 %; major component Phaeocystis pouchetii), Prasinophytes (15 % to 30 %; major components Micromonas pusilla and Bathycoccus pusilla) and Prymnesiophytes (13 % to 24 %; major component Emiliania huxleyi) (Giraudeau at al., 2016). Three other features are noteworthy (Fig. 2d): (i) the dominance of dinoflagellates (24 %) at the southernmost station of the transect (close to the Norwegian coast) contrasted with their total absence in the well-mixed NwAW, north of 74.5 N; (ii) the presence of diatoms (10 % to 20 %) in the surficial NwCW, but rare (<5 %) to the north; (iii) the constant increase in relative abundance of Cyanobacteria from <5 % to more than 15 % along the south-to-north transect.

Figure 2South–north, 0–100 m deep sections across the Barents Sea Opening, compiled from data of the 32 CTD casts: (a) temperature (C), (b) salinity, and (c) Chl a concentrations (mg m−3). Vertical dashed lines and associated numbers correspond to the location of the MultiNet hauls. The lower right panel (d) displays pigment concentrations translated into the relative abundances of the major phytoplankton species at seven sampling stations along the transect (Giraudeau et al., 2016).

3 Material and methods

In late summer 2014, from 8 August to 20 September, the SHOM (French Hydrographic Office) operated the oceanographic cruise MOCOSED 2014, on board the R/V Pourquoi pas?. Along a 700 km south-to-north transect from the Norwegian (68 N) to the Spitsbergen (76 N) coasts, investigations of hydrological processes at the BSO were carried out coupled with the exploration of the phytoplankton and foraminiferal communities using a total of 32 vertical casts deployed ≈20 km apart from each other (Figs. 1 and 2).

## 3.1 Living planktonic foraminifera from stratified plankton samples (MultiNet)

Living PF were collected at 7 of the 32 conductivity–temperature–depth (CTD) south-to-north transect stations (3 to 9) and at two stations (1 and 2) located west to east ≈20km apart from the central point of the main south-to-north CTD transect (Fig. 1; Table 1), using a stratified plankton tow (MultiNet Hydro-Bios type Midi, opening of 0.25 m2) equipped with five nets (mesh size 100 µm to avoid nets clogging in case of intense phytoplankton bloom). This sampling strategy was used in order to observe the potential effect of latitudinal changes but also of bathymetry, longitudinally, on PF distribution. At each station, one vertical haul sampled five successive water layers from the sea surface to 100 m depth. A second haul has been deployed for the stations 1, 3, 4, 5, 6, 7 and 8, to collect material down to 700 m depth (Table 1). For each of the depth intervals (0–20, 20–40, 40–60, 60–80, 80–100; 0–100, 100–200, 200–300, 300–500 and 500–700 m), the filtered water volume was measured by means of a flowmeter attached to the MultiNet mouth. Each MultiNet sample was preserved in a 250 mL vial with ethanol (90 %) buffered with hexamethylenetetramine until processing at the land-based laboratory. Back at the laboratory, MultiNet samples were washed over a 100 µm mesh, and all foraminifera were removed from the sample and dried in an oven at 50 C. All living PF, distinguished by their coloured cytoplasm visible through the shell, were picked, stored in counting cells and identified at the species level, following the SCOR WG138 taxonomy as implemented in Siccha and Kucera (2017). Empty tests, considered dead individuals, were counted separately. Correlations following a non-metric multidimensional scaling (NMDS) ordination were carried out with the R package Vegan (Oksanen et al., 2013). Using the Bray–Curtis distance these correlations were tested between PF species absolute concentrations, the latitude of the station and parameters of the ambient waters (temperature, salinity, Chl a concentration). Results are given in relative abundances (percent of the total, live or dead fauna) or in absolute abundances in number of individuals per cubic metre of filtered water ($\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$).

Table 1Location (latitude and longitude), sampling date and water depth, of the nine MultiNet and five multitube stations, incremented from south to north (stations 1 and 2 being positioned beside the main transect midpoint 6). Stations where phytoplankton analyses were performed are also indicated.

### Protein biomass and test size measurements

For protein extraction and measurement, a few living individuals (≈60) were picked on board from the shallowest water samples at stations 3 to 9 immediately after sampling. Only shells that were completely filled with cytoplasm were selected. After picking, individuals were immediately cleaned with a brush and filtered seawater to remove all particles, including organic matter, that were stuck to the test. The individuals were stored in a 1.5 mL Eppendorf vial and analysed on board, using the bicinchoninic acid (BCA) method as explained in Meilland et al., (2016). Morphometric analyses on single foraminiferal tests were carried out at the University of Angers with an automated incident light microscope (Bollmann et al., 2004; Clayton et al., 2009) at a resolution of 1.4 µm2. Images were analysed for their two-dimensional (silhouette) morphometry (Beer et al., 2010), including minimum test diameter, which is the shortest distance from wall to wall that passes through the centre of the proloculus (the initial chamber of a foraminifer). Protein-to-size relations were determined for the minimum diameter of each test providing size-normalized protein content (SNP) for data analyses. Foraminifera protein concentrations were linearly normalized to 1 µm minimum test diameter, being aware of any unavoidable errors related to non-linear increments of biomass at volumetric test growth (see Beer et al., 2010).

## 3.2 Fossil planktonic foraminifera assemblages from core tops (multitube)

At five sites of the main CTD transect, an interface corer (multitube type, Oktopus GmbH, Institut National des Sciences de l'Univers division of Brest, France) was used to simultaneously obtain eight short sediment cores (less than 1 m in length) (Fig. 1; Table 1). At each station, the core with the more even and undisturbed water–sediment interface was selected. The uppermost 0.5 cm of the core was sampled and fixed with 95 % ethanol. Samples were stained with Rose Bengal, reacting with cytoplasm, to distinguish PF still bearing cytoplasm (fresh or in degradation) and thus very recently deposited from fossil PF without cytoplasm. Stained shells of foraminifera probably reflect the spring and summer population of the year of sampling, even though the exact degradation time of cytoplasm is still poorly constrained (Schönfeld et al., 2013). The core-top sediments were wet sieved using a 100 µm mesh and the foraminifera were identified using the same SCOR WG138 taxonomy in order to be directly comparable to the plankton tow samples.

4 Results

## 4.1 Planktonic foraminifera diversity and distribution in the water column

Data from the seven stations of the south-to-north CTD transect with five depth values per station were compiled to display the repartition of PF absolute abundances (for total assemblage and species specific) in the upper 100 m of the vertical section across the BSO. Abundances of living specimens below 100 m depth were very low (total concentration $<\mathrm{5}\phantom{\rule{0.125em}{0ex}}\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$) and are therefore not discussed. Data from station 1 (western) and station 2 (eastern) located 20 km on either side of the S–N transect were compared to the data of station 6 at the middle of the transect.

Figure 3Distribution of planktonic foraminifera total abundances ($\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$) in the 0–100 m depth section across the Barents Sea Opening (a) and in the west–east transect (stations 1, 6 and 2; b). Station names are indicated above vertical alignments of five dots representing the middle points of the five net sampling intervals.

Total concentrations of living PF fauna varied between 0 and 400 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$ (Fig. 3). Along the south-to-north CTD transect, the highest concentrations were all observed above 20 m water depth, in the surface mixed layer of the well-stratified water area, i.e. at the edge of the NwCW. The two stations located at the south and north extremities of the transect (3, off the Norwegian coast, and 9, off the Spitsbergen coast) displayed low densities (10 to 50 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$). At station 1, located above the Barents Sea margin slope, the maximum abundance of 220 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$ was recorded in a deeper habitat (20–40 m, Fig. 3). Station 2, located inside the Barents Sea, was very poor in PF ($<\mathrm{10}\phantom{\rule{0.125em}{0ex}}\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$). In total, 10 species were observed. The studied area was characterized by high abundances of subpolar to polar species (Fig. 4), listed in descending order: Globigerinita uvula (45 %), Turborotalita quinqueloba (26 %), Neogloboquadrina incompta (15 %) and Neogloboquadrina pachyderma (9 %). There was a notable presence of the temperate water species Globigerina bulloides (3 %) and negligible percentages (<1 %) of Globigerinita glutinata, Neogloboquadrina dutertrei, Globigerinoides ruber, Globigerinoides sacculifer and Orcadia riedeli.

Figure 4Distribution of the four major species abundances ($\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$) in the 0–100 m depth section across the Barents Sea Opening. (a) For Globigerinita uvula with species abundances (z axes) from 0 to 250 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$; for (b) Turborotalita quinqueloba, (c) Neogloboquadrina incompta (d) and Neogloboquadrina pachyderma values range from 0 to 80 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$. Station names are indicated above vertical alignments of five dots representing the middle points of the five net sampling intervals.

In the surface waters of the south-to-north transect (0–20 m depth), except at the septentrional station 9, G. uvula was the most abundant species, reaching 64 % of the total fauna at station 7. The second most abundant species, T. quinqueloba, dominated the PF fauna only at station 9, with 45 % of the total assemblage and 26 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$. Both G. uvula and T. quinqueloba have a patchy repartition with two patches of maximum abundances located in the first 0–20 m, at stations 4 and 7 for G. uvula (175 and 245 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$, respectively) and at stations 5 and 7 for T. quinqueloba (53 and 80 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$, respectively). The northern and common patch for both species is the patch with a higher concentration. In the 20–40 m deep layer at station 1 (West of the transect), these two species also showed relatively high concentrations (121 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$ for G. uvula and 30 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$ for T. quinqueloba). The two other major species N. pachyderma and N. incompta had similarly low relative abundances (4 % to 22 % and 8 % to 25 %, respectively). For both species, maximum absolute abundance of about 40–45 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$ occurred in the central part of the transect between 72 and 74 N.

The NMDS analysis of species abundances with regard to environmental parameters (latitude of the station, temperature, salinity, Chl a concentration) indicates that none of the species-specific distribution displays a significant correlation to any of the tested variables (p values >0.1). NMDS documents distributional affinity (Fig. 5), with N. pachyderma and N. incompta plotting in the same area and T. quinqueloba and G. uvula plotting separately from each and also from the N. pachydermaN. incompta area.

Figure 5Standard non-metric multidimensional scaling (NMDS) ordination analysis of planktonic foraminifera species distribution (red dots) with temperature, salinity, Chl a and station location as factors (blue arrows).

## 4.2 Planktonic foraminifera protein biomass

Individual protein content (BCA method) and associated test minimum diameter were measured for a total of 272 specimens of the four major species, including 32 specimens of Neogloboquadrina pachyderma, 58 Neogloboquadrina incompta, 72 Globigerinita uvula and 110 Turborotalita quinqueloba. About 5 to 25 individuals per species were selected at each sampled depth interval of the seven stations along the S–N transect, paying careful attention to sample the whole size range of the populations. At station 7, the protein extraction was successful for only one specimen of N. pachyderma. Therefore, no value is displayed for this species at this station in Fig. 6.

Figure 6Boxplot of the size-normalized protein biomass (SNP, µg µm−1), along the south-to-north transect (from station 3 = 69.8 N to station 9 = 75.6 N) (a) for Globigerinita uvula (b), Turborotalita quinqueloba (c), Neogloboquadrina incompta (d) and Neogloboquadrina pachyderma. Blue dots highlight the data dispersion. Potential outliers such as data for N. pachyderma were removed at station 7 as analyses were only run on one individual. Thick lines indicate median, boxes extend to the interquartile range (IQR) and whiskers indicate 1.5×IQR.

Minimum diameters of the 272 selected tests cover a large size range, from 65 to 315 µm with a median value of 160 µm. N. incompta is the biggest species with a median of 200 µm, and G. uvula is the smallest with a median of 110 µm (Table 2). For each studied species, the mean size is close enough to the median size to say that the size distribution of the picked tests is symmetric, thus making us confident that our test selection properly represents the natural test size range of each studied species. The biomass of a single individual normalized by its test size (SNP) averages out about 0.0055 µg of protein per micrometre of foraminiferal shell diameter. It varies depending on species from 0.0004 (G. uvula) to 0.0426 µg µm−1 (T. quinqueloba). The SNP of all four species displays a northward decrease from 70 to 74 (Fig. 6). T. quinqueloba and G. uvula have slightly (but not significantly) higher relative protein concentrations than N. pachyderma and N. incompta (Fig. 6).

Table 2Descriptive statistics (minimum, maximum, average, and median values in micrometres) of the minimum (i.e. the shortest) diameter of all the 272 BCA measured specimens, per species.

## 4.3 Planktonic foraminifera diversity and distribution in surface sediments

Concentrations of planktonic foraminifera with colourless empty tests varied from a maximum of 6200 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{cm}}^{-\mathrm{3}}$ at station 4 (71.3 N) to a minimum of 200 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{cm}}^{-\mathrm{3}}$ at the septentrional station 9 (Fig. 7a). Neogloboquadrina pachyderma was the most abundant species (31 % to 59 %) along the entire transect. Assemblages were more mixed at the two ends of the transect where N. pachyderma was less abundant. The assemblage at the southernmost point also contained Turborotalita quinqueloba (33 %) and Neogloboquadrina incompta (24 %). At the northernmost point, station 9, T. quinqueloba (23 %) co-occurred with Globigerinita uvula (25 %).

Figure 7Relative species occurrence (percent of the total fauna) of planktonic foraminifera found in the upper 0.5 cm core-top sediment (histograms) and total concentration of individuals per cubic centimetre of dry sediment (logarithmic y axis on the right): (a) colourless empty tests (station numbers are in brackets) and (b) individuals bearing a Rose Bengal coloured cytoplasm.

Concentrations of planktonic foraminifera bearing a coloured cytoplasm (Fig. 7b) varied from 100 to 300 $\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{cm}}^{-\mathrm{3}}$. All along the transect, the relative abundance of coloured N. pachyderma remained between 10 % and 26 %. The species T. quinqueloba occurred everywhere above 20 % and up to 40 % south of 72 N. The central station 6 was dominated by G. uvula (38 %). North of 74, the fauna was balanced between N. incompta (33 % and 9 %) and G. uvula (8 % and 34 %).

5 Discussion

## Distribution pattern of living planktonic foraminifera at the Barents Sea Opening

In the late summer of 2014 the hydrology at the BSO was characterized by a strong water stratification with a 30 m thick Chl a enriched lens of NwCW that to the north overlapped with the NwAW (saltier and colder) from 69.8 to 74.5 N. Further north, a well-mixed water column with characteristics of the NwAW occupied the Storfjordrenna Trough. Here a coccolithophore bloom and the highest concentration of cyanobacteria were recorded in the upper water column (Giraudeau et al., 2016). Despite these marked features the pattern of planktonic foraminifera abundance did not correlate with any of the studied environmental parameters (Fig. 5). These observations confirm the low influence of commonly imputed parameters such as temperature, salinity and primary production on PF density (Schiebel et al., 2017). In accordance with the conclusions of Retailleau et al. (2018), there are multiple indications of a possible role of water turbidity in PF abundance variation. The highest densities of PF occurred in the 0–20 m upper water layer between 70.5 and 74.5 N. Their very low abundances (total concentration $<\mathrm{5}\phantom{\rule{0.125em}{0ex}}\mathrm{ind}.\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{3}}$) below 100 m depth suggest a shallow depth habitat for individuals in the region, especially for N. pachyderma, which was recently reported to live between 25 and 280 m depth in the north Atlantic Arctic region (Greco et al., 2019). Very low abundances were also recorded close to the Norwegian and Spitsbergen coasts. The low abundances at the two ends of the studied transect could reflect a planktonic foraminifera patchiness pattern of distribution (Meilland et al., 2019) or highlight the fact that waters under continental influences (nutrient-enriched, more turbid) likely hamper the foraminiferal production. In line with this, the abrupt decrease in abundances from west to east (stations 2 to 6 to 1) may be ascribed to the decrease in depth of the Bjømøyrenna Trough up to the Barents Sea shelf (from 1850 to 430 m), as foraminifera are suspected to avoid neritic waters over continental shelves (Schmuker, 2000).

The remarkable point of our results is the dominance of Globigerinita uvula. This species, described as a temperate to polar species (Schiebel and Hemleben, 2017), is known to account for less than 2 % of the assemblages in marginal Arctic seas based on material collected with a 63 µm plankton net mesh size (Volkmann, 2000). Neogloboquadrina pachyderma is considered the dominant species in polar regions, making up more than 90 % of the total planktonic foraminifera assemblages (e.g. Schiebel et al., 2017). The high densities of G. uvula recorded at the BSO in 2014 seem to contradict the former statements but are consistent with a recent study reporting G. uvula as one of the dominant species in southern high latitudes, south of the Polar Front (Meilland et al., 2017). A possible explanation could be the warming experienced by the western Barents Sea (sea surface temperature (SST) anomalies $\approx +\mathrm{2}\phantom{\rule{0.125em}{0ex}}{}^{\circ }\mathrm{C}$) and its increase in salinity (sea surface salinity (SSS) anomalies $\approx +\mathrm{0.3}$) over the last decades (Dobrynin and Pohlmann, 2015). These hydrological changes impact the plankton dynamics and biogeography, with a northward shift of the natural range of biological communities (Barton et al., 2016). Thus, the species distribution of planktonic foraminifera could be affected by an eventual expansion of subpolar/temperate species towards high latitudes, leading to phytoplankton composition changes in response to sea temperature warming under global climate change. Our observations from the north polar region support the shift of planktonic foraminifera assemblages to warmer conditions already asserted from North Atlantic data (Jonkers et al., 2019) and from southern Indian Ocean data (Meilland et al., 2017). However, a single observational dataset in the Barents Sea is not sufficient to robustly validate this assumption and a second hypothesis for the dominance of G. uvula in our sampling area could be a response to specific phytoplankton composition and ambient water conditions by pulsed reproduction events only in summer conditions. This seasonal pattern is known to occur in polar regions for Turborotalita quinqueloba (Schiebel and Hemleben, 2017). In fact, this species is the second most dominant one in our late summer 2014 samples. As observed in this study, T. quinqueloba is also known to display high concentrations in the Barents Sea and western Spitsbergen (Volkmann, 2000) and to co-occur with the typically polar species Neogloboquadrina pachyderma in the high-latitude cold-water assemblages (Volkmann, 2000; Eynaud, 2011).

A discrepancy between the species-specific distribution patterns was observed in late summer 2014 at the BSO. The low abundances of Neogloboquadrina pachyderma and Neogloboquadrina incompta consistent over the studied area versus the patchy distribution and high densities of Globigerinita uvula and Turborotalita quinqueloba suggest differences in the ecological strategy and behaviour between these two pairs of species. The patchy pattern of planktonic foraminifera distribution has been observed before (Boltovskoy, 1971; Siccha et al., 2012; Meilland et al., 2019), suggesting that high densities are not exclusively constrained by the physical structure of the (sub-)surface layers.

Potential differences in diet preferences could explain the observed species distribution in late summer 2014 at the BSO. Both G. uvula and T. quinqueloba are thought to follow food availability and primary production (Volkmann, 2000; Schiebel and Hemleben, 2017). However, we observed no correlation between their distribution and Chl a concentrations (Fig. 5). In late summer 2014, G. uvula and T. quinqueloba showed high concentrations, especially at station 7, located at the crossroads of the Atlantic (NwAW) and Arctic waters flowing out of Storfjordrenna (Fig. 1), at the edge of the Polar Front (Oziel et al., 2017). From this location to the northern station, the concentration of phytoplankton was relatively low and the phytoplankton community showed singular characteristics, in comparison to the southern part of the transect, fuco-flagellates became dominant and diatom concentrations decreased. The fuco-flagellate blooms (mainly Phaeocystis pouchetii in late summer 2014; Giraudeau et al., 2016) are well known to occur in the Barents Sea (Wassmann et al., 1990; Vaquer-Sunyer et al., 2013). Our hypothesis is thus that high densities of G. uvula and T. quinqueloba are due to food composition (quality) rather than food concentrations (quantity). This also implies that satellite-derived chlorophyll concentrations, considered a potential indicator of algal bloom, may not always be good indicators of foraminiferal concentration and distribution.

### Planktonic foraminifera protein concentration, potential marker of their metabolism

Proteins are the main component of zooplankton biomass (Corg) in all oceanographic regions, from the tropics to polar areas (e.g. Percy and Fife, 1981; Donnelly et al., 1994; Kumar et al., 2013; Yun et al., 2015). Their role is essential to organisms' growth and their concentration and composition likely reflect the environment individuals grew in and how well they adjust to it. Based on previous studies, the protein concentration of PF can be used as a proxy of its biomass (Corg), and foraminiferal biomass should remain the same for a given test size (Schiebel and Movellan, 2012). However, in our study, the SNP (size-normalized protein content) of PF decreases with higher latitude and hence with decreased Chl a concentration and temperature. The size of individuals picked for these analyses remains however constant. This observation suggests that foraminifera metabolisms (i.e. ability to consume and degrade food and to grow) are decreasing towards the north. This would be consistent with the observation of lower metabolism for zooplankton with decreasing salinity and temperature in the Arctic (Alcaraz et al., 2010). Proteins are the main component of zooplankton biomass (Corg), closely followed by lipids. Lipids in zooplankton organisms are very variable geographically, showing a latitudinal pattern with high percentages in polar areas and low percentages in warm tropical waters. Lipid percentages also display seasonal features, with higher values in summer than in winter (Falk-Petersen et al., 1999; Mayzaud et al., 2011; Kumar et al., 2013). It is thus possible that a part of the energy (biomass Corg) of the PF collected along the south-to-north transect shifts from being stocked as protein in warmer waters to being stocked as lipids in colder waters. This strategy would allow foraminifera to resist the cold to potentially overwinter. This hypothesis is supported by analogous observations made on different size fraction of zooplankton in the southern Indian Ocean showing variability in protein and lipid percentages among the 80 to 200 µm populations (Harmelin-Vivien et al., 2019) and also by observations made on pteropods in the Arctic (Kattner et al., 1998; Phleger et al., 2001; Böer et al., 2005). The fact that higher SNPs of foraminifera were observed where Chl a is higher is compatible with the fact that polar organisms rely on their protein catabolism when food is easily accessible rather than on their lipid storage (Brockington and Clarke, 2001). It has also been shown that a single organism in a cold environment is able to switch between predominantly protein or lipid catabolism across its life (Mayzaud, 1976). This suggests that individuals from the same species can display more or fewer proteins for the same biomass in different locations. With the reduction of PF protein concentration (and likely metabolism) going north, one could expect lower abundances. However, we observe no link between PF concentrations, which appeared to be species specific, and protein concentrations (evolving similarly for all four species), suggesting a decoupling between individuals' metabolism and densities.

### Discrepancy between upper water column and interface sediment samples

The PF species compositions recorded during the late summer 2014 in the water column and in surface sediments are similar while species relative abundances are drastically different. Indeed, the living fauna (collected by plankton net) displays large relative and absolute abundances of the two species Globigerinita uvula and Turborotalita quinqueloba, whereas the fossil assemblages (found in core tops) are largely dominated by Neogloboquadrina pachyderma or, at the southernmost station, co-dominated by T. quinqueloba and N. pachyderma (Fig. 7a). Affected by differential settling velocities (200 and 500 m d−1 in normal conditions), water depth and test sizes of different species, the foraminiferal fluxes exported from the upper productive surface and reaching the sea bottom depend on direction and intensity of currents (Takahashi and Bé, 1984). Lateral advection may transport shells over long distances >25 km for N. pachyderma and >50km for T. quinqueloba (Von Gyldenfeldt et al., 2000). Lateral advection of shells is also strengthened by water stratification that increases resident time at the shear boundary between water masses (Kuhnt at al., 2013). The BSO has a complex hydrography with the buoyant NwCW flowing northwards above the NwAW that is entering the Barents Sea eastwards when cold Barents Sea Atlantic Water (BSAW) is flowing westwards. In such areas, the PF settling velocities and extension of lateral advection are poorly constrained. Consequently, the sediment core records cannot match exactly with the place and/or the intensity of production (Von Gyldenfeldt et al., 2000; Jonkers et al., 2015). However, considering potential lateral advection of shells, Pados and Spielhagen (2014) observed from a study through the dynamic Fram Strait that the distribution pattern obtained by plankton tows was clearly reflected on the sediment surface and that the assemblage on the sediment surface can be used as an indicator for modern planktonic foraminiferal fauna. This suggests that the large discrepancy between upper water column and interface sediment samples collected at the BSO in late summer 2014 should be taken into consideration. As a sedimentation rate of $\mathrm{1.3}±\mathrm{0.6}\phantom{\rule{0.125em}{0ex}}\mathrm{mm}\phantom{\rule{0.125em}{0ex}}{\mathrm{yr}}^{-\mathrm{1}}$ has been recently measured in the Storfjordrenna outlet (76 N–17 E), close to our station 9 (Fossile et al., 2019), the core-top sediments may have recorded less than a decade. Even though sedimentation rates are likely to vary along the transect, we hypothesise that the sea surface–bottom differences in the foraminiferal assemblages along the south-to-north transect at the BSO might reflect a community change within a short period of time (about a decade).

Furthermore, the analysis of sediment from the five core tops demonstrated important differences between the assemblages of fossil fauna and recently settled tests (likely coming from surface spring–summer production), i.e. Rose Bengal-stained tests bearing not yet decomposed cytoplasm. For example, at 71.3 N, the percentages of coloured T. quinqueloba and G. uvula are twice as high as the ones observed for the fossil faunas (Fig. 7). At 72.9 N in the surface sediment, G. uvula reaches up to 38 % of the coloured assemblages (Fig. 7b), whereas it never exceeds 25 % in the non-coloured ones (Fig. 7a). The high abundance of the two species G. uvula and T. quinqueloba in the living fauna as well as in the recently settled shells but not in the fossil faunas suggest that they may present a seasonal character with a production period focussed in late summer as a response to environmental and trophic conditions. This is supported by previous studies in the Arctic where T. quinqueloba has been found to dominate assemblages sampled in August (Carstens et al., 1997; Volkmann, 2000) but not in June–early July (Pados and Spielhagen, 2014) and by sediment trap observations from the subpolar North Atlantic where T. quinqueloba reaches its maximum in autumn (Jonkers et al., 2010). The dominance of N. pachyderma in the fossil faunas collected at the BSO and its low but constant presence in the coloured shells of surface sediment and plankton tow sampled in late summer 2014 suggest that this species may demonstrate a regular production throughout the year. Neogloboquadrina pachyderma production appears to be sustained and constant year-round in the area, whereas other species clearly respond to a local seasonal signal.

6 Conclusions

The sampling and analytic approaches deployed during the MOCOSED14 cruise and combining the use of plankton net, core-top, molecular biology (protein measurement), environmental parameters and phytoplankton characterization provide us with a unique dataset to better constrain the distribution of planktonic foraminifera within the highly complex studied area of the western Barents Sea.

The observed abundances of PF in the studied area are high offshore and the lower densities were recorded close to the Norwegian and Spitsbergen coasts. These observations highlight the fact that waters under continental influences (nutrient-enriched, more turbid) are rather inhospitable for PF production. The PF species composition observed at the BSO is diverse, with more than 10 different species in the net samples including Globigerinita uvula (45 %), Turborotalita quinqueloba (26.2 %), Neogloboquadrina incompta (15 %) and Neogloboquadrina pachyderma (8.9 %). The two species G. uvula and T. quinqueloba dominate the living (water sample) population and display highly patchy abundances, suggesting they occur in late summer in response to physico-chemical conditions and related specific primary productivity. The dominance of G. uvula in water samples could also be a signal of the temperature increase experienced over the last decades in the Barents Sea and the North Atlantic Ocean. Further sampling in the area is thus needed to test this hypothesis. The species N. pachyderma and N. incompta show low densities but a continuous distribution pattern in the water samples. They also dominate the core-top assemblages, suggesting that both species present a more consistent production over the course of spring–summer seasons.

Unlike their species-specific abundance pattern of distribution, size-normalized protein concentrations of all four major species decrease with increasing latitude (and a decrease in temperature and Chl a concentration). This observation leads us to hypothesize that (1) PF abundance and metabolism are decoupled and (2) foraminifera metabolism in the north of the studied region is lower than in the south. It opens the following question: can individuals of the same species balance the ratio between their protein and lipid concentrations (major components of zooplankton Corg) in order to adapt to environmental conditions (e.g. temperature)? Further analyses on planktonic foraminifera lipid concentration and composition are thus needed and would help us to better understand the metabolism of these organisms and their fate in a context of climate change.

Data availability
Data availability.

Data will be made available on request to the main author until their online publication on https://doi.pangaea.de/10.1594/PANGAEA.913782 (Meilland and Howa, 2020).

Author contributions
Author contributions.

JM and HH designed the study. JM, VH, ID and JS generated the data and carried out the analyses. TG provided access to the MOCOSED 2014 cruise. All authors contributed to writing the manuscript.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

The authors are thankful to the crew and captain of R/V Pourquoi pas?, as well as the scientific participants from the SHOM (French Hydrographic Office) for their support in sampling during the cruise “MOCOSED 2014”. We fully acknowledge the efficient technical help provided by Sophie Sanchez at the LPG-BIAF laboratory. We thank the three anonymous reviewers and Lukas Jonkers for their help on improving the manuscript. A grant was kindly provided to Vivien Hulot by the SHOM to study phytoplankton and PF communities.

Financial support
Financial support.

The article processing charges for this open-access publication were covered by the University of Bremen.

Review statement
Review statement.

This paper was edited by Carol Robinson and reviewed by three anonymous referees.

References

Alcaraz, M., Almeda, R., Calbet, A., Saiz, E., Duarte, C. M., Lasternas, S., Agustí, S., Santiago, S., Movilla, J., and Alonso, A.: The role of arctic zooplankton in biogeochemical cycles: respiration and excretion of ammonia and phosphate during summer, Polar Biol., 33, 1719–1731, 2010.

AMAP Assessment 2018: Arctic Ocean Acidification, Arctic Monitoring and Assessment Programme (AMAP), Tromsø, Norway, available at: https://www.amap.no, last access: 10 October 2018.

Årthun, M., Eldevik, T., Smedsrud, L. H., Skagseth, Ø., and Ingvaldsen, R. B.: Quantifying the Influence of Atlantic Heat on Barents Sea Ice Variability and Retreat, J. Clim., 25, 4736–4743, https://doi.org/10.1175/JCLI-D-11-00466.1, 2012.

Barton, A. D., Irwin, A. J., Finkel, Z. V., and Stock, C. A.: Anthropogenic climate change drives shift and shuffle in North Atlantic phytoplankton communities, P. Natl. Acad. Sci. USA, 113, 2964–2969, https://doi.org/10.1073/pnas.1519080113, 2016.

Beaugrand, G., McQuatters-Gollop, A., Edwards, M., and Goberville, E.: Long-term responses of North Atlantic calcifying plankton to climate change, Nat. Clim. Change, 3, 263–267, 2013.

Beer, C. J., Schiebel, R., and Wilson, P. A.: Technical Note: On methodologies for determining the size-normalised weight of planktic foraminifera, Biogeosciences, 7, 2193–2198, https://doi.org/10.5194/bg-7-2193-2010, 2010.

Böer, M., Gannefors, C., Kattner, G., Graeve, M., Hop, H., and Falk-Petersen, S.: The Arctic pteropod Clione limacina: seasonal lipid dynamics and life-strategy, Mar. Biol., 147, 707–717, 2005.

Bollmann, J., Quinn, P. S., Vela, M., Brabec, B., Brechner, S., Cortés, M. Y., Hilbrecht, H., Schmidt, D. N., Schiebel, R., and Thierstein, H. R.: Automated particle analysis: calcareous microfossils, edited by: Francus, P., Image Analysis, Sediments and Paleoenvironments, Springer, Dordrecht, chap. 12, 229–252, 2004.

Boltovskoy, E.: Patchiness in the distribution of planktonic foraminifera, in: Proceedings of the Second Planktonic Conference, Rome 1970, edited by: Farinacci, A., Edizioni Technoscienza, Roma, 107–116, 1971.

Brockington, S. and Clarke, A.: The relative influence of temperature and food on the metabolism of a marine invertebrate, J. Exp. Mar. Biol. Ecol., 258, 87–99, 2001.

Burenkov, V. I., Kopelevich, O. V., Rat'kova, T. N., and Sheberstov, S. V.: Satellite observations of the coccolithophorid bloom in the Barents Sea, Oceanology, 51, 766–774, https://doi.org/10.1134/s0001437011050043, 2011.

Carstens, J., Hebbeln, D., and Wefer, G.: Distribution of planktic foraminifera at the ice margin in the Arctic (Fram Strait), Mar. Micropaleontol., 29, 257–269, 1997.

Clayton, C. R. I., Abbireddy, C. O. R., and Schiebel, R.: A method of estimating the form of coarse particulates, Geotechnique, 59, 493–501, https://doi.org/10.1680/geot.2007.00195, 2009.

Dai, A., Luo, D., Song, M., and Liu, J.: Arctic amplification is caused by sea-ice loss under increasing CO2, Nat. Commun., 10, 1–13, https://doi.org/10.1038/s41467-018-07954-9, 2019.

Dalpadado, P., Ingvaldsen, R. B., Stige, L. C., Bogstad, B., Knutsen, T., Ottersen, G., and Ellertsen, B.: Climate effects on Barents Sea ecosystem dynamics, ICES J. Mar. Sci., 69, 1303–1316, https://doi.org/10.1093/icesjms/fss063, 2012.

Dobrynin, M. and Pohlmann, T.: Anomalous hydrographic conditions in the western Barents Sea observed in March 2014, Cont. Shelf Res., 111, 69–82, https://doi.org/10.1016/j.csr.2015.10.020, 2015.

Donnelly, J., Torres, J. J., Hopkins, T. L., and Lancraft, T. M.: Chemical composition of Antarctic zooplankton during austral fall and winter, Polar Biol., 14, 171–183, 1994.

Dylmer, C. V., Giraudeau, J., Eynaud, F., Husum, K., and De Vernal, A.: Northward advection of Atlantic water in the eastern Nordic Seas over the last 3000 yr, Clim. Past, 9, 1505–1518, https://doi.org/10.5194/cp-9-1505-2013, 2013.

Eynaud, F.: Planktonic foraminifera in the Arctic: potentials and issues regarding modern and quaternary populations, IOP C Ser. Earth Env., 14, O12005, https://doi.org/10.1088/1755-1315/14/1/012005, 2011.

Falk-Petersen, S., Sargent, J. R., Lønne, O. J., and Timofeev, S.: Functional biodiversity of lipids in Antarctic zooplankton: Calanoides acutus, Calanus propinquus, Thysanoessa macrura and Euphausia crystallorophias, Polar Biol., 21, 37–47, 1999.

Fossheim, M., Primicerio, R., Johannesen, E., Ingvaldsen, R. B., Aschan, M. M., and Dolgov, A. V.: Recent warming leads to a rapid borealization of fish communities in the Arctic, Nat. Clim. Change, 5, 673–677, https://doi.org/10.1038/nclimate2647, 2015.

Fossile, E., Nardelli, M. P., Jouini, A., Lansard, B., Pusceddu, A., Moccia, D., Michel, E., Péron, O., Howa, H., and Mojtahid, M.: Benthic foraminifera as tracers of brine production in Storfjorden “sea ice factory”, Biogeosciences Discuss., https://doi.org/10.5194/bg-2019-405, in review, 2019.

Giraudeau, J., Hulot, V., Hanquiez, V., Devaux, L., Howa, H., and Garlan, T.: A survey of the summer coccolithophore community in the western Barents Sea, J. Mar. Syst., 158, 93–105, https://doi.org/10.1016/j.jmarsys.2016.02.012, 2016.

Greco, M., Jonkers, L., Kretschmer, K., Bijma, J., and Kucera, M.: Depth habitat of the planktonic foraminifera Neogloboquadrina pachyderma in the northern high latitudes explained by sea-ice and chlorophyll concentrations, Biogeosciences, 16, 3425–3437, https://doi.org/10.5194/bg-16-3425-2019, 2019.

Harmelin-Vivien, M., Bănaru, D., Dromard, C. R., Ourgaud, M., and Carlotti, F.: Biochemical composition and energy content of size-fractionated zooplankton east of the Kerguelen Islands, Polar Biol., 42, 603–617, 2019.

Holding, J. M., Duarte, C. M., Sanz-Martín, M., Mesa, E., Arrieta, J. M., Chierici, M., Hendriks, I. E., 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, https://doi.org/10.1038/nclimate2768, 2015.

Jonkers, L., Brummer, G.-J. A., Peeters, F. J. C., van Aken, H. M., and de Jong, M. F.: Seasonal stratification, shell flux, and oxygen isotope dynamics of left-coiling N. pachyderma and T. quinqueloba in the western subpolar North Atlantic, Paleoceanography, 25, PA2204, https://doi.org/10.1029/2009PA001849, 2010.

Jonkers, L., van Heuven, S., Zahn, R., and Peeters, F. J. C.: Seasonal patterns of shell flux, δ18O and δ13C of small and large N. pachyderma (s) and G. bulloides in the subpolar North Atlantic, Paleoceanogr. Paleoclimatol., 28, 164–174, https://doi.org/10.1002/palo.20018, 2013.

Jonkers, L., Reynolds, C. E., Richey, J., and Hall, I. R.: Lunar periodicity in the shell flux of planktonic foraminifera in the Gulf of Mexico, Biogeosciences, 12, 3061–3070, https://doi.org/10.5194/bg-12-3061-2015, 2015.

Jonkers, L., Hillebrand, H., and Kucera, M.: Global change drives modern plankton communities away from the pre-industrial state, Nature, 570, 372–375, https://doi.org/10.1038/s41586-019-1230-3, 2019.

Kattner, G., Hagen, W., Graeve, M., and Albers, C.: Exceptional lipids and fatty acids in the pteropod Clione limacina (Gastropoda) from both polar oceans, Mar. Chem., 61, 219–228, 1998.

Kucera, M., Weinelt, M., Kiefer, T., Pflaumann, U., Hayes, A., Weinelt, M., Chen, M.-T., Mix, A. C., Barrows, T. T., Cortijo, E., Duprat, J., Juggins, S., and Waelbroeck, C.: Reconstruction of sea-surface temperatures from assemblages of planktonic foraminifera: Multi-technique approach based on geographically constrained calibration data sets and its application to glacial Atlantic and Pacific Oceans, Quaternary Sci. Rev., 24, 951–998, 2005.

Kuhnt, T., Howa, H., Schmidt, S., Marié, L., and Schiebel, R.: Flux dynamics of planktic foraminifer tests at a hemipelagic site of the inner Bay of Biscay (Northeast Atlantic margin), J. Mar. Syst., 109–110, S169–S181, https://doi.org/10.1016/j.jmarsys.2011.11.026, 2013.

Kumar, M. A., Padmavati, G., and Anandavelu, I.: Biochemical composition and calorific value of zooplankton from the coastal waters of South Andaman, Proc. Int. Acad. Ecol. Environ. Sc., 3, 278–287, 2013.

Loeng, H.: Features of the physical oceanographic conditions of the Barents Sea, Polar Res., 10, 5–18, 1991.

Mackey, M. D., Mackey, D. J., Higgins, H. W., and Wright, S. W.: CHEMTAX – a program for estimating class abundances from chemical markers: application to HPLC measurement of phytoplankton, Mar. Ecol. Prog. Ser., 144, 265–283, 1996.

Mayzaud, P., Lacombre, S., and Boutoute, M.: Seasonal and growth stage changes in lipid and fatty acid composition in the multigeneration copepod Drepanopus pectinatus from Iles Kerguelen, Antarct Sci., 23, 3–17, 2011.

Mayzaud, P.: Respiration and nitrogen excretion of zooplankton. IV. The influence of starvation on the metabolism and the biochemical composition of some species, Mar. Biol., 37, 47–58, 1976.

Meilland, J. and Howa, H.: Barents Sea plankton tows and surface sediment planktonic foraminifera assemblages, PANGAEA, https://doi.org/10.1594/PANGAEA.913782, 2020.

Meilland, J., Howa, H., LoMonaco, C., and Schiebel, R.: Individual planktic foraminifer protein-biomass affected by trophic conditions in the Southwest Indian Ocean, 30 S–60 S, Mar. Micropaleontol., 124, 63–74, https://doi.org/10.1016/j.marmicro.2016.02.004, 2016.

Meilland, J., Schiebel, R., Monaco, C. L., Sanchez, S., and Howa, H.: Abundances and test weights of living planktic foraminifers across the Southwest Indian Ocean: Implications for carbon fluxes, Deep-Sea Res. Pt. I, 131, 27–40, 2018.

Meilland, J., Siccha, M., Weinkauf, M. F. G., Jonkers, L., Morard, R., Baranowski, U., Baumeister, A., Bertlich, J., Brummer, G. J., Debray, P., Fritz-Endres, T., Groeneveld, J., Magerl, L., Munz, P., Rillo, M. C., Schmidt, C., Takagi, H., Theara, G., and Kucera, M.: Highly replicated sampling reveals no diurnal vertical migration but stable species-specific vertical habitats in planktonic foraminifera, J. Plank. Res., 41, 127–141, 2019.

Neukermans, G., Oziel, L., and Babin, M.: Increased intrusion of warming Atlantic water leads to rapid expansion of temperate phytoplankton in the Arctic, Global Change Biol Wiley, 24, 2545–2553, https://doi.org/10.1111/gcb.14075, 2018.

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O'Hara, R. B., Simpson, G. L., Solymos, P., Stevens, M. H. H., and Wagner, H.: Vegan: Community Ecology Package, available at: http://www.CRAN.R-project.org/package=vegan, 2013.

Oziel, L., Neukermans, G., Ardyna, M., Lancelot, C., Tison, J-L., Wassmann, P., Sirven, J., Ruiz-Pino, D., and Gascard, J.-C.: Role for Atlantic inflows and sea ice loss on shifting phytoplankton blooms in the Barents Sea, J. Geophys. Res.-Ocean., 122, 5121–5139, https://doi.org/10.1002/2016JC012582, 2017.

Pados, T. and Spielhagen, R. F.: Species distribution and depth habitat of recent planktic foraminifera in Fram Strait, Arctic Ocean. Polar Res., 33, https://doi.org/10.3402/polar.v33.22483, 2014.

Percy, J. A. and Fife, F. J.: The biochemical composition and energy content of Arctic marine macrozooplankton, Arctic, 307–313, 1981.

Phleger, C. F., Nelson, M. M., Mooney, B. D., and Nichols, P. D.: Interannual variations in the lipids of the Antarctic pteropods Clione limacina and Cliopyramidata, Comp. Biochem. Phys. B, 128, 553–564, 2001.

Retailleau, S., Eynaud, F., Mary, Y., Schiebel, R., and Howa, H.: Canyon heads and river plumes: how might they influence neritic planktonic foraminifera communities in the SE Bay of Biscay?, J. Foramin. Res., 42, 257–269, 2012.

Schiebel, R. and Hemleben, C.: Planktic Foraminifers in the Modern Ocean, Springer-Verlag GmbH, Berlin, 357 pp., https://doi.org/10.1007/978-3-662-50297-6, 2017.

Schiebel, R., Spielhagen, R. F., Garnier, J., Hagemann, J., Howa, H., Jentzen, A., Martínez-Garcia, A., Meilland, J., Michel, E., Repschläger, J., Salter, I., Yamasakig, M., and Haug, G.: Modern planktic foraminifers in the high-latitude ocean, Mar. Micropaleontol., 136, 1–13, 2017.

Schiebel, R. and Movellan, A.: First-order estimate of the planktic foraminifer biomass in the modern ocean, Earth Syst. Sci. Data, 4, 75–89, https://doi.org/10.5194/essd-4-75-2012, 2012.

Schlitzer, R.: Ocean data view, available at: http://odv.awi.de, 2014.

Schmuker, B.: The influence of shelf vicinity on the distribution of planktic foraminifera south of Puerto Rico, Mar. Geol., 166, 125–143, 2000.

Schönfeld, J., Golikova, E., Korsun, S., and Spezzaferri, S.: The Helgoland Experiment – assessing the influence of methodologies on Recent benthic foraminiferal assemblage composition, J. Micropalaeontol., 32, 161–182, https://doi.org/10.1144/jmpaleo2012-022, 2013.

Shepherd, T. G.: Effects of a warming Arctic, Science, 353, 989–990, https://doi.org/10.1126/science.aag2349, 2016.

Siccha, M. and Kucera, M.: ForCenS, a curated database of planktonic foraminifera census counts in marine surface sediment samples, Sci. Data, 4, 170109, https://doi.org/10.1038/sdata.2017.109, 2017.

Siccha, M., Schiebel, R., Schmidt, S., and Howa, H.: Short-term and small scale variability in planktic foraminifera test flux in the Bay of Biscay, Deep-Sea Res. Pt. I, 64, 146–156, https://doi.org/10.1016/j.dsr.2012.02.004, 2012.

Skagseth, Ø., Furevik, T., Ingvaldsen, R., Loeng, H., Mork, K. A., Orvik, K. A., and Ozhigin, V.: Volume and heat transports to the Arctic Ocean via the Norwegian and Barents Seas, in: Arctic–Subarctic Ocean Fluxes, Springer, Dordrecht, 45–64, 2008.

Smedsrud, L. H., Esau, I., Ingvaldsen, R. B., Eldevik, T., Haugan, P. M., Li, C., Lien, V. S., Olsen, A., Omar, A. M., Otterå, O. H., Risebrobakken, B., Sandø, A. B., Semenov, V. A., and Sorokina, S. A.: The role of the Barents Sea in the Arctic climate system, Rev. Geophys., 51, 415–449, 2013.

Smyth, T.J., Tyrell, T., and Tarrant, B.: Time series of coccolithophore activity in the Barents Sea, from twenty years of satellite imagery, Geophys. Res. Lett., 31, L11302, https://doi.org/10.1029/2004GL019735, 2004.

Takahashi, K. and Bé, A. W. H.: Planktonic foraminifera: factors controlling sinking speeds, Deep-Sea Res., 31, 1477–1500, https://doi.org/10.1016/0198-0149(84)90083-9, 1984.

Troupin, C., Barth, A., Sirjacobs, D., Ouberdous, M., Brankart, J. M., Brasseur, P., Rixen, M., Alvera-Azcárate, A., Belounis, M., Capet, A., Lenartz, F., Toussaint, M. E., and Beckers, J. M.: Generation of analysis and consistent error fields using the Data Interpolating Variational Analysis (DIVA), Ocean Model., 52–53, 90–101, 2012.

Van Heukelem, L. and Thomas, C. S.: Computer-assisted high performance liquid chromatography method development with applications to the isolation and analysis of phytoplankton pigments, J. Chromatogr., 910, 31–49, 2001.

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.

Volkmann, R.: Planktic foraminifers in the outer Laptev Sea and the Fram Strait modern distribution and ecology, J. Foramini. Res., 30, 157–176, 2000.

Von Gyldenfeldt, A. B., Carstens, J., and Meincke, J.: Estimation of the catchment area of a sediment trap by means of current meters and foraminiferal tests, Deep-Sea Res Pt. II, 47, 1701–1717, 2000.

Wassmann, P., Vernet, M., Mitchell, B. G., and Reyh, F.: Mass sedimentation of Phaeocystis pouchettii in the Barents Sea, Mar. Ecol. Prog. Ser., 66, 183–195, 1990.

Yun, M. S., Lee, D. B., Kim, B. K., Kang, J. J., Lee, J. H., Yang, E. J., Park, W. G., Chung, K. H., and Lee, S. H.: Comparison of phytoplankton macromolecular compositions and zooplankton proximate compositions in the northern Chukchi Sea, Deep-Sea Res. Pt. II, 120, 82–90, 2015.