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.     

Stability of ice on Mars and the water vapor diurnal cycle: Experimental study of the sublimation of ice through a fine-grained basaltic regolith

In order to advance our understanding of the long-term stability of subsurface ice, the diurnal martian water cycle, and implications for liquid water, we determined diffusion coefficients and adsorption kinetics for the water vapor produced by the sublimation of ice buried beneath various layers of fine-grained (<63, 63-125, and 125-250 μm) basaltic powder under simulated martian conditions. Sublimation rates at shallower depths, <10 mm, were determined to be affected by mass transfer through the atmosphere in addition to the basalt layer. For greater depths, the measured diffusion coefficients for water vapor moving through basalt grains were 1.56±0.53×¹⁰⁻⁴, 2.05±0.82×¹⁰⁻⁴, and for the <63, 63-125, and 125-250 μm basaltic layers, respectively. Through the Brunauer, Emmett and Teller (BET) isotherm, which assumes multiple molecular layers of adsorbed water, we determined the adsorption constants of 52.6±8.3 at 270 K for <63 μm, 39.0±6.4 at 267 K for 63-125 μm, and 54.3±9.3 at 266 K for 125-250 μm, resulting in surface areas of 2.6±0.1×¹⁰⁴, 1.7±0.3×¹⁰⁴, , respectively. These results suggest that while diffusion is too rapid to explain the purported diurnal cycle in water content of the atmosphere, adsorption is efficient and rapid, and does provide an effective mechanism to explain such a cycle. The present diffusion data suggest that very thin, <50 pr μm, shallow, 10 mm, ice deposits would last for >10 h at ∼224 K, just above the freezing point of saturated CaCl₂. Temperatures can remain above ∼224 K over most of the planet, which means that water, even as saturated brine, will sublimate before the freezing point is reached and liquid could be formed. On the other hand, 1 m ice layers below 1 m of fine-grained basaltic regolith at 235 K and 10 Pa of atmospheric water could last 600 to 1300 years. At deeper depths and lower temperatures, ice could last since the last major obliquity change 400,000 years ago.     

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.     

Observations of C2H6 and C2H2 in the stratosphere of Saturn

We have performed high-resolution spectral observations at mid-infrared wavelengths of C₂H₆ (12.16 μm), and C₂H₂ (13.45 μm) on Saturn. These emission features probe the stratosphere of the planet and provide information on the hydrocarbon photochemical processes taking place in that region of the atmosphere. The observations were performed using our cryogenic echelle spectrometer Celeste, in conjunction with the McMath-Pierce 1.5-m solar telescope in November and December 1994. We used Voyager IRIS CH₄ observations (7.67 μm) to derive a temperature profile on the saturnian atmosphere for the region of the stratosphere. This profile was then used in conjunction with height-dependent volume mixing ratios of each hydrocarbon to determine global abundances for ethane and acetylene. Our ground-based measurements indicate abundances of for C₂H₆ (1.0 mbar pressure level), and for C₂H₂ (1.6 mbar pressure level). We also derived new mixing ratios from the Voyager mid-latitude IRIS observations; 8.6±0.9×¹⁰⁻⁶ for C₂H₆ (0.1-3.0 mbar pressure level), and 1.6±0.2×¹⁰⁻⁷ for C₂H₂ (2.0 mbar pressure level).     

N2 escape rates from Pluto's atmosphere

Hydrodynamic escape of N₂ molecules from Pluto's atmosphere is calculated under the assumption of a high density, slow outflow expansion driven by solar EUV heating by N₂ absorption, near-IR and UV heating by CH₄ absorption, and CO cooling by rotational line emission as a function of solar activity. At 30 AU, the N₂ escape rate varies from in the absence of heating, but driven by an upward thermal heat conduction flux from the stratosphere, for lower boundary temperatures varying from 70-100 K. With solar heating varying from solar minimum to solar maximum conditions and a calculated lower boundary temperature, 88.2 K, the N₂ escape rate range is , respectively. LTE rotational line emission by CO reduces the net solar heat input by at most 35% and plays a minor role in lowering the calculated escape rates, but ensures that the lower boundary temperature can be calculated by radiative equilibrium with near-IR CH₄ heating. While an upward thermal conduction heat flux at the lower boundary plays a fundamental role in the absence of heating, with solar heating it is downward at solar minimum, and is, at most, 13% of the integrated net heating rate over the range of solar activity. For the arrival of the New Horizons spacecraft at Pluto in July 2015, predictions are lower boundary temperature, T₀∼81 K, and N₂ escape rate , and peak thermospheric temperature ∼103 K at 1890 km, based on expected solar medium conditions.     

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.     

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.     

