Introduction

From the late Carboniferous to the early Permian, the transition from the final stage of the Variscan orogeny to the stable continental regime of Pangea in Central Europe led to a change from compressive regimes to extensional/transtensional ones (Ziegler 1990; van Wees et al. 2000; McCann et al. 2006; Timmerman 2008; Kroner et al. 2016, 2020). Paleogeographically, large marine areas of high sea level turned into shallow marine marginal seas, merging in fringes of extensive swamp forests which fossilized to coal measures of Namurian to Westphalian age. In the latest Carboniferous (Stephanian), the humid environments became drier and fluviolacustrine sedimentation took place (e.g., Roscher and Schneider 2006; Paul 2012). Similar sedimentary sequences developed in elongated troughs inside the orogenic belt as in the Saar–Nahe Basin (e.g., Boy et al. 2012) or in the Thuringian Forest Basin (TFB) (Lützner et al. 2012). During the early Permian, these intramontane basins continued to form and fill up due to exhumation and erosion of the Variscan orogen, increasingly under the control of fault tectonics and magmatism.

The TFB is situated in a central and important geographic position in the Variscan orogenic belt in the Mid-German Crystalline Rise and Saxothuringian Zone between the Saar–Nahe and Kraichgau Basins in the SW and the Saale Basin in the NE (Lützner and Kowalczyk 2012; Andreas 2014). These basins were thoroughly studied during the last > 70 years for their geological structure, volcano-sedimentary successions (e.g., Lützner et al. 2012), magnetostratigraphy (e.g., Menning and Bachtadse 2012), and biostratigraphy (e.g., Martens 2012; Schneider and Werneburg 2012). Therefore, they are considered type sections with which the late Carboniferous to Permian/Triassic time period can be studied in detail. In particular, the TFB includes the Carboniferous–Permian boundary and the transition from the late Permian to the Triassic. Besides numerous small gaps in sedimentation, there was a long-lasting middle Permian hiatus detected in previous studies that is present in diverse European basins, particularly in the central Southern Permian Basin in Northern Germany (e.g., van Wees et al. 2000; Geißler et al. 2008). In parts of mid-European basins, sedimentation restarted close to the major geomagnetic Illawarra Reversal (e.g., Hounslow and Balabanov 2018) in the Late Permian. Therefore, because it covers important intervals of Permian sedimentation and non-sedimentation (“intra-Rotliegend hiatus”), the TFB is exceptionally well suited to further our understanding of sediment fill and magmatism during late Paleozoic geodynamic evolution.

Geochronological studies have already been performed in the TFB (Goll and Lippolt 2001; Lippolt et al. 2009; Lützner et al. 2007; Zeh and Brätz 2000), in the Saar–Nahe Basin (von Seckendorff 2012), and in the Saale Basin (Breitkreuz et al. 2009). However, recent developments in zircon U–Pb CA–ID–TIMS (chemical abrasion–isotope dilution–thermal ionization mass spectrometry) dating allow precise determination of crystallization ages with an accuracy of 0.1% (Condon et al. 2015; Schaltegger et al. 2015). This method is, thus, particularly well suited to establish chronological sequences of sedimentary basins if zircon-bearing volcanic intercalations are available. We present new geochronological data from seven felsic volcanic rock samples of the TFB with zircons dated by the U–Pb CA–ID–TIMS method to:

  • 1. Define the age range of volcanic activity in the TFB and compare it to the existing biostratigraphic data,

  • 2. Constrain the gap between the end of major volcanic activity with associated intramontane sediments and the voluminous post-Illawarra paleomagnetic Reversal sediment accumulation of the late Permian stage, and,

  • 3. To compare the timing of the main volcanic activity across Permian basins in Central Europe.

Geological setting

Carboniferous–Permian volcanism in Central Europe

At the end of the Variscan orogeny, an extensive phase of Carboniferous–Permian magmatism led to the extrusion of thick volcanic successions across Europe and NW Africa (Fig. 1a, e.g., Oslo Rift, Scotland, England, Southern Permian Basin/N Germany, W and central Germany, Saxony, N Bohemia, NW part of the Polish Basin, SW Poland, NE Czech Republic, N Spain, Morocco, France, N Italy), with a maximum volume in the late Carboniferous (Stephanian) to early Permian. The eruption of volcanic rocks was accompanied by massive synextensional granitic intrusions and hydrothermal mineralisation (Lorenz and Nicholls 1976, 1984; Hoth et al. 1993; Marx 1995; Marx et al. 1995; Benek et al. 1996; Doblas et al. 1998; Breitkreuz and Kennedy 1999; Lorenz and Haneke 2004; Timmerman 2004, 2008; McCann et al. 2006; Breitkreuz et al. 2007, 2009; Geißler et al. 2008; von Seckendorff 2012; Hoffmann et al. 2013; Willcock et al. 2013; Schmiedel et al. 2015).

Fig. 1
figure 1

Simplified geological maps. a Geological map of the main tectonic elements of Central Europe in relation to the distribution of Late Carboniferous and Early Permian sedimentary basins and volcanic rocks (based on McCann et al. 2006; Mazur et al. 2010). Major tectonic elements are ACF Alpine–Carpathian front, VF Variscan deformation front, RH Rheno-Hercynian Zone, ST Saxothuringian Zone, MZ Moldanubian Zone, TB Tepla–Barrandian Zone, EEC East European Platform, FSS Fennoscandian Shield, TTZ Tornquist–Teisseyre Zone, OR Oslo Rift. Two sub-basins of the Southern Permian Basin are labeled as NGB North German Basin, PB Polish Basin, and the main sedimentary basins of the Saxothuringian Zone are SNB Saar–Nahe Basin, TFB Thuringian Forest Basin, SB Saale Basin. Thick black lines are tectonic elements; thin black lines represent the coastlines of Central Europe. b Geological Map of the Thuringian Forest Basin. Sample locations for new zircon U–Pb CA–ID–TIMS dating are shown by black circles with sample numbers in yellow boxes. (Lützner et al. 2012, amended). RCC Ruhla Crystalline Complex, THG Thüringer Hauptgranit

