Maxime ; Bertauts
email: bertauts.maxime@gmail.com

Emilie ; Janots

Magali ; Rossi

Isabelle ; Duhamel-Achin

M.-C ; Boiron

Laura ; Airaghi

Pierre ; Lanari

Philippe ; Lach

Peiffert

O Adrian Pfiffner

Maria Iole Spalla

Gianni Balestro

Eduard Kissling

Olivier Lacombe

Jesus Martinez-Frias

Marie-Christine Boiron

Chantal Peiffert

Valérie Magnin

A New Alpine Metallogenic Model for the Pb-Ag Orogenic Deposits of Macôt-la Plagne and Peisey-Nancroix (Western Alps, France)

Keywords: Pb-Ag orogenic deposits, Western Alps, Briançonnais, U-Th-Pb geochronology, syn-tectonic mineralization, hydrothermal circulation, trace elements in sulfides, chlorite-phengite thermobarometry

Understanding mass transfer associated with fluids circulation and deformation in the Alpine orogeny is often complex due to common multistage crystallization. For example, in two emblematic and historic Pb-Ag deposits of the French Alps, Macôt-la Plagne (MP) and Peisey-Nancroix (PN), a sedimentary or orogenic origin is still debated. To discriminate between the metallogenic models of the two deposits, an integrative methodology combining field, microstructural, mineralogical, thermobarometrical, and geochronological data was here applied for establishing detailed Pressure-Temperature-Time-Deformation (P-T-t-d) mineralization conditions. Both deposits are located in Permo-Triassic quartzite of the External Briançonnais domain along the Internal Briançonnais Front (Internal Western Alps). The ore mainly occurs as veins and disseminated textures containing galena, pyrite, and variable content of tetrahedrite-tennantite and chalcopyrite. Quartz porphyroclasts and sulfide microstructures indicate a dynamic recrystallization of the quartzite during the main fluid mineralization episode. Chlorites and K-white micas (phengite) chemical analysis and thermodynamic modeling from compositional maps indicate an onset of the mineralization at 280 • C, with a main precipitation stage at 315 ± 35 • C and 6.25 ± 0.75 kbar. In situ U-Pb dating on monazite, cogenetic with sulfides, gives ages around 35 Ma for both deposits. The integrative dataset converges for a cogenetic MP-PN Alpine Pb-Ag mineralization during deformation in relation to the thrusting of the "Nappe des Gypses" and the Internal Briançonnais at the metamorphic peak.

Introduction

