断裂带中方解石脉显微结构、地球化学和年代学研究：藏南亚东裂谷带古地震记录*

李春锐 李海兵 潘家伟 刘栋梁 白明坤 刘富财 张进江

1. 中国地质科学院地质研究所，自然资源部深地动力学重点实验室，北京 100037

2. 北京大学地球与空间科学学院，北京 100871

3. 南方海洋科学与工程广东省实验室(广州)，广州 511458

断层以及断层附近的破裂带作为流体运移的通道，在地震活动期间能够沟通不同层次流体，如大气水、地壳流体、幔源流体(Menziesetal., 2014；Barberioetal., 2017)。地震过程中富含CO2流体运移时，由于岩石破裂导致的流体压力降低，以及不同来源流体混合时导致的矿物溶解度变化，在断层的阶步以及破碎带裂缝中会产生矿物沉淀，形成方解石脉(Dolejš and Manning，2010；De Filippisetal., 2013；Brogietal., 2014)。方解石脉的显微结构可以指示其与地震事件的关系：同构造期沉淀脉体结构为块状、拉长块状和扩展状，其中同震阶段形成的脉体由于裂隙开启过程中流体压力快速下降导致的高矿物饱和度，一般具有块状等粒结构(Verhaertetal., 2003；Boullieretal., 2004；Uysaletal., 2007；Watanabeetal., 2008；Nurieletal., 2012)；震间期压溶-沉淀，常见结构为“压力影”和超压泥页岩层中发育的纤维脉(Eyaletal., 1992；Bonsetal., 2012; Zhangetal., 2015)；震间期开放裂缝沉淀，常见类型为浅表暴露面的泥质薄层(Flottéetal., 2001；Verhaertetal., 2004；Nurieletal., 2011)。沿断层发育的方解石脉作为地震期间断层活动和裂缝开放的记录(Boullieretal., 2004；Williamsetal., 2017)，其地球化学示踪剂如同位素，REY(稀土元素和钇)模式，可以用于推断发震时期以及地壳流体循环路径(Nolletetal., 2005；Uysaletal., 2011；Bergmanetal., 2013)。由于非地震活动成因的脉体流体来源为相邻围岩，因此脉体地球化学和同位素多具有与围岩相似的特征。地震活动相关脉体具有更大的外部流体来源贡献，地球化学特征与围岩常具有明显差异。因此，对断层方解石脉的研究有助于更好地理解地震行为、地震周期以及地壳内部流体循环模式。

青藏高原的东西向伸展产生了一系列南北向裂谷(图1)，是藏南地区最显著的活动构造(Molnar and Tapponnier，1978；Tapponnieretal., 1981；Armijoetal., 1986；Kapp and Guynn，2004)。亚东-古露裂谷是南北向裂谷系中规模最大的裂谷，地震活动性最强、发震频率最高(Wuetal., 2011)。自公元1264年以来，亚东-古露裂谷内部及周缘共发生了17次Mw>6.0的地震，其中14次发生在裂谷北段和中段，即雅鲁藏布江至古露之间(Chevalieretal., 2020)。裂谷南段第四纪地震活动研究同样缺乏，仅在帕里地堑有两处报道，分别为地震崩滑体以及地表破裂带(Pengetal., 2018；Haetal., 2019a)。作为藏南地区规模最大的裂谷系，亚东-谷露裂谷南段延伸约200km，在裂谷东侧的基岩内有广泛的断层出露(图2)。裂谷断层附近发育的方解石脉能够为同震断裂活动研究提供依托。

图1 藏南地区正断层分布和亚东-谷露裂谷位置

图2 亚东裂谷区域地质简图(据Zhang et al., 2012; 刘文灿等，2014)

本文对亚东裂谷断层涅如堆段和帕里段发育的方解石脉体进行了显微结构，C-O-Sr同位素，REY元素以及U-Th定年研究，确定了矿物沉淀与地震事件的关系，发生时间以及流体来源。结合野外调查，讨论了中更新世以来亚东裂谷的伸展速率变化。

1 地质背景

