Preliminary studies concerning subsurface probes for the exploration of icy planetary bodies

It is well-known that the permanent terrestrial ice sheets (glaciers and polar caps) contain a lot of information about the recent geological history and in particular about climatic changes. Extrapolating this fact to other ice sheets in the solar system (e.g. the Mars polar regions, the icy moons of the outer planets, etc.), we may expect a similar wealth of information. To obtain this information it is possible to drill holes or melt the ice by a heated probe, which in this way is able to penetrate the surface and investigate the deeper layers in situ. In the latter case the driving agent is the heating power and the weight of the probe. In this paper we consider the application of such “melting probes” for exploring the structure of ice sheets in extraterrestrial environments. We describe several laboratory experiments with simple melting probes performed under cryo-vacuum conditions and compare the results with tests in a terrestrial environment. The experiments revealed that under space conditions the downward motion of a heated probe in an ice sheet is characterized by intermittent periods of sublimation and melting of the surrounding ice, sometimes interrupted by periods where a part of the probe's outer surface is frozen to the surrounding ice. This leads to a temporary blocking of the probe's downward motion. A similar situation can occur when the trailing tether is frozen in behind the probe. During the periods of ice sublimation the penetration process is significantly more power consuming, due to the large difference between the latent heat of sublimation and the latent heat of melting for water ice.     

Huygens HASI servo accelerometer: A review and lessons learned

