Mercury's atmosphere: A perspective after Mariner 10

Measurements made during the Mariner 10 flybys of Mercury have shown that this planet has a tenuous atmosphere, somewhat similar to that of the Moon, which consists of at least helium and can be classified as an exosphere. The amount of helium observed can be supplied by either the accretion of only a fraction of the solar wind He²⁺ diffusing across the magnetopause, or from outgassing of radiogenic helium from the planetary crust. The role of solar wind in the maintenance and depletion of Mercury's atmosphere is discussed in view of the density upper limits established from Mariner 10. The argon supply rate on Mercury is probably not more than that on the Earth, but it is difficult to say whether Mercury is deficient in potassium or not on the basis of the present data. The global outgassing of CO₂ and H₂O from the planet interior is estimated to be at least four orders of magnitude smaller than for Earth which indicates that either Mercury is deficient in volatiles or that this planet is very inactive.     

Models of dense cores in translucent regions of low mass star formation

We have constructed models for a region of low mass star formation where stellar winds ablate material from dark dense cores and return it to a translucent intercore medium from which subsequent generations of cores condense. Depletion of gas phase species onto grains plays a major role in the chemistry. For reasonable agreement between model core chemical fractional abundances and measured TMC-1 fractional abundances to obtain, the core collapse, once started, must be relatively uninhibited by turbulence or magnetic fields and the core lifetime must fall in a limited range determined by the assumed depletion rates. In a core with the TMC-1 fractional abundances, CH, OH, C₂H, H₂CO, HCN, HNC, and CN are the only simple species that have been detected in TMC-1 at radio and millimeter wavelengths to have fractional abundances that are roughly constant or increasing with time; this result bears considerably on previous work concerned with searches for spectroscopic evidence for and the diagnosis of collapse during protostellar formation, but depends on the fractions of the OH and CH emissions that are associated with the core centre rather than more extended gas or a core-stellar wind boundary layer. Model results for the abundance ratios of H₂O, CH₄, and NH₃ ices are in good agreement with those inferred for Halley's Comet.     

The energy flux number and three types of planetary dynamo

It is proposed that the existence and nature of a planetary dynamo can be characterized by a dimensionless number Φ ≡ F_eR/ϱλ²ω, called the energy flux number, where F_e is the energy flux available for dynamo generation, R is the core radius (or thickness of the dynamo generating region), ϱ is the fluid density, λ is the magnetic diffusivity and ω is the angular velocity. For Φ ≲ 1, there is no dynamo. For 1 ≲ Φ ≲ 10^2.5 there is an “energy-limited dynamo”, in which F_e is insufficient to enable the dynamo to reach the dynamically desirable state A ≡ B²/8πϱλω ∼ 1, where B is a typical field amplitude (in Gauss). For 10^2.5 ≲ Φ ≲ 10⁵, there is a dynamically determined dynamo (Λ ∼ 1) in which the magnetic Reynolds number of turbulent eddies is small. For Φ ≳ 10⁵, there is a turbulent dynamo. Probable planetary examples of these three dynamo states are Mercury (Φ ∼ 10²-10³), Earth (Φ ∼ 10⁴) and Jupiter (Φ ∼ 10¹¹), respectively.     

The 3:2 spin-orbit resonant motion of Mercury

Our purpose is to build a model of rotation for a rigid Mercury, involving the planetary perturbations and the non-spherical shape of the planet. The approach is purely analytical, based on Hamiltonian formalism; we start with a first-order basic averaged resonant potential (including J ₂ and C ₂₂, and the first powers of the eccentricity and the inclination of Mercury). With this kernel model, we identify the present equilibrium of Mercury; we introduce local canonical variables, describing the motion around this 3:2 resonance. We perform a canonical untangling transformation, to generate three sets of action-angle variables, and identify the three basic frequencies associated to this motion. We show how to reintroduce the short-periodic terms, lost in the averaging process, thanks to the Lie generator; we also comment about the damping effects and the planetary perturbations. At any point of the development, we use the model SONYR to compare and check our calculations.     

2000—2049年行星天象（一）

