Glacier-Induced Alluvial Fan Development on the Northeast Tibetan Plateau Since the Late Pleistocene

The origin of alluvial fans at the drainage basin or the margins of the arid sedimentary basin on the northeastern Tibetan Plateau (NETP) has been a focus of debate among scientific communities. Extensive alluvial fans had developed in the mountain-basin systems of this region during the late Pleistocene. Based on geomorphic/stratigraphic studies and the optically stimulated luminescence (OSL) chronology, we investigated numerous alluvial sequences in the drainage basin of the coupled mountain-basin system on the NETP. Sedimentologic analyses showed that these alluvial sediments mainly comprise crudely meter and decimeter beds of fanglomerate and sandy lenticles, which occasionally contain boulders in the thinner layers. OSL dating results showed that the alluvial sediments were mainly developed during the late MIS 5, MIS 3, and Last Glaciation and Deglaciation (maximum aggradations), while little had occurred during the early Holocene. Comparing our ages and their probability density curves with those of glacial advances from the northeastern/eastern TP, the consistency of both records during different stages since the late Pleistocene suggested that the development of alluvial fans was driven by glacier activities, which yielded abundant outwash, feeding alluvial aggradations. Based on this finding, together with results of previous studies, we argued that the dynamics of alluvial geomorphic processes had a paleoclimatic origin rather than surface uplifts. Further studies indicated that the precipitation of glaciations originated from the ISM northward intrusion during the late MIS 5, MIS 3, and Holocene; however, it was mainly influenced by the westerlies during the Last Glaciation and Deglaciation. Therefore, the development of alluvial fans on the NETP represents coupling signals of the ISM and the westerlies on the interglacial–glacial cycle timescales. The coupling evolution of glacier activities and alluvial aggradations in the drainage basins of mountain-basin systems on the NETP provides a model for assessing the linkages between regional geomorphic processes and atmospheric circulations on hemispheric scales.

The origin of alluvial fans at the drainage basin or the margins of the arid sedimentary basin on the northeastern Tibetan Plateau (NETP) has been a focus of debate among scientific communities. Extensive alluvial fans had developed in the mountain-basin systems of this region during the late Pleistocene. Based on geomorphic/stratigraphic studies and the optically stimulated luminescence (OSL) chronology, we investigated numerous alluvial sequences in the drainage basin of the coupled mountain-basin system on the NETP. Sedimentologic analyses showed that these alluvial sediments mainly comprise crudely meter and decimeter beds of fanglomerate and sandy lenticles, which occasionally contain boulders in the thinner layers. OSL dating results showed that the alluvial sediments were mainly developed during the late MIS 5, MIS 3, and Last Glaciation and Deglaciation (maximum aggradations), while little had occurred during the early Holocene. Comparing our ages and their probability density curves with those of glacial advances from the northeastern/eastern TP, the consistency of both records during different stages since the late Pleistocene suggested that the development of alluvial fans was driven by glacier activities, which yielded abundant outwash, feeding alluvial aggradations. Based on this finding, together with results of previous studies, we argued that the dynamics of alluvial geomorphic processes had a paleoclimatic origin rather than surface uplifts. Further studies indicated that the precipitation of glaciations originated from the ISM northward intrusion during the late MIS 5, MIS 3, and Holocene; however, it was mainly influenced by the westerlies during the Last Glaciation and Deglaciation. Therefore, the development of alluvial fans on the NETP represents coupling signals of the ISM and the westerlies on the interglacial-glacial cycle timescales. The coupling evolution of glacier activities and alluvial aggradations in the drainage basins of mountain-basin systems on the NETP provides a model for assessing the linkages between regional geomorphic processes and atmospheric circulations on hemispheric scales.

