GRACE反演南极冰盖质量变化的Gauss的Wiener滤波比较.pdfVIP

GRACE 反演南极冰盖质量变化的Gauss 与W iener 滤波比较 1 2 1 1 2 3 李军海 ，文汉江 ，刘焕玲 ，朱广彬 ( 1 中国测绘科学研究院 北京 100830 ； 2 辽宁工程技术大学 辽宁 阜新 123000 ； 3 武汉大学测绘学院 湖北 武汉 430079) 摘要：利用UTCSR 2003 年 1 月到2008 年8 月间的GRACE Level-2 RL04 重力场模型估计 了南极冰盖质量变化。计算过程中分别采用Gauss 和Wiener 滤波两种平滑方法，分别采用 22 、43 和 65 个月时间段的重力场模型计算Wiener 滤波信号与噪声函数。发现其信号与噪 声函数与模型所选的时间段有关。因此在实际的计算过程中需要具体计算Wiener 滤波平滑 因子值。将65 个月GRACE 重力场模型计算得到的Wiener 滤波权值与Gauss 滤波权值比较， 发现其非常接近于平滑半径为 540km 高斯滤波权值。通过对西南极 Amundsen 和东南极 Enderby 地区冰盖质量等效水高变化率比较发现，采用两种不同的滤波方法所计算得到的相 同区域质量变化率基本一致。 关键词：GRACE，南极，Gauss 滤波，Wiener 滤波 1 引言 南极地区储存有地球上最大的冰川，因此了解南极冰盖质量变化的规律对研究全球平均 海平面变化、气候变化、水文变化以及大气动力学等问题非常重要。GRACE重力卫星计划由 NASA和DLR联合开发，旨在获取高精度地球重力场及其时变特征，它的误差分析表明，在季 节性或年际性时间尺度和几百公里或更大空间尺度上，GRACE可探测平均小于1cm的表层水变 [1] 化，这为我们探测南极冰盖质量的变化提供了可能 。 GRACE卫星对地球重力场的探测能力随着空间分辨率的提高逐渐降低，因此随着球谐 系数的增加，误差积累逐渐变快。为了得到更合理的时变重力场,应对时变重力场进行平滑 [2] 处理。自Jekeli1981年将平滑技术用于球谐函数的滤波器以来,该技术已广泛应用于地球 [3] 物理、大地测量等地学学科。Wahr等 在1998年提出了利用Gauss平滑函数处理该误差，不 [4] 过该方法需要对平滑半径进行预定义。I.SASGEN等 在2005年提出了利用Wiener平滑函数处 理该误差，与Gauss平滑不同的是该方法采用阶方差谱来分别建立信号与噪声模型，受信号 与噪声比的约束很小。 [5] 本文利用UTCSR发布的GRACE Level-2 RL04版本重力场位模型 ，使用2003年01月 -2008年08月共65个月 GRACE重力场位模型应用于南极冰盖质量变化反演，其中没有2004 [6] 年01月的模型。由于单独使用GRACE得到的C20 项精度较低，朱广彬等 利用主成分分析发 现SLR提供的C20 能够提高南极冰盖质量变化反演精度。因此本文采用SLR得到的C20 代替 GRACE的C20 以提高反演精度。 2 GRACE数据反演质量变化的基本原理 [3,7] 物质质量分布变化所引起的球谐系数变化为 ： 2p p r l cos ml C 1 2 r lm