基于空间域的GRACE时变重力场去噪方法的研究及应用.pdfVIP

·830· 中国地球物理2011 (23)地球重力场变化与在地学中应用 基于空间域的GRACE时变重力场去噪方法研究及应用 and ofthe Methodin Domain StudyapplicationFHtering Space BasedonGRACETimeVariableField Gravity 王秋革+陈超+杜劲松陈波 Chao+DuChenBo WangQiuge+Chen Jinsong 中国地质大学(武汉)地球物理与空间信息学院武汉430074 1 ’E-mail：qiuge06006 05@163．com；chenchao@cug．edu．cn 由于地球内部质量的重新分布导致了不同时间尺度的地球重力场变化，利用足够精度和时空分辨率的重力 观测可以了解地球时变重力场。自从2002年3月成功发射，GRACE卫星已经开始承担起了以月为时间尺度来 测量全球时变重力场的任务，目前测量的空间精度可以达到500千米。尤其是在陆地上，这种测量反应的是整 体的陆地水形态，包括陆表水(如湖泊、湿地以及河流等)，土壤湿度，地下水以及覆盖的积雪(GRamillien，2006)。 然而，GRACE每月观测资料得到的地球重力场模型本身会受卫星轨道误差、双星K波段测距误差、加速 度计测量误差以及卫星姿态测量误差等的影响，同时，GRACE月重力场模型给出的是截断到一定阶次的位系数， 导致了时变重力场产生截断误差，月重力场模型的高阶误差占主要成分，表现为随模型球谐系数阶数的增大而迅 速增加，高阶项参与计算给时变重力场带来了更大的误差。而且，科学家们感兴趣的通常不是某一点的表面密度 变化，而是某一区域乃至全球的质量变化。为了有效提取有效信号，必须对时变重力场进行去噪或者平滑。以往 的去噪方法均在球谐域内进行，本文尝试在空间域进行滤波器的设计，进而对GRACE时变重力场去噪。 1．方法简介 计算一个点的平滑值时，首先以该点为球冠极点，构建球冠坐标系，各个网格点的原始时变重力值根据球 坐标变换，利用已有GRACE时变重力场模型直接解算：其次，将球冠内各个节点的时变重力值与平滑核函数 相乘再累加即得球冠极点的平滑值。该方法的关键在于设计平滑滤波器，即平滑核函数。本文首先设计各向同 性高斯滤波器，与球谐域内的高斯滤波比较验证了本文方法的正确性，进而设计各向异性滤波器与自适应高斯 滤波器。该方法比较灵活，在研究局部区域时，可以设定平滑范围，以减小泄露效应，而且适用于全球。 2．应用及结果分析 最后，将该方法运用到美国堪萨斯州地下水储量变化的研究中，我q]q用GRACE卫星测得的时变重 储量的时间空间变化特征，然后再与堪萨斯西部实际原位井(2345口井)测量的水位变化数据进行对比。 研究结果呈现出一种明显的季节性特征：冬春季的地下水储量比较大而夏秋季的地下水储量达到最低值。 结果显示，由GRACE卫星数据计算得到的地下水的变化和实际的原位井测量的数据在发现两者在2003 年至2009年中符合较好。另外，计算结果中，两种方式得到的地下水储量变化很明显的均呈现出一种下 降趋势。结果体现了GRACE卫星用于干旱地区地下水量变化监测的能力，另一方面也说明了基于空间域 的GRACE时变重力场去噪方法的可行性和有效性。 学基金联合资助。 参考文献 J，eta1．Time oftheEarth’s andoceaniceffectsandtheir detection [1】Wahr variability gravity field：hydrological possible