Abstract

The Late Permian on the periphery of the Songliao Basin, eastern Inner Mongolia, is an important hydrocarbon source rock system. Its sedimentary environment plays an important role in the evaluation of hydrocarbon prospects in the area. Unfortunately, until now, the interpretation of the sedimentary environment of this area has been controversial. We investigated the Late Permian sedimentary environment by studying the sedimentary characteristics and geochemistry. Based on these investigations, we conclude that the Linxi Formation is mainly composed of clastic sediments, interbedded with limestone lenses, with bioherm limestone at the top of the formation. Inner-layer marine fossils (calcium algae, bryozoans, and sponges) and freshwater and blackish water microfossils (bivalves) are all present, indicative of a typical shallow water sedimentary environment with an open and concussion background. In terms of geochemistry, the formation is relatively light rare Earth enriched, with significant positive Eu anomaly, slight positive La and Y anomaly, weak positive Gd anomaly, and lack of Ce anomaly. The average B/Ga ratio of the mudstone is greater than 3.3, and the average Sr/Ba ratio of the limestone is greater than 1.0. The range of the 87Sr/86Sr ratio is from 0.707285 to 0.707953. The range of δ13C values is from −4.0‰ to 2.4‰. The sediment assemblages, rare Earth elements, trace elements, and 87Sr/86Sr and C isotopes of the formation indicate that the Linxi Formation formed in a marine sedimentary environment and occasional marine-terrestrial transitional facies. The formation can be further divided into littoral facies, neritic facies, bathyal facies, and delta front.

1. Introduction

The eastern part of Inner Mongolia in China is located on the junction of the Sino-Korea Plate and the Siberian Plate. During the Permian period, this area experienced an important tectonic transition, and a complex tectonic paleogeographic environment was formed within the sedimentary basin. Since the 20th century, different interpretations have been put forth concerning the paleogeographic environment of the Linxi area of eastern Inner Mongolia during the Late Permian. He et al. interpreted the sedimentary environment of the Linxi Formation, according to the lithology, sedimentary structure, and paleontology of an exposed section in Linxi County. They concluded that the middle-lower parts of the Linxi Formation deposited in the marine environment and the upper part of the Linxi Formation deposited in the terrestrial lake environment [1]. Zhang et al. and Tian et al. concluded that the Linxi Formation formed in a marine sedimentary environment, based on the presence of marine fauna and bioherms in the exposed sections of the Late Permian Linxi Formation in the Linxi area [2, 3]. The other interpretation was that the Late Permian Linxi Formation deposited in terrestrial lake sedimentary environment, according the development of the freshwater bivalve fossils and Angara flora of the Late Permian Linxi Formation in the Linxi area [46].

The Late Paleozoic Permian strata in eastern Inner Mongolia are a potential hydrocarbon (oil and gas) resource on the periphery of Songliao Basin. Clarifying the evolution of the sedimentary environment of the strata is important to locate an Upper Paleozoic hydrocarbon generation horizon and predict hydrocarbon prospects in the area. Meanwhile, the lithofacies of the Linxi Formation record the evolution from filling to closing of the Xing-Meng Trough, which is located on the boundary of the Sino-Korea and Siberian Plates. Although previous studies formed some conclusions concerning the Permian sedimentary environment of eastern Inner Mongolia, especially in the Linxi area, some issues still require clarification. This is mainly due to the fact that the fossil types are not consistent throughout the formation and that determining the sedimentary environment of a complete set of formations based on environmentally meaningful typical fossils, sedimentary structures, or bioconstructions has its limitations. This research study investigates the Late Permian Linxi Formation’s sedimentary environment, by analyzing the sedimentary characteristics and the geochemical characteristics of the Linxi Formation in the Linxi area, eastern Inner Mongolia.

2. Regional Geological Survey

2.1. Structure Distribution

The eastern part of Inner Mongolia is located in the Xing-meng orogenic zone. It is bordered in the west and north by the Mongolia-Okhotsk suture zone, which is adjacent to Mongolia and Russia. In the south, it is bordered by the Xar Moron River-Yanji suture zone, which is adjacent to the North China plate. In the east, it is bordered by the Central Sikhote subduction belt near the Pacific Plate. The Linxi area is located east of the Xing-meng orogenic belt and north of the Xar Moron River fault in the area where the Sino-Korea and Siberian Plates meet (Figure 1(a)). The Linxi Formation is mainly distributed in the Keshiketeng Banner and Linxi County. It is an important Permian lithostratigraphic unit in the Upper Permian strata of eastern Xing-meng. The underlying Middle Permian Wujiatun Formation is composed of marine carbonates, volcanic rocks, and clastic sediment, and the overlaying Lower Triassic Xingfulu Formation is composed of terrestrial clastic sediments [7].

