The fluvial history of the northern Lower Rhine Embayment shows interplay of three main river systems: Rhine, Meuse and smaller rivers draining the central and northern part of Belgium.

The Pliocene and Early Pleistocene (pre-)Rhine and Meuse river systems had their conjunction in the southern part of the Roer Valley Graben between Aachen and Jülich. Despite slight differences in the heavy-mineral assemblages the lithological composition of the Pliocene deposits of the three river systems shows close resemblance and therefore they cannot be mapped separately. However, due to a marked change of the petrographical composition the Upper Pliocene and Lower Pleistocene deposits of the Rhine are easily recognised and as a result Rhine and Meuse deposits can be mapped separately upstream of their confluence.

The Lower Pleistocene deposits of Rhine, Meuse and the Belgian rivers show a clear interrelationship. They are bounded by two regional well-mapable unconformities and are preserved in from west to east changing lithostratigraphical sequences. Revision of the lithostratigraphical schemes in Germany and the Netherlands and the better defined lithostratigraphical position of Meuse deposits in Germany now strongly constrain the correlation of the various fluvial deposits. As a result existing reconstructions of the fluvial deposition and tectonic history of the southern Roer Valley Graben can be evaluated and re-adjusted.

It is concluded that the main course of the Meuse was aligned through the so-called East Meuse valley during the larger part of the Early Pleistocene. Available pollen data do not conflict with this conclusion. At the same time the Rhine ceased to enter the southern part of the Roer Valley Graben. Instead, the Meuse accumulated here a series of deposits derived from the East-Meuse valley. Simultaneously, the Belgian rivers filled available accommodation space in the Roer Valley Graben of the southern Netherlands. The conclusions are based primarily on the revised lithostratigraphical framework. In general they simplify the picture of fluvial and tectonic behaviour of the area.

Geothermal energy is a viable alternative to gas for the heating of buildings, industrial areas and greenhouses, and can thus play an important role in making the transition to sustainable energy in the Netherlands. Heat is currently produced from the Dutch subsurface through circulation of water between two wells in deep (1.5–3 km) geothermal formations with temperature of up to ∼100 °C. As the number of these so-called doublets is expected to increase significantly over the next decades, and targeted depths and temperatures increase, it is important to assess potential show-stoppers related to geothermal operations. One of these potential hazards is the possibility of the occurrence of felt seismic events, which could potentially damage infrastructure and housing, and affect public support. Such events have been observed in several geothermal systems in other countries. Here we review the occurrence (or the lack) of felt seismic events in geothermal systems worldwide and identify key factors influencing the occurrence and magnitude of these events. Based on this review, we project the findings for seismicity in geothermal systems to typical geothermal formations and future geothermal developments in the Netherlands. The case study review shows that doublets that circulate fluids through relatively shallow, porous, sedimentary aquifers far from the crystalline basement are unlikely to generate felt seismic events. On the other hand, stimulations or circulations in or near competent, fractured, basement rocks and production and reinjection operations in high-temperature geothermal fields are more prone to induce felt events, occasionally with magnitudes of M > 5.0. Many of these operations are situated in tectonically active areas, and stress and temperature changes may be large. The presence of large, optimally oriented and critically stressed faults increases the potential for induced seismicity. The insights from the case study review suggest that the potential for the occurrence of M > 2.0 seismicity for geothermal operations in several of the sandstone target formations in the Netherlands is low, especially if faults can be avoided. The potential for induced seismicity may be moderate for operations in faulted carbonate rocks. Induced seismicity always remains a complex and site-specific process with large unknowns, and can never be excluded entirely. However, assessing the potential for inducing felt seismic events can be improved by considering the relevant (site-specific) geological and operational key factors discussed in this article.

A 3D geological raster model has been constructed of the onshore of the Netherlands. The model displays geological units for the upper 500 m in 3D in an internally consistent way. The units are based on the lithostratigraphical classification of the Netherlands. This classification is used to interpret a selection of boreholes from the national subsurface database. Additional geological information regarding faults, the areal extent of each unit and conceptual genetic models have been combined in an automated workflow to interpolate the basal surfaces of each unit on 100 × 100 metre (x,y dimensions) raster cells. The combination of all interpolated basal surfaces results in a 3D Digital Geological Model (DGM) of the subsurface. A measure of uncertainty of each of these surfaces is also given. The automated workflow ensures an easily updatable subsurface model. The outputs are available for end users through www.dinoloket.nl.

