Structure and density of cometary nuclei

Abstract— Understanding the nature of the cometary nucleus remains one of the major problems in solar system science. Whipple's (1950) icy conglomerate model has been very successful at explaining a range of cometary phenomena, including the source of cometary activity and the nongravitational orbital motion of the nuclei. However, the internal structure of the nuclei is still largely unknown. We review herein the evidence for cometary nuclei as fluffy aggregates or primordial rubble piles, as first proposed by Donn et al. (1985) and Weissman (1986). These models assume that cometary nuclei are weakly bonded aggregations of smaller, icy‐conglomerate planetesimals, possibly held together only by self‐gravity. Evidence for this model comes from studies of the accretion and subsequent evolution of material in the solar nebula, from observations of disrupted comets, and in particular comet Shoemaker‐Levy 9, from measurements of the ensemble rotational properties of observed cometary nuclei, and from recent spacecraft missions to comets. Although the evidence for rubble pile nuclei is growing, the eventual answer to this question will likely not come until we can place a spacecraft in orbit around a cometary nucleus and study it in detail over many months to years. ESA's Rosetta mission, now en route to comet 67P/Churyumov‐Gerasimenko, will provide that opportunity.     

Apollo 12 Lunar Module exhaust plume impingement on Lunar Surveyor III

Understanding plume impingement by retrorockets on the surface of the Moon is paramount for safe lunar outpost design in NASA’s planned return to the Moon for the Constellation Program. Visual inspection, Scanning Electron Microscopy, and surface scanned topology have been used to investigate the damage to the Lunar Surveyor III spacecraft that was caused by the Apollo 12 Lunar Module’s close proximity landing. Two parts of the Surveyor III craft returned by the Apollo 12 astronauts, Coupons 2050 and 2051, which faced the Apollo 12 landing site, show that a fine layer of lunar regolith coated the materials and was subsequently removed by the Apollo 12 Lunar Module landing rocket. The coupons were also pitted by the impact of larger soil particles with an average of 103 pits/cm². The average entry size of the pits was 83.7 μm (major diameter) × 74.5 μm (minor diameter) and the average estimated penetration depth was 88.4 μm. Pitting in the surface of the coupons correlates to removal of lunar fines and is likely a signature of lunar material imparting localized momentum/energy sufficient to cause cracking of the paint. Comparison with the lunar soil particle size distribution and the optical density of blowing soil during lunar landings indicates that the Surveyor III spacecraft was not exposed to the direct spray of the landing Lunar Module, but instead experienced only the fringes of the spray of soil. Had Surveyor III been exposed to the direct spray, the damage would have been orders of magnitude higher.     

Processing of cometary grains at the nucleus surface

Cometary material inevitably undergoes chemical changes before and on leaving the nucleus. In seeking to explain comets as the origin of many IDPs (interplanetary dust particles), an understanding of potential surface chemistry is vital. Grains are formed and transformed at the nucleus surface; much of the cometary volatiles may arise from the organic material. In cometary near-surface permafrost, one expects cryogenic chemistry with crystal growth and isotope. This could be the hydrous environment where IDPs form. Seasonal and geographic variations imply a range of environmental conditions and surface evolution. Interplanetary dust impacts and electrostatic forces also have roles in generating cometary dust. The absence of predicted cometary dust ‘envelopes’ is compatible with the wide range of particle structures and compositions. Study of IDPs would distinguish between this model and alternatives that see comets as aggregates of core-mantle grains built in interstellar clouds.     

Processing of cometary grains at the nucleus surface

Cometary material inevitably undergoes chemical changes before and on leaving the nucleus. In seeking to explain comets as the origin of many IDPs (interplanetary dust particles), an understanding of potential surface chemistry is vital. Grains are formed and transformed at the nucleus surface; much of the cometary volatiles may arise from the organic material. In cometary near-surface permafrost, one expects cryogenic chemistry with crystal growth and isotope. This could be the hydrous environment where IDPs form. Seasonal and geographic variations imply a range of environmental conditions and surface evolution. Interplanetary dust impacts and electrostatic forces also have roles in generating cometary dust. The absence of predicted cometary dust envelopes is compatible with the wide range of particle structures and compositions. Study of IDPs would distinguish between this model and alternatives that see comets as aggregates of core-mantle grains built in interstellar clouds.     

