Introduction

Tin and W deposits are commonly associated with highly evolved ilmenite-series (peraluminous) granitoid rocks in collisional zones (Lehmann 2020); they occur in a number of mineralization styles, including vein-type, skarn, carbonate replacement, and greisen deposits. Major Sn-W ore provinces include the Bolivian Sn belt, the Southeast Asian Sn belt, the Western Tasmanian Sn province, and also the Acadian-Variscan-Appalachian Sn-W belt of northeast America and Europe (Hong et al. 2017; Richards 2011; Romer and Kroner 2016). The timing of Sn and W mineralization within the orogenic cycle, however, remains insufficiently understood in many cases. The internal zones of the Variscan Orogenic Belt of central and western Europe are an excellent example as they are host to several prolific Sn-W provinces. It has been widely assumed that formation of magmatic-hydrothermal ore deposits in these provinces is constrained to a narrow interval of  ≤ 15 Ma (e.g., Bouchot et al. 2005; Cuney et al. 2002; Lefebvre et al. 2019; Romer et al. 2007; Zhang et al. 2017) coinciding with intense magmatism during post-collisional extension (Förster et al. 1999). This is in contrast to several recent studies suggesting that magmatic-hydrothermal Sn-W-polymetallic mineralization formed episodically over a time span of more than 40 Ma (Burisch et al. 2019b; Harlaux et al. 2018; Moscati and Neymark 2020; Romer and Kroner 2015, 2016; Timón-Sánchez et al. 2019; Zhang et al. 2019).

These contrasting views on the timing and intervals in which magmatic-hydrothermal systems develop in collisional zones can mainly be attributed to the challenges associated with the direct dating of hydrothermal mineralization; greisen, skarn, and vein-type mineralization does not generally contain minerals suitable for conventional isotope geochronology. Furthermore, extensive hydrothermal overprinting may lead to systematic underestimation of the age of ore formation.

Recent advances in U–Pb LA-ICP-MS geochronology of carbonates (Burisch et al. 2017, 2018; Roberts et al. 2020), cassiterite (Yuan et al. 2011; Zhang et al. 2017), and garnet (Burisch et al. 2019b; Deng et al. 2017; Gevedon et al. 2018; Seman et al. 2017; Wafforn et al. 2018) provide a time- and cost-efficient toolkit to age-date hydrothermal mineralization. Garnet in particular is an ubiquitous skarn mineral, usually occurring early in the paragenetic sequence (Meinert et al. 2005). Garnet crystallization, hence, can be assumed to capture the onset of hydrothermal activity. This information may be used to correlate skarn formation with geodynamic processes, such as discrete magmatic-hydrothermal events. Despite the increasing number of studies that apply U–Pb LA-ICP-MS garnet dating, the precise mechanism of U incorporation into andradite-grossular garnet (grandite) still remains insufficiently constrained (cf. Guo et al. 2016; Rák et al. 2011; Smith et al. 2004). However, andradite-rich zones in grandite tend to be favorable for U–Pb LA-ICP-MS geochronology (cf. Gevedon et al. 2018; Seman et al. 2017; Wafforn et al. 2018), which may be attributed to a coupled substitution mechanism promoting the incorporation of U into andradite (Rák et al. 2011).

In this contribution, we present in situ U–Pb LA-ICP-MS geochronological data of andradite-grossular garnet from four major skarn prospects in the Schwarzenberg District, Erzgebirge (Germany). Complemented with published geochronological data, the timing and episodes in which magmatic-hydrothermal activity in the Erzgebirge occurred is constrained. Results are matched and critically evaluated together with available age data for magmatic-hydrothermal Sn and W mineralization in other parts of the internal (and external) zones of the Variscan Orogenic Belt.

Geological setting

The European Variscides

The Variscan Orogenic Belt in Europe is part of the Acadian-Variscan-Appalachian belt of northeast America and Europe. The closure of the Rheic ocean and subsequent continental collision of Gondwana and Laurussia from 390 to 330 Ma (Kroner et al. 2008) was followed by post-kinematic extension and crustal reorganization from ~ 330 to 280 Ma (Dias et al. 1998; Kroner and Romer 2013). The external zones of the European Variscan Belt (e.g., Rhenohercynian Zone, Cornwall, Cantabrian Zone, South Portuguese Zone) form an extensive fold-and-thrust belt of Palaeozoic low-grade (meta-)sedimentary rocks. In contrast, the internal zones (e.g., Central Iberian Zone, Amorican Massif, French Massif Central, Moldanubian Zone, Saxothuringian Zone) are characterized by the juxtaposition of low- and high-strain domains and regional low- to high-grade metamorphism (Dallmeyer et al. 1997; Faure et al. 2009; Kroner et al. 2008; Kroner and Romer 2013).

Magmatic activity is documented throughout the entire orogenic cycle of the Variscan Orogenic Belt, but is especially prominent during the late- to post-collisional stages (330 to 280 Ma), forming mostly ilmenite-series granitoid intrusions (e.g., Dias et al. 1998, 2002; Finger et al. 1997; Förster et al. 1999; Kroner and Romer 2013). Tin-W-polymetallic greisen, skarn, and vein-type mineralization in the internal and external zones of the European Variscan Belt is vaguely constrained to intervals from ~ 335 to 310 and 300 to 275 Ma related to late- and post-collisional magmatism, respectively (Romer and Kroner 2015, 2016 and reference therein). This includes the economically most significant Sn-W provinces, e.g., Cornwall, the Iberian Massif, and the Erzgebirge (Elsner 2014), but almost certainly also Europe’s largest operating tungsten mine, the Felbertal deposit in Austria (Kozlik et al. 2016). However, the precise timing of Sn-W-polymetallic ore formation remains rather poorly constrained for most of the districts.

