The E ring in the vicinity of Enceladus: I. Spatial distribution and properties of the ring particles

Saturn's diffuse E ring is the largest ring of the Solar System and extends from about (Saturn radius RS=60,330 km) to at least encompassing the icy moons Mimas, Enceladus, Tethys, Dione, and Rhea. After Cassini's insertion into her saturnian orbit in July 2004, the spacecraft performed a number of equatorial as well as steep traversals through the E ring inside the orbit of the icy moon Dione. Here, we report about dust impact data we obtained during 2 shallow and 6 steep crossings of the orbit of the dominant ring source—the ice moon Enceladus. Based on impact data of grains exceeding 0.9 μm we conclude that Enceladus feeds a torus populated by grains of at least this size along its orbit. The vertical ring structure at agrees well with a Gaussian with a full-width-half-maximum (FWHM) of ∼4200 km. We show that the FWHM at is due to three-body interactions of dust grains ejected by Enceladus' recently discovered ice volcanoes with the moon during their first orbit. We find that particles with initial speeds between 225 and 235 m s⁻¹ relative to the moon's surface dominate the vertical distribution of dust. Particles with initial velocities exceeding the moon's escape speed of 207 m s⁻¹ but slower than 225 m s⁻¹ re-collide with Enceladus and do not contribute to the ring particle population. We find the peak number density to range between 16×¹⁰⁻² m⁻³ and 21×¹⁰⁻² m⁻³ for grains larger 0.9 μm, and 2.1×¹⁰⁻² m⁻³ and 7.6×¹⁰⁻² m⁻³ for grains larger than 1.6 μm. Our data imply that the densest point is displaced outwards by at least with respect of the Enceladus orbit. This finding provides direct evidence for plume particles dragged outwards by the ambient plasma. The differential size distribution for grains >0.9 μm is described best by a power law with slopes between 4 and 5. We also obtained dust data during ring plane crossings in the vicinity of the orbits of Mimas and Tethys. The vertical distribution of grains >0.8 μm at Mimas orbit is also well described by Gaussian with a FWHM of ∼5400 km and displaced southwards by ∼1200 km with respect to the geometrical equator. The vertical distribution of ring particles in the vicinity of Tethys, however, does not match a Gaussian. We use the FWHM values obtained from the vertical crossings to establish a 2-dimensional model for the ring particle distribution which matches our observations during vertical and equatorial traversals through the E ring.     

Gravity field and interior of Rhea from Cassini data analysis

The Cassini spacecraft encountered Rhea on November 26, 2005. Analysis of the Doppler data acquired at and around closest approach yields the mass of Rhea and the quadrupole moments of its gravity field with unprecedented accuracy. We obtained which corresponds to a density of . Our results for J₂ and C₂₂ are (7.947±0.892)×¹⁰⁻⁴ and (2.3526±0.0476)×¹⁰⁻⁴, respectively. These values are consistent with hydrostatic equilibrium. From the value of C₂₂, we infer the non-dimensional moment of inertia C/MR²=0.3721±0.0036. Our models of Rhea's interior based on the gravity data favor an almost undifferentiated satellite. A discontinuity between a core and a mantle is possible but not required by the data. Models with a constant silicate mass fraction throughout the body cannot account for the determined quadrupole coefficients. The data exclude fully differentiated models in which the core would be composed of unhydrated silicates and the mantle would be composed of pure ice. If the mantle contains 10% in mass of silicates, the core extends to 630 km in radius and has a silicate mass fraction of 40%. A continuous model in which the silicates are more concentrated toward the center of the body than in the outer layers is allowed by the gravity data but excluded by thermal evolution considerations. The one model that fits the gravity data and is self-consistent when energy transport and ice melting are qualitatively considered is an “almost undifferentiated” Rhea, in which a very large uniform core is surrounded by a relatively thin ice shell containing no rock at all.     

A ballistics analysis of the Deep Impact ejecta plume: Determining Comet Tempel 1's gravity, mass, and density

