Spectral study of the Eunomia asteroid family : I. Eunomia

We present color ratio curves of the S-Asteroid 15 Eunomia, which have been extracted from high-precision photometric lightcurves obtained in three different VNIR wavelength bands at the Bochum Telescope, La Silla. The measured color ratio curves and near infrared spectra were used to derive a detailed surface composition model whose shape has been computed by V-lightcurve inversions. According to this analysis, the asteroid shows on one hemisphere a higher concentration of pyroxene, which causes an increased 440/700 nm and a reduced 940/700 nm reflectance ratio as well as a pronounced 2-μm absorption band. The remaining surface shows a higher concentration of olivine, leading to a reduced 440/700 nm and slightly increased 940/700 nm color ratio. In addition, we found that the maximum of the 440/700 nm color ratio curve coincide with the minimum of the 940/700 nm color ratio curve and vice versa. We demonstrate on the basis of USGS laboratory spectra that this anti-cyclical behavior can be explained by choosing Fe-rich olivine and a pyroxene with moderate Fe content as varying mineral phases. Furthermore, our observations confirm that 15 Eunomia is an irregular elongated and at least partially differentiated body. Previous spectral investigations of several smaller fragments of the Eunomia asteroid family revealed that the amount of fragments showing an increased pyroxene content exceeds the amount of pyroxene-poor fragments (Nathues, 2000, DLR Forschungsbericht, ISSN 1434-8454). This finding together with the observation that the major fraction of Eunomia's surface is enriched in olivine let us claim that a large fraction of the original pyroxene-enriched crust layer has been lost due to a major collision that created the Eunomia asteroid family. Significant spectral evidences, consistent with high concentrations of metals have been found neither in the rotational resolved spectra of 15 Eunomia nor in its fragments. This led to the conclusion that either no core consisting mainly of metals exists or that an eventual one has not been unearthed by the impact.     

Non-gravitational force modeling of Comet 81P/Wild 2: II. Rotational evolution

In this paper, we have studied both the dynamical and the rotational evolution of an 81P/Wild 2-like comet under the effects of the outgassing-induced force and torque. The main aim is to study if it is possible to reproduce the non-gravitational orbital changes observed in this comet, and to establish the likely evolution of both orbital and rotational parameters. To perform this study, a simple thermophysical model has been used to estimate the torque acting on the nucleus. Once the torque is calculated, Euler equations are solved numerically considering a nucleus mass directly estimated from the changes in the orbital elements (as determined from astrometry). According to these simulations, when the water production rate and changes in orbital parameters for 1997, as well as observational rotational parameters for 2004 are imposed as constraints, the change in the orbital period of 81P/Wild 2, , will decrease so that to , which is similar to the actual tendency observed from 1988 up to 1997. This nearly constant decreasing can be explained as due to a slight drift of the spin axis orientation towards larger ecliptic longitudes. After studying the possible spin axis orientations proposed for 1997, simulations suggest that the spin obliquity and argument (I,Φ)=(56°,167°) is the most likely. As for rotational evolution, changes per orbit smaller than 10% of the actual spin velocity are probable, while the most likely value corresponds to a change between 2 and 7% of the spin velocity. Equally, net changes in the spin axis orientation of 4°-8° per orbit are highly expected.     

Steady-state size distributions for collisional populations:: analytical solution with size-dependent strength

The steady-state population of bodies resulting from a collisional cascade depends on how material strength varies with size. We find a simple expression for the power-law index of the population, given a power law that describes how material strength varies with size. This result is extended to the case relevant for the asteroid belt and Kuiper belt, in which the material strength is described by 2 separate power laws—one for small bodies and one for larger bodies. We find that the power-law index of the small body population is unaffected by the strength law for the large bodies, and vice versa. Simple analytical expressions describe a wave that is superimposed on the large body population because of the transition between the two power laws describing the strength. These analytical results yield excellent agreement with a numerical simulation of collisional evolution. These results will help to interpret observations of the asteroids and KBOs, and constrain the strength properties of those objects.     

Non-gravitational force modeling of Comet 81P/Wild 2: I. A nucleus bulk density estimate