INTRODUCTION
Geomorphic processes in the drainage basins of river catchments offer insight into the response of mountain-basin surface dynamics to paleoclimate fluctuations (Owen et al., 2006;Chen et al., 2011;An et al., 2018a;Lehmkuhl et al., 2018;Li et al., 2020a;Tao et al., 2020). Using geomorphic processes to reconstruct paleoclimate changes often focuses on specific sedimentary or geomorphological archives such as fluvial terraces, eolian sediments, glacial landforms, lake shorelines, and alluvial fans . Geomorphic processes on the Northeastern Tibetan Plateau (NETP) have been intensively studied in recent decades due to its typical location of mountainbasin coupling and the interactional zone of diversiform atmospheric circulations (Owen et al., 2006;Wang et al., 2009;Chang et al., 2017;Stauch et al., 2017;An et al., 2018aAn et al., , 2018bLiu et al., 2019;Li et al., 2020b;Tao et al., 2020). For example, Owen et al. (2006) carried out systemic investigations on the alluvial fans, river terraces, and moraines along the eastern Kunlun Mountain slopes that indicated links between climatic amelioration, deglaciation, lake desiccation, and terraced alluvial fan evolution. Wang et al. (2009) likewise carried out an extensive optically stimulated luminescence (OSL) study of alluvial valley fills and associated terraces along the Kunlun Mountains and found synchronous aggradations in different rivers during the last glacial stage.
Similar to geomorphic units in the drainage basin, an alluvial fan is an important dynamic landform process, which commonly develops in areas of high topographic relief, especially in the transition zone between mountain-basin systems (An et al., 2018a). An alluvial fan is a cone-shaped deposit of silt, sand, and gravel that accumulates where a stream flows from an eroding upland catchment and debouches into an adjacent basin (Zhang et al., 2014). Alluvial fans trap sediment delivered from mountain source areas and exert an important control on the delivery of sediment to downstream environments, to axial drainages, and to arid sedimentary basins, and they preserve a sensitive record of environmental fluctuations within the high mountain source area (Ashworth, 2006). Alluvial processes are mainly influenced by tectonic activity, climate changes, and base-level controls. The understanding of such processes has progressed considerably in the past few decades (Owen et al., 2003(Owen et al., , 2006Zhang et al., 2011;An et al., 2018a), whereas effects of climate factors appear to exert an overwhelming control on late Quaternary sequences and sedimentary styles, especially in arid drainage basins (Gao et al., 2018). Alluvial sediments are sensitive to changes in orographic glacial activities and drainage hydrological conditions. They provide valuable information on paleoclimatic fluctuations of upstream areas (Ritter et al., 1995;Vandenberghe, 1995;Bridgland, 2010). Likewise, an alluvial fan is a medium of sediment transportation and groundwater discharge. For instance, high sediment supply to an alluvial fan resulted in the temporary blocking of the river channel and influenced the lake level of the adjacent Heihai Lake. This geomorphic process presumably occurred during or shortly after large glaciations (Stauch et al., 2017). The groundwater in the Qilian Mountains was recharged by high piedmonts, flowed through the alluvial fans from south to north, and was finally discharged several decades later from springs located along the toe of the fan (Guo et al., 2017). Recent studies have found that glacier advance was a key factor for the development of alluvial fans due to the abundant detritus production at the piedmont belt in the eastern Kunlun Mountains (Owen et al., 2006;An et al., 2018b;Cui et al., 2020). These researchers argued that the alluvial Sanchahe Formation filled in the Golmud River valley, described as "a set of periglacial fanglomerates," originated from the piedmont outwash yielded in the high Kunlun Mountains (Van Der Woerd et al., 2002). Therefore, it is important to investigate the association between glacial activities, alluvial fan developments, and atmospheric circulation changes on the orbital timescale of the late Quaternary and the spatial scale of the entire drainage basin of the NETP, through extensive geomorphic investigations and abundant chronology controls. If tectonic uplift were negligible for the formation of alluvial fans since the late Pleistocene, studying the climatic driving mechanism for the geomorphic process from the periglacial zone to the drainage basin would help us come up with a hypothesis that drainage alluvial geomorphic processes respond to climatic systems, moisture transportation, and orographic glacial activity on a hemispheric scale.
In this study, we conducted a geomorphologic and stratigraphic investigation using 59 OSL dating ages from a series of alluvial fans in widespread drainage basins on the NETP. The aims of our study included the following: (i) to establish a chronological framework of these alluvial fans since the last interglaciation, (ii) to establish links between the orographic glacial activity and alluvial geomorphic development, and (iii) to explore their paleoclimatic implication on the atmospheric circulation scale.

Study Area and Environmental Parameters
The Northeastern TP, with a series of high mountains, basins, alpine steppes, and deep-cutting gorges, covers an area of ∼580,000 km 2 (Bureau of Geology and Mineral Resources of Qinghai Province, 1991). The extensively developed basinmountain systems are the distinguishing topographical feature in this region (Kapp et al., 2011;An et al., 2018b).
The diastrophism in this region shows large-scale crustal shortening and strike-slip faults during the late Cenozoic . These faults bring up westnorthwest-oriented orogenic belts, such as the Kunlun Mountains, the Qilian Mountains, Ela Mountain, and Riyue Mountain, which separate the several large Cenozoic basins including the Qaidam, Gonghe, Qinghai Lake, and Xining basins ( Figure 1) Zhang et al., 2014). The alluvial fans associated with strong orographic weathering developed on the piedmonts along both sides of these high mountains. These alluvial fans are almost continuous within the tributary valleys and eventually join with the main valley fill that stretches along the margins of the basins (Owen et al., 2006;An et al., 2018b;Gao et al., 2018). In the southern margin of the Qaidam Basin, the alluvial landforms contain three subgeomorphic units: the piedmont alluvial fans along the slopes of the Kunlun Mountains, the main valley fill in the drainages, and the distal alluvial delta fan that cascades into the basin (An et al., 2018b). The highly dissected alluvial fans stretch out to the river valleys, and these are many tens of meters of alluvial fans that form the main valley fill. In the vicinity of the oasis preceding the mountain, a typical valley alluvial fan radiating from a series of large tributary valleys coalesces with the alluvial fill of the main drainage and then drains into the salt lakes in the Qaidam Basin and forms the distal alluvial delta fan (Owen et al., 2006). Five strath terraces have been developed within the alluvial fill (i.e., the Sanchahe Formation) in the middle reaches of the Golmud River, forming the spectacular fluvial landform. The Chaka Basin is part of the Gonghe Basin in the northeastern part of the TP, and it is separated from the Qaidam Basin by Ela Mountain (Liu et al., 2008). The basin is adjacent to the Qinghai Lake basin to the north by the Qinghainan Mountains (Zhang et al., 2012). The alluvial fans in the Chaka-Gonghe-Qinghai Lake regions are divided into three independent landform units, which include the north of the Ela Shan and both sides of the piedmonts of the Qinghainan Mountains. The northern alluvial fans in the Chaka-Gonghe Basin are characterized by bigger area and short lengths, but the southern Qinghai Lake area is relatively small (Zhang et al., 2014). The alluvial fans in the mountains of the eastern Qinghai Lake basin basically cascade from the lateral tributary valleys and stretch out to the main river drainage basin. These fans usually extend 10-20 km to form cliff terraces 30-65 m above the riverbed, and they likewise become detrital aggradations and valley fill (Guo et al., 2017;Li et al., 2018). Therefore, a large number of alluvial fans constitute the most impressive geomorphic processes in the mountain-basin systems of the NETP. These alluvial fans provide the ideal materials to study paleoclimatic/tectonic history in this region. Although earthquakes occur frequently in the region and the average strike-slip rate of the Kunlun-Qinling Fault in the Quaternary was up to ∼10 mm/ y (Van Der Woerd et al., 2002;Li et al., 2005), no indicators of tectonic uplift activity during the late Pleistocene, such as mole tracks, have been observed in the Golmud catchment (Li et al., 2005;Wang et al., 2009;Stauch et al., 2017).
The elevations of the sampling sites range from 2,700 to 4,200 m a.s.l. These elevation gradients present a trend of high in the west and low in the east. The drainage basins of the Kunlun and Qilian Mountains have the average elevations of ∼4,000 m a.s.l. But it is approximately ∼3,000 m a.s.l. in the Chaka-Gonghe Basin, the Qinghai Lake basin, and even much lower in the further southeast. Many types of glacial landforms and deposits have been developed in the U-shaped valleys, which are etched in the slopes on both sides of the mountains around the basins. The paleo-moraines are preserved at the elevations of 5,000-4,500 m a.s.l. in the western mountains of this region. But these forms are located at lower elevations (3,900-3,600 m a.s.l.) in the east of the mountains (Owen et al., 2006;Wang et al., 2013). The regional climate types are divided into three subzones: the arid feature and lower mean annual temperature in the Kunlun Mountains and the Chaka-Gonghe region, the semi-arid and moderate temperature in the Qilian and Qinghai Lake area, and the relatively humid and warm climate condition in southeastern Qinghai Province. The mean annual precipitation of the regions mentioned above not only depends on the atmospheric moisture but is also related to the topographic relief in the mountains (Yu et al., 2013).
At the present day, annual precipitation in the western region of the NETP is 100-200 mm, with the majority of it occurring during the summer season. The mean annual temperature is 2.0°C; the mean January temperature at the Wudaoliang meteorological station (4,780 m a.s.l.), which is 50 km to the south, is-8.4°C, and the mean July temperature is 12.1°C . Owing to the low annual mean temperature, permafrost features are widespread in the study area. Vegetation in the drainage basins consists of dry alpine steppe communities (Stauch et al., 2017). In the eastern regions, the annual mean temperature is-0.3°C, the highest monthly mean temperature is 10.9°C (July), and the lowest monthly mean temperature is -13.5°C (January). The annual mean precipitation is 300-400 mm. The annual mean evaporation is about 1,300 mm (Liu et al., 2011). Alpine steppe and alpine meadow are the main vegetation types in the eastern region.