In Europe, Stephanian to Permian wrenching resulted in thermal destabilization of the lithosphere, intense fracturing of the crust, disruption and erosion of its sedimentary cover, and regional uplift in the area of the future Southern Permian Basin (van Wees et al. 2000; Kroner et al. 2016). Transtension led to the formation of more than 70 rift basins across Europe. In the North German Basin (as part of the Southern Permian Basin), volcanism between ca. 300 and 290 Ma is closely related in time and space to an episode of major tectonic activity (right-lateral WNW- to W-trending transtensional strike-slip movements; Plein 1993; Gast and Gundlach 2006; McCann et al. 2006). The extensional tectonics were important for supplying pathways for intrusive and extrusive magmas.

The origin of the extensive and voluminous Carboniferous–Permian magmatism in Central Europe has been hotly debated for more than 30 years with diverse theories arising (see, e.g., Doblas et al. 1998; Timmerman 2004; McCann et al. 2006; Geißler et al. 2008). This is partly due to the large geochemical and mineralogical variety of the plutonic and volcanic rocks, but likely also due to the structural and lithospheric complexity of this part of Variscan Europe (Lorenz and Nicholls 1976, 1984; Breitkreuz and Kennedy 1999; Breitkreuz and Mock 2004; von Seckendorff 2004, 2012; McCann et al. 2006; Breitkreuz et al. 2007, 2009). A superplume has been invoked by Doblas et al. (1998) as a background mechanism of this magmatic episode. Rifting, thinning of highly asymmetric lithosphere, local decompression, and lithospheric mantle as the major source of the middle European Carboniferous–Permian volcanic rocks were defined by Neumann et al. (2004). Mantle sources metasomatized by subduction-related material were inferred by Timmerman (2004) to explain the geochemical character of the volcanic rocks in NE Germany that does not agree with an intra-plate setting.

The widespread Stephanian–early Permian (305–285 Ma, Kasimovian–late Artinskian) alkaline to calc-alkaline intrusive and extrusive magmatism of the central European Variscan area and its northern foreland (e.g., the North German Volcanic Complex) is, based on trace element ratios, mantle-derived and shows evidence of moderate crustal contamination (e.g., Marx et al. 1995; Benek et al. 1996; Breitkreuz and Kennedy 1999; McCann et al. 2006). In the internal Variscides (South of the Variscan front, VF in Fig. 1a), mafic rocks are rare with felsic-intermediate compositions being dominant; whereas in the external areas, calc-alkaline SiO2-rich volcanic rocks are more common (e.g., Timmerman 2008).

The late Stephanian to earliest Permian sequence (volcano-sedimentary and volcanic rocks as well as clastic fill) reached its maximum thickness in the TFB (ca. 2000 m, a value comparable to maxima in N Germany.) The igneous rocks are lavas and pyroclastic rocks of mainly trachyandesitic to rhyolitic composition (Andreas et al. 1974; Obst et al. 1999; Obst and Katzung 2000; Lützner et al. 2007, 2012).

Thuringian forest basin (TFB)

Geological overview

The TFB is of exceptional value for the age determination of late Variscan Carboniferous–Permian basins because it is well investigated concerning lithostratigraphy (Lützner et al. 2012) and biostratigraphy (e.g., Schneider and Werneburg 2012; Schneider et al. 2020). The TFB also shows recurrent strong intra-basinal volcanism that produced piles of lava flows and tuffs hundreds of meters thick. Precise and accurate age dating of the TFB formations is, therefore, critical for the correct correlation and stratigraphic categorisation of a variety of other European late Carboniferous/early Permian basins.

Tectonically, the Thuringian Forest represents an uplifted block bordered by mainly NW–SE-striking faults and flexures that predominantly coincide with the morphological slopes of the range. The maximum altitudes are reached in the middle part (Großer Beerberg, 982 m) and in the northwest (Großer Inselsberg, 916 m), with decreasing altitudes (ca. 700 m) between both peaks and further to the southeast. The forelands are made up of Triassic sedimentary formations (Buntsandstein, Muschelkalk, Keuper), in part with a narrow strip of late Permian Zechstein deposits along marginal flexures; whereas, within the elevated range, Stephanian and early Permian successions and their Variscan basement are revealed (Fig. 1b).

The Stephanian–Rotliegend (late Carboniferous–early Permian) filling of the TFB disconformably overlies the Variscan basement. On the southeastern margin, the basement consists of Proterozoic and Paleozoic up to lower Carboniferous sedimentary units of the Saxothuringian Zone which were folded and underwent a low-grade Variscan metamorphism (Bankwitz and Bankwitz 2003; Meinel 2003).

To the northwest, the basin filling covers the “Thüringer Hauptgranit”, an early post-kinematic granite pluton of early Carboniferous age (330.3 ± 3.4 Ma to 337 ± 4 Ma from zircon Pb–Pb evaporation and dark mica Ar–Ar, and K–Ar dating; Zeh and Brätz 2000; Zeh et al. 2000; Goll and Lippold 2001) with a contact belt against the early Paleozoic in the Southeast (Fig. 1b; Bankwitz and Kaemmel 1957). The composition of the “Thüringer Hauptgranit” is heterogeneous and varies between quartzdioritic, granodioritic, and leucogranitic (Meinel 2003). These plutonic rocks had been eroded and deeply weathered before the development of the TFB (Lützner et al. 2012). The extension of the TFB spreads beyond the recent tectonic borders of the Thuringian Forest at least 12 km to the Northeast (borehole Gotha 1), and at least 20 km to the Southwest (boreholes Metzels, Thüringen Süd, Rannungen, not shown in Fig. 1). The extension to the NW is bordered by the Ruhla Crystalline Complex (RCC, Fig. 1b). Paleocurrent patterns indicate that the RCC was one of the source areas for the TFB sedimentary formations.

Lithostratigraphy

