利用GRACE反演长江流域水储量变化

利用GRACE重力卫星2003年1月～2006年9月共计43个月的时变重力场数据,反演了长江流域水储量的变化.结果显示,基于GRACE数据的反演结果和CPC模型符合得相当好;若考虑低阶项的影响,对GRACE的反演结果有很明显的改进.     

GRACE时变重力位系数扇形滤波算法分析

卫星重力探测技术为监测全球陆地水储量变化提供了新的技术手段。采用Level-2 Release-05版本GRACE时变重力场模型数据计算了2010年全球陆地水储量的月变化;着重研究了扇形滤波对反演结果的影响;并结合GLDAS水文模型数据对反演结果进行了验证分析。实验结果表明:GRACE反演结果 GLDAS水文模型结果在时空分布上符合较好;扇形滤波能够削弱GRACE时变重力场模型的高阶项误差影响,有效去除反演结果中的条带状噪声。     

利用GRACE重力卫星反演2003~2013年新疆天山地区陆地水储量时空变化

利用2003-1~2013-12的GRACE月重力场数据来计算新疆天山地区陆地水储量的变化,采用去相关滤波和300 km高斯滤波相结合的方法来滤除GRACE数据中的噪声,同时采用尺度因子的方法来减小GRACE后处理误差的影响,并利用Paulson模型结果来扣除由冰川均衡调整（GIA）对水储量反演结果的影响,把得到的最终结果与同期GLDAS、CPC水文模型进行了验证分析。结果表明,GRACE反演结果与两个水文模型的模拟结果变化趋势基本一致,在2003~2013年间,研究区域的水储量整体上呈现下降趋势,速率为-0.54±0.27 mm/a左右,但期间水储量变化波动较大,在2008年10月份左右,该区域水储量较同期呈现明显减小,这与该时段的干旱事件相一致。     

GRACE离散时变重力场的大尺度动态变化模型构建分析

针对GRACE月时变重力场的离散球谐系数,提出联合经验正交函数分解和多通道奇异谱分析的方法,在有效实现空间维数压缩和条带噪声去除的同时,提取、分析其动态变化特征并构建相应的变化模型。结果表明,GRACE离散球谐系数存在显著的周年项和长期趋势项变化成分,其中,周年项模态的方差贡献达到50.6%,扣除周年项成分后长期趋势项的累积方差贡献达到77.4%。以长期趋势项和周年项成分构建的非线性模型能较好地反映离散球谐系数及其反演结果的动态变化特征,其对全球陆表等效水柱高的拟合标准差介于0.3~14.0 cm,平均标准差为1.8 cm。研究结果可为分析地球重力场各参量变化特征及地球动力学变化提供重要参考。     

京津冀陆地水储量变化与区域降水相关性分析

针对京津冀地区水资源灾害频发、水储量监测手段存在局限性的问题,本文利用该地区2015年1月至2017年6月的重力场恢复与气候实验（GRACE）月重力场模型数据反演陆地水储量变化信息,结合区域降水、全球陆地资料同化系统（GLDAS）模型数据,对研究区域内水储量变化时空特征进行相关性分析。结果表明：GRACE监测的陆地水储量变化滞后降水2个月达到强相关性,监测到2016年春旱灾害的发生,扣除季节项影响后,GRACE&GLDAS相关系数为0.88,研究区陆地水储量上升速率为0.118 mm/a。研究结果可为该区域水储量时空调配、旱涝灾害防治提供理论依据。     

星星跟踪重力场测量任务优化设计方法

在卫星重力场测量中,星星跟踪是获取中高阶重力场模型的有效方式,是GRACE Follow-on、GRACE II等下一代国际重力卫星所采用的测量方式.星星跟踪重力卫星任务设计需要考虑轨道高度、星间距离、定轨误差、星间距离变化率测量误差、非引力干扰确定误差、任务测量时间和数据采样间隔等任务参数,这些参数共同决定了重力场测量的时间分辨率、空间分辨率及其精度等重力场测量性能.如何分析这些系统参数对重力场测量性能的复杂物理机理,进而提出合理、优化的任务参数设计方法,是星星跟踪重力场测量系统设计中的重要问题.为此,本文建立了星星跟踪重力场测量性能的解析计算模型,并利用GRACE重力卫星测量参数验证了该解析模型,进而提出了重力卫星系统参数设计方法,为实现星星跟踪重力场测量性能最大化奠定了理论基础.     

科尔沁沙地陆地水储量变化的GRACE时变重力场探测