The Erzgebirge

The Erzgebirge/Krušné hory Province is an erosional window stretching northeast-southwest over ~ 145 km along the northern margin of the Bohemian Massif. Located in the Saxothuringian Zone, it forms part of the internal zone of the European Variscan Belt (Fig. 1a; Kroner et al. 2008; Romer and Kroner 2015). It is predominantly composed of high-grade gneisses in the east and southeast of the province, whereas low- to medium-grade mica schists, metacarbonates, metavolcanics, phyllites, and quartzites dominate the western and northwestern section with peak metamorphic conditions reached at ~ 340 Ma (Kröner and Willner 1998). Following rapid exhumation and erosion (Kröner and Willner 1998; Schmädicke et al. 1995), the metamorphic basement was intruded by voluminous late-collisional peraluminous granites (~ 325 to 310 Ma; Förster et al. 1999; Förster and Romer 2010; Tichomirowa et al. 2019). However, there are also some smaller magmatic stocks (e.g., Aue-Schwarzenberg granite suite) that are less fractionated and locally show ages > 330 Ma (Tichomirowa and Leonhardt 2010). Post-collisional (sub-)volcanic rocks occur from ~ 305 to 275 Ma and are associated with crustal reorganization and the onset of Permian rifting (Hoffmann et al. 2013; Nasdala et al. 1998; Seckendorff et al. 2004).

Fig. 1
figure 1

A Schematic map of the Variscan Orogenic Belt in western and central Europe (modified from Harlaux et al. 2018). Symbols mark the location of significant Sn-W deposits in the internal and external zones of the European Variscan Belt. Study area marked by red square (not to scale). B Simplified geological map of the Schwarzenberg District (modified from GeoSN 2020, with additional information from Anglo Saxony Mining Ltd. 2020 and Schuppan and Hiller 2012). Shaded areas mark the extent of mineralized skarn bodies as a surface projection. Stars indicate sample localities for this study. Garnet U–Pb LA-ICP-MS ages are provided as black boxes (this study) and red boxes (Burisch et al. 2019b) for each skarn. Reference coordinate system is ETRS89 UTM zone 33N. C: Cornwall; AM: Armorican Massif; RHZ: Rhenohercynian Zone; STZ: Saxothuringian Zone; MZ: Moldanubian Zone; ERZ: Erzgebirge; BM: Bohemian Massif; FMC: French Massif Central; IM:, Iberian Massif; ASGS: Aue-Schwarzenberg granite suite; EIB: Eibenstock granite; GER: Germany; CZ: Czech Republic

The Erzgebirge is host to abundant greisen, skarn, and epithermal vein-type deposits genetically related to late- and post-collisional magmatic-hydrothermal activity in conjunction with the aftermath of the Variscan orogeny. Published geochronological data related to magmatic-hydrothermal mineralization spans ages from ~ 326 to 276 Ma (cf. Burisch et al. 2019b; Lefebvre et al. 2019; Ostendorf et al. 2019; Romer et al. 2007; Zhang et al. 2017). The time window in which Sn-W(± Mo) greisen deposits formed has previously been constrained to 326 to 318 Ma (Lefebvre et al. 2019; Romer et al. 2007; Zhang et al. 2017), which coincides with the emplacement of late-collisional highly fractionated peraluminous Li-mica granites (Förster et al. 1999; Tichomirowa et al. 2019).

A first generation of skarn is associated with the peak metamorphic event at ~ 340 Ma (Burisch et al. 2019b) and lacks significant metal endowment. Large and well-endowed Sn-Zn-In skarns have previously been associated with two distinct temporal stages. Skarns of the older stage formed at ~ 325 to 313 Ma, overlapping closely with the age range reported for greisen deposits of the Erzgebirge; these are thought to be affiliated with the intrusion of late-collisional granitoids (Burisch et al. 2019b). In contrast, a younger stage of skarn formation in the Erzgebirge was dated to ~ 308.1 ± 3.6 to 295.5 ± 4.4 Ma, coinciding with post-collisional magmatism (Burisch et al. 2019b). These ages overlap within error with the 276 ± 16 Ma epithermal Ag-Pb-Zn-veins of the Freiberg District that have tentatively been related to post-collisional magmatic activity (Bauer et al. 2019b; Burisch et al. 2019a; Ostendorf et al. 2019; Swinkels et al. 2021) and that represent the youngest record of magmatic-hydrothermal activity in the Erzgebirge.

