Late Quaternary Climate Reconstruction and Lead-Lag Relationships of Biotic and Sediment-Geochemical Indicators at Lake Bolshoe Toko, Siberia

Millennial-scale climate change history in eastern Siberia and relationships between diatom diversity, paleoclimate, and sediment-geochemical lake system trajectories are still poorly understood. This study investigates multi-proxy time series reaching back to the Late Pleistocene derived from radiocarbon dated Lake Bolshoe Toko sediment cores, southeastern Yakutia, Russia. We analyzed diatoms, elements (XRF), minerals (XRD), grain-size, organic carbon, and included chironomid analyses and published pollen-data for quantitative paleoclimate reconstruction. Changes in diatom species abundances reveal repeated episodes of thermal stratification indicated by shifts from euplanktonic Aulacoseira to Cyclotella species. Chironomid and pollen-inferred temperature reconstruction reveal that the main shift between these diatom species is related to the onset of Holocene Thermal Maximum (HTM) at 7.1 cal ka BP. Comparison to other paleoclimate records along a north-south transect through Yakutia shows that the HTM was delayed as far south as the Stanovoy mountains. Relationships between sediment-geochemistry, paleoclimate variability and diatom species richness (alpha diversity) was tested in a moving temporal offset approach to detect lead-lag relationships. Sediment-geochemical data, mainly uniform during the Holocene, revealed strongest positive or negative correlations ahead of species richness changes. Mean July air temperature (T_July) reconstructions correlate with both Hill numbers and relative assemblage changes indicated by sample scores of multidimensional scaling analysis (MDS) over the entire time series. We found that sediment organic carbon revealed distinct positive correlations, i.e., centennial-scale delay to increases in diatom effective richness (Hill numbers N0 and N2). We conclude that a lag of deposited organic carbon concentrations behind changes in diatom alpha diversity reveals that species richness can augment the production and thus sequestration of organic matter in comparable lake systems.

Introduction

Over the past few decades, boreal and high elevation regions warmed faster than elsewhere (Mountain-Research-Initiative et al., 2015) and was strongest in northern Russia (Biskaborn B. K. et al., 2019). Although the Siberian region is generally regarded as being underrepresented in paleoenvironmental research (Sundqvist et al., 2014; Mckay et al., 2018; Kaufman et al., 2020), there is an ever-expanding interest in the area due to expected environmental changes and associated strong socio-environmental implications (Amap, 2017). In particular, remote Arctic lakes in permafrost regions represent pristine early-warning systems of environmental change, as they rapidly respond to climatic perturbations and changes in geo-ecological boundary conditions in their catchments (Subetto et al., 2017; Biskaborn et al., 2021; Nazarova L. B. et al., 2021). In terms of environmental sustainability, they still provide an infinite resource of uncontaminated fresh water, which make them unique compared to lake systems in the highly populated mid-latitude areas (Serreze et al., 2006). In order to assess the dimensions of modern landscape changes in Siberia, information about ecosystem changes during comparable climate warming episodes of the past provide valuable clues. In the terrestrial part of Siberia, such information, for instance, is stored in lake sedimentary records. The last strong warming of the atmosphere in a comparable setting took place after the Pleistocene-Holocene transition, referred to as the Holocene Thermal Maximum (HTM) that progressed in a non-uniform spatiotemporal event (Kaufman et al., 2004; Renssen et al., 2012). Knowledge about trends within this event was successively gained from paleolimnological records across a north-south transect through Siberia (Biskaborn et al., 2012; Biskaborn et al., 2016). These and our new record rely on radiocarbon dating and sediment-geochemical parameters that provide the background information on the sedimentation history necessary to reliably compile and evaluate paleoenvironmental bioindicator records (Biskaborn et al., 2013b; Wang et al., 2015; Heinecke et al., 2017).

The response of lake ecosystems to climate forcing is well represented by biotic indicators (Cohen, 2003). One of the main lake ecosystem components is given by diatoms, which are unicellular, siliceous microalgae that appear ubiquitous with opaline valves being well preserved in the sedimentary record allowing identification up to sub-species level by high-resolution light microscope analysis (Battarbee et al., 2001). Diatom assemblages (species communities) respond to hydrochemical changes (pH, Oxygen, nutrients) as well as seasonal climate changes via changes in the length of the ice cover, but also to physical habitat changes (water depth, sediment input, turbulence) (Biskaborn et al., 2012; Pestryakova et al., 2012; Biskaborn et al., 2013a; Herzschuh et al., 2013; Hoff et al., 2015; Biskaborn et al., 2016; Palagushkina et al., 2017; Pestryakova et al., 2018).

The simplest straightforward measure to describe diversity in lakes is alpha diversity (species richness). Alpha diversities in wetland ecosystem are affected by past long-term and ongoing recent climate change (Engels et al., 2020). In particular increasing summer temperature at northern sites, as a main driver, leads to increases in diatom diversity (Hobbs et al., 2010). As diatoms represent the main aquatic primary producers that are used as bioindicators in boreal lakes (Smol and Stoermer, 2010), the understanding of diatom diversity is important because it can stabilize primary production and thus overall biomass production (Marzetz et al., 2017). This in turn relevant for the global carbon budget as lakes act as sink for organic matter (Vyse et al., 2021 accepted and close to publication).

Larval head capsules of Chironomidae preserved in lake sediments have proven to be among the best indicators of environmental variations (Battarbee, 2000; Smol et al., 2005). Inference models (Barley et al., 2006; Nazarova et al., 2015; Luoto et al., 2016), which link the present distribution and abundance of chironomids to contemporary climate enable quantitative reconstructions of fundamental environmental parameters, among which are air paleotemperature (Nazarova et al., 2015), water depth (Nazarova et al., 2011a), continentality (Self et al., 2011), late-winter hypolimnetic oxygen (Luoto and Salonen, 2010).

Sedimentary pollen and their modern vegetation analogues have widely been used to reconstruct terrestrial summer temperatures (Biskaborn et al., 2016; Rudaya et al., 2016; Zhilich et al., 2017; Andreev et al., 2021). In our study design, we include a published pollen record from a parallel sediment core (Courtin et al., 2021) as additional independent proxy and counterpart to chironomid based air temperature reconstruction. Sediment-geochemical proxies that were analyzed from Bolshoe Toko surface sediments (XRF-gained elements, XRD-gained minerals, laser-derived grain-size distributions, organic carbon characteristics) are also used in support to biotic indicators. As novelty in context to previous research on Bolshoe Toko our study presents time series of aquatic ecosystem indicators and numerical paleoclimate reconstructions that can be correlated to sediment-geochemical downcore data. The resolution of bioindicator samples allows new insights into environmental variability and relationships up to centennial scale over the Holocene, while in the Pleistocene the age model (due to necessary projection of modern reservoir assumptions) only allows reliable conclusions at millennial scales.

The overall objectives of this paper are to 1) reconstruct millennial-scale climate change history in the yet understudied region of southern Yakutia and, by setting multi-proxy data in a paleoclimate context, to 2) infer relationships between diatom diversity, paleotemperature, and sediment-geochemical lake system trajectories.

Study Site

Lake Bolshoe Toko (56°15′N, 130°30′E, 903 m.a.s.l) was reported as possibly the deepest lake in Yakutia (Zhirkov et al., 2016). The freshwater lake is oligotrophic and circumneutral. Biskaborn B. K. et al. (2019) provided a detailed description of the lake system and reported pH 6.8 and conductivity 35.1 µS cm⁻¹ at the coring site PG2208 (Figure 1). According to their findings, the characteristics of the lake as coring locality for paleolimnological study can be summarized as follows: 1) pristine lake (far from settlements), 2) potentially old (multiple surrounding moraines and fault systems), deep (max. ca. 72.5 m), lake ecosystem represented in sediments (diatoms and chironomids preserved among multiple other proxies).

FIGURE 1

FIGURE 1. Study site of Lake Bolshoe Toko. (A) Overview map. (B) Lake Bolshoe Toko with coring sites. Water depth isobaths based on Konstantinov (2000) slightly modified from Biskaborn B. K. et al. (2019). Moraine cartography and topographic features based on Kornilov (1962).

Today, the sedimentary basin of Lake Bolshoe Toko is 15 km long and 8 km wide, situated at the northern foot of the Stanovoy Mountains representing most of the lake’s catchment (Biskaborn et al., 2021). The Bolshoe Toko basin forms a surficial landform on the southeastern Siberian Platform, which is built up of cratonic basement rocks of Precambrian age and overlain by Paleozoic to Mesozoic sedimentary rock strata and Quaternary cover sediments, in places with intercalated volcanic rock units (Nikishin et al., 2010; Barnet and Steiner, 2021). During late Jurassic to early Cretaceous time, the area formed part of an intramontane depression within the Aldan-Stanovoy Block (Imaeva et al., 2012) filled up with continental sedimentary rocks and coal deposits of the southern Yakutian coalfields (Golovin et al., 2007). Today a fault line separates the mountain range of the Stanovoy Mountains from Lake Bolshoe Toko, which is deepened into the coal deposits of southern Yakutia. Regional mountain glaciation during the Quaternary shaped the lake basin and left behind three generations of morainic arcs and ancient lake terraces of yet unknown ages (Kornilov, 1962; Konstantinov, 2000) (Figure 1).

Modern climate in the study area reveals strong seasonal contrasts under the regime of the Siberian High during winter and low-pressure conditions during the summer months (Mock et al., 1998; Shahgedanova, 2002). Regional climate data are available from the meterological station Toko, situated about 8 km north of Lake Bolshoe Toko (www.worldmeteo.info/en/europe/russia/toko/weather-101939). There, warmest temperatures and highest precipitation occur in July, ranging on average between 6 and 22°C and rain falls around 95 mm per month. Coldest temperatures between −49 and −29°C and low precipitation around 10 mm per month characterize the winter conditions in January. Climate conditions control the extent of lake ice. The study of modern satellite images (https://earthdata.nasa.gov) revealed that annual ice cover usually began in the first week of November and lasted until the end of May to mid-June during the last decade. In earlier times, lake ice cover may have existed longer, as for example observed in the year 1971, when lake ice did not disappear before the mid-July (Konstantinov, 2000).

The study site belongs to the Aldan Floristic Region (Kuznetsova et al., 2010) dominated by northern taiga consisting mainly of Larix accompanied by Picea and Pinus (Courtin et al., 2021).

Materials and Methods

Expedition and Core Retrieval

We conducted field work at Lake Bolshoe Toko in the first week of April 2013 as part of the expedition “Yakutia 2013” as a joint campaign of the North Eastern Federal State University (NEFU, Russia, Yakutsk) and the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI, Germany, Potsdam). Coring equipment was brought by trucks onto the lake. We used a Jiffy ice auger (250 mm) to open holes through the 80 cm thick ice cover as conduits for scientific instruments and coring equipment to the unfrozen water column and bottom sediments. A portable Echo sounder (HONDEX PS-7 LCD) gave estimates of water depth and verified ambiguous values with a calibrated rope. We used UWITEC piston and gravity corer to retrieve sediment cores. Core PG2208 was retrieved by repeated coring of sub-sections at a water depth of 68.3 m. To avoid core loss due to the piston coring technique, parallel 3 m core sections were retrieved with an overlap of 0.5 m. Compaction was reduced by leaving the cores for tension release for several hours after retrieval. Sediment cores were transported in PVC tubes to the AWI laboratories in Potsdam and stored at 4°C in a dark room. Hydrochemical measurements performed during fieldwork are already published in detail in Biskaborn B. K. et al. (2019). The study in hand is based on core PG2208 and additionally involves pollen time series published by Courtin et al. (2021) from a core at 26 m water depth.

X-Ray Fluorescence Scanning

Core sections were opened alongside and cut into halves of which one was used for non-destructive XRF scanning and back-up material, while the other one was subsampled for single sample analyses. We used an AVAATECH XRF Core Scanner at AWI (Bremerhaven), equipped with Rh x-ray tube at 1.2 mA, 12 s count time for 10 kV (no filter) and 1,5 mA, 15 s count time for 30 kV (Pd filter) to scan the core in 10 mm steps for elements ranging from Al to Bi within a He atmosphere. The mean error (expressed as chi square χ²) for elements of interest was Al 2.0, Si 3.5, Ti 1.6, Mn 14.4, Fe 12.1, Br 0.6, Rb 0.7, Sr 1.0.

Core Splicing and Dating

To gain a composite core record for core PG2208 (Figure 1) we used XRF data, i.e., Sr, to correlate overlapping parts of retrieved sections resulting in 367 cm in total. For dating we scanned the sediment core for macro fossils and retrieved bulk sediment samples at depths without any reliable macro remains. A prior attempt to date samples at the Poznan radiocarbon laboratory resulted in a large deviation between humic acid fraction (soluble in NaOH) and the humin fraction (residual, i.e., alkali-insoluble fraction). Therefore, we kept only the Poznan date of the single macro fossil (terrestrial wood remain) that was present and analyzed with standard (reliable) dating routines. For creating a reliable model, we dated another set of bulk sediment samples at MICADAS radiocarbon laboratory at AWI using an acid treatment method outlined in Vyse et al. (2020). All dating results are shown in detail in Courtin et al. (2021) including the age-model for PG2133 used for pollen in their study.

To estimate a reservoir (old carbon) effect we dated an undisturbed surface sample PG2119 that was 1.6 km SW of PG2208. To estimate recent accumulation rates, we developed chronologies for a short core PG2203 that was 935 m SE from coring location PG2208 using the Lead-210 (²¹⁰Pb) and ¹³⁷Cs technique by direct gamma assay in the Liverpool University Environmental Radioactivity Laboratory, using Ortec HPGe GWL series well-type coaxial low background intrinsic germanium detectors (Appleby et al., 1986). Details on the methods including analyses of short cores from Bolshoe Toko are provided by Biskaborn et al. (2021). All MICADAS datings performed on core PG2208, the Poznan dating of the plant remain, and the Lead-210 samples are shown in Table 1. We applied the R package “Bacon” (Blaauw and Christen, 2011) and the IntCal20 calibration curve (Reimer et al., 2020) to model the age-depth relationship based on 11 radiocarbon dates from PG2208 and nine calibrated ages (Lead-210 and ¹³⁷Cs technique) from the short core PG2203 reported onto composite depths from PG2208 (Figure 2). We used a reservoir value of 762 ± 66 years based on PG2119 to correct bulk sample radiocarbon estimates (Table 1).

TABLE 1

