Dynamical history of coplanar two-satellite systems

One of the possible early states of the Earth-Moon system was a system of several large satellites around the Earth. The dynamical evolution of coplanar three-body systems is studied; a planet (Earth) and two massive satellites (proto-moons) with geocentric orbits of slightly different radii. Such configurations may arise in multiple satellite systems receding from a planet due to tidal friction. The numerical integration of the equations of motion shows that initially circular Keplerian orbits are soon transformed into disturbed elliptic orbits which are intersecting. The life-time of such a coplanar system between two probable physical collisions of satellites is roughly from one day to one year for satellite systems with radii less than 20R⊕, and may reach 100 yr for three-dimensional systems. This time-scale is short in comparison with the duration of the removal of satellites due to tides raised on the planet, which is estimated as 10⁶–10⁸ yr for the same orbital dimensions. Therefore, the life-time of a system of several proto-moons is mainly determined by their tidal interactions with the Earth. For conditions which we have considered, the most probable result of the evolution was coalescence of satellites as the consequence of the collisions.     

Lunar paleotides and the origin of the Earth-Moon system

A new method for determining the early history of the Earth-Moon system is described. Called the study of lunar paleotides, it describes a method for explaining features of the remnant lunar gravity field, and the generation of the lunar mascons. A method for the determination of Earth-Moon distances compared with the radiometric ages of the maria is developed. It is shown that the Moon underwent strong anomalous gravitational tidal forces, for a durationt<10⁶yr, prior to the formation of the mascon surfaces. As these tidal forces had not been present at the time of the formation of the Moon, this shows that the Moon could not have been formed in orbit about the Earth.There are tides in the affairs of men which, taken at the flood, lead on to fortune... William Shakespeare 1564–1616     

On the accelerations of the Moon and Sun,the constant of gravitation,and the origin of mountains

In a previous paper Lyttleton (1976) has shown that the apparent secular accelerations of the Sun and Moon, as given by de Sitter, can be largely explained if the Earth is contracting at the rate required by the phase-change hypothesis for the nature of the core. More reliable values for these accelerations have since become available which warrant a redetermination of the various effects concerned on the basis of constantG, and this is first carried out in the present paper. The lunar tidal couple, which is the same whetherG is changing or not, is found to be (4.74±0.38)×10²³ cgs, about three-quarters that yielded by the de Sitter values, while within the theory the Moon would take correspondingly longer to reach close proximity to the Earth at about 1.5×10⁹ years ago.The more accurate values of the accelerations enable examination to be made of the effects that a decreasingG would have, and it is shown that a valueG/G=–3×10^–11 yr^–1 can be weakly satisfied compared with the close agreement found on the basis of constantG, while a value as large numerically asG/G=–6×10^–11 yr^–1 seems to be definitely ruled out. On the iron-core model, an intrinsic positive component of acceleration of the angular velocity cannot be reconciled at all with the secular accelerations even for constantG, and far less so ifG is decreasing at a rate suggested by any recent cosmological theory.ItG=0, the amount of contraction available for mountain-building would correspond to a reduction of surface area of about 49×10⁶ km² and a volume to be redistributed of 160×10⁹ km³ if the time of collapse were 2.5×10⁹ years ago. For earlier times, the values are only slightly reduced. IfG/G=–3×10^–11 yr^–1, the corresponding values are 44×10⁶ km² and 138×10⁹ km³ for collapse at –2.5×10⁹ yr, and not importantly smaller at 38×10⁶ km² and 122×10⁹ km³ for collapse at –4.5×10⁹ yr. Any of these values would suffice to account in order of magnitude for all the eras of mountain-building. An intense brief period of mountain-building on an immense scale would result from the Ramsey-collapse at whatever time past it may have occurred.     

Absolute dimensions of the close binary systems V453 Monocerotis and V523 Cassiopeiae

We present multi-colour CCD observations of the low-temperature contact binaries, V453 Mon and V523 Cas. Their light curves are modelled to determine a new set of stellar and orbital parameters. Analysis of mid-eclipse times yields a new linear ephemeris for both systems. A period decrease (dP/dt=2.3×10⁻⁷ days/yr) in V453 Mon is discovered. V523 Cas, however, is detected to show a period increase (dP/dt=9.8×10⁻⁸ days/yr) because of the mass transfer of a rate of 1.1×10⁻⁷ M_⊙ yr⁻¹, from a less massive donor. Using these findings we can determine the physical parameters of the components of V523 Cas to be M ₁=0.76 (3)M _⊙, M ₂=0.39 (2)M _⊙, R ₁=0.74 (2)R _⊙, R ₂=0.55 (2)R _⊙, L ₁=0.19 (3)L _⊙, L ₂=0.14 (3)L _⊙, and the distance of system as 46(9) pc.     