The diversity of morphological features of Mesozoic inversion structures in NW Germany as representatives of inversion structures in northern Europe is presented and their origin analysed and geologically dated. The particular role of salt in inverted basins and the re-shaping of pre-existing salt structures during the inversion act is demonstrated and the term ‘salt wedge’, a Zechstein salt intrusion into salt layers within the Triassic sedimentary pile, introduced. The leading theories on inversion (continent-continent collision, re-activation Variscan features) are discussed and discarded, but no new comprehensive theory was developed. The impact of inversion on HC prospectivity of sedimentary basins is debated and proposals for future interdisciplinary research are made.

During the Pleistocene the drainage pattern in the Lower Rhine Basin changed twice, from a flooding of the whole basin by the river Rhine from SW to NE to an influence restricted to the NE only.

The first dominance of the river Rhine is documented from the Reuverian to the Tiglian, the second one in the Cromerian. In between this time, the Meuse River drained the central Lower Rhine Basin in NE direction. For the sediments of that river, the term ‘Holzweiler Formation’ is introduced. Since the Late Cromerian, the influence of the Rhine is again restricted to the NE of the Lower Rhine Basin. The central part of the basin is drained by small local rivers.

In this paper we present a probabilistic fast model for performance assessment of geothermal doublets for direct heat applications. It is a simple yet versatile and multipurpose tool. It can be well applied in better understanding the sensitivity of performance to key subsurface parameters and depth trends therein, and for assessing the probability of success for geothermal projects under technical and financial constraints.

The underlying algorithms deliver a sensible accuracy given the uncertainties associated with geothermal projects at exploration state. A public release of the software, available under the name of DoubletCalc, is easy to handle and requires a limited set of input parameters. Thanks to an open source code, DoubletCalc can be implemented in other software applications and extended as it has been implemented for the integration into the national geothermal information system in the Netherlands (ThermoGIS, 2011).

Apart from its application for site assessments, the tool can be integrated into automated workflows processing faster representations of key aquifer properties and capable to produce indicative maps for predicted doublet power, economic feasibility and prediction of cumulative amount of heat that can be recovered. These capabilities are specifically important for decision support for policymakers while assessing the effects of particular insurance schemes and funding mechanisms.

DoubletCalc cannot and is not intended to substitute geologic exploration approaches. As exploration measures, such as seismic surveys are cost intensive, DoubletCalc can be used to focus geothermal exploration on areas and sites where an enhanced probability of success can be expected.

The Southern Alps – Ligurian basin junction is one of the most active seismic areas in Western Europe countries. The topographic and the structural setting of this region is complex because of (i) its position between the high topography of the Southern Alps and the deep, narrow Ligurian oceanic basin, and (ii) the large number of structures inherited from the Alpine orogeny. Historical seismicity reveals about twenty moderate-size earthquakes (up to M=6.0), mostly distributed along the Ligurian coast and the Vésubie valley. A recent geodetic experiment shows a significant strain rate during the last 50 years in the area between the Argentera massif and the Mediterranean coastline. Results of this experiment suggest a N-S shortening of about 2-4 mm/yr over the network, this shortening direction is consistent with the seismological (P-axes of earthquakes) and the microtectonic data. The Pennic front (E-NE of the Argentera massif) and the northern Ligurian margin are the most seismically active areas. In the Nice arc and in the Argentera massif, some seismic lineaments correspond to faults identified in the field (such as theTaggia-Saorge fault or the Monaco-Sospel fault). In the western part of the Alpes Maritimes, no seismic activity is recorded in the Castellane arc. In the field, geological evidence, such as offsets of recent alluvial sediments, recent fault breccia, speleothem deformations, radon anomalies and others indicates recent deformation along these faults. Nevertheless, to this date active fault scarps have not been identified: this probably results from a relatively high erosion rate versus deformation rate and from the lack of Quaternary markers. We also suspect the presence of two hidden active faults, one in the lower Var valley (Nice city area) and the other one at the base of the Argentera crustal thrust-sheet. Offshore, along the northern Ligurian margin, the seismic reflection data shows traces of Quaternary extensional deformation, but the accuracy of the data does not yet allow the construction of a structural map nor does it allow the determination of the continuity between the offshore and onshore structures. From these data set we propose a preliminary map of 11 active faults and we discuss the questions which remain unsolved in the perspective of seismic hazard evaluations.

