Impact-induced N2 production from ammonium sulfate: Implications for the origin and evolution of N2 in Titan’s atmosphere

Chemical reactions and volatile supply through hypervelocity impacts may have played a key role for the origin and evolution of both planetary and satellite atmospheres. In this study, we evaluate the role of impact-induced N₂ production from reduced nitrogen-bearing solids proposed to be contained in Titan’s crust, ammonium sulfate ((NH₄)₂SO₄), for the replenishment of N₂ to the atmosphere in Titan’s history. To investigate the conversion of (NH₄)₂SO₄ into N₂ by hypervelocity impacts, we measured gases released from (NH₄)₂SO₄ that was exposed to hypervelocity impacts created by a laser gun. The sensitivity and accuracy of the measurements were enhanced by using an isotope labeling technique for the target. We obtained the efficiency of N₂ production from (NH₄)₂SO₄ as a function of peak shock pressure ranging from ∼8 to ∼45 GPa. Our results indicate that the initial and complete shock pressures for N₂ degassing from (NH₄)₂SO₄ are ∼10 and ∼25 GPa, respectively. These results suggest that cometary impacts on Titan (i.e., impact velocity v_i > ∼8 km/s) produce N₂ efficiently; whereas satellitesimal impacts during the accretion (i.e., v_i < 4 km/s) produce N₂ only inefficiently. Even when using the proposed small amount of (NH₄)₂SO₄ content in the crust (∼4 wt.%) (Fortes, A.D. et al., 2007. Icarus 188, 139-153), the total amount of N₂ provided through cometary impacts over 4.5 Ga reaches ∼2-6 times the present atmospheric N₂ (i.e., ∼7 × 10²⁰-2 × 10²¹ [mol]) based on the measured production efficiency and results of a hydrodynamic simulation of cometary impacts onto Titan. This implies that cometary impacts onto Titan’s crust have the potential to account for a large part of the present N₂ through the atmospheric replenishment after the accretion.     

Molecular hydrogen in Titan’s atmosphere: Implications of the measured tropospheric and thermospheric mole fractions

The third most abundant species in Titan’s atmosphere is molecular hydrogen with a tropospheric/lower stratospheric mole fraction of 0.001 derived from Voyager and Cassini infrared measurements. The globally averaged thermospheric H₂ mole fraction profile from the Cassini Ion Neutral Mass Spectrometer (INMS) measurements implies a small positive gradient in the H₂ mixing ratio from the tropopause region to the lower thermosphere (∼950-1000 km), which drives a downward H₂ flux into Titan’s surface comparable to the H₂ escape flux out of the atmosphere (∼2 × 10¹⁰ cm⁻² s⁻¹ referenced to the surface) and requires larger photochemical production rates of H₂ than obtained by previous photochemical models. From detailed model calculations based on known photochemistry with eddy, molecular, and thermal diffusion, the tropospheric and thermospheric H₂ mole fractions are incompatible by a factor of ∼2. The measurements imply that the downward H₂ surface flux is in substantial excess of the speculative threshold value for methanogenic life consumption of H₂ (McKay, C.P., Smith, H.D. [2005], Icarus 178, 274-276. doi:10.1016/j.icarus.2005.05.018), but without the extreme reduction in the surface H₂ mixing ratio.     

Revised vertical cloud structure of Uranus from UKIRT/UIST observations and changes seen during Uranus’ Northern Spring Equinox from 2006 to 2008: Application of new methane absorption data and comparison with Neptune