In July of 2005, the Deep Impact mission collided a 366 kg impactor with the nucleus of Comet 9P/Tempel 1, at a closing speed of 10.2 km s⁻¹. In this work, we develop a first-order, three-dimensional, forward model of the ejecta plume behavior resulting from this cratering event, and then adjust the model parameters to match the flyby-spacecraft observations of the actual ejecta plume, image by image. This modeling exercise indicates Deep Impact to have been a reasonably “well-behaved” oblique impact, in which the impactor-spacecraft apparently struck a small, westward-facing slope of roughly 1/3-1/2 the size of the final crater produced (determined from initial ejecta plume geometry), and possessing an effective strength of not more than . The resulting ejecta plume followed well-established scaling relationships for cratering in a medium-to-high porosity target, consistent with a transient crater of not more than 85-140 m diameter, formed in not more than 250-550 s, for the case of (gravity-dominated cratering); and not less than 22-26 m diameter, formed in not less than 1-3 s, for the case of (strength-dominated cratering). At , an upper limit to the total ejected mass of 1.8×¹⁰⁷ kg (1.5-2.2×¹⁰⁷ kg) is consistent with measurements made via long-range remote sensing, after taking into account that 90% of this mass would have stayed close to the surface and then landed within 45 min of the impact. However, at , a lower limit to the total ejected mass of 2.3×¹⁰⁵ kg (1.5-2.9×¹⁰⁵ kg) is also consistent with these measurements. The expansion rate of the ejecta plume imaged during the look-back phase of observations leads to an estimate of the comet's mean surface gravity of (0.17-0.90 mm s⁻²), which corresponds to a comet mass of mt=4.5×¹⁰¹³ kg (2.3-12.0×¹⁰¹³ kg) and a bulk density of (200-1000 kg m⁻³), where the large high-end error is due to uncertainties in the magnitude of coma gas pressure effects on the ejecta particles in flight.     

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.     

Submillimeter Wave Astronomy Satellite observations of Comet 9P/Tempel 1 and Deep Impact

On 4 July 2005 at 5:52 UT the Deep Impact mission successfully completed its goal to hit the nucleus of 9P/Tempel 1 with an impactor, forming a crater on the nucleus and ejecting material into the coma of the comet. NASA's Submillimeter Wave Astronomy Satellite (SWAS) observed the ₁₁₀-₁₀₁ ortho-water ground-state rotational transition in Comet 9P/Tempel 1 before, during, and after the impact. No excess emission from the impact was detected by SWAS and we derive an upper limit of 1.8×¹⁰⁷ kg on the water ice evaporated by the impact. However, the water production rate of the comet showed large natural variations of more than a factor of three during the weeks before and after the impact. Episodes of increased activity with alternated with periods with low outgassing (). We estimate that 9P/Tempel 1 vaporized a total of N∼4.5×¹⁰³⁴ water molecules (∼1.3×¹⁰⁹ kg) during June-September 2005. Our observations indicate that only a small fraction of the nucleus of Tempel 1 appears to be covered with active areas. Water vapor is expected to emanate predominantly from topographic features periodically facing the Sun as the comet rotates. We calculate that appreciable asymmetries of these features could lead to a spin-down or spin-up of the nucleus at observable rates.     

Meridional variations of temperature, C2H2 and C2H6 abundances in Saturn's stratosphere at southern summer solstice

