Climate change inferred from aeolian sediments in a lake shore environment in the central Tibetan Plateau during recent centuries

BenLi Liu , JianJun Qu , ShiChang Kang , Bing Liu

1. Dunhuang Gobi Desert Ecological and Engineering Research Station/Key Laboratory of Desert and Desertification,Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, Gansu 730000, China

2. State Key Laboratory of Cryospheric Sciences, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, Gansu 730000, China

1 Introduction

The Tibetan Plateau (TP) is the highest plateau on Earth (reaching a height of over 4,200 m and covering about 2.5×106km2), and strongly influences the climate environment of nearby regions, East Asia and the whole globe as well (Wuet al., 2015). TP is one of the most sensitive regions in respect to global warming, as there is an increase in warming with altitude (Holmeset al., 2009). Intensive land degradation, desertification and aeolian activity since the Little Ice Age (LIA, during about 1350–1850 A.D.)are widely observed over the TP (Jinet al., 2007; Yuet al., 2017). The effects of climate change, the status of land degradation, permafrost warming, glacier shrinkage, desertification, among others, on the TP has received increased attention by researchers and society (Cheng and Wu, 2007; Xueet al., 2009; Yaoet al., 2012).

Some important climate reconstruction projects over the TP were carried out, in which a limited number of proxy data were used to represent this vast area(Yanget al., 2003). Overall warming conditions over the TP has occurred over the past few centuries as found in several lakes over the TP, and the warming rate increases with elevation, with the highest rate between 4,800 and 6,200 m (including the mountain area in this study) (Liu and Chen, 2000; Qinet al.,2009). However, a complete understanding of past environmental and climate conditions on the TP is still lacking. The density of data points over the TP is currently too low for identifying spatial patterns of any changes, because previous studies were performed at various locations using different climate proxies and of different temporal resolution (Br?uning, 2006;Holmeset al., 2009). Large deviation may arise given the strong variation of altitude and local climate over the Plateau. Data with high-confidence and high temporal interval in the central TP area is generally harder to obtain compared to other areas because most studies were concentrated on the northeastern and southern TP (Stauch, 2015). A new ice-core record was used to reconstruct the temperature variation for Geladaindong Mountain during the LIA and 20th-century (Zhanget al., 2016), but studies in the central TP area are still under-represented compared to other TP regions and needs more reinforcement.

Sediment deposition contains traces of past climate, environmental evolution, and responses to global change (Mügleret al., 2010). Sedimentary information has been used in some studies to evaluate longperiod paleoclimate at different lakes including Bangong and Sumxi in the west (Gasseet al., 1991,1996), Nam, Siling, Cuoe and Ahung in the central area (Morinagaet al., 1993; Mügleret al., 2010; Morrillet al., 2006; Wuet al., 2006), Kuhai and Qinghai in the northeast (Shenet al., 2005; Mischkeet al.,2010; Anet al., 2012; Luet al., 2015), Xinmocun in the east (Liang and Jiang, 2017), and Sugan in the Qaidam Basin of northern TP (Chenet al., 2009).Also, aeolian deposits have been recognized as providing information about past environments (Liuet al., 2011, 2013, 2015; Ijmkeret al., 2012; Yu and Lai, 2014; Qianget al., 2016) (Figure 1a). However,sediment reconstruction studies are considered highly uncertain and of low resolution in chronology(Holmeset al., 2009).

Recently, active wind erosion and desertification to the east side of Lake Co Nag in the central TP area has received increased attention. This is due to increased aeolian activity, resulting in significant sand accumulation which could obstruct the Qinghai-Tibetan Railway (QTR) (Zhanget al., 2010). The sands could be mobilized under the effects of freeze-thaw and wind erosion (Xieet al., 2016). Recorded climate data and model projections suggest that temperature and precipitation in the area increased over the past six decades and the trend may continue in the near future (Liuet al., 2016a). Unfortunately, research is lacking on past climate before the recording of meteorological data in the 1950s in the region of Lake Co Nag. Evaluating historic climate conditions over a longer period will help in understanding the cause and driving force of the rapid change of land surface conditions in the fragile environment on the TP.

In this study, we examined an aeolian sand sequence at the shore of Lake Co Nag with optically stimulated luminescence (OSL) age dating, grain-size and geochemical analysis. Grain-size and geochemical parameters, and depositional rate of sand sediments were used to understand climate change in the central TP. Our purpose is to reconstruct the past climate and to reveal the response of land surface process (aeolian activity) to climatic change, as well as to understand the future development of desertification and aeolian sand deposition problems. We demonstrate that in this region that lacks a rich climate proxy,aeolian sediments can be used to study the climatic history and geomorphic processes of central TP.

2 Materials and methods

2.1 Study area and sampling

The study area is located near to the QTR and at the east shore of Lake Co Nag in Anduo County,Tibet (Figure 1a). Lake Co Nag is the head of the Nujiang River (also called the Salween or Thanlwin River), which runs through Southwest China to Burma and Thailand. The area is characterized by high wind and pervasive wind-blown sand activity.The average wind speed (over multiple years) is 4.1 m/s, and there are 150 days per year with high wind speeds (＞17.2 m/s) (Liuet al., 2016a). Peak wind speed is 38 m/s in this area which is significantly higher than any other regions along the QTR(Zhanget al., 2010).

