Design of lunar landing trajectory under constraints

The design of a lunar landing trajectory which satisfies certain constraints is considered and discussed. The constraints are of two kinds, kinetic constraints, which deal with the relative positions among the Sun, the Moon, the Earth, the spacecraft and tracking stations, and dynamic constraints, which deal with the orbital motion of the spacecraft. After a discussion of the characteristics of lunar flight trajectory, a method of designing standard flight trajectory is suggested that satisfies the constraints. This method is applied to the Chinese lunar landing flight and to the pre-design of the orbit of a lunar satellite.     

日本SELENE月球探测计划和卫星间多普勒跟踪的数学模型

日本月球探测计划（SELENE）定于2004年夏季利用HIIa火箭发射一组共3颗绕月人造卫星。他们是主卫星、跟踪中断卫星和空间VLBI电波源。其主要科学目标之一是利用对绕月卫星的多普勒跟踪数据精确测定月球重力场,研究月球的起源与演化。SELENE计划中实现这个科学目标的关键技术是引入中继卫星,目的在于当处于低轨道的主卫星飞行到月球背面地面观测站无法观测时,采用卫星间跟踪方法（SST）,建立地面站与主卫星之间的联系,以得到月球背面重力场的直接测量数据。介绍了几种典型的四程卫星间多普勒跟踪模式和相应的数学模型,并针对SELENE计划中采用的特殊四程多普勒跟踪模式建立了卫星相对观测站速度与跟踪信号多普勒频移之间的转换关系。提出了利用GEODYNⅡ定轨分析软件处理SELENE多普勒跟踪数据的流程。     

月球卫星轨道力学综述

月球探测器的运动通常可分为3个阶段，这3个阶段分别对应3种不同类型的轨道：近地停泊轨道、向月飞行的过渡轨道与环月飞行的月球卫星轨道。近地停泊轨道实为一种地球卫星轨道；过渡轨道则涉及不同的过渡方式(大推力或小推力等)；环月飞行的月球卫星轨道则与地球卫星轨道有很多不同之处，它决不是地球卫星轨道的简单克隆。针对这一点，全面阐述月球卫星的轨道力学问题，特别是环月飞行中的一些热点问题，如轨道摄动解的构造、近月点高度的下降及其涉及的卫星轨道寿命、各种特殊卫星(如太阳同步卫星和冻结轨道卫星等)的轨道特征、月球卫星定轨等。     

Determining parameters of Moon’s orbital and rotational motion from LLR observations using GRAIL and IERS-recommended models

The aim of this work is to combine the model of orbital and rotational motion of the Moon developed for DE430 with up-to-date astronomical, geodynamical, and geo- and selenophysical models. The parameters of the orbit and physical libration are determined in this work from lunar laser ranging (LLR) observations made at different observatories in 1970–2013. Parameters of other models are taken from solutions that were obtained independently from LLR. A new implementation of the DE430 lunar model, including the liquid core equations, was done within the EPM ephemeris. The postfit residuals of LLR observations make evident that the terrestrial models and solutions recommended by the IERS Conventions are compatible with the lunar theory. That includes: EGM2008 gravitational potential with conventional corrections and variations from solid and ocean tides; displacement of stations due to solid and ocean loading tides; and precession-nutation model. Usage of these models in the solution for LLR observations has allowed us to reduce the number of parameters to be fit. The fixed model of tidal variations of the geopotential has resulted in a lesser value of Moon’s extra eccentricity rate, as compared to the original DE430 model with two fit parameters. A mixed model of lunar gravitational potential was used, with some coefficients determined from LLR observations, and other taken from the GL660b solution obtained from the GRAIL spacecraft mission. Solutions obtain accurate positions for the ranging stations and the five retroreflectors. Station motion is derived for sites with long data spans. Dissipation is detected at the lunar fluid core-solid mantle boundary demonstrating that a fluid core is present. Tidal dissipation is strong at both Earth and Moon. Consequently, the lunar semimajor axis is expanding by 38.20 mm/yr, the tidal acceleration in mean longitude is \(-25.90 {{}^{\prime \prime }}/\mathrm{cy}^2\), and the eccentricity is increasing by \(1.48\times 10^{-11}\) each year.     

Lunar perturbations of artificial satellites of the earth