本文给出 2 0 0 0 - 2 0 4 9年行星天象的一部分 ,包括水星、金星的东、西大距 ,留及上、下合日 ;金星最亮和最接近地球 ;外行星冲日、最接近地球、合日、留 ;地球和外行星过远近日点计1 5个表     

A procedure for determining the nature of Mercury's core

Abstract— We review the assertion that the precise measurement of the second degree gravitational harmonic coefficients, the obliquity, and the amplitude of the physical libration in longitude, C₂₀, C₂₂, θ, and φ₀, for Mercury are sufficient to determine whether or not Mercury has a molten core (Peale, 1976). The conditions for detecting the signature of the molten core are that such a core not follow the 88‐day physical libration of the mantle induced by periodic solar torques, but that it does follow the 250 000‐year precession of the spin axis that tracks the orbit precession within a Cassini spin state. These conditions are easily satisfied if the coupling between the liquid core and solid mantle is viscous in nature. The alternative coupling mechanisms of pressure forces on irregularities in the core‐mantle boundary (CMB), gravitational torques between an axially asymmetric mantle and an assumed axially asymmetric solid inner core, and magnetic coupling between the conducting molten core and a conducting layer in the mantle at the CMB are shown for a reasonable range of assumptions not to frustrate the first condition while making the second condition more secure. Simulations have shown that the combination of spacecraft tracking and laser altimetry during the planned MESSENGER (MErcury Surface, Space ENvironment, GEochemistry, Ranging) orbiter mission to Mercury will determine C₂₀, C₂₂, and θ to better than 1% and φ₀ to better than 8%—sufficient precision to distinguish a molten core and constrain its size. The possible determination of the latter two parameters to 1% or less with Earth‐based radar experiments and MESSENGER determination of C₂₀ and C₂₂ to 0.1% would lead to a maximum uncertainty in the ratio of the moment of inertia of the mantle to that of the whole planet, C_m/C, of ?2% with comparable precision in characterizing the extent of the molten core.     

Recent advances in planetary magnetism

During the past decade, significant advances in thein situ measurements of planetary magnetic fields have been made. The U.S.A. and U.S.S.R. have conducted spacecraft investigations of all the planets, from innermost Mercury out to Jupiter. Unexpectedly, Mercury was found to possess a global magnetic field but neither the Moon nor Venus do. The results at Mars are incomplete butif a global field exists, it is clearly quite weak. The main magnetic field of Jupiter has been measured directly for the first time and confirms, as well as augments appreciably, the past 2 decades of groundbased radio astronomical studies which provided indirect evidence of the field. Progress in developing analytically complete models of the dynamo process suggests a possible common origin for Mercury, Earth and Jupiter.Paper dedicated to Professor Hannes Alfvén on the occasion of his 70th birthday, 30 May, 1978.     

Love numbers of the moon and of the terrestrial planets

In the IERS Standards (1989), for the Moon the adopted value of the tide Love number, k ₂, is equal to 0.0222. In this paper using the latest geodetic parameters of the Moon a group of internal structure models are constructed for this celestial body (see Table V), then the dependence of the Moon's core size on calculated value of k ₂ is explored. The obtained results indicate that the second degree Love number, k ₂ = 0.02664, of the lunar model 91–04 is near its observed value (0.027 ± 0.006). This implies that the Moon may possess an outer core of 660 km radius and of 300 kbar mean rigidity. With the same method the static Love numbers from degree 2 to 30 are computed for the terrestrial planets — Mercury, Venus, and Mars (see Table VII), and the influence of some parameters (such as the rigidity) of the outer core on low degree Love numbers is discussed. Finally, the likely range of the second degree Love numbers is determined for the terrestrial planets (see Table XI). It seems that if low degree Love numbers of a terrestrial planet can be detected in the future space explorations, there is some possibility to improve the planetary internal structure model. For example, as soon as space techniques yield an observed value of k ₂ > 0.10 for Mercury, there will be reason to anticipate that a partly melted iron core exists in this planet.     

Did 26Al and impact‐induced heating differentiate Mercury?