亚东-谷露裂谷作为藏南地区规模最大的一条裂谷，自南向北依次穿过高喜马拉雅、特提斯喜马拉雅、雅鲁藏布江缝合带、拉萨地块，终止于拉萨地块北缘，全长约500km，总体表现为半地堑和地堑(图1)(Armijoetal., 1986；Wuetal., 1998；Leloupetal., 2010)。亚东-谷露裂谷南段即亚东裂谷，位于雅鲁藏布江缝合带与印度板块之间，长约200km，由三个北北东走向的地堑(帕里地堑、涅如堆地堑和热隆地堑)右阶组成(图2)。每个地堑长约45～80km，宽约20km，边界断层为裂谷东侧北西向延伸的西倾正断层(Haetal., 2019b; Wangetal., 2020)。亚东裂谷断层行迹清晰，可以观察到新鲜的地表破裂带，以及山根处的断层三角面，代表裂谷发生了较为活跃的伸展变形，晚第四纪以来亚东裂谷的平均伸展速率为0.9±0.3mm/yr(Haetal., 2019b；Wangetal., 2020)。

亚东裂谷发育在前寒武系变质岩、下古生界、上古生界、中生界和新生界中(图2)。前寒武系为高喜马拉雅结晶岩带亚东岩群。下古生界主要为奥陶系灰岩、粉砂岩。上古生界为陆相沉积，其中石炭系多为碳酸盐台地相沉积，以灰岩为主，二叠系主要为砂岩、砾岩及灰岩(刘文灿等，2014)。中生界在裂谷带分布广泛，下、中三叠统主要为砂质板岩、炭质板岩，上三叠统为变质细砂岩和粉砂质板岩。侏罗系、白垩系以海相沉积为主，主要岩性为砂岩、页岩和硅质岩等。新生代地层主要为碳酸盐和浅海陆源碎屑沉积(万晓樵等，2002)。新近系沉积受地貌、湖泊、断裂影响明显，多见河湖相、洪积和坡积相、冰碛相沉积，均不整合于下覆地层之上(刘文灿等，2006；Liuetal., 2017)。研究区内岩浆岩主要沿断裂出露，基性侵入岩出露较少，中酸性侵入岩分布面积较大，主要由加里东期花岗岩和新生代淡色花岗岩组成(张宏飞等，2005；刘文灿等，2014)。

2 方解石样品采集与测试方法

本文通过卫星照片和野外调查确定断层展布、产状和方解石脉体的发育位置，并采集定向样品。在对方解石脉样品切割、抛光，超声波清洗以去除风化表面后，使用1～2mm微钻获取样品粉末。对于方解石脉，为了避免围岩污染，至少远离围岩或角砾1mm采集样本。对于围岩，同样在远离裂隙和脉体1mm处采集样品。最终将所得的脉和围岩样本用于稀土元素和C-O-Sr同位素分析以及U系列测年。所有抛光样品均制备了显微薄片，以进行光学和阴极发光(CL)显微镜的微观结构分析。

Sr同位素分析在中国科学院地质与地球物理研究所固体同位素地球化学实验室完成。化学前处理及同位素测试流程参考Lietal.(2019)，测试仪器为Thermofisher-Triton-Plus多接收热电离质谱仪。采用国际标准样品NBS987对数据采集期间的仪器稳定性进行了评价，测量平均值为87Sr/86Sr=0.710249±0.000012，与报告值一致(Lietal., 2019)。

稀土元素含量分析在中国科学院地质与地球物理研究所成矿年代学实验室进行。称取40mg样品粉末加入2%HNO3溶解，待溶解完成后加入1mL 500×10-9In内标，并用1%HNO3稀释至50g。测试仪器为ICP-MS(Element)。采用标准曲线法以In内标校准仪器漂移。以国家标准参考物质(花岗岩GSR1；安山岩GSR2；玄武岩GSR3)进行质量监控，绝大多数元素RSD≤10%。

铀系测年在西安交通大学同位素实验室使用MC-ICP-MS(Thermo-Finnigan Neptune-plus)进行。化学前处理方法参照Edwardsetal.(1987)。仪器、标准化、半衰期和测试流程见Chengetal.(2000, 2013)。