Long-slit spectroscopy observations of Uranus by the United Kingdom InfraRed Telescope UIST instrument in 2006, 2007 and 2008 have been used to monitor the change in Uranus’ vertical and latitudinal cloud structure through the planet’s Northern Spring Equinox in December 2007.These spectra were analysed and presented by Irwin et al. (Irwin, P.G.J., Teanby, N.A., Davis, G.R. [2009]. Icarus 203, 287-302), but since publication, a new set of methane absorption data has become available (Karkoschka, E., Tomasko, M. [2010]. Methane absorption coefficients for the jovian planets from laboratory, Huygens, and HST data. Icarus 205, 674-694.), which appears to be more reliable at the cold temperatures and high pressures of Uranus’ deep atmosphere. We have fitted k-coefficients to these new methane absorption data and we find that although the latitudinal variation and inter-annual changes reported by Irwin et al. (2009) stand, the new k-data place the main cloud deck at lower pressures (2-3 bars) than derived previously in the H-band of ∼3-4 bars and ∼3 bars compared with ∼6 bars in the J-band. Indeed, we find that using the new k-data it is possible to reproduce satisfactorily the entire observed centre-of-disc Uranus spectrum from 1 to 1.75 μm with a single cloud at 2-3 bars provided that we make the particles more back-scattering at wavelengths less than 1.2 μm by, for example, increasing the assumed single-scattering albedo from 0.75 (assumed in the J and H-bands) to near 1.0. In addition, we find that using a deep methane mole fraction of 4% in combination with the associated warm ‘F’ temperature profile of Lindal et al. (Lindal, G.F., Lyons, J.R., Sweetnam, D.N., Eshleman, V.R., Hinson, D.P. [1987]. J. Geophys. Res. 92, 14987-15001), the retrieved cloud deck using the new (Karkoschka and Tomasko, 2010) methane absorption data moves to between 1 and 2 bars.The same methane absorption data and retrieval algorithm were applied to observations of Neptune made during the same programme and we find that we can again fit the entire 1-1.75 μm centre-of-disc spectrum with a single cloud model, providing that we make the stratospheric haze particles (of much greater opacity than for Uranus) conservatively scattering (i.e. ω = 1) and we also make the deeper cloud particles, again at around the 2 bar level more reflective for wavelengths less than 1.2 μm. Hence, apart from the increased opacity of stratospheric hazes in Neptune’s atmosphere, the deeper cloud structure and cloud composition of Uranus and Neptune would appear to be very similar.     

Minimum estimates of the amount and timing of gases released into the martian atmosphere from volcanic eruptions

Volcanism has been a major process during most of the geologic history of Mars. Based on data collected from terrestrial basaltic eruptions, we assume that the volatile content of martian lavas was typically ∼0.5 wt.% water, ∼0.7 wt.% carbon dioxide, ∼0.14 wt.% sulfur dioxide, and contained several other important volatile constituents. From the geologic record of volcanism on Mars we find that during the late Noachian and through the Amazonian volcanic degassing contributed ∼0.8 bar to the martian atmosphere. Because most of the outgassing consisted of greenhouse gases (i.e., CO₂ and SO₂) warmer surface temperatures resulting from volcanic eruptions may have been possible. Our estimates suggest that ∼1.1 × 10²¹ g (∼8 ± 1 m m⁻²) of juvenile water were released by volcanism; slightly more than half the amount contained in the north polar cap and atmosphere. Estimates for released CO₂ (1.6 × 10²¹ g) suggests that a large reservoir of carbon dioxide is adsorbed in the martian regolith or alternatively ∼300 cm cm⁻² of carbonates may have formed, although these materials would not occur readily in the presence of excess SO₂. Up to ∼120 cm cm⁻² (2.2 × 10²⁰ g) of acid rain (H₂SO₄) may have precipitated onto the martian surface as the result of SO₂ degassing. The hydrogen flux resulting from volcanic outgassing may help explain the martian atmospheric D/H ratio. The amount of outgassed nitrogen (∼1.3 mbar) may also be capable of explaining the martian atmospheric ¹⁵N/¹⁴N ratio. Minor gas constituents (HF, HCl, and H₂S) could have formed hydroxyl salts on the surface resulting in the physical weathering of geologic materials. The amount of hydrogen fluoride emitted (1.82 × 10¹⁸ g) could be capable of dissolving a global layer of quartz sand ∼5 mm thick, possibly explaining why this mineral has not been positively identified in spectral observations. The estimates of volcanic outgassing presented here will be useful in understanding how the martian atmosphere evolved over time.     

Ground-based measurements of the methane distribution on Titan

Using spectra taken with NIRSPEC (Near Infrared Spectrometer) and adaptive optics on the Keck II telescope, we resolved the latitudinal variation of the 3ν₂ band of CH₃D at 1.56 μm. As CH₃D is less abundant than CH₄ by a factor of 50±10×10^-5, these CH₃D lines do not saturate in Titan’s atmosphere, and are well characterized by laboratory measurements. Thus they do not suffer from the large uncertainties of the CH₄ lines that are weak enough to be unsaturated in Titan. Our measurements of the methane abundance are confined to the latitude range of 32°S-18°N and longitudes sampled by a 0.04″ slit centered at ∼195°W. The methane abundance below 10 km is constant to within 20% in the tropical atmosphere sampled by our observations, consistent with the low surface insolation and lack of surface methane [Griffith, C.A., McKay, C.P., Ferri, F., 2008. Astrophys. J. 687, L41-L44].     

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⁻².     

