The Late Pleistocene surface paleohydrography of the New Jersey Pine Barrens consisted of a series of broad braided alluvial surfaces with meandering paleochannels. This drainage is best explained in terms of impermeable (i.e. frozen) substrate, high sediment load, variable or decreasing discharge, and eolian sedimentary dynamics. Evidence for eolian activity is provided by wind-abraded sand grains, coversand, dunes, ventifacts, deflation hollows and wind-polished boulders. In several places stream avulsion occurred due to channel infilling by locally-derived wind-blown sediment. The braided and meandering river systems that characterise the tundra and polar semi-desert lowlands of the Western Canadian Arctic are described as modern analogs.

In this study we apply an evidence-based approach to model population-size fluctuations and their corresponding impact on land use during the Roman and early-medieval periods in the Rhine–Meuse delta in the present-day Netherlands. Past-population numbers are reconstructed based on Roman and early-medieval settlement patterns. Corresponding impacts of these demographic fluctuations on potential land use are calculated by integrating the newly developed demographic overviews with archaeological and geoscientific data using a new land-use model termed ‘Past Land-Use Scanner’ (PLUS). The primary aims are to reconstruct first-millennium palaeodemographics and to explore the potential of simulation modelling for testing the feasibility of archaeological hypotheses regarding past land use. Results show that in the study area the first millennium AD was characterised by two periods during which major population growth occurred: the middle-Roman period (AD 70–270) and early-medieval period C (AD 725–950). A major demographic decline of 78–85% occurred during the late-Roman period (AD 270–450), after which first-millennium population numbers never again reached middle-Roman period levels. The modelling outcomes demonstrate that the impact of population fluctuations (growth vs decline) on the limits of the natural landscape during the first millennium in general was low. During these thousand years, the natural landscape almost without exception (only scenario D deviates) provided sufficient options for arable farming, meadows and pastures and was not a limiting factor for population growth. These results underline the added value of simulation modelling for testing the feasibility of archaeological hypotheses and analysing human–landscape interactions in the past.

Climate change is very likely to cause a global acceleration in sea-level rise (SLR). The projected acceleration of SLR will also affect the Wadden Sea. In addition to an accelerated SLR, gas and salt extraction will cause subsidence that adds to an increase in water depth in the tidal basins. This will have consequences for the sediment budget of the Wadden Sea and especially for the intertidal flats that have a high ecological value. This synthesis presents projections of the future state of the Dutch Wadden Sea for the years 2030, 2050 and 2100.

The projected changes in mean sea level by 2100 for Den Helder and Delfzijl are above the global mean projections, mainly due to the above-average ocean dynamics and glacio-isostatic adjustment contributions in the regional projections. The projected rise in mean sea level for 2100 with relation to 2018 in these locations is 0.41m, 0.52m and 0.76m for, respectively, the RCP2.6, RCP 4.5 and RCP8.5 climate scenarios.

When we combine the presented SLR scenarios with the subsidence estimates and compare these rates to the critical rates for ‘drowning’ of intertidal flats that were calculated for the individual tidal basins, we can determine the moment that the maximum imported sediment volume can no longer compensate the increase in accommodation space in a basin and the intertidal flats will start to diminish in surface area and/or height. In the RCP2.6 scenario, the projected rates of relative SLR will be below the critical rate for drowning of the inlet systems in the Dutch Wadden Sea. For the RCP4.5 scenario, the critical SLR rate will be exceeded for Vlie Inlet in 2030, and for the RCP8.5 scenario the critical SLR rate will be exceeded for Vlie Inlet in 2030, Texel Inlet in 2050 and Ameland Inlet in 2100. For the other basins the critical rate will not be exceeded until 2100 or later.

The way the intertidal flats in a basin will react to ‘drowning’ is not clear beforehand. It is highly possible that erosion of flats in one place will produce the sediment to maintain flats in other places. Tidal flats close to the sediment-delivering tidal inlet are not likely to disappear, because there the balance between supply and erosion is not likely to change.

A number of 45 presumed pingo remnants have been analysed on lithological criteria while measured parameters of pingo remnants can be compared on a spatial scale. Methods used were hand augering and radiocarbon dating. Out of 45 cored sites thirty sites proved to represent pingo remnants after having tested the results according to the set of descriptive criteria. Does a difference in filling types in a spatial analysis influence the choice of early settlers in this region? Results show that two groups of pingo remnants can be distinguished by the nature of the pingo filling. Type 1 pingo remnants with predominantly peat filling are found in the lower valleys between the higher glacial till ridges. Type 2 pingo remnants with equal loam and peat filling are predominantly found on the higher glacial till ridges. Type 2 pingo remnants generally have a greater depth, are more ellipsoid in shape, and have a steeper slope angle than type 1 pingo remnants on the lower terrain. Results indicate that the substrate character influenced the development of infilling of pingo remnants. Marshy environments may have existed for a prolonged period throughout the Holocene in pingo remnants on the glacial ridges, while pingo remnants on the lower terrain are filled in by peat bogs. Landscape evolution based on the results of a distinct dataset of pingo remnants can help to solve archaeological problems of settlement issues when combined with new palaeobotanic data and high resolution dates.

