Terraces of the Ohře River in Žatec Area, Czech Republic: When Models of Holocene Fluvial Development Are Not Sufficient

The development of the Ohře River near Žatec between the Late Glacial Maximum (LGM) and present time was studied to obtain its natural trajectory and thus to understand the role of human impact on floodplain development. The study was based on geomorphic research, sampling and dating fluvial sediments by optically stimulated luminescence (OSL) and anthropogenic contamination by Sn since 16th century. During the LGM and the Pleistocene/Holocene transition (P/H), the river valley was shaped by large, incising meanders creating the entire valley floor. The Holocene river, due to decreased channel forming discharges, only slowly continued the LGM and P/H trajectory by cut-offs of the large palaeomeander necks and channel belt narrowing and incision. The last Holocene incision step was likely triggered by hydrological extremes of the Little Ice Age. The LGM and P/H Ohře River consequently changed from meandering to low sinuous, incised in the valley floor; this incision was completed by river engineering in the 20th century. The Ohře River development can thus be rationalised neither by traditional terrace-step formation model nor by model of the Anthropogenic aggradation.


Introduction
One major aim of the recent floodplain studies has been to distinguish intrinsic (autocyclic and stochastic) and extrinsic (climatic, anthropogenic, and sometimes tectonic) triggers in the floodplain development, typically to isolate climatic and/or anthropogenic signals in fluvial sedimentary records [1]. The meta-analysis of 14 C ages, that is, the identification and dating of periods of elevated fluvial activity, is a prominent step in this direction. It is conducted by constructing probability density functions (PDFs) and searching for the relative cumulation of ages of certain floodplain sedimentary units/facies [2]. PDFs have been employed systematically by many teams to produce a considerable volume of age datasets, but the outputs are not always straightforward. Individual catchments and regions did not produce catchment-independent and supra-regional isochronous patterns attributable to the Holocene climate change. PDFs were challenged from a methodological (theoretical) aspect [3]. Further studies identified the mechanisms by which the expected straightforward climatic record in floodplain sediments is blurred, such as the catchment-or reach-specific impacts of microtopography, bedrock resistance, and fluvial (dis)connectivity [4][5][6][7]. Faust and Wolf [5] explicitly stated that meta-analyses using PDFs "are often accompanied by the release of the age data from their original (stratigraphic) context, and the data is often unreflectedly used and correlated across a variety of scales". Faust and Wolf [5] also noted that researchers commonly tend to limit the number of samples for dating, typically for economic reasons. However, under-sampling can produce biases, and resulting studies can confirm paradigms rather than reveal the real complexity of each studied fluvial system. Understanding river systems requires complete realization of their complexity, and conclusions must be based on careful field observations [1]. This requirement is contrasting with current under-representation of case studies in research publishing.
The need to respect the local settings to isolate climatic or anthropogenic impacts on floodplain patterns has been demonstrated in individual case studies [4,[7][8][9]. For instance, the actual proportion of anthropogenic and climatic impacts on the development of middle-sized rivers in southeast Czech Republic, particularly the Morava River, has not been unequivocally distinguished, despite systematic and lasting research for several decades [9,10]. The fluvial record was more complicated than expected. The impact of autogenic processes has been considered in theoretical models of fluvial development. These focused on the spatial probability of sediment reworking [11], autogenic formation of terraces without exogenic forcing [12], impact of the slope processes on channel migrations [13], and impacts of variable discharge and sediment delivery on near-channel topography [14]. All of these aimed to visualise the consequences of theoretically known phenomena on the real floodplain topography. Autogenic processes and 'stochastic' site-specific features in the studied river segments can substantially blur the climatic and anthropogenic signals occurring in floodplains. However, the fluvial archives are valuable for their research and applicability in reconstructing past environmental changes.
This study deals with the Lower Ohře River valley in the Czech Republic ( Figure 1). This portion of the Ohře River and its entire lower reach does not show considerable recent aggradation [15,16]. However, the 'floodplain siltation' and similar consequences of aggradation have frequently occurred in European rivers in catchments under anthropogenic pressure [17]. The Ohře River has remarkable palaeomeanders in the study area (Figure 1), which have attracted researchers for at least one century [18][19][20][21]. Local authors [19][20][21] assumed that the river has produced a staircase (strath) terrace system with exactly one incision step per glacial cycle [19,20]. It was later recognised that the real number of incision steps per glacial cycle could differ [22,23], and all those steps need not be preserved [13,23]. Numerous studies have documented departures from a simple historical model of the strath-terrace formation, caused by the sediment reworking in cut-and-fill terraces and palaeochannel belts in paraglacial river valleys [24][25][26], polycyclic terraces [27], impact of the bedrock on terrace distribution and parameters [28], and upstream propagation of terrace formations due to changes in the erosion base in the lower reach [29]. These findings should prompt researchers to revisit preceding case studies of fluvial terraces that were not supported by the sediment dating, as shown by several European rivers by [30]. Recent optically stimulated luminescence (OSL) dating showed that the palaeomeanders of the Ohře River downstream of Žatec were formed in the Holocene [16] rather than the Pleistocene, as expected earlier [21] and specified in actual maps by the Czech Geological Survey (CGS). We have continued the critical revisiting of models for the Ohře River development initiated by Elznicová et al. [16] by increase of the OSL dating points, more extensive use of contamination dating for sediments younger than five centuries [31][32][33], and identification of major changes in the fluvial patterns in the study area. The palaeomeanders of the Ohře River offer unexplored opportunities to understand the river behaviour during the Holocene, including periods of expected human inherence. We focused on the P/H transition and Holocene evolution of the Ohře River in the Žatec area. The research questions that our work seeks to answer are: 1) did site-specific phenomena act in the study area, or did the development follow known scenarios? and 2) has the Holocene river development changed in the in the period of anthropogenic impact?

Geological Settings
The Ohře River catchment has varied geology and geochemistry. The uppermost reach, partly in eastern Germany, is covered by metamorphic rocks with a chemical composition close to the mean upper continental crust (UCC; [15]). The middle reach includes several geological and geochemical anomalies. Granite plutons of the Late Palaeozoic age occur in the river catchment near Karlovy Vary and in Slavkovský Les Mountains. Ore regions have occurred in both these areas since the late Mediaeval times [15,32]. Greisen rocks with ore veins are frequent in the Krušné Hory Mountains [32]. Further downstream, the river is incised into the Oligocene and Miocene basaltic rocks of the Doupov Mountains [15]. The Doupov Mountains are a remnant of the Oligocene basalt floods of the České Středohoří Mountains. This was followed by the Miocene subsidence and fluviolacustrine deposition in the Cheb, Sokolov, and most basins of the Ohře Rift [34]. These Miocene sediments represent the bedrock of the Ohře River in the studied area near Žatec [16]. Moreover, the Ohře Rift development has also been responsible for the formation of the Ohře River network [19].

