我国东部海域高精度地转流的确定方法

针对目前带有高频信息的全球稳态海面地形模型的计算过程比较复杂的问题,该文采用移去恢复法,对由全球海面高模型得到稳态海面地形模型添加高频信息,并用新模型计算我国东部海域的地转流,分析其精度。结果表明:使用卫星重力场模型GOCE_TIM5,全球平均海面高模型CNES-CLS15,计算出分辨率为15′×15′我国东部海域的地转流,高频信息不明显;采用移去恢复法,利用GOCE_TIM5模型、CNES-CLS15模型、全球稳态海面地形模型CNES-CLS13,得到相同分辨率我国东部海域的地转流,与未采用移去恢复法相比,其结果中细节信息明显,与全球实测浮标数据更为接近,其精度与由全球稳态海面地形模型CNES-CLS13得到结果相当。     

利用参考重力场模型基于能量法确定GRACE加速度计校准参数

利用多个参考重力场模型分别对GRACE一个月的实测加速度计观测数据进行检校.数值计算结果的比较分析表明了利用参考重力场模型确定加速度计校准参数是有效的.     

ICESAT卫星确定南极冰盖高程模型研究

研究了星载激光测高的理论和技术,利用ICESAT测高数据确定了南极冰盖高程模型.将ICESAT测高数据与南极实测GPS数据进行比较,结果显示ICESAT数据能够很好地确定南极冰盖高程模型.     

利用GRACE反演昆士兰州地下水储量变化

昆士兰州作为澳大利亚的第二大州,开展针对该地区地下水储量变化的监测与分析对当地生态环境与用水管理政策制定具有重要现实意义. 利用德克萨斯大学空间研究中心(CSR)发布的重力反演与气候实验卫星(GRACE)时变重力场模型与全球陆地数据同化系统(GLDAS)地表同化模型,对昆士兰州地下水储量时空变化进行监测分析,并与实测水井数据和全球降水气候学项目(GPCP)降雨资料进行验证分析. 研究结果表明:时间上昆士兰州地下水在2003—2015年约以1.3±0.09 cm/a的速率递增,空间上呈现东增西减的显著空间差异;与GPCP降雨资料对比发现,降雨是影响地下水储量变化的主要因素;与水井实测数据对比发现地下反演结果与水井水位变化趋势基本一致.      

基于星载激光雷达和雷达高度计数据的南极冰盖表面高程制图

南极数字高程模型DEMs(Digital Elevation Models)是研究极区大气环流模式,南极冰盖动态变化和南极科学考察非常重要的基础数据。目前,科学家已经发布了五种不同的南极数字表面高程模型。这些数据都是由卫星雷达高度计,激光雷达和部分地面实测数据等制作而成。尽管如此,由于海洋与冰盖交接的南极冰盖边缘区随时间的快速变化,有必要根据新的卫星数据及时更新南极冰盖表面高程数据。因此,我们利用雷达高度计数据(Envisat RA-2)和激光雷达数据(ICESat/GLAS)制作了最新的南极冰盖高程数据。为提高ICESat/GLAS数据的精度,本文采用了五种不同的质量控制指标对GLAS数据进行处理,滤除了8.36%的不合格数据。这五种质量控制指标分别针对卫星定位误差、大气前向散射、饱和度及云的影响。同时,对Envisat RA-2数据进行干湿对流层纠正、电离层纠正、固体潮汐纠正和极潮纠正。针对两种不同的测高数据,提出了一种基于Envisat RA-2和GLAS数据光斑脚印几何相交的高程相对纠正方法,即通过分析GLAS脚印点与Envisat RA-2数据中心点重叠的点对,建立这些相交点对的高度差(GLAS-RA-2)与表征地形起伏的粗糙度之间的相关关系,对具有稳定相关关系的点对进行Envisat RA-2数据的相对纠正。通过分析南极冰盖不同区域的测高点密度,确定最终DEM的分辨率为1000 m。考虑到南极普里兹湾和内陆地区的差异性,将南极冰盖分为16个区,利用半方差分析确定最佳插值模型和参数,采用克吕金插值方法生成了1000 m分辨率的南极冰盖高程数据。利用两种机载激光雷达数据和我国多次南极科考实测的GPS数据对新的南极DEM进行了验证。结果显示,新的DEM与实测数据的差值范围为3.21—27.84 m,其误差分布与坡度密切关系。与国际上发布的南极DEM数据相比,新的DEM在坡度较大地区和快速变化的冰盖边缘地区精度有较大改进。     

