In this contribution, a case is made for the significance of sedimentation and sea-level change for groundwater management of coastal areas. In groundwater practice these geological processes are rarely considered. The role of sediment loading in causing anomalous fluid pressures and flow fields in relatively shallow aquifer systems is discussed and illustrated via both case studies and generic modelling studies. The role of sea-level changes in controlling current salinity distributions is discussed likewise. Central in the discussion is the concept of memory of groundwater systems, which provides the basic reason why processes that were operative in the geological past are still of relevance today. It is argued and shown that awareness and knowledge of the influence of sediment loading and sea level change on current hydrological conditions can lead to improved characterization of the distribution of hydraulic parameters and of the distribution of water quality in coastal areas. This improved characterization, in turn, serves to enhance the validity of impact assessment studies for the long-term development and management of those areas.

We present new local groundwater-level rise data from two Late Glacial aeolian dunes, located near Barendrecht and Oud-Alblas in the western Rhine-Meuse delta. These data are based on AMS radiocarbon dating of terrestrial macrofossils, collected from the base of peat formed on the slopes of these dunes. This method avoids contamination of bulk peat samples by old soil carbon or younger rootlets and rhizomes, as well as the hardwater effect. The new data are used to assess the reliability of previously published groundwater-level index data based on conventional radiocarbon dating of bulk basal peat samples from the slopes of the Late Glacial aeolian dunes at Barendrecht, Hillegersberg, Bolnes and Wijngaarden, all located in the western Rhine-Meuse delta.

Comparison of the new and published groundwater-level data shows no significant systematic difference between conventionally dated bulk peat samples and AMS-dated samples of terrestrial macrofossils. The new data from the dune at Barendrecht confirm the reliability of the younger than 6600 cal yr BP age-depth data from the dunes at Hillegersberg and near Bolnes. This result supports the validity of this part of the mean sea-level (MSL) curve for the western Netherlands. Consequently, the position of the groundwater-level curve for Flevoland (central Netherlands) below this MSL curve can most likely be attributed to differential land-level movement.

The available data show that the groundwater-gradient effect in the western Rhine-Meuse delta became less than 5 cm/km after 6600 cal yr BP. Finally, temporal correlation between temporary increases in local groundwater-level rise with known shifts of river courses in the delta plain suggests, that avulsions can explain sudden local deviations from the trend in groundwater-level rise. A general conclusion of this study is that a complex relationship exists between sea level and local delta-plain water levels.

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.

We have built a 3D lithological model of the Netherlands, for the purpose of mapping on-land aggregate resources down to 50 m below the surface. The model consists of voxel cells (1000 ⋅ 1000 ⋅ 1 m), with lithological composition and aggregate content estimates as primary attributes. These attributes were derived from ∼350,000 borehole descriptions. Overburdens and intercalations of cohesive or otherwise non-dredgeable materials were taken into account to define geologically exploitable aggregates within the total stock. We arrive at about 520 ⋅ 109 m3 of aggregates occurring in the depth range investigated. Some 50% of this amount is considered geologically exploitable and about 25% would in principle (but largely not in reality) be accessible. Most aggregates resources (∼98%) are coarse sand, which is processed for use in concrete, masonry mortars, drains, filters, etc. The total exploitable stock of coarse sand in the depth range investigated amounts to roughly 7500 times the current annual consumption level, and is virtually indepletable. The gravel stock, estimated at some 12 ⋅ 109 m3, is small by comparison, and impels a dependency on imports.

Exploitable aggregates occur in all but the coastal provinces. In accordance with current policy changes, the future may show a shift from concentrated production along the upstream Dutch Rhine and Meuse rivers towards a more even distribution of small-sized operations over the country. Fairly large aggregate stocks, that have not yet been exploited to significant extent, are available in the northern extent of the aggregates occurrences.