The Schwarzenberg District in the western Erzgebirge (Fig. 1b) hosts prolific Sn-W-polymetallic skarn mineralization associated with laterally extensive meta-carbonate horizons (Lorenz and Hoth 1967) over an area of about 180 km2. The main skarn bodies in the Schwarzenberg District include the Globenstein W-Sn skarn as well as the Sn-Zn-In skarns of Breitenbrunn, Hämmerlein, and Antonstal (Fig. 1b). Total inferred resources reported for the district are 59.75 Mt at 0.39 wt% Sn, 15.85 Mt at 0.42 wt% WO3, 60.62 Mt at 1.01 wt% Zn, and 49 Mt at 70.36 g/t In (cf. Anglo Saxony Mining Ltd. 2020; Hösel 2003; see Table ES1 for more details). The skarns vary significantly in size, metal endowment, and mineralogy, which has recently been attributed to differing spatial positions relative to their magmatic-hydrothermal fluid source(s) (Burisch et al. 2019b; Korges et al. 2019; Reinhardt et al. 2021). Skarns are mainly hosted by marble horizons within the meta-sedimentary sequences; however, metasomatism of gneisses and schists may also be observed (e.g., Hösel 2003; Hösel and Haake 1965; Lorenz and Hoth 1967; Meyer and Röthig 1967; Reinhardt et al. 2021). The skarns occur within distances of < 100 to > 1000 m to the closest granite contact (cf. Hösel 2003; Schuppan and Hiller 2012). In the northern sector of the district, they tend to be smaller in volume and are characterized by a pyroxene-dominated mineralogy with disseminated base metal sulfides (Hösel and Haake 1965; Hoth et al. 2010), typical for skarns in a distal position relative to their causative intrusion (Meinert et al. 2005; Reinhardt et al. 2021). In contrast, skarns of the southern sector tend to be significantly larger in volume, are more proximal to a possibly causative intrusion, and are characterized by a higher abundance of garnet. Skarns of the southern sector comprise disseminated to massive lenses of magnetite that may reach up to several meters in thickness (e.g., Breitenbrunn, Globenstein, Hämmerlein skarns; cf. Hösel 2003; Hösel and Haake 1965; Hoth et al. 2010; Schuppan and Hiller 2012), disseminated to massive base metal sulfides (e.g., near monomineralic sphalerite bands up to several dm thick; Kern et al. 2019), and disseminated cassiterite and/or scheelite (e.g., Globenstein; Hösel 2003).

The paragenetic sequence of the skarns commonly starts with clinopyroxene (diopside-hedenbergite-johannsenite), grandite garnet, and magnetite followed by amphibole, epidote, quartz, chlorite, fluorite, and carbonate. Base metal sulfides (dominantly sphalerite and lesser galena, arsenopyrite, chalcopyrite, pyrite, löllingite etc.) as well as cassiterite and scheelite typically occur late in the paragenetic sequence and are thus mainly accompanied by quartz, fluorite, and chlorite (e.g., Hösel 2003; Hösel and Haake 1965; Kern et al. 2019; Reinhardt et al. 2021; Schuppan and Hiller 2012). One exception is the Globenstein skarn where scheelite also occurs intergrown with garnet without, e.g., chlorite and fluorite (see more details below). The skarns in the Schwarzenberg District may contain several generations of garnet. Garnet color, composition, and garnet/pyroxene ratio may vary spatially for a single magmatic-hydrothermal event as well as for consecutive skarn stages. This complexity, sometimes missing stages, and often intense retrograde alteration (Korges et al. 2019; Reinhardt et al. 2021) renders it challenging and sometimes even impossible to precisely correlate each garnet generation across hand specimens.

Methods

In total, 28 skarn samples containing euhedral to subhedral garnet were cut to size and mounted in 1-inch epoxy rounds. Uranium-Pb garnet age data was subsequently collected on a ThermoFisher Element 2 sector field ICP-MS coupled with a Resonetics RESOLution S-155 193 nm ArF Excimer laser equipped with a two-volume Laurin Technic ablation cell at the Frankfurt Isotope & Element Research Center (FIERCE) of Goethe-Universität Frankfurt following the methodology described in Burisch et al. (2019b). Measurement points for each polished section derive from a small area (< 1 cm2) and were set after careful screening to identify growth-zones with higher 238U/206Pb ratios (Fig. 2) before each analytical session and to avoid superficial exposure of mineral inclusions by monitoring trace element ratios during screening. Samples were then analyzed in situ in fully automated mode overnight during two subsequent analytical sessions, sequence 1 and sequence 2. During the measurements, the signals of 206Pb, 207Pb, 232Th, and 238U were acquired by peak jumping in pulse counting mode with a total integration time of c. 0.1 s, resulting in 400 mass scans. The raw data was corrected offline deploying an in-house Microsoft Excel spreadsheet program (Gerdes and Zeh 2006, 2009). Mali grandite (dated by TIMS at 202.0 ± 1.2 Ma; Seman et al. 2017) was used as a primary U–Pb age standard; garnet standards Lake Jaco (Seman et al. 2017) and MaliGUF (in-house reference material) were employed as secondary quality control standards. Common Pb correction has not been applied to the data due to its high variability although common Pb content was determined using the 208Pb signal after subtracting the radiogenic 208Pb (Millonig et al. 2012). Garnet U–Pb ages are calculated using linear regression in Tera-Wasserburg Concordia diagrams (Tera and Wasserburg 1972) since the andradite-grossular solid solution incorporates a mixture of non-radiogenic and radiogenic Pb formed due to in situ decay of U. Garnet U–Pb dates are defined as the lower intercept with the Concordia curve as determined by linear regression of discordant arrays. Reported ages are based on multiple concordant analyses from one or more individual garnets and were calculated using Isoplot 3.71. All uncertainties are reported at the 2σ level. After LA-ICP-MS analyses, all analyzed garnet samples were carbon-coated and the garnet composition of each sample was determined semi-quantitatively by SEM–EDX using a FEI Quanta 600F scanning electron microscope (SEM) equipped with two Bruker X-Flash EDX detectors located at the Institute of Mineralogy of the TU Bergakademie Freiberg to confirm petrographic observations. An acceleration voltage of 25 kV, a beam current of 10 nA, and a beam diameter of 5.8 µm were used during semi-quantitative EDX measurements.

Fig. 2
figure 2

