Hostname: page-component-7b9c58cd5d-bslzr Total loading time: 0.001 Render date: 2025-03-15T12:08:08.204Z Has data issue: false hasContentIssue false

Palynological evidence of sea-surface conditions in the Barents Sea off northeast Svalbard during the postglacial period

Published online by Cambridge University Press:  08 April 2020

Camille Brice
Affiliation:
GEOTOP, Université du Québec à Montréal, CP 8888, Montréa, QC H3C3P8 Canada
Anne de Vernal*
Affiliation:
GEOTOP, Université du Québec à Montréal, CP 8888, Montréa, QC H3C3P8 Canada
Elena Ivanova
Affiliation:
Shirshov Institute of Oceanology, Russian Academy of Sciences, 117997 Moscow, Russia
Simon van Bellen
Affiliation:
GEOTOP, Université du Québec à Montréal, CP 8888, Montréa, QC H3C3P8 Canada
Nicolas Van Nieuwenhove
Affiliation:
GEOTOP, Université du Québec à Montréal, CP 8888, Montréa, QC H3C3P8 Canada Shirshov Institute of Oceanology, Russian Academy of Sciences, 117997 Moscow, Russia University of New Brunswick, 2 Bailey Drive, Fredericton, NB E3B5A3, Canada
*
*Corresponding author email : devernal.anne@uqam.ca
Rights & Permissions [Opens in a new window]

Abstract

Postglacial changes in sea-surface conditions, including sea-ice cover, summer temperature, salinity, and productivity were reconstructed from the analyses of dinocyst assemblages in core S2528 collected in the northwestern Barents Sea. The results show glaciomarine-type conditions until about 11,300 ± 300 cal yr BP and limited influence of Atlantic water at the surface into the Barents Sea possibly due to the proximity of the Svalbard-Barents Sea ice sheet. This was followed by a transitional period generally characterized by cold conditions with dense sea-ice cover and low-salinity pulses likely related to episodic freshwater or meltwater discharge, which lasted until 8700 ± 700 cal yr BP. The onset of “interglacial” conditions in surface waters was marked by a major change in dinocyst assemblages, from dominant heterotrophic to dominant phototrophic taxa. Until 4100 ± 150 cal yr BP, however, sea-surface conditions remained cold, while sea-surface salinity and sea-ice cover recorded large amplitude variations. By ~4000 cal yr BP optimum sea-surface temperature of up to 4°C in summer and maximum salinity of ~34 psu suggest enhanced influence of Atlantic water, and productivity reached up to 150 gC/m2/yr. After 2200 ± 1300 cal yr BP, a distinct cooling trend accompanied by sea-ice spreading characterized surface waters. Hence, during the Holocene, with exception of an interval spanning about 4000 to 2000 cal yr BP, the northern Barents Sea experienced harsh environments, relatively low productivity, and unstable conditions probably unsuitable for human settlements.

Type
Research Article
Copyright
Copyright © University of Washington. Published by Cambridge University Press, 2020

Introduction

Most of the circum-Arctic lands were marked by human occupation for millennia, even the Canadian Arctic and northern Greenland where early settlements of pre-Inuit date from about 4000 yr (e.g., Raghavan et al., Reference Raghavan, DeGiorgio, Albrechtsen, Moltke, Skoglund, Korneliussen and Grønnow2014). However, Svalbard apparently remained free of human occupation until historical times (Arlov, Reference Arlov2005). Artifacts dating from the Stone Age were tentatively attributed by Christiansson and Simonsen (Reference Christiansson and Simonsen1970) to early settlements, but the nature of the artifacts was questioned, and the evidence denied by some archaeologists (e.g., Bjerk, Reference Bjerk2000). Norsemen and Russian Pomors possibly made early stays in Svalbard by the end of the Middle Age but the first unquestioned settlements in Svalbard date from the late sixteenth century, when whalers and hunters from northwest European started to exploit the Svarbard coastal regions, at least seasonally (Albrethsen and Arlov, Reference Albrethsen and Arlov1988; Arlov, Reference Arlov2005). Whereas harsh climate and low marine productivity that would have sustained and/or attracted human populations could explain the late arrival of humans on Svalbard, little is known about climate and surface ocean conditions prior to the instrumental period, especially in the northeast, where the coast lines are bathed by the Barents Sea and Arctic Ocean waters. Here, we document postglacial sea-surface conditions off northeast Svalbard based on the palynological content of a core collected in Kvitøya Trough (Fig. 1) in the northern Barents Sea, which is the second main gateway between the Atlantic and Arctic oceans after Fram Strait (Årthun et al., Reference Årthun, Eldevik, Smedsrud, Skagseth and Ingvaldsen2012). Hence, the northern Barents Sea that makes the transition from Atlantic to Arctic waters is a region that might have experienced important climate changes.

Figure 1. Maps showing the location of core S2528 and modern surface currents in the Barents Sea. The red arrows illustrate the circulation pattern of Atlantic waters at the surface, including the Norwegian Atlantic Current (NwAC) and the Norwegian Atlantic Slope Current (NwASC), which divides into the West Spitsbergen Current (WSC) and the North Cape Current (NCC). North of Svalbard, the WSC becomes the subsurface current of the Svalbard Branch (SB) and the Yermak Branch (YB). The blue arrows show the Arctic East Spitsbergen Current (ESC) and East Greenland Current (EGC). The green arrows correspond to coastal currents. The polar front position is represented by the white line based on Loeng (Reference Loeng1991).

The inflow of Atlantic water northeastward through the Barents Sea Opening presently represents a mean flux of 2.3 Sv and constitutes an important source of heat in the Barents Sea, with an annual average of 49 TW as recorded between 1998 and 2008 (Årthun et al., Reference Årthun, Eldevik, Smedsrud, Skagseth and Ingvaldsen2012). Because Atlantic water loses most of its heat content in the Barents Sea (Årthun and Schrum, Reference Årthun and Schrum2010), remaining heat transport through the northern exit is relatively small (Gammelsrød et al., Reference Gammelsrød, Leikvin, Lien, Budgell, Loeng and Maslowski2009). In the Barents Sea, the polar front position at the winter limit of sea-ice cover extent is highly sensitive to variations in heat advected from the North Atlantic Ocean. The northwestern Barents Sea is therefore located in the center of action for what concerns exchanges between cold surface waters from the Arctic Ocean and the northernmost branches of the relatively warm North Atlantic water. In the Barents Sea, any change in Atlantic current intensity may impact the distribution of sea ice (Årthun et al., Reference Årthun, Eldevik, Smedsrud, Skagseth and Ingvaldsen2012) and ocean-atmosphere heat exchanges, thus playing an important role in regional climate and accounting for a part of the overall strength of the Atlantic Meridional Overturning circulation (AMOC; Rudels et al., Reference Rudels, Anderson and Jones1996; Dieckmann and Hellmer, Reference Dieckmann, Hellmer, Thomas and Dieckmann2008; Smedsrud et al., Reference Smedsrud, Esau, Ingvaldsen, Eldevik, Haugan, Li and Lien2013).

With the aim of reconstructing sea-surface conditions since the deglaciation and their relationship with variations in the relative flow of North Atlantic water into the Barents Sea, we analyzed the palynological content of a marine sediment core (S2528) collected in Kvitøya Trough in the northwestern part of the Barents Sea (Fig. 1). Among palynological tracers, assemblages of dinoflagellate cysts (hereafter dinocysts) deserve special attention. Dinoflagellates are unicellular organisms that include heterotrophic, phototrophic, and mixotrophic taxa (e.g., Jeong et al., Reference Jeong, Yoo, Kim, Seong, Kang and Kim2010), thus allowing to characterize productivity in the upper part of the water column. During their life cycle, in relation with their sexual reproduction, about 10–12% of the species form a cyst composed of refractory organic matter that fossilizes in sediment (e.g., Head, Reference Head, Jansonius and McGregor1996). The dinocyst fluxes to the sea floor yield a fragmentary picture of biological productivity in surface waters and their assemblage composition provide information on oceanographic conditions, including sea-surface temperature and salinity, seasonal sea-ice cover, and productivity (Rochon et al., Reference Rochon, de Vernal, Turon, Matthiessen and Head1999; de Vernal et al., Reference de Vernal, Henry, Matthiessen, Mudie, Rochon, Boessenkool and Eynaud2001, Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a, Reference de Vernal, Hillaire-Marcel, Rochon, Fréchette, Henry, Solignac and Bonnet2013b, Reference de Vernal, Radi, Zaragosi, Van Nieuwenhove, Rochon, Allan and De Schepper2019; Radi and de Vernal, Reference Radi and de Vernal2008). In this paper, we used dinocyst data to estimate conditions steered by the northward flow of Atlantic waters at the surface. Although the study core site is located offshore, far from potential human settlements, it provides insights into climate and marine environmental conditions for the northeast Svalbard region over the last millennia.

GEOGRAPHIC AND OCEANOGRAPHIC SETTING

The Barents Sea is a relatively shallow, continental shelf sea with a mean water depth of 230 m. Surface waters in the Barents Sea consist of three main water masses: the Coastal Water, the North Atlantic Water (NAW), and the Arctic Water (Fig. 2; Loeng Reference Loeng1991). The NAW is characterized by salinities higher than 35.0 psu. Between Norway and Bear Island, where the NAW remains close to the surface, sea-surface temperature (SST) ranges from 3.5 to 6.5°C in summer (Fig. 2). Both temperature and salinity decrease northeastward. The Coastal Water, formed from the mixing of North Atlantic water and terrestrial freshwater inputs (Sætre and Ljøen, Reference Sætre and Ljøen1972), has almost the same temperature as the NAW, but its salinity is lower than 34.7 psu. The Arctic Water originates from the north and is characterized by low salinity (34.3–34.8 psu) and SST below 0°C. During winter, in the northern part of the Barents Sea, the Arctic Water occupies the upper 150 m of the water column, while meltwater may result in a low-salinity (31–34.2 psu) surface-water layer 5–20 m thick in summer (Helland-Hansen and Nansen, Reference Helland-Hansen and Nansen1909). The NAW layer occupies most of the water column in the southern Barents Sea (Risebrobakken et al., Reference Risebrobakken, Moros, Ivanova, Chistyakova and Rosenberg2010, Reference Risebrobakken, Dokken, Smedsrud, Andersson, Jansen, Moros and Ivanova2011), but it lies mostly at depths between 120 and 200 m in the northern Barents Sea (Duplessy et al., Reference Duplessy, Ivanova, Murdmaa, Paterne and Labeyrie2001).

Figure 2. Map of mean summer (a) sea-surface temperature and (b) salinity over the Barents Sea between 1955 and 2012, after the World Ocean Atlas (Locarinini et al., Reference Locarnini, Mishonov, Antonov, Boyer, Garcia, Baranova and Zweng2013; Zweng et al., Reference Zweng, Reagan, Antonov, Locarnini, Mishonov, Boyer and Garcia2013). The white dashed lines separate the North Atlantic waters (NAW), the Arctic waters (ArcW) and the Coastal waters (CW), as defined by Loeng (Reference Loeng1991).