A model of the macro-features of the surface layer of a cometary nucleus

The macro-features of the surface layer of a 'fresh' cometary nucleus are modelled by assuming that the dust and the snow particles of which it consists both have a mass distribution index of 1.65, and that the dust/gas mass ratio is 0.45. Conclusions are drawn as to how this model helps us to understand the cometary sublimation process and the cometary surface layer. The latter most probably consists of weak, low-density, friable, slightly dusty snow. Its ability effectively to support even the small weight of, say, the Rosetta landing probe is in considerable doubt.     

Simulation of ICE-Giacobini-Zinner and VEGA/GIOTTO-Halley encounters

The formation of a weak shock at the interaction of the solar wind with the cometary neutral gas is confirmed in a non-stationary gasdynamic model. No shocks will be detectable along the ICE path and if the VEGA/GIOTTO spacecraft enter the cometary atmosphere. A shock can only be expected at an exit from the Halley atmosphere.     

Hydrocode modeling of oblique impacts: The fate of the projectile

Abstract— All impacts are oblique to some degree. Only rarely do projectiles strike a planetary surface (near) vertically. The effects of an oblique impact event on the target are well known, producing craters that appear circular even for low impact angles (>15° with respect to the surface). However, we still have much to learn about the fate of the projectile, especially in oblique impact events. This work investigates the effect of angle of impact on the projectile. Sandia National Laboratories' three‐dimensional hydrocode CTH was used for a series of high‐resolution simulations (50 cells per projectile radius) with varying angle of impact. Simulations were carried out for impacts at 90, 60, 45, 30, and 15° from the horizontal, while keeping projectile size (5 km in radius), type (dunite), and impact velocity (20 km/s) constant. The three‐dimensional hydrocode simulations presented here show that in oblique impacts the distribution of shock pressure inside the projectile (and in the target as well) is highly complex, possessing only bilateral symmetry, even for a spherical projectile. Available experimental data suggest that only the vertical component of the impact velocity plays a role in an impact. If this were correct, simple theoretical considerations indicate that shock pressure, temperature, and energy would depend on sin²θ, where θ is the angle of impact (measured from the horizontal). However, our numerical simulations show that the mean shock pressure in the projectile is better fit by a sin θ dependence, whereas shock temperature and energy depend on sin^3/2 θ. This demonstrates that in impact events the shock wave is the result of complex processes that cannot be described by simple empirical rules. The mass of shock melt or vapor in the projectile decreases drastically for low impact angles as a result of the weakening of the shock for decreasing impact angles. In particular, for asteroidal impacts the amount of projectile vaporized is always limited to a small fraction of the projectile mass. In cometary impacts, however, most of the projectile is vaporized even at low impact angles. In the oblique impact simulations a large fraction of the projectile material retains a net downrange motion. In agreement with experimental work, the simulations show that for low impact angles (30 and 15°), a downrange focusing of projectile material occurs, and a significant amount of it travels at velocities larger than the escape velocity of Earth.     

SMART‐1 end of life shallow regolith impact simulations

The SMART‐1 end‐of‐life impact with the lunar surface was simulated with impacts in a two stage light‐gas gun onto inclined basalt targets with a shallow surface layer of sand. This simulated the probable impact site, where a loose regolith will have overlaid a well consolidated basaltic layer of rock. The impact angles used were at 5° and 10° from the horizontal. The impact speed was ~2 km s^?1 and the projectiles were 2.03 mm diameter aluminum spheres. The sand depth was between approximately 0.8 and 1.8 times the projectile diameter, implying a loose lunar surface regolith of similar dimensions to the SMART‐1 spacecraft. A crater in the basement rock itself was only observed in the impact at 10° incidence, and where the depth of loose surface material was less than the projectile diameter, in which case the basement rock also contained a small pit‐like crater. In all cases, the projectile ricocheted away from the impact site at a shallow angle. This implies that at the SMART‐1 impact site the crater will have a complicated structure, with exposed basement rock and some excavated rock displaced nearby, and the main spacecraft body itself will not be present at the main crater.     

MASS FLUX OF ASTEROIDAL ORIGIN METEOROIDS ON PERIODIC COMET NUCLEI

