刘清超:南极中山站的重力和区域海洋潮汐特征

南极中山站的重力和区域海洋潮汐特征

刘清超^1,2 , 孙和平^1,2 , 王泽民³ , 陈晓东¹ , 张胜凯³ , 徐建桥¹     

1. 中国科学院测量与地球物理研究所大地测量与地球动力学国家重点实验室, 湖北 武汉 430077; 
2. 中国科学院大学, 北京 100049; 
3. 武汉大学中国南极测绘研究中心, 湖北 武汉 430079

之一作者简介：刘清超(1990-), 男, 博士生, 研究方向为重力场理论、 *** 及应用。E-mail:liuqch1990@163.com

微信号：MeetyXiao
添加微信好友, 获取更多信息
复制微信号

摘要：本文利用中山站弹簧重力仪记录的重力潮汐时间序列、验潮站数据、CATS2008区域和Eot11a全球海潮模型研究重力和海洋潮汐特征。结果表明，在周日频段，潮波O1的海潮振幅达到28 cm，4个主要潮波(Q1、O1、P1和K1)的全球模型与验潮站潮高差之和为4.2 cm，区域模型与验潮站潮高差之和为4.4 cm; 在半日频段，潮波M2的海潮振幅达到20 cm，4个主要潮波(N2、M2、S2和K2)的潮高差之和分别为7.7 cm和5.1 cm，说明利用区域模型修正全球模型的重要性。经区域模型修正的全球海潮负荷改正后，重力主波K1、M2和S2的最终残差振幅分别下降了9.84%、56.14%和37.08%，说明区域海潮模型更能反映海洋潮汐的真实特征，用区域模型修正全球海潮模型的有效性得到验证。

关键词：南极    重力潮汐    区域和全球海潮模型    

The Characteristics of Gravity and Regional Ocean Tides at Zhongshan Station in Antarctic

LIU Qingchao^1,2 , SUN Heping^1,2 , WANG Zemin³ , CHEN Xiaodong¹ , ZHANG Shengkai³ , XU Jianqiao¹     

Abstract: This study investigates the characteristics of gravity and ocean tides using time series recorded with the spring gravimeter and tide gauge at Zhongshan station, CATS2008 regional ocean tide model, and Eot11a global ocean tide model. The results show that the ocean tidal amplitude of O1 wave reaches 28 cm in the diurnal tidal bands. The sum of four tidals (Q1, O1, P1 and K1) difference between the global model and the tide gauge station is 4.2 cm, and that between the regional model and the tide gauge station is 4.4 cm. In contrast, the tidal amplitude of M2 wave reaches 20 cm in the semi-diurnal bands, and the sum of four tidals (N2, M2, S2 and K2) difference are 7.7 cm and 5.1 cm, respectively. This shows that correcting global ocean tide model using regional models is important. Subsequent to the correction of global oceanic loading using the regional model, the amplitudes of final residuals for the main gravity waves K1, M2, and S2 decreased by 9.84%, 56.14%, and 37.08%, respectively, which indicates that the regional ocean tide model can better reflect the true characteristics of ocean tides. Moreover, the validity of correcting global ocean tide model using the regional model is verified.

Key words: Antarctic     tidal gravity     regional and global ocean tide model    

南极冰盖质量的变化主要发生在浮冰区的边缘, 对全球海平面的变化起着至关重要的作用, 而且是全球气候变化的一个重要指标。利用测高卫星及GRACE重力卫星研究冰盖表面高或体积长期项变化时, 必须去掉观测资料中的海潮信号^[1-2]。

国际同行基于流体动力学方程、有限元计算 *** , 以验潮站、卫星测高观测资料等作为限制条件, 编制了一系列海潮模型^[3]。但南极地区缺少高精度的卫星测高和验潮站资料, 造成模型精度较低, 如文献[4]利用中山站重力资料测定的主要潮波振幅因子经全球海潮模型作负荷改正后与理论值的相对差异高达5%。验潮站可以提供高精度的局部海洋潮汐特征, 进而提高模型精度, 如文献[5]利用昭和验潮站资料代替全球海潮模型TPXO7.2中相应区域的潮高数据, 主要潮波的最终残差下降了5%。高精度的区域海潮模型可以用于修正全球海潮模型, 文献[6]基于11个全球海潮模型与中国东海及南海潮汐观测资料研究了近海潮汐资料对中国及邻区重力场的影响; 文献[7]发表了南极区域海潮模型CATS, 随着验潮站、GPS和ICESat等观测数据的积累, 模型精度进一步得到优化。

