Sequence Stratigraphy, Depositional Environment and Associated Lithofacies of Lacustrine Shale: A Case From the Upper Fourth Member of Shahejie Formation, Dongying Depression, Bohai Bay Basin

The lithofacies assemblages in the lacustrine shale of the Shahejie Formation in the Dongying Depression have different compositions due to the complex depositional climate, lake water conditions, and sediment supply. Hence, the lithofacies types and depositional conditions were analyzed based on parasequence division and lithofacies classification. According to gamma ray integrated prediction error filter analysis (INPEFA), handheld x-ray fluorescence (XRF) measurements, mineral composition, total organic carbon (TOC) content, and sedimentary structures, six parasequences (P1–P6 from bottom to top) were identified, and seven lithofacies are recognized: (1) organic-poor thin-bedded calcareous mudstone, (2) organic-poor laminated calcareous mudstone, (3) intermediate-organic laminated calcareous mudstone, (4) organic-rich laminated calcareous mudstone, (5) organic-rich calcareous shale, (6) organic-rich laminated carbonate-bearing mudstone, and (7) laminated silty mudstone. When the climate became warm and wet, the water level of the lake increased with a decrease in salinity, and the corresponding lithofacies assemblages in the P1 to P2 parasequences gradually changed from organic-poor thin-bedded and laminated calcareous mudstone to intermediate-organic and organic-rich calcareous shale. When the climate became humid and warmer, resulting in an increasing water level and terrigenous input, the lithofacies assemblages in the P3–P4 parasequences changed from organic-rich calcareous shale and laminated carbonate-bearing mudstone to laminated silty mudstone and calcareous mudstone. At the end, the climate became slightly drier and cooler with a decrease in the water level, and the P5 and P6 parasequences were dominated by intermediate-organic and organic-poor laminated calcareous mudstones. Through elemental analysis and lithofacies observation, we found that the redox conditions of the water may have changed over time. Among the parasequences, the water column of the P2 parasequence was mainly stratified by salinity, that of the P4 parasequence was mainly stratified by temperature, and that of the P3 parasequence was in transition. The synergistic effect of various factors played a key role in the enrichment of organic matter. This study highlights the depositional processes of the lithofacies assemblages in a lacustrine shale and contributes to the understanding of shale oil accumulation in lacustrine shale.

The lithofacies assemblages in the lacustrine shale of the Shahejie Formation in the Dongying Depression have different compositions due to the complex depositional climate, lake water conditions, and sediment supply. Hence, the lithofacies types and depositional conditions were analyzed based on parasequence division and lithofacies classification. According to gamma ray integrated prediction error filter analysis (INPEFA), handheld x-ray fluorescence (XRF) measurements, mineral composition, total organic carbon (TOC) content, and sedimentary structures, six parasequences (P1-P6 from bottom to top) were identified, and seven lithofacies are recognized: (1) organic-poor thin-bedded calcareous mudstone, (2) organic-poor laminated calcareous mudstone, (3) intermediate-organic laminated calcareous mudstone, (4) organic-rich laminated calcareous mudstone, (5) organic-rich calcareous shale, (6) organic-rich laminated carbonate-bearing mudstone, and (7) laminated silty mudstone. When the climate became warm and wet, the water level of the lake increased with a decrease in salinity, and the corresponding lithofacies assemblages in the P1 to P2 parasequences gradually changed from organic-poor thin-bedded and laminated calcareous mudstone to intermediate-organic and organic-rich calcareous shale. When the climate became humid and warmer, resulting in an increasing water level and terrigenous input, the lithofacies assemblages in the P3-P4 parasequences changed from organic-rich calcareous shale and laminated carbonate-bearing mudstone to laminated silty mudstone and calcareous mudstone. At the end, the climate became slightly drier and cooler with a decrease in the water level, and the P5 and P6 parasequences were dominated by intermediate-organic and organic-poor laminated calcareous mudstones. Through elemental analysis and lithofacies observation, we found that the redox conditions of the water may have changed over time. Among the parasequences, the water column of the P2 parasequence was mainly stratified by salinity, that of the P4 parasequence was

