Alexander GUSEV,孟治國(guó),平勁松,Natalia PETROVA,Hideo HANADA

(1.中國(guó)科學(xué)院國(guó)家天文臺(tái)月球與深空探測(cè)重點(diǎn)實(shí)驗(yàn)室,北京100012;2.俄羅斯喀山聯(lián)邦大學(xué)天文系,喀山420008;3.吉林大學(xué)地質(zhì)探索科學(xué)技術(shù)學(xué)院,長(zhǎng)春130012; 4.日本國(guó)立天文臺(tái),水澤181-8588)

多層月球模型的自轉(zhuǎn)動(dòng)力學(xué)與著陸探測(cè)

Alexander GUSEV1,2,孟治國(guó)1,3,平勁松1,Natalia PETROVA2,Hideo HANADA4

(1.中國(guó)科學(xué)院國(guó)家天文臺(tái)月球與深空探測(cè)重點(diǎn)實(shí)驗(yàn)室,北京100012;2.俄羅斯喀山聯(lián)邦大學(xué)天文系,喀山420008;3.吉林大學(xué)地質(zhì)探索科學(xué)技術(shù)學(xué)院,長(zhǎng)春130012; 4.日本國(guó)立天文臺(tái),水澤181-8588)

著重于尋找月核證據(jù)并將月核模型引入月球天平動(dòng)理論中,探討如何將多種空間探測(cè)技術(shù)應(yīng)用于毫角秒精度的月球天平動(dòng)觀測(cè),進(jìn)而測(cè)定液態(tài)和固態(tài)月核參量。通過(guò)討論多層月球的月球物理參數(shù)、流體核的幾何與動(dòng)力學(xué)橢率、松散黏滯的月幔,可獲得月球相關(guān)詳細(xì)信息與參量,這些參量對(duì)評(píng)估多層結(jié)構(gòu)的月球自由天平動(dòng)很重要。物理天平動(dòng)的解析理論還可應(yīng)用于未來(lái)多種月球工作中,期望在近代月球科學(xué)研究基礎(chǔ)上能有進(jìn)一步發(fā)展。

月球;內(nèi)部結(jié)構(gòu);動(dòng)力學(xué);物理天平動(dòng);月球任務(wù)

0 引 言

月球物理天平動(dòng)的觀測(cè)及理論解析可以提供大量的關(guān)于月球內(nèi)部構(gòu)造的信息,近年來(lái)有多篇討論相關(guān)問(wèn)題和結(jié)果的文章[1-5]發(fā)表。在剛體月球物理模型與觀測(cè)結(jié)果一致度降低的情況下,研究月球自轉(zhuǎn)有助于我們深入了解其復(fù)雜的內(nèi)部構(gòu)造。

美國(guó)宇航局1998年實(shí)施的Lunar Prospector計(jì)劃,主要通過(guò)監(jiān)測(cè)月球軌道變化獲得其重力場(chǎng)感應(yīng)情況,進(jìn)而測(cè)量其潮汐影響隨時(shí)間變化的情況。重力場(chǎng)變化與高精度的月球電磁場(chǎng)和潮汐表面耗散有關(guān),測(cè)量其變化可獲得關(guān)于月球流體核的大小和組成成分等基本信息[1-2]。月震學(xué)也可提供清晰的月球內(nèi)部構(gòu)造圖像,一般是通過(guò)月面活動(dòng)實(shí)驗(yàn)室或月面天文地質(zhì)物理望遠(yuǎn)鏡測(cè)量獲得相關(guān)數(shù)據(jù)(見(jiàn)表1)。故月面著陸器,是提供月震等數(shù)據(jù)的最佳工具。月面著陸器可將激光反射鏡帶上月面,為月球激光測(cè)距(LLR)提供反射點(diǎn),為精確監(jiān)測(cè)月球軌道提供可能。同時(shí),月面著陸器還可從地面攜帶月面望遠(yuǎn)鏡,為研究月球天文地質(zhì)狀態(tài)提供可能。

表1 月球物理參量Table 1 Geophysical lunar parameters

在LLR、GRAIL、ChangE、Luna-Glob、Luna-Resource、SELENE-2等任務(wù)積累的數(shù)據(jù)基礎(chǔ)上,我們有機(jī)會(huì)實(shí)現(xiàn)月球?qū)Ш綒v表這個(gè)設(shè)想。對(duì)于月面光學(xué)及無(wú)線電定位觀測(cè)、激光/無(wú)線電單程實(shí)驗(yàn),構(gòu)建月球天文、導(dǎo)航歷表有至關(guān)重要的作用。