Back-scattered electron (BSE) images of laser ablation spots in representative garnet samples of the Schwarzenberg District. Laser ablation pits related to sample pre-screening are highlighted in red color and laser ablation pits of actual measurements (including spot analysis number) are indicated in green color. A Laser ablation pits in andraditic (Fe-rich) zones of garnet aggregates from the Hämmerlein skarn (Pöh-33A). B Laser ablation spots in the relatively andradite-rich rim of a garnet crystal from the Globenstein skarn (LfULG-13B). Note that measurement attempts in the more grossular-rich core have not been successful for this particular sample. C Zoned euhedral grandite garnet from the Breitenbrunn skarn (9890B). D Laser ablation spots in andradite-rich rims of grandite from the Hämmerlein skarn (Pöh-3C)

Results

Eighteen garnet samples (609 individual spot analyses) were successfully dated by U–Pb LA-ICP-MS. Although sample pre-screening identified andradite-rich growth zones as most suitable for analysis, there seems to be no clear trend between andradite component and U content of garnet (Fig. 2). Semi-quantitative garnet composition, U and Pb concentrations, and calculated U–Pb ages are provided in Table 1. Detailed descriptions of the samples, sample localities, and associated GPS coordinates are summarized in electronic supplement Table ES2. Uncorrected isotopic data may be found in Table ES3 of the electronic supplement. A histogram of the garnet age distribution in the Schwarzenberg District, analytical errors of individual samples, and Tera-Wasserburg Concordia diagrams of selected samples are shown in Fig. 5. An overview of measurement points on all samples, histograms of the garnet age distribution of individual deposits, and an overview of all Tera-Wasserburg Concordia diagrams are provided in Figures ES4, ES5, and ES6 of the electronic supplement, respectively. The available garnet ages of the entire Erzgebirge are summarized in Fig. 6. Although several generations of garnet may occur on the deposit scale, each of the 1-inch round mounts, that were used for geochronology, does only contain one garnet generation.

Table 1 Laser-ablation U–Pb ages of garnet from the Schwarzenberg District

Antonsthal

The Antonsthal skarn is dominated by fine-grained diopside-rich clinopyroxene accompanied by pale red to brown subhedral to euhedral grossular-rich grandite garnet (Fig. 3a). Sphalerite, galena, cassiterite, and lesser chalcopyrite, pyrite, and arsenopyrite occur paragenetically late and are associated with chlorite, quartz, fluorite, and calcite (Fig. 4a). Uranium and Pb contents in grandite vary from 0.09 to 118.09 µg/g and 0.07 to 15.39 µg/g, respectively. Three garnet samples yield overlapping lower intercept ages of 312.6 ± 3.9 Ma, 319.6 ± 5.2 Ma, and 320.3 ± 3.2 Ma (Fig. 5b, c).

Fig. 3
figure 3

Representative hand specimens and related epoxy mounts (~ 2.5 cm diameter) of samples from the four analyzed skarn prospects of the Schwarzenberg District. A Fine-grained pyroxene skarn (sample 10622), drift 9 ½-1, level 9 ½, shaft 235, Antonsthal. B Coarse-grained garnet-pyroxene skarn (sample SC-3–2), “westlicher Überhaun,” west field, St.-Christoph mine, Breitenbrunn. C Coarse-grained garnet skarn (sample Pöh-3), cross-drift 2, 590 m level, adit Pöhla, Hämmerlein. D Coarse-grained garnet-pyroxene skarn (sample + 350–1197), drift 103A, 350 m level, shaft 24, Globenstein. Px: pyroxene; Grt: garnet; Amp: amphibole; Sp:, sphalerite; Sch: scheelite

Fig. 4
figure 4

Simplified paragenetic sequences for the A Antonsthal, B Breitenbrunn, C Hämmerlein, and D Globenstein skarns. The paragenetic positions of grandite garnet and ore minerals (cassiterite, scheelite, and sulfides) have been highlighted in green (garnet) and orange (cassiterite, scheelite, sulfides). The paragenetic sequences have been compiled from Gärtner (2019), Hösel (2003), Hösel and Haake (1965), Hoth et al. (2010), Schuppan and Hiller (2012) and have been complemented by petrographic observations conducted during this study. Note that the paragenetic sequences are strongly simplified and that the skarns may exhibit distinct local mineralogical variations. Px:, pyroxene; Grt:, garnet; Fsp: feldspar; Amp: amphibole; Ep: epidote; Ves:, vesuvianite; Mag:, magnetite; Hem: hematite; Chl: chlorite; Qtz: quartz; Fl: fluorite; Cal:, calcite; Prh: prehnite; Ser: sericite; St: stokesite; Lu: ludwigite; Cst: cassiterite; Ssch:, scheelite; Sulf: sulfides

Fig. 5
figure 5

A Histogram of the garnet age distribution of skarns in the Schwarzenberg District. Garnet ages cluster around three maxima at ~ 340–335 Ma, ~ 325–315 Ma, and ~ 305–300 Ma (this study; Burisch et al. 2019b). B Distribution of individual sample ages and associated uncertainties for the Schwarzenberg District (n = 25 samples). The range of geochronological data related to granite stocks of the Aue-Schwarzenberg suite (~ 334 to 322 Ma) and the Eibenstock Li-mica granite (326 to 311 Ma) are indicated as boxes above the garnet ages (Tichomirowa et al. 2019; Tichomirowa and Leonhardt 2010). CF Selected Tera-Wasserburg Concordia diagrams of skarn-related garnet LA-ICP-MS analyses from the C Antonsthal, D Breitenbrunn, E Hämmerlein, and F Hämmerlein Globenstein skarn prospects. All data-point error ellipses represent 2σ uncertainties. Calculated garnet ages are reported as lower intercept ages. MSWD: mean square of weighted deviates; n: number of analyses. PG:, Pöhla-Globenstein; AT:, Antonsthal; BB: Breitenbrunn; PH: Pöhla-Hämmerlein; GB: Gelbe Birke