TABLE 1. List of samples and dating methods used for age-depth modelling. All dates are converted to C14 ages relative to 1950 CE. All dating results from MICADAS and Liverpool are shown. Just one Poznan date on wood residuals (Courtin et al., 2021) was used in this study. The composite depth indicated as “Com. Depth (cm)” represents the total depth below the sediment surface; the specific depth (cm) within core sections is indicated in “Core sample ID”.

FIGURE 2

FIGURE 2. Bacon age-depth model based on radiocarbon data from sediment core PG2208 (367 cm) and ²¹⁰Pb and ¹³⁷Cs dating from the uppermost 4 cm of short core PG2203.

Grain Size

Before conducting grain-size analysis, we treated 36 samples with hydrogen peroxide (30%), set on a platform shaker for 4–5 weeks to remove organic matter, while keeping the pH neutral using ammonia and acetic acid. After washing, sediment samples were treated with sodium pyrophosphate for dispersion and homogenization of the sediment suspension for a day. At least three measurements per sample were carried out using a laser diffraction particle size analyzer MASTERSIZER 3000. Their mean was used to estimate sand, silt, and clay fractions as sums of volume percentages between 2 and 63 mm, 63 and 2 mm, and <2 mm grain-size classes, respectively.

Carbon and Nitrogen

To gain information on organic matter variability we analyzed 204 samples for total carbon (TC) and nitrogen (N), 190 samples for total organic carbon (TOC), 45 samples for stable carbon isotopes (δ¹³C), and determined the water content in 417 samples. We used an Elementar Vario EL III to estimate TC and N from two repeated measurements of each freeze-dried and milled sample. The measurement accuracy is 0.05% for TC and 0.1% for TN. We used a Vario MAX C analyser to determine TOC. Measurement accuracy was 0.1%. We analyzed δ¹³C from carbonate-freed samples using a Finnigan Delta-S mass spectrometer. Results are presented as δ¹³C values relative to the PDB reference standard (‰) with an error of ±0.15%.

X-Ray Diffractometry

To analyze the mineralogical composition of sediments we freeze-dried and milled 19 samples for X-ray diffractometry (XRD) using a PHILIPS PW1820 goniometer (40 kV, 40 mA, from 3 to 100°, step-rate 0.05°, Co ka radiation). Data processing was performed using the free software MacDiff 4.0.7. To avoid influence of opal and organic matter in the samples (Voigt, 2009), we used the ratio of the peak area intensities of each mineral to the sum of the peak area intensities of all detected minerals, and used this value to estimate percentages.

Diatoms

