多通道奇异谱组合滤波反演水储量变化The inversion of terrestrial water storage changes by means of multi-channel singular spectrum combnation filter
李武东,常晓涛,朱广彬,郭金运
摘要(Abstract)：
针对GRACE时变重力场球谐系数的高阶次项中含有的较高噪声会造成反潢结果产生误差的问题,该文采用多通道奇异谱分析和高斯滤波相结合的方案来滤除GRACE数据中的噪声,并与普遍使用的去相关滤波和高斯滤波相结合的方案进行了对比验证。首先利用这两种方法分别计算了2007年4、10月份的全球水储量变化,结果显示,这两种方法反演的结果基本相同,而多通道奇异谱分析与高斯滤波相组合的方法提高了信噪比,减弱了信号泄露误差,验证了该组合方案的可行性。然后又分别计算了亚马逊流域2004—2010年的水储量时间变化序列,并与GLDAS水文模型进行了对比验证,结果显示,这两种方法得到的水储量时间序列变化趋势基本相同,各自之间具有强相关性,进一步说明了该文方案在GRACE数据滤波中的可行性。
关键词(KeyWords)： GRACE;多通道奇异谱分析;去相关滤波;水文模型;陆地水储量变化
基金项目(Foundation): 国家自然科学基金项目(41204007);; 国家高技术研究发展计划项目(2013AA122502);; 国家重点基础研究发展计划项目(2013CB733302)
作者(Author): 李武东,常晓涛,朱广彬,郭金运
DOI: 10.16251/j.cnki.1009-2307.2015.12.004
参考文献(References)：
- [1]詹金刚,王勇,郝晓光.GRACE时变重力位系数误差的改进去相关算法[J].测绘学报,2011,40(4):442-446.
- [2]翟宁,王泽民,伍岳,等.利用GRACE反演长江流域水储量变化[J].武汉大学学报·信息科学版,2009,34(4):436-439.
- [3]CHEN J L,WILSON C R,SEO K W.Optimized smoothing of Gravity Recovery and Climate Experiment(GRACE)time-variable gravity observations[J].Journal of Geophysical Research,2006,111(B6):466-480.
- [4]SWENSON S,WAHR J.Post-processing removal of correlated errors in GRACE data[J].Geophysical Research Letters,2006,33(8).
- [5]RANGELOVA E,VAN DER WAL W,BRAUN A,et al.Analysis of Gravity Recovery and Climate Experiment time-variable MSSA redistribution signals over North America by means of principal component analysis[J].Journal of Geophysical Research,2007,112(F3):256-260.
- [6]穆大鹏,郭金运,孙中昶,等.基于主成分分析的GRACE重力场模型等效水高[J].地球物理学进展,2014,29(4):1512-1517.
- [7]GUO J,MU D,LIU X,et al.Equivalent water height extracted from GRACE gravity field model with robust independent component analysis[J].Acta Geophysica,2014,62(4):953-972.
- [8]WANG X,DE LINAGE C,FAMIGLIETTI J,et al.Gravity Recovery and Climate Experiment(GRACE)detection of water storage changes in the Three Gorges Reservoir of China and comparison with in situ measurements[J].Water Resources Research,2011,47(12):1091-1096.
- [9]OROPEZA V,SACCHI M.Simultaneous seismic data denoising and reconstruction via multichannel singular spectrum analysis[J].Geophysics,2011,76(3):V25-V32.
- [10]ZOTOV L V,SHUM C K.Multichannel singular spectrum analysis of the gravity field data from GRACE satellites[C]//ASTROPHYSICS AND COSMOLOGY AFTER GAMOW:Proceedings of the 4th Gamow International Conference on Astrophysics and Cosmology After Gamow and the 9th Gamow Summer School“Astronomy and Beyond:Astrophysics,Cosmology,Radio Astronomy,High Energy Physics and Astrobiology”.AIP Publishing,2010,1206(1):473-479.
- [11]RANGELOVA E,SIDERIS M G,KIM J W.On the capabilities of the multi-channel singular spectrum method for extracting the main periodic and non-periodic variability from weekly GRACE data[J].Journal of Geodynamics,2012,54(2):64-78.
- [12]RANGELOVA E,VAN DER WAL W,SIDERIS M G,et al.Spatiotemporal analysis of the GRACE-derived MSSA variations in North America by means of multichannel singular spectrum analysis[J].Geoid&Earth Observation,2010,135(7):539-546.
- [13]WAHR J,MOLENAAR M,BRYAN F.Time variability of the Earth’s gravity field:hydrological and oceanic effects and their possible detection using GRACE[J].Journal of Geophysical Research Atmospheres,1998,103(B12):30205-30230.
- [14]CHENG M,TAPLEY B D.Variations in the Earth’s oblateness during the past 28years[J].Journal of Geophysical Research,2004,109(B9).
- [15]PELTIER W R.Global glacial isostasy and the surface of the ice-age Earth:the ICE-5G(VM2)model and GRACE[J].Annu.Rev.Earth Planet.Sci.,2004,32(12):111-149.
- [16]RODELL M,HOUSER P R,JAMBOR U,et al.The global land data assimilation system[J].Bulletin of the American Meteorological Society,2004,85(3):381-394.
- [17]HASSANI H.Singular spectrum analysis:methodology and comparison[J].Journal of Data Science,2007,5(2):239-257.
- [18]GOLYANDINA N,NEKRUTKIN V,ZHIGLJAVSKY A A.Analysis of time series structure:SSA and related techniques[M].CRC Press,2001.