Horseshoe crabs (Xiphosura) have a long evolutionary history starting in the Ordovician, but they have rarely been reported from the Netherlands. We report on the first Triassic horseshoe crab from the Netherlands identifiable to the species level, a specimen of the limulid Limulitella bronnii. We provide the first diagnosis for this species and refigure the holotype. The new specimen was found in the Middle Triassic (Anisian) Muschelkalk sediments of the Vossenveld Formation, in the Illyrian part of the stratigraphic profile of the Winterswijk quarry complex. The Winterswijk specimen represents the youngest occurrence of L. bronnii. The inferred non-marine habitat of this horseshoe crab species elsewhere in conjunction with occurrences of plant and insect remains within the same layer at Winterswijk suggest the specimen herein most probably did not live in marine conditions either. This species has previously been found in non-marine sediments in France and Germany, expanding its geographic range northward. Several faunal elements from Winterswijk including L. bronnii show resemblance to the roughly co-eval non-marine components of the Anisian Grès à Voltzia Formation in NE France, suggesting a paleobiogeographic connection between these regions in Western Europe.

We report a marine Middle Pleistocene mollusc fauna from a borehole near Luxwoude (Friesland, northern Netherlands). The fauna contains several species, including Bittium species, Acanthocardia paucicostata and Polititapes senescens, that hitherto have been used as indicative for warm (lusitanian) conditions in the southern North Sea Basin during the Late Pleistocene Eemian (MIS5e) interglacial. However, the stratigraphic context of the Luxwoude marine fauna indicates a MIS11 or older age for this new fauna. This thermophilous fauna demonstrates very warm-temperate conditions and probably an open marine connection through the Dover Strait towards the south, during this Middle Pleistocene interglacial. In the North Sea basin, this distinctive lusitanian fauna with Bittium-dominated assemblages can therefore no longer be presumed to be of Eemian age without additional evidence.

The Central Graben is a Mesozoic sedimentary basin that is significantly influenced by rift and salt tectonics. Its southern part is located in the German and Dutch sectors of the North Sea. Even though studies exist on the tectonic and stratigraphic development of the Danish and Dutch Central Graben, the German Central Graben as an important link is less investigated. We aim to fill this gap and to investigate the sedimentary development from the Latest Triassic to the Early Cretaceous, the relative influence of salt and rift tectonics on subsidence and how our results fit into the existing studies of the Danish and Dutch Central Graben. Knowledge of the development of the graben and its sedimentation is critical for any possible economic use like hydrocarbon exploitation or carbon capture and storage. Therefore, we mapped nine laterally traceable horizons on 2D and 3D reflection seismic data from the Lower Jurassic to the Lower Cretaceous within the German Central Graben and adjacent Danish Salt Dome Province as well as the northern Dutch Central Graben. These horizons include the base horizons of four tectonostratigraphic mega-sequences of the southern Central Graben adopted from the current Dutch tectonostratigraphic concept. Based on the mapping results, we constructed subsidence, thickness and erosion maps of the tectonostratigraphic mega-sequences and their subdivisions. The tectonostratigraphic mega-sequences were then correlated with well logs to determine the lithology. The results show that the structural and stratigraphic architecture of the German Central Graben was consecutively dominated by either subsidence controlled by rifting, salt tectonics or by thermal uplift and subsidence. We suggest that the German Central Graben is divided by a large strike-slip fault zone, the Mid Central Graben Transverse Zone, into a northern part that geologically rather belongs to the Danish and a southern part that rather belongs to the Dutch Central Graben. We discuss how this division and the tectonics influenced the regional lithology.

