Waves in Cassini UVIS stellar occultations: 2. The C ring

We performed a complete wavelet analysis of Saturn’s C ring on 62 stellar occultation profiles. These profiles were obtained by Cassini’s Ultraviolet Imaging Spectrograph High Speed Photometer. We used a WWZ wavelet power transform to analyze them. With a co-adding process, we found evidence of 40 wavelike structures, 18 of which are reported here for the first time. Seventeen of these appear to be propagating waves (wavelength changing systematically with distance from Saturn). The longest new wavetrain in the C ring is a 52-km-long wave in a plateau at 86,397 km. We produced a complete map of resonances with external satellites and possible structures rotating with Saturn’s rotation period up to the eighth order, allowing us to associate a previously observed wave with the Atlas 2:1 inner Lindblad resonance (ILR) and newly detected waves with the Mimas 6:2 ILR and the Pandora 4:2 ILR. We derived surface mass densities and mass extinction coefficients, finding σ = 0.22(±0.03) g cm⁻² for the Atlas 2:1 ILR, σ = 1.31(±0.20) g cm⁻² for the Mimas 6:2 ILR, and σ = 1.42(±0.21) g cm⁻² for the Pandora 4:2 ILR. We determined a range of mass extinction coefficients (κ = τ/σ) for the waves associated with resonances with κ = 0.13 (±0.03) to 0.28(±0.06) cm² g⁻¹, where τ is the optical depth. These values are higher than the reported values for the A ring (0.01-0.02 cm² g⁻¹) and the Cassini Division (0.07-0.12 cm² g⁻¹ from Colwell et al. (Colwell, J.E., Cooney, J.H., Esposito, L.W., Srem?evi?, M. [2009]. Icarus 200, 574-580)). We also note that the mass extinction coefficient is probably not constant across the C ring (in contrast to the A ring and the Cassini Division): it is systematically higher in the plateaus than elsewhere, suggesting smaller particles in the plateaus. We present the results of our analysis of these waves in the C ring and estimate the mass of the C ring to be between3.7(±0.9) × 10¹⁶ kg and 7.9(±2.0) × 10¹⁶ kg (equivalent to an icy satellite of radius between 28.0(±2.3) km and 36.2(±3.0) km with a density of 400 kg m⁻³, close to that of Pan or Atlas). Using the ring viscosity derived from the wave damping length, we also estimate the vertical thickness of the C ring between 1.9(±0.4) m and 5.6(±1.4) m, comparable to the vertical thickness of the Cassini Division.     

The effect of Nix and Hydra on the putative Pluto–Charon dust cloud

Impact-generated dust clouds around airless bodies have been observed or suggested to be present throughout the solar system, including around the Martian, Galilean and Saturnian satellites. Simulations have assessed Pluto and Charon as sources of a possible dust cloud or torus and found that such a cloud would be dominated by Charon-produced ejecta and would have an optical depth of τ≈10⁻¹¹. These simulations were conducted before the discovery of two additional, small satellites of Pluto, Nix and Hydra. These small moons may yield impact-generated dust in excess of their larger counterparts due to their lower escape velocities, despite their smaller cross sections. In this paper, we extend a previous model of the Pluto–Charon dust cloud to include Nix and Hydra, both as sinks for Pluto- and Charon-generated dust and as sources of impact-generated dust. We find that Nix- and Hydra-generated dust grains outlive Pluto and Charon dust grains significantly and are the dominant contributors of dust in the Pluto–Charon system. Furthermore, we estimate the net geometric optical depth of grains between 0.1 and to be on the order of 10⁻⁷.     

A dust cloud around Pluto and Charon

We suggest that Pluto and Charon are immersed in a tenuous dust cloud. The cloud consists of ejecta from Pluto and—especially—Charon, released from their surfaces by impacts of micrometeoroids originating from Edgeworth-Kuiper belt objects. The motion of the ejected grains is dominated by the gravity of Pluto and Charon, which determines a pear-shape of the densest part of the cloud. While the production rates of escaping particles from both sides are comparable, the lifetimes of the Charon particles inside the Hill sphere of Pluto-Charon with respect to the Sun are much longer than of the Pluto ejecta, so that the cloud is composed predominantly of Charon grains. The dust cloud is dense enough to be detected with an in situ dust detector onboard a future space mission to Pluto. The cloud's maximum optical depth of τ≈3×10⁻¹¹ is, however, too low to allow remote sensing observations.     