To the east of QTR, a sandy area that occupies approximately 180 km2covers the hills that were carved by extensive seasonal gullies and streams flowing into the lake (Figure 1b). The land surface is covered by short grass in the summer, but the sandy layers are continually being stripped away by aeolian process,leaving wind erosive pits and isolated sedimentary piles that still has grass growing on top. As a result,land degradation is far worse than can be seen from satellite images. The dislocated materials can be found along gully banks and in the sand control sys-tems at the railway sides in the downhill direction towards the lake, clearly indicating a transportation of deposition from the mountain to the lake. However,statistics (sand drift potential) on the energy of surface winds in terms of sand movement shows that the westerly wind from the lake to the mountain dominates. This means a strong interplay of sand transportation from the mountain and the lake, that is, the prevailing westerly winds moves the sands from the lakeshore to the hilly areas, while fluvial erosion and secondary-prevailing winds from the east transport sands towards the lake following the terrain downhill.

Figure 1 Location map of lake sediments (a) (AC: Ahung Co (Morrill et al., 2006); BC: Bangong Co (Gasse et al., 1996); CL:Cuoe Lake (Wu et al., 2006); KL: Kuhai Lake (Mischke et al., 2010); NC: Nam Co (Mügler et al., 2010); QL: Qinghai Lake (Shen et al., 2005; An et al., 2012; Lu et al., 2015); SL: Sugan Lake (Chen et al., 2009); SLC: Siling Co (Morinaga et al., 1993); SUC: Sumxi Co (Gasse et al., 1991); XL: Xinmucun Lake (Liang and Jiang, 2017)) and aeolian sediments (DC: Dongji Cona (Ijmker et al., 2012));GB: Gonghe Basin (Liu et al., 2013; Qiang et al., 2016); QB: Qaidam Basin (Yu and Lai, 2014); QL: Qinghai Lake (Lu et al., 2011,2015) in previous studies on the TP; (b) ambience environment around Lake Co Nag; (c) sampling site in this study

The altitude is about 4,800 m and annual mean temperature is ?2.57 °C. Records from the meteorological station in Anduo County, 23 km northeast of Lake Co Nag, show that annual precipitation is 446 mm with 70% in summer, 20% in autumn, 9% in spring and only about 1% in winter. The temperature raising trend is 0.357 °C per decade from 1966–2013 and 0.491 °C per decade from 1979–2013, which is almost twice the reported global land surface air temperature increase from CRUTEM (Climatic Research Unit TEMperature) and Berkeley models (Christensenet al., 2013). The warming may lead to a temperature increase of top permafrost, the melting of soil ice, the drying and loosing of surface sand, and thus further surface degradation and desertification.

An isolated outcropping sandy mound approximately 4 m above the water surface was found on the eastern lakeshore sand field (91°32'22"E, 32°04'09''N)(Figure 1c). The surrounding material had been eroded to form part of the sediment along the east railway or carried down west to the low reach of the lake. In August 2014, a 2.7-m sampling profile beginning at the top of the mound was created; samples were collected using a hand shovel every 5 cm interval from the bottom of the profile. The profile was composed of purely sand sediments. In total, 57 sand samples weighing 500 g each were collected. Meanwhile, six OSL samples to determine depositional ages were collected from the profile (Figure 1c). The six samples were obtained by hammering the samplers (aluminum tubes of 30-cm long and 5 cm in diameter) into a freshly created vertical section and sealing them with black plastic tape to avoid light exposure and moisture loss.

2.2 Grain-size analysis

Grain-size and geochemical analyses are two convenient methods to study sediment features and their depositional environment (Liuet al., 2014; Liuet al.,2015). The grain-size distribution of all the samples was determined between 0.02 and 2,000 μm by a Malvern Mastersizer laser grain-size analyzer in Lanzhou University, China. The accuracy and reproducibility of the analyzer was better than 1%, respectively. Values were converted toΦ(phi) units using the equationΦ=?log2(d), wheredis the grain-size diameter in mm. Before each measurement, chemical pretreatment following the procedure in Appendix C of Konert and Vandenberghe (1997) was performed to isolate discrete particles, removal of organic matter and carbonates using H2O2and HCl, respectively, and ensure an evenly dispersed suspension.

Grain-size parameters including mean grain size(Mz), standard deviation or sorting (σ), skewness(SK1), and kurtosis (KG) were calculated based on the method of Folk and Ward (1957). Sediment depositional environment can be determined by using a statistical method based on grain-size parameters.

Linear discrimination using grain-size parameters is an effective way to distinguish aeolian sediments.The discriminant functions proposed by Sahu (1964)were applied by several authors to the sand samples taken from the Kumtagh Desert (Heet al., 2009;Donget al., 2011; Liuet al., 2014). The four functions used are as follows:

Y4is less than 9.843, 3 indicates turbidity current deposition; otherwise, it is fluvial deposition.

The four grain-size parameters were substituted into these functions. The results present clear differences among different sediments.

2.3 Geochemical and chronological analysis

A fully automated sequential wavelength dispersive X-ray fluorescence spectrometer (AXIOS, PANalytical B.V., the Netherlands) was used for elemental analysis at the Key Laboratory of Desert and Desertification, Chinese Academy of Sciences (KLDD-CAS).Accessories include a semi-automatic press machine,grinding mill with tungsten carbon, and closed circuit cooling unit (ZHY-401A, ZHM-1A, BLK2-8FF-R,respectively, Zhonghe Corporation, Beijing, China).The samples were crushed into a powder finer than 75 μm using a multipurpose grinder and dried in an oven at 105 °C. Next, 4 g of the dry, powdered sample was pressed into 32 mm-diameter pellets under 30 tons of pressure using the pressed powder pellet technique. The briquettes were then stored in desiccators and tested in the spectrometer. Thirty chemical elements served as tracers were analyzed including Cl, P, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Ga, As, Br,Rb, Sr, Y, Zr, Nb, Ba, La, Ce, Nd, Pb, SiO2, Al2O3,Fe2O3, MgO, CaO, Na2O, and K2O. The oxide concentrations were in percentage, while other elements were in units of μg/g.