Sample Collection
The coordinates and elevations of all sampling sections in this study are shown in Table 1, and all profile photos are shown in Figure 2, in which the sampling sediment layers and dating results are displayed. The sections HF4, HF8, HF17, HHC, and HR are located in the Golmud River drainage basin. The first three sections consist of decimeters-thick alluvial sands to coarse gravels with moderate sorting, sub-rounding, and horizontal bedding traits. The alluvial fans are highly dissected, and several bluffs are developed due to river erosion. The section heights are 3, 24, and 13 m, respectively. The latter two are characterized by weakly sorted fine gravels, with a grain diameter of 2-5 dm, and sand lenticles sandwiched between the gravel layers. The sections NMH1, ALK3, and ALK4 lie on the Nuomuhong River and the Alake Lake catchments in the eastern Kunlun Mountains. The section NMH1 is a ∼8-m-high alluvial stack and comprises sub-rounded and poorly sorted alluvial sand and gravel layers, and the upper and middle layers contain sparse boulders. The section ALK3 is 7 m thick and composed of weakly rounded and sorted gravels, which show horizontal layers, sandwiching irregular sand layers. The alluvial deposits are capped by a discontinuous, decimeter-to-meter-thick layer of loess sediment. The section ALK4 comprises a 2.5-m-thick alluvial toe deposit comprising rufous, tawny horizontally bedded clay, and silt layers in the upper half and sand to five gravel layers in the lower part. The sections CK2, GH1, and GH2 are on the piedmont alluvial fan cascades from Ela Mountain and the Qinghainan Mountains. The section CK2 is a 10-m-thick alluvial stock and sandwiches half-meter sand lenticles. Inversely, the sections GH1 and GH2 are predominantly alluvial gravel deposits with 10-dm-thick sand lenticles, and have a thickness of 6.0 and 2.5 m. The sections RYS3, RYS4, and XS1 are situated in the piedmonts of Riyue Mountain. Sections RYS3 and RYS4 mainly comprise crudely stratified meter and decimeter beds of fanglomerates and debris flows (section RYS4 contains abundant boulders), with a section thickness of 3.5 and 3.0 m. Section XS1 is a set of alluvial deposits that contains some boulders, with a thickness of 15 m. In the easternmost segment of the NETP, there are four sections, HL1, TD1, TD2, and GN1, which merge with the distinctive alluvial deposits, with stratification, interbedding, and clear sand and gravels, also containing a few sand lenticles. The four sections have a section thickness of 25, 20, 15, and 18 m, respectively. A total of 26 samples were collected for OSL dating with the stratigraphic order from the alluvial sections mentioned above.