0.1 月核成分與構(gòu)造

通過(guò)月核動(dòng)力學(xué)數(shù)據(jù)研究,可進(jìn)行月核成分分異分析,早期月球幔殼結(jié)構(gòu)細(xì)化,核幔邊界(CMB)演化過(guò)程研究,核幔相互作用機(jī)制重構(gòu),自由及受迫月核章動(dòng)計(jì)算等研究。評(píng)估自由天平動(dòng)周期,需要考慮雙層及三層月球模型,其中包含錢(qián)德拉極移、自由核章動(dòng)、內(nèi)核極移、自由內(nèi)核章動(dòng)模式。

0.2 月核橢率

核幔相互作用(此相互作用是指,由對(duì)流運(yùn)動(dòng)及核幔層間摩擦造成的轉(zhuǎn)動(dòng)慣量交換)將導(dǎo)致月球橢率重新分布。根據(jù)克萊勞定理,非均勻旋轉(zhuǎn)體的扁率越靠近旋轉(zhuǎn)中心越小,月核的橢率應(yīng)該小于月幔。另外月球早期公轉(zhuǎn)演化中月核凍結(jié),在演化后期月幔才固體化,因而月核橢率可能大于月幔橢率。

1 月球內(nèi)部結(jié)構(gòu):基本參量

有部分基本觀點(diǎn)支持存在小尺寸月核及熱黏性的下月幔,其中月核由金屬鐵與少量的硫或氧構(gòu)成合金組成(見(jiàn)表2和圖1)。

表2 月核的月球物理參量Table 2 Geophysical parameter of lunar core

1.1 月球重力場(chǎng)結(jié)構(gòu)

通過(guò)比較Lunar Prospector(1998—1999)的高精度軌道測(cè)量、Clementine獲得的激光測(cè)高結(jié)果,以及LLR(1970—2014)數(shù)據(jù),可研究獲得月球重力場(chǎng)結(jié)構(gòu)模型,其中三者均假設(shè)存在金屬核,故其極軸方向上的轉(zhuǎn)動(dòng)慣量分量,可由下式得到。

式中:金屬鐵核的半徑為220～350 km,或?yàn)榱蜩F合金核,半徑為330～490 km[6-7]。

1.2 極軸方向轉(zhuǎn)動(dòng)慣量值

該值決定于月殼、月幔的化學(xué)、熱學(xué)及密度模型,并受月核的質(zhì)量和尺寸大小影響。另外,由圖1可知,在約1 250 km處,密度從3.8 g/cm3陡然升至7.0 g/cm3[8-10]。

1.3 月球質(zhì)心與幾何中心的非一致性

月核形成于5億年前,晚于月球本身的形成(41億年前),由此可提供證據(jù)證明月球質(zhì)心與幾何中心的非一致性與旋轉(zhuǎn)軸的慣性軸再定位的確存在[11]。

圖1 月核內(nèi)部結(jié)構(gòu)Fig.1 Interior structure of the Moon

圖2 月球物理參數(shù)分布輪廓Fig.2 Geophysical profile of the Moon

1.4 Love數(shù)估計(jì)

Love數(shù)是測(cè)量行星在外部引力作用下彈性形變的無(wú)量綱量,也是表征行星由于潮汐作用而產(chǎn)生形變的響應(yīng),k2則是其由形變產(chǎn)生的二階附加勢(shì)的相對(duì)整體形變影響的比例。LP遠(yuǎn)程跟蹤數(shù)據(jù)給出了詳細(xì)分析結(jié)果,k2=0.024 16±0.000 22,說(shuō)明部分熔融的下月幔包圍著流體核[12],可能是由約有200～300 km的月震極低速區(qū)域組成[2],而Love數(shù)的k2值本質(zhì)上限制了月核的半徑尺寸及其構(gòu)造。約2×10-5精度量級(jí)的k2值測(cè)量能力,為更高精度分辨流體核與固態(tài)核的邊界提供可能性。

1.5 月球公轉(zhuǎn)耗散的LLR分析

