绿帘石记录俯冲带变质流体活动*

唐盼 郭顺

1. 中国科学院地质与地球物理研究所岩石圈演化国家重点实验室，北京 1000292. 中国科学院大学，北京 1000493. 中国科学院青藏高原地球科学卓越创新中心，北京 100101

俯冲带是地壳物质循环和地球内部能量交换的最主要场所，维持着地球内部和外部之间长期的物质-能量收支平衡，并最终影响地球表生圈层的宜居环境。大量研究表明板片在深俯冲-折返过程中释放的高压-超高压流体是物质和能量交换的关键介质(McCulloch and Gamble, 1991; Brenanetal., 1998; Manning, 2004; Hermannetal., 2006; Zhengetal., 2011)。这些流体由俯冲板片中含水矿物通过一系列变质反应分解形成，随后在俯冲板片内部发生运移和汇聚，并最终上升进入地幔楔(Peacock, 1990; Schmidt and Poli, 1998; Scambelluri and Philippot, 2001; Hackeretal., 2003; Spandleretal., 2003; Angiboustetal., 2014)。在流体运移的过程中，它将与流经的俯冲隧道内各类岩石发生不同程度的水/岩相互作用，并通过一系列的溶解-沉淀机制进行物质交换，从而不断改变岩石和流体自身的化学组成(Gaoetal., 2007; Zhangetal., 2008; Beinlichetal., 2010; Guoetal., 2012, 2015a, 2019; Johnetal., 2012; Ague and Nicolescu, 2014; Angiboustetal., 2014; Xiaoetal., 2015)。因此，俯冲带流体的形成、运移和交代作用不仅控制着板片/地幔楔部分熔融、中/深源地震活动和板片流变学性质(Davies, 1999; Hackeretal., 2003; Tatsumi and Eggins, 1995; Zhengetal., 2011; Johnetal., 2012; Songetal., 2014; Huetal., 2017; Nietal., 2017; Xiaetal., 2019)，也极大影响着地球内部的元素分配、化学分异和重要矿产资源形成(Scambelluri and Philippot, 2001; Sunetal., 2004; Hermannetal., 2006; Gaoetal., 2007; Zhangetal., 2008; Zheng and Hermann, 2014; Xiaoetal., 2015)。准确限定俯冲带变质流体的形成机制和条件、元素/同位素特征、运移方式/通道/尺度以及流体活动引发的物质迁移是正确认识俯冲带重要地质过程的关键。

由于形成构造环境和温压条件的特殊性，初始形成的深俯冲变质流体在地表条件下难以保存(Manning, 2004)。除了少量被封存在矿物中的原生流体包裹体之外(Andersenetal., 1993; Xiaoetal., 2000; Fuetal., 2001; Klemd, 2004; Zhangetal., 2008; 肖益林等, 2018)，人们很难获取俯冲带流体的直接样本。出露于高压-超高压变质地体的变质岩石是俯冲板片折返到地表的产物，这些变质岩石中发育的含水矿物为我们提供了一个了解俯冲隧道内流体活动的重要窗口。由于俯冲带变质岩石普遍经历了不同期次流体活动的叠加(Lietal., 2001; Zhangetal., 2002; Rubatto and Hermann, 2003; Zhengetal., 2007; Wuetal., 2009; 张泽明和沈昆, 2009; Guoetal., 2012, 2014; Shengetal., 2012; Rubatto and Angiboust, 2015)，选择能够完整记录不同阶段流体事件的含水矿物是准确限定俯冲带多期次流体起源、成分和行为的关键。

绿帘石族矿物是俯冲带变质岩石中一类常见的含水Ca-Al硅酸盐矿物(含水量约为2.0%)(Poli and Schmidt, 1995)。该类矿物既可以在基性岩、酸性岩、泥质岩、碳酸盐岩等多种成分体系中稳定(Brastad, 1985; Spandleretal., 2003; Mattinsonetal., 2004; 梁金龙等, 2006; Liuetal., 2006; Gaoetal., 2007; Johnetal., 2008; van der Straatenetal., 2008; Guoetal., 2012, 2016; Lietal., 2013)，也常常发育在上述岩石的交代岩和变质脉体中(Brunsmannetal., 2000; Enamietal., 2004; Zhangetal., 2005, 2008; Guoetal., 2012, 2014; Brovaroneetal., 2018)。在超高压榴辉岩中，柯石英包裹体的发现表明绿帘石可形成于超高压条件下(Hirajimaetal., 1992; Zhang, 1992; 张泽明等, 2005; Guoetal., 2012)。此外，绿帘石族矿物独特的晶体结构特征使其可以显著富集各类微量元素，特别是二价大离子亲石元素(Sr和Pb)、稀土元素(REE)和过渡族金属元素(TME)等(Franz and Liebscher, 2004; Freietal., 2004)。因此，绿帘石化学成分的变化可以非常敏锐地指示地质过程中流体成分的改变(Martinetal., 2011; Chenetal., 2012; Lietal., 2013; Guoetal., 2015a)。另外，绿帘石族矿物一般具有缓慢的体扩散速率，所以天然产出的绿帘石常常可以保留明显的成分环带(Choo, 2002; 梁金龙等, 2006; Guoetal., 2015b)，这为我们认识不同期次流体活动提供了关键信息。