3 断层与脉体特征

3.1 宏观特征

方解石脉和围岩样品获取自亚东裂谷的三个采样点，分别为涅如堆地堑北侧达凯研究点和南侧土热勒研究点，以及帕里地堑北部的冲巴研究点(图2)。

达凯剖面位于涅如堆地堑北部东侧达凯村东南2km处，剖面长约350m，围岩为上二叠统板岩，自东向西分别发育断层破碎带以及F1和 F2两条断层(图3a-c)。断层F1倾向292°、倾角69°，F2倾向308°、倾角51°，与区域内亚东裂谷断层产状一致，破碎带内发育一组与断层走向一致倾向相反的共轭破裂(图3c-f)。沿次级破裂和断层面的阶步位置充填方解石脉(图3d-k)，脉体颜色为浅灰色和黄色，厚度0.3～1cm长度小于10cm(图3d-g)。

图3 涅如堆地堑正断层构造特征(达凯研究点)

土热勒剖面位于涅如堆地堑南部东侧，长约300m，围岩为上二叠统灰岩，自东向西分别发育F3和 F4两条断层，断层F3倾向289°、倾角57°，F4倾向291°、倾角64°，卫星影像上两条断层迹线清晰(图4a-c)。靠近盆地的F3断层错断上下盘灰岩，形成典型的断层三角面，断层面底部距离顶部约150m，底部1.5m为白色新鲜面，为最近一次古地震的地震遗迹(图4d-e)。F3和F4断层之间为碎裂岩带，表现为方解石胶结的灰岩角砾或早期方解石角砾，断层两侧发育破碎带(图4f)。该点方解石脉或方解石胶结物颜色为白色和黄色(图4e, f)。

图4 涅如堆地堑正断层构造特征(土热勒研究点)

冲巴剖面位于帕里地堑冲巴芒错北部约2km处，岩性为奥陶系灰岩以及上覆第四系冰川沉积物(图5a)。该点发育F5、F6、F7三条断层(图5b, c)，其中F5为裂谷边界主断层，倾向320°、倾角60°，上盘发育第四系冲积物和钙华，下盘为奥陶系灰岩，沿断层发育约1m宽的灰黑色断层泥及下盘破碎带。F6、F7为次级断层，断层F6倾向240°、倾角66°，F7倾向300°、倾角54°，断层切穿上下盘奥陶系泥晶灰岩，两侧均为断层破碎带，发育断层角砾岩(图5b-g)，沿F6断层面可见沿擦痕充填的白色方解石脉(图5g)，沿F7断层面发育厚层方解石脉，宽度可达5cm，颜色为白色和灰白色(图5g、图6g-j)。

意识的培养不是一蹴而就的，长期潜移默化的影响才能形成稳定的意识形态。面对我国高校创新创业教育滞后和师资不足的现状，光靠几节创新创业课程是不能有效激发师范生创新创业意识的。只有在日常的教学过程中，教师有计划的创新创业意识渗透与灌输才能提高教育成效。

图5 帕里地堑正断层构造特征(冲巴研究点)

图6 断层活动相关的方解石脉体野外宏观构造特征及采样位置

3.2 脉体微观结构

涅如堆地堑北侧达凯研究点(图7a, c-e)和帕里地堑北部的冲巴研究点(图7f)阶步处充填的方解石脉均为块状结构，部分脉体中可见岩石破碎形成的不规则的围岩角砾(图7a)。方解石脉体边界不平直，呈弯曲状(图7c)或锯齿状(图7e)。达凯研究点处方解石脉表现为多期次特征，同一裂隙中充填了粒径不同的多条平行脉体(图7c, e)。粗颗粒方解石粒径可达1mm，发育双晶，细颗粒方解石粒径0.1～0.4mm。方解石阴极发光颜色均匀，晶体没有观察到分带，表明在方解石脉生长过程中没有显著化学成分变化，粗粒的块状脉阴极发光颜色为暗红色，细粒块状方解石阴极发光颜色为红色(图7e)。单层脉体内部可见平直的后期破裂，延伸方向与边界主破裂方向相同(图7d)。