We examine the potential contamination of cometary nuclei through impacts from asteroidal origin meteoroids. The paper uses a simple model and has the goal of determining whether asteroidal contamination is potentially significant. We assume a meteoroid power law mass distribution with index values in the range from s=1.83 to s=2.09. We used maximum and minimum models which we believe will bracket the true meteoroid mass distribution. We identify those comets which are expected to be most significantly contaminated, and find values of up to 3.6 kg of asteroidal meteoroid impact per square meter of the cometary surface per orbital revolution. This is less than the expected mass loss per perihelion passage for most comets. Therefore any remnant effects of the contamination will depend on the penetration depth of the meteoroids in the cometary nucleus, and possibly on the distribution of active and inactive areas on cometary nuclei. We present a simple model which suggests that even small meteoroids will embed relatively deeply into a cometary nucleus.     

CoMA — An advanced space experiment forin situ analysis of cometary matter

Exploring one of the most pristine bodies in our solar system — a comet — with a spacecraft will be a great step towards a deeper understanding of our solar system's beginnings. We here present the advanced space experiment CoMA (cometarymatteranalyzer), which will be flown on NASA's cometary rendezvous and asteroid flyby mission CRAF. CoMA is a high resolution time-of-flight secondary ion mass spectrometer. It will analyze m-sized cometary dust grains and cometary gases with an unprecedented mass resolution and will yield data about the elemental, isotopic, and molecular composition.     

Deep penetration of projectile material into the planet during the impact

A mechanism through which water could be buried inside the Moon is found. If an icy comet strikes the planetary surface and a thin natural crack exists at the site of the impact, some amount of cometary material can penetrate deep into the ground. This happens due to peculiar features of hydrodynamic flow along the crack. Numerical simulations based on the free-Lagrangian method show that the amount of water buried under the crater is several percent of the original mass of the projectile.     

Correction of the ionospheric distortion on the MARSIS surface sounding echoes

Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) on the Mars Express (MEX) spacecraft has made numerous measurements of the Martian surface and subsurface. However, all of these measurements are distorted by the ionosphere and must be compensated before any analysis. We have developed a technique to compensate for the ionospheric distortions. This technique provides a powerful tool to derive the total electron content (TEC) and other higher-order terms of the limited expansion of the plasma dispersion function that are related to overall shape of the electron column profile. The derived parameters are fitted by using a Chapman model to derive ionospheric parameters like n₀, electron density primary peak (maximum for solar zenith angle (SZA) equal 0), and the neutral height scale H.

Our estimated ionospheric parameters are in good agreement with Mars Global Surveyor (MGS) radio-occultation data. However, since MARSIS does not have the observation geometry limitations of the radio occultation measurements, our derived parameters extend over a large range of SZA for each MEX orbit.

The first results from our technique have been discussed by Safaeinili et al. [2007, Estimation of the total electron content of the Martian ionosphere using radar sounder surface echoes. Geophys. Res. Lett. 34, L23204, doi:10.1029/2007GL032154].     

Properties of filamentary sublimation residues from dispersions of clay in ice

During vacuum sublimation experiments on simulated Martian polar deposits and cometary dirty ices, a fluffy filamentary sublimate residue material with unique physical properties was produced. The silica-to-silica bonds that we believe join the particles together are the result of conditions that may exist in some Martian polar deposits and on some cometary surfaces. Submicron particles of montmorillonite clay thinly dispersed (1: 1000 clay/water) and not contacting one another in water ice can form very-low-density structures (density as low as 10⁻³ g cm⁻³) during sublimation of the ice. The lightweight constructs, when viewed in scanning electron microscopy micrographs, are composed of long network chains of the clay particles. The material is sufficiently electrically conductive to drain away the scanning electron microscopy charge. It is also resistant (no change in electronic properties are apparent) to scanning electron microscopy electron-beam heating for hours in vacuo. Infrared spectra and X-ray diffraction patterns of the sublimate residues show little difference from spectra and patterns of the original minerals. Heating in an oven, in air, to 370°C produces little change in the structure of the sublimate residual material. The particle bonding forces are strong and produce a resilient, elastic lightweight material. The particle bonding will allow vapors to diffuse through it, and its thermal conductivity is very low. These properties produce a high-performance vacuum insulation. This material may have applications for insulating ice bodies (solid cryogens) in space. The incoming heat is partially carried away by the out-flowing water vapor.     

Rosetta mission mutual impedance probe modelling: the short and long Debye length plasma cases

