# Discovery of Variscan orogenic peridotites in the Pelvoux Massif (Western Alps, France)

Jean-Baptiste Jacob<fnref n="1" />,<fnref n="2" />,* and Emilie Janots<fnref n="1" />, Carole Cordier<fnref n="1" /> ID and Stéphane Guillot<fnref n="1" />

Received: 17 August 2022 / Accepted: 14 December 2022 / Publishing online: 8 March 2023

**Abstract** – Small bodies of mantle-derived peridotites and other ultramafic rocks are commonly found in exhumed lower crustal units of collisional orogens. They provide a direct record of the complex evolution of the upper mantle before and during an orogeny, and are therefore key markers of the geodynamic evolution of an orogen. We report here the discovery of such mantle-derived peridotites, which occur as fragmented enclaves in migmatites of the high-grade Variscan lower crust exposed in the Pelvoux Massif (external Western Alps). A wide petrographic diversity has been observed, from very fertile, garnet-bearing lherzolites, to more depleted spinel/chromite-bearing harzburgites. Thermobarometric calculations on a garnet lherzolite indicate an initial stage at 3.0–4.0 GPa and 970–1140 °C, followed by exhumation to 0.8–1.3 GPa and 800–850 °C, while the harzburgites do not show any evidence of equilibration in the garnet field. Petrological observations, whole-rock geochemistry and *in situ* mineral compositions suggest the peridotites have undergone a complex history prior to their incorporation in the lower crust during the Variscan Orogeny. They derive from a refractory mantle, which has experienced variable degrees of melt depletion, and has then been extensively refertilized. Cryptic metasomatism is observed in all samples. It is characterized by an enrichment in large-ion lithophile elements (LILE, in particular Cs, Rb, U and Pb) relative to high field strength elements (HFSE), in particular Nb and Ta. This cryptic metasomatism is presumably related to percolation of subduction-related fluids or melts in the mantle. In addition, modal metasomatism occurred in some samples, where crystallization of phlogopite, pargasite, chromite and apatite has been observed. This modal metasomatism resulted in significant enrichment in K₂O, P₂O₅ and Cr₂O₃ of the bulk rock, together with a strong enrichment in incompatible LREE relative to HREE. These geochemical characteristics are strikingly similar to that of syn-collisional, Mg–Cr–LILE rich mantle-derived (ultra)-potassic magmas such as durbachites and vaugnerites, which are ubiquitous in the Variscan metamorphic allochthons of Massif Central, external Alps, Vosges and Bohemian Massif. We therefore suggest that this metasomatism results from dynamic percolation of the peridotites by K₂O–P₂O₅–Cr₂O₃-rich melts from which the durbachites and vaugnerites are primarily derived. These geochemical characteristics are in line with whole-rock Nd isotopic compositions, which indicate enrichment of the mantle by a continental crust component, presumably related to Variscan subductions. This evolution is consistent with that of other Variscan peridotites in the Eastern Alps (Ulten) and the Bohemian Massif, where multiple metasomatic episodes related to melts or fluids released in Variscan subduction zones have been documented.

**Keywords:** orogenic peridotites / cryptic metasomatism / modal metasomatism / mantle refertilization / melt-peridotite interaction / Variscan Belt

**Résumé** – Découverte de péridotites orogéniques varisques dans le massif du Pelvoux (Alpes occidentales, France). Les unités de croûte inférieure exhumées au sein des orogènes collisionnels contiennent fréquemment des petits corps de péridotite et d'autres roches ultramafiques. Ces roches sont des marqueurs clés de l'évolution géodynamique d'un orogène, car elles fournissent un enregistrement direct de l'évolution du manteau avant et au cours des phases orogéniques. Nous rapportons ici la découverte de péridotites orogéniques au sein d'unités de croûte inférieure varisque, dans le massif du Pelvoux (Alpes occidentales externes). Celles-ci se présentent sous forme d'enclaves fragmentées, emballées dans des

<footnote n="1">¹ Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, IRD, Univ. Gustave Eifel, ISTerre, 38000 Grenoble, France</footnote>

<footnote n="2">² The Njord Centre, Department of Geoscience, University of Oslo, Oslo, Norway</footnote>

<footnote>*Corresponding author: j.b.jacob@mn.uio.no</footnote>
migmatites. On observe une grande diversité pétrographique, allant de lherzolites à grenat très fertile à des harzburgites à spinelle/chromite plus appauvries. Les estimations thermobarométriques effectués sur une lherzolite à grenat indiquent un stade initial à 3,0-4,0 GPa et 970-1140 °C, suivi d'une exhumation à 0,8-1,3 GPa et 800-850 °C, tandis que les harzburgites ne montrent aucune trace d'équilibration dans le champ du grenat. Les observations pétrographiques, les données géochimiques sur roche totale et les compositions minérales *in situ* indiquent que les péridotites ont enregistré une histoire complexe, précédant leur incorporation dans la croûte inférieure au cours de l'orogenèse Varisque. Elles dérivent de domaines mantelliques réfractaires, qui ont subi différents degrés de fusion partielle avant d'être refertilisés. Un métasomatisme cryptique est observé dans tous les échantillons. Il est caractérisé par un enrichissement en éléments lithophiles (en particulier en Cs, Rb, U, Pb) par rapport aux HFSE (High Field Strength Elements), en particulier Nb et Ta. Ce métasomatisme est vraisemblablement lié à la percolation de fluides / liquides silicatés issu d'un zone de subduction. Par ailleurs, on observe dans certains échantillons un métasomatisme modal, caractérisé par la cristallisation de phlogopite, de pargasite, de chromite et d'apatite ainsi que par un fort enrichissement en LREE (Light Rare Earth Elements) incompatibles par rapport aux HREE (Heavy Rare Earth Elements). Ces caractéristiques géochimiques sont très similaires à celles des durbachites et des vaugnérites, des séries de roches magmatiques mantelliques (ultra)-potassiques riches en Mg-Cr-LILE, mises en place au Carbonifère moyen à supérieur dans les unités métamorphiques allochthones du massif central, des Vosges, du massif de Bohème, ainsi que dans le socle varisque des Alpes externes. Ceci suggère l'existence d'un lien génétique entre les durbachites et vaugnérites et le métasomatisme modal observé dans les enclaves. Ce dernier pourrait ainsi être lié à la percolation dans le manteau des magmas primaires riches en K₂O-P₂O₅-Cr₂O₃ dont sont issues les durbachites. Ces caractéristiques géochimiques sont en accord avec les compositions isotopiques en Nd, qui indiquent un enrichissement du manteau par du matériel issu de la croûte continentale, vraisemblablement lié aux subductions varisques. Cette évolution est cohérente avec celle d'autres péridotites varisques dans les Alpes orientales (Ulten) et dans le massif de Bohème, où des épisodes métasomatiques liés à la percolation des fluides / liquides magmatiques mantelliques issus des zones de subduction varisques est bien documentée.

**Mots clés :** péridotites orogéniques / métasomatisme cryptique / métasomatisme modal / refertilisation mantellique / interactions magma-péridotite / chaîne Varisque

# 1 Introduction

Occurrence of mantle-derived ultramafic slices, lenses or boudins is a common feature of high-grade lower crustal units exposed in collisional orogens worldwide (e.g. Brueckner and Medaris, 2000; Nimis and Trommsdorff, 2001; Medaris et al., 2005; Chen et al., 2015; Brueckner, 2018). Although they only represent a tiny fraction of these high-grade units, ultramafic bodies provide extremely valuable information about metasomatic and magmatic processes occurring in the upper mantle during and before an orogeny (Malaspina et al., 2006; Scambelluri et al., 2006; Sapienza et al., 2009; Brueckner et al., 2010; Kubeš et al., 2022). They may also retain pressure-temperature information and thus provide constraints about burial and exhumation of (ultra)-high-pressure (U-HP) units in collision zones. They are therefore crucial markers to constrain the geodynamic evolution of an orogen.

In the Variscan Belt (Fig. 1a), small ultramafic mantle bodies, including garnet peridotites, spinel peridotites and pyroxenites, are found in the metamorphic allochthonous units that form exhumed portions of the lower-to-mid-crustal orogenic root, mainly in the Bohemian Massif (Medaris et al., 2005, 2015; Kubeš et al., 2022), the Vosges-Black Forest Massif (Altherr and Kalt, 1996; Altherr, 2021) and the French Massif Central (Gardien et al., 1990; Godard, 1990). In detail, these mantle bodies are found in exhumed lower to mid-crustal units from different tectonic domains, which record different metamorphic evolutions (Medaris et al., 2005; Lardeaux et al., 2014; Altherr and Soder, 2018). In the

Bohemian Massif, where the largest number of these ultramafic bodies has been found, a distinction has been recognized between peridotite-bearing units exhumed close to the plate interface (Saxo-Thuringian suture) and those exhumed further south, between the former Devonian subduction arc and the Brunia continental backstop (Moldanubian zone, Schulmann et al., 2009). The former typically record colder peak-P conditions (< 800 °C) and shorter-lived peak-T than the latter, which record peak-P at UHT conditions (850-1050 °C), followed by near-isothermal decompression (Maierová et al., 2021, and references therein). Therefore, these different units correspond to contrasted exhumation mechanisms, and mark different processes within the Variscan collision. While the plate-interface probably correspond to a former subduction *mélange* exhumed in the subduction channel, exhumation of Moldanubian peridotites and other UHP rocks might involve less conventional mechanisms of lower crustal “relimination” (Schulmann, Lexa et al., 2014; Maierová et al., 2018, 2021).

Multiple studies of the mantle-derived fragments have been undertaken in the Bohemian Massif (e.g. Medaris et al., 2005, 2015; Ackerman et al., 2020; Kubeš et al., 2022), which may therefore be a good reference site to discuss these processes. In other Variscan Massifs, deep crustal and mantle processes have been comparatively less studied, mainly because of the scarcity of well-preserved mantle bodies. This is particularly the case in the southeastern Variscides, where there is virtually no study on mantle fragments, apart from the garnet peridotites of the Ulten Zone, in the austroalpine
<figure><img src="image_1.png" /></figure>

**Fig. 1.** (a) Location of the External Crystalline Massifs of the Alps in the Variscan Belt. The Mid-Variscan allochthons, which contain portions of the exhumed crustal root, are highlighted in orange. Modified after Ballèvre *et al.* (2014) and Martínez Catalán *et al.* (2021). BF: Black Forest; ECM: External Crystalline Massifs; FMC: French Massif Central; GTMZ: Galicia-Trás-os-Montes; MT: Maures-Tanneron; RHS: Rheinische schiefergebirge; SPZ: South Portugese Zone; ST: Saxo-Thuringian suture; Vo: Vosges. (b) Simplified geological map of the central part of the ECM (Belledonne, Pelvoux, Mont-Blanc and Aiguilles-Rouges), modified from Guillot and Ménot (2009). The investigated area (black rectangle) is located in the central part of Pelvoux Massif (Fig. 2), but two samples of serpentinites were collected in the southern part of Belledonne.

**Fig. 1.** (a) Localisation des Massifs Crystallins Externes dans la chaîne Varisque. Les unités métamorphiques allochthones (« Mid-Variscan allochthons »), dans lesquels se trouvent les unités de croûte inférieure, sont figurés en orange. Modifié d'après Ballèvre *et al.* (2014) et Martínez Catalán *et al.* (2021). BF: Forêt-Noire; ECM: Massifs Cristallins Externes; FMC: Massif Central; GTMZ: Galice-Trás-os-Montes; MT: Maures-Tanneron; RHS: Rheinische schiefergebirge; SPZ: Zone sud-portugaise; ST: Suture Saxo-Thuringienne; Vo: Vosges. (b) Carte géologique simplifiée de la partie centrale des MCE (Belledonne, Pelvoux, Mont-Blanc et Aiguilles-Rouges), modifiée d'après Guillot et Ménot (2009). La zone étudiée (rectangle noir) est située dans la partie centrale du massif de Pelvoux (Fig. 2), mais deux échantillons de serpentinites ont été collectés dans la partie sud de Belledonne.

basement (Godard *et al.*, 1996; Tumiati *et al.*, 2003; Scambelluri *et al.*, 2006). We document here the discovery of new outcrops of peridotites in the Variscan basement of the external Western Alps, in exhumed lower crustal units of the Pelvoux massif (Fig. 1b). These peridotites occur as fragmented enclaves in migmatites, and some of them have been moderately affected by serpentinization.

This contribution aims to provide field, petrological, thermo-barometric, and whole-rock geochemical and isotopic data of these enclaves, in order to discuss their metamorphic and metasomatic evolution. This evolution will then be compared with that of other Variscan mantle bodies, in particular from the Bohemian Massif, in order to provide some insights into the evolution of the Variscan mantle in the external Western Alps.
# 2 Geological setting

## 2.1 General setting

The Pelvoux Massif is a well-exposed portion of Variscan middle-to-lower crust located in the French Western Alps (Fig. 1b). It belongs to the External Crystalline Massifs (ECM), which represent the exposed Paleozoic basement in the external Western Alps. The Alpine overprint in these massifs is generally considered to be very mild, and consists mostly in localized brittle/ductile deformation, which occurred at greenschist-facies conditions (<350 °C, Bellahsen *et al.*, 2014; Bellanger *et al.*, 2015). The Paleozoic basement of the ECM has been strongly overprinted during the Variscan orogeny (Guillot and Ménot, 2009; Jacob *et al.*, 2022), and displays similarities with the metamorphic allochthonous units exposed in the French Massif Central, Vosges and Bohemian Massifs (Faure *et al.*, 2009; Skrzypek *et al.*, 2012; Lardeaux *et al.*, 2014; Vanderhaeghe *et al.*, 2020; Martínez Catalán *et al.*, 2021). This includes in particular: (i) the nature of protoliths, which mainly consist of Cambrian–Ordovician magmatic and sedimentary sequences with widespread tholeitic magmatism, marking a stage of rifting along the northern margin of Gondwana (Von Raumer and Stampfli, 2008); (ii) the presence of ophiolitic relics, the largest one being the ophiolite complex of Chamrousse (496 ± 6 Ma) in the massif of Belledonne (Pin and Carme, 1987; Ménot *et al.*, 1988), which lies on the top of the metamorphic nappe pile and is devoid of HP metamorphism (Fig. 1b); (iii) the widespread occurrence of (retrogressed) eclogite/ mafic HP granulite bodies in the high-grade lower crustal units. These eclogites record a Visean (ca. 340–330 Ma) HP stage with peak-P conditions estimated around 1.4–1.8 GPa and 650–750 °C (Rubatto *et al.*, 2010; Jouffray *et al.*, 2020; Jacob *et al.*, 2021, 2022; Vanardois *et al.*, 2022). This eclogitisation is interpreted to mark the closure of a short-lived marginal basin opened in a back-arc setting on the Gondwanian plate at ca. 350 Ma (Fréville *et al.*, 2018), which would have been squeezed and incorporated into the Variscan orogenic wedge by southeastward propagation of deformation during the lower to middle Carboniferous (ca. 340–330 Ma). Felsic protoliths do however not record HP conditions. In weakly to non-migmatized gneiss and micaschist of the Western Pelvoux and Belledonne Massif, a Visean (ca. 340–330 Ma) MP–MT stage (0.5–1.0 GPa and 550–680 °C) is recorded, with no relics of HP/UHP assemblages (Guillot and Ménot, 1999; Fréville *et al.*, 2018, 2022). In the more highly migmatized units of inner Pelvoux, this MP–MT stage is overprinted by a HT granulitic stage at ca. 800–850 °C and 0.6–0.9 GPa (Grandjean *et al.*, 1996; Jacob *et al.*, 2022).

As in the rest of the Variscan metamorphic allochthons, Devonian–Carboniferous magmatism is widespread in the ECM. First, the opening of Devonian–Tournaisian basins, presumably in a back-arc setting, was associated with intrusion of trondhjemite dykes and sills around 350 Ma (Fréville *et al.*, 2018). Tonalite plutons of the same age have also been dated in the Aar Massif (Ruiz *et al.*, 2022). This low-K magmatism was followed by emplacement of K-rich sub-alkaline granitoids from the Visean to the early Permian (ca. 345–295 Ma, Debon and Lemmet, 1999). A large part of this magmatism results from crustal anatexis during the syn- to post-collisional stages, as evidenced by the large proportion of migmatized lower crust

exposed in the ECM (Guillot and Ménot, 2009). However, a significant contribution of mantle sources is also evidenced by the presence of (ultra)-potassic (UHK), mafic to intermediate igneous rocks, locally referred to as durbachites or vaugnerites (Debon *et al.*, 1998; Von Raumer *et al.*, 2014). These rocks are particularly rich both in compatible (Mg, Ni, Cr) and incompatible elements including large-ion lithophile elements (LILE: Cs, Rb, Ba, K), Th, U and light rare earth elements (LREE: La, Ce, Nd, Sm). Thus, the most widely accepted view attributes this magmatism to the melting of an enriched, metasomatized mantle source contaminated by subducted crustal material (Janoušek and Holub, 2007; Soder and Romer, 2018).

## 2.2 The Pelvoux Massif and peridotite-bearing localities

The Pelvoux Massif lies in the eastern part of ECM and exposes high-grade metamorphic series composed of highly migmatized mafic to felsic protoliths (metasediments, orthogneiss, amphibolites), intruded by several granitoid plutons (Fig. 2). The metamorphic series display different sets of N-N150° subvertical and subhorizontal migmatitic foliations, which are consistent with longitudinal flow of partially molten material in a N-N150° sinistral transpressive setting (Fréville *et al.*, 2022; Jacob *et al.*, 2022). However, the early syn-migmatitic fabrics have been partly overprinted by later, sub-solidus deformation occurring in the same sinistral wrenching setting, and by alpine ductile/brittle tectonics (Fig. 2). Occurrence of ultramafic rocks, mainly serpentinites, has been acknowledged for long in the Pelvoux Massif (Pecher, 1970; Le Fort, 1971), as well as in surrounding massifs of Belledonne and Aiguilles Rouges (Barfety *et al.*, 2000; Von Raumer and Bussy, 2004). Inspection of scree in multiple areas indicates that these serpentinites are ubiquitous in the metamorphic basement of the Pelvoux Massif, although in very subordinate volume (Fig. 2). When observed in place, they mostly occur as meter to decameter-wide lenses, which can be up to a hundred meters long. Some of them are concordant with the surrounding metamorphic foliation, but they are commonly transposed in brittle fault zones, as they acted as a weak layer localizing deformation (Fig. 3e). Thus, most of these serpentinites have been strongly deformed and retrogressed, and do not preserve any relics of olivine, pyroxene or any other primary mineral that existed before serpentinization. A second type of peculiar ultramafic rocks was reported by Le Fort (1971) and Pecher (1970) as “amphibolite boulders with a concentric structure”. They occur as decimeter-size dark boulders in migmatites, and consist of a core composed of amphiboles (actinolite-tremolite, hornblende, anthophyllite), chlorite, iron oxides, phlogopite, calcite, and a few relics of spinel and augite, which is surrounded by centimeter-wide reaction coronas of tremolite and biotite. The rocks described by Le Fort (1971) were obviously heavily retrogressed and did not preserve much of their initial assemblages. Thus, their origin remained unclear. They were regarded either as relics of metasomatized calc-silicate boudins, or as ultramafic cumulates of gabbro-dolerite series similar to the Chamrousse ophiolite in the massif of Belledonne.
<figure><img src="image_2.png" /></figure>

**Fig. 2.** (a) Geological map of Central Pelvoux, modified from the BRGM 1:50 000 map of Saint-Christophe en Oisans (Barféty and Pêcher, 1984) with the position and nature of collected samples. Two sites were peridotites have been observed in place are shown with yellow stars. Areas were peridotite enclaves are relatively abundant in screes are highlighted by black dotted lines. (b) Interpretative cross-sections in the Olan-Rouies and Bans areas showing the general structuration of the massif, after Jacob *et al.* (2022) and Fréville *et al.* (2022). Ultramafic enclaves are found both in the flat-lying S1 and in the subvertical S2.

**Fig. 2.** (a) Carte géologique de la partie centrale du massif du Pelvoux, modifiée à partir de la feuille BRGM 1:50 000 de Saint-Christophe en Oisans (Barféty et Pêcher, 1984), avec la position et la nature des échantillons collectés. Deux sites où des péridotites peuvent être observées en place sont indiqués par des étoiles jaunes. Les zones où les enclaves de péridotites sont relativement abondantes dans les éboulis sont mises en évidence par des lignes pointillées noires. (b) Coupes géologiques interprétatives montrant la structure générale du massif, d'après Jacob *et al.* (2022) et Fréville *et al.* (2022). Les enclaves ultramafiques se trouvent à la fois dans les foliations S1 et S2.
<figure><img src="image_3.png" /></figure>