土热勒断层面上的方解石脉同样具有块状晶体特征，发育双晶。方解胶结物中含有围岩或先期方解石角砾，角砾呈棱角状，磨圆差(图7b)。

冲巴研究点沿F6断层阶步位置发育的方解石脉具有块状结构，靠近擦痕的暴露面以及脉体-围岩边界呈弯曲状(图7f)。沿F7断层发育一处厚层(5cm)方解石脉，脉体-围岩边界呈弯曲或锯齿状，方解石晶体表现为晶体拉长、长轴方向垂直于围岩岩壁(图7g)。晶体从脉-围岩界面向脉中心方向生长，同向生长导致脉-围岩界面处存在生长竞争。脉体中间夹杂的围岩条带以及再次出现的方解石晶体竞争生长表明脉体经历了再次开启-愈合过程(图7g)。靠近F7断层的上盘围岩内发育的方解石细脉同样具有多次沉淀的特征，脉体-围岩边界较平直，晶体竞争生长方向表明第二次破裂发生在先期脉-围岩界面处(图7h)。

图7 亚东裂谷与正断活动相关的方解石脉微观构造特征

4 地球化学特征

为确定形成方解石脉的流体来源，本文在达凯、土热勒和冲巴研究点采集了碳氧同位素样品33个，锶同位素样品29个以及稀土元素-钇样品28个。

4.1 碳氧同位素

δ13CV-PDB和δ18OV-SMOW分析结果示于表1和图8a。

表1 亚东裂谷方解石脉、围岩和钙华的氧-碳-锶同位素值

涅如堆地堑达凯研究点方解石脉δ13CV-PDB值为-9.4‰～-4.8‰、平均值为-8.2±0.3‰，δ18OV-SMOW值为-4.4‰～-3.3‰、平均值为-3.9±1.5‰。涅如堆地堑达土热勒研究点方解石脉δ13CV-PDB值为+0.9‰～+2.2‰、平均值为+1.7±0.5‰，δ18OV-SMOW值为-4.2‰～-2‰、平均值为-3.3±0.8‰；围岩δ13CV-PDB值为+3.5‰～+4.5‰、平均值为+4.2±0.3‰，δ18OV-SMOW值为+9.5‰～+17.4‰、平均值为+13.6±2.4‰；钙华δ13CV-PDB值为+5‰～+5.9‰、平均值为+5.5±0.5‰，δ18OV-SMOW值为+13.9‰～+17.4‰、平均值为+15.6±1.8‰。帕里地堑冲巴研究点方解石脉δ13CV-PDB值为+1.6‰～+3.5‰、平均值为+2.8±0.7‰，δ18OV-SMOW值为-2.6‰～+0.3‰、平均值为-0.8±1.0‰；围岩δ13CV-PDB值为+0.8‰～+3.2‰、平均值为+2±1.2‰，δ18OV-SMOW值为+0.7‰～+2.4‰、平均值为+1.6±0.9‰。

藏南地区报道的奥陶系碳酸盐岩δ13CV-PDB为 0.17‰～2.32‰、δ18OV-PDB为-15.33‰～-11.45‰(Yuetal., 2021)， 二叠-三叠系碳酸盐岩δ13CV-PDB为0至4‰、δ18OV-PDB为-14‰至-12‰(Shenetal., 2010)。达凯研究点方解石脉与区域内碳酸盐岩相比具有较低的δ13CV-PDB值，其他两处采样点方解石脉δ13CV-PDB值与区域值相符(图8a)。

图8 方解石脉、围岩和钙华的同位素特征

