李 軍 趙德軍 陳永祥，2

1 西安測繪信息技術(shù)總站，西安市西影路36號，710054

2 大地測量與地球動力學(xué)國家重點實驗室，武漢市徐東大街340號，430077

地球自轉(zhuǎn)參數(shù)（ERP）包括日長變化和極移變化［1－3］。大地測量學(xué)家對ERP的預(yù)報通常是建立一個包含線性趨勢項、周年擺動項、錢德勒擺動項的最小二乘（LS）諧波模型，然后利用此模型對其進行外推預(yù)測。后來將LS殘差序列當作一個平穩(wěn)隨機時間序列，利用自回歸（AR）模型對殘差序列進行建模和推估，提出了基于LS 模型和AR模型組合的方法［4］。經(jīng)過2005～2008年多個國家參與的國際地球定向參數(shù)預(yù)報大會戰(zhàn)，發(fā)現(xiàn)LS＋AR 模型對ERP 參數(shù)的短期預(yù)報來說精度很高［3］。許多學(xué)者對LS＋AR 模型進行了研究和應(yīng)用［5－6］，并在LS＋AR 模型基礎(chǔ)上進行部分改進：文獻［7］認為LS模型中的錢德勒項具有時變特征，因而提出了基于錢德勒項時變修正的LS＋AR 模型；文獻［8］提出了加權(quán)LS＋AR 的組合模型；文獻［9］將AR 模型拓展成求和自回歸移動平均模型ARIMA，提出了LS＋ARIMA 的算法；文獻［10］將LS殘差序列看作灰色系統(tǒng)模型GM，提出LS＋GM 的算法。

地球物理學(xué)家主要分析引起ERP 變化的激發(fā)源，包括天體的引力矩和地球內(nèi)部各層的質(zhì)量重新分布等，海洋激發(fā)了ERP 部分年際變化，大氣可以解釋ERP中90%的季節(jié)性變化［11］。大氣角動量AAM 及海洋角動量OAM 的赤道方向分量X1、X2與極移之間存在著較強的相關(guān)性，軸向分量X3與日長變化有關(guān)［12］。本文分別用兩種方法推估ERP參數(shù)，對比兩者的預(yù)報效果。

1 流體激發(fā)函數(shù)

大氣和海洋是極移的主要激發(fā)源之一，極移與角動量的關(guān)系通過下面的復(fù)數(shù)微分給出［9－10］：

式中，i為虛數(shù)標記；m（t）＝x（t）－iy（t）為要計算的極移坐標的復(fù)數(shù)形式；X（t）＝X1（t）＋iX2（t）為大氣或海洋角動量的激發(fā)函數(shù)，X1（t）和X2（t）為角動量的赤道分量；σch＝［1＋i／（2Q）］2π／Tch為錢德勒擺動的復(fù)數(shù)頻率；Tch為錢德勒擺動周期，根據(jù)IERS2010規(guī)范取433d；Q是阻尼因子，本文取170。可用數(shù)值積分方法解算此微分方程：

Δt為數(shù)據(jù)采樣間隔。

2 數(shù)值實驗

2.1 數(shù)據(jù)來源及預(yù)處理

極移和日長時間序列來自IERS EOP 05C04（http：／／hpiers.obspm.fr／iers／eop／eopc04＿05／），該產(chǎn)品包含1962年至今的自轉(zhuǎn)參數(shù)，采樣間隔為1d。對于極移預(yù)測，可以直接利用該產(chǎn)品的時間序列，而日長變化主要受固體潮的影響。為保持日長時間序列的連續(xù)性，首先應(yīng)移除協(xié)調(diào)世界時UTC中的跳秒，得到UT1－TAI時間序列，再剔除周期為5 ～18.6a的固體潮帶諧項對UT1的影響δUT1，得到平滑的UT1R－TAI序列，然后基于此平滑時間序列進行預(yù)測。IERS 2010規(guī)范［1］給出了潮汐計算δUT1的詳細算法和相應(yīng)的計算軟件。