The source of widespread 3-μm absorption in Jupiter’s clouds: Constraints from 2000 Cassini VIMS observations

The Cassini flyby of Jupiter in 2000 provided spatially resolved spectra of Jupiter’s atmosphere using the Visual and Infrared Mapping Spectrometer (VIMS). A prominent characteristic of these spectra is the presence of a strong absorption at wavelengths from about 2.9 μm to 3.1 μm, previously noticed in a 3-μm spectrum obtained by the Infrared Space Observatory (ISO) in 1996. While Brooke et al. (Brooke, T.Y., Knacke, R.F., Encrenaz, T., Drossart, P., Crisp, D., Feuchtgruber, H. [1998]. Icarus 136, 1-13) were able to fit the ISO spectrum very well using ammonia ice as the sole source of particulate absorption, Sromovsky and Fry (Sromovsky, L.A., Fry, P.M. [2010]. Icarus 210, 211-229), using significantly revised NH₃ gas absorption models, showed that ammonium hydrosulfide (NH₄SH) provided a better fit to the ISO spectrum than NH₃, but that the best fit was obtained when both NH₃ and NH₄SH were present in the clouds. Although the large FOV of the ISO instrument precluded identification of the spatial distribution of these two components, the VIMS spectra at low and intermediate phase angles show that 3-μm absorption is present in zones and belts, in every region investigated, and both low- and high-opacity samples are best fit with a combination of NH₄SH and NH₃ particles at all locations. The best fits are obtained with a layer of small ammonia-coated particles (r ∼ 0.3 μm) overlying but often close to an optically thicker but still modest layer of much larger NH₄SH particles (r ∼ 10 μm), with a deeper optically thicker layer, which might also be composed of NH₄SH. Although these fits put NH₃ ice at pressures less than 500 mb, this is not inconsistent with the lack of prominent NH₃ features in Jupiter’s longwave spectrum because the reflectivity of the core particles strongly suppresses the NH₃ absorption features, at both near-IR and thermal wavelengths. Unlike Jupiter, Saturn lacks the broad 3-μm absorption feature, but does exhibit a small absorption near 2.965 μm, which resembles a similar jovian feature and suggests that both planets contain upper tropospheric clouds of sub-micron particles containing ammonia as a minor fraction.     

Determination of the complex refractive indices of Titan haze analogs using photothermal deflection spectroscopy

The spectrometers of the Cassini mission to the Saturn system have detected haze layers reaching up to 800 km in Titan’s atmosphere. Knowledge of the complex refractive index (k) of the haze is important for modeling the surface and atmosphere of Titan and retrieving some information about the functional groups present in the aerosols. Plasma discharges or ultraviolet radiation are commonly used to drive the formation of solid organics assumed to be good analogs of the Titan aerosols. [Tran, B.N., Ferris, J.P., Chera, J.J., 2003a. The photochemical formation of a Titan haze analog. Structural analysis by X-ray photoelectron and infrared spectroscopy. Icarus 162, 114-124; Tran, B.N., Force, M., Briggs, R., Ferris J.P., Persans, P., Chera, J.J., 2008. Photochemical processes on Titan: Irradiation of mixtures of gases that simulate Titan’s atmosphere. Icarus 177, 106-115] reported the index of refraction of analogs synthesized by far ultraviolet irradiation of various gas mixtures. k was determined in the 200-800 nm wavelength range from transmission and reflection spectroscopy. However, this technique is limited by (i) uncertainties in the absorption values because of the small amounts of organics available, (ii) light scattering by the surface roughness and particulates in the sample. These limitations prompted us to perform new measurements using photothermal deflection spectroscopy (PDS), a technique based on the conversion of absorbed light into heat in the material of interest. By combining traditional spectroscopy (λ < 500 nm) and PDS (λ > 500 nm), we determined values of k over the 375-1550 nm range. k values as low as 10⁻⁴ above 1000 nm were determined. This is one order of magnitude lower than the measurements generally used as a reference for Titan’s aerosols analogs [Khare, B.N., Sagan, C., Arakawa, E.T., Suits, F., Callicott, T.A., Williams, M.W., 1984. Optical-constants of organic Tholins produced in a simulated Titanian atmosphere—from soft-X-ray to microwave-frequencies. Icarus 60(1), 127-137]. We recommend that these results were used in models to describe the optical properties of the aerosols produced in Titan’s stratosphere.     