Titan's hydrodynamically escaping atmosphere: Escape rates and the structure of the exobase region

In Strobel [Strobel, D.F., 2008. Icarus, 193, 588-594] a mass loss rate from Titan's upper atmosphere, , was calculated for a single constituent, N₂ atmosphere by hydrodynamic escape as a high density, slow outward expansion driven principally by solar UV heating due to CH₄ absorption. It was estimated, but not proven, that the hydrodynamic mass loss is essentially CH₄ and H₂ escape. Here the individual conservation of momentum equations for the three major components of the upper atmosphere (N₂, CH₄, H₂) are solved in the low Mach number limit and compared with Cassini Ion Neutral Mass Spectrometer (INMS) measurements to demonstrate that light gases (CH₄, H₂) preferentially escape over the heavy gas (N₂). The lightest gas (H₂) escapes with a flux 99% of its limiting flux, whereas CH₄ is restricted to ?75% of its limiting flux because there is insufficient solar power to support escape at the limiting rate. The respective calculated H₂ and CH₄ escape rates are 9.2×¹⁰²⁷ and 1.7×¹⁰²⁷ s−1, for a total of . From the calculated densities, mean free paths of N₂, CH₄, H₂, and macroscopic length scales, an extended region above the classic exobase is inferred where frequent collisions are still occurring and thermal heat conduction can deliver power to lift the escaping gas out of the gravitational potential well. In this region rapid acceleration of CH₄ outflow occurs. With the thermal structure of Titan's thermosphere inferred from INMS data by Müller-Wodarg et al. [Müller-Wodarg, I.C.F., Yelle, R.V., Cui, J., Waite Jr., J.H., 2008. J. Geophys. Res. 113, doi:10.1029/2007JE003033. E10005], in combination with calculated temperature profiles that include sputter induced plasma heating at the exobase, it is concluded that on average that the integrated, globally average, orbit-averaged, plasma heating rate during the Cassini epoch does not exceed ().     

The 2.9-4.2 micron spectrum of Saturn: Clouds and CH4, PH3, and NH3

We have used synthetic spectra to analyze a medium resolution 2.9-4.2 μm spectrum of Saturn's temperate region observed at UKIRT using CGS4. The synthetic spectra include CH₄, PH₃, and NH₃ lines, for which mixing ratios were adopted from recent Cassini results. The observed absorption features in the spectrum are well accounted for by lines of these molecular species formed 22 +/− 8 km above the 1 bar pressure level at ∼610 mbar. The influence of optically thin haze particles at higher altitudes on the spectrum is not pronounced, with higher spectral resolution probably required to constrain the effects of haze in this wavelength region. Fluorescent line emission by CH₄ in its ν₃ and ν₃+ν₄−ν₄ bands, detected in the 3.2-3.5 μm region, originates between 400 km (∼0.06 mbar) and 800 km (∼0.01 μbar) above the 1 bar level, with peak contributions from the two major contributing bands at 550 km (∼3 μbar) and 700 km (∼0.1 μbar), respectively.     

New Horizons Alice ultraviolet observations of a stellar occultation by Jupiter’s atmosphere