2.2. Lithologic Characteristics

The Late Permian Linxi Formation can be divided into 5 sections, which we call sections Lin1-5. Section Lin1 is composed of medium-thick gray and gray-green medium-fine-grained sandstone, thin-bedded siltstone, and interbedded siltstone (A1 and A2 in Figure 2), which is locally interbedded with thin-medium calcarenite or limestone lenses. Section Lin2 is composed of dark gray, grayish-black shale, silty mudstone, and silty slate interbedded with gray-green fine sandstone, and it locally contains visible slump breccia and limestone lenses (A3 and A4 in Figure 2). Section Lin3 is composed of medium-thick light brown medium-coarse grained sandstone interbedded with thin gray, gray-brown siltstone, argillaceous siltstone, and silty mudstone and local areas of medium-thick gravel coarse grained sandstone and limestone (A5 in Figure 2). The entire section from bottom to top evolved from sandstone and pebbled sandstone to middle-thick sandstone and upper interbedded sandstone and mudstone. Section Lin4 is composed of dark gray, gray-black shale, silty mudstone, and argillaceous siltstones interbedded with thick gray-green, gray fine grained sandstone, thin-medium limestone, and limestone lenses (A6 and A7 in Figure 2). The lower part of this section is primarily composed of argillaceous siltstone and silty mudstone. The middle part of the section is dominated by shale, while the upper part of the section is composed of interbedded sandstone and mudstone interbedded with limestone lenses. Section Lin5 is composed of gray argillaceous siltstone, sandy mudstone, light gray medium-fine grained sandstone, tuffaceous sandstone, massive gray bioherm limestone, and framestone (A8 in Figure 2).

2.3. Sedimentary Structures and Paleontology Fossils

The Linxi Formation contains many types of sedimentary structures and paleontology fossils. In the lower part of section Lin1, wedge-shaped cross bedding, tabular cross bedding, and symmetrical ripple bedding are present in the sandstones, while horizontal bedding (B1 in Figure 2), current bedding, and small-scale cross bedding are present in the mudstone-siltstones and argillaceous siltstones. The upper part of section Lin1 sandstones contain parallel bedding, oblique bedding, small-scale cross bedding, and current bedding (B2 in Figure 2), indicating a shallow water concussion environment, containing bivalves Palaeanodonta sp., Palaeomutella sp., Anthraconaia sp., and fragments of plant fossils. The lower part of section Lin2 contains wash-out structures (B3 in Figure 2), locally developed slump structures, deformation bedding (B4 in Figure 2), and small current bedding, indicating a relatively deep-water environment. The imbricated structure observed in the lower part of section Lin3 exhibits parallel bedding and tabular cross bedding (B6 in Figure 2) in the sandstones and current bedding, concussion ripples, bioturbation structures (B5 in Figure 2), and mud cracks in the mudstones (B7 in Figure 2), which is indicative of an intermittently exposed shallow water environment. The middle and upper parts of section Lin3 exhibit lenticular foreset bedding and local contemporaneous slump structures and contain bivalve fossils (Palaeodoronella soronensis, Palaeodonta sp., P. cf. trigonalis, P. cf. semilunulata, and P. hahaiensis) and plant fossil debris, which is indicative of a delta front sedimentary environment. In section Lin4, low angle oblique bedding, graded bedding, concretionary structure (B8 in Figure 2), and small-scale ripple bedding developed can be seen. The clastic rocks of section Lin4 contain conchostracas (Huanghestheria sp., Palaeolimnadia sp., and Pemphicyclus cf. trochoides), bivalve fossils (Palaeanodonta sp., P. cf. trigonalis, and Palaeomutela soronensis), plant fragments, and brachiopod fossil fragments. Bryozoans, sponges, and calcareous fossils, including Monostysisyrinx gen. and Archaeolithoporella gen., can be seen in the limestone, indicating an expansive coastal environment that may have mixed a terrestrial source. The lower part of section Lin5 exhibits wedge-shaped cross bedding and parallel bedding, while the middle and upper parts exhibit visible bioturbation structures, concretionary structures, low angle oblique bedding, and small-scale cross bedding. Bryozoans and sponge bone fossils are present in the limestone [7], indicating a neritic environment.

