443 Eros: Problems with the meteorite magnetism record in attempting an asteroid match

Abstract— The magnetometer experiment (MAG) onboard the Near‐Earth Asteroid Rendezvous (NEAR)‐Shoemaker spacecraft detected no global scale magnetization and established a maximum magnetization of 2.1 times 10^?6 Am² kg^?1 for asteroid 433 Eros. This is in sharp contrast with the estimated magnetization of other S‐class asteroids (Gaspra, ?2.4 times 10^?2 Am² kg^?1; Braille, ?2.8 times 10^?2 Am² kg^?1) and is below published values for all types of ordinary chondrites. This includes the L/LL types considered to most closely match 433 Eros based on preliminary interpretations of NEAR remote geochemical experiments. The ordinary chondrite meteorite magnetization intensity data was reviewed in order to assess the reasonableness of an asteroid‐meteorite match based on magnetic property measurements. Natural remanent magnetization (NRM) intensities for the ordinary chondrite meteorites show at least a 2 order of magnitude range within each of the H, L, and LL groups, all well above the 2.1 times 10^?6 Am² kg^?1 level for 433 Eros. The REM values (ratio of the NRM to the SIRM (saturation remanent magnetization)) range over 3 orders of magnitude for all chondrite groups indicating no clear relationship between NRM and the amount of magnetic material. Levels of magnetic noise in chondrite meteorites can be as much as 70% or more of the NRM. Consequently, published values of the NRM should be considered suspect unless careful evaluation of the noise sources is done. NASA Goddard SFC studies of per unit mass intensities in large (>10 000 g) and small (down to <1 g) samples from the same meteorite demonstrate magnetic intensity decreases as size increases. This would appear to be explained by demagnetization due to magnetic vector randomness at unknown scale sizes in the larger samples. This would then argue for some level of demagnetization of large objects such as an asteroid. The possibility that 433 Eros is an LL chondrite cannot be discounted.     

Five-color photoelectric photometry of asteroid 433 Eros

Five-color photoelectric lightcurves of asteroid 433 Eros were obtained on 9 nights during the 1974/75 apparition. Although color differences due to changing solar phase angle were detected, color differences during a single rotation of Eros are less than 1%. Amplitudes of up to 1^m.44 were measured, and there are indications that three reversals in the relative depths of the two minima occured between late December 1974 and late January 1975. The absolute visual magnitude at primary maximum, corrected to zero phase and to one AU from Earth and Sun, is about V₀(1,0) = 10^m.8.     

Photometric lightcurves and pole determination of 433 Eros

Fourteen photometric lightcurves of 433 Eros were made at the Astronomical Observatory of Torino during the 1974–75 close passage. The absolute magnitude of the primary maximum (10^m78), the phase coefficient (0.023 mag/degree), the synodic and sidereal period of rotation (0^d.21956 and 0^d.21959, respectively) and the ecliptic coordinates of the pole (λ = 17°, β = 10°) were deduced.     

Magnetic field investigations during ROSETTA's 2867 Šteins flyby

During the 2867 Šteins flyby of the ROSETTA spacecraft on September 5, 2008 magnetic field measurements have been made with both the RPC orbiter magnetometer and the ROMAP lander magnetometer. These combined magnetic field measurements allow a detailed examination of any magnetic signatures caused either directly by the asteroid or indirectly by Šteins’ different modes of interaction with the solar wind. Comparing measurements with simulation results show that Šteins does not posses a significant remanent magnetization. The magnetization is estimated at less than 10⁻³ A m²/kg. This is significantly different from results at 9969 Braille and 951 Gaspra.     

Simulations of impact ejecta and regolith accumulation on Asteroid Eros

We have carried out a set of Monte Carlo simulations of the placement of impact ejecta on Asteroid 433 Eros, with the aim of understanding the distribution and accumulation of regolith. The simulations consisted of two stages: (1) random distribution of primary impact sites derived from a uniform isotropic flux of impactors, and (2) integration of the orbits of test particle ejecta launched from primary impact points until their re-impact or escape. We integrated the orbits of a large number of test particles (typically 10⁶ per individual case). For those particles that did not escape we collected the location of their re-impact points to build up a distribution on the asteroid surface. We find that secondary impact density is mostly controlled by the overall topography of the asteroid. A gray-scale image of the density of secondary ejecta impact points looks, in general, like a reduced-scale negative of the topography of the asteroid's surface. In other words, regolith migration tends to fill in the topography of Eros over time, whereas topographic highs are denuded of free material. Thus, the irregular shape of Eros is not a steady-state configuration, but the result of larger stochastic events.     

