利用GNSS垂直位移研究长江流域陆地水储量变化的时空分布特征

连续运行的GNSS参考站网络能够实时监测地表水文负荷变迁引起的弹性地壳垂直变形.因此,利用GNSS观测的水文负荷垂直位移估计区域陆地水储量变化是一种行之有效的方案.本文收集了长江流域98个GNSS台站2011—2020年的垂直位移时间序列,并采用Slepian基函数方法将其转化为相应的位移谱,基于质量负荷理论估算长江流域陆地水储量变化,然后联合GRACE、GLDAS和降水数据分析其时空分布特征及水文驱动机制.GNSS反演结果与GRACE、GLDAS水储量变化均表现出明显的季节特征,周年振幅的空间相关系数分别为0.79和0.91.三种结果均表现出年水储量变化东西部大、中部小的空间模式,但GNSS反演结果显示整个长江流域水储量变化的最大周年振幅为～214 mm,明显大于GRACE (～121 mm)和GLDAS (～107 mm)的结果.此外,本文深入调查了长江流域三个子区域(金沙江流域、传统上游及长江中下游)水储量变化的空间格局与时序特征,发现基于三个数据集的水储量变化周年振幅在金沙江流域均呈现出西南向东北递减的趋势,且GNSS与GRACE结果显示长江中下游水储量变化较为显著.GNSS时...     

应用GRACE卫星重力数据计算陆地水变化的相关进展评述

GRACE重力卫星自2002年3月发射至今,已进行了十多年的连续观测,由此获得的重力场变化数据被广泛地应用于研究地表河流及地下水储量变化、南极和格林兰岛冰盖厚度以及全球海平面变化等.本文从GRACE重力卫星数据的处理方法入手,对其数据特点、限制条件和水文模型计算方法等问题进行了系统总结,并针对近年来利用GRACE卫星数据开展的相关研究,从估算全球水储量的变化、区域水储量变化、地下水变化以及陆地河流流域的水储量变化等方面对相关研究和应用进行了简要评述.最后,对使用GRACE卫星数据反演陆地水储量的验证方法和存在的问题进行讨论.本文对于全面了解近年来应用GRACE卫星数据研究陆地水储量变化方面的相关进展具有参考意义.     

GRACE/GRACE-FO重力卫星揭示近20年河西走廊陆地水储量时空变化

河西走廊由疏勒河流域、黑河流域和石羊河流域组成,水资源保护对河西走廊生态平衡和经济发展有着重要意义.本文利用JPL GRACE/GRACE-FO Mascon模型反演该区域陆地水储量的时空变化,结合GLDAS模型、实测地下水位和冰川水模型等数据对陆地水储量进行水平衡分析及时空特征变化分析,结果表明：(1)2002-04—2020-01间由于降水和冰川融水的补充,疏勒河流域南部和黑河大部分区域陆地水储量空间变化呈上升趋势,而蒸散消耗与农业扩张则导致疏勒河流域北部和石羊河流域陆地水储量下降;(2)通过水平衡研究发现人类耗水是疏勒河流域、黑河流域和石羊河流域陆地水储量变化的重要因素,平均贡献率分别为-24.49%、-47.20%和-43.29%;(3)河西走廊水资源治理政策的实施减少了农业灌溉耗水量、控制了耕地面积的扩张、抑制了地下水储量的消耗.     

利用GRACE时变重力场探测2010年中国西南干旱陆地水储量变化

GRACE(Gravity Recovery And Climate Experiment)卫星计划为监测陆地水储量变化提供了有效技术手段.本文采用2003至2010年共计8年的GRACE月重力场模型反演中国西南区域陆地水储量变化,与GLDAS(Global Land Data Assimilation System)全球水文模型进行对比分析,其结果在时空分布上均符合较好,同时在2009年秋至2010年春该区域陆地水储量均呈现明显减少,与该时段云贵川三省的干旱事件相一致;比较分析了2009年秋至2010年春GRACE反演陆地水储量变化与TRMM(Tropical Rainfall Measuring Mission)合成数据计算的月降雨量的时空分布,两组结果均与西南干旱事件对应时段与区域十分吻合;对近8年的陆地水储量变化与月降雨量数据进行相关性分析,其结果表明陆地水储量变化与降雨量强相关,即降雨量是导致陆地水储量变化的主要因素;分析该区域地表温度变化,结果显示2009年9月至2010年3月地表温度均比历史同期高,地表温度的升高加剧了陆地水储量的减少.     

