Sedimentology and Geochemistry of the Upper Permian Linghao Formation Marine Shale, Central Nanpanjiang Basin, SW China

The Upper Permian Linghao Formation marine shale and contemporaneous transitional shale are the most potential shale gas targets in the Nanpanjiang basin, which is characterized by considerable TOC content, wide distribution, and considerable shale thickness. On the basis of division in Linghao Formation, petrographic, mineralogical, and high-resolution geochemical analyses were integrated to reveal the sedimentary environment including paleoproductivity, paleoredox conditions, detrital influx, paleoclimate, and the paleosalinity. There are two organic-rich shale intervals in Linghao Formation, which are Ling 1 member and the lower Ling 3 member. The lower Ling 1 is dominated by deep-water shelf facies, which are characterized by high TOC value (0.93%–6.36%, avg. 2.43%), high detrital influx proxies (Zr, 746–1508 ppm, avg. 1093 ppm; Ti, 19278–128730 ppm, avg. 16091 ppm), relatively warm–humid paleoclimate condition (CIA*, 75.94-91.90, avg. 82.26), low paleosalinity proxies (Sr/Ba, 0.13-0.34, avg. 0.22), and high paleoproductivity (P/Al (10−2), 1.06-2.06, avg. 1.63; Mn/Ca (10−3), 27.37-291.69, avg. 128.07). Detrital influx including gravity flow plays a critical role in the enrichment of organic matter. The sedimentary environment of upper Ling 1 and lower Ling 3 is the same as that of lower Ling 1. Unlike lower Ling 1, these intervals are characterized by low detrital influx proxies, moderate weathering, and relatively high paleosalinity proxies. The volcanic ash of Emei volcanism and felsic volcanism in South China plays a critical role in the enrichment of organic matter in upper Ling 1 and lower Ling 3, respectively. The sedimentary models for Linghao Formation organic-rich shale can reveal factors controlling the enrichment of organic matter.


INTRODUCTION
Shale oil and gas resources have been a significant part of unconventional resources in recent years (Jia, 2017;Liang et al., 2018). Organic-rich marine shales are the main target rocks for the exploration and development of unconventional hydrocarbon resources (Clarkson et al., 2012;Zou et al., 2019;Qiu and Zou, 2020;Xi et al., 2022;Zhang et al., 2022). In China, a total of 12 important organic-rich marine shale units are recognized in Mesoproterozoic through Cenozoic strata (Zou et al., 2019).
Under the exemplary role of the United States' "shale gas revolution" (Chen et al., 2019a;Chen et al., 2019b;Chen et al., 2019c), China has achieved a breakthrough in shale gas industrialization since 2010 (Jia, 2017;Xi et al., 2017;Xi et al., 2018). Weiyuan, Changning, Jiaoshiba, and other shale gas fields have been built, and large-scale shale gas development has been realized in marine organic-rich shales of the Late Ordovician-early Silurian (Zhang et al., 2019a;Zhang et al., 2019b;Zou et al., 2019;Zhang et al., 2020a;Chen et al., 2021). The studies on this set of shale are abundant and various, which lays a good foundation for efficient exploration and development (Zhang et al., 2020b;Zhang et al., 2020c). Permian marine or transitional organic-rich shale in South China is characterized by extensive distribution, considerable TOC content, and shale thickness, which is the most promising alternative strata for future shale gas exploration in China (Liu et al., 2013;Han et al., 2017). However, this set of shale is still in the early stage of exploration, and the lack of relevant studies seriously restricts the progress of shale gas exploration. Scientific problems such as spatial distribution, microscopic pore structure, and genesis of Permian organic-rich shale need to be solved first.
Geochemical proxies have been widely used to illustrate the effect of paleoproductivity, paleoredox condition, paleosalinity, and terrigenous input on the deposition of organic-rich shales (Ding et al., 2018;He et al., 2020;Liang et al., 2020;Leifu Zhang et al., 2021). Based on the division of the Linghao Formation in central Nanpanjiang basin, macroscopic sedimentary characteristics, microscopic petrological characteristics, and high-resolution geochemical analyses are integrated into this study to reveal depositional conditions (including terrigenous input, paleoclimate condition, paleosalinity, paleoproductivity, and paleoredox) of different members in Linghao Formation. This study aims to reveal the distribution and genesis of marinehosted organic-rich shale in the Permian Linghao Formation and provide scientific guidance for shale gas exploration of Permian marine.

