In the past 10 years the mature hydrocarbon province the West Netherlands Basin has hosted rapidly expanding geothermal development. Upper Jurassic to Lower Cretaceous strata from which gas and oil had been produced since the 1950s became targets for geothermal exploitation. The extensive publicly available subsurface data including seismic surveys, several cores and logs from hundreds of hydrocarbon wells, combined with understanding of the geology after decades of hydrocarbon exploitation, facilitated the offtake of geothermal exploitation. Whilst the first geothermal projects proved the suitability of the permeable Upper Jurassic to Lower Cretaceous sandstones for geothermal heat production, they also made clear that much detail of the aquifer geology is not yet fully understood. The aquifer architecture varies significantly across the basin because of the syn-tectonic sedimentation. The graben fault blocks that contain the geothermal targets experienced a different tectonic history compared to the horst and pop-up structures that host the hydrocarbon fields from which most subsurface data are derived. Accurate prediction of the continuity and thickness of aquifers is a prerequisite for efficient geothermal well deployment that aims at increasing heat recovery while avoiding the risk of early cold-water breakthrough. The potential recoverable heat and the current challenges to enhance further expansion of heat exploitation from this basin are evident. This paper presents an overview of the current understanding and uncertainties of the aquifer geology of the Upper Jurassic to Lower Cretaceous strata and discusses new sequence-stratigraphic updates of the regional sedimentary aquifer architecture.

The Netherlands has ample geothermal resources. During the last decade, development of these resources has picked up fast. In 2007 one geothermal system had been realised; to date (1 January 2019), 24 have been. Total geothermal heat production in 2018 was 3.7 PJ from 18 geothermal systems. The geothermal sources are located in the same reservoirs/aquifers in which the oil and gas accumulations are hosted: Cenozoic, Upper Jurassic – Lower Cretaceous, Triassic and Rotliegend reservoirs. Additionally, the yet unproven hydrocarbon play in the Lower Carboniferous (Dinantian) Limestones delivered geothermal heat in two geothermal systems. This is in contrast to the Upper Cretaceous and Upper Carboniferous with no producing geothermal systems but producing hydrocarbon fields. Similar to hydrocarbon development, developing the geothermal source relies on fluid flow through the reservoir. For geothermal application a transmissivity of 10 Dm is presently thought to be a minimum value for a standard doublet system. Regional mapping of the geothermal plays, with subsequent resource mapping, by TNO discloses the areas with favourable transmissivity within play areas for geothermal development. The website www.ThermoGis.nl provides the tool to evaluate the geothermal plays on a sub-regional scale. The Dutch geothermal source and resource portfolio can be classified using geothermal play classification of, for example, Moeck (2014). An appropriate adjective for play classification for the Dutch situation would be the predominant permeability type: matrix, karst, fracture or fault permeability. The Dutch geothermal play is a matrix-permeability dominated ‘Hot Sedimentary Aquifer’, ‘Hydrothermal’ or ‘Intra-cratonic Conductive’ play. The Dutch ‘Hot Sedimentary Aquifer’ play is subdivided according to the lithostratigraphical annotation of the reservoir. The main geothermal plays are the Delft Sandstone and Slochteren Sandstone plays.

The Danish onshore subsurface contains very large geothermal resources that have the potential to make a significant contribution to transforming Danish energy consumption toward a more sustainable energy mix. Presently, only a minor fraction of this green energy is exploited in three small plants. The main factors that have hampered and delayed larger-scale deployment are related to uncertainties in the geological models, which inevitably lead to high economic risks that are difficult for smaller district heating companies to mitigate without support from a compensation scheme. To facilitate and stimulate much wider use of the Danish geothermal resources, the Geological Survey of Denmark and Greenland (GEUS) and other research institutes have conducted several regional research projects focusing on the geological and geochemical obstacles with the principal objective of reducing the exploration risks by selecting the best geological reservoirs.

