堪察加半岛下方地幔过渡带（MTZ）约束板片穿透与板片边缘地幔流：来自P-to-S接收函数成像的证据

《Geochemistry, Geophysics, Geosystems》：Mantle Transition Zone Constraints on Slab Penetration and Slab-Edge Mantle Flow Beneath the Kamchatka Peninsula

　　

大量地震学研究已报道在堪察加北部部分地区之下缺乏可分辨的俯冲板片，尽管该异常的机制仍存在争议。研究人员利用1,696条P-to-S接收函数（receiver functions, RFs）对堪察加半岛及周边地区的地幔过渡带（mantle transition zone, MTZ）进行成像。在堪察加西南部之下，增厚的MTZ和显著的520-km不连续面支持冷的白垩纪太平洋板片相关物质，并可能在MTZ内发生水化。相反，在堪察加北部之下观测到显著的MTZ减薄，而先前层析成像（tomography）研究也显示浅部上地幔中板片缺失。这反对在采样的MTZ中存在大型冷拆离板片残余。结合先前的（seismic anisotropy）、层析成像和地球动力学（geodynamic）约束，研究人员的成果支持一种模型：热浮力地幔物质绕板片边缘流动并上升进入上地幔，形成半环状地幔流（semi-toroidal mantle flow）。推断的板片边缘地幔流（slab-edge mantle flow）可能至少延伸至MTZ深度，并为堪察加北部之下板片边缘的渐进改造和弱化提供了合理背景，同时允许板片在更西南方向穿透。这些观测突出强调了板片–地幔相互作用（slab–mantle interaction）在控制汇聚边缘板片命运中的一级作用。

研究背景、问题与意义：堪察加俯冲带位于西北太平洋，是白垩纪太平洋板片以约8 cm/year平均速率俯冲到鄂霍次克板片之下的活跃汇聚边缘。长期汇聚形成堪察加半岛及三条平行海沟的火山带：Eastern Volcanic Front、Central Kamchatka Depression和Sredinny Range。地震成像显示太平洋板片俯冲似乎终止于阿留申–堪察加交界附近，在北堪察加不活动火山弧之下未检测到清晰板片残余。该板片缺失机制存在争议：Levin et al. 基于约10 Myr构造与火山演化归因于海洋岩石圈拆离（detachment）；其他研究未检测到拆离板片碎片预期的高波速异常，而归因于板片边缘随地幔环流的渐进改造。上地幔层析成像（tomography）与地震各向异性（seismic anisotropy）提供板片边缘地幔流独立证据，但深度范围受层析分辨率限制。地幔过渡带（mantle transition zone, MTZ）由d₄₁₀与d₆₆₀限定，d₄₁₀对应橄榄石到瓦兹利石相变、正Clapeyron斜率；d₆₆₀对应林伍德石到布里奇曼石加方镁石，负Clapeyron斜率。冷异常通常增厚MTZ，热异常减薄；但d₆₆₀在高温下受majorite到perovskite相变控制，正Clapeyron斜率，且水含量也可增厚MTZ。因此需考虑温度、相变体制与挥发分。近期SdP研究显示堪察加和北千岛之下d₆₆₀明显凹陷，支持板片进入下地幔，中央堪察加d₄₁₀抬升和d₆₆₀凹陷导致MTZ增厚，可能由冷太平洋板片及撕裂板片窗附近地幔流热改造。但SdP结果主要沿中深源地震照亮的弧段，北堪察加MTZ约束不足。因此需要更广空间覆盖的接收函数（receiver functions, RFs）研究。

研究方法概要：数据来自EarthScope Consortium Data Services和GFZ Helmholtz Centre for Geosciences的GEOFON data center，包括84个台站，1990—2023年，150°E—168°E、50°N—60°N。采用可变截止震级M_c；经去卷积垂直分量计算径向P-to-S接收函数，质量筛选后获得来自466个地震的1,696条高质量P波RF。使用IASP91模型进行走时校正，在535 km参考深度计算射线piercing points，以1°半径圆Bin内至少10条RF叠加，生成300—800 km深度序列；采用50次bootstrap重采样估计均值与标准差；用TX2019slab和Jiang2009Vp模型进行速度校正。