基于GRACE卫星数据的高精度全球静态重力场模型

应用GRACE卫星数据反演高精度静态地球重力场是大地测量学界的热点之一。考虑到经典动力学法线性化误差随弧长拉长而迅速增长,本文以GRACE卫星轨道观测值为初值的线性化方法,建立了应用GRACE卫星轨道和星间距离变率反演地球重力场的改进动力学法理论模型。利用2003年1月至2010年12月的GRACE卫星姿态、轨道、星间距离变率和非保守力加速度等观测数据,解算了一个180阶次的无约束全球静态重力场模型Tongji-Dyn01s和一个采用Kaula规则约束的全球重力场模型Tongji-Dyn01k。与国际不同机构最新发布的纯GRACE数据解算的重力场模型(包括AIUB-GRACE03S、GGM05S、ITSG-Grace2014k和Tongji-GRACE01)进行比较,并利用DTU13海洋重力异常和GPS/水准高程异常进行外部检核,结果表明,Tongji-Dyn01s与国际最新模型精度处于同一水平,然而Tongji-Dyn01k模型总体上更加靠近EIGEN6C2重力场模型。     

论高精度卫星重力场模型和厘米级区域大地水准面的确定及水文学时变重力效应

阐述了联全新一代卫星重力测量数据、卫星测高数据及全球陆地重力数据确定高精度180阶全球重力场模型、以全球重力场模型为框架参考场、利用我国地面重力数据、GPS水准资料、数值高程模型和地形密度信息确定高分辨率cm级区域大地水准面的思想。指出了一个重点发展方向:利用GRACE卫星每30天的重力位模型分析时变重力场,联系合卫星测高同时相平均海面以及水文、气候和海洋模型,分析我国黄河流域和海洋地区水储量分布和海流季节性变化,并解释与气候要素变化的相关性。     

基于GRACE/GRACE-FO反演华北平原2002—2021年水储量变化

针对华北平原水资源分布不均的问题，该文利用2002年4月—2021年12月共237个月的重力恢复和气候实验(GRACE)卫星、GRACE-FO卫星重力数据对华北平原的水储量变化情况进行了研究，采用奇异谱插值填补GRACE与GRACE-FO间的空白期数据，并联合3家Mascon产品、全球陆面同化系统水文模型与全球降水气候中心降水数据进行综合性比较和分析，实现了对华北地区水储量时空变化情况的探究。结果表明，华北地区陆地水储量整体呈下降趋势，亏损速率为-0.927 6 cm/a,且具有明显的季节特征，与降水量的变化情况一致，经过相关性分析，GRACE、GRACE-FO与3个Mascon产品反演的水储量变化存在较高的相关性，相关系数均在0.89以上。     

基于GRACE的中国区域水储量变化研究

水是人类生活中广泛使用的重要资源,一旦开发利用不当,不仅会制约国民经济发展、破坏生态环境,还会引起水位持续下降、水质恶化和地面沉降等现象.而运用GRACE卫星时变重力场模型可以有效计算出陆地总水储量变化,同时结合相应水文模型还可以反演出地下水储量的变化,其相对于传统监测地下水方法具有简便、可以大范围进行监测以及预防由水储量变化引起的灾害等优势.针对水储量变化问题,采用GRACE重力卫星数据进行反演监测,将处理后的GRACE数据与GLDAS水文模型相结合,以反演中国地区的地下水储量变化趋势,结果表明:地下水储量变化与陆地总水储量变化趋势较为一致,且中国一半地区的地下水储量处于下降趋势,其中下降严重的有河北南部、新疆北部和云南北部等地区,最大变化速率分别达到了-1.72 cm/a、-2.44 cm/a和-3.01 cm/a.最后以陕西省与四川省的实测地下水数据来检验反演地下水变化的准确性,并分析降水与地下水变化的关系.     

联合GRACE/GOCE重力场模型和GPS/水准数据确定我国85高程基准重力位

