月球重力場研究的論述

楊宏偉, 趙文津, 吳珍漢

中國地質(zhì)科學(xué)院, 北京 100037

月球重力場研究的論述

楊宏偉, 趙文津, 吳珍漢

中國地質(zhì)科學(xué)院, 北京 100037

月球重力場探測是了解月球內(nèi)部結(jié)構(gòu)構(gòu)造最有效的手段之一, 也是未來登陸點(diǎn)選擇的重要依據(jù)。在眾多方法中, 衛(wèi)星重力探測是進(jìn)行全月球覆蓋以及獲得月球深部信息的重要方法。本文敘述了當(dāng)前月球重力探測歷史和重力模型發(fā)展過程, 并對這些模型進(jìn)行比較分析。在了解前沿研究方法的基礎(chǔ)上討論了月球重力位場解算的基本原理、高精度月球重力場獲得的方法和各種月球重力歸算方法, 以及當(dāng)前國際上利用月球重力場研究其內(nèi)部結(jié)構(gòu)構(gòu)造的熱點(diǎn)問題。最后對月球的典型重力異常特征進(jìn)行了綜合性的分析。

月球重力場探測; 月球重力場模型; 典型重力異常特征

了解一個行星的內(nèi)部結(jié)構(gòu)構(gòu)造可以為該行星形成和演化的研究提供重要依據(jù)。由于月球遠(yuǎn)離地球?qū)е氯祟惖窃绿綔y受到局限, 人類只能在月球登陸點(diǎn)附近有限的區(qū)域內(nèi)獲得月球地質(zhì)信息, 或者只能通過衛(wèi)星遙感等手段對月球表層物質(zhì)進(jìn)行研究。在眾多的方法中, 只有衛(wèi)星重力測量能夠獲得全月范圍的甚至包括人類未能登陸地區(qū)的重要的月球深部信息。因此, 通過衛(wèi)星重力測量研究月球的重力場,是了解月球內(nèi)部結(jié)構(gòu)以及認(rèn)識月球形成和演化的重要手段, 也是當(dāng)前最為可行的方法。因此, 月球全球高精度重力場的獲得和研究是當(dāng)前認(rèn)識月球內(nèi)部信息的關(guān)鍵所在。再者, 高精度的月球重力場也是保證未來登月艙精準(zhǔn)、安全著陸的重要數(shù)據(jù)。因此, 如何獲得高精度月球重力場以及對重力場的研究是月球探測和月球研究的重要內(nèi)容之一。

1 月球重力場探測的歷史和月球重力模型的發(fā)展

1.1 月球重力場探測的歷史任務(wù)

回顧月球衛(wèi)星重力探測歷史。自1966年起, 蘇聯(lián)的Luna 10號探測器成為首顆月球軌道衛(wèi)星并率先建立了月球重力場。盡管只有56天, 219個探測數(shù)據(jù), 350 km(近月點(diǎn))至 1000 km(遠(yuǎn)月點(diǎn))的高度,但卻第一次建立了月球重力場(Akim, 1966)。1966年至1967年間, Lunar Orbiter系列軌道探測器總共獲得了32600個用來計算月球重力場的觀測數(shù)據(jù)。Lunar Orbiter I-III軌道探測器系列在月球正面的赤道范圍內(nèi)近月點(diǎn)飛行高度達(dá)到40~1800 km, 取得了月球正面南緯20°至北緯30°范圍內(nèi)較Luna 10更高精度的月球重力場的模型。其中, Lunar Orbiter 4號和5號探測器還首次獲得了一些月球兩極區(qū)域的軌道數(shù)據(jù)(Michael et al., 1972)(如圖1)。圖1顯示了不同高度軌道數(shù)據(jù)在月球表面的覆蓋情況。1971年Apollo 15號和1972年的Apollo 16號的子衛(wèi)星首次將探測器降低到100 km高度的近圓軌道, 雖然只覆蓋了南北緯 300附近的范圍, 但是卻大大改進(jìn)了月球這一區(qū)域的重力場信息(Bills et al., 1980)。盡管如此, 以上探測任務(wù)的軌道覆蓋率以及數(shù)據(jù)采樣率都比較低。至1994年, Clementine計劃首次獲了月球激光高度信息, 這對于改進(jìn)月球低階重力場模型尤為重要(Zuber et al., 1994)。1998年Lunar Prospector進(jìn)行了100 km低圓極軌的探測, 此次探測是有史以來覆蓋率和采樣率最高的一次, 獲得了 3648853個數(shù)據(jù)。盡管如此, 由于月球的背面總是背對地球, 以上所提到的所有重力探測任務(wù)都沒能直接探測到月球背面的重力場信息。這部分未能覆蓋的區(qū)域占月球總面積的33%(如圖2)。這一關(guān)鍵性的難題導(dǎo)致了即使利用 LP探測數(shù)據(jù)建立的在當(dāng)時最高精度的月球重力場數(shù)據(jù)在月球背面的重力數(shù)據(jù)及其解釋也是不可靠的。其中 LP150Q模型的月球背面重力場的誤差達(dá)到 200 mGal, 月球正面的誤差為 30 mGal,而且異常的位置相對真實(shí)位置還存在水平偏移(Wieczorek et al., 2006)。2007 年, 日本Kaguya(SELENE)探測器利用一顆中繼衛(wèi)星和主衛(wèi)星以及 4-way Doppler雷達(dá)探測技術(shù)首次探測到了這塊一直未被覆蓋的月球背面的重力場信息。這是有史以來利用直接的探測數(shù)據(jù)首次建立了全月球的重力場信息(Namiki et al., 2009), 如圖3所示。至此,全月球重力場才能夠?yàn)樵虑騼?nèi)部結(jié)構(gòu)特征的解釋提供可靠依據(jù)。

