《Geomatics, Natural Hazards and Risk》：An integrated framework for decoupling geodetic parameters in early seismic deformation monitoring: a case study of the 2011 MW 9.0 Tohoku earthquake
编辑推荐：
准确利用全球导航卫星系统（GNSS）进行地震早期形变的实时监测长期受制于多类大地测量参数之间的耦合效应，导致在地震初始阶段产生虚假的同震信号。本研究提出了一种由模型驱动的全解耦框架，可将卫星模糊度、电离层与对流层延迟、钟差以及位置估计进行分离。将该框架应用于地
准确利用全球导航卫星系统（GNSS）进行地震早期形变的实时监测长期受制于多类大地测量参数之间的耦合效应，导致在地震初始阶段产生虚假的同震信号。本研究提出了一种由模型驱动的全解耦框架，可将卫星模糊度、电离层与对流层延迟、钟差以及位置估计进行分离。将该框架应用于地震事件早期阶段的分析表明，同震形变中的厘米级波动误差被降低至毫米量级，精度接近地震加速度计测量结果。即便在应用区域增强（RA）及模糊度固定（AR）之后，研究人员仍在地震早期阶段发现幅度高达3 cm的虚假同震波形。这一现象源于显著的参数耦合效应，其中电离层延迟与垂直方向位置及接收机钟差参数之间的相关系数的中位数绝对偏差（MAD）值超过0.26。通过完全解耦这些参数，该解耦框架有效校正了由电离层误差引入的形变伪影，从而增强了主震期间地震与电离层过程的探测能力。
本研究发表于《Geomatics, Natural Hazards and Risk》期刊，旨在解决GNSS地震监测中参数耦合这一长期存在的难题。传统上，地震形变通过积分加速度计信号或宽频带地震仪速度来推断，但这些仪器在大震级地震（M > 7）中易发生饱和或削波，导致形变估计出现严重偏差。GNSS技术因其不受饱和效应影响而成为一种变革性手段，可直接估计地面形变，对地震和海啸早期预警至关重要。然而，Chanard等人曾警告需警惕对GNSS时间序列中季节性信号的过度解读，指出系统误差、环境影响及GNSS处理中的参数相关性可能导致非真实的周期性变化。特别是在地震发生后的最初数秒，载波相位模糊度、对流层延迟、电离层延迟和接收机钟差等参数之间的强相关性会显著降低推断形变信号的可靠性。尽管模糊度固定技术和区域增强策略可提升定位精度，但研究人员发现即便应用这些技术后，垂直方向位置与对流层延迟、接收机钟差之间仍存在显著耦合，限制了GNSS形变估计的精度与稳健性。因此，构建能够系统解耦所有关键大地测量参数的框架具有重要科学价值。
研究人员利用2011年M
W 9.0东日本大地震前地震平静期的数据，选取国际GNSS服务组织（IGS）测站USUD作为主要监测站开展研究。该站配备高稳定度氢原子钟（H-maser），且邻近区域分布有密集的GEONET参考站网络，同时设有同址地震台站NGN014可用于独立验证。研究设计并比较了四种处理策略：标准浮点模糊度解（PPP）、模糊度固定解（PPP-AR）、附加对流层和电离层区域增强的PPP-AR解（PPP-AR + Tro and Ion RA），以及进一步引入接收机钟差建模的完整解耦框架（PPP-AR + Tro and Ion RA + RCM）。
**模糊度固定对参数解耦的影响** 通过对比PPP浮点解与PPP-AR解，研究人员发现初始收敛阶段浮点模糊度与位置参数存在显著相关性，相关系数常超过0.5；而成功固定模糊度后，相关系数急剧收敛至接近零，实现了模糊度参数与位置估计的结构解耦。定位精度方面，PPP-AR解在N、E、U三个方向的均方根（RMS）误差分别从1.6 cm、2.6 cm、6.4 cm降至0.8 cm、0.9 cm、4.1 cm，E方向改善尤为显著。
**区域增强的参数相关性与解耦分析** 研究采用改进的线性插值方法进行区域增强，包括对流层区域增强（Tro RA）和电离层区域增强（Ion RA）。结果显示，Tro RA对U分量有局部去相关作用，但对N、E分量影响甚微；Ion RA显著提升了定位收敛速度，将首次固定时间（TTFF）从数十分钟缩短至数分钟，并将电离层参数与ENU分量的相关系数降至0.35以下。然而，即便应用RA后，垂直方向定位精度仍劣于水平方向。
**高精度接收机钟差建模以增强参数解耦** 针对垂直方向与接收机钟差之间超过0.5的相关性，研究提出基于随机游走、闪烁噪声和白频率噪声功率谱密度（s
?2、s
?1、s
0）的接收机钟差建模（RCM）方法。通过将RCM纳入前向卡尔曼滤波，N分量相关系数高达0.5的波动被完全消除，U方向相关性从0.5以上降至0.25以下。RCM未影响RA提供的收敛性能，但显著改善了收敛后的垂直定位精度。
**全参数解耦框架及其在地震早期形变监测中的应用** 研究建立了统一数学表述，将AR、RA和RCM集成于共同GNSS观测模型中。收敛性能方面，浮点PPP的U分量收敛时间为84.0分钟，PPP-AR降至56.3分钟，而加入RA后N、E、U分量初始收敛时间分别降至21.7分钟、0分钟和15.0分钟；对于重收敛情形，RA辅助策略在USUD站无需重新收敛。定位精度方面，四种策略在N、E、U方向的RMS依次为：浮点PPP（2.4 cm、4.0 cm、6.7 cm）、PPP-AR（1.1 cm、1.1 cm、5.1 cm）、PPP-AR + Tro and Ion RA（0.8 cm、0.9 cm、3.8 cm）、以及完整解耦框架（0.8 cm、0.9 cm、3.2 cm）。
在地震早期波形分析中，研究人员将1 Hz GNSS数据与NGN014台站100 Hz强震加速度数据经二次积分得到的地震仪波形进行对比。结果显示，PPP-AR解在震前平静期（05:47:30–05:48:00）呈现约+5 cm的虚假垂直波动；引入Ion RA后，虚假波形幅度降至+3 cm但仍具相似形态；而完整解耦框架（PPP-AR + Tro and Ion RA + RCM）显著消除了这些厘米级伪影。与地震仪波形的差异标准差（STD）表明，U分量精度从PPP-AR的1.36 cm，经RA改善至1.06 cm，最终通过RCM进一步降至0.68 cm，总计改善约50%。
针对关键相关系数的深入分析揭示，RA情形下电离层参数与垂直位置和接收机钟差的MAD值分别为0.31和0.26；引入RCM后降至0.13和0.03，证实全解耦框架有效削弱了电离层延迟与其他参数间的耦合。研究表明，当参数解耦不完全时，外部电离层增强信息会干扰GNSS解算的内部动力学，导致虚假垂直位移；而完整解耦后，即便存在增强数据，电离层延迟估计仍保持稳健。
研究讨论部分指出，该全解耦框架的实施需要严格的数据条件：密集区域GNSS网络支持RA、配备高稳定度原子钟的测站实现RCM、以及同址地震台站提供独立验证。尽管这些条件限制了框架的直接适用范围，但研究已证明RCM策略在日本和智利等多站点多事件中的可行性，并将此前研究中的后处理模式扩展至实时前向滤波处理。此外，研究回应了Chanard等人关于谨慎解读GNSS季节性信号的关切，建议将全参数解耦作为信号解释与分离的核心组件，通过结构性分离位置、钟差、硬件偏差和大气参数，为区分真实地表形变与处理伪影提供更符合物理一致性的基础。研究还建议对大气增强信息和未校准相位延迟（UPD）等外部输入实施更严格的质量控制，以避免非物理信号的传播。讨论最后展望，该框架不仅可用于隔离关键地球物理参数，还可评估此类外部信息增强的内部一致性，为探测微弱季节或瞬态信号提供可靠保障。
**研究结论** 本研究提出了一种集成模型驱动框架，有效解决了GNSS地震监测中长期存在的参数耦合问题。引入模糊度固定（AR）和区域增强（RA）可如预期改善GNSS定位精度和收敛性，但即便收敛后，垂直方向定位精度仍低于水平分量，主要源于垂直方向位置与接收机钟差之间的强相关性（>0.5）。为此，研究人员进一步应用接收机钟差建模（RCM）进行参数解耦。RCM使N分量相关系数高达0.5的波动完全消除，U方向相关性从0.5降至0.25。通过分离卫星模糊度、大气延迟、钟差和位置估计，该框架消除了通常会扭曲大地震早期阶段的虚假同震信号。特别地，当应用电离层RA时，GNSS导出的垂直位移波形在震前平静期和地震初始阶段呈现高达3 cm的错误幅度。引入RCM进入全解耦框架后，电离层参数与垂直位置分量和接收机钟差参数之间相关系数的MAD值分别从0.31和0.26降至0.13和0.03，表明参数耦合显著降低。由此，全解耦框架有效校正了强地面运动期间由电离层误差引起的伪影，从而实现了地震和电离层过程的稳健联合探测。