The Rosetta spacecraft (S/C), which is planned to meet comet 46P/Wirtanen in 2011, will carry a set of five wave and plasma instruments (i.e. the Rosetta Plasma Consortium). This is to measure the cometary plasma properties from the minimum value of activity of the comet to its maximum value at perihelion. The mutual impedance probe, MIP, is one of those (Trotignon et al., 1999) five. This revised version was published online in July 2006 with corrections to the Cover Date.     

Photometric properties of Titan''s surface from Cassini VIMS: Relevance to titan''s hemispherical albedo dichotomy and surface stability

The Visual and Infrared Mapping Spectrometer (VIMS) instrument on the Cassini Saturn Orbiter returned spectral imaging data as the spacecraft undertook six close encounters with Titan beginning 7 July, 2004. Three of these flybys each produced overlapping coverage of two distinct regions of Titan's surface. Twenty-four points were selected on approximately opposite hemispheres to serve as photometric controls. Six points were selected in each of four reflectance classes. On one hemisphere each control point was observed at three distinct phase angles. From the derived phase coefficients, preliminary normal reflectances were derived for each reflectance class. The normal reflectance of Titan's surface units at 2.0178 μm ranged from 0.079 to 0.185 for the most absorbing to the most reflective units assuming no contribution from absorbing haze. When a modest haze contribution of τ=0.1 is considered these numbers increase to 0.089–0.215. We find that the lowest three reflectance classes have comparable normal reflectance on either hemisphere. However, for the highest brightness class the normal reflectance is higher on the hemisphere encompassing longitude 14–65° compared to the same high brightness class for the hemisphere encompassing 122–156° longitude. We conclude that an albedo dichotomy observed in continental sized units on Titan is due not only to one unit having more areal coverage of reflective material than the other but the material on the brighter unit is intrinsically more reflective than the most reflective material on the other unit. This suggests that surface renewal processes are more widespread on Titan's more reflective units than on its less reflective units.

We note that one of our photometric control points has increased in reflectance by 12% relative to the surrounding terrain from July of 2004 to April and May of 2005. Possible causes of this effect include atmospheric processes such as ground fog or orographic clouds; the suggestion of active volcanism cannot be ruled out.

Several interesting circular features which resembled impact craters were identified on Titan's surface at the time of the initial Titan flyby in July of 2004. We traced photometric profiles through two of these candidate craters and attempted to fit these profiles to the photometric properties expected from model depressions. We find that the best-fit attempt to model these features as craters requires that they be unrealistically deep, approximately 70 km deep. We conclude that despite their appearance, these circular features are not craters, however, the possibility that they are palimpsests cannot be ruled out.

We used two methods to test for the presence of vast expanses of liquids on Titan's surface that had been suggested to resemble oceans. Specular reflection of sunlight would be indicative of widespread liquids on the surface; we found no evidence of this. A large liquid body should also show uniformity in photometric profile; we found the profiles to be highly variable. The lack of specular reflection and the high photometric variability in the profiles across candidate oceans is inconsistent with the presence of vast expanses of flat-lying liquids on Titan's surface. While liquid accumulation may be present as small, sub-pixel-sized bodies, or in areas of the surface which still remain to be observed by VIMS, the presence of large ocean-sized accumulations of liquids can be ruled out.

The Cassini orbital tour offers the opportunity for VIMS to image the same parts of Titan's surface repeatedly at many different illumination and observation geometries. This creates the possibility of understanding the properties of Titan's atmosphere and haze by iteratively adapting models to create a best fit to the surface reflectance properties.     

Impact cratering of granular mixture targets made of H2O Ice–CO2 Ice–pyrophylite

Experiments related to impacts onto three-component targets which could simulate cometary nucleus or planetary regolith cemented by ices are presented here. The impact velocities are from 133 to 632 m s⁻¹. The components are powdered mineral (pyrophylite), H₂O ice, and CO₂ ice mixed 1:1:0.74 by mass. The porosity of fresh samples is about 0.48. Two types of the samples were studied: nonheated samples and samples heated by thermal radiation. Within the samples a layered structure was formed. The cratering pattern strongly depended on the history of the samples. The craters formed in nonheated targets had regular shapes. The volume was easy to be determined and it was proportional to impact energy E. The crater depth scales as E^0.5. Impacts on the thermally stratified target led to ejection of a large amount of material from the loose sub-crustal layer. For some particular interval of impact velocity a cratering pattern can demonstrate unusual properties: small hole through the rigid crust and considerable mass transfer (radially, outward of the impact point) within sub-crustal layer.     