1.2 月球重力場球諧模型的發(fā)展

利用衛(wèi)星重力測量恢復(fù)月球重力場是當(dāng)前研究月球重力的重要方法。主要影響重力場精度的因素為: 軌道傾角決定覆蓋范圍, 軌道高度決定重力場精度, 采樣間隔決定月球重力場在某高度上實(shí)際恢復(fù)重力場階次的高低以及軌道覆蓋率也決定重力場位系數(shù)的精度?？v觀月球重力場模型的發(fā)展過程,都是圍繞這些方面進(jìn)行不斷改進(jìn)的。

圖1 Lunar Orbiter系列探測器不同高度觀測值的個數(shù)(左: 月球正面; 右: 月球背面)(據(jù)Michael et al., 1972)Fig. 1 Coverage of observations in different altitudes from Lunar Orbiter series (left: Nearside, right: Farside)(after Michael et al., 1972)

圖2 Lunar Prospector計劃建立的未加約束的50階次月球重力場, 紫色區(qū)域?yàn)楸趁嫖锤采w區(qū)域, 占月球總面積的33%(Konopliv et al., 2001)Fig. 2 Unconstrained gravity field with degree and order 50 from Lunar Prospector. No direct observations were measured in the purple area, which equals to 33% of total area of the moon. (after Konopliv et al, 2001)

圖3 Kaguya(SELENE)首次利用直接探測數(shù)據(jù)建立全月自由空氣重力異常(Namiki et al., 2009)Fig. 3 The 1st free gravity anomalies field of the moon in global with direct tracking data from Kaguya (SELENE)(after Namiki et al., 2009)