3. Samples and Analytical Procedures

The strata discussed in this paper and the samples used for geochemical analysis were primarily collected from the Guandi profile in Linxi County in eastern Inner Mongolia (Figure 1(b)). A total of 33 samples were analyzed. All the samples are fresh and uncontaminated rocks by stripping a weathered layer. Observed under the microscope, the samples were without epigenic dikes or strong recrystallization. Of these samples, 4 of them were limestone of Early Triassic, which were analyzed for Sr-C isotopic ratios. 24 of the samples were primarily composed of lenticular limestones from the Linxi Formation. They were analyzed for major, trace, and rare Earth elements and Sr-C isotopic ratios (Tables 13). Seven of the samples were mudstone, which were analyzed for major, trace, and rare Earth elements (REEs) (Table 1). The X-ray fluorescence spectroscopy (XRF) was used to determine the oxides of major elements. The LOI (loss on ignition) of the samples was recorded as weight loss after burning at 1000°C. Trace and rare Earth elements were determined by an Element 6000 inductively-coupled plasma mass spectrometer (ICP-MS). Both of the tests were performed in the Chinese Academy of Geological Sciences National Geological Experimental Testing Center following the method described by Ryu et al. [15]. The Sr isotopes were analyzed using a British VG354 isotope mass spectrometer following the method described at the Nanjing Institute of Soil Science and Technology Service Center, Chinese Academy of Sciences. The measured value of 87Sr/86Sr for American isotope standard NBS987Sr was 0.710236 ± 7 (n = 10) with a precision of ±0.005%. The C isotopes were analyzed using MAT-252 isotope mass spectrometer with a precision of ±0.2‰ at Isotope Research Laboratory, Institute of Mineral Resources, Chinese Academy of Geological Sciences.

4. Results

4.1. Rare Earth Elements and Trace Elements

Trace element concentrations are reported in Table 1. REE + Y concentrations were normalized to the Post-Archaean Australian Shale (PAAS) composite [16]. The total REE content (ΣREE) of the sediments in the Linxi Formation ranges from 10.07 μg/g to 195.69 μg/g with an average of 76.56 μg/g. These samples of mudstones and lime mudstones with high REE are clearly contaminated by a considerable terrigenous debris component and perhaps masking the anomalies distribution of REE in marine environments. Thus, they are not used any further for palaeoenvironmental reconstructions. We mainly use the trace elements of the carbonates and argillaceous carbonates to discuss the palaeosedimentary environment of the Linxi Formation. Shale-normalised (SN) elemental anomalies were calculated on a linear scale, assuming that differences in concentration between neighbouring pairs are constant, as follows: La/La = La/(3Pr − 2Nd)PAAS, Ce/Ce = Ce/(0.5La + 0.5Pr)PAAS, Eu/Eu = Eu/(0.5Sm + 0.5Gd)PAAS, and Gd/Gd = Gd/(2Tb − Dy)PAAS.

The ΣREEs of the limestones are between 10.07 μg/g and 91.22 μg/g. Abbreviations are applied to light REE (LREE) and heavy REE (HREE). As can be seen from the PAAS normalized REE distribution pattern (Figure 3(a) and Table 1), the limestone samples of Linxi Formation in the Guandi profile of Linxi area display relatively uniform shale-normalised REE + Y patterns (Figure 3(a)). The LREE depleted relative to the HREE with LaSN/YbSN from 0.48 to 0.91, and the limestones display a significant positive Eu anomaly ((Eu/Eu)SN = 1.78). But significant differences can be noted between the limestones at the bottom of Lin1 and the other limestones of Linxi Formation. The limestones at the bottom of Lin1 show a slight positive La anomaly ((La/La)SN = 1.16), lack of Ce, Gd, and Y anomaly ((Ce/Ce)SN = 0.99), (Gd/Gd)SN = 1.0), and (Y/Ho)SN = 1.07). The other limestones of Lin1-5 show a distinct positive La anomaly ((La/La)SN = 1.6), slight positive Y anomaly ((Y/Ho)SN = 1.25), weak positive Gd anomaly ((Gd/Gd)SN = 1.08), weak negative Ce anomaly ((Ce/Ce)SN = 0.93, except sample LG 40-2), and lack of Ce anomaly ((Ce/Ce)SN = 0.95) (Table 1).

