过去的几十年里，全球尺度的地震成像研究揭示了下地幔底部的复杂结构。其中，位于非洲和太平洋板块之下核幔边界（core-mantle boundary，简称cmb）处的两个大型低剪切波速省（large low-shear velocity provinces，简称llsvps），尤为吸引广大学者们的目光。llsvps横向跨度可达数千公里（he and wen, 2012; ni and helmberger, 2003），占据了cmb约30%的面积（图1，garnero et al., 2016; torsvik et al., 2010）。纵向形态上，既有研究认为llsvps是一整个大的低速异常体，并由cmb向上延伸几百甚至上千公里（tsekhmistrenko et al., 2021）；也有学者提出llsvps是由多条独立地幔柱聚成的簇状体，且这些地幔柱源于cmb上一个低速、可能是较高密度的薄层（小于200 km厚）（davaille and romanowicz, 2020）。此外，源自深部地幔llsvps等区域的地幔柱可以将下地幔物质输送至地表（hofmann, 1997; morgan, 1971）。
图1 全球大型低剪切波速省 (llsvps)、大火成岩省、金伯利岩及与深源地幔柱相关的活跃热点分布。红色实线勾画了非洲和太平洋llsvps的轮廓（torsvik et al., 2010）
但对地幔柱的形态，以及地幔柱作为热物质通道的动力学过程的认识还存在颇多争议。基于地球动力学模拟实验，griffiths and campbell（1990）提出了经典的蘑菇状柱头-细长柱尾（直径小于200 km）的经典热地幔柱模型。而新的层析成像结果表明，地幔柱在下地幔可能是粗柱状（直径约500-800 km）的结构（french and romanowicz, 2015）。经典模型要求下地幔底部存在热异常，而单纯的热异常却不能形成粗柱状结构，需要热-化学异常的参与。进一步，针对位于非洲板块下的地幔柱，英国牛津大学地球科学学院的tsekhmistrenko 学者等提出了树形结构模型（详细内容可参阅报道《》）。相比下地幔粗柱状的地幔柱模型，该模型更为细致地刻画了1000 km以上的结构，并指出次生地幔柱的结构。这一模型得到了前人研究不同程度的支持（chang and van der lee, 2011; emry et al., 2019; liu and leng, 2020）。
近期，wamba et al.（2023）通过全波形反演（full-waveform inversion, 简称fwi）约束了位于印度洋西部的reunion和comores等热点及其周边区域（图2a）的下地幔速度结构（约800 km深度至cmb），结合该区域已有的地壳及上地幔速度模型semindo（深至约1200 km，wamba et al., 2021），他们提出新的全地幔尺度三维速度结构模型（semindo-wm3）。模型的二维纵向剖面（图2b-图2e）显示，该区域自cmb至约1000 km的深度范围存在至少三条相互独立、近直立的低速柱状体；在约1000至660 km的深度范围，这些柱状体向水平方向弯折并延展，形成水平层状低速体；自660 km以上，可以观测到更为直立的低速通道与地表热点火山相连接。三维模型更好地展示了这些结构的完整形态（图3），可以看到三条独立的低速柱状结构在约1000 km的深度处与水平低速层状延展区（low-seismic velocity ponding zone，简称lsvpz）相接。lsvpz的形成可能与该层较低的粘滞度有关（liu and leng, 2020; rudolph et al., 2015）。
图2 （a）研究区主要热点（黑色实心圆）及地质构造分布。cr: carlsberg ridge（卡尔斯伯格洋中脊）; cir: central indian ridge（中印度洋中脊）; swir: southwest indian ridge（西南印度洋中脊）; seir: southeast indian ridge（东南印度洋中脊）; rtz: rodrigues triple junction（罗德里格斯三联点）; asz: andaman subduction zone（安达曼俯冲带）; ssz: sumatra subduction zone（苏门答腊俯冲带）; jsz: java subduction zone（爪哇俯冲带）。（b）沿reunion和crozet热点连线的地幔剪切波速度扰动剖面图。（c）沿reunion和kerguelen热点连线的剖面图。（d）沿reunion和marion热点连线的剖面图。（e）沿comores和marion热点连线的剖面图（dongmo wamba et al., 2023）
图3 reunion下方cmb至850 km深度范围三维剪切波速度模型。（a）自东向西视角；（b）自南向北视角；（c）自北向南视角。顶部横截面的白色粗线代表地表地质构造的投影，红色粗线代表印度洋中脊的投影，绿色锥形指示了reunion和comores热点的位置（wamba et al., 2023）
大体来说，热管道系统从深部地幔到地表热点分为五层（图4）：（1）自cmb至约1000 km深度的粗柱状通道，将llsvp的热物质向上输送；（2）自约1000 km至660 km深度的水平层状延展区，接收来自下方的热物质，形成供给源；（3）自约660 km至250 km深度的细柱状通道，将热物质继续向上输送；（4）自约250 km至100 km深度的水平层状延展区，接收下方热物质，形成新的供给源；（5）自约100 km深度至地表的细枝状通道，将热物质输送至地表。值得注意的是，位于层（1）、（3）和（5）的近垂直通道并不一定处于同一水平位置，事实上他们可能会在横向上发生较大的偏移。这种垂直柱状通道-水平层状延展区交错的结构（图4），类似于地壳内火山系统的岩墙与岩床交织的结构，只不过尺度更大，贯穿了整个地幔。该结构模型与树形结构模型（tsekhmistrenko et al., 2021）存在一定的相似性：两个模型均指出自cmb上涌的物质在下地幔顶部（约1000 km以上）出现水平方向的延展，并在约660 km附近由于粘滞度的跳变形成直径较小的次生柱状通道。不同的是，树形结构模型认为cmb以上数百至上千公里的深度范围是一个分叉、树根状的结构，而该模型则认为在这个深度范围存在多条独立的柱状通道；同时树形结构模型并没有揭示出层（4）的水平层状延展区和层（5）的直径更小的枝状通道。
图4 地幔热管道系统概念图（wamba et al., 2023）
基于地幔速度结构模型（semindo-wm3）提出的热管道系统模型，有助于解析区域性地震学问题。例如，已有地球化学研究提出，reunion热点火山与中印度洋中脊（central indian ridge，简称cir）可能存在内在关联（furi et al., 2011）。新的热管道系统模型指出，此二者可能在下地幔源于同一个地幔柱，进一步支持了以上观点。二维纵向剖面显示，在reunion下方约660至250 km深度范围（层（3））存在一个直立的柱状低速通道，与上方的lsvpz（层（4））相连（图2b-2d）。在200 km深度的横截面上，可以观测到层（4）的lsvpz向东侧的cir和西侧的comores热点（mascareignes盆地）方向延伸（图5a）。结合地球化学观测，这表明reunion热点和cir在下地幔可能由同一个地幔柱提供热物质，只不过由于lsvpz的水平延展性质，导致热物质在地表不同位置喷出。那么reunion和comores热点是否在下地幔相连、由同一地幔柱供给热物质呢？或者说向西延展的lsvpz是否与地表的comores热点相连，并向其输送深部热物质呢？地球化学研究显示，reunion和comores热点处的玄武岩在地球化学性质上存在一定的差异，反映向西延展的lsvpz可能并没有与comores热点源区连接，而comores热点可能与其西侧的东非裂谷系是同源的（class et al., 2005）。
图5 与图3类似，展示自南向北视角reunion下方中地幔至上地幔三维剪切波速度模型。（a）1700 km至200 km深度范围；（b）1700 km至300 km深度范围；（c）1700 km至400 km深度范围（wamba et al., 2023）
新的热管道系统模型也可为解析全球性的地震学焦点争议提供新的观测和证据。例如，该模型支持llsvps是由多条独立地幔柱组成的，而非一个整体的大型低速异常体；热点火山作用是由cmb处的热物质通过管道系统输送至地表而引起的，而不仅是上地幔对流的地表响应等。同时也应注意，受区域观测数据覆盖限制，该模型的普适性还有待进一步针对不同地区的研究加以验证。
主要参考文献
chang s j, van der lee s. mantle plumes and associated flow beneath arabia and east africa[j]. earth and planetary science letters, 2011, 302(3-4): 448-454.
class c, goldstein s l, stute m, et al. grand comore island: a well-constrained “low 3he/4he” mantle plume[j]. earth and planetary science letters, 2005, 233(3-4): 391-409.
davaille a, romanowicz b. deflating the llsvps: bundles of mantle thermochemical plumes rather than thick stagnant “piles”[j]. tectonics, 2020, 39(10): e2020tc006265.
emry e l, shen y, nyblade a a, et al. upper mantle earth structure in africa from full‐wave ambient noise tomography[j]. geochemistry, geophysics, geosystems, 2019, 20(1): 120-147.
french s w, romanowicz b. broad plumes rooted at the base of the earth's mantle beneath major hotspots[j]. nature, 2015, 525(7567): 95-99.
füri e, hilton d r, murton b j, et al. helium isotope variations between réunion island and the central indian ridge (17°-21°s): new evidence for ridge–hot spot interaction[j]. journal of geophysical research: solid earth, 2011, 116(b2).
garnero e j, mcnamara a k, shim s h. continent-sized anomalous zones with low seismic velocity at the base of earth's mantle[j]. nature geoscience, 2016, 9(7): 481-489.
griffiths r w, campbell i h. stirring and structure in mantle starting plumes[j]. earth and planetary science letters, 1990, 99(1-2): 66-78.
he y, wen l. geographic boundary of the “pacific anomaly” and its geometry and transitional structure in the north[j]. journal of geophysical research: solid earth, 2012, 117(b9).
hofmann a w. mantle geochemistry: the message from oceanic volcanism[j]. nature, 1997, 385(6613): 219-229.
liu h, leng w. plume‐tree structure induced by low‐viscosity layers in the upper mantle[j]. geophysical research letters, 2020, 47(1): e2019gl086508.
morgan w j. convection plumes in the lower mantle[j]. nature, 1971, 230(5288): 42-43.
ni s, helmberger d v. seismological constraints on the south african superplume; could be the oldest distinct structure on earth[j]. earth and planetary science letters, 2003, 206(1-2): 119-131.
rudolph m l, leki？ v, lithgow-bertelloni c. viscosity jump in earth’s mid-mantle[j]. science, 2015, 350(6266): 1349-1352.
torsvik t h, burke k, steinberger b, et al. diamonds sampled by plumes from the core–mantle boundary[j]. nature, 2010, 466(7304): 352-355.
tsekhmistrenko m, sigloch k, hosseini k, et al. a tree of indo-african mantle plumes imaged by seismic tomography[j]. nature geoscience, 2021, 14(8): 612-619.
wamba m d, montagner j p, romanowicz b. imaging deep-mantle plumbing beneath la réunion and comores hot spots: vertical plume conduits and horizontal ponding zones[j]. science advances, 2023, 9(4): eade3723.（）
wamba m d, montagner j p, romanowicz b, et al. multi‐mode waveform tomography of the indian ocean upper and mid‐mantle around the réunion hotspot[j]. journal of geophysical research: solid earth, 2021, 126(8): e2020jb021490.
（撰稿：王拓，，/岩石圈室）