Two genera of tylosaurine mosasaurs, Tylosaurus and Hainosaurus, are recorded for the first time from Germany. Tylosaurus sp. is represented by two isolated tooth crowns, originally described as Mosasaurus? alseni (here considered a nomen dubium) from the latest Santonian–Early Campanian, which are very similar to T. ivoensis and T. gaudryi. The material of Hainosaurus sp. comprises a maxillary with associated postorbitofrontal, two pterygoid teeth and several indeterminate cranial fragments. The specimen from the Late Campanian is slightly less derived than H. bernardi from the Maastrichtian in retaining labiolingually less compressed anterior maxillary teeth and unserrated pterygoid teeth with only very weak carinae. Despite only minor skeletal differences, the genus Hainosaurus is considered to be distinct from Tylosaurus because of its significant modification of the dental apparatus compared to the plesiomorphic condition in the latter. This dental morphology suggests a phylogenetic trend from a generalised-piercing marginal dentition in Tylosaurus towards the increasingly labiolingually compressed, symmetrical, strongly bicarinate cutting marginal teeth in Hainosaurus spp. from the Early through Late Campanian and Maastrichtian. A similar trend is also present in pterygoid teeth with very indistinct unserrated carinae in the Campanian Hainosaurus sp. towards serrated ones in the Maastrichtian H. bernardi. A short review indicates the presence of Hainosaurus in northern, central and western Europe (Sweden to Spain) since the Early Campanian, and the occurrence of Tylosaurus spp. in the same area until the Late Campanian. Hainosaurus persisted until the end of the Maastrichtian; outside Europe it may have been present in the Late Campanian of the USA and the Maastrichtian of the Democratic Republic of Congo. Judging from a simple, uni- to bicarinate, stoutly conical tooth morphology in aigialosaurs and very basal mosasaurs as well as phylogenetic patterns, the development of blade-like cutting tooth crowns appears to have been convergent in several clades of large-bodied Campanian–Maastrichtian mosasaurids. These include both mosasaurines ('Leiodon' mosasauroides, Prognathodon? sectorius, Prognathodon? kianda, Eremiasaurus heterodontus) and tylosaurines (Hainosaurus spp.).

Sequence stratigraphic correlations indicate that intermittent changes of the kinematic far-field stress-field regimes, and the associated geodynamic re-organisations at the plate-tectonic contacts of the African, Apulian, Iberian and European plates, affected the Tertiary palaeogeographic evolution of the West European Platform through a combination of intra-plate tectonics and fluctuations of relative sea level. A temporal sequence of first-order stages in structural, palaeotopographic and palaeohydrographic development of the platform can be distinguished from the Paleocene onwards. These formative stages are closely linked to major plate-boundary events involving the development of the Pyrenean and Alpine orogens, and can be traced throughout the West European Platform.

Five new dentaries, representing three individuals of the poorly known durophagous mosasaur Carinodens are described from the Maastrichtian phosphates of Morocco. Contrary to previous assumptions, the dentary of Carinodens holds 17-18 teeth, rather than 13. One pair of dentaries represents a new species, Carinodens minalmamar, characterised by laterally flattened teeth with two relatively pronounced sulci, a tooth count of 18, and a more slender general aspect of the dentary itself. The new material requires adjustment of previous reconstructions, assumptions on jaw mechanics, and palaeobiological interpretations of Carinodens.

In the Weichselian, the Lower Rhine in the Dutch-German border region has used three courses, dissecting ice-marginal topography inherited from the Saalian. In the Late Weichselian, the three courses functioned simultaneously, with the central one gaining importance and the outer ones abandoning. This study aims to reconstruct the fluvial development and forcings that culminated in abandonment of the northern branch ‘Oude IJssel-Rhine’, at the time of the Lateglacial to Holocene transition. The fluvial architecture is studied using a cored transect over the full width of the valley, detailed cross-sections over palaeochannels and geomorphological analysis using digital elevation and borehole data. Biostratigraphy, radiocarbon dating and OSL dating provide a timeframe to reconstruct the temporal fluvial development. In its phase of abandonment, the fluvial evolution of the Oude IJssel-Rhine course is controlled by the ameliorating climate and related vegetation and discharge changes, besides by intrinsic (autogenic) fluvial behaviour such as the competition for discharge with the winning central branch and the vicinity of the Lippe tributary confluence. The rapid climate warming at the start of the Late Glacial resulted in flow contraction as the initial response. Other fluvial geomorphic adjustments followed, with some delay. An aggrading braided or transitional system persisted until the start of the Allerød, when channel patterns finally changed to meandering. Floodplain incision occurred at the Allerød - Younger Dryas transition and a multi-channel system developed fed by Rhine discharge. At the start of the Holocene, this system transformed into a small-scale, local meandering system, which was abandoned shortly after the start of the Holocene.

