2021, Vol. 37
2. 中国科学院大学, 北京 100049
2. University of Chinese Academy of Sciences, Beijing 100049, China
金红石由于具有较好的稳定性，可作为重要的副矿物广泛分布在近地表、岩浆岩、低压-高压-超高压变质作用的各种岩石、甚至在月球和陨石中(Goldsmith and Force, 1978; Force, 1991; Banfield and Veblen, 1991; Hwang et al., 2000; Zhang et al., 1998; Luvizotto and Zack, 2009; Hanaor and Sorrell, 2011)。由于类质同像的存在，金红石中常常载入许多其它金属离子(Morris, 1973; Kalfoun et al., 2002; Zack et al., 2002; Carruzzo et al., 2006)，为此金红石作为副矿物出现在各类岩石中，是微量元素特别是Nb、Ta、Ti、Zr、Hf等HFSE的主要载体矿物之一(Rudnick et al., 2000; Kalfoun et al., 2002; Liang et al., 2009; Schmidt et al., 2009)，所以它们在金红石中的含量以及替代关系可以为变质条件和地质演化提供信息(Brenan et al., 1994; Murad et al., 1995; Smith and Perseil, 1997; Zack et al., 2002; Spandler et al., 2003)；此外由于金红石具有较强的H亲和性，它也是名义无水矿物中含水量最高的矿物之一(Force, 1991; Carruzzo et al., 2006; 梁金龙等, 2007)。因此，高温高压下金红石矿物的稳定及其演化对于了解地球深部的水、微量元素的迁移和循环具有重要的科学意义。然而，在此之前从岩石圈到地球深部的运动的基本判断是依据热力学相平衡计算和岩石学观测的结果以及其推论，缺少充分的高温高压条件下的物理化学性质的限定(张雪彤等, 2016)。高温高压条件下金红石晶体结构和水溶解度、OH与之相关的微量元素耦合的数据是严重缺乏的，因为我们目前采集的天然的金红石样品中H2O可能已经受到后期地质事件如退变质反应、热液蚀变等作用的影响而达到新的平衡，从而不能有效地反映其初始组成。
因此，本文旨在通过原位的高温高压实验来模拟俯冲带中金红石晶体结构及其OH的稳定性，从而了解深俯冲带流体中金红石在地质过程中的行为及作用。
1 实验方法 1.1 实验初始样品我们选取来自大别山超高压变质带硬玉石英岩中金红石晶体作为拉曼光谱(Raman)和红外光谱(Infrared Spectroscopy, IR)实验的初始样品。硬玉石英岩主要由石英、硬玉、石榴石及少量的闪石、钠长石、金红石、锆石、榍石等组成(苏文等, 1995)，硬玉中的柯石英表明硬玉石英岩曾俯冲到120km的深度(Su et al., 1996)。金红石常以细小的包体(大小为20~80μm)或粒间基质(100~400μm)产出(苏文等, 1995; 高晓英, 2019)。金红石的U-Pb年龄为224~235Ma表明其超高压变质年龄(Gao et al., 2015; 高晓英, 2019)。对硬玉石英岩样品进行单矿物分选，获得金红石单矿物的晶体，颗粒大小为110~340μm，表明这些晶体主要来源于硬玉石英岩基质中的金红石。这些金红石作为本次实验研究的初始材料。从中挑选成分均一的金红石单晶，将单晶磨制、抛光(厚度约为16~40μm)，待实验备用。
电子探针分析在中国科学院地质与地球物理研究所电子探针与扫描电镜实验室进行，使用的是法国CAMECA公司生产的SXFiveFE高分辨场发射电子探针。测试的条件：电压为20kV，束流为200nA，束斑直径为3μm，背景驻留时间为20s。实验使用的标样为R10金红石标样(王娟等, 2017)。分析的元素为Si、Al、Ti、V、Cr、Fe和Nb。
电子探针分析了金红石的主量、微量元素成分(表 1)，微量元素中FeO含量较高，其次为V2O3、Nb2O5。通过电子探针数据及理想化学式的配比计算法，获得了金红石的晶体化学式，为Ti0.9896Fe0.0048V0.0021Nb0.0014 Cr0.0014 Si0.0007Al0.0002O2。
|
表 1 金红石样品电子探针分析结果(wt%) Table 1 Concentrations of rutiles by EPMA (wt%) |
原位室温高压的Raman和IR实验使用的是短柱状金刚石对顶砧压机(short symmetric-type diamond anvil cell, DAC, 图 1a)。金刚石压砧台面的直径为250μm，预压至大约35μm厚的钢片(型号为304)上打孔作为样品腔，样品腔的直径约为120μm。样品腔中除了装载样品外，还装载了用于传递压力的传压介质和用于标定压强的红宝石球(图 1b)。Raman实验使用的传压介质是高温硅油，IR实验使用的传压介质则是KBr，使用传压介质可确保样品腔内的压力达到均一。实验压强的标定使用的是红宝石荧光测压法，即每次加压完毕，均使用HR-2000光谱仪对红宝石的R1荧光谱线进行测试，通过其波长的偏移量计算出样品腔内的压强(Mao et al., 1986)。将金红石样品抛光至约25μm后装入样品腔内进行原位的光谱实验(图 1b)。
|
图 1 短柱状金刚石对顶砧(a)和样品腔(b)照片 Fig. 1 The photos of short symmetric-type diamond anvil cell (a) and the sample chamber (b) |
原位室温高压下的Raman实验是在中国科学院地质与地球物理研究所拉曼与红外光谱实验室进行的，使用的是德国WITec公司生产的alpha 300R共聚焦激光拉曼光谱仪，测试条件: 激光器为532nm，测试之前使用单晶硅片对拉曼光谱进行校正，经校正使单晶硅片的拉曼峰对应520.74cm-1；物镜为三丰公司生产的50倍长工作距离物镜(NA=0.42)，采集所用的激光功率为15mW，光栅为300g/mm，束斑5μm，采集时间为20s，累积次数为20。
原位室温高压下的IR实验是在中国科学技术大学国家同步辐射实验室的红外谱学和显微成像线站(BL01B)进行，使用的是德国Bruker公司VERTEX-70v红外光谱仪与HYPERIONTM-3000傅里叶变换红外显微镜。测试时实验室的温度为23℃，湿度为23%，分束器是KBr，使用的是液氮冷却的MCT-A检测器，样品和背景的扫描次数均为128次，样品和背景计数分别为1600和3160，分辨率是4cm-1，束斑为5μm，测量的波数范围是4000至650cm-1，覆盖了金红石中O-H伸缩模式的振动特征。每个压力下的红外光谱的测试中，我们均记录了其背景光谱，并在最终的结果中减去了背景光谱所造成的影响。使用OPUS软件进行了包括分峰、拟合等的光谱处理和分析。
1.3 原位高温高压下的Raman和IR实验原位高温高压实验使用的是装配有Ⅱ型金刚石压砧的BX90型金刚石对顶砧压机(BX90-type DAC，图 2a)。金刚石压砧台面的直径为500μm，通过在铼片上打孔形成样品腔，该孔提前被预压至约40μm厚，样品腔直径约230μm。样品腔中，除了样品还装载了用于传递压力的传压介质(高温硅油用于拉曼光谱，氩气用于红外光谱实验)和用于标定压力的红宝石球。压强的标定与室温高压实验的方法一致。
|
图 2 BX90型DAC高温高压实验装置示意图(a)及电阻炉(b)和温控仪(c) Fig. 2 Schematic photo of BX90-type DAC external heating device (a) with external heater (b) and schematic of the temperature controller (c) |
高温高压实验中的高温使用的是DAC外加温电阻进行加温(图 2b)，由一种特殊的陶瓷和Cr-Ni电阻丝构成的炉子(图 2b)，将电阻炉安装在BX-90压机内、金刚石压砧的周围(图 2a, b)，样品腔内的温度是通过k型热电偶(Omega Chal-005)测量的，热电偶用耐高温水泥固定并放置在金刚石压砧侧壁近样品处，并通过控制加载电流来对温度进行精确地调控(图 2c, 苏文等, 2019)，以实现在原位实验过程中进行实时的控温和测温。高温高压实验中，温度的范围为25~700℃，升温速率为5℃/min，降温速率为10℃/min，停留时间5min后进行样品的测试。
原位高温高压下的Raman实验是在中国科学院地质与地球物理研究所拉曼与红外光谱实验室进行的，使用的是德国WITec公司生产的alpha 300R共聚焦激光拉曼光谱仪，测试条件同原位室温高压下的Raman实验。
原位高温高压下的IR实验是在中国科学院地质与地球物理研究所拉曼与红外光谱实验室进行的，使用的是在德国Bruker公司VERTEX-70v红外光谱仪基础上，所构建的、具有水平光路的傅里叶变换红外显微系统上进行实验(苏文等, 2019)。测试时实验室的温度为21℃，湿度为23%，分束器是KBr，使用的是液氮冷却MCT-A检测器，样品和背景的扫描次数均为512次，样品和背景计数分别为900和330，分辨率是4cm-1，束斑为5μm，测量的波数范围是4000~650cm-1。每个压力下的红外光谱的测试中，我们同样记录了其背景光谱，并在最终的结果中减去了背景光谱所造成的影响。使用OPUS软件进行了包括分峰、拟合等的光谱处理和分析。
2 实验结果 2.1 常温常压下金红石的振动模式 2.1.1 拉曼光谱大别硬玉石英岩中的金红石样品在常温常压下的Raman如图 3所示，显示有三个拉曼特征峰，分别为ν1(239cm-1)、ν2(450cm-1)、ν3(615cm-1)。
金红石型TiO2结构属于D4hl4=P42/mnm空间群，其光学振动模式表示为(Chaves et al., 1974)，
|
因此，金红石的Raman通常有4种振动模式，在测得的金红石样品Raman中，观察到了ν1(B1g)、ν2(Eg)和ν3(A1g)三个特征峰，而前人研究中(Zhang et al., 1998)观察到的826cm-1(B2g)在本次金红石样品中并未出现(图 3)。ν2通常被指派为O-Ti-O的弯曲振动，而ν3则被指派为O-Ti-O的反伸缩振动(肖萍等, 2007)。ν1通常认为是多光子散射的结果(Ma et al., 2007), Hara and Nicol (1979)认为，金红石的晶格是无序的，而ν1的出现正是晶格无序导致的。
|
图 3 常温常压下金红石样品的拉曼光谱 Fig. 3 Raman of rutile at atmospheric pressure and room temperature |
常温常压下大别硬玉石英岩中金红石样品的IR如图 4所示。
|
图 4 金红石样品在室温室压下的红外光谱 Fig. 4 IR of rutile at atmospheric pressure and room temperature |
使用PeakFit软件(高斯)进行拟合后发现，样品中只有一个OH的红外特征峰a，为3320~3325cm-1。该特征峰是由Al3++H+⇌Ti4+的取代引起的(Bromiley et al., 2004; Guo, 2017)。
2.2 原位室温高压下金红石的振动模式 2.2.1 拉曼光谱对大别硬玉石英岩样品中的金红石进行原位室温高压下Raman实验，对样品先加压至最高压力点、后卸压，观察其拉曼谱峰随压力的变化，实验结果如图 5所示，实验所达到的最高压力为43.67GPa(图 5、表 2)。
|
图 5 原位室温高压下金红石样品的拉曼光谱 (a)从0.0001Gpa到43.67GPa的加压过程，(b)从43.67GPa到0.0001GPa的降压过程 Fig. 5 Changes of Raman Spectrum of rutile sample under high pressure (a) increasing pressure process from 0.0001GPa to 43.67GPa; (b) decreasing pressure process from 43.67GPa to 0.0001GPa |
|
|
表 2 金红石样品室温高压下拉曼光谱 Table 2 Raman spectrum of Rutile at high-pressure and room temperature |
如前所述，样品在常温常压下具有ν1(239cm-1)、ν2(450cm-1)、ν3(615cm-1)三个拉曼特征谱峰，当开始加压时，三个特征谱峰随着压力的变化向高波数方向偏移(图 5a、表 2)，当压力增加至23.43GPa，拉曼谱峰发生变化，在(199cm-1, ν4)、(232cm-1, ν5)、(324cm-1, ν6)和(448cm-1, ν7)附近出现新的拉曼谱峰(图 5a、表 2)，表明压力导致了金红石的相变；继续加压至34.98GPa，再次出现新的拉曼谱峰，其位于(741cm-1, ν8)(图 5a、表 2)，意味着金红石的第2次相变的出现。随着压力持续增加至43.67GPa，金红石样品的拉曼谱峰没有发生改变(图 5a)，表明金红石样品在34.98GPa之后没有出现进一步的相变。
降压过程中，当压力降至15.50GPa时开始发生相变，以(141cm-1, ν9)(图 5b、表 2)的出现为标志；卸压至5.05GPa，新的拉曼谱峰(180cm-1，ν10)出现(图 5b、表 2)，当卸至室压，其拉曼谱峰未再发生改变(图 5b)。
2.2.2 红外光谱对大别硬玉石英岩中金红石样品进行了原位室温下的加压IR实验，实验的最高压力为37.28GPa，实验结果如图 6所示。
|
图 6 室温高压下金红石的红外光谱 (a)从0.0001GPa到37.28GPa的加压过程；(b)从37.28GPa到0.0001GPa的降压过程 Fig. 6 Changes of IR of rutile under high pressure (a) increasing pressure process from 0.0001GPa to 37.28GPa; (b) decreasing pressure process from 37.28GPa to 0.0001GPa |
加压过程中，金红石样品的OH特征谱峰a(3323cm-1)随着压力的增加缓慢地向低波数方向偏移，并且强度降低(图 6a、表 3)。当压力增加至32.69GPa时，红外特征峰a开始向高波数方向偏移，强度继续降低(图 6a、表 3)。继续加压至实验的最高压力37.28GPa，并未有新的谱峰的出现或消失。在这个加压过程中，金红石样品的红外特征峰仅波数发生偏移，并没有特征峰的出现或消失。
|
|
表 3 金红石样品在室温高压下样品的红外光谱峰位 Table 3 IR of rutile at high-pressure and room temperature |
降压过程中的红外特征峰a则开始向低波数方向偏移(图 6b、表 3)，且强度逐渐增加；当卸压至12.51GPa时，金红石OH的红外特征峰a开始向高波数方向偏移，其强度继续增加(图 6b、表 3)。直至卸压至室压，红外特征峰a也并未消失。样品从37.28GPa卸压至室压的过程中，仅有红外特征峰a的波数发生偏移，并未有特征峰的消失或出现。
2.3 原位高温高压下金红石的振动模式 2.3.1 拉曼光谱将金刚石对顶砧技术与共聚焦激光拉曼光谱仪联用，进行原位的高温高压实验。本次实验先将大别硬玉石英岩中金红石样品加压至4.08GPa后再进行原位的高温高压实验，实验结果如图 7所示。
|
图 7 原位高温高压下金红石的拉曼光谱 (a)从室温到700℃的升温过程；(b)从700℃到室温的降温过程 Fig. 7 Changes of Raman of rutile under high-temperature and high-pressure (a) heating process from room temperature to 700℃; (b) cooling process from 700℃ to room temperature |
样品的ν1、ν2、ν3三个拉曼谱峰在加压至4.08GPa后整体向高波数方向偏移(蓝移)，此特征与室温高压下样品的拉曼谱峰的变化趋势一致(图 5a)，而此时并未出现新的拉曼谱峰，故仍为金红石型TiO2结构。在此压力下，对BX 90压机中的样品腔进行加温。随着温度的升高，金红石的拉曼特征谱峰ν1在232~258cm-1附近波动(图 7a、表 4)，仅表现为强度降低但并无明显的蓝移或红移。而ν2和ν3两个拉曼特征谱峰均向低波数方向偏移(红移)，至实验的最高温度700℃时，也并未出现新的拉曼谱峰(图 7a)，表明金红石样品未发生相变。降温过程中，ν1在243~258cm-1附近波动(图 7b、表 4)，ν2和ν3则向高波数方向偏移(蓝移)，直至接近实验初始的拉曼峰位(图 7b、表 4)。
通过图 7可以发现，大别硬玉石英岩中的金红石样品的三个拉曼特征峰在温度影响下的行为基本一致，为了进一步观察其在高温高压下的行为，我们对三个拉曼特征谱峰峰位的具体偏移进行了统计和计算，如表 4所示。
|
|
表 4 原位高温高压下金红石样品的拉曼光谱峰位 Table 4 The peak position of Raman of rutile at high-pressure and high-temperature |
通过计算发现，在加温过程中，拉曼特征峰ν2的变化速率比ν1和ν3的变化速率要快，为-0.063cm-1/℃。在降温过程中，ν1的变化速率快于ν2和ν3，为-0.013cm-1/℃。
2.3.2 红外光谱将金刚石对顶砧与IR相结合，进行原位的高温高压实验(苏文等, 2019)。本次实验先将大别硬玉石英岩中金红石样品加压至4.01GPa后再进行原位的高温高压实验，以此模拟金红石中OH在高温高压下的行为。实验结果如图 8所示。
|
图 8 原位高温高压下金红石的红外光谱 (a)从室温到700℃的升温过程；(b)从700℃到室温的降温过程 Fig. 8 Changes of IR of rutile under high-temperature and high-pressure (a) heating process from room temperature to 700℃; (b) cooling process from 700℃ to room temperature |
通过实验发现，该样品的OH红外特征峰a会随着温度的增加向低波数方向偏移，且其强度会随着温度的增加而降低(图 8a)，当实验达到最高温度700℃时，该样品的OH红外特征峰也并未完全消失。当样品开始降温时，其OH红外特征峰又会随着温度的降低向高波数方向偏移(图 8b)，特征峰的强度随着温度的降低有一定程度的增加。在整个实验过程中均未发生OH红外特征峰的增加或消失。
3 讨论 3.1 金红石在温度和压力下的结构转变大别硬玉石英岩中金红石样品在室温、加压过程中经历了两次相变，其中第一次相变以ν4、ν5、ν6和ν7的出现为标志(图 5a、表 2)，说明样品已经由金红石型TiO2结构转变为斜锆石型TiO2结构(Mammone et al., 1980; Lagarec and Desgreniers, 1995; 张红等, 2006; 肖万生等, 2007)，其相变的压力为23.43GPa(指示的深度约为670km)。第二次相变的标志为新的ν8拉曼谱峰的出现(图 5a、表 2)，标志着金红石样品已经由斜锆石型TiO2结构转变为萤石型TiO2结构(Zhou and Hu, 2011)，其相变的压力为34.98GPa(深度超过870km)。在降压过程中，当压力降低至15.50GPa时，样品出现新的拉曼谱峰ν9，卸压至5.05GPa，谱峰ν10出现(图 5b、表 2)，表明样品由萤石型TiO2结构完全转变为ɑ-PbO2型TiO2结构(Mammone et al., 1980; Lagarec and Desgreniers, 1995; 张红等, 2006), 而且卸压完毕后ɑ-PbO2型TiO2结构是唯一能够稳定存在的相。因此金红石在经历俯冲(至250km深度)折返后可能会以ɑ-PbO2型TiO2结构存在。孟大维等(2008)在大别山石马和双河地区超高压变质岩矿物中发现的纳米级ɑ-PbO2型TiO2晶体，表明ɑ-PbO2型TiO2结构可以在常温常压下稳定存在，同时也指示大别超高压变质岩片曾俯冲到约为200~250km深度。
单一的压力变量无法真实地反映深俯冲过程中金红石晶体结构及其OH的转变，因此应同时考虑温度和压力的高温高压实验来模拟金红石在俯冲过程中的变化特征。
对大别硬玉石英岩中的金红石样品ν1、ν2、ν3三个拉曼特征峰在高温高压下拉曼谱峰位置的变化情况进行了对比(图 9)。
|
图 9 金红石样品在高温高压条件下拉曼谱峰的变化趋势 (a)从室温到700℃的升温过程；(b)从700℃到室温的降温过程，其中线表示变化趋势线 Fig. 9 Changes of Raman of rutile under high-temperature and high-pressure (a) heating process; (b) cooling process, the line represents the change trend |
通过图 9可以看出，在加温过程中，ν2和ν3谱峰总体向低波数方向偏移，但其变化趋势是不同的。为了更精确的反应变化的速率，我们对金红石的不同拉曼特征峰在高温高压条件下的变化趋势进行了拟合，用ki值来表征其变化(表 5)。
|
|
表 5 金红石样品在高温高压下拉曼光谱的变化 Table 5 The changes in Raman of Rutile at high-pressure and high-temperature |
从表 5可以看出，以450℃为分界，ν2和ν3在加温至450℃前的变化速率要小于温度大于450℃时。说明在室温到450℃以前ν2代表的O-Ti-O弯曲振动和ν3代表的O-Ti-O的反伸缩振动能量减弱速率要小于450~700℃时速率。
通过原位高温高压Raman实验可以证明，大别硬玉石英岩中的金红石从常温常压加压至4.08GPa的过程中ν1、ν2、ν3三个拉曼谱峰向高波数方向偏移(蓝移，图 7a、表 4)，但因其压力并未达到金红石向斜锆石型TiO2结构转变的压力(23.43GPa，图 5a、表 2)，故此时样品仍为金红石结构。在4.08GPa压力下开始对样品加温，从室温加温到700℃的过程中(P=4.08GPa)，ν2和ν3向低波数方向偏移(红移)，在此过程中并未出现新的拉曼特征峰，故其晶体结构并未发生变化(图 7a)。Raman红移代表组成金红石晶格的化学键能量变小(Hara and Nicol, 1979)，因此，ν2代表的O-Ti-O弯曲振动和ν3代表的O-Ti-O的反伸缩振动在温度升高的过程中，能量减弱，金红石样品晶格有序度降低。在降温过程中，ν2和ν3向高波数方向偏移(蓝移)，说明在此过程中，样品晶格中的O-Ti-O弯曲振动和O-Ti-O的反伸缩振动能量增强，样品有序度增加。
3.2 高温高压下金红石中OH的稳定性在金红石结构发生变化的同时，金红石中OH红外特征谱峰也随着压力的变化而变化。
根据Beer-Lambert公式，A=Kbc。
其中A代表吸光度，可以用OH红外吸收峰的面积表示；K是红外吸收系数，仅与吸收材料的性质和入射光的波长有光，在相同实验中，K为常数；b代表材料的厚度；c是OH含量。可以发现，样品的含水量c与OH红外吸收峰的峰面积A成正比。因此，我们使用OPUS软件对样品在室温高压实验和高温高压实验中OH红外吸收峰的谱峰面积进行积分，获得其变化趋势图(图 10)，以期反映其含水量的变化。
|
图 10 金红石样品红外谱峰面积与压力关系图解 (a)从0.0001GPa到37.28GPa的加压过程；(b)从37.28GPa到0.0001GPa的降压过程 Fig. 10 Diagram of the relationship between the peak area of IR and pressure of rutile (a) increasing pressure process from 0.0001GPa to 37.28GPa; (b) decreasing pressure process from 37.28GPa to 0.0001GPa |
图 10显示了金红石中OH红外吸收峰面积在室温高压实验中的变化趋势。OH含量随着压力的增加总体呈减少趋势，在23.97GPa和32.69GPa处出现拐点。其中，压力从室温增加至23.97GPa时，金红石中OH含量降低，此过程中样品的红外特征峰a向低波数方向偏移，并且强度降低(图 6a)，此压力与金红石转变为斜锆石型TiO2结构的压力接近(23.43GPa，图 5a)。因此在这个过程中OH含量的降低可能是由于结构改变造成的。当压力从23.97GPa增加到32.69GPa时，金红石OH含量有一定程度的增加，此时样品的红外特征峰a继续向低波数方向偏移，强度继续降低(图 6a)，此压力与金红石转变为萤石型TiO2结构的压力接近(34.98GPa，图 5a)，因此，OH含量增加同样也可能是由于结构的改变造成的。当压力继续升高至最高压力(37.28GPa)的过程中，金红石中OH含量减少，值得注意的是，在此过程中，金红石的红外特征峰a开始向高波数方向转变。
因此，在压力增加过程中，金红石中水含量总体呈减少趋势，但直至37.28GPa(> 600km深度)，金红石中的水也并未完全脱去，证明金红石可以将水携带至深部地幔；卸压过程中，金红石中水含量总体呈增加趋势，但卸至室压，金红石中的水也并未恢复实验前的含量，说明经历过快速俯冲折返的金红石样品的水含量并不能代表金红石形成时的水含量。
图 11显示了金红石中OH红外吸收峰面积在原位高温高压实验中的变化趋势图。大别硬玉石英岩中金红石中OH含量随温度的变化呈规律性变化。在加温过程中，随着温度的升高，金红石中OH含量整体呈减少的趋势，在此过程中，样品的红外特征谱峰a随着温度的升高，整体向低波数方向偏移，并且其强度逐渐降低(图 8a)，表明OH含量的减少。但直至实验的最高温度700℃，OH的红外特征峰也并未完全消失，说明金红石可以携带水至地球深部约120km处。在降温过程中，其红外特征峰的峰面积整体呈增加的趋势，在此过程中，金红石的红外特征峰a整体向高波数方向偏移，强度随着温度的降低升高(图 8b)，表明OH含量的增加。但增加的速率要远远小于升温过程中降低的速率，并且直到温度降低至室温，其OH含量也并未恢复至高温高压前的含量，说明金红石在经历俯冲折返过程后的水含量无法代表初始时的金红石中水含量。
|
图 11 金红石样品红外谱峰面积与温度关系图解 (a)从室温到700℃的升温过程；(b)从700℃到室温的降温过程 Fig. 11 Diagram of the relationship between the peak area of IR and temperature of rutile (a) heating process from room temperature to 700℃; (b) cooling process from 700℃ to room temperature |
通过室温高压和高温高压下的Raman和IR实验，我们可以得出以下结论：
(1) 大别硬玉石英岩中金红石的结构在室温加压的过程中经历了两次相变，其中斜锆石结构型TiO2结构可以稳定存在于23.43~34.98GPa(深度为670~870km)范围内，而萤石型TiO2结构则可以在压力超过34.98GPa(>870km深度)的条件下稳定存在。降压过程中则由萤石型TiO2结构转变为ɑ-PbO2型TiO2结构，并在常温常压下稳定。暗示样品可能是在经历超深俯冲、折返后唯一稳定存在的相为ɑ-PbO2型TiO2相，标志着俯冲深度超过250km，因此ɑ-PbO2型TiO2结构的存在对俯冲深度有一定的指示意义。
(2) 金红石样品中的OH含量随着压力的增加总体呈减少趋势。但至最高压力37.28GPa(> 870km深度)金红石中的水也并没有完全脱去，证明金红石可以将水携带至深部地幔，并在俯冲过程中，不断释放H。压力降低时，金红石内部的OH含量有一定程度的增加，但并未恢复其初始状态。
(3) 金红石中OH含量随着温度的增加而减少，至实验最高温度700℃，金红石中的OH也并未完全消失，证明金红石中的H也可以在P=4.01GPa和T=700℃的温压条件下稳定存在。在降温过程中，金红石OH含量有所增加，但并未恢复至初始的状态，说明大别硬玉石英岩在经历俯冲超过120km的超高压变质作用后，再折返的金红石中的OH含量无法代表金红石的初始含量。
致谢 感谢中国科学院地质与地球物理研究所的毛骞、中国科学技术大学国家同步辐射实验室的戚泽明分别在电子探针和红外光谱实验中的帮助。
Banfield JF and Veblen DR. 1991. The structure and origin of Fe-bearing platelets in metamorphic rutile. American Mineralogist, 76(1-2): 113-127
|
Brenan JM, Shaw HF, Phinney DL and Ryerson FJ. 1994. Rutile-aqueous fluid partitioning of Nb, Ta, Hf, Zr, U and Th: Implications for high field strength element depletions in island-arc basalts. Earth and Planetary Science Letters, 128(3-4): 327-339 DOI:10.1016/0012-821X(94)90154-6
|
Bromiley GD, Keppler H, McCammon C, Bromiley FA and Jacobsen SD. 2004. Hydrogen solubility and speciation in natural, gem-quality chromian diopside. American Mineralogist, 89(7): 941-949 DOI:10.2138/am-2004-0703
|
Carruzzo S, Clarke DB, Pelrine KM and MacDonald MA. 2006. Texture, composition, and origin of rutile in the South Mountain Batholith, Nova Scotia. The Canadian Mineralogist, 44(3): 715-729 DOI:10.2113/gscanmin.44.3.715
|
Chaves A, Katiyan KS and Porto SPS. 1974. Coupled modes with A1 symmetry in tetragonal BaTiO3. Physical Review B, 10(8): 3522-3533 DOI:10.1103/PhysRevB.10.3522
|
Force ER. 1991. Geology of titanium-mineral deposits. In: Force ER (ed. ). Geology of Titanium-Mineral Deposits. America: Geological Society of America, 259: 1-112
|
Gao XY, Zheng YF, Chen YX, Tang HL and Li WC. 2015. Zircon geochemistry records the action of metamorphic fluid on the formation of ultrahigh-pressure jadeite quartzite in the Dabie orogen. Chemical Geology, 419: 158-175 DOI:10.1016/j.chemgeo.2015.10.043
|
Gao XY. 2019. Large-scale flow of metamorphic fluids in a continental subduction zone: Evidence from coesite-bearing jadeite quartzite in the Dabie Orogen. Earth Science, 44(12): 4064-4017 (in Chinese with English abstract)
|
Goldsmith R and Force ER. 1978. Distribution of rutile in metamorphic rocks and implications for placer deposits. Mineralium Deposita, 13(3): 329-343
|
Guo HH. 2017. In-situ infrared spectra of OH in rutile up to 1000℃. Physics and Chemistry of Minerals, 44(8): 547-552 DOI:10.1007/s00269-017-0881-6
|
Hanaor DAH and Sorrell CC. 2011. Review of the anatase to rutile phase transformation. Journal of Materials Science, 46(4): 855-874 DOI:10.1007/s10853-010-5113-0
|
Hara Y and Nicol M. 1979. Raman spectra and the structure of rutile at high pressures. Physica Status Solidi (b), 94(1): 317-322 DOI:10.1002/pssb.2220940137
|
Hwang SL, Shen PY, Chu HT and Yui TF. 2000. Nanometer-size α-PbO2-type TiO2 in garnet: A thermobarometer for ultrahigh-pressure metamorphism. Science, 288(5464): 321-324 DOI:10.1126/science.288.5464.321
|
Kalfoun F, Ionov D and Merlet C. 2002. HFSE residence and Nb/Ta ratios in metasomatised, rutile-bearing mantle peridotites. Earth and Planetary Science Letters, 199(1-2): 49-65 DOI:10.1016/S0012-821X(02)00555-1
|
Lagarec K and Desgreniers S. 1995. Raman study of single crystal Anatase TiO2 up to 70GPa. Solid State Communications, 94(7): 519-524 DOI:10.1016/0038-1098(95)00129-8
|
Liang JL, Sun XM, Xu L, Zhai W, Tang Q and Liang YH. 2007. Structural hydroxyl of rutile in UHP metamorphic rocks from Chinese Continental Scientific Drilling Project (CCSD) and areas nearby and its implication for geodynamics. Geological Review, 53(2): 267-272 (in Chinese with English abstract)
|
Liang JL, Ding X, Sun XM, Zhang ZM, Zhang H and Sun WD. 2009. Nb/Ta Fractionation observed in Eclogites from the Chinese Continental Scientific Drilling Project. Chemical Geology, 268(1-2): 27-40 DOI:10.1016/j.chemgeo.2009.07.006
|
Luvizotto GL and Zack T. 2009. Nb and Zr behavior in rutile during high-grade metamorphism and retrogression: An example from the Ivrea-Verbano Zone. Chemical Geology, 261(3-4): 303-317 DOI:10.1016/j.chemgeo.2008.07.023
|
Ma HL, Yang JY, Dai Y, Zhang YB, Lu B and Ma GH. 2007. Raman study of phase transformation of TiO2 rutile single crystal irradiated by infrared femtosecond laser. Applied Surface Science, 253(18): 7497-7500 DOI:10.1016/j.apsusc.2007.03.047
|
Mammone JF, Sharma SK and Nicol M. 1980. Raman study of rutile (TiO2) at high pressure. Solid State Communications, 34(10): 799-802 DOI:10.1016/0038-1098(80)91055-8
|
Mao HK, Xu J and Bell PM. 1986. Calibration of the ruby pressure gauge to 800kbar under quasi-hydrostatic conditions. Journal of Geophysical Research: Solid Earth, 91(B5): 4673-4676 DOI:10.1029/JB091iB05p04673
|
Meng DW, Wu XL, Sun F, Liu F and Zhang JP. 2008. Identification of α-PbO2-type TiO2 in tadeite quartzite from Shuanghe, Dabie Mountains, China. Earth Science (Journal of China University of Geosciences), 33(5): 706-715 (in Chinese with English abstract) DOI:10.3799/dqkx.2008.087
|
Morris GJ. 1973. The selection, training and development of computer personnel. IFAC Proceedings Volumes, 6(1): 445-447 DOI:10.1016/S1474-6670(17)68292-2
|
Murad E, Cashion JD, Noble CJ and Pilbrow JR. 1995. The chemical state of Fe in rutile from an Albitite in Norway. Mineralogical Magazine, 59(396): 557-560 DOI:10.1180/minmag.1995.059.396.17
|
Rudnick RL, Barth M, Horn I and McDonough WF. 2000. Rutile-bearing refractory Eclogites: Missing link between continents and depleted mantle. Science, 287(5451): 278-281 DOI:10.1126/science.287.5451.278
|
Schmidt A, Weyer S, John T and Brey GP. 2009. HFSE systematics of rutile-bearing eclogites: New insights into Subduction zone processes and implications for the Earth's HFSE budget. Geochimica et Cosmochimica Acta, 73(2): 455-468 DOI:10.1016/j.gca.2008.10.028
|
Smith DC and Perseil EA. 1997. Sb-rich rutile in the manganese concentrations at St. Marcel-praborna, Aosta Valley, Italy: Petrology and crystal-chemistry. Mineralogical Magazine, 61(408): 655-669
|
Spandler C, Hermann J, Arculus R and Mavrogenes J. 2003. Redistribution of trace elements during prograde metamorphism from lawsonite blueschist to eclogite facies: Implications for deep subduction-zone processes. Contributions to Mineralogy and Petrology, 146(2): 205-222 DOI:10.1007/s00410-003-0495-5
|
Su W, Xu ST, Jiang LL and Liu YC. 1995. Characteristics of jadeite-quartzite and associated eclogite in Hanchengzhong-Miaozhuyiu area, Qianshan. Geology of Anhui, 5(2): 7-21 (in Chinese with English abstract)
|
Su W, Xu ST, Jiang LL and Liu YC. 1996. Coesite from quartz-jadeitite in the Dabie Mountains, eastern China. Mineralogical Magazine, 60(401): 659-662 DOI:10.1180/minmag.1996.060.401.12
|
Su W, Liu ZX, Chen F, Gao J and Li XG. 2019. Construction of FT-IR micro-spectroscopy system: In situ studies on reaction-dynamics of material. Acta Petrologica Sinica, 35(1): 252-260 (in Chinese with English abstract) DOI:10.18654/1000-0569/2019.01.20
|
Wang J, Chen Y, Mao Q, Li QL, Ma YG, Shi YH and Song CZ. 2017. Electron microprobe trace element analysis of futile. Acta Petrologica Sinica, 33(6): 1934-1946 (in Chinese with English abstract)
|
Xiao P, Zheng SB, You JL, Jiang GC, Chen H and Zeng H. 2007. Structure and Raman spectra of titanium oxides. Spectroscopy and Spectral Analysis, 27(5): 936-939 (in Chinese with English abstract)
|
Xiao WS, Zhang H, Tan DY, Weng KN, Li YC, Luo CJ, Liu J and Xie HS. 2007. Raman characterization of rutile phase transitions under high-pressure and high-temperature. Spectroscopy and Spectral Analysis, 27(7): 1340-1343 (in Chinese with English abstract)
|
Zack T, Kronz A, Foley SF and Rivers T. 2002. Trace element abundances in rutiles from eclogites and associated garnet mica schists. Chemical Geology, 184(1-2): 97-122 DOI:10.1016/S0009-2541(01)00357-6
|
Zhang H, Xiao WS, Chen JY, Weng KN and Lü GC. 2006. Investigation of structural transformations in rutile (TiO2) at high pressure by Raman spectroscopy. Chinese Journal of Spectroscopy Laboratory, 23(1): 1-4 (in Chinese with English abstract)
|
Zhang XT, Zhang RH and Hu SM. 2016. Fluids in extreme conditions of the earth interior: 1. In situ observation of the high temperature and ultra-high pressure fluids by using diamond anvil cell. Acta Geologica Sinica, 90(8): 1846-1859 (in Chinese with English abstract)
|
Zhang YH, Chan CK, Porter JF and Guo W. 1998. Micro-Raman spectroscopic characterization of nanosized TiO2 powders prepared by vapor hydrolysis. Journal of Materials Research, 13(9): 2602-2609 DOI:10.1557/JMR.1998.0363
|
Zhou ZJ and Hu YF. 2011. First-principles study of rutile-to-fluorite transition. Chinese Journal of Computational Physics, 28(1): 111-118
|
高晓英. 2019. 大陆俯冲带大规模的变质流体活动: 来自大别造山带超高压硬玉石英岩的记录. 地球科学, 44(12): 4064-4017. |
梁金龙, 孙晓明, 徐莉, 翟伟, 汤倩, 梁业恒. 2007. CCSD超高压变质岩中金红石结构水(OH-)的傅里叶变换红外光谱测定及其地球动力学意义. 地质论评, 53(2): 267-272. |
孟大维, 吴秀玲, 孙凡, 刘富, 郑建平. 2008. 大别山硬玉石英岩中发现α-PbO2型TiO2超高压相. 地球科学(中国地质大学学报), 33(5): 706-715. |
苏文, 徐树桐, 江来利, 刘贻灿. 1995. 安徽潜山韩长冲-苗竹园一带硬玉石英岩及其伴生榴辉岩特征. 安徽地质, 5(2): 7-21. |
苏文, 刘振先, 陈菲, 高静, 李晓光. 2019. 构建红外光谱显微系统实验平台——原位模拟物质反应动力学过程. 岩石学报, 35(1): 252-260. |
王娟, 陈意, 毛骞, 李秋立, 马玉光, 石永红, 宋传中. 2017. 金红石微量元素电子探针分析. 岩石学报, 33(6): 1934-1946. |
肖萍, 郑少波, 尤静林, 蒋国昌, 陈辉, 曾昊. 2007. 钛氧化物结构及其拉曼光谱表征. 光谱学与光谱分析, 27(5): 936-939. DOI:10.3321/j.issn:1000-0593.2007.05.026 |
肖万生, 张红, 谭大勇, 翁克难, 李延春, 罗崇举, 刘景, 谢鸿森. 2007. 金红石高温高压相变的Raman光谱特征. 光谱学与光谱分析, 27(7): 1340-1343. |
张红, 肖万生, 陈晋阳, 翁克难, 律广才. 2006. 高压下金红石的拉曼光谱分析. 光谱实验室, 23(1): 1-4. |
张雪彤, 张荣华, 胡书敏. 2016. 地球内部高温超高压流体的水热金刚石压砧原位观测. 地质学报, 90(8): 1846-1859. |