Late Pleistocene evolution of the Rhine-Meuse system in the southern North Sea basin: imprints of climate change, sea-level oscillation and glacio-isostacy

https://doi.org/10.1016/j.quascirev.2007.07.013Get rights and content

Abstract

High-resolution continuous core material, geophysical measurements, and hundreds of archived core descriptions enabled to identify 13 Late Pleistocene Rhine-Meuse sedimentary units in the infill of the southern part of the North Sea basin (the Netherlands, northwestern Europe). This sediment record and a large set of Optical Stimulated Luminescence dates, 14C dates and biostratigraphical data, allowed to establish detailed relationships between climate change, sea-level oscillation, glaciation history and the sedimentary development of the Rhine fluvial system during the last glacial cycle (Marine Isotope Stages 5e-2, Eemian-Weichselian). A well-preserved Eemian sediment record was encountered as the infill of a Late Saalian (MIS6) subglacial basin. Part of this record reflects groundwater rise controlled (fine-grained) sedimentation as a result of postglacial (early) Eemian sea-level rise. It shows strong analogy to developments known from the Holocene Rhine-Meuse delta. Outside of the glacial depressions near coastal deposits are only fragmentarily preserved. The Early Glacial Rhine sediment record is dominated by organic debris and peat layers, marking landscape stability and low fluvial activity. Part of this record may have been formed under near coastal conditions. Significant amounts of reworked marine biomarkers in the lag-deposits of Early Pleniglacial (MIS4) fluvial systems indicate that this period is characterized by extensive reworking of older (MIS5) near-coastal sediments. Despite the marked Early Pleniglacial climatic cooling, input of new sediment from the drainage basin was relatively low, a feature that is related to the presence of regolith protective relic soil complexes in the basin. During the early Middle Pleniglacial, a major Rhine avulsion indicates the system was in an aggrading mode and that sediment supply into the lower reaches of the Rhine had strongly increased. This increase in sediment supply coincided with the timing of major climate cooling that occurred from ∼50 to 45 ka onwards. The increase in sediment supply is related to final breakup of the soil complexes in the drainage basin. After ∼24 ka, a strong input of coarse-grained gravelly sediments was observed which indicates a strong increase in physical weathering processes and periglacial-controlled supply of bedload sediment in the catchment. A time delay between climate change (∼30 ka) and channel belt aggradation (<24 ka), is explained as a result of transport path length between source and sink and/or effects of higher continental runoff rates after 22 ka. The Late Middle Pleniglacial, Late Pleniglacial and Lateglacial Rhine-Meuse record testifies for strong influence of glacio-isostatic-controlled differential upwarping of the study area. Glacio-isostatic-controlled forebulge upwarping and lateral valley tilting is shown to have deflected Rhine-Meuse channel belts after 35 ka. Glacio-isostatic upwarping is seen as the main cause for strong incision during the first phase of the Late Pleniglacial (30–24 ka). At later stage glacio-isostatic-controlled incision was overruled due to high climate-controlled sediment input from the catchment and probably initial glacio-isostatic subsidence. Migration of channel belts towards the direction of the former centre of glacio-isostatic uplift indicates that glacio-isostacy influenced Rhine-Meuse paleogeography until far into the Lateglacial.

Introduction

Rivers are transfer systems between areas of sediment production and sediment storage and therefore have an important role in long-term sediment re-organization across the Earth's surface. Insights into the sensitivity of the fluvial exchange system under external forcings are essential for understanding present and past landscape stability, for improvement of knowledge on basin-marginal sediment records and for predicting future change. Late Pleistocene fluvial records in northwestern Europe constitute important archives for improving the understanding of the response of fluvial systems to external forcing. Well-dated ice-core, marine and terrestrial records (Guiot et al., 1989; Sidall et al., 2003; NGRIP Members, 2004; Ménot et al., 2006) show that the rivers from this period experienced periodical, high amplitude changes in climate and sea-level while the Scandinavian and British ice-masses expanded and disintegrated multiple times (Clark et al., 2004; Mangerud, 2004; Svendssen et al., 2004). A long history of fluvial response to external forcing is registered in the record of the Rhine river system in northwestern Europe (Fig. 1). Its Late Pleistocene deposits are relatively well preserved and within dating range of several techniques enabling correlation between records of response and independently dated (proxy) records of forcing. The Rhine is a ∼1300 km long, rain- and snowmelt-fed river system that drains 185,000 km2 of the central and northwest European foreland towards the North Sea basin, which saw the southern limit of Scandinavian and British ice multiple times in the Middle and Late Pleistocene.

