Seasonal Adaptive Analysis of Radon Data and Correlation with Regional Earthquake Events

技术摘要：拉东数据的季节性自适应分析及其与区域地震事件的相关性

问题陈述
本研究旨在解决识别大气中拉东（Radon）排放中瞬态地球化学异常的挑战，这些异常可能作为区域地震事件的前兆。该领域的主要难点在于如何将真实的构造信号与多尺度环境波动（包括日间气象噪声，如温度、湿度；以及中期的天气系统因素，如气压、土壤湿度）区分开来。本文认为，若缺乏稳健的统计框架来过滤这些人为干扰，地壳脱气作用与地震加速之间的相关性将被掩盖。

方法论
研究实现了一种利用“全局标准差（Global Standard Deviation, SD）基准”的“季节性自适应分析”框架。

• 数据采集： 在海拔540米处（北纬 35.30°，东经 25.64°）通过商用 Airthings 硅光电二极管检测器进行连续拉东监测。传感器被保护在专门的封装结构内，以防止水分侵入和直接的环境辐射干扰，并通过 $\alpha$ 能谱法进行工作。
• 数据处理： 为了抑制日间噪声，高保真原始时间序列数据被聚合为日平均值。
• 统计框架：
  • 计算了整个时间序列背景的全局标准差 ($\sigma_{global}$)，以建立一个能够过滤中短期气象调制效应的统计基准。
  • 定义了一个自适应警报阈值：$Threshold = \mu_{6-day} + 2.5 \times \sigma_{global}$。
  • 该配置利用 6 天过去均值 ($\mu_{6-day}$) 来追踪短期趋势，同时利用全局标准差确保对快速、应力诱导的地壳脱气事件具有敏感性，并最大限度地减少来自局部环境噪声的误报。
• 相关性： 本研究将此地球化学时间序列与 2026 年 6 月 24 日至 7 月 28 日期间的区域地震活动数据（震级 $M_L$）相结合。

关键结果
对 2026 年 6 月至 7 月数据集的分析得出了以下观察结果：

• 统计参数： 在 35 天内，拉东平均值为 20.3 Bq/m³，标准差为 11.3，范围从最小值 7 Bq/m³ 到最大值 56 Bq/m³。
• 已验证异常（震前）： 2026 年 6 月 29 日检测到一个显著异常，偏差约为 2.9 SD。该事件发生在 2026 年 7 月 2 日 M 5.2 大地震发生前约 72 小时。研究将其解释为由临界地壳应变加速和微裂隙作用引起的短期地球化学前兆。
• 震后响应： 在 M 5.2 地震后，7 月 3 日出现了一个明显的亚阈值峰值，这归因于破裂后的应力松弛以及迫使深层拉东进入大气的瞬态水文地质扰动。
• 亚阈值前兆： 识别出两个中间峰值，分别为 **7 月 9 日（37 Bq/m³）**和 **7 月 11 日（33 Bq/m³）**。虽然它们未突破严格的 +2.5 SD 阈值，但其在时间上靠近二次地震簇（7 月 16 日 M 3.8 和 7 月 18 日 M 4.4），表明区域构造应力呈现阶梯式累积。
• 当前警报： 2026 年 7 月 28 日记录了一个关键的 2.9 SD 异常。作者将其归类为未经验证但高优先级的警报，并指出其统计行为与经过验证的 6 月 29 日前兆高度相似。

意义与主张
本文声称，通过将低成本、商用硅光电二极管技术（特别是 Airthings 传感器）与稳健的统计处理及环境屏蔽相结合，可以实现高分辨率的地震前兆监测。

• 噪声解耦： 结合基于随机森林（Random Forest）的气象基准与独立的 6 日过去均值，成功实现了气象噪声与真实构造信号的解耦。
• 机制验证： 拉东脱气峰值与区域地震震级之间的时间对齐，支持了如下前提：局部构造应力积累会调节地震发生前的地下气体传输机制。
• 可扩展性： 研究断言，该方法验证了部署密集型、低成本区域预警网络的潜力，其异常解析能力可与昂贵的科学仪器相媲美。
• 即时应用： 作者得出结论，7 月 28 日的异常使得更广泛的克里特（Crete）地区处于“高度关键的监测窗口”，需要进行实时观察和紧急地球物理评估，以作为潜在的短期前兆信号。