Guillaume Barré
email: guillaume.barre@univ-pau.fr

Pierre Strzerzynski

Raymond Michels

Stephane Guillot

Pierre Cartigny

Emilie Thomassot

Catherine Lorgeoux

Nelly Assayag

Laurent Truche

Stéphane Guillot

Tectono-metamorphic evolution of an evaporitic décollement

Keywords: Alps, Evaporites, Décollement, Metamorphism, Fluid inclusions, Stable Isotopes analysis

come

Introduction

Evaporites are found all over the world in fold and thrust belts including among others: Alps (e.g., [START_REF] Malavieille | Mylonitic deformation of evaporites in décollements: examples from the Southern Alps, France[END_REF][START_REF] Graham | Allochthonous salt in the sub-Alpine fold-thrust belt of Haute Provence, France[END_REF], Himalaya (e.g., [START_REF] Baker | Development of the Himalayan frontal thrust zone: Salt Range, Pakistan[END_REF], Pyrénées (e.g., [START_REF] Canérot | Mesozoic diapirism in the Pyrenean orogen: Salt tectonics on a transform plate boundary[END_REF], Turkey (e.g., [START_REF] Kergaravat | Minibasins and salt canopy in foreland fold-and-thrust belts: The central Sivas Basin, Turkey[END_REF], Zagros (e.g., [START_REF] Callot | Pre-existing salt structures and the folding of the Zagros Mountains[END_REF]. Evaporites are well documented in many orogenic structures as remarkable décollement layers [START_REF] Davis | The role of salt in fold-and-thrust belts[END_REF][START_REF] Vendeville | How the presence of a salt décollement in the sedimentary cover influences the behavior of subsalt thrusts in fold-and-thrust belts[END_REF]. Geomechanical experiments and models both under compressive (e.g., [START_REF] Costa | Experimental insights on the geometry and kinematics of fold-and-thrust belts above weak, viscous evaporitic décollement[END_REF][START_REF] Gemmer | Dynamic modelling of passive margin salt tectonics: effects of water loading, sediment properties and sedimentation patterns[END_REF][START_REF] Dooley | The structure and evolution of sutures in allochthonous salt[END_REF] and extensive regimes (e.g., [START_REF] Vendeville | Scale models of salt tectonics during basementinvolved extension[END_REF] have brought a series of key information on its mechanism. For example, it is shown that excess of pore pressure induced by the dehydration of gypsum into bassanite and anhydrite (e.g., [START_REF] Ko | Generation and maintenance of pore pressure excess in a dehydrating system 1. Experimental and microstructural observations[END_REF][START_REF] Leclère | Permeability control on transient slip weakening during gypsum dehydration: Implications for earthquakes in subduction zones[END_REF], are related to occurrences of earthquakes. However, the overall role of evaporites during subduction followed by exhumation is not fully understood due to the lack of major silicates assemblage in evaporitic units. This study aims at unravelling the motion of such a décollement layer and the processes responsible for the weakening of evaporites in this peculiar context by taking the "Nappe des Gypses" evaporitic formation (Western Alps) as a case study.

In the Alps the "Nappe des Gypses" evaporitic formation is described as a major décollement layer (e.g., [START_REF] Malavieille | Mylonitic deformation of evaporites in décollements: examples from the Southern Alps, France[END_REF] between the external and the internal Briançonnais zones and is considered as having played a major role in nappe formation and hence structuration of major tectonic units. The "Nappe des Gypses" formation was specifically studied for its geotechnical properties [START_REF] Fabre | Proprietes geotechniques de gypses et anhydrites du trias des alpes de savoie (France)[END_REF], composition of fluid inclusions [START_REF] Grappin | Les variations des rapports Cl/Br, Na/Br et K/Br dans les inclusions fluides des quartz de la cicatrice évaporitique de Bramans-Termignon (Vanoise, Alpes françaises)[END_REF] and for saline water springs related to the evaporites [START_REF] Zuppi | Eaux fortement minéralisées et circulations profondes dans le socle. Exemple des Alpes franco-italiennes[END_REF]. This formation is either studied as an independent unit (e.g., [START_REF] Jaillard | Evolutions sédimentaire et paléotectonique de la zone briançonnaise de Vanoise occidentale (Alpes occidentales françaises)[END_REF] or integrated as part of the Briançonnais cover [START_REF] Gabalda | Thermal structure of a fossil subduction wedge in the Western Alps[END_REF], the Vanoise cover [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] or the Modane-Aussois unit [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF]. Whatever its tectonic association, this layer always occurs structurally below the Jurassic/Cretaceous Liguro-Piemontese units.

Even if the internal Alps have been widely studied during the late decades, only few results focused on the "Nappe des Gypses". Only one temperature estimates of 420°C has been reported so far for the metamorphic peak [START_REF] Gabalda | Thermal structure of a fossil subduction wedge in the Western Alps[END_REF] and no pressure estimate is presently available. This lack of data is probably due to the nature of dominant rock lithology (anhydrite and dolomite) which is not optimal for structural, metamorphic or geochronology studies. This study proposes for the first time a multi-technique investigation documenting the petrography, the deformation features, pressure and temperature history of the "Nappe des Gypses" in order to decipher its burial/exhumation evolution within the framework of the western Alps geological history.

Geological settings

General structural context of the studied area

The internal Alps consists of a stack of continental and oceanic-margin derived nappes emplaced and metamorphosed during the Cenozoic (Fig. 1; e.g. de [START_REF] De Graciansky | The Western Alps, from rift to passive margin to orogenic belt: an integrated geoscience overview[END_REF]. The "Nappe des Gypses" unit corresponds to hundreds of meters (from 100 to 400 m) of dominant massive anhydrite interbedded with dolomites and micaschists layers from centimetric to plurimetric thickness (Fig. 2). The "Nappe des Gypses" formation corresponds to a typical marine basin wide evaporites depositions [START_REF] Warren | Evaporites through time: Tectonic, climatic and eustatic controls in marine and nonmarine deposits[END_REF] characterized by large sulfates accumulation during an extension phase of late Triassic (Carnian) as evidenced by the presence of Equisetum fossils in schists layers embedded in the evaporites [START_REF] Fudral | Notice explicative de la feuille Lanslebourg-Mont d[END_REF]. The "Nappe des Gypses" unit is covered by Liguro-Piemontese units. It consists of ophiolitic assemblages derived from Upper Jurassic oceanic crust (metagabbro, metabasalts) and mantle remnants characterized by serpentinites embedded into Cretaceous metasediments called "Schistes Lustrés" (e.g., [START_REF] Agard | Exhumation of the Schistes Lustrés complex: in situ laser probe 40Ar/39Ar constraints and implications for the Western Alps[END_REF][START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[END_REF]. The contact between the "Nappe des Gypses" and the Liguro-Piemontese units is of tectonic origin induced by the presence of evaporitic rocks that favored strain localization (e.g., [START_REF] Davis | The role of salt in fold-and-thrust belts[END_REF].

The Briançonnais zone corresponds to a continental basement covered by Carboniferous to Eocene deposits (Strzerzynski et al., 2012 and references therein). To the east, lower to middle Triassic detritic deposits, schists and carbonates are the most abundant whereas to the west, carboniferous conglomerates, sandstones, black shales and anthracite form the "Zone Houillère" (e.g., [START_REF] Mercier | Révision du carbonifère briançonnais: Stratigraphie et évolution du bassin[END_REF][START_REF] Barféty | Notice explicative de la feuille Névache-Bardonecchia-Modane à 1/50,000[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]. The hundred meters thick Liassic carbonates and calcschists of the Dent Parrachée unit forms the syn-rift series announcing the opening of the Liguro-Piemontese oceanic domain (e.g., [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]. This unit has been sometime associated with the "Nappe des Gypses" unit [START_REF] Gabalda | Thermal structure of a fossil subduction wedge in the Western Alps[END_REF].

As the evaporites of the "Nappe des Gypses" unit have sedimented after the early to middle Triassic Briançonnais deposits, it has been often proposed that the "Nappe des Gypses" forms the direct continuation of the Briançonnais units [START_REF] Mascle | Salt tectonics, Tethyan rifting and Alpine folding in the French Alps[END_REF]. However, the nature of the contact between the "Nappe des Gypses" and the Briançonnais units is still debated.

Regional tectonic and metamorphic evolution

The internal part of the Western Alps experienced a complex polyphased tectonic evolution that takes place in a general context of subduction/exhumation. Three tectonic phases are distinguished and called thereafter the D1, D2 and D3 events. Even if these events present structural similarities, recent studies have highlighted several changes in the relationships between tectonic events and metamorphic P-T path (e.g., [START_REF] Schwartz | Late tectonic and metamorphic evolution of the Piedmont accretionary wedge (Queyras Schistes lustrés, western Alps): Evidences for tilting during Alpine collision[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]. In addition, precise datings highlight a younging of successives tectono-metamorphic events toward the west [START_REF] Agard | Exhumation of the Schistes Lustrés complex: in situ laser probe 40Ar/39Ar constraints and implications for the Western Alps[END_REF][START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[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][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF].

1) The D1 event consists of a top-to-the north to north-west thrusting phase characterized by intense deformations, isoclinal folding and formation of the D1 foliation oriented sub-parallel to S0 [START_REF] Agard | Exhumation of the Schistes Lustrés complex: in situ laser probe 40Ar/39Ar constraints and implications for the Western Alps[END_REF][START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[END_REF][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[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]. It is responsible of the formation of the internal nappe stack of the Liguro-Piemontese unit over the European margin. Within the Briançonnais units, it is also responsible for the formation of several tectonic slices such as the internal crystalline massifs, the internal and external Briançonnais units [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[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]. The D1 event occurred during the metamorphic peak and the early stage of exhumation. P-T conditions are up to 570°C, 22-24 kbar in the Liguro-Piemontese unit [START_REF] Agard | Tectonometamorphic evolution of the Schistes Lustres Complex; implications for the exhumation of HP and UHP rocks in the Western Alps[END_REF][START_REF] Plunder | How continuous and precise is the record of P-T paths? Insights from combined thermobarometry and thermodynamic modelling into subduction dynamics (Schistes Lustrés, W. Alps)[END_REF], 500°C and 15 kbar in the internal Briançonnais zone [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[END_REF][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF] and 270 -350°C and 8 to 10 kbar in the external Briançonnais 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][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF]. Time constrains of the D1 event indicate that it occurred first at ca. 66-55 Ma within the oceanic units [START_REF] Agard | Exhumation of the Schistes Lustrés complex: in situ laser probe 40Ar/39Ar constraints and implications for the Western Alps[END_REF], then into the internal crystalline massifs and the internal part of the Briançonnais units between 55 and 43 Ma [START_REF] Meffan-Main | A Tertiary age for Alpine high-pressure metamorphism in the Gran Paradiso massif, Western Alps: A Rb-Sr microsampling study[END_REF][START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[END_REF][START_REF] Rosenbaum | Dating deformation in the Gran Paradiso Massif (NW Italian Alps): Implications for the exhumation of high-pressure rocks in a collisional belt[END_REF] and finally into the external Briançonnais units between 45 and 35 Ma [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF].

2) The D2 event is characterized by the top-to-the east shear zone associated with large folds.

A S2 schistosity developed under greenschist facies and affected the whole D1 nappe stack [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[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, Western Alps)[END_REF]. The tectonic significance of the D2 event is still debated. In the most internal part of the belt, top-to-the east D2 shear is frequently parallel to the D1 foliation and interpreted as detachment [START_REF] Agard | Exhumation of the Schistes Lustrés complex: in situ laser probe 40Ar/39Ar constraints and implications for the Western Alps[END_REF][START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[END_REF]. In the external Briançonnais units, the top-to-the east event is related to backthrusting and backfolding in a context of thickening of the belt [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]. Both extensional D2 detachment and backthrusting occurred at the transition between internal and external Briançonnais units and are very difficult to distinguish [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF]. Time constrains suggest that extensional top-to-the east shear zone occurred earlier within the oceanic unit (50-40 Ma; [START_REF] Agard | Exhumation of the Schistes Lustrés complex: in situ laser probe 40Ar/39Ar constraints and implications for the Western Alps[END_REF], and then migrated into the internal crystalline massif Ma; e.g., [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[END_REF]. Ages of the backthrusting are comprised between 35 and 32 Ma in the Briançonnais units (e.g., [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF].

3) The D3 event corresponds to the formation of large dome into the nappe stack in association with faulting [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[END_REF]. It is responsible of the last stage of exhumation below 300°C accompanied by a top-to-the west tilting phase within the Liguro-Piemontese unit [START_REF] Schwartz | Late tectonic and metamorphic evolution of the Piedmont accretionary wedge (Queyras Schistes lustrés, western Alps): Evidences for tilting during Alpine collision[END_REF] in a context of polyphase brittle tectonics [START_REF] Champagnac | Miocene lateral extrusion in the inner western Alps revealed by dynamic fault analysis[END_REF].

Deformations occurred during the Miocene and probably are still active today [START_REF] Delacou | Present-day geodynamics in the bend of the western and central Alps as constrained by earthquake analysis[END_REF].

Sampling and analytical methods

Sampling

The "Nappe des Gypses" formation was sampled from 6 major outcrops in the internal part of the western French Alps. They were selected to cover most of the geographical extent along the tectonic contact with the Briançonnais units, near the villages of Modane, Bramans, and Sollières-l'Envers (Arc Valley), the Mont-Cenis lake, the Roubion stream (near Névache village) and Tignes lake (Fig. 1). Lithologies as well as structural features were recognized as to bring geological background to rocks facies and mineral assemblages (e.g. in fractures).

Petrography

Hand samples as well as thin and thick sections were studied for petrography, mineralogy and fluid inclusions. In addition, X-ray Diffraction (XRD) analyses were performed on micaschists layers sampled on each studied site to determine their mineralogy using a ®Bruker D2 Phaser diffractometer, with Cu-Kα radiation. Diffractograms were interpreted using EVA software and spectra comparison to the NIST database. Micas may not be distinguished from their clay alteration products using solely XRD spectra because of similar chemical and crystallographic properties. To ensure the identification of unaltered micas, elemental compositions of micaschists layers were obtained using a Cameca computer-controlled SX-100 Electron Microprobe (EMP) equipped with a wavelength dispersive spectrometer (WDS). The accelerating voltage was set at 15 kV, beam current at 12 nA and peak counting time varied from 10 to 20s for silicates. The EMP analyses were used to determine pressure in the "Nappe des Gypses" formation using the phengite thermobarometer of [START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: thermodynamic modelling and implications for thermobarometric estimates[END_REF]. This technique is based on the K-white micas hydration calculated on the basis of their elemental composition determined by EMP. As the host-rock corresponds to gypsum dehydration into anhydrite [START_REF] Grappin | Les variations des rapports Cl/Br, Na/Br et K/Br dans les inclusions fluides des quartz de la cicatrice évaporitique de Bramans-Termignon (Vanoise, Alpes françaises)[END_REF], we consider that water is available in abundance in the system. Therefore, we fixed the water activity at 1. This allows a simultaneous estimation of the pressure and nH2O at a given temperature. The method of [START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: thermodynamic modelling and implications for thermobarometric estimates[END_REF] to estimate pressure requires that the temperatures were already determined on chlorite associated with phengite. As no chlorite where found to be associated with the phengite, we used the temperatures previously determined by Raman Spectroscopy of Carbonaceous Material (RSCM), fluid inclusions and stable isotopes.

Raman Spectroscopy of Carbonaceous Material (RSCM)

The maximum temperature of metamorphism was deduced using the geothermometer of [START_REF] Beyssac | Raman spectra of carbonaceous material in metasediments: a new geothermometer[END_REF]. This technique is based on the variation of the organization of carbonaceous material and the integration of the Raman band area typical of graphite. Raman spectra were recorded on a Labram spectrometer ( ® Jobin-Yvon, Horiba) using a 514 nm (green) Ar + laser excitation with a laser power of 200 mW. A filter was used to reduce the laser power on the sample at around 1 mW to prevent samples degradation during acquisition. Each acquisition was performed using a grating of 1800 lines/mm, a slit width of 100 μm and a confocal hole of 500 μm providing a spectral resolution of 1 to 1.5 cm -1 .

Fluid inclusions analysis

Fluid inclusions were studied on double polished 200 μm-thick sections from samples of i) anhydrite ± quartz ± fluorite veins present in the dolomitic "boudins", ii) automorph albites mineral hosted in some dolomitic "boudins" and iii) quartz veins associated to micaschists layers (see section 4.2. below for petrological description).

Microthermometric measurements were performed on a Linkam THMSG600 heatingcooling stage connected to an Olympus BX51 microscope. The stage was calibrated using

Carbon and oxygen isotopes analysis of dolomite

Carbon and oxygen isotopic compositions of dolomite were determined by using an auto sampler Gasbench coupled to a Thermo Scientific MAT253 isotope ratio mass spectrometer (IRMS). Each mineral phase was collected as 2 to 1000 mg of powder using a Dremel tool with boreholes of 1 and 2 mm in diameter. This procedure allows better targeting the desired mineral phases of <1 mm diameter whatever the host environment. Samples powder was then reacted with 2 mL of supersaturated orthophosphoric acid at 70°C during 5h in a He atmosphere. 10 measurements cycles of the isotopic composition of produced CO2 were performed. All samples measurements were adjusted to an internal reference calibrated on the international standards IAEA CO-1, IAEA CO-8 [START_REF] Stichler | Interlaboratory comparison of new materials for carbon and oxygen isotope ratio measurements. Reference and intercomparison materials for stable isotopes of light elements[END_REF] and an in-house dolomite standard. C and O isotopic measurements were performed at CRPG laboratory (Vandoeuvre-Lès-Nancy, France) and systematically doubled in order to check homogeneity of analyses. Isotopic compositions are reported in δ notation (in ‰) relative to V-PDB for carbon and converted to V-SMOW for oxygen. Standard deviations were 0.1‰ and 0.05‰ for δ 18 O and δ 13 C respectively.

Carbon and oxygen isotope compositions allow to determine the isotopic fractionation (∆ 13 C and ∆ 18 O) between two adjacent mineral phases (here gray and white dolomite). It is assumed that the isotopic equilibrium is reached when the two adjacent phases are directly bonded in similar structural positions (e.g., tension gashes). C and O isotope fractionation equations between dolomite and water in one hand [START_REF] Golyshev | Fractionation of stable oxygen and carbon isotopes in carbonate systems[END_REF][START_REF] Horita | Oxygen and carbon isotope fractionation in the system dolomite-water-CO2 to elevated temperatures[END_REF], and dolomite and dissolved CO2 in the other hand [START_REF] Golyshev | Fractionation of stable oxygen and carbon isotopes in carbonate systems[END_REF][START_REF] Ohmoto | Isotopes of sulfur and carbon[END_REF] were used to determine equilibrium temperature between the two adjacent carbonates. At isotopic equilibrium, this temperature corresponds to the precipitation temperature of the newly precipitated white dolomite [START_REF] Valley | Stable isotope thermometry at high temperatures[END_REF]. Only samples from Névache outcrop do not allow equilibrium temperature determination due to the lack of associated primary and secondary dolomites samples.

Results

Rock formations and deformations of the "Nappe des Gypses"

The "Nappe des Gypses" unit corresponds to an evaporitic sedimentary unit consisting of evaporites (mainly anhydrites), marine carbonates (dolomite), shales and clay layers (micaschists). Anhydrite rock is the most abundant rock formation of the unit. It forms massive outcrops all along the Arc valley in the Bramans, Sollières, Ambin outcrops and ruiniform landscape marked by numerous sinkholes in the Mont-Cenis lake area. Sedimentary bedding is underlined in white anhydrite by numerous pale gray layers that are commonly folded. Due to massive recrystallization, alpine foliations are most of the time invisible.

Marine carbonates have been observed in all sampling area. They consist of centimeter to meter scale bed of gray dolomite where for instance fine stromatolitic lamines are preserved (Fig. 3A-B). In the Bramans area, marine carbonates beds present automorph to sub-automorph albite (Fig. 3F and Supplementary Fig. A1) that forms crystals with X-Carlsbad and twins called "Roc Tourné" [START_REF] Rais | Distribution of authigenic albites in a limestone succession of the Helvetic Domain, eastern Switzerland[END_REF]. As these carbonate beds have been strongly folded and stretched by alpine tectonics, they form various scale blocks more or less rounded and elongated within the main anhydrite layering, called "boudins". A secondary white dolomite, in which no sedimentary structures are visible, is frequently observed around marine carbonate blocks. White dolomite probably formed by dissolution of a part of the gray dolomite beds and recrystallization.

Sedimentary clay layers are observed in each studied outcrop and present variable local thickness ranging from the plurimetric (Fig. 2C-D-F) to the microscopic scale (Fig. 4). They have been transformed into micaschist by alpine metamorphism where white micas are abundant and underline alpine foliations. Shales are scarcer as they have only been observed in Ambin stream and in Roubion stream near Névache village (Supplementary Fig. A2). They contain organic matter that leads to the formation of graphite during the metamorphism which is then used for the RSCM technique.

D1 deformations

At the regional scale, the D1 deformation is related to the nappes emplacement onto the Briançonnais unit (e.g., [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[END_REF][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF]. On the studied outcrops, such a dominant simple shear deformation induces a very penetrative deformation of anhydrite layers. Marine carbonate beds, micaschists and shales are folded, thinned and stretched. This produce spectacular isoclinal recumbent folds hinges (Fig. 2F), and "boudins" of the gray dolomitic beds (Fig. 2A,G). White dolomite crystals of mm to cm size form pressure shadows located between gray dolomite "boudins". This leads us to propose that these white dolomites crystallize during the stretching of the dolomite bed during D1 folding (Fig. 2E). Within clays layers, micas crystallize and underline a first foliation thereafter called D1 foliation. As D1 folding is isoclinal, bedding use to be transposed to the D1 foliation (Fig. 2A). This tectonic sequence is recognized in all of the studied outcrops (Fig. 2).

D2 foliation (S2)

Large scale D2 deformations are evidenced by the evolution of the S0-S1 structures from a site to another. For example, between the Sollières and the Ambin stream outcrops (Fig. 1), the main beddings dip to the SW and to the NE forming a syncline structure. From these outcrops, we estimate a N30-N45° striking fold axis postdating D1 structure that is compatible with D2 structures as observed elsewhere (e.g., [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[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][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][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF].

At the outcrop scale (Sollières, Névache and Mont-Cenis lake outcrops), D2 folds affect the D1 foliation and a very discrete D2 foliation underline the D2 fold axis planes (Fig. 2D). Within micaschist layers, D2 foliation, when present, is more penetrative and forms a crenulation foliation associated with a second generation of white micas (Fig. 4). At the Mont-Cenis, the main layering of the anhydrite rocks is here sub-horizontal (Fig. 2C). This layering cut some micaschists and dolomite blocks displaying a vertical foliation. Very similar relationships are observed in the Ambin Massif with the association of an early vertical D1 and a late D2 foliation [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[END_REF]. This suggests that D2 deformation is strongly expressed in the Mont Cenis area (Fig. 2C and4A).

Thin sections from the Bramans outcrop illustrate the relative relationships between white micas and albite: observations reveal that albite crystallize between two white micas crystallization phases (Fig. 4B). We propose, according to the surrounding area [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[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][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF], that albite crystallization, including D1 white micas, is surrounded by D2 white micas and then marks the very beginning of the D2 deformation phase (Fig 4).

D3 tectonics

In all sampling sites, late brittle structures are very scarce and late faults in evaporite rocks have not been clearly identified. In many places, dolomitic beds are affected by an intense fracturation without preferred orientation rather suggesting an episode of hydraulic fracturing postdating the D1 event (Fig. 3E). From available evidences, temporal relationship between this episode and the D2 structure cannot been established. These observations are in agreement with the global deformation history of the Vanoise cover during the Alpine subduction along the Briançonnais and Liguro-Piemont suture zone (e.g., [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].

Temperatures determination

C and O isotopes signatures and equilibrium temperatures in dolomite

Carbon and oxygen isotope compositions (δ 13 C and δ 18 O respectively) were determined on 50 dolomite samples collected in all the 6 studied outcrops. We collected micritic gray dolomites of sedimentary (primary) origin and white dolomite that crystallize during alpine events (Fig. 2G).

We distinguished during sampling a first white dolomite generation that is associated with D1 structures from a second generation, postdating D1 deformations that is observed all around gray dolomitic beds without preferential orientation (Fig. 2E-H).

Sedimentary "boudins" of gray dolomite present δ 18 O values around 23.75‰V-SMOW ± 1.49 and δ 13 C values around 3.50‰PDB ± 0.63. The D1 well-crystallized white dolomite presents δ 18 O values around 20.63‰V-SMOW ± 2.44 and δ 13 C values around 2.85‰PDB ± 1.04. In Sollières, Bramans, Ambin and Mont-Cenis outcrops, obtained isotopic equilibrium temperatures for the D1 white dolomites present the same range of temperature between 330°C and 432°C with a mean at 383°C ± 30°C, whereas temperatures at the Tignes outcrop seems to be shifted to higher temperature between 385°C and 475°C with a mean at 436°C ± 46°C (Table 1).

The later white dolomite surrounding the gray dolomitic "boudins" displays δ 18 O values around 16.95‰V-SMOW ± 1.52 and δ 13 C values around 1.94‰PDB ± 0.83. All oxygen and carbon stable isotopes values are summarized in Supplementary Table A1. Obtained isotopic equilibrium temperatures on samples from Sollières, Bramans and Ambin stream, ranged between 218°C and 289°C with a mean at 259°C ± 24°C (Table 1).

Fluid inclusions characterization

Our microthermometric data are considered with those from [START_REF] Barré | First evidence of the trisulfur radical ion S3and other sulfur polymers in natural fluid inclusions[END_REF] to bring a general view. Textural observations highlight that mineral phases suitable for fluid inclusions studies are albite, quartz and fluorite. Anhydrite always occurs as micro-crystallized phases indicating continued recrystallization during all the metamorphic events implying scarce and small fluid inclusions. On the contrary, quartz and fluorite appear well crystallized as isolated rounded crystals within the anhydrite (Fig. 3C-D) or within fracture cross cutting gray dolomitic beds. Such an association of quartz and fluorite is observed at Bramans and Sollières outcrops, when only quartz is observed in other studied outcrops. In every case, we never observed quartz and fluorite elongated within the S0-S1 foliation indicating that these minerals postdate the D1 deformation. Albites observed in the Bramans area (Fig. 3F and Supplementary Fig. A1), are of peculiar interest as they probably underline the beginning of the D2 deformation phase as evidenced by their crystallization between D1 and D2 foliation [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[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][START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF].

Fluid inclusions are found either isolated or in fracture planes whatever the host mineral (quartz, fluorite or albite). All the fluid inclusions are at least biphasic (aqueous and gaseous phase) and often contain one to three solid phases. All homogenize into the liquid phase, but they present different range of Th depending on the nature of the host mineral and the sampling sites (Fig. 5; Supplementary Table A2). Considering Th and fluid inclusions planes intersects in a single grain, it is possible to discriminate different generations (up to 4) of fluid inclusions (Fig. 6). There is a decrease of temperature from the first generation to the last implying episodic recording of the fluid flows during the "Nappe des Gypses" with a cooling from 207°C ± 3°C to 144°C ± 3°C.

Albite records the highest temperatures measured in fluid inclusions with Th ranging between 330 and 340°C with a mean value at 334°C ± 5°C. In quartz and fluorite, different generations of colder homogenization temperatures are recorded from a site to another: 1) A first generation of fluid inclusions hosted in quartz from Ambin show Th ranging between 293 and 319°C with a mean at 301°C ± 10°C. 2) A second generation hosted in quartz from the Mont-Cenis lake, show a range of Th from 230 to 287°C with a mean at 260°C ± 18°C, representing a regional temperature record. 3) A third phase of fluid inclusions trapping is recorded by a second generation of quartz from Ambin and a first generation of fluorite from Sollières l'Envers. The associated Th range between 203 and 244°C with a mean at 219°C ± 17°C at Ambin, and between 188 and 273°C with a mean at 214°C ± 27°C at Sollières l'Envers. 4) Finally, the last fluid inclusions trapping phase is recorded whatever the hosting mineral (quartz or fluorite) at Bramans, Ambin and Sollières l'Envers outcrops. At Bramans, Th range between 100 and 185°C with a mean at 137°C ± 28°C. They range between 138 and 179°C with a mean at 153°C ± 15°C at Ambin, and between 98 and 181°C with a mean at 147°C ± 26°C from Sollières outcrops.

The four outcrops studied here present a very similar range of Th between 50 and 350°C suggesting a very similar thermal evolution. Th of fluid inclusions hosted by the associated secondary mineral phases (quartz and fluorite) present similar range, confirming the coeval formation of these two minerals.

Graphite geothermometer

Graphite-rich metamorphic shales were collected in Ambin stream and in Roubion stream near Névache village (Supplementary Fig. A2). A total of 51 spectra was recorded from 2 samples from Ambin and 1 sample from Névache (Supplementary Fig. A3). A minimum of 15 spectra per sample was used to determine the temperatures. Table 2 shows the temperature mean values and the 1σ dispersion obtained from the different spectra per sample. The temperatures were deduced by Raman Spectroscopy of Carbonaceous Material (RSCM) with ranges from 397 to 438°C with a mean at 418°C ± 29°C for Ambin outcrop and at 359°C ± 8°C for Névache outcrop (Table 2).

Temperature ranges for D1, D2 and D3 event

We interpret the estimated isotopic equilibrium temperature between D1 white and gray sedimentary dolomites as the range of temperature of the D1 event: between 330°C and 432°C in Sollières, Bramans, Ambin and Mont-Cenis outcrops and between 385 and 475°C in the Tignes outcrop (Table 1).

The transition between D1 and D2 deformation phase is constrained in one hand by the lowest D1 isotopic equilibrium temperatures at 356°C ± 21°C (mean of the minimum D1 temperatures per site) and on the other hand by the 334°C ± 5°C Th on albite from fluid inclusions of the Bramans outcrop. Regarding the proximity and continuity of the outcrops of the Arc Valley (Bramans, Sollières, Ambin) and that of the Mont-Cenis outcrop, we consider similar conditions between these 4 outcrops. This leads us to propose a 350°C ± 20°C temperature for the D1-D2 transition on the Sollières, Bramans, Ambin and Mont-Cenis area.

Temperatures deduced on post-D1 deformation minerals such as quartz and fluorite studied by fluid inclusions analyses and stable isotopes signatures on white dolomites surrounding the dolomitic "boudins", display temperatures below the D1-D2 transition temperature. This suggests that D2 and later temperatures are never higher than D1 temperature in the "Nappe des Gypses". From a textural and thermal point of view, early fluid inclusions present the higher Th than the subsequent generations, as expected during a thermal retrograde path.

The maximum temperature of 418°C ± 29°C for Ambin outcrop provided by RSCM is consistent with the maximum 431°C ± 28°C (mean of the maximum D1 temperatures per site) provided by isotopic equilibrium temperature on D1 structure.

Temperature results suggest that maximum metamorphic temperature in the "Nappe des Gypses" is reached during D1 event at 431°C ± 28°C for the Ambin, Sollières, Bramans, Mont-Cenis and Tignes outcrops. Assuming that heating during the early stage of exhumation is never observed in the internal part of western alps, we propose that the maximum temperature is representative of the "Nappe des Gypses" temperature at its maximum of burial. In addition, in absence of heating evidence, we propose that temperatures decrease progressively during exhumation, from 431°C ± 28°C to 356°C ± 21°C during D1 event and then again at lower temperature during events D2 and later, as observed in fluid inclusions (Fig. 6).

Pressure estimates

White micas are mostly present within micaschist layers of the "Nappe des Gypses". Samples from the Bramans, Ambin stream and Mont-Cenis Lake outcrops contained sufficiently well preserved micaschists layers to allow Electron Microprobe measurement on large enough white micas (Fig. 4). The chemical composition of white micas from these outcrops follows the range of the muscovite-celadonite solid solutions that define phengites (Fig 7B .). Si content of phengite that is related to the crystallization pressure [START_REF] Massonne | Stability field of the high-pressure assemblage, talc+ phengite and two new phengite barometers[END_REF][START_REF] Vidal | P-T-deformation-Fe 3+ /Fe 2+ 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][START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: thermodynamic modelling and implications for thermobarometric estimates[END_REF] range between 3.54 and 3.27 Si a.p.f.u., 3.44 and 3.16 Si a.p.f.u. and 3.36 and 3.18 Si a.p.f.u. respectively in Bramans, Ambin stream and Mont-Cenis Lake outcrops (Supplementary Table A3 shows typical EMP chemical compositions of phengites used for pressure calculations). As a whole, the wide composition range suggest that a large range of pressure is recorded on the three outcrops. Highest pressures are expected on the Bramans and Ambin outcrop where phengite are the most Si-rich whereas phengites of the Mont-Cenis Lake outcrop probably only record pressures of D2 deformation according to field observations. In absence of textural relationships between chlorite and phengite as expected for thermobarometric data [START_REF] Vidal | P-T-deformation-Fe 3+ /Fe 2+ 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][START_REF] Dubacq | Dehydration of dioctahedral aluminous phyllosilicates: thermodynamic modelling and implications for thermobarometric estimates[END_REF], continuous estimate of the P-T conditions is not possible on the "Nappe des Gypses". We only provide estimation for the peak condition and the D1-D2 transition considering, in each case, D1 and D2 phengites that give the highest pressures.

According to the field observations and Electron Microprobe measurements, we selected phengites from the Bramans and Ambin outcrops to estimate the maximum conditions at a D1 maximum temperature of 430°C ± 7°C (Table 1) and 418°C ± 29°C, respectively. Obtained pressure range between 15.0 and 18.2 kbar at the Ambin and Bramans outcrops (Fig. 8A-B). We deduce a 5.2 to 7.8 kbars D1-D2 pressure using the highest-pressure results on the Mont-Cenis Lake outcrop at 356°C ± 21°C D1-D2 transition temperature (Fig. 8C). These data allow us to determine a P-T path for the "Nappe des Gypses" formation (Fig. 8D). These values are based on mean temperatures determined from fluid inclusions, Raman spectroscopy and oxygen isotope fractionation calculations. It implies a maximum overestimation of 2 kbars for the low temperature phengites.

Discussion

The maximum P-T condition of the "Nappe des Gypses".

Our results from the Sollières, Bramans, Ambin stream and Mont-Cenis Lake outcrops are consistent with metamorphic peak conditions with maximum temperature at 431°C ± 28°C and pressure at 16.6 ± 2.3 kbars. Note that the presence of jadeite in the "Nappe des Gypses" formation (e.g., [START_REF] Saliot | La jadéite dans les Alpes françaises[END_REF][START_REF] Grappin | Les variations des rapports Cl/Br, Na/Br et K/Br dans les inclusions fluides des quartz de la cicatrice évaporitique de Bramans-Termignon (Vanoise, Alpes françaises)[END_REF], is very consistent with these P-T conditions.

The albite-jadeite-quartz line at 431°C is at 11.8 ± 0.5 kbars [START_REF] Holland | The reaction albite= jadeite+ quartz determined experimentally in the range 600-1200 C[END_REF], corresponding to a minimum pressure in the "Nappe des Gypses". From the North-East (Tignes outcrop) to the South-West (Névache outcrop), RSCM temperatures and these obtained on ∆ 13 C and ∆ 18 O fractionations in dolomite show that the maximum temperatures increase from 359°C ± 7°C to 475 ± 7°C. As no micaschist layer suitable for pressure estimate was found in the Névache and Tignes outcrops, peak pressure conditions are not constrained. Whatever the pressure, maximum temperatures related to D1 tectonics recorded by the "Nappe des Gypses" increase from the South-West (Névache) to the North-East (Tignes). These temperatures fit perfectly into the North-South and East-West gradients already described in the literature [START_REF] Beyssac | Raman spectra of carbonaceous material in metasediments: a new geothermometer[END_REF][START_REF] Gabalda | Thermal structure of a fossil subduction wedge in the Western Alps[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] Plunder | How continuous and precise is the record of P-T paths? Insights from combined thermobarometry and thermodynamic modelling into subduction dynamics (Schistes Lustrés, W. Alps)[END_REF]. Here we highlight for the first time a variation of the P-T conditions during the metamorphic peak within the "Nappe des Gypses".

Metamorphic peak conditions on the Sollières, Bramans, Ambin stream and Mont-Cenis Lake outcrops are in very good agreement (Fig. 9) with maximum conditions reached in the overlying Liguro-Piemontese units (units III and IV after [START_REF] Agard | Exhumation of the Schistes Lustrés complex: in situ laser probe 40Ar/39Ar constraints and implications for the Western Alps[END_REF] and Median Units after [START_REF] Plunder | How continuous and precise is the record of P-T paths? Insights from combined thermobarometry and thermodynamic modelling into subduction dynamics (Schistes Lustrés, W. Alps)[END_REF]. This is also the case with the ∆ 13 C and ∆ 18 O fractionation temperatures obtained on the Tignes area as the maximum 475°C is similar with temperatures ranging between 450°C and 475°C in the Liguro-Piemontese units (Median Units after [START_REF] Plunder | How continuous and precise is the record of P-T paths? Insights from combined thermobarometry and thermodynamic modelling into subduction dynamics (Schistes Lustrés, W. Alps)[END_REF]. Such a good agreement regarding maximum P-T conditions is not found when comparing the results on the "Nappe des Gypses" and the covered Briançonnais units. In the Ambin Massif and Southern Vanoise unit, even if pressures are the same (14-16 kbars), temperatures are always slightly higher (500°C; [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[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]. In the Modane Aussois unit, nor the pressure (1.1 kbars), nor the temperature (350°C; [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] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF] fit with those of the "Nappe des Gypses" and a 5 kbars pressure gap implies very different maximum burial depth. In absence of pressure constraints, comparison with Houiller unit is only based on temperature that decrease from 359°C ± 7°C in the "Nappe des Gypses" to 300°C in the Houiller unit [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].

Assuming these similarities and differences, we propose that the "Nappe des Gypses" recorded HP metamorphism in very similar conditions than the Liguro-Piemontese median units.

As Triassic evaporitic formations are older than the Liguro-Piemontese oceanic domain, we propose that the "Nappe des Gypses" has been buried below Liguro-Piemontese units at high pressure conditions during Alpine continental subduction. This implies an initial distal position of evaporitic formations along the European passive margin as observed along current passive margins: Gulf of Mexico (e.g., [START_REF] Jackson | Regional extension as a geologic trigger for diapirism[END_REF], western Mediterranean Sea (e.g., [START_REF] Lofi | Seismic atlas of the Messinian Salinity Crisis markers in the Mediterranean and Black Seas[END_REF] and the Atlantic margins (e.g., [START_REF] Demercian | Style and pattern of salt diapirs due to thinskinned gravitational gliding, Campos and Santos basins, offshore Brazil[END_REF]. Such a burial of the "Nappe des Gypses" may have occurred during the early stage of exhumation of the Liguro-Piemontese units as they reached the same condition than the "Nappe des Gypses". Time constrains suggest a 62-55 Ma to 51-45 Ma range [START_REF] Agard | Exhumation of the Schistes Lustrés complex: in situ laser probe 40Ar/39Ar constraints and implications for the Western Alps[END_REF] for the burial of "Nappe des Gypses" below Liguro-Piemontese units.

The differences between peak conditions of the "Nappe des Gypses" and the Briançonnais units confirm the thrusting of the Triassic evaporitic formations onto the European basements and covers that have recorded various metamorphic conditions. As the nappe stack is folded during D2 event, we proposed that nappe emplacement occurred during the D1 geological event.

Exhumation of the "Nappe des Gypses"

Exhumation path of the "Nappe des Gypses" occurred in three stages (D1, D2, D3). Our results provide P-T constrains on metamorphic peak (431°C ± 28°C and 16.6 ± 2.3 kbars) and D1-D2 transition (350°C ± 20°C and 6.5 ± 1.8 kbars) conditions. For the D1-D2 transition, as in the case of the metamorphic peak, a good pressure agreement is observed with the Ambin units, whereas temperatures remain higher (500°C; [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[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]. This confirms the peculiar thermal evolution of the Ambin massif where a possible heating is described during D2 event [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]. P-T conditions for the D1-D2 transition are consistent with conditions recorded within the Modane-Aussois unit (350°C; [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] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF]. This confirms our structural observations indicating that the whole nappe stack is folded during D2 event. Such an event has been dated between 43 Ma and 35 Ma [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] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF].

In the surrounding units, D3 phase is characterized by brittle structure [START_REF] Ganne | Structural and metamorphic evolution of the Ambin massif (western Alps): toward a new alternative exhumation model for the Brianconnais domain[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][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]. Fluid inclusions hosted in quartz from carbonate-quartzsulfide veins provide a maximum 260°C ± 18°C homogenization temperatures (Fig. 5C; Supplementary Table A2) that may represent the D2-D3 transition temperature. This is consistent with the 259°C ± 24°C temperatures obtained on ∆ 13 C and ∆ 18 O fractionations on the second type of white-dolomite post-dating the D1 crystallization on the same site (Table 1). Thus, our data are consistent with a temperature of 259°C ± 24°C for the D2-D3 transition. In absence of a more precise structural control on phengite, it is however not possible to deduce a D2-D3 transition pressure based on our own results. Considering this, we propose that D2-D3 transition within the "Nappe des Gypses" and the Modane-Aussois unit are the same at 2.0 ± 1.0 kbars occurred at circa 30 Ma [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] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF].

Our new P-T-d results provide a new highlight on the exhumation of the Liguro-Piemontese units as we demonstrate that its tectonic contact with the "Nappe des Gypses" occurred at pressure as high as 16.6 ± 2.3 kbars. About 10 kbars of exhumation of the Liguro-Piemontese units is related to the D1 phase of the "Nappe des Gypses" in association with NW to W thrusting direction over the Briançonnais unit. At the same time, the internal deformation of the Liguro-Piemontese units is characterized by top to the east direction of shearing in a context of extension (e.g., [START_REF] Tricart | From passive margin to continental collision; a tectonic scenario for the Western Alps[END_REF][START_REF] Agard | Exhumation of the Schistes Lustrés complex: in situ laser probe 40Ar/39Ar constraints and implications for the Western Alps[END_REF]. These opposite senses of shear that are active during exhumation, are in good agreement with extrusion tectonics model of exhumation (e.g., [START_REF] Godin | Channel flow, ductile extrusion and exhumation in continental collision zones: an introduction[END_REF]. Such a scenario implies that the "Nappe des Gypses" unit acted as an efficient "décollement" layer especially as it reaches its metamorphic peak. Figure 10 presents a possible sketch of the kinematic of the western Alps with the integration of the "Nappe des Gypses" tectonic evolution according to our results.

Phase transitions in the "Nappe des Gypses" and implications for exhumation processes

Evaporite layer such as the "Nappe des Gypses" unit are considered as some of the most efficient "décollement" layers within fold and thrust belts [START_REF] Davis | The role of salt in fold-and-thrust belts[END_REF]. Anhydrite which is the most abundant component of the "Nappe des Gypses" unit, forms together with other salts a very soft material regarding to the rheological properties of crustal rocks as it reaches his brittle-ductile transition at around 5 km depth assuming a normal geothermal gradient [START_REF] Dahlen | Critical taper model of fold-and-thrust belts and accretionary wedges[END_REF]. This is consistent with the intense plastic deformation of the "Nappe des Gypses" unit (Fig. 2) and suggests that it can act as an efficient "décollement" layer during the whole burial exhumation path.

Transition phase of gypsum into bassanite and then bassanite into anhydrite release an amount of water estimated at 37% of the initial volume of gypsum [START_REF] Ko | Generation and maintenance of pore pressure excess in a dehydrating system 1. Experimental and microstructural observations[END_REF]. Recent experimental investigations on gypsum, bassanite and anhydrite [START_REF] Yang | Pressure-induced structural phase transition and dehydration for gypsum investigated by Raman spectroscopy and electrical conductivity[END_REF] better constrain the P-T conditions of these phase transitions (Fig. 9). Dehydration reactions are known to control the slip stabilization of the shearing and to generate pore pressure excess (e.g., [START_REF] Leclère | Permeability control on transient slip weakening during gypsum dehydration: Implications for earthquakes in subduction zones[END_REF].

Gypsum dehydration into bassanite is particularly efficient to generate pore overpressure [START_REF] Ko | Generation and maintenance of pore pressure excess in a dehydrating system 1. Experimental and microstructural observations[END_REF]. During the prograde path, dehydration reactions may have first transform gypsum into bassanite at ca. 9 kbar and 200°C, then bassanite into anhydrite at ca. 15 kbar and 300°C (Fig. 9; [START_REF] Yang | Pressure-induced structural phase transition and dehydration for gypsum investigated by Raman spectroscopy and electrical conductivity[END_REF]. This suggests that water is released from anhydrite near the end of the exhumation path.

In addition, [START_REF] Grappin | Les variations des rapports Cl/Br, Na/Br et K/Br dans les inclusions fluides des quartz de la cicatrice évaporitique de Bramans-Termignon (Vanoise, Alpes françaises)[END_REF] show that the water traps in the "Nappe des Gypses" only come from gypsum dehydration and Triassic connate seawater. Coupled to the pore overpressure generated by gypsum dehydration, it implies that the fluids evolved in a closed system all along the prograde and retrograde paths. Such a fluid overpressure at the metamorphic peak can generate mechanical decoupling and softening effects [START_REF] Fossen | Structural Geology[END_REF]. We propose that water releases from gypsum to anhydrite transition phases between 9 and 15 kbar and 200 and 300°C, is responsible for the increase of the "décollement" effect of the "Nappe des Gypses" and may have enhanced its exhumation. Such an overpressure event is evidenced, on the field, by hydraulic breccias and hydrothermal veins observed within the dolomite layers (Fig. 3E).

The transition from oceanic subduction to continental subduction between 50 and 40 Ma (HP ophiolites are dated at 45 Ma; [START_REF] Rubatto | Zircon formation during fluid circulation in eclogites (Monviso, Western Alps): implications for Zr and Hf budget in subduction zones[END_REF] most likely marks the implication of the soft "Nappe des Gypses" in the subduction/exhumation dynamics. Consequently, the exhumation of HP units initially controlled by the internal dynamic of the accretionary wedge [START_REF] Agard | Exhumation of the Schistes Lustrés complex: in situ laser probe 40Ar/39Ar constraints and implications for the Western Alps[END_REF][START_REF] Guillot | Exhumation processes in oceanic and continental subduction contexts: a review[END_REF] evolved to a thrust and fold belt system where the activation of a décollement layer within Triassic evaporites (Fig. 10) related to dehydration processes and fluid over-pressure play a major role, which was not recognized so far.

Conclusions

This study is the first to couple a complete petrographic, geochemical and tectonic investigation from the micro-to the subduction-collision scales of an evaporitic formation. The P-T-d evolution of the "Nappe des Gypses" formation within the Alpine tectonic framework has been determined for the first time. Its main conclusions are:

1. Maximum P-T conditions reached by the "Nappe des Gypses" are 16.6 ± 2.3 kbars and 431°C ± 28°C. These are similar to those reached by the overlying median Liguro-Piemontese units [START_REF] Plunder | How continuous and precise is the record of P-T paths? Insights from combined thermobarometry and thermodynamic modelling into subduction dynamics (Schistes Lustrés, W. Alps)[END_REF] and support 1) an HP tectonic contact between these units and 2) an active role of the "Nappe des Gypses" during the exhumation of the Alpine oceanic complex.

2. Differences in maximum P-T conditions reached by the "Nappe des Gypses" and the Briançonnais units (Ambin, Modane Aussois and Houiller units) confirm allochthonous origin of the "Nappe des Gypses" that is probably originating from a more distal part of the European Tethyan passive margin.

3. Fluid overpressure generated by sulfates dehydration at the end of the prograde path has probably later enhanced exhumation processes at 15.0 to 16.6 kbars. This integrative study shows how complementary use of petrography, crystallochemistry, fluid inclusions, stable isotopes (C, O) analysis and structural geology allows to constrain the P-T-d path of an evaporitic unit and precise its structural role. This methodology could be applied to other contexts where evaporites play a key role in the geodynamic structuration of orogens (e.g., Pyrénées, Himalaya, Zagros). [START_REF] Beyssac | Raman spectra of carbonaceous material in metasediments: a new geothermometer[END_REF] on carbonaceous material from black shales layers incorporated in the "Nappe des Gypses" formation. A2. The albite = jadeite + quartz line is also represented [START_REF] Holland | The reaction albite= jadeite+ quartz determined experimentally in the range 600-1200 C[END_REF], with ab = albite; jd = jadeite; qtz = quartz. HZ = 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]; MA = Modane-Aussois unit [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF]; LPU = Liguro-Piemontese Upper unit [START_REF] Agard | Exhumation of the Schistes Lustrés complex: in situ laser probe 40Ar/39Ar constraints and implications for the Western Alps[END_REF]; LPM = Liguro-Piemontese Median unit [START_REF] Plunder | How continuous and precise is the record of P-T paths? Insights from combined thermobarometry and thermodynamic modelling into subduction dynamics (Schistes Lustrés, W. Alps)[END_REF]; LPL = Liguro-Piemontese Lower unit [START_REF] Plunder | How continuous and precise is the record of P-T paths? Insights from combined thermobarometry and thermodynamic modelling into subduction dynamics (Schistes Lustrés, W. Alps)[END_REF]; SV = South-Vanoise unit [START_REF] Strzerzynski | Tectono-metamorphic evolution of the Briançonnais zone (Modane-Aussois and Southern Vanoise units, Lyon Turin transect, Western Alps)[END_REF]; GP = Gran Paradiso massif [START_REF] Gabudianu Radulescu | The age of HP metamorphism in the Gran Paradiso Massif, Western Alps: a petrological and geochronological study of "silvery micaschists[END_REF]; MV = Monviso massif [START_REF] Angiboust | Eclogitization of the Monviso ophiolite (W. Alps) and implications on subduction dynamics[END_REF].

FIGURE CAPTIONS

  synthetic and natural fluid inclusions with the following phase transitions: the ice melting of pure H2O at 0.0°C and the critical homogenization temperatures (Th) of natural fluid inclusions at 374.0°C. The accuracy for total homogenization temperatures is about ±1°C, according to the calibration curves. Total homogenization temperatures were measured in inclusions hosted in quartz, albite and fluorite. Due to high deformation of the crystal and change of the internal volume during heating, Th measurements from fluid inclusions hosted in anhydrite and fluorite are given for indicative purpose and for comparison with Th obtained in syngenetic quartz-hosted fluid inclusions.

Fig. 1 .Fig. 2 .

 12 Fig. 1. (A) Simplified geological map of the studied area (Western French Alps; modified from

Fig. 3 .

 3 Fig. 3. Photographs of typical paragenesis observed in the "Nappe des Gypses" formation. (A

Fig. 4 .

 4 Fig. 4. Microphotographs of thin sections at crossed polarized light showing typical

Fig. 5 .

 5 Fig. 5. Histograms of homogenization temperatures (Th) obtained by microthermometry on

Fig. 6 .

 6 Fig. 6. A. Microphotograph of typical fluid inclusions planes hosted in fluorite from Sollières

Fig. 7 .

 7 Fig. 7. (A) Si vs Interlayer cation diagram for analyzed phengites. (B) K-white mica chemistry

Fig. 8 .

 8 Fig. 8. (A), (B) and (C) represent the pressure determination for the Bramans, Ambin and Mont-

Fig. 9 .

 9 Fig. 9. Comparison between "Nappe des Gypses" unit retrograde P-T path and these from the

Fig. 10 .

 10 Fig. 10. Possible kinematic reconstruction of the Western Alps with highlights on the "Nappe

TABLE CAPTIONS Table 1 .

 CAPTIONS1 Equilibrium temperatures determined from dolomite-H2O fractionation equations from[START_REF] Golyshev | Fractionation of stable oxygen and carbon isotopes in carbonate systems[END_REF] and[START_REF] Horita | Oxygen and carbon isotope fractionation in the system dolomite-water-CO2 to elevated temperatures[END_REF] of carbon and oxygen isotopes values.

Table 2 .

 2 Maximum temperatures obtained by the geothermometer described by

Acknowledgements

This work was funded by LABEX ANR-10-LABX-21-01 Ressources21 (Strategic metal resources of the 21 st Century) and the French Ministry of Higher Education and Research. The authors are extremely grateful to T. Rigaudier for the C and O stable isotope analysis performed at Centre de Recherches Pétrographique et Géochimiques (CRPG) of Nancy, and to Christophe Scheffer for the pressure estimations calculations. The paper benefited greatly from the review provided by A. Plunder and one anonymous reviewer. We thank P. Agard for his very relevant comments that have improved this paper and for his editorial support.

Appendix A. Supplementary material

Supplementary material related to this article can be found as attached file.