The nucleus of Comet 81P/Wild 2 is modeled by assuming various smooth triaxial ellipsoidal or irregular body shapes, having different rotational periods, spin axis orientations, and thermophysical properties. For these model nuclei, a large number of surface activity patterns (e.g., maps of active and inactive areas) are studied, and in each case the resulting water production rate and non-gravitational force vector versus time are calculated. By requiring that the model nuclei simultaneously reproduce certain properties of the empirical water production curve and non-gravitational changes of the orbit (focusing on changes of the orbital period and in the longitude of perihelion), constraints are placed on several properties of the nucleus. The simulations suggest that the mass of Comet 81P/Wild 2 is , resulting in a rather low bulk density, (depending on the assumed nucleus volume), and that the nucleus rotation is prograde rather than retrograde. The active area fraction is difficult to constrain, but at most 60% of the nucleus is likely to have near-surface ice.     

Loss of water from Mars:: Implications for the oxidation of the soil

The evolution of the martian atmosphere with regard to its H₂O inventory is influenced by thermal loss processes of H, H₂, nonthermal atmospheric loss processes of H⁺, H₂⁺, O, O⁺, CO₂, and O₂⁺ into space, as well as by chemical weathering of the surface soil. The evolution of thermal and nonthermal escape processes depend on the history of the intensity of the solar XUV radiation and the solar wind density. Thus, we use actual data from the observation of solar proxies with different ages from the Sun in Time program for reconstructing the Sun's radiation and particle environment from the present to 3.5 Gyr ago. The correlation between mass loss and X-ray surface flux of solar proxies follows a power law relationship, which indicates a solar wind density up to 1000 times higher at the beginning of the Sun's main sequence lifetime. For the study of various atmospheric escape processes we used a gas dynamic test particle model for the estimation of the pick up ion loss rates and considered pick up ion sputtering, as well as dissociative recombination. The loss of H₂O from Mars over the last 3.5 Gyr was estimated to be equivalent to a global martian H₂O ocean with a depth of about 12 m, which is smaller than the values reported by previous studies. If ion momentum transport, a process studied in detail by Mars Express is significant on Mars, the water loss may be enhanced by a factor of about 2. In our investigation we found that the sum of thermal and nonthermal atmospheric loss rates of H and all nonthermal escape processes of O to space are not compatible with a ratio of 2:1, and is currently close to about 20:1. Escape to space cannot therefore be the only sink for oxygen on Mars. Our results suggest that the missing oxygen (needed for the validation of the 2:1 ratio between H and O) can be explained by the incorporation into the martian surface by chemical weathering processes since the onset of intense oxidation about 2 Gyr ago. Based on the evolution of the atmosphere-surface-interaction on Mars, an overall global surface sink of about 2×10⁴² oxygen particles in the regolith can be expected. Because of the intense oxidation of inorganic matter, this process may have led to the formation of considerable amounts of sulfates and ferric oxides on Mars. To model this effect we consider several factors: (1) the amount of incorporated oxygen, (2) the inorganic composition of the martian soil and (3) meteoritic gardening. We show that the oxygen incorporation has also implications for the oxidant extinction depth, which is an important parameter to determine required sampling depths on Mars aimed at finding putative organic material. We found that the oxidant extinction depth is expected to lie in a range between 2 and 5 m for global mean values.     

Properties of the Hermean Regolith: II. Disk-Resolved Multicolor Photometry and Color Variations of the “Unknown” Hemisphere

Multicolor photometric observations of the “unknown” hemisphere of Mercury have been performed with the Swedish Vacuum Solar Telescope on La Palma at maximal elongations from the Sun in 1997 and 1998. A set of six interference filters with central wavelengths from 450 to 940 nm were used. Multicolor photometry of Mercury was performed on disk-resolved images of the unknown hemisphere (longitudes 160°-340°) with a highest resolution of ∼200 km (J. Warell and S. Limaye 2001, Planet. Space Sci.49, 1531-1552).Disk-integrated spectrophotometry shows that (1) the spectrum of Mercury displays a linear slope from 650 to 940 nm, indicating that the average mercurian regolith is considerably more mature than relatively immature pure anorthosite regions on the Moon; (2) there is negative evidence for the presence of the putative 1-μm absorption feature near 940 nm due to the presence of ferrous iron (Fe²⁺) in pyroxenes; and (3) no effect of phase reddening of the integrated disk is observed between phase angles of 63° and 84°.For the first time, disk-resolved spectrophotometry of Mercury's surface has been obtained, from which it is inferred that (4) the scattering properties of Mercury's regolith are more homogeneous than for the Moon and that there is no clear relation between reflectance and chemical properties at spatial scales of ∼300 km on the unknown hemisphere and (5) there exists an inverse relation of spectral slope with emission angle which is larger for Mercury than for the Moon, indicating that the average mercurian regolith is more backscattering and that this effect increases with wavelength.Finally, from filter ratio images of Mercury's disk it is found that (6) no color variations larger that 2% with respect to the surroundings are detected at a spatial resolution of ∼300 km.     