Eclipse reappearances of Io: Time-resolved spectroscopy (1.9-4.2 μm)

We obtained time-resolved, near-infrared spectra of Io during the 60-90 min following its reappearance from eclipse by Jupiter on five occasions in 2004. The purpose was to search for spectral changes, particularly in the well-known SO₂ frost absorption bands, that would indicate surface-atmosphere exchange of gaseous SO₂ induced by temperature changes during eclipse. These observations were a follow-on to eclipse spectroscopy observations in which Bellucci et al. [Bellucci et al., 2004. Icarus 172, 141-148] reported significant changes in the strengths of two strong SO₂ bands in data acquired with the VIMS instrument aboard the Cassini spacecraft. One of the bands (4.07 μm [ν₁ + ν₃]) observed by Bellucci et al. is visible from ground-based observatories and is included in our data. We detected no changes in Io’s spectrum at any of the five observed events during the approximately 60-90 min during which spectra were obtained following Io’s emergence from Jupiter’s shadow. The areas of the three strongest SO₂ bands in the region 3.5-4.15 μm were measured for each spectrum; the variation of the band areas with time does not exceed that which can be explained by the Io’s few degrees of axial rotation during the intervals of observation, and in no case does the change in band strength approach that seen in the Cassini VIMS data. Our data are of sufficient quality and resolution to show the weak 2.198 μm (4549.6 cm⁻¹) 4ν₁ band of SO₂ frost on Io for what we believe is the first time. At one of the events (June 22, 2004), we began the acquisition of spectra ∼6 min before Io reappeared from Jupiter’s shadow, during which time it was detected through its own thermal emission. No SO₂ bands were superimposed on the purely thermal spectrum on this occasion, suggesting that the upper limit to condensed SO₂ in the vertical column above Io’s surface was ∼4 × 10⁻⁵ g cm⁻².     

A search for ethane on Pluto and Triton

We present here a search for solid ethane, C₂H₆, on the surfaces of Pluto and Triton, based on near-infrared spectral observations in the H and K bands (1.4-2.45 μm) using the Very Large Telescope (VLT) and the United Kingdom Infrared Telescope (UKIRT). We model each surface using a radiative transfer model based on Hapke theory (Hapke, B. [1993]. Theory of Reflectance and Emittance Spectroscopy. Cambridge University Press, Cambridge, UK) with three basic models: without ethane, with pure ethane, and with ethane diluted in nitrogen. On Pluto we detect weak features near 2.27, 2.405, 2.457, and 2.461 μm that match the strongest features of pure ethane. An additional feature seen at 2.317 μm is shifted to longer wavelengths than ethane by at least 0.002 μm. The strength of the features seen in the models suggests that pure ethane is limited to no more than a few percent of the surface of Pluto. On Triton, features in the H band could potentially be explained by ethane diluted in N₂, however, the lack of corresponding features in the K band makes this unlikely (also noted by Quirico et al. (Quirico, E., Doute, S., Schmitt, B., de Bergh, C., Cruikshank, D.P., Owen, T.C., Geballe, T.R., Roush, T.L. [1999]. Icarus 139, 159-178)). While Cruikshank et al. (Cruikshank, D.P., Mason, R.E., Dalle Ore, C.M., Bernstein, M.P., Quirico, E., Mastrapa, R.M., Emery, J.P., Owen, T.C. [2006]. Bull. Am. Astron. Soc. 38, 518) find that the 2.406-μm feature on Triton could not be completely due to ¹³CO, our models show that it could not be accounted for entirely by ethane either. The multiple origin of this feature complicates constraints on the contribution of ethane for both bodies.     

New insights on the binary Asteroid 121 Hermione