验潮站资料可用于评价海潮模型的适定性, 重力固体潮的潮汐参数有利于研究地球内部构造、地壳的运动等地球动力学现象^[8]。南极中山固体潮站于1989年2月26日开始建造, 坐落于东南极大陆的拉斯曼丘陵。1991年3月起, 中国科学院测量与地球物理研究所使用LCR型弹簧重力仪ET20/21在南极中山站(69.371°S, 76.365°E)进行连续重力观测^[9]。中山验潮站(69.383°S, 76.380°E)始建于1999年, 位于东南极的Nella峡湾, 由于每年的3月初至11月底一直被冰覆盖^[10], 因此, 使用了压力式验潮仪。因验潮仪老化等导致验潮站无法记录海洋潮汐资料, 2010年2月, 中国第26次南极科考队综合南极已有验潮站的建站 *** 及中山站的实际情况, 采用新的海底压力式验潮仪成功建设了我国在南极的首个永久性自动验潮站系统^[11], 为研究南极海洋潮汐特征提供实测数据。

本文首先对中山站2009年至2013年观测资料进行调和分析处理, 获取该地区重力及海洋潮汐参数; 其次, 利用验潮站及重力潮汐资料检验南极区域海潮模型CATS2008修正全球海潮模型Eot11a在南极地区的有效性; 最后, 利用修正后的海潮模型对重力潮汐参数作负荷效应改正。

1 数据处理1.1 重力固体潮

本文使用的重力固体潮资料为安装在中山站的弹簧重力仪ET21记录所得, 区间为2009年12月1日至2013年12月12日, 原始采样间隔为1 s; 气压观测数据由中山气象站提供, 原始采样间隔为10 min。为了便于处理观测资料, 分别采用整点采样法和三次样条内插法将重力及气压数据的采样间隔转换为1 min。由于弹簧重力仪为相对重力仪, 其输出值是与重力变化对应的电压值, 根据格值因子将其转化为重力单位。然后, 利用国际地球潮汐中心推荐的重力潮汐观测数据预处理软件“T-Soft”^[12], 移去理论合成潮, 根据标准大气重力导纳值去掉大气效应, 利用修正器删掉突跳、尖峰和由地震等引起的异常信号, 对由仪器故障、断电等造成的间断进行插值拟合, 利用IERS(international earth rotation service)提供的极移资料计算出极移重力效应^[13], 从而扣除观测资料中的极移效应, 再加上上文移除的理论合成潮及大气效应即得到预处理后的重力潮汐数据。最后, 采用低通数字滤波器将高密度采样资料变换为每小时采样。基于Eterna3.30标准分析软件^[14]对预处理后的重力潮汐序列作调和分析, 并利用台站气压与重力潮汐残差的线性拟合去除大气效应, 采用切比雪夫多项式拟合漂移, 获取高精度的重力潮汐参数(振幅因子和相位滞后), 具体数值结果列于表 1。其中, δ_th为文献[15]给出的基于旋转椭球体、非流体静力平衡地球建立的理论潮汐模型值。研究表明, O1和M2潮波振幅因子观测的相对精度分别为0.06%和0.08%, 但与理论值的相对差异分别为13.63%和19.50%, 这主要由海潮负荷效应引起, 因为中山站距海岸线很近, 仅10 m左右。