Venus night airglow: Ground-based detection of OH, observations of O2 emissions, and photochemical model

Venus nightglow was observed at NASA IRTF using a high-resolution long-slit spectrograph CSHELL at LT = 21:30 and 4:00 on Venus. Variations of the O₂ airglow at 1.27 μm and its rotational temperature are extracted from the observed spectra. The mean O₂ nightglow is 0.57 MR at 21:30 at 35°S-35°N, and the temperature increases from 171 K near the equator to ∼200 K at ±35°. We have found a narrow window that covers the OH (1-0) P1(4.5) and (2-1) Q1(1.5) airglow lines. The detected line intensities are converted into the (1-0) and (2-1) band intensities of 7.2 ± 1.8 kR and <1.4 kR at 21:30 and 15.5 ± 2 kR and 4.7 ± 1 kR at 4:00. The f-component of the (1-0) P1(4.5) line has not been detected in either observation, possibly because of resonance quenching in CO₂. The observed Earth’s OH (1-0) and (2-1) bands were 400 and 90 kR at 19:30 and 250 and 65 kR at 9:40, respectively. A photochemical model for the nighttime atmosphere at 80-130 km has been made. The model involves 61 reactions of 24 species, including odd hydrogen and chlorine chemistries, with fluxes of O, N, and H at 130 km as input parameters. To fit the OH vibrational distribution observed by VEX, quenching of OH (v > 3) in CO₂ only to v ? 2 is assumed. According to the model, the nightside-mean O₂ emission of 0.52 MR from the VEX and our observations requires an O flux of 2.9 × 10¹² cm⁻² s⁻¹ which is 45% of the dayside production above 80 km. This makes questionable the nightside-mean O₂ intensities of ∼1 MR from some observations. Bright nightglow patches are not ruled out; however, the mean nightglow is ∼0.5 MR as observed by VEX and supported by the model. The NO nightglow of 425 R needs an N flux of 1.2 × 10⁹ cm⁻² s⁻¹, which is close to that from VTGCM at solar minimum. However, the dayside supply of N at solar maximum is half that required to explain the NO nightglow in the PV observations. The limited data on the OH nightglow variations from the VEX and our observations are in reasonable agreement with the model. The calculated intensities and peak altitudes of the O₂, NO, and OH nightglow agree with the observations. Relationships for the nightglow intensities as functions of the O, N, and H fluxes are derived.     

An investigation of Pluto’s troposphere using stellar occultation light curves and an atmospheric radiative-conductive-convective model

We use a radiative-conductive-convective model to assess the height of Pluto’s troposphere, as well as surface pressure and surface radius, from stellar occultation data from the years 1988, 2002, and 2006. The height of the troposphere, if it exists, is less than 1 km for all years analyzed. Pluto has at most a planetary boundary layer and not a troposphere. As in previous analyses of Pluto occultation light curves, we find that the surface pressure is increasing with time, assuming that latitude and longitude variations in Pluto’s atmosphere are negligible. The surface pressure is found to be slightly higher ( μbar in 1988,  μbar in 2002, and 18.5 ± 4.7 μbar in 2006) than in our previous analyses with the troposphere excluded. The surface radius is determined to be . Comparison of the minimum reduced chi-squared values between the best-fit radiative-conductive-convective (i.e., troposphere-included) model and best-fit radiative-conductive (i.e., troposphere-excluded) shows that the troposphere-included model is only a slightly better fit to the data for all 3 years. Uncertainties in the small-scale physical processes of Pluto’s lower atmosphere and consequently the functional form of the model troposphere lend more confidence to the troposphere-excluded results.     

Scale-dependent infrared radiative damping rates on Mars and their role in the deposition of gravity-wave momentum flux