Numerical models dealing with the planetary scale differentiation of Mercury are presented with the short‐lived nuclide, ²⁶Al, as the major heat source along with the impact‐induced heating during the accretion of planets. These two heat sources are considered to have caused differentiation of Mars, a planet with size comparable to Mercury. The chronological records and the thermal modeling of Mars indicate an early differentiation during the initial ~1 million years (Ma) of the formation of the solar system. We theorize that in case Mercury also accreted over an identical time scale, the two heat sources could have differentiated the planets. Although unlike Mars there is no chronological record of Mercury's differentiation, the proposed mechanism is worth investigation. We demonstrate distinct viable scenarios for a wide range of planetary compositions that could have produced the internal structure of Mercury as deduced by the MESSENGER mission, with a metallic iron (Fe‐Ni‐FeS) core of radius ~2000 km and a silicate mantle thickness of ~400 km. The initial compositions were derived from the enstatite and CB (Bencubbin) chondrites that were formed in the reducing environments of the early solar system. We have also considered distinct planetary accretion scenarios to understand their influence on thermal processing. The majority of our models would require impact‐induced mantle stripping of Mercury by hit and run mechanism with a protoplanet subsequent to its differentiation in order to produce the right size of mantle. However, this can be avoided if we increase the Fe‐Ni‐FeS contents to ~71% by weight. Finally, the models presented here can be used to understand the differentiation of Mercury‐like exoplanets and the planetary embryos of Venus and Earth.     

Magnetism and thermal evolution of the terrestrial planets

Of the terrestrial planets, Earth and probably Mercury possess substantial intrinsic magnetic fields generated by core dynamos, while Venus and Mars apparently lack such fields. Thermal histories are calculated for these planets and are found to admit several possible present states, including those which suggest simple explanations for the observations; whule the cores of Earth and Mercury are continuing to freeze, the cores of Venus and Mars may still be completely liquid. The models assume whole mantle convection, which is parameterized by a simple Nusselt-Rayleigh number relation and dictates the rate at which heat escapes from the core. It is found that completely fluid cores, devoid of intrinsic heat sources, are not likely to sustain thermal convection for the age of the solar system but cool to a subadiabatic, conductive state that can not maintain a dynamo. Planets which nucleate an inner core continue to sustain a dynamo because of the gravitational energy release and chemically driven convection that accompany inner core growth. The absence of a significant inner core can arise in Venus because of its slightly higher temperature and lower central pressure relative to Earth, while a Martian core avoids the onset of freezing if the abundance of sulfur in the core is ?15% by mass. All of the models presented assume that (I) core dynamos are driven by thermal and/or chemical convection; (ii) radiogenic heat production is confined to the mantle; (iii) mantle and core cool from initially hot states which are at the solidus and superliquidus, respectively; and (iv) any inner core excludes the light alloying material (sulfur or oxygen) which then mixes uniformly upward through the outer core. The models include realistic pressure and composition-dependent freezing curves for the core, and material parameters are chosen so that the correct present-day values of heat outflow, upper mantle temperature and viscosity, and inner core radius are obtained for the earth. It is found that Venus and Mars may have once had dynamos maintained by thermal convection alone. Earth may have had a completely fluid core and a dynamo maintained by thermal convection for the first 2 to 3 by, but an inner core nucleates and the dynamo energetics are subsequently dominated by gravitational energy release. Complete freezing of the Mercurian core is prohibited if it contains even a small amount of sulfur, and a dynamo can be maintained by chemical convection in a thin, fluid shell.     

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.     

On the oscillations in Mercury's obliquity

