Natural Geological Weighing Lysimeters: Calibration Tools for Satellite and Ground Surface Gravity Monitoring of Subsurface Water-Mass Change

Increased accuracy in measuring temporal variations in the Earth's gravity field allow inprinciple the use of gravity observations to deduce subsurface water-mass changes. This canbe with respect to a small area, or as a larger spatial average of water mass change usinggravity observations from low-altitude satellites, such as the forthcoming GRACE mission.At both scales, there is a need to validate gravity-based estimates against field recordings ofactual subsurface water-mass variations. In practice, this could prove difficult because thespatial integral of all water-storage change components can be subject to considerable fieldmeasurement error. An alternative approach to the validation process is proposed by whichsuitable geological formations are utilized as giant weighing devices to directly measure area-integratedwater-mass changes. The existence of such natural geological weighing lysimetersis demonstrated using observations from a replicated experimental site in New Zealand. Sitesof this type could be used to verify water-storage change estimates derived from sensitiveground surface gravity instrumentation. In addition, geological lysimeters could be used tomake local checks on the accuracy of any estimated regional water-mass time series, whichis proposed for satellite calibration. The land area weighed by a geological lysimeter increaseswith formation depth and it is speculated that recordings made at oil well depth may allowdirect monitoring of subsurface water mass changes at the regional scale.     

重力卫星五年运行对求定地球重力场模型的进展和展望

对近年升空的CHAMP和GRACE和将于2007年升空的GOCE卫星在测定地球重力场方面的技术特点和初步成果进行了回顾、比较和评估。并对它们今后在静态和动态的地球重力场构模方面可能的进展作一展望。现在只用一颗重力卫星的轨道摄动数据,就可以以前所未有的可靠性和精确性来求定地球重力场的长波和部分中波。如CHAMP重力卫星的33个月数据所求定的地球重力场模型,相对于曾利用多颗卫星资料所推算的GRlM5-S1重力场模型,在长波方面的精度和可靠性都有很大改善。而GRACE重力卫星的110天数据所导出成果的空间分辨率,又优于CHAMP的33个月的数据成果。GRACE卫星还有一个重要任务,就是测定重力场非潮汐的短期性或准实时的变化。还介绍了新发表的一个联合地球重力场模型EIGEN-CG03C,360完全阶次,分辨率约30′。CG03C同CHAMP／GRACE以前的重力场模型比较,在400km波长的精度方面改善了一个量级,大地水准面的精度改善了3cm,重力异常的精度改善了0．4mgal。     

On the use of iterative re-weighting least-squares and outlier detection for empirically modelling rates of vertical displacement

The proper identification and removal of outliers in the combination of rates of vertical displacements derived from GPS, tide gauges/satellite altimetry, and GRACE observations is presented. Outlier detection is a necessary pre-screening procedure in order to ensure reliable estimates of stochastic properties of the observations in the combined least-squares adjustment (via rescaling of covariance matrices) and to ensure that the final vertical motion model is not corrupted and/or distorted by erroneous data. Results from this study indicate that typical data snooping methods are inadequate in dealing with these heterogeneous data sets and their stochastic properties. Using simulated vertical displacement rates, it is demonstrated that a large variety of outliers (random scattered and adjacent, as well as jointly influential) can be dealt with if an iterative re-weighting least-squares adjustment is combined with a robust median estimator. Moreover, robust estimators are efficient in areas weakly constrained by the data, where even high quality observations may appear to be erroneous if their estimates are largely influenced by outliers. Four combined models for the vertical motion in the region of the Great Lakes are presented. The computed vertical displacements vary between  − 2 mm/year (subsidence) along the southern shores and 3 mm/year (uplift) along the northern shores. The derived models provide reliable empirical constraints and error bounds for postglacial rebound models in the region.     

冰盖消融的海平面指纹变化及其对GRACE监测结果的影响

