The tungsten-gold veins of Bonnac (French Massif central): new constraints for a Variscan granite-related genesis

– In the Brioude-Massiac district (French Massif Central: FMC), a network of W-As-Bi-Au quartz veins constitutes the Bonnac deposit, where tungsten is the major economic element, together with high-grade gold (up to 15 g/t Au). The evolution of this mineralization has been divided into 3 stages: (i) an early deep-seated wolframite-löllingite stage formed between 12 to 9 km, at up to 400 ° C; (ii) a ductile/brittle deformation stage associated with scheelite and arsenopyrite deposition, with an estimated temperature of 480 – 300 ° C; (iii) a late stage controlled by ﬂ uid-overpressure potentially triggered by fault-valve mechanism, at a depth of 7 to 5 km, and a temperature estimated between 266 to 240 ° C, is marked by micro-fracturing in ﬁ lled by native bismuth, bismuthinite, hedleyite, electrum, pyrite and base-metals. Structural analysis and apatite LA-ICP-MS U/Pb dating demonstrate a spatial and temporal link between the emplacement of the peraluminous leucogranitic dykes and the Bonnac mineralization. In more details, the mineralization was deposited between 321 – 316 Ma, during, or just after, the emplacement of the peraluminous dykes estimated around 329 – 315 Ma, suggesting a magmatic-hydrothermal transition for the ore-forming process. In the proposed model, the cooling of a hidden two-mica granitic pluton could have generated a magmatic ﬂ uid, and acted as the heat source responsible for ﬂ uid ﬂ ow towards inherited permeability zones. The magmatic ﬂ uid could have then re-equilibrated at high temperature by ﬂ uid-rocks interaction. The sharp changes in pressure, associated with the decrease of the temperature, and sul ﬁ de-fugacity generated by a late input of meteoric ﬂ uid were responsible for the deposition of the late gold-stage. At the regional scale, the tungsten-gold event is ascribed to an early hydrothermal stage, dissociated from the formation of the antimony event in the district. The leucogranitic dykes and Bonnac quartz veins are controlled by a NW-SE stretching direction, interpreted as an expression of the Serpukhovian-Bashkirian syn-orogenic extension (D4 event of the FMC). These new data provide evidence for an early tungsten and gold metallogenic event in the FMC, prior the “ Or300 ” event. The genetic classi ﬁ cation of the Bonnac mineralization is equivocal. The W-As-Bi-Au-quartz veins exhibit the features of both an “ orogenic gold ” deposit at a relatively deep emplacement level (mesozonal), and an Intrusion-Related-Gold-Deposit (IRGD) type with a spatial-temporal link with the peraluminous intrusion emplacement. We propose that the Bonnac deposits represent an intermediate type between a typical orogenic-gold deposit and an IRGD. We argue that the presence of economic high-grade gold content in tungsten vein-type, and more generally the IRGD deposits, have been underestimated in the Variscan French Massif Central. régionale, l ’ événement or-tungstène est attribué à un stade hydrothermal précoce, et dissocié de la formation de l ’ événement à antimoine du district. La mise en place des dykes leucogranitiques et des veines de quartz de Bonnac est contrôlée par une direction d ’ étirement NW-SE, interprétée comme une expression de l ’ extension syn-orogénique d ’ âge Serpukhovien-Bashkirien (événement D4 du Massif Central français). Ces nouvelles données fournissent des preuves d ’ une première période métallogénique à tungstène et or dans le FMC avant l ’ événement « Or 300 ». La classi ﬁ cation génétique des minéralisations de Bonnac est équivoque. Les veines de quartz à W-As-Bi-Au présentent à la fois les caractéristiques du modèle de l ’ « or orogénique », avec un niveau relativement profond de mise en place (mésozonal), et celles du modèle Intrusion-Related-Gold-Deposit (IRGD) avec un lien spatio-temporel entre la minéralisation et une intrusion peralumineuse. Nous proposons que Bonnac soit un type intermédiaire entre un modèle typique de l ’ or orogénique, et un IRGD. Nous soutenons que la présence de teneurs économiques en or dans les gisements de tungstène de type ﬁ lonien, et plus généralement ceux de type IRGD, ont été sous-estimés dans la partie varisque du Massif Central Français.