表 1 中山站重力潮汐参数Tab. 1 Gravity tidal parameters at Zhongshan station

      各潮波振幅是与纬度相关的, 长周期潮波振幅的更大值出现在两极地区, 振幅越大, 越容易从观测信号中分离出来^[16], 且求解精度越高。因此, 相对于其他纬度的台站, 当利用相同观测质量、相同观测区间的重力固体潮资料时, 中山站获得的长周期重力潮汐参数精度更高, 其中, MM(月潮)和MF(半月潮)为主要的长周期潮波, 表 1给出了对应的重力潮汐参数。数值结果表明, MM和MF潮波振幅因子观测的相对精度分别为9.05%和2.85%, 低于主要周日及半周日潮波的观测精度, 这是因为长周期潮波周期较长、受到干扰信号影响较大, 例如, 弹簧重力仪的零漂、南极冰盖的变化及其他环境因素等。由于MM和MF的观测相对精度较低, 在后续海潮模型修正的有效性检验中并未参与计算。

1.2 验潮站

本文分析的验潮资料观测区间为2010年2月21日至2012年12月13日, 采样间隔为1 h。压力式验潮仪给出的原始数据是压力值, 即验潮仪上方海水产生的压力和大气压的总和, 如图 1(a)所示。为了将其转化为真实的海面高, 首先在总压力资料中去除大气压得到海水产生的压力, 然后根据液压公式计算海面高h(如图 1(c)所示)
(1)

图 1 中山验潮站观测资料Fig. 1 Observation data at Zhongshan tide gauge station
 

式中, P是验潮仪记录的总压力; P_a为海面大气压(如图 1(b)所示), 由中山气象观测站提供; ρ为海水密度; G为中山站重力加速度。

某时刻的海面高(潮高)可以表示为^[17]

 (2)

式中, h(t_j)为t_j时刻的潮高; Z₀为常数; a为线性项系数; n为潮波数; A_k和g_k分别为潮波k的振幅与相位延迟; f_k(t_j)和u_k(t_j)分别为潮波k在t_j时刻振幅与相位的节点改正; V_k(t_j)为天文参数; R(t_j)为非潮汐残差信号; 节点改正及天文参数的计算见文献[18]。与传统调和分析相比, 该 *** 既消除了将天文参数假定为随时间线性变化而引起的误差, 又克服了将节点改正视为常数的问题, 且加入了随时间变化的线性趋势项。

令
 (3)

那么, 方程(2)转化为线性方程

 (4)

将海面高资料代入式(4), 基于versatile_tidana调和分析软件, 结合文献[19]给出的潮波选择顺序表, 利用奇异值分解(SVD)算法^[20]求出各分潮的调和常数(A_k和g_k), 可反映该地区的区域海洋潮汐特征, 具体数值结果列于表 2。在长周期频段, 潮波MF的海潮振幅更大, 数值为3.50 cm; 在周日频段, Q1、O1、P1和K1的振幅分别为6.55、28.21、8.82和27.37 cm, 其他潮波的振幅均未达到1.5 cm; 在半周日频段, N2、M2、S2和K2的振幅分别为5.68、19.99、17.92和5.03 cm, 而其他潮波均未超过1.7 cm。与文献[21]的研究结果相比, 主要潮波Q1、O1、P1、K1、N2、M2、S2和K2的振幅差异分别为0.13、0.21、0.28、0.23、0.20、0.13、0.28和0.09 cm, 这与观测时段有关, 也可归因于浅水非线性效应等。

表 2 中山验潮站各分潮的调和常数Tab. 2 Tidal harmonic constants of constituents at Zhongshan tide gauge station

2 全球海潮模型的修正

重力的海潮负荷效应可根据海潮潮高与重力负荷格林函数的褶积积分计算^[22-23], 当考虑目标台站的近海效应时, 用区域海潮模型修正全球海潮模型中相应区域的计算公式如下

 (5)

式中, L为负荷矢量; λ、φ为球坐标系下观测台站的经纬度; λ′、φ′分别是球坐标系下海洋潮汐负荷点的经纬度; H_L、H分别为近、远区海潮瞬时潮高; S_L和S_G-L分别代表台站的近和远海区域; ρ为海水密度; R为地球半径; G(ψ, A)为重力负荷格林函数。