AAM 數(shù)據(jù)（包括3個分量）來自美國環(huán)境預(yù)報中心NCEP 與美國大氣研究中心NCAR 合作發(fā)布的大氣角動量再分析數(shù)據(jù)（ftp：／／ftp.aer.com／pub／anon＿collaborations／sba／），該 產(chǎn) 品 包含1948年至今的AAM，采樣間隔為0.25d，本文將其平滑成1d的采樣間隔。OAM 采用ECCO 提供的ECCO＿kf066b.oam 文件（http：／／euler.jpl.nasa.gov／sbo），時間從1993年至今，采樣間隔為1d，但是只有赤道分量（X1、X2）。

2.2 LS＋AR 模型［4－10］結(jié)果分析

采用IERS EOP 快速預(yù)報中心推薦的平均絕對誤差MAE來評估模型的預(yù)報精度：

Pj為第j期的預(yù)報值，Qj為相應(yīng)時段的觀測值，N為總預(yù)報期數(shù)。

圖1為用LS和LS＋AR 計算的極移坐標預(yù)報值的誤差曲線。極移預(yù)報模型中考慮了趨勢項、周年項和錢德勒項（周期433d）。極移預(yù)報采用的樣本時間序列為1990 ～2009年，2010、2011年的數(shù)值作為檢測值。LS建模時分別采用10、6和4a的樣本時間序列，AR 建模使用800d的數(shù)據(jù)來擬合LS 殘差序列。從圖1 中看出，對于長期預(yù)報來說，LS＋AR 模型預(yù)報精度明顯高于LS模型，且隨著預(yù)報時間的拓展精度變化相對穩(wěn)定；LS＋AR 模型中，采用10a樣本數(shù)據(jù)建模的長期預(yù)報精度高于采用4a數(shù)據(jù)的精度。

圖2為用LS＋AR 模型預(yù)報日長UT1－UTC的誤差曲線。日長預(yù)報模型中考慮了趨勢項、半年項、周年項和錢德勒項。對預(yù)處理后的UT1R－TAI序列進行LS建模時，分別采用15、10 和5a的數(shù)據(jù)，用AR模型對LS殘差序列建模時，采用500d的數(shù)據(jù)。最后對預(yù)報值再加上潮汐改正、跳秒恢復(fù)成UT1－UTC序列，并統(tǒng)計平均絕對精度。從中看出，對于長期預(yù)報，采用10a已知時間序列建模的精度高于15、5a的精度。

2.3 流體激發(fā)函數(shù)結(jié)果分析

圖3為用LS＋AR 和流體激發(fā)函數(shù)計算的極移預(yù)報值誤差趨勢。LS＋AR 采用參數(shù)與圖1相同，仍采用10a的樣本數(shù)據(jù)。并分別采用大氣角動量、海洋角動量以及大氣角動量＋海洋角動量激發(fā)函數(shù)。從圖3看出，AAM 和OAM 激發(fā)的地球極移變化在100d內(nèi)的預(yù)報誤差精度處于同一量級，但是預(yù)報跨度在100d之后，OAM 精度迅速降低；AAM＋OAM 組合模型計算的精度反而低于單獨采用AAM 或OAM 計算的精度；AAM 長期預(yù)報精度明顯高于LS＋AR模型。

圖1 LS和LS＋AR 模型預(yù)測的極移誤差Fig.1 Prediction mean absolute errors of pole motion computed from LS and LS＋AR models

圖2 LS＋AR 模型預(yù)測的UT1－UTC誤差Fig.2 Prediction mean absolute errors of UT1－UTC computed from LS＋AR model

圖3 LS＋AR 模型和流體激發(fā)函數(shù)預(yù)測的極移誤差Fig.3 Predictionmean absolute errors of pole motion computed from LS＋AR modelsand fluid excitation functions

3 結(jié) 語

1）對ERP長期預(yù)測來說，LS＋AR 組合模型的精度明顯高于LS模型；

2）對于長期預(yù)報而言，采用10a的樣本時間序列預(yù)報精度高于5a，甚至15a樣本序列精度；

3）雖然當前對ERP 變化的地球物理解釋還不太明確，但是采用大氣角動量預(yù)報的精度明顯高于LS＋AR 模型精度，也高于海洋角動量的預(yù)報精度。因此，可考慮采用大氣角動量進行確定性建模，然后對確定性模型殘差進行AR 建模，即采用AAM＋AR 的組合模型。

