Reconstruction of Climate Changes on Based δ18Ocarb on the Northeastern Tibetan Plateau: A 16.1-cal kyr BP Record From Hurleg Lake

Hydroclimate evolution history and changes in the Tibetan Plateau play significant roles in depicting paleoclimate and evaluating climatic conditions in the coming future. However, the interaction of the westerlies and the Asian monsoon complicates our understanding of the mechanism of climate variation over the Tibetan Plateau. In this study, we assessed the paleoclimate of Hurleg Lake, which was previously located in the convergence area of the East Asian monsoon and westerly wind. We first reconstructed the climatic conditions based on fined-grained authigenic carbonate δ18O (δ18Ocarb), plant-derived proxies of C/N, and n-alkane-derived δ13C31. In the Hurleg Lake, δ18Ocarb was controlled by δ18O changes of the lake water and evaporation. The climate evolution since ∼16.1 cal kyr BP can be classified into three stages. The Lateglacial (16.1–11.0 cal kyr BP) was characterized by a warm-wet climate in the beginning, followed by a cold-dry climate since 12.0 cal kyr BP. Typical warm and cold phases occurred during 14.8–12.0 cal kyr BP and 12.0–11.1 cal kyr BP, which may correspond to the Bølling/Allerød (B/A) and Younger Dryas periods, respectively. The early to mid-Holocene was generally characterized by a warm-wet climate; however, notable cold-dry intervals occurred at ∼8.3 cal kyr BP. The Late Holocene (after 4.8 cal kyr BP) displayed a significantly cold-wet climate. Finally, we examined the possible mechanisms responsible for the climate variability in the study area. The results showed that the long-term warm trend in the Lateglacial and colder trend after early Holocene was controlled by insolation. The Asian summer monsoon and the westerlies played a significant role in determining moisture sources during the Lateglacial. The East Asian monsoon contributed greatly to the moisture variation from the early to mid-Holocene, whereas the westerly winds dominated during the late Holocene. Combined, our findings highlight the complex changes in hydroclimate conditions since the last glacial in the Tibetan Plateau and provide crucial implications for comprehending the hydroclimate pattern in the transition zone of westerlies and Asian monsoon.


INTRODUCTION
Since the 21st century, global environmental problems, such as sea level rise, extreme climate events, and ecological environmental damage, have increased in frequency and intensity in response to global warming (Klanderud, 2005;IPCC, 2014). These climate events have adversely impacted human and social development and attracted the attention of researchers across the world (Haug et al., 2003;Mondoni et al., 2012). Scientifically predicting future climate change is necessary for tackling global climatic problems; however, the scale of modern meteorological records is too short to reliably evaluate and simulate future climate change laws. Moreover, past climate can inform future climate variability (Tierney et al., 2020); therefore, determining the history and forcing mechanisms of past climate is necessary for estimating future climate trends (Lunt et al., 2013;Chen et al., 2016;Tierney et al., 2020).
As many large rivers originate from the extensive alpine plateau and its adjacent mountains, the Tibetan Plateau (TP) supports rapidly growing societies with a population of 3 billion people or almost half of the world's population . Due to its role in social development and its susceptibility to climate change, the TP is a hotspot for studying climate variability. The TP is located in a climatic junction where the Asian monsoon and westerlies interact intensely. The region has high elevation (average elevation >4,000 m above sea level) and cold climatic conditions (mean annual air temperature <4°C) (Liu and Chen, 2000;Chen et al., 2016;Luo et al., 2019), and the ecological environment of the TP is more sensitive to climatic variability than most regions on Earth (Liu et al., 2009;Diffenbaugh et al., 2013). In addition, thousands of lakes are distributed across the TP, where paleoclimate information is well preserved. These include the hydroclimate and temperature records derived from lakes Naleng, Nam Co., Qinghai, Hurleg, Qionghai, etc. (Kramer et al., 2010;An et al., 2012;Zhao et al., 2013;Günther et al., 2015;Wang et al., 2020). Analyses of past climate changes have advanced our understanding of the dynamic driving factors of hydroclimate and temperature, which are closely related to the monsoon and insolation, respectively, and are influenced by regional and global climate changes (Peterse et al., 2004(Peterse et al., , 2011Dykoski et al., 2005;Wang et al., 2005Wang et al., , 2008Zhao et al., 2011).
Hurleg Lake lies on the northeastern Tibetan Plateau, at the intersection area, which was previously affected by both the East Asian monsoon and westerly wind (Winkler and Wang, 1993;Zhao et al., 2011). Due to the particularity of its geographical location and its climatic variability and sensitivity, many studies have reconstructed the Holocene climate based on the results of lake cores derived from Hurleg Lake, including paleotemperature, paleomoisture, and paleovegetation (Zhao et al., , 2011Zhao et al., 2010;Zhao et al., 2013;Fan et al., 2014;Rao et al., 2014;He et al., 2016;Song et al., 2020). However, these climatic records from Hurleg Lake mostly focused on the Holocene, while studies on the Lateglacial were limited. Moreover, the records were often conflicting, due to uncertainties in sediment chronologies, the climatic significance of the proxies used, the detailed evolutionary model of climate and environmental change reconstructed by various environmental indexes from the same core, or the same proxies derived from different cores (Chen et al., 2016). In Hurleg Lake, pollen-based reconstruction of paleoclimate patterns displayed a relatively wet climate in the earliest Holocene (11.9-9.5 cal kyr BP), followed by a dry climate during the mid-Holocene (9.5-5.5 cal kyr BP), while a stable and less dry climate occurred after 5.5 cal kyr BP . The climate pattern was contradicted to most hydroclimate patterns on the northeastern Tibetan Plateau, such as the Qinghai Lake (Chen et al., 2016), Ximencuo Lake (Herzschuh et al., 2014), and Gahai Lake . What is more, Zhao et al. (2007) suggested that the position of the subtropical monsoon, the midlatitude westerlies, and the local topography all contributed to the climate pattern in the Hurleg Lake area, while Zhao et al. (2010) found that the generally dry climate between 7.8 and 1.0 cal kyr BP could be correlated with decreases in East Asian and Indian monsoon intensities, and Fan et al. (2014) hold that the Westerlies and local topography rather than Asian summer monsoon dominate the moisture availability in this region during the Holocene. Thus, the climate pattern of the northeastern Tibetan Plateau remains unclear and controversial. These records also suggest the need for a denser network of sites to better clarify the spatially complex pattern of climate change and atmospheric circulation in this region.
Herein, we present a well-dated and multi-proxy paleoclimate record from a sediment core in Hurleg Lake (northeastern TP) covering ∼16.1-∼0.2 cal kyr BP, based on the characteristics of authigenic carbonate oxygen isotopes, organic indexes (C/N), and n-alkane-derived indexes. Combined with multi-proxies including previously published organic and inorganic geochemical indexes, we reconstruct the detailed climate evolution in the study area and discuss the forcing mechanisms. Our results provide more accurate scientific data and references for paleoclimate reconstruction to strengthen predictions of future climate evolution on the Tibetan Plateau.