Introduction
The French Massif Central (FMC) represents a prominent metallogenic province for gold (Bouchot et al., 2005) but only a few tungsten deposits have been found, although resources of at least 45,000 t of WO 3 have been estimated (Audion and Labbé, 2012), this is probably the consequence of an underexploration for this metal in the FMC. The ubiquity of the quartz-wolframite hydrothermal veins makes it the main source of tungsten in the FMC. Some deposits with still a strong economic potential were formerly mined, such as Leucamp, Enguialès and Echassières (Marignac and Cuney, 1999; Fig. 1).
U-Pb TIMS dating on wolframite from the Limousin, Echassières, Montredon, and St-Mélany districts ( Fig. 1) has been recently performed (Harlaux et al., 2018). The obtained U-Pb ages define the existence of three tungsten mineralizing events closely related to magmatic activity: -Visean to Serpukhovian mineralization  coeval with the emplacement ages of large peraluminous two-mica plutons (ca. 335-325 Ma) during the synorogenic extension; -Serpukhovian to Bashkirian mineralization  synchronous with the emplacement of syn-tectonic granites (ca. 315-310 Ma) during the late-orogenic extension; -Gzhelian to Asselian mineralization (around 300-298 Ma), formed during the post-orogenic extension, coeval with the metallogenic peak of the Or300 event during which the main orogenic gold deposits were emplaced in the FMC (Bouchot et al., 2005). In the central part of the FMC, i.e. in the Pontgibaud, Brioude-Massiac, and Leucamp districts, the timing for the emplacement of the wolframite bearing veins remains poorly constrained. Furthermore, the genetic relationships between the ore deposits and the granitic emplacement is not yet demonstrated. The Bonnac area, in the central part of the Brioude-Massiac district is a good target to address these questions (Fig. 1). Indeed, in this area, a dozen of W-As-Bi-Au bearing quartz veins from the South of the Bonnac village to the Scoufour hamlet ( Fig. 3) further north (Périchaud, 1970;Sandras, 1988) are spatially associated with a leucogranitic dyke swarm dated at 322 ± 7 Ma (K/Ar on muscovite, Bril et al., 1991), while the W-As-Bi-Au mineralization yields late Permian age of about 250 Ma (K/Ar on muscovite from Bonnac; Bril et al., 1991). Variable gold grades ranging from 1 to 15 g/t have been reported in ore samples from the region (Périchaud, 1970;Sandras, 1988), but the relationship between gold and tungsten ore is not clearly documented.
The aim of this study is to better understand the W-As-Bi-Au quartz vein deposits in the FMC, and more particularly in the Bonnac area. It focuses on the structural control and the timing of these ore deposits, and their relationships with the Carboniferous tectonics and magmatism, and will be integrated within the hydrothermal framework of the Variscan FMC. Detailed field structural studies, and in situ apatite LA-ICP-MS U/Pb dating were carried out. Mineralogical and geochemical studies were also performed to decipher the tungsten-gold relationships and propose an updated metallogenic model. Located in the central part of the FMC, the 1000 km 2 Brioude-Massiac polymetallic district is covered by Cenozoic basaltic flows, and bounded to the northeast by the Oligocene Limagne graben. The district that belongs to the Variscan orogen, is a N120-130°E trending antiform that refolds the Fig. 2. Location of the study area within the geological and metallogenic framework of the Brioude-Massiac district. The simplified geology modified from Lasnier et al. (1982) and Thonat et al. (2014). Location and extension of mineralized veins are from Périchaud (1970). stack of metamorphic nappes recognized in the entire Massif Central (e.g. Burg and Matte, 1978;Ledru et al., 1989;Faure et al., 2005Faure et al., , 2009Figs. 1 and 2). The core of the area is occupied by the Lower Gneiss Unit (LGU) represented by biotite-sillimanite paragneiss, intruded by the Céroux orthogneiss derived from a Cambro-Ordovician alkaline granite (Burg and Matte, 1978;Lasnier et al., 1982;Mathonnat, 1983;Thonat et al., 2014). The LGU is tectonically overlain by the Upper Gneiss Unit (UGU), which is formed by migmatitic paragneiss with rare orthogneiss, and a felsic-mafic association called the "leptynite-amphibolite complex". The mafic rocks experienced an eclogite-facies or granulite-facies metamorphism during the eo-Variscan evolution (Marchand, 1974;Lasnier, 1977;Bernard-Griffiths et al., 1980). The nappe stacking in the Haut-Allier area was probably acquired during the early Carboniferous, around 350-340 Ma (monazite U/Pb dating on syntectonic granulite; Pin and Peucat, 1986). Then two extensional events, coeval with granitic intrusions accommodated the crustal thinning widely observed in the FMC (e.g. Malavieille et al., 1990;Faure and Pons, 1991;Faure, 1995): the syn-orogenic extension at ca 330-310 Ma, called D4 event , is characterized by a NW-SE maximum stretching recorded by the syn-tectonic biotite monzogranite and the two-mica granite plutons and their country rocks. In the Brioude-Massiac area, andalusite-leucogranitic dykes and their related two-mica granites (Sandras, 1988;Bril et al., 1991) belong to this D4 event; the late-orogenic extension (D5 event, between 310 to 299 Ma; Faure et al., 2009), with a NNE-SSW maximum stretching direction, was contemporaneous with the emplacement of highly differentiated peraluminous intrusions and pegmatites, and rare-metal granites (Cuney et al., 2002). The opening of the Late Carboniferous coal basins belongs to the D5 event (Faure, 1995).

The Bonnac W-As-Bi-Au deposits
The study area is located in the central part of the Brioude-Massiac polymetallic district which contains more than 200 veins-type deposits and occurrences with variable ore content: W, Sn or W-As-Bi-Au (Bonnac, Vèze-Bosberty veins), Sb (Au), Pb-Zn or F and Ba (Fig. 2). Antimony veins have been the main economic substance mined between 1850 and 1978 with an historical production of 43 000 t Sb metal (Périchaud, 1970;De Gramont et al., 1990). The Bonnac area was only mined for gold between 1880 and 1911 with a production of only 2 kg of gold, and never for tungsten (Périchaud, 1970;De Gramont et al., 1990). The Bonnac W-As-Bi-Au quartz-veins are spatially associated with Sb (Au) quartz veins (Fig. 2) and are located in the southwestern limb of a regional N120-130°E antiform. All veins seem to crosscut the foliation of the LGU paragneiss (Figs. 2 and 3). Some quartz veins develop immediately at the contact with andalusite leucogranitic dykes (Périchaud, 1970;Sandras, 1988). The ore mineralogy, is mainly represented by arsenopyrite-löllingite with wolframite, minor scheelite, associated with bismuthinite, native bismuth, tetradymite and visible gold in a quartz gangue with muscovite clusters (Périchaud, 1970;Bril, 1982;Bril and Beaufort, 1989).
Detailed studies on fluid inclusions (Bril, 1982), hydrothermal alteration (Bril and Beaufort, 1989), and lead isotopic signature (Marcoux and Bril, 1986) allow to constrain the P-T conditions of vein deposition. The Bonnac veins display two successive hydrothermal stages with an early high temperature (homogenization temperatures between 350-400°C) CO 2 -CH 4 -rich ± N 2 aqueous fluid with low-salinity (about 5 equiv wt % NaCl), interpreted as of metamorphic origin, and a late low temperature (homogenization temperatures between 200-300°C) aqueous fluid, with low-salinity, of meteoric origin, similar to the ore forming fluid from the Sb (Au) deposits. According to the same authors, the formation of the W-As-Bi-Au and Sb (Au) veins could be contemporaneous, the result of a single hydrothermal system coeval with the emplacement of an underlying, albeit hidden, Permian pluton. The ore veins were considered to show a zonal distribution, with W-As-Bi-Au occurring in a limited area, near the intrusion, the Sb (Au) veins being more distal.

