近日，中国科学技术大学地球和空间科学学院吴忠庆教授课题组与南京大学阮友谊副教授及中科院精密测量科学与技术创新研究院倪四道研究员合作，在约束地球下地幔的物质组成与温度分布方面取得重大突破。相关研究成果以“Compositional and thermal state of the lower mantle from joint 3D inversion with seismic tomography and mineral elasticity”为题，发表在国际综合学术期刊Proceedings of the National Academy of Sciences (PNAS) 上。

  地球内部可大致分为地壳、上地幔、下地幔和地核。其中，地球下地幔位于地表以下660~2890 km深度处，其体积和质量分别占地球整体的55%和52%，是地球最重要的圈层之一。地震学研究发现，下地幔并不均一，存在分布广泛、尺度不一的波速异常。尤其是位于非洲和太平洋下方的大型低剪切波速省（LLSVPs），其横向和径向尺度可达上千公里。这些异常体的性质、成因以及对地球演化的影响还不够清楚。因此，获取地球下地幔物质成分和温度的空间分布对了解地球的形成、演化与动力学至关重要。

  联合地震学层析成像和地球内部矿物的弹性性质是获得地幔物质成分和温度空间分布的关键手段。由于下地幔极高的温度和压强条件，通过实验测量处于下地幔条件下的矿物弹性性质极具挑战，而常规的矿物弹性第一性原理计算方法及其昂贵。为此，吴忠庆教授课题组提出并发展了计算量不到常规方法十分之一的弹性第一性原理计算方法。近年来，该课题组利用这一方法系统地研究了下地幔主要矿物的弹性性质，结果与相对低温和低压的实验数据很符合。

  该研究利用计算得到的下地幔主要矿物高温高压弹性数据，结合下地幔三维层析成像模型，采用马尔科夫链蒙特卡洛方法反演得到了整个下地幔的三维矿物成分与温度分布（图1），并进一步得到下地幔三维密度模型。

图1.不同深度下矿物成分、铁含量以及温度的空间分布。

  反演结果表明，下地幔的横向温度分布满足高斯分布，其横向展布在1600公里深度以内变化很小，然后随着深度的增加不断变宽。在下地幔最底部，其横向温度分布不再满足高斯分布（图2）。这表明下地幔底部的横向不均一性很强，可能和LLSVPs的存在有关。

图2.下地幔不同深度处温度的横向频率分布。

  进一步分析发现，下地幔上部的波速异常主要源于热异常，而最底部地幔中的波速异常则主要源于化学成分的变化（图3）。

图3.不同深度处温度和成分的变化对Vs异常的贡献。

  研究结果显示，LLSVPs与周围地幔相比在下地幔底部密度更高，在~2700公里深度以上密度更低。此外，LLSVPs相比周围地幔温度更高，更加富集铁和布里奇曼石（图4），这支持了LLSVPs可能源自地球早期形成的基底岩浆洋的假说。

图4.2680 km与2800 km深度下横向成分与温度分布。

  本研究的发现提供了关于下地幔物质组成和温度结构的关键见解，在很大程度上增进了我们对地球深部结构的认识，并将对有关地球的形成、演化与动力学等方面的研究产生重要影响。

  论文第一作者为博士生邓忻，通讯作者为吴忠庆教授。合作者包括南京大学阮友谊副教授、中科院精密测量科学与技术创新研究院倪四道研究员、中国科学技术大学王文忠研究员、徐印涵博士、博士生赵亚婕以及美国加州大学圣地亚哥分校博士生郝尚钦。中国科大为该论文第一单位。该项研究受到国家重点研发计划、国家自然科学基金项目以及中央高校基本科研业务费专项资金等项目的资助。计算工作主要在中国科大超级计算中心的超级计算机上进行。

论文链接：

  Deng, Xin; Xu, Yinhan; Hao, Shangqin; Ruan, Youyi; Zhao, Yajie; Wang, Wenzhong; Ni, Sidao; Wu, Zhongqing. (2023) "Compositional and thermal state of the lower mantle from joint 3D inversion with seismic tomography and mineral elasticity", Proceedings of the National Academy of Sciences, 120 (26). https://doi.org/10.1073/pnas.2220178120.