AMS 14C Dating Problem and High-Resolution Geochemical Record in Manzherok Lake Sediment Core From Siberia: Climatic and Environmental Reconstruction for Northwest Altai Over the Past 1,500 Years

This study presents high-resolution multi-proxy biological and geochemical records in an 82-cm sediment core from Manzherok Lake located in a forest–steppe zone on the western piedmonts of Altai Mountain, Russia. Based on 210Pb dating and 48 accelerator mass spectrometry 14C dates as well as pollen data and geochemical proxies, detailed lake history and local climatic changes over 1,500 years are obtained. Prior to the Medieval Warm Period (MWP), the lake had high productivity under stable moderately wet and warm conditions. During 1,150–1,070 year BP, strong surface runoff led to a high detritus input and an increasing lake level, reflecting the onset of MWP. The lake was deep and fresh under the warm and wet conditions of MWP (1,070–850 year BP). In this interval, more aquatic algae and submerged plants on the lake bottom were generated, which would use dissolved CO2 partially decomposed from organic matters in the deeper sediment layers. Consequently, many acid–base–acid-treated samples contain old carbon influence on their 14C dates. This calls for attention on the chronological construction of lake sediments. During 850–700 year BP, the lake level started to drop with reduced sediment load under cooling and drying conditions. Low sedimentation and lake productivity occurred due to cold and dry climates during 700–500 years BP. Very low sedimentation and hiatus were attributed to ice cover and weak water input between 500 and 50 years BP, corresponding to cold and dry Little Ice Age. Manzherok Lake has recovered productivity and deposition during the current warming century. Change in the total solar irradiance (TSI) is an important factor to influence the climate in the Altai Mountains. With decreased TSI, the Siberian High became strong, which led to the Westerly and the polar front being pushed away from this region, resulting in arid climates. The situation was reversed vice versa.


INTRODUCTION
A small endorheic or closed lake often leads to a high assimilation potential of living matters, which allows the lake to produce high organogenic masses and to form sapropel-type sediments (Korde, 1960;Popolzin, 1967;Lopatko, 1978). The source materials of sapropel formations in a lake normally include the remains of aquatic organisms such as plankton (phytoplankton and zooplankton), algae, macrophytes, and organic matters as well as mineral substances coming from the catchment area. As organic (humus) substances and inorganic mineral impurities (clay and sand) are introduced into the autochthonous organic matter (formed in the lake) under the influence of biochemical, microbiological, and mechanical processes, sapropel of allochthon origin is formed (Korde, 1960;Vine and Tourtelot, 1970;Lopatko, 1978;Poplavko et al., 1978;Gurari and Gavshin, 1981;Neruchev, 1982;Kuzin, 2007). When interpreting geochemical proxies in lake sediments, it is necessary to identify the source and the migration of those elements and isotopes. For instance, Leonova and Bobrov (2012) studied plankton development in reservoirs of the Siberian region. They revealed that, in small lakes, the organic matter of plankton detritus did not significantly change its micronutrient composition while sinking down to the lake bottom. They quantitatively calculated the supply of chemical elements directly through the "plankton channel" into the lake sediments. Vetrov and Kuznetsova (1997) and Granina (2008) studied the geochemistry of diatom pelagic silts and the microelement composition of plankton in Baikal Lake.
In recent years, interesting work has begun to be devoted to the study of the climatic indicator role of chemical elements in sediments (Syso and Yu, 2007;Kalugin et al., 2014;Darin et al., 2014). Syso and Yu (2007) proposed several geochemical indicators for determining the soil properties. It is believed that the enrichment of lake sapropel by mineral components reflects the enhancement of the surface flow of water . Indicators of this phenomena are considered to be "cluster" elements to which Si, Al, Ti, Fe, Mg, Ca, and K belong. The organic fraction of the sediment is more closely related to the increase of temperature at which the bio-productivity of the catchment basin and the water body increases. Some geochemical indices can be directly used for the paleo-reconstruction of the environment. For example, elements with variable valence can be specific indicators of external conditions . Fe is strictly connected with sulfur in sulfate under anoxic conditions. Additionally, Fe forms siderite in carbonate, whereas Mn-enriched layers mark long-term pauses of sedimentation in oxide systems. The Mo/Mn ratio is correlated well with the anoxic environment. According to Kalugin et al. (2014), the Sr/Rb ratio reflects a portion of clay fraction in annually laminated lake deposits of Shira Lake, which marks a stable sedimentation period in winter time. They divided the elements of the bottom sediments of Lake Shira into three clusters: carbonate (Sr, Sc, and Ca), clastic (Zn, Mo, As, Cu, Ni, U, Nb, Ba, Zr, Ga, Mn, Y, Co, Fe, Ti, Rb, and K), and organic pore water (Br, Cl, and S).
The trace metal elements (TME), heavy metal elements (HME), and rare earth elements (REE) in lake sediments are also used for environmental studies. Li (1991) and Leonova et al. (2010) investigated a number of elements that enrich the upper layers of modern lake deposits relatively to clay shale, including their sources and ways to migration in lake sediments. The groups of elements enriching the upper horizons of lake sapropel with respect to clayey shales are singled out as "highly sapropelophilic" elements (the ash concentration coefficients of which KK = 28-15)-P, Br, Mn, As, and Hg; "sapropelophilic" elements (KK = 7-8)-Fe, Mo, Zn, Cd, Cu, Pb, Ag, and Sb; "weakly sapropelophilic" elements (KK = 2-1)-Se and U; and "non-sapropelophilic" elements (KK <1)-alkaline, alkaline earth elements, and rare earth elements (Leonova et al., 2010). In general, it was found that, in small lakes, the plankton contribution significantly exceeds the supply of the microelements from the plant (macrophytic) source, with the exception of Mn. The share of Mn supply through vegetable detritus is about 15 vs. 1% through the plankton supply channel. Plankton is found as a biogeochemical barrier between the atmosphere and the water surface boundary, which intensively retards the chalcophile elements (Cd, As, Sb, Zn, Pb, Hg, Se, etc.) emanating from atmospheric precipitation (Leonova and Bobrov, 2012). Low-biophilic elements of Al, Ti, and Zr are used to characterize the role of clastic and clay material in suspended matter and precipitation (Perelman, 1979).
According to previous studies (Boyarkina et al., 1993;Kutsenogiy and Kutsenogiy, 2000;Raputa et al., 2000;Smolyakov, 2000;Moiseenko et al., 2006), an important additional source of micronutrient supply in small closed lakes is atmospheric mineral matter, and the role of the plankton biofilter is very strongly manifested in such environments (Leonova and Bobrov, 2012). Phyto-and zooplankton are mainly the first substance (biofilter) which involves chemical elements in the biological cycle at the atmosphere-water interface. On this barrier, phyto-and zooplankton get enriched by chalcophile elements which enter the water surface of lakes mainly as part of atmospheric precipitation. These elements are separated into a group of so-called "volatile" elements (Malakhov and Makhonko, 1990;Shotyk et al., 1996;Gavshin et al., 2004;Bobrov, 2007;Dauvalter et al., 2008;Bobrov et al., 2011). While precipitating on the earth's surface as part of atmospheric aerosols, chalcophile "volatile" elements get actively induced in the processes of biodifferentiation by living matter as well as by plankton in small lakes (Leonova et al., 2006;Bobrov et al., 2007). The further fate of these elements is determined by the hydrochemical characteristics of the lake water. Moreover, these elements enter the bottom sediments of the plankton detrital composition if the hydrocarbonate class of water shows a pH range of 7-8 (Leonova et al., 2008Leonova and Bobrov, 2012). In lakes with sulfate class of water with low pH ranging 5-6, these elements would be leached out from plankton detritus into an aqueous solution and get re-introduced into trophic chains. As, Se, Cd, Hg, and Pb are toxic elements with extremely low background concentrations in water. However, they are characterized by high accumulation coefficients in living organisms (Leonova and Bychinsky, 1998;Leonova, 2004;Leonova et al., 2007). Sc, Zr, Nb, Hf, and Th are included in the group of terrigenous elements. Sc is characterized by low solubility, and it has been suggested to use, for standardization, the basic elements to evaluate the accumulation of elements in various water bodies using the "enrichment factor" (Shotyk et al., 1996).
In Blyakharchuk et al. (2017), detailed palaeopollen and diatom data and brief elemental data with 21 accelerator mass spectrometry (AMS) 14 C dates in an 82-cm sediment core from Manzherok Lake located in a forest-steppe zone on the western piedmonts of Altai Mountain were used for climate and environmental reconstruction over the last 1,400 years. Among the 21 14 C dates, nine non-acid-base-acid (ABA)-treated samples are older than their ABA-treated pairs, even though the ABA-treated samples show a rather scattered chronological stratigraphy (Blyakharchuk et al., 2017). Thus, we have dated more layers with ABA treatment. In this study, a total of 48 AMS 14 C dates plus 210 Pb dating results from the 82-cm core yield better chronology. The present study reveals detailed TME, HME, and REE contents in the core. Together with physical, biological, and geochemical data sets, a detailed history of the lake and the environmental conditions under climate influence over the past 1,500 years has been discussed. In addition, archeological activity and human impact in the study area will be addressed.

