1 Introduction

The Western Carpathians in central Europe represent the northernmost branch of the Alps and form the westernmost segment of a curved Carpathian orogenic belt (Fig. 1a). Similarly to other sectors of the Alpine belt, the Western Carpathians have experienced a rather complex Alpine tectonothermal history, comprising Jurassic rifting and basin formation, Cretaceous collisional tectonics, extensional collapse and lateral escape of fragments of the Adriatic (Apulian) plate and their complex interaction with the European foreland in the Tertiary (e.g. Royden et al. 1982; Mahel’ 1986; Csontos 1995; Ratschbacher et al. 1991; Tari et al. 1992; Plašienka 1996, 1997; Frisch et al. 2000; Sperner et al. 2002). Although orogen dynamics in the case of the Western Carpathians is well described by conceptual models (see e.g. Kováč et al. 1994; Plašienka et al. 1997a), important aspects, such as the timing and grade of metamorphism or the timing and rate of exhumation and topography formation, are little understood or controversial. This is due to lack of empirical data and a sparse geological record.

Fig. 1
figure 1

a Tectonic sketch map of the Western Carpathians with exposures of Variscan crystalline complexes belonging to three principal units and occurrences of Paleogene sediments and Neogene volcanic rocks. Inset map shows the location a; WC Western Carpathians, EA Eastern Alps. b Digital elevation model of the Nízke Tatry Mts. (NT) with position of the samples and trace of the profile presented in the Fig. 2. c Geological sketch map of the NT and surrounding areas (modified after Biely et al. 1992; Lexa et al. 2000) with location of the samples and AFT ages reported by Král’ (1977). d Schematic profile after Biely et al. (1992) (for location, see dashed line in b)

In this study, we focus on the Nízke Tatry Mts. (NT; Lower Tatra Mountains in English) in the central part of the Western Carpathians (Fig. 1a), a key area for understanding the Alpine history of the region. Tectonothermal evolution of the NT is believed to be well understood (e.g. Kováč et al. 1994), but, this understanding is not supported by empirical data. Furthermore, inferring from our experiences in other mountain ranges in the Western Carpathians we suspect that some classic models and long accepted definitions might be misconceptions:

  1. 1.

    Similar to other ‘core’ mountains in the Western Carpathians (i.e. basement-cored mountains belonging to the Tatric superunit sensu Mahel’ 1986 and Plašienka et al. 1997a), the core of the NT is comprised of Variscan crystalline basement consisting of Variscan granitoids and metamorphic rocks. In Carpathian literature, the Tatric Variscan basement complexes including the NT are traditionally defined as completely lacking or being only weakly overprinted by an Alpine metamorphic event. Traditionally, the Alpine overprint is thought to be restricted to narrow shear zones, reaching P–T conditions of the anchizone or lower greenschist facies (Bujnovský and Lukáčik 1985; Krist et al. 1992; Plašienka et al. 1997a; Madarás et al. 1996; Plašienka 2003a, b). However, based on petrographical and thermochronological evidence, this widely accepted definition was challenged by some authors who argued that some Variscan granitic cores did experience very low-grade (anchizonal) Alpine metamorphic overprint (Faryad and Dianiška 2002; Danišík et al. 2008a, 2010).

  2. 2.

    The generally accepted model of exhumation for the Western Carpathians was presented by Kováč et al. (1994). Based on AFT data, the authors proposed a progressive exhumation of crystalline bodies, migrating from internal parts of the belt towards the orogenic front, with the NT being exhumed in the Eocene (37–52 Ma). Recently, work on other ‘core’ mountains has implied that this model is imprecise (Danišík et al. 2004, 2008a, 2010), however, no new data from the NT have as yet been reported.

  3. 3.

    The NT form an E–W oriented ridge, about 70 km long, 20–30 km wide, >2,000 m high (Fig. 1b), which is the second highest ridge in the Western Carpathians. It represents the most important watershed in the Western Carpathians, but its topographic evolution is poorly understood. In the Pleistocene, the topography of the NT was probably fairly similar to that of the present, dictating the pattern of local mountain glaciation (Louček et al. 1960), but the question of when the present-day high relief formed has not been addressed.

