Understanding how alluvial stratigraphy responds to sediment supply perturbations is critical for interpreting past environmental changes from the sedimentary record, characterizing subsurface reservoirs, and forecasting future landscape evolution. However, identifying and quantifying sediment supply signals preserved in the rock record remain challenging, leaving their stratigraphic imprint insufficiently understood. To help address this issue, we use a process-based numerical model to simulate alluvial stratigraphy under different sediment supply scenarios, independently testing the effects of supply magnitude and variability. Our results show that sediment supply variability has a stronger impact than magnitude: increased variability leads to much thicker channel-belt deposits and elevated yet alternating high and low down-valley slopes. In contrast, greater total sediment supply results in only slightly thicker channel-belt deposits and uniformly elevated down-valley slopes. These results reconcile diverse fluvial stratigraphic responses to sediment-supply changes across basins during climatic perturbations.

Alluvial stratigraphy builds up over geologic time under the complex interplay of external climatic and tectonic forces and internal stochastic processes. This complexity makes it challenging to attribute alluvial stratigraphic changes to specific factors. Geological records indicate pronounced and persistent climatic changes during the Phanerozoic, while the effects of these changes on alluvial stratigraphy remain insufficiently documented. We provide evidence for 405 k.y. long-eccentricity climate forcing of alluvial stratigraphy in the lower Eocene Willwood Formation of the Bighorn Basin, Wyoming (USA). Two ∼90-m-thick intervals, characterized by a relative paucity of sand, dominance of sinuous-river channels, and floodplain sediments with better-developed paleosols, coincide with eccentricity maxima as determined through integrated stratigraphic methods. These intervals are interspersed with three contrasting intervals, marked by relatively high sand content, prevalent braided-river channels, and less-developed paleosols, corresponding to eccentricity minima. A comprehensive genetic model that integrates climate, source-to-sink system, and alluvial dynamics to explain these findings remains to be elucidated. Given the consistent presence of the 405 k.y. eccentricity cycle throughout Earth’s history, it is plausible to infer that its influence may be discernible across a wide array of alluvial stratigraphic records.

Sedimentation on river floodplains is a complex process that involves overbank flooding, crevasse splaying, and river avulsion. The resulting floodplain stratigraphy often exhibits floodplain aggradation cycles with alternating fine-grained overbank flooding deposits that underwent significant petrogenesis, and coarser-grained, avulsion-belt deposits largely devoid of pedogenic impact. These cycles are linked to lateral migration and avulsion of channels driven by internal dynamics, external factors, or a combination of both. To better understand the spatial and vertical variability of such floodplain aggradation cycles, we map these in three dimensions using a photogrammetric model of the lower Eocene Willwood Formation in the northern Bighorn Basin, Wyoming, USA. This allows identifying 44 floodplain aggradation cycles in ∼300 m of strata with an average thickness of 6.8 m and a standard deviation of 2.0 m. All the cycles are traceable over the entire model, pointing to their spatial consistency over the 10 km2 study area. At the same time, rapid lateral thickness changes of the floodplain aggradation cycles occur with changes up to 4 m over a lateral distance of 400 m. Variogram analyses of both field and numerical-model results reveal stronger consistency of floodplain aggradation cycle thicknesses along the paleoflow direction compared to perpendicular to paleoflow. Strong compensational stacking occurs at the vertical scale of 2–3 floodplain aggradation cycles (14–20 m), while full compensational stacking occurs at larger scales of more than six floodplain aggradation cycles (>41 m). The lateral and vertical thickness variability of the floodplain aggradation cycles, as well as their compensational stacking behavior, are interpreted to be dominantly driven by autogenic processes such as crevasse splaying and avulsing that preferentially fill topographic lows. External climate forcing may have interacted with these autogenic processes, producing the laterally persistent and vertically repetitive floodplain aggradation cycles. The spatial variability of floodplain aggradation cycles demonstrated in this study highlights again the need for three-dimensional data collection in alluvial floodplain settings rather than depending on one-dimensional records.