The intrusive rocks
Leucogranitic dykes, striking N30°E, N60°E and sometimes N100°E, with a variable thickness (between 50 cm to 20 m), and small extension (less than 3 km), are cropping out in a few clusters in the district, the Bonnac area being one of the most important (Figs. 2 and 3). These dykes display a homogeneous mineralogical and chemical composition in the overall district (Sandras, 1988). Quartz is the most abundant mineral, associated with plagioclase, and sub-euhedral Kfeldspar (<5 mm) with biotite aggregates and muscovite. Andalusite is observed as small sized grains (0.6 to 0.8 mm), associated with the K-feldspar and quartz. A temperature emplacement of ca 600-650°C, and a pressure of 0.4 GPa, was estimated from the andalusite-quartz assemblage (Sandras, 1988). An emplacement age of 322 ± 7 Ma was documented by K/Ar radiometric dating on muscovite (Bril et al., 1991).
The dykes chemical compositions correspond to a peraluminous two-mica granite with 74 to 76% SiO 2, and a Al 2 O 3 /CaO þ Na 2 O þ K 2 O ratio up to 1.6 (see EDM. 1 for analytical details). This bulk chemistry is similar to those found in the Chatelet gold district, and other Visean to Serpukhovian leucogranites from the Limousin district (e.g. St-Sylvestre massif; Scaillet et al., 1996). A local partial melting of the LGU paragneiss could be responsible for the genesis of those intrusive rocks (Sandras, 1988). A spatial and potentially genetic link with the W-As-Bi-Au deposits of the Brioude-Massiac district is assumed. At Bonnac, the dykes are interpreted as the shallow expression of an underlying pluton (Sandras, 1988;Bril and Beaufort, 1989).

Materials and methods
The accessory mineral phases were identified using a Merlin compact Zeiss Scanning Electron Microscopy (SEM), co-operated by BRGM-CNRS-Orléans University, and equipped with Energy Dispersive System (EDS) for qualitative analyses. In addition, SEM-cathodoluminescence images have been made using MEB-CL MIRA-3 Tescan (BRGM-CNRS-Orléans University). Chemical analyses of arsenopyrite were obtained with the Cameca SX-Five Electron Probe Micro-Analyzer (EPMA) at the Institut des Sciences de la Terre d'Orléans (ISTO). Analyses were performed using an accelerating voltage of 20 kv, a beam current at 40 nA, a beam diameter of 3 mm, and a counting time of 30 s for Fe, Ni, Co, Sb, Pb, Cu, As, Zn, Ag, Bi, and 60 s for gold. Details on the method, parameters and analytical data are given in ESM2.
U/Pb dating on hydrothermal apatite from the mineralization was conducted in situ by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) at the GeOHeLiS analytical platform (Géosciences Rennes/OSUR, Univ. Rennes) using an ESI NWR193UC Excimer laser coupled to an Agilent 7700 Â quadrupole ICP-MS equipped with a dual pumping system to enhance sensitivity (Paquette and Piro, 2014). The analytical procedure followed the method described in Pochon et al. (2016). Details, standards analyses and analytical data are provided in ESM3. All errors are listed at 2s level. Concordia diagrams and calculated ages have been produced using IsoplotR (Vermeesch, 2018).
4 Field characterization of the W-As-Bi-Au quartz veins

Field and structural features from the Bonnac area
The strike of the mineralized quartz veins displays four different orientations similar to those of the leucogranitic dykes. The N30°E trend is represented by the Riolet, Androl, Rode, Borie veins), the N60°E by the Costillon one, while few veins strike N90°E and N160-170°E (Fig. 3). The vein dip varies between 30°to 45°to the SE whatever their strike. Mineralized quartz veins are located in the hangingwall or footwall of the leucogranitic dykes (Figs. 4 and 5). They display similar features all around the area, with a lenticular shape, and a narrow width of 1 to 20 cm. The host-rocks are weakly affected by the hydrothermal alteration. The geometric relationships show that the leucogranitic dykes crosscut the paragneiss foliation, and are cut by the mineralized veins (Figs. 4 and Figs. 4B and 5) with locally developed breccias (Fig. 4B).
Most of the vein wallrocks are fault planes with a down-dip striation, indicating a normal motion (Fig. 5). Some veins display also a striation with a NE 40-50°pitch, in which quartz steps indicate a normal component of displacement ( Fig. 5C, D). The observed kinematics of the Costillon vein ( Fig. 4A) also indicate a normal motion that complies with the drag folds asymmetry developed around the leucogranitic dykes (Fig. 5B). Late N60-65°E striking sub-vertical barren faults, with gouge, offset the quartz veins (Fig. 4A).

Field and structural features from the Scoufour area
The Scoufour mineralized quartz veins exhibit some differences with the Bonnac ones. Three N10-30°E striking low-angle veins with a dip between 10 and 20°to the SE are observed (Fig. 6). Their thickness is comprised between 30 and 50 cm, but may exceed 1 m. Here, leucogranitic dykes are absent, and the host-rock underwent a hydrothermal alteration related to the mineralization. A greisen-type alteration is thus observed up to 1 or 2 m from the vein, resulting in the bleaching of the paragneiss by transformation of biotite into muscovite with associated silicification. Drag folds observed near the Ferbert vein show a normal motion (Fig. 6B). A late fracturing event, represented by N20-25°E steeply NW dipping barren faults with gouge, is also recognized (Fig. 6A). All the structural features observed in the Bonnac and Scoufour areas are summarized in Figure 7. In both cases, the emplacement of the mineralized quartz veins occurred during a single deformation phase coeval with a hydrothermal event. The magmatic event responsible for dyke emplacement postdated the regional metamorphism but preceded the ore veins emplacement. Furthermore, normal faults rework the preexisting faults that controlled the emplacement of the leucogranitic dykes. Lastly, both dykes and mineralized veins are cut by NE-SW subvertical faults with subhorizontal striation indicating a sinistral motion.
5 Mineralogy and geochemistry of the W-As-Bi-Au ore veins

