地震b值(b-values)在台湾板块边界上的空间、时间和深度依赖性变化

《Tectonophysics》：Spatial, temporal, and depth-dependent variability of earthquake b-values across the Taiwan plate boundary

　　

• 在三十年间，b值(b-values)没有明显的长期时间趋势。
• 台湾地区b值(b-values)存在强烈的空间和深度依赖性异质性。
• b值(b-values)与地震聚类或生产力之间的耦合较弱。
• 通过深度分层b值映射(b-value mapping)解析了不同的地震成因机制。
• 深度依赖性和结构变化与台湾板块边界沿线的b值(b-value)变异性有关。

**论文解读：台湾板块边界地震b值(b-value)的时空与深度依赖性变异性**

**研究背景与问题**
古登堡-里克特（Gutenberg-Richter）频率-震级关系是观测地震学的基础，通过b值参数描述地震发生随震级增大呈指数衰减的规律（Gutenberg和Richter，1944）。b值长期以来被解释为地壳和岩石圈物理状态的代理指标，反映差异应力、断层异质性、材料性质和破裂过程的变化（Aki，1965；Scholz，1968，2015，2019）。实验室实验、数值模拟和野外观测一致表明，低b值与高应力水平和均匀断层系统相关，而高b值则倾向于表征异质性、热弱化或流体影响的环境。近期区域尺度研究进一步证明了空间分辨b值分析在调查地震构造变异性、应力异质性和地震潜力方面的有效性（Ahadov和?ztürk，2022；?ztürk和Alkan，2024；Alkan等，2025）。

尽管经过数十年的研究，b值变异性的物理解释仍存在争议。许多研究记录了俯冲带、大陆碰撞带和火山区域的显著空间和时间b值变化（Wiemer和Wyss，2000；Wyss等，1997；Scholz，2015）。台湾的区域研究同样报告了显著的空间、时间和深度依赖性b值变异性，突出了当地地震构造背景的复杂性（Wang等，2015）。Hui等（2020）利用1995–2019年的地震活动性，进一步记录了b值与多重分形谱宽度之间的强正相关，表明地震大小分布与台湾地震活动的时空复杂性共变。多项研究还探讨了台湾b值在大型地震前后的时空变化及其潜在前兆行为（Chan等，2012）。然而，三个主要问题仍限制着可靠解释：第一，b值估计强烈依赖于目录完整性和空间分辨率，难以区分真实物理信号与采样伪影（Woessner和Wiemer，2005）；第二，b值变化常被隐含地等同于地震聚类或余震产率，尽管这些量描述的是地震发生的不同统计性质，可能响应不同的物理控制（Helmstetter等，2007；Zhuang等，2002）；第三，b值的深度依赖性——特别是在脆-韧转换带和俯冲板块内部——很少在系统且内部一致的框架内利用均匀区域目录进行探索（?ztürk和Alkan，2024）。

台湾为解决这些挑战提供了独特的天然实验室。它位于欧亚板块与菲律宾海板块的汇聚边界，在短空间尺度上拥有广泛的构造体制，包括活跃的弧-陆碰撞、逆冲断层、走滑变形、伸展域和深部板块内地震活动性。更重要的是，该区域由中央气象署（Central Weather Administration，CWA）运营的全球最密集地震台网之一监测，自1990年代初以来提供了长期、内部一致的本地震级（ML）目录。这一独特组合使得深度分辨、空间自适应分析成为可能，这在其他板块边界环境中很少可行。

**研究内容与结论**
研究人员利用覆盖1991–2025年的CWA ML地震目录，在统一分析框架内调查了台湾b值的空间、时间和深度依赖性变异性。主要结论包括：b值在参考震级阈值下无显著长期时间漂移；空间上呈现显著侧向异质性，典型值范围约0.6至1.2；b值表现出显著的深度依赖性变化，浅层、中层和深层地震活动具有不同的b值分布和空间模式；b值与地震聚类或产率之间的关联较弱；基于剖面的分析揭示了b值与深度结构的协同变化，反映了台湾的三维构造结构。该研究为理解复杂板块边界的地震行为提供了基准，并强调了在将b值用作地震成因状态指标时需联合考虑深度、空间分辨率和触发指标的重要性。论文发表在《Tectonophysics》（Tectonophysics）。

