重力场恢复与气候实验的石羊河流域水储量变化

为了更好地从空间尺度显示石羊河流域水储量的变化情况,该文利用2003年1月至2014年12月GRACE RL06重力场数据计算其陆地水储量变化,对数据进行预处理以及去相关滤波和300km高斯平滑处理,并利用流域内4个气象站的降水平均数据与GRACE、GLDAS的结果进行对比分析。结果表明,GRACE反演结果与GLDAS水文模型变化趋势基本一致,时间上呈现明显季节性变化,研究区水储量整体呈现0.93mm/a的上升趋势,期间水储量变化波动较大主要与季节性降水有关。     

Seasonal global mean sea level change from satellite altimeter, GRACE, and geophysical models

We estimate seasonal global mean sea level changes using different data resources, including sea level anomalies from satellite radar altimetry, ocean temperature and salinity from the World Ocean Atlas 2001, time-variable gravity observations from the Gravity Recovery and Climate Experiment (GRACE) mission, and terrestrial water storage and atmospheric water vapor changes from the NASA global land data assimilation system and National Centers for Environmental Prediction reanalysis atmospheric model. The results from all estimates are consistent in amplitude and phase at the annual period, in some cases with remarkably good agreement. The results provide a good measure of average annual variation of water stored within atmospheric, land, and ocean reservoirs. We examine how varied treatments of degree-2 and degree-1 spherical harmonics from GRACE, laser ranging, and Earth rotation variations affect GRACE mean sea level change estimates. We also show that correcting the standard equilibrium ocean pole tide correction for mass conservation is needed when using satellite altimeter data in global mean sea level studies. These encouraging results indicate that is reasonable to consider estimating longer-term time series of water storage in these reservoirs, as a way of tracking climate change.     

GRACE反演南极冰盖质量变化的高斯与Wiener滤波比较

本文利用UTCSR 2003年1月到2008年8月间的GRACE Level-2 RL04重力场模型估计了南极冰盖质量变化。计算过程中分别采用高斯和Wiener滤波两种平滑方法,分别采用22、43和65个月重力场模型计算Wiener滤波信号与噪声函数,得出以下结论:在实际的计算过程中需要具体计算Wiener滤波平滑因子值,65个月GRACE重力场模型计算得到的Wiener滤波权值非常接近于平滑半径为540km高斯滤波权值;采用两种不同的滤波方法在相同区域质量变化率基本相同。     

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

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

利用卫星观测数据评估GLDAS与WGHM水文模型的适用性

利用2002—2012年的GLDAS和WGHM模型模拟水文产品,以及重力恢复与气候试验卫星（Gravity Recovery and Climate Experiment,GRACE）观测数据,计算了全球范围内30个主要流域的水储量变化时间序列,从模拟数据与观测数据的年周期振幅、长期趋势项及时空分布一致性等几个方面,对GLDAS和WGHM进行了评估。结果表明,GLDAS的4个子模型都表现出了明显的季节性变化,CLM年周期振幅输出最小,MOSAIC和VIC最大,NOAH居中,且最接近4个子模型的平均值。与GRACE结果相比,约80%流域的GLDAS与WGHM模型年周期振幅输出呈明显低估现象,且GLDAS的低估程度大于WGHM,但靠近北极高纬度地区的流域有相反的情况出现。在长期趋势项方面,三者结果差异较大,尤其是对于面积较小且人类活动影响较大的流域,GLDAS与WGHM模型不能充分反映人类活动的影响,模型输出表现较差,GRACE结果更接近实际情况。此外,还研究了流域水储量长期变化趋势与灌溉率的关系,发现呈现明显下降趋势的流域主要集中在高灌溉率（>10%）地区,而灌溉率是影响流域水储量变化的重要因素之一。     

日本Mw9.0地震前GRACE卫星重力变化