The final abandonment of the Oude IJssel-Rhine and Niers-Rhine courses can be attributed to deep incision of the Central Rhine course in the earliest Holocene and is considered to be controlled by flow contraction induced by climate and related vegetation and discharge changes.

The records of the last interglacial in several palaeogeographical zones and various sedimentary environments from Central Asia to the Black Sea shoreline are presented. The last interglacial is characterized in the study areas by a two- or three-unit stratigraphical complex with both terrestrial and marine formations. Finds of significant small mammalian remains, together with a characteristic mollusc fauna in the Karangatian marine sediments, provide key levels for biostratigraphical correlation between the marine and continental deposits. New U/Th dates allow the correlation of the Karangatian transgression with δ18O substage 5e.

The impact of oil and, in particular, gas fields discovered in the Dutch subsurface has been very significant. However, 50 years after the discovery of the giant Groningen gas field the Netherlands has become very mature for exploration of oil and gas, and the gas volume left to be discovered in conventional traps is insignificant compared to what has been found already. The total portfolio of conventional prospects held by the industry contains several 100s of billions of cubic metres (bcm), as reported by the Ministry of Economic Affairs, but many of these prospects are unattractive to drill because of their small size or other geologically unfavourable aspects. Hence, for planning purposes of future national gas production the risk should be taken into account that the size of the conventional portfolio is overestimated. The major E&P companies have reduced their exploration efforts and the number of wells drilled as well as the size and total volume of discovered gas reserves has seen a steady decline over the last 10 years. Some surprises may still be in store and can occasionally add a welcome addition of gas. But the follow-up potential of new play and trapping concepts has been disappointing for many years now, and it is concluded that this is unlikely to be different in the future. Remaining conventional discoveries will mainly be in small near-field targets that as a result of technological advances made in the last few decades can be drilled with high confidence, despite their small volumes.

This leaves the so-called unconventional gas (UG) resources for a real and significant increase in the exploration potential of the Netherlands. UG resources occur outside conventional structural or stratigraphic traps in tight (low permeability) rocks and are of regional or sub-regional extent, without well-defined hydrocarbon-water contacts. The potential for Basin Centred Gas, Shale Gas and Coal Bed Methane is reviewed. As, according to present-day technology, development of UG requires very dense drilling at low costs with well spacing of a few 100s of metres, only the onshore potential can be commercial, even in the longer term.

Recent geological uplift is a characteristic for all North American commercial UG developments. Uplift helps bringing the resources close to the surface and facilitates development of fractures, which are essential for achieving commercial flow rates. This significantly reduces the area where commercial UG resources may occur in the Netherlands. In addition, sweet spots, where commercial flow rates and ultimate recovery per well can be achieved, represent only a fraction of the total ‘play area’. The UG plays in the Dutch subsurface remain to be proven, and there is still a significant technical risk associated with these plays, on top of the commercial risk. Therefore, despite potentially enormous in-place gas volumes in these unconventional plays, recoverable volumes are much less. If UG resources can be proven and are commercially developable, their cumulative volume potential is estimated by us in the order of a few tens to one or two hundreds bcm of recoverable gas at best. Finally, as UG resources produce at very low rates and require large numbers of wells to develop, the environmental impact in a densely populated country like the Netherlands is enormous, and needs to be seriously considered, already in the exploration phase.