Bauxite, the raw material of aluminum, has been one of the economically vital natural resources for Suriname. Mining operations started about a century ago, and subsequent development of a refinery industry and hydro-electric power made Suriname one of the foremost bauxite and alumina producers worldwide for a long period of time. This paper presents a concise survey of the main geological attributes of its bauxite deposits and examines significant aspects in the development of mining in the country where alumina dominated the export revenues until a decade ago. The lateritic bauxite deposits are spread across the northern part of the country and developed on various parent rocks during Late Cretaceous–Early Tertiary times. Bauxites in the coastal lowlands formed on Cenozoic sedimentary deposits, whereas plateau bauxites originated on various crystalline rocks in inland regions of the Precambrian Guiana Shield. The composition of parent rocks and timing of bauxitisation point to a genetic correspondence with West African bauxites and a strong control of paleoclimatic conditions on the distribution and properties of bauxite in both regions. The more accessible bauxite deposits in the coastal lowlands are almost mined out, whereas the plateau bauxites have been extensively explored but have not been brought into production to date. For economic and environmental reasons, the future of the bauxite industry in Suriname is currently uncertain.

A brief overview is given of the history of geological mapping of theHolocene Rhine-Meuse delta. The first accurate map of the delta, based onfield observations, was made by Vink (1926). The geological map of theNetherlands, scale 1 : 50,000, made by the ‘Geologische Stichting’ (1927 –1938) under the supervision of P. Tesch totally neglected Vink’s work, andwas a step backwards with regard to the mapping of the Holocene delta.Between 1940 and 1965, the Wageningen group of soil scientists produceddetailed regional soil maps, that had a strong ‘geogenetic’ component. Inthe 1960’s a revolutionary ‘profile type legend’ was introduced by theNetherlands’ Geological Survey. This allowed to map not only the outcroppingsediments, but the whole Holocene succession, which gave more insight intothe geological history. Over the past 30 years, the Rhine-Meuse delta hasbeen studied extensively by students of physical geography at UtrechtUniversity. More than 250,000 borehole descriptions, 1500 14Cdates and over 36,000 archeological artifacts with associated ages(collected by the National Service for Archaeological Heritage) haveaccumulated, resulting in the largest database of a delta in the world. Theproduction of detailed maps has been crucial to the solution of manyscientific problems. The use of GIS has greatly enhanced geological andgeomorphological mapping, and subsequently, understanding of the evolutionof the Holocene Rhine-Meuse delta. A new detailed digital elevation map ofthe Netherlands, based on very accurate laser-altimetry data, will enable usto map larger areas in greater detail, with greater accuracy, and in a muchshorter period of time.

Accelerator mass spectrometry (AMS) radiocarbon (14C) dating of Cerastoderma edule (Linnaeus 1767) and Mytilus edulis (Linnaeus 1758) shells sampled in AD 1889 near the island of Wangerooge gave a new local correction factor ΔR of −85 ± 17 14C years for the Wadden Sea area. The value is considerably higher than the available scattered data from the North Sea, which were obtained from pre-bomb growth rings of living Arctica islandica (Linnaeus 1767). This can be explained by the incorporation of 14C-depleted terrestrial carbon into the shell material which compensates the intensified exchange of CO2 between atmosphere and shallow coastal water, e.g. by tidal currents. Additionally, two examples of application of the new ΔR value in coastal research give deeper insights into the dynamics of bivalve shell preservation in the Wadden Sea and the need for further research to clarify the Holocene reintroduction of Mya arenaria (Linnaeus 1758) into European waters.

We describe the occurrence in the Netherlands of three teeth of Mimomys hajnackensis (Arvicolidae, Rodentia) from three separate boreholes extending into matine deposits of the Maassluis Formation. These marine Pliocene/Early Pleistocene deposits at depth are overlain by, and interdigitate eastwards with, the freshwater deposits of the Waalre and Peize Formations (formerly known in part as the Tegelen Formation). Teeth identified as M. hajnackensis have not previously been described from the Netherlands although similar material has been reported from Pliocene deposits in Germany at Frechen and Hambach close to the Dutch border. The material described here comes from boreholes at the Polder de Biesbosch near Dordrecht, Diepenveen and Deventer. To our knowledge, these finds represent the oldest described rodent fossils from the Netherlands, corresponding to an absolute age of 2.9 to 3.3 Ma on the basis of comparison with other European localities, implying that the remains are reworked from older deposits.