文獻(xiàn)[4]分析指出,存在兩種可能的耗散源,月周期的固態(tài)流及流體核的影響,月核與低黏度月幔的較差轉(zhuǎn)動(dòng),這兩種效應(yīng)均已由月球公轉(zhuǎn)數(shù)值積分方程計(jì)算。另外,考慮純鐵核的半徑為Rc=352 km,Fe-FeS共晶成分的核半徑為Rc=374 km。有幾個(gè)實(shí)驗(yàn)室仿照月核的極高溫、高壓條件,得出液態(tài)鐵核的密度為ρс=7 g/cm3,固態(tài)鐵核密度為ρc= 7.8 g/cm3,共晶成分Fe-FeS組成的核密度為ρc= 5.3 g/cm3。

2 月球的共振公轉(zhuǎn)與物理天平動(dòng)

月球赤道平面與其公轉(zhuǎn)軌道面基本一致,轉(zhuǎn)動(dòng)慣量極軸方向分量繞公轉(zhuǎn)極軸方向做受迫振蕩,且與其他兩個(gè)分量存在一個(gè)有限夾角,此夾角存在微小變化量,此即為月球緯度方向上的受迫物理天平動(dòng)。

月球繞地公轉(zhuǎn)規(guī)律由1693年提出的卡西尼三定則描述,其具體內(nèi)容為:月球自轉(zhuǎn)周期等于其公轉(zhuǎn)周期;月球赤道與黃道相交成1.57°;月球自轉(zhuǎn)軸、月球軌道平面的法線及黃道面的法線三者共面,且第三者處于前兩者之間。

月球物理天平動(dòng)是由于月球的3條主慣性軸長(zhǎng)度不等,加上橢圓軌道造成的距離改變,在地球、太陽(yáng)及其他大天體的引力影響下,發(fā)生對(duì)平均位置的偏移。物理天平動(dòng)中的經(jīng)天平動(dòng)最大幅值約為120″,比視天平動(dòng)小230倍。

3 物理天平動(dòng)可用表格

物理天平動(dòng)根據(jù)其周期來(lái)源可分為三組,第一組主要涉及三體(日-地-月)問(wèn)題,其頻率是純德勞內(nèi)參數(shù)的組合;第二組主要為月球歲差作用和行星相互作用有關(guān)的參數(shù);第三組為月球自由天平動(dòng),周期分別約為2.89 a、80.9 a、74.6 a。表3中的頻率及幅值大小均選自物理天平動(dòng)數(shù)表[4-5]。

4 多層月球

在歷史上考慮了月幔和月核的雙層模型(見(jiàn)圖3(a))中,錢(qián)德拉極移是月球突起繞軸轉(zhuǎn)動(dòng)異常的動(dòng)力學(xué)結(jié)果,對(duì)于剛體月球來(lái)說(shuō),其周期約為74.6 a,是由月球激光測(cè)月觀測(cè)得來(lái)。

表3 三層月球自由天平動(dòng)的幅值和周期[4-5]Table 3 Amplitudes and periods of free librations of the three-layer Moon[4-5]

月球自由核章動(dòng)表征月幔與流體核的較差轉(zhuǎn)動(dòng)狀況,只有核為流體或者半流體時(shí)才會(huì)存在。通過(guò)分析解的方法,李文瀟等(2014)[13]對(duì)雙層模型的自由天平動(dòng),給出了一個(gè)分析解。另外,伴隨著月球激光測(cè)月技術(shù)的發(fā)展,測(cè)量精度不斷提高,可通過(guò)軌道數(shù)據(jù),估計(jì)月球內(nèi)部構(gòu)造情況。

在雙層模型基礎(chǔ)上,考慮了存在固態(tài)內(nèi)核及流體外核的月球的三層模型,見(jiàn)圖3(b),有4種標(biāo)準(zhǔn)公轉(zhuǎn)模式,若少了內(nèi)核,模式數(shù)則減為2,沒(méi)有了流體核則減為1。4種標(biāo)準(zhǔn)模式分別為錢(qián)德拉極移、自由核章動(dòng)、自由內(nèi)核章動(dòng)、內(nèi)核極移。所有模式均因內(nèi)外核及幔的旋轉(zhuǎn)非一致性造成。

圖3 多層月球模型Fig.3 Lunar multi-layer model

5 月球動(dòng)力學(xué)領(lǐng)域新的著陸探測(cè)任務(wù)

