Abstract

Sedimentologic and provenance analyses for the Qaidam Basin in the northern Tibetan Plateau help to elucidate the stratigraphic signatures of initial deformation and exhumation in basin-bounding ranges. The basin recorded sedimentary transitions in response to uplift and unroofing of several distinctive source regions. Along the NE basin margin, a detrital record of exhumation and basin isolation is preserved in the 6200-m-thick Cenozoic succession at the Dahonggou anticline. An up-section shift from axial fluvial and marginal lacustrine deposition to transverse fluvial sedimentation suggests progradation and increasingly proximal sediment sources, reflecting activation and advance of crustal deformation. Provenance results from sandstone petrology, U-Pb geochronology, and heavy mineral analyses indicate initial late Paleocene–early Eocene derivation from igneous, metamorphic, and sedimentary sources, consistent with Permian–Triassic arc rocks dominating the southern (Kunlun Shan) or southwestern (Qimen Tagh) basin margins. Up-section variations in sediment composition and detrital zircon U-Pb age distributions are attributed to Eocene–Oligocene derivation from lower Paleozoic and Mesozoic igneous and metamorphic rocks of the central to northern Qilian Shan–Nan Shan. Disappearance of igneous sources and persistence of metamorphic sources are consistent with derivation from the southern Qilian Shan–Nan Shan during early–middle Miocene shortening along the frontal Nan Shan–North Qaidam thrust belt. These results are supported by paleocurrent analyses revealing an Eocene shift from roughly E-directed (axial) to SW-directed (transverse) dispersal of sediment. Variations in lithofacies, composition, U-Pb ages, and paleoflow are consistent with late Paleocene–early Eocene exhumation in the Kunlun Shan followed by middle Eocene–middle Miocene exhumation in the Qilian Shan–Nan Shan. The up-section disappearance and reappearance of diagnostic U-Pb age populations can be associated with progressive unroofing of multiple thrust sheets, successive input of sedimentary and magmatic sources, and southward encroachment of Qilian Shan–Nan Shan shortening into the Qaidam Basin. The sedimentary record presented here indicates that during the Paleogene, the unified Qaidam-Tarim Basin was partitioned and uplifted as it was incorporated into the growing Tibetan Plateau. Comparison with basins on and surrounding the Tibetan Plateau suggests that basement strength and lateral homogeneity, and formation of syndepositional structural dams are among the primary controls on formation of giant sedimentary basins.

INTRODUCTION

The 5–14-km-thick Cenozoic succession in the Qaidam Basin (Fig. 1) provides an important detrital record of the intraplate response to the India-Asia collision. Growth of this internally drained basin at ∼2500 m elevation has been ascribed to northward advance of Neogene collisional deformation and accompanying flexural subsidence (Meyer et al., 1998; Métivier et al., 1998; Tapponnier et al., 2001). However, evidence for pre-Neogene shortening in the northern Tibetan Plateau contradicts the proposed systematic northward propagation of deformation. Specifically, rapid late Eocene–Oligocene sedimentation and tectonic rotation (Horton et al., 2004; Dupont-Nivet et al., 2004, 2008), exhumational cooling (Clark et al., 2010), and thrust faulting (Duvall et al., 2011) indicate early deformation in northeastern Tibet (see review by Yuan et al., 2013). These findings highlight the limited understanding of early Cenozoic deformation and plate-interior basin evolution within the Tibetan Plateau.

Contrasting reconstructions have been proposed for mountain ranges bounding the Qaidam Basin. To the south, the eastern Kunlun Shan records pre-Cenozoic deformation (Cheng et al., 2015), an early phase of Eocene shortening (Jolivet et al., 2003; Yin et al., 2007a; Duvall et al., 2013), and complex Oligocene to Quaternary transpressional deformation ranging from ca. 30–20 Ma in the central range (Mock et al., 1999; Wang et al., 2004; Yuan et al., 2006; Clark et al., 2010; Mao et al., 2014) to 12–8 Ma in the west and 8–5 Ma in the east (Duvall et al., 2013). Left-lateral slip on the Kunlun fault initiated by middle–late Miocene time (Kirby et al., 2007; Fu and Awata, 2007; Duvall et al., 2013). To the north, the Qilian Shan–Nan Shan may record early Cenozoic faulting and middle–late Cenozoic migration of deformation and sedimentation to the south and east (Yin et al., 2008a). In contrast, thermochronology studies in the northern Qilian Shan identify bedrock cooling at 12–6 Ma (Zheng et al., 2010; Craddock et al., 2011b).

Continuous sediment accumulation in the Qaidam Basin from the Paleogene to Holocene may record the onset of exhumation in basin-bounding ranges, and associated changes in sediment sources. However, the Qaidam Basin epitomizes the challenge posed by plate-interior basins, in which multiple sediment sources are characterized by similar lithologies, compositions, or thermal histories (e.g., Graham et al., 1993). Sediment provenance tracers are often consistent with multiple potential sources and do not point to systematic up-section changes in sediment sourcing (Rieser et al., 2005, 2006a, 2006b; Zhuang et al., 2011; Jian et al., 2013; Cheng et al., 2015). However, application of multiple provenance techniques has identified changes in sediment sources (e.g., Hanson, 1999).

This study addresses the evolution of the Qaidam Basin by integrating multiple sedimentologic and provenance records from the northeastern margin of the basin to evaluate Cenozoic exhumation patterns in major basin-bounding ranges. Uniting stratigraphic and sedimentologic relationships with multiple provenance techniques, including petrographic, U-Pb geochronologic, heavy mineral, and paleocurrent analyses, sheds light on the growth of flanking ranges and associated evolution of Qaidam depositional systems. Interpretations of detrital zircon U-Pb results are guided by multidimensional scaling, a visualization technique representing statistical similarity among samples. Collectively, these data sets identify the early intraplate response to India-Asia collision through identification of nonuniform exhumation along Qaidam Basin margins.

GEOLOGIC SETTING

The Qaidam Basin is surrounded by diverse igneous and metamorphic rocks constituting Precambrian basement, Paleozoic–Mesozoic magmatic arc, and accretionary prism assemblages overlapped by Mesozoic–Cenozoic basin fill. Flanking ranges (Figs. 1 and 2) include the Qilian Shan–Nan Shan system to the north, Kunlun Shan to the south, Altyn Shan to the west, and Ela Shan and Dulan-Chaka highland to the east.

The basin is commonly modeled as a rigid block composed of strong lithosphere and continental basement (Clark and Royden, 2000; Braitenberg et al., 2003; Dayem et al., 2009). Crystalline basement in northern Tibet is composed of the Qilian and Qaidam terranes: Proterozoic and Paleozoic assemblages accreted to the Tarim and Sino-Korean cratons during Paleozoic–early Mesozoic closure of the Tethys Ocean (Hsü et al., 1995; Sobel and Arnaud, 1999; Yin and Harrison, 2000; Yin et al., 2007b; Pullen et al., 2008; Gehrels et al., 2011).

To the north-northeast, the Qilian Shan–Nan Shan (Figs. 1 and 2) forms a major NW-trending fold-thrust belt overprinting a system originally formed during early Paleozoic collision and arc accretion and later reactivated during Mesozoic–Cenozoic extension, contraction, and strike-slip deformation (Ritts and Biffi, 2001; Jolivet et al., 2001; Xiao et al., 2009; Yan et al., 2010; Zheng et al., 2010). Synorogenic deposits in the adjacent Hexi Corridor document the onset of crustal shortening in the northern Qilian Shan by the Miocene (Wang and Coward, 1990; George et al., 2001; Bovet et al., 2009).

To the south, the E-trending Eastern Kunlun Shan (Figs. 1 and 2) is a transpressional fold-thrust system overprinting a Permian–Triassic magmatic arc (Dewey et al., 1988; Mock et al., 1999; Jolivet et al., 2001, 2003; Xiao et al., 2002; Mo et al., 2007; Yuan et al., 2009; Clark et al., 2010) and early Paleozoic arc in the westernmost Qimen Tagh range (W. Li et al., 2013). Left-lateral displacement along the Kunlun fault is thought to have initiated with late Miocene extension in southern Tibet and continues to the present (Van der Woerd et al., 1998; Lin et al., 2002; Fu and Awata, 2007; Kirby et al., 2007; Duvall et al., 2013).

To the northwest, the NE-striking Altyn Tagh fault (Figs. 1 and 2) persists for 1600 km, with left-lateral displacement absorbed by NE-SW shortening in the Qilian Shan–Nan Shan. Altyn Tagh deformation initiated between 49 and 30 Ma, with strike-slip motion dominant by Miocene time, and progressively decreasing strike-slip motion (Jolivet et al., 2001; Yin et al., 2002; Wang et al., 2005) related to a shift from plate-like lateral extrusion to distributed crustal shortening (Yue et al., 2004; Ritts et al., 2008; Yin et al., 2008a, 2008b).