GEOLOGICAL BACKGROUND
The Nanpanjiang basin, with an area of 380, 000 km 2 , is located at the junction of Yunnan, Guizhou, and Guangxi, and the geotectonic position is located on the southwestern margin of South China Block (Xia et al., 2018) (Figures 1A,B). Similar to the Carboniferous Datang Formation in the Nanpanjiang basin, the Linghao Formation is a special and representative stratigraphic unit formed by the deep-water sedimentary environment between isolated platforms during the late Permian. Observations from the studied well and outcrops show that the Linghao Formation is in conformable contact with underlying limestone or muddy limestone of the Maokou Formation, and the overlying Luolou Formation. According to the combination of lithology, logging characteristics, and lithology, the Linghao Formation can be divided into Ling 1 member, Ling 2 member, and Ling 3 member. The Ling 1 member consists of shale with Emeishan basalt interbed. Previous studies suggested that the eruption timing of Emeishan basalt is between 257 Ma and 259 Ma (Shellnutt et al., 2012), and its formation period belongs to Wuchiapingian. The Ling 2 member is composed of mudstone and limestone. The Ling 3 Member is composed of shale ( Figure 1C). According to previous studies, the northwest Nanpanjiang basin is dominated by swamp, tidal flat, and lagoon facies during Linghao period (Luo et al., 2018;He et al., 2020). The contemporaneous strata of the Linghao period comprise mudstones, silty mudstones, siltstones, limestones, and coals (He et al., 2020). The studied well is located in the central Nanpanjiang basin which is dominated by shelf facies during Linghao deposition ( Figure 1B).

Sedimentological Characterization
The study was well carried out continuous coring in Linghao Formation. Macroscopic observation of core and outcrops identified different rock types, including limestone, mudstone, silty mudstone, and shale, and record the lithology combination of different members. 150 core samples were made into thin sections and the core observation results were further verified by a Leica DM4 M optical microscope. Macroscopic characteristics of sedimentary texture or structure in typical outcrops were described in the study, with emphasis on the difference between each member of the Linghao Formation.

Organic Matter Content Analysis
TOC was measured by a Leco carbon/sulfur analyzer in the State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation of China, and the analytical precisions are ±0.5%. 116 samples were treated with 10% hydrochloric acid to remove carbonates. The acid-treated sample was washed with distilled water to neutral, then the sample was dried in an oven at 60°C-80°C. Dried samples were added to the cosolvent and fully burned in the high-temperature oxygen flow, ensuring that the organic carbon can be completely converted into carbon dioxide, and the content of total organic carbon was tested by the infrared detector.

Major Element Composition Analysis
The content of major elements (Al, Si, Ti, Ca, Fe, Mg, Mn, K, and Na) was determined by X-ray fluorescence (XRF) spectroscopy. 60 powder samples were used to eliminate the mineral and particle size effects and cast to suitable fused glass beads to fit the X-ray fluorescence spectrometer. The fluorescent X-ray intensity of the elements was measured. Based on the calibration curve or equation, the interference effect between elements was corrected and the element content was obtained. The accuracy of the XRF analysis is better than 1% for all major oxides.
Paleoclimate conditions were evaluated by the chemical index of alteration (CIA), which was calculated using the following formula (Nesbitt and Young, 1982;Price and Velbel, 2003;Liu et al., 2017): (1)

Trace Element Composition Analysis
Trace elemental contents were measured by using inductively coupled plasma-mass spectrometry (ICP-MS). 60 samples were dissolved with hydrofluoric acid and nitric acid in a closed container. The hydrofluoric acid was wiped out by evaporation on the electric heating plate and then dissolved by nitric acid. After dilution, the samples were directly measured by ICP-MS. The analytical uncertainties are estimated to be 5%. Standard rock reference material (GSD-9) was used to monitor the analytical accuracy and precision. Authigenic enrichment of redox-sensitive trace elements was widely used to indicate redox changes in the water column (Algeo and Tribovillard, 2009;Algeo and Liu, 2020). The enrichment factor (EF) is calculated by comparing the Al-normalized metal content to that of the continental crust (Wedepohl, 1995). The calculation formula is as follows: ( 2 ) In the formula, X represents the content of element X, and (X/ Al) CC represents the X/Al ratio of the continental crust (Wedepohl, 1995).

