季节性水文应力与地球潮汐驱动下的加州地震节律：探索地震孕育过程的摩擦物理与灾害评估新途径

《SCIENCE ADVANCES》：Seismic rhythms: Earthquake response to tectonic, hydrological, and tidal forcing in California

地震活动的季节性调制揭示了地球物理过程的微妙平衡，但地震响应应力扰动的内在物理机制尚未完全阐明。为此，研究人员基于2006-2021年加州地震目录，利用先进的去丛集技术分离背景地震活动，并应用Schuster频谱分析，揭示了由水文地表荷载（季节性）而非固体地球潮汐（半日）驱动的显著地震活动调制现象。研究发现，地震活动率峰值比应力率峰值滞后约半月，经典的瞬时库仑破裂模型无法解释观测到的振幅和时间延迟，而基于速率-状态摩擦（RSF）定律的非瞬时孕震模型则能很好地拟合观测数据。该研究通过地震活动对周期性应力扰动的响应，定量揭示了控制断层摩擦性质的应力参数（aσ ~ 1-10 kPa）和特征松弛时间（ta~ 0.05-1 年），为利用应力变化探究断层力学特性、改进地震危险性评估开辟了新途径。

　　

在我们脚下，地球的“脉搏”——地震，并非随机跳动。驱动它们的主要力量是板块构造运动带来的缓慢应变积累，但越来越多的研究表明，一些微小、周期性的外界“压力”变化也能“拨动”地震的“琴弦”。这些外界压力包括地表水体、冰雪的季节性增减（水文荷载），以及由月球和太阳引力引起的固体地球周期性形变（地球潮汐）。尽管这些变化产生的应力仅有几千帕斯卡（kPa），相比地壳内部的岩石静压（约100-150兆帕）微不足道，却能在多个地区引发地震活动率的显著季节性波动。然而，为什么在一些地区季节调制信号很强，而同等量级的潮汐应力调制却不显著？地震活动响应这些应力变化的具体机制是什么？回答这些问题，不仅能帮助我们理解地震是如何“被触发”的，更关键的是，能借此揭示控制地震孕育过程的关键物理参数，如断层的摩擦状态和地下的绝对应力水平，从而为改进地震预测和危险性评估提供物理基础。为此，一组研究人员聚焦于拥有高质量地震目录和显著季节性荷载变化的美国加州，开展了一项深入分析，相关成果发表在了《SCIENCE ADVANCES》上。

研究团队运用了几项关键技术方法。首先，他们利用改进的最近邻距离法，对加州2006-2021年间的完整地震目录进行了去丛集分析，有效分离了由构造应力直接驱动的背景地震和由余震、震群等构成的丛集地震。其次，采用Schuster谱分析方法，在空间网格点上评估了背景地震活动率在不同周期（特别是1年周期和半日潮汐周期）的调制振幅和相位，精确识别了具有统计显著性的季节信号。再者，基于GPS和GRACE卫星观测数据反演得到的水文表面荷载（等效水厚度）模型，结合弹性半空间点荷载解析解，计算了在5公里深度、不同预设断层（最优构造应力向断层、最优水力应力向断层、平行圣安德烈斯断层的垂直走滑断层）上引起的库仑应力变化。最后，将观测到的地震响应振幅和时间延迟，与两种孕震物理模型（CFM和RSM）的预测结果进行对比，反演出控制断层孕震过程的摩擦应力参数（aσ）和特征松弛时间（t_a）。

研究结果

1. 长期背景地震活动由构造加载驱动

去丛集后的背景地震活动集中分布在已知的活动断层带上，并与大地测量获取的构造应变率呈良好的线性正相关关系。这表明，在过滤掉余震和震群后，剩余的背景地震活动主要由长期构造应变加载所驱动，为后续分析周期性的水文和潮汐调制提供了纯净的样本。

2. 显著的季节性调制，但无显著潮汐调制

空间分析表明，加州多个区域（特别是北加州）的背景地震活动率存在显著的、空间一致的年度（季节性）调制。地震活动率的峰值多出现在7月至8月。调制振幅在局部区域可达背景地震活动率的15%。然而，在相同的空间尺度上，未检测到具有统计显著性的半日潮汐调制信号。这表明地震活动对谐波应力扰动的响应具有周期依赖性。

3. 地震活动率峰值滞后于应力率峰值

比较观测到的地震活动季节性峰值时间与水文荷载引起的库仑应力率峰值时间，发现地震活动响应存在明显的滞后，中位数时间延迟约为0.52个月（约半个月）。这种滞后是瞬时破裂模型无法解释的，暗示地震孕育过程并非瞬时完成。

4. 用物理模型解释地震活动调制

研究人员首先检验了库仑破裂模型（CFM）。该模型假设地震在应力超过断层强度时瞬时成核。结果表明，无论临界周期T_τ如何设定，CFM都严重高估了地震响应振幅，且无法解释观测到的时间延迟。随后，他们采用了基于速率-状态摩擦（rate-and-state friction, RSF）定律的孕震模型（RSM）。RSM考虑了摩擦过程的时间依赖性，其响应由两个关键参数控制：摩擦应力参数aσ（乘积）主要控制响应振幅（aσ越小，响应越大），特征松弛时间t_a主要控制响应持续时间。通过将观测到的地震调制振幅和相位滞后与RSM预测进行匹配，研究人员成功反演出了加州大部分地区的aσ值（主要集中在1-10 kPa之间）和t_a值（主要集中在0.05-1年之间）。

5. 地震敏感性随时间变化

以帕克菲尔德地区为例，研究发现该地区在2004年M6.0地震之前（1984-2003年）存在清晰的季节性调制，但在该地震之后（2006-2021年）的观测时段内，调制信号却消失了。这提示断层带的应力状态或摩擦属性可能被大地震所改变，从而影响了地震活动对应力扰动的敏感性，与实验室观测到的应力依赖性敏感性一致。

结论与意义

这项研究通过系统分析加州地震活动对水文季节性荷载和固体地球潮汐的响应，揭示了地震孕育过程的非瞬时特性，并成功利用观测数据约束了断层摩擦的物理参数。研究得出的aσ值（~1-10 kPa）与从慢滑移驱动震群、诱发地震以及余震序列中得到的估计值范围一致。若取典型的实验室摩擦参数a=0.001，则意味着断层有效正应力σ仅为几兆帕，远低于5公里深度约150 MPa的岩石静压，这可能表明天然断层的a值比实验室条件下低得多，和/或断层内流体压力极高，导致有效正应力极低。

该研究的核心意义在于，它提供了一种“被动”探测地球地壳应力状态和断层力学性质的新方法。通过分析地震活动对已知的、周期性的微小应力扰动（如季节性水文荷载、地球潮汐）的响应，可以在不进行主动干预的情况下，定量估算控制地震成核的关键摩擦参数。这为评估区域断层带的“脆弱性”和应力状态开辟了新途径。更重要的是，这种方法可以应用于与可持续能源生产（如地热开发、碳封存）相关的诱发地震危险性评估中。在那些断层应力和摩擦性质未知的区域，观测地震活动是否存在对潮汐或季节应力的响应，可以为建模和缓解由人类活动（流体注入/抽取）引发的地震风险提供宝贵的约束信息。尽管当前采用的RSM模型未考虑有限断层效应等复杂因素，但它作为一个高效的半解析工具，首次实现了在广阔空间范围内对孕震参数的一阶估计，为未来的精细化数值模拟和预测奠定了重要基础。