To the east, the Qaidam Basin is flanked by the Ela Shan range and Dulan-Chaka highland (Fig. 2), which form the boundary between Qaidam and externally drained Neogene basins of northeastern Tibet. These include the Miocene–Pliocene Gonghe, Xunhua, and Chaka basins, which occupy the Yellow River headwaters (Craddock et al., 2011b; Lease et al., 2011; Lu et al., 2012; H.-P. Zhang et al., 2012).

Ages of the Paleocene–Quaternary succession in the Qaidam Basin, based on nonmarine fossil assemblages and magnetostratigraphy (Fig. 3), have been extrapolated across the basin on the basis of seismic reflectors (QBGMR, 1991; Yang et al., 1992; Sun et al., 2005; Rieser et al., 2005, 2006a; Zhang, 2006; Fang et al., 2007; Wang et al., 2007; Gao et al., 2009; Lu and Xiong, 2009). However, ages from critical sedimentary sections are often only broadly constrained due to a lack of key marker beds. This study presents sedimentologic and multiproxy provenance results for the northern basin margin (Figs. 2 and 4), complementing magnetostratigraphic, paleoclimatologic, and paleontological data from the Dahonggou and Huaitoutala anticlines (Fang et al., 2007; Wang et al., 2007; Lu and Xiong, 2009; Zhuang et al., 2011; Song et al., 2013a, 2013b). Ages used here are consistent with the regional study conducted by Zhuang et al. (2011) and published magnetostratigraphic results from the two sections examined here (Fang et al., 2007; Lu and Xiong, 2009). A 6200-m-thick stratigraphic section was measured on the southern limb of the Dahonggou anticline west of Xitieshanzhen (37.5486°N, 95.1655°E). Continuous exposures include the upper Paleocene–Eocene Lulehe (∼1100 m), Eocene Xiaganchaigou (∼2400 m), Oligocene Shangganchaigou (∼1300 m), lower Miocene Xiayoushashan (∼800 m), middle–upper Miocene Shangyoushashan (∼300 m), upper Miocene–Pliocene Shizigou (300 m), and overlapping Quaternary Qigequan Formations (Fig. 3). Although the lower 3500 m section of the Cenozoic succession remains poorly dated, magnetostratigraphic and sedimentologic studies of neighboring sections have constrained the age of the upper part of the section to between 34 Ma and 8 Ma and have provided documentation of general lithologic trends of upper stratigraphic levels (Lu and Xiong, 2009; Zhuang et al., 2011).

POTENTIAL SEDIMENT SOURCES

Potential sediment sources include mountain ranges surrounding the Qaidam Basin on all flanks (Fig. 2). Buried sedimentary and metamorphic units form the Qaidam basement, which could also have provided sediment during early basin filling. Next, we highlight the characteristic crystallization ages, lithologies, and heavy mineral constituents of each of these sources.

Southern Sources

Along the southern basin margin, the Eastern Kunlun Shan east of ∼92°E (Figs. 1 and 2) is principally composed of Permian–Middle Triassic (260–210 Ma) plutons with isolated lower Paleozoic granitic rocks (450–380 Ma; Fig. 5C; Harris et al., 1988; Roger et al., 2003; Mo et al., 2007; Pullen et al., 2008; Yuan et al., 2009; McRivette, 2011; J.Y. Zhang et al., 2012) intruding late Proterozoic to Paleozoic metamorphic basement rock. Heavy minerals associated with magmatic and metamorphic basement (pyroxene, hornblende, epidote, garnet) are expected from the Kunlun Shan and Qimen Tagh (L. Li et al., 2015). In contrast to the East Kunlun Shan, early Paleozoic (513–420 Ma) subduction-related intrusions occupy the Qimen Tagh (W. Li et al., 2013). Farther south, gneissic basement rocks north of the Kunlun fault yield Paleoproterozoic ages (Zhang and Zheng, 1994). Upper Carboniferous to Lower Permian shallow-marine sedimentary and volcanic strata are unconformably overlain by Jurassic–Cenozoic nonmarine strata (Mock et al., 1999). Southeast of the Kunlun Shan, there is the Songpan-Ganzi Basin (Fig. 1), a belt of Middle–Upper Triassic deep-marine clastic deposits with a broad mixture of Permian–Triassic, middle Paleozoic, and rare Precambrian detrital zircon age populations (Fig. 5A) from the North China block, South China block, Kunlun, and Qiangtang terranes (Weislogel et al., 2006, 2010; Weislogel, 2008; Enkelmann et al., 2007; Ding et al., 2013).

Directly south of the Eastern Kunlun Shan, the Hoh Xil Basin contains Cenozoic siliciclastic rocks characterized by Neoproterozoic, Paleozoic, and Permian–Triassic sources (Fig. 5B) derived from the Lhasa and Qiangtang terranes, as well as the Eastern Kunlun Shan (Dai et al., 2012; Staisch et al., 2014).

Eastern Sources

To the east, Proterozoic metamorphic basement, Paleozoic–Mesozoic sedimentary and metasedimentary units, and Permian–Triassic granitoids are exposed in the Ela Shan and Dulan-Chaka highland (Fig. 2). These ranges divide the Qaidam Basin from Cenozoic intermontane basins of northeastern Tibet. Basement is composed of early Proterozoic high-grade metamorphic rocks, Paleozoic–Mesozoic sedimentary and low-grade metamorphic rocks, and extensive Permian–Triassic igneous bodies (Fig. 5G; Liu et al., 2004; Wu et al., 2004; Chen et al., 2012; Duvall et al., 2013).

Western Sources

The Altyn Shan to the west of Qaidam Basin includes Archean–lower Paleozoic magmatic arc, accretionary prism, and related sedimentary successions. The Altyn Tagh fault (Figs. 1 and 2), which bisects the Altyn Shan, is a 1600-km-long left-lateral fault. Sedimentary units exposed in the Altyn Shan include Neoproterozoic quartz arenite and marbleized stromatolitic limestone (Gehrels et al., 2003a). Paleoproterozoic (1.6–2.0 Ga), Neoproterozoic (922–928 Ma), Ordovician–Silurian (430–480 Ma), and Permian (260–280 Ma) zircon U-Pb ages are common in igneous and metamorphic basement of the Altyn Shan (Fig. 5F; Gehrels et al., 2003b).

Northern Sources

Finally, north of the basin, localized ultrahigh-pressure metamorphic complexes, magmatic arcs, and recycled sedimentary units comprise the northern Qaidam and Qilian Shan–Nan Shan region (QBGMR, 1991). These sedimentary and magmatic units are underlain by Proterozoic and Archean crystalline basement (Yin and Harrison, 2000; Gehrels et al., 2003a, 2003b). The Qilian Shan–Nan Shan (Figs. 1, 2, and 5D) includes the North Qaidam terrane of Proterozoic gneiss overlain by Paleozoic metasedimentary rocks (QBGMR, 1991). Early Paleozoic arcs developed along the southern margin of the North China craton (Yin and Harrison, 2000; Xiao et al., 2009), including marine sedimentary rocks, ophiolites, and ultrahigh-pressure metamorphic complexes (Menold et al., 2009; Xiao et al., 2009). Heavy mineral provinces of the southern Qilian Shan–Nan Shan can be expected to contain high volumes of epidote, hornblende, garnet, rutile, and metamorphic accessory minerals such as kyanite, andalusite, and sillimanite. The Qilian Shan–Nan Shan contains isolated Permian–Triassic plutons (Chen et al., 2012). Silurian–Ordovician strata are dominated by 450–490 Ma U-Pb ages from a local volcanic source and a minor Proterozoic component (Gehrels et al., 2011). In addition, recycling of Mesozoic and late Paleozoic sedimentary units would provide Permian–Triassic-aged zircons, such as the unimodal 270–280 Ma population observed in Jurassic sandstones associated with magmatic ages in the Kunlun arc (Yan et al., 2010; Gehrels et al., 2011).

Local Basement Sources

Neoproterozoic and lower Paleozoic metamorphic rocks, flysch, and carbonates form part of the Qaidam basement beneath Mesozoic–Cenozoic basin fill (Wang and Coward, 1990; Xia et al., 2001; Chen et al., 2012). Jurassic–Cretaceous nonmarine rocks exposed on the northern Qaidam margins and in well logs (QBGMR, 1991) were sourced from early exhumation in the Qilian and Altyn Shan (Ritts and Biffi, 2001; Sobel et al., 2001; Gehrels et al., 2003b).

SEDIMENTOLOGY AND DEPOSITIONAL ENVIRONMENTS

