A fragmentary and partially dissociated arm of an astropectinid starfish, assigned to the genus Astropecten Gray, 1840, is recorded from decalcified and limonitised, fine-grained sandstone concretions and slabs of late Miocene (Tortonian) age, collected from the former gravel and sand pit ‘de Hoogdonk’ near Liessel, province of Noord-Brabant. Added are two examples of the ichnofossil Asteriacites lumbricalis von Schlotheim, 1820, both preserved in convex hyporelief and co-occurring with other, indeterminate traces. Although astropectinids have been recorded previously from Miocene and Pliocene deposits in the Netherlands, the present is the first record of fossilised traces of starfish.

New petrographical and geochemical evidence of tourmalinite pebbles from two Lower and three Middle Devonian formations from Belgium is presented. Petrography, REE, transition metal and microprobe data of the studied rocks suggest it are (meta)sediment-derived tourmalinites formed by boron metasomatism (and hydrothermal brecciation) in an intrusive setting. Tourmaline mineralizations within eastern Avalonia are known in south-east Ireland, the English Lake District and East Anglia. Based on previously suggested relations between Early Palaeozoic igneous activity in last two mentioned regions and the Brabant Massif, it is presumed that the Brabant Massif also underwent granite-related tourmalinization and that this was the source of the studied pebbles. Petrologie differences between pebbles found in Middle Devonian formations and pebbles from Lower Devonian formations, suggest that fluid circulation occurred on a regional scale, possibly extending into the Stavelot massif.

The decline of domestic natural gas production, increasing dependency on gas imports and lagging development of renewable energy production may pose serious challenges to the current high standards of secure energy supply in the Netherlands. This paper examines synergy between hydrocarbon- and geothermal exploitation as a means to reinforce energy security. The Roden gas field is used as an example to demonstrate potential delay of water breakthrough in the gas well and a resulting increase of recovered gas (up to 19%), by positioning of a geothermal doublet in the water leg of the gas field. The reservoir simulations show that the total increase of gas production primarily depends on the amount of aquifer support. An optimal configuration of gas- and geothermal wells is key to maximise gas recovery and strongly depends on the distribution of reservoir properties. The study also reveals that this option can still be beneficial for gas fields in a late stage of production.

Net Present Value calculations show that the added value from the geothermal doublet on total gas production could lead to an early repayment of initial investments in the geothermal project, thereby reducing the overall financial risk. If no subsidies are taken into account, the additional profits can also be used to finance the geothermal project up to break-even level within 15 years. However, this comes with a cost as the additional profits from improved gas recovery are significantly reduced.

In the conservation of fragile fossil bone material, impregnation by solvent-borne consolidant is often required. Understanding the mode of penetration of consolidants into fossil bone is of crucial importance. It is governed by a variety of factors; neutron imaging is a powerful tool to monitor and visualise this penetration (non-destructively). The consolidation of a vertebrate fossil from the Maastrichtian of the southeast Netherlands was imaged at the High Flux Reactor facility at Petten, the Netherlands. The analysis shows current conservation practice to result in a sufficiently deep and isotropic penetration.

Until now hemispherical astylospongiid sponges were invariably referred to as Caryospongia juglans var. basiplana Rauff. Renewed investigations have now shown that part of the material should be assigned to a new genus and species, Tympanospongia vankempeni, which is characterised by a system of very irregular canals. These flat-based sponges originate from the Baltic region and occur in two assemblages of silicified Late Ordovician sponges known exclusively as erratics from The Netherlands and northern Germany. These fossils were transported by the River Eridanos, a former drainage system from the Baltic region that filled the North European Basin during the Miocene to Early Pleistocene. Specimens of Tympanospongia vankempeni gen. et sp. nov. also occur in the Upper Pleistocene of Gotland, Sweden. The new sponge described herein principally differs from other genera of the Astylospongiidae found frequently in the erratic sponge assemblages by its irregular system of apochetes which ramify and anastomose commonly.

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 ‘Wadden Sea Archive of landscape evolution, climate change and settlement history’ project (WASA) focuses on the analysis of marine sediment archives from the East Frisian Wadden Sea region. It aims at understanding the formation of palaeolandscapes since the end of the last ice age. One part of the project studies the possible correlation and shift of archaeological settlement patterns, climate change and sea-level rise through time in order to derive archaeological expectancy maps. In this paper we present our findings for a quantifiable set of Stone Age sites in the area of the prehistorical Dornumer tidal basin, discussing them against the background of coastal environmental factors and the applied methodology of our modelling. To enable spatial analysis of these sites, we developed a palaeographic elevation model, which was subsequently flooded at 2000-year intervals between the Boreal and early Subboreal periods. Particular challenges are posed by the dynamics of marine transgression, the related changes in the natural environment and their spatial extent. As a result of our GIS-based approach, the model can be extended geographically and provides a basis for future research.

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.