Lunisolar perturbations of an artificial satellite for general terms of the disturbing function were derived by Kaula (1962). However, his formulas use equatorial elements for the Moon and do not give a definite algorithm for computational procedures. As Kozai (1966, 1973) noted, both inclination and node of the Moon's orbit with respect to the equator of the Earth are not simple functions of time, while the same elements with respect to the ecliptic are well approximated by a constant and a linear function of time, respectively. In the present work, we obtain the disturbing function for the Lunar perturbations using ecliptic elements for the Moon and equatorial elements for the satellite. Secular, long-period, and short-period perturbations are then computed, with the expressions kept in closed form in both inclination and eccentricity of the satellite. Alternative expressions for short-period perturbations of high satellites are also given, assuming small values of the eccentricity. The Moon's position is specified by the inclination, node, argument of perigee, true (or mean) longitude, and its radius vector from the center of the Earth. We can then apply the results to numerical integration by using coordinates of the Moon from ephemeris tapes or to analytical representation by using results from lunar theory, with the Moon's motion represented by a precessing and rotating elliptical orbit.     

Electron content near the lunar surface using dual-frequency VLBI tracking data in a single lunar orbiter mission

In VLBI observations of Vstar, a subsatellite of the Japanese lunar mission SELENE, there were opportunities for lunar grazing occultation when Vstar was very close to the limb of the Moon. This kind of chance made it possible to probe the thin plasma layer above the Moon's surface as a meaningful by-product of VLBI,by using the radio occultation method with coherent radio waves from the S/X bands.The dual-frequency measurements were carried out at Earth-based VLBI stations. In the line-of-sight direction between the satellite and the ground-based tracking station where VLBI measurements were made, the effects of the terrestrial ionosphere, interplanetary plasma and the thin lunar ionosphere mixed together in the combined observables of dual-frequency Doppler shift and phase shift. To separate the variation of the ionospheric total electron content(TEC) near the surface of the Moon from the mixed signal, the influences of the terrestrial ionosphere and interplanetary plasma have been removed by using an extrapolation method based on a short-term trend. The lunar TEC is estimated from the dual-frequency observation for Vstar from UT 22:18to UT 22:20 on 2008 June 28 at several tracking stations. The TEC results obtained from VLBI sites are identical, however, they are not as remarkable as the result obtained at the Usuda deep space tracking station.     

A Proposed lunar orbiting gravity gradiometer experiment

Analysis of the gravity gradiometer developed by R. L. Forward and C. C. Bell at the Hughes Research Laboratories suggest than an accuracy, in the range 0.1 to 0.5 EU can be expected in a lunar orbiter application. This accuracy will allow gradient anomalies associated with mascons to be mapped with 1% accuracy and should reveal a great deal of new information about the lunar gravity field.The proposed experiment calls for putting such a gradiometer into a closely circular polar orbit at an average height of about 30 km above the lunar surface. This orbit allows the entire lunar surface to be covered in fourteen days, the gradiometer to be checked twice per revolution and results in successive passes above the lunar surface being spaced at about the resolution limit of about 30 km set both by the satellite altitude and instrumental integration time. Doppler tracking will be employed and the spacecraft will carry an electromagnetic altimeter. Gradient and altitude data from the far side of the Moon can be stored for replay when communication is re-established.     

On the deviation of the lunar center of mass to the East. Two possible mechanisms based on evolution of the orbit and rounding off the shape of the Moon

From the observations of the gravitational field and the figure of the Moon, it is known that its center of mass (briefly COM) does not coincide with the center of figure (COF), and the line “COF/COM” is not directed to the center of the Earth, but deviates from it to the South–East. Here we study the deviation of the lunar COM to the East from the mean direction to Earth.At first, we consider the optical libration of a satellite with synchronous rotation around the planet for an observer at a point on second (empty) orbit focus. It is found that the main axis of inertia of the satellite has asymmetric nonlinear oscillations with amplitude proportional to the square of the orbit eccentricity. Given this effect, a mechanism of tidal secular evolution of the Moon’s orbit is offered that explains up to \(20\%\) of the known displacement of the lunar COM to the East. It is concluded that from the alternative—evolution of the Moon’s orbit with a decrease or increase in eccentricity—only the scenario of evolution with a monotonous increase in orbit eccentricity agrees with the displacement of lunar COM to the East. The precise calculations available confirm that now the eccentricity of the lunar orbit is actually increasing and therefore in the past it was less than its modern value, \(e = 0.0549\).To fully explain the displacement of the Moon’s COM to the East was deduced a second mechanism, which is based on the reliable effect of tidal changes in the shape of the Moon. For this purpose the differential equation which governs the process of displacement of the Moon’s COM to the East with inevitable rounding off its form in the tidal increase process of the distance between the Earth and the Moon is derived. The second mechanism not only explains the Moon’s COM displacement to the East, but it also predicts that the elongation of the lunar figure in the early epoch was significant and could reach the value \(\varepsilon\approx0.31\). Applying the theory of tidal equilibrium figures, we can estimate how close to the Earth the Moon could have formed.     