Breitenbrunn

The Breitenbrunn skarn comprises two major skarn bodies that can be traced for more than 4 km along strike (Hösel and Haake 1965). The early-stage mineralogy is dominated by pale green clinopyroxene (diopside-hedenbergite), reddish-brown fine- to medium-grained grandite garnet (Fig. 3b, 4b), and magnetite, which may form massive lenses up to 5 m in thickness (Hösel and Haake 1965; Hoth et al. 2010). The intensity of retrograde alteration is relatively weak and may include variable proportions of amphibole, epidote, vesuvianite, quartz, chlorite, prehnite, and carbonate (Fig. 4b). Disseminated cassiterite and sulfides are invariably associated with the retrograde stage (Fig. 4b), with sulfides being dominated by sphalerite and galena. Dump sample 9890B (Table 1) from the southern sector of the Breitenbrunn skarn (Fortuna mine) consists of pale green euhedral coarse-grained garnet pre-dating magnetite and sphalerite. However, due to the lack of field evidence, the paragenetic position of this garnet relative to other garnet generations at the Breitenbrunn skarn remains unconstrained. The chemical composition of garnet varies from low to high andradite proportions within the grandite solid solution. Uranium and Pb contents range from 0.06 to 158.47 µg/g and 0.05 to 119.94 µg/g, respectively. Reddish-brown fine-to medium-grained garnet from pyroxene-garnet skarn shows overlapping lower intercept ages between 317.3 ± 6.3 Ma and 313.5 ± 4.7 Ma (Fig. 5b, d). The pale green garnet sample does not fall into this range and yields a considerably older lower intercept age of 337.5 ± 3.3 Ma.

Hämmerlein

The Hämmerlein-Tellerhäuser skarn complex comprises three skarn units (Hämmerlein, Dreiberg, and Zweibach, cf. Bauer et al. 2019a) that extend laterally for more than 5 km (see Fig. 1b). Samples studied here are invariably related to the Hämmerlein unit (cf. Schuppan and Hiller 2012). Skarn is mainly composed of clinopyroxene (diopside-hedenbergite), pale red to green euhedral grandite garnet and epidote, and fine-grained skarnoid domains that are composed of feldspar, clinopyroxene, and epidote (Fig. 4c). Intense retrograde alteration is recognized, which is characterized by early amphibole and magnetite, followed by late-stage chlorite, quartz, fluorite, and carbonate (Fig. 4c). Disseminated cassiterite and sphalerite as well as minor chalcopyrite, arsenopyrite, and galena accompany late-stage chlorite and quartz (Kern et al. 2019; Schuppan and Hiller 2012). The late-paragenetic position of ore minerals is furthermore supported by relatively low homogenization temperatures (185 to 220 °C) in sphalerite (Korges et al. 2019; Schuppan and Hiller 2012). Locally, sphalerite may form almost monomineralic bands that comprise late-stage coarse-grained euhedral pale brownish-orange to greenish grandite garnet (cf. Kern et al. 2019; Schuppan and Hiller 2012; Fig. 3c). Garnet varies strongly from andradite- to grossular-dominated compositions and has U and Pb contents that vary from 0.05 to 6.93 µg/g and 0.03 to 2.79 µg/g, respectively. Lower intercept ages of pale red to green garnet from pyroxene-garnet skarn with metasomatic or skarnoid textures are in the range of 327.3 ± 8.4 Ma to 324 ± 7.4 Ma, whereas pale brownish to greenish garnet intergrown with sphalerite yields younger ages of 304 ± 5.7 to 294 ± 8.3 Ma (Fig. 5b, e).

Globenstein

The Globenstein prospect consists of two marble horizons with variable degrees of metasomatism. The deeper marble horizon hosts two intensely metasomatized skarn layers (Lager 3 and Lager 4), whereas less intense metasomatism affects the hanging-wall marble horizon (Lager 5; Hösel 2003).

The prograde skarn assemblage at the Globenstein skarn is mainly composed of clinopyroxene (diopside-hedenbergite) and medium- to coarse-grained euhedral pale green to brown grandite garnet (Figs. 3d, 4d). Prograde skarn is often replaced by amphibole, epidote, vesuvianite, quartz, fluorite, and chlorite (Fig. 4d). Magnetite occupies a paragenetic position between prograde and retrograde skarn mineral assemblages and forms semi-massive to massive lenses of up to several meters in thickness (9 m maximum; Hösel 2003). Cassiterite, scheelite, and sulfide minerals are mainly associated to chlorite, fluorite, and quartz of the retrograde stage (Fig. 4d). However, scheelite may also occur as coarse-grained aggregates which are intergrown with garnet. Although scheelite seems to be paragenetically younger than garnet, the conspicuous absence of associated retrograde minerals such as chlorite and fluorite with this scheelite type suggests that it may pre-date the retrograde assemblage (cf. Hösel 2003). In the samples dated in this study, scheelite is associated with quartz and forms aggregates interstitial to garnet and pyroxene. Garnet compositions at Globenstein span a wide range of andradite-grossular proportions. Uranium and Pb contents vary from 0.05 to 6.15 µg/g and 0.03 to 18.06 µg/g, respectively. Lower intercept ages of all analyzed samples overlap within error and vary from 338.2 ± 2.5 Ma to 331.2 ± 3.0 Ma (Fig. 5b, f).

Discussion

Timing of skarn formation in the Schwarzenberg District