We report on the results of a 6-month photometric study of the main-belt binary C-type Asteroid 121 Hermione, performed during its 2007 opposition. We took advantage of the rare observational opportunity afforded by one of the annual equinoxes of Hermione occurring close to its opposition in June 2007. The equinox provides an edge-on aspect for an Earth-based observer, which is well suited to a thorough study of Hermione’s physical characteristics. The catalog of observations carried out with small telescopes is presented in this work, together with new adaptive optics (AO) imaging obtained between 2005 and 2008 with the Yepun 8-m VLT telescope and the 10-m Keck telescope. The most striking result is confirmation that Hermione is a bifurcated and elongated body, as suggested by Marchis, et al. [Marchis, F., Hestroffer, D., Descamps, P., Berthier, J., Laver, C., de Pater, I., 2005. Icarus 178, 450-464]. A new effective diameter of 187 ± 6 km was calculated from the combination of AO, photometric and thermal observations. The new diameter is some 10% smaller than the hitherto accepted radiometric diameter based on IRAS data. The reason for the discrepancy is that IRAS viewed the system almost pole-on. New thermal observations with the Spitzer Space Telescope agree with the diameter derived from AO and lightcurve observations. On the basis of the new AO astrometric observations of the small 32-km diameter satellite we have refined the orbit solution and derived a new value of the bulk density of Hermione of 1.4 + 0.5/−0.2 g cm⁻³. We infer a macroscopic porosity of ∼33 + 5/−20%.     

Meridional transport of HCN from SL9 impacts on Jupiter

In July 1994, the Shoemaker-Levy 9 (SL9) impacts introduced hydrogen cyanide (HCN) to Jupiter at a well confined latitude band around −44°, over a range of specific longitudes corresponding to each of the 21 fragments (Bézard et al. 1997, Icarus 125, 94-120). This newcomer to Jupiter's stratosphere traces jovian dynamics. HCN rapidly mixed with longitude, so that observations recorded later than several months after impact witnessed primarily the meridional transport of HCN north and south of the impact latitude band. We report spatially resolved spectroscopy of HCN emission 10 months and 6 years following the impacts. We detect a total mass of HCN in Jupiter's stratosphere of 1.5±0.7×10¹³ g in 1995 and 1.7±0.4×10¹³ g in 2000, comparable to that observed several days following the impacts (Bézard et al. 1997, Icarus 125, 94-120). In 1995, 10 months after impact, HCN spread to −30° and −65° latitude (half column masses), consistent with a horizontal eddy diffusion coefficient of K_yy=2-3×10¹⁰ cm² s⁻¹. Six years following impact HCN is observed in the northern hemisphere, while still being concentrated at 44° south latitude. Our meridional distribution of HCN suggests that mixing occurred rapidly north of the equator, with K_yy=2-5×10¹¹ cm² s⁻¹, consistent with the findings of Moreno et al. (2003, Planet. Space Sci. 51, 591-611) and Lellouch et al. (2002, Icarus 159, 112-131). These inferred eddy diffusion coefficients for Jupiter's stratosphere at 0.1-0.5 mbar generally exceed those that characterize transport on Earth. The low abundance of HCN detected at high latitudes suggests that, like on Earth, polar regions are dynamically isolated from lower latitudes.     

Thermal model of water and CO activity of Comet C/1995 O1 (Hale-Bopp)

