This special issue of the Netherlands Journal of Geosciences / Geologie en Mijnbouw contains research papers offering important new insights into aspects of the occurrence of gas-production-induced seismicity in Groningen, the Netherlands. The issue is a ‘work in progress’ and provides a snapshot of our knowledge at the end of 2016, which is when most of the contributions were written.

The up to 60 m thick Neurath Sand (Serravallian, late middle Miocene) is one of several marine sands in the Lower Rhine Basin which were deposited as a result of North Sea transgressive activity in Cenozoic times. The shallow-marine Neurath Sand is well exposed in the Garzweiler open-cast mine, which is located in the centre of the Lower Rhine Basin. Detailed examination of three sediment profiles extending from the underlying Frimmersdorf Seam via the Neurath Sand and through to the overlying Garzweiler Seam, integrating both sedimentological and palaeontological data, has enabled the depositional setting of the area to be reconstructed.

Six subenvironments are recognised in the Neurath Sand, commencing with the upper shoreface (1) sediments characterised by glauconite-rich sands and an extensive biota (Ophiomorpha ichnosp.). These are associated with the silt-rich sands of a transitional subenvironment (2), containing Skolithos linearis, Planolites ichnosp. and Teichichnus ichnosp. These silt-rich sands grade up to the upper shoreface subenvironment (1), which is indicative of an initial regressive trend. The overlying intertidal deposits can be subdivided into a lower breaker zone (3), characterised by ridge-and-runnel systems, and the swash zone (4) where the surge and backwash of waves resulted in the deposition of high-energy laminites. The intertidal deposits were capped by aeolian backshore sediments (5). Extensive root traces present in this latter subenvironment reflect the development of the overlying peatland (i.e. Garzweiler Seam). Within the Garzweiler Seam, restricted sand lenses indicate a lagoonal or estuarine depositional environment (6). Regional correlation with adjacent wells establishes that shallow-marine conditions were widespread across the Lower Rhine Basin in middle Serravallian times. The shoreline profile, characterised by both tidal and wave activity and influenced by fluvial input from the adjacent Rhenish Massif, is indicative of the complexity of the coastal depositional setting within the Lower Rhine Basin.

This special issue contains papers presented at the conference ‘The Meuse Valley: 15,000 years of (A)biotic Landscape Development and Human Impact’. The conference was part of the annual meeting (Palynologendagen) of the Palynologische Kring (Dutch Palynological Society) held on 9 and 10 October 2014 at the Castle De Berckt in Baarlo (province of Limburg).

Tidal inlets and the associated ebb-tidal deltas can significantly impact the coastal sediment budget due to their ability to store or release large quantities of sand. Nearly 300 million m3 (mcm) of sediments were eroded from Texel Inlet's ebb-tidal delta and the adjacent coasts following the closure of the Zuiderzee in 1932. This erosion continues even today as a net loss of 77 mcm was observed between 1986 and 2015. To compensate, over 30 mcm of sand has been placed on the adjacent coastlines since 1990, making maintenance of these beaches the most intensive of the entire Dutch coastal system.

Highly frequent and detailed observations of both the hydrodynamics and morphodynamics of Texel Inlet have resulted in a unique dataset of this largest inlet of the Wadden Sea, providing an opportunity to investigate inlet sediment dynamics under the influence of anthropogenic pressure. By linking detailed measurements of bathymetric change to direct observations of processes we were able to unravel the various components that have contributed to the supply of sediment to the basin, and develop a four-stage conceptual model describing the multi-decadal adaptation of the ebb-tidal delta.

Prior to closure of the Zuiderzee a dynamic equilibrium state (stage 1) existed with a stable ebb-tidal delta. The largest morphological changes occurred in roughly the first 40 years since the closure, and were dominated by the rotation and scouring of large tidal channels and landward retreat of the Noorderhaaks ebb shoal (stage 2; adaptation). Between 1975 and 2001 the general layout of main channels and shoals was stable, but large sediment losses continued to occur (stage 3; equilibrium erosional state). Since 2001, the erosion rates have significantly reduced to 2 mcm a−1 (stage 4; stabilisation).

