Earth’s Orbital Control on Organic Matter Enrichment in the Black Shales of the Wufeng–Longmaxi Formation in the Upper Yangtze Region, South China

China’s most successful horizon for shale-gas exploration and development is the Upper Ordovician Wufeng Formation through the Lower Silurian Longmaxi Formation in its Upper Yangtze Region. In this study, the Wufeng–Longmaxi Formation black shales in the Upper Yangtze Region are analysed to determine their microstructural characteristics, total organic carbon (TOC) content, and well log characteristics and to identify information on the parameters of Earth’s orbit from its gamma-ray (GR) data series. Then, paleoenvironmental modes are established over a long time scale. On this basis, the Earth’s orbital control on organic matter enrichment in the black shales is examined. The black shales of the Wufeng Formation and the lower Longmaxi Formation are rich in biogenic siliceous fossils and framboidal pyrite. Their TOC content, GR value, and uranium (U) content increase gradually from the bottom of the Wufeng Formation to the Guanyinqiao Member at its top, peak at the Guanyinqiao Member, and then decrease gradually up to the Longmaxi Formation. Approximately six and seven long eccentricity cycles can be identified from the GR curves of the Wufeng Formation and the lower Longmaxi Formation, respectively. During the long eccentricity minima, corresponding to Earth’s cold period, the black shales have a relatively high level of enrichment of organic matter. This study can provide an important reference for investigating the mechanism by which Earth’s orbits control the climate and sedimentary environment, as well as the mechanism of organic matter enrichment.


INTRODUCTION
Periodic changes in the parameters (i.e., eccentricity, obliquity, and precession) of Earth's orbit lead to periodic fluctuations in the amount of solar radiation and thus the surface climatic system on Earth. These periodic climatic changes are recorded in sediments sensitive to climate change, resulting in cyclic features in sedimentary formations, known as Milankovitch cycles (Schwarzacher, 1989;Schwarzacher, 1993;Strasser et al., 2006). Researchers have identified Milankovitch cycles (or Earth's orbital parameters) in sedimentary formations using metrics reflective of palaeoclimatic changes, such as carbon and oxygen isotopes, gamma-ray (GR) values, magnetic susceptibility, organic carbon, and element ratios in conjunction with methods such as wavelet and spectral analysis (Weedon et al., 2004;Sageman et al., 2006;Huang et al., 2010;Cheng et al., 2017;Jin S. et al., 2020;Ravier et al., 2020;Shi et al., 2020;Shi et al., 2021;Kochhanna et al., 2021;Zhang et al., 2021). Analysis of Earth's orbital parameters in sedimentary formations has produced a series of progress, e.g., the improvement of the International Geological Time Scale (Lourens et al., 2004), the estimation of the duration and rate of deposition of formations (Huang et al., 2010;Cheng et al., 2017), the estimation of the duration of major geological events and the determination of their trigger mechanisms (Sageman et al., 2006;Zhong et al., 2020), and the recognition of the effects of orbits on sedimentary environments (Weedon et al., 2004;Yao et al., 2015;Jin S. et al., 2020;Ravier et al., 2020;Shi et al., 2020;Shi et al., 2021;Kochhanna et al., 2021;Zhang et al., 2021). Rachold and Brumsack (2001) summarized some effects of Earth's orbits on climatic processes and sedimentary environments (including sea-level changes, weathering processes, organic productivity, and reduction-oxidation conditions). However, the mechanisms by which orbits control sedimentary environments have yet to be understood. For example, research on the mechanism by which orbits affected the sedimentary environment and organic matter enrichment in Palaeozoic shales stagnates in only two hypothetical models. 1) The eccentricity cycle controls primary productivity. This model states that during maximum eccentricity, intense weathering causes oceanward transport of terrestrial detritus and nutrients in large quantities and thereby aids in increasing oceanic primary productivity. That is, a high level of enrichment of organic matter occurred in shales during maximum eccentricity. 2) The eccentricity cycle controls organic matter preservation conditions. This model states that a change in seasons occurs less frequently during minimum eccentricity, thereby providing a stable sedimentary environment with relatively good preservation conditions for organic matter that allow for a high level of enrichment of organic matter in shales (Gambacorta et al., 2018;Jin S. et al., 2020).
Extraction of shale gas from the black shales of the Upper Ordovician Wufeng Formation and the Lower Silurian Longmaxi Formation in southern China's Yangtze Region has been a great success (Zou et al., 2018). Through relatively extensive investigations of the paleoenvironment during the deposition of the Wufeng-Longmaxi Formation black shales, previous researchers have noted that the Yangtze Sea had relatively high primary productivity and an anoxic bottom-water environment during this period Zhao et al., 2019;Jin C. et al., 2020;Yang S. et al., 2021;Zan et al., 2021). In addition, they have discussed the controlling factors of organic matter enrichment in the black shales of the Wufeng-Longmaxi Formation from multiple perspectives, including upwelling currents, volcanic ash, sea-level changes, and occlusive settings Zhao et al., 2019;Jin C. et al., 2020;Yang S. et al., 2021;Zan et al., 2021). Million-year-scale palaeo-ocean and palaeoclimate research findings show that while the shales in the Wufeng-Longmaxi Formation in the Yangtze Region are several hundred metres in thickness, organic matter-rich black shales are concentrated mainly within the bottom several tens of metres (Zou et al., 2018). Jin S. et al. (2020) and Zhang et al. (2021) identified the cycles of Earth's orbital parameters (e.g., eccentricity, obliquity, and precession) from the black shales of the Wufeng-Longmaxi Formation in the Upper Yangtze Region. The investigation of the Earth's orbital control on organic matterrich enrichment in the black shales of the Wufeng-Longmaxi Formation is still in the exploratory stage.
To address the above problem, this study analyses black shales of the Wufeng-Longmaxi Formation in well N211 and the Shuanghe profile to determine their microscopic characteristics, total organic carbon (TOC) content, and well log characteristics. Information for the parameters of Earth's orbit contained in the black shales of the Wufeng-Longmaxi Formation is identified using the GR curves of well N211. Then, paleoenvironmental modes are established, and the controlling effects of a long eccentricity cycle on organic matter enrichment in the black shales are examined. This study can become an important reference for investigating the mechanisms by which Earth's orbits control climate and sedimentary environments and the mechanism of organic matter enrichment in black shales.