The middle Eocene Dongying sag in the Bohai Bay Basin of China has an estimated shale oil resource of approximately 1.1 billion t (8.06 billion bbl); flows of shale oil have been produced in the succession from tens of wells, where the daily production of a single well generally varies between 10 and 100 t (73.3–733 bbl). Therein, the mudrock successions composed of meter-scale mudstone–limestone couplets are the most important shale oil-producing layers. The controls on the deposition of the meter-scale mudstone–limestone couplets, however, remain enigmatic, constraining the analysis of lithofacies and, therefore, sweet spot distributions. Here, we analyze three continuously cored organic-rich successions of mudstone–limestone couplets (371 m [1217 ft] in total) in the middle Eocene Dongying sag, accompanied by decimeter- to meter-scale sampling and testing of mineralogy, organic geochemistry, and paleontology of the rocks. Our integrated cyclostratigraphic analysis shows that the observed mudstone–limestone couplets occur at periods that coincide with Milankovitch periodicities; 21-k.y. precession cycles are the main driver of the meter-scale mudstone–limestone couplets, with additional imprints of 41-k.y. obliquity cycles. Specifically, precession minima are associated with high summer insolation and consequently high summer monsoonal precipitation, which increased river discharge and terrigenous input to the basin, resulting in the deposition of siliciclastic-rich mudstones. In the study, low summer insolation during precession maxima led to decreased summer monsoonal precipitation, lower river discharge and terrigenous input, and increased lake water salinity, resulting in the deposition of authigenic lime mudstones. The shale reservoir quality kept pace with the orbital climate changes; compared with lime mudstones deposited during precession maxima, mudstones deposited during precession minima had higher total organic carbon, porosity, and oil content, but lower brittleness.

Formation of alluvial stratigraphy is controlled by autogenic processes that mix their imprints with allogenic forcing. In some alluvial successions, sedimentary cycles have been linked to astronomically-driven, cyclic climate changes. However, it remains challenging to define how such cyclic allogenic forcing leads to sedimentary cycles when it continuously occurs in concert with autogenic forcing. Accordingly, we evaluate the impact of cyclic and non-cyclic upstream forcing on alluvial stratigraphy through a process-based alluvial architecture model, the Karssenberg and Bridge (2008) model (KB08). The KB08 model depicts diffusion-based sediment transport, erosion and deposition within a network of channel belts and associated floodplains, with river avulsion dependent on lateral floodplain gradient, flood magnitude and frequency, and stochastic components. We find cyclic alluvial stratigraphic patterns to occur when there is cyclicity in the ratio of sediment supply over water discharge (Q_s/Q_w ratio), in the precondition that the allogenic forcing has sufficiently large amplitudes and long, but not very long, wavelengths, depending on inherent properties of the modelled basin (e.g. basin subsidence, size, and slope). Each alluvial stratigraphic cycle consists of two phases: an aggradation phase characterized by rapid sedimentation due to frequent channel shifting and a non-deposition phase characterized by channel belt stability and, depending on Q_s/Q_w amplitudes, incision. Larger Q_s/Q_w ratio amplitudes contribute to weaker downstream signal shredding by stochastic components in the model. Floodplain topographic differences are found to be compensated by autogenic dynamics at certain compensational timescales in fully autogenic runs, while the presence of allogenic forcing clearly impacts the compensational stacking patterns.

Fossil fuel resources are invaluable to economic growth and social development. Understanding the formation and distribution of fossil fuel resources is critical for the search and exploration of them. Until now, the vertical distribution depth of fossil fuel resources has not been confirmed due to different understandings of their origins and the substantial variation in reservoir depths from basin to basin. Geological and geochemical data of 13 634 source rock samples from 1286 exploration wells in six representative petroliferous basins were examined to identify the maximum burial depth of active source rocks in each basin, which is referred to in this study as the active source rock depth limit (ASDL). Beyond the ASDL, source rocks no longer generate or expel hydrocarbons and become inactive. Therefore, the ASDL also sets the maximum depth for fossil fuel resources. The ASDLs of basins around the world are found to range from 3000 to 16 000 m, while the thermal maturities (Ro) of source rocks at the ASDLs are almost the same, with Ro ≈ 3:5±0:5 %. The Ro of 3.5% can be regarded as a general criterion to identify ASDLs. High heat flow and more oil-prone kerogen are associated with shallow ASDLs. In addition, tectonic uplift of source rocks can significantly affect ASDLs; 21.6 billion tons of reserves in six representative basins in China and 52 926 documented oil and gas reservoirs in 1186 basins around the world are all located above ASDLs, demonstrating the universal presence of ASDLs in petroliferous basins and their control on the vertical distribution of fossil fuel resources. The data used in this study are deposited in the repository of the PANGAEA database at: https://doi.org/10.1594/PANGAEA.900865 (Pang et al., 2019).