Valley Morphology: Terraces
The Ohře River terraces were first described scientifically about a century ago [18]. The age model for terraces was proposed by Balatka and Červinka [19]. Tyráček [20] attempted to use palaeomagnetic dating of porcelanite for supporting the terrace model; he used clay beds in coal seams that were fired by the self-ignition of the coal in the outcrops. The results were promoted in international literature [21]. The terrace stratigraphy by Balatka and Červinka [19] assumed a regular step-like incision in each subsequent climatic cycle. This reportedly produced as much as 24 terrace levels in a total elevation range of <100 m, that is, ca. 4 m per terrace. Balatka and Červinka [19] as well as Tyráček [20,21] assumed that the terraces represented the entire Quaternary period.
Previous studies [19,20] provide a valuable summary of preceding geological mapping and research, including the palaeopedology and archaeology. It indicated the middle Pleistocene age for the higher terrace levels; in particular, >250 m a.s.l. in the Číňov area and >230 m a.s.l. in the Staňkovice area. This is several tens of metres above the levels addressed in this work. However, the ages of the terraces in the lower levels (<210 m a.s.l.), which have been determined by Elznicová et al. [16], are much less than the previously assumed value.
In our work, the fluvial forms along the Lower Ohre River were studied in detail, downstream of the Nechranice Reservoir in the Žatec area ( Figure 1). Though the valley is 80 m deep, only the youngest fluvial features were studied, which are located within 30-40 m elevation above the valley floor. These lower surfaces were classified as development stages. In some cases, they were considered equal to conventionally defined strath terraces. Their numbering was taken from Elznicová et al. [16]. Their characteristic features and definitions are presented in Table 1 and shown in Figures 2 and 3, respectively.

Anthropogenic Impacts
The entire Ohře River catchment has been intensively used by humans since prehistorical times, from the river valley to the montane regions. The Žatec area has settled since the Neolithic period. The archaeological findings in the Žatec area have been summarised by Holodňák [35]; commonly found are graves, ceramic vessels and shards, and tools from the Neolithic, Bronze Age, and younger periods. These settlements were found in the city of Žatec and practically in all historical villages. In level 0 of the floodplain, traces of pre-historic settlements are absent. In level I, they are exceptional, as that observed in site Horkau (in the historical SK map) or "ovčárna" (sheep farm) situated north of the Nové Sedlo village at 219 m a.s.l., which is 6 m above the current level 0a. Moreover, artefacts from slightly above the surrounding floodplain level show occupation from the Neolithic period to the 19th century. Other long-term settled sites at relatively low elevations above the floodplain are in Libočany (ca. 215 m a.s.l. and ca. 8 m above level 0a, i.e. at level II). Another persistently occupied region was the Staňkovice village at the slope of the palaeomeander II core. Other archaeological findings were identified at higher terraces.
The extent and character of the mediaeval colonisation of the ridges of the Krušné Hory Mountains (Erzgebirge in German), charcoal and glass production, and ore exploitation can be inferred from both the complex reconstruction of a nearby German montane site, Ullersdorf [36], and the extensive contamination of the Ohře River channel sediments [32]. Large-scale deforestation for charcoal produced for glass making and metallurgy changed the local plant cover dominated by broad leaved trees in the Krušné Hory Mountains into secondary forests dominated by spruce.
The locations of the major ore regions in the Ohře River catchment are shown in Figure 1A. The ore veins in the Slavkovský Les and Krušné Hory mountains have likely been exploited from prehistoric times, but the most extensive exploitation occurred in the 16th century for ores bearing Ag, Cu, Pb, and Sn [32,33,37]. The Ohře River system has consequently been contaminated mainly by As, Cu, Pb, Sn, and Zn [32,33,37]. Among these elements, Sn is the least mobile after sediment deposition, and we thus consider it as the most reliable marker of sediments younger than the mining peak [16,31].

GIS Analysis
The geomorphology of the study areas was described using a field survey and a digital terrain model (DTM) derived from the laser scanning (LiDAR) dataset of 2010 (DMR 5G, Czech Office for Surveying, Mapping and Cadastre, ČÚZK). Data processing, analysis, and visualisation were conducted in ArcGIS Desktop 10.7, with the 3D Analyst extension (ESRI, USA). Geomorphological development phases were identified based on the study by Elznicová et al. [16]; these included the terraces, incised (palaeo)meanders, and palaeochannels occasionally preserved on the valley floor. The individual elevations/ages distinguished in the topography were denoted as development phases and numbered from the active channel belt (0) upwards (I to III). More channel generations were identified in stage 0 and indicated by letters, from actual (a) to older (b,c). Segments in development stage I were further classified for active (Im; circular, regularly developing meanders) and passive meandering (Ip; transverse and low-sinuosity channels).
The hydrological/geographical datasets on the channel and active floodplain (digital basis of the water management data, DIBAVOD) were purchased from the T. G. Masaryk Water Research Institute (Prague). Based on this dataset and DTM, the channel elevation and valley width of the Ohře River were calculated using ArcGIS. The main historical maps -the Imperial Obligatory Imprints of the Stable Cadastre (SK maps) produced in the 1840s -were purchased from the ČÚZK.

Geochemical Analyses
The elemental composition of the sediments, which were derived from the youngest floodplain in its development stages, was determined by X-ray fluorescence spectrometry (XRF). We conducted laboratory analyses of the sampled sediments and in-situ analyses of drill cores directly in the field. Sediments were obtained by manually operated soil corers (groove corer with an internal diameter 3 cm, from Eijkelkamp, Giesbeek, the Netherlands), which could be used in sediments ranging from coarse sand to clayey mud with negligible contents content of clasts with diameters of ca. 5 mm or large. In-situ analyses were conducted using either a Niton XL3t 950 GOLDD+ instrument (Thermo Scientific, , Waltham, MA) or a DELTA Premium (Olympus Innov-X, Waltham, MA). Both instruments comprise energy dispersive, Peltier cooled large-area Si drift detectors and a preinstalled software for the quantitative analysis of geological samples. The performance of the in-situ XRF analyses was demonstrated [38] for U and Zn contamination in another Czech floodplain. Tin concentrations below the limit of detection (LOD) of the two instruments in our study were estimated from the error values produced by the quantitation software extrapolated to low concentrations. The lowest values produced by XRF (LOD) were 16 and 38 mg kg −11 for the Niton and DELTA instruments, respectively. Whereas, the corresponding extrapolated error values were 10 and 13 mg kg −11 ; specifically, the LODs of the Niton and DELTA instruments showed approximately two-fold and three-fold errors, respectively. Sn concentrations greater than 50 mg kg -1 were considered to be potentially contaminated.
Selected sediment cores were brought to the laboratory in Institute of Inorganic Chemistry in Řež, and subsequently dried, powdered in a planetary micromill, and subjected to the XRF analysis using a desktop instrument Epsilon 3X (Panalytical, Almelo, the Netherlands) with a Peltier-cooled large-area Si detector (SDD, an energy dispersive detector). The device was calibrated using the certified reference materials, as described previously [33].
The mean concentration of Sn in the UCC is 2.1 mg kg −1 [39], but the local background for the Ohře River sediment is presumably higher due to ore veins and their weathering products in the river catchment of the Krušné Hory Mountains. The thresholds for contamination (Sn concentrations and thickness of the contaminated strata) were estimated using histograms and probability plots by following the generalised approach of Matschullat et al. [40], and applied to the contaminated sediments of the Ohře River [37].