The mineralized facies
All the mineralized veins of the study area have similar habitus, with a massive infill of coarse white-grey quartz gangue with parallel ribbons of a few centimeters, or aggregates of arsenopyrite and coarse wolframite (Figs. 4A and 8). Secondary scheelite is common in cm-scale pseudomorphic aggregates of early wolframite that may  appear as rare relicts and rarely observed as isolated, coarse grains or infilling fractures (Fig. 8).
On the basis of gangue and ore mineral textural observations, 3 stages have been recognized, namely: stage 1 with tungsten and arsenic; stage 2 with arsenic, and tungsten; stage 3 with bismuth, gold, and base metals.

Stage 1: early tungsten and arsenic
The first opening of the veins corresponds to the deposition of a wolframite and löllingite assemblage in a macrocristalline Qz1 quartz (crystal size >500 mm; Fig. 9A, B). Observed with cathodoluminescence, Qz1 displays a blocky habitus with regular growth bands without significant deformations (Fig. 9B), and numerous fluid inclusions. Coarse euhedral wolframite crystals (up to 1 cm) are associated with cm-size aggregates of löllingite (Fig. 9A). SEM analyses of wolframite exhibits a high iron grade (6-8 wt %) and a lower manganese one suggesting a ferberite. Large apatite crystals (>100 mm) have been observed within Qz1 in close association with wolframite (Fig. 9C). Their euhedral shape and straight grain boundaries suggest it is coeval with Qz1.

Stage 2: arsenic and late tungsten
The beginning of stage 2 is marked by an intense and heterogenous episode of deformation ( Fig. 9A-D), which has affected the previous stage 1 assemblage. Qz1 grains experienced a dynamic recrystallization (Stipp et al., 2002) shown by undulose extinction, subgrain rotation and recrystallization with golfed or jagged shape of the Qz1 boundary, bulging and quartz neograin (Qz2) formation at the expense of the early Qz1 (Figs. 9B, D and 10A and B). Qz2 displays undulose extinction and grain size reduction and is located along 10 to more than 500 mm shear bands that surround the edges of Qz1, or crosscut it.
The crystallization of arsenopyrite (Fig. 9B) and scheelite are coeval with the deformation episode. Arsenopyrite forms small rhombic crystals (up to 100 mm) and fills the fracture network of the stage 1 löllingite or wolframite and sometimes totally pseudomorphs after löllingite. Xenomorphic scheelite forms at this stage by partial or total pseudomorphosis of previous wolframite (Fig. 9A, C).The scheelites, show rare inclusions of niobium oxides rich in niobium as on the Ferbert vein (65 wt % Nb 2 0 5 ), iron (25 wt % FeO), manganese (4.5 wt % MnO), tungsten (2.4 wt % WO 3 ), titanium (2.2 wt % TiO 2 ) and scandium (0.9 wt % Sc 2 O 3 ). Those uncommon oxides have been observed in the Puy-les-Vignes tungsten deposit (Harlaux et al., 2015). The 56 EPMA analyses of arsenopyrite (see Tab. 1 and ESM 2 for detail) have been compared with previous study of Bril (1983) from the Bonnac area and from other W-As-Au regional veins (Vèze and Bosberty, see Fig. 2 for localization). Arsenopyrite has a very similar composition in all deposits with As % average values of 33.4%, 32.8%, and 33.0% for Bonnac, Vèze, and Bosberty veins, respectively. Minor elements, e.g. Co, Ni, Sb, Au are under the detection limit of the method. According to Kretschmar and Scott (1976), and Sharp et al. (1985), the use of the arsenopyrite chemical compositions to assess ore formation temperatures may be considered with a good confidence level, because of the absence of minor element. Assuming the arsenopyrite-pyrrhotite equilibrium, the calculated formation temperature was around 480°C. These values are in the same range than those deduced from the H 2 O-CO 2 L1-type fluid-inclusions (homogenization temperature between 300 and 420°C; Bril, 1982).

