文章快速检索    
  地震地磁观测与研究  2021, Vol. 42 Issue (2): 214-216  DOI: 10.3969/j.issn.1003-3246.2021.02.069
0

引用本文  

翟丽娜, 王松阳, 孔祥瑞. 海潮负荷对大连地震台重力固体潮的影响[J]. 地震地磁观测与研究, 2021, 42(2): 214-216. DOI: 10.3969/j.issn.1003-3246.2021.02.069.
ZHAI Lina, WANG Songyang, KONG Xiangrui. The influence of ocean tide loading on the gravity earth tide of Dalian Seismic Station[J]. Seismological and Geomagnetic Observation and Research, 2021, 42(2): 214-216. DOI: 10.3969/j.issn.1003-3246.2021.02.069.

基金项目

震情跟踪青年项目(项目编号: 2021010204);星火计划青年项目(项目编号: XH21013Y); 星火计划攻关项目(项目编号: XH19059);辽宁省地震局青年基金(项目编号: LZ201911)

作者简介

翟丽娜(1986-), 女, 工程师, 主要从事地壳形变与地震预测方面的研究工作。E-mail: 529436296@qq.com
海潮负荷对大连地震台重力固体潮的影响
翟丽娜 , 王松阳 , 孔祥瑞     
中国沈阳 110034 辽宁省地震局
关键词gPhone重力仪    固体潮    海潮负荷    海潮模型    调和分析    
The influence of ocean tide loading on the gravity earth tide of Dalian Seismic Station
ZHAI Lina , WANG Songyang , KONG Xiangrui     
Liaoning Earthquake Agency, Shenyang 110034, China
Key words: gPhone gravimeter    solid tide    ocean tide loading    regional OTL models    harmonic analysis    

海潮是海洋表面受到太阳和月亮引力后产生的周期性涨落现象,与重力固体潮具有同源性,会产生相同的频谱特征,因此,在重力固体潮观测资料中会夹杂海潮信息(张邵栋等,1989)。分析重力固体潮汐残差信号可以获得地球内部动力学信息,但其中夹杂的全球和局部海潮负荷对重力固体潮影响较大,仅采用滤波法难以剔除,尤其在沿海地区,海潮负荷在重力潮汐振幅中占比更大(杜文成等,2018)。为了更加有效地利用已知重力资料精密测定地球潮汐常数等地球物理学参数,研究地核运动等地球动力学问题,有必要对重力观测资料进行海潮负荷信号改正(孙和平等,2002)。

Farrell(1972)首次提出采用海潮潮高和负荷格林函数进行褶积积分,求取海潮负荷影响;Schwiderski(1980)采用流体动力学差值方法构建较高精度的SCW80海潮模型,被广泛认可并作为标准使用。伴随着现代高精度Topex/Poseidon(T/P)卫星测高技术的发展,结合地域特点,学者们提出多种高精度海潮模型(Sun,1992Neumeyer et al,2005Bos et al,2005周江存等, 2005, 2007, 2008孙和平等, 2005, 2006Florsch et al,2009Arnoso et al,2011Kim et al,2011Kim and Shibuya, 2013Zhou et al,2013;徐建桥等,2014;Ducarme et al,2014;刘清超等,2016;杨锦玲等,2016;杜文成等,2018翟丽娜等,2020)。辽宁省大连市国家基准台(下文简称大连台)是距海岸线较近的重力观测台站之一,海潮对重力固体潮观测资料的影响较大。文中搜集了现阶段应用较广泛的8个海潮模型:FES2004、GOT4.7、EOT11a、DTU10、NAO99b、HAMTIDE11a、TPXO72atlas和TPXO72,对大连基准地震台重力固体潮受到的海潮负荷影响进行分析。

根据Farrell等(1972)提出的理论,结合Agnew(1997)的积分格林函数计算法,将海潮与重力格林函数做褶积积分计算(刘清超等,2016;杨锦玲等,2016;杜文成等,2018),具体计算公式为

$ \begin{aligned} L(\theta, \lambda, t) &=a^{2} \iint \rho H\left(\theta^{\prime}, \lambda^{\prime}, t^{\prime}\right) \cdot G(\psi) S(A) \sin \theta^{\prime} \mathrm{d} \theta^{\prime} \mathrm{d} \lambda^{\prime} \\ &=a^{2} \iint \rho H(\psi, A, t) G(\psi) S(A) \sin \psi \mathrm{d} \psi \mathrm{d} A \end{aligned} $ (1)

其中,L(θ, λ, t)表示某时刻某点的海潮负荷量,θλ分别为计算点坐标,t为时间;ρ为海水密度;a为地球半径;H(θ′, λ′, t′)为负荷点(θ′, λ′)处的瞬时潮高;G(ψ)表示与角距相关的重力负荷格林函数(该函数与负荷勒夫数相关)。