Several geochemical parameters including chemical index of alteration (CIA), SiO2/TiO2and CaO/MgO were calculated based on the chemical analysis data. The CIA index is calculated as:

where the materials are expressed by their molar fraction (m). The value of CaO* is taken from m(Na2O)if m(CaO) ＞ m(Na2O). This is done to consider CaO only from silicates while removing the content in carbonate or phosphate, as the molar fractions of CaO and Na2O are usually equal in silicates.

The OSL measurements were carried out using an automated Ris?TL/OSL-DA-15 reader, also in the KLDD-CAS. The OSL signal from quartz was detected through one 5-mm thick Hoya U-340 filter. The IRSL signal from feldspar was detected by combining one BG-39 and one Coring-759 filters. Laboratory irradiation was carried out using90Sr/90Y sources mounted within the reader with a dose rate of 0.084 Gy/s.Details on the procedures of OSL measurements can be found in Zhaoet al. (2015). OSL results are presented in Table 1.

3 Results

3.1 Grain-size and parameters

The profile samples are coarse in particle size with particles ＞63 μm being at least 98.3% and an average of 99.7%. No samples have particles of ＜16 μm.No clay and silt were found in all samples, and the soil texture of all the sediments is classified as sand.The sediments are classified into fine sand group(72.9%±4.5%), followed by very fine sand group(14.6%±4.0%) and medium sand group (12.2%±3.5%).

Deposition environment discrimination analysis was applied to all the samples. Results indicate that all samples were formed under an aeolian environment except for the third from the top sample (No.55), which has a higher coarse particle ratio and tends to be shallow marine deposition. This outlier might be impacted by stronger lake waves in the near past compared to other layers, but overall deposition materials should still be transported by an aeolian process.

Table 1 OSL dating results and related parameters of aeolian sediments at Lake Co Nag

The uniform sediment environment and source are further indicated by particle parameters. The scatter plots ofσ,SK1andKGagainstMzall show close sorting, skewness, and kurtosis of the samples except for No. 55 sample outlier (Figure 2). Mean grain sizes are also concentrated in a narrow range of about 2.5 to 2.8Φ, corresponding to 0.17–0.14 mm. Theσvalue ranges from 0.50 to 0.65, representing relatively good sorting (Figure 2a).SK1values of close to 0 indicate symmetric distribution of the sample particles (Figure 2b).KGvalues of 0.90–1.05 show that the sediments have medium ratios of middle-size particles to the tail particles (Figure 2c). The concentrated grain size, sorting, skewness, and kurtosis all indicate that the sediments were separated by wind power. As a result, the variations of aeolian material characteristics can be used to reflect past aeolian activity conditions and perhaps project future desertification trends with additional information from climate data.

Figure 2 Scatter plot of statistical grain size parameters. (a): Mz vs. σ; (b): Mz vs. SK1, and (c): Mz vs. KG

3.2 Geochemical elements

For tracer elements, all the layers are rich in Co and leaching in all others (Figure 3a). The Co element concentration is six times higher than the Upper Continental Crust (UCC) (Rudnick and Gao, 2014),10 times higher than the background concentration on the TP (Zhanget al., 2002), 20 times the desert sands to the north of the TP (Su and Yang, 2008), and 9 times than the average of Chinese soils (Chenet al., 1991).The element richness characters together with small coefficient of variation suggest that lakeshore sand deposition is derived from a single and uniform protolith, ultrabasic igneous rock (Jowittet al., 2012). It can be inferred that the sediments were originally from the weathering of mountain rocks to the east of the lake, which were then transported by fluvial and wind processes toward the lake.

For major elements (expressed as oxides), the layers are significantly rich in Ca and slightly lower in Si while leaching in others (Figure 3b). The enrichment of Ca indicates a dry climatic condition with weak weathering, and Si is enriched in cold and dry conditions with less weathering and decreased leaching while its content is positively correlated with coarse particle content in sediments (Liuet al., 2013). However, there is no clear correlation of Si content with the coarse sediment content in the samples because the latter has uniform high values. Al is relatively stable and is separated out at the end of the chemical weathering process, and its content (Al2O3ranges at 5.8%–6.5%) is clearly lower than sands in other deserts (Zimbelman and Williams, 2002; Liuet al.,2015).

It is interesting to note that two distinct synchronous changes are found from the profiles of many elements, a stronger one at about 90 years before present(B.P.) and a weaker one at about 380 years B.P.(Figure 4). Si content at the two periods are found to be significantly lower than in adjacent samples, while other linked elements such as Mn, Co, Ca and Fe are clearly higher. This feature is not found in other tested elements. The corresponding changes of other oxides including MgO, CaO, and Fe2O3can be responses to SiO2. Because Si existed as oxides, leaching was easier in warm or humid conditions due to enhanced weathering. These changes may indicate significant warm or wet periods for the corresponding sedimentary layers. The temperature condition was generally known as colder, especially during the LIA, than present; thus, we can assume that the two strong synchronous profile changes were attributed to strong wet conditions, which may have occurred before the dated time because of lagged responses of geochemical elements.

Figure 3 UCC normalized chemical composition of the sediments

Figure 4 Profiles of several elements that show two synchronous variations at about 90 and 380 years B.P.

3.3 Geochemical parameters

