基于现代大地测量手段的青藏高原南缘地壳形变分析

基于2003~2013年GRACE卫星重力、GPS和绝对重力数据对青藏高原南缘地壳形变进行探讨分析。结果表明,拉萨测站1999~2015年垂向隆升速率为1.5 mm/a,显示出拉萨测站持续隆升的特点。结合南加州综合网络构建的远程连续GPS参考站,另外选择TPLJ、CHLM、JMSM测站数据进行分析。GPS结果表明,4个测站2003~2013年平均隆升速率为2.3±0.11 mm/a,GRACE垂向位移时间序列显示4个测站2003~2013年除季节性变化外也具有轻微上升趋势,平均速率为0.35 mm/a。联合绝对重力资料在考虑诸多因素(地壳隆升、侵蚀、GIA等)后得到拉萨测站近年重力变化约为-0.97 μGal/a,青藏高原拉萨测站地下地壳底部增厚速率约为4.4 cm/a,揭示了青藏高原南部地区的地壳形变特征--地壳隆升与底部增厚。     

基于插值理论的GRACE卫星精密轨道内插的研究

利用GRACE卫星高采样率精密轨道,进行了拉格朗日插值、牛顿插值和线性插值方法的研究,并对插值结果进行了比较和分析,探讨了插值方法的插值精度及其应用于低轨卫星轨道内插的适用性。     

�����Ǽ�΢�����ϵͳ�۲�����ģ����Ԥ������

??侫????????????GRACE????????????Ч?????????????????????????GRACE????????о???????????????????????????λ????????????????????λ???????????????????????????????????????????????????????????????λ???????????????????????????????????λ??????????????????????     

基于GRACE卫星重力数据确定地球重力场模型WHU-GM-05

基于卫星轨道运动的能量积分方程,可导出利用卫星跟踪卫星数据求解地球重力场的实用公式.本文在Jekeli给出的公式基础上导出了基于能量守恒方程利用两颗低-低卫星跟踪的扰动位差求解重力位系数的严密关系式.基于两颗GRACE卫星的观测数据,采用本文导出的严密能量积分方法求解得到120阶的GRACE地球重力场模型,命名为WHU-GM-05;将WHU-GM-05模型与国际上同类重力场模型EIGEN-GRACE系列和GGM02S分别在阶方差和大地水准面高等方面作了比较,并与美国和中国的部分地区GPS水准观测值进行了精度分析.结果表明基于本文推导的严密双星能量守恒方程得到的WHU-GM-05重力场模型精度与国际上同类重力场模型的精度相当.     

Intercomparison of the total storage deficit index (TSDI) over two Canadian Prairie catchments

Retrieval of the terrestrial moisture storage dataset from the Gravity Recovery and Climate Experiment (GRACE) satellite remote sensing system is possible when the catchment of interest is of large spatial scale. These dataset are of paramount importance for the estimation of the total storage deficit index (TSDI), which enables the characterization of a particular drought event from the perspective of the terrestrial moisture storage over that catchment. Incidentally, the GRACE gravity signal over the 13,000 km² Upper Assiniboine River Basin on the drought-prone Canadian Prairie is so poor therefore making the computation of the total storage deficit index for this basin infeasible. Consequently, the estimation of the terrestrial moisture storage from other reliable sources becomes imperative in order to enable the computation of the TSDI over this basin.This study explores the utilization of the Variable Infiltration Capacity (VIC) model, a physically based, spatially distributed hydrologic model to simulate the total moisture storage over the Upper Assiniboine River Basin which was then employed in the estimation of the TSDI over this basin for subsequent characterization of the recent Prairie-wide drought. Interestingly, the temporal patterns in the computed TSDI from the VIC model reveal a strong resemblance with the same drought characterization undertaken over the larger adjacent Saskatchewan River Basin, which was accomplished utilizing terrestrial moisture storage from the GRACE-based approach. Additionally, these independent techniques employed in the characterization of the last Prairie drought over the two adjacently situated basins resulted in similar drought severity classification from the standpoint of the total moisture storage deficits over these basins. This study has therefore shown that in the computation of the total storage deficit index over small-scale catchments during anomalous climatic conditions that propagate extreme dryness through the terrestrial hydrologic systems, simulations of the total water storage from a structurally sound model such as the VIC model could be resourceful for the computation of the monthly total storage deficit index if no constraint is placed on the availability of accurate meteorological forcing.     