INTRODUCTION
With the successful development of marine shale gas in North America, shale oil and gas have become important unconventional energy sources Jarvie et al., 2007;Aplin and Macquaker, 2011;Huang et al., 2020;Kang et al., 2019). Nevertheless, the evolution characteristics and distribution patterns of continental shale lithofacies types are still not well known due to their complexity and variability. Relatedly, the sedimentary environment and evolution process of continental shale remains unclear. In recent years, many scholars have performed studies on shale lithofacies division, sequence division, sedimentary environment analysis, diagenetic analysis and so on (Loucks and Ruppel, 2007;Rodriguez and Philp, 2010;Abouelresh and Slatt, 2012;Hao et al., 2013;Tian et al., 2021;Wu et al., 2019;Zeng et al., 2022). Although the sedimentary sequence and lithofacies development characteristics of shale are relatively complex, they are controlled by the evolution of the sedimentary environment, showing obvious segmented characteristics as a whole, and the lithofacies association changes regularly (Abouelresh and Slatt, 2012;Zhang et al., 2020a). Therefore, it is important to clarify the sedimentary evolution of mud shale to establish a high-precision sequence stratigraphic framework and then carry out the fine characterization of lithofacies types and analysis of lithofacies association characteristics.
For the division of continental shale lithofacies, the method involving four components and three endmembers has been mainly adopted at present. However, due to different angles of consideration, there are differences in lithofacies division in different research areas, as well as in the same research area (Dong et al., 2015;Liang et al., 2018a;Chen et al., 2008Wu et al., 2017Li et al., 2017;Zhang et al., 2018). It is considered the key to establish a reasonable lithofacies division scheme to consider the distribution pattern of rock components and carefully observe the characteristic differences in sedimentary structures through fine observation of core thin sections.
In addition, it is very difficult to study the stratigraphy of shale with a uniform structure deposited in relatively deep water, which prevents the application of traditional sequence stratigraphy in shale (Boulila et al., 2014). Some scholars have proposed solutions to the sequence division of shale, such as the gamma ray parasequences (GRP) method established by Slatt et al. (2008;Slatt et al., 2012), which is the most widely used method; this method mainly divides GRP cycles based on the superposition combination pattern of gamma ray (GR) curves. However, due to the frequent changes in GR curves, it is difficult to divide systems tracts and quasi-sequences and consider the existence of many subjective factors, leading to inconsistent results in the division of the same stratum and work area (Chen et al., 2008;Liang et al., 2018b). Therefore, finding a method that can objectively reflect the boundary of sequence division more intuitively and obviously is a difficult problem in the division of mud shale sequences.
The upper fourth member of the Eocene Shahejie Formation (Es4s) is a typical example of a Bohai Bay Basin lacustrine shale and is predominantly composed of carbonate, clay minerals and other mineral components that vary more widely than those of marine siliciclastic shales. Based on the division of a highprecision sequence stratigraphic framework, we clarify the characteristics of lithofacies assemblages at different stages, discuss the response of lithofacies to the evolution of the sedimentary environment, explore the development and distribution pattern of the dominant lithofacies assemblages and guide shale oil exploration (Yu et al., 2022;Gao, 2021;Wang et al., 2020).

GEOLOGICAL SETTING
The Dongying depression is located in the Jiyang Depression of the Bohai Bay Basin. It is a continental lake basin with an asymmetrical half-graben structure ( Figure 1). The Dongying depression can be divided into four depressions: the Boxing depression, Niuzhuang depression, Lijin depression and Minfeng depression. The Dongying depression is approximately 90 km long from east to west and 65 km wide from north to south. The exploration area in the Dongying depression is approximately 5,760 km 2 . The Paleogene strata are widely distributed; from bottom to top, they are the Kongdian Formation, Shahejie Formation, Dongying Formation, Guantao Formation and Minghuazhen Formation. The Shahejie Formation is subdivided into four members, Es4, Es3, Es2, and Es1 (from bottom to top). The Es4 strata are subdivided into the Upper Es4 (Es4s) interval and the Lower Es4 (Es4x) interval ( Figure 1). The Es4s lake basin was in the prime stage of rift expansion, with strong subsidence and rapid expansion, forming a broad semideep to deep lake environment, and in this lake, lacustrine organic-rich black shale up to 200 m thick was deposited in the center of the lake basin.

Core and Thin Section Observations
Data for this study were collected from nearly 160 m of cores from Well NY1 in the Dongying depression, which penetrates the entirety of Es4s in the Shahejie Formation ( Figure 1). Almost all publications about Es4s in the Shahejie Formation in the Dongying depression are based on this well. A total of 62 core samples were collected from the Es4s shale target interval for observation and testing. Thin sections with a thickness of 0.03 mm were prepared and observed under a polarizing Zeiss microscope (Axio Scope A1). These samples were also examined using a field emission scanning electron microscope (FE-SEM) from Carl Zeise (Supra 55).