2005—2017年两次强ENSO事件对中国区域陆地水储量变化影响的卫星重力观测

近年来极端气候事件的频发对全球和区域性水循环产生了重大影响,特别是2005—2017年间两次强ENSO(El Nino-Southern Oscillation)事件使得全球陆地水储量出现了较大的年际波动.GRACE(Gravity Recovery and Climate Experiment)重力卫星随着数据质量的提高、后处理方法的完善和超过十年的连续观测,捕捉陆地水储量异常的能力明显提高,这为研究2005—2017年间两次强ENSO事件对中国区域陆地水储量变化的影响提供了观测基础.本文综合利用GRACE卫星重力数据、GLDAS水文模型和实测降水资料分析了中国区域陆地水储量年际变化和与ENSO的关系.研究发现:长江流域中、下游地区和东南诸河流域与ENSO存在较高的相关性,与ENSO的相关系数最大值分别为0.55、0.78、0.70,较ENSO分别滞后约7个月、5个月和5个月.其中长江流域下游地区与ENSO的相关性最强,2010/11 La Nina和2015/16 El Nino两次强ENSO事件使得陆地水储量分别发生了约-24.1亿吨和27.9亿吨的波动.在2010/11 La Nina期间,长江流域下游地区和东南诸河流域陆地水储量异常约在2011年4—5月达到谷值,而长江流域中游地区晚1~2月达到谷值.在2015/16 El Nino期间,长江流域中、下游地区和东南诸河流域陆地水储量从2015年9月到2016年7月持续出现正异常信号.其中,2015年秋冬季(2015年9月至2016年1月)陆地水储量异常明显是受此次El Nino同期影响的结果;2016年春季(4—5月)陆地水异常是受到此次厄尔尼诺峰值的滞后影响所致;2016年7月的陆地水储量异常则与西北太平洋存在的异常反气旋环流有关.     

GNSS(GPS)观测研究地壳运动的新进展

介绍了近几年来GNSS(GPS)系统、观测技术和数据处理方法和技术的发展,数据格式标准化的进展,全球参考框架现状,参考框架对地壳运动研究的影响和IGS的产品。介绍利用GNSS观测研究地壳运动的成果,包括大地震前后及震时的地壳运动结果,时间序列观测结果的研究进展等。     

基于Forward-Modeling方法的黑河流域水储量变化特征研究

黑河流域陆地水储量变化对流域下游等周边区域水资源的合理利用以及经济和社会发展等有着重要的意义.本文利用2003年1月至2013年12月的GRACE RL05数据反演了黑河流域陆地水储量长时间序列的变化,并针对重力场模型和数据处理中产生的信号泄漏问题,采用Forward-Modeling方法进行了改正并恢复泄漏信号;将GRACE获得的泄漏信号恢复前后的黑河流域水储量变化结果与全球水文模型GLDAS和CPC进行比较分析,结果表明泄漏信号改正后的结果与水文模型结果的时间序列相关性均有明显提高,从其空间分布结果可以看出Forward-Modeling方法有效地恢复初始信号、增强被湮没的信号,泄漏信号误差减小;通过分析黑河流域水储量变化的长时间序列结果,发现其具有明显的阶段性变化特征,即2003—2006年呈明显下降趋势,约为-0.86cm·a-1,在2007—2010年趋于平衡状态,而2011—2013年则呈现缓慢上升趋势约为0.14cm·a-1;联合GRACE数据和GLDAS数据反演了黑河流域地下水储量变化,并与全球降雨数据GPCC进行了比较分析,两者相关性可达到0.88以上.     

基于北斗、GLONASS和GPS系统的中低纬电离层特性联合探测