本文总结了近年来关于俯冲带变质绿帘石研究的最新进展以及该矿物在揭示俯冲带流体活动方面的重要应用，特别是在示踪流体源区，指示流体化学组成和结晶历史，限定流体氧逸度条件以及揭示水/岩交代效应和元素运移方式等方面的应用。阐明了绿帘石族矿物是一个记录俯冲带流体演化和元素迁移的理想矿物，在流体研究方面具有巨大潜力。

图1 不同体系下绿帘石族矿物组合的P-T稳定域(据Poli and Schmidt, 2004修改)(a) CaO-Al2O3-SiO2-H2O体系下的相稳定性；(b)水饱和MORB组分的相稳定性.绿色粗线表示不同绿帘石矿物组合的转变线Fig.1 Stability of epidote-bearing mineral assemblages in P-T projection for different systems (modified after Poli and Schmidt, 2004)(a) phase stabilities in the CaO-Al2O3-SiO2-H2O system; (b) phase stabilities in hydrated MORB compositions. Thick green lines delimit the stability fields of different epidote-bearing assemblages

1 绿帘石晶体化学特征与稳定温压范围

绿帘石族矿物(epidote-group minerals)的名称源自希腊语“epidosis”，含义是“增长”，指的是该类矿物晶形呈(棱)柱状，且(棱)柱的一面比另一面长(Haüy, 1801)。绿帘石族矿物的结构通式为：A1A2M1M2M3(Si2O7)(SiO4)O(OH)，其中A1A2位置以配位数大于6的Ca2+为主(A1配位数是7～9，A2配位数是8～10)，Fe2+和Mn2+可以取代A1位的Ca2+，而Pb2+、REE3+、Sr2+或U4+等元素可以取代A2位的Ca2+；M1M2M3主要以6次配位的Al3+为主，其中Fe3+可以取代M3位的Al3+，Cr3+、V3+或Ti4+可以取代M1M2M3位的Al3+(Deeretal., 1986; Enamietal., 2004)。尽管绿帘石族矿物端元组分的划分和命名非常复杂(Franz and Liebscher, 2004; Gottschalk, 2004)，通常意义上的绿帘石是指单斜晶系的斜黝帘石[Ca2Al2FeSi3O12(OH)]-绿帘石[Ca2Al3Si3O12(OH)]固溶体以及斜黝帘石的斜方晶系同质异像体-黝帘石(Enamietal., 2004; Gattaetal., 2011)。在本文中，如不进行特殊说明，“绿帘石”一词都是指斜黝帘石-绿帘石固溶体和黝帘石。

大量实验岩石学表明，绿帘石的稳定温压范围非常宽广(Liou, 1973; Liouetal., 1983; Poli and Schmidt, 1998, 2004)。在CaO-Al2O3-SiO2-H2O(CASH)简单体系中，绿帘石的稳定温压范围为300～1200℃和0.1～7GPa。其稳定上限主要受控于三个变质反应，分别是：(a)Zo/Ep=Lws+Grt+Ky+Coe、(b)Zo/Ep=Grt+Ky+Coe+H2O和(c)Zo/Ep=Grt+Ky+melt+H2O，稳定下限主要受控于变质反应(d)An+Grt+Crn+H2O=Zo/Ep(图1a)(矿物缩写参见Whitney and Evans, 2010)。同时，相对于水饱和体系，无水体系更有利于绿帘石在高温高压条件下稳定(Naney, 1983; Schmidt and Thompson, 1996; Skjerlie and Patio Douce, 2002; Poli, 2016)。在含Fe体系(CaO-FeO-Al2O3-SiO2-H2O-O2，CFASHO)中，绿帘石的稳定性进一步受氧逸度条件的影响。高的氧逸度条件有利于绿帘石在更高的温压条件下保持稳定(Holdaway, 1972; Liou 1973; Liouetal., 1983)。在水饱和MORB体系中(fO2=NNO)，绿帘石稳定的最大压力可达3GPa，最大温度可达700℃，其形成反应通常与硬柱石分解相关，反应式可表达为Lws+Fe-silicate(如Grt、Cpx等)=Ep+Ky+Coe/Qz+H2O(Poli and Schmidt, 2004; Weietal., 2010; Guoetal., 2013)，随着氧逸度的升高，基性岩岩石中绿帘石的稳定压力上限会进一步提高(图1b)。在天然样品中，Zhang (1992)首次在苏鲁榴辉岩绿帘石中发现柯石英假象，推测绿帘石可在至少3.2GPa和850℃条件下稳定。随后大量的研究证明柯石英包裹体可存在绿帘石中(Hirajimaetal., 1992; Zhangetal., 1995; Rolfoetal., 2000; Yaoetal., 2000; Massonne and O’Brien, 2003; Mattinsonetal., 2004; Guoetal., 2012, 2015a; Liuetal., 2015)。Sobolev and Shatsky (1990)和Zhangetal. (1997)进一步发现绿帘石甚至可以与含金刚石包裹体的石榴石平衡共生，推测其可以在金刚石稳定域存在。因此，实验模拟和天然样品观察均说明在俯冲带变质岩石中，绿帘石可以从俯冲起始阶段(绿片岩相条件)至深度大于90km条件下一直稳定存在。