Using a Curtis-matrix model of 15 μm CO₂ radiative cooling rates for the martian atmosphere, we have computed vertical scale-dependent IR radiative damping rates from 0 to 200 km altitude over a broad band of vertical wavenumbers ∣m∣ = 2π(1-500 km)⁻¹ for representative meteorological conditions at 40°N and average levels of solar activity and dust loading. In the middle atmosphere, infrared (IR) radiative damping rates increase with decreasing vertical scale and peak in excess of 30 days⁻¹ at ∼50-80 km altitude, before gradually transitioning to scale-independent rates above ∼100 km due to breakdown of local thermodynamic equilibrium. We incorporate these computed IR radiative damping rates into a linear anelastic gravity-wave model to assess the impact of IR radiative damping, relative to wave breaking and molecular viscosity, in the dissipation of gravity-wave momentum flux. The model results indicate that IR radiative damping is the dominant process in dissipating gravity-wave momentum fluxes at ∼0-50 km altitude, and is the dominant process at all altitudes for gravity waves with vertical wavelengths ?10-15 km. Wave breaking becomes dominant at higher altitudes only for “fast” waves of short horizontal and long vertical wavelengths. Molecular viscosity plays a negligible role in overall momentum flux deposition. Our results provide compelling evidence that IR radiative damping is a major, and often dominant physical process controlling the dissipation of gravity-wave momentum fluxes on Mars, and therefore should be incorporated into future parameterizations of gravity-wave drag within Mars GCMs. Lookup tables for doing so, based on the current computations, are provided.     

Core-collapse supernovae and the Younger Dryas/terminal Rancholabrean extinctions

Early predictions that some supernovae release large quantities of prompt high energy photons are now corroborated by optical identification of core-collapse supernovae associated with extragalactic GRBS (beamed γ-ray bursts) and XRFS (beamed or un-beamed X-ray flashes). Given the in-galaxy supernova frequency and GRB and XRF recurrence statistics, significant Earth-incident events during the past several million years very likely occurred and nearby events should have affected the Earth and other planetary atmospheres, including terrestrial surface solar UV, the Earth’s climate, and its ecology. The Younger Dryas Stadial (∼12,900 to 11,550 calendar yr BP) began with sharply cooler temperatures in the Earth’s northern hemisphere, regional drought, paleoecological evidence compatible with increased UV, and abrupt increases in cosmogenic ¹⁴C and ¹⁰Be in ice and marine cores and tree rings. In North America, stratigraphic and faunal sequences indicate that a major pulse of mammalian extinctions (at least 23-31 genera) began very close to 12,830 calendar yr BP and was sudden: deposits one century younger are devoid of diverse extinct fauna remains. A 10 s beamed GRB within 2 kpc of the Earth delivers 100 kJ m⁻² fluence to the Earth’s atmosphere, where it causes spallation and catalytic reactions depleting 35-50% O³, and producing excess NO_x species (which favor cooling, drought, and surface fertility), ¹⁴C, and ¹⁰Be. An un-beamed, 10⁵⁰ erg hard photon impulse at ∼250 pc produces similar terrestrial atmospheric effects. A well-characterized massive star supernova, the unusually close Vela event (d = 250 ± 30 pc; total energy of 1-2 × 10⁵¹ erg; age constrained from remnant nebula shock velocities considerations at 13,000-16,000 yr and from the pulsar characteristic age at ∼11,400 yr) may have initiated the Younger Dryas climate change, and caused the extinction of the terminal Rancholabrean fauna.     

Combined experimental and theoretical studies on methane photolysis at 121.6 and 248 nm—implications on a program of laboratory simulations of Titan’s atmosphere

Methane is, together with N₂, the main precursor of Titan’s atmospheric chemistry. In our laboratory, we are currently developing a program of laboratory simulations of Titan’s atmosphere, where methane is intended to be dissociated by multiphotonic photolysis at 248 nm. A preliminary study has shown that multiphotonic absorption of methane at 248 nm is efficient and leads to the production of hydrocarbons such as C₂H₂ (Romanzin et al., 2008). Yet, at this wavelength, little is known about the branching ratios of the hydrocarbon radicals (CH₃, CH₂ and CH) and their following photochemistry. This paper thus aims at investigating methane photochemistry at 248 nm by comparing the chemical evolution observed after irradiation of CH₄ at 248 and at 121.6 nm (Ly-α). It is indeed important to see if the chemistry is driven the same way at both wavelengths in particular because, on Titan, methane photolysis mainly involves Ly-α photons. An approach combining experiments and theoretical analysis by means of a specifically adapted 0-D model has thus been developed and is presented in this paper. The results obtained clearly indicate that the chemistry is different depending on the wavelength. They also suggest that at 248 nm, methane dissociation is in competition with ionisation, which could occur through a three-photon absorption process. As a consequence, 248 nm photolysis appears to be unsuitable to study methane neutral photochemistry alone. The implications of this result on our laboratory simulation program and new experimental developments are discussed. Additional information on methane photochemistry at 121.6 nm are also obtained.     