OSL Dating
OSL sampling was conducted from shovel-dug pits that were 30-100 cm deep, with a thick black polyethylene foil cover over the pit and the sampling person. From coarse, gravelly sediments (typical for point bars of large and incising meanders), fine sand was obtained from the matrix among clasts. In finer lithologies deposited by laterally unstable channels, deposits containing sand were obtained from drill holes made by a soil auger with using PVC tubes.
The 150-220 μm grain size fraction was used for the OSL measurements in University of Szeged. The preparation of samples followed standard techniques [41,42]. Although multi-grain aliquots were used during the measurements, sample size was kept low (~20-30 grains/aliquot) in order to be able to capture the scatter of data in relation with partial bleaching and potential post-depositional mixing. Measurements were carried out with a RISØ TL/OSL-DA-20 automated luminescence reader, using blue LED simulation and a Hoya U340 optical filter. Radioactive irradiation was conducted using a calibrated 90 Sr/ 90 Y β-source. Equivalent doses (De) were determined using a single aliquot regeneration (SAR) protocol [43,44].
The OSL intensities were fitted using a single saturating exponential function. Data processing procedures, such as the rejection criteria for selecting acceptable aliquots, identification of feldspar contamination in quartz extracts, and the calculation of aliquot De values are described by Elznicová et al. [16]. Sample De was calculated using either the central age model (CAM), the minimum age model (MAM), or the finite mixture model (FMM) [45]. When selecting between CAM and MAM, decision was made on the basis of the skewness and overdispersion of De values. However, in terms of 3 samples the distribution of single aliquot De values showed two very distinct peaks, in these cases the FM model was applied and two separate sample De results were given.
The environmental dose rate (D*) was measured using a high-resolution, extendedrange gamma-ray spectrometer (Canberra XtRa Coaxial Ge detector). Dry dose rates were calculated from the conversion factors of Liritzis et al. [46]. The attenuation factors for the β dose rates were obtained from [47]. Wet dose rates were assessed based on the in-situ water content [48]. The rate of cosmic radiation was determined based on burial depth, following the equation of Prescott and Hutton [49].

