Upper Jurassic–Lower Cretaceous Source Rocks in the North of Western Siberia: Comprehensive Geochemical Characterization and Reconstruction of Paleo-Sedimentation Conditions

The present work is devoted to geochemical studies of the Bazhenov Formation in the north of the West Siberian Petroleum Basin. The object is the Upper Jurassic–Lower Cretaceous section, characterized by significant variations in total organic carbon content and petroleum generation potential of organic matter at the beginning of the oil window. The manuscript presents the integration of isotopic and geochemical analyses aimed at the evaluation of the genesis of the rocks in the peripheral part of the Bazhenov Sea and reconstruction of paleoenvironments that controlled the accumulation of organic matter in sediments, its composition and diagenetic alterations. According to the obtained data, the sediments were accumulated under marine conditions with a generally moderate and periodically increasing terrigenous influx. The variations in organic matter composition are determined by redox conditions and terrigenous input which correlate with the eustatic sea level changes during transgressive/regressive cycles and activation of currents. Transgression is associated with an intensive accumulation of organic matter under anoxic to euxinic conditions and insignificant influence of terrigenous sources, resulting in the formation of rocks with oil-generating properties. During the regression periods, the terrigenous sedimentation increased along with the dissolved oxygen concentration, and deposits with low organic matter content and gas-generating properties were formed.


