Introduction
The relation(s) between fluvial systems and climate and vegetation for the Late-glacial and Early Holocene have long been studied in the Netherlands (Pons, Reference Pons1957; Van den Broek & Maarleveld, Reference Van den Broek and Maarleveld1963; Berendsen et al., Reference Berendsen, Hoek, Schorn and Frenzel1995; Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995; Huisink, Reference Huisink1997; Tebbens et al., Reference Tebbens, Veldkamp, Westerhoff and Kroonenberg1999; Van Huissteden & Kasse, Reference Van Huissteden and Kasse2001; Cohen et al., Reference Cohen, Stouthamer and Berendsen2002; Cohen, Reference Cohen2003; Busschers et al., Reference Busschers, Kasse, Van Balen, Vandenberghe, Cohen, Weerts, Wallinga, Johns, Cleveringa and Bunnik2007; Hij-ma et al., Reference Hijma, Cohen, Hoffmann, Van der Spek and Stouthamer2009; Erkens et al., Reference Erkens, Hoffmann, Gerlach and Klostermann2011; Janssens et al., Reference Janssens, Kasse, Bohncke, Greaves, Cohen, Wallinga and Hoek2012; Zuidhoff & Huizer, Reference Zuidhoff and Huizer2015). In these studies braided systems are correlated to cold periods, whereas meandering systems are mostly connected to interstadial and interglacial conditions. The Meuse river and Meuse–Niers–Rhine confluence area between Roermond and Nijmegen revealed an especially strong connection between climate, channel planform changes and phases of instability as demonstrated by incision and terrace formation (Pons, Reference Pons1957; Berendsen et al., Reference Berendsen, Hoek, Schorn and Frenzel1995; Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995; Huisink, Reference Huisink1997; Tebbens et al., Reference Tebbens, Veldkamp, Westerhoff and Kroonenberg1999). Subsequently, river systems in Poland, Germany and Hungary have been investigated in order to establish whether the changes to those systems were similar and synchronous to those documented in the Netherlands, to determine local (vegetation, tectonics) or regional climate forcing (Vandenberghe et al., Reference Vandenberghe, Kasse, Bohncke and Kozarski1994; Kasse et al., Reference Kasse, Hoek, Bohncke, Konert, Weijers, Cassee and Van der Zee2005, Reference Kasse, Bohncke, Vandenberghe and Gábris2010; Kiss et al., Reference Kiss, Hernesz, Sümeghy, Györgyövics and Sipos2015; Panin & Matlakhova, Reference Panin and Matlakhova2015; Starkel et al., Reference Starkel, Michczyńska, Gębica, Kiss, Panin and Perşoiu2015). It was established that the major climate shift from the Pleniglacial to the Late-glacial was reflected, with variable delay, regionally in river systems by a change from braided to meandering. Several studies have identified multiple parallel functioning meandering channels during braided to meandering transitional stages at the Pleniglacial to Late-glacial and Younger Dryas to Holocene transitions (Vandenberghe et al., Reference Vandenberghe, Kasse, Bohncke and Kozarski1994; Erkens et al., Reference Erkens, Hoffmann, Gerlach and Klostermann2011). However, short-lived climate events were recorded differently. For example, the response of river systems to the Younger Dryas cooling differed regionally and was not always reflected by meandering to braided channel planform changes, as in the Meuse and Rhine valleys. The Warta river in Poland (Vandenberghe et al., Reference Vandenberghe, Kasse, Bohncke and Kozarski1994), the Niers in Germany (Kasse et al., Reference Kasse, Hoek, Bohncke, Konert, Weijers, Cassee and Van der Zee2005) and the Tisza in Hungary (Kasse et al., Reference Kasse, Bohncke, Vandenberghe and Gábris2010) persisted in a meandering mode during the Younger Dryas, and it was therefore concluded that the threshold towards braiding was not always crossed due to differences in gradient, stream power, lithology and bank stability. Even intra-catchment variable response to external forcing was established in the Wetter valley in Germany (Houben, Reference Houben2003).
In addition to climate-related planform changes of river systems, climate change may also be reflected by short phases of fluvial instability, with incision and terrace formation. For this reason former floodplains and palaeochannels may be preserved and their geomorphology, lithology and palynology can be studied. At both cold–warm and warm–cold climate transitions changes in precipitation, evapotranspiration, river discharge and soil cohesion result in changes in the water/sediment ratio which can lead to incision and terrace formation (Vandenberghe, Reference Vandenberghe1995, Reference Vandenberghe2008). However, in most fluvial terrace staircase records the products of cold–warm and warm–cold transitions are poorly preserved due to widespread lateral erosion and valley widening by braided systems in the next cold phase. In many cases, only the fluvial archive related to the last glacial–interglacial transition has been preserved, as is the case in the Roer valley, enabling reconstructions and dating of glacial–interglacial fluvial-systems changes. Although the main rivers Rhine and Meuse have received most attention in the Netherlands (Van den Broek & Maarleveld, Reference Van den Broek and Maarleveld1963; Berendsen et al., Reference Berendsen, Hoek, Schorn and Frenzel1995; Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995; Huisink, Reference Huisink1997; Tebbens et al., Reference Tebbens, Veldkamp, Westerhoff and Kroonenberg1999; Van Huissteden & Kasse, Reference Van Huissteden and Kasse2001; Cohen, Reference Cohen2003; Busschers et al., Reference Busschers, Kasse, Van Balen, Vandenberghe, Cohen, Weerts, Wallinga, Johns, Cleveringa and Bunnik2007; Hijma et al., Reference Hijma, Cohen, Hoffmann, Van der Spek and Stouthamer2009; Erkens et al., Reference Erkens, Hoffmann, Gerlach and Klostermann2011), their tributaries have been studied in some detail as well. From the reconstructed fluvial development of the Geul (De Moor et al., Reference De Moor, Kasse, Van Balen, Vandenberghe and Wallinga2008), Niers (Kasse et al., Reference Kasse, Hoek, Bohncke, Konert, Weijers, Cassee and Van der Zee2005), Lippe and Oude IJssel (Janssens et al., Reference Janssens, Kasse, Bohncke, Greaves, Cohen, Wallinga and Hoek2012) it was concluded that river system changes in the tributaries differ in some aspects from those seen in the main rivers. In this study we focus on the Roer river, a tributary of the Meuse. Despite its moderate catchment size, it has so far received little attention (Van den Broek & Van der Marel, Reference Van den Broek and Van der Marel1964; Schalich, Reference Schalich1968; Van Zuidam, Reference Van Zuidam1980; Van den Berg, Reference Van den Berg1989; Bakels, Reference Bakels2017). The age control on the Late-glacial and Holocene development is still poor (Hoek, Reference Hoek1997b: sites RMO3, RMO4, Melick). Attempts have been made to fill this knowledge gap in several bachelor and master research projects (Van Bostelen & Vreugdenhil, Reference Van Bostelen and Vreugdenhil2009; Janssens, Reference Janssens2011; Vreugdenhil, Reference Vreugdenhil2011).
The aims of this study are to: (1) reconstruct the Late Pleniglacial, Late-glacial and Early Holocene fluvial channel pattern changes, incision events and terrace formation in the lower Roer valley; (2) establish the impact of climate change on the vegetation and the fluvial system; and (3) explore the role of the Meuse as a local base level for the Roer.
Study site
The Roer (Rur in Germany) catchment (2354 km2) is situated in the southern Netherlands and adjacent parts of Germany and Belgium (Figs 1 and 2). The Roer has a length of 165 km and its source in the Hautes Fagnes in Belgium is at an altitude of 660 m a.s.l. The river debouches into the Meuse near Roermond at an altitude of c. 19 m a.s.l. The mean annual discharge (AD 1969–2009) at Stahe, c. 7 km upstream of the Dutch–German border in Germany, is 21.8 m3 s−1 and the discharge varies between 8 and 124 m3 s−1. The Roer river is located in the SE–NW-oriented Roer Valley Graben bounded by the Peel Boundary Fault zone in the northeast and the Feldbiss Fault zone in the southwest (Fig. 1). The current extension phase of the graben started at the Oligocene–Miocene transition (Van Balen et al., Reference Van Balen, Houtgast, Van der Wateren, Vandenberghe and Bogaart2000, Reference Van Balen, Houtgast and Cloetingh2005). The graben was filled with Rhine and Meuse river deposits during the Early and Middle Pleistocene until ~400 ka, and with local fluvial and aeolian sediments after that time (Schokker et al., Reference Schokker, Cleveringa, Murray, Wallinga and Westerhoff2005). The Roer Valley Graben is tectonically active, as demonstrated by c. 1 m displacement of Allerød-age terraces and Younger Dryas dunes along the Peel Boundary Fault (Michon & Van Balen, Reference Michon and Van Balen2005). Despite long-term subsidence of the Roer Valley Graben relative to adjacent areas, well-developed terraces and palaeochannels have formed in the confluence area of the Roer and Meuse during the glacial–interglacial transition, related to climate change.
Location map and tectonic setting of the study area in the southeastern Netherlands and western Germany. RVG = Roer Valley Graben.
Digital elevation map of the Roer river catchment. Coordinates according to the Dutch Ordnance System. For location see Figure 1.
Methods
The geomorphological analysis was based on topographical maps (1 : 25,000), digital elevation models of the Netherlands (AHN) and Germany and morphological mapping (Figs 3, 4 and 5). Morphogenetic units were distinguished based on surface morphology and elevation, location of terrace scarps, palaeochannel planform characteristics and cross-cutting relationships (Fig. 3b).
(A) Digital elevation map of the lower reaches of the Roer; (B) morphogenetic interpretation of the study area with the terraces and floodplain levels, aeolian deflation and accumulation areas; (C) locations mentioned in the text. For location see Figure 2.
Digital elevation map of Late-glacial level B and Holocene floodplain level C in the downstream part of the Roer valley. Position of cross section AAʹ, corings (open dots) and OSL-sample locations (closed dots) and ages in ka are indicated. For location see Figure 3.
Digital elevation map of Late-glacial level B in the Roer valley at the Dutch–German border. Position of cross section BBʹ, corings (open dots) and locations of pollen diagrams from abandoned channels (green dots) are indicated (a=Karken; b=Geraerds; c=Haaserdriesch; d=Kitscherholz; e= Kitscherholz b). For location see Figure 3.
Coring transects and lithological description
Two main cross sections were cored perpendicular to the valley axis to investigate the lithology and to reconstruct the fluvial architecture of the area (Fig. 6). Five detailed lithological cross sections were made across abandoned channels to find the deepest and oldest infill for pollen analysis, to investigate vegetation development and to date the timing of channel abandonment (Figs 5 and 7). Coring was performed with an Edelman auger, a gauge and Van der Staay suction-corer. Sediment description in the field was done according to the Dutch classification of unconsolidated soil samples (NEN-5104). In this paper, textural classes according to the United States Department of Agriculture (USDA) have been used. Sampling for pollen analyses from the fills of abandoned channels was done with a gauge or piston corer (diameter 6 cm).
Lithogenetic cross sections over abandoned channels of level B in the Dutch–German border area of the Roer valley. Asterisks indicate positions of pollen diagrams. For location see Figure 5.
Pollen analyses
Pollen analysis was performed to reconstruct the vegetation development (Figs 8, 9, 10 and 11). Local pollen assemblage zones, based on changes in species composition, were correlated to the existing biochronostratigraphic zonation for the Late-glacial and the Early Holocene in the Netherlands (Hoek, 1997a). This provides an indirect dating tool for fluvial system change. The cores from the palaeochannels were described in detail in the laboratory. A total of 48 samples from four cores were selected, based on lithology and lithological transitions. The pollen samples were prepared according to standard procedures (Faegri & Iversen, Reference Faegri and Iversen1989). The pollen slides were examined using a Zeiss Axioskop 50 light microscope with a magnification of 630× with and without phase contrast. Pollen determination was conducted using the pollen key of Moore et al. (Reference Moore, Webb and Collinson1991). The pollen sum varied between 100 and 200 and included tree taxa, shrubs and upland herbs. Wetland herbs, aquatics and Cyperaceae were excluded from the pollen sum. Pollen diagrams were constructed using TILIA software in which taxa are arranged in ecological groups (Grimm, Reference Grimm1992).
Pollen diagram Karken. For location see Figure 5.
Pollen diagram Geraerds. For location see Figure 5.
Pollen diagram Haaserdriesch. For location see Figure 5.
Pollen diagram Kitscherholz. For location see Figure 5.
Optically stimulated luminescence dating
Fluvial and aeolian sands were sampled for luminescence dating to reconstruct meandering channel dynamics, and to date incision events and aeolian dune formation (Fig. 4; Table 1). Boreholes were made with an Edelman corer (10 cm diameter), and the disturbed and bioturbated soil and fine-grained overbank sediments were removed down to the sampling depth. Sampling was done in the boreholes by vertically hammering 40 cm long stainless steel tubes with a diameter of 7 cm into the sand (Van der Horst coring equipment). The optically stimulated luminescence (OSL) samples were taken from homogeneous sand units, mostly with a median grain size of 210–300 µm, to exclude the effect of lithological variability. A total of 11 OSL samples were collected. Eight samples were collected from fluvial pointbar deposits at 1–2 m depth below the surface, to prevent the influence of later overbank deposition, soil formation and bioturbation. Sample Kraamweg 8 was taken from the sandy infill of a meander channel. Samples Bergerweg 14 and 16 were taken at a shallower depth from aeolian dune deposits overlying meandering level B.
Quartz optical dating results of pointbar, channel and aeolian sediments from level B. For location of the samples see Figure 4.
In the laboratory of the Netherlands Centre for Luminescence Dating (situated at Delft University of Technology at the time of analysis) samples were split into two parts. Outer ends of the sample tubes were exposed to light during sampling, and were used to determine the dose rate. Activity concentrations of radionuclides were determined with gamma spectrometry, and dose rates were calculated taking into account attenuation due to water and grain size, and a contribution from cosmic rays (for details see Sevink et al., Reference Sevink, Koster, Van Geel and Wallinga2013 and references therein).
Material from the inside of the sample tubes was prepared for equivalent dose determination. Quartz grains of size 180–212 µm were prepared for analysis, and purity of quartz extracts was tested by measuring the luminescence response to infrared stimulation (Duller, Reference Duller2003). Luminescence measurements were made using Risø TL/OSL-DA-15 and DA-20 instruments (Bøtter-Jensen et al., Reference Bøtter-Jensen, Andersen, Duller and Murray2003), equipped with an internal 1.48 GBq 90Sr/90Y source for irradiation, and blue light-emitting diodes (LEDs) for optical stimulation (470 nm, ~35 mW cm−2). Quartz luminescence signals were collected by a selected photomultiplier tube, using 7.5 mm of U340 filters to block the stimulating light and select the quartz emission around 340 nm. The single-aliquot regenerative (SAR) dose procedure (Murray & Wintle, Reference Murray and Wintle2000, Reference Murray and Wintle2003) was applied for equivalent dose determination, using small aliquots (2 mm; ~100 grains) to minimise within-aliquot averaging and allow detection of heterogeneous bleaching (Cunningham et al., Reference Cunningham, Wallinga and Minderhoud2011).
Based on preheat plateau tests for a subset of the samples, a preheat of 240°C and cutheat of 220°C were selected. Stimulation was for 20 s at 125°C; an early background approach was adopted to obtain net OSL signals (Cunningham & Wallinga, Reference Cunningham and Wallinga2010; initial signal 0–0.5 s; background 0.5–1.75 s), and at the end of each SAR cycle the signals were bleached by blue LEDs for 40 s at 260°C. The adopted protocol was tested using a dose recovery test, yielding excellent results (average 0.99 ± 0.02; n = 25). Each measurement included a recycling test (threshold 10% from unity), recuperation test (threshold 10% of first regenerative response) and feldspar test (threshold infrared signal <20% of blue signal, and infrared depletion <10%). Measurements were repeated for each sample until at least 22 aliquots passed rejection criteria; for one sample additional measurements were made (n = 40) to allow interpretation of the complex dose distributions. To obtain burial doses from the equivalent dose distributions, the Central Age Model (Galbraith et al., Reference Galbraith, Roberts, Laslett, Yoshida and Olley1999) was used for all but one sample (see discussion below). Burial ages were obtained by dividing the equivalent dose by the dose rate, taking into account all uncertainties in both. Uncertainties reflect one sigma, and ages are reported in ka relative to the year of sampling (2008). For sample NCL-6108166 an overdispersed and skewed dose distribution indicated heterogeneous bleaching. For this sample, a bootstrapped version of the minimum age model (BSMAM) was used (Cunningham & Wallinga, Reference Cunningham and Wallinga2012). Based on the overdispersion of the other fluvial samples in this study, an overdispersion of 15 ± 2% due to other sources than heterogeneous bleaching was used for the BSMAM analysis.
Results
Three (palaeo)floodplain levels, separated by terrace scarps, were distinguished in the lower reaches of the Roer valley (Fig. 3), based on elevation, differences in surface morphology and lithological characteristics.
Level A
Morphology
The elevation in the study area declines from 34 m a.s.l. near Karken to 28 m a.s.l. in the downstream area near Linne, resulting in a valley gradient of c. 55 cm km−1. Level A is characterised by gentle relief of 1.5 m and an undulating morphology; neither fluvial channel patterns nor pointbar features were observed at this level. Extensive aeolian dune fields with parabolic dune forms, up to 10 m high, are present at level A. These dune fields are situated at the northeastern side of circular to oval-shaped deflation hollows (Fig. 3).
Lithology
Few observations were done on this level. In the downstream area, fine to medium and locally gravelly coarse sands are overlain by 1–2 m loamy sand and sandy loam (Fig. 6: cross section AAʹ, unit A1) (after Van Bostelen & Vreugdenhil, Reference Van Bostelen and Vreugdenhil2009). The upper ~5 m mostly consist of slightly to moderately loamy, fine to medium sand (150–420 µm) with weakly developed fining-upward sequences of <2 m thickness (unit A2). In the more upstream part of the study area medium to coarse gravelly sand is overlain by 2–3 m of loamy, fine to medium sand (Fig. 6: cross section BBʹ, units A1 and A2) (after Vreugdenhil, Reference Vreugdenhil2011).
Chronology
One sample was taken from level A for OSL dating (NCL-6108159; Stella Maris 3) (Fig. 4; Table 1). The age of 15.5 ± 1.1 ka indicates a Late Pleniglacial age close to the Late-glacial boundary of 14.7 ka. The equivalent dose distribution for this sample was broader than for other samples, with an overdispersion of 23%. This may indicate mixing of grains of different age, or heterogeneous bleaching of the sample. However, there were no clear signs of poor-bleaching (e.g. skewed distribution), and hence the Central Age Model (CAM) was adopted for determining the burial dose.
Level B
Morphology
Level B is situated south of the present-day Roer meander belt as a palaeosystem that is incised in level A (Fig. 3). It dips from 32 m a.s.l. at the Dutch–German border near Karken to 25 m a.s.l. near the confluence with the Meuse (valley gradient of c. 62 cm km−1). The channel-belt width increases from c. 1 km in the southeast to 2 km in the northwest. Level B is characterised by well-developed high-sinuosity palaeomeanders of several generations, locally with neck-cut-off morphology, testifying to intense lateral channel migration. The present relief difference between palaeochannels and pointbar surfaces increases from c. 1 m in the more upstream part of the study area to c. 2 m in the downstream part. In the downstream part of level B, palaeomeanders with a larger amplitude and width are present at a lower elevation of c. 25 m a.s.l. (Figs 3a and 4). These larger meanders are dissected by a southeast–northwest-oriented small terrace scarp near Lerop, resulting in a local intermediate level between level B and the Holocene floodplain C (Fig. 3). Locally, northeast and north of Montfort, level B is covered by parabolic dune forms (Figs 3 and 4).
Lithology
The channel lag and pointbar deposits in the downstream part, at cross section AAʹ, are up to 5 m (smaller meanders) or 6 m (larger meanders) thick, with a lower boundary at c. 22 m a.s.l. (Van Bostelen & Vreugdenhil, Reference Van Bostelen and Vreugdenhil2009) (Fig. 6: cross section AAʹ). The lithological successions are characterised by a channel lag deposit of slightly gravelly, medium to coarse sand (300–850 µm) (unit B1), overlain by pointbar sediments with fining-upward sequences of coarse to medium sand (300–600 µm) grading to fine sand (150–300 µm) (unit B2). The upper part consists of loamy sand (105–210 µm) and sandy loam (unit B3). The lithology and surface morphology indicate that the upper part was formed by overbank deposition of suspended material on the pointbars. The channel lag and pointbar deposits in the more upstream area are c. 3 m thick with a lower boundary at 29 m a.s.l. (Vreugdenhil, Reference Vreugdenhil2011) (Fig. 6: cross section BBʹ). The lithological trends are similar, but the sedimentary sequences are thinner and gravel content of the channel lag deposits is higher than in the downstream area.
The residual channel-fill deposits in the downstream area are generally 2–3 m thick in the smaller abandoned meanders and up to 4 m in the large meanders and are dominated by clastic deposits of the Roer (Fig. 6: cross section AAʹ, unit B4) (Van Bostelen & Vreugdenhil, Reference Van Bostelen and Vreugdenhil2009). Most channel fills show a fining-upward sequence with medium to coarse sand (300–600 µm) grading to fine sand (150–300 µm) and loamy fine sand (150–210 µm). Locally the upper part of the fill consists of clay, loam and loamy sand. The residual channel fills in the upstream part of the study area are generally 1.5–2.5 m thick and contain an alternation of organic (gyttja and peat) and clastic units (loam and clay) deposited by the Roer and tributary brooks (Fig. 6: cross section BBʹ, unit B4) (Vreugdenhil, Reference Vreugdenhil2011).
Palynology
Because of the presence of fine-grained and organic units in the abandoned meanders in the upstream part of level B, five detailed cross sections were made over the palaeochannels (Figs 5 and 7). The deepest locations were sampled for pollen analysis. The general outline of the vegetation development is described below. A detailed description of the lithology, local pollen zones and biostratigraphical interpretation is provided in the Appendix. The local pollen zones are correlated with the radiocarbon-dated biostratigraphy of Hoek (Reference Hoek1997a,b).
Sites Karken, Geraerds and Haaserdriesch (Figs 8–10) demonstrate that the first infill of the abandoned channels is characterised by low arboreal pollen and high values in Poaceae and heliophilous herbs that can be correlated with the Older Dryas (PAZ 1c) biozone. The Kitscherholz palaeochannel (Fig. 11) shows a younger infill that started in the later part of the Younger Dryas biozone.
The vegetation shows an increase in trees, especially of Betula during the Allerød (PAZ 2a) (Figs 8–10). The increase of Pinus, characteristic for the second part of the Allerød (PAZ 2b), was only observed at site Karken (Fig. 8). Shielding effects by Salix stands near the Haaserdriesch palaeochannel may have reduced the influx of Pinus pollen (Fig. 10).
The Younger Dryas biozone (PAZ 3) is reflected by an increase in grasses and heliophilous herbs in the diagrams, indicating more open vegetation and unstable soil conditions (Figs 8–10). The lithology of the channel fills at the three sites changes from organic to clay deposits at the Allerød to Younger Dryas transition.
The pollen diagrams show a rapid increase in arboreal pollen marking the onset of the Early Holocene (Figs 8–11). Diagrams Geraerds and Karken reveal a hiatus, associated with soil formation, covering the Preboreal biozone (Figs 8 and 9). In the deeper palaeochannel Haaserdriesch, the Early Holocene record is complete. Juniperus and Betula show a peak during the Preboreal (Figs 10 and 11), followed by an increase in Pinus and Corylus during the Boreal (Figs 8 and 9) and Alnus and Tilia during the Atlanticum (Fig. 8: diagram Karken).
In contrast to the upstream area, the meander-fill sediments in the downstream part of level B are exclusively clastic with a fining-upward sequence of sand, sandy loam and loam. The absence of organic material in the abandoned meanders is caused by groundwater levels below the base of the palaeochannels. The deep groundwater level is related to the incision of the Late-glacial and Early Holocene Meuse and downstream part of the Roer, therefore inhibiting the formation of oxbow lakes and organic accumulation. Instead, brown soils (Dutch ooivaaggronden, vorstvaaggronden) developed in the palaeochannels during the Holocene (Bodemkaart van Nederland, 1972).
Chronology
In order to date the fluvial activity and abandonment of meandering level B, eight OSL samples (NCL-6108158 and -160 to -166), mostly from sandy pointbar deposits, and two samples from overlying dune sediments were dated (Table 1). It was hypothesised during sampling that spatial age differentiation could be present between the different meander generations of level B with younger meander sediments in the direction of the present-day floodplain (Fig. 4). However, the ages range from 16.0 ± 1.1 to 12.7 ± 0.7 ka, and many dates overlap within errors. No systematic age trend can be observed in the morphological development.
Two samples from aeolian sediments from a parabolic dune complex overlying level B (Fig. 4) (NCL-6108157-156; Bergerweg 16 and 14) were dated at 13.6 ± 0.9 and 0.85 ± 0.05 ka / AD 1158 ± 49. The dates provide a minimum age for level B. The date of Bergerweg 16 (0.8 m) is stratigraphically consistent with the date of 14.1 ± 0.8 ka from the underlying fluvial pointbar (NCL-6108158; Bergerweg 16: 2.4 m). The aeolian deposits and parabolic dunes can be assigned to the Late-glacial, and the OSL age agrees within errors with the Younger Dryas stadial (12.9–11.7 ka) (Kasse, Reference Kasse, Troelstra, Van Hinte and Ganssen1995b). The date of Bergerweg 14 (0.85 ± 0.05 ka / AD 1158 ± 49) indicates a Medieval age and is surprisingly recent as it was expected to be from the same lithological unit as Bergerweg 16. Apparently, the sand has been locally reworked by aeolian processes and redeposited as drift sand.
Level C (Holocene floodplain)
Morphology
The present-day Roer floodplain is situated at c. 28 m a.s.l. at the Dutch–German border and 19 m a.s.l. at the confluence with the Meuse and has a floodplain gradient of c. 75 cm km−1 (Fig. 3). This gradient is steeper than the valley gradient of meandering level B, and both gradient lines converge upstream c. 10 km southeast of the Dutch–German border. The Roer is an unconstrained meandering river and the floodplain is characterised by many meander scars, visible in the terrace scarp of c. 4–5 m between levels B and C in the downstream reach, and as remnants of oxbow lakes especially in the younger parts of the floodplain along the present-day river (Fig. 3a). Within the Holocene floodplain, two sublevels can be distinguished, especially downstream of the Dutch–German border (Janssens, Reference Janssens2011). In the downstream part of the floodplain, near the confluence with the Meuse, these two sublevels merge and only one floodplain can be distinguished. The morphological details and fluvial development of the Holocene floodplain level C are outside the scope of this paper, and only the main elements that constrain the Late-glacial to Early Holocene evolution are presented here.
Lithology
The lithology of the Holocene floodplain is characterised mostly by fining-upward sequences (Janssens, Reference Janssens2011, fig. 8). Near Karken the floodplain sequence is 3 m thick, near Herkenbosch 4 m and downstream near Melick the sequence is 5 m thick (after Van Bostelen & Vreugdenhil, Reference Van Bostelen and Vreugdenhil2009; Janssens, Reference Janssens2011) (Fig. 6: cross section AAʹ). At the base, channel lag deposits of gravel and coarse sand are overlain by medium to coarse-grained pointbar sands (300–850 µm) and loamy sand and sandy loam overbank deposits (c. 2 m) (units C1, C2, C3). The channel lag deposits of level C are coarser-grained and have a higher gravel content than those of level B, where gravel is scarce. In addition, loam and clay layers that are restricted at level B to channel fills are more extensive at level C. Channel fills at level C are up to 5 m thick in the downstream part near Melick (Fig. 6: cross section AA, unit C4). They consist of loamy fine sand at the base, grading upward into silty clay, clay loam and sandy loam deposited by the Roer. In general, the sediments at level B are dominated by pointbar sands, while level C is characterised by pointbar sands and floodplain loams.
Interpretation and discussion
The Late-glacial and Holocene fluvial development of the Roer is characterised by channel pattern changes, incision, and by terrace formation and abandonment of a channel belt.
Late Pleniglacial evolution (level A)
The undulating, irregular morphology and absence of fluvial morphology at the surface of level A, in combination with loamy fine sand in the upper 1–2 m of the sedimentary succession, indicates that level A may be covered by aeolian coversand deposits with limited relief (Fig. 3a). At the base of the succession thin fining-upward sequences, gravel and coarse sand have been observed pointing to fluvial channel deposition. The intermediate zone, characterised by an alternation of medium sand and loamy beds, may be interpreted as fluvial bar and overbank sediments (Fig. 6: cross section AAʹ, top unit A1). The absence of fluvial surface morphology hampers the genetic interpretation in terms of meandering or braided, as all fluvial systems are characterised by fining-upward sequences (Bridge, Reference Bridge1985). However, the absence of long fining-upward sequences, as has been found at level B, and the apparent absence of large channel fills, might indicate an origin in sandy braided systems. The sedimentary succession at level A resembles the fluvial to aeolian succession as described in the Meuse valley near Grubbenvorst (Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995, photo 2; Reference Kasse, Vandenberghe, De Corte and Van den haute2007). A partial co-functioning of the braided system and coversand deposition may have occurred. The absence of curved meander scars at the transition to older terraces might also point to deposition by straight or low-sinuosity channels during the formation of level A. Level A rapidly radiates in a downstream direction, which can be interpreted as an aggrading sandy alluvial fan that formed because of the gradient decrease in the confluence area of the Roer and the Meuse, along the eastern margin of the Meuse valley (Fig. 3b).
The single optical date of 15.5 ± 1.1 ka obtained from fluvial sediment (210–300 µm) (Table 1: NCL-6108159, Stella Maris 3) indicates formation of level A during the final Late Pleniglacial (Fig. 4). Level A dips to the northwest in the direction of the Meuse valley and merges with Meuse terraces near Roermond (Van den Broek & Maarleveld, Reference Van den Broek and Maarleveld1963; Van den Berg & Schwan, Reference Van den Berg, Schwan and Van den Berg1996). These Meuse terrace levels, covered by aeolian deposits, have been described previously (Mol et al., Reference Mol, Vandenberghe, Kasse and Stel1993; Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995: level 1; Van den Berg & Schwan, Reference Van den Berg, Schwan and Van den Berg1996: level 2; Van Huissteden & Kasse, Reference Van Huissteden and Kasse2001). At Grubbenvorst, 28 km north of Roermond, Kasse et al. (Reference Kasse, Vandenberghe, De Corte and Van den haute2007) optically dated the fluvial to aeolian succession of terrace level 1, including the Older Coversand I and Older Coversand II, separated by the Beuningen Gravel Bed. For the fluvial deposits, the obtained ages are between 29 and 25 ka, whereas overlying fluvio-aeolian deposits (Older Coversand I) were dated between 25 and 17 ka. West of the Meuse river, near Neer, this terrace was investigated by trenching to reconstruct the fault displacement rate along the Peel Boundary Fault (Van den Berg et al., Reference Van den Berg, Vanneste, Dost, Lokhorst, Van Eijk and Verbeeck2002). The fluvial to aeolian succession of the foot wall, below the Beuningen Gravel Bed, was infrared stimulated luminescence (IRSL) / thermoluminescence (TL)-dated between 36.1 ± 3.9 / 25.6 ± 2.8 ka and 16.4 ± 2.5 / 16.4 ± 1.9 ka (Frechen & Van den Berg, Reference Frechen and Van den Berg2002). Both studies indicate a Middle to Late Pleniglacial age. However, it cannot be excluded that level A is polygenetic. According to Van den Berg (Reference Van den Berg1989), different terraces separated by inconspicuous scarps are present, which were formed during the Saalian (Caberg and Eisden Lanklaar Terraces) and Weichselian (Mechelen aan de Maas Terrace).
The Late Pleniglacial was characterised by cold conditions with mean winter temperatures of −25 to −15°C and mean summer temperatures of 10 to 15°C (Renssen & Isarin, Reference Renssen and Isarin2001). In the Roer catchment and wider surroundings of the Lower Rhine Embayment and the Ardennes foothills the sparse vegetation (Hoek, Reference Hoek1997a) and continuous permafrost conditions (Vandenberghe & Pissart, Reference Vandenberghe and Pissart1993; Kasse, Reference Kasse2002) resulted in a low soil permeability and overland flow during snow melt (Woo & Winter, Reference Woo and Winter1993). These factors caused a high sediment supply from the upstream catchment of the Roer in the Eifel mountains, which downstream, in the Roer Valley Graben, resulted in aggradation by a braided river system forming a low-angle alluvial fan in the Roer–Meuse confluence area. A decline of fluvial activity followed by aeolian deposition is inferred in the Roer valley (Fig. 6: units A2, A3). A similar evolution was also determined for several western and central European valleys (Schwan, Reference Schwan1987; Vandenberghe & Van Huissteden, Reference Vandenberghe and Van Huissteden1988; Bohncke et al., Reference Bohncke, Kasse and Vandenberghe1995; Van Huissteden & Kasse, Reference Van Huissteden and Kasse2001; Vandenberghe et al., Reference Vandenberghe, Derese, Kasse and Van den haute2013; Zieliński et al., Reference Zieliński, Sokołowski, Fedorowicz and Zaleski2014; Kadlec et al., Reference Kadlec, Kocurek, Mohrig, Shinde, Murari, Varma, Stehlík, Beneš and Singhvi2015) and was related to increasing climatic aridity during the Late Pleniglacial (Kasse, Reference Kasse1997, Reference Kasse2002; Kasse et al., Reference Kasse, Vandenberghe, De Corte and Van den haute2007).
Pleniglacial – Late-glacial fluvial transition (level A to B)
The biostratigraphy of the basal infill of several abandoned meanders in the more upstream part of level B (Figs 5 and 7) shows that the infill occurred during the Older Dryas and Allerød biozones (Figs 8–11), and indicates that the Roer river was a meandering river during the Older Dryas or even the Bølling period. The start of the channel fills is in close agreement with palynological results from site Posterholt where an Allerød age was established at the base of a channel fill (Bakels, Reference Bakels2017). Therefore, the channel pattern change from an aggradational, probably braided, system (level A), blanketed by coversands, to an incising meandering system (level B; Fig. 3) possibly occurred during the early Late-glacial or at the Pleniglacial to Late-glacial transition (~14.7 ka).
The Roer fluvial development echoes that of the Meuse system in central and northern Limburg between Roermond and Cuijck, although the start of the meandering phase in the Roer may be earlier than in the Meuse valley. In the Meuse valley most Late-glacial palaeomeanders have an Allerød or Younger Dryas basal infill (Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995; Van den Berg & Schwan, Reference Van den Berg, Schwan and Van den Berg1996; Huisink, Reference Huisink1997) although older fills have also been reported (Tebbens et al., Reference Tebbens, Veldkamp, Westerhoff and Kroonenberg1999), indicating that the Meuse was a meandering system during and before the Allerød period. A transitional phase of pattern change from braided to meandering, as reported for the Meuse and Niers rivers (Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995, Reference Kasse, Hoek, Bohncke, Konert, Weijers, Cassee and Van der Zee2005 respectively), was not observed in the Roer valley. The Meuse and Niers had a gradual or delayed fluvial response to the rapid warming at the start of the Late-glacial at 14.7 ka of c. 700 radiocarbon years, despite the abrupt change in temperature and precipitation (Kasse et al., Reference Kasse, Hoek, Bohncke, Konert, Weijers, Cassee and Van der Zee2005). The delayed response can be explained by reach-to-reach and catchment-to-catchment differences in the balance between discharge and sediment supply (Vandenberghe, Reference Vandenberghe1995), which would be affected by the gradual migration of vegetation and associated stabilisation of the surface by vegetation and soil formation (Hoek, Reference Hoek2001). The pollen diagrams of the Roer valley (Figs 8–11) show that the vegetation was still rather open and characterised by heliophilous herbs, Poaceae and Betula.
Aggradation of the river valley by a braided system (level A) during the Late Pleniglacial was probably related to the snow-melt peak discharges, a frozen subsoil and high sediment supply (as a result of weathering and lack of vegetation). During the early Late-glacial, peak discharges probably decreased because of permafrost degradation and the resulting increased soil permeability. Because of higher winter temperatures and reduced snow cover, discharge became more regular, and as a consequence the braided system (level A) transformed into a meandering system (level B) because of reduced stream power (Van den Berg, Reference Van den Berg1995; Bogaart et al., Reference Bogaart, Van Balen, Kasse and Vandenberghe2003; Fig. 3) despite the floodplain slope increase. Lower floodplain gradients would also imply a reduction of stream power and channel-pattern changes from braided to meandering. However, the Roer valley gradient of meandering level B (c. 62 cm km−1) is higher than for level A (c. 55 cm km−1). Therefore, gradient changes cannot be responsible for the observed river system changes from braided to meandering. The gradual immigration of vegetation, start of soil formation and increasing cohesion of the surface resulted in a gradual decrease in sediment supply leading to incision by the Roer river, which is evident from the small terrace scarp (c. 2 m) between levels A and B (Fig. 3).
Similar fluvial system changes have been described in central and northwestern Europe during this time interval (Kozarski, Reference Kozarski1983; Kiden, Reference Kiden, Starkel, Gregory and Thornes1991; Bohncke et al., Reference Bohncke, Kasse and Vandenberghe1995; Rose, Reference Rose and Frenzel1995; Schirmer, Reference Schirmer1995; Pastre et al., Reference Pastre, Limondin-Lozouet, Leroyer, Ponel and Fontugne2003; Starkel et al., Reference Starkel, Gębica and Superson2007, Reference Starkel, Michczyńska, Gębica, Kiss, Panin and Perşoiu2015; Erkens et al., Reference Erkens, Hoffmann, Gerlach and Klostermann2011; Meylemans et al., Reference Meylemans, Bogemans, Stormea, Perdaen, Verdurmen and Deforce2013; Turner et al., Reference Turner, Tolksdorf, Viehberg, Schwalb, Kaiser, Bittmann, von Bramann, Pott, Staesche, Breest and Veil2013). It has been postulated that the timing of the fluvial system transformation from braided to meandering is not uniform for all river systems but that it depends on the discharge, grain size, gradient and surrounding vegetation, which all influence sediment supply and river transport capacity. Turner et al. (Reference Turner, Tolksdorf, Viehberg, Schwalb, Kaiser, Bittmann, von Bramann, Pott, Staesche, Breest and Veil2013) reported that changes in the Jeetzel river from a braided to an incised system even pre-dated the Late-glacial warming. They also stress the importance of the vegetation cover determining evapotranspiration, soil cohesion and sediment supply. According to Kasse et al. (Reference Kasse, Bohncke, Vandenberghe and Gábris2010), the earlier transition from braided to meandering in the Tisza river in Hungary, compared to northwestern and central European rivers, may be related to earlier vegetation development in southeastern Europe.
Late-glacial meandering system (level B)
The OSL dates from pointbar deposits and channel fills of level B range between 16.0 ± 1.1 ka (NCL-6108160) and 12.7 ± 0.7 ka (NCL-6108162), indicating that fluvial activity and lateral migration by a meandering system probably occurred during the Late-glacial period (14.7–11.7 ka) (Fig. 4; Table 1). Due to the range of ages and the precision of the individual estimates, it is not possible to assign specific parts of level B to the Bølling–Allerød interstadial or the Younger Dryas stadial without additional independent dating. The dates are in agreement with previous results obtained from channel fills on meandering level B (Hoek, Reference Hoek1997a, site RMO 3; Bakels, Reference Bakels2017). The sediments have been attributed to the Younger Dryas (biozone 3) supported by a radiocarbon date of 10,440 ± 160 BP (12,761–11,860 cal BP; Oxcal 4.2) (Hoek, Reference Hoek1997b, p. 53). Meandering level B is locally overlain by parabolic dunes OSL-dated at 13.6 ± 0.9 ka (NCL-6108157; Bergerweg 16) and 0.85 ± 0.05 ka (NCL-6108156; Bergerweg 14) (see discussion below) (Fig. 4). Parabolic dune formation in the Roer valley probably occurred during the Younger Dryas stadial (12.9–11.7 ka), similar to the Meuse valley (Bohncke et al., Reference Bohncke, Vandenberghe and Huijzer1993; Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995). Therefore, meandering level B at location Bergerweg is older than the Younger Dryas, as also indicated by the OSL date of 14.1 ± 0.8 ka (NCL-6108158) from the fluvial sediments underlying the parabolic dune. However, this does not necessarily mean that level B was inactive during the Younger Dryas since younger meandering channels are present outside the dune complex (Figs 3 and 4). The optical dates of the meander pointbars and palynological results from the youngest meander fill Kitscherholz indicate that the Roer was a meandering system during the entire Late-glacial period, starting in the Bølling and continuing up to the Younger Dryas – Holocene transition.
Pollen diagrams Karken, Haaserdriesch and Geraerds (Figs 8–10) show an increase in trees, especially Betula, attributed to the Allerød biozone. An increase of Pinus, characteristic for the second part of the Allerød (PAZ 2b), was only observed at Karken (Fig. 8). The rare presence of pine can be explained by the location of the sites in the wide and wet valley (Fig. 5). In addition, shielding effects by Salix stands near the abandoned channels may have reduced the influx of Pinus pollen (Fig. 10). An increase in forest cover and associated increase in evapotranspiration could have resulted in lower bank-full discharges and reduced channel dimensions (Ward et al., Reference Ward, Renssen, Aerts, Van Balen and Vandenberghe2008; Kasse et al., Reference Kasse, Bohncke, Vandenberghe and Gábris2010). In the more upstream part of level B, smaller meanders in the axis of the valley seem to indicate a decrease in channel dimension (Fig. 5). However, in the downstream part this change was not observed (Fig. 4). The increase in vegetation cover will have reduced the sediment supply to the river, resulting in erosion of the river channel. Successive meander generations of level B are situated in a slightly lower position in the landscape, illustrating the tendency for gradual incision during meander formation (Fig. 4). The incision of the Meuse and Roer during the Late-glacial can be explained by the climate and vegetation-controlled change in the discharge/sediment ratio, in combination with the long-term tectonic uplift of this southeastern part of the Roer Valley Graben (Van Balen et al., Reference Van Balen, Houtgast, Van der Wateren, Vandenberghe and Bogaart2000) comparable to other river systems (e.g. Maddy et al., Reference Maddy, Bridgland and Westaway2001). The Late-glacial incision and terrace formation cannot be explained by glacio-isostatic movements as this region was probably subsiding due to forebulge collapse (Busschers et al., Reference Busschers, Kasse, Van Balen, Vandenberghe, Cohen, Weerts, Wallinga, Johns, Cleveringa and Bunnik2007).
Level B is characterised by several generations of palaeomeanders that show cross-cutting relationships and a tendency towards incision (Fig. 4). The optical dates indicate rapid lateral channel migration, pointbar and channel cut-off formation during the Late-glacial, despite the Late-glacial vegetation increase and increased soil cohesion. The presence of the Usselo Soil of Allerød age in northwestern Europe indeed indicates landscape stability outside floodplains (Kasse, Reference Kasse2002; Kaiser et al., Reference Kaiser, Hilgers, Schlaak, Jankowski, Kühn, Bussemer and Przegiętka2009). The strong fluvial dynamics of the meandering system can be explained by (1) the sandy substrate and (2) increase of floodplain gradient due to downstream incision. The cut banks and the subsoil of the meanders consist mostly of fine to medium sand of level A (Fig. 6: cross section AAʹ). The non-cohesive character could have enhanced lateral migration of the meandering channels (Allen, Reference Allen1965; Berendsen & Stouthamer, Reference Berendsen and Stouthamer2001). The valley gradients of the successive levels A, B and C show an increase from c. 55 cm km−1 to 75 cm km−1 (Van Bostelen & Vreugdenhil, Reference Van Bostelen and Vreugdenhil2009). The upstream convergence of the valley-gradient lines indicates a downstream control governed by incision by the Meuse during formation of level B during the Late-glacial (Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995; Van den Berg & Schwan, Reference Van den Berg, Schwan and Van den Berg1996; Huisink, Reference Huisink1997). This resulted in a local base-level lowering at the confluence of the Roer and the Meuse, leading to Roer incision that propagated upstream from the river mouth in the downstream part of the Roer valley. In addition, lateral migration of the Meuse to the east, in the Meuse–Roer confluence area during this period, might have shortened the downstream part of the Roer, thereby increasing the gradient of the Roer. Following the reduction in stream power due to the decrease in peak discharge from the Pleniglacial (level A) to the Late-glacial (level B), the gradual steepening of the gradient by incision generated an increase in stream power that may have caused intensified lateral migration and meandering on level B and in the Holocene floodplain (Holbrook & Schumm, Reference Holbrook and Schumm1999). However, the increase in stream power apparently was insufficient to induce a transformation to a braided system. The headward incision since the start of the Late-glacial has resulted in upstream knickpoint migration of the present-day river up to c. 10 km southeast of the Dutch–German border. East of this point, clear terraces cannot be distinguished in the Roer valley.
Younger Dryas impact on the Roer
Although some delay may be present between the moment of channel abandonment and the first channel infill, the palynological results of the channel fills and the OSL dates of level B indicate that a meandering system was active in the southern valley during the Late-glacial. As the basal part of the youngest channel fill (Kitscherholz) has been interpreted as the Younger Dryas biozone (Fig. 11), it is concluded that the meandering channel belt of level B was still active during the Younger Dryas. The southern course was probably abandoned during the late Younger Dryas or at the Younger Dryas to Holocene transition (see Holocene evolution of level C). In the pollen diagrams Haaserdriesch, Geraerds and Karken (Figs 8–10) the increase in grasses and heliophilous herbs indicates more open vegetation and unstable soil conditions during the Younger Dryas. Betula is the dominant tree species while Pinus is low. Empetrum (diagram Haaserdriesch) and Ephedra (dia-gram Karken) are characteristic elements of the Younger Dryas vegetation.
The Younger Dryas was characterised by low winter temperatures, discontinuous permafrost conditions and opening of the vegetation cover (Kasse, Reference Kasse1995a; Hoek, Reference Hoek1997a; Renssen, Reference Renssen2001). Lower evapotranspiration, more snow accumulation in winter and lower soil permeability due to permafrost or deep seasonal frost will have resulted in a more nival discharge regime with higher peak flood discharges and stream power than during the preceding Allerød period (Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995; Huisink, Reference Huisink1997; Pastre et al., Reference Pastre, Limondin-Lozouet, Leroyer, Ponel and Fontugne2003; Dzieduszyńska et al., Reference Dzieduszyńska, Kittel, Petera-Zganiacz, Brooks, Korzeń, Krąpiec, Pawłowski, Płaza, Płóciennik, Stachowicz-Rybka and Twardy2014). The start of meander fill Kitscherholz seems to indicate that the Roer was a meandering system during the Younger Dryas period. No clear channel-pattern change from meandering to braided has been found in the Late-glacial southern course of the Roer valley, in contrast to the larger rivers Meuse and Rhine (Berendsen et al., Reference Berendsen, Hoek, Schorn and Frenzel1995; Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995; Huisink, Reference Huisink1997; Erkens et al., Reference Erkens, Hoffmann, Gerlach and Klostermann2011; Janssens et al., Reference Janssens, Kasse, Bohncke, Greaves, Cohen, Wallinga and Hoek2012). Intensive lateral meander migration in the northern course (level C) during the Holocene may have erased possible traces of the Younger Dryas course in that area. Apparently the threshold value in stream power change was not crossed. In this respect the Roer resembles the Niers in Germany (Kasse et al., Reference Kasse, Hoek, Bohncke, Konert, Weijers, Cassee and Van der Zee2005), the Scheldt in Belgium (Kiden, Reference Kiden, Starkel, Gregory and Thornes1991; Meylemans et al., Reference Meylemans, Bogemans, Stormea, Perdaen, Verdurmen and Deforce2013) and central European rivers like the Tisza in Hungary (Kasse et al., Reference Kasse, Bohncke, Vandenberghe and Gábris2010) and Warta in Poland (Vandenberghe et al., Reference Vandenberghe, Kasse, Bohncke and Kozarski1994).
Although a valley-wide channel pattern change was not observed, some fluvial changes in the Roer valley might be related to the Younger Dryas cooling. In the downstream part of level B larger and deeper palaeomeanders are present near Stekveld–Lerop (Figs 3 and 4) that indicate an increase in bank-full discharge, possibly during the Younger Dryas, which is comparable to the Niers valley (Kasse et al., Reference Kasse, Hoek, Bohncke, Konert, Weijers, Cassee and Van der Zee2005). The optical dates (15.8 ± 0.9 and 13.5 ± 0.8 ka) of the related meander pointbars, however, do not indicate a Younger Dryas age (Fig. 4; Table 1). The southeast–northwest-oriented straight lineament near Lerop, at an intermediate level between meandering levels B and C (Fig. 3), might represent a small terrace scarp formed by a straight, possibly braided, channel that could have been formed during the Younger Dryas. This interpretation seems to be supported by the optical date of 13.1 ± 0.7 ka (Lerop 9) indicating an Allerød or Younger Dryas age. However, the upstream continuation of this intermediate level could not be established. Both larger meanders and straight terrace fragment have lower elevations than the rest of level B (Figs 3 and 4). This indicates incision during their formation, possibly during the Younger Dryas, which is in accordance with the Younger Dryas fluvial incision in the Meuse valley (Kasse et al., Reference Kasse, Vandenberghe, Bohncke and Frenzel1995; Huisink, Reference Huisink1997).
The palaeochannel fills (Haaserdriesch, Geraerds and Karken) in the upstream part of level B show a lithological change from organic to clay deposition at the Allerød to Younger Dryas transition (Fig. 7). This probably indicates higher peak discharges in the Younger Dryas active channel and deposition of suspended material in the abandoned channels. A similar fluvial response with increased clay deposition during the Younger Dryas in abandoned channels of Allerød age and on terraces was described for the Meuse valley near Beugen and Bosscherheide (Bohncke et al., Reference Bohncke, Vandenberghe and Huijzer1993; Kasse, Reference Kasse1995a) and in the Paris Basin in France (Pastre et al., Reference Pastre, Limondin-Lozouet, Leroyer, Ponel and Fontugne2003).
Parabolic dunes near Montfort, up to 10 m high, and associated deflation hollows (Fig. 3) have been found superposed on the undulating morphology of level A, and they locally overlie meander scars of level B (Van den Berg, Reference Van den Berg1989). The close spatial link of deflation hollows and downwind dune forms indicates that the dunes were formed by local deflation of terrace level A by southwesterly winds (Fig. 3). The possible cause for deflation may be related to the incision of the Meuse and the Roer during the Allerød and Younger Dryas in combination with a reduced vegetation cover during the Younger Dryas. The incision will have resulted in groundwater level lowering in the adjacent terrace level A, making it vulnerable for deflation. In addition, the sandy character of level A (Fig. 6: cross section AAʹ) provided good source material for deflation.
The parabolic dune complex at Bergerweg (sample 16) has been dated at 13.6 ± 0.9 ka, indicating that dune formation occurred during the Late-glacial and possibly during the Younger Dryas stadial (12.9–11.7 ka) (Fig. 4; Table 1). Bergerweg 14 (0.85 ± 0.05 ka / 1158 ± 49 AD), taken from the same dune complex, can be explained by local remobilisation of the dune and drift-sand formation during the Middle Ages (Castel et al., Reference Castel, Koster and Slotboom1989; Vandenberghe et al., Reference Vandenberghe, Derese, Kasse and Van den haute2013). Younger Dryas dune formation is in accordance with previous findings along the east bank of the Meuse between Roermond and Nijmegen (Bohncke et al., Reference Bohncke, Vandenberghe and Huijzer1993; Kasse, Reference Kasse, Troelstra, Van Hinte and Ganssen1995b; Michon & Van Balen, Reference Michon and Van Balen2005) and along the Schelde river between Antwerp and Bergen op Zoom (Vandenberghe et al., Reference Vandenberghe, Kasse, Hossain, De Corte, Van den haute, Fuchs and Murray2004). However, the formation of parabolic dunes in the Roer valley by local deflation of higher terraces contrasts with the genesis of the extensive Younger Dryas dune fields along the Meuse and Scheldt, which are attributed to deflation of sand from the then active Younger Dryas river braid plain and deposition in source-bordering (river) dunes on higher terraces east of the rivers (Bohncke et al., Reference Bohncke, Vandenberghe and Huijzer1993; Kasse, Reference Kasse, Troelstra, Van Hinte and Ganssen1995b).
Holocene evolution (level C)
The youngest OSL date (12.7 ± 0.7 ka) from pointbar sediment and thus meander activity on level B (Fig. 4) in combination with the start of the channel fill (i.e. meander abandonment) at Kitscherholz during the Younger Dryas and Preboreal (Figs 5 and 11), indicates that level B (southern branch) was abandoned approximately at the Younger Dryas to Holocene transition.
Levels B and C are separated by a clear terrace scarp of c. 4–5 m formed by incision of the Roer river (Figs 3 and 4). The downstream incision of the Roer was partly controlled by the incision of the Meuse as indicated by the steepening of the floodplain gradients at levels A, B and C. In addition, incision can be explained by a change in the discharge/sediment supply ratio related to the climatic amelioration and rapid vegetation response at the start of the Holocene (Hoek, Reference Hoek2001; Renssen & Isarin, Reference Renssen and Isarin2001). The pollen diagrams from the abandoned channels at level B show a rapid increase in arboreal pollen during the Early Holocene. Juniperus and Betula show a peak during the Preboreal (Figs 10 and 11: diagrams Haaserdriesch and Kitscherholz), followed by an increase in Pinus and Corylus during the Boreal (Figs 8 and 9: diagrams Geraerds and Karken) and Alnus and Tilia during the Atlanticum (Fig. 8: diagram Karken). The incision event and related lowering of the groundwater table is registered in the channel fills Geraerds and Karken by soil formation and a higher amount of indeterminable pollen in the top of the clay. In the pollen diagrams a hiatus is present since the Preboreal biozone is missing (Figs 8 and 9). In the deeper palaeochannel Haaserdriesch, the regional groundwater table lowering resulted in a drop in the local water level that is reflected by a lithological change from clay and gyttja to moss peat (Fig. 10).
The incision and level B abandonment by the Roer can further be bracketed by comparing biostratigraphical and OSL-dating results from level B (see above) and the Holocene floodplain level C. According to Janssens (Reference Janssens2011, pollen diagrams Turfkoelen and Bennebroek), the oldest fills of abandoned meanders from level C date from the Preboreal and Boreal. These age estimates are in accordance with previous biostratigraphical results from Holocene floodplain sediments south of Melick where the Boreal (biozone 5) has been found (Hoek, Reference Hoek1997a, sites Melick and RMO 4). Therefore, floodplain incision, which had started during the formation of level B during the Allerød and Younger Dryas, intensified during the late Younger Dryas and early Preboreal. Since the Early Holocene sites occur along both valley borders, this indicates that both vertical incision and widespread lateral erosion occurred during this instability phase. The Early Holocene incision of the Roer is in agreement with previous results for northwest and central European rivers (Vandenberghe et al., Reference Vandenberghe, Kasse, Bohncke and Kozarski1994; Huisink, Reference Huisink1997; Starkel, Reference Starkel2002; Kasse et al., Reference Kasse, Hoek, Bohncke, Konert, Weijers, Cassee and Van der Zee2005, Reference Kasse, Bohncke, Vandenberghe and Gábris2010; Erkens et al., Reference Erkens, Hoffmann, Gerlach and Klostermann2011).
The abandonment of level B (southern branch) and incision of the Roer in a northern course (level C) indicates that probably two courses were simultaneously present during the Late-glacial (Fig. 3). It is unlikely that the Roer shifted from level B to level C by channel avulsion, due to the incised position of level B. However, little hard evidence is present for the Late-glacial northern course, except for a Late-glacial meander remnant west of Herkenbosch (Fig. 3). We suggest that most of the Late-glacial fluvial morphology in the northern branch was eroded by the lateral migration of the Holocene river.
Trends in number of active channels
The fluvial development over the studied time interval indicates a reduction of the number of channels and channel belts over the Pleniglacial with an aggrading braided system, to the Late-glacial with two incised meandering channel belts and finally to the Holocene with one deeply incised meandering channel belt. Similar channel-belt reductions have been described previously for the Niers, Oude IJssel and Rhine systems (Cohen, Reference Cohen2003; Kasse et al., Reference Kasse, Hoek, Bohncke, Konert, Weijers, Cassee and Van der Zee2005; Busschers et al., Reference Busschers, Kasse, Van Balen, Vandenberghe, Cohen, Weerts, Wallinga, Johns, Cleveringa and Bunnik2007; Hijma et al., Reference Hijma, Cohen, Hoffmann, Van der Spek and Stouthamer2009; Erkens et al., Reference Erkens, Hoffmann, Gerlach and Klostermann2011; Janssens et al., Reference Janssens, Kasse, Bohncke, Greaves, Cohen, Wallinga and Hoek2012). It is suggested that during an incision phase the channel with the highest discharge erodes more rapidly, making a deeper channel and therefore reducing the flow in the smaller channel belts at the bifurcating point (Cohen, Reference Cohen2003; Erkens et al., Reference Erkens, Dambeck, Volleberg, Bouman, Bos, Cohen, Wallinga and Hoek2009). The supposed mechanism would imply that the northern Roer branch was more important in terms of discharge. In addition, parabolic dune formation, probably during the Younger Dryas, may have obstructed the southern course of the Roer (Van Zuidam, Reference Van Zuidam1980; Van den Berg, Reference Van den Berg1989; Fig. 3). Differential tectonic movements along the faults of the Roer Valley Graben could also explain the northward channel-belt shift (Fig. 1), as along one of the Rhine branches (Cohen, Reference Cohen2003). Along the Peel Boundary Fault (bounding the study area in the northeast) a 1 m displacement since the Younger Dryas has been established (Michon & Van Balen, Reference Michon and Van Balen2005), while no displacement has been observed along the southern Beegden Fault. This should have resulted in a northeastward tilt of the Roer Valley Graben towards the Peel Boundary Fault, possibly leading to a northward shift of the channel to the area of maximum subsidence. This proposed mechanism of northward block tilting and channel-incision competition was previously also postulated to explain parts of the Late-glacial and Early Holocene fluvial evolution of the Rhine–Meuse area (Cohen, Reference Cohen2003; Hijma et al., Reference Hijma, Cohen, Hoffmann, Van der Spek and Stouthamer2009). Finally, the northward shift may be related to collapse of the glacial forebulge during this period (Busschers et al., Reference Busschers, Kasse, Van Balen, Vandenberghe, Cohen, Weerts, Wallinga, Johns, Cleveringa and Bunnik2007). As the centre of the forebulge was situated towards the north, its collapse will have had a similar effect to the tilt of the Roer Valley Graben.
Lithological trends
Levels B and C are both characterised by fining-upward sequences which are comparable in thickness, indicating that the sequences were formed by laterally migrating meandering channels of similar dimensions and depth (Fig. 6: cross section AAʹ). However, some differences in the sedimentary successions have been observed. Level C has a coarse-grained, gravelly channel lag deposit and an overlying sequence dominated by sand and loam, while level B has a channel lag deposit with slightly gravelly, coarse sand and pointbar sediments dominated by sand (unit B2). The differences in grain size cannot be explained by differences in gradient and transport capacity, as the finer-grained Holocene floodplain has a higher gradient than level B. The different character of the channel lag deposits may be related to the grain size of the substrate. The meanders of level B eroded the sandy units A2/A3 of level A, while the meanders of the more incised level C eroded into the gravelly unit A1 of level A (Fig. 6: cross section AAʹ). The different lithological successions at levels B and C can further be explained by the formational time span of the levels and human activity. According to the optical dates and biostratigraphical results, level B was formed in a geologically short period of 3000 years covering the Late-glacial (Fig. 4). Limited time was available for fine-grained overbank deposition on the pointbar systems and in the abandoned meander channels. The wealth of palaeomeanders and meander cut-offs on level B shows that lateral migration rates of the meandering channels were high, eroding previously deposited overbank fines, and therefore preservation of overbank fines was limited (Hobo et al., Reference Hobo, Makaske, Wallinga and Middelkoop2014). In the abandoned channels of level B, fine-grained channel fills are thin and for that reason the palaeochannels are easily visible as geomorphological features in the present-day landscape (Fig. 4). Level C, on the other hand, covers the longer time span of the Holocene, and therefore deposition of suspended material on the floodplain has formed thick fine-grained floodplain and channel-fill units (Fig. 6: cross section AAʹ, units C3/C4). In addition, human activity, deforestation and agriculture may have increased erosion in the catchment of the Roer, leading to a higher suspended load and increased rates of loam deposition in the floodplain (Janssens, Reference Janssens2011). Neolithic occupation was reported on the Aldenhoven loess plateau in the Roer catchment (Kalis & Bunnik, Reference Kalis, Bunnik and Schirmer1990). However, the impact on river systems was still low during the Neolithic, probably because of localised deforestation and limited connectivity of arable fields with the river floodplain as established for the Dijle river by Broothaerts et al. (Reference Broothaerts, Verstraeten, Kasse, Bohncke, Notebaert and Vandenberghe2014a). From the Bronze Age onwards, human impact on the Roer landscape increased (Bunnik, Reference Bunnik1999) and increased loam deposition is reported in river valleys in central Belgium (Broothaerts et al., Reference Broothaerts, Verstraeten, Kasse, Bohncke, Notebaert and Vandenberghe2014a). Especially during the Roman period and the Middle Ages the impact of man on the landscape was strong, which was reflected by high sedimentation rates in alluvial fans and floodplains in the southern Netherlands and central Belgium (De Moor et al., Reference De Moor, Kasse, Van Balen, Vandenberghe and Wallinga2008; Broothaerts et al., Reference Broothaerts, Verstraeten, Kasse, Bohncke, Notebaert and Vandenberghe2014b)
Conclusions
The last glacial to interglacial transition is characterised by rapid climate change and related environmental impact on the landscape and vegetation. In this paper, fluvial channel pattern changes, incision trends and vegetation development in the Roer valley in the southeastern Netherlands and western Germany have been presented for the Late Pleniglacial, Late-glacial and Early Holocene periods. We conclude that:
-
• two phases of terrace formation have occurred coinciding with the Pleniglacial to Late-glacial and Younger Dryas to Holocene transitions. The development is explained by changes in the discharge/sediment ratio, related to climate and vegetation change, and by the downstream control of the incising Meuse river.
-
• the number of channels reduced stepwise from a multi-channel braided system in the Pleniglacial, to a twofold meander-belt system in the Late-glacial and a single-channel meandering system in the Holocene.
-
• the Roer changed from an aggrading braided to an incising meandering system at the Pleniglacial to Late-glacial transition. During the Late-glacial a dynamic and gradually incising meandering system established in the Roer valley as indicated by OSL dating. The sandy substratum and higher river gradient in the downstream part of the Roer, related to the incision of the Meuse, promoted lateral meander migration.
-
• a fluvial system response related to the Younger Dryas cooling could not be established conclusively. The youngest Late-glacial meanders are larger, and a small, straight terrace fragment has been preserved locally. Holocene lateral meander migration and incision might have erased the Younger Dryas fluvial landforms. Incision, groundwater level lowering and opening of the vegetation cover resulted in localised deflation of the sandy Pleniglacial terrace and parabolic dune formation by southwesterly winds.
-
• the Younger Dryas to Early Holocene transition is reflected by enhanced incision, causing terrace formation and preservation of the glacial to interglacial fluvial record. Incision was induced by forest recovery, reduced sediment supply and the lowering base level of the Meuse.
-
• the forcing factors of fluvial system change in the Roer valley are climate (discharge and vegetation change) and downstream base-level control by the Meuse. The sandy subsurface lithology might have played a role in the rapid meander migration.
Acknowledgements
The authors thank Toine van Bostelen and Marlies Janssens for collecting the field data used in this study. Kim Cohen and Henk Weerts are thanked for constructive reviews that improved the paper.
Appendix: Description of the cross sections and pollen diagrams
Cross section Karken
This cross section is situated over an abandoned meander along the right bank of level B (Fig. 5). The channel fill consists of gravelly sand overlain by channel and pointbar sands (Fig. 7). Its morphology shows in-channel bars and two scour pools (cores 1, 2, 5). The fine-grained fill of the pools starts with loamy sand and clay followed by gyttja and oxidised peat. Overlying the peat, clay has been found with traces of soil formation at the top, testifying to a period of non-deposition. Peaty clays were deposited later over the width of the palaeochannel. Coring 5 has been sampled for pollen analysis (Fig. 8).
Zone 1 (145–123 cm; clay) shows low arboreal pollen and high values in Poaceae and heliophilous herbs, in particular Artemisia and Asteraceae. Combined with the presence of Empetrum and Saxifraga, this indicates a relatively cold and dry period. The presence of Corylus, Alnus, Quercus and Tertiary pollen in the clastic basal fill of the channel can be attributed to reworking of older sediments. This zone probably reflects the Older Dryas biozone (PAZ 1c; pollen zones according to Hoek, Reference Hoek1997a).
Zone 2 (123–91 cm; gyttja) is characterised by a rise in Betula up to 60%, wetland herbs and aquatics. Koenigia is present and, as a typical arctic–alpine species, it is associated with the Late-glacial. Therefore, this zone can be attributed to the Betula phase of the Allerød (PAZ 2a), rather than the start of the Holocene. Both vegetation and lithology indicate warmer and wetter conditions and higher lake level in the abandoned meander (Bohncke & Wijmstra, Reference Bohncke and Wijmstra1988; Hoek, Reference Hoek1997a).
Zone 3 (91–82 cm; oxidised peat) shows an increase in Pinus, a decrease in grasses and low values of aquatic species. The lithological change from gyttja to peat indicates lower water levels and a change from open water to a riparian environment. The increase in Pine is characteristic for the Pine-phase of the Allerød and might be a response to increased continentality (Hoek, Reference Hoek1997a; PAZ 2b). The presence of Typha latifolia points to summer temperatures above 13°C (Isarin & Bohncke, Reference Isarin and Bohncke1999).
Zone 4 (82–77 cm; clay) shows a decrease in arboreal pollen, and an increase in Poaceae and heliophilous herbs, like Potentilla, Brassicaceae and Asteraceae, indicating opening of the vegetation probably related to climate deterioration. Steppe elements, like Ephedra, are typical for the Younger Dryas (PAZ 3). The presence of Carya and Quercus and the clastic character of the sediment suggest reworking of sediment in the upstream catchment.
Zone 5 (77–70 cm; clay with soil formation) shows a decrease in Poaceae and upland herbs, a strong increase in Pinus and Corylus and low values of Quercus. The presence of Corylus is characteristic for the Boreal (PAZ 5), indicating that the Preboreal is missing in this diagram. This hiatus may reflect a water-level drop with associated soil formation in the Karken palaeochannel, as indicated by high values of indeterminable and corroded pollen. Groundwater level lowering was possibly related to incision of the active channel at the Younger Dryas to Holocene transition (see discussion). The upper part of zone 5 (70–62 cm; peat) is characterised by high values of Alnus, Corylus and Tilia, which indicates the Atlantic period of the Holocene. Low values of Carpinus may even suggest a Subatlantic age for the upper part of the diagram, although anthropogenic indicators have not been observed.
Cross section Geraerds
Cross section Geraerds has been made over a meander scar of level B, 1.5 km southeast of Posterholt (Fig. 5). A channel lag deposit at c. 29 m a.s.l. is overlain by medium-grained, channel bar and pointbar sands (Fig. 7). The channel was first filled with loam and clay after channel abandonment, gradually grading into gyttja representing shallow lake conditions (core 77, 78). The gyttja is sharply overlain by a clay layer representing increased flooding of the palaeomeander by a nearby active channel. As in Karken, the top of the clay has been influenced by soil formation and a hiatus can be assumed at that level. The upper part of the fill consists of peat and peaty clay. Coring 78 has been sampled for pollen analysis (Fig. 9).
Zone 1 (135–127 cm; clay) is characterised by relatively high non-arboreal pollen percentages and low Betula values. The presence of Juniperus, Empetrum, Plantago, Thalictrum and Asteraceae indicates cold and dry conditions and rather open vegetation that can be attributed to the Older Dryas (PAZ 1c). The presence of Alnus and Quercus can be explained by deposition of clastic sediment in the channel, derived from erosion of older deposits in the catchment.
Zone 2 (127–105 cm; humic clay and gyttja) shows an increase of Betula up to 60% and a decrease of grasses and heliophilous herbs like Helianthemum, Plantago and Empetrum. Wetland herbs and aquatics increase, indicating warmer and wetter conditions or higher lake levels that can be attributed to the Allerød (PAZ 2a). At the top of this zone Pinus increases slightly and aquatics decrease, which might be attributed to the Pinus-phase of the Allerød (PAZ 2b). The low pine values, compared to diagram Karken, can be related to site-specific conditions as site Geraerds is located more in the central wetter part of the palaeovalley. Pinus probably colonised the higher, freely drained surface of terrace level A (Hoek, Reference Hoek2000).
Zone 3 (105–86 cm; clay) is characterised by a decrease in arboreal pollen, probably due to climate cooling. Poaceae show an increase and heliophilous herbs and Juniperus are present, indicating rather open vegetation that can be attributed to the Younger Dryas (PAZ 3). Wetland herbs and aquatics show a maximum in this interval, indicating wetter conditions that may be related to channel reactivation, flooding and deposition of clastic material in this interval.
Zone 4 (86–76 cm; clay and peat) shows a strong increase in arboreal pollen and decrease in Poaceae and heliophilous herbs. Corylus content is high and Quercus is present in low values which is characteristic for the Boreal (PAZ 5). As in the diagram Karken, the Preboreal pollen zone is not observed. Traces of soil formation in the top of the clay and high values of indeterminable pollen indicate a water-level drop in the Geraerds palaeochannel, probably related to incision of the Roer at the Late-glacial to Early Holocene transition.
Cross section Haaserdriesch
Cross section Haaserdriesch is situated in the central part of level B (Fig. 5). The channel lag is found slightly below 29 m a.s.l. and the channel fill is thicker than in the previous cross sections, so a more complete succession might be present (Fig. 7). Gravelly coarse sand of the channel lag deposit is overlain by medium-grained channel and pointbar sands that formed during the active stage of the channel. Sandy loams were deposited at the margins of the channel (cores 98, 95) during the first stage of channel abandonment, thereby decreasing the width of the channel. In the remaining channel (core 97), carbonate-rich gyttja (lake marl) was deposited, abruptly followed by clay and moss peat. Unlike Karken and Geraerds, no signs of soil formation have been observed in the top of the clay. The upper part of the fill consists of clay, peat and loam. Coring 97 has been sampled for pollen analysis (Fig. 10).
Zone 1 (258–245 cm; clayey gyttja) is characterised by low tree values and high percentages of Poaceae (40%) and heliophilous herbs (Asteraceae). The presence of Helianthemum, Plantago and Rumex in the regional vegetation indicates unstable soil conditions and treeless vegetation that can be associated with the cold and dry Older Dryas (PAZ 1c). The high value of aquatics in the local vegetation, mainly Potamogeton, demonstrates the presence of lake conditions in the former channel. The clayey character of the gyttja can be explained by occasional flooding of the palaeochannel shortly after its abandonment.
Zone 2 (245–195 cm; calcareous gyttja) shows a strong increase in Salix up to 40% and a decrease in heliophilous herbs and Poaceae. Betula values are low, probably due to the local presence of Salix. Despite the low Betula values, this interval can be attributed to the Allerød (PAZ 2), also because of the high calcium carbonate content pointing to higher temperature conditions. The drop in Myriophyllum during this zone can be attributed to the general lowering of the groundwater table at the transition from PAZ 2a to 2b (Hoek Reference Hoek1997b). However, the Allerød Pinus-phase (PAZ 2b) is absent in the diagram, possibly due to the nearby presence of dense Salix stands colonising the channel margins and filtering out the Pinus pollen.
Zone 3 (195–158 cm; clay) starts with a brief increase in Poaceae, while Betula decreases. Non-arboreal pollen, like Artemisia, Potentilla, Chenopodiaceae, Asteraceae, Helianthemum and Rumex, increase, indicating unstable soils and open vegetation. The presence of Alnus and Quercus can be explained by reworking, which is supported by the siliciclastic character of the sediment. The presence of Empetrum and Ephedra is characteristic for the Younger Dryas (PAZ 3).
Zone 4 (158–140 cm; gyttja and moss peat) shows a Juniper peak and sudden increase of Filipendula, followed by a strong increase in Betula up to 80%, indicating the start of the Holocene (Preboreal, PAZ 4). The presence of Typha latifolia indicates summer temperatures above 13°C (Isarin & Bohncke, Reference Isarin and Bohncke1999). The transition from gyttja to peat and increase of the wetland herbs (Cyperaceae) can be explained by a decrease in water depth and hydroseral succession. In contrast to Karken and Geraerds, where a hiatus was observed in the Preboreal, at the Haaserdriesch site the Preboreal is present, probably because of the larger channel depth. The upper part of zone 4 (140–125 cm; clay) shows an increase of Pinus, the first occurrence of Corylus and low upland herb values. This zone can be interpreted as the early Boreal (PAZ 5).
Cross sections Kitscherholz and Kitscherholz b
Cross sections Kitscherholz and Kitscherholz b are located in the axis and lowest part of the palaeovalley (level B) and therefore might contain the youngest sedimentary succession (Fig. 5). At the base of the cross section around 29 m a.s.l. gravelly coarse sand and medium sands occur, which were deposited on the channel bed and possibly as pointbars (Fig. 7). Loams and fine sand were deposited at the margins of the channel under low flow velocities, and the wet cross-sectional area of the channel therefore decreased during or shortly after the abandonment (cores 99, 101, 46, 48). Following the final abandonment clayey gyttja, gyttja and oxidised peat were deposited in the remaining small channel under standing water conditions (cores 100, 45, 47). The peat is overlain with a sharp boundary by sandy clay, suggesting channel reactivation by floods with suspended material. Peaty clay, loam and sandy loam in the upper part possibly reflect young brook deposits. Coring 100 has been sampled for pollen analysis, but the sampling interval is large, hampering the biostratigraphic zonation (Fig. 11).
Zone 1 (148–140 cm; clay) has low Betula values and high values of Artemisia, Thalictrum and Asteraceae. The presence of Juniperus might indicate the final stage of the Younger Dryas or the Late-glacial to Holocene transition.
Zone 2 (140–115 cm; gyttja and peat) is characterised by high Betula values (40–60%), Pinus (10–20%) and Poaceae and this zone could reflect the Preboreal (PAZ 4). The changes in the Poaceae curve might reflect the climate oscillation of the Rammelbeek phase.
Zone 3 (115–112 cm; clay) shows the presence of Corylus, indicating the start of the Boreal.