CIA reflects the chemical weathering of silicate minerals. A higher CIA mirrors a humid climate and stronger geochemical weathering with more leaching of elements such as Na, K, and Ca. Results show that CIA values are low for all the samples and vary modestly from 49.8 to 53.5, indicating low level chemical weathering through the whole profile (Figure 5).However, they are still higher than the weathering intensity of sands collected from deserts (Donget al.,2011), which can be explained by the relatively high sediment water content due to their closeness to the lake. A decreasing tendency was found from the bottom to the top of the profile, indicating recent weather condition is drier or colder than the past. This slightly decreasing trend before about 150–160 years B.P. was caused by a cold and dry environment during the late LIA, while the more clearly decreasing after the end of LIA till recent time was due to a significant drying caused by warming climate (Yanget al., 2002;Jinet al., 2007; Ijmkeret al., 2012; Gouet al., 2015).

Figure 5 Profiles of element parameters with linear regressions for the past four centuries that show a cold period during the LIA and a drying change since about 160 years B.P.

Ti content is usually stable because it mainly exists in stable minerals. Si accumulates in cold and dry environments, as mentioned above. Therefore, an increased SiO2/TiO2ratio indicates a colder or drier climate, while a decreased ratio reflects enhanced weathering. The SiO2/TiO2values of the profile generally increases with fluctuation from the bottom of the profile, and reflects a drier recent climate compared to about 400 years B.P. (Figure 5). However the ratio was slightly reduced before 150–160 years B.P. and quickly increased afterwards, corresponding to a relative stable or lower dryness condition during the late LIA with a cold and rapidly drying trend after LIA with warming temperatures since about 150–160 years B.P. (Zhanget al., 2016).

The Ca2+/Mg2+ratio is also a good indicator of an arid climate, because these two elements are typical in arid and semi-arid environments, though Ca2+migrates easier than Mg2+. Therefore, an increasing CaO/MgO ratio represents a drier environmental changing trend. Variation of the CaO/MgO profile is synergetic to and more significant compared to that of the SiO2/TiO2profile, indicating the consistency of these geochemical parameters (Figure 5).

3.4 Change of deposition rate

The deposition rate of sediments can also be analyzed using the OSL dating and core depth information to trace past aeolian environments due to its sensitivity to alternations of aridity and wetness, and to related atmospheric circulations (Luet al., 2005; Longet al., 2016). The transportation distance was short in the study area so that the blown sand can reach the sedimentary site in one wind event or a short period.Thus, the deposition rate can imply the richness of loose surface material and general climate conditions at the time of sand accumulation. Using the six OSL dating values and sample depths, we obtained the time span and depth change between two adjacent samples;thereafter, the rate of sand deposition change (cm/a) at five different depositional periods can be estimated.

Results show that the deposition rate of sand sediments increased quickly in the past four centuries.The rate was 0.29 cm/a in the period of 380–240 years B.P., but with a nearly 6 times increase to 1.63 cm/a in the most recent period of 60–20 years B.P.. This significant change can be expressed by a power function with high confidence level (R2=0.90), in whichxis time (years B.P.) andyis the deposition rate(Figure 6). The change of rate indicates that aeolian activity became much more prevalent during the tested period with significantly increased sands available at the land surface. However, there was a stagnation of deposition rate change between 140 and 90 years B.P., which means the erosibility of surface materials was low although it was under a warming background. This can be explained by the wet period before 90 years B.P. as indicated in section 3.2, which may increase surface moisture and plant growth,while reducing aeolian activity over numerous years in mostly winter and spring dry seasons.

4 Discussion

4.1 Source of lakeshore depositions

Gain-size analysis with parameters such as mean grain size, sorting, skewness and kurtosis, together with special geochemical characters such as Co richness and absence of Zn, indicates that lakeshore deposition was transported by the aeolian process and their same origin was the weathering of local mountain rocks east of the lake. This is consistent with other studies that show local and small-scale source of the TP sands (Ijmkeret al., 2012). These results also agree with a previous study on the sources of lakeshore sands based on composite element fingerprinting analysis (Liuet al., 2016b).

Figure 6 Depositional rate change of sand sediments with time

Alpine meadows at different degraded levels can be seen throughout the mountains, and aeolian material is plentiful although the natural weathering process is slow in this cold area. The mountain surface was devoid of developed soil during the field survey, and the top several meters of sands may have been totally removed by wind (Liuet al., 2016a). The dominant wind direction was westerly; however, the downhill transportation of sands from the eastern mountain area should have a significant lower threshold, which means the movement of sands toward the lake is much easier with secondary prevailing easterly wind.The accumulated aeolian sands rest at or beside the gullies can also be carried by floods during the rainy season and deposited along the lakeshore.

A similar sand material transportation and the interplay of fluvial/aeolian processes may be common to the many other lakes over the TP, as a large part of this high plateau is influenced by the westerlies.

4.2 Recent climate change

Several geochemical parameters including CIA and ratios of SiO2/TiO2and CaO/MgO indicate that climate at Lake Co Nag during the past four centuries was impacted by LIA and significant drying after the end of LIA in the 20th century. LIA cold temperature resulted in low weathering which is reflected by a slightly reduced CIA. The frozen surface, low evaporation, and high soil moisture in the sand surface that resulted in a slight reduction in dryness is reflected by ratios of SiO2/TiO2and CaO/MgO. After LIA, temperatures started to increase over the TP with a significant warming in the 20th century (Kanget al., 2015;Zhanget al., 2016). Warming in the dry winter and spring seasons can increase evaporation and reduce surface soil moisture, which leads to higher erodibility of surface sands, while the change of precipitation in the wet season has a limited effect. This is reflected by geochemical elements with lower CIA, higher SiO2/TiO2and CaO/MgO values.