In a mature area like the Netherlands, industry focus should be on technology development to reduce risk, increase recovery, reduce cost and minimize surface impact. Cooperation between Operators to build multi-well campaigns is therefore strongly recommended to reduce mobilisation cost. In addition, government incentives should be targeted at the development phase, in order to increase economic attractiveness for difficult reservoirs, both conventional and unconventional. In this way State and industry will both be able to maximize their returns on the remaining potential for gas and oil in the next two to three decades.

Hydrogeological research in coastal areas has gained considerable attention over the last decades due to increasing stresses on fresh groundwater resources. Fundamental groundwater flow and solute transport analyses remain essential for a concise understanding of the governing processes that lead to salinisation of fresh groundwater resources. However, the challenge of modern research is the application and quantification of these processes in real world cases. In this context, deltaic areas are amongst the most difficult study areas as they often have a complex groundwater salinity distribution. The Wadden Sea area in the northern part of the Netherlands is an example of such an area.

We quantified salt water intrusion and salinisation of groundwater flow systems in two representative case studies in the Wadden Sea area, using the density dependent groundwater flow and transport code M0CDENS3D. The results indicate that sea-level rise and autonomous processes will cause severe salinisation in the future, especially in the low polder areas close to the sea. In addition, we show that enhanced land subsidence due to salt exploitation accelerates this process. Salinisation can be mitigated to some extent by raising surface water levels in polders and by creating saline groundwater collection areas that maintain a low controlled water level.

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.

Rising sea levels due to climate change can have severe consequences for coastal populations and ecosystems all around the world. Understanding and projecting sea-level rise is especially important for low-lying countries such as the Netherlands. It is of specific interest for vulnerable ecological and morphodynamic regions, such as the Wadden Sea UNESCO World Heritage region.

Here we provide an overview of sea-level projections for the 21st century for the Wadden Sea region and a condensed review of the scientific data, understanding and uncertainties underpinning the projections. The sea-level projections are formulated in the framework of the geological history of the Wadden Sea region and are based on the regional sea-level projections published in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR5). These IPCC AR5 projections are compared against updates derived from more recent literature and evaluated for the Wadden Sea region. The projections are further put into perspective by including interannual variability based on long-term tide-gauge records from observing stations at Den Helder and Delfzijl.

We consider three climate scenarios, following the Representative Concentration Pathways (RCPs), as defined in IPCC AR5: the RCP2.6 scenario assumes that greenhouse gas (GHG) emissions decline after 2020; the RCP4.5 scenario assumes that GHG emissions peak at 2040 and decline thereafter; and the RCP8.5 scenario represents a continued rise of GHG emissions throughout the 21st century. For RCP8.5, we also evaluate several scenarios from recent literature where the mass loss in Antarctica accelerates at rates exceeding those presented in IPCC AR5.

For the Dutch Wadden Sea, the IPCC AR5-based projected sea-level rise is 0.07±0.06m for the RCP4.5 scenario for the period 2018–30 (uncertainties representing 5–95%), with the RCP2.6 and RCP8.5 scenarios projecting 0.01m less and more, respectively. The projected rates of sea-level change in 2030 range between 2.6mma−1 for the 5th percentile of the RCP2.6 scenario to 9.1mma−1 for the 95th percentile of the RCP8.5 scenario. For the period 2018–50, the differences between the scenarios increase, with projected changes of 0.16±0.12m for RCP2.6, 0.19±0.11m for RCP4.5 and 0.23±0.12m for RCP8.5. The accompanying rates of change range between 2.3 and 12.4mma−1 in 2050. The differences between the scenarios amplify for the 2018–2100 period, with projected total changes of 0.41±0.25m for RCP2.6, 0.52±0.27m for RCP4.5 and 0.76±0.36m for RCP8.5. The projections for the RCP8.5 scenario are larger than the high-end projections presented in the 2008 Delta Commission Report (0.74m for 1990–2100) when the differences in time period are considered. The sea-level change rates range from 2.2 to 18.3mma−1 for the year 2100.

We also assess the effect of accelerated ice mass loss on the sea-level projections under the RCP8.5 scenario, as recent literature suggests that there may be a larger contribution from Antarctica than presented in IPCC AR5 (potentially exceeding 1m in 2100). Changes in episodic extreme events, such as storm surges, and periodic (tidal) contributions on (sub-)daily timescales, have not been included in these sea-level projections. However, the potential impacts of these processes on sea-level change rates have been assessed in the report.