This study presents the 3-D geometry of the Oligocene Rupel Clay Member and a review of its depositional environment based on new and published data. The Rupel Clay Member (RCM) is a clay layer in the Dutch subsurface which is informally known as the Boom Clay. New depth and thickness maps show that the member is present in nearly the whole subsurface of the onshore Netherlands to a depth of about 1500 m. The thickness of the member is variable but rarely exceeds 125 m. We identified three subunits: RCM-U, with a thickness of 15–40 m; RCM-M, with a thickness of 40–90 m; and RCM-L with a thickness of 25–50 m. The Rupel Clay Member is not a homogeneous clay; both vertical and lateral grain-size trends are present. These trends match with general palaeogeographic trends which foresee near-shore facies in the south and southeast of the onshore Netherlands and a more distal setting when moving northward. The three subunits are correlated with global climatic and regional tectonic events, which have mainly affected marginal facies. Faults are known to cut through the clay layer, but have not been mapped. A mismatch of lithostratigraphic nomenclature between the Netherlands and neighbouring countries is present. Part of the solution for this mismatch lies in the proper interpretation of well data in the Netherlands, both on- and offshore.

The geometry and depositional history of the Maassluis Formation is described from an East-West oriented transect located in the west-central Netherlands and P- and Q-blocks in the Dutch offshore area. The Late Pliocene to Early Pleistocene Maassluis Formation was deposited under near coastal marine conditions. Two distinct facies are distinguished: (1) medium grained sands with a blocky GR/SP-log pattern that were deposited under intertidal and aeolian conditions and (2) medium to coarse grained sands deposited in a subtidal environment showing a clear coarsening upward trend. The lower part of the Maassluis Formation is laterally equivalent to the shallow marine Oosterhout Formation and gets progressively younger towards the West. The upper part is lateral equivalent to estuarine and mudflat deposits (Balk Member, c.q. Peize Formation) in the central part of The Netherlands and to river deposits (Peize Formation) further to the East.

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.

Two closely spaced geothermal doublets were operated in the Californië geothermal field near Venlo, the Netherlands. The geothermal wells target the Dinantian Zeeland formation below 2 km depth. For several years, hot fluid was produced from the Tegelen fault, a regional fault in the Roer Valley rift system, until a felt M1.7 earthquake led to the suspension of geothermal activities. The Californië showcase provides a rare opportunity to retrospectively evaluate the assessment and the management of induced seismicity risks for a geothermal project. A seismic hazard assessment was conducted at several stages of the project, and seismicity was continuously monitored with a local station network.

In this paper, we report on the characteristics of the induced seismicity and evaluate the findings of the seismic hazard assessments conducted prior to the earthquakes. Seismic hazard assessments were based on numerical simulations of subsurface stress changes associated with geothermal operations. A geomechanical analysis indicated that the mapped faults in the subsurface are likely to be critically stressed. The largest hazard was inferred to result from thermo-elastic stresses, originating from cold water injection close to the Tegelen fault.

Subsequent earthquakes predominantly occurred near a production well after stopping or reducing production. We attributed this observation to a thermo-elastic stress load caused by cold water injection close to the Tegelen fault, combined with a counter-acting stabilisation of the fault due to pressure depletion during production. This mechanism was consistent with the dominating mechanism considered in the preceeding seismic hazard assessments. Although geothermal operations have not resumed yet, the geomechanical analysis indicates that re-locating one of the injection wells further away from the Tegelen fault could provide an efficient measure for mitigating induced seismicity risks at Californië.

Micromorphological studies relating the soil formation history (processes and timing) to activity events of the Peel Boundary Fault (PBF) showed rotation features (circular distribution pattern of sand grains) in mechanically displaced rounded fragments of Bt bands. These features are interpreted as being caused by ‘mudflow’ during active faulting event (PBF event F2). The micromorphological interpretation of Late Weichselian soil formation (clay illuviation, degradation features and offsetting of Bt bands) agrees with the hypothesized 3 PBF periods of fault activity events (Fl, F2 and F3).