1966年發(fā)射的Lunar Orbiter I、II和1967年發(fā)射的Lunar Orbiter III的傾角在 10o～20o內(nèi)變化, 而1967年后半年發(fā)射的Lunar Orbiter IV、V的全球覆蓋有效數(shù)據(jù)只有不到 40000個數(shù)據(jù)點(diǎn), 因此能夠解算的球諧階次較小, 且軌道覆蓋區(qū)域主要集中在月球正面的南緯 20o～北緯 30o范圍。1968年,Lorell和 Sjogren(1968)利用 Lunar Orbiters的5個探測任務(wù)建立了月球重力球諧場 8×4模型(階數(shù)為 8,級數(shù)為 4)(Lorell et al., 1968)。之后 Michael和Blackshear(Michael et al., 1972)又獲得了 13×13(階數(shù)和級數(shù)都為 13)的球諧模型, 此模型在月球南緯200至北緯300范圍較前一模型有所提高。Apollo15、16號子衛(wèi)星的軌道高度降低到了100 km, 但傾角為29°, 因此對比前面的重力模型不僅擴(kuò)大了重力場的覆蓋范圍而且在此區(qū)域的精度也有所提高。至1980年, Bills和Ferrari(1977, 1980)利用Lunar Orbiter共5次探測數(shù)據(jù)和Apollo軌道數(shù)據(jù)以及激光測距數(shù)據(jù)建立了 16階次模型, 這一模型一直沿用了 13年。隨著高性能計算機(jī)技術(shù)的提高, Konopliv等人(1993),將球諧模型提升到了60階次(Lun60d模型)。高精度的激光數(shù)據(jù)有助于改進(jìn)月球球諧重力場的低階系數(shù)。1994年Clementine任務(wù)獲得的大量激光測距數(shù)據(jù)改進(jìn)了低階重力場系數(shù)并將階數(shù)再一次提升到了70階次(GLGM-2模型)(Lemoine et al., 1997)。而GLGM-2模型在高緯度地區(qū)的精度比赤道附近的重力場精度要低。1998年, Konopliv等人(1999)利用100 km高度和30 km高度的Lunar Prospector數(shù)據(jù)逐次建立了 75階次月球位場模型(LP75D, LP75G),100階次月球位場模型(LP100J, LP100K), 150階次月球位場模型(LP150Q)以及165階次月球位場模型(LP65P)。其中, 150階模型是利用新的算法得到的。但是上述這些模型都因缺少背面直接的探測數(shù)據(jù),因此背面的精度只能達(dá)到 200 mGal(正面的數(shù)據(jù)能夠達(dá)到 30 mGal), 這對于地球物理解釋來說是非常不可靠的(Wieczorek et al., 2006)。2009年 Namiki等人(2009)在Kaguya任務(wù)中利用雙星技術(shù)直接探測得到了月球背面的重力數(shù)據(jù)解決了缺少背面數(shù)據(jù)的難題, 首次建立了全月球的球諧場模型——SGM90d, 即90階次模型。

圖4、5顯示出美國Lunar Prospetor任務(wù)解算的LP165P與日本SELENE解算的SGM90d這兩個模型在月球背面以及背面的Apollo坑(36.1°S, 151.8°W)的重力場的對比。雖然在沒有背面探測數(shù)據(jù)時利用Kaula約束同樣能夠解算高階次月球重力場, 但是在月球背面的精度上較直接探測建立的結(jié)果差距較大。圖4、5便顯示出在沒有背面直接探測數(shù)據(jù)時比有直接探測數(shù)據(jù)時的精度要低。

由此看出, 全球范圍、多次覆蓋、高密度采樣率、低軌道的探測數(shù)據(jù)是獲得高精度重力模型的基礎(chǔ)。

圖4 LP150Q(150th, a)和SGM90d(90th,b)月球背面自由空氣重力場的對比(據(jù)JAXA, 2010)Fig. 4 Comparison of free air gravity anomalies on lunar farside between LP165P (165th, a)and SGM90d (90th, b) (Provided by JAXA, 2010)

圖5 月球背面Apollo坑(36.1°S / 151.8°W)自由空氣重力異常(a圖來自LP150Q模型, b圖來自SGM90d模型)(據(jù)JAXA-Kaguya image gallery, 2008)Fig. 5 Free gravity anomaly in Apollo crater (36.1°S, 151.8°W) on lunar far side(They are retrieved respectively from LP 150Q (a) and SGM90d (b)) (Provided by JAXA-Kaguya image gallery, 2008)

表1 比較重要的月球重力場球諧模型Table 1 Some of fundamental lunar gravity coefficient models

2 高精度月球球諧重力場模型的解算

2.1 重力場解算原理

美國和日本的重力數(shù)據(jù)是以三種方式給出月球重力場的。首先是上面提到的SGM90d以及LP150Q等月球重力場球諧模型; 第二種是給出月球網(wǎng)格點(diǎn)上的相應(yīng)重力值。第三種是給出衛(wèi)星軌道數(shù)據(jù), 如衛(wèi)星位置或速度數(shù)據(jù)等。高精度月球重力場模型的解算就是利用衛(wèi)星軌道數(shù)據(jù)建立的。

以下將介紹利用衛(wèi)星軌道數(shù)據(jù)計算得到重力場球諧模型的基本原理和計算流程(圖6), 這是解算高精度重力場球諧模型的基礎(chǔ)。

圖6 利用軌道跟蹤數(shù)據(jù)計算球諧重力場模型流程圖Fig. 6 Flow chart for solution of gravity coefficient model from orbit tracking data