The servo accelerometer constituted a vital part of the Huygens Atmospheric Structure Instrument (HASI): flown aboard the Huygens probe, it operated successfully during the probe's entry, descent, and landing on Titan, on 14th January 2005. This paper reviews the Servo accelerometer, starting from its development/assembly in the mid-1990s, to monitoring its technical performance through its seven-year long in-flight (or cruise) journey, and finally its performance in measuring acceleration (or deceleration) upon encountering Titan's atmosphere.The aim of this article is to review the design, ground tests, in-flight tests and operational performance of the Huygens servo accelerometer. Techniques used for data analysis and lessons learned that may be useful for accelerometry payloads on future planetary missions are also addressed.The main finding of this review is that the conventional approach of having multiple channels to cover a very broad measurement range: from 10⁻⁶g to the order of 10g (where g=Earth's surface gravity, 9.8 m/s²), with on-board software deciding which of the channels to telemeter depending on the magnitude of the measured acceleration, works well. However, improvements in understanding the potential effects of the sensor drifts and ageing on the measurements can be achieved in future missions by monitoring the ‘scale factor’—a measure of such sensors’ sensitivity, along with the already implemented monitoring of the sensor's offset during the in-flight phase.     

Numerical modeling of protocore destabilization during planetary accretion: Methodology and results

We developed and tested an efficient 2D numerical methodology for modeling gravitational redistribution processes in a quasi spherical planetary body based on a simple Cartesian grid. This methodology allows one to implement large viscosity contrasts and to handle properly a free surface and self-gravitation. With this novel method we investigated in a simplified way the evolution of gravitationally unstable global three-layer structures in the interiors of large metal-silicate planetary bodies like those suggested by previous models of cold accretion [Sasaki, S., Nakazawa, K., 1986. J. Geophys. Res. 91, 9231-9238; Karato, S., Murthy, V.R., 1997. Phys. Earth Planet Interios 100, 61-79; Senshu, H., Kuramoto, K., Matsui, T., 2002. J. Geophys. Res. 107 (E12), 5118. 10.1029/2001JE001819]: an innermost solid protocore (either undifferentiated or partly differentiated), an intermediate metal-rich layer (either continuous or disrupted), and an outermost silicate-rich layer. Long-wavelength (degree-one) instability of this three-layer structure may strongly contribute to core formation dynamics by triggering planetary-scale gravitational redistribution processes. We studied possible geometrical modes of the resulting planetary reshaping using scaled 2D numerical experiments for self-gravitating planetary bodies with Mercury-, Mars- and Earth-size. In our simplified model the viscosity of each material remains constant during the experiment and rheological effects of gravitational energy dissipation are not taken into account. However, in contrast to a previously conducted numerical study [Honda, R., Mizutani, H., Yamamoto, T., 1993. J. Geophys. Res. 98, 2075-2089] we explored a freely deformable planetary surface and a broad range of viscosity ratios between the metallic layer and the protocore (0.001-1000) as well as between the silicate layer and the protocore (0.001-1000). An important new prediction from our study is that realistic modes of planetary reshaping characterized by a high viscosity protocore and low viscosity molten silicate and metal [Senshu, H., Kuramoto, K., Matsui, T., 2002. J. Geophys. Res. 107 (E12), 5118. 10.1029/2001JE001819] may result in the transient exposure of the protocore to the planetary surface and a strongly (up to 8% of the planetary diameter) aspherical deviation of the planetary shape during the early stages of core formation. Exposure of the protocore might happen in the early stages of iron core formation. This process may conceivably convert a large amount of potential energy into temperature increase and a transient strongly non-uniform depth of the magma ocean around the protoplanet. Our simplified model also predicts that the time for metallic core formation out of the metal-rich layer depends mainly on the dynamics of the deformation of the solid strong protocore. In nature this dynamics will be strongly dependent on the effective viscosity of the protocore, which should generally have non-Newtonian pressure-, temperature-, and stress-dependent rheology with strong thermomechanical feedbacks from gravitational energy dissipation.     

Spin rates of small moonlets embedded in planetary rings: I. Three-body calculations

We investigate the spin rates of moonlets embedded in planetary rings, subject to collisions with surrounding small particles, using three-body integrations including friction and spins. All successive impacts of the particle with the moonlet are followed, including a possible sliding phase after the initial inelastic rebounds. Two methods for treating impacts, (1) as instantaneous velocity changes and (2) using an impact force model, are applied after Salo (1995, Icarus 117, 287). Conducting a series of integrations with various initial summed spin velocity of the moonlet and the particle, we determine the equilibrium spin rate for which the averaged torque vanishes. This equilibrium spin rate corresponds to the final spin rate of the moonlet if the moonlet is much larger than the surrounding particles; it also corresponds to the mean spin rate for a ring composed of identical particles. We find that the equilibrium spin rate is enhanced by sliding orbits as compared with the spin rate determined by considering only the first impacts of the particles with the moonlet. If the random velocities of incident particles are small enough, the resulting equilibrium spin rate of the moonlet can be larger than the synchronous rotation rate, for r_p∼1, where r_p denotes the sum of radii of the colliding pair normalized by their mutual Hill radius. In this special case aggregates without internal strength may become rotationally unstable. However, the equilibrium spin rate decreases with increasing random velocity, and aggregates are always rotationally stable in the more likely case where the relative velocities are comparable to the mutual escape velocity. We also compare our results with the mean spin rates found in previous N-body simulations, and find a good agreement for optically thin rings; however the spin rates for optically thick rings are significantly larger than those predicted by our three-body calculations.     

Melting and Sublimation of Planetary Ices Under Low Pressure Conditions: Laboratory Experiments with a Melting Probe Prototype

One possibility to explore the subsurface layers of icy bodies is to use a probe with a “hot tip", which is able to penetrate ice layers by melting. Such probes have been built and used in the past for the exploration of terrestrial polar ice sheets and may also become useful tools to explore other icy layers in the Solar System. Examples for such layers are the polar areas of Mars or the icy crust of Jupiter’s moon Europa. However, while on Earth a heated probe launched into an ice sheet always causes melting with subsequent refreezing, the behaviour of such a probe in a low pressure environment is quite different. We report on the results of some experiments with a simple “melting probe" prototype with two different kinds of hot tips in a vacuum environment. For one of the tips the probe moved into two types of ice samples: (i) compact water ice and (ii) porous water ice with a snow (firn) like texture. It was also found that the penetration behaviour was basically different for the two sample types even when the same kind of tip was used. While in the porous sample the ice was only subliming, the phase changes occurring during the interaction of the tip with the compact ice are much more complex. Here alternating phases of melting and sublimation occur. The absence of the liquid phase has severe consequences on the performance of a “melting probe" under vacuum conditions: In this environment we find a high thermal resistance between the probe surface and the underlying ice. Therefore, only a low percentage of the heat that is generated in the tip is used to melt or sublime the ice, the bulk of the power is transferred towards the rear end of the probe. This is particularly a problem in the initial phases of an ice penetration experiment, when the probe has not yet penetrated the ice over its whole length. In the compact ice sample, phases could be observed, where a high enough gas pressure had built up locally underneath the probe, so that melting becomes possible. Only during these melting periods the thermal contact between the probe and the ice is good and in consequence the melting probe works effectively.     

Moonlet induced wakes in planetary rings: Analytical model including eccentric orbits of moon and ring particles

Saturn’s rings host two known moons, Pan and Daphnis, which are massive enough to clear circumferential gaps in the ring around their orbits. Both moons create wake patterns at the gap edges by gravitational deflection of the ring material (Cuzzi, J.N., Scargle, J.D. [1985]. Astrophys. J. 292, 276-290; Showalter, M.R., Cuzzi, J.N., Marouf, E.A., Esposito, L.W. [1986]. Icarus 66, 297-323). New Cassini observations revealed that these wavy edges deviate from the sinusoidal waveform, which one would expect from a theory that assumes a circular orbit of the perturbing moon and neglects particle interactions. Resonant perturbations of the edges by moons outside the ring system, as well as an eccentric orbit of the embedded moon, may partly explain this behavior (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S., Burns, J.A., Dones, L. [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S. [2009]. Astron. J. 138, 272-286). Here we present an extended non-collisional streamline model which accounts for both effects. We describe the resulting variations of the density structure and the modification of the nonlinearity parameter q. Furthermore, an estimate is given for the applicability of the model. We use the streamwire model introduced by Stewart (Stewart, G.R. [1991]. Icarus 94, 436-450) to plot the perturbed ring density at the gap edges.We apply our model to the Keeler gap edges undulated by Daphnis and to a faint ringlet in the Encke gap close to the orbit of Pan. The modulations of the latter ringlet, induced by the perturbations of Pan (Burns, J.A., Hedman, M.M., Tiscareno, M.S., Nicholson, P.D., Streetman, B.J., Colwell, J.E., Showalter, M.R., Murray, C.D., Cuzzi, J.N., Porco, C.C., and the Cassini ISS team [2005]. Bull. Am. Astron. Soc. 37, 766), can be well described by our analytical model. Our analysis yields a Hill radius of Pan of 17.5 km, which is 9% smaller than the value presented by Porco (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236), but fits well to the radial semi-axis of Pan of 17.4 km. This supports the idea that Pan has filled its Hill sphere with accreted material (Porco, C.C., Thomas, P.C., Weiss, J.W., Richardson, D.C. [2007]. Science 318, 1602-1607). A numerical solution of a streamline is used to estimate the parameters of the Daphnis-Keeler gap system, since the close proximity of the gap edge to the moon induces strong perturbations, not allowing an application of the analytic streamline model. We obtain a Hill radius of 5.1 km for Daphnis, an inner edge variation of 8 km, and an eccentricity for Daphnis of 1.5 × 10⁻⁵. The latter two quantities deviate by a factor of two from values gained by direct observations (Jacobson, R.A., Spitale, J., Porco, C.C., Beurle, K., Cooper, N.J., Evans, M.W., Murray, C.D. [2008]. Astron. J. 135, 261-263; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767), which might be attributed to the neglect of particle interactions and vertical motion in our model.     

Stability of outer planetary orbits around binary stars: A comparison of Hill's and Laplace's stability criteria

A comparison is made between the stability criteria of Hill and that of Laplace to determine the stability of outer planetary orbits encircling binary stars. The restricted, analytically determined results of Hill's method by Szebehely and co-workers and the general, numerically integrated results of Laplace's method by Graziani and Black are compared for varying values of the mass parameter =m ₂/(m ₁+m ₂). For 00.15, the closest orbit (lower limit of radius) an outer planet in a binary system can have and still remain stable is determined by Hill's stability criterion. For >0.15, the critical radius is determined by Laplace's stability criterion. It appears that the Graziani-Black stability criterion describes the critical orbit within a few percent for all values of .     

The full set of gas giant structures: I: On the origin of planetary masses and the planetary initial mass function

Context

Current planet search programs are detecting extrasolar planets at a rate of 60 planets per year. These planets show more diverse properties than was expected.

Aims

We try to get an overview of possible gas giant (proto-) planets for a full range of orbital periods and stellar masses. This allows the prediction of the full range of possible planetary properties which might be discovered in the near future.

Methods

We calculate the purely hydrostatic structure of the envelopes of proto-planets that are embedded in protoplanetary disks for all conceivable locations: combinations of different planetesimal accretion rates, host star masses, and orbital separations. At each location all hydrostatic equilibrium solutions to the planetary structure equations are determined by variation of core mass and pressure over many orders of magnitude. For each location we analyze the distribution of planetary masses.

Results

We get a wide spectrum of core-envelope structures. However, practically all calculated proto-planets are in the planetary mass range. Furthermore, the planet masses show a characteristic bimodal, sometimes trimodal, distribution. For the first time, we identify three physical processes that are responsible for the three characteristic planet masses: self-gravity in the Hill sphere, compact objects, and a region of very low adiabatic pressure gradient in the hydrogen equation of state. Using these processes, we can explain the dependence of the characteristic masses on the planet’s location: orbital period, host star mass, and planetesimal accretion rate (luminosity). The characteristic mass caused by the self-gravity effect at close proximity to the host star is typically one Neptune mass, thus producing the so-called hot Neptunes.

Conclusions

Our results suggest that hot Jupiters with orbital period less than 64 days (the exact location of the boundary depends on stellar type and accretion rate) have quite distinct properties which we expect to be reflected in a different mass distribution of these planets when compared to the “normal” planetary population. We use our theoretical survey to produce an upper mass limit for embedded planets: the maximum embedded equilibrium mass (MEEM). This naturally explains the lack of high mass planets between 3 and 64 days orbital period.     

Transformation of Trojans into quasi-satellites during planetary migration and their subsequent close-encounters with the host planet

We use numerical integrations to investigate the dynamical evolution of resonant Trojan and quasi-satellite companions during the late stages of migration of the giant planets Jupiter, Saturn, Uranus, and Neptune. Our migration simulations begin with Jupiter and Saturn on orbits already well separated from their mutual 2:1 mean-motion resonance. Neptune and Uranus are decoupled from each other and have orbital eccentricities damped to near their current values. From this point we adopt a planet migration model in which the migration speed decreases exponentially with a characteristic timescale τ (the e-folding time). We perform a series of numerical simulations, each involving the migrating giant planets plus test particle Trojans and quasi-satellites. We find that the libration frequencies of Trojans are similar to those of quasi-satellites. This similarity enables a dynamical exchange of objects back and forth between the Trojan and quasi-satellite resonances during planetary migration. This exchange is facilitated by secondary resonances that arise whenever there is more than one migrating planet. For example, secondary resonances may occur when the circulation frequencies, f, of critical arguments for the Uranus-Neptune 2:1 mean-motion near-resonance are commensurate with harmonics of the libration frequency of the critical argument for the Trojan and quasi-satellite 1:1 mean-motion resonance . Furthermore, under the influence of these secondary resonances quasi-satellites can have their libration amplitudes enlarged until they undergo a close-encounter with their host planet and escape from the resonance. High-resolution simulations of this escape process reveal that ≈80% of jovian quasi-satellites experience one or more close-encounters within Jupiter’s Hill radius (R_H) as they are forced out of the quasi-satellite resonance. As many as ≈20% come within R_H/4 and ≈2.5% come within R_H/10. Close-encounters of escaping quasi-satellites occur near or even below the 2-body escape velocity from the host planet. Finally, the exchange and escape of Trojans and quasi-satellites continues to as late as 6-9τ in some simulations. By this time the dynamical evolution of the planets is strongly dominated by distant gravitational perturbations between the planets rather than the migration force. This suggests that exchange and escape of Trojans and quasi-satellites may be a contemporary process associated with the present-day near-resonant configuration of some of the giant planets in our Solar System.     

Highly reducing conditions during core formation on Mercury: Implications for internal structure and the origin of a magnetic field

The high average density and low surface FeO content of the planet Mercury are shown to be consistent with very low oxygen fugacity during core segregation, in the range 3-6 log units below the iron-wüstite buffer. These low oxygen fugacities, and associated high metal content, are characteristic of high-iron enstatite (EH) and Bencubbinite (CB) chondrites, raising the possibility that such materials may have been important building blocks for this planet. With this idea in mind we have explored the internal structure of a Mercury sized planet of EH or CB bulk composition. Phase equilibria in the silicate mantle have been modeled using the thermodynamic calculator p-MELTS, and these simulations suggest that orthopyroxene will be the dominant mantle phase for both EH and CB compositions, with crystalline SiO₂ being an important minor phase at all pressures. Simulations for both compositions predict a plagioclase-bearing “crust” at low pressure, significant clinopyroxene also being calculated for the CB bulk composition. Concerning the core, comparison with recent high pressure and high temperature experiments relevant to the formation of enstatite meteorites, suggest that the core of Mercury may contain several wt.% silicon, in addition to sulfur. In light of the pressure of the core-mantle boundary on Mercury (∼7 GPa) and the pressure at which the immiscibility gap in the system Fe-S-Si closes (∼15 GPa) we suggest that Mercury’s core may have a complex shell structure comprising: (i) an outer layer of Fe-S liquid, poor in Si; (ii) a middle layer of Fe-Si liquid, poor in S; and (iii) an inner core of solid metal. The distribution of heat-producing elements between mantle and core, and within a layered core have been quantified. Available data for Th and K suggest that these elements will not enter the core in significant amounts. On the other hand, for the case of U both recently published metal/silicate partitioning data, as well as observations of U distribution in enstatite chondrites, suggest that this element behaves as a chalcophile element at low oxygen fugacity. Using these new data we predict that U will be concentrated in the outer layer of the mercurian core. Heat from the decay of U could thus act to maintain this part of Mercury’s core molten, potentially contributing to the origin of Mercury’s magnetic field. This result contrasts with the Earth where the radioactive decay of U represents a negligible contribution to core heating.     

Global contraction of planetary bodies due to despinning: Application to Mercury and Iapetus

We extend previous work on the global tectonic patterns generated by despinning with a self-consistent treatment of the isotropic despinning contraction that has been ignored. We provide simple analytic approximations that quantify the effect of the isotropic despinning contraction on the global shape and tectonic pattern. The isotropic despinning contraction of Mercury is ∼93 m (T/1 day)⁻², where T is the initial rotation period. If we take into account both the isotropic contraction and the degree-2 deformations associated with despinning, the preponderance of compressional tectonic features on Mercury’s surface requires an additional isotropic contraction ?1 km (T/1 day)⁻², presumably due to cooling of the interior and growth of the solid inner core. The isotropic despinning contraction of Iapetus is ∼9 m (T/16 h)⁻², and it is not sensitive to the presence of a core or the thickness of the elastic lithosphere. The tectonic pattern expected for despinning, including the isotropic contraction, does not explain Iapetus’ ridge. Furthermore, the ridge remains unexplained with the addition of any isotropic compressional stresses, including those generating by cooling.     

The deep interior of Venus, Mars, and the Earth: A brief review and the need for planetary surface-based measurements

A comparison of the internal structure of Earth-like planets is unavoidable to understand the formation and evolution of the solar system, and the differences between Earth’s, Mars’, and Venus’ atmospheres, surfaces and tectonic behaviors. Recent studies point at the role of core structure and dynamics in the evolution of the atmosphere, mantle and crust. On Earth, the crust thickness and the radius and physical state of the cores are known for almost one century, since the advent of seismological observations, but the lack of long-term surface-based geodetic, electromagnetic and seismological observations on the other planets, results in very large uncertainties on the crust thickness, on the temperature and composition of their mantle, and on the size and physical state of their cores. According to the currently available geodetic data, Mars’ dimensionless mean-moment-of-inertia ratio is equal to 0.3653±0.0008. When combined with geochemical observations and with the inputs of laboratory experiments on planetary materials at high pressure and high temperature, this result constrains a narrow range of density values for Mars’ mantle and favors a light [6200-6765 kg m⁻³] sulfur-rich core, but it still allows for a 1600-1750 km range for the core radius, i.e. an uncertainty at least ten times larger than the precision obtained in 1913 by Gutenberg for the Earth’s core. Mars’ mantle density distribution may be explained by a large range of temperatures and mineralogical compositions, either olivine- or pyroxene-rich. The unknown mean thickness of Mars’ crust makes necessary a number of working assumptions for the interpretation of gravimetric and magnetic data. The situation is worse for Venus, and the most conservative model of its deep interior is a transposition of the Earth’s structure scaled to Venus’ radius and mass. The temperature conditions at the surface of this planet hardly make possible long-term ground-based measurements, but this is indeed feasible at the surface of Mars. Precise measurements of Mars’ crust thickness, core radius and structure, and the proof of the existence or absence of an inner core, would put tight constraints on mantle dynamics and thermal evolution, and on possible scenarios leading to the extinction of Mars’ magnetic field about 4.0 Ga ago. Long-lasting surface-based geodetic, seismological and magnetic observations would provide this information, as well as the distributions as a function of depth of the density, elastic and anelastic parameters, and electrical conductivity. Current studies on the structure of Earth’s deep interior demonstrate that the latter data set, when constrained by laboratory experiments, may be inverted in terms of temperature, chemical, and mineralogical compositions.     

Formation of gas giant planets: core accretion models with fragmentation and planetary envelope

We have calculated formation of gas giant planets based on the standard core accretion model including effects of fragmentation and planetary envelope. The accretion process is found to proceed as follows. As a result of runaway growth of planetesimals with initial radii of ∼10 km, planetary embryos with a mass of ∼10²⁷ g (∼ Mars mass) are found to form in ∼10⁵ years at Jupiter's position (5.2 AU), assuming a large enough value of the surface density of solid material (25 g/cm²) in the accretion disk at that distance. Strong gravitational perturbations between the runaway planetary embryos and the remaining planetesimals cause the random velocities of the planetesimals to become large enough for collisions between small planetesimals to lead to their catastrophic disruption. This produces a large number of fragments. At the same time, the planetary embryos have envelopes, that reduce energies of fragments by gas drag and capture them. The large radius of the envelope increases the collision rate between them, resulting in rapid growth of the planetary embryos. By the combined effects of fragmentation and planetary envelope, the largest planetary embryo with 21M_⊕ forms at 5.2 AU in 3.8×10⁶ years. The planetary embryo is massive enough to start a rapid gas accretion and forms a gas giant planet.     

Gamma–hadron separation methods for the VERITAS array of four imaging atmospheric Cherenkov telescopes

Ground-based arrays of imaging atmospheric Cherenkov telescopes have emerged as the most sensitive γ-ray detectors in the energy range of about 100 GeV and above. The strengths of these arrays are a very large effective collection area on the order of 10⁵ m², combined with excellent single photon angular and energy resolutions. The sensitivity of such detectors is limited by statistical fluctuations in the number of Cosmic-ray initiated air showers that resemble γ-ray air showers in many ways. In this paper, we study the performance of simple event reconstruction methods when applied to simulated data of the Very Energetic Radiation Imaging Telescope Array System (VERITAS) experiment. We review methods for reconstructing the arrival direction and the energy of the primary photons, and examine means to improve on their performance. For a software threshold energy of 300 GeV (100 GeV), the methods achieve point source angular and energy resolutions of σ_63% = 0.1° (0.2°) and σ_68% = 15% (22%), respectively. The main emphasis of the paper is the discussion of γ–hadron separation methods for the VERITAS experiment. We find that the information from several methods can be combined based on a likelihood ratio approach and the resulting algorithm achieves a γ–hadron suppression with a quality factor that is substantially higher than that achieved with the standard methods used so far.     

Noble gases in planetary atmospheres: Implications for the origin and evolution of atmospheres

The radiogenic and primordial noble gas content of the atmospheres of Venus, Earth, and Mars are compared with one another and with the noble gas content of other extraterrestial samples, especially meteorites. The fourfold depletion of ⁴⁰Ar for Venus relative to the Earth is attributed to the outgassing rates and associated tectonics and volcanic styles for the two planets diverging significantly within the first billion or so years of their history, with the outgassing rate for Venus becoming much less than that for the Earth at subsequent times. This early divergence in the tectonic style of the two planets may be due to a corresponding early onset of the runaway greenhouse on Venus. The 16-fold depletion of ⁴⁰Ar for Mars relative to the Earth may be due to a combination of a mild K depletion for Mars, a smaller fraction of its interior being outgassed, and to an early reduction in its outgassing rate. Venus has lost virtually all of its primordial He and some of its radiogenic He. The escape flux of He may have been quite substantial in Venus' early history, but much diminished at later times, with this time variation being perhaps strongly influenced by massive losses of H₂ resulting from efficient H₂O loss processes.Key trends in the primordial noble gas content of terrestial planetary atmospheres include (1) a several orders of magnitude decrease in ²⁰Ne and ³⁶Ar from Venus to Earth to Mars; (2) a nearly constant ²⁰Ne/³⁶Ar ratio which is comparable to that found in the more primitive carbonaceous chondrites and which is two orders of magnitude smaller than the solar ratio; (3) a sizable fractionation of Ar, Kr, and Xe from their solar ratios, although the degree of fractionation, especially for ³⁶Ar/¹³²Xe, seems to decrease systematically from carbonaceous chondrites to Mars to Earth to Venus; and (4) large differences in Ne and Xe isotopic ratios among Earth, meteorites, and the Sun. Explaining trends (2), (2) and (4), and (1) pose the biggest problems for the solar-wind implantation, primitive atmosphere, and late veneer hypotheses, respectively. It is suggested that the grain-accretion hypothesis can explain all four trends, although the assumptions needed to achieve this agreement are far from proven. In particular, trends (1), (2), (3), and (4) are attributed to large pressure but small temperature differences in various regions of the inner solar system at the times of noble gas incorporation by host phases; similar proportions of the host phases that incorporated most of the He and Ne on the one hand (X) and Ar, Kr, and Xe on the other hand (Q); a decrease in the degree of fractionation with increasing noble-gas partial pressure; and the presence of interstellar carriers containing isotopically anomalous noble gases.Our analysis also suggests that primordial noble gases were incorporated throughout the interior of the outer terrestial planets, i.e., homogeneous accretion is favored over inhomogeneous accretion. In accord with meteorite data, we propose that carbonaceous materials were key hosts for the primordial noble gases incorporated into planets and that they provided a major source of the planets' CO₂ and N₂.     

Tungsten isotopes in ferroan anorthosites: Implications for the age of the Moon and lifetime of its magma ocean

New W isotope data for ferroan anorthosites 60025 and 62255 and low-Ti mare basalt 15555 show that these samples, contrary to previous reports [Lee, D.C., et al., 1997. Science 278, 1098-1103; Lee, D.C., et al., 2002. Earth Planet. Sci. Lett. 198, 267-274], have a W isotope composition that is indistinguishable from KREEP and other mare basalts. This requires crust extraction on the Moon later than ∼60 Myr after CAI formation, consistent with ¹⁴⁷Sm-¹⁴³Nd ages for ferroan anorthosites. The absence of ¹⁸²Hf-induced ¹⁸²W variations in the Moon is consistent with formation of the Moon at after CAI formation that has been inferred based on the indistinguishable ¹⁸²W/¹⁸⁴W ratios of the bulk Moon and the bulk silicate Earth. The uncertainties on the age of the Moon and the age of the oldest lunar samples currently hamper a precise determination of the duration of magma ocean solidification and are consistent with both an almost immediate crystallization and a more protracted timescale of ∼100 Myr.     

Basaltic glass as a habitat for microbial life: Implications for astrobiology and planetary exploration

Recent studies have demonstrated that terrestrial subaqueous basalts and hyaloclastites are suitable microbial habitats. During subaqueous basaltic volcanism, glass is produced by the quenching of basaltic magma upon contact with water. On Earth, microbes rapidly begin colonizing the glassy surfaces along fractures and cracks that have been exposed to water. Microbial colonization of basaltic glass leads to the alteration and modification of the rocks and produces characteristic granular and/or tubular bioalteration textures. Infilling of the alteration textures by minerals such as phyllosilicates, zeolites and titanite may enable their preservation through geologic time. Basaltic rocks are a major component of the Martian crust and are widespread on other solar system bodies. A variety of lines of evidence strongly suggests the long-term existence of abundant liquid water on ancient Mars. Recent orbiter, lander and rover missions have found evidence for the presence of transient liquid water on Mars, perhaps persisting to the present day. Many other solar system bodies, notably Europa, Enceladus and other icy satellites, may contain (or have once hosted) subaqueous basaltic glasses. The record of terrestrial glass bioalteration has been interpreted to extend as far back as ∼3.5 billion years ago and is widespread in oceanic crust and its metamorphic equivalents. The terrestrial record of glass bioalteration strongly suggests that glassy or formerly glassy basaltic rocks on extraterrestrial bodies that have interacted with liquid water are high-value targets for astrobiological exploration.     

The frequency-area distribution of volcanic units on Venus: Implications for planetary resurfacing

The areas of volcanic units on Venus have been measured on the 1:5000000 geological maps published by NASA/USGS. These data were used to obtain a frequency-area distribution. The cumulative frequency-area distribution of 1544 specific occurrence of units cover six orders of magnitude from the largest unit (30 × 10⁶ km²) to the smallest (20 km²). The probability distribution function has been calculated. The medium and large volcanic units correlate well with a power-law (fractal) relation for the dependence of frequency on area with a slope of −1.83. There are fewer small units than the expected values provided by the power-law relation. Our measurements cover 21.02% of the planetary surface, 3.59% of the study area was found to be tessera terrain and is excluded from this study of volcanism. The measurements were restricted to areas where geological maps have been published. The analysis was performed on two independent areas of the planet, with a complete coverage of published maps. In both areas the largest volcanic unit covers a significant portion of the surface (58.75% and 63.64%, respectively). For the total measured volcanic units (excluding tessera), these two largest units (that could correspond to the same unit or not) cover the 61.18% and they are stratigraphically superimposed on older volcanic units which cover 3.37% of the area. The remaining area (35.45%) is occupied by younger volcanic units stratigraphically superimposed on the large volcanic unit(s). These results are based on the independent mapping of a large number of geologists with different ideas about the geodynamical evolution of Venus and different criteria for geological mapping. Despite this fact, the presence of these very large units is incompatible with the equilibrium resurfacing models, because their generation at different ages would destroy the crater randomness. Our frequency-area distribution of the mapped volcanic units supports a catastrophic resurfacing due to the emplacement of the largest unit(s) followed by a decay of volcanism. Our data for the frequency-area distribution of volcanic units provide new support for catastrophic resurfacing models. It is difficult to make our observations compatible with equilibrium, steady-state resurfacing models.     

Effects of impacts on the atmospheric evolution: Comparison between Mars, Earth, and Venus

Classified as a terrestrial planet, Venus, Mars, and Earth are similar in several aspects such as bulk composition and density. Their atmospheres on the other hand have significant differences. Venus has the densest atmosphere, composed of CO₂ mainly, with atmospheric pressure at the planet's surface 92 times that of the Earth, while Mars has the thinnest atmosphere, composed also essentially of CO₂, with only several millibars of atmospheric surface pressure. In the past, both Mars and Venus could have possessed Earth-like climate permitting the presence of surface liquid water reservoirs. Impacts by asteroids and comets could have played a significant role in the evolution of the early atmospheres of the Earth, Mars, and Venus, not only by causing atmospheric erosion but also by delivering material and volatiles to the planets. Here we investigate the atmospheric loss and the delivery of volatiles for the three terrestrial planets using a parameterized model that takes into account the impact simulation results and the flux of impactors given in the literature. We show that the dimensions of the planets, the initial atmospheric surface pressures and the volatiles contents of the impactors are of high importance for the impact delivery and erosion, and that they might be responsible for the differences in the atmospheric evolution of Mars, Earth and Venus.     

On the stability of a planet between Mars and the asteroid belt: Implications for the Planet V hypothesis

The stability of an additional planet between the orbit of Mars and the asteroid belt is examined in the context of the Planet V hypothesis. In this model, the Solar System initially contained a fifth terrestrial planet, “Planet V,” which was removed after ∼700 Myr, a possible trigger for the late heavy bombardment on the inner planets. The model is investigated using 96 N-body integrations of the 8 major planets with an additional body between Mars and the asteroid belt. In more than 1/4 of simulations, Planet V survives for 1000 Myr. In most other cases, Planet V collides with the Sun or hits another planet after several hundred Myr, leaving 4 surviving terrestrial planets. In 24/96 simulations, Planet V is lost by ejection or collision with the Sun while the other four terrestrial planets survive without undergoing a collision. In 18 cases, Planet V is removed at least 200 Myr after the beginning of the simulation. The endstate depends sensitively on the mass of Planet V. Collision with the Sun is likely when Planet V's mass is 0.25 Mars masses or less. When Planet V is more massive than this, collisions involving it and/or other terrestrial planets become commonplace. In unstable systems, the times of first encounter and first collision/ejection depend on the initial aphelion distance of Mars. Reducing Mars's aphelion distance increases these times and also increases the fraction of systems surviving for 1000 Myr. When Mars's current orbit is used, the stability of Planet V increases when these two planets are widely separated initially. Planet V's aphelion distance Q typically begins to cross the asteroid belt within a few tens to a few hundred Myr, and its orbit last leaves the belt several hundred Myr later in most cases. The total time spent with Q>2.1 AU is typically less than 200 Myr.