Lightcurves and the axis of rotation of 433 Eros

Ten lightcurves and UBV photometry of 433 Eros were obtained between August 1972 and May 1975. The absolute magnitude of the lightcurve maximum is 10.75 and the phase coefficient is 0.025 mag/deg. There may be a small difference in B-V color between the northern and southern hemispheres. The pole of the axis of rotation is directed toward λ₀ = 16°, β₀ = 12°, ecliptic longitude and latitude, respectively, and the rotation is direct with a sidereal period of 0.^d219599 or 5^h16^m13^s4 ± 0.^s2. The dimensions derived from the polarimetric albedo and the lightcurve amplitudes are 12km × 12km × 31km for a smooth cylinder with hemispherical ends.     

An Estimate of Eros's Porosity and Implications for Internal Structure

Earth-based spectral measurements and NEAR Shoemaker magnetometer, X-ray, and near-infrared spectrometer data are all consistent with Eros having a bulk composition and mineralogy similar to ordinary chondrite meteorites (OC). By comparing the bulk density of 433 Eros (2.67±0.03 g/cm³) with that of OCs (3.40 g/cm³), we estimate the total porosity of the asteroid to be 21-33%. Macro (or structural) porosity, best estimated to be ∼20%, is constrained to be between 6 and 33%. We conclude that Eros is a heavily fractured body, but we find no evidence that it was ever catastrophically disrupted and reaccumulated into a rubble pile.     

Studies of microparticle impact phenomena leading to the development of a highly sensitive micrometeoroid detector

A highly sensitive device has been developed for the detection of micrometeoroids with masses extending down to and below the expected classical solar radiation pressure limit. This limit occurs at a mass of about 10⁻¹⁶ kg for a particle density of 3 × 10³ kg m⁻³. The detector operates on the principle that charge (Q) is released when a hypervelocity microparticle impacts on a solid surface. This process has previously been investigated at particle velocities (ν) greater than 1 km sec⁻¹ and the empirical relationship Q = Km^αν^β obtained where m is the particle mass and K, α and β are constants. In the present study the validity of this relationship has been demonstrated for iron particles with masses from 10⁻¹⁶ to 10⁻¹³ kg in the velocity range 0.05–1.4 km sec⁻¹ impacting on a molybdenum target. A value of α ≈ 1 is indicated and a value of β of 3.2 ± 0.4 has been obtained.     

A Global Solution for the Gravity Field, Rotation, Landmarks, and Ephemeris of Eros

As part of the NEAR Radio Science investigation, a global solution that includes both spherical and ellipsoidal harmonic gravity fields of Eros, Eros pole and rotation rate, Eros ephemeris, and landmark positions from the optical data was generated. This solution uses the entire one-year in orbit collection of X-band radiometric tracking (Doppler and range) from the Deep Space Network and landmark tracking observations generated from the NEAR spacecraft images of Eros. When compared to a constant density shape model, the gravity field shows a nearly homogeneous Eros. The Eros landmark solutions are in good agreement with the Eros shape model, and they reduce the center-of-mass and center-of-figure offset in the z direction to 13 m. Most of the NEAR spacecraft orbits are determined in all directions to an accuracy of several meters. The solution for the ephemeris of Eros constrains the mass of Vesta to 18.2±0.4 km³/s² and reduces the uncertainty in the Earth-Moon mass ratio.     

An experimental study of the formation of an ice crust and migration of water vapor in a comet's upper layers

From recent close encounters with Comets Wild-2 and Tempel 1 we learned that their surfaces are very rugged and no simple uniform layers model can be applied to them. Rather, a glaciological approach should be applied for describing their surface features and behavior. Such intrinsically rugged surface is formed in our large scale experiments, where an agglomerate of ∼200 μm gas-laden amorphous ice particles is accumulated to form a 20 cm diameter and few cm high ice sample. The density, tensile strength and thermal inertia of our ice sample were found to be very close to those found by Deep Impact for Comet Tempel 1: density 250-300 kg m⁻³ vs DI 350-400 kg m⁻³; tensile strength 2-4 kPa vs DI 1-10 kPa; thermal inertia 80 W K⁻¹ m⁻² s^1/2 vs <100 W K⁻¹ m⁻² s^1/2 and <50 W K⁻¹ m⁻² s^1/2. From the close agreement between the thermal inertias measured in our ice sample, which had no dust coverage and that of Comet Tempel 1, we deduce that the low thermal inertia is an intrinsic property of the fluffy structure of the ice as a result of its low density, with an addition by the broken terrain and not due to the formation of a dust layer. Upon warming up of the ice, water vapor migrates both outward into the coma and inward. Reaching cooler layers, the water vapor condenses, forming a denser ice crust, as we show experimentally. We also demonstrate the inward and outward flow of water vapor in the outer ice layers through the exchange between layers of D₂O ice and H₂O ice, to form HDO.     

