Distributions and Geodynamic Implications of High Field Strength Elements in Rutile from Ultrahigh Pressure Eclogites
-
摘要: 利用LA-ICP-MS对CCSD-MH超高压榴辉岩中金红石进行了详细的原位微区微量元素组成分析.金红石中高场强元素Nb和Ta含量主要受全岩Nb、Ta和TiO2含量控制, Zr、Hf含量比较稳定基本不受全岩含量影响.粒间金红石中, 同一颗粒金红石核部Zr含量系统高于边部, 而边部则出现了明显的Pb和Sr富集特征.CCSD-MH榴辉岩中金红石与全岩的Nb/Ta比值呈现明显的不一致性.全岩Nb/Ta比值明显低于金红石且与全岩TiO2含量负相关, 而金红石的Nb/Ta比值与全岩Nb、Ta含量和Nb/Ta比值没有明显的相关关系.金红石和全岩之间非完全耦合的Nb/Ta组成表明, 金红石并非形成于原岩的结晶过程中而是在超高压变质作用过程中形成, 尽管金红石是榴辉岩中Nb、Ta含量的主要载体矿物, 但金红石的Nb/Ta比值并不一定能完全代表全岩的特征, 而与全岩Nb、Ta和TiO2的含量有关.粒间金红石核部Zr含量所记录的温度与粒径之间具有明显的正相关性, 反映金红石中的Zr在其形成后没有封闭.粒间金红石所表现出的明显的边部富集Pb和Sr的特征, 反映了后期流体活动对金红石组成的影响.这些研究结果为金红石中Zr在高温下的扩散作用和后期流体活动的影响提供了重要证据, 这可能是利用金红石Zr含量地质温度计计算的苏鲁-大别榴辉岩变质温度(598~827℃) 偏低的主要原因.Abstract: Trace element compositions of rutiles in eclogites from the Chinese Continental Scientific Drilling (CCSD) main hole were analyzed using LA-ICP-MS.The results indicate that Nb and Ta contents of rutiles are significantly controlled by whole rock compositions, while Zr and Hf show no obvious dependence on the whole rock compositions.Coupled with enrichments of Pb and Sr at the rim of the interstitial rutiles, Zr contents decrease from the core to the rim.Nb/Ta ratios of rutiles are much higher than those of the parent eclogites.Moreover, Nb/Ta ratios of rutile show no correlation with Nb and Ta concentrations and Nb/Ta ratios of the parent eclogites, while Nb/Ta ratios of eclogites correlate negatively with the TiO2 contents.These observations suggest that rutiles were metamorphic products of the Triassic ultrahigh pressure (UHP) metamorphism, rather than formed by crystallization in equilibrium with melt.The partially decoupling between rutiles and parent eclogites in terms of Nb/Ta ratios implies that although rutiles are the dominant carriers of Nb and Ta in eclogites, they might not completely reflect the Nb/Ta ratios of the parent eclogites, especially for low-Ti eclogites (TiO2 < 0.5%).Nb/Ta ratios of rutiles are dominated by the combination of Nb, Ta and TiO2 contents of the parent eclogites.The positive correlations between temperatures recorded by Zr contents and the size of rutile indicate that Zr in rutile was not immediately closed after rutile formation.On the other hand, the enrichments of Pb and Sr at the rim of the interstitial rutiles suggest influence of subsequent fluids on the compositions of rutiles, which could have affected Zr of the rutiles as well.This study provides the evidence for the Zr diffusion at a high temperature and fluid activity that could have contributed to the lower temperatures (598-827℃) recorded by Zr-in-rutile for the eclogites from the Dabie-Sulu UHP terrane.
-
Key words:
- CCSD-MH /
- eclogite /
- rutile /
- trace element /
- Nb/Ta ratio
-
金红石作为副矿物可以出现在从低级变质岩(如绿片岩) 至高级变质岩(如榴辉岩) 的各种类型的岩石中, 是高场强元素Nb、Ta、Zr和Hf的主要载体矿物之一(Zack et al., 2002).金红石不仅会强烈富集Nb和Ta (DNb、Ta >> 1), 且会使Nb和Ta显著分异(DTa > DNb) (Green and Pearson, 1987; Brenan et al., 1994; Stalder et al., 1998; Foley et al., 2000, 2002; Green, 2000; Horng and Hess, 2000; Schmidt et al., 2004; Klemme et al., 2005; Xiong et al., 2005).因此, 金红石在解释大陆地壳和亏损地幔Nb、Ta含量和Nb/Ta比值不平衡特征、岛弧火山岩微量元组成上受到特别关注(Brophy and Marsh, 1986; Tatsumi and Nakamura, 1986; Morris et al., 1990; Ringwood, 1990; Ulmer and Trommsdorff, 1995; Stalder et al., 1998; Rudnick et al., 2000; Audetat and Keppler, 2005).另外, 金红石还被利用其中的Nb和Cr含量判别含金红石榴辉岩的原岩性质(Zack et al., 2002; 李晓峰等, 2005; 王汝成等, 2005; 陈振宇等, 2006a; 余金杰等, 2006), 利用Zr含量温度计(Zack et al., 2004; Watson et al., 2006) 计算金红石榴辉岩形成温度(王汝成等, 2005; 高天山等, 2006; Spear et al., 2006; 余金杰等, 2006; 陈振宇等, 2006b; Miller et al., 2007). Rudnick et al. (2000)基于金红石的Nb、Ta含量和Nb/Ta比值推测含金红石榴辉岩可能是补偿硅酸岩地球大陆地壳与亏损地幔之间Nb/Ta比值不平衡的储库.Xiao et al. (2006)对大别-苏鲁超高压榴辉岩中金红石的研究认为, Nb和Ta的分馏最可能发生在没有金红石形成条件下的蓝片岩向角闪榴辉岩相变过程中, 热板片释放的低Nb/Ta比值流体运移至冷板片被含金红石榴辉岩所保存导致了热板片区域榴辉岩具有高Nb/Ta比值.这些研究大部分都是基于某些元素(如Nb、Ta) 在金红石中的高度相容性, 假设这些元素在金红石中的含量或比值可以代表全岩特征而进行的.因此, 一个很重要的问题是金红石中这些元素的含量或比值是否一定能够代表全岩特征?什么条件下金红石中的Nb/Ta比值和全岩一致或不一致?
大别-苏鲁超高压变质地体是华北与扬子板块在三叠纪俯冲-碰撞作用的产物, 被认为是目前世界上已知的最大超高压变质带.中国大陆科学钻探工程主孔(CCSD-MH) 位于江苏省东海县东南17 km的毛北村附近, 构造上位于苏鲁超高压变质带的南部.CCSD-MH榴辉岩不同程度地出现了金红石或其他含钛矿物(许志琴, 2004; 张泽明等, 2004), 为系统研究超高压变质岩中金红石的形成机制及其对微量元素分馏的控制作用提供了良好的天然实验室. 陈晶等(2005)发现与石榴石和绿辉石共生的金红石曾遭受过强烈的应力作用, 是在大陆板块俯冲或折返过程中形成的, 而石榴石和绿辉石中的微粒金红石包裹体则形成于大陆板块俯冲前或折返后的构造地质作用相对稳定期. 王汝成等(2005)、李晓峰等(2005)、余金杰等(2006)以及高天山等(2006)曾对大别-苏鲁榴辉岩中金红石的部分微量元素进行过测定, 通过金红石中Zr含量温度计得出的温度(608~787 ℃) 整体低于用其他方法所估算的榴辉岩峰期变质温度(Wang et al., 1993; Zhang et al., 1994, 1995; Zhang et al., 1996; Banno et al., 2000; Yao et al., 2000), 并推测导致偏低的原因可能是金红石中的Zr与锆石之间没有达到平衡, 或者是金红石中的Zr在退变质过程中由于重结晶作用和流体活动而发生了扩散, 或者这些金红石形成于超高压峰期变质作用之前.那么, 导致上述温度估计值之间差异的原因究竟是什么?已有的研究工作仅限于推测, 而没有给出很好的证据.
根据锆石U-Pb同位素定年所估计的大别-苏鲁超高压变质峰期年龄范围较宽(220~245 Ma; Hacker et al., 1998; Ayers et al., 2002; Wan et al., 2005; 张泽明等, 2005; Liu et al., 2006b; Zhao et al., 2007).利用榴辉岩中石榴石、绿辉石和金红石Sm-Nd等时线得出的超高压峰期变质年龄则较低(226.3±3.2 Ma; Li et al., 2000), 而金红石U-Pb定年获得的年龄值(218±1.2 Ma; Li et al., 2003) 与普遍认为的角闪岩相退变质年龄一致(< 220 Ma; Zhao et al., 2006; Liu et al., 2008).上述年龄差别暗示至少Sm、Nd、U、Pb在金红石形成后并没有封闭.那么金红石中的Zr是否也存在类似特征?因此, 榴辉岩中金红石的形成作用及是否存在后期改造对于利用金红石准确判别榴辉岩原岩的性质和解释金红石所记录的温度都具有重要意义.
本文在利用LA-ICP-MS对CCSD-MH榴辉岩中金红石微量元素含量和Nb/Ta比值进行系统原位微区分析的基础上, 对上述问题进行了深入探讨.研究结果表明榴辉岩全岩Nb/Ta比值主要受控于其原岩的形成过程; 金红石的Nb/Ta比值并不一定代表全岩特征, 而取决于全岩Nb、Ta和Ti的含量; 尽管金红石主体形成于超高压变质作用峰期, 但由于高温下的扩散作用以及折返过程中可能的流体活动使金红石Zr含量温度计并不能完全记录超高压变质作用的温度.
1. 研究样品与分析方法
在详细的岩相学工作基础上, 对所选择的10个榴辉岩样品进行了ICP-MS和XRF分析(表 1).这些样品分别采自CCSD-MH 100~1 100 m不同深度和类型的榴辉岩中.其中SD001、SD004、SD022和SD033属于低Mg-Ti榴辉岩(TiO2 < 2%), SD006、SD009、SD014、SD017和SD018属于高Ti榴辉岩(TiO2 > 2%), SD019属于高Mg榴辉岩.主要矿物组合为石榴石+绿辉石±石英±多硅白云母, 副矿物以Fe-Ti氧化物、金红石、黄铁矿和磷灰石为主, 有时出现退变质矿物绿帘石和角闪石.榴辉岩中金红石分2种类型: (1) 位于石榴石或绿辉石内部的包裹体金红石; (2) 位于石榴石或绿辉石颗粒之间的粒间金红石.其中包裹体金红石含量相对较少, 绝大多数金红石分布于石榴石或绿辉石颗粒之间.
表 1 CCSD-MH榴辉岩全岩主、微量元素组成Table Supplementary Table Major and trace element contents of the CCSD-MH eclogites样品粉碎工作和全岩常量元素组成分析在西北大学地质系大陆动力学国家重点实验室完成. 常量元素分析结果的精密度和准确度同Rudnick et al. (2004).全岩微量元素分析用中国地质大学地质过程与矿产资源国家重点实验室电感耦合等离子体质谱仪(Aglient 7500a ICP-MS) 完成.用于微量元素分析的样品制备流程为: (1) 准确称取50 mg左右约200目的岩石粉末样品放入Teflon坩埚中, 加入1.5 mL高纯HNO3摇匀后, 再加入1.5 mL高纯HF摇匀; (2) 把盛有样品和混合酸的Teflon坩埚放入特制钢套中, 在高温(190 ℃) 密闭烘箱中溶解样品48 h; (3) 把Teflon坩埚从钢套中取出后放在110 ℃电热板上加热3~4 h直至样品蒸干; (4) 在已蒸干的样品中加入1 mL高纯HNO3摇匀后再次蒸干; (5) 经过2次蒸干的样品冷却后加入2 mL高纯去离子水和1 mL高纯HNO3, 摇匀后重新放回钢套中, 并在高温(190 ℃) 密闭烘箱中加热24 h, 提取样品; (6) 把盛有样品提取液的Teflon坩埚从钢套中取出, 并用2%HNO3将提取液在聚乙烯塑料瓶中定溶至100 g左右, 待测(为了避免Nb、Ta在稀硝酸介质中不稳定问题, 微量元素ICP-MS分析工作均在样品定容后12 h内完成).
金红石原位微区微量元素分析由中国地质大学地质过程与矿产资源国家重点实验室的LA-ICP-MS完成, 激光剥蚀系统是配备有193 nm ArF-excimer激光器的GeoLas 2005.分析采用的激光剥蚀孔径为32 μm, 激光脉冲为8 Hz, 能量为80 mJ.电子探针(EMP) 分析在中国地质大学地质过程与矿产资源国家重点实验室完成, 分析仪器JEOL733, 分析条件为: 加速电压15 kV, 束流20 nA.对于单颗粒金红石的微量元素测定, 首先在电子显微镜下对加厚(~100 μm) 探针片进行观察, 将欲进行分析的单颗粒金红石进行标记, 然后对做好标记位置的金红石利用LA-ICP-MS进行微量元素含量分析.LA-ICP-MS分析数据处理采用NIST610作外标、Ti作内标元素, 利用Glitter4.4软件完成.其中NIST610微量元素含量使用的是GeoReM (http://georem.mpch-mainz.gwdg.de/) 在2006年给出的最新推荐值, 该推荐值纠正了以前NIST610推荐值中Mn、Ta含量偏低12%和20%的问题.为了对应进行EMP分析, 对经LA-ICP-MS剥蚀的金红石样品进行了轻微抛光以除去表面溅落的剥蚀颗粒.然后在激光剥蚀坑的附近利用EMP对金红石中Ti含量进行分析.
2. 结果
CCSD-MH榴辉岩主量元素MgO和TiO2含量差别很大, 且高场强元素Nb、Ta、Zr和Hf的含量以及Nb/Ta比值与TiO2含量具有负相关关系(图 1a、1c、1e、1f), 而Nb/Ta比值与Nb、Ta含量明显正相关(图 1b、1d) (相关系数R分别为0.89和0.77).高Mg榴辉岩全岩MgO和Al2O3的含量相对很高, 分别为12.42%和13.23%, 但是TiO2的含量却很低(0.39%).Ni、Cr含量也明显高于其他类型的榴辉岩, Nb、Ta、Zr和Hf的含量以及Nb/Ta比值略高于高Ti榴辉岩, 而与低Mg-Ti榴辉岩接近(图 1).高Ti榴辉岩具有异常高的TiO2 (2.75%~4.52%) 和TFe2O3 (15.34%~24.47%) 含量, 而MgO的含量(5.29%~7.37%) 明显低于高Mg榴辉岩(12.42%), 却比低Mg-Ti榴辉岩(1.96%~5.76%) 高.但是高Ti榴辉岩中高场强元素Nb、Ta、Zr和Hf含量以及Nb/Ta比值却明显较其他类型榴辉岩低很多(图 1).低Mg-Ti榴辉岩TiO2和MgO含量分别介于高Ti榴辉岩和高Mg榴辉岩之间, 但是却具有最高的Nb、Ta、Zr和Hf含量以及Nb/Ta比值(图 1).
金红石中微量元素含量在不同类型榴辉岩、同一类型榴辉岩不同金红石颗粒以及同一颗粒金红石的不同部位之间均有明显变化.V和Fe是金红石中含量最高的微量元素(表 2).其中Fe的含量在不同样品中相对较稳定, 而V的含量在不同全岩类型的金红石之间变化范围较大.过渡元素Sc、Cr、Cu、Zn、W、Mo以及碱土金属元素Mg和Al在金红石中含量相对较高(1.2×10-6~350×10-6), 而Rb、Ba、Cs和稀土元素(REE) 等含量大部分在LA-ICP-MS检出限以下.Sr、Pb在金红石中的含量相对较低, 但金红石边部较核部明显富集Pb、Sr (图 3a、3b).
表 2 榴辉岩中金红石微量元素的LA-ICP-MS测定结果Table Supplementary Table Trace element compositions of rutiles in eclogites analyzed by LA-ICP-MS金红石中高场强元素Nb、Ta、Zr、Hf含量非常富集.其中Nb和Ta含量主要受全岩Nb、Ta和TiO2含量控制, 与全岩Nb、Ta含量正相关(图 2a, 2c) (R=0.87和0.86), 与全岩TiO2含量明显负相关(R=-0.67和-0.70) (图 2b, 2d).Nb/Ta比值与全岩Nb、Ta含量没有明显相关性(图 2i), 但是与全岩TiO2含量负相关(R=-0.61) (图 2f).不同类型榴辉岩中金红石的Nb/Ta比值与其Nb、Ta含量明显正相关(R=0.96和0.94) (图 2g, 2h).同一样品中不同颗粒金红石之间的Nb/Ta比值与Ta含量明显负相关(R=-0.50~-0.96) (图 2g), 但是Nb/Ta比值与Nb的关系比较复杂, 低Mg-Ti榴辉岩金红石的Nb/Ta比值与其Nb含量正相关(R=0.54~0.70) (图 2h), 而高Mg和高Ti榴辉岩中金红石的Nb/Ta比值与其Nb的含量呈弱负相关(R=-0.34~-0.66) (图 2h).同一颗粒金红石Nb、Ta含量以及Nb/Ta比值比较均一, 核边之间没有明显变化(图 3d, 3f).但是Zr含量具有系统的核、边差异, 核部Zr含量高于边部(图 3c).
图 2 金红石Nb、Ta含量及Nb/Ta比值与全岩Nb、Ta和TiO2含量及Nb/Ta比值的关系(图例同图 1)相关系数(R) 由各样品平均值计算所获得.图 2a中相关系数计算未包括高Mg榴辉岩SD0019 (该样品金红石中Nb含量极高)Fig. 2. Plots of Nb and Ta contents and Nb/Ta ratios of rutiles versus Nb, Ta and TiO2contents and Nb/Ta ratios of parenteclogites利用Zack et al. (2004)的金红石中Zr含量温度计(ZMK) 计算的结果表明, 除高Mg榴辉岩SD019外, 其他榴辉岩样品中粒间金红石的温度范围为623~796 ℃ (平均645 ℃), 包裹体金红石的温度范围598~639 ℃ (平均614 ℃) (表 3, 图 4a).高Mg榴辉岩SD019中金红石计算的温度明显高于其他样品(达827 ℃; 表 3).利用Watson et al. (2006)的金红石中Zr含量温度计(WWT) 所得到的温度总体上较利用Zack et al. (2004)金红石温度计的计算结果低(平均低30~100 ℃) (表 3).但是利用2种Zr含量温度计获得的温度都表现出粒间金红石温度较包裹体金红石温度高的特征(表 3, 图 4a).另外, 粒间金红石核部温度明显高于边部温度(图 4b), 且对于直径小于300 μm的金红石核部温度与金红石直径具有正相关关系(图 4c).
表 3 金红石Zr含量地质温度计计算结果Table Supplementary Table Temperatures calculated by the Zr in rutile thermometer图 4 CCSD榴辉岩中金红石温度变化特征a、b.不同产状金红石Zr含量温度计计算结果.包裹体金红石温度低于粒间金红石温度, 而粒间金红石边部温度低于核部温度; c.利用Zack et al. (2004)温度计获得的粒间金红石核部温度对颗粒直径的变化; d.粒径200μm的金红石在不同温度条件下其中心位置Zr发生扩散所需要的时间及金红石中Zr封闭温度和粒径的关系.温度计算是根据Dodson (1973)的封闭温度公式和Cherniak et al. (2007)金红石中Zr扩散系数进行的, 分别计算了冷却速率为1℃Ma-1 (点线)、10℃Ma-1 (实线) 和100℃Ma-1 (虚线) 条件下的封闭温度Fig. 4. Temperature variations of rutile in eclogites from CCSD3. 讨论
3.1 金红石Nb/Ta比值能否反映全岩特征?
Nb、Ta具有相同的电价(+5) 和相近的原子半径(Shannon, 1976), 被认为在没有金红石存在条件下的分离结晶过程中具有一致的地球化学行为.实验研究表明金红石不仅会强烈富集Nb和Ta (DNb、Ta >> 1, 最高达540), 且会使Nb和Ta显著分异(DTa > DNb) (Green and Pearson, 1987; Brenan et al., 1994; Stalder et al., 1998; Foley et al., 2000, 2002; Green, 2000; Horng and Hess, 2000; Schmidt et al., 2004; Klemme et al., 2005; Xiong et al., 2005).因此, 金红石存在条件下的部分熔融或变质分异作用会显著影响所在体系的Nb/Ta比值.Nb和Ta在金红石中的高度相容性特征, 使全岩Nb和Ta绝大部分集中在金红石中(Foley et al., 2000; Schmidt et al., 2004; Klemme et al., 2005; Xiong et al., 2005).因此, Rudnick et al. (2000)认为金红石的Nb/Ta比值可以代表榴辉岩全岩的Nb/Ta比值, 并基于金红石Nb/Ta比值得出榴辉岩全岩Nb/Ta比值高于球粒陨石值.笔者对CCSD-MH榴辉岩全岩以及其中金红石单矿物精确的微量元素含量分析发现, 榴辉岩全岩Nb/Ta比值(17.61~3.36) 整体较原始地幔值(17.5±2.0) (Sun and McDonough, 1989) 低, 而金红石Nb/Ta比值(平均21.8) 则主要表现为高于原始地幔值的特征, Nb/Ta比值在榴辉岩全岩与金红石之间明显不一致(图 5).
图 5 榴辉岩和金红石Nb/Ta比值与Nb的变化关系大别榴辉岩全岩据Chavagnac and Jahn (1996)、Zhai and Cong (1996)、Liou et al. (1997)、Jahn (1998)、刘贻灿等(2000)、吴元保等(2000); 苏鲁榴辉岩全岩据Jahn et al. (1996)、Jahn (1998)、Zhang et al. (2000)、孙景贵和胡受奚(2000)、王大志等(2006).CCSD榴辉岩数据包括本文及张泽明等(2003, 2004)、刘勇胜等(2005)、戚学祥等(2006).原始地幔(PM) 值据Sun and McDonough (1989)Fig. 5. Nb/Ta ratios versus Nb contents for eclogites and rutiles通过改变榴辉岩中石榴石、绿辉石和金红石的比例对全岩TiO2、Nb、Ta含量和Nb/Ta比值进行平衡计算结果表明, 榴辉岩全岩Nb/Ta比值受全岩TiO2、Nb、Ta含量及矿物组合共同控制(图 6a, 6c).在全岩TiO2含量很低的情况下(< 0.5%), 金红石含量极少, 尽管其中Nb、Ta含量很高但对全岩贡献有限(图 6b, 6d).因此, 全岩Nb/Ta比值更多地受石榴石和绿辉石控制, 而表现出低Nb/Ta比值的特征(图 6a).而当全岩TiO2的含量为0.5%~2.5%时, 全岩Nb、Ta含量和Nb/Ta比值虽然受到矿物组合的影响(图 6), 但是主要受金红石含量控制, 而表现出金红石Nb/Ta比值的特征(图 6a).在全岩TiO2含量 > ~2.5%时, 金红石对全岩Nb、Ta的含量起主导作用(图 6b, 6d), Nb/Ta比值基本上完全受金红石Nb/Ta比值控制且相对稳定(图 6a, 6c).在TiO2含量高于0.5%而低于~3%时, 实测榴辉岩组成和模型计算结果吻合较好.但高Ti榴辉岩(TiO2 > 3%) 全岩Nb/Ta比值明显低于假设TiO2全部以金红石形式存在时的模型计算结果(图 6a, 6c), 表明高Ti榴辉岩中Ti主体不是以金红石形式存在的, 这与该类样品中大量出现的钛赤铁矿、钛铁矿一致(徐珏等, 2004).其极低的Nb/Ta比值应形成于原岩的岩浆作用过程.所以, 尽管金红石是榴辉岩中Nb、Ta的主要载体矿物, 但是并不能简单地用金红石Nb/Ta比值来代替全岩的Nb/Ta比值, 低TiO2 (< 0.5%) 或异常高Ti榴辉岩中金红石Nb/Ta比值不能直接反映全岩Nb/Ta比值特征(图 6a, 6c).因此, 在利用金红石Nb/Ta比值代替全岩比值时不仅要考虑全岩TiO2含量, 更重要的是考虑有多少全岩TiO2是由金红石贡献的.
图 6 对榴辉岩Nb/Ta比值与TiO2、Nb和Ta含量关系的模拟计算(图例同图 1)计算中假设榴辉岩由石榴石+绿辉石+金红石组成, 全岩TiO2含量由石榴石+绿辉石中溶解的TiO2和金红石中的TiO2组成.金红石中TiO2含量取为99%, 并以SD006中金红石、石榴石和绿辉石Nb、Ta和TiO2含量平均值代表金红石、石榴石和绿辉石中Nb、Ta和TiO2含量Fig. 6. Model calculations of Nb/Ta ratios versus TiO2, Nb and Ta contents for eclogites3.2 榴辉岩Nb/Ta比值的控制因素及其成因意义
大别-苏鲁不同类型榴辉岩全岩Nb/Ta比值变化范围很大(1.6~64), 但是总体表现出低于原始地幔值(17.5±2.0) (Sun and McDonough, 1989) 的特征(图 5), 不支持与俯冲有关的榴辉岩可能是高Nb/Ta比值的隐藏储库的观点(Rudnick et al., 2000).这与Kamber et al. (2002)对Franciscan Complex、Baja、Samana Peninsula以及Alps等地榴辉岩Nb/Ta比值的统计结果(平均值33.7) 不同.CCSD-MH榴辉岩明显偏低的Nb/Ta比值可能是其原岩特征的反映. Liu et al. (2008)对CCSD-MH榴辉岩的研究认为高Mg榴辉岩原岩是在玄武质岩浆房中形成的橄榄辉石岩堆晶产物, 高Ti榴辉岩原岩则是由铁质玄武岩浆中堆晶形成的富钛磁铁矿辉长岩, 低Mg-Ti榴辉岩原岩则可能是由高Mg和高Ti榴辉岩原岩高度分离结晶后的演化熔体.矿物与玄武质熔体之间Nb、Ta分配系数实验研究表明Nb、Ta在橄榄石(DNb/DTa=0.001 6~0.006 5/0.065~0.17; Dunn and Sen, 1994)、单斜辉石(DNb/DTa=0.008 1/0.010 2; Hauri et al., 1994)、石榴石(DNb/DTa=0.015~0.04/0.051~0.15; Green et al., 1989; Jenner et al., 1993)、斜长石(DNb/DTa=0.008~0.045/0.027~0.11; Dunn and Sen, 1994) 和磁铁矿(DNb/DTa=0.01~1.8/0.01~2.53; Nielsen et al., 1994; Nielsen and Beard, 2000) 中均是不相容元素, 且Nb比Ta更不相容, 即DNb < DTa < 1.在玄武质岩浆演化早期的橄榄石/辉石分离结晶作用中, Ti、Nb、Ta的分配主要发生在橄榄石+辉石与熔体之间, 因此高Mg榴辉岩原岩不仅具有低Ti、Nb、Ta含量, 同时其Nb/Ta比值也低于原始熔体(图 1).伴随着玄武质岩浆演化, 熔体中FeO和TiO2含量增加, 导致钛磁铁矿分离结晶从而形成具有非耦合的高Ti、低Nb (Ta) 特征的富钛磁铁矿辉长质堆晶产物(刘勇胜等, 2005; Liu et al., 2008).但是Nb、Ta在磁铁矿中同样具有DNb < DTa < 1的特征(Nielsen et al., 1994; Nielsen and Beard, 2000).因此, 高Ti榴辉岩原岩(富钛磁铁矿辉长质堆晶岩) 具有比原始玄武质熔体更低的Nb/Ta比值和Nb、Ta含量(图 1).橄榄石、单斜辉石、斜长石和磁铁矿的分离结晶作用导致残余熔体中Nb、Ta含量和Nb/Ta比值随分离结晶作用程度的增加而增加.因此经历过橄榄石、单斜辉石、斜长石和磁铁矿分离结晶的演化熔体固结或高度分离结晶形成的低Mg-Ti榴辉岩的原岩Nb、Ta含量和Nb/Ta比值将会明显高于高Mg和高Ti榴辉岩原岩(图 1).CCSD-MH中榴辉岩所特有的非耦合的高Ti、低Nb (Ta) 特征(图 1a, 1c), 以及Nb/Ta比值与Nb、Ta含量正相关(图 1b, 1d) 而与TiO2含量负相关的特征(图 1f) 均表明全岩高度变化的Nb/Ta比值特征是其原岩特征的反映, 记录了原岩形成过程中的岩浆分离结晶作用.
金红石Nb、Ta含量与榴辉岩全岩Nb、Ta含量正相关(图 2a、2c), 与全岩TiO2含量负相关(图 2b、2d), 而且Nb/Ta比值明显高于全岩(图 5), 且与全岩Nb/Ta比值无明显相关性(图 2e), 表明金红石应形成于超高压变质作用, 而不是继承自原岩的结晶作用.继承自原岩和熔体平衡的金红石其Nb/Ta比值及Nb、Ta含量应主要受分离结晶程度和结晶产物组成控制.因为Nb和Ta在金红石中高度相容(DNb、Ta >> 1, 最高达540), 且Ta比Nb更相容(DTa > DNb) (Green and Pearson, 1987; Brenan et al., 1994; Stalder et al., 1998; Foley et al., 2000, 2002; Green, 2000; Horng and Hess, 2000; Schmidt et al., 2004; Klemme et al., 2005; Xiong et al., 2005), 分离结晶形成的金红石Nb/Ta比值与全岩Nb、Ta含量应具有明显的负相关性(图 7a), 且金红石Nb/Ta比值与全岩Nb/Ta比值基本一致(图 7b).因此, CCSD-MH榴辉岩中的金红石应是在超高压或退变质作用过程中相变形成的.
图 7 含金红石辉长岩分离结晶作用过程中Nb/Ta比值的变化计算中以大陆溢流玄武岩平均值(Nb=15.1 ppm, Ta=1.07 ppm) 作为原始熔体, 假设结晶辉长岩的矿物组成为60%~55%斜长石+25%单斜辉石+10%斜方辉石+4.5%橄榄石+0%~5%金红石.根据Nb、Ta在斜长石(Dunn and Sen, 1994)、单斜辉石(Hauri et al., 1994)、斜方辉石(Dunn and Sen, 1994)、橄榄石(Ulmer, 1989) 以及金红石(Jenner et al., 1993) 与熔体之间的分配系数, 分别计算了分离结晶程度为30%和70%时全岩和金红石的Nb、Ta含量和Nb/Ta比值.图中虚线代表不同金红石含量(增量为0.5%)Fig. 7. Nb/Ta variations during fractionation crystallization of rutile-bearing gabbro对于高Ti榴辉岩, 原岩变质作用前Ti主要赋存于钛磁铁矿中, 而Nb和Ta则在钛磁铁矿、单斜辉石、斜长石中都有分布.以钛磁铁矿形式出现的Ti相对Nb (Ta) 过剩.超高压变质作用导致钛磁铁矿分解形成“大量”金红石, “少量”的Nb和Ta从钛磁铁矿、单斜辉石、斜长石中析出而在“大量”金红石中重新分配, 因此金红石Nb、Ta含量与全岩TiO2含量负相关(图 2b、2d).高Mg榴辉岩原岩Ti、Nb、Ta均赋存在橄榄石、单斜辉石、斜长石中, 超高压变质形成金红石的数量则取决于在特定温度-压力条件下Ti在石榴石、绿辉石中的溶解度(溶解度越大, 形成的金红石数量则越少), 而Nb、Ta由于其在金红石中的高度相容性则会从石榴石、绿辉石中析出而进入少量的金红石中, 因此高Mg榴辉岩中金红石高度富集Nb、Ta (图 2a、2d).低Mg-Ti榴辉岩原岩是由高Mg和高Ti榴辉岩原岩分馏结晶后的残余熔体形成, 其中Ti、Nb、Ta等的分配整体符合稀溶液定律.因此, 尽管其原岩在超高压变质作用中形成金红石的数量同样与Ti在石榴石、绿辉石中的溶解度有关, 但由于Ti含量较高, 溶解于石榴石、绿辉石中的Ti远远低于全岩TiO2总含量.因此低Mg-Ti榴辉岩形成金红石时, Nb、Ta既不会相对Ti过饱和(如高Mg榴辉岩) 也不会相对Ti严重不足(如高Ti榴辉岩), 而与全岩的特征保持一致(图 2a、2d).因而, 低Mg-Ti榴辉岩实测全岩Nb/Ta比值与理论计算值相对一致, 而不同于高Mg和高Ti榴辉岩(图 6a、6c).
3.3 金红石中Zr含量温度计与Zr的扩散作用
金红石-锆石-石英体系中金红石Zr含量地质温度计近年来得到了较多应用(Zack et al., 2004; 王汝成等, 2005; Spear et al., 2006; Watson et al., 2006).相对于榴辉岩峰期变质温度(700~860 ℃) (Wang et al., 1993; Zhang et al., 1994, 1995; Zhang et al., 1996; Banno et al., 2000; Yao et al., 2000), 前人利用金红石中Zr含量地质温度计(Zack et al., 2004; Watson et al., 2006) 获得的大别-苏鲁榴辉岩中金红石的形成温度为608~787 ℃ (王汝成等, 2005; 陈振宇等, 2006a, 2006b; 高天山等, 2006), 明显低于通过其他温度计获得的超高压峰期变质温度.笔者对CCSD-MH金红石形成温度的计算结果表明, 除高Mg榴辉岩相对较高以外(最高827 ℃), 其他类型榴辉岩的变化范围在598 ℃至787 ℃之间, 同样呈现偏低的特征(图 4).金红石Zr含量计算的温度能否记录其形成温度, 要视其形成后Zr的封闭情况(影响因素包括温度、流体、扩散速率等) 而定(Zack et al., 2004; Watson et al., 2006).利用Zr含量地质温度计获得的温度偏低的可能原因包括: (1) 榴辉岩原岩中Zr-Ti“不平衡”, 即Zr相对Ti不足; (2) 金红石不是在超高压变质峰期形成; (3) 金红石形成后受后期作用影响导致Zr含量降低, 从而使得计算温度偏低.
CCSD-MH榴辉岩全岩呈现明显的非耦合高TiO2、低Zr特征(图 8), 对于TiO2含量异常高(> ~2.5%) 的榴辉岩, 在超高压变质作用形成大量金红石的过程中可能存在Zr-Ti“不平衡”, 即Zr含量相对其他类型榴辉岩可能是不饱和的(至少是局部不饱和).Zack et al. (2004)的经验公式以及Watson et al. (2006)给出的实验公式均未考虑这一现象.因此, 这可能是利用高Ti榴辉岩中金红石计算的温度偏低的原因之一.而高Mg榴辉岩在超高压变质作用中形成的金红石含量极少, Zr含量相对饱和, 以致有足够的Zr进入金红石.因此, 高Mg榴辉岩中的金红石记录的较高的温度可能更接近榴辉岩相变质温度.
图 8 榴辉岩及大洋中脊玄武岩中Zr对TiO2变化关系大别和苏鲁榴辉岩数据来源同图 5.阴影区域为大洋中脊玄武岩(MORB), 据http://petdb.columbia.petdb/Fig. 8. Plot of Zr versus TiO2for eclogites and the middle ocean ridge basalts (MORB)另外, 如果金红石并非形成于峰期变质阶段, 而是形成于进变质阶段或者折返过程中的退变质冷却阶段, 同时金红石中Zr又能够在形成后及时封闭, 那么金红石Zr含量所记录的温度同样会低于榴辉岩峰期变质温度.CCSD-MH榴辉岩中金红石多数以粒间金红石的形式存在, 包裹体金红石数量相对较少(王硕等, 2006).整体而言, 包裹体金红石所记录的温度(平均614 ℃) 低于粒间金红石(平均640 ℃) (图 4a).包裹体金红石由于受寄主矿物的保护, 很少遭受后期改造作用(Zack et al., 2004), 记录的温度应该是其形成温度的真实反映, 但其记录的温度与大别-苏鲁榴辉岩经历绿帘角闪岩相进变质阶段的温度一致(500~600 ℃; 张泽明等, 2005), 暗示包裹体金红石可能主要形成于进变质作用阶段, 记录了进变质阶段的温度变化.
尽管粒间金红石记录的温度高于包裹体金红石(图 4a), 但仍低于榴辉岩的峰期变质温度(700~860 ℃; Wang et al., 1993; Zhang et al., 1994, 1995; Zhang et al., 1996; Banno et al., 2000; Yao et al., 2000).这暗示金红石形成后其中的Zr并没有封闭, 或者受到了后期作用的改造.粒间金红石边部温度普遍低于核部温度(图 4b), 且核部温度与其粒径呈正相关关系(图 4c).这种核-边温度变化及其与粒径的关系充分反映了后期Zr扩散作用的影响.另一方面, 这些金红石边部明显的高Pb、Sr特征(图 3a、3b) 与退变质过程中流体活动的特征一致(张斌辉等, 2006).利用Dodson (1973)元素封闭温度公式以及Cherniak et al. (2007)给出的金红石Zr扩散系数的计算结果表明: 对于半径100~200 μm的金红石, 在冷却速率为1~100 ℃/Ma的条件下Zr的封闭温度范围为560~800 ℃ (图 4d).因此, 金红石中Zr “较容易”受到扩散作用改造.在900 ℃条件下, 半径为100 μm的金红石其中心位置Zr发生扩散所需要的时间仅为0.04 Ma (图 4d), 即便是在600 ℃条件下其中心位置Zr发生扩散所需要的时间也仅16.6 Ma (图 4d), 该时间接近或低于榴辉岩变质作用所持续的时间(10~20 Ma) (Zheng et al., 2003; 张泽明等, 2005; Liu et al., 2006b; 刘福来等, 2006a; 刘贻灿等, 2006; Wu et al., 2006).如果存在流体活动, 金红石中Zr发生扩散所需要的时间将更短(Watson et al., 2006).因此, 扩散作用和/或后期流体作用的影响可能是导致苏鲁-大别榴辉岩金红石Zr含量温度计计算结果偏低的重要原因.
4. 结论
对CCSD-MH 100~1 100 m榴辉岩全岩主、微量元素以及金红石的LA-ICP-MS分析研究表明, 榴辉岩全岩低Nb/Ta比值的特征反映了其原岩形成过程中与玄武质岩浆平衡分配的结果.金红石高Nb/Ta比值的特征则是Nb、Ta在超高压变质作用过程中在不同矿物之间重新分配形成.Nb/Ta比值在全岩和金红石之间的非耦合关系表明尽管金红石是榴辉岩中Nb、Ta的主要寄主矿物, 但并不一定能够完全表征榴辉岩全岩特征, 而与榴辉岩原岩中Ti的赋存形式及含量有关.
本研究表明金红石Zr含量地质温度计在实际应用中可能存在一定的局限性, 对于异常高Ti榴辉岩可能是不合适的.粒间金红石中与流体活动密切相关的元素Pb、Sr所具有的边部富集特征以及金红石核部温度与粒径之间的正相关性充分反映粒间金红石形成后其中的Zr并没有封闭, 这应是导致大别-苏鲁榴辉岩中金红石Zr含量所记录温度偏低的重要原因, 为金红石中Zr含量受到扩散作用和/或后期流体活动的影响提供了证据.
致谢: 感谢刘慧芳老师在EMP分析中所给予的帮助和指导. -
图 2 金红石Nb、Ta含量及Nb/Ta比值与全岩Nb、Ta和TiO2含量及Nb/Ta比值的关系(图例同图 1)
相关系数(R) 由各样品平均值计算所获得.图 2a中相关系数计算未包括高Mg榴辉岩SD0019 (该样品金红石中Nb含量极高)
Fig. 2. Plots of Nb and Ta contents and Nb/Ta ratios of rutiles versus Nb, Ta and TiO2contents and Nb/Ta ratios of parenteclogites
图 4 CCSD榴辉岩中金红石温度变化特征
a、b.不同产状金红石Zr含量温度计计算结果.包裹体金红石温度低于粒间金红石温度, 而粒间金红石边部温度低于核部温度; c.利用Zack et al. (2004)温度计获得的粒间金红石核部温度对颗粒直径的变化; d.粒径200μm的金红石在不同温度条件下其中心位置Zr发生扩散所需要的时间及金红石中Zr封闭温度和粒径的关系.温度计算是根据Dodson (1973)的封闭温度公式和Cherniak et al. (2007)金红石中Zr扩散系数进行的, 分别计算了冷却速率为1℃Ma-1 (点线)、10℃Ma-1 (实线) 和100℃Ma-1 (虚线) 条件下的封闭温度
Fig. 4. Temperature variations of rutile in eclogites from CCSD
图 5 榴辉岩和金红石Nb/Ta比值与Nb的变化关系
大别榴辉岩全岩据Chavagnac and Jahn (1996)、Zhai and Cong (1996)、Liou et al. (1997)、Jahn (1998)、刘贻灿等(2000)、吴元保等(2000); 苏鲁榴辉岩全岩据Jahn et al. (1996)、Jahn (1998)、Zhang et al. (2000)、孙景贵和胡受奚(2000)、王大志等(2006).CCSD榴辉岩数据包括本文及张泽明等(2003, 2004)、刘勇胜等(2005)、戚学祥等(2006).原始地幔(PM) 值据Sun and McDonough (1989)
Fig. 5. Nb/Ta ratios versus Nb contents for eclogites and rutiles
图 6 对榴辉岩Nb/Ta比值与TiO2、Nb和Ta含量关系的模拟计算(图例同图 1)
计算中假设榴辉岩由石榴石+绿辉石+金红石组成, 全岩TiO2含量由石榴石+绿辉石中溶解的TiO2和金红石中的TiO2组成.金红石中TiO2含量取为99%, 并以SD006中金红石、石榴石和绿辉石Nb、Ta和TiO2含量平均值代表金红石、石榴石和绿辉石中Nb、Ta和TiO2含量
Fig. 6. Model calculations of Nb/Ta ratios versus TiO2, Nb and Ta contents for eclogites
图 7 含金红石辉长岩分离结晶作用过程中Nb/Ta比值的变化
计算中以大陆溢流玄武岩平均值(Nb=15.1 ppm, Ta=1.07 ppm) 作为原始熔体, 假设结晶辉长岩的矿物组成为60%~55%斜长石+25%单斜辉石+10%斜方辉石+4.5%橄榄石+0%~5%金红石.根据Nb、Ta在斜长石(Dunn and Sen, 1994)、单斜辉石(Hauri et al., 1994)、斜方辉石(Dunn and Sen, 1994)、橄榄石(Ulmer, 1989) 以及金红石(Jenner et al., 1993) 与熔体之间的分配系数, 分别计算了分离结晶程度为30%和70%时全岩和金红石的Nb、Ta含量和Nb/Ta比值.图中虚线代表不同金红石含量(增量为0.5%)
Fig. 7. Nb/Ta variations during fractionation crystallization of rutile-bearing gabbro
图 8 榴辉岩及大洋中脊玄武岩中Zr对TiO2变化关系
大别和苏鲁榴辉岩数据来源同图 5.阴影区域为大洋中脊玄武岩(MORB), 据http://petdb.columbia.petdb/
Fig. 8. Plot of Zr versus TiO2for eclogites and the middle ocean ridge basalts (MORB)
表 1 CCSD-MH榴辉岩全岩主、微量元素组成
Table 1. Major and trace element contents of the CCSD-MH eclogites
表 2 榴辉岩中金红石微量元素的LA-ICP-MS测定结果
Table 2. Trace element compositions of rutiles in eclogites analyzed by LA-ICP-MS
表 3 金红石Zr含量地质温度计计算结果
Table 3. Temperatures calculated by the Zr in rutile thermometer
-
[1] Audetat, A., Keppler, H., 2005. Solubility of rutile in subduction zone fluids, as determined by experiments in the hydrothermal diamond anvil cell. Earth and Planetary Science Letters, 232 (3-4): 393-402. doi: 10.1016/j.epsl.2005.01.028 [2] Ayers, J. C., Dunkle, S., Gao, S., et al., 2002. Constraints on timing of peak and retrograde metamorphismin the Dabie Shan ultrahigh-pressure metamorphic belt, east-central China, using U-Th-Pb dating of zircon and monazite. Chemical Geology, 186 (3-4): 315-331. doi: 10.1016/S0009-2541(02)00008-6 [3] Banno, S., Enami, M., Hirajima, T., et al., 2000. Decompression P-T path of coesite eclogite to granulite from Weihai, eastern China. Lithos, 52 (1-4): 97-108. doi: 10.1016/S0024-4937(99)00086-9 [4] Brenan, J. M., Shaw, H. F., Phinney, D. L., et al., 1994. Rutile-aqueous fluid partitioning of Nb, Ta, Hf, Zr, U and Th: Implications for high field strength element depletions inisland-arc basalts. Earth and Planetary Science Letters, 128 (3-4): 327-339. doi: 10.1016/0012-821X(94)90154-6 [5] Brophy, J. G., Marsh, B. D., 1986. On the origin of high-alumina arc basalt and the mechanics of melt extraction. J. Petrol. , 27: 763-789. doi: 10.1093/petrology/27.4.763 [6] Chavagnac, V., Jahn, B. M., 1996. Coesite-bearing eclogites from the Bixiling complex, Dabie Mountains, China: SmNd ages, geochemical characteristics and tectonic implications. Chemical Geology, 133 (1-4): 29-51. doi: 10.1016/S0009-2541(96)00068-X [7] Chen, J., Xu, Z. Q., Li, X. P., 2005. The formation of nanometer twins of rutile and its textural characteristics in UHP eclogite. Acta Petrologica Sinica, 21 (2): 399-404 (in Chinese with English abstract). [8] Chen, Z. Y., Wang, D. H., Chen, Y. C., et al., 2006a. Mineral geochemistry of rutilein eclogite anditsimplications. Earth Science—Journal of China University of Geosciences, 31 (4): 533-538, 550 (in Chinese with English abstract). [9] Chen, Z. Y., Yu, J. J., Xu, J., et al., 2006b. EMP trace elements analysis of rutile in eclogite. J. Chin. Electr. Microsc. Soc. , 25 (Suppl.): 295-296 (in Chinese). [10] Cherniak, D. J., Manchester, J., Watson, E. B., 2007. Zr and Hf diffusionin rutile. Earthand Planetary Science Letters, 261 (1-2): 267-279. doi: 10.1016/j.epsl.2007.06.027 [11] Dodson, M. H., 1973. Closure temperature in cooling geochronological and petrological systems. Contrib. Mineral. Petrol. , 40: 259-274. doi: 10.1007/BF00373790 [12] Dunn, T., Sen, C., 1994. Mineral/matrix partition coefficients for orthopyroxene, plagioclase, and olivine in basaltic to andesitic systems: A combined analytical and experimental study. Geochimica et Cosmochimica Acta, 58 (2): 717-733. doi: 10.1016/0016-7037(94)90501-0 [13] Foley, S. F., Barth, M. G., Jenner, G. A., 2000. Rutile/melt partition coefficients for trace elements and an assessment of theinfluence of rutile onthe trace element characteristics of subduction zone magmas. Geochimica et Cosmochimica Acta, 64 (5): 933-938. doi: 10.1016/S0016-7037(99)00355-5 [14] Foley, S. F., Massimo, T., Riccardo, V., 2002. Growth of early continental crust controlled by melting of amphibolite in subduction zones. Nature, 417: 837-840. doi: 10.1038/nature00799 [15] Gao, T. S., Zheng, Y. F., Chen, R. X., 2006. Applications of the different types of mineral-pair geothermometers to eclogite at Huangzhenin the Dabie orogen. Acta Petrologica Sinica, 22 (7): 1957-1968 (in Chinese withEnglish abstract). [16] Green, T. H., 2000. New partition coefficient determinations pertinent to hydrous melting processes in subduction zones. State of the Arc 2000: Processes and timescales. Wellington, 92-95. [17] Green, T. H., Pearson, N. J., 1987. An experimental study of Nb and Ta partitioning between Ti-rich minerals and silicate liquids at high pressure and temperature. Geochimica et Cosmochimica Acta, 51 (1): 55-62. doi: 10.1016/0016-7037(87)90006-8 [18] Green, T. H., Sie, S. H., Ryan, C. G., et al., 1989. Proton microprobe-determined partitioning of Nb, Ta, Zr, Sr and Y between garnet, clinopyroxene and basaltic magma at high pressure and temperature. Chemical Geology, 74 (3-4): 201-216. doi: 10.1016/0009-2541(89)90032-6 [19] Hacker, B. R., Ratschbacher, L., Webb, L., et al., 1998. U/Pb zircon ages constrain the architecture of the ultrahigh-pressure Qinling-Dabie orogen, China. Earth and Planetary Science Letters, 161: 215-230. doi: 10.1016/S0012-821X(98)00152-6 [20] Hauri, E. H., Wagner, T. P., Grove, T. L., 1994. Experimental and natural partitioning of Th, U, Pb and other trace elements between garnet, clinopyroxene and basaltic melts. Chemical Geology, 117: 149-166. doi: 10.1016/0009-2541(94)90126-0 [21] Horng, W. S., Hess, P. C., 2000. Partition coefficients of Nb and Ta between rutile and anhydrous haplogranite melts. Contrib. Mineral. Petrol. , 138: 176-185. doi: 10.1007/s004100050016 [22] Jahn, B. M., 1998. Geochemical andisotopic characteristics of UHP eclogites and ultramafic rocks of the Dabie orogen: Implications for continental subduction and collisional tectonics. In: Hacker, B. R., Liou, J. G., eds., When continents collide: Geodynamics and geochemistry of ultrahigh-pressure rocks. Kluwer Academic Publishers, Netherlands, 203-239. [23] Jahn, B. M., Cornichet, J., Cong, B. L., et al., 1996. Ultrahigh-εNd eclogites froman ultrahigh-pressure metamorphic terrane of China. Chemical Geology, 127 (1-3): 61-79. doi: 10.1016/0009-2541(95)00108-5 [24] Jenner, G. A., Foley, S. F., Jackson, S. E., et al., 1993. Determination of partition coefficients for trace elements in high pressure-temperature experimental run products by laser ablation microprobe-inductively coupled plasma-mass spectrometry (LAM-ICP-MS). Geochimica et Cosmochimica Acta, 57 (23-24): 5099-5103. doi: 10.1016/0016-7037(93)90611-Y [25] Kamber, B. S., Ewart, A., Collerson, K. D., et al., 2002. Fluid-mobile trace element constraints on the role of slab melting and implications for Archaean crustal growth models. Contrib. Mineral. Petrol. , 144: 38-56. doi: 10.1007/s00410-002-0374-5 [26] Klemme, S., Prowatke, S., Hametner, K., et al., 2005. Partitioning of trace elements between rutile and silicate melts: Implications for subduction zones. Geochimica et Cosmochimica Acta, 69 (9): 2361-2371. doi: 10.1016/j.gca.2004.11.015 [27] Li, Q., Li, S. Q., Zheng, Y. F., et al., 2003. Ahigh precision U-Pb age of metamorphic rutile in coesite-bearing eclogite from the Dabie Mountains in Central China: A newconstraint on the cooling history. Chemical Geology, 200 (3-4): 255-265. doi: 10.1016/S0009-2541(03)00194-3 [28] Li, S. Q., Jagoutz, E., Chen, Y. I., et al., 2000. Sm-Nd and Rb-Sr isotopic chronology and cooling history of ultrahigh pressure metamorphic rocks and their country rocks at Shuanghe in the Dabie Mountains, Central China. Geochimicaet Cosmochimica Acta, 64 (6): 1077-1093. doi: 10.1016/S0016-7037(99)00319-1 [29] Li, X. F., Chen, Z. Y., Wang, R. C., et al., 2005. Mineralogical and geochemical features of hair-like rutilein rocked quartzin Donghai, Jiangsu Province, China. Acta Petrologica Sinica, 21 (2): 475-481 (in Chinese with English abstract). [30] Liou, J. G., Zhang, R. Y., Jahn, B. M., 1997. Petrology, geochemistry and isotope data on an ultrahigh-pressure jadeite quartzite from Shuanghe, Dabie Mountains, eastcentral China. Lithos, 41 (1-3): 59-78. doi: 10.1016/S0024-4937(97)82005-1 [31] Liu, D. Y., Jian, P., Kr ner, A., et al., 2006. Dating of prograde metamorphic events deciphered frome pisodic zircon growth in rocks of the Dabie-Sulu UHP complex, China. Earth and Planetary Science Letters, 250 (3-4): 650-666. doi: 10.1016/j.epsl.2006.07.043 [32] Liu, F. L., Xue, H. M., Xu, Z. Q., et al., 2006a. SHRI MP U-Pb zircon dating fromeclogite lens in marble, Shuanghe area, Dabie UHPterrane: Restriction on the prograde, UHP and retrograde metamorphic ages. Acta Petrologica Sinica, 22 (7): 1761-1778 (in Chinese with English abstract). [33] Liu, F. L., Gerdes, A., Liou, J. G., et al., 2006b. SHRI MP U/Pb zircon dating from Sulu-Dabie dolomitic marble, eastern China: Constraints on prograde, ultrahigh pressure and retrograde metamporphic ages. J. Metamorph. Geol. , 24 (7): 569-589. doi: 10.1111/j.1525-1314.2006.00655.x [34] Liu, F. L., Gerdes, A., Zeng, L. S., et al., 2008. SHRI MP U-Pb dating, trace elements and the Lu-Hf isotope system of coesite-bearing zircon from amphibolite in the SW Sulu UHPterrane, eastern China. Geochim. Cosmochim. Acta, doi: 10.1016/j.gca.2008.04.007. [35] Liu, Y. C., Li, S. G., Xu, S. T., et al., 2006. Ultrahigh-pressure metamorphic evidence for gneiss fromthe northern Dabie complex, Central China: Message fromzircon. Acta Petrologica Sinica, 22 (7): 1827-1832 (in Chinese with English abstract). [36] Liu, Y. C., Xu, S. T., Li, S. G., et al., 2000. Eclogites from the northern Dabie Mountains, eastern China: Geochemical characteristics, Sr-Nd isotopic compositions and tectonic implications. Sciencein China (Series D), 43: 178-188. doi: 10.1007/BF02911943 [37] Liu, Y. S., Iong, K. Q., Kelemen, P. B., et al., 2008. Geochemistry and magmatic history of eclogites and ultramafic rocks from the Chinese Continental Scientific Drill. Chem. Geol. , 247: 133-153. doi: 10.1016/j.chemgeo.2007.10.016 [38] Liu, Y. S., Zhang, Z. M., Lee, C. T., et al., 2005. Decoupled high-Ti fromlow-Nb (Zr) of eclogites fromthe CCSD: Implications for magnetite fractional crystallization in basalt chamber. Acta Petrologica Sinica, 21 (2): 339-346 (in Chiese with English abstract). [39] Miller, C., Zanetti, A., Thoni, M., et al., 2007. Eclogitisation of gabbroic rocks: Redistribution of trace elements and Zr in rutile thermometryin an Eo-Alpine subduction zone (eastern Alps). Chemical Geology, 239 (1-2): 96-123. doi: 10.1016/j.chemgeo.2007.01.001 [40] Morris, J. D., Leeman, W. P., Tera, F., 1990. The subducted component inisland arc lavas: Constraints from Be isotopes and B-Be systematics. Nature, 344: 31-36. doi: 10.1038/344031a0 [41] Nielsen, R. L., Beard, J. S., 2000. Magnetite-melt HFSE partitioning. Chemical Geology, 164 (1-2): 21-34. doi: 10.1016/S0009-2541(99)00139-4 [42] Nielsen, R. L., Forsythe, L. M., Gallahan, W. E., et al., 1994. Major and trace-element magnetite-melt equilibria. Chemical Geology, 117 (1-4): 167-191. doi: 10.1016/0009-2541(94)90127-9 [43] Qi, X. X., Tang, Z. M., Yan, L., 2006. Geochemical characteristics of garnet-amphibolite fromthe pre-pilot hole of the Chinese Continental Scientific Drilling (CCSD-PP2) and its geological significance. Earth Science—Journal of China University of Geosciences, 31 (4): 539-550 (in Chinese with English abstract). [44] Ringwood, A. E., 1990. Slab-mantle interactions 3. Petrogenesis of intraplate magmas and structure of the upper mantle. Chem. Geol. , 82: 187-207. doi: 10.1016/0009-2541(90)90081-H [45] Rudnick, R. L., Barth, M., Horn, I., et al., 2000. Rutilebearing refractory eclogites: Missinglink between continents and depleted mantle. Science, 287: 278-281. doi: 10.1126/science.287.5451.278 [46] Rudnick, R. L., Gao, S., Ling, W. L., et al., 2004. Petrology and geochemistry of spinel peridotite xenoliths from Hannuoba and Qixia, North China craton. Lithos, 77 (1-4): 609-637. doi: 10.1016/j.lithos.2004.03.033 [47] Schmidt, M. W., Dardon, A., Chazot, G., et al., 2004. The dependence of Nb and Ta rutile-melt partitioning on melt composition and Nb/Ta fractionation during subduction processes. Earth and Planetary Science Letters, 226 (3-4): 415-432. doi: 10.1016/j.epsl.2004.08.010 [48] Shannon, R. D., 1976. Revised effective ionicradii and systematic studies of interatomic distances in chalides and halcogenides. Acta Crystallographica Section A, 32 (5): 751-767. doi: 10.1107/S0567739476001551 [49] Spear, F. S., Wark, D. A., Cheney, J. T., et al., 2006. Zr-inrutile thermometry in blueschists from Sifnos, Greece. Contrib. Mineral. Petrol. , 152: 375-385. doi: 10.1007/s00410-006-0113-4 [50] Stalder, R., Foley, S. F., Brey, G. P., et al., 1998. Mineralaqueous fluid partitioning of trace elements at 900-1200℃ and 3.0-5.7GPa: Newexperimental data for garnet, clinopyroxene, and rutile, and implications for mantle metasomatism. Geochimicaet Cosmochimica Acta, 62 (10): 1781-1801. doi: 10.1016/S0016-7037(98)00101-X [51] Sun, J. G., Hu, S. X., 2000. Research on geochemical variation during the retrograde metamorphic processes of eclogite—From the eclogite bodiesinthe Rongcheng-Weihai rock belt of Sulu (Jiangsu-Shandong) UHMB. J. Mineral Petrol. , 20 (1): 17-23 (in Chinese with English abstract). [52] Sun, S. S., McDonough, W. F., 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. In: Saunders, A. D., Norry, M. J., eds., Magmatismin the ocean basins. Geological Society Special Pubications, London, 313-345. [53] Tatsumi, Y., Nakamura, N., 1986. Composition of aqueous fluid from serpentinite in the subducted lithosphere. Geochem. J. , 20: 191-196. doi: 10.2343/geochemj.20.191 [54] Ul mer, P., 1989. Partitioning of high field strength elements among olivine, pyroxenes, garnet and calc alkaline picrobasalt: Experimental results and an application. International Journal of Mass Spectrometry and Ion Physics, : 42-47. [55] Ul mer, P., Trommsdorff, V., 1995. Serpentine stability to mantle depths and subduction-related magmatism. Science, 268: 858-861. doi: 10.1126/science.268.5212.858 [56] Wan, Y., Li, R. W., Wilde, S. A., et al., 2005. UHP metamorphismand exhumation of the Dabie orogen, China: Evidence from SHRI MP dating of zircon and monazite froma UHP granitic gneiss cobble fromthe Hefei basin. Geochim. Cosmochim. Acta, 69 (17): 4333-4348. doi: 10.1016/j.gca.2005.03.055 [57] Wang, D. Z., Zhang, Z. M., Shen, K., et al., 2006. Origin of Ti-rich eclogite and rutile deposit in the Donghai area of the southern Sulu UHP metamorphic belt, China. Geological Bulletin of China, 25 (7): 839-849 (in Chinesewith English abstract). [58] Wang, Q. C., Ishi watari, A., Zhao, Z. Y., et al., 1993. Coesitebearing granulite retrograded fromeclogite in Weihai, eastern China. Eur. J. Mineral. , 5 (1): 141-152. doi: 10.1127/ejm/5/1/0141 [59] Wang, R. C., Wang, S., Qiu, J. S., et al., 2005. Rutile in the UHP eclogites fromthe CCSD main drillhole (Donghai, eastern China): Trace-element geochemistry and metallogenetic implications. Acta Petrologica Sinica, 21 (2): 465-474 (in Chinese with English abstract). [60] Wang, S., Wang, R. C., Qiu, J. S., et al., 2006. Mineral inclusions in rutile of Sulu UHP eclogites fromthe main drill hole of CCSD, Donghai. Acta Petrologica et Mineralogica, 25 (1): 61-70 (in Chinese with English abstract). [61] Watson, E. B., Wark, D. A., Thomas, J. B., 2006. Crystallization thermometers for zircon and rutile. Contrib. Mineral. Petrol. , 151 (4): 413-433. doi: 10.1007/s00410-006-0068-5 [62] Wu, Y. B., Chen, D. G., Cheng, H., et al., 2000. Geochemical characteristics of retrometamorphic eclogite in Raobazai, northern Dabie Mountains. Seismology and Geology, 22 (Suppl.): 99-103 (in Chinese with English abstract). [63] Wu, Y. B., Zheng, Y. F., Zhao, Z. F., et al., 2006. U-Pb, Hf and Oisotope evidence for two episodes of fluid-assisted zircon growthin marble-hosted eclogites fromthe Dabie orogen. Geochimicaet Cosmochimica Acta, 70 (14): 3743-3761. doi: 10.1016/j.gca.2006.05.011 [64] Xiao, Y. L., Sun, W. D., Hoefs, J., et al., 2006. Making continental crust through slab melting: Constraints from niobium-tantalum fractionation in UHP metamorphic rutile. Geochimicaet Cosmochimica Acta, 70 (18): 4770-4782. doi: 10.1016/j.gca.2006.07.010 [65] Xiong, X. L., Adam, J., Green, T. H., 2005. Rutile stability and rutile/melt HFSE partitioning during partial melting of hydrous basalt: Implications for TTG genesis. Chemical Geology, 218 (3-4): 339-359. doi: 10.1016/j.chemgeo.2005.01.014 [66] Xu, J., Chen, Y. C., Wang, D. H., et al., 2004. Titaniumm ineralization in the ultrahigh-pressure metamorphic rocks from Chinese Continental Scientific Drilling 100-2000m main hole. Acta Petrologica Sinica, 20 (1): 119-126 (in Chinese with English abstract). [67] Xu, Z. Q., 2004. The scientific goals and investigation progress of the Chinese Continental Scientific Drilling Project. Acta Petrologica Sinica, 20 (1): 1-8 (in Chinese with English abstract). [68] Yao, Y., Ye, K., Liu, J., et al., 2000. Atransitional eclogite to high pressure granulite-facies overprint on coesite-eclogite at Taohang in the Sulu ultrahigh-pressure terrane, eastern China. Lithos, 52 (1-4): 109-120. doi: 10.1016/S0024-4937(99)00087-0 [69] Yu, J. J., Xu, J., Chen, Z. Y., et al., 2006. Trace element geochemistry of rutiles in the eclogites fromthe Chinese Continental Scientific Drilling Project Main Hole. Acta Geologica Sinica, 80 (12): 1835-1841 (in Chinese with English abstract). [70] Zack, T., Kronz, A., Foley, S. F., et al., 2002. Trace element abundances in rutiles fromeclogites and associated garnet mica schists. Chemical Geology, 184: 97-122. doi: 10.1016/S0009-2541(01)00357-6 [71] Zack, T., Moraes, R., Kronz, A., 2004. Temperature dependence of Zr in rutile: Empirical calibration of a rutile thermometer. Contrib. Mineral. Petrol. , 148: 471-488. doi: 10.1007/s00410-004-0617-8 [72] Zhai, M. G., Cong, B. L., 1996. Major and trace element geochemistry of eclogites and related rocks. In: Cong, B. L., ed., Ultrahigh-pressure metamorphic rocks in the Dabieshan-Sulu region of China. Science Press and Kluwer Academic Publishers, Netherlands, 69-89. [73] Zhang, B. H., Liu, Y. S., Zong, K. Q., et al., 2006. Geochemical behavior of trace elements during retrograde metamorphismof eclogites: Study on the retrograde eclogites fromthe CCSD main hole. Acta Petrologica Sinica, 22 (T): 1833-1844 (in Chinese with English abstract). [74] Zhang, R. Y., Hirajima, T., Banno, S., et al., 1995. Petrology of ultrahigh-pressure rocks fromthe southern Su-Lu region, eastern China. J. Metamorph. Geol. , 13: 659-675. doi: 10.1111/j.1525-1314.1995.tb00250.x [75] Zhang, R. Y., Liou, J. G., Cong, B. L., 1994. Petrogenesis of garnet-bearing ultramafic rocks and associated eclogites in the Su-Lu ultrahigh-P metamorphic terrane, eastern China. J. Metamorph. Geol. , 12: 169-186. doi: 10.1111/j.1525-1314.1994.tb00012.x [76] Zhang, Z. M., Xu, Z. Q., Liu, F. L., et al., 2003. Petrology and petrochemistry of UHP metamorphic rocks fromthe Sulu orogenic belt, eastern central China. Acta Geologica Sinica, 77 (4): 478-491 (in Chinese with English abstract). [77] Zhang, Z. M., Xu, Z. Q., Liu, F. L., et al., 2004. Geochemistry of eclogites fromthe main hole (100-2050m) of the Chinese Continental Scientific Drilling Project. Acta Petrologica Sinica, 20 (1): 27-42 (in Chinese with English abstract). [78] Zhang, Z. M., Xu, Z. Q., Xu, H. F., 2000. Petrology of ultrahigh-pressure eclogites fromthe ZK703 drill hole in the Donghai, eastern China. Lithos, 52 (1-4): 35-50. doi: 10.1016/S0024-4937(99)00083-3 [79] Zhang, Z. M., You, Z. D., Han, Y. J., et al., 1996. Petrology, metamorphic process and genesis of the Dabie-Sulu eclogite belt, eastern-central China. Acta Geologica Sinica, 9: 134-156. [80] Zhang, Z. M., Zhang, J. F., You, Z. D., et al., 2005. Ultrahigh-pressure metamorphic P-T-t path of the Sulu orogenic belt, eastern central China. Acta Petrologica Sinica, 21 (2): 257-270 (in Chinese with English abstract). [81] Zhao, R., Zhang, R. Y., Liou, J. G., et al., 2007. Petrochemistry, oxygen isotopes and U-Pb SHRI MP geochronology of mafic-ultramafic bodies from the Sulu UHP terrane, China. J. Metamorph. Geol. , 25 (2): 207-224. doi: 10.1111/j.1525-1314.2007.00691.x [82] Zhao, Z. F., Zheng, Y. F., Gao, T. S., et al., 2006. Isotopic constraints on age and duration of fluid-assisted highpressure eclogite-facies recrystallization during exhumation of deeply subducted continental crust in the Sulu orogen. J. Metamorph. Geol. , 24: 687-702. doi: 10.1111/j.1525-1314.2006.00662.x [83] Zheng, Y. F., Fu, B., Gong, B., et al., 2003. Stable isotope geochemistry of ultrahigh pressure metamorphic rocks fromthe Dabie-Sulu orogen in China: Implications for geodynamics and fluid regime. Earth-Science Reviews, 62 (1-2): 105-161. doi: 10.1016/S0012-8252(02)00133-2 [84] 陈晶, 许志琴, 李旭平, 2005. 超高压榴辉岩中金红石的纳米级微构造特征及成因机制探讨. 岩石学报, 21 (2) 399-404. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200502013.htm [85] 陈振宇, 王登红, 陈毓川, 等, 2006a. 榴辉岩中金红石的矿物地球化学研究及其意义. 地球科学——中国地质大学学报, 31 (4): 533-538, 550. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200604010.htm [86] 陈振宇, 余金杰, 徐珏, 等, 2006b. 榴辉岩中金红石的微量元素电子探针分析. 电子显微学报, 25 (增刊): 295-296. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXV2006S1163.htm [87] 高天山, 郑永飞, 陈仁旭, 2006. 大别造山带黄镇榴辉岩矿物不同类型地质温度计应用和对比. 岩石学报, 22 (7): 1957-1968. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200607023.htm [88] 李晓峰, 陈振宇, 王汝成, 等, 2005. 江苏东海毛发状水晶中金红石矿物学、地球化学特征. 岩石学报, 21 (2): 475-481. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200502022.htm [89] 刘福来, 薛怀民, 许志琴, 等, 2006a. 大别超高压变质带的进变质、超高压和退变质时代的准确限定: 以双河大理岩中榴辉岩锆石SHRI MP U-Pb定年为例. 岩石学报, 22 (7): 1761-1778. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200607002.htm [90] 刘贻灿, 李曙光, 徐树桐, 等, 2006. 北大别片麻岩的超高压变质证据——来自锆石提供的信息. 岩石学报, 22 (7): 1827-1832. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200607008.htm [91] 刘贻灿, 徐树桐, 李曙光, 等, 2000. 大别山北部榴辉岩的地球化学特征和Sr, Nd同位素组成及其大地构造意义. 中国科学(D辑), 30 (增刊): 99-107. https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK2000S1012.htm [92] 刘勇胜, 张泽明, Lee, C. T., 等, 2005. CCSD主孔高Ti榴辉岩非耦合的高Ti、低Nb (Zr) 特征: 对玄武质岩浆房中磁铁矿分离结晶作用的指示. 岩石学报, 21 (2): 339-346. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201503017.htm [93] 戚学祥, 唐哲民, 闫玲, 2006. 中国大陆科学钻探预先导孔(CCSD-PP2) 榴闪岩的地球化学组成及其地质意义. 地球科学——中国地质大学学报, 31 (4): 539-550. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200604011.htm [94] 孙景贵, 胡受奚, 2000. 榴辉岩退变质过程的地球化学研究——以荣成-威海岩段的榴辉岩体为例. 矿物岩石, 20 (1): 17-23. https://www.cnki.com.cn/Article/CJFDTOTAL-KWYS200001003.htm [95] 王大志, 张泽明, 沈昆, 等, 2006. 南苏鲁超高压变质带东海地区富钛榴辉岩及金红石矿的成因. 地质通报, 25 (7): 839-849. doi: 10.3969/j.issn.1671-2552.2006.07.011 [96] 王汝成, 王硕, 邱检生, 等, 2005. CCSD主孔揭示的东海超高压榴辉岩中金红石: 微量元素地球化学及其成矿意义. 岩石学报, 21 (2): 465-474. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200502021.htm [97] 王硕, 王汝成, 邱检生, 等, 2006. CCSD主孔超高压榴辉岩金红石中的矿物包裹体研究. 岩石矿物学杂志, 25 (1): 61-70. doi: 10.3969/j.issn.1000-6524.2006.01.008 [98] 吴元保, 陈道公, 程昊, 等, 2000. 北大别饶拔寨退变质榴辉岩的地球化学特征. 地震地质, 22 (增刊): 99-103. https://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ2000S1013.htm [99] 徐珏, 陈毓川, 王登红, 等, 2004. 中国大陆科学钻探主孔100-2000m超高压变质岩中的钛矿化. 岩石学报, 20 (1): 119-126. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200401009.htm [100] 许志琴, 2004. 中国大陆科学钻探工程的科学目标及初步成果. 岩石学报, 20 (1): 1-8. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200401000.htm [101] 余金杰, 徐珏, 陈振宇, 等, 2006. 中国大陆科学钻探工程主孔榴辉岩中金红石微量元素地球化学特征. 地质学报, 80 (12): 1835-1841. doi: 10.3321/j.issn:0001-5717.2006.12.006 [102] 张斌辉, 刘勇胜, 宗克清, 等, 2006. 榴辉岩退变质过程中的微量元素地球化学行为: 对CCSD主孔退变质榴辉岩的研究. 岩石学报, 22 (T): 1833-1844. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200607009.htm [103] 张泽明, 许志琴, 刘福来, 等, 2003. 南苏鲁造山带的超高压变质岩及岩石化学研究. 地质学报, 77 (4): 478-491. doi: 10.3321/j.issn:0001-5717.2003.04.004 [104] 张泽明, 许志琴, 刘福来, 等, 2004. 中国大陆科学钻探工程主孔(100-2050m) 榴辉岩岩石化学研究. 岩石学报, 20 (1): 27-42. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200401002.htm [105] 张泽明, 张金凤, 游振东, 等, 2005. 苏鲁造山带超高压变质作用及其P-T-t轨迹. 岩石学报, 21 (2): 257-270. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200502001.htm -