大地震前，地球真的会“偷偷”活动吗？
地震预测，一直是人类与自然搏斗中最具挑战性的目标之一。几十年来，科学家们一直在努力寻找大地震发生前可能存在的蛛丝马迹，即所谓的“前兆信号”。其中，一个核心问题是：在断层最终破裂、引发灾难性地震的几小时前，它是否真的会以一种可以被仪器捕捉到的方式“偷偷”加速滑动，为后续的大爆发积蓄能量？这就是“预滑”或“前兆性断层滑移”假说。
2023年，一项发表在《科学》杂志上的研究（Bletery 和 Nocquet, 2023）在全球范围内投下了一枚重磅“炸弹”。研究者们分析了90个大地震发生前两天的GNSS（全球导航卫星系统）数据，声称通过一种叫做“点积叠加”的方法，在震前数小时普遍观测到了断层滑移呈现“指数加速”的迹象，就像一颗不断加速的、最终冲向终点的弹珠。这一发现让许多人看到了短期地震预测的新曙光。
然而，科学界并未被轻易说服。紧随其后的几项研究给出了不同的声音。例如，Hirose等人(2024)在分析日本东北大地震的倾斜仪数据时，并未发现类似的加速现象。更有研究指出，BN23的结果可能受到GNSS数据中普遍存在的“共模噪声”影响，一旦将这些背景噪音剔除，加速现象便消失了。这些争议使得“预滑”假说的真实性再度蒙上阴影。
更重要的是，BN23研究的样本存在一个显著的“盲区”：他们主要分析了发生在板块边界俯冲带的大型地震，但对于发生在板块内部（如中国大陆大部分地区）的地震，则缺乏足够的考察。那么，一个关键问题悬而未决：BN23发现的所谓“前兆加速”现象，是地震的普遍规律，还是仅仅是板块边界俯冲带地震的“专利”？对于地震频发的中国大陆而言，搞清楚这个问题，对于评估本地的地震风险和理解地震孕育机制具有至关重要的意义。
为了填补这一知识空白，西北工业大学的张成阳、贾柯、王浩浩、郝明组成的研究团队，将目光投向了中国大陆地区。他们决定用同样的“尺子”——即BN23提出的研究方法，来丈量我们脚下的这片土地。他们的核心问题是：在2010年至2024年间，发生在中国大陆内部的那些震级（M）不小于6.5的强烈地震，在发生前是否也普遍存在可检测到的地壳形变加速迹象？这项研究的结果，以一篇题为《Is there accelerated crustal deformation before large earthquakes? A case study of strong intra-plate earthquakes in China》的论文，发表在《Earthquake Science》期刊上，为我们带来了来自大陆内部的答案。
关键技术方法
为了回答上述问题，研究人员系统性地开展了一项回顾性案例研究。他们首先从中国地震台网中心选取了2010年至2024年间发生在中国大陆的15个M≥6.5地震（不包括2014年于田地震）。1 利用以震中为中心、半径500公里范围内的中国大陆构造环境监测网络(CMONOC)连续GNSS观测站数据，获取了每个地震前两天的观测序列。数据处理是本研究的关键环节，团队使用了PRIDE PPP-AR软件进行动态精密单点定位(PPP)，以获取高精度的测站位置时间序列。在剔除异常值后，他们将数据采样间隔从30秒降采样至5分钟，以匹配BN23的研究参数。核心的分析方法是严格复现BN23提出的S函数计算方法：首先计算每个GNSS测站记录的水平位移与基于断层模型计算出的格林函数(Green’s function)的点积，以量化观测信号与理论断层滑动方向的吻合程度；然后，通过对所有可用测站的归一化点积结果进行叠加，生成表征每个地震“前兆信号强度”随时间变化的S函数。通过对单个地震以及分组、整体地震的S函数进行拟合（包括正弦函数、指数函数拟合和移动平均分析），来判断是否存在周期性和震前加速现象。
研究结果：来自中国大陆的否定答案
1. 单个地震的S函数：大部分“静悄悄”
研究团队逐一计算了15个地震的S函数，并根据其震前趋势将它们分为三类：上升型、下降型和无明显趋势型。结果显示：
• 仅1例上升 ：2024年新疆乌什M7.1地震的S函数在震前4小时内表现出上升趋势，指数拟合时间常数为1.65小时，接近BN23报告的1.3小时。但进一步分析发现，这一结果严重依赖于距离震源最近的WUSH站，当剔除该站数据后，上升趋势消失。
• 仅1例下降 ：2022年四川泸定M6.8地震的S函数在震前约6小时开始显著下降。
• 13例无明显趋势 ：其余13个地震（包括2021年青海玛多M7.4、2017年四川九寨沟M7.0等）的S函数在震前没有显示出清晰的、可与其他时间段明显区分的上升或下降趋势，其波动相对平稳。例如，对2021年玛多地震的分析显示，尽管其1小时50分钟的移动平均值在震前8小时有微弱上浮，但数值仍低于更早时期的多数S函数值，指数拟合甚至得到了负系数，未发现有效的预滑加速证据。
2. 叠加结果：依然缺乏普遍加速证据
为了增强信号、探索普遍规律，研究人员将所有15个地震的S函数进行了叠加分析。结果显示：
• 无整体加速迹象 ：综合S函数在震前没有显示出加速上升的现象。相反，指数函数拟合表明叠加的S函数在震前有轻微下降趋势，拟合得到的时间常数为4.52小时，远大于BN23报告的1.3小时。
• 存在周期性 ：叠加S函数显示出了约11.49小时的周期性，这与BN23报告中观察到的12.9小时周期现象有相似之处。尽管数据已进行潮汐改正，但该周期性现象的根源尚不确定。
3. 按断层机制分类：加速现象仍未显现
为了探究不同断层类型的地震是否有不同表现，研究人员将地震按震源机制分为逆断层、走滑断层等类别，并分别进行叠加分析。结果显示：
• 逆断层地震 ：6个逆断层地震的叠加S函数表现出周期性特征，但未发现明显的预滑加速。
• 走滑断层地震 ：7个走滑断层地震的叠加S函数同样显示出周期性，但也缺乏加速证据。
由于其他类型地震数量太少，其叠加结果不具有统计意义。这些分类分析进一步表明，在板内地震中，不同断层类型的地震均未普遍表现出BN23所报告的指数加速前兆。
结论与讨论：前兆加速并非“标配”，方法敏感性需警惕
本研究的主要结论是：在2010-2024年间中国内陆发生的15个M≥6.5强板内地震中，没有发现普遍的、清晰的震前地壳形变加速（预滑）现象 。单个地震分析中，仅有个别案例显示出趋势，且易受单个临近台站数据的支配。综合所有地震或按断层类型分组叠加的结果，也未能证实BN23所报告的指数加速前兆的普适性。这提示，断层预滑加速可能并非大地震发生前的“标配”现象，其出现或许与特定的地震类型（如巨型俯冲带地震）、区域构造环境或数据条件有关。
研究在讨论部分深入剖析了导致与BN23结果差异的可能原因及S函数方法的局限性：
• 台站分布与数据质量的重大影响 ：分析表明，S函数结果对GNSS台站的分布极为敏感。靠近震中的台站由于信号强、格林函数计算的静态位移大，其点积结果会不成比例地主导整个地震的S函数。在中国大陆，GNSS台网密度和分布不均，且数据质量参差不齐，这导致单个地震的S函数容易被少数甚至单个台站“绑架”，结果的稳健性存疑。例如，乌什地震的“前兆”信号在剔除最近台站后即消失。
• 震源机制解的不确定性 ：格林函数的计算依赖于震源机制解（走向、倾角、滑动角）。测试表明，这些角度在20°的不确定度范围内变化时，虽然不改变S函数的总体上升或下降趋势，但会影响其具体形态，带来额外的不确定性。
• 板内与俯冲带地震的潜在差异 ：BN23的研究主要针对板块边界俯冲带上的大地震，而本研究聚焦于板内地震。两者在构造背景、应力加载方式、断层规模等方面可能存在本质差异，这可能是导致观测现象不同的根本原因之一。
• 周期性现象的谜团 ：一个有趣的发现是，尽管缺乏加速证据，但本研究和BN23的研究都在叠加的S函数中观察到了约12小时左右的周期性波动。在已进行潮汐改正的前提下，这种现象的来源仍不明确，其是否蕴含了与地震孕育相关的物理信息，值得在未来利用更庞大的数据集进行深入探究。
重要意义
这项研究的意义在于，它首次系统性地将引发广泛关注的“震前形变加速”假说检验场，从板块边界延伸到了中国大陆这样的板内地震活跃区，并得出了基本否定的答案。这为评估“预滑”作为地震普适性前兆的可靠性提供了关键且必要的补充案例。研究表明，地震前兆现象可能比想象中更为复杂和多样，不存在“放之四海而皆准”的单一模式。同时，研究也警示，基于GNSS点积叠加的S函数分析方法虽然有力，但其结果严重依赖于台网密度、数据质量和数据处理细节，在应用于台站稀疏或分布不均的区域时需要格外谨慎。该工作强调了在探讨地震前兆时，必须考虑地震类型、构造背景的差异性，并需要整合多种数据源进行交叉验证。未来，通过对更长时间尺度、更多地震案例，尤其是不同构造环境下的大地震进行归类分析，或许能更清晰地揭示前兆信号的真相，为最终理解地震孕育过程和提升灾害风险评估能力累积宝贵的知识。
打赏