基于GRACE与降水的干旱研究分析

针对干旱问题,采用GRACE重力卫星数据进行反演监测,并对该技术在监测应用中的问题进行了研究.以陕西省与四川省为例,分析了陆地总水储量及降水对干旱的影响,采用基于两者相结合进行计算的干旱指数来检测干旱的发生,并与实际干旱情况及其他干旱指数进行对比,结果表明:该指数能够较好地检测出实际干旱事件且与其他干旱指数有较强的相关...     

Use of GRACE determined secular gravity rates for glacial isostatic adjustment studies in North-America

Monthly geopotential spherical harmonic coefficients from the GRACE satellite mission are used to determine their usefulness and limitations for studying glacial isostatic adjustment (GIA) in North-America. Secular gravity rates are estimated by unweighted least-squares estimation using release 4 coefficients from August 2002 to August 2007 provided by the Center for Space Research (CSR), University of Texas. Smoothing is required to suppress short wavelength noise, in addition to filtering to diminish geographically correlated errors, as shown in previous studies. Optimal cut-off degrees and orders are determined for the destriping filter to maximize the signal to noise ratio. The halfwidth of the Gaussian filter is shown to significantly affect the sensitivity of the GRACE data (with respect to upper mantle viscosity and ice loading history). Therefore, the halfwidth should be selected based on the desired sensitivity.It is shown that increase in water storage in an area south west of Hudson Bay, from the summer of 2003 to the summer of 2006, contributes up to half of the maximum estimated gravity rate. Hydrology models differ in the predictions of the secular change in water storage, therefore even 4-year trend estimates are influenced by the uncertainty in water storage changes. Land ice melting in Greenland and Alaska has a non-negligible contribution, up to one-fourth of the maximum gravity rate.The estimated secular gravity rate shows two distinct peaks that can possibly be due to two domes in the former Pleistocene ice cover: west and south east of Hudson Bay. With a limited number of models, a better fit is obtained with models that use the ICE-3G model compared to the ICE-5G model. However, the uncertainty in interannual variations in hydrology models is too large to constrain the ice loading history with the current data span. For future work in which GRACE will be used to constrain ice loading history and the Earth's radial viscosity profile, it is important to include realistic uncertainty estimates for hydrology models and land ice melting in addition to the effects of lateral heterogeneity.     

Contributions of terrestrial and GRACE data to the study of the secular geoid changes in North America

This paper tests and discusses different statistical methods for modelling secular rates of change of the geoid in North America. In particular, we use the method of principal component/empirical orthogonal functions (PC/EOF) analysis to model the geoid rates from Gravity Recovery and Climate Experiment (GRACE) satellite data. As demonstrated, the PC/EOF analysis is useful for studying the contributions from different signals (mainly residual hydrology signals and leakage effects) to the GRACE-derived geoid rates. The PC/EOF analysis leads to smaller geoid rates compared to the conventional least-squares fitting of a trend and annual and semi-annual cycles to the time series of the spherical harmonic coefficients. This is because we filter out particular spatiotemporal modes of the regional geoid changes.We apply the method of least-squares collocation with parameters to combine terrestrial data (GPS vertical velocities from the Canadian Base Network and terrestrial gravity rates from the Canadian Gravity Standardization Net) with the GRACE-derived vertical motion to obtain again the geoid rates. The combined model has a peak geoid rate of 1.4 mm/year in the southeastern area of Hudson Bay contrary to the GRACE-derived geoid rates that show a large peak of 1.6–1.7 mm/year west of Hudson Bay. We demonstrate that the terrestrial data, which have a longer time span than the GRACE data, are important for constraining the GRACE-derived secular signal in the areas that are well sampled by the data.     