GEOLOGICAL SETTING
The Yangtze Block collided and merged with the Cathaysian Block during the Middle Ordovician to form the South China Block (Wang, 1985;Chen et al., 2004;Su et al., 2009). During the Late Ordovician through the early Silurian, the South China Block was situated in the equatorial tropical or subtropical region on the western margin of the Gondwana continent ( Figures 1A,B) (Metcalfe, 1994;Li et al., 2005;Huang et al., 2018). In this period, the Central Sichuan Submarine Uplift and the Central Guizhou Uplift were formed on the western margin and in the southern part of the Upper Yangtze Region, respectively, while the Xuefengshan Uplift and Western Hubei Submarine Uplift developed in the eastern part ( Figure 1B). These factors led to the formation of an expansive, low-energy, and undercompensated sedimentary environment (Chen et al., 2004;Yan et al., 2009;Chen et al., 2013), where organic matterrich black shales were deposited over a large area. Of these black shales, those with a TOC > 2% have a thickness of up to approximately 100 m ( Figure 1C) .
With a lithological composition consisting primarily of siliceous shales (Figure 2A), the Wufeng Formation and the lower Longmaxi Formation in the Upper Yangtze Region are rich in graptolite fossils ( Figure 2B) and are intercalated with thin pyrite layers ( Figure 2C) and thin potash bentonite layers ( Figure 2D). The Guanyinqiao Member at the top of the Wufeng Formation consists of three facies-siliceous shale facies, calcareous-siliceous shale facies, and limestone facies-in the Upper Yangtze Region (Wang et al., 2016). The siliceous and calcareous-siliceous shales contain an abundance of Hirnantian faunal fossils ( Figure 2E), suggesting the presence of a cold, regressive, relatively shallow-water, and oxygen-rich sedimentary environment (Rong et al., 2010). Chen et al. (2015) divided the Wufeng-Longmaxi Formation into 13 graptolite biozones, each with a corresponding isotopic age ( Figure 3A). The organic matter-rich black shales are located in the Wufeng Formation and the lower Longmaxi Formation, primarily corresponding to the WF1-LM5 graptolite zones (Zou et al., 2018). Therefore, only the black shales in the Block during the Late Ordovician through the Early Silurian. The South China Block was situated in the equatorial tropical or subtropical region on the western margin of the Gondwana continent at approximately 90°counterclockwise from its present-day location, according to Cocks and Torsvik (2002) and Huang et al. (2020). (B) Paleogeographic map of the Yangtze Block during the Late Ordovician through the Early Silurian. The present-day northern margin of the Yangtze Region was its western margin during the Late Ordovician through the Early Silurian. The Upper Yangtze Region is surrounded by the Central Sichuan Submarine Uplift, Central Guizhou Uplift, and Western Hubei Submarine Uplift, according to Zhou et al. (2015); Wang et al. (2017); Liu et al. (2017), and Tang et al. (2020). (C) Distribution of the organic matter-rich shales in the Wufeng-Longmaxi Formation in the Upper Yangtze Region. Those with a TOC > 2% have a thickness up to approximately 100 m, and the thickest organic-rich shales are near Shizhu and well Ys2, according to Wang et al. (2017

Samples and Data
Five samples of the black shales of the Wufeng Formation and the lower Longmaxi Formation were collected from the Shuanghe profile and then analysed to determine their microstructural characteristics. Forty black shale samples were retrieved from the Wufeng Formation and the lower Longmaxi Formation in well N211 (9 from the Wufeng Formation (including 1 from the Guanyinqiao Member) and 31 from the lower Longmaxi Formation) for TOC analysis. Well log data [GR values and uranium (U) content, data point interval: 0.1524 m] were obtained by Schlumberger Limited from well N211. The GR values and U content were measured using natural gamma-ray spectral logging.

Microscopic Characteristics of the Shales and the Total Organic Carbon Analytical Method
The microscopic characteristics of the black shales were analysed using a Leica DM4500P polarized-light microscope and a Quanta 250 FEG scanning electron microscope at the State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Chengdu University of Technology. The black shale samples were ultrasonically washed multiple times in deionized water and then dried. Then, they were ground in an agate mortar and passed through an 80-mesh screen for TOC analysis, which was performed using a CS230 C-S analyser at the Analytical Experimental Centre of the Research Institute of Petroleum Exploration and Development, PetroChina Southwest Oil and Gas Field Company.

Spectral Analytical Method
The black shales in the Wufeng Formation and the lower Longmaxi Formation were deposited at fairly different rates. Therefore, the GR values retrieved from the Wufeng Formation and the lower Longmaxi Formation were separately subjected to spectral analysis using the Acycle 2.3.1 software package based on the MATLAB platform . The main analytical methods involved in this study include multitaper-method spectral analysis and fast-Fourier transform sliding-window spectral analysis.

Microscopic Characteristics of the Shales
The black shales of the Wufeng Formation and the lower Longmaxi Formation are rich in biogenic siliceous fossils, including siliceous Radiolaria with diameters of generally less than 200 μm ( Figure 4A) and sponge spicules with diameters of approximately 50 µm and lengths of approximately 300 μm ( Figure 4B). Ma et al. (2018) suggested that sources of the siliceous matter in the Wufeng and lower Longmaxi Formation black shales include siliceous organisms (67-90%) and a small amount of terrestrial detritus. This finding, combined with the microscopic analysis conducted in this study, clearly shows that the siliceous components of the Wufeng and lower Longmaxi Formation black shales are primarily biogenic. Extensive framboidal pyrite with a particle diameter of less than 5 µm has developed in the Wufeng and lower Longmaxi Formation black shales ( Figures 4C,D), suggesting that the bottom seawater or the porewater in the sediments had an anoxic environment during the deposition of these black shales (Bond and Wignall, 2010).

Total Organic Carbon Content
The black shales of the Wufeng and lower Longmaxi Formations have a high TOC content ( Figure 5). The TOC content in the black shales of the Wufeng Formation ranges from 1.01% to 8.35% and increases gradually from the bottom up. Specifically, this TOC content is 1.01-5.75% (average: 3.66%) at the bottom-middle and 8.35% (the highest TOC content in the Wufeng and lower Longmaxi Formation black shales) in the Guanyinqiao Member at the top. In addition, the Guanyinqiao Member is rich in benthic fossils (i.e., Hirnantian fauna fossils), suggesting the presence of oxygen-rich bottom seawater during the deposition of the Guanyinqiao Member but the presence of anoxic porewater in the sediments in the Guanyinqiao Member. The TOC content in the lower Longmaxi Formation black shales ranges from 1.28% to 6.04% (average: 3.15%) and decreases gradually from the bottom up ( Figure 5). The TOC content increases gradually from the bottom of the Wufeng Formation to its Guanyinqiao Member, peaks in the Guanyinqiao Member, and then decreases gradually going further up to the Longmaxi Formation ( Figure 5).  Figure 5). In summary, the GR value and U content increase gradually from the bottom of the Wufeng Formation to its Guanyinqiao Member, peak in the Guanyinqiao Member, and then decrease gradually going further up to the Longmaxi Formation ( Figure 5).