Some effects of elasticity on lunar rotation

A general Hamiltonian for a rotating Moon in the field of the Earth is expanded in terms of parameters orienting the spin angular momentum relative to the pricipal axes of the Moon and relative to coordinate axes fixed in the orbital plane. The effects of elastic distortion are included as modifications of the moment of inertia tensor, where the magnitude of the distortion is parameterized by the Love numberk ₂. The principal periodic terms in the longitude of a point on the Moon due to variations of the tide caused by the Earth are shown to have amplitudes between 3.9 × 10^–3 and 1.6 × 10^–2 with a period of an anomalistic month, 3.0 × 10^–4 and 1.2 × 10^–3 with a period of one-half an anomalistic month and 2.4 × 10^–4 and 9.6 × 10^–4 with a period of one-half of a nodical month. The extremes in the amplitudes correspond to rigidities of 8 × 10¹¹ cgs and 2 × 10¹¹ cgs, respectively, the former rigidity being comparable to that of the Earth. Only the largest amplitude given above is comparable to that detectable by the projected precision of the laser ranging to the lunar retrorereflectors, and this amplitude corresponds to an improbably low rigidity for the Moon. A detailed derivation of the free wobble of the lunar spin axis about the axis of maximum moment of inertia is given, where it is shown that elasticity can alter the period of the free wobble of 75.3 yr by only 3 × 10^–4 to 10^–3 of this period. Also, the effect of elasticity on the period of free libration is completely negligible by many orders of magnitude. If the Moon's rigidity is close to that of the Earth there is no effect of elasticity on the rotation which can be measured with the laser ranging and, therefore, no elastic properties of the Moon can be determined from variations in the rotation.Currently on leave from the Dept. of Physics, University of California, Santa, Barbara, California.Communication presented at the conference on Lunar Dynamics and Observational Coordinate Systems held January 15–17, 1973 at the Lunar Science Institute, Houston, Tex., U.S.A.     

Photometric Stability of the Lunar Surface

The rate at which cratering events currently occur on the Moon is considered in light of their influence on the use of the Moon as a radiometric standard. The radiometric effect of small impact events is determined empirically from the study of Clementine images. Events that would change the integral brightness of the moon by 1% are expected once per 1.4 Gyr. Events that cause a 1% shift in one pixel for low Earth-orbiting instruments with a 1-km nadir field of view are expected approximately once each 43 Myr. Events discernible at 1% radiometric resolution with a 5 arc-sec telescope resolution correspond to crater diameters of approximately 210 m and are expected once every 200 years. These rates are uncertain by a factor of two. For a fixed illumination and observation geometry, the Moon can be considered photometrically stable to 1 × 10⁻⁸per annum for irradiance, and 1 × 10⁻⁷per annum for radiance at a resolution common for spacecraft imaging instruments, exceeding reasonable instrument goals by six orders of magnitude.     

Hypervelocity Impact Experiments in the Laboratory Relating to Lunar Astrobiology

The results of a set of laboratory impact experiments (speeds in the range 1–5 km s⁻¹) are reviewed. They are discussed in the context of terrestrial impact ejecta impacting the Moon and hence lunar astrobiology through using the Moon to learn about the history of life on Earth. A review of recent results indicates that survival of quite complex organic molecules can be expected in terrestrial meteorites impacting the lunar surface, but they may have undergone selective thermal processing both during ejection from the Earth and during lunar impact. Depending on the conditions of the lunar impact (speed, angle of impact etc.) the shock pressures generated can cause significant but not complete sterilisation of any microbial load on a meteorite (e.g. at a few GPa 1–0.1% of the microbial load can survive, but at 20 GPa this falls to typically 0.01–0.001%). For more sophisticated biological products such as seeds (trapped in rocks) the lunar impact speeds generate shock pressures that disrupt the seeds (experiments show this occurs at approximately 1 GPa or semi-equivalently 1 km s⁻¹). Overall, the delivery of terrestrial material of astrobiological interest to the Moon is supported by these experiments, although its long term survival on the Moon is a separate issue not discussed here.     

Pre-capture orbits of the moon

If the mass of the Earth was not considerably larger than at present, the pre-capture orbit of the Moon was in the range 0.9–1.1 A.U. Capture occurred within several 10⁸ years after formation of the Moon.     