多層月球模型自轉(zhuǎn)問(wèn)題目前依然是測(cè)月學(xué)、月球動(dòng)力學(xué)和月球物理學(xué)領(lǐng)域中的一個(gè)有待解決的關(guān)鍵問(wèn)題之一。伴隨著各類(lèi)觀測(cè)技術(shù)精度提高,這個(gè)問(wèn)題變得越發(fā)顯著。在觀測(cè)領(lǐng)域,測(cè)量月面反射鏡到地球距離的月球激光測(cè)距技術(shù)目前可以達(dá)到亞厘米和毫米精度;而最近和計(jì)劃中的GRAIL、SELENE-2中的ILOM計(jì)劃,通過(guò)測(cè)量10-12精度的月球重立場(chǎng)異?；蛘哂^測(cè)月球物理天平動(dòng)信息來(lái)實(shí)現(xiàn)對(duì)月球深部構(gòu)造的探測(cè)。

新的和計(jì)劃中的月球著陸探測(cè)為月球自轉(zhuǎn)觀測(cè)帶來(lái)了契機(jī)。這些探測(cè)計(jì)劃包括Chang'e-3/4、Luna-Glob、Luna-Resource、SELENE-2等。其中Chang'e-3/4著陸器在月球正面著陸,其攜帶的X波段無(wú)線電鎖相轉(zhuǎn)發(fā)器,為地面探測(cè)提供了往返地-月之間鏈路的高精度無(wú)線電測(cè)距機(jī)會(huì)[14]。Luna-Glob和Luna-Resource的著陸器將分別著陸在月球的南、北兩極接近永久的光照區(qū)域,并且都攜帶了超穩(wěn)定晶體振蕩器和X波段DOR信號(hào)發(fā)射設(shè)備,為月球歷表的轉(zhuǎn)動(dòng)的測(cè)量提供了機(jī)會(huì)[15]。計(jì)劃中的SELENE-2將攜帶角反射器和ILOM數(shù)字照相天頂筒,后者通過(guò)對(duì)恒星位置在CCD上的變化,直接用于測(cè)量月球自轉(zhuǎn)角速度變化[16]。

6 結(jié)束語(yǔ)

此文著重介紹了月核存在的證據(jù),并將其考慮入月球物理天平動(dòng)理論研究中,在此基礎(chǔ)上,評(píng)估了自由天平動(dòng)周期,討論了幾種不同模式的天平動(dòng)狀態(tài),并簡(jiǎn)單介紹了天平動(dòng)解析結(jié)果在獲得月核參量上的作用。根據(jù)月球內(nèi)部分層結(jié)構(gòu),討論了不同分層情況下的物理天平動(dòng),同時(shí)討論了不同的月核成分及密度對(duì)物理天平動(dòng)的影響。未來(lái)更多、更精確的觀測(cè)數(shù)據(jù),將促進(jìn)對(duì)月球,尤其是月球內(nèi)部構(gòu)造的深入研究。而月面著陸器則為以上研究提供數(shù)據(jù)來(lái)源,在更高觀測(cè)精度的需求基礎(chǔ)上,月面著陸器數(shù)量、分布以及工作需求,都還需要仔細(xì)研究探討。

致謝

作者感謝北京師范大學(xué)天文系李文瀟將本文由英文翻譯成中文。

[1]Gusev A,Kawano N,Petrova N.Investigation and exploration of the Moon[J].Actual problems of aerocosmic systems,2007(12):94-128.

[2]Gusev A,Petrova N.The Moon-2012+:Spin-orbit evolution,geophysics and selenodesy of the Moon[J]. Advances in Space Research,2008(42):289-293.

[3]Williams J M.A scheme for lunar inner core detection[J]. Geophys.Research Letters,2007(34):L03202.doi: 10.1029/2006GL028185.

[4]Williams J M,Boggs D H,Folkner W M.DE430 Lunar orbit,physical librations,and surface coordinates[P].JPL IOM335-JW,DB,WF-2013,0722-016,2013.

[5]Williams J G,Konopliv A S,Boggs D H,et al.Lunar interior properties from the GRAIL mission[J].J.Geophys. Res.2013,doi:10.1002/2013JE004559.

[6]Konopliv A S,Binder A B,Hood L L,et al.Improved gravity field of the Moon from lunar prospector[J].Science, 1998(281):1476-1478.

[7]Konopliv A S,Park R S,Yuan D N,et al.High-resolution lunar gravity fields from the GRAIL[J].Primary and Extended Missions,Geophys.Res.Lett.,2014(41): 1452-1458.

[8]Dickey J O,Bender P L,Faller J E,et al.Lunar laser ranging:a continuing legacy of the Apollo program[J]. Science,1994(265):482-490.

