Multistage origin of dunite in the Purang ophiolite, southern Tibet, documented by composition, exsolution and Li isotope characteristics of constituent minerals

The Purang ophiolite, which crops out over an area of about 650 km2 in the western Yarlung–Zangbo suture zone, consists chiefly of mantle peridotite, pyroxenite and gabbro. The mantle peridotite is comprised mainly harzburgite and minor dunite. Locally, the latter contains small pods of chromitite. Pyroxenite and gabbro occur as veins of variable size within the peridotite; most of them strike northwest, parallel to the main structure of the ophiolite. Three types of dunite occur in the Purang ophiolite: dunite that envelopes podiform chromitite (1) and lenses of dunite with either Cr-rich spinel (2) or Cr-poor spinel (3) in a harzburgite host. The constituent minerals of dunite envelopes around podiform chromitite are similar in composition to those of transition-zone dunite (Fo91.01−91.87 in olivine; Cr/(Cr+Al) (Cr#) = 41.5–47.0 and Mg/(Mg+Fe2+) (Mg#) = 58.9–63.0 in Cr-spinel). Forsterite contents in olivine decrease from type 2 lenses with Cr-rich spinel (91.9–93.0) to type 1 dunite enveloping chromitite (91.7–93.7) to type 3 lenses with Cr-poor spinel (95.3–96.0). Similarly, Cr# in spinel decreases from type 2 (66.9–67.9) to type 1 (41.5–47.0) to type 3 (19.8–20.6). In addition, Al2O3 in clinopyroxene is highest in type 2 (3.48–5.24 wt %) and decreases to type 1 (1.56–3.29 wt %) and type 3 (0.78–0.86 wt %). Olivine in type 1 dunite enveloping podiform chromitite has Li concentrations and δ7Li values of 1.48–1.71 ppm and 6.19 ‰– 7.98 ‰, respectively. Type 2 dunite lenses with Cr-rich spinel contain olivine with Li = 0.98–1.64 ppm and δ7Li = 6.77 ‰–10.99 ‰. The type 3 dunite lenses with Cr-poor spinel show the highest values of Li = 0.94– 1.40 ppm and δ7Li = 10.25 ‰–14.20 ‰. Exsolution lamellae of clinopyroxene and magnetite occur as oriented intergrowths in olivine of type 3 dunite lenses with Cr-poor spinel. We suggest that the Purang ophiolite developed during two main stages of formation. In the first stage, abyssal peridotites formed in a mid-ocean-ridge environment. During the second stage, hydrous high-Mg boninitic melts were produced by high degrees of partial melting in a supra-subduction zone mantle wedge, which reacted with peridotite to form type 2 dunite pods with high-Cr# spinel. At lower degrees of partial melting in the same mantle wedge, Al-rich melts were produced, which reacted with peridotite to form type 3 dunite pods that contain lowCr# spinel. These Al-rich melts were also relatively rich in Ti4+, Ca2+ and Fe3+, which were incorporated into the olivine structure by appropriate substitutions. During cooling, these elements exsolved as lamellae of magnetite and clinopyroxene. Published by Copernicus Publications on behalf of the European mineralogical societies DMG, SEM, SIMP & SFMC. 188 F. Xiong et al.: Multistage origin of dunite in the Purang ophiolite, southern Tibet


Introduction
The Li isotope system is of particular significance in distinguishing diverse sources of melt in subduction zones; it is also very sensitive to different types of dehydration reactions and metamorphic settings (Chen et al., 2000). The 7 Li can be preferentially released from the subducting plate, resulting in 7 Li isotope enrichment in the overlying mantle wedge and 6 Li isotope enrichment in the residue (Tomascak et al., 2002). Release from the slab to the mantle wedge causes changes in both the quantity and isotopic composition of Li. Although there is overlap, Li isotope compositions can also be of help in distinguishing different sources of formation. Eclogite and granulite, in particular, as well as altered mid-ocean-ridge basalt (MORB) and marine sediments, for instance, may show much lower values of 7 Li compared to ocean-island basalt (OIB), MORB or arc lava (Tomascak et al., 2002). However, Li isotope variations have been studied far less in residual mantle peridotite, although such information is essential for our understanding of partial melting processes in subduction zones . Previous work has shown that Li isotope compositions can be used to trace the process of rock melting reactions (Su et al., 2015Xiao et al., 2017). Such an approach is thus of particular interest for understanding the origin of the different types of dunite investigated in the current study. Based on the Li isotope results, we can expect important information for tracking and distinguishing various subduction-related processes and the interaction between melt and olivine to understand the extent of mantle heterogeneity in the Purang peridotite.
In this contribution, we describe and present mineralogical and Li isotope data on three types of dunite associated with chromitite deposits of the Purang ophiolite in Tibet to interpret their formation.

Geological setting
The Purang massif (Fig. 1) is the largest ophiolite sequence in the western part of the Yarlung-Zangbo suture zone and consists mainly of a large mafic to ultramafic body extending for ∼ 600 km 2 , Xiong et al., 2018. The ultramafic rocks include mainly harzburgite and minor lherzolite and dunite, all of which contain irregularly shaped and deformed bodies of chromitite (Xiong et al., 2018;Fig. 1b). The main fabric of the peridotites is directed northwestsoutheast. Individual dikes of pyroxenite, gabbronorite, dolerite and rare basalt locally crosscut in a northwesterly direction of the peridotitic mantle lithologies (Fig. 1b;Liu et al., 2015). The ophiolite tectonically overlies a late Cretaceous mélange along its southern boundary that includes some volcanic rocks and metamorphosed, shallow to bathypelagic carbonate rocks and shales of Triassic age (Xiong et al., 2013, Fig. 1). In the southwest, the ophiolite is bounded by a fault zone that contains rocks of a serpentinite mélange and is thrust over Triassic carbonaceous slates of the Tethyan Himalaya in the southeast (Fig. 1b).
U-Pb laser ablation, multi-collector, inductively coupled plasma mass spectrometry (LA-ICP-MS) dating of zircon crystals recovered from the dolerite and gabbro dikes within the peridotites has yielded crystallization ages ranging from 140 to 120 Ma (Miller et al., 2003;Li et al., 2008;Liu et al., 2011;Chan et al., 2015;Xiong et al., 2019b;Meng et al., 2016Meng et al., , 2019Zheng et al., 2019). Thus, the magmatic age of the Purang massif is considered to be Early Cretaceous. Major geotectonic sources of the Purang ophiolite are connected with a detachment fault at a slow-spreading ridge (Liu Figure 1. (a) Simplified tectonic map of the Tibetan Plateau, composed of several terranes and separated by well-defined suture zones. The (Indus-)Yarlung-Zangbo suture zone is the southernmost of these sutures, separating the Lhasa block to the north from the Indian plate to the south. Abbreviations: BNS, Bangong-Nujiang suture; IYZS, Indus-Yarlung-Zangbo suture. Locations of other UHP-mineral-bearing ophiolites along the YZSZ follow Xiong et al. (2019a). (b) Detailed geological map of the study area indicated in (a) (modified from Xiong et al., 2013;Dai et al., 2013;Liu et al., 2014Liu et al., ). et al., 2014, providing a continental-margin ophiolite influenced by plume-type magmatism , recycled subcontinental lithospheric mantle  and a sub-oceanic lithospheric mantle wedge from an underlying downgoing slab Su et al., 2015).

Sampling and analytical methods
Our regional-scale geological fieldwork on the Purang ophiolite focused on the different rock types of the massif. Harzburgite and three types of dunite could be distinguished based on their contrasting appearance in the field (Fig. 2 -169-4.2 L-169-4.3 L-169-4.4 L-169-4.7 L-169-4.9 L-169-4 Fig. 1b). The dunite lenses with Cr-rich spinel are randomly distributed over the entire massif, while the dunite forming envelopes around chromitite occurs both in the northwestern and southeastern parts of the massif (Fig. 1b).
The microtextural characteristics of the Purang peridotite were studied with the backscattered electron (BSE) technique using an FEI Nova Versa-3D scanning electron microscope (SEM) equipped with an energy-dispersive spectrometer (EDS) at the Key Laboratory of Deep-Earth Dynamics of the Ministry of Natural Resources, Institute of Geology, Chinese Academy of Geological Sciences. The operating conditions for the SEM were 20 kV excitation voltage and 15 nA beam current (Co metal used for calibration). Olivine, pyroxene and Cr-spinel were analyzed with a JEOL JXA-8100 electron microprobe (EMP) at the Key Laboratory of Deep-Earth Dynamics of Ministry of Natural Resources, Institute of Geology, Chinese Academy of Geological Sciences. All analyses were carried out in the wavelength-dispersive Xray spectroscopy (WDS) mode. Analytical conditions were 15 kV accelerating voltage, 20 nA beam current, and 2 µm spot diameter with a counting time of 20 s on the peaks and 10 s on the background.
Transmission electron microscopy (TEM) was used to determine the composition and symmetry of the mineral inclusions in olivine; sample preparation was performed using the focused ion beam (FIB) technology from MNR Key Laboratory of Metallogeny and Mineral Assessment, Institute of Mineral Resources, Chinese Academy of Geological Sciences. Typical TEM foils prepared had the dimensions 15 µm ×10µm ×0.20 µm. A FIB single beam device (FEI Versa 3D) operated at CAGS was used for sample preparation. A Ga ion beam (30 keV acceleration voltage) focused onto a selected location of the sample surface was used to sputter material from the sample. For details of the TEM method used here, see Wirth et al. (2009).
The in situ Li isotope measurement on olivine was performed using a Cameca IMS-1280 Secondary ion mass spectrometry (SIMS) from the Key Laboratory of Mineral Resources, Institute of Geology and Geophysics, Chinese Academy of Sciences. The primary O ion beam accelerated at 13 kV with an intensity of about 15 to 30 nA. The beam spot was elliptical, with a size of about 20 µm ×30 µm. Positive secondary ions were measured on the ion multiplier in pulse-counting mode with a mass resolution (M/DM) of 1500 and an energy slit opening of 40 eV without any energy offset. A total of 180 s of pre-sputter without grating was applied prior to analysis. The secondary ion beam position in the contrast aperture and the magnetic field and en-ergy offset were automatically centered prior to each measurement. The measurement was carried out for 30 cycles, with a counting time of 12, 2, and 4 s for 6 Li, a background sample of 6.5 mass, and 7 Li, respectively. Using an olivine sample 09XDTC1-24OL with a forsterite component of 94.2 as standard , the similar composition of the standard as compared to the analyzed olivine eliminates any possible contribution of the matrix to the data obtained . The measured δ 7 Li values are given as δ 7 Li ([( 7 Li/ 6 Li) sample /( 7 Li/ 6 Li) L-SVEC − 1] × 1000), relative to units of the standard National Institute of Standard and Technology (NIST) SRM 8545 (L-SVEC). The instrumental mass fractionation is expressed as δ 7 Li = δ 7 Li SIMS − δ 7 Li MC-ICPMS . A total of 22 analyses of the standard confirmed a homogeneous Li isotopic composition of i = 26.0±1.9 ‰ (2SD) (Su et al., 2015). The external 2σ errors of the isotope compositions for both standard and sample are less than 2.5 ‰.

Type 1 dunite envelopes around chromitite
Most of the podiform chromitites in the Purang ophiolite are surrounded by dunite envelopes that vary from a few centimeters to several meters in thickness (Fig. 2b). The thickest and best-developed ones occur around disseminated chromitite. This type of dunite is moderately fresh, light yellowish-brown in color and massive. Typical samples consist of 98 vol % olivine and 2 % magnesiochromite (Fig. 3b); neither clinopyroxene nor orthopyroxene has been identified in these rocks, not even as inclusions in olivine or chromite. Most olivine grains are subhedral, form granular crystals, and show grain sizes between 2 and 6 mm. They rarely exhibit deformation structures such as kink banding or undulatory extinction. Spinel occurs as small euhedral to subhedral grains (Fig. 3b) and as larger (up to 5 mm), irregular, amoeboid crystals, some of which contain small inclusions of euhedral olivine.

Type 2 dunite lenses with Cr-rich spinel
This type of dunite lens is enclosed in harzburgite, forming very fresh and mainly brown rock inclusions (Fig. 2c). In general, the lenses are oriented northwest-southeast and show a width of 1-5 m and a length of 5-10 m. The rock is mainly composed of olivine (95-99 modal %, Fig. 3) and contains small amounts of Cr-rich spinel, orthopyroxene and clinopyroxene. The olivine is large in grain size, mainly 0.5-3 mm in length. This type of olivine does not contain any inclusions or exsolution textures. Locally, the rims of olivine are replaced by serpentine. Chromian spinel has a particle size of between 0.2 and 1 mm.

Type 3 dunite lenses with Cr-poor spinel
Type 3 dunite lenses and pods in the Purang peridotites are mostly 0.1-2 m in width and commonly aligned in a northwest-southeast direction, roughly parallel to the main foliation of the host harzburgites (Fig. 2d). The dunite is fresh and nearly unaltered (< 2 vol % serpentine) with brownishblack to yellowish-brown colors. It consists mainly of olivine (95-99 modal %) and orthopyroxene (2 %-4 %) with minor amounts of clinopyroxene and Cr-poor spinel. Olivine grains are commonly 2.5-3 mm in diameter, but recrystal-lized grains are much smaller. Some of the largest grains show undulatory extinction and deformation twinning. Crpoor spinel, 0.2-1 mm in size, is mostly randomly distributed in the rocks but can locally concentrate to form small, disseminated clusters.

Magnetite and diopside exsolution lamellae in olivine
Magnetite and diopside exsolution lamellae were observed in thin sections of type 3 dunite lenses with Cr-poor spinel and confirmed by electron microprobe (EMP), SEM and TEM analysis. These lamellae (Figs. 3d, 4-6) can be oriented in different planes of the olivine structure (see below). Magnetite lamellae are generally 1-4 µm wide and show a length of 5-50 µm. Diopside lamellae have the same width but are typically a little shorter (maximum length 30 µm). The two minerals are commonly intergrown, as shown in Fig. 5, where a single exsolution lamella is composed of magnetite (lower part) and diopside (upper part). Volumetrically, diopside exsolution lamellae are much less abundant than magnetite lamellae.
Because they are very small and commonly superimposed, it is difficult to obtain precise compositions. However, microprobe analyses show that exsolution lamellae of diopside have compositions that are nearly identical to irregularly distributed inclusions of diopside in olivine grains (Table 4). Some relatively large diopside inclusions appear to   have grown contemporaneously with the magnetite exsolution lamellae. Figure 4 also shows elemental X-ray maps of olivine grains with exsolved diopside and magnetite.
6 Major element and isotope chemistry 6.1 Major element chemistry of olivine Figure 7a and b compare Fo vs. NiO and MnO, respectively. Figure 7a shows that olivine compositions from harzburgite and the three types of dunite overlap to a large degree. However, forsterite contents in olivine type 1 dunite envelopes and dunite lenses with Cr-rich spinel reach somewhat higher values of Fo 91.8 and Fo 92.2 , respectively, than harzburgite (Fo 91.2 ) or dunite lenses with Cr-poor spinel (Fo 90.7 ).
In Fig. 7b the olivine compositions of the different source rocks are presented in a MnO vs. Fo component diagram, again documenting similar compositions of olivines in harzburgites and type 3 dunite. Olivine of dunite that envelopes chromitite lenses shows a higher Fo component; the respective values of olivine from type 2 dunite lenses are even higher.

Li isotope chemistry of olivine
Lithium isotope compositions, the respective amounts of forsterite component of olivine, as well as selected chemical parameters of accompanying Cr-rich spinel of the harzburgite, and the three different types of dunite are listed in Table 2. Generally, the Li concentrations and δ 7 Li values of the olivines studied show large variations from 0.75 to 1.74 ppm and 6.19 ‰ to +14.72 ‰, respectively, and they are negatively correlated (Fig. 8). Note that the olivines in harzburgite exhibit the lowest amounts of Li (0.8 ppm on average; range of 0.75-0.89 ppm) and the highest concentrations of δ 7 Li (average of 13.8 ‰; range of 11.82-14.72 ppm), whereas those from type 1 dunite enveloping chromitite lenses have the highest amounts of Li (1.6 ppm on average; range of 1.48-1.74 ppm) and the lowest δ 7 Li concentrations (average of 7.1 ‰; range of 6.19 ‰-7.98 ‰). Li contents of olivines from type 2 and type 3 dunite lenses vary widely between 0.94 and 1.64 ppm, but some olivines of the type 2 dunite lenses have considerably higher amounts of Li (Fig. 8). The δ 7 Li values of olivine in the type 3 dunite lenses are in general higher (10.25 ‰-14.2 ‰) than those of the type 2 dunite lenses (6.77 ‰-10.76 ‰). In Fig. 9 a series of correlation diagrams display the amount of Li (in ppm) and the δ 7 Li values (in ‰) in olivine relative to the Fo component in olivine and the NiO and MnO content, as well as Cr# of the spinels in the different lithological types studied. These diagrams indicate that dunite formed in different environments and/or stages. Note, however, that olivine from harzburgite and type 3 dunite lenses with Cr-poor chromite exhibit quite similar chemical characteristics.

Chemical composition of the spinel supergroup minerals
Spinel supergroup minerals (Bosi et al., 2019) are significant indicators that allow us to distinguish the different host rocks and gain insight into possible geotectonic provenance. In Fig. 10, discrimination diagrams were chosen to characterize the spinel compositions; the respective analyses are given in Table 3. Figure 10c, a Cr# [100 × Cr/(Cr + Al)] vs. Mg# [100 × Mg/(Mg + Fe 2+ + Mn)] diagram, documents a compositional trend with increasing Cr# and Mg# from type 3 dunite, via harzburgite and type 1 dunite enveloping chromitite, to type 2 dunite. Magnesiochromite in type 2 dunite lenses, both in the matrix and as inclusions in olivine, is typically small and euhedral and shows Cr# between 66 and 67 and Mg# between 52 and 54, whereas spinel. in the type 3 dunite lenses yields Cr# between 19 and 20 and Mg# between 72 and 78. The diagram also documents that only the spinel supergroup minerals of the type 2 dunites are magnesiochromite, whereas those of the three other lithological types are spinel. The analyses of magnesiochromite from the type 2 dunite lenses fall into the compositional field of ophiolitic chromitites (Kamenetsky et al., 2001) and are located between the two fields of abyssal and forearc peridotites (Fig. 10c).
The Cr# vs. TiO 2 diagram (Fig.10b) documents low amounts of TiO 2 (< 0.1 wt %) in spinel supergroup minerals of the type 1 dunite that envelopes chromitite and in the harzburgite; the respective Cr# is about 50 and 30. The type 3 and type 2 dunite lenses contain spinel supergroup minerals with TiO 2 contents between 0.13 and 0.23. The correlation between the Cr# of spinel supergroup minerals and the Fo content of the adjacent olivine reveals a compositional trend within the olivine-spinel mantle array (OSMA; Arai, 1994; Fig. 10a).

Chemical composition of clinopyroxene
Clinopyroxene is essentially diopsidic with compositions between 31.2 and 50.1 mol % En end-member (Table 4). Correlation diagrams for selected oxides are shown in Fig. 11.
In general, the clinopyroxenes in the harzburgite have the lowest Mg#, between 90.7 and 93.0. In the type 1 dunite that envelopes chromitite lenses, the Mg# is between 91.7 and 93.7. In the type 3 dunite lenses the Mg# is similar (95.3-96.0), but in the type 2 dunite lenses Mg# is considerably higher and varies between 91.9 and 93.0. The Al 2 O 3 content of clinopyroxene is variable, lying in the range of 4.28 wt %-5.45 wt % in the harzburgites, 1.56 wt %-3.29 wt % in the type 1 dunite enveloping chromitite, 0.78 wt %-0.86 wt % in the type 2 dunite lenses and 3.48 wt %-5.24 wt % in the type 3 dunite lenses.
The contents of Cr 2 O 3 in clinopyroxene of the peridotites range from 0.1 wt % to 1.32 wt %, and Na 2 O contents are below 1.0 wt %. The clinopyroxenes of harzburgites show the highest amounts of Cr 2 O 3 , Na 2 O and TiO 2 (0.16 wt %-0.25 wt %) compared to the remaining peridotites. The amount of TiO 2 in clinopyroxene is also relatively high in the type 3 dunite lenses (0.12-0.19), whereas it is significantly lower in the type 1 dunite enveloping chromitite and in the type 2 dunite.

Origin of different dunite types
The bulk rock characteristics and Li isotope ratios of the different types of dunite from the Purang ophiolite, as well as their corresponding mineral compositions, are quite different. The type 3 dunite lenses contain chromian spinel with the lowest Cr# observed, which may have experienced interaction with an aluminum-rich melt. In addition, exsolution textures were produced in type 3 dunite during uplift Fig. 10a). Note that the composition of chromian spinel and olivine from the type 3 dunite is similar to that in abyssal peridotite (Zhou et al., 2005;Arai and Miura, 2016). Note also that the type 2 dunite lenses and the type 1 dunite enveloping chromitite differ from the type 3 dunite. The type 1 dunites were possibly more intensely affected by an interaction with melt. A similar behavior was reported by Zhou et al. (2005), who suggested that the melt in a supra-subduction zone setting would be of a boninite composition, whereas the type 1 dunite envelope is considered to be formed by the reaction of a subduction-related melt with lherzolite or wehrlite (Zhou et al., 2005;Xiong et al., 2015;Zhang et al., 2019a, b). In an olivine vs. chromian spinel mantle-array diagram (Arai, 1994), the type 3 dunites and the harzburgites plot entirely within the area defined by abyssal peridotite, with approximately 8 % and 10 % partial melting (Fig. 10a), indicating MOR peridotite. In contrast, type 1 dunite envelopes and type 2 dunites, which can be estimated to have experienced about 20 % and 35 % partial melting, respectively, show compositional features similar to subduction zone (SSZ) peridotite (Fig. 10a), which may be the result of a multistage development.
The formation of high-Cr chromitites, dunites and refractory harzburgites is thought to result from varying degrees of reaction between a boninite-like melt and Cpxbearing harzburgite (Zhou et al., 2005). In addition, most of the dunite envelopes are found only around disseminated chromitites, not the massive varieties (Fig. 2b). The prevailing view that podiform chromitites form only in SSZ settings by reaction between hydrous boninitic melts and mantle peridotites in the shallow mantle (Zhou et al., 1996(Zhou et al., , 2005Uysal et al., 2005;Xiong et al., 2015) is based on the study of large chromitite bodies. Kelemen et al. (1990) provided reaction experiments between a melt of olivine tholeiite and harzburgite at 5 kbar and 1050-1150 • C and demonstrated that the amount of SiO 2 in the solid phase gradually increases with the consumption of the harzburgite. The Mg/(Mg + Fe) ratio in the solid phase also decreases until the harzburgite is finally consumed. Their study showed that a similar process might be expected during interaction with depleted mantle peridotites, e.g., when an arc of magma rises through the mantle lithosphere of the overriding plate. This process involves extensive reaction events between melt, olivine, orthopyroxene and spinel. Sakae and Jun-Ichi (2007) pointed out that through this boninitic melt-rock reaction a large number of different dunite lenses can be produced in a peridotite massif in New Zealand. The type 2 dunite lenses with Cr-rich spinel form intercalations within harzburgite, and the compositions of olivine, clinopyroxene and chromian spinel show gradual changes from harzburgite. Therefore, the type 2 dunite lenses may be the remnants of a mantle peridotite after high partial melting, and the type 3 dunite with Cr-poor spinel could result from an in situ reaction between an aluminum-enriched melt and peridotite. The type 1 dunite envelopes are probably a result of a melt-peridotite reaction.
Several authors have described clinopyroxene and magnetite intergrowth structures in olivine and interpreted them as exsolution features (e.g., Arai, 1978;Moseley, 1984;     and Trinity (Lundstrom et al., 2005) ophiolites and the Gakkel Ridge peridotites  plotted for comparison. The δ 7 Li ranges of MORB, OIB, arc lava, eclogite + granulite, altered MORB and marine sediment are from Tomascak et al. (2016) and references therein. The two arrows, "partial melting" and "fractional crystallization" are defined based on the Trinity data (Lundstrom et al., 2005). Xiong et al., 2017;Liang et al., 2014Liang et al., , 2018. Ashworth and Chambers (2000) and Ren et al. (2008), however, interpreted them as symplectites. The structures observed in the Purang dunites of the present study are here interpreted as exsolution lamellae. Exsolution of magnetite from olivine could presumably represent either a decomposition product of a solid solution between olivine and laihunite [Fe 2+ Fe 3+ 2 (SiO 4 ) 2 ] or an oxidization product of Fe 2+ in the olivine (Moseley, 1984).
Olivine accommodates extremely small quantities of Ca 2+ and Na + ions in its crystal lattice due to the larger ionic radii of these elements compared to those of Fe 2+ , Mg 2+ and Mn 2+ . Olivine is thus not expected to exsolve clinopyroxene. Experimental studies on the olivine-melt system show that the higher the contents of CaO, Na 2 O and FeO in the melt, the higher the Ca content in the olivine (Libourel, 1999). In this study, olivine grains with exsolution microstructures point to a relatively high CaO content prior to exsolution (Fig. 4). The CaO content in the olivine cannot indicate fractional crystallization, nor can it be due to to the effects of oxygen fugacity, temperature, and pressure (Libourel, 1999). Thus, the variation in CaO content of olivine in the dunite from the Purang ophiolite (Fig. 4) may be related to later modification by melt or fluid. The clinopyroxene in the rodlike intergrowths of magnetite in the type 3 dunite olivine grains is aligned in the same direction as the needle-like intergrowths (Figs. 4,5), indicating their simultaneous exsolution. It has been documented that clinopyroxene can crystallize from hydrothermal fluids at temperatures higher than 750 • C (Python et al., 2007;Akizawa et al., 2011;Akizawa and Arai, 2014). The association of clinopyroxene and magnetite and the presence of interstitial clinopyroxene in the Purang ophiolite suggest the involvement of hydrothermal fluids during the formation of the dunite. Hence, the reactant Arrows represent the magma differentiation trend with high partial melting (Lundstrom et al., 2005;Su et al., 2016). Symbols are the same as in Fig. 7. fluids could have originated from hydrothermal fluids. Reequilibration between olivine and fluids at high temperatures could produce solid solutions with Ca 2+ and Na + in olivine. This feature has been well documented in dunites of the Luobusa ophiolite in the eastern part of the Yarlung-Zangbo suture zone (Ren et al., 2008). Because of the rapid uplift and cooling of the Purang peridotite (Liu et al., 2014), the exsolution lamellae in the olivine grains do not show plastic deformation; thus, their primary exsolution features are well preserved.

Li isotope fractionation of the different types of dunite and harzburgite
The isotopes of Li, 6 Li and 7 Li can be easily fractionated because of their large mass difference (∼ 15 %), resulting in large δ 7 Li variation in the terrestrial environment. For example, the δ 7 Li value of seawater is ∼ +32 ‰ and that of N-MORB and MORB sources ranges from ∼ +1.5 ‰ to +6.5 ‰, whereas altered MORB has a larger range from −2 ‰ to +14 ‰ (Chan et al., 1999;Decitre et al., 2002;Lundstrom et al., 2005;Vils et al., 2008). During subduction, Li in subducted slabs can be easily mobilized and carried into the overlying mantle wedge by fluids (Ishikawa et al., 2005;Savov et al., 2007). Li isotopic signatures developed in the crust and mantle can be used to study mantle heterogeneity. Abyssal peridotites have Li concentrations of < 3 ppm, which are lower than those of rocks formed by a melt-peridotite reaction. Based on the dissolution of olivine during the process of partial melting, Li isotopic diversities of the different peridotites studied can be seen as the result of different amounts of melt extraction and high-temperature hydrothermal alteration Xiao et al., 2017) (Fig. 8).
The different types of olivine-rich rocks discovered in the Purang ophiolite formed due to different mechanisms. In olivine the Fo end-member content correlates with the different Li isotope characteristics in the various peridotites (Fig. 9). The Purang peridotites have heavy Li isotopes with high δ 7 Li values of 6.19-14.72, which are similar to arc lavas but also overlap with values found in metamorphic products such as eclogites and granulites, as well as altered MORBs and marine sediments (Fig. 8). Note that the type 1 dunite enveloping chromitite has higher amounts of Li than the harzburgite but the lowest values of δ 7 Li. Since Li is highly mobile, it is sensitive to low-temperature fluid-rock interaction (Chan et al., 2002). Lithium and (preferentially) 6 Li move from olivine grains to the hydrothermal fluid during serpentinization, and as a result olivine becomes Li-depleted and isotopically heavier in Li (Decitre et al., 2002;Lundstrom et al., 2005;Wimpenny et al., 2010;Wunder et al., 2010;Su et al., 2016). In addition, crustal contamination may lead to an increase in the Li content in the magma and a change of its 7 Li isotopic composition, especially during magma ascent and ophiolite emplacement (Teng et al., 2004). The Purang peridotite is high in MgO and rare-earth elements (REE) Xiong et al., 2018), and the Li content in olivine is very low ( Table 2). This means that the process of crustal contamination is negligible. Teng et al. (2011), Xiao et al. (2016 and Bai et al. (2017) demonstrated the possibility of an intense exchange of chemical components between olivine and chromian spinel. It is widely accepted that Fe diffuses from olivine to chromite as the temperature decreases (Mg does the opposite), thus increasing the Fo component of olivine (Su et al., 2015;Xiong et al., 2015;. Previous studies have suggested (e.g., Lundstrom et al., 2005;Tomascak et al., 2016) that because 6 Li diffuses slightly faster than 7 Li inside minerals, a negative anomaly of about 10 ‰-20 ‰ should be produced in olivine. This diffu-  Ozawa, 1994;Pagé et al., 2008;Dick and Bullen, 1984).  Dick and Bullen (1984), Bonatti and Michael (1989), Hébert et al. (1990), Johnson et al. (1990), Ishii et al. (1992), Arai (1994), Pearce (1998), Pearce et al. (2000) and Parkinson et al. (2003). The degrees of melt extraction calculated based on spinel composition in (a) are based on the empirical formula of Hellebrand et al. (2001). The arrows and thick dashed line in Fig. 10b show interaction between MORB-like melt and a supra-subduction zone melt (Pearce et al., 2000). Symbols are the same as in Fig. 7. Figure 11. Compositional ranges of selected oxides for clinopyroxene from different lithologies of the Purang ophiolite. Data for abyssal peridotite and fore-arc peridotite are from Pagé et al. (2008). sion mechanism may explain the slightly negative correlation between Li and δ 7 Li observed in the Purang olivines of the different ultrabasic rocks (Fig. 8), but it does not explain the variability of the entire dataset of 10 ‰ of δ 7 Li difference and its correlation with magma differentiation indicators between the different dunites. The decreased δ 7 Li trend of olivine in dunite was explained as a result of Li diffusion from melt in a melt conduit to the host peridotite during melt-harzburgite interaction .
Chemically or isotopically driven intragranular diffusion may affect the Li concentration of a mineral, depending on whether its growth occurred from a fluid or a melt (Tomascak et al., 2016). Even if the Li isotopes take more time to achieve compositional uniformity between melt and crystals than the Li abundance (Richter et al., 2014), the Li content in Purang peridotites is similar, whereas the δ 7 Li values have some unclear variation (Fig. 9). The δ 7 Li values of forsterite of the type 2 dunite lenses and the type 1 dunite enveloping chromitite are lower than those of the harzburgite and the type 3 dunite lenses. The same features are evident in NiO and Fo of olivine and Cr# of spinel, respectively, which show a negative correlation with the δ 7 Li values (Fig. 9).
Lithium is a moderately incompatible element during mantle melting processes and magma differentiation (Brenan et al., 1998). The concentration of Li in the magma increases as the degree of differentiation increases (e.g., Hamelin et al., 2009;Marks et al., 2007;Weyer and Seitz, 2012;Ackerman et al., 2015). The Purang olivine shows a tendency towards increased Li with a decrease in Fo component (Fig. 9).
In addition, the Li content of olivine is negatively correlated with NiO and also negatively correlated with Cr# in spinel supergroup minerals (Fig. 9). These trends are consistent with trends that result from magma differentiation. Marks et al. (2007) found evidence that olivine fractionation drives the derived magma to move in the direction of high Li and that this elevated Li trend is thought to be related to the redox conditions of the magma. The diffusion process alone does not explain the preferred enrichment of δ 7 Li in the melt relative to olivine phenocrysts. The Purang peridotite acts as a quasi-stable aquifer system that can carry large amounts of fluids that were introduced during the SSZ subduction (Xiong et al., 2018). A large number of studies have demonstrated that Li, in general, is an element that can easily transport via a fluid phase (Brenan et al., 1998;. With respect to questions related to isotope fractionation processes, Chan et al. (1999) were able to show that 7 Li becomes preferentially enriched in the fluid phase. Su et al. (2015Su et al. ( , 2016 have demonstrated that during the subduction of the basaltic oceanic crust and overlying sediments, fluids may become released and percolate through the overlying mantle wedge; magma that is generated in such an environment can thus incorporate the Li isotope characteristics of these fluids.

Tectonic setting of the Purang ophiolite
The Purang ophiolite is one of the numerous ophiolite bodies located along the Yarlung-Zangbo suture zone in southern Tibet, which are considered to be either remnants of the Neo-Tethyan lithosphere emplaced during the collision between the Indian and Asian plates or as ancient subcontinental lithospheric mantle domains . It consists chiefly of mantle peridotites locally accompanied by gabbros, dikes and solidified lavas. The solidified lavas range from arc tholeiites to boninitic and OIB-like basalts. Most of the peridotites have "U-shaped", chondrite-normalized REE patterns, which indicate the re-fertilization of depleted peridotites in a supra-subduction zone environment (Xiong et al., 2018).
However, the presence of both type 2 and type 3 dunite pods in the Purang peridotites and the presence of magnetite + clinopyroxene exsolution lamellae in olivine of the type 3 dunite provide important clues to the tectonic environment in which this ophiolite formed. There is strong evidence for a multistage origin of the Purang peridotites. After having reached relatively shallow levels in the upper mantle, the peridotites were trapped in a supra-subduction zone wedge, where they underwent interaction with SSZ melts and fluids (Xiong et al., 2019c). Initially, the hydrous melting of the peridotite produced Mg-rich magmas, which reacted with the peridotites to form dunite pods characterized by Mg-rich olivine and high-Cr spinel. In some ophiolites, such as Luobusa, the dunite pods are closely associated with podiform chromitites. A second stage of melting in the Purang peri-dotites produced boninitic magmas via lower degrees of partial melting, which were relatively rich in Ti 4+ , Al 3+ , Ca 2+ , Na 1+ and Fe 3+ . During the reaction of these melts with the peridotites, small amounts of these elements may have become incorporated into the newly crystallized olivine. During subsequent cooling, they very likely became exsolved as lamellae of magnetite (containing Fe 3+ , Ti 4+ ) and clinopyroxene (containing Al 3+ , Ca 2+ and Na 1+ ). Thus, our preferred explanation for the origin and evolution of the Purang ophiolite is that the mantle peridotites formed in a MOR environment and were then modified by later-stage melts and fluids in a SSZ setting during intra-oceanic subduction along the Indus-Yarlung-Zangbo suture (see Zhou et al., 1996Zhou et al., , 2005Malpas et al., 2003;Hébert et al., 2003;Dubois-Côté et al., 2005;Bédard et al., 2009;Liu et al., 2010;Xu et al., 2011;Xiong et al., 2013;Su et al., 2019). The discovery of a large number of water-bearing minerals in the Purang chromitite is a further indicator of fluid exchange processes (Xiong et al., 2013(Xiong et al., , 2018. The above-described characteristics and discussed processes for the Purang ophiolite are consistent with the geochemical behavior of Li isotopes in a subduction-zone setting, as described by Marschall et al. (2007), Wunder et al. (2006) and Yamaji et al. (2001), who found evidence for a progressive release of fluid from the subducted plate that migrated into the overlying mantle wedge and caused a lowering of δ 7 Li.

Conclusions
1. We recognized three types of dunite in the Purang ophiolite: type 1 dunite, which envelops podiform chromitite and contains Cr-rich spinel with Cr# = 41.5-46.8; type 2 dunite lenses containing magnesiochromite with Cr# = 66.9-67.9; and type 3 dunite lenses with spinel that shows Cr# between 19.8 and 20.6. Olivine in the type 3 dunite exhibits exsolution lamellae of diopside and magnetite; such features are absent in type 2 dunite. The type 3 dunite probably formed by a reaction between harzburgite and an arc tholeiite magma that was relatively rich in Ti, Fe 3+ , Al, Ca and Na, whereas the type 2 variety formed from interaction with boninitic melts containing high amounts of Mg and Cr. Due to rapid cooling during tectonic emplacement of the rocks, the exsolution textures in olivine were produced and preserved; the compositional zoning observed in the spinel supergroup minerals also reflects a formation that occurred during uplift.
2. The Purang ophiolite contains rocks that have different amounts of Li and that also differ in δ 7 Li compositions. The Li isotopic variations in the various types of dunite are a result of different types of occurrence, i.e., as primary dunite at mid-ocean ridges or as dunite that experienced subsequent modification in a supra-subduction zone setting. Olivine from the harzburgite documents a negative correlation between the forsterite component of olivine and δ 7 Li values due to the diffusive influx of Li from interconnected melts. The δ 7 Li isotopic compositions of the dunites are similar to arc-like melts. However, there is also an overlap with MORB, OIB, or even Li that originates from sedimentary sources. The interacting melts leading to crystallization of magnesiochromite are inferred to have been depleted in Li content and enriched in 6 Li relative to 7 Li, thus producing values inferred for highly dehydrated slabs. These different types of melts continuously and progressively reacted with the oceanic lithospheric mantle, resulting in the formation of dunite and chromitite and, in addition, accounting for the observed Li isotope heterogeneity.
Data availability. The data of this paper can be accessed from reliable public data repositories, by assigning digital object identifiers (DOIs) https://doi.org/10.5194/ejm-32-1-2020 (Miyawaki, 2020) and properly citing datasets as individual contributions, such as the journal of European Journal of Mineralogy.
Author contributions. FX, JY, ZL and XX undertook the fieldwork, collected the rock samples and acquired the different analytical data. FX, JY and HPS shaped the work's conception, studied the microscopic data and SEM features, and interpreted the elemental and isotope data. FX and HPS wrote the manuscript and finalized the submission. All authors have approved the submitted version.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "To Christian Chopin, in recognition of 30 years of dedicated service to EJM". It is not associated with a conference.