2 俯冲带变质绿帘石的化学成分

2.1 主量元素

绿帘石的主量组分变化主要涉及M3结构位置的Fe3+-Al替换。因此，绿帘石主量元素变化通常用XFe[=Fe3+/(Al+Fe3+)]值(即Ps值)进行表征。通常情况下，绿帘石-斜黝帘石的XFe值变化范围为0.09至0.30，黝帘石的XFe值变化范围为0至0.06(Enamietal., 2004)。

全岩的Fe/Al比值是影响绿帘石中XFe值大小的一个关键因素。例如，变基性岩石中绿帘石XFe值明显高于变泥质岩石中绿帘石的XFe值。同时，两类岩石的全岩Fe/Al比值均与绿帘石XFe呈正相关(Miyashiro and Seki, 1958; Brunsmannetal., 2002; Poteletal., 2002)。此外，当岩石发生变质相转变时，其全岩Fe含量变化也直接影响绿帘石XFe值。例如在蓝片岩相向榴辉岩相转变过程中，榴辉岩和泥质岩中Fe含量的降低会导致其中发育的绿帘石XFe值明显降低(Ohetal., 1991; Spandleretal., 2003)。

此外，绿帘石XFe值还与温度、压力条件和氧逸度条件密切相关。大量研究表明，不同变质阶段生长的绿帘石具有明显不同的XFe值 (Eide and Liou, 2000; Enamietal., 2004; Mattinsonetal., 2004; 翟伟等, 2006; Guoetal., 2013, 2015a)。然而，不同样品中XFe值表现出不同的变化趋势。苏鲁青龙山和西大别红安榴辉岩中进变质绿帘石XFe值(0.21～0.30)明显高于超高压峰期绿帘石XFe值(0.03～0.19)(Eide and Liou, 2000; Mattinsonetal., 2004; 梁金龙等, 2006)。相反，东大别港河和花凉亭榴辉岩中进变质绿帘石XFe值(0.08～0.15)却明显低于超高压峰期绿帘石XFe值(0.21～0.24)(Guoetal., 2013, 2015a)。这种变化差异表明绿帘石的XFe值可能不仅仅受变质温度、压力因素控制，还可能受其他因素，如氧逸度等影响。Holdaway (1972)和Liou (1973)通过实验证明绿帘石XFe值与俯冲带氧逸度条件呈正相关，即岩石氧化程度越高，绿帘石XFe越大。Guoetal. (2017b)也发现绿片岩相高XFe绿帘石的形成与高氧逸度流体渗透有关。

2.2 微量元素

绿帘石可以显著富集多种微量元素，如大离子亲石元素(Sr和Pb)、锕系元素(Th和U)、稀土元素REE以及过渡族金属(Cr和V)等(Freietal., 2004)。已有的研究表明绿帘石是俯冲带变质岩石中Sr、Pb、Th、U和轻稀土元素(LREE)的最主要载体(图2a)(Brastad, 1985; Spandleretal., 2003; Gaoetal., 2007; van der Straatenetal., 2008; Guoetal., 2012, 2016; Lietal., 2013)。图2b, c分别展示典型榴辉岩和泥质片麻岩中造岩矿物对全岩微量元素的贡献(Spandleretal., 2003)。在榴辉岩中，超过95%的Sr、Pb、Th、U和LREE赋存在绿帘石中；泥质片麻岩中，LREE几乎全部赋存在绿帘石中(>90%)，此外，全岩50%的Sr和25%的Pb赋存在绿帘石中。因此，绿帘石控制着变质作用过程中这些微量元素的分配和迁移。

2.2.1 Sr和Pb

2.2.2 REE

REE能够以不同置换反应进入绿帘石晶格，包括：(1)REE3+和碱金属元素置换A位的Ca2+，如REE3++Na+=2Ca2+；(2)REE3+和Fe2+置换A位的Ca2+和M位的Al3+，如REE3++Fe2+=Ca2++Al3+；(3)REE3+置换A位的Ca2+和Th4+(或U4+)，如2REE3+=Ca2++Th4+/U4+；(4)REE3+和Fe3+置换A位的Fe2+和Th4+(或U4+)，如REE3++Fe3+=Fe2++Th4+/U4+；(5)REE3+与Fe2+和F-置换A位Ca2+和M位Fe3+以及附加阴离子O2-，如REE3++2Fe2++F-=Ca2++2Fe3++O2-(Freietal., 2004; Gieré and Sorensen, 2004)。