Measurements of the vertical and latitudinal variations of temperature and C₂H₂ and C₂H₆ abundances in the stratosphere of Saturn can be used as stringent constraints on seasonal climate models, photochemical models, and dynamics. The summertime photochemical loss timescale for C₂H₆ in Saturn's middle and lower stratosphere (∼40-10,000 years, depending on altitude and latitude) is much greater than the atmospheric transport timescale; ethane observations may therefore be used to trace stratospheric dynamics. The shorter chemical lifetime for C₂H₂ (∼1-7 years depending on altitude and latitude) makes the acetylene abundance less sensitive to transport effects and more sensitive to insolation and seasonal effects. To obtain information on the temperature and hydrocarbon abundance distributions in Saturn's stratosphere, high-resolution spectral observations were obtained on September 13-14, 2002 UT at NASA's IRTF using the mid-infrared TEXES grating spectrograph. At the time of the observations, Saturn was at a LS≈270°, corresponding to Saturn's southern summer solstice. The observed spectra exhibit a strong increase in the strength of methane emission at 1230 cm⁻¹ with increasing southern latitude. Line-by-line radiative transfer calculations indicate that a temperature increase in the stratosphere of ≈10 K from the equator to the south pole between 10 and 0.01 mbar is implied. Similar observations of acetylene and ethane were also recorded. We find the 1.16 mbar mixing ratio of C₂H₂ at −1° and −83° planetocentric latitude to be and , respectively. The C₂H₂ mixing ratio at 0.12 mbar is found to be at −1° planetocentric latitude and at −83° planetocentric latitude. The 2.3 mbar mixing ratio of C₂H₆ inferred from the data is and at −1° and −83° planetocentric latitude, respectively. Further observations, creating a time baseline, will be required to completely resolve the question of how much the latitudinal variations of C₂H₂ and C₂H₆ are affected by seasonal forcing and/or stratospheric circulation.     

Heavy ion escape from Mars, influence from solar wind conditions and crustal magnetic fields

We have used more than 4 years of Mars Express ion data to estimate the escape of heavy ions ( and ) from Mars. To take the limited field of view of the instrument into account, the data has been binned into spatial bins and angular bins to create average distribution functions for different positions in the near Mars space. The net escape flux for the studied low solar activity period, between May 2007 and May 2011, is 2.0 ± 0.2 × 10²⁴ s⁻¹. The escape has been calculated independently for four different quadrants in the Y_MSO − Z_MSO plane, south, dusk, north and dawn. Escape is highest from the northern and dusk quadrants, 0.6 ± 0.1 × 10²⁴ s⁻¹, and smallest from the south and dawn quadrants, 0.4 ± 0.1 × 10²⁴ s⁻¹. The flux ratio of molecular ( and ) to O⁺ ions is 0.9 ± 0.1, averaged over all quadrants. The flux difference between the north and south quadrants is statistically significant, and is presumed to be due to the presence of significant crustal magnetic fields in the southern hemisphere, reducing the outflow. The difference between the dawn and dusk quadrants is likely due to the magnetic tension associated with the nominal Parker angle spiral, which should lead to higher average magnetic tension on the dusk side. The escape increases during periods of high solar wind flux and during times when co-rotating interaction regions (CIR) affect Mars. In the latter case the increase is a factor 2.4-2.9 as compared to average conditions.     

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.     

Mars high resolution gravity fields from MRO, Mars seasonal gravity, and other dynamical parameters

With 2 years of tracking data collection from the MRO spacecraft, there is noticeable improvement in the high frequency portion of the spherical harmonic Mars gravity field. The new JPL Mars gravity fields, MRO110B and MRO110B2, show resolution near degree 90. Additional years of MGS and Mars Odyssey tracking data result in improvement for the seasonal gravity changes which compares well to global circulation models and Odyssey neutron data and Mars rotation and precession (). Once atmospheric dust is accounted for in the spacecraft solar pressure model, solutions for Mars solar tide are consistent between data sets and show slightly larger values (k₂ = 0.164 ± 0.009, after correction for atmospheric tide) compared to previous results, further constraining core models. An additional 4 years of Mars range data improves the Mars ephemeris, determines 21 asteroid masses and bounds solar mass loss (dGM_Sun/dt < 1.6 × 10⁻¹³ GM_Sun year⁻¹).     

Detection of C2HD and the D/H ratio on Titan