Fluids circulation associated with deformation in an orogenic context is at the origin of mass transfers, leading to the transport-deposition and remobilization of polymetallic deposits. Understanding metals' transfer and remobilization processes and scales is complex because deformation and fluids circulation are preferentially located at lithological contacts and on structures inherited from previous tectonic episodes. The French Alps present abundant polymetallic mineralizations that provide evidence of fluids circulation and significant mass transfer in the crust during successive orogenic cycles. Around 250 deposits and occurrences are identified by the BRGM (Bureau de Recherche Géologiques et Minières) in the northern French Alps ("Armada" the GIS database of French metallic occurrences; available online http://infoterre.brgm.fr; accessed on 19 July 2019), mainly as Pb-Zn-Ag (44%) and Fe-Cu (29%) mineralizations, hosted in the Variscan basement, the Carboniferous and Permo-Triassic cover, and in the Liassic limestones (Figure 1) [START_REF] Meloux | Bilan Synthétique de dix Années de Prospection[END_REF]. These polymetallic mineralizations were mainly investigated between the 1960s and 1980s, based mainly on the analysis of the structural context and petrographic analysis. Metallogenic models for most of these deposits are generally ambiguous, with no consensus on the geodynamic context that is proposed as late-Variscan (Sarrazins deposit; [START_REF] Caby | Etude Géologique et Métallogénique du Bord Interne de la Zone Houillère et de la Bordure des Schistes Lustrés Entre Modane et la Vallée étroite (Savoie, Haut Val de Susa[END_REF]), linked to Tethys Ocean opening (Macôt-la Plagne deposits; [START_REF] Rogel | Le Gisement de Plomb de La Plagne (Savoie)-Etude Géologique et Métallogénique-Vanoise[END_REF][START_REF] Rogel | Exemple de Gisement de Galène du Trias des Alpes[END_REF][START_REF] Bouladon | Contribution à une systématique des gisements de plomb et de zinc[END_REF]) and to hydrothermal circulation during the Alpine collision (Macôt-la Plagne and Le Reveyret deposits; [START_REF] Raguin | La mine de plomb argentifère de La Plagne (Savoie)[END_REF][START_REF] Waline | La mine de La Plagne[END_REF][START_REF] Cabrol | Etude Géologique et Minière du Massif des Hurtières-Chaine de Belledonne (Savoie)-Alpes Françaises[END_REF]). According to previous authors, e.g., [START_REF] Ypma | Rejuvenation of Ore Deposits as Exemplified by the Belledonne Metalliferous Province[END_REF][START_REF] Court | Etude Géologique et Métallogénique de la Basse-Maurienne (Savoie)-Alpes Françaises[END_REF][START_REF] Ochoa Alencastre | Etude Gîtologique des Minéralisations Fluorées du Massif de Rocheray, Savoie-Alpes Françaises Externes[END_REF][START_REF] Negga | Comportement de L'uranium lors des Métamorphismes tardi-Hercinien et Alpin Dans les Massifs des Aiguilles Rouges et de Belledonne (Vallorçine, Lauzière)[END_REF], several mineralizations of the external crystalline massifs could result from previous mineralizations remobilized during the Alpine orogeny. In this study, we propose to evaluate the metallogenic models of two Pb-Ag deposits of the French Alps: the Macôt-la Plagne (MP) and Peisey-Nancroix (PN). These two ore deposits are emblematic as being among the most economically important Pb-Ag ore deposits in France during the 19th century. They are located in the same litho-structural context in Permo-Triassic phyllitic quartzite of the External Briançonnais domain of the Internal Alps. In the absence of ages or mineralizing conditions (P, T, speciation), metallogenic models are still debated because of the stratiform or vein-like type mineralizations that could either result from sedimentary or orogenic processes [START_REF] Rogel | Le Gisement de Plomb de La Plagne (Savoie)-Etude Géologique et Métallogénique-Vanoise[END_REF][START_REF] Rogel | Exemple de Gisement de Galène du Trias des Alpes[END_REF][START_REF] Bouladon | Contribution à une systématique des gisements de plomb et de zinc[END_REF][START_REF] Raguin | La mine de plomb argentifère de La Plagne (Savoie)[END_REF][START_REF] Waline | La mine de La Plagne[END_REF][START_REF] Cabrol | Etude Géologique et Minière du Massif des Hurtières-Chaine de Belledonne (Savoie)-Alpes Françaises[END_REF][START_REF] Routhier | Géologie et Principes de Recherche[END_REF], as for other worldwide Pb-Zn(-Ag-Cu) deposits in several mountain belts (e.g., Pyrenees, Mount Isa region, and North China craton) mainly considered initially as sedimentary but revealed in fact to be orogen-driven [START_REF] Zhong | Ore-Forming Conditions and Genesis of the Huogeqi Cu-Pb-Zn-Fe Deposit in the Northern Margin of the North China Craton: Evidence from Ore Petrologic Characteristics[END_REF][START_REF] Gibson | Basin Inversion and Supercontinent Assembly as Drivers of Sediment-Hosted Pb-Zn Mineralization in the Mount Isa Region, Northern Australia[END_REF][START_REF] Cugerone | Structural Control on the Formation of Pb-Zn Deposits: An Example from the Pyrenean Axial Zone[END_REF][START_REF] Cugerone | Variscan U-Th-Pb Age for Stratabound Pb-Zn Mineralization in the Bossòst Dome (Pyrenean Axial Zone)[END_REF][START_REF] Cave | Linking Uplift and Mineralisation at the Mount Novit Zn-Pb-Ag Deposit, Northern Australia: Evidence from Geology, U-Pb Geochronology and Sphalerite Geochemistry[END_REF]. This study proposes for the first time a multitechnique investigation on these two deposits, coupling field observations, microstructure and mineral paragenesis identification, and thermobarometry and geochronology. Based on recent success in determining the ages of Alpine clefts by the in situ dating of monazite (e.g., [START_REF] Gasquet | Miocene to Messinian deformation and hydrothermal activity in a pre-Alpine basement massif of the French western Alps: New U-Th-Pb and argon ages from the Lauzière Massif[END_REF][START_REF] Grand'homme | Interpretation of U-Th-Pb in-situ ages of hydrothermal monazite-(Ce) and xenotime-(Y): Evidence from a large-scale regional study in clefts from the western alps[END_REF][START_REF] Gnos | Episodes of fissure formation in the Alps: Connecting quartz fluid inclusion, fissure monazite age, and fissure orientation data[END_REF]), a special interest was set on monazite characterization in all the investigated samples. Thermobarometry was based on detailed chlorite-phengite characterization using compositional maps.

Geological Setting and Ore Mineralization 2.1. Regional Geology

The Peisey-Nancroix and Macôt-la Plagne mining sites are located in the Western Alpine arc, in the Briançonnais domain (Internal Alps, Figure 1a). The arrangement of the nappe pile shows a subduction-type structure in which most of the tectonic units dip southeastward. The current architecture of the Internal Alps results from the closure of the Valais and Piedmont-Liguria oceans (branches of the Alpine Tethys) during the Late Cretaceous and Cenozoic convergence between the European and African (Apulian/Adrian) tectonic plates. The Briançonnais domain corresponds to the continental margin located between both oceans, which was subducted under the African tectonic plates and thrust over the European tectonic plate during the alpine convergence (Figure 1a) [START_REF] Schmid | Tectonic map and overall architecture of the Alpine Orogen[END_REF]. The Internal metamorphic Alps are separated from the "unmetamorphic" External Alps by a major thrust called the Penninic Front, which delimits the onset of the collision of the European margin at around 35 Ma (Figure 1a,b) [START_REF] Ceriani | Multi-stage thrusting at the "Penninic Front" in the Western Alps between Mont Blanc and Pelvoux massifs[END_REF]. The External Alps correspond to the European margin; they comprise the Dauphinois and the Ultra Dauphinois domains. The Internal Alps consist of a stack of oceanic (Valais and Piedmont Schistes Lustrés nappes and associated ophiolites) and continental-margin-derived nappes (the Briançonnais domain and some pieces of the Adrian plate such as the Margna-Sesia fragments). The Briançonnais domain is divided into two units (the External and Internal Briançonnais) separated by the Internal Briançonnais Front (IBF). The IBF is a major detachment zone highlighted by a Late Triassic (Carnian) formation called the "Nappe des Gypses", which is composed of hundreds of meters (from 100 to 400 m) of dominant massive evaporites (mainly anhydrites), marine carbonates (dolomites tectonically transformed into cargneules), and shales and clays (metamorphosed into the Schistes Lustrés micaschists) (Barré et al. [START_REF] Barré | Tectono-metamorphic evolution of an evaporitic décollement as recorded by mineral and fluid geochemistry: The "Nappe Des Gypses" (Western Alps) case study[END_REF] and references therein). To the east of the IBF, the Internal Briançonnais is composed of Mesozoic sediments lying on top of a pre-Carboniferous metamorphic basement locally called the Vanoise-Mont Pourri unit. To the west of the IBF, the External Briançonnais is composed of the Variscan Sapey gneiss [START_REF] Giorgis | The Mont-Mort metapelites:Variscan metamorphism and geodynamic context (Briançonnais basement, Western Alps, Switzerland)[END_REF][START_REF] Guillot | Zircon U-Pb geochronology of Ordovician magmatism in the polycyclic Ruitor Massif (Internal W Alps)[END_REF] and Carboniferous black schists containing anthracitic lenses, arkoses, and conglomerates, which are locally called the Briançonnais Houiller Zone (BHZ) [START_REF] Fabre | Contribution à l'étude de la Zone Houillère en Maurienne et en Tarentaise (Alpes de Savoie)[END_REF][START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF]. This basement rock is overlain by a Permo-Triassic sedimentary cover represented by a series of schists, quartzites, limestones, and dolomites from bottom to top [START_REF] Ellenberger | Etude géologique du pays de Vanoise[END_REF]. The External Briançonnais exhibits a low-grade alpine metamorphism, from lower greenschist facies to upper greenschist facies [START_REF] Goffé | Explanatory notes to the map: Metamorphic structure of the Alps Western and Ligurian Alps[END_REF][START_REF] Bousquet | Metamorphism of metasediments at the scale of an orogen: A key to the Tertiary geodynamic evolution of the Alps[END_REF].

Structural Evolution of the Briançonnais Zone

The internal part of the Western Alpine arc experienced a complex and polyphase tectonic history. Precise dating highlights a diachronic tectono-metamorphic evolution, with a rejuvenation of the ages towards the west. In the External Briançonnais, Bertrand et al. [START_REF] Bertrand | The Pennine Front zone in Savoie (Western Alps), a review and new interpretations from the Zone Houillère Briançonnaise[END_REF] proposed three compressional stages (D1, D2 and D3 events), followed by a late extensional one.

1.

The D1 event is characterized by an intense and penetrative schistosity (S1) oriented subparallel to S0 and by isoclinal folding. This early tectonic phase is interpreted as a top-to-the-northwest thrusting of the Internal Briançonnais over the External Briançonnais [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF] (Lanari et al. [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF] and references therein). The D1 deformation in the External Briançonnais occurred during the metamorphic peak. It is attributed to the exhumation of the Internal Briançonnais. There is no absolute age for the D1 deformation and metamorphic peak in the External Briançonnais. In the Internal Briançonnais, the D1 ranges from 55 and 43 Ma in the most internal parts [START_REF] Ganne | Les Dômes de Socles HP-BT dans le Domaine Pennique des Alpes Nord-Occidentales (Massifs d'Ambin et de Vanoise Sud): Modalités de Leur Exhumation[END_REF][START_REF] Gerber | Evolution Tectono-Métamorphique du Briançonnais Interne (Alpes Occidentales, Massifs de Vanoise Sud et d'Ambin): Comportement du Socle et de sa Couverture Dans un Contexte de Subduction Continentale Profonde[END_REF] and until 35 Ma in the most external part (in the Modane-Aussois unit) [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF]. The P-T conditions for D1 were constrained between 250-300 • C and 4-8 kbar in the BHZ [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF], and 320-370 • C and 9-11 kbar in the cover of the Internal Briançonnais [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF]. In the "Nappe des Gypses", the P-T conditions for the D1-D2 transition were estimated at 350 ± 20 • C and 6.5 ± 1.8 kbar using a δ 18 O thermometer on carbonates [START_REF] Barré | Tectono-metamorphic evolution of an evaporitic décollement as recorded by mineral and fluid geochemistry: The "Nappe Des Gypses" (Western Alps) case study[END_REF].

2.

The D2 event consists of a top-to-the-east shearing and thrusting highlighted by large-scale east-verging recumbent folds associated with the S2 schistosity, dipping to the west except into the eastern part of the BHZ. The D2 deformation occurred under greenschist facies conditions [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (Western Alps): Toward a new alternative exhumation model for the Briançonnais domain[END_REF]. In the most internal part of the Briançonnais domain, the large-scale D2 shear zones are often parallel to S1 and interpreted as detachment faults related to extension [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (Western Alps): Toward a new alternative exhumation model for the Briançonnais domain[END_REF]. The tectonic significance of the D2 event is still debated, as backthrusting [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF][START_REF] Platt | Structural history of high-pressure metamorphic rocks in the southern Vanoise massif, French Alps, and their relation to Alpine tectonic events[END_REF][START_REF] Freeman | Dating deformation using Rb-Sr in white mica: Greenschist facies deformation ages from the Entrelor shear zone, Italian Alps[END_REF][START_REF] Gabalda | Processus d'exhumation dans les Alpes Occidentales: Modélisation Géométrique et Reconstitution Géodynamique sur la Transversale Chartreuse-Maurienne, une Approche Multi-Échelle[END_REF] or as backfolding [START_REF] Bucher | Late-stage deformation in a collisional orogen (Western Alps): Nappe refolding, back-thrusting or normal faulting?[END_REF] during the continental collision, or as a detachment predating the continental collision [START_REF] Gerber | Evolution Tectono-Métamorphique du Briançonnais Interne (Alpes Occidentales, Massifs de Vanoise Sud et d'Ambin): Comportement du Socle et de sa Couverture Dans un Contexte de Subduction Continentale Profonde[END_REF][START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (Western Alps): Toward a new alternative exhumation model for the Briançonnais domain[END_REF]. Ages of backthrusting using 40 Ar/ 39 Ar and Rb-Sr dating on white mica are constrained between 35 and 32 Ma into the Internal Briançonnais [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF][START_REF] Freeman | Dating deformation using Rb-Sr in white mica: Greenschist facies deformation ages from the Entrelor shear zone, Italian Alps[END_REF]. The D2-D3 transition occurred at P-T conditions between 230-300 • C and 1-3 kbar in the "Nappe des Gypses" and Internal Alps sedimentary cover [START_REF] Barré | Tectono-metamorphic evolution of an evaporitic décollement as recorded by mineral and fluid geochemistry: The "Nappe Des Gypses" (Western Alps) case study[END_REF][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF][START_REF] Gerber | Evolution Tectono-Métamorphique du Briançonnais Interne (Alpes Occidentales, Massifs de Vanoise Sud et d'Ambin): Comportement du Socle et de sa Couverture Dans un Contexte de Subduction Continentale Profonde[END_REF].

3.

The D3 event is characterized by brittle structures and corresponds to the last stage of exhumation and to the formation of a large dome of crystalline basement associated with faulting in a brittle-ductile tectonic regime [START_REF] Gerber | Evolution Tectono-Métamorphique du Briançonnais Interne (Alpes Occidentales, Massifs de Vanoise Sud et d'Ambin): Comportement du Socle et de sa Couverture Dans un Contexte de Subduction Continentale Profonde[END_REF][START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (Western Alps): Toward a new alternative exhumation model for the Briançonnais domain[END_REF][START_REF] Lanari | P-T-t estimation of deformation in low-grade quartz-feldspar-bearing rocks using thermodynamic modelling and 40Ar/39Ar dating techniques: Example of the Plan-de-Phasy shear zone unit (Briançonnais Zone, Western Alps)[END_REF]. Bertrand et al. [START_REF] Bertrand | The Pennine Front zone in Savoie (Western Alps), a review and new interpretations from the Zone Houillère Briançonnaise[END_REF] suggest that the fan geometry of the BHZ results from the large-scale bending of the S2 schistosity due to D3 km-size folds. These D3 folds are associated with a sub-horizontal S3 schistosity (Lanari et al. [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF] and references therein). This last D3 event is related to the last stage of exhumation below 300 • C during Oligocene-Miocene times with a southward tilting around an east-west axis into the Briançonnais zone [START_REF] Barré | Tectono-metamorphic evolution of an evaporitic décollement as recorded by mineral and fluid geochemistry: The "Nappe Des Gypses" (Western Alps) case study[END_REF][START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF].

The Pb-Ag Mining District of Macôt-la Plagne (MP) and Peisey-Nancroix (PN)

The investigated Pb-Ag deposits of Macôt-la Plagne (MP) and Peisey-Nancroix (PN) are located in the Vanoise massif, at the contact of the Permo-Triassic External Briançonnais with cargneules (Figures 1 and2a). Cargneules are Triassic dolomites that have been brecciated and partly dissolved during the alpine orogeny from sulfate-rich fluids derived from Triassic gypsum. In the studied area, cargneules can be related to the "Nappe des Gypses", and therefore to the Internal Briançonnais Fault (IBF; Figure 1b). Both deposits are hosted in Triassic (Werfenian stage) quartzite (tQ on Figure 2a) and Permian sericite schists and phyllitic quartzite (rt on Figure 2a). The MP ore deposit is located within a recumbent dissymmetrical syncline, with a NNE/SSW axis (Figure 1c) [START_REF] Guillot | La Vallée du Ponturin (Vanoise Septentrionale, Alpes Françaises)[END_REF]. The PN deposit is bordered by two faults: the Mine Fault (M.F.) and the IBF (Figure 2a). The IBF is described as the thrusting of the Internal Briançonnais (Vanoise massif in the study zone) on top of the External Briançonnais, while the Mine Fault is a detachment level of the allochthonous External Briançonnais cover on the basement. In the study area, the IBF is characterized by two successive faulting activities, first with a top-to-the-west horizontal thrusting, followed by vertical normal faulting with a likely sinistral shearing component [START_REF] Guillot | La Vallée du Ponturin (Vanoise Septentrionale, Alpes Françaises)[END_REF].

The Pb-Ag mines of MP (main exploitation: 1810-1973) and PN (exploitation: 1734-1866) were the most economically important of the French Alps. The historical production is estimated at 125,600t of Pb and 360t of Ag (with an average content of 3.7% of Pb and 0.015% of Ag) and 22,000 t of Pb and 53 t of Ag, for MP and PN, respectively, with a peak of activity during the 19th century [START_REF] Meloux | Bilan Synthétique de dix Années de Prospection[END_REF][START_REF] Rogel | Le Gisement de Plomb de La Plagne (Savoie)-Etude Géologique et Métallogénique-Vanoise[END_REF][START_REF]Les Gisements de Pb-Zn Français[END_REF][START_REF] Gimard | Mine de Plomb Argentifère (Savoie)[END_REF]. In both deposits, the ore is characterized by Ag-rich galena, either massive or disseminated in a quartz-or barite-rich gangue. Tetrahedrite and pyrite are systematically associated with galena, with variable chalcopyrite, arsenopyrite, sphalerite, and bournonite [START_REF] Raguin | La mine de plomb argentifère de La Plagne (Savoie)[END_REF][START_REF] Routhier | Géologie et Principes de Recherche[END_REF][START_REF] Tachat | Mine de Plomb Argentifère[END_REF][START_REF] Gimard | Mines de Plomb Argentifère (Savoie)[END_REF]. Fluorite was never reported for these deposits.

Waline [START_REF] Waline | La mine de La Plagne[END_REF] describes three types of ores for the Macôt-la Plagne deposit (Figure 1c): (1) the Charles-Albert-type Pb-Ag veins, consisting of mineralized quartzite intercalated within the Permian schists, forming discontinuous lenses; (2) the Sarrazins-type veins, with banded fine-grained galena and barite; (3) pods-type ore, with well-crystallized galena hosted in mylonitic quartzite. The mineralization is well-aligned along the schistosity, with no particular relationship between the mineralization and the Permian schist-Triassic quartzite contact. The Pb-Ag mines of MP (main exploitation: 1810-1973) and PN (exploitation: 1734-1866) were the most economically important of the French Alps. The historical production is estimated at 125 600t of Pb and 360t of Ag (with an average content of 3.7% of Pb and 0.015% of Ag) and 22 000t of Pb and 53t of Ag, for MP and PN, respectively, with a peak of activity during the 19th century [START_REF] Meloux | Bilan Synthétique de dix Années de Prospection[END_REF][START_REF] Rogel | Le Gisement de Plomb de La Plagne (Savoie)-Etude Géologique et Métallogénique-Vanoise[END_REF][START_REF]Les Gisements de Pb-Zn Français[END_REF][START_REF] Gimard | Mine de Plomb Argentifère (Savoie)[END_REF]. In both deposits, the ore is characterized by Ag-rich galena, either massive or disseminated in a quartz-or barite-rich gangue. Tetrahedrite and pyrite are systematically associated with galena, with variable chalcopyrite, arsenopyrite, sphalerite, and bournonite [START_REF] Raguin | La mine de plomb argentifère de La Plagne (Savoie)[END_REF][START_REF] Routhier | Géologie et Principes de Recherche[END_REF][START_REF] Tachat | Mine de Plomb Argentifère[END_REF][START_REF] Gimard | Mines de Plomb Argentifère (Savoie)[END_REF]. Fluorite was never reported for these deposits.

Waline [START_REF] Waline | La mine de La Plagne[END_REF] describes three types of ores for the Macôt-la Plagne deposit (Figure 1c): (1) the Charles-Albert-type Pb-Ag veins, consisting of mineralized quartzite intercalated within the Permian schists, forming discontinuous lenses; (2) the Sarrazins-type veins, with banded fine-grained galena and barite; (3) pods-type ore, with well-crystallized galena hosted in mylonitic quartzite. The mineralization is well-aligned along the schistosity, with no particular relationship between the mineralization and the Permian schist-Triassic quartzite contact.

Field Description, Structural Analyses, and Sampling

Field Description and Structural Analyses

For the PN deposit, the main expressed foliation in the Permian and Triassic rocks (S1) is subhorizontal and dips to the east at the contact with the Mine Fault (M.F., Figure 2a,b). The coarse lineation is roughly NW-SE oriented. In the Permian mylonitic phyllitic quartzite (PN20-03, Figure 2c), the S1 foliation is highlighted by fine-grained white micas. Porphyroclasts of deformed quartz indicate a top-to-the-west kinematic. The C'-type shear bands are short, anastomosing, and wavy, and cross-cut the S1 mylonitic schistosity at relatively high angles (Passchier and Trouw [48] and references therein). The direction of deformation observed in the phyllitic quartzite is consistent with the C'-type shear bands structure dipping westward, shown by the Permian sericite schist (PN20-04, Figure 2d). S2 foliation was only observed to affect the pervasive S1 foliation as isoclinal folds on the other side of the Ponturin River.

Samples

The entrances of the two abandoned mines are sealed and no longer accessible for in situ structural observations or sampling. Four samples from the MP deposit were selected from Raguin's historical collection, which is kept by BRGM (Orléans, France). The samples come from two mineralized districts within Triassic quartzite: the Quartzite Est and the Sarrazins Districts. In the former district, a massive sulfide was selected from vein 0/8 (3501-9) and a disseminated ore from vein 0/10 (3501-4). Similarly, a massive sulfide (3501-10) and a more disseminated ore showing millimetric veins (3501-5) were selected from the Sarrazins District (Figure 1c). For the PN deposit, three samples were collected on a mine waste heap (Figure 1d). The mineralization occurs in veins or impregnation within mylonitic phyllitic Triassic quartzite (PN16-R05a, PN20-11, PN20-12).

Analytical Methods

Petrographical Observations

A detailed petrographic study was carried out at ISTerre using an optical microscope combined with back-scattered-electron imaging using a scanning electron microscope (SEM) Tescan Vega 3 in order to identify the mineralogical parageneses, including the sulfides-sulfosalts, the gangue, and accessory minerals. Particular attention was paid to micrometric minerals inclusions and possible compositional zoning in the back-scatteredelectron images.

Sulfides-Sulfosalts EPMA and LA-ICP-MS Trace Elements Protocols

Major and minor element analyses were performed using a JEOL JXA-8230 electron microprobe (EPMA) at ISTerre (Grenoble, France) equipped with five wavelengths dispersive spectrometers. Quantitative point analyses were made at 20 keV and 50 nA. The following 15 elements, standards, analytical lines, and monochromators were used for measurements: S (PbS, Ka, PETH), Cu (CuFeS2, Ka, LiFH), Fe (CuFeS2, Ka, LiFH), Pb (PbS, Ma, PETH), Zn (ZnS, Ka, LiFH), Sb (Sb, La, PETH), Ag (Ag, La, PETH), Cd (CdS, Lb, PETH), Sn (SnO2, La, PETH), Mn (MnSiO3, Ka, LiFH), Hg (HgS, La, LiFL), Co (CoO, Ka, LiFL), Ni (Ni, Ka, LiFL), As (FeAsS/InAs, La, TAP for galena and Ka, LiFL for chalcopyrite/tetrahedrite-tennantite), and Bi (Bi, La, LiFL for galena and Ma, PETH for chalcopyrite/tetrahedrite-tennantite). Peak counting times were scaled from 40 to 400 s for a total point analysis time of 400 s. Internal JEOL procedures were used for data reduction (i.e., quantification of element contents, uncertainties, and limits of detection). Measurements were collected on 3501-5a, 3501-9b, and PN16-R05a thin sections (30 µm). Among the 15 elements monitored, Co, Mn, Ni, and Bi were not measured, as well as As for galena, and Cd and Hg for chalcopyrite. Nonetheless, trace element analyses using EMPA are limited due to the compromise between the counting time, the number of measured elements, and the required limits of detection. Conditions for X-ray maps in the tetrahedrite grains from the PN16-R05a sample were fixed at 20 keV, 100 nA, dwell time of 200 ms, for a surface area of 300 × 360 µm with a step of 1 µm.

Trace element analyses in sulfides were carried out using Laser Ablation Inductively Coupled Plasma Mass Spectrometry (LA-ICP-MS) at the GeoRessources laboratory (Vandoeuvre-lès-Nancy, France), equipped with an ESI New Wave Excimer laser (193 nm) coupled with the Agilent 7900 quadrupole ICP-MS. Laser ablation conditions were fixed at 5Hz pulse rate, 6 mJ/cm 2 laser energy, 200 laser pulses. Selected zones of the minerals were ablated with 60 µm spot diameter from the 3501-9b, 3501-10a, PN20-11, and PN20-12 150 µm-thick sections. The ablated material was transported using a constant He flows (700 mL/min) and mixed with Ar (0.85l/min) in a cyclone coaxial mixer prior to entering the ICP torch and being ionized. Data were collected during a 30 ms dwell time per channel to enhance the number of counts for trace elements. The following isotopes were monitored: 34 S, 52 Cr, 55 Mn, 57 Fe, 59 Co, 60 Ni, 63 Cu, 66 Zn, 71 Ga, 74 Ge, 75 As, 82 Se, 95 Mo, 107 Ag, 111 Cd, 115 In, 118 Sn, 121 Sb, 125 Te, 197 Au, 202 Hg, 205 Tl, 208 Pb, and 209 Bi. All elements were measured with an integration time of 20 ms, except 34 S, 57 Fe (10 ms), 197 Au, and 205 Tl (30 ms). The 74 Ge was chosen over 72 Ge following the study of interferences from Bellissont et al. [START_REF] Belissont | LA-ICP-MS analyses of minor and trace elements and bulk Ge isotopes in zoned Ge-rich sphalerites from the Noailhac-Saint-Salvy deposit (France): Insights into incorporation mechanisms and ore deposition processes[END_REF]. External standard calibration was performed with the synthetic polymetallic sulfide material MASS-1 [START_REF] Wilson | Development of sulfide calibration standards for the laser ablation inductively-coupled plasma mass spectrometry technique[END_REF]. Data reduction was carried out using Iolite software (version 4) [START_REF] Paton | Iolite: Freeware for the visualisation and processing of mass spectrometric data[END_REF], following the standard methods of [START_REF] Longerich | Laser Ablation Inductively Coupled Plasma Mass Spectrometric Transient Signal Data Acquisition and Analyte Concentration Calculation[END_REF], and using Pb or Cu content-known from prior EPMA analyses (in same grains or averaged grains from the same location)-as an internal standard for galena, chalcopyrite, and tetrahedrite-tennantite solid solution, respectively.

Chlorite-Phengite EPMA Protocol and Thermometry Calculation

All investigated samples contain white mica (phengite). From the 7 rock samples investigated (Supplementary Table S1), the only MP sample containing both chlorite and phengite was selected for geothermobarometry (3501-5b). No chlorite was found in samples from the PN deposit. Comparison of the chemistry of phengite from the matrix with phengite in inclusions was performed on sample PN16-R05a, which was the most precisely investigated and the only one showing a poikilitic texture. Chemical compositions of chlorite and phengite were determined by point analyses and X-ray maps with electron microprobe analyses using a JEOL JXA-8230 at ISTerre. Quantitative point analyses were made at 15 keV and 10 nA. Conditions for X-ray maps in the selected sample 3501-5b were fixed at 15 keV, 100 nA, dwell time of 200 ms, for a surface area of 310 × 300 µm with a step of 0.5 µm. Raw X-ray maps were converted into oxide wt% maps using point analyses measured on the map area as internal standards [START_REF] De Andrade | Quantification of electron microprobe compositional maps of rock thin sections: An optimized method and examples[END_REF] with the software XMapTools 3.4.1 [START_REF] Lanari | XMapTools: A MATLAB©-based program for electron microprobe X-ray image processing and geothermobarometry[END_REF][START_REF] Lanari | Quantitative compositional mapping of mineral phases by electron probe micro-analyser[END_REF]. Structural formulae were then calculated for each pixel on a 14-oxygen basis and 11-oxygen basis for chlorite and phengite, respectively.

Chlorite analyses that did not satisfy the following criteria were rejected: (i) Alkali sum (K 2 O + Na 2 O + CaO) < 1 wt% [START_REF] Bourdelle | Low-temperature chlorite geothermometry: A graphical representation based on a T-R2+ -Si diagram[END_REF][START_REF] Vidal | Deciphering temperature, pressure and oxygen-activity conditions of chlorite formation[END_REF], and (ii) Octahedral sum (Al VI + Fe + Mn + Ca + Na + K) ≤ 6 apfu. Phengite analyses that did not satisfy the following criteria were rejected (following Parra et al. [START_REF] Parra | A thermodynamic model for Fe-Mg dioctahedral K white micas using data from phase-equilibrium experiments and natural pelitic assemblages[END_REF]): (i) 92 wt% ≤ Oxygen sum ≤ 96 wt%; (ii) TiO 2 + MnO + CaO < 0.5 wt% and (iii) Alkali sum ≥ 0.9 apfu. The K-mica-quartz-water thermodynamic model of Dubacq et al. [START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: Thermodynamic modelling and implications for thermobarometric estimates[END_REF] allows modeling K-mica compositional variation as a function of temperature and pressure.

The temperature range of chlorite growth was calculated for each pixel of the chlorite compositional map using the chlorite-quartz-water thermometer of [START_REF] Vidal | A thermodynamic model for Fe-Mg aluminous chlorite using data from phase equilibrium experiments and natural pelitic assemblages in the 100 to 600 • C, 1 to 25 kb range[END_REF][START_REF] Vidal | Thermodynamic properties of the Tschermak solid solution in Fe-Chlorite: Application to natural examples and possible role of oxidation[END_REF][START_REF] Vidal | P-T-Deformation-Fe3+/Fe2+ mapping at the thin section scale and comparison with XANES mapping: Application to a garnet-bearing metapelite from the Sambagawa metamorphic belt (Japan)[END_REF] with individual estimation of the Fe 3+ /Fe tot (XFe 3+ ). Pressure was fixed at 5 kbar, water activity at 1, and the minimum XFe 3+ at 10% to avoid the overestimation of temperature highlighted by Lanari [START_REF] Lanari | Micro-Cartographie P-T-ε Dans les roches Métamorphiques[END_REF], especially for HT (>350 • C)-LP (<5 kbar) chlorite. An average temperature was calculated from chlorite zones using a Monte Carlo technique (3000 iterations for each pixel) to test the temperature distribution. Conventional empiric thermometer based on the Al IV content of chlorite [START_REF] Cathelineau | Cation site occupancy in chlorites and illites as a function of emperature[END_REF] was also applied for comparison.

The temperature range of phengite growth was estimated using the phengite-quartzwater thermometer of Dubacq et al. [START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: Thermodynamic modelling and implications for thermobarometric estimates[END_REF]. The equilibrium conditions of K-mica-quartz-water were represented by a P-T line along which the water content varies. The thermodynamic model of Dubacq et al. [START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: Thermodynamic modelling and implications for thermobarometric estimates[END_REF] was also used for pressure estimates at fixed temperature values (see Lanari et al. [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF]).

Monazite Composition Analysis and Dating (EPMA and LA-ICP-MS)

Monazite chemical compositions were analyzed using a JEOL JXA-8230 electron microprobe at ISTerre to unravel intragrain zoning on the PN16-R05a and 3501-5b thin sections. The X-ray maps were acquired at 15keV and 100nA, and 1000 ms dwell time for a surface area of 170 × 55 µm and using 1 µm step on a PN16-R05a monazite grain. Quantitative point analyses of the monazite were obtained using a 15keV and 200nA beam modified after Scherrer et al. [START_REF] Tachat | Mine de Plomb Argentifère[END_REF]. From U, Th, and Pb concentrations, a chemical age was then calculated using NiLeDAM software [START_REF] Seydoux-Guillaume | Low-temperature alteration of monazite: Fluid mediated coupled dissolution-precipitation, irradiation damage, and disturbance of the U-Pb and Th-Pb chronometers[END_REF][START_REF] Villa-Vialaneix | Monazite Datation for the NiLeDAM Team[END_REF] with an age uncertainty at 2σ. All analyses with SiO 2 content > 1 wt% (contamination from surrounding silicates) or Pb detection below the detection limit of around 100 ppm were not taken into account.

In situ monazite U-Pb isotopic dating was done on one PN (PN16-R05a) and three and MP thin sections (3501-4a, 3501-5b, and 3501-10a) using the LA-ICP-MS at the BRGM (Orléans, France). The mass spectrometer is an ICP-MS X series II set up in Xs lenses coupled with a Cetac Excite 193 nm laser at short timing pulse. The sample was placed in a double-volume cell swept by a helium flux (analysis conditions are summarized in Supplementary Table S2). Analyses were normalized to the Trebilcock monazite (272 ± 2 Ma) [START_REF] Tomascak | U-Pb Monazite geochronology of granitic rocks from Maine: Implications for late Paleozoic tectonics in the northern Appalachians[END_REF] with a quality control on the Madmon standard monazite (523 ± 12 Ma) [START_REF] Schulz | In situ Th-Pb dating of monazite by 266 nm laser ablation and ICP-MS with a single collector, and its control by EMP analysis[END_REF]. Reduction data was done using the Glitter © software. The 235 U contents were recalculated from measured contents of 206 Pb, 207 Pb, and 238 U using the disintegration equation. The age calculation was performed using IsoplotR on R [START_REF] Vermeesch | IsoplotR: A free and open toolbox for geochronology[END_REF] with a standard deviation of 2σ, and the isochron was drawn using an ordinary least squares regression. A size spot of 10 or 15 µm was chosen to maximize the signal. For the PN deposit, the analyses which statistically diverge from the Discordia line have been removed from the age calculation. A SEM control was then realized for the 6 removed analyses, providing evidence of interference as inclusions or mineral interfaces.

Petrographic Results

The MP and PN Pb-Ag deposits show significant similarities in terms of microstructure, ore mineralization, and mineral assemblages. Consequently, a unique petrogenetic sequence is proposed for the two deposits.

Sample Microstructures

The MP and PN Pb-Ag deposits correspond to mineralized quartzite and to phyllitic quartzite (1 to 15% white mica; Figure 3) with intergranular disseminated mineralization along the foliation, millimetric to pluricentimetric veins, and massive sulfides (Figure 3).

Samples with a disseminated texture have a mylonitic fabric consisting of quartz porphyroclasts surrounded mainly by galena and phengite (Figures 3 and4). Submillimetric to millimetric anhedral quartz porphyroclasts (Qz1) are elongated along the stretching lineation (Figures 3c and4a,c). They exhibit undulose extinction and are replaced and surrounded by recrystallized equigranular and isometric grains with a micrometer size (Qz2). Relics of Qz1 porphyroclasts show pervasive intragranular fractures at high angles to the foliation that can be healed with trails of fluid inclusions or sealed by Qz2 with some minor sulfide-sulfosalts (mainly galena; Figure 4a,c). Newly formed Qz2 and sulfide-sulfosalts assemblages (galena, pyrite, tetrahedrite-tennantite, and chalcopyrite) can also occur in dilatation zones that crosscut the healed microcracks (Figure 4b), or in strain shadows of Qz1 (Figure 4a,b). In some places, Qz1 is completely replaced by Qz2 associated with sulfides, but the Qz1 porphyroclasts (or Qz1 porphyroclast aggregates) initial shape is preserved and molded by phengite and sulfides (Figures 3f and4c,d). Apart from quartz, phengite is the second main silicate associated with the mineralization. Phengite occurs as micrometric acicular grains oriented within the S1 foliation (Figures 3c and4), but it also sometimes appears as randomly oriented and folded aggregates up to 100 µm in size.

Paragenetic Sequence

The ore is mainly made of galena with variable amounts of pyrite, tetrahedritetennantite, and chalcopyrite, and rare arsenopyrite, cobaltite-gersdorffite, sphalerite, and bournonite (Figures 345). In samples with a massive texture, galena contains relics of Qz1 porphyroclasts recrystallized into Qz2, pyrite, and rare gypsum (Figure 3b,f) associated with fractured dolomite (Figure 3b). In disseminated texture, sulfides-sulfosalts occur mostly within the foliation between quartz and phengite (Figure 4c,d) or in millimetric to centimetric sulfide veins parallel or perpendicular to the foliation often intergrown with tetrahedrite-tennantite, pyrite, or chalcopyrite (Figures 3a,c,e and4e,f). In these samples, galena may occur as fine-grained intergrowth (displaying 120 • triple-junction grain boundaries; Figure 5f) with tetrahedrite-tennantite, chalcopyrite, pyrite, barite, and dolomite-ankerite (Dol-Ank), displaying sector zoning with variable Mg, Fe, and Mn contents (Figures 3e and8a), quartz, monazite, and phengite (Figures 3b, 4e,f and5). Apart from galena, pyrite can be abundant within some veins. Galena, monazite, barite, and dolomite can be found in inclusions within pyrite (Figure 5e,f). Generally, PN samples contain more abundant tetrahedrite-tennantite compared to MP, where this mineral is rare. In contrast, the MP deposit contains abundant barite, which is rarely represented in PN. One sample from PN shows peculiar features, such as (i) a poikilitic texture with abundant inclusions of quartz and phengite found in sulfides and sulfosalts, and (ii) the presence of microspheres of florencite (a REE-aluminophosphate) in sulfides-sulfosalts, Qz2 micrograins, and phengite (Figure 5c). The entire description above is resumed in the paragenesis presented in Figure 6.

Sulfides and Sulfosalts Trace Element Composition

Trace elements were investigated first using EPMA and then LA-ICP-MS. Sulfidessulfosalts trace element compositions are summarized in Table 1. For galena, the trace elements are generally below the EPMA detection limits, except for Ag and Sb. With LA-ICP-MS, variations in Ag, Sb, Bi, Tl, Cd, As, Sn, Mn, Fe, Hg, and Cr contents are detected, while Te, Au, Mo, Ge, and Ga contents remain below the detection limits. The MP galena (3501-5a, 3501-9b, and 3501-10a) shows a significant correlation between Ag (480-1665 ppm) and Sb (580-2080 ppm) concentrations (Figure 7), which attests for a coupled substitution of these elements (George et al. [START_REF] George | Trace and Minor Elements in Galena: A Reconnaissance LA-ICP-MS Study[END_REF] and references therein). The highest Ag and Sb contents are only observed within the massive galena texture of 3501-9b. In 3501-10a (Figure 3b), massive galena within the dolomite fractures (zone 1) exhibit slightly lower to similar Ag and Sb contents than the massive galena (3501-9b. Opposite, the fine-grained galena (zone 2) has much lower Ag-Sb contents. Galena from the PN samples (PN16-R05a, PN20-11 and PN20-12) displays low Ag (<540 ppm) and Sb (<700 ppm) contents, similar

Sulfides and Sulfosalts Trace Element Composition

Trace elements were investigated first using EPMA and then LA-ICP-MS. Sulfidessulfosalts trace element compositions are summarized in Table 1. For galena, the trace elements are generally below the EPMA detection limits, except for Ag and Sb. With LA-ICP-MS, variations in Ag, Sb, Bi, Tl, Cd, As, Sn, Mn, Fe, Hg, and Cr contents are detected, while Te, Au, Mo, Ge, and Ga contents remain below the detection limits. The MP galena (3501-5a, 3501-9b, and 3501-10a) shows a significant correlation between Ag (480-1665 ppm) and Sb (580-2080 ppm) concentrations (Figure 7), which attests for a coupled substitution of these elements (George et al. [START_REF] George | Trace and Minor Elements in Galena: A Reconnaissance LA-ICP-MS Study[END_REF] and references therein). The highest Ag and Sb contents are only observed within the massive galena texture of 3501-9b. In 3501-10a (Figure 3b), massive galena within the dolomite fractures (zone 1) exhibit slightly lower to similar Ag and Sb contents than the massive galena (3501-9b. Opposite, the fine-grained galena (zone 2) has much lower Ag-Sb contents. Galena from the PN samples (PN16-R05a, PN20-11 and PN20-12) displays low Ag (<540 ppm) and Sb (<700 ppm) contents, similar to the lower Ag-Sb contents from the disseminated MP samples (3501-5b and 3501-10a zone2; Figure 7). The Bi content in the PN galena (up to 60 ppm) is significantly higher than in the MP galena (<1 ppm), especially in vein-type galena. Both the MP and PN galena analyses are aligned along the main coupled substitution: Ag + + (Sb + Bi) 3+ → 2Pb 2+ . The microprobe analyses display the same correlation between the Ag-Sb-rich galena (3501-9b) and Ag-Sb-poor galena from the vein-type mineralization (3501-5a and PN16-R05a; Figure 7). They also show a slightly richer Ag-Sb galena in larger veins (PN16-R05a; Figure 3a) compared to millimetric veins (PN20-11; Figure 3c). Variations of Mn, Fe, and Hg (up to 3100 ppm for Fe in one analysis) are not correlated with other trace elements. Cd and Cu show content variation mainly ranging from 1 to 100 ppm, with rare concentrations up to 400 ppm. 3b) and PN samples. Circles and diamonds represent LA-ICP-MS and EPMA analyses, respectively. The correlation line was drawn using all the dataset; the correlation coefficient is given by R².

Thermobarometric Results

Chlorites

Rosette-shape chlorite aggregates show an oscillatory zoning pattern in BSE images (Figure 8) and X-ray maps, mainly caused by Fe and Mg content variations (XMg = Mg/(Mg + Fetot)). Based on variations in Mg, Si, Al, and Mn content, three different composition zones can be recognized from the core to the rim of the rosette (Figure 8; Table 2). Three representative composition domains (Figure 8) were selected from the chlorite zoned map (n = 36,616). The chlorite zones have a size of 138, 86, and 66 pixels from Chl1 1 to Chl3, respectively. From core (Chl1) to the rim (Chl3), the Mg and Al increase, while Si, Fe, and Mn decrease with an Mg-poor, Fe-rich intermediate rim (Chl2) (Figure 8). The Si and Al anticorrelation do not follow exactly the same features as the XMg and Mn contents. The Si and Al contents have a moderate chemical transition occurring in the XMgpoor chlorite (Chl2). Textures and chemical zoning, from Chl1 to Chl3, display chlorite growth without dissolution patterns. Except for FeO and MgO contents, other chemical variations are under 1 wt% (Table 2). The Mn content displays an oscillatory zoning anticorrelated with the XMg zoning (Figure 8). It reaches its maximum with 0.065 apfu in the Chl2 domain. Assuming that the Fe measured is purely divalent, chlorite compositions are very close to the calculation limits of the semiempirical and thermodynamic models, with an average vacancy value of 0.03 ± 0.03 apfu (Masci, [START_REF] Masci | Cristallochimie du fer dans les Chlorites Métamorphiques: Approche Analytique Multiéchelle, Expérimentale, et Implications Pétrologiques[END_REF] and references therein). So, to avoid meaningless temperatures (e.g., negative or > 1000 °C values using Inoue et al. [START_REF] Inoue | Application of chemical geothermometry to low-temperature trioctahedral chlorites[END_REF] and Bourdelle et al. [START_REF] Bourdelle | A new chlorite geothermometer for diagenetic to low-grade metamorphic conditions[END_REF]), it was necessary to evaluate the Xfe 3+ = Fe 3+ /Fetot value to approach the real chemical chlorite compositions and, by extension, to estimate the chlorite temperature formation. The Xfe 3+ estimated with the chlorite-quartz-water multiequilibria approach of Vidal et al. [START_REF] Vidal | Thermodynamic properties of the Tschermak solid solution in Fe-Chlorite: Application to natural examples and possible role of oxidation[END_REF] is 32 ± 5% for Chl1, 30 ± 4% for Chl2, and 27 ± 5% for the Chl3. Then, the ranges of temperature calculated for each zone are 285 ± 30 °C, 280 ± 30 °C, and 315 ± 35 °C, for Chl1, Chl2, and Chl3, respectively (Figure 8f).

The empirical thermometer of Cathelineau [START_REF] Cathelineau | Cation site occupancy in chlorites and illites as a function of emperature[END_REF] based on the Al IV content shows a lower temperature range for Chl1 and Chl2 (320 ± 10 °C) than for Chl3 (330 ± 10 °C). Thus, the Al IV content is similar for Chl1 and Chl2 (1.19 ± 0.03 apfu) and slightly higher for Chl3 (1.23 ± 0.03 apfu). Chl1 and Chl2 display a similar Al IV content (1.19 ± 0.03 apfu), while Chl3 shows slightly higher Al IV content (1.23 ± 0.03 apfu). The range (280-330 °C) and the 3b) and PN samples. Circles and diamonds represent LA-ICP-MS and EPMA analyses, respectively. The correlation line was drawn using all the dataset; the correlation coefficient is given by R 2 .

The tetrahedrite group is a copper sulfosalt solid solution between tetrahedrite (Sb), tennantite (As), and freibergite (Ag) end-members [START_REF] Moëlo | Sulfosalt systematics: A review. Report of the sulfosalt sub-committee of the IMA commission on ore mineralogy[END_REF][START_REF] Biagioni | The tetrahedrite group: Nomenclature and classification[END_REF]. The tetrahedrite group mineral is found as anhedral grains often in intergrowth with galena, chalcopyrite, phengite, and monazite (Figures 4 and5). In both deposits, the tetrahedrite group displays mainly between the tetrahedrite (Sb) and tennantite (As) end-member (Table 1), with mainly a tennantite core and tetrahedrite rim observed in PN (Figure S1). In PN16-R05a, both galena and tetrahedrite show a poikilitic texture with inclusions of microquartz (Qz2), phengite, and florencite opposite to the inclusion-free texture (Figure 5a,d). The tetrahedrite grains display significant intragranular variations for Sb (21-25 wt%), As (3.2-5.6 wt%), and Ag (5.3-6.7 wt%). These variations reflect a complex zoning with a good correlation between Sb and Ag and an anticorrelation between As and Sb-Ag. Chemical zoning in tetrahedrite is not readily correlated with poikilitic or inclusion-free textures. However, tetrahedrite shows intergrowth texture with a Pb-Ag-depleted galena, which is in agreement with the tetrahedrite-preferred host for Ag-Sb in case of cocrystallizing tetrahedrite-galena assemblage [START_REF] George | Minor and trace elements in natural tetrahedrite-tennantite: Effects on element partitioning among base metal sulphides[END_REF].

Pyrite is ubiquitous in both Pb-Ag deposits and mainly found as (i) inclusions within galena or (ii) as the main ore mineral in sample 3501-5, where it forms a vein-type structure. Pyrite is usually found as subeuhedral to euhedral cubic crystals that may contain inclusions of galena, barite, phengite, hematite, florencite, and monazite. The inclusions may locally highlight a cubic growth zoning pattern (Figure 5e,f).

Chalcopyrite was only, but systematically, identified in the PN samples, associated with galena in millimetric veins perpendicular to the foliation, and in inclusions or intergrowth with the galena and tetrahedrite-tennantite aligned within the foliation (Figures 3c and4e,f). All of the chalcopyrite trace elements contents are under 130 ppm, with no significant variation (Table 1). Other ore sulfides-sulfosalts, such as bournonite, cobaltite-gersdorffite, and sphalerite, were exceptionally observed in our samples.

Thermobarometric Results

Chlorites

Rosette-shape chlorite aggregates show an oscillatory zoning pattern in BSE images (Figure 8) and X-ray maps, mainly caused by Fe and Mg content variations (XMg = Mg/ (Mg + Fe tot )). Based on variations in Mg, Si, Al, and Mn content, three different composition zones can be recognized from the core to the rim of the rosette (Figure 8; Table 2). Three representative composition domains (Figure 8) were selected from the chlorite zoned map (n = 36,616). The chlorite zones have a size of 138, 86, and 66 pixels from Chl1 1 to Chl3, respectively. From core (Chl1) to the rim (Chl3), the Mg and Al increase, while Si, Fe, and Mn decrease with an Mg-poor, Fe-rich intermediate rim (Chl2) (Figure 8). The Si and Al anticorrelation do not follow exactly the same features as the XMg and Mn contents. The Si and Al contents have a moderate chemical transition occurring in the XMg-poor chlorite (Chl2). Textures and chemical zoning, from Chl1 to Chl3, display chlorite growth without dissolution patterns. Except for FeO and MgO contents, other chemical variations are under 1 wt% (Table 2). The Mn content displays an oscillatory zoning anticorrelated with the XMg zoning (Figure 8). It reaches its maximum with 0.065 apfu in the Chl2 domain. Assuming that the Fe measured is purely divalent, chlorite compositions are very close to the calculation limits of the semiempirical and thermodynamic models, with an average vacancy value of 0.03 ± 0.03 apfu (Masci,[START_REF] Masci | Cristallochimie du fer dans les Chlorites Métamorphiques: Approche Analytique Multiéchelle, Expérimentale, et Implications Pétrologiques[END_REF] and references therein). So, to avoid meaningless temperatures (e.g., negative or >1000 • C values using Inoue et al. [START_REF] Inoue | Application of chemical geothermometry to low-temperature trioctahedral chlorites[END_REF] and Bourdelle et al. it was necessary to evaluate the Xfe 3+ = Fe 3+ /Fe tot value to approach the real chemical chlorite compositions and, by extension, to estimate the chlorite temperature formation. The Xfe 3+ estimated with the chlorite-quartz-water multiequilibria approach of Vidal et al. [START_REF] Vidal | Thermodynamic properties of the Tschermak solid solution in Fe-Chlorite: Application to natural examples and possible role of oxidation[END_REF] is 32 ± 5% for Chl1, 30 ± 4% for Chl2, and 27 ± 5% for the Chl3. Then, the ranges of temperature calculated for each zone are 285 ± 30 • C, 280 ± 30 • C, and 315 ± 35 • C, for Chl1, Chl2, and Chl3, respectively (Figure 8f). The empirical thermometer of Cathelineau [START_REF] Cathelineau | Cation site occupancy in chlorites and illites as a function of emperature[END_REF] based on the Al IV content shows a lower temperature range for Chl1 and Chl2 (320 ± 10 • C) than for Chl3 (330 ± 10 • C). Thus, the Al IV content is similar for Chl1 and Chl2 (1.19 ± 0.03 apfu) and slightly higher for Chl3

(1.23 ± 0.03 apfu). Chl1 and Chl2 display a similar Al IV content (1.19 ± 0.03 apfu), while Chl3 shows slightly higher Al IV content (1.23 ± 0.03 apfu). The range (280-330 • C) and the tendency of increasing temperature between Chl1-Chl2 and Chl3 are similarly displayed by both empirical and thermodynamic chlorite thermometers.

Phengites

In the chlorite-bearing sample (3501-5b), euhedral acicular phengite is commonly found parallel and in contact with the individual platelet Chl3 crystals (Figure 8a) and is thus assumed to be cogenetic with Chl3, at least for its outer compositions. Variations in Si and Al contents allow two groups of phengite to be distinguished (Table 2; Figure 9a,b). The first group (Ph1: n = 212) has lower Si and higher Al contents (3.31-3.45 and 2.06-2.24 apfu, respectively) than the second group (Ph2: n = 120, Si of 3.37-3.52 apfu, Al of 2.00-2.18 apfu). The temperature of Ph1 growth estimated with the thermometer of Dubacq et al. [START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: Thermodynamic modelling and implications for thermobarometric estimates[END_REF] is slightly higher (275-360 • C at a fixed reference pressure of 5 kbar) than that for Ph2 (250-330 • C, Figure 9). However, the quantitative point analyses (used as internal standards for the Ph map) show intermediate temperature included within both Ph zones (individual Ph from MP: n = 18, 275-310 • C at 5 kbar, Figure 9). The Ph1 and Ph2 P-T dispersion may be explained by the Ph textures (acicular or larger crystals) and are not sustained by the individual Ph analyses from the same area. Thus, for more accuracy and representativeness, the individual Ph data were used to constrain the P-T deposit formation. Following the approach proposed by Lanari et al. [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF], approximative pressure conditions for the growth of phengite and chlorite aggregates are obtained by combining results of phengite-quartz-water (for Individual Ph from MP) and chlorite-quartz-water (for Chl3) models, respectively. They yield at 315 ± 35 • C and 6.25 ± 0.75 kbar (Figure 9).

The PN phengite shows no significant variations between microstructural sites, and overall mineralogical variations are moderate (n = 15, Table 2), with Si content ranging from 3.19 to 3.39 apfu. The PN phengite displays Al content ranging between 2.32 and 2.55 apfu and Fe content between 0.05 and 0.07 apfu. Despite the PN analyses of phengite being in the same SiO 2 content range as those from MP, it appears that PN phengite is closer to the muscovite + paragonite end-members (0.70 ± 0.04) in comparison to those from MP (0.54 ± 0.02). From both deposits, the maximum pyrophillitic ratio calculated in a single analysis is 0.10 (XPrl in Table 2). Hence, the Si content is controlled by the Tschermark substitution between celadonite and muscovite (respectively, XCel and XMs, as shown in Table 2) at the same bulk rock composition (Dubacq et al. [START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: Thermodynamic modelling and implications for thermobarometric estimates[END_REF] and references therein). However, at a fixed pressure (5 kbar), the PN phengite growth temperature using Dubacq et al. [START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: Thermodynamic modelling and implications for thermobarometric estimates[END_REF] ranges between 280 and 330 • C, which is in agreement with both the MP phengite and chlorite temperature formation (Figures 8 and9). The black boxes show the pixel zones analyzed. (c) Phengite P-T estimates (lines) using the method of Dubacq et al. [START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: Thermodynamic modelling and implications for thermobarometric estimates[END_REF]. The box is constructed using the result of Chl3 thermometry (Figure 8f) and the punctual data. For a given temperature, pressure has been estimated using the phengite P-T equilibrium lines of the punctual data from the map area. Boundaries of the Ph1 and Ph2 were determined using more than 95% of P-T lines for each zone. [START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: Thermodynamic modelling and implications for thermobarometric estimates[END_REF]. The box is constructed using the result of Chl3 thermometry (Figure 8f) and the punctual data. For a given temperature, pressure has been estimated using the phengite P-T equilibrium lines of the punctual data from the map area. Boundaries of the Ph1 and Ph2 were determined using more than 95% of P-T lines for each zone.

Monazite Characterization and U-Pb In Situ Dating

Despite the micrometric size (<10 µm) and similarity in contrast and brightness between the monazite and barite in the SEM-BSE mode (Figure 5e), few crystals of monazite were found in all rock samples, except for the massive sulfides. However, in only one sample from the PN deposit (PN16-R05a), monazite was abundant and large enough for detailed petrological characterization and further chemical and isotopic dating. For the MP deposit, monazite grains from three different samples were dated.

Monazite Characterization and Composition

In PN16-R05a, monazite crystals show an elongated subhedral shape with a length and width up to 100 µm and 20 µm. Monazite has a poikilitic texture with inclusions of quartz, phengite, and florencite. It occurs in the recrystallized Qz2 domain or inclusions in sulfide-sulfosalt. Elongated grains can share sharp straight boundary contact with sulfide-sulfosalt or phengite. In PN16-R05a, the monazite grains are zoned (Figure 10; Supplementary Table S3). The core corresponds to Nd-monazite (Nd > Ce wt% content) enriched in middle-REE, especially in Eu (up to 5 wt%), U, Th, Ca, Si, and Y. The rim corresponds to a Ce-monazite rim (Ce > Nd content, Th < 3 wt% and Ca < 0.8 wt%) enriched in La (light-REE). In the MP samples, monazite shows minor variations with a typical composition of monazite-(Ce) and a lower U and Th compared to those obtained in the core of the PN monazite-(Nd).

Monazite U-Th-Pb Chemical Dating

In PN16-R05a, U can reach up to 1.24 wt% in the core of monazite-(Nd). Uranium generally tends to decrease towards the rim. Similarly, Th is high but decreases to extremely low values towards the rim. In the Th-and U-rich core, it was possible to calculate a chemical U-Th-Pb age as young as 33 Ma (Figure 10), excluding the presence of an inherited core. The rim provides meaningless older ages and is to low contribution of radiogenic Pb in the U-and Th-poor rim (EPMA ages of Figure 10).

As regards MP monazite grains, the systematically low Pb content, below or close to the detection limit (Supplementary Table S3), could not allow calculating an age, but also suggests a young crystallization with no evidence of an older event.

Monazite U-Th-Pb Chemical Dating

In PN16-R05a, U can reach up to 1.24 wt% in the core of monazite-(Nd). Uranium generally tends to decrease towards the rim. Similarly, Th is high but decreases to extremely low values towards the rim. In the Th-and U-rich core, it was possible to calculate a chemical U-Th-Pb age as young as 33 Ma (Figure 10), excluding the presence of an inherited core. The rim provides meaningless older ages and is to low contribution of radiogenic Pb in the U-and Th-poor rim (EPMA ages of Figure 10).

As regards MP monazite grains, the systematically low Pb content, below or close to the detection limit (Supplementary Table S3), could not allow calculating an age, but also suggests a young crystallization with no evidence of an older event.

Monazite In Situ Isotopic U-Pb Dating

In total, 44 LA-ICP-MS analyses were performed on twenty-six monazite individual grains on PN16-R05a (PN ore), and seventeen analyses were performed on fourteen monazite grains from three MP samples (3501-4b, n = 3; 3501-5b, n = 11; 3501-10a, n = 3), as shown in Table 3.

Table 3. Summary of the LA-ICP-MS raw data for Peisey-Nancroix (PN16-R05a) and Macôt-la Plagne monazite (3501-4b, 3501-5b, and 3501-10a). Note that monazite grains are referenced by the first digit in the analysis n • after the sample names, and the spot numbers in each monazite grain are displayed by the third and last digit. In the Tera-Wasserburg diagram, the data align well on a single line for both deposits. The lower intercept of the Concordia line provides an estimated age (without correction) of 35.1 ± 0.6 Ma (MSWD = 2.1; n = 38) for PN (Figure 11a) and 34.2 ± 3.5 Ma (MSWD = 1.3; n = 17) for MP (Figure 11c). By fixing the value of the initial 207 Pb/ 206 Pb ratio at 0.835, using the bulk-Pb evolution model of Stacey and Kramers [START_REF] Stacey | Approximation of terrestrial lead isotope evolution by a two-stage model[END_REF] and Cumming and Richards [START_REF] Cumming | Ore lead isotope ratios in a continuously changing Earth[END_REF], it is possible to correct the calculated age from the common Pb contamination. Thus, the corrected age is 34.9 ± 0.5 Ma (MSWD = 2.1, n = 38/44) for PN monazite grains (Figure 11b) and 33.1 ± 2.9 Ma (MSWD = 1.3, n = 17) for MP (Figure 11d). Anchoring the Pb initial composition to the common Pb values (Stacey and Kramers model [START_REF] Stacey | Approximation of terrestrial lead isotope evolution by a two-stage model[END_REF]) has no significant impact on the monazite age, indicating that the isotopic Pb composition of the MP and PN mineralizing fluids are close to the predicted common Pb values, opposite to some mineralizing fluids in relation with the Penninic front [START_REF] Rossi | Geochemical Variations and Element Transfer during Shear-Zone Development and Related Episyenites at Middle Crust Depths: Insights from the Mont Blanc Granite (French-Italian Alps)[END_REF][START_REF] Cenki-Tok | Direct Dating of Mid-Crustal Shear Zones with Synkinematic Allanite: New in Situ U-Th-Pb Geochronological Approaches Applied to the Mont Blanc Massif[END_REF]. Moreover, even if monazite grains from three different rock samples are considered for the MP deposit, there is one sample (3501-5b) more represented and well dispersed along the Discordia line, which gives an age of 32.15 ± 3.25 Ma (n = 11) within the uncertainty of the age obtained by considering the entire MP dataset (Figure 11d). [START_REF] Stacey | Approximation of terrestrial lead isotope evolution by a two-stage model[END_REF] and Cumming and Richards [START_REF] Cumming | Ore lead isotope ratios in a continuously changing Earth[END_REF], to correct the calculated age from the common Pb contamination. Forty-four analyses from twenty-six monazite grains (white ellipses in (a,b)) were used for calculating the age of PN16-R05a (gray ellipses (n = 6) were removed from the age calculation after SEM control, §4.4.). Seventeen analyses from fourteen monazite grains were used for calculating the age (without removed data) from the Macôtla Plagne samples.

Discussion

The Pb-Ag mineralizations of the two MP and PN deposits are located in a similar litho-structural domain. Further, they show a main similar mineral paragenesis, deformation microstructure, and ages, suggesting that they formed under similar conditions. 3) for the Peisey-Nancroix (a,b) and the Macôt-la Plagne (c,d) deposits. (a,c) provides Discordia lines and ages calculated without any correction, whereas for (b,d) the initial 207 Pb/ 206 Pb ratio is fixed at 0.835, using the bulk-Pb evolution model of Stacey and Kramers [START_REF] Stacey | Approximation of terrestrial lead isotope evolution by a two-stage model[END_REF] and Cumming and Richards [START_REF] Cumming | Ore lead isotope ratios in a continuously changing Earth[END_REF], to correct the calculated age from the common Pb contamination. Forty-four analyses from twenty-six monazite grains (white ellipses in (a,b)) were used for calculating the age of PN16-R05a (gray ellipses (n = 6) were removed from the age calculation after SEM control, Section 4.4). Seventeen analyses from fourteen monazite grains were used for calculating the age (without removed data) from the Macôt-la Plagne samples.

Discussion

The Pb-Ag mineralizations of the MP and PN deposits are located in a similar litho-structural domain. Further, they show a main similar mineral paragenesis, deformation microstructure, and ages, suggesting that they formed under similar conditions. Amongst the mineralogical variations detected, we can notice: (i) between deposits, chalcopyrite exclusively occurs within PN; (ii) through the MP deposit (Figure 1c); higher amounts of barite are found within the Sarrazins-type veins contrary to the pod-type ores; (iii) depending on the sample, a high variability of the monazite distribution. Despite these paragenesis variations, the microstructural, thermobarometric, and geochronological similarities between the two deposits is legitimate, and the metallogenic model behind the PN and MP deposits will be further discussed indistinctly.

An Alpine Pb-Ag Mineralization

A main issue for determining metallogenic models of ore deposits located in orogenic belts resides in the possibility of inheritance or remobilization through time from the initial crystallization/deposition up to the different orogenic deformation stages (e.g., [START_REF] Cook | Sulphide Mineralogy, metamorphism and deformation in the matchless massive sulphide deposit, Namibia[END_REF][START_REF] Gu | Copper, gold and silver enrichment in ore mylonites within massive sulphide orebodies at Hongtoushan VHMS deposit[END_REF][START_REF] Turlin | Synmetamorphic Cu remobilization during the Pan-African orogeny: Microstructural, petrological and geochronological data on the kyanite-micaschists hosting the Cu(-U) Lumwana deposit in the Western Zambian copperbelt of the Lufilian belt[END_REF][START_REF] Zheng | Micro-textural and fluid inclusion data constraints on metallic remobilization of the Ashele VMS Cu-Zn deposit, Altay, NW China[END_REF]). In this study, major efforts were thus made to investigate mineral microstructures and compositional variations for discriminating crystallization steps. Whatever the disseminated or massive texture of the Pb-Ag mineralization, the mineral parageneses are very similar in all investigated samples. They display intergrown microstructures, suggesting cogenetic growth during a single hydrothermal event (Figures 4f and5). However, mineral textures are often not sufficient to correlate or decouple mineralization phases and/or fluid evolution during crystal growth between ore samples, especially without controlling the in situ structural context of the investigated samples within the deposit. Trace elements associated with precise textural observations are a powerful method to constrain the number of crystallization phases and/or fluid evolution in a mineralization stage. Galena is the only sulfide to occur in all rock samples. Galena is quite homogeneous in trace elements, except for Ag and Sb. Fine-grained galena disseminated between intergranular gangue grains or in vein-type was depleted in Ag-Sb compared to the massive galena. Intragranular zoning was observed in sulfide (pyrite), sulfosalt (tetrahedrite), silicates (chlorite and phengite, Figures 8 and9), carbonates (dolomite-ankerite group, Figure 8a), and accessory minerals (monazite, Figure 10). Whatever the minerals and their textures, the inclusions are very similar: they are generally represented by phengite and quartz, which commonly formed during the Alpine tectonics within the External Briançonnais domain [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF][START_REF] Goffé | Explanatory notes to the map: Metamorphic structure of the Alps Western and Ligurian Alps[END_REF]. More exceptional is the presence of microspheric inclusions of florencite in sample PN16-R05a. This unique feature may be due to higher and uncommon REE composition (see Supplementary Table S3) combined with the low-grade metamorphic conditions [START_REF] Janots | Evolution of the REE mineralogy in HP-LT metapelites of the Sebtide complex, rif, Morocco: Monazite stability and geochronology[END_REF]. In this PN sample, ubiquitous florencite microinclusions in recrystallized isometric and isogranular Qz2 (Figure 5c), sulfides (galena, Figure 5b), sulfosalt (tetrahedrite, Figure 5b), and monazite (Figure 5b) associated with intergrown microstructures suggest that the mineral paragenesis occurred during a single crystallization stage at relatively low-grade metamorphic conditions. Mineral zoning is thus suggested to result from the physicochemical evolution (P, T, speciation) of the fluid through time during the crystallization.

A single-stage crystallization during the Alpine orogeny is independently confirmed by monazite in situ dating. In this study, monazite is interpreted to be cogenetic of the ore mineralization based on their syngenetic texture with sulfides-sulfosalts: (i) monazite is included within galena and pyrite (Figure 5); (ii) it is intergrown with galena and tetrahedrite (Figure 5a,b); (iii) it can have a poikilitic texture and it contains the same inclusions (quartz, phengite, florencite; Figure 5a,b) as the sulfide-sulfosalts. Their elongated shape with fractionation of LREE/HREE from core to the rim is a relatively common feature of low-grade metamorphic monazite (T < 400 • C; [START_REF] Wing | Prograde destruction and formation of monazite and allanite during contact and regional metamorphism of pelites: Petrology and geochronology[END_REF][START_REF] Rasmussen | In situ U-Pb dating and element mapping of three generations of monazite: Unravelling cryptic tectonothermal events in low-grade terranes[END_REF][START_REF] Janots | Prograde metamorphic sequence of REE minerals in pelitic rocks of the central Alps: Implications for allanite-monazite-xenotime phase relations from 250 to 610[END_REF][START_REF] Seybold | A ductile extensional shear zone at the contact area between HP-LT metamorphic units in the Talea Ori, central Crete, Greece: Deformation during early stages of exhumation from peak metamorphic conditions[END_REF]). Monazite zoning is thus attributed to different stages of formation [START_REF] Rasmussen | In situ U-Pb dating and element mapping of three generations of monazite: Unravelling cryptic tectonothermal events in low-grade terranes[END_REF] or evolution in pore fluid chemistry from early diagenesis to very low-grade metamorphism [START_REF] Čopjaková | Formation of authigenic monazite-(Ce) to monazite-(Nd) from upper Carboniferous graywackes of the Drahany upland: Roles of the chemical composition of host rock and burial temperature[END_REF]. Prior to isotopic dating, preliminary EPMA U-Th-Pb dating at microscale excludes pre-Alpine ages in monazite. In the MREE-and Eu-rich core of monazite-(Nd), the exceptionally high (and Th) contents were high enough that Pb can be detected to determine an age around 33 Ma (Figure 10). For other compositions, radiogenic Pb was, however, insufficient, and common Pb too abundant for EPMA dating. The EPMA Alpine age was thus sharpened using in situ LA-ICP-MS dating, providing an age of about 35 Ma. Data are more accurate for the PN sample than for the MP samples, with a good alignment of the entire U-Pb dataset on the Tera-Wasserburg diagram, which intersects the concordia at 35.1 Ma ± 0.6 Ma. This indicates that the age of the different monazite zoning could not be distinguished using LA-ICP-MS, and that mineralization probably occurred during a time interval below the analytical uncertainties obtained from the calculated age (0.6 Ma). Based on the paragenetic sequence, this age and duration are thus assumed to be valid for the overall Pb-Ag mineralization.

Historically, textural variations in the Triassic quartzite are interpreted either as a depositional facies' variation during sedimentary deposition [START_REF] Rogel | Le Gisement de Plomb de La Plagne (Savoie)-Etude Géologique et Métallogénique-Vanoise[END_REF][START_REF] Rogel | Exemple de Gisement de Galène du Trias des Alpes[END_REF] or as a tectonic crushing [START_REF] Raguin | La mine de plomb argentifère de La Plagne (Savoie)[END_REF][START_REF] Waline | La mine de La Plagne[END_REF]. These two views have sparked the debate between two main metallogenic models for PN and MP deposits with either: a syn-diagenetic origin [START_REF] Rogel | Le Gisement de Plomb de La Plagne (Savoie)-Etude Géologique et Métallogénique-Vanoise[END_REF][START_REF] Rogel | Exemple de Gisement de Galène du Trias des Alpes[END_REF][START_REF] Bouladon | Contribution à une systématique des gisements de plomb et de zinc[END_REF][START_REF] Routhier | Géologie et Principes de Recherche[END_REF] or a tectonic cause associated with the Alpine collision [START_REF] Raguin | La mine de plomb argentifère de La Plagne (Savoie)[END_REF][START_REF] Waline | La mine de La Plagne[END_REF][START_REF] Ellenberger | Etude géologique du pays de Vanoise[END_REF][START_REF] Goffé | Explanatory notes to the map: Metamorphic structure of the Alps Western and Ligurian Alps[END_REF]]. Yet, a sedimentary/diagenetic origin can definitely be excluded based on the mineral, thermobarometry, and geochronology datasets. Such epigenetic models at relatively similar thermobarometric conditions are coherent with some recent reappreciations of other worldwide Pb-Zn(-Ag-Cu-Fe) deposits from Australian and Chinese belts [START_REF] Zhong | Ore-Forming Conditions and Genesis of the Huogeqi Cu-Pb-Zn-Fe Deposit in the Northern Margin of the North China Craton: Evidence from Ore Petrologic Characteristics[END_REF][START_REF] Cave | Linking Uplift and Mineralisation at the Mount Novit Zn-Pb-Ag Deposit, Northern Australia: Evidence from Geology, U-Pb Geochronology and Sphalerite Geochemistry[END_REF].

A Mineralizing Fluid Circulation during the D1 Deformation Stage

Field observations combined with quartz and sulfides microstructures allow reconciling the apparent stratiform texture with an ore mineralization occurring during the Alpine ductile deformation. All investigated samples show intense quartz recrystallization, up to quartzite mylonitization, in samples showing disseminated textures (Figures 3c and4a,c,d). In a single thin section (PN20-11a; Figure 4), relics of initial Qz1 porphyroclasts are unevenly preserved from recrystallization. The progressive evolution of quartz recrystallization in association with mineralization migration is schematized in Figure 12. From step 1 to step 3, the Qz1 porphyroclasts show different degrees of recrystallization by newly formed isometric and equigranular microquartz (Qz2). In the first step, Qz1 porphyroclasts are well preserved with abundant cracks generated by high-stress conditions in the presence of the mineralizing fluids, as attested by the presence of galena within the sealed cracks (Figure 4b). More generally, Qz1 recrystallization into Qz2 is mainly restricted to cracks and pressure shadows. At step 1, the sulfide-sulfosalts are mainly observed in the foliation, at porphyroclasts' grain boundaries, or associated with Qz2 micrograins in dilation zones, sealed cracks, and strain shadows. This attests to the migration of the mineralizing fluids in the foliation, at grain boundaries, and within the porosity generated by fracturing or further dissolution [START_REF] Boulter | The golden mile, Kalgoorlie; a giant gold deposit localized in ductile shear zones by structurally induced infiltration of an auriferous metamorphic fluid[END_REF]. Then, the (re)crystallization and mineralization propagate pervasively through the porphyroclast (step 2) up to reaching a complete porphyroclast recrystallization with no evidence of the initial Qz1 left (step 3). Concomitant Qz1 recrystallization with fluid propagation is underlined by the increasing sulfide and Qz2 proportion from steps 1 to 3. Truncated and elongated porphyroclast relics along the main foliation provide evidence of the action of the pressure solution in quartz deformation (S1) [START_REF] Gratier | The role of pressure solution creep in the ductility of the Earth's upper crust[END_REF]. The presence of white mica aligned parallel to the foliation may play the role of a coating and retard the pressure solution (step 1) (Renard et al. [START_REF] Renard | The different processes involved in the mechanism of pressure solution in quartz-rich rocks and their interactions[END_REF] and references therein), thus favoring porphyroclast preservation (Figure 4a). Fluid-assisted deformation is thus proposed as the main mechanism accounting for the initial fracturing and later porphyroclast recrystallization associated with sulfide precipitation following a scenario such as: Fluid-assisted deformation is thus proposed as the main mechanism accounting for the initial fracturing and later porphyroclast recrystallization associated with sulfide precipitation following a scenario such as:

(i) a stress increase, probably due to seismic pumping to reaching hydraulic fracturing [START_REF] Sibson | Implications of fault-valve behaviour for rupture nucleation and recurrence[END_REF][START_REF] Cox | The application of failure mode diagrams for exploring the roles of fluid pressure and stress states in controlling styles of fracture-controlled permeability enhancement in faults and shear zones[END_REF], which is illustrated by the abundance of cracks in the porphyroclast and the vein-type texture of the samples;

(ii) a stress drop during fluid migration through the cracks, grain boundaries, and porosity generated by quartz dissolution-reprecipitation reactions [START_REF] Renard | Kinetics of crack-sealing, intergranular pressure solution, and compaction around active faults[END_REF]. Lower stress conditions are attested by the isometric and equigranular microstructure of the newly formed Qz2 grains and the euhedral shape of the mineralization (Figures 4 and5c). It is in that stage that the sulfides-sulfosalts mineralization formed, through progressive precipitation from the mineralizing fluids propagating in the dynamically generated porosity.

Preferential orientation of the sulfides along the main foliation of the mylonite results from a syn-deformation ore genesis model associated with the Alpine collision.

P-T-t-d Correlation between the Genesis of the Pb-Ag Deposits and the Tectonics of the Briançonnais Zone

The Mine Fault (Figure 2) mylonitizes the Permian phyllitic quartzite and sericitoschist with top-to-the-NW deformation highlighted by deformed quartz porphyroclasts along the S1 foliation and the C'-type shear bands. In our mineralized samples, the mylonitic fabrics with deflection of the S1 foliations around the porphyroclast Qz1 relics is consistent with the mylonitic fabrics of nonmineralized phyllitic quartzite. Based on the similarities of the microstructures, we propose that the Pb-Ag mineralizations occur during the top-to-the-NW D1 deformation in the External Briançonnais [START_REF] Bertrand | The Pennine Front zone in Savoie (Western Alps), a review and new interpretations from the Zone Houillère Briançonnaise[END_REF].

This study constrains the main P-T-t-d conditions for ore genesis in PN and MP with an onset at 280 • C and a main stage of mineralization at 315 ± 35 • C, 6.25 ± 0.75 kbar and 35 ± 1 Ma. The P-T conditions are consistent with the Alpine peak metamorphic conditions recorded in the Briançonnais Houiller Zone (275 • C and 5.9 kbar) [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF] and with the D1-D2 transition estimation for the cover of the Internal Briançonnais (Modane-Aussois units) and the "Nappe des Gypses" (350 • C and 6.5 kbar) [START_REF] Barré | Tectono-metamorphic evolution of an evaporitic décollement as recorded by mineral and fluid geochemistry: The "Nappe Des Gypses" (Western Alps) case study[END_REF][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF]. These Alpine peak metamorphic conditions are estimated in the External Briançonnais between 34-35 MaMa by the D1-D2 transition recorded within the Modane-Aussois unit (Internal Briançonnais) [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF], during the exhumation of the Internal Briançonnais on top of the External Briançonnais through NW-directed thrusting (Figure 13). These new P-T-t-d constraints suggest that the ore formed probably due to mineralized fluids flowing along the contact between the "Nappe des Gypses" and the Permo-Triassic sediments of the External Briançonnais (Figure 13). Thus, in light of the new ore formation model, we propose that the foliation that is parallel to the Internal Briançonnais Front, which is highlighted by the contact with the "Nappe des Gypses" in MP (noted S3 on Figure 1c), must be a S1 foliation reoriented by the thrusting of the "Nappe des Gypses" on top of the External Briançonnais. The important ore concentration in both deposits illustrated by the Pb and Ag production involves an important amount of fluid. A large part of the fluid source may come from the "Nappe des Gypses". Barré et al. [START_REF] Barré | Tectono-metamorphic evolution of an evaporitic décollement as recorded by mineral and fluid geochemistry: The "Nappe Des Gypses" (Western Alps) case study[END_REF] (and references therein) indicate that the gypsum dehydration into anhydrite, probably at the end of the prograde path, might release an amount of water estimated at 37%. [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF] with associated P-T paths for Ambin basement and Vanoise cover [START_REF] Gerber | Evolution Tectono-Métamorphique du Briançonnais Interne (Alpes Occidentales, Massifs de Vanoise Sud et d'Ambin): Comportement du Socle et de sa Couverture Dans un Contexte de Subduction Continentale Profonde[END_REF], Vanoise basement, and Modane-Aussois units [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF], "Nappe des Gypses" [START_REF] Barré | Tectono-metamorphic evolution of an evaporitic décollement as recorded by mineral and fluid geochemistry: The "Nappe Des Gypses" (Western Alps) case study[END_REF] Briançonnais Houiller Zone [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF], and the PN-MP Pb-Ag deposits (this study). Geothermal gradients of 8 °C/km (a), 15 °C/km (b), and 40 °C/km (c) are also shown. [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF] with associated P-T paths for Ambin basement and Vanoise cover [START_REF] Gerber | Evolution Tectono-Métamorphique du Briançonnais Interne (Alpes Occidentales, Massifs de Vanoise Sud et d'Ambin): Comportement du Socle et de sa Couverture Dans un Contexte de Subduction Continentale Profonde[END_REF], Vanoise basement, and Modane-Aussois units [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF], "Nappe des Gypses" [START_REF] Barré | Tectono-metamorphic evolution of an evaporitic décollement as recorded by mineral and fluid geochemistry: The "Nappe Des Gypses" (Western Alps) case study[END_REF] Briançonnais Houiller Zone [START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF], and the PN-MP Pb-Ag deposits (this study). Geothermal gradients of 8

Conclusions

This integrative study shows how the complementary use of structural geology, petrology, trace element analyses in sulfides, chlorite-phengite thermobarometry, and in situ U-Pb monazite dating was successful in constraining the P-T-t-d conditions of the fluid circulation responsible for Pb-Ag deposits located in a complex and polyphase orogenic history (the Peisey-Nancroix and the Macôt-la Plagne deposits). The main conclusions are:

1.

Both Pb-Ag deposits are found in the same litho-structural context and show comparable paragenetic sequences attributed to a single mineralization episode with possible chemico-physical conditions evolution through progressive crystallization.

2.

Despite the high Pb contents, the two Pb-Ag deposits could be successfully dated using the in situ LA-ICP-MS isotopic U-Pb dating of monazite. The geochronological dataset reveals that Pb-Ag mineralization formed during the Alpine orogenesis at about 35 Ma. P-T conditions of the mineralization are estimated at 315 ± 35 • C and 6.25 ± 0.75 kbar. These correspond to the Alpine peak metamorphism of the External Briançonnais.

3.

P-T-t-d conditions, structural, and textural observation of deformed quartz porphyroclasts demonstrate a syn-kinematic Pb-Ag mineralization with the Alpine late-D1 deformation in the Internal Alps.

To our knowledge, this is the first time that an Alpine age is unambiguously proven for Pb-Ag deposits of the Western Alps. This methodology could be further applied to other low-grade metamorphic deposits from the Alps or other orogenic belts. In particular, understanding the transfer and source of the fluid, as here probably in association with the thrusting of the "Nappe des Gypses", will be of great interest to better characterize polymetallic mineralization during ductile deformation.

Figure 1 .Figure 1 .

 11 Figure 1. Geological map of (a) the Western Alps (after "Armada" France database) with square representing map b location; (b) the Briançonnais zone between the Lauzière and Vanoise massifs (modified after the Moûtiers geological map (1:50,000) by J. Debelmas). MP: Macôt-la Plagne; PN: Peisey-Nancroix; PF: Penninic Front; BF: Briançonnais Front; IBF: Internal Briançonnais Front. Additionally, sections modified after Guillot [22] of the (c) Macôt-la Plagne mine; and (d) Peisey-Figure 1. Geological map of (a) the Western Alps (after "Armada" France database) with square representing map b location; (b) the Briançonnais zone between the Lauzière and Vanoise massifs (modified after the Moûtiers geological map (1:50,000) by J. Debelmas). MP: Macôt-la Plagne; PN: Peisey-Nancroix; PF: Penninic Front; BF: Briançonnais Front; IBF: Internal Briançonnais Front. Additionally, sections modified after Guillot [22] of the (c) Macôt-la Plagne mine; and (d) Peisey-Nancroix valley. Sampling zones are given with the sample references. BHZ: Briançonnais Houiller Zone.

Figure 2 .

 2 Figure 2. (a) Photo of the Peisey-Nancroix valley (PN on the Figure 1b), Ponturin River right bank with outcrop location b (yellow star), and (b) outcrop of Permian schist and quartzite (rt) with the stereonet of the S1 schistosity and samples location (yellow stars). (c,d) Photo of the two samples from the PN outcrop showing a top-to-the-west kinematic of deformation: (c) Mylonitic quartzite with the schistosity highlighted by phyllite and the direction given by deformed porphyroclast of quartz (PN20-03); (d) Schist with shear planes (red line; C') transecting the main schistosity (S1) at an angle of about 30-60° (PN20-04). ζ: Sapey gneiss; rt: Permian schist and phyllitic quartzite; tQ: Triassic quartzite; tC: Triassic carbonates; tK: Triassic evaporites; M.F.: Mine Fault.

Figure 2 .

 2 Figure 2. (a) Photo of the Peisey-Nancroix valley (PN on the Figure 1b), Ponturin River right bank with outcrop location b (yellow star), and (b) outcrop of Permian schist and quartzite (rt) with the stereonet of the S1 schistosity and samples location (yellow stars). (c,d) Photo of the two samples from the PN outcrop showing a top-to-the-west kinematic of deformation: (c) Mylonitic quartzite with the schistosity highlighted by phyllite and the direction given by deformed porphyroclast of quartz (PN20-03); (d) Schist with shear planes (red line; C') transecting the main schistosity (S1) at an angle of about 30-60 • (PN20-04). ζ: Sapey gneiss; rt: Permian schist and phyllitic quartzite; tQ: Triassic quartzite; tC: Triassic carbonates; tK: Triassic evaporites; M.F.: Mine Fault.

36 Figure 3 .

 363 Figure 3. Macroscopic views of rock samples (a); BSE thin section image (b); thick sections under cross polars (c-f) of the various typologies of the Pb-Ag ore. (a) Massive galena in a vein crosscutting a phyllitic quartzite from a mining waste of the Peisey-Nancroix mine (PN16-R05a). (b) Fractured dolomite filled mainly by massive galena (zone 1) at the contact with fine-grained tetrahedrite-tennantite, galena, pyrite, and gangue (barite, quartz, and phengite) (zone 2) in MP sample (3501-10a). Note the deformed porphyroclasts of quartz with some chlorite in 3501-10. Images of thick sections from: (c) Thick section of quartzite with disseminated ore aligned in the schistosity (PN20-11a); (d) Quartzite with pyrite-carbonates veins (3501-5a); (e) Quartzite with massive ore (PN20-11b); (f) Massive sulfides with quartz porphyroclasts (3501-9b). Brt: barite; Ccp: chalcopyrite; Chl: chlorite; Dol: dolomite; Dol-Ank: dolomite-ankerite; Gn: galena; Py: pyrite; Qz: quartz (after Whitney and Evans [70]).

Figure 3 .

 3 Figure 3. Macroscopic views of rock samples (a); BSE thin section image (b); thick sections under cross polars (c-f) of the various typologies of the Pb-Ag ore. (a) Massive galena in a vein crosscutting a phyllitic quartzite from a mining waste of the Peisey-Nancroix mine (PN16-R05a). (b) Fractured dolomite filled mainly by massive galena (zone 1) at the contact with fine-grained tetrahedrite-tennantite, galena, pyrite, and gangue (barite, quartz, and phengite) (zone 2) in MP sample (3501-10a). Note the deformed porphyroclasts of quartz with some chlorite in 3501-10. Images of thick sections from: (c) Thick section of quartzite with disseminated ore aligned in the schistosity (PN20-11a); (d) Quartzite with pyrite-carbonates veins (3501-5a); (e) Quartzite with massive ore (PN20-11b); (f) Massive sulfides with quartz porphyroclasts (3501-9b). Brt: barite; Ccp: chalcopyrite; Chl: chlorite; Dol: dolomite; Dol-Ank: dolomite-ankerite; Gn: galena; Py: pyrite; Qz: quartz (after Whitney and Evans [70]).

Figure 4 . 36 Figure 5 .Figure 5 .

 43655 Figure 4. Examples of microstructures observed using optical microscopy. (a) Deformed Qz1 porphyroclast along the main foliation, surrounded and molded by micrometric phengite oriented in the foliation, and equigranular micrograins of Qz2 (thick section, PN20-11a, polarized light); (b) Microfractures healed and sealed by recrystallized Qz2 and sulfide-sulfosalts in a Qz1 porphyroclast (zoom of image a, polarized and reflected lights). Note that healed microcracks are cut off by newly formed Qz2 and sulfide-sulfosalts. (c,d) Qz1 porphyroclasts variably recrystallized in isometric and equigranular Qz2 micrograins and sulfides (galena and pyrite) in polarized light (c) and reflected light ((d); PN20-11a; thick section). Note that, when absent, the initial shape of Qz1 porphyroclast is

Figure 6 .

 6 Figure 6. A common paragenesis for the Macôt-la Plagne and the Peisey-Nancroix Pb-Ag deposits. The thickness of the line represents the proportion of each phase.

Figure 6 .

 6 Figure 6. A common paragenesis for the Macôt-la Plagne and the Peisey-Nancroix Pb-Ag deposits. The thickness of the line represents the proportion of each phase.

Figure 7 .

 7 Figure 7. Ag versus Sb + Bi plot for trace elements in galena from MP samples (3501-10a-1 and 3501-10a-2 correspond, respectively, to galena from zone 1 and 2 on Figure3b) and PN samples. Circles and diamonds represent LA-ICP-MS and EPMA analyses, respectively. The correlation line was drawn using all the dataset; the correlation coefficient is given by R².

Figure 7 .

 7 Figure 7. Ag versus Sb + Bi plot for trace elements in galena from MP samples (3501-10a-1 and 3501-10a-2 correspond, respectively, to galena from zone 1 and 2 on Figure3b) and PN samples. Circles and diamonds represent LA-ICP-MS and EPMA analyses, respectively. The correlation line was drawn using all the dataset; the correlation coefficient is given by R 2 .

Figure 8 .

 8 Figure 8. (a) Backscatter (BSE) electron image of euhedral rosette-shape chlorite aggregates with chemical zoning. The red frame displays the zone of the X-ray map. Quantified X-ray chlorite maps (in apfu) of: (b) Mn; (c) Si; (d) Al; (e) XMg content (XMg = Mg/(Mg + Fetot)) with black boxes showing the pixel zone analyzed of each representative chlorite zone, as discussed in the text. (f) Histogram of chlorite temperature formation calculated on each pixel from the map.

Figure 8 .

 8 Figure 8. (a) Backscatter (BSE) electron image of euhedral rosette-shape chlorite aggregates with chemical zoning. The red frame displays the zone of the X-ray map. Quantified X-ray chlorite maps (in apfu) of: (b) Mn; (c) Si; (d) Al; (e) XMg content (XMg = Mg/(Mg + Fe tot )) with black boxes showing the pixel zone analyzed of each representative chlorite zone, as discussed in the text. (f) Histogram of chlorite temperature formation calculated on each pixel from the map.

Figure 9 .

 9 Figure 9. (a,b) Quantified X-ray phengite maps (in apfu) of the Si content (a) and the Al content (b).The black boxes show the pixel zones analyzed. (c) Phengite P-T estimates (lines) using the method of Dubacq et al.[START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: Thermodynamic modelling and implications for thermobarometric estimates[END_REF]. The box is constructed using the result of Chl3 thermometry (Figure8f) and the punctual data. For a given temperature, pressure has been estimated using the phengite P-T equilibrium lines of the punctual data from the map area. Boundaries of the Ph1 and Ph2 were determined using more than 95% of P-T lines for each zone.

Figure 9 .

 9 Figure 9. (a,b) Quantified X-ray phengite maps (in apfu) of the Si content (a) and the Al content (b). The black boxes show the pixel zones analyzed. (c) Phengite P-T estimates (lines) using the method of Dubacq et al.[START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: Thermodynamic modelling and implications for thermobarometric estimates[END_REF]. The box is constructed using the result of Chl3 thermometry (Figure8f) and the punctual data. For a given temperature, pressure has been estimated using the phengite P-T equilibrium lines of the punctual data from the map area. Boundaries of the Ph1 and Ph2 were determined using more than 95% of P-T lines for each zone.

Geosciences 2022 ,

 2022 12, x FOR PEER REVIEW 24 of 36

Figure 10 .

 10 Figure 10. Elementary maps (La, Th, Ce, U, Nd, Si, Eu, Ca, Y) of a Peisey-Nancroix monazite surrounded by galena (white) and quartz (black; PN16-R05a), performed using electron microprobe and BSE image showing electron microprobe ages (white dots). The violet circles show LA-ICP-MS spots of U-Th-Pb dating.

Figure 10 .

 10 Figure 10. Elementary maps (La, Th, Ce, U, Nd, Si, Eu, Ca, Y) of a Peisey-Nancroix monazite surrounded by galena (white) and quartz (black; PN16-R05a), performed using electron microprobe and BSE image showing electron microprobe ages (white dots). The violet circles show LA-ICP-MS spots of U-Th-Pb dating.

  MP and PN mineralizing fluids are close to the predicted common Pb values, opposite to some mineralizing fluids in relation with the Penninic front[START_REF] Rossi | Geochemical Variations and Element Transfer during Shear-Zone Development and Related Episyenites at Middle Crust Depths: Insights from the Mont Blanc Granite (French-Italian Alps)[END_REF][START_REF] Cenki-Tok | Direct Dating of Mid-Crustal Shear Zones with Synkinematic Allanite: New in Situ U-Th-Pb Geochronological Approaches Applied to the Mont Blanc Massif[END_REF]. Moreover, even if monazite grains from three different rock samples are considered for the MP deposit, there is one sample (3501-5b) more represented and well dispersed along the Discordia line, which gives an age of 32.15 ± 3.25 Ma (n = 11) within the uncertainty of the age obtained by considering the entire MP dataset (Figure11d).

Figure 11 .

 11 Figure 11. Tera-Wasserburg diagrams of the monazite LA-ICP-MS analysis (Table 3) for the Peisey-Nancroix (a,b) and the Macôt-la Plagne (c,d) deposits. (a) and (c) provides Discordia lines and ages calculated without any correction, whereas for (b) and (d) the initial 207 Pb/ 206 Pb ratio is fixed at 0.835, using the bulk-Pb evolution model of Stacey and Kramers[START_REF] Stacey | Approximation of terrestrial lead isotope evolution by a two-stage model[END_REF] and Cumming and Richards[START_REF] Cumming | Ore lead isotope ratios in a continuously changing Earth[END_REF], to correct the calculated age from the common Pb contamination. Forty-four analyses from twenty-six monazite grains (white ellipses in (a,b)) were used for calculating the age of PN16-R05a (gray ellipses (n = 6) were removed from the age calculation after SEM control, §4.4.). Seventeen analyses from fourteen monazite grains were used for calculating the age (without removed data) from the Macôtla Plagne samples.

Figure 11 .

 11 Figure 11. Tera-Wasserburg diagrams of the monazite LA-ICP-MS analysis (Table3) for the Peisey-Nancroix (a,b) and the Macôt-la Plagne (c,d) deposits. (a,c) provides Discordia lines and ages calculated without any correction, whereas for (b,d) the initial 207 Pb/ 206 Pb ratio is fixed at 0.835, using the bulk-Pb evolution model of Stacey and Kramers[START_REF] Stacey | Approximation of terrestrial lead isotope evolution by a two-stage model[END_REF] and Cumming and Richards[START_REF] Cumming | Ore lead isotope ratios in a continuously changing Earth[END_REF], to correct the calculated age from the common Pb contamination. Forty-four analyses from twenty-six monazite grains (white ellipses in (a,b)) were used for calculating the age of PN16-R05a (gray ellipses (n = 6) were removed from the age calculation after SEM control, Section 4.4). Seventeen analyses from fourteen monazite grains were used for calculating the age (without removed data) from the Macôt-la Plagne samples.

Figure 12 .

 12 Figure 12. Scheme of quartz (Qz1) and textural evolution of the ore in mylonites due to increasing deformation (Step 1 to 3) and fluid circulation: from disseminated ore (step 1 and 2) to massive ore (step 3).

Figure 12 .

 12 Figure 12. Scheme of quartz porphyroclast (Qz1) and textural evolution of the ore in mylonites due to increasing deformation (Step 1 to 3) and fluid circulation: from disseminated ore (step 1 and 2) to massive ore (step 3).

Figure 13 .

 13 Figure 13. Possible kinematic reconstruction this part of the Western Alps from Eocene to Oligocene modified after Lanari et al.[START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF] with associated P-T paths for Ambin basement and Vanoise cover[START_REF] Gerber | Evolution Tectono-Métamorphique du Briançonnais Interne (Alpes Occidentales, Massifs de Vanoise Sud et d'Ambin): Comportement du Socle et de sa Couverture Dans un Contexte de Subduction Continentale Profonde[END_REF], Vanoise basement, and Modane-Aussois units[START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF], "Nappe des Gypses"[START_REF] Barré | Tectono-metamorphic evolution of an evaporitic décollement as recorded by mineral and fluid geochemistry: The "Nappe Des Gypses" (Western Alps) case study[END_REF] Briançonnais Houiller Zone[START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF], and the PN-MP Pb-Ag deposits (this study). Geothermal gradients of 8 °C/km (a), 15 °C/km (b), and 40 °C/km (c) are also shown.

Figure 13 .

 13 Figure 13. Possible kinematic reconstruction of this part of the Western Alps from Eocene to Oligocene modified after Lanari et al.[START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF] with associated P-T paths for Ambin basement and Vanoise cover[START_REF] Gerber | Evolution Tectono-Métamorphique du Briançonnais Interne (Alpes Occidentales, Massifs de Vanoise Sud et d'Ambin): Comportement du Socle et de sa Couverture Dans un Contexte de Subduction Continentale Profonde[END_REF], Vanoise basement, and Modane-Aussois units[START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF], "Nappe des Gypses"[START_REF] Barré | Tectono-metamorphic evolution of an evaporitic décollement as recorded by mineral and fluid geochemistry: The "Nappe Des Gypses" (Western Alps) case study[END_REF] Briançonnais Houiller Zone[START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF], and the PN-MP Pb-Ag deposits (this study). Geothermal gradients of 8 • C/km (a), 15 • C/km (b), and 40 • C/km (c) are also shown.

  Figure 13. Possible kinematic reconstruction of this part of the Western Alps from Eocene to Oligocene modified after Lanari et al.[START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF] with associated P-T paths for Ambin basement and Vanoise cover[START_REF] Gerber | Evolution Tectono-Métamorphique du Briançonnais Interne (Alpes Occidentales, Massifs de Vanoise Sud et d'Ambin): Comportement du Socle et de sa Couverture Dans un Contexte de Subduction Continentale Profonde[END_REF], Vanoise basement, and Modane-Aussois units[START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect[END_REF], "Nappe des Gypses"[START_REF] Barré | Tectono-metamorphic evolution of an evaporitic décollement as recorded by mineral and fluid geochemistry: The "Nappe Des Gypses" (Western Alps) case study[END_REF] Briançonnais Houiller Zone[START_REF] Lanari | Diachronous evolution of the alpine continental subduction wedge: Evidence from P-T estimates in the Briançonnais Zone Houillère (France-Western Alps)[END_REF], and the PN-MP Pb-Ag deposits (this study). Geothermal gradients of 8 • C/km (a), 15 • C/km (b), and 40 • C/km (c) are also shown.

Table 1 .

 1 Sulfides-sulfosalts major element concentrations determined by EPMA in wt%, and trace elements by LA-ICP-MS in ppm.

	Mineral	Sample	Pb (wt%)	Cu	Fe	Zn	As		Ag		Sn	Sb	Bi
	Galena	3501-5a *
		Min.	85.31	-	-	-	-		260		-	369	-
		Max.	86.93	-	179	-	-		1003		152	912	-
		Mean (17)	86.2	-	20	-	-		588		8.9	548	-
		1σ	0.52	-	56	-	-		175		37	160	-
		3501-9b
		Min.	86.52	1.3	0.2	1.1	16		1340		0.2	1670	0.1
		Max.	87.17	19	11	1.6	37		1665		0.4	2081	0.9
		Mean (13)	86.8	9.8	3.6	1.3	30		1507		0.2	1866	0.4
		1σ	0.33	4.3	4.3	0.3	5.8		92		0.1	125	0.2
		3501-10a
		Min.		0.6	1.1	0.4	3.0		481		0.2	585	0.1
		Max. Mean (11)	86.0	18 8.9	3113 493	5.3 1.9	16 13		1447 957		0.3 0.2	1729 1165	0.5 0.3
		1σ		6.8	1157	1.5	4.0		374		0.0	446	0.1
		PN16-R05a *
		Min.	85.49	-	-	-	-		369		-	277	-
		Max.	86.53	-	203	-	-		537		-	694	-
		Mean (5)	86.0	-	41	-	-		449		-	524	-
		1σ	0.39	-	91	-	-		68		-	156	-
		PN20-11
		Min.		229	0.5	0.1	38		246		0.3	245	9.2
		Max. Mean (6)	86.0	259 244	249 53	0.9 0.5	40 39		369 297		0.7 0.5	416 314	55 30
		1σ		22	110	0.5	1.0		48		0.2	62	22
		PN20-12
		Min.		-	1.5	0.7	17		174		0.9	111	12
		Max. Mean (14)	86.0	66 6.9	26 8.4	3.6 1.5	75 38		328 240		2.0 1.2	289 191	61 31
		1σ		21	9.7	0.9	22		47		0.3	47	19
	Mineral	Sample	Pb	Cu (wt%)	Fe (wt%)	Zn (wt%)	As (wt%)	Ag (wt%)	Sn	Sb (wt%)	Bi
	Tetrahedrite	PN16-R05a *
	-Tennantite	Min.	-	32	2.8	3.8	3.2		5.3		1132	21	-
		Max.	-	35	3.5	4.4	5.6		6.7		1370	25	-
		Mean (10)	-	34	3.1	4.1	4.5		6.1		1232	23	-
		1σ	-	0.9	0.2	0.2	1.0		0.4		100	1.5	-
		PN20-12
		Min.	0.3		2.2	6.2	20		1.2		0.1	3.8	0.8
		Max. Mean (5)	2.4 1.4	35.0	2.3 2.2	6.3 6.3	23 22		1.5 1.3		2.5 0.6	5.2 4.7	2.5 1.3
		1σ	0.8		0.1	0.1	1.5		0.1		1.1	0.6	0.7
	Mineral	Sample	Pb	Cu (wt%)	Mn	Fe (wt%)	Zn	As	Ag	Cd	Sn	Sb	Bi
	Chalcopyrite	PN20-11
		Min.	0.4		0.0	29	46	48	0.3	6.1	6.0	1.3	2.9
		Max. Mean (8)	21 7.6	34.0	10 2.1	30 29	129 62	92 67	101 20	20 9.8	8.9 7.1	89 24	2.9 2.9
		1σ	12		4.5	0.4	28	13	34	5.0	0.9	34	-
			Bold values of Pb and Cu correspond to internal standards obtained using EPMA. * Samples with trace element
			concentrations determined by EPMA.

Table 2 .

 2 Representative EPMA chlorite and phengite compositions for the MP and PN ores. Standard deviations (σ) are expressed within parentheses. Oxides are expressed in wt% and elementary atoms per formula unit (apfu). Sum VI corresponds to the octahedral sum.

	Samples						3501-5b (MP)						PN16-R05a
	Mineral	Chl1		Chl2		Chl3			Ph1		Ph2		Ph
	SiO 2	26.70	(0.60)	26.31	(0.54)	26.52	(0.66)		49.09	(1.11)	49.73	(0.72)	48.59	(1.13)
	Al 2 O 3	19.66	(0.31)	19.52	(0.33)	20.54	(0.31)		26.57	(0.52)	25.80	(0.47)	31.20	(0.72)
	FeO	27.89	(0.49)	28.61	(0.56)	25.84	(0.61)		4.70	(0.40)	4.46	(0.30)	1.15	(0.10)
	MgO	14.09	(0.29)	13.00	(0.21)	15.23	(0.55)		2.39	(0.11)	2.42	(0.08)	2.00	(0.13)
	MnO	0.35	(0.10)	0.56	(0.08)	0.35	(0.08)		0.04	(0.06)	0.04	(0.01)	-
	K 2 O	-		-		-			10.19	(0.22)	10.04	(0.21)	10.83	(0.20)
	Na 2 O	0.03	(0.01)	0.03	(0.01)	0.03	(0.01)		0.35	(0.22)	0.43	(0.18)	0.14	(0.06)
	CaO	0.03	(0.02)	0.02	(0.02)	0.04	(0.03)		0.04	(0.06)	0.08	(0.08)	-
	TiO 2	0.01	(0.03)	0.01	(0.03)	0.04	(0.06)		0.02	(0.01)	0.02	(0.00)	0.14	(0.06)
	Total	88.76	(0.84)	88.04	(0.82)	88.58	(1.08)		93.38	(1.19)	93.03	(0.74)	94.11	(1.21)
	apfu			(14 anhydrous oxygen basis)					(11 anhydrous oxygen basis)
	Si	2.81	(0.04)	2.81	(0.04)	2.77	(0.04)		3.38	(0.04)	3.43	(0.03)	3.26	(0.04)
	Al	2.44	(0.05)	2.46	(0.03)	2.53	(0.05)		2.16	(0.05)	2.10	(0.04)	2.47	(0.05)
	Fe	2.46	(0.04)	2.56	(0.05)	2.26	(0.06)		0.27	(0.02)	0.26	(0.02)	0.06	(0.01)
	Mg	2.21	(0.04)	2.07	(0.03)	2.37	(0.07)		0.25	(0.01)	0.25	(0.01)	0.20	(0.02)
	Mn	0.03	(0.01)	0.05	(0.01)	0.03	(0.01)		-		-		-
	K	-		-		-			0.90	(0.02)	0.88	(0.02)	0.93	(0.02)
	Na	-		-		-			0.05	(0.03)	0.06	(0.02)	0.02	(0.01)
	Ca	-		-		-			-		0.01	(0.01)	-
	Ti	-		-		-			-		-		0.01	(0.00)
	∑VI	5.97	(0.03)	5.96	(0.03)	5.97	(0.03)	Xcel	0.40	(0.04)	0.43	(0.02)	0.25	(0.02)
	XMg	0.47	(0.01)	0.45	(0.01)	0.51	(0.01)	XMs	0.50	(0.04)	0.47	(0.03)	0.68	(0.03)
	Xfe 3+	32	(5)	30	(4)	27	(5)	XPrl	0.05	(0.03)	0.05	(0.02)	0.05	(0.03)
	T ( • C) *	285	(30)	280	(30)	315	(35)	XPg	0.04	(0.02)	0.05	(0.02)	0.02	(0.01)
				* Vidal et al. [40,41]; Pinit = 5 kbar.

Geosciences 2022, 12, x FOR PEER REVIEW

Acknowledgments:

The authors gratefully acknowledge Roger De Ascenção Guedes and Serge Fudral for their precious knowledge of the Savoie base metal occurrences and on the regional geology. We thank Laurent Truche for his precious advice, and Ahmed Nouibat for his help on MatLab software. The three anonymous reviewers helped clarifying several aspects of the paper and are sincerely acknowledged.

Funding: This work was first supported by the "Référentiel Géologique de la France" RGF-Alps Program by BRGM, which funds the doctoral grant of M. Bertauts. This work was further supported by the TelluS Program of CNRS/INSU, the IDEX-IRGA 2021, and local BQR funding. The work was mainly done in ISTerre (Université Grenoble Alpes), which is part of Labex OSUG@2020 (ANR10 LABX56). The sulfide analyses were carried out at the LA-ICP-MS laboratory of GeoRessources in Nancy, which is funded by the Labex Ressources 21 (ANR-10-LABX-21-RESSOURCES21), the Région Lorraine, and the European Community through the FEDER program.

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.

Conflicts of Interest:

The authors declare no conflict of interest.