Geomorphic Description
The actual channel slope upstream of Žatec ( Figure 1B) varies between 1.5 and 2.5 m km −1 (table 1C, Figure 2). However, the historical channel slope was lower as the result of the high sinuosity of palaeomeanders. Downstream of Žatec, the active channel slope decreases to ~0.5 m km -1 in average. The slope remains low further downstream [16].
The river valley floor (elevation up to ca. 10 m above the active floodplain 0a, Figure  2 Table 1. Based on the topography (e.g. Figures 1B and 1C) and elevation transects drawn roughly perpendicularly to the valley axis ( Figure 3), characteristic surfaces and form assemblages were identified. These represent fluvial stages labelled from 0 to III according to Elznicová et al. [16], and also presented in Table 1. Table 1. Fluvial forms and surface sediments in development stages or terrace levels. Number of OSL ages in the stages in parentheses, and elevation of the surfaces relative to the lowest surface of the active floodplain 0a. The number of OSL ages includes both this work and work by Elznicová et al. [16]. Channel widths are in form of mean (error), n.p. stands for genuine meanders not present. Footnotes: 1 upstream Žatec, 2 downstream Žatec. The uppermost studied surface (development stage) III is flat and longitudinally descending. It is located 20-30 m above the current channel (0a), running parallel to it. Surface III represents a possible fragmented surface of the strath-terrace level in a conventional sense ( Figure 2A). Fragments of surface III are partly represented by small top parts of circular meander cores of stage II (Číňov, Nové Sedlo, Staňkovice, and Levonice), partly by irregular elongated hills scattered above the floodplain along the edges of the valley (Zlatník, Figures 3B and 4), and by flat remnants of a former floodplain without meander scars (Žatec, transect 10 on Figure 3C, and near Bezděkov).
Stage II surfaces are mostly represented by abandoned incising (ingrown) palaeomeanders. They were formed by simultaneous lateral and vertical erosion, as discussed by Kiss et al. [50] for the Tisza River floodplain. The incising character of meanders in stage II is obvious from the descending meander core surfaces in the elevation transects shown in Figures 3 and 6. Consistent with the incising character of stage II meanders, there are surfaces at elevations between levels II and I ( Figure 2B).
The floodplain surface of fluvial stage I upstream of Žatec is flat and hosts the meander scars remaining after laterally unstable channels, frequently with abandoned loops (Nové Sedlo I, Alte Eger, and Libočany; Figures 1B, 2B, 3, and 4). These channels have a nearly transverse direction relative to the valley axis. The Alte Eger (Old Ohře in German) was the active channel in historical times, as it is represented by the SK map (Figures 4  and 5). Its meander core is only about 1 m above the current floodplain 0a (Figure 4). The Alte Eger meander scar is still not filled by sediments, indicating channel activity until the recent past. The channel slope in this fluvial stage was much lower (0.5-1.5 m km −1 ) than the present (ca. 2 m km −1 ; Figure 2C). Neither of the boundaries between stages II, I and 0 are sharp, but consist of a continuum of individual levels ( Figure 2B).  Tables 2 and 3.
On the active floodplain (0) three channel belts were identified. The most recent channel belt (0a) has been shaped by river channel engineering, particularly by channel straightening and shortening in the 20th century. The engineering is particularly obvious in the channel topography upstream of Libočany ( Figure 5B) and near Tvršice ( Figure 5C). The original (pre-engineered historical) multi-channel system (0b) is clearly discernible in the DTM upstream of Libočany ( Figure 4). The deepest meander scars correspond to the channel position in the SK map from the first half of the 19th century. Dikes were constructed around the artificial 20th century channel ( Figures 5D and 5E), mainly downstream of Žatec, to protect the villages, crop fields, and hop gardens in former palaeomeanders. Lateral channel movement near settlements and towns are prevented by revetment ( Figure 5F). Based on their geomorphological setting and the applied dating methods, the fragments of an older channel belt (0c) were also identified. These palaeomeanders have been preserved at the distal parts of the active floodplain.

OSL Ages
From the study area some OSL ages were already published [16] and also used hereafter ( Table 2). The new OSL ages of the Ohře River sediments are presented in Table 3. Based on the OSL data the timing of the formation of the levels described in Chapter 3.1. was determined. While the fluvial stage 0a represents the modern channel belt, the palaeochannel belts of 0b and 0c that are a few metres above 0a were formed in historical times. Stage I includes Middle to Late Holocene meandering channels and their floodplains; stage II denotes the large meanders formed around the P/H transition; and stage III includes flat surfaces developed during the Late Pleistocene. To understand the development history of the valley, the OSL samples were collected from at least two points of the ingrown meanders ( Figures 1B,1C,6). The uppermost sample was taken from the rim of the meander core: here the point bar sediments represent the beginning of the incision phase. The lowermost sample usually represents the youngest (last) point bar of the palaeochannel, thus the end of the incision of a given form (or stage) could be determined. The timeframe for the development stages summarised in Table 1 was complemented by the information derived from the Sn contamination discussed in sections 3.3 and 3.4.  Table 1.
The plot of elevation against OSL ages at the sampling sites ( Figure 7) shows a general incision tendency of the river valley floor since (at least) the Last Glacial Maximum (LGM; terminated 19-20 ka). The oldest age yet obtained, 36.50±2.89 ka at the top of the Číňov meander core (CIN-1), fits the same trend. Figure 7. The OSL ages plot against elevation above the lowest floodplain level 0a, or 0b. OSL ages were taken from previous ( Table 2) and actual works (Table 3).

Sediment Geochemistry
The geochemistry of the sediments in the valley floor was investigated for 1) Sn contamination by historical mining and 2) markers of the sediment provenance such as K/Rb known from studying the Miocene sediments of the Most Basin [51] and recent Ohře River sediments [15]. Contamination by Sn was found in the lowermost floodplain levels (0a and 0b), and only exceptionally in the 0c meander scars (Table 1).
A histogram of the Sn concentrations in the studied sediments obtained by laboratory XRF is shown in Figure 8A. The major Sn concentration mode centred at ca. 8 mg kg −1 is ascribed to uncontaminated sediments (slightly above the UCC), while values above ca. 20 mg kg −1 are attributed to contamination. The main Sn-contamination mode was centred at 250 mg kg −1 . In majority of the examined Ohře River sediments, Sn concentrations were indirectly proportional to the K/Rb ratio ( Figure 8B). However, several sediment samples showed increased Sn concentrations at a relatively high K/Rb ratio, i.e. the Sn concentrations and K/Rb ratios were probably not directly interrelated, but their change in the Ohře River sediments was temporally close. Three groups of sediments can thus be distinguished ( Figure 8B): uncontaminated (low Sn), Sn-contaminated with K/Rb<120, and Sncontaminated with K/Rb>120. The variability of K/Rb was not associated with the sediment lithology, as no relationship between K/Rb and Al/Si was found; the latter is a versatile proxy for the sediment grain size. The K/Rb histogram shows two main modes: one centred at 100 and another at 140. Both values are substantially lower than the global average (K/Rb in the UCC is 280; [39]). The low K/Rb values in the studied sediments are a consequence of geochemical anomalies in the Ohře River catchment. Low values of K/Rb have also been found in two sediment sources of the Ohře River in the study area: the Miocene Most Basin sediments (K/Rb of 100-150, [52]) and the Rolava River sediments (100-140, [15]). The Rolava River is a tributary of the Ohře River with granitic rocks in its catchments [15]. In the lower reach of the Ohře River, where the geogenic anomalies in the upper and middle reaches of the river catchment are compensated by mixing with the Cretaceous and Quaternary sediments, the K/Rb values range from 250 to 300 [15], which is much closer to the UCC value. The indirect proportionality of Sn and K/Rb in the Žatec area ( Figure 8B) is interpreted to be due to the mixing of Sn-rich and extremely low K/Rb materials, presumably those from the Krušné Hory ore regions ( Figure 1A), and some local uncontaminated material with low K/Rb ratios. The minority pattern of the sediments labelled "Sn-contaminated, unchanged K/Rb" in Figure 8B  Concentrations of Sn are elevated in conjunction with those of Cu, Pb, and Zn (only Cu/Fe vs. Sn is shown in Figure 8C), which supports their attribution to ore mining and processing in the Ohře River catchment. The Cu/Fe ratio is plotted in Figure 8C instead of the raw Cu concentration to minimise the grain-size control of the sediment contamination by potentially mobile divalent metals [53]. Mining of Sn and other metal ores in the Ohře River catchment peaked in approximately the same period (16th century) for all the mined metals. Their concentration depth profiles are similar, although each of the exploited metals (Cu, Sn, and Pb) was present in specific minerals and distinct ore regions [32,37]. For instance, Cu was mined near Kraslice, Pb near Oloví, and Sn in all other ore regions, as shown in Figure 1A.
Typical depth profiles of Sn concentrations in the fine Ohře floodplain sediments sampled by the grove corer are shown in Figure 9. The entire thickness of these sediments was contaminated by Sn in the lowest-lying floodplain segments 0a and 0b ( Figure 9C), particularly in the direct vicinity of the depressions of the abandoned channels. Weaker Sn contamination was limited to the surface sediment strata at level 0b ( Figure 9B), probably due to contaminant dilution during floods, which levelled distal floodplain. In several exceptional cases, thick Sn-contaminated strata were found in abandoned channels of the development stages Ip (sediment fill of Alte Eger, drill cores Břežany-1, Břežany-2, and NS-16), and IIm (HR-3).
The Sn contamination facilitated the identification of OSL ages departing from the general trend shown in Figure 7. The Sn contamination was not identified in the CIN-4 sediment core ( Figure 9A) in the fill of the Číňov palaeomeander scar, which is stratigraphically below the OSL-dated samples CIN-1, CIN-2, and CIN-3. We consider the lack of Sn contamination in CIN-4 to evidence its pre-16th century age, including the last channel fill (meander scar). This supports the young ages of CIN-2 and CIN-3 (Table 3, Figure 7) as evaluated by OSL is not primary, probably due to subsequent sediment reworking (scenario 1a, Table 4).  Such mechanism has not been described yet The study area has settings resembling alluvial fan, where authigenic cycles (avulsions, reworking) are common Further study is needed