The contribution of icy grains to the activity of comets: I. Grain lifetime and distribution

We have developed a computer code (GEM—grain evolution model) to simulate the behavior of ice grains in a comet coma. The grains are assumed to be composed of water-ice with an admixture of dark material (“dirt”). An initial size distribution of grains is assumed to be ejected from the nucleus. The ejected mass is taken to be proportional to the rate of gas production by the nucleus. The efficiency for absorption and re-radiation of sunlight is computed from Mie scattering theory. The grain temperature and sublimation rate at a given heliocentric distance is then derived from energy balance considerations. The evolution of the grain size distribution is followed as a function of distance from the nucleus.     

On the dynamics of the jovian ionosphere and thermosphere.: IV. Ion-neutral coupling

We use the fully coupled, three-dimensional, global circulation Jovian Ionospheric Model (JIM) to calculate the coupling between ions in the jovian auroral ovals and the co-existing neutral atmosphere. The model shows that ions subject to drift motion around the auroral oval, as a result of the E×B coupling between a meridional, equatorward electric field and the jovian magnetic field, generate neutral winds in the planetary frame of reference. Unconstrained by the magnetic field, these neutral winds have a greater latitudinal extent than the corresponding ion drifts. Values of the coupling coefficient, k(h), are presented as a function of altitude and cross-auroral electric field strength, for different incoming electron fluxes and energies. The results show that, with ion velocities of several hundred metres per second to over 1 km s⁻¹, k(h) can attain values greater than 0.5 at the ion production peak. This parameter is key to calculating the effective conductivities required to model magnetosphere-ionosphere coupling correctly. The extent to which angular momentum (and therefore energy) is transported vertically in JIM is much more limited than earlier, one-dimensional, studies have predicted.     

“Methane on Mars: A product of H2O photolysis in the presence of CO”: Response to V.A. Krasnopolsky

The detection of CH₄ in the martian atmosphere, at a mixing ratio of about 10 ppb, prompted Krasnopolsky et al. [Krasnopolsky, V.A., Maillard, J.P., Owen, T.C., 2004. Icarus 172, 537-547] and Krasnopolsky [Krasnopolsky, V.A., 2006. Icarus 180, 359-367] to propose that the CH₄ is of biogenic origin. Bar-Nun and Dimitrov [Bar-Nun, A., Dimitrov, V., 2006. Icarus 181, 320-322] proposed that CH₄ can be formed in the martian atmosphere by photolysis of H₂O in the presence of CO. We based our arguments on a clear demonstration that CH₄ is formed in our experiments, and on thermodynamic equilibrium calculations, which show that CH₄ formation is favored even in the presence of oxygen at a mixing ratio 1.3×¹⁰⁻³, as observed on Mars. In the present comment, Krasnopolsky [Krasnopolsky, V.A., 2007. Icarus, in press (this issue)] presents his arguments against the suggestion of Bar-Nun and Dimitrov [Bar-Nun, A., Dimitrov, V., 2006. Icarus 181, 320-322], based on the effect of O₂ on CH₄ formation, the absence of kinetic pathways for CH₄ formation and on the inadequacy of thermodynamic equilibrium calculations to describe the martian atmosphere. In this rebuttal we demonstrate that experiments with molecular oxygen at a ratio of O₂/CO₂=(8.9-17)×¹⁰⁻³, exceeding the martian ratio, still form CH₄. Thermodynamic equilibrium calculations replicate the experimental CH₄ mixing ratio to within a factor of 1.9 and demonstrate that CH₄ production is favored in the martian atmosphere, which is obviously not in thermodynamic equilibrium. Consequently, we do not find the presence of methane to be a sign of biological activity on Mars.     