During the Late Pleistocene, considerable quantities of primarily coarse-grained sediment were produced in the upland areas of the catchment that includes the Alps, Vosges, Black Forest and uplands of the Rhenish Shield (Fig. 1). A major part of that sediment was trapped in the south-eastern part of the North Sea sedimentary basin (Fig. 1). It forms a record consisting of 10–25 m of primarily coarse-grained gravely sands in the subsurface of the Netherlands (Fig. 2, Kreftenheye Formation, cf. Zonneveld, 1947; Doppert et al., 1975; Westerhoff et al., 2003), that within the North Sea basin is well-preserved owing to burial by the Holocene coastal prism (e.g. Berendsen and Stouthamer, 2000) and owing to the fact that Late Pleistocene fluvial units are much wider than Holocene (deltaic) channel belts in the area, that only locally caused erosion.

This record reflects the infill of the southern North Sea basin of the last interglacial–glacial cycle and is well preserved due to background tectonic subsidence, relatively unconstrained valley widths and occurrence of deeply scoured glacial basins that were formed during the Penultimate (Saalian) glaciation (Fig. 2). The record was formed under strongly varying climate (Fig. 3) including periods of full interglacial conditions (Eemian) and times of alternating stadial–interstadial conditions (Weichselian) with superimposed rapid oscillations at shifts between colder and warmer states (Voelker, 2002).

Eustatic sea-level oscillations of tens of meters in amplitude occurred repeatedly during the Late Pleistocene, but were interspersed with periods of relative sea-level stability (Fig. 3). The ultimate low stand was over 100 m below modern MSL some 25–20 ka BP. Sea-level oscillation led to repetitive extension and truncation of the river Rhine across the North Sea floor and beyond through the Dover Strait (Bridgland and D’Olier, 1995; Gibbard, 1995; Antoine et al., 2003; Fig. 1). In addition, the lower reaches of the Rhine experienced significant glacio-isostasy-driven crustal warping (forebulge buildup and collapse) because of its particular distance to centers of ice-sheet formation in Scandinavia, the Baltic and Britain (Mörner, 1979; Peltier, 1994, Peltier, 2004; Lambeck, 1995; Kiden et al., 2002; Lambeck et al., 2006).

Recently, a series of publications have presented new data on the sedimentary buildup and chronology of the Late Pleistocene Rhine system (Törnqvist et al., 2000, Törnqvist et al., 2003; Wallinga et al., 2004; Busschers et al., 2005). These papers documented the first applications of OSL-dating on the Rhine sediments, carried out in a limited study area in the west-central Netherlands. Their discussion sections focused on the role of sea-level change and related sequence stratigraphic interpretations (Törnqvist et al., 2000, Törnqvist et al., 2003) and later also on evaluating the relative roles of coeval forcing factors during the last interglacial–glacial cycle (Wallinga et al., 2004; Busschers et al., 2005). The major part of the Late Pleistocene fluvial deposits, however, remained to be studied in similar detail. Such work was a necessity to validate the many ideas on the bearings of the Rhine-record, and the accumulated data presented in this paper first allows doing so.

The present paper evaluates the imprints of northwest European climate-change, sea-level oscillations and glaciation history on the sedimentary record of the Rhine fluvial system in the southern part of the North Sea basin. We focus on the Rhine record of the last glacial cycle (Eemian-Weichselian-MIS5e-2, cf. Bassinot et al., 1994) that was deposited in the Netherlands and adjacent offshore area. We make use of new continuous core–material and numerous digitally archived core descriptions from the DINO-database of TNO Built Environment and Geosciences, Geological Survey of the Netherlands (TNO, 2006). This data is used for facies analysis (grain-size, sedimentary structures and mineralogy) and construction of geological transects and palaeogeographical maps. Chronological control is obtained by 94 quartz OSL-dates (optically stimulated luminescence) and by biostratigraphy (pollen zonation schemes). This enabled well-constrained reconstruction of sedimentation in the downstream reaches of the Late Pleistocene Rhine valley(s) and allows for explaining sediment composition variations, shifting loci of sediment preservation, variation in channel belt geometry, bounding surface formation and palaeo-geographical evolution by causal linkage to records of external forcing.

