Early Jurassic climate and atmospheric CO2 concentration in the Sichuan paleobasin, southwestern China

Climatic oscillations have been developed through the (Early) Jurassic from marine sedimentary archives but remain unclear from terrestrial records. This work presents investigation of climate-sensitive sediments and carbon and oxygen isotope analyses of lacustrine and pedogenic carbonates for the Early Jurassic Ziliujing Formation taken from the Basin in southwestern China. Sedimentary and stable isotope proxies manifest that an overall secular (semi)arid climate dominated the Sichuan Basin during the Early Jurassic, except for the Hettangian. This climate pattern is similar to the arid climate in the Colorado Plateau region in western North America but is distinct from the relatively warm and humid climate in northern China and at high latitudes in the Southern Hemisphere. The estimated atmospheric CO2 concentration (pCO2) from carbon isotopes of pedogenic carbonates shows a range of 980–2610 ppmv (∼ 3.5– 10 times the pre-industrial value) with a mean of 1660 ppmv. Three phases of pCO2 (the Sinemurian 1500–2000 ppmv, the Pliensbachian 1000–1500 ppmv, and the early Toarcian 1094–2610 ppmv) and two events of pCO2 rapidly falling by ∼ 1000–1300 ppmv are observed, illustrating the pCO2 perturbation in the Early Jurassic. The perturbation of pCO2 is compatible with seawater temperature and carbon cycle from the coeval marine sediments, suggesting a positive feedback of climate to pCO2 through the Early Jurassic.

Data from the terrestrial realm also provide important details of environmental and climatic change (e.g., Hesselbo et al., 2000;Suan et al., 2010;Jenkyns, 2010;Philippe et al., 2017), from which the oscillated climate could be observed and revealed as well. Terrestrial proxies, such as flora Published by Copernicus Publications on behalf of the European Geosciences Union. 2056 X. Li et al.: Early Jurassic climate and atmospheric CO 2 concentration (e.g., Riding et al., 2013;Deng et al., 2017;Philippe et al., 2017), vegetation (Pole, 2009), and geochemistry (e.g., Riding et al., 2013;Kenny, 2015;Tramoy et al., 2016), as well as the atmospheric CO 2 concentration (pCO 2 ) record (e.g., Retallack, 2001a;Beerling and Royer, 2002;McElwain et al., 2005;Berner, 2006;Steinthorsdottir and Vajda, 2015), provide an emerging record of the Early Jurassic terrestrial climate and environment changes. Correspondingly, the proxy application of terrestrial sedimentary archives could play a key role in the global Early Jurassic correlation of the marine and terrestrial climate.
Proxies for pCO 2 are the important linkage between the marine and terrestrial climatic condition. Studies of the terrestrial pCO 2 record have focused on the Triassic-Jurassic boundary (e.g., Tanner et al., 2001;Cleveland et al., 2008;Schaller et al., 2011;Steinthorsdottir and Vajda, 2015) and the Toarcian oceanic anoxic event (McElwain et al., 2005), where pCO 2 estimates range from 1000 ppm to 4000 ppmv (e.g., Tanner et al., 2001;Cleveland et al., 2008;Schaller et al., 2011). Few relatively continuous pCO 2 records and coupled terrestrial climate changes have been documented for the Early Jurassic.
There are several large Triassic-Jurassic terrestrial basins in western China, in which the Sichuan Basin has a relatively complete and continuous continental sedimentary sequence of the Upper Triassic-Paleogene (e.g., SBGM, 1991SBGM, , 1997Wang et al., 2010). During the Early Jurassic, the Sichuan Basin was in the Boreotropical climate zone, a suggestion based on climate-sensitive sediments ( Fig. 1a; Boucot et al., 2013), or a warm and temperate climate, suggested based on clay mineralogy and phytogeography (e.g., Dera et al., 2009). In this work, we present a field investigation, including lithofacies and paleosol interpretation, and carbon and oxygen isotope analyses of both lacustrine and pedogenic carbonates in the Sichuan Basin. New results allow us to reconstruct the paleoclimate and relatively consecutive pCO 2 record through the Early Jurassic, which we use to compare stable isotopes of marine sediments and estimated sea water temperature.