Twenty-five years of data on ‘Dynamic Preservation’ prove that sand nourishments are well capable of keeping the coastlines adjacent to the Texel Inlet in place. Moreover, the abundant supply of sediment may also have compensated for the sediment losses on the larger scale of the southern part of the ebb-tidal delta, resulting in a recent stabilisation of its volume. This response illustrates the potential benefits of Dynamic Preservation not only for coastline resilience but also on the larger scale of the inlet system. Such knowledge is essential for future preservation, management and maintenance of inlet systems in the scope of climate change and accelerated sea-level rise.

The primary challenge for efficient geothermal doublet design and deployment is the adequate prediction of the size, shape, lateral extent and thickness (or aquifer architecture) of aquifers. In the West Netherlands Basin, fluvial Lower Cretaceous sandstone-rich successions form the main aquifers for geothermal heat exploitation. Large variations in the thickness of these successions are recognised in currently active doublet systems that cannot be explained. This creates an uncertainty in aquifer thickness prediction, which increases the uncertainty in doublet lifetime prediction as it has an impact on net aquifer volume. The goal of this study was to improve our understanding of the thickness variations and regional aquifer architecture of the Nieuwerkerk Formation geothermal aquifers. For this purpose, new palynological data were evaluated to correlate aquifers in currently active doublet systems based on their chronostratigraphic position and regional Maximum Flooding Surfaces. Based on the palynological cuttings analysis, the fluvial interval of the Nieuwerkerk Formation was subdivided into two successions: a Late Ryazanian to Early Valanginian succession and a Valanginian succession. Within these successions trends were identified in sandstone content. In combination with seismic interpretation, maps were constructed that predict aquifer thickness and their lateral extent in the basin. The study emphasises the value of palynological analyses to reduce the uncertainty of fluvial hot sedimentary aquifer exploitation.

The fluvial development of the Roer river in the southeastern Netherlands and western Germany is presented for the Late Pleniglacial, Late-glacial and Early Holocene periods. Reconstruction of fluvial-style changes is based on geomorphological and sedimentological analysis. Time control comes from correlation to the pollen-based biochronostratigraphic framework of the Netherlands combined with independent optically stimulated luminescence (OSL) ages. At the Pleniglacial to Late-glacial transition a system and channel pattern change occurred from an aggrading braided to an incising meandering system. Rapid rates of meander migration, as established for the Late-glacial by optical dating, were likely related to the sandy nature of the substratum and the Late-glacial incision of the Meuse that resulted in a higher river gradient in the downstream part of the Roer. In the Roer valley the Younger Dryas cooling is not clearly reflected by a fluvial system response, but this may also be related to Holocene erosion of Younger Dryas fluvial forms. An important incision and terrace formation was established at the Younger Dryas to Early Holocene transition, probably related to forest recovery, reduced sediment supply and base-level lowering of the Meuse. The results of this study show a stepwise reduction in the number of channel courses from a multi-channel braided system in the Pleniglacial, to a double meander-belt system in the Late-glacial and a single-channel meandering system in the Early Holocene. The results show that the forcing factors of fluvial-system change in the Roer valley are climate change (precipitation, permafrost and vegetation) and downstream base-level control by the Meuse.

A recent find of a Middle Villafranchian (c. 2.35–2.10 Ma) Canis cf. C. etruscus in the trawlings from the Oosterschelde concerns the oldest dog known from the Netherlands and is the first appearance record of this canid in the North Sea Basin. It shows that the tribe Canini was dispersed beyond south central Europe up to the northwestern edge of the continent. The find confirms the lack of synchroneity and usefulness of the so-called ‘Wolf Event’.