Energy-Dispersive X-Ray Fluorescence Analysis
The elemental compositions (23 elements per sampling site) of core samples were measured on clean and smooth core surfaces using a Thermo Scientific Niton XL handheld energy-dispersive Frontiers in Earth Science | www.frontiersin.org April 2022 | Volume 10 | Article 906987 3 X-ray fluorescence (ED-XRF) analyzer from the Reservoir Geology Key Laboratory of Shandong Province (East China), Qingdao, China, with a spatial resolution of 3-5 cm. The ED-XRF instrument has a built-in calibration system that measures in ppm and has been shown to provide a reliable first-hand view of the chemical stratigraphy of borehole cores (Rowe et al., 2012).

X-Ray Diffraction Analysis
X-ray diffraction (XRD) analysis was performed using an X'pert-Multipurpose diffractometer (MPD) from the Shengli Oilfield Geological Research Center, which consists of a copper butt, 30 kV pipe pressure, 40 mA conduit flow, and 2θ/min scanning speed.

Chemical Composition of Organic Matter
Total organic carbon (TOC) and hydrocarbon potential analyses were also performed on the samples. TOC was measured with a CS344 carbon and sulfur detector (Leco, United States), and free hydrocarbon analysis (e.g., S1; Mg HC/g rock) and kerogen cracking (S2; Mg HC/g rock) were performed using a Rock-Eval pyrolysis system (Rock-EVAL6). All analyses were carried out in the Shengli Oilfield Geological Research Center.

Integrated Prediction Error Filter Analysis
Spectral trend attribute analysis is a technique for identifying stratigraphic cycles by means of spectral analysis of logging or time series data (Nio et al., 2005). Based on the theory of cycle stratigraphy, the logging curve was converted into an integrated prediction error filter analysis (INPEFA) curve using digital signal processing, which reflects the hidden periods in the formation characteristic logging data. According to the principle of logging geology, different logging curves reflect different geological characteristics. Compared with other logging curves, the natural GR curve is the most sensitive to changes in mud content and is not sensitive to borehole conditions. It is a good index for analysis of the sedimentary environment. Thus, it is the most effective index to use for stratigraphic sequence research.

Mineralogical Characteristics
The Es4s shale are mainly composed of high percentages of calcite, clay minerals and quartz, while the contents of dolomite, pyrite, plagioclase and other minerals are low (Liang et al., 2018a). The mineral compositions of shale from Es4s in the Shahejie Formation in Well NY1 are shown in Figure 2 and Table 1. The content of carbonate minerals is high in Es4s: the content of calcite is low and that of dolomite is high at the bottom of the section, and the content of calcite increases and that of dolomite decreases upward. Clay minerals and quartz are abundant in the range of 3,360-3,390 m, with average contents of 35 and 25%, respectively. The quartz content increases significantly at 3,390 m, with an average content of 33%. The feldspar content is low in Es4s, with an average content of 6%.
The pyrite content is generally low, with an average content of 3%. The anhydrite content is low, and anhydrite is found only at the bottom of the section.

Organic Characteristics
The organic matter types were classified using rock pyrolysis data (HI versus Tmax values, Van Krevelen, 1961;Tissot et al., 1974). The diagram shows that the organic matter in Es4s samples is dominated by type I kerogen and type II kerogen ( Figure 3A). The vitrinite reflectance values range from 0.65% to 0.77%, and most samples are mature and plot in the oil generation window.
The TOC content of Es4s samples ranges from 0.58% to 1.57% (average 2.4%). S1 ranges from 0.54 mg/g to 9.62 mg/g, with an average of 4.2 mg/g, while S2 ranges from 1.99 mg/g to 31.52 mg/ g, with an average of 16 mg/g. The HI values range from 350-700 mg HC/g TOC, with an average of 584 mg HC/g TOC, and the Tmax values range from 418 to 447°C, with an average of 442°C. Proportional to water depth and reducing conditions (Andrea, 2010;Slatt et al., 2012), the relative hydrocarbon generation potential (RHP; that is, S1 + S2/TOC) ranges from 4.36 to 8.95 mg HC/g TOC (with an average of 7.7 mg HC/g TOC) (Figures 2, 3). The results show that the shale samples are rich in organic matter and high in TOC. S1 is positively correlated with TOC ( Figure 3B). The cross plot of S1 and TOC shows that most of the samples have indigenous hydrocarbons and a small part has nonindigenous hydrocarbons (Jarvie and Baker, 1984;Li et al., 2017b). The relationship between S2 and TOC is a useful indicator for evaluating the hydrocarbon generation potential of source rocks (Peters, 1986;Peters and Cassa, 1994;Alias et al., 2012;Shalaby et al., 2012;Zhao et al., 2014). According to the ratio diagram of S2 to TOC, the NY1 shale sample has good hydrocarbon generation potential ( Figure 3C).