The North Cape Current, a North Atlantic branch flowing along the Scandinavian coasts, transports most of the heat to the Barents Sea through the Bear Island Trough (Fig. 1). It controls the polar front position at the sea-ice limit (Loeng, Reference Loeng1991). The inflow of the North Cape Current mainly follows a counterclockwise trajectory before exiting the Barents Sea between Novaya Zemlya and Franz Josef Land (Schauer et al., Reference Schauer, Loeng, Rudels, Ozhigin and Dieck2002). Relatively warm (+0.5 to 1°C) and saline (34.7–35.1 psu) modified Atlantic water subducting in Fram Strait spreads as a subsurface current along the continental slope of the Nansen Basin and enters the Barents Sea from the north via the relatively deep (400–600 m) Kvitøya and Franz Victoria Troughs, and in the Kara Sea via the St. Anna Trough (Lind and Invaldsen, Reference Lind and Ingvaldsen2012).

The surface water at our study site is presently under the dominant influence of the Arctic Water with SST <0°C and salinity between 34.3 and 34.7 psu in summer (Fig. 2; Pfirman et al., Reference Pfirman, Bauch, Gammelsrød, Johannessen, Muench and Overland1994). The Arctic Water is carried into the Barents Sea with the East Spitsbergen Current in the Franz Victoria and Kvitøya troughs, from where it flows south and west along the eastern Svalbard coast. Conductivity, temperature and density data show that the Arctic Water occupies the upper 50 m of the water column in the coring site area, with fresher water at the top (Ivanova et al., Reference Ivanova, Murdmaa, de Vernal, Risebrobakken, Peyve, Brice, Seitkalieva and Pisarev2019). Below 50 m, the subsurface to intermediate layer of the water column is occupied by the Atlantic-derived water coming from the Svalbard Branch that flows into the Kvitøya Trough from the northwest (Pfirman et al., Reference Pfirman, Bauch, Gammelsrød, Johannessen, Muench and Overland1994).

METHODS

Sediment core and chronology

The 385-cm-long sediment core S2528 (80°40.80N; 29°36.70E) was retrieved at a depth of 428 m in the Kvitøya Trough, between Nordaustlandet and the Kvitøya islands during Cruise 25 of R/V Akademik Nikolaj Strakhov in 2007.

The chronology of S2528 was established from accelerator mass spectrometry (AMS) 14C dates of mixed benthic foraminifers and bivalve shells collected in sediment slices of 3 cm or more in order to obtain enough material (Table 1). The AMS 14C dates are reported using the Libby half-life of 5568 yr and calibrated using the Marine13 calibration curve (Reimer et al., Reference Reimer, Bard, Bayliss, Beck, Blackwell, Bronk Ramsey and Buck2013). An additional correction (ΔR) of 71 ± 21 yr was used to account for the air–sea reservoir difference in the Barents Sea (Mangerud et al., Reference Mangerud, Bondevik, Gulliksen, Hufthammer and Høisæter2006). To create a chronological framework, we used Bacon (Blaauw and Christen, Reference Blaauw and Christen2011) and Clam (Blaauw, Reference Blaauw2010) software, both run in R. Bacon relies on a Bayesian approach and applies a high number of iterations for creating the “best model” within a confidence interval. The Clam software is simpler than Bacon but provides the option to define the type of interpolation and presents the advantage of allowing the sample thickness to be included. Because of the sparse dated intervals and because thick samples had to be sieved to collect enough biogenic carbonate for dating, we first applied Clam for interpolation and we developed the final age model with Bacon for calculating the confidence intervals or uncertainties.

Table 1. 14C dates of core S2528 and corresponding calibrated ages (Ivanova et al. [2019] with additions). Laboratory codes: Poz, Poznań Radiocarbon Laboratory; OS, NOSAMS facility at the Woods Hole Oceanographic Institution; CAMS, Lawrence Livermore National Laboratory.

Palynological analysis

Sediments from core S2528 were subsampled over 3 cm at 3-cm intervals for palynological analysis. Each 3-cm-thick sample was processed according to a standardized protocol (de Vernal et al., Reference de Vernal, Bilodeau and Henry2010). In brief, 1.5 to 5.0 cm3 was weighted and sieved at 106 and 10 μm. The dried fraction >106 μm, mainly composed of detrital material, was weighed and used as a proxy for ice rafting deposition and for hand-picking of foraminifers and other biogenic remains. The fraction between 10 and 106 μm was used for palynological preparations, which consisted in chemical treatment with room temperature HCl (10%) and HF (49%) to dissolve carbonate and silica particles. The residue was concentrated on a 10-μm nylon mesh and mounted on microscope slides with glycerol gelatin. Prior to the chemical treatment, Lycopodium clavatum spore tablets were added to the samples to allow estimating palynomorph concentrations (e.g., Mertens et al., Reference Mertens, Verhoeven, Verleye, Louwye, Amorim, Ribeiro and Deaf2009).

The palynological analysis was conducted by identifying and counting common palynomorphs, including contemporaneous and reworked pollen grains, spores, and dinocysts, Halodinium and other acritarchs, and foraminifer organic linings, with an optical microscope in transmitted light at magnifications of 400 × to 1000 × . A minimum of 300 dinocysts were identified and counted in most sample using the nomenclature summarized in Rochon et al. (Reference Rochon, de Vernal, Turon, Matthiessen and Head1999) and Radi et al. (Reference Radi, Bonnet, Cormier, de Vernal, Durantou, Faubert, Head, Henry, Pospelova, Rochon and Van Nieuwenhove2013). Relative abundance of dinocyst taxa was expressed as percentages versus counted dinocysts. Pollen grains and spores were counted to quantify inputs from land vegetation, using identification keys (e.g., McAndrews et al., Reference McAndrews, Berti and Norris1973). Reworked palynomorphs were counted as indicators of past sediment erosion, subsequent transport, and resedimentation. The concentration of palynomorphs was expressed as the number of specimens per cm3 of wet sediment. Palynomorph concentrations and dinocyst percentage data are presented in Supplementary Tables 1a and 1b.

Special attention was paid to the dinocyst data. Multivariate analyses using CANOCO 5 software (ter Braak and Smilauer, Reference ter Braak and Smilauer2012) were performed to quantitatively describe the assemblage. We used the modern analogue technique (MAT; Guiot and de Vernal, Reference Guiot and de Vernal2007) to quantitatively reconstruct the sea-surface conditions. In paleoceanography, the analogue approach is commonly used for quantitative estimation of past SST (Guiot and de Vernal, Reference Guiot and de Vernal2007). In the case of dinocyst assemblages for which we have developed a very large reference database (de Vernal et al., Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a, Reference de Vernal, Radi, Zaragosi, Van Nieuwenhove, Rochon, Allan and De Schepper2020), it is the preferred technique because it relies on similarities between fossil and modern assemblages and does not assume linear or unimodal relationships for calibrations based on regression. Moreover, because dinocyst assemblages depend upon temperature in addition to salinity, sea-ice cover, and trophic level, the MAT has the advantage of allowing simultaneous reconstruction of many parameters. Furthermore, MAT is appropriate in the present case since the reference database used here (n = 1492) includes a large number of data points in the polar and subpolar domains with a wide range of hydrographical conditions regarding salinity, sea-ice cover and temperatures (de Vernal et al., Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a, Reference de Vernal, Hillaire-Marcel, Rochon, Fréchette, Henry, Solignac and Bonnet2013b).