全球海潮模型Eot11a是德国DGFI(Deutsches Geodätisches Forschungs Institut)以FES2004模型为基础, 依据1992年至2010年的TOPEX/Poseidon(T/P)、Jason-1/2、ERS-2和Envisat等测高卫星的调和分析资料作经验性改正而得到的模型, 分辨率为1/8°×1/8°(见文献[24])。CATS2008是包含整个南大洋的高精度南极区域海潮模型, 范围是56°S以南的整个区域, 网格分辨率为4 km×4 km, 以天文引潮位和全球海潮模型TPXO7.1的海潮高作为求解线性浅水方程的限制条件, 其同化数据包括:无冰区的T/P雷达测高资料, 50个高精度潮汐记录数据集(底部气压计、沿岸验潮站和GPS记录等, 具体信息参考Antarctic Tide Gauge Database)以及ICESat激光测高卫星获取的罗斯和龙尼-菲尔希纳冰架的表面高资料^[25]。其中, 卫星T/P和Jason-1/2可观测的最南端为66.2°S, ERS-2和Envisat为81.5°S, 而ICESat为86.0°S。而南极更大的两个冰架(龙尼-菲尔希纳和罗斯冰架)的浮冰区分别延伸到83.0°S和85.4°S, 因此只有ICESat测高卫星的观测范围可以覆盖整个南极冰盖的浮冰区, 为精化海潮模型提供必要的观测资料。

海潮模型Eot11a和CATS2008均给出了主要潮波在网格节点上的潮高信息(振幅及相位), 考虑到长周期潮波振幅较小且容易受区域环境因素的影响, 因此, 计算了4个周日(Q1、O1、P1和K1)及4个半周日(M2、N2、S2和K2)潮波潮高模型振幅与验潮站观测振幅之间的差异diff_i, 具体数值结果见表 3。此外, 还计算了周日潮与半日潮几个主要潮波潮高模型振幅与验潮站观测振幅之差的和, 分别用diff_D和diff_SD表示。由表 3可知, 在周日频段, 对于Q1和O1潮波, Eot11a模型值比CATS2008更接近潮汐的实际观测值, P1与K1潮波则相反, 但两模型的diff_D值基本相当。而在半周日频段中, 除N2潮波外, 其他潮波的CATS2008模型值比Eot11a更接近潮汐的实际观测值, 其中M2波的差值从3.7 cm(Eot11a-TGOTP)下降到1.6 cm(CATS2008-TGOTP)。

表 3 潮高模型振幅与验潮站观测振幅之差Tab. 3 The tidal amplitude's difference between model values and observations of the tide gauge

    根据式(5), 利用SPOTL软件^[26], 计算了南极中山站主要潮波的重力海潮负荷效应, 为方便后文讨论, 将Eot11a全球海潮模型记为模型1, 将Eot11a全球海潮模型经CATS2008区域海潮模型修正后记为模型2, 将验潮站潮高资料替换模型2中以验潮站为中心的1/16°×1/16°网格区域后记为模型3, 具体数值结果列于表 4。可以看出, 利用模型2和模型3计算得到的重力海潮负荷效应相差很小(小于0.001 μGal)。这是因为模型CATS2008同化了中山站附近的戴维斯验潮站的观测资料, 而模型1和模型2计算得到的海潮对K1、M2和S2潮波的负荷效应的振幅之差达到0.1 μGal, 其他主要潮波振幅之间的差异小于0.1 μGal, 相位差异更大值为3.8°(K2)。

表 4 中山站主要潮波的重力海潮负荷矢量Tab. 4 Gravity loading vectors for main waves at Zhongshan station