To tackle these issues, samples from a transect across the NT ridge were analyzed using a combination of four thermochronometers: zircon fission track (ZFT), zircon (U–Th)/He (ZHe), apatite fission track (AFT) and apatite (U–Th–[Sm])/He (AHe) dating methods. Sensitivities of these methods cover the temperature range of ~270–40°C (Hurford 1986; Wagner and Van den haute 1992; Brandon et al. 1998; Wolf et al. 1998; Carlson et al. 1999; Farley 2000). The multi-dating approach facilitates a reconstruction of the full thermal evolution at mid- to shallow crustal levels and addresses the issues of metamorphism, burial and exhumation, and topographic evolution raised above.

2 Geological setting

In this study we focus on the western part of the NT, which is dominated by Variscan crystalline basement with sparsely preserved Mesozoic sedimentary cover and two Mesozoic superficial nappes (Fatric and Hronic, respectively) flanking the basement on the north and south (Fig. 1c, d).

The basement consists of Variscan granitoids (zircon U/Pb ages: 343 ± 4 and 330 ± 10 Ma; Poller et al. 2001; Putiš et al. 2003) and Variscan metamorphic rocks (mostly amphibolite-facies paragneiss, orthogneiss, amphibolite and migmatite; U/Pb ages ~360 Ma; e.g. Putiš et al. 2009a). Variscan cooling of the basement is constrained by K/Ar and Ar/Ar dating of muscovite and biotite (ages ~330–300 Ma; Kantor et al. 1984; Maluski et al. 1993; Dallmeyer et al. 1996). The studied area is classified in the Carpathian literature as a part of the Tatric superunit (or Tatricum; Mahel’ 1986; Plašienka et al. 1997a) and according to the traditional definition, as the part of the Tatricum, the basement should lack an Alpine metamorphic overprint (Mahel’ 1986; Plašienka et al. 1997a). To the southeast, the Tatric basement is bordered by a distinct thrust zone (Čertovica line) composed of mylonitic and phyllonitic rocks (Biely and Fusán 1969). This thrust zone separates Tatric basement and cover in the footwall from a hanging wall basement unit (Veporicum) in the southeast and is considered one of the most important Alpine shortening zones in the Western Carpathians (Biely and Fusán 1969; Plašienka et al. 1997b). Although the age of activity of the Čertovica line has not been radiometrically determined, it is accepted that thrusting of basement units and emplacement of superficial nappes occurred during the Alpine (Eo-Alpine) collision between the European and the Adriatic plate from ~110 to ~85 Ma. This has been inferred from paleontological evidence from the nappe sequences, and both K/Ar and Ar/Ar data from similar shear zones (e.g. Dallmeyer et al. 1996; Plašienka 1997; Putiš et al. 2009a).

Post-thrusting evolution of the NT remains unclear because there are no post-Cretaceous sediments preserved directly on the basement of the NT, but it can be traced in the sedimentary record of the surrounding depressions (Fig. 1c, d). The first post-Cretaceous record is represented by deposits of the Central Carpathian Palaeogene Basin (CCPB) preserved north and south of the range (Fig. 1c, d; Gross 1978; Gross 2008; Gross et al. 1984). The CCPB sediments reach a thickness of up to 3.5–4 km in the deepest part of the basin (Soták 1998; Janočko and Jacko 2001) and consist of shallow to deep marine sediments of Lutetian to Aquitanian age (Gross et al. 1984; Gross 2008). On the northern slopes of the NT, thickness of these clastic and organo-detritic sediments reaches 500–600 m and increases towards the north, reaching a maximum of 1,600 m (Gross et al. 1980). South of the NT, a maximum of 550 m is reached (Filo in Bezák et al. 2009). It is noteworthy that clastic material from the NT crystalline basement was found in CCPB sequences on both sides of the NT mountains ridge. The position of the NT during the CCPB sedimentation is controversial. Gross et al. (1984) and Gross (2008) argued that the crystalline basement was eroded and has supplied material to the CCPB since the Late Eocene (~36 Ma), whereas Kázmér et al. (2003) argued for a burial of the NT by CCPB sediments from the Middle to Late Eocene (~40–38 Ma) until the Late Oligocene (~28 Ma).