Seismic hazard in low seismicity areas of Europe has traditionally been considered insignificant. However, in the light of the recently conducted paleoseismic studies along the Rhine Graben, a revision is required. Previously applied standard probabilistic seismic hazard assessment (PSHA) methods, using Poissonian approach for the earthquake occurrence, can now be substituted by renewal models where fault parameters such as the maximum magnitude, recurrence interval and the elapsed time since the last occurrence of a large earthquake, can be utilized. In this study, the application and the influence of the available paleoseismic data in the Lower Rhine Graben to seismic hazard analysis is demonstrated. The resulting hazard maps, when compared to the standard PSHA using Poissonian approach, indicate a more precise geographical distribution of the estimated seismic hazard levels. The influence of the paleoseismic data seem to be less important for return periods less than a 1000 years. Among the different input models, the highest values reach to 170 cm/sec2 for a 1000 year return period using a combination of Poissonian and renewal models.

Numerous archaeological investigations have been performed along the river Meuse in the Netherlands’ southeastern province of Limburg as part of the major ‘Maaswerken’ infrastructural project. To improve flood risk management and navigability, and for the purpose of gravel production and nature development, several areas of land covering a total of almost 2000 ha are being excavated to a great depth. In anticipation of this, archaeological research was performed for the purposes of recording and documenting archaeological remains in the most important areas and locations. From 1998 to 2015 the Cultural Heritage Agency of the Netherlands (Amersfoort) was in charge of the investigations, and acted as adviser to national public works agency Rijkswaterstaat.

The archaeological research connected with the Maaswerken project differed from regular, site-based investigations in terms of the landscape archaeology perspective on which it was based. The research themes and principles associated with this perspective were published in several documents, including a scientific policy plan published in 2004, and presented in further detail in area programmes and project briefs. The policy plan assigned each project area to one of five value assessment categories, based on the intactness of the landscape and the archaeological potential for addressing the research questions. In areas of high landscape intactness and great archaeological potential (category 1) the Agency selected zones to be surveyed and assessed, and for archaeological excavation. Though most of the fieldwork, including specialist analysis, was performed in these zones, other category project areas have also been the subject of archaeological fieldwork, including borehole surveys, site-oriented research and watching briefs, but on a more incidental basis. Observations were also made in the river Meuse itself and in the river's winter bed.

The archaeological investigations resulted in a large number of standard reports of desk studies and fieldwork, including reports of specialist analyses. A considerable proportion of these refer to the large-scale investigations at Borgharen and Itteren to the north of Maastricht, and at Lomm and Well–Aijen to the north of Venlo. The results of the investigations suggest the archaeological record here is rich and varied, with a time depth of c. 11,500 years, and traces of occupation and land use ranging from the Early Mesolithic (Well–Aijen, Borgharen) to the Second World War (Lomm).

This paper reflects on almost 20 years of archaeological research in the project areas of the Maaswerken and on the principles and methods used in the field research. The common thread is the results of landscape and archaeological studies and the relationship between them. Examples are used to illustrate results that can be regarded as important from a national perspective, and in terms of archaeological heritage management.

For more than 120 years since the name Caprovis savinii was founded on the basis of a bovid horn core from the ‘Forest Bed’ at Overstrand (Norfolk, UK), only one other find has been assigned to this species. Recently, two horn cores and one horn core fragment of the species have been discovered at two sand and gravel sorting centres in the Netherlands. The bone-bearing sediments are exposed on the North Sea floor, just off the coast of Great Yarmouth (Norfolk, UK). By comparing fossil and extant taxa, a classification within the bovid tribe Antilopini is proposed. Based on the accompanying fauna from the specific dredging area and the litho- and chronostratigraphy at Overstrand, C. savinii is placed within the Early Pleistocene.

In the Netherlands, mineral extraction by means of dredging or quarrying meets with considerable societal resistance. Land surface lowering prior to large land reconstruction projects may raise fewer objections. We have calculated the potential yields of sand and gravel from land surface lowering embedded in planned building and construction projects, and in nature, farmland and recreation area development. Our primary data sets were a compilation of spatial plans for the period 1995 – 2005 and about 95,000 borehole descriptions. Even if embedded consistently, land surface lowering would contribute modestly (up to 5.4 Mio m3/a) to the filling sand provision (annual demand 45 – 50 Mio m3/a).

