留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

利用2006—2015年VLBI数据进行地球定向参数解算与分析

李斐 邵先远 曲春凯 鄢建国 郝卫峰 叶茂 杨轩 金炜桐

李斐, 邵先远, 曲春凯, 鄢建国, 郝卫峰, 叶茂, 杨轩, 金炜桐. 利用2006—2015年VLBI数据进行地球定向参数解算与分析[J]. 武汉大学学报 ● 信息科学版, 2019, 44(11): 1581-1587. doi: 10.13203/j.whugis20170365
引用本文: 李斐, 邵先远, 曲春凯, 鄢建国, 郝卫峰, 叶茂, 杨轩, 金炜桐. 利用2006—2015年VLBI数据进行地球定向参数解算与分析[J]. 武汉大学学报 ● 信息科学版, 2019, 44(11): 1581-1587. doi: 10.13203/j.whugis20170365
LI Fei, SHAO Xianyuan, QU Chunkai, YAN Jianguo, HAO Weifeng, YE Mao, YANG Xuan, JIN Weitong. Solution and Analysis of Earth Orientation Parameters with 2006-2015 VLBI Observation[J]. Geomatics and Information Science of Wuhan University, 2019, 44(11): 1581-1587. doi: 10.13203/j.whugis20170365
Citation: LI Fei, SHAO Xianyuan, QU Chunkai, YAN Jianguo, HAO Weifeng, YE Mao, YANG Xuan, JIN Weitong. Solution and Analysis of Earth Orientation Parameters with 2006-2015 VLBI Observation[J]. Geomatics and Information Science of Wuhan University, 2019, 44(11): 1581-1587. doi: 10.13203/j.whugis20170365

利用2006—2015年VLBI数据进行地球定向参数解算与分析

doi: 10.13203/j.whugis20170365
基金项目: 

国家自然科学基金 41374024

国家自然科学基金 41604004

国家自然科学基金 41174019

湖北省自然科学基金重点项目创新群体类项目 2015CFA011

宇航动力学国家重点实验室开放基金 2016ADL-DW0103

中国博士后科学基金 2016M602360

详细信息
    作者简介:

    李斐, 博士, 教授, 研究方向为重力场理论及应用。fli@whu.edu.cn

    通讯作者: 曲春凯, 硕士生。quchunkai999@163.com
  • 中图分类号: P228;P223

Solution and Analysis of Earth Orientation Parameters with 2006-2015 VLBI Observation

Funds: 

The National Natural Science Foundation of China 41374024

The National Natural Science Foundation of China 41604004

The National Natural Science Foundation of China 41174019

Hubei Province Natural Science Foundation Innovation Group Project 2015CFA011

the Open Research Fund of State Key Laboratory of Astronautic Dynamics 2016ADL-DW0103

China Postdoctoral Science Foundation 2016M602360