A total of 37 samples from sediment core PG2208 were analyses for diatoms, of which we found valves in 32 samples. Average resolutions for diatoms and other indicators over the Holocene and Pleistocene was generally centennial and millennial, respectively (listed in the Supplementary Table; Holocene 264 years, Pleistocene ca. 1,400 years). For species analysis we prepared diatom slides for light microscopy following the common procedure developed by Battarbee et al. (2001). About 0.1 g of freeze-dried sample material was treated with hydrogen peroxide (30%) for several hours and washed before adding 5 × 10⁶ microspheres to estimate diatom valve concentration (DVC). Diatom valves were identified using a Zeiss AXIO Scope.A1 light microscope with a Plan-Apochromat 100x/1.4 Oil Ph3 objective at 1,000× magnification. A minimum of 300 diatom valves (Wolfe, 1997) were counted (mean 316) in each sample. Only 23 valves were found in the deepest sample at 310.5 cm. In parallel to diatom analysis, Chrysophyte cysts and Mallomonas scales were counted but not further identified. Diatom species identification was supported by images from a scanning electron microscope (SEM) at GeoForschungZentrum Potsdam. The diatom species identification was carried to lowest possible taxonomic level. To guarantee correct identification literature was used including Hofmann et al. (2011) and Krammer and Lange-Bertalot (1986–1991) as well as online databases (i.e., http://www.algaebase.org).

Chironomids

A total of 57 samples were studied for chironomid analysis. Sample resolution was 317 years. in the Holocene and ca. 1800 years in the Pleistocene (Supplementary Table). However, the lower (291–366 cm) and the upper (0–40 cm) parts of the core contained less than 50 chironomid head capsule remains per sample, that is necessary for a reliable quantitative paleotemperature reconstruction (Heiri and Lotter, 2001; Quinlan and Smol, 2001). In order to reach the adequate number of head capsules for a reliable estimate of inferred temperature we merged adjacent samples from these two core intervals. Treatment of sediment samples for chironomid analysis followed standard techniques described in Brooks et al. (2007). Subsamples of wet sediments were deflocculated in 10% KOH, heated to 70°C for 10 min by adding boiling water and left for another 20 min. The sediment was passed through stacked 225 and 90 µm sieves. Chironomid larval head capsules were picked out of a grooved Bogorov sorting tray under a stereomicroscope and were mounted in Hydromatrix two at a time, ventral side up, under a 6 mm diameter cover slip. Chironomids were identified to the highest taxonomic resolution possible with reference to Wiederholm (1983) and Brooks et al. (2007). Information on the ecology of chironomid taxa was taken from Brooks et al. (2007), Nazarova et al. (2008), Nazarova et al. (2011a), Moller Pillot (2013), Nazarova et al. (2015) and Nazarova L. B. et al. (2017).

Data Processing and Statistical Techniques

For statistical analysis of diatoms and sediment-geochemical variables we used the R environment (R Core Team, 2016).

XRF data was cleaned, i.e., negative values (<1%) resulting from the model were excluded. Therefore, some graphs reveal minor data gaps. For better comparison with lower resolution proxies and to reduce noise, we applied 7-point running mean shown on top of unsmoothed data. We transformed the data to account for compositional data using centered log ratios (CLR) for single elements and additive log ratios (ALR) for ratios of two elements.

Diatom zones were established guided by a constrained incremental sums-of-squares clustering (CONISS) applied on all species counts while using Euclidean dissimilarity after log transformation to downweigh abundant species (Grimm, 1987). Before choosing a suitable ordination technique for diatom data we performed a detrended correspondence analysis (DCA) to calculate gradient length (in standard deviation SD units) using the “decorana” function in the package “vegan” (Oksanen et al., 2020). The DCA gradient length for the original count data was (SD = 3.1) which is above the threshold suggested by Birks (2010). Therefore, to reveal major downcore trends in diatom assemblage data, multidimensional scaling (MDS) technique based on Euclidean distance was chosen (Mardia, 1978) using the “cmdscale” function from the package “stats”. Before MDS analysis we filtered the diatom dataset to include only species present with ≥2.5% in ≥2 samples, excluded the bottom sample at 310 cm because of too scarce valve counts, and applied square-root transformation to downweigh abundant species (but less aggressively than log transformation on unfiltered data applied for CONISS).

We estimated diatom alpha diversity (species richness) based on Hill’s N0, and N2 diversity estimates (Hill, 1973) using all diatom counts excluding the bottom sample that had too low number of counts. We applied rarefaction to correct diversity estimates for sampling differences (Birks et al., 2016) using the minimum base sum of all samples (n = 228). Hill’s diversity number N0 estimates species richness, while N2 estimates the inverse of the Simpson concentration that amplifies abundant species but discounts rare species (Chao et al., 2014). Rarefaction and diversity indices were calculated using the “vegan” R package (Oksanen et al., 2020).

The quantitative chironomid-based reconstruction of mean July air temperatures (T_July) is based on calibration chironomid data sets of the lakes from northern Russia (Nazarova et al., 2015). T_July was inferred by using a North Russian (NR) chironomid-based temperature inference model (WA-PLS, 2 component; r2 boot = 0.81; RMSEP boot = 1.43°C) based on a modern calibration data set of 193 lakes and 162 taxa from East and West Siberia (61–75°N, 50–140°E, T_July range 1.8–18.8°C) (Nazarova et al., 2015). Modern mean July air temperature for the lakes from the calibration data set was derived from New et al. (2002). The T_July NR model was previously applied for paleoclimatic inferences in Europe, Eurasia with emphasis on Arctic Russia (Nazarova et al., 2013; Syrykh et al., 2017; Wetterich et al., 2018; Plikk et al., 2019; Nazarova L. et al., 2021; Nazarova L. B. et al., 2021), and demonstrated high reliability of the reconstructed parameters. The chironomid-inferred T_July were corrected to 0 m a.s.l. using a modern July air temperature lapse rate of 6°C km⁻¹ (Livingstone et al., 1999; Heiri et al., 2014). Chironomid-based reconstructions were performed and results were displayed in stratigraphic diagrams produced using the software C2 version 1.7.7 (Juggins, 2007). The chironomid data were square-root transformed to stabilize species variance. Zonation of stratigraphy was accomplished using the optimal sum-of-squares partitioning method (Birks and Gordon, 1985) and the program ZONE (Lotter and Juggins, 1991). The number of significant zones was assessed by a broken stick model (Bennett, 1996) using program BSTICK (Birks H.J.B. and Line J.M., unpublished).

We used the pollen dataset from Courtin et al. (2021) to perform July temperature reconstruction. Climate reconstructions are based on harmonized pollen data in accordance with Cao et al. (2013) i.e., woody taxa and common herb/shrub taxa are harmonized on genus level and all other taxa on family level; only terrestrial taxa were included. To minimize the impact of false-positive analogues (Cao et al., 2017), the used modern pollen training dataset was selected from sites within a 2000 km radius around the lake Bolshoe Toko from the Eurasian (Davis et al., 2020), a North American (Whitmore et al., 2005) and East Asian (Cao et al., 2014) modern pollen compilations. T_July and Pann values were extracted for 1701 modern pollen sites from WorldClim 2 data (Fick and Hijmans, 2017; https://www.worldclim.org/) having a T_July range from 3.1 to 24.3°C. The climate reconstructions were performed applying by using the MAT-function in the rioja package version 0.9-21 (Juggins, 2014) for R software (RCore Team, 2016). Cross-validation of modern data set yielded a root mean squared-error of prediction of 1.88°C and a r2 of 0.77 for observed vs. predicted T_July values. In addition, a statistical significance test (Telford and Birks, 2011) was performed for the reconstruction using the randomTF-function in the palaeoSig package (version 2.0-3).

To investigate lead-lag correlations we developed a moving temporal offset approach in the R environment (R Core Team, 2016). To make data sets comparable, new data points were predicted by applying local polynomial regression fitting using the “loess” function in the standard R package “stats”, based on Cleveland and Devlin (1988), and then linear interpolated the newly generated points to a 10-year resolution. This procedure helped us to minimize fluctuations and thus reduced noise in the integrated data set. It further ensured that data of different time series can be compared on a 10-year resolution for correlation.

We calculated the correlation coefficient for each 10-year interval step with a maximum offset of ±1,000 year while keeping one parameter fixed (either Hill’s N0 or N2 diatom diversity). Thus, 200 correlations were calculated for each involved parameter with the selected diatom alpha diversity, using the “cor” function of the package “compositions”, which is based on the Spearman correlation (Van Den Boogaart and Tolosana-Delgado, 2008). Parameters with a correlation coefficient above 0.7 or below −0.7 at any point are included.

Supportive and additional statistical analyses are provided in the Supplementary Material to this paper. To visualize multiproxy data structures, we performed a MDS and a Pearson correlation matrix based on the ensemble of indicator data (interpolated as described above for lead-lag correlation) presented in the study. These include chironomid- and pollen-based temperature reconstructions, sediment-geochemical variables (elements, grain-size, organic carbon estimates), and diatom-inferred variables (Supplementary Figures S1, S2). To compare our lead-lag analysis with other methods we also tested for cross-correlation. We used the same data pretreatment as described above but used 100-year interval steps instead. This change in intervals was necessary to make the data suitable for the correlation calculation in the “ccf” R function, based on Brockwell et al. (2016), for estimating the lead-lag relationships of two selected parameters. We used the outcome of our moving temporal offset approach to test the highest/lowest correlations of interest within the “ccf” function (Supplementary Figures S3–S5). 95% confidence intervals were used from the error margin of the standard error of the “loess” function.

We summarized the basic statistics of parameters used for this study in an overview, which is accessible in the Supplementary Material.

Results

Chronology

Age-depth modelling performed using the R package “Bacon” revealed a continuous age-depth relationship without any reversals (Figure 2). Inferred sedimentation rates from the applied age model show that accumulation during the Holocene (0.02 cm a⁻¹) was higher than in the Pleistocene (0.01 cm a⁻¹). Mean 95% confidence of age estimates over the Holocene were 842 years (±418 years), and 2,595 years (±1,298 years) over the Pleistocene core section. The maximum confidence range was 4,220 years (±2,110 years) at the bottom of the core (367 cm). The full set (100% of the dates) overlap with 95% envelope ranges developed by the Bacon model.

Sedimentology and Biogeochemistry

The Pleistocene part of the sediment core starting at 367 cm (ca. 23.3 cal ka BP), is characterized by medium grey, homogeneous, silty clay. In a lithological transition zone, it passes over to light brown finely laminated silty clay between 288 and 268 cm including the Pleistocene-Holocene boundary. The Holocene part, and especially the upper part starting at 200 cm (7.1 cal ka BP) to the top of the core, consists of homogenous but (partly visible) finely laminated brown organic mud (Figure 3).

FIGURE 3

FIGURE 3. Variability of grain size and organic geochemistry (TOC, total organic carbon; TC, total carbon; TN, total nitrogen; ∂13C, stable carbon isotopes) in Bolshoe Toko (sediment core PG2208). A stretched core photograph is added below the mean grain-size data at the left side of the figure showing the lithological development and fine lamination.

The generally silty grain-size composition revealed increasing but overall low mean grain size between the bottom and 268 cm (min. 2.4, max. 10.3, mean 4.9 µm). In this section sand is almost absent (mean 2.0 vol%) but clay revealed highest values (mean 31.2 vol%). In samples above 268 cm mean grain size revealed higher, but also strongly fluctuating values (min. 7.9, max. 38.7, mean 17.8 µm) and highest proportions of sand (mean 18.7 vol%), while clay decreased to 9.5 vol%. Mean proportion of silt over the entire core was 70.3 vol%.

The mean difference between TC and TOC is 0.3 wt%, with an exception of two samples at 0.25 cm (4.2 wt%) and 202.5 cm (3.7 wt%), which, regarding the high number of samples measured in both devices (n = 204), is assumed to represent outliers. In the Pleistocene, carbon is almost absent in the core. From the core bottom upward to ca. 288 cm mean carbon values remain below 1%. The transition zone between ca. 288 and 268 cm reveals rising carbon concentrations up to ca. 5 wt%. Carbon concentrations in the above Early Holocene between 268 and 200 cm are slightly but constantly lower (mean 5.5 wt%) than when compared to the mid- and Late Holocene part above 200 cm (mean 7.5 wt%). Fluctuations visible in (Figure 3) are likely an expression of the fine laminations. TOC/TN_atomic ratios generally follow the carbon concentrations increasing from ca. 6 to 16 by 268 cm and subsequently remaining continuously high between 16 and 20 until the core top.

The water content throughout the core increases slowly with lower values between the bottom and 268 cm (min. 30.1 wt%, mean 52.7 wt%) and higher values starting from 268 cm until the top (max. 90.5 wt%, mean 79.9 wt%).

δ¹³C values range from −27.4 to −26.7‰ (mean 27.0‰). While values decline in the Late Holocene, there is an increase to maximum values in the uppermost part of the core which corroborates findings from a nearby short core (PG2203, Figure 1), reported by Biskaborn et al. (2021).

Elements and Minerals

CLR and ALR transformed XRF data suggest a general shift in the sediment core at 280 cm (ca. 14.0 cal ka BP) (Figure 4). At that depth Sr/Rb ratios, Br, Mn/Fe, and Si/Al strongly increase, while Ti decrease and stay approximately uniform with minor but high-frequent fluctuations in the Holocene. The transition is rather abrupt with the exception of Br which starts to increase gradually already at 320 cm. The transition at 280 cm is in accordance with XRD inferred mineral composition which shows a general decline of plagioclase, K-feldspar, and hornblende at that depth, while quartz rises. The maximum (measured) quartz peak was found even earlier at 300 cm.

FIGURE 4

FIGURE 4. Sediment-geochemical downcore variability in Bolshoe Toko (sediment core PG2208). Elemental data from XRF scanning given as additive log ratios (ALR) or centered log ratios (CLR) with 5-point running averages (black lines). Mineral data from XRD given as semi-quantitative percentage data.

Diatom Species Assemblages

Diatoms occurred between 311 and 0 cm depth (ca. 18.0 and −0.063 cal ka BP) in the sediment core PG2208 (Figure 5). In total, 158 diatom species were identified, dominated by a variety of Cyclotella, Aulacoseira and small achnanthoid taxa. In addition, Pliocaenicus bolshetokoensis appeared in considerable quantities between 5.1 and 34.4% (mean 15.0%) in the entire core record. Diatom valve concentration varied between 1.1 and 9.2 10⁸ g⁻¹ (mean 3.4 10⁸ g⁻¹), started to increase at 200 cm (7.1 cal ka BP) towards the top of the core, revealed a low value at 101.5 cm and peaked at 35.25 cm. The valve preservation index (F-Index) was generally high, ranging between 0.79 and 0.94 in all samples except the lowermost two samples at 300 and 310 cm, that had 0.51 and 0.09, respectively. The lower part of the core between 290 and 311 cm was dominated by Cyclotella sp. C. intermedia and C. ocellata species reaching up to 47.8, 21.7, and 36.3% in some of the samples, respectively. Aulacoseira ambigua started to appear in high frequency between 290 and 250 cm (35.8–10.5%), remained at intermediate values (mean 9.0%) between 250 and 100 cm, and decreased until the top (mean 2.9%). A. subarctica had intermediate values (mean 7.1%) but peaked several times reaching maximum abundance of 23.5% at 250 cm, while A. distans and A. valida had rather minor occurrences with highest values between 290 and 240 cm. C. atomus and C. tripartita started to appear with minor occurrences in the lower part of the core but appeared in higher frequencies above ca. 200 cm, reaching 23.3 and 18.4%, respectively. C. comensis (max. 23.6%) started to appear in high frequencies at 200 cm, and similar to C. cyclopuncta, increased distinctly again between 70 cm and the top. Achnanthidium minutissimum was abundant throughout the record with medium values between 1.9 and 18.5% (mean 10.0%). Other noteworthy taxa occurred throughout the record including Tabellaria (0–10, mean 2.8%), Fragilaria gracilis (0–5.5, mean 2.0%), and Rossithidium pussilum (0–6.2, mean 1.4%), of which the latter decreased in the uppermost part of the core. Mallomonas scales appeared from 270 cm onward and peaked at 170.5, 110.5, 80.5, and 12.5 cm before the Mallomonas index (0–38.3, mean 12.8) reached a constant maximum in the uppermost 3.25 cm of the core. The chrysophyte index peaked between 270 and 210, 140 and 80 cm, and revealed minor peaks in the uppermost core section.

FIGURE 5

FIGURE 5. Relative abundance of diatom species, diatom valve concentration, Mallomonas and Crysophyte indices, and sums of selected planktonic Aulacoseira and Cyclotella (sensu lato) taxa. Zonation guided by CONISS. Values are plotted next to composite core depth and median age values from the Bacon model. The species shown represent taxa present with ≥2.5% in ≥2 samples including the bottom sample.

MDS analysis depicts grouping of key-planktonic and benthic species along the primary axis and in particular a distinct clustering Cyclotella versus Aulacoseira together with Pliocaenicus and Achnanthidium minutissimum along the secondary axis (Figure 6). Both axes together explain 76.1% of the data variance.

FIGURE 6

FIGURE 6. Biplot of the first two dimensions (Dim1, Dim2) generated by multi-dimensional scaling (MDS) analysis of diatom species data. Colored circles indicate eco-taxonomical clusters with similar environmental preferences. Explained variance of each dimension is given at the axis label.

4.5 Chironomid Species Assemblages

In total, 65 chironomid taxa were identified. In the lowermost samples (below 290 cm; 15.7 cal ka BP, Figure 7) we found low concentration of chironomid remains represented only by profundal taxa from the genera Heterotrissocladius (H. subpilosus-type, H. grimschawi-type, H marcidus-type, H. maeaeri-type 2). Diversity of chironomid communities remains low: Hill`s N2 varies between 1 and 3.7. Reconstructed T_July range from 7.1 to 7.8°C.

FIGURE 7

FIGURE 7. Relative proportions of the most abundant chironomid taxa (%), concentration of chironomid head capsules (HC), PCA axes 1 scores for chironomid data, chironomid-inferred mean July air temperature (T_July), Hill’s N2 diversity and relative abundances (%) of littoral, profundal, lotic and semiterrestrial taxa in the sediment core PG2208 of Lake Bolshoe Toko. Grey vertical line at T_July represents modern mean July air T (°C).

From the core depth of 280 cm (14.0 cal ka BP) upwards, the diversity and concentration of head capsules rise. Between 280 and 200 cm (14.0–7.1 cal ka BP) though still at low abundancies, several littoral taxa (Cladotanytarsus mancus-type, Corynoneura-type, Psectrocladius taxa, Orthocladius type S) and semiterrestrial Smittia-Parasmittia appear in the lake. Reconstructed T_July rise from 7.8°C at the core depth 291 cm (15.85 cal ka BP) to 10.0°C at the depth 211 cm (7.6 cal ka BP). Characteristic taxa for flowing waters (lotic) appear after the core depth of 231 cm (8.8 cal ka BP).

The strongest taxonomic shift in the chironomid assemblages, reflected as well by the decline of the PCA 1 below 0, was found at 201 cm depth (7.1 cal ka BP; Figure 7). At this depth we observe rise of species representing littoral, lotic and semiterrestrial taxa. Taxonomic diversity of chironomid communities (Hill’s N2) increases up to 14.6 at the core depth 121 cm (3.9 cal ka BP). Among the profundal taxa we observe a shift towards dominance of acidophobic Micropsectra insignilobus-type. Abundances of semiterrestrial and lotic taxa are the highest in the record and phytophilic taxa from Cricotopus and Orthocladius genera became more abundant. Reconstructed T_July between 201 (7.1 cal ka BP) and 111 cm (3.4 cal ka BP) of the core depth vary between 9 and 13°C.

At a depth of 161 cm (5.6 cal ka BP) we observe a sharp decrease of the head capsules concentration, decrease of the chironomid diversity and a strong increase in the abundance of the semiterrestrial taxon Limnophies-Paralimnophies.

From 151 cm (5.2 cal ka BP) to 111 cm (3.4 cal ka BP) the acidophobic Micropsectra insignilobus-type and lotic and semi-terrestrial fauna declined while profundal acidophilic Heterotrissocladius taxa increase again in line with the rise of Corynoneura arctica-type and several phytophilic taxa from the genera Cricotopus, Orthocladius, Paratanytarsus, Cladotanytarsus. The reconstructed T_July remain slightly above the modern temperature (10.5–11.9°C).

After 111 cm there is a break in the development of chironomid communities. No chironomids have been found till 86 cm, when the chironomids re-appear in the record. The proportion of littoral taxa rises up reaching 66% at 86 cm (2.4 cal ka BP) alongside with the latter increasing abundances of lotic fauna that reach 22% at the depth 71 cm (1.9 cal ka BP) which could indicate lake deepening and new expansion of littoral areas. At the 71 cm Hydrobaenus lugubris-type appears. Reconstructed T_July decreases from 10.5°C at 111 cm to 8.1°C at 71 cm.

Above 71 cm the concentrations of lotic taxa decrease and above 9 cm (ca. 188 cal years BP) they are not found in the chironomid communities anymore. The relative proportion of profundal to littoral taxa remain at ca. 2.2 on average until ca 9 cm. Between 9 and 7.5 cm (ca. 147 cal years BP) the concentration of littoral taxa starts to increase and the concentration of the profundal taxa decreases again. The communities are dominated by profundal acidophobic M. insignilobus-type. Reconstructed T_July has an increasing trend towards 15 cm (11.7°C; 394 cal years BP) then a decrease with a min. of 9.4°C at ca 8 cm and a rise again towards in modern time up to 11.2°C.

Discussion

Our following discussion will start with inferring the paleoclimate history from Lake Bolshoe Toko data and synthesize the findings with adjacent paleoenvironmental records, i.e., along a north-south profile through Yakutia. Subsequently we discuss the ecological developments of the lake and inter-proxy relationships, such as lead-lag correlations between temperature reconstruction, bioindicator indices, and sediment-geochemical trends.

Paleoclimate Variability in Northeastern Siberia

Evaluation of Climate Reconstructions From Bolshoe Toko Sediment Cores

Here we discuss pollen-based and chironomid-based T_July reconstructions performed on samples from two parallel sediment cores from Lake Bolshoe Toko. Chironomid-based reconstruction in core PG2208 provides information derived from aquatic bioindicators; pollen-based reconstruction taken from Courtin et al. (2021, core PG2133) provides terrestrial-based information from the catchment including aeolian pollen transport (Figure 1). Deviations between both reconstructions (Figure 8) can be related to different conceptual proxy and ecosystem settings. Distributions of chironomids are mainly driven by air temperature (Brooks et al., 2007; Self et al., 2011; Nazarova et al., 2015) but can be influenced by changes in limnic habitat preferences (e.g., water depth, substrate, ice cover, oxygen) (Stief et al., 2005; Nazarova et al., 2011b; Engels et al., 2012), while terrestrial vegetation compositions are limited by water supply (e.g., Antonsson et al., 2006). A major palynological constraint is that long-range transport of pollen obscures vegetation-based temperature reconstruction at the single lake site, that is especially true in mountain regions (Rosén et al., 2003) and amplifying uncertainties caused by complex vegetation-landscape interactions and long-term lagged climate response, i.e., in the Holocene (Cao et al., 2019). Therefore, in our study we use T_July (pollen) mainly for inferring climate trends at the Pleistocene/Holocene boundary where pollen shifts are well established (Velichko et al., 2002), and T_July (chironomid) within the Holocene interglacial, when limnoecological conditions were relatively stable. Potential deviations caused by age control or different sampling strategies cannot be ruled out, but are likely neglectable as common features for both reconstructions. One major feature for both data sets is a temperature rise at the Pleistocene-Holocene boundary. Here, pollen data reveal strong short-term fluctuations. Artemisia and Poaceae increased while Betula and Alnus decreased around 12.5 cal ka BP, suggesting evidence for a Younger Dryas cooling (Courtin et al., 2021). Climate cooling and increased continentality during the Younger dryas was reported for other sites in southern (Harding et al., 2020) and eastern Siberia (Velichko et al., 2002; Müller et al., 2009). However, tree-growth (i.e., Larix) persisted during this cold phase at multiple sites in eastern Siberia (Müller et al., 2009; Werner et al., 2010).

FIGURE 8

FIGURE 8. Graph showing synchronized variables used for lead-lag correlations. Temperature reconstructions (T_July, °C) by pollen are from core PG2133. All other variables are from core PG2208: chironomids, diatom alpha diversity estimates (Hill’s N2, N0), MDS first dimension Dim 1 (multiplied by −1), organic carbon concentration, selected elements (CLR values), and total sand concentrations. Variables are drawn using their ages from PG2208 and PG2133 for pollen (Courtin et al., 2021). Smoothing and confidence intervals were gained using LOESS. The less reliable (lowermost four) samples were excluded from the chironomid data due to the low number of head capsules. HTM, Holocene thermal maximum; P/H, Pleistocene/Holocene transition.

Our pollen-based reconstruction of summer temperatures from Bolshoe Toko suggests that T_July (pollen) fluctuated around 11°C during the late Glacial with slightly higher values during 13 to 12 cal ka BP and a reversal shortly before the early Holocene warming. Holocene temperature ranges around 13°C with a slight decline toward the late Holocene. The T_July reconstruction was significant (p < 0.028) (according to Telford and Birks, 2011). In detail, pollen reconstructions revealed a distinct but ambiguous climatic event between ca. 12.9 and 11.7 cal ka BP (Figure 8) indicating cold temperatures at the beginning and the end of the Younger Dryas period, but a warm interim period (Figure 8). Such a short warming event in the middle Younger Dryas in Siberian records was already considered by Velichko et al. (2002). The lower resolution of our T_July (chironomid) reconstruction in this period can somewhat validate Younger Dryas cooling but not a middle Younger Dryas short warming episode.

The chironomid taxa from the genera Heterotrissocladius that are dominant until 13.0 cal ka BP are profundal, oligotrophic, acidophilic, cold stenotherm (H. subpilosus-type, H. grimschawi-type) to cold tolerant (H. marcidus-type, H. maeaeri-type 2) (Brooks et al., 2007). Such composition of chironomid communities indicates cold and dystrophic lake conditions unfavorable for development of rich chironomid fauna. Head capsule concentration was too low to enable high resolution quantitative temperature reconstructions from samples deeper than 290 cm (15.8 cal ka BP). Nevertheless, the reconstructed T_July (chironomid) was ca 2–3°C below modern.

After ca. 14.0 cal ka BP rising concentration of chironomid remains, diversity of chironomid communities and emergence of littoral and semiterrestrial taxa indicate beginning of development of suitable littoral habitats at the lake shores and slowly rising temperatures (Figure 8). Emergence of the lotic taxa after 8.8 cal ka BP reflects strengthening of the surface water inflow into the lake that can be caused by warming and eventually thawing of cryotic (mountain glacier) catchment features leading to rise of the lake water level.

The strongest change in the chironomid assemblage occurred at 7.1 cal ka BP and is associated with a replacement of acidophilic profundal cold-stenotherm Heterotrissocladius taxa by acidophobic Micropsectra insignilobus-type and a number of littoral, semiterrestrial, phytophilic and lotic taxa. The dominating acidophobic Micropsectra insignilobus-type, though mostly attributed to profundal conditions, is known for thriving as well in the littoral zone of mountain lakes (Lindegaard, 1992; Bitušık and Kubovcık, 1999). We assume that water depth remained high with well-developed littoral zones overgrown by macrophytes. This major shift supports chironomid-reconstructed warm temperatures that reached about 2–3°C above todays values (Figure 8) indicating a mid-Holocene Thermal Maximum (HTM) starting at 7.1 cal ka BP.

However, sharp decrease of the head capsule concentration, decrease of chironomid diversity, and rising abundance of semiterrestrial Limnophies-Paralimnophies possibly indicate unstable lake levels (Massaferro and Brooks, 2002; Brooks et al., 2007) at around 5,600 cal years BP.

Between 5.2 and 2.0 cal ka BP Micropsectra insignilobus-type and lotic and semi-terrestrial fauna declined while cold stenotherm Heterotrissocladius taxa, Corynoneura arctica-type and several phytophilic taxa increased. Corynoneuraarctica-type is littoral, possibly preferring floating leaved macrophytes, especially duckweeds and can be acidophilic because species of this genus were found in lakes with a pH around 4 (Moller Pillot and Buskens, 1990). This taxonomic shift indicates Neoglacial cooling climate conditions at Bolshoe Toko and lowering of the water pH. The pH decrease can be caused by respiration and decomposition processes in shrinking and shallowing littoral zones (Nazarova L. et al., 2017; Nazarova et al., 2020).

Between ca. 3.4 and 2.4 cal ka BP, the temperatures dropped. Chironomids were not found in the sediments during this period, probably as a result of lake level and productivity lowering in response to Neoglacial cooling. This cooling trend was also observed in Kamchatka (Brooks et al., 2015) consistent with the strengthening of the Siberian High (Hammarlund et al., 2015). Furthermore it is prominent in East Eurasian records as it was amplified by thermohaline circulation (Wanner et al., 2011).

At ca. 2.3 cal ka BP the rise of littoral taxa and recovery of Micropsectra insignilobus-type with increasing abundancies of lotic fauna indicate increasing water flow peaking at ca. 1.9 cal ka BP under a relatively cool climatic condition (chironomid T_July drops up to 2°C below modern). Here, emergence of Hydrobaenus lugubris-type confirms strengthening of the water flow. This taxon is specially adopted to the floodplains of large rivers (Moller Pillot, 2013) and might have been brought to the lake with inflowing rivers.

During the last ca. 2000 years the relatively stable ratio of profundal to littoral taxa indicates stable water levels during this interval. Smaller shifts between habitat preferences after ca. 188 cal BP could potentially be associated to water level changes related to centennial/decadal climate disturbance during the little ice age. However, the low number of chironomid head capsules in core PG2208 and the mean Holocene resolution of pollen in core PG2133 of 750 years did not allow sufficient resolution of the reconstructions of ecological conditions during the last millennium to reach centennial scale.

Bolshoe Toko as the Southern End of a North-South Paleoclimate Transect in Siberia

Our study completes the southern end of a N-S transect of paleoclimate records in Yakutia (Figure 9) established by multiple studies along the Lena River over the last decades (Laing et al., 1999; Andreev et al., 2004; Fradkina et al., 2005; Popp, 2006; Müller et al., 2009; Biskaborn et al., 2012; Klemm et al., 2013; Nazarova et al., 2013; Tarasov et al., 2013; Biskaborn et al., 2016; Kostrova et al., 2021). The Holocene warming trend in Siberia is a common and very clear feature reflecting the glacial-interglacial transition in the Northern Hemisphere (Kaufman et al., 2020). However, previous investigations of lake records along the transect showed that the onset of most pronounced warming within the interglacial lagged by up to ∼3 kyrs in southern direction (Biskaborn et al., 2012), contradicting the assumption that warming was primarily related to solar insolation (Bond et al., 2001; Miller et al., 2010). The global compilation of Holocene temperature reconstructions by Marcott et al. (2013) suggested that the Early to mid-Holocene warming was modulated by deglaciation of Northern Hemisphere ice sheets. Biskaborn et al. (2012) found that the spatiotemporal pattern along a latitudinal gradient in eastern Siberia matches with global climate models (Renssen et al., 2009; Renssen et al., 2012) reconstructing deglaciation processes of the Laurentide ice-sheets and associated supply of cold air to Central Siberia until ∼7 cal ka BP. An exception to this latitudinal trend in eastern Siberia is provided by the Lake Emanda record, where diatom opal ∂18O and species assemblage changes suggested slightly earlier warming at 7.9 cal ka BP compared to sites along the Lena transect (Kostrova et al., 2021). An explanation might be the location of Lake Emanda to the east beyond the Verkhoyansk Mountains. As a consequence, the mountain range could have acted as a climate barrier for cold air supply by westerlies.

FIGURE 9

FIGURE 9. Spatiotemporal variability of the Holocene thermal Maximum (HTM) in Siberia compiled from paleolimnological studies (Laing et al., 1999; Andreev et al., 2004; Fradkina et al., 2005; Popp, 2006; Müller et al., 2009; Biskaborn et al., 2012; Klemm et al., 2013; Nazarova et al., 2013; Tarasov et al., 2013; Biskaborn et al., 2016; Kostrova et al., 2021). Horizontal dashed orange lines indicate the onset of the HTM in each record. This graph was adopted from Biskaborn et al. (2016) and further developed in eastern and southern directions with two new sites.

The new Lake Bolshoe Toko record presented here demonstrates that the cooling mechanism in the Early Holocene could possibly also be applicable for terrestrial aquatic settings close to the Sea of Okhotsk, as Bolshoe Toko is ca. 300 km far from its coast inland. Renssen’s model show that there is a tendency of delayed warming towards latitudes near Bolshoe Toko (Renssen et al., 2009). However, the relative proximity of our study site to the marine coast possibly exceeds the resolution of their model.

On the other hand, pollen and chironomid inferred summer temperatures differ (Figure 8) suggesting, apart from the mechanisms explained above, a bias caused by comparisons of summer versus winter sensitive proxies as highlighted by Baker et al. (2017). Their study demonstrates that winter climate dominated Holocene temperature dynamics is associated with decay of ice sheets in the Northern Hemisphere in continental Eurasia in contrast to regions more proximal to marine coasts.

Several short-term cold events disrupting the overall Holocene climate development were reported by Wanner et al. (2011), which, in most cases, are not found in the Bolshoe Toko record, possibly due to low temporal resolution of the study. A prominent cooling signal that only was found in the chironomid-based T_July reconstruction was dated to 5.6 cal ka BP (5.1–6.2 ka 2σ range) but could not be correlated to common cold events. There is, however, a cooling signal found at around 2.5 cal ka BP in both chironomid- and, to a lesser degree, in pollen-based reconstructions (Figure 8). If taken into account the confidence range of the age model at this core section (2.2–2.7 ka 2σ range) this cold phase could match to cold events between 3.3 and 2.5 ka BP that were found in other Eurasian records associated to Bond events (Wanner et al., 2008).

Paleolimnological Variability and Inter-proxy Relationships

Development of the Sedimentary Depositional Environment

The sediment-geochemical variability in core PG2208 reveals a major change at the Pleistocene/Holocene boundary from siliciclastic erosional regime, represented by high proportions of minerogenic clastic input and low organic matter contents, to higher lake and catchment primary productivity. This is indicated by rising organic matter concentrations and plant remains present in the Holocene (Table 1; Figures 3, 4). This change is accompanied by an increase in sedimentation rates, possibly slightly overestimated because of the low consolidation of the pore-water-rich sediments.

Abiotic proxy data exhibit a decline of detrital elements and minerals (Ti, Feldspar, Hornblende) after ca. 15 cal ka BP (Figure 4). Mn/Fe as oxygenation proxy (Naeher et al., 2013) and Si/Al as biogenic opal and thus diatom proxy (Vyse et al., 2020) increased distinctly and more abrupt at ca. 15 cal ka BP, while production of organic matter (TOC, TOC/TN_atomic) increased more gradually between ca. 15 and 11 cal ka BP together with coarsening grain size. Changes in the mineralogy of detrital components, basically reflected by an increase in quartz at the expense of hornblende, point to modified routes of exogenic sediment supply. Sediments enriched in hornblende indicate sources from the crystalline rocks of the Stanovoy Ranges, while quartz-rich detritus typically points to shore erosion of older Quaternary sediments around the northern fringe of the lake (Biskaborn B. K. et al., 2019). This sedimentary pattern implies a successive delay of bio-related responses such as bioproductivity after physico-chemical changes already occurred. We therefore assume that the sediment-geochemical depositional setting of the lake system first experienced a rather qualitative change driven by provenance changes in a shrinking deglacial mountain catchment. Successive changes, i.e. decreases in winter ice-cover already led to stronger oxygenation and diatom growth, before increased summer warming initiated a fully established Holocene lake ecosystem. As a consequence, the beginning of the interglacial period was accompanied by a decrease of terrigenous sediment supply to the benefit of biogenic sedimentation.

Paleoecological Interpretation of Diatom Taxa

To understand the downcore diatom assemblage trends in a paleolimnological context we attribute all cyclotelloid and Aulacoseira taxa found in the core to principle ecological groups summarized in Figure 5. Lindavia comensis is described as opportunistic and fast-growing (Rühland et al., 2015), oligotrophic reference species that occurs in August in the Great Lakes (Wunsam et al., 1995) and is positively correlated to temperature and thermocline depth (Reavie et al., 2017). Cyclotella cf. atomus occurs in spring in the Great Lakes positively correlated to thermocline depth (Reavie et al., 2017), but not clearly with temperature. Lindavia cyclopuncta is reported as an oligotrophic species (Scussolini et al., 2011). Lindavia intermedia is reported as oligotrophic “bodanicoid” (Novis et al., 2017) and thus not necessarily related to thermal stratification (Rühland et al., 2015), growing in late summer and competing well under low light under P-limitation (Malik, 2016). Lindavia ocellata is described as planktonic and periphytic, ultraoligotrophic to mesotrophic taxon (Wunsam et al., 1995). Lindavia tripartite is correlated with temperature in the Great Lakes (Reavie et al., 2017). Cyclotella cf. iris was reported from fossil records (Khursevich et al., 2001) and thus we believe that a sparce occurrence in the upper part of the core could be misinterpreted. Other Cyclotella-like valves in the record could not be identified to the species level, just a few counts were similar to Cyclotella schumannii, also found in Lake Immandra, Kola Penninsula (Denisov and Genkal, 2018). Given the above assessment of ecological preferences we defined C. comensis, cyclopuncta, tripartite and cf. atomus as planktonic and lightly silicified forms (Figure 5) that dwell under thermally stratified conditions.

The genus Aulacoseira generally builds strongly silicified, heavy, and thus rapidly sinking frustules that are abundant in deep boreal lakes (Laing and Smol, 2003). Euplanktonic forms A. ambigua and A. subarctica, however, differ in their lifestyles to tychoplanktonic taxa A. distans and A. valida. We therefore interpret A. ambigua and A. subarctica as pelagic forms which require water turbulence to remain suspended in the photic zone (Gibson et al., 2003; Rühland et al., 2008) for comparison with abundancies of light cyclotelloid forms preferring stratified water.

A special feature of Lake Bolshoe Toko is the dominance of newly described Pliocaenicus bolshetokoensis (Genkal et al., 2018; Genkal and Gabyshev, 2018), which is, as an endemic planktonic species, ecologically not yet fully understood. Biskaborn B. K. et al. (2019) found that its spatial within-lake variability might relate this taxon to increased nutrient availability at the Utuk river mouth and other small inflow streams. However, this rather is contradicted by the general appearance of Pliocaenicus costatus (sensu lato) forms which are found in cold, oligotrophic, circumneutral to weakly acidic mountain lakes in East Siberia (Cremer and Van De Vijver, 2006).

Benthic species reveal only minor long-term trends, i.e., the most frequent benthic/tychoplanktonic form Achnanthidium minutissimum reveals relative uniform downcore abundances. The main ecological signal is visible in the decrease of Rossithidium pussillum in the last 1.5 cal ka BP parallel to the maximum of DVC, which, because R. pussillum is an oligotrophy indicator can be an indicator that DVC is driven by nutrient supply.

The MDS analysis revealed key-planktonic and benthic species along the first axis and in particular a distinct clustering of light planktonic Cyclotella species versus heavy and pelagic Aulacoseira frustules together with planktonic Pliocaenicus and the tychoplanktonic Achnanthidium minutissimum along the secondary axis (Figure 6). The data structure thus supports that changes in these key species, with ecologies explained above, are i.e., driven by stratification. Sample scores of the first MDS dimension (Dim1, Figures 5, 6, 8) therefore may indicate changes from turbulent to stratified system tendencies.

Late Quaternary Diatom-Assemblage Trends Driven by Paleoclimate Change

The earliest preserved diatoms are represented by only a few (23) valves of Cyclotella and Pliocaenicus species dated to ca. 17.8 cal ka BP. While C. ocellata was reported from a high diversity of cold lakes with differing preferences from ultra-oligotrophic (Cremer and Wagner, 2003) to eutrophic (Schlegel and Scheffler, 1999), C. intermedia is more indicative for low light and nutrient conditions and thus both forms are able to withstand harsh environmental conditions (Malik, 2016). Fragilarioid forms are typical pioneering forms in other boreal shallow lakes (Valiranta et al., 2011) but are missing in this central and deep part of the lake basin.

Valve concentration increased and preservation improved at ca. 15.7 cal ka BP when Aulacoseira taxa and also Achnanthidium joined the assemblage. At the same time July temperatures reconstructed from pollen (and less distinctly from chironomids) increased (Figure 8). The Pleistocene-Holocene boundary is represented by strong increase of euplanktonic Aulacoseira taxa (i.e., A. subarctica) and thus an increase in turbulent lake conditions related to water mixing caused by strong summer winds. Fluctuating and rather cool reconstructed July temperatures support this assumption.

The first strong increase of euplanktonic and light Cyclotella taxa is dated to ca. 7.1 cal ka BP indicating a distinct increase in thermal summer stratification (Rühland et al., 2015; Biskaborn et al., 2021). This shift is accompanied by increases in DVC as a warming indicator (Biskaborn et al., 2012) and in accordance with the onset of the Holocene Thermal Maximum (HTM) indicated by chironomid-inferred July temperature. Increasing temperatures are however not ratified by pollen-based reconstructions (see “T_July (pollen)”, Figure 8).

At ca. 3.4 cal ka BP low values of DVC, TOC, Cyclotella, Crysophytes and Mallomonas and a small clay peak point to a short-term event with low lake bioproductivity related to cooling inferred by chironomid data. The second major increase in light planktonic Cyclotella taxa, together with high DVC values, indicates thermal summer stratification starting at ca. 2.0 cal ka BP. Until recent decades, summers with stratification tend to be frequent. Although this would be contradictory to decreasing insolation over the Late Holocene, this development of the aquatic proxies (partially including chironomid-inferred T_July) can be explained by long-term trends of warming Siberian winters (Meyer et al., 2015) possibly leading to prolonged open water seasons. In this context, increases in the proportion of planktonic taxa after ca. 2.0 cal ka BP (Figure 5) can also indicate water depth increase as a result of increased snow melt in the mountain catchment leading to strengthened river water input. This could partly also have contributed to the formation of the eastern lake terraces (Figure 1) discussed in Biskaborn B. K. et al. (2019) after retreat of the lake level.

For the most recent history, Biskaborn et al. (2021) showed, that there is a distinct increase of summer stratification in Lake Bolshoe Toko during the industrial era as response to recent northern hemispheric warming, accompanied with acidification trends represented by Mallomonas.

Relationships Between Diatom Diversity, Paleoclimate and Sediment-Geochemistry

Reconstructed temperatures and sediment-geochemical parameters were correlated with diatom diversity (Hill’s N0 and N2, Figures 10A,B), revealing strong positive correlations to Hill’s N2 (Figure 10A, Supplementary Figures S1, S2). Reconstructed temperatures from both pollen and chironomids fit well, i.e., chironomid-inferred temperatures show strongest correlation with neglectable time offset. Direct linkage between diatoms and climate due to their fast reproduction rates (Julius and Theriot, 2010), and a general lag of terrestrial vegetation (and thus pollen-inferred temperatures) have been reported by Liao et al. (2020). The good fit between T_July (chironomids) and diatom diversity is supported by near synchronous compositional changes of both chironomids and diatoms over the Holocene reported by Solovieva et al. (2015).

FIGURE 10

FIGURE 10. Lead-lag correlations of diatom alpha diversity are based on (A) Hill’s N2 and (B) Hill’s N0, and include all involved proxies correlating >0.7 (solid lines indicate positive correlation) or <−0.7 (dashed lines indicate negative correlation). Before analysis the data were LOESS smoothed and re-sampled into an integrated decadal resolution time series. Then, spearman-based correlation was applied in a successive bi-directional 10-year moving temporal offset to maxima of ±1,000 year. Cross-correlation tests based on auto correlation function are provided in the supplementary material.

However, our analyses generally showed that limiting factors for detection of reliable lead-lag relationships is sampling strategy and resolution of the proxies involved. Therefore, to account for deviations in the proxy sampling strategies (Supplementary Material), we only infer directional lead or lag behaviors of proxies and omit illusive absolute temporal durations in this study.

XRF-derived elemental variability revealed positive correlations of redox-related elements (Fe, Mn) and Cl related to water content, and negative correlations of detrital elements (Al, Ti, K), all showing that these sediment-geochemical proxies shifted before diatom alpha diversity responded. Positive correlations to sand and negative correlations to silt and clay confirm this trend. Interestingly, negative correlations ahead of species richness appears slightly stronger for Hill’s N0 index, suggesting that sedimentary changes, i.e., the reduction of very fine (clay, Al) and coarse (gravel) components somewhat allows the initiation of rare taxa growth. The strongest positive correlations are found between measurements of organic proxies (TOC, TC, N) and Hill’s N2, and, to a lesser degree Hill’s N0. The observed centennial-scale delay of organic carbon behind changes in diatom alpha diversity demonstrates a strong causal relationship between species richness and production of organic matter (Marzetz et al., 2017). This finding may motivate further research on diversity-carbon trajectories as understanding of nature-based solutions for carbon sequestration will increasingly important for future climate change management.

Conclusion

The study of radiocarbon dated sediment cores from the remote Siberian Lake Bolshoe Toko enabled new paleoclimate reconstructions that were used to extend a north-south transect in Yakutia. The multiproxy data were further analyzed to gain understanding of the relationships between different paleolimnological system components, i.e., lead-lag relationships between diatom alpha diversity versus sediment-geochemical and paleoenvironmental variability. We summarize our conclusions as follows:

• The PG2208 sediment core revealed 23.3 cal ka BP over a length of 367 cm with higher sedimentation rates in the Holocene (0.02 cm a⁻¹) than in the Pleistocene section (0.01 cm a⁻¹).

• High resolution of pollen-based summer temperature reconstruction T_July (pollen) provides evidence for Younger Dryas cooling between ca. 12.9 and 11.7 cal ka BP and possibly an interim warm period.

• Chironomid-inferred summer temperature reconstruction T_July (chironomid) reveals an onset of Holocene Thermal Maximum (HTM) at 7.1 cal ka BP.

• Extension of a Siberian north-south paleoclimate transect shows that the delay of the HTM reaches far south close to the Sea of Okhotsk.

• Relative shifts in abundances of euplanktonic diatom species Aulacoseira and Cyclotella reveal repeated episodes of thermal stratification in response to mid-Holocene paleoclimate changes.

• Lead-lag correlations between diatom alpha diversity, involved environmental proxy data and temperature reconstructions revealed a lead of geochemical parameters, quasi-synchronous paleoclimate responses, and a lag of sedimentary organic carbon. Centennial-scale delay of organic carbon behind changes in diatom Hill’s N0 and N2 reveals the strongest causal relationship, which implies that species richness can augment the production and accumulation of organic matter in comparable lake systems. This may contribute to assessments of future carbon sequestration during ongoing global warming.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

BB led the sediment-geochemical and diatom analyses and writing of the manuscript. LN and LS led the chironomid analyses including chironomid-based reconstructions. UH led the pollen-based reconstruction. BD and LP led the field work and contributed to the environmental interpretation of the data. TK and GP assisted with data processing and statistical analyses. All authors contributed to analysis of the results and revision of the manuscript.

Funding

This research was supported by the BMBF project PALMOD #01LP1510D (Germany).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The handling Editor declared a past co-authorship with one of the authors (LN).

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank Kim Funck for diatom identification and counting, and Vanessa Warmke for help in the laboratory for high resolution organic carbon and grain-size measurements.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2021.737353/full#supplementary-material

References

Amap (2017). Snow, Water, Ice and Permafrost in the Arctic (SWIPA) 2017. Oslo, Norway: Arctic Monitoring and Assessment Programme (AMAP).

Andreev, A. A., Raschke, E., Biskaborn, B. K., Vyse, S. A., Courtin, J., Böhmer, T., et al. (2021). Late Pleistocene to Holocene Vegetation and Climate Changes in Northwestern Chukotka (Far East Russia) Deduced from Lakes Ilirney and Rauchuagytgyn Pollen Records. Boreas 50, 652. doi:10.1111/bor.12521

CrossRef Full Text | Google Scholar

Andreev, A., Tarasov, P., Schwamborn, G., Ilyashuk, B., Ilyashuk, E., Bobrov, A., et al. (2004). Holocene Paleoenvironmental Records from Nikolay Lake, Lena River Delta, Arctic Russia. Palaeogeogr. Palaeoclimatol. Palaeoecol. 209, 197–217. doi:10.1016/j.palaeo.2004.02.010

CrossRef Full Text | Google Scholar

Antonsson, K., Brooks, S. J., Seppä, H., Telford, R. J., and Birks, H. J. B. (2006). Quantitative Palaeotemperature Records Inferred from Fossil Pollen and Chironomid Assemblages from Lake Gilltjärnen, Northern central Sweden. J. Quat. Sci. 21, 831–841. doi:10.1002/jqs.1004

CrossRef Full Text | Google Scholar

Appleby, P. G., Nolan, P. J., Gifford, D. W., Godfrey, M. J., Oldfield, F., Anderson, N. J., et al. (1986). 210Pb Dating by Low Background Gamma Counting. Hydrobiologia 143, 21–27. doi:10.1007/bf00026640

CrossRef Full Text | Google Scholar

Baker, J. L., Lachniet, M. S., Chervyatsova, O., Asmerom, Y., and Polyak, V. J. (2017). Holocene Warming in Western continental Eurasia Driven by Glacial Retreat and Greenhouse Forcing. Nat. Geosci. 10, 430–435. doi:10.1038/ngeo2953

CrossRef Full Text | Google Scholar

Barley, E. M., Walker, I. R., Kurek, J., Cwynar, L. C., Mathewes, R. W., Gajewski, K., et al. (2006). A Northwest North American Training Set: Distribution of Freshwater Midges in Relation to Air Temperature and lake Depth. J. Paleolimnol. 36, 295–314. doi:10.1007/s10933-006-0014-6

CrossRef Full Text | Google Scholar

Barnet, J. S. K., and Steiner, B. M. (2021). Unravelling the Complex Geological Evolution of One of Earth's Final Remaining Frontiers: East Siberia. Geology Today 37, 12–17. doi:10.1111/gto.12336

CrossRef Full Text | Google Scholar

Battarbee, R. W., Jones, V. J., Flower, R. J., Cameron, N. G., Bennion, H., Carvalho, L., et al. (2001). “Diatoms,” in Tracking Environmental Change Using Lake Sediments, Vol. 3. Terrestrial, Algal, and Siliceous Indicators. Editors J. P. Smol, H. J. B. Birks, and W. M. Last (Dordrecht, Netherlands: Kluwer Academic Publishers), 155–202.

Google Scholar

Battarbee, R. W. (2000). Palaeolimnological Approaches to Climate Change, with Special Regard to the Biological Record. Quat. Sci. Rev. 19, 107–124. doi:10.1016/s0277-3791(99)00057-8

CrossRef Full Text | Google Scholar

Bennett, K. D. (1996). Determination of the Number of Zones in a Biostratigraphical Sequence. New Phytol. 132, 155–170. doi:10.1111/j.1469-8137.1996.tb04521.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Birks, H., and Gordon, A. (1985). “The Analysis of Pollen Stratigraphical Data. Zonation,” in Numerical Methods in Quaternary Pollen Analysis (London: Academic Press), 47–90.

Google Scholar

Birks, H. J. B. (2010). “Numerical Methods for the Analysis of Diatom Assemblage Data,” in The Diatoms: Applications for the Environmental and Earth Science. Editors J. P. Smol, and E. F. Stoermer (Cambridge: Cambridge University Press), 23–54.

Google Scholar

Birks, H. J. B, Felde, V. A., Bjune, A. E., Grytnes, J.-A., Seppa, H., and Giesecke, T. (2016). Does Pollenassemblage Richness Reflect Floristic Richness? A Review of Recent Developments and Future Challenges. Rev. Palaeobot. Palynol. 228, 1–25.

Biskaborn, B. K., Herzschuh, U., Bolshiyanov, D., Savelieva, L., and Diekmann, B. (2012). Environmental Variability in Northeastern Siberia during the Last ∼ 13,300 Yr Inferred from lake Diatoms and Sediment-Geochemical Parameters. Palaeogeogr. Palaeoclimatol. Palaeoecol. 329-330, 22–36. doi:10.1016/j.palaeo.2012.02.003

CrossRef Full Text | Google Scholar

Biskaborn, B. K., Herzschuh, U., Bolshiyanov, D., Savelieva, L., Zibulski, R., and Diekmann, B. (2013a). Late Holocene Thermokarst Variability Inferred from Diatoms in a lake Sediment Record from the Lena Delta, Siberian Arctic. J. Paleolimnol. 49, 155–170. doi:10.1007/s10933-012-9650-1

CrossRef Full Text | Google Scholar

Biskaborn, B. K., Herzschuh, U., Bolshiyanov, D. Y., Schwamborn, G., and Diekmann, B. (2013b). Thermokarst Processes and Depositional Events in a Tundra Lake, Northeastern Siberia. Permafrost Periglac. Process. 24, 160–174. doi:10.1002/ppp.1769

CrossRef Full Text | Google Scholar

Biskaborn, B. K., Narancic, B., Stoof-Leichsenring, K. R., Pestryakova, L. A., Appleby, P. G., Piliposian, G. T., et al. (2021). Effects of Climate Change and Industrialization on Lake Bolshoe Toko, Eastern Siberia. J. Paleolimnology 65, 335. doi:10.1007/s10933-021-00175-z

CrossRef Full Text | Google Scholar

Biskaborn, B. K., Nazarova, L., Pestryakova, L. A., Syrykh, L., Funck, K., Meyer, H., et al. (2019a). Spatial Distribution of Environmental Indicators in Surface Sediments of Lake Bolshoe Toko, Yakutia, Russia. Biogeosciences 16, 4023–4049. doi:10.5194/bg-16-4023-2019

CrossRef Full Text | Google Scholar

Biskaborn, B. K., Smith, S. L., Noetzli, J., Matthes, H., Vieira, G., Streletskiy, D. A., et al. (2019b). Permafrost Is Warming at a Global Scale. Nat. Commun. 10, 264. doi:10.1038/s41467-018-08240-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Biskaborn, B. K., Subetto, D. A., Savelieva, L. A., Vakhrameeva, P. S., Hansche, A., Herzschuh, U., et al. (2016). Late Quaternary Vegetation and lake System Dynamics in north-eastern Siberia: Implications for Seasonal Climate Variability. Quat. Sci. Rev. 147, 406–421. doi:10.1016/j.quascirev.2015.08.014

CrossRef Full Text | Google Scholar

Bitušık, P., and Kubovcık, V. (1999). Sub-fossil Chironomids (Diptera: Chironomidae) from the Sediments of the Nizné Terianske Pleso (High Tatra Mts., Slovakia). Dipterologica bohemoslovaca 9, 11–20.

Google Scholar

Blaauw, M., and Christen, J. A. (2011). Flexible Paleoclimate Age-Depth Models Using an Autoregressive Gamma Process. Bayesian Anal. 6, 457–474. doi:10.1214/ba/1339616472

CrossRef Full Text | Google Scholar

Bond, G., Kromer, B., Beer, J., Muscheler, R., Evans, M., Showers, W., et al. (2001). Persistent Solar Influence on north Atlantic Climate during the Holocene. Science 294, 2130–2136. doi:10.1126/science.1065680

PubMed Abstract | CrossRef Full Text | Google Scholar

Brockwell, P. J., Brockwell, P. J., Davis, R. A., and Davis, R. A. (2016). Introduction to Time Series and Forecasting. Berlin, Germany: Springer.

Brooks, S. J., Diekmann, B., Jones, V. J., and Hammarlund, D. (2015). Holocene Environmental Change in Kamchatka: A Synopsis. Glob. Planet. Change 134, 166–174. doi:10.1016/j.gloplacha.2015.09.004

CrossRef Full Text | Google Scholar

Brooks, S. J., Langdon, P. G., and Heiri, O. (2007). The Identification and Use of Palaearctic Chironomidae Larvae in Palaeoecology. London: Quaternary Research Association.

Cao, X.-Y., Herzschuh, U., Telford, R. J., and Ni, J. (2014). A Modern Pollen-Climate Dataset from China and Mongolia: Assessing its Potential for Climate Reconstruction. Rev. Palaeobotany Palynology 211, 87–96. doi:10.1016/j.revpalbo.2014.08.007

CrossRef Full Text | Google Scholar

Cao, X.-Y., Ni, J., Herzschuh, U., Wang, Y.-B., and Zhao, Y. (2013). A Late Quaternary Pollen Dataset from Eastern continental Asia for Vegetation and Climate Reconstructions: Set up and Evaluation. Rev. Palaeobotany Palynology 194, 21–37. doi:10.1016/j.revpalbo.2013.02.003

CrossRef Full Text | Google Scholar

Cao, X., Tian, F., Li, F., Gaillard, M. J., Rudaya, N., Xu, Q., et al. (2019). Pollen-based Quantitative Land-Cover Reconstruction for Northern Asia Covering the Last 40 ka cal BP. Clim. Past 15, 1503–1536. doi:10.5194/cp-15-1503-2019

CrossRef Full Text | Google Scholar

Cao, X., Tian, F., Telford, R. J., Ni, J., Xu, Q., Chen, F., et al. (2017). Impacts of the Spatial Extent of Pollen-Climate Calibration-Set on the Absolute Values, Range and Trends of Reconstructed Holocene Precipitation. Quat. Sci. Rev. 178, 37–53. doi:10.1016/j.quascirev.2017.10.030

CrossRef Full Text | Google Scholar

Chao, A., Gotelli, N. J., Hsieh, T. C., Sander, E. L., Ma, K. H., Colwell, R. K., et al. (2014). Rarefaction and Extrapolation with Hill Numbers: a Framework for Sampling and Estimation in Species Diversity Studies. Ecol. Monogr. 84, 45–67. doi:10.1890/13-0133.1

CrossRef Full Text | Google Scholar

Cleveland, W. S., and Devlin, S. J. (1988). Locally Weighted Regression: an Approach to Regression Analysis by Local Fitting. J. Am. Stat. Assoc. 83, 596–610. doi:10.1080/01621459.1988.10478639

CrossRef Full Text | Google Scholar

Cohen, A. S. (2003). Palaeolimnology - the History and Evolution of Lake Systems. Oxford: Oxford University Press. doi:10.1093/oso/9780195133530.001.0001

CrossRef Full Text

Courtin, J., Andreev, A. A., Raschke, E., Bala, S., Biskaborn, B. K., Liu, S., et al. (2021). Vegetation Changes in Southeastern Siberia during the Late Pleistocene and the Holocene. Front. Ecol. Evol. 9, 625096. doi:10.3389/fevo.2021.625096

CrossRef Full Text | Google Scholar

Cremer, H., and Van De Vijver, B. (2006). OnPliocaenicus costatus(Bacillariophyceae) in Lake El'gygytgyn, East Siberia. Eur. J. Phycology 41, 169–178. doi:10.1080/09670260600621932

CrossRef Full Text | Google Scholar

Cremer, H., and Wagner, B. (2003). The Diatom flora in the Ultra-oligotrophic Lake El'gygytgyn, Chukotka. Polar Biol. 26, 105–114. doi:10.1007/s00300-002-0445-0

CrossRef Full Text | Google Scholar

Davis, B. A. S., Chevalier, M., Sommer, P., Carter, V. A., Finsinger, W., Mauri, A., et al. (2020). The Eurasian Modern Pollen Database (EMPD), Version 2. Earth Syst. Sci. Data 12, 2423–2445. doi:10.5194/essd-12-2423-2020

CrossRef Full Text | Google Scholar

Denisov, D., and Genkal, S. (2018). Centric Diatom Algae of Lake Imandra (Kola Peninsula, Russia). Int. J. Algae 20. doi:10.1615/interjalgae.v20.i1.30

CrossRef Full Text | Google Scholar

Engels, S., Cwynar, L. C., Rees, A. B. H., and Shuman, B. N. (2012). Chironomid-based Water Depth Reconstructions: an Independent Evaluation of Site-specific and Local Inference Models. J. Paleolimnol. 48, 693–709. doi:10.1007/s10933-012-9638-x

CrossRef Full Text | Google Scholar

Engels, S., Medeiros, A. S., Axford, Y., Brooks, S. J., Heiri, O., Luoto, T. P., et al. (2020). Temperature Change as a Driver of Spatial Patterns and Long‐term Trends in Chironomid (Insecta: Diptera) Diversity. Glob. Change Biol. 26, 1155–1169. doi:10.1111/gcb.14862

PubMed Abstract | CrossRef Full Text | Google Scholar

Fick, S. E., and Hijmans, R. J. (2017). WorldClim 2: New 1‐km Spatial Resolution Climate Surfaces for Global Land Areas. Int. J. Climatol. 37, 4302–4315. doi:10.1002/joc.5086

CrossRef Full Text | Google Scholar

Fradkina, A. F., Alekseev, M. N., Andreev, A. A., and Klimanov, V. A. (2005). “East Siberia (Based on Data Obtained Mainly in Central Yakutia),” in Cenozoic Climatic and Environmental Changes in Russia. Editors A. A. Velichko, and V. P. Nechaev (Boulder, Colorado: Geological Society of America), 89–103.

CrossRef Full Text | Google Scholar

Genkal, S., Gabyshev, V., Kulilovskiy, M., and Kuznetsova, I. (2018). Pliocaenicus Bolshetokoensis – a New Species from Lake Bolshoe Toko (Yakutia, Eastern Siberia, Russia). Diatom Res. 33 (2), 1–9. doi:10.1080/0269249x.2018.1477690

CrossRef Full Text | Google Scholar

Genkal, S., and Gabyshev, V. (2018). New Records of Centric Diatoms from Yakutia (Bolshoe Toko Lake): SEM Morphology, Ecology and Distribution. Novosti sistematiki nizshikh rastenii 52, 245. doi:10.31111/nsnr/2018.52.2.245

CrossRef Full Text | Google Scholar

Gibson, C. E., Anderson, N. J., and Haworth, E. Y. (2003). Aulacoseira Subarctica: Taxonomy, Physiology, Ecology and Palaeoecology. Eur. J. Phycology 38, 83–101. doi:10.1080/0967026031000094102

CrossRef Full Text | Google Scholar

Golovin, G. S., Bychev, M. I., Moskalenko, T. V., Petrova, G. I., and Mikheev, V. A. (2007). Coal Reserves of Republic of Sakha (Yakutia) and Principal Lines of Their Use. Soil Fuel Chem. 41, 65–70. doi:10.3103/s0361521907020012

CrossRef Full Text | Google Scholar

Grimm, E. C. (1987). CONISS: a FORTRAN 77 Program for Stratigraphically Constrained Cluster Analysis by the Method of Incremental Sum of Squares. Comput. Geosciences 13, 13–35. doi:10.1016/0098-3004(87)90022-7

CrossRef Full Text | Google Scholar

Hammarlund, D., Klimaschewski, A., St. Amour, N. A. N. A., Andrén, E., Self, A. E., Solovieva, N., et al. (2015). Late Holocene Expansion of Siberian dwarf pine (Pinus Pumila) in Kamchatka in Response to Increased Snow Cover as Inferred from Lacustrine Oxygen-Isotope Records. Glob. Planet. Change 134, 91–100. doi:10.1016/j.gloplacha.2015.04.004

CrossRef Full Text | Google Scholar

Harding, P., Bezrukova, E. V., Kostrova, S. S., Lacey, J. H., Leng, M. J., Meyer, H., et al. (2020). Hydrological (In)stability in Southern Siberia during the Younger Dryas and Early Holocene. Glob. Planet. Change 195, 103333. doi:10.1016/j.gloplacha.2020.103333

CrossRef Full Text | Google Scholar

Heinecke, L., Mischke, S., Adler, K., Barth, A., Biskaborn, B. K., Plessen, B., et al. (2017). Climatic and Limnological Changes at Lake Karakul (Tajikistan) during the Last ∼29 Cal Ka. J. Paleolimnol. 58, 317–334. doi:10.1007/s10933-017-9980-0

CrossRef Full Text | Google Scholar

Heiri, O., Brooks, S. J., Renssen, H., Bedford, A., Hazekamp, M., Ilyashuk, B., et al. (2014). Validation of Climate Model-Inferred Regional Temperature Change for Late-Glacial Europe. Nat. Commun. 5, 4914. doi:10.1038/ncomms5914

PubMed Abstract | CrossRef Full Text | Google Scholar

Heiri, O., and Lotter, A. F. (2001). Effect of Low Count Sums on Quantitative Environmental Reconstructions: an Example Using Subfossil Chironomids. J. Paleolimnology 26, 343–350. doi:10.1023/a:1017568913302

CrossRef Full Text | Google Scholar

Herzschuh, U., Pestryakova, L. A., Savelieva, L. A., Heinecke, L., Böhmer, T., Biskaborn, B., et al. (2013). Siberian Larch Forests and the Ion Content of Thaw Lakes Form a Geochemically Functional Entity. Nat. Commun. 4, 2408. doi:10.1038/ncomms3408

PubMed Abstract | CrossRef Full Text | Google Scholar

Hill, M. O. (1973). Diversity and evenness ‐ unifying notation and its consequences. Ecology 54, 427–432.

Hobbs, W. O., Telford, R. J., Birks, H. J. B., Saros, J. E., Hazewinkel, R. R. O., Perren, B. B., et al. (2010). Quantifying Recent Ecological Changes in Remote Lakes of North America and Greenland Using Sediment Diatom Assemblages. PLOS ONE 5, e10026. doi:10.1371/journal.pone.0010026

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoff, U., Biskaborn, B. K., Dirksen, V. G., Dirksen, O., Kuhn, G., Meyer, H., et al. (2015). Holocene Environment of Central Kamchatka, Russia: Implications from a Multi-Proxy Record of Two-Yurts Lake. Glob. Planet. Change 134, 101–117. doi:10.1016/j.gloplacha.2015.07.011

CrossRef Full Text | Google Scholar

Hofmann, G., Lange-Bertalot, H., and Werum, M. (2011). Diatomeen im Süßwasser - Benthos von Mitteleuropa. Waldkirch, Germany: Ganter Verlag.

Imaeva, L. P., Imaev, V. S., and Koz’min, B. M. (2012). Seismogeodynamics of the Aldan-Stanovoi Block. Russ. J. Pac. Geol. 6, 1–12. doi:10.1134/s1819714012010071

CrossRef Full Text | Google Scholar

Juggins, S. (2014). “rioja: Analysis of Quaternary Science Data, R Package”. 0.9.3 Ed.

Juggins, S. (2007). “C2 Version 1.5 User Guide,” in Software for Ecological and Palaeoecological Data Analysis and Visualisation (Newcastle upon Tyne, UK: Newcastle University), 73.

Google Scholar

Julius, M. L., and Theriot, E. C. (2010). “The Diatoms: a Primer,” in The Diatoms: Applications for the Environmental and Earth Sciences. Editors J. P. Smol, and E. F. Stoermer (Cambridge: Cambridge University Press), 8e22.

Google Scholar

Kaufman, D., Ager, T. A., Anderson, N. J., Anderson, P. M., Andrews, J. T., Bartlein, P. J., et al. (2004). Holocene thermal Maximum in the Western Arctic (0-180°W). Quat. Sci. Rev. 23, 529–560. doi:10.1016/j.quascirev.2003.09.007

CrossRef Full Text | Google Scholar

Kaufman, D., Mckay, N., Routson, C., Erb, M., Davis, B., Heiri, O., et al. (2020). A Global Database of Holocene Paleotemperature Records. Scientific Data 7, 115. doi:10.1038/s41597-020-0515-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Khursevich, G. K., Karabanov, E. B., Prokopenko, A. A., Williams, D. F., Kuzmin, M. I., and Fedenya, S. A. (2001). Biostratigraphic Significance of New Fossil Species of the Diatom Genera Stephanodiscus and Cyclotella from Upper Cenozoic Deposits of Lake Baikal, Siberia. Micropaleontology 47, 47–71. doi:10.2113/47.1.47

CrossRef Full Text | Google Scholar

Klemm, J., Herzschuh, U., Pisaric, M. F. J., Telford, R. J., Heim, B., and Pestryakova, L. A. (2013). A Pollen-Climate Transfer Function from the Tundra and Taiga Vegetation in Arctic Siberia and its Applicability to a Holocene Record. Palaeogeogr. Palaeoclimatol. Palaeoecol. 386, 702–713. doi:10.1016/j.palaeo.2013.06.033

CrossRef Full Text | Google Scholar

Konstantinov, A. F. (2000). “Environmental Problems of lake Bolshoe Toko (In Russian),” in Lakes of Cold Environments, Part V: Resource Study, Resource Use, Ecology and Nature Protection Issue. Editors I. I. Zhirkov, L. A. Pestryakova, and G. N. Maksimov (Yakutsk, Russia: Yakutsk State University Publishing House), 85–93.

Google Scholar

Kornilov, B. A. (1962). Relief: The Southeast Suburbs of Aldan Mountains (In Russian). Moscow: Publishing House of Academy of Sciences of the USSR.

Kostrova, S. S., Biskaborn, B. K., Pestryakova, L. A., Fernandoy, F., Lenz, M. M., and Meyer, H. (2021). Climate and Environmental Changes of the Lateglacial Transition and Holocene in Northeastern Siberia: Evidence from Diatom Oxygen Isotopes and Assemblage Composition at Lake Emanda. Quat. Sci. Rev. 259, 106905. doi:10.1016/j.quascirev.2021.106905

CrossRef Full Text | Google Scholar

Krammer, K., and Lange-Bertalot, H. (1986-1991). Bacillariophyceae Band 2/2. Stuttgart: Gustav Fischer Verlag.

Kuznetsova, L. V., Zakharova, V. I., Sosina, N. K., Nikolin, E. G., Ivanova, E. I., Sofronova, E. V., et al. (2010). “Flora of Yakutia: Composition and Ecological Structure,” in The Far North (Berlin, Germany: Springer), 24–140. doi:10.1007/978-90-481-3774-9_2

CrossRef Full Text | Google Scholar

Laing, T. E., Rühland, K. M., and Smol, J. P. (1999). Past Environmental and Climatic Changes Related to Tree-Line Shifts Inferred from Fossil Diatoms from a lake Near the Lena River Delta, Siberia. The Holocene 9, 547–557. doi:10.1191/095968399675614733

CrossRef Full Text | Google Scholar

Laing, T. E., and Smol, J. P. (2003). Late Holocene Environmental Changes Inferred from Diatoms in a lake on the Western Taimyr Peninsula, Northern Russia. J. Paleolimnology 30, 231–247. doi:10.1023/a:1025561905506

CrossRef Full Text | Google Scholar

Liao, M., Herzschuh, U., Wang, Y., Liu, X., Ni, J., and Li, K. (2020). Lake Diatom Response to Climate Change and Sedimentary Events on the southeastern Tibetan Plateau during the Last Millennium. Quat. Sci. Rev. 241, 106409. doi:10.1016/j.quascirev.2020.106409

CrossRef Full Text | Google Scholar

Lindegaard, C. (1992). Zoobenthos Ecology of Thingvallavatn: Vertical Distribution, Abundance, Population Dynamics and Production. Oikos 64, 257–304. doi:10.2307/3545055

CrossRef Full Text | Google Scholar

Livingstone, D. M., Lotter, A. F., and Walkery, I. R. (1999). The Decrease in Summer Surface Water Temperature with Altitude in Swiss Alpine Lakes: a Comparison with Air Temperature Lapse Rates. Arctic, Antarctic, Alpine Res. 31, 341–352. doi:10.1080/15230430.1999.12003319

CrossRef Full Text | Google Scholar

Lotter, A., and Juggins, S. (1991). “POLPROF, TRAN and ZONE: Programs for Plotting, Editing and Zoning Pollen and Diatom Data,” in INQUA-Subcommission for the study of the Holocene Working Group on Data-Handling Methods, Newsletter 6, 4–6.

Google Scholar

Luoto, T. P., Rantala, M. V., Galkin, A., Rautio, M., and Nevalainen, L. (2016). Environmental Determinants of Chironomid Communities in Remote Northern Lakes across the Treeline - Implications for Climate Change Assessments. Ecol. Indicators 61, 991–999. doi:10.1016/j.ecolind.2015.10.057

CrossRef Full Text | Google Scholar

Luoto, T. P., and Salonen, V.-P. (2010). Fossil Midge Larvae (Diptera: Chironomidae) as Quantitative Indicators of Late-winter Hypolimnetic Oxygen in Southern Finland: a Calibration Model, Case Studies and Potentialities. Boreal Environ. Res. 15 (1), 1–18.

Google Scholar

Malik, H. I. (2016). Deciphering the Mechanisms behind Climate-Driven Changes in the Relative Abundances of the Diatom Cyclotella. Electron. Theses Dissertations 2545, 7845085. doi:10.1109/intech.2016.7845085

CrossRef Full Text | Google Scholar

Marcott, S. A., Shakun, J. D., Clark, P. U., and Mix, A. C. (2013). A Reconstruction of Regional and Global Temperature for the Past 11,300 Years. Science 339, 1198–1201. doi:10.1126/science.1228026

PubMed Abstract | CrossRef Full Text | Google Scholar

Mardia, K. V. (1978). Some Properties of Clasical Multi-Dimesional Scaling. Commun. Stat. - Theor. Methods 7, 1233–1241. doi:10.1080/03610927808827707

CrossRef Full Text | Google Scholar

Marzetz, V., Koussoroplis, A.-M., Martin-Creuzburg, D., Striebel, M., and Wacker, A. (2017). Linking Primary Producer Diversity and Food Quality Effects on Herbivores: A Biochemical Perspective. Scientific Rep. 7, 11035. doi:10.1038/s41598-017-11183-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Massaferro, J., and Brooks, S. J. (2002). Response of Chironomids to Late Quaternary Environmental Change in the Taitao Peninsula, Southern Chile. J. Quater. Sci. 17, 101–111.

Mckay, N. P., Kaufman, D. S., Routson, C. C., Erb, M. P., and Zander, P. D. (2018). The Onset and Rate of Holocene Neoglacial Cooling in the Arctic. Geophys. Res. Lett. 45, 12487. doi:10.1029/2018GL079773

CrossRef Full Text | Google Scholar

Meyer, H., Opel, T., Laepple, T., Dereviagin, A. Y., Hoffmann, K., and Werner, M. (2015). Long-term winter Warming Trend in the Siberian Arctic during the Mid- to Late Holocene. Nat. Geosci. 8, 122–125. doi:10.1038/ngeo2349

CrossRef Full Text | Google Scholar

Miller, G. H., Alley, R. B., Brigham-Grette, J., Fitzpatrick, J. J., Polyak, L., Serreze, M. C., et al. (2010). Arctic Amplification: Can the Past Constrain the Future? Quat. Sci. Rev. 29, 1779–1790. doi:10.1016/j.quascirev.2010.02.008

CrossRef Full Text | Google Scholar

Mock, C. J., Bartlein, P. J., and Anderson, P. M. (1998). Atmospheric Circulation Patterns and Spatial Climatic Variations in Beringia. Int. J. Climatol. 18, 1085–1104. doi:10.1002/(sici)1097-0088(199808)18:10<1085::aid-joc305>3.0.co;2-k

CrossRef Full Text | Google Scholar

Moller Pillot, H., and Buskens, R. (1990). De larven der Nederlandse Chironomidae (Diptera): autoekologie en verspreiding. Nederland: 1C. Nederlandse faunistische mededelingen.

Moller Pillot, H. (2013). Chironomidae Larvae of the Netherlands and Adjacent Lowlands, Biology and Ecology of the Aquatic Orthocladiinae. Zeist: KNNV Publishing. doi:10.1163/9789004278059

CrossRef Full Text

Müller, S., Tarasov, P. E., Andreev, A. A., and Diekmann, B. (2009). Late Glacial to Holocene Environments in the Present-Day Coldest Region of the Northern Hemisphere Inferred from a Pollen Record of Lake Billyakh, Verkhoyansk Mts, NE Siberia. Clim. Past 5, 73–84. doi:10.5194/cp-5-73-2009

CrossRef Full Text | Google Scholar

Naeher, S., Gilli, A., North, R. P., Hamann, Y., and Schubert, C. J. (2013). Tracing Bottom Water Oxygenation with Sedimentary Mn/Fe Ratios in Lake Zurich, Switzerland. Chem. Geology. 352, 125–133. doi:10.1016/j.chemgeo.2013.06.006

CrossRef Full Text | Google Scholar

Nazarova, L. B., Frolova, L. A., Palagushkina, O. V., Rudaya, N. A., Syrykh, L. S., Grekov, I. M., et al. (2021b). Recent Shift in Biological Communities: A Case Study from the Eastern European Russian Arctic (Bol`shezemelskaya Tundra). Polar Biol. 44, 1107–1125. doi:10.1007/s00300-021-02876-7

CrossRef Full Text | Google Scholar

Nazarova, L., Bleibtreu, A., Hoff, U., Dirksen, V., and Diekmann, B. (2017a). Changes in Temperature and Water Depth of a Small mountain lake during the Past 3000 Years in Central Kamchatka Reflected by a Chironomid Record. Quat. Int. 447, 46–58. doi:10.1016/j.quaint.2016.10.008

CrossRef Full Text | Google Scholar

Nazarova, L. B., Pestryakova, L. A., Ushnitskaya, L. A., and Hubberten, H.-W. (2008). Chironomids (Diptera: Chironomidae) in Lakes of central Yakutia and Their Indicative Potential for Paleoclimatic Research. Contemp. Probl. Ecol. 1, 335–345. doi:10.1134/s1995425508030089

CrossRef Full Text | Google Scholar

Nazarova, L. B., Self, A. E., Brooks, S. J., Solovieva, N., Syrykh, L. S., and Dauvalter, V. A. (2017b). Chironomid Fauna of the Lakes from the Pechora River basin (East of European Part of Russian Arctic): Ecology and Reconstruction of Recent Ecological Changes in the Region. Contemp. Probl. Ecol. 10, 350–362. doi:10.1134/s1995425517040059

CrossRef Full Text | Google Scholar

Nazarova, L., Herzschuh, U., Wetterich, S., Kumke, T., and Pestryakova, L. (2011a). Chironomid-based Inference Models for Estimating Mean July Air Temperature and Water Depth from Lakes in Yakutia, Northeastern Russia. J. Paleolimnol. 45, 57–71. doi:10.1007/s10933-010-9479-4

CrossRef Full Text | Google Scholar

Nazarova, L., Lüpfert, H., Subetto, D., Pestryakova, L., and Diekmann, B. (2013). Holocene Climate Conditions in central Yakutia (Eastern Siberia) Inferred from Sediment Composition and Fossil Chironomids of Lake Temje. Quat. Int. 290-291, 264–274. doi:10.1016/j.quaint.2012.11.006

CrossRef Full Text | Google Scholar

Nazarova, L., Razjigaeva, N. G., Ganzey, L. A., Makarova, T. R., Lyashevskaya, M. S., Biskaborn, B. K., et al. (2021a). The Middle to Late Holocene Environment on the Iturup Island (Kurils, North Western Pacific). Quat. Int. doi:10.1016/j.quaint.2021.05.003

CrossRef Full Text | Google Scholar

Nazarova, L., Self, A. E., Brooks, S. J., Van Hardenbroek, M., Herzschuh, U., and Diekmann, B. (2015). Northern Russian Chironomid-Based Modern Summer Temperature Data Set and Inference Models. Glob. Planet. Change 134, 10–25. doi:10.1016/j.gloplacha.2014.11.015

CrossRef Full Text | Google Scholar

Nazarova, L., Syrykh, L. S., Mayfield, R. J., Frolova, L. A., Ibragimova, A. G., Grekov, I. M., et al. (2020). Palaeoecological and Palaeoclimatic Conditions on the Karelian Isthmus (Northwestern Russia) during the Holocene. Quat. Res. 95, 65–83. doi:10.1017/qua.2019.88

CrossRef Full Text | Google Scholar

New, M., Lister, D., Hulme, M., and Makin, I. (2002). A High-Resolution Data Set of Surface Climate over Global Land Areas. Clim. Res. 21, 1–25. doi:10.3354/cr021001

CrossRef Full Text | Google Scholar

Nikishin, A. M., Sobornov, K. O., Prokopiev, A. V., and Frolov, S. V. (2010). Tectonic Evolution of the Siberian Platform during the Vendian and Phanerozoic. Mosc. Univ. Geol. Bull. 65, 1–16. doi:10.3103/s0145875210010011

CrossRef Full Text | Google Scholar

Novis, P., Schallenberg, M., Saulnier-Talbot, É., Kilroy, C., and Reid, M. (2017). The Diatom Lindavia Intermedia Identified as the Producer of Nuisance Pelagic Mucilage in Lakes. New Zealand J. Bot. 55, 479–495. doi:10.1080/0028825x.2017.1377263

CrossRef Full Text | Google Scholar

Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., Mcglinn, D., et al. (2020). Vegan: Community Ecology Package. R Package Version 2.5-7. Available at: https://CRAN.R-project.org/package=vegan.

Google Scholar

Palagushkina, O., Wetterich, S., Biskaborn, B. K., Nazarova, L., Schirrmeister, L., Lenz, J., et al. (2017). Diatom Records and Tephra Mineralogy in Pingo Deposits of Seward Peninsula, Alaska. Palaeogeogr. Palaeoclimatol. Palaeoecol. 479, 1–15. doi:10.1016/j.palaeo.2017.04.006

CrossRef Full Text | Google Scholar

Mountain-Research-Initiative Pepin, N., Bradley, R. S., Diaz, H. F., Baraer, M., Caceres, E. B., Forsythe, N., et al. (2015). Elevation-dependent Warming in Mountain Regions of the World. Nat. Clim. Change 5, 424. doi:10.1038/nclimate2563

CrossRef Full Text | Google Scholar

Pestryakova, L. A., Herzschuh, U., Gorodnichev, R., and Wetterich, S. (2018). The Sensitivity of Diatom Taxa from Yakutian Lakes (north-eastern Siberia) to Electrical Conductivity and Other Environmental Variables. Polar Res. 37, 1485625. doi:10.1080/17518369.2018.1485625

CrossRef Full Text | Google Scholar

Pestryakova, L. A., Herzschuh, U., Wetterich, S., and Ulrich, M. (2012). Present-day Variability and Holocene Dynamics of Permafrost-Affected Lakes in central Yakutia (Eastern Siberia) Inferred from Diatom Records. Quat. Sci. Rev. 51, 56–70. doi:10.1016/j.quascirev.2012.06.020

CrossRef Full Text | Google Scholar

Plikk, A., Engels, S., Luoto, T. P., Nazarova, L., Salonen, J. S., and Helmens, K. F. (2019). Chironomid-based temperature reconstruction for the Eemian Interglacial (MIS 5e) at Sokli, northeast Finland. J. Paleolimnol. 61, 355–371. doi:10.1007/s10933-018-00064-y

CrossRef Full Text | Google Scholar

Popp, S. (2006). “Late Quaternary Environment of Central Yakutia (NE Siberia): Signals in Frozen Ground and Terrestrial Sediments,” (Potsdam, Germany: University Potsdam). Unpublished PhD thesis.

Google Scholar

Quinlan, R., and Smol, J. P. (2001). Setting Minimum Head Capsule Abundance and Taxa Deletion Criteria in Chironomid-Based Inference Models. J. Paleolimnology 26, 327–342. doi:10.1023/a:1017546821591

CrossRef Full Text | Google Scholar

R Core Team (2016). R: A Language and Environment for Statistical Computing.

Reavie, E. D., Sgro, G. V., Estepp, L. R., Bramburger, A. J., Shaw Chraïbi, V. L., Pillsbury, R. W., et al. (2017). Climate Warming and Changes in Cyclotella Sensu Lato in the Laurentian Great Lakes. Limnol. Oceanogr. 62, 768–783. doi:10.1002/lno.10459

CrossRef Full Text | Google Scholar

Reimer, P. J., Austin, W. E. N., Bard, E., Bayliss, A., Blackwell, P. G., Bronk Ramsey, C., et al. (2020). The IntCal20 Northern Hemisphere Radiocarbon Age Calibration Curve (0-55 Cal kBP). Radiocarbon 62, 725–757. doi:10.1017/rdc.2020.41

CrossRef Full Text | Google Scholar

Renssen, H., Seppä, H., Crosta, X., Goosse, H., and Roche, D. M. (2012). Global Characterization of the Holocene Thermal Maximum. Quat. Sci. Rev. 48, 7–19. doi:10.1016/j.quascirev.2012.05.022

CrossRef Full Text | Google Scholar

Renssen, H., Seppa, H., Heiri, O., Roche, D., Goosse, H., and Fichefet, T. (2009). The Spatial and Temporal Complexity of the Holocene thermal Maximum. Nat. Geosci. 2, 410–413. doi:10.1038/ngeo513

CrossRef Full Text | Google Scholar

Rosén, P., Segerström, U., Eriksson, L., and Renberg, I. (2003). Do Diatom, Chironomid, and Pollen Records Consistently Infer Holocene July Air Temperature? A Comparison Using Sediment Cores from Four Alpine Lakes in Northern Sweden. Arctic, Antarctic, Alpine Res. 35, 279–290. doi:10.1657/1523-0430(2003)035[0279:ddcapr]2.0.co;2

CrossRef Full Text | Google Scholar

Rudaya, N., Nazarova, L., Novenko, E., Andreev, A., Kalugin, I., Daryin, A., et al. (2016). Quantitative Reconstructions of Mid- to Late Holocene Climate and Vegetation in the north-eastern Altai Mountains Recorded in lake Teletskoye. Glob. Planet. Change 141, 12–24. doi:10.1016/j.gloplacha.2016.04.002

CrossRef Full Text | Google Scholar

Rühland, K. M., Paterson, A. M., and Smol, J. P. (2015). Lake Diatom Responses to Warming: Reviewing the Evidence. J. Paleolimnol. 54, 1–35. doi:10.1007/s10933-015-9837-3

CrossRef Full Text | Google Scholar

Rühland, K., Paterson, A. M., and Smol, J. P. (2008). Hemispheric-scale Patterns of Climate-Related Shifts in Planktonic Diatoms from North American and European Lakes. Glob. Change Biol. 14, 2740–2754. doi:10.1111/j.1365-2486.2008.01670.x

CrossRef Full Text | Google Scholar

Schlegel, I., and Scheffler, W. (1999). Seasonal Development and Morphological Variability of Cyclotella Ocellata (Bacillariophyceae) in the Eutrophic Lake Dagow (Germany). Int. Rev. hydrobiology 84, 469–478.

Google Scholar

Scussolini, P., Vegas-Vilarrúbia, T., Rull, V., Corella, J. P., Valero-Garcés, B., and Gomà, J. (2011). Middle and Late Holocene Climate Change and Human Impact Inferred from Diatoms, Algae and Aquatic Macrophyte Pollen in Sediments from Lake Montcortès (NE Iberian Peninsula). J. Paleolimnol. 46, 369–385. doi:10.1007/s10933-011-9524-y

CrossRef Full Text | Google Scholar

Self, A. E., Brooks, S. J., Birks, H. J. B., Nazarova, L., Porinchu, D., Odland, A., et al. (2011). The Distribution and Abundance of Chironomids in High-Latitude Eurasian Lakes with Respect to Temperature and Continentality: Development and Application of New Chironomid-Based Climate-Inference Models in Northern Russia. Quat. Sci. Rev. 30, 1122–1141. doi:10.1016/j.quascirev.2011.01.022

CrossRef Full Text | Google Scholar

Serreze, M. C., Barrett, A. P., Slater, A. G., Woodgate, R. A., Aagaard, K., Lammers, R. B., et al. (2006). The Large-Scale Freshwater Cycle of the Arctic. J. Geophys. Res. 111, C11010. doi:10.1029/2005jc003424

CrossRef Full Text | Google Scholar

Shahgedanova, M. (2002). “Climate at Present and in the Historical Past,” in The Physical Geography of Northern Eurasia. Editor M. Shahgedanova (Oxford: Oxford University Press), 70–102.

Google Scholar

Smol, J. P., and Stoermer, E. F. (2010). The Diatoms: Applications for the Environmental and Earth Sciences. Cambridge: Cambridge Univ Pr.

Smol, J. P., Wolfe, A. P., Birks, H. J. B., Douglas, M. S. V., Jones, V. J., Korhola, A., et al. (2005). Climate-driven Regime Shifts in the Biological Communities of Arctic Lakes. Proc. Natl. Acad. Sci. 102, 4397–4402. doi:10.1073/pnas.0500245102

PubMed Abstract | CrossRef Full Text | Google Scholar

Solovieva, N., Klimaschewski, A., Self, A. E., Jones, V. J., Andrén, E., Andreev, A. A., et al. (2015). The Holocene Environmental History of a Small Coastal lake on the north-eastern Kamchatka Peninsula. Glob. Planet. Change 134, 55–66. doi:10.1016/j.gloplacha.2015.06.010

CrossRef Full Text | Google Scholar

Stief, P., Nazarova, L., and De Beer, D. (2005). Chimney Construction byChironomus Ripariuslarvae in Response to Hypoxia: Microbial Implications for Freshwater Sediments. J. North Am. Benthological Soc. 24, 858–871. doi:10.1899/04-137.1

CrossRef Full Text | Google Scholar

Subetto, D. A., Nazarova, L. B., Pestryakova, L. A., Syrykh, L. S., Andronikov, A. V., Biskaborn, B., et al. (2017). Paleolimnological Studies in Russian Northern Eurasia: A Review. Contemp. Probl. Ecol. 10, 327–335. doi:10.1134/s1995425517040102

CrossRef Full Text | Google Scholar

Sundqvist, H. S., Kaufman, D. S., Mckay, N. P., Balascio, N. L., Briner, J. P., Cwynar, L. C., et al. (2014). Arctic Holocene Proxy Climate Database - New Approaches to Assessing Geochronological Accuracy and Encoding Climate Variables. Clim. Past 10, 1605–1631. doi:10.5194/cp-10-1605-2014

CrossRef Full Text | Google Scholar

Syrykh, L. S., Nazarova, L. B., Herzschuh, U., Subetto, D. A., and Grekov, I. M. (2017). Reconstruction of Palaeoecological and Palaeoclimatic Conditions of the Holocene in the South of the Taimyr According to an Analysis of lake Sediments. Contemp. Probl. Ecol. 10, 363–369. doi:10.1134/s1995425517040114

CrossRef Full Text | Google Scholar

Tarasov, P. E., Müller, S., Zech, M., Andreeva, D., Diekmann, B., and Leipe, C. (2013). Last Glacial Vegetation Reconstructions in the Extreme-continental Eastern Asia: Potentials of Pollen and N-Alkane Biomarker Analyses. Quat. Int. 290-291, 253–263. doi:10.1016/j.quaint.2012.04.007

CrossRef Full Text | Google Scholar

Telford, R. J., and Birks, H. J. B. (2011). A Novel Method for Assessing the Statistical Significance of Quantitative Reconstructions Inferred from Biotic Assemblages. Quat. Sci. Rev. 30, 1272–1278. doi:10.1016/j.quascirev.2011.03.002

CrossRef Full Text | Google Scholar

Väliranta, M., Weckström, J., Siitonen, S., Seppä, H., Alkio, J., Juutinen, S., et al. (2011). Holocene Aquatic Ecosystem Change in the Boreal Vegetation Zone of Northern Finland. J. Paleolimnol. 45, 339–352. doi:10.1007/s10933-011-9501-5

CrossRef Full Text | Google Scholar

Van Den Boogaart, K. G., and Tolosana-Delgado, R. (2008). “Compositions”: a Unified R Package to Analyze Compositional Data. Comput. Geosciences 34, 320–338. doi:10.1016/j.cageo.2006.11.017

CrossRef Full Text | Google Scholar

Velichko, A. A., Catto, N., Drenova, A. N., Klimanov, V. A., Kremenetski, K. V., and Nechaev, V. P. (2002). Climate Changes in East Europe and Siberia at the Late Glacial-Holocene Transition. Quat. Int. 91, 75–99. doi:10.1016/s1040-6182(01)00104-5

CrossRef Full Text | Google Scholar

Voigt, C. (2009). “Data Report: Semiquantitative Determination of Detrital Input to ACEX Sites Based on Bulk Sample X-ray Diffraction Data,” in Proceedings of the Integrated Ocean Drilling Program, Volume 302. Editors J. Backman, K. Moran, D. B. Mcinroy, and L. A. Mayer (Edinburgh: Expedition 302 Scientists).

Google Scholar

Vyse, S. A., Herzschuh, U., Andreev, A. A., Pestryakova, L. A., Diekmann, B., Armitage, S. J., et al. (2020). Geochemical and Sedimentological Responses of Arctic Glacial Lake Ilirney, Chukotka (Far East Russia) to Palaeoenvironmental Change since ∼51.8 Ka BP. Quat. Sci. Rev. 247, 106607. doi:10.1016/j.quascirev.2020.106607

CrossRef Full Text | Google Scholar

Vyse, S., Herzschuh, U., Pfalz, G., Pestryakova, L. A., Diekmann, B., Nowaczyk, N., et al. (2021). Sediment and Carbon Accumulation in a Glacial Lake in Chukotka (Arctic Siberia) During the Late Pleistocene and Holocene: Combining Hydroacoustic Profiling and Down-Core Analyses. Biogeosciences. doi:10.5194/bg-2021-39

Wang, R., Zhang, Y., Wünnemann, B., Biskaborn, B. K., Yin, H., Xia, F., et al. (2015). Linkages between Quaternary Climate Change and Sedimentary Processes in Hala Lake, Northern Tibetan Plateau, China. J. Asian Earth Sci. 107, 140–150. doi:10.1016/j.jseaes.2015.04.008

CrossRef Full Text | Google Scholar

Wanner, H., Beer, J., Bütikofer, J., Crowley, T. J., Cubasch, U., Flückiger, J., et al. (2008). Mid- to Late Holocene Climate Change: an Overview. Quat. Sci. Rev. 27, 1791–1828. doi:10.1016/j.quascirev.2008.06.013

CrossRef Full Text | Google Scholar

Wanner, H., Solomina, O., Grosjean, M., Ritz, S. P., and Jetel, M. (2011). Structure and Origin of Holocene Cold Events. Quat. Sci. Rev. 30, 3109–3123. doi:10.1016/j.quascirev.2011.07.010

CrossRef Full Text | Google Scholar

Werner, K., Tarasov, P. E., Andreev, A. A., MãœLler, S., Kienast, F., Zech, M., et al. (2010). A 12.5-kyr History of Vegetation Dynamics and Mire Development with Evidence of Younger Dryas Larch Presence in the Verkhoyansk Mountains, East Siberia, Russia. Boreas 39, 56–68. doi:10.1111/j.1502-3885.2009.00116.x

CrossRef Full Text | Google Scholar

Wetterich, S., Schirrmeister, L., Nazarova, L., Palagushkina, O., Bobrov, A., Pogosyan, L., et al. (2018). Holocene Thermokarst and Pingo Development in the Kolyma Lowland (NE Siberia). Permafrost and Periglac Process 29, 182–198. doi:10.1002/ppp.1979

CrossRef Full Text | Google Scholar

Whitmore, J., Gajewski, K., Sawada, M., Williams, J. W., Shuman, B., Bartlein, P. J., et al. (2005). Modern Pollen Data from North America and Greenland for Multi-Scale Paleoenvironmental Applications. Quat. Sci. Rev. 24, 1828–1848. doi:10.1016/j.quascirev.2005.03.005

CrossRef Full Text | Google Scholar

Wiederholm, T. (1983). Chironomidae of Holarctic Region: Keys and Diagnoses. Part 1. Larvae Entomol. Scand. Suppl. 19, 1–457.

Google Scholar

Wolfe, A. P. (1997). On Diatom Concentrations in lake Sediments: Results from an Inter-laboratory Comparison and Other Tests Performed on a Uniform Sample. J. Paleolimnology 18, 261–268. doi:10.1023/a:1007937300347

CrossRef Full Text | Google Scholar

Wunsam, S., Schmidt, R., and Klee, R. (1995). Cyclotella-taxa (Bacillariophyceae) in Lakes of the Alpine Region and Their Relationship to Environmental Variables. Aquat. Sci. 57, 360–386. doi:10.1007/bf00878399

CrossRef Full Text | Google Scholar

Zhilich, S., Rudaya, N., Krivonogov, S., Nazarova, L., and Pozdnyakov, D. (2017). Environmental Dynamics of the Baraba forest-steppe (Siberia) over the Last 8000 Years and Their Impact on the Types of Economic Life of the Population. Quat. Sci. Rev. 163, 152–161. doi:10.1016/j.quascirev.2017.03.022

CrossRef Full Text | Google Scholar

Zhirkov, I., Trofimova, T., Zhirkov, K., Pestryakova, L., Sobakina, I., and Ivanov, K. (2016). Current Geoecological State of Lake Bolshoe Toko. Int. J. Appl. Fundam. Res. 8, 208–213.

Google Scholar