We recognize major depositional shifts based on facies distributions of six formations (Fig. 6). Initial coarsening (Lulehe) was followed by fining (Xiaganchaigou), continued fining and transition from fluvial to marginal lacustrine deposits (Shangganchaigou), and final coarsening (Xiayoushashan, Shangyoushashan, and Shizigou).

Lulehe Formation

Description

Fine-grained, distinctively red deposits of the lower Lulehe Formation (0–220 m) are characterized by trough cross-stratified and horizontally bedded sandstone (Fig. 7A) with laminated mudstone lenses and pebble lags above sharp bases (Fig. 8F). Fine- to medium-grained sandstone packages contain downstream accretion surfaces and are typically 3–7 m thick and 10–20 m wide with lenticular geometries. These fine upward into mudstone intervals up to 3 m thick, which feature mottling and variable destruction of primary sedimentary structures.

The upper Lulehe Formation (220–1100 m) contains 1–5-m-thick lenticular beds of polymict, poorly sorted, clast-supported, granule-cobble conglomerate interbedded with medium- to coarse-grained, trough cross-stratified sandstone (Figs. 7B and 8A). Conglomerate beds have sharp erosional surfaces incising up to 1 m into underlying beds and a matrix of poorly sorted, calcite-cemented sandstone. These beds pinch out laterally, fine upward to coarse-grained sandstone, and include horizontally stratified, imbricated, and structureless varieties. Overlying conglomerate beds are horizontally laminated or structureless red mudstone beds 2–5 m thick (Fig. 8H).

Interpretation

Lower Lulehe levels represent a moderate-energy fluvial system with stable floodplains. The sharp bases, gradational upper surfaces, and low width-to-thickness ratios of the sandstones are similar to sand bodies typical of channel-fill deposits (e.g., Friend et al., 1979; Tornqvist, 1993). Structureless red mudstones are interpreted as overbank deposits with subaerial exposure and pedogenesis. Trough cross-stratified sandstone bodies, stable overbank deposits, and lack of laterally accreting bar forms are consistent with a low-sinuosity channel belt in a sandy, low- to moderate-gradient braided or anastomosing fluvial system with multiple channels, downstream accreting bar forms, and vertical aggradation (Cant and Walker, 1978; Ethridge et al., 1999; Makaske, 2001). We interpret the increase in grain size and sand-to-mud ratio in the upper 20 m of this interval as a transition to a higher-competence fluvial system.

Upper Lulehe conglomerate beds preserve channel thicknesses 1–5 m thick (Fig. 7B). Clast support, imbrication, poor sorting, coarse grain size, and scoured bases indicate deposition by a higher-energy system (Nemec and Steel, 1984; Smith, 1990). Structureless and horizontally stratified conglomerate beds are considered components of downstream-migrating compound bars (Miall, 1977). This interval is consistent with a transition to a higher-gradient gravelly braided fluvial system, although the presence of mudstone beds indicates highly variable flows.

Xiaganchaigou Formation

Description

The intermediate succession (1100–3550 m) contains finer material with few gravel clasts (Fig. 7C). The deposits consist of multistory sand bodies in 5–10-m-thick packages with sharp bases, commonly fining upward in the basal 50 cm; these are composed of 5–10-cm-thick beds of horizontally and ripple cross-stratified, fine- to very fine-grained silty sandstone (Fig. 8B). The packages contain mudstone drapes and interbeds and are capped by 1–3-m-thick mudstone intervals displaying blocky soil aggregates (peds), pedogenic slickensides, carbonate nodules, red and green mottling, rhizoliths, and moderate bioturbation.

Above 1250 m, interbedded sandstone and mudrock deposits dominate, featuring upward-coarsening and upward-thickening intervals of laminated mudstone and massive to horizontally stratified, laterally continuous sandstone sheets 0.5–2 m thick (Fig. 8C). In sharp contact with these interbedded deposits, trough cross-stratified, medium-grained sandstone packages reach 20–40 m (up to 75 m) in thickness (Fig. 8I). These thick sandstone packages persist to ∼3550 m, with mudstone intervals increasing above the ∼3000 m level.

Interpretation

Upward-fining mudstone, cross-bedded sandstone, and interbedded paleosols represent a fine-grained fluvial system with pedogenically modified overbank settings. Trough cross-stratified sandstones indicate dune migration, with laminated mudstones resulting from low-flow isolation of dunes and mud deposition (Smith, 1990). Stability of the surrounding floodplain is suggested by ubiquitous root traces, mottling, and soil structures. The lack of distinctive soil profiles points to aggradational alluvial soils (Besly and Fielding, 1989). Thick sandstone bodies are typical of amalgamated channels in a high-sinuosity channel belt; upward-coarsening and upward-thickening mudstone and sandstone intervals represent crevasse splays. These overbank, channel margin, and in-channel deposits are considered representative of a meandering fluvial system (e.g., Nanson, 1980; Kraus, 1996).

Shangganchaigou Formation

Description

The Shangganchaigou (3550–4805 m) shows continued fining and is dominated by laminated mudstone and ripple-laminated fine-grained sandstone (Fig. 7D). Thin beds (0.5–2 m) of mottled, nodular mudstone coarsen upward into fine-grained structureless or cross-laminated sandstone (Fig. 8K). These thinly bedded heterolithic deposits are commonly incised by upward-fining, lenticular cross-laminated sandstone bodies. Cross-laminae are generally asymmetric, with occasional symmetrical cross-laminae (Fig. 8C). Rare lateral accretion surfaces in sandstone beds are less than 1 m thick.

Interpretation

The thin-bedded sandstone and laminated to structureless mudstone of this association are attributed to floodplain deposition, based on upward-coarsening packages and association with channel sandstones (Porter and Gallois, 2008). Mudstone packages capping channel sandstones indicate phases of channel abandonment. Lack of lacustrine ichnofacies, freshwater fish fauna (Wang et al., 2007), soft sediment deformation, and lacustrine carbonate or beach facies helps rule out an open lacustrine setting, though such facies are identified in the distal basin (Zhuang et al., 2011). These deposits likely represent an alluvial plain with splays into shallow ephemeral lakes, as expressed by bidirectional ripples and thick mudstone packages.

Xiayoushashan, Shangyoushashan, and Shizigou Formations

Description

Upper stratigraphic levels (4805–6200 m) comprise the second major coarse-grained interval (Fig. 7E). Deposits are dominated by normally graded, trough cross-stratified conglomerate and medium- to very coarse-grained sandstone (Fig. 8E). Beds are 1–4 m thick, with sharp bases and lenticular geometries (Fig. 8L). At 4800–5800 m, <1-m-thick beds of mottled, structureless siltstone are interbedded with fine- to very fine-grained, structureless to ripple cross-laminated sandstone. The frequency and thickness of granule-pebble conglomerate beds increase up section (4805–5787 m), with capping deposits (above 5787 m) dominated by clast-supported, planar to trough cross-stratified, imbricated conglomerate containing sandstone lenses and clasts up to 1 m in diameter.

Interpretation

Lithofacies and stratal geometries are characteristic of a sandy to gravelly braided fluvial system (Rust, 1972; Cant and Walker, 1978). Normally graded, trough cross-stratified, clast-supported conglomerates are consistent with dune migration in fluvial channels during flood and waning flow conditions (Miall, 1977). The absence of matrix-supported debris-flow conglomerates and sandy sheetflood deposits precludes a proximal alluvial-fan setting. Coarse grain size, cross-stratification, erosive bases, and lenticular bed geometries are representative of proximal deposition in an aggradational, low-sinuosity channel belt. The upward-thickening and upward-coarsening pattern likely indicates a progradational system in which coarse-grained transverse fluvial systems advanced into the basin.

Regional Evolution of Depositional Systems

Shifting depositional systems are evidenced by facies variations and grain-size trends in the Dahonggou section, consistent with other studies of Cenozoic strata across the basin (e.g., Bally et al., 1986; Hanson et al., 2001; Zhuang et al., 2011). Coarse intervals represent basinward progradation, potentially linked to climate change and/or deformation shifts in source regions, or reduced accommodation in proximal basin settings. Fine intervals may represent changing basin architecture in response to differential subsidence, cutoff of basin outlets and resulting internal drainage, and/or varying sediment discharge in response to changes in climate, source lithologies, and drainage divide position.

The first major shift is represented by the Eocene transition from a low-gradient fluvial system (lower Lulehe Formation) to coarse-grained braided river system (upper Lulehe Formation). The red color, palynologic data (Miao et al., 2008), and weathering indices (chemical index of alteration, plagioclase index of alteration, chemical index of weathering, and modified chemical index of weathering from Song et al., 2013b) are consistent with strong chemical weathering linked to a warm, wet climate. Similar deposits are observed along strike to the northwest at Lake Mahai and Lulehe (Fig. 2; Zhuang et al., 2011), suggesting regional drainage reorganization of a braided fluvial system.

