Straight to Low-Sinuosity Drainage Systems in a Variscan-Type Orogen—Constraints from Tectonics, Lithology and Climate

A holistic-modular approach has been taken to study the evolution of three straight to low-sinuosity drainage systems (=SSS) in an uplifted basement block of the Central European Variscides. The development of the SSS is described by means of a quadripartite model. (1) The geological framework of the SSS: Forming the lithological and structural features in the bedrock as a result of different temperature, pressure and dynamic-metamorphic processes. (2) Prestage of SSS: Forming the paleo-landscape with a stable fluvial regime as a starting point for the SSS. (3) Proto-SSS: Transition into the metastable fluvial regime of the SSS. (4) Modern SSS: Operation of the metastable fluvial regime Tectonics plays a dual role. Late Paleozoic fold tectonic creates the basis for the studied SSS and has a guiding effect on the development of morphotectonic units during the Neogene and Quaternary. Late Cenozoic fault tectonics triggered the SSS to incise into the Paleozoic basement. The change in the bedrock lithology has an impact on the fluvial and colluvial sediments as well as their landforms. The latter reflects a conspicuous modification: straight drainage system ⇒ higher sinuosity and paired terraces ⇒ hillwash plains. Climate change has an indirect effect controlling via the bedrock the intensity of mechanical and chemical weathering. The impact on the development of the SSS can be assessed as follows: Tectonics >> climate ≅ bedrock lithology. The three parameters cause a facies zonation: (1) wide-and-shallow valley (Miocene), (2) wide-angle V-shaped valley (Plio-Pleistocene), (3) acute-angle V-shaped valley (Pleistocene), (4) V-shaped to U-shaped valleys (Pleistocene-Holocene). Numerical data relevant for the hydrographic studies of the SSS are determined in each reference area: (1) Quantification of fluvial and colluvial deposits along the drainage system, (2) slope angles, (3) degree of sinuosity as a function of river facies, (4) grain size distribution, (5) grain morphological categorization, (6) grain orientation (“situmetry”), (7) channel density, (8) channel/floodplain ratios. Thermodynamic computations (Eh, pH, concentration of solubles) are made to constrain the paleoclimatic regime during formation of the SSS. The current model of the SSS is restricted in its application to the basement of the Variscan-Type orogens, to an intermediate crustal maturity state.


