An overview of induced seismicity in the Netherlands

We present an overview of induced seismicity due to subsurface engineering in the Netherlands. Our overview includes events induced by gas extraction, underground gas storage, geothermal heat extraction, salt solution mining and post-mining water ingress. Compared to natural seismicity, induced events are usually small (magnitudes ≤ 4.0). However, due to the soft topsoils in combination with shallow hypocentres, in the Netherlands events exceeding magnitude 1.5–2.0 may be felt by the public. These events can potentially damage houses and infrastructure, and undermine public acceptance. Felt events were induced by gas production in the north of the Netherlands and by post-mining water ingress in the south-east. Notorious examples are the earthquakes induced by gas production from the large Groningen gas field with magnitudes up to 3.6. Here, extensive non-structural damage incurred and public support was revoked. As a consequence, production will be terminated in 2022 leaving approximately 800 billion cubic metres of gas unexploited. The magnitudes of the events observed at underground gas storage, geothermal heat production and salt solution mining projects have so far been very limited (magnitudes ≤ 1.7). However, in the future larger events cannot be excluded. Project- or industry-specific risk governance protocols, extensive gathering of subsurface data and adequate seismic monitoring are therefore essential to allow sustainable use of the Dutch subsurface now and over the decades to come.

Palynological indications for Silurian – earliest Devonian age strata in the Netherlands

Knowledge of the stratigraphic development of pre-Carboniferous strata in the subsurface of the Netherlands is very limited, leaving the lithostratigraphic nomenclature for this time interval informal. In two wells from the southwestern Netherlands, Silurian strata have repeatedly been reported, suggesting that these are the oldest ever recovered in the Netherlands. The hypothesised presence of Silurian-aged strata has not been tested by biostratigraphic analysis. A similar lack of biostratigraphic control applies to the overlying Devonian succession. We present the results of a palynological study of core material from wells KTG-01 and S05-01. Relatively low-diversity and poorly preserved miospore associations were recorded. These, nonetheless, provide new insights into the regional stratigraphic development of the pre-Carboniferous of the SW Netherlands. The lower two cores from well KTG-01 are of a late Silurian (Ludlow–Pridoli Epoch) to earliest Devonian (Lochkovian) age, confirming that these are the oldest sedimentary strata ever recovered in the Netherlands. The results from the upper cored section from the pre-Carboniferous succession in well KTG-01 and the cored sections from the pre-Carboniferous succession in well S05-01 are more ambiguous. This inferred Devonian succession is, in the current informal lithostratigraphy of the Netherlands, assigned to the Banjaard group and its subordinate Bollen Claystone formation, of presumed Frasnian (i.e. early Late Devonian) age. Age-indicative Middle to Late Devonian palynomorphs were, however, not recorded, and the overall character of the poorly preserved palynological associations in wells KTG-01 and S05-01 may also suggest an Early Devonian age. In terms of lithofacies, however, the cores in well S05-01 can be correlated to the upper Frasnian – lower Famennian Falisolle Formation in the Campine Basin in Belgium. Hence, it remains plausible that an unconformity separates Silurian to Lower Devonian strata from Upper Devonian (Frasnian–Famennian) strata in the SW Netherlands. In general, the abundance of miospore associations points to the presence of a vegetated hinterland and a relatively proximal yet relatively deep marine setting during late Silurian and Early Devonian times. This differs markedly from the open marine depositional settings reported from the Brabant Massif area to the south in present-day Belgium, suggesting a sediment source to the north. The episodic presence of reworked (marine) acritarchs of Ordovician age suggests the influx of sedimentary material from uplifted elements on the present-day Brabant Massif to the south, possibly in relation to the activation of a Brabant Arch system.

The impact of natural fractures on heat extraction from tight Triassic sandstones in the West Netherlands Basin: a case study combining well, seismic and numerical data