The Xiaganchaigou Formation is characterized by thick multistory channel stacking within a lower-gradient, sandy fluvial system. These deposits can be correlated with coarser meandering fluvial deposits at the Lake Mahai and Lulehe sections (Fig. 2; Yin et al., 2008a, 2008b; Zhuang et al., 2011).

The Shangganchaigou Formation is composed of interbedded paleosol, channel belt, and marginal lacustrine deposits of a low-gradient fluvial system with small lake basins (Fielding, 1986; Miall, 2013). During this phase, the depocenter was situated in the western basin (Ganchaigou and Hongshanhan sections of Zhuang et al., 2011; Fig. 2), where laminated mudstones and evaporites indicate profundal and hypersaline lacustrine environments (Bally et al., 1986). Adjacent to the depocenter, there are marginal lacustrine clastic and calcareous mudstones, including subaerial mud cracks and subaqueous algal stromatolites, ostracods, and pisoliths (Hanson et al., 2001). During endorheic periods, a central evaporative lake was likely quite shallow, resulting in significant shoreline advance and retreat.

Miocene–Pliocene deposits are clast-supported conglomerate and coarse-grained trough cross-bedded sandstones of braided fluvial systems. The interval at Dahonggou exceeds 1000 m, but it thins to ∼600 m at Lake Mahai (Zhuang et al., 2011). Upward coarsening and thickening of this conglomeratic succession are attributed to progradation of a proximal fluvial system, possibly a fluvial megafan (Horton and DeCelles, 2001) derived from a large hinterland catchment area in the southern Qilian Shan fold-thrust belt. This interval is coincident with a regional tectonic reorganization proposed for northeastern Tibet.

SEDIMENT PROVENANCE

A multiproxy provenance approach (paleocurrent, sandstone petrographic, detrital zircon U-Pb geochronologic, and heavy mineral analyses) for the 6200-m-thick Dahonggou section (Figs. 4 and 6) helps to constrain topographic development of several ranges around the basin. To constrain lateral variability, U-Pb results are also presented for the 1400-m-thick Huaitoutala section, ∼100 km to the E (Figs. 3 and 4), where Miocene and younger deposits (magnetostratigraphy from Fang et al., 2007) are exposed in the hanging wall of the Eastern Luliang, Xitie Shan, and Lenghu thrust system of the southern Qilian Shan–Nan Shan.

Paleocurrent Analysis

Paleocurrents from 30 localities help to constrain dispersal patterns and indicate a major switch in paleoflow. At each locality, paleoflow was determined from measurements of at least 10 imbricated clasts or 15 trough limbs. Data (Fig. 9) show a major switch from E-directed flow for the lower Lulehe Formation (four sites, 220–500 m; Fig. 6) to principally SW-directed flow for upper sites (29 sites, 1100–6200 m; Fig. 6). Paleocurrents at the base of the section point to an axial fluvial system flowing from the west, eliminating the northeastern region as a primary source during the early Cenozoic. In contrast, SW-directed flow for upper levels is comparable to the modern configuration, revealing transverse dispersal away from the Qilian Shan–Nan Shan since the Eocene.

Sandstone Petrography

Sixty-seven petrographic thin sections of medium-grained sandstone samples were examined, with 14 selected for point counting using the Gazzi-Dickinson method (Gazzi, 1966; Dickinson, 1970, 1985). These submature, moderately sorted feldspathic litharenites–litharenites (mean composition Q63Ft10Lt27; Folk, 1980; Fig. 10) preserve unstable metamorphic rock fragments (range 3–18%, mean 11%) and heavy minerals including kyanite, sillimanite, garnet, zircon, and chlorite. Polycrystalline quartz with extensive grain-boundary migration is common. Sedimentary rock fragments are also common throughout the section (range 2–14%, mean 7%).

Petrographic data are comparable to sandstones from the same units in the western Qaidam Basin, though our samples contain a higher ratio of metamorphic to volcanic rock fragments (Rieser et al., 2005). Mesozoic sandstones from the northeastern basin are quartz- and feldspar-rich sandstones (Qm62F24Lt14; Fig. 10), with fewer rock fragments (Ritts, 1998). Compositional data indicate recycled orogenic and collisional suture or fold-thrust belt sources (Fig. 10; Dickinson and Suczek, 1979). Sandstone compositions show a minor increase in the abundance of lithic grains between the Shangganchaigou and Xiaganchaigou Formations, followed by a return to slightly higher quartz content in the overlying Shangyoushashan Formation (Fig. 10). Ubiquitous and diverse rock fragments attest to the proximity and/or complex geology of the catchment, but they do not allow for identification of a specific source.

Detrital Zircon Geochronology

Detrital zircon U-Pb geochronologic analyses were employed to constrain Cenozoic sediment provenance for the Qaidam Basin. Shifts in detrital U-Pb ages allow discrimination among potential source areas and may be linked to changes in deformation patterns, shifting upstream drainage configurations, or unroofing of progressively deeper rock units.

Methods

Detrital zircon U-Pb ages were obtained for 19 sandstone samples: 11 from Dahonggou and 8 from Huaitoutala. For each sample, ∼5 kg of medium-grained sandstone were collected and processed according to standard mineral separation techniques. 120 zircon grains were randomly selected from the nonmagnetic heavy mineral fraction of each sample and analyzed by laser ablation–multicollector–inductively coupled plasma–mass spectrometry (LA-MC-ICP-MS) at the Arizona LaserChron Center at the University of Arizona, Tucson, Arizona. Analyses were conducted using a laser spot size of 30 μm, avoiding impurities, and using a Sri Lankan zircon standard for calibration according to Arizona LaserChron Center protocols (Gehrels et al., 2008, 2011). Errors in determining 206Pb/238U, 206Pb/204Pb, and 206Pb/207Pb, result in a measurement error of ∼1–2% (at 2σ level), but this can be substantially larger for grains younger than 1.0 Ga due to low intensity of the 207Pb signal. The 206Pb/238U ages are reported for grains with 206Pb/238U ages younger than 900 Ma, whereas 206Pb/207Pb ages are reported for older grains. Only ages with <20% discordance or <5% reverse discordance were considered robust enough for interpretation. Analytical results (Table DR11) consist of 1999 analyses for 19 samples and are presented as age histograms and kernel density estimates (Fig. 11A) using the DensityPlotter routine outlined by Vermeesch (2012).

Results

U-Pb results help to constrain Cenozoic sediment provenance. Major populations in the Qaidam Basin samples include: (1) Permian–Triassic (299–200 Ma), (2) Late Cambrian to Early Devonian (500–400 Ma), and (3) Paleoproterozoic (2.5–1.6 Ga).

At the base of the section, two lower Lulehe Formation samples (1DH0 and 2DH0) display a bimodal distribution, with peaks at 250 Ma and 436 Ma, and rare Proterozoic grains. Up section, two upper Lulehe Formation samples (2DH224 and 220611-01) yield a unimodal early Paleozoic peak of ca. 440 Ma, a departure from samples below (2DH0) and above (2DH1184). Though rare, there is a significant increase in Proterozoic zircons in upper Lulehe Formation samples (2DH224 and 220611-01). Three Xiaganchaigou Formation samples (2DH1184, 2DH1770, 2DH3040) show diverse zircon age spectra with peaks corresponding to the three aforementioned age populations, including the initial appearance of Late Triassic crystallization ages. Contrasting detrital zircon age spectra are revealed for two Shangganchaigou Formation samples (2DH3774, 2DH4580). Minor Precambrian peaks are present in both samples, though their frequency increases significantly at 2DH4580, continuing to the top of the section. The remaining 10 samples from the Dahonggou and Huaitoutala sections, from the Xiayoushashan, Shangyoushashan, and Shizigou Formations (from Dahonggou: 2DH5484, 2DH6126; from Huaitoutala: 1QUA1, 2QUA144, 2QUA416, 160611-01, 160611-03, 160611-02, 160611-01, 190611-01), contain all three major age populations.

The D value of the Kolmogorov-Smirnov test (K-S test) was used to compare detrital zircon U-Pb age cumulative distribution functions, and to determine the probability that samples are derived from the same parent population. The K-S test D value is calculated as the maximum absolute difference between two samples’ empirical cumulative distribution function and is a measure of the dissimilarity between two samples. Multidimensional scaling was then used to create a spatial visualization or “map” (Fig. 11B) of the dissimilarity measure (Vermeesch, 2013). Three groups of similarly sourced samples were identified using this metric: (1) 1DH0 and 2DH0, both displaying a zircon assemblage consisting of both early Paleozoic and Permian–Triassic ages, (2) 2DH224, 220611-01, and 2DH3774, which are dominated by a single early Paleozoic population, and (3) the remaining samples, which contain a cosmopolitan assemblage of all major detrital zircon populations (Fig. 11B).