Optically Stimulated Luminescence Dating
The OSL samples were collected by driving iron tubes (25 cm in length and 5 cm in diameter) into newly cleaned vertical sections. Each tube was covered with black plastic bags and was sealed with black packaging tape to avoid light exposure and moisture losses. In the laboratory, 3-cm sections of the sediments at each end of the cylinders were scraped away and were used for water content and dose rate measurements. The unexposed part in the middle of the tube was used for equivalent dose (D e ) determination. We followed the procedures of Lai (2010) to chemically pretreat OSL samples to remove feldspar contamination that might lead to age underestimation . We checked the purity of the retrieved quartz by infrared light stimulation (830 nm), and no obvious IRSL was observed. The pure quartz samples were mounted on the center (0.7 cm in diameter) of stainless steel discs with a diameter of 0.97 cm with silicone oil.
We used an automated Risø TL/OSL-DA-20 reader for OSL dating, which was equipped with a 90 Sr/ 90 Y beta source for OSL measurements at the Qinghai Institute of Salt Lakes, Chinese Academy of Sciences. The procedure involved stimulation by blue LEDs at 130°C for 40 s and signal detection using 7.5-mm-thick U-340 filters. We used neutron activation analysis to determine the U, Th, and K concentration at the China Institute of Atomic Energy in Beijing. For grains in the range of 38-63 μm, the α efficiency was taken as 0.035 ± 0.003 . The cosmic ray dose rate was estimated for each sample as a function of depth, altitude, and geomagnetic latitude (Prescott and Hutton, Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 702340 5 1994). Water content was calculated from moisture mass/dry mass (Aitken, 1985), and most of the measured results were near 0%, which is attributed to the arid modern environment. However, considering the uncertainties associated with burial and sediment types, values of 10-15% were estimated for the water content and dose rate calculations ( Table 1).
The equivalent dose (D e ) was determined using a combination of the single-aliquot regenerative dose (SAR) Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 702340 6 protocol Murray and Wintle (2000) and a standardized growth curve (SGC) method (Lai and Wintle, 2006), that is, the SAR-SGC method was applied. For each sample, six aliquots were measured using the SAR to construct the SGC, and then twelve aliquots were measured using the SGC. The final D e for a sample was given by the average of all eighteen D e values. Figure 3 shows the typical OSL decay curves and growth curves for ALK4-1 and GH1-1, respectively, suggesting that the OSL signals were from the fast component. Preheat plateau tests were conducted on samples ALK4-3 and RYS3-2, and a preheat plateau was clearly identified from 240 to 280°C in both samples (Figure 4). Based on the results, a preheat of 260°C for 10 s and a second heat of 220°C for 10 s were used for D e measurements. Signals of the first 0.64 s of stimulation were integrated for growth curve construction after back ground subtraction.
Dose recovery tests (Murray and Wintle, 2003), which can validate the SAR procedure, were performed on samples ALK4-3 and RYS3-2, and six aliquots of each sample were tested. The ratios of the average measured dose (54.4 and 106.6 Gy) to the given dose (50 and 105 Gy), which fall into the acceptable range of 0.9-1.1, suggested that the SAR protocol was suitable for D e determination in this study (Figure 4). Recuperation was calculated by comparing the sensitivity-corrected OSL signal of 0 Gy to the sensitivity-corrected natural signal to check the thermo-transferred signals. In this study, recuperations for different preheat temperatures were less than 5% of the natural signal, which was negligible (Figure 4). The "recycling ratio" was introduced to check for sensitivity change correction Murray and Wintle (2000), and for most aliquots, the recycling ratios fell into the acceptable range of 0.9-1.1 (Figure 4).

