Thermochronological constraints on the post-Variscan exhumation history of the southeastern Bohemian Massif (Waldviertel and Weinsberg Forest, Austria): palaeogeographic and geomorphologic implications

Resolving the Mesozoic and Cenozoic palaeogeography and geomorphologic development of outcropping Variscan basement is a pre-condition for the understanding of central European geodynamics. For our study, we have applied apatite fission-track (AFT) and apatite (U–Th)/He (AHe) thermochronology to surface rocks of the southeastern Bohemian Massif. 46 samples were examined by the AFT method. Additional AHe dating was applied to 8 of them. The AFT ages range from 251 ± 46 to 60.2 ± 4.8 Ma. AHe ages range from 25 to 525 Ma with rather high intra-sample scatter. On a regional scale, the AFT ages generally decrease from mainly late Variscan in the NE to Late Cretaceous and Paleocene in the SW. This regional age asymmetry relative to the NW–SE trending watershed of the Weinsberg Forest is neither compatible with regional uplift of a single block nor with large-scale lithospheric updoming. The lack of age breaks along late Variscan faults demonstrates that strong vertical offset cannot have occurred in Cretaceous and Cenozoic times. Inverse modeling of thermochronological data indicates regional Early Cretaceous cooling and subsequent reheating during the Late Cretaceous. Rocks of the present-day surface were heated up to a temperature of ca. 80 °C without full reset of the AFT system. This thermal history is compatible with the existence of a large mainland in Early Cretaceous times and a subsequent sedimentary reburial until the Campanian on the order of up to 1 km overburden. Parts of the exhumed weathering basal relief to the N and NE of the Weinsberg Forest are inherited as ‘sealed relief’ from Middle Cretaceous time. The observed regional asymmetry of AFT data is best explained by the development of a continental escarpment adjacent to the North Penninic Ocean in latest Cretaceous to Paleogene times. A final episode of accelerated cooling after ca. 20 Ma, as indicated by thermochronological modeling, is tentatively ascribed to either collisional coupling of the Alpine-Carpathian nappe pile with its northern foreland or to East-Alpine slab detachment.