Mercury's capture into the 3:2 spin–orbit resonance can be explained as a result of its chaotic orbital dynamics. One major objective of MESSENGER and BepiColombo spatial missions is to accurately measure Mercury's rotation and its obliquity in order to obtain constraints on internal structure of the planet. Analytical approaches at the first-order level using the Cassini state assumptions give the obliquity constant or quasi-constant. Which is the obliquity's dynamical behavior deriving from a complete spin–orbit motion of Mercury simultaneously integrated with planetary interactions? We have used our SONYR model (acronym of Spin–Orbit N-bodY Relativistic model) integrating the spin–orbit N-body problem applied to the Solar System (Sun and planets). For lack of current accurate observations or ephemerides of Mercury's rotation, and therefore for lack of valid initial conditions for a numerical integration, we have built an original method for finding the libration center of the spin–orbit system and, as a consequence, for avoiding arbitrary amplitudes in librations of the spin–orbit motion as well as in Mercury's obliquity. The method has been carried out in two cases: (1) the spin–orbit motion of Mercury in the 2-body problem case (Sun–Mercury) where an uniform precession of the Keplerian orbital plane is kinematically added at a fixed inclination (S_2K case), (2) the spin–orbit motion of Mercury in the N-body problem case (Sun and planets) (S_n case). We find that the remaining amplitude of the oscillations in the S_n case is one order of magnitude larger than in the S_2K case, namely 4 versus 0.4 arcseconds (peak-to-peak). The mean obliquity is also larger, namely 1.98 versus 1.80 arcminutes, for a difference of 10.8 arcseconds. These theoretical results are in a good agreement with recent radar observations but it is not excluded that it should be possible to push farther the convergence process by drawing nearer still more precisely to the libration center. We note that the dynamically driven spin precession, which occurs when the planetary interactions are included, is more complex than the purely kinematic case. Nevertheless, in such a N-body problem, we find that the 3:2 spin–orbit resonance is really combined to a synchronism where the spin and orbit poles on average precess at the same rate while the orbit inclination and the spin axis orientation on average decrease at the same rate. As a consequence and whether it would turn out that there exists an irreducible minimum of the oscillation amplitude, quasi-periodic oscillations found in Mercury's obliquity should be to geometrically understood as librations related to these synchronisms that both follow a Cassini state. Whatever the open question on the minimal amplitude in the obliquity's oscillations and in spite of the planetary interactions indirectly acting by the solar torque on Mercury's rotation, Mercury remains therefore in a stable equilibrium state that proceeds from a 2-body Cassini state.     

Numerical theory of rotation of the deformable Earth with the two-layer fluid core. Part 1: Mathematical model

Improved differential equations of the rotation of the deformable Earth with the two-layer fluid core are developed. The equations describe both the precession-nutational motion and the axial rotation (i.e. variations of the Universal Time UT). Poincaré’s method of modeling the dynamical effects of the fluid core, and Sasao’s approach for calculating the tidal interaction between the core and mantle in terms of the dynamical Love number are generalized for the case of the two-layer fluid core. Some important perturbations ignored in the currently adopted theory of the Earth’s rotation are considered. In particular, these are the perturbing torques induced by redistribution of the density within the Earth due to the tidal deformations of the Earth and its core (including the effects of the dissipative cross interaction of the lunar tides with the Sun and the solar tides with the Moon). Perturbations of this kind could not be accounted for in the adopted Nutation IAU 2000, in which the tidal variations of the moments of inertia of the mantle and core are the only body tide effects taken into consideration. The equations explicitly depend on the three tidal phase lags δ, δ _c, δ _i responsible for dissipation of energy in the Earth as a whole, and in its external and inner cores, respectively. Apart from the tidal effects, the differential equations account for the non-tidal interaction between the mantle and external core near their boundary. The equations are presented in a simple close form suitable for numerical integration. Such integration has been carried out with subsequent fitting the constructed numerical theory to the VLBI-based Celestial Pole positions and variations of UT for the time span 1984–2005. Details of the fitting are given in the second part of this work presented as a separate paper (Krasinsky and Vasilyev 2006) hereafter referred to as Paper 2. The resulting Weighted Root Mean Square (WRMS) errors of the residuals dθ, sin θd for the angles of nutation θ and precession are 0.136 mas and 0.129 mas, respectively. They are significantly less than the corresponding values 0.172 and 0.165 mas for IAU 2000 theory. The WRMS error of the UT residuals is 18 ms.     

A Postulated Planetary Collision, the Terrestrial Planets, the Moon and Smaller Solar-System Bodies

