从全球一维地震学模型看地球的下地幔,绝大部分似乎是“平淡无奇”的区域。近年来,随着地震观测与反演技术的进步,人们逐渐“看”到下地幔其实有着不同寻常的精彩。其中最吸引人的,是人们通过不同观测均发现核幔边界以上存在厚达上千公里的结构体,它们集中分布于非洲和太平洋以下,横向分布可达数千公里,这些结构体有较为清晰的边界,呈现与周围地幔相比偏低的横波波速(因而被称为大型横波低速带,Large Low Shear Velocity Provinces,简称为LLSVPs),并且具有较高的横波纵波波速比。显然,研究LLSVPs的特征及成因对人们理解地幔对流机制具有十分重要的意义。
什么原因造成LLSVPs的横波波速异常?其主因究竟是温度异常还是由于物质成分异常,以及进一步地造成横波波速异常的成分又来自于什么地方?人们对此有很多解读,也存在很大争议。近日,来自英国伦敦大学学院的Thomson等在Nature上发表了题为“Seismic velocities of CaSiO3 perovskite can explain LLSVPs in Earth’s lower mantle”的最新成果,为解释LLSVPs提供了新的约束。
Thomson等的工作源于矿物学的“传统”思路:利用热力学模型计算矿物组合的波速,并将其与实际地震观测数据相比较,从中推算出造成地震波速异常的矿物学成因。如图1所示,用现有的下地幔矿物数据库计算出的玄武岩波速(图1中红色虚线)始终明显大于平均地幔(图1中的地幔岩pyrolite模型)。他们注意到,该数据库中默认CaSiO3钙钛矿(以下简称为CaPv)为立方晶体,并且参数主要来自于对第一原理模拟数据的拟合,但最近的高温高压实验表明,立方结构CaPv的实际横波波速也许要远低于理论模拟结果(Gréaux et al., 2019;参见前沿报道《立方结构CaSiO3钙钛矿的声速原位测量及对下地幔顶部物质成分的启示》)。由于CaPv在俯冲带玄武岩矿物组合中的占比可达20%以上(如图1),因此现有热力学模型显然高估了洋壳深俯冲物质的横波波速。
图1 利用热力学模型计算出的玄武岩模型在地幔条件下的矿物组合与横波波速(Stixrude and Lithgow-Bertelloni, 2012)
和Gréaux等的工作类似,Thomson等利用高温高压实验方法原位测定了CaPv随温压变化的结构相变及其波速,不同的是他们在样品中还加入了少量Ti,并且特别注意了高温下样品的温度梯度控制。如图1所示,Thomson等发现他们在常温常压下测得的CaPv波速与Gréaux等的结果吻合得很好,但是在高温下存在较大偏差,他们把这些偏差归咎于Gréaux等的实验在高温下可能存在较大的温度梯度(样品实际温度可能要低于热电偶测量值约250℃-300℃)。与此同时,新的实验发现Ti的加入将CaPv从四方相转变到立方的温度提高了近800 K,而高压相变会带来波速的较大异常变化(纵波减小4%-14%,横波减小8%-20%,如图3所示)。
图2 Thomson等测得的CaPv横波波速与前人研究结果的对比。绿色粗实线为沿1500K地幔绝热曲线的立方CaPv的横波波速
图3 CaPv样品在高温高压条件下的波速变化
有了新的实验数据,再精选前人的相关实验数据,Thomson等重新拟合了立方和四方结构CaPv的热力学状态方程参数,然后重新计算了各种矿物学模型的波速。如图4所示,沿着1500 K的地幔绝热线MORB的波速要明显低于平均地幔,考虑到MORB的密度要大于平均地幔,因此洋壳俯冲至下地幔堆积成为LLSVPs的合理解释。这解决了前人热力学模型预测MORB波速高于平均地幔的矛盾(如图1所示)。Thomson等进一步提出,如果CaPv已经相转变为立方结构,则100 GPa下64% MORB+36%harzburgite模型或者125 GPa下48% MORB+52%harzburgite模型可以较好地解释LLSVPs的横波异常以及横波纵波波速比偏大等地震学特征。与此同时,联系图3所反映的CaPv相变带来的波速异常,Thomson等提出1000公里附近的地幔中部的散射体和异常也许也和CaPv密切相关。
图4 MORB模型的波速与平均地幔波速的对比。图中红色为沿着1500 K地幔绝热线的波速,而蓝色为沿着1000 K俯冲带绝热线的波速
回到CaPv矿物本身,Thomson等的工作是不是就宣告了其弹性波速研究的胜利和终结呢?答案显然是否定的:首先,正如Thomson等在其文中述及,他们的实验温压区间还非常有限(大约12GPa、300-1500 K),从如此有限的数据回归热力学模型推测整个地幔深处(与LLSVPs相关的温压可高达130GPa、3000 K)显然有很大的不确定性;其次,不同实验之间有差别,Thomson等虽然做了分析和推断,但评价数据的可靠性显然还有待更多研究的交叉验证;再有,什么原因造成第一原理模拟得到的波速要远高于新的实验结果?这仍然是困扰高温高压领域的热点和难点,同时也是更加深入理解CaPv这类具有明显非谐效应矿物的关键。
主要参考文献
Gréaux S, Irifune T, Higo Y, et al. Sound velocity of CaSiO3 perovskite suggests the presence of basaltic crust in the Earth’s lower mantle[J]. Nature, 2019, 565:218-221.(链接)
Stixrude and Lithgow-Bertelloni. Geophysics of chemical heterogeneity in the mantle[J]. Annual Review of Earth and Planetary Sciences, 2012, 40: 569-595.(链接)
Thomson A R, et al. Seismic velocities of CaSiO3 perovskite can explain LLSVPs in Earth’s lower mantle[J]. Nature, 2019, 572(7771): 643-647.(链接)
(撰稿:张志明、张志刚/地星室)
删除或更新信息,请邮件至freekaoyan#163.com(#换成@)
【前沿报道】Nature:CaSiO3钙钛矿的地震波速实验测定及其对下地幔LLSVPs的启示
本站小编 Free考研/2020-05-30
相关话题/实验 结构
【前沿报道】Nature: Mohns超慢速扩张脊的深部电性结构成像
洋脊是地球系统中火山和地震活动极其频繁、岩浆大规模上涌及新生洋壳形成的巨型活动构造带,是大陆漂移和海底扩张的中心,关乎板块构造学说的理论根基,对于研究板块动力、熔体形成和地幔对流等基本地球科学问题具有十分重要的意义。根据扩张速率大小,洋脊可划分为快速、中速、慢速和超慢速扩张脊等四种类型(图1a)。其 ...中科院地质与地球物理研究所 本站小编 Free考研 2020-05-30【前沿论坛】钮凤林:全波地震层析成像研究东亚上地幔结构
报告人:钮凤林∣编写:蒋梦凡,塔力哈尔·哈帕尔,孙伟家(地球与行星室) 深俯冲岩石圈动力学研究近年来取得很大进展,但完整理解该俯冲过程依然存在困难。地震层析成像显示,一些板块俯冲并停滞在地幔过渡带[1],一些板块则俯冲至下地幔深处[2]。不同俯冲带可表现出迥异的俯冲过程,如海沟后撤或前进;俯冲带的 ...中科院地质与地球物理研究所 本站小编 Free考研 2020-05-30【前沿报道】Nature:立方结构CaSiO3钙钛矿的声速原位测量及对下地幔顶部物质成分的启示
CaSiO3钙钛矿(以下简称CaPv)是地幔转换带(MTR,410-660km)和下地幔(660-2900 km)中的重要矿物组分。前人的计算模拟与高温高压实验均表明,CaPv在高压下会随温度的升高而从四方结构转变至立方结构,但对于该相变的温压边界以及立方结构CaPv的弹性波速还存在很大争议。从高温 ...中科院地质与地球物理研究所 本站小编 Free考研 2020-05-30【前沿论坛】Windley:太古代板块结构——来自热液通道、OPS、异剥钙榴岩和块状硫化物沉积的启示
报告人:Brian F. Windley |整编:李咏晨、陈振宇(岩石圈室) 板块构造理论是20世纪地质学最重要的突破,1912年德国科学家魏格纳提出大陆漂移假说,约半个世纪后美国海洋地质学家赫斯和海洋地球物理学家迪茨几乎同时提出海底扩张说, 1965年加拿大地球物理学家威尔逊首次使用“板块”一词 ...中科院地质与地球物理研究所 本站小编 Free考研 2020-05-30郑天愉等-NC:印度大陆俯冲到缅甸下方的结构证据
印度大陆俯冲(subduction)或底冲(underthrusting)是解释印度与亚洲碰撞以来青藏高原新生代地壳变形、岩浆活动和隆升的一个强有力的构造模式。在喜马拉雅地区,地震成像揭示了包括中下地壳的印度板块向青藏高原之下的底冲构造(Nábelek et al., 2009)。在印度-欧亚大陆碰 ...中科院地质与地球物理研究所 本站小编 Free考研 2020-05-30臧明东等-NHESS:考虑结构面粗糙度影响的区域地震滑坡预测模型
地震滑坡是由地震触发的边坡失稳事件,常常造成巨大的生命和财产损失。近年来,全球强震频发,关于地震滑坡的研究已经成为国际上的热点。我国是一个多山多地震的国家,长期面临着极高的地震滑坡风险,准确预测地震滑坡的空间分布是一项极富挑战性的研究,对于地震防灾减灾具有非常重要的意义。 大量野外调查发现,自然界 ...中科院地质与地球物理研究所 本站小编 Free考研 2020-05-30吕彦等-JGR:西北太平洋地区复杂的上地幔结构和变形
西北太平洋地区在新生代期间经历了广泛的构造变形,是研究板块相互作用动力学过程的重要区域。太平洋板块和菲律宾海板块在鄂霍次克板块和欧亚板块下俯冲,形成边缘海、扩张带和广泛分布的岛弧。然而,由于数据和方法的限制,该区域已有的层析成像结果仍然缺乏共识,了解西北太平洋地区的地球动力学特征是一项具有挑战性的课 ...中科院地质与地球物理研究所 本站小编 Free考研 2020-05-30段武辉等-EPSL:石笋记录的华北倒二冰消期的时间和结构特征
20世纪50年代,Emiliani首次将δ18O分析方法应用到海洋沉积物中,并根据海洋沉积物有孔虫δ18O记录率先提出第四纪存在多次冰期旋回。随后,越来越多的深海岩芯、黄土沉积以及冰芯等古气候记录证实了这种“锯齿状”冰期-间冰期旋回:从间冰期向冰期转换的过程是缓慢的,而从冰期向间冰期转换的过程却是迅 ...中科院地质与地球物理研究所 本站小编 Free考研 2020-05-30范兴利等-JGR:长白山火山及其邻区的高精度三维剪切波速度结构
长白山是一座规模巨大的板内层状复式休眠活火山,自形成以来曾发生多次剧烈喷发,发生在公元946年的“千年大喷发”被认为是有人类历史记载以来的最大的火山喷发事件之一,其飘散的火山灰在日本北海道地区沉降了数厘米厚,甚至在北极格林兰的冰芯中都有记录。此后,长白山又在1413年、1668年、1702年和190 ...中科院地质与地球物理研究所 本站小编 Free考研 2020-05-30宋鹏汉等-JGR:秦岭造山带中下地壳存在北东—南西向低速结构
秦岭造山带位于华北板块与华南板块之间,在中国的中央造山带中心,西连昆仑、祁连造山带,东连大别苏鲁造山带 (图1)。秦岭形成于三叠纪华北板块、华南板块和秦岭微板块之间的汇聚、俯冲和碰撞过程。此后,太平洋板块西向俯冲,秦岭的构造活动在中生代由东向西被重新激活。新生代,秦岭被部分卷入青藏高原东北缘的生长和 ...中科院地质与地球物理研究所 本站小编 Free考研 2020-05-30