Holocene Paleoenvironmental Implications of Diatom, Non-Pollen Palynomorph, and Organic Carbon Records from the Kandalaksha Bay of the White Sea (European Arctic)

: Variations in sea surface conditions and sea level through the Holocene in the Kandalaksha Bay, the White Sea, were reconstructed based on the study of core sediments from the outer Kan-dalaksha Bay, using the modern analog technique applied to dinocysts in addition to diatoms, TOC, δ 13 C org , CaCO 3 , and grain size data. The chronostratigraphy of the core sediments was deﬁned from accelerator mass spectrometry 14 C dates on mollusk shells. The results indicated an increase in water depth in the outer Kandalaksha Bay and in the central Dvina Bay until the late Holocene. From about 9.5 to 7.5 cal kyr BP, the data suggested a general trend of increasing sea surface temperatures (up to 14 ◦ C), at least in areas with inﬂow of Atlantic waters. The last 2.5 kyr were characterized by increased freshwater runoff to the White Sea.


Introduction
The amplified response to the recent global climate change in the Arctic that has been widely documented relates to increased heat transfer from the North Atlantic, atmospheric and surface ocean temperature rise, and reduction of the sea ice cover [1][2][3][4][5][6].Therefore, it is particularly relevant to study past and ongoing changes in Arctic marine ecosystems and to reconstruct the variations with the aim to understand and quantify the natural environmental variability and to predict its further evolution under anthropogenic stresses.
The White Sea belongs to the Arctic Ocean as a small and isolated water body (Figure 1a).Its hydrological regime is primarily governed by the inflow of the Barents Sea waters, which are formed from both the Arctic water mass flowing from the north and the Atlantic waters entering the Arctic Ocean from the southwest via the North Atlantic Current (NAC).In the southwestern Barents Sea, the NAC is divided into two main branches: the West Spitsbergen Current turning north along the continental slope and the North Cape stream heading into the Barents Sea (Figure 1a,b).Therefore, the transformed Barents Sea waters entering the White Sea contain an oceanographic signal from the North Atlantic [7][8][9][10].
Despite several research projects conducted on White Sea sediment cores [20,[24][25][26][27], the chronostratigraphy of the area remained poorly defined due to low content in biogenic carbonate remains in Holocene sequences making it difficult to obtain the 14 C date.The rare dates obtained through accelerator mass spectrometry (AMS) 14 C measurements were published in our previous studies [21,28].
Here we present new paleoenvironmental data based on the analyses of the sediment core PS-6066, in which an age model based on new AMS 14 C dates was defined.The core PS-6066 was collected in the Kandalaksha Bay, the deepest estuary in the White Sea that reaches a maximum water depth of 343 m [29].Our aim was to reconstruct the sea surface salinity and temperature, seasonal sea ice cover, and other environmental parameters based on multiproxy records.These included micropaleontological tracers (diatoms and dinocysts), in addition to lithological and geochemical properties (grain size, CaCO3 and organic carbon concentration, and δ 13 C of organic carbon) of the sediment.Diatom and non-pollen palynomorph (NPP) assemblages are among the most useful micropaleontological groups for paleoenvironmental reconstructions in the White Sea region [20][21][22][30][31][32][33].Diatoms dominate in the phytoplankton communities [34][35][36][37][38], and dinoflagellates are very important components of the White Sea ecosystem.Our interpretations of fossil assemblages of diatoms and NPP in core sediments were based primarily on the previously established patterns of their distribution in the surface sediments and their relationship with hydrological and sedimentation processes [21,31,32,39,40].We used the modern analog technique [41,42] for the reconstruction of the salinity, temperature, and duration of the sea ice cover.

Modern Setting
The White Sea is the most enclosed part of the Eurasian sector of the Arctic Ocean and is almost entirely located south of the polar circle.The origin of the White Sea basin is associated with the development of a rift system that arose in the Middle-Late Riphean on the Early Precambrian consolidated basement [43].The configuration of the modern basin and the physiography of the surrounding areas have formed under the Neogene-Quaternary tectonic and Quaternary glacial activity [11,44].[20], PS-6050, PS-6042 by [21,22], and PL 96-112 by [23].
Despite several research projects conducted on White Sea sediment cores [20,[24][25][26][27], the chronostratigraphy of the area remained poorly defined due to low content in biogenic carbonate remains in Holocene sequences making it difficult to obtain the 14 C date.The rare dates obtained through accelerator mass spectrometry (AMS) 14 C measurements were published in our previous studies [21,28].
Here we present new paleoenvironmental data based on the analyses of the sediment core PS-6066, in which an age model based on new AMS 14 C dates was defined.The core PS-6066 was collected in the Kandalaksha Bay, the deepest estuary in the White Sea that reaches a maximum water depth of 343 m [29].Our aim was to reconstruct the sea surface salinity and temperature, seasonal sea ice cover, and other environmental parameters based on multiproxy records.These included micropaleontological tracers (diatoms and dinocysts), in addition to lithological and geochemical properties (grain size, CaCO 3 and organic carbon concentration, and δ 13 C of organic carbon) of the sediment.Diatom and nonpollen palynomorph (NPP) assemblages are among the most useful micropaleontological groups for paleoenvironmental reconstructions in the White Sea region [20][21][22][30][31][32][33].Diatoms dominate in the phytoplankton communities [34][35][36][37][38], and dinoflagellates are very important components of the White Sea ecosystem.Our interpretations of fossil assemblages of diatoms and NPP in core sediments were based primarily on the previously established patterns of their distribution in the surface sediments and their relationship with hydrological and sedimentation processes [21,31,32,39,40].We used the modern analog technique [41,42] for the reconstruction of the salinity, temperature, and duration of the sea ice cover.

Modern Setting
The White Sea is the most enclosed part of the Eurasian sector of the Arctic Ocean and is almost entirely located south of the polar circle.The origin of the White Sea basin is associated with the development of a rift system that arose in the Middle-Late Riphean on the Early Precambrian consolidated basement [43].The configuration of the modern basin and the physiography of the surrounding areas have formed under the Neogene-Quaternary tectonic and Quaternary glacial activity [11,44].
The average seawater depth of the White Sea is about 67 m.The maximum depths are recorded in the central depression named the Basin (down to 350 m water depth) and in the Kandalaksha Bay, where it reaches 343 m [29].With respect to the structure and geomorphology, the following major structural elements are distinguished: the Voronka, the Gorlo, the Basin, and the bays of Onega, Dvina, Mezensky, and Kandalaksha, into which large rivers of the same name flow (Figure 1b).These bays differ in the freshwater supply, amplitude of tides, salinity gradients, sea ice conditions, and biota [45][46][47][48].
The present-day hydrology of the White Sea is mostly under the influence of the freshwater runoff, water exchange with the Barents Sea, and strong tidal currents.Freshwater runoff is carried out mainly by the largest rivers, the Dvina, Mezen, Onega, and others, the total flow of which exceeds 225 km 3 /year, which defines the regularities and particular features of the hydrological and biogeochemical processes within its area [7,29,46].Tides in the White Sea reach their maximal heights of up to 7.7 m in the Mezensky Bay and up to 1.6 m in the Dvina Bay, while they do not exceed 2.5 m in the Kandalaksha Bay.Under the conditions of extensive shallow waters in bays, tides have a profound impact on the mixing of water masses in the White Sea.
Water masses entering the White Sea from the Barents Sea have an Atlantic origin [7,10].They are characterized by a consistently high salinity of 34-35 psu and a relatively high temperature (Figure 2a,b) as compared with other Arctic seas.The surface water temperatures in the White Sea range from +7 to +15 • C in summer and fall down to a freezing point below −1.6 • C in winter [29].There is a general cyclonic circulation in the sea, homogenizing the water masses.The hydrographic system is characterized by uniform currents due to the continental runoff, concentrated in the upstream part of the bays.In the outer Dvina Bay, the cold deep-sea layers are much closer to the surface than in other parts of the basin.The temperature of 0 • C is observed here in summer only 12-15 m below the water surface.The existence of this unique "hydrological pole of cold" [49] is due to the cyclonic circulation of surface waters with deep-water upwelling in the center.The average seawater depth of the White Sea is about 67 m.The maximum depths are recorded in the central depression named the Basin (down to 350 m water depth) and in the Kandalaksha Bay, where it reaches 343 m [29].With respect to the structure and geomorphology, the following major structural elements are distinguished: the Voronka, the Gorlo, the Basin, and the bays of Onega, Dvina, Mezensky, and Kandalaksha, into which large rivers of the same name flow (Figure 1b).These bays differ in the freshwater supply, amplitude of tides, salinity gradients, sea ice conditions, and biota [45][46][47][48].
The present-day hydrology of the White Sea is mostly under the influence of the freshwater runoff, water exchange with the Barents Sea, and strong tidal currents.Freshwater runoff is carried out mainly by the largest rivers, the Dvina, Mezen, Onega, and others, the total flow of which exceeds 225 km 3 /year, which defines the regularities and particular features of the hydrological and biogeochemical processes within its area [7,29,46].Tides in the White Sea reach their maximal heights of up to 7.7 m in the Mezensky Bay and up to 1.6 m in the Dvina Bay, while they do not exceed 2.5 m in the Kandalaksha Bay.Under the conditions of extensive shallow waters in bays, tides have a profound impact on the mixing of water masses in the White Sea.
Water masses entering the White Sea from the Barents Sea have an Atlantic origin [7,10].They are characterized by a consistently high salinity of 34-35 psu and a relatively high temperature (Figure 2 a,b) as compared with other Arctic seas.The surface water temperatures in the White Sea range from +7 to +15 °C in summer and fall down to a freezing point below −1.6 °C in winter [29].There is a general cyclonic circulation in the sea, homogenizing the water masses.The hydrographic system is characterized by uniform currents due to the continental runoff, concentrated in the upstream part of the bays.In the outer Dvina Bay, the cold deep-sea layers are much closer to the surface than in other parts of the basin.The temperature of 0 °C is observed here in summer only 12-15 m below the water surface.The existence of this unique "hydrological pole of cold" [49] is due to the cyclonic circulation of surface waters with deep-water upwelling in the center.

Core Location and Processing
The gravity core PS-6066 (65°59.691′'N, 35°32,387′ E; water depth, 266 m) was collected in the Kandalaksha Bay (Figure 1a,b), the NW White Sea, during the R/V Professor

Core Location and Processing
The gravity core PS-6066 (65 • 59.691 N, 35 • 32,387 E; water depth, 266 m) was collected in the Kandalaksha Bay (Figure 1a,b), the NW White Sea, during the R/V Professor Shtokman Cruise 80 in 2006.The core was sampled at 1 cm intervals, and the sediment samples were placed in plastic boxes and stored in a cold room (<5 • C) until freeze-drying.

Organic Matter Analysis
The total organic carbon (TOC) content was determined by the automatic coulometric method using an AN7529M Carbon Analyzer M (GZIP, Gomel, Republic of Belarus) that yielded results with an accuracy of ±0.1%.The CaCO 3 content was calculated from inorganic carbon with a coefficient of 8.33 according to molar ratio.The results were given in weight percent (wt.%).
The δ 13 C org was determined from carbon dioxide obtained through the burning of the samples previously treated with hydrochloric acid to remove carbonate compounds.The δ 13 C was measured using a Delta Plus mass spectrometer (Thermo Electron Corporation, Germany) at the Winogradov Institute of Microbiology, Russian Academy of Sciences, Federal Research Centre "Fundamentals of Biotechnology" of the Russian Academy of Sciences.The reproducibility was better than ±0.1‰ for δ 13 C.The results were presented as δ 13 C versus the international standard Vienna Pee Dee Belemnite (VPDB).

Chronology and Sedimentation Rates
Two AMS 14 C dates on bivalves and brachiopods were obtained from the core PS-6066 (Table 1) and calibrated to calendar ages using the Bacon software [55] and the Marine 20 calibration data set that included corrections for the air-sea reservoir difference [56].We further applied a delta R of −132 ± 48.0 yr calculated using the method of Bevington [57] and the regional data reported in Zaretskaya et al. [58].In the White Sea, Zaretskaya et al. [58] obtained marine reservoir values ranging from 435 to 535 14 C years for mollusks inhabiting depths of 20-30 m, thus practically leveled out by the correction for isotopic fractionation.The authors, however, mentioned that some mollusks inhabiting deeper than 30-40 m settings may record much larger reservoir effects [58].It should be noted that the White Sea sediments do not contain much material for dating, most of the foraminifer shells being agglutinated.Because carbonate material could not be recovered in the 39−211 cm interval from core PS-6066, we tried to constrain the chronology from correlations with the lithology of the previously studied core PS-6050 [22].On this basis, we set an interpolated age (see also Section 4.1.Lithology) that was added in the Bacon input file to generate a depth-age model as illustrated in Figure 3.
The ages are given hereafter in cal kyr BP (thousands of calendar years before 1950 CE).The estimated age-depth relationship suggested that average sedimentation rates in the Kandalaksha Bay were about 18-20 cm yr −1 , with maximum values of up to 27 cm yr −1 for the last ~1.5 cal kyr BP (Figure 3).Thus, the sedimentation rates seemed mostly uniform, which was consistent with continuous marine sedimentation during the Holocene.
The age model allowed us to make ecostratigraphical correlations based on diatom ecozones and pollen and spore assemblages of other cores from the White Sea area [20,21,25,28].The uncertainties in the absolute ages are larger than 200 years, as is often the case in marine sequences of subarctic settings [59,60].Figure 3. Age-depth model for core PS-6066 constrained from the Bacon software [55] to calibrated accelerator mass spectrometry 14 C dates from Table 1.Numbers are the calibrated ages (yr BP), and the one in italics is the added age from interpolation with PS-6050.
We also used the classification of phytogeographical groups developed by Beklemishev and Semina [86], with some clarifications of the species affiliation [34,87].The diatom species were further divided into the following phytogeographical groups: cosmopoli-tan, boreal, arcto-boreal, and bipolar.Among the neritic and panthalassic species, we distinguished relatively warm-water species such as Coscinodiscus radiatus, C. perforatus, Shionodiscus oestrupii, Actinoptychus senarius, and others, which were all indicators of Atlantic Water advection into the Arctic seas [84][85][86].Particular attention was paid to the sea ice species; Attheya septentrionalis, Melosira arctica, and Nitzschia frigida, among others, are associated with sea ice at least for a part of their life cycle [85,88,89].We also noted the presence of ice-neritic species, which developed in the phytoplankton of the marginal ice zone, at low surface water temperatures close to or below 0 • C, such as Porosira glacialis, Thalassiosira nordenskioeldii, Bacterosira bathyomphala, Rhizosolenia hebetata f. hebetata, and Chaetoceros furcellatus [85,[89][90][91][92].
The freshwater diatom genera Aulacoseira and Cyclotella and green algae appeared on the shelf with a riverine outflow, and most of them were deposited within the coagulationsorption stage of the marginal filter of rivers [32, 73,93].
According to the change in total diatom concentration in sediments, taxonomic diversity, and distribution of the main ecological groups, as well as the dominant and subdominant taxa, five diatom ecozones (DZ) were distinguished.

Non-Pollen Palynomorphs, "Pollen and Spores" Studies
Marine dinoflagellate cysts (dinocysts) and other non-pollen palynomorphs (mainly freshwater green algae, acritarchs, and foraminifer organic linings) as well as pollen and spores were studied in 34 samples using a standard palynological technique [32,94].The sediment samples were treated with cold 10% HCl for 1 h to remove carbonates and with cold 48% HF for 5−10 days to remove siliceous particles.The residue was rinsed twice with distilled water and centrifuged after each step.The samples were then sonicated for 0.2 and 0.8 min and sieved through 15 µm mesh sieves to eliminate pelite.No oxidation was conducted to prevent the loss of the more fragile Protoperidinium cysts [95].A total of 100-300 dinocysts were identified in each slide using an Axiostar Plus microscope (Carl Zeiss, magnification 400× or 1000×) following the nomenclature provided in the literature [96][97][98][99][100].Pollen grains and spores were counted to quantify inputs from land vegetation using Kupriyanova and Aleshina [101,102] Kupriyanova [103], and Reille [104].Pollen and spores of terrestrial plants in the bottom sediments of the White Sea related to inputs from river runoff and aeolian transport [25].NPP, pollen, and spore concentrations were calculated based on the marker grain method using tablets of Lycopodium clavatum spores, cf.Stockmarr [105].
Special attention was paid to dinocyst assemblages as the percentages of the dinocyst taxa were determined for reconstructions of past sea surface conditions.The dinocyst assemblages were related to a combination of hydrographic parameters, among which sea surface temperature (SST), sea surface salinity (SSS), sea ice, and productivity appeared the most determinant [106,107].We used the modern analog technique (MAT) [41,106,107] to quantitatively reconstruct the SST ( • C) and SSS (PSU) in winter and summer, sea ice cover (months yr −1 > 50%), and primary productivity (gC/m 2 yr −1 ).The reconstructions were based on the dinocyst reference database updated in 2020 [42,107] that includes data from 1968 sites documenting the modern distribution of dinocysts and a wide range of hydrographical conditions (Table S2).The reconstructions were created from the five best modern analogs as identified from the distance (inversely proportional to the similarity) calculated after log-transformation of percentage data.In the case of the n = 1968 database, the threshold distance for good analogs was 1.2 [107].

Lithology and Organic Matter Content in Core PS-6066
In core PS-6066, which encompassed the entire Holocene according to the AMS 14 C ages (Table 1 and Figure 3), the sediments consisted of olive-gray silty pelite (2.42-1.40m), dark olive silty pelite (1.40-0.04m), and dark brown clayey silt (0.04-0.00 m) (Figure 4).The upper unit of core PS-6066 (1.40-0 m) correlated with the 2.20-0 m interval in core PS-6050, which both showed distinct boundaries from grayish sediments to olive sediments dated at about 7.0 cal kyr BP in core PS-6050 [22].According to the results of the grain size analysis, pelite (>95%) dominated the sequence (Figure 4).The sand content varied between 0.6 and 3.9% with abundance peaks at 195, 89, and 31 cm.
The proportions of organic carbon in the core sediments ranged from 1 to 1.5%, while they varied in surface sediments from 0 in the Gorlo Strait to 3.07% in the deepwater part of the White Sea [108].The δ 13 C org suggested changes in the nature of organic fluxes.From 2.42 to 1.60 m, δ 13 C org ranged from -24.41 to -24.57‰, which is typical of coastal marine sediments with high terrestrial material inputs.From 1.60 to 0.60 m, the δ 13 C org increased to values ranging between -24.31 and −22.68‰, indicating a decrease in terrestrial input relative to marine fluxes.The upper part of the core was characterized by lower δ 13 C org , down to -26‰, suggesting a high proportion of organic matter from a terrestrial origin [109].
The CaCO 3 varied between 1.1 and 4.5%, which corresponded to the contents in the surface sediments of about 1-5% in the central part of the White Sea [24].

Diatom Assemblages in Core PS-6066
In core PS-6066, the concentration of diatoms was high and varied from 0.49 × 10 6 to 3.79 × 10 6 valves/g.The overall taxonomic diversity of diatoms in the assemblages reached 81 taxa, among which 17 were freshwater and 64 were marine or brackish.The marine diatoms dominated (89-100%), with abundant planktonic neritic and panthalassic species.The percentage of panthalassic taxa that increased significantly toward the top of the core (from 29 to 87%) indicated an increase in water depth in the outer Kandalaksha Bay.Coscinodiscus radiatus, which is an indicator of Atlantic waters in the western Arctic seas [83][84][85], was a common panthalassic species of the assemblages.Marine sublittoral diatoms were mostly represented by the meroplanktonic Paralia sulcata, typical for the fouling of higher aquatic plants in the White Sea [110,111].The abundance of Paralia sulcata was maximal in the lower part of the core, reaching up to 54%, and gradually decreased by half up core, confirming an increase in water depth.It was also worth noting the occurrence of freshwater diatoms in core sediments (0.008-0.124 × 10 6 valves/g).The proportions of the freshwater species (Aulacoseira ambigua, Cyclotella meneghiniana, Stephanodiscus neoastrea, etc.) commonly transported to the shelf seas with river runoff varied from 1 to 8%.The main sources of riverine diatoms in the study area were the Varzuga and Umba rivers, along with numerous streams draining the Tersky and Karel'sky coasts.In the core, six diatom assemblage ecozones (DZ) were distinguished based on changes in the ecological affinities of taxa and the composition of dominant and subdominant species (Figure 5).
3.79 × 10 6 valves/g.The overall taxonomic diversity of diatoms in the assemblages reached 81 taxa, among which 17 were freshwater and 64 were marine or brackish.The marine diatoms dominated (89-100%), with abundant planktonic neritic and panthalassic species.The percentage of panthalassic taxa that increased significantly toward the top of the core (from 29 to 87%) indicated an increase in water depth in the outer Kandalaksha Bay.Coscinodiscus radiatus, which is an indicator of Atlantic waters in the western Arctic seas [83][84][85], was a common panthalassic species of the assemblages.Marine sublittoral diatoms were mostly represented by the meroplanktonic Paralia sulcatа, typical for the fouling of higher aquatic plants in the White Sea [110,111].The abundance of Paralia sulcatа was maximal in the lower part of the core, reaching up to 54%, and gradually decreased by half up core, confirming an increase in water depth.It was also worth noting the occurrence of freshwater diatoms in core sediments (0.008-0.124 × 10 6 valves/g).The proportions of the freshwater species (Aulacoseira ambigua, Cyclotella meneghiniana, Stephanodiscus neoastrea, etc.) commonly transported to the shelf seas with river runoff varied from 1 to 8%.The main sources of riverine diatoms in the study area were the Varzuga and Umba rivers, along with numerous streams draining the Tersky and Karel'sky coasts.In the core, six diatom assemblage ecozones (DZ) were distinguished based on changes in the ecological affinities of taxa and the composition of dominant and subdominant species (Figure 5).DZ VI (2.28-2.10m, ~11.7-10.9cal kyr BP) corresponded to the Preboreal.The concentrations of diatoms in the olive-gray pelite at the base of core were low and did not exceed 0.9 × 10 6 valves/g.DZ VI was characterized by the dominance (62-70%) of sublittoral species (Delphineis surirella, Grammatophora hamulifera, and G. marina), with maximum meroplanktonic Paralia sulcata (up to 54%).Most of the marine planktonic neritic diatoms were represented by relatively warm-water species (Coscinodiscus radiatus, C. perforatus, and Shionodiscus oestrupii; up to 32%), indicating advection of Atlantic waters in the Kandalaksha Bay.DZ V (2.10-1.85m) from the olive-gray pelite corresponded to the beginning of the Boreal (ca.10.9-9.5 cal kyr BP).The concentrations of diatoms slightly increased but did not exceed 1.3 × 10 6 valves/g.The marine sublittoral diatoms remained predominant (up to 62%) with a high proportion of the meroplanktonic species Paralia sulcata (up to 48%).The abundances of the Atlantic water species such as Shionodiscus oestrupii, Coscinodiscus radiatus, C. perforatus, and Actinoptychus senarius increased up to 42%.

Palynomorph Assemblages in Core PS-6066
In the modern sediments of the White Sea and Barents Sea, the influence of Atlantic water masses is marked by the occurrence of the phototrophic dinocyst species, Operculodinium centrocarpum [32].Assemblages with common O. centrocarpum are typical of the western shelf edge of the Barents Sea and West Spitsbergen [112,113], as well as of the White Sea [32].Islandinium minutum subsp.minutum often dominate in sediments of areas marked by Arctic surface waters mass distribution area, together with other heterotrophic taxa, notably Brigantedinium spp.[112,113] in the Barents Sea, and Echinidinium karaense and Selenopemphix quanta in the White Sea [32].
Micropaleontological analysis in core PS-6066 showed the presence of common palynomorphs, including dinoflagellate cysts (dinocyst), freshwater green algae, acritarchs, foraminifer organic linings, pollen grains, and spores.The most abundant palynomorphs were dinocysts (up to 271 × 10 3 cysts/g), pollen and spores of terrestrial plants (up to ~2 × 10 6 grains/g), and green algae (up to 7.4 × 10 3 sp./g).Dinocysts were represented in the core by 14 taxa.In general, cosmopolitan species such as O. centrocarpum were associated with the penetration of Atlantic waters [114,115].The palynomorphs concentrations, percentages of the main dinocyst taxa, and grain size led to the definition of three main assemblage zones in the PS-6066 core (Figure 6).

Palynomorph Assemblages in Core PS-6066
In the modern sediments of the White Sea and Barents Sea, the influence of Atlantic water masses is marked by the occurrence of the phototrophic dinocyst species, Operculodinium centrocarpum [32].Assemblages with common O. centrocarpum are typical of the western shelf edge of the Barents Sea and West Spitsbergen [112,113], as well as of the White Sea [32].Islandinium minutum subsp.minutum often dominate in sediments of areas marked by Arctic surface waters mass distribution area, together with other heterotrophic taxa, notably Brigantedinium spp.[112,113] in the Barents Sea, and Echinidinium karaense and Selenopemphix quanta in the White Sea [32].
Micropaleontological analysis in core PS-6066 showed the presence of common palynomorphs, including dinoflagellate cysts (dinocyst), freshwater green algae, acritarchs, foraminifer organic linings, pollen grains, and spores.The most abundant palynomorphs were dinocysts (up to 271 × 10 3 cysts/g), pollen and spores of terrestrial plants (up to ~2 × 10 6 grains/g), and green algae (up to 7.4 × 10 3 sp./g).Dinocysts were represented in the core by 14 taxa.In general, cosmopolitan species such as O. centrocarpum were associated with the penetration of Atlantic waters [114,115].The palynomorphs concentrations, percentages of the main dinocyst taxa, and grain size led to the definition of three main assemblage zones in the PS-6066 core (Figure 6).Zone III at the base of the core (2.42-1.65 m; 11.7-8.4cal kyr BP) was subdivided into three subzones.Subzone IIIc (2.42-2.10m; 11.7-10.9cal kyr BP) was characterized by the dominance of Operculodinium centrocarpum that includes various morphological types (up to 71%) and Spiniferites spp.(up to 23%).The heterotrophic species Selenopemphix quanta, Islandinium minutum subsp.minutum, and Echinidinium karaense do not exceed 6.9%.Dinocyst concentrations are 30-42 × 10 3 cysts/g, and the content of freshwater green algae is low, not exceeding 0.7 sp./g.The concentrations of pollen grains and spores were 167-280 × 10 3 grains/g and were much higher than those of marine NPP (Figure 6).
In subzone IIIb (2.10-1.92m; 10.9-9.9 cal kyr BP), the species diversity of dinocysts increased, and the Impagidinium pallidum, which is typical for fully marine conditions with high salinity, appeared singly [116].At a depth of about 2 m, the sediment became more fine-grained (pelite ~98%), but then the sandy fraction reached up to 1.7%.
A slight increase in coarse-grained sediment in subzone IIIa (1.92-1.65 m; 9.9-8.4cal kyr BP) was accompanied by a sharp short-term increase in heterotrophic species (I.minutum, E. karaense, S. quanta, and Brigantedinium spp.) up to 8-10%, although autotrophic cosmopolitan species generally dominated [32].Heterotrophic dinocysts were associated with an environment with dense seasonal sea ice cover [96,97,106,107].In addition, the Nematosphaeropsis labyrinthus (about 1%), found exclusively in completely marine conditions, was encountered for the first time, and the frequency of the findings of Impagidinium pallidum (1-2%) also increased.Further, the concentrations of dinocysts reached up to 270 × 10 3 cysts/g, and the freshwater green algae (up to 7.4 sp./g) and the pollen and spores (up to 2 × 10 6 grains/g) sharply increased (Figure 6).Subzone IIIa was characterized also by the diversity of spores and pollen and the predominance of up to 99% of tree species (mainly Pinus and Picea genus) with single grains of Betula nana and Alnus.However, it should be noted that there was a sharp maximum in the proportion of Sphagnum and Polypodiaceae spores (almost up to 4% in total) at the zone boundary at 1.92 m.
Zone II (IIc, 1.65-1.38m; 8.4-6.9 cal kyr BP) was characterized by a decrease in the proportion of the O. centrocarpum species to 42% and was accompanied by an increase in other autotrophic species: Pentapharsodinium dalei cysts up to 9% and Bitectatodinium tepikiense up to 16%.The presence of the latter two species in modern sediments of the Arctic seas was often associated with low-duration (<4 months a year) sea ice cover, as well as with the passage of the polar front and/or increased productivity [97,116].
Subzone IIb (1.38-1.22m; 6.9-6 cal kyr BP) was characterized by a slight decrease in the proportion of heterotrophic species down to 2%, while autotrophic species reached their maximum.In addition to the general predominance of cysts of the genus Operculodinium (up to 61%), the cysts of S. ramosus and S. elongatus reached 44% by the end of the period.In addition, almost all samples contained single cysts of Impagidinium pallidum up to 1%.
Subzone IIa (1.22-0.92m; 6-4.4 cal kyr BP) was characterized by increases in the proportion of heterotrophic relatively cold-water species up to 7.5% and was noted especially at a depth of 1.16-1.02m.In addition, the proportion of the Operculodinium cyst, which currently occupies about 50% in the associations, was reduced.
The palynomorph Zone I (Ib, 0.92-0.60m; 4.4-2.7 cal kyr BP) was marked at its base by an increase in the concentrations of aquatic and terrestrial palynomorphs.The dinocysts content reached a second maximum with values of 260 × 10 3 cysts/g.The concentrations of green algae and pollen and spores reached 6 and 754 × 10 3 grains/g, respectively.Toward the top of the zone (Ia, 0.60-0.00m; 2.7 cal ka BP to the present), the palynomorph concentrations decreased.
In Zone I, the dinocyst assemblages were characterized by high percentages of Operculodinium centrocarpum and a decreasing trend of Spiniferites spp.Rare cysts of Nematosphaeropsis labyrinthus and Impagidinium pallidum, which are typical for subpolar environments of the northeast North Atlantic [106,107,117], were recovered.

Changes in P/B Ratio Values
The P/B ratio (Figure 7) expresses the proportion of diatom valves belonging to planktonic species to the sum vs. benthic and planktonic species [118][119][120] following Wang et al. [119].
In core PS-6066 from the outer Kandalaksha Bay, the P/B ratio ranges between 0.29 and 0.86, increasing from the base to the top, which suggests deepening of the water depth.The increase in the P/B ratio from 0.29 to 0.63 occurred before 5.0 cal kyr BP.Two low P/B value peaks (0.62 and 0.60) at about 8.4 and 6.8 cal kyr BP were noted.In core PS-6066 from the outer Kandalaksha Bay, the P/B ratio ranges between 0.29 and 0.86, increasing from the base to the top, which suggests deepening of the water depth.The increase in the P/B ratio from 0.29 to 0.63 occurred before 5.0 cal kyr BP.Two low P/B value peaks (0.62 and 0.60) at about 8.4 and 6.8 cal kyr BP were noted.
We calculated the P/B ratios from the data of the previously studied cores PS-6050 and PS-6042 (Figures 1b and 7) collected in the outer and central Dvina Bay, at 101 and 61 meters of water depth, respectively [21,22,28].The P/B ratio of core PS-6050 showed a trend from 0.3 to 0.8 over the past 8.6 cal kyr BP, which also suggested a temporary increase in water depth in the outer Dvina Bay.Core PS-6042, from the central part of the Dvina Bay, differed significantly from the deep-water areas of the White Sea (Figure 7).The P/B ratio increased from 8.8 to 1.5 cal kyr BP, with the most pronounced change occurring from 8.6 to 5.6 cal kyr BP.The last 1.6 kyr was, however, marked by a decrease in the P/B ratio by almost half, which could have been due to local changes in the water depths most likely caused by local block movements in the Dvinsky graben.

Reconstruction of Sea Surface Water Conditions
The sea surface conditions were estimated based on the modern analog technique (MAT) using the n = 1968 dinocyst database [42,109] that was applied to assemblages of cores PS-6066 (this study), PS-6042, and PS-6050 [21,22,28].The distance between fossil and modern analogs mostly varied between 0.3 and 0.6 and never exceeded 0.9 for the first analog, which indicated good modern analogs exist for most samples.The MAT results showed the variation of hydrographic parameters in time and space (Figure 7).
From >11.4 to ~9.5 cal kyr BP, we reconstructed the harsh conditions in the central Dvina Bay with the sea ice cover of 8-10 months yr −1 , summer SST of 2-8 °C, and SSS fluctuating between 27 and 32 PSU.At the deepest site of the Kandalaksha Bay, the sea surface conditions were less severe with 4-6 months yr −1 of sea ice cover, SSS around 30  [23]).The Kandalaksha Bay is additionally highlighted by gray.
We calculated the P/B ratios from the data of the previously studied cores PS-6050 and PS-6042 (Figures 1b and 7) collected in the outer and central Dvina Bay, at 101 and 61 meters of water depth, respectively [21,22,28].The P/B ratio of core PS-6050 showed a trend from 0.3 to 0.8 over the past 8.6 cal kyr BP, which also suggested a temporary increase in water depth in the outer Dvina Bay.Core PS-6042, from the central part of the Dvina Bay, differed significantly from the deep-water areas of the White Sea (Figure 7).The P/B ratio increased from 8.8 to 1.5 cal kyr BP, with the most pronounced change occurring from 8.6 to 5.6 cal kyr BP.The last 1.6 kyr was, however, marked by a decrease in the P/B ratio by almost half, which could have been due to local changes in the water depths most likely caused by local block movements in the Dvinsky graben.

Reconstruction of Sea Surface Water Conditions
The sea surface conditions were estimated based on the modern analog technique (MAT) using the n = 1968 dinocyst database [42,109] that was applied to assemblages of cores PS-6066 (this study), PS-6042, and PS-6050 [21,22,28].The distance between fossil and modern analogs mostly varied between 0.3 and 0.6 and never exceeded 0.9 for the first analog, which indicated good modern analogs exist for most samples.The MAT results showed the variation of hydrographic parameters in time and space (Figure 7).
From >11.4 to ~9.5 cal kyr BP, we reconstructed the harsh conditions in the central Dvina Bay with the sea ice cover of 8-10 months yr −1 , summer SST of 2-8 • C, and SSS fluctuating between 27 and 32 PSU.At the deepest site of the Kandalaksha Bay, the sea surface conditions were less severe with 4-6 months yr −1 of sea ice cover, SSS around 30 PSU, and SST 8-10 • C (Figure 7).From 9.5 to ~7 cal kyr BP a warming trend was reconstructed.The outer part of the Kandalaksha Bay became almost free of sea ice, and the central Dvina Bay recorded the minimum sea ice (~2 months yr −1 ), while SST reached the maximum of 12-14 • C, and the annual primary productivity increased up 500 gC m −2 yr −1 .In the second half of the Middle Holocene at 7-6 cal kyr BP, there was an increase in the duration of seasonal ice cover and a decrease in SSS (to 28-30 PSU) and SST (to about 8 • C).After 6 cal kyr BP, there were few fluctuations in the sea surface parameters.The last notable change in surface water occurred at about 3.5 cal kyr BP when SSS decreased to 28 PSU in the central Dvina Bay, and SST increased up to 2-3 • C. Around 0.5 cal kyr BP, high temperatures were reconstructed in Kandalaksha and the outer part of Dvina Bay (up to 14-16 • C), while they were 2-4 degrees lower in the central Dvina Bay (Figure 7).
It should be noted that the milder conditions in the Kandalaksha Bay than in the Dvina Bay may be associated with the transformed Atlantic water distributions that flow into the Kandalaksha Bay along the Tersky Coast, prior to reaching the Dvina Bay, thus resulting in local particularity in the "hydrological pole of cold" and the "hydrological thermal pole".

Holocene Sea Level Changes in the White Sea Area
Reconstruction of the sea level in the White Sea area during the Holocene is a complicated task since the factors that operate in this region include eustatic sea level rise, glacioisostatic uplift, and regional neotectonic movements.During the last glacial maximum around 21,000 years ago, the global sea level was >120 m lower than at present [12,[16][17][18][121][122][123][124].The average rates of eustatic sea level rise (SLR) were uneven [15,18,125].From 16.5 to ~8.2 cal kyr BP, it was about 12 mm yr −1 with faster rates of ~40 mm yr −1 at ~14.5-14.0cal kyr BP and 11.3 cal kyr BP and slower rates during the Younger Dryas, 12.5-11.5cal kyr BP [18,19].The modern eustatic sea level was reached by the end of the Middle Holocene [15,18].
Elevated shorelines evidence isostatic uplift along the White Sea coasts, notably along the Kandalaksha Bay at altitudes up to 100 meters [14,[126][127][128][129] and in the Dvina Bay at 5−20 m above the present sea level [129].The discrepancy in the elevation of the shorelines results from the glacioisostatic processes [14,[127][128][129][130][131][132].Our records from the outer Kandalaksha Bay indicated local deepening of the water depth since the beginning of the Preboreal.The evidence came from the consistent increase in planktonic diatoms from 29 to 87% in core PS-6066 (Figure 5) relative to the decrease in sublittoral diatoms as indicated by the P/B ratio increasing from 0.29 to 0.86 (Figure 7).The record of core PS-6066 was consistent with data from the central White Sea [20] (Figure 1b, core 253) and the inner part of the Kandalaksha Bay [14,128,129].
The deepening of the White Sea was primarily caused by eustatic sea level rise, which was about 50 meters between 11.7 and 5 cal kyr BP [15].However, the diatom data suggested that the basin deepening continued until recently, which points to neotectonics as the leading factor.The Kandalaksha graben of the White Sea is one of the most active tectonic zones of the eastern part of the Fennoscandian crystal shield.This is confirmed by a system of Holocene active faults, earthquake, and paleoseismic dislocations [133,134].
The Dvina Bay is located within the southeastern tip of the Kandalaksha-Dvinsky Rift, which extends in a northwesterly direction from the Kandalaksha Bay through the central part of the White Sea to the Dvina Bay.It is assumed that this rift experienced continuous subsidence during the Holocene.
Previous diatom investigations of the Dvina Bay sediments [21,22] indicated a gradual increase in marine planktonic taxa.It increased from 29 to 80%, and the P/B ratio increased (from 0.3 to 0.8, Figure 7) toward the top of core PS-6050 from the outer Dvina Bay, indicating deepening over the past 8.6 cal kyr BP [21].Similarly, in the sediments of the central part of the Dvina Bay (core PS-6042), the planktonic diatoms consistently increased from 9% to 79%, and the proportion of marine sublittoral species decreased from 91 to 21% [22].From these data, we may conclude that the water depth of the Dvina Bay increased in its outer part and at least until past 1.6 cal kyr in its central part.
The decrease in the P/B ratio in the central Dvina Bay sediments since about 1500 years ago might be caused by decreased water depths locally, which could be due to local tectonic block movements [135,136].A decreased bathymetry in the central Dvina Bay during the middle of the Subatlantic was also evidenced by the formation of the Northern Dvina River delta at about 2 cal kyr BP and by coastal ridges on the nearby islands [137].
Our records indicated that while the coasts of the Kandalaksha Bay continued to uplift in the Holocene, the deep-water areas of its outer part and the adjacent Basin, as well as the outer part of the Dvina Bay, experienced sea level rise.

Holocene Environments in the White Sea
The main environmental changes of the late glacial and postglacial history of the White Sea were related to water exchange with the Barents Sea, which started between 11.7 and 11.2 cal kyr BP [14] when sea level reached about the 40-42 m deep Gorlo Strait sill [11][12][13].According to previous studies, there was permanently a strong current from the Barents Sea, carrying transformed Atlantic waters into the White Sea basin along the Tersky coast to the Kandalaksha Bay (Figure 1b), as at present [11,20].
In the outer Kandalaksha Bay, at the location of core PS 6066, fully marine conditions were already established by the beginning of the Preboreal at ~11.7 cal kyr BP as indicated from diatom (Figure 5) and palynomorph assemblages (Figure 6, subzone IIIc) and δ 13 C org values (Figure 4).The dominant sublittoral diatoms and low P/B ratios in the olive-gray pelite with hydrotroilite, however, suggested lower sea level in the outer Kandalaksha Bay than at present (Figure 7).
During the Preboreal, the marine waters of the Kandalaksha Bay penetrated into the land, as evidenced by the abundance of neritic euryhaline diatom species in sediment cores from the upstream part of bay [20,24].The high abundances of planktonic diatoms in these sediments may indicate the likelihood of high water productivity all over the bay.
The compositions of diatom and aquatic palynomorph assemblages were similar to those recovered near surface sediments (Figures 5 and 6, core PS-6066), which suggested an early onset of modern-like hydrological and hydrobiological conditions.Of particular interest was the confirmation that the penetration of transformed Atlantic waters into the Kandalaksha Bay started at ~11.7 cal kyr BP.The high concentrations and proportions of relatively warm-water planktonic diatoms (Coscinodiscus radiatus and Shionodiscus oestrupii) and dinocysts (O.centrocarpum, cyst of P. dalei), which are indicators of Atlantic-sourced waters in the western Arctic [84][85][86]114,115], proved the wide water exchange between the White and the Barents Seas and the intensive advection of Atlantic waters during the Preboreal.The data led to the inference of an increase in SST, from 2 to 12 • C, and phytoplanktonic productivity up to 730 gC m −2 g −1 .
Our reconstructions are in a good agreement with the diatom and dinocyst data from the northwestern continental slope of the Barents Sea [136].Based on the records from of the well-dated core SV-04, it was shown that the modern-like conditions characterized by the intensive advection of Atlantic waters were established early during the Holocene (11.2 cal kyr BP) [138].
At the beginning of the Preboreal (~11.7-11.4cal kyr BP), transformed Atlantic waters also entered the central part of the Dvina Bay, although with lower fluxes as evidenced by diatom and dinocyst of core PS-6042 [22].Based on low diatom and NPP concentrations, maximum freshwater green algae input, and low δ 13 C org (-27.04-28.7),we deduced intensive input of terrigenous organic matter at ~11.7 cal kyr BP into the central Dvina Bay [22,28].This was consistent with dinocyst data that led to the reconstruction of relatively low salinity (29-30 psu) and productivity (82-457 gC m −2 yr −1 ) in the surface waters.This was likely caused by the proximity of study site PS6042 to the Northern Dvina River mouth at a time of lower sea level, by about −40 to −35 m compared with the present [16,139].
During the Boreal (~10.9-8.4 cal kyr BP), a gradual deepening of the outer water depth in the Kandalaksha Bay (PS-6066) was evidenced by a decrease in the relative abundances of sublittoral diatoms and a rise of the P/B ratios (Figure 7).The species compositions of planktonic diatom and dinocyst assemblages indicated changes in the hydrological conditions.
The first half of the Boreal (Figure 5, ~10.9-9.5 cal kyr BP) was characterized by a slight reduction in the influence of transformed Atlantic waters in the outer Kandalaksha Bay, as inferred from the general decrease in the concentrations of diatoms and dinocysts.An increase in the duration of seasonal sea ice cover (2-4 months yr −1 ) and a decrease in sea surface temperature (down to 2-8 • C) and annual productivity (100-200 gC m −2 yr −1 ) were reconstructed.Nevertheless, the eustatic sea level rise was reflected in the increasing sea surface salinity of the White Sea record.In the second half of the Boreal (~9.2-8.4 cal ka BP), there was a sharp increase in concentrations of palynomorphs in sediments.A shorter duration of seasonal sea ice cover (0-2 months yr −1 ) may have led to an increase in productivity (189-489 gC m −2 yr −1 ) in surface waters.In the outer part of the Kandalaksha Bay, the second half of the Boreal period 9.5-8.4cal ka BP was marked by an increase in SST and SSS (11-12 • C and 32-34 PSU).
In the central Dvina Bay, core PS 6042 contains sediments dating back to the Boreal period.In the first half of the Boreal (10.9-9.5 cal kyr BP), the sea level rise led to a shift of the Northern Dvina River mouth toward its present position.This was marked by a diminution of the green algae concentrations and changes of δ 13 C org to more positive values.Surface water temperatures (0-2 • C) during this period were slightly lower than in the Kandalaksha Bay despite the occurrence of Atlantic water indicator species in the diatom and aquatic palynomorph assemblages.During the second half of the Boreal (9.4-8.2 cal kyr BP), the SST and SSS increased (5-10 • C and 28-32 PSU), and the seasonal sea ice cover decreased (from 5-10 to 1-5 months yr −1 ).
Sea ice formation in the White Sea starts first at the mouths of bays, where large river fluxes result in low-salinity surface waters, thus fostering early freezing, 1-3 months earlier than more offshore.The extent of seasonal sea ice cover is also governed by wind.By April-May, the "spring bloom" of sea ice diatoms and the active growth of phytoplankton start in the Arctic seas when 50-80% of the sea ice is dispersed to the Barents Sea through the Gorlo Strait by strong winds [29].The sea ice distribution in the White Sea is reflected in diatom and dinocyst assemblages recovered from bottom sediments as shown from previous studies [31,32,83].Therefore, we hypothesized that the extensive sea ice and minimal surface water temperatures characterizing the first half of Boreal in the central part of the Dvina Bay were caused by freshening in surface waters due to the relatively close location of the Northern Dvina River mouth [29].
The warming of the end of the Boreal after 9.5 cal kyr BP in the Kandalaksha and Dvina Bays was interrupted by a cooling outer Dvina Bay and outer Kandalaksha Bay (Figure 7).Relatively high percentages of the heterotrophic species I. minutum subsp.minutum and E. karaense suggest more extensive sea ice cover (up to 10 months/year in the outer Dvina Bay) and a decrease in surface water temperatures in the Kandalaksha and Dvina Bays of the White Sea (Figure 7).In the first part of Atlantic period (~8.4-7 cal kyr BP) the outer Kandalaksha Bay recorded an increase in SST, up to 14 • C which may coincide with the thermal optimum the Holocene recorded in the Barents Sea to occur from 7.8 ka B.P. to 6.8 ka [140].This interval was also marked by the deposition of olive-gray pelite in the outer Kandalaksha Bay, the occurrence of relatively warm-water diatom and NNP species (Figure 5), high carbonate content (CaCO 3 ~3-4%), and low terrigenous input (δ 13 C org ranging from -22.68 to −24.31‰).The SST rise reached its maximum at about 8-7.5 cal kyr BP (Figure 7).In the central Dvina Bay, SST slightly decreased (down to 6-8 • C) during the first half of the Atlantic, whereas it increased from 3 to 12 • C in the outer Dvina Bay.
We hypothesized that higher SST in the Kandalaksha Bay and the outer Dvina Bay in comparison with the central Dvina Bay during the first half of the Atlantic were related to transformed Atlantic waters entering into the White Sea, first into the Kandalaksha Bay along the Tersky Coast and then into the Dvina Bay.The regional differences may reflect the local structures of the "hydrological pole of cold" and the "hydrological thermal pole".The warming in the Kandalaksha Bay and the outer Dvina Bay is in good agreement with the SST changes in the first half of the Atlantic in the Barents Sea (Figure 7) [23,140].
In the outer Kandalaksha Bay, there was a decrease in SST and SSS (5 • C and 29-31 PSU) and productivity  gC m −2 yr −1 ), while sea ice cover extended to 6 months yr −1 during the second part of the Atlantic (about 7-5.8 cal kyr BP).Such a mid-Holocene cooling was also recorded in the Barents Sea after 6.7 kyr BP [140].In the White Sea the cooling was also evidenced from the increasing number of cold-water diatom species Porosira gracialis, Bacterosira bathyomphala, and Rhizosolenia hebetata f. hebetata and the occurrence of the sea ice species Nitzschia frigida and Chaetoceros furcellatus (Figure 5).
In the outer Dvina Bay, SST, SSS, and productivity decreased down to 6 • C, 25 PSU, and 149 gC m −2 yr −1 , respectively, whereas sea ice cover increased up to 6 months yr −1 .In the central Dvina Bay, the conditions remained similar to those of the first half of the Atlantic.
The maximum diatom concentrations in the outer Kandalaksha Bay (core PS-6066, the upper part of DZ III, Figure 5) and in the outer and central Dvina Bay [21,22] occurred at about 6.7-6 cal kyr BP.Thus, the White Sea sediments recorded a maximum concentration of diatoms, which may correspond to a change in the sedimentary regime, as well as an increase in the productivity of diatoms.We interpreted the increased diatoms and NPP concentrations as the result of the input of biophilic trace elements from seasonal sea ice.Apparently, the mechanism of sea ice melt and opening differed from the modern one and was primarily associated with wind-driven upwelling.This was consistent with the occurrence of sea ice and ice-neritic diatom species in the outer Kandalaksha Bay and the Dvina Bay [21,22] sediments.The high Si content could alternatively derive from river discharge and the underground runoff, which increased due to relict permafrost thawing [141].
In the outer Kandalaksha Bay, the beginning of the Subboreal (after ~5.8 cal kyr BP) was marked by increased SST (12 • C) and SSS (34 PSU).In the outer and central Dvina Bay, SST and SSS also increased up to 8-14 and 8-13 • C and up to 34 and 30 PSU, respectively.Nonetheless, the decrease in diatom concentration, the increase in the number of relatively cold-water species, and the occurrence of the diatom sea ice species Attheya septentrionalis in the Kandalaksha (Figure 5, ~6-4.4 cal kyr BP) and Dvina Bays [21,22] indicated a cooling in the first half of the Subboreal.We also assumed that the presence of sea ice species in sediments may have been associated with the seasonal sea ice cover boundary's closer position to the central parts of the bays during the growing season.This corresponded well with two high sea ice concentration events (~5 and 3.5 cal kyr BP) from the Kandalaksha Bay.In the outer Kandalaksha Bay, the dinocyst data indicate the SSS was about 31-34 PSU and SST-8-10 • C during the second half of the Subboreal.In the outer Dvina Bay, SST increased from 8 to 14 • C, whereas in the central Dvina Bay, it decreased from 13 to 8 • C, and SSS reached 31-34 PSU.The slight increase in the number of warm-water species and diatom concentrations in the outer Kandalaksha might have been due to the increase in the intensity of transformed Atlantic waters.After about 3.5 cal kyr BP, there was a decrease in salinity in the study areas, which could be linked with a decrease in Atlantic water fluxes, as also evidenced in the south of the Barents Sea (Figure 7) [23].
At the beginning of the Subatlantic, about 2000 years ago, there was a short-term increase in SSS up to 34 PSU, probably related to an increase in the intensity of transformed Atlantic waters.In the second part of the Subatlantic, after ~2.0 cal kyr BP, a consistent increase in the number of relatively warm-water diatom and dinocyst species was noted in the outer Kandalaksha Bay, outer Dvina Bay, and central Dvina Bay [21,22,28].SST values increased during the Subatlantic up to 13 • C in the outer Kandalaksha Bay, to 16 • C in the outer Dvina Bay, and to 11 • C in the central Dvina Bay.The trend toward lower surface temperatures in the outer Dvina Bay that was particularly pronounced after 2.0 cal kyr BP was probably due to upwelling of cold bottom waters.
The increase in freshwater species (up to 8%) in the outer Kandalaksha Bay core (Figure 5, DZ I) was associated with an increase in the river discharge after ~2.0 cal kyr BP, which was also noted in the central Dvina Bay (up to 40%) and the outer Dvina Bay (up to 8%) cores [21,22].It was also related to the lowering of sea level along the coastline of the Kandalaksha Bay [14,128,129,142,143] that fostered the influence of the distribution of river fresh waters.The freshwater diatom distributions in the Dvina Bay, up to 40% in the central part and 8% in the outer part, were similar to the modern ones in surface sediments [31].The increase in river surface runoff was also evidenced from the isotopically light δ 13 C org of -26‰ in the upper part of core PS-6066 [109], in agreement with studies of the White Sea drainage basin [111,139].According to the regional physiography and radiocarbon dating, an increase in river runoff accompanied with the formation of new channels was established in the Subatlantic time period [141,144,145].

Conclusions
Diatom and palynomorph assemblages from the 14 C-dated sediment core PS-6066 provided new information on the postglacial history in the Kandalaksha Bay, the White Sea, over the Holocene.Data indicated a deepening of the Kandalaksha Bay throughout the Holocene.
The penetration of transformed Atlantic waters into the Kandalaksha Bay started at least 11,700 years ago and between ~11.7 and 11.4 cal kyr BP into the central part of the Dvina Bay (core PS-6042).The data also indicated that the highest SST prevailed in the White Sea from 8 to 7.5 cal ka BP.
In the Kandalaksha (PS-6066) and Dvina (PS-6042, PS-6050) Bays, there was a cooling event in the surface waters marked by relatively low temperatures and salinity, a decrease in productivity, and an increase in sea ice for more than 6 months of the year at the end of the middle Holocene at about ~7-5.8 cal kyr BP.
In the central Dvina Bay, a decrease in surface temperatures more pronounced than in the outer Kandalaksha and outer Dvina Bays was recorded after 2.0 cal kyr BP.The upwelling of cold and nutrient-rich bottom waters could be invoked.Furthermore, the increase in the freshwater diatom species in the outer Kandalaksha Bay and the central and outer Dvina Bay could be associated with increasing river discharge after ~2.0-1.0 cal kyr BP.
The overall records of the White Sea evidence the impact of sea level and hydrological changes in sea surface conditions, including sea ice cover and primary productivity.Our local study of Holocene natural changes in the White Sea contributes to a broader understanding of climate change in the Barents Sea region as a whole and will be considered in further research.

Figure 1 .
Figure 1.Map showing the locations of the studied cores and modern surface currents in the North Atlantic (a) and the White Sea (b).Arrows indicate the direction of surface currents: ESC-East Spitsbergen Current; WSC-West Spitsbergen Current; NAC-North Atlantic Current; NCC-Norwegian Coastal Current.Sediment core sites from this and previous studies are indicated by red and blue circles, respectively: Core 253 [20], PS-6050, PS-6042 by[21,22], and PL 96-112 by[23].

Figure 1 .
Figure 1.Map showing the locations of the studied cores and modern surface currents in the North Atlantic (a) and the White Sea (b).Arrows indicate the direction of surface currents: ESC-East Spitsbergen Current; WSC-West Spitsbergen Current; NAC-North Atlantic Current; NCC-Norwegian Coastal Current.Sediment core sites from this and previous studies are indicated by red and blue circles, respectively: Core 253 [20], PS-6050, PS-6042 by[21,22], and PL 96-112 by[23].

Figure 2 .
Figure 2. Map of mean summer (a) sea surface salinity and (b) temperature in the White and Barents Seas between 1955 and 2012, after the World Ocean Atlas [50,51].

Figure 2 .
Figure 2. Map of mean summer (a) sea surface salinity and (b) temperature in the White and Barents Seas between 1955 and 2012, after the World Ocean Atlas [50,51].

Figure 4 .
Figure 4. Lithological units; grain size; total organic carbon; calcium carbonate; carbon isotope composition of organic carbon, marked in blue; concentrations of diatoms, marked in gray; dinocysts; green algae; acritarchs; organic foraminifera linings; and pollen and spores (marked in green) expressed in the number of specimens per gram of sediment in core PS-6066, the Kandalaksha Bay.The age model dates in cal kyr BP are marked in red, and the boundary dates in cal kyr BP are black.

Figure 6 .
Figure 6.Dinocyst concentrations and relative abundance (percentages of the main taxa in core PS-6066.Dark gray and light gray mark heterotrophic and autotrophic species, respectively.The Holocene stages and dinocyst ecozones are marked on the right.Age model dates in cal kyr BP are marked in red, and boundary dates in cal kyr BP are in black.Zone III at the base of the core (2.42-1.65 m; 11.7-8.4cal kyr BP) was subdivided into three subzones.Subzone IIIc (2.42-2.10m; 11.7-10.9cal kyr BP) was characterized by the dominance of Operculodinium centrocarpum that includes various morphological types (up to 71%) and Spiniferites spp.(up to 23%).The heterotrophic species Selenopemphix quanta,

Figure 6 .
Figure 6.Dinocyst concentrations and relative abundance (percentages of the main taxa in core PS-6066.Dark gray and light gray mark heterotrophic and autotrophic species, respectively.The Holocene stages and dinocyst ecozones are marked on the right.Age model dates in cal kyr BP are marked in red, and boundary dates in cal kyr BP are in black.

Figure 7 .
Figure 7. Diatom P/B ratio and quantitative estimates of sea surface parameters in the White Sea (this study) in and the Barents Sea (PL 96-112 by[23]).The Kandalaksha Bay is additionally highlighted by gray.

Figure 7 .
Figure 7. Diatom P/B ratio and quantitative estimates of sea surface parameters in the White Sea (this study) in and the Barents Sea (PL 96-112 by[23]).The Kandalaksha Bay is additionally highlighted by gray.

Funding:
The sedimentological and micropaleontological study was supported by the Russian Science Foundation (project No 21-17-00235).The diatom researches of the Dvina Bay sediments were carried out within the framework of the State Assignments of Ministry of Science and High Education, Russia (Lomonosov Moscow State University, theme no.121051100135-0).The studies by A. de V are supported by the Natural Sciences and Engineering Council of Canada.Data Availability Statement: Not applicable.