Provenance analyses
Paleoweathering and sedimentary recycling
The geochemical features of clastic rocks are strongly influenced by the existence and degree of chemical weathering and sedimentary recycling (Nesbitt and Young 1982; Krzeszowska 2019; Bokanda et al. 2021; Omietimi et al. 2022), thus a variety of different weathering conditions can be used to evaluate chemical weathering history of the source area (Nesbitt and Young 1982; McLennan et al. 1993; Fedo et al. 1995; Lewin et al. 2018). CIA values can be calculated by the formula [Al2O3/(Al2O3 + CaO* + Na2O + K2O)] * 100 (molar proportion; Nesbitt and Young 1982), and the CaO* contents were calculated based on the method described by McLennan et al. (1993). Generally, intense weathering in the source area may result in the increase of CIA values (80–100) in sediments, whereas weak weathering may cause the sediments to have relatively low CIA values (50–70) (Yan et al. 2010; Fig. 6).
The CIA values for all collected mudstones (46.0-68.9, avg. 59.1) are lower than those of the PAAS (70.36), which indicate weak chemical weathering in the source areas (Table 1, Fig. 6). The A-CN-K ternary plot is widely applied to evaluate the degree of weathering in the source areas as well (molar proportion; Fedo et al. 1995). On the A–CN–K plot, all collected mudstone samples are plotted near to the plagioclase–K-feldspar join line and clustered between granodiorite and granite average compositions (Fig. 6). In addition, the linear weathering trend of these mudstones on the A–CN–K plot reflects that the source areas were stable (Fig. 6). These results indicate that the source areas of the collected mudstones were affected by a weak chemical weathering.
The index of compositional variability (ICV= (TFe2O3 + K2O + Na2O + CaO + TiO2) / Al2O3) is widely used to evaluate sediment recycling and maturity (Cox et al. 1995; Awasthi 2017; Sahariah and Bhattacharyya 2019; Li et al. 2022). The high ICV value (> 1) indicate the first-cycle products in the tectonically active area, while the low ICV value (< 0.84) suggest intense weathering and multiple sedimentary cycles (Van de Kamp and Leake 1985; Cox et al. 1995; Chen et al. 2014). The ICV values of the Liwaxia to Naijiahe formations are between 0.85 and 13.85 (Table 1), with an average of 3.60, which is significantly greater than that of the PAAS (0.80). These values indicate that the Liwaxia to Naijiahe formations are first-cycle and compositionally immature. The Th/Sc and Zr/Sc ratios are widely used to deduce the sorting degree, compositional maturity, and heavy mineral accumulation of clastic sediments (McLennan et al. 1993; Dypvik and Harris 2001; ArmstrongAltrin et al. 2012; Wang et al. 2018). The Th/Sc and Zr/Sc ratios of all collected samples are 0.45–1.26 and 5.18–13.80, respectively. These data combined with Th/Sc–Zr/Sc diagram indicate that the LPS mudstones were little affected by sedimentary recycling (Fig. 7).
Type of source rocks
The detrital components and elemental composition of clastic rocks are controlled by its provenance (Verma and Armstrong-Altrin 2013; Xie et al. 2018; Li et al. 2022). Some trace elements (e.g., La, Zr, Th, Sc, Hf, Co) are non-migrating and provide a reliable indication of provenance (Li et al. 2022). In the provenance discrimination diagrams of these trace elements (Zr/Sc versus Th/Sc, Hf versus La/Th, and La/Sc versus Co/Th), most samples fall into the field of felsic rocks (Fig. 7). The mixed felsic/basic source was suggested in the Hf versus La/Th diagram (Fig. 7). However, the low MgO and TiO2 contents indicate the depletion of Mg and Ti, and less mafic minerals, further excluded an abundant basic provenance (McCann 1991; Tao et al. 2016; Zhang et al. 2020; Li et al. 2022). The Al2O3/TiO2 ratio is sensitive to the change of parent rocks: high ratios (19–28) indicate a felsic parent rock and low ratios (< 14) indicate mafic parent rocks (Table 1, Girty et al. 1996). As illustrated on Al2O3 versus TiO2 diagram, the Al2O3/TiO2 ratios (22.84–33.44) of all samples show a felsic rocks feature (Fig. 7). The TiO2/Zr ratios increase obviously from felsic to mafic source, thus the TiO2 versus Zr diagram can distinguish mafic, intermediate, and felsic rocks (Hayashi et al. 1997; Armstrong-Altrin et al. 2015a, b; Moradi et al. 2016; Wang et al. 2017c;Wang et al. 2018). This diagram suggests that the collected samples originated from felsic and intermediate rocks (Fig. 7). An acidic rocks source was also suggested by TiO2 versus Ni diagram (Fig. 7). In addition, the REE’s features can be used to infer the source of fine-grained sedimentary rocks: felsic source commonly displays higher LREE/HREE ratios and negative Eu anomalies, whereas mafic source demonstrates low LREE/HREE ratios and no pronounced Eu anomalies (Taylor and McLennan 1985; Roddaz et al. 2006; Kasanzu et al 2008). All samples show the characteristics of LREE enrichment and HREE depletion with obvious negative Eu anomalies, further denoting a felsic provenance (Fig. 5). Abundant intermediate-acidic rocks exposed in eastern Qilian and western Qinling areas (southwestern to the study area), which has been confirmed to be the provenance of the Sanqiao-Heshangpu formations of the Liupanshan Basin (Zhao et al. 2020). In addition, the detrital zircon ages distributions of Liwaxia and Madongshan formations are consistent with the those of the eastern Qilian and western Qinling orogens (Ning 2017). Hence, the provenance of the Liwaxia-Naijiahe Formations is generally inherited from the Sanqiao-Heshangpu period.
Based on the analysis above, it can be concluded that the provenance of the sediments of the Liwaxia-Naijiahe formations were mainly felsic acidic rocks of the eastern Qilian and western Qinling areas.
Paleoenvironment conditions
Paleoclimate condition
The Sr/Cu ratio is an important indicator to reveal the paleoclimate conditions: low Sr/Cu ratios (< 5.0) suggest a humid climate condition, and high Sr/Cu ratios (> 5.0) indicate an arid climate condition(Wang et al. 2017;Yu et al. 2021; Lerman 1989; Meng et al. 2012; Cao et al. 2015). The Sr/Cu ratios of the Liwaxia to Naijiahe formations are mostly > 5.0 (1.91–47.36, avg., 12.44; 7.62–33.84, avg. 23.58; and 1.71–284.52, avg. 76.02 respectively) (Fig. 8), indicating a semiarid-arid climate condition. The Sr/Cu values of some samples (e.g., NX-20-41 NX-20-43) of the Naijiahe formation are extremely high, with extremely high Sr and Ba contents, which probably indicate extreme arid and evaporative environments (e.g., Dai et al. 2021). The Sr/Ba ratio is not only used for reconstruction of the paleosalinity, but paleoclimate conditions as well (Wang et al. 2018). Generally, high Sr/Ba ratios (> 1.0) indicate an arid climate conditions and low Sr/Ba ratios (< 0.50) reflect a humid climate condition (Meng et al. 2012; Fu et al. 2016; Wang et al. 2018). The high Sr/Ba ratio of the Liwaxia to Naijiahe formations (0.12–3.15, avg. 0.74; 0.80–3.22, avg. 2.16, and 0.31–15.13, avg. 4.71) suggest a general arid climate conditions as well. The Rb/Sr ratio also an important indicator to reveal the paleoclimate conditions due to different geochemical properties (Chen et al. 2022). The low Rb/Sr ratios suggest an arid climate condition, and high Rb/Sr ratios indicate a humid climate condition (Zheng et al. 2015; Ma et al. 2019; Chen et al. 2022). The low Rb/Sr ratio of the Liwaxia to Naijiahe formations (0.07–1.75, avg. 1.06; 0.10–0.45, avg. 0.22, and 0.01–1.32, avg. 0.22) also suggest an arid climate condition (Fig. 8). The carbon isotope values of the Pseudofrenelopsis leaf collected from the Naijiahe Formation also indicate an arid or semi-arid climate during the late Albian (Du et al. 2018). In addition, the discovery of Pseudofrenelopsis and Caddisfly (Du et al. 2014; He et al. 2014), the dominant Classopollis (Li and Du 2006; Zhang et al. 2012), high pCO2 estimates (Du et al. 2014), occurrence of gypsum layers and evidence of the carbon and oxygen isotopes (Li et al. 2013; Du et al. 2014) also support an arid climate condition during the late Early Cretaceous.
Paleoredox condition
Trace element ratios such as V/Cr, V/Ni, and V/Ni + V are considered to be credible redox indicators (Tonger et al. 2004; Liu et al. 2007; Li et al. 2020). Scheffler (2006) demonstrated that V/Cr ratio < 2 indicate an oxic condition, 2-4.25 reflect a dysoxic condition, and > 4.25 reveal an anoxic environment. V/Ni ratios will increase under reducing environments (> 1) and decrease under oxic conditions (< 1) (Tan et al. 2013). V/(V + Ni) ratios less than 0.45 indicate oxic conditions, and > 0.50 indicate anoxic (Liu et al. 2007). The low V/Cr ratios of the collected mudstone samples (0.48–2.15) fluctuate somewhat around the reference material, but all are in the oxic field except one (Fig. 8). The V/Ni ratios of all samples are greater than 1 with an average value of 2.87, indicating a reductive deposition (Fig. 8). The V/(V + Ni) ratios of all samples, ranging 0.57 to 0.83, plot into the anoxic (Fig. 8), and indicate moderate water stratification (Peng et al. 2012; Zheng et al. 2015).
The rare element U in water is oxidized to soluble U6+ under an oxidizing condition and results in the loss of U, while Th is generally present as insoluble Th4+ and stable under redox conditions (Morford et al. 2009). Thus, the U/Th ratio and δU value (δU = 2U/(U + Th/3)) are widely used to reconstruct the paleoredox condition. And U/Th ratios < 0.27 indicate an oxic condition, 0.27–0.50 reflect a dysoxic condition, and > 0.50 reveal an anoxic environment (Wignall and Twitchett 1966). The δU value < 1 indicate an oxic condition, and > 1 reveal an anoxic condition (Tonger et al. 2004; Wang et al. 2017). The U/Th ratios of the Liwaxia to Naijiahe formations are mostly > 0.27 (0.13–3.55, avg., 0.90; 0.77–1.99, avg. 1.50; and 0.31–2.24, avg. 0.92, respectively), indicating a dysoxic-anoxic condition (Fig. 8). This conclusion is also supported by high δU value (mostly > 1) (Fig. 8).
Paleosalinity
Paleosalinity is a significant indicator that is used to reflect the sedimentary environment of the water column in the geologic history (Cheng et al. 2021). The value of 100MgO/Al2O3 can be used as an indicator for the paleosalinity (Zhang 1988; Lin et al. 2020; Stanistreet et al. 2020; Li et al. 2022). The 100MgO/Al2O3 values of all samples ranging 15.45-861.64, suggest a high salinity condition. Generally, Sr is derived from a saline water column, whereas Ba is accumulated on the fine-grained clastic sediments (Wang et al. 2018, and reference therein). Sr/Ba ratio therefore is extensively used for reconstruction of the paleosalinity, and high Sr/Ba ratios (> 1.0) and low ratios (< 0.50) indicate a high-salinity or low-salinity column water, respectively (Meng et al. 2012; Fu et al. 2016; Wang et al. 2017). The relative lower Sr/Ba ratio of the Liwaxia formation (0.12–3.15, avg. 0.74) suggest a low-moderate salinity condition (Fig. 8). In contrast, the high Sr/Ba ratio of the Madongshan and Naijiahe formations (0.80–3.22, avg. 2.16, and 0.31–15.13, avg. 4.71), indicating a high-salinity condition. Upwards through the stratigraphy, the values gradually increase from the Liwaxia to Naijiahe formation, which record the highest paleosalinity values. Thus, this feature depicts the transition from freshwater to brackish water. The high salinity of the Madongshan and Naijiahe formations is also supported by appearance of gypsum crystals. In addition, the discovery of Oncolites, Caddisfly cases in the Madongshan-Naijiahe Formation further confirms high-salinity condition during this period (Zhong et al. 2010; He et al. 2014). Water salinity is closely related to paleoclimate: hot and arid climates commonly have high evaporation rates, resulting in high salinity, whereas warm and humid climates have lower evaporation rates, resulting in lower salinity (Wei et al. 2021). The discriminant parameters of paleoclimate and paleo-salinity show a similar trend (Fig. 8), indicating that the paleoclimate play a significant role in the Liupanshan basin's salinity fluctuation.
Reconstruction of the paleoclimate evolution model of the Liupanshan Basin and its implications
Based on the above geochemical analyses, combined with previous lithology, sedimentary, paleontology, and carbon burial characteristics of the Liwaxia-Naijiahe formations, the paleoclimate evolution model of the Liwaxia-Naijiahe period in the Liupanshan Basin was established (Fig. 9), so as to further discuss the impact of regional paleoclimatic changes on the paleontological and sedimentary evolution in the northern China during the Early Cretaceous.
During the Liwaxia period, the grain size of sediments became finer obviously than early stage, featured by predominant purple-red and gray-green sand-mudstone. The oil shale deposits developed, suggesting that the paleowater depth has obviously increased compared with the early stage. The variable colors of fine-grain sediments imply climatic fluctuation during this period. Previous studies also proposed different views of paleoenvironment, such as high salinity and arid (Zhong et al. 2014), or freshwater and warm environment (Jin et al. 2006; Cai 2021). Based on investigated samples, we prefer a semiarid-arid, anoxic, and low-moderate salinity condition during the Liwaxia period (Fig. 9). The different interpretation of paleoenvironment may be related to the climatic fluctuation, which is resulted by some geologically abrupt event. For instance, the mass mortality of Lycoptera during the Liwaxia period, was caused by a rapid redox change of the water body and toxicity of H2S, which were interpreted as the results of eruption of LIPs and ocean anoxic events (Liang et al. 2022).
During the Madongshan-Naijiahe period, the Liupanshan basin was under an arid, anoxic, and high salinity condition (Fig. 9). The abundant Cypridea-Mongolocypris-Liupanshannia, Lycopera-Huashia, and Kuntulunia- Tongxinichthys in the Madongshan formation appeared, and dominating oil shale deposits developed, suggesting that the paleowater depth has further increased and reached the highest (Fig. 9). Semi-deep lake/deep lake facies deposits are mainly developed in the whole region (Zhao et al. 2020). In addition, the Madongshan period was the main stage of organic carbon burial in the Liupanshan basin. The grey-black oil shale deposits are developed in the stratum, and the TOC abundance reached its highest value (Cai 2021). The sporopollen (e.g., Schizaeoisporites, Ephedripites, and Jugella) and fossil plant taxa, and the leafy shoot morphological and epidermal structures of the present Pseudofrenelopsis in the Naijiahe formation indicate an arid climate as well (Du et al. 2014). The increasing of gypsum sequestration in the Naijiahe Fm. supports a gradual enhancement of salinization (Fig. 9). However, the lake level declined in the Naijiahe period, indicated by sedimentary characteristics (Qu et al. 2003), which is probably caused by extremely dry with high evaporation condition or regional uplift.
The Cretaceous is known to be a time of “hothouse climate” (Wang and Hu 2005; Zhang et al. 2020). In this study, a process of gradual drought and increasing salinity is suggested in the Liupanshan basin. Especially, the arid and anoxic climate conditions during the Madongshan-Naijiahe period probably have a close relationship with the global “hothouse climate” (Fig. 10). On a large scale, the Lower Cretaceous strata are widely distributed in northern China (Cao 2013; Xi et al. 2019), and most Early Cretaceous terrestrial basins display similar sedimentary features and basin evolution (Fig. 10). These basins were mainly filled under the rift basins system and featured by lacustrine sedimentary environments in the Aptian and deposited rich organic shale or coal, indicating a climate favorable for hydrocarbon generation. From NE to NW China, it displays a trend of intensified aridity. Although the development of coal in the NE China indicates a humid climate, the climatic fluctuation also led to the weakening or stopping of coal accumulation (Wang 2018). The enrichment of organic matter in the lake water, featured by the development of black shales, is considered caused by hot and dry continental paleoclimate and lake water salinity under the hothouse climate during the OAE 1 period (Zhang et al. 2021). The geochemical and organic carbon isotope analysis of the lacustrine sedimentary strata in the norther China basins, e.g., Liupanshan, Jiuquan, Liaoning, Jiaolai basins, revealed a good compatibility with the marine sedimentary strata during OAE 1 (Fig. 10) (Yang et al. 2007; Dai et al. 2012; Li et al. 2013; Suarez et al. 2013; Zhang et al. 2016; Zhang et al. 2021b). These evidences further support the hypothesis that OAE 1 has extensive responses in terrestrial basins in northern China. In addition, the gypsum deposits from Naijiahe formation (K1n) of Liupanshan Basin, Zhonggou Formation (K1z) of Jiuquan Basin, and Lianmuqin formation (K1l) of the Junggar and Tuha basins (Cao 2010) are reliable evidence for this regional strong evaporation climate event. A semiarid climate interrupted by arid evaporation alternations was interpreted for the salinization in some lakes in NW China during the Early Cretaceous (Li et al. 2013; Zhang et al. 2021). The high contents of Classopollis were also observed in Lanzhou-Hekou, Yin’e, and Liupanshan basins (Zhang 2011), indicating that a hot, dry, and high-evaporation climate dominated the NW China during the middle-late Early Cretaceous. Therefore, the hothouse climate played a significant role in the development of terrestrial sediments (especially lacustrine hydrocarbon source rocks) during the late Early Cretaceous in northern China. This strengthens understanding of the paleoclimatic patterns of terrestrial lacustrine system in northern China under the background of a greenhouse climate. Also noteworthy is that the hothouse climate pattern was not stable and it was frequently interrupted by short-term cooling events, which has been confirmed by a large number of studies, such as, carbon-oxygen isotope values (Li et al. 2013; Zhang et al. 2021).