In a scenario produced by the Capture Theory of planetary formation, a collision between erstwhile solar-system giant planets, of masses 798.75 and 598.37 M _⊕, is simulated using smoothed-particle hydrodynamics. Due to grain-surface chemistry that takes place in star-forming clouds, molecular species containing hydrogen, with a high D/H ratio taken as 0.01, form a layer around each planetary core. Temperatures generated by the collision initiate D–D reactions in these layers that, in their turn, trigger a reaction chain involving heavier elements. The nuclear explosion shatters and disperses both planets, leaving iron-plus-silicate stable residues identified as a proto-Venus and proto-Earth. A satellite of one of the colliding planets, captured or retained by the proto-Earth core, gave the Moon; two massive satellites released into heliocentric orbits became Mercury and Mars. For the Moon and Mars, abrasion of their surfaces exposed to collision debris results in hemispherical asymmetry. Mercury, having lost a large part of its mantle due to massive abrasion, reformed to give the present high-density body. Debris from the collision gave rise to asteroids and comets, much of the latter forming an inner reservoir stretching outwards from the inner Kuiper Belt that replenishes the Oort Cloud when it is depleted by a severe perturbation. Other features resulting from the outcome of the planetary collision are the relationship of Pluto and Triton to Neptune, the presence of dwarf planets and light-atom isotopic anomalies in meteorites.     

Search for sulfates on the surface of Ceres

The formation of hydrated salts is an expected consequence of aqueous alteration of Main Belt objects, particularly for large, volatile‐rich protoplanets like Ceres. Sulfates, present on water‐bearing planetary bodies (e.g., Earth, Mars, and carbonaceous chondrite parent bodies) across the inner solar system, may contribute to Ceres’ UV and IR spectral signature along with phyllosilicates and carbonates. We investigate the presence and stability of hydrated sulfates under Ceres’ cryogenic, low‐pressure environment and the consequent spectral effects, using UV–Vis–IR reflectance spectroscopy. H₂O loss begins instantaneously with vacuum exposure, measured by the attenuation of spectral water absorption bands, and a phase transition from crystalline to amorphous is observed for MgSO₄·6H₂O by X‐ray powder diffraction. Long‐term (>40 h), continuous exposure of MgSO₄·nH₂O (n = 0, 6, 7) to low pressure (10^?3–10^?6 Torr) causes material decomposition and strong UV absorption below 0.5 μm. Our measurements suggest that MgSO₄·6H₂O grains (45–83 μm) dehydrate to 2% of the original 1.9 μm water band area over ~0.3 Ma at 200 K on Ceres and after ~42 Ma for 147 K. These rates, inferred from an Avrami dehydration model, preclude MgSO₄·6H₂O as a component of Ceres’ surface, although anhydrous and minimally hydrated sulfates may be present. A comparison between Ceres emissivity spectra and laboratory reflectance measurements over the infrared range (5–17 μm) suggests sulfates cannot be excluded from Ceres’ mineralogy.     

Clues on the importance of comets in the origin and evolution of the atmospheres of Titan and Earth

Earth and Titan are two planetary bodies formed far from each other. Nevertheless the chemical composition of their atmospheres exhibits common indications of being produced by the accretion, plus ulterior in-situ processing of cometary materials. This is remarkable because while the Earth formed in the inner part of the disk, presumably from the accretion of rocky planetesimals depleted in oxygen and exhibiting a chemical similitude with enstatite chondrites, Titan formed within Saturn's sub-nebula from oxygen- and volatile-rich bodies, called cometesimals. From a cosmochemical and astrobiological perspective, the study of the H, C, N, and O isotopes on Earth and Titan could be the key to decipher the processes occurred in the early stages of formation of both planetary bodies. The main goal of this paper is to quantify the presumable ways of chemical evolution of both planetary bodies, in particular the abundance of CO and N₂ in their early atmospheres. In order to do that the primeval atmospheres and evolution of Titan and Earth have been analyzed from a thermodynamic point of view. The most relevant chemical reactions involving these species and presumably important at their early stages are discussed. Then, we have interpreted the results of this study in light of the results obtained by the Cassini–Huygens mission on these species and their isotopes. Given that H, C, N, and O were preferentially depleted from inner disk materials that formed our planet, the observed similitude of their isotopic fractionation, and subsequent close evolution of Earth's and Titan's atmospheres points towards a cometary origin of Earth atmosphere. Consequently, our scenario also supports the key role of late veneers (comets and water-rich carbonaceous asteroids) enriching the volatile content of the Earth at the time of the late heavy bombardment of terrestrial planets.     

Oxidation of CO on surface hematite in high CO2 atmospheres