From the 40s A.D. onwards a dense military system was established in the Lower Rhine delta in the Netherlands. Long since, it is questioned why this system was established in a wetland area and even turned into the northwest frontier of the Roman Empire, the Limes. A new detailed palaeogeographical map, based on a digital elevation model (LIDAR), soil maps and excavation results, was constructed. This reconstruction provides insight and understanding of the interactions between the natural environment in this part of the delta on the one hand and the establishment of this part of the Limes along the Old Rhine between Utrecht and Katwijk on the other. This study shows that the distinctive landscape of the western Rhine-Meuse delta, with an exceptionally large number of tributaries, determined the spatial pattern of the military structures. All forts (castella) were erected on the southern natural levees of the river Rhine, directly alongside the river, regardless of height and composition of the subsoil and alongside or opposite routes that provided natural access to the river. We conclude that their aim was to guard all waterways that gave access to the river Rhine from the Germanic residential areas further north and from/to the Meuse tributary further south in the delta. In addition, a system of small military structures, mostly watchtowers, was erected between the forts to watch over the river Rhine and its river traffic. Furthermore, at least two canals were established to create shorter and safely navigable transport routes to the river Meuse. At first, this integrated system of castella and watchtowers probably aimed to protect against Germanic invasions and to create a safe corridor for transport and built up of army supplies for the British invasion in 43 A.D. Only later on, probably by the end of the first century, this corridor turned into a frontier zone.

This paper describes the background of five palaeogeographical maps between the Marsdiep and the Weser River, and discusses the natural and anthropogenic processes driving the coastal changes during the last part of the Holocene. Before 2500 BC, during the first half of the Holocene, tidal basins were formed in the lower lying Pleistocene valley system as a result of the Holocene sea-level rise. The tidal basins were filled during the second half of the Holocene and on the deposits from the Pleistocene in the hinterland large coastal peat bogs developed. These peat bogs were vulnerable and sensitive to marine ingressions when the peat surface subsided due to drainage, compaction and erosion. During the Subatlantic (450 BC to present), the different ingression systems in the coastal area between Marsdiep and Weser had their own histories in timing and evolution. The ingressions were naturally caused by lateral migration of coastal barrier and tidal-inlet systems or by changes in the natural drainage system in the hinterland. From the Late Iron Age onwards, humans started to be the major cause of ingressions. By reclaiming and cultivating the seaward margins of coastal peat bogs, these areas subsided significantly and were flooded by high storm surges. When coastal areas were embanked during the historical period, the situation for the lower lying peat lands became more dramatic. When the sea dikes breached, the peat land was flooded, leading to casualties and huge material damages and loss of land. Drowning of the peat lands of the Jade and Dollard in the 14th and 15th centuries are examples of such catastrophes.

The Pliocene and Quaternary unconsolidated sediments of the Upper Rhine Graben (URG) were petrographically analysed in numerous high quality drill cores. The heavy mineral composition of the Graben sediments was compared to those from the Graben margins. In addition, the sedimentary lithofacies were investigated. The chronological classification of the sedimentary successions was established by the interpretation of sporadic palaeontological and palaeomagnetic data.

Within the Pliocene sediments, two distinguishable heavy mineral assemblages indicate different source areas of the Graben fill. At first, a heavy mineral assemblage of stable minerals (turmaline, zircon and anatase) indicates a contribution of debris supplied from Buntsandstein areas at the Graben margins. Secondly, a mixed association of stable minerals in combination with unstable (garnet, hornblende, less epidote) and distinct rare minerals (e.g. monazite, xenotime) can be traced back to debris derived from the crystalline rocks of the southern Graben margins (Black Forest, Vosges). The distribution of sediments with this mixed heavy mineral assemblage proves the fluvial sediment transport from south to north and therefore the course of the Pliocene proto-Rhine along the Graben.

The correlation between the Quaternary sediment successions in the south and the north of the Graben is problematic due to their variable thicknesses as well as their changing lithofacies.

In the southern URG, the Quaternary strata could be subdivided into the older Breisgau Formation and the younger Neuenburg Formation based on characteristic lithofacies. Within this succession, the lower part of the Breisgau Formation (lower Breisgau beds) can be distinguished by noticeable lower contents of hornblende, which probably reflects the effects of weathering and solution of these unstable minerals due to repeated discontinuities during the sediment accumulation. The sediments of the upper part of the Breisgau Formation (upper Breisgau beds) and of the Neuenburg Formation contain a heavy mineral assemblage of garnet, epidote and hornblende, which is typical for Rhine deposits with Alpine contribution. This probably unaffected composition indicates a more unvaried and rapid accumulation of predominantly Alpine debris. In the northern URG, the Quaternary strata are subdivided into three aquifers and intercalated fine-grained horizons. Here, the Quaternary sediments can be petrographically classified into Rhine deposits (garnet, epidote and hornblende) and local accumulations contributed from the Graben margins (turmaline, zircon and anatase) without contemporaneous influence of the Rhine. The analytical results obtained from several drill cores in the northern URG provide evidence for the spatial and temporal variability of the course of the Rhine during the Quaternary.