Figure 2 illustrates the lithostratigraphic subdivisions of the TFB. The lowermost formations (Georgenthal Formation in the NW and Möhrenbach Formation in the SE, thought to be coeval) reflect a first culmination of magmatic activity and are both dominated by intermediate volcanic rocks with additional rhyolites present in the Möhrenbach Formation. Thin sedimentary intercalations were used for internal correlations (Lützner et al. 2012: p. 430; Andreas 2014: p. 96). The overlying Ilmenau Formation displays local sedimentary beds alternating with volcanic and pyroclastic rocks of bimodal basaltic and rhyolitic composition (Lützner et al. 2012).

Fig. 2
figure 2

Lithostratigraphic subdivision of the Carboniferous–Permian sediments of the Thuringian Forest Basin, presented as a schematic NW–SE cross–section. Circled letters label biostratigraphically important fossiliferous horizons: A—Öhrenkammer horizon (Georgenthal Fm.), B—Ilmtal horizon (Möhrenbach Fm.), C—Lindenberg horizon (base of Ilmenau Fm.), D—Sembachtal horizon (top of Ilmenau Fm.), E—Manebach fossiliferous horizons (Manebach Fm.), F—Sperbersbach horizon (lower Goldlauter Fm.), G—Gottlob horizon (upper Goldlauter Fm.), H—Lochbrunnen horizon (lower Oberhof Fm.), I—Wintersbrunn horizon (upper Oberhof Fm.), J—Bromacker horizon (Tambach Fm.). Numbers in yellow boxes represent new sample locations for zircon U–Pb CA–ID–TIMS dating

Extrusive volcanic rocks are extremely scarce in the two subsequent formations, i.e., in the coal-bearing Manebach Formation and in the fluvial–lacustrine Goldlauter Formation. Volcanic events are restricted to several tuff layers, which are used as tephrostratigraphic marker horizons (Andreas and Haubold 1975).

In the TFB, a second culmination of Rotliegend volcanism occurred during the following Oberhof Formation where volcanic rocks alternate with volcaniclastics and subordinately with alluvial red beds and lacustrine gray beds. The basal group of volcanic rocks within the Oberhof Formation consists of an association of basic, intermediate, and acidic volcanic rocks that is intercalated within the lower sedimentary beds of fluvial and lacustrine origin. The middle and upper parts of the Oberhof Formation are dominated by rhyolites and associated volcaniclastics. So-called Older Rhyolites prevail in the middle part. They are characterized by medium to large feldspar phenocrysts. Finally, the so-called Younger Rhyolites constitute the upper part of the Oberhof Formation. These rocks have typically small qtz–fsp phenocrysts of about 1 mm size. However, a few rhyolite domes intruding the Younger Rhyolite level also have large qtz–fsp phenocrysts.

The Rotterode Formation consists of red sandstones with conglomeratic intercalations that surround several rhyolitic bodies and tuff layers. This formation was deposited on the Oberhof Formation with a low-angle unconformity separating the two formations. The angular unconformity was likely related to the tectonically induced uplift of the Ruhla crystalline complex, which caused rejuvenation of the relief and resulted in deposition of granite-derived sands and gravels (Thomson and Zeh 2000; Andreas and Lützner, 2009).

The Tambach Formation starts with coarse conglomerates that cover an eroded surface. The conglomerates grade upwards into red sandstones and mudstones that contain the youngest fossil beds of stratigraphic significance in the TFB. The fossil-bearing red beds at the Bromacker locality are famous for tetrapod footprints and skeletal remains of terrestrial amphibians and reptiles. Paleontological projects were performed by T. Martens in cooperation with experts from the USA (Martens 2001; Eberth et al. 2000 and others). The state of biostratigraphic correlation was summarized by Schneider et al. (2020). They consider for the Bromacker Horizon an Artinskian age, which is not proven by radioisotopic studies. Volcanic rocks are missing in the Tambach Formation; volcaniclastic rocks are also presently unknown.

The Eisenach Formation only occurs on the western flank of the Ruhla Crystalline Complex and the adjacent Werra basin. Along the basin fringes, marginal alluvial fans are interfingered with sandy mudstones.

The Neuenhof Formation is the uppermost Permian formation of continental red beds. It consists of a 5–10 m thin-bedded alternation with fine gravels, sandstones, and siltstones of a mainly fluvial environment. Locally intercalated cross-bedded sandstones are probably of aeolian origin. At the top, the red-colored sandstones gradually pass over into gray sediments that were reworked by the invading Zechstein Sea. A scarce settling with brachiopods has been proven in places. The Kupferschiefer black shale follows with a sharp boundary.

Magnetostratigraphically, all formations of the TFB, including the Eisenach Formation, belong to the Kiaman Reversed Superchron. Only the red beds of the Neuenhof Formation proved to be part of the Permian–Triassic Mixed Superchron (Isozaki 2009; Menning et al. 1988). Thus, the Illawarra Reversal (266.66 ± 0.76 Ma, Hounslow and Balabanov 2018) took place during a long hiatus (break of sedimentation, gap due to onlapping cover on morphology, and/or erosion). Nevertheless, the Neuenhof Formation has the appearance of the Rotliegend facies (Lützner et al. 2012), which is directly covered by the basal conglomerate and black shale (Kupferschiefer) of the Zechstein sea.

Samples and methods

Samples

We aimed to sample the whole Stephanian and early Permian sequence of the TFB by sampling various rhyolitic and rhyodacitic volcanic rocks (Figs. 1b, 2). From each formation, one or two samples were collected (Table 1). We preferred to sample acidic rocks with medium- to large-sized phenocrysts of quartz and feldspar because these rocks generally yielded sufficient zircon crystals. We briefly comment on the sample sites in the lithostratigraphic order (Fig. 2).

Table 1 List of samples

In the lowest section, rhyolites from the Möhrenbach and the Ilmenau Formations of the Southeastern TFB region were sampled. Samples were taken from two rhyolite outcrops within the Möhrenbach Formation: (i) from road cuts near Stützerbach (sample 7) and (ii) from the quarry at the Hundskopf hillside (sample 6). According to several authors (Andreas et al. 1974; Michael 1972; Lützner et al. 2012), the Stützerbach Rhyolite is situated in the middle section of the Möhrenbach Formation and thus has a somewhat lower position than the Hundskopf Rhyolite. The latter is considered to be a part of the Kienberg Rhyolite that occurs as a dissipated system of dykes, sill-like bodies, vent fillings and lava sheets. Both samples are good representatives of the Möhrenbach Formation.