为了实现科尔沁沙地的水资源监测,该文采用2003年7月至2010年12月的GRACE月重力场模型,经过去相关滤波、高斯平滑滤波和GIA改正,利用尺度因子法恢复重力场信号,反演科尔沁沙地的陆地水储量变化,与CPC水文模型反演结果进行对比分析。研究结果表明:由GRACE Release-05Level-2数据反演得到的2003年7月至2010年12月科尔沁沙地陆地水储量下降速率为-13.2±2.6(mm·a^-1);CPC水文模型反演的该地区陆地水储量变化曲线与GRACE反演结果呈现出良好的一致性,具有相似的季节变化;2009年秋至2010年春该地区陆地水储量呈现直线减少,并达到最低点,这一现象与该时段中国北方的干旱事件相一致。     

GRACE卫星时变重力场的小波多尺度分解——以2008年汶川Ms8.0大地震为例

根据UTCSR发布的GRACE卫星资料,计算了汶川及其周边地区的卫星时变重力场,并利用小波多尺度分析方法对卫星重力场进行分解,得到了反映不同深度的重力场细节和逼近。通过研究各个细节和逼近成分的重力场变化,结合区域构造运动的特点,对GRACE观测的时变重力场及汶川地震的动力机制进行了初步解释和讨论。结果表明,4阶细节和4阶逼近较好地揭示了汶川地区构造形变和深部物质流动引起的重力变化,而且这种震前重力变化特征满足地震孕育发生的条件。     

基于改进短弧积分法的GRACE重力反演理论、方法及应用

正CHAMP、GRACE和GOCE等卫星重力任务的成功实施,为大地测量学、冰川学、海洋学、水文学等学科提供了诸多高时空分辨率的地球重力场模型。由于GRACE对地球重力场的长波段信号十分敏感,且能以较高的精度恢复中波段重力场信号,因此应用GRACE重力数据恢复时变与静态地球重力场,一直以来备受大地测量学者关注。本文在经典短弧积分法的基础上,对重力场反演理论和方法作进一步的探讨和改进,并用GRACE实测数据解算了静态和时变重力场模型,主要研究成果     

不同重力场模型对青藏高原莫霍面反演的影响

针对当前少有学者研究不同重力场模型对莫霍面地形反演的影响,该文利用Parker-Olderburg方法深入讨论纯GRACE重力场模型、纯GOCE重力场模型以及GRACE和GOCE混合重力场模型对青藏高原地区莫霍面反演的影响。实验结果表明,不同重力场模型所反演的莫霍面地形整体空间分布上比较一致,但也存在10 m量级的细微差异;该地区莫霍面地形分布呈顺时针旋转模式,可以较好反映主要构造结构的特点;最深的莫霍面约为70 km,远超过了其他地区正常地壳的厚度。     

中国西部区域GRACE卫星重力变化

针对研究中国西部长期重力变化的问题,该文利用GRACE月重力场模型,基于Slepian变换构建中国西部区域局部重力场,并采用GIA和GLDAS等模型扣除相应误差,通过时间序列分析方法得到该区域2003—2013年卫星重力变化。结果表明,东天山呈现负重力变化,可能由区域冰川消融引起的;西天山的正重力变化,应与冰川小幅增加和构造运动引起的壳幔物质积累有关;青藏高原内陆正重力变化可能由区域地壳抬升、粘弹性地壳的构造应变以及壳幔物质的质量迁移与积累引起;青藏高原边界区域的负重力变化,应与冰川加速消融以及地下水抽取有关。Slepian方法较好地克服了滤波平滑处理带来的重力变化信号压制及细节平滑,可为中国西部的地壳运动、地震活动和气候变化等研究提供数据支撑。     

联合卫星重力和卫星测高数据确定稳态海洋动力地形

新一代卫星重力探测任务GRACE大大提高了地球重力场模型中长波分量的精度,使得联合卫星测高平均海面分离更精细稳态海洋动力地形成为可能。本文利用T/P(1994年~2003年)和JASON-1(2003年~2005年)卫星测高数据确定了全球30′×30′平均海面高;基于重力场模型WHU-GM-05,计算得到对应于海面高的GRACE海洋大地水准面格网值;利用“移去-恢复法”和高斯滤波求得全球稳态海面地形。与EGM96、R io05、ECCO和GGM02模型进行比较,检验结果表明GRACE任务有效的改善了海洋大地水准面的精度,使得稳态海洋动力地形能够呈现更多细部。     

多源卫星重力场模型精度与频谱分析