利用GRACE重力卫星月重力场数据,通过去相关与高斯滤波等方法获取日本大地震震前震源区周缘年度、季度和年度差分重力变化,以及若干点位时间序列重力变化;分析结果表明:日本大地震前5年内在震源区周边出现了比较明显的卫星重力异常正负交替和迁移现象,至震前1－2年,震区周边形成了明显正负异常区,正重力异常区重力增加现象明显;点位重力时间序列分布指出了日本MW9.0级地震前存在与1976年唐山地震类似的重力变化现象。这些结果反映了日本大地震震前震源区周边地下物质运动、质量迁移和能量积累等问题,为研究该地震的孕育过程提供了证据。     

Simulated estimation of hydrological loads from GRACE

Four different basin functions are developed to estimate water storage variations within individual river basins from time variations in the Stokes coefficients now available from the GRACE mission. The four basin functions are evaluated using simulated data. Basin functions differ in how they minimize effects of three major error sources: measurement error; leakage of signal from one region to another; and errors in the atmospheric pressure field removed during GRACE data processing. Three of the basin functions are constant in time, while the fourth changes monthly using information about the signal (hydrologic and oceanic load variations). To test basin functions performance, Stokes coefficient variations from land and ocean models are synthesized, and error levels 50 and 100 times greater than pre-launch GRACE error estimate are used to corrupt them. Errors at 50 times pre-launch estimates approximately simulate current GRACE data. GRACE recovery of water storage variations is attempted for five different river basins (Amazon, Mississippi, Lena, Huang He and Oranje), representing a variety of sizes, locations, and signal variance. In the large basins (Amazon, Mississippi and Lena), water storage variations are recovered successfully at both error levels. As the error level increases from 50 to 100 times, basin functions change their shape, yielding less atmospheric pressure error and more leakage error. Amplitude spectra of measurement and atmospheric pressure errors have different shapes, but the best results are obtained when both are used in basin function design. When high-quality information about the signal is available, for example from climate and ocean models, changing the basin function each month can reduce leakage error and improve estimates of time variable water storage within basins.     

利用GPS和GRACE研究澳大利亚地壳垂向季节性变化

为探究重力场恢复与气候实验(gravity recovery and climate experiment,GRACE)卫星与全球定位系统(global positioning system,GPS)两种独立技术获取的因陆地水储量变化引起的地壳垂向季节性位移的一致性,选取澳大利亚27个GPS站点5～10 a的高程时间序...     

Annual variations in water storage and precipitation in the Amazon Basin

We combine satellite gravity data from the gravity recovery and climate experiment (GRACE) and precipitation measurements from the National Oceanic and Atmospheric Administration (NOAA) Climate Prediction Center’s (CPC) Merged Analysis of Precipitation (CMAP) and the Tropical Rainfall Measuring Mission (TRMM), over the period from mid-2002 to mid-2006, to investigate the relative importance of sink (runoff and evaporation) and source (precipitation) terms in the hydrological balance of the Amazon Basin. When linear and quadratic terms are removed, the time-series of land water storage variations estimated from GRACE exhibits a dominant annual signal of 250 mm peak-to-peak, which is equivalent to a water volume change of ~1,800 km³. A comparison of this trend with accumulated (i.e., integrated) precipitation shows excellent agreement and no evidence of basin saturation. The agreement indicates that the net runoff and evaporation contributes significantly less than precipitation to the annual hydrological mass balance. Indeed, raw residuals between the de-trended water storage and precipitation anomalies range from ±40 mm. This range is consistent with stream-flow measurements from the region, although the latter are characterized by a stronger annual signal than our residuals, suggesting that runoff and evaporation may act to partially cancel each other.     

利用GPS垂直位移反演云南省陆地水储量变化