The estuaries in the SW Netherlands, a series of distributaries of the rivers Rhine, Meuse and Scheldt known as the Dutch Delta, have been engineered to a large extent as part of the Delta Project. The Voordelta, a coalescing system of the ebb-tidal deltas of these estuaries, extends c.10 km offshore and covers c.90 km of the coast. The complete or partial damming of the estuaries had an enormous impact on the ebb-tidal deltas. The strong reduction of the cross-shore directed tidal flow triggered a series of morphological changes that continue until today. This paper aims to give a concise overview of half a century of morphological changes and a sediment budget, both for the individual ebb-tidal deltas and the Voordelta as a whole, based on the analysis of a unique series of frequent bathymetric surveys. The well-monitored changes in the Voordelta, showing the differences in responses of the ebb-tidal deltas, provide clear insight into the underlying processes. Despite anthropogenic dominance, knowledge based on natural inlets can still explain the observed developments. Complete damming of the three northern estuaries Brielse Maas, Haringvliet and Grevelingen resulted in a regime shift, from mixed-energy to wave-dominated, and sediments are transported in landward and downdrift direction. This results in large morphodynamic changes – sediments are redistributed from the delta front landward – but small net volume changes – a 0.1–0.2 × 106 m3 a−1 increase in volume over the period 1965–2010 – since the dams block sediment transport into the estuaries. Sediment volume losses of 106 m3 a−1 are observed on the ebb-tidal delta of the partially closed Eastern Scheldt and still open Western Scheldt estuary. As a result of a reduction of the estuarine tide in the mouth of the Eastern Scheldt, the north–south-running North Sea tidal wave has gained impact on its ebb-tidal delta, which causes morphological adjustments and erosion of the Banjaard shoal area. Moreover, the Eastern Scheldt ebb-tidal delta delivers sediment to its neighbours. The stable ebb-tidal delta configuration in the Western Scheldt, despite major dredging activities, illustrates that these large inlet systems are robust and resilient to significant anthropogenic change, as long as the balance between the dominant hydrodynamic processes (tides and waves) does not alter significantly.

After more than half a century of production and with some 350 wells, the Groningen gas field must be one of the best-studied gas fields in the world. Initially, it was considered to be relatively simple and behaving like one big tank. Now that it is entering a phase of declining production it has become clear that many subtleties are not fully understood yet. Prediction and management of subsidence and induced earth tremors require a detailed understanding of the field geology. In addition, an optimum gas recovery is only possible if details of, for example, reservoir quality distribution and faulting, that did not appear relevant before, are well understood.

The large Groningen field comprises a structurally high block during much of its history, probably already from Devonian times onwards. The desert sandstones of the Rotliegend reservoir exhibit a strong south-to-north proximal–distal relationship. Whilst diagenesis has in many fields led to deterioration of reservoir properties, this effect is small in the Groningen field. The field is dipping to the north, and bounded by a series of normal faults in the west, south and east. Almost all faults are normal extensional faults, but locally inverse reactivation has led to small pop-up structures. Reactivation of older faults must have resulted in oblique movements along most faults.

The challenges for future development of the Groningen field are immense. Managing the risks associated with induced seismicity and recovery of the remaining gas will continue to require an increasingly detailed knowledge and understanding of its geology.

The artefact size of the Early Middle Palaeolithic (EMP) assemblages in ice-pushed Rhine–Meuse deposits in the central Netherlands decreases northwestward. This trend correlates to the downstream fining direction of the Rhine–Meuse fluvial system, the source of the rock material, showing that locally available material was used. Furthermore, also in line with the fluvial trend, the gravel and cobble trends indicate that EMP artefacts could be present in the buried part of the ice-pushed ridges in the northern part of the central Netherlands (southwest Flevoland). Based on the gravel and cobble data, combined with the literature study, we argue that the artefacts date to marine isotope stages (MIS) 7–6. The review of published data shows that during the time of deposition of the artefact-bearing layer (early MIS 6), the Rhine–Meuse study area was part of a braidplain located in a slightly incised valley. To the north a polar desert was present.

