Corundum-bearing and spinel-bearing symplectites in ultrahigh-pressure eclogites record high-temperature overprint and partial melting during slab exhumation

. Corundum-and spinel-bearing symplectites after muscovite were found in ultrahigh-pressure (UHP) eclogites from the Dabie terrane, China. Three types of symplectites were recognized based on their mineral assemblages: (1) symplectitic intergrowths of corundum + plagioclase + biotite after phengite (CPB), (2) sym-plectitic intergrowths of spinel + plagioclase + biotite after phengite (SPB), and (3) symplectitic intergrowths of spinel + plagioclase after paragonite (SP). The microtextures and mineral assemblages of the symplectites, in combination with the results of thermodynamic modeling on local regions, indicate that these symplectites formed by the breakdown of phengite and paragonite during the granulite-facies metamorphic overprint (770– 850 ◦ C) of the eclogite at pressures of 0.8–0.9 GPa. Dehydration partial melting reactions occurred during the breakdown of muscovite, which leads to the formation of thin plagioclase ﬁlms (silicate melts) along grain (garnet, rutile, quartz) boundaries. Mass balance calculations indicate that the development of CPB and SPB symplectites after phengite requires the introduction of Al, Ca, Na, and Fe and loss of Si, Mg, and K. However, the formation of SP symplectites after paragonite requires the input of Mg, Ca, and Fe and removal of Si, Al, and Na. By summarizing the occurrence and growth mechanism of corundum-and spinel-bearing symplectites in global UHP terranes, we ﬁnd that such symplectites can form by both the subsolidus replacement of an Al-rich anhydrous mineral (kyanite) and the dehydration melting of an Al-rich hydrous phase during high-temperature metamorphism. This study reveals that muscovite-bearing eclogites may experience multiple episodes of partial melting during the slab exhumation, not only at the high-pressure (HP) exhumation stage but also at the lower-pressure metamorphic overprinting stage. Kyanite is a reaction product during the HP partial melting, whereas the low-pressure (LP) melting will consume kyanite. We propose that the occurrence of corundum-and spinel-bearing symplectites after muscovite in eclogites is a potential mineralogical indicator of LP melting in exhumed slabs.


Introduction
Partial melting of subducted crust is an important petrologic process, which could modify the rheology of the slabs (e.g., Rosenberg and Handy, 2005;Labrousse et al., 2011), produce adakitic and tonalite-trondhjemitegranodiorite (TTG) magmas (e.g., Goss and Kay, 2006;Palin et al., 2016), and trigger element transfer and chemical differentiation at convergent plate margins (e.g., Hermann et al., 2006;Brown et al., 2011;Stepanov et al., 2014;Ferrero et al., 2016;Guo et al., 2022a). Thus, documenting the partial melting of subduction-zone metamorphic rocks is of significant importance for understanding the lithospheric rheology, thermal structure, and process of material cycling between the shallow and deep regions of the Earth. However, these rocks generally have experienced complex metamorphic evolution histories during subduction-exhumation processes, involving the prograde stage, the HP-UHP (high-pressure and ultrahigh-pressure) peak stage, and multiple stages of retrogression (e.g., Malaspina et al., 2006;Groppo et al., 2015). Thus, it is essential to precisely identify the potential melting events at different stages and determine the mechanisms and P -T (pressure and temperature) conditions of the partial melting reactions.
Numerous studies have reported the occurrence of partial melting of eclogites in global UHP metamorphic belts, such as those in the Dabie-Sulu orogen, the Kazakhstan Kokshetau massif, and the Western Gneiss Region of Norway (Gao et al., 2017, and references therein). The partial melting events of eclogites are identified based on both the microtextural evidence, such as the felsic films and cuspate minerals around eclogitic matrix grain boundaries, the multiphase solid inclusions (MSIs) in eclogite minerals, and the glassfilled pores in refractory minerals (e.g., Holness and Sawyer, 2008;Zeng et al., 2009;Chen et al., 2014;Gao et al., 2017;Wang et al., 2017), and also the outcrop-scale phenomena of melt flow, such as felsic segregations, eclogitic migmatites, and pegmatites veins (e.g., Liebscher et al., 2007;Song et al., 2014;Wang et al., 2014). Investigations on the above textures and migmatites reveal that the partial melting of the eclogites mostly occurs during the initial exhumation stage at HP conditions (1.5-2.0 GPa) (e.g., Liu et al., 2009;Gao et al., 2012;Tang et al., 2021). However, the temperatures may further increase when these eclogites exhumed to the crustal depth (pressure < 1.5 GPa) (e.g., Gordon et al., 2012;Deng et al., 2019;Cao et al., 2021). Whether the partial melting occurs at lower-pressure conditions during slab exhumation is still largely unknown due to a lack of diagnostic petrological records.
In this study, we report distinct corundum-and spinelbearing symplectites in the UHP eclogites from the Dabie terrane, China. We present comprehensive investigations of petrology, Raman analyses, mass balance calculations, and thermodynamic modeling to constrain the formation mechanism and P -T conditions of these symplectites. Our results reveal that the corundum-and spinel-bearing symplectites formed by the partial melting of phengite and paragonite during the granulite-facies overprinting. This study indicates that the UHP eclogites have experienced multiple episodes (HP vs. LP) of partial melting during exhumation. We also summarize and compare the occurrence, condition, and formation mechanism of corundum-and spinelbearing symplectites in eclogites from global UHP terranes and stress that the occurrence of such symplectites replacing muscovite in eclogites provides important implications for fluid or melt actions in exhumed slabs. Mineral abbreviations used in the text, figures, and tables are based on Whitney and Evans (2010).