Over the last decades luminescence dating techniques have been developedthat allow earth scientists to determine the time of deposition ofsediments. In this contribution we review: 1) the development of themethodology; 2) tests of the reliability of luminescence dating onNetherlands’ sediments; and 3) geological applications of the method in theNetherlands. Our review shows that optically stimulated luminescence datingof quartz grains using the single aliquot regenerative dose method yieldsresults in agreement with independent age control for deposits ranging inage from a few years up to 125 ka. Optical dating of quartz has successfullybeen applied to sediments from a wide range of depositional environmentssuch as coastal dunes, cover sands, fluvial channel deposits, colluvialdeposits and fimic soils. These results demonstrate that optical dating is apowerful tool to explore the natural archive of the Netherlands’subsurface.

The presence of permafrost in Poland north of the line indicating the maximum ice extent of the Vistulian (Weichselian) glaciation after retreat of the land-ice cap has been a subject of debate for a long time. Investigations in an area at the line of the maximum ice extent of the Pomeranian phase prove that permafrost existed, indeed, after the ice retreat. This conclusion is drawn on the basis of morphological data (the presence of oriented kettle holes), sedimentological data (the nature of the infilling of the kettle holes) and pedological data (permafrost-affected horizons in soil profiles). It appears that the permafrost mostly developed in the ice-free zone that appeared after the retreat of the land-ice cap, but it is likely that some relict permafrost that had originated earlier in front of advancing ice was also still present. The landscape of northern Poland owes its relief partly to the Late Glacial permafrost.

A new site with Lateglacial palaeosols covered by 0.8 - 2.4 m thick aeolian sands is presented. The buried soils were subjected to multidisciplinary analyses (pedology, micromorphology, geochronology, dendrology, palynology, macrofossils). The buried soil cover comprises a catena from relatively dry (’Nano’-Podzol, Arenosol) via moist (Histic Gleysol, Gleysol) to wet conditions (Histosol). Dry soils are similar to the so-called Usselo soil, as described from sites in NW Europe and central Poland. The buried soil surface covers ca. 3.4 km2. Pollen analyses date this surface into the late Allerød. Due to a possible contamination by younger carbon, radiocarbon dates are too young. OSL dates indicate that the covering by aeolian sands most probably occurred during the Younger Dryas. Botanical analyses enables the reconstruction of a vegetation pattern typical for the late Allerød. Large wooden remains of pine and birch were recorded.

An inventory of the available surficial data on active faults in Italy has been compiled by gathering all the available information on peninsular Italy (project by CNR, National Group for the Defense against Earthquakes), the central-eastern Alps and the Po Plain (EC ‘PALEOSIS’ project). Such information has been summarised in maps (reporting surficial expressions of faults with length L≥11 km) and in a table where fault parameters relevant for seismic hazard assessment (e.g. slip rates, recurrence intervals for surface faulting events, etc..) have been reported. Based on the geological characteristics of the Italian territory, a fault has been considered as active if it shows evidence of Late Pleistocene-Holocene displacements. Active faults in Italy are distributed throughout the entire Apennine chain, in the Sicilian and Calabrian regions and in some Alpine sectors, but knowledge is not homogeneously distributed through the territory. The largest amount of data is related to the central Apennines. In contrast, fault geometries and parameters are less well defined in the southern Apennines, Sicily and Calabria, where investigations have started more recently. Knowledge is sparse in the northern Apeninnes, where data necessary to define fault parameters are lacking and also the chronology of the activity has to be considered cautiously. Abundant blind faulting in the Po Plain hinders the detection of active faults by means of the classical surficial investigations and therefore the present knowledge is limited to the Mantova fault. Blind faults and the peculiar recent geological history of the Alpine areas, which is strongly conditioned by the erosional and depositional activity during and after the last glacial maximum, also hinder the identification of active faults in the central-eastern Alps. Some faults in this Alpine sector are believed to be active, but data on their segmentation are still missing. Available information indicates that Italian active faults are usually characterised by slip rates lower than 1 mm/yr. Recurrence intervals for surface faulting events are longer than 1,000 years in the central and southern Apennines. This review on the Italian active faults represents the first step to produce a map of the major seismic sources in Italy, which in turn will result from the merge of surficial data with seismological and geological subsurficial data. The available knowledge gathered in this paper indicates those areas where data are presently sparse. It should be, therefore, possible to better plan future geomorphological and paleoseismological investigations.