Occurrences of Neogene sediments can be found only south of the range and are represented by Oligocene(?)–Lower Miocene conglomerates, consisting of pebbles of dolomites, limestones and the crystalline basement (Biely and Samuel 1982). It is widely accepted that during Neogene times, the NT were subjected to erosion and karstification but did not yet form a pronounced mountain range as there are no post-Lower Miocene sediments preserved in the area (Biely et al. 1992; Kováč 2000 and references therein). The mountain range must have existed in the Pleistocene as evidenced by the presence of glaciers at that time (Louček et al. 1960).

The exhumation history of the NT has been investigated by means of thermochronology by Král’ (1977). The author reported two AFT ages (52 ± 9 and 37 ± 5 Ma) measured by the now obsolete population dating method and argued for Eocene uplift although no track length data justifying this statement were provided. This conclusion was then widely accepted by the geological community and has served as the main support for the currently accepted exhumation model for the Western Carpathians (Kováč et al. 1994).

3 Samples and methods

Five samples of the Variscan granite were collected for thermochronological investigation from a tectonically undisturbed, ~8 km long profile, crossing the ridge of the NT through the highest peak (Figs. 1b, c, 2a). The maximum vertical difference between the samples is 887 m; sample NT-1 was taken from the tectonic contact of the basement and Triassic limestone, representing the basal part of the Fatric nappe thrust upon the crystalline basement in the Turonian (Biely et al. 1992); sample NT-4 was taken from the highest peak of the NT (Ďumbier, 2,045 m a.s.l.).

Fig. 2
figure 2

a Topographic profile with position of the samples and measured thermochronological data. AFT apatite fission track, AHe apatite (U–Th–[Sm])/He, ZFT zircon fission track, ZHe zircon (U–Th)/He. Large spread of AHe ages hampers any attempts of correlation with topography. b Apatite fission track length data. Explanation in histograms (from top): sample code; mean track length ± SD (both in μm); number of measured tracks. c All data displayed on a radial plot (Galbraith 1988; Vermeesch 2009), clearly showing minimum difference between AFT and ZHe ages, pointing to a fast cooling in the Eocene. ZFT and AHe form well separated clusters related to different thermal events

Apatite and zircon grains were separated using conventional magnetic and heavy liquid separation techniques. Fission track analysis was carried out using standard procedures described in Danišík et al. (2007). The external detector method (Gleadow 1981) was applied with the etching protocols of Donelick et al. (1999) for apatite (5.5 M HNO3 for 20 s at 21°C) and Zaun and Wagner (1985) for zircon (eutectic mixture of KOH and NaOH at 215°C for 12 h). The zeta calibration approach (Hurford and Green 1983) was adopted to determine the ages. FT ages were calculated using TrackKey 4.2 g (Dunkl 2002). In apatite, horizontal confined tracks in tracks were measured in c-axis parallel surfaces and were normalized for crystallographic angle using a c-axis projection (Donelick et al. 1999; Ketcham et al. 2007a). The annealing properties of apatite were assessed by measuring Dpar values (Dpar, the mean etch pit diameter of fission tracks measured on the apatite polished surface parallel to the crystallographic c-axis; e.g. Burtner et al. 1994).

For (U–Th–[Sm])/He analysis, apatite and zircon crystals were hand-picked following strict selection criteria (Farley 2002; Reiners 2005), then photographed and measured. Apatite was loaded in Pt tubes, degassed at ~960°C under vacuum using laser-heating or furnace and analyzed for 4He using a Pfeiffer Prisma QMS-200 mass spectrometer. Following He measurements, the apatite was spiked with 233U and 230Th, dissolved in nitric acid and analyzed by isotope dilution inductively coupled mass spectrometry for U, Th and in one batch also for Sm on a Perkin Elmer (ELAN DRC II) ICP-MS. Zircon was loaded in Nb tubes, degassed at ~1,250°C and analyzed for 4He using the same facility as for apatite. Degassed zircon was dissolved following the procedure of Evans et al. (2005) and analyzed by isotope dilution inductively coupled mass spectrometry for U and Th on an Agilent 7500 ICP-MS. For more details on analytical procedures, the reader is referred to Evans et al. (2005) and Danišík et al. (2008c).

