Mikhail Vasilyevich VASILYEV,Eleonora Ivanovna YAGUDINA

(俄羅斯科學(xué)院應(yīng)用天文研究所,圣彼得堡,197110)

俄羅斯應(yīng)用天文研究所月球歷表研究現(xiàn)狀

Mikhail Vasilyevich VASILYEV,Eleonora Ivanovna YAGUDINA

(俄羅斯科學(xué)院應(yīng)用天文研究所,圣彼得堡,197110)

空間探測需要越來越高精度的觀測和歷表支持。俄羅斯科學(xué)院應(yīng)用天文研究所研發(fā)了滿足精度要求的月球歷表EPM-ERA,它與同期的法國INPOP歷表和美國DE歷表比較,精度方面有微小差距。分析了進一步提高歷表精度的方法。然而EPM-ERA僅僅依賴于法國和美國臺站提供的月球激光測距LLR數(shù)據(jù),為了克服這個限制,應(yīng)用天文研究所在推進俄羅斯西伯利亞阿爾塔伊LLR地面站的重建任務(wù),以改善LLR測站分布的幾何構(gòu)型。俄羅斯的月球-全球和月球-資源探月計劃為改進月球歷表提供了機會,參與了多次國際深空探測的俄羅斯VLBI網(wǎng)“QUASAR”也將參與這些俄羅斯月球探測計劃的精密測軌跟蹤。

月球激光測距;歷表;VLBI

0 引 言

當(dāng)前和未來的太陽系研究的發(fā)展需要越來越精密的行星歷表來支持航天計劃?；拘行菤v表也與衛(wèi)星導(dǎo)航系統(tǒng)如GPS和GLONASS等的發(fā)展緊密關(guān)聯(lián)[1]。例如,當(dāng)導(dǎo)航衛(wèi)星星歷是由星載導(dǎo)航處理器通過對衛(wèi)星運動方程數(shù)值積分來直接實現(xiàn)時, GNSS衛(wèi)星星歷自動預(yù)報的可能性使得行星星歷成為了現(xiàn)代GNSS系統(tǒng)用戶設(shè)備的基本要求。目前美國給出的DE/LE歷表系列被用于大行星、太陽和月球的標(biāo)準(zhǔn)歷表模型[2]。這些行星歷表的精度一直在改進提高著,以滿足現(xiàn)代太空探測任務(wù)的需求。遺憾的是,生成這些歷表的軟件代碼不是開源的,為此限制了科學(xué)研究者通過分析探測任務(wù)數(shù)據(jù)和歷表數(shù)據(jù)獲得科學(xué)成果的途徑,也不可能通過對這類數(shù)據(jù)的分析來提高天文常數(shù)的精度。為了克服這個障礙,歐洲空間局ESA支持了一個法國項目來構(gòu)建與DE/LE歷表精度相同的獨立的行星歷表INPOP[3]。而俄羅斯版本的月球和行星歷表EPM也由俄羅斯科學(xué)院應(yīng)用天文研究所獨立完成[4-5]。月球歷表EPM-ERA是EPM歷表的一部分[6]。用于構(gòu)建EPM歷表的軟件是ERA[7]。

1 數(shù)學(xué)模型

EPM-ERA歷表的數(shù)學(xué)模型是在G A Krasinsky的主持下構(gòu)建的。文獻(xiàn)[8-9]詳細(xì)描述了行星天體的軌道運動方程、彈性體的自轉(zhuǎn)運動特性,可以方便地用于數(shù)值積分和分析解研究,甚至用于對潮汐效應(yīng)的分析研究。為了構(gòu)建EPM-ERA 2013模型,軌道積分和天體自轉(zhuǎn)運算時,同時對太陽、月球、大行星、小行星以及TNO天體和環(huán)進行了積分。構(gòu)建EPM-ERA歷表的ERA軟件使用到的詳細(xì)數(shù)學(xué)模型在文獻(xiàn)[10]中給出。該文獻(xiàn)在J2000.0歷元太陽系質(zhì)心直角坐標(biāo)無旋轉(zhuǎn)參考框架中給出了n-體問題的軌道運動方程,并考慮了PPN參數(shù)的影響。為了構(gòu)建高精度的月球運動理論,在數(shù)學(xué)模型中還需要考慮到以下影響:

1)太陽引力勢的二階帶諧項產(chǎn)生的加速度效應(yīng);

