📄 地球科学

The Magnitude–Frequency Kink in the 2020–2024 Noto Peninsula Earthquake Swarm

本研究分析了 2020–2024 年能登半岛地震群，揭示了深部地震活动中存在具有统计学意义的震级–频率拐点，该现象很可能是震级偏差造成的伪影而非物理现象，从而警示在不考虑深度和震级范围的情况下，不应将该目录中的单一 b 值进行解释。

原作者： RamaKrishna Pasupuleti

发布于 2026-10-13

📖 1 分钟阅读☕ 轻松阅读

原作者： RamaKrishna Pasupuleti

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

地震很少孤立发生；它们成群出现，其规模从仅能被灵敏仪器察觉的微小震颤，到足以重塑地貌的巨大破裂不等。几十年来，地震学家一直依赖一条简单的规则来描述这些事件的分布：对于每一次规模的提升，地震的数量会按可预测的幅度下降。这种被称为古腾堡-里希特定律（Gutenberg-Richter law）的模式表明，如果你统计了小规模地震的数量，就可以可靠地估算出可能发生的几次大规模地震。这种下降的陡峭程度由一个被称为 b 值（b-value）的数值来衡量，科学家利用该值来评估地壳内部的应力水平。当这条规则成立时，它为理解地震风险提供了稳定的基础。然而，自然界往往比直线更为复杂，当这种模式发生断裂时，可能预示着地下正在发生不寻常的变化，例如流体的移动或重大事件发生前的应力积累。

在2024年1月1日袭击日本能登半岛的里氏7.5级毁灭性地震发生前的几年里，该地区经历了一场由数千次较小震颤组成的巨大地震群。这些事件发生在2020年至2024年之间，与从地球地幔上升并导致地壳滑动（但未产生大型冲击）的深层流体有关。一位独立科学家 RamaKrishna Pasupuleti 决定检查这些数千次震颤的目录，以观察它们是否遵循标准的地震统计规则。利用一份经过高度精确重新映射、包含超过31,000个事件的列表，研究人员专门针对发生在2021年至2023年间、即大灾难发生前的地震群活动进行了研究。他们根据深度对地震进行了分类，将靠近地表的浅层地震与埋藏在地下13公里以下的地震分离开来。

研究人员发现，浅层地震的表现完全符合预期，遵循着小规模事件频繁而大规模事件罕见的标准直线模式。然而，深层地震却讲述了一个不同的故事。在13公里深度以下，模式发生了剧烈弯曲。根据微小震颤的数量，深层中等规模地震的数量远高于标准规则所预测的数量。具体而言，数据显示，对于深度超过13公里的事件，震级在3级或以上的地震数量大约是通常统计模型计算值的五倍。这并非微小的波动；统计证据如此有力，以至于研究人员可以自信地拒绝“这些深层事件遵循单一简单规则”的假设。相反，深层地震似乎属于两个不同的群体：一个是由极小且频繁的震颤组成的大群体，以及一个不符合常规模式的规模稍大且更频繁的事件群体。

为了确保这一奇特的发现并非源于数据错误或计数方式的偏差，科学家对其发现进行了严格测试。他们检查了在剔除相关的余震簇后模式是否消失，结果并未消失。他们验证了即使仅观察特定类型仪器读数测量的地震，结果依然成立，从而排除了因更换测量工具而导致结果偏差的可能性。他们还确认，无论是在浅层还是深层，模式的断裂点都发生在大约震级1.9左右，尽管深层的异常现象更为显著。这些检查证实，这种异常模式是地震群的一个真实特征，而非目录中的错误。

尽管证明了该模式的存在，研究人员并未止步于宣布其物理意义究竟为何。他们指出，这种弯曲在所有深度都发生在非常特定的震级处，这是当我们测量微小地震大小的方式发生变化时经常出现的特征。有可能用于计算这些微小深层震级大小的仪器在与它们实际释放的能量进行缩放时并不完美，从而在统计中制造了一个人为的弯曲。作者指出，如果这是一个物理现象，例如深层地下断层滑动类型的变化，则需要通过将这些测量值与另一种不易受微小事件误差影响的地震大小测量方法进行对比来予以确认。在进行此类对比之前，这种弯曲虽然是数据中一个稳健的事实，但其真实原因仍是一个悬而未决的问题。

这些发现为科学家未来如何解读地震统计数据提供了至关重要的警示。这项研究表明，如果数据包含混合的不同行为，或者模式在特定规模处发生变化，那么仅仅计算一个数字来描述某一区域的地震频率可能会产生误导。对于能登半岛地震群而言，深层地震并未遵循标准规则，忽视这种断裂可能会导致对可能发生的较大规模事件的错误估计。研究人员总结道，虽然深层地震群的行为显然与浅层不同，但我们必须谨慎，在能够排除“测量工具对最小事件表现异常”的可能性之前，不要假设我们已经掌握了其物理原因。能登地震群的故事是一个清晰、可测量的异常现象，它挑战了我们的标准模型，提醒我们地球深层的过程所产生的模式远比一条简单的直线要复杂得多。

技术摘要：2020–2024年能登半岛地震群中的震级–频率转折点（Kink）

