On the mean oxygen isotope composition of the Solar System

Since the first discovery of extraordinary oxygen isotope compositions in carbonaceous meteorites by Clayton et al. [Clayton, R.N., Grossman, L., Mayeda, T.K., 1973. Science 182, 485-488], numerous studies have been done to explain the unusual mass-independent isotope fractionation, but the problem is still unresolved to this day. Clayton's latest interpretation [Clayton, R.N., 2002. Nature 415, 860-861] sheds new light on the problem, and possible hypotheses now seem to be fairly well defined. A key issue is to resolve whether the oxygen isotopes in the Solar System represented by the Sun (solar oxygen) are the same as oxygen isotopes in planetary objects such as bulk meteorites, Mars, Earth, and Moon, or whether the solar oxygen is more similar to the lightest oxygen isotopes observed in CAIs (Calcium Aluminum-rich Inclusions) in primitive meteorites. Here, we examined the problem using oxygen isotope analytical data of about 400 bulk meteorite samples of various classes or types (data compiled by K. Lodders). We used in our discussion exclusively the parameter , a direct measure of the degree of mass-independent isotope fractionation of oxygen isotopes. When is arranged according to a characteristic size of their host planetary object, it shows a systematic trend: (1) values scatter around zero; (2) the scatter from the mean () decreases with increasing representative size of the respective host planetary object. This systematic trend is easily understood on the basis of a hierarchical scenario of planetary formation, that is, larger planetary objects have formed by progressive accretion of planetesimals by random sampling over a wide spectrum of proto-solar materials. If this progressive random sampling of planetesimals were the essential process of planetary formation, the isotopic composition of planetary oxygen should approach that of the solar oxygen. To test this random sampling hypothesis, we applied a multiscale, multistep bootstrap statistical method [Shimodaira, H., 2004. Ann. Statist. 32, 2616-2641] to the meteorite oxygen isotope data, and deduced a σ-N relation, where σ is the standard deviation of , and N is the representative size of a host planetary object. If we assign 200 and 500 km as a representative sizes of the chondrite and achondrite parent bodies, the observed σ of agree well with the values predicted by the σ-N relation. A common mean value of for all planetary objects also agrees with the progressive random sampling process. Therefore, we conclude that the solar oxygen is the same as planetary oxygen, but differs from CAI oxygen. The conclusion implies that a massive enrichment in ¹⁷O and ¹⁸O resulting from CO self-shielding, a current influential interpretation of CAI-O, did not occur.     

Compositional and physical results for Rosetta's new target Comet 67P/Churyumov-Gerasimenko from narrowband photometry and imaging