More Information
    Author Bio:

    LI Fei, PhD, professor, specializes in the theory and application of gravity field. E-mail:fli@whu.edu.cn

    Corresponding author: QU Chunkai, postgraduate.E-mail: quchunkai999@163.com
  • 摘要: 目前正处在下一代甚长基线干涉测量(very long baseline interferometry,VLBI)系统的建设时期。利用维也纳VLBI与卫星软件(Vienna VLBI and satellite software,VieVS)解算了2006—2015年的VLBI数据,得到了10 a的地球定向参数(Earth orientation parameters,EOP)时间序列,并与国际地球自转服务机构的结果进行了对比。利用解算结果得到了10 a的日长变化时间序列,通过傅里叶分析得出了日长变化的短周期、半月周期、月周期、半年周期和周年周期,同时还分析得到了极移序列中的周年项和张德勒周期项以及章动改正序列中的自由核章动项。此次解算工作可为武汉大学卫星台站日后的VLBI数据解析积累一定的经验。
  • 图  1  2006—2015年参与观测台站数目

    Figure  1.  Number of Observing Stations from 2006 to 2015

    图  2  EOP时间序列

    Figure  2.  Time Series of EOP

    图  3  ERP解算结果相对于IERS 08C04序列的改正

    Figure  3.  ERP Corrections to IERS 08C04

    图  4  去趋势项后的日长变化

    Figure  4.  ΔLOD Without the Trend Term

    图  5  日长变化幅频谱

    Figure  5.  Magnitude Spectrum of ΔLOD

    图  6  极移参数幅频谱

    Figure  6.  Magnitude Spectrum of Polar Motion Paremeters

    图  7  章动参数幅频谱

    Figure  7.  Magnitude Spectrum of Nutation Paremeters

    表  1  不同台站的观测次数

    Table  1.   Observation Times of Different Stations

    站名 观测次数
    AIRA 20
    ALGOPARK 27
    BADARY 262
    FORTLEZA 599
    HART15M 167
    HARTRAO 180
    HOBART12 279
    HOBART26 253
    ISHIOKA 7
    KASHIM11 9
    KASHIM34 6
    KATH12M 253
    KOKEE 774
    MATERA 408
    MEDICINA 142
    NYALES20 719
    ONSALA60 170
    PARKES 4
    RAEGYEB 36
    SEJONG 41
    SESHAN25 103
    SVETLOE 302
    TIANMA65 1
    TIGOCONC 749
    TSUKUB32 381
    URUMQI 16
    WARK12M 86
    WESTFORD 364
    WETTZ13N 28
    WETTZELL 950
    YARRA12M 249
    YEBES40M 143
    ZELENCHK 340
    下载: 导出CSV

    表  2  2006‒2015年EOP参数的内符合精度

    Table  2.   Internal Accord Accuracy of EOP During 2006-2015

    年份 Xpol/mas Ypol/mas UT1-UTC/ms dX/mas dY/mas
    2006 0.065 0.065 0.003 0.048 0.049
    2007 0.051 0.057 0.002 0.040 0.042
    2008 0.069 0.071 0.003 0.051 0.051
    2009 0.064 0.070 0.003 0.052 0.050
    2010 0.083 0.073 0.003 0.050 0.050
    2011 0.074 0.068 0.003 0.046 0.046
    2012 0.066 0.060 0.003 0.040 0.040
    2013 0.073 0.068 0.003 0.046 0.045
    2014 0.060 0.054 0.003 0.041 0.040
    2015 0.045 0.048 0.003 0.028 0.028
    下载: 导出CSV

    表  3  2006-2015年EOP参数的外符合精度

    Table  3.   External Accord Accuracy of EOP During 2006-2015

    年份 Xpol/mas Ypol/mas UT1-UTC/ms dX/mas dY/mas
    2006 0.115 0.146 0.007 0.149 0.216
    2007 0.129 0.147 0.007 0.102 0.125
    2008 0.159 0.137 0.008 0.104 0.149
    2009 0.195 0.143 0.008 0.094 0.093
    2010 0.189 0.163 0.009 0.103 0.095
    2011 0.161 0.192 0.010 0.098 0.148
    2012 0.165 0.181 0.009 0.120 0.113
    2013 0.207 0.252 0.011 0.112 0.106
    2014 0.202 0.209 0.014 0.089 0.128
    2015 0.140 0.175 0.011 0.082 0.088
    下载: 导出CSV

    表  4  2006-2008年ERP参数的外符合精度

    Table  4.   External Accord Accuracy of ERP During 2006—2008

    方法 年份 Xpol/mas Ypol/mas UT1-UTC/ms
    文献[9-10] 2006 0.338 0.320 0.176
    2007 0.338 0.320 0.176
    2008 0.450 0.290 1.608
    本文方法 2006 0.092 0.112 0.005
    2007 0.097 0.100 0.005
    2008 0.117 0.109 0.006
    下载: 导出CSV
  • [1] 王广利, 钱志瀚.1979—1999全球天测与测地VLBI观测资料处理结果与讨论[J].天文学报, 2000, 41 (3): 320-326 doi:  10.3321/j.issn:0001-5245.2000.03.011

    Wang Guangli, Qian Zhihan. Solution and Discussion of Global Astrometric and Geodetic VLBI Data (1979—1999) Reduction[J]. Acta Astronomica Sinica, 2000, 41(3):320-326 doi:  10.3321/j.issn:0001-5245.2000.03.011
    [2] 李金岭, 王广利. VLBI观测得到的极移序列分析[J].科学通报, 2000, 45(12): 1 260-1 262 doi:  10.3321/j.issn:0023-074X.2000.12.005

    Li Jinling, Wang Guangli. An Analysis of the Polar Motion Series from VLBI Observations[J]. Chinese Sci Bull, 2000, 45(12): 1 260-1 262 doi:  10.3321/j.issn:0023-074X.2000.12.005
    [3] Böhm J, Böhm S, Nilsson T, et al. The New Vienna VLBI Software VieVS[C]. IAG 2009 Scientific Assembly, Buenos Aires, Argentina, 2009
    [4] Wei Erhu, Tang Shenquan. On Module Framework of VieVS and Data Processing for Teaching and Research[J]. Coordinates, 2015, 11(2): 50-54
    [5] 马高峰. VLBI2010与GNSS联合数据分析理论及方法研究[D].郑州: 信息工程大学, 2011

    Ma Gaofeng. Research on the Theory and Method of Data Analysis Combined with VLBI2010 and GNSS[D]. Zhengzhou: Information Engineering University, 2011
    [6] Belda S, Heinkelmann R, Ferrandiz J M, et al. On the Consistency of the Current Conventional EOP Series and the Celestial and Terrestrial Reference Frames[J]. Journal of Geodesy, 2017, 91(2): 135-149 doi:  10.1007/s00190-016-0944-3
    [7] 张波.基于VLBI资料的地球定向参数高频变化研究[D].上海: 中国科学院上海天文台, 2004

    Zhang Bo. A Study of the High Frequency Variation of Earth Rotation Parameters Based on VLBI Obervation[D]. Shanghai: Shanghai Astronomical Observatory, Chinese Academy of Sciences, 2004
    [8] Wei Erhu, Liu Jingnan, Shi Chuang. Estimation of Geodetic Parameters with VLBI Data of Last 5 Years[J]. Geo⁃spatial Information Science, 2007, 10(1): 12-16 http://d.old.wanfangdata.com.cn/Periodical/dqkjxxkxxb-e200701003
    [9] 魏二虎, 刘经南, 潘培鲸.近三年VLBI观测的数据处理及分析[J].武汉大学学报∙信息科学版, 2008, 33(12): 1 275-1 278 http://ch.whu.edu.cn/CN/abstract/abstract1777.shtml

    Wei Erhu, Liu Jingnan, Pan Peijing. On Data Processing with Last Three Years of VLBI Observation[J]. Geomatics and Information Science of Wuhan University, 2008, 33(12): 1 275-1 278 http://ch.whu.edu.cn/CN/abstract/abstract1777.shtml
    [10] 魏二虎, 李雪川, 易慧, 等.利用2008年VLBI观测数据进行EOP参数的解算与分析[J].武汉大学学报∙信息科学版, 2010, 35(8): 988-990 http://ch.whu.edu.cn/CN/abstract/abstract1035.shtml

    Wei Erhu, Li Xuechuan, Yi Hui, et al. On EOP Partameter Based on 2008 VLBI Observation Data[J]. Geomatics and Information Science of Wuhan University, 2010, 35(8): 988-990 http://ch.whu.edu.cn/CN/abstract/abstract1035.shtml
    [11] 张志斌, 王广利, 刘祥, 等.中国VLBI网观测地球定向参数能力分析[J].武汉大学学报∙信息科学版, 2013, 38(8): 911-915 http://ch.whu.edu.cn/CN/abstract/abstract2711.shtml

    Zhang Zhibin, Wang Guangli, Liu Xiang, et al. Analysis of EOP Determination via Chinese VLBI Network[J]. Geomatics and Information Science of Wuhan University, 2013, 38(8): 911-915 http://ch.whu.edu.cn/CN/abstract/abstract2711.shtml
    [12] 魏二虎, 张奇.利用1985—2010年VLBI数据进行日长变化分析[J].武汉大学学报∙信息科学版, 2012, 37(9): 1 020-1 023 http://ch.whu.edu.cn/CN/abstract/abstract305.shtml

    Wei Erhu, Zhang Qi.On Length of Day with 1985—2010 VLBI Observation[J]. Geomatics and Information Science of Wuhan University, 2012, 37(9):1 020-1 023 http://ch.whu.edu.cn/CN/abstract/abstract305.shtml
    [13] 韩永刚, 李志安.地球自转速率变化主要中短周期的分析[J].地球物理学进展, 2002, 17(2): 349-352 doi:  10.3969/j.issn.1004-2903.2002.02.024

    Han Yonggang, Li Zhian. Analysis of Main Middle and Short Periods in the Variation of the Earth's Ratation[J]. Progress in Geophysics, 2002, 17(2): 349-352 doi:  10.3969/j.issn.1004-2903.2002.02.024
    [14] 徐建桥, 孙和平, 罗少聪.利用国际超导重力仪观测资料研究地球自由核章动[J].中国科学(D辑), 2001, 31(9): 719-726 http://d.old.wanfangdata.com.cn/Periodical/zgkx-cd200109003

    Xu Jianqiao, Sun Heping, Luo Shaocong. Study of the Earth's Free Core Nutation by Tidal Gravity Data Recorded with International Superconducting Gravimeters[J]. Science in China (Ser D), 2001, 31(9): 719-726 http://d.old.wanfangdata.com.cn/Periodical/zgkx-cd200109003
  • [1] 张宇, 郭丽, 王广利, 李金岭, 黄勇, 李培佳, 郑为民, 王文彬, 舒逢春, 刘庆会.  定位归算在嫦娥四号中继星变轨阶段的应用 . 武汉大学学报 ● 信息科学版, 2020, 45(2): 189-195. doi: 10.13203/j.whugis20180365
    [2] 李金岭, 赵达新, 柳聪, 张会, 虞林峰, 赵融冰, 王锦清, 范庆元.  佘山13 m口径射电望远镜指向扫描数据解析 . 武汉大学学报 ● 信息科学版, 2020, 45(2): 159-166. doi: 10.13203/j.whugis20180415
    [3] 张志斌, 王广利, 宋世泽, 张阿丽, 王宏.  归心测量中轮轨式VLBI天线的水平钢轨形变建模 . 武汉大学学报 ● 信息科学版, 2019, 44(10): 1498-1504. doi: 10.13203/j.whugis20180050
    [4] 孙和平, 刘清超, 崔小明, 徐建桥.  利用重力和VLBI技术检测地球液态核的动力学效应 . 武汉大学学报 ● 信息科学版, 2018, 43(12): 2058-2063. doi: 10.13203/j.whugis20180110
    [5] 魏二虎, 汤深权, 刘建栋, 董翠军.  基于联邦卡尔曼滤波的月球车联合定位方法 . 武汉大学学报 ● 信息科学版, 2017, 42(10): 1431-1437. doi: 10.13203/j.whugis20160137
    [6] 李金岭, 黄飞, 孙中苗.  Ka波段天测与测地VLBI观测的利弊分析 . 武汉大学学报 ● 信息科学版, 2017, 42(2): 257-262. doi: 10.13203/j.whugis20140723
    [7] 魏二虎, 刘文杰, WEI Jianan, 金双根, 刘经南.  VLBI和GPS观测联合解算地球自转参数和日长变化 . 武汉大学学报 ● 信息科学版, 2016, 41(1): 66-71,92. doi: 10.13203/j.whugis20130435
    [8] 刘鹂, 李金岭.  中国VLBI网测站参数解析 . 武汉大学学报 ● 信息科学版, 2014, 39(3): 262-266. doi: 10.13203/j.whugis20120005
    [9] 李金岭, 张津维.  利用GPS测量监测VLBI天线参考点的仿真分析 . 武汉大学学报 ● 信息科学版, 2013, 38(12): 1387-1391.
    [10] 李金岭, 郭丽, 孙中苗.  关于脉冲星参数的我国VLBI网测量 . 武汉大学学报 ● 信息科学版, 2013, 38(4): 386-389.
    [11] 张志斌, 王广利, 刘祥, 唐正宏.  中国VLBI网观测地球定向参数能力分析 . 武汉大学学报 ● 信息科学版, 2013, 38(8): 911-915.
    [12] 魏二虎, 张奇.  利用1985~2010年VLBI数据进行日长变化分析 . 武汉大学学报 ● 信息科学版, 2012, 37(9): 1020-1023.
    [13] 李金岭, 魏二虎, 孙中苗, 钱志瀚.  关于我国天测与测地VLBI网络未来建设的讨论 . 武汉大学学报 ● 信息科学版, 2010, 35(6): 670-673.
    [14] 魏二虎, 周培来, 刘翠芝, 阿里根.  利用RADIOASTRON模拟估计的EOP对相关参数精度的灵敏度分析 . 武汉大学学报 ● 信息科学版, 2010, 35(6): 674-678.
    [15] 魏二虎, 李雪川, 易慧, 刘经南.  利用2008年VLBI观测数据进行EOP参数的解算与分析 . 武汉大学学报 ● 信息科学版, 2010, 35(8): 988-990.
    [16] 魏二虎, 刘经南, 潘培鲸.  近三年VLBI观测的数据处理及分析 . 武汉大学学报 ● 信息科学版, 2008, 33(12): 1275-1278.
    [17] 朱新慧, 孙付平, 秦勇.  组合GPS和VLBI数据建立板块运动模型 . 武汉大学学报 ● 信息科学版, 2005, 30(7): 604-608.
    [18] 张捍卫, 郑勇, 杜兰.  河外射电源地面VLBI观测的相对论时延模型 . 武汉大学学报 ● 信息科学版, 2003, 28(5): 585-588.
    [19] 郭丽, 李金岭, 王广利, 黄逸丹, 王文彬, 郑为民.  嫦娥五号探测器月面定位与精度分析 . 武汉大学学报 ● 信息科学版, 0, 0(0): -. doi: 10.13203/j.whugis20210226
    [20] 钟胜坚, 舒逢春, 李金岭, 吴徳, 何旋, 童锋贤.  利用中国VGOS站开展国际联测的EOP精度仿真分析 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20210204
  • 加载中