Recorded meteorological data for the past half century show that precipitation and temperature increased at the same time; however, the increased rainfall occurred mainly in the summer rainy season,while in the dry winter and spring with strong wind activity, dry and bare surfaces received little additional water (Liuet al., 2016a; Zhanget al., 2016). Another hypothesis is that the added plants that responded to the increase of summer rainfall may intercept more shifting sands and lead to an increase of sand depositional rate. However, aeolian activity was mainly dominant in winter and spring, while the ephemeral short grass in summer has little affect as no biological residual remained in the studied sand profile. Even increasing precipitation in summer might not be consistent in other areas, and it remains difficult to evaluate the precipitation change trend in the face of limited observational data in the TP (Conroyet al., 2017). This is because the summer was dominated by the Indian monsoon and the winter and spring were controlled by the westerlies, and the two climate systems have different patterns.

A similar conclusion in the southern TP was found where temperature instead of precipitation is the primary factor on local evapotranspiration and effective moisture (Klinge and Lehmkuhl, 2015). As a result, the dry season became warmer and drier, resulting in more available sands and stronger aeolian activity over the eastern mountain area and more sands transported west to the lake. This increasing richness of surface sands caused by the warming in the dry seasons was verified by the change of deposition rate of sands at the lakeshore of Lake Co Nag,which increased nearly six times at about 60–20 years ago compared to 380–240 years ago. This trend might have ceased in some wet periods with higher surface moisture conditions in the dry seasons.

The increase in the deposition rate of sands with the warming trend in this area suggests that the strengthening of surface sand erodibility is sensitive to climate change. Significant surface aeolian activity and desertification may be triggered by climatic forces in the central TP, resulting in problems including strong sand depositional at the two sides of the QTR. This condition may be worse in the future as there is high confidence that temperatures are projected to increase (Christensenet al., 2013). Thus, more attention and action is needed to address the challenging future of increased sand and dust.

4.3 Evidence of significant wet perio ds

The wet/dry episodes over the TP are explained by the weakening of monsoons from East Asia since 4.2 cal ka B.P. (Mügleret al., 2010; Chenet al.,2016). This weakening of Asian summer monsoon can lead to strengthening of aeolian activity on the TP andvice versa(Stauchet al., 2016, 2017). The coordinated changes of deposition rate and geochemical profile indicate a significant wet period between 90 and 140 years B.P., which led to a cessation of deposition rate increase and abrupt changes of many elemental contents. There was another older wet period at about 380 years B.P. that can be inferred from a similar geochemical profile change, although we have no deposition rate data to verify this because no older age dating values were available. Thus, mutual verification of physical and chemical information indicates that geochemical data, together with supporting data from grain size and age dating data, can be used to detect ancient climate events.

The recent wet period between 90 and 140 years B.P. was confirmed by a wet phase during 1883–1906 according to historical archives (Lin and Wu, 1986).Some studies revealed a similar pronounced wet condition, for example during 1890–1900 in the southeast (Liet al., 2017), 1891–1913 in south-central (Heet al., 2013), and 1886–1913 in north (Wanget al.,2013) of TP. In addition, sometimes a relative dry period was also shown at various locations over TP(Chenet al., 2009; Griessingeret al., 2011). This assumption can hardly be verified by other high-precision precipitation reconstruction studies, because proxy data with high temporal resolution and reliability like those from tree rings are not available in central TP close to our study area. However, results in this study still raise the possibility of understanding the historical climate and environmental evolution in the remote plateau area. Future evidence may still come from soil or lake sediments, as they can provide accurate long term records of past climate (Holmeset al., 2009).

5 Conclusions

The TP is sensitive to climate change and is therefore a good indicator of climate change, but available paleoclimate information with high temporal precision is rare in the remote and harsh central TP. In this study, grain-size analysis with size parameters and deposition environment discriminant, geochemical element analysis with deposition profile change and chemical parameters, and sediment age dating with deposition rate estimation of a sand profile were combined to examine the environment and climate change with seasonal differences at the shore of Lake Co Nag in central TP during the past four centuries.

The physical and chemical factors support each other in showing changes of aeolian sediments and climatic forcing. The results suggest an agreement with the broader picture of LIA and 20th century environmental evolution. This study adds data to the overall picture of past climate and geomorphic history of the TP in showing that when other proxies are not available, aeolian sediments can reflect climatic variation and surface process with multi information including physical, chemical, and dating data.

Some facts in this study can be drawn as follows:

(1) The depositional rate of aeolian sands increased by six times during the past four centuries resulting from increasingly active surface sands, and problems related to aeolian activity such as surface erosion, blown-sand deposition and desertification.This type of activity may worsen in the future under a projected warming climate which may occur in other cold and high altitude areas. The source of the east lakeshore sediments of Co Nag was the mountain area to the east from ultrabasic igneous rock.

(2) Collaborative analysis using chemical, depositional and meteorological data indicates different seasonal effects of climate variation on geomorphic processes. Aeolian activity is very sensitive to warming in dry seasons in the study area. The dry seasons of winter and spring are warmer and drier compared to 400 years ago although the summer season might have received more precipitation. Climatic conditions before and after the LIA were clearly different as indicated by geochemical parameters including CIA and ratios of SiO2/TiO2and CaO/MgO.

(3) Two significant wet periods might have happened between 90–140 years B.P. and about 380 years B.P. in the central TP, which led to abrupt changes of the geochemical profile and pauses of deposition rate increase.

We thank Pro. Hui Zhao and Dr. JiangLin Wang at the Northwest Institute of Eco-Environment and Engineering, Chinese Academy of Sciences for discussions on past climatic events in the Tibetan Plateau.This study was supported by the National Science Fund of China (41501008), the China Postdoctoral Science Foundation (2014M550518), the Youth Innovation Promotion Association (2016373), and the"Light of West China" Program of the Chinese Academy of Sciences.

An ZS, Colman SM, Zhou WJ,et al., 2012. Interplay between the Westerlies and Asian monsoon recorded in Lake Qinghai sediments since 32 ka. Scientific Reports, 2: 619. DOI: 10.1038/srep00619.