One of the most important geological factors causing uncertainty is the quality of the reservoirs and their ability to produce the expected volume of warm geothermal brine. Thus, great emphasis has been placed on investigating and understanding the relationships between reservoir sandstone, porosity, permeability, petrography, diagenetic processes and alterations related to variable sediment sources, basin entry points, depositional systems and climate, burial and thermal history. Mesozoic sandstones comprise the most important geothermal reservoirs in Denmark. Details concerning the reservoir quality are compiled and compared for the Lower Triassic Bunter Sandstone, Triassic Skagerrak, Upper Triassic – Lower Jurassic Gassum and Middle Jurassic Haldager Sand formations. The Bunter Sandstone Formation contains extensive aeolian and more confined fluvial sandstones with high porosity and permeability. However, highly saline formation water could be unfavourable. The Skagerrak Formation comprises well-sorted braided stream sandstones in the centre of the basin, and is otherwise characterised by muddy sandstones and alluvial fan conglomerates. An immature mineralogical composition has caused intensive diagenetic changes in the deepest buried parts of the basin. The Gassum Formation consists of shoreface, fluvial and estuarine sandstones interbedded with marine and lacustrine mudstones. In the upper part of the formation, the sandstone beds pinch out into mudstones towards the basin centre. Pervasive siderite- and calcite cement occurs locally in shallowly buried sandstones, and with burial depth the maximum abundances of quartz and ankerite cement increase. Sandstones of shallow burial represent excellent reservoirs. The relatively coarse grain size of the Haldager Sand Formation results in high porosity and permeability even at deep burial, so the formation comprises a high-quality geothermal reservoir.

Substantial progress has been made, and a well-established regional geological model combined with reservoir quality is now available for areas with cored wells. This has enabled an improved estimation of reservoir quality between wells for exploration of geothermal reservoirs.

Over the past decade in the Netherlands, most operators have only developed a single doublet. The learning effect from these single events is suboptimal, and operators have only been capable of developing doublets in areas with relatively low exploration risk. This ‘stand-alone’ approach can be significantly improved by a collective approach to derisk regions with similar subsurface characteristics. Such a play-based portfolio approach, which is common in the oil and gas industry, can help to accelerate the development of the geothermal industry through unlocking resource potential in areas marked by high upfront geological risk, effectively helping reduce costs for the development. The basis of the methodology is to deploy new information to the play portfolio by trading off with the risk of the first wells, resulting in a strong geological risk reduction.

The added value of the portfolio approach is demonstrated for the Netherlands in this paper through a comparison with a ‘stand-alone’ development. In the stand-alone approach, each new project will be equally risky, and therefore relatively unprofitable. In the case of a portfolio approach, all experience about the play is used optimally for derisking. In case of success, subsequent projects will have a higher chance of being successful, due to the experience gained in previous projects. Even if a project fails, this may help in increasing the probability of success for subsequent projects. For plays that are initially considered too risky for the market to start developing, the value of information (VoI) of a play-based portfolio approach will help by derisking the play to such an extent that it becomes attractive for the market to develop, even at high initial risk. It can be demonstrated for several geothermal plays in the Netherlands that by adopting the portfolio approach, the probability of a play being developed becomes higher, the number of successfully developed projects increases and the average profitability of the project will also be higher. Five more advantages are: (1) continuous improvement by integrated project development, (2) cost reduction through synergy, efficiency and standardisation, (3) optimisation of the surface heat demand and infrastructure, (4) the possibility of structural research and development (R&D) and innovation, and (5) financing advantages. The advantages reinforce each other.

A preliminary estimate of the geothermal potential of the Netherlands adopting the portfolio approach is between 90 and 275 Petajoules (PJ). For about 350 doublets being developed, producing about 70 PJ, the value of the advantage of the play-based portfolio approach is €2 billion for the three main plays: Rotliegend, Triassic and Jurassic/Cretaceous. The learning effects of synergy, efficiency and standardisation are expected to be significant.

The main objective of this paper is to give an overview of the risks seen in the exploration and production of geothermal energy from the viewpoint of the regulator. The risks are categorised as conventional risks, ultra-deep risks and enhancing factors. These risks are similar to those seen in the oil and gas industry, but the maturity of the geothermal sector in terms of managing such risks is much lower.

Another objective of this paper is to discuss how these risks are managed and mitigated by the sector and the supervisor, State Supervision of Mines (SodM). Portfolio operators developing multiple projects, using skilled employees and embracing continuous improvement are seen as the way forward for the sector to grow safely and sustainably.

This paper concludes that positive developments have started, but a lot of work still needs to be done to ensure safe growth of the geothermal energy sector.