In addition to REE analysis, other trace element (Sr, Ba, B, and Ga) analyses were conducted on the limestones and mudstones from the Guandi profile in the Linxi area (Table 1). In the Linxi Formation, the Sr content ranges from 58.9 μg/g to 3052 μg/g with an average of 918.43 μg/g. The Ba content ranges from 22.2 μg/g to 591 μg/g with an average of 238.18 μg/g. The Sr/Ba ratios of the limestone range from 1.76 to 23.68 with an average of 8.18. The B content of the mudstones ranges from 18 μg/g to 150 μg/g with an average of 57.21 μg/g, which is enriched compared to the average B content of the crust (7.6 μg/g). The Ga content ranges from 4.34 μg/g to 22.1 μg/g with an average of 12.95 μg/g. The B/Ga ratios of the mudstone range from 1.90 to 7.94 with an average of 4.42.

4.2. Sr and C Isotope Distribution

The whole rock total Sr content and Sr isotopic data for the Linlin Formation in the Linxi area are presented in Table 2. The C isotopic data for the Linlin Formation are presented in Table 3. The results show that the Sr contents of all nine limestone samples in the Linxi Formation are greater than 850.0 μg/g with a maximum value of 2762.0 μg/g and an average of 1435.1 μg/g. The 87Sr/86Sr values of the limestone samples range from 0.707285 to 0.707953 with an average of 0.707624. The δ13C values of the limestone samples vary from −4.0‰ to 2.4‰, and the average is −1.3‰. As shown in Figure 4, there is no correlation between Sr-C isotope and Mn/Sr, and the Mn/Sr ratio reaches the upper limit of 2-3 for qualified samples recommended by Kaufman et al. [27]. Seven samples have low Mg/Ca ratios (<0.1) and may be presumed to be unaltered. The Sr contents of all of the limestone samples are significantly larger than the minimum Sr content (400 μg/g) suggested by Korte et al. in their Sr isotope deposition studies [28]. In addition, there is no correlation between 87Sr/86Sr and Si (representing SiO2 content) (Figure 4(b)), indicating that the Sr isotopic composition is completely derived from the endogenous sedimentary limestone and is not affected by other Si-bearing terrestrial debris and smaller diagenesis events. This indicates that the Sr–C isotopes in the limestone represent the Sr-C isotopes of the sedimentary fluids.

5. Discussion

5.1. REE Distribution and Sedimentary Environment

The REEs of different water bodies are quite different. The chemical composition of the authigenic minerals and rocks formed by these water bodies can be used to determine the geochemical characteristics of modern and ancient sedimentary media [2932]. The water involved in the formation of carbonate rocks mainly fall into three categories: seawater, rivers or atmospheric precipitation, and hydrothermal fluid. The behavior of REE in a given type of water body results in distinguishing characteristics. In seawater, the LREE is being depleted relative to the HREE and results in La enrichment [17]. In addition, the REE of seawater in oxygen environment is resulting in Ce depletion [33]. In rivers, lakes, and other freshwater medium, Ce usually exhibits a slight anomaly or no anomaly. The infiltration of weathered leaching water near Eu-rich rocks or Eu bonded to iron hydroxide in water also causes slight positive Eu anomalies [22]. Beyond that, an influx of high-temperature midocean ridge oceanic hydrothermal fluids can also cause positive Eu anomalies.