Geological background and sample description
The Dabie-Sulu metamorphic terrane, located in the eastern portion of central China collision belt (Fig. 1a), formed by the subduction of the South China Block (SCB) beneath the North China Block (NCB) in the Triassic (e.g., Zhang et al., 2009a). The Dabie terrane is located in the southwestern part of the Dabie-Sulu UHP terrane, which is separated by ∼ 500 km of the sinistral Tan-Lu strike-slip fault. The occurrence of micro-diamond and coesite in metamorphic minerals suggests that the South China Block had been subducted to 120 km depths at least (e.g., Okay et al., 1989;Xu et al., 1992). The Dabie terrane is generally subdivided into six metamorphic zones from north to south ( Fig. 1a): (1) the north Huaiyang low-grade metamorphic zone, (2) the north Dabie metamorphic core complex, (3) the coesite eclogitefacies metamorphic zone, (4) the quartz eclogite-facies metamorphic zone, (5) the amphibolite-facies metamorphic zone, and (6) the greenschist-to blueschist-facies metamorphic zone (Zheng et al., 2005). The eclogite-facies metamorphic zones mainly consist of various types of gneisses, which are intruded by numerous Cretaceous granites (Fig. 1a) (Faure et al., 2003;Zhang et al., 2009a). Variable-size blocks of eclogite, garnet peridotite, marble, jadeite quartzite, peridotite, and marble can be found in the gneisses.
The present study area is located in the Yanwo area in the Dabie UHP terrane (Fig. 1b), where UHP eclogites occur as various sizes of lenses and blocks (a few centimeters to over hundreds of meters in length) within granitic gneisses (Fig. 1b). A new petrological study (Tang et al., 2021) indicated that the Yanwo UHP eclogites had experienced a clockwise P -T metamorphic path with peak metamorphic conditions of 2.8-3.1 GPa and 630-700 • C. Moreover, both metamorphic dehydration (∼ 2.2-2.7 GPa and 620-750 • C) and partial melting events (∼ 1.5-2.0 GPa and 700-750 • C) during exhumation were recognized in the Yanwo eclogites and HP veins (Tang et al., 2021).