Geological setting and stratigraphy
Southwestern China, including the provinces of Yunnan, Sichuan, Chongqing, and Guizhou, had been the main part of the upper Yangtze Plate since the Proterozoic and possibly since the Neoarchean. With the amalgamation of the Cathaysia and Yangtze plates, it became the western South China plate or cratonic basin since the Neoproterozoic (Sinian) and continued to the late Middle Triassic. Via the Indosinian orogeny, new foreland basins were formed since the Late Triassic (e.g., He and Liao, 1985;, recording the Mesozoic and Cenozoic evolution of tectonics, environment, and climate in southwestern China. The Mesozoic Sichuan paleobasin was confined by the Longmenshan thrust belt in the northwest, the Micangshan-Dabashan arcuate thrust belt in the northeast (Fig. 1b), and the northern hilly topography boundary of the Yunnan-Guizhou plateau in the south and east. It was mainly developed during the Late Triassic-Jurassic and includes provincial areas of eastern Sichuan, northern Guizhou, western Hubei, northwestern Hunan, and the entirety of Chongqing. This Triassic-Jurassic Sichuan foreland basin was much larger than the present Sichuan Basin in eastern Sichuan province. We estimate the size of the Sichuan paleobasin is roughly 480 000 km 2 based on lithofacies paleogeography (Fig. 1b;Ma et al., 2009;Li and He, 2014) and suggest naming this the grand Sichuan paleobasin (GSB).
The Mesozoic terrestrial sediments accumulated up to ∼ 9 km (Guo et al., 1996) in the GSB, and the Jurassic part can be up to 3-3.5 km thick (SBGM, 1991). Two types of Lower Jurassic deposits have been distinguished (Table 1): the Baitianba Formation (Fm) in the north (∼ 10 %) and the Ziliujing Fm (e.g., SBGM, 1991;Wang et al., 2010) in the south (over 90 % of the basin).
The Baitianba Fm was deposited unconformably on the Upper Triassic Xujiahe Fm and is overlain conformably by the Middle Jurassic Xintiangou Fm and Qianfoyan Fm (Table 1). It is mainly composed of grayish shales and sandstones with coal layers and massive conglomerates. Abundant plant fossils, sporopollens, conchostracans, bivalves, and gastropods indicate it is of the Early Jurassic (SBGM, 1991(SBGM, , 1997. Sporopollen assemblages of the Hettangian-Sinemurian age were found in the lower part (Zhang and Meng, 1987) and the Pliensbachian-Toarcian assemblages were reported in the upper part (Wang et al., 2010).
The Da'anzhai Member is characterized by dark gray to black shales and bioclastic limestones with a southward increase of reddish mudrocks (SBGM, 1991(SBGM, , 1997Wang et al., 2010) and is regarded as the sediment in a grand Sichuan paleolake (e.g., Ma et al., 2009;Li and He, 2014). Ostracod assemblages indicate it is the late Early Jurassic (e.g., Wei, 1982;Wang et al., 2010). A Re-Os isochron age of 180.3 ± 3.2 Ma associated with an organic carbon isotope  (Boucot et al., 2013) laid over an Early Jurassic (∼ 193 Ma, Sinemurian) paleogeographic map (Scotese, 2014). (b) Lithofacies paleogeographic sketch of the grand Sichuan paleobasin (GSB) in the early Early Jurassic (Zhenzhuchong and Dongyuemiao members) showing locations of the observed and analyzed sections and climate-sensitive sediments. The lithofacies paleogeographic map was composed and modified from Ma et al. (2009) and Li and He (2014). The blue area is the extent of the paleolake, estimated as ∼ 380 000 km 2 . The blue-gray region is the basin shape, estimated as ∼ 480 000 km 2 . The dotted red line encircles the deeper lake area in the late Early Jurassic (Ma'anshan and Da'anzhai members). The black dashed line is the northern edge of Calcisol occurrence, which may separate the climate of the GSB into the northern and southern types. Triangles with numbers are locations of observed and analyzed sections: A1, Xiasi section, Jian'ge; A2, Puji section, Wangcang; A3, Shiguansi section, Wanyuan; A4, Shaping section, Ya'an (bed and thickness from Wen and Zhao, 2010); A6, Tanba and Maliuping section, Hechuan (bed and thickness from Wang et al., 2010); A7, Wenquan section, Kaixian (thickness from Wang et al., 2010). The location and source data of sections A5 and A8-A17 (climate-sensitive sediments) refer to Supplement Table S1. 2058 X. Li et al.: Early Jurassic climate and atmospheric CO 2 concentration Table 1. Stratigraphic framework of the Lower Jurassic Ziliujing Fm in Sichuan and adjacent area (GSB), southwestern China.
The Ma'anshan Member is comprised of violet and red mudrocks with a few greyish and greenish thin-bedded fine sandstones and siltstones in which floral fossils are common . The Dongyuemiao Member consists of greenish and reddish mudrocks and siltstones with greyish bioclastic limestone and marlstone, of which abundant bivalve and plant fossils were reported from eastern Sichuan and Chongqing Meng et al., 2003;Wang et al., 2010). The Zhenzhuchong Member is dominated by violet red mudrocks and shales intercalated with thin-bedded sandstones and/or siltstones and numerous plant fossils of the Early Jurassic affinity (e.g., Duan and Chen, 1982;Ye et al., 1986). Taken together, fossil associations suggest that the three members were deposited in the middle to late Early Jurassic. The age limitation of the overlying Da'anzhai Member and the correlation to the Lufeng dinosaur fauna places these members in the Sinemurian-Pliensbachian, and the Zhenzhuchong and Dongyuemiao Fms are suggested to be from the Sinemurian (Table 1).
The Qijiang Member is composed of quartz arenite interbedded and intercalated with dark shales. Coal seams are often seen in the middle of the Qijiang Member. This member mainly occurs in the central part of the GSB. It is likely the earliest Jurassic, possibly Hettangian age, but plant fossils cannot precisely indicate the age (Wang et al., 2010).

Materials and methods
We have measured sections and made detailed observations and descriptions of sedimentary characteristics for lithofacies analysis at six outcrop sections (Locations A1-A4, A6, and A7, Fig. 1). Published descriptions for other sections (Locations A5, A8, and A9, Fig. 1) are integrated into our observations. Details of microscopic examination of sedimentary rocks and analysis of sedimentary facies underpinning the climate analysis are attached as the Supplement Note S1. Below we state climate-sensitive sediment observations, carbon and oxygen isotope analyses, and estimates of pCO 2 .

Observation of climate-sensitive sediments
Climate-sensitive sediments are mainly dolomites, gypsum, and paleosols, which are used to analyze the climate in this work (Table S1).
Dolomites and gypsum are relatively easy to recognize both in the field and under a microscope. We distinguish dolomites from limestones following Tucker (2011) andFlügel (2004). As Flügel (2004) stated, field distinctions of limestone and dolomite can also be made, although detailed differentiation of carbonate rocks is best performed in the laboratory. In the field, we recognize gypsum via particular structures, such as a chicken-wire cage, gypsum pseudomorph, and clusters of (0.5-1 cm) pore.
There are multiple classifications of paleosols (e.g., Wright, 1992;Mack et al., 1993;Retallack, 2001b;Imbellone, 2011), mostly based on the US Soil Taxonomy. We recognized paleosols in the field based on color, structures, horizonation, root traces, and textures, and followed the general classification paleosols by Mack et al. (1993) and Retallack (2001b). In this paper, paleosols are described following the procedures of the Soil Survey Manual and classified according to Soil Survey Staff (1998).
Within the measured and observed sections, paleosol profiles were mainly identified from the two main sections A4 and A6 (Figs. S1 and S2, and Table S2). Horizonation, Bk horizon thickness, boundaries, structures, trace fossils, rootlets, carbonate accumulations (calcretes), etc. were recorded (Table S2). Paleosols interpreted in other cited sections ( Fig. 1) rely on the description of lithology, structure, and calcrete in the original references. Based upon a modification of the Retallack (1998) categorization of paleosol maturity, the relative paleosol development (maturity) was assigned.

Analyses of carbon and oxygen isotopes
The 10 lacustrine carbonate samples were analyzed for carbon and oxygen isotopes from the Da'anzhai Member at the Shaping section, Ya'an (Location A4; Fig. S1 and Table S3). A total of 26 pedogenic carbonate samples were analyzed for carbon and oxygen isotopes from 32 paleosols of the Ziliujing Fm at the same section ( Fig. S1 and Table S4). Two or three microdrilling powder samples (columns 7 and 8 in Table S4) were taken from the same individual calcrete for stable isotope analysis, and then a mean value for each calcrete sample was calculated (columns 9 and 10 in Table S4).
At the field scale, calcretes are ginger-like and sporadically spaced within the soil horizon. We observed no linear and planar calcretes that would indicate precipitation at or below the water table. Before drilling, thin sections were petrographically studied using polarized light microscopy and cathodoluminescence imaging. Micritic calcite is predominant in both lacustrine and pedogenic carbonate samples, with no evidence for carbonate detritus in calcretes ( Fig. 2a and b). The micritic calcites used for stable isotope analyses are chiefly null-to non-luminescent, with <10 % light orange and brownish luminescence, indicating genesis primarily in the vadose zone. While luminescent calcretes indicate a high possibility of hydrological influence (e.g., Mintz, et al., 2011), we sampled to avoid this. Based on petrography and cathodoluminescence (CL) imaging together with the field observations, the dense micritic zones sampled for the stable isotope composition should give pristine δ 13 C values that can be used to estimate pCO 2 .
Microsampling of lacustrine and pedogenic carbonates focused on only micrites, avoiding diagenetic spar from cracks, veins, and vug spaces. Powder samples were obtained using a dental drill (aiguille diameter φ = 1-2 mm).
Isotopic analyses were conducted on 0.3-0.5 mg powder samples. Powder samples were dried in an oven at 60 • C for 10 h before being moved to the instrument. Carbon dioxide for isotopic analysis was released using orthophosphoric acid at 70 • C and analyzed online in a DELTA-Plus XP (CF-IRMS) mass spectrometer at the State Key Laboratory for Mineral Deposits Research, Nanjing University. The precision of the measurements was regularly checked with a Chinese national carbonate standard (GBW04405) and the international standard (NBS19) and the standard deviation of δ 13 C was ± 0.1 ‰ over the period of analysis. Calibration to the international Pee Dee Belemnite (PDB) scale was performed using NBS19 and NBS18 standards.

Calculation of atmospheric CO 2 concentration
The Cerling (1991Cerling ( , 1999 equation was used to calculate the pCO 2 using the carbon isotope of pedogenic carbonates as below: where C a is pCO 2 ; δ 13 C s , δ 13 C r , and δ 13 C a are the isotopic compositions (‰) of soil CO 2 , soil-respired CO 2 , and atmospheric CO 2 , respectively; and S (z) is the CO 2 contributed by soil respiration (ppmv).
The value of δ 13 C r represents the carbon isotope ratio of average bulk C3 vascular tissue (Arens et al., 2000), reflecting atmospheric δ 13 C (Jahren et al., 2008). The δ 13 C om of organic matter within paleosols based on the range of modern C3 ecosystem fractionations (Buchmann et al., 1998;Ekart et al., 1999) is commonly used for δ 13 C r . However, the δ 13 C r could be compromised in fossil soils due to oxidation and metabolism of organic matter after burial (Nadelhofer and Fry, 1988). In this paper, we use the δ 13 C om from the Paris Basin Peti et al., 2017) for the Sinemurian-Pliensbachian δ 13 C r and from Cardigan Bay, UK (Xu et al., 2018), for the Toarcian.
The δ 13 C a value, the carbon isotopic composition of the atmosphere, was about −8 ‰ in the 1980s, being depleted relative to the pre-industrial atmosphere, which was around −6.5 ‰ (Friedli et al., 1986). The average value of −6.5 ‰ has been chosen as the δ 13 C a for acquiring δ 13 C r and S (z) (e.g., Ekart et al., 1999;Robinson et al., 2002), and the δ 13 C a was generally calibrated as δ 13 C ac from δ 13 C r using the formula (δ 13 C r + 18.67)/1.1 (Arens et al., 2000). Herein we used both calibrations to calculate the δ 13 C a (Table S4).
S (z) is the largest source of uncertainty in pCO 2 estimates (Breecker, 2013) and the uncertainty arises primarily from the sensitivity to soil-respired CO 2 (S (z) (Montañez, 2013). It is a function of depth and effectively constant below 50 cm (e.g., Cerling, 1991). S (z) = 2500 ppmv is suggested for the subhumid temperate and tropical climates (Breecker et al., 2010), 2500-5000 ppmv is suggested for higher moisture and productivity soil (Montañez, 2013), 2000 ppmv is suggested for semi-arid areas (Breecker et al., 2009), 1500-2000 ppmv is suggested for aridisols and alfisols (calcisolargillisol), and 2000 ± 1000 is suggested for paleo-vertisol (Montañez, 2013), with 1000 ppmv in desert areas (Breecker et al., 2010) or 400 ± 200 ppmv for immature soil (Montañez, 2013). In this context, we chose the S (z) = 2000 ppmv for calculating pCO 2 at 25 • C as the calcisols are reddishbrownish aridisols, and we also compared the results with that by S (z) = 2500 ppmv (Table S4). Additionally, we took samples at the middle and lower Bk horizon (often > ∼ 20-30 cm to the Bk top). This means that the depth of calcrete samples in the examined paleosols was generally deeper than 50 cm below the paleosol surface, meeting the requirement for a constant value of S (z) .

Results
Based on the investigation of cross sections (locations A1-A4, and A6-A7; Fig. 1), we have classified six sedimentary facies units in the Ziliujing Fm. They are alluvial fan, fluvial river, floodplain, lake, lake delta, and swamp facies. Details of description and interpretation are in the Supplement Note S1. Below are results of climate-sensitive sediment observation, stable isotope analyses, and pCO 2 calculation.

Climate-sensitive sediments
Field observation combined with published calcrete materials shows that paleosols widely occur in the Lower Jurassic Ziliujing Fm of the GSB (Figs. 1, 3, and 4). A total of 32 paleosols were observed and described at the Shaping section, Ya'an, and five paleosols were found at the Tanba section, Hechuan (Table S2).
Dolomites were found at seven locations in the central and southern GSB (Figs. 1, 4, and Table S1). The dolomites chiefly occur in the Toarcian Da'anzhai Member and a few in the Sinemurian-Pliensbachian Dongyuemiao and Ma'anshan members (Fig. 4). They are often massive whitish (Figs. 3f and S3e) and micritic ( Fig. S4b and d), likely indicating an authigenic origin.

Carbon and oxygen isotope values
The δ 13 C values of lacustrine carbonate samples range from −2.02 ‰ to −4.07 ‰, and δ 18 O values range from −9.91 ‰ to −12.28 ‰ (Table S3 and Fig. 5). An increasing trend of both carbon and oxygen isotope ratios is observed from lower to upper horizons across a 45 m stratal interval of the lower Da'anzhai Member (Fig. 6).
Pedogenic carbonate samples have δ 13 C values from −3.52 ‰ to −8.10 ‰, which fall in the typical stable isotope range for pedogenic carbonates. Values of −6 ‰ to −8.0 ‰ characterize the sequence of the Zhenzhuchong Member and the main Ma'anshan Member, with an abrupt increase to −5.5 ‰ to −3.5 ‰ at the top of the Ma'anshan Member (samples J1z-16-01 and J1z-18-01; Fig. 6). The δ 18 O values are mainly from −11.3 ‰ to −13.10 ‰ in the interval of the Zhenzhuchong Member and Ma'anshan Member. The δ 18 O follows δ 13 C with a sudden increase to −5.5 ‰ at the top of the Ma'anshan Member (Fig. 6). Large and frequent variations of both carbon and oxygen isotope ratios can be observed in the lower Da'anzhai Member ( Fig. 6 and Table S4).

CO 2 concentrations
pCO 2 values based on paleobarometer modeling of paleosol calcite (Cerling, 1999) vary depending on the parameters used for the calculation.
Results show that pCO 2 values at S (z) = 2500 ppmv are larger than at S (z) = 2000 ppmv and that the S (z) is the largest uncertainty leading to different pCO 2 values. The highest difference (columns 16 and 18 in Table S4) of pCO 2 is ∼ 650 ppmv (3260 − 2610 ppmv) (Sample J1z-22-01 at depth 324.5 m), while the lowest difference is ∼ 245 ppmv (1226−981 ppmv) (Sample J1z-12-01 at depth 148.9 m), and the mean is ∼ 370 ppmv. In addition, when S (z) is the same, the pCO 2 values are close even if other parameters are different (comparison between columns 15 and 16 and 17 and 18 in Table S4; Fig. 6).
Whatever parameters are used, the trend of pCO 2 over the epoch is quite similar (Fig. 6). We chose S (z) = 2000 ppmv (column 18 in Table S4) to illustrate the nature of the Early Jurassic pCO 2 in the GSB.
pCO 2 values mostly range between 980 and 2610 ppmv, and the mean 1660 ppmv is ∼ 6 times the pre-industrial 275 ppmv. Most of the pCO 2 values are 1000-2000 ppmv, with the mean 1580 ppmv in the Zhenzhuchong and Figure 6. Diagram of the Lower Jurassic strata and lithological log at the Shaping section, Ya'an, with carbon and oxygen isotope values of pedogenic and lacustrine carbonates and the pCO 2 curve. Three phases and two events can be observed for both stable isotope values of pedogenic carbonates and pCO 2 estimate. A legend of lithology in log refers to Supplement Figs. S1 and S2. T-OAE stands for the Toarcian oceanic anoxic event. The 1E CO 2 and 2E CO 2 are events of rapid falling pCO 2. Numbers 15 to 18 are the curves of pCO 2 in different parameters, and details refer to Supplement Table S4. Discrepancies (errors) are produced from the pCO 2 subtraction of column 16 (S (z) =2500 ppmv) from column 18 (S (z) = 2000 ppmv) in Table S4, indicating that the largest uncertainty for the estimate of pCO 2 is the S (z) . The highest difference (columns 16 and 18 in Table S4) of pCO 2 is 652 ppmv (3262 − 2610 ppmv) (Sample J1z-22-01 at depth 324.5 m), the lowest difference is 245 ppmv (1226 − 981 ppmv) (Sample J1z-12-01 at depth 148.9 m), and the mean is ∼ 370 ppmv.
It is noted that the Gaussian errors of pCO 2 range from 384 to 1017 ppmv with a mean 647 ppmv (Table S5), leading to a large uncertainty of the mean ∼ 39 %. The largest source of the uncertainty is the standard error (766 ppmv) of modern soil carbonate (Breecker and Retallack, 2014). The pCO 2 uncertainty decreases by ∼ 20 % if half (383 ppmv) of the standard error of soil carbonate is selected, and decreases to ∼ 12 % if 0.25 (∼ 191 ppmv) standard error is used. The second largest source of error in the pCO 2 is the S (z) estimate. The uncertainty of pCO 2 becomes much smaller when the S (z) is larger, e.g., it will fall from ∼ 39 % to ∼ 17 % if S (z) = 5000 ppmv instead of 2000 ppmv. Other parameters such as temperature, δ 13 C r , δ 13 C a , and δ 13 C s contribute very little to the calculated pCO 2 uncertainty. The uncertainty of pCO 2 is the same when δ 13 C a is determined by the transfer from δ 13 C c of marine planktonic fossil carbonates (Table S8) or from organic matter (Table S5).

Discussion
The Jurassic marine record shows climatic and environmental oscillations (e.g., van de Schootbrugge et al., 2005;Dera et al., 2011;Gómez et al., 2016;Arabas et al., 2017), including sea water temperature fluctuation and carbon cycle reorganization recorded in both carbonate and organic matter. The climate changes and events recorded in the marine realm have been mainly attributed to Karoo-Ferrar volcanism (e.g., Hesselbo et al., 2000;Caruthers et a., 2013), sea level change (e.g., Hesselbo and Jenkyns, 1998;Hallam and Wignall, 1999), orbital forcing (e.g., Kemp et al., 2005;Huang and Hesselbo, 2014;Storm et al., 2020), and/or the opening of the Hispanic corridor (e.g., van de Schootbrugge et al., 2005;Arias, 2009). Eruption of the Karoo-Ferrar and Central Atlantic magma is thought to have released large amounts of CO 2 into the atmosphere in a short amount of time, resulting in rising temperatures of both marine and continental realms. The nearly continuous record of the Jurassic strata in the GSB provides an excellent test of this hypothesis in the terrestrial realm. We compare the climate and pCO 2 record from the GSB in relationship to the marine temperature records.

Paleoclimate variation
During the Late Triassic, southwestern China was warm to hot and humid and occupied a tropical and/or subtropical zone, as demonstrated by palynoflora, coals, and perennial riverine and lacustrine lithofacies in the Xujiahe Fm (e.g., Huang, 1995;. However, the climate became dry through the Early Jurassic, manifesting as climatesensitive sediments and stable isotopes, although there are two lithofacies packages reflecting two major lake stages (for details refer to the Supplement Note S1) in the GSB.

The Hettangian age
In the Hettangian, the climate was warm and humid like the Late Triassic in the GSB. The Qijiang Member is comprised of mainly mature quartz arenites and siltstones with coals ( Fig. 7), as well as siderite concretions, indicating a stable tectonic setting and warm and humid climate in the eastern and southern GSB. Climate was similar across the whole region, because multiple coal layers occur in the lower Baitianba Fm. The alluvial fan system of the lower Baitianba Fm (Figs. 7 and S6) is characterized by moderate to good roundness and sorting of gravels with sandy matrix (Fig. S3a; e.g., Liu et al., 2016;Qian et al., 2016;and this work). In the Newark basin of eastern North America, climate-sensitive sediments such as nodules of carbonate and gypsum (pseudomorph), as well as mudcrack in mudflat facies, indicate an arid climate in the fifth cycle of the Hettangian (> 199 Ma) Passaic Fm (Kent et al., 2017). The more widespread eolian Navajo Sandstone, dated as Hettangian-Sinemurian (200-195 Ma;Parrish et al., 2019), indicates arid in the Colorado Plateau ( Fig. 1a; Boucot et al., 2013).

The Sinemurian age
The early Sinemurian Zhenzhuchong Member is a combination of riverine floodplain and lacustrine facies (Supplement Note S1). The lithology is dominated by violet red mudrocks with few thin greyish and greenish fine sandstones and siltstones. The reddish color of the rocks may indicate a change of climate. Differences in the color appearance show that the reddish color started in the middle member in the central basin (Location A6; Fig. S2) but almost developed through the whole member in the western basin (Location A4; Fig. 6).
This climate change, interpreted from reddish mudrocks and paleosols, is consistent with the floral fossils (e.g., Huang, 2001;Wang et al., 2010), suggesting the decreasing humidity and increasing temperature from the Late Triassic epoch and the Hettangian age into the Sinemurian age in the southern GSB. However, in the northern GSB there are few proxies for climate change, and alluvial fan and lacustrine delta facies common in the middle Baitianba Fm (Fig. S6; e.g., Qian et al., 2016) do not give us information on climate.
The Sinemurian (semi)arid climate interpreted from reddish mudrocks and calcisols is supported by the flora  and the mudrock geochemistry (Guo et al., 2017). Few records of coeval terrestrial climate are documented from other continents or regions in the literature. The Whitmore Point Member of the Moenave Fm deposited in dryland lakes (Tanner and Lucas, 2008) and the upper part of eolian Navajo Sandstone (Blakey et al., 1988) could represent the coeval strata with sediments of similar climate in the Colorado Plateau. It might be an interruption of relatively cool (∼ 9 to 18 • C) continental climate inferred from oxygen and hydrogen isotope composition of chert precipitated in inter-dune, freshwater lakes in the Navajo Sandstone (Kenny, 2015). In eastern England, the co-occurrence of the acmes of thermophilic pollens Classopollis classoides and Liasidium variabile indicates a warm and humid climate in the late Sinemurian (Riding et al., 2013), likely illustrating a different climate scenario in Europe.

The Pliensbachian age
The Ma'anshan Member of the Pliensbachian displays a prominent change in the distribution and extent of red color sediment and pedogenesis. The reddish sediments extend through the entire member (cf. Figs. 6 and S2) and can be ob- served across most of the GSB. Calcisols are documented in both the western and central GSB (Figs. 6, 7, S1, and S2). A total of 10 calcisol horizons were recognized at the Shaping section, Ya'an (Figs. 6 and S1). Strongly leached pedogenic structures and mudcracks are seen in Bed H8 of the Tanba section, Hechuan ( Fig. 3a and d). Abundant calcretes within terrestrial red mudrocks are widely described at Gaoxian of Dafang (Location A8; Zhang et al., 2016), Hulukou of Weiyuan (Location A10; SBG, 1980), Geyaoguan of Gulin (Location A13; SBG, 1976), Taiyuan of Fengdu (Location A16;SBG, 1975), and Yaxi of Zunyi (Location A17; Yang, 2015). The widespread distribution of redbeds and calcisols (Figs. 4 and 7) denotes an intensification of the (semi)arid climate.
Plant and sporopollen fossils also show a change to drier climate in the Pliensbachian. Compared to the Sinemurian members, more plant fossils are reported in this member (e.g., Meng and Chen, 1997;Wang et al., 2010). The Pliensbachian-Toarcian sporopollen assemblages are dominated by a sporomorph genera assemblage of Dictyophyllidites, Cyathidites, and Classopollis, in which the dry-type gymnosperm spore Classopollis is more prevalent than in the Hettangian-Sinemurian (Zhang and Meng, 1987).
A similar dry temperate or subtropical climate is interpreted for the upland coniferous forest in Qaidam Basin in northwestern China (Wang et al., 2005) and by inter-dune playa mudstones of the Kayenta Fm in the Colorado Plateau (e.g., Bromley, 1992), although there was a cool and humid climate in southern Kazakhstan, central Asia (Tramoy et al., 2016).

The Toarcian age
In spite the fact that the Da'anzhai Member was deposited in the largest lacustrine transgression period ( Fig. 7; for details see Supplement Note S1), abundant evidence for arid conditions, including backshore reddish mudrocks with calcisols, lacustrine micritic dolomites and/or gypsum, and stable isotopic geochemistry of lacustrine carbonate, indicate that the Toarcian aridification could be the most intensive of the late Early Jurassic in the GSB.
Gypsum and micritic dolomites are reported in the western and southern GSB (SBG, 1980;Mo and Yu, 1987;Peng, 2009; and this work) (Figs. 1, 4, and 7). Though there are a number of hypotheses on the dolomite formation in deep time, such as authigenic origin, diagenetic replacement, microbial mediation (e.g., Vasconcelos et al., 1995;Mckenzie and Asconcelos, 2009;Petrash et al., 2017), a high abundance of dolomite was interpreted as forming during greenhouse periods, characterized by warm climates, probably reflecting favorable conditions for evaporite deposition and dolomitization via hypersaline reflux (Warren, 2000). Dolomites are also thought to be the result of interplay of climate with sea level and base level change (e.g., Newport et al., 2017) or of interaction with climatic regimes ( Vandeginste et al., 2012). The widespread micritic dolomites in the Da'anzhai Member, which are associated with gypsum (Fig. 3f), likely indicate an arid climate in the central and western GSB (Fig. 1b). Gypsum occasionally occurs at Maliuping of Hechuan (Fig. 3f) and Wujiaba of Zigong (SBG, 1980), showing a possible evaporitic climate in the early Toarcian in the central GSB.
Carbon and oxygen isotopes of lacustrine carbonates further support the interpretation of an arid climate in the Toarcian age in the GSB. The mainly positive δ 13 C values, 0 to 2 ‰ (Fig. 5), from Hechuan  indicate the lakes were brackish or even saline. The relatively heavy negative δ 13 C values, −1 ‰ to −3.5 ‰ (Fig. 5), from Zigong  and Ya'an (this work) denote low depletions of 13 C during calcite and aragonite precipitation and mean that the lakes were possibly brackish. Lightly negative δ 18 O values, −5 ‰ to −12 ‰ (Fig. 5), of the lacustrine carbonates suggest closed lacustrine, palustrine, and pond systems formed in a regional arid or semiarid climate with evaporation exceeding precipitation.
The covariance of δ 13 C and δ 18 O is a criterion to distinguish closed or open lakes (e.g., Talbot, 1990;Li and Ku, 1997). Pronounced positive covariances (R 2 = 0.44-0.96) between carbon and oxygen isotopes (Fig. 5) indicate a typical arid-semiarid pattern of lakes in the central and western GSB.
The Da'anzhai Member has the same palynofloral assemblage with the Ma'anshan Member, in which the dry-type gymnosperm spore Classopollis is more abundant than in underlying strata (e.g., Zhang and Meng, 1987;Wang et al., 2010), supporting the aridification indicated by climatesensitive sediments and stable isotope ratios of lacustrine carbonates mentioned above.
Coastal Cheirolepidiacean (gymnosperm) forests indicate (temperate to subtropical) warm and humid climate punctuated by locally dry and/or arid events in the Toarcian in Qaidam Basin in northwestern China (Wang et al., 2005). In Inner Mongolia in northern China, thermophilous plants, such as the dipteridaceous fern Hausmannia, Bennettitales Ptilophyllum, display similar warm and humid climate interrupted by hot and even arid conditions in a short intervals of the Toarcian (Deng et al., 2017). The warm and wet climate was also indicated by assemblages of sporomorph and vegetation in the late Early Jurassic in Jurong of Jiangsu, in the Lower Yangtze area . In southern Kazakhstan, central Asia, paleoflora and δ 2 H values suggest slightly less humid and warmer conditions starting from the early Toarcian (Tramoy et al., 2016).
Climate-sensitive sediments, carbon and oxygen isotope values and covariance, and palynoflora together indicate that an overall (semi)arid climate dominated the GSB during the Early Jurassic, possibly accompanied by occasional evaporitic climate. Relatively abundant calcisols suggest that the GSB was in a subtropical arid zone based on the paleoclimatic zonation model of paleosols (Mack and James, 1994) during the middle-late Early Jurassic. Through the Early Jurassic, this (semi)arid climate in GSB is thoroughly comparable with the simultaneous arid climate recorded in dryland lacustrine and eolian facies in the Colorado Plateau (e.g., Blakey et al., 1988;Bromley, 1992;Tanner and Lucas, 2008;Parrish et al., 2017) but distinct from the relatively warm and humid climate indicated by sedimentological and floral characteristics in northern China (e.g., Wang et al., 2005;Deng et al., 2017) and in relatively high latitudes of the Southern Hemisphere (Pole, 2009).
In summary, the increasing aridity and warming in the GSB and arid climate in the Colorado Plateau could have been consecutive through the Early Jurassic and seems to not be harmonized with the global fluctuated climate that could be imprinted by two large volcanic eruptions of the Central Atlantic magmatic province and Karro-Ferrar large igneous province. The secular arid climate in the two areas is more probably constrained by paleotopography, where both were laid in the relatively low latitudes 15-30 • N (Fig. 1a).

pCO 2 perturbations and events
Pedogenic carbonates found in various continental settings precipitate in direct contact with soil atmosphere and bedrock and hold a meaningful signature of past climate (Alonso-Zarza and Tanner, 2006). There are few high age resolution pCO 2 reconstructions for the Early Jurassic. The pCO 2 estimates have a focus on the event horizons, such as the transition from the Triassic to Jurassic (e.g., Tanner et al., 2001;Schaller et al., 2011). Herein we present a pCO 2 estimate based on data from the GSB at ∼ 1.0 Myr age resolution for a ∼ 20 Myr  interval of the Early Jurassic (Figs. 6 and 8a).

pCO 2 perturbation
Results of model estimates show that the pCO 2 values range from 980 to 2610 ppmv with a mean 1660 ppmv in  Tables S6 and S7). Vertical bars are errors (1σ ) of pCO 2 (Table S5). Errors are propagated using the Gaussian approach (Breecker and Retallack, 2014). Note that (1) pCO 2 = 4027 ppmv (black solid diamond, sample J1z-20-01) if the δ 13 C r = −29.0 ‰ at 181 Ma, taken from Xu et al. (2018), in the case of other constant parameters. (2) The early published pCO 2 values from both carbon isotope of pedogenic carbonates and stomatal index of fossil plants (data refer to Tables S6 and S7) were very roughly dated using the average age of a lithostratigraphic formation or group, with which the uncertainty can be up to 10 Myr, leading to the difficulty of precise and accurate pCO 2 correlation in pace, frequency, and event. (b) The δ 18 O and seawater temperature (dotted black line) of marine invertebrate fossils compiled from Rosales et al. (2001), Jenkyns et al. (2002, Bailey et al. (2003), van de Schootbrugge et al. (2005), Gómez et al. (2008), Metodiev and Koleva-Rekalova (2008), Suan et al. (2008), Korte et al. (2009), Dera et al. (2011), and Gómez et al. (2016. (c) The dotted red line shows δ 13 C of marine carbonates in the western Tethys Ocean, composed from Clayton (1986, 1997), Hesselbo et al. (2000), Dera et al. (2011), Arabas et al., 2017; the dotted and solid black line δ 13 C of organic matters from Paris Basin, France (Peti et al., 2017). Smoothed δ 18 O and seawater temperature (red curves) in (b) and (c) are following Dera et al. (2011). (d) The δ 13 C of organic matter from the North Atlantic taken from the Mochras borehole, Cardigan Bay Basin, UK (Xu et al., 2018;Storm et al., 2020), with a seven-point average smoothing against depth (mbs).
The evolution and level of pCO 2 estimated by carbon isotope ratios of the pedogenic carbonates from the GSB compare favorably with the global composite based on the plant stomata method (for data of the composite curve, see Table S6) but show significant differences relative to the global composite pCO 2 based on paleosols ( Fig. 8a; Suchechi et al., 1988;Cerling, 1991;Ekart et al., 1999), which may be attributed to the shortage (<4 samples) of global data and large age uncertainties ( Fig. 8a and Tables S5 and S6). The changes in pCO 2 from the GSB have a similar pattern to coeval seawater temperature estimates through the Early Jurassic, although there are some discrepancies in pace and in detail (cf. Fig. 8a and b). Therefore, the relatively high pCO 2 1500-2000 ppmv approximately corresponds to the relatively high seawater mean temperature −2 to +2 • C in the Sinemurian, low pCO 2 1000-1500 ppmv corresponds to low seawater mean temperature −5 to −2 • C in the Pliensbachian, and quick rising pCO 2 of 1200 ppmv to ∼ 2500 ppmv corresponds to the rapidly increased seawater temperature of −4 to +4 • C in the late Pliensbachian-early Toarcian.
The pCO 2 record roughly trends with the carbon isotope records of marine carbonates and organic matter (cf. Fig. 8a to c and d), suggesting a possible linkage of the pCO 2 record in the GSB to the global carbon cycle (see Sect. 5.2.2). Nevertheless, it is difficult for the proxies to compare in a higher detail, making it difficult to relate the record to orbital forcing of the global carbon cycle in the Sinemurian-Pliensbachian (Storm et al., 2020).
As a greenhouse gas, atmospheric CO 2 has a strong control over global temperatures for much of the Phanerozoic (e.g., Crowley and Berner, 2001;Royer, 2006;Price et al., 2013), but a decoupling of CO 2 and temperature has also been suggested (e.g., Veizer et al., 2000;Dera et al., 2011;Schaller et al., 2011). The pattern of the Early Jurassic pCO 2 reconstructed from the carbon isotope of pedogenic carbonates in the GSB in southwestern China supports the coupled relationship of CO 2 with temperature. Models of the coupling and decoupling of CO 2 -temperature and the CO 2carbon cycle have to consider (1) the age order of CO 2temperature and carbon cycle relevance, i.e. they should be related in the same age (long-term or short-term) hierarchy; (2) precise age constraints of individual CO 2 and temperature data; (3) methods of CO 2 and temperature estimates, depending on precondition, presumptions, parameters, uncertainty, sample diagenesis, etc.; and (4) controls or influences of key factors such ice sheets, tectonics, paleogeography, cosmic ray fluxes, biota, volcanic eruptions, and so on.
While the corresponding early Pliensbachian climatic and isotopic-shifting events are not observed in the smoothed curves of the Early Jurassic seawater temperature and carbon cycle (Dera et al., 2011), the rapid-falling event 1E CO 2 is well correlated to the nearly coeval excursion events of carbon and oxygen isotopes recorded in western Tethys and the North Atlantic (Fig. 8). The 1E CO 2 event compares well to (1) the rapid carbon isotope negative excursion of (oysters, belemnites, and brachiopods) shells from the Cleveland Basin, UK (Korte and Hesselbo, 2011), and northwest Algeria (Baghli et al., 2020); (2) that of organic matter and marine carbonates from southern Pairs Basin Peti et al., 2017) and Cardigan Bay Basin, UK (Storm et al., 2020); and (3) rapid oxygen isotope negative excursion (seawater warming) of belemnites from northern Spain (van de Schootbrugge et al., 2005). The rapid change of the stable isotope record has been called the Sinemurian-Pliensbachian boundary event (SPBE) and dated in the ammonite of the upper raricostatum and lower jamesoni zones .
Multiple hypotheses have been proposed to interpret the 5-6 • C decrease of sea surface temperatures in the late Pliensbachian (Bailey et al., 2003;van de Schootbrugge et al., 2005;Suan et al., 2010) and warming of ∼ 8 • C in the early Toarcian (Bailey et al., 2003;Suan et al., 2010), such as the sea level falling and rising, methane release, a Karoo-Ferrar eruption, the Hispanic corridor opening, etc. Perhaps these hypotheses somewhat explain the rapid change of sea surface temperatures but might not link to drastic falling of pCO 2 . As we know, atmospheric CO 2 is controlled by volcanism, weathering, vegetation on land and phytoplankton in the ocean, and orbital forcing. The Sr isotope curve shows a rapid change in the early Toarcian but does not show a change in the early Pliensachian (e.g., Jones et al., 1994), indicating a distinct transfer of weathering took place on the land only at the time of the T-OAE. No robust evidence shows the rapid changes of terrestrial vegetation and marine primary productivity for the two intervals, except for the floral change in western Tethys during the T-OAE (Slater et al., 2019). The Karoo-Ferrar eruption could be responsible for the rapid rising of pCO 2 but not for the falling. The orbital forcing might be an alternative for explaining this.
To sum up, the rapid falling events of the Early Jurassic pCO 2 values in the GSB are compatible with the response of stable isotopes (carbon cycle) and seawater temperature from coeval marine sediments in a total tendency and eventful change but are not harmonized at a high-resolution timescale. Whatever caused the rapid variations of sea surface temperatures, stable isotopes, and pCO 2 , their near concordance suggests that it is a positive feedback of the sea surface temperature and carbon cycle to the pCO 2 in trend and event through the Early Jurassic, whereas the higher frequency changes in the Sinemurian-Pliensbachian may support other causal driving of the climate, such as orbital forcing (Storm et al., 2020).
1. An overall warm and hot and (semi)arid climate dominated the GSB during the Early Jurassic, possibly accompanied by an occasionally evaporitic climate in the Toarcian. This (semi)arid climate in GSB is comparable with that in the Colorado Plateau, western North America, but distinct from the relatively warm and humid terrestrial climate recognized in other places on the Chinese mainland (e.g., Qaidam, Inner Mongolia, and the Lower Yangtze) and the high latitudes of the Southern Hemisphere.
2. The Early Jurassic pCO 2 values show that a range between 980 and 2610 ppmv is ∼ 3.5-10 times the preindustrial value 275 ppmv and that the mean 1720 ppmv is ∼ 6 times the pre-industrial value. Three phases of pCO 2 values were distinguished: 1500-2000 ppmv (mean ∼ 1700 ppmv) in the Sinemurian age, 1000-1500 ppmv (mean ∼ 1300 ppmv) in the Pliensbachian age, and 1094-2610 ppmv (mean ∼ 1980 ppmv) in the early Toarcian. Two events of rapidly falling pCO 2 were also recognized: a ∼ 1000-1300 ppmv drop at the Sinemurian-Pliensbachian boundary and a quickly falling (rising) event of ∼ 1500 ppmv in the early Toarcian. The phases and events manifest the perturbation of pCO 2 in the Early Jurassic.
3. The perturbation and rapid falling events of the Early Jurassic pCO 2 from the GSB are compatible with the carbon cycle and seawater temperature from coeval marine sediments in the North Atlantic and western Tethys in terms of total tendency and eventful change. The compatibility suggests that there is a positive linkage of the sea surface temperature and carbon cycle to the pCO 2 through the Early Jurassic. On the contrary, differences at a high-resolution timescale imply additional climate drivers, such as orbital forcing, are important in the Sinemurian-Pliensbachian record.
Data availability. All the data are available in the Supplement.
Author contributions. XL designed the concept and working idea and created the paper text, tables, and figures. JW took part in the field work and contributed to writing the paper and analyzing carbon and oxygen isotopes. TR polished the language and advised on some ideas. MZ, ZW, and CZ took part in field work and took some samples.
Review statement. This paper was edited by Appy Sluijs and reviewed by Dan Breecker and Helmut Weissert.