GRACE、GOCE卫星重力计划的实施,对确定高精度重力场模型具有重要贡献。联合GRACE、GOCE卫星数据建立的重力场模型和我国均匀分布的649个GPS/水准数据可以确定我国高程基准重力位,但我国高程基准对应的参考面为似大地水准面,是非等位面,将似大地水准面转化为大地水准面后确定的大地水准面重力位为62 636 854.395 3m~2s~(-2),为提高高阶项对确定大地水准面的贡献,利用高分辨率重力场模型EGM2008扩展GRACE/GOCE模型至2190阶,同时将重力场模型和GPS/水准数据统一到同一参考框架和潮汐系统,最后利用扩展后的模型确定的我国大地水准面重力位为62 636 852.751 8m~2s~(-2)。其中组合模型TIM_R4+EGM2008确定的我国85高程基准重力位值62 636 852.704 5m~2s~(-2)精度最高。重力场模型截断误差对确定我国大地水准面的影响约16cm,潮汐系统影响约4~6cm。     

联合GOCE卫星数据和GRACE法方程确定SWJTU-GOGR01S全球重力场模型

联合地球重力场和海洋环流探测器（Gravity Field and Steady-State Ocean Circulation Explorer,GOCE）和重力恢复与气候实验（Gravity Recovery and Climate Experiment,GRACE）卫星观测数据确定全球静态重力场模型是当前大地测量学的研究热点之一。联合近3 a的GOCE卫星梯度数据和7 a左右的GRACE星间距离变率数据计算的ITG-GRACE2010S模型的法方程恢复了210阶次的重力场模型SWJTU-GOGR01S。采用带通数字滤波方法处理GOCE卫星的4个高精度梯度观测分量,利用梯度数据恢复重力场模型的观测方程直接建立在梯度仪坐标系中,可以避免坐标转换过程中高精度的梯度观测分量受低精度分量的影响;联合法方程解的最优权采用方差分量估计迭代计算,GOCE数据的两极空白引起的病态问题采用Kaula正则化方法进行约束。基于EIGEN-6C2模型和北美地区的GPS水准网观测数据,对SWJTU-GOGR01S模型进行内外符合精度分析,结果表明,SWJTU-GOGR01S模型在210阶次的大地水准面误差和累计误差分别为1.3 cm和5.7 cm,精度与欧洲空间局公布的第四代时域法模型相当,略优于GOCO02S和GOCO03S模型的精度。     

利用GRACE、GPS和绝对重力数据监测斯堪的纳维亚陆地水储量变化

重力卫星GRACE(gravity recovery and climate experiment)监测斯堪的纳维亚半岛陆地水储量变化会受到冰川均衡调整(GIA)信号的严重影响。首先根据该地区绝对重力和GPS并址观测数据计算了GIA重力和垂直位移的实测线性比值,利用该比值和GPS网观测的垂直位移速度场得到了GIA重力。然后,对GRACE观测的重力变化速率进行GIA重力改正,进而可分离陆地水储量变化趋势,避免了使用GIA模型所带来的巨大不确定性,并根据观测数据完整估计了所得结果的不确定性。最后与水文模型作对比分析。结果表明,实测的GIA重力-垂直位移线性比值为0.148±0.020μGal/mm(1Gal=10-2 m/s2),该结果检验了Wahr的理论近似值且与北美实测的结果非常接近。2003年1月至2011年3月期间,斯堪的纳维亚半岛陆地水储量存在明显的增加趋势,信号的主体位于半岛南端的维纳恩湖附近,总的水量增加速率为4.6±2.1km3/a,数据观测期间的累积增加水量为38±17km3。研究结果与WGHM水文模型的结果有较好的一致性,相关系数达到0.69,而与GLDAS水文模型的相关性略小。     

激光测距对静态重力场恢复精度的影响分析

鉴于下一代重力卫星设计中利用更高精度的激光测距技术代替微波测距,对激光测距提高地球重力场探测精度的问题进行了讨论。通过高精度的动力学重力场模型反演方法,推导了线性化的星间变率公式,并以一定的权融合卫星精密轨道与星间变率数据。通过模拟计算结果可知,当卫星轨道、定轨精度、加速度计精度与Gravity Recovery and Climate Experiment(GRACE)相同,星间变率精度依次从1.0×10~(–6) m/s提高到5.0×10~(–7) m/s、1.0×10~(–7) m/s、5.0×10~(–8) m/s、1.0×10~(–8) m/s时,累积大地水准面误差则在120阶时依次从85.14 cm降低为33.09 cm、7.33 cm、3.70 cm、3.59 cm。结果表明当采用高精度的激光测距后,采用低低卫卫跟踪模式,地球静态重力场模型的探测精度有望比GRACE提高1个量级。当激光测距精度提高至10 nm/s时,计算结果与精度为50 nm/s的计算结果无明显差别,这表明过高的星间变率测量精度相对于其他指标而言有冗余。本文建议在其他测量技术精度未有提高的前提下,只需将激光测距的精度提高至50 nm/s即可。     