Introduction
Post-Variscan thermal evolution and exhumation history of the Bohemian Massif has already been investigated in several areas. Vercoutere (1994), Wagner et al. (1994) and Hejl et al. (1997) have applied apatite fission-track (AFT) dating to the Naab mountains and to the basement adjacent to the Franconian Fault in Bavaria. Glasmacher et al. (2002) investigated the thermal history of Upper Proterozoic to Upper Carboniferous formations of the Teplá-Barrandium in the Czech Republic based on modeling of AFT data. Such data and corresponding modeling results are also available for the Ruhla crystalline complex in Thuringia (Thomson and Zeh 2000), for the Erzgebirge (Ventura and Lisker 2003) and the Lusatian Block in Saxony (Ventura et al. 2009), for the Sudetes (Aramowicz et al. 2006;Danišik et al. 2012;Sobczyk et al. 2020), for the south-eastern Bohemian Massif in Austria (Hejl et al. 2003) and for the Bavarian Forest (Siebel et al. 2009;Vamvaka et al. 2014). Both zircon (U-Th)/ He thermochronology and AFT analysis were applied to Lower Carboniferous greywackes of the Moravo-Silesian Culm Basin (Botor 2014;Botor et al. 2017). In a recently published, comprehensive study, von Eynatten et al. (2021) have presented new apatite fission-track and (U-Th)/He data for more than 100 samples from central Germany (eastern Rhenish Massif, Harz Mountains, Thuringian Forest, western Bohemian Massif, Franconian Platform, Buntsandstein uplands in northern Hesse and southern Lower Saxony). They interpreted these data together with several hundred, formerly published AFT data from the Rhenish Massif, the Vosges, the Black Forest, the entire Bohemian Massif and the Holy Cross Mountains in Poland. Von Eynatten et al. (2021) concluded that Late Cretaceous to Paleocene exhumation and cooling predominates in many areas and that the spatial exhumation pattern is governed by thrust-bordered basement uplift and superimposed regional domal uplift. In addition to crustal thickening by intraplate compressional stresses, also upwelling asthenosphere is supposed to have contributed to the total amount of uplift and exhumation.
The present investigation aims at a better understanding of the post-Variscan history of the Waldviertel and the Weinsberg Forest (Austria), including its regional exhumation and its tectonic and geomorphologic development. The pilot study of Hejl et al. (2003) has already detected some regional variations in the area's cooling history, but sampling was restricted to an E-W profile along the Kamp valley and a small region further south, in the vicinity of the towns of Krems and Spitz. Triassic AFT ages were found for the Moravian Thaya window, and mainly Cretaceous ages further to the W and S. No discernible age breaks were observed at the Diendorf and Vitis Faults. Thus, significant vertical offset due to these faults is precluded for Cretaceous and Cenozoic times. Furthermore, an integrated study of Eberts et al. (2021) based on Bouguer gravity anomaly, topographic data, low-temperature thermochronology and metamorphic grades in the SW Bohemian Massif has revealed a mainly late-to early post-Variscan age of basement segmentation and differential uplift along the NW-SE striking Pfahl and Danube shear zones in Bavaria. Evidence for Cretaceous and/or Cenozoic vertical offset was found for the NNW-SSE striking Cham Fault.
Thermochronological modeling based on AFT ages and track length data has shown that a regional sedimentary reburial of the Waldviertel basement during the Late Cretaceous is consistent with the data set (Hejl et al. 2003).
However, this area was not well covered by thermochronological data-especially samples from the eastern Weinsberg Forest and the Danube tributaries to its S and SW were missing. With 46 instead of 19 samples, the data set of the present investigation covers a larger area and yields a better spatial resolution of the regional cooling trends.
Besides a better resolution of post-Variscan exhumation trends, our study intends to set up well-defined boundary conditions for the thermal evolution, vertical movements, long-term denudation and burial of the Bohemian Massif in a wider, central European context. The post-Variscan cooling history of crustal segments depends on heat flow and geothermal gradients, as well as on exhumation produced by either large-scale domal uplift or by more localized processes like thrusting or fault-scarp tectonics (or a combination of both). Nielsen et al. (2005) have suggested that plate-wide stress relaxation could explain Paleocene basin inversions in Europe. Regional updoming may be due to tectonic thickening by compression, to lithospheric folding (Cloetingh and Burov 2011;Cloetingh et al. 1999), to underplating by basaltic melts (Brodie and White 1994), to thinning of the lithospheric mantle (Meier et al. 2016) or to dynamic topography caused by mantle upwelling (Braun et al. 2013). Von Eynatten et al. (2021) have argued that crustal shortening by compression would have been insufficient to explain the full amplitude of the Late Cretaceous and Paleogene exhumation (often in excess of 4 km) in central Europe, because crustal thickening associated with inversion tectonics can only contribute to half of the doming signal. Therefore, they suggested additional uplift caused by upwelling asthenosphere. We suppose that a denser net of thermochronological data will provide a better resolution and understanding of such processes.
Our new thermochronological age constraints are combined with compiled data from regional stratigraphy, structural geology and geomorphology to deduce the cooling and exhumation history of this former passive margin and its hinterland. Regarding the area under discussion (Waldviertel and the area to the east of Linz), the reader is referred to paleogeographic syntheses of Ziegler (1990), Faupl (1997 and Park (2014). The state of geomorphologic research in the granite domains of the Waldviertel has been reviewed in detail by Huber (1996) who specially emphasized features of the exhumed weathering basal relief. Regional age distributions and the modeled cooling histories that were found by the present investigation, partly confirm former concepts but also depict some trends of the macro-relief development, which were unknown so far. The most important finding is the discovery of an ancient escarpment at the southwestern edge of the Weinsberg Forest. This escarpment has developed adjacent to the passive continental margin of the North Penninic Ocean in Paleogene times. It formed the transition zone between the shelf area (underneath the present-day 1 3 North Alpine Molasse) and the exhumed weathering basal relief of a Mesozoic peneplain on top of the Weinsberg Forest and further north, in the Kamp and Thaya catchment area.

Basement geology and geomorphological features
The Variscan basement is mainly composed of Moldanubian nappes, Moravian units of the Thaya Window (adjacent to the Diendorf Fault), and plutonic rocks of the South Bohemian batholith ( Fig. 1; cf. also Schnabel et al. 2002). A Grenvillian cratonic provenance of the original Moravo-Silesian crust was suggested by SHRIMP dating of inherited zircons in metagranitoid rocks which has revealed Mesoproterozoic and late Paleoproterozoic tectono-magmatic events (Friedl et al. 2000). On the other hand, the Moldanubian units are characterized by Neoproterozoic and Eburnian (middle Paleoproterozoic) ages of inherited zircon populations (Friedl et al. 2000(Friedl et al. , 2004Finger et al. 2007).
The Upper Moravian nappe comprises the Bittesch orthogneiss with U/Pb age 580 ± 7 Ma (Friedl 1997). Eastward thrusting of Moldanubian units onto Moravo-Silesian crust occurred in the Early Carboniferous and produced a distinct metamorphic zonation in the lower plate (Höck 1975(Höck , 1994. The Moldanubian zone of the Waldviertel consists of a Variscan nappe pile. Remarkable formations of this nappe pile are the Dobra gneiss (Drosendorf nappe) with a protolith age of 1.38 ± 0.01 Ga (Gebauer and Friedl 1993;Lindner et al. 2020), and the granulites on top of the Gföhl nappe, which were produced by a high-pressure, high temperature (HPHT) metamorphism at ~ 340 Ma (Kröner et al. 1988;Friedl et al. 1994). Intra-Moldanubian nappe stacking slightly postdates this metamorphic event but was finished before 335 Ma (Friedl et al. 2011;Finger et al. 2007).
Several generations of granitoid intrusions build up the South Bohemian Batholith (Frasl and Finger 1991). The volumetrically dominant, coarse-grained Weinsberg granite intruded between ca. 331 and 322 Ma (± ca. 3 Ma each; Finger et al. 2022). The Eisgarn granite with S-type composition and some fine-grained granites with I-type composition postdate the emplacement of the Weinsberg granite and occupy smaller areas of the batholith (Frasl and Finger 1991;Finger et al. 2022).
Prominent late Variscan faults dissect the Bohemian Massif (Fig. 1a). The area of investigation (Fig. 1b) is affected by the SW-NE trending Diendorf, Vitis, and Rodl Faults. Their shear zones that are exposed at the present-day surface, were formed at different synkinematic paleotemperatures. Mylonites of the Rodl Fault were formed under ductile greenschist facies conditions, but rather cool brittle deformation features are observed along the Diendorf Fault (Brandmayr et al. 1995;Tollmann 1985). The latter has a length of more than 80 km extending from the Alpine foreland in the SW to the southern end of the Boskovice furrow in Moravia. It produced a sinistral strike-slip displacement of almost 60 km (Schermann 1966;Finger and Riegler 1999). Field evidence for Quaternary faulting ca. 25 km to the SW of Brno (Moravia, Czech Republic) was reported by Leichmann and Hejl (1996).
The post-Variscan sediments covering the margins of the Bohemian Massif will be considered in a following section (Mesozoic and Cenozoic paleogeography) with regard to potential reburial and regional erosion. Intensive tropical weathering, strong erosion and peneplanation affected the Bohemian Massif in Meozoic and Cenozoic times. Several stages of peneplanation interacted with shallow-marine ingressions during the Late Jurassic, the Middle to Late Cretaceous and the Miocene (cf. Hejl et al. 2003).
Woolsack weathering (Pötzsch 1803), exhumed inselbergs, exfoliation domes (Gilbert 1904) and castle koppies are common surface features of the South Bohemian batholith (cf. Huber 1996, and references cited therein). Giant boulders and wobbling stones (so-called 'Wackelsteine') occur in the vicinity of Gmünd (Kieslinger 1964;Köhler 1976) and on many other granite exposures of the Waldviertel. The woolsacks result from blocky dissection by tension cracks (Cloos 1925) and sub-horizontal exfoliation (Gilbert 1904), and from subsequent hydrolysis of feldspar that starts at the joints, softens the rock and produces granular disintegration (Wilhelmy 1981). Strong hydrolysis mainly takes place in the zone of groundwater saturation, below the water table. Resulting deep weathering profiles of tropical peneplains can be tens or even hundreds of meters thick (Ollier 1988). It is generally assumed that such deep weathering and peneplanation of Variscan granite in central Europe has developed in a humid tropical climate (Borger 1992). Thus, the present-day weathering basal relief (Büdel 1957) should be inherited from Early Cenozoic (pre-Miocene) or even earlier times. Preservation of Mesozoic weathering profiles until present is indeed known from Australia (Ollier 1988). On the other hand, geomorphological investigations in the Aegean archipelago (Greece) have shown that peneplanation can proceed rather fast (within less than 10 Ma) when climatic conditions are appropriate (Riedl 2004).
Relict planation levels are widespread in the southern Bohemian Massif. Fischer (1979) interpreted them in terms of a peneplain staircase ('Rumpfflächentreppe') that corresponds to the concept of a piedmont staircase (Spreitzer 1951;cf. also Riedl 2007, and references cited therein). He distinguished four palaeoplanation levels that he ascribed to post-Cretaceous to pre-Eggenburgian piedmont steps. The Weinsberg peneplain staircase ('Rumpfflächentreppe im Weinsberger Wald') comprises the oldest planation relics on top of the Weinsberg Forest (> 940 m, 'Hochfläche von Bärnkopf'), and three lower planation levels which become progressively younger with decreasing altitude. These lower stages of peneplanation occur between 850 and 880 m, between 760 and 810 m, and between 720 and 750 m. These palaeoplanation levels are supposed to occur at quite similar altitudes on both sides of the NW-SE trending watershed of the Weinsberg Forest. Younger levels of pedimentation between 550 and 400 m altitude, in the vicinity of the river Danube (Donau) and its tributaries, are supposed to have a Neogene age.  Glasmacher et al. 2002) and geological map (b) of the southeastern Bohemian Massif and its surroundings (simplified after Schnabel et al. 2002 Quaternary terraces of the Danube are located below 400 m altitude. Two essential aspects of the geomorphological concept of Fischer (1979) should be born in mind. First, the NW-SE trending ridge of Weinsberg Forest is not a first-order continental watershed. The continental divide between the North Sea catchment (Vltava-Elbe) and the Black Sea catchment (Donau) is farther to the N, close to the national border between the Czech Republic and Austria. This continental divide is unspectacular from a geomorphological point of view. The Weinsberg Forest is even not a second-order watershed, like those between the catchments of major rivers as Rhine, Elbe or Oder. It just constitutes a local divide between the small rivers Krems and Kamp to the NE, and even smaller rivers to the SW (Gusen, Aist, Naarn, Isper, etc.). The horizontal distance between their headwater sources and their confluence into the Danube is less than 80 km. A second important fact is that the most spectacular giant woolsacks, mushroom rocks, wobbling stones and castle koppies do not occur on top of the Weinsberg Forest ('Zentrales Bergland' above 1000 m, and 'Hochfläche von Bärnkopf' at 940-980 m altitude) but farther to the N at altitudes below 600 m-especially in the nature park 'Blockheide Eibenstein' close to the town of Gmünd (Kieslinger 1964;Huber 1996).

AFT dating and track length measurements
46 samples were collected from the South Bohemian batholith and from metamorphosed rocks of the Moldanubicum and the Moravo-Silesicum. Except for one paragneiss sample (BM 7), all samples originate from plutonic rocks or metamorphosed plutonic rocks. Locations and rock types are indicated in Supplement file 1. For a better comparison with the already published dataset of Hejl et al. (2003), we have merged the 19 old and the 27 new samples in common tables (Tables 1, 2, and 3) and in Fig. 2.
After crushing and dry sieving, the fractions with grain sizes between 100 and 250 µm were chosen for mineral separation. Apatite was enriched by chemical flotation (Hejl and Ney 1994), by magnetic and heavy liquid separation.
The grain-population method (cf. Gleadow 1981; Wagner and Van den haute 1992) was applied for dating, because it is well suitable for apatite with fairly homogeneous uranium concentrations, as it can be assumed for apatite of magmatic origin. An aliquot of each apatite concentrate was heated for 24 h at 550 °C in a muffle furnace and then irradiated with thermal neutrons. Samples BM 47 to BM 58 were irradiated in channel 8 of the Thetis nuclear reactor at the University of Ghent (Belgium). The other samples were irradiated in the Astra reactor of the Forschungszentrum Seibersdorf (Austria). Both irradiation facilities had a good thermalization. Thermal neutron fluences were between 3.7 × 10 14 and 2.2 × 10 15 cm −2 (Table 1), and the epithermal neutron flux amounted to less than 0.1% of the total flux. Thermal neutron fluences were determined by γ-spectrometry of activated Co and Au monitors (Van den haute et al. 1988). Natural and irradiated grains of each sample were embedded separately in epoxy resin, ground, polished and etched with 5% HNO 3 for 60 s at 20 °C. The 238 U spontaneous fission decay constant λ f = 8.46 × 10 -7 a −1 (Galliker et al. 1970) was used for age calculations. System calibration was repeatedly tested with Durango apatite (Jonckheere et al. 1993(Jonckheere et al. , 2015 and apatite from Mt. Dromedary banatite (Green 1985). Dating results are reported in Table 1.
Polished apatite mounts of 27 unirradiated samples with spontaneous track densities above 7 × 10 5 cm −2 were overetched with an aqueous solution of 5.5 M HNO 3 for 20 s at 21 °C, and then examined for fully confined, horizontal tracks (Gleadow et al. 1986a, b). Length measurements were performed with a drawing tube attachment to the microscope and a digitizing tablet (cf. Wagner and Van den haute 1992). Individual track length measurements were performed with a precision of ± 0.2 µm. Results of confined track lengths measurements-including mean lengths, standard deviations and Dpar-are reported in Table 2. Frequency distributions (histogram classes) are reported in Table 3. A mean induced track length of 16.2 µm was determined on a sample of irradiated Durango apatite age standard (Jonckheere et al. 1993(Jonckheere et al. , 2015. This value is similar to that reported by Gleadow et al. (1986b).

Apatite (U-Th)/He dating
AHe analysis was conducted on a subset of samples that had been previously dated with AFT. Three samples (BM 1,5,11) were analyzed at the University of Göttingen (Germany). Five samples (BM 14,40,43,46,54) were dated at Dalhousie University (Canada). Results are given in Tables 4  and 5. For a detailed description of AHe dating, the reader is referred to a recent extensive review in Flowers et al. (2022a, b). Details on the analytical procedure of Göttingen and Dalhouse laboratories are provided in supplementary file 2. Here, only a brief outline is given.
AHe thermochronology is based on the production of α-particles (i.e., He atoms) during radioactive decay of 238 U, 235 U, 232 Th and 147 Sm. Retention of these 4 He nuclides within the crystal lattice versus their loss by diffusion is temperature dependent and closure temperatures may range between ∼ 40 and ∼115 °C.
He diffusivity is impacted by various factors, among which radiation damage (e.g., Shuster et al. 2006;Flowers et al. 2009;Gautheron et al. 2009), and grain size (e.g., Farley 2000) are considered as the most important ones. Grain size determines how much He nuclei are ejected from the crystal, since the α-particle travels through the grain upon its generation (the alpha stopping distance in apatite is on average 18 µm, Flowers et al. (2022a). He atoms produced by parent nuclides that are located in the proximity of the crystal's rim are more likely to get ejected. Therefore, AHe ages are generally corrected for He loss using a geometric correction factor (Ft), which takes the size and the shape of each dated crystal into account (Farley et al. 1996). The kinetic differences in He diffusion between individual grains are particularly important for samples experiencing protracted cooling, i.e., with a sufficiently long residence time in the He partial retention zone to induce partial diffusive He loss. For datasets with high intra-sample dispersion, such as the one presented within this study, single-grain ages are plotted vs. equivalent spherical radii (ESR; proxy for grain size) and effective uranium (eU; proxy for radiation damage calculated as eU = U + 0.235·Th + 0.0046·Sm). If correlations are observed, even highly dispersed datasets can enhance thermal history models, as kinetic models incorporate He diffusion as a function of damage accumulation and annealing (e.g., Flowers et al. 2009;Gautheron et al. 2009). If no such correlation exists, we here follow the approach of not including AHe data into the interpretive process of thermal history modeling and only report the dataset. We define a threshold level of 20% dispersion, below which we calculate and report an unweighted mean age and the standard deviation. Following these criteria, it turned out that only one sample remained to be integrated into thermal history modeling implementing the RDAAM annealing code of Flowers et al. (2009).
Analysts: Ewald Hejl (EH) and Gert Sekyra (GS). The age is given as an error-weighted average for those samples that have been counted by both analysts. N s , N i , total number of counted tracks (spontaneous and induced, respectively); ρ s , ρ i , areal track densities (spontaneous and induced, respectively)

Regional pattern and frequency distribution of measured AFT ages
The AFT ages range from 251 ± 46 Ma (Latest Permian) to 60.2 ± 4.8 Ma (Paleocene). We observe a decrease of measured ages from NE to SW, but no significant age breaks (i.e., exceeding the error limits) at the Late Variscan Vitis and Diendorf faults (Fig. 2). Thus, strong vertical offset (> ca. 500 m) on these is precluded for Cretaceous and Cenozoic times. This has already been stated by Hejl et al. (2003). The highest AFT age was that of sample BM 53 which originates from a Gföhl gneiss between the towns of Horn and Gföhl. Seven other AFT ages above 200 Ma were measured for surrounding locations (BM 12, 47 and 54) and for the area farther to the N (BIT 1, BM 49, 57 and 58). The samples of this highest age group originate from Gföhl gneiss (3 samples), granulite (2 samples), Dobra gneiss (1 sample) and Bittesch gneiss (2 samples). No age break occurs at the thrust front of the Moldanubian onto the Moravian units. Within the area with highest AFT ages (> 200 Ma), also lower ages were measured. The lowest age of this area was found for sample BM 10 (109.2 ± 10.6 Ma).
Because of the small variation in sampling altitudes (from 280 to 551 m) compared to the much larger horizontal distances between sampling locations (up to ca. 40 km), a plot of AFT ages versus altitude has not been made because it could be misleading for the area under discussion.
A quite homogeneous group of Early Cretaceous AFT ages (between 142 ± 17 and 104 ± 7.6 Ma) dominates the summit region and the northeastern slopes of the Weinsberg Forest, down to the Kamp valley (Fig. 2) to the E of the town of Zwettl. Two samples to the W and S of the town of Groß-Gerungs (BM 2 and 3) yielded slightly lower AFT ages of 91.9 ± 13.5 Ma and 98.5 ± 5 Ma, respectively (Fig. 2). Thus, the whole area exhibits only moderate age variations over a horizontal distance of ca. 60 km from NW to SE.
Distinctly lower AFT ages were measured for the samples from the SW of the Weinsberg Forest, i.e., from an elevation of 750 m (BM 39) down to the slopes of the river Danube at ca. 300 m altitude (BM 43). Ages of this group range from 104.9 ± 13.6 Ma to 60.2 ± 4.8 Ma (Middle Cretaceous to Paleocene). They are significantly lower than the AFT ages of the samples from similar elevation on the northeastern side of the Weinsberg Forest.
The measured AFT ages exhibit an Early Cretaceous frequency maximum between ca. 145 and 100 Ma (Fig. 3). Most AFT ages of samples from the northeastern slopes of the Weinsberg Forest belong to this group.

AFT length data
Frequency distributions of confined track lengths (Gleadow et al. 1986a, b) as well as corresponding mean track lengths and standard deviations were determined for 27 samples (Figs. 4, 5, 6 and 7). Only those samples with at least 80 confined track lengths measurements were chosen for thermochronological modeling. To facilitate their interpretation, the histograms are grouped according to the AFT age classes of Figs. 2, 3 and 4 (i.e. > 200 Ma;< 100 Ma).
Arithmetic means of all length distributions range between 11.99 µm and 13.40 µm. According to Gleadow et al. (1986b), values around 12 µm or 13 µm are typical for basement terrains, which have never been significantly disturbed thermally since their original magmatic emplacement or metamorphism. Standard deviations of all length distributions range from 1.38 to 2.07 µm.

Apatite (U-Th)/He results
All analyzed samples delivered inclusion-free, fracturefree, euhedral grains of sufficient size with sphere radii well above 50 µm. For the Göttingen sample subset, 3-4 grains per sample were measured (Table 4), for the Dalhousie sample subset 5 crystals per sample were measured (Table 5). Intra-sample data variation is high for almost all samples with individual dates ranging from 25 ± 0.1 to 525 ± 42 Ma and dispersion varying from 16 to 64%. Supplementary file 3 displays age-eU and size-eU plots for all samples. Since we neither find a positive correlation between age and eU nor between age and grain size, we cannot pinpoint the firstorder kinetic controls for date dispersion and thus refrain from including these samples into thermal history modeling. We report a mean age for sample BM 5, since its dispersion is 16% only. The AHe mean age of 110 ± 18 Ma is younger than the AFT age of 122 ± 7 Ma. There is no standard protocol on when and how to combine individual dates, but we here apply a threshold of 20%. In the absence of positive relationships, He implantation or uncharacterized eU zonation may be responsible for data scatter, since He ejection correction assumes a homogeneous distribution of the parent nuclides. However, our AFT analysis revealed very little internal zoning.

Mesozoic and Cenozoic palaeogeography
For a better understanding of the central European paleogeography, we have added the sketch maps of Fig. 8, which roughly display the regional extent of mainland, marine and coastal areas for selected periods of the Mesozoic and Cenozoic. These sketch maps were compiled after Ziegler (1990), Faupl (1997) and Park (2014).
In the following, we briefly outline the waxing and waning of seas in our study area as derived from the sedimentary record. Since discordant sedimentary horizons pin surface (re-) exhumation of basement rocks, they serve as first-order constraints for thermal modeling of thermochronology data, as will be explained in the next chapter.
Sediments of the Boskovice furrow in Moravia (Fig. 1a) are of Late Carboniferous to Permian age. Similar sediments occur farther south at the village of Zöbing, where a tectonically isolated block of arkoses, sand-, and siltstones is squeezed between two segments of the Diendorf fault (cf. Fig. 1b). The Late Carboniferous to Early Permian Zöbing formation has a thickness of ca. 1000 m (Vasicek 1991). This isolated occurrence of thick pre-Mesozoic sediments demonstrates that late-Variscan exhumation of basement rocks was not regionally uniform, but was complicated by a mosaic of mountain ranges and intermontane basins (Fazlikhani et al. 2022;Ziegler 1990).
In Triassic times, the Bohemian Massif was part of the Vindelicic-Bohemian mainland, which separated the shallow-marine, epicontinental German basin from the northern Tethys (Fig. 8a). Coastal fans, fluvial and lacustrine sedimentation were restricted to the western margin of the Bohemian Massif (Schröder et al. 1997), while most parts of the regions farther east were mainland. A Jurassic transgression of the northern Tethys inundated the southern shore of the Vindelicic-Bohemian mainland and culminated in Late Jurassic time. Marginal parts of the Bohemian Massif became submerged, and a marine channel separated its major part from a Sudetic island in the N (Fig. 8b). The preserved thickness of Jurassic formations in the area of the later North Alpine molasse zone locally exceeds 500 m (Malzer et al. 1993). Significant reburial of the southeastern Bohemian Massif by Late Jurassic sediments may have occurred too, because gravels of cherty limestones of presumable Late Jurassic age were found in Neogene sediments to the NE of the town of Krems (Steininger and . However, position and shape of the Late Jurassic shoreline are unknown. Because of a strong regression at the end of the Jurassic, many regions of western and central Europe became part of a large mainland. In the Early Cretaceous at ca. 130 Ma, it included the Brabant Massif, middle and northern Germany, Bohemia and Moravia (Fig. 8c). Limnic and brackish sedimentation-including coal measures-occurred along the coast, while central parts of the mainland were subjected to erosion under a humid tropical climate (cf . Ziegler 1990;Faupl 1997). This period of intensive tropical weathering on a broadly emerged mainland temporally coincides with the frequency maximum of AFT ages between ca. 145 and 100 Ma that we interpret as cooling ages produced by strong erosion and continuous cooling through the entire partial annealing zone (PAZ) for apatite fission tracks.
Extensive transgressions occurred during the Middle Cretaceous (Niebuhr et al. 2009). Major parts of the former mainland became progressively inundated. The sea-level rise culminated about 75 Ma ago, when southern Bohemia and the Rhenish Massif were still connected as a narrow elongated island (Fig. 8d). Several hundred meters of Late Cretaceous marine sediments accumulated in central and northern Bohemia (Ulčný 2001;Ulčný et al. 2009), as well as in the shelf area of the later North Alpine Molasse Zone (Malzer et al. 1993). The sedimentary basins of Ťrebon and Česke Budejovice to the NW of the town of Gmünd contain non-marine Late Cretaceous sediments with a thickness of 450 m (Suk 1984). Spectacular giant woolsacks of the nature park 'Blockheide Eibenstein' are located directly at the sedimentary base of Late Cretaceous sediments. Thus, they may belong to a 'sealed relief', i.e., a tropical weathering basal relief that could have developed during the Early Cretaceous and afterward was protected from further erosion by a thick cover of Late Cretaceous sediments. Such Late Cretaceous reburial has probably affected major parts of the Waldviertel (Hejl et al. 2003).
Paleocene and Eocene sediments are unknown from the Waldviertel. Most likely, the whole region was emerged und subjected to erosion . A Late Eocene transgression occurred farther south. Corresponding shallow marine and adjacent deltaic sediments of up to 120 m thickness are known from boreholes in the Upper Austrian molasses basin (Malzer et al. 1993). The shoreline of a Late Eocene sea was close to the southern margin of the present-day outcropping basement (Fig. 8e).
Strong subsidence of the North Alpine molasses basin started in Oligocene time, after the closure of the North Penninic Ocean. Consequently, the sea advanced toward the N and the shoreline touched the present-day border of the outcropping Variscan basement in Upper and Lower Austria (cf. Kuhlemann and Kempf 2002). Coastal sands of the Linz and Melk formations lie directly on top of the Variscan basement. They comprise brackish to shallow-marine sands and clays, as well as coal-bearing freshwater deposits (Malzer et al. 1993). Fluvial and lacustrine sediments of Oligocene age are even preserved deep inside the Waldviertel basement area. These silty sands, gravels and clays of the St. Marein-Freischling formation (cf. Roetzel and Steininger 1996) occur between the towns of Gmünd and Horn, in close vicinity to the headwaters of the Thaya and Kamp rivers. These sediments temporally constrain emergence of the basement units, but due to their limited thickness of less than 40 m (cf. Tollmann 1985), deep reburial is ruled out.
Few Badenian and Sarmatian marine to brackish sediments are preserved in the vicinity of the town of Krems, but in general, the Miocene sea did not submerge the presentday basement area and no transgressions occurred after the Sarmatian . Ongoing uplift since the Pannonian (about 10 Ma ago) created an eastward-orientated fluvial drainage of a Palaeo-Danube from Upper Austria to the Vienna basin (Unger 1989;Gyurits and Kurzweil 1976). A general drainage of the North Alpine Molasse Zone (including southern Germany) from W to E was established during the Late Pannonian or at the beginning of the Pliocene, i.e., about 7-5 Ma ago (Unger 1989). Quaternary gravel terraces of the Danube represent climatically and uplift controlled steps of fluvial incision. The sedimentary record of preserved formations in the vicinity of the Bohemian Massif demonstrates that substantial reburial of the Variscan basement may have occurred during the Late Jurassic and/or Late Cretaceous. But the Early Cretaceous must have been a period of persistent regional erosion under a tropical climate. However, the total thicknesses of sedimentary overburden and precise delimitations of erosion and sedimentation areas are unknown. Thus, reconstructions of Mesozoic paleogeography yield only 'soft constraints' for thermochronological modeling. A rather 'strong constraint' is given for those samples, which were taken close to the base of the St. Marein-Freischling Formation (Oligocene). These rocks must have reached very shallow depth about 30 Ma ago, and significant reburial is excluded due to the formation's limited thickness.

Thermochronological modeling and general discussion
According to Gleadow et al. (1986b), mean confined track lengths around 12 and 13 µm and standard deviations between 1.2 and 2.0 µm are typical for undisturbed basement that underwent steady denudation.
Thus, our AFT ages, mean confined fission-track lengths (from 11.99 to 13.40 µm) and standard deviations (from 1.38 to 2.07 µm) do not contradict the assumption of slow steady cooling through the partial annealing zone (PAZ). However, the shape of some histograms is not negatively skewed, and some histograms are even bimodal (Figs. 5, 6, 7). Therefore, more complex cooling paths must be considered for some samples and age groups.
In a first interpretation step, we assume that the distinct frequency maximum of AFT ages between ca. 145 and 100 Ma represents cooling ages, i.e., that the samples underwent rather fast and continuous cooling through the whole AFT partial annealing zone between ca. 120 and 60 °C (cf. Wagner and Van den haute 1992). The AFT ages above 200 Ma probably represent post-Variscan cooling ages, while the AFT ages between 200 and 145 Ma may be mixed ages that are either the result of stepwise cooling or the result of moderate reheating between a late-Variscan and an Early Cretaceous cooling event. Such moderate reheating could have been caused by a Late Jurassic sedimentary reburial, as has been outlined in the previous chapter. Figures 4 and 6 show that the confined track lengths of most samples of the age group between 145 and 100 Ma have rather high arithmetic means above 12.8 µm, while those of the other age groups (> 145 Ma and < 100 Ma) are mainly lower than 12.8 µm. This underlines the foregoing assumption that the frequency maximum of AFT ages between ca. 145 and 100 Ma coincides with a regional cooling event during the Early Cretaceous. Samples of the lowermost age group (< 100 Ma) have mean lengths of less than 12.75 µm and-except for sample BM 36-rather high standard deviations (> 1.8 µm). This indicates that the samples with ages of less than 100 Ma have not cooled continuously through the whole PAZ in Cenozoic times but did already reside at temperatures below 120 °C during the Middle and Late Cretaceous.
The outlined Mesozoic and Cenozoic palaeogeography helps to formulate rather 'strong' and 'weak' constraints for modeling of thermal histories. The fact that a thin cover of Oligocene freshwater sediments occurs along the Thaya and Kamp rivers on top of the Variscan basement yields a strong constraint for the potential cooling paths of adjacent basement rocks. The basement underneath or close to the St. Marein-Freischling formation (mainly between the towns of Gmünd and Horn) must have been close to the surface at about 30 Ma ago. On the other hand, the Mesozoic palaeogeography of already eroded formations on top of the Bohemian Massif is not well known. An Upper Jurassic and/ or Middle Cretaceous to Campanian sedimentary cover may have buried marginal parts of the Bohemian Massif, but both sedimentary thicknesses and the position of ancient coastlines remain doubtful. Thus, potential reburial and reheating between ca. 100 to 80 Ma (Cenomanian to Campanian) may be considered as weak constraint. At least strong erosion of the hinterland is rather unlikely for the quoted time interval.
We used the program HeFTy (version 2.0.9, 2022 provided by Richard A. Ketcham) for inverse modeling of AFT data and AHe ages with dispersion below 20%. This program calculates the range of thermal histories (temperature versus time) that are potentially consistent with measured thermochronometric data (AFT age, frequency distribution of confined track lengths and AHe age) and with constraints that are imposed to the t-T path by use of independent geological evidence. Annealing kinetics of the AFT system rely on the investigations of Carlson et al. (1999), Donelick et al. (1999) and Ketcham et al. (1999), and on the improved model of Ketcham et al. (2007). Inverse modeling of AHe data uses the He diffusion kinetic model of Flowers et al. (2009). Random paths were generated by a Monte Carlo scheme  Ziegler (1990), Faupl (1997 and Park (2014). The hatched area corresponds to the presentday outcropping basement shown in Fig. 1 with random spacing of sub-segments. We have imposed a maximum |dT/dt| of 5 °C/Ma because exhumation and/ or sedimentary burial rates of more than ca. 200 m/Ma are improbable for a non-orogenic setting. We have estimated the apatite's resistance to track annealing by including the diameters of etch pits parallel to the crystallographic c-axis (Dpar values in Table 2).

Fig. 9
Inversion models (temperature versus time) generated with HeFTy software (version 2.0.9, Richard A. Ketcham) for samples BM 53,6,5,7,14,8,1,4,36,and 43 (in descending order of AFT ages). Technical specifications: annealing model of Ketcham et al. (2007); goodness-of-fit (GOF) by Kolmogorov-Smirnov test; Monte Carlo method of search; random subsegment spacing between constraints; |dT/dt|< 5 °C/Ma for each subsegment; cf. text for geological constraints. The inverse model of BM 5 includes the concordant (U-Th)/He age as additional thermochronometer (calibration after Flowers et al. 2009) We have not measured the angle between confined tracks and the crystallographic c-axis because the apatite grains were embedded in epoxy in random orientation. A systematic assessment of such angles would have restricted the length measurements to those grains viewed perpendicular to the c-axis. In our case, most grains have a low length-width ratio. The right orientation is rare, and only few confined tracks would have been suitable for both length and angle measurement. Thus, we opted for a modeling procedure without taking angles into account. In this case, the HeFTy modeling routine is based on the assumption that confined tracks were chosen randomly for length measurement-regardless of track and grain orientation-and that no systematic bias has occurred.
For each path, the resulting AFT age and length distribution were calculated, and the goodness-of-fit (G.O.F.) between calculated and measured data was evaluated by a Kolmogorov-Smirnov test. In the case of sample BM 5, concordant AHe data have been included in the inverse model, i.e., in the estimation of G.O.F. values, but the other inverse models include AFT data only. The program HeFTy maps out time-temperature regions with "good" and "acceptable" fit, corresponding to G.O.F. values from 0.5 to 1 and from 0.05 to 0.5, respectively. Inverse modeling was terminated when at least 100 "good paths" had been found (inversion models of Fig. 9).
Depending on the geological setting of the sample, constraints were imposed on the modeled thermal histories. These constraints are based on geological observations, such as sedimentary horizons discordantly overlying basement rocks, thus proving subaerial unroofing previous to the time of deposition. They are defined by a series of boxes in the time-temperature graph. The starting and end point of each path segment must be located in these boxes. For all model runs, a present-day surface temperature of 8 ± 5 °C was imposed as constraint. This temperature is an estimation of the mean annual temperature during the last ca. 10 ka.
We imposed 20-30 °C at 35-25 Ma for samples from locations close to discordant contact with the St. 7 and 8). Thermal histories of samples of the AFT age group between 145 and 100 Ma (BM 1,4,5,6,7,8 and 14) were allowed to reach low temperatures in the Early to Middle Cretaceous (between 145 and 95 Ma) and to be moderately reheated in the Late Cretaceous with a temperature climax between 80 and 65 Ma. These constraints were rather soft because of a broad temperature spread and a large time interval of the boxes (25-70 °C between 145 and 95 Ma;. In other words, Late Cretaceous reheating was enabled but not enforced. Stepwise cooling or even continuous cooling would have been compatible with these weak constraints. Two samples of the AFT age group with the lowest ages (BM 36 and 46) were modeled without the option of cooling to less than 60 °C in Cretaceous times. Inverse modeling with a constraining box of 90-150 °C from 145 to 65 Ma produced a much higher percentage of cooling histories with "good" fit than unconstrained modeling from 300 Ma to the present.
Inverse modeling of sample BM 53 ( Fig. 9) indicates fast post-Variscan cooling during the Late Carboniferous and the Early Permian. The sample resided at ca. 60 °C from 250 to 70 Ma, before cooling became faster again. Good fit was obtained without Late Cretaceous reburial, but the inverse model is compatible with stepwise cooling after ca. 70 Ma.
Inverse modeling of samples with Early Cretaceous AFT ages (BM 6,5,7,14,8,1 and 4, in descending order of the AFT ages) was performed with the weak constraints mentioned above that allow Early Cretaceous surface approach and subsequent reheating during the Late Cretaceous (Fig. 9). These models indicate that most samples did not cool to less than 100 °C before ca. 145 Ma ago. Fast cooling occurred during the Early Cretaceous. Afterward, the samples were either moderately reheated or resided at a fairly constant temperature until ca. 65 Ma ago. The inverse model of BM 5 with included AHe data most narrowly confines good fit and depicts distinct Late Cretaceous reheating. Cenozoic cooling was again faster. Some samples (BM 5, 7, 14 and 8) show a last step of accelerated cooling after ca. 30 Ma ago.
Inverse modeling of samples with Paleogene AFT ages (BM 36 and 46) was performed without Cenozoic time-temperature constraints, but cooling to less than 90 °C was prevented for the whole Cretaceous ( Fig. 9) because test runs without this constraint produced a lesser percentage of "good" fit. This supports the assumption that the samples resided in the PAZ until the end of the Cretaceous. The resulting models indicate stepwise cooling after ca. 80 Ma. Fast cooling from 80 to 60 Ma is followed by a period of moderate cooling from 60 to 20 Ma, and a second period of fast cooling after 20 Ma. The area of good fit would allow for a moderate reheating between 100 and 80 Ma but does not imply it.
Accelerated basement exhumation in latest Cretaceous or earliest Paleogene times is also indicated by independent geological evidence: Heavy mineral spectra in Late Cretaceous and Paleogene sediments of the Rhenodanubian flysch exhibit a distinct compositional change at the Cretaceous/ Paleogene boundary. Associations rich in garnet dominate in Campanian and Maastrichtian flysch sandstones (Woletz 1963). The source rocks of these garnets can be found in the metamorphosed units of the Austroalpine hinterland. This supply of garnet-rich detritus terminated at the end of the Cretaceous and was replaced by sandy material rich in zircon, typical for Paleogene sandstones of the Flysch zone and the Helveticum (Woletz 1963). This sudden change in detritus supply is best explained by the progressive closure of the North Penninic Ocean and the increasing terrigeneous input from its northern margin, i.e., from the European shelf and the adjacent crystalline-mainly granitic-basement of the Bohemian Massif. This scenario is compatible with the development of a continental escarpment along the NW-SE trending Weinsberg ridge. The southern margin of the Bohemian Massif was buried under a thick cover of marine and coastal sediments during the Late Cretaceous. The transgression culminated in the Campanian. Afterward, the sea level dropped, and the sedimentary overburden on the southern Bohemian Massif was removed by erosion. In large areas, Weinsberg and other Variscan granites were exposed during the Paleogene and delivered zircon-bearing sand to the European shelf, the Helveticum and turbiditic fans of the North Penninic Ocean. The continental escarpment on the northern (passive) margin of the Flysch Ocean developed from about 60 to 35 Ma ago, i.e., during a time span of ca. 35 Ma. This situation ended in the Late Eocene, when the Austroalpine nappe pile collided with the European shelf.
A last phase of accelerated cooling at ca. 20 Ma is implied by thermal simulations of samples 5, 36 and 43. Interestingly, the sedimentary archive of the eastern part of the Molasse basin records a change in the sedimentation pattern and onset of rapid basin infill at ca. 19 Ma (Kuhlemann 2007;Hülscher et al. 2019). Concurrent exhumation of the nearby Alpine source area is indicated by fully reset AFT ages ranging from ca. 18 to 25 Ma from the easternmost flysch zone W of Vienna (Trautwein et al. 2001). New AHe data from clastic formations of the nappe pile of the Northern Calcareous Alps in the easternmost Eastern Alps confirm such a young cooling pulse (Heberer et al. 2022). It may be correlated with thrust activity of the advancing nappe stack (e.g., Beidinger and Decker 2014), but it may also, at least partly, be the surficial expression of more deep-seated processes. Based on the AlpArray seismological campaign, Handy et al. (2021) proposed a new model of slab configuration beneath the Alps, including a completely detached slab beneath the Eastern Alps.
With regard to domal uplift on a lithospheric scale, as it has been proposed by von Eynatten et al. (2021) for major areas of central Europe, we have to consider the regional size and amplitude of the exhumed areas. The NW-SE trending ridge of the Weinsberg Forest has a length of about 80 km, and its width in NE-SW direction is even smaller (Fig. 2b). Thus, the horizontal size of this topographic feature is only twice as much as the crustal thickness. It is rather unlikely that exhumation produced by lithospheric folding (Cloetingh and Burov 2011;Cloetingh et al. 1999) or dynamic topography caused by mantle upwelling (Braun et al. 2013) may generate a distinct differential exhumation on such a small scale. The observed asymmetry of measured AFT ages on both sides of the Weinsberg Forest is better explained by the development of an ancient escarpment from the latest Cretaceous (Maastrichtian) to the Eocene. During this time span of ca. 35 Ma, the southern Bohemian Massif and the adjacent shelf area formed a passive continental margin at the northeastern side of the North Penninic Ocean.
Prominent examples of escarpments on passive continental margins are those of the Drakensberg in South Africa (cf. Brown et al. 2002), the Transantarctic Mountains (cf. Fitzgerald and Gleadow 1990) and the Snowy Mountains in New South Wales (Australia; cf. Kohn et al. 1999). Their surface uplift, denudation rates, drainage evolution and erosional retreat are produced by fluvial incisions from the coast backward and by progressive isostatic rebound of the crust underneath the eroded coastal slopes. This process creates an asymmetric exhumation pattern at both sides of the crest. By means of backstacking AFT topography, Brown et al. (2002) have estimated 4.5 km of denudation for the coastal zone adjacent to the Drakensberg Escarpment in South Africa, which rises up to a present-day altitude of 3482 m. Maximum altitudes of the Weinsberg Forest are much smaller (up to 1053 m). But a cumulative exhumation of about 2000 m (including basement erosion and removal of a putative overburden of Upper Cretaceous sediments) since the latest Cretaceous would have been sufficient to bring rocks from a Late Cretaceous PAZ to the surface. With an assumed geothermal gradient of 30 °C/km, latest Cretaceous and Cenozoic erosion could have produced a regional cooling in the order of 60 °C on the sloping area to the SW of the crest. In addition to a Neogene surface temperature of ca. 15 °C, this cooling would correspond to initial rock temperatures of ca. 75 °C previous to latest Cretaceous and Cenozoic exhumation.
We propose the name of 'Weinsberg escarpment' for the Paleogene geomorphic situation on the southwestern slopes of the Weinsberg Forest. The term 'Danube escarpment' would be misleading, because the river Danube did not yet exist in the Paleogene and the original escarpment was not shaped by the present-day fluvial catchment.
After the closure of the North Penninic Ocean, during the Late Eocene, either collisional coupling of the advancing Alpine-Carpathian nappe stack with its northern foreland (Ziegler and Dèzes 2007) or-more likely-East-Alpine slab detachment (Handy et al. 2021) may have caused accelerated uplift, exhumation and cooling during the latest 20 Ma.

Conclusions
The regional distribution of AFT ages, the synopsis of already known paleogeographical reconstructions, as well as inverse modeling of both AFT and concordant AHe data lead to the following conclusions: and Vitis faults. Important vertical offset along these faults during the Cretaceous and Cenozoic is thus precluded. 3. Shallow-marine and terrestrial sediments have reburied large areas of the southern Bohemian Massif during the Late Cretaceous, until the Campanian. This cover protected a former weathering basal relief from later erosion. Morphological remnants of this Mesozoic weathering basal relief are preserved as 'sealed relief' until present. The Late Cretaceous reburial and reheating caused broad and even bimodal frequency distributions of confined track lengths in apatite, but did not produce a full reset of the AFT and AHe systems. 4. The regional asymmetry of measured AFT ages relative to the NW-SE trending ridge of the Weinsberg Forest is not compatible with the geomorphological model of a piedmont staircase, which would require uplift en bloc, similar exhumation rates on both sides of the ridge, and similar cooling ages at similar altitude on both sides of the ridge. 5. Large-scale domal uplift of the Variscan basement may have produced Early Cretaceous exhumation but is ruled out for Cenozoic time because of the quoted regional asymmetry around the highest summits of the Weinsberg Forest. 6. The sloping area between the Weinsberg Forest and the Danube started to develop as an escarpment on a passive continental margin adjacent to the North Penninic Ocean. In latest Cretaceous (Maastrichtian) and Paleogene times, this 'Weinsberg escarpment' formed the transition between the shelf area and the elevated peneplain on top of the Weinsberg Forest. Present-day surface rocks of this sloping area did already reside in the apatite PAZ before their final Cenozoic exhumation. 7. Either collisional coupling of the East-Alpine and Carpathian orogens with their northern foreland or slab detachment underneath the Eastern Alps caused a final accelerated exhumation that started ca. 20 Ma agoespecially in the sloping area to the NE of the Danube.