STUDY AREA
Hurleg Lake (37°14′-37°20′ N, 96°51′-96°57′ E, elevation 2,817 m a.s.l) lies east of the Qaidam Basin, in the northeastern Tibetan Plateau (Figure 1). The lake covers an area of approximately 57.9 km 2 , with an average water depth of 2.0 m and a maximum water depth of 9.6 m. The salinity of Hurleg Lake, which belongs to the freshwater lake of the sodium sulfate subtype, is about 1 g/l . The Hurleg Lake region is outside the East Asian summer monsoon (EASM) edge currently and generally has a continental arid climate (Tian et al., 2001(Tian et al., , 2003. The region is characterized by a large difference in the mean temperature between summer and winter, ranging from −9.7°C to 15.6°C, and the average annual temperature was about 4°C (1.6°C-5.2°C) over the past 48 years , Delingha climate station) (Yi et al., 1992). The annual precipitation is approximately 160 mm, and annual evaporation is ca. 2,000 mm (Yi et al., 1992;Zhao et al., 2013). The water in Hurleg Lake is mainly recharged by surface runoff, from the Bayin River and Barogen River, which both eventually flow into Toson Lake through a small stream (Wang and Dou, 1998;Zhao et al., 2010Zhao et al., , 2011. Located in the temperate desert zone, the main vegetation in the lake area is desert grassland (gravel desert) and desert shrub (vial desert) (Zhou et al., 1990). The coastal area of the lake is surrounded by a large portion of marsh vegetation including common reeds (Phragmites communis), Hippuris spp., and Scirpus spp., .
Six sediment organic matter samples and two plants residue samples derived from the HL18-3 core ( Table 1) were collected for accelerator mass spectrometry (AMS) 14 C dating (Ma et al., 2021). The 14 C age was calibrated with the IntCAL13 curve (Reimer et al., 2013), and then the calibrated 14 C ages were used to reconstruct the age-depth model in the R software package bacon v.2.3.7 (Figure 2) (Blaauw and Christen, 2011;Blaauw et al., 2020). Considering the "old carbon" effect, the upper 3.5-m core was analyzed through the method of geomagnetic relative paleointensity (RPI). Before the RPI test, the rock magnetic and environmental magnetic experiments were conducted. Using NRM30mT/SIRM30mT as the index of RPI change, the upper 167 cm generally conforms to the established criteria for RPI reconstruction (Ma X. et al., 2021). By matching our RPI results with published RPI records, three credible RPI ages were gotten in the upper 167 cm of the HL18-3 core ( Table 1) (Ma et al., 2021). The AMS 14 C ages were determined on the AMS at the Beta Analytic, Miami, USA. The magnetic experiments were carried out in the Palaeomagnetism and Geochronology Laboratory of the Institute of Geology and Geophysics, Chinese Academy of Sciences. Data of AMS 14 C ages as well as the rock magnetic and environmental magnetic were published on the journal "Quaternary Geochronology" (Ma et al., 2021).

Authigenic Carbonates
The analyses of δ 13 C and δ 18 O of carbonates were conducted on the gas stable isotope mass spectrometer (Thermo Scientific Delta V-GasBench II, Thermo Fisher Scientific, Waltham, MA, USA). Prior to measurement, plant residues and microbiological shells were picked out under a magnifying glass. After that, each sample was sieved less than 200 mesh. The treated sediment samples were loaded into glass vials and placed in an aluminum heating block of GasBench II; samples (∼0.1 g) were purged by helium for 7 min, followed by treating with supersaturated phosphoric acid (98%) for 2 h, and the reaction temperature is 72°C. Ultimately the CO 2 generated by the reaction was carried into the Delta V by helium for detection. The international reference material is NBS-18, the national reference material is GBW04405, and the delta values are reported in per mil relative to the Vienna Pee Dee Belemnite (VPDB), with error less than 0.1‰.

Organic Proxies
The organic elements were determined by a vario MICRO cube (Elementar Company, Langenselbold, Germany). After combustion, the sample was further catalyzed and C and N elements were converted into detectable gases within a redox process. Pure oxygen was introduced into the system to support combustion. After combustion, the sample was further catalyzed and redox catalyzed, in which C and N elements were all converted into various detectable gases. After further separation by a chromatographic column, these gases were finally analyzed by the detector.

X-ray Diffraction Method
The X-ray diffraction (XRD) experiment was conducted on Ultima IV. Samples were crushed to <200 mesh in an agate mortar, and then the powder sample was loaded into the aluminum circular sample tank with the sample spoon. Finally, the sample was transferred into the instrument for FIGURE 2 | Age-depth model of the HL18-3 core based on bacon. The red solid circles are the RPI age, the blue marks are AMS 14 C ages dated by organic matter, and the green mark is AMS 14 C age dated by plant remains; the shaded part between the gray dotted lines represents the error interval.

RESULTS
A total of 12 samples were collected at intervals of ∼50 cm for the analysis of the XRD results to determine the mineral types. The main carbonate mineral is calcite, ranging from 17.4% to 86.3%, with an average value of 52% (n 12). The total percentage of calcite and aragonite varied from 74.0% to 92.2%, with an average value of 82.4% (n 12). Other minerals, such as dolomite and gypsum, are rare ( Table 2).
A total of 65 subsamples were collected at 8-cm intervals from the HL18-3 core and analyzed on the Delta V (Thermo Scientific). The values of δ 13 C carb ranged from −2.29% to 5.19‰, with an average value of 1.88‰ (n 65), and those of δ 18 O carb ranged from −8.39% to −1.52‰, averaging −5.16‰ (n 65) ( Table 3). δ 13 C carb displayed a generally positive trend during the Lateglacial and became negatively correlated with fluctuation from Early to Late Holocene ( Figure 5A). C/N ratios were generally utilized to distinguish organic matter derived from various sources. In the HL18-3 core, 64 samples were analyzed for C and N elements (Table 3), and the C/N ratios varied from 29.3 to 0.1, with an average value of 12.7 (n 64). The C/N ratios were relatively high during the early to mid-Holocene, while an uptrend occurred in the Lateglacial and displayed a decreased tendency in the late Holocene interval ( Figure 5A).

Age-Depth Model
Using linear extrapolation, the calculated RPI age at 9 cm of the core was 160 cal yr BP, and that at 50 cm of the core was 1,000 cal yr BP; the corresponding AMS 14 C ages were 2,595 cal yr BP and 3,773 cal yr BP, respectively. Thus, the "old carbon" reservoir ages at 9 and 50 cm were 2,435 cal yr BP and 2,773 cal yr BP, respectively, which is close to the reservoir age (2,758 cal yr BP) of core HL06-2 from Hurleg Lake (Zhao et al., 2009;Zhao et al., 2010).
However, owing to the spatiotemporal variability of the carbon reservoir, the carbon reservoir could differ between the sediment  cores from the same area (Reimer and Reimer, 2006;Mischke et al., 2013;An et al., 2012;Zhou et al., 2014). Moreover, the carbon reservoirs within the same core are not always constant (Schneider et al., 2019;Soulet et al., 2019), and the commonly used method of subtracting the carbon reservoir ages from the whole section may be unreasonable. Nevertheless, the occurrence of a linear relationship or a quadratic regression suggests that the carbon reservoirs are constant because the hydrochemical conditions in the lake are generally stable (Xiao et al., 2004;Herzschuh et al., 2006;An et al., 2012;Zhang et al., 2016). According to the analysis of the relationship between ages (including three RPI ages in the upper 110 cm and five AMS 14 C ages under 110 cm) and depths, we obtained a significant linear relationship (R 2 0.93, p < 0.001, n 8) and quadratic regression (R 2 0.99, p < 0.001, n 8) (Figure 3). Three ages obtained by the geomagnetic RPI method are reliable and not affected by the reservoirs. The linear and quadratic regression relationship suggests that there is little or no carbon reservoir effect below 110 cm (Ma et al., 2021). To support this finding, the radiocarbon dating of plant remains in a sample at 380 cm returned a date practically identical to that dated by the organic matter at the same depth ( Table 1).

Origin of Carbonates
Previous studies have shown that the carbonate in lake sediments with grain sizes less than 200-300 mesh is authigenic carbonate (Jiang and Wu, 2003;Zhang et al., 2016). In order to minimize the possible effect of exogenous carbonate, samples <200 mesh were selected for testing carbonate isotopes. In addition, the XRD results displayed that calcite was the dominated carbonate mineral, with the calcite and aragonite contents accounting for 82.3% of the total carbonate. In contract, the contents of dolomite and gypsum are relatively low (Table 2). Moreover, we found no large exposed carbonatite in the study area. Consequently, finegrained carbonates can represent lake authigenic carbonate.
Furthermore, several studies have reported positive or negative correlations between δ 18 O carb and δ 13 C carb . Talbot and Livingstone (1989) and Talbot (1990) found that δ 18 O carb was positively correlated with δ 13 C carb across a large number of lakes in central Asian and Africa; they also identified a stronger correlation with longer water retention time but no correlation in open lakes. In addition, Leng and Marshall (2004) revealed that δ 18 O carb correlated positively with δ 13 C carb in the sediment of freshwater and saltwater lakes, but not in groundwater recharge lakes. The positive correlation between δ 13 C carb and δ 18 O carb in the lakes of northwestern China reported by Zhang et al. (2016) indicated a closed lake, with a long retention time and obvious salinization characteristics. The δ 18 O carb and δ 13 C carb results of carbonate in the HL18-3 core showed consistent variability patterns, exhibiting a positive correlation (Figure 4, n 65, R 2 0.4071, p < 0.001). The Hurleg Lake is a freshwater lake in which the sodium sulfate subtype is about 1 g/l . Although Hurleg Lake is currently an open freshwater lake, its salinity was generally higher in the past ( Figure 6E) (He et al., 2016). In addition, according to the annual total river inflow and lake volume, the mean residence time of the Hurleg Lake is approximately 1 year .
The variability in carbonate content under different water depths at Hurleg Lake was found to be controlled by the dilution of endogenic carbonate by fine-grained silicates delivered from surrounding deserts with sparse vegetation cover . Thus, the change in the carbonate percentages was used to reflect the variation in water depth. Prior to the study, Zhao et al. (2010) observed a flash flood event in 2005, in which a large amount of alluvial sediment was transported from the area west of the lake. However, more detrital material was deposited in shallow water because the dense vegetation covering the lake bottom trapped the materials in the littoral zone. Similarly, the inflowing rivers are surrounded by dense marshes, and river channels are unlikely to transport clastics into the lake . The HL18-3 core is located in the southwestern part of the lake, where it is far away from the river mouth ( Figure 1B). Although exogenous clastic minerals might have an effect on sedimentary carbonate, the effect is likely to be slight. Accordingly, we suggest that the fine-grained carbonate in Hurleg Lake is mainly derived from authigenic carbonates, and the variabilities of carbonate δ 18 O carb are capable of indicating the climate change in the study area.

Implications of Climatic Proxies
Climatic Implications of δ 18 O carb The carbonate δ 18 O carb values in lake sediments are generally controlled by the δ 18 O of lake water, water temperature, and disequilibrium effects (Bowen, 1990;Gasse et al., 1991;Yu and Eicher, 1998;Leng and Marshall, 2004). Tian et al. (2003) analyzed the δ 18 O of modern precipitation water derived from the Delingha, where it is nearby the Hurleg Lake; the results revealed a strong temperature influence on precipitation isotopic composition in the Qaidam Basin, with higher temperature corresponding to higher isotopic value in precipitation. However, water temperature affects the δ 18 O fractionation from lake water into carbonate, which could induce an equilibrium fractionation coefficient of −0.25‰ per°C during calcification (Epstein et al., 1953). The changes in δ 18 O values of the HL18-3 core sediments reached 6.87‰, which led to an approximately 27.5°C temperature difference, and there is no evidence for such a high-amplitude change. Therefore, the variabilities of δ 18 O carb cannot account for the variation in the water temperature. The disequilibrium effects were proven to be negligible at Gengahai Lake on the northeastern Tibetan Plateau (Qiang et al., 2017) and were rarely discussed when evaluating the paleoclimate. Therefore, the δ 18 O carb records of the Hurleg Lake could have been controlled primarily by variations in the δ 18 O of lake water (δ 18 O Lw ). There is another fact that largely parallel first-order variations in the δ 18 O values derived from various carbonate species suggest that all the δ 18 O records provide a reasonable approximation of lake water δ 18 O (Apolinarska and Hammarlund, 2009;Clegg and Hu, 2010;Qiang et al., 2017). In Hurleg Lake, water samples have an average δ 18 O value of −4.9‰ (n 15) (He et al., 2016); besides, water isotope measurements of 13 samples taken from Hurleg Lake between 2004 and 2007 have a range from −6.6% to −3.3‰ for δ 18 O . The average value of δ 18 O carb in the HL18-3 core is −5.16‰ (n 65) ( Table 3), and the average value of δ 18 O carb from the upper 50 cm of the HL18-3 core is −6.20‰ (n 6); the δ 18 O carb values are similar to the δ 18 O in the lake water. Therefore, it is considered that the changes in δ 18 O carb can represent the variation of δ 18 O in lake water. Hence, factors contributing to δ 18 O Lw will in turn affect the composition of δ 18 O carb . Among these factors, the δ 18 O values of the regional rainfall and lake hydrology are the most significant ones (Zhang et al., 2011). Meteorological data from the Delingha weather station show that the lake area of Hurleg Lake was significantly correlated with the annual runoff of the Bayin River (R 0.78), while the latter is closely related to the mean annual precipitation (Fu et al., 2008;Fan et al., 2014). On the other hand, compared to the annual runoff of the Bayin River, meteorological analyses show that melting ice and snow from glaciers are in low proportion (Fu et al., 2008;Fan et al., 2014). Therefore, the lake area and lake level of Hurleg Lake are dominated by the annual precipitation of the lake catchment (Fan et al., 2014). Correspondingly, the variations in moisture sources and precipitation amounts were considered as significant controls on δ 18 O carb .
In general, moisture transported by the Westerly wind to the region is much more 18 O-depleted than that from the summer monsoon (Araguás-Araguás et al., 1998;Tian et al., 2007；;Zhang et al., 2011). Hurleg Lake is located in the transition zone of Westerly wind and Asian monsoon in the past (Winkler and Wang, 1993;Zhao et al., 2011); the variation of moisture source has excellent potential for affecting δ 18 O in precipitation. Besides, there is a negative correlation arising from the "amount effect," whereby the heaviest rainfall is most depleted in δ 18 O (Dansgaard, 1964;Zhang et al., 2011). In the monsoon region, significant correlations between δ 18 O of precipitation and rainfall amount were observed from the meteorological station (Johnson and Ingram, 2004;Dayem et al., 2010). Although the precipitation amount in the Hurleg Lake region now is low, in the past, the precipitation amounts probably are in large value. This is due to the elevation of three palaeoshorelines in the Holocene which are ∼4, ∼7, and ∼9 m above the present lake level (Fan et al., 2014), and there was a mega-paleolake during the early Pleistocene to late Pleistocene (∼2,800-2,700 m a.l.s) in the north-eastern Qaidam Basin (Huang and Cai, 1987). Thereby, the rainfall amounts could also contribute to the variabilities of δ 18 O in precipitation, in turn affecting δ 18 O carb in lake sediments.
In dryland regions, for closed lakes, the water loss is dominated by the evaporation; δ 18 O Lw is therefore strongly influenced by the precipitation/evaporation (P/E) ratio. In previous studies, variations in the δ 18 O carb of lake sediments derived from arid to semiarid regions have always been interpreted as reflecting the balance between regional precipitation and evaporation (P/E) (Lister et al., 1991;Wei and Gasse, 1999;Liu et al., 2007;Zhang et al., 2011). In south Tibet, Tian et al. (2008) suggested that the δ 18 O Lw values in several closed lakes are very sensitive to relative humidity. However, the Hurleg Lake is an open lake, the rapid discharge of a lake system always reduces the lake water residence time and thus decreases the effect of evaporation on δ 18 O Lw (Develle et al., 2010). The Donggi Cona Lake in the Tibetan Plateau provided good evidence, due to the δ 18 O carb values becoming more negative after ∼4.0 cal kyr BP, when the lake became hydrologically open . Nevertheless, a recent study on isotopes of Hurleg Lake water displayed that the δ 18 O values suffered from a cumulative evaporation effect in which it became stronger from source water to hydrologically open lakes (e.g., Hurleg Lake), then to closed lake systems (He et al., 2016). Further, the results of 2 H and 18 O derived from Hurleg Lake water samples fall on a local evaporation line, indicating that the lake water isotopes were influenced by evaporation (He et al., 2016). Therefore, it seems that the impact of evaporation on open lakes cannot be ignored.
Overall, the shifts of δ 18 O carb values from negative to more positive in the Hurleg Lake region over the episodes since the Lateglacial can be interpreted primarily as a decrease in the amounts of precipitation coupled with an increase in the regional P/E ratios, indicating the decrease in regional effective precipitation (Zhang et al., 2011) and representing a dryer climate. Besides, the changes in the moisture source between Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 745972 8 westerly wind and Asian monsoon may also have played a dramatic role in controlling the δ 18 O carb values.

Climatic Implications of Organic Proxies
Organic matter in the lake sediments originating from terrestrial higher plants contains a large number of nitrogen-free biomacromolecules, such as lignin, cellulose, and hemicellulose, which are characterized by carbon enrichment. In contrast, aquatic plants, such as phytoplankton, are characterized by nitrogen enrichment because they are rich in proteins and other substances (Meyers, 1994b;Meyers, 1997;Hedges, et al., 1997). Therefore, the C/N ratio can be utilized to distinguish organic matter derived from various sources, where the C/N ratio ranges of 4-10, 10-20, and >20 correspond to aquatic plants or algae, submerged/phytoplankton or mixed sources of terrestrial and emergent plants, and terrestrial higher plants, respectively (Meyers, 1994b;Jia and Peng, 2003;Hu et al., 2009).
Atmospheric CO 2 exists in three isotopic forms ( 12 CO 2 , 13 CO 2 , and 14 CO 2 ); owing to the different activities of light and heavy isotopes in thermal gradients or biochemical reactions, carbon isotope fractionation occurs in plants during transpiration and biosynthesis. The mean δ 13 C TOC value for C 3 plants is −27‰ (range from −33‰ to −22‰), and for C 4 it is −13‰ (range, −16-9‰) (Cerling et al., 1997). For individual leaf wax n-alkanes, δ 13 C values are ∼6‰-8‰ lower than for bulk tissues (Rieley et al., 1993;Huang et al., 2001); thus, the variabilities of individual leaf wax n-alkanes can reflect changes in C 3 and C 4 plants. C 4 plants prefer to grow in areas with sufficient sunshine, higher temperatures, and drier conditions, primarily composed of warm and dry grassland vegetation (Edwards and Smith, 2010;Strömberg, 2011). C 4 plants have a higher water-use efficiency than C 3 plants and enriched in δ 13 C 31 (Edwards and Smith, 2010;Strömberg, 2011). Several recent studies suggested that submerged aquatic plants with enriched δ 13 C may contribute to the concentrations of long-chain n-alkanes (nC 27 , nC 29 ), while nC 31 is a reliable proxy for estimating the variabilities of C 4 plants (Liu and Yang, 2015;Liu and Liu, 2016). Thus, δ 13 C 31 was used to reconstruct the evolutionary history of C 4 vegetation (Ma et al., submitted for publication). Generally, higher δ 13 C 31 suggests an expansion of C 4 plants, corresponding to a warmer climate.

Climate Changes on the Northeastern Tibetan Plateau Since the Lateglacial
Lateglacial (Ca. 16.1-11.1 cal kyr BP) The generally enriched δ 18 O carb value of fine-grained carbonates suggested a decreased effective precipitation and thus a dryer climate during this period ( Figure 5A). However, during this period, the fluctuating upward trend of mean annual precipitation in the Ximencuo Lake northeastern TP coincided with a wetter climatic condition ( Figure 7E) (Herzschuh et al., 2014), and the depleted δ 18 O of stalagmite calcite from Dongge Cave indicates a strong intensity of Asian monsoon (Dykoski et al., 2005). In addition, the increased C/N ratio revealed that the contribution of terrestrial higher plants to sediments increased ( Figure 5B), and the C 4 plants expanded as indicated by more positive δ 13 C 31 values ( Figure 5C). These results indicated a flourishing of terrestrial and aquatic plants, corresponding to a warm-wet climatic environment. Typically, moisture transported by the summer monsoon to the Tibetan Plateau is more δ 18 O-enriched than that of the westerly wind, because of the higher condensation temperature and/or weaker "rain-out" after a shorter air-mass transport distance (Araguás-Araguás et al., 1998;Tian et al., 2007;Qiang et al., 2017). In addition, the δ 18 O values of fresh snow samples derived from the Qaidam Basin in January 2000 ranged from −29.4‰ to 20.8‰ (average of −26.2‰), which is significantly more negative than that of Hurleg Lake water (Qiang, 2002;Zhao et al., 2010). Furthermore, a previous study on the northeastern Tibetan Plateau revealed that the hydroclimate was controlled by monsoon-derived, westerly derived, and/or locally recycled moisture in several durations of the Lateglacial (Qiang et al., 2017). Henderson et al. (2010) attributed the depleted δ 18 O carb during the Little Ice Age at Lake Qinghai to the enhanced input of westerly wind. Considering the rapid discharge of the Hurleg Lake system, the δ 18 O of lake water was significantly correlated with the annual runoff of Bayin River (R 0.78), (Fu et al., 2008;Fan et al., 2014). The dramatic hydro-variation of the δ 18 O signal induced by the transformation of westerly wind and monsoon would be recorded by the δ 18 O carb at Hurleg Lake. Therefore, the positive shift in δ 18 O carb during 16.1-14.8 cal kyr BP in the Hurleg Lake may be dominated by the strengthened monsoon, and after ∼14.8 cal kyr BP, δ 18 O carb represents multiple changes in monsoon intensity ( Figure 5A). The δ 18 O carb derived from the lake sediments in the adjacent areas of Geanggahai Lake in the Qaidam Basin, northeastern Tibetan Plateau, displayed a similar tendency with that of Hurleg Lake ( Figure 7B). Combined with other evidence derived from Geanggahai Lake, the data suggest that the Asian summer monsoon may start to dominate the study region at ∼15.3 cal kyr BP (Qiang et al., 2017).
Accordingly, the depleted δ 18 O carb indicated a wet climate from ∼14.8 to ∼12.0 cal kyr BP and enriched δ 18 O carb values during 12.0-11.1 cal kyr BP ( Figure 5A), indicating a switch to a drier climate. These findings are in good agreement with the climate reconstructions interpreted by C/N and δ 13 C 31 ( Figure 5B). The wet-warm conditions at 14.8-12.0 cal kyr BP were likely a response to the wet-warm Bølling/Allerød (B/A) event in the northeast Atlantic Ocean (Bard et al., 2000). At that time, pollen records from Lake Ximencuo on the northeastern Tibetan Plateau inferred that a Cyperaceae-rich high alpine meadow replaced the dry glacial flora at ∼15.5 cal kyr BP, which was then further replaced by alpine meadow and shrubland at 14.0-12.5 cal kyr BP, indicating a humid climate (Herzschuh et al., 2014). The corresponding pollen-based mean annual temperature and precipitation showed a gradual increasing trend in 15.5-12.5 cal kyr BP and showed lower values during 12.5-10.8 cal kyr BP ( Figures 7D,E) (Herzschuh et al., 2014). In addition, the U K 37 record at Lake Qinghai inferred a significantly increase in temperature since 15.1 cal kyr BP, which remained at relatively high temperature until 12.6 cal kyr BP. This Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 745972 was followed by a lower temperature during 12.6-11.7 cal kyr BP ( Figure 7C) (Hou et al., 2016). After that, the 12.0-11.1-cal kyr BP dry interval at the Hurleg Lake may be linked to the Younger Dryas period (Shakun and Carlson, 2010), which markedly affected most regions of the Northern Hemisphere, including the Tibetan Plateau. Lake Dongxiong Co. in the southwestern Tibetan Plateau experienced a 3.8°C drop in temperature during 12.3-11.4 cal kyr BP (Ling et al., 2017). Other lake records on the Tibetan Plateau, such as Qinghai Lake (Ji et al., 2005), Lake Sumxi Co. (Gasse et al., 1991), and Lake Nam Co. (Günther et al., 2015), revealed a commonly cold-dry/wet climate lasting from ∼12.9 to ∼10.0 cal kyr BP.
Early to Mid-Holocene (Ca. 11.1-4.8 cal kyr BP) Soon after 11.1 cal kyr BP, the δ 18 O carb values rapidly decreased and reached a relatively low value at ∼10 cal kyr BP and then increased ( Figure 5A), inferring a wet climate in the early Holocene, which may be induced by the enhanced Asian monsoon ( Figure 8E) (Dykoski et al., 2005). Since ∼8.3 cal kyr BP, δ 18 O carb quickly depleted and then gradually enriched, reflecting a wetter climate than the previous stage, and a drier climate in the late mid-Holocene ( Figure 6A). During this interval, the δ 13 C 31 values were generally high in 11.1-9.2 cal kyr BP and then decreased to low values during 9.2-8.3 cal kyr BP. This corresponds to a warm climate at the beginning of the Holocene, which was interrupted by a cold interval around 8.3 cal kyr BP ( Figure 5C). The ∼8.3-cal kyr BP cold event has been commonly acknowledged in large numbers of studies on the Tibetan Plateau (Herzschuh et al., 2006(Herzschuh et al., , 2014Mischke and Zhang, 2010;Zhang and Mischke, 2009) and in the Hurleg Lake Zhao et al., 2010;Zhao et al., 2013;He et al., 2016). The precipitated carbonate δ 18 O of the HL05-2 core showed a similar variation trend to the δ 18 O carb in the HL18-3 core during 9.2-8.3 cal kyr BP ( Figure 6B), which signified a dry climate . Moreover, the lake level reflected by the carbonate content suggested that of the lake contracted notably at ∼7.9 cal kyr BP ( Figure 6C) . The reconstructed salinity via alkenone-based %C 37:4 exhibited relatively high values at 7.7-8.2 cal kyr BP ( Figure 6E), indicating a dry climate (He et al., 2016). Besides, U k' 37 -derived ΔT fluctuated rapidly at low values during this period, which then increased distinctly after ∼8.0 cal kyr BP ( Figure 6D) . Furthermore, pollen concentrations were low during intervals ( Figure 6F) Zhao et al., 2013), and the C 4 plant indicator was significantly reduced ( Figure 5C, Ma et al., submitted for publication); correspondingly, the organic matter percentage declined during 8.3-8.7 cal kyr BP ( Figure 6G) . When considering chronological errors, we determined the occurrence of a vital cold-dry event, which may be related to the frequently recorded 8.2-cal kyr BP cold event in the North Atlantic and many parts of the world (Johnsen et al., 1992;Alley et al., 1997;Thomas et al., 2007;Hede et al., 2010).
The enriched δ 13 C 31 value after ∼8.3 cal kyr BP indicates a recovery of C 4 plants ( Figure 5C). From ∼8.0 to ∼6.2 cal kyr BP, the low δ 18 O carb value suggested a generally wet climate, corresponding to a high lake level, reported by Fan et al. inferred temperature record at Lake Qinghai for the past 16 ka (Hou et al., 2016). (D) Pollen-based T ann (mean annual temperature) from Lake Ximen Cuo on the northeastern TP (Herzschuh et al., 2014). (E) Pollen-based P ann (mean annual precipitation) from Lake Ximen Cuo on the northeastern TP (Herzschuh et al., 2014). (F) Arctic sea ice extent (de Vernal et al., 2013).
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 745972 12 (2014), The excessively humid climate did not favor C 4 plant growth; thus, the δ 13 C 31 value was depleted, representing a decrease in C 4 plant biomass ( Figures 5A,C). Afterward, the higher δ 18 O carb values during 6.2-4.8 cal kyr BP suggest a gradual drying climate, which is beneficial to the growth and expansion of C 4 plants, as reflected by the enriched δ 13 C 31 value ( Figures  5A,C). The changes in precipitated carbonate δ 18 O of the HL05-2 core ( Figure 6B)  and Genggahai Lake ( Figure 7A) (Qiang et al., 2017) were consistent with the δ 18 O carb of the Hl18-3core during this period ( Figure 6B) .
Late Holocene (Ca. 4.8-0.2 cal kyr BP) In the late Holocene, the δ 18 O carb value showed a general depleted trend, similar to that of other climate indices, such as C/N and δ 13 C 31 (Figures 5A-C). The δ 18 O carb value was enriched during 4.8-4.3 cal kyr BP, while δ 13 C 31 decreased rapidly, indicating a contraction of C 4 plants. All of these proxies signified a strong abrupt cold-dry event ( Figures  5A,C). The temperature indicator U k' 37 -derived ΔT showed low values during 4.9-4.3 cal kyr BP. Moreover, the carbonate content declined revealing a low lake level, and a higher salinity appeared at ∼4.6 cal kyr BP, suggesting a dry climate ( Figures 6C-E) Zhao et al., 2010;He et al., 2016). The pollen concentration and organic matter percentage were in low value from ∼4.2 to 3.7 cal kyr BP ( Figures 6F,G) Zhao et al., 2013), which demonstrated a dry and deteriorated climate. These results generally signify a cold-dry event, correlated with a 4.2cal kyr BP event (Bond et al., 1997;Mischke and Zhang, 2010;Hou et al., 2016). The positive δ 18 O carb values observed in the Genggahai Lake during ∼5.0-∼3.5 cal kyr BP, as well as the decreased precipitation in Ximencuo Lake on northeastern Tibetan Plateau during 4.9-3.8 cal kyr BP, indicated that a hydroclimate was dry, while the temperature in the adjacent area of Hurleg Lake, such as the Qinghai and Ximencuo lakes, decreased dramatically ( Figures 6C,D) (Herzschuh et al., 2014;Hou et al., 2016). In addition, the sea ice coverage in the Barents Sea increased significantly ( Figure 7F) (De Vernal et al., 2013), sea surface temperature decreased (Lea et al., 2003), and Greenland ice core δ 18 O (Stuiver et al., 1995) notably decreased during this period, indicating a global cold event at ∼4.2 cal kyr BP.
After the rapid cold-dry event, the climate conditions improved slightly around 4.2-3.0 cal kyr BP, as shown by a negative δ 18 O carb (Figure 5A), indicating a higher input of effective precipitation. However, the C/N ratio revealed a higher input of terrestrial plants, while the depleted δ 13 C 31 pointed to less C 4 plants ( Figure 5C). Here we argued that although the C 4 plants decreased, with the more rapid expansion in C 3 plants, it led to a high value of the C/N ratio. In the meantime, the decrease of C 4 plants also suggests that the climate is gradually turning cold ( Figure 5C). Subsequently, the climate displayed a generally cold-wet trend, as revealed by δ 18 O carb , C/N, and δ 13 C 31 (Figures 5A-C). Fan et al. (2014) reported that a high lake level occurred during 2.2-1.4 cal kyr, and Song et al. (2020) also revealed a moderately wet climate after 2.5 cal kyr BP by sediment grain size and δ 13 C of authigenic carbonate in the Hurleg Lake.

Implications Between EASM-Westerly Interaction and Climate Changes
Multiple climate systems including the EASM, the Indian Summer Monsoon (ISM), and the westerly jet An et al., 2012;Hou et al., 2017) together contributed to the complex climate condition of the Tibetan Plateau. The Hurleg Lake region now is located outside of the northern edge of the EASM ( Figure 1A), deduced from the stable isotope analysis of local precipitation collected nearby the Delingha station located ∼30 km northeast of the lake (Tian et al., 2001(Tian et al., , 2003. The region is located in the transition zone between the westerlies and the EASM, and owing to the instability of the transition zone between the westerly jet and the EASM, the Hurleg Lake area was influenced by both the westerlies and the EASM in the past (Bryson, 1986;Winkler and Wang, 1993;Zhao et al., 2007). However, in contrast to the acknowledged synchronous variation between precipitation and temperature (Zhang et al., 2011) controlled by the EASM, regimes dominated by the westerly jets are multivariable in moisture and other climate conditions (Chen et al., 2019). Neither our temperature indicator (δ 13 C 31 ) nor our moisture index (δ 18 O carb ) showed a consistent pattern with EASM intensity ( Figure 8E) or the westerlies; thus, we assume that both of them affected the study area.
In general, the n-alkane-derived δ 13 C 31 values ( Figure 8A, green line)-the temperature-related proxy-showed a similar trend to that of maximum insolation at 35°N since ∼16.1 cal kyr BP ( Figure 8A, black line), which suggests that the long-term temperature variability was controlled by insolation. However, the variation in humidity in the northeastern TP is much more complex than that of temperature. At ∼14.8 cal kyr BP, the westerly index (WI) recorded in Qinghai Lake decreased significantly ( Figure 8C) , while the Asian monsoon intensified gradually ( Figure 8E) (Dykoski et al., 2005). This result coincides with our δ 18 O carb values, which suggests that the enriched δ 18 O carb at ∼14.8 cal kyr BP was caused by the enhanced Asian monsoon. Although we attributed the δ 18 O carb variabilities in the HL18-3 core to changes in the Asian monsoon precipitation during 14.8-12.0 cal kyr BP, there may be other controlling factors. We noticed another turning point at ∼13.3 cal kyr BP, in which the westerly wind intensified and the Asian monsoon weakened during 13.3-12.0 cal kyr BP. The opposite change occurred during 12.0-11.0 cal kyr BP ( Figures 8C,E) Dykoski et al., 2005). The negative δ 18 O carb value indicates that a wetter climate could have been caused by both the Asian monsoon and westerly wind during 13.3-12.0 cal kyr BP. According to the reconstructed δ 18 O and δ 2 H values of precipitation in the Qinghai, Bangong Co., and Linggo Co. Lakes in the northern TP, the positive δ 2 H values in the three lakes during 10.0-12.0 cal kyr BP support a weakened monsoon in the TP during the late Pleistocene (Yao et al., 2013;Hou et al., 2017). The abrupt decrease in δ 2 H in all three lakes during the early Holocene signified an increased monsoon intensity (Yao, Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 745972 1997; Hou et al., 2017). Thus, it is considered that since ∼11 cal kyr BP, the variation of moisture in the Hurleg Lake area was predominantly controlled by Asian monsoon, until ∼4.0 cal kyr BP, in which the intensity of monsoon declined gradually ( Figure 8E) (Dykoski et al., 2005). Nevertheless, the intensity of the monsoon has been decreasing over the past ∼4.0 cal kyr BP (Dykoski et al., 2005), while the depleted δ 18 O carb value reflects a wetter climate condition ( Figures 8B,E). Thus, the Asian monsoon is unlikely to be the dominant moisture source after ∼4.0 cal kyr BP, while the synthetic model of humidity from the entire area of westerly Asia (ACA) displayed an increasing trend ( Figure 8D) (Chen et al., 2019), which likely contributed to the increased humidity in the Hurleg Lake region. Thomas et al. (2016) suggested that the hydroclimate of Qinghai Lake during the late Quaternary was mainly dominated by westerlyderived moisture, and the low δ 18 O carb values in the lake during the Little Ice Age was likely caused by westerlysourced moisture (Henderson et al., 2010). Therefore, the negative shift in the δ 18 O carb value during the late Holocene could have resulted from strengthened westerlies, which provided more depleted δ 18 O precipitation to the lake (Qiang et al., 2017). The pattern of oxygen isotope fluctuations at Hurleg Lake strongly suggests that the northeastern TP was alternately dominated by the westerlies in the late Holocene.
In summary, the hydroclimatic evolution of Hurleg Lake, as revealed by multiple proxies, including δ 18 O carb , C/N of HL18-3 core, carbonate content and organic matter percentage , and pollen concentration , was generally controlled by both the westerlies and the Asian monsoon during the Lateglacial, the Asian monsoon during the early to mid-Holocene, and the westerlies during the late Holocene.
2) 11.1-4.8 cal kyr BP (early to mid-Holocene): warm-wet climate conditions occurred in the early Holocene, which was interrupted by the 8.3 cal kyr BP cold-dry interval. After then, the climate became wet and warm during the mid-Holocene. 3) 4.8-0.2 cal kyr BP (late Holocene): a cold-dry event occurred during 4.8-4.3 cal kyr BP, which was followed by a cooler climate with increasing humidity.
Comparisons of regional climate records suggest that although several global cold events affected the climate condition of the northeastern TP in the short term, the long-term alternation of cold-warm environmental conditions was controlled by insolation. The hydroclimatic evolution reflected by δ 18 O carb documented both westerly and Asian monsoon effective precipitation. By combining multiple hydroclimate proxies, we found that the precipitation in the study region was firstly controlled by both westerly wind and Asian monsoon during the Lateglacial, the Asian monsoon during the early to mid-Holocene, and the westerly winds during the late Holocene.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.