Spectral Analysis
The TOC content and GR values in shales are strongly positively correlated ( Figure 5). Organic matter-rich shales generally have relatively high GR values. Therefore, GR data series obtained from dense sampling sites are used in this study to perform spectral analysis and examine the orbital control on organic matter enrichment. GR well logs may contain interference signals unrelated to geological factors. To better identify orbital cycles from the GR well logs, we removed the effects of the deformation of the lowfrequency portions of the well log data for the Wufeng Formation and the lower Longmaxi Formation using the "Detrending" module in the Acycle 2.3.1 software package before spectral analysis ( Figures 6A,B).
We subjected the detrended GR well log data for the Wufeng and lower Longmaxi Formations to multitaper-method spectral analysis. Analysis of the multitaper-method spectrum for the Wufeng Formation reveals that the sedimentary cycle thickness corresponding to the peak power is 1.60 m (90% confidence interval, 1.03-2.29 m) ( Figure 7A). Analysis of the multitaper- method spectrum for the lower Longmaxi Formation reveals that the sedimentary cycle thicknesses corresponding to the peak powers are 5.64, 4.83, 3.26, 1.47, and 1.16 m, respectively, each with a confidence level of greater than 90% ( Figure 7B).
To identify how the principal frequency varies with depth (i.e., the principal cycle) in the Wufeng and lower Longmaxi Formations and to ascertain the principal cycle stability with respect to depth, we subjected the GR curves for the Wufeng and lower Longmaxi Formations in well N211 to fast-Fourier transform sliding-window spectral analysis at a window step size of 0.1524 m and sliding-window lengths of 2 and 8 m, respectively. Analysis of the fast-Fourier transform sliding-window spectra reveals that the frequencies corresponding to high powers are consistent with those corresponding to the peak powers in the multitaper-method spectra FIGURE 7 | Multitaper-method spectra and fast-Fourier transform sliding-window spectra produced from the GR data series for the Wufeng Formation (A and C, respectively) and the lower Longmaxi Formation (B and D, respectively) in well N211 of the Upper Yangtze Region.
Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 938323 9 ( Figures 7C,D), suggesting that principal frequency signals are preserved in each depth section of the Wufeng and lower Longmaxi Formations. Early Silurian (440 Ma) orbital parameters-obliquity: 33.4 ± 3.8 kyr; precession: 20.9 ± 1.5, 19.9 ± 1.3, 17.1 ± 1.0, and 17.2 ± 1.0 kyr. Hence, 405:125:95: 33.4:20.4:17.15 is treated as the ratio of Earth's orbital parameter cycle in the Late Ordovician-Early Silurian. The ratio of the sedimentary cycle thickness corresponding to the peak power (at a high level of confidence) in spectral analysis is compared with this theoretical ratio of the orbital parameter cycle to ascertain the presence of orbital records in the GR data series for the Wufeng and lower Longmaxi Formations in well N211.