In order to facilitate a more direct comparison of the REE contents of the different water bodies, PAAS was used to normalize the REE of the ancient marine carbonates and minerals, modern seawater, lakes, lacustrine carbonate rocks, and hydrothermal fluids [17, 18, 2025] (Figures 3(b)3(h)). Among the rare Earth patterns of the different sediments, water, and minerals, normal seawater has enriched La, negative Ce anomalies, and slightly enriched Gd (Figures 3(b) and 3(c)). The hydrothermal fluids from the median ridge have obvious positive Eu anomalies (Figure 3(g)). The freshwater and the lacustrine carbonate rocks are characterized by a relatively flat REE patterns without obvious anomaly distribution (Figures 3(e) and 3(h)). The gulf and coastal waters have positive Ce anomalies and small negative Eu anomalies and show relatively gentle HREE distribution patterns (Figure 3(d)). The marine carbonate rocks and minerals have enriched La, small negative Ce anomalies, and positive Eu anomalies (Figure 3(f)). The samples from a distal section of near-shore marine carbonates display a distinct positive Y anomaly, positive La anomaly, and slight positive Gd anomaly. In contrast, the samples from the proximal section of near-shore marine carbonates display a slight positive Y anomaly, weak positive Gd anomaly, and a lack of La anomaly (Figure 3(i)). It can be noted by the distribution characteristics of REE, except of the Eu, the limestones of Linxi Formation in the Guandi profile of Linxi area display relatively uniform shale-normalised REE + Y patterns (Figure 3(a)). The limestones at the bottom of Lin1 display a slight positive La anomaly and lack of Ce, Gd, and Y anomaly, which are closed to the REE distribution patterns of lake water, no-marine carbonates, and the proximal section of near-shore marine carbonates (Figures 3(e), 3(h), and 3(i) and 5). The middle-upper limestone samples of Lin1 and the limestones of Lin2-5 display a distinct positive La anomaly, slight positive Y anomaly, weak positive Gd anomaly, and lack of Ce anomaly, and its REE patterns resemble those deposits from the distal section of littoral environments, which is dominated by the ocean but with freshwater inputs from continental weathering (Figures 3(f) and 3(i)) or the neritic environments with a terrestrial source mixing. In addition to the above REE of La, Ce, Gd, and Y, almost all the limestones of the Linxi Formation display a distinct positive Eu anomaly. The Eu positive anomalies have no obvious indication of the sedimentary environment. With the lack of Eu positive anomaly in mudstones and no obvious correlation between the concentration of Eu and iron oxide (not shown), we infer that the reason for the positive Eu anomaly is related to the hydrothermal fluid mixing.

5.2. Trace Element Distribution and Sedimentary Environment

In addition to REE, other trace element ratios of sedimentary rocks can also be used to determine their sedimentary environments. The Sr in fine-grained marine sediments tends to be enriched relative to Ba, causing the higher Sr/Ba value in the fine-grained sediments. Therefore, the relative contents of Sr and Ba recorded by the fine-grained sediments, especially in mudstones and shales, i.e., Sr/Ba, are positively correlated with the paleosalinity [3436]. Wei and Algeo conducted a comprehensive study of the Sr/Ba ratios as paleosalinity proxies in ancient shales and mudrocks, and they concluded that the Sr/Ba ratio of freshwater is less than 0.2, the Sr/Ba ratio from 0.2 to 0.5 indicate brackish environments, and ratios greater than 0.5 represent the marine environments. When water depth increases in sedimentary basins, the energy of seawater decreases with the increased clay minerals and fine-grained material. The fine-grained sediment enhances the adsorption of Ba ions and causes a lower Sr/Ba ratio. In general, the Sr/Ba ratios is lower in deep sea or bathyal sediments [36]. In addition, the incorporation of terrigenous debris in the offshore area can also result in a decrease in Sr/Ba ratios of argillaceous sediments. We compared the Sr and Ba contents of the Permian sedimentary rocks in the other areas of China (Figure 6 and Table 4) [34, 37, 38]. The Sr/Ba ratios of the marine mudstones of Ordovician Pingliang Formation in southern Ordos Basin range from 0.19 to 2.59 with an average of 0.96. The Sr/Ba ratios of the marine mudstones of the lower Permian Liangshan Formation in Dushan area of Guizhou Province range from 0.14 to 2.79 with an average of 0.91. The Sr/Ba ratios of the mudstones of marine-continental transitional facie of Permian Shanxi Formation in Ordos Basin range from 0.13 to 0.89 with an average of 0.36. The Sr/Ba ratios of the mudstones and lime mudstones in the Linxi Formation of the study area range from 0.1 to 3.09 with an average of 0.93. The Sr/Ba ratios of the mudstones and lime mudstones of Linxi Formation correspond well with the Sr/Ba ratios of the marine mudstones and the mudstones of marine-continental transitional facie in the other areas of China. The Sr/Ba ratios indicate that the Linxi Formation limestone mainly formed in a salt water environment and develop bathyal sea sedimentary environment or have been affected by terrestrial debris (the Sr/Ba ratios of sample LG 29-1H <0.2).

