地幔柱动力学模拟.docxVIP

PAGE1/NUMPAGES1

地幔柱动力学模拟

TOC\o1-3\h\z\u

第一部分地幔柱形成机制 2

第二部分热-化学耦合模型 4

第三部分数值模拟方法综述 8

第四部分地幔对流边界条件 13

第五部分柱头羽状结构演化 17

第六部分板块运动相互作用 21

第七部分地表火山响应特征 25

第八部分模拟结果验证策略 29

第一部分地幔柱形成机制

地幔柱形成机制是地球深部动力学研究中的核心科学问题之一，其本质涉及地球内部热-化学对流过程、相变效应、地核-地幔边界（Core-MantleBoundary,CMB）热结构演化以及岩石圈-软流圈相互作用等多个物理与地球化学过程的耦合。地幔柱通常被定义为起源于地核-地幔边界或下地幔底部、以热异常和/或化学异常为特征、向上穿透整个地幔并最终影响地表地质活动的柱状上涌流体结构。其形成机制主要可归纳为热驱动型、化学浮力驱动型以及二者耦合驱动型三种模式。

首先，热驱动型地幔柱模型认为，地幔柱源于CMB附近局部区域的显著热积累。由于地核向地幔持续释放热量，若局部区域热传导效率较低或存在热边界层增厚现象，则会在CMB上方形成热边界层（ThermalBoundaryLayer,TBL）。当该边界层内温度梯度超过临界瑞利数（Rayleighnumber）所对应的对流阈值时，热物质将因浮力失稳而发生上涌，形成热地幔柱。数值模拟研究表明，在典型地幔黏度（1021–1023Pa·s）和热扩散系数（~10??m2/s）条件下，CMB处热异常幅度需达到200–300K以上，才足以触发稳定地幔柱的生成。此外，地幔相变界面（如660km不连续面）对地幔柱上升具有调制作用：部分热柱可能在此界面受阻并横向扩展，形成“蘑菇头”结构；若热浮力足够强，则可穿透相变带继续上涌至浅部地幔。

其次，化学浮力驱动机制强调地幔中存在密度低于周围地幔的化学异常体（如再循环洋壳、亏损地幔残留物或富集组分），其低密度特性可提供额外浮力，促进上涌。地球化学证据显示，热点火山岩（如夏威夷、冰岛）常具有不同于正常洋中脊玄武岩（MORB）的同位素特征（如高3He/?He比值、特定的Sr-Nd-Pb同位素组成），暗示其源区包含古老、未充分混合的化学储库。此类化学异常体可能源于早期地幔分异残留、俯冲板片再循环产物或地核-地幔相互作用形成的超低速带（Ultra-LowVelocityZones,ULVZs）物质。在数值模拟中，引入化学组分差异后，即使温差较小（100K），只要密度差达到0.5%–1%，亦可驱动柱状上涌。值得注意的是，纯化学浮力驱动的地幔柱上升速度通常较慢，且形态更弥散，难以解释快速、集中式热点活动。

第三，热-化学耦合机制被认为是更符合实际地球条件的综合模型。在此框架下，CMB附近的热异常与化学异常共存并相互增强：热异常降低局部黏度，促进化学异常体的聚集与上涌；而化学异常体提供的额外浮力又可降低对热异常幅度的要求，使地幔柱在更广泛的参数空间内得以形成。高分辨率三维球坐标对流模拟表明，当同时考虑热与化学浮力时，地幔柱可呈现周期性脉冲式喷发特征，其头部直径可达数百公里，尾部直径约50–100km，上升速率约为5–15cm/yr，与地震层析成像观测到的深部低速异常体尺度及位置具有较好一致性。例如，太平洋下方LLSVPs（LargeLowShearVelocityProvinces）边缘区域常被视为地幔柱的主要发源地，其既具有显著的热异常（ΔT≈200–400K），又富含铁、钛等重元素导致的化学密度异常。

此外，地幔柱的形成还受到全球板块构造格局的调控。俯冲带通过冷却上覆地幔、改变应力场分布及诱发地幔回流，可间接影响地幔柱的定位与稳定性。部分研究指出，超级地幔柱（如白垩纪非洲下方）的爆发可能与同期全球板块重组事件相关。同时，地幔黏度的径向与横向非均匀性（如660km相变导致的黏度跃升、岩石圈刚性盖层的抑制效应）亦显著调控地幔柱的上升路径与地表表现形式。

综上所述，地幔柱的形成机制是一个多因素

第二部分热-化学耦合模型

关键词

关键要点

热-化学耦合模型的基本理论框架

1.热-化学耦合模型整合了地幔中热对流与化学组分演化两个物理过程，通过Navier-Stokes方程、能量守恒方程及物质输运方程构建多场耦合控制体系。该模型强调温度场与成分场之间的非线性反馈机制，例如放射性生热率依赖于特定元素（如U、Th、K）的浓度分布，而黏度则同时受温度和化学成分（如铁镁硅酸盐比例）调控。

2.在数值实现上，通常采用有限元或有限体积方法离散偏微分方程组，并引入追踪粒