In fully developed evaporite cycles, effective viscosity contrasts of up to five orders of magnitude are possible between different layers, but the structures and mechanics in evaporites with such extreme mechanical stratification are not well understood. The Zechstein 3 unit in the Veendam salt pillow in the Netherlands contains anhydrite, halite, carnallite and bischofite, showing this extreme mechanical stratification. The Veendam Pillow has a complex multiphase salt tectonic history as shown by seismic reflection data: salt withdrawal followed by convergent flow into the salt pillow produced ruptures and folds in the underlying Z3-anhydrite–carbonate stringer and deformed the soft Z3-1b layer

We analysed a unique carnallite- and bischofite-rich drill core from the soft Z3-1b layer by macroscale photography, bulk chemical methods, X-ray diffraction and optical microscopy. Results show high strain in the weaker bischofite- and carnallite-rich layers, with associated dynamic recrystallisation at very low differential stress, completely overprinting the original texture. Stronger layers formed by alternating beds of halite and carnallite show complex recumbent folding on different scales commonly interrupted by sub-horizontal shear zones with brittle deformation, veins and boudinage. We attribute this tectonic fragmentation to be associated with a softening of the complete Z3-1b subunit during its deformation. The result is a tectonic mélange with cm- to 10 m-size blocks with frequent folds and boudinage. We infer that these structures and processes are common in deformed, rheologically strongly stratified evaporites.

During the Pleistocene, the coastal marine bivalve mollusc Mya arenaria became extinct in northwest Europe. The species remained present in North America. Datings of Mya shells found in northern Denmark and the southern Baltic Sea suggest that repopulation of northwest European coasts already occurred before Columbus’ discovery of America (1492), possibly facilitated by Viking (Norse) settlers at Greenland and northeast North America. In this paper we report on findings of M. arenaria at five locations in the coastal landscape of the Netherlands: polders reclaimed from the Wadden Sea and the former estuaries of Oer-IJ and Old Rhine. The shells from four of these locations also date before 1492 AD.

This study provides an overview of existing palynological and chronological data of the northern Meuse valley which have been collected over recent decades. The palynological data were used to make a vegetation reconstruction in time and space for the Lateglacial and early Holocene. The vegetation development is strongly influenced by the rapid climate changes that occurred during this time period. It is shown that the biostratigraphy can be used to provide better age estimates for the abandoned channel fills, which have been dated in most cases using conventional bulk 14C dates. Furthermore, the combination with a geomorphological reconstruction based on AHN (actual height model of the Netherlands) lidar data has been used to evaluate the interactions between fluvial terrace formation and vegetation development. It appears that, although the vegetation development is comparable to the general vegetation development in the Netherlands, slight differences occur, in particular of the vegetation composition, presumably linked to the dynamic geomorphological activity in the Meuse valley. Finally, the spatial distribution of sites may give indications for the migration routes of pine (Pinus) and poplar (Populus) during this period of rapid vegetation development over the Last Glacial–Interglacial Transition.

Prediction of gas-production-induced subsidence and seismicity is much more difficult and uncertain than generally recognised in the past. It is now widely accepted that uncertainties in predicted subsidence and seismicity are large prior to and during the initial stages of production. At later stages, predictions remain highly uncertain for periods more than three to five years into the future. This requires a different regulatory framework to ensure that associated risks remain within accepted boundaries. Previously, single-scenario operator predictions were checked against field measurements. When subsidence or seismicity started to deviate beyond claimed uncertainties, the operator was asked to provide prediction updates. The practice was long considered acceptable, as structural damage to buildings and infrastructure or personal risk to people was not expected. This all changed following the 2012 Huizinge seismic event, necessitating better identification, assessment and ranking of risks, the use of scenarios, probabilistic forecasting and a much intensified field monitoring and control loop. It requires that the regulator becomes actively involved in assuring the integrated control loop of risk identification, predictions, monitoring, updating, mitigation measures and the closing of knowledge gaps, to ensure that subsidence (rate) and induced seismicity remain within acceptable limits. And it requires that this increased involvement of the regulator is supported in the mining law and by appropriate conditions in the Production Plan assent.