Elemental composition of 433 Eros: New calibration of the NEAR-Shoemaker XRS data

We present a new calibration of the elemental-abundance data for Asteroid 433 Eros taken by the X-ray spectrometer (XRS) aboard the NEAR-Shoemaker spacecraft. (NEAR is an acronym for “Near-Earth Asteroid Rendezvous.”) Quantification of the asteroid surface elemental abundance ratios depends critically on accurate knowledge of the incident solar X-ray spectrum, which was monitored simultaneously with asteroid observations. Previously published results suffered from incompletely characterized systematic uncertainties due to an imperfect ground calibration of the NEAR gas solar monitor. The solar monitor response function and associated uncertainties have now been characterized by cross-calibration of a large sample of NEAR solar monitor flight data against contemporary broadband solar X-ray data from the Earth-orbiting GOES-8 (Geostationary Operational Environmental Satellite). The results have been used to analyze XRS spectra acquired from Eros during eight major solar flares (including three that have not previously been reported). The end product of this analysis is a revised set of Eros surface elemental abundance ratios with new error estimates that more accurately reflect the remaining uncertainties in the solar flare spectra: Mg/Si=0.753+0.078/−0.055, Al/Si=0.069±0.055, S/Si=0.005±0.008, Ca/Si=0.060+0.023/−0.024, and Fe/Si=1.678+0.338/−0.320. These revised abundance ratios are consistent within cited uncertainties with the results of Nittler et al. [Nittler, L.R., and 14 colleagues, 2001. Meteorit. Planet. Sci. 36, 1673-1695] and thus support the prior conclusions that 433 Eros has a major-element composition similar to ordinary chondrites with the exception of a strong depletion in sulfur, most likely caused by space weathering.     

A precise modeling of Eros 433 rotation

Thanks to the recent data obtained from the NEAR space probe, we calculate in this paper, with a precision never reached so far for an asteroid, the precession and the nutation of Eros 433. In a preliminary step, we show that Eros obliquity has a remarkable value of 89.0° which tends to align its figure axis along the orbital plane. This very specific obliquity has some consequences on the motion of the axis of figure: one is the very small amplitude of the precession in longitude, for which we get the value . Moreover, we calculate Eros nutation for the figure axis due to the Sun, after developing the perturbing potential at the 4th order of the eccentricity. We show that the figure axis undergoes very large oscillations in the direction perpendicular to Eros orbital plane, due to the nutation in obliquity. Peak to peak, these oscillations reach 55″, which is far larger than the amplitudes of the nutations of the Earth due to the Sun (of the order of 2″). Moreover, we give the analytical developments of Δψ and Δε, both for the axis of angular momentum, and the axis of figure.     

Inflight Calibration of the NEAR Multispectral Imager: II. Results from Eros Approach and Orbit

During the Near-Earth Asteroid Rendezvous (NEAR) spacecraft's investigation of asteroid 433 Eros, inflight calibration measurements from the multispectral imager (MSI) have provided refined knowledge of the camera's radiometric performance, pointing, and light-scattering characteristics. Measurements while at Eros corroborate most earlier calibration results, although there appears to be a small, gradual change in instrument dark current and flat field due to effects of aging in the space environment. The most pronounced change in instrument behavior, however, is a dramatic increase in scattered light due to contaminants accumulated on the optics during unscheduled fuel usage in December 1998. Procedures to accurately quantify and to remediate the scattered light are described in a companion paper (Li et al. 2002, Icarus155, 00-00). Acquisition of Eros measurements has clarified the relative, filter-to-filter, radiometric performance of the MSI. Absolute radiometric calibration appears very well constrained from flight measurements, with an accuracy of ∼5%. Pointing relative to the spacecraft coordinate system can be determined from the temperature of the spacecraft deck with an accuracy of ∼1 pixel.     

In situ measurements of particle load and transport in dust devils

In situ (mobile) sampling of 33 natural dust devil vortices reveals very high total suspended particle (TSP) mean values of 296 mg m⁻³ and fine dust loadings (PM10) mean values ranging from 15.1 to 43.8 mg m⁻³ (milligrams per cubic meter). Concurrent three-dimensional wind profiles show mean tangential rotation of 12.3 m s⁻¹ and vertical uplift of 2.7 m s⁻¹ driving mean vertical TSP flux of 1689 mg m⁻³ s⁻¹ and fine particle flux of ∼1.0 to ∼50 mg m⁻³ s⁻¹. Peak PM10 dust loading and flux within the dust column are three times greater than mean values, suggesting previous estimates of dust devil flux might be too high. We find that deflation rates caused by dust devil erosion are ∼2.5-50 μm per year in dust devil active zones on Earth. Similar values are expected for Mars, and may be more significant there where competing erosional mechanisms are less likely.     