RESULTS
Considerable success has been obtained from applying OSL dating to dammed fluvial-lacustrine sediments in the central southern TP, with well-bleached samples that yielded accurate age estimates Chen et al., 2016;Liu et al., 2018;Wang et al., 2019). Consistent ages of eolian sand lenticels and adjacent alluvial sedimentary layers in the alluvial delta fan supported the reliability of OSL dating of alluvial sediments in the eastern Kunlun Mountains (An et al., 2018a). A modern fluvial sediment sample from the Kunlun River catchment (proximity to the HL3), which produced a D e of 1.20 ± 0.20 Gy (0.4 ± 0.1 ka) at a depth of 0.3 m, indicated a negligible residual dose in samples from the drainage basin .
The OSL ages determined in the current study are presented in Table 1, and the quoted chronology data from previous studies are listed in Supplementary Table S1. The entire data on alluvial fans cover from the most western part of the eastern Kunlun Mountains to the most eastern end of Qinghai Lake and Riyue Mountain, occupying the main monsoon-dominated region on the NETP. The initiating of the alluvial fans involved in this investigation occurred since the age of 99.2 ± 13.1 ka. Inferences suggested that there are older alluvial sediments than these. However, the older geomorphic processes could have greater odds of tectonic influence, especially the surface uplift; thus, we mainly focused on the alluvial formation since the Last Interglaciation. The youngest ages for a debris flow that is covered by loess have a date of 8.1 ± 0.6 ka. Analysis of the entire forming ages of the alluvial fans on the NETP suggests that the Last Glaciation and Deglaciation were the main developing stages, and the alluvial fans of different scales filled diffusely in the drainage basins of this region. The MIS 3 period is another relative activity stage of alluvial geomorphic processes. From 83.9.1 ± 9.36 to 99.2 ± 13.1 ka, the strong hydrological processes promoted the development of the alluvial fans. This stage corresponded to the late period of the MIS 5. Analyses of these dating results from regionalism showed that the alluvial Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 702340 fans in the drainage basins of the western high mountains yielded the broad-spectrum ages ranging from the Last Interglaciation to early Holocene periods. Inversely, in the eastern lower mountains, the age data were mainly concentrated in the MIS 3 period, and other climate stage ages were relatively scarce. These results could be ascribed to the high mountain topography (altitude variation from ∼4,000 to 6,000 a.s.l), which was more sensitive to the moisture fluctuation dominated by the hemispherical atmospheric circulations.