利用GRACE资料构造大尺度时变重力场统一模型

基于2003~2012年的GRACE卫星重力资料,采用最小二乘拟合的方法,构建了时变重力场统一模型IGG-TVG2013。该模型以球谐系数的形式表达,在考虑趋势项和周期项等经验参数的基础上,还考虑了加速度项和潮汐模型误差、大地震等因素的影响。将IGG-TVG2013模型与GRACE资料进行了比较分析,在全球92%以上的区域二者符合精度优于±1ugal;利用该模型外推预测了2013年1~6月的重力场变化,结果与GRACE实测数据符合较好。这表明IGG-TVG2013模型不但能较好地描述重力场的连续时空变化,而且具有一定的短期预测能力。     

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

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

海洋重力异常数据的EGM2008模型插补方法

针对加权平均模型中反映信号区域性变化特征的权系数不易确定的问题,提出了利用EGM2008模型计算的重力异常数据表示重力场的趋势性变化以改进传统的加权平均法。基本步骤是利用实测海洋重力异常数据,借助EGM2008全球重力位模型,在移去-恢复技术框架下内插出残差重力异常数据,进而得到内插点的海洋空间重力异常。海区实测数据分析结果表明,本文所述方法计算的内插值,无论是分布情况还是与实测值差值的大小,都明显优于传统的加权平均法。     

重力卫星的星载GPS精密定轨

利用CHAMP和GRACE卫星的实测数据,研究了重力卫星的精密定轨问题,并针对几何法精密定轨方法给出了一种有效的星载数据编辑策略;在PANDA软件的基础上,处理了101 d的实测数据;通过与不同机构卫星轨道的比较、激光测距观测值检验以及重力场模型恢复等外部检核的方式,分析了卫星轨道的精度.结果显示,本文的简化动力学轨道的精度为2～3 cm;几何学轨道的定轨精度为3～4 cm,适用于重力场模型的解算.     

联合Jason-1与GRACE卫星数据研究全球海平面变化

海平面变化主要由海水质量变化和比容海平面变化组成。联合GRACE时变重力场与Jason-1卫星测高数据分别研究了全球平均海平面及其质量分量与比容分量的季节性变化(2002-09—2008-04),联合两种卫星数据所得平均海平面变化在周年振幅和相位上与WOA05海洋模型结果具有一致性。估算出研究时间段内全球平均海平面变化及其两个主要分量的长期性趋势,在超过6.5a的时间尺度上,全球平均海平面以+2.0±0.4mm/a的速率上升,其中海水质量变化分量的贡献为+1.4±0.4mm/a,比容变化分量的贡献为+0.5±0.3mm/a。     

GPS卫星轨道及钟差对GRACE卫星定轨精度影响

针对GPS卫星精密轨道和钟差插值对GRACE卫星定轨精度影响进行了分析,分别使用IGS(International GNSS Service)30 s间隔钟差、CODE(the Center for Orbit Determination in Europe)30和5 s间隔钟差以及15 min精密星历进行GRACE卫星定轨实验。结果表明:GPS轨道插值精度可以达到cm级,将15 min GPS轨道插值为30 s间隔利用9阶拉格朗日插值定轨结果精度最高,继续增加阶数定轨精度不会增加;利用CODE钟差计算GRACE非差运动学轨道,码伪距结果精度较IGS产品提高6%,载波相位运动学定轨结果和约化动力学定轨结果精度都提高10%左右;5 s间隔卫星钟数据对定轨结果改进并不明显。采用CODE间隔为30 s钟差进行GRACE运动学定轨的计算精度能满足cm级轨道的应用需求。     

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

针对京津冀地区水资源灾害频发、水储量监测手段存在局限性的问题,本文利用该地区2015年1月至2017年6月的重力场恢复与气候实验(GRACE)月重力场模型数据反演陆地水储量变化信息,结合区域降水、全球陆地资料同化系统(GLDAS)模型数据,对研究区域内水储量变化时空特征进行相关性分析.结果表明:GRACE监测的陆地水储...