图2 绿帘石微量元素特征(a)不同岩性中绿帘石微量元素地球化学特征汇总(数据来源：Brastad, 1985; Gao et al., 2007; Guo et al., 2012, 2016; Li et al., 2013; Spandler et al., 2003; van der Straaten et al., 2008)；基性榴辉岩中(b)和泥质片岩中(c)主要造岩矿物微量元素分布简图(据Spandler et al., 2003修改)Fig.2 Trace elements characters of epidote(a) review of epidote trace elements geochemistry from different rocks (data from Brastad, 1985; Gao et al., 2007; Guo et al., 2012, 2016; Li et al., 2013; Spandler et al., 2003; van der Straaten et al., 2008); trace element distribution or budget from mafic eclogites (b) and from pelitic schists (c) (modified after Spandler et al., 2003)

2.2.3 Cr和V

Cr3+与Al3+电价相同，离子半径相似，因此Cr3+易于取代M1位置的Al3+进入绿帘石晶格(Freietal., 2004; Nagashimaetal., 2019)。绿帘石一般包含10n×10-6～100n×10-6的Cr(Nagasaki and Enami, 1998; Brunsmannetal., 2000)，但在一些特殊地质过程中，绿帘石也可以显著富集Cr2O3(最高可达16%)(Eskola, 1933; Grapes, 1981; Treloar, 1987a; Freietal., 2004; Nagashimaetal., 2006)。这些富Cr绿帘石的形成主要与基性/超基性岩蚀变释放的富Cr流体交代或早期富Cr矿物的分解有关(Grapes, 1981; Yang and Enami, 2003; Tumiatietal., 2005; Nagashimaetal., 2006)。另外，即使在同一岩石样品中，绿帘石Cr含量也可以存在明显差异。如变基性岩中生长在铬铁矿周边的绿帘石比远离铬铁矿的绿帘石具有更高的Cr含量(Nagashimaetal., 2006)。导致这一现象的原因可能是Cr在流体活动中迁移能力有限，难以在岩石中达到含量均一(Treloar, 1987b; Watenphuletal., 2014)。

V3+(0.64Å)的离子半径稍大于Cr3+(0.62Å)，因此V3+主要取代Al3+进入配位数较大的M3位置(Nagashimaetal., 2019)。V在绿帘石中具有很大的相容性(Freietal., 2004)，可以显著富集在绿帘石中。绿帘石V含量与原岩性质有关，如变花岗岩中绿帘石V含量为155×10-6～624×10-6(Spiegeletal., 2002)，而变基性岩中绿帘石V含量为846×10-6～985×10-6(Spandleretal., 2003)。此外，后期富V流体的交代作用也会促使绿帘石V含量增加(Uheretal., 2008)。如出露于坦桑尼亚东北部片麻岩岩石中的 “坦桑石”是一种由富V流体交代形成的宝石级绿帘石，其V2O3含量可以达到0.26%～0.53%(Malisa, 2003; Harrisetal., 2014)。

表1绿帘石-流体稀土元素分配系数

Table 1 REE partition coefficients of epidote-fluid

样品Zo-QzsegregationMF-60MF-61MF-62MF-63Zo30-700Zo30-850P-T条件1.5GPa/500～550℃2GPa/750℃2GPa/800℃2GPa/900℃2GPa/800℃3GPa/700℃3GPa/850℃文献Brunsmannetal.,2001Feinemanetal.,2007Martinetal.,2011分配系数DDσDσDσDσDσDσLa0.081.710.514.961.9515.90.506.611.150.090.052.290.99Ce0.111.970.406.302.6044.44.908.761.620.200.125.061.16Nd0.191.720.395.782.3129.12.307.621.540.060.042.511.14Sm0.501.690.296.622.5125.52.208.092.050.090.063.071.23Eu0.802.480.2510.54.1049.95.4010.62.400.120.073.681.34Ho1.830.1911.13.2024.92.6010.41.900.270.134.340.92Tm2.450.2914.74.4022.02.4011.61.400.420.185.480.83Yb3.640.4025.29.8046.110.5019.82.200.490.185.360.61Lu3.130.4120.26.2025.73.2015.11.400.480.165.530.80

图3 不同P-T条件下绿帘石-流体间稀土元素分配系数1.5±0.2GPa/500±50℃数据来自Brunsmann et al., 2001；2GPa/750～900℃数据来自Feineman et al., 2007；3GPa/700～850℃数据来自Martin et al., 2011Fig.3 REE partition coefficients of epidote-fluid in different P-T conditionsData of 1.5±0.2GPa/500±50℃ from Brunsmann et al., 2001; data of 2GPa/750～900℃ from Feineman et al., 2007; data of 3GPa/700～850℃ from Martin et al., 2011

2.2.4 其他重要元素

除了上述元素，绿帘石同样富集锕系元素Th和U以及变价元素Mn。

Th4+和U4+电价较高，因此Th4+和U4+取代A位置的Ca2+进入绿帘石晶格时，通常伴随其他二价离子取代M位的Al3+。绿帘石Th和U含量主要与原岩性质有关。已有数据显示：变花岗岩和变泥质岩中绿帘石的Th和U含量较高，约是变基性岩中二者含量的10～20倍(Sorensen and Grossman, 1989; Fingeretal., 1998; Spandleretal., 2003; Freietal., 2004; Zacketal., 2004)。另外，后期熔/流体的交代作用以及富Th和U矿物(如独居石、褐帘石和磷钇矿等)的分解也可以导致绿帘石中Th和U含量发生显著改变(Fingeretal., 1998; Guoetal., 2015a)。