Cyclostratigraphy of the Wufeng-Longmaxi Formation Black Shales
The sedimentary cycle thickness corresponding to a cycle of orbital parameters can be identified only if it is greater than four times the data point interval (Weedon, 2003;Wu et al., 2011). The deposition rate of the Wufeng Formation is calculated based on its isotopic age to be 2.56-3.04 m/Myr. Multiplying this deposition rate by the duration of each of the abovementioned orbital parameters yields the possible corresponding formation thickness. A comparison of the results with four times the data point interval (0.1542 m) reveals that only long eccentricity cycles can be identified from the GR data for the Wufeng Formation. Similarly, long eccentricity and short eccentricity cycles can be identified from the GR data of the lower Longmaxi Formation in well N211, whereas obliquity and precession cycles cannot.
Assuming that 1.03-to 2.29-m-thick sedimentary cycles in the Wufeng Formation correspond to the 405-kyr long eccentricity cycle, the deposition rate of the shales is estimated to be 2.6-5.7 m/Myr, which is similar to that (2.56-3.04 m/Myr) calculated from the isotopic age of the formation. We can conclude that the 1.03-to 2.29-m-thick sedimentary cycles correspond to the long eccentricity cycle during the deposition of the Wufeng Formation in well N211 (Figures 7A,C) Figures 7B,D). Based on the orbital cycles corresponding to these sedimentary cycles, the shale deposition rate is estimated at 11.76-14.00 m/Myr, which is similar to that (9.56-11.01 m/Myr) calculated from the isotopic age of the formation, which validates the above assumption.
Based on the spectral analysis, the sedimentary cycles (with a thickness of 4.83 m) corresponding to the 405-kyr long eccentricity cycle recorded in the lower Longmaxi Formation were extracted using a Gaussian bandpass filter with a frequency of 0.21 ± 0.04 cycles/m. The filtering results show that approximately seven long eccentricity cycles are recorded in the lower Longmaxi Formation in well N211 ( Figure 6A). Similarly, the sedimentary cycles (with a thickness of 1.60 m) corresponding to the 405-kyr long eccentricity cycles recorded in the Wufeng Formation were extracted using a Gaussian bandpass filter with a frequency of 0.63 ± 0.12 cycles/m. The filtering results show that approximately six long eccentricity cycles are recorded in the Wufeng Formation in well N211 ( Figure 6B). Following these filters, we can construct a floating astronomical time scale (ATS) using the interpreted 405-kyr long eccentricity cycles of the GR series.  6A,B). This value is in the range of the duration (447.62 ± 1.1 to 440.77 ± 1.2) estimates provided by Chen et al. (2015) and Gradstein et al. (2012). The result we obtained from this study provides an important refinement of the Late Katian-Rhuddanian stratigraphy.