We report here the first detection of mono-deuterated acetylene (acetylene-d1, C₂HD) in Titan's atmosphere from the presence of two of its emission bands at 678 and 519 cm⁻¹ as observed in CIRS spectral averages of nadir and limb observations taken between July 2004 and mid-2007. By using new laboratory spectra for this molecule, we were able to derive its abundance at different locations over Titan's disk. We find the C₂HD value () to be roughly constant with latitude from the South to about 45° N and then to increase slightly in the North, as is the case for C₂H₂. Fitting the 678 cm⁻¹ν₅ band simultaneously with the nearby C₂H₂ 729 cm⁻¹ν₅ band, allows us to infer a D/H ratio in acetylene on Titan with an average of the modal values of 2.09±0.45×¹⁰⁻⁴ from the nadir observations, the uncertainties being mainly due to the vertical profile used for the fit of the acetylene band. Although still subject to significant uncertainty, this D/H ratio appears to be significantly larger than the one derived in methane from the CH₃D band (upper limit of 1.5×¹⁰⁻⁴; Bézard, B., Nixon, C.A., Kleiner, I., Jennings, D.E., 2007. Icarus, 191, 397-400; Coustenis, A., Achterberg, R., Conrath, B., Jennings, D., Marten, A., Gautier, D., Bjoraker, G., Nixon, C., Romani, P., Carlson, R., Flasar, M., Samuelson, R.E., Teanby, N., Irwin, P., Bézard, B., Orton, G., Kunde, V., Abbas, M., Courtin, R., Fouchet, Th., Hubert, A., Lellouch, E., Mondellini, J., Taylor, F.W., Vinatier, S., 2007. Icarus 189, 35-62). From the analysis of limb data we infer D/H values of (at 54° S), (at 15° S), (at 54° N) and (at 80° N), which average to a mean value of 1.63±0.27×¹⁰⁻⁴.     

Jupiter's hydrocarbon polar brightening: Discovery of 3-micron line emission from south polar CH4, C2H2, and C2H6

Jupiter exhibits bright H⁺₃ auroral arcs at 3-4 microns that cool the hot (>1000 K) ionosphere above the ∼¹⁰−7 bar level through the infrared bands of this trace constituent. Below the ¹⁰−7 bar level significant cooling proceeds through infrared active bands of CH₄, C₂H₂, and C₂H₆. We report the discovery of 3-micron line emission from these hydrocarbon species in spectra of the jovian south polar region obtained on April 18 and 20, 2006 (UT) with CGS4 on the United Kingdom Infrared Telescope. Estimated cooling rates through these molecules are 7.5×¹⁰−3, 1.4×¹⁰−3, and , respectively, for a total nearly half that of H⁺₃. We derive a temperature of 450 ± 50 K in the ¹⁰−7-¹⁰−5 bar region from the C₂H₂ lines.     

Prebiotic organic synthesis in early Earth and Mars atmospheres: Laboratory experiments with quantitative determination of products formed in a cold plasma flow reactor

The goal of this study was to explore prebiotic chemistry in a range of plausible early Earth and Mars atmospheres. To achieve this laboratory continuous flow plasma irradiation experiments were performed on N₂/H₂/CO/CO₂ gas mixtures chosen to represent mildly reducing early Earth and Mars atmospheres derived from a secondary volcanic outgassing of volatiles in chemical equilibrium with magmas near present day oxidation state. Under mildly reducing conditions (91.79% N₂, 5.89% H₂, 2.21% CO, and 0.11% CO₂), simple nitriles are produced in the gas phase with yield (G in molecules per 100 eV), for the key prebiotic marker molecule HCN at G∼1×¹⁰⁻³ (0.1 nmol J⁻¹). In this atmosphere localized HCN concentrations possibly could approach the 10⁻² M needed for HCN oligomerization. Yields under mildly oxidizing conditions (45.5% N₂, 0.1% H₂, 27.2% CO, 27.2% CO₂) are significantly less as expected, with HCN at G∼3×¹⁰⁻⁵ (). Yields in a Triton atmosphere which can be plausibly extrapolated to represent what might be produced in trace CH₄ conditions (99.9% N₂, 0.1% CH₄) are significant with HCN at G∼1×¹⁰⁻² (1 nmol J⁻¹) and tholins produced. Recently higher methane abundance atmospheres have been examined for their greenhouse warming potential, and higher abundance hydrogen atmospheres have been proposed based on a low early Earth exosphere temperature. A reducing (64.04% N₂, 28.8% H₂, 3.60% CO₂, and 3.56% CH₄), representing a high CH₄ and H₂ abundance early Earth atmosphere had HCN yields of G∼5×¹⁰⁻³ (0.5 nmol J⁻¹). Tholins generated in high methane hydrogen gas mixtures is much less than in a similar mixture without hydrogen. The same mixture with the oxidizing component CO₂ removed (66.43% N₂, 29.88% H₂, 0% CO₂, and 3.69% CH₄) had HCN yields of G∼1×¹⁰⁻³ (0.1 nmol J⁻¹) but more significant tholin yields.     