［1］Gerard P，Brian L.IERS Conventions 2010［EB／OL］.http：／／www.iers.org／TN36／，2010

［2］張衛(wèi)星，劉萬科，龔曉穎.EOP預(yù)報誤差對自主定軌結(jié)果影響分析［J］.大地測量與地球動力學(xué)，2011（5）：106－110（Zhang Weixing，Liu Wanke，Gong Xiaoying.Analysis of Influence of EOP Prediction Error on Autonomous Orbit Determination［J］.Journal of Geodesy and Geodynamics，2011（5）：106－110）

［3］Kalarus M，Schch H，Kosek W，et al.Achievements of Earth Orientation Parameters Prediction Comparison Campaign［J］.Journal of Geodesy，2010，84：587－596

［4］Kosek W，Kalarus M，Niedzielski T.Forecasting of the Earth Orientation Parameters－Comparison of Different Algorithms［R］.Paris，2007

［5］嚴鳳，姚宜斌.地球自轉(zhuǎn)參數(shù)短期預(yù)報方法及其實現(xiàn)［J］.大地測量與地球動力學(xué)，2012（4）：71－75（Yan Feng，Yao Yibin.Short－Term Prediction Methods and Realization of Earth Rotation Parameters［J］.Journal of Geodesy and Geodynamics，2012（4）：71－75）

［6］孫張振，徐天河.基于LS＋AR 模型的LOD 預(yù)報及其周期性影響分析［J］.大地測量與地球動力學(xué)，2013，33（4）：69－72（Sun Zhangzhen，Xu Tianhe.Prediction of Length of Day Based on LS＋AR Model and Accuracy Analysis of Different Cycle Term［J］.Journal of Geodesy and Geodynamics，2013，33（4）：69－72）

［7］張昊，王琪潔，朱建軍.對錢德勒參數(shù)進行時變修正的CLS＋AR 模型在極移預(yù)測中的應(yīng)用［J］.武漢大學(xué)學(xué)報：信息科學(xué)版，2012，37（3）：286－289（Zhang Hao，Wang Qijie，Zhu Jianjun.Application of CLS＋AR Model Polar Motion to Prediction Based on Time－Varying Parameters Correction of Chandler Wobble［J］.Geomatics and Information Science of Wuhan University，2012，37（3）：286－289）

［8］張昊，王琪潔，朱建軍.加權(quán)最小二乘法與AR 組合模型在極移預(yù)測中的應(yīng)用研究［J］.天文學(xué)進展，2011，29（3）：343－352（Zhang Hao，Wang Qijie，Zhu Jianjun.Joint Model of Weighted Least－Squares and AR in Prediction of Polar Motion［J］.Progress in Astronomy，2011，29（3）：343－352）

［9］張昊，王琪潔，朱建軍.最小二乘外推與ARIMA（p，1，0）組合模型在極移短期預(yù)報中的應(yīng)用［J］.大地測量與地球動力學(xué)，2011（3）：106－109（Zhang Hao，Wang Qijie，Zhu Jianjun.Application of Combinational Model of Least Squares Extrapolation and ARIMA（p，1，0）in Short－Term Prediction of Polar Motion［J］.Journal of Geodesy and Geodynamics，2011（3）：106－109）

［10］徐君毅.基于截距修正的GM（1，1）模型在極移實時預(yù)測中的應(yīng)用［J］.大地測量與地球動力學(xué)，2010（1）：88－91（Xu Junyi.GM（1，1）Model Based on Intercept Correction and Application to Polar Motion Prediction［J］.Journal of Geodesy and Geodynamics，2010（1）：88－91）

［11］Gross R S，F(xiàn)ukumori I，Menemenlis D，et al.Atmospheric and Oceanic Excitation of Length－of－Day Variations during 1980－2000［J］.J Geophys Res，2004，109：61－72

［12］Niedzielski T，Kosek W.Prediction of UT1－UTC，LOD and AAM χ3by Combination of Least－Squares and Multivariate Stochastic Methods［J］.Journal of Geodesy，2008，82：83－92