基于GNSS(Global Navigation Satellite Systems)的发展,我们利用具有北斗、GLONASS和GPS三系统信号接收功能的接收机观测的数据,结合电离层总电子含量(Total Electron Content, TEC)的反演算法,提取出GNSS三系统观测的电离层TEC;同时,将GNSS三系统获取的TEC应用到电离层TEC地图、行进式扰动、不规则体结构和电离层的太阳耀斑响应等方面的研究中,这也是首次使用三种GNSS系统数据对电离层进行联合探测研究.研究结果表明,增加了北斗系统的GNSS三系统在研究中国地区电离层TEC地图、周日变化、逐日变化,行进式扰动以及电离层的实时监测等方面较单系统的GPS具有明显的优势.     

基于随机森林模型的GRACE数据3种空间降尺度对比

陆地水储量是赋存在陆地上各种形式水的综合体现,研究其时空变化对认识区域水循环过程和水资源调控等具有重要意义。然而现有陆地水储量变化数据实际分辨率较低,限制了其在中小流域或地区中的应用。针对这一问题,本文基于GRACE重力卫星和其后续卫星GRACE-FO反演的陆地水储量变化数据,首先采用随机森林模型,分别基于格点、区域(流域)和区域(全国)3种空间降尺度思路将GRACE数据降尺度至0.25°×0.25°,后结合GLDAS模型数据,基于水量平衡原理计算得到地下水储量变化数据,最后基于降尺度模型模拟效果和实测地下水位数据评估3种降尺度思路在全国的适用性。结果表明:随机森林模型能够较好地模拟驱动数据(降水、气温、植被条件指数和土壤水储量)与GRACE数据的统计关系,验证期格点降尺度思路的平均相关系数总体在0.6左右,区域降尺度思路的平均纳什效率系数、相关系数和均方根误差分别>0.5、>0.75和<6.6 cm,3种空间降尺度思路的模拟精度均满足基本要求;2003—2021年间,GRACE数据、格点降尺度、区域降尺度(流域)和区域降尺度(全国)得到的我国陆地水储量亏缺量分别约为...     

地球模型对季节性负荷形变计算的影响——以川滇地区为例

利用川滇地区的GNSS和GRACE数据,结合不同地球模型和负荷理论,研究了地球模型对地表季节性负荷形变计算的影响,该工作对于选取合适的地球模型开展负荷形变研究具有一定的参考价值。研究表明：(1)川滇地区GNSS观测的地壳垂向季节性形变振幅为20 mm左右,GRACE反演的垂向形变与GNSS的结果相位一致,振幅存在差异。(2)区域地球模型的负荷勒夫数与其他地球模型的差异较大,且负荷勒夫数h_n对地球结构的变化较为敏感。(3)区域地球模型可以改善GRACE反演的负荷形变结果,从而减小与GNSS观测结果的差异。(4)川滇地区大部分GNSS测站的加权均方根比值减小量呈现由东北向西南方向逐渐增加的变化趋势。     

利用GRACE卫星重力数据监测关中地区地下水储量变化

关中地区作为一带一路重要的工农业发达地区之一,开展针对该地区地下水储量变化的监测和分析工作对揭示地下水储量变化特征与经济社会发展具有重要现实意义.本文基于2003—2014年GRACE卫星重力场模型数据,采用组合滤波及单一尺度因子方法反演了关中地区陆地水储量变化,扣除GLDAS地表水平均结果,对关中地区地下水储量变化进行了监测分析.将陆地水储量变化与GLDAS进行相关性分析,将地下水储量变化与WGHM地下水模型及实测地下水位结果进行对比分析.研究结果表明：①关中地区陆地水变化与GLDAS模型结果具有较强的相关性,相关系数多数大于0.7,其中与模型平均结果的相关系数可达0.8.② 2003—2008年关中地区地下水呈正增长趋势,增加速率为0.25 cm·a^-1,与同期实测数据变化趋势一致;但2003—2013年地下水存在长期亏损,亏损速率为-0.37 cm·a^-1等效水高,这与同时期WGHM估算结果-0.35 cm·a^-1十分吻合.③关中地区地下水存在明显的年变化特征,在2003—2014年期间地下水减少速率为-0.44 cm·a^-1,与该地区降雨量有较好的一致性,在降雨偏少的2008、2012和2013年,地下水也显著减少.     

利用GPS和GRACE分析四川地表垂向位移变化