The Ilmenau Formation (occurring above the Möhrenbach Formation) contains the Kickelhahn Rhyolite that typically has plenty of medium-sized quartz and feldspar phenocrysts and is well exposed at the type locality “High Hill Kickelhahn” near Ilmenau where sample 5 was taken. This formation is associated with thick pyroclastics.

The predominantly sedimentary Manebach and Goldlauter Formations were not sampled because these units lack extrusive volcanic rocks.

Extensive magmatic activity (rhyolites, rhyodacites, felsic intrusions, and various pyroclastic deposits) occurred in the Oberhof Formation. Two localities in the uppermost lava flow of the “Older Rhyolite Series” were selected near Luisenthal to sample a representative rhyolite in the middle section of the Oberhof Formation (based on the revised Geological Map; Enderlein and Ziegenhardt 1975; scale of 1: 25,000). Sample 2 was taken from the Klotzegrund valley and sample 3 from outcrops at the mouth of the Lehmgrund valley. Both samples represent different parts of a uniform rhyolite body.

The Rotterode Formation comprises two fault-bound rhyolite bodies that did not contain zircon. In addition, several rhyolite dykes penetrate the red beds of the Rotterode Formation. They have been mapped in detail by Patzelt (1966) and were adopted in a revised map of the Rotterode Formation by Andreas and Lützner (2009). We sampled one of the rhyolite dykes in an outcrop in the village Struth-Helmershof (sample 1). It should indicate the median or minimum age of the Rotterode Formation.

The red beds overlying the Rotterode Formation, i.e., the Tambach and Eisenach Formations, do not contain volcanic rocks. Only the Elgersburg Formation, which is considered as an equivalent of the Tambach Formation, contains two rhyolite bodies and a basaltic layer in the lower part (Schwalbenstein Conglomerate Member). The Elgersburg Rhyolite represents the uppermost volcanic rock in the TFB. Unfortunately, separated zircons from this rhyolite had xenomorphic forms; idiomorphic grains were lacking. Therefore, we excluded this sample from dating and instead used the Wolfstein Rhyolite, which appears near the base of Elgersburg conglomerates (No. 4). It was first described as “Preußenshöhe Porphyr” by Zimmermann (in Loretz et al. 1908, p. 122) and Deubel (1960, p. 429). The Wolfstein Rhyolite is characterized by medium to large phenocrysts of quartz and feldspar and contains idiomorphic zircon crystals suitable for high-precision dating.

Methods