The area ‘Het Gooi’ in the Netherlands is part of a Pleistocene ice-pushed ridge system that partially drowned during the Holocene upon sea level and associated groundwater rise. As a result, the ridge system was gradually encroached by peat. From the late Middle Ages onward, man reclaimed the peatlands surrounding Het Gooi, heavily reducing their extension and lowering the regional groundwater level by increasingly intensive drainage. Based on historical and archaeological arguments, several authors assume that the Holocene peat cover in the border zones of ‘Het Gooi’ formed the extension of large raised peat bogs that formed further to the west and east, respectively. They presume that in the late Middle Ages these extensions reached ‘upslope’ to a maximum altitude of 3 m + NAP (Dutch Ordnance Datum – approximating mean sea level). However, the original extension is difficult to reconstruct, as this peat has disappeared as a result of its exploitation and oxidation, if having been present at all.

In this study, the maximum extension of the Holocene peat cover on the ice-pushed ridge system was reconstructed based on soil characteristics. Used soil characteristics concerned the presence of iron coatings around sand grains and the upper boundary of gleyic features, because these are indicators for the mean highest groundwater level (MHG). For peat to form, this MHG needs to be at or just above the ground surface for most of the year. Based on study of a number of soil transects, we reconstructed to what maximum altitude peat encroachment may have occurred. This ‘maximum extension’ can alternatively be described as the maximum altitude of the bottom of the peat onlapping the ridge system.

In the western border zone, this peat cover was found to have reached to c. NAP or just above, near Hilversum. No indications were found for the occurrence of raised bogs. We conclude that the phreatic groundwater level in this zone was controlled by the sea level and associated lake levels (Naardermeer and Horstermeer), a dominant role being played by the shallow presence of Pleistocene formations with a high hydraulic conductivity. In the eastern border zone, altitudes were more variable and in places reached 2 m + NAP. Peat at this higher elevation probably formed under the influence of a higher phreatic groundwater level, induced by the presence of a clayey Eemian fill with low hydraulic conductivity in the adjacent glacial basin (the Eem valley).

This study demonstrates the value of detailed soil transect studies for palaeogeographical reconstructions of the former Holocene peat cover in Pleistocene landscapes of NW Europe. It also provides independent data for validation of geohydrological models for such landscapes.

Recently drilled geothermal boreholes in Mol, northern Belgium, provide new information on the Lower Carboniferous carbonates in the Campine–Brabant Basin. Because of low primary porosity, fractures in these limestones and dolostones are of major importance for reservoir permeability. The Fullbore Formation MicroImager (FMI) log of the MOL-GT-01 borehole enabled interpretation of bed boundaries and fractures in the subsurface. Relationships between the frequency of these fractures and bed thickness, lithology variations and the presence of faults were explored. The results show that thick beds contain relatively few fractures and thin beds relatively many. Except for lower values in shaly intervals, the fracture frequency (number per metre) is largely independent of lithology. Zones with substantial changes in the structural dip (called a cusp) and/or azimuth of bed boundaries were identified. The clearest cusp is present at a depth of 3284 m. Since the presence of a normal fault is most likely regarding the local and regional geology, this cusp likely resembles a fault-tip fold of a WSW-ward dipping normal fault with an inclination of at least 54°. It is uncertain whether the borehole crossed the fault itself or only a monocline on top of it. Fracture frequency is increased in the vicinity of the interpreted possible faults. Up to a vertical distance of c.45 m from the faults, the mean fracture frequency is higher than in a non-faulted zone with similar lithology. However, frequency differences between these faulted and non-faulted zones are mostly insignificant, so no clear damage zones are present.

Geological, geomorphological and soil maps provide important information on the substrate as well as on the past and present physical landscape. For the intensely studied Netherlands coastal plain and Rhine–Meuse delta, many such map datasets have been compiled over the last two centuries. These mapping materials comprise older and younger legacy datasets, often fragmented over regions. They have been compiled within various research traditions and by various parties, involving geologists, soil scientists, geomorphologists and landscape archaeologists. The maps and datasets summarise overwhelming amounts of underlying data accumulated over the last few centuries, and are therefore valuable for reconstructing past landscapes.