Heavy Mineral Analysis

Further discrimination of sources was obtained through heavy mineral (density > 2.89 g/cm3) analyses for sandstones from the Dahonggou section. Heavy minerals generally constitute a minor volume, but they are a useful tracer of source evolution in orogenic belts and stratigraphic shifts in provenance (Mange and Maurer, 1992). Distinct petrologic characteristics of the Qilian Shan–Nan Shan and Kunlun Shan make heavy mineral analysis an important tool for understanding Qaidam Basin provenance.

Methods

Ten heavy mineral fractions were separated from the same Dahonggou samples used for detrital zircon geochronology and analyzed using automated scanning electron microscopy (QEMSCAN) at the Colorado School of Mines, Golden, Colorado. Four energy dispersive X-ray (EDX) spectrometers acquired spectra from each particle with a beam stepping interval of 15 µm, an accelerating voltage of 25 keV, and a beam current of 5 nA; results were analyzed using the control program (iDiscover, FEI). Backscattered electron (BSE) values and EDX spectra were compared with previously determined values and spectra for a representative catalog of diverse mineral phases, facilitating a compositional assignment for each acquisition point and generation of a compositional map for each mount. This approach reflects the instrument design, which is optimized for rapid mineral identification rather than detailed intragranular compositional analysis (for which other techniques are more appropriate).

Results

Results of QEMSCAN heavy mineral analyses revealed three groups with distinct assemblages (Fig. 12). In order to compare heavy mineral assemblages, the index ratios ZTR (zircon-rutile-tourmaline; Hubert, 1962), ATi (apatite-tourmaline index), RZi (rutile-zircon index), and MZi (monazite-zircon index; Morton and Hallsworth, 1994) were used, comparing minerals of similar weathering and diagenetic stability, and hydraulic behavior. In addition, volume percent of each heavy mineral species was calculated for each sample.

Lulehe Formation samples (1DH0, 2DH0, and 2DH224) contain a high volume percent of epidote (average 23.29) with increasing hornblende (average 18.95, increasing from 14.68 to 26.70) and tourmaline (average 0.12, increasing from 0.06 to 0.12). The highest GZi (98.07) and lowest ATi (80.93) are in this group. Whereas alumino-silicate metamorphic accessory minerals andalusite, kyanite, and sillimanite are absent from this sample, these minerals are susceptible to dissolution and/or replacement due to both metamorphic and sedimentary processes. The association of epidote and hornblende, low volume percent of sphene (average 3.29), orthopyroxene (average 0.02), and monazite (average 0.02), and high GZi/ATi ratio (82.5) suggest a source terrane dominated by greenschist-facies metamorphic rocks.

Xiaganchaigou samples (2DH1184, 2DH1770, and 2DH3040) show a significant departure from Lulehe heavy mineral assemblages. There is a switch to low epidote (average 7.78) and hornblende (average 4.93), and decreasing tourmaline (average 0.15, decreasing from 0.20 to 0.12), as well as high volumes of apatite (average 9.37), TiO2 polymorphs (average 30.95), siderite (average 25.67), zircon (average 6.75), monazite (average 0.07), and sphene (average 6.54). This interval has the highest mean ZTR index (37.36), and lowest mean GZi (55.46), RZi (82.39), and MZi (1.20). Zircon, apatite, and sphene are all associated with igneous rocks, and the decrease in epidote and hornblende suggests either dilution or increasing removal of the low-grade metamorphic source. High ZTR indices characteristic of these samples are associated with enhanced recycling (Hubert, 1962).

Upper Xiaganchaigou and lower Shangganchaigou samples (2DH3040 and 2DH3774) have the highest volumes of hornblende (average 25.40), epidote (average 26.38), clinopyroxene (average 0.65), orthopyroxene (average 0.21), sphene (average 6.71), and andalusite/kyanite/sillimanite (average 0.02), yet they have the lowest volumes of apatite (average 0.03), TiO2 polymorphs (average 14.90), and tourmaline (average 0.03). The lowest ZTR index (16.69) and highest ATi, RZi, and MZi are in this zone. The high abundance of metamorphic minerals (hornblende, epidote, pyroxene, andalusite/kyanite/sillimanite), low abundance of typically magmatic minerals (apatite, tourmaline), and high indices correspond to a shift from igneous (2DH1184, 2DH1770, and 2DH3040) to low- to high-grade metamorphic sources (2DH3774, 2DH4580, 2DH5484, and 2DH6126).

Provenance Interpretations

We identified four shifts in sediment provenance based on major changes in paleocurrent orientation, heavy mineral compositions, and detrital zircon age spectra. Consistent with previous research, changes in sandstone framework grain composition at Dahonggou are minor and provide secondary constraints on provenance reconstructions (Rieser et al., 2005). The first shift involves a change from an early Paleozoic and Permian–Triassic zircon U-Pb assemblage to a largely unimodal early Paleozoic assemblage (Figs. 11A and 11B) above 220 m in the Lulehe Formation. The second shift occurs between 770 m and 1100 m in the Lulehe Formation and is characterized by a change in all three major proxies. At 770 m, paleoflow changes from eastward to southwest-directed flow (Fig. 6). At ∼1100 m, there is a shift in U-Pb age spectra from an early Paleozoic assemblage to a cosmopolitan assemblage (Figs. 11A and 11B), and the heavy mineral assemblage shifts toward higher ZTR and lower GZi indices (Fig. 12). This shift is also accompanied by a minor increase in abundance of lithic grains in sandstone (Fig. 10). Between 3040 m and 3774 m, there is a return to a single early Paleozoic–dominated zircon U-Pb age spectrum coupled with a decrease in ZTR and GZi indices in the upper Xiaganchaigou–lower Shangganchaigou Formations. There is also a return to more compositionally mature sandstones. Finally, the cosmopolitan zircon U-Pb age assemblage reemerges between 4580 m and 5484 m (uppermost Shangganchaigou Formation–lower Xiayoushashan Formation) in the Dahonggou section and synchronously dominates all samples in the Huaitoutala section (Fig. 11A).

E-directed paleocurrent orientations in the lower Lulehe Formation indicate sediment sources to the south or west of the Dahonggou location. Potential sources include (1) the Hoh Xil Basin, (2) Altyn Shan, (3) Western Kunlun Shan, and (4) Eastern Kunlun Shan.

(1) Detrital zircon age spectra (1DH0 and 2DH0) lack the large Neoproterozoic, Paleoproterozoic, and Cretaceous populations of contemporaneous Lhasa-Qiangtang—derived deposits in the Hoh Xil Basin to the south (Fig. 5; Dai et al., 2012; Staisch et al., 2014), suggesting the Paleogene appearance of a Kunlun Shan topographic barrier between the Qaidam and Hoh Xil Basins.

(2) Initial displacement on the Altyn Tagh fault and uplift of the Altyn Shan are usually ascribed to Oligocene–early Miocene time, with possible Mesozoic–early Cenozoic precursor faulting (Jolivet et al., 2001; Sobel et al., 2001; Yin et al., 2002, 2007a; Cowgill et al., 2003; Darby et al., 2005; Wang et al., 2005; Ritts et al., 2008). Because the Altyn Shan contains a major 1800–2000 Ma source (Yue et al., 2005; Gehrels et al., 2003a, 2003b) that is not observed in lower Lulehe results, it is eliminated as a possible source.

(3) On the basis of basin architecture and apatite fission-track data, Wang et al. (2006) interpreted an ∼2000 km “paleo–Kunlun River” sourced in the Western Kunlun Shan to northeastern Pamir, following the incipient Altyn Tagh and Karakash faults. The expected provenance signature, however, would include early Paleozoic and Cretaceous zircons but would be dominated by Permian–Late Triassic ages (Fig. 5E; Cowgill et al., 2003; Zhang et al., 2007; Robinson et al., 2012; Carrapa et al., 2014). This association of ages is inconsistent with basal Lulehe results, which are dominated by early Paleozoic–aged zircons and include no Cretaceous ages (Fig. 11). However, exclusion of the Cretaceous-aged population would make the Western Kunlun Shan a viable secondary sediment source.