**Fig. 3.** Field pictures of peridotite and serpentinite outcrops in the Pelvoux Massif. (a) Decameter-wide boulder of diatexitic migmatite containing multiple enclaves of peridotite and amphibolite, eastern side of Sellar Pass (N44.8344°; E06.3440°). Enclaves are surrounded by a thin hornblende-bearing coronas at the contact with migmatite. (b) Same kind of peridotite enclaves in diatexitic migmatite on the southern slopes of Olan, Valgaudemar (N44.8509°; E06.2031°). Reaction with the migmatite is very limited here, and enclaves display sharp border devoid of wide coronas. (c, d) Peridotite enclaves in place in a diatexitic migmatite, eastern side of Sellar pass (N44.8361°; E06.3391°). Enclaves display 1–3 cm-wide coronas composed of radially-oriented tremolite surrounded by a second rim of biotite, and are cut by centimeter-wide veins of orthopyroxenite. (e) Serpentinite outcrop at the Sellettes pass. The serpentinites are transposed in an alpine brittle fault oblique to the Variscan foliation (N44.8540°; 06.2181°). (f) Hand-size sample of retrogressed garnet lherzolite (JB-19-33) found in a scree above the Olan hut (see Tab. 1 for coordinates). (g) Close-up view of a phlogopite-bearing peridotite enclave on the large block of picture (a).

**Fig. 3.** Photos des principaux affleurements de péridotite et de serpentinite dans le massif du Pelvoux. (a) Bloc décamétrique dans des migmatite (diatexite), contenant de multiples enclaves de péridotite et d'amphibolite; versant est du col du Sellar (N44,8344°; E06,3440°). Une fine couronne de hornblende entoure les enclaves au contact de la migmatite. (b) Même type d'enclaves de péridotite dans une migmatite diatexitique sur les pentes sud de l'Olan, Valgaudemar (N44,8509°; E06,2031°). La réaction avec la migmatite est ici très limitée, et les enclaves présentent par endroits une bordure nette dépourvue couronnes réactionnelles. (c, d) Enclaves de péridotite en place dans une migmatite diatexitique, côté est du col de Sellar (N44,8361°; E06,3391°). Les enclaves présentent des couronnes réactionnelles de 1 à 3 cm de large, composées de trémolite se développant de façon radiaire, elles-mêmes entourée d'une seconde couronne de biotite. Le coeur des enclaves est recoupé par des veines d'orthopyroxénite de largeur centimétrique. (e) Affleurement de serpentinite au niveau du col des Sellettes. Ces serpentinites sont transposées dans une faille cassante alpine, qui recoupe la foliation varisque (N44,8540°; 06,2181°). (f) Échantillon de lherzolite à grenat rétromorphosée (JB-19-33) trouvé en éboulis au-dessus du refuge d'Olan (coordonnées en Tab. 1). (g) Gros plan sur une enclave de péridotite à phlogopite sur le gros bloc de la photo (a).
**Table 1.** List of collected samples with name, localization (WGS84), petrographic type and bulk composition in major elements recalculated on an anhydrous basis and normalized to 100 wt%. < d.l: values below detection limit. Grt Lrz: Garnet Lherzolite; Spl Hrz: Spinel Harzburgite; Phl Hrz: Phlogopite-chromite Harzburgite; Srp: Serpentinite.

**Tableau 1.** Liste des échantillons collectés avec nom, localisation (WGS84), type pétrographique et composition roche totale éléments majeurs, recalculée sur une base anhydre et normalisée à 100% (en masse). < d.l: valeurs inférieures à la limite de détection. Grt Lrz: Lherzolite à grenat; Spl Hrz: Harzburgite à spinelle; Phl Hrz: Harzburgite à phlogopite-chromite; Srp: Serpentinite.

<table><thead><tr><th colspan="2">Sample</th><th colspan="2">Position (WGS84)</th><th colspan="10">wt% oxides</th></tr><tr><th>Name</th><th>Type</th><th>Long°</th><th>Lat°</th><th>SiO<sub>2</sub></th><th>TiO<sub>2</sub></th><th>Al<sub>2</sub>O<sub>3</sub></th><th>FeO</th><th>MnO</th><th>MgO</th><th>CaO</th><th>Na<sub>2</sub>O</th><th>K<sub>2</sub>O</th><th>P<sub>2</sub>O<sub>5</sub></th><th>Mg#</th><th>Cr#</th></tr></thead><tbody><tr><td>JB-19-20</td><td>Grt Lrz</td><td>6.35604</td><td>44.83521</td><td>45.11</td><td>0.11</td><td>3.29</td><td>8.14</td><td>0.12</td><td>40.70</td><td>2.33</td><td>0.20</td><td>&lt;d.l</td><td>&lt;d.l</td><td>0.90</td><td>0.07</td></tr><tr><td>JB-19-23</td><td>Grt Lrz</td><td>6.34905</td><td>44.83562</td><td>48.19</td><td>0.06</td><td>2.20</td><td>8.67</td><td>0.12</td><td>39.16</td><td>1.41</td><td>0.14</td><td>0.05</td><td>&lt;d.l</td><td>0.89</td><td>0.11</td></tr><tr><td>JB-19-31</td><td>Grt Lrz</td><td>6.22290</td><td>44.84258</td><td>46.52</td><td>0.15</td><td>3.72</td><td>8.15</td><td>0.10</td><td>38.11</td><td>2.69</td><td>0.26</td><td>0.29</td><td>&lt;d.l</td><td>0.89</td><td>0.06</td></tr><tr><td>JB-19-33</td><td>Grt Lrz</td><td>6.20450</td><td>44.84665</td><td>45.33</td><td>0.15</td><td>3.62</td><td>8.52</td><td>0.11</td><td>39.46</td><td>2.58</td><td>0.23</td><td>&lt;d.l</td><td>&lt;d.l</td><td>0.89</td><td>0.06</td></tr><tr><td>JB-19-34</td><td>Grt Lrz</td><td>6.20450</td><td>44.84665</td><td>45.71</td><td>0.13</td><td>3.41</td><td>8.29</td><td>0.08</td><td>40.23</td><td>2.02</td><td>0.13</td><td>&lt;d.l</td><td>&lt;d.l</td><td>0.90</td><td>0.07</td></tr><tr><td>JB-20-06</td><td>Grt Lrz</td><td>6.20309</td><td>44.84967</td><td>46.06</td><td>0.11</td><td>2.78</td><td>8.30</td><td>0.08</td><td>40.54</td><td>1.93</td><td>0.16</td><td>0.04</td><td>&lt;d.l</td><td>0.90</td><td>0.08</td></tr><tr><td>JB-20-12</td><td>Grt Lrz</td><td>6.20473</td><td>44.84611</td><td>45.98</td><td>0.15</td><td>3.80</td><td>8.33</td><td>0.08</td><td>38.71</td><td>2.65</td><td>0.27</td><td>0.04</td><td>&lt;d.l</td><td>0.89</td><td>0.07</td></tr><tr><td>JB-20-13</td><td>Grt Lrz</td><td>6.20319</td><td>44.84945</td><td>45.22</td><td>0.17</td><td>4.49</td><td>8.29</td><td>0.11</td><td>38.60</td><td>2.84</td><td>0.22</td><td>0.06</td><td>&lt;d.l</td><td>0.89</td><td>0.06</td></tr><tr><td>JB-20-14</td><td>Grt Lrz</td><td>6.20274</td><td>44.84982</td><td>45.72</td><td>0.05</td><td>2.20</td><td>8.09</td><td>0.11</td><td>41.98</td><td>1.62</td><td>0.08</td><td>0.15</td><td>&lt;d.l</td><td>0.90</td><td>0.11</td></tr><tr><td>JB-19-21</td><td>Phl Hrz</td><td>6.33910</td><td>44.83612</td><td>45.98</td><td>0.10</td><td>1.03</td><td>8.45</td><td>0.10</td><td>43.21</td><td>0.36</td><td>&lt;d.l</td><td>0.64</td><td>0.13</td><td>0.90</td><td>0.35</td></tr><tr><td>JB-19-22</td><td>Spl.Hrz</td><td>6.33910</td><td>44.83612</td><td>44.69</td><td>0.03</td><td>1.63</td><td>8.79</td><td>0.10</td><td>42.94</td><td>1.42</td><td>0.08</td><td>0.33</td><td>&lt;d.l</td><td>0.90</td><td>0.14</td></tr><tr><td>JB-18-27</td><td>Srp</td><td>6.39256</td><td>44.87250</td><td>45.01</td><td>&lt;d.l</td><td>1.22</td><td>9.11</td><td>0.11</td><td>44.15</td><td>0.39</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>0.90</td><td>0.20</td></tr><tr><td>JB-19-45</td><td>Srp</td><td>6.23450</td><td>44.86690</td><td>46.66</td><td>&lt;d.l</td><td>1.26</td><td>8.17</td><td>0.08</td><td>42.75</td><td>0.79</td><td>&lt;d.l</td><td>0.29</td><td>&lt;d.l</td><td>0.90</td><td>0.18</td></tr><tr><td>JB-19-48</td><td>Srp</td><td>6.10814</td><td>45.23891</td><td>46.20</td><td>0.04</td><td>1.59</td><td>8.98</td><td>0.07</td><td>42.86</td><td>0.28</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>0.89</td><td>0.15</td></tr><tr><td>JB-19-49</td><td>Srp</td><td>6.11245</td><td>45.23940</td><td>46.20</td><td>&lt;d.l</td><td>0.84</td><td>8.84</td><td>0.05</td><td>43.63</td><td>0.43</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>0.90</td><td>0.24</td></tr></tbody></table>

Extensive investigation of the screes in the area covered by the map in Figure 2 has led us to find well-preserved specimens of these boulders, with limited reactions with the migmatites, which contain a primary assemblages dominated by olivine, orthopyroxene and clinopyroxene, with spinel and/or garnet. It is therefore clear that these rocks are retrogressed peridotites. The main occurrences have been found in screes on the southern slopes of Valgaudemar Valley, between the summits of Olan and Les Rouies, and in the Bans Valley. Two sites where they are in place are shown with yellow stars in Figure 2. The first one is located on the eastern side of the Sellar pass, in the Bans Valley (GPS: N44.83612°, E6.33910°, alt. 2770 m), and the second one above Refuge de l'Olan (GPS: N44.85089°, E6.20314°, alt. 2744 m). These peridotites form rounded to sub-angular 5 cm to 1 m wide enclaves. They generally occur in highly-mobilized part of migmatites (diatexites), and commonly appear as "nest" of enclaves, without any particular orientation (Figs. 3a-3c). This disposition suggests the enclaves represent pieces of larger peridotite bodies, possibly up to decimeter-size, which have been fractured during injection of granitic magmas and dismembered into smaller blocks. The enclaves are neither particularly elongated (aspect ratios between 1 and 1.8) nor oriented in any preferential direction, which suggests they did not experience significant ductile deformation in the migmatite. Reaction with the hosting migmatite is evidenced by the development of coronas around the peridotite enclaves. However, they are not systematic. Some enclaves are almost devoid of these coronas, and display instead sharp contacts with the migmatite (Fig. 3b), while hornblende ± biotite reaction coronas may develop in other cases (Figs. 3a and 3c). The most spectacular ones have been observed on peridotites from the Sellar pass, and consist

of a 1–3 cm wide coronas of radially-oriented tremolite, surrounded by a ca. 0.5–1 cm concentric layer of biotite (Figs. 3c and 3d). However, these coronitic reactions only affected the enclaves at their outer rims or along fracturation pathways, leaving the inner parts of largest enclaves (> 10 cm) relatively well-preserved (Figs. 3f and 3g).

## 3 Samples and methods

Eleven samples of peridotite and 4 samples of serpentinite were collected for petrological study and whole-rock geochemical and isotopic analyses. All the samples come from the Pelvoux Massif, excepted two samples of serpentinite collected in a large hectometer-size lense in the Belledonne Massif. Sampling coordinates and petrological classification of each sample are summarized in Table 1.

### 3.1 Petrography and mineral analysis

Backscattered electron (BSE) images were acquired at ISTERRE Grenoble with a Tescan Vega 3 scanning electron microscope equipped with an energy dispersive spectrometer (EDS) for semi-quantitative analysis. It was operated with an accelerating voltage of 16 kV and a beam current of 10 nA. Inclusions in garnet were identified by Raman spectroscopy, using a Horiba Labram Soleil instrument installed at ISTERRE Grenoble, with a 532 nm laser beam. Mineral compositions were measured using a JEOL JXA-8230 electron probe microanalyser (EPMA) at ISTERRE Grenoble, which is equipped with five wavelength dispersive spectrometers (WDS) and an additional EDS. All silicates except olivine
were analyzed using standard EPMA protocols, with an accelerating voltage of 15 kV (amphibole, phlogopite) to 20 kV (garnet, pyroxenes), a beam current of 10 nA (amphibole, phlogopite) to 20 nA (garnet, pyroxenes), and a counting time per element of 30 to 60 s for peak and both background positions. Natural minerals and synthetic glass were used for standardization, and the ZAF matrix procedure was applied for data reduction. The spot size was set to 1–3 µm depending on the size of minerals and the presence of volatile elements.

The composition of olivine was analyzed following high-precision method for major (Si, Mg, Fe), minor (Ni, Mn) and trace (Al, Ca, Cr, Ti, P) elements (Batanova et al., 2015) at an acceleration voltage of 20 kV, beam current of 900 nA and beam diameter of 2 microns. Total analysis time of each point was 7.2 minutes. MongOl sh11-2 olivine reference material (Batanova et al., 2019) was used as primary standard to verify precision and accuracy (e.g. Sobolev et al., 2007). Deviation from the reference measured during the analytical session was between 2 and 13 ppm for trace elements, between 10 and 85 ppm for minor elements and between 0.01 and 0.06 wt% for major oxides, which is in any cases below the 2-sigma uncertainty reported in Batanova et al. (2019). The composition of Cr-spinel was measured at 20 kV and 50 nA beam current. Counting time was adapted to element content and was between 60 to 120 s for each element (for peak and both background positions). The standardization was made using a synthetic oxide standard set (P&H Developments Ltd., Calibration Standards for Electron Probe Microanalysis, Standard Block GEO) for all elements except Mn (on rhodonite). Ferric iron in spinel was calculated assuming perfect stoichiometry. Repeated measurements of the chromite USNM 117075 (Jarosewich et al., 1980) standard and spinel samples Bar 8601-10 and Dar 8502-2, whose $Fe^{3+}/Fe_{total}$ ratio had been measured by Mossbauer spectroscopy (Ionov and Wood, 1992) have shown that the selected method provides $Fe^{3+}/Fe_{total}$ ratios accurate to within the measurement error. Spinel Bar 8601-10 (Ionov and Wood, 1992) was used as a primary standard and was analyzed multiple times during the session to verify precision and accuracy. Deviation from the reference value measured during the session was below 50 ppm for most of trace elements (Ti, V, Mn, Ni) and bit higher for Zn (66 ppm) and Si (86 ppm). For major oxides, deviation from the reference value was about 0.64 wt% for $Al_2O_3$ and 0.06 wt% for FeO and MgO.

## 3.2 Whole-rock composition

Samples were cut to remove the coronas formed at the contact with the migmatites, and domains devoid of large veins were carefully selected to study only the best-preserved domains. They were then crushed in a jaw crusher and pulverized into a 70–80 µm powder using an agate crusher.

whole-rock geochemical analyses were performed at the SARM in Nancy, using a Thermo Fischer iCap6500 ICP-OES for major oxides and an iCapQ ICP-MS for minor and trace elements. Details about the analytical procedures for whole-rock geochemical analyses are described in Carignan et al. (2001). Samples display high loss on ignition between 7 and 13% due to serpentinization, and therefore the analyses were recalculated on an anhydrous basis to allow comparison

between samples. Raw analyses are provided in the Supplementary Material. The R-based GCDKit software (Janoušek et al., 2016) was used for statistical analysis and plotting of whole-rock geochemical data.

## 3.3 Sr-Nd isotopes

Nine samples were selected for Sr-Nd isotope analyses. Due to low Nd concentration, only five of them were analyzed for Nd isotopes. Isotopic measurements were also carried out at the SARM in Nancy. 100 to 200 mg of powdered samples were dissolved into a mixture of concentrated $HNO_3$ and HF heated to 115 °C during 24 to 48 h, followed by concentrated HCl at 125 °C during 24 h.

Sr and Nd were isolated via ion-exchange chromatography using Sr.Spec, Tru.Spec and Ln.Spec resins, following the procedure described in Pin et al. (1994) and Pin and Zalduegui (1997). Sr and Nd isotopes were analyzed by thermal ionization mass spectrometry (TIMS), using respectively a Triton Plus and a Neptune Plus instrument operated in static multi-collection mode. The $^{87}\text{Sr}/^{86}\text{Sr}$ ratios were corrected for mass bias assuming $^{87}\text{Sr}/^{88}\text{Sr} = 0.119400$, and $^{143}\text{Nd}/^{144}\text{Nd}$ ratios were corrected to $^{146}\text{Nd}/^{144}\text{Nd} = 0.721900$ (Luais et al., 1997). Internal standards NBS 987 and JNdi-1 (Tanaka et al., 2000) were used for Sr and Nd, respectively. The decay constants applied to age-correct the isotopic ratios are from Steiger and Jäger (1977) for Sr and Lugmair and Marti (1978) for Nd. The $\epsilon_{Nd_i}$ values were obtained using composition of the chondritic uniform reservoir (CHUR) of Bouvier et al. (2008).

## 4 Petrography

All the samples collected in the ultramafic enclaves have been retrogressed to various degrees and contain a high amount of volatiles (loss on ignition between 7 and 13 wt%), which reflects fluid addition during retrogression. Typical serpentinization textures are observed, with the development of a mesh of serpentine + Fe-oxide in cracks and at olivine grain boundaries. Minor amount of carbonates are also observed in the most retrogressed samples. Four different types of ultramafic rocks have been identified based on petrological observations: serpentinites (4 samples), retrogressed garnet lherzolite (9 samples), retrogressed spinel harzburgite (1 sample) and retrogressed phlogopite-chromite harzburgite (1 sample; Fig. 4). The serpentinites present whole-rock major element compositions similar to the harzburgites (see Sect. 6, Geochemistry Section below), but rarely preserve relics of their primary minerals. They are mostly composed of HT serpentine (antigorite) and Fe-oxides, with a few relics of spinel. Petrological descriptions are therefore focused on the peridotite samples.

### 4.1 The garnet lherzolites

The garnet lherzolites are the most distinctive type in the field. Macroscopically, they are composed of a dark green, fine-grained matrix containing millimeter to centimeter-size reddish spots (former garnet) and millimeter-size grains of
<figure><img src="image_4.png" /></figure>

**Fig. 4.** (a) Thin section scans of representative samples from the three petrographic groups identified in the field (garnet lherzolites, spinel harzburgites, phlogopite-chromite harzburgites). Note the presence of Opx-rich layers in the harzburgites, which may be related to melt-rock interaction. (b) Polarizing microscope images showing the successive stages of recrystallization in the kelyphite formed by destabilization of garnet. Stage 1: Thin kelyphite composed of micrometer-scale intergrowth of spinel + enstatite ± diopside, containing a few relics of garnet. Stage 2: Recrystallization into a coarser-grained symplectite of spinel + enstatite ± diopside. Amphibole may occur in some samples, due to external fluid input during garnet recrystallization. Stage 3: complete annealing of the symplectites, which have recrystallized into a polygonal assemblage of enstatite + diopside, with relics of spinel.