In addition, B is relatively enriched in marine shales, whereas Ga is relatively enriched in terrestrial shales. Therefore, the B/Ga ratio can be used to distinguish between marine and terrestrial sedimentary environments to a certain degree [36]. Yan et al. conducted a comprehensive study of the Mesozoic biological, mineral, and geochemical characteristics of the Paleocene-Pliocene strata of the Northern Jiangsu basin [39]. Based on a statistical comparison between ancient American marine sediments, modern marine sediments, ancient terrestrial sediments, modern freshwater sediments, modern Chinese marine sediments, and modern lacustrine sediments, they concluded that the B/Ga ratio of terrestrial sediment is less than 3.3, and ratios greater than 3.3 represent sea-land transitional zones and marine sedimentary environments. To more intuitively contrast the content of B and Ga and assess the deposition environment, the data of B and G contents of shale in different regions of the world have been counted and are presented in Table 5 and shown in Figure 6. The analysis indicates that the B and Ga content of the Changxingjie marine shale in southern China varies from 34.7 μg/g to 141.0 μg/g and from 6.8 μg/g to 25.4 μg/g, respectively, and the average value of B/Ga is 4.79 [40]. The B and Ga content of the Permian terrestrial mudstone in Zhungeer Basins of China varies from 14.2 μg/g to 20.0 μg/g and from 4.4 μg/g to 18.5 μg/g, respectively, and with the average values (B/Ga) of 2.27 [37]. The B and Ga content of the ancient freshwater argillaceous sediments in the United States varies from 8.0 μg/g to 55.0 μg/g and from 9.0 μg/g to 25.0 μg/g, respectively, and with the average values (B/Ga) of 2.02 [41]. The B and Ga content of the ancient marine argillaceous sediment in the United States varies from 65.0 μg/g to 187.0 μg/g and from 10.0 μg/g to 39.0 μg/g, respectively, and with the average values (B/Ga) of 5.77 [41]. The B and Ga content of the inner shelf mudstone off the east coast of Inida varies from 80.0 μg/g to 130.0 μg/g and from 11.0 μg/g to 33.0 μg/g, respectively, and with the average values (B/Ga) of 4.84. The B and Ga content of the outer shelf mudstone off the east coast of Inida varies from 90.0 μg/g to 160.0 μg/g and from 16.0 μg/g to 34.0 μg/g, respectively, and with the average values (B/Ga) of 4.81 [42]. In the study area, the B and Ga content of Linxi Formation mudstone varies from 33.4 μg/g to 150.0 μg/g and from 7.7 μg/g to 22.1 μg/g, respectively. The B/Ga ratios of Linxi Formation mudstone varies from 1.90 to 7.94 with an average of 4.62 being greater than 3.3, except for the mudstone (Sample LG 75-4, B/Ga, 1.9) at the top of section Lin5, and this corresponds well with the marine mudstone from the other area of the world (Figure 7). The relative B and Ga contents of the mudstone from the Linxi Formation indicate a predominantly marine sedimentary environment during deposition of the Linxi Formation and that the formation was affected by terrestrial debris at the end of the deposition of section Lin5.

5.3. Strontium Isotope Distribution and Sedimentary Environment