Lunar gravity derived from long-period satellite motion — A proposed method

A new method has been devised to determine the spherical harmonic coefficients of the lunar gravity field. This method consists of a two-step data reduction and estimation process. In the first step, a weighted least-squares empirical orbit determination scheme is applied to Doppler tracking data from lunar orbits to estimate longpperiod Kepler elements and rates. Each of the Kepler elements is represented by an independent function of time. The long-period perturbing effects of the Earth, Sun, and solar radiation are explicitly modeled in this scheme. Kepler element variations estimated by this empirical processor are then ascribed to the non-central lunar gravitation features. Doppler data are reduced in this manner for as many orbits as are available. In the second step, the Kepler element rates are used as input to a second least-squares processor that estimates lunar gravity coefficients using the long-period Lagrange perturbation equations.Pseudo Doppler data have been generated simulating two different lunar orbits. This analysis included the perturbing effects of the L1 lunar gravity field, the Earth, the Sun, and solar radiation pressure. Orbit determinations were performed on these data and long-period orbital elements obtained. The Kepler element rates from these solutions were used to recover L1 lunar gravity coefficients. Overall results of this controlled experiment show that lunar gravity coefficients can be accurately determined and that the method is dynamically consistent with long-period perturbation theory.     

Cratering rate over the surface of a synchronous satellite

Equations have been derived for the asymmetries of crater frequency over the surface of a synchronously rotating satellite, when the orbital velocities of projectiles about the parent planet are always larger than the satellite's circular orbital velocity. If the projectiles orbit the planet in moderately eccentric ellipses, no marked apex-antapex asymmetries of crater frequency distribution are expected. Theoretical values of apex-antapex asymmetries are presented for the Earth, the Moon, and some Jovian and Saturnian satellites, and are compared with available observational and theoretical results.     

Study of the transfer between libration point orbits and lunar orbits in Earth–Moon system

This paper is devoted to the study of the transfer problem from a libration point orbit of the Earth–Moon system to an orbit around the Moon. The transfer procedure analysed has two legs: the first one is an orbit of the unstable manifold of the libration orbit and the second one is a transfer orbit between a certain point on the manifold and the final lunar orbit. There are only two manoeuvres involved in the method and they are applied at the beginning and at the end of the second leg. Although the numerical results given in this paper correspond to transfers between halo orbits around the \(L_1\) point (of several amplitudes) and lunar polar orbits with altitudes varying between 100 and 500 km, the procedure we develop can be applied to any kind of lunar orbits, libration orbits around the \(L_1\) or \(L_2\) points of the Earth–Moon system, or to other similar cases with different values of the mass ratio.     

Origin of the moon by tidal capture and some geophysical consequences

Of the many proposed modes of origin of the Moon, some violate physical laws; many are in conflict with observations; all are improbable. Perhaps the least improbable - based on recent tidal theory calculations and on the interpretation of lunar rock data - is capture of the Moon as it passed near the Earth in adirect (prograde) orbit, shortly after the formation of Moon and Earth, about 4.5 billion years ago. (Capture of the Moon from an initiallyretrograde orbit which had been proposed some years ago, leads to physically unacceptable consequences.) The effects of capture on the Earth would have been cataclysmic, leading to intensive heating of its interior, to volcanism, and to the immediate formation of an atmosphere and hydrosphere. Thus capture of a Moon may have given rise to the unique properties of the Earth (in the Solar System) and to the early evolution of life, about 3.5 billion years ago.Presented at the NATO Advanced Study Institute on Lunar Studies in Patras, Greece, September, 1971.     

Direct and indirect capture of near-Earth asteroids in the Earth–Moon system

