Monitoring of Crustal Vertical Deformation and Gravity Change Caused by Environmental Load in Wenzhou-Lishui Region Using CORS Network
-
摘要: 基于CORS(continuously operating reference stations)站网综合解算,采用负荷场移去恢复技术,研究了温州-丽水区域环境负荷引起的地壳垂直形变与地面重力时空变化,并与GRACE (gravity recovery and climate experiment)重力卫星结果进行比较。研究发现:①环境负荷对地壳垂直形变与重力变化的影响分别达到厘米和十几微伽量级,季节性变化特征显著;②冬季地壳垂直形变与重力变化分别呈现由西至东依次递减与增加趋势;③与GRACE结果相比,除局部区域存在差异,整体上时空变化趋势具有较高一致性;④CORS站网能够监测地表环境负荷引起的地壳垂直形变与重力场时空变化。该研究结果可为环境动力学研究与地质灾害监测预警提供重要参考。Abstract: Based on the comprehensive calculation of CORS network, this paper uses the load field removal-recovery technology to study the crustal vertical deformation and ground gravity spatio-temporal changes caused by environmental load in Wenzhou-Lishui region and compares it with GRACE gravity satellite results. Some conclusions are drown as follows: ① The effects of environmental load on the vertical deformation and gravity change of crust are up to the centimeter and ten micro-gauge levels, respectively, and the seasonal variation characteristics are significant. ② The vertical deformation and gravity change in winter show a decreasing and increasing trend from west to east, respectively. ③ Compared with GRACE results, in addition to local region differences, the overall trend of spatiotemporal changes has a higher consistency. ④The CORS network can monitor the spatiotemporal changes of crustal vertical deformation and gravity field caused by total surface environmental load. The research methods and results of this paper can provide an important reference for environmental dynamics research and geological disaster monitoring and early warning.
-
Keywords:
- CORS network /
- environmental load /
- vertical deformation /
- gravity change
-
火星是人类深空探测器到访最多的行星。火星快车(Mars Express,MEX)是欧空局(European Space Agency,ESA)的首颗火星探测器[1],在轨运行期间传回的大量火星地表影像等资料表明,火星大气层内存在甲烷[2],火星南极冠存在水冰。MEX还多次飞越火卫一(Phobos),测量了其质量和低阶重力场等[3]。国内外开发了许多火星探测器定轨软件[4-6],如武汉大学深空探测器精密定轨与重力场解算软件系统中的火星重力场解算和分析系统模块(Mars Gravity Recovery and Analysis Software/System,MAGREAS)对MEX的定轨结果已经达到了比利时皇家天文台发布的精密星历的精度水平[7-8],本文将选用MAGREAS作为MEX的定轨软件。
目前定轨大多采用原始数据,其中含有较强的噪声,严重影响定轨的精度。噪声也是影响各类接收机性能的主要因素之一[9]。噪声可以分为热噪声、散弹噪声和闪烁噪声。热噪声是电子设备中一种十分常见的白噪声[10],也是本文主要考虑消除的噪声成分。深空通讯设备中的热噪声主要来源于地面发射器/接收机和星载转发器。由MEX 1 s采样间隔获取的X波段原始数据提取出的航天器速度的噪声为0.05~0.37 mm/s,具体值取决于太阳-火星-地球的夹角,平均为0.13 mm/s[3]。为达到更好的定轨效果,有研究者采用滤波器处理原始数据,提高数据信噪比。Andert采用Kaiser窗低通滤波器处理了MEX 2006年3月和2008年7月两次飞掠Phobos的频率数据,结果表明,两次飞掠Phobos的频率数据的残差标准差分别从6.56 mHz减少至1.71 mHz、从7.32 mHz减少至1.97 mHz[11]。为了除去原位相位特征的影响,并精确保留月球重力场的信息,Liu等利用Kaiser窗设计了一个低通滤波器, 对日本“月亮女神”月球探测器的两颗子卫星的原始多普勒数据进行滤波处理[12]。滤波后,Vstar的双程多普勒测量残差均方根(root mean square, RMS)从0.355 Hz减小到0.001 6 Hz; Rstar的双程多普勒测量残差RMS从0.293 Hz减小到0.002 1 Hz,四程(Rstar-主卫星)多普勒测量残差RMS从1.028 Hz减小到0.025 Hz。
基于零相位分析,本文比较了FRR(forward-filter reverse-filter reverse-output)、RRF(reverse-filter reverse-filter forward-output)和Matlab中的filtfilt这3种滤波器的优劣,并设计了一种有效的零相位Kaiser窗低通滤波器。
1 零相位滤波算法
1.1 滤波器的基本概念
滤波器的传递函数H(ejω)可以用极坐标表示为:
$$ H({{\rm{e}}^{j\omega }}) = \left| {H\left( \omega \right)} \right|\cdot{{\rm{e}}^{j\varphi (\omega )}} $$ (1) 式中,|H(ω)|、φ(ω)分别为滤波器的振幅响应和相位响应, 计算公式分别为:
$$ \left| {H\left( \omega \right)} \right| = \sqrt {{\rm{R}}{{\rm{e}}^2}\left[ {H({{\rm{e}}^{j\omega }})} \right] + {\rm{I}}{{\rm{m}}^2}[H({{\rm{e}}^{j\omega }})]} $$ (2) $$ \varphi \left( \omega \right) = {\rm{arctan}}\frac{{{\rm{Im}}[H({{\rm{e}}^{j\omega }})]}}{{{\rm{Re}}[H({{\rm{e}}^{j\omega }})]}} $$ (3) 式中,Re(·)为函数实部;Im(·)为函数虚部;j为虚数单位;ω为数字频率。
根据滤波器的幅度响应,可以将滤波器分为低通、高通、通带和阻带4类滤波器[13],图 1为不同类型滤波器的理想幅度响应。
图 1所示的滤波器幅度响应是一种理想状态,在现实中不可能实现,只能尽可能地逼近这种状态。在逼近的过程中,滤波器的不同性能指标通常不可能同时达到最优。因此在设计滤波器时,一般根据实际需要,允许滤波器在通带和阻带内与理想状态有一定的偏差,通带与阻带之间也允许有一个过渡带。以低通滤波器为例,滤波器的实际幅频特性如图 2所示。
图 2中,δp、δs分别称为通带波纹和阻带波纹;ωp、ωs分别称为通带截止频率和阻带截止频率;Δω=ωs-ωp称为过渡带。
设计滤波器时,一般情况下,振幅特性由给定的通带和阻带衰减确定,衰减A(ω)用反映功率增益的幅度平方函数(或称作模平方函数)|H(ω)|2来定义:
$$ A\left( \omega \right) = - 10{\rm{lg}}{\left| {H\left( \omega \right)} \right|^2} = - 20{\rm{lg}}\left| {H\left( \omega \right)} \right| $$ (4) 所以通带衰减Ap和阻带衰减As可以表示为:
$$ {A_p} = - 20{\rm{lg}}(1 - {\delta _p}) $$ (5) $$ {A_s} = - 20{\rm{lg}}{\delta _s} $$ (6) 1.2 零相位滤波
一个信号经过一个滤波器系统后,会将信号每个频率分量的振幅乘上系统振幅响应的模,以改变信号不同频率成分的能量,实现噪声的滤除。滤波器系统在改变信号幅频性质的同时也会在原信号相位上附加一个相位,称为系统的相移。如果这种相位的改变不是所预期的,就会造成相位的失真,影响数据质量[14]。
对多普勒数据进行滤波时不希望相位发生变化,零相位滤波器具有零相位系统特性,可以获得精确零相位失真的信号[15]。零相位滤波可以采用FRR滤波方法,首先将输入序列按顺序滤波,然后将得到的结果逆转后再滤波,最后将所得结果逆转后输出,即可得到精确零相位失真的序列。
FRR滤波的时域描述可以表示为:
$$ \left\{ \begin{array}{l} {y_1}\left( n \right) = x\left( n \right)\cdot h\left( n \right)\\ {y_2}\left( n \right) = {y_1}({L_1} - n + 1)\\ {y_3}\left( n \right) = {y_2}\left( n \right)\cdot h\left( n \right)\\ {y_4}\left( n \right) = {y_3}({L_1} - n + 1) \end{array} \right. $$ (7) 式中,L1为序列长度;n∈[1, L1];x(n)表示输入序列;h(n)表示数据滤波器冲击响应序列;y(n)表示滤波或者序列逆转后的结果。
FRR滤波的频率描述即为式(7)相应的频域表示:
$$ \left\{ \begin{array}{l} {Y_1}({{\rm{e}}^{j\omega }}) = X({{\rm{e}}^{j\omega }})\cdot H({{\rm{e}}^{j\omega }})\\ {Y_2}({{\rm{e}}^{j\omega }}) = {{\rm{e}}^{ - j\omega (N + 1)}}\cdot{Y_1}({{\rm{e}}^{ - j\omega }})\\ {Y_3}({{\rm{e}}^{j\omega }}) = {Y_2}({{\rm{e}}^{j\omega }})\cdot H({{\rm{e}}^{j\omega }})\\ {Y_4}({{\rm{e}}^{j\omega }}) = {{\rm{e}}^{ - j\omega (N + 1)}}\cdot{Y_3}({{\rm{e}}^{ - j\omega }}) \end{array} \right. $$ (8) 式中,X(ejω)是x(n)的频率描述;Y(ejω)是y(n)的频率描述。由式(8)可得:
$$ Y({{\rm{e}}^{j\omega }}) = X({{\rm{e}}^{j\omega }})\cdot{\left| {H({{\rm{e}}^{j\omega }})} \right|^2} $$ (9) 由式(9)可知,输出Y(ejω)与输入X(ejω)之间不存在附加相位,FRR滤波实现了精确零相位失真。
Matlab软件中提供了一种零相位滤波的函数filtfilt,其本质上也是FRR滤波,但为了减小数字滤波都会遇到的边界效应问题,filtfilt函数在进行滤波前,在数据首尾两个方向上各进行了与滤波器节数相同的延拓[16]。设原数据为x(n),数据量大小为L2,滤波器节数为M,则拓展后的数据用Matlab语言可表达为:
$$ \begin{array}{l} \left[ {2\cdot x\left( 1 \right) - x\left( {M + 1: - 1:2} \right);x;} \right.\\ \left. {2\cdot x({L_2}) - x({L_2} - 1: - 1:{L_2} - M)} \right] \end{array} $$ 图 3为使用常规滤波方法和3种零相位滤波方法对某一固定周期的正弦信号进行滤波处理的结果,正弦信号的表达式为:x(t)=3sin(2π·2 000·t),采样频率为10 000 Hz,采样时间为0.025 s,滤波器设计为全通滤波器。从图 3可以看出,常规滤波与理论值存在一定的相位偏移,而FRR、RRF和filtfilt滤波与理论值相位保持一致,但FRR和RRF方法存在较明显的边界效应问题,在数据两端出现异常扰动。鉴于filtfilt函数的零相位滤波性能和对于边界效应的改善效果较好,本文将采用此方法进行零相位滤波。
2 Kaiser窗低通滤波器
本文需要一个低通滤波器滤除原始信号中的高频热噪声,提高信号的信噪比,以期改善定轨精度。常见的数字滤波器分为两种:无限脉冲响应(infinite impulse response,IIR)数字滤波器和有限脉冲响应(finite impulse response,FIR)数字滤波器。考虑到FIR数字滤波器能做到严格线性相位,且在滤波实现时可以采用快速傅里叶变换,因此本文使用窗函数法设计了一种低通FIR数字滤波器。
常见的窗函数有三角形窗、Hanning窗、Hamming窗、Blackman窗和Kaiser窗等。由于Kaiser窗具有可调性(即可选择不同参数α以适应各种不同需要),所以选用适应性较大的Kaiser窗作为低通滤波器的窗函数[17-19]。Kaiser窗的定义为:
$$ \omega \left( n \right) = \frac{{{I_0}\left( {\alpha \sqrt {1 - {{\left( {1 - \frac{{2n}}{{N - 1}}} \right)}^2}} } \right)}}{{{I_0}\left( \alpha \right)}} $$ (10) 式中,α为Kaiser窗参数;N为滤波器节数;n∈[0, N-1];I0(·)是第一类修正零阶贝塞尔函数。I0(x)可用下述级数来计算:
$$ {I_0}\left( x \right) = 1 + {\sum\limits_{k = 1}^\infty {\left[ {\frac{1}{{k!}}{{\left( {\frac{x}{2}} \right)}^k}} \right]} ^2} $$ (11) α是一个可调的参数,与主瓣宽度和旁瓣衰减有关。一般来说,α越大,过渡带越宽,阻带越小,衰减越大。图 4为不同α值对应的Kaiser窗的形状。若阻带最小衰减表示为As=-20lgδs,则α的确定可采用经验公式:
$$ \alpha = \left\{ \begin{array}{l} 0, {A_s} \le 21\\ 0.5842{\left( {{A_s} - 21} \right)^{0.4}} + 0.07886\left( {{A_s} - 21} \right), \\ 21 < {A_s} \le 50\\ 0.1102\left( {{A_s} - 8.7} \right), {A_s} > 50 \end{array} \right. $$ (12) 若滤波器通带和阻带波纹相等,即δp=δs时,则滤波器的节数N可以通过式(13)确定:
$$ N \approx \frac{{ - 20{\rm{lg}}{\delta _p} - 7.95}}{{\frac{{14.36\Delta \omega }}{{2{\rm{ \mathsf{ π} }}}}}} = \frac{{{A_s} - 7.95}}{{\frac{{14.36\Delta \omega }}{{2{\rm{ \mathsf{ π} }}}}}} $$ (13) Kaiser窗低通滤波器最基本的参数是截断频率ωc、滤波器节数N和Kaiser窗参数α。但在实际设计滤波器时,Kaiser窗参数α不够直观,通常选用滤波器衰减As,且As和α可通过式(12)转换。
3 数据处理与分析
3.1 仿真实验
为了验证零相位Kaiser窗滤波器滤除MEX多普勒数据噪声后对定轨的改善效果,本文首先进行仿真实验。使用MAGREAS对2010-06-26 T07:17:21至T12:06:09时段内的MEX双程多普勒测量过程进行模拟,得到精确的频率数据Data1(数据量为17 329,采样间隔为1 s)。在Data1上混入均值为0、方差为0.005 6 Hz的白噪数据(Data_noise),得到仿真数据Data2。将Data_noise通过零相位Kaiser窗低通滤波器(截断频率ωc=0.01 Hz,滤波器节数N=5 775,滤波器衰减As=-28.90 dB,滤波器最优参数通过反复测试得到)后叠加到Data1上得到Data3,将Data3作为仿真数据Data2滤除白噪后的结果。利用MAGREAS软件对模拟数据进行定轨,并对初轨在X、Y、Z方向各添加100 m的偏移。表 1为利用仿真数据(Data2和Data3)定轨后对初轨的修正结果。
仿真定轨数据 测速残差RMS
/(mm·s-1)定轨误差
/mΔX=-14.036 滤波前数据Data2 0.107 ΔY=-30.058 ΔZ=-25.156 ΔX=-4.594 滤波后数据Data3 0.038 ΔY=-10.406 ΔZ=-8.592 从表 1中可以看出,滤波后的数据比滤波前的数据所确定的轨道更为精确,双程多普勒测速残差RMS从0.107 mm/s减少至0.038 mm/s,定轨误差在X、Y、Z方向上也减小为原来的1/3左右。
3.2 MEX实测数据分析
通过仿真数据的测试可以初步确定,零相位Kaiser窗滤波器滤除多普勒数据中的噪声后,明显改善了定轨精度。为验证Kaiser窗滤波器在实际应用中的效果,选择MEX 2010-06-26 T07:17:21至T12:06:09弧段的多普勒实测数据进行测试。在该弧段内,数据按1 s间隔连续采样,数据量为17 329,数据文件为ESA 2010-06-26 MEX观测文件。滤波器的输入为残余频率,是观测数据和预测数据的差值; 滤波后的值叠加到预测数据上作为观测数据滤除白噪声后的结果。零相位Kaiser窗低通滤波器的参数设置为:截断频率ωc=0.002 5 Hz,滤波器节数N=5 774,滤波器衰减As=-100 dB,滤波器的最优参数通过反复测试得到。图 5为实测数据滤波前后的波形,是定轨之前的情况。图 6和图 7分别为利用原始数据和滤波后数据解算的轨道外推3 h后与ESA精密轨道在径向(R)、切向(T)和法向(N)3个方向上的差异。表 2为利用滤波前后数据解算的轨道对精密初轨的修正结果。其中,VX、VY、VZ分别表示航天器的初始速度在X、Y、Z方向上的分量;ΔVX、ΔVY、ΔVZ分别是对应的速度差。
MEX实测
定轨数据MAGREAS 测速残差RMS
/(mm·s-1)初轨修正结果 原始数据 0.107 ΔX=-62.761 m ΔY=-138.464 m ΔZ=-114.972 m ΔVX =-3.010 mm/s ΔVY= 12.698 mm/s ΔVZ=-44.553 mm/s 滤波后数据 0.031 ΔX=-3.260 m ΔY=-8.124 m ΔZ=-6.726 m ΔVX=-0.222 mm/s ΔVY= 0.859 mm/s ΔVZ=-2.145 mm/s 对比定轨结果发现,滤除实测数据中的白噪声后,定轨的改善效果明显。与仿真分析情况类似,双程多普勒测速残差RMS从0.107 mm/s减少至0.031 mm/s,本弧段的轨道位置和速度与ESA精密轨道的差异也明显变小。
由图 6和图 7可以看出,利用滤波后数据解算的轨道与ESA精密轨道的差异要明显小于原始数据的差异。图 6中原始数据对应的位置差异沿切向最大约为300 m,而图 7中滤波后数据对应的位置差异沿切向最大仅为50 m左右。此外,滤波后数据对应的速度差异也明显小于原始数据对应的速度差异。这一结果与表 2中给出的两种数据的测量精度相一致。
4 结语
通过仿真数据和实测数据的测试可以确定,零相位Kaiser窗低通滤波器滤除MEX多普勒频率数据中的白噪声后,对定轨精度有很大的改善。滤波后,双程测速数据残差RMS达到了0.031 mm/s,减小为原来的1/3左右;轨道位置和速度与ESA精密轨道的差异也明显变小。该滤波算法作为定轨前的数据预处理可以提高定轨的精度,从而为中国火星探测器的轨道数据处理提供一定的参考。
-
表 1 温州-丽水区域2015-2017年3种负荷引起的地壳垂直形变与重力变化统计
Table 1 Statistics of Crustal Vertical Deformation and Gravity Change Caused by Three Loads in Wenzhou-Lishui Region in 2015-2017
大地测量要素非潮汐变化 大气负荷影响 土壤水负荷影响 海平面变化影响 最大/最小值 均值/标准差 年变化幅度 最大/最小值 均值/标准差 年变化幅度 最大/最小值 均值/标准差 垂直形变/mm 6.86/-6.85 -0.00/4.33 12.4 2.72/-2.22 -0.00/0.92 3.5 1.26/-8.77 -2.23/1.38 重力变化/μGal 6.57/-6.64 0.00/3.82 12.1 4.24/-6.06 0.00/1.32 5.2 15.23/-6.16 -0.54/1.72 注:年变化幅度为每年中最大值与最小值之差的平均值 表 2 温州-丽水区域2015-2017年地下水负荷垂直形变与重力变化统计
Table 2 Statistics of Crustal Vertical Deformation and Gravity Change in Wenzhou-Lishui Region in 2015-2017
影响类型 最大值 最小值 年变化幅度 地壳垂直形变/mm 20.50 -15.43 28.60 地面重力变化/μGal 0.44 -0.59 0.80 表 3 温州-丽水区域2015-2017年地表环境负荷影响统计
Table 3 Statistics of the Impact of Total Surface Environmental Load in Wenzhou-Lishui Region in 2015-2017
影响类型 最大值 最小值 年变化幅度 地壳垂直形变/mm 13.12 -20.39 26.6 地面重力变化/μGal 9.00 -10.78 17.9 -
[1] 姜卫平, 李昭, 刘鸿飞, 等.中国区域IGS基准站坐标时间序列非线性变化的成因分析[J].地球物理学报, 2013, 56(7):2 228-2 237 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb201307010 Jiang Weiping, Li Zhao, Liu Hongfei, et al. Cause Analysis of the Non-Linear Variation of the IGS Reference Station Coordinate Time Series Inside China[J]. Chinese J Geophys, 2013, 56(7): 2 228-2 237 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb201307010
[2] 居恢扬, 顾仁杰.地质灾害研究的基础理论:环境因素控制论[J].自然灾害学报, 1994, 3(4):15-22 Ju Huiyang, Gu Renjie. Basic Theory on Study of Geological Calamitimes-Theory of Environmental Factors Control[J]. Journal of Natural Disasters, 1994, 3(4):15-22
[3] Tapley B D, Bettadpur S, Watkins M, et al. The Gravity Recovery and Climate Experiment: Mission Overview and Early Results[J]. Geophysical Research Letters, 2004, 31(9):4-5 doi: 10.1029/2004GL019920
[4] Dam V, Wahr J, Lavallée D. A Comparison of Annual Vertical Crustal Displacements from GPS and Gravity Recovery and Climate Experiment (GRACE) over Europe[J]. Journal of Geophysical Research: Solid Earth, 2007, 112(B3):200-203 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1029/2006JB004335
[5] Khan S A, Wahr J, Bevis M, et al. Spread of Ice Mass Loss into Northwest Greenland Observed by GRACE and GPS[J]. Geophysical Research Letters, 2010, 37(6):460-472 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1029/2010GL042460
[6] Bevis M, Alsdorf D, Kendrick E, et al. Seasonal Fluctations in the Mass of the Amazon River System and Earth's Elastic Response[J]. Geophys Res Lett, 2005, 32(16): L16308 doi: 10.1029/2005GL023491
[7] 王林松, 陈超, 邹蓉, 等.利用GPS与GRACE监测陆地水负荷导致的季节性水平形变:以喜马拉雅山地区为例[J].地球物理学报, 2014, 57(6):1 792-1 804 http://www.cnki.com.cn/Article/CJFDTotal-DQWX201406011.htm Wang Linsong, Chen Chao, Zou Rong, et al. Using GPS and GRACE to Detect Horizontal Deformation Caused by Loading Terrestrial Water: A Case Study in the Himalayas[J]. Chinese J Geophys, 2014, 57(6):1 792-1 804 http://www.cnki.com.cn/Article/CJFDTotal-DQWX201406011.htm
[8] Dong D, Fang P, Bock Y, et al. Anatomy of Apparent Seasonal Variations from GPS‐Derived Site Position Time Series[J]. Journal of Geophysical Research: Solid Earth, 2002, 107(B4):ETG-1-ETG9-16 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1029/2001JB000573
[9] Yan H M, Chen W, Zhu, Y Z, et al. Contributions of Thermal Expansion of Monuments and Nearby Bedrock to Observed GPS Height Changes[J]. Geophysical Research Letters, 2009, 36(13):88-97 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1029/2009GL038152
[10] 周东旭, 周兴华, 张化疑, 等.利用GPS连续观测进行中国沿海验潮站地壳垂直形变分析[J].武汉大学学报·信息科学版, 2016, 41(4):516-522 doi: 10.13203/j.whugis20140714 Zhou Dongxu, Zhou Xinghua, Zhang Huayi, et al. Analysis of the Vertical Deformation of China Coastal Tide Stations Using GPS Continuous Observations[J]. Geomatics and Information Science of Wuhan University, 2016, 41(4):516-522 doi: 10.13203/j.whugis20140714
[11] 李少魁.长江三角洲地区极端气候事件及其成因分析[D].南京: 南京信息工程大学, 2014 Li Shaokui. Analysis of Extreme Climate Events and Its Causes in Yangtz River Delta Region[D]. Nanjing: Nanjing Universiy of Information Science and Technology, 2014
[12] Wu X P, Heflin M B, Ivins E R, et al. Seasonal and Internnual Global Surface Mass Variations from Multisatellite Geodetic Data[J]. J Geophys Res, 2006, 111(B9): B09401 doi: 10.1029/2005JB004100/pdf
[13] Petit G, Luzum B. IERS Conventions 2010. International Earth Rotation and Reference Systems Service[R]. IERS Technical Note No. 36, Verlag des Bundesamtesfür Kartographie und Geodäsie, Frankfurt am Main, 2010
[14] 王伟, 章传银, 杨强, 等.大气负荷对区域地壳形变和重力变化的影响分析[J].武汉大学学报·信息科学版, 2018, 43(9): 1 302-1 308 doi: 10.13203/j.whugis20160392 Wang Wei, Zhang Chuanyin, Yang Qiang, et al. Impact of Atmospheric Pressure Loading on Regional Crustal Deformation and Gravity Change[J]. Geomatics and Information Science of Wuhan University, 2018, 43(9):1 302-1 308 doi: 10.13203/j.whugis20160392
[15] Rodell M, Houser P R, Jambor U, et al. The Global Land Data Assimilation System[R]. Bull Amer Meteor SOC, 2004, 85(3): 381-394
[16] Ablain M, Cazenave A, Valladeau G, et al. A New Assessment of the Error Budget of Global Mean Sea Level Rate Estimated by Satellite Altimetry over 1993-2008[J]. Ocean Science, 2009, 5(2):193-201 doi: 10.5194/os-5-193-2009
[17] Save H, CSR Level-2 Team. "GRACE RL06 Reprocessing and Results from CSR" EGU2018-10697[R]. EGU General Assembly 2012, Vienna, Austria, 2018
[18] 章传银, 王伟, 甘卫军, 等.利用CORS站网监测三峡地区环境负荷引起的地壳形变与重力场时空变化[J].武汉大学学报·信息科学版, 2018, 43(9): 1 287-1 294 doi: 10.13203/j.whugis20160419 Zhang Chuanyin, Wang Wei, Gan Weijun, et al. Monitoring Temporal and Spatial Changes of Crustal Deformation and Gravity Field Caused by Environmental Load in the Three Gorges Reservoir Region Based on CORS[J]. Geomatics and Information Science of Wuhan University, 2018, 43(9):1 287- 1 294 doi: 10.13203/j.whugis20160419
[19] Argus D F, Fu Y, Landerer F W. Seasonal Variation in Total Water Storage in California Inferred from GPS Observations of Vertical Land Motion[J]. Geophysical Research Letters, 2014, 41(6); 1 971-1 980 doi: 10.1002/2014GL059570
[20] Han S C, Shum C K, Jekeli C, et al. Non-Isotropic Filtering of GRACE Temporal Gravity for Geophysical Signal Enhancement[J]. Geophysical Journal International, 2005, 163(1): 18-25 doi: 10.1111/j.1365-246X.2005.02756.x
[21] Zhang Z, Chao B F, Yang L, et al. An Effective Filtering for GRACE Time-Variable Gravity: Fan Filter[J]. Geophysical Research Letters, 2009, 36(17): 1 397-1 413 http://www.onacademic.com/detail/journal_1000035772156710_8210.html
[22] Cheng M, Tapley B D. Variations in the Earth's Oblateness During the Past 28 Years[J]. Journal of Geophysical Research Solid Earth, 2005, 110(B9): 1 404-1 406 http://gji.oxfordjournals.org/external-ref?access_num=10.1029/2004JB003028&link_type=DOI
[23] 王伟, 党亚民, 章传银, 等.基于CORS站网监测三峡地区陆地水负荷对地表形变和重力变化的影响[J].地球物理学报, 2017, 60(3):962-971 http://www.cnki.com.cn/Article/CJFDTotal-DQWX201703012.htm Wang Wei, Dang Yaming, Zhang Chuanyin, et al. Monitoring of Crustal Deformation and Gravity Variation from Terrestrial Water Loading in the Three Gorges Region by the CORS Network[J]. Chinese J Geophys, 2017, 60(3):962-971 http://www.cnki.com.cn/Article/CJFDTotal-DQWX201703012.htm
-
期刊类型引用(1)
1. 郭斐,邱耀东,王正涛. 基于FIR滤波的磁静期电离层-磁层磁场信号分离. 武汉大学学报(信息科学版). 2019(06): 844-850 . 百度学术 其他类型引用(1)