(4) Permian–Triassic and early Paleozoic age signatures from the lower Lulehe Formation (1DH0 and 2DH0) are compatible with the Kunlun arc terrane of the Eastern Kunlun Shan (west of Golmud) and southern Altyn Shan (Cowgill et al., 2003) and Triassic sedimentary units of the Kunlun Shan (Ding et al., 2013), and they are incompatible with Late Triassic plutons of the easternmost Eastern Kunlun Shan (Mo et al., 2007). The high abundance of epidote, hornblende, tourmaline, and garnet is consistent with a metamorphic basement source. Whereas low-temperature thermochronologic results from the easternmost Eastern Kunlun Shan support Oligocene–Miocene cooling (Mock et al., 1999; Jolivet et al., 2003; Yuan et al., 2006), this exhumation may be related to S-directed backthrusting following earlier Paleocene–Eocene N-directed thrusting in the Eastern Kunlun Shan (e.g., Yin et al., 2007a). Therefore, the westernmost Eastern Kunlun Shan is favored as the principal sediment source.

The switch to an early Paleozoic detrital zircon signature above 220 m occurs in the absence of a change in paleocurrent orientation, suggesting consistent derivation from the southern Qaidam Basin margin. The unimodal early Paleozoic zircon U-Pb age spectrum is consistent with a westward shift in active exhumation to the northwestern Qimen Tagh region of the Eastern Kunlun Shan, adjacent to the Altyn Shan (W. Li et al., 2013) feeding the E-directed axial fluvial system.

There is a clear switch from variable axial paleocurrents of the lower–middle Lulehe Formation to SW-directed, transverse paleocurrents at 770 m in the uppermost Lulehe and Xiaganchaigou Formations (Fig. 6). At the Dahonggou section, this interval is characterized by a cosmopolitan detrital zircon signature with Permian–Triassic, early Paleozoic, and Proterozoic age peaks. Samples from the Xiaganchaigou Formation (2DH1184, 2DH1770, 2DH3040) include Paleozoic and Paleoproterozoic grains matching the Qilian Shan–Nan Shan basement (Fig. 5; Cowgill et al., 2003; Gehrels et al., 2003a, 2003b; Menold, 2006; Yan et al., 2010). The origin of Permian–Triassic grains is more complex, because such plutons are less common in the Qilian Shan–Nan Shan than the Eastern Kunlun Shan. Given SW-directed paleocurrents and high ZTR heavy mineral index, the Permian–Triassic peak is attributed to grains recycled from Jurassic deposits situated to the north in the Qilian Shan–Nan Shan, but originally shed from the Permian–Triassic Kunlun arc (Fig. 5). The heavy mineral assemblage represents dominantly igneous and recycled sedimentary sources, consistent with populations in the Qilian Shan–Nan Shan, where a well-established catchment eroded diverse source rocks. The absence of a significant metamorphic component, common in the North Qaidam terrane, suggests that high topography was confined to the northern Qilian Shan during Xiaganchaigou deposition.

A pronounced 400–500 Ma detrital zircon age peak and corresponding reduction or disappearance of other age peaks in the two Shangganchaigou Formation samples (2DH3774, 2DH4580) could be explained by several mechanisms. First, units containing recycled Permian–Triassic and Paleoproterozoic grains may have been preferentially eroded away, leaving only lower Paleozoic rocks. Second, new drainage networks may have contained local Paleozoic sources, before drainage expansion and elaboration provided more diverse populations. Neither option, however, is likely to have completely removed two source populations while maintaining early Paleozoic sources. A third explanation involves in-sequence deformation and growth of a topographic barrier along the Nan Shan–North Qaidam thrust belt (NQF in Fig. 1). At the same level, the increase in epidote and hornblende abundance and switch to high GZi and RZi heavy mineral indices require deeper unroofing and/or exposures of metamorphic rocks. Paleoflow was to the southwest, consistent with a north/northeast source in the Qilian Shan. Possible sources include ophiolitic and ultrahigh-pressure metamorphic assemblages in the Nan Shan–North Qaidam thrust belt along the Luliang and Xitieshan thrusts (Yin et al., 2008a), or possibly farther north in the Qilian Shan–Nan Shan.

Detrital zircon ages for the uppermost Xiayoushashan–Shizigou Formations at both sections contain three major populations, including the early Paleozoic and reintroduction of Permian–Triassic and Proterozoic populations (Fig. 11). Metamorphic minerals continue to dominate heavy mineral assemblages. Lithofacies, sedimentation rates, and provenance shifts indicate continued in-sequence advance of the Luliang-Xitieshan thrust system along the southern front of the Qilian Shan–Nan Shan. Age populations can be traced to the southern Qilian Shan–Nan Shan by erosional unroofing of Paleoproterozoic basement, lower Paleozoic magmatic rocks, isolated Permian–Triassic plutons related to the Kunlun arc, as well as Mesozoic sedimentary units, collectively resulting in well-mixed detrital zircon age spectra.

DISCUSSION

Basin Reconstruction

Sedimentologic, provenance, and geochronologic data support a reconstruction of Qaidam Basin evolution (Fig. 13) involving multiple source regions and variations in drainage patterns dictated by shifts in the tectonic development of flanking mountain ranges. Results presented here document the Cenozoic depositional record in the northeastern Qaidam Basin, and its relation to the broader northern Tibetan Plateau. For this reconstruction, the chronology of the Dahonggou section is roughly correlated to nearby magnetostratigraphic sections (e.g., Lu and Xiong, 2009; Ke et al., 2013).

Lulehe Formation (Paleocene)

The composition and sediment dispersal pathways of Paleocene Lulehe Formation sediments are consistent with derivation from the westernmost Eastern Kunlun Shan (Fig. 2). This observation suggests that the Eastern Kunlun Shan had developed into a topographic edifice, largely or totally blocking drainage from the south (Fig. 13A). The drainage divide would have shed Permian–Triassic and early Paleozoic detritus to both the Qaidam and Hoh Xil Basins, and blocked the majority of Neoproterozoic detritus from the Qiangtang and other southern terranes (Fig. 5). Detrital zircon U-Pb data and carbonate rock fragments in the Lulehe Formation in western Qaidam confirm early exhumation of the Eastern Kunlun Shan (Cheng et al., 2015). Finally, E- and NE-directed paleocurrent orientations at the Dahonggou location on the northern margin of the Qaidam Basin also strongly suggest that the Qilian Shan was not a significant topographic barrier in the Paleocene.

Lulehe Formation (Paleocene–Eocene Transition)

The increased fluvial competence indicated by the upward-coarsening trend may be attributed to increased precipitation and runoff or to enlargement of the catchment (Fig. 13B). However, the shift to a unimodal early Paleozoic detrital zircon peak in samples 2DH224 and 220611-01 (Fig. 11) is inconsistent with enlargement of the fluvial catchment. Oxygen isotopic values from sedimentary cements and paleosol carbonates from around the Qaidam Basin decrease from the lower to upper Lulehe Formation (Kent-Corson et al., 2009), consistent with enhanced precipitation in the late Paleocene–early Eocene. Low ZTR indices through this interval (Fig. 12), despite enhanced weathering in the warm Paleocene–Eocene climate, suggest rapid denudation and sediment transport to the Dahonggou site consistent with the increased grain size and bed-load proportions. Continued E-oriented paleoflow during this interval suggests quiescence in the Qilian Shan, continued sourcing of sediment from the SW-W basin margin, and a Qaidam Basin that remained open during this time interval.

Uppermost Lulehe–Xiaganchaigou Formations (Early–Late Eocene)

Evidence for Eocene shortening in the Qilian Shan–Nan Shan (Figs. 13C and 13D) is revealed in the Dahonggou section by shifts in lithofacies and sediment dispersal. Results indicate Eocene establishment of the northern basin margin along the southern Qilian Shan–Nan Shan north of Daqaidam (Fig. 1), with the transition to the cosmopolitan zircon assemblage at sample DH1184. Provenance data presented here are complementary to the stratigraphic synthesis of Zhuang et al. (2011), in which they recognized Eocene accumulation of coarse sediment along the southern front of the Qilian Shan and Nan Shan–North Qaidam fold-thrust belt.

Shangganchaigou Formation (Oligocene)

The transition to a dominant early Paleozoic peak during the Oligocene–early Miocene with continued SW-directed paleoflow of a meandering fluvial system requires the southward advance of the Nan Shan–North Qaidam thrust belt within the Qilian Shan (Fig. 13E). This predicts isolation from earlier sources and a huge early Paleozoic signal from plutons along the Qilian Shan–Nan Shan thrust front. A feature consistent with this interpretation is the drastic shift in heavy minerals, from a recycled sedimentary to metamorphic source. The Nan Shan–North Qaidam thrust belt contains an early Paleozoic suture zone (Mattinson et al., 2007; Xiao et al., 2009) that underwent high- and ultrahigh-pressure metamorphism, and that is exposed along the Luliang and Xitie Shan faults (Fig. 4; Yang et al., 2001; Yin et al., 2008a; Menold et al., 2009).