地表质量的重分布会引起固体地球的弹性形变,GPS连续运行观测站能够精确测定地表负荷引起的地壳形变。本文通过模拟数据对利用云南省及其周边47个中国大陆构造环境监测网（陆态网）台站反演云南地区陆地水储量的可行性进行分析：以水文模型周年振幅为真值,计算47个台站点的负荷形变,同时加入随机误差构成模拟观测数据,最后采用模型反演陆地水储量变化;1000次的随机模拟试验表明利用当前GPS台站数据可有效地反演云南地区陆地水储量变化。基于上述结论,笔者反演了云南省2010—2014年陆地水储量变化,GPS反演结果表明：云南省陆地水变化呈现明显的地域分布特征,西南部高山地区的水储量周年变化高于东部平原地区;在时间尺度上,云南省大部分地区水储量在10月（夏季末）达到最大值,在4月（冬季末）达到最小值;云南省2010—2014年陆地水呈缓慢增长趋势,约为20 mm/a。通过GPS陆地水储量反演结果与GRACE、GLDAS以及TRMM数据综合对比分析,表明利用云南地区当前GPS台站可以作为独立观测量用于GRACE与GRACE Follow-on衔接期间的陆地水储量变化监测。     

不同模型的地表质量异常一阶项、二阶项估计

为了合理补充重力场恢复与气候试验卫星（Gravity Recovery and Climate Experiment,GRACE）时变重力场的一阶斯托克斯系数（C₁₀、C₁₁、S₁₁）和替换二阶斯托克斯系数（C₂₀）,介绍了相关GRACE-OBP算法及其改进的算法,比较了相应的Chamber Model和4个Sun Model的一阶系数及其计算的地表质量异常,同时比较了基于卫星激光测距观测的Cheng Model与4个Sun Model的C₂₀及其地表质量异常。结果表明,GRACE-OBP算法的一阶系数、卫星激光测距观测的C₂₀及其地表质量异常与改进的GRACE-OBP算法在趋势项上有很大差异,但周年项差异相对较小。利用不同截断阶数和不同机构的GRACE时变重力场模型,对其趋势项和周年项都有一定影响,且对趋势项影响更大。因此,在计算陆地水储量变化时,建议使用改进的GRACE-OBP算法的估计结果,使用较理想的、截断阶数较高的GRACE时变重力模型。     

利用ICA方法提取奥卡万戈三角洲水储量变化信号

利用重力恢复与气候实验(gravity recovery and climate experiment, GRACE)时变地球重力场模型计算得到非洲奥卡万戈三角洲地区2003-01—2014-12的陆地水储量变化信息,分别采用主成分分析(principal component analysis, PCA)和独立成分分析(independent component analysis, ICA)提取质量变化信号,并与全球陆地数据同化系统(global land data assimilation system, GLDAS)的水文模型进行对比。结果显示,在奥卡万戈河流域东北部,水储量表现出很强的周期性变化,两种数据空间特征分布的信号出现在相同位置的成分GRACE-IC1和GLDAS-IC1对应的时间序列的相关系数达到0.85。奥卡万戈三角洲地区水储量从2003-01—2011-10呈现上升趋势,两种数据空间特征分布的信号出现在相同位置的成分GRACE-IC2和GLDAS-IC3对应的时间序列的相关系数达到0.81,说明GRACE反演结果与GLDAS水文模型反演结果在研究区域内具有很强的一致...     

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

华北地区长期重力变化监测

利用GRACE卫星重力资料,计算了华北地区的长期重力变化结果,利用6个测站的绝对重力观测资料,获取了测站的重力变化时间序列,同时获取了北京、泰安测站的GRACE卫星月重力变化时间序列。卫星重力观测结果显示华北地区地下水流失严重,绝对重力观测结果表明地面沉降严重。     

Approximate decorrelation and non-isotropic smoothing of time-variable GRACE-type gravity field models