The vertical structure of Titan's upper atmosphere from Cassini Ion Neutral Mass Spectrometer measurements

Data acquired by the Ion Neutral Mass Spectrometer (INMS) on the Cassini spacecraft during its close encounter with Titan on 26 October 2004 reveal the structure of its upper atmosphere. Altitude profiles of N₂, CH₄, and H₂, inferred from INMS measurements, determine the temperature, vertical mixing rate, and escape flux from the upper atmosphere. The mean atmospheric temperature in the region sampled by the INMS is 149±3 K, where the variance is a consequence of local time variations in temperature. The CH₄ mole fraction at 1174 km is 2.71±0.1%. The effects of diffusive separation are clearly seen in the data that we interpret as an eddy diffusion coefficient of , that, along with the measured CH₄ mole fraction, implies a mole fraction in the stratosphere of 2.2±0.2%. The H₂ distribution is affected primarily by upward flow and atmospheric escape. The H₂ mole fraction at 1200 km is 4±1×¹⁰⁻³ and analysis of the altitude profile indicates an upward flux of , referred to the surface. If horizontal variations in temperature and H₂ density are small, this upward flux also represents the escape flux from the atmosphere. The CH₄ density exhibits significant horizontal variations that are likely an indication of dynamical processes in the upper atmosphere.     

A sensitive search for nitric oxide in the lower atmospheres of Venus and Mars: Detection on Venus and upper limit for Mars

High-resolution spectra of Venus and Mars at the NO fundamental band at 5.3 μm with resolving power ν/δν=76,000 were acquired using the TEXES spectrograph at NASA IRTF on Mauna Kea, Hawaii. The observed spectrum of Venus covered three NO lines of the P-branch. One of the lines is strongly contaminated, and the other two lines reveal NO in the lower atmosphere at a detection level of 9 sigma. A simple photochemical model for NO and N at 50-112 km was coupled with a radiative transfer code to simulate the observed equivalent widths of the NO and some CO₂ lines. The derived NO mixing ratio is 5.5±1.5 ppb below 60 km and its flux is . Predissociation of NO at the (0-0) 191 nm and (1-0) 183 nm bands of the δ-system and the reaction with N are the only important loss processes for NO in the lower atmosphere of Venus. The photochemical impact of the measured NO abundance is significant and should be taken into account in photochemical modeling of the Venus atmosphere. Lightning is the only known source of NO in the lower atmosphere of Venus, and the detection of NO is a convincing and independent proof of lightning on Venus. The required flux of NO is corrected for the production of NO and N by the cosmic ray ionization and corresponds to the lightning energy deposition of . For a flash energy on Venus similar to that on the Earth (∼¹⁰⁹ J), the global flashing rate is ∼90 s⁻¹ and ∼6 km⁻² y⁻¹ which is in reasonable agreement with the existing optical observations. The observed spectrum of Mars covered three NO lines of the R-branch. Two of these lines are contaminated by CO₂ lines, and the line at 1900.076 cm⁻¹ is clean and shows some excess over the continuum. Some photochemical reactions may result in a significant excitation of NO (v=1) in the lowest 20 km on Mars. However, quenching of NO (v=1) by CO₂ is very effective below 40 km. Excitation of NO (v=1) in the collisions with atomic oxygen is weak because of the low temperature in the martian atmosphere, and we do not see any explanation of a possible emission of NO at 5.3 μm. Therefore the data are treated as the lack of absorption with a 2 sigma upper limit of 1.7 ppb to the NO abundance in the lower atmosphere of Mars. This limit is above the predictions of photochemical models by a factor of 3.     

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]).     

The orbit, mass, and albedo of transneptunian binary (66652) 1999 RZ253

