Sandstones from the Main Buntsandstein Subgroup represent a promising deep geothermal target in the subsurface of the Netherlands considering their widespread distribution and temperatures locally reaching 140-150 °C at depths of ~ 3 to 5 km. The Main Buntsandstein Subgroup is a sand-prone interval, but the reservoir quality of these sandstones is known to be heterogeneous as result of an interplay between depositional and diagenetic processes. This makes the Buntsandstein sediments an uncertain and risky geothermal play. In this project, we assess the syn- and post-depositional history of these sediments. The aim is to define structural, sedimentary, and diagenetic heterogeneities within the Main Buntsandstein sediments and assess their impact on reservoir quality. This will help reduce uncertainties for geothermal operations in the Triassic in the southern Netherlands and beyond. The structural analysis of the study area using seismic and well data reveals that the Main Buntsandstein sediments represent an early syn-rift sequence and that their present-day distribution is strongly controlled by faulting. In parallel, the study of the sedimentology and stratigraphy conducted on core and wireline data indicates that the depositional environment evolves through the Buntsandstein stratigraphy, resulting in the development of different reservoir architectures. Diagenesis has largely altered the primary relationship between sedimentary facies and porosity and permeability. Overall cementation seems to have a larger impact on reducing reservoir quality than compaction, with quartz, dolomite, and illite representing the most abundant types of cement. The analysis of fractures using core and image logs suggests that the fracture density is driven by the lithological variability within the Main Buntsandstein and that fracture joints and stylolites locally may contribute to enhancing the system permeability. The integrated assessment of the results allows the development of prospect play maps for the Buntsandstein in the southern Netherlands, addressing uncertainties and providing future recommendations for further exploration and optimizing geothermal operations in the Triassic. @en

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.

@en