Near-Earth asteroids have attracted attention for both scientific and commercial mission applications. Due to the fact that the Earth–Moon \(\hbox {L}_{1}\) and \(\hbox {L}_{2}\) points are candidates for gateway stations for lunar exploration, and an ideal location for space science, capturing asteroids and inserting them into periodic orbits around these points is of significant interest for the future. In this paper, we define a new type of lunar asteroid capture, termed direct capture. In this capture strategy, the candidate asteroid leaves its heliocentric orbit after an initial impulse, with its dynamics modeled using the Sun–Earth–Moon restricted four-body problem until its insertion, with a second impulse, onto the \(\hbox {L}_{2}\) stable manifold in the Earth–Moon circular restricted three-body problem. A Lambert arc in the Sun-asteroid two-body problem is used as an initial guess and a differential corrector used to generate the transfer trajectory from the asteroid’s initial obit to the stable manifold associated with Earth–Moon \(\hbox {L}_{2}\) point. Results show that the direct asteroid capture strategy needs a shorter flight time compared to an indirect asteroid capture, which couples capture in the Sun–Earth circular restricted three-body problem and subsequent transfer to the Earth–Moon circular restricted three-body problem. Finally, the direct and indirect asteroid capture strategies are also applied to consider capture of asteroids at the triangular libration points in the Earth–Moon system.     

A new variant of the large-impact hypothesis for the origin of the Moon

The hypothesis is advanced that after collision of a Mars-sized impact with the Earth, collisions between debris particles themselves are able to place enough material into Earth orbit, to form the Moon. Collision probability estimates show that the collision frequency is high enough to place about one lunar mass into Earth orbit, if the average semimajor axis is of order of the Earth's Roche limit of 18 500 km.     

Numerical simulations of a Mars geodesy network experiment: Effect of orbiter angular momentum desaturation on Mars' rotation estimation

The scientific objectives of a geodetic experiment based on a network of landers, such as NEIGE (NEtlander Ionosphere and Geodesy Experiment) are to improve the current knowledge of Mars' interior and atmosphere dynamics. Such a network science experiment allows monitoring the motions of the Martian rotation axis with a precision of a few centimeters (or milli-arc-seconds (mas)) over annual and sub-annual periods. Thereto, besides radio tracking of a Mars orbiter from the Earth, radio Doppler shifts between this orbiter and several landers at the planet's surface will be performed. From the analysis of these radio Doppler data, it is possible to reconstruct the orbiter motion and Mars' orientation in space. The errors on the orbit determination (position and velocity of the orbiter) have an impact on the geodetic parameters determination from the Doppler shifts and must be removed from the signal in order to achieve a high enough accuracy. In this paper, we perform numerical simulations of the two Doppler signals involved in such an experiment to estimate the impact of the spacecraft angular momentum desaturations on the determination of Mars' orientation variations. The attitude control of the orbiter needs such desaturation maneuvers regularly repeated. They produce velocity variations that must be taken into account when determining the orbit. For our simulations, we use a priori models of the Martian rotation and introduce the spacecraft velocity variations induced by each desaturation event. By a least-squares adjustment of the simulated Doppler signals, we then estimate the orbiter velocity variations and the parameters of the Mars' rotation model. We show that these velocity variations are ill resolved when the spacecraft is not tracked, therefore requiring a near-continuous tracking from the Earth to accurately determine the orbit. In such conditions we show that only 15- of lander-orbiter tracking per week allows recovering Mars' orientation parameters with a precision of a few mas over a period of 1 Martian year.     

Luni-solar perturbations of an Earth satellite

Luni-solar perturbations of the orbit of an artificial Earth satellite are given by modifying the analytical theory of an artificial lunar satellite derived by the author in recent papers. Expressions for the first-order changes, both secular and periodic, in the elements of the geocentric Keplerian orbit of the earth satellite are given, the moon's geocentric orbit, including solar perturbations in it, being found by using Brown's lunar theory.The effects of Sun and Moon on the satellite orbit are described to a high order of accuracy so that the theory may be used for distant earth satellites.     

On the possibility of investigating the galactic center region with the integral observatory by the method of lunar and terrestrial occultations

The possibility of investigating the sky region near the Galactic center with instruments of the INTEGRAL orbital astrophysical gamma-ray observatory by the method of its occultation by the Earth and the Moon is considered. Existing engineering constraints on the observing conditions, such as the admissible orientation of the INTEGRAL satellite relative to the direction to the Sun and the performance of measurements only outside the Earth??s radiation belts, are taken into account. Long time intervals during which the lunar occultation center passes at angular distances of less than 2° from the Galactic center have been found. Such events occur under the adopted constraints two or three times per year without any correction of the INTEGRAL satellite orbit. The orbit can be corrected to reduce the angular distance between the Moon and the Galactic center in occultation events. The required velocity impulses do not exceed several meters per second. The possibility of the Galactic center being occulted by the Earth has been analyzed. In this case, to perform measurements, the admissible (in radiation exposure) height of the working segment of the orbit should be reduced to 25 000 km, which can be problematic. At the same time, part of the Galaxy??s equatorial region is shadowed by the Earth for a time long enough to carry out the corresponding experiments.     