The Climatic/Tectonic Influences on the Alluvial Fans Since 200 Ka
The valley fills or the piedmont cascade of alluvial detritus in the arid and semi-arid drainage basins are commonly ascribed to climate change or tectonic uplift (Chen et al., 2011;Zhang et al., 2014;An et al., 2018a;Tao et al., 2020). In recent decades, various geomorphic units in the northeastern margin of the NETP, such as alluvial fans, fluvial terraces, high-stand lacustrine sediments, sand dunes, and moraines, have been intensively studied. Almost all previous studies have focused on the effect of climate changes for the surface geomorphic dynamics since the penultimate Interglaciation, although the influence of tectonic uplift could not be ruled out in some local places (for the specific statistics of the references, see Table 2). As shown in Table 2, the paleoclimate change is one of the major factors controlling depositional processes, and it played a prominent role in the formation of various geomorphic units on the NETP. Most landform processes in the Golmud River of the eastern Kunlun Mountains, including impressive piedmont alluvial fans and delta alluvial fans, stretch into the arid basin and are related to climate change rather than surface uplift (Owen et al., 2006;An et al., 2018a, b). The formation of the alluvial fan deposits is related to the hydrological variation in the mountains at the southern border of the Heihai Lake, which is located in the headwater area of the Golmud River, and it influenced the variation of the lake level through blockage of the lake outflow (Stauch et al., 2017). The synthetic wetness index from various geomorphic units in the mountain valleys of the eastern Kunlun Mountains showed a good correspondence to moisture records that were controlled by the ISM, showing their paleoclimatic cause (An et al., 2018a). In a study of the 10 Be flux of lacustrine sediments collected at Kunlun Pass, Chen et al. (2021) suggested that there was a shift in the dominant influence in the region from monsoons to westerlies after 4 ka BP. In the Qinghai Lake watershed, alluvial sediments mainly deposited during the middle-to-late MIS 3, and they had a paleoclimatic origin . A stream-dominated fan in the Xining Basin, which lies on the margin of the NETP, was reported by Gao et al. (2018), Their results indicated that a climate-dominated transition of the sedimentary environment from high-to low-energy fans occurred from the penultimate glacial to the last interglacial. A database of the terrace ages of Quaternary fluvial terraces on the TP and surrounding orogenic belts was integrated by Tao et al. (2020), and they grouped terraces into the three regional climate systems, the ISM, EASM, and westerlies, and found that the formation of terraces younger than 200 ka was mainly controlled by climate changes. However, the incision of fluvial terraces in the Golmud River drainage basin was ascribed to the transitory surface uplift by Chen et al. (2011). They argued that the valley fills occurred during the tectonic stable stage of the Last Glaciation. From the analysis mentioned above, we can draw a Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 702340 conclusion that the alluvial fans developed in the studied area responded to the paleoclimatic fluctuations rather than the tectonic uplifts during the late Quaternary.