Further research into seismicity caused by natural gas production from the Groningen field is necessary to improve the assessment of seismic risk and develop means to control and reduce it. Research into subsurface aspects is primarily of relevance to assess the seismic hazard component in the cause-and-effect chain that governs the seismic risk. It requires a wide range of research activities that can be broadly classified as follows:

• • Increasing understanding of the physical mechanisms that govern production-induced seismicity, in particular source mechanisms, compaction behaviour, propagation of energy to the surface, and the effects of fluctuating production.

• • Reducing uncertainty by acquiring additional field data to improve statistical inference, and developing statistical methods and procedures that can cope with the non-stationary nature of the process.

• • Developing tools and techniques to improve risk management, and support operational control and policy measures under uncertainty.

An essential requirement for further research will be the possibility of developing competing theories for many aspects of the modelling chain. This requires an overall hazard and risk assessment methodology that can accommodate multiple models, and an organisational structure that facilitates the comparison of competing approaches while safeguarding their independent development. This will have to be supported by the availability of reliable data via shared databases. Finally, the scientific community should be prepared to make a major effort to translate their research results into popular scientific versions in order to keep stakeholders abreast of progressive insight into the origin, predictability and prevention of induced seismicity.

The presence of salt in dilatant normal faults may have a strong influence on fault mechanics in the Groningen field and on the related induced seismicity. At present, little is known of the structure of these fault zones. This study starts with the geological evolution of the Groningen area, where, during tectonic faulting, rock salt may have migrated downwards into dilatant faults. These fault zones therefore may contain inclusions of rock salt. Because of its rate-dependent mechanical properties, the presence of salt in a fault may introduce a loading-rate dependency into fault movement and affect the distribution of magnitudes of seismic events. We present a first-look study showing how these processes can be investigated using a combination of analogue and numerical modelling. Full scaling of the models and quantification of implications for induced seismicity in Groningen require further, more detailed studies: an understanding of fault zone structure in the Groningen field is required for improved predictions of induced seismicity. The analogue experiments are based on a simplified stratigraphy of the Groningen area, where it is generally thought that most of the Rotliegend faulting has taken place in the Jurassic, after deposition of the Zechstein. This suggests that, at the time of faulting, the sulphates were already transformed into brittle anhydrite. If these layers were sufficiently brittle to fault in a dilatant fashion, rock salt was able to flow downwards into the dilatant fractures. To test this hypothesis, we use sandbox experiments where we combine cohesive powder as analogue for brittle anhydrites and carbonates with viscous salt analogues to explore the developing fault geometry and the resulting distribution of salt in the faults. Using the observations from analogue models as input, numerical models investigate the stick-slip behaviour of fault zones containing ductile material qualitatively with the discrete element method (DEM). Results show that the DEM approach is suitable for modelling the seismicity of faults containing salt. The stick-slip motion of the fault becomes dependent on shear loading rate with a modification of the frequency–magnitude distribution of the generated seismic events.

During the 2007 ‘Kor en Bot’ collecting trip across the Oosterschelde (province of Zeeland, southwest Netherlands), on board trawler cutter ZZ10, a stem fragment of a fossil isocrinid was recognised amongst the contents of the nets pulled on deck. This specimen is here interpreted to be of Early Jurassic age and assignable to the genus Isocrinus. However, because only internodals are preserved in this pluricolumnal, specific identification cannot be but approximate (Isocrinus (Chladocrinus) cf. tuberculatus). In the absence of any outcrop of Jurassic deposits in Zeeland and adjacent Dutch and Belgian territory, the most likely explanation is that this crinoid represents erratic material transported by precursors of the present-day River Maas (Meuse). Between the Langres Plateau and Sedan (northeast France), this river cuts through several occurrences of Lower Jurassic strata from which the present isocrinid might have originated. A less likely explanation is that it stems from boulders used for coastal reinforcement or from a Roman limestone votive altarpiece put up at the temple complex for the goddess Nehalennia, formerly present at Colijnsplaat, near Domburg (Noord-Beveland, Zeeland). Transportation from either northwest France or the southern or eastern United Kingdom, where there are coastal exposures of Jurassic strata, via the North Sea, is another option which, however, is also considered less feasible in view of the good state of preservation of the crinoid.

Pleistocene glacial stages were implemented into a 3D basin and petroleum systems model of the northeastern Netherlands to address the influence of low surface temperatures and the mechanical loading of ice sheets on the subsurface. Two ice sheet thickness scenarios were used based on published data. Overall, Quaternary glacial stages have a substantial impact on the temperature and pressure distribution in the subsurface. Subsurface temperatures are significantly reduced during glacial stages, leading to lowered present-day temperatures and a low geothermal gradient in the shallow subsurface. In deeply buried sedimentary formations, pressures build up with every glacial advance resulting in overpressures at the present day. Glacial stages do not directly influence the petroleum generation of petroleum source rocks in the area, but high pressures during loading might have impacted petroleum expulsion of the early mature Coevorden Formation. Hydrocarbon accumulations in the Lower Saxony Basin were simulated to investigate the possible effects of mechanical ice loading and unloading on hydrocarbon migration. A loss of Coevorden Formation-sourced hydrocarbons to the surface was calculated in the Lower Saxony Basin during the glacial stages, indicating an influence of glacial loading on the Mesozoic petroleum system.