针对新一代卫星重力探测技术对地球重力场的频谱贡献问题,该文提出了一种基于GPS/水准数据获取多源卫星重力场模型频谱变化特征的方法。采用GPS/水准外符合检验,有效分析评估了多源卫星重力场模型在中国东、西部地区的精度水平。研究结果表明,以CHAMP、GRACE和GOCE卫星为代表的高-低卫星跟踪卫星、低-低卫星跟踪卫星和卫星重力梯度技术,对地球重力场的频谱贡献分别集中在600km以上的长波和中长波、300km以上的中波、200~350km之间的中短波部分。     

GRACE时变地球重力场模型反演全球质量变化

根据GRACE卫星所得到的2008年1～12月份的重力场模型数据反演各种质量变化的分析结果,求出等效水高。对全球和局部地区的等效水高进行分析,得出关于时变重力场模型能够在较大的时间尺度范围和地理范围内很好地反演全球质量变化的结论。     

EIGEN-GL04C重力场模型的精度分析

本文对GFZ发布的GRACE EIGEN-GL04C重力场模型从球谐系数分析、误差阶方差分析两方面进行了精度评价。研究表明,相较于以往的重力场模型,EIGEN-GL04C重力场模型精度对120阶以下(未包含J2项)的地球中长波部分具有明显的改善。GRACE重力卫星测量已经成为获取地球重力场信息的重要手段。     

利用GRACE数据检测日本Mw9.0地震同震和震后重力变化

基于德国地学中心（GFZ）发布的GRACE RL05月重力场模型数据,考虑全球陆地同化系统陆地水储量的影响,采用300 km的扇形滤波,利用叠积法提取了日本Mw9.0地震的同震和震后重力时变信号,并利用最小二乘拟合的方法计算了两个同震重力变化极值点在地震前后85个月的重力年变率。结果表明：土壤水分和雪水引起的重力变化为-0.34~0.09 μgal;利用GRACE数据检测到的同震重力变化为-5.3~4.2 μgal,与基于PSGRN/PSCMP模型计算的结果在空间分布和量级上具有较好的一致性;震后5 a震中附近区域的重力整体上呈现增加趋势,断层上盘所在的日本海域与下盘所在的太平洋区域重力增加的最大值分别为2.6、4.5 μgal,下盘重力增加较大,可能与断层下盘所在地层的黏滞性相对较低有关。     

青藏高原及其周边地区水储量变化的独立成分分析

采用2003年1月至2013年12月共132个月的GRACE卫星数据反演得到的水储量变化信息,利用独立成分分析方法将水储量变化信息分解成多个信号成分,然后与NOAH和WGHM两个水文模型的分析结果进行比较。成分对比结果表明,GRACE反演得到的水储量变化与NOAH和WGHM水文模型在第一个主成分方面符合很好,相关系数分别是0.884和0.877。说明GRACE反演的水储量变化和水文模型在青藏高原及其周边区域具有很强的一致性。从空间模式上看,GRACE反演水储量变化信号的强度比水文模型信号要大,可能与GRACE反演的水储量变化还包括地下水的变化情况等有关。     

GRACE时变重力场滤波方法研究北大核心CSCD

用3种不同的滤波方法获得了2007年相对于2005年的卫星重力变化图像,并与同期地面重力测量的结果进行比较,对比分析GRACE月重力场滤波方法的优缺点。结果表明,去相关平滑滤波算法优于高斯滤波和直接截断法,且去相关平滑滤波DDK5处理得到的卫星重力动态变化图像与地面观测结果符合最好,表明GRACE卫星时变重力场可以用来分析大区域重力动态变化。     

GRACE时变重力场滤波方法研究

用3种不同的滤波方法获得了2007年相对于2005年的卫星重力变化图像,并与同期地面重力测量的结果进行比较,对比分析GRACE月重力场滤波方法的优缺点。结果表明,去相关平滑滤波算法优于高斯滤波和直接截断法,且去相关平滑滤波DDK5处理得到的卫星重力动态变化图像与地面观测结果符合最好,表明GRACE卫星时变重力场可以用来分析大区域重力动态变化。     

利用GRACE RL05模型监测长江流域水储量变化

本文基于GRACE最新重力场模型RL05序列研究了高斯滤波、Wiener滤波、各向异性滤波三种方法在长江流域水储量变化监测中的适用性,计算了应用三种方法得到的水储量变化速率。通过与长江流域水文模型的比较,高斯滤波平滑半径为430km时所得的结果与Wiener滤波基本一致,但各向异性滤波反演的结果与水文模型更为接近,并且优于前面两种方法。研究结果表明GRACE RL05时变重力场球谐系数误差存在各向异性的分布特征,因此各向异性滤波更适用于GRACE区域水储量变化的研究。