2)地球和月球形狀帶來的擾動;

3)在地月系統(tǒng)中地球潮汐形變的滯后效應(yīng);

4)地球形狀與點質(zhì)量太陽、月亮之間的相互作用;

5)月球形狀與點質(zhì)量太陽、地球、木星、金星之間的相互作用。

2 月球自轉(zhuǎn)模型

以月球中心為參考的月固坐標(biāo)系與慣性參考系的連接可以通過3個歐拉角實現(xiàn):章動角θ,歲差角φ,以及進動旋轉(zhuǎn)角ψ。自轉(zhuǎn)角速度分量ω=(ωx, ωy,ωz)可以展開為歐拉角的導(dǎo)數(shù)組合。在旋轉(zhuǎn)坐標(biāo)系里,歐拉運動方程描述了角動量矢量的變化特征,同時需要考慮月球的特性特征。根據(jù)文獻(xiàn)[11],月球自轉(zhuǎn)的潮汐擾動是由以下因素引起的:

1)地球和太陽的引潮力導(dǎo)致月球引力勢和轉(zhuǎn)動慣量的變化;

2)月球自身向心力產(chǎn)生的擠壓導(dǎo)致月球引力勢和轉(zhuǎn)動慣量的變化;

3)月球本體滯彈性導(dǎo)致的潮汐耗散。

3個包括二階項的關(guān)于歐拉角的微分方程可以與太陽系天體的軌道運動方程同步進行數(shù)值積分運算。文獻(xiàn)[10]中給出了詳細(xì)的模型介紹。

3 EPM-ERA:實際演算

EPM-ERA歷表的全部數(shù)學(xué)模型都包含在ERA軟件包中。軟件包通過對太陽系大天體、小行星、月球等運動方程和月球自轉(zhuǎn)的同步數(shù)值積分實現(xiàn)。其中月球自轉(zhuǎn)的計算使用了文獻(xiàn)[11]中的模型。月球的引力勢考慮到4階球諧函數(shù)模型。地球的引力勢考慮到5階球諧函數(shù)模型。這參考了IERS對高軌道地球衛(wèi)星軌道積分推薦使用的要求。

帶有滯后參數(shù)τ的微分方程組的數(shù)值積分通過等步長23階Everhart方法實現(xiàn)。為了保證計算精度與LLR觀測的匹配,計算時使用了80位浮點處理器。對月球軌道和自轉(zhuǎn)動力學(xué)參數(shù)的偏導(dǎo)數(shù)通過積分可變方程組得到。特別地,這種導(dǎo)數(shù)算法可以用來計算月球的初始坐標(biāo)和速度,用以支持計算歐拉角和時間導(dǎo)數(shù)、月球引力勢系數(shù)以及地球潮汐延遲角的計算。比如月球潮汐因子K2M,其導(dǎo)數(shù)計算就是通過對方程嚴(yán)格系統(tǒng)的變參數(shù)估計值積分得到。

通過上述演算,得到了1970—2013年之間月球和行星歷表的切比雪夫內(nèi)插多項式系數(shù)文件。基于動力學(xué)和運動學(xué)理論,這些文件用于計算LLR觀測量對待估計參數(shù)的殘差和偏導(dǎo)數(shù);這里使用了IERS標(biāo)準(zhǔn)規(guī)范來獲得這類殘差。進一步,使用最小二乘法迭代解算條件方程,并且使用nσ來評估殘差水平。當(dāng)所有的參數(shù)都被擬合后,進行最后的數(shù)值積分。切比雪夫內(nèi)插多項式系數(shù)結(jié)果文件包含了1800—2200年之間的大行星、太陽、月球的位置和運動速度,以及月球的天平動,這個歷表被稱為EPM-ERA 2013。外部用戶可以使用文獻(xiàn)[12]推薦的算法流程實現(xiàn),代碼包括C、Fortran和Java版本。

4 LLR觀測數(shù)據(jù)