**Fig. 4.** (a) Scans de lames minces d'échantillons représentatifs des trois groupes pétrographiques identifiés sur le terrain (lherzolites à grenat, harzburgites à spinelle, harzburgites à phlogopite-chromite). Notez la présence de niveaux riches en Opx dans les harzburgites, qui semblent résulter d'interactions avec un magma silicaté. (b) Images prises au microscope polarisant, montrant les étapes successives de recristallisation des kelyphites formées par la déstabilisation du grenat. Stade 1 : kelyphites fines, composée d'intercroissance à l'échelle micrométrique de spinelle + enstatite ± diopside, contenant quelques reliques de grenat. Stade 2 : Recristallisation en symplectites plus grossières à spinelle + enstatite ± diopside. De l'amphibole peut être présente dans certains échantillons, en raison d'un apport extérieur de fluide pendant la recristallisation du grenat. Stade 3 : texture de recuit des symplectites, qui ont recristallisé en un assemblage polygonal à enstatite + diopside, avec quelques reliques de spinelle.
diopside (Figs. 3f and 4a). The dark matrix contains 100–200 µm grains of enstatite (labeled Opx-I in Fig. 5a), diospide (Cpx-I) and olivine (Fig. 5a). Garnet is almost completely destabilized, and only a few millimeter-size relics were identified (Fig. 4). Most of the former grains (reddish spots) have been replaced by kelyphites mainly composed of micrometer-size intergrowth of spinel and enstatite (Opx-II in Fig. 5b) with subordinate amounts of larger grains (10–20 µm) of diopside (Cpx-II) and amphibole. These kelyphites are surrounded by a coarse-grained corona of enstatite (Opx-II) + diopside (Cpx-II) (Figs. 5a and 5c). In most samples, the fine-grained kelyphite has partially recrystallized into a coarse-grained globular symplectite of orthopyroxene + spinel ± clinopyroxene ± amphibole and eventually evolved toward a coarse-grained polygonal assemblage of the same minerals with relics of vermicular spinel (Figs. 4b and 5c). The rare garnet relics are directly surrounded by a second generation of symplectite composed of plagioclase + amphibole + Cr-rich spinel (Fig. 5b). The garnet relics contain small inclusions (ca. 20 µm wide) of orthopyroxene, which themselves contain tiny dark inclusions of graphite, identified by Raman spectroscopy (Fig. 5d).

## 4.2 The spinel-pargasite-(phlogopite) harzburgite (JB-19-22)

The harzburgite sample JB-19-22 is mostly composed of olivine, enstatite, pargasitic amphibole and Mg-spinel, with rare diopside, phlogopite, and accessory Fe-oxides (Figs. 4a, 6a and 6b). Spinel forms 200–500 µm grains that contain inclusions of pargasitic amphibole (Fig. 6b), but never occurs in symplectites or kelyphites like in the garnet lherzolites, which suggests that garnet was absent from the primary assemblage. Olivine forms 50–200 µm grains, which have been partly serpentinized (Fig. 6a). Enstatite is abundant and forms large grains up to 500 µm wide. The largest grains are preferentially localized in discontinuous, millimeter-wide bands (Fig. 4a). By contrast, diopside is extremely rare and only a few grains < 100 µm, rimmed by amphibole, were identified (Fig. 6c). Pargasite forms 100–300 µm grains and occurs as smaller inclusions in spinel. Phlogopite is scarce and generally found in association with pargasite.

## 4.3 The phlogopite-chromite harzburgite (JB-19-21)

The phlogopite-chromite harzburgite contains olivine, enstatite, phlogopite, chromian spinel (chromite) and accessory apatite, and is completely devoid of clinopyroxene (Figs. 4a and 6d). The sample investigated in detail (JB-19-21) is an enclave surrounded by a thick corona of biotite and tremolite (Figs. 3c, 3d and 4a), but large grains of phlogopite are also visible in some enclaves devoid of large reaction coronas (Fig. 3g). Olivine forms partly serpentinized grains 200–500 µm in size (Fig. 5d). Enstatite and phlogopite form large grains up to 1 mm-wide. Chromite is ubiquitous. It forms subeuhedral grains ranging from 10 µm to > 100 µm, which occur either in the serpentinized matrix or as inclusions in phlogopite, and contain itself small inclusions of phlogopite (Fig. 6e). Apatite is also closely associated with phlogopite, and generally occurs

in contact with this mineral (Fig. 6d). Harzburgitic domains are cut by centimeter-wide sub-parallel veins of orthopyroxenite (Fig. 3a), which are mainly composed of large grains of enstatite (up to 3 mm wide), with subordinate olivine, phlogopite and chromite (Fig. 6f). Partial replacement of Mg-rich primary olivine by, BSE-brighter Fe-rich secondary olivine is commonly observed in olivine grains located close to or within these orthopyroxenite veins (Fig. 6f). This replacement has mostly occurred in the rims and along cracks, leaving the olivine cores unmodified (Fig. 6f).

# 5 Mineral composition

Mineral composition has been investigated by EPMA in three samples: a garnet lherzolite (JB-19-33), a spinel-pargasite-(phlogopite) harzburgite (JB-19-22) and a phlogopite-chromite harzburgite (JB-19-21). A set of selected representative analyses is provided in Tables 2 and 3. The complete set of electron microprobe analyses is provided in the Supplementary Material.

In all samples, olivine displays high Mg# (molar ratio Mg/Mg + Fe²⁺) between 0.88 and 0.91, high Ni (2700–4200 ppm) and moderate Mn (1000–1300 ppm) content, with Ni/Mn between 2.1 and 3.3 (Fig. 7a), which is indicative of mantle peridotite olivine and precludes any cumulative origin for these samples (Wang et al., 2021). This is consistent with the low concentrations in trace elements: below 50 ppm for Ti and Al and below 200 ppm for Ca, Cr and P. In sample JB-19-21 (phlogopite-chromite bearing harzburgite), the rims of olivine grains within or close to the orthopyroxenite layers have lower Mg# (0.81–0.86) and higher Mn content (Mn = 2000–5800 ppm) than the cores, with low Ni/Mn between 0.32 and 1.74. This suggests this olivine population formed by reaction between mantle olivine and melts. Composition variations are observed between samples. In particular, there is a consistent increase in Ni/Mn correlated with increasing Mg# (Fig. 7a), from the least magnesian olivine in clinopyroxene-garnet lherzolites (Ni/Mn = 2.1–2.3) to the more magnesian olivine in clinopyroxene-free phlogopite-chromite harzburgite (Ni/Mn = 3.0–3.3).

Enstatite displays high Mg# between 0.89 and 0.91. As for olivine, enstatite displays an increase in Mg# correlated with decreasing Mn content from the garnet lherzolite to the phlogopite-chromite harzburgite (Fig. 7b). Al content of enstatite is very low in the phlogopite-chromite harzburgite (ca. 0.01 p.f.u) and higher in the spinel harzurgite (0.09–0.13 p.f.u) and the garnet lherzolite (Opx-I: 0.05–0.08 p.f.u). In the latter sample, secondary enstatite that crystallized in the coronas around garnet (Opx-II) displays higher Al-content up to 0.15 p.f.u, but does not show significant difference in Mg# with Opx-I (Fig. 7c).

Diopside and garnet were analyzed only in the garnet lherzolite. Diopside has high Mg# (0.92–0.94 p.f.u) and variable Al (0.14–0.17 p.f.u Al). Diopside in the coronas (Cpx-II) contains slightly more Al and slightly less Cr than primary diopside Cpx-I, with again no difference in Mg# (Fig. 7d).