Is there methane on Mars?

There have been several reports of methane on Mars at the 10-60 ppbv level. Most suggest that methane is both seasonally and latitudinally variable. Here we review why variable methane on Mars is physically and chemically implausible, and then we critically review the published reports. There is no known mechanism for destroying methane chemically on Mars. But if there is one, methane oxidation would deplete the O₂ in Mars’s atmosphere in less than 10,000 years unless balanced by an equally large unknown source of oxidizing power. Physical sequestration does not raise these questions, but adsorption in the regolith or condensation in clathrates ignore competition for adsorption sites or are inconsistent with clathrate stability, respectively. Furthermore, any mechanism that relies on methane’s van der Waals’ attraction is inconsistent with the continued presence of Xe in the atmosphere at the 60 ppbv level. We then use the HITRAN database and transmission calculations to identify and characterize the absorption lines that would be present on Earth or Mars at the wavelengths of the published observations. These reveal strong competing telluric absorption that is most problematic at just those wavelengths where methane’s signature seems most clearly seen from Earth. The competing telluric lines must be removed with models. The best case for martian methane was made for the ¹²CH₄ν₃ R0 and R1 lines seen in blueshift when Mars was approaching Earth in early 2003 (Mumma, M.J., Villanueva, G.L., Novak, R.E., Hewagama, T., Bonev, B.P., DiSanti, M.A., Mandell, A.M., Smith, M.D. [2009]. Science 323, 1041-1045). For these the Doppler shift moves the two martian lines into near coincidence with telluric ¹³CH₄ν₃ R1 and R2 lines that are 10-50× stronger than the inferred martian lines. By contrast, the ¹²CH₄ν₃ R0 and R1 lines when observed in redshift do not contend with telluric ¹³CH₄. For these lines, Mumma et al.’s observations and analyses are consistent with an upper limit on the order of 3 ppbv.     

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.     

Monte Carlo modeling of Io’s [OI] 6300 Å and [SII] 6716 Å auroral emission in eclipse

We present a Monte Carlo (MC) model of [OI] 6300 Å and [SII] 6716 Å emission from Io entering eclipse. The simulation accounts for the 3-D distribution of SO₂, O, SO, S, and O₂ in Io’s atmosphere, several volcanic plumes, and the magnetic field around Io. Thermal electrons from the jovian plasma torus are input along the simulation domain boundaries and move along the magnetic field lines distorted by Io, occasionally participating in collisions with neutrals. We find that the atmospheric asymmetry resulting from varying degrees of atmospheric collapse across Io (due to eclipse ingress) and the presence of volcanoes contributes significantly to the unique morphology of the [OI] 6300 Å emission. The [OI] radiation lifetime of ∼134 s limits the emission to regions that have a sufficiently low neutral density so that intermolecular collisions are rare. We find that at low altitudes (typically <40 km) and in volcanic plumes (Pele, Prometheus, etc.) the number density is large enough (>4 × 10⁹ cm⁻³) to collisionally quench nearly all (>95%) of the excited oxygen for reasonable quenching efficiencies. Upstream (relative to the plasma flow), Io’s perturbation of the jovian magnetic field mirrors electrons with high pitch angles, while downstream collisions can trap the electrons. This magnetic field perturbation is one of the main physical mechanisms that results in the upstream/downstream brightness asymmetry in [OI] emission seen in the observation by Trauger et al. (Trauger, J.T., Stapelfeldt, K.R., Ballester, G.E., Clarke, J.I., 1997. HST observations of [OI] emissions from Io in eclipse. AAS-DPS Abstract (1997DPS29.1802T)). There are two other main causes for the observed brightness asymmetry. First, the observation’s viewing geometry of the wake spot crosses the dayside atmosphere and therefore the wake’s observational field of view includes higher oxygen column density than the upstream side. Second, the phased entry into eclipse results in less atmospheric collapse and thus higher collisional quenching on the upstream side relative to the wake. We compute a location (both in altitude and latitude) for the intense wake emission feature that agrees reasonably well with this observation. Furthermore, the peak intensity of the simulated wake feature is less than that observed by a factor of ∼3, most likely because our model does not include direct dissociation-excitation of SO₂ and SO. We find that the latitudinal location of the emission feature depends not so much on the tilt of the magnetic field as on the relative north/south flux tube depletion that occurs due to Io’s changing magnetic latitude in the plasma torus. From 1-D simulations, we also find that the intensity of [SII] 6716 and 6731 Å emission is much weaker than that of [OI] even if the [SII] excitation cross section is 10³ times larger than excitation to [OI]. This is because the density of S⁺ is much less than that of O and because the Einstein-A coefficient of the [SII] emission is a factor of ∼10 smaller than that of [OI].     