**主要关键技术方法**
研究人员采用以下关键技术方法：（1）最大似然法（Aki，1965）估计古登堡-里克特b值，使用0.1的震级宽度和500次自助重采样（bootstrap）获取95%置信区间；（2）自适应常数邻域法（Wiemer和Wyss，2000）进行空间b值映射，设定每个网格节点最近N=500个事件，并计算有效平滑半径R_n；（3）深度分层分析：将目录分为浅层（0–30 km）、中层（30–80 km）和深层（80–300 km），分别进行自适应映射（N=300）；（4）最近邻时空-震级聚类框架（Baiesi和Paczuski，2004；Zaliapin等，2008）量化聚类分数和产率；（5）沿剖面分析：定义东西向和南北向两条剖面，宽度±30 km，在25 km距离箱内计算b值。数据来源为台湾中央气象署（CWA）1991–2025年本地震级（ML）目录，仅使用ML≥3.0的事件。

**研究结果**
**4.1 时间变化**
通过重叠365天窗口（步长90天）估计b值时间序列，结果显示1991–2025年期间无明确长期时间趋势，短期波动通常处于95%自助置信区间内。b值普遍低于1.0，符合台湾活跃板块汇聚、高差异应力的构造背景。尽管存在局部下降（如目录末期），但未形成统计显著的长期趋势。

**4.2 空间变化**
自适应常数-N方法（N=500）得到的空间b值图（图3a）显示显著侧向异质性，典型值范围约0.6至1.2。有效平滑半径R_n（图3b）空间变化反映了地震密度，异常b值区域通常对应中等至较小的R_n，表明观测模式并非由数据密度主导。对ML≥3.5的稳健性测试保留了主要空间特征（补充图S2）。

**4.3 深度依赖性**
深度依赖性分析（图4）显示，浅层地壳b值约0.75–0.80，中间深度局部较高，更深层b值变异性增大且置信区间扩大。使用ML≥3.5的重复分析（补充图S1）保持了整体b(z)结构，表明深度趋势并非由低阈值事件控制。深度分层自适应图（图5）显示三个深度区间具有不同的空间模式，无法简单缩放全深度分布。

**4.4 b值与地震聚类/产率的关系**
通过最近邻聚类框架导出的聚类分数和产率代理（图6），与b值的散点图显示弱相关，斯皮尔曼秩相关系数仅约0.05–0.06，表明b值变化与短期聚类或触发效率在空间上解耦。

**4.5 沿剖面分析**
东西向和南北向剖面（图7）显示b值在数十至数百公里尺度上空间连贯，并伴有中位深度和深度离散度的系统变化。东西向剖面捕捉了从西部碰撞域到东部俯冲系统的转变，南北向剖面反映了沿造山带走向的板块几何、地壳厚度和变形风格变化。

**讨论与结论**
讨论部分总结：b值在台湾板块边界的变异性受深度、构造结构和地震体制的共同影响。时间分析无长期漂移，空间分布呈现显著侧向异质性，深度依赖性表明不同深度存在不同的地震成因机制。b值与聚类/产率之间的弱关联支持了b值与触发效率受不同物理控制的观点。剖面分析强调了构造结构对b值变异性的控制。研究局限性包括仅使用本地震级（ML）和目录截断在ML=3.0，但深度依赖性稳健性测试支持主要结果。最终结论翻译如下：

**结论**
利用1991–2025年台湾均匀本地震级地震目录，本研究对复杂板块边界上古登堡-里克特b值的空间、时间和深度依赖性变异性进行了全面评估。主要结论可概括如下：
第一，b值在参考震级阈值下无明确长期时间漂移，表明目录完整性或台网性能的长期变化不主导观测变异性，支持使用聚合数据集进行空间和深度分析。
第二，b值在台湾区域呈现显著侧向异质性。自适应空间映射表明这些模式对空间分辨率稳健，并非数据密度伪影。不同低b值和高b值域（典型值约0.6至1.2）在可比空间尺度上被分辨，突出了地震行为的真实区域差异。
第三，b值表现出显著的深度依赖性变化。浅层、中层和深层地震活动具有不同的b值分布和空间模式，表明地震成因体积物理状态的深度依赖性变化起主要作用。这些深度效应无法通过简单缩放单一空间模式解释，而指向不同深度运行的根本不同的地震成因机制。阈值敏感性测试（ML≥3.5）支持这一结论。
第四，b值与地震聚类或产率之间的关系较弱。低b值区域（通常~0.6–0.8）并不系统对应增强聚类或更高触发效率，表明b值变异性与地震触发在本文考虑的尺度上基本解耦。
第五，基于剖面的分析揭示了b值与深度结构的协同剖面变化，反映台湾的三维构造结构。东西向和南北向剖面的差异强调了板块边界几何、深度结构和变形风格对b值变异性的联合影响。
总体而言，这些结果表明复杂板块边界的b值变异性与深度依赖性和结构变化密切相关，而具体物理机制有待通过独立地球物理和构造观测进一步约束。本文采用的综合框架为解释其他构造活跃区域的b值变异性提供了基准，并强调了在将b值用作地震成因状态指标时需联合考虑深度、空间分辨率和触发指标。