Br?uning A, 2006. Tree-ring evidence of 'Little Ice Age' glacier advances in southern Tibet. The Holocene, 16(3): 369–380. DOI:10.1191/0959683606hl922rp.

Chen FH, Wu D, Chen JH,et al., 2016. Holocene moisture and East Asian summer monsoon evolution in the northeastern Tibetan Plateau recorded by Lake Qinghai and its environs: A review of conflicting proxies. Quaternary Science Reviews, 154: 111–129. DOI:10.1016/j.quascirev.2016.10.021.

Chen JH, Chen FH, Zhang EL,et al., 2009. A 1000-year chironomidbased salinity reconstruction from varved sediments of Sugan Lake, Qaidam Basin, arid Northwest China, and its palaeoclimatic significance. Chinese Science Bulletin, 54(20): 3749–3759. DOI:10.1007/s11434-009-0201-8.

Chen JS, Wei FS, Zheng CJ,et al., 1991. Background concentrations of elements in soils of China. Water, Air, & Soil Pollution,57–58(1): 699–712. DOI: 10.1007/BF00282934.

Cheng GD, Wu TH, 2007. Responses of permafrost to climate change and their environmental significance, Qinghai-Tibet Plateau. Journal of Geophysical Research: Earth Surface, 112(F2): F02S03. DOI:10.1029/2006JF000631.

Christensen JH, Kanikicharla KK, Marshall G,et al., 2013. Climate phenomena and their relevance for future regional climate change.In: Stocker TF, Qin DH, Plattner GK (eds.). Climate Change 2013:The Physical Science Basis. Cambridge: Cambridge University Press, pp. 1217–1308.

Conroy JL, Hudson AM, Overpeck JT,et al., 2017. The primacy of multidecadal to centennial variability over late-Holocene forced change of the Asian Monsoon on the southern Tibetan Plateau.Earth and Planetary Science Letters, 458: 337–348. DOI:10.1016/j.epsl.2016.10.044.

Dong Z, Su ZZ, Qian GQ,et al., 2011. Aeolian Geomorphology in the Kumtagh Desert. Beijing: Science Press.

Folk RL, Ward WC, 1957. Brazos River Bar: A study in the significance of grain size parameters. Journal of Sedimentary Research,27(1): 3–26. DOI: 10.1306/74D70646-2B21-11D7-8648000102 C1865D.

Gasse F, Arnold M, Fontes JC,et al., 1991. A 13,000-year climate record from western Tibet. Nature, 353(6346): 742–745. DOI:10.1038/353742a0.

Gasse F, Fontes JC, Van Campo E,et al., 1996. Holocene environmental changes in Bangong Co basin (Western Tibet). Part 4: Discussion and conclusions. Palaeogeography, Palaeoclimatology, Palaeoecology, 120(1–2): 79–92. DOI: 10.1016/0031-0182(95)00035-6.

Gou XH, Gao LL, Deng Y,et al., 2015. An 850-year tree-ring-based reconstruction of drought history in the western Qilian Mountains of northwestern China. International Journal of Climatology,35(11): 3308–3319. DOI: 10.1002/joc.4208.

Griessinger J, Brauning A, Helle G,et al., 2011. Late Holocene Asian summer monsoon variability reflected by δ18O in tree-rings from Tibetan junipers. Geophysical Research Letters, 38(3): L03701.DOI: 10.1029/2010GL045988.

He MH, Yang B, Br?uning A,et al., 2013. Tree-ring derived millennial precipitation record for the south-central Tibetan Plateau and its possible driving mechanism. The Holocene, 23(1): 36–45.

He Q, Yang XH, Huo W,et al., 2009. Characteristics of sand granularity from kumtag desert and its environmental significance. Journal of Desert Research, 29(1): 18–22.

Holmes JA, Cook ER, Yang B, 2007. Climate change over the past 2000 years in Western China. Quaternary International, 194(1–2):91–107. DOI: 10.1016/j.quaint.2007.10.013.

Ijmker J, Stauch G, P?tsch S,et al., 2012. Dry periods on the NE Tibetan Plateau during the late Quaternary. Palaeogeography, Palaeoclimatology, Palaeoecology, 346–347: 108–119. DOI:10.1016/j.palaeo.2012.06.005.

Jin HJ, Chang XL, Wang SL, 2007. Evolution of permafrost on the Qinghai-Xizang (Tibet) Plateau since the end of the late Pleistocene. Journal of Geophysical Research: Earth Surface, 112(F2):F02S09. DOI: 10.1029/2006JF000521.

Jowitt SM, Jenkin GRT, Coogan LA,et al., 2012. Quantifying the release of base metals from source rocks for volcanogenic massive sulfide deposits: Effects of protolith composition and alteration mineralogy. Journal of Geochemical Exploration, 118: 47–59.DOI: 10.1016/j.gexplo.2012.04.005.

Kang S, Wang F, Morgenstern U,et al., 2015. Dramatic loss of glacier accumulation area on the Tibetan Plateau revealed by ice core tritium and mercury records. The Cryosphere, 9(3): 1213–1222. DOI:10.5194/tc-9-1213-2015.

Klinge M, Lehmkuhl F, 2015. Holocene aeolian mantles and inter-bedded paleosols on the southern Tibetan Plateau. Quaternary International, 372: 33–44. DOI: 10.1016/j.quaint.2014.06.068.

Konert M, Vandenberghe J, 1997. Comparison of laser grain size analysis with pipette and sieve analysis: a solution for the underestimation of the clay fraction. Sedimentology, 44(3): 523–535. DOI:10.1046/j.1365-3091.1997.d01-38.x.