研究结果：可靠测量在55个Bin获得d₄₁₀、54个d₆₆₀、54个同时包含二者。表观MTZ厚度197—293 km，平均248±18 km；表观d₄₁₀ 416±13 km，d₆₆₀ 664±11 km。TX2019slab校正后均值d₄₁₀ 415±13、d₆₆₀ 666±12、MTZ 251±18 km；Jiang2009Vp校正后d₄₁₀ 417±13、d₆₆₀ 666±17、MTZ 249±21 km。划分为A、B、C三个子区。

子区A：西东向，无俯冲相关中深源地震，穿过Sredinny Range、Central Kamchatka Depression至堪察加–阿留申交界。d₄₁₀平均凹陷约17 km；表观d₆₆₀ 656±11 km，比参考浅约4 km；MTZ比250 km薄约21 km。TX2019slab校正后d₄₁₀ 423±8、d₆₆₀ 655±11、MTZ 232±12 km；Jiang2009Vp校正后422±8、651±13、229±12 km。两种校正一致，确认A区MTZ减薄。

子区B：东北–西南向，平行俯冲带，包括堪察加半岛最南部。MTZ 235±5 km，d₄₁₀平均凹陷10 km，d₆₆₀平均抬升5 km。TX2019slab校正：d₄₁₀ 423±4、d₆₆₀ 661±7、MTZ 238±5 km；Jiang2009Vp：425±5、655±4、230±5 km。

子区C：MTZ平均278±12 km，d₄₁₀抬升约22 km，d₆₆₀凹陷约6 km。TX2019slab校正：d₄₁₀ 388±6、d₆₆₀ 671±11、MTZ 283±13 km；Jiang2009Vp：397±6、684±19、287±17 km。C区与A、B形成鲜明对比。

讨论总结：
4.1 半环状板片边缘地幔流与MTZ减薄：A区MTZ减薄18–21 km，速度校正后仍存在，最可能反映MTZ高温而非速度结构假象。按dV_p/dT=?4.8×10^?4 km s^?1K^?1，d₄₁₀ Clapeyron +2.9 MPa K^?1，d₆₆₀ ?1.3 MPa K^?1，减薄对应高温异常约150–180 K，不足以使660 km主相变转为majorite控制，对应MTZ内平均负P波速度异常约0.7–0.9%，与区域远震层析低速异常大体一致。深地幔柱起源不太可能，因北堪察加无年龄递进热点轨迹、下地幔热通道和plume-type地球化学特征。该温度异常与全球俯冲相关板片边缘上涌环境范围一致。东部MTZ减薄由d₄₁₀凹陷和d₆₆₀抬升，热效应延伸至MTZ大部分；西部由d₄₁₀凹陷主导，d₆₆₀近参考，热异常主要限于上MTZ。地球动力学模型表明海沟后撤、板片撕裂或板片丢失可诱导板片边缘局部地幔上涌。独立地震各向异性研究记录北堪察加推断板片边缘周围半环状快方向，通常解释为浅上地幔板片边缘驱动地幔流。其与MTZ减薄空间对应，表明热异常可能至少延伸至MTZ。热浮力地幔物质绕板片边缘循环并上升进入上地幔，与各向异性及板片边缘地幔环流模型一致。Cui et al. 推断热地幔流绕撕裂板片窗可能热改造MTZ深度太平洋板片；但其研究记录冷板片主导的增厚MTZ，而本研究北板片边缘显示减薄，指示相对较暖热状态。

