Middle Holocene Climate Oscillations Recorded in the Western Dvina Lakeland

: Although extensive archeological research works have been conducted in the Serteya region in recent years, the Holocene climate history in the Western Dvina Lakeland in Western Russia is still poorly understood. The Neolithic human occupation of the Serteyka lake–river system responded to climate oscillations, resulting in the development of a pile-dwelling settlement between 5.9 and 4.2 ka cal BP. In this paper, we present the quantitative paleoclimatic reconstructions of the Northgrippian stage (8.2–4.2 ka cal BP) from the Great Serteya Palaeolake Basin. The reconstructions were created based on a multiproxy (Chironomidae, pollen and Cladocera) approach. The mean July air temperature remained at 17–20 ◦ C, which is similar to the present temperature in the Smolensk Upland. The summer temperature revealed only weak oscillations during 5.9 and 4.2 ka cal BP. A more remarkable feature during those events was an increase in continentality, manifested by a lower winter temperature and lower annual precipitation. During the third, intermediate oscillation in 5.0–4.7 ka cal BP, a rise in summer temperature and stronger shifts in continental air masses were recorded. It is still unclear if the above-described climate ﬂuctuations are linked to the North


Introduction
Neolithization processes in Central and Northern Europe were mainly caused by climate fluctuations, especially during the period between 8.2 and 6 ka cal BP [1]. The cold event in 8.2 ka cal BP caused the migration of Neolithic communities northward from the steppe area of the Black Sea Lowland region toward the coniferous and mixed forest belt of the East European Plain [2]. The rise of temperature in Northern Europe in ca. 6 ka cal BP favored agriculture dispersion and an increase of population in the circum-Baltic Region [3]. In the Central-Western European territories, the expansion of Linear Pottery Culture communities appears to coincide with wet and relatively warm winters and cool summers [1], which was also supported by quantitative paleoclimatic studies. Several paleotemperature records from the E Scandinavia and NW Russia were cited by Kaufman et al. [4], e.g., two quantitative reconstructions based on pollen data [5,6] and one based on Chironomidae [7,8].
Paleoclimatic research is less developed in Eastern Europe compared to the western part of the continent [2,[9][10][11]. Thus, the paleoclimatic background for archeological discoveries on Mesolithic and Neolithic cultures in Eastern Europe is still insufficient. Good qualitative data, unfortunately, the temperature reconstructions are not available for the following publication: Tallinn region (Estonia) [12], Myshetskoe-Dolgoe Lake in Moscow [13], Lake Kenozero [14], Ukraine and Belarus [15,16]. However, plant communities in the area reveal a lagged response to short-lived temperature oscillations that might have been crucial for human settlement, which made it possible to adopt different nutritional strategies at different times.
Aquatic invertebrates are more reactive and sensitive to short-term changes in temperature than plants [17]. Chironomids (Chironomidae), a midge group within the Culicomorpha (Diptera) infraorder, have short life cycles and respond sensitively to climate changes with a decadal resolution [18]. Their development is strongly influenced by environmental conditions [19,20]. The characteristics of Chironomids, namely species diversity, short life cycles, and sensitivity to environmental conditions, make them one of the most reliable proxies for summer air reconstructions [18,21].
Although the summer temperature is crucial for vegetation in temperate to boreal biomes, and influences the economy of the Neolithic communities, precipitation also determines primitive agriculture as well as a hunter-gatherer subsistence strategy. The effective long-term precipitation may be estimated from changes in lake level [22]. Alongside chironomids, small crustaceans of the order Cladocera are very sensitive to lake level fluctuations, as many of the Cladocera taxa are planktonic [23,24]. Therefore, cladocerans can be successfully used for water level reconstructions [25,26]. These microcrustaceans also indicate changes in the trophic status of lakes, including human-environment relationships [27], lake levels [28] including overbank episodes in floodplain areas [29,30], and climate change [31,32]. The Holocene lake level changes are very well documented for Eastern Europe. In the case of the East Peribaltic Region, lake levels and climate wetness were recorded from Northern Poland [33][34][35] and Baltic States [36][37][38][39][40]. Currently, the entire Holocene sequence from the East European Plain can be found in as many as 50 lake records [41,42]. From Eastern Europe quantitative reconstruction have been made from pollen-inferred paleotemperature and precipitation at Mshinskoye raised bog and Lembolovskoye Lake [43], Bobrovichskoe and Oltushskoe (Polesye) [44], Upper Don [45] and Myshetskoe-Dolgoe Lake [13]. Also pollen-inferred Holocene paleotemperature data from Nikokolsko-Lutinskoye bog, Shirinsky Mokh bog, Lammin-Suo bog, Vishnevskoye Lake and Sakkala bog [46] are available. Chironomid-based mean July air temperature reconstructions are available from Nikolay Lake [47], Lake Medvedevskoe [48] and Lake Glubokoye [49,50]. One of the important regions in terms of East European paleoclimatology is the Western Dvina Lakeland (Western Russia), where the terrain relief was formed after the recession of the Valdai Ice Sheet [51]. This region was considered crucial for studying the development of human settlement [52][53][54][55]. The Neolithic occupation in this area began at ca. 8.3 ka cal BP, and the development of Neolithic communities lasted up to around 3.8 ka cal BP [53,54,56,57]. The pile-dwelling settlement existed in the region during the Middle and Late Neolithic period (between 5.9 and 4.2 ka cal BP) [53,58]. The existence of Neolithic communities in this part of Europe was based on a hunter-gatherer subsistence strategy. The elements of agriculture were adopted gradually only in the Late Neolithic (from ca. 3.8 ka cal BP) [59][60][61][62]. However, the local production of biocenoses was influenced by climate variability, particularly mean July air temperature (T Jul), precipitation, and duration of the growing season. Hunter-gatherer communities were also dependent on climate conditions [55]. Thus, the sudden climate fluctuations linked to the Bond events, during 5.9 and 4.2 ka cal BP, might have triggered cultural changes in the Serteya region, forcing a lifestyle change due to the lower productivity of lake and terrestrial ecosystems, lake level fluctuations, and changes in snow cover.
Bond events are climate fluctuations in the Holocene which are based on petrological traces of drifting ice in the North Atlantic over a 1000-year cycle [63,64]. The possible causes for these events include: (1) Orbital insolation-changes in the Earth's orbit are, however, a process too long to cause a 1000-year fluctuation [65]. Solar irradiance-it is contested at event 4.2, where the peak values of 14 C and 10 Be were very low [66], and is speculated as a possible cause of a decrease in the amount of ultraviolet radiation that may have cooled the ozone layer. These changes could ignite the negative North Atlantic oscillation (NAO) phase, which transported cold and dry air masses over Northern Europe [67,68]. The NAO phenomena in the Holocene correspond chronologically to a decrease in solar activity (increase in 14 C and 10 Be) [63]. NAO is a circulation system that determines mostly autumn and winter weather conditions in Europe. It is the transition between the positive and negative phases which are distinguished by the amount of atmospheric pressure. Based on the observations of the fluctuations connected with global air circulation and oceanic water, it is known that the NAO is characterized by two phases: positive and negative. The positive phase involves the transfer of warm and humid air masses from the Atlantic Ocean in the direction of Northeastern Europe, while in the negative phase warm and humid air masses move toward the Mediterranean area, causing drought due to decreased snowfall in the northeastern part of Europe [69][70][71]. (2) Volcanic activity-it is possible only with huge eruptions, but according to Cronin et al. [72] even the largest ones may cause a reduction of the global temperature by only 0.2-0.3 • C for one to several years. There were 12 volcanic eruptions during the Little Ice Age [73]. (3) Ice-sheet dynamics-the melting of icebergs lowers salinity, while decreasing THC (thermohaline circulation) causes negative feedback leading to ice growth in Laurentide [72,74,75]. This paper presents quantitative reconstructions of the main climate factors (mean summer and winter air temperature, and precipitation) influencing the neolithization processes in Eastern Europe, especially the activities related to human occupation in the Western Dvina Lakeland from the 5.9 ka cal BP event (at 6.0-5.75 ka cal BP) to the 4.2-ka cal BP event (at 4.43-3.97 ka cal BP). The T Jul, lake water level related to effective precipitation, mean temperature of the coldest month (TCM), and annual precipitation (AnP) were inferred from the analyses of Chironomidae, Cladocera, and pollen assemblages [18,25,76] and Nazarova et al. unpublished data). We aimed to identify the scale of temperature fluctuations in the region and climate wetness phases. According to the null hypothesis, during the Northgrippian stage (Middle Holocene), the mean July temperature and climate humidity did not vary significantly in the Eastern European Plain and were not decisive for the late emergence of agriculture in the Serteya area. Alternatively, we hypothesized that significant climate changes, coincident with global events, took place during the mid-Holocene. These changes might influence the long-term development of the Serteya II settlement, e.g.: a switch to the Neolithic culture at 8.2 ka cal BP, the development of the pile-dwelling settlement at 4.2 ka cal BP, the prolonged existence of a hunter-gatherer subsistence strategy, and the late emergence of agriculture after 4.2 ka cal BP [57,77].

Study Area
The Serteya region, including the Serteyka River valley, is situated in Western Russia (the Smolensk Oblast of the Russian Federation), in the vicinity of the European watershed of the Western Dvina (Daugava) River and the Dnieper River Basin (Figure 1). The Serteyka River is a left-bank tributary of the Western Dvina River. According to the physicogeographical division, the Serteya region belongs to the Western Dvina Lakeland [78] or the Vitebsk Lakeland [79]. It lies in the temperate continental climate zone characterized by high-dynamic changes in modern meteorological conditions which strongly influence the natural environment of the area [80]. Winters are moderately cold, with a stable snow cover from November to April, and summers are moderately cool and humid. Much of the AnP (ca. 60%) falls between the months of May and October. The climate data of Velizh, situated 20 km west of the site, has an average annual temperature of 3.6 • C to 8.4 • C and average annual atmospheric precipitation varying from 488 to 1296 mm per year for the years 1955-2017. The mean temperature of the warmest month (TWM) (July) varies between 14.8 • C and 23.0 • C, while that of the coldest month (January) ranges from −17.7 • C to −0.8 • C [81]. tively, we hypothesized that significant climate changes, coincident with global events, took place during the mid-Holocene. These changes might influence the long-term development of the Serteya II settlement, e.g.: a switch to the Neolithic culture at 8.2 ka cal BP, the development of the pile-dwelling settlement at 4.2 ka cal BP, the prolonged existence of a hunter-gatherer subsistence strategy, and the late emergence of agriculture after 4.2 ka cal BP [57,77].

Study Area
The Serteya region, including the Serteyka River valley, is situated in Western Russia (the Smolensk Oblast of the Russian Federation), in the vicinity of the European watershed of the Western Dvina (Daugava) River and the Dnieper River Basin (Figure 1). The Serteyka River is a left-bank tributary of the Western Dvina River. According to the physico-geographical division, the Serteya region belongs to the Western Dvina Lakeland [78] or the Vitebsk Lakeland [79]. It lies in the temperate continental climate zone characterized by high-dynamic changes in modern meteorological conditions which strongly influence the natural environment of the area [80]. Winters are moderately cold, with a stable snow cover from November to April, and summers are moderately cool and humid. Much of the AnP (ca. 60%) falls between the months of May and October. The climate data of Velizh, situated 20 km west of the site, has an average annual temperature of 3.6 °C to 8.4 °C and average annual atmospheric precipitation varying from 488 to 1296 mm per year for the years 1955-2017. The mean temperature of the warmest month (TWM) (July) varies between 14.8 °C and 23.0 °C, while that of the coldest month (January) ranges from −17.7 °C to −0.8 °C [81]. At the end of the Valdai Glaciation, geomorphological processes led to the formation of the land relief framework in the area [82,83]. Within the tunnel valley, which is occupied in modern times by the lower Serteyka River valley, a few lake basins existed in the Late Vistulian and during the Holocene. Previous reports indicate that, in the Late Holocene, the lake landscape was gradually replaced by a fluvial system due to the headward erosion and subsequent draining of the lake basins during the development of the Serteyka River valley [84,85]. This scheme is similar to that described by Błaszkiewicz [86,87] for the Kociewskie Lakeland, Poland. At present, the Great Serteya Palaeolake Basin is the largest extension of the valley, with a length of ca. 1.2 km and a width of over 0.5 km [85].
During the last 30 years of intense archaeological survey, about 60 archaeological sites, originating from different periods since the Late Paleolithic to the Middle Ages, have been discovered in the surrounding area. The results of the long-term archeological research formed the basis for the reconstruction of the human settlement development At the end of the Valdai Glaciation, geomorphological processes led to the formation of the land relief framework in the area [82,83]. Within the tunnel valley, which is occupied in modern times by the lower Serteyka River valley, a few lake basins existed in the Late Vistulian and during the Holocene. Previous reports indicate that, in the Late Holocene, the lake landscape was gradually replaced by a fluvial system due to the headward erosion and subsequent draining of the lake basins during the development of the Serteyka River valley [84,85]. This scheme is similar to that described by Błaszkiewicz [86,87] for the Kociewskie Lakeland, Poland. At present, the Great Serteya Palaeolake Basin is the largest extension of the valley, with a length of ca. 1.2 km and a width of over 0.5 km [85].
During the last 30 years of intense archaeological survey, about 60 archaeological sites, originating from different periods since the Late Paleolithic to the Middle Ages, have been discovered in the surrounding area. The results of the long-term archeological research formed the basis for the reconstruction of the human settlement development and economic foundations of communities in the subsequent periods. The human settlement in this area was strongly affected by paleoenvironmental conditions. Previous palynological, diatom, and geochemical research of biogenic sediment cores M25 from the Serteya region provided an insight into the general paleoclimatic situation of the study area. Kittel et al. [57] explained the paleoecological changes of the Great Serteya Palaeolake Basin in the last 6.3 ky BP. Based on previous research of the Neolithic environmental con-ditions between ca. 4300 and 1600 cal BC, the period 4150-3250 cal BC is distinguished as a time of increased settlement activity of Neolithic hunter-gatherers on the paleolake shore. On the basis of ecological proxies, fluctuations in the lake water level and the increase in the lake trophic state caused by human influence were also reconstructed. The results show that the Neolithic community and its functioning were strongly influenced by changing environmental conditions, including climate change. Łuców et al. [88] described the last 300-year ecological history of the Serteya region, based on analysis of the organic sediments of an ombrotrophic fen, which revealed the influence of local anthropogenic disturbances on the functioning of the fen. This work reconstructed climate and lake level change based on three proxies (testate amoeba and pollen, plant macrofossils). In the present study, a complete chronology of the sediment was obtained from AMS radiocarbon dating based on the macrofossils of terrestrial plants [cf. 57]. The chronology used follows Kittel et al. [57]. The results presented below are a further analysis of the data obtained from the deposits of the STII M25 core.

Archeological Investigation
Two main areas of the Serteya II have been explored so far [52,57]. The eastern part of the site contains the remains of the pile-dwelling structures, which are dated to 5-4 ka cal BP. The western part of the site was first inhabited in the Late Mesolithic period. During the Early Neolithic, periodic settlements were established at 9-6 ka cal BP. At 6-5 ka cal BP, people of Eneolithic steppe cultures were active in the area. Numerous wooden and bone objects, remnants of residential structures, and fishing devices were found at these complexes. In addition, finds from ca. 5-4 ka cal BP belonging to the Usviaty Culture and Zhizhitsa Culture were recorded, similar to materials from pile dwellings found in the eastern part of the site. At ca. 4.5 ka cal BP, the bodies of two young women were deposited in the paleolake shore zone of site's western part [52,57,89,90].

Field Studies and Coring
The geological and geomorphological survey made in 2013-2019 enabled a detailed reconstruction of the geology of the Serteya II site and its vicinity [52,85]. In 2016, a 1.4-m core (STII M25) was taken from the wall of the archeological outcrop in close proximity to the human bones discovered in 2015 [57] covering deposits from 23 to 160 cm below ground level (b.g.l.). The core was placed in three metal boxes measuring 50 cm × 10 cm × 10 cm in size. It consists of the following deposits [57] The age of the deposits was determined by Kittel et al. [57] from selected macrofossils of terrestrial plants through 14C AMS (accelerator mass spectrometry) analysis. The calibration and chronology were estimated using OxCal v. 4.4.3 with IntCal 20 calibration curve and P_Sequence model [91,92]. For the 1.4-m-long core, a total of six AMS dates, from depths 34, 85, 116, 138, 148, and 158 cm b.g.l., were used for the construction of the depth-age model (not shown) [57].

Chironomidae-Laboratory Techniques, Identification, and Statistics
The subfossil Chironomidae assemblage analysis was carried out on the STII M25 core with a resolution of 2 cm. For this purpose, the sediment samples were sieved on 90 µm sieves. The laboratory methods described by Brooks et al. [18] were used for the analysis. In total, 3503 specimens were identified and classified into three Chironomidae subfamilies and undetermined Ceratopogonidae. The subfossils were identified to the lowest possible taxonomic level, mainly using the keys of Brooks et al. [18] and Andersen et al. [93]. The sequence zonation was determined using the optimal sum-of-squares partitioning method described by Birks [94], Birks and Gordon [95], and Bennett [96]. The statistically significant number of zones was calculated based on the broken-stick model [97]. The number of zones was determined using ZONE software [98] and BSTICK [99], and the stratigraphic diagram was made using C2 software [100].

Cladocera-Laboratory Techniques, Identification, and Statistics
The sediments the Serteya STII M25 core were examined at 2-cm intervals. Each 1-cm 3 sample of the deposit was processed according to the standard procedure [101]. A minimum of 200 cladoceran remains per sample were identified. The taxonomy follows Szeroczyńska and Sarmaja-Korjonen [102] and Van Damm et al. [103]. The ecological preferences of cladoceran taxa were determined using Bjerring et al. [104]. The cladoceran zones were determined by a stratigraphically constrained cluster analysis (CONISS) using POLPAL software [105].

Pollen-Laboratory Techniques, Identification, and Statistics
A total of 55 one cm 3 samples (2-cm interval) were treated following the standard pollen preparation procedures (treated with HCl and NaOH, sieved using a 250-µm sieve, treated with HF, acetolysis) [106]. Pollen and NPPs were identified at a magnification of ×400 and ×630 with reference to published illustrations, morphological keys-including those of Faegri and Iversen [106], Reille [107], van Geel [108], van Geel and Aptroot [109], and Beug [110], and a laboratory reference collection. Pollen percentage and influx diagrams were constructed using Tilia software [111,112]. Aquatic and hygrophilous plants, undetermined pollen, ferns, moss, and NPPs were excluded from the pollen sum. The percentages of excluded taxa were calculated using the total pollen sum [113].

Temperature, Precipitation, and Water Depth Reconstructions
Climate parameters are inferred from chironomids, cladocerans, and pollen assemblages, and water level fluctuations are estimated using chironomids and cladoceran data. The characteristics of the training sets used are presented in Table 1. The pollen-inferred (P-I) reconstructions of temperature and precipitation at Serteya across the Holocene were performed using the Modern Analogue Technique (MAT) [116]. This method has been utilized at a global scale for samples from Europe [117,118] and at a more regional scale for samples from Mediterranean pollen sequences [119,120]. The MAT is not a transfer function (cf. the WA-PLS) but an "assemblage approach," and is based on the primary assumption that pollen samples sharing a similar composition are a by-product of comparable vegetation assemblages. Using MAT, the composition and abundance of fossil pollen samples are compared with those of modern pollen surface samples, and the similarity between each fossil sample and modern pollen assemblage is evaluated by a squared chord distance.
Once a set of modern pollen samples has been selected as analogues (samples with the smallest distance), climate parameters are assigned to each fossil sample as the weighted average of the climate parameters of the modern samples which are the best analogues. In this study, the number of analogues chosen was five, which were selected using a leave-oneout cross-validation test. The training set of modern pollen samples contains more than 3000 samples from Europe and the Mediterranean area [76]. We also calculated and plotted the chord distance of the first best analogue (Distmin1) and the last analogue selected (Distmin2) (not illustrated) based on a threshold defined by a Monte Carlo method, in order to check the accuracy of the pollen-based climatic reconstruction. The minimum distance (Distmin1) calculated was low for P-I TCM, P-I AnP, indicating that the analogues selected are consistent. Based on previous studies, we reconstructed the climate for precipitation and mean TCM so that the results would be comparable across the region. These reconstructions were carried out using the package Rioja with R software [121].
A Detrended Correspondence Analysis (DCA) was carried out on the percentage data for the Chironomidae, Cladocera, and pollen datasets, with detrending by segments and down-weighting of rare species using Canoco 4.5 [122]. The temperature, water depth, and precipitation reconstructions used C2 and NeuroGenetic Optimizer software (NGO, version 2.6.130, © BioComp Systems, Inc.) [100,123].

Chironomidae, Cladocera, and Pollen Assemblage Stratigraphy
The detailed ecological description of the zones can be found in Kittel et al. [54], here we present only a summary of the ecological interpretation and diagrams ( Figure S1). Based on the similarity analysis using CONISS and BSTICK [99,124], the zones are divided as follows for each proxy. The Chironomidae STII M25 sequence is divided into five significant zones ( Figure S1

Detrended Correspondence Analysis
The first DCA axis of the Chironomidae sequence explained 23.6% of the cumulative variance of species data and had a gradient length of 1.742, proving the linear distribution of the data. The second DCA axis explained 7.7% of the cumulative variance of species data. The DCA Ax1 of chironomid assemblages revealed a similar trend to the estimations of mean July air temperature from the Fn TS and SNP TS WA-PLS ( Figure 5), which indicates that the summer temperature was the main driver of the midge communities. Furthermore, the estimations of mean July air temperature from the SNP TS artificial neural network (ANN) were similar to DCA Ax2 of Chironomidae assemblages, but it is a completely different trend than that marked on Ax1. The reconstruction based on the SNP of ANN reflects the influence of another unknown factor. Regarding Cladocera communities, the first DCA axis explained 27.3% of the cumulative variance of species data, with a gradient length of 1.339, proving the linear distribution of the data. The second DCA axis explained 9.9% of the cumulative variance. Similar to chironomid assemblages, the DCA Ax1 of cladoceran communities also reflected the mean July air temperature estimations from the Fn TS ( Figure 5), indicating that the summer temperature was the main driver of the water flea communities, whereas Ax2 revealed a different trend.

The Chironomidae-Inferred and Cladocera-Inferred Mean Summer Air Temperature Reconstructions
According to the CH-I T Jul of the 31 Fn TS samples, seven samples remained below the 2 percentile threshold (minDC < 7.40328), representing that they are samples with very good modern analogues, 14 remained below the 5 percentile threshold (minDC < 8.85913), representing good modern analogues, and only two had values over the 10 percentile threshold (minDC > 10.1025), representing poor or very poor modern analogues. Referring to the SNP TS WA-PLS reconstruction, among the 31 samples, only one was below the 2 percentile threshold (minDC < 6.98983) and three were below the 5 percentile threshold (minDC < 8.57757), representing that they are samples with very good and good modern analogues, respectively. As many as 16   AnP. The cold-winter arid phase Ic fell on the 5.9-ka cal BP event and phase IIIc on the 4.2-ka cal BP event. Phase IIc coincided with a distinct rise of CH-I and CL-I summer temperature reconstructions derived from all three TSs. The minimum distance (Distmin1) calculated was low, indicating that the selected analogues are consistent and that our reconstruction is acceptable considering the possible bias due to the human impact on the pollen data. However, no indicator of farming activities was identified in the pollen diagram [57].

The Cladocera-Inferred and Chironomidae-Inferred Lake Water Level Reconstructions
All the samples of CL-I water depth reconstruction (CL-I depth) represented very poor modern analogues, which were over the 20 percentile threshold (minDC > 3.68913). According to the CH-I water depth reconstruction (CH-I depth), only one sample from 124 cm remained below the 2 percentile threshold (minDC < 8.57652), representing very good modern analogues, and six remained within the 5-10 percentile threshold (9.81504 < minDC < 10.8182), representing moderate modern analogues. All the other 24 samples remained over the 10 percentile threshold (minDC > 10.8182), representing that they are samples with poor and very poor modern analogues.
The Fn TS CL-I depth reconstruction showed a constantly shallow littoral depth ranging from 1.2 to 2 m. However, this reconstruction indicated a relative increase in water level at 142 cm/ca. 5.9 ka cal BP (from 1.

Discussion
The CH-I T Jul reconstructions were based on three different models. As the training sets originated from different regions-Western and Central Europe (SNP TS), Central and Northern Europe (Fn TS), and Eastern Eurasia (Rn TS), we obtained different results. Nonetheless, the reconstructions were consistent, giving similar values and showing more or less similar trends. It is not surprising that the SNP TS and Fn TS showed similar trends compared to the more dissimilar Rn TS, because they shared ca. 100 sites from Poland. The SNP TS revealed high temperature variability because it combines three different climatic regions representing separate midge fauna, with a wide range of temperatures. The Fn TS gave the closest modern analogues as it uses samples from the nearby geographical region. The results of the Fn TS reconstruction are strongly flattened. The Rn TS gave the highest temperature values, in the range of continental climates of Central Siberia. This reconstruction indicates an overall slightly decreasing trend, but is erratic. The DCA and N2 showed similar trends to the SNP TS reconstruction indicating that the summer temperature was a dominant driver for Chironomidae communities. The cold oscillations took place at about 5.9 and 4.2 ka. Generally, these oscillations are weak and are apparent in a few consecutive samples (1)(2)(3)(4) and indicate a drop in T Jul of ca. 1-2 • C. This exceeds the RMSEP of the models, which suggests this is not noise. It can be seen from the reconstructions that during the Northgrippian stage the summer air temperature was high, mostly in the range of 17-20 • C. These values are similar to the current July air temperatures in the Serteya region [80,125]. A comparison with other sites in Eastern Europe indicates that the values of P-I TWM and CH-I T Jul were usually 2 • C lower in the circum-Baltic region, whereas on the Central Polish Plain (e.g.,Żabieniec, Błędowo) and the Central Russian Plain (Staroselsky Moch, Peatland Klukva) the values were similar to those at the Serteya STII M25 (Table S1, Figure S2). The data from the Leningrad Oblast shows a cold oscillation between 5.9 and 4.2 ka cal BP [43]. This result indicates that in more continental regions, distant from the Baltic influence, the summer temperatures were slightly higher.
The SNP TS reconstructions showed a increasing temperature trend from 5.9 to 4.2 ka cal BP, punctuated by two cool oscillations indicated at 5.9 and 4.2 ka. However, the Rn TS showed a weakly decreasing trend. An increasing trend is also indicated from Upper Don [45]. The global temperature trend of terrestrial records from 30 • to 60 • latitude is decreasing after 6 ka cal BP [126,127]. In addition, the simulations of Renssen et al. [128] indicated that the Holocene Thermal Maximum (HTM) occurred in Eastern Europe earlier than 6 ka cal BP. A comparison of the temperature of both cooling events in Table S1 indicated a similar trend in the East European Lowland. At some sites in the southern part of the region, the temperatures were slightly higher at 5.9 ka cal BP. However, in the northern part, some sites had slightly higher TWM during 4.2 ka cal BP. Warden et al. [3] clearly indicated that the HTM in the pari-Baltic region occurred at ca. 5.9-5.5 ka cal BP, with a cooling trend to 4.2 ka cal BP. The nearby P-I TWM reconstruction from Staroselsky Moch revealed more or less the same values through 6-4 ka cal BP. Tarasov et al. [62] also suggested that in the Western Dvina Lakeland HTM occurred in ca 8.6-6.9 ka cal BP and in the Karelian Isthmus in ca 7.8 cal BP [40]. It started and ended some centuries earlier than suggested for the neighboring Baltic region-ca 8.1-5.6 ka cal BP. Based on pollen data, Tarasov et al. [62] reported that the climate in Middle Holocene was dryer and warmer, whereas after 5.4 ka cal BP the climate was cooler and wetter. Nonetheless, the authors did not perform quantitative reconstructions from the investigated sites, but their results seem to be more consistent with the Rn TS reconstruction.
The SNP TS and Fn TS revealed a cold oscillation from the Mid-to Late Northgrippian (5.9-4.2 ka cal BP) stage. It may be a microclimate pattern that is dependent on local factors such as lake transgressions. The warm-stenotherm chironomid taxa usually have preferences for a higher trophic state [17]. Human-induced eutrophication might have overridden the influence of July air temperature resulting in higher reconstructed temperatures [57].
Juggins [129] suggested that only one quantitative reconstruction should be obtained from a particular proxy. Here, we decided to illustrate both reconstructions from the pollen data. Both revealed very similar trends and were consistent with Distmin1 ( Figure 5). This suggests that plant communities were influenced by one climatic factor: continentality fluctuations. Such a coincidence of decreased precipitation and TCM indicates that continentality is the main driving factor of plant communities in the region. At about 4.9 ka occurs sudden drop in precipitation and TCM. The values of TCM, TWM/T Jul, and AnP were compared for a wide group of 22 sites from Eastern Europe (Table S1), no clear signal of the 5.9-and 4.2-ka cal BP oscillations in the region was noted.
Knowledge about the cold episode at ca 5.9 ka BP is still insufficient. There is no consensus on the magnitude and timing of the 5.9 Bond event oscillation in eastern Europe [36,[128][129][130][131]. This raises a doubt whether there are any well-documented quantitative reconstructions of these Bond events in Eastern Europe. Where there are some records indicating oscillations in a few consecutive samples, the signal is weak and unclear. Hence, climatic records of much better resolution are needed from the region to verify the hypothesis of the presence of these Bond events.
Because the Bond events are linked to the NAO fluctuations, they might be more significant in Western Europe than in the East. The European territory along the North Atlantic is the main area to be affected by this climatic phenomenon. The STII M25 site is located in the Western Dvina Lakeland in the contact zone of transitional and continental climate, which could have made it sensitive to the paleoclimatic regime changes (NAO circulation). Based on their studies, Yao and Luo [130] and Rousi et al. [131] concluded that the eastern and western types of NAO differently influence the temperature and precipitation regimes in Europe. To study the magnitude of the NAO influence and the relationship between NAO and the temperate climate zone, the domain size must also be considered [132]. For this zone, the optimal domain sizes are lower in summer compared to winter, and for continental rather than oceanic regions. The NAO will likely have a stronger impact on the climate of Northeast Europe in winter in continental regions. From the results ( Figure 5), one can speculate that the positive NAO phases during the 4.2 and 5.9 ka events was characterized by cool summers and severely cold winters, as well as intense precipitation mostly in autumn. The differences between the dynamics of precipitation during these events may result from changes in the activity of ultraviolet radiation and glacier size. For further conclusions, it is necessary to study more sites in Eastern Europe at high temporal resolution. However, for the LIA, correlated with the Bond Event no. 0, a predominance of negative NAO phases is postulated by Trouet [133]. Moreover, a positive NAO phase is postulated for the MCA-LIA period [134,135]. It must be also stressed that Moreno-Chamarro et al. [136] emphasis a reduction in the heat transport by the subpolar gyre rather than a NAO impact on the LIA cooling, mostly during the winter seasons.
The boundary of boreal plant assemblages and temperate deciduous forest in the studied area also reflect the climate zones. Thus, P-I reconstructions should reflect the change in the range of vegetation composition in Eastern Europe. Pleskot et al. [137] questions the presence of the 4.2-ka cal BP oscillation in Poland. At the STII M25 P-I, the higher continentality signal is clear at both events and has also been recorded in a few subsequent samples. Pleskot et al. [137] performed only CH-I T Jul reconstruction. The STII M25 CH-I T Jul reconstruction also gave a weak and ambiguous signal of 4.2-and 5.9-ka cal BP oscillations. P-I reconstructions revealed a much clearer pattern at the STII M25 site, but not at the nearby Staroselsky Moch and Klukva Peatland [138,139].
The P-I TCM/AnP and CH-I T Jul reconstructions from STII M25 indicated the third intermediate oscillation falling on 5.0-4.7 ka cal BP. This oscillation manifests by a small increase of T Jul at 5.1 ka and distinct increase in continentality. Its maximum effect falls on 4.9-4.8 ka cal BP, and its magnitude is comparable with the 5.9-ka cal BP signal and stronger than the 4.2-ka cal BP one. Modern maps indicate that this microregion differs from the wider Smolensk region. The oscillation may have had a more local character as there are warmer and colder microregions within the Western Dvina Lakeland. Besides P-I and CH-I, the present study performed CL-I summer temperature and lake depth reconstructions. Although water level fluctuations are regarded as the primary factor influencing Cladocera communities in the long term [28], the DCA ( Figure 5) results were consistent with the CL-I summer temperature reconstruction, indicating it as the primary driver of cladoceran communities. The values of temperature are lower than CH-I T Jul because midges are more linked to air temperature while water fleas are linked to water temperature [6,24,140,141]. Thus, the lower values of CL-I temperature are not surprising; however, the decrease of CL-I T Jul also differs from CH-I. From 136 to 104 cm b.g.l., the temperature increasing was observed, similar to the SNP and Fn TS CH-I T Jul reconstructions, although the earlier (158-134 cm b.g.l.) oscillation was opposite. In addition, the reconstruction from the three topmost samples revealed a dissimilar pattern to all CH-I reconstructions. This can also be explained by the differences and complexity of proxy responses to temperature which merges with water trophic state and other environmental factors [142]. Moreover, it needs to be mentioned that CL-I reconstructions from STII M25 gave poor and very poor modern analogues.
Finally, the question of CH-I water level reconstruction remains. Because the site is located in the lakeshore zone [57,85], water level fluctuations at STII M25 are crucial for explaining the environmental causes of Neolithic settlements. The results showed by CL-I and CH-I lake depth reconstructions are consistent with general observations that STII M25 is located in the shallow littoral zone [57]. The CH-I lake depth suggests that the water level was unstable and may have even been seasonal before the 5.9-ka cal BP event, and then there were very shallow water conditions up to ca. 4.5 ka cal BP. According to the CH-I lake depth, the bodies of women recovered at ca. 115 cm b.g.l. should have been placed when the water level was increasing before the 4.2-ka cal BP event. During this event, the lake revealed a higher water level than during the 5.9-ka cal BP event. Nonetheless, the CH-I depth gave generally better modern analogues between the 5.9-and 4.2-ka cal BP event in comparison to the CL-I depth, and the fluctuations of both reconstructions remained within the RMSEP range. The mean July temperature was the leading factor in the case of both Chironomidae and Cladocera, overriding water level fluctuations in the permanently shallow littoral zone. Therefore, the constructed quantitative reconstructions of water level fluctuations are doubtful, and in this case, qualitative reconstructions of the habitat change presented in a previous study [57] are more valid.
The above-presented reconstructions prove that climate conditions during the Northgrippian stage were suitable for human communities. The local microclimate in the Serteya could have been buffered by the lake-river system, weakening the fluctuations of summer temperature during 5.9 and 4.2 ka cal BP. Summer temperature was relatively high during the whole Northgrippian stage, due to which tribes from the Black Sea Lowland settled here after 8.2 ka cal BP [1,55]. Stable climate conditions would imply constant primary and secondary production of an ecosystem supporting a hunter-gatherer economy for a long time after agricultural development in localities in the West and the North [3]. However, plant communities here consist of those that clearly reacted negatively to the 5.9-and 4.2-ka cal BP event. This suggests that the Bond events might have indeed influenced the ecology of local settlers triggering a change in their culture to altered, more difficult (higher continentality) conditions. The intermediate climate oscillation falling on 5.0-4.7 ka cal BP coincides with the appearance of relicts of the Globular Amphora culture at the Serteyka River Valley. It indicates that local settlers participated in spreading trails of amber to the East. Hence, the Western Dvina Lakeland must have had high economic status in the region [57,85] which can be attributed to productive local ecosystems, high geodiversity, and favorable [142] climate conditions.

Conclusions
In recent years, quantitative reconstructions have impacted paleoecology mainly by providing precise quantitative records, though they have their own limitations. One of them is the reconstructed chronology of non-laminated lake sediments, another is sedimentation rate in lakes, and the lag in the reaction of plant biota to climate change that limits reconstruction to decadal resolution [143]. Moreover, the methods used for quantitative reconstructions have their own statistical limitations which can be accessed from the R2, RMSEP, and modern analogues. The climatic records may be obscured by local impacts from non-climatic influences on the proxies [129,142]. The present study has included several reconstructions using three different proxies: chironomids, cladocerans, and pollen. Summer mean air temperature appeared to have a stronger influence than the lake depth on the midges and water fleas, which is in line with the general observations of their ecology [144,145]. This indicates that local water level fluctuations had only a secondary effect on the littoral lake ecosystem. The CH-I mean air temperature and CL-I mean water temperature of the summer season are similar to the recent, Late Holocene values. Although STII M25 is located in the western part of the archaeological site, not within a major dwelling area of the pile-dwelling settlement, human-induced local eutrophication might have still reinforced the effect of temperature, to produce a false signal of temperature increase inferred from Chironomidae as well as Cladocera. The pollen analysis indicated an increase in continentality during the investigated period. Furthermore, P-I TCM and AnP showed clearer oscillations at the 5.9-and 4.2-ka cal BP events than the chironomids and Cladocera. Hence, the P-I reconstructions may be considered as the primary result of the presented studies. It is difficult to recognize mid-Holocene (Northgrippian) oscillations, especially in Eastern Europe where the influence of NAO is weaker [24,69,137]. In order to speculate about their time shift (5.9 ka cal BP and 4.2 cal BP), it would be useful to increase the number of quantitative reconstructions for the area of Eastern Europe. The comparison of the results obtained with broad-spectrum quantitative reconstructions from the region (Table S1) shows that the SNP TS results do agree with other global reconstructions. Most reconstructions from Eastern Europe have insufficient resolution to capture Bond oscillations in the Northgrippian stage. Considering all the methodological cautions in the interpretation of the presented results, we are inclined to support the alternative hypothesis presented in the Introduction, though with some reservations. Nevertheless, we cannot prove that the oscillations are at a regional or global scale and result from NAO. We assume that these climate changes might have influenced the progression of the Serteya II settlement over time, including the possible collapse of the pile-dwelling settlement at 4.2 ka cal BP and the prolonged existence of a hunter-fisher-gatherer subsistence strategy as the conditions were favorable for the productivity of local ecosystems.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.