Distribution of Sn-Contamination in Stage 0 Sediments
Detailed examination of the Sn contamination was conducted in the channel belt 0 upstream of Žatec in the floodplain segments ( Figures 1B and 10) by in-situ XRF analysis. The sampling provided insight into the thickness of the total and contaminated sediment strata of the abandoned channel fills or overbank fines, as only such fine sediments were retrieved by the manual grove corer.  Figure 1B.
In Stroupeček ( Figure 10B), meander scars were identified in the DTM and the field, where 70-120 cm of fine sediments was retrieved. All sites in the depressions showed Sn contamination at depths of at least 70 cm, while the higher bar surfaces (STR-12 and STR- 14) were not contaminated, apart from a thin surface layer of flood deposits. At this site, elevated parts of floodplain 0 were thus cut and laterally filled with sediments before the Sn mining peak.
Approximately 100 cm of fine sediments was retrieved in Stroupeč ( Figure 10C), with most profiles strongly contaminated by Sn to depths more than 50 cm. Sediments in shallow depressions (meander scars: STR-3 and STR-7) were contaminated to depths of more than 120 cm, i.e. the entire fine stratum. Small loops of meander scars and lack of scroll bars in both Stroupeček and Stroupeč indicate a braided system similar to the Libočany channel represented on historical maps ( Figure 5B) rather than active meandering, which would have produced scroll bars or larger circular meander scars with cut-off necks. Their fills in STR-3 and STR-7 were all younger than 16th century.
In Záhoří ( Figure 10D), typically only ca. 60 cm of fine sediments were found, wherein moderate contamination by Sn occurs in the top ~40 cm. Several profiles, such as ZAH-4, ZAH-5, and ZAH-12, were contaminated by Sn at depths greater than 90 cm, indicating a channel that was abandoned during the Sn mining peak or soon afterwards. Only one core (ZAH-11) showed nearly no Sn contamination, which clearly represented a remnant of an older surface, similar to that at higher levels in the 0 channel belt in the Stroupeček site.
In Libočany, the channel belt around the multi-channel system depicted in the digital elevation model and on historical SK maps (Figures 4 and 5B) was strongly contaminated by Sn. Floodplain topography resembled that in the left bank sites (Stroupeček, Stroupeč, and Záhoří) depressions remained after the former low-sinuous channels with large bars in the channel. The distribution of Sn contamination showed lateral channel shifts that started before the mining peak and continued during that peak. The same pattern was found in the lower reach of the Ohře River in the Písty area [16]. In the Písty area, remarkable scroll bar structures documented the enhanced lateral floodplain reworking, which started just before the mining peak and continued during or soon after this peak, as evaluated by the Sn contamination pattern [16].
In valley floor below the site Horkau, located in the 0b channel belt upstream of Libočany and on the right side of channel belt 0 opposite to Stroupeč ( Figure 1B), showed only a thin layer of floodplain fines (a few dm) over a cobbly and pebbly channel bottom deposits, while all the fine deposits were contaminated by Sn. This observation documents the limited deposition (or subsequent erosion) of finer sediments in the low levels of the 0b channel belt.

Channel Pattern of the Ohře River Since the Late Pleistocene
The geomorphic evaluation (Table 1) and OSL ages (Tables 2 and 3) show that the Ohře River had meandering pattern during the LGM, P/H transition, and a major part of the Holocene. The palaeomeanders from the LGM and P/H transition (stage II-III) were large, similarly to some other European rivers over the time periods [54][55][56]. However, the channel-forming discharges of the LGM rivers were likely achieved only for a few days or weeks of ice or snow melting in their catchments [57]. Later, during the Holocene, the channels of the Ohře River (stages I and 0) became much narrower (40-90 m) than those of the large P/H meanders (stage II, 130-200 m, Table 1), and the mean radii decreased substantially (>500 m in stage II and 250 M in stage I, Table 1). These trends have been found universally in European rivers around the P/H transition [54,55,58]. Thus, during the Middle and Late Holocene the discharge of the river dropped (stage I), though it reached the smallest values during historical times (stage 0). Finally, in historical times the sinuosity of the channel dropped, and it incised into the valley floor. This incision was completed by river engineering in the 20th century. Before the 20 th century multiple channel sections developed, however as the result of engineering works the active channel became simple with low-sinuosity.