Mn有两种价态：Mn2+和Mn3+，其中Mn2+与Ca2+电价相等，且Mn2+离子半径小于Ca2+，导致Mn2+易于取代Ca2+进入绿帘石中配位数较小的A1位置，而Mn3+优先取代绿帘石M3位置的Al3+(Freietal., 2004)。因此Mn3+/Mn2+可能反映体系氧逸度条件(Armbrusteretal., 2006)。Freietal. (2004) 总结发现不同岩石体系中绿帘石的Mn含量差异不明显，如变基性岩中绿帘石Mn2O3含量为0.02%～0.87%，而变质花岗岩中绿帘石Mn2O3含量为0.11%～0.43%。

3 绿帘石在俯冲带流体活动研究中的应用

3.1 示踪俯冲带(多阶段)流体源区

图4 大别山港河地区全岩和不同期次绿帘石的初始Sr同位素投图(据Guo et al., 2014修改)t1=230Ma和t2=210Ma分别代表大别山UHP岩石峰期变质年龄和角闪岩相退变质年龄Fig.4 Initial whole-rock Sr isotopic compositions of Ganghe eclogites, veins, and gneisses as well as the initial in situ Sr isotopic compositions of the various stages of epidote (modified after Guo et al., 2014)t1=230Ma and t2=210Ma represent the ages of UHP metamorphism and amphibolite-facies retrogression for the Dabie UHP rocks, respectively

准确限定俯冲带流体的源区是正确认识这些流体性质和行为的前提。然而，俯冲带岩石中变质流体的源区示踪一直是个难题。最主要的原因是这些岩石中的流体活动往往是幕式的，具有不同成因机制、形成时代和表现形式，如含水矿物生长、变质成脉和退变质阶段的水化作用等等(Santoshetal., 2010; Chenetal., 2012; Guoetal., 2012; Massonne, 2012; Wangetal., 2017)。这些变质流体既可以来自寄主岩石内部脱水(Lietal., 2004; Zhengetal., 2007, 2011; Wuetal., 2009; Guoetal., 2012; Zhaoetal., 2016)，也可以来自寄主岩石之外的其他源区(Austrheim, 1987; Franzetal., 2001; Svensenetal., 2001; Johnetal., 2008; Vrijmoedetal., 2013; Guoetal., 2019)，还可以由内部流体和外部流体叠加形成(Spandleretal., 2011; Lietal., 2013; Guoetal., 2014, 2015b)。要了解俯冲带岩石中流体活动过程，就需要在俯冲带岩石中识别出各个期次的变质流体，并逐一确定这些不同期次流体的来源。

绿帘石通常具有较高的Sr(普遍高于1000×10-6)，有利于进行原位微区的Sr同位素分析。同时绿帘石极低的Rb/Sr比值(通常<0.0001)，使得实测的绿帘石87Sr/86Sr比值与其变质生长时的初始87Sr/86Sr比值十分接近，无须再进行初始比值的计算。Campos-Alvarezetal. (2010)利用飞秒LA-MC-ICP-MS对Sudbury火成杂岩体中的热液绿帘石进行了原位Sr同位素分析尝试，揭示了该岩体中至少存在三期不同来源的岩浆-热液流体活动。Guoetal. (2014)利用LA-MC-ICP-MS对超高压榴辉岩中不同期次绿帘石进行了原位微区Sr同位素测定(图4)，发现榴辉岩相绿帘石具有相同且均一的初始Sr同位素组成(87Sr/86Sr=0.70693～0.70720)，并与全岩Sr同位素组成接近；相反，角闪岩相绿帘石具有更高且明显变化的初始Sr同位素组成(87Sr/86Sr=0.70894～0.71172)，与榴辉岩的围岩片麻岩更接近。因此表明高压-超高压成脉流体来源于榴辉岩本身，而角闪岩相流体来源于外部片麻岩。绿帘石原位Sr同位素分析也被应用于多期复合高压脉体的流体源区识别(Guoetal., 2015a)和中-低压退变流体的源区识别(Guoetal., 2016, 2017b)。Bieseleretal. (2018)在研究洋壳经典剖面阿曼蛇绿岩时，发现位于蛇绿岩序列下部的辉长岩中发育新鲜绿帘石脉体，其利用TIMS测得脉体绿帘石的原位87Sr/86Sr比值为0.70429～0.70512，介于幔源岩浆(87Sr/86Sr=0.7032)和海水(87Sr/86Sr=0.7073)之间，表明绿帘石的形成与海底热液蚀变作用有关，为海洋-地壳相互作用提供可靠的岩石学证据。这些研究都表明绿帘石的原位Sr同位素分析是示踪流体源区的一个有效手段。

3.2 指示含水矿物(硬柱石)分解和流体释放过程