Total analytical uncertainty was calculated as the square root of the sum of the squares of weighted uncertainties on U, Th, Sm and He measurements. Total analytical uncertainty was less than ~5% in all cases and was used to calculate the error of raw (U–Th–[Sm])/He ages. The raw AHe and ZHe ages were corrected for alpha ejection (Ft correction) after Farley et al. (1996) and Hourigan et al. (2005), respectively. A value of 5% was adopted as the uncertainty on the Ft correction, and was used to calculate errors for the corrected AHe and ZHe ages.

The low-temperature thermal history based on fission track and (U–Th–[Sm])/He data was modelled using the HeFTy modelling program (Ketcham 2005) operated with the multikinetic fission track annealing model of Ketcham et al. (2007b) and the diffusion kinetics of the Durango apatite after Farley (2000) and zircon after Reiners et al. (2004). Inferring from fission track mounts, we assumed homogeneous distribution of U and Th in the apatite and zircon.

4 Results

The results of thermochronological analyses are summarized in Tables 1 and 2 and shown in Fig. 2. All five samples yielded fairly consistent results, clearly defining distinct age groups. Higher temperature chronometers revealed older ages than chronometers with lower temperature sensitivity, which is in agreement with the closure temperature concept (Dodson 1973).

Table 1 Fission track results
Table 2 (U–Th–[Sm])/He results

ZFT age spectra are fairly broad, but passed the Chi-square test and thus represent a single age population. The ZFT ages form a loose cluster of Late Jurassic–Early Cretaceous age (132.1 ± 8.3; 155.1 ± 12.9; 146.8 ± 8.6; 144.9 ± 11.0 Ma).

ZHe data ages reproduce well: in most of the cases four out of five replicates reproduce within one sigma error and one replicate is a ‘flier’ considerably older than the rest. Single grain ZHe ages of reproducible replicates form a tight cluster of Early–Middle Eocene age ranging from 53.4 ± 3.8 to 39.4 ± 2.6 Ma (15 replicates). The ‘fliers’ range from 77.8 ± 5.2 to 58.3 ± 4.1 Ma (5 replicates). AFT ages are slightly younger than ZHe ages and form a distinct cluster of Middle Eocene age (40.3 ± 1.7; 42.1 ± 1.9; 43.9 ± 2.4; 41.5 ± 2.8; 40.0 ± 2.3 Ma; Fig. 2a, c), indicating rapid cooling through the zircon ZHe partial retention zone (~200–160°C; Reiners et al. 2004) and apatite partial annealing zone (~120–60°C; Wagner and Van den haute 1992) during that period. All AFT samples passed the Chi-square test and are therefore considered to form one age population. The average Dpar value for all samples is ~1.6 μm, indicating fluorine-rich apatite, characterized by relatively low annealing temperature (~60–120°C; e.g. Wagner and Van den haute 1992; Ketcham et al. 1999). Four of five track length distributions are bimodal, with short mean track lengths (MTL 12.0–13.1 μm) and relatively large standard deviations (Fig. 2b), pointing to a complex thermal evolution with a long residence or a phase of reheating within the apatite partial annealing zone (Gleadow et al. 1986a, b).

Due to the unacceptable quality of apatite grains for purposes of (U–Th–[Sm])/He dating and frequent occurrence of fluid and mineral inclusions, we could not always measure five replicates per sample as in the case of zircon. The presence of microscopically undetected inclusions became obvious after the second gas extraction (so called ‘re-extract’). In such cases we did not proceed with U–Th measurements but rather tried to select another crystal. Despite significant effort, in only one sample (NT-2) were five replicates successfully measured, whereas the rest of the samples were run in duplicate and triplicate. Nevertheless, AHe ages are consistent, with the majority clustering around the Early Miocene. This suggests a Miocene residence of the samples within AHe partial retention zone and corroborates the assumption of complex cooling or reheating inferred from the track length distributions.