Incision of the Ohře River Since the Late Pleistocene
Two main types of fluvial surfaces were identified in the area. On higher elevations, along the rims of the valley patches of flat surfaces remained (stage III). For these surfaces the typical terrace development scheme could be applied [12,13,[24][25][26][27]: they were formed by fluvial aggradation during an aggradational terrace development stage, and they bear some remnant of large palaeomeanders. However, younger surfaces (stage I-II) are related to incised meanders: their meander cores reflect the elevation of that flat surface, which existed before the slow incision of the river, while the related palaeomeanders reflect the termination of a given incision phase.
The incision phase in the Ohře River valley occurred between the end of the LGM and the Early Holocene (i.e. around the P/H transition). The OSL data from various palaeomeanders reflect, that erosion was a dominant process during the entire studied period, only very short (yet undated) aggradational periods could appear. Panin et al. [59] dated the last major incision phase in the Seym (Seim) River (Ukraine/Russian transboundary river) to a period between the LGM and the beginning of the Holocene and supported their findings with reports on other European rivers; the overview of the central European rivers [60,61] confirms that model. Kasse et al. [55] concluded that the major incision phase of the Roer River occurred between the Younger Dryas and the Early Holocene. The ingrown (incising) character of the LGM meanders of the Ohře River is obvious from the declining surfaces of meander cores (Figures 4 and 6), and can be found in all big palaeomeanders of stage II or the remnants of their cores in stage III, particularly in Číňov, Nové Sedlo II, Rybňany-Zálužice, and Hradiště-Strkovice. Occasionally, only the erosion remnants of meander cores have been preserved, while their lowest channel fills were eroded (Zlatník Hill, Figure 4; Rybňany-Zálužice palaeomeander, Figure 6).
The considerable Holocene incision and formation of a novel valley floor is indirectly supported by the distribution of bones and teeth of large Pleistocene mammals, such as rhinoceros and mammoths, attributed to the Palaeolithic hunters. They were found mainly on edges of the gravel terraces south of Žatec [35], i.e. terraces IV or higher based on notation by Elznicová et al. [16]. These levels are much higher than those addressed in this paper. One notable exception was a single locality above the former Dreher brewery ( [35], pp. 113 and 114), just above the former entry to the Žatec palaeomeander. This locality has an elevation between 200 and 210 m a.s.l. (i.e. between levels I and II), and is thus assigned to the P/H transition ( Figure 2B).
The developed surfaces could be considered as strath terraces, as the river incised into much older valley-floor deposits. The bedrock of the well-preserved ingrown meanders of the Ohře River are very cohesive Miocene sediments, mainly clays and coaly clays. The mechanical resistance of these rocks made the outer slopes of the meander steep and smooth ( Figure 11F). The meander cores are covered by cobbly, gravelly, and sandy lateral deposits with a typical thickness of ca. 3-4 m according to the CGS database of drill cores [16]. The Pleistocene sediment layers are therefore thinner than the entire elevation difference of the declining surface of the meander cores (5-10 m, Figure 5). In particular, the stage II palaeomeanders of the Ohře River are genetically closer to the strath terraces than to cut-and-fill terraces, the latter being easier to date due to a larger amount of deposits and their protection to aggradation.
The event-driven large-scale reworking of sediments is typical for alluvial fans [54]. The settings of the Žatec area seem favorable for the formation of an aggradational body, as the V-shaped valley of the narrow gorge at Libočany controls the slope and sedimentation, and the propagation of fluvial processes especially during extreme floods. Upstream of the study area (upstream Kadaň, Figure 1A) the high slope and V-shaped river valley provided a limited accommodation space for sediments [15,32,33], in contrary downstream of the "bottleneck" the decreased slope conditions had favoured sediment deposition. Therefore, the Žatec area represents the sediment accommodation zone immediately downstream of the sediment transfer zone. The Libočany gorge was probably even relatively narrower to the upstream valley segment during the formation of large palaeomeanders of development stage II. Assuming that the palaeomeanders II preserved on the right side of the valley (Číňov and Nové Sedlo) formerly had their symmetric counterparts on the left side of the valley that are currently missing, the valley width was ca. 2000 m upstream of Libočany and only 700 m in the Libočany gorge, with the former valley width being estimated from the current topography at a 30 m elevation above the lowest level 0. The existence of the Libočany gorge could have resulted from the mouth of a side valley of the Liboc Creek, which springs in the Doupov Mountains. The inputs from side valleys can contribute to the sediment plug formation in the Ohře River channel [32]. This effect could have favoured Late Glacial sediment deposition upstream Žatec.
Loess deposition could also partake in the late glacial development, as has been shown for other European river valleys during the Late Glacial [55,59,62]. Loess deposits were found in higher terraces of the Ohře River ( [20], CGS maps) and probably in levelled surfaces of stage III and higher, where meander scars are sparsely visible [16]. Sediment accumulation in the Late Glacial followed by erosion and channel incision before the P/H transition and Holocene was a phenomenon typical of periglacial rivers [22,55,57,61] and presumably also for other rivers globally (e.g. Portugal, [63], Australia, [64]). The expected presence of voluminous sediment bodies in the Žatec area could depict periods of repetitive lateral and vertical deposition and reworking in scenario 3.  lines in B and C). Staňkovice (C) and Hradiště-Strkovice meander cores (D) have nearly triangular shape rather than expected circular meander core. Edges of palaeomeanders cut into Oligocene to Miocene loose sandy sediments in Stroupeček are irregular (E), while in Selibice they are smooth (F) as the result of more cohesive Miocene clayey strata. Position of areas A-F are shown in inset (G).

