《1. 引言》
1. 引言
了解地球内部碳元素分布对于理解全球碳循环[包括大气圈和水圈中的二氧化碳（CO2 ）]至关重要[1,2]。全球碳循环中涉及固体地球的部分，包括沿洋中脊和岛弧火山的CO2 脱气过程，以及通过玄武岩的热液蚀变和地幔橄榄岩的蛇纹石化，在洋壳和地幔中CO2 的储存过程[3,4]。前者将挥发性成分释放到地表和地下储层，而海洋岩石圈俯冲作用是向地幔储层提供挥发性成分的主要入口[5,6]。在海洋沉积物中，进入俯冲带的全球平均CO2 含量约为3.01%±1.44%（质量分数）[7]。据估计，新的洋中脊玄武岩（MORB）包含约0.15%（质量分数）的CO2 ，但在交代和海底风化期间，洋壳600 m以上的CO2 含量可能增加到3%（质量分数），并且整个地壳的平均CO2 含量大约增加0.4%（质量分数）[8]。在最上面的火山岩部分，俯冲岩石圈平均含有约2%（质量分数）的CO2 ，在其余约7 km的地壳中则含有100~5000 ppm（1 ppm=10−6 ）的CO2 [9]。在此过程中，变质再平衡的俯冲板片碳和弧下条件下蛇纹石脱水过程中产生的氧化流体，将有利于碳酸盐化反应，这将使得大部分俯冲板片碳转变为碳酸盐矿物[10]。碳酸盐矿物也可形成于蛇纹石热液碳酸盐化和风化过程中[11]。俯冲板片脱水发生在40~200 km之间，这与蛇纹石、绿泥石、角闪石、辉长岩、云母等各种含水相的稳定性有关[12−15]。
相当大数量的碳可能在俯冲板块脱水后保留下来，并被俯冲至地幔深处[6,16−18]。地幔可能是地球上最大的碳库，并可能超过地表或地下环境中其他所有碳库。菱镁矿（MgCO3 ）被认为是地幔中碳的主要存储者之一[19]。
蛇纹石矿物是含水上地幔大洋板块和俯冲带以上地幔楔的主要成分[20,21]。因此，它们在上地幔水和挥发物元素循环[20,22]，以及构造和地震过程中发挥着重要作用[23−26]。蛇纹石矿物是由1:1的八面体和四面体的层状结构组成。在蛇纹石矿物中，纤蛇纹石[Mg3 Si2 O5 (OH)4 ]，也称为白石棉，其特征为圆柱形纤维晶体，这种结构也赋予石棉相应的结构特性（图1）。该纤维是非手性的，有常规的绕着α轴卷起的纤蛇纹石，或不太常见的绕着b轴卷起的副纤蛇纹石。这种氢氧化镁硅酸盐展现出由四面体[SiO4 ]4– 硅酸盐层（T）组成的一种卷轴状双层结构，沿着与三八面体[Mg3 O2 (OH)4 ]2– 布氏体层（O）相互接触的表面铺展[图1（a）]。当俯冲作用后，纤蛇纹石出现在从地表或近地表环境到一定深度的区域[27]。据报道，沿冷俯冲带纤蛇纹石的温度稳定区域可达600 ℃和（或）100 km深[28]。
《图1》
图1.（a）常温常压条件下沿c轴（上）和a轴（下）观察纤蛇纹石[Mg3 Si2 O5 (OH)4 ]的晶体结构模型。在常温常压（b）及高温高压（5.32 GPa、220 ℃）下与CO2 反应后（c）纤蛇纹石的扫描电镜（SEM）图像。
为理解CO2 与俯冲矿物之间的相互作用，我们以CO2 为传压介质，完成纤蛇纹石的高温高压X射线衍射（XRD）实验，温度、压力分别为500 ℃和5.0(1) GPa。
《2. 材料与方法》
2. 材料与方法
纤蛇纹石（采集自加拿大魁北克，来自史密森研究所、编号为NMMH 156082的样品）的高温高压实验在韩国浦项市加速器实验室（Pohang Accelerator Laboratory, PAL）10C 线站和美国阿贡国家实验室先进光子源（Advanced Photon Source, APS）16BMD 线站进行。在16BMD线站，使用光束波长和直径分别为0.4959(1)Å和90 μm的单色X射线束，以及用来测量XRD模式的MAR345探测器，样品到探测器的距离约为299 mm，是为了在300 s的曝光时间下2θ的角度可以覆盖到大约30°。
改进的MerrillBassett型金刚石压腔（diamond anvil cell, DAC）配备700 石以及碳化钨垫块用于高压实验。初始厚度为700 μm的不锈钢片被预压至130 μm，并通过电火花打孔。纤蛇纹石粉末样品与红宝石压标一起被放入垫片孔中。位于金刚石压腔里面的干燥样品在常温常压下被首次采集。随后，在液氮低温冷却下，固体CO2 被装入已填充好粉末样品的金刚石压腔样品腔中，并密封至第一个压力点。对于非原位加热，金刚石压腔被放置在烘箱中1 h，然后冷却到室温。在原位加热时，我们使用了由环形陶瓷加热器组成的活塞圆筒金刚石压腔。直径为0.3 mm、电阻为0.6 Ω的铂丝被用作加热元件。
通过检测红宝石片R1发射谱线上的位移来测量样品在金刚石压腔中的压力（精度为±0.05 GPa）[29]。在每一压力测量点，平均大约需要10 min，样品在金刚石压腔中的压力才能够稳定下来。在收集每个粉末衍射图谱后，压力大约增加0.5 GPa。
利用高压数据和通用结构分析系统（GSAS）软件进行拟合，通过一系列全剖面拟合程序，得出了晶胞长度和体积随压力变化的关系[30,31]。背景通过一个系数≤20的Chebyshev多项式拟合，采用Thompson等[32]提出的pseudoVoigt剖面函数对观测到的布拉格峰进行拟合。然后使用在2.5(1) GPa 和170 ℃加热1 h后降至常温常压下测得的数据进行精修。通过对四面体和八面体原子和氧原子构架分别进行分组，得到了结构模型的各向同性位移因子[33]。在Si−O四面体和Mg−O八面体中，T−O键和O−O键距离分别用目标值（1.640±0.001）Å和（2.692±0.005）Å、（2.070±0.001）Å和（2.001±0.005）Å进行几何软约束（表1）。纤蛇纹石的（等温）体积压缩用体弹模量K0 来描述（K0 = 1/β = −V·∂P/∂V，其中，β是等温压缩系数），K0 通过三阶BirchMurnaghan状态方程拟合得出[34]，利用的是EOS-fit 7.0程序[35]和基于压力（P）、体积（V）不确定性加权后的数据。
采用日本延世大学（Yonsei University）JSM7001F （JEOL, Japan）型场发射扫描电子显微镜（field emission scanning electron microscope, FESEM）对压力温度处理前后试样的形貌进行了观察。
《表1》
表1 常温常压和高温高压下纤蛇纹石的内部原子距离和角度
（续表）
《3. 结果与讨论》
3. 结果与讨论
以CO2 作为压力传压介质，在2.5(1) GPa下于170 ℃加热1 h后，同步辐射XRD图谱显示纤蛇纹石的结构开始发生剧烈变化（图2）。这些XRD图谱的变化包括固体CO2 峰的消失和菱镁矿峰的出现，这表明压力和温度导致了CO2 传压介质与样品发生了（部分）反应[图2（a）]。有趣的是，在部分碳酸盐反应后，初始的纤维状纤蛇纹石转变为土状[图1（b）、（c）]。
《图2》
图2. 在非原位（a）和原位（b）高温高压下，以二氧化碳为压力传压介质，在PAL 10C线站和APS 16BMD线站所测得的纤蛇纹石同步辐射XRD图谱中的压力引起的变化。Mag：菱镁矿（MgCO3 ）。
利用常温常压和2.5 GPa 和170 ℃的温度下加热1 h后采集到的粉末XRD数据精修得到纤蛇纹石结构模型。对于合成样品和天然样品，纤蛇纹石的空间群结构分别为Cc和C2/m [36,37]。我们的XRD图谱表明纤蛇纹石的空间群为Cc，常温下纤蛇纹石的组成为Mg3 Si2 O5 (OH) 4 ，晶胞常数α = 5.368(1) Å，b = 9.201(2) Å，c = 14.747(4) Å 且 β = 96.25(5)°，最终的一致性因子为Rp = 1.64%，Rwp = 2.66%，且 x2 = 13.35（表2、表3）。沿α轴，SiT1和SiT2四面体的角度为65.72(1)°，在M1位置，Mg八面体中OH2和OH3羟基群的距离为2.18(6) Å[图3（a）]。由于OH3不含硅（Si）四面体，因而容易发生脱羟基或碳酸盐化反应。OH1、OH2和OH4之间的夹角为63.91(1)°，所有这些位点的占有率都设置为1.0[图3（a）]。当M1位点Mg发生碳酸盐化反应时，M1位点周围的八面体受到影响。在2.5 GPa和170 ℃加热1 h后，纤蛇纹石的组成变为Mg2.4(1) Si2 O5 (OH) 2.4(1) ，空间群依旧为Cc，晶胞常数α = 5.355(2) Å，b = 9.178(2) Å，c= 14.535(2) Å，且 β = 96.25(5)°，最终的一致性因子为Rp = 1.98%，Rwp = 3.09%和x2 = 15.59。
《表2》
表2 常温常压和高温高压下纤蛇纹石的结构参数
《表3》
表3 常温常压和高温高压下纤蛇纹石的原子位置、占位、位移参数
（续表）
（续表）
x, y, and z indicate relative position a specific atom in the unit cell; x is fraction of α axis, y is fraction of b axis, and z is fraction of c axis. Occu. refers to occupancy, which indicates what fraction of a site is occupied by a specific atom. U iso indicates isotropic thermal parameter.
相比于常温下的模型，M1和 OH3占位减少了大约50(1)%，反映了菱镁矿的形成。结果在Mg八面体中M1占位OH2与OH3的距离减小至1.99(1) Å，同时OH1、OH2、OH4之间的夹角沿a轴增大至65.08(1)° [图3 （b）]。我们怀疑这是因为M1占位结合OH3占位，而M2和M3占位连接到硅四面体（图3）。在菱镁矿形成时，MgM1八面体占位减少到0.49(1)，并且羟基群占位分别减少OH1 = 0.2(1)，OH2 = 0.4 (2)，OH3 = 0.5(1)和OH4 =0.8(4) [图3（b）]。
《图3》
图3. Rietveld拟合并推导出纤蛇纹石在常温常压条件下（a）和在2.5 GPa、170 ℃加热后（b）的结构细节。
根据纤蛇纹石脱羟基和碳酸盐化的机理模型[38,39]，我们推测在2.5 GPa 和 170 ℃加热1 h后亚稳态纤蛇纹石和菱镁矿形成时发生了质子转移反应。初始的结构水通过质子转移最外层羟基群（OH1和OH2 [40]），从纤蛇纹石中T−O单元释放出来。这个反应的发生是以消耗最内层羟基群OH3的质子为代价的。然后质子会不断向上移动，这样就更容易填满位于脱氢氧原子左边的空位。在此中间态，晶胞重新排布将可能产生亚稳态单脱羟基纤蛇纹石中间态，从而导致晶胞化学式变为Mg2.4(1) Si2 O5 (OH)2.4(1) ，如反应式（1）所示。
这个反应持续至520 ℃，但660 ℃纤蛇纹石开始出现热分解。根据前人的研究[41]，我们认为亚稳态纤蛇纹石将在520 ℃以上分解。
我们的实验条件与苏门答腊岛南部和琉球俯冲带的冷俯冲环境相当[42]（图4），纤蛇纹石出现至2 GPa，而在2.5 GPa压力下170 ℃加热后菱镁矿和亚稳态纤蛇纹石开始形成，并持续至5 GPa。我们知道，在600 ℃以上纤蛇纹石在5~8 GPa（150 ~ 250 km）之间脱水，并释放约13%（质量分数）的水[43]。根据最近的调查，CO2 在80~120 km的深度流失，然后通过导弧火山作用释放到大气中[44,45]。我们的观察结果表明，纤蛇纹石碳酸盐化作用与这个深度范围非常吻合。
《图4》
图4. 南苏门答腊岛和琉球俯冲带的板块表面（实线）和地表下7 km（莫霍板片，虚线）压力温度路径。正方形符号表示热处理后的高压数据。钻石符号表示原位高压高温实验；实心符号表示亚稳态温石棉和菱镁矿的混合物。
《4. 结论》
4. 结论
蛇纹石是一种相对稳定的矿物，它存在于较广的温度和压力范围，从地球表面到俯冲带深度[46]。在本次工作中，我们报道了CO2 与纤蛇纹石的相互作用，研究压力、温度最高至5.0(1) GPa和520 ℃。在80 km深度冷俯冲条件下，常温常压下的纤蛇纹石经部分碳酸盐化作用转变为亚稳态纤蛇纹石和菱镁矿，这与引起岛弧火山作用的碳循环过程相一致。有趣的是，纤蛇纹石特有的纤维形态在碳酸盐化反应后变成了土状。这一发现为石棉在中等压力和温度下解毒提供了可能。进一步的研究正在进行中，以了解亚稳态温石棉在更高的压力和温度条件下的稳定性。
《Acknowledgements》
Acknowledgements
This research was supported by the project Crustal Evolution of Victoria Land, Antarctica and Formative Process of Planets (20140409 and PM18030) funded by the Ministry of Ocean and Fisheries, Korea. The authors also thank the partial supports by the Leader Researcher program (NRF2018R1A3B1052042) of the Korean Ministry of Science and ICT and the NRF grant 2016K1A4A3914691.
《Compliance with ethics guidelines》
Compliance with ethics guidelines
Mihye Kong and Yongjae Lee declare that they have no conflict of interest or financial conflicts to disclose.
京公网安备 11010502051620号