Besides the Oligocene Boom Clay, the Ypresian clays – part of the Eocene Ieper Group (Kortrijk Formation and Kortemark Member) – are currently being investigated as an alternative host rock for the deep geological disposal of radioactive waste in Belgium and the Netherlands. In this study, broad-ion-beam milling and high resolution scanning electron microscopy (BIB-SEM) analyses were carried out to analyse the mineral fabrics and microstructures of representative Ypresian clay samples from different depths of the ON-Kallo-1 borehole (Kallo, Belgium). Qualitative microstructural observations indicate that mineral fabrics and pore morphologies in fine-grained samples are comparable to those found for fine-grained Boom Clay, but most of the Ypresian clay samples analysed also contained a significant silt fraction, which is associated with larger inter-aggregate pores, coated by a thin, very low porous clay layer. Quantitative pore-shape analysis shows lower axial ratios and elongations, as well as higher roundness and circularities for pores in the clay matrix of the more coarse-grained samples, compared to the fine-grained samples. The contribution of large pores (>1 × 107 nm² pore area) to the total BIB-SEM observed porosity was found to correlate with the non-clay mineral (NCM) content of a sample. Frequencies of pore sizes within the clay matrix follow a power-law distribution, hinting towards the possibility of up-scaling of the nanometre-scale observations to larger scale (micro-) structural features of the material. Power-law exponents are comparable to values found for power-law pore-size distributions within the clay matrix of the Boom Clay, which indicates similarity of the pore-space morphologies within the clay matrix of the Boom Clay and the Ypresian clays. Wood's metal injection, followed by (cooled BIB)-SEM analysis shows that all visible pores are connected via pore throats of diameter down to ~10 nm.

Ever since the first mosasaur restorations were published, these extinct marine reptiles have been pictured with either notched, forked or undivided tongues. Here, we present an overview of existing iconography, a review of the previous literature, and we discuss how best to reconstruct tongue form in mosasaurs. Despite disagreement about their precise phylogenetic position, most authors consider mosasaurs members of the Varanoidea, derived anguimorphans including Helodermatidae, Varanidae, Lanthanotus and probably snakes. All anguimorphans share a diploglossan (two-part) tongue, in which the foretongue is derived and modified into a highly protrusible chemosensor, while the hindtongue is plesiomorphic, retaining well-developed papillae, mucocytes and robust posterior lobes. We suggest that mosasaurs had a diploglossan tongue that remained in a relatively underived state. The form of the tongue would probably have been most like modern Heloderma or Lanthanotus with a protrusible chemosensory foretongue and a plesiomorphic, papillose hindtongue. Such a tongue is consistent with well-developed vomeronasal chemoreception through tongue-flicking, with the retention of the ancestral function of hyolingual food transport and swallowing following jaw-prehension of prey. The presence of paired fenestrae in the palate associated with the vomers, as well as the presence of pterygoid teeth are in accordance with such a tongue form in mosasaurs.

Presented are the first formal descriptions of mosasaur remains of Maastrichtian age from the Danish mainland (Jylland, Sjælland); just two taxa are recognised, viz. Mosasaurus cf. hoffmanni and Plioplatecarpus sp. Recent finds include an association of skeletal fragments, inclusive of a single marginal tooth crown, from within 10 metres or so of the K/T boundary near Holtug, Stevns Klint (Sjælland, eastern Denmark), here assigned to Plioplatecarpus sp. A brief review of previous records of Danish mosasaurs is added.

Two new records from the Amersfoort glacial basin are investigated by means of pollen analysis. The cores are situated in the deeper part, close to the original Eemian stratotype Amersfoort 1 (Zagwijn, 1961) and at the margin of the basin. The aim is to reconstruct the Eemian and Early Weichselian vegetation development and to explore the impact of accommodation space, influx of allochthonous pollen and geomorphology on the vegetation composition. The results of the Amersfoort Basin are compared to the current Eemian stratotype in the Amsterdam Basin and other Eemian sites in the Netherlands. An almost complete Eemian to Early Weichselian sequence (E2-EWII) was retrieved from the deeper part of the Amersfoort basin. The late Saalian (LS) to early Eemian transition is not recorded in the Amersfoort basin, in contrast to the deeper Amsterdam Basin. The basin marginal core Den Treek reveals a condensed late Eemian (E5-6) and Early Weichselian (EW I-II) succession showing the importance of accommodation space. The first impact of the Eemian transgression is registered at the E3 to E4a boundary in the Amersfoort and Amsterdam basins, and highest sea level is proposed at the end of pollen zone E5. Upstream in the Eemian delta, in the palaeo-Vecht valley and IJssel Basin, the transgression is recorded later.

The influx of reworked (allochthonous) pollen in clastic sediment units hampers vegetation and climatic reconstructions during the LS and Eemian. The early appearance of Picea in zone E4 and Abies in zone E5 in clastic sediment intervals can be related to long-distance transport by the river Rhine and redistribution in the Eemian delta.

