The perigranitic W-Au Salau deposit (Pyrenees, France): polyphase genesis of a late Variscan intrusion related deposit

– A ﬁ eld study combined with a laboratory study and 3D modeling have been performed in order to decipher the genesis of the Salau deposit W-Au mineralization (Pyrenees, France), one of the most important for tungsten in Europe. Results show the existence of two superimposed ore types, emplaced ca. 10 km depth and within decreasing temperature conditions: a calcic silicates skarn with rare scheelite and disseminated sulphides followed by a mineralized breccia with massive sulphides (pyrrhotite and chalcopyrite dominant), coarse-grained scheelite and gold, representing the main part of the ore mined in the past. This breccia is localized in ductile-brittle shear-zones which crosscut the granodiorite. U/Pb dating on zircon, apatite and scheelite, previously realized, con ﬁ rmed this polyphase evolution. These two types of mineralization, linked to the emplacement of two successive intrusions as con ﬁ rmed by sulphur isotopic analysis, granodioritic then leucogranitic, can be classi ﬁ ed as belonging to the Intrusion-Related Gold Deposit type (IRGD). The emplacement of the high-grade gold and scheelite breccia was initiated by the progressive localization of the regional deformation in the Axial Zone of the Pyrenees during the Permian within E-W dextral-reverse faults.


