文章快速检索     高级检索
  大地测量与地球动力学  2019, Vol. 39 Issue (2): 189-194  DOI: 10.14075/j.jgg.2019.02.015

引用本文  

王海涛, 王斌, 王伟, 等. 大气负荷对新疆地区地壳形变和地面重力变化的影响[J]. 大地测量与地球动力学, 2019, 39(2): 189-194.
WANG Haitao, WANG Bin, WANG Wei, et al. Research on the Effect of Atmospheric Pressure Load on Crust Deformation and Gravity Change in Xinjiang Region[J]. Journal of Geodesy and Geodynamics, 2019, 39(2): 189-194.

项目来源

国家自然科学基金(41304009,41374081);国家重点研发计划(2016YFB0501702)。

Foundation support

National Natural Science Foundation of China, No. 41304009, 41374081; National Key Research and Development Program of China, No. 2016YFB0501702.

第一作者简介

王海涛,硕士生,主要研究方向为大地测量数据处理,E-mail:1280181825@qq.com

About the first author

WANG Haitao, postgraduate, majors in geodetic data processing, E-mail: 1280181825@qq.com.

文章历史

收稿日期:2018-03-07
大气负荷对新疆地区地壳形变和地面重力变化的影响
王海涛1,2     王斌1     王伟2     章传银2     
1. 国家测绘地理信息局大地测量数据处理中心,西安市友谊东路334号,710054;
2. 中国测绘科学研究院,北京市莲花池西路28号,100830
摘要:基于ECMWF全球大气模型和中国地面气候资料日值数据集中新疆及周边的气象站气压数据,根据弹性地壳形变理论,采用移去-恢复法,利用负荷格林函数法和球谐函数法,计算新疆及周边区域2011~2015年大气负荷变化对地壳垂直、水平形变的影响,同时采用大气导纳方法计算大气变化对地面重力的影响。结果表明,大气负荷对新疆地区垂直形变的影响达到cm级,对地面重力的影响可达10 μGal;大气负荷影响具有明显的季节性。
关键词大气负荷球谐函数移去-恢复法地壳形变地面重力

现代大地测量和地球物理观测技术得到不断突破,超导重力仪的观测精度可达±0.1 μGal,甚长基线干涉测量(VLBI)、卫星激光测距(SLR)、全球定位系统(GPS)基线解算精度可达±0.1~±1 cm(10-8~10-9)量级[1]。要获得高精度的大地测量结果,必须进行严密的数据处理和严格的干扰信号剔除,而干扰信号中大气变化是主要干扰噪声。大气压周期性变化引起地面点地壳垂直形变可达±25 mm,对地面重力影响可达±20 μGal[2]。显然,对于高精度的大地测量,如参考框架的建立、基准点的确定,大气影响改正不可忽略,因此精确的大气负荷变化研究不仅可用于GPS基线解算和精密重力测量等其中的大气改正,而且对提高不连续观测站的观测精度也具有重要的实用价值[3]

大气负荷影响的计算方法主要有负荷格林积分法和球谐函数法。负荷格林积分法是Farrell[4-5]基于Longman地球负荷理论,从理论上推出球对称地球在单位点质量负荷下的静态格林函数,解决其计算问题,即结合气压变化数据和格林函数的褶积积分,可以得到大气变化产生的位移、重力、倾向和应变效应。罗少聪等[1]利用弹性地球负荷理论和流体静力学平衡,计算大气压作用下地表位移重力和格林倾斜函数,表明大气效应弹性影响主要来自于远离台站的区域,引力项对重力影响的主要贡献来自于台站近区。许多学者[6-8]又根据Farrell理论给出应用球谐函数直接计算地表流体负荷(大气、海洋、陆地水、冰雪圈等)质量变化产生的负荷响应公式。张诗玉等[3]利用球谐系数方法研究大气负荷效应对GPS基准站的影响表明,大气负荷引起的地壳垂直形变可达20~30 mm,表现为明显的季节性变化。沈迎春等[9]比较上述两种计算方法的结果精度和运算速度表明,两种方法在数学本质上是一致的,计算精度和运算速度在相同分辨率下差异不大。

本文利用全球大气变化模型资料,结合研究区域内气象站实测气压变化数据,采用移去-恢复法计算新疆区域大气负荷变化对地壳形变和重力的影响,并分析其影响特征。

1 理论方法 1.1 基本理论和公式

根据气压原理将全球大气模型转换为等效水高变化,地面一点的等效水高变化量Δhw($φ$, λ)通过式(1)规格化展开为球谐系数[10-11],进而计算相应的负荷形变:

$ \begin{array}{l} \mathit{\Delta }{h_{\rm{w}}}(\varphi, \lambda ) = R\sum\limits_{n = 1}^N {\sum\limits_{m = 0}^n {\left[{\mathit{\Delta }C_{nm}^q{\rm{cos}}m\lambda + } \right.} } \\ \left. {\mathit{\Delta }S_{nm}^q{\rm{sin}}m\lambda } \right]{{\bar P}_{nm}}({\rm{sin}}\varphi ) \end{array} $ (1)

根据地球负荷形变和重力场理论,地面等效水高变化引起一点的大地高(垂直位移形变)负荷变化量为:

$ \begin{array}{l} \;\;\mathit{\Delta }r(\varphi, \lambda ) =-3\frac{{\rho {}_{\rm{w}}}}{{\rho {}_{\rm{e}}}}\frac{{GM}}{{R2}}\sum\limits_{n = 1}^N {\frac{{{{h'}_n}}}{{2n + 1}} \cdot } \\ \sum\limits_{m = 0}^n {(\mathit{\Delta }C_{nm}^q{\rm{cos}}m\lambda + \mathit{\Delta }S_{nm}^q{\rm{sin}}m\lambda )} {{\bar P}_{nm}}({\rm{sin}}\varphi ) \end{array} $ (2)

引起的地面重力变化为:

$ \begin{array}{l} {g_t}(\varphi, \lambda ) =-3\frac{{\rho {}_{\rm{w}}}}{{\rho {}_{\rm{e}}}}\frac{{GM}}{{R2}}\sum\limits_{n = 1}^N {\frac{{n + 2{{h'}_n}-(n + 1){{k'}_n}}}{{2n + 1}} \cdot } \\ \;\;\;\;\;\;\;\sum\limits_{m = 0}^n {(\mathit{\Delta }C_{nm}^{\rm{q}}{\rm{cos}}m\lambda + \mathit{\Delta }S_{nm}^q{\rm{sin}}m\lambda )} {{\bar P}_{nm}}({\rm{sin}}\varphi ) \end{array} $ (3)

式(1)~(3)中,$φ$λ为地面点的地心纬度和经度,R为地球平均半径,ΔCnmqΔSnmqnm次规格化球谐系数,Pnm(·)为nm次规格化勒让德函数,h′为n阶径向负荷勒夫数,knn阶负荷勒夫数,ρw≈1.0×103 kg/m3为水的密度,ρe≈5.5×103 kg/m3为固体地球的平均密度,G为万有引力常数,M为地球的总质量。

同样,大气负荷引起的地面重力变化和地表大地高变化也可以用格林函数积分法进行计算。负荷格林函数积分通用公式为:

$ \mathit{\Delta \Theta }\left( {\varphi, \lambda } \right) = G{\rho _{\rm{w}}}\mathop \int\!\!\!\int \limits_S^{} \frac{{{h_{\rm{w}}}\left( {\varphi ', \lambda '} \right)}}{L}G\left( \psi \right){\rm{d}}S $ (4)

式中,($φ$′, λ′)为地面积分流动点,($φ$, λ)为地面计算点,L为流动点到计算点的空间距离,ψ是球面角距,S为地球表面,G(ψ)为通用的负荷格林函数。

重力格林函数计算公式为:

$ \begin{array}{l} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{G^g}\left( {\varphi, \lambda } \right){\rm{ = }}\\ - \frac{\gamma }{M}\sum\limits_{n = 1}^\infty {\left[{n-(n + 1){{k'}_n} + 2{{h'}_n}} \right]{P_n}{\rm{cos}}(\psi )} \end{array} $ (5)

式中,Pn(·)为n阶勒让德函数。

大地高变化格林函数Gr为:

$ {G^r}\left( {\varphi, \lambda } \right){\rm{ = }}\frac{R}{M}\sum\limits_{n = 1}^\infty {{{h'}_n}{P_n}{\rm{cos}}(\psi )} $ (6)

将式(6)代入式(4),得到计算任意点($φ$, λ)大地高变化公式:

$ \begin{array}{l} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\mathit{\Delta \Theta }\left( {\varphi, \lambda } \right) = \\ G{\rho _{\rm{w}}}\mathop \int\!\!\!\int \limits_S^{} \frac{{{h_{\rm{w}}}\left( {\varphi ', \lambda '} \right)}}{L}\frac{R}{M}\sum\limits_{n = 1}^\infty {{{h'}_n}{P_n}{\rm{cos}}(\psi )} {\rm{d}}S \end{array} $ (7)
1.2 移去-恢复法理论

移去-恢复方法广泛应用于拟合大地水准面和重力场研究[12]。采用移去-恢复法,可以发挥区域实测气压变化优势,显现区域差异,提高模型的区域适用性,准确表现区域负荷变化特征。在移去过程中,需要将式(7)进行改进,将式中的积分域S分成近区S0和远区S-S0,则近区的等效水高积分公式变为:

$ \begin{array}{l} \mathit{\Delta \Theta }(\varphi, \lambda ){|_{{S_0}}} = G{\rho _{\rm{w}}}\mathop \int\!\!\!\int \limits_S^{} \frac{{{h_{\rm{w}}}\left( {\varphi ', \lambda '} \right)}}{L}\frac{R}{M} \cdot \\ \;\;\;\sum\limits_{n = 1}^\infty {{{h'}_n}{P_n}{\rm{cos}}(\psi )} {\rm{d}}S + \mathop \int\!\!\!\int \limits_{{S_0}}^{} {\rm{d}}{h_{\rm{w}}}G({\psi _0}){\rm{d}}S \end{array} $ (8)

式中,dhw=Δhw($φ$′, λ′)-ΔhwM($φ$′, λ′),称之为残差等效水高;Δhw($φ$′, λ′)为区域实测等效水高;ΔhwM($φ$′, λ′)为全球大气模型对应的等效水高,称为球谐系数等效水高。

移去-恢复法主要流程为:

1) 将全球大气模型数据扣除基准值,转换为等效水高值,得到每月全球大气等效水高相对于基准值的变化量,根据球谐分析法展开为等效水高球谐系数。

2) 利用等效水高球谐系数计算气象站位置处的模型等效水高,将气象站气压值同样扣除基准值后转换为等效水高。移去球谐系数等效水高,得到气象站位置处模型与实测的差异,称之为残差等效水高。这一步称之为“移去”。

3) 将站点残差等效水高采用克里金插值法格网化[2],基于格林函数积分法计算残差等效水高对区域地壳形变和地面重力变化的影响。

4) 将全球大气等效水高球谐系数计算的地壳形变或地面重力变化与残差等效水高计算的地壳形变或地面重力变化相加,得到大气负荷的间接影响。这一步称之“恢复”。

5) 由于大气对地面重力变化的影响还包括负荷的直接影响,即大气负荷对重力变化的影响为间接影响与直接影响之和。

2 研究区域和数据处理 2.1 研究区域

本文研究区及区域内气象站分布如图 1所示。

图 1 研究区域和气象站分布 Fig. 1 Distribution of meteorological stations in the studied area
2.2 数据处理 2.2.1 全球大气模型处理

将下载自欧洲中期天气预报中心ECMWF(European centre for medium-range weather forecasts)的0.25°×0.25°月平均值全球大气模型数据扣除一个基准值,作为气压变化量[13]

根据气压原理,按照1 hPa气压值等于10.197 2 mm等效水高的关系,将全球大气数据资料转换为相应的等效水高。根据式(1)将其展开为360阶球谐系数,并根据区域内已知的气象站点位置坐标求得每月等效水高,作为全球大气模型等效水高数据,在后续将进行“移除”操作。利用区域内0.25°×0.25°分辨率的坐标数据,根据式(2)求得区域内各点的垂直形变量。

2.2.2 气象站观测气压数据处理

气象站气压观测资料来自于国家气象科学数据共享服务平台网提供的中国地面气候资料日值数据集。通过将气象站观测数据剔除缺省值和粗差后绘图发现,该区域内低海拔气象站气压年变化表现为:一般区域夏季气压低、冬季气压高,如图 2中51243站海拔449.5 m;而位于昆仑山脉高海拔地区的气象站气压变化呈现夏季气压高、冬季气压低的反转现象,且年变化量较小,如图 2中52908站海拔4 612.2 m,呈现区域变化差异的复杂性。之后将气象站实测日值数据转换为月平均值,并扣除同一月份的基准值,根据气压原理将气压值转换为等效水高,作为区域实测气压值变化。减去全球大气等效水高球谐系数计算的等效水高,即移去模型等效水高值,得到的值称为残差等效水高,并采用surfer中的克里金插值法将残差等效水高进行区域格网化[2]。根据等效水高格林积分法(式(7))计算垂直形变的影响,高程数据采用等分辨率的SRTM数据。

图 2 气象站气压日值变化 Fig. 2 Variation of daily data of atmospheric pressure at the meteorological stations
3 形变结果分析 3.1 地壳形变

本文研究时间跨度为2011~2015年,在基准选择上根据研究时间跨度,选择2011-05的月平均值作为基准值扣除。将球谐系数计算的全球大气模型地壳形变与格林函数计算的残差等效水高相应的地壳形变相加,得到最终大气负荷对新疆区域地壳形变的总影响。

根据计算的大气负荷对地壳形变的影响,表现为明显的年周期变化。图 3图 4为2011年和2013年中4个月的地壳垂直形变,由图可见,垂直形变最大年变化量出现在吐鲁番盆地,大约为25 mm,这与该地区气压年变化量较大相符。图 5图 6为2012年和2014年中4个月的地壳水平形变变化情况,由图可见,水平形变量较小,在1 mm左右。

图 3 2011年大气负荷对地壳垂直形变的影响 Fig. 3 Effect of atmospheric load on vertical deformation of the crust in 2011

图 4 2013年大气负荷对地壳垂直形变的影响 Fig. 4 Effect of atmospheric load on vertical deformation of the crust in 2013

图 5 2012年大气负荷对地壳水平形变的影响 Fig. 5 Effect of atmospheric loading on horizontal deformation of the crust in 2012

图 6 2014年大气负荷对地壳水平形变的影响 Fig. 6 Effect of atmospheric load on horizontal deformation of the crust in 2014
3.2 地面重力影响

大气负荷对地面重力的影响分为直接影响和间接影响,直接影响就是气压变化直接引起的重力变化;间接影响是指大气负荷变化导致的固体地球内部密度质量重新分布引起的重力变化[9],直接影响占主要部分。

采用大气导纳方法计算大气负荷对地面重力变化的影响,理论大气导纳值为-0.303 7 μGal/hPa[2],包含直接影响和间接影响。由全球大气模型中提取区域大气数据,扣除参考基准值,得到大气变化数据,然后乘以理论大气导纳值,得到区域大气负荷对地面重力的影响。2015-03、2015-06、2015-09、2015-12大气对地面重力的影响变化如图 7所示,由图可见,重力变化与气压值变化成反比,年变化幅度在10 μGal左右。

图 7 2015年大气负荷对地面重力的影响 Fig. 7 Effect of atmospheric load on ground gravity in 2015
3.3 大气负荷形变统计

为了全面了解大气负荷对新疆区域大地测量要素的影响,对其进行统计分析,如表 1所示。大气负荷对地壳形变和地面重力变化的影响具有年周期变化特征,通过统计其年变化幅度,了解其影响量级和变化规律。由表 1看出,大气负荷对垂直形变的影响最大为12.00 mm、最小为-17.97 mm,对水平影响较小(在1 mm以内);大气对地面重力的影响与大气变化具有明显的负相关性,年变化幅度在11.00 μGal以内。

表 1 大气负荷对大地测量要素影响统计 Tab. 1 Influence statistics of atmospheric load on geodetic factors
4 结语

本文采用移去-恢复法,将全球大气模型资料和气象站实测气压数据结合起来,利用球谐函数法和格林函数法,计算新疆区域的大气负荷月变化对区域地壳形变和地面重力的影响。结果表明,大气负荷影响具有明显的年周期变化特征,在新疆地区大气负荷对地壳垂直形变影响的年变化量超过20 mm,对水平形变影响的年变化量在1.2 mm以内;对地面重力的影响超过10 μGal,与区域气压变化量呈现正相关。

参考文献
[1]
罗少聪, 孙和平, 徐建桥. 大气负荷引起的重力与位移变化理论计算的精度估计[J]. 测绘学报, 2001, 30(4): 309-315 (Luo Shaocong, Sun Heping, Xu Jianqiao. Theoretical Computation in the Precision Estimation of the Change in Gravity and Displacements Due to Atmospheric Loading[J]. Acta Geodaetica et Cartographica Sinica, 2001, 30(4): 309-315 DOI:10.3321/j.issn:1001-1595.2001.04.007) (0)
[2]
罗少聪.大气负荷效应问题研究[D].武汉: 中国科学院测量与地球物理研究所, 2003 (Luo Shaocong. Study on the Loading Effects on the Atmospheric Pressure[D]. Wuhan: Institute of Geodesy and Geophysics, CAS, 2003) (0)
[3]
张诗玉, 钟敏, 唐诗华. 我国GPS基准站地壳垂直形变的大气负荷效应[J]. 武汉大学学报:信息科学版, 2006, 31(12): 1090-1093 (Zhang Shiyu, Zhong Min, Tang Shihua. Vertical Crustal Displacements Due to Atmospheric Loading Effects at GPS Fiducial Stations in China[J]. Geomatics and Information Science of Wuhan University, 2006, 31(12): 1090-1093) (0)
[4]
Farrell W E. Deformation of the Earth by Surface Loads[J]. Reviews of Geophysics, 1972, 10(3): 761-797 (0)
[5]
Longman I M. A Green's Function for Determining the Deformation of the Earth under Surface Mass Loads[J]. Journal of Geophysical Research, 1962, 67(2): 845-850 (0)
[6]
Spratt R S. Modeling the Effect of Atmospheric Pressure Variations on Gravity[J]. Geophysical Journal International, 2010, 71(1): 173-186 (0)
[7]
Mangiarotti S, Cazenave A, Soudarin L, et al. Annual Vertical Crustal Motions Predicted from Surface Mass Redistribution and Observed by Space Geodesy[J]. Journal of Geophysical Research:Solid Earth, 2001, 106(B3): 4277-4292 (0)
[8]
Mitrovica J X, Davis J L, Shapiro I I. A Spectral Formalism for Computing Three-Dimensional Deformations Due to Surface Loads: 2. Present-Day Glacial Isostatic Adjustment[J]. Journal of Geophysical Research :Solid Earth, 1994, 99(B4): 7075-7101 (0)
[9]
沈迎春, 闫昊明, 彭鹏, 等. 质量负荷引起地表形变的格林函数和球谐函数方法对比研究[J]. 武汉大学学报:信息科学版, 2017, 42(7): 1008-1014 (Shen Yingchun, Yan Haoming, Peng Peng, et al. Comparative Study of Green's Function and Spherical Harmonic Function Methods on Surface Deformation Caused by Mass Loading[J]. Geomatics and Information Science of Wuhan University, 2017, 42(7): 1008-1014) (0)
[10]
王伟, 党亚民, 章传银, 等. 基于CORS站网监测三峡地区陆地水负荷对地壳形变和重力变化的影响[J]. 地球物理学报, 2017, 60(3): 962-971 (Wang Wei, Dang Yamin, 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 Journal of Geophysics, 2017, 60(3): 962-971) (0)
[11]
党亚民. 现代大地测量基准[M]. 北京: 测绘出版社, 2015 (Dang Yamin. Modern Geodetic Datum[M]. Beijing: Surveying and Mapping Press, 2015) (0)
[12]
郭春喜, 宁津生, 陈俊勇, 等. 珠峰地区似大地水准面精化与珠峰顶正高的确定[J]. 地球物理学报, 2008, 51(1): 101-107 (Guo Chunxi, Ning Jinsheng, Chen Junyong, et al. Improvement of Regional Quasi-Geoid in Qomolangma(Mt. Everest) and Determination of Orthometric Elevation[J]. Chinese Journal of Geophysics, 2008, 51(1): 101-107 DOI:10.3321/j.issn:0001-5733.2008.01.013) (0)
[13]
Ray R D, Ponte R M. Barometric Tides from ECMWF Operational Analyses[J]. Ann Geophys, 2003, 21(8): 1897-1910 (0)
Research on the Effect of Atmospheric Pressure Load on Crust Deformation and Gravity Change in Xinjiang Region
WANG Haitao1,2     WANG Bin1     WANG Wei2     ZHANG Chuanyin2     
1. Geodetic Data Processing Centre, NASMG, 334 East-Youyi Road, Xi'an 710054, China;
2. Chinese Academy of Surveying & Mapping, 28 West-Lianhuachi Road, Beijing 100830, China
Abstract: Based on the ECMWF global atmospheric model and the daily data of the surface climatological data in China, the pressure data of the meteorological stations in Xinjiang and surrounding areas are concentrated. According to the theory of elastic crustal deformation, the method of removal and recovery is adopted, and Green's function and the spherical harmonic function method are used to study and calculate its effect on atmospheric load on the vertical and horizontal deformation of the crust, during the period from 2011 to 2015 in Xinjiang and the surrounding area. At the same time, the atmospheric admittance method is used to calculate the influence of atmospheric variance on ground gravity. The results indicate: 1) the influence of atmospheric load on vertical deformation in Xinjiang reaches the centimeter level, and the influence on ground gravity can reach 10 μGal; 2) the influence of atmospheric load has significant seasonality.
Key words: atmospheric pressure load; spherical harmonic function; remove-restore method; crustal deformation; gravity over the earth surface