全球变暖背景下的冰盖消融以及由此带来海平面上升日益明显,直接影响地球表面的陆地水质量平衡,以及固体地球瞬间弹性响应,研究冰盖质量变化的海平面指纹能够帮助深入了解未来海平面区域变化的驱动因素.本文基于海平面变化方程并考虑负荷自吸效应（SAL）与地球极移反馈的影响,借助美国德克萨斯大学空间研究中心（Center for Space Research,CSR）发布的2003年到2012年十年期间的GRACE重力场月模型数据（RL05）,结合加权高斯平滑的区域核函数,反演得到格陵兰与南极地区冰盖质量变化的时空分布,并利用海平面变化方程计算得到了相对海平面的空间变化,结果表明：格陵兰与南极冰盖质量整体呈明显的消融趋势,变化速率分别为-273.31 Gt/a及-155.56 Gt/a,由此导致整个北极圈相对海平面降低,最高可达约-0.6 cm·a^-1;而南极地区冰盖质量变化趋势分布不一,导致西南极近海相对海平面下降,而东南极地区近海相对海平面上升,最高可达约0.2 cm·a^-1.远离质量负荷区域的全球海平面以上升趋势为主,平均全球相对海平面上升0.71 mm·a^-1,部分远海地区相对海平面上升更加突出（例如北美与澳大利亚）,高出全球平均海平面上升速率将近30%.此外,本文也重点探讨了GRACE监测冰盖消融结果中由于极地近海海平面变化导致的泄漏影响,经此项影响校正后的结果表明：海平面指纹效应对GRACE监测格陵兰与南极地区2003-2012期间整体冰盖消融速率的贡献分别为约3%与9%,建议在后期利用GRACE更精确地估算研究区冰盖质量变化时,应考虑海平面指纹效应的渗透影响.

     

GGM02 – An improved Earth gravity field model from GRACE

A new generation of Earth gravity field models called GGM02 are derived using approximately 14 months of data spanning from April 2002 to December 2003 from the Gravity Recovery And Climate Experiment (GRACE). Relative to the preceding generation, GGM01, there have been improvements to the data products, the gravity estimation methods and the background models. Based on the calibrated covariances, GGM02 (both the GRACE-only model GGM02S and the combination model GGM02C) represents an improvement greater than a factor of two over the previous GGM01 models. Error estimates indicate a cumulative error less than 1 cm geoid height to spherical harmonic degree 70, which can be said to have met the GRACE minimum mission goals. Electronic Supplementary Material Supplementary material is available in the online version of this article at     

2004年苏门答腊地震前后海面高趋势的变化

利用苏门答腊附近海域T/P、Jason-1测高卫星近20年的海面高连续观测资料,分别计算了该区域(80°E~105°E,5°S~20°N)在2004年苏门答腊大地震前后的海面高变化趋势,并与该区域对应时间段GRACE重力卫星反演的地表质量迁移结果进行比较。研究结果表明,由卫星测高观测资料估算的震区地震前后海平面趋势的变化与卫星重力反演结果基本一致。由于卫星测高沿轨观测具有高精度、高时空分辨率的特性,卫星测高资料估算的海平面趋势变化可以更为准确地反映震区海平面变化的局部特征。还对该区域震后形变特征进行了初步分析。     

GRACE重力卫星在陆地水储量变化监测中的应用

介绍了GRACE重力卫星,并对GRACE重力卫星数据在陆地水储量变化中的应用现状进行分析,总结了GRACE重力卫星数据在陆地水储量变化检测中的数据获取、计算方法和精度分析,以及Grace数据在不同区域尺度陆地水储量变化估算中的应用情况,最后,指出GRACE在水储量应用中的不足和未来的研究方向。     

利用GOCE和GRACE卫星观测数据确定静态重力场模型

高精度静态卫星重力场模型在全球海洋环流研究、全球/区域数字高程基准面确定等领域有重要应用,本文研究仅利用GOCE卫星和联合GRACE卫星观测数据确定高精度高阶次静态重力场模型.利用GOCE卫星全周期高精度引力梯度分量(V_xx、V_yy、V_zz和V_xz)观测值基于直接最小二乘法构建300阶次的SGG(Satellite Gravity Gradiometry)法方程,并利用卫星跟踪卫星观测值基于点域加速度法构建130阶SST(Satellite-to-Satellite Tracking)法方程,然后利用方差分量估计联合SGG和SST法方程确定300阶次纯GOCE卫星重力场模型GOSG02S.利用全周期GRACE观测数据由动力学方法解算了180阶次的SWPU-GRACE2021S模型,并将其对应法方程与GOCE卫星法方程联合解算了GRACE和GOCE的联合模型WHU-SWPU-GOGR2022S.分别基于XGM2019模型和GPS水准数据对本文解算的三个模型GOSG02S、SWPU-GRACE2021S和WHU-SWPU-GOGR2022S在频域和空域进行了精度分析,结果表明,GOSG02S和WHU-SWPU-GOGR2022S模型与GO_CONS_GCF_2_DIR_R6、GO_CONS_GCF_2_TIM_R6、GO_CONS_GCF_2_SPW_R5、GOCO06s和Tongji-GMMG2021S等使用了GOCE卫星全周期数据的模型精度相当,精度差异基本都在毫米量级;SWPU-GRACE2021S模型在160阶次之前与国际主流GRACE卫星重力场模型ITSG-Grace2018s和Tongji-Grace02s精度相当.

     

利用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水文模型的相关性略小。