Abstract
In high northern latitudes, the Middle to Late-Holocene was a time of orbitally-induced atmospheric cooling. This led to increased sea-ice production in the Arctic Ocean and its export southward, a decrease in sea surface temperatures (SST), and glacier advances at least since 5–4 ka BP. However, the response of the ocean-climate system to decreasing insolation was not uniform. Our research shows that the sea-ice cover in the northwestern Barents Sea experienced a late response to Neoglacial cooling. We analyzed dinoflagellate cyst assemblages from a sediment core from Storfjordrenna, south of Svalbard. We found that the area experienced ice-free conditions throughout most of the Mid- and Late-Holocene. It was only after 2.3 ka BP that the study site became covered with winter drift ice and primary productivity decreased subsequently. Other regional data support the decrease in SST, the expansion of the sea-ice cover, and the deterioration of the environmental conditions around that time. Our findings indicate that the sea-ice cover in the northwestern Barents Sea required a significant amount of time to respond to the general cooling trend in the region. These results have important implications for present-day environmental changes. Even if the current warming trend is revoked in the future, the observed sea-ice loss in the Barents Sea may be incredibly challenging to reverse.
Keywords
Introduction
The Middle to Late-Holocene in high northern latitudes was a time of decreasing temperatures (McKay et al., 2018) caused by declining boreal summer insolation (Laskar et al., 2004) and referred to as Neoglacial cooling (Wanner, 2021). Around 5 ka BP the modern sea level was reached and the postglacial flooding of the Laptev Sea shelves was finalized (Bauch et al., 2001), allowing the Arctic sea-ice production, which predominantly takes place on the shallow Arctic shelf areas, to reach its modern magnitude (Werner et al., 2013). This resulted in the onset of modern-like conditions with perennial sea ice in the Arctic Ocean (Cronin et al., 2010). Enhanced Mid-Holocene Siberian river runoff caused an eastward shift of Transpolar Drift and increased sea-ice export through the Fram Strait (Dyke et al., 1997; Prange and Lohmann, 2003). Data from the Greenland Ice Core Project ice core show the onset of atmospheric cooling between 5 and 4 ka BP (Dahl-Jensen et al., 1998). Although in the northern North Atlantic region there is no compelling evidence for any significant and widespread climatic anomaly (Bradley and Bakke, 2019) associated with the 4.2 ka BP event (Renssen, 2022) that marks the onset of the Late-Holocene (Walker et al., 2019), a compilation of glacier records indicates a significant climatic transition referred to as “the Holocene Turnover” occurred ~4 ka BP. It represents a dynamical adjustment that subsequently resulted in the establishment of a new climate regime or mode rather than being a multidecadal or centennial deviation from mean conditions (Paasche et al., 2004; Paasche and Bakke, 2009). On Svalbard, Late-Holocene glacier re-advances started around 4 ka BP (Farnsworth et al., 2020). Tidewater glaciers, which were largely absent in Svalbard during the early and Mid-Holocene, reappeared 3–4 ka BP (Jang et al., 2023; Svendsen and Mangerud, 1997). Similarly, in northern Norway, glaciers reappeared ~4 ka BP (Bakke et al., 2005). This coincides with the onset of sea surface temperature (SST) decrease in the northwestern Barents Sea continental slope (Rigual-Hernández et al., 2017; Risebrobakken et al., 2010). In the subpolar North Atlantic, pronounced SST cooling was observed between 4 and 2 ka BP (Orme et al., 2018). A significant sea-ice advance accompanied by distinct sea-ice fluctuations occurred in the eastern Fram Strait after 3 ka BP (Müller et al., 2012). By 2 ka BP, the SST in the Norwegian Sea had already decreased to its Holocene low (Andersen et al., 2004). However, the response of the ocean environment to decreasing boreal summer insolation was not uniform (Andersen et al., 2004, 2004; Wanner, 2021). Around 2 ka BP, the North Atlantic Oscillation, the dominant mode of atmospheric variability at mid-latitudes in the North Atlantic region, changed from variable, intermittently negative to generally positive conditions (Olsen et al., 2012), indicating stronger westerlies. This caused increased Atlantic Water (AW) advection into the Nordic Seas (e.g. Giraudeau et al., 2010; Spielhagen et al., 2011; Telesiński et al., 2014, 2015; Werner et al., 2013). As a result, ocean (e.g. Andersen et al., 2004; Sarnthein et al., 2003), as well as atmospheric (e.g. Johnsen et al., 2001; McDermott et al., 2001) warming, occurred. In Europe and the North Atlantic region, it is recognized as the Roman Warm Period (e.g. Bianchi and McCave, 1999; Matul et al., 2018; Wang et al., 2012).
Here we reconstruct paleoenvironmental changes in the northwestern Barents Sea during the Late-Holocene to estimate how fast the sea-ice cover reacted to changes in the atmosphere and the ocean. We analyze proxy data from a marine sediment core from Storfjordrenna, south of Svalbard, including dinoflagellate cyst (dinocyst) assemblages, as well as previously published XRF data, stable isotope, and alkenone-based SST records. Based on the abundance of indicator species of dinocysts, we reconstruct the reappearance of winter drift ice, as well as the changing influence of AW on the study site. We also compare our data with biomarker-based reconstruction of sea ice from a core from the Olga Basin, northern Barents Sea (Berben et al., 2017).
Oceanographic setting
The Barents Sea is an Arctic shelf sea between the Nansen Basin of the Arctic Ocean to the north, Novaya Zemlya to the east, Scandinavia to the south, the Norwegian Sea to the west, and the Svalbard archipelago to the north-west (Figure 1). The Barents Sea is influenced by several water masses. For this reason, strong environmental gradients can be observed here, making it a great area for studying paleoceanographic changes (e.g. Berben et al., 2017; Knies et al., 2017; Łącka et al., 2015, 2019).