Garnet is a slightly chromian pyrope (Cr# = 0.04–0.05). Composition profile (Fig. 7e) reveals a slight Mg-decrease and Fe-increase from core to rim (from Prp₇₃Alm₁₅Grs₁₁ to Prp₆₅Alm₂₃Grs₁₁)
<figure><img src="image_5.png" /><figcaption>Grt Lherzolite (JB-19-33)</figcaption></figure>

**Fig. 5.** Backscattered electron (BSE) and optical microscope images of the garnet lherzolites (sample JB-19-33). (a) BSE wide view on the lherzolite matrix, with a former grain of garnet on the bottom left. Note the abundance of diopside (Cpx-I) in the matrix, typical of a lherzolite. Garnet has been replaced by a kelyphite surrounded by a coarse-grained corona of clinopyroxene + orthopyroxene (Cpx-II + Opx-II). (b) Zoom on an Opx + Spl ± Cpx kelyphite surrounding a relic of garnet (top-right). A thin 20–50 µm-wide symplectite of amphibole + plagioclase ± Cr-spinel has developed at the contact with the garnet relics during a later stage of retrogression. (c) Advanced stage of annealing of a former kelyphite, which has been replaced by a coarse-grained assemblage of amphibole, orthopyroxene, clinopyroxene and Cr-spinel. Relics of a globular spinel + orthopyroxene symplectite is visible in the centre. (d) Close-up view on orthopyroxene inclusions in a garnet relics of sample JB-19-33 (reflected light optical microscopy). Raman spectroscopy on the tiny dark inclusions in orthopyroxene reveals the presence of graphite, indicated by a strong peak around 1566–1575 cm⁻¹ (Beyssac *et al.*, 2002).

**Fig. 5.** Images en microscopie électronique (électrons rétrodiffusés, BSE) et en microscopie optique (lumière réfléchie) d'une lherzolites à grenat (échantillon JB-19-33). (a) Image en BSE de la matrice de la lherzolite, avec un ancien grain de grenat en bas à gauche. Notez l'abondance du diopside (Cpx-I) dans la matrice, typique d'une lherzolite. Le grenat a été remplacé par une kelyphite, entourée d'une couronne à gros grains de clinopyroxène + orthopyroxène (Cpx-II + Opx-II). (b) Zoom sur une kelyphite à Opx + Spl ± Cpx entourant une relique de grenat (en haut à droite). Une fine symplectite de 20–50 µm de large à amphibole + plagioclase ± spinelle chromifère s'est développée au contact des reliques de grenat pendant un stade ultérieur de rétromorphose. (c) Stade avancé de recristallisation d'une ancienne kelyphite, qui a été remplacée par un assemblage à gros grains d'amphibole, d'orthopyroxène, de clinopyroxène et de spinelle chromifère. Les restes d'une symplectite globulaire à spinelle + orthopyroxène sont visibles au centre. (d) Gros plan sur des inclusions d'orthopyroxène dans les reliques de grenat de l'échantillon JB-19-33 (microscopie optique en lumière réfléchie). Les spectres Raman obtenus sur les petites inclusions sombres dans les orthopyroxènes révèle la présence de graphite, indiquée par un pic à environ 1566–1575 cm⁻¹ (Beyssac *et al.*, 2002).
Harzburgites (JB-19-21; JB-19-22)

<figure><img src="image_6.png" /></figure>

**Fig. 6.** Backscattered electron (BSE) images of harzburgites (samples JB-19-21 and JB-19-22). (a) Olivine-dominated assemblage in sample JB-19-22 (Spinel harzburgite), which is almost devoid of clinopyroxene. (b) Large grain of spinel in the spinel harzburgite (JB-19-22). This sample is rich in pargasite, which occurs as 100–200 µm grains and as inclusions in spinel. (c) Rare grain of clinopyroxene surrounded by pargasite in the spinel harzburgite. (d) Mineral assemblage in the phlogopite-chromite harzburgite (JB-19-21). Chromite is ubiquitous and can be found as inclusion in phlogopite. Apatite and dolomite commonly crystallize close to phlogopite. (e) Close-up view on a grain of chromite, which contains a 5-µm wide inclusion of phlogopite. (f) Zoned olivine close to an orthopyroxenite layer, with a Mg-rich core surrounded by a BSE-brighter rim of Fe-Mn-rich olivine. Fe X-ray mapping reveals an irregular boundary between Mg-rich and Fe-rich domains. Note that this zoning is cut by the serpentine network, which indicates that partial replacement of Mg-rich by Fe-Mn-rich olivine predates serpentinization.

**Fig. 6.** Images en microscopie électronique (BSE) des harzburgites (échantillons JB-19-21 et JB-19-22). (a) Assemblage dominé par l'olivine dans l'échantillon JB-19-22 (harzburgite à spinelle), qui est presque dépourvu de clinopyroxène. (b) Gros grain de spinelle dans la harzburgite à spinelle (JB-19-22). Cet échantillon est riche en pargasite, qui se présente sous forme de grains de 100–200 µm et d'inclusions dans les spinelles. (c) Rare grain de clinopyroxène entouré de pargasite dans la harzburgite à spinelle. (d) Assemblage minéral dans la harzburgite à phlogopite-chromite (JB-19-21). La chromite est omniprésente et peut se présenter sous forme d'inclusion dans la phlogopite. L'apatite et la dolomite cristallisent à proximité de la phlogopite. (e) Zoom sur un grain de chromite, qui contient une inclusion de phlogopite de 5 µm de large. (f) Olivine zonée près d'un niveau d'orthopyroxénite, avec un coeur riche en Mg entouré d'une bordure plus claire riche en Fe et Mn. La cartographie du Fe révèle une limite irrégulière entre les domaines riches en Mg et en Fe. On peut noter que la zonation est recoupée par des veines de serpentine, ce qui indique que le remplacement de l'olivine riche en Mg par l'olivine riche en Fe–Mn est antérieur à la serpentinisation.
<caption>**Table 2.** Representative electron microprobe analyses of garnet, orthopyroxene, clinopyroxene, amphibole and phlogopite. na: not analyzed. The complete set of analyzes with end-member recalculations is provided in the Supplementary Material.</caption>

<table><caption>**Tableau 2.** Sélection représentative d'analyses à la microsonde électronique du grenat, de l'orthopyroxène, du clinopyroxène, de l'amphibole et de la phlogopite. na : non analysé. Le jeu de données complet est disponible en matériel supplémentaire.</caption><thead><tr><th>Sample</th><th>JB-19-33</th><th>JB-19-33</th><th>JB-19-21</th><th>JB-19-22</th><th>JB-19-33</th><th>JB-19-33</th><th>JB-19-33</th><th>JB-19-33</th><th>JB-19-21</th><th>JB-19-22</th><th>JB-19-21</th><th>JB-19-21</th></tr></thead><tbody><tr><td>Mineral</td><td>Grt</td><td>Grt</td><td>Opx</td><td>Opx</td><td>Opx</td><td>Opx</td><td>Cpx</td><td>Cpx</td><td>Amp</td><td>Amp</td><td>Phl</td><td>Bt</td></tr><tr><td>LABEL</td><td>Grt-12</td><td>Grt-01</td><td>Opx-02</td><td>Opx-04</td><td>Opx-07</td><td>Opx-08</td><td>Cpx-11</td><td>Cpx-06</td><td>Amp-5</td><td>Amp-1</td><td>Phl-04</td><td>Phl-11</td></tr><tr><td>Position</td><td>core</td><td>rim</td><td>core</td><td>core</td><td>Matrix Core</td><td>Grt Corona</td><td>Grt Corona</td><td>Matrix Core</td><td>Corona</td><td>core</td><td>core</td><td>Corona</td></tr><tr><td>wt% oxides</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SiO<sub>2</sub></td><td>41.32</td><td>41.45</td><td>57.57</td><td>56.02</td><td>55.71</td><td>55.01</td><td>52.99</td><td>53.49</td><td>55.29</td><td>44.74</td><td>41.27</td><td>38.32</td></tr><tr><td>TiO<sub>2</sub></td><td>0.63</td><td>0.61</td><td>0.06</td><td>0.03</td><td>0.09</td><td>0.07</td><td>0.34</td><td>0.31</td><td>0.03</td><td>0.39</td><td>1.15</td><td>2.02</td></tr><tr><td>Al<sub>2</sub>O<sub>3</sub></td><td>21.45</td><td>21.64</td><td>0.15</td><td>2.12</td><td>1.91</td><td>3.57</td><td>3.99</td><td>3.42</td><td>1.23</td><td>12.73</td><td>12.96</td><td>14.50</td></tr><tr><td>Cr<sub>2</sub>O<sub>3</sub></td><td>1.33</td><td>1.30</td><td>0.05</td><td>0.32</td><td>0.28</td><td>0.09</td><td>0.76</td><td>0.87</td><td>0.24</td><td>1.49</td><td>1.03</td><td>0.19</td></tr><tr><td>FeO</td><td>8.65</td><td>11.69</td><td>6.06</td><td>6.78</td><td>7.27</td><td>7.26</td><td>2.04</td><td>2.75</td><td>7.58</td><td>3.63</td><td>2.45</td><td>15.93</td></tr><tr><td>MnO</td><td>0.28</td><td>0.70</td><td>0.16</td><td>0.17</td><td>0.20</td><td>0.18</td><td>0.09</td><td>0.11</td><td>0.42</td><td>0.06</td><td>0.00</td><td>0.21</td></tr><tr><td>MgO</td><td>20.57</td><td>18.25</td><td>35.94</td><td>34.01</td><td>33.73</td><td>33.35</td><td>16.33</td><td>18.05</td><td>19.96</td><td>18.25</td><td>25.39</td><td>15.15</td></tr><tr><td>CaO</td><td>4.41</td><td>4.31</td><td>0.12</td><td>0.21</td><td>0.43</td><td>0.25</td><td>22.39</td><td>19.55</td><td>11.46</td><td>12.34</td><td>0.00</td><td>0.00</td></tr><tr><td>Na<sub>2</sub>O</td><td>0.00</td><td>0.02</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>1.04</td><td>1.03</td><td>0.26</td><td>1.76</td><td>0.15</td><td>0.20</td></tr><tr><td>K<sub>2</sub>O</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>0.13</td><td>0.99</td><td>9.60</td><td>9.36</td></tr><tr><td>P<sub>2</sub>O<sub>5</sub></td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td></tr><tr><td>NiO</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td></tr><tr><td>ZnO</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td></tr><tr><td>Total</td><td>98.63</td><td>99.96</td><td>100.10</td><td>99.67</td><td>99.61</td><td>99.78</td><td>99.95</td><td>99.57</td><td>96.61</td><td>96.39</td><td>93.99</td><td>95.87</td></tr><tr><td>Oxygen<br>Cation</td><td>12</td><td>12</td><td>6</td><td>6</td><td>6</td><td>6</td><td>6</td><td>6</td><td>23</td><td>23</td><td>11</td><td>11</td></tr><tr><td>Si</td><td>2.988</td><td>3.000</td><td>1.980</td><td>1.943</td><td>1.941</td><td>1.911</td><td>1.920</td><td>1.936</td><td>7.695</td><td>6.364</td><td>2.937</td><td>2.852</td></tr><tr><td>Ti</td><td>0.034</td><td>0.033</td><td>0.001</td><td>0.001</td><td>0.002</td><td>0.002</td><td>0.009</td><td>0.008</td><td>0.004</td><td>0.042</td><td>0.062</td><td>0.113</td></tr><tr><td>Al</td><td>1.828</td><td>1.845</td><td>0.006</td><td>0.087</td><td>0.078</td><td>0.146</td><td>0.170</td><td>0.146</td><td>0.201</td><td>2.133</td><td>1.087</td><td>1.272</td></tr><tr><td>Cr</td><td>0.076</td><td>0.074</td><td>0.001</td><td>0.009</td><td>0.008</td><td>0.002</td><td>0.022</td><td>0.025</td><td>0.027</td><td>0.168</td><td>0.058</td><td>0.011</td></tr><tr><td>Fe<sup>2+</sup></td><td>0.523</td><td>0.707</td><td>0.174</td><td>0.197</td><td>0.212</td><td>0.211</td><td>0.062</td><td>0.083</td><td>0.864</td><td>0.432</td><td>0.146</td><td>0.992</td></tr><tr><td>Fe<sup>3+</sup></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.018</td><td>0.000</td><td>—</td><td>—</td></tr><tr><td>Mn</td><td>0.017</td><td>0.043</td><td>0.005</td><td>0.005</td><td>0.006</td><td>0.005</td><td>0.003</td><td>0.003</td><td>0.050</td><td>0.007</td><td>0.000</td><td>0.014</td></tr><tr><td>Mg</td><td>2.217</td><td>1.969</td><td>1.843</td><td>1.759</td><td>1.752</td><td>1.727</td><td>0.882</td><td>0.974</td><td>4.142</td><td>3.869</td><td>2.694</td><td>1.681</td></tr><tr><td>Ca</td><td>0.342</td><td>0.334</td><td>0.005</td><td>0.008</td><td>0.016</td><td>0.009</td><td>0.869</td><td>0.758</td><td>1.709</td><td>1.880</td><td>0.000</td><td>0.000</td></tr><tr><td>Na</td><td>0.000</td><td>0.003</td><td>0.000</td><td>0.000</td><td>0.000</td><td>0.000</td><td>0.073</td><td>0.073</td><td>0.070</td><td>0.485</td><td>0.020</td><td>0.029</td></tr><tr><td>K</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.023</td><td>0.179</td><td>0.871</td><td>0.888</td></tr><tr><td>Total cation</td><td>8.03</td><td>8.01</td><td>4.01</td><td>4.01</td><td>4.01</td><td>4.01</td><td>4.01</td><td>4.01</td><td>16.80</td><td>17.56</td><td>7.87</td><td>7.85</td></tr><tr><td>Mg#</td><td>0.81</td><td>0.74</td><td>0.91</td><td>0.90</td><td>0.89</td><td>0.89</td><td>0.93</td><td>0.92</td><td>1.00</td><td>1.00</td><td>0.95</td><td>0.63</td></tr><tr><td>Cr#</td><td>2.83</td><td>2.85</td><td>0.20</td><td>0.09</td><td>0.09</td><td>0.02</td><td>0.11</td><td>0.15</td><td>0.12</td><td>0.07</td><td>0.05</td><td>0.01</td></tr><tr><td>(Na + K)<sub>tot</sub></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.09</td><td>0.66</td><td>—</td><td>—</td></tr></tbody></table>

probably related to diffusion, while no zoning is observed neither in Ca nor in Cr.

In the well-preserved part of the spinel harzburgite sample JB-19-22, amphibole is a slightly chromian pargasite (Cr=0.15–0.19 p.f.u, Mg# > 0.94 and Al<sup>iv</sup>=1.40–1.76 p.f.u). In sample JB-19-21 (phlogopite-chromite harzburgite), amphibole was only observed in the corona that rims the peridotite enclave. It is a nearly pure tremolite with low Cr-content (Cr=0.02–0.05 p.f.u, Mg# > 0.95, Al<sup>iv</sup>= 0.20–0.33, Fig. 7f).

Spinel in sample JB-19-22 is a slightly chromian Mg-spinel (Mg#=0.63–0.66, Cr#=0.29–0.35, Fig. 7g) with high Ni-content (1250–1500 ppm, Fig. 7h) and low Mn

(1240–1365 ppm), Zn (2700–3200 ppm), V (820–950 ppm) and Ti (300–400 ppm) content. By contrast, sample JB-19-21 contains ferroan chromite (Mg#=0.15–0.40, Cr#=0.81–0.91) with much lower Ni-content (500–740 ppm) but much higher Zn (up to 7400 ppm), V (1000–1250 ppm), Mn (2600–10 050 ppm) and Ti (4500–6200 ppm) content than spinel of sample JB-19-22.

Finally, mica was only analyzed in sample JB-19-21 (Fig. 7i). The harzburgite core of the sample contains nearly pure, Cr-rich–Ti-poor phlogopite (Mg#=0.95, Ti=0.05–0.06 p.f.u, Cr=0.06 p.f.u). In contrast, the coronas around the enclave contain less magnesian, Ti-rich–Cr-poor biotite (Mg#=0.62–0.63, Cr < 0.02 p.f.u, Ti=0.11).
<caption>**Table 3.** Representative EPMA analyses of spinel, chromite and olivine, including trace elements measured with a high intensity electron beam (Batanova et al., 2015). na: not analyzed. The complete set of analyzes is provided in the Supplementary Material.</caption>

<table><caption>**Tableau 3.** Sélection représentative d'analyses à la microsonde électronique du spinelle, de la chromite et de l'olivine, incluant les éléments traces mesurés avec une forte intensité du faisceau d'électrons (Batanova et al., 2015). na : non analysé. Le jeu de données complet est disponible en matériel supplémentaire.</caption><thead><tr><th>Sample</th><th>JB-19-21</th><th>JB-19-22</th><th>JB-19-21</th><th>JB-19-21</th><th>JB-19-22</th><th>JB-19-33</th></tr></thead><tbody><tr><td>Mineral</td><td>Chr</td><td>Spl</td><td>Olivine</td><td>Olivine</td><td>Olivine</td><td>Olivine</td></tr><tr><td>LABEL</td><td>Spl-09</td><td>Spl-01</td><td>Ol-07</td><td>Ol-18</td><td>Ol-06</td><td>Ol-07</td></tr><tr><td>Position</td><td>core</td><td>core</td><td>core</td><td>rim</td><td>core</td><td>core</td></tr><tr><td>wt% oxides</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SiO<sub>2</sub></td><td>0.02</td><td>0.00</td><td>40.83</td><td>39.29</td><td>40.80</td><td>40.58</td></tr><tr><td>TiO<sub>2</sub></td><td>0.93</td><td>0.06</td><td>0.0035</td><td>0.0016</td><td>0.0010</td><td>0.0028</td></tr><tr><td>Al<sub>2</sub>O<sub>3</sub></td><td>4.38</td><td>37.75</td><td>0.0000</td><td>0.0000</td><td>0.0000</td><td>0.0018</td></tr><tr><td>Cr<sub>2</sub>O<sub>3</sub></td><td>58.96</td><td>28.98</td><td>0.0268</td><td>0.0110</td><td>0.0022</td><td>0.0028</td></tr><tr><td>V<sub>2</sub>O<sub>3</sub></td><td>0.17</td><td>0.12</td><td>na</td><td>na</td><td>na</td><td>na</td></tr><tr><td>FeO</td><td>29.47</td><td>18.10</td><td>8.79</td><td>17.27</td><td>9.66</td><td>10.78</td></tr><tr><td>MnO</td><td>0.87</td><td>0.18</td><td>0.1308</td><td>0.7278</td><td>0.1455</td><td>0.1669</td></tr><tr><td>MgO</td><td>3.77</td><td>14.75</td><td>49.54</td><td>42.36</td><td>48.98</td><td>48.11</td></tr><tr><td>CaO</td><td>na</td><td>na</td><td>0.0091</td><td>0.0052</td><td>0.0133</td><td>0.0133</td></tr><tr><td>P<sub>2</sub>O<sub>5</sub></td><td>na</td><td>na</td><td>0.0214</td><td>0.0779</td><td>0.0136</td><td>0.0100</td></tr><tr><td>NiO</td><td>0.07</td><td>0.16</td><td>0.4256</td><td>0.2446</td><td>0.3731</td><td>0.3590</td></tr><tr><td>ZnO</td><td>0.64</td><td>0.40</td><td>na</td><td>na</td><td>na</td><td>na</td></tr><tr><td>Total</td><td>99.29</td><td>100.50</td><td>99.78</td><td>99.99</td><td>100.00</td><td>100.02</td></tr><tr><td>Oxygen</td><td>4</td><td>4</td><td>4</td><td>4</td><td>4</td><td>4</td></tr><tr><td>Cation</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Si</td><td>—</td><td>—</td><td>1.02</td><td>0.97</td><td>1.01</td><td>1.01</td></tr><tr><td>Ti</td><td>0.025</td><td>0.001</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Al</td><td>0.184</td><td>1.268</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Cr</td><td>1.661</td><td>0.653</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Fe<sup>2+</sup></td><td>0.752</td><td>0.343</td><td>0.36</td><td>0.71</td><td>0.40</td><td>0.44</td></tr><tr><td>Fe<sup>3+</sup></td><td>0.127</td><td>0.089</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Mg</td><td>0.200</td><td>0.627</td><td>1.60</td><td>1.35</td><td>1.57</td><td>1.54</td></tr><tr><td>Total cation</td><td>2.95</td><td>2.98</td><td>2.98</td><td>3.03</td><td>2.99</td><td>2.99</td></tr><tr><td>Mg#</td><td>0.90</td><td>0.34</td><td>0.91</td><td>0.81</td><td>0.90</td><td>0.89</td></tr><tr><td>Cr#</td><td>0.21</td><td>0.65</td><td>—</td><td>—</td><td>—</td><td>0.61</td></tr><tr><td>Trace ppm</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Ni</td><td>586</td><td>1259</td><td>3344</td><td>1922</td><td>2932</td><td>2821</td></tr><tr><td>Mn</td><td>6766</td><td>1365</td><td>1013</td><td>5637</td><td>1127</td><td>1293</td></tr><tr><td>Cr</td><td>—</td><td>—</td><td>183</td><td>75</td><td>15</td><td>19</td></tr><tr><td>V</td><td>1155</td><td>840</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Ca</td><td>—</td><td>—</td><td>65</td><td>37</td><td>95</td><td>95</td></tr><tr><td>Ti</td><td>5570</td><td>357</td><td>21</td><td>10</td><td>6</td><td>17</td></tr><tr><td>Si</td><td>90</td><td>0</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Al</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>10</td></tr><tr><td>P</td><td>—</td><td>—</td><td>93</td><td>340</td><td>59</td><td>44</td></tr><tr><td>Zn</td><td>5158</td><td>3213</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Ni/Mn</td><td>0.09</td><td>0.92</td><td>3.30</td><td>0.34</td><td>2.60</td><td>2.18</td></tr></tbody></table>

# 6 Whole-rock geochemistry and Nd–Sr isotopes

## 6.1 Major elements

All the samples have been serpentinized to various degrees during retrogression, but this process is generally considered to

induce no or very minor changes in major element composition, and preserves in particular Al₂O₃/SiO₂ and MgO/SiO₂ ratios (Deschamps et al., 2013). Formation of talc during serpentinization may nevertheless be associated with strong chemical modification, notably SiO₂ enrichment (Paulick et al., 2006). Apart from highly retrogressed or altered domains which have been removed before crushing, no
<figure><img src="image_7.png" /></figure>

**Fig. 7.** Mineral compositions in selected representative samples from the three different peridotite groups. Note in particular the negative correlation between Mg# and Ni–Mn content in olivine and orthopyroxene on panel (a and b), and the significant difference in trace element composition between high Mg# and low-Mg# olivine in the phlogopite-harzburgite (JB-19-21). Distinction between mantle and magmatic olivine is after Wang *et al.* (2021). Note also the sharp difference in composition between the Mg–Cr-rich mantle-phlogopite in the core of the harzburgite enclave and biotite formed by melt-rock interaction in the coronas (g). Mineral compositions are expressed in atom per formula unit, and Mg# is expressed as the molar ratio Mg/Mg + Fe²⁺.

**Fig. 7.** Compositions minérales dans des échantillons représentatifs trois différents groupes de péridotites. Noter en particulier la corrélation négative entre le Mg# et la teneur en Ni-Mn dans l'olivine et l'orthopyroxène (a et b), et la différence significative de composition en éléments traces entre l'olivine riche en Mg et l'olivine pauvre en Mg dans la harzburgite à phlogopite (JB-19-21). Distinction entre olivine mantellique et magmatique d'après Wang *et al.* (2021). Notez également la forte différence de composition entre les phlogopites mantelliques riches en Mg–Cr dans le cœur de l'enclave de harzburgite et les biotites formées par interaction avec la migmatite dans la couronne réactionelle (g). Les compositions minérales sont exprimées en atome par formule structurale, et Mg# est exprimé comme le rapport molaire Mg/Mg + Fe²⁺.
<figure><img src="image_8.png" /></figure>

**Fig. 8.** Whole-rock composition of collected samples recalculated to an anhydrous basis and plotted on the Al₂O₃/SiO₂–MgO/SiO₂ diagram of Paulick *et al.* (2006). All the samples align slightly below the terrestrial mantle melting array (Jagoutz *et al.*, 1979), possibly due to SiO₂ addition during metasomatism (Paulick *et al.*, 2006) or MgO loss during serpentinization. Larger symbols correspond to samples analyzed by electron microprobe.

**Fig. 8.** Composition roche totale des échantillons collectés, recalculée sur une base anhydre et reportée sur le diagramme Al₂O₃/SiO₂–MgO/SiO₂ de Paulick *et al.* (2006). Tous les échantillons s'alignent légèrement en dessous de la ligne « terrestre » (Jagoutz *et al.*, 1979), probablement en raison d'un ajout de SiO₂ pendant le métasomatisme (Paulick *et al.*, 2006) ou d'une perte de MgO pendant la serpentinisation. Les échantillons analysés par microsonde électronique sont mis en évidence par des symboles plus gros.

talc has been observed in the samples. Therefore, the measured compositions re-normalized to an anhydrous basis (Tab. 1) are considered as representative of the whole-rock compositions before serpentinization. Normalizing to dry composition may induce a bias for the phlogopite/pargasite-bearing samples, which were not completely anhydrous before serpentinization. This bias is however assumed to be small considering the very high loss on ignition (8–10 wt%), which is dominated by serpentine and other volatile-bearing phases formed during retrogression rather than pre-existing metasomatic minerals.

In the Al₂O₃/SiO₂–MgO/SiO₂ diagram (Fig. 8), all peridotite and serpentinite analyses plot slightly below the terrestrial array (melting trend of primitive mantle), which could result from addition of silicic melt or metasomatic SiO₂ addition (Paulick *et al.*, 2006). The clinopyroxene-rich garnet lherzolites are particularly fertile and are characterized by high Al₂O₃/SiO₂ > 0.04, low MgO/SiO₂ (0.81–0.91 wt%), and high CaO (1.41–2.69 wt%). By contrast, the serpentinites, the spinel and the phlogopite-chromite harzburgites display low Al₂O₃/SiO₂ < 0.04, high MgO/SiO₂ (0.91–0.96) and low CaO (0.28–1.42 wt%wt%), which is more typical of residual harzburgitic mantle that has experienced partial melting (Fig. 8). The phlogopite–chromite harzburgite JB-19-21 is also enriched in K₂O (0.63 wt%) and P₂O₅ (0.12 wt%) relative to the other samples, which all contain less than 0.35 wt% K₂O and P₂O₅ below detection limit (< 0.1 wt%).

Variation in major element composition reflects the mineralogical composition of the peridotites. In particular, the trend in the Al₂O₃/SiO₂–MgO/SiO₂ diagram is consistent with composition of olivine and orthopyroxene, which are more magnesian in the most refractory harzburgites than in the most fertile garnet lherzolite (Figs. 7a and 7b). Furthermore, higher CaO and Al₂O₃ in the lherzolites is consistent with the large abundance of –retrogressed– garnet and clinopyroxene, and higher K₂O and P₂O₅ in sample JB-19-21 is consistent with the presence of phlogopite and apatite.

## 6.2 Trace elements

Trace element analyses recast on an anhydrous basis are given in Table 4. Large variability is observed in trace element composition, and three main groups are identified based on REE content (Fig. 9a). The first group includes most of the garnet lherzolites, which present LREE-depleted patterns (La/N/Yb_N = 0.46–0.96), with HREE content close to or slightly below that of the primitive mantle (Yb_N = 0.5–1.1). Samples from this group display no or slightly positive Eu anomaly (Eu/Eu* = 1.0–1.2, where Eu/Eu* = Eu_N/ (Sm_N.Gd_N)¹/²), and are not particularly enriched in Th (Th_N/Yb_N = 0.48–2.12), in contrast to the other groups. The second group includes the spinel and the phlogopite harzburgites, plus three samples of
<caption>**Table 4.** Whole-rock trace element composition recalculated on an anhydrous basis. Key race element ratios normalized to the primitive mantle composition (McDonough and Sun, 1995) have also been calculated. < d.l: values below detection limit. Grt Lrz: Garnet Lherzolite; Spl Hrz: Spinel Harzburgite; Phl Hrz: Phlogopite-chromite Harzburgite; Srp: Serpentinite.</caption>

<table><caption>**Tableau 4.** Analyses roche totale en éléments traces, recalculées sur une base anhydre. < d.l. : valeurs inférieures à la limite de détection. Grt Lrz : Lherzolite à grenat; Spl Hrz : Harzburgite à spinelle; Phl Hrz : Harzburgite à phlogopite-chromite; Srp : Serpentinite.</caption><thead><tr><th>Sample</th><th>JB-19-20</th><th>JB-19-23</th><th>JB-19-31</th><th>JB-19-33</th><th>JB-19-34</th><th>JB-20-06</th><th>JB-20-12</th><th>JB-20-13</th><th>JB-20-14</th><th>JB-19-21</th><th>JB-19-22</th><th>JB-18-27</th><th>JB-19-45</th><th>JB-19-48</th><th>JB-19-49</th></tr><tr><th>Type</th><th>Grt Lrz</th><th>Grt Lrz</th><th>Grt Lrz</th><th>Grt Lrz</th><th>Grt Lrz</th><th>Grt Lrz</th><th>Grt Lrz</th><th>Grt Lrz</th><th>Grt Lrz</th><th>Phl Hrz</th><th>Spl.Hrz</th><th>Srp</th><th>Srp</th><th>Srp</th><th>Srp</th></tr></thead><tbody><tr><td>As</td><td>0.92</td><td>90.21</td><td>2.12</td><td>3.88</td><td>3.84</td><td>5.81</td><td>5.16</td><td>1.08</td><td>4.07</td><td>12.30</td><td>2.76</td><td>26.93</td><td>78.75</td><td>496.29</td><td>65.25</td></tr><tr><td>Ba</td><td>22.1</td><td>&lt;d.l</td><td>21.6</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>50.0</td><td>171.5</td><td>26.9</td><td>&lt;d.l</td><td>9.6</td><td>&lt;d.l</td><td>&lt;d.l</td></tr><tr><td>Be</td><td>&lt;d.l</td><td>0.91</td><td>0.55</td><td>&lt;d.l</td><td>0.71</td><td>0.26</td><td>0.95</td><td>0.15</td><td>0.49</td><td>1.18</td><td>0.32</td><td>0.16</td><td>0.50</td><td>0.21</td><td>0.09</td></tr><tr><td>Bi</td><td>&lt;d.l</td><td>0.07</td><td>0.75</td><td>0.09</td><td>0.10</td><td>0.06</td><td>&lt;d.l</td><td>0.08</td><td>0.26</td><td>0.31</td><td>1.21</td><td>0.12</td><td>0.42</td><td>1.09</td><td>0.13</td></tr><tr><td>Cd</td><td>0.03</td><td>0.03</td><td>0.08</td><td>0.04</td><td>&lt;d.l</td><td>0.03</td><td>0.03</td><td>0.06</td><td>0.06</td><td>0.07</td><td>&lt;d.l</td><td>1.89</td><td>0.09</td><td>&lt;d.l</td><td>&lt;d.l</td></tr><tr><td>Co</td><td>101.9</td><td>110.8</td><td>97.9</td><td>103.5</td><td>105.1</td><td>101.4</td><td>100.3</td><td>99.7</td><td>111.8</td><td>129.2</td><td>111.8</td><td>121.2</td><td>113.6</td><td>119.8</td><td>81.0</td></tr><tr><td>Cr</td><td>2718</td><td>2803</td><td>2426</td><td>2499</td><td>2563</td><td>2598</td><td>2710</td><td>2871</td><td>2723</td><td>5674</td><td>2705</td><td>3038</td><td>2815</td><td>2772</td><td>2675</td></tr><tr><td>Cs</td><td>1.52</td><td>7.90</td><td>3.17</td><td>1.30</td><td>3.65</td><td>2.63</td><td>2.56</td><td>0.66</td><td>1.51</td><td>4.02</td><td>5.96</td><td>1.46</td><td>3.56</td><td>0.47</td><td>0.20</td></tr><tr><td>Cu</td><td>7.33</td><td>5.95</td><td>8.85</td><td>8.51</td><td>3.49</td><td>6.55</td><td>5.69</td><td>13.74</td><td>14.13</td><td></td><td>5.59</td><td>5.81</td><td>&lt;d.l</td><td>3.35</td><td>&lt;d.l</td></tr><tr><td>Ga</td><td>3.0</td><td>2.2</td><td>4.0</td><td>3.4</td><td>3.1</td><td>2.6</td><td>3.6</td><td>3.6</td><td>2.5</td><td>2.1</td><td>1.9</td><td>1.1</td><td>1.4</td><td>1.5</td><td>0.8</td></tr><tr><td>Ge</td><td>0.85</td><td>2.18</td><td>1.02</td><td>0.88</td><td>0.93</td><td>0.95</td><td>1.01</td><td>0.92</td><td>0.96</td><td>2.13</td><td>0.88</td><td>1.03</td><td>1.20</td><td>1.42</td><td>1.02</td></tr><tr><td>Hf</td><td>0.14</td><td>0.09</td><td>0.32</td><td>0.20</td><td>0.18</td><td>0.19</td><td>0.20</td><td>0.25</td><td>0.31</td><td>0.55</td><td>0.18</td><td>0.05</td><td>&lt;d.l</td><td>0.13</td><td>&lt;d.l</td></tr><tr><td>Nb</td><td>0.08</td><td>0.05</td><td>0.37</td><td>0.07</td><td>0.09</td><td>0.06</td><td>0.03</td><td>0.08</td><td>0.37</td><td>0.88</td><td>0.63</td><td>0.06</td><td>0.19</td><td>0.18</td><td>0.08</td></tr><tr><td>Ni</td><td>1979</td><td>2129</td><td>1891</td><td>1909</td><td>2038</td><td>2064</td><td>2005</td><td>1865</td><td>2157</td><td>2788</td><td>2104</td><td>2088</td><td>2171</td><td>2201</td><td>2322</td></tr><tr><td>Pb</td><td>1.57</td><td>2.42</td><td>1.57</td><td>0.58</td><td>0.60</td><td>11.88</td><td>5.97</td><td>0.59</td><td>4.91</td><td>1.84</td><td>1.44</td><td>3.15</td><td>2.62</td><td>1.35</td><td>&lt;d.l</td></tr><tr><td>Rb</td><td>1.09</td><td>5.68</td><td>14.81</td><td>2.38</td><td>4.15</td><td>3.36</td><td>4.24</td><td>2.52</td><td>7.76</td><td>34.93</td><td>17.31</td><td>0.99</td><td>17.70</td><td>0.29</td><td>0.28</td></tr><tr><td>Sb</td><td>0.07</td><td>2.00</td><td>0.13</td><td>0.12</td><td>3.09</td><td>2.07</td><td>1.07</td><td>0.35</td><td>0.41</td><td>3.78</td><td>0.63</td><td>0.81</td><td>4.66</td><td>7.69</td><td>9.51</td></tr><tr><td>Sc</td><td>15.0</td><td>11.8</td><td>15.2</td><td>15.5</td><td>15.0</td><td>13.1</td><td>16.4</td><td>19.2</td><td>13.0</td><td>3.2</td><td>12.4</td><td>12.4</td><td>12.3</td><td>10.1</td><td>7.5</td></tr><tr><td>Sn</td><td>&lt;d.l</td><td>&lt;d.l</td><td>1.64</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>0.93</td><td>&lt;d.l</td><td>0.34</td><td>0.66</td><td>0.46</td><td>0.43</td><td>&lt;d.l</td></tr><tr><td>Sr</td><td>9.7</td><td>11.1</td><td>8.0</td><td>10.0</td><td>17.8</td><td>12.3</td><td>11.6</td><td>9.6</td><td>29.9</td><td>7.0</td><td>15.2</td><td>13.8</td><td>9.0</td><td>5.6</td><td>6.6</td></tr><tr><td>Ta</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.1</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>Th</td><td>0.50</td><td>0.08</td><td>0.65</td><td>0.03</td><td>0.03</td><td>0.08</td><td>0.05</td><td>0.48</td><td>1.24</td><td>0.55</td><td>0.08</td><td>0.04</td><td>0.37</td><td>0.08</td><td>0.08</td></tr><tr><td>U</td><td>0.11</td><td>0.05</td><td>0.38</td><td>0.03</td><td>0.02</td><td>0.04</td><td>0.06</td><td>0.11</td><td>0.37</td><td>0.80</td><td>0.32</td><td>0.37</td><td>0.11</td><td>0.10</td><td>0.06</td></tr><tr><td>V</td><td>63.3</td><td>46.5</td><td>69.4</td><td>70.4</td><td>65.6</td><td>55.9</td><td>74.8</td><td>79.0</td><td>53.4</td><td>15.9</td><td>41.1</td><td>34.0</td><td>42.3</td><td>35.5</td><td>27.0</td></tr><tr><td>W</td><td>&lt;d.l</td><td>3.3</td><td>4.3</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>&lt;d.l</td><td>1.3</td><td>11.7</td><td>1.0</td><td>3.0</td><td>6.2</td></tr><tr><td>Y</td><td>2.70</td><td>1.92</td><td>4.90</td><td>3.19</td><td>3.16</td><td>2.28</td><td>3.54</td><td>4.36</td><td>1.57</td><td>2.47</td><td>1.28</td><td>0.66</td><td>0.31</td><td>1.16</td><td>0.33</td></tr><tr><td>Zn</td><td>56.8</td><td>61.0</td><td>68.4</td><td>54.4</td><td>53.9</td><td>59.7</td><td>71.5</td><td>62.1</td><td>72.4</td><td>99.9</td><td>54.1</td><td>87.4</td><td>49.3</td><td>67.0</td><td>62.2</td></tr><tr><td>Zr</td><td>4.1</td><td>2.8</td><td>10.8</td><td>7.2</td><td>6.2</td><td>6.7</td><td>6.7</td><td>8.5</td><td>10.4</td><td>17.7</td><td>6.5</td><td>2.2</td><td>&lt;d.l</td><td>4.8</td><td>&lt;d.l</td></tr><tr><td>La</td><td>0.930</td><td>0.293</td><td>1.053</td><td>0.235</td><td>0.210</td><td>0.186</td><td>0.434</td><td>0.329</td><td>1.356</td><td>1.771</td><td>1.949</td><td>0.389</td><td>0.152</td><td>0.750</td><td>0.503</td></tr><tr><td>Ce</td><td>1.866</td><td>0.696</td><td>2.499</td><td>0.762</td><td>0.665</td><td>0.601</td><td>0.999</td><td>0.986</td><td>4.099</td><td>4.206</td><td>4.936</td><td>0.581</td><td>0.397</td><td>1.543</td><td>0.936</td></tr><tr><td>Pr</td><td>0.234</td><td>0.106</td><td>0.367</td><td>0.141</td><td>0.117</td><td>0.110</td><td>0.159</td><td>0.172</td><td>0.658</td><td>0.589</td><td>0.620</td><td>0.088</td><td>0.052</td><td>0.200</td><td>0.118</td></tr><tr><td>Nd</td><td>1.039</td><td>0.526</td><td>1.810</td><td>0.807</td><td>0.665</td><td>0.599</td><td>0.839</td><td>0.895</td><td>3.067</td><td>2.577</td><td>2.266</td><td>0.359</td><td>0.190</td><td>0.833</td><td>0.464</td></tr><tr><td>Sm</td><td>0.291</td><td>0.202</td><td>0.608</td><td>0.298</td><td>0.250</td><td>0.232</td><td>0.317</td><td>0.353</td><td>0.859</td><td>0.615</td><td>0.402</td><td>0.075</td><td>0.043</td><td>0.204</td><td>0.101</td></tr><tr><td>Eu</td><td>0.090</td><td>0.084</td><td>0.158</td><td>0.121</td><td>0.120</td><td>0.102</td><td>0.146</td><td>0.139</td><td>0.176</td><td>0.085</td><td>0.084</td><td>0.018</td><td>0.038</td><td>0.111</td><td>0.082</td></tr><tr><td>Gd</td><td>0.356</td><td>0.252</td><td>0.707</td><td>0.405</td><td>0.362</td><td>0.299</td><td>0.427</td><td>0.497</td><td>0.604</td><td>0.525</td><td>0.270</td><td>0.083</td><td>0.031</td><td>0.196</td><td>0.082</td></tr><tr><td>Tb</td><td>0.064</td><td>0.047</td><td>0.128</td><td>0.075</td><td>0.066</td><td>0.056</td><td>0.086</td><td>0.098</td><td>0.069</td><td>0.079</td><td>0.039</td><td>0.014</td><td>0.005</td><td>0.032</td><td>0.011</td></tr><tr><td>Dy</td><td>0.465</td><td>0.326</td><td>0.839</td><td>0.505</td><td>0.473</td><td>0.404</td><td>0.596</td><td>0.719</td><td>0.311</td><td>0.473</td><td>0.227</td><td>0.091</td><td>0.046</td><td>0.207</td><td>0.062</td></tr><tr><td>Ho</td><td>0.107</td><td>0.074</td><td>0.185</td><td>0.117</td><td>0.111</td><td>0.089</td><td>0.139</td><td>0.166</td><td>0.058</td><td>0.094</td><td>0.047</td><td>0.023</td><td>0.012</td><td>0.046</td><td>0.013</td></tr><tr><td>Er</td><td>0.294</td><td>0.216</td><td>0.509</td><td>0.330</td><td>0.317</td><td>0.254</td><td>0.396</td><td>0.477</td><td>0.156</td><td>0.242</td><td>0.131</td><td>0.066</td><td>0.040</td><td>0.126</td><td>0.036</td></tr><tr><td>Tm</td><td>0.049</td><td>0.033</td><td>0.080</td><td>0.053</td><td>0.050</td><td>0.040</td><td>0.061</td><td>0.074</td><td>0.023</td><td>0.037</td><td>0.022</td><td>0.012</td><td>0.007</td><td>0.020</td><td>0.005</td></tr><tr><td>Yb</td><td>0.312</td><td>0.219</td><td>0.532</td><td>0.344</td><td>0.327</td><td>0.267</td><td>0.393</td><td>0.482</td><td>0.164</td><td>0.257</td><td>0.153</td><td>0.085</td><td>0.059</td><td>0.123</td><td>0.035</td></tr><tr><td>Lu</td><td>0.049</td><td>0.036</td><td>0.083</td><td>0.051</td><td>0.055</td><td>0.041</td><td>0.061</td><td>0.076</td><td>0.025</td><td>0.042</td><td>0.027</td><td>0.015</td><td>0.011</td><td>0.019</td><td>0.007</td></tr><tr><td>(Th/Nb)<sub>N</sub></td><td>54.857</td><td>13.843</td><td>14.605</td><td></td><td></td><td>2.671</td><td>4.917</td><td>21.938</td><td>5.354</td><td>10.796</td><td>11.709</td><td>7.288</td><td>11.542</td><td>1.711</td><td>16.961</td><td>7.748</td></tr><tr><td>(Th/Yb)<sub>N</sub></td><td>9.328</td><td>2.119</td><td>7.099</td><td></td><td></td><td>0.484</td><td>0.737</td><td>1.191</td><td>0.644</td><td>16.988</td><td>27.870</td><td>20.825</td><td>5.306</td><td>3.768</td><td>17.542</td><td>12.749</td></tr><tr><td>(Th/La)<sub>N</sub></td><td>4.358</td><td>2.204</td><td>5.000</td><td></td><td></td><td>1.051</td><td>1.474</td><td>1.500</td><td>1.315</td><td>2.869</td><td>5.635</td><td>2.273</td><td>1.619</td><td>2.030</td><td>3.996</td><td>1.240</td></tr><tr><td>(La/Yb)<sub>N</sub></td><td>2.140</td><td>0.962</td><td>1.420</td><td>0.490</td><td>0.460</td><td>0.500</td><td>0.794</td><td>0.490</td><td>5.921</td><td>4.946</td><td>9.162</td><td>3.278</td><td>1.856</td><td>4.390</td><td>10.279</td></tr></tbody></table>
<figure><img src="image_9.png" /><figcaption>**Fig. 9.** Trace element composition of collected samples, normalized to the composition of primitive mantle after McDonough and Sun (1995). Three groups are identified based on their trends in rare earth elements (see description in text). See Figure 8 for the legend.</figcaption></figure>

<caption>**Fig. 9.** Composition en éléments traces des échantillons collectés, normalisée à la composition du manteau primitif d'après McDonough et Sun (1995). Trois groupes sont identifiés à partir des spectres en terres rares (voir description dans le texte). Légende identique à la figure 8.</caption>

garnet lherzolite, which present similar REE patterns. They display a HREE content similar to samples of the first group (Yb$_N$=0.4–1.2), but they are strongly enriched in LREE and Th (La$_N$/Yb$_N$=1.42–9.16, Th$_N$/Yb$_N$=7.10–27.87) and present a negative Eu anomaly (Eu/Eu* = 0.46–0.85). The third group consists of serpentinites. They are all strongly depleted in HREE relative to the primitive mantle (Yb$_N$ = 0.08–0.28), and display either U-shaped REE patterns, with depletion in MREE relative to HREE and LREE, or continuous enrichment from HREE to LREE (La$_N$/Yb$_N$ = 1.86–10.28, Th$_N$/Yb$_N$ = 3.77–17.54). With one exception, all the analyzed serpentinites have a positive Eu anomaly (Eu/Eu* up to 3.2). These three groups share in common a strong enrichment in LILE and other fluid-mobile incompatible elements, in particular Cs, Rb, U and Pb, and strongly negative anomalies in HFSE, in particular in Nb and Ta, independently of their mineral composition or REE and Th content (Fig. 9b). Th$_N$/Nb$_N$ ratio, which quantifies this anomaly, ranges between 1.71 and 54.86.

## 6.3 Sr–Nd isotopes

Nine samples have been selected for bulk-rock Sr–Nd isotopic analyses (Tab. 5). However, due to very low Nd content, no Nd isotopic data could be acquired for samples from the LREE-depleted group, and the only available Nd data are from LREE-enriched samples. Initial εNd$_i$ have been recalculated at 330 Ma, which corresponds to the climax of the Variscan collision in the ECM. Variations of ±20 Ma on the initial age have little effect on recalculated εNd$_i$ (±0.07–0.23) compared to the observed variation range. εNd$_{330}$ values cover a large range from slightly positive to strongly negative values (–8.12; +0.59). Slight contamination by granitic melt in the migmatite cannot be excluded, which would shift εNd$_{330}$ toward lower values. However, only the freshest part of the enclaves, which did not present macroscopic evidences of melt-rock interaction in the migmatite, has been selected for whole-rock isotopic analyses. Therefore, we consider that the measured whole-rock isotopic compositions are dominated by
**Table 5.** Whole-rock Sr–Nd isotopic compositions. Initial εNd have been recalculated at 330 Ma, which corresponds to the climax of the Variscan collision in the ECM. Sensitivity of εNdi to variation in initial age is shown by recalculating isotopic compositions at 350 Ma and 300 Ma (values in bracket). Recalculation of initial $^{87}$Sr/$^{86}$Sr is hampered by high uncertainty related to high Rb/Sr ratios in certain samples, and possible contamination by fluids or melts in the crust.

**Tableau 5.** Compositions isotopiques en Sr–Nd sur la roche totale. Les εNd initiaux ont été recalculés à 330 Ma, ce qui correspond au pic de la collision varisque dans les Massifs Cristallins Externes. La sensibilité de εNdi à la variation de l’âge initial est testée en recalculant les compositions isotopiques à 350 Ma et 300 Ma (valeurs entre parenthèses). Le calcul des rapports $^{87}$Sr/$^{86}$Sr initiaux est entravé par une forte incertitude, liée à des rapports Rb/Sr élevés dans certains échantillons, et à une contamination possible par des fluides ou des liquides silicatés dans la croûte.

<table><thead><tr><th rowspan="2">Sample</th><th rowspan="2">Type</th><th colspan="4">Composition (ppm)</th><th colspan="2">Sr isotopes</th><th colspan="3">Nd isotopes</th></tr><tr><th>Rb</th><th>Sr</th><th>Nd</th><th>Sm</th><th><sup>87</sup>Rb/<sup>86</sup>Sr</th><th><sup>87</sup>Sr/<sup>86</sup>Sr</th><th><sup>147</sup>Sm/<sup>144</sup>Nd</th><th><sup>143</sup>Nd/<sup>144</sup>Nd</th><th>εNd<sub>330</sub></th><th>εNd<sub>[350-310]</sub></th></tr></thead><tbody><tr><td>JB-19-20</td><td>Grt Lrz</td><td>1.09</td><td>9.7</td><td>1.039</td><td>0.291</td><td>0.32445</td><td>0.70895</td><td>0.16938</td><td>0.51261</td><td>0.59</td><td>[+0.78; +0.64]</td></tr><tr><td>JB-19-33</td><td>Grt Lrz</td><td>2.38</td><td>10.0</td><td>0.807</td><td>0.298</td><td>0.68983</td><td>0.70440</td><td></td><td></td><td></td><td></td></tr><tr><td>JB-20-12</td><td>Grt Lrz</td><td>4.24</td><td>11.6</td><td>0.839</td><td>0.317</td><td>1.05645</td><td>0.71017</td><td></td><td></td><td></td><td></td></tr><tr><td>JB-20-13</td><td>Grt Lrz</td><td>2.52</td><td>9.6</td><td>0.895</td><td>0.353</td><td>0.75968</td><td>0.70572</td><td></td><td></td><td></td><td></td></tr><tr><td>JB-20-14</td><td>Grt Lrz</td><td>7.76</td><td>29.9</td><td>3.067</td><td>0.859</td><td>0.75138</td><td>0.71442</td><td>0.16924</td><td>0.51216</td><td>-8.12</td><td>[-7.92; -8.06]</td></tr><tr><td>JB-19-21</td><td>Phl Hrz</td><td>34.93</td><td>7.0</td><td>2.577</td><td>0.615</td><td>14.4460</td><td>0.76771</td><td>0.14426</td><td>0.51224</td><td>-5.64</td><td>[-5.38; -5.65]</td></tr><tr><td>JB-19-22</td><td>Spl.Hrz</td><td>17.31</td><td>15.2</td><td>2.266</td><td>0.402</td><td>3.30482</td><td>0.72185</td><td>0.10713</td><td>0.51230</td><td>-2.78</td><td>[-2.44; -2.89]</td></tr><tr><td>JB-19-48</td><td>Srp</td><td>0.29</td><td>5.6</td><td>0.833</td><td>0.204</td><td>0.15114</td><td>0.71727</td><td>0.14816</td><td>0.51246</td><td>-1.37</td><td>[-1.13; -1.37]</td></tr><tr><td>JB-19-49</td><td>Srp</td><td>0.28</td><td>6.6</td><td>0.464</td><td>0.101</td><td>0.12105</td><td>0.71783</td><td></td><td></td><td></td><td></td></tr></tbody></table>

<figure><img src="image_10.png" /></figure>

**Fig. 10.** (a) Isochemical phase diagram computed with PerpleX 6.9.0 (Connolly, 2009). P–T conditions are inferred using the composition quality factor $Q_{\text{cmp}}$ (Duesterhoeft and Lanari, 2020), which is computed over the P–T grid for the primary Grt + Opx-I + Cpx-I assemblage (b). Best fit between modeled and measured composition ($Q_{\text{cmp}}$ close to 1, in yellow) is obtained around 940–980 °C and 2–4 GPa and is consistent with results obtained by inverse thermobarometry (red lines). Another diagram with P–T estimations obtained from standard composition isopleths is provided in the Supplementary Material. Inverse thermobarometry on Opx-II + Cpx-II coronas around garnet indicate lower temperature of equilibration around 800–850 °C (blue lines). TA98: Taylor (1998); BK90: Brey and Köhler (1990); NG85: Nickel and Green (1985).

**Fig. 10.** (a) Diagramme de phase (pseudosection) calculé avec PerpleX 6.9.0 (Connolly, 2009). Les conditions P–T sont estimées à l’aide du critère de composition $Q_{\text{cmp}}$ (Duesterhoeft et Lanari, 2020), calculé sur l’ensemble sur la grille P–T pour l’assemblage primaire à Grt + Opx-I + Cpx-I (b). Le meilleur ajustement entre la composition modélisée et celle mesurée (indiqué par une valeur de $Q_{\text{cmp}}$ proche de 1, en jaune) est obtenu autour de 940–980 °C et 2–4 GPa, et est cohérent avec les résultats obtenus par thermobarométrie inverse (lignes rouges). Un diagramme montrant les estimations P–T obtenues avec la méthode standard (intersection des isoplèthes de composition) est disponible en matériel supplémentaire. La thermobarométrie inverse sur les couronnes secondaires autour du grenat (Opx-II + Cpx-II) indique une température d’équilibre plus basse, autour de 800–850 °C (lignes bleues). TA98: Taylor (1998); BK90: Brey et Köhler (1990); NG85: Nickel et Green (1985).
**Table 6.** Summary of P–T estimates in sample JB-19-33. For inverse thermobarometry, reported P and T are average of several measurements. Reported uncertainties of ±50 °C and ±0.5 GPa are conservative estimates taking into account composition variation and internal calibration uncertainty (Nimis and Grütter, 2010).

**Tableau 6.** Résumé des estimations P–T sur l’échantillon JB-19-33. Pour la thermobarométrie inverse, les valeurs de P et T reportées, sont des moyennes obtenues à partir de plusieurs mesures dans chaque domaine. Les incertitudes rapportées de ±50°C et ±0.5 GPa sont des estimations conservatrices prenant en compte les variations de composition et les incertitudes internes des différents thermobaromètres (Nimis et Grütter, 2010).

<table><thead><tr><th rowspan="2">Assemblage</th><th colspan="3">Inverse thermobarometry</th><th colspan="2">PerpleX</th></tr><tr><th>P (GPa)</th><th>T (°C) TA98</th><th>T (°C) BK90</th><th>T (°C)</th><th>P (GPa)</th></tr></thead><tbody><tr><td>Opx-I + Cpx-I + Grt</td><td>3.0 ± 0.5 (NG85)</td><td>973 ± 50</td><td>1137 ± 50</td><td>940–980</td><td>2.0–4.0</td></tr><tr><td>Opx-II + Cpx-II (coronas)</td><td>1.0 (assumed)</td><td>827 ± 50</td><td>803 ± 50</td><td>nd</td><td>0.8–1.3</td></tr></tbody></table>

the peridotite component, and are thus representative of the Variscan mantle composition.

In contrast, Sr isotope data are more difficult to interpret. Measured ⁸⁷Sr/⁸⁶Sr ratios range from 0.7057 to 0.7677, with highest values measured in the harzburgites (0.7218–0.7677) and in the serpentinites (0.7173–0.7178). Recalculation of initial isotopic composition at 330 Ma may be flawed by late alteration, in particular in the serpentinites, which present high ⁸⁷Sr/⁸⁶Sr ratios. Furthermore, high Rb/Sr ratios (17 to 35, Tab. 5) in the phlogopite-bearing samples make recalculation of initial ⁸⁷Sr/⁸⁶Sr very sensitive to even small analytical errors on Rb or Sr contents or age variations. Thus, the Sr isotopic compositions cannot be easily interpreted in terms of mantle processes, and we will not discuss these results further.

## 7 Thermobarometry

The P–T evolution of the garnet lherzolites has been constrained using inverse thermobarometers based on equilibrium between garnet, clinopyroxene and orthopyroxene. Among the collected samples, only JB-19-33 preserved a complete assemblage with these three phases. P–T conditions have therefore been estimated only on this sample. Two metamorphic stages can be constrained: (i) a first stage in the garnet stability field, which we constrained using the composition of garnet, Opx-I and Cpx-I; (ii) a lower P stage in the spinel stability field, which we constrained using the composition of Opx-II–Cpx-II pairs in the coarse-grained coronas around the kelyphites. We only used composition data acquired in the core of largest grains (ca. 200 µm), which are the least affected by post-growth diffusion. As pointed out by Nimis and Grütter (2010), many formulations of thermobarometric equations commonly used for mantle rocks are not internally consistent and may yield biased results, even in well-equilibrated samples. We followed their recommendations and used the Al-in-Opx barometer of Nickel and Green (1985) (NG85) in combination with the two-pyroxene thermometer of Taylor (1998) (TA98) and the Ca-in-Opx thermometer of Brey and Köhler (1990) (BK90) which have proven to yield the most robust results. Furthermore, these thermobarometers do not rely on Fe–Mg partitioning between co-existing minerals, which is particularly sensitive to diffusion resetting (Carlson, 2006; Cherniak and Dimanov, 2010).

Inverse thermobarometry has been completed by thermodynamic modeling with PerpleX 6.9.0 (Connolly, 2009), using the thermodynamic database of Holland and Powell (2011) and the set of solution models from Jennings and Holland (2015) developed for phase relations in peridotite in the composition space NaO–CaO–FeO–MgO–Al₂O₃–SiO₂–Cr₂O₃. All iron was assumed to be ferrous (no Fe³⁺), and the bulk composition of sample JB-19-33 recalculated on an anhydrous basis was used in input. All P–T results are summarized in Table 6. The phase diagram computed with PerpleX is shown in Figure 10a. The primary assemblage composed of olivine, garnet, clinopyroxene and orthoproxene is stable over a large P–T range above 1.3–2.0 GPa, while spinel appears at lower P between 0.8 and 1.3 GPa. Equilibration conditions of garnet, Cpx-I and Opx-I has been assessed by computing the composition quality factor Q<sub>cmp</sub> of Duesterhoef and Lanari (2020) over the whole P–T grid (Fig. 10b). Q<sub>cmp</sub> gives the match between modeled and measured mineral composition in function of P and T, using the same mineral compositions as those used for inverse thermobarometry. The closer Q<sub>cmp</sub> is to 1 (yellow shades in Fig. 10b), the better the match. The best composition match (Q<sub>cmp</sub> > 0.95) is obtained for P between 2.0 and 4.0 GPa and T between 940 and 980 °C (Fig. 10b). Inverse thermobarometry using the NG85 and TA98 calibrations yields a comparable estimate for the primary assemblage at P=3.0 GPa and T=973 °C. Internal calibration uncertainty of these thermobarometers can realistically be assumed to be around ±0.5 GPa for P and ±50 °C for T (Nimis and Grütter, 2010). However, assuming the same pressure (3.0 GPa), the Ca-in-Opx thermometer (BK90) yields higher T at 1136 °C. In contrast, Opx-II–Cpx-II pairs in the coronas yield consistent T of 827 ± 50 °C for the TA98 two-pyroxene thermometer and 803 ± 50 °C for the BK90 Ca-in-Opx thermometer, assuming P = 1.0 GPa.

By contrast with the lherzolites, spinel/phlogopite-chromite harzburgites lack of garnet and clinopyroxene, which prevents the use of classical thermobarometers. Moreover, these samples display clear evidence of modal metasomatism, which makes the identification of equilibrium assemblages for thermobarometry more challenging. Therefore, no P–T estimates have been obtained on these samples. Lack of garnet nevertheless indicates equilibration below the spinel/garnet transition. However, according to Ziberna *et al.* (2013), this transition may be significantly shifted up to > 5 GPa in
Cr-rich depleted harzburgites. Therefore, lack of garnet in the
metasomatized harzburgites does not provide any information
about initial P, which may well be in the same range as that of
the garnet lherzolites.

8 Discussion

**8.1 Chemical overprint of enclaves in the crust**

The relatively small (< 1 m) peridotite enclaves hosted in
migmatites have experienced significant chemical overprint,
first at HT conditions in the partially molten crust, and later by
low-T fluid addition during retrogression. Precise characteri-
zation of this late overprint is needed before discussing more
extensively about the mantle composition and mantle
processes recorded by these enclaves. Chemical interaction
with the migmatite host, either at supra-solidus or sub-solidus
conditions, is evidenced by the development of reaction
coronas at the contact of the enclaves. However, many
enclaves display sharp contact with the migmatites, or are
surrounded only by a thin rim of biotite and/or hornblende
(Figs. 3a, 3b and 3g), which suggests that chemical interactions
with the host have been restricted to a thin outer rim, with little
effect in the core of enclaves. Migmatite-peridotite interactions
are however particularly visible around some harzburgite
enclaves, especially on sample JB-19-21, which is surrounded
by spectacular radial coronas composed of biotite + tremolite
(Figs. 3c and 3d). Furthermore, Opx-rich veins observed in this
sample (Fig. 4a) are interpreted as a product of the reaction Ol
+ Liq 1 → Opx + Liq 2, in which primary olivine is dissolved
by a relatively SiO₂-rich melt (Liq 1), producing orthopyrox-
ene and a more mafic melt (Liq 2). These veins could therefore
indicate reaction of the peridotite with granitic liquid in the
migmatite. However, such reactions only affected centimeter-
wide domains, as shown by the distribution of olivine grains
with Fe-rich rims, which are restrained to the immediate
vicinity of Opx-rich veins. Since both the coronas and the
metasomatic veins have been carefully removed by cutting
before doing geochemical analyses, they should only
marginally affect the results of the geochemical analyses.

Serpentinization and other low T retrogression reactions
occurred to various degrees in all the collected samples, and
resulted in significant addition of volatile components (Loss on
Ignition (LOI) up to 13 wt%). However, it does not seem that
these reactions have significantly affected the composition of
minerals. In sample JB-19-21 (phlogopite harzburgite), Fe-
Mg-Mn zoning previously developed in olivine was not
modified by the development of the serpentine mesh (Fig. 6f),
which suggests little or no effect of serpentinization on olivine
composition. Fluid-addition does not seem to be associated
either with strong modification of the whole-rock trace element
and isotopic composition. Positive anomalies in some large ion
elements mobile in crustal fluids (U, Pb, Cs) are observed in all
samples, independently of the volatile content in these
samples, as shown with U/Th plotted against LOI in Figure 11.
Two serpentinite outliers with high LOI however display high
U/Th relative to other samples, which might be the result of U
addition by high U/Th crustal fluids. Negative HFSE
anomalies (Th/Nb > 1), LREE enrichment relative to HREE
(La/N/YbN) and Nd content are also not correlated with the LOI

(Fig. 11), which suggest that the main geochemical features of
the peridotites and serpentinites were acquired before
serpentinization. Moreover, the strongly negative anomalies
in HFSE correspond to Nb and Ta contents well below the
primitive mantle composition for most of the samples
(Fig. 9b). Thus, even accounting for possible addition of
La, Th and U by crustal fluids, these anomalies likely existed
before the chemical overprint in the crust.

8.2 P–T evolution of the garnet lherzolites

The primary Grt + Cpx-I + Opx-I assemblage in sample
JB-19-33 correspond to an initial equilibration stage of the
lherzolites in the garnet field, estimated at 973 ± 50 °C and
3.0 ± 0.5 GPa using the NG85 + TA98 thermobarometers,
which is consistent with estimation obtained by forward
thermodynamic modeling (2–4 GPa, 940–980 °C). However,
there is quite a large discrepancy between the TA98 two-
pyroxene and the BK90Ca-in-Opx thermometer, which yields
T ca. 150 °C difference. This discrepancy cannot be explained
solely by internal uncertainties of the respective calibrations,
but more probably results also from composition resetting by
diffusion after mineral growth. This may explain for instance
the relatively large spread in Mg# and Ca observed for
clinopyroxene (Fig. 7d, Table S1). In contrast, Al content in
Opx presents little dispersion (Fig. 7c), which suggests little or
no diffusion has occurred. We therefore consider the Al-in-
Opx barometry is quite robust. Taking the TA98 and BK90 as
respectively lower and upper estimates for T, and considering
the dependence in T of the NG85 barometer, we estimate that
the initial equilibration stage occurred around 970–1140 °C
and 3.0–4.0 GPa.

Formation of kelyphites is indicative of a second
metamorphic recrystallization stage at lower pressure.
Kelyphitization results from reaction between garnet and
olivine in the spinel stability field, following the reaction
Grt + Ol (+H₂O) → Opx + Cpx + Spl (+Amp) (Godard and
Martin, 2000; Obata, 2011), which occurred between 1.3
and 0.8 GPa according to our thermodynamic calculations
(Fig. 10). The coarse-grained Opx-II + Cpx-II coronas
developed concomitantly with the kelyphite (Obata, 2011),
and constrain T between 800 and 850 °C during decompression
using the TA98 and BK90 thermometers. Finally, the
formation of plagioclase-bearing symplectites around garnet
relics indicates further decompression in the plagioclase
stability field below ca. 0.8 GPa. Transformation of the
kelyphites into coarser symplectites and eventually into a
polygonal assemblage of orthopyroxene + clinopyroxene ±
spinel ± amphibole (Fig. 4b) is a process indicative of static
recrystallization at relatively high T conditions > 800 °C
(Obata, 2011). Resulting textures are very similar to annealing
textures observed in garnet-clinopyroxene HP mafic granulites
of the Oisans–Pelvoux Massif (Jacob et al., 2022), which have
experienced an overprint at similar P–T conditions (0.6–
0.9 GPa; 800–870 °C) during the late Variscan collisional
stages (ca. 310–295 Ma). These textures therefore suggest that
the lherzolite enclaves have remained in the lower crust at
T > 800 °C, following the first stage of decompression that led
to garnet breakdown.
<figure><img src="image_11.png" /><figcaption>**Fig. 11.** Whole-rock (U/Th)N, (Th/Nb)N, (La/Yb)N and NdN normalized to the primitive mantle composition in function of the loss on ignition (LOI). LOI is considered as a proxy for fluid addition during retrogression in the crust (mostly serpentinization). Lack of correlation with the LOI ((U/Th)N, (Th/Nb)N, (La/Yb)N) or slightly negative correlation (Nd) suggests that late fluid-addition during retrogression had little effect on whole-rock trace element compositions. Two serpentinite outliers with high (U/Th)N (contoured spots) might have nevertheless experienced significant contamination by fluid. See caption of Figure 8 for the legend. Primitive mantle composition after McDonough and Sun (1995).</figcaption></figure>

<caption>**Fig. 11. Sélections de paramètres géochimiques sur roche totale (U/Th)N, (Th/Nb)N, (La/Yb)N et NdN (composition normalisée au manteau primitif) en fonction de la perte au feu (LOI).** Celle-ci est utilisée comme un indicateur de l'intensité des réactions de rétrogression tardives (serpentinisation principalement), associées à un fort apport de fluide. L'absence de corrélation avec la LOI ((U/Th)N, (Th/Nb)N, (La/Yb)N) ou la présence de corrélation légèrement négative (Nd) suggère que l'apport tardif de fluides pendant la rétrogression a eu peu d'effets sur la composition en éléments traces de la roche totale. Légende identique à la figure 8. Composition du manteau primitif d'après McDonough and Sun (1995).</caption>
## 8.3 A refractory mantle refertilized by melts and fluids

Geochemical and mineralogical investigations of the peridotites and serpentinites reveal significant decoupling between major, trace elements and mineral compositions, which we interpret as the result of metasomatic overprint of variably melt-depleted refractory mantle. Initial melt depletion is inferred from whole-rock mineralogy, Al₂O₃/SiO₂ and MgO/SiO₂ compositions, which present large variations from very fertile (high Al₂O₃/SiO₂, cpx-rich) garnet lherzolites to refractory (low Al₂O₃/SiO₂, cpx-poor to cpx-free) harzburgites and serpentinites (Fig. 8). These variations are furthermore consistent with the composition of olivine and orthopyroxene, which display increasing Mg# and Ni and decreasing Mn from the more fertile to the more refractory samples. Finally, even in the most fertile lherzolites, LREE depletion suggests they originally derive from a refractory mantle.

The variable degrees of melt-depletion are however completely decorrelated from the whole-rock trace element compositions. Indeed, largest enrichment in incompatible elements, in particular LREE, Rb, Ba and Th are observed in the refractory harzburgites, while the more fertile lherzolites are depleted in LREE relative to HREE (Fig. 7). Selective enrichment in some LILE (Cs, Rb, U, Pb) relative to HFSE (Nb–Ta) is also observed independently of the mineral and whole-rock composition of the samples, and thus marks a cryptic metasomatic overprint. Some samples also record modal metasomatism, as evidenced by the presence of volatile-bearing phases like phlogopite and pargasite in the refractory harzburgites. In the lherzolites, hydrous minerals (pargasite) are also found, but they are mostly located in former garnet sites (kelyphites or polygonal mosaic assemblages), and may be related to retrogression in the crust rather than mantle metasomatic processes. However, the tiny graphite-bearing inclusions found in garnet relics (Fig. 5) suggest the presence of reduced C–O–H fluids during garnet crystallization. Therefore, it is possible that crystallization of garnet-clinopyroxene-rich assemblages in the LREE-depleted lherzolites was catalyzed by the presence of these reduced fluids.

Similar decoupling between major and trace elements content in orogenic mantle rocks has been widely documented throughout the Variscan Belt (Scambelluri *et al.*, 2006; Medaris *et al.*, 2015; Kubeš *et al.*, 2022). It has commonly been interpreted as the result of an early phase of partial melting, followed by one or multiple phases of metasomatic overprint that led to selective enrichment of the bulk rock in some melt or fluid-mobile trace elements (Medaris *et al.*, 2015; Kubeš *et al.*, 2022). The initial refractory character of the mantle rocks may either indicate derivation from a subducted oceanic lithosphere (Kubeš *et al.*, 2022), or point toward old (possibly Proterozoic) partial melting events in the sub-continental mantle (Medaris *et al.*, 2015). Decoupling between major and trace element composition may result from limited metasomatic influx, which would only produce selective enrichment in specific elements mobile in the metasomatic phase, leaving the initial major element characteristics of the bulk rock largely unmodified. Such a case has been documented by Borghini *et al.* (2018) in mantle eclogites from the Granulitgebirge (NW Bohemian Massif). In these rocks, percolation by felsic melts rich in Cs, Li, B, Pb, Rb, Th, and U is evidenced by the presence of tiny crystallized melt inclusions in garnet, while

the bulk rock only recorded very limited enrichment in these elements.

Two different metasomatic trends (subsequently referred to as type-I and type-II) are recorded in the Pelvoux mantle enclaves, which variably affected the different samples: all of them present variable enrichment in LILE (Cs, Rb, U, Pb) relative to Nb and Ta (type-I), while only the harzburgites and some lherzolites are enriched in incompatible LREE and Th (type-II). This dichotomy appears more clearly when plotting Th/Nb (amplitude of Nb anomaly) against Th/Yb (enrichment in incompatible elements relative to HREE). The LREE-depleted lherzolite samples are aligned along a very steep trend (Fig. 12), which indicates increasing Nb anomalies at nearly constant Th/Yb (type-I). In contrast, the LREE-enriched samples are generally more scattered and define a much shallower trend (type-II), which indicates significant enrichment in incompatible elements (Th and LREE) relative to more incompatible ones (HREE), without additional fractionation of LILE relative to Nb and Ta. We interpret these two trends as the result of interactions with different types of fluids or melts, which present different affinities with LREE relative to HREE and LILE relative to HFSE.

The signature of type-I metasomatism is very similar to subduction-related metasomatism, which as been extensively discussed in the literature (e.g. Baier *et al.*, 2008; Hermann and Rubatto, 2009; Zheng, 2019). Strong enrichment in LILE relative to Nb and Ta in the metasomatic phase is probably related to melting and/or dehydration at UHP conditions (> 2.5 GPa) of subducted metasediments, in which the stability at HP of phengite, rutile, allanite and monazite exerts a strong control on LILE, Th, LREE and HFSE content of the fluid or melt produced (Hermann and Rubatto, 2009). This type of signature is very common in Variscan mantle (U)HP rocks (Scambelluri *et al.*, 2006; Medaris *et al.*, 2015; Borghini *et al.*, 2018, 2020; Kubeš *et al.*, 2022), which suggests that metasomatism mainly resulted from percolation by silicic melts released during Variscan subduction. This is furthermore confirmed by Lu–Hf and Sm–Nd garnet ages of some metasomatized garnet peridotites which point toward Devonian–Carboniferous age of metasomatism and UHP metamorphism (Tumiati *et al.*, 2003; Ackerman *et al.*, 2020; Kubeš *et al.*, 2022).

Type-II metasomatism presents distinctive characteristics, which can be investigated with more details using mineral compositions and petrographic observations in the phlogopite-chromite harzburgite (sample JB-19-21), one of the most LREE-enriched samples (Fig. 9). Petrographic relations in this sample show that metasomatism is associated with coeval crystallization of phlogopite, chromite and apatite. Composition of low-Mg# chromite, which is poor in compatible Ni, and comparatively richer in more incompatible Ti and Mn (Figs. 7g and 7h, Tab. 3), suggests it crystallized from a melt phase. We therefore suggest that type-II metasomatism was driven by percolation of the harzburgite by a hydrous silicate melt rich in K, P, Cr and LREE. This signature, which is characterized by simultaneous enrichment in compatible (Cr) and incompatible elements, is typical of ultrapotassic magmas like lamproites or lamprophyres (Foley *et al.*, 1987). It appears that similar ultrapotassic magmas are quite widespread in the Variscan basement of ECM (Debon *et al.*, 1998) and other Variscan Massifs of Western and Central Europe. With their more
<figure><img src="image_12.png" /></figure>

**Fig. 12.** (Th/Nb)N vs. (Th/Yb)N diagram (a), and NbN vs. LaN diagram (b) showing the two metasomatic trends (type-I and type-II) observed in the samples. Compositions are normalized to the primitive mantle of McDonough and Sun (1995), to show the relative evolution of the different groups from this reference. Type-I metasomatism is characterized by increasing HFSE anomaly with little or no fractionation of incompatible elements (LREE + Th) relative to HREE. This corresponds to an increase in (Th/Nb)N at constant (Th/Yb)N in Figure 12a, and to an increase in La at constant Nb in Figure 12b. In constrast, type-II metasomatism is characterized by strong enrichment in incompatible elements without fractionating fluid-mobile LREE and fluid-immobile HFSE. This corresponds to an increase in (Th/Yb)N in Figure 12a, and to a correlated increase in La and Nb in Figure 12b. According to Paulick *et al.* (2006), these two trends are indicative of fluid-rock and melt-rock mantle interaction for type-I and type-II metasomatism, respectively. See caption of Figure 8 for the legend.

**Fig. 12.** Diagramme (Th/Nb)N vs. (Th/Yb)N (a) et NbN vs. LaN (b), montrant les deux tendances métasomatiques observées dans les échantillons (type-I et type-II). Les compositions sont normalisées par rapport au manteau primitif, afin de montrer l'évolution relative des différents groupes à partir de cette référence. Le métasomatisme de type-I est caractérisé par une augmentation des anomalies en HFSE avec peu ou pas de fractionnement des éléments incompatibles (LREE + Th) par rapport aux HREE. Cela se traduit par une augmentation du rapport (Th/Nb)N sans modification majeure du rapport (Th/Yb)N (Fig. 12a), et par une augmentation de la teneur en La sans évolution notable de la teneur en Nb (Fig. 12b). À l'inverse, le métasomatisme de type II est caractérisé par un fort enrichissement en éléments incompatibles, sans fractionnement notable en fonction de la mobilité dans les fluides aqueux (LREE vs. HFSE). Cela correspond à une augmentation du rapport (Th/Yb)N (Fig. 12a), et à une augmentation corrélée de teneurs en La et en Nb (Fig. 12b). D'après Paulick *et al.* (2006), ces deux tendances sont les marqueurs respectifs d'interaction du manteau avec un fluide aqueux (type-I) et avec un liquide magmatique (type II). Légende identique à la figure 8.

differentiated counterparts (Mg–K-rich monzodiorite to qtz-
monzonite and mela-syenites), they collectively form a
magmatic series referred to as durbachites (Von Raumer
et al., 2014; Moyen et al., 2017; Soder and Romer, 2018).
These durbachites presumably derived from melting of a
metasomatized mantle source, enriched by melts derived from
a subducted mature continental crust (Janoušek and Holub,
2007; Janoušek et al., 2019, 2020). We therefore speculate that
type-II metasomatism is related to percolation in the mantle of
the primary magmas from which the durbachites are derived.
This hypothesis is supported by experimental work of Foley
and Pertermann (2021), who produced a phlogopite harzbur-
gite composition similar to our sample by dynamic percolation
of lamproite melt into a garnet peridotite at ca. 3.0 GPa.
Finally, the very low εNdᵢ measured in the phlogopite
harzburgite sample is also similar to that of durbachites
(Janoušek et al., 2019) and confirms involvement of mature
continental material in the source of the metasomatic melts.
Variable Nd isotopic signatures (−8.12 < εNdᵢ < +0.59) in the
LREE-enriched peridotites suggest that the extent of type-II
metasomatism was spatially very heterogeneous, and probably
focused along localized melt migration pathways.

8.4 Geodynamic considerations

The ultramafic enclaves from the Pelvoux Massif display
a large petrographic variety, and thus probably sample
different mantle levels in the garnet and spinel stability field,
which were variably affected by metasomatic overprint. UHP
conditions recorded in the garnet lherzolite (3.0–4.0 GPa;
970–1140 °C) exceed by far the peak-P conditions of crustal
eclogites and HP granulites from the ECM, which are
typically between 1.4 and 1.8 GPa (Ferrando et al., 2008;
Jouffray et al., 2020; Jacob et al., 2021, 2022; Vanardois
et al., 2022). The MP–HT conditions associated with
decompression in the spinel field (0.8–1.3 GPa and 800–
850 °C) are however quite consistent with the granulite-facies
overprint (0.6–0.9 GPa; 800–870 °C) recorded in mafic HP
granulites from the Pelvoux Massif, which presumably
occurred during the late collisional stages around 325–
300 Ma (Jacob et al., 2022). Because of the contrasted P–T
evolution, we suggest that the peridotite enclaves were
incorporated in the lower crust only during the late stages of
the Variscan collision in the ECM (ca. 325–300 Ma), and are
not pre-orogenic enclaves.
Decompression of the garnet lherzolites from > 3 GPa to
0.8–1.3 GPa requires significant vertical uplift of at least 50–
70 km, assuming lithostatic pressure and a density of
3300 kg.m⁻³ for the lithospheric mantle. However, the exact
mechanism for vertical transport of the mantle enclaves
remains unclear. Buoyancy-driven exhumation in a subduction
channel, or incorporation as tectonic slices in the downgoing
continental slab is commonly suggested to explain the
presence of mantle enclaves in crustal units (Guillot et al.,
2009; Brueckner et al., 2010; Malatesta et al., 2012). However,
the Pelvoux Massif lies hundreds of kilometers south of the
main Variscan (Saxo-Thuringian/Rheic) subduction zones
(Fig. 1), and there is no clear evidence of a former subduction
channel preserved in the ECM (Jacob et al., 2021, 2022).
Either a subduction was active in the ECM during the early
orogenic stages, and was then dismembered by subsequent
tectonics and magmatism, forming a “cryptic suture zone”
(Schulmann, Lexa et al., 2014; Jouffray et al., 2020), or such a
subduction never existed, and alternative mechanisms have to
be involved. In the Moldanubian Zone of the Bohemian
Massif, “relamination” mechanisms have been suggested to
explain the presence of (U)HP units far from the plate
interface. These (U)HP units consist of felsic garnet-kyanite-
bearing granulites with relics of coesite and/or diamond
(Kotková et al., 2011; Perraki and Faryad, 2014), which
contain numerous bodies of garnet peridotites and other UHP
mantle rocks (Medaris et al., 2005). It is suggested that these
units represent buoyant *mélanges* derived from the slab/mantle
wedge interface, which were exhumed by diapiric ascent far
behind the subduction front, and eventually incorporated into
the overlying lower crust. This relaminated lower crust would
have then been locally extruded by up-doming (Janoušek and
Holub, 2007; Lexa et al., 2011; Schulmann, Catalán et al.,
2014). This mechanism has been corroborated by recent
thermomechanical modeling (Maierová et al., 2018, 2021),
and could be a viable explanation to the presence of mantle
enclaves in the lower crust of ECM, far from the subduction
trench. Furthermore, it is rather consistent with (i) the extent of
metasomatism observed in the enclaves, and (ii) the sampling
of different mantle domains at different depths. However, as
opposed to the Bohemian Moldanubian Zone, no evidences of
UHP conditions in the felsic crust have been found so far. It is
therefore still not clear whether the felsic migmatitic host has
experienced such extreme conditions. Recent tomographic
reconstruction of the crust and upper mantle of the Western and
Central Alps indicates the presence of a low-velocity anomaly
in the lower crust below the ECM (Nouibat et al., 2022), which
might sign the presence of relaminated felsic material.
However, extrusion of this lower crust may have been less
efficient than in the Bohemian Massif, which would explain
the scarcity of mantle enclaves and the lack of large UHP felsic
granulite units in the ECM.

Extensive metasomatic interactions with hydrous melts
and/or fluids recorded in the ultramafic enclaves are very likely
related to slab/mantle interactions during the Variscan orogeny,
although we lack of geochronological data to constrain
precisely the age. Indeed, older mantle-derived magmatism in
the ECM does not present any evidence of mantle enrichment.
In particular, Cambrian–Ordovician tholeiitic magmatism,
which is widespread in the ECM, consists only of N/E-MORBs
with positive εNdᵢ between +6 and +8 indicative of a depleted

mantle source (Paquette et al., 1989). Thus, the mantle below
the ECM prior to the Variscan orogeny was probably mostly
composed of depleted, refractory harzburgitic material. The
initial melt-depletion event inferred from whole-rock and
mineral compositions of the peridotites may be related to this
widespread magmatic event, although older depletion events
are not excluded. The isotopic and trace element enrichment
observed in the peridotite enclaves was necessarily acquired
later, more probably during Variscan subduction and collision.
This is in good agreement with the metamorphic and
metasomatic evolution of other Variscan orogenic peridotites.
Dating of garnet in lherzolites from the Bohemian Massif
(Kubeš et al., 2022) and from the Ulten Zone (Tumiati et al.,
2003) yield indeed Carboniferous ages (ca. 340–330 Ma),
which points toward a widespread metasomatic event in the
mantle occurring during the Variscan collisional stages. This
age is very similar to the age of emplacement of the durbachite
magmas (Debon et al., 1998), which clearly suggests a
connection between this peculiar magmatism and metasoma-
tism observed in mantle rocks.

9 Conclusion

Enclaves of peridotites and serpentinites discovered in
exhumed lower and mid-crustal domains in the Oisans-
Pelvoux and Belledonne Massif provide insight into the
composition and evolution of the Variscan mantle in the ECM.
Petrological and geochemical diversity among the enclaves
mark the important heterogeneity of the mantle domain from
which they derive. Thermobarometric estimations indicate
provenance of these enclaves in the upper mantle wedge or the
overlying lithospheric mantle. Refractory compositions,
characterized in particular by high MgO/SiO₂ and low
Al₂O₃/SiO₂ in the whole-rock and high Mg# and Ni-content
in olivine, mark an early melt-depletion event, presumably
related to the widespread Cambrian–Ordovician magmatism
recorded all across the Variscan Belt. Depleted mantle domains
were variably metasomatized by a LILE-rich phase,
probably derived from a Variscan subduction zone, which
produced selective enrichment in Cs, Rb, U, Pb relative to
less mobile Nb and Ta. Some samples display in addition
modal metasomatism (crystallization of phlogopite and/or
pargasite with accessory chromite and apatite), associated
with a strong enrichment in LREE relative to HREE. This
metasomatism is attributed to percolation by a hydrous
LILE–K–P–Cr-rich melt, which may be genetically related
with the ultrapotassic magmas of the durbachite series.
These geochemical characteristics are in line with whole-
rock Nd isotopic compositions, which indicate enrichment
of the mantle by a continental crust component. Further
work is however needed to constrain the precise timing of
metasomatic events, understand the exhumation mecha-
nisms of the mantle enclaves and establish more precise
correlations with the Variscan magmatic and tectono-
metamorphic evolution in the ECM.

**Supplementary Material**

**Figure S1.** P–T phase diagram with selected composition isopleths in sample JB-19-33 (garnet lherzolite).
**Table S1.** Complete set of *in situ* electron microprobe analyses of minerals.

**Table S2.** Complete set of whole-rock geochemical data, including raw data and data recalculated to an anhydrous basis.

The Supplementary Material is available at http://www.bsgf.fr/10.1051/bsgf/2022021/olm.

*Acknowledgements.* This work was supported by the BRGM through the Référentiel Géologique de la France program (RGF). Jean-Baptiste Jacob was funded by a doctoral grant from the ENS de Lyon. We greatly thank Laurent Jolivet and Olivier Vanderhaeghe for editorial handling, and two anonymous reviewers for their useful suggestions and comments, in particular regarding the comparison with Moldanubian mantle rocks. We would like to warmly thank Arnaud Pêcher, who worked about fifty years ago during his PhD and guided Jean-Baptiste Jacob to some key outcrops. We also would like to thank Marianna Jagercikova, who assisted Jean-Baptiste Jacob in the field during the summer 2019, and actually found some of the best-preserved samples. We wish to thank all the technical staff at ISTerre in Grenoble for the preparation of samples and technical help during the analytical sessions. Most of the field investigations were carried out in the central zone of the Ecrins National park, in which collecting rock samples is usually strictly forbidden. We therefore wish to thank the administration staff of the park, who gave us authorization to collect the samples.

## References

Ackerman L, Kotková J, Čopjaková R, Sláma J, Trubač J, Dillingerová V. 2020. Petrogenesis and Lu-Hf dating of (ultra) mafic rocks from the Kutná Hora Crystalline Complex: Implications for the Devonian evolution of the Bohemian Massif. *Journal of Petrology*. https://doi.org/10.1093/petrology/egaa075.

Altherr R. 2021. Retrograded garnet peridotites from Col des Bagenelles and Crébimont in the Variscan Vosges Mountains (NE France). *Contributions to Mineralogy and Petrology* 176(7): 53. https://doi.org/10.1007/s00410-021-01811-7.

Altherr R, Kalt A. 1996. Metamorphic evolution of ultrahigh-pressure garnet peridotites from the Variscan Vosges Mountains (France). *Chemical Geology* 134(1): 27–47. https://doi.org/10.1016/S0009-2541(96)00088-5.

Altherr R, Soder CG. 2018. Sapphirine as a breakdown product of garnet in a Variscan UHP/HT Peridotite from the Vosges Mountains (France)— An indication of near-isothermal decompression. *Journal of Petrology* 59(11): 2221–2243. https://doi.org/10.1093/petrology/egy096.

Ballèvre M, Martínez Catalán JR, López-Carmona A, Pitra P, Abati J, Díez Fernández R, et al. 2014. Correlation of the nappe stack in the Ibero-Armoricán arc across the Bay of Biscay: a joint French-Spanish project. *Geological Society, London, Special Publications* 405: 77–113. https://doi.org/10.1144/SP405.13.

Baier J, Audôtat A, Keppler H. 2008. The origin of the negative niobium tantalum anomaly in subduction zone magmas. *Earth and Planetary Science Letters* 267(1): 290–300. https://doi.org/10.1016/j.epsl.2007.11.032.

Barféty JC, Pêcher A. 1984. Geological map and explanatory text of the sheet n 822 "St Christophe en Oisans", Scale: 1:50 000. Ed. BRGM.

Barféty JC, Gidon M, Ménot RP, Debon F, Pêcher S, Guillot S. 2000. Notice de la carte géologique de la France, feuille Domène (773), Scale 1: 50 000. Orléans : BRGM.

Batanova VG, Sobolev AV, Kuzmin DV. 2015. Trace element analysis of olivine: High precision analytical method for JEOL JXA-8230 electron probe microanalyser. *Chemical Geology* 419: 149–157. https://doi.org/10.1016/j.chemgeo.2015.10.0.

Batanova VG, Thompson JM, Danyushevsky LV, Portnyagin MV, Garbe-Schönberg D, Hauri E, et al. 2019. New olivine reference material for *in situ* microanalysis. *Geostandards and Geoanalytical Research* 43(3): 453–473. https://doi.org/10.1111/ggr.12266.

Bellahsen N, Mouthereau F, Boutoux A, Bellanger M, Lacombe O, Jolivet L, et al. 2014. Collision kinematics in the western external Alps. *Tectonics* 33(6): 1055–1088. https://doi.org/10.1002/2013TC003453.

Bellanger M, Augier R, Bellahsen N, Jolivet L, Monié P, Baudin T, et al. 2015. Shortening of the European Dauphinois Margin (Oisans Massif, Western Alps): New insights from RSCM maximum temperature estimates and $^{40}$Ar/$^{39}$Ar *in situ* dating. *Journal of Geodynamics* 83: 37–64. https://doi.org/10.1016/j.jog.2014.09.004.

Beyssac O, Goffé B, Chopin C, Rouzaud JN. 2002. Raman spectra of carbonaceous material in metasediments: A new geothermometer. *Journal of Metamorphic Geology* 20: 859–871.

Borghini A, Ferrero S, Wunder B, Laurent O, O’Brien PJ, Ziemann MA. 2018. Granitoid melt inclusions in orogenic peridotite and the origin of garnet clinopyroxenite. *Geology* 46(11): 1007–1010. https://doi.org/10.1130/G45316.1.

Borghini A, Ferrero S, O’Brien PJ, Laurent O, Günter C, Ziemann MA. 2020. Cryptic metasomatic agent measured *in situ* in Variscan mantle rocks: Melt inclusions in garnet of eclogite, Granulitgebirge, Germany. *Journal of Metamorphic Geology* 38(3): 207–234. https://doi.org/10.1111/jmg.12519.

Bouvier A, Vervoort JD, Patchett PJ. 2008. The Lu-Hf and Sm-Nd isotopic composition of CHUR: Constraints from unequilibrated chondrites and implications for the bulk composition of terrestrial planets. *Earth and Planetary Science Letters* 273(1-2): 48–57.

Brey GP, Köhler T. 1990. Geothermobarometry in four-phase Lherzolites II. New thermobarometers, and practical assessment of existing thermobarometers. *Journal of Petrology* 31(6): 1353–1378. https://doi.org/10.1093/petrology/31.6.1353.

Brueckner H. 2018. The great eclogite debate of the Western Gneiss Region, Norwegian Caledonides: The *in situ* crustal vs. exotic mantle origin controversy. *Journal of Metamorphic Geology* 36(5): 517–527. https://doi.org/10.1111/jmg.12314.

Brueckner H, Medaris LG. 2000. A general model for the intrusion and evolution of “mantle” garnet peridotites in high-pressure and ultra-high-pressure metamorphic terranes: Intrusion and evolution of ‘mantle’ garnet peridotites. *Journal of Metamorphic Geology* 18 (2): 123–133. https://doi.org/10.1046/j.1525-1314.2000.00250.x.

Brueckner HK, Carswell DA, Griffin WL, Medaris LG, Van Roermund HLM, Cuthbert SJ. 2010. The mantle and crustal evolution of two garnet peridotite suites from the Western Gneiss Region, Norwegian Caledonides: An isotopic investigation. *Lithos* 117(1): 1–19. https://doi.org/10.1016/j.lithos.2010.01.011.

Carignan J, Hild P, Mevelle G, Morel J, Yeghicheyan D. 2001. Routine analyses of trace elements in geological samples using flow injection and low pressure on-line liquid chromatography coupled to ICP-MS: A study of geochemical reference materials BR, DR-N, UB-N, AN-G and GH. *Geostandards Newsletter* 25 (2-3): 187–198. https://doi.org/10.1111/j.1751-908X.2001.tb00595.x.
Carlson WD. 2006. Rates of Fe, Mg, Mn, and Ca diffusion in garnet. *American Mineralogist* 91(1): 1–11. https://doi.org/10.2138/am.2006.2043.

Chen Y, Su B, Guo S. 2015. The Dabie–Sulu orogenic peridotites: Progress and key issues. *Science China Earth Sciences* 58(10): 1679–1699. https://doi.org/10.1007/s11430-015-5148-9.

Cherniak DJ, Dimanov A. 2010. Diffusion in pyroxene, mica and amphibole. *Reviews in Mineralogy and Geochemistry* 72(1): 641–690. https://doi.org/10.2138/rmg.2010.72.14.

Connolly JAD. 2009. The geodynamic equation of state: What and how. *Geochemistry, Geophysics, Geosystems* 10(10). https://doi.org/10.1029/2009GC002540.

Debon F, Lemmet M. 1999. Evolution of Mg/Fe ratios in Late Variscan Plutonic rocks from the External Crystalline Massifs of the Alps (France, Italy, Switzerland). *Journal of Petrology* 40(7): 1151–1185. https://doi.org/10.1093/petroj/40.7.1151.

Debon F, Guerrot C, Ménot R-P. 1998. Late Variscan granites of the Belledonne Massif (French Western Alps): An Early Visean magnesian plutonism. *Schweizerische Mineralogische Und Petrographische Mitteilungen* 78: 67–85. https://doi.org/10.5169/seals-59275.

Deschamps F, Godard M, Guillot S, Hattori K. 2013. Geochemistry of subduction zone serpentinites: A review. *Lithos* 178: 96–127. https://doi.org/10.1016/j.lithos.2013.05.019.

Duesterhoeff E, Lanari P. 2020. Iterative thermodynamic modelling – Part 1: A theoretical scoring technique and a computer program (Bingo-Antidote). *Journal of Metamorphic Geology* 38(5): 527–551. https://doi.org/10.1111/jmg.12538.

Faure M, Lardeaux J-M, Ledru P. 2009. A review of the pre-Permian geology of the Variscan French Massif Central. *Comptes Rendus Geoscience* 341(2-3): 202–213.

Ferrando S, Lombardo B, Compagnoni R. 2008. Metamorphic history of HP mafic granulites from the Gesso-Stura terrain (Argentera Massif, Western Alps, Italy). *European Journal of Mineralogy* 20 (5): 777–790.

Foley SF, Pertermann M. 2021. Dynamic metasomatism experiments investigating the interaction between migrating potassic melt and garnet peridotite. *Geosciences* 11(10): article 10. https://doi.org/10.3390/geosciences11100432.

Foley SF, Venturelli G, Green DH, Toscani L. 1987. The ultrapotassic rocks: Characteristics, classification, and constraints for petrologic- netic models. *Earth-Science Reviews* 24(2): 81–134. https://doi.org/10.1016/0012-8252(87)90001-8.

Fréville K, Trap P, Faure M, Melleton J, Li X-H, Lin W, et al. 2018. Structural, metamorphic and geochronological insights on the Variscan evolution of the Alpine basement in the Belledonne Massif (France). *Tectonophysics* 726: 14–42.

Fréville K, Trap P, Vanardois J, Melleton J, Faure M, Bruguier O, et al. 2022. Carboniferous-Permian tectono-metamorphic evolution of the Pelvoux Massif (External Crystalline Massif, Western Alps), with discussion on flow kinematics of the Eastern-Variscan Shear Zone. *Bulletin de la Société Géologique de France* 193(1): 13. https://doi.org/10.1051/bsgf/2022008.

Gardien V, Tegyey M, Lardeaux JM, Misseri M, Dufour E. 1990. Crust-mantle relationships in the French Variscan chain: The example of the Southern Monts du Lyonnais unit (Eastern French Massif Central). *Journal of Metamorphic Geology* 8(5): 477–492. https://doi.org/10.1111/j.1525-1314.1990.tb00481.x.

Godard G. 1990. Découverte d'éclogites, de péridotites à spinelle et d'amphibolites à anorthite, spinelle et corindon. *Comptes Rendus de l'Académie Des Sciences. Série 2, Mécanique, Physique, Chimie, Sciences de l'univers, Sciences de la Terre* 310(3): 227–232.

Godard G, Martin S. 2000. Petrogenesis of kelyphites in garnet peridotites: A case study from the Ulten Zone, Italian Alps. *Journal of Geodynamics* 30(1): 117–145. https://doi.org/10.1016/S0264-3707(99)00030-7.

Godard G, Martin S, Prosser G, Kienast JR, Morten L. 1996. Variscan migmatites, eclogites and garnet-peridotites of the Ulten Zone, Eastern Austroalpine system. *Tectonophysics* 259(4): 313–341. https://doi.org/10.1016/0040-1951(95)00145-X.

Grandjean V, Guillot S, Pecher A. 1996. A new record of the LP-HT Late Variscan metamorphism: The Peyre-Arguet unit (Haut- Dauphine). *Comptes Rendus de l'Académie des sciences* 322: 189– 195.

Guillot S, Ménot R. 1999. Nappe stacking and first evidence of Late Variscan extension in the Belledonne Massif (External Crystalline Massifs, French Alps). *Geodinamica Acta* 12: 97–111. https://doi.org/10.1016/S0985-3111(99)80026-6.

Guillot S, Ménot R-P. 2009. Paleozoic evolution of the External Crystalline Massifs of the Western Alps. *Comptes Rendus Geoscience* 341(2-3): 253–265. https://doi.org/10.1016/j.crte.2008.11.010.

Guillot S, Hattori K, Agard P, Schwartz S, Vidal O. 2009. Exhumation processes in oceanic and continental subduction contexts: A review. *Subduction Zone Geodynamics*: 175–205.

Hermann J, Rubatto D. 2009. Accessory phase control on the trace element signature of sediment melts in subduction zones. *Chemical Geology* 265(3): 512–526. https://doi.org/10.1016/j.chemgeo.2009.05.018.

Holland TJB, Powell R. 2011. An improved and extended internally consistent thermodynamic dataset for phases of petrological interest, involving a new equation of state for solids. *Journal of Metamorphic Geology* 29(3): 333–383. https://doi.org/10.1111/j.1525-1314.2010.00923.x.

Ionov DA, Wood BJ. 1992. The oxidation state of subcontinental mantle: Oxygen thermobarometry of mantle xenoliths from Central Asia. *Contributions to Mineralogy and Petrology* 111 (2): 179–193.

Jacob J-B, Guillot S, Rubatto D, Janots E, Melleton J, Faure M. 2021. Carboniferous high pressure metamorphism and deformation in the Belledonne Massif (Western Alps). *Journal of Metamorphic Geology* n/a(n/a). https://doi.org/10.1111/jmg.12600.

Jacob J-B, Janots E, Guillot S, Rubatto D, Fréville K, Melleton J, et al. 2022. HT overprint of HP granulites in the Oisans-Pelvoux Massif: Implications for the dynamics of the Variscan collision in the external Western Alps. *Lithos* 416-417: 106650. https://doi.org/10.1016/j.lithos.2022.106650.

Jagoutz E, Palme H, Baddenhausen H, et al. 1979. The abundances of major, minor and trace elements in the earth's mantle as derived from primitive ultramafic nodules. In: *Lunar and Planetary Science Conference Proceedings*, pp. 2031–2050.

Janoušek V, Holub FV. 2007. The causal link between HP-HT metamorphism and ultrapotassic magmatism in collisional orogens: Case study from the Moldanubian Zone of the Bohemian Massif. *Proceedings of the Geologists' Association* 118(1): 75–86. https://doi.org/10.1016/S0016-7878(07)80049-6.

Janoušek V, Moyen J-F, Martin H, Erban V, Farrow C. 2016. Classical plots. In: Janoušek V, Moyen J-F, Martin H, Erban V, Farrow C, eds. *Geochemical Modelling of Igneous Processes – Principles And Recipes in R Language: Bringing the Power of R to a Geochemical Community* (pp. 27–43). Springer. https://doi.org/10.1007/978-3-662-46792-3_3.

Janoušek V, Bonin B, Collins WJ, Farina F, Bowden P. 2019. Post-Archean granitic rocks: Contrasting petrogenetic processes and tectonic environments. Geological Society of London.
Janoušek V, Hanžl P, Svojtka M, Hora JM, Kochergina YVE, Gadas P, et al. 2020. Ultrapotassic magmatism in the heyday of the Variscan Orogeny: The story of the Třebič Pluton, the largest durbachitic body in the Bohemian Massif. *International Journal of Earth Sciences* 109 (5): 1767–1810. https://doi.org/10.1007/s00531-020-01872-2.

Jarosewich E, Nelen JA, Norberg JA. 1980. Reference samples for electron microprobe analysis. *Geostandards Newsletter* 4(1): 43–47.

Jennings ES, Holland TJB. 2015. A simple thermodynamic model for melting of peridotite in the system NCFMASOcr. *Journal of Petrology* 56(5): 869–892. https://doi.org/10.1093/petrology/egv020.

Jouffray F, Spalla MI, Lardeaux JM, Filippi M, Rebay G, Corsini M, et al. 2020. Variscan eclogites from the Argentera–Mercantour Massif (External Crystalline Massifs, SW Alps): A dismembered cryptic suture zone. *International Journal of Earth Sciences*. https://doi.org/10.1007/s00531-020-01848-2.

Kotková J, O'Brien PJ, Ziemann MA. 2011. Diamond and coesite discovered in Saxony-type granulite: Solution to the Variscan garnet peridotite enigma. *Geology* 39(7): 667–670. https://doi.org/10.1130/G31971.1.

Kubeš M, Leichmann J, Kotková J, Čopjaková R, Holá M, Sláma J. 2022. Diversity of origin and geodynamic evolution of the mantle beneath the Variscan Orogen indicating rapid exhumation within subduction-related mélange (Moldanubian Zone, Bohemian Massif). *Lithos* 422-423: 106726. https://doi.org/10.1016/j.lithos.2022.106726.

Lardeaux JM, Schulmann K, Faure M, Janoušek V, Lexa O, Skrzypek E, et al. 2014. The moldanubian zone in the French Massif Central, Vosges/Schwarzwald and Bohemian Massif revisited: Differences and similarities. *Geological Society, London, Special Publications* 405(1): 7–44.

Le Fort P. 1971. Géologie du Haut-Dauphiné cristallin (Alpes Française): Études pétrologique et structurale de la partie occidentale. PhD Thesis, Université Nancy I.

Lexa O, Schulmann K, Janoušek V, Štípská P, Guy A, Racek M. 2011. Heat sources and trigger mechanisms of exhumation of HP granulites in Variscan orogenic root: Heat sources and exhumation mechanisms. *Journal of Metamorphic Geology* 29(1): 79–102. https://doi.org/10.1111/j.1525-1314.2010.00906.x.

Luais B, Telouk P, Albaréde F. 1997. Precise and accurate neodymium isotopic measurements by plasma-source mass spectrometry. *Geochimica et Cosmochimica Acta* 61(22): 4847–4854.

Lugmair GW, Marti K. 1978. Lunar initial ¹⁴³Nd/¹⁴⁴Nd: Differential evolution of the lunar crust and mantle. *Earth and Planetary Science Letters* 39(3): 349–357.

Maierová P, Schulmann K, Gerya T. 2018. Relamination styles in collisional orogens. *Tectonics* 37(1): 224–250. https://doi.org/10.1002/2017TC004677.

Maierová P, Schulmann K, Štípská P, Gerya T, Lexa O. 2021. Translithospheric diapirism explains the presence of ultra-high pressure rocks in the European Variscides. *Communications Earth & Environment* 2(1): 56. https://doi.org/10.1038/s43247-021-00122-w.

Malaspina N, Hermann J, Scambelluri M, Compagnoni R. 2006. Polyphase inclusions in garnet–orthopyroxenite (Dabie Shan, China) as monitors for metasomatism and fluid-related trace element transfer in subduction zone peridotite. *Earth and Planetary Science Letters* 249(3): 173–187. https://doi.org/10.1016/j.epsl.2006.07.017.

Malatesta C, Crispini L, Federico L, Capponi G, Scambelluri M. 2012. The exhumation of high pressure ophiolites (Voltri Massif, Western Alps): Insights from structural and petrologic data on metagabbro bodies. *Tectonophysics* 568-569: 102–123. https://doi.org/10.1016/j.tecto.2011.08.024.

Martínez Catalán JR, Schulmann K, Ghienne J-F. 2021. The Mid-Variscan Allochthon: Keys from correlation, partial retrodeformation and plate-tectonic reconstruction to unlock the geometry of a non-cylindrical belt. *Earth-Science Reviews* 220: 103700. https://doi.org/10.1016/j.earscirev.2021.103700.

McDonough WF, Sun S-S. 1995. The composition of the Earth. *Chemical Geology* 120: 223–253. https://doi.org/10.1016/0009-2541(94)00140-4.

Medaris G, Wang H, Jelínek E, Mihaljevič M, Jakeš P. 2005. Characteristics and origins of diverse Variscan peridotites in the Gföhl Nappe, Bohemian Massif, Czech Republic. *Lithos* 82(1): 1–23. https://doi.org/10.1016/j.lithos.2004.12.004.

Medaris G, Ackerman L, Jelínek E, Michels ZD, Erban V, Kotková J. 2015. Depletion, cryptic metasomatism, and modal metasomatism (refertilization) of Variscan lithospheric mantle: Evidence from major elements, trace elements, and Sr–Nd–Os isotopes in a Saxothuringian garnet peridotite. *Lithos* 226: 81–97. https://doi.org/10.1016/j.lithos.2014.10.007.

Ménot R-P, Peucat JJ, Scarenzi D, Piboule M. 1988. 496 My age of plagiogranites in the Chamrousse ophiolite complex (external crystalline massifs in the French Alps): Evidence of a Lower Paleozoic oceanization. *Earth and Planetary Science Letters* 88(1-2): 82–92.

Moyen J-F, Laurent O, Chelle-Michou C, Couzinié S, Vanderhaeghe O, Zeh A, et al. 2017. Collision vs. subduction-related magmatism: Two contrasting ways of granite formation and implications for crustal growth. *Lithos* 277: 154–177. https://doi.org/10.1016/j.lithos.2016.09.018.

Nickel KG, Green DH. 1985. Empirical geothermobarometry for garnet peridotites and implications for the nature of the lithosphere, kimberlites and diamonds. *Earth and Planetary Science Letters* 73 (1): 158–170. https://doi.org/10.1016/0012-821X(85)90043-3.

Nimis P, Trommsdorff V. 2001. Revised thermobarometry of Alpe Arami and other Garnet Peridotites from the Central Alps. *Journal of Petrology* 42(1): 103–115. https://doi.org/10.1093/petrology/42.1.103.

Nimis P, Grütter H. 2010. Internally consistent geothermometers for garnet peridotites and pyroxenites. *Contributions to Mineralogy and Petrology* 159(3): 411–427. https://doi.org/10.1007/s00410-009-0455-9.

Nouibat A, Stehly L, Paul A, Schwartz S, Bodin T, Dumont T, et al. 2022. Lithospheric transdimensional ambient-noise tomography of W- Europe: Implications for crustal-scale geometry of the W-Alps. *Geophysical Journal International* 229(2): 862–879. https://doi.org/10.1093/gji/ggab520.

Obata M. 2011. Kelyphite and symplectite: Textural and mineralogical diversities and universality, and a new dynamic view of their structural formation. In: Sharkov E, ed. *New Frontiers in Tectonic Research–General Problems, Sedimentary Basins and Island Arcs*. InTech. https://doi.org/10.5772/20265.

Paquette JL, Ménot RP, Peucat JJ. 1989. REE, Sm–Nd and U–Pb zircon study of eclogites from the Alpine External Massifs (Western Alps): Evidence for crustal contamination. *Earth and Planetary Science Letters* 96(1-2): 181–198.

Paulick H, Bach W, Godard M, De Hoog JCM, Suhr G, Harvey J. 2006. Geochemistry of abyssal peridotites (Mid-Atlantic Ridge, 15.20'N, ODP Leg 209): Implications for fluid/rock interaction in slow spreading environments. *Chemical Geology* 234(3-4): 179–210. https://doi.org/10.1016/j.chemgeo.2006.04.011.

Pecher A. 1970. Étude pétrographique de la partie orientale du massif des Écrins–Pelvoux : Le socle ancien – Alpes françaises. https://tel.archives-ouvertes.fr/tel-00576930.
Perraki M, Faryad SW. 2014. First finding of microdiamond, coesite and other UHP phases in felsic granulites in the Moldanubian Zone: Implications for deep subduction and a revised geodynamic model for Variscan Orogeny in the Bohemian Massif. *Lithos* 202-203: 157–166. https://doi.org/ 10.1016/j.lithos.2014.05.025.

Pin C, Carme F. 1987. A Sm-Nd isotopic study of 500 Ma old oceanic crust in the Variscan Belt of Western Europe: The Chamrousse ophiolite complex, Western Alps (France). *Contributions to Mineralogy and Petrology* 96(3): 406–413.

Pin C, Zalduegui JS. 1997. Sequential separation of light rare-earth elements, thorium and uranium by miniaturized extraction chromatography: Application to isotopic analyses of silicate rocks. *Analytica Chimica Acta* 339(1-2): 79-89.

Pin C, Briot D, Bassin C, Poitrasson F. 1994. Concomitant separation of strontium and samarium-neodymium for isotopic analysis in silicate samples, based on specific extraction chromatography. *Analytica Chimica Acta* 298(2): 209–217.

Rubatto D, Ferrando S, Compagnoni R, Lombardo B. 2010. Carboniferous high-pressure metamorphism of Ordovician protoliths in the Argentera Massif (Italy), Southern European Variscan Belt. *Lithos* 116(1-2): 65–76.

Ruiz M, Schaltegger U, Gaynor SP, Chiaradia M, Abrecht J, Gisler C, et al. 2022. Reassessing the intrusive tempo and magma genesis of the late Variscan Aar batholith: U–Pb geochronology, trace element and initial Hf isotope composition of zircon. *Swiss Journal of Geosciences* 115(1): 20. https://doi.org/10.1186/s00015-022-00420-1.

Sapienza GT, Scambelluri M, Braga R. 2009. Dolomite-bearing orogenic garnet peridotites witness fluid-mediated carbon recyc-
ling in a mantle wedge (Ulten Zone, Eastern Alps, Italy).
*Contributions to Mineralogy and Petrology* 158(3): 401–420.
https://doi.org/10.1007/s00410-009-0389-2.

Scambelluri M, Hermann J, Morten L, Rampone E. 2006. Melt-versus fluid-induced metasomatism in spinel to garnet wedge peridotites (Ulten Zone, Eastern Italian Alps): Clues from trace element and Li abundances. *Contributions to Mineralogy and Petrology* 151(4): 372–394. https://doi.org/10.1007/s00410-006-0064-9.

Schulmann K, Konopásek J, Janoušek V, Lexa O, Lardeaux J-M, Edel J-B, et al. 2009. An Andean type Palaeozoic convergence in the Bohemian Massif. *Comptes Rendus Geoscience* 341: 266–286. https://doi.org/10.1016/j.crte.2008.12.006.

Schulmann K, Catalán J-RM, Lardeaux J-M, Janoušek V, Oggiano G. 2014. The Variscan orogeny: Extent, timescale and the formation of the European crust. *Geological Society, London, Special Publications* 405(1): 1–6. https://doi.org/10.1144/SP405.15.

Schulmann K, Lexa O, Janoušek V, Lardeaux J-M, Edel J-B. 2014. Anatomy of a diffuse cryptic suture zone: An example from the Bohemian Massif, European Variscides. *Geology* 42(4): 275–278. https://doi.org/10.1130/G35290.1.

Skrzypek E, Štípská P, Cocherie A. 2012. The origin of zircon and the significance of U–Pb ages in high-grade metamorphic rocks: A case study from the Variscan orogenic root (Vosges Mountains, NE France). *Contributions to Mineralogy and Petrology* 164(6): 935–957. https://doi.org/10.1007/s00410-012-0781-1.

Sobolev AV, Hofmann AW, Kuzmin DV, Yaxley GM, Arndt NT, Chung S-L, et al. 2007. The amount of recycled crust in sources of mantle-derived melts. *Science* 316(5823): 412–417. https://doi.org/10.1126/science.1138113.

Soder CG, Romer RL. 2018. Post-collisional Potassic–Ultrapotassic Magmatism of the Variscan Orogen: Implications for Mantle Metasomatism during Continental Subduction. *Journal of Petrology* 59(6): 1007–1034. https://doi.org/10.1093/petrology/egy053.

Steiger RH, Jäger E. 1977. Subcommission on geochronology: Convention on the use of decay constants in geo- and cosmochronology. *Earth and Planetary Science Letters* 36(3): 359–362. https://doi.org/10.1016/0012-821X(77)90060-7.

Tanaka T, Togashi S, Kamioka H, Amakawa H, Kagami H, Hamamoto T, et al. 2000. JNdi-1: A neodymium isotopic reference in consistency with LaJolla neodymium. *Chemical Geology* 168 (3-4): 279-281.

Taylor WR. 1998. An experimental test of some geothermometer and geobaro-meter formulations for upper mantle peridotites with application to the ther-mobarometry of fertile lherzolite and garnet websterite. *Neues Jahrbuch Für Mineralogie–Abhandlungen*: 381–408. https://doi.org/10.1127/njma/172/1998/381.

Tumiati S, Thöni M, Nimis P, Martin S, Mair V. 2003. Mantle-crust interactions during Variscan subduction in the Eastern Alps (Nonsberg-Ulten Zone): Geochronology and new petrological constraints. *Earth and Planetary Science Letters* 210(3-4): 509–526.

Vanardois J, Roger F, Trap P, Goncalves P, Lanari P, Paquette J-L, et al. 2022. Exhumation of deep continental crust in a transpressive regime: The example of Variscan eclogites from the Aiguilles-Rouges Massif (Western Alps). *Journal of Metamorphic Geology*.

Vanderhaeghe O, Laurent O, Gardien V, Moyen J-F, Gébelin A, Chelle-Michou C, et al. 2020. Flow of partially molten crust controlling construction, growth and collapse of the Variscan orogenic belt: The geologic record of the French Massif Central. *Bulletin de la Société Géologique de France* 191(25). https://doi.org/10.1051/bsgf/2020013.

Von Raumer JF, Bussy F. 2004. Mont Blanc and Aiguilles Rouges geology of their polymetamorphic basement (external massifs, Westerns Alps, France-Switzerland). *Mémoires de Géologie (Lausanne)* 42: 1–210.

Von Raumer JF, Stampfli GM. 2008. The birth of the Rheic Ocean–Early Palaeozoic subsidence patterns and subsequent tectonic plate scenarios. *Tectonophysics* 461(1-4): 9–20.

Von Raumer JF, Finger F, Veselá P, Stampfli GM. 2014. Durbachites–Vaugnerites–A geodynamic marker in the Central European Variscan orogen. *Terra Nova* 26(2): 85–95.

Wang J-M, Lanari P, Wu F-Y, Zhang J-J, Khanal GP, Yang L. 2021. First evidence of eclogites overprinted by ultrahigh temperature metamorphism in Everest East, Himalaya: Implications for collisional tectonics on early Earth. *Earth and Planetary Science Letters* 558: 116760. https://doi.org/10.1016/j.epsl.2021.116760.

Zheng Y-F. 2019. Subduction zone geochemistry. *Geoscience Frontiers* 10(4): 1223–1254. https://doi.org/10.1016/j.gsf.2019.02.003.

Ziberna L, Klemme S, Nimis P. 2013. Garnet and spinel in fertile and depleted mantle: Insights from thermodynamic modelling. *Contributions to Mineralogy and Petrology* 166(2): 411–421. https://doi.org/10.1007/s00410-013-0882-5.

Cite this article as: Jacob J-B, Janots E, Cordier C, Guillot S. 2023. Discovery of Variscan orogenic peridotites in the Pelvoux Massif (Western Alps, France), BSGF - *Earth Sciences Bulletin* 194: 2.