陆地水储量的季节性变化是导致地表周期性负荷形变位移的主要因素,有效地剔除地表位移中的陆地水储量影响,是获取地壳构造垂向运动的必要过程.四川地处青藏高原东边缘,地形分区明显,境内以长江水系为主,水资源丰富,研究四川地区地表负荷形变位移,有助于分析陆地水储量的时空分布特性及地壳构造形变信息.本文利用研究区域内59个CORS站的GPS观测数据,计算了CORS站点的垂向位移,并将其与GRACE所得相应结果进行对比分析.结果显示,GPS和GRACE所得垂向位移时间序列的振幅大小整体相符,但存在明显的相位差.GPS站点振幅最大值为12.7 mm,对应HANY站,最小值为1.5 mm,对应SCMX站.GRACE所得的地表垂向位移振幅大小均为3~4 mm,且最大位移集中出现在7-9月份;而GPS站点出现最大位移的月份和地形相关,东部盆地、西北部高原和南部山地分别出现在7-8月份、10-11月份和10月份.GPS站点时间序列中的周年项与陆地水的季节性变化强相关,为了讨论陆地水储量对GPS站点位移的影响,本文利用改进的总体经验模态分解方法（MEEMD：Modified Ensemble Empirical Mode Decomposition）,从GPS垂向位移时间序列中提取出周年项及约2年的年际变化项.发现利用MEEMD获取的周年项改正原始GPS时间序列时可使其WRMS（Weight Root Mean Square）减少量减小约26%,结果优于最小二乘拟合方法提取的GPS周年项改正效果,验证了MEEMD方法在GPS坐标时间序列处理中的可行性及有效性.     

GRACE估算陆地水储量季节和年际变化

利用最新公布的GRACE GFZ RL04数据,分析了2003年1月～2007年12月全球27条流域和陆地水储量的季节性和年际变化.结果表明,相近流域季节性变化相位接近.2003年1月～2007年12月陆地水储量季节性变化为1572.4 km3,其中变化最大流域为亚马逊河,其次分别为鄂毕河、尼罗河和尼日尔河等流域.5年来 GRACE陆地水储量的年际变化为-75.4±40.3 km3/a,其中亚马逊河、勒拿河和马更些河等流域的年际变化呈现正增长,而刚果河、密西西比河、恒河、育空河和雅鲁藏布江等流域则相反.GRACE与GLDAS数据均表明2006年后陆地水储量年际变化存在明显增加.     

联合高频GNSS和强震数据快速反演四川泸定6.8级地震断层破裂特征

北京时间2022年9月5日12时52分,四川甘孜藏族自治州泸定县发生6.8级地震.利用震中附近1 Hz高频GNSS观测数据获取了同震速度和位移波形,并快速测定了泸定地震的震中和震级.实验结果表明：高频GNSS反演的震中与美国地质调查局(USGS)发布的震中相差32 km,与中国地震台网中心发布值相差16 km;高频GNSS反演的震级,与两个机构均仅差0.1个震级单位.针对地震预警、震后快速响应等时效性应用,提出了一种联合高频GNSS和强震数据的线源破裂特征快速反演方法.泸定地震实验结果表明：在震后20 s时可获得稳定的线源模型,破裂长度、方向和破裂模式值分别为33.3 km、151°和0.6,破裂方向与USGS震源机制解断层走向相差14°,反演的断层破裂模式为双侧破裂.提出的地震断层破裂特征快速反演方法可用于地震预警、震后灾害快速评估以及紧急响应,同时可为今后联合高频GNSS和强震数据快速测定地震破裂特征提供参考.     

基于GRACE和GLDAS水文模型反演安徽省地下水储量变化

地下水储量的有效监测是实现区域水资源管理的重要依据,传统监测方法存在各自的局限性导致其实现较为困难.本文提供一种使用GRACE卫星重力数据与GLDAS水文模型数据反演得到安徽省区域地下水储量变化的方法.利用2002年4月-2017年6月不同机构GRACE卫星重力数据的综合解,经过DDK去相关光滑滤波与退卷积法分别消除或削弱南北条带误差、改正信号泄露反演得到安徽省陆地水储量变化,扣除由GLDAS水文模型数据获取的地表水储量变化,得到安徽省地下水储量变化时间序列,并结合国家统计局官方发布的安徽省地下水资源量进行初步验证分析.研究结果表明:安徽省地下水储量变化长期表现为波动下降趋势,其年变化率约为-5.35 mm·a-1,且呈现出自东南方向至西北方向逐次递减的显著空间差异;地下水储量表现出明显的季节性变化,夏季和冬季地下水储量呈现回升趋势,春季和秋季呈现出下降趋势.除去反演过程存在较大干扰因素的情况,反演结果与国家统计局官方数据的相关系数达到89.62％,因此本文反演得到的地下水储量变化的结果是相对可靠的.     