The Alice ultraviolet spectrograph onboard the New Horizons spacecraft observed two occultations of the bright star χ Ophiucus by Jupiter’s atmosphere on February 22 and 23, 2007 during the approach phase of the Jupiter flyby. The ingress occultation probed the atmosphere at 32°N latitude near the dawn terminator, while egress probed 18°N latitude near the dusk terminator. A detailed analysis of both the ingress and egress occultations, including the effects of molecular hydrogen, methane, acetylene, ethylene, and ethane absorptions in the far ultraviolet (FUV), constrains the eddy diffusion coefficient at the homopause level to be  cm² s⁻¹, consistent with Voyager measurements and other analyses (Festou, M.C., Atreya, S.K., Donahue, T.M., Sandel, B.R., Shemansky, D.E., Broadfoot, A.L. [1981]. J. Geophys. Res. 86, 5717-5725; Vervack Jr., R.J., Sandel, B.R., Gladstone, G.R., McConnell, J.C., Parkinson, C.D. [1995]. Icarus 114, 163-173; Yelle, R.V., Young, L.A., Vervack Jr., R.J., Young, R., Pfister, L., Sandel, B.R. [1996]. J. Geophys. Res. 101 (E1), 2149-2162). However, the actual derived pressure level of the methane homopause for both occultations differs from that derived by [Festou et al., 1981] and [Yelle et al., 1996] from the Voyager ultraviolet occultations, suggesting possible changes in the strength of atmospheric mixing with time. We find that at 32°N latitude, the methane concentration is  cm⁻³ at 70,397 km, the methane concentration is  cm⁻³ at 70,383 km, the acetylene concentration is  cm⁻³ at 70,364 km, and the ethane concentration is  cm⁻³ at 70,360 km. At 18°N latitude, the methane concentration is  cm⁻³ at 71,345 km, the methane concentration is  cm⁻³ at 71,332 km, the acetylene concentration is cm⁻³ at 71,318 km, and the ethane concentration is  cm⁻³ at 71,315 km. We also find that the H₂ occultation light curve is best reproduced if the atmosphere remains cold in the microbar region such that the base of the thermosphere is located at a lower pressure level than that determined by in situ instruments aboard the Galileo probe (Seiff, A., Kirk, D.B., Knight, T.C.D., Young, R.E., Mihalov, J.D., Young, L.A., Milos, F.S., Schubert, G., Blanchard, R.C., Atkinson, D. [1998]. J. Geophys. Res. 103 (E10), 22857-22889) - the Sieff et al. temperature profile leads to too much absorption from H₂ at high altitudes. However, this result is highly model dependent and non-unique. The observations and analysis help constrain photochemical models of Jupiter’s atmosphere.     

Impact crater formed on sintered snow surface simulating porous icy bodies

To improve the scaling parameter controlling the impact crater formation in the strength regime, we conducted impact experiments on sintered snow targets with the dynamic strength continuously changed from 20 to 200 kPa, and the largest crater size formed on small icy satellites was considered by using the revised scaling parameter. Ice and snow projectiles were impacted on a snow surface with 36% porosity at an impact velocity from 31 m s⁻¹ to 150 m s⁻¹. The snow target was sintered at the temperature from −5 °C to −18 °C, and the snow dynamic strength was changed with the sintering duration at each temperature. We found that the mass ejected from the crater normalized by the projectile mass, π_V, was related to the ratio of the dynamic strength to the impact pressure, , as follows: , where the impact pressure was indicated by P = ρ_tC_0tv_i/2 with the target density of ρ_t, when the impact velocity, v_i, was much smaller than the bulk sound velocity C_0t (typically 1.8 km s⁻¹ in our targets). The ratio of the largest crater diameter to the diameter of the target body, d_max/D, was estimated by calculating the crater diameter at the impact condition for catastrophic disruption and then compared to the observed d_max/D of jovian and saturnian small satellites, in order to discuss the formation condition of these large d_max/D in the strength regime.     

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.     

Titan's corona: The contribution of exothermic chemistry

The contribution of exothermic ion and neutral chemistry to Titan's corona is studied. The production rates for fast neutrals N₂, CH₄, H, H₂, ³CH₂, CH₃, C₂H₄, C₂H₅, C₂H₆, N(⁴S), NH, and HCN are determined using a coupled ion and neutral model of Titan's upper atmosphere. After production, the formation of the suprathermal particles is modeled using a two-stream simulation, as they travel simultaneously through a thermal mixture of N₂, CH₄, and H₂. The resulting suprathermal fluxes, hot density profiles, and energy distributions are compared to the N₂ and CH₄ INMS exospheric data presented in [De La Haye, V., Waite Jr., J.H., Johnson, R.E., Yelle, R.V., Cravens, T.E., Luhmann, J.G., Kasprzak, W.T., Gell, D.A., Magee, B., Leblanc, F., Michael, M., Jurac, S., Robertson, I.P., 2007. J. Geophys. Res., doi:10.1029/2006JA012222, in press], and are found insufficient for producing the suprathermal populations measured. Global losses of nitrogen atoms and carbon atoms in all forms due to exothermic chemistry are estimated to be and .     

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×¹⁰⁻⁴.     