Scour holes are common features in deltaic rivers which can destabilise embankments through oversteepening of the river bed. Their development has been studied extensively from the hydraulic perspective, but another important control is the erodibility of the river bed which varies considerably due to thickening of heterogeneous deltaic substrate towards the coast. Therefore, we assessed the influence of delta-scale geological heterogeneity and local subsurface architecture on scour hole formation in addition to the hydrodynamic controls. We (1) created an inventory of 165 scour hole locations in the Rhine–Meuse delta, (2) assessed the hydrodynamic conditions at the locations, (3) extracted geometric characteristics and (4) determined the subsurface architecture from geological data. Central and lower delta branches have 0.6–0.7 scours per km while upper delta branches have less than 0.2. Downstream, 58% of scour holes were related to architectural elements, notably sand bodies from former Holocene channel belts and Early Holocene cohesive beds. These scours have steeper slopes due to higher proportions of cohesive sediments near the river bed. Furthermore, scours related to channel belt sand bodies are limited in downstream length and depth, up to maximum of approximately two times the water depth. From our results, we provide a delta-scale explanatory framework that relates the position of present-day river channels with respect to Pleistocene river deposits and Holocene fluvio-deltaic deposits to scour hole formation. Upstream rivers are incised in Pleistocene deposits showing less local variation in erodibility. The majority of scour holes here relate to engineering works. In central and lower delta branches, geologically inherited heterogeneity of the Holocene substrate at critical depths near the channel bottom adds to anthropogenic induced scours and results in high abundances. This demonstrates that downstream variation in subsurface architecture should be considered as a key control on scour locations and characteristics for management purposes.

A mandible of a Late Cretaceous sea turtle with affinities to Ctenochelys is reported from the Maastrichtian type area of the Netherlands. The triangular mandible has a well-developed symphyseal ridge surrounded on both sides by large, concave areas on the triturating surface. It represents the first potential occurrence of Ctenochelys from the Maastrichtian type area. This finding increases the diversity of the turtle fauna known from the Maastrichtian type area.

The Netherlands is in the midst of an energy transition with hydrocarbon production gradually declining, whereas the role of sustainable energy technologies is on the rise. One of these technologies is geothermal energy production from porous reservoirs at 1.5–3 km depth. As the number of geothermal projects increases, there is a growing concern that felt and/or damaging induced seismic events could occur as a result of geothermal operations. Over the last two decades, such events have occurred in the Netherlands due to gas production, notably in the Groningen gas field. However, the occurrence of felt events is limited to hydrocarbon fields in certain regions or reservoirs. Understanding where and for which plays these events are observed helps to estimate seismogenic potential for geothermal operations and other sustainable subsurface activities. Here, we summarise and review the main similarities and differences in terms of geological and geomechanical characteristics between the hydrocarbon and geothermal plays in the Netherlands, and we consider the differences in pressure and temperature changes. By doing so, we provide better insights into the factors that could play a role for fault reactivation and induced seismicity, and how these differ for hydrocarbon production and geothermal operations in the Netherlands. The review shows that geological characteristics for most geothermal target reservoirs are similar to those of hydrocarbon, albeit geothermal projects so far target higher porosity rocks than hydrocarbon reservoirs. On the other hand, pressure and temperature changes are very different, with significant depletion for hydrocarbon fields vs significant cooling around geothermal injection wells. The different operations result not only in different expected stress change magnitudes but also in a distinct spatio-temporal stress build-up on faults, which has implications for seismogenic potential and monitoring of these different operations.

The utilisation of geothermal energy in the Netherlands is primarily focused on deep sedimentary aquifers, which are often intersected by major faults. Geothermal operations (i.e. fluid production and injection) may alter the effective stress state along these faults and trigger induced seismic events. Pore pressure perturbations have been generally considered the main driver of injection-induced seismicity. However, thermal stresses caused by temperature gradients between the re-injected cold fluid and the reservoir rock may also contribute to the triggering of earthquakes in geothermal reservoirs. While existing geothermal power plants operating in sandstone reservoirs did not produce any major induced seismicity, it is a matter of debate whether a reduction in the temperature of the re-injected fluid could increase the seismic hazard potential. In this study, we applied modified Gutenberg–Richter statistics based on frictional Coulomb stress variations implemented in a coupled thermo-hydro-mechanical model to estimate the seismic hazard caused by the operation of a geothermal doublet. We conducted a systematic parametric study to assess and rank the impact of different intrinsic (geological) and extrinsic (operational) parameters on the induced seismic hazard potential. We identified a competing mechanism between induced variations in pore pressure and thermal stress within the reservoir in controlling induced seismicity. We found that stress changes induced by pore pressure variations are the main cause of seismic hazard, although thermally induced stresses also contribute significantly. The results indicate that by optimising the operational parameters it is possible to increase production efficiency while maintaining a long-term control over the fluid injection-induced seismicity.

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.