The Qaidam Basin was fully closed at this stage, with organic-rich lacustrine deposits accumulating in the basin center (Hanson et al., 2001). At the Dahonggou site, marginal lacustrine deposits are interbedded with paleosols and meandering fluvial deposits throughout this interval, as the depocenter shifted to the southeast due to flexural loading by the Nan Shan–North Qaidam thrust belt (Yin et al., 2008a, 2008b).

Xiayoushashan, Shangyoushashan, and Shizigou Formations (Miocene–Pliocene)

During Miocene deposition of the Xiayoushashan and Shangyoushashan Formations, the Qaidam depocenter shifted eastward (Wang et al., 2006); the deposits begin coarsening upward at the 4800 m level in the Dahonggou section, and sedimentation initiated at the Huaitoutala section (Fang et al., 2007) in a piggyback hanging-wall setting on the S-directed Luliang and Xitieshan thrust system (Fig. 13F).

Growth strata, regional unconformities, and shifts in depositional systems and accumulation rates support a deformation phase commencing at 12–8 Ma and accelerating after ca. 3.6 Ma (Clark et al., 2010; Zheng et al., 2010; Lease et al., 2011; Yuan et al., 2013). This reorganization includes diminished slip along the Altyn Tagh fault and deformation migration in the southern Qilian Shan–Nan Shan to the south and east with the development of the Eastern Luliang, Xitie Shan, and Lenghu thrusts (Yue et al., 2005; Yin et al., 2008a). In the northwestern Qaidam Basin, a shift from magmatic to metamorphic sources (Jian et al., 2013) is attributed to unroofing of the Nan Shan–North Qaidam and South Qilian Shan belts. At the Dahonggou section, diverse detrital zircon and heavy mineral assemblages imply drainage integration and source unroofing, producing diverse grains from Paleoproterozoic basement, lower Paleozoic magmatic rocks, and Permian–Triassic grains from recycled sedimentary, igneous, and metamorphic sources.

Miocene and younger deposits preserve >2000 m of coarse-grained fill, with an upward-coarsening and upward-thickening pattern reflecting progradation of large distributive fluvial lobes from range-front outlets. This succession is suggestive of the large area (>1000 km2), low slope, and waterlain facies of fluvial megafans (Gupta, 1997; Horton and DeCelles, 2001; Hartley et al., 2010). During middle Miocene reorganization of deformation, a basinwide grain-size increase and initiation of deposition in the Delingha depression (subbasin containing the Huaitoutala section) represent significant topographic growth in surrounding ranges.

Modern drainage configuration and source areas were established during late Miocene–Pliocene deposition of the Shangyoushashan and Shizigou Formations (Pullen et al., 2011). Coarse-grained sediment accumulation and growth strata formation continue to characterize the modern margins and interior of the Qaidam Basin (Wang and Coward, 1990).

Regional Implications

Basin Isolation

The pattern of basin isolation for the Qaidam region (Fig. 13) is consistent with models for plateau expansion involving structural incorporation and partitioning of once-giant peripheral basins (Métivier et al., 1998; Meyer et al., 1998; Tapponnier et al., 2001; Horton et al., 2004; Carroll et al., 2010; Craddock et al., 2011b; Horton, 2012). In this model, growth of new basin-bounding ranges forms critical topographic barriers to surface flow and atmospheric circulation, promoting sediment ponding and arid rain-shadow conditions within the basin (Sobel et al., 2003; Strecker et al., 2007). Following surface uplift, basin-margin topographic barriers and associated arid conditions inhibit erosional evacuation. Progressive partitioning of large intraplate basins is supported by studies suggesting a shared Mesozoic history for the northwestern Qaidam Basin and the southeastern Tarim Basin (Ritts and Biffi, 2000, 2001; Meng and Fang, 2008; Craddock et al., 2011a; Johnson and Ritts, 2012). Tectonic subsidence waned in the Early Cretaceous and was followed by regional burial (Jolivet et al., 2001) and potential formation of a giant peripheral basin on the northern margin of the Tibetan Plateau. The sedimentary record presented here indicates that the early Cenozoic Kunlun Shan largely isolated the Hoh Xil and Qaidam Basins by the Paleocene and would have constituted the southern margin of a formerly unified Qaidam-Tarim Basin (Ritts and Biffi, 2000; Figs. 13A and 13B).

Partitioning of this giant basin and progressive isolation of the Qaidam Basin from the Tarim Basin (to the WNW) and the Hexi Corridor (to the NE) began in the Paleogene. Paleocene to early Eocene isolation from the Hexi Corridor is documented by the provenance shift and SW-directed paleoflow in the uppermost Lulehe Formation, and displacement along thrusts bounding the NE margin of the Qaidam Basin (Fig. 13C; Yin et al., 2008a). Late Eocene–early Oligocene isolation from the Tarim Basin via movement on the Altyn Tagh fault is documented by thermochronology and basin paleogeography for the NW Qaidam region (e.g., Yue et al., 2004, 2005; Wang et al., 2005; Meng and Fang, 2008; Ritts et al., 2008). By the middle Oligocene, the Qaidam Basin was tectonically isolated from basins to the north and west. SW-directed paleoflow at Dahonggou and endorheic lacustrine deposition in western Qaidam Basin (Hanson et al., 2001) indicate that by the Oligocene, the Qaidam Basin was also hydrologically closed. Southward encroachment of Qilian Shan fold-thrust deformation into the Qaidam Basin focused sedimentation in transverse fluvial systems along basin margins (Figs. 13C–13E). Connections to basins of the northeastern Tibetan Plateau are less clear. Increased exhumation is not documented in the Dulan-Chaka highland E of Qaidam Basin until the early Miocene (Duvall et al., 2013), raising the possibility that the Qaidam Basin remained tectonically integrated with regions to the east until at least the early Miocene.

The early Eocene onset of basin isolation is consistent with previous research suggesting that strain from India-Asia collision was transmitted rapidly to the northern edge of the modern Tibetan Plateau (e.g., Horton et al., 2004; Dupont-Nivet et al., 2004; Yin et al., 2008a; Clark et al., 2010). This may have been facilitated by high topography in the southern Tibetan Plateau (Ding et al., 2014), strong lithospheric blocks underlying the Qaidam-Tarim Basin (Braitenberg et al., 2003), and preexisting zones of weakened lithosphere in basin-bounding mountain ranges due to pre-Cenozoic deformation (Ritts and Biffi, 2001; Wang et al., 2005; Johnson and Ritts, 2012). Oxygen isotopic data from basins in the vicinity of Qaidam Basin indicate topographic growth from Eocene to Oligocene, coincident with the timing of basin isolation (Graham et al., 2005; Kent-Corson et al., 2009). Surface uplift of Qaidam Basin and its bounding ranges may have stopped or slowed due to a potential early–middle Miocene shift from largely N-S shortening to a deformation regime dominated by eastward translation of crustal blocks along continental-scale strike-slip faults (Zhang et al., 2004; Gan et al., 2007; Duvall and Clark, 2010; Zheng et al., 2010; Lease et al., 2011; Yuan et al., 2013).

Factors Favoring Development of Giant Basins

The Qaidam Basin contrasts with basins in the southern Tibetan Plateau, as well as the mosaic of basins to its east in the northeastern Tibetan Plateau. Middle Miocene–Quaternary basins in the southern and central Tibetan Plateau are commonly related to extension or strike-slip faulting (e.g., Garzione et al., 2003; Spurlin et al., 2005; Kapp et al., 2008; Saylor et al., 2010; Sanchez et al., 2010; Horton, 2012; Taylor et al., 2012; Woodruff et al., 2013). Older basins in central Tibet such as Lunpola, Nima, Linzhou, and Hoh Xil (Fig. 1) formed in response to, or coeval with, shortening, and often contain only an early Cenozoic record of sedimentation, after which they were exhumed during ongoing deformation (Xu, 1984; Leeder et al., 1988; Harrison et al., 1992; Liu and Wang, 2001; Liu et al., 2001; Horton et al., 2002; DeCelles et al., 2007; He et al., 2007; Kapp et al., 2007). However, the southern and central Tibetan basins are considerably smaller, thinner, and have shorter life spans than the Qaidam Basin. Similarly, a collection of subbasins in the northeastern Tibetan Plateau have thinner stratigraphic sequences and discontinuous sedimentation, reflecting late Eocene to late Miocene segregation events (Horton et al., 2004; Dupont-Nivet et al., 2008; Lease et al., 2011), and ultimately drainage integration and evacuation in the Pleistocene (Li et al., 1997; Fang et al., 2005; Craddock et al., 2010; Nie et al., 2015).