An investigation of the activity of Comet C/1995 O1 (Hale-Bopp) with a thermophysical nucleus model that does not rely on the existence of amorphous ice is presented. Our approach incorporates recent observations allowing to constrain important parameters that control cometary activity. The model accounts for heat conduction, heat advection, gas diffusion, sublimation, and condensation in a porous ice-dust matrix with moving boundaries. Erosion due to surface sublimation of water ice leads to a moving boundary. The movement of the boundary is modeled by applying a temperature remapping technique which allows us to account for the loss in the internal energy of the eroded surface material. These kind of problems are commonly referred to as Stefan problems. The model takes into account the diurnal rotation of the nucleus and seasonal effects due to the strong obliquity of Hale-Bopp as reported by Jorda et al. (Jorda, L., Rembor, K., Lecacheux, J., Colom, P., Colas, F., Frappa, E., Lara, L.M. [1997]. Earth Moon Planets 77, 167-180). Only bulk sublimation of water and CO ice are considered without further assumptions such as amorphous ices with certain amount of occluded CO gas. Confined and localized activity patterns are investigated following the reports of Lederer and Campins (Lederer, S.M., Campins, H. [2002]. Earth Moon Planets 90, 381-389) about the chemical heterogeneity of Hale-Bopp and of Bockelée-Morvan et al. (Bockelée-Morvan, D., Henry, F., Biver, N., Boissier, J., Colom, P., Crovisier, J., Despois, D., Moreno, R., Wink, J. [2009]. Astron. Astrophys. 505, 825-843) about a strong CO source at a latitude of 20°. The best fit to the observations of Biver et al. (Biver, N. et al. [2002]. Earth Moon Planets 90, 5-14) is obtained with a low thermal conductivity of 0.01 W m⁻¹ K⁻¹. This is in agreement with recent results of the Deep Impact mission to 9P/Tempel 1 (Groussin, O., A’Hearn, M.F., Li, J.-Y., Thomas, P.C., Sunshine, J.M., Lisse, C.M., Meech, K.J., Farnham, T.L., Feaga, L.M., Delamere, W.A. [2007]. Icarus 187, 16-25) and with previous thermal simulations (Kührt, E. [1999]. Space Sci. Rev. 90, 75-82). The water production curve matches the production rates well from −4 AU pre-perihelion to the outgoing leg while the model does not reproduce so well the water production beyond 4 AU pre-perihelion. The CO production curve is a good fit to the measurements of Biver et al. (2002) over the whole measured heliocentric range from −7 AU pre- to 15 AU post-perihelion.     

In situ observation of penetration process in silica aerogel: Deceleration mechanism of hard spherical projectiles

A large number of cometary dust particles were captured with low-density silica aerogels by NASA’s Stardust Mission. Knowledge of the details of the capture mechanism of hypervelocity particles in silica aerogel is needed in order to correctly derive the original particle features from impact tracks. However, the mechanism has not been fully understood yet. We shot hard spherical projectiles of several different materials into silica aerogel of density 60 mg cm⁻³ and observed their penetration processes using an image converter or a high-speed video camera. In order to observe the deceleration of projectiles clearly, we carried out impact experiments at two velocity ranges; ∼4 km s⁻¹ and ∼200 m s⁻¹. From the movies we took, it was indicated that the projectiles were decelerated by hydrodynamic force which was proportional to v² (v: projectile velocity) during the faster penetration process (∼4 km s⁻¹) and they were merely overcoming the aerogel crushing strength during the slower penetration process (∼200 m s⁻¹). We applied these deceleration mechanisms for whole capture process to calculate the track length. Our model well explains the track length in the experimental data set by Burchell et al. (Burchell, M.J., Creighton, J.A., Cole, M.J., Mann, J., Kearsley, A.T. [2001]. Meteorit. Planet. Sci. 36, 209-221).     

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.     

Formation of Jupiter using opacities based on detailed grain physics

Numerical simulations, based on the core-nucleated accretion model, are presented for the formation of Jupiter at 5.2 AU in three primordial disks with three different assumed values of the surface density of solid particles. The grain opacities in the envelope of the protoplanet are computed using a detailed model that includes settling and coagulation of grains and that incorporates a recalculation of the grain size distribution at each point in time and space. We generally find lower opacities than the 2% of interstellar values used in previous calculations (Hubickyj, O., Bodenheimer, P., Lissauer, J.J. [2005]. Icarus 179, 415-431; Lissauer, J.J., Hubickyj, O., D’Angelo, G., Bodenheimer, P. [2009]. Icarus 199, 338-350). These lower opacities result in more rapid heat loss from and more rapid contraction of the protoplanetary envelope. For a given surface density of solids, the new calculations result in a substantial speedup in formation time as compared with those previous calculations. Formation times are calculated to be 1.0, 1.9, and 4.0 Myr, and solid core masses are found to be 16.8, 8.9, and 4.7 M_⊕, for solid surface densities, σ, of 10, 6, and 4 g cm⁻², respectively. For σ = 10 and σ = 6 g cm⁻², respectively, these formation times are reduced by more than 50% and more than 80% compared with those in a previously published calculation with the old approximation to the opacity.     