文献[27]利用超导重力仪器的潮汐观测资料研究了海潮模型的适定性问题, 为比较海潮模型CATS2008和Eot11a在南极地区的精度, 计算了主要潮波重力潮汐参数的最终残差, 其振幅越小, 说明模型精度越高。本文利用与验潮站相同观测区间的中山站固体潮观测资料获得重力潮汐参数, 根据矢量叠加原理, 减掉理论重力潮汐参数, 并作海潮负荷效应改正, 即得到各潮波的最终残差, 具体数值结果列于表 5, 振幅如图 2所示。结果表明, 在周日频段内, K1波经模型2作海潮负荷改正后, 最终振幅残差比经模型1作负荷改正时下降了0.8 nm/s², 下降幅度为9.84%, 其他潮波的最终残差变化不明显; 在半周日频段, N2、M2、S2和K2潮波的最终残差振幅均减小, 其中, M2和S2最终残差振幅分别下降了1.6和1.2 nm/s², 下降幅度分别为56.14%和37.08%。综上, 可证明南极区域海潮模型CATS2008计算的重力海潮负荷效应比全球海潮模型Eot11a的计算值更接近实际观测值, 即模型CATS2008修正模型Eot11a在南极地区是有效的。这是因为前者同化了南极地区多个高精度验潮站资料以及ICESat测高卫星获取的罗斯和龙尼-菲尔希纳冰架浮冰区的数据信息。但O1、K1潮波残差矢量仍然较大, 原因是其对应的理论振幅大, 分别达到了205和288 nm/s², 经模型2作海潮负荷效应改正得到的最终残差振幅与理论振幅的比值分别为3.24%和2.60%, 相应的, Q1、P1、N2、M2、S2和K2的比值分别为3.83%、3.24%、2.08%、1.29%、4.54%和2.45%。这表明经海潮负荷效应改正后, 各主要潮波观测振幅与理论值的相对差异在同一个量级。

表 5 中山站主要潮波最终残差矢量Tab. 5 The final residual vectors for main waves at Zhongshan station

图 2 中山站主要潮波最终残差振幅Fig. 2 Amplitudes of the final residual vectors for main waves at Zhongshan station
 

最后, 利用模型2对本文1.1节中相应的主要重力潮汐参数作海潮负荷效应改正, 得到主要长周期(MM和MF)、周日(O1、K1、P1和Q1)和半周日(M2、S2、N2和K2)潮波的振幅因子, 具体数值结果见表 1。相应的, 与文献[15]给出的理论振幅因子的相对差异分别为0.38%、4.28%、5.12%、4.63%、2.27%、3.31%、2.43%、3.39%、4.00%和1.77%, 这主要由海潮模型的不确定性、仪器观测误差及区域环境等引起。

3 结语

基于南极中山站重力潮汐及海洋潮汐观测序列, 采用调和分析 *** 精密确定了重力潮汐参数及潮高值, 分别反映了该地区的重力和海洋潮汐特征。计算结果表明, O1和M2潮波的重力振幅因子观测相对精度分别为0.06%和0.08%, 而MM和MF潮波的对应值分别为9.05%和2.85%, 说明虽然在两极地区长周期潮波振幅大, 但由于容易受到干扰信号(弹簧重力仪的漂移、冰盖的消融和凝固等)的影响, 造成观测精度不高。计算了全球海潮模型Eot11a和区域海潮模型CATS2008在中山验潮站的潮高模型值, 进而计算了潮高模型振幅与验潮站实际观测振幅的差值, 使用CATS2008对Eot11a相应区域作了修正(记为模型2), 并对其进行了有效性检验。研究表明, 在周日频段, 模型Eot11a和CATS2008提供的4个潮波潮高模型振幅与验潮站观测振幅之差的和基本相当, 而在半周日频段, 对应值分别为7.7 cm和5.1 cm, 与利用模型Eot11a作海潮负荷效应改正相比, 使用模型2时, Q1、K1、N2、M2、S2和K2潮波的最终残差振幅均减小, 其中, 潮波K1、M2和S2分别下降了0.82、1.55和1.16 nm/s², 证明了在南极地区用区域模型修正全球海潮模型的有效性, 因此, 今后在研究该地区的区域动力学问题时, 必须使用区域海潮模型修正后的全球模型计算海潮负荷效应。
 

【引文格式】刘清超, 孙和平, 王泽民, 等. 南极中山站的重力和区域海洋潮汐特征[J]. 测绘学报，2018，47(3)：316-323. DOI: 10.11947/j.AGCS.2018.20170123