Time-stratigraphic interpretations of Late Pliocene to Early Pleistocene sediments from onshore locations and from marginal marine settings of the North Sea Basin often refer to the subdivision of the Dutch and British ’Quaternary’ regional stratigraphic stages. Since age control for these stages and their stage boundaries are based on relative dating methods, in this study pollen, dinoflagellate cysts and foraminiferal assemblages were investigated to correlate the regional stratigraphic stages independently to the global chronostratigraphy and the paleomagnetic timescale. The data were obtained from eight boreholes located in the depocentre setting of the Late Pliocene North Sea Basin comprising a 1000 m thick sedimentary succession. The British Gedgravian and Waltonian stages, the Dutch Reuverian to Brunssumian as well as published foraminiferal zones (NSB 14, FB and the lower part of the FA2 zone) fall within the Zanclean and Piacenzian. The lower boundaries of the Pre-Ludhamian and Pretiglian stages and of the NSB 14 to 15 zones are close to the paleomagnetic Gauss-Matuyama boundary. The Pre-Ludhamian, Ludhamian, Thurnian and the Pretiglian, Tiglian A and Tiglian B stages presumably cover the marine isotope stages 103 to 95. It is proposed that the Ludhamian, Thurnian and the Tiglian A were short lasting, warm, periods during which sea level highstand facilitated sedimentary deposition at the marginal areas of the North Sea Basin. The lower boundary of the paleomagnetic Olduvai subchron is situated in the Tiglian Cl-4b stage while the TC4c stage is found within the Olduvai subchron. Foraminiferal NSB 15 and NSB 16 zone as well as the upper part of the FA2 and FA1 zone fall within the Gelasian and cover the Matuyama chron as well as the lower part of the Olduvai subchron. Comparison with formerly dated North Sea sediments shows a good agreement between foraminiferal zonations on a broader scale but significant differences in absolute ages occur. Strontium isotope values indicate approximately 1 Ma younger ages as expected from our chronostratigraphic model. This discrepancy is explained by the dominance of freshwater from river discharge contributing high amounts of eroded material to the basin, leading to an increase of the 87Sr/86Sr ratio in the shelf-sea water.

The north-central Texas Cretaceous section and its contained fossils, as compared to sections at ‘Ein Yabrud in the eastern Mediterranean region, demonstrate that dolichosaurs and primitive mosasaurs inhabited marine environments at least from the intertidal zone to <100 m in depth. The small body size of dolichosaurs and primitive mosasaurs, and the association of Haasiasaurus with marine snakes at ‘Ein Yabrud, suggest similar temperature requirements to modern sea snakes and marine iguanas (Amblyrhynchus cristatus), specifically sea-surface temperatures between approximately 20° and 30° С The Cenomanian and Turonian stages are characterised by widespread shallow seas with relatively high sea-surface temperatures extending to at least 45° N and 65° S latitude. The distribution of dolichosaurs, mosasaurs and snakes during this interval contradicts palaeoclimate models that predict high (>30° C) equatorial sea-surface temperatures and a steep latitudinal temperature gradient in the Northern Hemisphere.

Approximately 200,000 lithological borehole descriptions, 1200 14C dates, 36,000 dated archaeological artifacts, and gradients of palaeochannels were used to reconstruct the Holocene evolution of the fluvial part of the Rhine-Meuse delta. Ages of all Holocene channel belts were stored in a Geographical Information System database that enables generation of palaeogeographic maps for any time during the Holocene. The time resolution of the palaeogeographic reconstruction is about 200 years.