Ling 1 Member
Based on observation of typical outcrop and drilling core from the studied well, the lower Ling 1 consist of black shale with interbed of gravity flow sandstone ( Figure 2A). Scour surface can be observed between gravity flow sandstone and black shale in the core ( Figure 3A). Above Emeishan basalt belongs to the upper part of Ling 1 Member black shale with multiple interbeds of bentonite ( Figures 2B, 3B). Combined with previous study results (Shellnutt et al., 2012), the bentonite interbeds are derived from Emeishan volcanic plume. Under plane-polarized light, thin sections show that black shale with abundant bioclasts including sponge spicule ( Figure 4A), radiolarian, and mollusk (Figures 4B,C) in upper Ling 1 suggest high paleoproductivity.

Ling 2 Member
The Ling 2 member consists of multiple rock types including micritic limestone, mudstone ( Figure 2C), silty mudstone ( Figure 2D), siltstone ( Figure 3C), and bioclastic limestone ( Figure 3D). Frequent interbeds between different types of sedimentary rock can be observed in core and outcrops ( Figure 3E). Under planepolarized light, thin sections show that mudstone in Ling 2 Member is rich in bioclasts including branchiopod, arthropoda, and mollusk ( Figure 4D). Silty mudstone is rich in silt-sized detrital quartz ( Figure 4E). Bioclastic limestone in Ling 2 is mainly composed of shallow benthic bioclasts ( Figure 4F).

Ling 3 Member
The lower Ling 3 Member is composed of black shale with horizontal bedding (Figures 2E,F). In the core, bentonite interbeds are developed in black shale ( Figures 3F,G). In certain intervals of black shale in lower Ling 3 ( Figure 3H), black shale is characterized by massive bedding without any

TOC Content Characteristics
Based on the TOC measurement results of marine shale or mudstone ( Figure 5), the TOC of Ling 1 marine shale varies widely, ranging from 0.14% to 6.36%, with an average of 2.36% (n = 14). The TOC range of marine shale in Ling 2 is relatively small, ranging from 0.50% to 3.08%, and the average value is only 1.30% (n = 19). For lower Ling 3, the range of TOC is 0.58%-4.36%, with an average of 2.43% (n = 9). In upper Ling 3, the range of TOC is 0.79%-3.24%, with an average value of 1.87% (n = 4).

Major Element Content
The major element content of different members in the Linghao Formation is different ( Figure 6). The Ca content of Ling 2 member is obviously higher, ranging from 2.25% to 19.81%, with an average of 7.40% (n = 31). The content of Ca in Ling 1 is 0.31%-13.90%, with average value of 3.04% (n = 14). The average content of Ca in a high TOC interval (3,756 m-3805 m) is 2.83%. The content of Ca in lower Ling 3 member is 1.17%-10.31%, with an average of 4.26% (n = 10). The Ca content of upper Ling 3 is 3.18%-12.93%, with an average of 8.04% (n = 7).