We present compositional and physical results of Comet 67P/Churyumov-Gerasimenko, the new target of ESA's Rosetta mission. A total of 16 nights of narrowband photometry were obtained at Lowell Observatory during the 1982/83 and 1995/96 apparitions, along with one night of imaging near perihelion in 1996. These data encompass an interval of −61 to +118 days from perihelion, corresponding to a range of heliocentric distances before perihelion from 1.48 to 1.34 AU, and an outbound range from 1.30 to 1.86 AU. Production rates were determined for OH, NH, CN, C₃, and C₂, along with A(θ)fρ, a proxy of the dust production. Water production, based on OH, has a steep () power-law rH-dependence post-perihelion and the minor species are somewhat less steep ( to −4), while the dust is quite shallow (), possibly due to a lingering population of large, slow-moving grains. All species exhibit larger production rates after perihelion, with water having a ∼2×pre/post-perihelion asymmetry, while minor species and dust have larger asymmetries. These asymmetries imply a strong seasonal effect and probable high obliquity of the rotational axis, along with one or more isolated source regions coming into sunlight near perihelion. Peak water production (which occurred about 1 month after perihelion) was and, when combined with a standard water vaporization model, implies an effective active area on the surface of the nucleus of ∼1.5-2.2 km² or an active fraction of only about 3-4%. Abundances of carbon-chain molecules yield a classification of slightly “depleted” in the A'Hearn et al. [A'Hearn, M.F., Millis, R.L., Schleicher, D.G., Osip, D.J., Birch, P.V., 1995. Icarus 118, 223-270] database. The peak dust production (as measured by A(θ)fρ, and uncorrected for phase angle) was ∼450 cm, while the color of the dust is moderately reddened, and the mean radial profile has a power-law slope of −1.3. Large night-to-night variability is also present, presumably due to the source region(s) rotating in and out of sunlight along with effects due to the use of differently sized apertures. A strong sunward radial feature was detected in images obtained near perihelion, along with a significant asymmetry between the two perpendicular directions from the Sun/tail line. These features may be the result of a mid-latitude source region sweeping out a cone with each rotation, which we are viewing from the side and where the sunward radial feature is one edge of the cone seen in projection. When combined with other constraints on the pole orientation, a possible pole solution is found having an obliquity of about 134° at an RA of about 223° and a Dec of −65°, with a source region located near +50° and in overall agreement with the photometric results. In comparison to the original Rosetta target Comet 46P/Wirtanen, Comet Churyumov-Gerasimenko has essentially the same peak water production but a peak dust production about 3 times greater than does Wirtanen based on A(θ)fρ (i.e., if one assumes that the properties of the dust grains are similar) (cf. Farnham and Schleicher [1998. Astron. Astrophys. 335, L50-L55]).     

Parametric analysis of modeled ion escape from Mars

We develop a parametric fit to the results of a detailed magnetohydrodynamic (MHD) study of the response of ion escape rates (O⁺, and ) to strongly varied solar forcing factors, as a way to efficiently extend the MHD results to different conditions. We then use this to develop a second, evolutionary model of solar forced ion escape. We treat the escape fluxes of ion species at Mars as proportional to the product of power laws of four factors - that of the EUV flux R_euv, the solar wind particle density R_ρ, its velocity (squared) R_v², and the interplanetary magnetic field pressure R_B², where forcing factors are expressed in units of the current epoch-averaged values. Our parametric model is: , where ?(i) is the escape flux of ion i. We base our study on the results of just six provided MHD model runs employing large forcing factor variations, and thus construct a successful, first-order parametric model of the MHD program. We perform a five-dimensional least squares fit of this power law model to the MHD results to derive the flux normalizations and the indices of the solar forcing factors. For O⁺, we obtain the values, 1.73 × 10²⁴ s⁻¹, 0.782, 0.251, 0.382, and 0.214, for ?₀, α, β, γ, and δ, respectively. For , the corresponding values are 1.68 × 10²⁴ s⁻¹, −0.393, 0.798, 0.967, and 0.533. For , they are 8.66 × 10²² s⁻¹, −0.427, 1.083, 1.214, and 0.690. The fit reproduces the MHD results to an average error of about 5%, suggesting that the power laws are broadly representative of the MHD model results. Our analysis of the MHD model shows that by itself an increase in R_EUV enhances O⁺ loss, but suppresses the escape of and , whereas increases in solar wind (i.e., in , and R_B², with R_euv constant) favors the escape of heavier ions more than light ions. The ratios of escaping ions detectable at Mars today can be predicted by this parametric fit as a function of the solar forcing factors. We also use the parametric model to compute escape rates over martian history. This second parametric model expresses ion escape functions of one variable (per ion), ?(i) = ?₀(i)(t/t₀)^−ξ(i). The ξ(i) are linear combinations of the epoch-averaged ion escape sensitivities, which are seen to increase with ion mass. We integrate the and oxygen ion escape rates over time, and find that in the last 3.85 Gyr, Mars would have lost about mbars of , and of water (from O⁺ and ) from ion escape.     

Titania's radius and an upper limit on its atmosphere from the September 8, 2001 stellar occultation

On September 8, 2001 around 2 h UT, the largest uranian moon, Titania, occulted Hipparcos star 106829 (alias SAO 164538, a V=7.2, K0 III star). This was the first-ever observed occultation by this satellite, a rare event as Titania subtends only 0.11 arcsec on the sky. The star's unusual brightness allowed many observers, both amateurs or professionals, to monitor this unique event, providing fifty-seven occultations chords over three continents, all reported here. Selecting the best 27 occultation chords, and assuming a circular limb, we derive Titania's radius: (1-σ error bar). This implies a density of using the value derived by Taylor [Taylor, D.B., 1998. Astron. Astrophys. 330, 362-374]. We do not detect any significant difference between equatorial and polar radii, in the limit , in agreement with Voyager limb image retrieval during the 1986 flyby. Titania's offset with respect to the DE405 + URA027 (based on GUST86 theory) ephemeris is derived: ΔαTcos(δT)=−108±13 mas and ΔδT=−62±7 mas (ICRF J2000.0 system). Most of this offset is attributable to a Uranus' barycentric offset with respect to DE405, that we estimate to be: and ΔδU=−85±25 mas at the moment of occultation. This offset is confirmed by another Titania stellar occultation observed on August 1st, 2003, which provides an offset of ΔαTcos(δT)=−127±20 mas and ΔδT=−97±13 mas for the satellite. The combined ingress and egress data do not show any significant hint for atmospheric refraction, allowing us to set surface pressure limits at the level of 10-20 nbar. More specifically, we find an upper limit of 13 nbar (1-σ level) at 70 K and 17 nbar at 80 K, for a putative isothermal CO₂ atmosphere. We also provide an upper limit of 8 nbar for a possible CH₄ atmosphere, and 22 nbar for pure N₂, again at the 1-σ level. We finally constrain the stellar size using the time-resolved star disappearance and reappearance at ingress and egress. We find an angular diameter of 0.54±0.03 mas (corresponding to projected at Titania). With a distance of 170±25 parsecs, this corresponds to a radius of 9.8±0.2 solar radii for HIP 106829, typical of a K0 III giant.     

Fragment-collision model for compound chondrule formation: Estimation of collision probability

We propose a new scenario for compound chondrule formation named as “fragment-collision model,” in the framework of the shock-wave heating model. A molten cm-sized dust particle (parent) is disrupted in the high-velocity gas flow. The extracted fragments (ejectors) are scattered behind the parent and the mutual collisions between them will occur. We modeled the disruption event by analytic considerations in order to estimate the probability of the mutual collisions assuming that all ejectors have the same radius. In the typical case, the molten thin () layer of the parent surface will be stripped by the gas flow. The stripped layer is divided into about 200 molten ejectors (assuming that the radius of ejectors is 300 μm) and then they are blown away by the gas flow in a short period of time (). The stripped layer is leaving from the parent with the velocity of depending on the viscosity, and we assumed that the extracted ejectors have a random velocity Δv of the same order of magnitude. Using above values, we can estimate the number density of ejectors behind the parent as . These ejectors occupy ∼9% of the space behind the parent in volume. Considering that the collision rate (number of collisions per unit time experienced by an ejector) is given by Rcoll=σcollneΔv, where σcoll is the cross-section of collision [e.g., Gooding, J.K., Keil, K., 1981. Meteoritics 16, 17-43], we obtain by substituting above values. Since most collisions occur within the short duration () before the ejectors are blown away, we obtain the collision probability of Pcoll∼0.36, which is the probability of collisions experienced by an ejector in one disruption event. The estimated collision probability is about one order of magnitude larger than the observed fraction of compound chondrules. In addition, the model predictions are qualitatively consistent with other observational data (oxygen isotopic composition, textural types, and size ratios of constituents). Based on these results, we concluded that this new model can be one of the strongest candidates for the compound chondrule formation. It should be noted that all collisions do not necessarily lead to the compound chondrule formation. The formation efficiency and the future works which should be investigated in the forthcoming paper are also discussed.     

A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris

With the collection of six years of MGS tracking data and three years of Mars Odyssey tracking data, there has been a continual improvement in the JPL Mars gravity field determination. This includes the measurement of the seasonal changes in the gravity coefficients (e.g., , , , , , ) caused by the mass exchange between the polar ice caps and atmosphere. This paper describes the latest gravity field MGS95J to degree and order 95. The improvement comes from additional tracking data and the adoption of a more complete Mars orientation model with nutation, instead of the IAU 2000 model. Free wobble of the Mars' spin axis, i.e. polar motion, has been constrained to be less than 10 mas by looking at the temporal history of and . A strong annual signature is observed in , and this is a mixture of polar motion and ice mass redistribution. The Love number solution with a subset of Odyssey tracking data is consistent with the previous liquid outer core determination from MGS tracking data [Yoder et al., 2003. Science 300, 299-303], giving a combined solution of k₂=0.152±0.009 using MGS and Odyssey tracking data. The solutions for the masses of the Mars' moons show consistency between MGS, Odyssey, and Viking data sets; Phobos GM=(7.16±0.005)×¹⁰⁻⁴ km³/s² and Deimos GM=(0.98±0.07)×¹⁰⁻⁴ km³/s². Average MGS orbit errors, determined from differences in the overlaps of orbit solutions, have been reduced to 10-cm in the radial direction and 1.5 m along the spacecraft velocity and normal to the orbit plane. Hence, the ranging to the MGS and Odyssey spacecraft has resulted in position measurements of the Mars system center-of-mass relative to the Earth to an accuracy of one meter, greatly reducing the Mars ephemeris errors by several orders of magnitude, and providing mass estimates for Asteroids 1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, and 324 Bamberga.     

Evolution of the dusty rings of Uranus

We present high quality images of the uranian ring system, obtained in August 2002, October 2003, and July 2004 at 2.2 μm with the adaptive optics camera NIRC2 on the Keck II telescope. Using these data, we report the first detection in backscattered light of a ring (which we refer to as the ζ ring) interior to Uranus' known rings. This ring consists of a generally uniform sheet of dust between 37,850 and 41,350 km with an equivalent width (in 2004; or ), and extends inward to 32,600 km at a gradually decreasing brightness. This ring might be related to the Voyager ring R/1986 U 2, although both its location and extent differ. This could be attributed to a difference in observing wavelength and/or solar phase angle, or perhaps to temporal variations in the ring. Through careful modeling of the I/F of the individual rings at each ansa, we reveal the presence of narrow (few 100 km wide) sheets of dust between the δ and ε rings, and between rings 4 and α. We derived a typical anisotropy factor g≈0.7 in the scattering behavior of these particles. The spatial distribution and relative intensity of these dust sheets is different than that seen in Voyager images taken in forward scattered light, due either to a difference in observing wavelength, and/or solar phase angle or to changes over time. We may have detected the λ ring in one scan at , but other scans provided upper limits below this value. A single detection, however, would be consistent with azimuthal asymmetries known to exist in this ring. We further demonstrate the presence of azimuthal asymmetries in all rings. We confirm the eccentricity of ∼0.001 in rings 4, 5, 6, which in 2004 are ∼70 km closer to Uranus in the north (near periapse; lower I/F) than in the south. We find a global optical depth of τ∼0.3 in the main rings, and of τ=0.25±0.05 in the ε ring.     

Gravitational tidal waves in Titan's upper atmosphere

Tidal waves driven by Titan's orbital eccentricity through the time-dependent component of Saturn's gravitational potential attain nonlinear, saturation amplitudes (|T^′|>10 K, , and ) in the upper atmosphere (?500 km) due to the approximate exponential growth as the inverse square root of pressure. The gravitational tides, with vertical wavelengths of ∼100-150 km above 500 km altitude, carry energy fluxes sufficient in magnitude to affect the energy balance of the upper atmosphere with heating rates in the altitude range of 500-900 km.     

The capture of interstellar dust: the Lorentz force case

The capture of arbitrarily shaped interstellar dust in the Solar System is investigated. Electromagnetic radiation and gravitational forces of the Sun and Lorentz force generated by interplanetary magnetic field are considered. The capture conditions appear to be very sensitive to the particle shape. Non-spherical particles as well as their spherical equivalents are captured only when they are moving initially in the vicinity of ecliptic plane. Capture of non-charged non-spherical dust typically occurs in the region , where R_Sun is solar radius and impact parameter b is defined as the smallest distance between the particle and the Sun if no forces existed. In contrast, charged particles are typically captured at b>150 R_Sun. The total amount of captured non-spherical sub-micron particles differs significantly from the corresponding amount of spherical dust grains. However, both amounts are comparable in the micron-sized range. It is shown that a certain mass of captured non-spherical particles may survive in the Solar System, while captured spherical ones hit the Sun or sublimate in its vicinity. Only a negligible amount of spherical particles can survive. Consideration of solar wind within around of yields that 20% of the captured non-spherical particles of the effective radius survive; the corresponding percentage for particles of the radius is 7%. The total mass of the surviving charged particles is about two orders of magnitude larger than the mass of the surviving non-charged particles. As a result, the sub-micron-sized particles are candidates to contribute to the density increase of the circumsolar dust cloud.