Xiaoyan Gu

Jannick Ingrin

Etienne Deloule

Lydéric France

Qunke Xia

Mantle metasomatism in the sub-continental lithospheric mantle beneath the South French Massif Central: Constraints from trace elements, Li and H in peridotite minerals

Keywords: French Massif Central, Mantle metasomatism, Peridotite xenoliths, Li concentrations, δ 7 Li, Water

Mantle metasomatism by percolating melts/fluids can significantly modify the geochemical and mineralogical compositions of the sub-continental lithospheric mantle (SCLM). We present a detailed study of these modifications in Li concentrations, δ 7 Li and water contents in mantle minerals from a suite of peridotite xenoliths entrained by a Cenozoic strombolian volcano in the south of the French Massif Central (FMC). The trace element compositions in clinopyroxenes and amphiboles were first investigated to track metasomatism. The wide spectrum of trace element distribution in clinopyroxenes (e.g., (La/Yb) N from 0.25 to 22.21; Ti/Eu ratios from 453 to 4892) suggests that the SCLM have undergone metasomatism by carbonatitic melts/fluids or melts/fluids related to subducted materials. But two amphibole-bearing samples exhibit equilibrated trace element partitioning between clinopyroxenes and coexisting amphiboles, with depletion of LREE ((La/Yb) N =0.26 and 0.30), indicating that amphiboles have grown in another modally metasomatic process predating to the cryptic metasomatism accounting for LREE enrichment and HFSE negative anomalies in other samples. Li concentrations in minerals are similar to those in the normal mantle, with inter-mineral Li partitioning nearly equilibrated and intragranular Li distribution nearly homogeneous. However, the negative δ 7 Li values of clinopyroxenes and orthopyroxenes in some samples (as low as -8.8 ‰ in clinopyroxene of sample MC38) can be attributed to exchange with a small-volume melt of moderate Li concentration and light Li isotopic compositions, originally associated with a recycled component. The preservation of currently observed inter-mineral large Li isotopic difference indicates that the percolation of the melt should have occurred shortly before the entrainment of peridotite xenoliths by the host magma. Water contents in minerals vary from 41 to 428 ppm in clinopyroxenes, from 28 to 152 ppm in orthopyroxenes, and their roughly negative co-variation with Fo contents in coexisting olivines implies that partial melting has mainly controlled their variations among most samples. The different behaviors of OH in different metasomatized samples with LREE enrichment point to the involvement of metasomatic agents of different origins. The aqueous agent responsible for generation of amphiboles in two samples has not lead to a notable increase of water content.

Introduction

