《Geosciences》：Comparative Statistical Detection of Ionospheric GPS-TEC Anomalies Associated with the 2021 Haiti and 2022 Cyprus Earthquakes
Sanjoy Kumar Pal,
Kousik Nanda,
Soumen Sarkar,
Stelios M. Potirakis,
Masashi Hayakawa and
Sudipta Sasmal
编辑推荐：
为解决地震前兆信号难以从背景电离层波动中可靠分离的问题，研究人员对2021年海地Mw 7.2和2022年塞浦路斯Mw 6.6地震，利用位于地震孕育区（EPZ）内的IGS站GPS数据，提取了垂直总电子含量（VTEC）序列。通过高斯、稳健中位数绝对偏差（MAD）和分位数三类统计模型构建“发生次数”和“能量”指数，并经归一化和95%分位数阈值滤波后发现，两次地震前数周均存在显著的VTEC异常增强，验证了归一化结合稳健统计方法对于从背景波动中甄别地震相关离子层TEC异常的实用性。该研究为岩石圈-大气层-电离层耦合（LAIC）机制提供了新的定量证据，对深入理解地震-电离层效应及提升地震监测能力具有重要意义。
地震，这种破坏力巨大的自然灾害，其可靠的短期预报至今仍是悬而未决的科学难题。在探索地震预警信号的道路上，科学家们将目光投向了地球上空几百公里处的电离层。一个名为“岩石圈-大气层-电离层耦合（LAIC）”的理论框架认为，地震孕育过程中地下岩石的应力积累、微破裂和流体运移等过程，可能通过电磁、化学和动力学等复杂通道，将扰动向上传播，最终影响到高层大气的电离层，使其中的电子密度发生异常变化。这为通过监测电离层状态来探寻地震前兆提供了可能。然而，这条路充满挑战：电离层本身就像一个“情绪”极其不稳定的高空海洋，太阳辐射、地磁活动等“天气”变化随时可能激起惊涛骇浪，其扰动幅度常常不亚于甚至超过地震可能引发的变化。因此，如何从这片充满“噪音”的海洋中，准确识别出与地震孕育相关的微弱“信号”，成为该领域研究的核心瓶颈。
为了攻克这一难题，并深入探究特定地震事件前的电离层响应，一篇发表在《Geosciences》上的研究，对两次中强震——2021年海地Mw 7.2地震和2022年塞浦路斯Mw 6.6地震——展开了详尽的对比分析。研究人员巧妙地选择了这两次事件，因为它们发生前数周没有其他强震干扰，且地磁活动相对平静，这为“纯净”地分析可能的地震-电离层信号创造了条件。研究的目标是，利用全球定位系统（GPS）衍生的垂直总电子含量（VTEC）数据，运用多种统计“武器”，去捕捉和验证震前是否存在具有统计显著性的电离层异常。
这项研究的技术路线清晰而系统。首先，数据来源于国际GNSS服务（IGS）网络中位于两次地震震中约300公里内的测站（海地：SCUB站；塞浦路斯：MERS站），这些站点均位于根据多贝罗夫斯基公式估算的地震孕育区（EPZ）内。研究人员处理了每个震例前后共53天（震前45天，震后7天）的双频GPS观测数据，使用Gopi Seemala软件计算VTEC，并通过扣除日变化分量得到了残差序列，以隔离短期的电离层变化。核心的异常检测采用了三套统计框架进行对比：基于均值和标准差的高斯模型、基于中位数和中位数绝对偏差（MAD）的稳健模型，以及无需分布假设的分位数模型。为了量化异常，他们构建了每日“发生次数”和“能量”指数，但发现直接使用这些指数会检测到过多异常，缺乏选择性。因此，研究引入了关键的后处理步骤：将所有指数除以其中位数进行归一化，然后仅保留超过95%分位数的极端值。这一步骤极大地压制了背景波动，凸显了少数统计上显著的异常峰值。
研究结果
1. 数据处理与背景表征
研究首先展示了地震位置、IGS测站及EPZ的范围。对地磁指数Dst和Kp的分析表明，两次地震前数周均处于地磁宁静期，这降低了空间天气效应主导观测到异常的可能性。
2. 多方法异常检测结果
无论采用高斯、MAD还是分位数方法，计算出的“发生次数”和“能量”指数在时间演变上均表现出相似的模式。然而，未经处理的指数普遍显示异常频发，统计学选择性低。在应用了归一化和95%分位数阈值过滤后，背景波动被大幅抑制，研究成功分离出了少量的显著性异常峰值。
3. 两次地震的异常特征对比
对于两次地震，增强的异常活动均在震前数周被识别出来。相比之下，震后的异常峰值则与升高的气象和地磁活动时期相吻合。这一时间关联性的差异，暗示震前检测到的异常可能具有不同的起源。结果证明，结合归一化的稳健统计方法，对于从背景变率中辨别出显著的离子层TEC异常至关重要。
结论与讨论
本项研究通过多统计方法对比分析，为2021年海地Mw 7.2和2022年塞浦路斯Mw 6.6地震前存在电离层VTEC异常提供了新的证据。研究表明，在严格的数据筛选（选择地磁宁静、无强震干扰的震例）和精心的数据处理（去除日变化、使用残差序列）基础上，结合归一化与高分位数阈值过滤的稳健统计方法，能够有效抑制背景电离层波动，提取出可能与地震孕育过程相关的显著性异常信号。
该研究的核心意义在于其方法论上的贡献。它强调了在复杂多变的电离层环境中，单一或简单的异常检测标准极易产生误报。通过引入并对比多种统计模型，并实施“归一化+极端值过滤”的后处理策略，研究建立了一套更严格、选择性更强的分析流程，提升了从背景“噪声”中识别地震相关“信号”的可靠性。这为未来基于GPS-TEC或其他电离层参数进行地震前兆的统计研究提供了重要的方法学参考。尽管两次震例的分析支持了LAIC机制的存在，但作者也谨慎指出，要确立电离层异常作为普适性的地震前兆，仍需在更大规模的地震目录上进行检验，并更严格地控制太阳和地磁活动的混杂影响。这项工作推动地震-电离层耦合研究向更定量、更稳健的统计验证方向迈进了一步。