The positive impact that natural fractures can have on geothermal heat production from low-permeability reservoirs has become increasingly recognised and proven by subsurface case studies. In this study, we assess the potential impact of natural fractures on heat extraction from the tight Lower Buntsandstein Subgroup targeted by the recently drilled NLW-GT-01 well (West Netherlands Basin (WNB)). We integrate: (1) reservoir property characterisation using petrophysical analysis and geostatistical inversion, (2) image-log and core interpretation, (3) large-scale seismic fault extraction and characterisation, (4) Discrete Fracture Network (DFN) modelling and permeability upscaling, and (5) fluid-flow and temperature modelling. First, the results of the petrophysical analysis and geostatistical inversion indicate that the Volpriehausen has almost no intrinsic porosity or permeability in the rock volume surrounding the NLW-GT-01 well. The Detfurth and Hardegsen sandstones show better reservoir properties. Second, the image-log interpretation shows predominately NW–SE-orientated fractures, which are hydraulically conductive and show log-normal and negative-power-law behaviour for their length and aperture, respectively. Third, the faults extracted from the seismic data have four different orientations: NW–SE, N–S, NE–SW and E–W, with faults in proximity to the NLW-GT-01 having a similar strike to the observed fractures. Fourth, inspection of the reservoir-scale 2D DFNs, upscaled permeability models and fluid-flow/temperature simulations indicates that these potentially open natural fractures significantly enhance the effective permeability and heat production of the normally tight reservoir volume. However, our modelling results also show that when the natural fractures are closed, production values are negligible. Furthermore, because active well tests were not performed prior to the abandonment of the Triassic formations targeted by the NLW-GT-01, no conclusive data exist on whether the observed natural fractures are connected and hydraulically conductive under subsurface conditions. Therefore, based on the presented findings and remaining uncertainties, we propose that measures which can test the potential of fracture-enhanced permeability under subsurface conditions should become standard procedure in projects targeting deep and potentially fractured geothermal reservoirs.

Luminescence dating of a late Middle Pleistocene glacial advance in eastern England

Previous investigation of isolated landforms, on the eastern margin of the East Anglian Fenland, England, has demonstrated that they represent an ice-marginal delta and alluvial fan complex deposited at the margin of an ice lobe that entered the Fenland during the ‘Tottenhill glaciation’ (termed the ‘Skertchly Line’). They have been attributed, based on regional correlations, to a glaciation during the Late Wolstonian (i.e. Late Saalian) Substage (Drenthe Stadial, early Marine Isotope Stage (MIS) 6). This paper aimed to test this correlation by directly optically luminescence dating, for the first time, sediments found within the Skertchly Line at Shouldham Thorpe, Norfolk, and Maidscross Hill, Suffolk, together with those in associated kame terrace deposits at Watlington, Norfolk. Ages ranged from 244 ± 10 ka to 12.8 ± 0.46 ka, all the results being younger than MIS 8 with some clearly showing the landforms have been subsequently subjected to periglacial processes, particularly during the Late Devensian Substage (∼MIS 2). Most of the remainder fall within the range 169–212 ka and could be assigned to MIS 6, thus confirming the previously proposed age of the glaciation. The local and regional implications of these conclusions are discussed, the maximum ice limit being linked to that of the Amersfoort–Nijmegen glaciotectonic ridge limit in the central Netherlands.

Holocene relative sea-level data for the East Frisian barrier coast, NW Germany, southern North Sea

Collecting sea-level data from restricted coastal areas is essential for understanding local effects on relative sea level. Here, a revised relative mean sea-level curve for the area of the East Frisian island Langeoog, northwestern Germany, for the time period from 7200 cal BP until Recent is presented. The revision is based on the reinterpretation of previously published and unpublished data following the HOLSEA standardisation of data handling. Altogether 68 sea-level data taken from 32 cores and outcrops from Langeoog, its back-barrier and the adjacent mainland, which have been collected since the 1950s for mapping and landscape reconstruction purposes, are presented. The age constraints, derived from radiocarbon ages of basal peat, intercalated peat and molluscs and optical dating of tidal deposits, were evaluated in terms of the HOLSEA sea-level protocol and their stratigraphic context. For 7200 cal BP until modern times, 30 sea-level index points with different uncertainty ranges were defined. Additionally, a factor of decompaction was estimated for the remaining basal peat samples as well as for the underlying sediments of intercalated peat samples. The comparison of the Langeoog relative sea-level curve with the relative sea-level curve from the western Netherlands shows that the Langeoog curve lies up to 0.80 m lower than the Dutch curve and diverges for the time before 6000 cal BP. Though the offset coincides with the overall predicted trend of glacial-isostatic adjustment, it is less than predicted. Our study provides a useful assessment of legacy data and contributes to an improved sea-level index dataset for the southern North Sea coast.

Landscape changes and human–landscape interaction during the first millennium AD in the Netherlands