EPM-ERA2013歷表構(gòu)建過程中使用了1970—2013年之間18 700個LLR觀測數(shù)據(jù)用以擬合月球的軌道和自轉(zhuǎn)參數(shù)。目前,月面有5個激光反射器,分別是Apollo-11,Apollo-14,Apollo-15,以及前蘇聯(lián)與法國合作的Lunokhod1和Lunokhod2。截止到2013年為止,這些反射器的觀測數(shù)據(jù)分別是1990、1961、14 092、503和105。在整個的LLR常規(guī)觀測歷史中,主要有4個站參加了觀測,分別是美國德克薩斯州McDonald、夏威夷Haleakala、新墨西哥州Apache和法國CERGA。截止到1985年Mc Donald天文臺LLR使用了一個通用的SLR望遠(yuǎn)鏡設(shè)備,之后更新為MLRS專用LLR設(shè)備。1988年MLRS被重置在天文臺內(nèi)的一個新位置。

從2006年開始,新建設(shè)完成的Apache臺站實現(xiàn)了毫米精度月球測距觀測,使得月球軌道和動力學(xué)自轉(zhuǎn)的測量精度達(dá)到了前所未有的精度。現(xiàn)有的歷表模型DE、INPOP、和EPM-ERA都不滿足Apache臺站觀測數(shù)據(jù)的測量和解析精度的要求。

1998年之前,觀測臺站應(yīng)用戶的要求提供觀測數(shù)據(jù)。之后借助SLR/IERS服務(wù)平臺開始通過FTP方式網(wǎng)絡(luò)發(fā)布觀測數(shù)據(jù),用戶可以從http:// oca.eugemini.donnees.las_lunes以及http:// ccdisa.gsfc.nasa.gov/pub/slr獲得LLR數(shù)據(jù)。

5 擬合參數(shù)與殘差

構(gòu)建EPM-ERA2013歷表模型時總共擬合了65個參數(shù),列在表1中。由于LLR觀測量融合了地-月系統(tǒng)的運動和自轉(zhuǎn),參數(shù)之間相關(guān)性較強,無法同時從LLR觀測量中擬合這個系統(tǒng)的指向參數(shù)?？紤]到Apollo-5有著最多的觀測、位置定位精度最高,在擬合中我們固定了它的坐標(biāo)。然后使用給定的月球天平動,初步擬合一組簡單的精度較低的參數(shù)解。

LLR測站的動力學(xué)模型對地球引力場數(shù)的J2項很敏感,分析表明J2不能與月球x坐標(biāo)同時精確擬合估計,為此表1中沒有包括J2項參數(shù)。在迭代過程中,殘差按照4σ剔出。圖1是擬合殘差分布及其歷史沿革情況。

表1 擬合估計參數(shù)列表Table 1 List of estimated parameters

續(xù)表1

圖1 EPM-ERA2013歷表LLR擬合殘差Fig.1 EPM-ERA2013 fitting residuals for LLR

我們注意到與歷史模型如EPM-ERA 2010比較,EPM-ERA 2013版本的模型由于使用了新的高精度LLR數(shù)據(jù),并且軟件中優(yōu)化了數(shù)學(xué)模型,計算機也從64位浮點運算改為80位浮點運算,使得模型有了很大的改進。

6 俄羅斯LLR臺站重建

前蘇聯(lián)的克里米亞天體物理臺CRAO在1969—1983年之間成功地開展了LLR觀測實驗。并于1973年開始了常規(guī)的LLR觀測,測量精度達(dá)到了0.9 m。從1978年開始,CRAO安裝了新的激光設(shè)備,使得測量精度達(dá)到0.25 m。這個站總共得到了Apollo-15和Lunachod2反射器的1 400個LLR測距觀測數(shù)據(jù)。1983年觀測計劃中止后,俄羅斯月球歷表的觀測量主要依賴美國和法國LLR臺站提供。為改變這種狀況,包括俄羅斯科學(xué)院應(yīng)用天文研究所在內(nèi)的多家研究機構(gòu)在考慮重新啟動LLR觀測。新的臺站設(shè)置將改進當(dāng)前LLR測量網(wǎng)的幾何構(gòu)型分布。臺站目前景觀見圖2,該站主要參數(shù)如下。

1)位置:西伯利亞阿爾塔伊(50°N,82°E);

2)望遠(yuǎn)鏡口徑:3.12 m;

3)LLR觀測精度目標(biāo):3 mm;

4)新研制的激光器參數(shù)還在進一步調(diào)整中,將滿足標(biāo)稱的測量精度需求;

5)氣象條件:每年1 400晴夜小時,240夜晚適合LLR觀測。

圖2 俄羅斯阿爾塔伊天文臺Fig.2 Russian Altay station:panorama