Recent development in ID-TIMS dating substantially reduced the analytical error of zircon dating, especially through the application of chemical abrasion (CA) on single zircon grains prior to their analysis. It was demonstrated in several studies that this pre-treatment considerably reduces Pb loss that—if not detected—would result in exceedingly young ages (Mattinson 2005; Schaltegger et al. 2015; Tichomirowa et al. 2019). In addition, due to efforts of the EARTHTIME initiative (http://www.earthtimetestsite.com/), the inter-laboratory bias of zircon ages has been reduced to 0.1% through the use of precisely and accurately calibrated tracers (Condon et al. 2015) and the standardization of raw data correction by the development and supply of open source programs (McLean et al. 2015). As a result, the accuracy and external reproducibility of zircon ages from participating laboratories can be ≤ 0.1% and can be proved by their reported international zircon standard ages. Consequently, the CA–ID–TIMS method has the ability to produce the most precise and accurate ages (Schaltegger et al. 2015).

U–Pb CA–ID–TIMS dating was performed in the TIMS lab at the TU Bergakademie Freiberg (Germany). Selected zircon grains with the most idiomorphic forms (ca. 15–30 per sample) were annealed for 72–96 h at 850 °C and subsequently chemically abraded for 12–15 h at 180 °C with concentrated HF and HNO3. This procedure dissolves crystal domains with the most intense radiation damage, which are suspected to have experienced post-crystallization lead loss (Mattinson 2005). After this chemical leaching procedure, single idiomorphic grains are transformed into fragments that represent those zircon parts with the lowest uranium content (Mattinson 2005). Afterwards, the acid, together with dissolved zircon material, was completely pipetted out. 3.5 N HNO3 was added to the remaining zircon fragments, and the mixture was left for 30 min at 50 °C to remove surface lead. Several cleaning cycles with water, combined with repeated ultrasonic treatment were conducted before single zircon fragments were selected for further processing.

Single zircon fragments were washed with 3 N HNO3 and transferred into cleaned microcapsules with a small drop of this fluid and 3 drops of concentrated HF. Samples were spiked with 6–10 mg of a 205Pb–233U–235U tracer solution (ET535 at TU Bergakademie Freiberg; Condon et al. 2015). For final dissolution, the microcapsules were placed in Parr vessels and heated to 200° C for 48 h followed by drying at 130 °C and then re-dissolution in 6 N HCl for 24 h at 200 °C to transfer them into chlorides. After repeated drying, the samples were dissolved in 10 drops (ca. 50–60 mg) of 3.1 N HNO3 and transferred into 50 µl micro-columns for column chemistry. U and Pb were separated from the rest of the sample by anion exchange chromatography (AG 1-X8, 200–400 mesh) using HCl and H2O. The U- and Pb-containing fraction was loaded on degassed rhenium filaments with a drop of silica gel and measured with a Finnigan Mass Spectrometer MAT 262 using a secondary electron multiplier (SEM). The common Pb was corrected using the model of Stacey and Kramers (1975) instead of applying a specific lab blank value because this correction resulted in the most concordant values for reference materials (two international standards—see below).

Ages of two international zircon standards (Temora 2, 91500) were acquired contemporaneously with those of the Thuringian Forest volcanic rocks to ensure their accuracy. Similar to the samples, both standards were pre-treated with chemical abrasion (CA) to get comparable results. Data are given as 206Pb/238U age-ranked plots in Fig. 3. The Temora 2 standard has an accepted age of 416.8 ± 0.3 Ma (Black et al. 2003). The majority of our single grain ages have a spread of <  ± 0.3% with a few slightly older ages (Table 2, Fig. 3). We obtained a weighted mean age of 417.34 ± 0.28/0.33/0.56 Ma for this standard. The standard 91500 was dated at 1062.4 ± 0.4 Ma (Wiedenbeck et al. 1995) and 1063.6 ± 0.3 Ma (Schoene et al. 2006). In our data, a few ages are slightly (< 1%) younger than the main age cluster, indicating the presence of small lead loss despite CA pre-treatment (Table 2). The remaining data used for the weighted mean age calculation have a spread of < 0.3% and yielded a weighted mean age of 1064.6 ± 0.55/0.60/1.30 Ma (Fig. 3). Altogether, our results are within 0.1% of the accepted values, assuring the accuracy of our CA–ID–TIMS ages at the 0.1% level.

Fig. 3
figure 3

206Pb/238U ages of two international standards obtained by CA–ID–TIMS. Single ages are shown by bars representing 2σ errors (including analytical, tracer, and decay constant uncertainties); unfilled bars are considered outliers. All uncertainties associated with weighted mean 206Pb/238U ages are reported at the 95% confidence level and given as ± x/y/z, with x as analytical uncertainty, y including tracer calibration uncertainty, and z including 238U decay constant uncertainty. The accepted mean ages, together with their 2σ uncertainties for these standards (see text for values) are given in green for comparison

Table 2 U–Pb CA–ID–TIMS analyses of standard and sample zircons; data used for 206Pb/238U-weighted mean age calculations are given in bold

Results

The U–Pb isotopic results of seven samples are given in Table 2 and presented in Fig. 4 as 206Pb/238U age-ranked plots for each individual sample. Eight to fifteen grains were analyzed for each sample. Mean sample ages were calculated from age clusters that are interpreted as crystallization ages. The error of weighted mean ages is given as ± x/y/z, where x is the internal 2σ measurement error, y is the internal 2σ error plus tracer calibration uncertainty, and z additionally includes the decay constant uncertainty (Schoene et al. 2006).

Fig. 4
figure 4

New single grain zircon analyses and 206Pb/238U-weighted mean ages for the Thuringian Forest Section. Each single zircon grain/fragment analysis including its 2σ uncertainty is represented by a horizontal bar. Grey bars are not included in the calculation of weighted mean ages. The weighted mean ages are shown by vertical lines. Uncertainty in the weighted mean ages are reported as x/y/z, where x is the internal 2σ measurement error, y is the internal 2σ error plus tracer calibration uncertainty, and z additionally includes the uncertainty of the decay constant (Schoene et al. 2006)

Nine zircons were dated from sample 1. The youngest age was discarded and interpreted as slight lead loss because of its offset from the remaining eight results. The other eight analyses yielded ages from 295.7 ± 0.2 Ma to 297.6 ± 0.9 Ma. The youngest three of these eight dates define a cluster of identical ages with a weighted mean age of 295.8 ± 0.1/0.2/0.4 Ma [mean square of weighted deviates (MSWD) = 0.1], which we interpret as the crystallization age. The five remaining dates show a slight increase from 296.3 ± 1.5 to 297.6 ± 0.9 Ma. These deviating ages are interpreted as zircons incorporated during intrusion from slightly older magma batches (antecrystic grains, Miller et al. 2007).

Analysis of six individual zircons from sample 2 yielded a statistically significant cluster with a weighted mean age of 296.8 ± 0.2/0.2/0.4 Ma (MSWD = 0.8). Two slightly older ages and one slightly younger age deviate from this cluster and have been discarded from the mean age calculation and are interpreted as antecrystic grains and grains affected by lead loss, respectively.

The weighted mean age of sample 3 was calculated from the youngest age cluster of three zircons yielding 296.9 ± 0.1/0.2/0.4 Ma (MSWD = 0.9, n = 3). Four grains are distinctly older, varying from 297.1 ± 0.3 Ma to 301.2 ± 1.2 Ma. These ages do not overlap with the calculated mean age and are interpreted as antecrystic (Miller et al. 2007) and/or inherited ages.

Ten single zircon grains were analyzed from sample 4. Three of them are more than 5 Ma older than the rest and are interpreted as inherited grains. Four of the remaining dates show a tendency towards younger ages, while the other three grains define an age cluster at 299.4 ± 0.3/0.3/0.5 Ma (MSWD = 1.3), interpreted as the crystallization age of this sample. If this interpretation is valid, severe lead loss has to be assumed for the remaining four grains (the youngest date results in 296.5 ± 0.4 Ma, which is 1% below the defined mean age). Therefore, lead loss is most pronounced in zircons from sample 4 but also present in zircons from samples 1, 2, 5, 6, and 7.

The analysis of 15 zircon grains from sample 5 yielded ages that range from 298.4 ± 0.2 Ma to 302.8 ± 0.1 Ma, with one inherited age at 319.3 ± 0.2 Ma. Three other grains are significantly older than the majority of ages and are interpreted as antecrystic zircons. The remaining eleven zircons reveal a much smaller scatter from 298.4 ± 0.2 Ma to 299.3 ± 0.1 Ma. We interpret the younger dates as affected by lead loss and calculated a weighted mean age of 299.3 ± 0.1/0.1/0.3 Ma (MSWD = 1.2) using the five oldest grains from this population.

Nine zircon grains were dated from sample 6. The statistically significant age cluster includes seven single grain ages and yields a weighted mean age of 300.0 ± 0.3/0.3/0.5 Ma (MSWD = 0.9). One age is significantly younger than this cluster and is interpreted as having been affected by lead loss; another age is significantly older and is interpreted as an antecrystic grain.

Ten zircon grains were analyzed from sample 7 and show ages from 298.0 ± 0.6 Ma to 299.9 ± 0.7 Ma. The youngest ages are not identical but show a tendency toward lower ages that is interpreted as having been affected by lead loss. A well-defined age cluster is formed by the four oldest grains of this sample at 299.7 ± 0.3/0.3/0.4 Ma (MSWD = 0.2).

Discussion

New radiometric zircon ages

We discuss the results in stratigraphic order, first with regard to the lithostratigraphic subdivision of the TFB, then to the biostratigraphic aspects. In the following, we use the error expressed as z that also includes the error in tracer calibration and that of the decay constant (see above). This allows us to compare our ages to those obtained with other tracers and other dating systems.

Samples 7 and 6 yielded identical (within errors) 206Pb/238U mean ages (299.7 ± 0.4 Ma vs. 300.0 ± 0.5 Ma) representing the age of the Möhrenbach Formation. Considering the uncertain lithostratigraphic relation between the two sample sites (compare with Section “Geological setting”), the results of samples 7 and 6 clearly indicate an age in the range between 300.5 Ma and 299.3 Ma for the middle to upper part of the Möhrenbach Formation. Consequently, the Möhrenbach Formation was deposited during the latest Carboniferous, i.e., the Gzhelian (Standard Global Chronostratigraphic Scale) or Stephanian (West European regional scale).

Our new age of the Ilmenau Formation (sample 3) is the first precise depositional age for this unit. The sample dates the Kickelhahn Rhyolite near Manebach to 299.3 ± 0.3 Ma. This age overlaps with the recorded age of the Carboniferous–Permian boundary (298.90 ± 0.15 Ma; International Chronostratigraphic Chart 2019/05; the age was obtained by the same U–Pb CA–ID–TIMS method using the 535 spike; Ramezani et al. 2007). Thus, the age of the Kickelhahn Rhyolite provides important information for the stratigraphic position of the C–P boundary in the Thuringian Forest Basin, which has long been discussed as not exactly known.

After mapping the TFB on a scale of 1: 25,000, Beyschlag (1895: p. 599) concluded that “late” Carboniferous deposits are missing in the Thuringian Forest and that the entire sequence belongs to the Rotliegend (synonym for “Lower” Permian). Later, Gothan (1928) described from the basal, coal-seam-bearing sediments at the Öhrenkammer near Ruhla (letter A in Fig. 2) a peat-forming flora with Stephanian characteristic forms, which (supposedly) no longer occur in the Rotliegend. The Stephanian age of the Öhrenkammer and Ilmtal horizons (letters A and B, Fig. 2) was also confirmed by the fauna of insects and amphibians (Schneider 1996; Schneider and Werneburg 2006, 2012). From then on, the C–P boundary was assumed to be undefined in the overlying volcanic sequence. Barthel and Rößler (2010) described plant remains from a silicified tuff horizon at the base of the Ilmenau Formation on the Lindenberg hill near Ilmenau (letter C in Fig. 2) as a typical Rotliegend flora. According to this assessment, the C–P boundary was "tentatively" assumed at the base of the Ilmenau Formation (Lützner et al. 2012; Gebhardt et al. 2018; Schneider et al. 2020).

The Kickelhahn Rhyolite (near Manebach village) is only covered by 10–15 m of tuffaceous conglomerates, above which the Manebach Formation begins with grey fluviatile sandstones. In the northwestern part of the Thuringian Forest, above an approx. 100 m thick rhyolite, 30 m of grey to black, fossil-rich lake sediments with tuff layers and subeffusive mafic rocks follow (formerly called “Melaphyr”). These rocks are covered by gray sandstones of the Manebach Formation. Considering the low thickness of post-rhyolitic volcanoclastic sedimentation, we estimate that the C–P boundary should be positioned approximately at the boundary between the Ilmenau and Manebach Formations (cp. Fig. 5).

Fig. 5
figure 5

Compilation of age data from the Rotliegend section of the Thuringian Forest. Radiometric age data are from a Goll and Lippolt (2001), b Zeh and Brätz (2000), c Lützner et al. (2007). The blue bars show the tentative assessment of the Thuringian Forest formations to the International Stratigraphic Chart 2019 after d Schneider and Werneburg (2012), e Menning et al. (2016; explanatory notes in Gebhardt et al. 2018) and f Schneider et al. (2020). Our new CA–ID–TIMS ages are shown in red. Numbers in yellow boxes are sample numbers for new zircon U–Pb CA–ID–TIMS dating. Obviously, most of the previous age data and assessments are remarkably younger than the new zircon ages

The U–Pb zircon age of a silicified tuff layer from the middle part of the Manebach Formation was determined to be 297.8 ± 2.0 Ma by laser ICP–MS (Trümper et al. 2020). This finding fits well with the result of sample 5 (Fig. 6).

Fig. 6
figure 6

CA–ID–TIMS U–Pb zircon age data of the TFB and surrounding basins in relation to biostratigraphic zoning based on amphibians and insects

Our new data suggest an age difference of about 2.5 Myr between sample 5 und sample 3. During this time, the Manebach and Goldlauter Formations filled a pull-apart basin with paludal, lacustrine, fluvial, and alluvial sediments. The thickness of both formations is nearly 1000 m, i.e., the rate of subsidence (uncorrected) corresponds to about 400 m/Myr. Such subsidence rates are common in grabens/rift settings (e.g., Einsele 1992; Doglioni et al. 1998) and can reach extreme values in pull-apart basins like in the Ridge Basin in California (e.g., Yeats 1978: 2000 m/Myr).

Samples 3 and 2 reflect the culmination of the second volcanic phase within the TFB, which occurred in the range between 296 and 297 Ma. In future studies, it will be possible to better constrain the age of the Oberhof Formation by dating the rhyolites at the base (Heuberg- and Inselsberg Rhyolite) and the coarse crystalline Greifenberg rhyolite in dome structures near the upper boundary at the site of Großer Buchenberg near Tambach-Dietharz.

The age of the Komberg Rhyolite intrusion is provided by sample 1 and dates the red beds of the Rotterode Formation. As this formation clearly overlies the Oberhof Formation, the youngest zircon age of this sample (295.8 ± 0.4 Ma) is in agreement with the stratigraphy.

Sample 4 was dated with the idea of determining a representative age for the Elgersburg Formation. However, the obtained age is much older than expected; similar to that of the Möhrenbach Formation. The Wolfstein Rhyolite is directly bordered by the Schöffenhaus fault that dislocates the Elgersburg Formation against members of the Ilmenau Formation and the Oberhof or (in question) Rotterode Formations. Thus, we suggest that the Wolfstein Rhyolite could be interpreted as a small, tectonically uplifted block. Similar small uplifted blocks are well known from faults along the southwestern margin of the Thuringian Forest. The rock type of the Wolfstein Rhyolite with especially large phenocrysts of quartz and feldspar occurs within the Möhrenbach Formation as a widely ramified rhyolite dyke system. It has been described as Meyersgrund Rhyolite by Zimmermann (in Loretz et al. 1908: p. 100, 126), named after a tributary valley of the Ilm river.

In most occurrences of this Meyersgrund Rhyolite, the large quartz and feldspar phenocrysts are very fresh. In contrast, the Wolfstein Rhyolite is found in a highly altered state. In the small outcrop area, the rock is often so strongly weathered that the phenocrysts can only be recognized as hollow forms. It is possible that the fluid movements on the Schöffenhaus fault zone led to a particularly strong alteration. The result of sample 4 is also difficult to interpret due to the wide dispersion of the single grain ages and the low cluster density. For example, taking only the two youngest dates as clusters, the resulting age would be almost identical with sample 1, close to the expected age of the Elgersburg Formation.

The new zircon data for the profile of Carboniferous–Permian volcanic rocks in the TFB show older ages compared to the previously published age data (Fig. 5; data in summarized form in Lützner et al. 2012). This reflects the fact that the pre-treatment of the dated zircons by CA removed those components that suffered severe Pb loss. For most dated samples, the still-existing slight Pb loss resulted in few single ages that are 1–2 Ma younger than the calculated mean age (Fig. 4). Exceptions are sample 4 (which still has pronounced Pb loss; many single grain ages are up to 4 Ma younger) and one measurement of sample 6 (that is 5 Ma younger than the mean age). In general, the new zircon ages (single zircon ages as well as mean sample ages) show less age scatter than those of earlier publications. Consequently, the total duration of the volcanic activity in the TFB was considerably shorter. Whereas previously an age interval between 295 and 275 Ma was expected for the interval between the Möhrenbach Formation and the Elgersburg Formation (Lützner et al. 2007), our new data suggest a range between ca. 300 Ma for the Möhrenbach Formation and ca. 296 Ma for the Rotterode Formation. Even if the age of the youngest rhyolite (Elgersburg) has not yet been determined, the total duration of volcanic activity has to be reduced from ca. 20 Myr to ca. 4 Myr. The sedimentary gap up to the Illawarra Reversal (ca. 265 Ma after Isozaki 2009; 266.66 ± 0.76 Ma after Hounslow and Balabanov 2018; base of Upper Rotliegend II) therefore has to be extended to > 25 Myr. This result requires a new view on the early evolution of the Thuringian Forest Basin. The sedimentary hiatus, which can extend locally to 30 Myr in Northern Germany (e.g., between Permian volcanics and the Wustrow Member as the oldest Permian sediment), is at least partially caused by volcano-topographic features (morphology) besides intra-basinal tectonics according to Geißler et al. (2008). Further work is required to evaluate this hiatus in sedimentation for the TFB region as well.

Comparison to biostratigraphic data

The modern state of biostratigraphy in the continental basins of Pangea has been summarized by Schneider et al. (2020). The authors describe “biostratigraphic timelines” which are calibrated for Europe and North America by marine zone fossils and are dated for certain levels by radioisotopic zircon ages (including U–Pb CA–ID–TIMS zircon ages) from the Bashkirian (early Pennsylvanian: ca. 319 Ma) to the late Asselian (early Cisuralian: ca. 295 Ma). On this basis, stratigraphic correlations can be done at a local and regional scale. However, there are only a few continental basins at the Pangea margins that are well linked to the global marine standard biostratigraphy. Accordingly, more chronostratigraphic data from Pangea continental basins are urgently needed.

In the TFB, biostratigraphy is well investigated in the sedimentary record. Collecting biostratigraphic data (occurrences) and evaluating biostratigraphic events (first and last appearances) have been used to establish biostratigraphic zoning for insects, especially cockroach wings, for amphibians, vertebrate foot prints and conchostraca as well. However, the correlation with biostratigraphic data in continental Carboniferous and Permian basins is partially problematic (e.g., Schneider et al. 2020). Comparison of Fig. 2 of Schneider et al. (2020) with our new age data for the formations near the Carboniferous–Permian boundary (Möhrenbach Fm. plus Ilmenau Fm. = Gehren subgroup) indicates that both datasets are compatible. However, contradictory interpretations are obtained for the Rotterode Formation: According to our new ages, the Rotterode Formation corresponds to the Asselian Stage while Schneider et al. (2020) place the Rotterode Formation in the Artinskian Stage.

Schneider and Werneburg (2006, 2012) biostratigraphically correlated the Ilmtal fossiliferous horizon in the Möhrenbach Formation (letter B in Fig. 2 and Fig. 6) with the Breitenbach Formation (Saar–Nahe Basin) and the Wettin Subformation (Saale Basin) on the basis of index fossils of the Branchierpeton saalensis–Apateon intermedius assemblage zone (amphibians) and of the Sysciophlebia euglyptica–Syscioblatta dohrni assemblage zone (insect wings). Schneider et al. (2020) have inferred the late Stephanian (late Gzhelian) time from a zircon SHRIMP age of an intrusive rhyolite near Halle (Breitkreuz et al. 2009). The Stephanian age is also confirmed here by the age for the Kickelhahn Rhyolite (Ilmenau Formation).

The next amphibian biozone following upward, i.e., Apateon dracyensis–Melanerpeton sembachense assemblage zone comes from the Sembach fossiliferous horizon near the top of the Ilmenau Formation (letter D in Figs. 2 and 6) and is also reported from the Muse Formation near Dracy St. Loup (Autun basin, France). The upper part of the Muse oil shale is dated between 298.6 ± 0.4 Ma and 298.0 ± 0.4 Ma by the U–Pb CA–ID–TIMS zircon age method (Pellenard et al. 2017; Fig. 6). This fits well to the 299.3 ± 0.3 Ma age of the Kickelhahn Rhyolite when taking into account that the Sembach Horizon is situated at the top of the Ilmenau Formation and that the Kickelhahn Rhyolite has a somewhat lower position.

Furthermore, the Sysciophlebia balteata–Spiloblattina homigtalensis a.z. (insect zone 7), known from the lower Goldlauter Formation in the TFB (letter F in Fig. 2 and Fig. 6), was also discovered in the Ričany Formation, Boskovice Furrow (Czech Republic) (Schneider and Werneburg 2012). There, at the top of the underlying formation, a tuff yielded a CA–ID–TIMS U–Pb zircon age of 298.9 ± 0.4 Ma (Opluštil et al. 2017; Fig. 6). This age stands approximately for the Sysciophlebia balteata–Spiloblattina homigtalensis assemblage zone, which had been established in the TFB section by Schneider and Werneburg (2012).

The Goldlauter Formation comprises amphibian zone 5, i.e., Branchierpeton reinholdiApateon flagrifera flagrifera–Melanerpeton eisfeldi assemblage zone of Werneburg in Schneider et al. (2020) (letters F and G in Figs. 2 and 6). There are no precise radiometric ages for this biostratigraphic zone (Fig. 6).

An important amphibian biozone occurs in the uppermost part of the Oberhof Formation (amphibian zone 7 in Schneider et al. 2020, formerly zone 8 in Schneider and Werneburg 2012, Fig. 8) termed the Melanerpeton pusillum–Melanerpeton gracile assemblage zone, found in the Wintersbrunn fossiliferous horizon [syn. Upper Protriton Horizon (letter I in Figs. 2 and 6)]. The evidence of this biozone is widespread. It is reported from the Niederhäslich Formation (Döhlen Basin, Saxony), the Schöneck Formation (Hessian Basins), the Humberg Horizon (top of Meisenheim Formation, Saar–Nahe Basin), from several French basins and from the Olevětín Member (Broumov Formation, Intrasudetic Basin) as well (Fig. 6). In the Intrasudetic Basin, this biozone is found in the Ruprechtice limestone horizon in the middle part of the Olevětín Member. Below, at the top of the Nowaruda Member, there occurs the Vraní hory Rhyolite, which is dated by the CA–ID–TIMS U–Pb method at 297.1 ± 0.3 Ma (Opluštil et al. 2016). This is in good agreement with samples 2 and 3 from the Oberhof Formation and sample 1 from the Rotterode Formation of the TFB (Fig. 6). Approximately, 70 m of clastic deposits was accumulated between the Vraní hory Rhyolite and the Ruprechtice Horizon.

The amphibian zone 7 is situated in a median level between the uppermost Older Rhyolite of the Oberhof Formation represented by samples 2 and 3 and the Komberg Rhyolite of the Rotterode Formation (sample 1). Hence, the age of this biozone can be estimated at about 296 Ma from our data.

Conclusions

We report the first investigation by the U–Pb CA–ID–TIMS method of a Stephanian to early Permian (lower Rotliegend) profile in the Thuringian Forest Basin. We draw the following conclusions.

  1. 1.

    First, the results fulfill the expectation that the ages change in the same order as the lithostratigraphic sequence. That is, formations following one another in the order given by lithological mapping yielded zircon ages from older to younger in the same order. Sample 4 is an exception that can be explained by the local tectonics of this sample site.

  2. 2.

    New geochronological results for the Ilmenau Formation (sample 5) and the Oberhof Formation (samples 2 and 3) indicate an age difference of 2.5 Myr between both formations, indicating a ca. 2.5 Myr quiescent magmatic period for the predominantly sedimentary Manebach and Goldlauter Formations. This is in accordance with different biostratigraphic zones. However, still-existing discrepancies between biostratigraphy and geochronology have to be resolved by further geochronological dating and additional paleobiological data.

  3. 3.

    The zircon ages of sample 5 (Kickelhahn Rhyolite, Ilmenau Formation) and the two rhyolite samples (samples 6 and 7, Möhrenbach Formation) are very close to the internationally valid chronostratigraphic Carboniferous-Permian boundary (International Chronostratigraphic Chart 2019) at 298.90 ± 0.15 Ma. Based on the new Kickelhahn Rhyolite zircon age, we assume that the position of the C–P boundary should be located at the base of the Manebach Formation. Consequently, the position of the C–P boundary in the TFB could be established more precisely compared to previous estimates.

  4. 4.

    Compared to earlier concepts, the new zircon ages indicate that the total duration of the volcanic activity in the TFB was considerably shorter. Whereas previously an age interval between 295 and 275 Ma was expected for the whole volcanic period in the TFB, now the range is reduced to a period from ca. 300 Ma (Möhrenbach Formation) to ca. 296 Ma (Rotterode Formation). Consequently, the total duration of volcanic activity is shortened from ca. 20 Myr to ca. 4 Myr. We consider the pre-treatment of the dated zircons by chemical abrasion an important procedure in removing zircon parts with strong Pb loss, e.g., due to late magmatic–hydrothermal activities, that resulted in exceedingly young ages in previous studies.

  5. 5.

    Consequently, the possible time span from the early Permian volcanic rocks (top Rotterode Formation) up to the Illawarra Reversal has to be extended to ca. 30 Myr. This result requires a new view on the evolution of the TFB. This stratigraphic interval can cover 25–30 Myr in parts of the external Variscides.

  6. 6.

    Our newly obtained ages of the dated TFB formations that contain important biostratigraphic fossils can be used as time marks for certain biostratigraphic zones in other continental late Carboniferous–early Permian basins.

  7. 7.

    In the future, the chronostratigraphic sequence in the Thuringian Forest could be even more precisely defined by dating the Inselsberg and/or Heuberg Rhyolite as a mark for the base of the Oberhof Formation, and the Elgersburg Rhyolite as a mark for the youngest rhyolite. In addition, new age data from the other continental basins are urgently needed, in particular from the Saar–Nahe Basin, the Saale Basin near Halle and the Chemnitz Basin. Such new data will allow the definition of volcanic–subvolcanic–plutonic complexes as cogenetic systems and the assessment of their geochemical signatures, thermal history, and magmatic–hydrothermal evolution in an integrated way.