doi: 10.1111/ter.12221 Late Miocene increasing exhumation rates in the eastern part of the Alps – implications from low temperature thermochronology
Andreas W€olfler,1 Walter Kurz,2 Harald Fritz,2 Christoph Glotzbach1 and Martin Dani s ık3 1Institut fur€ Geologie, Leibniz Universitat€ Hannover, Callinstraße 30, Hannover D-30167, Germany; 2Institut fur€ Erdwissenschaften, Karl-Franzens Universitat€ Graz, Graz A-8010, Austria; 3John de Laeter Centre for Isotope Research, Department of Applied Geology, The Institute for Geoscience Research (TIGeR), Curtin University, GPO Box U1987, Perth, WA 6845, Australia
ABSTRACT A new set of apatite fission-track and apatite (U–Th)/He data operated simultaneously during lateral extrusion of the East- reveals a hitherto undated late Miocene exhumation pulse in ern Alps. As the higher late Miocene/Pliocene exhumation the eastern part of the Eastern Alps. While distinct parts of rates are restricted to a single tectonic block, namely the Nie- the study area, including the Seckauer Tauern, have been at dere Tauern, we infer a tectonic trigger that is probably near surface conditions (<100 °C) since the Eocene, the neigh- related to a change in the external stress field that affected bouring Niedere Tauern experienced enhanced cooling and the Alps during this time. exhumation in the middle Miocene and again at the late Mio- cene/Pliocene boundary. Middle Miocene exhumation is inter- Terra Nova, 00: 1–9, 2016 preted as a result of tectonic escape and convergence that
thermochronological point of view, and perturbation of isotherms, which Introduction this part of the Alps is less well may result in false conclusions about The existing low-temperature ther- investigated than other areas such as exhumation rates in the upper crust mochronological datasets from the the Tauern Window or the Western (Stuwe€ et al., 1994; Braun, 2002). To European Alps have been success- Alps (e.g. Foeken et al., 2007; Luth estimate the influence of topography fully used to derive an orogen-wide and Willingshofer, 2008; Vernon on perturbation of isotherms and erosional history, which indicates an et al., 2008). In particular, AHe data, thermochronological data, we correct increase in erosion rates during the which have the lowest closure tem- the derived exhumation rates follow- late Miocene and Pliocene (Vernon perature (~60 °C) among all ther- ing the recently proposed approach of et al., 2008; Herman et al., 2013). mochronometers, are rare, and thus Glotzbach et al. (2015). In addition, Although the findings from ther- the late stages of final exhumation we use conventional thermal-history mochronological data agree with from shallow crustal depths with modelling to constrain the near sur- those from studies based on indepen- temperatures <100 °C are not well face cooling history of the study area. dent approaches such as sediment- constrained. Accordingly, this region budget analysis or cosmogenic may bear an undiscovered thermal Geological setting nuclides (Kuhlemann et al., 2002; signal that can provide valuable Wittmann et al., 2007; Wagner et al., information about the late Neogene The European Alps (Figs 1a,b) are 2010; Legrain et al., 2014a,b), the exhumation history of the Eastern the result of the convergence between discussion about possible drivers (cli- Alps. the European and Adriatic plates. In mate change vs. tectonics) of the Low-temperature thermochronol- the course of their collision during late-stage uplift is ongoing (Ceder- ogy has been proven to be a powerful the Palaeogene, the Penninic domain bom et al., 2004, 2011; Willett, 2010; tool to infer the exhumation history was overridden by the Austroalpine Herman et al., 2013; Baran et al., of rocks in qualitative as well as quan- nappes, formerly part of the Adriatic 2014). titative ways (Reiners and Ehlers, plate (Schmid et al., 2004). Subse- To contribute to this debate, we 2005). However, thermochronology quent crustal-scale folding, orogen- provide new apatite fission-track has distinct limitations for estimating parallel extension and lateral extru- (AFT) and apatite (U–Th)/He (AHe) exhumation rates. In particular, the sion led to the formation of the data to infer exhumation rates from palaeo-geothermal gradient is usually Tauern Window (Fig. 1b; Ratschba- the eastern part of the Eastern poorly constrained so that the depth cher et al., 1991; Frisch et al., 1998; Alps (Figs 1 and 2). From the of the closure temperature can be Rosenberg et al., 2007). Although specified only vaguely. Therefore, we the timing of rapid cooling and collected samples from different ele- footwall exhumation in the Tauern € Correspondence: Andreas Wolfler, Institut vations to calculate exhumation rates Window is well constrained for the for Geology, Leibniz University of from the slope of the age–elevation early and middle Miocene (Fugen-€ Hannover, Callinstraße 30, Hannover relationship (AER). However, such schuh et al., 1997; Luth and Willing- D-30167, Germany. Tel.: +0049 511 762 interpretations are often hampered by shofer, 2008; Scharf et al., 2013; 3871; fax: +0049 511 762 2172; e-mail: inhomogeneities in the thermal field Favaro et al., 2015), the relative woelfl[email protected]
© 2016 John Wiley & Sons Ltd 1 Early to Late Miocene exhumation rates in the Eastern Alps • A. W€olfler et al. Terra Nova, Vol 0, No. 0, 1–9 ......
17° 16° 15° (a)
17° 5° 7° 9° 11° 13° 15° Neogene and Quarternyry sediments Donau WienWien München 48° 48° Nappes derived from the Tethys Wien Zürich Adria derived units Maribor 48° 46° 46° Penninic nappes Milano Europe derived units 17° 11° 13° 15° 44°
(b) SEMP
Tauern Window Graz 47° 47° BF KF DAV Pannonian MV Fig. 2 basin PF MariborDrava
46° 17° 16° Adamello 15° 46° 14° 10° 11° 12° 13° Neogene sedimentary and magmatic rocks Upper Austroalpine nappe system: Penninic nappe system: European derived units
Upper Cretaceous to Paleogen sediments Drauzug-Gurktal and Ötztal- Piemontais Helvetic and of the Gosau Group Bundschuh nappe system Ultrahelvetic nappes Units derived from the Meliata-Vardar domain Koralpe-Wölz nappe system Brianconnais undeformed European continent Adriatic plate Silvretta-Seckau nappe system Valais
Dinarides Subpenninic nappes
Southalpine unit Lower Austroalpine subunit
Fig. 1 (a) Simplified tectonic map of the Alps. (b) Tectonic map of the Eastern Alps following the nomenclature of Schmid et al. (2004) (modified after Froitzheim et al., 2008). BF, Brenner normal fault; KF, Katschberg normal fault; DAV, Defer- eggen–Antholz–Vals fault; SEMP, Salzach–Ennstal–Mariazell–Puchberg fault; PF, Pustertal fault as part of the Periadriatic fault system; MV, Moll€ valley fault. contributions of orogen-parallel are sparse; the existing data indicate kilometers during the whole Ceno- extension by normal faulting, strike- segmentation of the Austroalpine zoic. The Niedere and Seckauer slip faulting and compression are still upper crust into tectonic blocks with Tauern are bordered by the Salzach- a matter of discussion (Lammerer, different cooling histories (Frisch tal–Mariazell–Puchberg fault 1988; Fugenschuh€ et al., 1997, 2012; et al., 1998; Fig. 2). The Niedere (SEMP), the Paltental-Liesingtal Kuhlemann et al., 2001; Linzer et al., Tauern revealed middle Miocene fault (PLF), the Pols-Lavanttal€ fault 2002; Rosenberg and Garcia, 2011). AFT ages (Hejl, 1997; Reinecker, (PoLF),€ the Mur-Murz€ fault (MMZ) Our study area is situated within 2000) documenting cooling that is and the Niedere Tauern South fault the Austroalpine nappe complex to contemporaneous with the cooling of system (NTSFS; Fig. 2). the east of the Tauern Window Penninic units of the eastern Tauern The Niedere Tauern and the (Fig. 1b), in the so-called Niedere Window (Reinecker, 2000; Wolfler€ Northern Calcareous Alps are char- and Seckauer Tauern (Fig. 2). It et al., 2011). The Austroalpine units acterised by high relief with steep comprises polymetamorphic rocks to the south and east of the Niedere slopes and an average elevation of that experienced the most recent Tauern show Cretaceous zircon fis- 1.67 km; on the other hand, the metamorphism during the late Creta- sion track ages (Kurz et al., 2011), Seckauer Tauern and the Sau and ceous (Schmid et al., 2004) and sub- Eocene to Oligocene AFT ages (Hejl, Koralpe are characterised by lower sequent cooling, as indicated by 1997; Reinecker, 2000; Wolfler€ et al., relief and a distinctly smoother muscovite K–Ar ages of 95–70 Ma 2010) and Oligocene to early Mio- topography with an average eleva- (Frank et al., 1987). Low-tempera- cene AHe ages (Legrain et al., tion between 1.44 and 1.12 km ture thermochronological data from 2014a), showing that this region (Figs 1 and 2; Frisch et al., 2000a). the eastern part of the Eastern Alps resided at a depth of only a few In the late Oligocene and early
2 © 2016 John Wiley & Sons Ltd Terra Nova, Vol 0, No. 0, 1–9 A. W€olfler et al. • Early to Late Miocene exhumation rates in the Eastern Alps ......
PLF Northern Calcaraous Alps
SEMP Niedere Tauern Seckauer Tauern Ju9: 48.5±7.4 16.0±1.6 19.4±1.5 43.7±2.5 38.6±1.4 14.4±1.1 15.4±1.3 Ju10: 46.3±4.0 17.8±2.6 Ju7: 44.4±6.4 So9: 16.9±2.6 19.6±1.6 PöLF Ju8: 42.4±4.2 6.6±0.5 22.9±4.0 23.8±3.5 23.8±3.5 So7: 15.9±2.3 B Ju11: 40.9±5.9 6.3±0.3 So11: 15.9±2.5 38.9±5.2 15.3±2.0 A So13: 15.0±2.6 So8: 15.2±2.1 47.1±4.3 So14: 14.3±2.5 FoB 6.1±0.3 15.6±1.8 So15: 13.9±2.4 8.6±2.3 So6: 14.1±2.1 NTSFSTB So5: 14.1±1.0 Koralpe 5.7±0.9
Tauern Window MMZ
30.9±4.0 GöF Saualpe 32.3±3.3 32.9±2.7 17.3±1.2 Gurktal Alps 23.0±4.2 35.1±3.1 29.0±2.4 KF Ju11: 40.9±2.3 AFT this study 43.9±9.5 AHe this study
12.9±1.4 43.7±2.5 AFT Hejl (1997)
22.9±4.0 AFT Reinecker (2000)
MV 8.6±2.3 AFT Bertrand (2013) 03.5 7 14 21 28
N Kilometers
Fig. 2 Digital elevation model of the study area with new and previously published AFT and AHe data. The Niedere and Seck- auer Tauern are bordered by the sinistral Salzachtal–Ennstal–Mariazell–Puchberg fault (SEMP) and the Paltental-Liesingtal fault (PLF) to the north and the sinistral Mur-Murz€ fault (MMZ) and the Niedere Tauern South fault system (NTSFS) to the south. The dextral Pols-Lavanttal€ fault (PoLF)€ separates the Niedere and the Seckauer Tauern. These fault zones were active during Middle Miocene lateral extrusion (Ratschbacher et al., 1991). We note that the complex geometry and the possible tem- poral and structural interplay between the NTSF and the MMZ are still not fully understood and require further investigations (e.g. Reinecker, 2000; Wolfler€ et al., 2011). Go,€ Gortschitz€ fault; MMZ, Mur-Murz€ fault; MV, Moll€ valley fault; TB, Tamsweg basin; FoB, Fohnsdorf basin. AFT and AHe data from the eastern Tauern Window range from 25 to 3 Ma and from 15 to 8 Ma respectively (Foeken et al., 2007; Staufenberg, 1987; Wolfler€ et al., 2008, 2012; Bertrand, 2013).
Miocene, the study area was charac- enhanced block movement in the et al., 1998). The modern topogra- terised by lowlands and hilly areas course of lateral extrusion led to the phy of the eastern Tauern Window with a northward-directed drainage formation of intramontane basins and adjacent eastern Austroalpine system (Frisch et al., 1998). At the (Fig. 2) and differential exhumation units formed in the late Miocene, early/middle Miocene boundary, and uplift of crustal blocks (Frisch contemporaneous with a switch from
© 2016 John Wiley & Sons Ltd 3 Early to Late Miocene exhumation rates in the Eastern Alps • A. W€olfler et al. Terra Nova, Vol 0, No. 0, 1–9 ...... aN–S- to an E–W-directed parameters can be varied to test, for are given in the supporting informa- drainage system (Frisch et al., 1998; instance, possible variations in the tion. Foeken et al., 2007; Robl et al., palaeo-topography by applying a ‘to- 2008, 2015). pographic amplification factor’. This Results and discussion factor linearly scales the present topography: a factor of 0.5 reduces Samples and methods Data and modelling results relief by 50% and a factor of 1.5 For AFT and AHe analysis, we col- increases relief by 50% (Glotzbach The samples from the Seckauer and lected six orthogneisses and eight et al., 2015). Niedere Tauern display Eocene and micaschists, which were treated using For thermal-history modelling, we middle to late Miocene AFT and standard magnetic and heavy liquid used the HeFTy software (Ketcham, AHe ages respectively (Fig. 2). The separation techniques. Details of 2005) with the annealing algorithm apparent exhumation rate deduced AFT and AHe analytical methods of Ketcham et al. (1999). More from the AER is 0.14 mm/yr for the are given in the supporting informa- details on the modelling approaches Seckauer Tauern (Fig. 3a), indicating tion. The slopes of the AER were corrected for the effect of isotherm 3000 perturbation caused by topography. (a) Seckauer T auern The age perturbations at sample 2500 localities were estimated by taking into account the dominant wave- AFT: 0.14 mm year–1 2000 lengths and amplitudes of the topog- R 2 = 0.98 raphy using a Fast Fourier 1500 Transform (FFT) algorithm (Glotz- Elevation (m) bach et al., 2015). The age perturba- 1000 tion was calculated assuming a steady-state topography, but also for 500 Cret. Palaeogene possible palaeo-topographies with Maas. Palaeocene Eocene Olig. higher or lower relief. In general, a higher relief results in a lower AER- 70 60 50 40 30 derived exhumation rate and vice Age (Ma) versa. The steps involved in the (b) Niedere Tauern: Profile A 3000 quantification of the age peraturba- tion are (Glotzbach et al., 2015): AFT: AHe: 2500 0.54 mm year–1 1.0 mm year–1 1 Decomposing topography (digital R 2 = 0.89 R 2 = 0.98 2000 elevation model (DEM) with 1 km
resolution) using a FFT algorithm 1500
along N–S, E–W, NW–SE and Elevation (m) – NE SW profiles. 1000 2 For each combination of wave- length, amplitude, exhumation rate Early M. Middle M. Late Miocene Pli. P. and geothermal gradient, the age 20 15 10 5 0 perturbation was previously esti- Age (Ma) mated with PECUBE (Braun et al., 2012). (c) Niedere Tauern: Profile B 3000 3 The age perturbation along each AFT: 2500 decomposed profile and Fourier 0.52 mm year–1 series is calculated based on the R 2 = 0.97 2000 PECUBE-derived age perturbation dataset. 1500
4 Combining the age perturbations Elevation (m) of all Fourier series along all pro- 1000 files yields an age perturbation map the same size as the input Early M. Middle M. Late Miocene Pli. P. DEM. 20 15 10 5 0 Age (Ma) The resulting age perturbation map can be used to estimate the site- Fig. 3 Age–elevation relationships of the data, and apparent exhumation rates specific potential perturbation of derived from them. Exhumation rates (in mm/yr) are calculated as error-weighted thermochronological ages caused by linear best fit (solid red lines) with 95% confidence level (dashed lines). Black topography and to correct ages, squares are AFT data; white circles are AHe data. Cret., Cretaceous; Maas., Maas- age–elevation profiles and the result- trichtian; Olig., Oligocene; Early M, Middle M., Early and Middle Miocene; ing exhumation rates. Input model Pli., Pliocene; P., Pleistocene.
4 © 2016 John Wiley & Sons Ltd Terra Nova, Vol 0, No. 0, 1–9 A. W€olfler et al. • Early to Late Miocene exhumation rates in the Eastern Alps ...... that the samples cooled slowly with contradiction probably arises from slowly through the APAZ and little exhumation and the ages have limitations caused by the unknown HePRZ and were at near surface probably not been affected by the geothermal gradient and/or the conditions by Eocene to Oligocene development of topography. Slow chemical composition. Although times (Table 2). cooling rates can also be inferred Dpar values are between 1.57 and The middle Miocene AFT ages of from the relatively small track length 1.58 lm, typical of a fluorine-apatite the Niedere Tauern yield apparent distribution and thermal-history composition, little is known about exhumation rates of 0.54 and modelling (Fig. 4c, Table 1). It the role of ions such as OH, Mn, Sr 0.52 mm a 1 (Fig. 3b,c). In contrast, should be noted that sample Ju11 or the REE content (Barbarand the AHe data yield an apparent yields AFT and AHe ages that over- et al., 2003; Spiegel et al., 2007). In exhumation rate of 1 mm a 1, thus lap within their error range, which summary, we suggest that the Seck- indicating increasing exhumation usually implies fast exhumation. This auer Tauern had already cooled rates from the middle to the late Miocene. Regardless of the topo- graphical amplification factor, the (a) AFT (b) AHe corrected exhumation rates of the AFT data are slightly higher than ~0.5 mm a 1 (Fig. 4a) and thus equal to those calculated from the AER. The AHe data from the same profile show a relationship between topographical amplification and exhumation rate (Fig. 4b), indicating that the exhumation rate calculated from the AER is an overestimate. Assuming a topographical steady state (amplification factor = 1), the corrected exhumation is ~0.8 mm a 1 (star in Fig. 4b). We suggest that the lower temperature Modelled AFT age (Ma) Modelled AHe age (Ma) Exhumation rate (km/Ma) Exhumation rate (km/Ma) isotherms were warped sufficiently to affect the AHe ages but not the AFT data, which may be attributed to the topography. Stuwe€ et al. (1994) suggested a significant influ- Topographic amplification factor Topographic amplification factor ence on the 110 °C isotherm of a wavelength of ~20 km, a relief of (c) Seckauer tauern (d) Niedere tauern ~3 km and denudation rates of ~ 1 0 1kmMa . However, these condi- sample Ju9 sample So6 tions are clearly not met in the Nie-
50 HePRZ HePRZ dere Tauern, which are characterised
APAZ APAZ by a topographical wavelength of 100 ~3–6 km and a relief that hardly 150 exceeds 2 km. Independent of the chosen palaeo-topography (ampli- 200 fication factor), corrected exhuma- 50 40 30 20 10 0 25 20 15 10 5 0 tion rates from the AHe system Time (Ma) Time (Ma) exceed those from the AFT system (Fig. 4a,b). This observation implies Fig. 4 Topographical correction of (a) AFT and (b) AHe data from Profile A that there might have been an epi- ° 1 (Figs 2 and 3b) assuming an initial crustal geothermal gradient of 20 Ckm (for sode of faster cooling at the late details see Glotzbach et al., 2015). We have modelled, for a large variety of param- Miocene/Pliocene boundary, which eter combinations (exhumation rates and topographical amplification factors), ther- can also be seen with conventional mochronological ages and resulting apparent exhumation rates from the AER. The modelled and observed AER have been compared to find the best fitting parameter thermal-history modelling (Fig. 4d). combination, which is marked as a thick red line. The thin red lines are the 95% confidence intervals. The thin black lines are modelled AFT and AHe ages assum- Implications for the study area ing continuous exhumation. The small stars are the modelled parameter combina- tions (exhumation rate and topographical amplification factor). (c,d) Thermal- During the early and middle Mio- history inverse models for (c) the Seckauer and (d) the Niedere Tauern using the cene, the Eastern Alps were affected HeFTy programme (Ketcham, 2005). Input parameters: central AFT age with 1r by large-scale extension that resulted error, corrected AHe ages, track length distributions and Dpar values as kinematic in block segmentation and extrusion parameter. Light grey paths: acceptable fit; dark grey paths: good fit; black line: towards the Pannonian basin in best fit. the east (Frisch et al., 2000b). This
© 2016 John Wiley & Sons Ltd 5 Early to Late Miocene exhumation rates in the Eastern Alps • A. W€olfler et al. Terra Nova, Vol 0, No. 0, 1–9 ......
m) process was accompanied by fast : l
D cooling and exhumation of footwall q units exposed in the Tauern Win- dow (Frisch et al., 2000b). However,
(L) Dpar ( our results document that exhuma- N tion can also be observed in dis- tinct parts of the hangingwall. m) l While the Niedere Tauern experi- enced exhumation at a rate of
and longitude coordinates 1 degrees of freedom (where ~0.5 mm a , the Seckauer Tauern n m) SD (
l and other areas such as the Gurktal Alps and the Koralpe (Fig. 2) did ) for 2 not experience substantial exhuma- v tion and surface uplift (Frisch et al.,
(Ma) MTL ( 2001). We attribute this to the r
1 influence of the north-directed com-