Direct Monte Carlo and multifluid modeling of the circumnuclear dust coma: Spherical grain dynamics revisited

This paper describes the first computations of dust distributions in the vicinity of an active cometary nucleus, using a multidimensional Direct Simulation Monte Carlo Method (DSMC). The physical model is simplistic: spherical grains of a broad range of sizes are liberated by H₂O sublimation from a selection of nonrotating sunlit spherical nuclei, and submitted to the nucleus gravity, the gas drag, and the solar radiation pressure. The results are compared to those obtained by the previously described Dust Multi-Fluid Method (DMF) and demonstrate an excellent agreement in the regions where the DMF is usable. Most importantly, the DSMC allows the discovery of hitherto unsuspected dust coma properties in those cases which cannot be treated by the DMF. This leads to a thorough reconsideration of the properties of the near-nucleus dust dynamics. In particular, the results show that (1) none of the three forces considered here can be neglected a priori, in particular not the radiation pressure; (2) hitherto unsuspected new families of grain trajectories exist, for instance trajectories leading from the nightside surface to the dayside coma; (3) a wealth of balistic-like trajectories leading from one point of the surface to another point exist; on the dayside, such trajectories lead to the formation of “mini-volcanoes.” The present model and results are discussed carefully. It is shown that (1) the neglected forces (inertia associated with a nucleus rotation, solar tidal force) are, in general, not negligible everywhere, and (2) when allowing for these additional forces, a time-dependent model will, in general, have to be used. The future steps of development of the model are outlined.     

Latitudinal transport by barotropic waves in Titan's stratosphere.: II. Results from a coupled dynamics-microphysics-photochemistry GCM

We present a 2D general circulation model of Titan's atmosphere, coupling axisymmetric dynamics with haze microphysics, a simplified photochemistry and eddy mixing. We develop a parameterization of latitudinal eddy mixing by barotropic waves based on a shallow-water, longitude-latitude model. The parameterization acts locally and in real time both on passive tracers and momentum. The mixing coefficient varies exponentially with a measure of the barotropic instability of the mean zonal flow. The coupled GCM approximately reproduces the Voyager temperature measurements and the latitudinal contrasts in the distributions of HCN and C₂H₂, as well as the main features of the zonal wind retrieved from the 1989 stellar occultation. Wind velocities are consistent with the observed reversal time of the North-South albedo asymmetry of 5 terrestrial years. Model results support the hypothesis of a non-uniform distribution of infrared opacity as the cause of the Voyager temperature asymmetry. Transport by the mean meridional circulation, combined with polar vortex isolation may be at the origin of the latitudinal contrasts of trace species, with eddy mixing remaining restricted to low latitudes most of the Titan year. We interpret the contrasts as a signature of non-axisymmetric motions.     

Optical alteration of complex organics induced by ion irradiation:: 1. Laboratory experiments suggest unusual space weathering trend

Most ion irradiation experiments relevant to primitive outer Solar System objects have been performed on ice and silicate targets. Here we present the first ion irradiation experiments performed on natural complex hydrocarbons (asphaltite and kerite). These materials are very dark in the visible and have red-sloped spectra in the visible and near-infrared. They may be comparable in composition and structure to refractory organic solids on the surfaces of primitive outer Solar System objects. We irradiated the samples with 15-400 keV H⁺, N⁺, Ar⁺⁺, and He⁺ ions and measured their reflectance spectra in the range of 0.3-2.5 μm before ion implantation and after each irradiation step. The results show that irradiation-induced carbonization gradually neutralizes the spectral slopes of these red organic solids. This implies a similar space weathering trend for the surfaces of airless bodies optically dominated by spectrally red organic components. The reduction of spectral slope was observed in all experiments. Irradiation with 30 keV protons, which transfers energy to the target mostly via electronic (inelastic) collisions, showed lower efficiency than the heavier ions. We found that spectral alteration in our experiments increased with increasing contribution of nuclear versus electronic energy loss. This implies that nuclear (elastic) energy deposition plays an important role in changing the optical properties of irradiated refractory complex hydrocarbon materials. Finally, our results indicated that temperature variations from 40 K to room temperature did not influence the spectral properties of these complex hydrocarbon solids.     