The Glacial Forcing Mechanism on Alluvial Fans' Development
Enhanced evidence suggested that the dynamic of alluvial fans on the NETP showed a periglacial origin, especially the high relief and elevation regions (Van Der Woerd et al., 2002;Owen et al., 2006;Stauch et al., 2017;An et al., 2018a;Cui et al., 2020). The glacier activities in the high mountains produced abundant glacial debris in the piedmont during glaciation periods (Owen et al., 2006). Subsequently, the glacial debris cascaded onto alluvial fans and was transported into tributary and main valleys by abundant glacier meltwater and finally completed valley filling (An et al., 2018b). The process of aggradations from the piedmont to the alluvial delta fan was probably multistage or leapt forward to the downstream drainage basin, as there are several reservoir segments in the main valley through which detrital materials are transported step by step (Wang et al., 2009;An et al., 2018b).
Therefore, the multistage alluvial fans or valley fills in drainage basins and the margin area of arid basins are the main linkage between the orographic glaciers and the distal lower basins. The dating age for the moraine sediments is limited because there are few suitable chronology methods or dating materials (Cui et al., 2020). Thus, the dynamics of the alluvial fans reflect environmental/climatic information of the provenance area, namely, glacial activity processes. Both the geomorphic dynamics showed more consistent feedback to each other on the orbital timescale. In view of the foregoing, we compared the developing ages of alluvial fans in different areas on the NETP with those of the glacial advances from the eastern Kunlun Mountains, the Dalijia Peak of eastern Riyue Mountain, the Queer Shan ranges of the eastern TP, and Yulong Mountain of the southeastern TP (Owen et al., 2006;Wang et al., 2013;Ou et al., 2014;He et al., 2016) (Figures 5A,B). There is broad consistency in the timing of both geomorphic dynamics, especially in the early Holocene, Last Deglaciation, and MIS 2. The stages MIS 3 and late MIS 5 were also relatively active periods for the alluvial fans and glaciers. Even in the records of both geomorphic processes, the Younger Dryas (YD) event can be distinguished during ∼10.5-12.0 ka ( Figures 5A,B). The good correspondence between forming ages of alluvial fans and those of glacial advances indicates that the latter controls the former's extension. Meanwhile, strong hydrographic discharge also probably promoted the stretches of alluvial fans and finished most valley fills of drainage basins; the hypothesis can be tested by the mass development of alluvial fans during the Last Deglaciation/Last Interglaciation periods. The melting glaciers yielded abundant high-load runoffs and transported surplus sediments, which boosted the aggradations toward the fan toes (Ou et al., 2014). The strengthened precipitation together with moderate glacier activities during the Late Interglaciation was probably somewhat responsible for the development of the alluvial fans, corresponding to the slight enhancement for the current probability density function (PDF) curves ( Figures  5A,B). Previous studies provided enhanced evidence that the glacial advances occurred during the penultimate glacial cycle, MIS 3, Last Glacial Maximum (LGM), Late Glacial, and Holocene periods in the northeastern or southeastern mountain regions of the TP (Owen et al., 2003(Owen et al., , 2006Ou et al., 2014;Wang et al., 2014;He et al., 2016;Cui et al., 2021). The expansion of the extent of glaciers was affected by the change in precipitation, and the extent of glaciers during the MIS 3 was more extensive than that during the LGM (Owen et al., 2003;Cui et al., 2021). The alluvial fans originated from orographic outwash formed on the southern side of the Kunlun Mountains prior to 200 ka (Owen et al., 2006) and  Alivernini et al., 2018). (D) δ 18 O record of the stalagmites from the Xiaobailong cave, southwestern TP (Cai et al., 2015). (E) Arabian Sea summer monsoon stack (Caley et al., 2011;Bolton et al., 2013). (F) Qinghai Lake westerlies climate index . The gray frames represent the prominent increase stages of alluvial aggradations. MIS, Marne Isotope Stage; LDG, Last Deglaciation; YD, Younger Dryas event.
Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 702340 later frequently occurred in different paleoclimatic stages (Xu et al., 2010). Meanwhile, alluvial sequences in the eastern part of the NETP, such as Riyue Mountain and the Qinghainan Mountains, intensively formed during the MIS 3. This phenomenon is probably ascribed to the lower temperature and more abundant precipitation in this stage, including glacier activities combining the sufficient rainfall that were responsible for the alluvial geomorphic evolutions. Cui et al. (2021) argued that the coupling effect of lower temperature and increased precipitation on the NETP resulted in the expansion of the orographic glaciers during the MIS 3 rather than the LGM. A simulation result of the CAM4 climate model suggests that an expansion of glaciers would lead to a decrease in surface air temperatures during both the summer and winter over most regions of the TP, and coupled with increased summer precipitation and increase in cold glacier meltwater supply, it would be in favor of alluvial advance and high lake levels (Liu et al., 2020).