The broad range of garnet ages in the Schwarzenberg District indicates that skarn formation during the Variscan orogenic cycle is not constrained to a narrow window, but occurs episodically over a broad interval from ~ 338 to 294 Ma (Fig. 5, Table 1). This result is in excellent agreement with the regional study of Burisch et al. (2019b). The four major skarn bodies of this district — Globenstein, Breitenbrunn, Anthonsthal, and Hämmerlein — show distinctively different ages and age ranges, suggesting that they do not originate from the same fluid-source/intrusion.

Five garnet samples from the Globenstein skarn have, within their errors, very similar ages (338.2 ± 2.5 Ma to 331.2 ± 3.0 Ma) suggesting that at least the prograde evolution of this W-dominated skarn is related to a single magmatic-hydrothermal event occurring soon after the regional peak metamorphic event (341 ± 1 Ma; Kröner and Willner 1998). Textural, geochemical, and mineralogical arguments, including the coarse-grained isometric habit of the garnet and the observation that garnet overprints metamorphic features of the host rock (e.g., Hösel 2003; Lorenz and Hoth 1967), strongly support a magmatic-hydrothermal origin of the Globenstein skarn. Although scheelite and cassiterite post-date the prograde skarn stage, it is likely that they are related to the retrograde stage of the same magmatic-hydrothermal event, since available geochronological data and petrographic observations do not suggest a polyphase evolution of the Globenstein skarn. Nevertheless, the presented garnet ages constrain the prograde evolution at Globenstein and, hence, have to be regarded as a maximum age for the retrograde stage and associated W-Sn mineralization (see Fig. 4d).

The age obtained for the Globenstein W-Sn skarn can tentatively be linked to granitic stocks of the Aue-Schwarzenberg suite nearby, with ages of 334 ± 12 to 329 ± 3 Ma (U–Pb zircon ages; Tichomirowa and Leonhardt 2010). The garnet ages of the Globenstein skarn presented here thus reflect the oldest record of magmatic-hydrothermal fluid flow in conjunction with the Variscan orogeny in the Erzgebirge (cf. Burisch et al. 2019b; Förster et al. 1999; Zhang et al. 2017), yet the significance of these ages for the timing of Sn-W mineralization in the area is not yet fully explored.

Although located in the immediate vicinity of the Globenstein skarn, the Breitenbrunn, Antonsthal, and Hämmerlein skarn have a distinctly different metal tenor (Zn, Sn, In) and show a more complex distribution of garnet ages. This observation is used here as evidence that these three skarns may record a series of successive and temporally distinct (prograde) magmatic-hydrothermal events.

Available garnet ages suggest the presence of at least three distinct hydrothermal events at the Breitenbrunn skarn. Most of the garnet ages related to this skarn fall into the range of 317 to 313 Ma, which correlates well with U–Pb zircon ages of the proximal Eibenstock granite (~ 315 to 314 Ma; Tichomirowa et al. 2019; Tichomirowa and Leonhardt 2010) and also agrees well with garnet ages recorded for the Antonsthal skarn. Nevertheless, one garnet age of 337.5 ± 3.3 Ma records an older magmatic-hydrothermal event that coincides with the age of the Globenstein skarn. Two previously published garnet ages of 308.1 ± 3.6 and 302 ± 2.7 Ma for the Breitenbrunn skarn (cf. Burisch et al. 2019b) provide evidence for a third, younger magmatic-hydrothermal event, which post-dates the emplacement of all known granitic intrusions in the area.

Garnet ages obtained for the Hämmerlein skarn prospect are again distinctly different. The paragenetically oldest garnet, which pre-dates the ore stage (Kern et al. 2019), has ages ranging from 327.3 ± 8.4 to 324.3 ± 7.4 Ma, overlapping both with the youngest U–Pb zircon ages of the Aue-Schwarzenberg granite suite (334 to 322 Ma; Tichomirowa et al. 2019; Tichomirowa and Leonhardt 2010) and the oldest magmatic records of the emplacement of the Eibenstock granite (~ 315 Ma). One single garnet age of 312.8 ± 3.4 Ma (Burisch et al. 2019b) at Hämmerlein may record magmatic-hydrothermal activity related to the main emplacement stage of the Eibenstock granite (Tichomirowa et al. 2019). The majority of garnet ages from this deposit (304 ± 5.7 to 294 ± 8.3 Ma), however, are related to paragenetically late garnet associated with massive sphalerite bands. These garnet ages are significantly younger than all known granite intrusions in the Schwarzenberg District. Yet, they coincide with post-collisional volcanic lithologies in the larger region that lack exposed plutonic equivalents (c.f. Burisch et al. 2019b; Hoffmann et al. 2013; Kroner and Romer 2013; Seckendorff et al. 2004). Therefore, the main episode of fertile skarn formation at Hämmerlein is tentatively related to this late, post-collisional magmatic event.

The broad range of skarn-related garnet ages in the Schwarzenberg District indicates that skarn formation in the district has been polyphase and not restricted to a narrow time window. The formation of metal-endowed skarns is by no means related exclusively to the emplacement of late-collisional Li-mica granites as often assumed in previous studies (cf. Förster et al. 1999; Lefebvre et al. 2019; Romer et al. 2007; Zhang et al. 2017). Instead, skarn formation in the Schwarzenberg District occurred episodically within a broad time interval of up to ~ 40 Ma.

Timing of skarn formation in the Erzgebirge

