IERS EOP 文件的解读

一、预备知识

IERS（International Earth Rotation Service）即国际地球自转服务，于1988年由国际大地测量学和地球物理学联合会（International Union of Geodesy and Geophysics， IUGG)与国际天文学联合会（International Astronomical Union，IAU）共同建立，用以取代国际时间局（BIH）的地球自转部分和原有的国际极移服务（IPMS），于1988年1月1日正式开始运作。

IERS 通过分布在全球各地的观测网获得各种观测数据数据，IERS 中心局根据各分析中心的处理结果进行综合分析，得出 ICRF、ITRF 和 EOP 的最终结果，由 IERS 年度报告和技术备忘录向全世界发布。

• ICRF: The International Celestial Reference Frame (ICRF) consists of equatorial coordinates of extragalactic radio sources observed with VLBI, realizing an ideal reference system, the International Celestial Reference System (ICRS). （即国际天体参考框架（ICRF）由VLBI观察到的赤道式无线电源组成，实现了理想的参考系统，即国际天体参考系统（ICRS）。）
• The International Terrestrial Reference Frame (ITRF) is a set of points with their 3-dimensional cartesian coordinates which realize an ideal reference system, the International Terrestrial Reference System (ITRS).（国际地面参考框架（ITRF）是一组具有三维笛卡尔坐标的点，实现了理想的参考系统，国际陆地参考系统（ITRS）。）
• EOP is Earth Orientation Parameter(地球定位定向参数), IERS provides rapid, monthly and long term earth orientation data as well as leap second announcements and announcements of DUT1 in the form of bulletins or data files.（IERS以公告或数据文件的形式提供了快速，每月和长期地球定向数据以及DUT1的LEAP第二公告和DUT1的公告。）

1、数据下载

官网：https://www.iers.org/
数据下载：https://datacenter.iers.org/products/

dir: products/eop/rapid/standard/

下载到的数据样例：

2、一些概念

• 极移：地球自转轴相对于地球本体的位置是变化的，这种运动称地极移动，简称极移。
• 章动：月球绕地球旋转的轨道称为白道，月球运行的轨道与地球的之间距离是不断变化的，使得月球引力产生的大小和方向不断变化，从而导致北天极在天球上绕黄极旋转的轨道不是平滑的小圆，而是类似圆的波浪曲线运动，即地球旋转轴在岁差的基础上叠加周期为18.6年，且振幅为9.21″的短周期运动。这种现象称为章动。
• 岁差：由于日、月等天体的影响，地球的旋转轴在空间围绕黄极发生缓慢旋转，类似于旋转陀螺，形成一个倒圆锥体，其锥角等于黄赤交角ε=23.5°，旋转周期为26000年，这种运动称为岁差，是地轴方向相对于空间的长周期运动。岁差使春分点每年向西移动50.3″。

3、单位换算

• 极移经常用到的单位是角秒（sec. of arc，arc sec）
• 角分：arc minute；1 角分 = 1/60 度 = 60 角秒
• $180°=\pi rad$，那么 $1rad\approx 57.3°$，$1°\approx 0.01745rad$
• 1 arcsecond = 1 / 60 arcminute = 1 / 3600 degree
• 1 毫角秒 = $10^{-3}$ 角秒；1 微($\mu$)角秒 = $10^{-6}$ 角秒
• 1 arcsec = 4.8481368 μrad

Since one degree is defined as one three hundred sixtieth (1/360) of a circle, 1 minute of arc is 1/21,600 of the amount of arc in a closed circle.
1 arcsecond = 1 / 60 arcminute = 1 / 3600 degree
—— 摘自维基百科
ps: arc sec 直译为弧秒，实际上和弧度半毛钱关系都没有，之前笔者也是被这个翻译所迷惑了，现在笔者觉得确切的翻译应该为角秒。关于弧度，没有 60 进制的小单位转换！

二、数据解读

官方说明文档：https://www.iers.org/IERS/EN/DataProducts/EarthOrientationData/eop.html
https://maia.usno.navy.mil/ser7/readme.finals2000A