4.2 北堪察加板片缺失的可能动力学机制：地质和火山记录表明太平洋板片曾存在于北堪察加之下，Sredinny Range上新世火山作用为证。后续火山重组归因于板片拆离或板片边缘被地幔环流渐进改造。若拆离，拆离碎片应因浮力屏障和d₆₆₀附近黏度增加停滞在MTZ，产生MTZ增厚。全球综合显示俯冲板片可在MTZ停滞较长时间后进入下地幔。因此大型冷拆离板片残余在采样MTZ应导致增厚。观测却为减薄，不支持拆离主导。小或横向位移碎片不能完全排除，但本研究MTZ不支持。MTZ减薄、板片边缘各向异性和adakitic火山作用综合更符合周围地幔流渐进改造板片边缘。物理机制是板片边缘与周围地幔流持续热机械相互作用。Kamchatka corner热模拟表明先前减薄和撕裂板片边缘易再加热。三维地球动力学模型显示板片边缘可发展快速地幔环流。Calabrian slab也提出类似热机械破坏。RF不直接约束该过程速率或历史细节。

4.3 西南堪察加板片穿透d₆₆₀的证据：C区MTZ显著厚于250 km参考，d₄₁₀抬升、d₆₆₀凹陷，与Cui et al. 观测的大d₆₆₀凹陷一致。层析也显示高波速异常延伸进MTZ，表明白垩纪太平洋板片保留足够密度穿透MTZ并在西南堪察加到达其底部，导致局部冷地幔。但RF主要约束MTZ和d₆₆₀附近结构，不直接成像更深下地幔板片几何；层析提供西北太平洋d₆₆₀以下板片相关结构独立证据，SS precursor stacks成像鄂霍次克海以西最上部下地幔中巨型加厚俯冲洋壳。C区观察到明确P₅₂₀s与d₅₂₀相关，区别于A、B。d₅₂₀对应瓦兹利石到林伍德石相变，正Clapeyron斜率大于d₄₁₀，对中MTZ热扰动敏感。无水条件下该相变阻抗差小、两相共存区间宽，难以检测；C区P₅₂₀s振幅大，局部接近P₄₁₀s和P₆₆₀s，难以仅用相变解释。叠加深度序列显示P₅₂₀s清晰一致，非叠加假象。两个互补机制：太平洋板片通过俯冲前洋岩石圈水化和下降过程脱水将水带入MTZ；实验和模拟表明板片来源水缩窄瓦兹利石–林伍德石两相共存区间，使速度梯度增至三倍，增强d₅₂₀地震表现；MTZ储水也可降低d₄₁₀阻抗差并促进MTZ增厚。整体MTZ组成变化也可改变瓦兹利石–林伍德石地震对比，因为520-km不连续面强度对橄榄石多型丰度和组成敏感。总之C区增厚MTZ、抬升d₄₁₀、凹陷d₆₆₀、异常强d₅₂₀与冷且水化的太平洋板片延伸穿过大部分MTZ一致。相反A、B区d₅₂₀弱或缺失；实验表明高温使瓦兹利石–林伍德石相变移向更深并降低地震可检测性，支持A、B区MTZ高温异常，并与不存在类似C区的冷、水化太平洋板片一致。

结论翻译：基于接收函数的地幔过渡带（MTZ）成像揭示堪察加半岛之下板片行为沿走向的显著差异，这些差异与板片–地幔相互作用紧密相关。北堪察加之下一致的18–21 km MTZ减薄表明温度升高与板片边缘地幔流有关，该地幔流可能向下至少延伸至MTZ深度。这种绕板片边缘的热浮力地幔环流与半环状地震各向异性型式一致，并为板片边缘渐进改造和弱化提供合理背景，从而可在不要求采样MTZ中存在大型冷板片残余的情况下，解释不活动北部火山弧之下缺乏可检测板片残余。形成鲜明对比的是，西南堪察加以显著增厚的MTZ（278±12 km）和显著的520-km地震不连续面为特征，二者共同支持冷太平洋板片相关物质延伸穿过大部分MTZ。异常强的d₅₂₀还可能反映与俯冲板片相关的水化或组成非均质性。这些观测强调板片–地幔相互作用在控制汇聚边缘板片命运中的关键作用，为改进板片边缘俯冲和深部地幔过程的地球动力学模型提供关键约束。