Astronomical evidence concerning non-gravitational forces in the Earth-Moon system

The most conspicuous effects of non-gravitational forces in the Earth-Moon system are the accelerations of the Earth's spin and of the Moon's mean angular velocity. Evidence indicates that the present acceleration of the Moon is between –20 and –52 s of arc per century per century and that the present average acceleration of the Earth is between –5 and –23 parts in 10⁹ per century. Over the past 2000 yr, the average for the Moon has been about –42 s per century per century and for the Earth has been about –28 parts in 10⁹ per century; these values are probably correct within 10%. Evidence that does not involve any assumptions about the present values shows strongly that there was a square wave in the accelerations that lasted from about 700–1300, and that the accelerations were different by a factor of perhaps 5 during the time of this wave from what they were at neighboring times.An effect that seems to be changing the obliquity of the ecliptic has been reported in recent literature, on the basis of data obtained within the past century. The effect amounts to about 1/4 s of arc per century if it is real. Older data are not accurate enough to give information about an effect this small.There are no satisfactory explanations of the accelerations. Existing theories of tidal friction are quite inadequate.Paper presented at the AAAS Symposium on the Early History of the Earth and Moon in Philadelphia on 28 December 1971.     

The large-number hypothesis and the Earth's expansion-II

Careful study by Owen of the displacements of the continents renders it likely that the Earth radius increased by almost 20% during the past 200 million years. Ramseys proposal that the Earth liquid core is a high-pressure phase of the mantle material is adopted to calculate the Earth radius under variableG. It is shown that ifG varies according to the large-numbers hypothesis of Dirac, the primordial Earth would have had radius 700 km less than the modern value. If it is assumed that the Earth began to respond and, hence, expand only 2 × 10⁸ yr B.P., the expansion as large as 700 km is shown to be energetically plausible.     

Satellite motion in a non-singular gravitational potential

We study the effects of a non-singular gravitational potential on satellite orbits by deriving the corresponding time rates of change of its orbital elements. This is achieved by expanding the non-singular potential into power series up to second order. This series contains three terms, the first been the Newtonian potential and the other two, here R ₁ (first order term) and R ₂ (second order term), express deviations of the singular potential from the Newtonian. These deviations from the Newtonian potential are taken as disturbing potential terms in the Lagrange planetary equations that provide the time rates of change of the orbital elements of a satellite in a non-singular gravitational field. We split these effects into secular, low and high frequency components and we evaluate them numerically using the low Earth orbiting mission Gravity Recovery and Climate Experiment (GRACE). We show that the secular effect of the second-order disturbing term R ₂ on the perigee and the mean anomaly are 4″.307×10⁻⁹/a, and −2″.533×10⁻¹⁵/a, respectively. These effects are far too small and most likely cannot easily be observed with today’s technology. Numerical evaluation of the low and high frequency effects of the disturbing term R ₂ on low Earth orbiters like GRACE are very small and undetectable by current observational means.     

Absolute parameters of the close binary BW Boo with a superasynchronous A2m primary component

Based on two high-dispersion spectra of the close binary BW Boo, we have detected lines of the secondary component whose contribution to the combined spectrum does not exceed 2%. We have determined the rotation velocities of the components and spectroscopic orbital elements. Numerous lines of neutral and ionized iron have been used to determine the effective temperature and surface gravity for the primary component. The photometric light curves for this binary have been solved for the first time. Its primary component is an A2Vm star with a mass of 2 ± 0.1M _⊙ and a radius of 1.9 ± 0.4R _⊙. Its rotation velocity is 2 km s⁻¹, which is a factor of 18 lower than the pseudo-synchronous velocity for this component. The G6 secondary component, a T Tau star, has a rotation velocity of 17 km s⁻¹, amass of 1.1M _⊙, and a radius of 1 R _⊙. The age of the binary has been estimated to be 10⁷ yr.     

Bias-corrected population, size distribution, and impact hazard for the near-Earth objects