在了解流体源区基础上，准确限定这些流体形成的岩石学过程对理解流体的性质和行为同样重要。已有的研究表明俯冲板片中含水矿物的脱水分解是俯冲带自由流体的主要来源(Pawley, 1994; Schmidt and Poli, 1998; Hermannetal., 2006; Zhengetal., 2009, 2011)。硬柱石[CaAl2Si2O7(OH)2·H2O]是深俯冲板片中一个重要的含水矿物(含水量高达11.5%)，其形成和分解将极大影响俯冲带流体活动、中/深源地震和板片流变学性质(Zacketal., 2004; Davis and Whitney, 2006; Teyssieretal., 2010; Brovarone and Beyssac, 2014)。准确判断硬柱石的出现对限定岩石变质P-T轨迹、流体活动和元素迁移具有重要指示意义。然而，硬柱石极易受折返阶段温度升高和流体渗透的影响而发生脱水分解(Whitney and Davis, 2006)。P-T-X(H2O)相图模拟表明在水不饱和的前提下，随着岩石体系水含量增加，硬柱石会逐渐失稳发生分解(魏春景和崔莹, 2011)。这些因素导致在自然界出露的岩石中硬柱石难以保存(Clarkeetal., 2006; Wei and Clarke, 2011; Guoetal., 2013)。准确识别硬柱石早期是否存在是岩石学家面临的一大挑战。

图5 绿帘石的显微结构和成分特征(a-e)大别山港河榴辉岩中两期绿帘石(Ep-In和Ep-P)结构特征(据Guo et al., 2013)；(f)大别山港河榴辉岩中两期绿帘石(Ep-In和Ep-P)微量元素La-Sr投图(据Guo et al., 2013)；(g、h)西南天山蓝片岩交代岩中绿帘石的溶蚀结构(据John et al., 2008)Fig.5 Textural and compositional characters of epidotes(a-e) textural characters of two-stage epidotes (Ep-In and Ep-P) in Ganghe eclogite (after Guo et al., 2013); (f) plot of La vs. Sr of two-stage epidotes (Ep-In and Ep-P) in Ganghe eclogite (after Guo et al., 2013); (g, h) corrosion structures of epidote in Tianshan blueschist-alteration zone (after John et al., 2008)

在高压-超高压岩石中，绿帘石是硬柱石分解的一个常见产物(Schmidt and Poli, 1994, 1998; Castellietal., 1998; Lietal., 2004; Weietal., 2010)，因此对疑似硬柱石分解形成的绿帘石研究是检验硬柱石出现的一个潜在突破口。El Korhetal. (2009)发现与硬柱石平衡共生的其他矿物(如绿帘石、榍石、磷灰石)具有较低的Sr和LREE含量，而由硬柱石分解形成的绿帘石的Sr和LREE含量明显升高。基于对不同期次绿帘石微量元素变化的认识，Guoetal. (2013)进一步提出质量平衡关系能够更准确地判断榴辉岩中早期是否存在硬柱石。他们在大别港河榴辉岩中识别出两期绿帘石：早期高压绿帘石包体和峰期超高压绿帘石变斑晶(图5a-e)，矿物成分分析显示高压绿帘石的Sr(990×10-6～1890×10-6)和LREE(La为60×10-6～110×10-6)含量远低于超高压绿帘石的Sr(7200×10-6～10300×10-6)和LREE(La为160×10-6～1300×10-6)含量(图5f)。根据超高压绿帘石的体积含量，要求早期绿帘石的含量超过100%，因此，早期高压矿物组合中必然存在除绿帘石之外的一个更加富集Sr和LREE的矿物。根据基性岩体系相平衡模拟结果，在早期可能存在的所有(超)高压矿物中，只有硬柱石能够富集Sr和LREE。这一结果说明早期硬柱石一定存在。因此进变质低Sr和LREE绿帘石的出现可以作为榴辉岩中早期存在硬柱石的一个重要指示标志(Guoetal., 2013)。Lietal. (2013)在西南天山含硬柱石榴辉岩中发现与硬柱石平衡共生的蓝片岩相绿帘石包体Sr和LREE含量明显低于榴辉岩相绿帘石变斑晶Sr和LREE含量，指示硬柱石的分解的确可以引起新生绿帘石的Sr和LREE含量增加。

3.3 记录俯冲带水/岩反应过程及微量元素迁移

俯冲板片中含水矿物脱水释放的初始流体通常具有较低微量元素含量(Gao and Klemd, 2001; Spandleretal., 2003; Volkovaetal., 2004)。这些流体在脱水位置不断的汇聚和迁移，并在运移过程中，与俯冲带隧道内流经的各类岩石发生不同程度的交代作用。在此过程中，流体将通过一系列的溶解-再沉淀机制(ICDP)不断改变岩石和流体自身的化学组成(Bebout, 2007; Gaoetal., 2007; Zhangetal., 2008; Beinlichetal., 2010; Xiaoetal., 2011, 2015; Guoetal., 2012; Klemd, 2012; Angiboustetal., 2014)。因此，水/岩相互作用是控制和影响流体化学特征的关键过程。

图6 大别山港河榴辉岩(样品05GH10、05GH12A-3、05GH12A-2和05GH12A-1)中绿帘石含量、形态和尺寸的空间变化规律(据Guo et al., 2015b修改)Fig.6 A summary sketch showing the spatial varation in abundance, shape, and size of epidote from the Ganghe eclogites (samples 05GH10, 05GH12A-3, 05GH12A-2 and 05GH12A-1) (modified after Guo et al., 2015b)