The stable evolutionary characteristics of Sr isotopes provide a basis for the comparison of Sr isotopes in seawater of the same geological age around the world [43]. In terms of the various sources of Sr, due to the weathering and leaching of crustal Si-Al rocks, lakes, rivers, and groundwater from terrestrial sources transport a large amount of radiogenic 87Sr, which is generated in their catchments (basin), and thus, they have high 87Sr/86Sr values (about 0.7115) [44]. However, seawater transports Sr into the mantle at subduction zones, and Sr enters the seawater at midoceanic ridge hydrothermal systems [45]. The seawater also interacts with terrestrial waters that drain into the oceans, resulting in the 87Sr/86Sr values of modern seawater being approximately 0.7092 [46]. Carbonates and some fossils can record the composition of Sr isotopes in sedimentary media. The composition and evolution of Sr isotopes (87Sr/86Sr values) throughout geologic history have been well documented and reported. Therefore, by comparing the Sr isotopic compositions of endogenous carbonate rocks that have not been altered to the Sr isotopic composition of seawater during that geologic period, we can determine whether the sedimentary strata was deposited in a marine sedimentary environment or a terrestrial sedimentary environment. As shown in Figure 7, Korte et al. established the background value (0.7073–0.7074) of 87Sr/86Sr during the Permian and the Triassic by the analysis of Sr isotopes from brachiopods [9, 10, 28]. Tian and Zheng found that 87Sr/86Sr values of Permian marine limestone range from 0.7072 to 0.7077 in South China [8]. Huang et al. found that 87Sr/86Sr values of the Late Permian seawater range from 0.7067 to 0.7076 by analyzing the Sr isotopes of marine carbonate rocks from the Upper Permian in the Upper Yangtze region [11]. Martin and Macdougall researched the contents of Sr isotopes in conodonts of Permian/Triassic from the United States and Pakistan, and they found that 87Sr/86Sr values of seawater in Permian/Triassic range from 0.7069 to 0.7089 (Table 2) [12]. The results of Sr isotopic analysis show that the 87Sr/86Sr values of the Linxi Formation limestone in the Linxi Guandi profile range from 0.70728 to 0.70795 and distribute within the range of Sr isotope composition of Permian seawater (Figure 8). The 87Sr/86Sr values of the four limestone from section Lin1 to the bottom of section Lin3 range from 0.70753 to 0.70776 and fall within the Sr isotope range of the Late Permian Wujiaping and Changxing Period marine carbonates in southern China, indicating a shallow water carbonate sedimentary environment. The 87Sr/86Sr values of five limestone samples from sections Lin4 and Lin5 range from 0.70729 to 0.70795. Comparing with the contents of Sr isotope of the Permian marine limestone in other area of China, the 87Sr/86Sr values of the lower and upper limestone of section Lin4 have abnormally high values of 0.70795 and 0.70780, respectively, reflecting the fact that terrestrial water or terrestrial debris mixed with the main sediment supply. In addition, the 87Sr/86Sr values of one limestone sample each (sample numbers LG72-13T and LG73-2CT) from the top of section Lin4 and the bottom of section Lin5 are consistent with the 87Sr/86Sr background values (0.7073–0.7074) of the global Late Permian paleocean and exhibit relatively low 87Sr/86Sr values. The interpretation of this phenomenon is that the upper part of the Linxi Formation was deposited during a stage of rift closing and uplift accompanied by intense volcanic eruptions. The eruptions of magmas from a mantle source affected the composition of the Sr isotopes, which resulted in a decrease in the 87Sr/86Sr value of the sediments. The tuff interlayer in section Lin5 also supports the occurrence of volcanic activity.

5.4. C Isotope Distribution and Sedimentary Environment

The compositions of carbon isotopes in carbonate rocks and their changes can effectively trace sedimentary environments, diagenetic environments, and other geological processes [14, 4752]. The δ13C values of all nine limestone samples in the Linxi Formation range from −4.0‰ to 2.4‰ with a variation range of 6.4‰, and 67% of them are in the range of δ13C values (0 ± 2‰) of normal seawater, indicating a normal marine sedimentary environment in the Linxi area. Compared with the δ13C of marine Permian carbonate rocks (e.g., eastern Inner Mongolia, Changxing District of Zhejiang and Khorat Basin of Thailand) and Triassic lacustrine limestone in Linxi area (Table 3), the δ13C values of limestone in the Linxi Formation are in the range of marine limestone but significantly higher than those of the lacustrine limestone (Figure 9). However, it should be noted that individual sample with a low δ13C value indicates a further differentiation in the sedimentary environment under normal marine background. For instance, the difference of seawater circulation, atmospheric circulation, and the accordingly production rate of organic carbon may cause changes in carbon isotope composition in limestone. The extensive development of plant debris and bioturbation structures in the Linxi mudstone indicate a shallow water environment with sufficient biomass during the deposition period of the Late Permian Linxi Formation. And so, we infer that the 12C produced by the biological respiration, the oxidation of organic matter, and the freshwater recharge of low-carbon isotopic compositions all contribute to the low carbon isotope value of diagenetic fluids. Moreover, frequent volcanic activity was triggered by strong tectonics at the Late Permian in the study area and the volcanism produced a lot of CO2. The carbon isotopic compositions of CO2 (−5‰) may also be one of the reasons for the lower δ13C value of limestone in the study area.