Comet nucleus size distributions from HST and Keck telescopes

The Wide Field Camera (WFC) on the Hubble Space Telescope and the Low Resolution Imaging Spectrograph (LRIS) on the Keck II telescope have been used to image 21 distant dynamically new, long-period (LP) and short-period (SP) Jupiter-family (JF) comet nuclei (near aphelion), as part of a long-term program to search for physical differences between short-period comets and Oort cloud comets. WFC data were obtained on Comets C/1987 H1 (Shoemaker) and C/1984 K1 (Shoemaker) during Cycle 5 (1995 December) and on C/1988 B1 (Shoemaker), C/1987 F1 (Torres), and C/1983 O1 (?ernis) during Cycle 6 (1997 April, May, and June). The HST comets were at heliocentric distances 20.4 < r[AU] < 29.5. Each comet observation was allocated 7 orbits, for ≈3.6 hrs of integration. The most difficult part of the image reduction was the removal of cosmic rays. We present our scheme for cosmic ray removal. None of the HST comet nuclei was detected to the 3-σ level at m_R∼27. The inferred upper limits to the nucleus radii are . The SP comets range in radius between , with a median value of R_N∼1.61 km. The LP comets ranged in size between <4.0-56 km. Over a range of radii between 1-10 km, the nuclei can be fit with a cumulative distribution N(>R_N)∝R_N^−α with α=1.45±0.05, and for nuclei in the range 2-5 km, α=1.91±0.06. Statistical analysis and modeling shows that the slopes of the observed TNO and JF comet distributions are not compatible, suggesting that the intrinsic distribution of JF comet nuclei is a differential a^−3.5 power law truncated at small nucleus radii between 0.3 and 2.0 km.     

Mass and interior of Enceladus from Cassini data analysis

Gravity results are available from radio Doppler data acquired by the Deep Space Network during the encounter of the Cassini spacecraft with Enceladus in February 2005. We report the mass of Enceladus to be (1.0798±0.0016)×¹⁰²⁰ kg, which implies a density of . For a core made of hydrated silicates with a density of 2500 kg m⁻³ the core radius is ∼190 km and the quadrupole moment C₂₂∼1.4×¹⁰⁻³. If Enceladus is in hydrostatic equilibrium, the larger than previously anticipated density implies that the recently proposed secondary spin-orbit resonance cannot be present. Therefore, the source of endogenic activity of Enceladus remains unexplained.     

Comet Hale-Bopp in outburst: Imaging the dynamics of icy particles with HST/NICMOS

Comet Hale-Bopp was imaged at wavelengths from 1.87 to 2.22 μm by HST/NICMOS in post-perihelion observations starting on UT 1997 August 27.95. Diffraction-limited (∼02) images were obtained at high signal-to-noise (∼1500) to probe the composition and dynamics of the inner coma and also the size and activity of the nucleus. The velocities of several unusual morphological features over a 1.7 h period, indicate that a significant outburst occurred 7.4 h prior to these images while the comet was at a heliocentric distance of 2.49 AU. Similar features are also apparent after re-analysis of pre-perihelion ground-based images. The inner coma (radius ?2500 km) is dominated by an “arc” feature, which expanded and became more diffuse with time. This feature can be modeled as the bright central portion of a “jet of outburst” from a near-equatorial region of the nucleus. Less prominent, time-variable linear and circular morphologies are also apparent. The expansion rates of both the arc feature and the circular morphologies imply a common origin and also suggest a grain size distribution with two broad maxima. In addition, several static linear features extend to the edge of the field of view (21,100 km). Radial brightness profiles are highly asymmetric and only approach a ρ⁻¹ decline at distances ?15,000 km. Images in a narrow-band filter at 2.04 μm exhibit a ∼4% absorption feature relative to nearly simultaneous images at wavelengths of 2.22, 1.90, and 1.87 μm. This absorption is attributed to H₂O ice in the coma grains. The spatial distribution and expansion velocity of the absorption at 2.04 μm indicate that these grains are associated with the outburst. The constancy of the absorption feature indicates no appreciable sublimation over 1.7 h. The unresolved nucleus has a flux density consistent with a 40±10 km diameter assuming a 4% geometric albedo.     

