Timing of magmatic-hydrothermal activity in the Variscan Orogenic Belt: LA-ICP-MS U–Pb geochronology of skarn-related garnet from the Schwarzenberg District, Erzgebirge

Here, we present in situ U–Pb laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) ages of andradite-grossular garnet from four magmatic-hydrothermal polymetallic skarn prospects in the Schwarzenberg District, Erzgebirge (Germany), located in the internal zone of the Variscan Orogenic Belt. Within the geochronological framework of igneous rocks and hydrothermal mineralization in the Erzgebirge, the obtained garnet ages define three distinct episodes of Variscan skarn formation: (I) early late-collisional mineralization (338–331 Ma) recording the onset of magmatic-hydrothermal fluid flow shortly after the peak metamorphic event, (II) late-collisional mineralization (~ 327–310 Ma) related to the emplacement of large peraluminous granites following large-scale extension caused by orogenic collapse and (III) post-collisional mineralization (~ 310–295 Ma) contemporaneous with widespread volcanism associated with Permian crustal reorganization. Our results demonstrate that the formation of skarns in the Schwarzenberg District occurred episodically in all sub-stages of the Variscan orogenic cycle over a time range of at least 40 Ma. This observation is consistent with the age range of available geochronological data related to magmatic-hydrothermal ore deposits from other internal zones of the Variscan Orogenic Belt in central and western Europe. In analogy to the time–space relationship of major porphyry-Cu belts in South America, the congruent magmatic-hydrothermal evolution in the internal zones and the distinctly later (by ~ 30 Ma) occurrence of magmatic-hydrothermal ore deposits in the external zones of the Variscan Orogenic Belt may be interpreted as a function of their tectonic position relative to the Variscan collisional front.


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 . 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 Harlaux et al. 2018; Moscati and Neymark 2020;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(Burisch et al. , 2018Roberts et al. 2020), cassiterite (Yuan et al. 2011;Zhang et al. 2017), and garnet 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.

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. 1998Dias et al. , 2002Finger 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 Kroner 2015, 2016 andreference 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 Fig. 1 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).
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 Tichomirowa et al. 2019).
A first generation of skarn is associated with the peak metamorphic event at ~ 340 Ma  and lacks significant metal endowment. Large and wellendowed 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 latecollisional granitoids . 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 ). 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 km 2 . 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) 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) 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 cm 2 ) and were set after careful screening to identify growth-zones with higher 238 U/ 206 Pb 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 206 Pb, 207 Pb, 232 Th, and 238 U 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 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 208 Pb signal after subtracting the radiogenic 208 Pb (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 semiquantitatively 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 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 grossularrich 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) 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.

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). Semiquantitative 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.

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), reddishbrown 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  (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 coarsegrained 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).

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 magmatichydrothermal 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 magmatichydrothermal 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  at Hämmerlein may record magmatichydrothermal 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 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 ). 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 magmatichydrothermal 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.

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 skarnforming hydrothermal activity in the Erzgebirge ). 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 magmatichydrothermal 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 oreformation 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).
In contrast to the internal zones, the onset of magmatichydrothermal 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 postcollisional 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).