AREA OF STUDY
Manzherok Lake [51.822 • N, 85.811 • E, 373 m above sea level (a.s.l.)] is a small, freshwater lake located in the western periphery of the Altai Mountains of Southern Siberia on a terrace of the right-hand bend of the Katun' River, 18 km southwest of Gorno-Altaisk City (51 • 49 ′ 15.5" N, 85 • 48 ′ 35.7" E, 423 m a.s.l.) (Figure 1). The lake has an elevation of 373 m, which is much higher than that of Katun' River (297 m a.s.l.). Therefore, although the lake is about 2 km away from the river, there is no evidence that the river has invaded the lake over the late Holocene. Geologically, Manzherok Lake lies within the tectonic unit of Bijsk-Katun' anticlinorium-exhibiting more ancient geological structures of the Altai Early Caledonian fold system. Consequently, in the area of Manzherok Lake, a positive element of relief is presented by ancient carbonate rocks of the Baratalsk series (R3-V), formed in the Proterozoik marine environment (Geology of the USSA, 1997). The low mountains around the lake are covered by forest-steppe vegetation with Betula pendula (BP) and Pinus sylvestris, with patches of Pinus sibirica. The climate is moist with a relatively mild temperature even though it is classified as of the West Siberian continental type (Ogureeva, 1980). A more detailed regional setting with detailed pollen and diatom data of Manzherok Lake is presented in an earlier publication (Blyakharchuk et al., 2017). Although the origin and the history of Manzherok is still under study, Rusanov and Vazhov (2017) provided the most comprehensive summary of the Manzherok Lake studies. The origin of Lake Manzherok was probably caused by a catastrophic breakthrough of the large glacial-dammed lakes of Altai at the end of the Ice Age. Fed by streams flowing from the surrounding slopes, precipitation, and groundwater, the lake is elliptical in shape and elongated from southwest to northeast. Manzherok Lake is a freshwater open lake, with an outlet at the southwest end through a swampy runoff hollow. The lake water has low pH of 6.2-7.2 and contains low ionic concentrations. The salinity and the alkalinity of the lake are quite low, e.g., CO 2− 3 and HCO − 3 measurements in the lake in 1972 were 0.0 (non-detectable) and 0.28 mg/100 g water, respectively (Ilyin, 1982). The modern Manzherok Lake is a popular recreation area. Human activities around the lake, such as fishing, farming, and grazing, may have impact on the sedimentary processes to the lake.

Elemental Geochemical Analyses of Lake Sediment Samples
About 0.05 g of dry powder lake sediment from every 1-cm horizon of the upper 31-cm section and every 2-cm interval of the 32-82-cm section of the core was precisely weighted and then digested with a mixed acid solution (HF + HNO 3 + HCl) in a polytetrafluoroethylene beaker using a microwave system (SpeedWave 3, Berghof, Germany). The microwave program ran for 15 min at about 145 • C, for 20 min at 200 • C, and a holding time of 15 min, followed by a cool-down procedure of 10 min. After the residue was completely diluted, the evaporate was dissolved repeatedly with concentrated HCl to convert the solution in chloride medium. Then, the solution was diluted to 100 ml volume with 0.5 N HCl. The acid solution was filtered through an acetate filter membrane with 0.45-µm pore size to remove any undissolved particles. The filtered solution was analyzed for 45 elements using inductively coupled plasma-mass spectrometry (Agilent Technologies 7700×, Japan). Multi-element standards (Agilent) were routinely analyzed as quality control. The analytical performance was assessed through related materials including lake Baikal sediment reference material BIL-1 (Russia) and LKSD-1 (Canada), which all have certified or recommended values. The relative standard deviation was <8%. The recovery of the reference materials ranges 91-112%. A total of 54 samples have been analyzed for TME, HME, and REE.