Lithofacies Description
Based on previous works on lithofacies classification (Dong et al., 2015;Liang et al., 2018b;Chen et al., 2017;Zhang et al., 2018), the mineral compositions (three end members of clay, carbonates, and quartz and feldspar), TOC contents (<2.0% as organic-poor, 2.0-3.5% as intermediate-organic, and >3.5% as organic-rich) and sedimentary structure characteristics (thin-bedded, laminated, and shaly structures), were considered in this study, with seven different lithofacies types identified in the Es4s shale samples. The detailed characteristics of these lithofacies are shown in Figure 4 and Table 2.
Organic-rich calcareous shale: The hand samples are dark gray, rich in carbonate clay minerals and organic matter, and have a foliation structure ( Figure 4A). The organic matter content is greater than 3.5%, the carbonate content is generally greater than 75%, and calcite veins are present. The laminae are well developed and clearly bounded, consisting of light and dark laminae. The light laminae are calcite veins, and the dark laminae are organicrich clay minerals ( Figures 4A,B), with abundant framboidal pyrite ( Figure 4N). This lithofacies is generally developed in the middle and lower parts of Es4s.
Organic-rich laminated calcareous mudstone: The hand samples are dark gray but lighter in color than the organic-rich calcareous Frontiers in Earth Science | www.frontiersin.org April 2022 | Volume 10 | Article 906987 shale and rich in carbonate clay minerals and organic matter ( Figure 4C). The content of organic matter is greater than 3.5%, and the content of carbonate is greater than 75%, but it is generally lower than that in organic-rich calcareous shale. Granular calcite grains of different sizes are present ( Figure 4D) with obvious laminae and boundaries, and framboidal pyrite is developed. This lithofacies is mainly developed in the middle and lower parts of the upper Fourth Member of the Shahejie Formation. Organic-rich laminated carbonate-bearing mudstone: The hand samples are dark gray and rich in clay minerals ( Figure 4E). The content of organic matter is greater than 3.5%, the content of clay minerals is generally between 60 and 70%, and the carbonate mineral content is low. Granular calcite laminae are mainly developed in the form of interbedded organic matter and clay minerals, and framboidal pyrite is present ( Figure 4F). This lithofacies is mainly developed in the middle of the upper Fourth Member of the Shahejie Formation.
Intermediate-organic laminated calcareous mudstone: The hand samples are gray in color and are rich in carbonate and clay minerals ( Figure 4G). The organic matter content varies between 2 and 3.5%, carbonate content is greater than 65%, and micritic calcite laminae are well developed but variable ( Figure 4H).
Organic-poor laminated calcareous mudstone: The hand samples are gray, but the color is lighter than that of the intermediate-organic laminated calcareous mudstone. These samples are rich in carbonate and clay minerals ( Figure 4I). The content of organic matter is less than 2%, and the carbonate mineral content is generally greater than 50%. The laminae are slightly developed, but the boundaries are not obvious ( Figure 4J).
Organic-poor thin-bedded calcareous mudstone: The hand samples are gray, lighter in color than all other lithofacies, and rich in carbonate minerals and clay minerals ( Figure 4K). The content of organic matter is less than 2%, and the content of carbonate minerals is greater than 50%, but the content of dolomite is significantly higher than that of calcite. The dolomite is distributed in clay minerals in the form of micrite ( Figure 4O), and laminae are not developed ( Figure 4L).
Laminated silty mudstone: This lithofacies is rich in felsic minerals and clay minerals and generally has high organic matter contents (3.5% on average). Compared with the other lithofacies, this lithofacies has the highest contents of felsic minerals (up to 45%) and pyrite. This lithofacies is well laminated and characterized by interbedded organic matter and clay minerals ( Figure 4M). Granular quartz is usually surrounded by clay minerals and mainly exists in the form of terrigenous detrital material. This lithofacies is mainly developed in the middle and upper parts of Es4s.