Recent results on the mass,gravitational field and moments of inertia of the moon

Doppler tracking data from the Lunar Orbiter series of spacecraft have been used in a more complete analysis of the spherical harmonic coefficients of the lunar gravitational field through thirteenth degree and order. The value obtained for the mass of the Moon,GM = 4902.84 km³ s^–2, is in good agreement with previous results and with results obtained by alternate procedures. Acceleration contour plots, derived from the gravitational coefficients, show correlations with surface features on the near side of the Moon, but are of questionable validity for the far side because of the lack of direct tracking data on the far side. Based on the most recent gravitational field data, the current estimate for the polar moment of inertia of the Moon isC/Ma ² = 0.4019_-0.002 ^+0.004. This value indicates that the interior of the Moon can be homogeneous, but some results presented strongly suggest that the Moon is differentiated, with an excess of mass in the direction toward the Earth.Paper presented at the NATO Advanced Study Institute on Lunar Studies, Patras, Greece, September, 1971.     

Formation of Embryos of the Earth and the Moon from a Common Rarefied Condensation and Their Subsequent Growth

Embryos of the Moon and the Earth may have formed as a result of contraction of a common parental rarefied condensation. The required angular momentum of this condensation could largely be acquired in a collision of two rarefied condensations producing the parental condensation. With the subsequent growth of embryos of the Moon and the Earth taken into account, the total mass of as-formed embryos needed to reach the current angular momentum of the Earth–Moon system could be below 0.01 of the Earth mass. For the low lunar iron abundance to be reproduced with the growth of originally iron-depleted embryos of the Moon and the Earth just by the accretion of planetesimals, the mass of the lunar embryo should have increased by a factor of 1.3 at the most. The maximum increase in the mass of the Earth embryo due to the accumulation of planetesimals in a gas-free medium is then threefold, and the current terrestrial iron abundance is not attained. If the embryos are assumed to have grown just by accumulating solid planetesimals (without the ejection of matter from the embryos), it is hard to reproduce the current lunar and terrestrial iron abundances at any initial abundance in the embryos. For the current lunar iron abundance to be reproduced, the amount of matter ejected from the Earth embryo and infalling onto the Moon embryo should have been an order of magnitude larger than the sum of the overall mass of planetesimals infalling directly on the Moon embryo and the initial mass of the Moon embryo, which had formed from the parental condensation, if the original embryo had the same iron abundance as the planetesimals. The greater part of matter incorporated into the Moon embryo could be ejected from the Earth in its multiple collisions with planetesimals (and smaller bodies).     

Zonal harmonics of azimuthally averaged potential of rotating inhomogeneous ellipsoids

The distant retrograde orbits (DROs) can serve as the parking orbits for a long-term cis-lunar space station. This paper gives a comprehensive study on the transfer problem from DROs to Earth orbits, including low Earth orbits (LEOs), medium Earth orbits (MEOs), and geosynchronous orbits (GSOs), in the bicircular restricted four-body problem (BR4BP) via optimizations within a large solution space. The planar transfer problem is firstly solved by grid search and optimization techniques, and two types of transfer orbits, direct ones and low-energy ones, are both constructed. Then, the nonplanar transfer problem to Earth orbits with inclinations between 0 and 90 degrees are solved via sequential optimizations based on the planar transfers. The transfer characteristics in the cases of different destination orbit inclinations are discussed for both the direct and the low-energy transfer orbits. The important role of the lunar gravity in the low-energy transfers is also discussed, which can overcome the increase of transfer cost caused by the high inclination of Earth orbits. The distinct features of different transfer scenarios, including multiple revolutions around the Earth and Moon, the exterior phase, and the lunar flyby, are discovered. The energy of transfer orbits is exploited to discuss the effects of close lunar flybys. The results will be helpful for the transfer design in future manned or unmanned return missions, and can also provide valuable information for selecting proper parking DROs for cis-lunar space stations.