Influence of high abundances of aerosols on the electrical conductivity of the Titan atmosphere

Observations of optical depth and scattering by instrumentation onboard the Huygens probe have been used by Tomasko et al. [Tomasko et al., 2005. Rain, winds and haze during Huygens probe's descent to Titan's surface. Nature 438 (8), 765-778] to deduce that the size and abundance of Titan aerosols could be nearly independent of altitude. Here we show that by assuming a constant mass flux with altitude and using the measured optical depth as a constraint, we obtain more realistic size and abundance distributions. In particular, the calculated abundance decreases from 3.5×10⁷ m⁻³ at 100 km to 8×10⁶ m⁻³ near the surface while the particle radius varies from 0.25 μm at 150 km to 1.1 μm at the surface. These distributions are consistent with the reported measurements for these quantities. Our results are then employed to compute electron and ion densities and conductivities for various solar UV photoelectron emission thresholds. Our model shows that to get agreement with the published (preliminary) conductivity measurements, photoemission cannot be an important source of electrons and ions. To get agreement with the electron and ion conductivity observations, both an additional population of aerosol embryos above 50 km and a trace amount of an electrophillic molecular species below 50 km are needed.     

Transport of ionizing radiation in terrestrial-like exoplanet atmospheres

The propagation of ionizing radiation through model atmospheres of terrestrial-like exoplanets is studied for a large range of column densities and incident photon energies using a Monte Carlo code we have developed to treat Compton scattering and photoabsorption. Incident spectra from parent star flares, supernovae, and gamma-ray bursts are modeled and compared to energetic particles in importance. Large irradiation events with fluences of 10⁶-10⁹ erg cm⁻² at the conventional habitable zone can occur at a rate from many per day (flares from young low-mass parent stars) to ∼100 per Gyr (supernovae and gamma-ray bursts). We find that terrestrial-like exoplanets with atmospheres thinner than about 100 g cm⁻² block nearly all X-rays, but transmit and reprocess a significant fraction of incident γ-rays, producing a characteristic, flat surficial spectrum. Thick atmospheres (?100 g cm⁻²) efficiently block even γ-rays, but nearly all the incident energy is redistributed into diffuse UV and visible aurora-like emission, increasing the effective atmospheric transmission by many orders of magnitude. Depending on the presence of molecular UV absorbers and atmospheric thickness, up to 10% of the incident energy can reach the surface as UV reemission. For the Earth, between 2×10⁻³ and 4×10⁻² of the incident flux reaches the ground in the biologically effective 200-320 nm range, depending on O₂/O₃ shielding. For atmospheres thicker than ∼50 g cm⁻² in the case of pure Rayleigh scattering and ∼100 g cm⁻² in the case of O₂/O₃ absorption, the UV reemission exceeds the surficial transmitted ionizing radiation. We also discuss the effects of angle of incidence and derive a modified two-stream approximation solution for the UV transfer. Finally, we suggest that transient atmospheric ionization layers can be frequently created at altitudes lower than the equilibrium layers that result from steady irradiation and winds from the parent star. We suggest that these events can produce frequent fluctuations in atmospheric ionization levels and surficial UV fluxes on terrestrial-like planets.     

Vertical profile of H2SO4 vapor at 70-110 km on Venus and some related problems