Properties of the dust cloud caused by the Deep Impact experiment

We present an analysis of the observations of the Deep Impact event performed by the OSIRIS narrow angle camera aboard the Rosetta spacecraft over two weeks, in an effort to characterize the cometary dust grains ejected from the nucleus of Comet 9P/Tempel 1. We adopt a Monte Carlo approach to generate calibrated synthetic images, and a linear combination of them is fitted to the calibrated images so as to determine the physical parameters of the dust cloud. Our model considers spherical olivine particles with a density of 3780 kg m⁻³. It incorporates constraints on the direction of the cone of emission coming from additional images obtained at Pic du Midi observatory, and constraints on the dust terminal velocities coming from the physics of the impact. We find that the slope of the differential dust size distribution of grains with radii <20 μm (β>0.008) is 3.1±0.3, a value typical of cometary dust tails. This shows that there is no evidence in our data for an enhancement in sub-micron particles in the ejecta compared to the typical dust distribution of active comets. We estimate the mass of particles with radii <1.4 μm (β>0.14) to be 1.5±0.2×¹⁰⁵ kg. These particles represent more than 80% of the cross-section of the observed dust cloud. The mass carried by larger particles depends whether the gas significantly increases the kinetic energy of the grains in the inner coma; it lies in the range 1-14×¹⁰⁶ kg for particles with radii <100 μm (β>0.002). We obtain the distribution of terminal velocities reached by the dust after the dust-gas interaction which is very well constrained between 10 and 600 m s⁻¹. It is characterized by Gaussian with a maximum at about 190 m s⁻¹ and a width at half maximum of 150 m s⁻¹.     

UBV photometry of asteroid 433 Eros

UBV observations of asteroid 433 Eros were conducted on 17 nights during the winter of 1974/75. The peak-to-peak amplitude of the lightcurve varied from about 0.3 mag to nearly 1.4mmag. The absolute V mag at maximum light, extrapolated to zero phase, is 10.85. Phase coefficients of 0.0233 mag/degree, 0.0009 mag/degree and 0.0004 mag/degree were derived for V, B-V, and U-B, respectively. The zero-phase color of Eros (B?V = 0.88, U?B = 0.50) is representative of an S (silicaceous) compositional type asteroid. The color does not vary with rotation. The photometric behavior of Eros can be modeled by a cylinder with rounded ends having an axial ratio of about 2.3:1. The asteroid is rotating about a short axis with the north pole at λ₀ = 15° and β₀ = 9°.     

Revisiting the thermal inertia of Iapetus: Clues to the thickness of the dark material

The energy balance at the surface of an airless planetary body is strongly influenced by the bolometric Bond albedo and the surface thermal inertia. Both of these values may be calculated through the application of a thermal model to measured surface temperatures. The accuracy of either, though, increases if the value of the other is better constrained. In this study, we used the improved global bolometric Bond albedo map of Iapetus derived from Cassini VIMS and ISS and Voyager ISS data in conjunction with Cassini CIRS temperature data to reevaluate surface thermal inertia across Iapetus. Results showed the thermal inertia of the dark terrain varies between 11 and 14.8 J m⁻² K⁻¹ s^−1/2 while the light material varies between 15 and 25 J m⁻² K⁻¹ s^−1/2. Using an approximation to the thermal properties of the dark overburden derived from our thermal inertia results, we can implement our thermal model to provide estimates on the dark material thickness, which was found to lie between 7 cm and 16 cm. In order to develop an accurate global thermal model, a weighted function that approximates the surface thermal inertia across Iapetus was developed and verified via our measurements. The global bolometric Bond albedo map, surface thermal inertia map, and the thermal model are then used to synthesize global temperature maps that may be used to study the stability of volatiles.     

Physical properties of morphological units on Comet 9P/Tempel 1 derived from near-IR Deep Impact spectra