Peculiarities of the Ingrown Meanders in Terrace Development
In the study area the remnants of meander cores (Figures 11A and 11B) are common, they belong to the development stages I and II. They could be identified as 'floodplain islands' [50], 'hillocks' [65], or 'cliffs' in the Ohře River floodplain, which are elongated or nearly polygonal elevations with gently descending surfaces, scattered on the present valley floor, such as Zlatník hill (Figures 4 and 12). These 'cliffs' can be remnants of meander cores created by neck-cutoffs, the shape of which can be inferred by extrapolating the 'cliffs' to opposite sides of the valley. An example of such an extrapolation is shown in Figure 13; the hypothetic meander neck shown here would explain the existence of transverse meanders of Nové Sedlo I and Alte Eger (Figure 4). Its cut-off would rationalise the incision of the channel belt stage 0 north of Nové Sedlo.
The meander cores usually have gentle slopes towards their palaeomeander scars, but their distal sides could be heavily eroded by subsequent lateral erosion ( Figure 11C). These features were formed in the material of an older valley floor. Although the oldest incised meanders belong to stage II, their meander cores were probably developed in the alluvium of level III, deposited by the Late Pleistocene meandering Ohře River or in the alluvium of an older, probably braided river. Our hypothesis on braided pattern is based on common fragments of low-sinuosity and relatively narrow palaeomeander scars on tops of meander cores (indicated by straight black lines in Figures 11B and 11C). Braided channel pattern was more frequently reported in European rivers around the LGM, compared to meandering [23,[58][59][60], thus, we hypothesise that the Late Glacial meandering of the Ohře River started on the braided plain and was associated with incision, as proposed by the general model from Vandenberghe [23].

Peculiarities of the Holocene Evolution of the Ohře River Valley
The Holocene channels laterally reworked some of the palaeomeander scars, such as in the Hradiště-Strkovice meander, of which core was shaped during the LGM (HR-1 and HR-2, Table 1). However, its lowermost point bar and meander scar were filled during the last millennium (HR-3, Table 3), when it was finally abandoned by the formation of the current channel belt 0. The scars subsequent to the narrower wandering channels of the middle and late Holocene can be distinguished at the bottom of the last palaeochannels of development stage II (Figures 11B-D).
Level Ip upstream of Žatec hosts numerous scars after abandonment of river channels ( Figure 4) with no visible traces of active meandering, such as scroll bars or oxbow lakes after the meander neck breakthrough and cut off. The Holocene channels upstream of Žatec, for example, Alte Eger, is a loop elongated nearly perpendicularly to the valley axis (transverse meanders or loops). Transverse meanders typically occur for rivers incised in rocks that are too resistant with respect to the actual channel forming discharge. Thus, softer material in lithological heterogeneities or faults in resistant rocks is preferentially eroded, while the meander neck cut-offs are exceptional [6]. In the case of the Ohře River, the Miocene clays and coaly clays in the Late Glacial meander cores represented highly resistant rocks, while the Holocene fluvial sediments in the meander scars represented the softer material that was easier eroded or laterally reworked. Sandy Oligocene to Miocene deposits were also much less cohesive. They occasionally outcrop in the left (northern) slope of the river valley (Stroupeček and surroundings). Meanders with irregular slopes and gully erosion features and landslides in the meander edges (DTM in Figure 11E) were cut into Miocene sandy clays and sands according to the CGS maps. It suggests, that bedrock heterogeneities affected the fluvial development of the Ohře River in the study area during the Holocene, similar to the models discussed by Malatesta et al. [13].
The elevation profiles (Figures 3 and 6) show that the Holocene development, particularly upstream of Žatec, was controlled mainly by the occasional abandonment of transverse channels of stage I by breaking through the remnants of the cores of the P/H meanders (stage II). These breakthroughs resulted in the formation of the incised belt around the current low-sinuosity channel of stage 0a, with the channel slope increasing due to the natural channel shortening ( Figure 2C). Similar mechanisms (channel avulsion and shortening) resulted in a channel belt incision of the Vistula River in the area studied by Dobrzańska and Kalicki [58].
It is tempting to attribute the cut-offs of the P/H incised meander necks to the Late Holocene climatic extremes, as the average channel forming discharges in the Holocene were considerably smaller than those around the P/H transition and usual Holocene discharges were too low to shape the valley floor. Thus, the Holocene discharges were insufficient to 'finish' the P/H meander cut-offs. The major climatic extremes of the last millennium were definitely the three extreme periods of the Little Ice Age (LIA). The general challenge is in dating these assumed incision events as they mainly caused erosion. Alte Eger is a candidate for the youngest abandoned palaeochannel upstream of Žatec, because it is still not completely filled with sediments yet and the channel has even obtained its historical name. The Sn contamination of the channel fill (Břežany-1, Břežany-2) and the proximal floodplain (NS-16) shows that Alte Eger was probably abandoned during the mining peak in the 16th century, or not long before that peak. The abandonment occurred even before finishing the incision of the channel belt 0 down to the present-day level, as the contaminated sediment would not have been deposited in the Ip floodplain level. The OSL datum of the point bar of Alte Eger (NS-16, Table 3) confirms that channel developed in the 1st millennium AD.
The Hradiště-Strkovice meander ( Figure 11D) was perhaps the youngest abandoned one downstream of Žatec, as inferred from the OSL age of HR-3 (Table 3) and thick Sn contamination in the meander scar [16]. The temporal coincidence of the abandonment of the Alte Eger arm and Sn mining peak would explain the nearly simultaneous changes in the Sn concentration and K/Rb ratio. The K/Rb change could be attributed to the input of large volumes of Miocene sediments from the hypothetical remnant of the Late Glacial meander neck (Figure 13) to the Ohře River system. Then the Sn contamination and K/Rb changes could be temporally close but independent processes, which would explain the minority group of sediments in Figure 8B.
In a recent study [16], we found two remarkable palaeomeanders with clearly visible scroll bars in the Písty site 60 km downstream of Žatec, developed and abandoned in the time period shortly preceding or coeval with the Sn mining peak. These processes followed a period of stable floodplain before ca. 2 ka [16]. Similarly, as in the Žatec area, after the abandonment of these palaeomeanders in Písty, the river channel stopped meandering and turned into a low-sinuous narrow channel belt.
The latest identified steps of cut-offs/channel belt incisions in the Žatec area, Alte Eger, Hradiště-Strkovice palaeomeander, and the two scroll bar systems in the Písty area in the lower river reach [16] could thus have occurred in the first or second period of the LIA, i.e. in the middle of 15th or 17th centuries. The last reorganisation of the Morava River channel system (eastern Czech Republic) in the Strážnice area was also attributed to the LIA [11]. The possible relationship between the elevated flood activities and LIA has been discussed by several researchers in other European [9,66] or US temperate rivers [67]. The Weiße Elster River flows from the Fichtel Mountains to the west of the Krušné Hory Mountains; that river eroded loams and replaced them by gravelly sediments during the LIA (H. Von Suchodoletz, Dresden University, personal communication). The particular impact of LIA on European rivers should be addressed in further case studies, where the climatic effect can be blurred or interrelated with the changing human impact in river catchments. The general problem of distinguishing human and climatic forcing in fluvial systems stems here from co-incidence of the end of Middle Ages, which opened new technological perspectives also in ore mining and metallurgy, and beginning of LIA, that both could affect budget of fluvially transported solids.

The Anthropogenic Impact on the Ohře River Valley
The Ohře River development in the middle and late Holocene produced relatively flat, low gradient surfaces (Ip and Im) in the Ohře River floodplain, which hosted laterally shifting or meandering channels, with active meandering being limited to some areas downstream of Žatec. Further incision of these Holocene channels occurred by the occasional breaks through the meander core remnants. The term "break through" here stand for avulsion into meander neck to distinguish it from meander cut off by laterally shifting channel of actively meandering river. The Holocene narrow, low-sinuous channels of stage I (Nové Sedlo I, Alte Eger, Libočany I, Figure 4) have thus been abandoned in a stepwise manner by ongoing channel shortening and weak incision ( Figure 2B). According to the CGS core database, the current valley floor is covered by a 3-4 m channel and overbank sediments on the Oligocene to Miocene bedrock [16]. Human action could contribute to this natural process (channel shortening and weak incision) by artificial meander cutoffs, one of which was documented in the Žatec city chronicle [68], and by the construction of levees (DTM, Figure 5), weirs, and revetments that limit lateral channel erosion.
The most recent channel shortening in the Žatec area occurred between the time of drawing of the SK maps (the 1840s) and the present ( Figure 5). In 1960s and 1970s, the Kadaň and Nechranice Reservoirs ( Figure 1A) were established, which trap most particulates transported formerly from the middle reach of the Ohře River [33,69]. River damming has resulted in 'hungry waters' and downstream vertical channel erosion [51], which are in conjunction with the natural long-term incision of the Ohře River channel belt. The Ohře River has thus followed general incision trends, typical for many European rivers: loss of connectivity of the river channel and floodplain [70]. Presumably, both natural and long-term Holocene trajectories together with the engineering measures ( Figure 5) could have caused a lack of alluviation and vertical aggradation in the Ohře River valley floor, prevented the current formation of natural levees and deposition of fine silty overbank sediments, irrespective of a strong human impact in the river catchment.

Methodological Notes
The OSL ages of the surfaces of the big Ohře palaeomeanders are not always ordered according to the elevation of the sampled strata relative to the level of stage 0a. Some samples (Figures 4. and 6) show reverse orders. For example, the younger ST-1 was found at higher elevations than the older HR-2. Other samples seem too young, such as those from the slope of the Číňov meander core (CIN-2 and CIN-3), or too old with respect to their positions relative to the recent channel belt (HR-2 and RY-3). We also found nonmonomodal distribution of ages in the OSL results (Table 3, samples RY-1 and RY-3).
The elevation-age misfits can be interpreted by the scenarios presented in Table 4. Among them, particular attention should be given to scenario 3, wherein repetitive erosion and accumulation of sediments occur in stage II during the LGM, and between the LGM and P/H transition, i.e. erosion/deposition/remodelling of the meander cores. The current surface of the floodplain in scenario 3 would be an envelope of a time-series of erosion bases of several P/H channels, with the vast majority of temporarily deposited sediments eroded before the ultimate abandonment of the last channels of such large meanders. This envelope would show (in current topography) the deepest local erosion of the cohesive Miocene bedrock (covered by the Quaternary channel armour) rather than the unidirectional downslope formation of scroll bars. Previous studies have attributed distinct parts (elevation intervals) of the slopes of these meander cores to distinct glacial/stadial periods [19,20]. We assigned these distinct elevation intervals to a set of erosion bases representing a much shorter time interval during the Late Glacial. This is presumably attributable to the environmental disruptions between the LGM (ca. 22-18 ka) and the Holocene, such as the Bølling-Allerød interstadial and Younger Drays stadial.

Outlines of Future Research
The OSL dating of the Ohře River sediments provided a new perspective for fluvial research in this valley as older local studies require complete revision. Detailed insight showed that the P/H incision was not a unidirectional vertical process, resulting in plain strath terraces with one step per a glacial cycle. The meanders, formerly declared to be of the Pleistocene age and assigned so in the Czech geological maps, although no direct sediment dating has yet been available for them), were formed and abandoned during the Late Glacial and Holocene, as stated by Elznicová et al. [16] and confirmed and extended in this work. Natural processes caused channel shortening by very slow breakthroughs through the P/H meander cores during the Holocene, postponed process initiated during the LGM and P/H transition ( Figure 13). Some terrace levels of the Ohře River have not been formed as a direct consequence of major climate change, but presumably due to the lateral channel movement combined with persistent [12] and stepwise incisions of the channel belt, caused by its shrinking on advancing the vertical incision [13]. The extensive floodplain sediment reworking could have produced fragments of elevated palaeosurfaces that agree with models of the mosaic structure of floodplains due to the spatially diverse probability of sediment reworking by meandering rivers [11]. The realistic results available for the Ohře River can be rationalised in two ways: either by a systematic general incision trend, expressed as a straight line in elevation vs. age (line in Figure 7), or event driven incision, where the major abrupt events were the LGM and P/H transition, and the late Holocene abandonment of the meandering pattern (grey rectangles in Figure 7). The decision as to which of these models is more realistic would require further sediment sampling and dating or process modelling.
We hypothesise that the Holocene climatic extremes could be responsible for the most recent cut-offs of the Late Glacial meanders, in particular Nové Sedlo I and Alte Eger, Staňkovice, Stekník, and Hradiště-Strkovice palaeomeanders. In particular, the general incision trend in Figure 7 is only a general tendency, and is an envelope of the past events in reality, as discussed above. The pattern of Sn contamination indicates the last temporal cumulation of the meander cut-offs in the case of Alte Eger and Hradiště-Strkovice; the corresponding channel belt incision occurred during the LIA. Detailed work is required to improve the age constraints of these episodes to evaluate the hypothetic event-driven model for the Ohře River valley development. If proven, the Ohře River will require an interpretation that is substantially different from that of the PDF discussed in the Introduction. The case study in the Ohře River documents inherent complexity of the fluvial system, that prevented simple strath terrace formation during the global environmental change (P/H transition) and resulted in the P/H incision fragmented to the entire Holocene by inherited LGM topography, heterogeneous bedrock geology, and decreased channel forming discharge. These individual factors should all be considered in data interpretation, rather than attempting to "automatically" summarise available data in the PDF.
The results presented above could impact the decision making for the future restoration of the Ohře River. Brown et al. [17] declared a need to understand the river behaviour before restoration activities. For European rivers, Brown et al. [17] considered the trajectory from sediment-starving swampy floodplains with laterally stable river channels during the early to middle Holocene to clastic floodplains aggraded by anthropogenic siltation and channel meandering during the late Holocene as valid. The Ohře River has followed an opposite scenario, with persistent Holocene incision and change in the naturally meandering river to a narrow, incised channel belt, with a lack of siltation despite mining activities, large-scale deforestation in montane parts of the catchment, and intensive agriculture in the entire catchment. The Ohře River development underpins the idea by Cordier et al. [1] which states that understanding the current river behaviour, its trajectory must be followed not only for a few centuries, but for the entire Holocene and beyond. If the river restoration practices in Europe would mean the renewal of connectivity between the floodplains and channels, as discussed by Maaß and Schüttrumpf [69], this goal would contradict to the general Holocene evolution in the Ohře River.