Pb Dating and Sedimentation Models
Lead-210 measurement can be achieved by either directly using low-background gamma spectrometry or measuring its decay product ( 210 Po) via alpha spectrometry. Despite being proven efficient for low-density samples (e.g., peat, lake sediment), gamma spectrometry is limited in small samples (e.g., Ebaid and Khater, 2006), while the measurement of low-energy 210 Pb gamma photons (E 210Pb = 46.5 keV) is virtually impossible. The 210 Pb dating of the Manzherok Lake core was determined by the 210 Po method (Pawlyta et al., 2004;Cooke et al., 2007;Baskaran et al., 2014).
Dry sediment samples of 1.0 g for the top core section and 2.0 g for the bottom position with a known amount of 209 Po spike solution were placed into Teflon vessels and digested at a temperature of 120 • C using concentrated HNO 3 + HClO 4 + HF. After 12 h of digestion, the solution was centrifuged. The supernatant was transferred into a Teflon beaker and evaporated with 6 M HCl until dryness. The evaporate in the Teflon beaker was dissolved in 10 ml 0.5 N HCl and transferred into a 50ml centrifuge tube. Ascorbic acid powder was added into the solution to form Fe complex ions in order to prevent the Fe ions to be co-precipitated with Po ions. A silver disk, 8 mm in diameter, was placed in the solution of the centrifuge tube. Then, the centrifuge tube was placed in a water bath at 60 • C. Polonium isotopes were spontaneously deposited within 4 h on the Ag disk. The activities of 209 Po and 210 Po were measured by an alpha spectrometer. Two blank samples were analyzed with each sample batch to verify the quality of the chemical extraction. The procedure of acid extraction and deposition of polonium was outlined in Suriyanarayanan et al. (2008) and Ugur et al. (2003). The total 210 Pb activity reached a relatively constant below 13 cm in depth (Figure 2). We use the average 210 Pb activity below the 13-cm depth as the supported 210 Pb activity. The excessive 210 Pb ( 210 Pb ex ) in each sample was obtained by the measured total activity of 210 Pb subtracted by the supported 210 Pb activity (Figure 2). The Constant Rate of Supply model (Appleby, 2001) was applied to the 210 Pb inventories calculated from the 210 Pb ex data to generate ages.

AMS Radiocarbon Dating
In Blyakharchuk et al. (2017), we reported 21 AMS 14 C dates from the core. In that study, we have found that the lake sediments need ABA treatment (Brock et al., 2010). Otherwise, the 14 C dates became older than their true ages. Therefore, in the present study, about 0.1 g of bulk gyttja sample from different depths of the core was treated with 0.5 N HCl, 0.5 M NaOH, and 0.5 N HCl sequentially. The dry-treated sample was then placed into a 9-mm quartz tube with pre-combusted CuO powder and a piece of silver and then placed on a vacuum line. The quartz tube was sealed under vacuum of 1e-5 mbar and then combusted for 6 h at 850 • C in a muffle furnace. The CO 2 produced by oxidation was transferred and purified cryogenically on the vacuum line and sealed into combination tubes, which included a 9-mm glass tube with Zn and TiH 2 powders and an inner 6-mm center tube containing Fe powder (Xu et al., 2007). Graphitization of the CO 2 in the tube took place in the muffle furnace at 550 • C for 6 h.
Sample graphite was pressed into a target holder and measured for its 14 C/ 12 C and 13 C/ 12 C ratios with a 1.0 MV Tandetron Model 4110 BO-Accelerator Mass Spectrometer in the NTUAMS Lab. Every sample batch contains at least three international standards (OXII, 4,900 C), three backgrounds (BKG), and two inter-comparison samples (IRIs). The measurement mode was 14 C 3+ to avoid 2Li + interference with 14 C 2+ . Using the 14 C/ 12 C and the 13 C/ 12 C ratios of OXII, BKG, and the samples, the percentages of modern carbon (pMC), D 14 C (= pMC/100-1) × 1,000, and 14 C conventional age were calculated with the Libby half-life of 5,568 years (Stuiver and Polach, 1977). δ( 13 C) is calculated from the 13 C/ 12 C ratios of OXII (δ 13 C = −18‰) and sample, which is used for the correction of carbon isotopic fractionation (against −25‰) during natural and AMS dating laboratory processes. Hence, the δ( 13 C) value is not only determined by the carbon isotopic composition of the sample but also strongly influenced by the ABA treatment and the AMS measurement. Measured by AMS, the δ( 13 C) is different from δ 13 C (measured by isotope-ratio mass spectrometry) and cannot be used for stable isotope interpretation. The conventional age is transferred into calibrated 14 C ages with 1σ error using the calibration curve IntCal13 (Stuiver and Reimer, 1993;Niu et al., 2013;Reimer et al., 2013). The 14 C ages expressed in this paper are calibrated 14 C ages in years BP (0 year BP = 1,950 CE). A total of 48 AMS 14 C dates from the core are listed in Table 1.

