Abstract
The Arctic marine biome, shrinking with increasing temperature and receding sea-ice cover, is tightly connected to lower latitudes through the North Atlantic. By flowing northward through the European Arctic Corridor (the main Arctic gateway where 80% of in- and outflow takes place), the North Atlantic Waters transport most of the ocean heat, but also nutrients and planktonic organisms toward the Arctic Ocean. Using satellite-derived altimetry observations, we reveal an increase, up to two-fold, in North Atlantic current surface velocities over the last 24 years. More importantly, we show evidence that the North Atlantic current and its variability shape the spatial distribution of the coccolithophore Emiliania huxleyi (Ehux), a tracer for temperate ecosystems. We further demonstrate that bio-advection, rather than water temperature as previously assumed, is a major mechanism responsible for the recent poleward intrusions of southern species like Ehux. Our findings confirm the biological and physical “Atlantification” of the Arctic Ocean with potential alterations of the Arctic marine food web and biogeochemical cycles.
Similar content being viewed by others
Introduction
The European Arctic Corridor (EAC, Fig. 1) is the main Arctic gateway where 80% of both in- and outflow takes place1. One of the most prominent synoptic feature of the EAC is the northward flowing North Atlantic Waters (NAW) entering the Arctic Ocean through two main branches (indicated in red in Fig. 1), which separate around 70°N. One branch keeps flowing northward toward Fram Strait while the second one turns eastward into the Barents Sea (BS). Over the last three decades, in a context of amplified warming2 and sea-ice loss3, substantial changes have been documented in the NAW inflow with a two-fold increase of its volume occupation and a northward shift of the polar front structure in the BS4. The NAW inflow largely controls the physical and sea-ice conditions of the region5,6,7,8. In addition, almost 50% of the total Arctic primary production takes place in the EAC1,9,10,11, which is of major importance for fisheries12.
Recent warming of the EAC related to the NAW inflow has been suspected to trigger poleward intrusions of temperate phytoplankton13,14,15,16 and species from higher trophic levels17,18,19,20,21. By carrying biomass and nutrients produced elsewhere, bio-advection has recently been proposed as an “essential mechanism” for ecosystem dynamics in the Arctic Ocean22,23. Although the actual role of advection has already been identified as a potential driver altering zooplankton dynamics24,25,26, it has never been assessed quantitatively from observations for phytoplankton in the EAC, or more generally in the Arctic Ocean22,23,27,28,29.
Here we investigate how ocean currents control the spatial dynamics of a specific coccolithophore bloom-forming species, Emiliania huxleyi (hereafter referred to as Ehux). Ehux is usually associated with the surface layer of the temperate NAW in summer, typically in a post-spring bloom context characterized by low nutrients, low light and strong stratification30,31,32. Since Ehux does not form winter resting spores, this tracer of Atlantic ecosystem is generally considered to be a summer visitor in the BS, unlike neritic diatoms species13. A combination of bottom-up (i.e., winter darkness, cold temperatures and intense vertical mixing; see Supplementary Note 2.6) and top-down controls (i.e., grazing, viral lysis33,34,35) prevent Ehux from year-to-year survival in the north-easternmost parts of the BS. Because of their high abundance in Ehux (i.e., 115 × 106 cells L−1 36), coastal regions and fjords of the Norwegian Sea have been suspected to be the source of Ehux for the whole EAC37.
Using a Lagrangian tracking approach based on satellite-derived current fields, we explore how the advection of Ehux cells from these upstream coastal temperate regions shapes the distribution of the massive Ehux blooms in the BS (i.e., via the “seeding effect”38,39). By combining satellite-derived altimetry with ocean-color Particulate Inorganic Carbon (PIC, a coccolithophore biomass proxy) estimates, we demonstrate a major role of bio-advection in phytoplankton transport along the EAC.
Results and discussion
Climatic modes and increasing NAW surface current velocities
To document the interannual and decadal variability of surface NAW currents, we first performed an Empirical Orthogonal Function (EOF) analysis of the Sea Level Anomalies (1993–2016, see Material and Methods), on which the seasonal signal was removed. The EOF analysis showed that the two first modes of variability accounted for more than 80% of the non-seasonal variability (Fig. 2a; Supplementary Note 2.3). When summing the two first modes of variability, a dipole structure emerged, centered on the NAW path between the Barents and Norwegian Sea shelves (about east of 5°E) and the center of the Norwegian Sea (about west of 5°E, Fig. 2a). The associated time-series showed a linear positive trend of about +9 (normalized units, Fig. 2b) over the last 24 years with high energy at the interannual and decadal time scale (Supplementary Fig. 7). This trend toward a more positive phase corresponds to an increase of the sea level gradient across the NAW path, hence intensifying NAW currents, while negative phases are associated with weaker NAW currents.
To reveal the impact of changing sea level on surface velocity fields, absolute surface geostrophic velocities, derived from Absolute Dynamic Topography fields, are shown in Fig. 3a, b during the extreme negative (i.e., 1993) and positive (i.e., 2015) phases of the EOF. These two contrasting years illustrate the ongoing strengthening of the ocean surface circulation in the EAC associated with the Atlantic inflow. Surface geostrophic current linear trends showed that surface currents changed in most of the Atlantic pathway (Fig. 3c). In the Norwegian Sea, the increase in the Atlantic surface current speed reached +2 mm s−1 yr−1, corresponding to a relative increase ranging from 30% on the shelf to 100% in the basin. The surface currents in the eastern BS Atlantic corridor also significantly increased by about +0.6 mm s−1 yr−1 (+14%). The positive EOF phase associated with stronger north-eastward NAW surface currents since year 2000 reflects an increase in advection at or near the surface. This result is consistent with the recent trends in NAW current velocities observed upstream in the Nordic Seas40 or modelled in the BS41.
The interannual fluctuations as well as the long-term trend observed in the time-series of surface geostrophic velocities are mainly attributed to the dynamics of the North Atlantic Subpolar Gyre (Pearson’s correlation coefficient r = 0.58; p = 0.0028 two-sided t-test), and to the atmospheric forcing of the North Atlantic Oscillation (Pearson’s correlation coefficient r = −0.60, p = 0.003 two-sided t-test), which are both tightly coupled42. This suggests that the increase in surface advection in the EAC is likely due to a natural multidecadal oscillation related to the upstream synoptic oceanic circulation and atmospheric forcing, which could, in turn, drive long-term climatic change in the BS43.
Increasing NAW surface currents control the spatial distribution of Ehux
The first hypothesis that ocean currents control the summer spatial distribution of Ehux in the BS was tested using a Lagrangian experiment (EXP1). We advected virtual particles (here considered as the inoculum of Ehux cells), using observations of surface geostrophic velocities, from their expected overwintering location in March (defined by sea surface temperature >4 °C and distance from coast ≤180 km: see material and methods and Supplementary Note 2.5 and 2.6 for more details on the areal distribution and timing of the inoculum). EXP1 revealed that the spatial distribution of the virtual particles at the end of the advection period matched the extent of Ehux blooms (evidenced by satellite-derived Particulate Inorganic Concentration, PIC), with 80% of tracked particles ending up within 50 km of an Ehux blooming location (see Supplementary Note 2.7). The spatial distribution of particles in 1998 (beginning of the satellite ocean-color time-series) and 2015 are shown in Fig. 4, while the entire 1998–2016 period is displayed in Supplementary Fig. 14. The year-to-year robustness of the matchup between virtual particles and PIC clearly supports the fact that the NAW surface currents shape the location and extent of Ehux blooms in the EAC.
Increasing NAW surface currents control the poleward expansion of Ehux
It is noteworthy that the particles reached further north and east in the EAC in 2015 than in 1998, in agreement with Ehux blooms (Fig. 4). The north-eastward expansion of the Ehux bloom, expressed here as the distance reached by the leading-edge (see the Supplementary Method 1.1), increased on average by 325 km during the last 19 years as indicated by ocean-color observations (Fig. 5a), in close agreement with previous estimations (324 km, 40–50°E, 1989–201616;) and with the estimates from EXP1 (424 km, Fig. 5b).
The 4 °C surface isotherm (Fig. 4) is considered for Ehux as the lowest temperature required for sufficient growth to allow bloom formation (Supplementary Note 2.5). However, the Ehux blooms do not seem to follow this isotherm in summer and appear to be constrained by another factor. For example, the north-eastward expansion of the winter 4 °C isotherm, which delimits the areal extent of the inoculum in our Lagrangian method, could contribute to reduce the distance between the Ehux winter location and the Arctic domain. To test this hypothesis, we performed two additional Lagrangian experiments to examine both the role of currents and winter temperature on the poleward expansion of Ehux. In EXP2 we constantly initialized the virtual particles at the same inoculation region for all years, using a climatological mean temperature field to determine the inoculation region. Hence, EXP2 exclusively reflected the role of currents on the interannual variability of the advected particles. In contrast, in EXP3, the constant inoculation region varied from one year to another, but we used a climatological mean current to advect the particles. In this way, the interannual variability of the advected particles due to interannual variability of the inoculation region was quantified. The interannual PIC leading-edge location was highly correlated with the leading-edge from EXP2 (Pearson’s correlation coefficient r = 0.71, p < 0.01 two-sided t-test), suggesting a stronger control of the bloom expansion by currents than by winter temperature (EXP3, Pearson’s correlation coefficient r = 0.43, p = 0.06 two-sided t-test). On the decadal scale, currents (EXP2) were found responsible for the 56% (240 km) increase in the long-term leading-edge expansion against 44% (186 km) for winter temperature (EXP3), when compared to the reference EXP1. This significant increasing trend in current velocities (see Supplementary Note 2.8) was also revealed by the greater distance covered by the virtual particles reaching the BS, which increased by 110 km on average since 1993 (EXP2, Supplementary Note 2.9).
Increasing water temperature has previously been assumed to be the main driver of the spatial distribution of Ehux bloom in the BS15,16,44. Our results demonstrate that the primary driver of the Ehux dynamics (i.e., spatial distribution and timing) could be, in fact, stronger surface currents, which in turn intrinsically shape the temperature field and frontal structures. We show that oceanic currents (i.e. their intensity and fluctuations) drive the spatial distribution of the bloom, its interannual variability and more than 50% of the long-term poleward expansion of Ehux bloom in the BS. More importantly, from 2006 and onward, the contribution of water temperature to the expansion of Ehux bloom becomes negligible (Supplementary Table 3), and its poleward expansion is entirely due to the accelerating currents.
A newly “Atlantified” and retreating marine Arctic ecosystem
Ehux is largely studied given its significant role in marine geochemical cycles45,46 as illustrated by its sensitivity to ocean acidification47, its effect on pCO2 and CO2 uptake48 and its role on carbon export by providing calcite ballast effect49. Here, Ehux was used as an indicator of Atlantic ecosystems32,50. By expanding poleward and doubling its areal extent in the BS (Supplementary Note 2.10), Ehux attests to the ongoing “Atlantification” of the Arctic Ocean. Both arctic and Atlantic domains have distinct ecological signatures51, and the latter is undeniably “invading” the former4,52. Advected with the surface currents, Ehux will have to survive during the travel (e.g., avoid grazing53,54, subduction under the polar mixed layer in Fram Strait55) until finding more favourable blooming conditions in the BS in summer13. The fate of Ehux in the BS is therefore of major importance as it determines the potential “seeding effect” of Ehux in the Arctic regions. Ehux seems to be adapted to the low light, low nutrient, oligotrophic and highly stratified conditions of the NAW in summer30 such that its expansion, growth and blooming in the Arctic Ocean will be limited at some point by those constraints (bottom-up) but also by the grazing pressure (top-down).
Despite its adaptation to low light conditions, Ehux still requires sufficient light levels to sustain the energy-demanding calcification of its coccoliths. Such conditions are met in highly stratified oceans where Ehux can accumulate in the surface layer56. At high latitudes (>81°N), even with sea surface temperature above 3 or 4 °C, the survival of Ehux would require adaptation to rapidly decreasing solar radiations in late summer57. Ehux fate thus mainly relies on its ability to drift, with the appropriate timing, to highly stratified and temperate areas that allow it to stay in the surface euphotic layer. These conditions58 would likely be met in the Eurasian interior shelves of the Arctic Ocean28 (i.e., Kara, Laptev, and Siberian seas) where surface waters are warming and freshening59. If the increase in advection along the shelf slope continues in the future, we can expect Ehux to become a summer resident of the newly “Atlantified” Eurasian interior shelves, as previously revealed during the last interglacial50,60.
By driving such a poleward expansion, advective processes could affect the entire marine ecosystems by shifting species distribution61 and modifying interactions at higher trophic levels62. The concomitant decline of silicate concentrations in NAW63 may also contribute to the success of non-silicifying and small phytoplankton such as Ehux64. In addition, a change toward temperate pelagic species could have an impact on energy transfer to higher trophic levels65 including marine mammals and commercial fish stocks66. Considering the role of “bio-advection” in ecological models (i.e., trait-based and niche-based approaches) must improve predictions of community shifts67. The comparable increase in poleward advection of Pacific waters occurring in the Bering Strait68 suggests that the shrinking polar domain of the Arctic Ocean may be prone to intrusions of temperate species at a pan-Arctic scale.
Methods
Satellite altimetry
The satellite altimetry product was generated by the Sea Level Thematic Assembly Center (SL-TAC) and distributed by the Copernicus Marine and Environment Monitoring Service (CMEMS, http://marine.copernicus.eu/). This is a daily gridded global product in delayed-time (re-analysis product ID: SEALEVEL_GLO_PHY_L4_REP_OBSERVATIONS_008_047) with a spatial resolution of 1/4° available from 1993 to present time. A 1/4° Mercator projection grid roughly corresponds to a resolution of 27.8 km in latitude × 9.5 km in longitude at 70°N. This product is generated using Sentinel 3 A, Jason 3, HY2A, Saral/AltiKa, Cryosat 2, OSTM/Jason 2, Jason-1, Topex/Poseidon, Envisat, GFO, ERS-1/2. The SL-TAC produced this dataset with a special care for synchronization, homogenization, cross-calibration, correction of large-scale biases, and cross-validation. Before merging missions by optimal interpolation, the SL-TAC filtered the dataset to remove residual noise. The reader is referred to the product documentation for more details. The satellite altimetry product provides Sea Level Anomaly that estimates the variation of the Sea Surface height around a vertical reference (Mean Sea Surface). The product also provides Absolute Dynamic Topography (=Sea Level Anomaly+Mean Dynamic Topography), geostrophic velocity anomalies and absolute geostrophic velocities (derived respectively from Sea Level Anomaly and Absolute Dynamic Topography). The Mean Dynamic Topography employed by the SL-TAC was the MDT-CNES-CLS13 produced by the Collecte Localisation Satellites Space Oceanography division and distributed separately by Aviso (http://www.aviso.altimetry.fr).
SLA has been successfully used over the last 25 years to study the variability of sea level and surface geostrophic circulation at different spatial and temporal scales. However, most of the validation studies focused on the areas below the polar circle (<66°N) such as in the Norwegian Sea where derived geostrophic currents were successfully evaluated against in situ current meters69. The use of altimetry in Arctic regions has been avoided in the past due to (1) persistent sea-ice cover and (2) limitation of some radar altimeters (TOPEX/POSEIDON, Jason-1, OSTM/Jason-2 and Jason-3 missions), which do not acquire measurements above the polar circle. First, sea ice was not an issue for this study because virtual particles were only evolving in the year-round ice-free Atlantic domain. For the trend and EOF calculation, data associated with sea ice were systematically precluded (see below). Second, measurements above 66°N can rely on several other satellites: ERS-1 and -2, Envisat, SARAL/AltiKa, GFO (<72°N), HY-2A, ICESat, Cryosat-2, Sentinel-3. A comprehensive validation of the gridded product of geostrophic velocities was recently provided for ERS-1/2 and Envisat satellites using surface drifters and tide gauges in the whole EAC, with a particular focus on the Barents Sea70,71. The authors demonstrated that RMS differences between the drifter (corrected for Ekman currents) and altimetry velocities ranged within 7–15 cm s−1, which was comparable to previous estimates at lower latitudes. More importantly, they also concluded that the dataset provided sufficient spatial and temporal coverage to consistently resolve the synoptic and large mesoscale variability. Moreover, the product was recently found appropriate to monitor the long-term northward Atlantic inflow in the Norwegian Sea40 with even lower RMS differences between drifters and satellites velocities (0.7–8 cm s−1).
The dataset used in this study was enhanced from additional satellites (i.e., CryoSat-2 and Sentinel-3A) that offered a better data coverage at high latitudes for recent years (since 2011). Surface currents inferred from altimetry were found to reflect the variability of the deeper layers in the slope current in the Nordic Seas40. For this study, the derived surface geostrophic velocities were evaluated again toward five long-term in situ current meters in the southwestern Barents Sea and showed reasonably good Pearson’s correlation coefficients (r = 0.43–0.71, Supplementary Note 2.1). In addition, by evaluating the effect of steric and mass related effect on the SLA, we demonstrated that the increasing long-term trend was driven by change in currents (Supplementary Note 2.2).
Satellite sea surface temperature (SST)
We used the NOAA 1/4° degree and daily Optimum Interpolation Sea Surface Temperature version 2 (OISSTv2, provided by the NOAA/OAR/ESRL PSD, Boulder, Colorado, USA, https://www.esrl.noaa.gov/psd/data/gridded/data.noaa.oisst.v2.highres.html72) to extract the 4 °C isotherms during the month of March between 1998 and 2016. This dataset was constructed by combining satellites observations from the Advanced Very High-Resolution Radiometer (AVHRR) and the Advanced Microwave Scanning Radiometer on the Earth Observing System (AMSR-E) with in situ platforms (ships, buoys) on a regular global spatial grid. Missing observations in SST were optimally interpolated. The dataset was delivered with error fields providing a measure of confidence. This allowed us to build monthly composites maps excluding values with estimated errors higher than 50%, following the same method used for the altimetry dataset. Note that the dataset also provided daily sea-ice concentration maps, which were used as masks for the altimetry dataset. Additional information can be found here: https://www.ncdc.noaa.gov/oisst.
Satellite particulate inorganic carbon (PIC)
Ocean colour satellites detect coccolithophore blooms because they produce calcite plates (coccoliths) that shed light and produce ‘turquoise’ colored waters. Remotely sensed ocean-color offers a 20-year time-series from 1998 to present. In this study, we used the GlobColour database (http://hermes.acri.fr), which provides a continuous dataset of ocean-color satellite data products derived from the merging of five ocean-color sensors MEdium Resolution Imaging Spectrometer (MERIS), Moderate Resolution Imaging Spectroradiometer (MODIS), Sea-viewing Wide Field-of-view Sensor (Sea-WiFS), Visible Infrared Imaging Radiometer Suite (VIIRS) and Ocean and Land Colour Imager (OLCI-A). These sensors measure visible and near-infrared solar radiation reflected back from the ocean surface layer. Such remotely sensed information is only available during the daytime and in the absence of ice and clouds.
Level 3 PIC products were derived from NASA’s standard PIC algorithm73,74. This algorithm is based on a robust relationship between the light backscattering coefficient and the concentration of coccoliths, calcite plates forming the coccosphere of the coccolithophores30,75. Some coccolithophore species, including Ehux, produce and release coccoliths into the water during the later stages of a bloom76,77,78, creating large patches of highly reflective waters, which can be observed from space56. Satellite-derived PIC concentration provides a good proxy for Ehux coccolith concentration because Ehux dominates the coccolithophore population in the Norwegian Sea79 and in the Barents Sea80. Moreover, satellite PIC products have been successfully validated during Ehux blooms in the Barents Sea13,80,81, as well as in the global ocean56,74,82.
The leading-edge detection
The leading-edge was the main metric used in this study to follow the poleward (i.e., northward and eastward) expansion of the Ehux bloom. The leading-edge was derived as the 95th percentile of the longitude and latitude positions of the virtual particles entering the eastern Barents Sea. See Supplementary Method 1.1 for a more detailed description.
Empirical Orthogonal Function analysis of sea level and trends in surface geostrophic velocities
Empirical Orthogonal Function (EOF) and trends analysis were performed between 1993 and 2016, with monthly composite maps of non-seasonal Sea Level Anomalies (SLA) and surface absolute geostrophic velocities. Monthly composite maps were derived excluding pixels with error estimates (percentage of the SLA signal variance) higher than 50% or sea-ice concentrations higher than 15%. The seasonal cycle, known to be highly variable for SLA and the velocity fields, was removed because the analysis focused on the interannual and decadal time scale. The maximum amplitude in SLA generally occurs in fall83, whereas the velocity fields are maximum during the winter84. More details on the EOF analysis are available in Supplementary Note 2.3. Trends in surface absolute geostrophic velocities were only derived for pixels associated with more than 50% of data points.
Lagrangian calculation of particles trajectories
The Lagrangian experiments needed daily temporal resolution and directly employed daily absolute surface geostrophic velocities delivered from CMEMS. Three different Lagrangian simulations were used to address specific objectives. In each simulation, particles were advected with a 4th order Runge-Kutta integration scheme, with a time step of 3h. The simulations showed good robustness to the addition of noise, which was representative of the velocity fields uncertainties (Supplementary Note 2.4). The different simulations are described below.
-
(1)
Variable temperature and currents. The first simulation (EXP1) explored the influence of seasonal current variability on Ehux distribution from 1998 to 2016. On the 1st of March of each year, we selected a patch of particles representing the theoretical inoculum area of Ehux (i.e., SST ≥ 4 °C and distance from coast < 180 km, see Supplementary Note 2.5). These particles were advected, using satellite-derived velocities, over the entire growing season, when sufficient light availability is expected to allow phytoplankton growth (i.e., from the March 1st to September 1st). The choice of March 1st as the starting date is discussed in Supplementary Note 2.6.
-
(2)
Constant temperature and variable currents. The second simulation (EXP2) investigated the long-term increase in particle advection in the BS (1993–2016). In this simulation, instead of having a distinct ‘inoculum patch’ for each year, we inoculated the same patch every year. This patch was built using the March SST climatology (1993–2016) and the same criterion than EXP1 (SST ≥ 4 °C and distance from coast < 180 km). The particles were then also advected during the entire growing season (i.e., from March 1st to September 1st) as in EXP.1.
-
(3)
Variable temperature and constant currents. The third simulation (EXP3) aimed at assessing the role of temperature in the long-term poleward expansion of Ehux. EXP3 is similar to EXP1 with the exception that the particles are advected each year on the same climatological velocity field (varying ‘inoculum patch’).
Statistical analysis
Trends were estimated using linear least-squares regression. Correlations were calculated using Pearson’s coefficient correlation r. Two-sided t-tests with N–2 degrees of freedom were used to find significant levels (p-values) for trends and correlations. The significance level for this study was considered as 95% (p < 0.05).
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Data availability
All the data used in this research are freely available to the public and may be downloaded through the links detailed in the Methods section. Analysis outputs data for the main figures are made available at https://gitlab.com/loloziel/oziel_et_al_nc_2020.
Code availability
The code for Lagrangian experiments is based on Lamta 0.2 Copyright (C) 2007/12/5 Francesco d’Ovidio (francesco.dovidio@locean-ipsl.upmc.fr). This is a free Octave software under the terms of the GNU General Public License (http://www.gnu.org/licenses/) and can be distributed and modified. Outputs were analyzed using ©Matlab custom scripts that are accessible at https://gitlab.com/loloziel/oziel_et_al_nc_2020.
References
Wassmann, P., Slagstad, D. & Ellingsen, I. Primary production and climatic variability in the European sector of the Arctic Ocean prior to 2007: preliminary results. Polar Biol. 33, 1641–1650 (2010).
Huang, J. et al. Recently amplified arctic warming has contributed to a continual global warming trend. Nat. Clim. Chang 7, 875–879 (2017).
Stroeve, J. & Notz, D. Changing state of Arctic sea ice across all seasons. Environ. Res. Lett. 13, 103001 (2018).
Oziel, L., Sirven, J. & Gascard, J. C. The Barents Sea frontal zones and water masses variability (1980–2011). Ocean Sci. 12, 169–184 (2016).
Årthun, M., Eldevik, T., Smedsrud, L. H., Skagseth, Øystein & Ingvaldsen, R. B. Quantifying the influence of atlantic heat on barents sea ice variability and retreat. J. Clim. 25, 4736–4743 (2012).
Smedsrud, L. H. et al. The role of the Barents Sea in the Arctic Climate System. Rev. Geophys. 51, 415–449 (2013).
Lien, V. S., Schlichtholz, P., Skagseth, Ø. & Vikebø, F. B. Wind-driven Atlantic water flow as a direct mode for reduced Barents Sea ice cover. J. Clim. 30, 803–812 (2016).
Polyakov, I. V. et al. Greater role for Atlantic inflows on sea-ice loss in the Eurasian Basin of the Arctic Ocean. Science 356, 285–291 (2017).
Wassmann, P., Slagstad, D., Riser, C. W. & Reigstad, M. Modelling the ecosystem dynamics of the Barents Sea including the marginal ice zone. J. Mar. Syst. 59, 1–24 (2006).
Ardyna, M. et al. Parameterization of vertical chlorophyll a in the Arctic Ocean: impact of the subsurface chlorophyll maximum on regional, seasonal, and annual primary production estimates. Biogeosciences 10, 4383–4404 (2013).
Bélanger, S., Babin, M. & Tremblay, J.-É. Increasing cloudiness in Arctic damps the increase in phytoplankton primary production due to sea ice receding. Biogeosciences 10, 4087–4101 (2013).
Dalpadado, P. et al. Productivity in the Barents Sea - Response to Recent Climate Variability. PLoS ONE 9, e95273 (2014).
Hegseth, E. N. & Sundfjord, A. Intrusion and blooming of Atlantic phytoplankton species in the high Arctic. J. Mar. Syst. 74, 108–119 (2008).
Winter, A., Henderiks, J., Beaufort, L., Rickaby, R. E. M. & Brown, C. W. Poleward expansion of the coccolithophore Emiliania huxleyi. J. Plankton Res. 36, 316–325 (2014).
Oziel, L. et al. Role for Atlantic inflows and sea ice loss on shifting phytoplankton blooms in the Barents Sea. J. Geophys. Res. Ocean 122, 5121–5139 (2017).
Neukermans, G., Oziel, L. & Babin, M. Increased intrusion of warming Atlantic water leads to rapid expansion of temperate phytoplankton in the Arctic. Glob. Chang. Biol. 24, 2545–2553 (2018).
Fossheim, M. et al. Recent warming leads to a rapid borealization of fish communities in the Arctic. Nat. Clim. Chang. 5, 673–677 (2015).
Ingvaldsen, R. B. & Gjøsæter, H. Responses in spatial distribution of Barents Sea capelin to changes in stock size, ocean temperature and ice cover. Mar. Biol. Res. 9, 867–877 (2013).
Eriksen, E., Skjoldal, H. R., Gjøsæter, H. & Primicerio, R. Spatial and temporal changes in the Barents Sea pelagic compartment during the recent warming. Prog. Oceanogr. 151, 206–226 (2017).
Olszewska, A. et al. Interannual zooplankton variability in the main pathways of the Atlantic water flow into the Arctic Ocean (Fram Strait and Barents Sea branches). ICES J. Mar. Sci. 74, 1921–1936 (2017).
Edvardsen, A., Slagstad, D., Tande, K. S. & Jaccard, P. Assessing zooplankton advection in the Barents Sea using underway measurements and modelling. Fish. Oceanogr. 12, 61–74 (2003).
Hunt, G. L. et al. Advection in polar and sub-polar environments: Impacts on high latitude marine ecosystems. Prog. Oceanogr. 149, 40–81 (2016).
Wassmann, P. et al. The contiguous domains of Arctic Ocean advection: trails of life and death. Prog. Oceanogr. 139, 42–65 (2015).
Sundfjord, A. et al. Seasonal variation in transport of zooplankton Into the Arctic basin through the Atlantic gateway, Fram Strait. Front. Mar. Sci. 5, 1–22 (2018).
Wassmann, P., Slagstad, D., Ellingsen, I. & Ross, R. M. Advection of Mesozooplankton Into the Northern Svalbard Shelf Region. Front. Mar. Sci. 6, 1–10 (2019).
Slagstad, D., Wassmann, P. F. J. & Ellingsen, I. Physical constrains and productivity in the future Arctic Ocean. Front. Mar. Sci. 2, 1–23 (2015).
Popova, E. E., Yool, A., Aksenov, Y. & Coward, A. C. Role of advection in Arctic Ocean lower trophic dynamics: a modeling perspective. J. Geophys. Res. Ocean 118, 1571–1586 (2013).
Carmack, E. & Wassmann, P. Food webs and physical-biological coupling on pan-Arctic shelves: unifying concepts and comprehensive perspectives. Prog. Oceanogr. 71, 446–477 (2006).
Vernet, M. et al. Influence of phytoplankton advection on the productivity along the Atlantic Water Inflow to the Arctic Ocean. Front. Mar. Sci. 6, 583 (2019).
Balch, W. M. Re-evaluation of the physiological ecology of coccolithophores. Coccolithophores 165–190 (Springer Berlin Heidelberg, 2004).
Paasche, E. A review of the coccolithophorid Emiliania huxleyi (Prymnesiophyceae), with particular reference to growth, coccolith formation, and calcification-photosynthesis interactions. Phycologia 40, 503–529 (2001).
Gafar, N. A. & Schulz, K. G. A three-dimensional niche comparison of Emiliania huxleyi and Gephyrocapsa oceanica: Reconciling observations with projections. Biogeosciences 15, 3541–3560 (2018).
Bratbak, G., Wilson, W. & Heldal, M. Viral control of Emiliania huxleyi blooms? J. Mar. Syst. 9, 75–81 (1996).
Vardi, A. et al. Host-virus dynamics and subcellular controls of cell fate in a natural coccolithophore population. Proc. Natl Acad. Sci. USA 109, 19327–19332 (2012).
Highfield, A., Evans, C., Walne, A., Miller, P. I. & Schroeder, D. C. How many Coccolithovirus genotypes does it take to terminate an Emiliania huxleyi bloom? Virology 466–467, 138–145 (2014).
Berge, G. Discoloration of the sea due to coccolithus huxleyi “bloom”. Sarsia 6, 27–40 (1962).
Tyrrell, T. & Merico, A. in Coccolithophores (eds. Thierstein, H. R. et al.) 75–97 (Springer, 2004).
Birkenes, E. Phytoplankton in the Oslo Fjord during a ‘Coccolithus huxleyi-summer’. Avh. Nor. Vidensk. Akad. Oslo I. Mat. Nat. Kl. 2, 1–23 (1952).
Garcia, V. M. T. et al. Environmental factors controlling the phytoplankton blooms at the Patagonia shelf-break in spring. Deep. Res. Part I Oceanogr. Res. Pap. 55, 1150–1166 (2008).
Raj, R. P. et al. Quantifying Atlantic Water transport to the Nordic Seas by remote sensing. Remote Sens. Environ. 216, 758–769 (2018).
Lien, V. S., Vikebø, F. B. & Skagseth, Ø. One mechanism contributing to co-variability of the Atlantic inflow branches to the Arctic. Nat. Commun. 4, 1488 (2013).
Hátún, H., Sande, A. B., Drange, H., Hansen, B. & Valdimarsson, H. Ocean science: influence of the atlantic subpolar gyre on the thermohaline circulation. Science 309, 1841–1844 (2005).
Zhang, L. et al. The North Atlantic Oscillation as a driver of rapid climate change in the Northern Hemisphere. Nat. Geosci. 9, 509–512 (2016).
Signorini, S. R. & McClain, C. R. Environmental factors controlling the Barents Sea spring-summer phytoplankton blooms. Geophys. Res. Lett. 36, 1–5 (2009).
Riebesell, U. et al. Competitive fitness of a predominant pelagic calcifier impaired by ocean acidification. Nat. Geosci. 10, 19–23 (2017).
Taylor, A. R., Brownlee, C. & Wheeler, G. Coccolithophore cell biology: chalking up progress. Ann. Rev. Mar. Sci. 9, 283–310 (2016).
Beaufort, L. et al. Sensitivity of coccolithophores to carbonate chemistry and ocean acidification. Nature 476, 80–83 (2011).
Delille, B. et al. Response of primary production and calcification to changes of pCO2 during experimental blooms of the coccolithophorid Emiliania huxleyi. Glob. Biogeochem. Cycles 19, 1–14 (2005).
Klaas, C. & Archer, D. E. Association of sinking organic matter with various types of mineral ballast in the deep sea: Implications for the rain ratio. Glob. Biogeochem. Cycles 16, 63-1–63–14 (2002).
Baumann, M. in Geological History of the Polar Oceans: Arctic versus Antarcti pp. 437–445 (Springer, Netherlands, 1990).
Dalpadado, P. et al. Climate effects on Barents Sea ecosystem dynamics. ICES J. Mar. Sci. 69, 1303–1316 (2012).
Lind, S., Ingvaldsen, R. & Furevik, T. Arctic warming hotspot in the northern Barents Sea linked to declining sea-ice import. Nat. Clim. Change 8, 634–639 (2018).
Valiela, I. Coastal altimetry. Marine Ecological Processes 189-416 (Springer, 1995).
Olli, K. et al. The fate of production in the central Arctic Ocean—top-down regulation by zooplankton expatriates? Prog. Oceanogr. 72, 84–113 (2007).
Rudels, B., Jones, E. P., Schauer, U. & Eriksson, P. Atlantic sources of the Arctic Ocean halocline. Polar Res. 23, 10767 (2003).
Holligan, P. M. et al. A biogeochemical study of the coccolithophore, Emiliania huxleyi, in the North Atlantic. Glob. Biogeochem. Cycles 7, 879–900 (1993).
Kaartvedt, S. Photoperiod may constrain the effect of global warming in arctic marine systems. J. Plankton Res. 30, 1203–1206 (2008).
Krumhardt, K. M., Lovenduski, N. S., Iglesias-Rodriguez, M. D. & Kleypas, J. A. Coccolithophore growth and calcification in a changing ocean. Prog. Oceanogr. 159, 276–295 (2017).
Williams, W. J. & Carmack, E. C. The ‘interior’ shelves of the Arctic Ocean: physical oceanographic setting, climatology and effects of sea-ice retreat on cross-shelf exchange. Prog. Oceanogr. 139, 24–41 (2015).
Backman, J., Fornaciari, E. & Rio, D. Biochronology and paleoceanography of late Pleistocene and Holocene calcareous nannofossil abundances across the Arctic Basin. Mar. Micropaleontol. 72, 86–98 (2009).
Kortsch, S., Primicerio, R., Fossheim, M., Dolgov, A. V. & Aschan, M. Climate change alters the structure of arctic marine food webs due to poleward shifts of boreal generalists. Proc. R. Soc. B Biol. Sci. 282, 20151546 (2015).
Frainer, A. et al. Climate-driven changes in functional biogeography of Arctic marine fish communities. Proc. Natl Acad. Sci. USA 114, 12202–12207 (2017).
Rey, F. Declining silicate concentrations in the Norwegian and Barents Seas. ICES J. Mar. Sci. 69, 208–212 (2012).
Hátún, H. et al. The subpolar gyre regulates silicate concentrations in the North Atlantic. Sci. Rep. 7, 1–9 (2017).
Renaud, P. E. et al. Pelagic food-webs in a changing Arctic: a trait-based perspective suggests a mode of resilience. ICES J. Mar. Sci. 75, 1871–1881 (2018).
Bogstad, B., Gjøsæter, H., Haug, T. & Lindstrøm, U. A review of the battle for food in the Barents Sea: Cod vs. marine mammals. Front. Ecol. Evol. 3, 29 (2015).
Beaugrand, G. et al. Prediction of unprecedented biological shifts in the global ocean. Nat. Clim. Chang 9, 237–243 (2019).
Woodgate, R. A. Increases in the Pacific inflow to the Arctic from 1990 to 2015, and insights into seasonal trends and driving mechanisms from year-round Bering Strait mooring data. Prog. Oceanogr. 160, 124–154 (2018).
Mork, K. A. & Skagseth, Øystein. A quantitative description of the Norwegian Atlantic current by combining altimetry and hydrography. Ocean Sci. 6, 901–911 (2010).
Vignudelli, S., Kostianoy, A. G., Cipollini, P. & Benveniste, J. Coastal altimetry. 389–413 (Springer 2011).
Volkov, D. L. & Pujol, M.-I. I. Quality assessment of a satellite altimetry data product in the Nordic, Barents, and Kara seas. J. Geophys. Res. Ocean 117, 1–18 (2012).
Reynolds, R. W. et al. Daily high-resolution-blended analyses for sea surface temperature. J. Clim. 20, 5473–5496 (2007).
Gordon, H. R. et al. Retrieval of coccolithophore calcite concentration from SeaWiFS imagery. Geophys. Res. Lett. 28, 1587–1590 (2001).
Balch, W. M., Gordon, H. R., Bowler, B. C., Drapeau, D. T. & Booth, E. S. Calcium carbonate measurements in the surface global ocean based on Moderate-Resolution Imaging Spectroradiometer data. J. Geophys. Res. Ocean 110, 1–21 (2005).
Paasche, E. Roles of nitrogen and phosphorus in coccolith formation in Emiliania Huxleyi (prymnesiophyceae). Eur. J. Phycol. 33, 33–42 (1998).
Borman, A. H., De Jong, E. W., Huizinga, M. & Westbroek, P. in Biomineralization and Biological Metal Accumulation 303–305 (Springer, Netherlands, 1983).
Westbroek, P., Young, J. R. & Linshooten, K. Coccolith production (Biomineralization) in the Marine Alga Emiliania huxleyi. J. Protozool. 36, 368–373 (1989).
Feng, Y. et al. Interactive effects of increased pCO2, temperature and irradiance on the marine coccolithophore Emiliania huxleyi (Prymnesiophyceae). Eur. J. Phycol. 43, 87–98 (2008).
Dylmer, C. V., Giraudeau, J., Hanquiez, V. & Husum, K. The coccolithophores Emiliania huxleyi and Coccolithus pelagicus: Extant populations from the Norwegian-Iceland Seas and Fram Strait. Deep. Res. Part I Oceanogr. Res. Pap. 98, 1–9 (2015).
Giraudeau, J. et al. A survey of the summer coccolithophore community in the western Barents Sea. J. Mar. Syst. 158, 93–105 (2016).
Hovland, E. K. et al. Optical impact of an Emiliania huxleyi bloom in the frontal region of the Barents Sea. J. Mar. Syst. 130, 228–240 (2014).
Hopkins, J., Henson, S. A., Painter, S. C., Tyrrell, T. & Poulton, A. J. Phenological characteristics of global coccolithophore blooms. Glob. Biogeochem. Cycles 29, 239–253 (2015).
Volkov, D. L., Landerer, F. W. & Kirillov, S. A. The genesis of sea level variability in the Barents Sea. Cont. Shelf Res. 66, 92–104 (2013).
Ingvaldsen, R. B., Asplin, L. & Loeng, H. Velocity field of the western entrance to the Barents Sea. J. Geophys. Res. Ocean 109, 1–12 (2004).
Acknowledgements
L.O. was funded by an NSERC/CRSNG “Visiting Fellowship” in Canadian laboratories in the framework of the DFO-International Governance Strategy (IGS). L.O. was also supported by the joint international laboratory Takuvik (ULaval/CNRS), the Canada Excellence Research Chair in remote sensing of Canada’s new Arctic frontier. This research was supported by the Sentinel North program of Université Laval, made possible, in part, thanks to funding from the Canada First Research Excellence Fund. M.A. was supported by a European Union’s Horizon 2020 Marie Sklodowska-Curie grant (no. 746748). The authors would like to thank the international agencies and programs from the Copernicus Marine Environment Monitoring Service (CMEMS), the National Oceanic and Atmospheric Administration (NOAA), and the ACRI-ST company for freely providing access to their satellite altimetry, Sea Surface Temperature and Ocean-color data. This work heavily relies on the following toolbox for mapping: Pawlowicz, R., 2018. “M_Map: A mapping package for MATLAB”, version 1.4j, [Computer software], available online at www.eoas.ubc.ca/~rich/map.html as well as environmental toolbox developed by Chad Greene (2019): The Climate Data Toolbox for Matlab (https://www.github.com/chadagreene/CDT), GitHub. Retrieved March 11, 2019. The authors thank Francesco D’Ovidio for the Lagrangian tools, Martine Lizotte for language support and Marie-Isabelle Pujol for her advices on satellite velocity uncertainties.
Author information
Authors and Affiliations
Contributions
L.O. conceptualized the study, conducted data analysis, and wrote the manuscript. A.B. performed the Lagrangian experiments. M.A. helped in the writing and in the conceptualization of the manuscript. P.M. and R.B.I. helped in the analysis of the steric contribution on SLA. R.I. provided current meters data and helped in the validation analysis. A.R. derived the virtual particles match-ups. All authors (L.O., A.B., M.A., P.M., A.R., J-B.S., R.B.I, E.D., and M.B.) contributed to the ideas and edited the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Luc Beaufort, Paul Wassmann, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Oziel, L., Baudena, A., Ardyna, M. et al. Faster Atlantic currents drive poleward expansion of temperate phytoplankton in the Arctic Ocean. Nat Commun 11, 1705 (2020). https://doi.org/10.1038/s41467-020-15485-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-15485-5
This article is cited by
-
Biogeographic gradients of picoplankton diversity indicate increasing dominance of prokaryotes in warmer Arctic fjords
Communications Biology (2024)
-
Algal blooms in the ocean: hot spots for chemically mediated microbial interactions
Nature Reviews Microbiology (2024)
-
The distribution and diversity of eukaryotic phytoplankton in the Icelandic marine environment
Scientific Reports (2023)
-
Atlantic water influx and sea-ice cover drive taxonomic and functional shifts in Arctic marine bacterial communities
The ISME Journal (2023)
-
Atlantic-origin water extension into the Pacific Arctic induced an anomalous biogeochemical event
Nature Communications (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.