Due to their potential as a petroleum or geothermal system, the Dinantian carbonates of the Netherlands have recently attracted renewed interest because of the identified presence of excellent reservoir properties. This notion contrasts with the general assumption that these carbonates are tight. Therefore, in order to give the current knowledge state, this paper re-examines the sparse publicly available well and seismic data and literature to assess the distribution and reservoir properties of the Dinantian carbonates.

Dinantian carbonate deposition occurred throughout the study area (southern onshore and offshore of the Netherlands and northern Belgium), which is situated on the northern margin of the London–Brabant Massif, progressively onlapping the latter structure. This study confirms the presence of three carbonate facies types in the study area: a Tournaisian low-gradient carbonate ramp system, succeeded by a succession in which the carbonate ramp system evolved to a rimmed shelf setting. Subsidence of the northern margin of the London–Brabant Massif resulted in a landward shift of the shallow-marine facies belts, while the formation of normal faults resulted in a ‘staircase’-shaped shallow-water platform–slope–basin profile, associated with large-scale resedimentation processes. After deposition, the limestone deposits were frequently exhumed and reburied. A first period of regional exhumation occurred at the end of the Dinantian, which seems to be associated with porosity-enhancing meteoric karstification, possibly limited to the palaeo-shelf edge. The most intense alterations seem to be present as a deep leached horizon below the Cretaceous unconformity at the top of the Dinantian sequences. In addition, clear evidence for hydrothermal fluid migration is found locally, enhancing reservoir properties at some places while occluding porosity at others. The timing of these phases of hydrothermal fluid circulation is poorly understood.

Whereas in the United Kingdom hydrocarbons have been produced from karstified Dinantian carbonates, this petroleum play has received little attention in the Netherlands. This paper shows that, also for the Netherlands, a karstic reservoir probably existed before the start of hydrocarbon generation from the onlapping basal Namurian shales. The hydrocarbon prospectivity of these sediments, however, is primarily controlled by the presence of both a karst-related reservoir and migration routes from a decent-quality source rock. Two geothermal projects producing from this reservoir in the southern onshore Netherlands have shown the potential of the Dinantian carbonates for ultra-deep geothermal projects. To conclude, the findings presented herein are relevant for studies of the hydrocarbon prospectivity and studies of the geothermal potential of Dinantian carbonates in the Dutch on- and offshore.

Depletion of the Groningen gas field has induced earthquakes, although the north of the Netherlands is a tectonically inactive region. Increased seismic activity raised public concern which led the government to initiate a number of studies with the aim of understanding the cause(s) of the earthquakes. If the relationship between production and seismicity were understood then production could be optimized in such a way that the risk of induced seismicity would be minimal. The main question remains how production is correlated with induced seismicity. The Minister of Economic Affairs of the Netherlands decided to reduce production starting from 17 January 2014, specifically in the centre of the gas field as it has the highest rates of seismicity, the largest-magnitude events and the highest compaction values of the field.

A reduction in production could possibly lead to a reduced rate of compaction. Additionally a reduction of production rate could lead to a reduced stress rate increase on the existing faults and consequently fewer seismic events per year. One might envisage a ‘bonus effect’ in the events reduction in the sense that the total number of events will be less, with the same total production smeared out over a longer period. This is as yet unclear.

In this paper we apply different statistical methods to look for evidence supporting or disproving a decrease in the number of seismic events due to production reduction.

High-resolution Early Holocene palynological records from the middle Meuse River valley were missing until recently. In order to investigate environmental and inferred climate changes during the Preboreal, sediments from a former residual channel of the Meuse River near Haelen were studied. Detailed multi-proxy analyses, including microfossils, macroremains and loss-on-ignition measurements, were carried out at a high temporal resolution. An accurate chronology of the >1000-year-long record was provided by accelerator mass spectrometry (AMS) 14C wiggle-match dating.