We found no obvious correlation between thermochronological data and sample elevation (Fig. 2a), suggesting that all samples experienced a fairly similar cooling history.

5 Interpretation and discussion

5.1 Alpine metamorphism

The Late Jurassic–Early Cretaceous ZFT ages show that the NT basement must have reached temperatures >210°C during post-Variscan times (assuming an effective closure temperature of 240 ± 30°C; Brandon et al. 1998) as the ‘Variscan memory’ in the ZFT system was fully reset. The maximum temperature must have been lower than ~300°C, since the Ar/Ar system retains Variscan ages (Maluski et al. 1993; closure temperature of Ar/Ar system in biotite is ~300°C; Harrison et al. 1985; McDougall and Harrison 1988). Given these constraints, the basement must have experienced a very low-grade or low-grade Alpine metamorphic overprint, in contrast to the traditionally accepted notion of no overprint (see Sect. 2; Mahel’ 1986; Plašienka et al. 1997a). These data thus clearly prove a very low-grade (anchizonal) Alpine overprint in the Tatric superunit (Faryad and Dianiška 2002; Danišík et al. 2008a, 2010).

Although the presence of Alpine metamorphism is documented by the ZFT data, its timing and geodynamic context remains uncertain. This is because no track lengths in zircons were measured and the thermal history could not be modelled, thus the meaning of the data is less clear. From the perspective of regional geology it is important to note that similar Jurassic to Early Cretaceous ZFT ages have been found in other parts of the Western Carpathians (Kováč et al. 1994; Plašienka et al. 2007). Furthermore, Danišík et al. (2010) reported almost identical Late Jurassic–Early Cretaceous ZFT ages from an adjacent Tatric crystalline basement (Malá Fatra Mts.) that should have had a common evolution to the NT throughout the Mesozoic (Plašienka et al. 1997a). The authors provide an extensive discussion of the possible meaning of these ZFT ages and we speculate that the ZFT ages are either (1) cooling ages recording a Jurassic/Cretaceous thermal event, perhaps related to the rifting (Plašienka 2003b; Putiš et al. 2009b) associated with increased heat flow and magmatic activity (Hovorka and Spišiak 1988; Spišiak and Hovorka 1997; Spišiak and Balogh 2002) or (2) apparent ages resulting from partial rejuvenation of the ZFT system during the Eo-Alpine collision (~110–85 Ma; Danišík et al. 2010).

5.2 Eocene exhumation

A tight cluster of Early–Middle Eocene ZHe ages and Middle Eocene AFT ages, backed up by the thermal modelling results (see Sect. 5.3) give strong evidence for a rapid cooling event between ~50 and ~40 Ma, when the NT basement cooled from ZHe partial retention zone to apatite partial annealing zone. Exhumation rates between 50 and 40 Ma determined by the mineral pair method and by conversion of modelled trajectories are in the range of ~250–900 m/Ma (assuming a geothermal gradient of 30°C/km; closure temperatures of ~190 and ~110°C for ZHe and AFT system, respectively; Reiners et al. 2004; Wagner and Van den haute 1992).

We interpret the cooling event as being related to collapse of the Carpathian orogenic wedge and explain it as follows: during the Eo-Alpine collision in the mid-Cretaceous, the Tatric superunit was overthrust by superficial nappes and its interior, including the NT crystalline core, was perhaps partly overridden by the Veporic basement/cover complex. This compressional event is well documented in the stratigraphy and by Ar/Ar ages from shear zones ranging between ~100 and 85 Ma (Dallmeyer et al. 1996; Putiš et al. 2009a). As a result of thrusting, the NT basement must have been buried to depths where the temperature was elevated enough to fully reset the ZHe system (>190°C) and even partially reset the ZFT system (>210°C).