For the first time, geomorphology and archaeology are combined for a 165 km long stretch of the Meuse river, resulting in a geomorphogenetic map (GKM) and a series of archaeological predictive maps (AVM). The maps cover the central part the Meuse valley, located in the province of Limburg between Mook in the north and Eijsden in the south. The area consists of fluvial and aeolian landforms of the Holocene Meuse floodplain and Younger Dryas aged terraces along it, spanning a period of approximately 15,000 years of landscape genesis and human habitation. The GKM more clearly discriminates between map units of Younger Dryas and early Holocene age than in previous mappings of the Meuse valley. The AVM series provide predictive information on the location of sites for four distinct consecutive archaeological periods and four main cultural themes. The maps contribute to a better understanding of landscape processes (fluvial and aeolian geomorphology and the impact of man on river behaviour), and the possibilities for human habitation and land use in prehistoric and historic times.

A micropalaeontological analysis (foraminifera, ostracoda and bolboforma) has been undertaken for the interval between 299.7 and 201.5m of the Maaseik borehole. This interval includes marine glauconitic sands, which are generally not or not strongly decalcified. Independent from each other, the observed taxa of the different microfossil groups indicate a Late (range down to Middle) Miocene age for the interval. On the base of the biozonations of the various microfossil groups, some correlations can be presented with Late Miocene deposits of other areas in Belgium (Antwerp and Campine area), and of other places in NW Europe.

Hydraulic fracturing is a long-established method of stimulating a well to improve the inflow or outflow potential. Hydraulic fracturing is the most successful stimulation method used by the oil and gas industry, and is also used for water injection and production wells around the world, even for drinking-water wells. Hydraulic fracturing creates a crack in the earth that is then filled with a highly conductive material (proppant). This fracture has a large inflow area compared to an unstimulated wellbore and provides a high-permeability path for the fluid to flow in or out of the reservoir.

Hydraulic fracturing has a long history of being used in hot dry rock (HDR) geothermal applications since the 1980s (Murphy & Fehler, 1986). In those often very tight reservoirs, the aim is to create fracture networks that generate the reservoir flow capacity. In high-permeability formations, fracturing can potentially double the productivity of a well. In low-permeability formations, well performance can be increased by a factor of 5–10 in most cases.

In this paper, we focus on two different scenarios of geothermal stimulation. The first is for permeable, porous formations where the heat exchange happens through the perfect contact between the fluid and the porous reservoir. Stimulation may then be necessary to create a small fracture if the pressure drop near the well is too large due to insufficient reservoir permeability. The other scenario is a formation at great depth, where the formation permeability is so extremely small that very long propped fractures would be needed to obtain sufficient flow or even where the porous system does not provide sufficient heat exchange but the heat exchange has to be facilitated by an artificial or stimulated fracture network: a so-called Enhanced Geothermal System.

For porous, permeable formations we will present examples of fracture treatments that can increase the flow rate so that the economics of the project is improved. In some formations, stimulation is then a contingency in case of poorer than expected reservoir quality. A worst-case well with a large skin value of 20 can perform with stimulation like a base-case unstimulated well. In other formations, stimulation will be integral to well design in order to optimise the project performance. For those cases the Coefficient of Performance can be improved from 7 to 25 with the aid of stimulation.

In Ultra-Deep Geothermal (UDG) recovery, the targets are reservoirs below 4000 m, because industrial heat demand requires a minimum temperature of 120°C up to 250°C. For an economic business case, the rate over a period of 15 to 25 years should be from 150 to 450 m3 h−1, depending on the boundary conditions.

Shallower reservoirs in the Netherlands often show very high permeability, but at great depth the target layers could have very low permeability (Veldkamp et al., 2018). Several stimulation methods can be used, of which hydraulic fracture stimulation with water (proppantless) is the primary candidate. Other stimulation methods are propped fracturing in sandstone, acid fracturing in carbonates and thermal stimulation.

For a geological play that is attractive for UDG in the Netherlands, the most likely stimulation method is with water fracturing, because propped fracturing would require a huge amount of proppant that is very costly. Based on analogues and conceptual designs, the expected flow rate is estimated under selected boundary conditions.