該站址的優(yōu)勢是改進了LLR臺站網(wǎng)的經(jīng)度分布,劣勢是緯度太高,與Apache臺站比較可能會失去15%的觀測機會。新臺站的業(yè)務(wù)中,應(yīng)用天文研究所將主要負(fù)責(zé)軟件業(yè)務(wù)支持和以下工作:

1)LLR觀測計劃制定;

2)LLR觀測的歷表支持;

3)原始觀測數(shù)據(jù)的預(yù)處理以及NP計算;

4)LLR觀測量數(shù)據(jù)庫的管理;

5)EPM-ERA參數(shù)擬合改進;

6)EPM-ERA月球歷表的更新。

項目當(dāng)前的狀態(tài)是由OJC研究生產(chǎn)聯(lián)合體精密系統(tǒng)與儀器公司贏得了政府支持,牽頭負(fù)責(zé)該項目的推進。應(yīng)用天文研究所是參加機構(gòu)之一。

7 Svetloe臺站的深空探測器VLBI觀測

俄羅斯“QUASAR”VLBI網(wǎng)包括3個臺站: Svetloe、Zelenchukskaya和Badary。它們組成的基線長度分別是2 015、4 282和4 404 km。網(wǎng)絡(luò)控制和數(shù)據(jù)處理中心在圣彼得堡[13]?！癚UASAR”網(wǎng)主要開展天體測量、天體物理和大地測量觀測研究?？臻g探測器和著陸器也可以用該網(wǎng)進行高精度定位測量。由于該網(wǎng)觀測任務(wù)負(fù)擔(dān)很重,目前正在進行升級改造,在Badary和Zelenchukskaya增設(shè)了13 m天線,主要用于大地測量和天體測量觀測。預(yù)計2015年開始可以進行空間探測器常規(guī)觀測。而使用目前“QUASAR”網(wǎng)也可以開展空間探測器的VLBI觀測,觀測目的主要包括火星快車、金星快車、空間VLBI衛(wèi)星Radio Astron的國際VLBI和D-VLBI聯(lián)測,以及GLONASS、北斗等導(dǎo)航衛(wèi)星的觀測。

試驗觀測的目的是為了研發(fā)VLBI和D-VLBI技術(shù)和積累相關(guān)經(jīng)驗,用于未來對俄羅斯月球-全球和月球-資源探測計劃。已經(jīng)完成的初步試驗包括技術(shù)開發(fā)和觀測流程驗證:

1)開發(fā)用于空間衛(wèi)星VLBI觀測的軟件;

2)開發(fā)歷表服務(wù)于衛(wèi)星VLBI觀測和數(shù)據(jù)處理;

3)寬帶和點頻信號的VLBI相關(guān)處理設(shè)計研發(fā);

4)寬帶納秒和點頻皮秒精度的VLBI技術(shù)驗證。

8 結(jié) 論

俄羅斯應(yīng)用天文研究所歷表研究進展表現(xiàn)為:

1)長期維護EPM-ERA月球歷表,與同期國際上發(fā)布的月球歷表精度相當(dāng);

2)當(dāng)前EPM-ERA歷表表明,為了與LLR的測量精度需求匹配,在構(gòu)建歷表的數(shù)學(xué)模型中應(yīng)該改進月球自轉(zhuǎn)模型,達(dá)到足夠的精度,并應(yīng)用于LLR數(shù)據(jù)分析全過程;

3)最新的EPM-ERA歷表主要使用了美國和法國LLR站提供的測量數(shù)據(jù);

4)積極參與俄羅斯國內(nèi)LLR站的重建;

5)以其主導(dǎo)的“QUASAR”VLBI網(wǎng)為依托,積極推進俄羅斯空間探測器的VLBI和D-VLBI研發(fā)和驗證觀測,以改進月球歷表為目標(biāo),為觀測“嫦娥3號”/“嫦娥4號”、月球-全球和月球-資源做充分的準(zhǔn)備。

[1]Finkel'shtein A M,Gayazov I S,Gubanov V S,et al. Fundamental coordinate time supply[J].Trudy Inst.Prikl. Astron.Russ.Akad.Nauk,2005(13):1-42.

[2]Standish E M.JPL planetary and lunar ephemerides[R]. DE405/LE405Interroficememorandum,1998,no.312.F98048.

[3]Fienga A,Laskar J,Kuchynka P,et al.The INPOP10 a planetary ephemeris and its applications in fundamental physics[J].Celest.Mech.Dyn.Astron.,2011,111(3): 363-385.