Digital-infrastructure developments have enhanced possibilities for recombining existing and new data over the last few decades, e.g. through GIS solutions such as palaeogeographical base maps, from which multiple derived map products can be generated. Integration of thematic information from various source maps and underlying data is needed to use the accumulated data diversity to its full potential and to answer applied and fundamental scientific questions. Using diverse information to compile or update maps, however, requires awareness of legacy surveying strategies and the state of knowledge at the time the original data and maps were produced. This paper reviews the soil, geological and geomorphological mapping traditions. We evaluate their products, underlying data and the reasoning behind their compilation, focusing on their use in conventional and digital palaeogeographical mapping. This helps get the most out of large quantities of legacy and modern data, a major challenge for surface and substrate digital mapping in the big-data era.

This study presents an annually resolved dendrochronological reconstruction of hydrological impacts on the Roman and early-medieval landscape in the Low Countries of northwestern Europe. Around 600 hydrologically sensitive ring-width patterns, mostly oak (Quercus robur/petraea) as well as some ash (Fraxinus excelsior) and elm (Ulmus sp.), were selected from an initial dataset of >5000 and compiled into two chronologies that span the first millennium AD. Their content and (dis)similarities to established tree-ring chronologies from this and surrounding regions were used to assess their provenance, which in both cases is in the area where the majority of the wood was recovered. Instances of high groundwater levels and/or inundation were catalogued by identifying multi-year intervals of strongly reduced annual growth that occurred simultaneously throughout the research area. The resulting record contains 164 events dated between AD 1 and 1000, of which 21 have a recurrence frequency ≥50 years. One-third of the ≥50-yr events date between AD 185 and 282, making this the most flood-intense interval of the first millennium. The severest reconstructed impact of the first millennium dates to AD 602. A comparison to historically documented river floods/sea breaches and drought/heat spells shows that the predominant cause of the inferred impacts in the research area was river overflow. Synchronous inundation responses of oaks preserved in former bogs in Lower Saxony (NW Germany) indicate that half of the reconstructed events occurred on a supra-regional level, pointing to regional precipitation as a main forcing. River floods documented in written sources do not seem to have affected tree growth in Lower Saxony in a significant manner, indicating that the majority of documented floods most likely were caused by hydrological circumstances upstream of the catchments of the Rhine and/or Meuse. Reconstructed flood impacts during the Early Middle Ages coincide remarkably well with construction and repair of Rhine revetments at the early-medieval site of Leiderdorp-Plantage in the western Netherlands.

The results from analyses of botanical remains (pollen, wood, charcoal, seeds) from several archaeological features excavated in Kluizen (northern Belgium) are presented. The region was largely uninhabited until the Iron Age and Roman period when a rural settlement was established, resulting in small-scale woodland clearance. The site was subsequently abandoned from c. AD 270 till the High Middle Ages. The results of the archaeological and archaeobotanical analyses provide information on changes in land use and resulting dynamics of woodland cover and composition between c.600 BC and AD 1200, with a spatial and temporal resolution unrivalled in northern Belgium. Especially the long period of woodland regeneration following abandonment of the site around AD 270, covering the Late Roman and Early Medieval period, could be reconstructed in detail. Abandoned fields were first covered with pioneer woodland (Salix, Corylus and Betula), then Quercus-dominated secondary forest and finally a late-successional forest with Fagus sylvatica, Carpinus betulus and Ilex aquifolium, an evolution that took over 300 years. The results also indicate that the observed increase of Fagus during the Early Middle Ages, which was never an important element in the woodland vegetation in northern Belgium before, was related to climatic changes rather than anthropogenic factors.

The Younger Dryas cold period caused major changes in vegetation and depositional environments. This study focuses on the aeolian river-connected dunes along the former, Weichselian Late Glacial, course of the Scheldt River in the southern Netherlands. Aeolian dunes along the Scheldt have received little attention, as they are partly covered by Holocene peat and marine deposits. The spatial distribution of the dunes is reconstructed by digital elevation model analysis and coring transects. Dunes are present on the high eastern bank of the Scheldt and in the subsurface of the polder area west of the Brabantse Wal escarpment. A reach-specific higher channel gradient probably caused a channel pattern change from meandering to braiding during the Younger Dryas. This enabled deflation from the braid plain and accumulation in source-bordering river dunes east of the incised and terraced, subsurface Late Pleistocene Scheldt valley. The age of the dune formation is established by pollen analysis and radiocarbon dating of underlying and overlying peat beds. The peat layer below the dune at Zomerbaan is attributed to the Allerød and early Younger Dryas periods. Dune formation occurred predominantly during the second part of the Younger Dryas stadial, both on and in front (west) of the Brabantse Wal escarpment. Wind direction was reconstructed by geomorphic analysis and sedimentary structures on lacquer peels. A southwesterly wind direction is demonstrated by the parabolic dune morphology. For the first time, Younger Dryas wind direction is reconstructed based on adhesion ripple cross-laminated sets on lacquer peels. Sand-transporting south-southwesterly winds were dominant during the Younger Dryas, most likely during summer.