An analysis of Pluto occultation light curves using an atmospheric radiative-conductive model

We use a radiative-conductive model to least-squares fit Pluto stellar occultation light curve data. This model predicts atmospheric temperature based on surface temperature, surface pressure, surface radius, and CH₄ and CO mixing ratios, from which model light curves are to be calculated. The model improves upon previous techniques for deriving Pluto’s atmospheric thermal structure from stellar occultation light curves by calculating temperature (as a function of height) caused by heating and cooling by species in Pluto’s atmosphere, instead of a general assumption that temperature follows a power law with height or some other idealized function. We are able to fit for model surface radius, surface pressure, and CH₄ mixing ratio with one of the 2006 datasets and for surface pressure and CH₄ mixing ratio for other datasets from the years 1988, 2002, 2006, and 2008. It was not possible to fit for CO mixing ratio and surface temperature because the light curves are not sensitive to these parameters. We determine that the model surface radius, under the assumption of a stratosphere only (i.e. no troposphere) model in radiative-conductive balance, is . The CH₄ mixing ratio results are more scattered with time and are in the range of 1.8-9.4 × 10⁻³. The surface pressure results show an increasing trend from 1988 to 2002, although it is not as dramatic as the factor of 2 from previous studies.     

Analysis of Titan's neutral upper atmosphere from Cassini Ion Neutral Mass Spectrometer measurements

In this paper we present an in-depth study of the distributions of various neutral species in Titan's upper atmosphere, between 950 and 1500 km for abundant species (N₂, CH₄, H₂) and between 950 and 1200 km for other minor species. Our analysis is based on a large sample of Cassini/INMS (Ion Neutral Mass Spectrometer) measurements in the CSN (Closed Source Neutral) mode, obtained during 15 close flybys of Titan. To untangle the overlapping cracking patterns, we adopt Singular Value Decomposition (SVD) to determine simultaneously the densities of different species. Except for N₂, CH₄, H₂ and ⁴⁰Ar (as well as their isotopes), all species present density enhancements measured during the outbound legs. This can be interpreted as a result of wall effects, which could be either adsorption/desorption of these molecules or heterogeneous surface chemistry of the associated radicals on the chamber walls. In this paper, we provide both direct inbound measurements assuming ram pressure enhancement only and abundances corrected for wall adsorption/desorption based on a simple model to reproduce the observed time behavior. Among all minor species of photochemical interest, we have firm detections of C₂H₂, C₂H₄, C₂H₆, CH₃C₂H, C₄H₂, C₆H₆, CH₃CN, HC₃N, C₂N₂ and NH₃ in Titan's upper atmosphere. Upper limits are given for other minor species.The globally averaged distributions of N₂, CH₄ and H₂ are each modeled with the diffusion approximation. The N₂ profile suggests an average thermospheric temperature of 151 K. The CH₄ and H₂ profiles constrain their fluxes to be and , referred to Titan's surface. Both fluxes are significantly higher than the Jeans escape values. The INMS data also suggest horizontal/diurnal variations of temperature and neutral gas distribution in Titan's thermosphere. The equatorial region, the ramside, as well as the nightside hemisphere of Titan appear to be warmer and present some evidence for the depletion of light species such as CH₄. Meridional variations of some heavy species are also observed, with a trend of depletion toward the north pole. Though some of the above variations might be interpreted by either the solar-driven models or auroral-driven models, a physical scenario that reconciles all the observed horizontal/diurnal variations in a consistent way is still missing. With a careful evaluation of the effect of restricted sampling, some of the features shown in the INMS data are more likely to be observational biases.