Trace Element Content
The content of indicator trace elements or their ratios has been extensively applied to reveal the paleoredox conditions, paleosalinity, paleoclimate, and paleoproductivity of the The variation trend of Sr/Cu is opposite to that of CIA* (Figure 8). The Sr/Cu of lower Ling 1 is the lowest, ranging from 0.34 to 20.54, and the average value is only 3.01 (n = 9). The average Sr/Cu ratios for upper Ling 1, Ling 2 and upper Ling 3 is 5.26 (4.02-6.98, n = 5), 8.35 (1.10-24.88, n = 31) and 6.88 (4.39-10.43, n = 6), respectively. The Sr/Cu of lower Ling 3 is relatively low, ranging from 1.36 to 11.77, with an average of only 3.60 (n = 10).
The variation of Sr/Ba is similar to that of Sr/Cu (Figure 8). The Sr/Ba of lower Ling 1 is the lowest, ranging from 0.13 to 5.48, and the average value is only 0.85 (n = 9). The average Sr/Ba ratios for upper Ling 1, Ling 2, and upper Ling 3 is 1.04 (0.59-1.69, n =

Detrital Influx Proxies
As important elements in sedimentary rocks, the contents of Al, Ti, and Zr in rocks are almost not affected by weathering or diagenesis, so these elements are applied for evaluation of the influence degree of terrestrial input (Li et al., 2017;Liu et al., 2021). Al only exists in the clay minerals of fine-grained sedimentary rocks, while Ti and Zr are mainly assigned to clay, sand, and silt particles composed of ilmenite, rutile, and augite (Caplan and Bustin, 1998;Murphy et al., 2000). The detrital influx proxies represented by Al, Ti, and Zr, are generally high for lower Ling 1 (Figure 6), especially in the high TOC interval (3,800 m-3807 m). It is suggested that the gravity flow delivered more terrestrial organic matters into the water during the depositional stage of the lower Ling 1. This terrigenous organic matter, on the one hand, provides a lot of nutrients for microbes or directly increases the organic carbon content of shales (Khripounoff et al., 2009). Zr/Al and Ti/Al ratios are thought to closely relate to the coarser part of the sediments (Bertrand et al., 1996;Caplan and Bustin, 1998). Good Ti-Al correlation suggests that Ti comes from the lattice of clay minerals or stable terrigenous clastic materials (Leifu Zhang et al., 2021). Ling 2 and Ling 3 are characterized by good Ti-Al correlation, while Ling 1 is not ( Figure 9A). It is suggested that the detrital influx of Ling 2 and Ling 3 is relatively stable, while that of Ling 1 is not. Zr usually exists in clay minerals or heavy minerals of silt size (e.g., zircons) (Rachold and Brumsack, 2001;Liu et al., 2017). The correlation between Al and Zr is relatively good for Ling 2 and Ling 3, while Ling 1 is not ( Figure 9B). The results also show that the detrital influx of Ling 1 is strongly unstable and heterogeneous.

Paleoclimate Conditions
The warm-humid climate is favorable for the atmospheric water cycle, accelerating the chemical weathering intensity, and causing the transportation of nutrients to marine or lakes. Moreover, a warm-humid climate benefits the prosperity of plankton in surface water, and the burial amount of organic matter.
Since the paleoclimate variations will alter the chemical composition of rocks to a certain degree, the paleoclimate of the depositional period can be reflected by the change in the chemical composition of rocks in turn. Sr and Cu are very sensitive to climate change, so the concentrations of these two elements can be applied to illustrate the paleoclimate conditions. Sr is likely to be lost due to weathering or leaching, while Cu is relatively stable (Leifu Zhang et al., 2021). Under the warm conditions, due to strong chemical weathering, Sr is more easily lost, resulting in increased Sr/Cu in sediments. Under dry and hot conditions, the weathering intensity is relatively low, and more Sr elements remain in the parent rock, resulting in lower Sr/Cu values in sediments. For dry and hot climates, Sr/Cu ratio is between 1-5, while for warm and humid climates, Sr/Cu ratio is greater than 5 (Yandoka et al., 2015;Xie et al., 2018). The average Sr/Cu ratio of Ling 1 is only 3.81, especially the Sr/Cu ratio for lower Ling 1 interval is the lowest (Figure 8). The average Sr/Cu ratio of Ling 3 is 4.83. The average Sr/Cu ratio of Ling 2 is 8.35, which is significantly higher than that of other members. It is suggested that a relatively dry climate during the depositional stage of Ling 1 and Ling 3 leads to less terrigenous input, which benefits the deposition of organic-rich shale. The relatively high Sr/Cu ratio for Ling 2 indicates that the paleoclimate tends to be warm and humid.
CIA* (Chemical Index of Alteration) has been extensively utilized to reflect the chemical weathering intensity (Nesbitt and Young, 1982;Bai et al., 2015). In general, high CIA* values suggested warm, humid paleoclimate and strong chemical weathering. Low CIA* values reflect dry and cold paleoclimate, and chemical weathering is weak. Bai et al. (2015) proposed that when CIA* is between 50 and 65, it reflected the cold-dry climate under the background of low chemical weathering. When CIA* is between 65 and 85, it represents the warm-humid climate under the background of moderate chemical weathering. When CIA* is between 85 and 100, it represents the hot and humid paleoclimate under strong chemical weathering. The CIA* for lower Ling 1 is the highest, which is in the middle range (65-85) of chemical weathering. The CIA* values of the Ling 2 and Ling 3 are similar (69-79) (Figure 8), which indicates that the climate tends to be warm and humid. It seems that the climate of lower Ling 1 tends to be the warmest and humid interval during the deposition of the Linghao Formation. According to the Sr/Cu ratio and the field outcrop observation results, it is most likely that the high detrital influx represented by gravity flow leads to the high CIA* value in lower Ling 1.

Paleosalinity
Paleosalinity is an important index for restoring the paleosedimentary environment. Ba and Sr were proved to be two indicator elements sensitive to paleosalinity (Wang et al., 1979;Wang, 1996). Sr and Ba showed different enrichment under different salinity. Under low salinity, Ba and Sr existed in soluble bicarbonate. With the increase in salinity, Ba gradually precipitated in the form of BaSO 4 , and Ba content in water decreased compared to Sr (Wolgemuth and Broecker, 1970). Under higher salinity, Sr is precipitated only in the form of SrSO 4 (Wang et al., 1979). The enrichment of Sr in sediments is associated with an increase in salinity, so high Sr/Ba ratios represent an increase in seawater depth and salinity (Chegrouche et al., 2009;Wei and Algeo, 2020). Therefore, Sr/ Ba ratio has been extensively applied to qualitatively restore paleosalinity Wei and Algeo, 2020). Sr content and Sr/Ba ratio in sedimentary rocks have a strong positive linear relationship with paleosalinity, while Ba content is negatively correlated with paleosalinity. Since Sr and Ca are similar in atomic radius, Sr content may be much higher in sediments containing carbonate components. Before restoring paleosalinity by Sr/Ba ratio, the interference of carbonate rocks on Sr concentration must be excluded. CaO content was used as a proxy for carbonate content and the positive Sr/Ba-CaO relationships are present in Ling 1 and Ling2 samples with CaO content larger than 10% ( Figure 10A). We selected CaO = 10% as the maximum carbonate threshold for Linghao Formation sediments in this study. In Figure 10B, the CaO threshold (CaO = 10%) was applied to eliminate samples that probably contained carbonate-hosted Sr Wei and Algeo (2020) proposed that sedimentary Sr/Ba ratios of <0.2, 0.2-0.5, and >0.5 are indicative of freshwater, brackish, and marine water, respectively. Sr/Ba ratios are different in different members (Figure 8). Ling 2 and Ling 3 sediments are deposited in marine environments, and the average Sr/Ba ratio is 1.86 and 0.96, respectively. The Sr/Ba ratio of Ling 1 is higher than that of the other two members, with an average value of 0.40, which represents brackish water-fresh water ( Figure 10B). Combined with outcrop observation (Figure 2A), it is suggested that gravity flow may be the direct cause of paleo-water salinity closer to freshwater during the deposition of Ling 1.