The Eocene cooling so well documented by ZHe and AFT data is most likely related to exhumation of the NT basement, pointing to an extensional collapse of the Carpathian orogenic wedge during the Eocene. Although predominantly erosion-driven exhumation can be inferred from extensive occurrences of upper Middle to Upper Eocene detritus from the NT in the surrounding depressions, relatively high cooling rates are difficult to explain solely by erosion. Thus, we speculate that the exhumation could have an important tectonic component. Several studies have mentioned that in the Tertiary, the Čertovica line was reactivated as a transtensional or extensional shear zone (e.g. Šefara et al. 1998; Plašienka 2003a; Putiš et al. 2009a). It is likely that this process tectonically exhumed the footwall Tatric basement from below the hanging wall Veporic basement sheet.

5.3 Post-Eocene thermal event?

The post-tectonic evolution of the NT is difficult to track from the limited geological record, so the Miocene AHe ages may represent an important contribution to this topic. In order to elucidate the meaning of the AHe data, we first review relevant facts and currently accepted models, which we then test by thermal modelling, before presenting our conclusion.

Bartonian to Rupelian flysch sediments of the CCPB basin gave rise to two different models for the Palaeogene period (for current distribution of CCPB sediments see Fig. 1a):

  1. 1.

    Gross et al. (1984) and Gross (1978, 2008) suggest that already in the early Middle Eocene (Lutetian), prior to sedimentation of CCPB, the emerged NT formed a topographic high which was prone to erosion and was never covered by Palaeogene sediments. The first CCPB sediments of Bartonian age were derived from the Mesozoic cover nappes and Tatric sedimentary cover. Crystalline basement was first eroded and supplied material to the CCPB at the Bartonian/Priabonian boundary (~37 Ma), as documented by submarine fan deposits with crystalline clasts preserved north of the range. Material supply from the NT, however, quickly ceased by the Priabonian, as inferred from transport direction as well as lithological data. Nevertheless, the NT should have formed an island in the Palaeogene sea during Priabonian–Rupelian times (Gross 2008).

  2. 2.

    Kázmér et al. (2003), in contrast, proposed a large scale model for the Eastern Alps and Western Carpathians, whereby, from Bartonian to Egerian times (~40–28 Ma), the NT were completely buried by a thick pile of CCPB sediments.

Neogene tectonothermal history is even more difficult to track. Given the lack of sediments, it is widely accepted that the NT were not affected by Neogene transgression (Kováč 2000), so no reheating induced by burial should be expected. On the other hand, several studies have shown that the Tatric crystalline complexes were affected by a Miocene regional thermal event related to increased heat flow and volcanic activity in the Western Carpathian region (Danišík et al. 2008a, b; for current distribution of Neogene volcanic rocks see Fig. 1a).

In order to test these hypotheses, we used the HeFTy modelling program (Ketcham 2005) to calculate thermal histories that reconcile ZHe, AFT and AHe data. According to the hypotheses presented above, we tested four plausible cooling scenarios: simple cooling in the Eocene (after Gross 2008), Palaeogene burial (after Kázmér et al. 2003), Miocene reheating (after Danišík et al. 2008a, b), and superimposed Palaeogene burial and Miocene reheating (after Kázmér et al. 2003; Danišík et al. 2008a, b combined).

We found that the thermal trajectories, which best reconcile the measured data as given by the highest ‘goodness of fit’ values, are characterized by a fast cooling from >190 to <60°C at ~55–40 Ma, and require a reheating to ~55–90°C before cooling to the present-day temperature conditions (Fig. 3a–d). While the interpretation of the fast cooling is straight forward in terms of post-collisional collapse and exhumation, the meaning of the reheating episode is ambiguous. Unfortunately, the modelling could not resolve exactly when the heating pulse occurred and the suggested time span of 30–5 Ma is too broad to favour any of the above-mentioned scenarios. The only conclusion which can be drawn with confidence from the modelling results is that the NT must have experienced a thermal event between ~30 and 5 Ma with minimum peak temperatures of ~55–90°C. This event could be related either to the CCPB burial in the Oligocene or magmatic activity in the Miocene or Miocene sedimentary burial, where all options are equally plausible given the lack of geological constraints and considering the modelling results only. Nevertheless, inferring from the Early Miocene AHe age cluster, we are inclined to prefer a Miocene thermal event, related to magmatic activity and/or sedimentary burial.