RESULTS AND DISCUSSIONS
Problems of 14 C Dates and New Chronology Table 1 and Figure 3 show all the 14 C dates. First of all, the core contains very high organic content (>20 wt.%) ( Figure 3B). The core materials have black color with swamp odor, which look totally different from the river sediments and the surrounding debris. The latter ones contain much lower organic matter. The materials in the core are very suitable for 14 C dating, especially the lower part which contains organic matter like peat as high as 35 wt.%. However, these organic matters are well-decomposed and difficult to be isolated from the bulk sediments. In the earlier publication (Blyakharchuk et al., 2017), several plant remain samples were AMS 14 C dated, even though those plant remains still show old carbon influence on the non-ABA treated 14 C dates ( Table 1). In fact, all non-ABA-treated samples have older ages than those of ABA-treated samples from the same bulk sample pairs. Since the remaining core materials could not pick up more plant remains, in this study we have dated bulk gyttja (sapropel) samples from multiple horizons above 35 cm in depth, with ABA pre-treatment. As the lake sediments were deposited stratigraphically, their ages should become older from top to bottom. In Table 1, the first criterion for selecting the corrected dates is to remove all the non-ABA treated samples. The second criterion is to remove an older date (ABA-treated) in the upper layer as the younger dates contain least old carbon influence. If an age in the upper layer is slightly older that in the level below but is within uncertainties, we will keep it.
Although the lake has low pH (6.2-7.2), the carbonate content in the bulk sample still has about 2-5% (Figure 3). Those carbonates contain normally old carbon from detritus carried by surface runoff. With ABA treatment, those carbonates can be removed. Another older carbon source comes from organic carbon such as humic acids and/or CO 2 decomposed from deeper sediments. The humic acids decomposed from older organic matter in the deeper sediments may migrate upward and enter the bulk gyttja. The function of base treatment in the ABA procedure is to remove such mobile humic acids. Thus, the ABA-treated samples should be able to eliminate the influence of carbonate and mobile humic acid. However, there are many newly measured ABA-treated samples showing reversed 14 C ages (Table 1 and Figure 3). Our hypothesis of this old carbon influence is that dissolved CO 2 in the lake water contained CO 2 decomposed from organic matters in the deeper sediments that was uptaken by submerged plants and aquatic algae to form organic compounds which could not be removed by ABA treatment.
Since the CO 2 decomposed from organic matters in the deeper sediments is a function of lake chemistry such as redox condition, water temperature, pH, and bacterial activity, etc., the later ones are affected by water depth and climatic condition. Thus, this old carbon influence is different from "hard water effect" or "reservoir effect." "Hard water effect" is caused by high concentrations of HCO − 3 and CO 2− 3 in lake water, which commonly exist in closed alkaline lakes. Manzherok Lake does not belong to such a lake. "Reservoir effect" including hard water effect sometimes can be caused by organic carbon input from inflow sediments, for example, Teletskoye Lake in Altai has reservoir age in the 14 C ages of bulk organic carbon in the sediments (Rudaya et al., 2016). This lake is also a freshwater lake, but its sediments contain a low percentage of total organic carbon. A total of 16 AMS 14 C dates (from two labs) on the TOC of bulk sediments from the lake core show about 2,400-3,000 years of "reservoir age" throughout the sediment core. This is not the case for Manzherok Lake core. In the Manzherok Lake core, there are several 14 C dates from plant remains. Even on the bulk gyttja samples from the upper 5-cm part, four ABAtreated samples contain nuclear bomb 14 C signal, indicating no reservoir effect (Table 1). Thus, the old carbon influence on the 14 C age of the Manzherok Lake sediments is not considered as "hard water effect" or "reservoir effect." This influence varies with lake condition with time. Figure 3A shows that the old carbon influence was negligible above 19 cm in depth and was relatively weak below 46 cm in depth.
In general, deeper water (higher lake level) is in favor of the anoxic condition in the bottom of the lake. Warmer lake water provides a stronger bacterial activity. The reduced environment would generate CH 4 , which would be oxidized into CO 2 during degassing in the lake bottom. In Figure 3, most ABA-treated samples with rejected 14 C ages are in 40-20-cm depths. This interval belongs to the Medieval Warm Period (MWP). In later  sections, we will interpret that this interval had a higher lake level under warm and wet climatic conditions. Since the lake sediments should be deposited in age sequence, those reversed ages could not be used for chronology construction. Based on the criteria described above, we remain as many reasonable 14 C dates as possible for the Bacon age model (Blaauw and Christen, 2011). Figure 4 exhibits the Bacon age model of the Manzherok Lake core, which shows the core to contain deposition from 1,440 years BP to the present. Our findings about the problems of 14 C dating on the Manzherok Lake sediments call for attention about chronological construction: (1) it is necessary to perform ABA treatment for lake sediments and (2) ABA treatment cannot remove old carbon influence if uptake is older dissolved CO 2 (non-equilibrium with the atmospheric CO 2 ) by organisms in the lake. In freshwater lakes, terrestrial plant remains in lake cores are often difficult to be obtained for dating. Old carbon influence and age reversion are common, e.g., in Rudaya et al. (2016). Sometimes, with a few 14 C dates, the chronology may be not good even if the dates are in a stratigraphic order.
Although the new Bacon age model provides a better chronology than that in Blyakharchuk et al. (2017), it has a problem in the upper 5-cm part. The Bacon age model results yielded ages in the upper 5-cm part of 0-305 years BP. However, four 14 C dates in the upper 5 cm show nuclear bomb 14 C, reflecting that the lake sediments in this interval must be deposited after 1,950 CE. The 14 C ages from 5-to 9-cm depths quickly jumped to 390-580 years BP, indicating that the lake deposition between 5-and 7-cm depths could have hiatus.
Because the Bacon age model does not allow sedimentary hiatus, the modeled result forces the ages in the upper 5 cm to become older. Besides that, the 210 Pb dating result also indicates that the upper 10-cm part has an apparent 210 Pb ex decay trend, reflecting modern deposition (Figure 2). Therefore, the results of the Bacon age model for the upper 5-cm part were not used. We instead use a linear sedimentation rate given by the 210 Pb dating for the upper 5-cm part. The reason that we do not use the 210 Pb ages shown in Figure 2 is that 210 Pb in the surface sediments of a freshwater lake can diffuse downward (Benoit and Hemond, 1991), whereas the majority of the organic carbon (plant remains) in the lake sediments are not mobile. Hence, we use the agedepth relationship of the Bacon age model below the 5-cm depth. In the new chronology, there is only one data point (1-cm layer) between 0 year BP (4.5-cm depth) and 465 years BP (6.5-cm depth). This interval represents the depositional period of the Little Ice Age. Nevertheless, the new chronology is significantly improved by more 14 C dates and 210 Pb. We shall re-interpret the Manzherok Lake record based on the new chronology.