In this paper we analyze near-infrared thermal emission spectra of the spatially resolved nucleus of Comet 9P/Tempel 1 obtained by the NASA spacecraft Deep Impact. Maps of spectral reddening, the product X^′ between the beaming function and directional emissivity, as well as surface temperature are constructed. Thermophysical modeling is used to estimate the degree of small scale surface roughness and thermal inertia by detailed reproduction of the empirical temperature map. Mie and Hapke theories are used in combination with numerically calculated beaming functions to analyze the X^′ map and place constraints on composition and grain size of the surface material. We show that it is absolutely mandatory to include small scale surface roughness in thermophysical modeling of this object, since the resulting self heating is vital for reproducing the measured temperatures. A small scale self heating parameter in the range 0.6?ξ?0.75 is common, but smoother areas where 0.2?ξ?0.3 are also found. Contrary to models neglecting small scale surface roughness, we find that the thermal inertia of Comet 9P/Tempel 1 generally is high (1000-3000 J m⁻² K⁻¹ s^−1/2), although it may be substantially lower (40-380 J m⁻² K⁻¹ s^−1/2) in specific areas. We obtain a disk-averaged reddening of 3.5% kÅ⁻¹, with statistically significant local variations around that value on a ±1.0% kÅ−1 level. Vast regions appear covered by small (∼0.1 μm) highly absorbing grains such as carbon or iron-rich silicates. Other regions appear dominated by somewhat larger (∼0.5 μm) and/or less absorbing grains such as troilite or magnesium-rich silicates. Surface variations in reddening, roughness, thermal inertia, composition and/or grain size are moderately to strongly correlated to the locations of morphological units on the surface. The existence of morphological units with differing physical properties may be primordial, hence reflecting a diversity in the building block cometesimals, or resulting from evolutionary processes.     

Stability and regions of motion in the restricted three-body problem when both the primaries are finite straight segments

In this paper we have studied the locations and stability of the Lagrangian equilibrium points in the restricted three-body problem under the assumption that both the primaries are finite straight segments. We have found that the triangular equilibrium points are conditional stable for 0<μ<μ _c , and unstable in the range μ _c <μ≤1/2, where μ is the mass ratio. The critical mass ratio μ _c depends on the lengths of the segments and it is observed that the range of μ _c increases when compared with the classical case. The collinear equilibrium points are unstable for all values of μ. We have also studied the regions of motion of the infinitesimal mass. It has been observed that the Jacobian constant decreases when compared with the classical restricted three-body problem for a fixed value of μ and lengths l ₁ and l ₂ of the segments. Beside this we have found the numerical values for the position of the collinear and triangular equilibrium points in the case of some asteroids systems: (i) 216 Kleopatra-951 Gaspara, (ii) 9 Metis-433 Eros, (iii) 22 Kalliope-243 Ida and checked the linear stability of stationary solutions of these asteroids systems.     

Carbon Monoxide on Jupiter: Evidence for Both Internal and External Sources

Thirteen lines of the CO band near 4.7 μm have been observed on a jovian hot spot at a resolution of 0.045 cm⁻¹. The measured line profiles indicate that the CO mole fraction is 1.0±0.2 ppb around the 6-bar level and is larger in the upper troposphere and/or stratosphere. An external source of CO providing an abundance of 4⁺³₋₂×10¹⁶ molecules cm⁻² is implied by the observations in addition to the amount deposited at high altitude by the Shoemaker-Levy 9 collision. From a simple diffusion model, we estimate that the CO production rate is (1.5-10)×10⁶ molecules cm⁻² s⁻¹ assuming an eddy diffusion coefficient around the tropopause between 300 and 1500 cm² s⁻¹. Precipitation of oxygen atoms from the jovian magnetosphere or photochemistry of water vapor from meteoroidal material can only provide a negligible contribution to this amount. A significant fraction of the CO in Jupiter's upper atmosphere may be formed by shock chemistry due to the infall of kilometer- to subkilometer-size Jupiter family comets. Using the impact rate from Levison et al. (2000, Icarus143, 415-420) rescaled by Bottke et al. (2002, Icarus156, 399-433), this source can provide the observed stratospheric CO only if the eddy diffusion coefficient around the tropopause is 100-300 cm² s⁻¹. Higher values, ∼700 cm² s⁻¹, would require an impact rate larger by a factor of 5-10, which cannot be excluded considering uncertainties in the distribution of Jupiter family comets. Such a large rate is indeed consistent with the observed cratering record of the Galilean satellites (Zahnle et al. 1998, Icarus136, 202-222). On the other hand, the ∼1 ppb concentration in the lower troposphere requires an internal source. Revisiting the disequilibrium chemistry of CO in Jupiter, we conclude that rapid vertical mixing can provide the required amount of CO at ∼6 bar for a global oxygen abundance of 0.2-9 times the solar value considering the uncertainties in the convective mixing rate and in the chemical constants.