During the last decades a considerable amount of data on mammals from Mikulino (=Eemian) deposits of the central and southern parts of the Russian Plain has been produced. Mammuthus primigenius (the early type), Palaeoloxodon antiquus (the advanced form), Arvicola ex gr. terrestris, Eolagurus cf. luteus and Lagurus cf. lagurus characterize this period. The so-called ‘Shkurlatian mammal assemblage’ was distinguished on the basis of a number of Eemian faunas.

The age of the mammal localities was established by using various geological and palaeontological evidence, together with the results of palaeomagnetic studies and absolute dating. The evolutionary level of the diagnostic species allows correlations between localities from different parts of the Russian Plain that have varying taphonomical conditions. The mammal- and malacofauna of the Mikulino Interglacial can also be used as a basis for correlations between the Karangat marine strata of the Black Sea and the continental deposits (alluvial sediments of the second terrace of the Sudost’ River and the Salyn phase of the Mezin fossil pedocomplex (see Dodonov et al., this volume). Forest, forest-steppe and steppe zone landscapes have been reconstructed for the central and southern parts of the Russian Plain on the basis of the Eemian faunas.

A unique complex of multiphased Holocene drift sands and paleosols, with at least two lacustrine phases, was discovered during a major sanitation project in the Laarder Wasmeren area near Hilversum, the Netherlands. The complex could be studied in detail, highly facilitated by the excellent and large-scale exposure of the various deposits and soils. OSL dating was used to establish ages of the phases, where possible differentiating between time of deposition and time of burial by taking into account the potential effects of bioturbation. Pollen analysis served to reconstruct the vegetation during the various phases.

A first minor phase of aeolian activity already started before approx. 5,000 BC, followed around 4,000 BC by a second phase and a rather massive third phase around 3,000 BC. After a long phase of soil formation, the latest, massive drift sand phase started around the 14th to 15th century. It clearly represents the classic drift sand phase that started in the Late Middle Ages in the Netherlands. Sand drifting followed on soil forming phases during which the vegetation became increasingly dominated by ericaceous plants and culminated in heathlands. The first three aeolian phases and associated heathlands are much older than generally assumed for heathland and drift sand to occur in the Netherlands. Moreover, podzolisation was found to have started very early, true podzols already occurring before 4,000 BC.

Around 3,000 BC groundwater in the area reached a maximum altitude of about 230 cm +NAP, resulting in local open water in the area. This rise is probably linked to the development of the Dutch coastal area, where at that time peat accumulated and drainage was poor, inducing a rise of the groundwater level in ‘het Gooi’. This groundwater level fell later on, to never reach this altitude again. The Groot Wasmeer was formed by local stagnation on a slowly permeable podzol and already reached a level of 320-325 cm +NAP by 400 BC, which more or less equals its 20th century level.

The results demonstrate that earlier concepts on the occurrence and age of aeolian phases, podzols and heathland vegetations in the Netherlands are far too schematic, and that early, pre-agricultural cultures may already have had an impact on the stability of fragile cover sand landscapes, e.g. through burning. Results are in line with those from several contemporary studies on early prehistoric cultures and their impact in river dune areas in the Central and Eastern Netherlands.

The carbon isotope composition of tooth enamel from all five mosasaur taxa known to date from the type Maastrichtian (southeast Netherlands, northeast Belgium) has been analysed. Differences in enamel δ13C values between taxa suggest resource partitioning. Body size and δ13C value correlate surprisingly well, with larger taxa having δ13C values more depleted in the heavier 13C isotope.