The channel was abandoned during the late Younger Dryas, when accumulation started with gyttja. This period was characterised by an open landscape with herbaceous vegetation and dwarf shrubs. Patches of birch were present on the floodplains around depressions and (oxbow) lakes. Some pines survived the cold in sheltered locations. In the residual channel the water was flowing temporarily and aquatic plant communities developed with predominantly submerged taxa and algae. The shores were fringed by willows and sedges and were probably used as a watering place by large herbivores.

Following the Late-glacial/Holocene climate warming, dated in the Haelen record around 11,520 cal BP, birch woodlands expanded on the river floodplains and slopes of terraces during the Friesland Phase. Open vegetation with herbs and juniper remained present on the nearby terraces. An increase in the water level of the oxbow lake and seepage of groundwater occurred. Along the shores herbaceous vegetation was present. Around 11,420 cal BP, birch expansion was interrupted by the dry continental Rammelbeek Phase. On the river floodplain and terrace slopes, open grassland vegetation developed and on the terraces, grasslands and open grounds were abundant. In the residual channel the water became stagnant and floating-leaved vegetation developed. At the start of the Late Preboreal, around 11,270 cal BP, a sudden shift to a more humid climate took place and birch forests expanded again on the river floodplains and terrace slopes. Poplar became more abundant in these forests, and birch and poplar swamp forests were present near the site. Pine expanded at c. 11,160 cal BP on the higher sandy and gravelly terraces. During the Late Preboreal a reed swamp developed on the shores of the residual channel.

At the onset of the Boreal, around 10,710 cal BP, woodlands, initially with hazel, but later also with oak, elm and lime, started to develop, while pine forest remained present on the higher terraces. Hazel shrubs were growing on the terrace slopes. Birch and poplar forests occurred on moist parts of the floodplains. Around the residual channel they formed a zone behind the reed swamps surrounding the oxbow lake. Vegetation with water lilies was present in open water.

The Haelen record shows, despite a lack of archaeological evidence, indications for the presence of Mesolithic people in the area during the Preboreal. These include the occurrence of (natural or man-made) fires, in combination with the presence of trampled areas and disturbed grounds and possibly consumption of Nymphaeaceae seeds and tubers.

We describe a tetrapod swimming traceway from the Middle Triassic Vossenveld Formation of the Netherlands. Forty-five individual traces, each consisting of two parallel claw drag marks, were followed over 9 m in a roughly east–west direction. The asymmetry of the traceway geometry indicates the trace maker negotiated a lateral current. The trace maker could not be identified, but the traces described here are markedly different from Dikoposichnus traces attributed to swimming nothosaurs.

We studied a Holocene peat fill of a small depression in Pleistocene coversand, in the western border zone of Het Gooi, to assess the early local and regional vegetation history in relation to sea-level rise, soil development and potential human impact. In the fourth millenium BC, a podzol which had formed in the depression became stagnative, leading to the development of a moorland pool. The local vegetation changed from dry heathland, through an amphibic vegetation type with, among others, Littorella uniflora and Lycopodiella inundata, to a permanently moist Sphagnum-dominated vegetation. The existence of moorland pools and the development of such habitats into Sphagnum-dominated vegetation are known from Late-Holocene anthropogenic, more or less open landscapes that were formed on a podzolising sandy soil under ericaceous vegetation. However, the recorded vegetation succession did not show any recognisable local human impact and therefore is attributed to natural succession. In the period concerned, sea level was still about 4 m below the land surface in the depression, implying that water logging occurred independent from a rise in sea level and associated groundwater level. It took until the Late Middle Ages before such rise led to significant water logging and peat growth in this border zone, but the mean groundwater level never reached to above NAP (Dutch Ordnance Datum).