During several archaeological excavations on a river terrace of the river Meuse near the village of Lomm (southeast Netherlands) information was gathered for a reconstruction of the sedimentation and vegetation history during the Holocene. Various geoarchaeological methods – geomorphological, micromorphological and botanical analyses – were applied, while accelerator mass spectrometry (AMS) 14C and optically stimulated luminescence (OSL) dating provided an accurate chronology for the sediments.

During the Early Holocene, many former braided river channels were deepened due to climate amelioration. Later, river flow concentrated in one main river channel to the west, at the location of the modern Meuse. The other channels were only active during floods, and infilling continued until the Bronze Age. Because of the higher setting of the Lomm terrace, it was only occasionally flooded and therefore formed an excellent location for habitation. Humans adapted to the changing landscape, as most remains were found on the higher river terraces or their slopes, a short distance from the Maas river. The Lomm terrace was more or less continuously inhabited from the Mesolithic onwards.

During the Early Holocene, river terraces were initially densely forested with birch and pine. From the Boreal (Mesolithic) onwards, dense mixed forests with deciduous shrubs and trees such as hazel, oak, elm and lime developed. During the Atlantic (Meso/Neolithic), the deciduous forests became dominated by oak. Due to human activities from the Late Subboreal (Late Bronze Age) onwards, forests slowly became more open, yet remained relatively dense in comparison to other Dutch areas. The botanical data, however, show that within the Lomm study area there was a large difference in the composition, distribution and openness of the vegetation. The spatial variation in openness came into existence during the Late Bronze Age, as soon as the higher areas started to be used for human activities (i.e. habitation, agriculture and livestock herding). Due to human activities, the northern part of the study area became very open during the Early Roman period. In the lower-situated areas of the southern part, however, forests remained present much longer, until the Early Middle Ages. Due to large-scale deforestation in the Lomm area and hinterland during the Roman period and Middle Ages, the sediment load of the river increased, large floods occurred and overbank sediments were deposited, burying the archaeological remains. The largest increase in sedimentation occurred after the Middle Ages.

Remains of elasmosaurid plesiosaurs are exceedingly rare in the type-Maastrichtian strata (Late Cretaceous, southeast Netherlands and northeast Belgium), in stark contrast to relatively common skeletal remains of mosasaurs. Here, we present an analysis of δ13C stable isotope values for tooth enamel of two elasmosaur teeth from the type Maastrichtian. The δ13C signal is a proxy for foraging area, trophic level and diving behaviour, the net value of which in these rare elasmosaurs turns out to be not noticeably different from that for the much commoner mosasaurs in the type Maastrichtian. Therefore, the rarity of elasmosaurs in the area probably reflects a primary near-absence of such reptiles during the latest Cretaceous, rather than a taphonomic artefact.

In the Monnikenberg estate, built in the mid-19th century in a Pleistocene coversand area to the south of Hilversum (the Netherlands), a depression with a small man-made mire was found to contain the remains of a complex peaty fill of a larger mire. Microfossil study and radiocarbon dating showed that its development started at the beginning of the third millennium BC through water stagnation on a podzol in combination with a high groundwater level, and about two millennia ago the vegetation had progressed from a mesotrophic mire to an ombrotrophic bog. Such origin, development and age are in conformity with the results for some other recently studied mires in the Het Gooi region (Gijzenveen; Van Geel et al., 2016). In the Monnikenberg mire, a discrepancy was found between the radiocarbon age of the upper part of the peaty fill and its clearly medieval or younger age as evidenced by its microfossil content, notably the occurrence of pollen of buckwheat (Fagopyrum) and cornflower (Centaurea cyanus). These species were first introduced in the Netherlands during the late Middle Ages. From further observations in a large pit in the former border zone of the mire, it is concluded that this discrepancy and the unique sedimentary structures, observed in this pit, result from peat exploitation and contemporary reworking of its sandy overburden during the Late Middle Ages or Early Modern Times. This is the first evidence for local early peat exploitation for such mires in Het Gooi. The results underpin the importance of microfossil analyses for proper interpretation of radiocarbon data, when studying mire fills. Additionally, they provide further proof for the development of podzols and ericaceous vegetation prior to and independent of early prehistoric agriculture in the Dutch Pleistocene coversand areas, and for the early (pre-medieval) start of wind erosion and deposition of drift sands.

The assessment of the seismic hazard and risk associated with the extraction of gas from the Groningen field involves a chain of modelling efforts. The first step is a description of the 3D distribution of reservoir properties in the reservoir – the static reservoir model – and is the subject of this paper. Consecutive steps in the chain of models are described elsewhere in this volume. The construction of a static reservoir model is not strictly a scientific endeavour, but many of the applied modelling techniques are underpinned by extensive scientific research. This paper aims to give a general introduction to the approach followed by NAM to build static models for the Groningen field. More detailed accounts of the applied modelling techniques, the assessment of associated uncertainties or the usage of multiple modelling scenarios are beyond the scope of the current paper, but are referenced in the text.

In the above mentioned article (Janssen et al., 2011) the caption of Figure 1 was not complete and missed a reference citation.

The vegetation history of the area around the confluence of the rivers Meuse and Swalm (the Netherlands) during the Middle Ages is covered by two pollen diagrams. The diagram Swalmen reveals a large-scale deforestation as a result of the foundation of a nobleman's homestead around 950. The diagram Syperhof shows a period during which the forest partly returns after a long history of unremitting anthropogenic stress. This temporary phenomenon is ascribed to the onslaught of the Black Death in 1349. Both diagrams provide evidence of the start of buckwheat growing.