Introduction
During the last twenty years, new metallogenic models have been proposed, such as the Intrusion Related Gold Deposit model (Lang and Baker, 2001;Hart, 2007) or the Carlin type gold deposit model (Kesler et al., 2005), in particular thanks to the use of innovative techniques which enable to re-examine known deposits: isotopic geochemistry (Raju et al., 2015;Mao et al., 2016), geochronology (Charles et al., 2013;Ballouard et al., 2017;Harlaux et al., 2018), thermicity (Delchini et al., 2016;Poitrenaud, 2018) or 3D modeling (Cowan, 2020). The complexity of the deposits often reveals a polyphase of the mineralization, previously often described as resulting from a single mineralizing event (Lefebvre et al., 2018). Although comprising numerous historically exploited deposits and a significant mining potential, the Pyrenean segment of the Western European Variscan chain has only been the subject of very few recent studies (Link et al., 2019;Cugerone et al., 2020). This consideration is particularly true for tungsten mineralization, but also for gold, which represents the most widespread deposits of interest in the Pyrenees (Poitrenaud, 2018;Poitrenaud et al., 2019).
It is described that the Salau mineralization was formed during two stages, a prograde stage where metasomatism at the contact with the intrusion transformed carbonate rocks in calcic silicate skarns and a retrograde stage associated with "hydrothermal alteration" where the previously formed skarns are impregnated with massive sulphides (Derré et al., 1980;Fonteilles et al., 1989). Most of the scheelite appears during this retrograde stage (Lecouffe, 1987). During previous mining exploitation, the cut off grade of the ore was 0.8% WO 3 and only the "hydrothermally altered" retrograde skarn (Derré et al., 1980;Fonteilles et al., 1989), with abundant scheelite, was exploited. The exploitation of the deposit changed dramatically in the last few years of production, when mining was focused on a very rich ore type, with an average grade of 1.5% WO 3 . Gold has also been found in significant grades (5-10 g/t). Despite still significant proven resources, the mine was closed in 1986, due to the drop of tungsten prices. In comparison with other skarn mineralizations, the Salau deposit shows exceptional tungsten contents, among the highest in the world (Derré et al., 1980;Fonteilles et al., 1988). The existence of gold in tungsten mineralization, also with high grade, is a rare phenomenon in most other deposits of the same type (Meinert et al., 2005) and raises the question of the superposition of mineralizing events.
Recently, a U/Pb dating study was performed to provide time constraints on the genesis of the Salau deposit and to understand the setting of this complex mineralization (Poitrenaud et al., 2019). U/Pb datings performed on zircon, apatite and scheelite show that the granodioritic intrusion and the mineralization associated with the skarn are contemporary at 295 ± 2 Ma, while the massive sulphides and gold ore was formed later, at 289 ± 2 Ma suggesting the superposition of two mineralizing events.
In this study, and in order to complement the geochronological data recently acquired, we bring new observations and constraints for the tungsten mineralization genesis combining field survey, 3D modeling, mineralogy, microchemistry and isotopic geochemistry, at the scale of the orebodies. It appears that contrary to what was described by several authors as a simple skarn (Derré et al., 1984;Lecouffe, 1987;Fonteilles et al., 1989), the Salau deposit corresponds to the superposition of a vein breccia with massive sulphides (pyrrhotite and chalcopyrite dominant), coarse-grained scheelite and gold located in ductile-brittle shear zones intersecting a calcic silicate skarn, with rare scheelite and disseminated sulphides. The metallogenic implications of these results are discussed by repositioning the deposit in the regional late Variscan geodynamic, magmatic and hydrothermal context of the Pyrenees area.
2 Geological context 2.1 Regional geology The Axial Zone of the Pyrenees was affected by the Variscan orogeny between 330 and 290 Ma. This zone consists essentially of Precambrian to Paleozoic metasediments intruded by numerous late Variscan granitic plutons characterizing the suprastructure, and of vast gneiss and migmatites domes forming the infrastructure (Denèle et al., 2014;Fig. 1a). The Variscan deformation seems to relate to a single event during which the Axial Zone is characterized by a strong structural contrast between the suprastructure and the infrastructure. Indeed, the suprastructure mainly records a NS shortening in upper to medium crust, and the roof of the infrastructure outcropping in HT metamorphic domes, was exhumed in a context of EW stretching (de Sitter and Zwart, 1960;Zwart, 1979;Carreras and Capellà, 1994;Cochelin et al., 2018). Recent studies suggest that the deformation of these two domains appears to be largely synchronous in a dextral transpressive regime (Gleizes et al., 1998b;Mezger, 2009;Denèle et al., 2014;Cochelin, 2016;Cochelin et al., 2018). This episode is associated with a dominant calc-alkaline magmatism and the intrusion of numerous plutons, mainly of  Denèle et al. (2011); (b) Regional structural diagram of the Salau zone modified according to Colchen et al. (1997). granodioritic type (Debon and Guitard, 1996;Denèle et al., 2014).
The central part of the Axial Zone is structured by several plurikilometric shear zones, generally considered of Variscan age and possibly reactivated during the Pyrenean collision ( Fig. 1a) (McCaig and Miller, 1986;Barnolas and Chiron, 1996). However, there is no consensus about the age and the relative importance of the reactivation of these different shear zones, in particular due to the scarcity of post Variscan stratigraphic markers (Cochelin et al., 2018). Synkinematic white micas of some ductile shear zones have been dated between the Eocene and Jurassic by the 40 Ar/ 39 Ar method (McCaig and Miller, 1986;Monié et al., 1994;Wayne and McCaig, 1998;Vissers et al., 2017); their ages suggesting both a mica neoformation during Pyrenean or even Mesozoic tectonometamophic episodes but also possibly a partial reopening of the K/Ar system of inherited Variscan minerals. However, the 40 Ar/ 39 Ar ages obtained on mica cannot be used with confidence to determine the age of these faults. On one hand, the closing temperature of the 40 Ar/ 39 Ar system in mica is between 300 and 425°C (Harrison et al., 1985(Harrison et al., , 2009Villa, 1998;Cassata et al., 2009). On the other hand, the Mesozoic ages obtained in these shear zones show a partial to total rejuvenation of the isotopic systems linked to metasomatic events (250-550°C) of Permian to Cretaceous age (Boulvais et al., 2007;Poujol et al., 2010;Fallourd et al., 2014;Boutin et al., 2016). This problem about the age rejuvenation of synkinematic mica has been raised in particular in the plutons of Maladeta, Néouvielle and Canigou (Jolivet et al., 2007;Maurel et al., 2008;Metcalf et al., 2009). On the Salau deposit, a U/Pb age on apatite of 289 ± 2 Ma was recently obtained on one of these shear zones known as the Veronique fault, where most of the sulphide mineralization is hosted (Poitrenaud et al., 2019).
Late Variscan magmatism in the Pyrenees is represented by approximatively thirty plutons, mainly of calc-alkaline nature, distributed over the whole Axial Zone (Autran et al., 1970;Barnolas and Chiron, 1996). The magnetic structures of these plutons are consistent with their synkinematic emplacement during a D2b dextral transpressive episode (Evans et al., 1997;Gleizes et al., 1998bGleizes et al., , 2001Gleizes et al., , 2006Auréjac et al., 2004;Denèle et al., 2014). A close relationship between the activation of shear zones and pluton emplacement is proposed as an ascent mechanism of magma in the crust (Denèle et al., 2008). U/Pb ages on zircon dates the emplacement of most of these plutons from the Upper Carboniferous to the Lower Permian (Romer and Soler 1995;Paquette et al., 1997;Roberts et al., 2000;Maurel et al., 2008;Olivier et al., 2008;Aguilar et al., 2013). Most of the late Variscan plutonic massifs of the Pyrenees show a wide compositional variation, with a mafic external part (gabbros, diorites, granodiorite), while their core typically consist of monzogranites (Barnolas and Chiron, 1996). All of these facies are then cross cut by small leucogranitic bodies, as dyke swarms (Paquette et al., 1997;Metcalf et al., 2009). These granitoids, with an increasing differentiation, are emplaced according to several pulses, during a short period of a few My (Gleizes et al., 2006).
A single crustal thickening subdivided into 3 phases, is proposed for the Variscan segment of the Pyrenees (Denèle et al., 2009a(Denèle et al., , 2014Laumonier et al., 2010). The early D1 phase is characterized by crustal thickening during an episode of N-S to NE-SW shortening in the upper and middle crust with thrusts and folds development. This phase is possibly associated with a low to medium degree Barrovian metamorphism. According to these authors, it induces the development of a penetrative major schistosity in the middle crustal levels. A D2a syn-convergence phase of crustal scale lateral flow of the lower crust with an EW to NW-SE stretching direction, contemporary with the metamorphic peak of HT/LP reaching partial melting (Deloule et al., 2002;Maurel et al., 2008). A D2b crustal-scale folding phase leading to the formation of tight sub-vertical, multi-scale, south-verging E-W folds in the upper crust and open, kilometer-wide folds in the middle crust. The open folds of more competent Ordovician laccoliths in the middle crust led to the formation of gneissic domes (Denèle et al., 2007(Denèle et al., , 2009b. The intrusion of calc-alkaline plutons is considered contemporaneous with this dextral transpressive event (Bouchez and Gleizes, 1995;Evans et al., 1997;Gleizes et al., 1998a). The last phase of deformation D3 corresponds to the localization of the deformation during the cooling of the Variscan crust in a still NS to NE-SW general shortening context leading to the formation of EW to NW-SE multikilometric shear zones with a dextral reverse kinematic (Mezger et al., 2012;Van den Eeckhout and de Bresser, 2014).

The Salau deposit
The Salau sector, located on the northern flank of the Pallaresa anticlinorium, within Paleozoic metasediments, is bordered to the south by the Ordovician terranes of the Pallaresa zone and to the north by the North-Pyrenean fault (Fig. 1b;Colchen et al., 1997). This fault separates the Paleozoic terranes of the axial zone with the Mesozoic Aulus rift. The Salau deposit is hosted in a limestone and shale formation of the Lower Devonian known as the Port of Salau formation, with south verging recumbent folds (Bodin and Ledru, 1986). This unit is intruded by late Variscan calcalkaline granitoids (e.g. Marimaña, Bassiès or La Fourque; Fig. 1b). These plutons developed contact metamorphic aureoles associated with scheelite skarn mineralization such as Salau (Derré, 1982).
The Salau deposit is located at the contact between the small granodioritic massif of La Fourque and carbonate metasedimentary formations of Devonian age (Fig. 2). La Fourque massif includes two main intrusive facies, the apical La Fourque granodiorite on the periphery, and a porphyroid monzogranite in the core of the massif. The apical granodiorite is predominantly gray in color, with a fine-grained equigranular texture (0.5-1 mm) and forms an external halo 100 to 250 m wide, surrounding the monzogranitic core. It consists essentially of interstitial xenomorphic quartz, zoned plagioclase feldspar, sparse microcline, abundant biotite and rare zircon and apatite (Soler, 1977). This apical granodiorite thins and disappears at depth. In underground mining works, porphyroid granodiorite is the dominant facies. This facies contains microcline crystals (20%) up to 1 cm in length. Geochemical studies of Guy (1979), Raimbault (1981) and Toulhoat (1982), demonstrated the chemical composition differences between these 2 petrographic facies.
The lithological interface between Devonian marbles and the plutonic rocks is partly reactivated by shear zones, the main one is the Veronique fault, emplaced at the southern boundary of the La Fourque pluton (Lecouffe, 1987). The emplacement of La Fourque pluton initiated important metasomatic reactions at its contact. A calcic silicate skarn with scheelite mineralization (grossular, hedenbergite) was developed. The particularity of the ore mined at Salau is its high content in sulphides, in particular of massive pyrrhotite (Derré et al., 1980). According to Derré et al. (1984), the origin of the massive sulphides in the ore is linked to the circulation of hydrothermal fluids rich in tungsten, gold and copper during the retrograde stage of contact metamorphism. According to Fonteilles et al. (1989); this "hydrothermal alteration" reacted at high temperature with the previously metasomatized carbonate host rocks during the prograde stage. This episode precipitates a mineralization rich in scheelite, pyrrhotite and chalcopyrite. These authors interpreted the source of this "hydrothermal alteration" as linked to the formation of skarns and therefore to La Fourque granodiorite. Considering crosscuting relationships between La Fourque pluton and the shear zones, other authors (Ledru, 1986;Ledru and Autran, 1987), attributed the precipitation of late abundant sulphides to the development of the shear zones. These authors also proposed a different source of the fluids from that of La Fourque granodioritic pluton. These hypotheses were recently confirmed by the study of Poitrenaud et al. (2019), which shows that there is a telescoping of two mineralization types. The activation of the shear zones was initiated ca. 6 My after the emplacement of La Fourque pluton in the metasedimentary series and the development of the calcic silicate skarn. Other similar mineralizations to the Salau deposit have been described across the axial zone of the Pyrenees (Derré et al., 1984). The location of the Aurenère skarn occurrence at the west of Salau (Palau i Ramirez et al., 1997) and Bassiès skarn occurrence (Colchen et al., 1997;Poitrenaud, 2018) at the east of Salau, from which some samples were collected for comparison with the results of this study have been reported in Figure 1b. In addition, a few samples were also collected, for comparison, on the Pb-Zn SEDEX mineralization of Saube and Hoque Rabe, located east of Salau (Derré et al., 1980;Pouit, 1984;Fig. 1b).

Description of the different types of orebodies
The Salau deposit was mined from two main ore bodies, the Bois d'Anglade and the Véronique sectors. The Bois d́Anglade is located at the south-east of La Fourque pluton.

The Bois d'Anglade sector
The mineralized occurrences of the Bois d'Anglade sector are located in an embayment of the granodioritic massif (Derré et al., 1980) in the carbonate metasediments (Fig. 2). The host rocks are structured in a sub-vertical N110°E trending antiform. Most of the Bois d'Anglade mineralization is located to the east and to the west of a sub-outcropping granodiorite stock, cross cut by a of N120°E 80°N shear zones system. Exoskarns (0-30 m thickness) are developed in carbonate metasediments along the contact with the intrusion. The contact metamorphism isograds intersect the regional folds that precede the establishment of the Fourque intrusion in a largescale antiform. The massive sulphide ore is located along the Véronique fault on the southern contact of the granodiorite. Later, late faults E-W 70°N move the ore more than 100 in a reverse movement. Fig. 2. Coupe N-S du massif de la Fourque réalisée à partir des observations de terrain et des données de sondages historiques. La localisation des échantillons étudiés est reportée en bleu. Les isogrades de métamorphisme de contact recoupent les plis régionaux qui précèdent la mise en place de l'intrusion de la Fourque dans un antiforme à grande échelle. Le minerai à sulfures massifs est localisé le long de la faille Véronique sur le contact sud de la granodiorite. Postérieurement, des failles tardives E-W 70°N déplacent le minerai de plus de 100 selon un mouvement inverse.
Endoskarns are also observed at the contact with the pluton, in a relatively thin zone (< 50 cm). At the surface, the Bois d'Anglade area corresponds to an open pit of more than 50 m in height, mined by the Société Minière d'Anglade (SMA). It consists of three main levels that allow to observe the contact between the granodiorite and the host rock. The Veronique fault cannot be observed directly from the open pit, however, a few developments oriented to the south allow to identify figures of intense deformation. For example, on some outcrops in the southern part of the open pit, the granodiorite is strongly altered (silicified and invaded with massive sulphides) and fracturated according to networks of EW brittle faults steeply dipping towards the north. Moreover, on this same wall, exoskarns in contact with the intrusion are generally rich in massive sulphides (dominant pyrrhotite, but also arsenopyrite, chalcopyrite and sphalerite). The eastern part of this open pit does not show sheared granodiorite but the skarns are rich in massive sulphides. The northern part shows an undeformed contact and the skarn is not rich in massive sulphides (Fig. 3a). Accordingly, the presence of these massive sulphides seems to be strongly correlated with the deformation observed. The surrounding area mainly consists of marble with a sub-vertical N100°E stratification. At contact with the host rocks, the granodiorite is altered in endoskarn which frequently shows massive epidotization in parallel bands of few decimeters. Scheelite and pyrrhotite are rare.
In the Bois d'Anglade area, exoskarns are generally poor in scheelite and pyrrhotite (Fig. 3a). They consist of monomineral ribbons, generally dominated by pyroxenes (hedenbergite or diopside) or by the decimetric alternation of ribbons rich in pyroxene, epidote and titanite (light skarns), and of garnet, pyroxene and epidote ribbons (dark skarns). Skarns are locally cross cut by quartz, calcite and chlorite veins precipitated in tension jogs near shear zones in the granodiorite. Scheelite is commonly concentrated in the hedenbergite zone of exoskarn ( Fig. 3a), but can also be found at the border of the endoskarn. It is generally fine (< 100-200 mm) and forms ribbons in the bedding of metasomatized marbles. When the skarns have a low content in sulphides (disseminated), the scheelite is fine and also disseminated.

The Veronique fault
In underground mining works, the Véronique fault is mainly outcropping close to the orebodies of the Véronique district, at the south of La Fourque pluton. The fault has a N120°E average strike with variable dip towards the north (40°N to 80°N). The orebodies have a relatively irregular morphology and distribution in several lenses, from 1490 m to 1080 m in the area explored by mining works and drilling. The ore mined by the SMA in these orebodies had the particularity of containing very high scheelite average contents (> 1.5% WO 3 ) sometimes reaching locally more than 20% WO 3 (Poitrenaud, 2018). In the Véronique district, important gold anomalies were discovered with contents of 5 to 10 g/t Au (Fonteilles et al., 1989) with analyzes up to 16 g/t Au. Mineralization mainly occurs as, steeply dipping to the north, lenses of massive sulphides with coarse-grained scheelite hosted in skarns, marbles or in altered and mylonitized granodiorite (Figs. 3b and 3c). The surrounding granodiorite close to these lenses is laminated, silicified and contains scheelite and arsenopyrite veins (Fig. 3g).
The dextral reverse movement of the shear zone has been determined with sigmoid-shape microlithons between shear zones in the granodiorite and in the skarns intersected by the Veronique fault (Figs. 3e and 3f). Folds show the posteriority of this massive sulphides-bearing fault compared to the calcic silicate skarns (Fig. 3c). The extension of the mineralized bodies at the footwall of the Veronique fault, intersected by drilling to the northwest down-dip extension, show a reverse dextral pluri-hectometric shift between two skarn compartments, which is in agreement with observations realized by Lecouffe (1987). In La Fourque granodiorite, close to openings along faults or in breccias, subhorizontal quartz, calcite, adularia, pyrrhotite and powdery chlorite veins precipitated in tension jogs highlights the reverse movement of the Veronique fault.
The ore from the Veronique fault has a breccia texture ( Fig. 3b) composed of host rock elements (shreds of intensely folded graphitic marbles, skarns, blocks of altered mylonitic granodiorite) cemented by massive pyrrhotite and coarsegrained scheelite (100 mm-1 cm). Shear zones at contact with the Veronique fault are represented by a several meter thick band in the pluton and the host rocks. Decimetric to metric veins of laminated white quartz are commonly found in contact with granodiorite (Fig. 3c). Massive sulphides and coarsegrained scheelite are discordant on the skarns, which remain as relics (Fig. 3b). Massive pyrrhotite locally occur as a stockwork in the skarn. Quartz and calcite are in the interstices between large and ovoid scheelite crystals (> 200 mm). This assemblage is also found in endoskarns and in fractures of altered granodiorite (quartz, muscovite, epidote, calcite).
The Veronique fault has also been observed far from the La Fourque pluton, in the area of the Plagnau long ravine (1 km to the west of the deposit) and in the Pallaresa valley (10 km to the west of the deposit). In the first case, it is located at the interface between the Ordovician shales (Mont Rouch formation; Fig. 1b) and the Devonian marbles (Salau marble formation; Fig. 1b) and is marked by a 3m thick EW fault breccia steeply dipping to the north. The breccia consists of schist shreds or multi-metric marble blocks and is cemented by a fine quartz and clay matrix. The breccia is also rich in pyrite, often as decimetric lenses. At surface, the drags of parallel satellite faults close to the breccia are consistent with a dextral kinematic in the N100°E 85°S schistosity of the Ordovician schists. Despite the presence of this fault, which appears to be the lateral extension of the Veronique fault, no massive sulphide and scheelite mineralization has been observed. In the Pallaresa valley, the fault intersects the Ordovician schists along a N110°E 85°S orientation. The fault zone is marked by several veins of white quartz, rich in pyrite, separating several schist compartments, sheared and brecciated by a dextral reverse movement. Analyzes on areas rich in pyrite did not show significant gold (< 0.1 g/t) or tungsten contents (< 0.01% WO 3 ).

Relationships between mineralizations and magmatic rocks
The R1 cuddy of the Salau mine at level 1230 (Quer de l'Aigle; Lecouffe, 1987) shows boudins of dark diorite enclaves in La Fourque granodiorite and late aplitic dykes, which cross cut both La Fourque granodiorite and these dioritic enclaves. Aplite dykes are frequent at the base of the Salau deposit and near the Veronique fault where they outcrop as swarms. A cluster of coarse-grained scheelite and massive sulphide ore occurs at the border of an aplitic stockwork hosted by a granodioritic breccia (Fig. 3d). Aplite mainly consists of quartz and feldspar, with epidote alteration. Its contact with massive sulphides, marked by a border poor in feldspar and rich in quartz, is progressive and well zoned, showing their contemporaneity, with multi-millimeter ribbons of sphalerite at contact with the dyke, then galena and finally pyrrhotite and massive chalcopyrite (Fig. 3d).
Near the Veronique fault, centimetric veins of black quartz and tourmaline cross cut La Fourque granodiorite (Fig. 3h) in which they constitute a network of parallel veins resulting in a ribboned texture. They also outcrop locally as decimetric veins near the aplitic veins of the R1 cuddy at level 1230. These veins are grouped into two sub-orthogonal populations with relatively constant orientations (N80E 45°N and N170E 70°S). The tourmaline veins are parallel to the Veronique fault when they are in its surroundings. At contact with these veins, the granodiorite is silicified, exempt of biotite, and altered in a quartz, muscovite and titanite greisen (Fig. 3h). Arsenopyrite and more rarely wolframite (Lecouffe, 1987) are neoformed in silicified granodiorite (Fig. 3g) with small discordant calcite and pyrrhotite veins. Tourmaline veins evolve laterally towards the massive sulphide ore lenses of the Veronique fault.

3D modeling
To constrain the orebodies geometry at the deposit scale and the spatial relationships between skarn and massive sulphides, a 3D model of the mineralizations was constructed using the software Micromine TM (Fig. 4a). The data set used was extracted from logs of diamond and destructive drill holes realized during historical exploration and exploitation of the Salau deposit (Guiraudie et al., 1964;Passaqui and Costargent, 1965;Prouhet et al., 1966;Michard and Bouquet, 1986; internal SMA archives). In total, 677 logs of diamond drilling, 14 logs of destructive drilling and 141 air-blast holes from BRGM and SMA were collected (Fig. 4a). In addition to this drilling information, data from 2373 samples collected along exploitation faces were digitized. The length of the diamond drills varies from 25 to 492 m. Most of the analysis data relate to tungsten, expressed in % WO 3 . Only a hundred analyzes provide informations for gold and silver contents. Most of the high tungsten contents are located in the surroundings of the Veronique fault. Their use, combined with harmonized lithological logs, gives an indication of the geological structure of the deposit, which seems complex, in particular for the correlation between neighbouring drillholes. This is explained by the very irregular morphology between the granodioritic intrusion and the host rocks. 3D modeling of the mineralization geometry, required the realization of numerous serial sections to best constrain their extent.
The Veronique fault connects the Veronique district and the Bois d'Anglade, as observed in the field. To the east of the deposit, the fault cross cuts the Quer de l'Aigle area where a breccia-type mineralization with massive sulphides is outcropping. In this zone, the granodiorite has a very irregular contact with the host rocks, which explains the numerous embayments with skarn mineralization. The fault extends eastward under the Quer de l'Aigle massif. A E-W trending satellite fault crosscutting by the granodiorite was identified 50 m to the north of the Veronique fault. This satellite fault is also mineralized in massive sulphides and coarse-grained scheelite, but with a restricted thickness of 0.50 m to 1 m. At the scale of the deposit, the 3D model allows to visualize the shift generated by the fault between the orebodies of the Veronique district (Figs. 4b and 4c). The two skarn compartments identified by mining work and drilling are shifted with a dextral reverse movement of more than 250 m on each side of the Veronique fault. The fault has an average N100°E 70°N direction and is represented by a breccia cemented with massive sulphides and coarse-grained scheelite at the contact between granodiorite and host rocks.

Paragenetic synthesis
The mineralogical assemblages of the Salau deposit are complex due to the spatial and temporal superposition of two distinct mineralizations. Observation of more than 250 thin and polished sections allow to distinguish the main characteristics for each type of ore (Fig. 5).
The Bois d'Anglade ore type shows a calcic silicate and disseminated scheelite skarn, hosted by Devonian marbles. The prograde stage consists of a grossular, hedenbergite and diopside assemblage and the retrograde stage of a wollastonite, vesuvianite and fine scheelite assemblage (Fig. 6a). During retrograde hydroxylation, garnets and pyroxenes are respectively retromorphosed in epidote (clinozoisite) and amphibole (ferroactinote; Fig. 6f). Disseminated sulphides (arsenopyrite, pyrrhotite) filled the grain boundaries of calcic silicates (Fig. 6b). Crystallization of sulphides show the end of the retrograde stage.
The ore from the Veronique fault is a vein breccia which cross cuts granodiorite, skarns and marbles. Black quartz and tourmaline veins are also discordant on the granodiorite. They develop a greisen-type alteration and are locally associated with arsenopyrite. An assemblage of abundant coarse-grained scheelite and massive sulphides with dominant pyrrhotite, chalcopyrite, sphalerite, galena, electrum, native bismuth, laitakarite, native silver (Figs. 6c, 6d and 6e) systematically takes place close to N100°E 70°N fault zones (see Fig. 3c). These sulphides constitute the cement of clasts from adjacent rocks (granodiorite, marbles, calcic hornfels, skarn, skarnoids). Numerous euhedral hydrothermal apatite are associated with these massive sulphides ore. The composition of the sulphide assemblage varies from place to place and is correlated to the nature of the host rock. Graphitic levels are the richest in electrum. Quartz, chlorite and calcite sub-vertical tension jogs mark the vertical stretch induced by the reverse movement of the Veronique fault near the massive sulphide mineralization. Late assemblages with adularia, titanite, laumontite and apophyllite have also been observed, mostly in geodes.

Thermobarometric constraints
To estimate the temperature and pressure conditions during the emplacement of the La Fourque granodiorite, microprobe analyzes of black and white mica have been performed. All the samples come from La Fourque granodiorite but white mica samples, neoformed from plagioclases, come from E-W shear zones.
Twenty-seven analyzes for black mica were compiled from Lecouffe, 1987 (ESM Tab. 1). The structural formulas were recalculated at 22 oxygen and show that these black micas are biotites, with 55% phlogopite and 45% annite (K 0.99 Mg 1.38 Fe 1.16 Ti 0.04 Al 1.10 Si 2.90 O 10 (OH) 1.85 F 0.27 ). In these conditions, it is possible to use the geobarometer [P (kbar) = 3.03 Â (Al tot ) À 6.53] developed by Uchida et al. (2007). The range of values in Al tot varies from 3.07 to 3.37 and the range of results obtained varies from 2.75 to 3.63 kbar. The average pressure is around 3.3 kbar, indicating an emplacement depth of La Fourque granodiorite between 10 and 11 km.
The white micas analyzed (4 samples) come from the Veronique district at level 1430 (TP492) and from satellite faults in La Fourque ravine (TP627, TP652a and TP752a; ESM Fig. 1). All analyzes were performed with an ISTO CAMECA SX 100 electronic microprobe, set at 15 kV and 6 nA for silicates. The count time was 10 s for Si, Al, Ti, Fe, Mn, Mg, Ca, Na, K, Cl and F. The structural formulas were recalculated at 22 oxygen and show that these white micas are phengites (K 0.94 Na 0.03 Fe 0.12 Al 1.76 Mg 0.15 Al 0.80 Si 3.19 O 10 (OH) 2 poor in Fluorine (0.00 to 0.01% F; ESM Tab. 2). The results obtained on biotites and phengites were projected in the thermometric diagram of Monier and Robert (1986) to estimate their crystallization temperatures from Mg and Fe contents (Fig. 7a). Biotites show temperatures above 800°C. whereas the phengites only reach temperature values between 400 and 500°C; lower values in agreement with the secondary nature of these white micas, as attested by their structuraltextural position. The results on phengites have been projected in the diagram of Massone and Schreyer (1987), to determine a pressure range for the phengite equilibration or minimum pressure if they are as a secondary phase (Fig. 7b). For a range of temperatures between 400 and 500°C and an average phengitic substitution of Si = 3.2, the phengites of the EW shear zones show minimum formation pressures between 3.2 to 5.9 kbar, so an approximately depth of 10 km.
The composition of amphiboles also allows to estimate the pressure conditions of magma crystallization. Two amphibole analyzes from the La Fourque granodiorite were performed with an electron microprobe by Soler (1977). One analysis was realized on the rim and another on the core of the amphibole (ESM Tab . In these conditions, it is possible to use several geobarometers on the amphiboles developed by Hammarstrom and Zen (1986), Hollister et al. (1987), Johnson and Rutherford (1989), Thomas and Ernst (1990), Schmidt (1992) and Mutch et al. (2016). The calibration of Mutch et al. (2016) seems to be the most suitable due to the wide range of P-T conditions covered for igneous rocks. The range of values in Al tot varies from 1.26 for the rim to 1.96 for the core of the hornblende and results obtained are respectively of 2.50 and 4.97 kbar (ESM Tab. 3). The pressure estimated at the core of the hornblende at around 4.97 kbar, indicates that the magma started to crystallize around 16-17 km during their ascent in the crust before reaching the emplacement depth of La Fourque granodiorite at 2.5 ± 0.25 kbar (8 to 10 km).

Sulfur isotopes
To determine the sulfur source of the different sulphides in the two ore types (skarn and massive sulphide breccia), sulfur isotopic analyzes were performed in collaboration with the LEHNA laboratory at the Claude Bernard University of Lyon (Tab. 1). For comparison, 5 sulphide samples from the adjacent Aurenère and Bassiès tungsten occurrences as well as Saube and Hoque Rabe lead/zinc occurrences were included in the study (Fig. 1b). 4 samples containing pyrite crystals from different sedimentary formations of the surrounding area were also analyzed to assess their possible role as a source. The various sulphides analyzed (pyrrhotite, arsenopyrite, chalcopyrite, sphalerite, pyrite) were sorted manually under a binocular and approximately 5 g were crushed in an agate mortar to obtain a powder with a particle size of 100 to 300 mm. 16 samples were then analyzed by mass spectrometry (IRMS), online and in continuous flow mode with an elemental  analyzer (EA) at the LEHNA laboratory according to the protocol developed by Fourel et al. (2014). The mass spectrometer was of the IsoPime100 type and the elemental analyzer of the Pyrocube type in NCS combustion mode, marketed by the company Elementar. The particularity of the elemental analyzer used is to implement the technique of separating the gases to be analyzed, of the "purge and trap" type, particularly suitable for measuring the isotope ratios of sulfur (Fourel et al., 2014). Each sample was analyzed two to three times for quality control. For the calibration, 3 silver sulphide reference materials provided by the IAEA: IAEA-S1 (d 34 S CDT = À0.3‰), IAEA-S2 (d 34 S CDT = þ22.7‰) and IAEA-S3 (d 34 S CDT = À32.3‰) were used to frame the values of the unknown samples. The values of the isotope ratios are expressed in d 34 S relative to the CDT standard (Tab. 1). Sulfur concentrations were calculated from elemental analyzer TCD data. d 34 S ratios obtained define a wide range of values, from À0.42‰ to 18.14‰. Pyrites of the metasedimentary host rocks (Fig. 2) have the widest range with low values for the Silurian schists (À0.42‰) or the Salau marbles (4.86‰), and high values for the Mont Rouch schists (18.14‰). These results are consistent with the domain of sedimentary sulphides (Fig. 8). In the same figure, Saube (4.18‰) and Hoque Rabe (6.36‰) sphalerites fit with the volcanic sulfur source field, which is in agreement with their origin interpreted as volcano-sedimentary (Bois and Pouit, 1976;Pouit, 1984;Cugerone et al., 2018). In the Salau deposit, disseminated pyrrhotite of the calcic silicate skarn reach a value of 5.01‰. Massive sulphides from the Veronique fault show spread but close values: 4.56 to 5.12‰ (arsenopyrite), 1.22 to 4.85‰ (pyrrhotite) and 2.83 to 5.22‰ (chalcopyrite). These very similar values show that the sulfur from the sulphides of the two ore types is likely of magmatic origin, a conclusion that can be extended to the massive sulphides of Aurenère (arsenopyrite 1.36‰ and pyrrhotite 1.26‰). Finally, the massive pyrrhotite from Bassiès and the pyrrhotite from Plabidaous show higher values (7.63‰ and 7.29‰) but which still correspond to the domain of magmatic sulfur, with probably a participation of a sedimentary or metamorphic sulfur source (Rye and Ohmoto, 1974).

Polyphase evolution of the Salau deposit
In the Salau region, the mineralized occurrences studied belong to two distinct types of mineralization, a skarn and a vein breccia, with different geometry, mineralogy and grades. They occur as several mineralized bodies and depend on the host rock nature and its proximity to the Veronique fault.
The skarn mineralization corresponds to the orebodies of the Bois d'Anglade sector, located in embayments of the complex geometry between the granodiorite and metasedimentary host rocks. They consist of calcic silicates (garnet, hedenbergite), rare fine-grained scheelite and disseminated sulphides (arsenopyrite, pyrrhotite). Grades are generally below 0.5% WO 3 and are close to the detection threshold for gold (0.1 g/t). The prograde stage is marked by an assemblage of grossular, diopside, hedenbergite and fine-grained scheelite (< 100 mm). The retrograde stage is then marked by an epidotization of the grossular, which is also marked by alteration of the adjacent granodiorite and by the replacement of the clinopyroxene by ferroactinote. Disseminated sulphides (pyrrhotite, arsenopyrite) also occur in the grain boundaries of calcic silicates.
The vein breccia mineralization located along the E-W 70°N trending Veronique fault, has no calcic silicates and completely differs from the skarns sensu stricto (Meinert et al., 2005). This vein breccia consists mainly of massive sulphides with abundant coarse-grained scheelite in cement of host rocks clasts. Field criteria show that this mineralization is clearly posterior to the skarns, and cross cut both granodiorite, marble host rocks and skarns, also confirmed by U-Pb geochronology (Poitrenaud et al., 2019). Tungsten and gold contents of this vein breccia are high (> 1.5% WO 3 and 5-10 g/t Au) and correlated with massive sulphides. From a mineralogical point of view, the assemblage consists of coarse-grained scheelite (> 200 mm) and massive arsenopyrite, cross cut by a procession of dominant massive pyrrhotite, chalcopyrite, sphalerite, galena, euhedral apatite, native bismuth, electrum, bismuthinite, hessite, laitakarite, tetrahedrite and native silver.  Fig. 6. Photographies de sections polies montrant les éléments suivants : (a) pyrrhotite disséminée associée avec de la scheelite fine et des grenats grossulaires dans le skarn à silicates calciques du Bois d'Anglade ; (b) paragenèse à sulfures massifs avec un assemblage à pyrrhotitechalcopyrite-sphalérite recoupant la scheelite grossière ; (c) scheelite grossière fragmentée et cimentée par de la pyrrhotite massive ; (d) cristaux d'arsénopyrite fragmentés et cimentés par de la chalcopyrite et de la pyrrhotite massive ; (e) assemblage à bismuth natif À laïtakarite À electrum typique de la paragenèse aurifère de Salau ; (f) hédenbergite rétromorphosée en amphibole avec intercroissance de pyrrhotite disséminée dans les skarns du Bois d'Anglade ; (g) image MEB d'une apatite magmatique fragmentée contemporaine des biotites de la granodiorite de la Fourque ; (h) image MEB d'une apatite hydrothermale automorphe contemporaine de la pyrrhotite massive.  Monier and Robert (1986) diagram to evaluate the temperature of crystallization of black and white micas from their content in Fe and Mg; (b) Velde (1965) diagram modified by Massone and Schreyer (1987) to evaluate the H 2 O pressure conditions (kbar) during white mica crystallization. 3D modeling of the exploitation drillholes shows that the two skarn compartments of the Veronique district (developments at levels 1450 and 1230) are sectioned by the Veronique fault with a reverse dextral movement and an offset of more than 200 m. The kinematic, deduced from the shear criteria observed at the surface or in mining works, is consistent with the offset of the mineralized bodies in the 3D model. The Veronique fault is regionally connected to multikilometric shear corridors described throughout the Axial Zone (Cochelin et al., 2018). However, far from the Salau mine area, the fault does not contain any trace of mineralization. Thus, away from La Fourque pluton, as in Plagnau long ravine and in the Pallaresa valley, the extension of the Veronique fault is only represented by a barren pyritic breccia. The mineralization is therefore not only controlled by the regional faults but also by the proximity of the intrusion. The faults potentially localize fluid circulation while the heat provided by the La Fourque pluton at the time of emplacement is probably the source of hydrothermal circulation (Marcoux et al., 2015).
In the Salau deposit, massive sulphides are closely related to aplite dyke swarms and tourmaline veins which cross cut the base of La Fourque granodiorite, as shown by their progressive lateral transition to massive sulphides. We interpret them as constituting the apical part of a second intrusion, more differentiated and underlying the La Fourque granodiorite. Such relationships between aplitic dykes and sulphide mineralization have been described at a regional scale, in particular at Bassiès (Gerbal and Ledru, 1985) and Cauterets (Ledru, 1986), but also in Canada (Cantung; Rasmussen, 2004;Rasmussen et al., 2011) and in Russia (Lermontovskoe deposit; Soloviev et al., 2017a;Vostok-2 deposit, Soloviev et al., 2017b;Skrytoe deposit, Soloviev and Kryazhev, 2017c). The assemblage of native bismuth, bismuthinite and Pb-Bi and Cu-Bi sulfosalts shows a stronger affinity with late evolved magmas than with intermediate magmas (Derré et al., 1980;Ciobanu et al., 2010;Cepedal et al., 2013;Botros, 2015).
Datings combined with the composition of fluid inclusions (Krier-Schellen, 1988;Poitrenaud et al., 2019), indicate two different mineralizing episodes which take place successively under conditions of decreasing temperatures and changes in the fluid composition (Fig. 9). The results on magmatic zircon show that the emplacement of granodiorite, at the origin of the scheelite skarn, occurred at 295 ± 2 Ma (Poitrenaud et al., 2019). The P-T conditions of emplacement for the granodioritic intrusion of the Fourque, obtained from microprobe analyzes on black mica and amphibole in this study show a temperature above 800°C for a pressure of 2.5 to 3 kbar (8 to 10 km depth). These results for the granodiorite intrusion are consistent with the estimations for skarn formation temperature obtained by thermodynamics (505-540°C; Soler, 1977), fluid inclusions (microthermometric data: 470-530°C; Krier-Schellen, 1988) and RSCM (600-650°C; Poitrenaud, 2018). This first generation of fluid, linked to the formation of skarns, essentially consists of a highly saline H 2 O-NaCl solution (25 wt.% Eq. NaCl), with halite crystals, typical of a magmatic origin (Krier-Schellen, 1988). The high temperatures and the geological context confirmed that it is a magmatic fluid. Datings of hydrothermal apatite grains from massive sulphides ore show an age of 289.1 ± 1.9 Ma, younger than the formation of the skarn (Poitrenaud et al., 2019). Therefore, the two types of ore are separated by a minimum of about 2 Ma and a maximum of about 9 Ma. The second generation of fluid inclusions observed in the coarse-grained scheelites and quartz associated with the massive sulphides hosted by the Veronique fault show two types of solutions, one with H 2 O-NaCl identical to the fluid found in the skarns and the other with H 2 O-CO 2 -CH 4 interpreted of metamorphic origin by Krier-Schellen (1988) and Fonteilles et al. (1989). The temperatures estimated by different methods are relatively consistent: thermochemistry (450°C, Soler, 1977;327-366°C, Zahm, 1987), fluid inclusions (335-363°C for 2.4 kbar, Krier-Schellen, 1988) or microchemistry on phengites in shear zones for this study (400°C to 500°C for a pressure of 3.2 to 5.9 kbar). The pressure interval obtained on phengites is intrinsically explained by the degree of precision provided with the thermometric diagram of Monier and Robert (1986). Phengites were collected from the same type of shear zones and show temperatures of 400°C to 500°C. The impact of this temperature interval is important in the diagram of Massone and Schreyer (1987) due to the high-pressure variation for a small temperature variation, which explains the wide pressure range obtained (3.2 to 5.9 kbar).
In mineralized systems, a method to identify the source of the fluids is to evaluate the isotopic fractionation of stable elements which enter into their composition, such as sulfur, carbon or oxygen (Rye and Ohmoto, 1974;Höll et al., 1987;Ishihara and Sasaki, 1989;Coplen and Krouse, 1998;Seal, 2006;Yang and Lentz, 2010;Sharman et al., 2014;Tanner et al., 2016). The sulfur isotopic values obtained in this study have been compared with a compilation of all others bibliographical results obtained on the sulphides of the axial zone of the Pyrenees (Fig. 8). This comparison show that the sulphides from the metasedimentary host rocks (Vall de Ribes, Salau and San Bernat) have a d 34 S spread over a wide range of values (À15 to þ20‰), which fit the composition domain of sedimentary sulfur. Ag-Zn-Cu mineralizations of the Silurian black shales (Roca de Ponent, Coma Fosca, San Miquel (Pura et al., 2002) also show d 34 S values, typical of sulfur from sedimentary origin (À15‰ to À5‰), while most lead-zinc SEDEX in the axial zone (deposits of Liat, Solitaria, Font dels Lladres, Victoria, Reparadora, Preciosa Pepita, Saube, Hoque Rabe, Yenefrito, Tebarray), with a few exceptions (Margalida deposit þ10 to þ12‰, Preciosa Pepita deposit þ16‰; Cardellach et al., 1996), show values between À3‰ and þ10‰ in agreement with sulfur from volcanic origin. Pyrrhotite and pyrite of the barren skarns developed around the intrusions of Quérigut, Lacourt, Counozouls and Puyvalador (Guy, 1979;Toulhoat, 1982) have very low values (À13‰ to þ2‰). The origin of the latter can result from a mix between magmatic and sedimentary sulfur. Auriferous veins with quartz-arsenopyrite from Andorra and Vall de Ribes show values between þ5‰ and þ15‰. Their signature shows a possible metamorphic influence. Some authors described them as retrograde shear zones Ayora and Casas, 1995). Finally, the sulphides of mineralized skarns and massive sulphide breccias from the deposits and occurrences of Andorra, Aurenere, Sarraera, Escunyau, Arties, Bassiès, Salau, Roc Jalère, Balbonne, Boutadiol and Costabonne (Guy, 1979;Soler, 1990;Soler et al., 1990;Cardellach et al., 1992) show relatively low values, between À2‰ and þ8‰, which contribute to consider a sulfur from magmatic origin.

Metallogenic model
All these results converge to propose a two-stage model for the formation of the Salau deposit (Figs. 9 and 10). This model is consistent with the regional geological context of the Axial Zone and fits perfectly into the polyphase evolution described by Denèle et al. (2014) and Cochelin et al. (2018). The intrusion of La Fourque granodiorite is synchronous (295 ± 2 Ma) with most of the calc-alkaline plutons emplaced between 305-295 Ma during an N-S isoclinal folding episode of the Axial Zone. The emplacement P-T conditions are 800°C and 3.3 kbar. Calcic silicate skarns, with rare fine-grained scheelite and disseminated sulphide, were developed on its border at contact with the Devonian marble host rocks, during cooling. The tungsten contents are low (< 0.5% WO 3 ) and gold is absent. During the general cooling of the Variscan crust, the regional deformation is localized at the level of rheological contrasts (granitoid contactsurrounding and at periphery of metamorphic domes) along dextral reverse shear zones. These shear zones appear later, at 289 ± 2 Ma in Salau and Aurenère. A second mineralization (massive sulphides, abundant coarsegrained scheelite and gold) with high tungsten and gold contents (> 1.5% WO 3 and 10 g/t Au) is hosted by these faults and deposited at about 450°C and 3.2 to 5.9 kbar. Associated fluids come from a swarm of aplite dykes which constitute the apical part of a second intrusion, leucogranitic and more differentiated.
The Salau deposit therefore reflects the Carboniferous-Lower Permian evolution of the Variscan crust in the Axial Fig. 8. Compilation of d 34 S CDT assays realized for this study and for sulphides from different mineralizations through the axial zone of the Pyrenees (data compiled from Guy, 1979;Soler, 1990;Soler et al., 1990;Cardellach et al., 1992). The orange place names correspond to the sites from which the samples analyzed in this study were collected, while the black place names correspond to the sites from which the data was extracted in the existing literature. The result of the analyzes performed in this study is represented by a red box in the figure. The d 34 S CDT compositional domains of magmatic, volcanic and sedimentary sulfur were compiled from Lancelot et al. (1971), Rye and Ohmoto (1974) and Seal (2006). Fig. 8. Compilation de mesures de d 34 S CDT réalisées dans cette étude et pour les sulfures de différentes minéralisations à travers la zone axiale des Pyrénées (données compilées d'après Guy, 1979 ;Soler, 1990 ;Soler et al., 1990 ;Cardellach et al., 1992). Les noms de lieu en orange correspondent aux sites d'où proviennent les échantillons analysés dans cette étude alors que les noms de lieu en noir correspondent à des sites d'où ont été extraites les données de la littérature existante. Le résultat des analyses réalisées dans cette étude est représenté par un encadré rouge sur la figure. Les domaines de composition du d 34 S CDT du soufre magmatique, volcanique et sédimentaires sont compilés de Lancelot et al. (1971), Rye andOhmoto (1974) et Seal (2006).
Zone and show the transition from a homogeneous thickening stage of the upper crust, while the lower crust is partially molten, towards a system where the deformation is localized at the level of EW reverse dextral faults which drain W-Au fluids of more evolved magmas. At the scale of the Axial Zone of the Pyrenees, many calc-alkaline intrusive massifs are cross cut by swarms of dykes or even small stocks of leucogranites (Gleizes, 1992;Barnolas and Chiron, 1996;Gleizes et al., 1997;Roberts et al., 2000;Aleksandrov and Troneva, 2011;Olivier et al., 2016). Similar occurrences to the Salau deposit with the spatial and temporal superposition of two mineralization types have also been described for the intrusions of Bassiès (Poitrenaud, 2018) or Aurenère (Spain; Palau i Ramirez et al., 1997;Poitrenaud et al., 2019).
The Salau deposit has many similarities with the Mactung and Cantung deposits in Canada (Dick and Hodgson, 1982;Rasmussen, 2004;Rasmussen et al., 2011). The geological context is very similar between these Canadian deposits and Salau with polyphase intrusions, granodioritic and then leucogranitic in a carbonate metasedimentary series. Their sulphide sequence with pyrrhotite, chalcopyrite, arsenopyrite and sphalerite, posterior to the skarns, are identical to Salau with in particular native bismuth associated with gold (Palmer et al., 2013;Lentz et al., 2015). Their genesis is similar with an early metasomatic stage which develops skarns with a low tungsten content (0.2 to 0.5% WO 3 ) followed by a late hydrothermal stage, associated to a sulphide mineralization rich in scheelite (> 1% WO 3 ; Hodgson, 2000), linked to faults and to a more differentiated late granitic intrusion. Only the age and the geodynamic context of the intrusion differ: late Variscan for Salau and Cretaceous for the Yukon skarns.
The massive sulphide, coarse-grained scheelite and gold mineralization hosted by the Veronique fault, share numerous similarities with the Intrusion Related Gold Deposit model (IRGD; Hart, 2007). This recently proposed model groups a large diversity of mineralization morphology: quartz veins, sheeted or in stockwork in the Tighza deposit (Morocco; Marcoux et al., 2015), veins or disseminated in granite (Timbarra deposit, Autralia; Mustard, 2001), breccia (Kidston deposit, Australia; Baker and Andrew, 1991), greisen (Salave deposit, Spain; Harris, 1980) or skarn in contact with granite (Dublin Gulch deposit, Alaska; Maloof et al., 2001). The ore morphology of the Veronique fault (massive sulphide breccia) is compatible with the type of fault described in the Ryan Lode deposit (Bakke et al., 2000). The common characteristic of all these mineralizations is their magmatic origin, frequently marked by a mineral assemblage which contains metals and lithophilic elements such as arsenic, tungsten, tin, molybdenum, bismuth, tellurium, gold, silver and antimony (Thompson et al., 1999;Thompson and Newberry, 2000;Lang and Baker, 2001;Hart, 2007). The massive sulphide ore also shows very good correlations between tungsten, bismuth, tellurium and gold contents. The second generation of fluid inclusions observed in coarse-grained scheelites and quartz in massive sulphides ore from the Veronique fault is a solution with H 2 O-CO 2 -CH 4 (Krier-Schellen, 1988;Fonteilles et al., 1989), characteristic of fluids rich in CH 4 , resulting from reduced plutons of the ilmenite series (Zhang et al., 2019). Finally, Fig. 9. Chronological diagram of the genesis of the Salau ore deposit at the scale of the upper and medium crust of the axial zone of the Pyrenees. A N-S compression stage led to the formation of N110°E folds, the intrusion of the Fourque granodiorite and the genesis of the associated scheelite skarn ca. 295 Ma. A second stage, ca. 289 Ma corresponds to the intrusion of a late evolved granite and to the appearance of the deformation in shear zones during the gradual cooling of the crust leading to the formation of coarse scheelite, massive sulphides and gold mineralization in the brecciated veins. Fig. 9. Schéma chronologique de la genèse du gisement de Salau à l'échelle de la croûte supérieure et moyenne dans la zone axiale des Pyrénées. Un stade de compression N-S conduit à la formation de plis N110°E, à l'intrusion de granodiorite de la Fourque et à la genèse du skarn à scheelite associé vers 295 Ma. Un deuxième stade, vers 289 Ma, correspond à l'intrusion d'un granite évolué tardif et à l'apparition de la déformation dans des zones de cisaillement durant le refroidissement progressif de la croûte conduisant à la formation de scheelite grossière, de sulfures massifs et de la minéralisation aurifère dans les brèches filoniennes. mineralization is not only controlled by regional faults but essentially by the distance to the intrusion. All these criteria converge to classify the Salau deposit in the Intrusion Related Gold Deposit (IRGD) model.

Conclusions
Main results of our work show that: 1 two different types of mineralization coexist at Salau (Fig. 10), a skarn with calcic silicates, fine-grained scheelite and disseminated sulphides (arsenopyrite, pyrrhotite) developed around La Fourque intrusion emplaced at 800°C and 3.3 kbar at 295 ± 2 Ma, as well as a vein breccia with massive sulphide, coarse-grained scheelite and gold hosted by the Veronique fault, emplaced at 450°C and 3.2 to 5.9 kbar at 289 ± 2 Ma. These two mineralizations have mineralogical characteristics, peri-intrusion locations, morphology and emplacement in decreasing temperature conditions which link them to an IRGD model; 2 the coexistence of two types of mineralization also results from a magmatic polyphase, like the copper porphyry systems, the source of the vein breccia mineralization with high tungsten and gold content being linked to the emplacement of an evolved leucogranite which ascent in the upper crust through regional EW faults; 3 the evolution of mineralizing phenomenon in the Salau deposit reflects the dynamics of the differentiation at depth of late Variscan magmatism in the axial zone of the Pyrenees.
The polyphase evolution of the Salau deposit is fully in line with the late Variscan evolution of the Paleozoic basement of the Pyrenees and reflects the differentiation dynamics of lateorogenic Variscan, granodioritic and then leucogranitic magmatism. The granodioritic intrusion of Salau emplaced within an antiform at 3.3 kbar and 800°C during a general shortening stage due to the closure of the Ibero-Armorican arc at ca. 295 Ma (Cochelin et al., 2018). This period corresponds to the end of the metamorphic domes exhumation and is characterized by a general state of abnormally warm crust. Skarns were developed at the contact of the granodiorite during its progressive cooling. The temperature of the Paleozoic basement then gradually decrease and the deformation regime changes. The deformations previously homogeneously distributed, are then focalized within regional E-W dextral reverse shear zones developed in rheological contrasts. The periphery of metamorphic domes and calc-alkaline granitoids are particularly suitable. In these shear zones, differentiated intrusions, like leucogranites, rise in the upper crust at ca. 289 Ma as the last residual melt extracted from a deep magma chamber. A massive sulphide, coarsegrained scheelite and gold mineralization was developed in the faults which act as a drain and cross cut the earlier skarns with a plurikilometric offset under P-T conditions of 450°C and 3.2 to 5.9 kbar. The Salau deposit is a skarn over which a massive sulphide breccia, rich in tungsten and gold, is superimposed. According to numerous geochemical, genetic and structural criteria, it reflects the evolution of an Intrusion Related Gold Deposit (IRGD) model.   Lecouffe, 1987). Tab. 2. Result of microprobe analyzes of white mica from E-W trending shear zones (Veronique and satellite faults) intersecting the La Fourque granodiorite. N correspond to the number of analyzes. Tab. 3. Table of analytical results by electron microprobe of the hornblendes from la Fourque granodiorite (Soler, 1977) and the corresponding pressures calculated from the equations of Hammarstrom and Zen (1986);Hollister et al. (1987); Johnson and Rutherford (1989); Thomas and Ernst (1990); Schmidt (1992) and Mutch et al. (2016).