杨帆^１，２，许厚泽^１，２，钟敏^２，王长青^２，周泽兵^１
１ 华中科技大学物理学院地球物理研究所，武汉　430074
２ 中国科学院测量与地球物理研究所大地测量与地球动力学国家重点实验室，武汉　430077

摘要　本文利用GRACE（Gravity Recovery And Climate Experiment）level 1b数据和径向基函数RBF（radial basis function）方法解算了全球时变地球重力场．RBF基函数相比传统球谐（spherical harmonic）基函数，其高度的空域局部特性使得正则化过程易于添加先验协方差信息，从而可能揭示更加准确的重力场信号．本文研究表明，RBF 基函数算法在精化现有的GRACE 全球时变重力场模型，如提升部分区域信号幅度等方面具有一定优势．本文通过将RBF的尺度因子作为待解参数，基于GRACE 卫星的level 1b数据和变分方程法，成功获取了2009-2010 年90阶无约束全球时变重力场RBF 模型Hust-IGG03，以及正则化全球时变重力场RBF 模型Hust-IGG04．通过与GRACE官方数据处理中心GFZ发布的最新90阶球谐基时变模型RL05a进行对比，结果表明：（１）无约束RBF模型Hust-IGG03和GFZ RL05a在空域和频域表现基本一致；（２）正则化RBF模型Hust-IGG04无需进行后处理滤波已经显示较高信噪比，噪音水平接近于球谐基模型GFZ RL05a经400km 高斯滤波后的效果；（３）Hust-IGG04相比400km 高斯滤波GFZ RL05a在周年振幅图和趋势图上显示出更多的细节信息，并且呈现出更强的信号幅度，如在格陵兰冰川融化趋势估计上Hust-IGG04比GFZ RL05a提高了24.2%．以上结果均显示RBF方法有助于进一步挖掘GRACE 观测值所包含的时变重力场信息．