The vertical profile of H₂SO₄ vapor is calculated using current atmospheric and thermodynamic data. The atmospheric data include the H₂O profiles observed at 70-112 km by the SOIR solar occultations, the SPICAV-UV profiles of the haze extinction at 220 nm, the VeRa temperature profiles, and a typical profile of eddy diffusion. The thermodynamic data are the saturated vapor pressures of H₂O and H₂SO₄ and chemical potentials of these species in sulfuric acid solutions. The calculated concentration of sulfuric acid in the cloud droplets varies from 85% at 70 km to a minimum of 70% at 90 km and then gradually increasing to 90-100% at 110 km. The H₂SO₄ vapor mixing ratio is ∼10⁻¹² at 70 and 110 km with a deep minimum of 3 × 10⁻¹⁸ at 88 km. The H₂O-H₂SO₄ system matches the local thermodynamic equilibrium conditions up to 87 km. The column photolysis rate of H₂SO₄ is 1.6 × 10⁵ cm⁻² s⁻¹ at 70 km and 23 cm⁻² s⁻¹ at 90 km. The calculated abundance of H₂SO₄ vapor at 90-110 km and its photolysis rate are smaller than those presented in the recent model by Zhang et al. (Zhang, X., Liang, M.C., Montmessin, F., Bertaux, J.L., Parkinson, C., Yung, Y.L. [2010]. Nat. Geosci. 3, 834-837) by factors of 10⁶ and 10⁹, respectively. Assumptions of 100% sulfuric acid, local thermodynamic equilibrium, too warm atmosphere, supersaturation of H₂SO₄ (impossible for a source of SO_X), and cross sections for H₂SO₄·H₂O (impossible above the pure H₂SO₄) are the main reasons of this huge difference. Significant differences and contradictions between the SPICAV-UV, SOIR, and ground-based submillimeter observations of SO_X at 70-110 km are briefly discussed and some weaknesses are outlined. The possible source of high altitude SO_X on Venus remains unclear and probably does not exist.     

NN detectability in Pluto’s atmosphere

Based on the vapor pressure behavior of Pluto’s surface ices, Pluto’s atmosphere is expected to be predominantly composed of N₂ gas. Measurement of the N₂ isotopologue ¹⁵N/¹⁴N ratio within Pluto’s atmosphere would provide important clues to the evolution of Pluto’s atmosphere from the time of formation to its present state. The most straightforward way of determining the N₂ isotopologue ¹⁵N/¹⁴N ratio in Pluto’s atmosphere is via spectroscopic observation of the ¹⁴N¹⁵N gas species. Recent calculations of the 80–100 nm absorption behavior of the ¹⁴N₂ and ¹⁴N¹⁵N isotopologues by Heays et al. (Heays, A.N. et al. [2011]. J. Chem. Phys. 135, 244301), Lewis et al. (Lewis, B.R., Heays, A.N., Gibson, S.T., Lefebvre-Brion, H., Lefebvre, R. [2008]. J. Chem. Phys. 129, 164306); Lewis et al. (Lewis, B.R., Gibson, S.T., Zhang, W., Lefebvre-Brion, H., Robbe, J.-M. [2005]. J. Chem. Phys. 122, 144302), and Haverd et al. (Haverd, V.E., Lewis, B.R., Gibson, S.T., Stark, G. [2005]. J. Chem. Phys. 123, 214304) show that the peak magnitudes of the ¹⁴N₂ and ¹⁴N¹⁵N absorption bandhead cross-sections are similar, but the locations of the bandhead peaks are offset in wavelength by ∼0.05–0.1 nm. These offsets make the segregation of the ¹⁴N₂ and ¹⁴N¹⁵N absorption signatures possible. We use the most recent N₂ isotopologue absorption cross-section calculations and the atmospheric density profiles resulting from photochemical models developed by Krasnopolsky and Cruickshank (Krasnopolsky, V.A., Cruickshank, D.P. [1999]. J. Geophys. Res. 104, 21979–21996) to predict the level of solar light that will be transmitted through Pluto’s atmosphere as a function of altitude during a Pluto solar occultation. We characterize the detectability of the isotopic absorption signature per altitude assuming ¹⁴N¹⁵N concentrations ranging from 0.1% to 2% of the ¹⁴N₂ density and instrumental spectral resolutions ranging from 0.01 to 0.3 nm. Our simulations indicate that optical depth of unity is attained in the key ¹⁴N¹⁵N absorption bands located between 85 and 90 nm at altitudes ∼1100–1600 km above Pluto’s surface. Additionally, an ¹⁴N¹⁵N isotope absorption depth ∼4–15% is predicted for observations obtained at these altitudes at a spectral resolution of ∼0.2–0.3 nm, if the N₂ isotopologue ¹⁵N/¹⁴N percent ratio is comparable to the 0.37–0.6% ratio observed at Earth, Titan and Mars. If we presume that the predicted absorption depth must be at least 25% greater than the expected observational uncertainty, then it follows that a statistically significant detection of these signatures and constraint of the N₂ isotopologue ¹⁴N/¹⁵N ratio within Pluto’s atmosphere will be possible if the attainable observational signal-to noise (S/N) ratio is ?9. The New Horizons (NH) Mission will be able to obtain high S/N, 0.27–0.35 nm full-width half-max 80–100 nm spectral observations of Pluto using the Alice spectrograph. Based on the NH/Alice specifications we have simulated 0.3 nm spectral resolution solar occultation spectra for the 1100–1600 km altitude range, assuming 30 s integration times. These simulations indicate that NH/Alice will obtain spectral observations within this altitude range with a S/N ratio ∼25–50, and should be able to reliably detect the ¹⁴N¹⁵N gas absorption signature between 85 and 90 nm if the ¹⁴N¹⁵N concentration is ∼0.3% or greater. This, additionally, implies that the non-detection of the ¹⁴N¹⁵N species in the 1100–1600 km range by NH/Alice may be used to reliably establish an upper limit to the N₂ isotopologue ¹⁵N/¹⁴N ratio within Pluto’s atmosphere. Similar results may be derived from 0.2 to 0.3 nm spectral resolution observations of any other N₂-rich Solar System or exoplanet atmosphere, provided the observations are attained with similar S/N levels.     