The middle Eocene was a key period of global climate change from a warm “greenhouse” to a cooling “doubthouse”. At the middle Eocene, an extensive arid climate in China was recorded in the mudstones and shales interbedded with salinastone which formed good hydrocarbon source rocks, especially in the Bohai Bay Basin. The sedimentary paleoenvironment and its variations affected the source rock quality and distribution. Based on successive and dense sampling of the middle Eocene Shahejie Formation shale (MES shale) in Nanpu Sag of Bohai Bay Basin, petro-mineralogical tests, geochemical methods and spectrum analysis were performed to analyse the palaeoenvironmental fluctuations and the responses to Milankovitch cycles. Element analysis on MES shale shows that the middle Eocene climate was arid and the sediments were deposited under dysoxic/suboxic conditions with relatively high salinity. Based on the Δlog R method, the organic matter (OM) abundance was evaluated and ranged from 0.18 to 3.55 wt % for MES shale. The mean primary productivity of 545 g C m⁻² yr⁻¹ shows a eutrophic environment in middle Eocene Nanpu paleolake. Wavelet analysis performed on gamma logging and palaeoenvironmental proxies revealed clear Milankovitch cycles for MES shale, and the cyclostratigraphic and palaeoenvironmental fluctuations were effected by astronomical oscillation periods of the short eccentricity cycle (100 k.y.), obliquity cycle (39 k.y.) and precession cycle (19 k.y.). The depositional rate of MES shale was corrected to 0.43 m/k.y. The sedimental process of MES shale was associated with palaeoenvironmental variation, which have been divided into four stages affecting OM accumulation. A hypothesis that the vast carbon of “greenhouse” gases in the Palaeocene and early Eocene was fixed in organic matters during middle Eocene and enriched in the hydrocarbon source rocks of the sedimentary basins formed at early Eocene was proposed to explain the “doubthouse” world in the Eocene.

In order to determine the source of the light hydrocarbons in the Kuqa Depression of the Tarim Basin, Northwest China, their geochemical characteristics are documented, and possible influencing factors are discussed in a broad context on the basis of geochemical analyses of natural gas samples from 19 wells in combination with previously reported natural gas components and carbon isotope values from the same depression. The biomarker analytical results show that (a) methyl cyclohexane (37.5–60.3%) of natural gas predominates the C₇ light hydrocarbon series (methyl cyclohexane, heptanes, and dimethyl cyclopentane), and (b) the heptane content (25.0–54.8%) is also relatively high. The C_6–7 light hydrocarbon series (aromatics, normal alkanes, and cycloalkanes) have abnormally high aromatics (benzene and toluene) content (26.2–83.8%). Heptane value (13.31–37.77%) and isoheptane value (2.10–7.64%) indicate that the Kuqa natural gases are of high maturity and typical of coal-derived gas (gas generated by sapropelic source rock). Their light hydrocarbons, however, appear to have some characteristics of mixed-source organic matter. Specifically, the light hydrocarbons in the Dina 2 gas field show characteristics of the coal-derived gas; those in the Dabei gas field show some features of the oil-associated gas (gas generated by humic source rock), whereas those in the Kela 2 and Keshen 2 gas fields and Well Yangtake-101 display characteristics of the mixed source organic matter. The abnormally high content of the aromatics in the Kuqa Depression may have resulted mainly from their high productivity during the thermal maturity of humic type organic matters and the gases accumulation of late stage. The different relative contents of the aromatics in various oil and gas fields may have been accounted by different organic matter types, maturity, and distribution patterns of source rocks.

The effects of pore throat structure on gas permeability in tight sand reservoir were investigated using helium-measured porosity, pulse decay permeability, casting thin sections, scanning electron microscopy (SEM), nuclear magnetic resonance (NMR), and constant-rate mercury injection for five sandstone samples from the Upper Triassic Yanchang Formation in the Western Ordos Basin, China. Results show that permeability values are in the range of 0.474-2.290 mD. The strong variation may be due to the effect of mineral composition and pore-throat structure. Generally, sandstone samples of higher content of quartz, chlorite film, and slit-shaped throat tend to have higher permeability. The movable fluid porosity numerically equals to the product of porosity and the movable fluid saturation based on the film bound water, which can comprehensively reflect the influence of fluid flow channels on porosity while neglecting the influence of bound fluids in the fine pores of clay minerals in the tight reservoirs. Larger movable fluid porosity values mean larger pore spaces for the movable fluids, which contribute to larger permeability. It is more reasonable to analyze the movable fluid pore-throat radius as a range. The homogeneous pore throat structure lead to a narrow lower limit range of the movable fluid pore-throat radius, which is favorable for the fluid flow. Meanwhile, if the lower bound of the lower limit range of the movable fluid pore-throat radius is small, pores and throats are less plugged, and thus make the conditions favorable for fluid flow. The contribution of throats to the pore volume is critical in influencing reservoir permeability as it can directly control the lower limit of the movable fluid pore-throat size. The pore-throat radius where throats account for the smallest amount in the pore space increases with the permeability, which indicates that larger proportions of large throats result in larger throat radius and higher pore connectivity, leading to higher permeability.