Utilizing the largest available data sets for the observed taxonomic (Binzel et al., 2004, Icarus 170, 259-294) and albedo (Delbo et al., 2003, Icarus 166, 116-130) distributions of the near-Earth object population, we model the bias-corrected population. Diameter-limited fractional abundances of the taxonomic complexes are A-0.2%; C-10%, D-17%, O-0.5%, Q-14%, R-0.1%, S-22%, U-0.4%, V-1%, X-34%. In a diameter-limited sample, ∼30% of the NEO population has jovian Tisserand parameter less than 3, where the D-types and X-types dominate. The large contribution from the X-types is surprising and highlights the need to better understand this group with more albedo measurements. Combining the C, D, and X complexes into a “dark” group and the others into a “bright” group yields a debiased dark-to-bright ratio of ∼1.6. Overall, the bias-corrected mean albedo for the NEO population is 0.14±0.02, for which an H magnitude of 17.8±0.1 translates to a diameter of 1 km, in close agreement with Morbidelli et al. (2002, Icarus 158 (2), 329-342). Coupling this bias corrected taxonomic and albedo model with the H magnitude dependent size distribution of (Stuart, 2001, Science 294, 1691-1693) yields a diameter distribution with 1090±180 NEOs with diameters larger than 1 km. As of 2004 June, the Spaceguard Survey has discovered 56% of the NEOs larger than 1 km. Using our size distribution model, and orbital distribution of (Stuart, 2001, Science 294, 1691-1693) we calculate the frequency of impacts into the Earth and the Moon. Globally destructive collisions (∼10²¹ J) of asteroids 1 km or larger strike the Earth once every 0.60±0.1 Myr on average. Regionally destructive collisions with impact energy greater than 4×10¹⁸ J (∼200 m diameter) strike the Earth every 56,000±6000 yr. Collisions in the range of the Tunguska event (4-8×10¹⁶ J) occur every 2000-3000 yr. These values represent the average time between randomly spaced impacts; actual impacts could occur more or less closely spaced solely by chance. As a verification of these impact rates, the crater production function of Shoemaker et al. (1990, Geological Society of American Special Paper 247) has been updated by combining this new population model with a crater formation model to find that the observed crater production function on both the Earth and Moon agrees with the rate of crater production expected from the current population of NEOs.     

Effects of a physical librations of the moon caused by a liquid core,and determination of the fourth mode of a free libration

An analytical theory of lunar physical librations based on its two-layer model consisting of a non-spherical solid mantle and ellipsoidal liquid core is developed. The Moon moves on a high-precision orbit in the gravitational field of the Earth and other celestial bodies. The defined fourth mode of a free libration is caused by the influence of the liquid core, with a long period of 205.7 yr, with amplitude S = 0″0395 and with an initial phase Π₀ = ?134° (for the initial epoch 2000.0). Estimates of dynamic (meridional) oblatenesses of a liquid core of the Moon have been estimated: ?_D = 4.42 × 10^?4, μ_D = 2.83 × 10^?4 (?_D + μ_D = 7.24 × 10^?4). These results have been obtained as a result of comparison of the developed analytical theory of physical librations of the Moon with the empirical theory of librations of the Moon constructed on the basis of laser observations.     

Chemical composition of the Earth after the giant impact

The giant impact hypothesis for the origin of the Moon has been widely accepted. One of the most important features of this hypothesis is that the impactor's metallic core was incorporated in the Earth after impact. If the mass of the impactor is 0.82 × 10²⁷ g, the mass of the impactor core was estimated to be 0.19 × 10²⁷ g, which is about 1/10 of present Earth's core. Liu (1982) derived the bulk composition of the Earth from CI chondrites, and concluded that the Fe content of his model appears to be low in comparison with the present Earth, which, however, can be rationalized by the addition of impactor core into the proto-Earth developed by Liu (1982). If the impactor's mantle contains 14 wt% FeO as suggested, the mass ratio of impactor/proto-Earth should not exceed 0.22. The same ratio is not likely to exceed 0.30, if a giant blowoff did not occur during impact.     

On the origin of the moon by rotational fission