Paleoenvironmental Modes During the Deposition of the Wufeng-Longmaxi Formation Black Shales
In the Late Ordovician through the Early Silurian, the Yangtze Block was located near the equator on the eastern margin of the Proto-Tethys Ocean, was situated within the southeasterly or northeasterly trade-wind zone on the west side of Gondwana, and was a semioccluded sea with a westward opening (i.e., the present-day northern margin of the Upper Yangtze Region was its western margin in the Late Ordovician through the Early Silurian) (Figures 1A-C) (Metcalfe, 1994;Li et al., 2005;Huang et al., 2018). The palaeogeographic location of the Yangtze Block is comparable to the location of the present-day coastal regions near the equator on the west side of the continent (i.e., the east side of the ocean). The oceanic and climatic conditions and paleoenvironmental mode during the deposition of the black shales in the Wufeng-Longmaxi Formation may have been similar to those in the present-day coastal regions near the equator on the west side of the continent. These coastal regions are located within the southeasterly or northeasterly trade-wind zone where offshore winds are formed. Offshore winds cause a loss of surface seawater by generating westward surface offshore currents in these coastal regions. Intermediate-Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 938323 depth and deep water in the ocean rise to compensate for the lost water in the coastal regions, forming upwelling currents (Schweigger, 1958;Huyer et al., 1987;Blanco et al., 2002;Capet et al., 2004;Chavez and Messié, 2009). The Wufeng Formation and the lower Longmaxi Formation deposited in shallow water in the Upper Yangtze Region contain extensive records of storm deposits, including rip-up clasts deposited by proximal storms, hummocky cross-beds deposited by distal storms, and storm-scoured surfaces (Cheng and Wang, 1991;Li, 1997;Zhou et al., 2014), suggesting that storms were prevalent during the deposition of the black shales in the Wufeng-Longmaxi Formation in the Upper Yangtze Region. Multiple layers of potash bentonite are present in the black shales of the Wufeng-Longmaxi Formation in the Yangtze Block ( Figure 2C). They are a product of the deposition and alteration of volcanic ash (Hu et al., 2009;Su et al., 2009;Du et al., 2021;Shen et al., 2021) generated by volcanic activity on the southeastern margin of the Yangtze Block (Su et al., 2009). The thickness of the potash bentonite in the black shales shows a gradual decrease from the southern margin to the northern margin of the Yangtze Block (Su et al., 2009). The main driver of the drift of the volcanic ash through the air is winds. Therefore, the origin of the volcanic ash and the distribution of the potash bentonite in the Yangtze Block suggest the presence of offshore winds travelling from the southern margin of the Yangtze Block to the northern margin during the deposition of the black shales in the Wufeng-Longmaxi Formation.
Siliceous organisms are an important indicator of upwelling zones (Leinen et al., 1986). The profusion of biogenic siliceous fossils in the black shales of the Wufeng-Longmaxi Formation suggests the presence of upwelling currents in the Upper Yangtze Region FIGURE 8 | (A) Vertical variations in the GR value, TOC content, and productivity metric (i.e., the P 2 O 5 /TiO 2 ratio) of the Wufeng-Longmaxi Formation black shales in well N211 and the upwelling intensity [data for the P 2 O 5 /TiO 2 ratio were obtained from Wang et al. (2017), and the upwelling intensity is from Zhao et al. (2019)]. (B) Glacier development and global sea-level change during the Late Ordovician-early Silurian, according to Haq and Schutter (2008). (C) Paleoenvironmental mode during the deposition of the Wufeng-Longmaxi Formation black shales in the Upper Yangtze Region.
Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 938323 ( Figures 4A,B). Based on palaeontological, mineralogical, and geochemical evidence, previous researchers have noted that upwelling occurred during the deposition of the black shales of the Wufeng-Longmaxi Formation in the Upper Yangtze Region (Zhao et al., 2019;Yang S. et al., 2021;Zan et al., 2021). The upwelling currents brought nutrients to the upper Yangtze Sea and therefore considerably improved its sea-surface primary productivity. The flourishing, death, and decay of organisms consumed oxygen and thereby promoted the formation of an anoxic environment in the bottom seawater or the porewater in the sediments. High productivity and the anoxic environment together contributed to organic matter enrichment in the black shales of the Wufeng-Longmaxi Formation. The upwelling grew gradually more intensely from the early and middle stages of deposition of the Wufeng Formation to the depositional period of the Guanyinqiao Member and peaked during the deposition of the Guanyinqiao Member, followed by a gradual decrease during the deposition of the black shales in the Longmaxi Formation ( Figure 8A) (Zhao et al., 2019;Yang S. et al., 2021). The vertical variation in the TOC content, GR value, U content, and productivity metric (i.e., the P 2 O 5 /TiO 2 ratio) reveals that these attributes increase gradually from the bottom of the Wufeng Formation to its top, peak at the Guanyinqiao Member, and decrease gradually going further up to the Longmaxi Formation ( Figures 5, 8A), which shows a similar pattern with upwelling intensity characterized by an initial increase followed by a decrease ( Figure 8A). upper Yangtze Sea via the western margin of the Upper Yangtze Region (i.e., its present-day northern margin), forming upwelling currents. This mechanism raised the surface primary productivity of the sea and promoted organic matter enrichment in the sediments ( Figure 8C). During glacial periods, trade winds over the upwelling zone on the east side of the ocean intensified, thereby increasing the upwelling intensity (Abrantes, 1991;Weber et al., 1995;Wells and Okada, 1996;Kim et al., 2003;Romero et al., 2008;Matsuzaki et al., 2011). For example, the prevalent La Niña conditions and intensified southeasterly trade winds near the equator during the last glacial period resulted in stronger upwelling currents on the eastern margin of the ocean near the equator (Patrick and Thunell, 1997;Pisias and Mix, 1997;Mix et al., 1999;Feldberg and Mix, 2002;Kim et al., 2003;Martinez et al., 2003;Dubois et al., 2009), which might be related to the increase in the pole-equator temperature gradient or the intensification of the Walker circulation during the glacial period (Andreasen and Ravelo, 1997). The deposition of the black shales in the Wufeng-Longmaxi Formation corresponded to the glacial development-to-shrink process during the Late Ordovician to the early Silurian ( Figure 8B) (Haq and Schutter, 2008). The Guanyinqiao Member in the Upper Yangtze Region was deposited in concurrence with the lowest sea level during the deposition of the black shales in the Wufeng-Longmaxi Formation, the most intensive Antarctic glacial development (Brenchley et al., 2003;Haq and Schutter, 2008;Rong et al., 2010). Therefore, during the deposition of the Guanyinqiao Member, the Upper Yangtze Region saw its strongest offshore winds and upwelling currents and peak marine productivity. The TOC content in the Guanyinqiao Member also shows a peak value ( Figure 8C)

