《Frontiers in Earth Science》：Seismic anisotropy reveals stress rotation caused by magmatic-tectonic interactions in the Natron rift
编辑推荐：
Oldoinyo Lengai火山位于纳特龙裂谷内，这是东非裂谷系统（East African Rift System, EARS）中一个年轻的岩浆活动段，此处构造伸展与岩浆过程紧密耦合。理解应力、变形和岩浆运移在此类环境中的相互作用仍是一项关键挑战。本研究利
Oldoinyo Lengai火山位于纳特龙裂谷内，这是东非裂谷系统（East African Rift System, EARS）中一个年轻的岩浆活动段，此处构造伸展与岩浆过程紧密耦合。理解应力、变形和岩浆运移在此类环境中的相互作用仍是一项关键挑战。本研究利用局部剪切波分裂测量来刻画地震各向异性，并约束纳特龙裂谷的应力与变形分布。研究人员分析了密集SEISVOL地震台网的约4400个台站-事件对，识别出平均延迟时间约0.12秒、优势快轴方向约21°N的各向异性，与区域伸展应力场总体一致。然而，显著的横向变化揭示了纳特龙盆地、Gelai以东地区以及Oldoinyo Lengai附近的岩浆活动区之间存在不同的各向异性域。在所有区域，各向异性均由约5–6公里的浅部地壳层主导，归因于应力定向排列的裂隙系统，并且各向异性分层和横向非均质性提供了额外复杂性。在岩浆活动区，观测到与区域模式的显著偏离，包括快轴方向的局部旋转和延迟时间增大。随深度分析表明，在中地壳深度（约11–15公里）存在额外的各向异性贡献，与先前推断的岩浆存储区一致。分裂参数的时间变化进一步与火山活动增强期相关，表明碳酸岩与硅质岩浆管道系统强烈耦合，相应的熔体和流体运移动态地改变应力场和裂隙网络。研究人员将这些观测解释为受岩浆影响的转换带的表达，其中应力旋转、剪切调节和集中岩浆上升相互作用。结果表明，地震各向异性既能捕捉区域应力，也能捕捉岩浆过程对其的局部改造，突显了岩浆-构造反馈在年轻大陆裂谷变形与分段控制中的关键作用。
研究背景与意义
大陆裂谷作用的核心问题之一在于区分岩浆与构造的相对贡献，以及这一贡献如何随裂谷演化而变化。纳特龙裂谷位于东非裂谷系统（East African Rift System, EARS）北坦桑尼亚分歧区，是一个年轻的岩浆型裂谷段，构造伸展与岩浆活动紧密耦合。该区域拥有全球唯一现生碳酸岩（carbonatite）火山Oldoinyo Lengai，以及Gelai、Kerimasi等火山。已有地震层析成像、衰减成像和大地测量研究揭示了复杂的多级岩浆管道系统，也观察到2007年岩墙侵入与纳特龙边界断层滑动之间的力学联系。然而，局部应力场如何响应断层分段、岩浆储库和流体运移仍缺乏高分辨率直接约束。为此，研究人员利用SEISVOL台网的局部剪切波分裂（shear-wave splitting）测量，对纳特龙裂谷的应力与变形分布进行成像，以理解岩浆与构造反馈在年轻大陆裂谷中的作用。该研究发表在《Frontiers in Earth Science》。
研究人员使用由Reiss与Rümpker构建的SEISVOL台网（编码9J）记录的2019年3月至2020年5月连续地震数据，共26个宽带台站，数据取自GEOFON存档；地震事件目录和P波、S波到时沿用Reiss等人先前定位结果。主要技术方法包括：应用S变换和协方差矩阵特征分解进行时频域极化分析，以分离S波；利用自动SplitRacer工具包，结合特征值最小化与能量最小化方法反演分裂参数；依据分裂强度、椭圆率、能量降低等参数对测量进行空、好、一般质量分类；对台站开展双层各向异性联合分析，并采用bootstrap统计与聚类分析评估模型稳定性；同时通过网格搜索拟合视各向异性随深度分布模型，约束各向异性层的深度范围。
研究结果
3.1 单次测量结果
研究人员从13528个台站-事件对中筛选出4226个可观测分裂事件，其中1676个被归类为良好质量。区域平均延迟时间为0.12秒，良好测量的平均值为0.10秒；405个全部测量及96个良好测量显示大于0.25秒的大延迟时间。优势快轴方向为北北东向，但表现出系统性空间变化。据此将研究区划分为三个各向异性域：A为纳特龙盆地，B为Gelai以东区域，C为Oldoinyo Lengai与岩浆带。
3.2 纳特龙盆地
纳特龙盆地的多数测量显示近南北向快轴，从西部裂谷肩附近的约5°N向东至Gelai附近旋转到约-10°N，并出现约70°N至130°N的近东西向快轴。研究区地震活动相对较低且震源更深，导致采样不足。近边界断层台站的双层各向异性分析显示，浅层快轴约-20°N，与边界断层走向近平行；下层快轴近东西向，平行于伸展方向，指示地壳内存在各向异性分层。
3.3 Gelai以东区域
Gelai以东区域的优势快轴方向垂直于伸展方向，并与Gelai火山的长轴近平行。在约2.62°S、2.675°S和2.71°S处，出现间隔规则的、与Gelai长轴垂直的窄带快轴方向，研究人员将其解释为局部剪切带，与区域内少量走滑型震源机制一致。双层联合分析显示，上层快轴垂直于Gelai地形轮廓，下层平行于裂谷轴，指示地形载荷对浅层应力场的局部改造。
3.4 Oldoinyo Lengai与岩浆带
在Gelai以南的转换区域，快轴方向发生大尺度旋转，出现垂直于裂谷轴的取向，即平行于伸展方向；该取向在约36°E、-2.75°N附近最为集中，该处曾发生岩浆侵入。Oldoinyo Lengai附近快轴以平行裂谷方向为主，但整体复杂性更高，反映了岩浆管道系统对应力场的干扰。
3.5 深度变化
3.5.1 视各向异性深度估计
通过计算视各向异性（apparent anisotropy）随深度的分布，研究人员发现浅部约5–6公里范围内各向异性最强，估计值介于约6%至26%之间；模型拟合给出浅层各向异性层厚度约5–9公里，峰值各向异性约12%–31%。在岩浆活动区，最大视各向异性估计在11–12公里深度、平均估计在14–15公里深度出现偏离模型衰减趋势的增大，与先前推断的中地壳岩浆存储区位置一致。
3.5.2 双层各向异性
24个台站中，有7个台站采用双层模型比单层模型的横向能量最小化改进超过10%，6个台站达到超过40%的总能量降低，表明部分台站下方存在明显的各向异性分层。Gelai以东区域的双层结果较为稳健，但总体能量降低相对较低，暗示该区域存在更复杂的横向各向异性结构。
讨论与结论
讨论部分指出，纳特龙裂谷以约1.67毫米/年的速率伸展，伸展方向约为101°。在最小主应力σ
3水平、裂隙垂直伸展方向开启的背景下，区域应力场预期产生约11°N的快轴方向，与观测的优势取向一致。浅部裂隙系统是各向异性的主要来源。纳特龙盆地近边界断层的分层各向异性可能与深部断裂网络、地形载荷引起的重力应力有关。Gelai以东的窄带快轴旋转指示分段正断层之间的走滑剪切调节。Oldoinyo Lengai区域的时间变化显示，火山活动增强期与中地壳岩浆房上方的延迟时间增大相对应，说明碳酸岩与硅质岩浆管道系统强烈耦合，熔体和流体运移动态改变裂隙网络和局部应力场。最突出的观测结果是在Oldoinyo Lengai与Gelai以南之间形成局部转换带（transfer zone）：由于边界断层和Gelai下方活动断层沿走向错列，差异运动导致转换带内张应力轴旋转至近平行于裂谷方向，从而产生垂直裂谷的快轴方向，并使岩浆趋于在该带汇聚。地形应力模型虽能解释局部快轴取向，但预测的应力幅值与观测到的强各向异性不匹配，因此需要叠加岩浆与流体过程的贡献。
研究结论翻译如下：本研究为东非裂谷系统纳特龙裂谷段中构造变形、应力传递与岩浆作用之间的相互作用提供了新约束。利用SEISVOL台网的局部剪切波分裂测量，研究人员获得了普遍存在但空间可变的 seismic anisotropy（地震各向异性），它同时反映了区域伸展应力场及其局部重组。各向异性主要由约5–6公里深的浅部地壳层主导，符合区域伸展形成的应力定向裂隙系统；但系统性横向变化表明，断层分段引起的应力传递叠加在区域模式之上。从边界断层主导变形向Gelai下方活动增强的过渡，表现为快轴方向旋转和局部剪切带存在，指示应变在错列断层段之间再分配。在Gelai以东和Oldoinyo Lengai岩浆系统方向，各向异性记录了从伸展主导到剪切影响的变形转换。该调节带具有增强的各向异性变化，反映分布式破裂和应力条件演化。随深度分析进一步显示中地壳深度存在额外各向异性贡献，与岩浆存储区一致；分裂参数的时间变化与火山活动增强期相关，证明熔体和流体运移动态改变裂隙网络和局部应力场。总体而言，研究结果证明地震各向异性能够以高时空分辨率捕捉应力再分配、断层分段和岩浆过程的耦合效应。在纳特龙裂谷中，岩浆-构造相互作用在地壳变形塑造中发挥积极作用，应力转换带成为应变调节和潜在岩浆汇聚的关键位置，在早期大陆裂谷阶段尤为重要。