gPhone型连续重力观测仪置于观测山洞,为了更加真实反映海潮影响程度,剔除重力干扰数据,估算重力固体潮各个波群的潮汐参数及精度,选取适合大连地区的8个海潮模型(FES2004、GOT4.7、EOT11a、DTU10、NAO99b、HAMTIDE11a、TPXO72atlas和TPXO72),对主要潮波S2、M2、K2、N2、K1、O1、Q1和P1进行海潮负荷改正计算,对负荷改正效果进行综合评价,发现3个主潮波的残差负荷改正有效性为43%—53%,其中对于K1波和O1波,TPXO7.2atlas模型海潮改正效果最好,对于M2波,TPXO7.2模型潮波改正效果最好(图 1)。

图 1 K1、O1和M2波海潮负荷改正有效性 Fig.1 Effectiveness of K1, O1 and M2 wave sea tide loading correction

总体来说,采用不同模型,主潮波的剩余残差值仍较大,可能与大连海湾位置的特殊性有关,近海模型存在误差,若采用近海高精度海潮模型进行修正,并在预处理方法中尝试改进,可能有助于计算效果的提升。

参考文献
杜文成, 袁林果, 张宁宁, 等. 青岛台站重力固体潮和海潮负荷特征研究[J]. 大地测量与地球动力学, 2018, 38(9): 913-916.
孙和平, Ducarme B, 许厚泽, 等. 基于全球超导重力仪观测研究海潮和固体潮模型的适定性[J]. 中国科学(D辑), 2005a, 35(7): 649-657.
孙和平, 许厚泽, 陈武, 等. 香港地区重力固体潮和海潮负荷特征研究[J]. 地球物理学报, 2006, 49(3): 724-734. DOI:10.3321/j.issn:0001-5733.2006.03.016
孙和平, 许厚泽, 周江存, 等. 武汉超导重力仪观测最新结果和海潮模型研究[J]. 地球物理学报, 2005b, 48(2): 299-307.
孙和平, 周江存. 中国地壳运动观测网络基准站重力场变化的海潮负荷信号改正问题[J]. 地球科学进展, 2002, 17(1): 39-43. DOI:10.3321/j.issn:1001-8166.2002.01.007
翟丽娜, 孔祥瑞, 杨锦玲, 等. 大连台重力固体潮受海潮负荷影响特征研究[J]. 中国地震, 2020, 36(4): 953-960.
周江存, 孙和平. 海潮对卫星重力场恢复的影响[J]. 地球物理学报, 2007, 50(1): 115-121. DOI:10.3321/j.issn:0001-5733.2007.01.016
周江存, 孙和平. 用东海和南海潮汐资料修正全球海潮模型对中国及邻区重力场负荷计算的影响[J]. 地震学报, 2005, 27(3): 332-338. DOI:10.3321/j.issn:0253-3782.2005.03.012
周江存. 固体潮和卫星重力场恢复中的负荷问题及其应用研究[D]. 武汉: 中国科学院测量与地球物理研究所, 2008: 39-43.
Arnoso J, Benavent M, Bos M S, et al. Verifying the body tide at the Canary Islands using tidal gravimeter observations[J]. J Geodyn, 2011, 51(5): 358-365. DOI:10.1016/j.jog.2010.10.004
Bos M S, Baker T F. An estimate of the errors in gravity ocean tide loading computations[J]. J Geod, 2005, 79(1/2/3): 50-63. DOI:10.1007/s00190-005-0442-5
Farrell D E, Huebener R P, Kampwirth R T. Direct observation of growth of the landau domain structure[J]. Solid State Communication, 1972, 12: 1647-1649.
Florsch N, Llubes M, Wöppelmann G, et al. Oceanic loading monitored by ground-based tiltmeters at Cherbourg (France)[J]. J Geodyn, 2009, 48(3/4/5): 211-218.
Kim T H, Shibuya K, Doi K, et al. Validation of global ocean tide models using the superconducting gravimeter data at Syowa Station, Antarctica, and in situ tide gauge and bottom-pressure observations[J]. Polar Sci, 2011, 5(1): 21-39. DOI:10.1016/j.polar.2010.11.001
Kim T H, Shibuya K. Verification of the ellipsoidal earth model with an inelastic and convective mantle using tidal gravity variations revisited[J]. Geophys J Int, 2013, 194(1): 230-248. DOI:10.1093/gji/ggt108
Neumeyer J, del Pino J, Dierks O, et al. Improvement of ocean loading correction on gravity data with additional tide gauge measurements[J]. J Geodyn, 2005, 40(1): 104-111. DOI:10.1016/j.jog.2005.07.012
Schwiderski E W. On charting global ocean tides[J]. Rev Geophys, 1980, 18(1): 243-268. DOI:10.1029/RG018i001p00243
Sun H P. Comprehensive researches for the effect of the ocean loading on gravity observations in the western Parific area. Bull Int[J]. Marees Terrestres, 1992, 113: 8271-8292.
Zhou J C, Hwang C W, Sun H P, et al. Precise determination of ocean tide loading gravity effect for absolute gravity stations in coastal area of China: effects of land-sea boundary and station coordinate[J]. J Geodyn, 2013, 68: 29-36. DOI:10.1016/j.jog.2013.03.003