Section snippets

Geographical and geological setting

The Rhine fluvial system at present drains ∼185,000 km2 of Alpine, central and northwestern Europe before it enters the present North Sea near Rotterdam. Along its pathway from headwaters in the Swiss Alps to the North Sea, the Rhine is joined by several large tributaries of which the Aare, Mosel, Neckar, Main and Meuse are most important (Fig. 1). Modern Rhine discharges measured at the Dutch–German border vary from 2300 m3/s (mean annual discharge) to a maximum of 12,000 m3/s following extreme

Data collection and sampling strategies

The sequence is described by means of twenty 20–40-m deep continuous cores (Busschers et al., 2005; this study) and hundreds of archived core descriptions (of counter-flush and air-lift drillings), lacquer peels and associated material archived by the Geological Survey (TNO, 2006). Four valley-wide transects perpendicular to the main palaeo flow direction of the Rhine (Fig. 2) document complete stratigraphic sequences produced by the Late Pleistocene Rhine-Meuse river system. For this study, 11

Facies descriptions and interpretations of Area A

The sequence of Area A (10–15 m thick) consists of five amalgamated sedimentary units that comprise fluvial channel-belts and floodbasin deposits. The geometry and stratigraphy of the units of Area A is depicted for a valley-wide transect (Fig. 4). That sequence buries coarse-grained fluvial sands and glacio-lacustrine clays of Saalian age and is subdivided into five sedimentary units that vary 12.5–35 km in width and 2–8 m in thickness (Table 1). Fig. 5 shows the sedimentological logs of cores in

Facies descriptions and interpretations of Area B

The sequence in Area B (10–20 m thick) consists of seven sedimentary units that largely represent amalgamated fluvial channel-belt deposits. The stratigraphical position and geometry of the units in Area B are depicted in three valley-wide cross-sections for the upstream (Fig. 6(A)), central (Fig. 6(B)) and downstream part (Fig. 6(C)) of Area B respectively. The sediments cover coarse-grained deposits of Early and Middle Pleistocene age, including deposits associated to the maximum Saalian ice

Chronostratigraphy

Age control on the oldest part of the Rhine record, i.e. the Eemian and Early Weichselian sequence in Area A, was obtained by pollen stratigraphy complemented by two OSL dates. Twelve OSL-dates provide chronological control for the Pleniglacial and younger sequence in Area A. In Area B, chronological control is based on 80 OSL-dates and some 14C dates, while no (useful) pollen records were encountered here.

Eemian—Early Glacial/MIS5 (Units A1, A2, A3, B1)

During the disintegration of the Saalian ice-sheet, the Rhine established a course through the IJssel glacial basin (Area A) and the northern part of the Netherlands (Fig. 9(A)). Subsequent deposition of glacio-lacustrine and fluvial sediments caused the basin to be largely filled up at the onset of the Eemian (Van de Meene and Zagwijn, 1978). In Area A, we observed multiple deep channels at the top of this Saalian sequence. Pollen-assemblage registration of the Zeifen-Kattegat climate

Fluvial evolution in relation to glacio-isostacy

Glacial reconstructions of Scandinavia and Britain show that during the Weichselian, multiple large-scale ice-sheet extensions occurred of which the Late Pleniglacial glaciation is the most recent one and of best-known extent (Fig. 1, Fig. 2). Relative sea-level rise reconstructions for Lateglacial to Middle Holocene times testify for the occurrence of anomalous high non-linear subsidence rates along the southern North Sea coasts while the deglaciated Scottish and Norwegian coasts show strong

Conclusions

Our large study area and multidisciplinary approach has yielded new data and advanced insight in the Late Pleistocene infill of the Southern North Sea basin. Integration of high-resolution core material, archived core data, OSL-dating and pollen analysis has delivered a detailed reconstruction that documents the development of a major river system over a full glacial cycle. The reconstruction enabled to determine relative controls of climate change, base-level variation and glaciation on

