泄漏校正用于修正球谐截断和滤波引起的信号扩散、振幅衰减和周边污染。对于湖泊、封闭流域、冰盖边缘和海岸带等边界陡变区域，泄漏误差会显著影响 GRACE 反演的季节振幅和长期趋势。因此，滤波结果不能只看图像是否平滑，还必须通过尺度因子、周边信号扣除、流域统计和外部产品验证进行综合评估。
背景
GRACE Level-2 球谐解通常需要截断到有限阶数，并经过滤波处理。这一过程会改变真实信号的空间谱结构，使目标区信号向外扩散，同时也可能让周边区域信号混入目标区。前者称为外泄漏，后者称为内泄漏。
对于里海这类封闭水体，真实信号主要集中在水体范围内，边界清晰。如果不进行泄漏校正，区域平均 TWSA 可能被低估，长期趋势和季节振幅也会产生系统偏差。
基本原理
泄漏校正可以采用两步法：第一步，利用 GLDAS 等陆面水文模型估算目标区周边陆地水储量变化，并施加与 GRACE 相同的截断和滤波处理，得到周边信号对目标区的污染贡献。第二步，构造目标区内均匀 1 cm 的合成信号，经过同样处理链计算尺度因子。
\[SF=\frac{h_{input}}{h_{filtered}},\quad TWSA_{corrected}=SF\cdot(TWSA_{GRACE}-TWSA_{leakage\text{-}in})\]
处理流程
- 准备目标区掩膜，例如里海水体范围或流域边界。
- 计算原始 GRACE 滤波结果的区域平均序列。
- 利用 GLDAS 或其他水文模型估算周边陆地污染信号。
- 对模型场施加与 GRACE 完全一致的处理流程。
- 从 GRACE 区域序列中扣除周边泄入信号。
- 构造合成信号并计算尺度因子。
- 对区域序列进行尺度因子恢复。
- 与 Hydroweb 卫星测高、CSR Mascon 或其他独立产品进行一致性验证。
产品标签与配置示例
该页面中的方法不应被统一保存为 filtered。配置文件至少应保留方法名、产品标签、处理域、上游输入和输出对象。
method: leakage_correction
tag: CASPIAN_HSAF_SF
target_region: Caspian_Sea
background_model: GLDAS_Noah_2_1
scale_factor: 1.22
validation: Hydroweb_altimetry
metrics: annual_amplitude, trend, correlation, rmse效果与局限
能够明显改善小尺度目标区的振幅恢复和趋势估计。
依赖外部模型和区域掩膜；目标区边界、尺度因子构造方式和滤波流程一致性都会影响最终结果。
页面图件建议
内泄漏与外泄漏示意图。
两步法泄漏校正流程图。
里海校正前后空间分布对比图。
GRACE 与 Hydroweb 测高时间序列对比图。
年振幅、趋势、RMSE 指标表。
滤波解决的是噪声问题，泄漏校正解决的是区域信号恢复问题；二者应作为连续处理链而不是互相替代的方法。