A geothermal well doublet, designed with two primary aims; one of research and the second of commercial thermal energy supply, is currently being installed on the campus of Delft University of Technology, with the wells being drilled in the second half of 2023. The project includes a comprehensive research program, involving the installation of a wide range of instruments alongside an extensive logging and coring program and monitoring network. The doublet has been cored, with continuous samples from the heterogenous reservoir being complimented with more distributed side-wall cores, alongside a large suite of open-hole well logs in the reservoir section of both wells. Such investigation is rarely undertaken in geothermal projects. A fiber optic cable will monitor the production well, and will be installed all-the-way down to the reservoir section when the well completion is installed, at approximately 2300m depth. The reservoir is the fluvial Lower Cretaceous Delft Sandstone that is used as a geothermal reservoir in a series of existing and planned doublets in the West Netherlands Basin. A local seismic monitoring network has been installed in the surrounding area with the aim of monitoring very low-magnitude natural or induced seismicity. A vertical observation well with electromagnetic sensors will be drilled in a few y ears’ time between the injector and producer to monitor cold-front propagation. The total project is targeted to supply around 25 MW of thermal energy at peak conditions, next to this project a thermal energy storage system is planned to provide a seasonal buffer. The project is a key national research infrastructure and is being incorporated into the European infrastructure EPOS (European Plate Observing System, https://www.epos-eu.org/), such that accessibility and data availability will be as wide as possible. All observations will be included in a digital-twin framework that will allow better decisions to be made in future geothermal projects. This paper presents the implementation and initial data collection from the project, including an initial evaluation of the logging and coring campaigns.@en

Aquatic biodiversity hotspots often emerge in regions with active tectonism, diverse climate conditions and complex basin configurations enabling episodic biotic isolation and exchange. The Anatolian microcontinent, located between the Mediterranean and Pontocaspian regions, has been considered a cradle of biodiversity for continental aquatic organisms. The Denizli Basin succession of SW Anatolia contains a “Didacna” mollusc fauna that could be the precursor of the modern Pontocaspian mollusc faunas of the Black Sea-Caspian Sea regions. However, the appearance of Pontocaspian faunas in the Denizli Basin and constraints upon their ages and dispersal pathways remain enigmatic. Moreover, the emergence of the Pontocaspian biota far into the Anatolian continental interior raises questions regarding the connectivity history and tectonic evolution of the Anatolian, Aegean and Pontocaspian realms. Here, we present an integrated stratigraphy of the ∼1 km thick succession of the Kolankaya Formation of the Denizli Basin, previously assigned to the Late Miocene. To date the first occurrence of Pontocaspian fauna in the Denizli Basin and to characterise accompanying palaeoenvironmental/palaeohydrological changes, we focus on three sets of approaches: dating (magnetostratigraphy and ⁴⁰Ar/³⁹Ar), biotic record (molluscs, ostracods and dinoflagellates) and hydrological connectivity (O- and C-isotopes and ⁸⁷Sr/⁸⁶Sr). We date the studied section as Early Pleistocene, spanning a time range of 2.6 Ma to 0.7 Ma. During that time, the Denizli Basin hosted an isolated to partially hydrologically open oligo-to mesohaline lake. The biotic record shows a drastic turnover of mollusc fauna from endemic Aegean-Anatolian and Pannonian/Paratethyan to Pontocaspian affinity at ∼1.8 Ma. The palaeogeographic evolution of the region, along with the geographically limited appearance of the Pontocaspian faunas, suggests a dispersal pathway from the Black Sea Basin via the Aegean Basin. Subsequently, a short incursion into the Denizli Basin may have occurred via a series of graben-type basins: either via the Söke-Milet Basin – Büyük Menderes Graben or via Izmir Bay – Gediz Graben. Our study shows that the Denizli Basin was not a cradle but rather a sink of the Pontocaspian biota during the Early Pleistocene. The new Early Pleistocene age assignment for the Pontocaspian fauna and the Kolankaya Formation in Denizli calls for a major reappraisal of models for the tectonic and stratigraphic evolution of SW Anatolia, including the regional interbasinal connectivity history.

@en

The lower Triassic Main Buntsandstein Subgroup represents a promising, but high-risk geothermal play in the Netherlands. Although the gross thickness in boreholes locally exceeds 200 m, the spatial distribution, geometries and preservation of these sedimentary units remained uncertain due to the lack of seismic data with sufficient resolution and the sparse borehole network. This creates uncertainty in the quantification of the aquifer dimensions that is essential for the planning of geothermal operations.

In this study, seismic interpretation and 2D palinspastic restoration of new and reprocessed seismic data were conducted and combined with borehole data to assess the tectonic evolution of the Roer Valley Graben in the southeastern Netherlands and its control on the spatial distribution of the Main Buntsandstein Subgroup sediments. Our results show that the central and southern parts of the Roer Valley Graben were active depocenters in the Early to Middle Triassic times dominated by fluvial sandstone deposition, providing important play elements for prospective leads on geothermal exploration. The northern part of the basin was a more marginal area where mostly fine-grained sediments were deposited. To the northwest, differential subsidence resulted in the development of areas where the Buntsandstein thickness is reduced to ∼150 m.

After deposition, the Main Buntsandstein sediments were compartmentalised by faulting related to post-depositional tectonic activity, locally reducing the lateral extent of the geothermal target areas down to 1–2 km in a ∼NE–SW direction. On the platform areas adjacent to the Roer Valley Graben and to the southeast, Jurassic sediments are largely absent and the Main Buntsandstein sediments are present at depths shallower than 2 km. These platforms are promising targets for further investigation, as the relatively shallow burial depths, compared to the central part of the Graben, may have contributed to the preservation of more favourable reservoir properties.@en

Nearly half of the Netherlands’ natural gas consump tion is allocated to heating, with direct -use geothermal heating being one of the available low-carbon energy solutions. A geothermal well doublet, designed with the two primary aims of research and commercial heat supply, is currently being installed on the campus of Delft University of Technology. The project is a key national research infrastructure and is being incorporated into the European sustainable and distributed infrastructure (EPOS: European Plate Observing System, https://www.epos-eu.org/), such that accessibility and data availability will be as wide as possible. All observations will be included in a digital-twin framework, which will allow us to make better decisions in future geothermal projects. The project includes a comprehens ive research program, involving the installation of a wide range of instruments alongside an extensive logging and coring program and monitoring network. The doublet has been cored, with substantial continuous samples from the heterogeneous reservoir, alongside a large suite of well logs in both the reservoir and overlying geological units. Such investigation is rarely undertaken in geothermal projects. A fiber-optic cable will monitor the producer well all the way down to the reservoir section, at approximately 2300m depth, in the Lower Cretaceous Delft Sandstone that is used as a geothermal reservoir in a series of existing and planned doublets in the West Netherlands Basin. A local seismic monitoring network has been installed in the surrounding area with the aim of monitoring very low-magnitude natural or induced seismicity. A vertical observation well with electromagnetic sensors will be drilled in the near future between the injector and producer to monitor cold-front propagation. This paper presents the initial modeling for the project and steps towards the production of a digital twin. Two modeling examples in the paper will emp hasize current operational challenges relevant to the project.@en

A geothermal doublet has been installed in a sedimentary reservoir for direct-use heating on the TU Delft campus, targeted to supply around 25 MW of thermal energy at peak conditions. This contribution presents the implementation and initial data collection from the doublet, including an initial evaluation of the logging and coring campaign. Nearly half of Netherlands natural gas consumption is allocated to heating, and the on-campus CO2 emissions from heating exceed 50%. The doublet has been designed with two primary aims of research and commercial heat supply, with the wells being completed in December 2023. The project will be operated by a commercial entity, and built into a larger thermal energy system including a high temperature underground storage system, with the first energy production planned in 2025. The research questions relate to field-scale geothermal operations, e.g. how reliable is the long-term energy production?, how do materials perform in the long-term? and how can geothermal projects be best monitored? The research programme involves the installation of a wide range of instruments alongside an extensive logging and coring program and monitoring network. The doublet has been cored, with substantial continuous samples from the heterogenous reservoir, alongside a large suite of open hole well logs in the reservoir and through casing logs in overlying geological units. A fiber-optic cable will monitor distributed pressure throughout the producer reservoir section, at approximately 2300m depth, which will be installed during commissioning. A local seismic monitoring network has been installed in the surrounding area with the aim of monitoring very low-magnitude natural or induced seismicity. The project is a key national research infrastructure and is being incorporated into the European EPOS (European Plate Observing System, https://www.epos-eu.org/), such that accessibility and data availability will be as wide as possible. All observations will be included in a digital-twin framework that will allow to make better decisions in future geothermal projects.@en

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.

@en

The Lower Triassic Main Buntsandstein Subgroup represents one of the most promising deep geothermal plays in the Netherlands. Reservoir zones consist of sandstone units deposited in an arid to semi-arid alluvial plain reaching total thicknesses of over 200 m in the Roer Valley Graben. However, the lack of a comprehensive reservoir model integrating stratigraphy, sedimentology, and petrography makes the Buntsandstein a high-risk target. Here, we present a reservoir architecture model of the spatial and temporal variation of sedimentary facies and inherent permeability barriers and baffles, based on the integration of data from 33 wells from different stratigraphic levels and different areas of the Roer Valley Graben. Core samples and thin sections analysis revealed two major reservoir facies that were interpreted as the products of transport and depositional processes in braided and sinuous river settings. The identified facies were then coupled to wireline logs to assess the spatial and temporal variation in sedimentary architecture. The lower part of the stratigraphy is dominated by braided river reservoir facies with a high degree of connectivity, where regional lacustrine-playa lake sediments represent the main potential permeability barriers. By contrast, the upper part of the stratigraphy is characterized by an increase in the proportion of sinuous river complexes. These latter yield a lower degree of connectivity with different types of baffles such as intercalated fine-grained overbank sediments, abandonment plugs, bar-draping fines, and cemented dolocrete scour fills. These are much more localized compared to the braided complexes-related barriers, making the prediction of the upper stratigraphy architecture uncertain. @en

The Buntsandstein subgroup in the southeastern part of the Netherlands represents one of the most promising, but risky, geothermal plays. To understand the main controls on Buntsandstein reservoir quality, we combine petrophysical (porosity and permeability) and petrographic (point counting) data derived from different wells and different depth levels. Results show that porosity ranges from 2 to 18.5 and permeability from 0.001 to 285 mD. Dolomite represents the most abundant cement and show an inverse correlation with porosity. Illite occurs in higher concentrations in samples with values of permeability below 20 mD, while kaolinite becomes the most dominant phyllosilicate cement in samples with higher permeability. By looking at the main cement distribution over the sedimentary facies, it appears that dolomite is strongly related to depositional facies and has a positive correlation with grain size, while illite and kaolinite yield a negative correlation with grain size. Pedogenic dolomite nodules are often reworked as detrital grain into the channel scour deposits and are the main source for dolomite cementation. The current study has shown how diagenesis makes Buntsandstein reservoir complex and heterogeneous, and how reservoir quality is strongly related to the depositional environment.@en

In this study, a re-evaluation is performed of the well data of the NLW-GT-01 and VAL-01 wells in the Lower Triassic sandstones in the West Netherlands Basin. Core, geophysical and image logs, are compared to document the characteristics of natural fractures distribution, and investigate the geological parameters influencing their development. The main control on the distribution is identified. Natural fractures are concentrated in heterolithic lithological intervals of the VAL-01 and NLW-GT-01 wells. The identification of more fractures in VAL-01 compared to NLW-GT-01 can be explained by the difference in basin location. VAL-01 was located in the centre of the basin where distal playa environments produced fine-grained material alternating with coarser sands. The more proximal NLW-GT-01 was dominated by fluvial sands. The lithological variability produces Young’s modulus variability that seems to be driving increased fracture density rather than that the absolute value of the Young’s modulus. These fluctuations could be the result of concentrated compressional stress, which is supported by the existence of stylolites, indicative of high compressional stresses in the same intervals. The identification of controls on fracture distribution in the targeted stratigraphic interval can be used to optimize the planning of future geothermal doublets and to de-risk upcoming operations.@en

Orbital driven climate control on sedimentation produces regional, stratigraphically repetitive characters and so cyclostratigraphic correlation can improve correlation and identify stratigraphic trends in borehole sections. This concept is commonly used to correlate marine and lacustrine strata. However, in the alluvial domain, its use is more challenging because internal, local dynamics controlling sedimentation may interfere with the expression of cyclic climate forcing. Intervals of low net-to-gross may be important for successful application in this domain as they tend to better document regional changes. This study applies climate-based stratigraphic correlation concepts to improve well correlations, characterise vertical sand distribution, and identify potential reservoir targets in a generally low net-to-gross interval. Coarsening upward sedimentary repetitions (cyclothems) are identified and correlated with high certainty in nineteen well sections in the upper Carboniferous Westoe and Cleaver formations of the Silverpit Basin. Local sedimentary dynamics provide variability in the character of the cyclothems and several types of cyclothem are classified. Correlation of sections using cyclothems recognised on wireline logs is done twice: once manually and once semi-automatically. The semi-automated correlation is based on calculation of deviation curves which depict stratigraphic changes that are less dependent on absolute wireline values and follow vertical trends more clearly. The correlations provide composite stratigraphies that are analysed using vertical proportions curves. Both approaches yield similar results in terms of stratigraphic trends. However, for detailed correlation of wells, the manual correlation is better at accounting for any local variability within the system. The same two zones of higher net-to-gross ratios are found using both correlation methods. These are linked to palaeoclimatic changes driven by long eccentricity and the proposed climate stratigraphic model has predictive value for identifying sandstone occurrence. The climate-based stratigraphic correlation improves the assessment reservoir distribution and properties on small (10–20 m thickness) and large (100–200 m thickness) stratigraphical scales.@en

The Paleocene-Eocene Thermal Maximum (PETM) global warming event at ∼56 million years before present changed catchment weathering and erosion. Increased chemical weathering of silicate minerals is thought to be an important process removing CO₂ from the atmosphere. However, changes in clay mineralogy can often be explained by enhanced erosion of catchment laterites during the event. Here, we investigate chemical and physical weathering and erosive flux changes through the PETM interval in the Bighorn Basin, Wyoming, a Laramide foreland basin, in a proximal continental-interior alluvial setting. These show an increase of detrital smectite with a lag time of 20-kyr after the main onset the PETM. The smectite increase continued for at least 50-kyr after the event. In-situ, post-depositional pedogenic clay mineral formation is similar between pre-PETM and PETM soil profiles, despite large macroscopic differences between soils that formed before and during the event. Drier, hotter summers during the PETM probably caused decreased vegetation cover that, in concert with more frequent and heavier rainstorms, intensified the erosion of smectite-rich Cretaceous bentonites on the margins of the catchment, which exceeded changes in chemical weathering within the catchment. The lagged response in reaching full PETM clay mineral values can be explained by the time required for upstream sediment to reach the catchment basin floodplain. The prolonged nature of smectite enhancement after the PETM event may again relate to signal propagation times that are now even longer due to lower fluvial recycling rates. Our results indicate that chemical weathering changes were probably superceded by enhanced physical weathering and clay-mineral transport from basin margins at this continental-interior study site.

@en

The multi-purpose research borehole at the Delftse Hout is the third of four seismic monitoring locations of the seismic monitoring network for the geothermal research project on the TU Delft campus (Geothermal Delft GTD, also known as DAPwell, https://geothermiedelft.nl/). For the geothermal research project, two deep wells (“a doublet” consisting of an injector and a producer) for geothermal energy extraction will be installed on the TU Delft campus next to the combined heat and power plant (“warmtekrachtcentrale - WKC”). The system will produce geothermal heat to supply the campus of TU Delft and part of the city of Delft.
The herein presented borehole describes the installation of a multi-purpose research borehole (called DAPGEO-02), which was installed in the period February - May 2022. DAPGEO-02 is part of a seismic monitoring system for the shallow and deeper subsurface in the vicinity of the planned geothermal doublet. The locations of all four stations are given in Figure 1. The monitoring network and the related research gathers knowledge about the current status of the subsurface on the basis of periodic data measurements, and possible seasonal effects.
Within the seismic monitoring network, three seismic monitoring stations have already been installed, respectively DAPGEO-01 on the proposed location of the geothermal project near the Leeghwaterstraat in Delft, DAPGEO-03 on the Kerkpolderweg in Delft, and ZH03 in on the Ackersdijkseweg in Pijnacker-Nootdorp (installed and equipped by KNMI).@en

ABSTRACT The lower Eocene Willwood Formation of the Bighorn Basin, Wyoming, USA, is an alluvial succession with a sand content varying around 25 palaeoenvironments and palaeoclimates, as well as sedimentological and stratigraphic analysis. Channel dynamics were studied at a relatively low resolution throughout the basin over the geological time from late Palaeocene to early Eocene. Here, a high-resolution study is reported to complement previous research at the basin scale. Efforts are made to document the characteristics and river planform styles of most sandstone bodies encountered through ca 300 m of alluvial stratigraphy in a 10 km2 area of the Deer Creek part of the McCullough Peaks area situated in the basin axis of northern Bighorn Basin. Four channel facies associations are recognized and ascribed to four river planform styles: crevasse channel, trunk channel, braided-like channel and sinuous-like channel, with the latter two types dominant. Braided-like and sinuous-like channel sandstone bodies differ significantly in thicknesses, being on average 6.1 m versus 9.0 m, but they have similar palaeoflow–perpendicular widths of on average 231 m and palaeoflow directions of on average N 003°. Braided-like and sinuous-like river planform styles show no spatial dependency in the 10 km2 study area. Results of this study are in line with existing basin-scale depositional models that are composed of a single axial system fed by several transverse systems dominantly from the west. The feeding of these systems could be influenced by palaeoclimate changes possibly controlling their contribution over time, thereby impacting river planform styles. At the same time, changing water discharge hydrograph, sediment load, and overbank cohesiveness may have equally driven the observed river planform style changes within the basin without a major role of catchments.@en

The ambition is to significantly enhance the knowledge, mapping, and prediction of the geological character of Dutch geothermal plays, focussing on thermal and geomechanical properties. The aim is to improve fundamental understanding of the geological causes behind rock properties, the understanding and quantification of the conversion of downhole log responses into rock properties, the statistical approaches to study these properties, and the implementation of rock properties into reservoir models. The study aims at leading to lower research costs for geothermal operators, lower uncertainties concerning production prediction and risk assessment, and improved productivity due to more optimal well placement and production strategies. With that, our fully open-access results will be applicable to all geothermal targets in the Netherlands and thereby be able to calibrate reservoirs, geomechanical and thermal reservoir models with the ultimate goal to optimise the exploitation of geothermal heat in the Netherlands in a sustainable and safe way. The project has a measuring component producing data and geological understanding, a correlating component, linking rock properties to petrophysical log data by various innovative means, and an implementation component setting the findings into geomodel application. Novel microstructural scanning data enables to link the nano-scale rock composition to decimetre scale bedding and logs. The focus will be on pore throats, detrital, authigenic, and new-grown cement and their types and how these relate to flow, and in particular thermal, and mechanical properties of the rocks. Next, petrophysical downhole logging data will be analysed using statistical and machine learning techniques to produce a much enhanced methodology to relate petrophysical log responses to different rock properties. Improved correlations will be produced per play investigated and can be applied to clastic geothermal reservoirs in general. This will allow for quick and improved screening of rock properties through different wells and finally beyond those into white spots. @en

A combined study of depositional facies and diagenesis variation was carried out to understand the main controls on aquifer quality of the Middle Buntsandstein in the southeastern part of the Netherlands. Heterogeneities in continental sandstone bodies occur at different spatial scales, ranging from micrometers to hundreds of meters. Commonly, such heterogeneities result from the interaction of depositional processes at various spatial and time scales. These processes partially also influence subsequent diagenetic evolution, hence present-day aquifer properties. Understanding the role of the resulting architectural heterogeneities in controlling the dynamic reservoir behavior is key in determining aquifer properties and improving pre-drilling prediction. The sandstones of the Main Buntsandstein subgroup in the southeastern part of the Netherlands provide an excellent example where different detrital compositions, internal sedimentary architectures, and diverse burial histories have resulted in a wide range of present-day aquifer properties. In the study area, the aquifers are composed of stacked heterogenous alluvial sandstones bodies intercalated with mud-prone intervals deposited in arid to semi-arid conditions. Differences in sediment sources, transport mechanisms, and intrabasinal conditions resulted in a wide distribution of composition and texture. Additionally, the effect of post-depositional burial diagenesis in a basin with complex tectonic history created diverse burial histories across the basin. The study aims to investigate the variation of present-day aquifer hydraulic parameters about changes in aquifer facies and architecture, detrital composition, as well as compaction and cementation during burial. Core sample analysis unfolded a diverse spectrum of sedimentary facies and lithic fragments, which differ between formations. Thin section analysis provides insights about mechanical compaction, cementations, and authigenic phases. By combining these results with petrophysical data on permeability and porosity of core samples, the major controls on present-day aquifer quality can be assessed. @en

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.

@en

Smectites affecting water and soil exchangeable cations, especially potassium and ammonium related to plants, are common swelling clay minerals in sediments and paleosols that can record paleoenvironmental and paleoclimatic information. Paleoenvironmental and paleoclimatic interpretations are potentially ambiguous, because smectite minerals may have either detrital or pedogenic origins. This study analyzes smectites in Lower Eocene paleosols of the Qaidam Basin, China, with the goal of determining their origin and paleoclimatic significance. Twelve smectite-rich vertisol samples were analyzed by X-ray diffraction (XRD), Fourier transform infrared spectroscopy (FTIR), inductively coupled plasma-mass spectrometry (ICP-MS), scanning electron microscopy (SEM), and differential scanning calorimetry (DSC) in the Dahonggou (DHG) section of the Lulehe Formation. The samples from DHG section contained predominant dioctahedral smectites (60–65%) in the <2 μm clay fractions. Minor dioctahedral illites, dioctahedral interstratified illite-smectite (I-Sm), kaolinite, and traces of dioctahedral chlorites and palygorskites were present. The chemical data and FTIR spectra revealed both tetrahedral substitutions (i.e., Al for Si) and octahedral substitutions (i.e., Mg and minor Fe for Al) in the dioctahedral smectites. Montmorillonite comprised the dominance of the fine-clay fractions (< 0.2 μm) with minor contributions of beidellite, as determined by the Greene-Kelly test (GK-test), geochemical and FTIR analyses. SEM observations revealed two morphological types of dioctahedral smectites, the less abundant one with typical ‘honeycomb’ structure, interpreted as neoformed, and the other occurring as thin, well-defined plates and irregular masses considered to be detrital. DSC revealed two dehydroxylation temperatures (500–550 °C and 600–700 °C), indicating the presence of trans-vacant (tv) and cis-vacant (cv) layers, respectively. The cv layers were attributed to montmorillonite originating as catchment-delivered detritus, as the study Lulehe Formation exhibited no evidence of volcanic ash layers or diagenesis that would have promoted formation of montmorillonite with high dehydroxylation temperatures. The tv layers represented pedogenic beidellite formed via solution crystallization in a floodplain pedogenic environment. The accumulation of abundant inherited montmorillonite and neoformed beidellite in paleosols of the Lulehe Formation was probably promoted by a relatively warm climate with alternating wet-dry seasons during the Early Eocene. Admixture of detrital montmorillonite makes paleoclimatic reconstructions based on bulk-smectites concentrations and isotopic compositions unreliable, showing that integrated clay mineralogical analyses of the type undertaken in this study are necessary to identify the provenance and paleoclimatic significance of smectites in sediments.

@en

Asian mineral dust has been studied extensively for its role in affecting regional-to global-scale climate and for its deposits, which enable reconstructing Asian atmospheric circulation in the past. However, the timing and origin of the dust deposits remain debated. Numerous loess records have been reported across the Asian continent with ages varying from the Miocene to the Eocene and linked to various mechanisms including global cooling, Tibetan Plateau uplift and retreat of the inland proto-Paratethys Sea. Here, we study the Eocene terrestrial mudrocks of the Xining Basin in central China and use nonparametric end-member analysis of grain-size distributions to identify a loess-like dust component appearing in the record at 40 Ma. This is coeval with the onset of high-latitude orbital cycles and a shift to predominant steppe-desert vegetation as recognized by previous studies in the same record. Furthermore, we derive wind directions from eolian dune deposits which suggest northwesterly winds, similar to the modern-day winter monsoon which is driven by a high pressure system developing over Siberia. We propose that the observed shifts at 40 Ma reflect the onset of the Siberian High interacting with westerly derived moisture at obliquity timescales and promoting dust storms and aridification in central China. The timing suggests that the onset may have been triggered by increased continentality due to the retreating proto-Paratethys Sea.

@en