Acknowledgements

This paper is part of the Ph.D. study of Freek S. Busschers at the Vrije Universiteit Amsterdam (VU). The project was financed by TNO Built Environment and Geosciences—Geological Survey of the Netherlands (TNO-B&O) and by the Dutch Ministry of Transport, Public Works and Water Management, Directorate General for Water Affairs (Grant DWW-1804). The Netherlands Centre for Luminescence dating (NCL) is supported by the Netherlands Organization for Scientific research (NWO; Grants 863.03.006 and

References (121)

  • M. Houmark-Nielsen

    The last interglacial/glacial cycle in Denmark

    Quaternary International

    (1989)
  • J. Holbrook et al.

    Geomorphic and sedimentary response of rivers to tectonic deformation: a brief review and critique of a tool for recognizing subtle epeirogenic deformation in modern and ancient settings

    Tectonophysics

    (1999)
  • R.F.B. Isarin et al.

    Mean July temperatures during the younger dryas in Northern and Central Europe as inferred from climate indicator plant species

    Quaternary Research

    (1999)
  • H. Kooi et al.

    Geological causes of recent (∼100 yr) vertical land movement in the Netherlands

    Tectonophysics

    (1998)
  • J. Mangerud

    Ice sheet limits on Norway and the Norwegian continental shelf

  • A.S. Murray et al.

    Luminescence dating of quartz using an improved single-aliquot regenerative-dose protocol

    Radiation Measurements

    (2000)
  • A.S. Murray et al.

    The single aliquot regenerative dose protocol: potential for improvements in reliability

    Radiation Measurements

    (2003)
  • H.S. Schmincke et al.

    Evolution and environmental impacts of the eruption of Laacher See Volcano (Germany) 12,900 a BP

    Quaternary International

    (1999)
  • J. Schokker et al.

    An OSL dated Middle and Late Quaternary sedimentary record in the Roer Valley Graben (southeastern Netherlands)

    Quaternary Science Reviews

    (2005)
  • M.-S. Seidenkrantz et al.

    Two-step deglaciation at the Oxygen Isotope Stage 6/5E transition: the Zeifen-Kattegat climate scillation

    Quaternary Science Reviews

    (1996)
  • R.L. Skelly et al.

    Architecture of channel-belt deposits in an aggrading shallow sandbed braided river: the lower Niobrara River, northeast Nebraska

    Sedimentary Geology

    (2003)
  • T.E. Törnqvist et al.

    Response of the Rhine-Meuse system (west-central Netherlands) to the last Quaternary glacio-eustatic cycles: a first assessment

    Global and Planetary Change

    (2000)
  • R.T. Van Balen et al.

    Neotectonics of The Netherlands: a review

    Quaternary Science Reviews

    (2005)
  • J. Vandenberghe

    Paleoenvironment and stratigraphy during the Last Glacial in the Belgian–Dutch border region

    Quaternary Research

    (1985)
  • J. Vandenberghe

    Timescales, climate and river development

    Quaternary Science Reviews

    (1995)
  • G. Aalbersberg et al.

    Multiproxy climate reconstruction for the Eemian and Early Weichselian

    Journal of Quaternary Science

    (1998)
  • L. Ahorner

    Untersuchungen zur quartären Bruchschollentektonik der Niederrheinischen Bucht

    Eiszeitalter und Gegenwart

    (1962)
  • P. Antoine et al.

    The Pleistocene rivers of the English Channel region

    Journal of Quaternary Science

    (2003)
  • M. Ballarini et al.

    Spatial variation of dose rate from beta sources as measured using single grains

    Ancient TL

    (2006)
  • M.D. Bateman et al.

    The timing of Lastglacial periglacial and aeolian events, Twente, eastern Netherlands

    Journal of Quaternary Science

    (1999)
  • H.J.A. Berendsen et al.

    Late Weichselian and Holocene river channel changes of the rivers Rhine and Meuse in the central Netherlands

    Paläoklimaforschung

    (1995)
  • Berendsen, H.J.A., Stouthamer, E., 2001. Palaeogeographic Development of the Rhine–Meuse Delta, The Netherlands, Van...
  • S.J. Blott et al.

    Technical communication: GRADISTAT: a grain size distribution and statistics package for the analysis of unconsolidated sediments

    Earth Surface Processes and Landforms

    (2001)
  • W. Boenigk

    Die flußgeschichtliche Entwicklung der Niederrheinischen Bucht im Jungtertiär und Altquartär

    Eiszeitalter und Gegenwart

    (1978)
  • W. Boenigk

    The Pleistocene drainage pattern in the Lower Rhine basin

    Netherlands Journal of Geosciences

    (2002)
  • J.A.A. Bos et al.

    Vegetation and climate during the Weichselian Early Glacial and Pleniglacial in the Niederlausitz, eastern Germany—macrofossil and pollen evidence

    Journal of Quaternary Science

    (2001)
  • J. Bridge

    Rivers and Floodplains—Forms, Processes and Sedimentary Record

    (2003)
  • Bridgland, D.R., D’Olier, B., 1995. The Pleistocene evolution of the Thames and Rhine drainage systems in the southern...
  • F.S. Busschers et al.

    Sedimentary architecture and optical dating of Middle and Late Pleistocene Rhine-Meuse deposits-fluvial response to climate change, sea-level fluctuation and glaciation

    Netherlands Journal of Geosciences

    (2005)
  • F.S. Busschers

    Unraveling the Rhine, response of a fluvial system to climate change, sea-level oscillation and glaciation. Published Ph.D. dissertation, Vrije Universiteit Amsterdam/TNO Built Environment and Geosciences

    Geology of the Netherlands

    (2008)
  • S.J. Carr et al.

    The Last Glacial Maximum in the North Sea basin: micromorphological evidence of extensive glaciation

    Journal of Quaternary Science

    (2006)
  • G. Caspers et al.

    Vegetation and climate in the Early and Pleni-Weichselian in northern central Europe

    Journal of Quaternary Science

    (2001)
  • P. Cleveringa et al.

    The Eemian stratotype locality at Amersfoort in the central Netherlands: a re-evaluation of old and new data

    Geologie en Mijnbouw

    (2000)
  • Cohen, K.M., 2003. Differential subsidence within a coastal prism. Lateglacial–Holocene tectonics in the Rhine-Meuse...
  • K.M. Cohen

    3D Geostatistical interpolation and geological interpretation of palaeo-groundwater rise in the Holocene coastal prism in the Netherlands

  • K.M. Cohen et al.

    Fluvio-deltaic floodbasin deposits recording differential subsidence within a coastal prism (central Rhine-Meuse delta, the Netherlands)

  • W. De Gans et al.

    Late Saalian and Eemian deposits in the Amsterdam glacial basin

    Geologie en Mijnbouw

    (2000)
  • J.W.Chr. Doppert et al.

    Formaties van het Kwartair en Boven-Tertiair in Nederland

  • Ehlers, J., Gibbard, P.L., 2004. (Eds.) Quaternary Glaciations: Extent and Chronology, vol. 1, Europe. Developments in...
  • D. Ellwanger

    Eine landschaftsübergreifende Lockergesteinsgliederung vom Alpenrand zum Oberrhein

  • Cited by (227)

    • Distribution of rare earth elements and yttrium in water, suspended matter and bottom sediments in Lake Onego: Evidence of the watershed transformation in the Late Pleistocene

      2023, Quaternary International
      Citation Excerpt :

      Glacial-isostatic changes in the local topography caused by glacial retreat had a substantial impact on the inclination of the terrain, which after the glacial retreat accounted for 30 m/100 km comparing with contemporary (Demidov, 2006). The gradual uplift of the territory after the last glaciation could also affect the drainage routes of rivers (Busschers et al., 2007; Brandes et al., 2011; Panin et al., 2020), which in turn could be reflected in the watershed transformation and changes in sedimentation environment of Lake Onego. Since the composition and geochemistry of the bottom sediments of Lake Onego are generally controlled by two major sources: the Karelian craton (the western and northwestern, eastern parts of the lake drainage areas) and the southern part - the rocks of RP, it is possible to use the data on the REE + Y to obtain independent additional information of the changes in sedimentation environment of the lake in the Late Pleistocene and Holocene.

    View all citing articles on Scopus
    View full text