图(7) / 表(4)
计量
  • 文章访问数:  1164
  • HTML全文浏览量:  76
  • PDF下载量:  332
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-04-25
  • 刊出日期:  2019-11-05

利用2006—2015年VLBI数据进行地球定向参数解算与分析

doi: 10.13203/j.whugis20170365
    基金项目:

    国家自然科学基金 41374024

    国家自然科学基金 41604004

    国家自然科学基金 41174019

    湖北省自然科学基金重点项目创新群体类项目 2015CFA011

    宇航动力学国家重点实验室开放基金 2016ADL-DW0103

    中国博士后科学基金 2016M602360

    作者简介:

    李斐, 博士, 教授, 研究方向为重力场理论及应用。fli@whu.edu.cn

    通讯作者: 曲春凯, 硕士生。quchunkai999@163.com
  • 中图分类号: P228;P223

摘要: 目前正处在下一代甚长基线干涉测量(very long baseline interferometry,VLBI)系统的建设时期。利用维也纳VLBI与卫星软件(Vienna VLBI and satellite software,VieVS)解算了2006—2015年的VLBI数据,得到了10 a的地球定向参数(Earth orientation parameters,EOP)时间序列,并与国际地球自转服务机构的结果进行了对比。利用解算结果得到了10 a的日长变化时间序列,通过傅里叶分析得出了日长变化的短周期、半月周期、月周期、半年周期和周年周期,同时还分析得到了极移序列中的周年项和张德勒周期项以及章动改正序列中的自由核章动项。此次解算工作可为武汉大学卫星台站日后的VLBI数据解析积累一定的经验。