Schematic map showing present-day surface water circulation in the Barents Sea. Red arrows indicate Atlantic Water, light blue arrows – Polar/Arctic Water, and white dashed line – Arctic Front (AF). The location of core JM09-020 is marked with an asterisk. The location of core NP05-11-70GC (Berben et al., 2017) also discussed in the paper is marked with a dot. BIC – Bear Island Current, ESC – East Spitsbergen Current, NCaC – North Cape Current, NwAC – Norwegian Atlantic Current, WSC – West Spitsbergen Current.
The relatively warm and saline AW (T > 3°C, S > 35.0; Loeng, 1991) is carried northward by the Norwegian Atlantic Current (Hopkins, 1991). The current is divided into the West Spitsbergen Current (WSC) and the North Cape Current. The North Cape Current enters the Barents Sea directly from the south-west, through Bjørnøyrenna, while the WSC continues northward along the shelf break, encircles Svalbard (Manley, 1995), and enters the Barents Sea from the north as a subsurface current, through Franz Victoria Trough (Abrahamsen et al., 2006; Rudels et al., 2015). Subsequently, AW is advected southwestward into the Olga Basin, where it has been observed year-round (Abrahamsen et al., 2006). After mixing and heat loss, AW exits the Barents Sea and reaches the Arctic Ocean via the St. Anna Trough (e.g. Rudels et al., 2015; Schauer et al., 2002).
The Polar Water (PW) is brought from the Arctic Ocean into the Barents Sea through the Franz Victoria and St. Anna Troughs, via the East Spitsbergen Current and the Bear Island Current, respectively. Arctic Water (ArW) is formed when relatively warm AW mixes with cold, less saline, and ice-loaded PW (Hopkins, 1991). Hence, surface water in the north-eastern Barents Sea is dominated by ArW, characterized by reduced temperature and salinity, as well as seasonal sea ice conditions (Hopkins, 1991).
The main oceanographic features of the near-surface waters of the Barents Sea are the oceanic fronts (Pfirman et al., 1994). Defined as sharp gradients in terms of temperature, salinity, and sea ice, the Polar and Arctic fronts are the respective boundaries between PW/ArW and ArW/AW. The positions of the Polar and Arctic fronts are closely related to the overall sea ice conditions and, in particular, align with the average summer and winter sea ice margins, respectively (Vinje, 1977). Although sea ice advection from the Arctic Ocean occurs, sea ice within the Barents Sea is formed mainly locally during autumn and winter (Loeng, 1991). The southward extent of the oceanic fronts and the sea-ice conditions are regulated by the inflow of AW into the western Barents Sea (Årthun et al., 2012), though in the west the PF is topographically controlled and therefore rather stable (Lien et al., 2017). On the contrary, the north-eastern Barents Sea experiences large changes in seasonal sea-ice conditions (Sorteberg and Kvingedal, 2006; Vinje, 2001) with maximum sea-ice conditions during March/April and minimum occurring throughout August/September.
The interplay between water masses determines the position of the marginal ice zone (MIZ) (Divine and Dick, 2006), an area characterized by high surface productivity during the summer season (e.g. Smith and Sakshaug, 1990). Within the Barents Sea, enhanced primary production results from a peak algal bloom along the MIZ as sea ice retreats in late spring (Hebbeln and Wefer, 1991; Ramseier et al., 1999; Sakshaug, 2004). Additionally, AW advection contributes to longer productive seasons, compared to other Arctic areas (Wassmann, 2011). Consequently, the Barents Sea is one of the most productive areas of the Arctic seas (Wassmann, 2011; Wassmann et al., 2006).
Material and methods
Sediment gravity core JM09-020 was retrieved from Storfjordrenna, northwestern Barents Sea (Figure 1, 76°19′N, 19°42′E, 253 m water depth, Łącka et al., 2015) and has been successfully used to reconstruct paleoceanographic conditions in the area over the last 14 kyr (Łącka et al., 2015, 2019, 2020).
For dinocyst analysis, the core was sampled every 4–6 cm. Each 1-cm-thick slab of sediment was collected into a zip bag and stored at a temperature of −20°C. After thawing, 3–4 cm3 of well-mixed sediment was put in a polypropylene test tube, dried at >40°C, and weighed with an analytical balance. Samples were subsequently soaked with distilled water for 12 h, centrifuged at 3600 rpm for 6 min, and then processed using a standard palynological technique (e.g. Pospelova et al., 2005, 2010).
Marker grains of a known number of Lycopodium clavatum spores (e.g. Mertens et al., 2009, 2012b) were added to allow quantitative estimates of the absolute concentrations of dinocysts. At room temperature, about 7 ml of hydrochloric acid (HCl, 10%) was slowly added to samples to dissolve the L. clavatum spore tablets and remove carbonates. After 30 min samples were centrifuged and decanted. Subsequently, ~9 ml of distilled water was added and samples were centrifuged and decanted again. The procedure was repeated until the pH of the supernatant reached a neutral level. Afterward, the samples were wet-sieved through 125 and 15 µm mesh to remove fractions of sediment above and below the maximal and minimal size of dinocysts.
After sieving, centrifuging, and decanting, ~7 ml of room-temperature hydrofluoric acid (HF, 48%) was added to the sediment to remove silicate. Samples were left in a fume hood for 72 h, with regular digestion checking and stirring. After silicate dissolution, samples were once again centrifuged and decanted and ~7 ml of hydrochloric acid (HCl, 10%, at room temperature) was added. Samples were rinsed with distilled water as described above and sieved through a 15 µm mesh. Aliquots of a few drops of sample residue were placed on a glass slide and left for 24 h at room temperature to dry. Glycerine gel was used to mount a cover slide to the glass slide.
Approximately 300 dinocyst specimens (min 201, max 341) were counted from each sample. Dinocysts were identified to the lowest possible taxonomical level. The paleontological taxonomy system used throughout this paper follows Zonneveld (1997), Kunz-Pirrung (1998), Montresor et al. (1999), Rochon et al. (1999), Head et al. (2001), Pospelova and Head (2002), Moestrup et al. (2009), Mertens et al. (2013, 2015, 2012a), and Zonneveld and Pospelova (2015). Cysts with unknown taxonomic affinity were classified into one of four groups: unidentified 1 – round transparent cyst, unidentified 2 – spiny transparent cyst, RBC – round brown cyst, and SBC – spiny brown cyst. Cysts of Biecheleria cf. baltica are mostly very small (~5–10 µm) and were partly lost during sample preparation (sieving). Therefore, we excluded them from the total cyst concentrations statistical analyses. Furthermore, it cannot be excluded that some thin-walled transparent Impagidinium spp. cysts have been missed during the counting (Telesiński et al., 2023). Dinocyst fluxes [cysts cm−2 yr−1] were calculated using absolute dinocyst abundances [cysts g−1], sedimentation rates [cm yr−1], and dry bulk density [g cm−3].
The chronology of core JM09-020 was based on radiocarbon dating (Łącka et al., 2015). We recalibrated the AMS 14C dates using CALIB 14C age calibration software (rev 8.1.0; Stuiver and Reimer, 1993) and the Marine20 calibration curve (Heaton et al., 2020). A regional correction of ΔR = −53 ± 36 14C years was applied. This value was calculated with the Marine Reservoir Correction database (Reimer and Reimer, 2001) and the Marine20 curve (Heaton et al., 2020, 2022) using the same mollusk samples as those used by Mangerud et al. (2006) for Svalbard. The difference between the resulting and the original age model (Łącka et al., 2015) is less than 50 years within the Holocene, which is insignificant for the present study, allowing for a direct comparison with previous studies of the core (Łącka et al., 2015, 2019).
Results
Here we present only selected parameters of the dinocyst assemblage analysis that are important for the current study. Complete results can be found in the Supplemental Material. Dinocyst flux was low (<100 cysts cm−2 yr−1) in the earliest Holocene (Figure 2a). It increased to 100–400 cysts cm−2 yr−1 between 11 and 9 ka BP. Subsequently, the flux decreased again to <100 cysts cm−2 yr−1. After 4 ka BP, the flux increased gradually to reach a maximum of ~900 cysts cm−2 yr−1 at 2.2 ka BP and then decreased again to reach ~300 cysts cm−2 yr−1 at 1.3 ka BP. The dinocyst assemblage was generally dominated by heterotrophic species (Figure 2b). However, the percentage of autotrophic dinocysts gradually increased from <5% in the Early Holocene to a maximum of 51% at 2.3 ka BP. Subsequently, the percentage of autotrophic cysts decreased again to 21% at the end of the record. The abundance of Echinidinium karaense was relatively high in the Early Holocene, though it never exceeded a relative abundance of 5% or a flux of 5 cysts cm−2 yr−1 (Figure 2c). After 8 ka BP the species disappeared completely from the record and reappeared only around 2.1 ka BP, reaching a flux of around 4 cysts cm−2 yr−1 at the end of the record, though its relative abundances were lower (up to 1.3%) than in the Early Holocene. The abundance of Operculodinium centrocarpum s.l. was extremely low (<5 cysts cm−2 yr−1 and <2%, respectively) throughout the Early Holocene (Figure 2d). Starting from 7.5 ka BP, its relative abundance increased gradually to reach a maximum of 44% around 2.3 ka BP. The flux remained low until ~4 ka BP but later also increased to reach a maximum of 356 cysts cm−2 yr−1 around 2.2 ka BP. Subsequently, both the relative abundance and the flux decreased (to 17% and 58 cysts cm−2 yr−1, respectively) toward the end of the record. The relative abundance of Islandinium minutum (Figure 2e) in the earliest Holocene was high (10–20%), though its flux remained relatively low (<20 cysts cm−2 yr−1). Between 11 and 2.1 ka BP, both relative abundance and flux were low. Only around 2.1 ka BP both relative abundance and the flux of this species increased rapidly to approximately 20% and 40–100 cysts cm−2 yr−1, respectively, and remained high until the end of the record.

