杨帅， 郑伟,3*， 尹文杰， 刘杰

1 河南理工大学测绘与国土信息工程学院， 焦作河南 454003 2 中国空间技术研究院钱学森空间技术实验室， 北京 100094 3 辽宁工程技术大学测绘与地理科学学院， 阜新辽宁 123000

0 引言

自2002年重力恢复与气候实验(Gravity Recovery and Climate Experiment，GRACE)卫星发射以来，GRACE重力卫星提供了高精度观测地球时变重力场信号的新手段(Tapley et al., 2004).截止到目前，GRACE重力卫星已被广泛应用于陆地水储量(Terrestrial Water Storage，TWS)变化(Rodell et al., 2009; Seyoum and Milewski, 2016)、冰川变化及其质量平衡(Velicogna and Wahr, 2006a,b; Mohajerani et al., 2018)、全球平均海洋质量变化(Chen et al., 2019)、地震(Han et al., 2006; Chen et al., 2007; Zhou et al., 2018)等多领域研究.GRACE重力卫星在2017年10月停止运行，随着GRACE Follow-on重力卫星的成功发射，该卫星反演的地球重力场模型继续在地球科学领域发挥着极其重要的作用(郑伟等, 2014).

基于Wahr等(1998)的理论，地球表面的质量变化可以利用GRACE球谐系数解反演得到.然而，球谐系数中既含有随系数阶数而增加的随机误差，又包含与系统相关的系统误差.因此，许多研究者相继提出了空间平滑滤波(如高斯平滑滤波(Wahr et al., 1998)、维纳滤波(Sasgen et al., 2007)等)和去条带滤波(Swenson和Wahr(2006)的去条带方法、Chen等(2009)的P4M6方法等).这些滤波方法在消除误差影响的同时也会改变真实的物理信号.Chen等(2015)通过实验模拟得出截断及平滑300 km产生的泄露误差在南极洲可达80.95 Gt·a-1(实测值-150.00 Gt·a-1，模拟值-69.05 Gt·a-1)；Jin和Zou(2015)研究表明格陵兰岛冰盖质量损失的信号中，42.4%是由泄漏误差造成，对泄漏误差进行校正后，格陵兰岛的冰川融化速度从-15.58±0.93 Gt·a-2增至-26.19±1.67 Gt·a-2.因此，泄露误差是GRACE反演质量变化中的主要限制因素，校正泄露误差对提高反演结果的准确性至关重要.

目前，GRACE数据泄露误差常用的……