Li JB, Shi JF, Zhang DD,et al., 2017. Moisture increase in response to high-altitude warming evidenced by tree-rings on the southeastern Tibetan Plateau. Climate Dynamics, 48(1–2): 649–660. DOI:10.1007/s00382-016-3101-z.

Liang LJ, Jiang HC, 2017. Geochemical composition of the last deglacial lacustrine sediments in East Tibet and implications for provenance, weathering, and earthquake events. Quaternary International,430: 41–51. DOI: 10.1016/j.quaint.2015.07.037.

Lin ZY, Wu XD, 1986. A preliminary analysis of the regularity in flood, drought and snowstorm in Tibetan Plateau during historical times. Acta Meteorologica Sinica, 44(3): 257–264.

Liu B, Jin HL, Sun LY,et al., 2013. Holocene climatic change revealed by aeolian deposits from the Gonghe Basin, northeastern Qinghai–Tibetan Plateau. Quaternary International, 296: 231–240.DOI: 10.1016/j.quaint.2012.05.003.

Liu B, Jin HL, Sun LY,et al., 2015. Geochemical evidence for Holocene millennial-scale climatic and environmental changes in the south-eastern Mu Us Desert, northern China. International Journal of Earth Sciences, 104(7): 1889–1900. DOI: 10.1007/s00531-015-1161-7.

Liu BL, Qu JJ, Ning DH,et al., 2014. Grain-size study of aeolian sediments found east of Kumtagh Desert. Aeolian Research, 13: 1–6.DOI: 10.1016/j.aeolia.2014.01.001.

Liu BL, Qu JJ, Kang SC, 2016a. Response of dune activity on the Tibetan Plateau to near future climate change. Climate Research,69(1): 1–8. DOI: 10.3354/cr01385.

Liu BL, Niu QH, Qu JJ,et al., 2016b. Quantifying the provenance of aeolian sediments using multiple composite fingerprints. Aeolian Research, 22: 117–122. DOI: 10.1016/j.aeolia.2016.08.002.

Liu JJ, Yang B, Qin C, 2011. Tree-ring based annual precipitation reconstruction since AD 1480 in south central Tibet. Quaternary International, 236(1–2): 75–81. DOI: 10.1016/j.quaint.2010.03.020.

Liu XD, Chen BD, 2000. Climatic warming in the Tibetan plateau during recent decades. International Journal of Climatology, 20(14):1729–1742. DOI: 10.1002/1097-0088(20001130)20:14＜1729::AID-JOC556＞3.0.CO;2-Y.

Long H, Fuchs M, Yang LH,et al., 2016. Abrupt sand-dune accumulation at the northeastern margin of the Tibetan Plateau challenges the wet MIS3a inferred from numerous lake-highstands. Scientific Reports, 6: 25820. DOI: 10.1038/srep25820.

Lu HY, Miao XD, Zhou YL,et al., 2005. Late Quaternary aeolian activity in the Mu Us and Otindag dune fields (north China) and lagged response to insolation forcing. Geophysical Research Letters, 32(21): 2465–2475. DOI: 10.1029/2005GL024560.

Lu HY, Zhao CF, Mason J,et al., 2011. Holocene climatic changes revealed by aeolian deposits from the Qinghai Lake area (northeastern Qinghai-Tibetan Plateau) and possible forcing mechanisms.The Holocene, 21(2): 297–304. DOI: 10.1177/0959683610378884.Lu RJ, Jia FF, Gao SY,et al., 2015. Holocene aeolian activity and climatic change in Qinghai Lake basin, northeastern Qinghai-Tibetan Plateau. Palaeogeography Palaeoclimatology, Palaeoecology, 430:1–10. DOI: 10.1016/j.palaeo.2015.03.044.

Mischke S, Zhang CJ, Borner A,et al., 2010. Lateglacial and Holocene variation in aeolian sediment flux over the northeastern Tibetan Plateau recorded by laminated sediments of a saline meromictic lake. Journal of Quaternary Science, 25(2): 162–177.DOI: 10.1002/jqs.1288.

Morinaga H, Itota C, Isezaki N,et al., 1993. Oxygen-18 and carbon-13 records for the last 14, 000 years from Lacustrine carbonates of Siling-Co (Lake) in the Qinghai-Tibetan Plateau. Geophysical Research Letters, 20(24): 2909–2912. DOI: 10.1029/93GL02982.

Morrill C, Overpeck JT, Cole JE,et al., 2006. Holocene variations in the Asian monsoon inferred from the geochemistry of lake sediments in central Tibet. Quaternary Research, 65(2): 232–243. DOI:10.1016/j.yqres.2005.02.014.

Mügler I, Gleixner G, Günther F,et al., 2010. A multi-proxy approach to reconstruct hydrological changes and Holocene climate development of Nam Co, Central Tibet. Journal of Paleolimnology, 43(4):625–648. DOI: 10.1007/s10933-009-9357-0.

Qiang MR, Jin YX, Liu XX,et al., 2016. Late Pleistocene and Holocene aeolian sedimentation in Gonghe Basin, northeastern Qinghai-Tibetan Plateau: Variability, processes, and climatic implications.Quaternary Science Reviews, 132: 57–73. DOI: 10.1016/j.quascirev.2015.11.010.

Qin J, Yang K, Liang SL,et al., 2009. The altitudinal dependence of recent rapid warming over the Tibetan Plateau. Climatic Change,97(1–2): 321–327. DOI: 10.1007/s10584-009-9733-9.

Rudnick RL, Gao S, 2014. 4.1-Composition of the continental crust A2-Holland, Heinrich D. In: Turekian KK (ed.). Treatise on Geochemistry (2nd ed.). Oxford: Elsevier, pp. 1–51. DOI: 10.1016/B978-0-08-095975-7.00301-6.