Paleoclimatic Implications for Glacier-Induced Alluvial Fans
The NETP is situated in a sensitive arid and semi-arid zone between the ISM-controlled and the westerly-influenced regions of Asia, and the dominated atmospheric circulations were commonly altered from the interglacial to glacial periods . Thus, construction of the paleoclimate changes from the dynamics of the glacier and alluvial fan geomorphic processes is an intricate procedure. But, this allows us to distinguish the different climatic systems on the interglacial-to-glacial or interstadialto-stadial cycle timescales. The formation of the alluvial fans involved in the current study began since the late MIS 5 and ended during the early Holocene (Table 2 and Figure 5A). These results referred to various climatic stages, including the MIS 5, MIS 3, MIS 2, Last Deglaciation, and YD periods. Ramisch et al. (2016) suggested that the moisture status of the eastern Kunlun Mountain region is mainly dominated by the ISM and the mountains form the northern boundary of the monsoon system, which triggered orographic rainfall or snowfall. Cui et al. (2020) argued that the temperature drops ranged from -0.5 to 0.1°C with precipitation being 140 ∼ 200% during the MIS 3 in the Kunlun Mountains. The lower temperature and high moisture were in favor of glacial advances (Cui et al., 2021), which contributed to abundant outwash supply and alluvial fill in the drainage basin. Recent studies have indicated that the glaciations of most regions of the TP were influenced by the ISM or the Asian monsoon during the MIS 5 and MIS 3, and the increased south Asian monsoon supplied more precipitation as snow at high altitudes (Owen et al., 2003;Xu et al., 2010;Stauch et al., 2017;Cui et al., 2021). However, a remarkable feature of the age distribution for alluvial fans indicates that the periods of Last Glaciation and Deglaciation are the main developing stages for the alluvial aggradations. The results are well consistent with the ages of the glacial advance ( Figures  5A,B). An et al. (2012) found an antiphase relationship between the westerlies and the ASM for glacial millennial Frontiers in Earth Science | www.frontiersin.org July 2021 | Volume 9 | Article 702340 timescales in the Qinghai Lake catchment during the Last Glaciation. Li et al. (2020) argued that the Asian monsoon and the westerlies showed positive correlation on millennial timescales but anticorrelation during suborbital cycles. Therefore, we synthesized a record comparison between our alluvial age data with the proxies for ISM-dominated regions and the westerlies index from Qinghai Lake ( Figures  5C-F). As shown in the figure, the alluvial age record closely matches that of glacier activities and the record of high lake levels on the southern TP, such as those of the Taro Co lake group ( Figure 5C). Alluvial aggradations are also consistent with the δ 18 O record of the stalagmites from the southwestern TP and the summer monsoon stack from the Arabian Sea ( Figures 5D,E), suggesting the influence of the ISM during the MIS 5, MIS 3, and early Holocene. But the prominent increase in the alluvial ages during the Last Glaciation and Deglaciation showed more consistency with the westerlies index. Thus, we considered that strengthened westerlies combined with slight monsoon dominated the intensive glacial activities and alluvial aggradations in the drainage basin on the NETP during the Last Glaciation. This glacial expansion stage was verified using a simulation. The decrease in surface air temperatures and increase in precipitation rates during both the summer and winter boosted the orographic glacier activities (Figure 6), and they also showed inter-feedback to each other, especially the western mountains of the region. Subsequently, massive glaciers melted with the warming climate, finally facilitating a large scale of development for the alluvial fans during the Last Deglaciation.

CONCLUSION
Investigation on the climatically induced alluvial fans that cascaded from mountains in the semi-arid regions provides insight into the geomorphic and hydrographic processes driven by paleoclimatic fluctuations. We investigated a series of alluvial fans on the NETP region using the morphostratigraphic study and OSL dating. These alluvial sequences cascaded from the high piedmonts in the eastern Kunlun Mountains, Ela Mountain, the Qinghainan Mountains, and Riyue Mountain, ultimately filling the drainage basins or stretching into the arid sedimentary basins. Integrated analysis indicated that the alluvial geomorphic processes were related to the paleoclimate changes rather than the surface uplifts. These are basically composed of the coarse-to-fine, weak, and angular gravels and sandy lenticles, which occasionally contained boulders in specific layers. The OSL dating results showed that the alluvial fans had mainly developed during the late MIS 5, MIS3, and Last Glaciation and Deglaciation, little occurring during the early Holocene. The findings are firmly consistent with the orographic glacial advance periods in the northeastern and eastern regions of the TP. Considering that the alluvial fan is a production of downward outwash transportation according to previous studies, we argued that the dynamics of the alluvial fans and alluvial valley fills should be ascribed to the paleo-glacier activities during different climatic stages since the late Pleistocene. The comparisons of wet indexes between our age data and records from the ISMand westerly-dominated regions showed that both the climatic systems had coupling impacts on the glacial advances and subsequent alluvial aggradations on the supply of moisture in different paleoclimatic stages. The simulation also suggested that lower temperature and enhanced precipitation during the LGM resulted in the glacier expansion. This result matches with our data that indicate that the alluvial aggradations mainly occurred during the MIS 2 and the sequent postglacial stage.

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.

ETHICS STATEMENT
Written informed consent was obtained from the relevant individual(s) for the publication of any potentially identifiable images or data included in this article.

AUTHOR CONTRIBUTIONS
FA and SL designed the research. FA, BN, DG, and JZ performed the research. LC and XC made the figures and analyzed the data. FA and BQ wrote the manuscript. TC analyzed the OSL dating results with age models.

FUNDING
This study was supported by the Natural Science Foundation of Qinghai Province (Grant no. 2020-ZJ-727) and the China NSF (Grant no. 41961014).