Spore-Pollen Data for Manzherok Lake
A detailed spore-pollen diagram, with descriptions of the pollen zones and the phases of vegetation development in the study area obtained from the Lake Manzherok core, was published earlier (Blyakharchuk et al., 2017). According to the new chronology, the time boundaries of the phases in vegetation development reconstructed based on the pollen zones have changed slightly (Figure 5), which can be summarized as follows: Based on the Manzherok pollen data, we reconstruct climate changes in terms of humidity and compare with neighboring palaeoecological data of Teletskoye Lake record (Rudaya et al., 2016) (Figure 6). Although the new chronology practically has no significant influence on the interpretations which were made in Blyakharchuk et al. (2017), some differences are addressed below: 1. According to the modified chronology, the period of humid climates during medieval episode became shorter than that in the previous chronology (200 years compared with 400 years). Consequently, it can be concluded that the climate of the medieval climatic anomaly in the western foothills of the Altai was not uniform in terms of the degree of wetness. The early Middle Ages (1,050-900 years BP) was wet, which is marked by the spread of dark coniferous  and pine forests in the vicinity of the Manzherok Lake. The late Middle Ages (900-700 years BP) was characterized by relatively dry climate and domination of the birch forest-steppe.
2. The Mongolian invasion (1,236-1,242 years CE) in the territory of Altai to the new chronology is synchronous with the layer of lake sediment at the depth of 11-13 cm or ∼710 years BP (25-30 cm in the old version). During this period, the climate became less arid than before, and the area occupied by wormwood steppes in the foothills of the Altai decreased. Below this layer of sediment (850-750 years BP), a new series of geochemical data shows two high maxima of Zn and an increase of Cu and Fe content. In the pollen spectrum of this layer, the abundance of pollen of anthropogenic indicators and charcoal was slightly decreased. These indicators confirm our earlier interpretation of the influence of Mongolian invasion on the Altai region (Blyakharchuk et al., 2017). 3. During the cold Little Ice Age (LIA), according to the modified chronology, the rate of sediment accumulation in Manzherok Lake was very low or even absent. The cold and dry climates during LIA could lead to a longer period of frozen lake surface and very low biological productivity in the lake. Low surface runoff under the dry climates reduced the sediment input into the lake. Unfortunately, based on the new chronology, there is only a 2-cm layer falling into LIA. Thus, it is not possible to trace the detailed effect of LIA climates on the lake deposition by biological and geochemical proxies. Besides the very low sedimentation which indicates cold/dry climatic condition, our interpretation on the biological and geochemical results of the 5-7-cm sediment sample deposited during the LIA should be with caution.

Results of High-Resolution Geochemical Analysis
In this study, a new and detailed series of geochemical data ( Table 2) was obtained for samples of the same core for which spore-pollen, diatom, and preliminary geochemical analyses were made earlier (Blyakharchuk et al., 2017). A total of 54 samples (every 1 cm in the upper 31-cm part and every 2 cm below the 31-cm depth) were used for high-resolution geochemical analysis. A total of 46 elements, including TME, HME, and REE, were analyzed for every sample. It is important to point out that these elements are from the bulk sediment samples which were completely dissolved in HF + HNO 3 + HCl acids by microwave digestion. Therefore, the majority of the elements should come from detritus, mainly from exogenous origins. Their variations probably are mainly due to input change instead of lake chemistry, except for some elements such as alkaline and alkaline earth elements. Figure 7 plots the variations of selected elements in the core. Among the measured elements, Al, Fe, and Ti are the dominant elements and have high concentrations in the order of mg/g, whereas the other elements are on the orders of mg/kg. Al, Fe, and Ti mainly come from detritus and have a strong positive correlation (R 2 > 0.9) with the ash content of the sediment which was determined by loss of ignition at 550 • C (Blyakharchuk et al., 2017). Some trace metal and heavy metal elements such as Li, Rb, Th, Nb, Zr, Cr, V, Hf, Sc, and Y have similar patterns as those of Al, Fe, and Ti, indicating that they have the same origins. Together Al, Fe, Ti, Li, Rb, Th, Nb, Zr, Cr, V, Hf, Sc, and Y are considered as lithophilic elements, and their variations are chiefly controlled by detrital input. Furthermore, the REE contents are also strongly correlated with the ash content (R 2 = 0.88) (Figure 8). Thus, the detrital content of the sediments is the main dominant factor to control the variations of the abovementioned elements (Dobrovolsky, 1999). Hence, we can classify the geochemical variation of the core as follows: Zone I: 1,440-1,150 years BP (82-50 cm): high organic (OM) and CaCO 3 contents, low ash content, and low TME, HME, and REE contents. Zone II: 1,150-1,070 years BP (50-40 cm): strongly increased ash content but decreased OM and CaCO 3 contents and strongly increased TME, HME, and REE contents. Zone III: 1,070-850 years BP (40-17 cm): maximum values of ash, lithophilic element, and REE contents, but low OM and CaCO 3 contents in the beginning. Then, the lithophilic element and the REE contents gradually decreased. Zone IV: 850-500 years BP (17-7 cm): the lithophilic element and the REE contents were stable at moderate levels. The sedimentation rate had a significant decrease. Zone V: 500-50 years BP (7-5 cm): very low sedimentation rate; CaCO 3 content increased. Except As and Pb, most elemental concentrations were relatively low. Zone VI: The last century (5-0 cm): OM and CaCO 3 contents increase. Most elemental concentrations are reduced to low values.
The geochemical records shown in Figures 7, 8 provide us the lake history under climate influence. During zone I, the lake had relatively high biological productivity as reflected by high OM and CaCO 3 . The lake was a freshwater lake with a moderate lake level because the organic matter was generally high (28-35%) and the CaCO 3 content was only 2-5% in the bulk sediments. With relatively low detritus content, the lake sediments contained low TME, HME, and REE. During zone II, as the surface runoff sharply increased, the OM and the CaCO 3 contents rapidly decreased and the ash contents strongly increased, which resulted in strong increases of TME, HME, and REE. During zone III, the surface runoff could be still high. However, as the lake became deeper and larger, sediment loading in the coring site became finer and detritus/OM ratio gradually decreased. The lake level was still high in the beginning of zone IV. Then, as the climate became cooling and drying, the lake input decreased, resulting to lowering of the sedimentation rate during zone IV. The lake had very low sedimentation during zone V due to cold and dry climatic conditions. Until about 100 years ago, the lake started to recover its biological productivity and water input. Rare earth elements include La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Tu, Yb, and Lu. Light REEs involve La to Gd, and heavy REEs contain Tb to Lu. Many previous studies have used REE patterns to trace the depositional processes and identify the sediment provenances since they behave conservatively during sediment formation (Prajith et al., 2015;Xu et al., 2017). However, most REE studies of sediment provenances are used in marine environments. For the Manzherok Lake core, the sediment source is relatively simple and the depositional time is relatively short. Therefore, we only use the REE pattern to discuss the sediment feature throughout the core. Figure 9 shows the REE normalized to North American Shale Composite (NASC) and Figure 10 reveals the REE normalized to chondrite in different sediment layers. We averaged the patterns    Frontiers in Earth Science | www.frontiersin.org of 51-82 cm as their ratios in the different layers of this zone are quite similar. In the two figures, although the patterns are all similar, the values of the normalized ratios can be grouped into three categories: (1) 51-82 cm, (2) above 29-33-cm layer, and (3) between 34-35-cm layer and 49-50-cm layer. In the first category which is corresponding to zone I, both NASC and chondrite normalized ratios are the lowest, and the ratios of Ho, Er, Tm, Yb, and Lu to NASC are the same (Figure 9). The feature of the first category can be considered as the organic-rich sediments with a relatively constant input, which can be regarded as an endmember of the mixed lake sediments with another end-member from the surface runoff. In the second category, the values of both NASC and chondrite normalized patterns increased from the 49-50-cm layer to the 38-39-cm layer and then slightly decreased to the 34-35-cm layer (Figures 9, 10). In principle, the REEs in the lake sediments are mainly from detrital input carried by surface runoff. When the detrital input increases with the surface runoff under wet climates, the REE contents should increase. The patterns at different layers are basically the same, indicating that the detrital source, which is from the surrounding area of the lake, remains relatively constant. The variations in the NASC and chondrite normalized ratios shown in Figures 9, 10 depend on the percentage of the input detritus in the lake sediments. The higher ratios reflect a higher percentage of the input detritus, which is confirmed by the increasing ash content from the 49-50-cm layer to 38-39-cm layer. Therefore, the changes in the geochemical proxies during the depositional period of 50-40 cm in the lake core (zone II, 1,150-1,070 years BP) reflected strongly the increased detrital input due to the enhanced surface runoff under wet climates. The third category shows that the patterns and the ratios are similar, especially for the chondrite normalized ratios (Figure 10). This means that the amount of detrital input to the lake sediments at the coring site from the surface runoff reduced after deposition at 33-cm depth (1,005 years BP). Such a situation could be caused by both increased lake level, so that it was difficult for the sediments from the surface runoff to reach the depo-center where the core was retrieved, and reduced surface runoff. The variation trends of As, Cd (not shown), and Pb are different from those of lithophilic elements (Figures 7, 8). Except  Frontiers in Earth Science | www.frontiersin.org at the surface 5 cm, these elements have a general increasing trend from the bottom to the upper section. Figure 11 exhibits that those toxic elements have some positive correlations (R 2 = 0.39 for Pb vs. Cd and R 2 = 0.46 for Pb vs. As). Normally, these elements have low concentration in the natural background so that their origins should not come from detrital materials. The accumulation of these elements in the lake sediments may be related to the metabolic processes of living organisms and anthropogenic activities.
Sr, Ba, and Rb are IIA elements in the periodic table and have similar chemical characteristics. However, Sr is more influenced by lake water chemistry, whereas Ba and Rb are mainly from detritus. Figure 11 shows the significant correlation (R 2 = 0.735) between Rb and Ba, but weak or no correlation (R 2 = 0.15) between Rb and Sr. Using the Sr/Ba and the Sr/Rb ratios, one may eliminate detrital influence so that the effect of lake chemistry can be seen. Sr prefers to co-precipitate with carbonate and favors high alkalinity and lake productivity. The variations of Sr/Ba and Sr/Rb in the core show peaks at 1,410, 1,335, 1,175, and 725 year BP, which also correspond to peaks in CaCO 3 (Figure 8). During those periods, the lake had relatively high productivity and alkalinity. A slight increase in Sr/Ba, Sr/Rb, and CaCO 3 during LIA might not reflect an increase in lake productivity. Under dry climatic condition during LIA, the lake level or volume decreased so that the concentrations of Ca 2+ and CO 2− 3 in the lake might be enhanced to raise alkalinity and pH, resulting in carbonate precipitation.
The Mo/Mn ratio is correlated well with the anoxic environment , whereas the Mn-enriched layer marks a long-term pause of sedimentation in the oxidized systems. In Figure 8, the Mo/Mn peaks appeared at 1,335, 1,270, 1,140, 930, and 630 years BP, indicating the anoxic conditions in the lake, which reflected high lake productivity under relatively warm and moist conditions on decadal time scales in the study area.
These periods also coincide with the mass bloom of algae and cyanobacteria Anabaena. Thus, the geochemical proxies and the biological proxies agree with each other, confirming the observations.

