PNAS3He4He研究揭示印度板块在青藏高原下方的俯冲模式

　　自50Ma以来印度与欧亚的陆陆碰撞和持续汇聚，造就了世界最高峰喜马拉雅山脉和世界屋脊青藏高原，这里是研究陆陆碰撞绝佳的天然实验室。在陆陆碰撞中，印度岩石圈在青藏高原下方是如何俯冲的一直存在争议。目前两种典型的模型有两种：第一种是印度岩石圈地幔以平板俯冲方式直接底垫在青藏高原地壳之下（McKenzieetal。，2019）；第二种是印度岩石圈板片在欧亚陆陆碰撞边界向高原下方地幔中俯冲，角度由西向东变陡（Lietal。，2008）。这两种模式的本质区别是青藏高原地壳与俯冲的印度岩石圈地幔之间是否存在热的软流圈地幔（图1）。
　　图1大陆碰撞的几何模式图和位置图（Klempereretal。，2022）。（A、B）定义平板俯冲、俯冲和地幔缝合线（印度岩石圈与（A）青藏岩石圈或（B）软流圈的边界）的卡通图。蓝色：地幔岩石圈（冷）；红色：软流圈（热）；（C）研究区域构造简图。黄色条带：图4中的横截面；红色条带：体波层析成像给出的地幔缝合线；绿色条带：接收函数显示的印度地壳前缘；蓝色条带：全波形成像给出的240km厚印度岩石圈的北部边界；细灰线：缝合线主要断裂带：MFT（主前缘断层）；STD（藏南拆离系）；YZS（雅鲁藏布缝合线）；BNS（班公怒江缝合线）；JRS（金沙江缝合线）；主要走滑断层ATF（阿尔金断裂）、KF（昆仑断裂）和KKF（喀拉昆仑断裂）。粗黑线为活动裂谷：YR（亚热裂谷）；TG（塔口拉裂谷）；LG（隆格尔裂谷）；TY（当惹雍错裂谷）；PX（朋曲申扎裂谷）；YG（亚东谷露裂谷）；CS（错那裂谷）
　　Klempereretal。（2022）通过多年的努力，测量了青藏高原上225个温泉的3He4He，使用该比值分析了青藏高原下方是否存在热地幔楔，由此限定了印度欧亚板块之间的地幔缝合线的位置，进而揭示了印欧陆陆碰撞的深部过程，研究成果于2022年3月发表在PNAS上。
　　3He4He比值分析方法基于地壳流体中存在不同成分的惰性气体，通过测量不同惰性气体在地壳流体的含量，解析流体的来源和运移过程。4He是由地壳中U、Th等放射性元素衰变产生的；而3He来自地幔，主要通过地幔熔融或交代流体流动释放，并通过断层、地壳流扩散以及岩浆侵入等方式运移到地表，主要发生在大陆伸展或岩浆活动区域内（Ballentineetal。，2002）。印度岩石圈是古老而冷的克拉通，4He丰度较高，但3He无法释放至地表。而位于印度俯冲板片上方的地幔楔和青藏高原Sn不发育区的上地幔都是热的，并且容易发生交代作用，因此这些区域的热泉中富集3He。
　　Klempereretal。（2022）对测量的3He4He进行统计分析，RA表示空气中的3He4He，R表示温泉样本中的3He4He，RC表示对R进行大气贡献校正。标定古老而寒冷的印度岩石圈RCRA0。06（地壳样本），初始熔融且含CO2和H2O的热青藏高原地幔RCRA0。1（地幔样本），地幔缝合线在地表的位置则为两者之间的过渡（图2）。结果给出了亚欧板块之间地幔缝合线位置（图3）。该边界与YZS近乎平行，东西跨越青藏高原1000多公里，边界以南分布中源地震（65100km）、上地幔顶部地震波速度高，以北分布第四纪火山，上地幔顶部地震波速度低（图3），地幔熔融或有流体流动。
　　图2青藏高原块体同位素数据（Klempereretal。，2022）。纵轴为RRA，横轴为（HeNe）sample（HeNe）ASW，FieldT。（HeNe）sample（HeNe）ASW，FieldT10竖线分隔本研究中的高质量数据（圆圈）和不确定数据（六边形）；三角形表示其他作者的高质量数据。采样位置：塔里木、柴达木盆地（蓝色）；羌塘块体（橙色）；拉萨块体（白色）；特提斯喜马拉雅（灰色）；大、小、亚喜马拉雅（黑色）。十字表示水样。、中间，地幔标定分别对应0。5和1的地幔含量。红色、黄色、蓝色箭头和字母，以及绿色字母与图三中样本的颜色对应
　　图3地幔样本和样本分布图（Klempereretal。，2022）。（A）3He4He比值（RCRA）分布。白色细实线：首选的10km宽的氦边界；白色透明虚线：氦边界的两种可能；白色小方块：莫霍面错断和上地幔速度突变的位置（连接方块的线）；黑线：裂谷；蓝色绿色红色地球物理边界如图1C所示；（B）RCRA是距离首选氦边界距离的函数（使用A中的北侧氦边界虚线，并根据地体亲和度进行颜色编码，如图1C和图2所示：蓝色（印度样本）；黄色（地幔样本）；绿色（过渡样本））；（C）测量的〔3He〕和〔4He〕的丰度如主地图中的彩色编码，显示出明显的双峰性。（D）地震体波层析成像结果在100km深度的P波速度扰动。（E）伴随层析成像（全波形拟合）结果在80km深度的S波速度扰动
　　与青藏高原区域的多种地球物理研究成果进行对比分析，Klempereretal。（2022）对地幔缝合线的研究在一定程度上解决了地球物理研究结果中的部分争议，地幔缝合线为研究印度板片的俯冲撕裂提供了新证据。在青藏高原南部分布有多条近南北走向的裂谷，裂谷的形成与印度板片撕裂、软流圈上涌有着密切的关系（Chenetal。，2015；Liuetal。，2020）。地幔缝合线在裂谷处发生分段，裂谷中富集的地幔样本与裂谷下方的莫霍偏移以及板片撕裂一致，意味着印度岩石圈在俯冲同时存在板片撕裂（图3）。比较特殊的是在CS裂谷，地幔缝合线相比地表缝合线YZS向南还偏离了100km，说明印度岩石圈地幔从地壳拆离（Shietal。，2015）。
　　地幔缝合线只在局部（85E附近）与印度地壳前缘重合，而在高原西部80E以西地幔缝合线向北西侧延伸至阿尔金断裂，标志印度岩石圈俯冲的北部边界。而在85E以东，地幔缝合线与印度地壳前缘逐渐分离，在92E至93E向南撤退了200km至300km（图3），标志印度俯冲板片在地幔缝合线处从印度地壳剥离并向深部俯冲，从西至东俯冲角度越来越陡。结合体波走时速度结构成像结果（Lietal。，2008）、全波形拟合成像结果（Chenetal。，2017）与剪切波传播特征（BarronandPriestley，2009）等，Klempereretal。（2022）认为青藏高原下方岩石圈地幔虽然已经被印度岩石圈地幔所取代，但青藏地壳和印度岩石圈地幔之间仍存在一层薄的、热的软流圈地幔（图4）。
　　Klempereretal。（2022）对印欧板块之间的地幔缝合线的定位，否定了印度岩石圈直接底垫在青藏高原地壳下方的平板俯冲模式，给出了印度冷的岩石圈地幔与热的或流动的软流圈地幔在高原下方叠置的模式（图4），同时有利于直接的研究大陆俯冲板片撕裂。这种陆陆碰撞模式与现代纳斯卡海脊的平板俯冲有一定的几何相似性（Antonijevicetal。，2015）。
　　图4青藏高原横纵等比例的横截面（沿图1C中的黄色条带）（Klempereretal。，2022）。所有实灰线都来自已发表的地球物理结果。地幔缝合线位于在地壳氦和地幔氦样本之间的边界正下方。黑色箭头显示印度地壳（绿色）和岩石圈（蓝色）相对于亚洲板块（灰色）的相对运动。俯冲板块的倾角不受氦同位素比值研究的限制，可能沿走向变化，如图所示，根据地震衰减和层析成像研究，倾角为5至30。印度下地壳的前缘位置与地幔缝合线的位置沿东西向变化，由东向西，分别在地幔缝合线的北部、与地幔缝合线重合、或可能在其南侧。MHT：主喜马拉雅逆冲断层
　　主要参考文献（上下滑动查看）
　　AntonijevicSK，WagnerLS，KumarA，etal。Theroleofridgesintheformationandlongevityofflatslabs〔J〕。Nature，2015，524（7564）：212215。
　　BallentineCJ，BurgessR，MartyB。Tracingfluidorigin，transportandinteractioninthecrust〔J〕。ReviewsinMineralogyandGeochemistry，2002，47。
　　BarronJ，PriestleyK。ObservationsoffrequencydependentSnpropagationinNorthernTibet〔J〕。GeophysicalJournalInternational，2009，179（1）：475488。
　　BarronJ，PriestleyK。ObservationsoffrequencydependentSnpropagationinNorthernTibet〔J〕。GeophysicalJournalInternational，2009，179（1）：475488。
　　ChenM，NiuF，TrompJ，etal。LithosphericfounderingandunderthrustingimagedbeneathTibet〔J〕。NatureCommunications，2017，8（1）：15659。。
　　ChenY，LiW，YuanX，etal。TearingoftheIndianlithosphericslabbeneathsouthernTibetrevealedbySKSwavesplittingmeasurements〔J〕。EarthandPlanetaryScienceLetters，2015，413：1324。
　　KlempererSL，ZhaoP，WhyteCJ，etal。LimitedunderthrustingofIndiabelowTibet：3He4Heanalysisofthermalspringslocatesthemantlesutureincontinentalcollision〔J〕。ProceedingsoftheNationalAcademyofSciences，2022，119（12）：e2113877119。
　　LiC，VanderHilstRD，MeltzerAS，etal。SubductionoftheIndianlithospherebeneaththeTibetanPlateauandBurma〔J〕。EarthandPlanetaryScienceLetters，2008，274（12）：157168。
　　LiuZ，TianX，YuanX，etal。ComplexstructureofuppermantlebeneaththeYadongGuluriftinTibetrevealedbyStoPconvertedwaves〔J〕。EarthandPlanetaryScienceLetters，2020，531：115954。
　　McKenzieD，JacksonJ，PriestleyK。Continentalcollisionsandtheoriginofsubcrustalcontinentalearthquakes〔J〕。CanadianJournalofEarthSciences，2019，56（11）：11011118。
　　ShiD，WuZ，KlempererSL，etal。Receiverfunctionimagingofcrustalsuture，steepsubduction，andmantlewedgeintheeasternIndiaTibetcontinentalcollisionzone〔J〕。EarthandPlanetaryScienceLetters，2015，414：615。
　　撰稿：李聃聃，田小波岩石圈室
　　美编：陈菲菲