[4]Krasinsky G A,Pitjeva E V,Sveshnikov M L,Chu-najeva L I.The motion of major planets from observations 17691988 and some astronomical constants[J].Celest. Mech.Dyn.Astron.,1993(55):1-23.

[5]Pitjeva E V.High precision ephemerides of planets-EPM and determination of some astronomical constants[J].SolarSyst.Res.,2005,39(3):176-186.

[6]Krasinsky G A,Vasiliev M V.ERA:knowledge base for ephemeris and dynamical stronomy[C]∥in Dynamics and Astrometry of Natural and Artificial Celestial Bodies,IAU Coll.165,Wytrzyszczak,I.M.,Lieske,J.H.,and Feldman,R.A.,Eds.,Dortrecht:Kluwer Acad.Publ., 1997:239-244.

[7]Krasinsky G A,Prokhorenko S O,Yagudina E I,New version of EPM ERA Lunar theory[R].Journees2010, Capitaine,N.,Ed.,Paris,2010:61-64.

[8]Krasinsky G A.Tidal effects in Earth and Moon rotations.1. Mathematical model[J].Trudy Inst.Prikl.Astron.Russ. Akad.Nauk,Astrometr.Geodinam.,1998(3):96-124.

[9]Krasinsky G A.Tidal effects in Earth and Moon rotations. 2.Evolution of rotating motion of the Earth and Moon orbits [J].Trudy Inst.Prikl.Astron.Russ.Akad.Nauk, Astrometr.,Geo-dynam.,Nebesn.Mekh.,1999a(4): 174-186.

[10]Vasilyev M V,Yagudina E I.Russian Lunar ephemeris EPMERA 2012[J].Solar System Research,2014,48(2):158-165.

[11]Krasinsky G A.Tidal effects in the Earth-Moon system and the Earth rotation[J].Celest.Mech.Dyn.Astron.,1999b, 75(1):40-66.

[12]Brattseva O A,Novikov F A,Pit'eva EV,et al.Access procedure to high-accuracy numerical ephemeredes of the Moon and planets[J].Trudy Inst.Prikl.Astron.Russ. Akad.Nauk,2010(21):201-204.

[13]Gayazov I,Mitryaev V,Smolentsev S,Rahimov I,Diyakov A,Shpilevski V,Pshenkin V,Rets Ya.SLR Observa-tions at the“Quasar”network stations[C]∥International Scientific-Technical Conference WPLTN-2012.Book of abstracts,41-42,2012.

Mikhail Vasilyevich VASILYEV(1958—),男,副教授,博士生導(dǎo)師,主要研究方向:衛(wèi)星導(dǎo)航,天體力學(xué),行星歷表。

通信地址:俄羅斯科學(xué)院應(yīng)用天文研究所,圣彼得堡(197110)

E-mail:mvv@ipa.nw

[責(zé)任編輯:宋宏]

Current State and Prospects of Lunar Ephemeris in IAA RAS

Mikhail Vasilyevich VASILYEV,Eleonora Ivanovna YAGUDINA
(Institute of Applied Astronomy Russian Academy of Sciences,St.Petersburg 197110,Russia)

More accurate observations and ephemerides to support the space programs are required.Currently, lunar ephemerides of adequate accuracy are being developed and maintained by IAA RAS(Institute of Applied Astronomy Russian Academy of Sciences)under the name EPM-ERA.The accuracy of EPM-ERA ephemeris is still worse a little bit than that of competitive French INPOP and American DE ephemerides.The ways to improve the quality of EPM-ERA ephemeris are presented and discussed.However,EPM-ERA ephemeris is based on only measurements kindly provided by American and French LLR observatories.To overcome this problem the project of Russian LLR station was initiated by several organizations including IAA RAS.The proposed location of this station is Siberia,Altai region that promised visible improvement of existing LLR network geometry.Perspectives of the project,current problems and progress are presented and discussed.The impact on improving ephemeris of the Moon is activity of IAA RAS connected with Russian Luna-Glob and Luna-Resource programs.VLBI network“Quasar”is prepared now to observe space vehicles and lunar landers.

lunar laser ranging;ephemerides;VLBI

P228;P691

:A

:2095-7777(2014)03-0187-05

10.15982/j.issn.2095-7777.2014.03.004

2014-07-20

2014-08-01