问题陈述
古登堡–里希特（Gutenberg–Richter, GR）定律假设地震震级与频率之间存在线性关系，这是地震学的基础假设，常用于推断微分应力（通过 $b$ 值）并区分地震阶段。然而，观测中经常出现偏离单一线性 GR 定律的情况，这可能源于目录不完整、混合种群或小震震级标度发生变化。2020–2024年发生在能登半岛东北部下方的地震群（在 2024 年 1 月 $M_w$ 7.5 主震之前发生）提供了一个高精度数据集，用于研究这些偏差。虽然流体驱动的地震群以复杂的震级统计特征著称，但本研究旨在探讨能登地震群的震级–频率分布（FMD）是遵循单一 GR 定律，还是表现出随深度变化的统计显著“转折点”（kink），以及这一特征是物理现象还是震级计算产生的伪影。

研究方法
本研究利用了 Yoshida 等人 (2024) 的双差定位目录，该目录包含 31,158 个事件（2003–2024年）。分析重点在于主要地震群阶段（2021–2023年），包含 14,810 个事件， $M_{JMA} \ge 1.2$ ，并将其分为三个深度带：浅层（< 10 km）、中层（10–13 km）和深层（≥ 13 km）。

关键分析协议包括：

完备性估计： 通过最大曲率和拟合优度法确定，在所有深度带中建立统一的完备震级（ $M_c$ ）为 1.2。
模型比较： 测试了三种 FMD 模型：单一 GR 定律、破碎 GR 定理（具有两个斜率和一个转折震级 $M_k$ ）以及两类种群混合模型。使用贝叶斯信息准则（BIC）对模型进行比较。
统计显著性： 通过参数化自助法（200 次模拟）生成似然比统计量的零分布，以评估破碎模型的显著性。
稳健性检查：
- 去簇化（Declustering）： 应用最近邻去簇法（Zaliapin & Ben-Zion, 2013）以移除聚集事件，并测试该转折是否由余震序列驱动。
- 震级类型审计： 审计了 JMA 震级类型（位移型 vs. 速度型），以排除震级计算方法的切换是否为导致转折的原因。
- 离散度分析： 计算了 GR 条件下的震级离散度残差（ $Z_C$ ），以检测超出 $b$ 值所暗示的震级分组之外的偏差。

主要结果

深度依赖的转折： 浅层事件（< 10 km）符合单一 GR 定律（ $b \approx 0.92$ ）。相比之下，深层事件（≥ 13 km）表现出对单一 GR 定律的统计显著偏离。其分布在小震级时较陡（ $M_k$ 以下 $b_1 \approx 1.57$ ），而在大震级时较平缓（ $M_k$ 以上 $b_2 \approx 0.75$ ）。
转折震级与强度： 转折震级（ $M_k$ ）在不同深度保持一致，出现在 $M_{JMA} \approx 1.9$ 附近（95% 区间为 1.8–2.3）。然而，偏离的强度（以 $\Delta$ BIC 衡量）随深度增加而显著增强。对于深层带，破碎模型优于单一 GR 定律， $\Delta$ BIC $\approx -107$ ，且参数化自助法 $p < 0.005$ 。
大震级事件过剩： 深层带包含 37 个 $M \ge 3$ 的事件，而单一 GR 拟合仅预测 7 个事件（95% 预测范围为 2–14）。同样，观测到 12 个 $M \ge 4$ 的事件，而预期值为 0。
稳健性： 在进行最近邻去簇化处理后，以及在仅限于速度型震级（Type V，占相关震级范围内事件的约 98%）进行分析时，该转折依然存在。审计确认，在转折震级附近不存在会导致解释转折的显著震级类型组成变化。
离散度： 深层事件的 $Z_C$ 指标为负值，表明其离散度高于单一 GR 分布所预期的水平，这与存在混合种群或破碎分布的情况一致。

意义与主张
本文确立了深层能登地震群的震级–频率分布在统计上与单一古登堡–里希特定律不一致。这种偏差是目录的一个稳健特征，通过了完备性、去簇化和震级类型切换的测试。

作者明确指出，该转折的物理起源仍是一个开放性问题。作者指出，在不同深度出现共同转折震级是已知现象，即小震级事件中基于振幅的震级与地震矩之间的非线性标度关系。然而，作者认为，已知的 JMA 小事件偏差（低估微震）在理论上会使小震级斜率变平，而观测到的能登深层 FMD 在小震级处反而更陡。

因此，本研究得出结论：虽然该转折是一个记录在案的统计事实，但其解释究竟是物理过程（例如，两种不同的断层斑块或流体驱动机制）还是震级标度伪影，仅凭目前的目录无法判定。作者提出，决定性的测试是与小规模地震群的矩震级（ $M_w$ ）或谱震级进行比较：如果转折在矩震级中依然存在，则物理性质的解释是合理的；如果转折消失，则它是震级标度的伪影。

启示
研究结果提醒，在解读该目录的 $b$ 值时，如果不指定震级范围和深度，应保持谨慎，因为有效 $b$ 值在转折点上下会发生剧烈变化。研究还表明，像 $Z_C$ 这样的离散度指标可以作为 FMD 形状变化的探测器，尽管在本例中，时间上的变化被解释为地震群在不同深度种群之间的迁移，而非前兆信号。

类似论文