On the Light-Absorbing Surface Layer of Cometary Nuclei: II. Thermal Modeling

The classical way to treat absorption of solar light in thermophysical modeling of cometary nuclei (and other ice-rich bodies such as jovian satellites) has been to assume complete opaqueness of the surface material. However, as shown by Davidsson and Skorov (2002, Icarus156, 223-248), substantial light penetration can occur in porous ice even if it is very dusty, implying that gradual absorption of energy in a surface layer should be accounted for.We present a thorough comparison between a surface energy absorption model and a layer energy absorption model, for various combinations of heliocentric distances, conductivities, opacities, pore sizes, and rotational periods relevant for cometary nuclei, by fully solving the coupled differential equations of heat transfer and gas diffusion. We find substantial differences between the models in terms of gas production rate, thermal lag angle, surface temperature, and the origin of coma molecules. For example, the surface energy absorption model overestimates the total gas production by a factor of 2-7, underestimates the lag angle by a factor of 2-3, and places the origin of coma molecules at the surface, instead of the near-surface interior.     

Influence of the magnetic field on the density distribution of solar wind protons and cometary ions in the shock layer ahead of cometary ionospheres

The “mass loading” of the solar wind by cometary ions produced by the photoionization of neutral molecules outflowing from the cometary nucleus plays a major role in the interaction of the solar wind with cometary atmospheres. In particular, this process leads to a decrease in the solar wind velocity with a transition from supersonic velocities to subsonic ones through the bow shock. The so-called single-fluid approximation, in which the interacting plasma flows are considered as a single fluid, is commonly used in modeling such an interaction. However, it is occasionally necessary to know the distribution of parameters for the components of the interacting plasma flows. For example, when the flow of the cometary dust component in the interplanetary magnetic field is considered, the dust particle charge, which depends significantly on the composition of the surrounding plasma, needs to be known. In this paper, within the framework of a three-dimensional magnetohydrodynamic model of the solar wind flow around cometary ionospheres, we have managed to separately obtain the density distributions of solar wind protons and cometary ions between the bow shock and the cometary ionopause (in the shock layer). The influence of the interplanetary magnetic field on the position of the point of intersection between the densities with the formation of a region near the ionopause where the proton density is essentially negligible compared to the density of cometary ions is investigated. Such a region was experimentally detected by the Vega-2 spacecraft when investigating Comet Halley in March 1986. The results of the model considered below are compared with some experimental data obtained by the Giotto spacecraft under the conditions of flow around Comets Halley and Grigg–Skjellerup in 1986 and 1992, respectively. Unfortunately, our results of calculations on Comet Churyumov–Gerasimenko are only predictive in character, because the trajectory of the Rosetta spacecraft, which manoeuvred near its surface for several months, is complex.     

Physico-Chemistry of Comets: Models and Laboratory Experiments

Until cometary matter can be studied in-situ or cometary samples are brought back to Earth for analysis, theoretical models and laboratory studies remain a crucial tool for revealing the nature of cometary matter. Constraints on the nature of the primordial material available for incorporation into comets and other solar system material comes from analysis of data from space-based and ground-based observatories. The structure of the nuclear ice component, which may have coexisting amorphous/crystalline phases and include clathrates and other trapped guest molecules, strongly influences the cometary outgassing properties. This paper reviews laboratory work on ice and carbon aceous compounds and discusses their significance for cometary chemistry. Special emphasis will be given to studies on the thermal processing of ices and their implications for the structure changes and subsequent release of volatiles. We also describe the preliminary results of a model of nuclear outgassing, and discuss how such model scan be used to infer the chemical structure of the nuclearices. Furthermore, we confront cometary data with the analysis of carbonaceous meteorites. Recent laboratory results on volatile compounds and the macro molecular structure of carbonaceous meteorites allow us to investigate the link of small bodies in the Solar System. Until ROSETTA will land on comet Wirtanen and study directly the nuclear composition, laboratory measurements of ice and refractory analogs will — together with the analysis of meteorites —significantly improve our knowledge on the origin and structure ofcomets.