During the Holocene, avulsion was an important process, resulting in frequent shifts of areas of clastic sedimentation. Palaeogeographic evolution and avulsion history of the Rhine-Meuse delta are governed by complex interactions among several factors. These are: (1) Location and shape of the Late Weichselian palaeovalley. In the Early Holocene, rivers were confined to the LateWeichselian valley. When aggradation shifted upstream, the margins of the valley were crossed by newly formed channel belts. (2) Sealevel rise, which resulted in back-filling of the palaeovalley. (3) River channel pattern. In the central-western part of the delta, a straight anastomosed channel pattern with large-scale crevassing developed as a result of sealevel rise and the associated decrease of stream power. (4) Neotectonics. Differential tectonic movements of the Peel Horst and Roer Valley Graben seem to have influenced river behaviour (formation of an asymmetrical meander belt, location of avulsion nodes in fault zones), especially from 4500–2800 14C yr BP when the rate of sealevel rise had decreased. After 2800 14C yr BP sealevel rise further decreased, and tectonic influence still may have influenced avulsions, but from then on other factors became dominant. (5) Increased discharge, sediment load and/or within-channel sedimentation. After 2800 14C yr BP, meander wavelenghts increased, which is interpreted as a result of increased bankfull discharge and/or within channel sedimentation. After 2000 14C yr BP both discharge and sediment load increased as a result of human influence. (6) Coastal configuration. The limited number of tidal inlets and extensive peat formation restricted the number of avulsions in the western part of the delta, and enhanced channel reoccupation. (7) Composition of the substrate and river banks. Meandering river channels tended to adhere to the sandy margins of the LateWeichselian palaeovalley, and high channel sinuosity is found in areas where river banks consisted of sand. Peat formation, which was most extensive in the western part of the back-barrier area especially between 4000 and 3000 14C yr BP, more or less fixed the river pattern at that time, hampering avulsions. (8) Marine ingressions, e.g. the 1421 AD St. Elizabeth’s flood caused large-scale erosion in the southwestern part of the fluvial deltaic plain, resulting in a shift of the main drainage to the SW. (9) Human influence. Since about 1100 AD human influence dominated the palaeogeographic evolution. Rivers were embanked and natural avulsions did no longer occur.

The Wadden Sea is a unique coastal wetland containing an uninterrupted stretch of tidal flats that span a distance of nearly 500km along the North Sea coast from the Netherlands to Denmark. The development of this system is under pressure of climate change and especially the associated acceleration in sea-level rise (SLR). Sustainable management of the system to ensure safety against flooding of the hinterland, to protect the environmental value and to optimise the economic activities in the area requires predictions of the future morphological development.

The Dutch Wadden Sea has been accreting by importing sediment from the ebb-tidal deltas and the North Sea coasts of the barrier islands. The average accretion rate since 1926 has been higher than that of the local relative SLR. The large sediment imports are predominantly caused by the damming of the Zuiderzee and Lauwerszee rather than due to response to this rise in sea level. The intertidal flats in all tidal basins increased in height to compensate for SLR.

The barrier islands, the ebb-tidal deltas and the tidal basins that comprise tidal channels and flats together form a sediment-sharing system. The residual sediment transport between a tidal basin and its ebb-tidal delta through the tidal inlet is influenced by different processes and mechanisms. In the Dutch Wadden Sea, residual flow, tidal asymmetry and dispersion are dominant. The interaction between tidal channels and tidal flats is governed by both tides and waves. The height of the tidal flats is the result of the balance between sand supply by the tide and resuspension by waves.

At present, long-term modelling for evaluating the effects of accelerated SLR mainly relies on aggregated models. These models are used to evaluate the maximum rates of sediment import into the tidal basins in the Dutch Wadden Sea. These maximum rates are compared to the combined scenarios of SLR and extraction-induced subsidence, in order to explore the future state of the Dutch Wadden Sea.

For the near future, up to 2030, the effect of accelerated SLR will be limited and hardly noticeable. Over the long term, by the year 2100, the effect depends on the SLR scenarios. According to the low-end scenario, there will be hardly any effect due to SLR until 2100, whereas according to the high-end scenario the effect will be noticeable already in 2050.

Fluvial sequences, particularly major terrace staircases, represent potential archives of palaeoclimatic fluctuation during the Quaternary. Such sequences can span much if not all of the Quaternary and, provided that dating is possible, can serve as stratigraphical frameworks for correlation with evidence from other depositional environments. In particular, they can provide a terrestrial lithostratigraphical framework that can be correlated with the global marine oxygen isotope record. Fluvial lithos-tratigraphical frameworks also provide important contexts for records of faunal evolution and human occupation, the latter largely determined from the occurrence of artefacts in fluviatile sediments.

This paper announces a new project within the International Geological Correlation Programme, devoted to fluvial sequences (IGCP 449 ‘Global Correlation of Late Cenozoic fluvial deposits’). It attempts to summarize existing baseline knowledge at the outset of the project and outlines the proposed methods and criteria for establishing a database of fluvial sequences.