Major and Trace Element Geochemistry
Suitable environmental indicators can better reflect the evolution of the paleoenvironment. The chemical index of alteration (CIA) is often used as an index to evaluate the intensity of weathering and the climate (Price and velbel, 2003). The CIA data in this paper are from Chao, Liang (2018). Generally, the element enrichment coefficient is used to indicate whether the element content in sediments is enriched or depleted relative to the standard content (Taylor and Mclennan, 1985;Walsh et al., 1988;Murray et al., 1992;Saito, 1992;Algeo and Liu, 2020). In addition, since the TOC content is usually closely related to  redox conditions, TOC can also be used as an indicator to reconstruct the redox conditions (Dong et al., 2018). The strontium abundance and Sr/Ba ratio can be used as sensitive indicators of paleosalinity. Therefore, the Sr abundance and Sr/Ba ratio recorded in sediments are significantly positively correlated with paleosalinity (Zheng and Liu, 1999). Siliceous organisms and authigenic quartz are scarce in the saltwater basin, so the content of felsic minerals can be used as a good indicator to evaluate terrigenous input. The Fe/Mn ratio can be used as a marker of offshore distance (or water depth) (Chen et al., 2008). It has been confirmed that Ba has a positive correlation with surface productivity and organic matter enrichment in water and has been widely used in evaluations of paleoproductivity (Dymond et al., 1992;Li et al.,. 2015;Teng, 2006). The results show that the Sr/Ba ratio of the tested shale samples varies between 0.25 and 13.23, with an average of 2.88. The Mn/Fe ratio varies between 0.02 and 0.14, with an average of 0.06. In addition, Mn/Fe and Sr/Ba ratios vary greatly. The Sr/Ba ratio of the middle and lower strata of Es4s (mean value: 7.33) is much higher than that of the upper strata of Es4s (mean value: 0.55). The Mn/Fe ratio (mean value: 0.04) is the lowest in the lower strata of Es4s, followed by the Mn/Fe ratio (mean value: 0.055) in the upper strata and the Mn/Fe ratio (mean value: 0.075) in the middle strata ( Figure 5 and Table 3). Figure 5 shows that the Mn/Fe and Sr/Ba ratios in the upper submember of the Fourth Member of the Shahejie Formation are negatively correlated; from bottom to top, the Mn/Fe ratios tend to increase gradually, although the ratios tend to decrease slightly in the upper strata. The Sr/Ba ratio gradually decreases.