In their regional study across the Erzgebirge, Burisch et al. (2019b) identified skarns of regional metamorphic origin associated with peak regional metamorphism at ~ 340 Ma. However, these skarns are typically barren and are, hence, not discussed further here. Burisch et al. (2019b) also delineated two main windows of fertile skarn formation across the Erzgebirge, at ~ 325 to 313 Ma and ~ 308 to 295 Ma. These were related to late- and post-collisional magmatic activity (Burisch et al. 2019b). Most of the geochronological data presented in this study for fertile skarn bodies of the Schwarzenberg District are in excellent agreement with these results. However, 338 to 331 Ma garnet dates related to the Globenstein W-Sn skarn provide evidence for a previously unrecognized magmatic-hydrothermal episode. In the following, we use available garnet age data and compare it to published geochronological data of igneous rocks and hydrothermal mineralization from the Erzgebirge to define temporal windows and geological environments in which fertile skarn formation occurred during the Variscan orogenic cycle.

Stage I

Garnet ages of > 335 Ma in the Schwarzenberg District are restricted to the Globenstein and Breitenbrunn skarns (Fig. 5b, Table 1) and reflect the oldest record of magmatic-hydrothermal fluid flow associated with the Variscan orogeny in the Erzgebirge (cf. Burisch et al. 2019b; Lefebvre et al. 2019; Romer et al. 2007; Zhang et al. 2017). Oldest garnet ages from the Globenstein and Breitenbrunn skarns overlap within error with the peak of regional metamorphism in the Erzgebirge at 341 ± 1 Ma (Fig. 6d; Kröner and Willner 1998), but textural evidence as well as the close proximity to granites of similar age support that magmatic-hydrothermal skarns locally already formed in the immediate aftermath of peak metamorphism. Despite the remaining uncertainty regarding the age of scheelite in these skarns, the narrow range of garnet ages and the anomalous W-endowment of the Globenstein prospect compared to all other skarns known in the Erzgebirge may reflect an early fertile skarn stage with particular exploration significance for W.

Fig. 6
figure 6

Skarn formation in the context of the geotectonic and magmatic evolution of the Erzgebirge. Relative frequency of A late-collisional granite intrusions, B post-collisional igneous rocks, and C skarns in the Erzgebirge. D Summary of skarn-related garnet ages within the geotectonic framework of the Erzgebirge. The age ranges of skarn stages I to III are indicated. n:, number of analyzed samples. Sources: Burisch et al. 2019b; Förster 1999; Förster and Romer 2010; Gehmlich et al. 2001; Hoffmann et al. 2013; Hofmann et al. 2009; Hübner et al. 2020; Kováříková et al. 2007; Kröner and Willner 1998; Lefebvre et al. 2019; Luthardt et al. 2018; Nasdala et al. 1998; Ostendorf et al. 2019; Quadt 1993; Romer et al. 2007; Romer et al. 2010; Seckendorff et al. 2004; Tichomirowa et al. 2019; Tichomirowa and Leonhardt 2010; Zhang et al. 2017

Stage II

Most ages of garnets in skarns of the Schwarzenberg District and across the Erzgebirge fall into the interval from ~ 327 to 310 Ma (Table 1, Fig. 6c; this study; Burisch et al. 2019b); this age range coincides with the emplacement of abundant Li-mica granites (Fig. 6a; Förster et al. 1999; Förster and Romer 2010; Hofmann et al. 2009; Kováříková et al. 2007; Tichomirowa et al. 2019; Tichomirowa and Leonhardt 2010) as well as the age of greisen-hosted Sn mineralization (Fig. 6d). Notably, the oldest garnet ages at Hämmerlein (327.3 ± 8.4 and 324.3 ± 7.4 Ma) overlap, due to their large errors, with stages I and II. However, based on petrographic evidence these ages have been related to stage II. The skarns in this age range have a Sn-Zn-(In) tenor (e.g., Bauer et al. 2019a; Hösel 1994; Kern et al. 2019; Schuppan and Hiller 2012). Due to their size and regional abundance, these skarns represent the most important targets of mineral exploration in the Erzgebirge (e.g., Anglo Saxony Mining Ltd. 2020; Burisch et al. 2019b; Hösel 1994; Schuppan and Hiller 2012).

Stage III

A significant number of garnet ages across the Erzgebirge are in the range of ~ 310 to 295 Ma (Table 1, Fig. 6c; this study; Burisch et al. 2019b) and reflect the youngest record of skarn-forming hydrothermal activity in the Erzgebirge (Burisch et al. 2019b). This skarn-forming event significantly post-dates the emplacement of voluminous granitic intrusions (cf. Förster et al. 1999; Förster and Romer 2010; Tichomirowa et al. 2019; Tichomirowa and Leonhardt 2010), but coincides with widespread post-collisional (sub)-volcanic units (Fig. 6b; this study; Burisch et al. 2019b). It is tentatively suggested that these garnet ages point to yet undiscovered (small) intrusions and related mineralization. Common trait of skarns that yield garnets of this age range, such as the Hämmerlein skarn, is their polyphase origin. They can be of considerable size and grade (e.g., Hämmerlein skarn; Anglo Saxony Mining Ltd. 2020; Hösel 1994; Schuppan and Hiller 2012), but it remains uncertain to which stage(s) of skarn formation they owe their metal endowment.

Timing of magmatic-hydrothermal mineral systems in the Variscan Belt