Conclusions
The palaeomeanders in the Ohře River valley are much younger than previously expected and actually presented in the maps of the Czech Geological Survey that prompts for revisiting preceding fluvial geomorphic studies not based on state-of-the-art knowledge and instrumental dating. The palaeomeanders document dynamic and spatially heterogenous development of the river valley between the Last Glacial Maximum and Holocene. The lack of anthropogenic aggradation in the Ohře River floodplain is caused by persistent river incisions that occurred since the Last Glacial Maximum or the Pleistocene/Holocene transition. The termination of this incision was postponed to the present time by the decrease in channel-forming discharges during the early Holocene. The incision was slowed down by much lower erosive power of the Holocene river compared to the Late Glacial; the last steps of that were postponed to the Mediaeval and modern periods, when the channel belt shrunk to its narrowest mode in the last 30 ka and the river has stopped active meandering, presumably due to the co-incident natural forcing and perhaps also anthropogenic impacts. The Holocene evolution of the Ohře River valley floor consisted of a series of incision events, which were supposedly triggered by climatic extremes and reflected inherited Late Glacial morphology and bedrock heterogeneity. The valley floor development was far from simple models of geomorphic record of the climate changes.

Data Availability Statement:
The data presented in this study are available in Table 3.

Acknowledgments:
The authors thank to M. Hošek, M. Faměra (both Institute of Inorganic Chemistry CAS), J. Elznic, and F. Elznic, who assisted in some fieldworks. Laboratory sample processing and XRF analyses were performed by P. Vorm and M. Maříková (Institute of Inorganic Chemistry CAS).

Conflicts of Interest:
The authors declare no conflict of interest.