We propose a mechanism for the oxidation of gaseous CO into CO₂ occurring on the surface mineral hematite (Fe₂O₃(s)) in hot, CO₂-rich planetary atmospheres, such as Venus. This mechanism is likely to constitute an important source of tropospheric CO₂ on Venus and could at least partly address the CO₂ stability problem in Venus’ stratosphere, since our results suggest that atmospheric CO₂ is produced from CO oxidation via surface hematite at a rate of 0.4 petagrammes (Pg) CO₂ per (Earth) year on Venus which is about 45% of the mass loss of CO₂ via photolysis in the Venusian stratosphere. We also investigated CO oxidation via the hematite mechanism for a range of planetary scenarios and found that modern Earth and Mars are probably too cold for the mechanism to be important because the rate-limiting step, involving CO(g) reacting onto the hematite surface, proceeds much slower at lower temperatures. The mechanism may feature on extrasolar planets such as Gliese 581c or CoRoT-7b assuming they can maintain solid surface hematite which, e.g. starts to melt above about 1200 K. The mechanism may also be important for hot Hadean-type environments and for the emerging class of hot Super-Earths with planetary surface temperatures between about 600 and 900 K.     

The abundance and isotopic composition of water in eucrites

Volatile elements play a key role in the dynamics of planetary evolution. Extensive work has been carried out to determine the abundance, distribution, and source(s) of volatiles in planetary bodies such as the Earth, Moon, and Mars. A recent study showed that the water in apatite from eucrites has similar hydrogen isotopic compositions compared to water in terrestrial rocks and carbonaceous chondrites, suggesting that water accreted very early in the inner solar system given the ancient crystallization ages (~4.5 Ga) of eucrites. Here, the measurements of water (reported as equivalent H₂O abundances) and the hydrogen isotopic composition (δD) of apatite from five basaltic eucrites and one cumulate eucrite are reported. Apatite H₂O abundances range from ~30 to ~3500 ppm and are associated with a weighted average δD value of ?34 ± 67‰. No systematic variations or correlations are observed in H₂O abundance or δD value with eucrite geochemical trend or metamorphic grade. These results extend the range of previously published hydrogen isotope data for eucrites and confirm the striking homogeneity in the H‐isotopic composition of water in eucrites, which is consistent with a common source for water in the inner solar system.     

Stability field and phase transition pathways of hydrous ferric sulfates in the temperature range 50 °C to 5 °C: Implication for martian ferric sulfates

We report the results from a systematic laboratory investigation on the fundamental properties of hydrous ferric sulfates. The study involves 150 experiments with duration of over 4 years on the stability field and phase transition pathways under Mars relevant environmental conditions for five ferric sulfates: ferricopiapite [Fe_4.67(SO₄)₆(OH)₂·20H₂O], kornelite [Fe₂(SO₄)₃·7H₂O], a crystalline and an amorphous pentahydrated ferric sulfate [Fe₂(SO₄)₃·5H₂O], and rhomboclase [FeH(SO₄)₂·4H₂O]. During the processes of phase transitions, we observed the phenomena that reflect fundamental properties of these species and the occurrence of other common hydrous ferric sulfates, e.g. paracoquimbite [Fe₂(SO₄)₃·9H₂O]. Based on the results of this set of experiments, we have drown the boundaries of deliquescence zone of five hydrous ferric sulfates and estimated the regions of their stability field in temperature (T) – relative humidity (RH) space. Furthermore, we selected the experimental parameters for a next step investigation, which is to determine the location of the phase boundary between two solid ferric sulfates, kornelite [Fe₂(SO₄)₃·7H₂O] and pentahydrated ferric sulfate [Fe₂(SO₄)₃·5H₂O]. The experimental observations in ferricopiapite dehydration processes were used to interpret the observed spectral change of Fe-sulfate-rich subsurface soils on Mars after their exposure by the Spirit rover to current martian atmospheric conditions.     

Planetary core size: A seismological approach

A key parameter for understanding the geodynamics of a terrestrial planet is the size of its core. Numerical evaluation of 28 different interior structure models of Mercury, Venus, Earth, the Moon, and Mars suggests that there is an almost linear relationship between the core radius and the extent of the seismic P-wave core shadow. A scaling law is derived from a simple mantle density and velocity model that permits the interpretation of respective seismic measurements on terrestrial planetary bodies.