Simulation of Io’s auroral emission: Constraints on the atmosphere in eclipse

We study the morphology of Io’s aurora by comparing simulation results of a three-dimensional (3D) two-fluid plasma model to observations by the high-resolution Long-Range Reconnaissance Imager (LORRI) on-board the New Horizons spacecraft and by the Hubble Space Telescope Advanced Camera for Surveys (HST/ACS). In 2007, Io’s auroral emission in eclipse has been observed simultaneously by LORRI and ACS and the observations revealed detailed features of the aurora, such as a huge glowing plume at the Tvashtar paterae close to the North pole. The auroral radiation is generated in Io’s atmosphere by collisions between impinging magnetospheric electrons and various neutral gas components. We calculate the interaction of the magnetospheric plasma with Io’s atmosphere-ionosphere and simulate the auroral emission. Our aurora model takes into account not only the direct influence of the atmospheric distribution on the morphology and intensity of the emission, but also the indirect influence of the atmosphere on the plasma environment and thus on the exciting electrons. We find that the observed morphology in eclipse can be explained by a smooth (non-patchy) equatorial atmosphere with a vertical column density that corresponds to ∼10% of the column density of the sunlit atmosphere. The atmosphere is asymmetric with two times higher density and extension on the downstream hemisphere. The auroral emission from the Tvashtar volcano enables us to constrain the plume gas content for the first time. According to our model, the observed intensity of the Tvashtar plume implies a mean column density of ∼5 × 10¹⁵ cm⁻² for the plume region.     

Negative ion chemistry in Titan's upper atmosphere

The Electron Spectrometer (ELS), one of the sensors making up the Cassini Plasma Spectrometer (CAPS) revealed the existence of numerous negative ions in Titan's upper atmosphere. The observations at closest approach (∼1000 km) show evidence for negatively charged ions up to ∼10,000 amu/q, as well as two distinct peaks at 22±4 and 44±8 amu/q, and maybe a third one at 82±14 amu/q. We present the first ionospheric model of Titan including negative ion chemistry. We find that dissociative electron attachment to neutral molecules (mostly HCN) initiates the formation of negative ions. The negative charge is then transferred to more acidic molecules such as HC₃N, HC₅N or C₄H₂. Loss occurs through associative detachment with radicals (H and CH₃). We attribute the three low mass peaks observed by ELS to CN⁻, C₃N⁻/C₄H⁻ and C₅N⁻. These species are the first intermediates in the formation of the even larger negative ions observed by ELS, which are most likely the precursors to the aerosols observed at lower altitudes.     

Radar measurements of Mercury’s north pole at 70 cm wavelength

We present radar imaging of Mercury using the Arecibo Observatory’s 70-cm wavelength radar system during the inferior conjunction of July 1999. At that time the sub-Earth latitude was ∼11°N and the highly reflective region at Mercury’s north pole that was first identified in radar images at the shorter wavelengths of 3.6 cm [Slade, M.A., Butler, B.J., Muhleman, D.O., 1992. Science 258, 635-640] and 13 cm [Harmon, J.K., Slade, M.A., 1992. Science 258, 640-643] was again clearly detected. The reflectivity averaged over a 75,000 km² region including the pole is similar to that measured at the other wavelengths over a comparable area, and the 70 cm circular polarization ratio of μ_C∼0.87 is possibly slightly lower. If this strong backscattering results from volume scattering in low absorption layers, the persistence of this effect over more than an order of magnitude change in wavelength scale has implications for the depth and thickness of the deposits responsible. The resolution of the radar maps at this wavelength is not sufficient to resolve individual craters, nor to discern features at other latitudes, but the planet’s total reflectivity is consistent with previous work and the scattering function suggests a surface roughness at this wavelength similar to the lunar highlands.