各研究点方解石脉的δ18OV-PDB值均轻于围岩以及藏南地区碳酸盐岩，值得注意的是，达凯研究点方解石脉具有极低的δ18OV-PDB值。方解石脉中的δ18OV-PDB特征取决于沉淀流体的来源和沉淀温度(Faure, 1998)。假设大气水在全新世没有显着变化，并且代表了过去的地下水，参考海拔 3000m 以上的高喜马拉雅河流已发布的δ18OV-SMOW值(-16‰～-19‰)(Gajurel al., 2006；Li and Garzione, 2017)，使用碳酸盐-流体古温度计 (O’Neiletal., 1969)，可以计算出，如果达凯研究点方解石脉来源为大气水，则沉淀温度在122～252℃之间(图8c)，若来源为大气水以及相邻区域碳酸盐岩，则需要更高的沉淀温度。

4.2 87Sr/86Sr同位素和REY元素模式

涅如堆地堑达凯和土热勒研究点方解石脉87Sr/86Sr同位素比值介于0.7224～0.7340之间(图8b, 表1)。明显高于土热勒灰岩围岩的87Sr/86Sr同位素比值(0.7139～0.7189)。帕里地堑冲巴研究点方解石脉87Sr/86Sr同位素比值介于0.7132～0.7159之间。与奥陶系灰岩围岩87Sr/86Sr同位素(0.7112～0.7148)一致。

方解石脉和围岩样品的稀土元素(包括钇)含量见表2。球粒陨石标准化的REY图中显示，达凯研究点方解石脉表现为富LREE和Eu正异常，土热勒研究点方解石脉表现为富LREE，Eu负异常以及Y正异常，冲巴研究点表现为富LREE，Eu负异常(图9)。

表2 亚东裂谷方解石脉和围岩REY元素分析结果(×10-6)

图9 达凯研究点(a)、土热勒研究点(b)和冲巴研究点(c)方解脉和围岩的球粒陨石标准化REY模式(标准化值据Haskin et al., 1968)

5 U系年代学

我们对达凯研究点3个擦痕方解石脉，土热勒研究点2个断层方解石脉，和冲巴研究点3个断层方解石脉进行U-Th年代学测定(表3)。

表3 亚东裂谷断层方解石脉U系测年结果

达凯研究点发育擦痕的方解石脉(YD112、YD132和YD164样品)均具有块状晶体结构，测得其年龄超过测试上限(大于600ka)。土热勒研究点F3方解石脉(NFA0样品)具有块状结构，测得年龄为679±129ka，F4方解石脉230Th/232Th值过低(4×10-6～5×10-6)，测得年龄超过测试上限(大于600ka)。冲巴研究点F6断层方解石脉YD410样品表面发与擦痕，且具有块状内部结构，测得年龄为544±26ka，F7断层方解石脉具有拉长的块状结构，在YD314和YD348样品中(图6i, j)测得年龄为408±13ka和444±15ka。微观结构观察显示这些方解石脉未重结晶，表明以上年龄是可靠的。

6 讨论

6.1 方解石脉结构与成因

断层作用各时期形成的方解石脉类型包括：(1)水力压裂产生的同震结晶：晶体形态为块状，非自形等粒结构，形成于流体压力快速下降或流体突然停止上升导致的高矿物饱和度(图10a)(Phillips, 1972；Bons，2001；Boullieretal., 2004；Uysaletal., 2007)。(2)断层拉张部位的同构造结晶：晶体形态为拉长块状和扩展状(Bonsetal., 2012)。拉长块状脉晶体具有明显的竞争生长特征，平均晶粒宽度沿生长方向增加，通常沿围岩表面向裂缝内生长(图10b)(Durney and Ramsay, 1973；Oliver and Bons, 2001)。扩展状脉具有横跨整个脉宽的晶体，是通过对在先存晶粒的不同位置多次破裂和愈合而形成的(Durney and Ramsay, 1973；Ramsay, 1980)。(3)震间期沉淀，如纤维脉(图10c)(Hilgersetal., 2001；Nolletetal., 2005)和暴露面的薄层泥晶涂层(Flottéetal., 2001；Nurieletal., 2012)。

图10 方解石脉内部特征结构(据Woodcock et al., 2007)