Orbital Control on Organic Matter Enrichment in the Wufeng-Longmaxi Formation Black Shales
In this study, the effects of a long eccentricity cycle on organic matter enrichment in black shales are analysed based on the long eccentricity cycles derived from the Wufeng Formation and the lower Longmaxi Formation in well N211. The established paleoenvironmental mode of the Upper Yangtze Region provides a basis for discussing the effects of a long eccentricity cycle on organic matter enrichment. As mentioned in Section 5.2, the paleoenvironmental mode during the deposition of the black shales in the Wufeng-Longmaxi Formation in the Upper Yangtze Region was related to the glacial period. The development of upwelling in the Upper Yangtze Region promoted the enrichment of organic matter in shales. Strong upwelling occurred in the colder period of the Earth in the Upper Yangtze Region (Figures 8A-C). Therefore, the black shales of the Wufeng-Longmaxi Formation deposited when Earth was cold are richer in organic matter ( Figures 8A-C).
During the deposition of the Wufeng-Longmaxi Formation black shales, periodic changes in long eccentricity led to periodic changes in the amount of solar radiation received by Earth, which in turn led to changes in cold and warm conditions (or changes in glacier size) in long eccentricity cycles (Zou et al., 2016;Cao et al., 2021). According to the paleoenvironmental model during the deposition of the Wufeng-Longmaxi Formation black shales in the Upper Yangtze Region, the changes in cold and warm conditions in long eccentricity cycles lead to periodic changes in the strength of upwelling in the Upper Yangtze Region and then affect the enrichment of organic matter in the Wufeng-Longmaxi Formation black shales. During long eccentricity minima, glaciers were formed at a relatively large scale (Zou et al., 2016;Cao et al., 2021); therefore, the offshore winds in the Upper Yangtze Region near the equator were relatively strong, resulting in the formation of relatively strong upwelling currents and a relatively high TOC content in the black shales, which is reflected by their relatively high GR values (Figures 6, 9A). During long eccentricity maxima, glaciers were formed at a relatively small scale (Zou et al., 2016;Cao et al., 2021); therefore, the offshore winds and upwelling currents in the Upper Yangtze Region were relatively weak, resulting in relatively low TOC contents in the black shales, which is reflected by relatively low GR values (Figures 6, 9B). Analysis of the controlling effects of the long eccentricity cycle on organic matter enrichment in the black shales of the Wufeng-Longmaxi Formation reveals a higher level of organic matter enrichment in the black shales deposited during long eccentricity minima. This understanding can provide a basis for predicting the distribution of shale with high organic matter abundance.