首先利用 Doppler和測軌數(shù)據(jù)計算出月球坐標(biāo)系中衛(wèi)星的軌道坐標(biāo), 將軌道坐標(biāo)利用 Taylor級數(shù)展開成與六軌道根數(shù)表示的級數(shù)形式。而Kelper定律又將軌道根數(shù)和重力場球諧模型的擾動重力位建立了聯(lián)系, 同時擾動重力位又可利用Laplace方程的系數(shù)展開表示, 由此便建立了衛(wèi)星坐標(biāo)與球諧系數(shù)之間的聯(lián)系。為了能夠獲得更高精度的重力場球諧系數(shù)模型主要從以下幾個方面考慮: 提高定軌精度,如利用VLBI技術(shù)等; 降低軌道高度; 降低采樣間隔;增加軌道覆蓋率。在最后的方程解算當(dāng)中, Konopliv等人曾提出利用一個完整矩陣來解算可以有效提高解算精度(Konopliv et al., 2001)。

重力位場球諧模型表達(dá)式:

上面的公式就是由 Laplace方程推導(dǎo)得出并進(jìn)行規(guī)一化處理和參數(shù)提取得到的(Hofmann- Wellenhof, 2005)。利用上面的重力位場球諧系數(shù)模型和Matlab軟件本人編寫了部分的重力位場計算代碼,并得到全月重力位場圖(圖7)。這個模型是今后解算各種月球重力異常場的基礎(chǔ)。

2.2 高精度月球重力場的解算

要獲得高精度月球重力場需要考慮下列因素:軌道覆蓋率、軌道采集率、軌道高度、觀測數(shù)據(jù)的精度等(陳俊勇等, 2005)。

(1)軌道覆蓋率

圖7 SGM90d (90階次)全月球重力位場圖, 25 km間隔(0.833°)(據(jù)JAXA, 2010)Fig. 7 Lunar gravity coefficient field in global from SGM90d (90th), spacing of 25 kilometers (equal to 0.833°)(Provided by JAXA, 2010)

一般來說, 當(dāng)前的月球軌道探測器都是圓極軌,因此能夠保證全球覆蓋以及全球等精度探測。即便如此, 如果能夠增加不同衛(wèi)星數(shù)據(jù)聯(lián)合求解就可以大幅實(shí)現(xiàn)全球的多次覆蓋, 使探測精度增加。除此以外, 由于月球背面總是背向地球, 就需要利用中繼衛(wèi)星來探測背面數(shù)據(jù)。日本于 2007年發(fā)射的SELENE衛(wèi)星利用一顆中繼衛(wèi)星首次直接探測到了月球背面重力場, 建立了月球全球重力場模型(Namiki et al., 2009)。

(2)軌道采集率和軌道高度

結(jié)合Hofmann-Wellenhof等人(2005)和陳俊勇等人(2005)提出的階次與精度的關(guān)系, 可用圖 8來表示,

圖8表示與緯線平行的分割是與解算的月球重力球諧場中的階數(shù)(degree)有關(guān)的, 而與精線平行的分割則與月球重力球諧場中的級數(shù)(order)有關(guān)。

其中重力場精度與階次的關(guān)系為:

這里D就是解算出來的精度,r是月球半徑, N是對應(yīng)階次;

當(dāng)近似認(rèn)為軌道為圓軌道時, 利用軌道高度就可以解算出對應(yīng)高度處的速度值, 隨之便可以解算出對應(yīng)精度(對應(yīng)星下點(diǎn)間距)的時間采樣間隔:

因此, 在某一軌道高度上不同的時間采樣間隔決定了探測精度, 見表2,

由此看來, 階次最終決定了重力場的精度, 時間采樣間隔決定在某一高度上所能達(dá)到探測精度的最小采樣時間間隔。

根據(jù)上面的理論并結(jié)合實(shí)際數(shù)據(jù)處理還發(fā)現(xiàn)了另外一些規(guī)律。日本JAXA公布的軌道數(shù)據(jù)的采樣間隔是60 s, 軌道高度為100 km, 結(jié)合上面的理論了解到此數(shù)據(jù)能夠解算出 60階次的月球重力場模型。實(shí)際解算結(jié)果如下:

如圖9~13所示, 隨著解算階次的升高, 精度逐漸增加, 但是從55階開始出現(xiàn)較大的噪聲, 高階次的結(jié)果的有效性所有降低。因此解算結(jié)果首先說明上面的理論是可靠的, 但是解算出的最后幾個階次會產(chǎn)生很大噪聲(圖11~13紅框部分)。

圖8 階次與精度的關(guān)系(Hoffmann-wellenhof et al.,2005)Fig. 8 Relationship between degree and accuracy(after Hoffmann-wellenhof et al., 2005)

(3)觀測數(shù)據(jù)的精度

當(dāng)前國際上普遍采用 Doppler雷達(dá)探測技術(shù)。當(dāng)前國際上最好的跟蹤精度為 0.1～0.3 mm/s(陳俊勇等, 2005)。為增加探測精度VLBI是未來增加精度的重要探測技術(shù)。

3 月球重力場的改正

為了能夠利用重力位場模型來認(rèn)識月球內(nèi)部結(jié)構(gòu), 需要對重力位場模型進(jìn)行相應(yīng)重力改正, 從而建立月球自由空氣重力異常場、月球布格重力異常場以及月球均衡重力異常場。

自由空氣重力異?？赏ㄟ^將月球重力值歸算到某一高度處, 并減去正常重力值。布格重力改正主要涉及地形和中間層改正。月球的布格改正的處理方法眾多, 而只有利用合適的計算方法得到的布格重力異常才能準(zhǔn)確了解地下物質(zhì)分布的信息?，F(xiàn)階段較重要的月球布格校正方法主要分為3種:

Prism作為區(qū)域布格校正是一種不含行星曲率的校正; 如需全球布格校正, 還需要引入曲率校正的改進(jìn);

Wild-Pfeiffer提出的Tessoid網(wǎng)格法可進(jìn)行含有曲率校正的全球布格校正, 它的計算效率比 Prism曲率校正改進(jìn)算法的效率要高一個數(shù)量級(Heck et al., 2007; Wild-Pfeiffer, 2008);

表2 軌道高度、階次與時間采樣間隔的關(guān)系Table 2 Relationship of latitude of satellite, degree, and interval of samplings

圖9 50階次月球自由空氣重力異常場(據(jù)JAXA, 2010)Fig. 9 50th degree free gravity anomaly field of the moon (Provide by JAXA, 2010)

圖10 53階次月球自由空氣重力異常場(據(jù)JAXA, 2010)Fig. 10 53th degree free gravity anomaly field of the moon (Provide by JAXA, 2010)

Bill和Ferrari(Bill et al., 1977)提出另一個計算思路:

當(dāng)計算一個行星的布格重力場時, 需要給定該行星的大地水準(zhǔn)面的位置信息。月球大地水準(zhǔn)面是月球布格重力場計算的參考面。對月球來說, Zuber等人(1994)把月球大地水準(zhǔn)面定為扁率為1/3234.93,半徑為1738 km的橢球體; 而Konopliv等人(2001)則將大地水準(zhǔn)面設(shè)定為1738 km的一個正球體。

利用JAXA提供的軌道數(shù)據(jù)計算了月球自由空氣重力異常圖, 如圖14所示。結(jié)合月形數(shù)據(jù)建立月球布格重力異常圖, 如圖15所示。

4 月球典型重力異常的研究

為了能夠正確地分析月球重力場特征, 需要利用已經(jīng)建立的月球地質(zhì)信息將不同巖性特征或顯示出不同構(gòu)造特征的地區(qū)加以分類。根據(jù)Apollo月巖樣品以及遙感衛(wèi)星得到的月球表面礦物和元素特征(圖16)以及由激光高度計探測得到的月球地形特征,Jolliff等人(2000)將月球分為三個不同區(qū)域。

圖11 55階次月球自由空氣重力異常場(據(jù)JAXA, 2010)Fig. 11 55th degree free gravity anomaly field of the moon (Provide by JAXA, 2010)

圖12 58階次月球自由空氣重力異常場(據(jù)JAXA, 2010)Fig. 12 58th degree free gravity anomaly field of the moon (Provide by JAXA, 2010)