Dinocyst record of core JM09-020. (a) Dinocyst flux [cysts cm−2 yr−1]. (b) Relative abundance of autotrophic versus heterotrophic species. (c) Flux [cysts cm−2 yr−1] and relative abundance [%] of Echinidinium karaense. (d) Flux [cysts cm−2 yr−1] and relative abundance [%] of Operculoidinium centrocarpum s.l. (e) Flux [cysts cm−2 yr−1] and relative abundance [%] of Islandinium minutum.
Discussion
The dinocyst species Echnidinium karaense, together with cysts of Polarella glacialis, has recently been identified as a winter drift ice indicator in waters around Svalbard (Telesiński et al., 2023). As the latter species is virtually absent in core JM09-020 (Supplemental Material), Echinidinium karaense remains the only available dinocyst sea-ice indicator. It was present in the western Barents Sea over the Early Holocene but it disappeared around 8 ka BP (Figure 2c), indicating ice-free conditions. Its reappearance at around 2.1 ka BP, after almost 6 thousand years of absence, clearly indicates a return of sea-ice conditions comparable to those in the Early Holocene. This is further supported by other dinocyst data from the same core. The peak in autotrophic dinocyst abundance at 2.3 ka BP (Figure 2b), followed by a peak in total dinocyst flux (Figure 2a) shortly thereafter, indicates increased primary productivity, which might suggest that the core site was reached by the MIZ (e.g. Barber et al., 2015; Ramseier et al., 1999; Sakshaug, 2004). After ~2 ka BP the total and autotrophic dinocyst abundance decreased, suggesting deteriorating surface water conditions, possibly due to the thickening of the sea-ice cover. Similarly, the abundance of Operculoidinium centrocarpum s.l., a cosmopolitan species whose high abundances in high northern latitudes are associated with AW dominance (Grøsfjeld et al., 2009; Rochon et al., 1999; Telesiński et al., 2023) reached a maximum around 2.3 ka BP but decreased sharply shortly thereafter, though remained higher than in the first half of the Holocene (Figure 2d). Furthermore, the relative percentage of O. centrocarpum s.l. versus I. minutum, which was relatively high throughout most of the Late-Holocene (Figure 3c), decreased after 2.1 ka BP. The relative percentage of these two species may be used to indicate whether warm AW flows at the surface or as a subsurface water mass (Grøsfjeld et al., 2009). This suggests that until ~2.1 ka BP, AW remained at the surface in the northwestern Barents Sea, while it subducted below ArW thereafter.