Peridotite xenoliths, commonly entrained by alkali basalts or kimberlites, can provide direct informations of the sub-continental lithospheric mantle (SCLM; e.g., Pearson et al., 2014). Studies on them have demonstrated that the SCLM, depleted by partial melting in various degrees, frequently experiences modal or cryptical metasomatic modifications by melts/fluids (known as mantle metasomatism; e.g., Menzies and Hawkesworth, 1987;[START_REF] O'reilly | Mantle Metasomatism. Metasomatism and the Chemical Transformation of Rock[END_REF].The latter processes are the source of significant mineralogical and geochemical variations of the SCLM (e.g., [START_REF] Frey | The mineralogy, geochemistry and origin of Iherzolite inclusions in Victorian basanites[END_REF][START_REF] Frey | Ultramafic inclusions from San Carlos, Arizona: petrologic and geochemical data bearing on their petrogenesis[END_REF][START_REF] O'reilly | Mantle metasomatism beneath western victoria, australia: i. metasomatic processes in cr-diopside lherzolites[END_REF][START_REF] Rudnick | Carbonatite metasomatism in the northern Tanzanian mantle: petrographic and geochemical characteristics[END_REF][START_REF] Rudnick | Petrology and geochemistry of spinel peridotite xenoliths from Hannuoba and Qixia, North China craton[END_REF][START_REF] Coltorti | Carbonatite metasomatism of the oceanic upper mantle: evidence from clinopyroxenes and glasses in ultramafic xenoliths of Grande Comore, Indian Ocean[END_REF][START_REF] Ionov | Mechanisms and sources of mantle metasomatism: major and trace element compositions of peridotite xenoliths from Spitsbergen in the context of numerical modelling[END_REF].

The fact that hydrogen can be incorporated into crystal lattices of mantle nominally anhydrous minerals (NAMs) has been largely recognized (e.g., [START_REF] Bell | Water in Earth's mantle: the role of nominally anhydrous minerals[END_REF][START_REF] Ingrin | Hydrogen in nominally anhydrous upper-mantle minerals concentration levels and implications[END_REF][START_REF] Hirschmann | Water, melting, and the deep Earth H 2 O cycle[END_REF][START_REF] Demouchy | Distribution and transport of hydrogen in the lithospheric mantle: A review[END_REF]. As a trace element in NAMs, hydrogen also behaves incompatibly during partial melting and fractional crystallization, with the incompatibility usually compared to Ce [START_REF] Michael | Regionally distinctive sources of depleted MORB: evidence from trace elements and H 2 O[END_REF]. Meanwhile, since water is a common component of metasomatic agents at mantle depths, hydrogen contents in mantle NAMs (conventionally measured in ppm wt. H 2 O) is potentially an effective tracer of mantle metasomatism [START_REF] Peslier | Metasomatic control of water contents in the Kaapvaal cratonic mantle[END_REF][START_REF] Doucet | High water contents in the Siberian cratonic mantle linked to metasomatism: An FTIR study of Udachnaya peridotite xenoliths[END_REF][START_REF] Demouchy | Characterization of hydration in the mantle lithosphere: Peridotite xenoliths from the Ontong Java Plateau as an example[END_REF][START_REF] Denis | Water content and hydrogen behaviour during metasomatism in the uppermost mantle beneath Ray Pic volcano (Massif Central, France)[END_REF].

Li incorporates into silicate mantle minerals by substituting Mg and Fe [START_REF] Seitz | The distribution of lithium in peridotitic and pyroxenitic mantle lithologies-an indicator of magmatic and metasomatic processes[END_REF][START_REF] Zhang | Coupled (Li + , Al 3+ ) substitutions in hydrous forsterite[END_REF] and it has an intermediate incompatibility similar to Yb during partial melting processes [START_REF] Ryan | The systematics of lithium abundances in young volcanic rocks[END_REF][START_REF] Brenan | Behaviour of boron, beryllium, and lithium during melting and crystallization: constraints from mineral-melt partitioning experiments[END_REF][START_REF] Ottolini | New experimental determination of Li and B partition coefficients during upper mantle partial melting[END_REF]. The large mass difference between 6 Li and 7 Li (-~16%) results in significant isotopic fractionation during various geochemical processes and produces large Li isotopic variations (see Tomascak et al. (2016) for a review).

Melt/fluid-rock interactions strongly affect both Li contents and Li isotopic compositions of mantle minerals (e.g., [START_REF] Seitz | The distribution of lithium in peridotitic and pyroxenitic mantle lithologies-an indicator of magmatic and metasomatic processes[END_REF][START_REF] Seitz | Lithium isotopic signatures of peridotite xenoliths and isotopic fractionation at high temperature between olivine and pyroxenes[END_REF][START_REF] Woodland | Varying behaviour of Li in metasomatised spinel peridotite xenoliths from western Victoria, Australia[END_REF][START_REF] Tang | Lithium isotopic systematics of peridotite xenoliths from Hannuoba, North China Craton: implications for melt-rock interaction in the considerably thinned lithospheric mantle[END_REF][START_REF] Tang | Multistage melt/fluid-peridotite interactions in the refertilized lithospheric mantle beneath the North China Craton: constraints from the Li-Sr-Nd isotopic disequilibrium between minerals of peridotite xenoliths[END_REF]Albach and Rudnick, 2009;Zhang et al., 2009;[START_REF] Su | Extremely high Li and low δ 7 Li signatures in the lithospheric mantle[END_REF][START_REF] Tang | Abnormal lithium isotope composition from the ancient lithospheric mantle beneath the North China Craton[END_REF][START_REF] Gu | Multi-stage metasomatism revealed by trace element and Li isotope distributions in minerals of peridotite xenoliths from Allègre volcano (French Massif Central)[END_REF].

The basement of the French Massif Central (FMC) was mainly formed during the Variscan orogeny (e.g., [START_REF] Matte | Tectonics and plate tectonics model for the Variscan belt of Europe[END_REF]. Based on petrological and geochemical observations on large quantities of peridotite xenoliths, it was proposed that the underlying SCLM is composed of two distinct domains to the north and south of the boundary at ~45°30'N (Fig. 1; Lenoire et al., 2000). The northern domain is more refractory than the southern domain. The Sr-Nd isotopic compositions of the former overlaps the compositional field of the European Asthenospheric Reservoir and those of the latter are more similar to MORB [START_REF] Granet | Imaging a mantle plume beneath the French Massif Central[END_REF][START_REF] Lenoir | Contrasting lithospheric mantle domains beneath the Massif Central (France) revealed by geochemistry of peridotite xenoliths[END_REF][START_REF] Downes | Mantle domains in the lithosphere beneath the French Massif Central: trace element and isotopic evidence from mantle clinopyroxenes[END_REF]. Several studies have documented that remarkable negative anomalies of HFSE were mainly observed in xenoliths from the northern domain but less common in xenoliths from the southern domain. Thus, the northern and southern domains of the FMC were considered to be metasomatized by different types of agents, fluids/carbonatitic melts and silicate melts, respectively [START_REF] Lenoir | Contrasting lithospheric mantle domains beneath the Massif Central (France) revealed by geochemistry of peridotite xenoliths[END_REF][START_REF] Downes | Mantle domains in the lithosphere beneath the French Massif Central: trace element and isotopic evidence from mantle clinopyroxenes[END_REF]Witting et al., 2007). Some other studies pointed out the imprints of melts/fluids likely related to recycled materials during Variscan subduction in some peridotites and pyroxenites [START_REF] Downes | Textural, isotopic and REE variations in spinel peridotite xenoliths, Massif Central, France[END_REF][START_REF] Deloule | Hydrogen isotope heterogeneities in the mantle from ion probe analysis of amphiboles from ultramafic rocks[END_REF][START_REF] Touron | Characterization of the metasomatic agent in mantle xenoliths from Deves, Massif Central (France) using coupled in situ trace-element and O, Sr and Nd isotopic compositions[END_REF][START_REF] Yoshikawa | Geochemical and Sr-Nd isotopic characteristics and pressure-temperature estimates of mantle xenoliths from the French Massif Central: evidence for melting and multiple metasomatism by silicate-rich carbonatite and asthenospheric melts[END_REF][START_REF] Gu | Multi-stage metasomatism revealed by trace element and Li isotope distributions in minerals of peridotite xenoliths from Allègre volcano (French Massif Central)[END_REF].

In this study, the major and trace element compositions, water contents, and Li concentrations and isotopic compositions of minerals in peridotite xenoliths transported by a Cenozoic volcano in the FMC were analyzed to deeply investigate the evolution of the SCLM beneath the FMC, and the behaviors of volatiles like H and Li during mantle metasomatism.

Geological background

Extensive magmatic activities occurred across the FMC since the early Cenozoic.

They have been sub-divided into three magmatic phases: pre-rift volcanism due to lithospheric bending, rift-related volcanism due to crustal thinning, and a later (15 Ma ~ 6.9 ka ago) magmatic event resulting from asthenospheric upwelling (Michon and Merle, 2001). They spatially distribute in a few regions (Fig. 1). Large amounts of mantle xenoliths have been transported to the surface by these volcanoes (Lenoire et al., 2000). Low-velocity anomalies beneath the FMC have been observed by seismic tomography [START_REF] Granet | Imaging a mantle plume beneath the French Massif Central[END_REF][START_REF] Goes | A lower mantle source for central European volcanism[END_REF][START_REF] Fichtner | Crust and upper mantle of the western Mediterranean-Constraints from full-waveform inversion[END_REF].

They have been attributed to either a mantle plume [START_REF] Granet | Imaging a mantle plume beneath the French Massif Central[END_REF][START_REF] Goes | A lower mantle source for central European volcanism[END_REF], or a shallow thermal anomaly related to the Alpine Orogeny [START_REF] Fichtner | Crust and upper mantle of the western Mediterranean-Constraints from full-waveform inversion[END_REF].

The sampling locality was the Razas Grand quarry in Mont Coupet, situated in the Devès in the southern FMC domain (Fig. 1). The volcano erupted in strombolian mode 2 Ma ago.

Analytical methods

Sample preparations

More than twenty peridotite xenoliths with diameters of 3~12 cm wrapped in volcanic bombs were collected in the Razas Grand quarry. The outer basaltic shells were sawed off and the cores were cut into billets for polishing thin sections. Then, detailed petrographic observations were performed on thin sections under microscope.

According to modal compositions, textures, the colors of spinel (Spl), eight representative fresh samples, free of visible intrusive veins and of any weathering feature, were selected for electron micro-probe analysis. After that analysis, a subset of six samples, covering the total major element composition range, was used for trace elements, Li, Li isotopic composition and water content measurements.

Electron micro-probe analysis

Mineral major element compositions were measured on thin sections using a Cameca SX100 electron microprobe at the Service CAMPARIS, in the Pierre and Marie Curie University (France). The operating conditions were set at an accelerating voltage of 15 kV and a beam current of 10 nA with a focused beam. Na and K were analyzed first with counting times of 10s on peak and 5s for backgrounds to minimize errors due to element migration. For other elements, counting times were 20s on peak and 10s for backgrounds. A series of natural and synthetic minerals was used as standards. At least three grains were analyzed for each mineral phase in a single thin section; at least three points were analyzed on the core and rims of each crystal.

Laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS)

Trace element compositions in clinopyroxene (Cpx) and amphibole (Amp) were determined using LA-ICP-MS at GeoRessources Laboratory, University of Lorraine (France). Double-polished thick sections (~0.15 mm thick) or thin sections were ablated in-situ using a nanosecond excimer laser (GEOLAS Pro;193 nm wavelength).

The spot size was 44 µm, the ablation frequency 5 Hz and the energy density 10 J/cm 2 per pulse. The ablation products were transported in a helium flow, mixed with argon gas, and then analyzed with an Agilent 7500 ICP-MS. A complete analysis comprises 30s for background acquisition, and 50s for sample acquisition. The raw intensities of ions were recorded as a function of time. SiO 2 contents obtained from electron micro-probe analysis were used as an internal standard; NIST 612 and 614, analyzed at the beginning and end of the analytical session, were used as external standards.

The concentrations of some minor elements in olivine (Ol), such as Al, Ca, Ti, Cr and Ni, were also analyzed following the same procedure.

Secondary ion mass spectrometry (SIMS)

Li concentrations and isotopic compositions were measured with a Cameca IMS1270 ion microprobe at the Centre de Recherches Pétrographiques et Géochimiques (CRPG) in France; thin sections were gold-coated before analysis. At least three points (the number of analyses depending on grain size) were measured from core to rim on each grain, in areas devoid of cracks and inclusions.

To sputter the secondary positive ions ( 7 Li + and 6 Li + ), a primary 16 O -beam was used at intensity ranging from 8 to 15nA, with spot diameters of 20-40 µm, slightly varying with the different instrumental settings in individual analytical sessions. The 10kV secondary positive ions were counted without energy offset at a mass resolution of 1100 (M/ΔM) to separate the 6 LiH from 7 Li, in mono-collection mode. Counting rates for 7 Li ranged from 1×10 4 to 1×10 5 cps, depending on the Li content and the primary beam intensity. After a 120s pre-sputtering, twenty-five to forty cycles were accumulated with counting times of 4s, 12s and 4s for background at mass 5.7, 6 Li and 7 Li, respectively. Prior to every analytical session, the deadtime of the counting system was calculated according to [START_REF] Deloule | Instrumental limitations for isotope measurements with a Caméca® ims-3f ion microprobe: Example of H, B, S and Sr[END_REF]. The Li isotopic ratios are expressed as δ 7 Li (δ 7 Li = [( 7 Li/ 6 Li) sample /( 7 Li/ 6 Li) L-SVEC -1]×1000, with /( 7 Li/ 6 Li) L-SVEC = 12.0192 [START_REF] Flesch | A secondary isotopic standard for 6 Li/7 Li determinations[END_REF]). A series of mineral standards established by [START_REF] Su | Potential Orthopyroxene, Clinopyroxene and Olivine Reference Materials for In Situ Lithium Isotope Determination[END_REF] was applied to calibrate the instrumental mass fractionation, defined as Δ i = δ 7 Li SIMS -δ 7 Li known [START_REF] Decitre | Behavior of Li and its isotopes during serpentinization of oceanic peridotites[END_REF]. These reference minerals are from mantle xenoliths, with Mg# comparable to the equivalent mineral phases in our samples (Supplementary table 1). Throughout the whole analytical session, the standards were analyzed to monitor the fluctuation of Δ i with temporal variations of instrumental parameters, including the primary beam intensity [START_REF] Fitzsimons | SIMS stable isotope measurement: counting statistics and analytical precision[END_REF]. Unlike previous studies using an IMS 3f ion probe (e.g., [START_REF] Decitre | Behavior of Li and its isotopes during serpentinization of oceanic peridotites[END_REF][START_REF] Beck | Li isotopic variations in single pyroxenes from the Northwest Africa 480 shergottite (NWA 480): a record of degassing of Martian magmas?[END_REF], reference pyroxenes and Ol do not produce the same Δ i ; pyroxenes and Ol have co-variations of -0.39 and -0.52 δ unit per Mg#, respectively (Fig. 2a).Thus, due to the narrow and similar range of Mg# in our samples relative to the standards, no large inaccuracy (±1‰) is introduced when applying an averaged Δ i of individual standard minerals to calculate the δ 7 Li values of the corresponding sample minerals. The absence of significant compositional variation across the analyzed grains, as well as the limited compositional range of the samples, dictate that variations in our results are not the consequence of matrix effects.

The lithium ionization yield ([Li ion counts / primary beam intensity] / Li concentration)) was determined on the reference minerals (Fig. 2b), and Li concentrations of samples were calculated using average ionization yields of the standards. The estimated precision of this procedure is better than 20%. Li concentrations and isotopic composition in Amp in samples MC36 and MC53 were also measured. Because of the similar crystalline structure of Amp and Cpx, the instrumental fractionation for measuring Li isotopic compositions of Amp was calibrated with Cpx standards.

Fourier Transform Infrared (FTIR) Spectrometer

Spectrum collection

Infrared spectra were acquired on double polished wafers of 0.1 mm to 0.2 mm in thickness, using a Bruker Hyperion 3000 FTIR-microscope attached to a Bruker Vertex 70 spectrometer equipped with a liquid nitrogen cooled MCT detector and a KBr beam splitter in the laboratory of Infrared and Raman spectrochemistry (LASIR, Université Lille 1). The entire instrument was continuously flushed with dry air during analysis. More than twenty non-oriented grains of each mineral phase (Ol, Opx and Cpx) in a single sample were analyzed under unpolarized light in the transmission mode, with the exception of sample MC38, which is too depleted to contain enough available Cpx grains for measurement. The spectra were collected close to the cores of the grains, through optically clean pathway, void of inclusion and crack. 256 scans were accumulated at a resolution of 4 cm -1 . A squared aperture (30×30 to 100×100 µm 2 ) was used and adjusted depending on the mineral grain size and quality. To check the intragranular water distribution homogeneity, several large Ol and Opx grains were chosen to make profile measurements with a step size of 60 µm. FTIR analyses by polarized light were also performed on Ol grains, the unpolarized spectra of which display hydroxyl absorption bands. The orientation of Ol crystals are identified through their Si-O overtone vibrations (1200-2200 cm -1 ; [START_REF] Lemaire | The effect of silica activity on the incorporation mechanisms of water in synthetic forsterite: a polarized infrared spectroscopic study[END_REF]. The crystals with at least one principal axis in the plane of the thin sections were selected and the polarized spectrum was acquired for this particular orientation.

The polarized light was obtained using a ZnSe grid polarizer. In every sample, two or three spectra were obtained on different grains for light parallel to each principle axis.

Calculation from infrared spectra to water contents

Water contents in minerals have been calculated from FTIR spectra based on the modified Beer-Lambert law:

2 =A/εtγ (1)
in which, 2 is OH concentrations (H 2 O in wt. ppm), A is the integral absorbance which is represented here by the integrated area beneath the hydroxyl absorption bands (cm -2 ), ε is the molar absorption coefficient (ppm -1 cm -2 ), t is the thickness of the sample (cm). The protocol followed to calculate A is to integrate the area beneath the OH absorption bands of the IR spectra averaged from more than 10 grains; this approach is reliable when linear unpolarized absorbance is below 0.3 (deviation < 10%; [START_REF] Kovács | Quantitative absorbance spectroscopy with unpolarized light: part II. Experimental evaluation and development of a protocol for quantitative analysis of mineral IR spectra[END_REF][START_REF] Withers | On the use of unpolarized infrared spectroscopy for quantitative analysis of absorbing species in birefringent crystals[END_REF]. Given that OH spectra in Ol are remarkably more anisotropic and have much lower absorbance than those in coexisting pyroxenes, FTIR measurements under polarized light on Ol grains were also made to compare and confirm the reliability of results from unpolarized absorbance. The bands located at wavenumbers between 3700 and 2800 cm -1 for Opx, at wavenumbers between 3700 and 3000 cm -1 for Cpx were considered when integrating. The absorption coefficients determined by [START_REF] Bell | Quantitative analysis of trace OH in garnet and pyroxenes[END_REF] are the most generally used for pyroxenes and it is the ones we used here: 14.84 ppm -1 cm -2 for Opx and 7.09 ppm -1 cm -2 for Cpx. In our samples, OH in Ol from mantle xenoliths is characterized by bands related to three types of defects located in three different wavenumber parts of OH absorption spectra: the group I located at wavenumbers between the two bands related to the Ti-H defectes at 3572 and 3525 cm -1 , the OH bands related to substitution of trivalent cations in the octahedral site in the range of 3300-3400 cm -1 , and the OH bands associated with octahedral vacancies in the range of 3100-3300 cm -1 [START_REF] Berry | Fingerprinting the water site in mantle olivine[END_REF][START_REF] Berry | The infrared signature of water associated with trivalent cations in olivine[END_REF][START_REF] Kovács | Site-specific infrared OH absorption coefficients for water substitution into olivine[END_REF]Blanchard et al., 2016). The spectra of Ol to estimate the absorption coefficient in [START_REF] Bell | Hydroxide in olivine: a quantitative determination of the absolute amount and calibration of the IR spectrum[END_REF] solely have bands related to Ti-H defects. Thus, when calculating the water contents in Ol, the bands will be integrated separately and equation 1 was applied along with the absorption coefficient (0.188 ppm -1 cm -2 ) from [START_REF] Bell | Hydroxide in olivine: a quantitative determination of the absolute amount and calibration of the IR spectrum[END_REF] for bands at 3572 and 3525 cm -1 , whereas for the bands at lower wavenumbers associated witg trivalent substitution and octahedral vacancies the calculation method established by [START_REF] Libowitzky | An IR absorption calibration for water in minerals[END_REF] was used:

2 =18000×3A/ε i tD (2) ε i =246.6×(3753-ν) (3)
where ε i is the integrated molar absorption coefficient (Lcm -2 mol -1 ) varying on the hydrogen bond strength (i.e., the wavenumbers (ν in equation 3; cm -1 )), D is the density of concerned mineral (gcm -3 ; here 3.3 gcm -3 for Ol). Here, for unpolarized average spectra the integrated polarized absorbance (requested in the initial formula of [START_REF] Libowitzky | An IR absorption calibration for water in minerals[END_REF]) is replaced by three times the average integrated absorbance (3×A). The thickness of samples was measured using a digital micrometer.

The baseline was fitted by a 3 rd order polynomial using KaleidaGraph software.

As for the polarized spectra, the integrated absorbances when E is parallel to three principle axes (x ([100]), y ([010]) and z ([001]); Fig. 4) were added together to take place of A in equation 1 and 3A in equation 2. The estimated errors from unpolarized measurements are less than 30% and those from polarized measurements are less than 10%.

Petrological description

According to the petrographical observations, the eight samples can be divided into two groups based on the Cpx modal proportions, spinel lherzolite (Cpx up to 15%) and spinel harzburgite (Cpx < 5%; Supplementary table 2). Two samples (MC36 and MC53) contain amphibole (Amp) in small modal proportions (less than 1%).

The samples display two types of texture, protogranular and porphyroclastic textures. The protogranular texture, more common in our samples, is regarded as an original texture that evolved (through plastic flow and recrystallization) to porphyroclastic [START_REF] Mercier | Textures and fabrics of upper-mantle peridotites as illustrated by xenoliths from basalts[END_REF][START_REF] Downes | Textural, isotopic and REE variations in spinel peridotite xenoliths, Massif Central, France[END_REF]. In protogranular samples, Ol and Opx are coarse-grained and have similar grain sizes

(1.5-2.5mm), whereas Cpx scattered interstitially with similar-sized or larger coexisting Ol and Opx (Fig. 3a). In samples with the porphyroclastic texture (Fig. 3b), large Ol and Opx porphyroclasts (up to 4mm across) are embedded in smaller Ol, Opx, Cpx and Spl granuloblasts (<1mm; Fig. 3b). Some Ol grains have kink-banded texture (Fig. 3c). Many Opx grain boundaries are cribriform-shaped (Fig. 3d) and fine Cpx exsolutions have been observed in large Opx grains under polarized light (Fig. 3b).

Amp grains in samples MC36 and MC53 are small with diameter of <1mm distributed in interstitial positions, usually adjacent to Cpx and Spl.

Results

Mineral major element compositions

The major element compositions of minerals in Mont Coupet samples are reported in Table 1. All the values are the averages from all the measured points on the same phase in a single sample. No significant core-rim variation was observed in any sample.

Ol forsterite (Fo) content (or Mg# = 100×molarMg/(Mg+Fe)) in samples from Mont Coupet varies from 89.6 to 91.3, falling in the range of Fo content in off-craton peridotites worldwide [START_REF] Rudnick | Petrology and geochemistry of spinel peridotite xenoliths from Hannuoba and Qixia, North China craton[END_REF]. No large inter-grain variation has been observed in individual samples (e.g., the standard deviation of MgO contents within every sample is less than 0.48). The Al, Ti, Cr and Ni concentrations in Ol have the range of 13.5-92.5 ppm, 3.7-14.0 ppm, 3.7-46.8 ppm and 3085-3267 ppm, respectively (Supplementary table 3). Most of these values fall in the mantle peridotite Ol ranges [START_REF] De Hoog | Trace-element geochemistry of mantle olivine and application to mantle petrogenesis and geothermobarometry[END_REF][START_REF] Foley | Minor and trace elements in olivines as probes into early igneous and mantle melting processes[END_REF].

Opx Mg# varies from 89.8 to 91.5, and co-existing Opx and Ol have similarMg#(Table 1), which suggests a chemical equilibrium between Ol and Opx after partial melting in individual samples [START_REF] Lee | Compositionally stratified cratonic lithosphere: petrology and geochemistry of peridotite xenoliths from the Labait volcano, Tanzania[END_REF]Ottoloni et al., 2004). Al 2 O 3 and CaO contents range from 1.72 to 3.91 wt. % and from 0.31 to 0.50 wt. %, respectively. TiO 2 contents (0.01 to 0.12 wt. %) are slightly lower than in other peridotites worldwide, but similar low TiO 2 contents have been reported in peridotites metasomatized by subduction-related melts (e.g., [START_REF] Zanetti | The Finero phlogopite-peridotite massif: an example of subduction-related metasomatism[END_REF][START_REF] Grégoire | Nature and evolution of the lithospheric mantle beneath the passive margin of East Oman: evidence from mantle xenoliths sampled by Cenozoic alkaline lavas[END_REF][START_REF] Woo | Silica-and LREE-enriched spinel peridotite xenoliths from the Quaternary intraplate alkali basalt, Jeju Island, South Korea: Old subarc fragments?[END_REF].

Cpx are Cr-diopsides with Cr 2 O 3 contents ranging from 0.57 to 1.34 wt. %. Mg# of Cpx (90.2 to 93.5) is higher and shows larger variations than Fo content of coexisting Ol and Mg# of coexisting Opx (Table 3). Al 2 O 3 and CaO contents in Cpx range from 1.88 to 6.74 wt. % and from 20.72 to 23.21wt. %, respectively.

Spl Cr# (=100×molarCr/(Cr+Al)) varies from 8.3 to 45.2, which falls in the Cr# range of Spl from subcontinental peridotite xenoliths [START_REF] Arai | Characterization of spinel peridotites by olivine-spinel compositional relationships: review and interpretation[END_REF]. The ratio of Fe 3+ /ΣFe, as an indicator of the redox state, has been calculated from 0.07 to 0.24 in studied samples. Some of these values are much lower than the range of the Spl in the FMC (0.16 to 0.46, [START_REF] Uenver-Thiele | Oxidation state of the lithospheric mantle below the Massif Central, France[END_REF].

Eight Amp grains (four for each sample) were analyzed using electron microprobe.

According to [START_REF] Leake | Nomenclature of amphiboles; report of the subcommittee on amphiboles of the International Mineralogical Association, Commission on New Minerals and Mineral Names[END_REF], Amp in both samples are pargasite. Most major elements have similar contents in the two samples, with the exception of TiO 2 , 2.76 and 3.46 wt. % in MC36 and MC53, respectively (Table 3). Compared with compositions of Amp documented in previous studies [START_REF] Touron | Characterization of the metasomatic agent in mantle xenoliths from Deves, Massif Central (France) using coupled in situ trace-element and O, Sr and Nd isotopic compositions[END_REF]Uenve-Thiele et al., 2014), no obvious distinction occurs in samples from Mont Coupet.

Trace element compositions in Cpx and Amp

Cpx and Amp trace element compositions are reported in Table 2. In each sample, more than 7 Cpx grains were analyzed at the cores, and the rims of 1 to 4 large grains were analyzed to monitor intragranular trace element homogeneity. No significant differences between cores and rims have been observed. For Amp, measurements were made at the cores, because Amp grains are too small to allow separate respectively). This similarity between REE content in Cpx and Amp attests for the equilibrium between these two phases (e.g., [START_REF] Chazot | Determination of partition coefficients between apatite, clinopyroxene, amphibole, and melt in natural spinel lherzolites from Yemen: implications for wet melting of the lithospheric mantle[END_REF].

According to the extended trace element spider diagrams, Amp in these two samples have a slightly different patterns (Fig. 5b). Amp in sample MC36 is enriched in Nb relative to neighboring elements (U and La) whereas Amp in sample MC53 is Nb depleted. A similar difference also occurs for Ba. Compared with the distribution patterns of coexisting Cpx, those of Amp have significant distinctions: Sr, Ti enrichment, and Zr, Hf depletion.

Li concentrations and isotopic compositions

Li concentrations and isotopic compositions of minerals in peridotite xenoliths from Mont Coupet are reported in Table 3.

Li concentrations in Ol range from 1.1 to 2.8 ppm, with small variation of averages in individual samples (2.1 to 2.5 ppm; Table 3). In samples MC36 and MC53, Li contents in Amp cover a narrow range of 0.7 to 1.4 ppm. In spite of the very small grain size of Amp, profiles have also been measured on each grain and no core-rim heterogeneous Li distribution was observed.

Amp in both samples have δ 7 Li values lower than the coexisting mineral phases, with δ 7 Li values varying from -11.4 to -2.1‰. The intragranular variations are regarded within the analytical errors.

Water contents in minerals

Infrared spectra

The averaged Ol spectra collected under unpolarized light have bands similar to those observed in previous studies (Fig. 6a; e.g., [START_REF] Miller | The natural occurrence of hydroxide in olivine[END_REF][START_REF] Bell | Hydroxide in olivine: a quantitative determination of the absolute amount and calibration of the IR spectrum[END_REF][START_REF] Matsyuk | Hydroxyl in olivines from mantle xenoliths in kimberlites of the Siberian platform[END_REF][START_REF] Berry | Fingerprinting the water site in mantle olivine[END_REF][START_REF] Peslier | Hydrogen loss from olivines in mantle xenoliths from Simcoe (USA) and Mexico: mafic alkalic magma ascent rates and water budget of the sub-continental lithosphere[END_REF][START_REF] Grant | Water partitioning between mantle minerals from peridotite xenoliths[END_REF][START_REF] Walker | Three water sites in upper mantle olivine and the role of titanium in the water weakening mechanism[END_REF][START_REF] Peslier | Fast kimberlite ascent rates estimated from hydrogen diffusion profiles in xenolithic mantle olivines from southern Africa[END_REF][START_REF] Peslier | Water disequilibrium in olivines from Hawaiian peridotites: Recent metasomatism, H diffusion and magma ascent rates[END_REF][START_REF] Denis | Evidence of dehydration in peridotites from Eifel Volcanic Field and estimates of the rate of magma ascent[END_REF][START_REF] Denis | Water content and hydrogen behaviour during metasomatism in the uppermost mantle beneath Ray Pic volcano (Massif Central, France)[END_REF]Schämdicke et al., 2013;[START_REF] Doucet | High water contents in the Siberian cratonic mantle linked to metasomatism: An FTIR study of Udachnaya peridotite xenoliths[END_REF][START_REF] Demouchy | Characterization of hydration in the mantle lithosphere: Peridotite xenoliths from the Ontong Java Plateau as an example[END_REF]. Only bands at wavenumber 3572 and 3520 cm -1 are present at higher frequencies (Fig. 6a). Bands at lower frequencies are present in Ol spectra of all but one sample (MC38). A broad band centered at 3225 cm -1 occurs and in some samples (MC34, MC36 and MC53) this band is even stronger than the band at 3572 cm -1 (Fig. 6a). Two other bands centered around 3350 and 3325 cm -1 are present in the Ol spectra from sample MC34, but not in spectra of Ol from other samples. Core-rim profiles have been made on some large Ol grains, but no significant changes have been observed in the spectra collected from grain cores to rims, including band position, band width, and peak height (Fig. 7a).

The spectra collected under polarized light with E parallel to X, Y and Z have similar behavior to those reported in [START_REF] Denis | Evidence of dehydration in peridotites from Eifel Volcanic Field and estimates of the rate of magma ascent[END_REF]Fig. 6b). The bands at high frequency are dominant for E parallel to x axis, while the band centered at 3225 cm -1 is dominant in the spectra collected along z axis.

The averaged Opx spectra are shown in Fig. 6c. The band positions are similar to those reported in previous studies (e.g., [START_REF] Skogby | Hydroxide in pyroxene: variations in the natural environment[END_REF][START_REF] Grant | Water partitioning between mantle minerals from peridotite xenoliths[END_REF][START_REF] Xia | Low water content of the Cenozoic lithospheric mantle beneath the eastern part of the North China Craton[END_REF][START_REF] Sundvall | Water in upper mantle pyroxene megacrysts and xenocrysts: A survey study[END_REF][START_REF] Baptiste | Deformation and hydration of the lithospheric mantle beneath the Kaapvaal craton, South Africa[END_REF]Mosenfelder and Rossman, 2013a;[START_REF] Demouchy | Characterization of hydration in the mantle lithosphere: Peridotite xenoliths from the Ontong Java Plateau as an example[END_REF][START_REF] Denis | Evidence of dehydration in peridotites from Eifel Volcanic Field and estimates of the rate of magma ascent[END_REF][START_REF] Denis | Water content and hydrogen behaviour during metasomatism in the uppermost mantle beneath Ray Pic volcano (Massif Central, France)[END_REF][START_REF] Hao | Regional heterogeneity in the water content of the Cenozoic lithospheric mantle of Eastern China[END_REF]. The bands are divided into several groups, between 3605-3590 cm -1 , 3525-3515 cm -1 , 3420-3390 cm -1 , 3315-3300 cm -1 , and 3080-3060 cm -1 . In this study, the averaged spectra of Opx from sample MC34, MC43 and MC49 are very similar, with the presence of all the aforementioned bands and the similar relative strengths of different bands (Fig. 6c). In Opx spectrum of sample MC38, an additional band at wavenumbers between 3578 to 3560 cm -1 appears and it overwhelm the band 3605-3590 cm -1 . In contrast, the bands centered at 3315-3300 cm -1 and 3080-3060 cm -1 are absent in the spectrum. The spectra of Opx from two hydrous samples MC36 and MC53 resemble with each other; bands centered at 3315-3300 cm -1 and 3080-3060 cm -1 also can't be observed; and the band centered at 3525-3515 cm -1 contributes less to the total OH absorbance when compared with the same band in sample MC34, MC43 and MC49. Core-rim profiles were also measured on some large Opx grains, and the infrared spectra do not display any changes from grain cores to rims (Fig. 7b).

The structural OH absorption bands in the averaged Cpx spectra of Mont Coupet are located in three narrow wavenumber ranges, between 3650-3630 cm -1 , 3540-3520 cm -1 , and 3465-3450 cm -1 (Fig. 6d), comparable to literature (e.g., [START_REF] Skogby | Hydroxide in pyroxene: variations in the natural environment[END_REF][START_REF] Grant | Water partitioning between mantle minerals from peridotite xenoliths[END_REF][START_REF] Sundvall | Water in upper mantle pyroxene megacrysts and xenocrysts: A survey study[END_REF][START_REF] Baptiste | Deformation and hydration of the lithospheric mantle beneath the Kaapvaal craton, South Africa[END_REF]Mosenfelder and Rossman, 2013b). The three band types are present together in the averaged spectra of Cpx from all the Mont Coupet samples.

Water contents

The results of water contents in Ol calculated from IR spectra are shown in Table 4.

No bands induced by structural OH vibration can be observed in the averaged spectrum of Ol in sample MC38 (Fig. 6a). The averaging is carried out on spectra collected from more than twenty grains, so it is acceptable to assume that water content of Ol in these samples is 0 ppm (at <<1 ppm). The deviation of total water contents calculated using unpolarized spectra relative to the values estimated from polarized spectra is within the analytical errors (<30%; Table 6), so in the following sections the values from unpolarized spectra will be adopted to represent Ol water contents. Ol in Mont Coupet samples (except sample MC38) contain 0.6 to 2.1 ppm water. The maximum value occurs in sample MC34.

Opx and Cpx water contents estimated from the averaged spectra are shown in Water contents of Ol, Opx and Cpx in Amp-bearing samples (MC36 and MC53)

are not higher than in other samples (Table 5). This suggests that although samples MC36 and MC53 have once experienced hydrous metasomatism, water contents of primary phases may not have been affected at that time or the imprint of water addition has been erased by subsequent metasomatic events.

Discussion

Mantle metasomatism after depletion by extraction of melts

During melting, incompatible elements (e.g., Ca, Al) preferentially enter the melt, whereas other elements (e.g., Mg) selectively reside in the solid phase. Cr has a higher compatibility in solid phases than Al, and in dry and low-pressure environment, more Mg remains in residuals than Fe. With the incremental degrees of partial melting of peridotites, Fo content in Ol and Cr# of Spl increase whereas Al 2 O 3 contents in Cpx decrease (e.g., [START_REF] Arai | Characterization of spinel peridotites by olivine-spinel compositional relationships: review and interpretation[END_REF][START_REF] Hellebrand | Coupled major and trace elements as indicators of the extent of melting in mid-ocean-ridge peridotites[END_REF]. Thus, the positive co-variation between Fo contents and Cr# of Spl and the negative co-variation between Al 2 O 3 contents in Cpx and Cr# of Spl suggest that peridotites from Mont

Coupet have experienced partial melting in varying degrees (Fig. 8a and8b).

Mildly incompatible elements in Cpx (HREE + Y) are commonly unaffected by later metasomatic events and negatively co-relate to Cr# of co-existing Spl, reflecting the partial melting event (Fig. 8c). In contrast, LREE contents (Fig. 8d) seem to be affected by later metasomatic events, and do not display such negative correlations. [START_REF] Werling | Thermal evolution of the lithosphere beneath the French Massif Central as deduced from geothermobarometry on mantle xenoliths[END_REF][START_REF] Xu | Melt percolation and reaction atop a plume: evidence from the poikiloblastic peridotite xenoliths from Borée (Massif Central, France)[END_REF][START_REF] Touron | Characterization of the metasomatic agent in mantle xenoliths from Deves, Massif Central (France) using coupled in situ trace-element and O, Sr and Nd isotopic compositions[END_REF][START_REF] Yoshikawa | Geochemical and Sr-Nd isotopic characteristics and pressure-temperature estimates of mantle xenoliths from the French Massif Central: evidence for melting and multiple metasomatism by silicate-rich carbonatite and asthenospheric melts[END_REF] and falls outside the available range of the equation.

Cpx trace element distributions in Mont Coupet samples display transitional variation of REE from LREE depletion to LREE high enrichment, and HFSE negative anomalies in varying degrees. Up to present, no exclusive opinions have been generally accepted on the type of metasomatic agents to induce such characteristic of trace elements in primary peridotites. HFSE negative anomalies associated with LREE enrichment have been considered as a signature of either carbonatite-bearing melts/fluids (e.g., [START_REF] Rudnick | Carbonatite metasomatism in the northern Tanzanian mantle: petrographic and geochemical characteristics[END_REF][START_REF] Yaxley | Carbonatite metasomatism in the southeastern Australian lithosphere[END_REF][START_REF] Coltorti | Carbonatite metasomatism of the oceanic upper mantle: evidence from clinopyroxenes and glasses in ultramafic xenoliths of Grande Comore, Indian Ocean[END_REF][START_REF] Ionov | Mechanisms and sources of mantle metasomatism: major and trace element compositions of peridotite xenoliths from Spitsbergen in the context of numerical modelling[END_REF], or melts/fluids released during subduction (e.g., [START_REF] Zanetti | The Finero phlogopite-peridotite massif: an example of subduction-related metasomatism[END_REF][START_REF] Downes | Ultramafic xenoliths from the Bearpaw Mountains, Montana, USA: Evidence for multiple metasomatic events in the lithospheric mantle beneath the Wyoming craton[END_REF]Rivalent et al., 2007).

When melts/fluids percolate into a peridotite suite in depth, the crystallization of hydrous accessory minerals commonly occurs at the early stage [START_REF] Ionov | Mechanisms and sources of mantle metasomatism: major and trace element compositions of peridotite xenoliths from Spitsbergen in the context of numerical modelling[END_REF].

Thus, if the melt responsible for generation of Amp in peridotites is the same as the one induced LREE enrichment in primary Cpx, Amp should keep the trace element distribution pattern characterized by LREE enrichment. But in samples MC36 and MC53, Amp are depleted in LREE (Fig. 5a) and REE partitioning between Cpx and Amp have achieved equilibrium (Table 2). Therefore, it can be considered that prior to the percolation of melts/fluids leading to LREE enrichment and HFSE negative anomalies in Cpx, another metasomatic event should have occurred to account for the generation of Amp. Meanwhile, Zr-Hf negative anomalies and Sr-Ti positive anomalies are present in Amp of both samples, totally different from coexisting Cpx, where weak Ti negative anomalies and no Sr-Zr-Hf anomalies occur. But K d Amp/Cpx of each element in both samples is close to its equilibrium values (Table 2). The different distribution patterns of Zr, Hf and Ti in coexisting Amp and Cpx from samples MC36 and MC53 can therefore be ascribed to equilibrated redistribution of each element between the two phases at mantle depths [START_REF] Vannucci | Origin of LREE-depleted amphiboles in the subcontinental mantle[END_REF]. [START_REF] Coltorti | Amphiboles from suprasubduction and intraplate lithospheric mantle[END_REF] have compiled major and trace element compositions of Amp in mantle xenoliths from many localities of various tectonic backgrounds; they suggested that Zr/Nb ratio plotted against Ti/Nb ratios and Ti/Zr in Amp could trace the origin of Amp (Fig. 9a andb; [START_REF] Coltorti | Amphiboles from suprasubduction and intraplate lithospheric mantle[END_REF]. Considering Amp in both samples plus a larger dataset from [START_REF] Touron | Geochemical Fingerprints of Mantle Metasomatism Beneath the Massif Central[END_REF], we propose that the majority of Amp in peridotites from the FMC may crystallize from suprasubduction fluids or from two end-member mixed fluids (supersubduction and intraplate; Fig. 9a andb; [START_REF] Coltorti | Amphiboles from suprasubduction and intraplate lithospheric mantle[END_REF].

Large inter-mineral and inter-sample Li isotopic variation

Since the volcano erupted in the strombolian mode, peridotite xenoliths from Mont

Coupet should have been transported to the earth surface quite fast and cooled to the ambient temperatures in a very short period. Thus, the original information of Li isotopic systems at mantle depths should be preserved in xenoliths.

Li concentrations in minerals in Mont Coupet samples range from 0.4 to 2.8 ppm, in the range of normal mantle (Fig. 10), and the Li partitioning among different coexisting phases are very close to equilibrium (Fig. 10; Table 3). Meanwhile, in the measured profiles, Li concentrations display almost no variation across most of analyzed grains of Ol, Opx, Cpx and Amp (Supplementary figure 1). In some grains (mainly Cpx), small decreases towards rims (<1 ppm) are limited within narrow zones close to grain rims; in fact, they are just a little larger than analytical errors.

In contrast with Li inter-mineral equilibrated partitioning, Li isotopic compositions show larger variation among different phases (Table 3).The averaged δ 7 Li values of different mineral phases in most of Mont Coupet samples follow the common sequence: Ol>Opx>Cpx, except for sample MC53, where Li isotopic composition of Cpx is higher than those of coexisting Ol and Opx (Table 3). The averaged δ 7 Li values of Ol range from +0.5 to +7.7 ‰, close to the estimated isotopic compositions of the upper mantle (+2 to +5 ‰); some Opx and Cpx in some samples have negative averaged δ 7 Li values (Table 3).

Given the large difference of diffusive rates between 6 Li and 7 Li, with a diffusion coefficient ratio D 7 /D 6 =(m 6 /m 7 ) β =1.0424 (m 6 and m 7 are the masses of 6 Li and 7 Li, respectively; exponent β=0.27 according to [START_REF] Richter | Lithium isotope fractionation by diffusion in minerals. Part 1: Pyroxenes[END_REF]), the Li isotopic composition will be altered during Li exchange by diffusion (e.g., [START_REF] Lundstrom | Observations of Li isotopic variations in the Trinity Ophiolite: evidence for isotopic fractionation by diffusion during mantle melting[END_REF][START_REF] Wagner | Behaviour of Li and its isotopes during metasomatism of French Massif Central lherzolites[END_REF][START_REF] Gu | Multi-stage metasomatism revealed by trace element and Li isotope distributions in minerals of peridotite xenoliths from Allègre volcano (French Massif Central)[END_REF]. During diffusion, the minerals as the receptors first become enriched in 6 Li owing to its faster diffusion rate, and the isotopic compositions of minerals evolve towards negative values, producing a 'trough' in isotopic compositions [START_REF] Parkinson | High-temperature lithium isotope fractionation: insights from lithium isotope diffusion in magmatic systems[END_REF] The negative δ 7 Li values, measured recently in peridotite xenoliths and mantle-derived magmas by more and more studies (e.g., [START_REF] Li | Anomalous lithium isotopic compositions of the Cenozoic lithospheric mantle beneath Penglai, Shandong Province: the ion probe analyses of peridotite xenoliths[END_REF][START_REF] Tang | Slab-derived lithium isotopic signatures in mantle xenoliths from northeastern North China Craton[END_REF][START_REF] Tang | Abnormal lithium isotope composition from the ancient lithospheric mantle beneath the North China Craton[END_REF][START_REF] Ackerman | Alkaline and carbonate-rich melt metasomatism and melting of subcontinental lithospheric mantle: evidence from mantle xenoliths, NE Bavaria, Bohemian Massif[END_REF]Adelfadil et al., 2014;[START_REF] Tian | The anomalous lithium isotopic signature of Himalayan collisional zone carbonatites in western Sichuan, SW China: Enriched mantle source and petrogenesis[END_REF][START_REF] Gu | Multi-stage metasomatism revealed by trace element and Li isotope distributions in minerals of peridotite xenoliths from Allègre volcano (French Massif Central)[END_REF], have been originally related to the recycled materials transported by subducted slabs into depths, forming a low-δ 7 Li reservoir in the deep mantle [START_REF] Zack | Extremely light Li in orogenic eclogites: the role of isotope fractionation during dehydration in subducted oceanic crust[END_REF].

Here, melts/fluids stemming from the subducted materials during Variscan orogeny may have metasomatized the lithospheric mantle beneath the FMC, and brought the light Li isotopic compositions [START_REF] Gu | Multi-stage metasomatism revealed by trace element and Li isotope distributions in minerals of peridotite xenoliths from Allègre volcano (French Massif Central)[END_REF].

Water content variation controlled by partial melting and metasomatism

Preservation of original water contents in minerals

Due to the high diffusivity of hydrogen at magmatic temperatures (e.g., Mackwell and Kohstedt, 1990;[START_REF] Hercule | Hydrogen in diopside: Diffusion, kinetics of extraction-incorporation, and solubility[END_REF][START_REF] Stalder | Hydrogen diffusion in natural and synthetic orthopyroxene[END_REF][START_REF] Demouchy | Mechanisms of hydrogen incorporation and diffusion in iron-bearing olivine[END_REF] and water solubility in NAMs decreasing with diminishing pressures (e.g., [START_REF] Kohlstedt | Solubility of water in the α, β and γ phases of (Mg, Fe) 2 SiO 4[END_REF]Fé 19 and Bolfan-Casanova, 2012), hydrogen incorporated in lattices of NAMs may be lost during xenolith ascent to the surface (Demouchy et al., 2006;[START_REF] Peslier | Hydrogen loss from olivines in mantle xenoliths from Simcoe (USA) and Mexico: mafic alkalic magma ascent rates and water budget of the sub-continental lithosphere[END_REF]. Although experimentally determined H diffusivity does not vary largely in Ol and pyroxenes (e.g., [START_REF] Hercule | Hydrogen in diopside: Diffusion, kinetics of extraction-incorporation, and solubility[END_REF][START_REF] Stalder | Hydrogen diffusion in natural and synthetic orthopyroxene[END_REF][START_REF] Demouchy | Mechanisms of hydrogen incorporation and diffusion in iron-bearing olivine[END_REF][START_REF] Ingrin | Diffusion of hydrogen in minerals[END_REF], profiles formed by H diffusive loss have been observed in xenolithic Ol but not in coexisting pyroxenes (e.g., [START_REF] Peslier | Hydrogen loss from olivines in mantle xenoliths from Simcoe (USA) and Mexico: mafic alkalic magma ascent rates and water budget of the sub-continental lithosphere[END_REF] except in Opx in one recent specific case [START_REF] Tian | Water concentration profiles in natural mantle orthopyroxenes: a geochronometer for long annealing of xenoliths within magma[END_REF]. In several recent studies, a certain amount of structural OH have been detected in Ol grains in peridotite xenoliths transported by alkali basalts and have been considered to represent the original amount prior to entrapment by magmas [START_REF] Demouchy | Characterization of hydration in the mantle lithosphere: Peridotite xenoliths from the Ontong Java Plateau as an example[END_REF][START_REF] Denis | Water content and hydrogen behaviour during metasomatism in the uppermost mantle beneath Ray Pic volcano (Massif Central, France)[END_REF]. There is no reason to believe that the host magmas of xenoliths from Mont Coupet have risen much slower than these alkali basalts and that main minerals have lost water through decompression-induced outward diffusion.

The homogeneous intragranular water distribution in Opx and Cpx grains from Mont Coupet xenoliths has been testified by the absence of variation in spectra collected from grain cores to rims (e.g., Fig. 7b). And the concentration ratios between Opx and Cpx in most samples from Mont Coupet, except two Amp-bearing samples MC36 and MC53, fall in the range restricted by mantle xenoliths from worldwide localities (Fig. 12a). They indicate that Opx and Cpx from Mont Coupet xenoliths should have preserved their initial water contents.

Ol in xenoliths from Mont Coupet does not show evidence of water loss during ascent to the surface. On the one hand, the FTIR profile analyses reveal no variation in spectra across Ol grains. Especially, for the profile in one grain from sample MC34, both the bandwidth and the peak height of the band at 3225 cm -1 do notshow any obvious changes (Fig. 7a). This band has been assigned to the mechanism of two H substituting Mg 2+ , by which hydrogen incorporated into Ol lattices [START_REF] Balan | Theoretical infrared spectrum of OH-defects in forsterite[END_REF]. It has been experimentally shown that this type of defect is one of fastest among the four mechanisms of H incorporation into Ol [START_REF] Padrón-Navarta | Site-specific hydrogen diffusion rates in forsterite[END_REF].

If hydrogen loss would have occurred in Ol, it can be expected that the band would display variation in spectra collected from core to rim, or have totally disappeared, which is in contradiction with the observations in Fig. 7a. On the other hand, water contents of Ol in Mont Coupet xenoliths have a roughly positive co-variation with those of co-existing Cpx (Fig. 12b).

Ol in xenoliths from Mont Coupet contain rather less hydrogen compared with Ol in cratonic peridotite xenoliths (e.g., [START_REF] Grant | Water partitioning between mantle minerals from peridotite xenoliths[END_REF][START_REF] Peslier | Fast kimberlite ascent rates estimated from hydrogen diffusion profiles in xenolithic mantle olivines from southern Africa[END_REF][START_REF] Doucet | High water contents in the Siberian cratonic mantle linked to metasomatism: An FTIR study of Udachnaya peridotite xenoliths[END_REF], but it is consistent with observations from other off-craton mantle xenoliths (e.g., [START_REF] Grant | Water partitioning between mantle minerals from peridotite xenoliths[END_REF][START_REF] Peslier | A review of water contents of nominally anhydrous natural minerals in the mantles of Earth, Mars and the Moon[END_REF] and especially compared to off-craton xenoliths also from the FMC [START_REF] Denis | Water content and hydrogen behaviour during metasomatism in the uppermost mantle beneath Ray Pic volcano (Massif Central, France)[END_REF]. It is in agreement with the general observation that Ol in the off-craton lithospheric mantle might indeed have lower water content than Ol from cratonic lithospheric mantle.

Water content variation controlled by different processes

Among xenoliths from Mont Coupet, Opx and Cpx in sample MC38 contain the least water compared with corresponding minerals in other samples, and Ol infrared spectra from sample MC38 do not display any OH-vibration bands at wavenumbers between 3700-3000 cm -1 (Fig. 6a), indicating that water content in Ol of this sample is lower than the detect limit of FTIR. According to major element compositions and Cpx HREE concentrations, sample MC38 had been depleted strongly by partial melting, with calculated partial melting degree of 16%, the highest among samples from Mont Coupet. Considering H as an incompatible element, the high-degree of partial melting of sample MC38 is coherent with a lower water content compared with other samples (Fig. 13a andb). Although LREE enrichment and HFSE negative anomalies in Cpx of sample MC38 points to a metasomatism subsequent to partial melting, the result of water content suggests that this metasomatism has at least not added water to the sample. This metasomatism is probably responsible for the alteration of the nature of main H defects in Opx and Cpx since their spectra are markedly different from those in other samples from Mont Coupet (Fig. 6c andd).

The band position shift is intrinsically related to the change of OH-stretching modes, corresponding to the change of P-T conditions or chemical compositions of host minerals (e.g., [START_REF] Yang | OH in natural orthopyroxene: an in situ FTIR investigation at varying temperatures[END_REF][START_REF] Ingrin | Low-temperature evolution of OH bands in synthetic forsterite, implication for the nature of H defects at high pressure[END_REF], which may have been induced

by the percolation of the metasomatic melt. Furthermore, sample MC38 also has the highest calculated Spl Fe 3+ /ΣFe value (Fig. 13c), consistent with its low water content.

This consistency has been observed in natural samples (e.g., [START_REF] Peslier | Low water contents in pyroxenes from spinel-peridotites of the oxidized, sub-arc mantle wedge[END_REF] and an experiment by [START_REF] Bromiley | Hydrogen solubility and speciation in natural, gem-quality chromian diopside[END_REF]. Everything is pointing toward a very low water activity of the fluid/melt involved in this particular metasomatism event.

Sample MC34 severely deviates from the roughly negative co-variation between Cpx/Opx water contents and Fo contents in Ol for samples from Mont Coupet (Fig. 13a andb). The metasomatism, characterized by LREE enrichment in Cpx of sample MC34, might have added some water to sample MC34. As mentioned above, sample MC38 display strongly enriched characteristics of trace element compositions in its

Cpx, but its water contents in minerals are the lowest among all the samples from Mont Coupet. The different effects of metasomatism on water contents point out that different percolating agents may have carried out the metasomatism in samples MC34 and MC38.

The whole rock water contents in samples MC43 and MC49, with values of 50 and 88 ppm, fall in the range of water content in depleted MORB source mantle (50-200 ppm). Their water contents in minerals may just be involved in partial melting but not major postdating metasomatism (Fig. 13a, b andd).

For the two Amp-bearing samples MC36 and MC53, water contents in Ol, Opx and Cpx are not remarkably higher than those in equivalent minerals from other anhydrous samples from Mont Coupet (Table 5; Fig. 12a andb). In fact, the highest water contents of Cpx and Opx are observed in anhydrous sample MC34 (Table 5).

The same observations have been reported in previous studies [START_REF] Yang | Water in minerals of the continental lithospheric mantle and overlying lower crust: a comparative study of peridotite and granulite xenoliths from the North China Craton[END_REF][START_REF] Schmädicke | Olivine from spinel peridotite xenoliths: Hydroxyl incorporation and mineral composition[END_REF][START_REF] Bonadiman | Crystal chemistry of amphiboles: implications for oxygen fugacity and water activity in lithospheric mantle beneath victoria land, antarctica[END_REF][START_REF] Denis | Water content and hydrogen behaviour during metasomatism in the uppermost mantle beneath Ray Pic volcano (Massif Central, France)[END_REF]. Bonadiman et al.

(2014) have calculated a very low water activity in Amp-bearing mantle peridotites, accompanied by low oxygen fugacity. They proposed that Amp could be regarded as a repository of water and have prevented circulation of water in the local environment composed of NAMs. From another viewpoint, the generation of Amp consumes water in the system and largely decreases the water activity in mantle peridotites. Thus, water contents in primary NAMs is not remarkably increased. In addition, as discussed in section 6.1, the formation of Amp may have occurred in another metasomatic event prior to the metasomatism characterized by LREE enrichment and HFSE negative anomalies in Cpx. The imprints from the foregoing modal metasomatism might have been modified by the lately percolating melts/fluids.

However, it is noteworthy that these two Amp-bearing samples have higher / ratios than the / ratio range restricted by mantle xenoliths from worldwide localities (Table 5; Fig. 12a).

Relationship between water contents and Ti concentrations in Ol

Although H can incorporate into Ol lattice through occupying Mg vacancies charge-balanced by trivalent ions (e.g., Al 3+ Cr 3+ , Fe 3+ ), its corresponding bands are only observed in Ol infrared spectra of sample MC34 but not in those of other

C H 2 O Cpx C H 2 O Opx C H 2 O Cpx C H 2 O Opx
samples. Thus, we will not discuss the relationship between concentrations of trivalent ions and concentrations of H responsible for bands at between 3400 and 3300 cm -1 .

The maximum concentration of H, entering Ol accompanied by Ti, is calculated from Ti concentrations in Ol according to the expression Ti 4+ + 2H + =Mg 2+ + Si 4+ , supposing that all Ti contributes to the Ti-2H defect. marked by a blue circle; a photo of the Razas Grand quarry from where the samples come is also shown on the lower right corner. Another locality (Allègre), from which peridotite xenoliths have also been studied [START_REF] Gu | Multi-stage metasomatism revealed by trace element and Li isotope distributions in minerals of peridotite xenoliths from Allègre volcano (French Massif Central)[END_REF], is also marked by a red circle for comparison. The boundary between the northern and southern domains (the horizontal black line; [START_REF] Lenoir | Contrasting lithospheric mantle domains beneath the Massif Central (France) revealed by geochemistry of peridotite xenoliths[END_REF] is also drawn. (d). [START_REF] Mcdonough | The composition of the Earth[END_REF].

C H 2 O Cpx C H 2 O Opx
Ti*=2×Ti n /(Eu n +Gd n ); Zr * =2×Zr n /(Nd n +Sm n ) where the subscript 'n' represents normalization to the primitive mantle. K d in the last row lists the trace element partition coefficients between Amp and Cpx calculated from the partition coefficients between Cpx (and Amp) and melts [START_REF] Ionov | Mechanisms and sources of mantle metasomatism: major and trace element compositions of peridotite xenoliths from Spitsbergen in the context of numerical modelling[END_REF]. spectra. The subscripts "B" and "L" denote the calculation based on the methods from [START_REF] Bell | Hydroxide in olivine: a quantitative determination of the absolute amount and calibration of the IR spectrum[END_REF] 3

Cpx

and [START_REF] Libowitzky | An IR absorption calibration for water in minerals[END_REF], respectively. And the subscripts "un" and "p" represent the 4 results calculated from the spectra under unpolarized and polarized light, respectively.

5

Deviation=(OH un -OH p )/OH p ×100, OH un (OH p ) denotes the total water contents (Group I + Group II) 6 calculated based on unpolarized (polarized) infrared spectra; A T = A X +A Y +A Z .

  measurements at rims. No significant variation was observed among different Cpx or among different Amp grains in any given sample; thus, we will use the average contents of individual samples in the following discussion. The chondrite-normalized REE distribution of Cpx displays different patterns among Mont Coupet samples (Fig.4a). Cpx in sample MC43 is slightly depleted in LREE relative to middle MREE, and has a nearly flat transition from MREE to HREE ((La/Yb) N =0.41; (Sm/Yb) N =0.96; the subscript 'N' represents normalization to chrondrite; McDonough and Sun, 1995). Cpx chondrite-normalized REE distribution patterns in Amp-bearing sample MC36 and MC53 ((La/Yb) N =0.25 and 0.32, respectively; (Sm/Yb) N =0.80 and 0.91, respectively), resemble with that in sample MC43. In sample MC49, the normalized pattern shows obvious enrichment in La compared with other elements and nearly no fluctuation from Ce N to Lu N with (La/Yb) N and (Sm/Yb) N values of 1.39 and 0.87, respectively. The distribution pattern of sample MC34 exhibit upward inflection from Nd to La, and flat transition from MREE to HREE ((La/Sm) N =3.20; (Sm/Yb) N =1.44). Sample MC38 is totally different. The sample is extremely depleted in MREE and HREE with Gd N and Er N values of 1.03 and 0.67, respectively, whereas it is highly enriched in LREE with La N value of 35.47. The pattern shows linear increase from Er to La. The extended trace element compositions of Cpx in xenoliths from Mont Coupet are characterized by HFSE negative anomalies when element abundances are normalized to the primitive mantle (Fig. 4b). Ti and Nb anomaly occurs in all the Mont Coupet samples. The magnitude of Ti negative anomaly gradually decreases from sample MC38 (Ti*=0.10) to sample MC53 (Ti*=0.63; Ti*=2×Ti n /(Eu n +Gd n ); the subscript 'n' represents normalization to the primitive mantle). A very strong Zr negative anomaly is shown by sample MC38. Sample MC34 displays a slight Zr negative anomaly and in other samples no Zr anomaly is observed. U and Thin Cpx of all the samples display marked enrichment. Amp REE have similar chondrite-normalized distribution patterns in sample MC36 and MC53 (Fig. 5a). They are very close to the patterns of coexisting Cpx, exhibiting LREE depletion and nearly flat transition from MREE to HREE in spider diagrams ((La/Yb) N =0.26 and 0.30; (Sm/Yb) N =0.61 and 0.90 in sample MC36 and MC53,

  Fig.11a). And, there are no negative co-variations between Cpx δ 7 Li values and the Li concentrations when taking all the Mont Coupet samples into consideration (Fig.11a). Since Li partitioning is at equilibrium among minerals and no signatures of Li diffusive addition into minerals have been observed (Fig.10and Supplementary figure1), the signatures of Li isotopic compositions imply that the light isotopic compositions of Opx and Cpx in some samples have been modified by an exotic small-volume melt with a light isotopic composition, and low or moderate Li content. In light of small or no core-rim concentration variations, Opx and Cpx grains in Mont Coupet samples may have achieved Li isotopic equilibrium with the percolating melt in a short time. Due to the lower diffusive rate of Li in Ol than Cpx, and higher Li content of Ol, Li isotopic compositions of Ol cannot have been largely affected.In fact, Li isotopic compositions in Cpx from Mont Coupet xenoliths display a gradual increase from sample MC38 to MC53 (Fig.11b). Here, an interpretation similar to the chromatographic effect described for REE distribution could be invoked to explain the Cpx Li isotopic composition of the different samples. Once the melt has percolated into the lithospheric mantle beneath Mont Coupet, harzburgites represented by MC38 were first passed through and the Li isotopic compositions in their Cpx (also Opx) grains were assimilated in a short time whereas the Li isotopic compositions of Ol were influenced in a smaller degree. As the melt left and continued advancing towards other parts of the mantle, δ 7 Li values of Cpx in sample MC38 was no longer changed by the percolating melt. With the melt moving among different parts of the mantle, more 6 Li than 7 Li diffused into ambient peridotites along the melt passage, and its light Li isotopic composition was gradually diluted. Thus, the samples far from the melt path could not record negative δ 7 Li signatures. Considering the high diffusive rate of Li at high temperatures, in a single sample the inter-mineral Li isotopic disequilibrium is a transient state. The preservation of inter-mineral and inter-sample large Li isotopic fractionation currently observed in Mont Coupet implies that the percolation of the melt should have occurred shortly prior to the entrainment of the peridotite xenoliths by the host magmas, and/or through focused channels.

Fig. 1

 1 Fig. 1 Simplified geological map of the FMC and the spatial distribution of Cenozoic volcanoes. The sampling locality of Mont Coupet, belonging to Dèves (region 2), is

Fig. 2

 2 Fig. 2 (a) Δ i plotted against Mg# of standards. We used the averaged Δ i of individual standard minerals (pyroxenes in blue and Ol in red) to calculate δ 7 Li values of the equivalent minerals in studied samples. (b) Lithium ion yields ([Li]) plotted against Mg# of standards. No co-variation exists. For calculation, we use the averages of individual phases (i.e., Ol, Opx and Cpx, respectively).

Fig. 3

 3 Fig. 3 (a) Protogranular-textured sample MC50. Small-grain Cpx and Spl dispersed interstitially with larger coexisting Ol and Opx. (b) Porphyroclastic sample MC30 with fine Cpx exsolutions in a large Opx porphyroclast. (c) Kink-banded texture in an Ol grain from sample MC30. (d) Cribriform-shaped Opx grain boundaries in sample MC53.

Fig. 5

 5 Fig. 5 Amp REE distribution patterns normalized to Chondrite (a) and extended trace element distribution patterns normalized to the primitive mantle (b). The dashed lines in graph (b) represent two samples from another locality in the southern domain of the

Fig. 6

 6 Fig. 6 Averaged unpolarized infrared spectra of Ol (a), Opx (c) and Cpx (d); representative polarized Ol infrared spectra (b) collected in sample MC36 with E parallel to principle axes X, Y and Z. The dashed lines represent the peak positions with corresponding wavenumbers.

Fig. 9

 9 Fig.9Zr/Nb ratios vs. Ti/Nb ratios (a) and Ti/Zr (b) in Amp from samples MC36 and MC53. S-Amp and I-Amp denote Amp originally related to suprasubduction and intraplate environments, respectively. The fields of different Amp types are defined according to[START_REF] Coltorti | Amphiboles from suprasubduction and intraplate lithospheric mantle[END_REF].

Fig. 10

 10 Fig. 10 Li partitioning between Ol and Cpx (a), between Opx and Cpx (b) in peridotite xenoliths from Mont Coupet. The equilibrated lines are plotted according to the equilibrated partition coefficients from Seitz and Woodland (2000), D Ol/Cpx =1.5-2.0 and D Opx/Cpx =0.5-1.1. Li contents of individual mineral phases in the normal upper mantle is delimitated by the red square. The metasomatic tendency is defined according to[START_REF] Seitz | The distribution of lithium in peridotitic and pyroxenitic mantle lithologies-an indicator of magmatic and metasomatic processes[END_REF] and[START_REF] Tang | Abnormal lithium isotope composition from the ancient lithospheric mantle beneath the North China Craton[END_REF].

Fig. 11 δ 7

 7 Fig. 11 δ 7 Li of Cpx plotted against Li concentrations in Cpx (a), and Mg# of coexisting Ol (b). In graph (a), along the direction of the dashed line, δ 7 Li of Cpx increase with incremental Li concentrations, in contrary to the view that Li uptake by

Fig. 13

 13 Fig. 13 Water contents in Cpx (a) and in Opx (b) plotted against Mg# of Ol; water contents in Cpx versus Fe 3+ /ΣFe values of coexisting Spl (c), and Sm contents in Cpx

  Fig. 1 17 18 19

  Fig. 3 27 28 29 30

  Fig. 6 51 52 53 54

  Fig. 10 91

  Fig. 12 101

Table 5 ,

 5 where Ol water contents are also included for comparison. Opx water contents in Mont Coupet samples range from 28 to 152 ppm, the maximum value measured in sample MC34 and the minimum one in sample MC38. Cpx in sample

	MC38 have water content of 41 ppm, almost one order of magnitude lower than other
	samples (304-428 ppm). The water content ratio between Cpx and coexisting Opx (
	/	) in each sample varies from 1.5 to 4.5.

  . Inversely, δ 7 Li values of the donors (percolating melts or fluids), which have lost Li, increase gradually. In receptor minerals, the 'trough' propagates into the interiors of grains until a new equilibrium is achieved, or the system cools below the closure temperature, freezing the disequilibrium state. Accordingly, the disequilibrated Li diffusive uptake can decrease the Li isotopic compositions of host minerals in certain degrees. In xenoliths from Mont Coupet, the lowest Cpx δ 7 Li value was measured in sample MC38, which means that Cpx in this sample should contain much more Li than other samples if the light Li isotopic composition is interpreted as the result of dis-equilibrated kinetic diffusion. However, in contrary to this assertion, Li concentration of Cpx in sample MC38 is the lowest among all in Mont Coupet (

  Li value, and low to moderate Li concentration, originally related to recycled materials. This metasomatic process should occur shortly before the xenolith entrainment by the host magmas. Minerals in Mont Coupet xenoliths, including Ol, have apparently retained their original water contents. The distinct OH behaviors in individual samples with marked Cpx LREE enrichment suggest that agents of different origins have conducted the metasomatic processes characterized by variations of trace element compositions inCpx. Ol, Opx and Cpx in two Amp-bearing samples show no remarkable difference in water content than other xenoliths. This point may be interpreted as the generation of Amp, as a repository of water, has significantly lowered the water activity in the whole system, so that the modal metasomatism have not imprinted the water content in NAMs with high values, or the effect of this modal metasomatism has perhaps been erased by some late metasomatism (e.g., the metasomatism characterized by LREE enrichment in Cpx).

	The calculated results are
	expressed in ppm wt. H 2 O ranging from 1.39 ppm in sample MC38 to 5.25 ppm in
	sample MC43, much higher than the total Ol water contents measured in
	corresponding samples, let alone the content of water associated with bands at 3572
	and 3525 cm -1 (Table 6; Schämdicke et al., 2013). It therefore suggests that Ti content
	was not a limiting factor in the incorporation of water in the mantle beneath Mont
	Coupet.
	Conclusion
	Different geochemical indicators in peridotite xenoliths from Mont Coupet
	uniformly illustrate that the SCLM beneath the FMC has undergone multi-stage
	metasomatism by melts/fluids of different origins. The H and Li homogeneous
	distributions at the mineral scale in each sample point out that the host magmas have
	erupted fast enough to preserve the original information from Li and H. Li
	concentrations and isotopic compositions, and water contents in minerals have been
	modified in different degrees by metasomatism at mantle depths.
	LREE enrichment and negative HFSE negative anomalies in Cpx point out a
	metasomatic process by carbonatitic or subduction-related melts/fluids. Prior to this
	event, another metasomatic event accounts for the generation of Amp in two xenoliths,
	because the trace element distribution in Amp shows characteristics of LREE
	depletion and equilibrium with co-existing Cpx.
	Minerals in peridotite xenoliths from Mont Coupet have Li concentrations similar
	to the normal mantle; but Li isotopic compositions display significant inter-mineral
	fractionation and large variations among different samples, which should be a result
	of percolation of a small-volume melt with a low δ 7

Table . 1

 . Major element compositions of minerals in peridotite xenoliths from Mont Coupet. The values are the averages of all the analyzed points ('n' denotes the number) on every mineral phase in a single sample. Abbreviations: Fe T O, total Fe; Mg#=molar (100×Mg/(Mg+Fe T )); Cr#= molar (100×Cr/(Cr+Al)).

			0.05 59.59	0.08	0.00	0.01	0.00	0.00 20.01 10.95	8.00 98.70	76.5	8.3
	Amp	n	SiO 2 Al 2 O 3 TiO 2	CaO	Na 2 O K 2 O MnO MgO Fe T O Cr 2 O 3 Total	Mg#
	MC36	12 42.75 14.30	2.76 11.39	3.62 0.12	0.02 16.71	4.00	0.78 96.46	88.1
	Ol MC53	n 15 42.87 14.53 SiO 2 Al 2 O 3 TiO 2 3.46 11.33 CaO	Na 2 O K 2 O MnO MgO Fe T O Cr 2 O 3 3.48 0.20 0.03 16.61 4.00 0.78 97.28 Total	Mg# 88.1
	MC30	23 40.65	0.02	0.02	0.03	0.01 0.00	0.16 48.42	9.64	0.01	98.96	89.9
	MC34	18 41.37	0.01	0.01	0.03	0.01 0.00	0.12 49.84	8.46	0.02	99.88	91.3
	MC36	19 41.21	0.00	0.02	0.02	0.01 0.01	0.14 48.51 10.01	0.04	99.96	89.6
	MC38	32 41.58	0.01	0.02	0.03	0.01 0.01	0.12 49.64	8.57	0.04	100.01	91.1
	MC43	21 41.66	0.02	0.01	0.06	0.01 0.00	0.13 49.16	9.63	0.00	100.68	90.1
	MC49	33 41.98	0.01	0.00	0.05	0.01 0.00	0.13 48.66 10.01	0.00	100.84	89.6
	MC50	15 40.92	0.00	0.00	0.02	0.00 0.00	0.13 48.54	9.61	0.00	99.21	90.0
	MC53	18 40.79	0.01	0.01	0.03	0.01 0.00	0.14 48.39 10.00	0.01	99.39	89.6
	Opx	n	SiO 2 Al 2 O 3 TiO 2	CaO	Na 2 O K 2 O MnO MgO Fe T O Cr 2 O 3	Total	Mg#
	MC30	15 54.99	4.01	0.11	0.41	0.05 0.00	0.14 32.53	6.37	0.36	98.97 90.1
	MC34	23 56.18	2.87	0.06	0.43	0.04 0.00	0.10 34.25	5.67	0.46	100.05 91.5
	MC36	19 55.74	3.85	0.10	0.43	0.04 0.00	0.13 32.98	6.64	0.27	100.19 89.8
	MC38	14 57.24	2.00	0.01	0.50	0.04 0.01	0.14 34.31	5.78	0.41	100.44 91.3
	MC43	21 56.70	3.64	0.10	0.45	0.06 0.00	0.12 33.51	6.36	0.32	101.27 90.4
	MC49	12 57.11	3.44	0.12	0.43	0.05 0.00	0.15 33.23	6.56	0.27	101.36 90.0
	MC50	18 56.50	3.57	0.06	0.31	0.03 0.01	0.13 33.21	6.33	0.39	100.55 90.3
	MC53	24 54.98	4.00	0.10	0.40	0.05 0.01	0.14 32.92	6.61	0.27	99.48 89.9
	Cpx	n	SiO 2 Al 2 O 3 TiO 2	CaO	Na 2 O K 2 O MnO MgO Fe T O Cr 2 O 3	Total	Mg#
	MC30	32 51.64	6.66	0.62 20.67	1.91 0.00	0.05 13.93	2.34	0.77	98.61	91.4
	MC34	17 52.40	5.21	0.23 21.23	1.65 0.01	0.01 15.07	2.29	1.38	99.48	92.1
	MC36	22 52.34	6.40	0.55 21.31	1.64 0.00	0.05 14.29	2.33	0.69	99.62	91.6
	MC38	11 54.41	1.94	0.06 23.05	0.62 0.01	0.04 16.75	2.15	0.61	99.64	93.3
	MC43	23 52.83	6.44	0.56 20.76	1.60 0.00	0.03 14.93	2.82	0.86	100.84	90.4
	MC49	30 53.41	6.46	0.56 20.89	1.78 0.00	0.07 14.65	2.64	0.82	101.28	90.8
	MC50	18 52.74	6.02	0.37 21.41	1.81 0.00	0.05 14.44	2.30	1.00	100.13	91.8
	MC53	20 51.45	6.99	0.69 20.94	1.89 0.00	0.06 14.21	2.45	0.66	99.33	91.2
	Spl	n SiO 2 Al 2 O 3 TiO 2	CaO	Na 2 O	K 2 O	MnO MgO Fe T O Cr 2 O 3 Total	Mg# Cr#
	MC30	8	0.03 55.54	0.05	0.00	0.00	0.00	0.00 19.32 11.89	9.48 96.31	74.3 10.3
	MC34	6	0.02 44.33	0.05	0.01	0.01	0.01	0.00 18.23 13.97 22.85 99.48	69.9 25.7
	MC36	6	0.04 57.78	0.04	0.00	0.00	0.00	0.00 19.77 12.02	8.54 98.19	74.5	9.0
	MC38	6	0.03 28.30	0.03	0.00	0.00	0.00	0.00 14.48 19.59 34.80 97.23	56.8 45.2
	MC43	7	0.04 54.43	0.05	0.00	0.00	0.00	0.00 20.38 11.88 11.25 97.66	75.3 12.2
	MC49	7	0.03 55.38	0.05	0.01	0.00	0.00	0.00 19.84 12.75 10.47 98.15	73.5 11.3
	MC50	7 -0.01 54.26	0.03	0.00	0.00	0.01	0.00 19.41 11.90 12.44 97.62	74.4 13.3

Table 2 :

 2 Trace element abundances in ppm for Cpx of Mont Coupet peridotite xenoliths. The subscript 'N' represents normalization to the compositions of Chondrite (McDonough

Table 4 :

 4 Ol water contents of minerals in Mont Coupet peridotite xenoliths calculated from infrared 2

	Amp	K d	Cpx/Amp	K d

Table 5 :

 5 Water contents of Opx and Cpx in Mont Coupet peridotite xenoliths calculated from 9 unpolarized infrared spectra. "n" denotes the numbers of analyzed grains whose spectra have been 10 selected to calculate the averaged spectra . In that series, the numbers bracketed are the total numbers 11 of analyzed grains. Ol water content is also shown for comparison. O in GI represent the total Ol water content measured by FTIR and the measured water content just 14 associated with bands in Group I region, respectively.

		7
	Sample t (cm)	n	Ol Group I un (cm -1 ) (ppm) Integral OH B	Ol Group II un (cm -1 ) (ppm) Integral OH L	Ol Group I p (cm -1 ) (ppm) A T OH B	Ol Group II p (cm -1 ) (ppm) A T OH L	OH un (ppm)	OH p (ppm)	Deviation
	MC34	0.015	38	0.01	0.4	0.123	1 0.024	0.3	0.29	0.7	1.4	1.1	0.27
	MC36	0.016	32	0.04	1.4	0.09	0.7 0.107	1.3 0.215	0.6	2.1	1.9	0.11
	MC38	0.016	34	-	-	-	-	-	-	-	-	-	-	-
	MC43	0.017	17	0.006	0.2	0.061	0.4	-	-	-	-	0.6	-	-
	MC49	0.015	26	0.01	0.4	0.09	0.7 0.025	0.3 0.243	0.7	1.1	1	0.1
	MC53	0.016	29	0.04	1.4	0.075	0.6	-	-	-	-	2	-	-

Acknowledgements

This research was supported by the INSU SYSTER program and the State Scholarship Fund of China. We thank Michel FIALIN for access to EMP analysis, Marie-Christine BOIRON and Chantal PEIFFERT for help in analyzing trace elements using LA-ICPMS. GU X.Y. also appreciates Marc CHAUSSIDON, Andrey GURENKO for direction in using SIMS at CRPG. Special gratitude is expressed to Nordine BOUDEN for assistance in guaranteeing the operation of IMS 1270 when GU X.Y. used it. Isabelle De Waele is thanked for her patient direction and aide for the FTIR-microscope to analyze the water content in Université Lille 1.