The parameters we reconstruct here include SST and °C in winter and summer, sea-surface salinity (SSS, psu) in winter and summer, sea-ice cover (months/yr > 50%) and primary productivity (gC/m2/yr). These parameters show a wide range of combinations in the reference database and explain the distribution of dinocyst assemblages as it was shown from canonical correspondence analyses (Radi and de Vernal, Reference Radi and de Vernal2008; de Vernal et al., Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a). The estimates are made from the standardized reference database of the Northern Hemisphere that includes 1492 reference sites and 66 taxa (http://www.geotop.ca; de Vernal et al., Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a). Reconstructions were made from the five best modern analogue sites found by the calculation of the distance (inversely proportional to the similarity) after log-transformation of percentage data, in order to give more weight to accompanying taxa that often have narrower ecological affinities than dominant taxa (e.g., de Vernal et al., Reference de Vernal, Henry, Matthiessen, Mudie, Rochon, Boessenkool and Eynaud2001). In the case of the n = 1492 database, the threshold distance for good analogues is 1.2. Modern oceanographical values at the five selected modern analogue sites were used to define the maximum and minimum possible and the most probable value, which corresponds to the average weighted inversely to the distance.

The validation tests to evaluate the accuracy of the approach indicated that the root mean square errors of prediction are ± 1.4 months/yr for the sea-ice cover, ±1.1°C and ±1.6°C for winter and summer SSTs, ±2.6 psu for summer SSS, and ± 55 gC/m2 for productivity (de Vernal et al., Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a). MAT has been criticized because it possibly underestimates the error of prediction due to spatial autocorrelation (Telford, Reference Telford2006; Telford and Birks, Reference Telford and Birks2009). However, different tests performed to overcome the effect of spatial autocorrelation illustrate that the performance of MAT are comparable or better than those of transfer function techniques (Guiot and de Vernal, Reference Guiot and de Vernal2011; Hohmann et al., Reference Hohmann, Kucera and de Vernal2020).

In order to assess the reliability of estimates beyond the error of prediction, a reliability index can be calculated, as proposed by de Vernal et al. (Reference de Vernal, Hillaire-Marcel and Darby2005, Reference de Vernal, Hillaire-Marcel, Rochon, Fréchette, Henry, Solignac and Bonnet2013b). The index takes into consideration the (1) dinocyst sum, (2) similarity of the best modern analogues, and (3) number of analogues used for reconstructions, based on the cumulative score system as described hereafter. (l) Dinocyst sum > 300 = 4; 300–200 = 3; 200–100 = 2; 100–50 = 1. (2) Distance of the first modern analogue < 0.3 = 4; 0.3–0.6 = 3; 0.6–0.9 = 2; 0.9–1.2 = 1. (3) Distance of the fifth modern analogue < 0.3 = 4; 0.3–0.6 = 3; 0.6–0.9 = 2; 0.9–1.2 = 1. (4) Number of analogues below threshold of 1.2 : 5 = 4; 4 = 3; 3 = 2; 2-1= 1. The cumulative scores then permit to define a reliability index, ranging from excellent (A = 14–16), very good (B = 12–13), good (C = 10–11), and weak (D < 10).

The results from MAT, including the most probable, maximum and minimum estimates, the distance and location of the first analogue, and the index of reliability are presented in Supplementary Table 1c.

RESULTS

Chronology of the core

Ten radiocarbon dates were obtained from intervals containing enough biogenic carbonates for AMS measurements (>1 mg). Between 170 and 300 cm, seven dates were obtained allowing to create a relatively high-resolution chronology spanning 13,400 to 9800 14C yr BP (Fig. 3, Table 1). Below this interval, no biogenic material to date could be recovered. At the top of the core, dates were obtained between 10 and 16 cm, yielding relatively recent ages of <850 14C yr BP. Between 15 and 170 cm, biogenic carbonate material was very rare. Nevertheless, it was possible to pick enough foraminifers in the 67–76 cm interval for an AMS 14C measurement that yielded an age of about 6470 14C yr BP. Hence, the setting of a chronology for the Holocene is hampered by the rarity of dating control points. Initial calibration and age modeling using the Bacon software suggested high sedimentation rates in the 300–170 cm interval and lower accumulation rates above, with a step decrease at 70 cm, which corresponds to the dating point. However, in Bacon, through its Bayesian approach, the potential for shifts in sedimentation rates between dated levels is limited. Therefore, the depth of changing sedimentation rates is heavily influenced by the position of dated levels. Hence, whereas Bacon-based age-depth modeling is suitable in the case of relatively constant sedimentation rates or when there is a high density of dated levels, it is not necessarily appropriate here. There is no sedimentary evidence for a change in sedimentary regime at the dated level at about 70 cm depth and the shift in sedimentation rate could have occurred gradually or abruptly at any point between 170 and 70 cm. We therefore created an age-depth relationship with Clam by applying a smooth spline interpolation, which resulted in progressive change in sedimentation rates. From the curve based on Clam, we retained an artificial tie point that was added in the input file used with Bacon. Specifically, from the Clam model, we retained the age and confidence interval obtained at 100 cm, which corresponds to a major change in dinocyst assemblages (Fig. 4 and 5). This age was added as a “datum” in the Bacon input file. As we wanted the final chronology to include the thickness of dated samples, we did two runs with Bacon including this datum: one with the upper level of the dated sample, and one with the lower level. Besides the different “depth” values, the settings for these runs were identical. By applying the whole depth interval range for each sample, we were able to more appropriately quantify the overall uncertainty associated with the age-depth model. The single best age estimation was defined by the mean value of the two medians of both confidence intervals, while the 95% uncertainty interval was defined by the youngest and oldest extremes of the two combined 95% uncertainties. The resulting chronology (Fig. 3) is used here to assess the timing, with uncertainties, of the main transitions identified from the dinocyst assemblages.

Figure 3. Age versus depth relationship and sedimentation rates of core S2528 based on 14C dates and the use of an approach combining the Clam and Bacon softwares (Blaauw, Reference Blaauw2010; Blaauw and Christen, Reference Blaauw and Christen2011). Numbers next to polygons are the calibrated ages (yr BP) and the one in italics is the added age from interpolation with Clam (see section Chronology of the core).

Figure 4. Core lithological units, concentrations of dinocysts, pollen and spores, organic foraminifer linings, and reworked palynomorphs expressed in number of specimens/cm3 of sediment, and percentage of the coarse fraction (> 106 μm) in core S2528. The coarse fraction is expressed on a logarithmic scale due to large amplitude variations throughout the core.

Figure 5. Relative abundance (percentages) of the main dinocyst taxa in core S2528, concentration of phototrophic and heterotrophic dinocyst taxa, and results of the principal component (PC) analyses. The PC1 (black line) and PC2 (gray dashed line) axes respectively explain 45.2 and 24.5%, respectively, of the variance. The gray bands correspond to the main zones as determined from the assemblages and the PC analyses. The ages of the transitions were estimated from the age versus depth relationship shown in Figure 3.

Lithology

Visual sediment description allowed to distinguish three lithological units described by Ivanova et al. (Reference Ivanova, Murdmaa, de Vernal, Risebrobakken, Peyve, Brice, Seitkalieva and Pisarev2019; Fig. 4). Unit I spans the upper 223 cm of the core and consists of dark gray mud with a greenish shade with the appearance of authigenic iron sulfide. Unit II, between 223 and 326 cm, consists of a homogeneous brown mud. Unit III, from 326 to 385 cm, is represented by a diamicton with a dusky brown silty-sandy muddy matrix containing unsorted sand to pebble-size grains and rock fragments.

The coarse fraction (>106 μm), which was examined under a binocular microscope, is mostly composed of mineral particles with very rare microfossils. The dry-weight percentage of this fraction (Fig. 4) ranges between 0.3 and 2%, with exceptionally high values, up to 50%, at the base of the core (330–375 cm).

Palynomorph assemblages

Palynomorphs are generally well-preserved, but their abundance is highly variable. In most samples, dinocysts are the most abundant palynomorphs. Their concentration ranges from 1000 to 10,000 cysts/cm3 (Fig. 4), except in the lower part of the core (>213 cm; > ~11,300 cal yr BP), where the concentrations are very low (80–300 cysts/cm3).

Pollen and spore concentrations are low throughout the core, with values of 100–1000 grains/cm3. Foraminifer lining concentrations range from 500 to 2000 linings/cm3, except at the bottom of the core (> 310 cm) where they are close to zero. The interval below 310 cm is also characterized by the highest percentages of the coarse fraction. Below 213 cm, concentrations of reworked palynomorphs are high (1100–2700/cm3) thus reflecting erosion from adjacent land.

The dinocysts assemblages are dominated by taxa typical of subarctic environments, notably Islandinium minutum, Brigantedinium spp., Operculodinium centrocarpum, Nematosphaeropsis labyrinthus, and Impagidinium pallidum (Matthiessen et al., Reference Matthiessen, de Vernal, Head, Okolodkov, Zonneveld and Harland2005; de Vernal et al., Reference de Vernal, Henry, Matthiessen, Mudie, Rochon, Boessenkool and Eynaud2001, Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a). They are accompanied by Spiniferites elongatus, the cyst of Pentapharsodinium dalei, and Islandinium? cezare (Fig. 5).

The analysis performed with a matrix of 22 species and 73 samples indicate that the response data have a gradient length of 1.9 SD units, suggesting a linear relationship and that a principal component (PC) analysis is the appropriate ordination method. The results from PC analysis show that the first and second components explain 45.2 and 24.5% of the variance, respectively. The first component (PC1) shows an opposition between heterotrophic taxa (I. minutum and Brigantedinium spp., notably) on the negative side and phototrophic taxa (notably O. centrocarpum, N. labyrinthus, cyst of P. dalei, I. pallidum, and S. elongatus) on the positive side. The second component (PC2) shows an opposition between I. minutum and Brigantedinium spp., which respectively record negative and positive scores (Fig. 5).

The dinocyst concentrations, the relative abundance of the main taxa, and multivariate analysis led us to define three main zones.

The base of the core (215–375 cm), which is older than ca. 11,300 cal yr BP, is characterized by very low dinocyst concentrations, the dominance of Brigantedinium spp., and negative PC1 suggesting very harsh environments with dense sea-ice cover and only occasional and short-lived periods of open water favorable to primary productivity (Rochon et al., Reference Rochon, de Vernal, Turon, Matthiessen and Head1999; de Vernal et al., Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a, Reference de Vernal, Hillaire-Marcel, Rochon, Fréchette, Henry, Solignac and Bonnet2013b).

The interval between 215 and 100 cm, spanning from ca. 11,300 to 8700 cal yr BP is characterized by high concentrations of heterotrophic taxa and the large dominance of I. minutum, and by negative scores of both PC1 and PC2 (Fig. 5). Although I. minutum is not necessarily produced in sea ice (Heikkilä et al., Reference Heikkilä, Pospelova, Hochheim, Kuzyk, Stern, Barber and Macdonald2014, Reference Heikkilä, Pospelova, Forest, Stern, Fortier and Macdonald2016), it may live in ice (Potvin et al., Reference Potvin, Rochon and Lovejoy2013) and its high dominance in a context of low species diversity together with occurrence of cold-water taxa I.? cezare reflect an environment with dense seasonal sea-ice cover (Rochon et al., Reference Rochon, de Vernal, Turon, Matthiessen and Head1999; de Vernal et al., Reference de Vernal, Henry, Matthiessen, Mudie, Rochon, Boessenkool and Eynaud2001, Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a, Reference de Vernal, Radi, Zaragosi, Van Nieuwenhove, Rochon, Allan and De Schepper2020).

At ~100 cm, there is a sharp change with a transition from heterotrophic to phototrophic dominant assemblages, as also reflected by the switch from negative to positive PC1 scores. This change is dated at about 8700 cal yr BP. The cosmopolitan taxon O. centrocarpum (de Vernal and Marret, Reference de Vernal and Marret2007; Zonneveld et al., Reference Zonneveld, Marret, Versteegh, Bonnet, Bouimetarhan, Crouch and de Vernal2013; de Vernal et al., Reference de Vernal, Radi, Zaragosi, Van Nieuwenhove, Rochon, Allan and De Schepper2020) that occurs throughout the core largely dominates in the upper 100 cm where it represents 40–60% of the assemblages. N. labyrinthus, S. elongatus, I. pallidum, and the cyst of P. dalei, which are typical of subpolar environments of the northeast North Atlantic (Bonnet et al., Reference Bonnet, de Vernal, Hillaire-Marcel, Radi and Husum2010; de Vernal et al., Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a, Reference de Vernal, Radi, Zaragosi, Van Nieuwenhove, Rochon, Allan and De Schepper2020) are common accompanying species above the transition, whereas there is a near disappearance of the cold-water I.? cezare and a sharp decrease of I. minutum (Fig. 5). N. labyrinthus, S. elongatus, I. pallidum, and the cyst of P. dalei together point to relatively high productivity in surface waters characterized by seasonal sea-ice-free conditions. N. labyrinthus and I. pallidum further suggest open ocean conditions and relatively high salinity (Bonnet et al., Reference Bonnet, de Vernal, Hillaire-Marcel, Radi and Husum2010; de Vernal et al., Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a, Reference de Vernal, Radi, Zaragosi, Van Nieuwenhove, Rochon, Allan and De Schepper2020), whereas the cyst of P. dalei, which is associated with the spring bloom, corresponds to a relatively long sea-ice-free season and mild summer temperature in the surface waters (e.g., Dale, Reference Dale2001; Harland et al., Reference Harland, Nordberg and Filipsson2004; Bonnet et al., Reference Bonnet, de Vernal, Hillaire-Marcel, Radi and Husum2010; de Vernal et al., Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a, Reference de Vernal, Radi, Zaragosi, Van Nieuwenhove, Rochon, Allan and De Schepper2020). Between about 4000 and 2200 cal yr BP, maximum absolute and relative abundances of phototrophic taxa, and relatively high diversity of species, probably illustrate optimal conditions.

Reconstruction of sea-surface conditions

MAT was applied to all samples, but we only report on results from the upper two assemblage zones (<11,300 cal yr BP), which contain abundant dinocysts (concentration >1000 cysts/cm3). The results are based on close analogues for all samples. The distance between fossil and modern analogues mostly varied from 0.1 to 0.4 and never exceeded 0.6, which is less than half that of the threshold for poor analogues. The reliability index yielded values ranging 12 to 16 for most samples (Fig. 6). Thus, the reconstructions from the MAT are considered very reliable, especially in the upper part of the core above 80 cm.

Figure 6. Sea-surface reconstruction from the application of MAT to the dinocyst assemblages of core S2528. Sea-surface temperature and salinity in summer are represented by red curves, and sea-ice cover and annual productivity are represented by purple and green lines, respectively. The lighter areas correspond to maximum and minimum possible temperatures as calculated from the set of five modern analogues. The gray dashed lines correspond to smoothed values according to a five-points window. The horizontal rectangles at the top of the graphs correspond to the mean ± 1 SD of summer temperature and salinity (World Ocean Atlas) and sea-ice cover (National Snow and Ice Data Center data). For productivity, no modern direct data are available for statistics. The reliability index, which is calculated based on counts and the statistical distance between the fossil assemblage and modern analogues, indicates high to very high reliability for most samples. The horizontal dashed lines represent the transitions discussed in the text, with tentative chronology estimated from the age versus depth relationship of Figure 3.

The reconstruction from MAT revealed important changes in sea-surface conditions (Fig. 6). From ca. 11,300 to 8700 cal yr BP, very cold conditions with an average of 8 months/yr of sea-ice cover, summer SST around 1–2°C and SSS fluctuating between 28 and 33.5 psu are reconstructed. Such conditions indicate polar climate and episodic meltwater discharge.

A warming trend is observed after ca. 8700 cal yr BP. Summer SST and SSS increased progressively to a maximum of 3.5°C and 33.5 psu, respectively, whereas seasonal sea ice cover decreased to a minimum of 2–4 months/yr, which is much lower than the modern range of values. High-frequency variations in sea-ice cover are observed until ca. 4000 cal yr BP and highest temperature and salinity were reached between 3000 and 2200 cal yr BP. After this optimum, a cooling trend is recorded with a summer SST drop of ~2°C, and a sea-ice cover increase up to 8 months/yr, with a trend towards modern values.

DISCUSSION

The glacial termination (> 11,300 cal yr BP)

The timing and structure of the last deglaciation of the northern Eurasian ice sheet has been documented by many studies (see the review by Hughes et al. [2016] and references therein; Murdmaa and Ivanova, Reference Murdmaa, Ivanova and Boone2017), notably in the northern Barents Sea sector (Hogan et al., Reference Hogan, Dowdeswell, Hillenbrand, Ehrmann, Noormets and Wacker2017). During the last glacial maximum, the Svalbard-Barents-Kara Ice Sheet occupied the northern Barents Sea including the present study area (e.g., Hughes et al., Reference Hughes, Gyllencreutz, Lohne, Mangerud and Svendsen2016). Along the western and northern Spitsbergen coasts, marine sedimentological data have shown that the ice retreat started around 15,000 14C y BP, which corresponds to about 17,700 cal yr BP (Mangerud et al., Reference Mangerud, Bolstad, Elgersma, Helliksen, Landvik, Lønne, Lycke, Salvigsen, Sandahl and Svendsen1992; Svendsen and Mangerud, Reference Svendsen and Mangerud1997; Landvik et al., Reference Landvik, Bondevik, Elverhøi, Fjeldskaar, Mangerud, Salvigsen and Siegert1998; Rasmussen et al., Reference Rasmussen, Thomsen, Ślubowska, Jessen, Solheim and Koç2007; Ślubowska-Woldengen et al., Reference Ślubowska-Woldengen, Rasmussen, Koç, Klitgaard-Kristensen, Nilsen and Solheim2007). High-resolution seismic profiles and lithostratigraphical data from the northern Barents Sea indicated that the grounded ice in the Franz Victoria Trough started retreating as early as 15,400 14C yr BP (= 18,100 cal yr BP; Kleiber et al., Reference Kleiber, Knies and Niessen2000). The Hinlopen Trough grounded ice retreated between 13,900 and 13,700 14C yr BP (~16,000–15,900 cal yr BP; Koç et al., Reference Koç, Klitgaard-Kristensen, Hasle, Forsberg and Solheim2002) and the area south of Kvitøya was ice-free at ca. 14,900 cal yr BP (Klitgaard Kristensen et al., Reference Klitgaard Kristensen, Rasmussen and Koç2013). Most of the Barents Sea was deglaciated by 13,000 14C yr BP (= 14,800 cal yr BP; Elverhøi et al., Reference Elverhøi, Fjeldskaar, Solheim, Nyland-Berg and Russwurm1993) with the Atlantic water flowing to the north via the Svalbard Branch and the North Cape Current. According to Hogan et al. (Reference Hogan, Dowdeswell, Hillenbrand, Ehrmann, Noormets and Wacker2017), the completely ice-free Franz Victoria Trough at 16,500 cal yr BP set off the deglaciation westward towards Spitsbergen via the Erik Eriksen Trough but grounded ice persisted in the southern part of Kvitøya Trough until ca. 13,500 cal yr BP.

The lower part of the core record, which is discussed in detail by Ivanova et al. (Reference Ivanova, Murdmaa, de Vernal, Risebrobakken, Peyve, Brice, Seitkalieva and Pisarev2019), indicates that the Kvitøya Trough was deglaciated prior to 15,500 cal yr BP. However, the interval older than 11,300 cal yr BP was characterized by extremely low concentrations of dinocysts, suggesting very low productivity, which we tentatively associate with harsh conditions and very dense sea-ice cover. The rare dinocyst specimens recovered mostly belong to Brigantedinium spp. and I. minutum, which presently occur together in marine environments with dense seasonal ice cover (Rochon et al., Reference Rochon, de Vernal, Turon, Matthiessen and Head1999; de Vernal et al., Reference de Vernal, Henry, Matthiessen, Mudie, Rochon, Boessenkool and Eynaud2001, Reference de Vernal, Hillaire-Marcel and Darby2005, Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a). The high concentrations of reworked palynomorphs and IRD in this interval suggest glaciomarine-type environments during a phase that might correspond to the deglaciation of the Svalbard margins (Knies et al., Reference Knies, Vogt and Stein1998). Whereas cold conditions prevailed in the surface waters of the northern Barents Sea, relatively warm Atlantic waters may have characterized the subsurface layers as suggested by foraminifer assemblages and the isotopic composition of benthic species (Lubinski et al., Reference Lubinski, Korsun, Polyak, Forman, Lehman, Herlihy and Miller1996, Reference Lubinski, Polyak and Forman2001; Koç et al., Reference Koç, Klitgaard-Kristensen, Hasle, Forsberg and Solheim2002; Klitgaard Kristensen et al., Reference Klitgaard Kristensen, Rasmussen and Koç2013) including those at the study site (Ivanova et al., Reference Ivanova, Murdmaa, de Vernal, Risebrobakken, Peyve, Brice, Seitkalieva and Pisarev2019). Hence, this interval was likely characterized by strongly stratified upper waters with dense sea-ice cover, overlying relatively warm Atlantic water subducted below the surface layer that may have contributed to the deglaciation of the Svalbard margins.

Transition towards “interglacial” conditions at about 8700 ± 700 cal yr BP

The occurrence of S. elongatus, cyst of P. dalei, and I.? cezare, the increased percentages of I. pallidum, N. labyrinthus, and O. centrocarpum together with occurrence of Brigantedinium spp. and the maximum abundance of I. minutum mark the end of glaciomarine-type conditions at ca. 11,300 cal yr BP. However, the large dominance of I. minutum and relatively high percentages of I.? cezare until 8700 cal yr BP suggest the persistence of cold conditions and dense sea-ice cover nearby (Rochon et al., Reference Rochon, de Vernal, Turon, Matthiessen and Head1999; de Vernal et al., Reference de Vernal, Henry, Matthiessen, Mudie, Rochon, Boessenkool and Eynaud2001, Reference de Vernal, Rochon, Fréchette, Henry, Radi and Solignac2013a). The application of MAT led to reconstruct low SSTs and 8–10 months/yr sea-ice cover (Fig. 6). It also suggests low salinity (<32 psu in summer), which could relate to dilution with freshwater/meltwater inputs from sea ice and/or remotely located ice or glaciers. Micropaleontological and biomarker data from northern and eastern Svalbard fjords also suggest high meltwater fluxes, which were attributed to glacier retreat fostered by high summer insolation during the early Holocene (Bartels et al., Reference Bartels, Titschack, Fahl, Stein, Seidenkrantz, Hillaire-Marcel and Hebbeln2017, Reference Bartels, Titschack, Fahl, Stein and Hebbeln2018). Although the timing can be challenged, the interval encompassing 8500–7500 cal yr BP was characterized by an important change related to the freshwater budget of the high latitude ocean. It was the time of the postglacial eustatic sea level rise deceleration marking the end of the main deglaciation (e.g., Peltier and Fairbanks, Reference Peltier and Fairbank2006), which corresponds to the final dislocation of the Laurentide Ice Sheet (Dyke et al., Reference Dyke, Moore and Robertson2003) and diminution of the Greenland Ice Sheet volume (Larsen et al., Reference Larsen, Kjær, Lecavalier, Bjørk, Colding, Huybrechts and Jakobsen2015). In our core located in the northern Barents Sea, the surface ocean transition recorded at 8700 ± 700 cal yr BP might thus correspond to large scale change in the circum-Arctic freshwater budget at the end of the glaciation.

Attainment of optimal conditions in surface waters during the Holocene

After ca. 8700 cal yr BP, the dinocyst assemblages dominated by heterotrophic taxa (notably I. minutum) were replaced by assemblages dominated by phototrophic taxa, including O. centrocarpum, N. labyrinthus, S. elongatus, I. pallidum, and the cyst of P. dalei (Fig. 5). This corresponds to the establishment of interglacial conditions in the surface waters at the coring site. From 8700 to about 4000 cal BP, there is a gradual decrease of I. minutum relative to phototrophic taxa (Fig. 5). Whereas concentrations of heterotrophic taxa remain similar, there is an increase of absolute abundance of phototrophic taxa, pointing to overall increase of productivity. Furthermore, the augmentation of I. pallidum, which is characteristic of offshore conditions in the subpolar Nordic Seas, suggests enhanced salinity (Bonnet et al., Reference Bonnet, de Vernal, Hillaire-Marcel, Radi and Husum2010; de Vernal et al., Reference de Vernal, Radi, Zaragosi, Van Nieuwenhove, Rochon, Allan and De Schepper2020). MAT led to reconstruct a gradual increase of SSTs and salinity and decrease of sea-ice cover (Fig. 6), with optimal conditions in surface waters being attained around 4100 ± 1500 yr BP, with low sea-ice cover for a couple of months per yr and summer temperatures reaching up to 4°C, which is warmer than modern conditions. The timing of events is poorly constrained from our core. Nevertheless, attainment of optimal conditions in surface waters during the mid-late Holocene was recorded from dinocyst data of many sites in the subpolar North Atlantic, notably in the eastern Nordic Seas (Van Nieuwenhove et al., Reference Van Nieuwenhove, Baumann, Matthiessen, Bonnet and de Vernal2016) and possibly in the Baffin Bay (Gibb et al., Reference Gibb, Steinhauer, Fréchette, de Vernal and Hillaire-Marcel2015). Other paleoceanographic records from the southwest Barents Sea and eastern Fram Strait rather suggest that the thermal optimum corresponds to an interval dated from the early Holocene to about 7500 cal yr BP (Baumann and Matthiessen, Reference Baumann and Matthiessen1992; Bauch et al., Reference Bauch, Erlenkeuser, Spielhagen, Struck, Matthiessen, Thiede and Heinemeier2001; Sarnthein et al., Reference Sarnthein, Kreveld, Erlenkeuser, Grootes, Kucera, Pflaumann and Schulz2003; Hald et al., Reference Hald, Ebbesen, Forwick, Godtliebsen, Khomenko, Korsun, Ringstad Olsen and Vorren2004, Reference Hald, Andersson, Ebbesen, Jansen, Klitgaard-Kristensen, Risebrobakken, Salomonsen, Sarnthein, Sejrup and Telford2007; Risebrobakken et al., Reference Risebrobakken, Dokken, Smedsrud, Andersson, Jansen, Moros and Ivanova2011; Müller et al., Reference Müller, Werner, Stein, Fahl, Moros and Jansen2012; Rigual-Hernàndez et al., Reference Rigual-Hernández, Colmenero-Hidalgo, Martrat, Bárcena, de Vernal, Sierro, Flores, Grimalt, Henry and Lucchi2017; Falardeau et al., Reference Falardeau, de Vernal and Spielhagen2018). In the Franz Victoria Trough, northeast Barents Sea, isotope data from benthic and planktic foraminifers indicate that optimal conditions in the bottom waters occurred at 7800–6800 cal yr BP, whereas warming in the subsurface waters was recorded until they reached present-day values at ca. 3700 cal yr BP (Duplessy et al., Reference Duplessy, Ivanova, Murdmaa, Paterne and Labeyrie2001).

The timing of the Holocene climate optimum has been the matter of much debate. In the circum-Arctic, it appears diachronous from one region to another (e.g., Kaufman et al., Reference Kaufman, Ager, Anderson, Anderson, Andrews, Bartlein and Brubaker2004). In Scandinavia, most terrestrial information available tends to indicate warmer than average conditions from the early Holocene to about 4000 yr ago (Sejrup et al., Reference Sejrup, McKay, Kaufman, Geirsdottir, de Vernal, Renssen, Husum, Jennings and Andrews2016). The marine data from the Nordic Seas, however, show variations with discrepancies from one site to another and from one proxy to another (Sejrup et al., Reference Sejrup, McKay, Kaufman, Geirsdottir, de Vernal, Renssen, Husum, Jennings and Andrews2016). For example, in the subpolar North Atlantic, surface ocean records based on alkenones show maximum surface temperature in the early–mid Holocene generally followed by decreasing trends after about 7000 cal yr BP, whereas records based on planktonic foraminifers show variations that make it difficult to identify clear trends (Leduc et al., Reference Leduc, Schneider, Kim and Lohmann2010; Risebrobakken et al., Reference Risebrobakken, Dokken, Smedsrud, Andersson, Jansen, Moros and Ivanova2011). In the northeast subpolar North Atlantic, the general discrepancies between temperature estimates from alkenones and those from planktic foraminifers has been interpreted as the differential signal between the surface ocean under direct influence of solar forcing and the subsurface ocean marked by heat advection and redistribution through currents (Risebrobakken et al., Reference Risebrobakken, Dokken, Smedsrud, Andersson, Jansen, Moros and Ivanova2011).

Studies conducted around Svalbard agree on the optimal conditions in subsurface water masses starting early, at ca. 10,000 cal yr BP if not earlier, as suggested from benthic and planktic foraminifer data (Hald et al., Reference Hald, Ebbesen, Forwick, Godtliebsen, Khomenko, Korsun, Ringstad Olsen and Vorren2004, Reference Hald, Andersson, Ebbesen, Jansen, Klitgaard-Kristensen, Risebrobakken, Salomonsen, Sarnthein, Sejrup and Telford2007; Ślubowska et al., Reference Ślubowska, Koç, Rasmussen and Klitgaard-Kristensen2005; Ebbesen et al., Reference Ebbesen, Hald and Eplet2007; Rasmussen et al., Reference Rasmussen, Thomsen, Ślubowska, Jessen, Solheim and Koç2007, Reference Rasmussen, Thomsen, Skirbekk, Slubowska-Woldengen, Klitgaard Kristensen and Koç2014; Ślubowska-Woldengen et al., Reference Ślubowska-Woldengen, Rasmussen, Koç, Klitgaard-Kristensen, Nilsen and Solheim2007, Reference Ślubowska-Woldengen, Koç, Rasmussen, Klitgaard-Kristensen, Hald and Jennings2008; Werner et al., Reference Werner, Spielhagen, Bauch, Hass and Kandiano2013; Bartels et al., Reference Bartels, Titschack, Fahl, Stein, Seidenkrantz, Hillaire-Marcel and Hebbeln2017). This is consistent with strengthening of the North Atlantic Current and a change in the regional ocean circulation in the North Atlantic (e.g., Jansen et al., Reference Jansen, Andersson, Moros, Nisancioglu, Nyland, Telford, Battarbee and Binney2009) together with maximum ocean heat transport following reorganization of the AMOC, as also inferred from paleoceanographic reconstruction in the Nordic Seas and Western Barents Sea (Risebrobakken et al., Reference Risebrobakken, Dokken, Smedsrud, Andersson, Jansen, Moros and Ivanova2011).

In the northern Barents Sea area, the late attainment of optimal conditions in the surface waters as recorded at about 4000 cal yr BP from the dinocyst data in core S2528 is accompanied by maximum salinity, which suggests minimum dilution with freshwater and meltwater. This optimum corresponds with maximum occurrences of phototrophic taxa including S. ramosus in association with O. centrocarpum and N. labyrinthus, which depicts an oceanic character to the assemblage and suggests dominant Atlantic influence. Hence, during the regional optimum, likely encompassing the 4000 to 2200 cal yr BP interval, the Atlantic water may have flowed close to the surface as far northeast as the core site. If our interpretation is correct, the apparent delay in optimal conditions could be an indirect consequence of insolation forcing during the early–middle Holocene, which contributed to enhanced ice melt in the northern Barents Sea, thus upper water mass stratification restricting the warm but saline Atlantic water to subsurface layers until about 4000 yr ago.

The late Holocene climate deterioration

After the relatively warm conditions that lasted until about 2200 cal yr BP, there was a cooling trend as estimated from the MAT applied to dinocysts. Sea-ice cover increased to >8 months/yr and summer SST decreased from 3.5 to ~1°C. Other records from the Nordic and the Barents Seas indicate a late Holocene deterioration of temperatures in both surface and bottom waters, which was associated with the diminution of summer insolation and the reduction in the northeastward flux of NAW (Hald et al., Reference Hald, Danielsen and Lorentzen1989, Reference Hald, Ebbesen, Forwick, Godtliebsen, Khomenko, Korsun, Ringstad Olsen and Vorren2004; Koç et al., Reference Koç, Jansen and Haidason1993; Polyak and Mikhailov, Reference Polyak, Mikhailov, Andrews, Austin, Bergsten and Jennings1996; Hald and Aspeli, Reference Hald and Aspeli1997; Jennings et al., Reference Jennings, Knudsen, Hald, Hansen and Andrews2002; Rasmussen et al., Reference Rasmussen, Bäckström, Heinemeier, Klitgaard-Kristensen, Knutz, Kuijpers, Lassen, Thomsen, Troelstra and van Weering2002; Sarnthein et al., Reference Sarnthein, Kreveld, Erlenkeuser, Grootes, Kucera, Pflaumann and Schulz2003a; Andersen et al., Reference Andersen, Koç, Jennings and Andrews2004; Bauch et al., Reference Bauch, Erlenkeuser, Spielhagen, Struck, Matthiessen, Thiede and Heinemeier2001; Ślubowska et al., Reference Ślubowska, Koç, Rasmussen and Klitgaard-Kristensen2005, Reference Ślubowska-Woldengen, Rasmussen, Koç, Klitgaard-Kristensen, Nilsen and Solheim2007; Rasmussen et al., Reference Rasmussen, Thomsen, Ślubowska, Jessen, Solheim and Koç2007; Ślubowska-Woldengen, Reference Ślubowska-Woldengen, Koç, Rasmussen, Klitgaard-Kristensen, Hald and Jennings2008; Risebrobakken et al., Reference Risebrobakken, Moros, Ivanova, Chistyakova and Rosenberg2010, Reference Risebrobakken, Dokken, Smedsrud, Andersson, Jansen, Moros and Ivanova2011). Hence, almost all data from the Nordic Seas and Scandinavia seem consistent both on land and at sea, showing that the last two millennia have been regionally the coldest of the Holocene (Sejrup et al., Reference Sejrup, McKay, Kaufman, Geirsdottir, de Vernal, Renssen, Husum, Jennings and Andrews2016), excluding the recent anthropogenic warming.

CONCLUSIONS

Palynological analysis from core S2528 allowed us to document surface conditions and the interactions between North Atlantic and Arctic waters northeast of Svalbard, which permit us to discuss their respective role in the regional hydrography from the deglaciation until present-day. Cold and harsh conditions lasted until 11,300 cal yr BP with dense sea-ice cover and strong stratification of water masses, restricting inflow of the warm and saline Atlantic water to subsurface layers. While the Nordic Seas and Fram Strait were ice-free since the last glacial maximum, the retreat of the Svalbard-Barents Sea ice sheet during the deglaciation and early Holocene was accompanied by large meltwater discharge that resulted in strong stratification and too harsh conditions in the surface waters of the northern Barents Sea for high dinoflagellate productivity. In the study area, the deglaciation accelerated at about 11,300 cal yr BP, with marked inflow of Atlantic water from the Svalbard Branch and pulses of meltwater from the Svalbard ice sheet resulting in generally low and highly variable surface salinity.

The first transition towards interglacial conditions in surface waters were recorded at about 8700 cal yr BP, with a very abrupt change in dinocysts assemblages, from dominant I. minutum to assemblages composed mainly of O. centrocarpum and N. labyrinthus. This transition suggests an increase of SSTs, salinity and productivity, and a significant decrease of sea-ice cover. After this main transition and until about 4000 cal yr BP, large amplitude variations are recorded. Sea-surface conditions finally stabilize to reach optimal conditions with respect to salinity, temperature, and productivity between ca. 3000 and 2200 cal yr BP. Although our core was collected offshore, far from human settlement, the data compiled here and the reconstruction we propose provide some clues on the regional climate and marine environment history. Hence, our study yields some information about the fact that the interval around 3000–2000 cal yr BP was probably the most favorable time window of the Holocene for human occupation in the area of northern Svalbard. Such a finding does not support directly the hypothesis of Stone Age settlements in Svalbard proposed by Christiansson and Simonsen (Reference Christiansson and Simonsen1970) but provides a climatological and oceanographical framework that would make it feasible. The trend towards colder conditions that is inferred over the last 2000 yr of the Holocene, prior to the ongoing anthropogenic warming, in the northwestern Barents Sea was probably an important limitation as in most other areas of the Arctic-subarctic sector of the North Atlantic.

ACKNOWLEDGMENTS

We thank the scientific party of Cruise 25 of R/V Akademik Nikolaj Strakhov for offering the samples from core S2528. This study was supported by the Fonds de Recherche du Québec sur la Nature et les Technologies (FRQNT) and the Natural Sciences and Engineering Research Council (NSERC) of Canada (grants to AdV) and by the Russian Federation state (assignment No. 0149-2019-0007 to EI). The laboratory analyses have been made at Geotop. We are grateful to the guest editor, Nicole Mirsati, Manuel Bringué, and one anonymous reviewer of the journal for their critical and constructive comments, which helped to prepare the revised version of the manuscript.

SUPPLEMENTARY MATERIAL

The supplementary material for this article can be found at https://doi.org/10.1017/qua.2020.2.

References

REFERENCES

Albrethsen, S.E., Arlov, T.B., 1988. The discovery of Svalbard—A problem reconsidered. Fennoscandia Archaeologica 5, 105110.Google Scholar
Andersen, C., Koç, N., Jennings, A., Andrews, J.T., 2004. Non-uniform response of the major surface currents in the Nordic Seas to insolation forcing: Implications for the Holocene climate variability. Paleoceanography 19, PA2003. http://dx.doi.org/10.1029/2002PA000873.CrossRefGoogle Scholar
Arlov, T.B., 2005. The discovery and early exploitation of Svalbard: Some historiographical notes. Acta Borealia 22, 319.CrossRefGoogle Scholar
Årthun, M., Eldevik, T., Smedsrud, L. H., Skagseth, Ø., Ingvaldsen, R.B., 2012. Quantifying the influence of Atlantic heat on Barents Sea ice variability and retreat. Journal of Climate 25, 47364743.CrossRefGoogle Scholar
Årthun, M., Schrum, C., 2010. Ocean surface heat flux variability in the Barents Sea. Journal Marine Systems 83: 8898.CrossRefGoogle Scholar
Bartels, M., Titschack, J., Fahl, K., Stein, R., Hebbeln, D., 2018. Wahlenbergfjord, eastern Svalbard: A glacier-surrounded fjord reflecting regional hydrographic variability during the Holocene? Boreas 47, 10031021CrossRefGoogle Scholar
Bartels, M., Titschack, J., Fahl, K., Stein, R., Seidenkrantz, M-S., Hillaire-Marcel, C., Hebbeln, D., 2017. Atlantic water advection vs. glacier dynamics in northern Spitsbergen since early deglaciation. Climate of the Past 13, 17171749.CrossRefGoogle Scholar
Bauch, H.A., Erlenkeuser, H., Spielhagen, R.F., Struck, U., Matthiessen, J., Thiede, J., Heinemeier, J., 2001. A multiproxy reconstruction of the evolution of deep and surface waters in the subarctic Nordic seas over the last 30,000 yr. Quaternary Science Reviews 20, 659678.CrossRefGoogle Scholar
Baumann, K.-H., Matthiessen, J., 1992. Variations in surface water mass conditions in the Norwegian Sea: evidence from Holocene coccolith and dinoflagellate cyst assemblages. Marine Micropaleontology 20, 129146.CrossRefGoogle Scholar
Bjerk, H.B., 2000. Stone Age settlement on Svalbard? A re-evaluation of previous finds and the results of a recent field survey. Polar Record 36, 97112.CrossRefGoogle Scholar
Blaauw, M., 2010. “Methods and code for ’classical’ age-modelling of radiocarbon sequences.” Quaternary Geochronology 5, 512518.CrossRefGoogle Scholar
Blaauw, M., Christen, J.A., 2011. Flexible paleoclimate age-depth models using an autoregressive gamma process. Bayesian Analysis 6, 457474.CrossRefGoogle Scholar
Bonnet, S., de Vernal, A., Hillaire-Marcel, C., Radi, T., Husum, K., 2010. Variability of sea-surface temperature and sea-ice cover in the Fram Strait over the last two millennia. Marine Micropaleontology 74, 5974.CrossRefGoogle Scholar
Christiansson, H., Simonsen, P., 1970. Stone Age finds from Spitsbergen. Acta Borealia, B. Humaniora 11, 524.Google Scholar
Dale, B., 2001. The sedimentary record of dinoflagellate cysts: Looking back into the future of phytoplankton blooms. Scientia Marina 65, 257272.CrossRefGoogle Scholar
de Vernal, A., Bilodeau, G., Henry, M., 2010. Micropaleontological Preparation Techniques and Analyses. Cahier du Geotop n° 3, Geotop, Montréal.Google Scholar
de Vernal, A., Henry, M., Matthiessen, J., Mudie, P.J., Rochon, A., Boessenkool, K.P., Eynaud, F. et al. , 2001. Dinoflagellate cyst assemblages as tracers of sea-surface conditions in the northern North Atlantic, Arctic and sub-Arctic seas: The new ‘n= 677’data base and its application for quantitative palaeoceanographic reconstruction. Journal of Quaternary Science 16, 681698.CrossRefGoogle Scholar
de Vernal, A., Hillaire-Marcel, C., Darby, D.A., 2005. Variability of sea ice cover in the Chukchi Sea (western Arctic Ocean) during the Holocene. Paleoceanography 20, PA4018. http://dx.doi.org/10.1029/2005PA001157.CrossRefGoogle Scholar
de Vernal, A., Hillaire-Marcel, C., Rochon, A., Fréchette, B., Henry, M., Solignac, S., Bonnet, S., 2013b. Dinocyst-based reconstructions of sea ice cover concentration during the Holocene in the Arctic Ocean, the northern North Atlantic Ocean and its adjacent seas. Quaternary Science Reviews 79, 111121.CrossRefGoogle Scholar
de Vernal, A., Marret, F., 2007. Organic-walled dinoflagellate cysts: Tracers of sea-surface conditions. Proxies in Late Cenozoic Paleoceanography 1, 371408.CrossRefGoogle Scholar
de Vernal, A., Radi, T., Zaragosi, S., Van Nieuwenhove, N., Rochon, A., Allan, E., De Schepper, S., et al. , 2020. Distribution of common modern dinoflagellate cyst taxa in surface sediments of the Northern Hemisphere in relation to environmental parameters: The new n = 1968 database. Marine Micropaleontology in press.CrossRefGoogle Scholar
de Vernal, A., Rochon, A., Fréchette, B., Henry, M., Radi, T., Solignac, S., 2013a. Reconstructing past sea ice cover of the Northern hemisphere from dinocyst assemblages: status of the approach. Quaternary Science Reviews 79, 122134.CrossRefGoogle Scholar
Dieckmann, G.S., Hellmer, H.H., 2008. The importance of sea ice: An overview. In: Thomas, D.N., Dieckmann, G.S (Eds.), Sea Ice: An Introduction to its Physics, Chemistry, Biology and Geology. Blackwell Science, Oxford. http://dx.doi.org/10.1002/9780470757161.ch1.Google Scholar
Duplessy, J.C., Ivanova, E., Murdmaa, I., Paterne, M., Labeyrie, L., 2001. Holocene paleoceanography of the northern Barents Sea and variations of the northward heat transport by the Atlantic Ocean. Boreas 30, 216.CrossRefGoogle Scholar
Dyke, A.S., Moore, A., Robertson, L., 2003. Deglaciation of North America. Geological Survey of Canada, Open File 1574. https://doi.org/10.4095/214399CrossRefGoogle Scholar
Ebbesen, H., Hald, M., Eplet, T.H., 2007. Late glacial and early Holocene climatic oscillations on the western Svalbard margin, European Arctic. Quaternary Science Review 26, 19992011.CrossRefGoogle Scholar
Elverhøi, A., Fjeldskaar, W., Solheim, A., Nyland-Berg, M., Russwurm, L., 1993. The Barents Sea Ice Sheet—A model of its growth and decay during the last ice maximum. Quaternary Science Reviews 12, 863873.CrossRefGoogle Scholar
Falardeau, J., de Vernal, A., Spielhagen, R., 2018. Paleoceanography of northeastern Fram Strait since the last glacial maximum: palynological evidence of large amplitude change. Quaternary Science Reviews, 195: 133152.CrossRefGoogle Scholar
Gammelsrød, T., Leikvin, Ø., Lien, V., Budgell, W.P., Loeng, H., Maslowski, W., 2009. Mass and heat transports in the NE Barents Sea: observations and models. Journal of Marine Systems 75: 5669.CrossRefGoogle Scholar
Gibb, O. T., Steinhauer, S., Fréchette, B., de Vernal, A., Hillaire-Marcel, C., 2015. Diachronous evolution of sea surface conditions in the Labrador Sea and Baffin Bay since the last deglaciation. The Holocene 25, 18821897.CrossRefGoogle Scholar
Guiot, J., de Vernal, A., 2007. Transfer functions: Methods for quantitative paleoceanography based on microfossils. Developments in Marine Geology 1, 523563.CrossRefGoogle Scholar
Guiot, J., de Vernal, A., 2011. Is spatial autocorrelation introducing biases in the apparent accuracy of paleoclimatic reconstructions? Quaternary Science Reviews 30, 19651972.CrossRefGoogle Scholar
Hald, M., Andersson, C., Ebbesen, H., Jansen, E., Klitgaard-Kristensen, D., Risebrobakken, B., Salomonsen, G.R., Sarnthein, M., Sejrup, H.P., Telford, R.J., 2007. Variations in temperature and extent of Atlantic Water in the northern North Atlantic during the Holocene. Quaternary Science Reviews 26, 34233440.CrossRefGoogle Scholar
Hald, M., Andersson, C., Ebbesen, H., Jansen, E., Klitgaard-Kristensen, D., Risebrobakken, B., Salomonsen, G.R., Sarnthein, M., Sejrup, H.P., Telford, R.J., 2007. Variations in temperature and extent of Atlantic water in the northern North Atlantic during the Holocene. Quaternary Science Reviews 26, 34233440.CrossRefGoogle Scholar
Hald, M., Aspeli, R., 1997. Rapid climatic shifts on the northern Norwegian Sea during the last deglaciation and the Holocene. Boreas 26, 1528.CrossRefGoogle Scholar
Hald, M., Danielsen, T.K., Lorentzen, S., 1989. Late Pleistocene-Holocene benthic foraminiferal distribution in the southwestern Barents Sea: Paleoenvironmental implications. Boreas 18, 368388.Google Scholar
Hald, M., Ebbesen, H., Forwick, M., Godtliebsen, F., Khomenko, L., Korsun, S., Ringstad Olsen, L., Vorren, T.O., 2004. Holocene paleoceanography and glacial history of the West Spitsbergen area, Euro-Arctic margin. Quaternary Science Reviews 23, 20752088.CrossRefGoogle Scholar
Harland, R., Nordberg, K., Filipsson, H.L., 2004. A high-resolution dinoflagellate cyst record from latest Holocene sediments in Koljö Fjord, Sweden. Review of Palaeobotany and Palynology 128, 119141.CrossRefGoogle Scholar
Head, M.J., 1996. Modern dinoflagellate cysts and their biological affinities. In: Jansonius, J., McGregor, D.C. (Eds.), Palynology: Principles and Applications. Vol. 3. American Association of Stratigraphic Palynologists Foundation, Dallas, pp. 11971248.Google Scholar
Heikkilä, M., Pospelova, V., Forest, A., Stern, G.A., Fortier, L., Macdonald, R.W., 2016. Dinoflagellate cyst production over an annual cycle in seasonally ice-covered Hudson Bay. Marine Micropaleontology 125, 124.CrossRefGoogle Scholar
Heikkilä, M., Pospelova, V., Hochheim, K.P., Kuzyk, Z.Z.A., Stern, G.A., Barber, D.G., Macdonald, R.W., 2014. Surface sediment dinoflagellate cysts from the Hudson Bay system and their relation to freshwater and nutrient cycling. Marine Micropaleontology 106, 79109.CrossRefGoogle Scholar
Helland-Hansen, B., Nansen, F., 1909. The Norwegian Sea. Fiskerhdirektoratets Skrifter Series Hauunders 2, 1360.Google Scholar
Hogan, K.A., Dowdeswell, J.A., Hillenbrand, C.-D., Ehrmann, W., Noormets, R., Wacker, L., 2017. Subglacial sediment pathways and deglacial chronology of the northern Barents Sea Ice Sheet. Boreas 46, 750771.CrossRefGoogle Scholar
Hohmann, S., Kucera, M., de Vernal, A., 2020. Identifying the signature of sea-surface properties in dinocyst assemblages: Implications for quantitative palaeoceanographical reconstructions by transfer functions and analogue techniques. Marine Micropaleontology in press.CrossRefGoogle Scholar
Hughes, A.L.C., Gyllencreutz, R., Lohne, Ø.S., Mangerud, J., Svendsen, J.I., 2016. The last Eurasian ice sheets—A chronological database and time-slice reconstruction Boreas 45, 145.CrossRefGoogle Scholar
Ivanova, E., Murdmaa, I., de Vernal, A., Risebrobakken, B., Peyve, A., Brice, C., Seitkalieva, E., Pisarev, S., 2019. Postglacial variations in Atlantic Water inflow, iceberg calving and sea ice conditions in the northwestern Barents Sea. Quaternary Research 120. https://doi.org/10.1017/qua.2019.18Google Scholar
Jansen, E., Andersson, C., Moros, M., Nisancioglu, K.H., Nyland, B.F., Telford, R.J., 2009. The early to mid-Holocene thermal optimum in the North Atlantic. In: Natural Climate Variability and Global Warming: A Holocene Perspective. Battarbee, R.W. and Binney, H. A, Eds. Blackwell Publishing Ltd, pp. 123137.Google Scholar
Jennings, A.E., Knudsen, K.L., Hald, M., Hansen, C.V., Andrews, J.T., 2002. A mid-Holocene shift in Arctic sea-ice variability on the East Greenland Shelf. The Holocene 12, 4958.CrossRefGoogle Scholar
Jeong, H.J., Yoo, Y.D., Kim, J.S., Seong, K.A., Kang, N.S., Kim, T.H., 2010. Growth, feeding and ecological roles of the mixotrophic and heterotrophic dinoflagellates in marine planktonic food webs. Ocean Science Journal 45, 6591.CrossRefGoogle Scholar
Kaufman, D.S., Ager, T.A., Anderson, N.J., Anderson, P.M., Andrews, J.T., Bartlein, P. J., Brubaker, L.B. et al. , 2004. Holocene thermal maximum in the western Arctic (0–180 W). Quaternary Science Reviews 23, 529560.CrossRefGoogle Scholar
Kleiber, H.P., Knies, J., Niessen, F., 2000. The Late Weichselian glaciation of the Franz Victoria Trough, northern Barents Sea: Ice sheet extent and timing. Marine Geology 168, 2544.CrossRefGoogle Scholar
Klitgaard Kristensen, D.K., Rasmussen, T.L., Koç, N., 2013. Palaeoceanographic changes in the northern Barents Sea during the last 16 000 years—New constraints on the last deglaciation of the Svalbard-Barents Sea Ice Sheet. Boreas 42, 798813.CrossRefGoogle Scholar
Knies, J., Vogt, C., Stein, R., 1998. Late Quaternary growth and decay of the Svalbard/Barents Sea ice sheet and paleoceanographic evolution in the adjacent Arctic Ocean. Geo-Marine Letters 18, 195202.CrossRefGoogle Scholar
Koç, N., Jansen, E., Haidason, H., 1993. Paleoceanographic reconstructions of surface ocean conditions in the Greenland, Iceland and Norwegian seas through the last 14 ka based on diatoms. Quaternary Science Reviews 12, 115140.CrossRefGoogle Scholar
Koç, N., Klitgaard-Kristensen, D., Hasle, K., Forsberg, C.F., Solheim, A., 2002. Late glacial palaeoceanography of Hinlopen Strait, northern Svalbard. Polar Research 21, 307314.CrossRefGoogle Scholar
Landvik, J.Y., Bondevik, S., Elverhøi, A., Fjeldskaar, W., Mangerud, J.A.N., Salvigsen, , Siegert, M.J. et al. , 1998. The last glacial maximum of Svalbard and the Barents Sea area: Ice sheet extent and configuration. Quaternary Science Reviews 17, 4375.CrossRefGoogle Scholar
Larsen, N.K., Kjær, K.H., Lecavalier, B., Bjørk, A.A., Colding, S., Huybrechts, P., Jakobsen, K.E., et al. , 2015., The response of the southern Greenland ice sheet to the Holocene thermal maximum. Geology 43, 291294.CrossRefGoogle Scholar
Leduc, G., Schneider, R., Kim, J.-H., Lohmann, G., 2010. Holocene and Eemian sea surface temperature trends as revealed by alkenone and Mg/Ca paleothermometry. Quaternary Science Reviews 29, 9891004.CrossRefGoogle Scholar
Lind, S., Ingvaldsen, R., 2012. Variability and impacts of Atlantic Water entering the Barents Sea from the north. Deep Sea Research Part I 62, 7088.CrossRefGoogle Scholar
Locarnini, R.A., Mishonov, A.V., Antonov, J.I., Boyer, T.P., Garcia, H.E., Baranova, O.K., Zweng, M.M., et al. , 2013. In: Levitus, S., Mishonov, A. (Eds.), World Ocean Atlas 2013. Volume 1: Temperature. S. Levitus, Ed., A. Mishonov Technical Ed.; NOAA Atlas NESDIS 73, 40 pp.Google Scholar
Loeng, H., 1991. Features of the physical oceanographic conditions of the Barents Sea. Polar Research 10, 518.CrossRefGoogle Scholar
Lubinski, D.J., Korsun, S., Polyak, L., Forman, S.L., Lehman, S.J., Herlihy, F.A., Miller, G.H., 1996. The last deglaciation of the Franz Victoria Trough, northern Barents Sea. Boreas 25, 89100.CrossRefGoogle Scholar
Lubinski, D.J., Polyak, L., Forman, S.L., 2001. Freshwater and Atlantic water inflows to the deep northern Barents and Kara seas since ca 13 14C ka: Foraminifera and stable isotopes. Quaternary Science Reviews 20, 18511879.CrossRefGoogle Scholar
Mangerud, J., Bolstad, M., Elgersma, A., Helliksen, D., Landvik, J.Y., Lønne, I., Lycke, A.K., Salvigsen, O., Sandahl, T., Svendsen, J.I., 1992. The last glacial maximum on western Svalbard. Quaternary Research 38, 131.CrossRefGoogle Scholar
Mangerud, J., Bondevik, S., Gulliksen, S., Hufthammer, A.K., Høisæter, T., 2006. Marine 14 C reservoir ages for 19th century whales and molluscs from the North Atlantic. Quaternary Science Reviews 25, 32283245.CrossRefGoogle Scholar
Matthiessen, J., de Vernal, A., Head, M., Okolodkov, Y., Zonneveld, K., Harland, R., 2005. Modern organic-walled dinoflagellate cysts in Arctic marine environments and their (paleo-) environmental significance. Paläontologische Zeitschrift 79/1, 351.CrossRefGoogle Scholar
McAndrews, J.H., Berti, A.A., Norris, G., 1973. Key to the Quaternary Pollen and Spores of the Great Lakes Region. Royal Ontario Museum Life Sciences Miscellaneous Publication. University of Toronto Press, Toronto.CrossRefGoogle Scholar
Mertens, K.N., Verhoeven, K., Verleye, T., Louwye, S., Amorim, A., Ribeiro, S., Deaf, A.S. et al. , 2009. Determining the absolute abundance of dinoflagellate cysts in recent marine sediments: The Lycopodium marker-grain method put to the test. Review of Palaeobotany and Palynology 157, 238252.CrossRefGoogle Scholar
Müller, J., Werner, K., Stein, R., Fahl, K., Moros, M., Jansen, E., 2012. Holocene cooling culminates in sea ice oscillations in Fram Strait. Quaternary Science Reviews 47, 114.CrossRefGoogle Scholar
Murdmaa, I., Ivanova, E., 2017. Deglaciation of the Late Weichselian Barents Sea ice sheet. In: Boone, M. (Ed.), Deglaciation Processes, Causes and Consequences. Nova Science Publishers, Hauppauge NY, pp. 141171.Google Scholar
Peltier, W.R., Fairbank, R.G., 2006. Global glacial ice volume and Last Glacial Maximum duration from an extended Barbados sea level record. Quaternary Science Reviews 25, 33223337.CrossRefGoogle Scholar
Pfirman, S.L., Bauch, D., Gammelsrød, T., 1994. The northern Barents Sea: water mass distribution and modification. In: Johannessen, O.M., Muench, R.D., Overland, J.E. (Eds.), The Polar Oceans and Their Role in Shaping the Global Environment: The Nansen Centennial Volume. Geophysical Monograph Series, Vol. 85. American Geophysical Union, Washington, DC, pp. 7794.Google Scholar
Polyak, L., Mikhailov, V., 1996. Post-glacial environments of the southeastern Barents Sea: foraminiferal evidence. In: Andrews, J.T., Austin, W.E.N., Bergsten, H., Jennings, A.E. (Eds.), Late Quaternary Paleoceanography of the North Atlantic Margins. Geological Society, London, Special Publication No. 111., pp. 323337.Google Scholar
Potvin, E., Rochon, A., Lovejoy, C., 2013. Cyst–theca relationship of the arctic dinoflagellate cyst Islandinium minutum (Dinophyceae) and phylogenetic position based on SSU rDNA and LSU rDNA. Journal of Phycology 49, 848866.CrossRefGoogle ScholarPubMed
Radi, T., Bonnet, S., Cormier, M.-A., de Vernal, A., Durantou, L., Faubert, E., Head, M.J., Henry, M., Pospelova, V., Rochon, A., and Van Nieuwenhove, N., 2013. Operational taxonomy and (paleo-)autecology of round, brown, spiny dinoflagellate cysts from the Quaternary of high northern latitudes. Marine Micropaleontology 98, 4157.CrossRefGoogle Scholar
Radi, T., de Vernal, A., 2008. Dinocysts as proxy of primary productivity in mid–high latitudes of the Northern Hemisphere. Marine Micropaleontology 68, 84114.CrossRefGoogle Scholar
Raghavan, M., DeGiorgio, M., Albrechtsen, A., Moltke, I., Skoglund, P., Korneliussen, T.S., Grønnow, B. et al. , 2014. The genetic prehistory of the New World Arctic. Science 345, 1255832.CrossRefGoogle ScholarPubMed
Rasmussen, T.L., Bäckström, D., Heinemeier, J., Klitgaard-Kristensen, D., Knutz, P.C., Kuijpers, A., Lassen, S., Thomsen, E., Troelstra, S.R., van Weering, T.C.E., 2002. The Faroe–Shetland Gateway: Late Quaternary water mass exchange between the Nordic Seas and the northeastern Atlantic. Marine Geology 188, 165192.CrossRefGoogle Scholar
Rasmussen, T.L., Thomsen, E., Skirbekk, K., Slubowska-Woldengen, M., Klitgaard Kristensen, D., Koç, N., 2014. Spatial and temporal distribution of Holocene temperature maxima in the northern Nordic seas: Interplay of Atlantic-, Arctic- and polar water masses. Quaternary Science Reviews 92, 280291.CrossRefGoogle Scholar
Rasmussen, T.L., Thomsen, E., Ślubowska, M.A., Jessen, S., Solheim, A., Koç, N., 2007. Paleoceanographic evolution of the SW Svalbard margin (76°N) since 20,000 14C yr BP. Quaternary Research 67, 100114.CrossRefGoogle Scholar
Reimer, P.J., Bard, E., Bayliss, A., Beck, J.W., Blackwell, P.G., Bronk Ramsey, C., Buck, C.E., et al. , 2013. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. Radiocarbon 55, 18691887.CrossRefGoogle Scholar
Rigual-Hernández, A.S., Colmenero-Hidalgo, E., Martrat, B., Bárcena, M.A., de Vernal, A., Sierro, F.J., Flores, J.A., Grimalt, J.O., Henry, M., Lucchi, R.G., 2017. Svalbard ice-sheet decay after the Last Glacial Maximum: new insights from micropalaeontological and organic biomarker paleoceanographical reconstructions. Palaeogeography, Palaeoclimatology, Palaeoecology 465, 225236.CrossRefGoogle Scholar
Risebrobakken, B., Dokken, T., Smedsrud, L.H., Andersson, C., Jansen, E., Moros, M., Ivanova, E.V., 2011. Early Holocene temperature variability in the Nordic Seas: The role of oceanic heat advection versus changes in orbital forcing. Paleoceanography 26, PA4206. http://dx.doi.org/10.1029/2011PA002117.CrossRefGoogle Scholar
Risebrobakken, B., Moros, M., Ivanova, E., Chistyakova, N., Rosenberg, R., 2010. Climate and oceanographic variability in the SW Barents Sea during the Holocene. The Holocene 20, 609621.CrossRefGoogle Scholar
Rochon, A., de Vernal, A., Turon, J.-L., Matthiessen, J., Head, M.J., 1999. Distribution of Dinoflagellate Cyst Assemblages in Surface Sediments from the North Atlantic Ocean and Adjacent Basins and Quantitative Reconstruction of Sea Surface Parameters. American Association of Stratigraphic Palynologists Contribution Series, No. 35. American Association of Stratigraphic Palynologists Foundation, Dallas.Google Scholar
Rudels, B., Anderson, L.G., Jones, E.P., 1996. Formation and evolution of the surface mixed layer and halocline of the Arctic Ocean. Journal of Geophysical Research: Oceans 101, 88078821.CrossRefGoogle Scholar
Sætre, R., Ljøen, R. 1972. The Norwegian Coastal Current. Proceedings from the First International Conference on Port and Ocean Engineering under Arctic Conditions 1, 514535.Google Scholar
Sarnthein, M., Kreveld, S., Erlenkeuser, H., Grootes, P., Kucera, M., Pflaumann, U., Schulz, M., 2003. Centennial-to-millennial-scale periodicities of Holocene climate and sediment injections off the western Barents shelf, 75 N. Boreas 32, 447461.CrossRefGoogle Scholar
Schauer, U., Loeng, H., Rudels, B., Ozhigin, V.K., Dieck, W., 2002. Atlantic Water flow through the Barents and Kara Seas. Deep-Sea Research I 49, 22812298.CrossRefGoogle Scholar
Sejrup, H.-P., McKay, N.P., Kaufman, D.S, Geirsdottir, A., de Vernal, A., Renssen, H., Husum, K., Jennings, A., Andrews, J.T., 2016. North Atlantic-Fennoscandian Holocene climate trends and Mechanisms. Quaternary Science Reviews 147, 365378.CrossRefGoogle Scholar
Ślubowska, M.A., Koç, N., Rasmussen, T.L., Klitgaard-Kristensen, D., 2005. Changes in the flow of Atlantic water into the Arctic Ocean since the last deglaciation: evidence from the northern Svalbard continental margin, 80 N. Paleoceanography 20. http://dx.doir.org/10.1029/2005PA001141.CrossRefGoogle Scholar
Ślubowska-Woldengen, M., Koç, N., Rasmussen, T.L., Klitgaard-Kristensen, D., Hald, M., Jennings, A.E., 2008. Time-slice reconstructions of ocean circulation changes on the continental shelf in the Nordic and Barents Seas during the last 16,000 cal yr BP. Quaternary Science Reviews 27, 14761492.CrossRefGoogle Scholar
Ślubowska-Woldengen, M., Rasmussen, T.L., Koç, N., Klitgaard-Kristensen, D., Nilsen, F., Solheim, A., 2007. Advection of Atlantic Water to the western and northern Svalbard shelf since 17,500 cal yr BP. Quaternary Science Reviews 26, 463478.CrossRefGoogle Scholar
Smedsrud, L.H., Esau, I., Ingvaldsen, R.B., Eldevik, T., Haugan, P.M., Li, C., Lien, V.S. et al. , 2013. The role of the Barents Sea in the Arctic climate system. Reviews of Geophysics 51, 415449.CrossRefGoogle Scholar
Svendsen, J.I., Mangerud, J., 1997. Holocene glacial and climatic variations on Spitsbergen, Svalbard. The Holocene 7, 4557.CrossRefGoogle Scholar
Telford, R., Birks, H., 2009. Evaluation of transfer functions in spatially structured environments. Quaternary Science Reviews 28, 13091316.CrossRefGoogle Scholar
Telford, R.J., 2006. Limitations of dinoflagellate cyst transfer functions. Quaternary Science Reviews 25, 13751382.CrossRefGoogle Scholar
ter Braak, C.J.F., Smilauer, P., 2012. Canoco Reference Manual and User's Guide: Software for Ordination (Version 5). Microcomputer Power, Ithaca.Google Scholar
Van Nieuwenhove, N., Baumann, A., Matthiessen, J., Bonnet, S., de Vernal, A., 2016. Sea surface conditions in the southern Nordic Seas during the Holocene based on dinoflagellate cyst assemblages. The Holocene 26, 722735.CrossRefGoogle Scholar
Werner, K., Spielhagen, R.F., Bauch, D., Hass, H.C., Kandiano, E., 2013. Atlantic Water advection versus sea-ice advances in the eastern Fram Strait during the last 9 ka: Multiproxy evidence for a two-phase Holocene. Paleoceanography 28, 283295.CrossRefGoogle Scholar
Zonneveld, K.A.F., Marret, F., Versteegh, G.J.M., Bonnet, S., Bouimetarhan, I., Crouch, E., de Vernal, A. et al. , 2013. Atlas of modern dinoflagellate cyst distribution based on 2405 datapoints. Review of Palaeobotany and Palynology 191, 1197.CrossRefGoogle Scholar
Zweng, M.M, Reagan, J.R., Antonov, J.I., Locarnini, R.A., Mishonov, A.V., Boyer, T.P., Garcia, H.E., et al. , 2013. World Ocean Atlas 2013. Volume 2: Salinity S. Levitus, Ed., A Mishonov, Technical Ed. NOAA Atlas NESDIS 74, 39 p.Google Scholar
Figure 0

Figure 1. Maps showing the location of core S2528 and modern surface currents in the Barents Sea. The red arrows illustrate the circulation pattern of Atlantic waters at the surface, including the Norwegian Atlantic Current (NwAC) and the Norwegian Atlantic Slope Current (NwASC), which divides into the West Spitsbergen Current (WSC) and the North Cape Current (NCC). North of Svalbard, the WSC becomes the subsurface current of the Svalbard Branch (SB) and the Yermak Branch (YB). The blue arrows show the Arctic East Spitsbergen Current (ESC) and East Greenland Current (EGC). The green arrows correspond to coastal currents. The polar front position is represented by the white line based on Loeng (1991).

Figure 1

Figure 2. Map of mean summer (a) sea-surface temperature and (b) salinity over the Barents Sea between 1955 and 2012, after the World Ocean Atlas (Locarinini et al., 2013; Zweng et al., 2013). The white dashed lines separate the North Atlantic waters (NAW), the Arctic waters (ArcW) and the Coastal waters (CW), as defined by Loeng (1991).

Figure 2

Table 1. 14C dates of core S2528 and corresponding calibrated ages (Ivanova et al. [2019] with additions). Laboratory codes: Poz, Poznań Radiocarbon Laboratory; OS, NOSAMS facility at the Woods Hole Oceanographic Institution; CAMS, Lawrence Livermore National Laboratory.

Figure 3

Figure 3. Age versus depth relationship and sedimentation rates of core S2528 based on 14C dates and the use of an approach combining the Clam and Bacon softwares (Blaauw, 2010; Blaauw and Christen, 2011). Numbers next to polygons are the calibrated ages (yr BP) and the one in italics is the added age from interpolation with Clam (see section Chronology of the core).

Figure 4

Figure 4. Core lithological units, concentrations of dinocysts, pollen and spores, organic foraminifer linings, and reworked palynomorphs expressed in number of specimens/cm3 of sediment, and percentage of the coarse fraction (> 106 μm) in core S2528. The coarse fraction is expressed on a logarithmic scale due to large amplitude variations throughout the core.

Figure 5

Figure 5. Relative abundance (percentages) of the main dinocyst taxa in core S2528, concentration of phototrophic and heterotrophic dinocyst taxa, and results of the principal component (PC) analyses. The PC1 (black line) and PC2 (gray dashed line) axes respectively explain 45.2 and 24.5%, respectively, of the variance. The gray bands correspond to the main zones as determined from the assemblages and the PC analyses. The ages of the transitions were estimated from the age versus depth relationship shown in Figure 3.

Figure 6

Figure 6. Sea-surface reconstruction from the application of MAT to the dinocyst assemblages of core S2528. Sea-surface temperature and salinity in summer are represented by red curves, and sea-ice cover and annual productivity are represented by purple and green lines, respectively. The lighter areas correspond to maximum and minimum possible temperatures as calculated from the set of five modern analogues. The gray dashed lines correspond to smoothed values according to a five-points window. The horizontal rectangles at the top of the graphs correspond to the mean ± 1 SD of summer temperature and salinity (World Ocean Atlas) and sea-ice cover (National Snow and Ice Data Center data). For productivity, no modern direct data are available for statistics. The reliability index, which is calculated based on counts and the statistical distance between the fossil assemblage and modern analogues, indicates high to very high reliability for most samples. The horizontal dashed lines represent the transitions discussed in the text, with tentative chronology estimated from the age versus depth relationship of Figure 3.

Supplementary material: PDF

Brice et al. supplementary material

Brice et al. supplementary material 1

Download Brice et al. supplementary material(PDF)
PDF 108.8 KB
Supplementary material: File

Brice et al. supplementary material

Brice et al. supplementary material 2

Download Brice et al. supplementary material(File)
File 11.7 KB
Supplementary material: PDF

Brice et al. supplementary material

Brice et al. supplementary material 3

Download Brice et al. supplementary material(PDF)
PDF 65.7 KB
Supplementary material: PDF

Brice et al. supplementary material

Brice et al. supplementary material 4

Download Brice et al. supplementary material(PDF)
PDF 97.2 KB