Sahu B K, 1964. Depositional mechanisms from the size analysis of clastic sediments. Journal of Sedimentary Research, 34(1): 73–83.

Shen J, Liu XQ, Matsumoto R,et al., 2005. A high-resolution climatic change since the Late Glacial Age inferred from multi-proxy of sediments in Qinghai Lake. Science in China Series D: Earth Sciences, 48(6): 742–751. DOI: 10.1360/03yd0148.

Stauch G, 2015. Geomorphological and palaeoclimate dynamics recorded by the formation of aeolian archives on the Tibetan Plateau.Earth-Science Reviews, 150: 393–408. DOI: 10.1016/j.earscirev.2015.08.009.

Stauch G, Lai Z, Lehmkuhl F,et al., 2016. Environmental changes during the late Pleistocene and the Holocene in the Gonghe Basin,north-eastern Tibetan Plateau. Palaeogeography Palaeoclimatology Palaeoecology, in press. DOI: 10.1016/j.palaeo.2016.12.032.

Stauch G, Schulte P, Ramisch A,et al., 2017. Landscape and climate on the northern Tibetan Plateau during the late Quaternary. Geomorphology, 286: 78–92. DOI: 10.1016/j.geomorph.2017.03.008.

Su YZ, Yang R, 2008. Background concentrations of elements in surface soils and their changes as affected by agriculture use in the desert-oasis ecotone in the middle of Heihe River Basin, Northwest China. Journal of Geochemical Exploration, 98(3): 57–64.DOI: 10.1016/j.gexplo.2007.12.001.

Wang WZ, Liu XH, Xu GB,et al., 2013. Moisture variations over the past millennium characterized by Qaidam Basin tree-ring δ18O.Chinese Science Bulletin, 58(32): 3956–3961. DOI: 10.1007/s 11434-013-5913-0.

Wu GX, Duan AM, Liu YM,et al., 2015. Tibetan Plateau climate dynamics: recent research progress and outlook. National Science Review, 2(1): 100–116. DOI: 10.1093/nsr/nwu045.

Wu YH, Lücke A, Jin ZD,et al., 2006. Holocene climate development on the central Tibetan Plateau: A sedimentary record from Cuoe Lake. Palaeogeography, Palaeoclimatology, Palaeoecology,234(2–4): 328–340. DOI: 10.1016/j.palaeo.2005.09.017.

Xie SB, Qu JJ, Wang T, 2016. Wind tunnel simulation of the effects of freeze-thaw cycles on soil erosion in the Qinghai-Tibet Plateau.Sciences in Cold and Arid Regions, 8(3): 187–195. DOI:10.3724/SP.J.1226.2016.00187.

Xue X, Guo J, Han BS,et al., 2009. The effect of climate warming and permafrost thaw on desertification in the Qinghai–Tibetan Plateau.Geomorphology, 108(3–4): 182–190. DOI: 10.1016/j.geomorph.2009.01.004.

Yang B, Braeuning A, Johnson KR,et al., 2002. General characteristics of temperature variation in China during the last two millennia.Geophysical Research Letters, 29(9): 1324. DOI: 10.1029/2001 GL014485.

Yang B, Br?uning A, Shi YF, 2003. Late Holocene temperature fluctuations on the Tibetan Plateau. Quaternary Science Reviews,22(21–22): 2335–2344. DOI: 10.1016/S0277-3791(03)00132-X.

Yao TD, Thompson L, Yang W,et al., 2012. Different glacier status with atmospheric circulations in Tibetan Plateau and surroundings.Nature Climate Change, 2(9): 663–667. DOI: 10.1038/nclimate 1580.

Yu KF, Lehmkuhl F, Falk D, 2017. Quantifying land degradation in the Zoige Basin, NE Tibetan Plateau using satellite remote sensing data. Journal of Mountain Science, 14(1): 77–93. DOI: 10.1007/s11629-016-3929-z.

Yu LP, Lai ZP, 2014. Holocene climate change inferred from stratigraphy and OSL chronology of aeolian sediments in the Qaidam Basin, northeastern Qinghai-Tibetan Plateau. Quaternary Research,81(3): 488–499. DOI: 10.1016/j.yqres.2013.09.006.

Zhang KC, Qu JJ, Liao KT,et al., 2010. Damage by wind-blown sand and its control along Qinghai-Tibet Railway in China. Aeolian Research, 1(3–4): 143–146. DOI: 10.1016/j.aeolia.2009.10.001.

Zhang XP, Deng W, Yang XM, 2002. The background concentrations of 13 soil trace elements and their relationships to parent materials and vegetation in Xizang (Tibet), China. Journal of Asian Earth Sciences, 21(2): 167–174. DOI: 10.1016/S1367-9120(02)00026-3.

Zhang YL, Kang SC, Grigholm B,et al., 2016. Twentieth-century warming preserved in a Geladaindong mountain ice core, central Tibetan Plateau. Annals of Glaciology, 57(71): 70–80. DOI:10.3189/2016AoG71A001.

Zhao H, Li SH, Li B,et al., 2015. Holocene climate changes in westerly-dominated areas of central Asia: Evidence from optical dating of two loess sections in Tianshan Mountain, China. Quaternary Geochronology, 30: 188–193. DOI: 10.1016/j.quageo.2015.04.002.

Zimbelman JR, Williams SH, 2002. Geochemical indicators of separate sources for eolian sands in the eastern Mojave Desert, California, and western Arizona. Geological Society of America Bulletin,114(4): 490–496. DOI: 10.1130/0016-7606(2002)114＜0490:GIOSSF＞2.0.CO;2.