如圖16、17所示, 三個主要月球地質(zhì)構(gòu)造區(qū)域包括月球正面的風(fēng)暴洋 KREEP地區(qū)(PKT: Procellarum KREEP Terrane), 背面的長石質(zhì)月球高地(FHT: Feldspathic Highland Terrane)以及南極艾肯盆地(SPA: South Pole Aitken Basin)。這些典型的地質(zhì)構(gòu)造區(qū)塊不僅與月形、礦物和元素分布有關(guān)還和月球重力異常有關(guān)。如圖 16所示, 在 PKT地區(qū) FeO和Th的含量較高; 在SPA地區(qū)其FeO和Th的含量也很高, 但Th的含量卻沒有PKT地區(qū)高; 而在FHT地區(qū)二者都很低。如圖17所示, 這些典型地質(zhì)特征還與重力異常區(qū)域有著重要的關(guān)系。圖 17表明在PKT地區(qū)其地形低矮, 但有很高的布格重力異常;在FHT地區(qū)是地形很高的月球高地, 而布格重力異常卻很低; 在SPA地區(qū)其特征又表現(xiàn)出高布格重力異常和低地形。

在月球重力場中另一個最顯著的特點(diǎn)就是高重力異常的盆地。Muller和 Sjogren稱其為 Mascon(Mass Concentration的簡稱)盆地(Muller et al.,1968)。在此需要注意的是Mascon的地質(zhì)含義并不與它的名字一樣是與“質(zhì)量瘤”有關(guān), 只不過是與高重力異常有關(guān)。至于其是否是由于高密度體所引起的并沒有得到證實(shí), 而且有資料表明這種異常與殼幔界面的起伏有關(guān)(后文有敘述)。因此, 下文即使使用Mascon, 也只是作為名稱引用。

圖13 60階次月球自由空氣重力異常場(據(jù)JAXA, 2010)Fig. 13 60th degree free gravity anomaly field of the moon (Provide by JAXA, 2010)

圖14 SGM90d (90階次)建立的月球自由空氣重力異常場(據(jù)JAXA, 2010)Fig. 14 Free air gravity anomaly field of the moon from SGM90d (90th) (Provide by JAXA, 2010)

圖15 利用重力位場(GLGM-1)和地形模型(GLTM-1)計算月球布格重力異常圖(Zuber et al., 1994)Fig. 15 Bouguer anomaly field of the moon from gravity coefficient model (GLGM-1) and topography model(GLTM-1) (Zuber et al., 1994)

Mascon主要表現(xiàn)為具有較高的正的重力異常和較低的地形特征(如圖18、19)。

利用 Kaguya衛(wèi)星得到的新的高精度重力球諧模型新發(fā)現(xiàn)了許多更小的Mascon盆地。至此, 月球上主要的Mascon盆地, 見表3。

圖16 月表FeO含量分布圖(a)和月表Th含量分布圖(b)(Jolliff et al., 2000)Fig. 16 Distribution of FeO (a) and Th (b) on lunar surface (Jolliff et al., 2000)

圖17 月球三個重點(diǎn)地區(qū)地形圖(a)與對應(yīng)的布格重力異常圖(b)(據(jù)JAXA, 2010)Fig. 17 Toporaphy(a) and Bouguer gravity anomaly(b) of the three typical areas on the moon(Provide by JAXA, 2010)

人們根據(jù) Mascon是否被月海玄武巖所覆蓋的原則, 將Mascon主要分為Mare Mascon Basin(月海Mascon盆地)和 Non-mare Mascon Basin(非月海Mascon盆地)。Konopliv等人(Konopliv et al., 2001;Konopliv et al., 1998)指出月海Mascon盆地包括月球正面的五個主要Mascon盆地, 即Imbrium(雨海),Serenitatis(澄海), Crisium(危海), Humorum(濕海),Smythii(史密斯海)。非月海 Mascon盆地包括有Schiller-Zucchius(席勒坑), Bailly(貝利坑),Schr?dinger(施羅丁格爾坑), Lorenta(勞倫茨坑),Hertzsprung(赫茲斯朋坑)。隨著新的更小的 Mascon的發(fā)現(xiàn), 這兩種類型的 Mascon盆地也在不斷更新中。

圖18 東方?！^低的地形和較高的布格重力異常特征(據(jù)JAXA, 2010)Fig. 18 Mare Orientale – low elevation and high Bouguer gravity anomaly(Provide by JAXA, 2010)

圖19 Mascons在月球表面分布圖(據(jù)JAXA, 2010)Fig. 19 Distribution of Mascons on lunar surface(Provide by JAXA, 2010)