Stage 3: Bi-Te-Au and base metals
This stage is marked by the deposition of a new discrete quartz generation (Qz3) in the fractures of arsenopyrite or löllingite (Fig. 10C, D). The medium grain-size (around 100 mm) Qz3 displays a subeuhedral shape or elongate-fibrous habitus. Its hyaline aspect is due to the scarcity of fluid and solid inclusions. This Qz3 never experienced dynamic recrystallization (Fig. 10C, D) and is coeval with a bismuth-telluride-gold-base metals assemblage (Fig.10C-E). Rhombic arsenopyrite crystals are commonly crosscut by large mm-size, subeuhedral pyrite, associated with bismuth-telluride and base-metals minerals (Fig. 11D, E). The EPMA chemical analysis of the pyrite from Bonnac veins does not show minor element.
Bismuthinite and native bismuth are the dominant minerals species of this stage 3. They occur as 2 to 80 mm inclusions aligned along the cleavage plane of löllingite (Fig. 11A), also within arsenopyrite, or wolframite (Figs. 10C, D and 11A-C), and as small-size aggregates (<500 mm) or in fractures crosscutting löllingite or arsenopyrite (Fig. 10E, F). Native bismuth is frequently observed as ovoidal inclusions within bismuthinite (Figs. 10F and 11C). Xenomorphic hedleyite (Bi 7 Te 3 ) is the main tellurium mineral, always closely associated with the bismuthinite and native bismuth. Tetrahedrite [(Cu,Fe) 12 Sb 4 S 13 ] and a bismuth-antimony-sulfide, close to the chemical composition of the horobetsuite [(Bi,Sb) 2 S 3 , have been identified as inclusions in a hedleyite aggregate (Fig. 10F). In the Ferbert vein, a late deposition of Cu-Zn-Pb minerals overprints and brecciating all the previous minerals. Chalcopyrite is the main mineral, with small tetrahedrite inclusions (50 mm). Sphalerite and minor galena complete this assemblage. This late deposition stage is associated with a Qz4 quartz generation displaying undeformed pyramidal grains (Fig. 11E).
Visible gold has been commonly observed in thin-sections from all veins of the Bonnac area. Most of the gold intergrowths with the Bi-Te minerals (Figs. 10B and 11A-C), or more rarely appears within sulfides (mainly chalcopyrite, Fig. 11D, E), or sometimes pyrite or isolated between Qz2 grains. Gold patches range between 1 mm to 25 mm in size but can reach 50 mm within sulfides. All display a xenomorphous rounded shape (Fig. 11). Rare maldonite (Au 2 Bi) has been observed by SEM analysis as a few microns inclusions intergrowths with native bismuth.
Ten EPMA analyses of gold grains intergrowths with native bismuth and within chalcopyrite from the Ferbert vein are provided in Table 1 (see EDM 2 for detailed analyses). The small size of the gold grains explains why the analyses do not reach 100%. Silver is the only element present with gold (average of 16.7 wt % Ag). "Gold" of Bonnac area is thus an electrum with a gold fineness of 830 (Morisson et al., 1991).
The three economic-stages described above are summarized in Figure 12, however, this evolution represents increments of a single progressive crystallization-deformation event.
It should be noted that, a poorly developed and uneconomic post-ore stage, macroscopically highlighted by a dark-gray quartz facies, has been only observed in the Ferbert vein. All the previous stages are brecciated and cemented by a chalcedonic quartz gangue with minor carbonates. The associated metallic-minerals are pyrite-marcasite aggregates with pyrrhotite relics, with a new generation of chalcopyrite, and sphalerite. This does not greatly affect the previous assemblages.

LA-ICP-MS U/Pb geochronology on apatite
Hydrothermal apatite from the early tungsten-arsenic stage 1 (Figs. 9C and 12) of the Androl (Bonnac area) and Ferbert veins (Scoufour area) have been carefully characterized and 15 grains have been selected for dating. Namely, 9 from the Ferbert vein (mineralized quartz vein sample BM46D) and 6 from Androl vein (mineralized quartz vein sample BM39A; ESM 3 for analytical details). The cathodoluminescence image shows growth zones without evidence complex zoning (Fig. 13A, B). Between 2 and 6 U/Pb analysis have been performed on each grain (Fig. 13C, D). Apatite grains from the Ferbert vein have variable 207 Pb/ 206 Pb ratios, ranging between 0.1849 and 0.7010. The discordant points plot along a discordia line that yields a lower intercept date of 315.6 ± 2.8 Ma (MSWD = 3.7; Fig. 13C). The data plot along a discordia line yielding a lower intercept date of 323.4 ± 3.2 Ma (MSWD = 2.4, Fig. 13D). The intersect point of the discordia  give the isotopic composition of the initial lead of the dated apatite (Fig. 13C, D). Their composition shows some difference between the Ferbert vein and the Androl vein with respectively 0.823 ± 0.005 and 0.853 ± 0.013. Those results are consistent with initial 207 Pb/ 206 Pb ratios measured in different nearby mineralized veins, and main lithological formations of the Brioude-Massiac area (Marcoux and Bril, 1986). Thus, the Fournial (Pb-Ag-Sn-As veins) have 207 Pb/ 206 Pb ratios different 0.86234 ± 0.00027 (n = 14), but that of the Sb-(Au) vein of the Céroux is closer 0.8598 ± 0.0014 (n = 1). On the other Sb-(Au) veins, the isotopic compositions of the lead remain variable within de district. The aplitic dykes close to a leucogranitic composition from the la Margeride pluton, in the southern part of the district, show the strongest correlation with the Androl vein with a 207 Pb/ 206 Pb ratio of 0.8569 (n = 1, composition recalculated at 320 Ma).

Mechanisms of ore-vein formation and chemicalthermal evolution
During stage 1, the initial vein opening, coeval with the massive crystallization of the Qz1 under low-strain, was contemporaneous with a high temperature (up to 400°C) paragenesis with wolframite and löllingite from an As and W enriched fluid with low-salinity (about 5 equiv wt % NaCl) CO 2 -CH 4 -rich ± N 2 aqueous fluid (V-type fluid inclusions in the classification of Bril, 1982).
A ductile/brittle deformation episode marks the beginning of stage 2 (Fig. 12). The first Qz1 experienced a dynamic recrystallization shown by undulose extinction, subgrain rotation recrystallization, bulging and quartz neograin (Qz2) formed at the expense of the early Qz1. Ductile/brittle deformation features argue for a deep-seated formation between 12 to 9 km (Scholz, 1988). This is supported by the range of temperature for this stage, estimated around 480°C, based on the arsenopyrite thermometry and, in agreement with the estimated temperature range (400 to 500°C) deduced from the quartz textural observations (e.g. Stipp et al., 2002). Those depth and temperature estimation are in good agreement with the estimated emplacement conditions of the leucogranitic dykes of Bonnac (see Sect. 2.3). The transformations of wolframite into scheelite, and löllingite into arsenopyrite suggest a change in fluid chemistry, still rich in As and W, and marked by an increase in sulfur fugacity and in the calcium activity. The pseudomorphic textures argue for a remobilization of the tungsten from the previous ferberite. The increase of the Ca/Fe ratio in the fluid, associated with a slightly decrease of the temperature, as predicted by Wood and Samson (2000) could possibly correspond to higher fluid interaction with the surrounding metamorphic formations. This fluid could correspond to the L1-type fluid inclusions defined by Bril (1982), with low-salinity (3.5-0 equiv wt % NaCl), lower CO 2 aqueous fluid, and trapping temperature estimated around 300 to 400°C (Fig. 14).
The onset of stage 3 (Fig. 12) was marked by brittle fracturing. The previous assemblage was micro-brecciated, and invaded by a hyaline, low-strain Quartz (Qz3). This new quartz generation was coeval with the deposition of native bismuth, bismuthinite, hedleyite, electrum, pyrite, and base-metals minerals (chalcopyrite, sphalerite, galena). The microcrack infill as well as fluid and solid-inclusions bands (Fig.10A, B) suggest fluid-assisted fracturation under a hydrostatic regime at a depth comprised between 7 and <5 km (Scholz, 1988). This stage 3 might be interpreted as a consequence of a crack-and-seal trigged by a fluid overpressure mechanism (Ramsay, 1980). The emplacement of the previous massive quartz generations could have sealed the porosity in the fault zone and act like an impermeable barrier to fluid flow allowing an increase in the fluid pressure. The crack-and-seal process has triggered pressure drop from lithostatic to hydrostatic regime, and might be interpreted as a consequence of a fault-valve mechanism (Sibson et al., 1988;Cox et al., 1991;Robert et al., 1995;Gaboury and Daigneault, 2000). The textural evolution reflects the transition from a lithostatic pressure context during stages 1 and 2 with a limited permeability to pervasive fluid circulation in connected fractures, therefore implying a transition from, an aseismic to, a seismic environment (Cox, 1987;Cox et al., 1991).
The deposition temperature of gold during stage 3 can be deduced from the paragenetic association of native bismuth and bismuthinite. According to Barton and Skinner (1979), a range of temperature between 320 to 240°C can be proposed. However, the association of native bismuth and hedleyite can only be stable below a temperature of 266°C (Dimitrova and Kerestedjian, 2006). Therefore, we propose a temperature range between 266 to 240°C for the deposition of stage 3, which is in good agreement with a ca 200°C homogenization temperature deduced from the L2-type secondary aqueous fluid inclusions (Bril, 1982). This relatively low temperature accounts for the deposition of antimony-rich Bi-Te minerals like the horobetsuite or the tetrahedrite.

Possible mechanisms for gold deposition
The common occurrence of electrum within Bi-Te minerals might emphasize the role of a polymetallic melt/liquid for the concentration of precious metal, as already reported for similar mineralization (Ciobanu et al., 2006;Cook et al., 2007;Ciobanu et al., 2010;Zachariá s et al., 2014). The polyphase character of the ore-inclusions (electrum intergrowth in native bismuth, or electrum, native bismuth intergrowths in hedleyite), the ovoidal shape of the electrum inside native bismuth, and the presence of minerals such as maldonite or horobetsuite might support this mechanism (Ciobanu et al., 2006;Zachariá s et al., 2014). But our deposition temperature estimate, between 266 to 240°C, is too close to, or even below, the bismuth solidus temperature, and the eutectic temperature of the Au-Bi system (241°C). This leads us to suggest that a large part of the gold and Bi-Te phases probably precipitated directly form a hydrothermal fluid enriched in Bi-Te-Au and were not scavenged by a bismuthrich liquid/melt.
Between stages 1 and 3, a temperature decrease is observed (Fig. 14). In addition, the stage 3 underwent a significant pressure drop associated with a change in tectonic regime, which might have triggered the input of a meteoritic fluid (Bril, 1982;Fig. 14). Furthermore, the decrease in the sulfide activity in the ore-fluid could be responsible for the gold deposition by destabilization of gold-sulfide complexes Table 1. Summary of the EPMA assays obtained on arsenopyrite from the Brioude-Massiac W-As-Bi-Au veins and gold from Bonnac. Estimated temperatures of crystallization with the arsenopyrite geothermometer are compared with the fluid-inclusions data (Bril, 1982).

Arsenopyrite composition
Number of analysis S (wt %) Std Fe Std As Std Sb Std Au Std Total As (at %) Std T (°C) (Kretschmar and Scott, 1976) Trapping temperature (°C) estimated by fluid inclusions (Bril, 1982) Asp from Bonnac W-As-Bi-Au veins (data from this study and Bril, 1983)

Tectonic control and timing of ore deposits within the Variscan evolution framework for the French Massif Central
The vein geometry observed in the field is controlled by the inherited structures. The medium dip (30-50°SE) of the Bonnac veins reused the previous fault system coeval with dyke emplacement, whereas the low-angle dip Scoufour veins are controlled by the host-rock foliation. These anisotropy surfaces can be easily reactivated during faulting, and channeled the ore-bearing fluids.
Kinematic data on the striated fault planes coeval with W-As-Bi-Au veins allowed us to derive the stress tensor by an optimized inversion method (Delvaux and Sperner, 2003). The result indicates an extensional tectonic regime with the maximum (s 1 ), intermediate (s 2 ), and minimum (s 3 ) principal stress axes oriented 39°E/86°, 37°E/04°and 127°E/0°, respectively (Fig. 15). The W-As-Bi-Au quartz veins deposition was controlled by a NW-SE horizontal extension direction, consistent with the NW-SE stretching direction documented for the D4 tectono-magmatic event in the FMC (Faure, 1995). This widespread event, related to the synorogenic Variscan extension in the FMC, dated between 325 and 310 Ma, is contemporaneous with the emplacement of the leucogranitic and monzogranitic syn-tectonic plutons (Figs. 1 and 16;Faure and Pons, 1991;Faure, 1995;Faure et al., 2009;Talbot et al., 2005aTalbot et al., , 2005bJoly et al., 2007Joly et al., , 2009. The hydrothermal apatite grains from two localities in the Bonnac area yield LA-ICP-MS U/Pb dates of 315.6 ± 2.8 for the Ferbert vein (Scoufour area) and 323.4 ± 3.2 Ma for the Androl vein (Bonnac area), respectively interpreted as the emplacement age of the W-As-Bi-Au veins. This age of ca. 320 Ma is consistent with the D4 tectonic, which is also coeval with the Serpukhovian-Bashkirian generation of tungsten ore deposits with minor gold in the FMC as the Puy-les-Vignes or St-Mélany ( Fig. 16; Harlaux et al., 2018). Closer to the Brioude-Massiac district, the Leucamp and Enguialès ore deposits in the Châtaigneraie district yielded the most important tungsten resources of the FMC (Derré, 1983;Bogdanoff et al., 1987;Demange et al., 1988;Lerouge and Bouchot, 2009) are other examples of this tungsten ore deposition event.
All previous works (Périchaud, 1970;Bril and Beaufort, 1989) in the Brioude-Massiac district, proposed a single hydrothermal event, responsible for the formation of both antimony and tungsten-gold from a perigranitic metal zonation. This model seems to be challenged by our results. Despite their spatial link with the antimony veins, the W-As-Bi-Au ones yield a distinct age, since the antimony veins of the Brioude-Massiac district belong to a younger hydrothermal event referred to as the "Or 300" event (Bouchot et al., 2005;Cheval-Garabédian, 2019). We have shown that the Bonnac W-As-Bi-Au mineralizations have a polyphased history in a continuum of time (Fig. 12), which fits with a unique deformation event, decorrelated with the Or300 event. The three mineralized stage are all trapped in quartz-veins with same and distinct structural control (D4 tectonic event), same kinematics, a same infilling chronology with a continuity from a ductile/brittle to only-brittle regime, associated with a progressive cooling of the ore-fluid (Figs. 12 and 14), and by the absence of late remobilization textures during the gold stage 3. This, is also attested by the obvious age difference.
Our results, together with the emplacement ages obtained on the gold-bearing quartz veins from the Cévennes district (Chauvet et al., 2012) argue for an earlier Bashkirian-Moscovian gold-hydrothermal event, sometimes associated with tungsten, during the D4 event. Our data allow us to propose a polyphase history in the FMC characterized by two distinct hydrothermal events carrying gold, namely: (i) a first Serpukhovian to Bashkirian, and (ii) a second, Gzhelian-Asselian, as the Or300 event, that might have re-used the structural discontinuities that controlled the early event.
Recent new data obtain on the orogenic-gold deposit of the la Bellière district (Armorican Massif) argue also for a tectonic-hydrothermal gold event earlier than the Or300 event around 340-325 Ma (Cheval-Garabédian et al., 2020).

Genetic link with the Serpukhovian leucogranites and ore deposit model for the Bonnac area
Based on our results, the genetic relationships between the mineralization and the leucogranitic dykes can also be re-assessed. The leucogranitic dykes from Brioude-Massiac, dated at 322 ± 7 Ma (Bril et al., 1991) were emplaced contemporaneously with the D4 event and are spatially related to the W-As-Bi-Au deposits. The mineralizing event that followed the intrusive rocks emplacement could be interpreted as a magmatic-to-hydrothermal evolution ( Fig. 17; Linnen and Cuney, 2005;Gloaguen, 2006) with a relatively long-time duration, around 329-316 Ma (dating of the leucogranitic magmatic event). Several lines of evidence support this interpretation as: the spatial and temporal link between dykes and mineralized veins; their formation during the D4 tectonic event; the high temperature paragenesis of stage 1 at 400-500°C, close to the 600-650°C emplacement temperature estimated for the andalusite-leucogranitic dykes; their close emplacement depths estimated around 10 km (Sandras, 1988); the reworking of previous structural discontinuities.
Throughout the Variscan Orogen, there are examples of W-Au deposits that were emplaced within a time gap of ca. 10 Ma after the intrusions of peraluminous magmas (Romer and Kroner, 2016). In France, the Puy-les-Vignes deposit (Harlaux et al., 2018), and the Salau deposit in the French Pyrenees (Poitrenaud et al., 2019), or in the Bohemian Massif, the Mokrsko deposit (Zachariá s et al., 2014), are good examples of this relation.
A two-stage evolutive model for the magmatic-hydrothermal system in the Bonnac area is proposed in a continuum of time (Fig. 17). Between 329 and 315 Ma, the emplacement of leucogranitic dykes, related to an underlying pluton, was controlled by normal faulting in the regional NW-SE extensional regime. Slightly younger from the dyke crystallizations, between 321 to 316 Ma, the deformation became more localized and contemporaneous with quartz veins emplacement that recorded a magmatic-to-hydrothermal transition (Fig. 17, stage 2). These hydrothermal fluids have reused the previous discontinuities as feeder channels.
The metamorphic or magmatic origin of the ore-forming fluid for peri-granitic ore deposits is still a matter of debate (Lang and Baker, 2001;Boiron et al., 2001;Hart et al., 2002;Bouchot et al., 2005;Hart, 2007;Zachariá s et al., 2014;Marcoux et al., 2015). In the Bonnac area, the magmatic origin of the fluid is not fully demonstrated by the existing fluid inclusion study (Bril, 1982) as they rather display a metamorphic origin (low salinity, complex chemistry enriched in CO 2 ). However, a magmatic fluid affinity is pointed out by: (i) the presence of magmatic affinity element such as the W and the Bi-Te-Au association (Lang and Backer, 2001;Blevin, 2004) and by (ii) the spatial and temporal link with leucogranitic dykes. In addition, the lead isotopic compositions measured for different mineralized veins and main lithological formations of the Brioude-Massiac area support such magmatic origin of the metals (Marcoux and Bril, 1986). The 207 Pb/ 206 Pb ratios estimated on our dating are close to those obtained on the aplitic dykes of leucogranitic composition of the la Margeride granite (see Sect. 6 for the data).
To take into account all those features, we propose a genetic link between those mineralized veins and the leucogranitic dykes, where the initial magmatic fluids were derived from the crystallization of an underlying granitic pluton. The heat source supplied by the pluton was responsible for fluid migration toward the high permeability zones. The high temperature (400-500°C) interaction between the fluid and the host-rocks (greisen type alteration in the vein hanging walls) could have re-equilibrated the fluid, overprinting its initial magmatic signature. Then faultvalve mechanisms channeled the ore-fluids into structural discontinuities (Sibson et al., 1988), and fluid overpressure associated with potential fault-valve mechanism could be responsible for vein opening and the gold-stage 3 deposition, in a shallow crustal environment.
The chemical composition of the Brioude-Massiac andalusite-leucogranitic S-type magma dykes (EDM 1) is characterized by an oxidizing redox state (Fig. 18), although a  Bril (1982), estimated temperature for quartz textural deformation are from Stipp et al. (2002). Fig. 15. Stereographic projection of faults planes contemporaneous of mineralized quartz veins with their kinematics and associated striae. s1, s2, s3: orientation of principal stress axe calculated from fault-slip data inversion using WinTensor software in Schmidt's lower hemisphere equal-area projection (Delvaux and Sperner, 2003).  (Béziat et al., 1980;Bouchot et al., 2005;Harlaux et al., 2018). White boxes correspond to U-Pb dating on wolframite by TIMS method, from Harlaux et al. (2018). Black boxes correspond to Ar/Ar dating on muscovite from Monié et al. (1999Monié et al. ( , 2000. For the orogenic-gold deposits, le Châtelet dating (K/Ar on illites) is from Bouchot et al. (2005); from Nicaud (2001) for Saint-Yrieix district (K/Ar on illite); from Chauvet et al. (2012) for the Cévennes district (Ar/Ar on muscovite). District name abbreviations are the same as in Figure 1. Geological map modified after Chantraine et al. (1996). (B) The Bonnac hydrothermal event is replaced in the different geodynamic events of the Variscan orogen in the French Massif Central  reduce redox state is more commonly observed. These characteristics are considered as more favorable than a reduce redox state to form a tungsten ore deposit. Indeed the early magnetite precipitation allows the preservation of Au in the melt fraction of granitic magmas (Blevin, 2004).
7.5 New ore deposit model for the Bonnac deposits?
Because of the metamorphic origin of fluids, the Bonnac deposits was previously classified as an orogenic gold deposit (Bouchot et al., 2005), although they do not fulfill all the necessary criteria (Groves et al., 1998).
The redox state of the leucogranitic dykes that spans the boundary of the ilmenite series (Fig. 18) is also a favorable context for IRGD deposits (Lang and Baker, 2001;Blevin, 2004).
On the other hand, the structural evolution from a ductile/ brittle deformation to a brittle one, fault-valve mechanism features, and the presence of a late meteoric fluid input, are classical features found for the orogenic gold type deposits (Boiron et al., 1990(Boiron et al., , 2003Groves et al., 1998;Bouchot et al., 2005). Therefore, the Bonnac ore-zone represents a mineralization type intermediate between the orogenic-gold and IRGD models as already demonstrated for the Mokrsko gold deposit in the Bohemian-Massif for example (Zachariá s et al., 2014).

Conclusions
The Bonnac veins show a W-As-Bi-Au paragenesis that results from polyphase hydrothermal and magmatic events. Three stages of ore deposition have been identified starting with wolframite and löllingite in a deep-seated setting, between 12 to 9 km and up to 400°C. The second ductile/ brittle deformation stage was responsible for scheelite and arsenopyrite deposition at 480-300°C. The late stage was marked by brittle micro-fracturing infilled with native bismuth, bismuthinite, hedleyite, electrum, pyrite and Cu-Zn-Pb minerals at 266 to 240°C and 7 to <5 km depth. This could be the result of fault-valve mechanism. The decrease of the temperature-pressure and the sulfide activity caused by ongoing sulfide precipitation seems to be the major factors responsible for the gold deposition during the late stage of deposition.
The mineralized veins are spatially and temporally linked to the emplacement of the andalusite-leucogranitic dykes of the Brioude-Massiac district around 329-315 Ma. Our new LA-ICP-MS U/Pb dating on hydrothermal apatite from the deposit documents an emplacement age around 320 Ma. This genetic link suggests a magmatic-to-hydrothermal transition for the Bonnac ore-forming process. In the newly proposed model, the cooling of a peraluminous pluton might have supplied magmatic fluids, and heat source for mineralized fluid flows toward high permeability zones such as previous discontinuities. The magmatic origin of the fluid has been partly overprinted by fluid-rocks interaction.
The emplacement of both the andalusite-leucogranitic dykes and quartz veins was controlled at the regional scale by a NW-SE extensional direction, consistent with the tectonomagmatic D4 event corresponding to the syn-orogenic extension in the FMC Variscan orogen. This geodynamic event is older than the "Or 300" one, suggesting therefore an earlier fertile metallogenic period at ca. 320 Ma for tungsten and gold deposition in the FMC.
The genetic classification of the Bonnac mineralization is equivocal. The W-As-Bi-Au-quartz veins exhibit the features of both orogenic gold (relatively deep level of emplacement) and IRGD model (spatial-temporal link between the mineralization and the peraluminous intrusion). We propose that it represents an intermediate type between the typical orogenicgold and IRGD model.
The Serpukhovian to Bashkirian peraluminous andalusiteleucogranitic dyke swarms appear as a favorable metallotect for future mining exploration in the FMC. The widespread occurrences related to this type of magmatism in the FMC suggest that other Bonnac-like, or IRGD-like, deposits probably exist. IRGD deposits, and the economic high-grade gold content in tungsten vein-type, might have been underestimated in the FMC.

Supplementary material
ESM1: Geochemical compositions of the Brioude-Massiac leucogranitic dykes replaced in the A) A-B plot, and B) Q-P plot from Debon and Le Fort (1988) diagrams. Data from Sandras (1988). ESM2: Electron microprobe analyses of arsenopyrite and the gold grain from the Bonnac veins. ESM3: Operating conditions for the LA-ICP-MS equipment and LA-ICP-MS U/Pb analyses for the apatites from the studied samples.