To gain an overview of the timing of magmatic-hydrothermal ore formation on the scale of the entire European Variscan Belt, we incorporated published geochronological data (e.g., this study; Burisch et al. 2019b; Harlaux et al. 2018; Moscati and Neymark 2020; Zhang et al. 2019) into the existing geochronological framework of the main Variscan Sn-W provinces (Fig. 7). In all internal zones of the Variscan Orogenic Belt (e.g., the French Massif Central/France, the Iberian Massif/Spain, and the Erzgebirge/Germany), the onset of magmatic-hydrothermal activity occurred simultaneously at approximately ~ 335 Ma (Fig. 7; this study; Harlaux et al. 2018; Zhang et al. 2019). This onset of magmatic-hydrothermal activity occurred less than 5 Ma after the peak metamorphic event. The youngest magmatic-hydrothermal record in the internal zones also occurred at a similar age around 295 to 275 Ma (including outliers; Fig. 7; Harlaux et al. 2018; Ostendorf et al. 2019). At least for the Erzgebirge and the Iberian Massif, the youngest ages are supported by late magmatic activity (Fig. 7). Geochronological data related to magmatic-hydrothermal deposits in the internal zones of the Variscides range from ~ 335 to 295 Ma (excluding outliers, see Fig. 7), suggesting that ore-formation occurred episodically within a time interval of at least 40 Ma. This argues against some previous models that proposed a temporally restricted occurrence of Variscan magmatic-hydrothermal ore deposits in a narrow time interval of ≤ 15 Ma related to late- and post-collisional magmatism/volcanism (cf. Bouchot et al. 2005; Cuney et al. 2002; Lefebvre et al. 2019; Romer et al. 2007; Zhang et al. 2017).

Fig. 7
figure 7

Comparison of peak metamorphism, ages of igneous units and Sn-W-polymetallic mineralization related to the Variscan orogenic cycle in internal and external zones of the European Variscan Belt. FMC: French Massif Central; blue:, peak metamorphism; pink: granitoid magmatism (and Permian volcanism); green: Sn-W-polymetallic mineralization; Eugaliès and Leucamp, post-Variscan hydrothermal overprint; Boden, metamorphic skarn (Burisch et al. 2019b); FG epithermal, Freiberg epithermal vein system (Ostendorf et al. 2019). Age data compiled from Chen et al. 1993; Chesley et al. 1993; Moscati and Neymark 2020; Neymark et al. 2018; Shail and Leveridge 2009; Warr et al. 1991 (Cornwall); Alexandrov et al. 2000; Alexandrov et al. 2002; Bertrand et al. 2001; Cheilletz et al. 1992; Cuney et al. 2002; Gebauer et al. 1981; Harlaux et al. 2018; Le Carlier de Veslud et al. 2013; Melleton et al. 2015; Scaillet et al. 1996 (French Massif Central); Carocci 2019; Chicharro et al. 2014; Chicharro et al. 2016; Dias et al. 1998; Díaz Alvarado et al. 2013; Escuder Viruete et al. 1998; Gloaguen 2006; Melleton et al. 2015; Moura et al. 2014; Orejana et al. 2012; Priem and Tex 1984; Snee et al. 1988; Timón-Sánchez et al. 2019; Zhang et al. 2019 (Iberian Massif); Burisch et al. 2019b; Förster 1999; Förster and Romer 2010; Gehmlich et al. 2001; Hoffmann et al. 2013; Hofmann et al. 2009; Hübner et al. 2020; Kováříková et al. 2007; Kröner and Willner 1998; Lefebvre et al. 2019; Luthardt et al. 2018; Nasdala et al. 1998; Ostendorf et al. 2019; Quadt 1993; Romer et al. 2007; Romer et al. 2010; Seckendorff et al. 2004; Tichomirowa et al. 2019; Tichomirowa and Leonhardt 2010; Zhang et al. 2017 (Erzgebirge)

The almost identical time span in which episodic Sn-W-polymetallic ore deposits occur in all internal zones of the European Variscan Belt suggests a broadly similar evolution of the metallogenic system during the Variscan orogenic cycle in central Europe. Episodic ore formation is associated with syn- to late-collisional, late-collisional, and post-collisional magmatism and/or volcanism.

In contrast to the internal zones, the onset of magmatic-hydrothermal activity in the external Variscan zones (i.e., Cornwall) occurred significantly later at ~ 295 to 274 Ma (Fig. 7) and appears to have a shorter duration (~ 21 Ma; Chesley et al. 1993; Moscati and Neymark 2020). The magmatism and associated ore-formation in Cornwall, thus, post-date the peak of regional metamorphism by at least 30 Ma and coincide with the youngest records of magmatic and magmatic-hydrothermal activity in the internal zones of the European Variscan Belt at 295 to 275 Ma.

Conclusions

Uranium-Pb LA-ICP-MS garnet ages from the Schwarzenberg District suggest that magmatic-hydrothermal skarn formation in the Erzgebirge occurred repeatedly and episodically. Skarns with significant (W-) metal endowment may have already formed in the immediate aftermath of peak metamorphism. Polyphase skarns, particularly those which are affected by late- and post-collisional hydrothermal overprinting, represent the most promising exploration targets in the Schwarzenberg District, since they tend to be large and exceptionally well-endowed (this study; Anglo Saxony Mining Ltd. 2020; Bauer et al. 2019a; Burisch et al. 2019b; Schuppan and Hiller 2012). From the results obtained in the Schwarzenberg District, the presence of laterally consistent (meta-) carbonate units appears to be a more important constraint for exploration for such skarn deposits in the Erzgebirge than the proximity of a known granite intrusion of a particular age or geotectonic position.

Reasons for the congruent magmatic-hydrothermal evolution that connects all the internal zones of the Variscan Orogenic Belt across central and western Europe — and the distinct later occurrence (by approximately 30 Ma) of magmatic-hydrothermal ore deposits in the external zones of the European Variscan Belt — will demand an explanation at the continental scale. We tentatively propose that the distinct difference may be related to the relative position and distance of each zone to the Variscan collisional front. In analogy to the time–space relationships documented for major porphyry-Cu belts of South America (Sillitoe and Perelló 2005), the timing of ore formation in the Variscan Sn-W provinces may be interpreted as a function of their tectonic position relative to the Rheic suture (Romer and Kroner 2015).