Paleoceanographic proxies of Late-Holocene changes in the northwestern Barents Sea from cores JM09-020 (a–f) and NP05-11-70GC (g). (a) Dinocyst flux [cysts cm−2 yr−1]. (b) Flux [cysts cm−2 yr−1] and relative abundance [%] of Echinidinium karaense. (c) Relative [%] abundance of Operculoidinium centrocarpum s.l. and Islandinium minutum. (d) Alkenone-based sea-surface temperature reconstruction (Łącka et al., 2019). Thin line – raw data, thick line – 3 pt moving average. e) Stable carbon isotope ratios [‰ vs VPDB] of benthic foraminifera Elphidium clavatum (Łącka et al., 2015). (f) Ba/Ti elemental ratios obtained from XRF core scanning (Łącka et al., 2015). Thin line – raw data, thick line – 5 pt moving average. (g) PIIIIP25 index and spring sea-ice concentration (SpSIC) [%] calculated from it (Berben et al., 2017). The horizontal dashed line marks the 0.8 PIIIIP25 threshold, indicating >5% summer sea-ice concentration.
Additional evidence from core JM09-020 corroborates the dinocyst data. An SST reconstruction based on alkenones (Łącka et al., 2019) shows a clear cooling trend between 2.3 and 2 ka BP (Figure 3d), which could be attributed to the expansion of sea ice. Similarly, the stable carbon isotope values of benthic foraminifera (Łącka et al., 2015) indicate increased variability in environmental conditions after 2 ka BP on the sea bottom (Figure 3e), which could be linked to enhanced sea-ice cover, variable productivity at the sea surface, and the amount of organic matter reaching the sea floor. Further details are provided by XRF data (Łącka et al., 2015). The Ba/Ti ratio exhibits a stepwise decrease around 2.3 ka BP (Figure 3f). Since the Ba/Ti ratio is believed to be broadly proportional to the organic carbon content in sediment (Thomson et al., 2006), such a decrease could indicate declining productivity (Croudace et al., 2006).
Based on the available data, it is evident that the northwestern Barents Sea witnessed a period of maximum productivity around 2.3 ka BP. This was mainly due to the inflow of warm surface AW from the west and the migration of the MIZ from the east. However, after 2.3–2.1 ka BP, the study site was covered with sea ice, and the AW submerged below surface ArW, which resulted in a decrease in SST and productivity.
Our results are further confirmed by data from core NP05-11-70GC from the Olga Basin, east of Svalbard (Berben et al., 2017). The PIIIIP25 index combines concentrations of tri-unsaturated highly branched isoprenoid (HBI) lipid (HBI III), a phytoplankton-derived biomarker, with IP25, a sea-ice proxy, to investigate past sea-ice conditions more quantitatively (Belt et al., 2007; Berben et al., 2017; Müller et al., 2011). The Olga Basin record indicates a constant increase in sea-ice concentration over the Holocene (Berben et al., 2017). Around 2.8 ka BP, the PIIIIP25 index crossed the 0.8 threshold (Figure 3g), indicating >5% summer sea-ice concentration (Smik et al., 2016). Around 2.5 ka BP, spring sea-ice concentration in the Olga Basin, derived from the PIIIIP25 index (Berben et al., 2017), reached 70%. Finally, around 1.9 ka BP, another stepwise increase of the PIIIIP25 index (approximately equal to eight of the total Holocene increase) occurred. The data from the Olga Basin confirm that a strong environmental gradient characterized the Barents Sea also in the past. While in the northern part of the basin, a dense spring sea-ice concentration was reached already around 2.5 ka BP, in the western part winter drift ice only appeared around 2.1 ka BP (Figure 3b and g). Nevertheless, data from both records confirm that the sea-ice cover reacted slowly to the Neoglacial cooling. Similarly, oxygen isotope data from core NP05-71GC from south of Kvitøya (Klitgaard-Kristensen et al., 2013) suggest that only after c. 2.5 ka BP the northwestern Barents Sea experienced cooling and/or increased brine formation, most probably related to the sea-ice expansion.
All the presented data suggest that the sea-ice expansion in the north-western Barents Sea occurred around 2.5–2.1 ka BP. Such a late response of the sea-ice cover to Neoglacial cooling is surprising. Firstly, the Arctic sea-ice production on the Siberian shelves has reached its modern magnitude already around 5 ka BP (Bauch et al., 2001; Werner et al., 2013) causing perennial sea-ice cover in the Arctic Ocean (Cronin et al., 2010). The fact that the sea-ice cover in the Barents Sea did not respond to this increase can be explained by the minor influence of advected (as opposed to locally formed) sea ice on the basin’s ice budget during the Late-Holocene (Loeng, 1991). On the other hand, terrestrial data indicate that the glacier re-advance in Svalbard began as early as ~4 ka BP (Farnsworth et al., 2020; Jang et al., 2023; Svendsen and Mangerud, 1997), suggesting that the atmospheric cooling in the north-western Barents Sea region required for the glacier growth was already achieved at the beginning of the Late-Holocene.
The approximately 2 kyr delay of the sea-ice expansion relative to the onset of atmospheric cooling and glacier advance in the region indicates that sea-ice cover in the Barents Sea needed significant time to recover, even in favorable climatic conditions. This was probably caused by the strong influence of AW, whose intrusions into the Barents Sea were frequent during the Middle and Late-Holocene (e.g. Pawłowska et al., 2020; Risebrobakken et al., 2010). It is worth noting that in core JM09-020, the abundance of the AW-indicating species Operculoidinium centrocarpum s.l. reached its maximum only ~2.3 ka BP (Figure 3c) and SST as well as productivity remained high until that time (Figure 3d and f), suggesting that in the western Barents Sea, the influence of AW was still increasing well into the Late-Holocene, despite the ongoing expansion of the sea-ice cover in the northern and eastern parts of the basin (Berben et al., 2017).
In the western Barents Sea, the PF is currently mainly topographically controlled (Lien et al., 2017). However, during the warm middle Holocene, the PF most probably decoupled from the bottom topography, allowing AW to reach much farther to the northeast (e.g. Berben et al., 2017). As a result, when orbitally forced Neoglacial cooling began, time was needed to push surface AW out of the central part of the Barents Sea, whereas in the west AW could have even increased its inflow on the surface. Even when atmospheric cooling in the region was advanced enough ~4 ka BP to allow the advance of Svalbard glaciers, another ~2 kyr was needed for the sea surface to cool enough to allow the sea ice to expand into the western Barents Sea. On the other hand, open water in the vicinity of Svalbard must have been an important source of moisture that allowed the growth of the glaciers (e.g. Hebbeln et al., 1994).
Over the last decades, AW intrusions on the Barents Sea shelf have become increasingly common (Kujawa et al., 2021; Telesiński et al., 2023; Walczowski and Piechura, 2011) as a result of enhanced northward heat transfer by the North Atlantic Drift (e.g. Spielhagen et al., 2011; Walczowski and Piechura, 2007), a phenomenon referred to as “Atlantification” of the Barents Sea (e.g. Årthun et al., 2012; Tesi et al., 2021). The delayed response of the sea-ice cover in the Barents Sea to Late-Holocene cooling demonstrated in this study suggests that even if the ongoing global warming is reversed in the future, which in itself is a highly challenging task, many centuries might be required for the sea ice to recover to its preindustrial extent. Taking into account that shrinking sea ice is one of the main drivers of the Arctic amplification (Serreze and Francis, 2006) as it reduces surface albedo, leading to greater surface solar absorption, amplifying warming, and further melt (e.g. Curry et al., 1995; Thackeray and Hall, 2019), the currently observed rapid sea-ice loss (e.g. Overland and Wang, 2013) might be an incredibly slow and long process to reverse.
Summary and conclusions
Reconstructing the paleoceanographic evolution of the northwestern Barents Sea during the Late-Holocene has been made possible by analyzing dinocyst assemblage data from sediment core JM09-020 from Storfjordrenna, south of Svalbard. The dinocyst data has been supplemented by stable carbon isotope, alkenone-based SST, and XRF data that have been previously published (Łącka et al., 2015, 2019). Furthermore, the dinocyst data has been compared with biomarker-based data from core NP05-11-70GC from the Olga Basin, east of Svalbard (Berben et al., 2017).
Based on the data, it appears that despite the ongoing Neoglacial cooling that began around 5 ka BP in high northern latitudes, the northeastern Barents Sea experienced a period of maximum productivity around 2.3 ka BP. This was due to two factors: the dominance of warm AW on the surface, and the proximity of the MIZ. Only after 2.3–2.1 ka BP did winter drift ice begin to cover the northwestern Barents Sea, resulting in a decrease in SST and productivity due to the subduction of AW below ArW.
Our findings have important implications for the current and future environmental changes. The presented data show that the recovery of the sea ice in the Barents Sea is a slow process. Even if the ongoing global warming can be halted or even revoked in the future, the reversing of the present sea-ice loss in the Barents Sea may be an incredibly long process.
Footnotes
Acknowledgements
We would like to express our gratitude to an anonymous reviewer for their helpful feedback and valuable suggestions, which greatly improved the quality of the manuscript.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by grant no. 2020/39/B/ST10/01698 funded by the National Science Centre, Poland.
Supplemental Material
Supplementary Material, containing all the dinocyst data from core JM09-020, is available on Zenodo (https://doi.org/10.5281/zenodo.8322505, Telesiński et al., 2023).