The Qaidam region shares elements with other large basins on the periphery of the Tibetan Plateau, including the Tarim, Sichuan, and Hexi Corridor regions (Fig. 1), but it also differs from them in several fundamental ways. The Qaidam and Tarim Basins preserve the thickest and most complete Cenozoic sedimentary sequence of any of these basins, with up to 15 km of Cenozoic sediment accumulation. In the Hexi Corridor, 2.5–5 km of sediment accumulated in a restricted broken foreland basin above the basal Cretaceous–Oligocene unconformity (Bovet et al., 2009, and references therein). The Sichuan Basin contains thick sequences of Mesozoic strata, but only a thin (<1 km) veneer of Cenozoic deposits (e.g., Guo et al., 1996; Meng et al., 2005).

High-elevation basins of the southern and central Tibetan Plateau and those on the periphery are also distinguished by the strength of their crustal basement relative to that of the surrounding mountain ranges, and development of syndepositional structural dams suggesting that these factors may control their evolution and accumulation histories. Low effective elastic thicknesses (typically <20 km) of basement underlying basins in the southern and central Tibetan Plateau (Masek et al., 1994; Braitenberg et al., 2003; Jordan and Watts, 2005), and a weak layer in the middle crust–upper mantle (Jin et al., 1994; Owens and Zandt, 1997; Klemperer, 2006) allow limited evolution of overlying basins before basins are disrupted and recycled (e.g., Kapp et al., 2008). The northeastern Tibetan region is characterized by lithosphere with moderately higher but laterally heterogeneous effective elastic thickness (Braitenberg et al., 2003; Jordan and Watts, 2005; L. Li et al., 2013), resulting in strain focusing in zones of weakness, disruption of once-contiguous basins during ongoing deformation, and development of the mosaic of basins. The Hexi Corridor is characterized by ∼20–30 km effective elastic thicknesses, but it is separated from the Qaidam Basin to the south by a low (<10 km) effective elastic thickness zone in the Qilian Shan–Nan Shan.

The Qaidam, Tarim, and Sichuan Basins are characterized by large regions with high effective elastic thicknesses, and high basement seismic velocities extending through the lithosphere (Braitenberg et al., 2003; Jordan and Watts, 2005; Liang and Song, 2006; Li et al., 2008a, 2008b; Li et al., 2012; Yang et al., 2012). Whereas the Tarim and Qaidam Basins developed thick Cenozoic sedimentary sequences by ponding due to Cenozoic deformation focused on the basin margins, no similar syndepositional structural dam developed in the Sichuan Basin, resulting in widespread erosion of basin strata (Richardson et al., 2008). The onset of isolation and uplift of the Qaidam Basin relative to the Tarim Basin may have been due either to injection of Tibetan lower crust beneath the Qaidam Basin (Yin et al., 2007a; Karplus et al., 2011) or structural transformation of the Qaidam Basin into a crustal-scale piggyback basin (Yin et al., 2007a; Horton, 2012).

Effect of Climate on Basin Filling and Evacuation

Finally, late Eocene–early Oligocene isolation of the Qaidam Basin predates a regional increase in aridity starting in the early–middle Miocene (Kent-Corson et al., 2009; Miao et al., 2011). Endorheic conditions in Qaidam Basin also predate the increase in aridity and persist through the modern era. This sequence of events points to early defeat of range-traversing rivers by tectonic uplift and development of rain-shadow–related aridity within an isolated Qaidam Basin prior to Miocene regional aridification. On the other hand, Cenozoic strata in the multiple subbasins of northeastern Tibet point to periods of both sediment accumulation and large-scale evacuation of basin fill following regional aridification (e.g., Fang et al., 2003; Craddock et al., 2010, 2011b; Hough et al., 2011, 2014). As in the Qaidam Basin, local rain-shadow aridity is superimposed on the regional signal of aridification (Hough et al., 2014). However, excavation of the subbasins following regional aridification indicates that despite increased regional aridity in the late Cenozoic, rivers maintained sufficient stream power to overcome structurally produced topographic barriers and effectively evacuate sediments that ponded upstream of those topographic barriers. The decoupling of regional climate changes and sediment accumulation or evacuation implies that regional climate change was not the primary factor in either the development or destruction of giant plateau-margin basins.

CONCLUSIONS

Constraints on the evolution of the Qaidam Basin provide insights into the Cenozoic deformation history of the northern Tibetan Plateau, which broadly include Paleocene–Eocene exhumation of the Kunlun Shan followed by activation and then outward growth of the Qilian Shan–Nan Shan. Supporting evidence consists of distinctive patterns and temporal shifts in depositional systems, detrital provenance, and sediment dispersal in the northeastern Qaidam Basin.

(1) Initial Cenozoic sedimentation in the northeastern Qaidam Basin involved eastward transport of sediment likely derived from the westernmost Eastern Kunlun Shan. Paleocene deposits (lower Lulehe Formation) contain major Permian–Triassic and early Paleozoic detrital zircon populations, with a heavy mineral assemblage of principally metamorphic origin, consistent with the Paleozoic basement, Permian–Triassic magmatic arc, and older metamorphic signatures along the Kunlun suture. These U-Pb detrital zircon ages support the existence of a drainage divide between the Qaidam Basin and the Hoh Xil Basin, preventing large volumes of Neoproterozoic grains of the Qiangtang terrane from reaching the Qaidam Basin. In the upper Lulehe Formation, a shift from sandstone- to conglomerate-dominated facies and a transition to early Paleozoic U-Pb detrital zircon ages indicate continued exhumation of the Eastern Kunlun Shan.

(2) The late Paleocene–Eocene shift is followed by an early Eocene switch to SW-directed paleoflow, igneous and recycled sedimentary heavy mineral assemblage, and diverse Permian–Triassic, early Paleozoic, and Proterozoic detrital zircon populations. This provenance analysis is consistent with the onset of extensive exhumation in the Qilian Shan–Nan Shan.

(3) In-sequence advance of the Qilian Shan–Nan Shan fold-thrust belt occurred during the Oligocene–early Miocene, with exhumation of metamorphic and ophiolitic complexes along the frontal Eastern Luliang, Xitie Shan, and Lenghu thrust faults and initial sedimentation at the Huaitoutala section located in the hanging wall of this frontal thrust belt. The Oligocene Shangganchaigou Formation preserves a significant provenance shift at the Dahonggou section, to a metamorphic heavy mineral assemblage and a unimodal early Paleozoic detrital zircon U-Pb peak.

(4) Since the middle Miocene, exhumation of the Qilian Shan–Nan Shan is recognizable in the detrital record preserved within the Qaidam Basin. Provenance data from both the Dahonggou and Huaitoutala sections show unroofing of Permian–Triassic, lower Paleozoic, and Paleoproterozoic sources, with metamorphic heavy mineral assemblages, in the Qilian Shan–Nan Shan through the Xiayoushashan and Shangyoushashan Formations, continuing through the Pliocene Shizigou Formation.

(5) During the Paleogene, the Qaidam Basin was sequentially isolated from the south, the west, and finally the north, becoming hydrologically closed by the Oligocene. The basin was eventually uplifted as a part of the Tibetan Plateau and further isolated from the basins of the northeastern Tibetan Plateau. This pattern of basin isolation and subsequent uplift demonstrates the critical role of basin formation and filling as a component in the growth of the Tibetan Plateau.

(6) This research joins a growing number of studies demonstrating the utility of an integrated approach to sediment provenance and quantitative comparison of detrital geochronology data. This is particularly true in regions where multiple sources have similar lithologies, compositions, or thermal histories or where recycling extensively masks provenance signals.

This research was supported by the University of Texas at Austin Jackson School of the Geosciences Ronald K. DeFord Field Scholarship, the American Association of Petroleum Geologists L. Austin Weeks Memorial Grant, a Geological Society of America graduate student grant, and an International Association of Sedimentologists Post-Graduate Grant Scheme awarded to Bush, and a subaward to Saylor and Horton from the National Science Foundation (NSF) grant EAR-0958704, awarded to Xiaoming Wang. Field work was facilitated by awards to Nie from the National Basic Research Program of China (grant 2013CB956400) and National Natural Science Foundation of China (grant 41422204) and logistical assistance from Yunfa Miao, Kuang He, Wenbin Peng, and Shunchuan Ji of Lanzhou University. We thank the staff at the Arizona LaserChron Center (NSF grant EAR-1032156) and Katharina Pfaff at the Colorado School of Mines QEMSCAN facility for laboratory assistance. The manuscript was improved by careful reviews by William Craddock, an anonymous reviewer, and editor Arlo Weil.

1GSA Data Repository Item 2015362, Table DR1. Sample locations, U-Pb isotopic results with 1σ uncertainties (rejected data in red italicized font), heavy mineral abundance data, sandstone petrographic point count data and paleocurrent measurement; Figure DR1. Concordia diagrams for all samples used in this study plotted using the Isoplot routine of Ludwig (2003), is available at www.geosociety.org/pubs/ft2015.htm, or on request from editing@geosociety.org.