BIOLOGICAL AND GEOCHEMICAL PROXIES IN THE MANZHEROK LAKE CORE AS INDICATORS OF CLIMATIC AND ENVIRONMENTAL CHANGE OVER THE PAST 1,500 YEARS
The biological and the geochemical records of the Manzherok Lake core with the new chronology provide us with details on lake history and vegetation change under climate regimes and human impacts since 1,440 years BP. In Figure 12, we summarize the major geochemical proxies in the core and compare them with the total solar irradiance (Bard et al., 2000;Steinhilber et al., 2009). The OM, ash, and CaCO 3 contents reflect the sediment feature which briefly dominate the lake geochemical and productivity signals. The concentrations of total REE, Al, and Ti represent detrital inputs. Sr/Ba and Sr/Rb can be considered as indicators of lake alkalinity and productivity. Mo/Mn may illustrate anoxic conditions. The shortage of this record is that the sedimentation rate was too low during the Little Ice Age so that the interpretation for this period is brief. In the following, we shall describe the record for each lake and climate stages one by one: 1,440-1,150 years BP (zone I for geochemical proxies and zone 1 for pollen records): This period is prior to the Medieval Warm Period. During this period, the lake had a moderately high level, with relatively high productivity and alkalinity under a frequently anoxic condition, which was shown by high organic matter and carbonate contents, peaks of Sr/Ba, Sr/Rb, and Mo/Mn, and mass bloom of algae and cyanobacteria Anabaena. Vegetation in the area was a phase of birch forest-steppe with a high abundance FIGURE 11 | Left: correlations among toxic elements of Pb, Cd, and As. Right: relationships among Sr, Ba, and Rb.
Frontiers in Earth Science | www.frontiersin.org FIGURE 12 | Comparison of important geochemical indices in the Manzherok Lake core with the total solar irradiance (TSI). The TSI curve is a five-point running average of the 5-years resolution TSI data from Steinhilber et al. (2009). Numerical numbers 1, 2, and 3 for the Sr/Ba and the TSI plots denote their correlations.
of Nymphaea and B. pendula and relatively low P. sylvestris. The climatic conditions were moderately warm and wet. Surface runoff to the lake was not strong so that detrital input from the surrounding lake was low. Decadal scales of warm but dry episodes appeared at 1,335, 1,270, and 1,140 years BP.
1,150-1,070 years BP (zone II for geochemical proxies and zone 2 for pollen records): Manzherok Lake had a sudden change due to wet climates. The surface runoff strongly increased and brought high sediment input to the lake, shown by strongly decreased OM and CaCO 3 and strongly increased ash content, lithophilic elements, and REE. The lake level increased so that lake productivity and alkalinity decreased (Sr/Ba and Sr/Rb strongly dropped) and the anoxic condition (Mo/Mn peak) was further enhanced. The REE concentrations continuously increased and reached a maximum at 1,070 years BP. The vegetation in the area became transitional phase from birch forest-steppe to birch and pine forests. The abundance of B. pendula sharply decreased and Nymphaea disappeared, while P. sylvestris increased. In the lake, algae and Anabaena were abundant, reflecting that the lake was still productive. The wet and warm climates marked the onset of MWP.
Anthropogenic influence on the landscape began. Approximately from 1,130 years BP or slightly earlier, one can see permanent findings of the pollen of cultivated plant Triticum supplemented by the pollen of field weeds. Simultaneously, we count abundant microcharcoal in pollen slides. It is evidence of spreading of fire-cutting agriculture on the western piedmonts of the Altai Mountains. According to these evidences, fire-cutting agriculture began to spread here a few centuries before the Russian colonization. It is worth to mention that human activity could cause the loosening of the earth surface, which might introduce more sediment input to the lake.
1,070-850 years BP (zone III for geochemical proxies and zone 3 for pollen records): This interval had warm and wet climates corresponding to MWP. Manzherok Lake was the largest and deepest during this period over the past 1,500 years. Due to the deepening of the lake throughout high and frequent surface runoffs in the previous stage, the sediments at the depo-center, where the core was retrieved, contained less detritus as shown by reduced contents of ash, REE, and lithophilic elements. In the early half of this period, the lake was very fresh and had low productivity (very low Sr/Ba and Sr/Rb). However, during the late half of the period, the lake productivity and the alkalinity (increased Sr/Ba, Sr/Rb, and OM content) increased, and the anoxic condition (highest Mo/Mn peak) at the bottom of the lake became very strong. The vegetation in the area was a phase of birch and pine forests during this period, reflecting warm and moist conditions. The abundance of P. sylvestris reached its maximum (ca. 40%), B. pendula dropped to its minimum (∼20%), and Nymphaea disappeared, showing the dynamic change of the local vegetation. The enhanced Bryales/algae and anabaena abundances and the REE and lithophilic elements before the end of this stage indicate that the lake started to drop its level (Figure 8). The humidity index, as reconstructed by the pollen data using the equation of Wang et al. (2010), provided a semi-quantitative moisture change (Blyakharchuk et al., 2017), showing the highest humidity at around 950-900 years BP (Figure 6). After 900 years BP, the climate became dry.
Owing to the deep and anoxic conditions of Manzherok Lake during MWP, the submerged plants and algae would uptake dissolved CO 2 which could not fully exchange with the atmospheric CO 2 . Because the dissolved CO 2 contained partial CO 2 which was decomposed from old organic matter in the deeper sediments, its 14 C had an initial age. Therefore, even though many samples from this interval were gone through ABA treatment, their ages are still older than their true ages because the ABA treatment could not remove the old carbon influence caused by the uptake of dissolved CO 2 .
The anthropogenic influence during this period became more significant as shown by the accumulation of toxic elements including Pb, As, and Cd. Metal elements such as Ni, Cu, and Zn commonly exit in the copper ore. The joint finding of Zn and Cu in the bones of people was considered as an indicator of their involvement in the metallurgical process (Aleksandrovskaya and Aleksandrovsky, 2007). The same authors pointed out that during the smelting of copper from copper and polymetallic ores and even during the remelting of copper objects, Zn could evaporate together with metallurgical gases. When the metallurgical gases contacted with cold air, Zn could condense and fall on the ground. The maximum Cu concentration appeared at 980 years BP, probably reflected as ancient smelting in the local area.
850-500 years BP (zone IV for geochemical proxies and lower zone 4 for pollen records): The vegetation after 850 years BP was determined as a phase of pine-birch foreststeppe. The lake sedimentation rate began to drop rapidly due to drying and cooling climatic conditions. The humidity index based on the pollen data dropped significantly and was kept at a low level (Figure 6). The concentrations of REEs and lithophilic elements in this interval was kept relatively constant at lower values than those in zone III, indicating a lower detrital input corresponding to reduced sedimentation rate under a dry climate. The lake level declined gradually so that carbonate precipitation slightly increased at 800 years BP (Figure 12). A slight increase in REE and lithophilic elements at 750 years BP probably indicated that the lake reached a low steady level with moderate productivity. The lake enhanced its alkalinity at 725 years BP as reflected by the peaks of Sr/Ba and Sr/Rb due to the further shrinking of the lake volume. As the lake volume shrank, its nutrient concentration increased so that the lake productivity also slightly increased, as shown by the enhanced abundance of algae and Anabaena between 700 and 600 years BP (Figure 6). A weak but apparent anoxic condition of the lake occurred at around 630 years BP. After 600 years BP, the climatic condition in the study area entered into cold LIA.
According to history record, the first Mongolian invasion to Siberia was 1,207 CE (743 years BP). The time period that Mongolia conquered the Altai region was 1,236-1,242 years CE (714-708 years BP). Thus, the two Zn peaks, together with a slight increase in the Cu and the Ba contents at 850-830 and 770-750 years BP, in our record might reflect the metallurgical practice of local human tribes who lived on the western piedmonts of the Altai Mountains in XII century CE. The humidity index based on the pollen record and the geochemical proxies of the Manzherok core illustrate that the climate was moist during the beginning of the Mongolian invasion, although the climate during this period was not as wet as that during MWP. The Mongol Empire can be divided in two parts: 1,206 CE (744 years BP) to 1,270 CE (680 years BP) before the Yuan Dynasty and 1,271 CE (679 years BP) to 1,368 CE (582 years BP) of the Yuan Dynasty. The first part, especially during the Chenghis Khan reign, was mainly active in the north regions including Altai, Siberia, whereas the second part was mainly active in Mongolia and northern China. According to the Manzherok Lake record, the climate in the study area began cooling and drying at around 750 years BP and became much worse after 600 years BP. The climatic conditions certainly affected human activity in the study area. 500-50 years BP (zone V for geochemical proxies and middle zone 4 for pollen records): The lake depositional record in this interval corresponds to LIA. The sedimentation rate was very low and often had depositional hiatuses due to cold and dry conditions during this period. For only the 2-cm sediments in this interval, the pollen and geochemical data were not good enough to interpret the detailed lake history. Nevertheless, the extreme sedimentation was attributed not only to low surface runoff under dry climates but also to longer frozen time under cold conditions.
After 1,950 CE (zone VI for geochemical proxies and upper zone 4 for pollen records): The top 5-cm sediments containing nuclear bomb 14 C indicate that they were deposited after 1,950 CE. In the 5-cm layer, one can still see clear trends of increasing OM, decreasing ash content, and declining concentrations of REE and lithophilic elements. These trends reflect the recovery of lake productivity and deposition in the current warming centenary.
The above description of Manzherok Lake history and the local climatic conditions over the past 1,500 years refine our findings in Blyakharchuk et al. (2017). The detailed pollen data and the geochemical proxies consistently illustrate the climate changes well. We have compared this record with the pollen and geochemical records from Teletskoye Lake which is 100 km east during the past 1,000 years Kalugin et al., 2007). The climate reconstruction-based pollen record (Figure 6) in Teletskoye Lake during the last 1,000 years  are summarized as follows: (1) humid climate of late MWP 1,000-880 years BP and dry climate 880-750 years BP, (2) humid climate 750-550 years BP, and (3) dry climate of LIA 520-110 years BP. Our Manzherok Lake record agrees well with the above-mentioned climate reconstructions. We have also found a good correlation of Sr/Rb index between the Manzherok Lake record and the Teletskoye Lake record . However, we should mention that Manzherok Lake is much smaller and shallower than Teletskoye Lake so that the sedimentations in the two lakes are different. In the Teletskoye Lake sediments, annual laminations can be found, whereas the Manzherok Lake sediments do not show laminations; they instead contain much more plant remains (OM > 20%). In addition, because Manzherok Lake is small and shallow, it has a longer frozen time than Teletskoye Lake. Hence, the sedimentary hiatus during LIA did not occur in Teletskoye Lake. Under such circumstances, Manzherok Lake is more sensitive to climate change than Teletskoye Lake in terms of geochemical proxies.
In Figure 12, a comparison of the Manzherok Lake record to the TSI record shows good correlations not only on the long-term major changes of the lake but also on decadal events corresponding to the climate changes. In general, a large, deep, and fresh lake with high sedimentation exists under warm and wet climates such as MWP. A small and shallow lake stage with low sedimentation appears under cold and dry climates such as LIA. During a long relatively warm and wet period, a TSI minimum resulting in a dry episode would lead to lake shrinkage that will cause a rise in lake alkalinity and productivity so that increased Sr/Ba and Sr/Rb could result. In Figure 12, numbers 1, 2, and 3 denote the correlations between Sr/Ba (and Sr/Rb) and TSI minimum within age uncertainties. Thus, we conclude that the climate of the study area is strongly influenced by total solar irradiance, with higher TSI resulting in warmer and wetter conditions. Since the westerly and the polar easterly are the two major moisture jet streams to the area, when lower TSI causes a stronger Siberian High, the latter probably pushes both the westerly and the polar front away from the study area, resulting in arid climates. The situation is reversed vice versa.