finals2000A.all contains the values from IERS Bulletin A for x/y pole, UT1-UTC, LOD, dX, dY, their errors and predictions for next 365 days for x/y pole, UT1-UTC, dX and dY as well as the values from IERS Bulletin B for x/y pole, UT1-UTC, dX, dY at daily intervals since 02. January 1973.
Celestial pole offsets (dX, dY) are related to the IAU2000A precession/nutation theory.

The format of the finals2000A.data, finals2000A.daily, and finals2000A.all files is:

Col.#    Format  Quantity
-------  ------  -------------------------------------------------------------
1-2      I2      year (to get true calendar year, add 1900 for MJD<=51543 or add 2000 for MJD>=51544)
3-4      I2      month number
5-6      I2      day of month
7        X       [blank]
8-15     F8.2    fractional Modified Julian Date (MJD UTC)
16       X       [blank]
17       A1      IERS (I) or Prediction (P) flag for Bull. A polar motion values
18       X       [blank]
19-27    F9.6    Bull. A PM-x (sec. of arc)
28-36    F9.6    error in PM-x (sec. of arc)
37       X       [blank]
38-46    F9.6    Bull. A PM-y (sec. of arc)
47-55    F9.6    error in PM-y (sec. of arc)
56-57    2X      [blanks]
58       A1      IERS (I) or Prediction (P) flag for Bull. A UT1-UTC values
59-68    F10.7   Bull. A UT1-UTC (sec. of time)
69-78    F10.7   error in UT1-UTC (sec. of time)
79       X       [blank]
80-86    F7.4    Bull. A LOD (msec. of time) -- NOT ALWAYS FILLED
87-93    F7.4    error in LOD (msec. of time) -- NOT ALWAYS FILLED
94-95    2X      [blanks]
96       A1      IERS (I) or Prediction (P) flag for Bull. A nutation values
97       X       [blank]
98-106   F9.3    Bull. A dX wrt IAU2000A Nutation (msec. of arc), Free Core Nutation NOT Removed
107-115  F9.3    error in dX (msec. of arc)
116      X       [blank]
117-125  F9.3    Bull. A dY wrt IAU2000A Nutation (msec. of arc), Free Core Nutation NOT Removed
126-134  F9.3    error in dY (msec. of arc)
135-144  F10.6   Bull. B PM-x (sec. of arc)
145-154  F10.6   Bull. B PM-y (sec. of arc)
155-165  F11.7   Bull. B UT1-UTC (sec. of time)
166-175  F10.3   Bull. B dX wrt IAU2000A Nutation (msec. of arc)
176-185  F10.3   Bull. B dY wrt IAU2000A Nutation (msec. of arc)

下面提出来一行，由于太长，将其分为几行如下所示：

73 1 2 41684.00 
I  0.120733 0.009786  0.136966 0.015902  I 0.8084178 0.0002710  0.0000 0.1916  
P    -0.766    0.199    -0.720    0.300   .143000   .137000   .8075000   -18.637    -3.667 

它们的含义依次为：

Item	含义
73 1 2	1973年1月2日
41684.00	儒列日
I	下面极移值的标签，I代表IERS，P代表预测
0.120733, 0.009786	x 向的极移值和精度，单位为 弧秒
0.136966 0.015902	y 向的极移值和精度，单位为 弧秒
I	下面 UT1-UTC 值的标签，I代表IERS，P代表预测
0.8084178 0.0002710	UT1-UTC 的值和精度，单位为秒
0.0000 0.1916	LOD(Length Of Day，日长变化)的值和精度，单位是毫秒，并不总是存在
P	下面章动值的标签，I代表IERS，P代表预测
-0.766 0.199	x 向的章动和精度，单位是毫弧秒
-0.720 0.300	y 向的章动和精度，单位是毫弧秒
	上面的都是Bull. A值，下面的是Bull. B值
.143000 .137000	x, y 向的极移值，单位为弧秒
.8075000	UT1-UTC 的值，单位为秒
-18.637 -3.667	x, y 向的章动值，单位为毫弧秒

注：Bull. A 和 Bull. B 是啥意思呢？Bull.是 Bulletin（公告，不知翻译是否准确得体）的缩写，换言之，它发两个公告，一个叫A，一个叫B。为什么要发两个公告呢？评论区大佬说：Bulletin A 是每星期发布的，发布方为美国海军在华盛顿的天文台， Bulletin B 是每个月发布的，发布方是法国巴黎的天文台。供参考~