Analytical methods
Backscattered electron (BSE) images were acquired using the Nova Nano scanning electron microscope (SEM) 450 at the Institute of Geology and Geophysics, Chinese Academy of Sciences (IGGCAS). The operation conditions are as follows: accelerating voltage of 15 kV, a beam current of 6 nA, and working distance of 6 mm.
Fine-grained minerals (< 1 µm) in symplectite assemblages were identified using a WITec alpha 300R confocal Raman microscope with a 532 nm excitation laser at the IG-GCAS. The detail positions of peaks in Raman spectra, referring to the RRuff online database (https://rruff.info/, last access: 15 June 2023), are as follows: 377, 419, 573, and 642 cm −1 for corundum; 406 and 666 cm −1 for spinel; 402, 487, 901, and 947 cm −1 for kyanite; 197, 331, 509, 677, and 1001 cm −1 for orthopyroxene; 180, 285, 481, and 508 cm −1 for plagioclase; and 221, 285, 400, 599, and 1350 cm −1 for hematite. For single spectrum analyses, the integration time and accumulations are set as 10 and 2 s, respectively. The size of the laser spot was 0.4 µm, and the laser power was 2-10 mW. Raman mapping was set at the step of 0.3-0.5 µm with an acquisition time of 3-5 s and one accumulation. The accuracy of Raman line shifts was checked by regular measurements of the 521 cm −1 signal for a synthetic monocrystalline Si standard.
Compositions of minerals were analyzed using a CAMECA SXFive electron probe microanalyzer (EPMA) at the IGGCAS. The accelerating voltage was 20 kV, and the beam current was 100 nA for fine-grained minerals (orthopyroxene, corundum, spinel, and kyanite) with a 1 µm beam diameter. The operating conditions for other minerals (plagioclase, muscovite, biotite, garnet, clinopyroxene, and amphibole) were 15 kV accelerating voltage, 18 nA beam current, and 5 µm beam diameter. The natural minerals and synthetic oxides used for calibration are diopside (Si, Ca, and Mg), albite (Na), orthoclase (K), rutile (Ti), Al 2 O 3 (Al), FeS 2 (Fe), Cr 2 O 3 (Cr), MnO (Mn), ZnS (Zn), and NiO (Ni). A program based on the ZAF procedure was used for data corrections of minerals (Armstrong, 1995). The analytical uncertainties for most major elements are less than 2 % relative. The X-ray maps (Si, Al, Fe, Ca, Mg, Na, and K) of the symplectite regions were acquired by wavelength dispersive spectroscopy (WDS) analyses using the CAMECA SXFive EMPA at the IGGCAS. The operation conditions were 15 kV accelerating voltage, 150 nA beam current, and a dwell time per point of 45 ms. The X-ray maps were produced at a 1024 × 1024 resolution.

Petrography
Three types of symplectites after muscovite (phengite and paragonite) were recognized according to the mineral assemblages  in the studied eclogites.
The first type of symplectite is the intergrowth assemblages of corundum + plagioclase + biotite after phengite (CPB). Figure 4 shows a typical CPB assemblage with 10 vol %-20 vol % of corundum, 55 vol %-65 vol % plagioclase, and 20 vol %-25 vol % of biotite . The corundum occurs as needlelike crystals ranging from 5 to 100 µm in length and has an intergrowth direction the same as the matrix of plagioclase ( Fig. 4c-g). The biotite typically shows a tabular shape ( Fig. 4c-k). Anhedral kyanite is occasionally found near the CPB symplectite, and the crystal size of kyanite is much bigger than the symplectitic minerals (k) Sketch map of (j).
( Fig. 4a-c). The CPB symplectites not only distribute around the phengite but also develop in the interior regions of the phengite where small fractures develop (Fig. 4j, k). Moreover, the CPB symplectites in the phengite display an internal porosity with sizes of irregular voids down to a few tens of nanometers (Fig. 4j, k).
The second type of symplectite is the intergrowth assemblages of spinel (15 vol %-20 vol %) + plagioclase (60 vol %-65 vol %) + biotite (15 vol %-20 vol %) after phengite (SPB). The minerals in the SPB symplectite typically show subhedral to anhedral shapes. Spinel is generally vermicular with a size of 1-5 µm in the matrix of plagioclase ( Fig. 5a-c). The growth of spinel is perpendicular to the phengite. Biotite forms tabular and subhedral crystal and has its long axes parallel or subparallel to the orientation of spinel. Several quartz grains appear in the symplectite surroundings (Fig. 5d, e). However, the spinel grains are never found to be in direct contact with quartz. Rather, spinel is separated from quartz by plagioclase. Kyanite is not observed in the SPB symplectites. In addition, we find that a single phengite grain can be replaced by both CPB and SPB symplectites (Fig. 4a). Both types of symplectites coexist with each other (Fig. 5f).
The third type of symplectite is the intergrowth assemblages of spinel + plagioclase after paragonite (SP). Figure 6a shows a typical SP symplectite, which is consisted of 25 vol %-35 vol % spinel and 60 vol %-70 vol % plagioclase. This type of symplectite develops after both matrix paragonite ( Fig. 6a, b) and paragonite inclusions in omphacite (Fig. 6c,d). Amphibole occasionally occurs in the SP symplectites (Fig. 6c). The Raman mapping (Fig. 3f), BSE images ( Fig. 6b, d), and X-ray maps (Fig. 7) show that the spinel grows perpendicularly to the paragonite. The grain size of spinel (up to 20 µm) in the SP symplectites after paragonite is relatively larger compared to that in SPB symplectites after phengite. Omphacite in contact with the SP symplectites shows a compositional variation (Fig. 6a, c). In addition, the orthopyroxene (Opx) grains after omphacite are found to occur adjacent to the SP symplectite (Fig. 6d).
Interestingly, a distinct texture of plagioclase films was found near the three types of symplectites (Fig. 8). In the CPB ( Fig. 8a-d), SPB (Fig. 8e, f), and SP ( Fig. 8g) symplectite domains, many fine films of plagioclase can be observed to distribute along the boundaries of the symplectites and other minerals (e.g., Qz and Rt). The films show finely tapered ends, suggesting consistently small dihedral angles (Fig. 8b,d,e,g). In the symplectite domain, plagioclase with clear cuspate margins extends along the symplectite boundaries, and the surrounding mineral grains develop an embayed shape ( Fig. 8e-g). Locally, the plagioclase, intergrown with biotite, occurs at the junction between garnet grains and forms the films separating these garnet grains (Fig. 8h).

Mineral chemistry
Most of symplectitic minerals have a very small grain size. In this study, the minerals whose grain sizes are larger than 1 µm were analyzed due to instrument limitation. Representative compositions of minerals are listed in Tables 1 and 2 and Tables S1-S3 in the Supplement.

Amphibole and kyanite
The matrix amphibole is classified as edenite according to Leake et al. (1997) with X Mg of 0.68-0.75 and (Na + K) A of 0.72-0.88. The amphibole in the SP symplectites is classified as pargasite or magnesiokatophorite with X Mg of 0.70-0.75 and (Na + K) A of 0.65-0.97 (Table S4). Kyanite is almost a pure Al 2 SiO 5 endmember with minor FeO contents (< 0.73 wt %).

Petrological interpretation of symplectites
The first type of symplectite (CPB) typically develops at the external boundary of coarse-grained phengite or occurs along the internal fractures of phengite (Fig. 4b, c). The mineral assemblages are generally distributed aligning parallel or subparallel to the long axis of the original phengite. This orientation relationship is consistent with the topotactic nucleation of the reaction products within the dehydroxylated phengite lattice (e.g., Brearley and Rubie, 1990;Rodriguez-Navarro et al., 2003;Devineau et al., 2006). Experimental investigations have shown that the thermal decomposition of phengite at high-T (≥ 750 • C) and low-P (< 1.0 GPa) conditions could generate Al-rich oxides as a peritectic phase with characteristic dendritic to skeletal morphologies (e.g., Brearley and Rubie, 1990;Devineau et al., 2006;Masotta et al., 2018), which is consistent with the textures of the skeletal and needle corundum crystals in the CPB symplectites. Similar textures have been also reported in phengite-bearing eclogites from the Black Forest, Germany, and were interpreted as the products of phengite melting by the reaction Ph → Crn + Pl + Bt + other Al silicates + melt during the transformation from a granulite-facies to a HT-LP amphibolitefacies stage (Wimmenauer and Stenger, 1989). Moreover, based on the petrological observation and P -T estimation, Dokukina et al. (2017Dokukina et al. ( , 2021 revealed the polymineral assemblages consisting of Crn + Kfs + Bt + newly formed Ms + dumortierite in phengite-bearing quartz-feldspar rocks from the Belomorian eclogite province (Russia) as a dehydration melting product of phengite by the reaction Ph → Crn + Kfs + Bt + melt under HT granulite-facies conditions. Furthermore, the compositions of corundum in CPB symplectites from the studied eclogites (0.52 wt %-1.01 wt % for FeO) are similar to those of the corundum (< 2.70 wt % for FeO) produced by the experimental melting of phengite (e.g., Brearley and Rubie, 1990;Rodriguez-Navarro et al., 2003;Devineau et al., 2006;Masotta et al., 2018). We thus consider that the corundum (Crn)-bearing symplectites in the Dabie eclogites are not formed by the subsolidus reactions of phengite. Rather, they are most likely the partial melting products of phengite. The anhedral shape and serrated grain boundary of the residual phengite ( Fig. 4a) are also consistent with those produced by the partial melting of phengite (Lang and Gilotti, 2007).
Interestingly, kyanite was found to always coexist with the CPB symplectites (Fig. 4b, c) but is never found near the other two types of symplectites, indicating that kyanite is also involved in the reaction of CPB symplectite formation. The grains of kyanite are much coarser than those of the symplectite minerals (Fig. 4c), suggesting that the kyanite is a reactant rather than a product of reaction. The interpretation is supported by many previous studies (e.g., Nakamura and Hirajima, 2000;Janák et al., 2012;Moulas et al., 2015), which reveal that kyanite could be replaced by the symplectites of Crn + Pl during HT metamorphism (≥ 700 • C). Thus, we conclude that the reaction responsible for the CPB symplectite formation is as follows: Ph + Ky → Crn + Pl + Bt + melt.
The second type of symplectite (SPB) also formed by the consumption of phengite (Fig. 5a, b). Vermicular spinel, rather than acicular corundum, occurs in the SPB symplectites. Several experimental studies reveal that phengite can be replaced by spinel, biotite, and other K-Al phases (mullite, feldspar, and K-feldspar) during the decomposition of phengite (e.g., Patiño Douce and Harris, 1998;Rushmer, 2001;Michaud et al., 2021). Importantly, silicate melts and glass were generally observed surrounding the spinel grains in these experiments (e.g., Patiño Douce and Harris, 1998;Rushmer, 2001;Devineau et al., 2006), which clearly indicates the formation of spinel + biotite assemblages by partial melting. The compositions of spinel (25.1 wt %-27.0 wt % for FeO; 9.66 wt %-10.8 wt % for MgO) in the SPB symplectites from the Yanwo eclogites are similar to those of spinel (23.7 wt %-32.8 wt % for FeO; 7.28 wt %-12.1 wt % for MgO) formed by the partial melting of phengite (e.g., Brearley and Rubie, 1990;Patiño Douce and Harris, 1998;Spicer et al., 2004), supporting the interpretation that the SPB symplectites are also the products of phengite melting. Similar symplectites of SPB in eclogites are also reported in other UHP terranes, such as the Bohemian Massif (Faryad et al., 2006); Papua New Guinea (Faryad et al., 2019); and Tromsø Nappe, Norway (Janák et al., 2012), and were proved to form at high temperatures of 700-850 • C at 0.5-1.5 GPa, which are beyond the fluid-absent solidus in metabasalt system (Holland and Powell, 2001). Different from the CPB symplectites, no kyanite is found near the regions of the SPB symplectites, implying that the involvement of kyanite during the breakdown of phengite facilitates the formation of corundum instead of spinel. Thus, the replacement reaction of SPB symplectites can be concluded as follows: Ph→ Spl + Pl + Bt + melt. The third type of symplectite (SP) develops around both the matrix paragonite and paragonite inclusions in omphacite. The symplectitic assemblage of spinel + plagioclase ± amphibole after paragonite has rarely been found in eclogites so far (Ferrando et al., 2005). Both the matrix paragonite and paragonite inclusions show obvious corroded structures (Fig. 6b-d). The Raman (Fig. 3f) and X-ray maps (Fig. 7) show that fine-grained vermicular spinel (up to 10 µm) grows perpendicular to the original paragonite grain   boundary, indicating the growth of spinel at the expense of paragonite. The omphacite in contact with paragonite commonly shows embayment textures and serrated grain boundaries, suggesting that omphacite grains were partly consumed and may be involved in the reaction of SP symplectite formation (Fig. 6a-c). For example, the lower Mg, Ca, and Fe compositions of omphacite adjacent to the SP suggest that MgO, CaO, and FeO from omphacite were added into the products (Fig. 9e, f; Table S3). Thus, the reaction of SP formation can be concluded as follows: Pg + Omp→ Spl + Pl ± Amp + melt.  Table 5.
In addition to these symplectites, several textures of partial melting and melt flow were observed in the eclogites (Fig. 8): (1) abundant fine films of plagioclase are observed to be distributed along the boundaries of the CPB and SPB symplectites with other minerals (e.g., quartz and rutile) (Fig. 8a-f). The films show typical tapered ends, which suggest consistently small dihedral angles (Fig. 8b, d). This textural feature is consistent with the presence of former silicate melts in HT metamorphic rocks (e.g., Brearley and Rubie, 1990;Holness et al., 2011;Vernon, 2011;Bartoli et al., 2013;Dokukina et al., 2021). The plagioclase moats are interpreted as textures of hydrous melt passing along the grain boundaries (Holness, 2011; Bartoli et al., 2013;Dokukina et al., 2021). (2) In the SP symplectite domain, plagioclase with clear cuspate margins extends along the symplectite boundaries (Fig. 8g). The cuspate-shaped areas of plagioclase have been regarded as the pools of crystallized melt that has penetrated along the surrounding minerals (e.g., Sawyer, 2001;Holness and Sawyer, 2008). The embayed shape of surrounding minerals (e.g., Grt and Amp) could be caused by the corrosion during partial melting, according to Vernon (2011). (3) In other regions, the intergrowth of plagioclase with fine-grained biotite distributes at the junction between garnet grains (Fig. 8h). This texture is similar to the interstitial pockets of plagioclase (crystallized from a silicate melt) at quartz triple junctions in partially melted quartz-plagioclase rock (Holness, 2011). Thus, these textures, in line with the three types of symplectites, again support the interpretation that hydrous melt has been released during the breakdown of a hydrous phase (phengite or paragonite) into an anhydrous assemblage (corundum, spinel, and plagioclase) during slab exhumation.
Considering that the three types of symplectites formed by the partial melting of muscovite, the occurrence of corundum + plagioclase + biotite + melt defines the P -T conditions of CPB symplectites after phengite at > 0.7 GPa and > 770 • C according to the pseudosection (Fig. 10a). The formation of SPB symplectites after phengite was constrained in the P -T range of < 0.95 GPa and > 770 • C (Fig. 10b). The assemblage of SP is stable in the P -T range of < 1.0 GPa and > 770 • C (Fig. 11). In addition, the isopleths of X Mg in spinel and X Ca (= Ca/(Ca + Na)) in plagioclase were calculated to further constrain the P -T conditions of symplectites. Based on the CPB assemblages and the X Ca values in plagioclase Figure 10. P -T pseudosection for the symplectites after phengite calculated in the NCKFMASHO system. The rock composition is given as molecule percentage (mol %) oxide. The green, red, yellow, and orange lines represent the stability curves of kyanite, phengite, paragonite, and solidus, respectively. X Mg isopleths of spinel and X Ca isopleths of plagioclase are shown. The blue regions show the P -T fields of the symplectites. (a) P -T pseudosection for the CPB symplectites after phengite. (b) P -T pseudosection for the SPB symplectites after phengite.  Figure 11. P -T pseudosection for the symplectite after paragonite calculated in the NCKFMASHO system. The rock composition is given as molecule percentage (mol %) oxide. The green, yellow, and orange lines represent the stability curves of kyanite, paragonite, and solidus, respectively. X Mg isopleths of spinel and X Ca isopleths of plagioclase are shown. The blue regions show the P -T field of the symplectites.
In summary, the three types of symplectites formed at similar HT conditions (0.8-0.9 GPa and 770-850 • C). It is thus reasonable to consider that these symplectites after phengite and paragonite developed at a same metamorphic stage during the exhumation of eclogites. Moreover, Opx-bearing symplectites after omphacite were found to coexist with the SP symplectites (Fig. 6d). Therefore, the P -T conditions of Opx-bearing symplectites were calculated using the Cpx-Opx solvus thermobarometers of Brey and Köhler (1990) and Taylor (1998). The lower limit of pressures for the eclogitefacies conditions in the Yanwo eclogites is ∼ 1.5 GPa (Tang et al., 2021). Thus, temperatures at pressures ranging from 0.5 to 1.5 GPa were estimated. In terms of the compositions of clinopyroxene (Di 56.1−59.5 Jd 3.8−6.7 Ae 0.9−5.3 ) and orthopyroxene (X Mg values of 0.67-0.76), five groups of Cpx-Opx pairs yielded temperature ranges of 772-837 • C (Brey and Köhler, 1990) and 763-879 • C (Taylor, 1998), respectively (Table S5), which are in good agreement with those constrained by the phase equilibrium modeling. Therefore, we conclude that the formation of Crn-bearing and Splbearing symplectites after muscovite, as well as the development of Di + Opx + Pl symplectites after omphacite, all occur at the granulite-facies stage during exhumation.

Element exchange during symplectite formation
The Crn-bearing and Spl-bearing symplectites were formed by the replacement of phengite or paragonite. If muscovite is regarded as a reference, such a replacement process can be considered as occurring in open systems and involving element exchange or transfer. Quantifying the element exchange during the formation of symplectites is essential to understand the replacement processes. In this study, a mass balance method of Godard and Mabit (1998) was used to constrain the element gains or losses from the system.
The formation of the phengite pseudomorph (CPB and SPB) can be expressed as follows: CPB: 1 Ph +X i → a 1 Crn +b 1 Pl +c 1 Bt; SPB: 1 Ph +X i → a 2 Spl +b 2 Pl +c 2 Bt. Paragonite pseudomorph can be expressed as SP, 1 Pg +X i → a 3 Spl + b 3 Pl, where a 1 , a 2 , a 3 , b 1 , b 2 , b 3 , c 1 , and c 2 are stoichiometric coefficients and X i is a vector (X Si , X Al , . . . , X Na ), representing the quantity of the element i (Si, Al,..., Na) transferred between the muscovite system and its surroundings. If X i = 0, the system is closed for the considered element i; otherwise, X i represents the gain (X i > 0) or the loss (X i < 0) of the element.
According to Godard and Mabit (1998), X i can be expressed by the following equation for each considered element i (taking the equation of the CPB symplectite for an example): where Fv is the volumetric factor defined as ( V/V) +1, and V M (e.g., V Crn , V Spl , V Pl , V Bt , and V Ph ) represents the mo-  land and Powell (1998). X Crn vol , X Spl vol , X Pl vol , and X Bt vol are the volumetric proportion of the products, which are obtained by estimating the mineral distributions (Table 4). N M i is the quantity of the element i in the mineral M, which is obtained by EPMA. In our calculation of mass transfer, eight elements (Si, Al, Mg, Fe, Ca, K, Na, and O) are considered, and then eight independent equations (α) with nine unknowns variables (X Si , X Al , X Mg , X Fe , X Ca , X K , X Na , X O , and Fv) are presented. There is probably an infinity of solutions because of the under-determined eight independent equations. It is commonly accepted that mass transfer needs to consume energy (Nakamura and Hirajima, 2000). Thus the most likely solution can be decided by minimizing the total mass transfer i |X i |. After gaining the optimal Fv, the X i can be calculated successively. The detailed calculation process has been described in Godard and Mabit (1998). The optimal Fv, X i , and volumetric proportion of the products are listed in Table 4. A similar calculation process is applied to the SPB and SP symplectites. The results are listed in Table 4 and shown in Fig. 12 ( The results indicate that the formations of both the CPB and SPB symplectites require the gain of Al, Ca, and Na from its surroundings and loss of Si, Mg, and K ( Fig. 12a, b). Al, Ca, and Na are most likely derived from the adjacent garnet and omphacite near the symplectites. In contrast, the replacement of paragonite by the SP symplectites needs the external input of Mg, Ca, and Fe and removal of Si, Al, and Na to the surroundings (Fig. 12c). Omphacite is the most likely reactant to provide Mg, Ca, and Fe for the formation of SP symplectites because the omphacite near the SP symplectites shows lower Mg, Ca, and Fe contents than that far away from SP (Fig. 9e, f). The replacement of muscovite by the symplectites and the element transfer are mainly by a dissolution-precipitation metasomatic process (e.g., Anderson and Moecher, 2007;Putnis et al., 2009;Martin and Duchêne, 2015;Guo et al., 2012Guo et al., , 2019Guo et al., , 2022b.  Table 5.

Comparison with Crn-and Spl-bearing symplectites in global UHP terranes
In fact, Crn-and Spl-bearing symplectites after eclogitefacies minerals have been reported in many other UHP metamorphic terranes. Figure 13 and Table 5 summarize the UHP terranes where Crn-and Spl-bearing symplectites have been found, including the Sulu terrane (Nakamura and Hirajima, 2000), the north Dabie complex (Groppo et al., 2015), and the north Qaidam (Zhang et al., 2009b) of China; the Armorican Massif in the Alps (Godard and Mabit, 1998); the Bohemian Massif of Germany (Faryad et al., 2006;Tajčmanová et al., 2007;Štípská et al., 2010); Sveconorwegian (Möller, 1998) and Saskatchewan of Canada (Baldwin et al., 2015); the Rhodope Massif of Greece (e.g., Moulas et al., 2013); the Caledonides of Greenland (Elvevold and Gilotti, 2000); and the Western Gneiss Region of Norway (Janák et al., 2012). Petrological studies indicate that most of these symplectites formed at relatively low-pressure (0.6-1.0 GPa) and high-temperature conditions (> 750 • C) (Fig. 13b), which is in agreement with our estimations. Such a result also demonstrates that the Crn-and Spl-bearing assemblages exclusively formed at the granulite-facies or HT amphibolite-facies stage during the exhumation of eclogites. The widespread occurrence of Crn-and Spl-bearing symplectites in global UHP terranes indicates that most exhumed continental slabs may have experienced a heating stage during decompression, especially for those with large sizes (e.g., Kylander-Clark et al., 2012;Liu et al., 2015;Faryad and Cuthbert, 2020;Wang et al., 2023). The compositions of spinel in the symplectites of SPB and SP fall in the composition range of symplectitic spinel from these previously reported UHP terranes (Fig. 9a). Mean-while, the composition of plagioclase in the three types of symplectite also falls in the composition range of plagioclase in symplectites from these reported UHP terranes (Fig. 9b). However, different from the symplectites in this study, the previously reported symplectites from other UHP terranes were typically formed by the breakdown of Al-rich anhydrous minerals (e.g., kyanite, omphacite, and garnet) (Table 5). This process always produces the anhydrous symplectite assemblages. In contrast, our new results on the Dabie eclogites reveal that the breakdown of Al-rich hydrous minerals (muscovite) can also produce the Crn-bearing and Splbearing assemblages. Importantly, the formation of the symplectites is caused by partial melting, which generated hydrous phases (e.g., biotite) in the products and led to the production and flow of silicate melts at HT conditions. Thus, reactants may play a key role in the formation of different mineral assemblages and the mechanism of metamorphic reactions during the formation of Crn-and Spl-bearing symplectites.
6 Implications for multiple episodes of partial melting of exhumed slab Partial melting of HP-UHP metamorphic rocks plays a key role in the geodynamic evolution and chemical exchange of subduction zones. Our new results demonstrate that the development of Crn-and Spl-bearing symplectites in the UHP eclogites is related to the HT metamorphic overprint and partial melting at LP stages. This result, in combination with Tang et al. (2021), reveals at least two episodes of melting in the Yanwo eclogites during the retrograde process ( Fig. 13a): (1) the first episode of partial melting occurred at the HP stage (∼ 1.5-2.0 GPa and 700-750 • C) and was induced by the breakdown of phengite via the reaction Ph + Qz + H 2 O → Ky + Kfs + melt. Water, which is derived from the breakdown of lawsonite (e.g., Guo et al., 2013Guo et al., , 2015, was involved in the reaction. This process resulted in the development of MSIs in eclogite minerals (e.g., garnet) and some veinlets of Ep + Pl + Kfs in the vein (Tang et al., 2021). (2) As the pressure further decreases to 0.8-0.9 GPa, the second episode of partial melting of phengite and paragonite took place at 770-850 • C, via the following dehydration melting reactions: Ph + Ky→ Crn + Pl + Bt + melt, Ph→ Spl + Pl + Bt + melt, and Pg + Omp→ Spl + Pl ± Amp + melt. This process formed low-pressure Crn-and Spl-bearing symplectites. We thus stress that the eclogites may experience multiple episodes of partial melting (HP vs. LP) during the exhumation process. Multiple episodes of partial melting of HP-UHP rocks were also found in Sulu terrane in China (Wang et al., 2017(Wang et al., , 2023, the Kokshetau in Kazakhstan , and the Caledonides in Greenland (Gilotti et al., 2014). Recognizing and characterizing these different melting events in the eclogites is the key to understand the thermal evolution history of exhumed slabs.  (1994), Moulas et al. (2013) in Greece The occurrence of Crn-and Spl-bearing symplectites replacing muscovite is a potential mineralogical indicator of LP partial melting in exhumed UHP eclogites.
Clearly, kyanite plays a different role in the two episodes of melting. During the HP stage, kyanite can be regarded as a peritectic product of partial melting (Tang et al., 2021), whereas during the LP stage, kyanite is a reactant that is involved in the melting reactions (Fig. 4a-c). This conclusion is also consistent with the experimental results (e.g., Skjerlie and Patiño Douce, 2002;Liu et al., 2009), which indicate that HP partial melting of eclogites (> 1.8-2.0 GPa) typically produces kyanite, whereas the LP (< 1.8 GPa) partial melting will consume kyanite. Thus, clarifying the behavior of kyanite during partial melting, based on the textural characteristics, can be used to distinguish the different episodes of the melting process.

Conclusions
Three types of Crn-and Spl-bearing symplectites in the Dabie UHP eclogites, including CPB (Crn + Pl + Bt), SPB (Spl + Pl + Bt), and SP (Spl + Pl) symplectites, were recognized and studied. The combination of petrology, mineral compositions, and thermodynamic modeling reveals that all the three types of symplectites formed by the breakdown of muscovite via the partial melting at the HT metamorphic stage during the exhumation process (0.8-0.9 GPa) of the UHP eclogites. The CPB and SPB symplectites formed by the partial melting of phengite at 780-850 • C, whereas the SP symplectite formed by the partial melting of paragonite at 770-800 • C. Mass balance calculations indicate that the formations of CPB and SPB symplectites require the gain of extra Al, Ca, and Na into the phengite regions and the loss of Si, Mg, and K from phengite. In contrast, the formation of SP symplectites is accompanied by the gain of Mg, Ca, and Fe into paragonite regions and by the loss of Si, Al, and Na to surrounding minerals. The integration of present results and investigations on Crn-and Spl-bearing symplectites from other UHP terranes demonstrates that such symplectites can be generated by the breakdown of both anhydrous and hydrous minerals. However, the breakdown of hydrous minerals typically produces hydrous phases and releases silicate melts. This study indicates that multiple episodes of partial melting events (HP vs. LP) may occur during the exhumation of a subducted UHP slab. By summarizing the P -T conditions of the Crn-and Spl-bearing symplectites in global UHP terranes, we stress that most of exhumed UHP slabs may experience a common heating stage during ascension to crustal depth.
Data availability. All data derived from this research are presented in the enclosed tables, figures, and Supplement.
Author contributions. SG collected the samples and designed the study. PT carried out detailed petrography, mineral chemistry, phase equilibrium modeling, and mass balance calculations. SG and PT interpreted the data and prepared the manuscript.
Competing interests. The contact author has declared that neither of the authors has any competing interests.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.