Local vegetation development can complicate regional biostratigraphic correlations. Alnus, considered characteristic for the late Eemian (E5-6), shows large differences over short distances in the Amersfoort Basin, related to local alder growth since Eemian E3. Carpinus, diagnostic for pollen zone E5, shows high values in the basins adjacent to higher, well-drained ice-pushed ridges, but low values in low-relief environments. Salt- to brackish-water marshes were present during high sea level in zone E5 in the Amsterdam and Amersfoort basins, while further upstream in the Rhine delta brackish to fresh-water tidal conditions dominated.

In line with Zagwijn (1961), the E6 to EWI boundary is defined at the first opening of the vegetation cover with Calluna, Poaceae and Artemisia increase, often coinciding with a lithological change from organic to clastic deposition, reflecting increased landscape instability.

The cores from the Amersfoort basin reveal a complete Eemian to Early Weichselian record. It is suggested to define the boundary stratotype for the base of the Weichselian Stage in the Amersfoort Basin. The current stratotype Amsterdam-Terminal reveals a fully developed LS to Eemian transition and contains the boundary stratotype for the base of the Eemian Stage.

Alluvial deposits of the Loire/Arroux trunk/tributary system record distinct, synchronous episodes of regional fluvial adjustment. Changes in facies and depositional style through time can be interpreted with a modern analogue model that relates vegetative cover/human influence with sediment supply, and modes of atmospheric circulation with the paths and styles of storms that drive variable discharge regimes across western Europe.

Zonal atmospheric circulation results in a Mediterranean style climate over southern Burgundy, producing dry conditions punctuated by infrequent, large floods. Episodic overbank sedimentation and the burial of thin paleosols in sandy overbank facies is indicative of this style of fluvial activity, ca 1300 years BP. Humans may have increased the available volume of fine grained sediment at this times through increased agricultural activity along valley axes, however facies match that expected from a ‘flashy’ discharge regime.

In contrast, meridional circulation patterns result in a maritime style climate over southern Burgundy, with the intrusion of storms, moist conditions and frequent, moderate magnitude discharges. Wide, deep channels, thick channel facies and thin overbank facies are indicative of this style of fluvial activity, recorded in deposits dating to ca 4050 to 3200 years BP. Strong meridional conditions and extreme climatic variability during the Little Ice Age resulted in very large discharges that straightened and widened channels, while scouring and obscuring older terraces (ca 500 years BP). Deposition over the last two centuries is related to increasingly zonal circulation and infrequent, large (over-bank) floods. Changes in fluvial dynamics over the last 300 years can be attributed primarily to climatic control, as there has been very little change in land-use over that period.

The Suriname–Guyana and Ghana–Ivorian Basins exhibit strong geological similarities which are of interest from a petroleum exploration point of view. These include (1) a well-developed system of Late Cretaceous erosional canyons allowing coarse-grained shallow-water clastics to enter the deep marine basin to form attractive turbidite exploration targets; (2) a broad shelf with strong longshore currents which sort and transport coarse clastics into the canyon heads; (3) an organic-rich Cenomanian–Turonian hydrocarbon source rock which is thermally mature in the centre of the basin; and (4) a series of extensional fault-networks along the shelf margin that extend upwards from the Rift Sequence into the overlying Drift Sequence and which, along with the canyon geometries, enable migrating oil and gas to accumulate in combination structural and stratigraphic traps. In June 2007, Tullow and its partners made an important discovery in offshore Ghana, at Mahogany, which subsequently became the giant Jubilee field. Tullow is currently applying the same geoscientific technologies in offshore Suriname in the search for analogous subtle combination traps.

Part of the gravelly deposits of a combined Rhine-Meuse river of Middle Pleistocene age in the central Netherlands contains early Middle Palaeolithic artefacts. Although not in their original position, a significant part of these artefacts is hardly abraded, indicating limited fluvial transport. The artefacts have mainly been made from fluvial flint gravel boulders, originating from the Meuse catchment. Thus far, inferences for the age of the artefacts are based on the stratigraphic context and floral and faunal remains, which suggest a MIS 7 age. In this paper, OSL dating carried-out in the framework of a research aimed at the paleogeographical reconstruction of the Rhine-Meuse fluvial system in the central Netherlands and a review of published data are used to provide absolute age constraints for the artefact-bearing deposits. It is argued that the deposits were formed during the glacial phase directly preceding the Drenthe substage of the late Saalian (early MIS 6), and that at least a part of the artefacts has approximately the same age.