水质量分布变化对重力观测的影响

利用全球7000多个陆地气象台站观测资料、TOPEX/Poseidon卫星测高资料、多层海水温盐度数据和负荷格林函数,采用数值积分法,计算了陆地水储量及海水质量分布变化引起我国20个测站的重力变化.结果表明,我国一些测站观测的重力受水质量迁移影响超过3μGal,可见在应用重力仪观测资料研究地球动力学时,有必要作水负荷改正.另外,文中还对武汉测站的超导重力仪观测的残差和水负荷的计算结果作了比较,并讨论了不同的地球模型及气象数据对计算测站重力的影响.     

全球水质量迁移对海平面空间模式周年变化的影响

大气、陆地水和海洋之间的水质量迁移对海平面的影响一般假定为均匀薄层分布.但实际上,水质量负荷重新分布一方面会使地壳产生形变,另一方面会引起重力位势场变化(引力位和离心力位),这都会对海平面时空变化特征产生影响,两方面之和称为负荷自吸引效应(SAL).海洋模式模拟的时变洋底压力结果一般符合Boussinesq假设即体积守恒,忽略了大气、陆地水和海洋之间水质量交换的影响.本文基于海平面变化方程,联合陆地水模型、大气地表气压模型、海洋洋底压力模型和GRACE反演的冰川质量变化,详细讨论了2003-2010年SAL对海平面周年变化的影响.主要结论有:(1)SAL对全球海平面周年变化有显著影响,振幅在1.3~19 mm.其中近海岸和低纬度区域受影响较大.(2)在SAL引起的海平面周年振幅变化的因素中,陆地水储量变化因素最大,大气因素次之,非潮汐海洋影响最小.但非潮汐海洋对海平面周年相位空间变化的影响最为复杂.(3)通过与国际长期验潮站观测数据结果比较,在ECCO海洋模式估计的洋底压力结果中引入SAL,能多解释约5.3%观测信号方差.     

利用断层应力积累探讨大地震发生的地点——以日本2011年东北9级大地震为例

大地震在哪里发生是地震预报首先要解决的问题.利用反演GNSS观测数据得到的2011年日本东北9级大地震前7年(2004-2010年)断层上的应力变化,我们发现了这次地震断层的孕震区.为了进一步研究该孕震区的演化过程,本文继续反演这次大地震在1997-2003年间的断层应力变化过程.通过这两期的反演工作,我们看到,在这1...     

利用精密单点定位求解电离层延迟

近年来,高时空分辨率的全球导航卫星系统(GNSS)观测信号已成为电离层研究的重要资源.利用GNSS研究电离层,需首先将观测资料转换成包含电离层信息的可观测量(Ionospheric 0bservables,称之为"电离层观测值").目前,最常用的电离层观测值一般采用联合无几何影响组合的码和相位观测,利用相位平滑伪距方法...     

联合GNSS和InSAR观测位移反演2008年汶川大地震断层位错模型参数

利用现代空间大地测量技术,尤其是卫星合成孔径雷达干涉测量,能够获取高精度、高空间分辨率的同震和孕震形变,为地震断层形变和破裂机制研究提供了前所未有的机遇。本文介绍了利用大地测量观测数据反演地震断层位错模型参数的贝叶斯反演方法。联合运用2008汶川大地震前后GNSS和InSAR技术观测获得的同震位移,反演了地震断层的几何参数和滑动位错分布。研究结果表明,汶川地震的断层滑动主要集中在倾角较陡的浅部,同时包含逆冲和右旋走滑,其中最大逆冲6.1m,最大右旋6.5m。根据断层滑动分布正演计算得到的上盘同震位移明显小于下盘,预示该断层两侧孕震形变可能存在较大的不对称性。