CONCLUSIONS
The high-resolution multiproxy analyses including pollen, diatom, contents of total organic carbon, carbonate and ash, concentrations of trace metal elements, heavy metal elements, and rare earth elements in an 82-cm sediment core from Manzherok Lake reveal detailed changes in vegetation and climate in the western foothills of the Altai Mountains over the past 1,500 years. A total of 48 AMS 14 C dates combined with 210 Pb dating indicate that the lake had very low sedimentation (only ∼2 cm) during LIA due to cold and dry climates. There are many ABA-treated samples that show old carbon influence on their 14 C ages because ABA treatment cannot remove organic compounds which used old dissolved CO 2 in the lake at high and anoxic stages. High-resolution 14 C dating on such lake cores should apply for solving such a problem.
Vegetation development in this region can be classified as (1) phase of birch forest-steppe during 1,440-1,150 years BP, (2) transitional phase from birch forest-steppe to birch and pine forests dated 1,150-1,070 years BP, (3) phase of birch and pine forests during 1,070-850 years BP, and (4) phase of pine-birch forest-steppe after 850 years BP. Aquatic pollen species such as algae, Nimphae cell, and cianobacterium Anabaena provide details on the changes in the lake ecosystem.
Changes in the contents of organic matter, ash, and CaCO 3 (dominant sediment feature), concentrations of lithophilic and REE elements (detrital inputs), Sr/Ba and Sr/Rb (indicators of lake alkalinity and productivity), and Mo/Mn (anoxic conditions) identify six stages in the lake history corresponding to the climate changes: (1) 1,440-1,150 years BP. The lake was moderately deep with high productivity under relatively warm and wet climates: (2) 1,150-1,070 years BP. Strong surface runoff into the lake resulting from wet climates marked the onset of MWP: (3) 1,070-850 years BP. Large, deep, and fresh lake exited due to warm and wet climates during MWP. Decomposed old organic matter in the anoxic bottom caused old carbon influence on the 14 C ages: (4) 850-500 years BP. In the early stage, the lake dropped its level slowly until about 700 years BP. Then, the lake shrank its size quickly corresponding to the cooling and drying climates: (5) 500-50 years BP. Small and shallow lake with very low sediment deposition corresponded to cold and dry LIA. (6) The lake has been recovered during the current warming century.
The Manzherok record compares well with the total solar irradiance record, indicating that changes in TSI is an important factor to influence climate in the Altai Mountains. The Siberian High became strong during the TSI minima. Consequently, the westerly and the polar front would be pushed away from this region, resulting in arid climates.
The pollen and charcoal records and the metal elemental concentrations indicate that fire-cutting agriculture began to spread in the Altai Mountains several centuries before the Russian colonization, and the metallurgical practice of local human tribes who lived on the western piedmonts of the Altai Mountains took place in XII century AD before the Mongolian invasion.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

AUTHOR CONTRIBUTIONS
TB and H-CL designed Research Topic and wrote the manuscript. H-CL and S-CK performed 14 C dating. TB and VU conducted pollen and geochemical analyses. All authors contributed to the article and approved the submitted version.

ACKNOWLEDGMENTS
We are thankful to T. S. Papina and P. A. Blyakharchuk for assistance in the field drilling of the Manzherok Lake sediments.