English Abstract

李斐, 邵先远, 曲春凯, 鄢建国, 郝卫峰, 叶茂, 杨轩, 金炜桐. 利用2006—2015年VLBI数据进行地球定向参数解算与分析[J]. 武汉大学学报 ● 信息科学版, 2019, 44(11): 1581-1587. doi: 10.13203/j.whugis20170365
引用本文: 李斐, 邵先远, 曲春凯, 鄢建国, 郝卫峰, 叶茂, 杨轩, 金炜桐. 利用2006—2015年VLBI数据进行地球定向参数解算与分析[J]. 武汉大学学报 ● 信息科学版, 2019, 44(11): 1581-1587. doi: 10.13203/j.whugis20170365
LI Fei, SHAO Xianyuan, QU Chunkai, YAN Jianguo, HAO Weifeng, YE Mao, YANG Xuan, JIN Weitong. Solution and Analysis of Earth Orientation Parameters with 2006-2015 VLBI Observation[J]. Geomatics and Information Science of Wuhan University, 2019, 44(11): 1581-1587. doi: 10.13203/j.whugis20170365
Citation: LI Fei, SHAO Xianyuan, QU Chunkai, YAN Jianguo, HAO Weifeng, YE Mao, YANG Xuan, JIN Weitong. Solution and Analysis of Earth Orientation Parameters with 2006-2015 VLBI Observation[J]. Geomatics and Information Science of Wuhan University, 2019, 44(11): 1581-1587. doi: 10.13203/j.whugis20170365
  • 甚长基线干涉测量(very long baseline interferometry,VLBI)是至今唯一能够同时提供天球参考框架(celestial reference frame, CRF)、地球参考框架(terrestrial reference frame, TRF)以及地球定向参数(Earth orientation parameters, EOP)的空间测地技术,在天体测量、大地测量、地球物理等相关领域应用较为广泛[1]。自20世纪80年代至今,由于VLBI的高稳定性与高精度等特点,其一直是测定EOP的主要支撑技术[2]。武汉大学诗琳通地球空间信息科学国际研究中心卫星台站(简称武汉大学卫星台站)已经立项,将观测能力拓展至测地VLBI观测。现已建成并投入使用的13 m射电望远镜的转速高于10°/s,目前已完成制冷接收机改造,并配备氢原子钟和水汽辐射计,基本符合VLBI2010全球观测系统(the VLBI2010 global observing system, VGOS)台站的规格要求。武汉大学卫星台站有望参与中国区域VGOS组网,本文将为该站的VLBI数据解析积累经验。

    维也纳VLBI与卫星软件(Vienna VLBI and satellite software,VieVS)作为专业的VLBI数据解析软件,采用国际地球自转服务机构(International Earth Rotation Service,IERS)的最新会议标准,包含会议推荐的最新计算模型,具有良好的人机交互界面,使用方便[3-4]。本文主要利用VieVS软件解算了2006—2015年的VLBI数据,得到了10 a的EOP时间序列,并与IERS的官方结果进行了精度对比,利用频谱分析方法对3次样条插值后的等间隔EOP时间序列进行了周期分析。

    • 本文采用国际地球参考框架ITRF 2008和国际天球参考框架ICRF2作为解算基准,选用的物理解算模型与IERS标准(2010)相一致。具体模型参数如下:采用喷气推进实验室DE421星历计算太阳系天体位置;从NGS文件中读取计算对流层时延所用的压力和温度参数,缺失值采用GPT2模型计算,对流层映射函数选取VMF1;电离层时延信息从NGS文件中读取;台站坐标改正考虑固体潮、极潮、海洋极潮、海洋潮负荷、大气潮负荷、大气非潮汐负荷、天线热变形等因素的影响;EOP的先验时间序列采用IERS 08C04,岁差章动模型选用IAU2006/2000A,插值方法选择拉格朗日插值。

    • 本文使用2006—2015年的观测数据进行解算,数据来源于国际大地测量与天体测量VLBI服务组织(International VLBI Service for Geodesy and Astrometry,IVS)的网站(ftp://cddisa.gsfc.nasa.gov),数据格式为NGS格式。此数据源于监测地球自转的两个24 h的定期观测计划,分别是每周一和周四观测的IVS-R1和IVS-R4计划[5]。2006—2015年,10 a一共进行了1 029次R1/R4观测,其中,2006-12-07的观测(其观测文件为06DEC07XE_N005)只有两个台站的数据,解算结果变差,因此本文解算时未包含此次观测。

      参与台站的数量随时间的变化如图 1所示,可以看出,2010年平均参与台站数最少,可能会影响解算的精度。

      图  1  2006—2015年参与观测台站数目

      Figure 1.  Number of Observing Stations from 2006 to 2015

      解算涉及全球范围内的33个台站,表 1统计了所有台站参与观测的次数。

      表 1  不同台站的观测次数

      Table 1.  Observation Times of Different Stations

      站名 观测次数
      AIRA 20
      ALGOPARK 27
      BADARY 262
      FORTLEZA 599
      HART15M 167
      HARTRAO 180
      HOBART12 279
      HOBART26 253
      ISHIOKA 7
      KASHIM11 9
      KASHIM34 6
      KATH12M 253
      KOKEE 774
      MATERA 408
      MEDICINA 142
      NYALES20 719
      ONSALA60 170
      PARKES 4
      RAEGYEB 36
      SEJONG 41
      SESHAN25 103
      SVETLOE 302
      TIANMA65 1
      TIGOCONC 749
      TSUKUB32 381
      URUMQI 16
      WARK12M 86
      WESTFORD 364
      WETTZ13N 28
      WETTZELL 950
      YARRA12M 249
      YEBES40M 143
      ZELENCHK 340

      由于IERS 08C04是在ITRF 2008框架下解算的,当解算台站坐标时,站坐标与ITRF 2008中的定义会存在差异,这将导致解算得到的EOP序列与IERS 08C04序列的差异变大[6]。考虑到本文研究的主要目的是解算高精度的EOP,因此采用单次解[7],即解算时固定台站坐标和射电源坐标[8],再对地球定向参数、大气参数和测站钟差等进行解算。

      解算过程中采取的主要策略有:首先,通过解算的时延残差分布发现台站的钟跳问题,选取某稳定台站(例如WETTZELL站)的时钟作为参考时钟,在OPT文件中设置问题台站的时钟间断后重新解算;然后,剔除残差异常的台站和基线后重新解算,例如TIGOCONC台站在2010年2月发生地震,则需剔除该站;最后,通过OUT文件剔除观测数据中的坏点后重新解算。

      解算共得到了2 879 d的EOP时间序列,如图 2所示。图 2中,dX与dY为两个章动参数, ${X_{{\rm{pol}}}}$和${Y_{{\rm{pol}}}}$为两个极移参数, UT1为世界时1(universal time 1), UTC为协调世界时(coordinated universal time), UT1-UTC表示UT1和UTC之差。下面从两个方面对解算结果进行分析讨论。

      图  2  EOP时间序列

      Figure 2.  Time Series of EOP

      1)内符合精度。将VieVS解算结果的形式误差取绝对值后再取平均,作为数据解析的内符合精度,5个EOP参数的内符合精度见表 2。由表 2可以看出,不同年份解算结果的内符合精度之间存在轻微差异,但都在同一个数量级。轻微的差异可能与当年观测资料质量有关,另外OPT文件中的解算策略也会对结果造成影响。对各个EOP参数2006—2015年的内符合精度分别取平均值,得到${X_{{\rm{pol}}}}$和${Y_{{\rm{pol}}}}$的平均精度分别为0.065 mas和0.063 mas,UT1-UTC的平均精度为0.003 ms,dX和dY的平均精度均为0.044 mas。由表 2可知,2010年的解算精度最差,这与当年平均参与台站数最少的情况相吻合;相反,2015年的解算精度最好,这与当年平均参与台站数最多的情况相一致。

      表 2  2006‒2015年EOP参数的内符合精度

      Table 2.  Internal Accord Accuracy of EOP During 2006-2015

      年份 Xpol/mas Ypol/mas UT1-UTC/ms dX/mas dY/mas
      2006 0.065 0.065 0.003 0.048 0.049
      2007 0.051 0.057 0.002 0.040 0.042
      2008 0.069 0.071 0.003 0.051 0.051
      2009 0.064 0.070 0.003 0.052 0.050
      2010 0.083 0.073 0.003 0.050 0.050
      2011 0.074 0.068 0.003 0.046 0.046
      2012 0.066 0.060 0.003 0.040 0.040
      2013 0.073 0.068 0.003 0.046 0.045
      2014 0.060 0.054 0.003 0.041 0.040
      2015 0.045 0.048 0.003 0.028 0.028

      2)外符合精度。将解算的EOP时间序列与IERS网站公布的EOP 08C04(IAU2000)时间序列作差,差值的均方根作为外符合精度,各参数的外符合精度见表 3。由表 3可见,本文解算的EOP结果与IERS官方结果接近。通过计算可得:XpolYpol的平均精度分别为0.169、0.178 mas,UT1-UTC的平均精度为0.010 ms,dX和dY的平均精度分别为0.107、0.131 mas。

      表 3  2006-2015年EOP参数的外符合精度

      Table 3.  External Accord Accuracy of EOP During 2006-2015

      年份 Xpol/mas Ypol/mas UT1-UTC/ms dX/mas dY/mas
      2006 0.115 0.146 0.007 0.149 0.216
      2007 0.129 0.147 0.007 0.102 0.125
      2008 0.159 0.137 0.008 0.104 0.149
      2009 0.195 0.143 0.008 0.094 0.093
      2010 0.189 0.163 0.009 0.103 0.095
      2011 0.161 0.192 0.010 0.098 0.148
      2012 0.165 0.181 0.009 0.120 0.113
      2013 0.207 0.252 0.011 0.112 0.106
      2014 0.202 0.209 0.014 0.089 0.128
      2015 0.140 0.175 0.011 0.082 0.088

      EOP 08C04时间序列是由IERS综合全球导航卫星系统、甚长基线干涉测量、卫星激光测距、星基多普勒轨道确定和无线电定位组合系统这4种技术的结果得到的。文献[9-10]的EOP时间序列仅来自VLBI数据,并与本文结果存在3 a(2006—2008年)的时间重合,因此将2006—2008年的结果与文献[9-10]作对比,验证本文解算的可靠性。考虑到本文采用了IAU2000A章动模型,而文献[9-10]采用的是IAU1980章动模型,这两个章动模型之间存在较大差异,不具有可比性,因此仅对比地球自转参数(Earth rotation parameters,ERP)的解算结果,即XpolYpol、UT1-UTC。为方便对比,3 a解算结果的外符合精度采用文献[9-10]中的定义方式,即解算结果与官方结果相减后取绝对值再取平均,解算的3 a ERP参数的外符合精度如表 4所示。

      表 4  2006-2008年ERP参数的外符合精度

      Table 4.  External Accord Accuracy of ERP During 2006—2008

      方法 年份 Xpol/mas Ypol/mas UT1-UTC/ms
      文献[9-10] 2006 0.338 0.320 0.176
      2007 0.338 0.320 0.176
      2008 0.450 0.290 1.608
      本文方法 2006 0.092 0.112 0.005
      2007 0.097 0.100 0.005
      2008 0.117 0.109 0.006

      表 4表明,本文解算的ERP参数明显优于参考文献[9-10]的结果。文献[11]使用了VieVS软件解算了佘山站和南山站在2011年之前参与的所有IVS 24 h观测数据,其解算的极移改正误差一般在1 mas内,UT1-UTC误差一般在20 µs以内,本文解算精度与其相当。图 3为本文ERP解算结果相对于IERS 08C04序列的改正。

      图  3  ERP解算结果相对于IERS 08C04序列的改正

      Figure 3.  ERP Corrections to IERS 08C04

    • 地球自转频率的变化通常以日长(length of day,LOD)的观测量相对于标准日长(86 400 s)的变化ΔLOD来表示。根据本文VLBI解算得到的UT1 - UTC,通过式(1)计算得到ΔLOD,其中UTC与国际原子时(international atomic time,TAI)仅存在跳秒的差异。当ΔLOD>0时,日长变短,地球自转速度加快;反之,日长变长,地球自转速度减慢[12]

      $$ {\Delta _{{\rm{LOD}}}} = \frac{{{{\left( {{T_{{\rm{UT}}1}} - {T_{{\rm{TAI}}}}} \right)}_t} - {{\left( {{T_{{\rm{UT}}1}} - {T_{{\rm{TAI}}}}} \right)}_{t + {\rm{\Delta }}t}}}}{{{\rm{\Delta }}t}}$$ (1)

      式中,ΔLOD表示日长的观测量相对于标准日长的变化;TUT1表示世界时;TTAI表示国际原子时;t表示某一时刻;Δt表示时间间隔。

      由于傅里叶变换只能处理等间隔采样的数据,而且数据的趋势项会影响周期项在幅频谱中的振幅,因此需要将得到的非等间隔ΔLOD数据进行3次样条插值,插值后的时间间隔为1 d,去除序列的平均值,以便通过频谱分析提取ΔLOD的周期特征[12]。去趋势项后的ΔLOD序列如图 4所示,对该序列进行傅里叶变换得到幅频谱如图 5所示。图 5中,X表示与横坐标频率对应的周期,Y表示纵坐标振幅。

      图  4  去趋势项后的日长变化

      Figure 4.  ΔLOD Without the Trend Term

      图  5  日长变化幅频谱

      Figure 5.  Magnitude Spectrum of ΔLOD

      根据频率与周期的倒数关系,图 5的幅频谱中,从左至右6个振幅峰值对应的周期分别为3 650.00 d、365.00 d、182.50 d、27.65 d、13.67 d、9.13 d。将本文所得周期结果与文献[12]的结果作比较,发现周年周期、半年周期、月周期和半月周期的符合度很高,两者结果最大差值仅在0.05 d以内;9.13 d的最小周期与文献[13]通过最大熵谱分析法得到的周期相一致;本文所得最大周期为3 650.00 d,约10 a,此周期与所选时间序列的长度有关,并非日长变化的真实周期。

    • 采用与日长变化相同的方式,对XpolYpol插值并去除趋势项后进行频谱分析,其幅频谱分别如图 6(a)图 6(b)所示。

      图  6  极移参数幅频谱

      Figure 6.  Magnitude Spectrum of Polar Motion Paremeters

      图 6中可以发现,极移参数XpolYpol存在3个明显的频率峰值,对应的周期分别为456.25 d、405.56 d、365.00 d,前两个周期与427 d的张德勒周期项有关,第3个周期为周年项,其中周年项的能量最为显著。

    • 采用与日长变化相同的方式,对dX和dY插值并去除趋势项后进行频谱分析,其幅频谱分别如图 7(a)图 7(b)所示。从图 7中可以看到,章动参数dX和dY存在两个明显的频率峰值,对应的周期分别为456.25 d和405.56 d,与429 d的自由核章动本征周期相吻合[14],此周期项的存在是因为IAU2000A章动模型未包含自由核章动。

      图  7  章动参数幅频谱

      Figure 7.  Magnitude Spectrum of Nutation Paremeters

    • 本文解算了2006—2015年的VLBI数据,得到了10 a跨度的EOP时间序列,并与IERS 08C04的结果进行了比较。结果表明,内符合精度和外符合精度均在误差范围之内。同时,利用解算的UT1-UTC结果得到了10 a的ΔLOD时间序列,并通过频谱分析得出了ΔLOD时间序列中的多个周期,包括9.13 d的短周期、13.67 d的半月周期、27.65 d的月周期、182.50 d的半年周期以及365.00 d的周年周期。此外,还分析出了极移序列中的周年项和张德勒周期项以及章动改正序列中的自由核章动项。本文VLBI解算工作将为武汉大学卫星台站加入VGOS网后的VLBI数据解析积累一定经验。

      未来希望通过解析中国区域VLBI数据,获取中国独立自主的EOP时间序列,从而摆脱对国外观测数据的依赖。

参考文献 (14)

目录

    /

    返回文章
    返回