表3 月球主要Mascon盆地Table 3 Main Mascon basins in the moon

許多人都對此特殊構(gòu)造進(jìn)行了重力反演和解釋。有人解釋這些異常特征是由未被補(bǔ)償?shù)男鋷r和殼幔界面起伏共同引起的(Wise et al., 1970)。也有人認(rèn)為是由月幔超均衡隆起引起的(Konopliv et al.,1998)。也有人認(rèn)為月海玄武巖是產(chǎn)生這種高重力異常的影響因素。

5 小結(jié)和討論

(1)月球重力場研究是了解月球內(nèi)部信息的重要手段, 也是未來登陸點(diǎn)選擇和登陸艙精準(zhǔn)著陸的重要保證;

(2)高精度重力場解算是當(dāng)前月球重力場研究的重點(diǎn), 急需低軌道、高采樣率和高覆蓋率的探測數(shù)據(jù); 盡管當(dāng)前國際上已經(jīng)有大量的數(shù)據(jù)公布, 但是數(shù)據(jù)仍然是有限的。為了能夠達(dá)到地球物理解釋還急需大量高精度數(shù)據(jù)的支持, 其中預(yù)計在2011年9月發(fā)射的 GRAIL衛(wèi)星其主要目的就是要獲得更高精度的月球重力場數(shù)據(jù);

(3)月球重力場的解釋還需要更多其他資料的結(jié)合, 如Apollo地震數(shù)據(jù)以及月表礦物和元素的分布特征和地質(zhì)信息, 并對重點(diǎn)地區(qū)進(jìn)行全面地分析和總結(jié)。

陳俊勇, 寧津生, 章傳銀, 羅佳. 2005. 在嫦娥一號探月工程中求定月球重力場[J]. 地球物理學(xué)報, 48(2): 275-281.

AKIM E L. 1966. Determination of the Gravitational Field of the Moon from the Motion of the Artificial Lunar Satellite ‘Lunar-10’[J]. Doklady Akademii Nauk SSSR, 170: 799-802.

BILLS B G, FERRARI A J. 1977. A Lunar Density Model Consistent with Topographic, Gravitational, Librational, and Seismic Data[J]. Journal of Geophysical Research, 82: 1306-1314.

BILLS B G, FERRARI A J. 1980. A Harmonic Analysis of Lunar Gravity[J]. Journal of Geophysical Research, 85: 1013-1025.

CHEN Jun-yong, NING Jin-sheng, ZHANG Chuan-yin, LUO Jia.2005. On the Determination of Lunar Gravity Field in the Chinese First Lunar Prospector Mission[J]. Chinese Journal of Geophysics, 48(2): 275-281(in Chinese).

HECK B, SEITZ K. 2007. A Comparison of the Tesseroid, Prism and Point-mass Approaches for Mass Reductions in Gravity Field Modeling[J]. Journal of Geodesy, 81: 121-136.

HOFMANN-WELLENHOF B, MORITZ H. 2005. Physical Geodesy[M]. US: Springer Wien New York.

Japan Aerospace Exploration Agency. 2008. Japan: Tokyo. Kaguya Image Gallery [EB/OL]. [2009-02-25]. http://wms.selene.jaxa.jp/selene_viewer/en/observation_mission/-rsat/-rsat_001.html.

Japan Aerospace Exploration Agency. 2010. Japan: Tokyo. Kaguya(SELENE) [EB/OL]. [2009-02-25]. http://www.kaguya.jaxa.jp/en/index.htm.

JOLLIFF B L, GILLIS J J, HASKIN L A, KOROTEV R.L,WIECZOREK M A. 2000. Major Lunar Crustal Terranes: Surface Expressions and Crust-mantle Origins[J]. Journal of Geophysical Research, 105: 4197-4216.

KONOPLIV A S, SJOGREN W L, WIMBERLY R N, COOK R A,VIJAYARAGHAVAN A. 1993. A High Resolution Lunar Gravity Field and Predicted Orbit Behavior[C]//AAS/AIAA Astrodynamics Specialist Conference, Victoria, B C, Canada,August, AAS Paper 93-622.

KONOPLIV A S, BINDER A B, HOOD L L, KUCINSKAS A B,SJOGREN W L, WILLIAMS J G. 1998. Improved Gravity Field of the Moon from Lunar Prospector[J]. Association for the Advancement of Science, 281: 1476-1480.