We discuss a new method for approximately decorrelating and non-isotropically filtering the monthly gravity fields provided by the gravity recovery and climate experiment (GRACE) twin-satellite mission. The procedure is more efficient than conventional Gaussian-type isotropic filters in reducing stripes and spurious patterns, while retaining the signal magnitudes. One of the problems that users of GRACE level 2 monthly gravity field solutions fight is the effect of increasing noise in higher frequencies. Simply truncating the spherical harmonic solution at low degrees causes the loss of a significant portion of signal, which is not an option if one is interested in geophysical phenomena on a scale of few hundred to few thousand km. The common approach is to filter the published solutions, that is to convolve them with an isotropic kernel that allows an interpretation as smoothed averaging. The downside of this approach is an amplitude bias and the fact that it neither accounts for the variable data density that increases towards the poles where the orbits converge nor for the anisotropic error correlation structure that the solutions exhibit. Here a relatively simple regularization procedure will be outlined, which allows one to take the latter two effects into account, on the basis of published level 2 products. This leads to a series of approximate decorrelation transformations applied to the monthly solutions, which enable a successive smoothing to reduce the noise in the higher frequencies. This smoothing effect may be used to generate solutions that behave, on average over all possible directions, very close to Gaussian-type filtered ones. The localizing and smoothing properties of our non-isotropic kernels are compared with Gaussian kernels in terms of the kernel variance and the resulting amplitude bias for a standard signal. Examples involving real GRACE level 2 fields as well as geophysical models are used to demonstrate the techniques. With the new method, we find that the characteristic striping pattern in the GRACE solutions are much more reduced than Gaussian-filtered solutions of comparable signal amplitude and root mean square.     

顾及GRACE季节影响的华北平原水储量变化反演

基于移去-恢复方法,利用GRACE卫星重力数据研究了华北平原2003年1月—2014年6月的陆地水变化,提出了一种顾及季节影响尺度因子计算方法,用于减小GRACE后处理误差。将本文方法所得结果与水文模式和降水模型结合进行比较分析,结果表明：2008年之前华北平原发生比较明显的干旱现象,陆地水和地表水分别以（-7.9±2.4）mm/a和（-7.3±2.8）mm/a的速度下降,但地下水仅以（-0.6±1.4）mm/a的速度减少;2008年之后,陆地水和地表水分别以（4.3±1.3）mm/a和（10.9±2.1）mm/a的速度上升,期间地下水的超采严重,平均以（-6.5±1.2）mm/a的速度下降,但下降速度以0.9 mm/a²的趋势减缓。研究期间内,华北平原的陆地水、地表水和地下水整体分别呈现（-2.0±0.6）、（2.9±0.7）和（-4.8±0.7）mm/a的变化趋势。最后,利用交叉小波谱分析了GRACE滞后于降水的现象。研究结论表明,降水和地下水开采是影响华北平原水储量变化的两大决定性因素。     

温度变化对中国大陆三维周年位移的影响

除了地表质量重新分布外,地表温度变化是影响地表周年变化的另一重要成因。本文利用全球温度变化数据,基于三维全空间热弹性形变模型,计算温度变化在中国大陆引起的地表热弹性形变,并讨论它对中国大陆三维周年位移的影响。结果表明,温度变化引起的地表周年变化振幅在毫米量级。中国大陆构造环境监测网络(简称:陆态网络)GPS台站受地表温度变化影响最大的台站是HLAR(海拉尔),东向、北向以及垂向的周年振幅矢量和为~2.293mm;影响最小的台站是HIYS(永暑礁),东向、北向以及垂向的周年振幅矢量和为~0.177mm。为了说明温度变化对地表周年形变的影响,本文联合GRACE以及物质负荷模型(MODEL)研究中国大陆地表三维周年位移。考虑温度变化后的MODEL、GRACE获取的地表周年形变在东向、北向、垂向的周年信号分别改进了6%、6%、2%;16%、5%、15%。结果表明,温度变化是物质负荷以外引起大陆地表形变的重要因素。     

根据人卫激光测距、GRACE和地球物理模型求解地球低阶重力场季节变化