Numerical modelling of heating in porous planetesimal collisions

Collisions between planetesimals at speeds of several kilometres per second were common during the early evolution of our Solar System. However, the collateral effects of these collisions are not well understood. In this paper, we quantify the efficiency of heating during high-velocity collisions between planetesimals using hydrocode modelling. We conducted a series of simulations to test the effect on shock heating of the initial porosity and temperature of the planetesimals, the relative velocity of the collision and the relative size of the two colliding bodies. Our results show that while heating is minor in collisions between non-porous planetesimals at impact velocities below 10 km s⁻¹, in agreement with previous work, much higher temperatures are reached in collisions between porous planetesimals. For example, collisions between nearly equal-sized, porous planetesimals can melt all, or nearly all, of the mass of the bodies at collision velocities below 7 km s⁻¹. For collisions of small bodies into larger ones, such as those with an impactor-to-target mass ratio below 0.1, significant localised heating occurs in the target body. At impact velocities as low as 5 km s⁻¹, the mass of melt will be nearly double the mass of the impactor, and the mass of material shock heated by 100 K will be nearly 10 times the mass of the impactor. We present a first-order estimate of the cumulative effects of impact heating on a porous planetesimal parent body by simulating the impact of a population of small bodies until a disruptive event occurs. Before disruption, impact heating is volumetrically minor and highly localised; in no case was more than about 3% of the parent body heated by more than 100 K. However, heating during the final disruptive collision can be significant; in about 10% of cases, almost all of the parent body is heated to 700 K (from an initial temperature of ∼300 K) and more than a tenth of the parent body mass is melted. Hence, energetic collisions between planetesimals could have had important effects on the thermal evolution of primitive materials in the early Solar System.     

Morphology, temperature, and eruption dynamics at Pele