The first millennium AD encompasses the Roman period (12 BC to AD 450) and the Early Middle Ages (AD 450 to 1050). In the Netherlands, this millennium saw population growth, steep decline and subsequent revival. In addition, many changes occurred in the physical landscape, marking a transition from a mainly natural prehistorical lowland landscape to an increasingly human-affected landscape. This paper synthesises the main landscape changes and human–landscape interactions in the Netherlands during this dynamic period. The degree of landscape change is compared between the coastal plain, the delta and the Pleistocene sand area. Human activities caused major often unintended geomorphological changes in all studied landscapes. Landscape sensitivity to human impact, however, strongly varied as a function of preceding landscape state. The most prominent changes took place in many parts of the coastal plain, where back-barrier peatlands transformed into open tidal basins. Presence of compaction-prone peat and intensified artificial drainage led to subsidence. This precondition and slow anthropogenic forcing combined, made the area more sensitive to stochastically occurring storms, which could serve as tipping points for large-scale drowning. Eventually, major peatlands turned into tidal areas that for many centuries would remain unsuitable for habitation. Human-induced peatland subsidence also led to the formation of the new Hollandse IJssel and Lek river branches. This marked a major reorganisation of the river network in the lower Rhine–Meuse delta. In the middle and upstream parts of the delta, the landscape was more stable. Yet, settlements on the natural levees show adaption to increasing flooding frequency from the Late Roman period onwards. The settlements shifted towards higher positions, while route networks between them largely remained intact. Smaller-scale landscape changes were found in the Pleistocene sand area. Here, local sand drifting occurred, most frequently occurring close to human movement corridors. Drift sand intensity became larger as population density increased after c.AD 900.

From dust till drowned: the Holocene landscape development at Norderney, East Frisian Islands

Within the multidisciplinary WASA project, 160 cores up to 5 m long have been obtained from the back-barrier area and off the coast of the East Frisian island of Norderney. Thirty-seven contained basal peats on top of Pleistocene sands of the former Geest and 10 of them also had intercalated peats. Based on 100 acclerator mass spectrometry (AMS) 14C dates and analyses of botanical as well as zoological remains from the peats, lagoonal sediments and the underlying sands, a variety of distinct habitats have been reconstructed. On the relatively steep slopes north of the present island, a swampy vegetation fringe several kilometres wide with carrs of alder (Alnus glutinosa) moved in front of the rising sea upwards of the Geest as it existed then until roughly 6 ka, when the sea level reached the current back-barrier region of Norderney at around −6 m NHN (German ordnance datum). From then on for nearly 4000 years a changing landscape with a mosaic of freshwater lakes and fens existed within this area. It was characterised by various stands of Cladium mariscus (fen sedge), alternating with brackish reed beds with Phragmites australis (common reed) and salt meadows with Aster tripolium (sea aster), Triglochin maritima (sea arrowgrass), Juncus gerardii (saltmarsh rush) as well as mudflats with Salicornia europaea (common glasswort). As far as shown by our cores, this highly diverse, and for humans potentially attractive landscape was at least some 4 km wide and followed the coast for about 10 km. Before the rising sea caused diversification of habitats, wet heath as well as dry and dusty sand areas existed. In the course of time, parts of the wet heath turned into raised Sphagnum bogs under an oceanic precipitation regime before this diverse landscape was drowned by the rising sea and finally covered by marine sediments, while the earlier sediments and peats were partly eroded and redeposited.

Facies characterisation of sediments from the East Frisian Wadden Sea (Germany): new insights from down-core scanning techniques

Sediment facies provide fundamental information to interpret palaeoenvironments, climatic variation, archaeological aspects and natural resource potentials since they are summary products of depositional processes, environmental conditions and biological activities for a given time and location. The conventional method of facies discrimination relies on macroscopic and/or microscopic determination of sediment structures combined with basic physical, chemical and biological information. It is a qualitative measure, depending on observer-dependent sedimentological descriptions, which cannot be reanalysed readily by further studies. Quantitative laboratory measurements can overcome this disadvantage, but are in need of large sample numbers and/or high temporal resolution, and are time-, labour- and cost-intensive. In order to facilitate an observer-independent and efficient method of facies classification, our study evaluates the potential of combining four non-destructive down-core scanning techniques: magnetic susceptibility (MS), X-ray computed tomography (CT), X-ray fluorescence (XRF) and digital photography. These techniques were applied on selected sections of sediment cores recovered around the island of Norderney (East Frisian Wadden Sea, Germany). We process and integrate the acquired scanning measurements of XRF elemental intensities, represented by principal components, MS, CT density and lightness of eight sediment facies previously recognised by conventional facies analysis: moraine, eolian/fluvial, soil, peat, lagoonal, sand flat, channel fill and beach-foreshore. A novel type of density plot is introduced to visualise the digitised sediment information that allows an observer-independent differentiation of these facies types. Thus, the presented methodology provides the first step towards automated supervised facies classification with the potential to reproduce human assessments in a fully reproducible and quantitative manner.