Paleoredox Conditions
Trace element proxies including U/Th, U EF , and Mo EF were extensively applied to reveal the redox conditions of paleowater, and the smaller these proxies represent the higher theoxidation degree, and the larger the ratio represents the stronger the reduction degree (Jian Liu et al., 2018;Zhang et al., 2019a;He et al., 2020). Previous studies point out that high U/Th ratios were consistent with highly reducing conditions, with U/Th > 1.25 representing dysoxic to anoxic conditions and U/Th < 0.75 representing normal oxic conditions (Hatch and Leventhal, 1992;Jones and Manning, 1994;Tribovillard et al., 2006). Authigenic Mo, authigenic U enrichment, and Mo-U covariant models have been applied to illustrate redox conditions and water mass limitation Algeo and Tribovillard, 2009;Algeo and Liu, 2020). In general, the oxic condition showed little or no enrichment of authigenic U and Mo, while the anoxic conditions showed strong enrichment of authigenic U and Mo (Algeo and Tribovillard, 2009). The U/Th values of each member are less than 0.75, suggesting that the environment condition of the Linghao Formation is generally inclined to oxic environments ( Figure 8). The data points of each member are all located between the oxic end and the dysoxic end of unrestricted marine trend ( Figure 11). It is suggested that the organic-rich shale of Ling 1 and Ling 3 were deposited in an oxic/dysoxic environment. It seems that the redox conditions of Ling 1 and Ling 3 are not favorable for the preservation of organic matter.