This paper focuses on unravelling the 1st millennium AD in the present-day Netherlands and the applicability of modelling when studying the past. By presenting the results of several studies analysing changes (or persistence) in connectivity and habitation patterns, the significance of these findings for (spatial) modelling is derived. The transition between the Roman and early-medieval periods is particularly interesting in this respect as it is characterised by severe pan-European political, socio-economic and demographic changes. Additionally, recent studies in geosciences increasingly point to marked climatic and landscape changes, such as river avulsions and floods, occurring at the same time. The extent to which these environmental and cultural dynamics were entwined and mutually influential is generally unknown, especially on larger-scale levels. Lowlands, such as the Netherlands, are especially suited to study these complex interactions since boundary conditions, i.e. the set of conditions required for maintaining the existing equilibrium in a region, in such areas are particularly sensitive to change.

In this paper the combined results of several recently developed landscape-archaeological models are presented. These models spatially analyse natural and cultural dynamics in five manifestations: route networks, long-distance transport, settlement patterns, palaeodemographics and land-use systems. Combined, these manifestations provide information on connectivity, persistence and habitation, key concepts for the cultural landscape as a whole. Results show that only by integrating these modelling outcomes is it possible to reconstruct boundary conditions and high-resolution spatio-temporal frameworks for cultural-landscape change. Equally, these models invite reflection on their applicability and, as such, point to the need for new theoretical framing and the development of more multi-proxy, evidence-based and transdisciplinary research approaches in archaeology. The evident interrelationship between cultural and natural-landscape dynamics necessitates a more integrated and transparent research attitude, covering multiple scales and studying the cultural landscape as a whole. Only then can models reflect historical reality as closely as possible.

The 1717 Christmas flood is one of the most catastrophic storm surges the Frisian coast (Netherlands and Germany) has ever experienced. With more than 13,700 casualties it is the last severe storm surge with a death toll of this order. At the same time, little is known about the hydrodynamic conditions and the morphological effects associated with this storm surge.

In this study, 41 potential dyke failures in the Province of Groningen (Netherlands) associated with the 1717 Christmas flood were systematically reconstructed and mapped by using historical maps and literature and by analysing the recent topography in search of typical pothole structures and sediment fans. The dimensions of the sediment fans as derived from the topography show a good accordance with the dimensions documented by vibracore profiles, direct push tests and electrical resistivity tomography data taken at three fieldwork sites. Moreover, the fan dimensions closely agree with the dimensions as simulated using a process-based morphodynamic numerical model for one of the three sites, the village of Wierhuizen. Consequently, the recent topography is still indicative for the locations and dimensions of dyke failures and sediment fans associated with the 1717 Christmas flood. Considering the large number of detected dyke failures (41) and the large dimensions of the potholes and particularly of the sediment fans up to a few hundred metres wide and up to 0.7 m thick, this study proves significant morphological effects of the 1717 Christmas flood on the mainland of the Province of Groningen.

Based on the numerical simulation approach and the comparison with field data and field observations, a maximum seaward water level of 5 m NAP for the dyke failure at Wierhuizen during the Christmas flood can be derived. A similar maximum water level is indicated for the two other fieldwork sites Vierhuizen and Kohol, which is in good agreement with the maximum storm surge level of 4.62 m NAP historically documented for the city of Emden located almost 50 km to the east of Wierhuizen.

The results of the current study demonstrate that the reconstruction of historical dyke failures based on (i) historical sources, (ii) recent lidar/high-resolution topographical data, (iii) multi-proxy sedimentary field data and (iv) hydro- and morphodynamic numerical simulations is a highly promising approach to derive hydrodynamic conditions and the morphological onshore response of the 1717 Christmas flood in the Province of Groningen. This knowledge is essential to improve our understanding of extreme storm surge dynamics, their influence on the coastal landscape and the associated hazards for the coastal population.