Chondrichthyan and actinopterygian fish remains from Rhaetian (c. 208.05–201.36 Ma) or perhaps Late Norian deposits in the Winterswijk quarry are described. The most abundant taxon is the actinopterygian Gyrolepis albertii, followed by the chondrichthyan Lissodus minimus. Furthermore, the palaeopterygian actinopterygians Saurichthys longidens and Birgeria acuminata, and some teeth of neopterygians Sargodon tomicus, ‘Lepidotes’ sp. and indeterminate pycnodontiforms are recorded in addition to the chondrichthyans Rhomphaiodon minor, Parascylloides turnerae and some ‘Hybodus’ cf. cuspidatus (senior synonym of H. cloacinus). Chondrichthyan dermal denticles, actinopterygian scales and gill rakers, tooth plates, and some fish bones were also found. There is considerable faunal resemblance to the various localities from the Rhaetian of the British Penarth Group, although it depends on the location as to whether chondrichthyans or actinopterygians prevail in the samples. On average, there are more chondrichthyan teeth present in the British samples than actinopterygian teeth, which is opposite to the situation in Winterswijk. That might be explained by different ecological circumstances, such as lower oxygen levels in bottom waters in Winterswijk and freshwater input and/or changes in salinity in the UK.

From dynamic rupture simulations, we reveal under which conditions a rupture in the Groningen gas field stops along fault dip or along fault strike after it starts on a fault in the reservoir. The simulations focus on the capabilities of fault plane irregularities to arrest ruptures. Such irregularities can be recognised in sandstone outcrops. Fault planes in the Groningen field, extracted from the 3D seismic data by seismic attribute extraction methods, show similar irregularities. A detailed surface of a major fault plane in the field indicates that steps and jogs of tenths of metres are possible. Although these irregularities are close to seismic resolution and could be partially artificial, we investigated their effect on rupture arrest.

For typical current stresses in the Groningen field, jogs and steps of this length scale are found to be remarkably effective to stop ruptures in the reservoir. Also, a significant increase in the fault dip along fault strike can stop these ruptures but a kink in the fault under a constant fault dip not.

Including non-planar fault features and pressure diffusion in the Carboniferous, the simulations in this paper follow trends of previous simulations in the literature using 2D planar faults. In particular, the horizontal stress in this formation and the strength of the Carboniferous fault zone are important for rupture propagation. If the fault would have been reactivated in the Neogene or Quaternary and poorly healed in clay-rich parts, rupture propagation into the Carboniferous can only be prevented by jogs of sufficient size and lateral continuity under the present estimate of the horizontal field stress.

In the Western Scheldt estuary, like in many estuaries, safe navigation, flood protection, and ecological targets require a balanced and sustainable sediment management. A thorough understanding of the morphodynamic functioning of the estuary and its response to changes in hydrodynamics (natural sediment transport) and large-scale interventions is imperative. This paper presents a detailed overview of over 65 years of morphological changes and a comprehensive sediment budget of the Western Scheldt estuary that is based on analysis of a unique series of frequent bathymetric surveys and historical data on human–sediment interactions of dredging, dredge disposal and sand mining. Solving the sediment budget reveals an annual sediment import of 2.2 million m3. The highest sediment imports occurred between 1980–1994 and 2005–2020 (2.9 and 3.7 million m3/year). A 1.3 million m3/year net export prevailed between 1994 and 2005. Natural variations in the hydrodynamics (e.g., tidal asymmetry and amplification) and sediment transports cannot explain the derived temporal variations in sediment import rates. Anthropogenic activities play a dominant role. Capital dredging of the main navigation channel has led to channel deepening and significantly increased dredge and disposal volumes. Disposal on tidal flats and in the secondary channel leads to a long-term storage of sand and, consequently, a local increase in bed level and a sand deficit in the remainder of the system that results in increased sediment imports. Large-scale disposal in the western part of the estuary can (temporarily) reverse the sediment exchange between the estuary and its mouth area, as observed between 1994 and 2005. Apparently, large-scale anthropogenic reallocation of sediment by dredging and/or disposal as part of navigation channel improvement, sand mining or nourishment essentially influences the morphological development of the Western Scheldt estuary. This reveals responsibilities as well as opportunities of sediment management for the Western Scheldt and similar estuaries worldwide.