本文利用2002年4月至2010年10月的Lageos1和Lageos2两颗激光卫星观测数据、GRACE以及地球物理模型三种独立的方法计算地球低阶重力场系数J 的变化,根据大气压强数据计算 J 时分别按反变气压计（IB）和非反变气压计(NIB)两种假设进行计算。通过分析 J 的季节特性表明,大气在NIB假设下得到的周年振幅比在IB假设下得到的振幅大3倍左右,相位相差47°;大气和陆地水的质量变化对 J 周年变化的贡献占主导地位,海洋的影响最小;大气、海洋和陆地水得到 J 半年振幅和相位值与SLR得到的振幅和相位值吻合较差,尤其是在IB假设下大气得到的结果与SLR结果相差最大; SLR、GRACE和地球物理模型三种独立方法得到的 J 周年项之间吻合相对较好,GRACE得到的周年振幅比SLR得到的周年振幅大50%左右, SLR观测得到的 J 周年振幅介于在NIB和IB两种假设下地球物理模型得到的结果之间;GRACE与SLR得到的 J 半年项的振幅相同,在IB假设下AOW得到的 J 半年振幅和相位与SLR结果差异最大。     

联合GRACE和GLDAS数据的河南省地下水储量反演与时空变化分析

针对利用传统监测手段难以高效获取地下水储量观测数据的问题,基于GRACE重力卫星的大尺度水资源储量反演已成为当前水资源调查的研究热点。本文利用2012-2016年CSR机构发布的GRACE RL06月解数据,通过等效水高反演得到河南省陆地水储量时序结果,扣除由同期GLDAS水文模型计算得到的地表水储量时序数据,从而得到河南省地下水储量时序变化数据结果。经与地下水位监测井实测数据进行对比验证,相关系数显著性水平达0.01,表明本文算法流程具有较高的可靠性。进一步的统计分析结果表明,河南省北部地区的地下水储量呈亏损态势,最大变化率超过26 mm/a;河南省中部和东部地区地下水储量有一定盈余,最大变化率超过16 mm/a,相关结果数据与河南省水利局公布的全省主要地下水超采区范围吻合。本文旨在利用GRACE重力卫星数据与GLDAS水文模型反演获取河南省地下水储量空间分布差异及演变趋势,相关算法流程可为广域地下水储量调查监测提供技术支撑;研究数据可为该区域地下水资源的合理利用与保护提供参考。     

Combination of temporal gravity variations resulting from superconducting gravimeter (SG) recordings, GRACE satellite observations and global hydrology models

Gravity recovery and climate experiment (GRACE)-derived temporal gravity variations can be resolved within the μgal (10^?8 m/s ²) range, if we restrict the spatial resolution to a half-wavelength of about 1,500 km and the temporal resolution to 1 month. For independent validations, a comparison with ground gravity measurements is of fundamental interest. For this purpose, data from selected superconducting gravimeter (SG) stations forming the Global Geodynamics Project (GGP) network are used. For comparison, GRACE and SG data sets are reduced for the same known gravity effects due to Earth and ocean tides, pole tide and atmosphere. In contrast to GRACE, the SG also measures gravity changes due to load-induced height variations, whereas the satellite-derived models do not contain this effect. For a solid spherical harmonic decomposition of the gravity field, this load effect can be modelled using degree-dependent load Love numbers, and this effect is added to the satellite-derived models. After reduction of the known gravity effects from both data sets, the remaining part can mainly be assumed to represent mass changes in terrestrial water storage. Therefore, gravity variations derived from global hydrological models are applied to verify the SG and GRACE results. Conversely, the hydrology models can be checked by gravity variations determined from GRACE and SG observations. Such a comparison shows quite a good agreement between gravity variation derived from SG, GRACE and hydrology models, which lie within their estimated error limits for most of the studied SG locations. It is shown that the SG gravity variations (point measurements) are representative for a large area within the accuracy, if local gravity effects are removed. The individual discrepancies between SG, GRACE and hydrology models may give hints for further investigations of each data series.