图7 大别山港河榴辉岩中四个不同空间位置绿帘石(Ep-PI、Ep-PII、Ep-PIII和Ep-PIV)的BSE图、X-ray面扫描和成分剖面图(据Guo et al., 2015b修改)Fig.7 BSE, X-ray mapping and compositional zonation profile showing the compositional character of four types epidotes (Ep-PI, Ep-PII, Ep-PIII and Ep-PIV) in the Ganghe eclogites (modified after Guo et al., 2015b)

已有的研究表明绿帘石的形态、环带结构和微量元素成分可以很好地记录水/岩相互作用过程。Johnetal. (2008)发现西南天山蓝片岩中绿帘石具有相对均一的矿物成分，而榴辉岩相反应带中绿帘石显示了明显的溶蚀结构(corrosion structures)(图5g, h)和变化的矿物成分，表明高压水/岩反应过程中流体部分溶解了绿帘石。Guoetal. (2012)通过对大别山港河含脉体榴辉岩剖面(约20cm)的研究，对超高压水/岩反应过程中绿帘石的行为进行了更加系统的调查。结果表明从远离脉体到靠近脉体，榴辉岩中绿帘石颗粒逐渐变小；绿帘石形态从自形逐渐变为他形；绿帘石含量逐渐降低直至在最靠近的脉体处(约1～2cm)完全消失(图6、图7a-h)。成分分析表明，远离脉体绿帘石的核部成分均匀，逐渐靠近脉体时，绿帘石核部显示补丁状结构(patchy texture)，且绿帘石的Sr、Pb、Th、U和LREE含量逐渐降低(图7)(Guoetal., 2015b)。这些结构形态和成分变化指示超高压水/岩反应过程中，流体溶解了绿帘石并导致微量元素进入流体相。

岛弧岩浆岩相对于亏损地幔楔具有明显高的Sr、Pb、Th、U和LREE，导致这一现象的主要原因是板片来源的流体携带上述元素进入了岛弧地幔源区(McCulloch and Gamble, 1991; Johnetal., 2004)。因此，高压-超高压变质条件下流体-绿帘石反应和溶解作用引发的微量元素迁移，对俯冲带关键元素的循环具有重要意义。

3.4 限定流体的氧逸度条件

流体的氧逸度条件是控制流体性质和行为的一个重要物理化学参数。例如，在俯冲带背景下，较高的氧逸度条件可以显著提高流体对难/微溶元素(如TME)的溶解程度和活化迁移能力，进而影响岛弧玄武岩的地球化学特征(Tatsumi and Eggins, 1995; Franketal., 2002; Mungall, 2002; Sunetal., 2004; Mengasonetal., 2011)。此外，流体的氧逸度还深刻影响着成岩矿物的变质相关系、元素赋存状态和同位素分馏行为(Arculus, 1985; Franketal., 2002; Mattinsonetal., 2004; Caoetal., 2018; Guoetal., 2017b)。因此俯冲带流体氧逸度变化为理解深俯冲变质作用与岛弧岩浆活动具有重要意义。

绿帘石晶格的M3位置主要由Fe3+和Al3+占据，因此绿帘石Fe含量能够灵敏地反映与之平衡流体的氧逸度条件，特别是不同期次绿帘石Fe含量的变化可能记录变质演化过程中流体氧逸度的变化。Donohue and Essene (2000)根据反应式2Ca2FeAl2Si3O12(OH)(绿帘石)=2Ca2FeAl2Si3O12(石榴石)+H2O+1/2O2提出了一个普遍适用于石榴石-绿帘石平衡体系的氧逸度计。Mattinsonetal. (2004)和Caoetal. (2011)根据绿帘石和石榴石平衡关系利用THERMOCALC程序分别对苏鲁榴辉岩中峰期变质流体氧逸度条件和北祁连榴辉岩中高压退变质流体氧逸度条件进行了约束(分别为ΔHM+2.5和ΔFMQ+(2-4))。Guoetal. (2017b)在大别变质地体的退变质脉体中发现铁钛氧化物被金红石+绿帘石+绿泥石组合取代结构，确定取代反应式为24Ilm+8An+8Qz+18H2O+O2=4Ep+4Dph (in Chl)+24Rt，认为该反应涉及铁钛氧化物中的Fe2+被氧化成Fe3+进入绿帘石中，游离出来的Ti结晶形成金红石。通过热力学模拟计算得出引发上述取代反应的绿片岩相退变质流体具有非常高的氧逸度条件(ΔFMQ+(2.5～4.5))。上述这些研究表明利用平衡矿物组合中绿帘石Fe的含量限定流体的氧逸度条件非常有效。