Evolution of Sedimentary Environments Under a High-Precision Sequence Stratigraphic Framework
The development of mineral characteristics and sequence stratigraphy is controlled by the sedimentary environment. Therefore, it is very important to analyze the evolution of the FIGURE 3 | Organic geochemical characteristics of the Es4s shale. (A) A plot of the HI of the sample relative to the maximum pyrolysis temperature (Tmax) of the pyrolysis product, reflecting the kerogen mass and thermal maturity stage (boundary modified by Mukhopadhyay et al., 1995). (B) Free hydrocarbons (S1) versus total organic carbon (TOC). (C) Pyrolytic hydrocarbon (S2) and TOC plots after kerogen cracking, showing source rock potential (boundaries revised by Peters and Cassa, 1994  Frontiers in Earth Science | www.frontiersin.org April 2022 | Volume 10 | Article 906987 8  with GRP theory to divide the sequences of the Barnett shale in the United States and proposed that the RHP could also be used as an index for sequence division. The INPEFA method proposed by a research team in the Netherlands mainly uses GR curves to analyze spectrum attributes and obtain an INPEFA curve (Nio et al., 2005), thus realizing the division of sequence stratigraphy (Liu et al., 2019;Ma et al., 2016;Ma et al., 2021). Based on the INPEFA of the GR curve and combined with the analysis of environmental indicators, we divided the upper submember of the Fourth Member of the Shahejie Formation from bottom to top into six parasequences ( Figure 5). In P1 (3,444.17-3,470.15 m depth range), the GR curve has a low INPEFA curve, and the CIA index is low, ranging from 73.715-76.346 (average 75.035), indicating that the climate was cold and dry, and the water body was shallow, corresponding to a low Mn/Fe ratio. The Sr/Ba ratio is high, with an average value of 8.43, indicating a high salinity, which may be related to the abundant gypsum at the bottom. The low TOC value in this section is 1.75% on average, and the RHP value is 2.87 mg HC/g TOC, indicating that the reducing conditions in this section were weak, which is not conducive to the preservation of organic matter. In addition, this result is consistent with the low enrichment coefficient of U and the lowest Bio-Ba value, indicating that the primary productivity was very low. In this period, the climate was relatively dry and cold, with little precipitation and shallow water, leading to weak reducing conditions, high salinity, a lack of biological development, low primary productivity and a low organic matter content.
In P2 (depth range 3,423.80-3,444.17 m), the INPEFA curve of the GR curve shifted significantly to the right, the CIA index increased significantly, and the increase in the Mn/Fe ratio indicates that the climate became warm and humid during this period. The water body became significantly deeper, and the Sr/Ba ratio decreased greatly, with an average value of 5.19, indicating that the salinity decreased. The development of framboidal pyrite in this period and the high TOC (average 2.5%) and high RHP (average 7.76 mg HC/g TOC) in this section indicate that the water conditions in this period were highly reducing, conducive to the preservation of organic matter, and that the Bio-Ba values were low, indicating low primary productivity. In this period, the climate was wet and warm, and the precipitation increased, leading to the deepening of the lake and a low salinity.
In P3 (depth range 3,393.34-3,423.80 m), the INPEFA curve of the GR curve continued to shift to the right, and the CIA index and Mn/Fe ratio continued to increase, indicating that the climate continued to become warmer and more humid and that the water depth further increased. The Sr/Ba ratio decreased to an average of 4.26, and the salinity continued to decrease. The average TOC was 2.7%, and the average RHP was 8.27 mg HC/g TOC, indicating that the water conditions in this period were highly reducing and conducive to the preservation of organic matter. In addition, the Bio-Ba value increased, indicating that the primary productivity increased in this period.
In P4 (depth range 3,362.70-3,393.34 m), the INPEFA curve of the GR curve has the largest rightward deviation, and the CIA index is the highest, with an average of 80.603, indicating that the water body reached its deepest depth at this time. The Sr/Ba ratio further decreases, with an average value of 2, indicating that the salinity of water at this time was significantly lower than that in P3. The enrichment coefficient of U is lower than that in P3. The mean TOC is 2.53%, and the RHP is 7.45 mg HC/g TOC, indicating strong reducing conditions in this period but weaker than those in P3. The maximum value of Bio-Ba and the maximum value of felsic minerals indicate the maximum primary productivity and terrigenous detrital input during this period. During this period, the climate was very wet and warm, and the precipitation was abundant, which allowed the water depth to reach the maximum; the salinity was further reduced, terrigenous debris was abundant, the primary productivity was high, and the water conditions were strongly reducing, which led to enrichment of organic matter in the sediment.
In P5 and P6 (depth range 3,336.81-3,336.81 m), the INPEFA curve of the GR curve began to deviate to the left, the CIA index decreased slightly, with an average of 80.60, and Mn/Fe decreased slightly, indicating that the water depth decreased at this time. The Sr/Ba ratio decreased with an average value of 0.97, indicating that the salinity of water further reduced the TOC by 2.00% and RHP by 6.17 mg HC/g TOC, indicating that the reducing conditions of the water weakened further in this period, but the reduction was stronger than in P1. During this period, the Bio-Ba value decreased, and the felsic mineral content decreased overall, indicating that primary productivity decreased and terrigenous detrital input decreased. During this period, the climate became slightly drier and cooler, and the water body decreased in size. However, the input of terrigenous debris and the existence of primary productivity during this period reduced the salinity of the water body, weakened the reducing conditions of the water body, and led to a decrease in organic matter abundance.
In general, the climate during the deposition of Es4s shifted from relatively dry to wet to a slightly drier and cooler. The salinity gradually decreased, the TOC and U element enrichment coefficients first decreased and then increased, and the reducing conditions of the water first decreased and then increased.

Corresponding Characteristics of Lithofacies Assemblages in Different Sedimentary Environments
The regular change in sedimentary environment controls the regular change in mineral composition and sedimentary structure, which leads to the regular change in the vertical lithofacies assemblage (Liang et al., 2018a;Zeng et al., 2021;Wu et al., 2020;Zhang et al., 2022b). The climate of the P1 parasequence was dry and cold, with little precipitation, resulting in shallow water, easy destruction of laminar structure and weak reducing conditions, which were not conducive to the preservation of organic matter. In addition, gypsum minerals developed at the bottom, indicating a higher salinity and more carbonate minerals in the area, with an average of approximately 40%. However, the content of dolomite is generally higher than the calcite content, which is caused by a dry and cold climate. In summary, the influence of the above factors led to the development of a lithofacies combination involving organic-poor thin-bedded calcareous mudstone and organic-poor laminated calcareous mudstone. During the P2 parasequence period, the climate became very wet and warm, and the rainfall increased. As a result, the water deepened rapidly, the salinity decreased, the reducing conditions increased, the organic matter was well preserved, and the laminar structure was well developed. In addition, carbonate minerals were dominant in this period, and the calcite content was higher than the dolomite content, so the lithofacies assemblage of organic-rich calcareous shale and organic-rich laminated calcareous mudstone was predominant in this period. As the climate became wetter and warmer, the water became deeper, leading to a further decrease in salinity, stronger reducing conditions, increased input of terrigenous detritus, and increased primary productivity. The organic matter was preserved and abundant in the P3 parasequence. The lithofacies assemblage is mainly composed of organic-rich laminated carbonate-bearing mudstone and organic-rich laminated calcareous mudstone. In the P4 period, the climate became humid and warm, and the precipitation increased further, which further deepened the water body, reduced the salinity, and strengthened the reducing conditions. Terrigenous detritus, primary productivity and organic matter were abundant in this period. However, terrigenous detritus also disrupted the reducing conditions to a certain extent. The lithofacies in this stage were mainly composed of laminated silty mudstone and organic-rich laminated calcareous mudstone. During the parasequence periods of P5 and P6, the climate became slightly drier and cooler, and the water depth and reducing conditions decreased. Moreover, terrigenous clastics were more abundant during this period, diluting the water salinity to a certain depth. The lithofacies were mainly laminated calcareous mudstones with intermediate and low amounts of organic matter ( Figure 6).

Effects of Paleoclimate on Lake Water Properties at the Parasequence Level
The paleoclimate plays a decisive role in the enrichment and preservation of organic matter by controlling the depth of the water body, influencing the supply of terrestrial matter, and controlling the primary productivity of the water body, all of which influence the redox conditions (Liang et al., 2018b;Zhang et al., 2022c). The data show that the CIA is positively correlated with the Mn/Fe ratio and negatively correlated with the Sr/Ba ratio ( Figures 7A,B), indicating that the paleoclimate had a significant influence on the paleosalinity by controlling the paleowater depth. A dry and cold climate leads to stronger evaporation and less precipitation, which increases the salinity of lake water. As the climate becomes wetter and warmer, evaporation weakened, and precipitation increases, leading to a decrease in salinity. In addition, a high salinity increases the ion concentration of the water, which promotes carbonate precipitation, and a humid climate leads to a greater supply of terrigenous material. The data show that felsic minerals are positively correlated with carbonate minerals ( Figure 7C), which indicates that terrigenous input can dilute the salinity of water and indirectly proves that a humid climate can reduce the Frontiers in Earth Science | www.frontiersin.org April 2022 | Volume 10 | Article 906987 13 salinity of water. The data also show that the Mn/Fe ratio is positively correlated with TOC during the P1 and P2 periods ( Figure 7D), and the Sr/Ba ratio is positively correlated with TOC ( Figure 7E), indicating that salinity controlled the stratification of the water and enhanced the reducing conditions of the water, which was conducive to the preservation of organic matter. Moreover, the P2 parasequence mainly developed unique vein calcite lithofacies, indicating that the content of carbonate deposited at that time was higher than that of the other parasequences, and the reducing conditions were also stronger. In the P3 parasequence period, the Mn/Fe ratio is negatively correlated with TOC ( Figure 7F), but the Sr/Ba ratio is positively correlated with TOC ( Figure 7G), indicating that during this period, due to the humid and warm climate, the lake water deepened, and the salinity further decreased, resulting in a low salinity and slightly weak reducing conditions; however, water stratification still existed, the salinity was still high, and the reducing conditions remained strong. As the climate gradually became wetter and warmer, precipitation continued to increase, and the salinity further decreased. In the P4 parasequence period, the TOC is still high, but the Sr/Ba ratio has a very weak positive correlation with TOC ( Figure 7K), while the CIA value has an obvious positive correlation with TOC ( Figure 7I), indicating that the salinity concentration of water at this time was low, and it may not have been sufficient to stratify the lake water. At this time, the lake was mainly stratified by temperature. There is also a significant positive correlation between felsic minerals and TOC ( Figure 7H), indicating that the wet and warm climate led to abundant terrigenous input, which had a significant dilution effect on water salinity and disrupted the salinity-driven stratification. However, there was no correlation between the Bio-Ba value and the felsic mineral content at this time ( Figure 7J), indicating that the strong primary productivity at this time was due to algae blooms rather than terrigenous input. Moreover, the crossplot in Figure 4 also proves this point. During the P5 and P6 periods, the climate became slightly drier and cooler, the depth of the water decreased slightly, and the stratification of water gradually weakened, which was not conducive to the enrichment of organic matter. The ratio of Sr/Ba also gradually decreased, and there was a negative correlation between felsic minerals and the Sr/Ba ratio ( Figure 7L), indicating that although terrigenous input was not large at this time, it had a dilution effect on water salinity, leading to a decrease in water salinity. In summary, the paleoclimate influenced the salinity and redox environment of the water by controlling the water depth and input of terrigenous debris and finally caused corresponding regular changes in the water properties.

Effects on Organic Matter Enrichment
The factors affecting the enrichment of organic matter in lake water mainly include the redox conditions, primary productivity and input of terrigenous debris (Liang et al., 2017b;Zhang et al., 2022d;Zhang et al., 2020b). During the evolution process of sedimentary environments, these factors jointly affect the enrichment and preservation of organic matter, but they often play different roles, with one or more factors playing a major role.
When the climate is dry and cold, the water is shallow and less reducing, which is not conducive to the preservation of organic matter. When the climate is warm and humid, the water is deep, and if the salinity is higher, then salinity-driven stratification will provide stronger reducing conditions. If terrigenous detritus is abundant and the water salinity is low, then temperature-driven stratification will provide strong reducing conditions. In addition, if algal blooms are common, the high primary productivity will also be conducive to the enrichment of organic matter. Therefore, the level of organic matter content is the result of many interrelated factors rather than being affected by a single factor. Therefore, for different periods, a detailed discussion of the changes and relationships of various factors is the main process of analyzing the causes of organic matter enrichment in a given period.
During deposition of the P1 parasequence, the TOC was low, with an average value of 1.75%, and the RHP was 2.87 mg HC/g TOC. At this time, the climate was dry and cold, and the water was shallow, leading to weak reducing conditions and high water salinity, which were not suitable for the survival of plankton. Moreover, the living space of plankton was limited, resulting in a low initial productivity. These results indicate that the enrichment of organic matter in this period was mainly controlled by the redox conditions. The TOC was high (mean 2.5%), and the RHP was also high (mean 7.76 mg HC/g TOC) in the P2 parasequence period. During this period, the climate became very wet and warm, the depth of lake water increased greatly, and the salinity was high, forming salinity-driven stratification, which strengthened the reduction of organic matter and was conducive to the preservation of organic matter. The enrichment of organic matter in this period was mainly controlled by the strong reducing conditions resulting from the water stratification caused by salinity.
In the P3 parasequence, the TOC and RHP averaged 2.7% and 8.27 mg HC/g TOC, respectively, indicating the enrichment of organic matter. At this time, the climate became more humid and warmer, the depth of the lake further increased, and the salinity was further reduced. The water body still retained strong reducing conditions. In addition, planktonic algae were more common in this period, and the primary productivity was higher, which provided a source of organic matter to a certain extent. These phenomena indicate that the organic matter enrichment in this period was mainly controlled by the strong reducing conditions caused by salinity-driven stratification, followed by primary productivity.
In the P4 parasequence period, the degree of humidity and warmth reached the maximum, the water depth also reached the maximum, the average TOC was 2.53%, the average RHP was 7.45 mg HC/g TOC, and the organic matter was relatively abundant. At this time, terrigenous detrital input was abundant, which diluted the water salinity to a certain extent and greatly reduced the water salinity concentration, and the salinity-driven stratification almost disappeared. Due to the very warm and humid climate, high water temperature and great depth of water, temperature-driven stratification occurred in this period, resulting in strong reducing properties. In addition, planktonic algae blooms occurred in this period, and the primary productivity was very high, so the organic matter was relatively abundant. Organic matter enrichment in this period was mainly controlled by the strong reducing conditions and primary productivity due to stratification caused by temperature, followed by terrigenous detrital input. Moreover, the TOC and RHP values in this period were lower than those of P2 and P3, indicating that salinity-driven stratification was better than temperature-driven stratification in the enrichment and preservation of organic matter, which may be related to the fact that an appropriate saline environment can effectively promote organic matter flocculation and thus improve the capture efficiency of organic matter (Hu et al., 2021).
In the P5 and P6 parasequence periods, the average TOC and RHP values were approximately 1.8% and 6.60 mg HC/g, respectively, which were both low. At this time, the climate became dry and cold, the depth of water decreased, the reducing conditions of the water were weakened, and the living space of plankton was restricted. During this period, some terrigenous debris was imported, which further diluted the salinity of the water and further reduced its salinity.

CONCLUSION
The upper fourth member of the Eocene Shahejie Formation is a typical lacustrine shale in the Bohai Bay Basin. It is mainly composed of carbonate and clay minerals, and the mineral composition is quite different from that of marine siliceous shale. Seven lithofacies types were classified based on core thin section and scanning electron microscope observations, combined with X-ray diffraction and geochemical analysis data. The organic matter in the shale in this area is mainly type I kerogen, with TOC values ranging from 0.58% to 4.57%, with an average of 2.41%, and RHP values ranging from 4.36 mg HC/g TOC to 8.95 mg HC/g TOC, with an average of 7.21 mg HC/g TOC. We introduced the INPEFA curve of the GR curve to divide the studied strata into parasequences, and a total of six parasequences were identified. According to the paleoenvironment indices and lithofacies characteristics, we establish an evolution model diagram of the lacustrine mud shale sedimentary environment and lithofacies assemblage. The results show that the regular change in sedimentary environment controlled the regular change in mineral composition and sedimentary structure, which led to the regular change in vertical lithofacies assemblage. The paleoclimate played a decisive role in the enrichment and preservation of organic matter by controlling the water body depth and subsequently influenced the supply of terrestrial matter, the primary productivity of the water body, and the redox conditions. The enrichment of organic matter was influenced by many interrelated factors rather than a single factor.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.