CONCLUSION
1) The black shales in the Wufeng Formation and the lower Longmaxi Formation in the Upper Yangtze Region are abundant in biogenic siliceous matter and framboidal pyrite. 2) TOC content, GR value, and U content increase gradually from the bottom of the Wufeng Formation to the Guanyinqiao Member, peak at the Guanyinqiao Member, and then decrease gradually going further up to the lower Longmaxi Formation. 3) Long eccentricity cycles can be identified from the GR curve for the Wufeng Formation in well N211 in the Upper Yangtze Region, and they correspond to sedimentary cycle thicknesses of 1.03-2.29 m. Approximately six long eccentricity cycles are recorded in the Wufeng Formation. Long eccentricity and short eccentricity cycles can be identified from the GR curve for the lower Longmaxi Formation. The long eccentricity cycle corresponds to sedimentary cycle thicknesses of 4.83-5.64 m, while the 125-and 95-kyr short eccentricity cycles correspond to sedimentary cycle thicknesses of 1.47 and 1.16 m, respectively. Approximately seven long eccentricity cycles are recorded in the lower Longmaxi Formation.
4) The black shales of the Wufeng-Longmaxi Formation in the Upper Yangtze Region deposited when Earth was relatively cold are richer in organic matter. At long eccentricity minima, the offshore winds in the Upper Yangtze Region were relatively strong, resulting in the formation of relatively strong upwelling currents, which led to a relatively high TOC content in the black shales. At long eccentricity maxima, the upwelling currents were weaker, resulting in lower TOC content in the black shales.

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.

AUTHOR CONTRIBUTIONS
SZ completed the writing and editing of the article. YL designed the research and focused on the scientific issues involved in the article. ZN provided the shale samples. YX, LC, and LZ performed the data analyses and prepared the figures and tables. All authors contributed to the article and approved the submitted version.