利用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月份左右,该区域水储量较同期呈现明显减小,这与该时段的干旱事件相一致。     

月球重力场模型GL0660B特征分析及其对绕月卫星轨道的影响

月球重力场可用来研究月球演化过程和内部结构,是影响绕月卫星精密定轨的重要因素。基于GRAIL任务数据解算的GL0660B重力场模型,极大提高了月球重力场空间频谱信号的强度和范围。本文首先通过计算相应重力场的阶方差和地形相关性分析,对GL0660B模型进行了精度分析;其次,利用GL0660B模型和其他几个月球重力场模型进行比较,对月球重力场的特征进行了分析;然后通过绘制GL0660B模型和LP150Q模型在月球外部不同高度处的重力异常图,分析比较了月球重力场模型在不同高度上所反映的月球重力场的特征和差异;最后,利用GEODYN软件模拟计算了不同高度卫星的轨道变化。可以看出绕月卫星离心率随时间的变化,以及周期性变化趋势,而且不同高度卫星轨道处质量瘤的摄动影响不同,远月点、近月点和偏心率的变化也存在差异。     

空间直线拟合的混合结构总体最小二乘算法

针对空间直线拟合不能直接采用最小二乘或总体最小二乘方法进行求解的问题,该文提出了一种基于混合结构总体最小二乘法的空间直线拟合算法。首先,由空间直线的参数方程转化成空间直线的通用函数模型;然后,根据得到的函数模型求解参数平差值,求解过程考虑了系数矩阵中常数列对平差结果的影响,顾及了系数矩阵中不同位置的重复元素得到相同改正数的事实,符合实际情况;最后,通过实例验证了算法的有效性。     

利用时变重力数据反演地下水储量变化

为了反映中国陆地区域地下水储量的变化情况,该文利用2003—2019年间GRACE、GRACE-FO重力卫星数据,对8个典型区域地下水储量变化情况进行了研究,并结合气象资料从相关性上分析各区域地下水储量显著变化的原因。结果表明,中国东南大部分地区地下水储量逐年增加,地下水主要靠降水补给;华北平原等人口稠密区地下水亏损严重,研究时段内持续呈下降趋势,降水仅能缓解地下水储量的亏损速度;天山山脉、念青唐古拉山脉等冰川区质量变化和温度异常的相关性较好,这些地区的质量亏损可能是冰川消融引起的。     

郯庐断裂与华北地块扬子地块的拼合

根据水准测量的结果以及地层,岩石,矿物方面的地质证据对华北与扬子地块的拼合界线和郯庐断裂的性质进行了讨论。     

基于多颗低轨卫星的SLR台站坐标几何解算方法

卫星激光测距（SLR）技术作为卫星精密定轨手段和轨道检核重要方法,激光反射器已经成为重力卫星和测高卫星等低轨卫星的基本载荷.经典的SLR台站坐标是使用动力学方法计算的,本文根据多颗低轨卫星（LEO）多历元的激光观测数据,采用几何方法开展地面SLR测站坐标计算.通过组建低轨卫星群实现对全球激光站的动态观测,为了合理配置不同低轨卫星间观测值权重,削弱低轨卫星群可能存在的系统性偏差,提出采用方差分量估计组合的最小二乘法进行解算.实测结果显示,解算出SLR台站坐标框架解与SLRF2014差异平均值在25.1 mm,外符合精度达到1~2 cm.该方法避免了复杂动力学模型,SLR台站坐标的几何计算方法既可以作为激光测站框架解算手段之一,同时将LEO卫星群作为空间并址站实现不同技术地球参考框架间的融合.     

基于高精度重力/地形模型研究亚速尔群岛弹性厚度

亚速尔群岛是大西洋中脊一个扩散中心,研究其内部的动力学作用机制是当前的研究热点.基于高精度的地球重力场模型计算得到的自由空气重力异常与最新的全球地形模型,进行频谱分析,克服传统频谱分析方法中频谱泄露的缺陷,采用多窗口的算法,计算了亚速尔群岛的弹性厚度,研究结果表明该地区海底岩石圈是比较年轻的,地壳厚度为12km,弹性厚度为5km.     

海量卫星重力梯度观测数据确定地球重力位模型的数值方法

基于空域最小二乘法,对卫星重力梯度数据确定地球重力场中的Cholesky分解法、预条件共轭梯度法以及OpenMP并行算法3种数值方法进行比较与分析。研究表明,在计算机硬件资源有限的情况下,传统的Cholesky分解法已经无法满足求解要求;预条件共轭梯度法的求解效率较之Cholesky分解法有改进,但其以损失小量精度为代价;OpenMP并行算法在不损失求解精度的条件下,可提高求解的效率。     

基于卫星测高交叉点的海洋表面地转流速度

在流体静力平衡状态下,海洋Coriolis力和压力梯度平衡就形成地转流,世界上大多数海流都近似为地转流.本文利用卫星测高交叉点方法计算海洋表面地转流速度,分析了利用测高交叉点计算地转流速度的不确定性,上升和下降弧段的海面倾斜在分辨率50 km上可以达到10^-7量级,才可能获得优于10 cm/s的地转流速度.在低纬度或者纬度接近卫星轨道倾角的地区,由交叉点方法计算的地转流速度精度低于中纬度地区.以中国台湾东部黑潮为试验区,利用最新的中国台湾周边海域大地水准面模型参考场计算高精度的大地水准面高,利用TOPEX/Poseidon和Jason-1的GDR数据(2002~2005年)计算海面高,然后计算交叉点的动力高,确定交叉点的地转流速度,结果与中国台湾NCOR(National Center for Ocean Research)的流速基本一致.     

近岸海域EnviSat卫星测高波形重定的Threshold优化算法

波形重定是改善近岸海域卫星测高数据精度的一种有效方法.分析了EnviSat雷达测高波形重定算法,并对Threshold算法进行了优化,开发了波形重定程序;以地中海为实验区,利用几种波形重定方法,对近岸海域的EnviSat测高波形进行重定,认为优化的Threshold算法是几种波形重定算法中最稳健的重定算法,适合于地中海近岸海域EnviSat测高波形的重定,对其他近岸海域的应用也有一定的参考意义.