将OpenMP并行运算库和Intel Math Kernel Library10.2科学计算库运用到高阶地球重力场反演,显著提高了计算效率。模拟结果表明:1)在以单历元为解算单元形成子法方程系数矩阵时,OpenMP可下三角压缩存储,其内存销耗和时间销耗均比较小;2)当利用多个历元组成高维系数矩阵,然后再形成法方程时,MKL算法才能体现出高效性能;3)MKL求逆算法效率远高于OpenMP算法;4)综合利用OpenMP和MKL算法的优势,可显著提高高阶地球重力场反演的效率。
Monitoring glacier mass balance is crucial to managing water resources and also to understanding climate change for the arid and semi-arid regions of Central Asia. This study extracted the inter-annual oscillations of glacier mass over Central Asia from the first ten principal components(S-PCs) of filtered variability via multichannel singular spectral analysis(MSSA), based on gridded data of glacier mass inferred from Gravity Recovery and Climate Experiment(GRACE) data obtained from July 2002 to March 2015. Two significant cycles of glacier mass balance oscillations were identified. The first cycle with a period of 6.1-year accounted for 54.5% of the total variance and the second with a period of 2.3-year accounted for 4.3%. The 6.1-year oscillation exhibited a stronger variability compared with the 2.3-year oscillation. For the 6.1-year oscillation, the results from lagged cross-correlation function suggested that there were significant correlations between glacier mass balances and precipitation variations with the precipitation variations leading the response of glacier mass balances by 9–16 months.
We use the average crustal structure of the CRUST1.0 model for the Tibetan Plateau to establish a realistic earth model termed as TC1 P, and data from the Global Land Data Assimilation System(GLDAS) hydrology model and Gravity Recovery and Climate Experiment(GRACE) data, to generate the hydrology signals assumed in this study. Modeling of surface radial displacements and gravity variation is performed using both TC1 P and the global Preliminary Reference Earth Model(PREM). Furthermore, inversions of the hydrology signals based on simulated Global Positioning System(GPS) and GRACE data are performed using PREM. Results show that crust in TC1 P is harder and softer than that in PREM above and below a depth of 15 km, respectively, causing larger differences in the computed load Love numbers and loading Green’s functions. When annual hydrology signals are assumed,the differences of the radial displacements are found to be as large as approximately0.6 mm for the truncated degree of 180; while for hydrology-trend signals the differences are very small. When annual hydrology signals and the trends are assumed, the differences in the surface gravity variation are very small. It is considered that TC1 P can be used to efficiently remove the hydrological effects on the monitoring of crustal movement. It was also found that when PREM is used inappropriately, the inversion of the hydrology signals from simulated annual GPS signals can only recover approximately 88.0% of the annual hydrology signals for the truncated degree of 180, and the inversion of hydrology signals from the simulated trend GPS signals can recover approximately 92.5% for the truncated degree of 90. However, when using the simulated GRACE data, it is possible to recover almost 100%. Therefore, in future, the TC1 P model can be used in the inversions ofhydrology signals based on GPS network data. PREM is also valid for use with inversions of hydrology signals from GRACE data at resolutions of approximately 220 km and larger.
Wang HanshengXiang LongweiWu PatrickJia LuluJiang LimingShen QiangSteffen Holger
As global warming continues,the monitoring of changes in terrestrial water storage becomes increasingly important since it plays a critical role in understanding global change and water resource management.In North America as elsewhere in the world,changes in water resources strongly impact agriculture and animal husbandry.From a combination of Gravity Recovery and Climate Experiment(GRACE) gravity and Global Positioning System(GPS) data,it is recently found that water storage from August,2002 to March,2011 recovered after the extreme Canadian Prairies drought between 1999 and 2005.In this paper,we use GRACE monthly gravity data of Release 5 to track the water storage change from August,2002 to June,2014.In Canadian Prairies and the Great Lakes areas,the total water storage is found to have increased during the last decade by a rate of 73.8 ± 14.5 Gt/a,which is larger than that found in the previous study due to the longer time span of GRACE observations used and the reduction of the leakage error.We also find a long term decrease of water storage at a rate of-12.0 ± 4.2 Gt/a in Ungava Peninsula,possibly due to permafrost degradation and less snow accumulation during the winter in the region.In addition,the effect of total mass gain in the surveyed area,on present-day sea level,amounts to-0.18 mm/a,and thus should be taken into account in studies of global sea level change.
Wang HanshengXiang LongweiJia LuluWu PatrickSteffen HolgerJiang LimingShen Qiang