Fig. 3
figure 3

Thermal modelling results of thermochronological data displayed in a time–temperature diagram modelled with the HeFTy program (Ketcham, 2005). The best fit is shown as a solid black line, shaded polygons (yellow acceptable fit, blue good fit) show the values of peak temperatures after Eocene cooling and during post-Eocene reheating. ZHePRZ ZHe partial retention zone, APAZ apatite partial annealing zone, AHePRZ AHe partial retention zone, MTL mean track length in μm, SD standard deviation in μm, GOF goodness of fit (statistical comparison of the measured input data and modelled output data, where a “good” result corresponds to value 0.5 or higher, “best” result corresponds to value 1). Note that according to the modelling results, reheating could occur anytime from the Oligocene to the Late Miocene and thus is not in favour of any of the suggested geological scenarios

5.4 Failed approach to palaeo-topography reconstruction

One of the objectives of this study was to obtain a better understanding of the topographic evolution of the NT. Our strategy was based on AHe dating of a ridge-valley profile crossing the NT ridge in the highest point. According to theory, the spatial pattern of AHe ages should reflect the amount of erosion or the shape of the subsurface isotherms, thereby providing constraints on palaeotopography and its change through time (Reiners 2007).

Despite demonstrated feasibility by several studies (e.g. House et al. 1998, 2001; Foeken et al. 2007), in our case this approach failed for two reasons: Firstly, the key requirement of this approach is that AHe ages are cooling ages, which allow straightforward calculation of erosion. Here, however, as proven by modelling results, AHe ages are not cooling ages related to a distinct cooling event, but are apparent ages resulting from a complex thermal history with a phase of reheating. In addition, we do not even know when this reheating occurred or its duration. Secondly, the resolution of the AHe dataset is too low (or the spread of single grain AHe ages is too high) to be correlated with present day topography (Fig. 2c). This is probably a consequence of a combination of several factors such as complex thermal history, quality of dated apatite grains, less pronounced relief (with too short wavelength and too low amplitude), or slow exhumation rates during the post-collisional stage. Therefore, we cannot draw any conclusion on the palaeotopographic evolution of the NT from our AHe data.

6 Conclusions

New ZFT, ZHe, AFT and AHe data enabled us to constrain the thermotectonic evolution of the NT and provide important constraints on the thermal and geodynamic evolution of the Western Carpathians. The most important results can be summarized as follows:

  • The Variscan crystalline basement of the NT was heated to temperatures above ~210°C and experienced a very low-grade to low-grade Alpine metamorphic overprint as documented by the Late Jurassic–Early Cretaceous ZFT ages ranging from 155.1 ± 12.9 to 132.1 ± 8.3 Ma. The timing and source of the heating remain unclear. We propose two alternative explanations: (1) either the ZFT ages are cooling ages related to a Jurassic/Cretaceous thermal event or (2) ZFT ages are apparent ages resulting from a partially reset during Eo-Alpine collision and thrusting in the mid-Cretaceous. Despite this ambiguity, ZFT data clearly disprove the widely accepted notion of non-existent Alpine metamorphism in the crystalline core complexes of the Tatric superunit.

  • ZHe and AFT data constrain a distinct fast cooling event in the Eocene between ~55 and 40 Ma, which we interpret in terms of erosional and tectonic exhumation of the basement, related to the collapse of the Carpathian orogenic wedge. This is the first thermochronological evidence univocally documenting cooling of the Tatric crystalline core in the Eocene.

  • AHe data revealed a thermal event in the Oligocene–Miocene, when the crystalline basement was heated to temperatures of ~55–90°C. Given the lack of geological information and insufficient resolution of the modelling results, this thermal event may be related to Oligocene/Miocene sedimentary burial, Miocene mantle upwelling, magmatic activity and/or increased heat flow in the Carpathian realm. Regardless of the cause, this finding disproves the widely accepted concept of thermal stability of the NT during the post-Eocene period.

  • Complex thermal evolution resulting in apparent AHe ages and limited resolution of thermochronological data did not allow us to place solid constraints on the topographic evolution of the NT.