本次研究中亚东裂谷涅如堆段发育的方解石胶结物和擦痕脉的微观结构具有等粒块状结构(图7c, d)，且具有不平直的脉体边界，代表了与地震断层滑动相关的拉张部位的方解石快速沉淀(Gratier and Gamond，1990；Nurieletal., 2019)。亚东裂谷帕里段一处沿断层发育的方解石脉具有拉长的块状结构(图7e, f)代表同构造时期的裂缝愈合。没有观察到与溶解沉淀相关的裂缝缓慢愈合的证据，如纤维脉和泥质涂层(Gratier and Gamond，1990；Hilgersetal., 2001；Nurieletal., 2012; Zhangetal., 2015)。CL观察表明，方解石擦痕样品CL颜色均匀(图7d)，表明在单个沉淀事件期间结晶。由于每个破裂-愈合事件都意味着裂缝打开，块状或拉长块状晶体从充满流体的开放裂缝中析出(Hilgersetal., 2001)，多条平行发育的方解石脉代表多次裂缝密封事件(图7g, h、图11)。本文中样品观测显示，不同期次的层状脉被围岩包体条带(图7e, f)分隔，或被锋利的晶体间破裂切割(图7c, d)，表明新的破裂事件均在早期裂缝完全密封后才发生。

图11 亚东裂谷地震相关方解石脉生长模式

6.2 流体来源

根据图8a,δ13CV-PDB和δ18OV-SMOW两者具有比较明显的线性关系，指示流体具有二元混合特征。本文的三个研究点中，达凯点具有异常轻的碳氧同位素特征。碳同位素组成表明，方解石中的碳可能来自地球深处(即地幔或下地壳)，也可能来自碳酸盐或有机质。然而，δ13C值较集中，排除了有机来源的可能性，位于幔源碳(-5.5±0.5‰)(Sakai and Matsuhisa，1996)和区域内碳酸盐岩之间，表明具有幔源贡献。极低的δ18O值反映了大气水来源，以及高于120℃的沉淀温度。方解石脉的高Sr同位素比值超过区域内奥陶系，二叠-三叠系灰岩的同位素值。因此排除了方解石仅来自溶解了碳酸钙的大气水。高Sr同位素比值表明，在流体的运移过程中通过了放射成因的基岩，如高Sr比值的硅质碎屑和花岗岩(McNutt, 2000)。因此区域内喜马拉雅淡色花岗岩体(Sr比值0.756～0.775)(杨晓松和金振民，2001；李开玉等，2020)以及二叠系板岩围岩为流体提供了高Sr比值特征。深层渗透的大气水与变质/花岗质基岩的相互作用很可能发生在温度高于122～252℃ 的条件下，Eu正异常同样指示了热液环境。因此该点的方解石脉形成机制可能为幔源脱气与深层大气水混合成因(图12)。外部来源的流体反映了脉体形成于开放系统，即断层活动时的裂隙连通状态，这与块状结构方解石脉代表的同震时期相匹配。

图12 地震相关的断裂带流体循环模式图

土热勒和冲巴研究点方解石脉具有相似的地球化学特征，δ13CV-PDB值与区域碳酸盐岩一致(Shenetal., 2010；Yuetal., 2021)，δ18OV-SMOW值低于碳酸盐围岩，表明具有大气水贡献。Ce和Eu负异常以及Y的正异常反映了低温氧化环境(Bauetal., 1997；Bolharetal., 2004)。考虑到三个研究点所处位置邻近，脉体记录的温度特征能够对应相对埋深，与达凯研究点方解石具有的高温条件相比，土热勒和冲巴研究点方解石形成于更浅的区域。分析的样品还显示了与灰岩围岩密切相关的方解石REY特征和87Sr/86Sr值，指示了流体来自大气水以及溶解的碳酸盐岩(图12)。

6.3 方解石脉与地震活动