Latitudinal transport by barotropic waves in Titan's stratosphere.: I. General properties from a horizontal shallow-water model

We present a numerical study of barotropic waves in Titan's stratosphere based on a shallow-water model. The forcing of the zonal flow by the mean meridional circulation is represented by a relaxation towards a barotropically unstable wind profile. The relaxation profile is consistent with observations and with previous results from a 3D general circulation model. The time constant of the forcing that best matches the northward eddy-transport of zonal momentum from the 3D model is τ∼5 Titan days. The eddy wind field is a zonal wavenumber-2 wave with a peak amplitude about 10% of the mean wind speed. The latitudinal transport of angular momentum by the wave tends to keep the flow close to marginal stability by carrying momentum upgradient, from the core of the jets into the low latitudes. Although the strongest eddy motions occur at the latitudes of the wind maxima, the strongest mixing takes place at the barotropically unstable regions, close to ±30° and spanning about 30° in latitude. An eddy-mixing time constant of the order of 1 Titan day is inferred within these regions, and of a few tens of days within regions of stable flow. Horizontal gradients in transient tracer fields are less than 10% of the latitudinal gradient of the meridional tracer profile. Cassini's detection of such waves could provide a direct observation of wind speeds at stratospheric levels.     

Photometric modeling of Saturn's rings: I. Monte Carlo method and the effect of nonzero volume filling factor

The scattering properties of particulate rings with volume filling factors in the interval D=0.001-0.3 are studied, with photometric Monte Carlo ray tracing simulations combining the advantages of direct (photons followed from the source) and indirect methods (brightness as seen from the observing direction). Besides vertically homogeneous models, ranging from monolayers to classical many-particle thick rings, particle distributions obtained from dynamical simulations are studied, possessing a nonuniform vertical profile and a power law distribution of particle sizes. Self-gravity is not included to assure homogeneity in planar directions. Our main goal is to check whether the moderately flattened ring models predicted by dynamical simulations (with central plane D>0.1) are consistent with the basic photometric properties of Saturn's rings seen in ground-based observations, including the brightening near zero phase angle (opposition effect), and the brightening of the B-ring with increasing elevation angle (tilt effect). Our photometric simulations indicate that dense rings are typically brighter in reflected light than those with D→0, due to enhanced single scattering. For a vertically illuminated layer of identical particles this enhancement amounts at intermediate viewing elevations to roughly 1+2D. Increased single scattering is also obtained for low elevation illumination, further augmented at low phase angles α by the opposition brightening when D increases: the simulated opposition effect agrees very well with the Lumme and Bowell (1981, Astron. J. 86, 1694-1704) theoretical formula. For large α the total intensity may also decrease, due to reduced amount of multiple scattering. For the low (α=13°) and high (α=155°) phase angle geometries analyzed in Dones et al. (1993, Icarus 105, 184-215) the brightness change for D=0.1 amounts to 20% and −17%, respectively. In the case of an extended size distribution, dynamical simulations indicate that the smallest particles typically occupy a layer several times thicker than the largest particles. Even if the large particles form a dynamically dense system, a narrow opposition peak can arise due to mutual shadowing among the small particles: for example, a size distribution extending about two decades can account for the observed about 1° wide opposition peak, solely in terms of mutual shadowing. The reduced width of the opposition peak for extended size distribution is in accordance with Hapke's (1986, Icarus 67, 264-280) treatment for semi-infinite layers. Due to vertical profile and particle size distribution, the photometric behavior is sensitive to the viewing elevation: this can account for the tilt-effect of the B-ring, as dense and thus bright central parts of the ring become better visible for larger elevation, whereas in the case of smaller elevation, mainly low volume density upper layers are visible. Since multiple scattering is not involved, the explanation works also for albedo well below unity. Inclusion of nonzero volume density helps also to model some of the Voyager observations. For example, the discrepancy between predicted and observed brightness at large phase angles for much of the A-ring (Dones et al., 1993, Icarus 105, 184-215) is removed when the enhanced low α single scattering and reduced large α multiple scattering is allowed for. Also, a model with vertical thickness increasing with saturnocentric distance offers at least a qualitative explanation for the observed contrast reversal between the inner and outer A-ring in low and high phase Voyager images. Differences in local size distribution and thus on the effective D may also account for the contrast reversal in resonance sites.