Introduction
Fluvial drainage systems have been in the focus of geomorphologists from different backgrounds, such as geology, geography and others. Their studies cover a wide range of different targets, from morphological studies of the fluvial drainage patterns to the degree of sinuosity and the measurements of the fluvial discharge [1][2][3][4][5][6][7][8][9][10] Studies of the straight and low-sinuosity channel systems (=SSS) are often eclipsed by more prominent braided, anastomosing, and meandering drainage systems. This lack of emphasis is also true for the studies on the control of endogenous parameters which have always been sidelined in favor of exogenous ones, such as climate change [11]. Structural analyses, particularly in neotectonics, are not anything out of the ordinary in geomorphology but frequently related to issues of civil engineering and sister drainage systems of the SSS [12][13][14][15][16][17][18][19][20][21][22][23][24][25]. A broader view has been taken on the evolution of fluvial drainage system as a function of tectonics in some comprehensive studies [17,[26][27][28][29][30][31][32][33][34]. Similar comprehensive studies dealing with the mineralogy and petrography of straight to low-sinuosity drainage systems cannot be quoted for the Variscan basement in Central Europe. Therefore, in the current terrain analysis emphasis is placed on a combination of geomorphology with tectonic, petrographic and climatic issues. The study area selected has one basement region with similar tectonic history and lithologies where the geological and geomorphological influences on the evolution of the 3 SSS can be investigated. Such an ideal study area is offered by the NE-Bavarian Basement, at the NW edge of the Bohemian Massif (Figure 1a,b).
One can find there three geodynamic units, which are drained by rivers heading towards the W and tributary to a common trunk river system, the Rhein-Main fluvial network. The hydrographic setting enables them to take a catena-like approach from the source (watershed) to the depocenter in the adjacent foreland of the drainage system. The three reference sites, albeit situated side-by-side with each other in the basement, are markedly different from each other with regard to their histories of the P (pressure) and the T (temperature) regimes both of which are decisive for the way composition and deformation of fluvial bedrocks can impact on the SSS. The changing physical parameters during the Paleozoic orogeny provoked a change of the bedrock lithology and the type of structural geology that influenced the formation of the landscape during the Cenozoic (see Sections 5.2 and 5.3). Consequently, these physical parameters form a good basis for a tripartite subdivision: (1) Frankenwald (very-low grade, VLG-type), (2) Münchberg Gneiss Complex with marginal facies (low-to high-grade, LHG-type), (3) Fichtelgebirge (low-grade with granitic intrusions LGG-type) [35] (Figure 1b,c, Table 1).
The uplifted basement block is part of the Central European Variscides framed by two sets of planes. It is bounded towards the Mesozoic foreland by a deep-seated subvertical lineamentary fault zone called the "Franconian Line Fault Zone" which herein is abbreviated to FLFZ and a concealed NW-SE structure forming the divide between the Rhein-Main and Elbe/Labe drainage basins and, hence, demarcating the study area to the NE. In vertical the crustal section is truncated by subhorizontal paleosurfaces arranged in a staircase-like way that play a major role in the fluvial evolution [38,39] (Figure 1b). The "stairs" demarcate the lithological outcrops within and the height of "stairs" control the gradient and longitudinal profile of the fluvial drainage systems (Section 5.2.1).

Figure 1.
Position and the geological setting of the study area in SE Germany with three type localities framed by stippled line: 1. VLG-type area, 2. LHG-type area, 3. LGG-type area. For features of each reference site see Table 1. (a) Position of the study area in the Central European Variscides in SE Germany at the NW edge of the uplifted basement of the Bohemian Massif (BM) and at the NE margin of the South German Block (SGB) representative of the Mesozoic and Cenozoic foreland basin. Rivers, towns and tectonic structures as well as geology are shown on the right-hand side, (b) The geological setting of the three type localities at the NW edge of the Bohemian Massif. Geological basis [35] (c) The position of the study areas referred to in the text in relation to the geomorphological-geodynamic setting and the distribution of Variscan-type orogens Maturity of crustal section: green: pre-mature, orange: mature, red: super-mature [36]. Table 1. Synoptical overview of the three straight to low-sinuosity reference drainage types as a function of the metamorphic grade and presence of intrusive rocks in the bedrock lithology. The three type were categorized according to the bedrock lithology the minerals of which decrease from left to right and the tectonic structures get obliterated in the same direction. both of which are besides climate the key elements of the origin of the SSS and form the framework stage. The rows describe the dual approach taken in the current study (see Figure 2). Blue: Geoscientific anamnesis and visual examination. Yellow: Geoscientific measurements, analysis, and computations. The "Basement region" denotes the three areas at the western edge of the Bohemian Massif ( Figure 1). The "Basement and foreland lithologies" provide an overview of the lithologies shown on the geological maps 1:25,000, a compilation of which is given in Figure 1. The landform series are described according to the proposals of Easterbrook (1999) and Bridge (2003) [8,37]. Approx. 300 samples/station points have been analyzed lithologically and mineralogically. Degree of sinuosity has been calculated according to the equation referred to in the text for each river section in the maps following the recommendations of Charlton (2014) [10] and using the ArcGis software 10 Single-channel to antler-shaped See Table 2 Dentritic to parallel See Table 2 Dentritic to antler-shaped See Table 2 Sediment Petrography Grain size distribution (%) Clay + silt: 12.  Table 3b See Table 3a See Table 3a Grain orientation ("situmetry") See Table 3b See Table 3a See Table 3a The study deals with a drainage system characterized by SSS and it is aimed to provide answers to the following questions: • What is the impact of the different types of structural geology on the fluvial systems? • What is the impact of the varied bedrock lithologies on the hydrography and the hydrodynamic evolution? • What is the impact of Neogene-Quaternary climate on the drainage systems? • To what extent can the GMS gravel (granulometry-morphology-situmetry) analysis contribute to the hydrodynamics of coarse-grained fluvial deposits? • Is the elaborated genetic model of the fluvial drainage system applicable to all landform systems in orogenic realms or geodynamically restricted to a special state of crustal maturity?
The current study is designed to account for the genesis of the SSS in a basement terrain and to provide a manual for a holistic-modular geomorphological approach. The modular workflow chart of Figure 2 forms the platform for the investigation and interlocks methods, results, and discussion.  Flow charts of the holistic-modular approach taken to study the SSS in a basement landscape. It is a dual-use approach flow chart placing emphasis on the methodology (a) and the origin of the SSS (b) with the land-forming processes linked by color codes. (a) The flow chart to link methods (descriptive diagnosis -upper part, processes center part, numerical methods lower part): Multicolored boxes: processes relevant for the SSS (see also color code in (b)); White boxes: Methods; Green boxes: Tools and products; Yellow boxes/framed in dark green with Arabic numerals reference to "Results". (b) The genetic flow chart to link genetic parameters as a function of evolutionary stages (framework, prestage, proto-SSS/CA, modern SSS/CA, modern SSS/TA, modern SSS/DA): Multicolored boxes: processes relevant for the SSS (see also color code in (a)); Blue boxes framed in red: Tectonic drivers and their age of operation; Beige boxes framed in red: Climatic drivers and the age of operation; Yellow boxes/ framed in dark green with Arabic numerals reference to "Discussion".

Strategy of the Flow Chart
The current paper takes a multidisciplinary approach, encompassing methods from geomorphology, sedimentology, structural geology, mineralogy, geophysics and radiometric age dating (Figure 2a). They are focused on landforms that are defined as physical recognizable forms or features, e.g., plains, valleys, slope, and hills of the Earth's surface, having a characteristic shape and produced by natural causes [40] (Figures 3-5).  Table 1-VLG-type) in digital terrain models (DTM) The V-shaped valleys are highlighted by the slope angle shaded in red (see angle bar = midslope). Rectangles framed in red denote landforms of fluvial relevance. (a) Geology and the geomorphological setting in the upper reaches (headwaters = (CA)) transitional into the proximal transport area (TA)); (b) Excerpt of the topography and the level of paleosurfaces of (a); (c) Geology and the geomorphological setting in the lower reaches (distal TA transitional into the depositional area (DA) within the basement and proximal foreland); (d) Excerpt of the topography and the level of paleosurfaces of (c).  Table 1-LHG-type). For legend see Figure 3; (a) Geology and the geomorphological setting. Blue stippled line marks the fan onto the hillwash plain; (b) Excerpt of the topography and the level of paleosurfaces of (a).  (Table 1-LGG-type): (a) Geology and the geomorphological setting; (b) Excerpt of the topography and the level of paleosurfaces of (a).
The classical geomorphology places a strong emphasis on the description of the Earth s surface [3], whereas more recently emphasis has also been placed on measurements using digital terrain models (DTM) and different types of maps and x-y plots to achieve solutions for issues of shaping the landscape [41]. Both approaches are likewise important in providing solutions to geomorphological issues but they should not cast in the role of stand-alone-procedures. That is the reason why a holistic approach has been taken in the current study to fill in the gap between descriptive and numerical methods relevant for landform evolution. It is a sort of a "semi-automated" technique used by the senior author in applied geomorphology and economic geology in central E Africa [42].
A holistic-multidisciplinary study cannot be performed in a simple way and presented as a unidirectional workflow. The various parameters, disciplines, and methods as well as processes and landforms are closely interlocked with each other so that the links between them have to be illustrated at the beginning in a special diagram. The qualitative-quantitative approach taken in the current study is illustrated in a binary flow chart with multiple links on display and the same codes numbering the (sub)sections in the text are used in both flow charts, too (Figure 2a

Geoscientific Mapping and Visual Examination
To create a geological and petrological mapping platform on the scale 1:25,000 is a prerequisite of any field study and of further investigations. In combination with a topographic map a composite digital terrain model (DTM) is created for each of the three study areas at the same scale (Figures 3-5). The XYZ data were obtained from "Geoportal Bayern", at a 50 m resolution, for the whole of Bavaria, in a GeoTIFF format. The GeoTIFF data were then cropped to fit the area of interest using QGIS and subsequently processed in Golden Software's Surfer to obtain a virtual 3 DTM by using the kriging method of gridding and data interpolation.
It can be used for all further geological-petrological-geomorphological terrain analyses throughout the entire study and assists in visual examination of geomorphological landforms in the field, their mapping (e.g., depositional terraces, floodplains) and the establishment of different relief generations (Figures 3-7). Furthermore, this DTM and the pertinent geological and petrological maps aid in delineating the supergene alteration zone and subdividing the regolith zone into duricrusts and argillitization ( Figure 2). The first principal part named "Geoscientific anamnesis and visual examination" results in the subdivision of the three reference landform series colluvial, cryogenic and fluvial. The latter landform series is categorized into master rivers (1st order) and tributary rivers (2nd order). The 1st order rivers, Rodach (VLG-Type), Steinach (LGG-Type), Schorgast-Coserbach drainage systems (LHG-Type) are in the focus of the current investigations and used in the presentation of the results and the discussion unless stated otherwise.  (Table 1). For peneplains, erosional and depositional terraces, and floodplains see Figure 3.

Geoscientific Analysis, Measurements, and Computations
The second part of the workflow chart denominated as "Geoscientific measurements, analysis, and computations" is split into two subsections, one dealing with indoor-only (laboratory) methods and another one working at the passage from field to laboratory.

Laboratory Methods
Dependent on the practical experience and knowledge a great deal of geoscientific identification and classification work in the field can be performed with the unarmed eye. For fine-grained, earthy and complex mineral aggregates, however, supplementary laboratory work is necessary (Table 1).
Petrographic microscopy: Routine thin section analysis was carried out in combination with optical microscopy of heavy and light minerals with 200 to 300 grains counted per sample. Heavy mineral separation was carried out for each sample to the grain size fractions richest in heavy minerals (the 3-4 ø fraction) using sodium polytungstate of a density of 2.9 g/cm 3 [43]. Approximately 70 samples were subjected to this mineralogical analysis.
X-ray analysis: X-ray diffraction (XRD) patterns are a supplement to the afore-mentioned method and were recorded using a diffractometer (Cu-Kα radiation generated at 40 kV and 40 mA), equipped with a 1 • divergence slit, a secondary monochromator, a point detector and a sample changer (sample diameter 28 mm). The samples were investigated from 2 • to 80 • 2Θ with a step size of 0.02 • 2Θ and a measuring time of 3 s per step. For specimen preparation the top loading technique was used.
Electron microscopy: Scanning-electron analyses (SEM-EDX) were carried out under low-vac-chamber conditions (1 to 10 mbar), no sputter coater was used prior to analysis. It is run in combination with the XRD analysis.
Raman spectroscopy: The micro-Raman spectra were acquired by using the Renishaw InVia Raman Microscope system, equipped with two laser sources with wavelengths of 532 nm and 785 nm and full power of 500 mW, with a CCD detector. The measurements were performed with the 532 nm wavelength laser with 1800 L/mm holographic grid, at a power of 10%, exposure time of 1-5 s, 10-20 acquisitions, using the zoom objective 50× from the Leica microscope attached to the spectrometer. The interpretation of the micro-Raman spectra was done using the LabSpec software by means of baseline correction and peakfitting. It is the most advanced level and used for fine-tuning of the mineralogical results.
The mineralogical analyses are applied to the bedrock of the channel system during hydrographic analysis together with measurements of the elevation and morphology of the valley of the SSS (Figure 2a). The current data obtained are also of assistance during petrological fine-tuning of gravel-sized debris in the course of the hydrodynamic analyses as well as useful for thermodynamic computations of the Eh, pH and temperature, all of which are essential to constrain the paleoclimate regime operative during the Late Cenozoic when these SSS came into being [44]. The latter goal can only be achieved in context with precise age dating mentioned in the succeeding paragraph.
Age dating: For the chronological studies of geomorphological processes radiocarbon dating and the OSL dating is common practice [45,46]. Both methods do not cover the period of time from the late Pleistocene through the Neogene relevant for the current investigation. Therefore, only the classical U/Pb and K/Ar or Ar/Ar age dating methods appropriate for K-bearing manganese oxide-hydrate and uriniferous silicates and phosphates can successfully been used [47][48][49][50]. In a few places where syngenetic organic matter has been identified in Neogene sediments palynology may be useful [51].

Field and Desk(top) Methods
Structural analysis: To investigate the spatial relation between structural elements mapped in the bedrock and the orientation of the SSS, the strike of faults, joints, cleavage planes and fold axes is determined and statistically treated in rose diagrams in plain view. In the course of this study, the areas under consideration are subdivided into litho-and morphotectonic units, each of which is characterized by a full-circle mirror rose diagram with 10 • sectors displaying the frequency distribution by percentage values (Figures 8-10). By contrast, the overview diagram of the entire VLG-Type area is plotted as a semicircle rose diagram with 20 • sectors, a technique which has been tested statistically to be sufficient for the different scales 1:25,000 and 1:200,000 ( Figure 8e).  Hydrographic analyses: The hydrographic analysis makes use of x-y plots (1-D), drainage and contour maps (2-D) and in one case also of DTM (3-D). The longitudinal profiles of the SSS reference drainage systems are plotted as function of the bedrock lithology in x-y graphs displaying the elevation vs. distance from the watershed in meters ( Figure 11). The data have been collected from the topographic map and in the field and the lithology is determined by mapping and petrological analyses ( Figure 2a). The slope angle diagram is identical as far as the x-axis is concerned with the host graph (Figure 11a-c). The y-axis, however, is different and yields the maximum slope angle measured with a tilt sensor in the field from the divide to the footslope and calculated from the numerical data which were also used for the construction of the 3-D diagrams (Figures 3-5). The slope angle can be deduced from this sub-diagram and correlated with the three facies zones CA (zone of headwaters), TA (transport area), and DA (deposition area).     The drainage maps (2-D) are simplified geological maps only revealing the geological and hydrographic relevant items such as the bedrock which is subdivided into basement and foreland (Figures 8-10). The hydrographic system is split up into the recent channel system/ floodplains and ancient floodplains/depositional terraces. To allow for a direct comparison between the direction of the talwegs of the 1st and 2nd order channels and the tectonic features, rose diagrams are added up to these maps (see previous paragraph). The degree of the sinuosity is also integrated into these composite maps. Numerical parameters to describe the sinuosity of fluvial channel systems were proposed by Horacio (2014) and Turowski (2018) [52,53]. In the current study a modified approach has been taken based upon the ArcGis software 10.2.2 and the sinuosity is calculated as follows: ratio S = channel length/valley length ( Table 2). An ideal straight channel yields a sinuosity degree of S = 1.
The density of drainage systems has been calculated according to the equation recommended by Ahnert (2015) [54]: density (D) = total length of drainage systems in kilometers divided by the size of the area measured in square kilometers. The data have been obtained from hydrographic maps which were subdivided into a grid of squares measuring 1 km × 1 km and 0.5 km × 0.5 km. For each square the density of aquatic systems was calculated (software: Adobe Photoshop CS3). The interpolation maps were created using the IDW interpolation method (Inverse distance weighting) and contour lines were generated after the resultant interpolation raster map was handled by means of ArcGIS 10.2.2. The density is shown by a color key together with the sinuosity and the cross-sectional shape of the valley in Figure 12. In addition to this river density the ratio of between the size of the floodplain and the size of the stream has been calculated and plotted side-by-side with each other. The database of the slope angle measurements consists of field measurements measuring the inclination in some places and a topographic map on a scale of 1:25,000 with a resolution of a 5-m spacing in height of the contour lines. It enables us to achieve a high precision in x-y plots such as Figure 11 and in plan view ( Figure 13). The determination of the slope angles was calculated for the same XYZ data by using Surfer's Terrain Slope Calculus gridding function and presented in two different ways. A contour map was designed for each drainage system as exemplified in Figure 13 yielding an inclination from 0 to 40 • . The contour map generated by this slope data was then overlaid on the 3 DTM to show areas of relevant slope angle intervals, in this case for the sake of readability, only for the maximum slope angle greater than 25 • (Figure 3c,d). Hydrodynamic analyses: The hydrodynamic analyses focus on the unconsolidated gravel-sized lithoclasts of the three landform types (fluvial, colluvial, cryogenic) identified in the three reference areas (VGL-, LGG-and LHG types) and counts on the GMS tool (Granulometry-Morphology-Situmetry), listed for the VGL-and LGG types ( Table 3). The preparatory work in the field is applicable to each of the three methods. A horizontal square measuring 1 m × 1 m was delineated in the transverse and longitudinal bars of the channels and in some places on the delta plain of some micro-delta wedges produced by 2nd order rivulets debouching the lithoclasts into the 1st order river. These sampling and measuring squares can also be form an overlay for mass wasting and cryogenic products. In outcrops of cryogenic material, the same grid was used in vertical position with sectors measuring 10 • to distinguish cryoturbated from fluvial deposits overriding these pockets filled up with lithoclasts ( Figure 7g). The number of gravel clasts totals approx. 100 to 150 samples, with the same specimen used for the each of the three methods. The hydrodynamic analysis always begins with the in situ measurements of the clast orientation with a hand-held Brunton compass using as reference the longest axis of the clasts. Its orientation is plotted relative to the width of the fluvial channel. In the case of mass wasting processes, the strike of the ridge of the hill is taken as a baseline for the gravity-driven motion along the slope. There are a few studies mainly devoted to colluvial deposits and to a lesser extent to fluvial or glacial processes [55][56][57][58]. The cited researchers plotted the obtained fabric data on stereonets, e.g., on Schmidt nets (Lambert equal-area projection). According to the senior author s practical experience gathered in the field of structural geology and sedimentology as an economic geologist, the stereonet projection is a valuable tool for a wide range of applications to interpret the distortion of beds by faulting and folding particularly where rather complex polycyclic processes have to be disentangled [59]. In the current situation the rose diagram, however, is the method of choice, since vertical inclination, e.g., imbrication, is negligible. To investigate the mode of transport and deposition of clasts, the use of the stereonet method is time consuming and not economic. The resultant measurements of length and angle are treated with the software ImageJ to plot semicircles with sectors in the intervals of 20 • . In case of tectonic structures mirror rose diagrams are used for distinction. For each interval, the segments length was summed up and the values were normalized to 100% (Software: OriginPro 2016). It is sufficient to determine the clast orientation in plain view and if necessary, in vertical cross sections and illustrate the orientation in a rose diagram ( Figure 14). The base line of the diagrams is positioned at right angle to the gradient of the fluvial channel or dip direction of the slope measured in the field.
The grain size of gravel is determined by measuring the maximum length of the clast while the individual lengths are plotted into classes of histograms revealing the frequency distribution of each size interval. Most of the samples show a lognormal distribution the graph of which is represented based on a Smooth-Kernel-type distribution function using the software OriginPro 2016.
The grain morphology of gravel cannot be done by automated procedures such as CAMSIZER (only <2000 µm) and therefore, visual inspection using morphology charts subdivided into six classes by Illenberger (1991) [60] is the most suitable method.
Some sites have been surveyed by means of geoelectrical methods using DC (lowfrequency AC) geoelectric resistivity measurements according to the "Schlumberger-Array " [61]. This special four-point method enables the operator to determine the thickness and predict the lithological built-up of the overburden by the aid of the resistivities of the clay, sand and encrustations and helps to assess the mass wasting and fluvial deposits.  (Table 3a) for more detail. The subfigures (diagrams) and their numbers are referred to in the text. Table 2. A classification scheme of fluvial drainage systems based on the degree of sinuosity for environment analysis. The numerical data of the fluvial reference types developing within the foreland are taken from Dill et al. (2020a) [62]. Column 1: Classification of the channel system based on visual examination in the field and using maps and satellite images; Column 2: Geological setting in the basement and the adjacent foreland; Column 3: Degree of sinuosity maximum; Column 4: Degree of sinuosity minimum; Column 5: Bedrock lithology in the basement and in the foreland.

Chanel System
Geological Dip-, anti-dip and strike-stream in the catchment (CA) and transport (TA) areas of rivers Meta-quartzites, phyllites, Meta-felsic rocks (epigneiss), felsic intrusive rocks (low-grade metamorphic rocks) Meta-arenites, slates, greywackes (very low-grade metamorphic rocks) Table 3. The GMS (Granulometry-Morphology-Situmetry) tool applied to the fluvial drainage systems under study in the basement as a function of bedrock lithology aimed at environment and landform analysis. ø = size of gravel specimen. In the situmetric columns the 1st and the 2nd maximum are given in degrees. The values denote the angle between the longest axis of the gravel relative to the reference line 0-180 • which is equal the orientation of the width of the channels at the station of measurement. See also semi-circle rose diagrams in Figure 14.

Geological Setting
The Saxothuringian Zone is a geodynamic unit of the Central European Variscides which subsided during the Cambro-Ordovician as a rift basin NW of the Moldanubian Zone [38]. It is located at the NW edge of the Bohemian Massif, a basement block which has been shared in by Austria, the Czech Republic, Poland, and Germany where it is uplifted along the FLFZ by as much as 1000 m relative to the Permo-Mesozoic foreland (Figure 1b,c). From the Precambrian through the Lower Carboniferous, a great variety of sedimentary and volcanic rocks formed in this basin until the Visean when tectonic disturbances put an end to the basin subsidence during the main phase of the Variscan Orogeny. This scenario is true for the autochthonous Frankenwald (FW) and Fichtelgebirge (FG) (Figure 1b,c). The Münchberg Gneiss Complex (MGC) is, however, an allochthonous unit that was shifted towards the NW over the FG and emplaced by thrustal nappe movements into its present position. In the MGC, we are faced with an allochthonous synform with a dual litho-stratigraphic subdivision into the "Liegend-Serie" (footwall series) and the "Hangend-Serie" (hanging-wall series) (

Geomorphological Setting
The NE Bavarian Basement forms part of the morphostratigraphic and morphotectonic crustal units in Central Europe described by Zöller et al. (2017) [63]. The geomorphological evolution of this basement at the western edge of the Bohemian Massif covers a timespan from the Neogene to the Recent. The basic post-war geomorphological study in the course of which a large area of the basement was mapped on a small-scale topographic map and various levels of relief generations were singled out was published by Thauer (1954) [64]. This purely descriptive study was backed then by paleopedological studies conducted by Strebel (1955) [65]. The timing of the paleosurfaces and terraces encountered in the field at various heights, however, has only been approximated by morphostratigraphic correlation and lacks any radiometric dating still today. There have no chronological data published during the last 70 years relevant for geomorphological processes in the basement others than those used in the current paper. By and large, the basement was taken as a provenance area for marker lithoclasts to create models describing the evolution of the foreland river systems Main and Donau during the Late Cenozoic [66][67][68][69]. For the pertinent fluvial landforms no radiometric age data have been provided either, excluding the post-glacial period [68]. The 1st order classification scheme of the drainage systems is based upon the metamorphic grade of the bedrocks according to Bucher and Grapes (2011) [70] and the presence of granitic rocks ( Table 1). The regional metamorphism of the Variscan Orogeny during the Paleozoic created a great variety of lithologies that can be held accountable for the marked differences in the built-up of landforms during the Neogene and Quaternary. The higher the metamorphic overprinting of the bedrocks, the higher the amount of high-T silicate minerals which are highly vulnerable to chemical weathering. Thereby, the impact of the bedrock on shaping the Cenozoic landscape diminishes because it signifies a loss of structural guidelines for the drainage systems due to the reduced rock strength of the crystalline rocks flanking the valleys-see also Section 5.2.1.
On a 2nd order level, each of the three compartments is subdivided with regard to the two principal groups of methods described in the flow chart of Figure 2a (geoscientific anamnesis and visual examination, geoscientific measurements, analysis, and computations).
The 3rd order classification in the hierarchy is focused on the prevailing channel processes responsible for the tripartite zonation along the talweg. It starts off with the headwaters (=CA) in the hinderland and the provenance area of the fluvial sediments. The three drainage systems are bounded upstream by a watershed forming the natural boundary. Erosion is operative in the upper reaches of the drainage system causing a retreat of the slopes and the incision into the subhorizontal paleosurfaces. It is denoted by a sudden change in the gradient and the valley slope ( Figure 11). The CA passes downstream into the area where transportation (transportation area = TA) of the sediment load prevails over erosion, albeit intermediate sediment traps can still be delineated along the talweg. The final section is called the depositional area (DA) where as a consequence of a sudden reduction in the transport capacity of a flood flow gives rise to flood plains (Table 1) (Figures 9c, 10b, 11a-c and 12). Prior to the quantitative treatment of the hydrographic issues, the basic landforms need to be identified, interpreted, and mapped relative to their topographic position in the DTM (Figures 3-5).

Fluvial Processes and Their Landforms
Objects and purpose (CA): Unless the fluvial landforms in the study area are identified, mapped relative to its topographic position and documented in 2D or 3D plots, a sound platform to build upon does not exist for further quantitative treatment of the hydrographic and hydrological issues (Figures 3-5).
The headwaters and relic landforms: The SSS have their beginning in the CA, where a stacked pattern of paleosurfaces truncated a great variety of basement lithologies at different altitudes, the maximum height of which is reached at 940 m a.m.s.l. in the Fichtelgebirge study area [71] (Figures 3-5). The term peneplain is used here for a subhorizontal planar landform which is still sparking a lot of genetic investigations [72][73][74][75]. There is common consensus that these paleosurfaces are irregularly shaped, uneven plains, with low relief. In places, they may show a considerable relief similar to what has been named etch surfaces in the Sudetes Mts. at the eastern edge of the uplifted Variscan basement block by Migón (1997) [76], updating "Büdel s Doppelte Einebnungsfläche" [77]. The paleosurface with the largest extension of approx. 2500 km 2 has largely been preserved in the study area at an altitude between 580 to 620 m and cuts across different lithologies exposed in the FW, MGC and the FG (Figures 3-5). The fluvial dynamics during the prestage of the SSS has not been precisely proven by sedimentological features or dated [64]. Along the boundary between the FW and MGM a relic wide-and-shallow drainage system has been preserved [50] (Figure 6a). The NE-SW trending alluvial-fluvial channel system has been filled up with silcretes and pebbles and cobbles of quartz (Figure 6a). These shallow depressions passing into the CA of the SSS gave birth to their embryonic rivulets of the fluvial system under study. The channel pattern across these SSS varies from straight non-alluvial to alluvial single channels in the CA of the VLG-type into alluvial multiple channels of the LGG-type. It is called dentritic to parallel with cascades and step-pools according to Easterbrook (1999) [37]. It grades into a wide-angle V-shaped valley (Figure 6c).
The transport area (TA): The transition downstream into the TA is marked by change in the valley profile as V-shaped gorge-like valleys incise into the relic channels of their predecessor fluvial systems (Figures 3-5 and 6d,e). Albeit denominated TA, there are also river sectors with conspicuous sedimentation. In the LGG-type an intermediate sediment trap named herein "Warmensteinach Bowl" can be recognized by its wide flood plain and drastic lowering of the slope angle prior to the Steinach River's incision into the rigid parametamorphic country rocks ( Figure 5). The upstream trap of sediments, results in an intensification of the channel erosion and greater downstream sediment supply and transport [78] ( Figure 5).
The depositional area (DA): In a cross-sectional view of the VLG-type the V-shaped valley converts downstream into a U-shaped valley. It maintains its steep hill slopes constituted by quartzite and chert while the floodplain is considerably increasing its width (Figure 6f,g). In the VLG-and LGG types, the depositional behavior of the rivers only differs from each other by the absence of terraces within the basement, e.g., the LGG. Its acute-angle V-shaped valleys abruptly turn into a vast floodplain bordered by a tripartite set of terraces along the low scarp of the valley. The most elevated one is located as many as 40 to 50 m above the present-day talweg [79][80][81]. The sediment load is in both cases dominated by gravel and abrupt changes in the valley morphology as well as flow direction as the rivers during their incision made contact with faults and lithologies which put up sufficient resistance to change the course of the fluvial system. The LHG-type channel systems lack these landforms owing to the missing bedrock characteristics and a considerably high river gradient (Figures 9c and 11d). In essence, these paired and unpaired terraces formed as a result of two different modes, a fact which is especially well-documented in the LGG-type, where cut-and-fill structures exist with each other in space but different in time (see discussion).
The grain size in the DA of the LHG-type, however, is in stark contrast to the equivalent deposits mentioned previously. The MGC s gravelly and argillaceous bimodal sediments eventually are spread across a fan-shaped large plain (Table 1, Figure 4). Emmert and Weinelt (1962) [82] reported quartz, black chert, diabase, prasinite, amphibolite, gneiss and late Triassic sandstone to contribute to the gravel assemblage. The sediments gain a thickness of as much as 6 m and an outward appearance of a fining-upward grain size gradation typical of point bar units in the lower part. This sedimentary feature fades out towards the top. Similar sediments have been described from the badland sediments where confined drainage systems (fluvial deposits) give way to unconfined ramifying drainage systems with sediments called hill wash sediments spread between buttes [83,84].

Mass Wasting and Cryogenic Processes and Their Landforms
Objects and purpose: In combination with the processes mentioned under 4.1.2 mass wasting is a decisive modifier of the valley slope leveling it to different intensities, under different degrees of wetness and at different speed of motion by slide, flow, fall or creep, while static cryogenic features are markers of the glacial to periglacial climate regime related in time and space with the evolution of the SSS.
Mass wasting processes: Mass wasting processes are operative in all study areas and they are responsible for the poorly sorted sediments mantling the gently inclined foot slopes of the SSS valleys as well as forming patchy blankets on the subhorizontal paleosurfaces below and above the steeply dipping midslope (Figures 3-5-see steeply dipping midslope marked in red). In some valleys, additional granulometric, morphological and situmetric analyses, however, are necessary to distinguish the various gravelly deposits during mapping (Section 4.2.5). Numerical representation has been achieved and given in Table 1 for the CA, TA and DA of each drainage type.
The strength of the mass wasting processes can be deduced from the residual landforms exhumed on top of the most elevated hills as the debris have been washed out downstream into the TA and DA (Figures 3-5 and 7a-f). From the VLG to the LGG types the number and size of monadnocks erecting from the various paleosurfaces increases accompanied by an enhancement of the roundness of the bare bedrock at outcrop as well as their fragments scattered in the apron around (block field). Pinnacles of very angular Ordovician quartzites and siliceous slates are the most common remnants of confined soil and talus creep of the VGL-type (Figure 7a). In the study area of the LHG-type even moderate variations in the chemical and mineralogical composition of ortho-metamorphic bedrocks provoke different mass wasting-induced landforms and grain morphologies (Figure 7b-d). Metabasic rocks (eclogite, amphibolite) are found surrounded by small blockfields resembling those recorded from many periglacial study sites proximal to the outcrop of metamorphic rocks [85,86]. By contrast, outcrops of metaultrabasic (serpentinite) rocks are depleted in gravel-sized clasts but surrounded by an apron of fine-grained colluvial material instead. When undergoing chemical weathering features of a silicate karst with karren and rillen appear. These fine-grained mass wasting products clog the fluvial drainage systems of the SSS downstream in the TA with a muddy "slurry" and contribute eventually to the hillwash planes in the DA. Mass wasting in the LGG-Type area results in the well-known tors derived from syeno-and monzo-granites in a way similar to many granitic terrains elsewhere when deprived of their regolith (Figure 7e) [87,88]. These landforms and their dislodged blocks present in the upper reaches of the Steinach River attained the highest degree of roundness ever (woolsack) while being embedded into a thick layer of grus sensu [89]. Judging from the landform series observed in the CA and upper reaches of the TA the mass wasting processes is mainly due to talus and soil creep, block sliding, topple and rockfall. Earth flows are subordinate in size and intensity in the study area relative to the afore-mentioned types and only play a significant part where metaultrabasic rocks such as bronzite serpentinites are exposed in the CA (Figure 7c,d). These silica-deficient magmatic rocks are surrounded by an apron of fine-grained weathering loam which increases the rivers' suspended load when washed downslope.
Cryogenic processes: The impact of cryogenic processes on the evolution of the fluvial drainage system can be verified in the CA and DA of the LGG-type where pockets in the Lower Triassic bedrock of the Steinach River are glutted with matrix-supported gravel that has undergone intense cryoturbation (Figure 7g). In the CA of the Steinach River between 720 and 940 cirques are arranged in a staircase-like way ( Figure 5) [90,91]. They are filled today with granitic grus and dislodged boulders, separated by lips. They contribute to the bed load of the Steinach River. The concave surface expressions are interpreted in the current setting as nivation cirques rather than relic cirque glacier [92][93][94][95][96][97]. In conclusions, relic periglacial landforms exist at the beginning (CA) and end of SSS (DA).

Tectonics and Bedrock Lithology in the Study Areas
Objects and purpose: A tectonic analysis is a pre-requisite for the geomorphological follow-up study of the drainage system. It involves statistical spatial treatments of cleavage, joint and fault planes as well as fold axes and handled in the same way as it is performed for the drainage system itself to assure equivalent datasets. The tectonic part is conducted according to the basics published by Fossen (2010) and Lisle et al. (2011) [98,99].
The structural geology at the western edge of the Bohemian Massif is subdivided into two different phases both of which have strong implications on the drainage network under study (Section 4.2.2). The younger Saxonian tectonic deformation or post-Variscan tectonic, a fracture-and-fold tectonic, comprises a protracted history of extension and graben evolution beginning during the Triassic with phases to be observed still during the Neogene and Pleistocene [100][101][102]. This fault tectonic is syngenetic with the SSS at least during its latest phases. The older structural deformation in the Variscan Basement ended during the Late Carboniferous. The tectonic and lithology of the bedrock constitute the endogenous framework within which the SSS develop. There is one prominent linear fold element common to each of the three study areas striking NE-SW and plunging towards the SW homoclinal with the course of the SSS. From the VGL-type through the LGG-type to the LHG-type the fold axes, however, become more and more overprinted as a morphotectonic element (Figure 1)

Channel Network and Structural Geology
Objects and purpose: The morphotectonic analysis involves a direct comparison between the basement tectonic and spatial elements of the drainage system to figure out the marker guidelines for drainage patterns across the SSS.
Similar morphotectonic studies on the most recent tectonic movements provoking changes in the drainage network have been made at the NE edge of the Bohemian Massif, Czech Republic, along the deep-seated lineamentary Sudetic Marginal Fault [20,108]. This fault is a pendant to the FLFZ at the NW edge of the Bohemian Massif. Unlike many papers devoted to the neotectonics sensu stricto, we applied in the current investigation this structural term on a much wider scale, covering the entire Neogene and Quaternary [24,109].
The numbers of structural elements considered throughout each study have been mentioned in the captions of Figures 8-10. The length and orientation of the trunk rivers-1st order-Rodach, Schorgast and Steinach can be deduced from the length and orientation in the pertinent maps (Figures 8-10). The same holds true for the 2nd order tributary rivers considered in the current study. River course and structural elements are compared with each other using the same type of mirror rose diagram, excluding the synoptical semicircle (Figure 8e).
Based upon the tectonic planar elements the FW study area was compartmentalized into four morphotectonic units (MTU 1 to 4), all of which reveal a marked NE-SW trend striking parallel to the fold axis of the anti-and synclines, e.g., Berga Anticline and Teuschnitz Syncline (Figures 1 and 8a-e). The channel patterns which are subdivided into 1st (Rodach, Wilde Rodach) and 2nd order streams differ, in places, markedly from the afore-mentioned structural patterns, being endowed with a second maximum displaying a NW-SE trend (Figure 8f-i). While in the Teuschnitz Syncline (MTU 1) the channel pattern accords well with the structural one, in the MTU 2 and MTU 3 the second NW-SE trend stands out irrespective of stream order. The situation becomes more complex upstream where both maxima are present. The closer to the NW-SE striking FLFZ and more downstream the NW trend prevails in the rose diagrams of the drainage pattern. At a more distal position of this highland boundary fault, the Variscan NE-SW arrangement becomes more pronounced. The latter structural pattern gives rise to a special drainage pattern called "antler-type" (Section 5.4). (Figures 1 and 8e). Artificial channels indicated in Figure 8a trenched to feed saw mills and facilitate the work of raftsmen are mentioned here but they are left unaddressed during any talweg measurements. By and large the SSS reflects a pristine channel system.
In the study area of the MGC, the NW-SE trend of the FLFZ and its parallel parasite faults can also be identified in the pattern of the 1st order stream represented by the Schorgast River (Figure 9a,b). The 2nd maximum straddling around ENE-WSW recorded from the 1st order rivers can be traced back to the Variscan fold axis which rotated anticlockwise from the NE-SW towards ENE-WSW. The "antler pattern" visualized by two equivalent maxima running NNE and NE find only representation in the drainage pattern of the 2nd order streams associated with a minor NW-SE maximum.
The south-eastern-most study area in the FG provides a different picture far more conspicuous than its NW counterparts. Its major anticline still stands out as a NE-SW striking domal structure, reaching an elevation of as much as 1052 m with the Schneeberg summit. The "antler pattern" of the river drainage system got almost obliterated ( Figure 10). The Steinach River no longer runs parallel to the 1st order fold axis or hinge line draining down towards SW but runs along its WNW limb as a straight strike stream parallel to the outcrop of the low-grade metamorphic rock series ( Figure 9). The 2nd and 3rd order fold axes run perpendicularly to the tributaries on the right side of the valley.

Degree of Sinuosity and Drainage Classification
Objects and purpose The hydrographic parameter channel sinuosity (=S value) is the principal tool to categorize the fluvial drainage systems, in general, and to draw a more detailed picture of the current SSS, in particular.
After Schumm s (1963) [110] classical studies of the sinuosity of alluvial channels, an progress report of investigations has been given by Church (2002) and Beechie et al. (2006) [111,112]. The S values obtained by the method put forward in Section 2 directly translate into the classification scheme proposed by Charlton (2014) [10]: Straight < 1.1, sinuous 1.1 to 1.5, meandering > 1.5. An amendment has been made based upon the datasets obtained from our previous studies of the proximal foreland rivers enabling us to subject the different streams to a more subtle subdivision and genetic interpretation. This newly introduced classification scheme can be applied to the basin-and-range topography in a sedimentary foreland and the upstream basement, alike (Figures 8-10, Table 2): Straight: <1.1, low-sinuosity: 1.1 to 1.3, high-sinuosity: 1.3 to 1.5, moderately meandering: 1.5 to 2.0, strongly meandering: 2 to 5, very strongly meandering with sedimentary trap sites > 5 [62]. According to this classification scheme the majority of channel systems under examination plot into the pigeonhole low-sinuosity, a channel type mostly found in the TA of the study areas (Table 2). True straight rivers are restricted to the CA and TA of low-grade and very-low-grade metamorphic rocks (Table 2, Figure 11a-c).

Channel Density and Channel Floodplain Ratio
Objects and purpose: The analyses of the channel density and channel-floodplain ratio are markers for the paleo-gradient, the flow rate and can be used as flow confinement indices. They should always be handled side-by-side with each other to find out which scope can be achieved most successfully by this approach (Figure 12).
The contour maps of the channel density of the VLG and LGG-types are rather homogeneous with a distribution of blue anomalies aligned like pearls on a string in a characteristic zig-tag pattern (Figure 12a,b). This regular patterns mirror a downstream trend of increasing sinuosity. The contour map illustrating the channel/flood plain ratio shows two eye-catching anomalies in each contour map (Figure 12c,d). They are remarkable in two different ways. It is the river section where older floodplains in the TA manifested by a tripartite set of intra-basement-paired and unpaired depositional terraces of the 440-500 m-level came into existence and they mark the transition from the TA into the DA (Figures 3b and 8a). Similar features in the LGG-type witness a stacked pattern of depositional terraces in the Fichtelgebirge Foreland of the "Weidenberg Bowl" (Figure 5).

Slope Morphology vs. Bedrock-Lithology and Tectonic
Objects and purpose: The x-y plot of elevation vs. distance illustrates the direct response of the talweg to the bedrock lithology from source to the point of confluence. It is not the first attempt of highlighting this fluvial relationship as shown by Hack (1973) [113] who put forward a stream gradient index. Substituting the slope angle in the x-y plot for the elevation gives an additional insight into the variation of valley morphology as a function of the fluvial facies (CA, TA, DA): wide-and-shallow valley ⇒ wide-angle ⇒ acute-angle V-shaped ⇒ U-shaped valleys (Figure 11a-c). Combining it with the sinuosity allows for a direct comparison of 1-D and 2-D patterns.
Longitudinal profiles of the rivers have been studied by Gomez et al. (1996), Zuchiewicz (1995) and Demoulin (1998) [114][115][116]. The latter study refers to a basement resembling, as far as the lithology and metamorphic regime are concerned, the VLG-type ( Table 1). The VLG and LGG types are described as concave with some knickpoints highlighting 2nd order concave riverbed sectors. The gradient of the LHG is made up of convex subsections within the basement and warrant an overall classification as convex. In the x-y plot illustrating the steepness of the three gradients the LHG type stands out by the steepest inclination ever in the area under study (Figure 11d). All longitudinal profiles start at the watershed and end on the floodplains in front of the basement edge so as enable a direct comparison of the SSS.
There is a third parameter related to the slope angle; it is the linear fold tectonic elements and the dip of slope in the Steinach drainage system ( Figure 13-see also section A to F). The Steinach River runs at an acute angle along the WNW limb of the Fichtelgebirge Anticline. The most steeply dipping slopes in the Steinach drainage system are within the 2nd order channels debouching into the 1st order river from the right-hand side ( Figure 13). The 2nd order channels are governed by NW-SE trending tectonic structures perpendicular to the 1st order fold axis of the Fichtelgebirge Anticline. The maximum of slope steepness is attained where the SW plunging fold axes of the B 1 and B 2 folding stage on the NW limb are perpendicular to the slope of the 2nd order channels.

Hydrodynamics and Grain Parameters (GMS Tool)
Objects and purpose: Reliable data on the most recent hydrodynamic regime can only be obtained by the GMS analysis of particles eroded, transported and deposited within these three SSS which mostly consist of gravel-sized siliciclasts (Table 1). It is an aid to the precise differentiate the landforms (Section 4.1) and determine its proportion to mixed-type depositional environments from land to sea. This is common in the SSS where fluvial, colluvial and cryogenic processes contribute to shaping the valleys.
A tripartite methodological approach has successfully been taken in previous projects investigating coarse-grained glacial, fluvial, marine sediments and different sorts of mass wasting products and named as GMS gravel analysis: Granulometry-morphologysitumetry (orientation of gravelly clasts) [62,84,117]. The gravel is subjected to a visual morphological examination while subdivided into six classes of roundness from very angular to well-rounded according to Illenberger (1991) [60] and discriminated granulometrically using the suggestions put forward by Pettijohn et al. (1987) [118]. The 3rd method is called situmetry which sedimentologists and geomorphologists have been losing sight of during the recent past, excluding applied geoscientists who use clast orientation in glacial till for mineral exploration [56,57,119]. To exemplify the tripartite approach of the GMS tool, an in-depth treatment of the LGG-Type has been conducted (Table 3, Figure 14).
The clast orientations suggest both colluvial and alluvial/fluvial deposits exist in the CA as indicated by the orientation of the 1st and 2nd order maxima and the angle between them ( Table 3). As noted in Table 3 and illustrated in Figure 14, the four major parameters are (1) the orientation parallel ( Figure 14 diagram 11a), or (2) perpendicular to the paleocurrent/slope direction ( Figure 14 diagram 12), (3) the angle between the maxima, and (4) the merging of maxima producing a fan-shape wide-angle sector. A trend to alluvial/fluvial processes is indicated by acute-angle maxima conspicuously well separated from each other, whereas colluvial and deltaic processes tend to form obtuseangle maxima and an amalgamation to fans. The method ought to be run in combination with grain size and morphology of clasts (Table 3).
In the TA, the channel system evolves from an intermediate sediment trap into a channel system which at its end converts into the foreland DA. The minor straight creeks of 2nd order debouch into the trunk river via coarse-grained deltas. The small channels in the CA of the VLG-type resemble its counterparts in the LGG-type and based upon that enable a morpho-stratigraphic correlation of the tributaries (Table 3-Stein 2, 11b, 11c). The LHG-type is a mixture of the two end-member types LGG and VLG.

The Composition of Autochthonous and Allochthonous Mineralization of the SSS
Objects and purpose: There are two types of mineralization. The autochthonous overburden contains the minerals useful for thermodynamical computations of the physicalchemical regime of the chemical weathering. The parameters are the redox conditions (Eh), the acidity/alkalinity of the meteoric fluids (Eh) and the near-ambient temperature regime. Radioactive age dating, e.g., U/Pb, and organic methods such as palynology allow a determination of the paleoclimate ( Table 5). The gravel-size debris are representative of the allochthonous mineral aggregates. Thermodynamic computations and age dating encompass a wider field of geosciences and therefore, the results are discussed in Section 5. This is also true for the composition of the debris (Table 4).

Composition and Age of Formation of the Mineral Association of the Supergene Alteration
The supergene alteration zone is subdivided into two alteration zones, one called zone of argillitization, due to the abundance in phyllosilicates and the other zone of duricrusts/orecretes rife with chemical residues made up mainly of Fe, Mn, Ti and Al [120] ( Figure 15, Table 1). A comparison of the mineral associations reveals great similarities between the VLG-and LHG types ( Figure 15). Both types only differ from each other by the absence of orecretes in the LHG-type (Mn orecretes, uraniferous orecretes and aluminumphosphate-sulfate minerals (APS)). The LGG-type is underlain by crystalline bedrocks abundant in alkaline feldspar and primary U-bearing host minerals rendering its bedrocks highly susceptible to chemical weathering. It is abundant in newly formed smectite, mainly nontronite, kaolinite and a varied spectrum of uranyl phosphates (torbernite, saléeite) and hydrosilicates (uranophane), all of which are amenable to radiometric age dating as it is the case with the Fe-Mn counterparts from the VLG-type ( Figure 15). Both sorts of duricrusts reflect paleo-aquifers that had direct implications as redox indicators and geo-acidometer on the chemical weathering of the bedrock in the CA and the incision of the drainage systems in the TA (Figures 11a,c and 15). Uraniferous phoscretes and silcretes and manganiferous ferricretes are two sides of the same coin, representatives of the continental aquatic regime. In one deposit of cryptomelane, a direct comparison has been conducted and yielded an age of formation of 4.20 ± 0.33 Ma for cryptomelane while torbernite gave an age of 4.55 ± 0.02 Ma. Thereby they furnish convincing evidence to belong to the same Pliocene groundwater system and cast in to the role of environment markers [48,49].

Thickness and Variation of Autochthonous Mineralization
In the LHG-type, reworked regolith is second to none compared to the adjacent study areas in terms of aerial distribution and thickness (Figures 8a,c, 9a,c and 10a,b, Table 1). Its loamy weathering mantle attains as much as 27 m on gneissic bedrock, with mean thicknesses around 11 m and minimum thicknesses of around 1.5 m on steeper slopes while it goes down to zero on bluffs [104,121]. Along the talweg of the TA the thickness diminishes to almost zero in the non-alluvial sectors of the SSS. By contrast in the DA the thickness increases again to values to more than 6 m [82]. Despite of this considerable thickness no regolith profiles have been well preserved due to its cryogenic reworking. Irrespective of the age of formation and its redeposition, the fine-grained weathering mantle is crucial as to the degree of sinuosity in the CA (Section 5.3.2). Over time as a channel incises through the regolith the sinuosity increases until the channel hits the rigid basement rocks underneath. There the channels had to succumb again to the guidance of the straight tectonic elements and the rigid lithology of the basement (Figures 8-10).
The size of deposits produced by mass wasting in the various sites well accord with the size of monadnocks and block fields and testifies the genetic relation between the two landform types ( Table 1). The quantification of fluvial and colluvial deposits is given by area percentage in Table 1 and displayed in the maps of Figures 8-10. Measuring the thickness of the overburden on top of the entire basement cannot be achieved without geophysical surveys across the entire study area supported by shallow drill holes used for reference and calibration. Geoelectrical methods have been applied selectively to get an insight into the thickness and the composition of the colluvial and fluvial deposits throughout the assessment of mass wasting [61]. The apparent resistivity values are as follows and arranged in order of increasing depth for comparison: Fluvial bedload deposits: 120-150 Ωm, solifluction creep: approx. 95 Ωm, mud and debris flow: 50-60 Ωm, regolith: 90-130 Ωm, basement: 340-750 Ωm. Fluvial deposits, mass wasting products, regolith and bedrocks can precisely be defined and the thickness determined. Therefore, based on the mapping results in Figures 8-10 and some thickness data from outcrops, geophysical logs and shallow drill holes the following zonation can be reported: CA very high thickness (max. 27 m), TA: considerably reduced thickness (max. 1.5 m), DA thickness increases (max. > 6 m).

The Physical-Chemical Conditions and the Paleo-Climatic Regime
The origin of the variable landforms is the result of endogenous (structural landforms) as well as the exogenous forces (sculptural landforms). The landforms developed at the interface between lithosphere, atmosphere, and hydrosphere leading to a special landscape, in the current study to these three SSS in a basement landscape (Table 1, Figure 16).
Weathering/supergene alteration mediating at the interface of the different spheres where the SSS formed [122][123][124][125][126][127][128]. It is a series of physical-chemical processes that is governed by the climate during the Neogene and Quaternary established on a Paleozoic and Mesozoic-Paleogene platform ( Table 5, Figures 2b and 16a,b). The thermo-dynamic calculations to constrain this physical-chemical regime have been extensively discussed in Dill et al. (2010a, b) and in Dill and Wemmer (2012) [48][49][50] and the results of which are portrayed for the region in Figure 16c-f. The translation of these physico-chemical data into the Cenozoic paleoclimate is discussed in context with the age data obtained from uranyl silicates and phosphates as well as K-Mn-Fe oxide hydroxides, unless stated otherwise (Figure 2b). The latter chemical systems U-Si-Mg-Ca-Cu-P and K-Mn-Fe are the only ones which can successfully be applied to the Cenozoic land-forming processes because both parameters, the "temporal" and "environmental" ones are inherent to each of the afore-mentioned mineral series.
Prestage of the SSS: The prestage of these SSS covers the Early to Middle Miocene (23 to 14 Ma) where lignite seams came into existence immediately SE of the study area of the LGG-Type in the Eger Rift Graben [51]. Forest peat swamps developed under oscillating water levels in a fluvial-lacustrine environment on the Neogene regolith. More basinward black shales yielded a palaeoflora of early Miocene age [129].  (Figure 16a,c,d)  Proto-SSS: The Pliocene to Early Pleistocene (4.0-1.8 Ma) evolution of the fluvial drainage systems is based on uranyl phosphates and -hydrosilicates/uranophane [Ca[UO 2 |SiO 3 OH] 2 ·5H 2 O] which evolved under a subtropical climate but slightly drier conditions as far as the Pliocene part of the fluvial landscape is concerned when silica-cemented gravelly deposits developed (Figure 16b,e,f) [131,132].
Modern SSS: The Early to Middle Pleistocene age (Early Pleistocene < 1.8 Ma) has been constrained based on the Mg uranyl phosphate saléeite [Mg[UO 2 |PO 4 ] 2 ·10-12H 2 O] referring to a period of time which is synonymous in Central Europe with the early phases of the glacial period. In the study area it is only represented by periglacial landforms within the CA and DA of the SSS such as nivation cirques and pockets undergoing cryoturbation (Tables 4 and 5) [92,93,97] (Figure 7g).
The youngest floodplains in the DA of the SSS are Holocene in age, the humid midlatitude temperate climate of which has been investigated by (Zöller et al., 2007;Dill et al., 2020a) [62,68].

Maintenance of a Metastable Hydrodynamic State Needs the Constraints of a Structural-Lithological Framework and Effective Drivers
The implementation of a SSS in nature needs a framework-driver couple. Both are linked with each other like "hardware" and "software".
The framework: Due to the metastable character of the SSS, its channels always tend to convert into a drainage system of higher sinuosity and try to achieve a higher level of hydrodynamic stability [11,111,133] (Table 2, Figure 6e). To maintain a SSS over a longer period of time and distance the course of the channels must be coerced into a rigid framework. For the case of the three SSS under study, such a setting is well established in basement terrains composed of sturdy bedrock lithologies and subjected to a tectonic overprinting that led to a perfectly well-structured fold and fault pattern with a pronounced morphotectonic outward appearance (Tables 1 and 5, Figures 3-5). A Paleozoic through Paleogene framework is mandatory to create a SSS (Figure 16a,b).
The drivers: Additionally, an SSS terrain needs a driver to assure a constant high gradient and to ensure continual erosion, transport, and deposition by fluvial and mass wasting processes ( Figure 11). In the present study areas, this driver is guaranteed by syngenetic tectonic pulses during the Cenozoic and fostered by contrasting climate regimes controlling the state of the agent water and impacting on the supergene alteration (Figures 2b and 15, Tables 4 and 5). The current study of the sinuosity over the full length of the drainage system illustrates that its development is not unidirectional in nature.
Stable drainage systems with a fairly high sinuosity degree of S > 2.5 occurred at the beginning but they developed into a metastable type with sinuosity degrees approximating S = 1; farther downstream they end up in the foreland in an ultrastable drainage systems as the S values reach values of as much as 11 at very low gradient (Tables 2 and 5) [62]. Such abnormally high values are indicative of the marginal facies of a fluvial-lacustrine environment. In the fluvial regime upstream of the entrance into this sediment trap sinuosity S values in the range 2 to 5 are common ( Table 2). They may episodically turn into more stable drainage systems within intermediate sediment traps of the SSS. This bimodal sinuosity distribution has been elaborated for the three facies zones CA, TA, and DA of the SSS-see Section 4.4.2. It is a consequence of climate change and caliber of sediment load delivered by fluvial and mass wasting processes.

The Paleozoic Tectonic Framework
Each of the three study areas has a south-westerly plunging 1st order anticline, the backbone of the framework (Figure 1). Their outward expression when exposed to surface; however, strongly differs as a function of the metamorphic-dynamic evolution during the Variscan orogeny and so does their impact several hundred million years later on the morphotectonic expression. The reference sites are arranged in order of increasing grade of dynamo-metamorphism, the lithological-structural features of which provoke a different morphotectonic expression in the landscape of the basement and thereby endow every type of SSS with peculiar features (Tables 1 and 5).
The deformation in the study area of the VLG-type took place between 360 Ma and 320 Ma. The orientation of the SSS can still be traced back to this folding because of the very low grade of metamorphic overprinting which left the structural features relevant for the SSS almost unaffected (Figure 8). The planar tectonic elements governed the river course of the Rodach drainage system, guided the progradation of the terraces into the VLG basement and the subsidence of the U-shaped valleys in the lower reaches of the Rodach drainage system (Figures 1, 3b, 6f,g and 8a). Upstream of the Wilde Rodach the stepwise piedmont staircase system evolved in accordance with the axial-plane parallel fault zones, and the litho-variance of the Paleozoic basement rocks (Figures 3a and 8b-e).
These NE-SW striking Variscan structural elements also contribute to the several riffles, cascades and rough bedforms particularly when chert forms the bedrock in low-sinuosity channels (Figures 6a and 8c). They are the sites where straight channels deviate into sinuous ones as the pebbles and boulders tend to increase in the sediment load (Figure 8c). Slight changes, be it an increase or decrease in the input of gravelly sediments, give rise to moderately meandering streams via wandering channels (Figures 6e and 8c-sinuosity:  1.125 ⇒ 1.580).
The NE-SW-, draining rivers display a ramification towards NE into two, rarely three, rivulets like the paired or triple-pointed ends of the antler of a roebuck ("antler drainage pattern"). This newly introduced channel pattern is intermediate between the dentritic and parallel ones [10]. The 1st order rivers, such as the Rodach River mark the 1st order fold axis of the large anticlines such as the Berga Anticline which plunges towards the SW. The 2nd and 3rd order axial plane parallel folds homoclinal with the major one are represented by the ramification of the channel network observed in the NE-SW and SW-NE draining 2nd order streams. Being homoclinal with the 1st order anticline provokes the "antler" open towards the NE. If the plunge of the 2nd and 3rd order folds reverses, the "antler" will be open towards the SW (Figure 8).
The "antler" pattern is vaguely expressed in the 2nd order streams of the MGC (Figure 8), because the Central Anticline in the MGC is rather poorly represented and only recognized in the lower reaches of the 1st order pattern. The middle-to high-temperature regional metamorphism between 430 and 360 Ma largely eradicated the fold structures [134].
The asymmetric development of the Steinach River, occurring as a consequence of an anti-clockwise rotation from the hinge-line parallel Proto-Steinach to the current position, did not allow for such antler patterns to survive undistorted. The rotation point is situated where the modern Steinach leaves the basement to debouch its bed load into the "Weidenberg Bowl" where it triggered a set of at least three depositional terraces arranged in a staircase-like manner-see examples from the South Swedish Dome [135] ( Figures 5, 10 and 13).
The low-pressure regional metamorphism in the FG began around 330 Ma and was immediately succeeded by the intrusion of granitic magmas which ceased around 290 Ma. Approximately 270 million years later, the Steinach fluvial system came into existing when doming along the FG Anticline forced the river to incise along the NE-SW hinge line from its culmination towards SW (1st order fold axis) and ended up on the WNW limp. As a result of this asymmetry the fluvial undercutting gradually shifted from the trunk river to the right-hand tributaries (2nd order stream) the steepest slopes at an almost right angle to the 2nd and 3rd order Variscan fold axes ( Figure 13).
During the Mesozoic uplift processes have been recognized in the area and proved by apatite-fission track studies [136]. They are, however, without any relevance for the evolution of the SSS. No proven chronologically tested evidence of landforms has been preserved.

The Cenozoic Tectonic Driver
The effect of these tectonic drivers impacting on the SSS need to be split up into three classes. They are destructive or non-destructive with regard to the bedrock of the SSS, they follow a tectonic "blue print" newly established during the Cenozoic or rejuvenate ancient structures of the framework, and can be named as 1st order or 2nd order knickpoint makers. It has to be noted that some methodological restrictions are imposed on the timing of the neotectonic processes. The interval between 0. 4 Ma (lower limit of U/Pb dating) and the time covered by radiocarbon data provided by Zöller et al. (2007) and Dill et al. (2020a) [62,68] is not covered by absolute age data but based only on lithostratigraphic correlation of terraces and a more precise chronological assignment cannot be given here.
Frameworkvs. bedrock-related tectonic processes: The NE-SW striking structural elements are prevalently caused by rejuvenation of the pre-existing Variscan-type near and planar tectonic elements (see Section 5.2.1). The NW-SE striking structural elements are provoked by neotectonic processes. Its most conspicuous representative is the FLFZ (Figure 1).
Non-destructive vs. destructive processes: Non-destructive processes reactivated the Variscan structures following the same NE-SW trend as their Paleozoic predecessors and normally resulting elongated domelike structures. The role model is the "Fichtelgebirge Anticline", the host structure of the LGG-type drainage pattern the onset of which can be traced back into the Mid-to Late Mesozoic (see Section 5.2.1). Similar to the SSS under study in Germany, domelike uplifts in a passive-margin domain strongly impacted in Sweden on the development of Neogene drainage systems such as the Eridano River of the Baltic Basin during Late Jurassic or Early Cretaceous also [135,137]. A newly formed non-destructive element, a swell running NW-SE terminates the study area towards the NE and is herein referred to as "Franconian Traverse Swell" (=FTS) (Figure 1). Its hinge line act as the watershed between the Rivers Rhein and Elbe/Labe and two subhorizontal paleosurfaces tilted away from it. The western paleosurface is under study in this project and homoclinal with the SSS. The destructive counterpart to it is the FLFZ bounding the basement towards the SW against the foreland (Figure 1). These planar architectural elements constitute the Cenozoic or neotectonic structural framework for the study area and come into effect in the SSS on a regional scale as 1st order knickpoints and on fluvial scale as 2nd order knickpoints (Figures 11 and 12).
1st and 2nd order knickpoints: When it comes to decipher the evolution of the Cenozoic tectonic processes within the lower reaches of the SSS, three parameters are important; the absolute height of vertical displacement, the type of displacement (e.g., non-destructive), and the positioning of the fault plane relative to the basement-foreland topography and the 1st order knickpoint of the SSS. The NW-SE trending FLFZ measures different heights of vertical displacement from NW towards the SE. The actual difference between the lowest erosional terrace on the fault s basement shoulder and the fault s lowest foreland floodplain measures in the VLG 140 m, in the LHG 200 m and the LGG types 150 m (Figures 3-5). The lowest erosional terraces are genetically different from the Quaternary depositional terraces situated at a height of 10 to 25 m above the present floodplain.
Faults, which are the leverage to rise the gradient of the SSS, reach its maximum in gradient and vertical displacement in the LHG-type (Figures 4, 9, 11d and 17). The basement-foreland topography coincides in this site with the FLFZ (Figure 9). In the LHGtype a sharp contact exists between the TA and DA, accentuated by the steepest escarpment ever observed in the region. The original river course started to follow the Variscan plan of construction being oriented NE-SW, while ongoing deepening during the Cenozoic redirect its course into the NW-SE direction of the FLFZ (Figure 9a,b). In the LHG-type, the NW-SE trend of the FLFZ can be identified in the pattern of the 1st order streams where it is represented by the upper reaches of the Schorgast River (Figures 9a,b and 11d). The geomorphological evolution of these SSS coincides with the tectonic evolution and the SSS development is held to be synkinematic-destructive (Figure 17b). Vertical movements along the FLFZ are also responsible for the development and preservation of paired and unpaired depositional terraces (Figures 8a, 10b and 12c,d). In the VLG-type area, the knickpoint of the gradient as well as the boundary between the DA-TA provoking the widening of the V-shaped valley to convert into a U-shaped one is located inside the basement (Figures 8a and 12c). In the VLG-type the DA invades along U-shaped valleys the basement while the FLFZ is located approx. 3 km to the W of the knickpoint which is denoted by a hydrographic anomaly (Figure 12c). The U-shaped valley is the outward expression of a tectonically abandoned part of the SSS and the system, consequently categorized as postkinematic-destructive (Figure 17a).
In the LGG-type this relation is reversed, the foreland sediments unconformably rest upon the basement forming a non-destructive flexure bulge. The FLFZ cuts through the Mesozoic foreland sediments in the "Weidenberg Bowl" which is largely occupied by fluvial terraces of Pleistocene to Holocene age located outside the basement, as it is with the 1st knickpoint (Figures 10b and 12d, Table 4).
The staircase-like bedsets of terraces known from the "Weidenberg Bowl" can be taken as reference for the fluvial incision and deposition adjacent to the FLFZ along the basement. These terraces evolved from ancient fluvial terraces and fluvial to colluvial pediment-glacis deposits of Neogene age laid down in the foreland of the basement by torrential rivers debouching into the Mesozoic Foreland Basin. They were subsequently almost completely "eaten away" during the Quaternary by the ongoing fluvial incision. The earlier terraces developed following a cycle of deposition and entrenchment in the formerly deposited sediments or as a secondary "inner" valley by entrenchment into older deposits provoked by morpho-tectonic processes under climatic regime strongly contrasting with the more recent ones (see climatic evolution this study). During the Pleistocene and Holocene a constant but stepwise adaptation to the changing climate shifted more and more to a forward carving out of the valley prior to the fluvial deposition. The geomorphological system is categorized as prekinematic-non-destructive (Figure 17c).
The knickpoint zones in the VLG-and LGG reference types are not only marked by a striking hydrographic anomaly but also accentuated by an "elbow" in the channel system (Figure 12c,d). These anomalies furnish evidence for a "bottleneck" attesting to a rising flow rate in the river course (highly confined flow) before relaxation into the flat-floor-shaped valley (unconfined flow) inside the basement or in the adjacent foreland according to the hydrodynamic Bernoulli Law. The 1st order tectonic knickpoints separate the DA from the TA of the SSS and help to constrain the relation between shaping the landscape by the SSS and the tectonic events. In connection with the terraces, each of which are assigned a post Middle Pleistocene age, these morpho-tectonic processes in the SSS are of younger than 0.4 Ma.
Parasite faults running parallel to the FLFZ within the basement created knickpoints in the talweg sparking different fluvial-colluvial environments to develop in the upper reaches of the SSS (Figures 9, 10 and 12c,d). They control the internal structuring of the TA and are denominated as 2nd order knickpoints. The latter are accountable for intermediate sediment traps within the LGG-and VLG types but not in the multiple-convex profile of the LHG-type which misses a tectonic pulse during this period of time. As a consequence of uplift in the VLG-and LGG-type the incision of both reference rivers could not keep pace and their upper reaches were dammed giving rise to the intermediate sediment traps while downstream a knickpoint developed in the river run (Figure 11a,c). In the LHG type, a series of staircase-like convex increments evolved instead (Figures 9b and 11b). In the LHG-type a convex knickpoint exists between 465 and 500 m a.m.s.l. marking the transition from a poorly developed CA into the TA accompanied by an erosional unconformity that "has eaten away" the wide-angle V-shaped valleys present in the neighboring reference study area (Table 5, Figure 11a-c).

The Lithological Impact on the SSS
The interaction of the bedrock lithology with the hydrographic (Figures 3-7, 11 and 15) and hydrodynamic processes operative in these SSS (Figures 1 and 14, Tables 1, 3a,b and 4) are a consequence of the mediating climatic regime leading to different cross-sectional morphologies of the valley and longitudinal profiles of the talweg (Figure 16).

Bedrock Geology-Strongly Altered by the Climate
The uplifted bedrock is truncated by a subhorizontal paleosurface tilted away from the watershed, with a relic landscape carved out of it-see Section 5.2. The paleosurface is interpreted as a peneplain or etch surfaces, the formation of which is still subject to controversial debate on a regional scale [76,77] and international level [138] (1) peneplanation (Davis' theory), (2) pediplanation (Penck' s theory), (3) etchplanation (Wayland's theory). What is beyond any debate is the strike of the paleochannels and the ensuing modern fluvial equivalents of the prestage and proto channels of the SSS all of which show a NE-SW trend and a drainage, homoclinal with the dip of relic landscape/peneplain ( Figure 16, Table 5). A great deal of the drainage channels is oriented along the Variscan anti-and synclines as to strike and dip which experienced a Neogene reactivation as the Eger Graben (Oherský příkop) subsided [139,140]. The latter tectonic process has already been heralded during the Late Cretaceous to Paleocene (68 to 59 Ma) ( Table 5) [141]. This period of time was constrained geochronologically at the interval from 29 to 19 Ma by Huckenholz and Kunzmann (1993) [142] and marks the period when the wide-and-shallow drainage system with a sinuosity > 2.5 and coarse-grained channel fill came into existence ( Table 5). The widespread occurrence of kaolinite-group clay minerals supports the idea of a humid tropical climate during this period which is also attested to by the presence of Fe phosphates. The variable contents of bivalent and trivalent Fe accommodated into the lattice of phosphates point to oscillating redox conditions and acidic meteoric fluids (Figure 16c). The wide-and-shallow valley of the VLG-type are filled with siliceous gravel and cemented by Fe-Mn orecretes under an oxidizing redox regime (Figures 6a, 11a and 16a,d). Gravel beds with encrustations of Fe-Mn oxides occur also in the Eger Rift Graben similar to those found in the wide valleys at Stegenwaldhaus. The lateral facies change in a fluviallacustrine regime reveals a trend to more allitic-sialitic weathering conditions [143]. This is also corroborated by the silicate karst observed in ultrabasic igneous rocks of the LHG -Type [144] (Figure 7c).
The ancient wide-and-shallow valleys are also the protore-collectors [subeconomic mineralized matter which was concentrated at a lower tenor than the subsequent ore that has been derived from this prestage mineral assemblage] for the V-shaped valleys which host the true ore deposits [145,146]. The genetic link between the Miocene predecessor drainage system and the VLG-type river is proved by the accumulation of cassiterite and the formation of Sn placers along the river run of the Wilde Rodach (Figure 11a).
During the waning prestage another tectonic phase took effect along the SW boundary of the basement under consideration in the South German Block [102] (Figure 1, Table 5). Uplift along the NW-SE running FLFZ are crucial for the height difference and the different fluvial gradients (Figures 1, 3-5 and 11). Tilting, uplift and disintegration of the South German Block that is geomorphologically synonymous with the foreland of the South German Basin, kickstarted the development of the Proto-SSS during the Plio-Pleistocene with the creation of the headwaters of the SSS [102] (Table 5, Figure 16b).

Bedrock Geology-Unaltered to Slightly Altered by the Climate
The headwaters (CA): The onset of the CA is marked by a steepening of the slope angle along with a diminishing degree of sinuosity (Figure 11a-d, Table 5). The drainage sector resembles the Modern SSS as far as the sinuosity is concerned but still possesses all hallmarks of the relic Neogene fluvial system as to the slope angle and the valley fill (Figure 11a-c) (Section 5.3.1). In the CA of the LGG-and LHG-types the slope angles are far below 20 • . Palisades and tors arise from hilltops the subhorizontal surface of which has an inclination of between 5 • to 10 • . Due to the joint action of fluvial and mass wasting (slide, creep) erosion at slope angles greater than 10 • these landforms with almost subvertical walls tend to disintegrate into felsenmeers, blockmeers, blockfields and block streams (Figures 7b,c,e and 16b, Table 5). The final result is a skeletal highly permeable landform series surrounded by deposits completely depleted as to their argillaceous matter and undergoing reworking and sorting during cold climates and periglacial conditions [86,[147][148][149][150][151][152][153][154][155].
Landforms and slope angles correlate with the F/M ratios (ratio of fluvial to mass wasting products) which are below 1 in contrast to the VLG-type where the ratio exceeds 2 ( Table 5). The less permeable lithologies of the slate mountains of the VLG exceed the LHG and LGG types with regard to the aerial extension of floodplain deposits; they are characterized by the highest drainage density, the highest degree of sinuosity and a high rate of slate pebbles and gravel in the sediment load. Their smooth hilltops are blanketed all across with topsoil and densely forested with softwood trees. This is also corroborated by the mass wasting products which rank lowest among the reference types ( Table 1). The slates that were subjected to very-low-grade regional metamorphism still feature the low permeability of their argillaceous parent material yet do no longer show the weathering behavior of their predecessors, grey and black shales [95,[156][157][158].
The LHG-and LGG-types preserved most conspicuously the Neogene supergene alteration with a thick regolith that was reworked under colder climatic condition into the CA by mass wasting processes, leaving as residues the different monadnocks which became instable on gentle inclined slopes. It has to be noted, that etch planation is held to be the initial process accountable for the built-up of the afore-mentioned landforms on hilltops but subsequent cryoplanation cannot be sidelined to have contributed to the formation and preservation of these rocky hilltops. On the other hand, the sinuosity increases relative to that of the VLG, where the fluvial deposits are more widespread around the hilltop devoid of monadnocks. This triple of numerical parameters sinuosity, slope angle and the F/M ratio is a diagnostic dataset to assess the supergene alteration or prestage impact on the development of the CA of the SSS.
The transport area (TA): The Steinach River is running subparallel to the strike of the metamorphic host rocks (Figures 10a and 13). Sectors of the drainage system framed by epigneisses of high rock strength tend to form fluvial deposits in a narrow gorge with a steepening of the slope angle accompanied by outcrops of the bare rocks and a decrease in the sinuosity between 1.023 and 1.095. In those sectors where phylliticquartzitic metasediments constitute the country rocks colluvial sediments prevail over fluvial deposits and the sinuosity increases to 1.085 to 1.116 while the valley widens, and the slope angle diminishes (Figures 10 and 13). Where the fresh bedrocks are largely exposed in the valley, irrespective of the varied lithology the slope angles are quite similar (Table 1) However, there are also marked differences in the lack of the V-shaped wide-angle valleys in the LHG-Type. The VLG-Type shows the steepest slopes, exceptionally high F/M ratios and the greatest variety in the cross-sectional valley morphology from wide-angle V-shaped valleys to flat-floor-shaped/U-shaped valleys. This fact closely correlated with elevated rock strength and low vulnerability of the bedrock against chemical weathering. The TA was shaped under oscillating cold and warmer phases of the glacial period ( Table 5).
The depositional area (DA): Only in the DA of the VLG-Type, with its intra-basement terraces and floodplains, occur exceptionally steep slopes in a fluvial landform denominated as flat-floor-shaped U-shaped valley (basement) (15 • -32.1 • -37 • ) (Figures 8a and 12a). Its rectangular cross-sectional shape originated from quartzitic bedrocks and chert of the high rock strength ever reported from this study area. (Figures 8a and 9c). The sinuosity of the SSS conspicuously rises when leaving the basement whereas the slope angles decrease (Tables 1 and 5). The size of the large plains of the DA in front of the basement is a function of the gradient attaining a maximum around the fan delta of the Schorgast-Coserbach Drainage System (LHG-Type) (Figure 11b,d). It is the only site where terraces were completely eradicated leaving behind some butte-like gravel piles (Figure 9c). On the hillwash plain boulders and gravel do not form part of the alluvial-fluvial flow processes. They are either supplied to the channels and rills (confined drainage system) by talus creep or respond only passively in the apron at the foot of triangular slopes as they got deprived of their fine-grained matrix by hillwash [84]. The stacked pattern of terraces within the DA formed during the glacial period between the Middle Pleistocene and the Holocene [81,108,159]. Štěpančíková et al. (2008) [160] recorded a similar set of terraces with the oldest one forming during the Elsterian 2 (0.40-0.46 Ma) from the East Sudeten Mountains. Younger terraces aged Saalian 1 (240-280 ka), Saalian 2 (130-180 ka) and Weichselian (10-80 ka), respectively, might be correlated with the three younger floodplains in front of the FLFZ. Fresh to slightly altered (limonitic rims) rock fragments are selected during gravel analysis using the GMS tool (Table 3). This sedimentological approach taken in different depositional environments from Antarctica to the Arctic is an additional means to decipher the mode of slope and valley formation in the statu nascendi at the beaches, in fluvial channels, mass wasting lobes and wedges and the periglacial (para)autochthonous landforms [62,84,117,161]. Situmetric measurements of the angle between the 1st and 2nd maximum in bimodal distributions of the clast orientation widen as the mass wasting on the valley fill increases relative to the fluvial deposition in the river sector (Table 3).
Cryogenic influence can be figured out using the situmetric reference plain in vertical position as exemplified by the sedimentology and sedimentary petrography of fluvial and cryogenic gravel deposits of the terrace T 1 of the LGG-Type in the "Weidenberg Bowl" (Figures 5a,b and 10b). Fluvial gravel deposits override randomly distributed lithoclasts in a cryoturbated cavity which has been dragged downstream. The vertical situgram displays fluvial gravel beds mixed with colluvial sediments underneath. Both maxima are aligned at acute angle to the floodplain. The vertical situgram of the cryoturbated sediments trapped in a pocket and fluvially reworked still shows a subvertical maximum between 40 to 80 • . The subvertical positioning of gravel underneath terraces and in fluvial sediments is indicative of periglacial landforms such as ice wedges. There are fluvial-(peri)glacial gravelly deposits in the DA of VLG and of the LGG types SSS in which based upon the situmetric measurements the relative influence of fluvial, cryogenic and mass wasting processes can be assessed (see Figure 8a, Table 1).
One can assess and describe the impact of cryogenic processes in a model as follows. Narrow V-shaped valleys carved out of the crystalline bedrock by ideal straight channels as at LGG-type are depleted in bed load deposits (Bernoulli Effect). There exists an intermittent but unimpeded flow regime at rather high speed in the straight channels so that gravel plains can easily be piled up and triggered by outburst and flash floods upstream in the nival zone. Conclusively, a low degree of sinuosity and a high gradient (see LHG-type) are both instrumental for the emplacement of such vast plains at the downstream end of the river run. Albeit no glaciers are covering the mountain tops, periglacially-related landslide triggered outburst floods, attested to by the 1st and 2nd order valleys clogged with mass wasting deposits work in the same way as even better due to the higher frequency [162,163]. Such indirect catastrophic lake outburst floods need no true glacier on the retreat but are typical proglacial lakes and especially destructive as the straight channel system does not slow down movement of the water mass but speed the process up [164].
The passage from the CA into the TA and from the TA into the DA is marked by an abrupt grain size increase caused by changes in the gradient and the presence of knickpoints. It is accompanied by a change of grain morphology as a function of facies with the roundness max fluvial is significantly greater than the roundness max colluvial . Considering the gravel size interval as a single entity, the gravel analyses of the SSS shows a rather homogeneous picture making up around 70% of the total grain size of the fluvial sediments (Table 1). A more subtle analysis discriminating the particle size into boulder, cobble, pebble and sand yields a characteristic split of the grain size range in the VLG-type, with the finer-grained lithoclasts to prevail over the coarser-grained ones: lithoclasts pebble > lithoclasts cobble+boulder (Table 1). This behavior of the clast community correlates with the metamorphic grade of the bedrock. Even though having been derived from the same pelitic parent material (mudstone, shales) they show up as different metapelites and with different ways of disintegration on erosion and transport in the SSS: VLG: slate (and flagstone), LGG: phyllite, LHG: paragneiss ( Table 1). The paragneisses of the LHG are most susceptible to chemical weathering during the Neogene and therefore the bedload got reduced at the expense of clay and silt (Section 5.4.2).

Lithology and Suspended Load
The mineralization in the regolith is controlled by the basement lithology and by the climate as exemplified in Figure 2b. The climate has a strong influence on the physical-chemical regime of aquatic system provoking a drop of the pH far below pH 7 and led to oscillating redox conditions around 0 (Section 5.1). It affected the Paleozoic bedrock (Section 5.3) to a different degree resultant in an argillaceous regolith of particular widespread occurrence in LHG (Section 5.3.1) (Tables 1 and 5, Figure 2b). The clay mineral assemblages in the reference types closely resemble each other. Some striking difference occur in ultrabasic rocks which create serpentine and talc when undergoing chemical weathering in the LHG-type bedrocks as a result of the high Mg content in the parent material. Smectite and nontronite are neither accounted for by a different climate nor by any pyroclastic input but caused by the preponderance of feldspar in the source rock such as the granites of the LGG type. Even at a high uplift rate a constant supply with feldspar is guaranteed in the LGG-type so as to maintain a considerable high pH in the saprock. Expandable clay minerals are very common within the proximal alteration zone of feldspar-enriched basic magmatic rocks [165,166]. The end product of this supergene argillitization/regolith is kaolinite, a marker mineral for acidic meteoric fluids in the LHG-and the LGG-types ( Figure 15).
In the VLG-type kaolinite co-exists with illite-muscovite, chlorite and vermiculite as a primary phyllosilicate in the slates. The primary metamorphic mineral assemblage is in equilibrium with the secondary supergene neomorphic mineral assemblage. Flakes of unaltered slate only need to get chipped off mechanically and transported by the SSS in the pebble grain-size fraction of the bedload ( Table 1). The same holds true also for the metapelites as far as the chloritic and micaceous phyllosilicates are concerned which are reworked as chlorite-mica-quartz phyllite gravel ( Table 1). The situation completely differs from that of the LHG-type where the high-grade metamorphic rocks are strongly susceptible to the pervasive chemical weathering and the adjustment of the clay mineral association by supergene alteration brought about loamy substrate and a complete eradication of the primary structures (Table 1, Figure 11b). In the CA, large areas are covered by sediments derived from argillaceous fluvial and mass wasting processes attaining in places a regolith as much as 23 m thick and unaltered gneissic gravel are the exception than the rule. Therefore, fluvial gravel in the TA are rare and a turbid slurry comes into existence which at the very end spreads across the hillwash plain in the DA (Figures 2b, 4, 9 and 11b,d). Unlike in many regolith blankets elsewhere, the results of supergene alteration in these basement regions are no longer in situ but have to be denominated as parautochthonous due to the reworking by periglacial land-forming processes [167][168][169][170][171].

Lithology and Bedload
Besides clay and silt, gravel-sized clasts form the second maximum within the bimodal grain-size distribution of the fluvial and colluvial sediments (Tables 1 and 4). This bimodal grain-size distribution is conspicuously expressed numerically by the three reference types all of which are undernourished with regard to sand relative to clay+silt and gravel ( Table 1). The coarse-grained fraction is especially widespread in the clast community of the VLG-and LGG types (Table 4). This is especially true for the FG Anticline which forms one of the most significant watersheds between the drainage systems of the Rhein and Donau Rivers (Figure 11c). Considering the DA of the Steinach River, a conspicuous impoverishment of the gravel assemblage in granite clasts is evident. They are very rare constituents in terrace No 2 and No1-Marginal Granite G2 and Core Granite G 3 (Table 4). At the transition CA-TA, the fresh granite exposed in the non-alluvial channels, abruptly converts via large blocks into grus at the edge of the granite. It disappears from the clast size amenable for differentiation with the unarmed eye in the field into a phyllosilicate spectrum sampled at the same site of observation made up of Illite/muscovite < 70%, vermiculite < 30%, kaolinite < 10% (Figure 11c). These mineralogical/ lithological observations agree with those made by Migoń and Lidmar-Bergström (2002) [172] who recorded increased amounts of 7-Å phyllosilicates from the Mesozoic and Paleogene saprolite but not from the Neogene-Quaternary interval. The phyllosilicate assemblage recorded above from the SSS does not only corroborate the attribution in time to the Quaternary but also delineates the metamorphic host rocks of the granite as the major provenance for the bedload. Creeks flowing into the Steinach originate from rills upstream, through confined debris flows grade into fan delta deposits that do not differ from the Pleistocene fluvial deposits of the master river with regard to the gravel assemblage ( Figure 14, Stein 11b, 11c) [173][174][175][176][177]. The bedload and suspended load representative of a bimodal distribution is in terms of size and composition typical of the SSS are another combination mirroring the bimodality during evolution of the drainage system (Table 5, Figure 2b).

The SSS in View of Geodynamics and Crustal Maturity
The SSS reflects a metastable hydrodynamic system which whenever it is possible tries to approximate the stable meandering final state. To force the channel to maintain its straight to low-sinuosity pattern needs a stringent control by the geology and petrography of the country rocks over a significant period of time which measures millions rather than thousands of years. The link between the endogenous processes litho-genesis and tectonics and the hydrodynamic evolution of the drainage systems is provided by the supergene alteration/ weathering [122,[124][125][126][127][128]178,179]. The joint action of the aforementioned processes related to the maturity of the geodynamic units, in the Central European Variscides, which can be subdivided on a global basis into three categories: (1) supermature, (2) mature, (3) pre-mature [120] (Table 6, Figure 1c). Due to the outstanding geodynamic position of the Central European Variscides (Figure 18), in general, and the Bohemian Massif, forming the core zone of Mesoeurope, in particular, the current study area can be taken as the reference type of the mature state of an orogenic crustal section. The supermature stage is today well exposed in the lowlands and peneplains of the deeply eroded Precambrian cratons. It is likewise different as to the landforms and earth surface processes from the mature state as it is the immature state represented, by modern fold belts, both of which need a separate treatment with regard to these drainage regimes and it would go far beyond the present topic to discuss the interlocking of them in more detail ( Table 6). This maturity concept applicable to drainage systems has been tested by Dill et al. (2019) [180] for the marginal zones of the Variscan Orogen of Mesoeuropean landscape in Morocco as well as by Dill et al. (2020c) [84] for the Mesoeuropean landscape of the Gondwana Continent in the South Atlantic Ocean (Figure 1c). Occurrence of weathering (see Figure 15) Absent (only slightly altered parent material) moderate argillitization Argillitization >> duricrusts Duricrusts-argillitization

The SSS in View of Higher and Lower Latitude Climate Zones
The modern SSS which in a Variscan-type orogen today is active in a humid mid-latitude climate zone in Central Europe ( Figure 18). What is the drainage system like in a Variscantype orogen located higher and lower latitudes than the archetypes (Figures 1c and 18)? For reference, the SSS equivalent drainage systems at higher latitude under dry continental to periglacial climatic conditions an area near Lake Hovsgul, Mongolia, has been considered. The landscape formed part of geomorphological transect from the taiga to the desert [42,183]. The drainage pattern of the very coarse-grained streams develops in a transition from the mountain forest steppe into a grass steppe. The channel systems is trellis-like with coarse grained -meandering streams, wide flood plains and terraces on both valley sides but of ephemeral type. Low-sinuosity to straight channel systems are consequent dip-stream systems. The higher sinuosity values are encountered in subsequent strike stream systems. In the Moroccan Meseta reference sites at lower latitude which evolved under tropical semiarid climatic conditions strike streams and anti-dip streams that break through the steeply dipping carbonate ridges evolved in some kind of a badland landscape [180,184,185]. Similar to the archetype in Germany and irrespective of the climate and latitude the typical SSS evolved where the tectonic drivers be it folding or faulting create the highest gradient and where a marked contrast in rock strength occurs. The sediment load, however, is controlled by the country rocks and the (paleo)climate which impacts on the lithosphere via the pedosphere and supergene alteration.

•
The "key message" to the evolution of the SSS is as follows (Table 5, Figure 2b): The SSS analyzed in this study are metastable drainage systems which are maintained over a longer period of time (approx. 25 Ma) and distance from source to confluence (approx. 30 km) by a rigid framework provided by the tectonics and bedrock lithology. The SSS is kept operating by a climate change from humid tropical through periglacial into humid temperate. In terms of chronology, it is most suitably described by a four-stage succession.

1.
The geological framework of the SSS: Forming the lithological and structural features in the bedrock as a result of different temperature, pressure and dynamic metamorphic processes.

2.
Prestage of SSS: Forming the paleo-landscape with a stable fluvial regime as starting point for the SSS 3.
Proto-SSS: Transition into the metastable fluvial regime of the SSS 4.
Modern SSS: Operation of the metastable fluvial regime with the tendency to regain a stable hydrodynamic state in the adjacent foreland • Tectonics plays a dual role. Late Paleozoic fold tectonic affects the Paleozoic bedrock and creates the basis for the drainage pattern of these SSS. Its influence diminishes from the very-low-grade to the high-grade stage of metamorphism (VLG⇒LGG⇒LHG). So does its guiding effect on the emplacement of the drainage systems during the Neogene and Quaternary. Late Cenozoic fault tectonic triggered the onset of the SSS to incise into the basement. The impact of fault tectonic on the SSS depends on the vertical displacement that correlates with the fluvial gradient of the SSS, the development of knickpoints, the occurrence of intermediate sediment traps and the cross-sectional shape of valleys.

•
The change in the bedrock lithology has an impact on the fluvial and colluvial sediments and their pertinent landforms. From slate to mica-gneiss (VLG⇒LGG⇒LHG) the vulnerability to chemical weathering increases and as a consequence of this a change among the process-related sediments takes place: gravel ⇒ clay + silt. The change of landforms follows suit: drainage system of straight ⇒ higher sinuosity and paired terraces ⇒ hillwash plains.

•
The climate change has a mediating effect via the bedrock controlling the intensity of mechanical and chemical weathering. The first type of supergene alteration increases towards younger intervals, whereas the second one shows a decreasing trend. Together with the Cenozoic tectonic, climate governs the intensity of fluvial, colluvial and cryogenic land-forming processes.

•
The elaborated genetic model of the SSS is restricted to the basement of the Variscan -Type orogens which are representative of an intermediate maturity state between the pre-mature Cenozoic fold belts and the super-mature Precambrian cratons. • Geomorphological-geological mapping of the fluvial, colluvial and cryogenic landform types is the descriptive part of the study while capturing numerical data relevant for the hydrographic studies of the SSS each reference area forms the geomorphometric one: (1) Quantification of fluvial and colluvial deposits along the drainage system, (2) slope angles, (3) degree of sinuosity as a function of river facies, (4) grain size distribution, (5) grain morphological categorization, (6) grain orientation ("situmetry"), (7) channel density, (8) channel/floodplain ratios.

•
The compositional approach involves new mineralogical for the regolith which has been subdivided into a zone of argillitization and duricrusts. Their mineral assemblages are the basis for the thermodynamic computations (Eh, pH, concentration of solubles) to constrain the paleoclimatic regime during formation of the SSS.

•
Comparing Variscan-type drainage systems from higher and lower-latitude climate zones relative to the archetypes reveals, the physical part (landform) of the SSS is controlled by tectonic and bedrock, whereas its chemical part (sediment composition) by the climate/weathering interacting with the country rocks.