基于GRACE RL06数据监测和分析南极冰盖27个流域质量变化

本文基于CSR最新公布的GRACE RL06版本数据,采用Slepian空域反演法估算了南极冰盖27个流域的质量变化.Slepian空域反演法结合了Slepian空间谱集中法和空域反演法的技术优势,能够有效降低GRACE在小区域反演时信号出现的严重泄漏和衰减,进而精确获得南极冰盖在每个流域的质量变化.相对于GRACE RL05版本数据,RL06在条带误差的控制上要更加优化,获得的南极冰盖质量变化时间序列也更加平滑,但在趋势估算上差别并不明显(小于10Gt/a).本文的估算结果显示:在2002年4月至2016年8月期间,整个南极冰盖质量变化速率为-118.6±16.3Gt/a,其中西南极为-142.4±10.5Gt/a,南极半岛为-29.2±2.1Gt/a,东南极则为52.9±8.6Gt/a.南极冰盖损失最大的区域集中在西南极Amundsen Sea Embayment(流域20-23),该地区质量变化速率为-203.5±4.1Gt/a,其次为南极半岛(流域24-27)以及东南极Victoria-Wilkes Land (流域13-15),质量变化速率分别为-29.2±2.1Gt/a和-19.0±4.7Gt/a,其中Amundsen Sea Embayment和南极半岛南部两个地区的冰排放呈现加速状态.南极冰盖质量显著增加的区域主要有西南极的Ellsworth Land(流域1)和Siple Coast(流域18和19)以及东南极的Coats-Queen Maud-Enderby Land (流域3-8),三个地区质量变化速率分别为17.2±2.4Gt/a、43.9±1.9Gt/a和62.7±3.8Gt/a,质量增加大多来自降雪累积,比如:Coats-Queen Maud-Enderby Land在2009年和2011年发生的大规模降雪事件,但也有来自冰川的增厚,如:Siple Coast地区Kamb冰流的持续加厚.此外,对GRACE估算的南极冰盖质量变化年际信号进行初步分析发现,GRACE年际信号与气候模型估算的冰盖表面质量平衡年际信号存在显著的线性相关关系,但与主要影响南极气候年际变化的气候事件之间却不存在线性相关关系,这说明南极冰盖质量变化的年际信号主要受冰盖表面质量平衡的支配,而气候事件对冰盖表面质量平衡的影响可能是复杂的非线性耦合过程.     

利用广义三角帽方法评估GRACE反演中国大陆地区水储量变化的不确定性

在无真实观测值的情况下,本文利用广义三角帽方法评估了五种GRACE时变重力场模型（CSR、GFZ、GRGS、HUST发布的球谐系数解和JPL发布的Mascon解）反演中国大陆地区2003-2013年水储量变化的不确定性.研究结果表明,CSR、GFZ、JPL、HUST和GRGS反演月水储量变化不确定性的区域平均RMS分别为14.4 mm、26.3 mm、25.3 mm、26.6 mm和56.1 mm,其中GRGS的结果未恢复泄漏信号;在季和年尺度上,模型的不确定性均小于月尺度;扣除周期和趋势信号后,各模型反演结果更为一致.除长江流域外,CSR在13个流域的不确定性均小于其他模型,GRGS反演各流域水储量变化的不确定性通常较大,且可能高估了温带大陆性气候地区水储量的波动;CSR和JPL的不确定性受流域周边水文特征、气候类型、流域面积和形状的影响相对较小,不确定性变化范围分别为2.3~17.1 mm和5.6~22.5 mm,GFZ和HUST受影响较大,不确定性变化范围分别为5.5~35.1 mm和4.0~40.6 mm.本文的研究结果为GRACE产品不确定性评估提供了新的途径,为GRACE时变重力场模型的选取提供参考.