The Pele region of Io has been the site of vigorous volcanic activity from the time of the first Voyager I observations in 1979 up through the final Galileo ones in 2001. There is high-temperature thermal emission from a visibly dark area that is thought to be a rapidly overturning lava lake, and is also the source of a large sulfur-rich plume. We present a new analysis of Voyager I visible wavelength images, and Galileo Solid State Imager (SSI) and Near Infrared Mapping Spectrometer (NIMS) thermal emission observations which better define the morphology of the region and the intensity of the emission. The observations show remarkable correlations between the locations of the emission and the features seen in the Voyager images, which provide insight into eruption mechanisms and constrain the longevity of the activity. We also analyze an additional wavelength channel of NIMS data (1.87 μm) which paradoxically, because of reduced sensitivity, allows us to estimate temperatures at the peak locations of emission. Measurements of eruption temperatures on Io are crucial because they provide our best clues to the composition of the magma. High color temperatures indicative of ultramafic composition have been reported for the Pillan hot spot and possibly for Pele, although recent work has called into question the requirement for magma temperatures above those expected for ordinary basalts. Our new analysis of the Pele emission near the peak of the hot spot shows color temperatures near the upper end of the basalt range during the I27 and I32 encounters. In order to analyze the observed color temperatures we also present an analytical model for the thermal emission from fire-fountains, which should prove generally useful for analyzing similar data. This is a modification of the lava flow emission model presented in Howell (Howell, R.R. [1997]. Icarus 127, 394-407), adapted to the fire-fountain cooling curves first discussed in Keszthelyi et al. (Keszthelyi, L., Jaeger, W., Milazzo, M., Radebaugh, J., Davies, A.G., Mitchell, K.L. [2007]. Icarus 192, 491-502). When applied to the I32 observations we obtain a fire-fountain mass eruption rate of 5.1 × 10⁵ kg s⁻¹ for the main vent area and 1.4 × 10⁴ kg s⁻¹ for each of two smaller vent regions to the west. These fire-fountain rates suggest a solution to the puzzling lack of extensive lava flows in the Pele region. Much of the erupted lava may be ejected at high speed into the fire-fountains and plumes, creating dispersed pyroclastic deposits rather than flows. We compare gas and silicate mass eruption rates and discuss briefly the dynamics of this ejection model and the observational evidence.     

Al decay: Heat production and a revised age for Iapetus

We revisit the appropriate energies to be used for computing heat production from ²⁶Al decay. Due to the complexity of the decay scheme of this radioisotope, previous geophysical studies have used values ranging from 1.2 to 4 MeV per decay. The upper bound corresponds to the difference in mass energy between the ²⁶Al and ²⁶Mg ground states. This includes energy carried away by neutrinos, which does not contribute to heating planetary material. The lower bound does not account for the heating caused by the absorption of the γ rays from the excited ²⁶Mg, or for the annihilation energy deposited in the material if the decay occurs inside even small planetesimals. Based on the calculations described by Schramm et al. [Schramm, D., Tera, F., Wasserburg, G.J., 1970. The isotopic abundance of ²⁶Mg and limits on ²⁶Al in the early Solar System. Earth Planet. Sci. Lett. 10, 44-59] updated with the most recent nuclear constants, we recommend using a heat production value of 3.12 MeV per decay, which is the total energy of disintegration minus the energy carried off by the neutrinos. This heat production value is higher than the value used in the modeling of Iapetus by Castillo-Rogez et al. [Castillo-Rogez, J., Matson, D.L., Sotin, C., Johnson, T.V., Lunine, J.I., Thomas, P.C., 2007. Iapetus’ geophysics: Rotation rate, shape, and equatorial ridge. Icarus 190, 179-202] by about a factor 2.5. The resulting estimate of the time of formation of Iapetus is shifted by about 1 Myr, to between ∼3.4 and 5.4 Myr after the production of the calcium-aluminum inclusions (CAIs).     

A Martian PFS average spectrum: Comparison with ISO SWS

The evaluation of the planetary Fourier spectrometer performance at Mars is presented by comparing an average spectrum with the ISO spectrum published by Lellouch et al. [2000. Planet. Space Sci. 48, 1393.]. First, the average conditions of Mars atmosphere are compared, then the mixing ratios of the major gases are evaluated. Major and minor bands of CO₂ are compared, from the point of view of features characteristics and bands depth. The spectral resolution is also compared using several solar lines. The result indicates that PFS radiance is valid to better than 1% in the wavenumber range 1800-4200 cm⁻¹ for the average spectrum considered (1680 measurements). The PFS monochromatic transfer function generates an overshooting on the left-hand side of strong narrow lines (solar or atmospheric). The spectral resolution of PFS is of the order of 1.3 cm⁻¹ or better. A large number of narrow features to be identified are discovered.