月球重力场可用来研究月球演化过程和内部结构，是影响绕月卫星精密定轨的重要因素。基于GRAIL任务数据解算的GL0660B重力场模型，极大提高了月球重力场空间频谱信号的强度和范围。本文首先通过计算相应重力场的阶方差和地形相关性分析，对GL0660B模型进行了精度分析；其次，利用GL0660B模型和其他几个月球重力场模型进行比较，对月球重力场的特征进行了分析；然后通过绘制GL0660B模型和LP150Q模型在月球外部不同高度处的重力异常图，分析比较了月球重力场模型在不同高度上所反映的月球重力场的特征和差异；最后，利用GEODYN软件模拟计算了不同高度卫星的轨道变化。可以看出绕月卫星离心率随时间的变化，以及周期性变化趋势，而且不同高度卫星轨道处质量瘤的摄动影响不同，远月点、近月点和偏心率的变化也存在差异。
The lunar gravity field provides a way to research moon's evolution and probes the interior structure of the moon. It is an important factor influencing the lunar satellite precise orbit determination as well. The new lunar gravity model GL0660B from GRAIL mission dramatically improves the gravity spectrum and spectral ranges. Using the model GL0660B, it can be computed that the corresponding degree-wise RMS and correlation of topography, with which the quality of model GL0660B can be analyzed. Then different characters of the lunar gravity field comparing with other lunar gravity fields are analyzed. Besides, gravity anomaly distribution figures at different height of the models are given, and the character and difference of the lunar gravity models at different height are compared. In addition, lunar satellite orbit revolutionary at different height are modeled by GEODYN. The result shows that the trend of lunar satellite eccentricity changes is a complex and long cycle of change trend. It is different affected by the perturbation of the mascons of different height, which causes different changes of apolune, perilune and eccentricity.
[1] KONOPLIV A S, ASMAR S W, CARRANZA E, et al. Recent Gravity Models as a Result of the Lunar Prospector Mission[J]. Icarus, 2001, 150(1):1-18.
[2] KONOPLIV A S, BINDER A B, HOOD L L, et al. Improved Gravity Field of the Moon from Lunar Prospector[J]. Science, 1998, 281(5382):1476-1480.
[3] NAMIKI N, IWATA T, MATSUMOTO K, et al. Farside Gravity Field of the Moon from Four-way Doppler Measurements of SELENE (Kaguya)[J]. Science, 2009, 323(5916):900-905.
[4] GOOSSENS S, MATSUMOTO K, LIU Q, et al. Lunar Gravity Field Determination Using SELENE Same-beam Differential VLBI Tracking Data[J]. Journal of Geodesy, 2011, 85(4):205-228.
[5] ROWLANDS D D, MARSHALL J A, MCCARTHY J, et al. GEODYN Ⅱ System Description. Vols.1-5[R]. Contractor Report, Hughes STX Corp. Greenbelt, MD:[s.n.], 1997.
[6] KAULA W M. Theory of Satellite Geodesy:Applications of Satellites to Geodesy[M]. London:Blaisdell, 1966.
[7] LEMOINE F G, GOOSSENS S, SABAKA T J, et al. High-degree Gravity Models from GRAIL Primary Mission Data[J]. Journal of Geophysical Research:Planets, 2013, 118(8):1676-1698. DOI:10.1002/jgre.20118.
[8] NEUMANN G A, ZUBER M T, SMITH D E, et al. The Lunar Crust:Global Structure and Signature of Major Basins[J]. Journal of Geophysical Research:Planets, 1996, 101(E7):16841-16863.
[9] MULLER P M, SJOGREN W L. Mascons:Lunar Mass Concentrations[J]. Science, 1968, 161(3842):680-684.
[10] HEISKANEN W A, MORITZ H. Physical Geodesy[M]. San Francisco:Freeman, 1967.
[11] LIU Lin, WANG Jiasong. An Analytic Solution of the Orbital Variation of Lunar Satellites[J]. Chinese Astronomy and Astrophysics, 1998, 22(3):328-351.
[12] KONOPLIV A S, PARK R S, YUAN D N, et al. The JPL Lunar Gravity Field to Spherical Harmonic Degree 660 from the GRAIL Primary Mission[J]. Journal of Geophysical Research:Planets, 2013, 118(7):1415-1434.
[13] ANDREWS-HANNA J C, ASMAR S W, HEAD Ⅲ J W, et al. Ancient Igneous Intrusions and Early Expansion of the Moon Revealed by GRAIL Gravity Gradiometry[J]. Science, 2013, 339(6120):675-678.
[14] ASMAR S W, KONOPLIV A,S,WATKINS M M, et al. The Scientific Measurement System of the Gravity Recovery and Interior Laboratory (GRAIL) Mission[J]. Space Science Reviews, 2013, 178(1):25-55. DOI:10.1007/s11214-013-9962-0.
[15] SMITH D E, ZUBER M T, NEUMANN G A, et al. Initial Observations from the Lunar Orbiter Laser Altimeter (LOLA)[J]. Geophysical Research Letters, 2010, 37(18):L18204. DOI:10.1029/2010GL043751.
[16] ARNOLD D, BERTONE S, JÄGGI A, et al. GRAIL Gravity Field Determination Using the Celestial Mechanics Approach[J]. Icarus, 2015, 261:182-192. DOI:10.1016/j.icarus.2015.08.015.
[17] 孙玉, 常晓涛, 郭金运, 等. 由SGM100i质量分析看SELENE的贡献[J]. 测绘科学, 2012, 37(2):176-178. SUN Yu, CHANG Xiaotao, GUO Jinyun, et al. Contribution of SELENE for Lunar Gravity Field Recovery from Quality Analysis of SGM100i[J]. Science of Surveying and Mapping, 2012, 37(2):176-178.
[18] LI Peijia, HU Xiaogong, HUANG Yong, et al. Orbit Determination for Chang'E-2 Lunar Probe and Evaluation of Lunar Gravity Models[J]. Science China Physics, Mechanics and Astronomy, 2012, 55(3):514-522.
[19] JIN Shuanggen, ARIVAZHAGAN S, ARAKI H. New Results and Questions of Lunar Exploration from SELENE, Chang'E-1, Chandrayaan-1 and LRO/LCROSS[J]. Advances in Space Research, 2013, 52(2):285-305. DOI:10.1016/j.asr.2012.11.022.
[20] 丰海, 李建成, 李大炜, 等. 月球Airy均衡状态与月壳厚度估计[J]. 测绘学报, 2012, 41(4):543-548. FENG Hai, LI Jiancheng, LI Dawei, et al. Airy Isostasy on the Moon:Implications for Crustal Thickness[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(4):543-548.
[21] 孙雪梅, 李斐, 鄢建国, 等. Slepian函数在月球局部重力场分析中的适用性分析[J]. 测绘学报, 2015, 44(3):264-273. DOI:10.11947/j.AGCS.2015.20130728. SUN Xuemei, LI Fei, YAN Jianguo, et al. An Analysis of the Applicability of Slepian Function in Analyzing Lunar Local Gravity Field[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(3):264-273. DOI:10.11947/j.AGCS.2015.20130728.