5.5. Sedimentary Environment and Its Evolution

These features of trace elements, isotopes, and sedimentary characteristics indicate that the environment of Linxi Formation is closed to the element and isotope compositions characteristic of littoral facies and neritic facies but with freshwater inputs from continental weathering.

Although the characteristics of limestone or shale of Linxi Formation yielded very consistent results, significant differences can be noted between 5 sections of Linxi Formation (Lin1–5) (Figures 59). Through the minor differences of trace compositions, lithologic association, sedimentary structures, and paleontology fossils of the Linxi Formation, we made the sedimentary environment dividing in detail. The lithology of section Lin1 of the Linxi Formation is mainly composed of sandy and silty sediments that exhibit symmetrical ripple bedding, horizontal bedding, and other sedimentary structures with directional significance for hydrodynamic conditions. The limestones at the lower section of Lin1 display a slight positive La anomaly, lack of Ce, Gd, and Y anomaly, and the upper section evolves into the characteristic of distinct positive La anomaly and slight positive Y anomaly. These distribution features of REE combining with Sr-C isotope and trace element indicate that the sedimentary environment of section Lin1 evolved from the proximal section of littoral facies with terrigenous clastic mixing to the distal section of littoral facies. Compare with the Lin1 section, the mud content of section Lin2 increased obviously which composed of mudstone, silt mudstone, and shale. The lower part of section Lin2 exhibited slump structures, convolute bedding, and flat bedding; the upper parts developed flat bedding, bioturbation structures, ripple mark structures, and current bedding. These sedimentary structures combining with the distribution of Sr-C isotope and Sr/Ba and B/Ga ratio indicate that the water body deepens during the depositional phase of section Lin2, and the sedimentary environments are interpreted to be bathyal slope to neritic facies. The sandstone content of the section Lin3 increases relatively. The lower parts of section Lin3 exhibited tabular cross bedding, wedge-shaped cross bedding, flat bedding, and ripple mark structures, and these sedimentary structures combining with the distributions of REE, C isotope, and Sr/Ba ratios of the limestones indicated that the sedimentary environment was close to the distal section of littoral facies. The middle and upper parts of section Lin3 developed foreset bedding, desiccation crack, bioturbation structure, slump structure, and current bedding, which indicated the sedimentary environment of delta front and littoral facies. The section Lin4 mainly composed of shale, silty mudstone, argillaceous siltstones, and limestone lenses. The geochemical characteristics of limestones with lower 87Sr/86Sr values, slight positive La and Y anomaly, and the sedimentary structures composing of flat bedding, tabular cross bedding, current bedding, bioturbation structure, and the typical marine fossils observed in this section indicate that the sedimentary environment of section Lin4 is mainly neritic facies. Section Lin5 is composed of sandstones and bioherm limestones and develops bryozoans and sponge bone fossils. The geochemical characteristics of the limestones in section Lin5 display a marine distribution pattern, which is displaying as low 87Sr/86Sr values and distinct positive La and Y anomaly. Combined with the sedimentary and geochemical characteristics, the sedimentary environment of section Lin5 is interpreted to be neritic facies.

6. Conclusions

The lithologic association of the Linxi Formation is mainly composed of clastic rocks, with interbedded carbonates and lenticular limestone bodies. The upper and lower units of the Linxi Formation limestones consistently show initial δ13C values, 87Sr/86Sr ratios, and the similar REE + Y patterns that are characterized by depletion of LREE relative to HREE, distinct positive La anomaly, minor positive Y anomaly, weak positive Gd anomaly, lack of Ce anomaly, and the same range of Sr/Ba and B/Ga ratio. In general, the sedimentary environment of the Linxi Formation of Late Permian is still dominated by an open marine environment and has freshwater inputs from continental weathering. The marine environments can be divided into 3 major types which are littoral facies, neritic facies, and bathyal facies (Figure 10). In addition, due to the regional uplift and the accompanying shallow water environment, a delta front environment of sea-land transition zone sediments also developed in the study area in addition to the marine sedimentary environment. The sedimentary environment of the Late Permian Linxi Formation underwent a process from littoral facies to neritic facies (bathyal facies) and then through littoral facies, delta front to neritic facies.

Data Availability

The raw data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was partially supported by the National Key R & D Program of China under grant 2017YFC0603103 and the National Natural Science Foundation of China (no. 41572098).