Based on simple CIPW norms for the proposed terrestrial upper mantle material, it is shown that if the Moon fissioned from the Earth and gravitationally differentiated, it could have a 72 km thick anorthosite (An₉₇) crust, a calcium poor (3.8% by weight) pyroxenite upper mantle 100 Mg/Mg + Fe = 75 to 80) ending at a depth of 313 km and a dunite (Fo_{93_95}) lower mantle below a depth of 313 km. Refinements of these simple norm models, based on the cooling history, crystallization sequence and the variations of the 100 Mg/Mg + Fe ratio of the liquid and crystals during the crystallization sequence, indicate that the final form of such a Moon could have the following properties: (1) a primitive, cumulate anorthosite - minor troctolite crust with intrusive and extrusive feldspathic basalts and KREEP rich norites; the thickness of this crust would be 75 km; (2) a zone in the bottom of the crust and the top of the upper mantle which is rich in KREEP, the incompatible elements, silica, and possibly voltiles; this zone would be the source area for the upland feldspathic basalts, KREEP rich norites and KREEP and silica rich fluids; (3) an upper mantle between the depths of 75 km and 350 to 400 km which consists of peridotite containing 80–85% pyroxene (Wo₁₀En_{68_72}Fs_{18_22}) and 15–20% olivine (Fo_{75_80}); the Al₂O₃ content of the upper mantle is 3%; the peridotite layer would be the source area for mare basalts and; (4) a lower mantle below a depth of 350–400 km which consists of dunite (Fo_{93_97}).The cooling history of such a moon indicates that the primitive anorthosite crust would have been completely formed within 10⁸ yr after fission. The extrusion and intrusion of upland basalts and KREEP rich norites and the metamorphism of the crustal rocks via KREEP and silica rich fluids would have ended about 4 × 10⁹ yr ago when cooling well below the solidus reached a depth of 150 km. As cooling continied, the only source of magmas after 4 × 10⁹ yr ago would have been the peridotite upper mantle, i.e. the source area of the mare basalts. Extrusion of mare basalts ended when cooling below the solidus reached the top of the refractory dunite lower mantle 3-3.3 × 10⁹ yr ago.Thus, it is shown that the chemistry, primary lithology, structure and developmental history of a fissioned Moon readily match those known for the real Moon. As such, the models presented in this paper strongly support the fission origin of the Moon.Guest Scientist, supported by the Alexander von Humboldt-Stiftung.Permanent Address.     

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.     

Hypothesis of variable g and the secular accelerations of the Sun and the Moon

It is known that the observed secular accelerations of the Sun and Moon are not consistent with the tidal interactions of the Earth with the Sun and Moon. Following Dicke, the hypothesis of variable constant of gravity is adopted and expressions for the accelerations are derived. It is shown that if the theoretical ratio of the acceleration is equated the observed one, a unique value for —/G can be calculated. Adopting the accelerations obtained by Fotheringham, Newton, Muller and Stephenson, and Stephenson, it is found that — /G ranges from 1.4 × 10^–11 to 3.3 × 10^–11 yr^–1. This estimate is consistent with the one based upon the comparison of the lunar accelerations measured with respect to atomic and ephemis times.     

The vector approach to the problem of physical libration of the Moon: the linearized problem

A flexible and informative vector approach to the problem of physical libration of the rigid Moon has been developed in which three Euler differential equations are supplemented by 12 kinematic ones. A linearized system of equations can be split into an even and odd systems with respect to the reflection in the plane of the lunar equator, and rotational oscillations of the Moon are presented by superposition of librations in longitude and latitude. The former is described by three equations and consists of unrestricted oscillations with a period of T ₁ = 2.878 Julian years (amplitude of 1.855″) and forced oscillations with periods of T ₂ = 27.201 days (15.304″), one stellar year (0.008″), half a year (0.115″), and the third of a year (0.0003″) (five harmonics altogether). A zero frequency solution has also been obtained. The effect of the Sun on these oscillations is two orders of magnitude less than that of the Earth. The libration in latitude is presented by five equations and, at pertrubations from the Earth, is described by two harmonics of unrestricted oscillations (T ₅ ≈ 74.180 Julian years, T ₆ ≈ 27.347 days) and one harmonic of forced oscillations (T ₃ = 27.212 days). The motion of the true pole is presented by the same harmonics, with the maximum deviation from the Cassini pole being 45.3″. The fifth (zero) frequency yields a stationary solution with a conic precession of the rotation axis (previously unknown). The third Cassini law has been proved. The amplitudes of unrestricted oscillations have been determined from comparison with observations. For the ratio $ \frac{{\sin I}} {{\sin \left( {I + i} \right)}} \approx 0.2311 $ \frac{{\sin I}} {{\sin \left( {I + i} \right)}} \approx 0.2311 , the theory gives 0.2319, which confirms the adequacy of the approach. Some statements of the previous theory are revised. Poinsot’s method is shown to be irrelevant in describing librations of the Moon. The Moon does not have free (Euler) oscillations; it has oscillations with a period of T ₅ ≈ 74.180 Julian years rather than T ≈ 148.167 Julian years.     

The Earth-figure perturbations in the Lunar theory

We compute the perturbations on the motion of the Moon due to the shape of the Earth. The zonal terms inJ ₂,J ₃, andJ ₄ are considered. The accuracy is estimated at 3×10^–5 and the results compared with previous theories.