We have observed (66652) 1999 RZ₂₅₃ with the Hubble Space Telescope at seven separate epochs and have fit an orbit to the observed relative positions of this binary. Two orbital solutions have been identified that differ primarily in the inclination of the orbit plane. The best fit corresponds to an orbital period, days, semimajor axis a=4660±170 km and orbital eccentricity e=0.460±0.013 corresponding to a system mass m=3.7±0.4×¹⁰¹⁸ kg. For a density of the albedo at 477 nm is p₄₇₇=0.12±0.01, significantly higher than has been commonly assumed for objects in the Kuiper belt. Multicolor, multiepoch photometry shows this pair to have colors typical for the Kuiper belt with a spectral gradient of 0.35 per 100 nm in the range between 475 and 775 nm. Photometric variations at the four epochs we observed were as large as 12±3% but the sampling is insufficient to confirm the existence of a lightcurve.     

On the Vertical Thermal Structure of Pluto's Atmosphere

A radiative–conductive model for the vertical thermal structure of Pluto's atmosphere is developed with a non-LTE treatment of solar heating in the CH₄3.3 μm and 2.3 μm bands, non-LTE radiative exchange and cooling in the CH₄7.6 μm band, and LTE cooling by CO rotational line emission. The model includes the effects of opacity and vibrational energy transfer in the CH₄molecule. Partial thermalization of absorbed solar radiation in the CH₄3.3 and 2.3 μm bands by rapid vibrational energy transfer from the stretch modes to the bending modes generates high altitude heating at sub-microbar pressures. Heating in the 2.3 μm bands exceeds heating in 3.3 μm bands by approximately a factor of 6 and occurs predominantly at microbar pressures to generate steep temperature gradients ∼10–20 K km⁻¹forp> 2 μbar when the surface or tropopause pressure is ∼3 μbar and the CH₄mixing ratio is a constant 3%. This calculated structure may account for the “knee” in the stellar occultation lightcurve. The vertical temperature structure in the first 100 km above the surface is similar for atmospheres with Ar, CO, and N₂individually as the major constituent. If a steep temperature gradient ∼20 K km⁻¹is required near the surface or above the tropopause, then the preferred major constituent is Ar with 3% CH₄mixing ratio to attain a calculated ratio ofT/(= 3.5 K amu⁻¹) in agreement with inferred values from stellar occultation data. However, pure Ar and N₂ices at the same temperature yield an Ar vapor pressure of only ∼0.04 times the N₂vapor pressure. Alternative scenarios are discussed that may yield acceptable fits with N₂as the dominant constituent. One possibility is a 3 μbar N₂atmosphere with 0.3% CH₄that has 106 K isothermal region (T/= 3.8 K amu⁻¹) and ∼8 K km⁻¹surface/tropopause temperature gradient. Another possibility would be a higher surface pressure ∼10 μbar with a scattering haze forp> 2 μbar. Our model with appropriate adjustments in the CH₄density profile to Triton's inferred profile yields a temperature profile consistent with the UVS solar occultation data (Krasnopolsky, V. A., B. R. Sandel, and F. Herbert 1992.J. Geophys. Res.98, 3065–3078.) and ground-based stellar occultation data (Elliot, J. L., E. W. Dunham, and C. B. Olkin 1993.Bull. Am. Astron. Soc.25, 1106.).     

Atomic oxygen distribution in the Venus mesosphere from observations of O2 infrared airglow by VIRTIS-Venus Express

