利用GRACE资料构造大尺度时变重力场统一模型_涂弋_陆洋_张子占_史红岭_杜宗.pdfVIP

第 ４１ 卷 第 ８ 期 ２０１６ 年 ８ 月 武 汉 大 学 学 报 · 信 息 科 学 版 Ｇｅｏｍａｔｉｃｓ 　 ａｎｄ 　 Ｉｎｆｏｒｍａｔｉｏｎ 　 Ｓｃｉｅｎｃｅ 　 ｏｆ 　 Ｗｕｈａｎ 　 Ｕｎｉｖｅｒｓｉｔ ｙ Ｖｏｌ．４１Ｎｏ．８ Ａｕ ｇ ．２０１６ 收稿日期： ２０１６ － ０１ － ３０ 项目资助：国家 ９７３ 计划（ ２０１３ＣＢ７３３３０１ ）；国家自然科学基金（，）。 第一作者：涂弋，博士生，主要从事 ＧＲＡＣＥ 卫星应用研究。 ｔｕ ｙ ｉ０１２６ ＠ １６３．ｃｏｍ ＤＯＩ ： １０．１３２０３ ／ ｊ ．ｗｈｕ ｇ ｉ文章编号： １６７１ － ８８６０ （ ２０１６ ） ０８ － １１００ － ０７ 利用 ＧＲＡＣＥ 资料构造大尺度时变重力场统一模型 涂 　 弋１，２ 　 陆 　 洋１，２ 　 张子占１，３ 　 史红岭１，４，５ 杜宗亮１，２ 　 高春春１，２ 　 朱传东１，２ １ 　 中国科学院测量与地球物理研究所大地测量与地球动力学国家重点实验室，湖北 武汉， ４３００７７ ２ 　 中国科学院大学，北京， １０００４９ ３ 　 美国得克萨斯大学空间研究中心，美国 得克萨斯州 奥斯汀， ７８７５９ ４ 　 地理空间信息工程国家测绘地理信息局重点实验室，北京， １００８３０ ５ 　 极地测绘科学国家测绘地理信息局重点实验室，湖北 武汉， ４３００７９ 摘 　 要：基于 ２００３ ～ ２０１２ 年的 ＧＲＡＣＥ 卫星重力资料，采用最小二 乘 拟 合 的 方 法，构 建 了 时 变 重 力 场 统 一 模 型 ＩＧＧ － ＴＶＧ２０１３ 。该模型以球谐系数的形式表达，在考虑趋势项和周期项等经验参数的基础上，还考虑了加 速度项和潮汐模型误差、大地 震 等 因 素 的 影 响。将 ＩＧＧ － ＴＶＧ２０１３ 模 型 与 ＧＲＡＣＥ 资 料 进 行 了 比 较 分 析，在 全球 ９２％ 以上的区域二者符合精度优于 ±１ｕ ｇ ａｌ ；利 用 该 模 型 外 推 预 测 了 ２０１３ 年 １ ～ ６ 月 的 重 力 场 变 化，结 果与 ＧＲＡＣＥ 实测数据符合较好。这表明 ＩＧＧ － ＴＶＧ２０１３ 模型不但能较好地描述重力场的连续时空变化，而 且具有一定的短期预测能力。 关键词： ＧＲＡＣＥ ；时变重力场统一模型；最小二乘拟合；球谐系数；加速度；潮汐模型误差；大地震 中图法分类号： Ｐ２３１．５ 　　　　　 文献标志码： Ａ 　　 ＧＲＡＣＥ （ Ｇｒａｖｉｔ ｙ 　 Ｒｅｃｏｖｅｒ ｙ 　 ａｎｄ 　 Ｃｌｉｍａｔｅ 　 Ｅｘ ｐ ｅｒｉ － ｍｅｎｔ ）卫星自 ２００２ 年 ３ 月发射至今，已对全球重力 场变化已进行了十多年的监测。目前， ＧＲＡＣＥ 卫 星已经广泛应用于陆地水循环及水储量变化，两极 冰盖和山岳冰川消融，海水质量变化，海平面变化 和大地震检测等研究领域［１］。相关机构如美国德 克萨 斯 大 学 空 间 研 究 中 心（ Ｃｅｎｔｅｒ 　 ｆｏｒ 　 Ｓ ｐ ａｃｅ 　 Ｒｅ － ｓｅａｒｃｈ ， Ｕｎｉｖｅｒｓｉｔ ｙ 　 ｏｆ 　 Ｔｅｘａｓ 　 ａｔ 　 Ａｕｓｔｉｎ ， ＣＳＲ ）、德国 地学中心（ Ｇｅｏ 　 Ｆｏｒｓｃｈｕｎ ｇ ｓ 　 Ｚｅｎｔｒｕｍ ， ＧＦＺ ）和美国 喷气动力实验室（ Ｊｅｔ 　 Ｐｒｏ ｐ ｕｌｓｉｏｎ 　 Ｌａｂｏｒａｔｏｒ ｙ ， ＪＰＬ ） 等发布的卫星重力时变资料具有较高的观测精度， 误差分析表明在周年时间尺度和几百 ｋｍ 的空间 尺度上， ＧＲＡＣＥ 可探测出小于 １ｃｍ 的陆地水和 海水储量 变 化［２］。总 之， ＧＲＡＣＥ 卫 星 能 够 提 供 大尺度的、连续的时变重力资料，这对于相关地学 研究具有重要的科学意义。 目前，重力场时变资料是基于引力位 球 谐 系 数表达的，一般按月间隔提供的离散模型。离散 模型仅能表示一段时间内的平均变化，不能连续 地描述地球重力场变化过程，也不能直接用于多 源动态数据融合。为了连续、动态地描述重力场 随时间的变化，需要用一种统一的函数形式将离 散的球谐系数联系起来，即在传统的描述三维空 间变化的重力场模型基础上，引入时间变量，构造 出能够同时描述重力场时间变化和空间变化的四 维结构模型［３］。该模型具有以下几个特点， ① 能 方便计算出任意时刻的重力场动态数据，并为多 源数据融合工作提供有效的输入值； ② 能同时计 算大地水准面、重力异常、垂线偏差和梯度等多种 重力场基本参量，对这些参量进行联合分析，方便 地提取出重力场时变信息，如趋势、季节变化和年 际变 化 等； ③ 可 用 于 模 拟 以 前 的 重 力 场 变 化 情 形，也