除Fe之外，绿帘石中的微量元素，如Mn、Ce和Eu，也能为流体的氧逸度条件提供重要信息(Grapes and Hoskin, 2004; Guoetal., 2016)。高的氧逸度条件有利于将流体中Mn2+、Eu2+和Ce3+分别氧化为Mn3+、Eu3+和Ce4+，其中Mn2+、Eu3+和Ce3+更容易取代绿帘石中A位置的Ca2+。因此，高氧逸度条件下结晶的绿帘石应具有高的Eu3+含量和较低的Mn3+和Ce4+含量。这一判别标准也与自然界的观察相吻合，如Guoetal. (2014, 2016)发现角闪岩相富Fe绿帘石(形成于高氧逸度条件下)具有明显的负Ce异常。

3.5 揭示流体结晶过程和脉体沉淀次序

流体在运移和汇聚过程中，随着温度压力等物理条件的改变，一部分矿物相会从流体中结晶出来形成脉体。大量的研究发现同一岩体中可以存在矿物组合和全岩成分均不同的脉体，表明成脉过程可能具有多期次性(Zhangetal., 2008; El Korhetal., 2009; Verlaguetetal., 2011; Chenetal., 2012; Guoetal., 2015a)。因此，准确限定脉体的结晶过程，特别是查明脉体的沉淀次序对于理解流体的初始化学组分具有重要意义。

大别造山带花凉亭榴辉岩中发育的复合高压脉体是流体多期次结晶的典型范例。从脉体-榴辉岩的边界到脉体内部，依次出现绿帘石-绿辉石脉、绿帘石-石英脉和绿帘石-蓝晶石-石英脉。成分分析显示，上述脉体中绿帘石的Cr、V和δEu值逐渐降低(分别从115×10-6、363×10-6、1.11降低至7.1×10-6、38.1×10-6、0.51)，而La含量逐渐升高(从2.72×10-6升高至1187×10-6)，反映了绿帘石-绿辉石脉最早结晶，绿帘石-石英脉次之，绿帘石-蓝晶石-石英脉结晶最晚，这一推断与流体成分模拟计算出的不同演化阶段结晶矿物的成分变化一致(Guoetal., 2015a)。同样，Chenetal. (2012)在野寨榴辉岩同一脉体中发现两种不同结构的绿帘石，其中高LREE的绿帘石具有负Eu异常(0.73～0.95)，而低LREE的绿帘石具有正Eu异常(1.03～1.42)，这一化学成分特征暗示前者可能形成于更早期的结晶过程。

4 结语及展望

本文系统论述和总结了俯冲带变质绿帘石在揭示俯冲带流体活动应用中的主要优势：(1)绿帘石稳定温压范围十分宽广，因此能记录板片在深俯冲-折返过程中多阶段的流体演化历史；(2)绿帘石在俯冲带多种岩石体系(特别是交代岩石和脉体)中可以广泛存在；(3)绿帘石的主量元素替代关系简单(Fe3+-Al替代)，能够指示平衡流体的氧逸度条件；(4)绿帘石是微量元素(特别是Sr、Pb、Th、U、Cr、V和LREE)的良好“容器”，其成分变化可以记录不同地质过程中(如水/岩反应过程和结晶过程)流体成分变化和元素迁移；(5)绿帘石具有非常低的Rb/Sr比值，适宜开展原位微区Sr同位素分析和流体源区示踪。

尽管俯冲带变质岩中绿帘石研究已经取得很多进展，但仍存在以下一些问题亟待解决：(1)目前绿帘石的研究主要集中在俯冲带富水流体活动，对于绿帘石在部分熔融和超临界流体活动中的行为认识非常有限。最近的研究发现含绿帘石的多相固相包裹体(熔体/超临界流体的结晶产物)广泛存在于高压-超高压岩石中(Korsakov and Hermann, 2006; Gaoetal., 2012; Chenetal., 2014; Liuetal., 2018)，暗示绿帘石的形成与岩石部分熔融和随后的结晶分异过程密切相关。因此有必要开展熔体和超临界流体中绿帘石行为的研究；(2)虽然绿帘石相对富集Nd，但其Nd的绝对含量通常低于激光原位同位素分析的最低含量要求(约500×10-6)。因此，目前对绿帘石Nd原位微区同位素分析及应用仅限于REE含量较高的绿帘石(如褐帘石)(Romer and Xiao, 2005; Guoetal., 2017a)。随着分析技术的提高，如果能够实现绿帘石激光原位Nd同位素分析，那么绿帘石原位Sr-Nd-(Pb)多元同位素分析的联合应用可能会提供更丰富和准确的流体源区信息；(3)除本文主要论述的绿帘石外，绿帘石族矿物还存在许多复杂的类质同象和同质异像体，如褐帘石、红帘石、钒铬绿帘石等等，它们也常常产出于俯冲带高压-超高压岩石中(Liuetal., 1999, 2003; Hermann, 2002; Bonazzi and Menchetti, 2004; Romer and Xiao, 2005; Uheretal., 2008; Zhangetal., 2008)。对这些特殊绿帘石的系统研究，可能为了解俯冲带不同种类元素循环提供重要信息。

致谢感谢杨岳衡研究员在绿帘石地球化学分析方面给予的大力支持。感谢陈仁旭教授和张聪副研究员审阅了全文并提出建设性的意见。

谨以此文祝贺叶大年院士八十华诞，愿先生健康长寿！