地震引起岩石破裂后，流体发生突然的压力降低和平衡浓度改变(Putnisetal., 1995；Robertetal., 1995；Barnes and Rose，1998；Putnis and Mauthe, 2001)，因此形成充填于裂缝或角砾间的矿物沉淀。受到矿物过饱和程度，裂隙张开宽度，以及先存晶体取向的不同，沉淀的矿物晶体具有块状或拉长块状等结构(Oliver and Bons, 2001；Nolletetal., 2005；Woodcocketal., 2007)。活动断层的实验和观察表明，在愈合作用影响下，同震期间形成的张裂隙将在此后很快(几天到几个月)关闭(Hiramatsuetal., 2005; Lietal., 2006)。本次研究中沿断层阶步部位发育的方解石脉均具有块状结构，代表同震时期方解石沉淀，部分样品发育由平直破裂面分隔的多层脉体，代表了流体多次注入。土热勒研究点沿断层面发育的方解石脉体内部发育棱角状围岩角砾，代表断层滑动期间岩石破裂伴随的流体注入。冲巴研究点沿断层面发育一套厚层方解石脉，具有拉长的块状结构，代表震后同构造期的破裂愈合。因此本文中选取的方解石脉能够代表同震或同构造的断层活动。即研究区自北向南，达凯断层记录了超过两期的地震活动，活动时间距今均超过600ka；土热勒研究点F3断层在距今679±129ka发生一期地震活动；冲巴研究点记录了NNE向断层544±26ka的地震活动以及，NE向F7断层408±13ka和444±15ka的两期地震活动。在以上三个研究点中，涅如堆地堑南段土热勒点发育的断层三角面较为完整，新鲜断层面与三角面顶部高差约150m，断层倾角为57°。若本研究中方解石脉记录的地震年龄为初始活动时间，则裂谷该段中更新世以来的E-W向伸展速率为0.17～0.24mm/yr，由于方解石脉的沉淀晚于最早期的断层活动，因此断层实际活动时期更长，因此该速率为最大值，即中更新世以来的伸展速率小于晚更新世以来0.8±0.3mm/yr的伸展速率(Wangetal., 2020),以及GPS观测的现今伸展速率6mm/yr(哈广浩，2019)。即中更新世以来，亚东裂谷涅如堆段断层活动速率增大。

7 结论

本文通过对藏南亚东裂谷地区断裂带中的方解石脉体成因、地球化学以及年代学分析，得到以下认识：

(1)亚东裂谷涅如堆段和帕里段沿裂谷断层发育方解石脉。三个采样点中沿断层面阶步产出的方解石脉具有块状结构，土热勒点沿断层面发育的方解石脉体中包含角砾状围岩，以上两种结构代表了流体在同震阶段的快速沉淀；冲巴点沿断层发育一处厚层方解石脉，具有拉长的块状结构，代表震后同构造阶段的流体沉淀。

(2)涅如堆地堑北部方解石脉REY模式反映了热液特征，异常轻的δ13CV-PDB值具有幔源特征，而极低的δ18OV-SMOW值显示出大气水贡献。涅如堆地堑南部和帕里地堑北部方解石脉REY模式反映了浅成的低温特征，δ13CV-PDB和87Sr/86Sr值与区域内灰岩相近，代表了围岩对方解脉的贡献，偏负的δ18OV-SMOW表明方解石脉同样存在大气水贡献。

(3)地球化学数据显示，达凯点方解石脉显示了来自大气水和地幔流体贡献，沉淀温度约为122～252℃，土热勒和冲巴点的流体来源主要为灰岩围岩和大气水，方解石在低温条件沉淀。地震相关的挤压和破裂导致地幔流体延断层上升，在深部区域与下渗或深埋的大气水共同导致方解石沉淀。在浅部区域，由于地震破裂产生的压力和温度变化，灰岩中过饱和流体和浅部地下水混合产生方解石沉淀。

(4)方解石脉的显微结构和地球化学数据显示，脉体形成与地震活动相关，其中裂谷的涅如堆段记录了一期679±129ka的地震活动，且该段断层活动速率从中更新世明显增大；帕里段北部记录了544±26ka、444±15ka和408±13ka三期地震活动。

致谢铀系测年工作得到西安交通大学同位素实验室宁有丰研究员和贾雪雪的悉心指导，碳氧同位素数据分析得到了中国地质科学院矿产资源研究所高建飞副研究员的宝贵建议和指导。刘俊来教授、张波副教授和本刊编辑提出了宝贵意见和建议。在此一并表示感谢。