Introduction
The Bazhenov Shale Formation of the West Siberian Petroleum Basin is one of the world's largest hydrocarbon (HC) reservoirs, and in recent decades it has attracted attention as an unconventional hydrocarbon source. The oil and gas production in the region is accompanied by the expansion of the fundamental knowledge on the genesis of the Bazhenov Formation, which is the main oil source of the Basin.
According to different estimations, the epicontinental Bazhenov Sea in the Upper Jurassic period occupied an area of up to 2 million km 2 [1,2], and its northern part was connected to the Boreal Sea through a relatively narrow channel [3]. Differences in sedimentation conditions and initial organic matter type in the central and peripheral areas led to the formation of the Bazhenov deposits with different lithological characteristics, total organic carbon (TOC) distribution, composition, and petroleum generation potential. The results of the studies of the Bazhenov Formation rocks and organic matter are reported in numerous publications, including the marine sediments of the central part (distal areas of Figure 1. Location of the studied area (modified after Zanin [4]) and stratigraphy of Mesozoic deposits (modified after Afanasenkov [14]); the investigated stratigraphic interval is marked by frame.

Geological Setting
The Bazhenov Formation is Upper Jurassic-Lower Cretaceous organic-rich shale, which mainly consists of biogenic silicites with clay components and carbonate admixture [24]. The Formation has an area of almost the whole of the West Siberian Basin, and is considered the major source rock in the region [1,25]. Siliceous mudstones of the Formation were deposited in the sea basin during the Tithonian-Lower Valanginian in the deep water depression in the territory of modern western Siberia ( Figure 1). The average thickness of the deposits varies between 20 and 50 m. The Bazhenov source rocks are underlayed by Abalak Formation (Callovian-Lower Tithonian) clay deposits also containing shale, carbonate rocks formed in shallow seas [26,27].
The Bazhenov Formation and its stratigraphic equivalents formed in the distal parts of the paleo-sea consist of siliceous rocks enriched in organic matter that Since the area was located at the border between the Bazhenov epicontinental and the Boreal Seas, the deposits were formed under frequently changing depositional conditions and the continuous influence of terrigenous influx. The sequence consists of alternating darker colored organic-rich and lighter organic-poor deposits, that differ from the Bazhenov Formation rocks of the central part of the Basin [13,23]. The mineral composition of rocks is characterized by high clay content, low carbonate content, and significant TOC variations.
The primary focus of this study is to investigate the key features of the different intervals, their genesis and to perform reconstruction of the sedimentation conditions. The conclusions are based on the results of a number of geochemical and isotope analyses of Geosciences 2021, 11, 320 3 of 28 organic and inorganic components of the rocks. We also investigated the relationships between paleo-sedimentation conditions and the present characteristics of the rocks and organic matter, which determine the prospects for oil and gas production in the region.

Geological Setting
The Bazhenov Formation is Upper Jurassic-Lower Cretaceous organic-rich shale, which mainly consists of biogenic silicites with clay components and carbonate admixture [24]. The Formation has an area of almost the whole of the West Siberian Basin, and is considered the major source rock in the region [1,25]. Siliceous mudstones of the Formation were deposited in the sea basin during the Tithonian-Lower Valanginian in the deep water depression in the territory of modern western Siberia (Figure 1). The average thickness of the deposits varies between 20 and 50 m. The Bazhenov source rocks are underlayed by Abalak Formation (Callovian-Lower Tithonian) clay deposits also containing shale, carbonate rocks formed in shallow seas [26,27].
The Bazhenov Formation and its stratigraphic equivalents formed in the distal parts of the paleo-sea consist of siliceous rocks enriched in organic matter that accumulated under anoxic marine conditions. Similar depositional conditions were present over the central part of the Basin, whereas the marginal paleo-sea areas were influenced by terrigenous influx and sea-level alterations. Therefore, the rocks in the central part of the Bazhenov are lithologically homogeneous over the large territory and form the largest oil shale formation in the world, both in thickness and area. Despite the high amount of organic matter in the deposits, no large oil fields are found in its central part due to the low permeability of rocks and preservation of liquid hydrocarbons in its pore space. The productive intervals are thus attributed to porous and fractured carbonates or diagenetic carbonates in radiolarites, which contain low amounts of syngenetic organic matter (kerogen and/or asphaltenes) [28]. During recent years, the Bazhenov Formation has been considered not only as a source rock but also as an unconventional reservoir. In contrast, the deposits in the peripheral parts were accumulated under the influence of various factors including the river and periodical flows, tectonic setting, connection with ocean waters, etc. Therefore, the deposits formed in the peripheral paleo-sea areas are depleted in organic matter and have higher terrigenous material and clay mineral content compared to that of the central part of the Basin.
The area investigated is located in the northern part of the West Siberian Basin close to the Arctic region. In terms of tectonics, the area is in the Yamal-Gydan megasaddle, at the boundary between the Yamal and Gydan regions ( Figure 1). There are widely known gas and oil fields in the region, such as Stormovoe, Gydanskoye, Salmanovskoye, and Tambeyskoye, and the Bazhenov Formation is considered the main source rock [17].

Samples and Methods
In the current study, a core from 1 well was examined, using a varying number of samples according to the applied technique. Details are given in the corresponding sections of the results.
The lithological study was performed at different scales, from core description to optical microscopy and scanning electron microscopy (SEM).
For the detailed study of micromorphological rock features, scanning electron microscopy (SEM) was used. Quattro S (Thermo Fisher Scientific, Waltham, MA USA) allowed us to analyze small (3-5 mm) rock samples with an electron beam current range from 1 pA to 200 nA with an accelerating voltage of 5 to 10 kV. Samples were coated with a thin (10 nm) gold layer to discharge samples.
Spectral gamma-ray logging of full-size core was performed on the USGC-1 gammaray spectrometer instrument. The vertical resolution was 5 cm.
The chemical composition of rocks was studied by X-ray fluorescence (XRF) using a portable X-MET 8000 Expert Geo spectrometer (Oxford Instruments Analytical GmbH, Uedem, Germany). The core samples were preliminary powdered. Determinations were made according to the settings of a universal calibration program used for a wide range of diagnostic chemical elements (from Mg to U), with a duration of 60 s. As a result, the quantitative analysis of petrogenic chemical elements (Mg, Al, Si, P, S, K, Ca, Ti, Mn, Fe) and of some trace elements (V, Ni, Cu, Zn, Rb, Sr, Zr, Mo, Ba, Th, U, etc.) was performed. This approach allows one to perform the continuous profiling of chemical composition along the interval, obtaining a large amount of geological data without special sample preparation and with an accuracy comparable with individual sample measurements (XRF, ICP-MS, carbonate content measurements). The obtained trends are reliable, since the maxima and minima found are consistent with stationary XRF and ICP-MS measurements performed for individual samples [30].
For Sm, Nd content and isotope ratio measurements, 20-30 mg of rock was crushed and dissolved in the mixture of concentrated acids HF and HNO 3 (5:1, 1 mL). According to the isotope dilution method, the tracer amount was calculated and added. The insoluble organic residue was additionally washed with concentrated HNO 3 and evaporated. The dilution with concentrated HCl and evaporation was repeated three times to obtain the solid matter, containing elements in the form of chlorides. The sample was then dissolved in 2.2 N HCl, and the insoluble residue was removed by centrifugation. The separation of Sm and Nd for isotope composition analysis was carried out by ion-exchange column chromatography. Isotope composition was measured using the "Triton" thermal ionization mass spectrometer (ThermoFisher Scientific, Waltham, MA, USA). To calculate the systematic error and normalize the obtained values, La Jolla standard was used, to identify the degree of contamination of a blank sample with a tracer was analyzed. The obtained 143 Nd/ 144 Nd ratios were corrected to the age of sedimentation (~145 Ma) using the formula:  [31,32]. We determined the main pyrolysis parameters: S0, S1, S2, S3, S4, S5, temperature Tmax ( • C), and calculated the values of total organic carbon (TOC, wt.%), generative organic carbon (GOC, wt.%) and non-generative organic carbon (NGOC, wt.%), Kgoc (GOC/TOC, %), hydrogen index (HI, mg HC/g TOC), oxygen index (OI, mg CO 2 /g TOC), index of productivity (PI), and carbonate carbon content (CC, wt.%).
The bulk kinetic studies of kerogen thermal decomposition were performed using the HAWK Resource Workstation. An extracted rock sample was subjected to nonisothermal pyrolysis at temperatures from 300 to 650 • C at three constant heating rates (3, 10, 30 • C/min) for calculation of kinetic parameters: frequency factor A and activation energy distribution. Hydrocarbon products were transferred from the oven to FID by a constant helium flow (50 mL/min). The obtained pyrograms were computed to discrete activation energy distributions using KINETICS2015 software version 6.020 (GeoIsoChem Corporation) for a fixed frequency factor A = 1 × 10 14 s −1 .
For Pyro(TD)-GC×GC-TOFMS/FID analysis, a powdered rock sample was placed in a quartz glass vial with a closed round end with a plug of quartz wool. Analyses were carried out using Pegasus 4D (LECO Corporation, St. Joseph, MI, USA) equipped with the pyrolysis injection unit (TDU2 and PYRO, Gerstel). The temperature of thermal extraction was 350 • C, and pyrolysis was performed at 500 • C. For thermal extraction analysis, the distribution of n-alkanes and a set of biomarker parameters were calculated, and for kerogen decomposition products, the temperature 500 • C and exposure time of 30 s were chosen and the main outcome was organic class distribution ( Figure 2).
Geosciences 2021, 11, x FOR PEER REVIEW 6 of 31 TOFMS electron ionization mode was at 70 eV, the ion source temperature was 250 °C, the voltage was 1600 V, the acquisition rate was 100 Hz and the mass range was 35-450 amu. The FID operation temperature was 350 °C, and the FID acquisition rate was 100 Hz. Data acquisition and processing were carried out using ChromaTOF™ software version 4.70.7.0 (LECO Corporation) and MS Office 2013. Individual component identification was based on the comparison of analysis results with the NIST database and literature data [33][34][35]. For elemental and isotope composition, DELTA V Plus mass spectrometer (Thermo) equipped with a Flash 2000 elemental analyzer was applied for isotope and elemental analysis of carbon, nitrogen, and sulfur. The accuracy of analytical procedures was established using international standards: oil NBS 22, strontium carbonate NBS 23, ammonium sulfate IAEA-N-2, and barium sulfate NBS 127 for carbon, nitrogen, and sulfur, respectively, as well as of the laboratory standards, whose elemental and isotope composition was determined during interlaboratory comparative measurements. The content and isotope composition of organic carbon (Corg) and carbonate carbon (Ccarb) were analyzed separately. The carbon and oxygen isotope composition of carbonates (δ 13 Ccarb and δ 18 Ocarb) were analyzed in gaseous CO2 released during the decomposition of carbonates after 100% orthophosphoric acid treatment. The amount of CO2 was quantified to obtain Ccarb content. The residual solution was filtered and the residue containing Corg. was analyzed with elemental analyzer using the combustion method. The relative error of element content measurements is ±10%. The measured C, H, N, and S isotopic ratios were recalculated and given in δ values in ‰, giving the deviation of the measured isotopic ratio relative to the isotopic ratio of the corresponding standard: PDB for carbon, SMOW for hydrogen and oxygen, AIR for nitrogen and CDT for sulfur. The absolute error of measurement for δ values was ±0.3 ‰ for carbon, and ±0.5 ‰ for sulfur and nitrogen.

Results and Discussion
This section presents the data on various inorganic and organic geochemical analyses: the description of the experimental results and their interpretation based on the published data related to the technique and object of the study. The final subsection is devoted to the integration of the various geochemical data and the reconstruction of sedimentation conditions to explain the alteration of rock composition, TOC distribution, and hydrocarbon generation potential within the studied section. Two-dimensional gas chromatography was performed using reverse column configuration: Rxi-17Sil MS 30 m and Rxi-5HT 1.5 m (Restek). Helium was used as the carrier gas at a flow rate of 1.4 mL/min. The oven temperature program started from 40 • C then ramped to 340 • C. Two detectors operated simultaneously with 150 s delay. TOFMS electron ionization mode was at 70 eV, the ion source temperature was 250 • C, the voltage was 1600 V, the acquisition rate was 100 Hz and the mass range was 35-450 amu. The FID operation temperature was 350 • C, and the FID acquisition rate was 100 Hz. Data acquisition and processing were carried out using ChromaTOF™ software version 4.70.7.0 (LECO Corporation) and MS Office 2013. Individual component identification was based on the comparison of analysis results with the NIST database and literature data [33][34][35].
For elemental and isotope composition, DELTA V Plus mass spectrometer (Thermo) equipped with a Flash 2000 elemental analyzer was applied for isotope and elemental analysis of carbon, nitrogen, and sulfur. The accuracy of analytical procedures was established using international standards: oil NBS 22, strontium carbonate NBS 23, ammonium sulfate IAEA-N-2, and barium sulfate NBS 127 for carbon, nitrogen, and sulfur, respectively, as well as of the laboratory standards, whose elemental and isotope composition was determined during interlaboratory comparative measurements. The content and isotope composition of organic carbon (C org ) and carbonate carbon (C carb ) were analyzed separately. The carbon and oxygen isotope composition of carbonates (δ 13 C carb and δ 18 O carb ) were analyzed in gaseous CO 2 released during the decomposition of carbonates after 100% orthophosphoric acid treatment. The amount of CO 2 was quantified to obtain C carb content. The residual solution was filtered and the residue containing C org. was analyzed with elemental analyzer using the combustion method. The relative error of element content measurements is ±10%. The measured C, H, N, and S isotopic ratios were recalculated and given in δ values in ‰, giving the deviation of the measured isotopic ratio relative to the isotopic ratio of the corresponding standard: PDB for carbon, SMOW for hydrogen and oxygen, AIR for nitrogen and CDT for sulfur. The absolute error of measurement for δ values was ±0.3 ‰ for carbon, and ±0.5 ‰ for sulfur and nitrogen.

Results and Discussion
This section presents the data on various inorganic and organic geochemical analyses: the description of the experimental results and their interpretation based on the published data related to the technique and object of the study. The final subsection is devoted to the integration of the various geochemical data and the reconstruction of sedimentation Geosciences 2021, 11, 320 6 of 28 conditions to explain the alteration of rock composition, TOC distribution, and hydrocarbon generation potential within the studied section.

Lithological Description
The samples are composed of argillaceous-siliceous rocks and argillites. The texture is pelitic, cryptocrystalline, and microcrystalline. The rocks have a relatively homogeneous mineral composition, mainly composed of clay minerals (average volume fraction is 42%) and of siliceous components (38% on average). The presence of feldspars (up to 10%), finely dispersed pyrite (up to 8%, in the form of framboids and pseudomorphs on organogenic detritus), and carbonates (mostly lower than 3%, with the exception of three layers with more than 40% carbonate) in the form of calcite and dolomite is detected. Little poorly preserved organogenic detritus of various origin is present (foraminifera shells, radiolarians, and bivalve mollusks).
Since rocks are characterized by relatively stable mineral composition, the main lithological types were distinguished according to rock texture and structure ( Figure 3). All established lithological types should be classified as argillaceous-siliceous mudstone according to the classification for fine-grained sedimentary rocks proposed by Lazar et al. [36].

Major and Trace Element Composition, Sm and Nd Isotopes
The elemental composition was determined for 131 samples taken at regular intervals along the section (the average distance between the samples is 15-25 cm) using portable XRF. In the resulting data set (Tables A1 and A2), significant variations in the elemental composition relative to the average values for the Abalak and Bazhenov Formations were identified. Therefore, we specified geochemical levels within the section.
Considering the characteristic paleo-geographic features of black shale formations [30,38], we paid special attention to the lithological composition and genesis of rocks: fluctuations in redox conditions, changes in the input of terrigenous and biogenic material, and bioproductivity of the paleo-sea.
According to the methods described in [39][40][41], we have calculated a set of lithochemical proxies, which were used to identify the genesis of rocks ( Figure 4). The most informative components to assess the contribution of terrigenous clastic components are Al, Ti, Zr, K, and the authigenic and biogenic components of the rock: S, P, Mo, Ca, Mg, Ni, Cu, Zn, V, etc. [41,42]. The carbonate content was studied via the ratios of Ca, Mg, Sr, Mn, and Fe. The predominant type of iron minerals (oxide, sulfide forms, the Fe affinity with S, Mn, and Ni) was revealed. Some of the proxies allowed us to assess paleo-facial conditions: dissolved oxygen concentration, depth fluctuations, distance to the source of terrigenous flux.  Argillaceous-siliceous organic-rich mudstone with planar microscale lamination (L-3): A distinctive feature is a pronounced thin lamination, detected only by microscope and coming from the layer-by-layer distribution of organic matter, clay minerals, and siliceous components. The average TOC content is 7-8 wt.%.

•
Argillaceous-siliceous mudstone with thin interlayering of organic-rich and organicpoor laminas and rare bioturbation (L-4): The rocks are characterized by rhythmic formationtal lamination due to the presence of dark gray (enriched in organic matter) and light gray layers, as well as bioturbation. The average TOC content is 2 wt.%.
For the investigated samples, the rock structure is related to organic matter particles shape, orientation and distribution, and its volume fraction is 1-10%. According to our results, L-3 is characterized by a strong dominance of thin laminas and lenticular inclusions of kerogen ( Figure 3C,G), whereas other lithological types host tiny dispersed or patchy organic matter, with a minor role of lenticular inclusions ( Figure 3). Variations of organic matter shape cause differences in petrophysical properties of rocks, such as mechanical anisotropy, wettability [37], and many others. Different structure and distribution patterns are determined by organic matter type and origin, which is further investigated in the paper.

Major and Trace Element Composition, Sm and Nd Isotopes
The elemental composition was determined for 131 samples taken at regular intervals along the section (the average distance between the samples is 15-25 cm) using portable XRF. In the resulting data set (Tables A1 and A2), significant variations in the elemental composition relative to the average values for the Abalak and Bazhenov Formations were identified. Therefore, we specified geochemical levels within the section.
Considering the characteristic paleo-geographic features of black shale formations [30,38], we paid special attention to the lithological composition and genesis of rocks: fluctuations in redox conditions, changes in the input of terrigenous and biogenic material, and bioproductivity of the paleo-sea.
According to the methods described in [39][40][41], we have calculated a set of lithochemical proxies, which were used to identify the genesis of rocks ( Figure 4). The most informative components to assess the contribution of terrigenous clastic components are Al, Ti, Zr, K, and the authigenic and biogenic components of the rock: S, P, Mo, Ca, Mg, Ni, Cu, Zn, V, etc. [41,42]. The carbonate content was studied via the ratios of Ca, Mg, Sr, Mn, and Fe. The predominant type of iron minerals (oxide, sulfide forms, the Fe affinity with S, Mn, and Ni) was revealed. Some of the proxies allowed us to assess paleo-facial conditions: dissolved oxygen concentration, depth fluctuations, distance to the source of terrigenous flux.
The total terrigenous contribution was estimated by the parameter (Al + K)/(Si + Ca + Mg), which was chosen specifically for siliceous black shales [30]. Unchanged ratios of chemical elements of terrigenous origin (Al, K) indicate the stable composition of the incoming terrigenous suspension. Mg is the main component of many clays, and Si can be of different origin: terrigenous (quartz grains, alumosilicates), biogenic (radiolarians), secondary (silicification). Variations in Ti/Zr, Ti/Al, Mg/Al, K/Al, and Rb/K are associated with the evolution of the terrigenous source or the contribution of several sources [40,41]. According to the increase in Zr (we consider the Ti/Zr ratio), one can reconstruct the distance to the terrigenous source, since Zr is the most stable component in the geochemical environment [40]. The absence of correlations between the contents of Si, Mg, and Fe and Al, Zr, Ti, and K, as well as a disproportionate excess of (Si + Mg + Fe) relative to (Al + Zr + Ti + K), indicates their authigenic origin and biogenic source.
The considered rocks contain a low amount of carbonate minerals, and the absolute carbonate content was determined for individual samples only (see Section 4.6). Relative carbonate content variations were evaluated by the ratio (Ca + Mg)/(Al + Si), where Ca and Mg are from the carbonate fraction, and Al and Si are attributed to the siliciclastic fraction and have terrigenous and mixed terrigenous-biogenic origin, respectively. Qualitative changes in the carbonate composition are caused by the source of carbonate matter (Sr/Ca [39]), depth and oxygen concentration (Mn/Ca), and partly by the total Mg content of the rocks (Mg/Al). Member #1 (3.5 m) is subdivided into three geochemical levels (1a, 1b, 1c), and the interval corresponds to the upper part of the Abalak Formation and is characterized by relatively high contents of terrigenous Fe, Ni, and Mn, low carbonate content, low contents of S and pyrite, and low concentrations of biophilic elements (Mo, V, and other transition metals), which refer to the presence of oxygen during sedimentation. Geochemical level 1a is characterized by the highest contents of non-biogenic Ni and Mn, by high K/Al, and by low Sr/Ba ratio. The sedimentation conditions of these rocks correspond to a shallow and well-aerated sea basin with a moderate terrigenous input, which is consistent with previously reported data [26,27]. For the above laying levels 1b and 1c, we observe slight fluctuations in sea level detected by Sr/(Ba+Sr), clay content (Al+K)/(Ca+Mg+Si), Ti/Zr, and a stepwise decrease in Ni concentrations as Ni/S and Ni/(Si+Ca+Mg) and Mn as Mn/Al and Mn/Ca. Geochemical member #2 (3 m) corresponds to the beginning of sedimentation of the Bazhenov Formation, which was established by the simultaneous prompt increase in the content of sulfur and sulfides, organic carbon, and Mo. Its lower level (2a) corresponds to a hydrologically active facies, which follows from the abrupt increase in Mn and P as well as fluctuations in clay content. The upper level (2b) has lower Mn, increased biogenic silica (Si/Al), and a decreased clay content, which can be associated with a slow rise in the sea level, as estimated from an increase in Sr/(Ba+Sr).
The bottom of geochemical member #3 (3 m) is the most distinct level in the section, characterized by the growth of biophilic S, Mo, Ni, Cu, Zn, changes in the concentrations of P, Ca, Sr, and Mn, with the lowest clay content within the section. The observed changes are in good agreement with the increase in U, S, content, total radioactivity, and TOC. Here, we interpret a rapid rise in sea level, which led to the weakening of the existing sources of terrigenous flux and subsequently to the change in the sedimentation to the biogenic regime. The transgressive period is marked near the boundary of levels 3a and 3b by the excess of P and Mn and minimum content of terrigenous components (Al, K, Ti, Zr), which all together indicates an active The contents of the most significant biophilic elements normalized to Al and some terrigenous elements were used as indicators of paleo-bioproductivity: P/Al, Si/Al, Mo/Al, S/(Rb + Zr), Zn/Al, Cu/Al, and Ni/Al. The amount of P and S in sedimentary rocks is predominantly controlled by biogenic material. Integration of the proxies is necessary, since the redox conditions were not stable during the accumulation of the studied deposits.
Water chemistry and element influx into a basin associated with terrigenous sediments can also impact the enrichment pattern of redox-sensitive trace elements. This controls the applicability of trace element-based redox proxies. In general, Mo, V, Ni, Cu, and Zn concentrate into the sediment under oxygen-depleted and anoxic environments together with organic matter [38,39,41,43,44]. The most sensitive indicators of the oxygen level for the rocks of the Bazhenov Formation are the ratios Mo/(Mn + Mo), V/Al, and Mn/Al [38][39][40][41]. Mn mainly exists in oxide mineral form, and therefore is not typical for sediments accumulated in anoxic conditions. Mn is one of the most sensitive indicators of oxygen concentration, and thus we consider the ratio Mn/Al.
The sulfur and sulfide content in the section is determined by several parameters: bioproductivity, anoxic conditions, and sedimentation rate. The sulphidity in the section is analyzed using the ratios of S and Fe, as well as S and Ni. In anoxic conditions, Fe and redox-sensitive metals exist mainly in the form of sulfide [38][39][40][41]. The total content of authigenic sulfur and its predominance over the terrigenous fraction was estimated by the S/(Rb + Zr) ratio, in which characteristic trace elements were used as a terrigenous indicator. Ni in rocks allows us to assess the authigenic biogenic contribution (by the direct relationship of Ni with S and other biophilic elements) or to reveal the introduced components (positive correlations of Ni with Fe, Mn, and terrigenous components).
One of the possible indicators of eustatic fluctuations with the relative stability of the drift sources is the ratio of Sr and Ba. The behavior of these elements in the marine geochemical environment is different, and the Sr content increases with increasing distance from the watershed (Sr/(Ba + Sr) ratio increases) [40,41], while Ba-based proxies should be used with caution in case of the alteration of redox conditions. During analysis of the elemental proxies, the frequent changes in rock composition along the interval were detected, and the section could be subdivided into 16 geochemical levels ( Figure 4). We combined the levels into eight geochemical members based on the data on gamma-spectral core profiling (U, Th, K), lithology, and TOC ( Figure 4, Table 1); members were considered to study rather general trends in paleoenvironment changes rather than the individual events. Member #1 (3.5 m) is subdivided into three geochemical levels (1a, 1b, 1c), and the interval corresponds to the upper part of the Abalak Formation and is characterized by relatively high contents of terrigenous Fe, Ni, and Mn, low carbonate content, low contents of S and pyrite, and low concentrations of biophilic elements (Mo, V, and other transition metals), which refer to the presence of oxygen during sedimentation. Geochemical level 1a is characterized by the highest contents of non-biogenic Ni and Mn, by high K/Al, and by low Sr/Ba ratio. The sedimentation conditions of these rocks correspond to a shallow and well-aerated sea basin with a moderate terrigenous input, which is consistent with previously reported data [26,27]. For the above laying levels 1b and 1c, we observe slight fluctuations in sea level detected by Sr/(Ba + Sr), clay content (Al + K)/(Ca + Mg + Si), Ti/Zr, and a stepwise decrease in Ni concentrations as Ni/S and Ni/(Si + Ca + Mg) and Mn as Mn/Al and Mn/Ca.
Geochemical member #2 (3 m) corresponds to the beginning of sedimentation of the Bazhenov Formation, which was established by the simultaneous prompt increase in the content of sulfur and sulfides, organic carbon, and Mo. Its lower level (2a) corresponds to a hydrologically active facies, which follows from the abrupt increase in Mn and P as well as fluctuations in clay content. The upper level (2b) has lower Mn, increased biogenic silica (Si/Al), and a decreased clay content, which can be associated with a slow rise in the sea level, as estimated from an increase in Sr/(Ba+Sr).
The bottom of geochemical member #3 (3 m) is the most distinct level in the section, characterized by the growth of biophilic S, Mo, Ni, Cu, Zn, changes in the concentrations of P, Ca, Sr, and Mn, with the lowest clay content within the section. The observed changes are in good agreement with the increase in U, S, content, total radioactivity, and TOC. Here, we interpret a rapid rise in sea level, which led to the weakening of the existing sources of terrigenous flux and subsequently to the change in the sedimentation to the biogenic regime. The transgressive period is marked near the boundary of levels 3a and 3b by the excess of P and Mn and minimum content of terrigenous components (Al, K, Ti, Zr), which all together indicates an active hydrological regime. Geochemical level 3b is distinguished by a short-term increase in clay content and a high content of redox-sensitive metals. Level 3c is characterized by the appearance of a stable anoxic environment with the accumulation of abnormally high organic matter content and the associated concentration of biophilic metals (Mo, V, Ni, U).
Rocks of member #4 (1.5 m) have an increased content of terrigenous material and decreased organic matter content, which coincides with lower gamma values and must result from moderate or periodical exposure to oxygen. The Ti/Zr, K/Al, and Si/Al ratios demonstrate a different composition of terrigenous material supplied to the area of sedimentation with respect to lower members #1 and #2. From this depth, biogenic silica accumulation is inhibited by the increasing volumes of clay. Member #4 is interpreted as an episode of the appearance of a new source of terrigenous material, the influence of which will gradually increase until member #8.
The composition of rocks from member #5 (2 m) reflects an episode of anoxia with intensive biogenic sedimentation, resulting in the accumulation of a large amount of TOC. Increased Mo/Al, Mo/(Mo + Mn), V/Al maximum and low Mn/Al ratios indicate that the anoxic environment expands from the bottom to the upper part of the member. The event of elevated bioproductivity is characterized by an increase in (Ca + Mg)/(Al + Si), P/Al, Cu/Al, Ni/(Si + Ca + Mg), Mo/Al, and Mo/(Mo+Mn), with a simultaneous decrease in terrigenous components (Ti/Zr, Mg/Al and others), and marks the boundary between two geochemical levels in member #5: the upper one (5b) corresponds to a more pronounced period of euxinia and organic matter accumulation, as it follows from TOC content, Mo and Ni proxies. The composition of terrigenous material is stable, according to qualitative proxies (K/Al, Rb/K, Si/Al, (Al + K))/(Ca + Mg + Si).
In member #6 (1 m), an intensive terrigenous inflow and reduced amounts of organic matter in the sediment are reconstructed by an increase in (Al + K)/(Ca + Mg + Si), Ti/Zr, and K/Al ratios. Such sedimentation conditions are associated with hydrodynamic activation of bottom currents and water aeration, which can be derived from increased Mn/Ca, Mn/Al ratios. The resulting accumulation of organic matter is reduced both by oxidative degradation of organic matter and dilution by significant amounts of clay material.
Member #7 (1 m) is the next stage of marine bioproductivity, with a reduced role of terrigenous sedimentation. An enriched in TOC member #7 was formed under a stable deep anoxic environment and anomalous enrichment in redox-sensitive Ni, V, and Mo was detected (the maxima for U, Mo/Al, Mo/(Mo + Mn), V/Al, Ni/(Si + Ca + Mg) and elevated Cu/Al are found). At the same time, terrigenous influx is moderate and its composition is stable, as it follows from (Al + K)/(Ca + Mg + Si), K/Al, and Ti/Zr. A thin carbonate layer is present in the upper part of the member, formed during the transition to oxic conditions (refer to Mn/Al), resulting in the drop in TOC accumulation at the boundary of members #7 and #8.
Member #8 (5 m) shows the most significant variations of composition and exhibits fluctuating enrichment in terrigenous components, with an occasional occurrence of biophilic components and organic matter in the rocks. It could be subdivided into three geochemical levels (8a, 8b, 8c). The depositional environments for 8a and 8c are close to that of member #6, with oxic conditions and intensive terrigenous influx, whereas for 8b, a suboxic environment is assumed. The rocks of 8b show increased clay content (Al + K)/ (Ca + Mg + Si), and Ti/Zr; in 8c, P, Ca, K/Al and Si/Al are increased. In general, member #8 is the less stable one in terms of sedimentation conditions, while the predominant regime is terrigenous sedimentation with low degree TOC accumulation.
Thus, lithological and geochemical data of the studied section indicate the homogeneity of the upper part of the Abalak deposits (member #1), transition to mixed terrigenous and marine sedimentation in #2 and #3 and alternation of the members with different genesis in the upper members #3-#8 ( Figure 4, Table 1). This alternation is caused by periodically increasing terrigenous influx and water oxygenation (members #4, #6, #8), and in between the sedimentation conditions were favorable for the concentration and preservation of organic matter (#3, #5, and #7) due to a low sedimentation rate, a high proportion of the biogenic component, and predominantly anoxic conditions. Several parameters indicate that the studied marine basin gradually expanded and deepened. An increase in the volume and depth of the basin started during member #2 accumulation led to the decreasing hydrodynamic activity of the bottom layers, which contributed to a gradual increase in the proportion of clay sedimentation (from member #4 to #8).
For 14 samples from the members accumulated under different sedimentation conditions, we used TIMS method to measure the content of Sm, Nd and isotope ratios 147 Sm/ 144 Nd, 143 Nd/ 144 Nd (Table 2, Figure 5). Variations in the isotopic composition of Nd along the section indicate the existence of at least two different terrigenous sources and their variable role in sedimentation. There are members characterized by mixed terrigenous and biogenic sedimentation #3, #4, #5, and #7, and their characteristic 143 Nd/ 144 Nd values do not exceed 0.512130 ( Figure 5). The increase in Nd isotope ratio for the members #1, #6, and #8 reflects an increase in terrigenous influx from the sources with higher values of 143 Nd/ 144 Nd: 0.512152, 0.512170, and 0.512149-0.512177, respectively.
According to [45,46], the material in the modern estuarine sediments of the Ob River (weathering crust of the Urals) has a lower isotopic ratio of Nd as compared to the Yenisei river deposits (weathering crust of the Siberian platform).   According to [45,46], the material in the modern estuarine sedim River (weathering crust of the Urals) has a lower isotopic ratio of Nd a the Yenisei river deposits (weathering crust of the Siberian platform). T values for the estuarine sediments of the Yenisei River are 0.512132-0.51 for the time of the Bazhenov Formation sedimentation), which corre measured values for members #1, #6, and #8 in the studied section (

Source Rock Characterization
We performed source rock characterization and TOC measure samples within the studied section. The results are shown in Figures 6 a A3.

Source Rock Characterization
We performed source rock characterization and TOC measurements for 150 samples within the studied section. The results are shown in Figures 6 and 7 and Table A3. According to the average TOC values for eight geochemical members (Tables 1  and A3), intervals with TOC < 5 wt.% can be considered as being depleted of organic matter to different extents (members #1, # 2, # 4, #6, and #8), and three members with TOC > 5 wt.% (#3, #5, and #7) are enriched in ТОС, although individual samples within each member can show higher or lower organic matter content (Figure 7).
The maturity of the organic matter in the studied section corresponds to the beginning of the oil window. The values of Tmax vary between 434 and 442 °C, with an average of 438 °C (Figure 6a).
The members with the increased role of biogenic sedimentation, member #3, #5, and #7, are enriched in TOC, and the rocks contain type II kerogen. They are potential oil source rocks. The hydrogen index HI varies from 249 to 511 mg HC/g TOC, and OI is below 6 mg CO2/g TOC with the only one exception (Figure 6b).
The organic matter from members #1, #2, #4, #6, #8 are mixed type II/III kerogen, with HI varying in a range from 44 to 300 mg HC/g TOC with a few exceptions, and moderate OI is from 5 to 30 mg CO2/g TOC for most of the samples. The significant admixture of humic organic matter with higher OI and lower HI is observed for the rocks in members #1, #2, #6, #8, whereas the individual samples from #4 contain type II kerogen.
These results are further confirmed by analysis of the composition and kinetics of kerogen thermal decomposition products.
The share of pyrolyzable carbon in TOC (GOC/TOC, %) is 10-30% for the samples from members #1, #2, #4, # 6, and #8 and 35-45% for the rocks from TOC-enriched #3, #5, and #7 with type II kerogen, with minor exceptions (Figure 7). Differences in GOC/TOC with the same thermal history and maturity within the studied 20 m originate from the composition of the initial organic matter and the differences in the kinetic characteristics of kerogen transformation to hydrocarbons.

Kerogen Thermal Decomposition Kinetics and Products
Analysis of kerogen decomposition products and bulk kinetics studies were used to obtain additional data on kerogen type, maturity, and petroleum generation characteristics for each member in accordance with the depositional environment.
The maturity of the organic matter in the studied section corresponds to the beginning of the oil window. The values of Tmax vary between 434 and 442 • C, with an average of 438 • C (Figure 6a).
The members with the increased role of biogenic sedimentation, member #3, #5, and #7, are enriched in TOC, and the rocks contain type II kerogen. They are potential oil source rocks. The hydrogen index HI varies from 249 to 511 mg HC/g TOC, and OI is below 6 mg CO 2 /g TOC with the only one exception (Figure 6b).
The organic matter from members #1, #2, #4, #6, #8 are mixed type II/III kerogen, with HI varying in a range from 44 to 300 mg HC/g TOC with a few exceptions, and moderate OI is from 5 to 30 mg CO 2 /g TOC for most of the samples. The significant admixture of humic organic matter with higher OI and lower HI is observed for the rocks in members #1, #2, #6, #8, whereas the individual samples from #4 contain type II kerogen.
These results are further confirmed by analysis of the composition and kinetics of kerogen thermal decomposition products.
The share of pyrolyzable carbon in TOC (GOC/TOC, %) is 10-30% for the samples from members #1, #2, #4, # 6, and #8 and 35-45% for the rocks from TOC-enriched #3, #5, and #7 with type II kerogen, with minor exceptions (Figure 7). Differences in GOC/TOC with the same thermal history and maturity within the studied 20 m originate from the composition of the initial organic matter and the differences in the kinetic characteristics of kerogen transformation to hydrocarbons.

Kerogen Thermal Decomposition Kinetics and Products
Analysis of kerogen decomposition products and bulk kinetics studies were used to obtain additional data on kerogen type, maturity, and petroleum generation characteristics for each member in accordance with the depositional environment.
Using flame ionization detector (FID) response in pyro-GC analysis, we calculated the distribution of different classes of compounds: gas and low-boiling liquid hydrocarbons (C1-C7), aliphatic (alkanes, cycloalkanes, and unsaturates), and aromatic components (mono-, di-, tricyclic) ( Table 3). Rock samples collected from intervals #1 and #8 contain gas-prone kerogen, which generates 36 and 44% percent of gas at elevated temperature, respectively. The aliphatic/aromatic compounds ratio is useful for oil hydrocarbons characterization since it correlates with oil density and viscosity. For intervals #1 and #8, the ratio is 17/83 and 11/89, respectively. Kerogen in rocks enriched with TOC (members #3, #5, and #7) is oil-prone, giving less than 10% of gas and oil with an aliphatic/aromatic ratio higher than 30/70. Intervals #2, #4, and #6 are in between two distinctive groups. There is 20% of gas in pyrolysis products for #4 and #6, and the aliphatic/aromatic ratio is around 15/85, whereas for #2, the composition is close to oil-prone intervals: the gas share is 16% and the aliphatic/aromatic ratio is 39/61. Therefore, kerogen in #2, #4, and #6 is classified as mixed oil-and gas-prone.

Thermal Extraction and Biomarker Analysis
Based on source rock characterization and analyses of kerogen decomposition products, nine samples were selected for biomarker analysis. These samples roughly represent common characteristics of organic matter and mobile hydrocarbon composition for each member, although the variations of composition within each interval are more significant. Since thermal extraction TD-GC×GC-MS (at 350 • C) is an express alternative to organic solvent extraction [47,48], we applied this approach for the analysis. The results are presented in Figures 8 and 9, and in Table 4.
The intervals enriched in TOC #3, #5, and #7 contain hydrocarbons similar to rock extracts and oils from the Bazhenov core in the central part of the West Siberian Basin, which were formed in deep-sea anoxic conditions from marine organic matter [9,49,50]. The main components of the extracts are petroleum hydrocarbons with a proportion of aliphatic components of about 80%, and that of aromatic components of no more than 10%. The maximum of the n-alkanes distribution corresponds to n-C 11 -C 13 (Figure 8). The extracts are characterized by molecular parameters and biomarkers of marine organic matter accumulated in suboxic conditions: Pr/Ph is in the range 2.3-2.6, the full distribution of homohopanes C31-C35, high C35 homohopane index, Ts/Tm = 0.9-1.5, C29/C30Hop = 0.5-0.7, the predominance of C27 in the distribution of regular steranes [20][21][22]51] and the samples are attributed to the area of mixed origin oils at the plot considered in [50] for the Bazhenov oils (Figure 9b). The regular steranes distribution is the only parameter that indicates the increasing terrigenous organic matter input in TOC of members #2, #3, #5, and #7: the share of C29 increases from #2 and #3 to #5 and #7 (19,19,21,26%, respectively) with the simultaneous decrease in C27.
The extracts from the depleted-in-TOC intervals #1, #2, #4, #6, and #8 differ from each other due to different depositional environments, and below we discuss the key features for each of the members.
Hydrocarbons in the rock from member #1 are condensate that contain alkanes, monoand di-aromatic components. The distribution of n-alkanes is narrow, with a pronounced maximum at n-C 15 . The extract contains no hopanes, steranes, or triaromatic steroids (TAS_), which is the evidence for oxic depositional environment and humic organic matter input [51].  The extractable organic matter in the sample from member #4 is represented by alkanes, mono-and di-aromatic components. The distribution of n-alkanes is broadened, and the maximum corresponds to n-C13-C16. The terpane ratios are Ts/Tm = 1, C29/C30Hop = 1.2. A distinctive feature of the extract is the presence of short C20 and C21 triaromatic steroids without higher homologs, TAS = 1. All parameters indicate that organic matter has a predominantly marine genesis, with a significant  The extractable organic matter in the sample from member #4 is represented by alkanes, mono-and di-aromatic components. The distribution of n-alkanes is broadened, and the maximum corresponds to n-C13-C16. The terpane ratios are Ts/Tm = 1, C29/C30Hop = 1.2. A distinctive feature of the extract is the presence of short C20 and C21 triaromatic steroids without higher homologs, TAS = 1. All parameters indicate that organic matter has a predominantly marine genesis, with a significant admixture of terrigenous as compared to members #3, #5, and #7.
Member #6 is an interval with easily recognizable lenses with different TOC. We performed two analyses of the extracts: for lighter and darker colored rocks. In the  The extract from member #2 combines features of both marine and terrigenous organics. n-alkanes distribution is wide, and there is an elevated concentration of long-chain alkanes from n-C 22 to C 30 . At the same time, the composition of sterane and hopane biomarkers is close to that of #3, which is the result of the transition from terrigenous in #1 to marine sedimentation that took place during member #3 accumulation.
For the studied deposits starting from member #3, accumulation of predominantly marine sapropelic organic matter in anoxic conditions periodically was interrupted by terrigenous material influx and water oxygenation which resulted in the dilution of organic matter by siliclastic material, the significant increase in sedimentation rate and lower degree of organic matter preservation.
The extractable organic matter in the sample from member #4 is represented by alkanes, mono-and di-aromatic components. The distribution of n-alkanes is broadened, and the maximum corresponds to n-C 13 -C 16 . The terpane ratios are Ts/Tm = 1, C29/C30Hop = 1.2. A distinctive feature of the extract is the presence of short C20 and C21 triaromatic steroids without higher homologs, TAS = 1. All parameters indicate that organic matter has a predominantly marine genesis, with a significant admixture of terrigenous as compared to members #3, #5, and #7.
Member #6 is an interval with easily recognizable lenses with different TOC. We performed two analyses of the extracts: for lighter and darker colored rocks. In the dark-colored sample, the distribution of n-alkanes is identical to that of the intervals #5 and #7 with high TOC (Table 4, Figure 8, see sample 6*). At the same time, in the lighter colored rock (sample 6), the distribution is wide, in a range from n-C12 to C20, without a pronounced maximum. There are no steranes and TAS. The terpane ratios are Ts/Tm = 2-3, C29/C30Hop = 1.0-1.5. Member #8 is a heterogeneous interval. For the selected sample, the ratio Ts/Tm = 0.8, C29/C30Hop = 1.1.
The extract composition for the members with TOC < 5 wt.% is determined by the admixture of terrigenous organic matter and oxidative accumulation conditions. Regular steranes composition indicates the increasing influx of terrigenous organic during predominant biogenic sedimentation: the admixture increases from #2 and #3 to #5 and #7.

Elemental and Isotope Composition of Carbon, Nitrogen, and Sulfur
The isotopic composition and content of organic carbon, nitrogen, sulfur and carbonate carbon, and oxygen were determined in 18 selected samples from all members (Figures 10 and 11 and Table A4). The carbonate carbon content is generally low in the studied sedimentary sequence, and Ccarb is less than 1 wt.% for most of the samples. Only three layers in investigated section demonstrate higher carbonate contents: one layer in the lower part (#1, sideritized argillaceous-siliceous rock) and two in the upper part (#7 and #8, calcified radiolarite), with Ccarb 8.1, 5.1, and 9.9 wt.%, respectively (Table A4). Samples from the middle part of the section (members #2-7) are weakly carbonated siliceous- The carbonate carbon content is generally low in the studied sedimentary sequence, and Ccarb is less than 1 wt.% for most of the samples. Only three layers in investigated section demonstrate higher carbonate contents: one layer in the lower part (#1, sideritized argillaceous-siliceous rock) and two in the upper part (#7 and #8, calcified radiolarite), with Ccarb 8.1, 5.1, and 9.9 wt.%, respectively (Table A4). Samples from the middle part of the section (members #2-7) are weakly carbonated siliceous- Figure 11. (a) Cross-plot of the elemental and isotopic composition of nitrogen. (b) Cross-plot of the elemental and isotopic composition of sulfur. The two series of samples represent rocks with C org < 5 wt.% (blue markers) and C org > 5 wt.% (red markers), numbers correspond to geochemical members.
The organic carbon concentrations measured using elemental analyzer (C org ) are consistent with the TOC values obtained by pyrolysis and vary from 0.5 to 11.1 wt.%. The isotopic composition of organic carbon takes the range of δ 13 C org from −32.4 to −28.5‰ (PDB) (Figure 10a). Samples with an organic carbon content higher than 5% are depleted in 13 C with δ 13 C org from −31.3 to −32.4‰, as compared to samples with organic carbon content below 5 wt.% with δ 13 C org in the range from −31.1 to −28.5‰.
The difference in isotope composition correlates with the difference in the environment of organic matter accumulation and its origin, which results in high and low TOC. According to [2,52], for the Bazhenov Formation, the δ 13 C org below −30‰ is characteristic for the marine polymer-lipid type of organic matter (#3, #5, and #7), while higher and more positive values reflect the input of organic matter of terrigenous origin, as in #1, #2, #4, #6, and #8.
The concentrations of nitrogen in the investigated samples vary from 0.05 to 0.38 wt.% and demonstrate a strong positive correlation with organic carbon concentrations, indicating the organic source of nitrogen (Figure 10b). The δ 15 N is in the range from 2.5 to 9.8‰ (AIR) (Figure 11a), and samples with C org > 5% are depleted in 15 N, with δ 15 N values from 2.5 to 5.9‰; meanwhile, samples with low C org are generally enriched in 15 N, and the δ 15 N value range is from 3.8 to 9.8‰ and tend to increase with decreasing organic carbon concentration. Observed C/N ratios are affected by both organic matter source and diagenetic alterations, and for high-TOC members, C/N varies in a narrow range, 24-31, whereas samples from low-TOC members show lower ratios, [8][9][10][11][12][13][14][15][16][17][18][19][20][21][22], that decrease with decreasing C org . The only exception is the sample with C/N 32, which has C org 0.6 wt.% and C carb 9.9 wt.%. These variations also speak in favor of the presence of organic matter of different origin in layers with high and low TOC.
Total sulfur content in the investigated rocks takes a wide range of values from 0.2 to 8.2 wt.%, exhibits a positive correlation with iron content, and no correlation with organic carbon concentration (Figure 11b). The main chemical form of sulfur is sulfide sulfur fixed as pyrite. Organic sulfur is present in a subordinate amount, since its content in organic matter in sedimentary rocks of the Bazhenov Formation is generally low and does not exceed 5 wt.% [2]. In contrast to sulfur concentration, the isotopic composition of sulfur demonstrates correlation with organic carbon concentration, the δ 34 S values from −33.7 to −7.8‰ are significantly higher for the samples with TOC < 5 wt.% than for the samples enriched in TOC with the δ 34 S from −37.3 to −25.8‰ (Figure 11b). The isotopic composition of sulfide (S 2− ) in sediments is controlled by the degree of bacterial reduction of sulfate (SO4 2− ) dissolved in water since the produced H 2 S is then used for pyrite formation and organic matter sulfurization. The process is affected by the redox environment in the bottom water layer and sediment column [53]. Under suboxic to euxinic conditions, sulfate reduction occurs both in bottom water and upper sediment layers (beneath the water-sediment interface). The system is open with respect to sulfate, and even at high sulfate reduction rates, the hydrogen sulfide and formed pyrite express very negative δ 34 S values due to the kinetic isotope effect between SO4 2− and H 2 S, which can reach −66‰ [54]. Under oxic conditions, there is no sulfate reduction in water and upper sediment layers. It only takes place in pore water of deeper buried sediments. Thus, the system is closed with respect to sulfate, and in the case of sulfate exhaustion, the hydrogen sulfide and pyrite can reach positive δ 34 S values. According to these mechanisms, low δ 34 S values and high organic carbon content in the investigated samples reflect suboxic to euxinic conditions in the paleo-sea, low organic carbon content (#3, #4, #5, #7, and #8) and higher δ 34 S values corresponding to oxic conditions during sedimentation (#1, #2, and #6).
The carbonate carbon content is generally low in the studied sedimentary sequence, and C carb is less than 1 wt.% for most of the samples. Only three layers in investigated section demonstrate higher carbonate contents: one layer in the lower part (#1, sideritized argillaceous-siliceous rock) and two in the upper part (#7 and #8, calcified radiolarite), with C carb 8.1, 5.1, and 9.9 wt.%, respectively (Table A4). Samples from the middle part of the section (members #2-7) are weakly carbonated siliceous-argillaceous, argillaceous-siliceous and argillaceous rocks.
The carbonate carbon and oxygen isotopic composition vary considerably within the studied section: δ 18 O carb varies from −14.6 to −0.3‰VPDB and δ 13 C carb from −26 to −0.4‰VPDB. The three abovementioned samples with carbonate content more than 5 wt.% demonstrate the highest δ 18 O values within the section, close to 0‰, which is typical for sedimentary carbonates. Other samples with C carb not exceeding 1 wt.% are essentially depleted in 18 O, and the δ 13 C carb values of these samples vary from −0.4 to −16.6‰. We suggest that the original isotopic composition, characteristic for sea water carbonates, was changed in carbonates of the investigated section in diagenesis and catagenesis, when the secondary carbonatization could take place at elevated temperatures in the presence of isotopically light CO 2 derived from the organic matter [55]. The layers with low carbonate content are mostly affected by secondary processes. One exception is the sample from member #8 with C carb 9.9 wt.%, which, due to the low δ 13 C value, −26‰, and high δ 18 O value, −0.3‰, corresponds to sedimentary-diagenetic methane-derived limestones [56,57]. This type of carbonate with isotopically light carbon could be formed on the seabed in the areas of focused methane discharge under anaerobic conditions when methane, entering the bottom water, undergoes oxidation by microorganisms.

Integration of the Geochemical Data and Reconstruction of Sedimentation Conditions
Two genetically distinct groups of rocks were identified according to their geochemical features, they correspond to two sedimentation regimes. The first regime is characterized by noticeable input of biogenic material and pronounced marine organic matter accumulation in sediments, and another regime is associated with terrigenous sedimentation and mixed marine and terrigenous organic matter in rocks (Table 1, Figure A3).
Episodes of the increased role of biogenic sedimentation and organic matter accumulation occurred during the sea level transgression, the beginning of which was detected at the boundary of members #2 and #3 by the minimum content of terrigenous particles and the highest contribution of biogenic Si, S, P, Mg, and redox-sensitive metals. The peaks of biogenic sedimentation occur in members #3, #5, and #7. Their presence in the section in most cases correlates with episodes of anoxic environments characterized by Mo/(Mo+Mn) and V/Al, U, Th/U, Mo, and Mo/TOC [43,[58][59][60][61][62][63]. The intervals are characterized by high TOC content, since a weak input of terrigenous material leads to the concentration of organic matter without its dilution.
The second series of rocks contain a significantly higher amount of the silt-clay component and could be distinguished by higher ratios of terrigenous elements Ti and Zi, Ti and Si, and Ti and K ( Figure A3). Such rocks were formed during the enrichment of sediments with terrigenous material. Fluctuations in the terrigenous influx have determined the rate of organic matter accumulation and composition. In this section, we can assume the influence of more than two terrigenous sources. Starting from member #4, a new terrigenous source appears (different from that influenced members #1 and #2), and its influence is detected within the section up to member #8 by Ti/Zr and Rb/K, and 143 Nd/ 144 Nd. Nd isotope ratio proves that the mineral part of members #1, #6, and #8 was significantly impacted by mineral component close in composition to the weathering crust of the Siberian platform.
Integrating all the factors, the studied section corresponds to an area with eustatic sea level changes during transgressive/regressive cycles and periodic activation of currents, supplying the sediments with terrigenous material and improving aeration of the bottom layers of the water column. The conclusion is consistent with the reported data on the paleogeography of the area in the Tithonian [3,64] (Figure 12). It was located in the isthmus between the Boreal Sea and the epicontinental Bazhenov Sea surrounded by Laurasia (modern Polar Urals and the Siberian Craton). Geosciences 2021, 11, x FOR PEER REVIEW 22 Figure 12. Paleogeography of the studied area (dark rhomb) in the Tithonian modified after Kontorovich [3].
Under the influence of the above discussed factors, including terrestrial influx redox conditions, the amount and composition of organic matter in the members va significantly (Figures 13 and 14).
In members with mixed terrigenous and biogenic sedimentation (#3, #5, and the accumulated TOC is higher than 5 wt.% in the vast majority of the samples. rocks contain marine organic matter, represented by type II oil-prone kerogen. hydrocarbon fluid characteristics are close to that of the sediments of the Bazhe Formation in the central part of the Basin in terms of molecular, biomarker and isot composition (Figures 9b, 13 and 14) [9,49,50].
Organic matter in member #4 with average TOC = 4.3 wt.% is predominantl marine origin, accumulated in more oxic conditions with a noticeable inpu terrigenous material. Its transformation leads to the formation of oil and condensate, but the amount is lower compared to members #3, #5, and #7 ( Figure   (a) (b) Figure 12. Paleogeography of the studied area (dark rhomb) in the Tithonian modified after Kontorovich [3].
Under the influence of the above discussed factors, including terrestrial influx and redox conditions, the amount and composition of organic matter in the members varies significantly (Figures 13 and 14). Under the influence of the above discussed factors, including terrestrial influx and redox conditions, the amount and composition of organic matter in the members varies significantly (Figures 13 and 14).
In members with mixed terrigenous and biogenic sedimentation (#3, #5, and #7) the accumulated TOC is higher than 5 wt.% in the vast majority of the samples. The rocks contain marine organic matter, represented by type II oil-prone kerogen. The hydrocarbon fluid characteristics are close to that of the sediments of the Bazhenov Formation in the central part of the Basin in terms of molecular, biomarker and isotopic composition (Figures 9b, 13 and 14) [9,49,50].
Organic matter in member #4 with average TOC = 4.3 wt.% is predominantly of marine origin, accumulated in more oxic conditions with a noticeable input of terrigenous material. Its transformation leads to the formation of oil and gas condensate, but the amount is lower compared to members #3, #5, and #7 ( Figure 14). Figure 13. Relations between the Th/U ratio (logarithmic scale) and the isotope composition of organic matter (a) and Nd (b). The two series of samples represent TOC < 5 wt.% (blue) and TOC > 5 wt.% (red markers), and numbers correspond to geochemical members.

(a) (b)
Geosciences 2021, 11, x FOR PEER REVIEW 23 of 31 Figure 13. Relations between the Th/U ratio (logarithmic scale) and the isotope composition of organic matter (a) and Nd (b). The two series of samples represent TOC < 5 wt.% (blue) and TOC >5 wt.% (red markers), and numbers correspond to geochemical members. The other members with average TOC < 5 wt.%, #1, #6, #8, were accumulated during periods of increased terrigenous influx and under oxic conditions. The terrigenous organic matter admixture is unambiguously identified by heavier Corg isotopic composition, only trace amounts of steranes, triaromatic steroids, and homohopanes, and C29/C30hop > 1 (in the case of interval #1, hopanes and steranes are absent). TOC in most samples does not exceed 4 wt.%, organic matter shows low GOC/TOC and HI, and high oxygen index (OI) = 19-26 mg CO2/g TOC. Kerogen thermal decomposition products are mainly gas and heavy polyaromatic components.
For member #6, rocks are characterized by higher PI and extract composition, which is similar to the extracts of nearby kerogen-saturated rocks #5 and #7. These observations indicate better reservoir properties of the rocks from interval #6. The improved accumulating properties of the rock and the maturity of organic matter at the beginning of oil generation suggest the prospects for petroleum development. For other members with low TOC and a moderate to high proportion of terrigenous components (members #1, #2, #4, and #8), no reservoir properties were observed.
The main arguments in favor of the promising development of the studied area are high content of organic matter in the rock, a sufficient degree of maturity for kerogen conversion into mobile hydrocarbons, the presence of syngenetic mobile hydrocarbons in the section, and the alternation of intervals enriched with kerogen (and thus showing generating properties) with potentially accumulating intervals with increased reservoir properties.

Conclusions
In the Upper Jurassic-Lower Cretaceous section in the north of the West Siberian Petroleum Basin, we have identified eight geochemical members, which differ according to the elemental chemistry, amount and composition of organic matter. Using the geochemical and isotope data, we have performed a reconstruction of the In members with mixed terrigenous and biogenic sedimentation (#3, #5, and #7) the accumulated TOC is higher than 5 wt.% in the vast majority of the samples. The rocks contain marine organic matter, represented by type II oil-prone kerogen. The hydrocarbon fluid characteristics are close to that of the sediments of the Bazhenov Formation in the central part of the Basin in terms of molecular, biomarker and isotopic composition (Figures 9b, 13 and 14) [9,49,50].
Organic matter in member #4 with average TOC = 4.3 wt.% is predominantly of marine origin, accumulated in more oxic conditions with a noticeable input of terrigenous material. Its transformation leads to the formation of oil and gas condensate, but the amount is lower compared to members #3, #5, and #7 ( Figure 14).
The other members with average TOC < 5 wt.%, #1, #6, #8, were accumulated during periods of increased terrigenous influx and under oxic conditions. The terrigenous organic matter admixture is unambiguously identified by heavier C org isotopic composition, only trace amounts of steranes, triaromatic steroids, and homohopanes, and C29/C30hop > 1 (in the case of interval #1, hopanes and steranes are absent). TOC in most samples does not exceed 4 wt.%, organic matter shows low GOC/TOC and HI, and high oxygen index (OI) = 19-26 mg CO 2 /g TOC. Kerogen thermal decomposition products are mainly gas and heavy polyaromatic components.
For member #6, rocks are characterized by higher PI and extract composition, which is similar to the extracts of nearby kerogen-saturated rocks #5 and #7. These observations indicate better reservoir properties of the rocks from interval #6. The improved accumulating properties of the rock and the maturity of organic matter at the beginning of oil generation suggest the prospects for petroleum development. For other members with low TOC and a moderate to high proportion of terrigenous components (members #1, #2, #4, and #8), no reservoir properties were observed.
The main arguments in favor of the promising development of the studied area are high content of organic matter in the rock, a sufficient degree of maturity for kerogen conversion into mobile hydrocarbons, the presence of syngenetic mobile hydrocarbons in the section, and the alternation of intervals enriched with kerogen (and thus showing generating properties) with potentially accumulating intervals with increased reservoir properties.

Conclusions
In the Upper Jurassic-Lower Cretaceous section in the north of the West Siberian Petroleum Basin, we have identified eight geochemical members, which differ according to the elemental chemistry, amount and composition of organic matter. Using the geochemical and isotope data, we have performed a reconstruction of the depositional environment and evaluated the factors that have affected the distribution and composition of organic matter in the source rocks.
It was shown that the TOC content in the rocks is mainly controlled by the redox conditions and the intensity of terrigenous influx during sedimentation.
We attribute a simultaneous change in redox conditions and the intensity of terrigenous influx in the studied area to eustatic sea level changes during transgressive/regressive cycles and activation of currents. Transgression periods are associated with a weakened influence of terrigenous influx, low oxygenation with periodically expanded anoxia and intense accumulation and preservation of marine organic matter with oil-generating characteristics. The members #3, #5, and #7 were developed under such a depositional environment.
During eustatic regression periods, an increase in terrigenous influx and activation of currents in the studied area occurred simultaneously with an increase in oxygen, dissolved in water, and the proportion of terrigenous organic matter. In the lower part of the considered section, members #1 and #2 reflect these conditions. Starting from members #2 and #3, the periods of an intense accumulation of marine organic matter under reducing conditions were periodically changed to oxic environments with a more abundant terrigenous influx from another weathering source, which resulted in the fluctuations in organic matter content and composition. The generally moderate terrigenous inflow along the section is found increased for member #4, and to a greater extent for members #6 and #8. Under oxic conditions with the increased proportion of the terrigenous material, the formation of rocks depleted in organic matter, containing gas-generating kerogen, occurred (members #1, #2, #6, and #8).
The obtained results allowed us to evaluate the relationships between the sedimentation conditions and the petroleum generation characteristics of the current deposits. The paleoenvironment patterns, identified using a set of geochemical and isotopic methods, form the basis for the studies of different marine basins.
Using a set of geological, geochemical and isotopic methods, we have evaluated the relationships between paleo-sedimentation conditions, the composition of the source rocks and the petroleum generation characteristics of organic matter for the section of Bazhenov Formation in the north of the West Siberian Petroleum Basin. The identified paleoenvironment patterns contribute to the understanding of the Bazhenov unconventional resource formation and the prediction of hydrocarbon source and reservoir character in the time equivalent deposits in the region.   Activation energy distribution for low TOC intervals (<5%) Figure A1. Representative bulk kinetic spectra (A = 1·10 14 s −1 ) of kerogen decomposition for low TOC members #1, #2, #4, #6, #8.