Paleoproductivity Proxies
The primary productivity in marine settings is thought to be a significant factor in the formation of organic-rich shales Zhang et al., 2019a). As an essential nutrient for plankton, P is widely used to evaluate the primary productivity of the paleomarine environment Jian Cao et al., 2018;He et al., 2020). To avoid the influence of the abiotic P from terrigenous input, the P/Al ratio can be used as a proxy for paleoproductivity in a marine environment, using the ratio of Mn/Ca as a reference value (Jian . The P/Al ratio for the lower Ling 1 ranges from 1.06 to 2.06 (average 1.63), and that for upper Ling 1 varies from 1.06 to 2.39 (average 1.80). The P/Al ratio for Ling 2 ranges from 1.23 to 3.17 (average 2.07), and that for lower Ling 3 varies from 1.30 to 2.46 (average 2.14). The P/Al ratio for the upper Ling 3 ranges from 0.78 to 2.44 (average 1.65). It seems that the low P/Al ratio of the Linghao Formation reflects low productivity and the vertical variation trend between TOC and P/Al ratio is not very consistent (Figure 8). In contrast, the vertical variation trend between TOC and Mn/Ca ratio is consistent. This phenomenon can be attributed to the effect of redox conditions on enrichment P in sediments (Algeo and Ingall, 2007;Westermann et al., 2013). In reducing environments, P dissolves from sediments into the water, while in oxidizing environments, P is easily adsorbed in Fe compounds (Rimmer et al., 2004), and this could have led to the high TOC and the low ratio of P/Al in the Linghao Formation. Therefore, the low P/Al value in the high TOC intervals of the Linghao Formation does not represent the low paleoproductivity, but on the contrary, represents the high productivity.

SEDIMENTARY MODEL OF ORGANIC-RICH SHALE
Based on the abovementioned discussion, a sedimentary model of organic-rich shale in the Linghao Formation is proposed. The sedimentary environment during deposition of the Maokou Formation is dominated by shallow-water shelves. During lower Ling 1 deposition, the sedimentary environment is dominated by deep-water shelves due to coeval basementrooted faulting under extension stress background (Changqun Wang et al., 2022). Strong weathering increases the detrital influx, and active basement-rooted faulting induces frequent gravity flow ( Figure 12A). Meanwhile, on the one hand, gravity flow brings a large amount of terrestrial organic matter, on the other hand, it provides rich nutrients for plankton and FIGURE 11 | Authigenic Mo-U covariation patterns for Linghao Formation marine shale. The diagonal dashed lines represent the Mo/U molar ratios of the seawater (sw). Grey area represents the "unrestricted marine" trend, characteristic of depositional systems with no limited trace metal renewal. The yellow area represents the "particulate shuttle" trend, characteristic of environments in which intense redox cycling of metals occurs within the water column (Algeo and Tribovillard, 2009;Tribovillard et al., 2012) Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 914426 improves paleoproductivity ( Figure 12A). Emeishan volcanic activity formed a set of basalt on lower Ling 1 but did not change the sedimentary environment. The bentonite layers of volcanic ash origin in the organic-rich shale of upper Ling 1 indicate that the volcanic ash benefits the enrichment of organic matter. Firstly, volcanic ash provides biological nutrients for plankton. Secondly, the volcanic ash layer has a strong oxygen insulation ability ( Figure 12B). During deposition of Ling 2, the sedimentary environment is dominated by shallow-water shelves when basement-rooted faulting entered the dormant period, and this kind of environment is not favorable for organic matter enrichment ( Figure 12C). With the reactivation of basementrooted faulting, the sedimentary environment changed into deepwater shelves again during the deposition of lower Ling 3. In Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 914426 13 general, the sedimentary model of organic-rich shale in lower Ling 3 is very similar to that in upper Ling 1 ( Figure 12D), but the difference is that the bentonite in lower Ling 3 originates from felsic volcanism in South China Hua Zhang et al., 2021).

CONCLUSIONS
1) There are two organic-rich shale intervals in the Linghao Formation, which are Ling 1 Member and lower Ling 3 Member. Organic-rich shales of these intervals are deposited in the environment of deep-water shelf triggered by basement-rooted faulting under extension stress background.
2) The changes of TOC and its correlation with geochemical proxies suggest that paleoproductivity is the critical factor that controls the enrichment of organic matter during the deposition of Linghao Formation organicrich shales. For lower Ling 1, detrital influx including gravity flow benefits paleoproductivity, thus controlling the enrichment of organic matter. For upper Ling 1 and lower Ling 3, the volcanic ash of Emei volcanism from Southwest China and felsic volcanism from South China plays a critical role in paleoproductivity, respectively, thus controlling the enrichment of organic matter.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.