This paper describes the landscape evolution of the Waasland Scheldt polders in the north of Belgium from the Late Glacial – early Holocene to the present time, and the effects of this changing landscape on the human settlement. The regional landscape evolution has been visualised in a series of palaeogeographical maps for successive time frames. Two different map series were produced: a series of Holocene palaeogeographical reconstructions (11,000–950 cal BP) based on geotechnical, geological and archaeological data, and a series of post-Medieval landscape reconstructions (16th- to 19th-century) based on historical maps, land registers and soil data. Additional palaeoenvironmental information from fossil pollen and plant remains allowed reconstruction of the vegetation and wetland changes, particularly for the middle to late Holocene. Peat growth was the main key to understanding the landscape evolution of the Waasland Scheldt polders. Whereas the landscape evolution during the Holocene was mainly sea-level driven, the transformation of the landscape during the last millennium was largely due to human interventions.

This paper presents the method applied to history-match the Groningen field dynamic model to conventional data (pressure data and water influx data) and to subsidence data, which is a novelty. Modelled subsidence is matched to subsidence data based on a simplified geomechanical model, which was built into the dynamic simulator.

A two-tier method was used to first achieve a match on a field-wide scale using field-average history-match quantifiers, which was then further improved at a regional/well level using regional history-match quantifiers. The history match was assisted by a space-filling experimental design. The resulting model has a field-average match to pressure data of ±2.17bar with a measurement uncertainty of ±0.4bar, to water influx data of ±2m with a measurement uncertainty of ±0.5m, and to subsidence data of ±4cm with a measurement uncertainty of ±1cm.

The output from this model is used as input for compaction, subsidence and production forecasts feeding into the hazard and risk assessment completed by NAM for the Groningen Winningsplan 2016.

An isolated, completely ossified vertebra tentatively ascribed to the non-cyamodontid placodont Pararcus diepenbroeki is described from the Anisian Vossenveld Formation in Winterswijk, the Netherlands, and compared to other material from the same locality. This fossil is the first completely ossified vertebra of the taxon and most likely originates from an adult specimen. It was recovered c. 16 m deeper in the stratigraphy than previously described material of the species, which is thus far known only from Winterswijk. Based on the slanting angle of the transverse process, the vertebra is interpreted to originate from the dorsal region. Besides the overall agreements in morphology that warrant a tentative identification as Pararcus diepenbroeki, the newly described vertebra deviates from other known Pararcus vertebrae in the presence of a longer, well-ossified neural spine and a strongly constricted, less pachyostotic and ovaloid vertebral centrum. General agreement in morphology with previously described vertebrae suggests this novel condition indicates a different anatomical position and perhaps a varied ossification pattern.

The infill of an abandoned branch of the small river Vlootbeek, a tributary of the Meuse, revealed a record of the vegetation history from the Allerød to the Middle Ages, albeit with at least two hiatuses in its upper reach. Striking observations are the continuous presence of Pinus in a period long after this tree was absent elsewhere in the Dutch landscape, and the high percentages of Tilia pollen in the mid-Holocene. The former is explained by the sandy–gravelly subsoil in the valley of the Vlootbeek. The latter is ascribed to the short distance between the coring location and the plateau on which Tilia must have been the dominant tree. Anthropogenic influence is apparent from the Early Neolithic Linearbandkeramik farmers onwards.

This paper describes a research programme recently initiated at Utrecht University that aims to contribute new, fundamental physical understanding and quantitative descriptions of rock and fault behaviour needed to advance understanding of reservoir compaction and fault behaviour in the context of induced seismicity and subsidence in the Groningen gas field. The NAM-funded programme involves experimental rock and fault mechanics work, microscale observational studies to determine the processes that control reservoir rock deformation and fault slip, modelling and experimental work aimed at establishing upscaling rules between laboratory and field scales, and geomechanical modelling of fault rupture and earthquake generation at the reservoir scale. Here, we focus on describing the programme and its intended contribution to understanding the response of the Groningen field to gas production. The key knowledge gaps that drive the programme are discussed and the approaches employed to address them are highlighted. Some of the first results emerging from the work in progress are also reported briefly and are providing important new insights.

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.

The Groningen field is the largest onshore gas field in Europe. The gas-bearing section comprises aeolian and fluvial Rotliegend sandstones of Permian age and fluvial sandstones of Carboniferous age. Continuous production since 1963 has led to induced seismicity starting in the early 1990s.

Faults at reservoir level play a major role in the seismicity in the Groningen field. Fault slip is expected when shear traction is sufficient to overcome frictional resistance on the fault surface. Clear insights into which faults and fault segments are most susceptible to seismicity could be used to optimise production and minimise the seismic risk. To gain these insights, a detailed and realistic fault model is required as input to both statistical analyses on seismicity and deterministic geomechanical modelling of seismogenic behaviour along fault planes. Geometrical seismic attributes and, subsequently, fault planes were extracted from a reprocessed and depth-imaged 3D seismic volume. This resulted in a detailed visualisation of the faults at reservoir level, with extension into the deeper strata below the reservoir in many cases. They represent fault planes with realistic dimensions and shapes. The fault map based on seismic attributes suggests the presence of faults that have not been included in studies on Groningen seismicity before. The improved fault definition correlates with recent earthquake hypocentres. We conclude that a detailed fault model of the Groningen field can be created using 3D seismic attributes and that detailed 3D fault planes can be extracted from these attributes. The results can be used as input to statistical and geomechanical analyses on 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.

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.