Over the past 80 years, the Rhine-Meuse delta has been mapped extensively by drilling boreholes. These maps are compared with a new detailed digital elevation map of the Netherlands (AHN) that became available in 2004 and is based on very accurate, (sub-decimeter) laser-altimetry data. Examples show, that existing maps can be significantly improved. However, field checks remain a necessity. Therefore, the AHN-data need to be confronted with existing borehole descriptions using a GIS. Geomorphological analysis and geological interpretation of surface elevation patterns now enable us to map larger areas, in greater detail, with greater accuracy, and much faster. It is argued that detailed (≥1 : 50,000) geological mapping of the Netherlands should be resumed, combining the databases of Utrecht University, the National Geological Survey and Alterra.

Information about the morphology and alluvial sediments of the Dnieper Valley is reviewed. The Dnieper Valley originated in the Late Miocene. The Middle Dnieper Valley is an intercontinental alluvial basin and the Lower Dnieper Valley is a shallow canyon that ends with a delta. Identification of the alluvial dynamic facies (channel, overbank, abandoned channel) is crucial for stratigraphical analysis. The dynamic facies form regular sequences - alluvial suites that combine into series. Individual suites and series are characterized by their mode of occurrence, facies composition, lithological features and expression in the modern landscape. Their stratigraphie position is established with reference to index beds and palaeontological, geochrono-logical and archaeological research, allowing them to be correlated along the valley. Correlation between different parts of the Dnieper system uses a combination of facies and geomorphological analyses, whereas correlation with other river systems makes use of mammalian and molluscan biostratigraphy. Global events (marine regressions and glaciations) that brought about similar reactions in different fluvial systems can be also used for correlation.

The marine Eemian (marine oxygen-isotope substage 5e: MIS 5e) is represented by shallow-water deposits in southern and western Denmark, while relatively deep-water environments occurred to the north and north-east, where complete interglacial successions seem to be present. We present an overview of the marine Eemian deposits in Denmark, and discuss in more detail indications of climate variability, both for the late Saalian and within the Eemian.

An explanation is presented for the origin of brackish to saline groundwater in the coastal area of the Netherlands based on geological, chemical (chlorinity), isotopic and geophysical data. A critical review of all possible salinization mechanisms shows that the origin of the brackish water is related to former transgressions. Both the vertical salinity distribution and the carbon-14 activity of the groundwater indicate that connate sea water from the Pliocene to Early Pleistocene is not the source of the brackish to saline waters in the overlying Pleistocene fluvial aquifers. Instead, it derives from Holocene transgressions. The salinization mechanism is discussed in relation to the paleogeographical development during the Holocene and the occurrence of low-permeability strata. Finally, freshening of the aquifers following retreat of the sea is briefly considered.

We present a revised relative mean sea-level (MSL) curve for the Rhine-Meuse delta, western Netherlands, for the period 7900-5300 cal yr BP. The revision is based on a series of new and previously unpublished local groundwater-level index data from buried Late Glacial aeolian dunes in the lower Rhine-Meuse delta, and reinterpretation of existing data.

The new index data consist of (AMS and conventional) radiocarbon dates of samples, collected from the base of peat formed on dune slopes, near Vlaardingen (21 index points), Hillegersberg (one index point), and Hardinxveld-Giessendam (10 index points). The Vlaardingen data represent the coast-nearest Rhine-Meuse delta local water-level record, which therefore is highly indicative for sea-level change. Pollen and macrofossil analysis, and dating of paired samples was carried out to assess the reliability of the groundwater-level index data.

The revision of the MSL curve involves: (1) a significant (0 to >1 m) upward adjustment for the period 7900-7300 cal yr BP; (2) a downward adjustment of <0.25 m for the period 6650-5300 cal yr BP. The new data indirectly support the reliability of the part of the curve for the period 7300-6650 cal yr BP. A longitudinally fairly uniform river gradient of 2.5-3.0 cm/km in the lower Rhine-Meuse delta during the period 6650-5600 cal yr BP can be inferred from the data sets. A significant river gradient extended further towards the coastline than previously thought and it may be that also the revised MSL curve reflects river-gradient effects. An increased floodbasin effect (stronger intra-coastal tidal damping) seems to have developed in the lower Rhine-Meuse delta in de period 7500-6600 cal yr BP, and was probably a complex response to a major avulsion of the Rhine.