[9]Kuskov O L,Kronrad V A.Constitution of the Moon,5, constraints on composition,density,temperature,and radius of a core[J].PEPI,1998(107):285-306.

[10]Kuskov O L,Kronrad V A.Core sizes and internal structure of Earth's and Jupiter's satellites[J].Icarus,2001(151): 204-227.

[11]Runcorn S K.The formation of the lunar core[J]. Geochimica et Cosmochimica Acta 60,1996(7):1205-1208.

[12]Yoder D F,Konopliv A S,Yuan D N,et al.Fluid core size of Mars from detection of the solar tide[J].Science,2003 (300):299-303.

[13]李文瀟,Gusev A,平勁松,等.雙層結(jié)構(gòu)月球的自由天平動(dòng)與受迫天平動(dòng)[J].深空探測(cè)學(xué)報(bào),2014,1(3):205-209.[Li W X,Gusev A,Ping J S,et al.Free and forced librations of tww-layer moon[J].Journal of Deep Space Exploration, 2014,1(3):205-209.]

[14]平勁松,王明遠(yuǎn),張素君,等.嫦娥系列探月衛(wèi)星無(wú)線電科學(xué)實(shí)驗(yàn)簡(jiǎn)介[J].深空探測(cè)學(xué)報(bào),2014,1(3):192-199.[Ping J S,Wang M Y,Zhang S J,et al.Radio science experiments in Chang'E series missions[J].Journal of Deep Space Exploration,2014,1(3):192-199.]

[15]Kosov A S,Gotlib V M,Korogod V V,等.俄月球-全球等探測(cè)任務(wù)中的無(wú)線電科學(xué)試驗(yàn)[J].深空探測(cè)學(xué)報(bào),2014,1 (3):181-186.[Kosov A S,Gotlib V M,Korogod V V,et al.Radio science experiments in russian Luna-Recourse and Luna-Glob projects[J].Journal of Deep Space Exploration, 2014,1(3):181-186.]

[16]Hanada H,Tsuruta S,Asari K,等.用于未來(lái)月球探測(cè)的ILOM技術(shù)現(xiàn)狀[J].深空探測(cè)學(xué)報(bào),2014,1(3):200-204. [Hanada H,Tsuruta S,Asari K,et al.Present status of ILOM for future lunar missions[J].Journal of Deep Space Exploration,2014,1(3):200-204.]

Alexander GUSEV(1954—),男,教授,主要研究方向:月球物理天平動(dòng)理論解析,月球內(nèi)部結(jié)構(gòu)分析,月球動(dòng)力學(xué)研究。

通信地址:中國(guó)科學(xué)院國(guó)家天文臺(tái)(100012)

電話(huà):(010)64807839

E-mail:Alexander.Gusev@mail.ru

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

Spin-Orbit Dynamics of the Multi-Layer Moon for Landing Missions

Alexander GUSEV1,2,MENG Zhiguo1,3,PING Jinsong1,Natalia PETROVA2,Hideo HANADA4
(1.Key Laboratory of Lunar and Deep Space Exploration,National Astronomical Observatories China,China Academy of Sciences,Beijing 100012,China;2.Kazan Federal University,Kazan 420008,Russia; 3.College of Geoexploration Science and Technology,Jilin University,Changchun 130012,China; 4.National Astronomical Observatory of Japan,Mizusawa 181-8588,Japan)

The emphasis will be put on the evidences of lunar core existence and on the necessity to take this fact into account in the lunar librations theory.Our goal is to show how the millisecond precision observations of lunar physical librations in the projects.We discuss geophysical parameters,geometrical and dynamic ellipticity of liquid core and viscose-elastic mantel of the multilayered Moon.The research results include the survey of internal structure of the Moon,tabulated values of geophysical parameters and geophysical profile of the Moon,including liquid lunar core,values of geometrical compression on the Moon's surface,densities,as well as graphic representations of the received data.The new prospects for horizons of research and development of the Moon for the nearest decade will be opened.

Moon;interior;dynamics;physical libration;lunar missions

P691

:A

:2095-7777(2014)03-0175-06

10.15982/j.issn.2095-7777.2014.03.002

2014-07-25

2014-08-19

中科院外籍專(zhuān)家項(xiàng)目(2013T2J0016);中科院俄烏白國(guó)際合作項(xiàng)目(2013-2014);國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)項(xiàng)目(2015CB857101);俄日科技合作項(xiàng)目RFBR(N-14-02-91212)