KONOPLIV A S, YUAN D N. 1999. Lunar Prospector 100th Degree Gravity Model Development[R]. Houston: Lunar and Planetary Science Conference XXX: 1067.

KONOPLIV A S, ASMAR S W, CARRANZA E, SJOGREN W L,YUAN D N. 2001. Recent Gravity Models as a Result of the Lunar Prospector Mission[J]. Icarus, 150: 1-18.

LEMOINE F G, SMITH D E, ZUBER M T, NEUMANN G A,ROWLANDS D D. 1997. A 70th Degree Lunar Gravity Model(GLGM-2) form Clementine and Other Tracking Data[J].Journal of Geophysical Research, 102: 16339-16359.

LORELL J, SJOGREN W L. 1968. Lunar Gravity: Preliminary Estimates from Lunar Orbiter[J]. Science, 159: 625-628.

MICHAEL W H, BLACKSHEAR W T. 1972. Recent Results on the Mass, Gravitational Field and Moments of Inertia of the Moon[J]. Earth moon and Planets, 3: 388-402.

MULLER P M, SJOGREN W L. 1968. Mascons: Lunar Mass Concentrions[J]. Science, 161: 680-684.

NAMIKI N, IWATA T, MATSUMOTO K, HANADA H, NODA H,GOOSSENS S, OGAWA M, KAWANO N, ASARI K,TSURUTA S, ISHIHARA Y, LIU Q H, KIKUCHI F,ISHIKAWA T, SASAKI S, AOSHIMA C, KUROSAWA K,SUGITA S, TAKANO T. 2009. Farside Gravity Field of the Moon from Four-way Doppler Measurements of SELENE(Kaguya)[J]. Science, 323: 900-905.

WIECZOREK M A, JOLLIFF B L, KHAN A, PRITCHARD M E,WEISS B P, WILLIAMS J G, HOOD L L, RIGHTER K, NEAL C R, SHEARER C K, Mccallum I S, TOMPKINS S, HAWKE B R, PETERSON C, GILLIS J J, BUSSEY B. 2006. The Constitution and Structure of the Lunar Interior[J]. Review in Mineralogy and Geochemistry, 60: 221-364.

WILD-PFEIFFER F. 2008. A Comparison of Different Mass Elements for Use in Gravity Gradiometry[J]. Journal of Geodesy,82: 637-653.

WISE D U, YATES M T. 1970. Mascons as Structural Relief on a Lunar ‘Moho’[J]. Journal of Geophysical Research, 75:261-268.

ZUBER M T, SMITH D E, LEMOINE F G, NEUMANN G A. 1994.The Shape and Internal Structure of the Moon from the Clementine Mission[J]. Science, 266: 1839-1843.

Discussions on the Study of Lunar Gravity Field

YANG Hong-wei1), ZHAO Wen-jin1), WU Zhen-han1)
Chinese Academy of Geological Sciences, Beijing100037

Lunar gravimetry is one of the most useful approaches to understanding the interior constitution and structure of the moon and also serves as a fundamental basis for choice of landing sites in future. Of varied methods, satellite gravimetry is an important means for covering the whole moon and obtaining deep information.This paper has described the history of lunar gravity exploration and the development of lunar gravity models, and analyzed and comparatively studied these models. Based on knowledge of forefront research methods, the authors deal with the principle for solution of lunar gravity coefficient model, approaches to acquiring lunar gravity model with high accuracy, varied tools for gravity field reductions, and current heated topics in the study of the interior structure of the moon. In addition, characteristics of typical gravity anomalies in the moon are comprehensively analyzed.

lunar gravity exploration; gravity model of the moon; characters of typical gravity anomalies

P184.57; P223.4

A

10.3975/cagsb.2011.03.09

本文由中國地質(zhì)調(diào)查局地質(zhì)調(diào)查工作項目(編號: 1212010811050)資助。

2011-03-15; 改回日期: 2011-04-11。責(zé)任編輯: 魏樂軍。

楊宏偉, 男, 1982年生。博士研究生。現(xiàn)主要從事利用月球重力場研究月球內(nèi)部結(jié)構(gòu)特征。通訊地址: 100037, 北京市西城區(qū)百萬莊大街26號。電話: 010-68999652。E-mail: yhw1106@163.com。