This VIRTIS instrument on board Venus Express has collected spectrally resolved images of the Venus nightside limb that show the presence of the (0,0) band of the infrared atmospheric system of O₂ at 1.27 μm. The emission is produced by three-body recombination of oxygen atoms created by photodissociation of CO₂ on the dayside. It is consistently bright so that emission limb profiles can be extracted from the images. The vertical distribution of O₂() may be derived following Abel inversion of the radiance limb profiles. Assuming photochemical equilibrium, it is combined with the CO₂ vertical distribution to determine the atomic oxygen density. The uncertainties on the O density caused by the Abel inversion reach a few percent at the peak, increasing to about 50% near 120 km. We first analyze a case when the CO₂ density was derived from a stellar occultation observed with the SPICAV spectrometer simultaneously with an image of the O₂ limb airglow. In other cases, an average CO₂ profile deduced from a series of ultraviolet stellar occultations is used to derive the O profile, leading to uncertainties on the O density less than 30%. It is found that the maximum O density is generally located between 94 and 115 km with a mean value of 104 km. It ranges from less than 1×¹⁰¹¹ to about 5×¹⁰¹¹ cm⁻³ with a global mean of 2.2×¹⁰¹¹ cm⁻³. These values are in reasonable agreement with the VIRA midnight oxygen profile. The vertical O distribution is generally in good agreement with the oxygen profile calculated with a one-dimensional chemical-diffusive model. No statistical latitudinal dependence of the altitude of the oxygen peak is observed, but the maximum O density tends to decrease with increasing northern latitudes. The latitudinal distribution at a given time exhibits large variations in the O density profile and its vertical structure. The vertical oxygen distribution frequently shows multiple peaks possibly caused by waves or variations in the structure of turbulent transport. It is concluded that the O₂ infrared night airglow is a powerful tool to map the distribution of atomic oxygen in the mesosphere between 90 and 115 km and improve future Venus reference atmosphere models.     

Accretion of the gaseous envelope of Jupiter around a 5-10 Earth-mass core

New numerical simulations of the formation and evolution of Jupiter are presented. The formation model assumes that first a solid core of several M_⊕ accretes from the planetesimals in the protoplanetary disk, and then the core captures a massive gaseous envelope from the protoplanetary disk. Earlier studies of the core accretion-gas capture model [Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y., 1996. Icarus 124, 62-85] demonstrated that it was possible for Jupiter to accrete with a solid core of 10-30 M_⊕ in a total formation time comparable to the observed lifetime of protoplanetary disks. Recent interior models of Jupiter and Saturn that agree with all observational constraints suggest that Jupiter's core mass is 0-11 M_⊕ and Saturn's is 9-22 M_⊕ [Saumon, G., Guillot, T., 2004. Astrophys. J. 609, 1170-1180]. We have computed simulations of the growth of Jupiter using various values for the opacity produced by grains in the protoplanet's atmosphere and for the initial planetesimal surface density, σinit,Z, in the protoplanetary disk. We also explore the implications of halting the solid accretion at selected core mass values during the protoplanet's growth. Halting planetesimal accretion at low core mass simulates the presence of a competing embryo, and decreasing the atmospheric opacity due to grains emulates the settling and coagulation of grains within the protoplanet's atmosphere. We examine the effects of adjusting these parameters to determine whether or not gas runaway can occur for small mass cores on a reasonable timescale. We compute four series of simulations with the latest version of our code, which contains updated equation of state and opacity tables as well as other improvements. Each series consists of a run without a cutoff in planetesimal accretion, plus up to three runs with a cutoff at a particular core mass. The first series of runs is computed with an atmospheric opacity due to grains (hereafter referred to as ‘grain opacity’) that is 2% of the interstellar value and . Cutoff runs are computed for core masses of 10, 5, and 3 M_⊕. The second series of Jupiter models is computed with the grain opacity at the full interstellar value and . Cutoff runs are computed for core masses of 10 and 5 M_⊕. The third series of runs is computed with the grain opacity at 2% of the interstellar value and . One cutoff run is computed with a core mass of 5 M_⊕. The final series consists of one run, without a cutoff, which is computed with a temperature dependent grain opacity (i.e., 2% of the interstellar value for ramping up to the full interstellar value for ) and . Our results demonstrate that reducing grain opacities results in formation times less than half of those for models computed with full interstellar grain opacity values. The reduction of opacity due to grains in the upper portion of the envelope with has the largest effect on the lowering of the formation time. If the accretion of planetesimals is not cut off prior to the accretion of gas, then decreasing the surface density of planetesimals lowers the final core mass of the protoplanet, but increases the formation timescale considerably. Finally, a core mass cutoff results in a reduction of the time needed for a protoplanet to evolve to the stage of runaway gas accretion, provided the cutoff mass is sufficiently large. The overall results indicate that, with reasonable parameters, it is possible that Jupiter formed at 5 AU via the core accretion process in 1 Myr with a core of 10 M_⊕ or in 5 Myr with a core of 5 M_⊕.