Impact of electron chemistry on the structure and composition of Io's atmosphere

Two-dimensional model calculations (altitude and solar zenith angle) are performed to investigate the impact of electron chemistry on the composition and structure of Io's atmosphere. The calculations are based upon the model of Wong and Smyth (2000, Icarus 146, 60-74) for Io's SO₂ sublimation atmosphere with the addition of new electron chemistry, where the interactions of the electrons and neutrals are treated in a simple fashion. The model calculations are presented for Io's atmosphere at western elongation (dusk ansa) for both a low-density case (subsolar temperature of 113 K) and a high-density case (subsolar temperature of 120 K). The impact of electron-neutral chemistry on the composition and structure of Io's atmosphere is confined primarily to an interaction layer. The penetration depth of the interaction layer is limited to high altitudes in the thicker dayside atmosphere but reaches the surface in the thinner dayside and/or nightside atmosphere at larger solar zenith angles. Within most of the thicker dayside atmosphere, the column density of SO₂ is not significantly altered by electrons, but in the interaction layer all number densities are significantly altered: SO₂ is reduced, O, SO, S, and O₂ are greatly enhanced, and O, SO, and S become comparable to SO₂ at high altitudes. For the thinner nightside atmosphere, the species number densities are dramatically altered: SO₂ is drastically reduced to the least abundant species of the SO₂ family, SO and O₂ are significantly reduced at all altitudes, and O and S are dramatically enhanced and become the dominant species at all altitudes except near the surface. The interaction layer also defines the location of the emission layer for neutrals excited by electron impact and hence determines the fraction of the total neutral column density that is visible in remote observation. Electron chemistry may also impact the ratio of the equatorial to polar SO₂ column density deduced from Lyman-α images and the north-south alternating and System III longitude-dependent asymmetry observed in polar O and S emissions.     

Cassini observations of Io's visible aurorae

More than 500 images of Io in eclipse were acquired by the Cassini spacecraft in late 2000 and early 2001 as it passed through the jovian system en route to Saturn (Porco et al., 2003, Science 299, 1541-1547). Io's bright equatorial glows were detected in Cassini's near-ultraviolet filters, supporting the interpretation that the visible emissions are predominantly due to molecular SO₂. Detailed comparisons of laboratory SO₂ spectra with the Cassini observations indicate that a mixture of gases contribute to the equatorial emissions. Potassium is suggested by new detections of the equatorial glows at near-infrared wavelengths from 730 to 800 nm. Neutral atomic oxygen and sodium are required to explain the brightness of the glows at visible wavelengths. The molecule S₂ is postulated to emit most of the glow intensity in the wavelength interval from 390 to 500 nm. The locations of the visible emissions vary in response to the changing orientation of the external magnetic field, tracking the tangent points of the jovian magnetic field lines. Limb glows distinct from the equatorial emissions were observed at visible to near-infrared wavelengths from 500 to 850 nm, indicating that atomic O, Na, and K are distributed across Io's surface. Stratification of the atmosphere is demonstrated by differences in the altitudes of emissions at various wavelengths: SO₂ emissions are confined to a region close to Io's surface, whereas neutral oxygen emissions are seen at altitudes that reach up to 900 km, or half the radius of the satellite. Pre-egress brightening demonstrates that light scattered into Jupiter's shadow by gases or aerosols in the giant planet's upper atmosphere contaminates images of Io taken within 13 minutes of entry into or emergence from Jupiter's umbra. Although partial atmospheric collapse is suggested by the longer timescale for post-ingress dimming than pre-egress brightening, Io's atmosphere must be substantially supported by volcanism to retain auroral emissions throughout the duration of eclipse.     

1-D DSMC simulation of Io's atmospheric collapse and reformation during and after eclipse

A one-dimensional Direct Simulation Monte Carlo (DSMC) model is used to examine the effects of a non-condensable species on Io's sulfur dioxide sublimation atmosphere during eclipse and just after egress. Since the vapor pressure of SO₂ is extremely sensitive to temperature, the frost-supported dayside sublimation atmosphere had generally been expected to collapse during eclipse as the surface temperature dropped. For a pure SO₂ atmosphere, however, it was found that during the first 10 min of eclipse, essentially no change in the atmospheric properties occurs at altitudes above ∼100 km due to the finite ballistic/acoustic time. Hence immediately after ingress the auroral emission morphology above 100 km should resemble that of the immediate pre-eclipse state. Furthermore, the collapse dynamics are found to be greatly altered by the presence of even a small amount of a non-condensable species which forms a diffusion layer near the surface that prevents rapid collapse. It is found that after 10 min essentially no collapse has occurred at altitudes above ∼20 km when a nominal mole fraction of non-condensable gas is present. Collapse near the surface occurs relatively quickly until a static diffusion layer many mean free paths thick of the non-condensable gas builds up which then retards further collapse of the SO₂ atmosphere. For example, for an initial surface temperature of 110 K and 35% non-condensable mole-fraction, the ratio of the SO₂ column density to the initial column density was found to be 0.73 after 10 min, 0.50 after 30 min, and 0.18 at the end of eclipse. However, real gas species (SO, O₂) may not be perfectly non-condensable at Io's surface temperatures. If the gas species was even weakly condensable (non-zero sticking/reaction coefficient) then the effect of the diffusion layer on the dynamics was dramatically reduced. In fact, if the sticking coefficient of the non-condensable exceeds ∼0.25, the collapse dynamics are effectively the same as if there were no non-condensable present. This sensitivity results because the loss of non-condensable to the surface reduces the effective diffusion layer size, and the formation of an effective diffusion layer requires that the layer be stationary; this does not occur if the surface is a sink. Upon egress, vertical stratification of the condensable and non-condensable species occurs, with the non-condensable species being lifted (or pushed) to higher altitudes by the sublimating SO₂ after the sublimating atmosphere becomes collisional. Stratification should affect the morphology and intensity of auroral glows shortly after egress.     

Mid-infrared detection of large longitudinal asymmetries in Io's SO2 atmosphere

We have observed about 16 absorption lines of the ν₂ SO₂ vibrational band on Io, in disk-integrated 19-μm spectra taken with the TEXES high spectral resolution mid-infrared spectrograph at the NASA Infrared Telescope Facility in November 2001, December 2002, and January 2004. These are the first ground-based infrared observations of Io's sunlit atmosphere, and provide a new window on the atmosphere that allows better longitudinal and temporal monitoring than previous techniques. Dramatic variations in band strength with longitude are seen that are stable over at least a 2 year period. The depth of the strongest feature, a blend of lines centered at 530.42 cm⁻¹, varies from about 7% near longitude 180° to about 1% near longitude 315° W, as measured at a spectral resolution of 57,000. Interpretation of the spectra requires modeling of surface temperatures and atmospheric density across Io's disk, and the variation in non-LTE ν₂ vibrational temperature with altitude, and depends on the assumed atmospheric and surface temperature structure. About half of Io's 19-μm radiation comes from the Sun-heated surface, and half from volcanic hot spots with temperatures primarily between 150 and 200 K, which occupy about 8% of the surface. The observations are thus weighted towards the atmosphere over these low-temperature hot spots. If we assume that the atmosphere over the hot spots is representative of the atmosphere elsewhere, and that the atmospheric density is a function of latitude, the most plausible interpretation of the data is that the equatorial atmospheric column density varies from about 1.5×¹⁰¹⁷ cm⁻² near longitude 180° W to about 1.5×¹⁰¹⁶ cm⁻² near longitude 300° W, roughly consistent with HST UV spectroscopy and Lyman-α imaging. The inferred atmospheric kinetic temperature is less than about 150 K, at least on the anti-Jupiter hemisphere where the bands are strongest, somewhat colder than inferred from HST UV spectroscopy and millimeter-wavelength spectroscopy. This longitudinal variability in atmospheric density correlates with the longitudinal variability in the abundance of optically thick, near-UV bright SO₂ frost. However it is not clear whether the correlation results from volcanic control (regions of large frost abundance result from greater condensation of atmospheric gases supported by more vigorous volcanic activity in these regions) or sublimation control (regions of large frost abundance produce a more extensive atmosphere due to more extensive sublimation). Comparison of data taken in 2001, 2002, and 2004 shows that with the possible exception of longitudes near 180° W between 2001 and 2002, Io's atmospheric density does not appear to decrease as Io recedes from the Sun, as would be expected if the atmosphere were supported by the sublimation of surface frost, suggesting that the atmosphere is dominantly supported by direct volcanic supply rather than by frost sublimation. However, other evidence such as the smooth variation in atmospheric abundance with latitude, and atmospheric changes during eclipse, suggest that sublimation support is more important than volcanic support, leaving the question of the dominant atmospheric support mechanism still unresolved.     

A comprehensive numerical simulation of Io’s sublimation-driven atmosphere

Io’s sublimation-driven atmosphere is modeled using the direct simulation Monte Carlo (DSMC) method. These rarefied gas dynamics simulations improve upon earlier models by using a three-dimensional domain encompassing the entire planet computed in parallel. The effects of plasma heating, planetary rotation, inhomogeneous surface frost, molecular residence time of SO₂ on the exposed (non-volatile) rocky surface, and surface temperature distribution are investigated. Circumplanetary flow is predicted to develop from the warm dayside toward the cooler nightside. Io’s rotation leads to a highly asymmetric frost surface temperature distribution (due to the frost’s high thermal inertia) which results in circumplanetary flow that is not axi-symmetric about the subsolar point. The non-equilibrium thermal structure of the atmosphere, specifically vibrational and rotational temperatures, is also examined. Plasma heating is found to significantly inflate the atmosphere on both the dayside and nightside. The plasma energy flux causes high temperatures at high altitudes but plasma energy depletion through the dense gas column above the warmest frost permits gas temperatures cooler than the surface at low altitudes. A frost map (Douté, S., Schmitt, B., Lopes-Gautier, R., Carlson, R., Soderblom, L., Shirley, J., and the Galileo NIMS Team [2001]. Icarus 149, 107-132) is used to control the sublimated flux of SO₂ which can result in inhomogeneous column densities that vary by nearly a factor of four for the same surface temperature. A short residence time for SO₂ molecules on the “rock” component is found to smooth lateral atmospheric inhomogeneities caused by variations in the surface frost distribution, creating an atmosphere that looks nearly identical to one with uniform frost coverage. A longer residence time is found to agree better with mid-infrared observations (Spencer, J.R., Lellouch, E., Richter, M.J., López-Valverde, M.A., Jessup, K.L, Greathouse, T.K., Flaud, J. [2005]. Icarus 176, 283-304) and reproduce the observed anti-jovian/sub-jovian column density asymmetry. The computed peak dayside column density for Io assuming a surface frost temperature of 115 K agrees with those suggested by Lyman-α observations (Feaga, L.M., McGrath, M., Feldman, P.D. [2009]. Icarus 201, 570-584). On the other hand, the peak dayside column density at 120 K is a factor of five larger and is higher than the upper range of observations (Jessup, K.L., Spencer, J.R., Ballester, G.E., Howell, R.R., Roesler, F., Vigel, M., Yelle, R. [2004]. Icarus 169, 197-215; Spencer et al., 2005).     

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.     

Temporal behavior of the SO 1.707 μm ro-vibronic emission band in Io's atmosphere

We report observations of the ro-vibronic a¹Δ→X³Σ⁻ transition of SO at 1.707 μm on Io. These data were taken while Io was eclipsed by Jupiter, on four nights between July 2000 and March 2003. We analyze these results in conjunction with a previously published night to investigate the temporal behavior of these emissions. The observations were all conducted using the near-infrared spectrometer NIRSPEC on the W.M. Keck II telescope. The integrated emitted intensity for this band varies from 0.8×¹⁰²⁷ to 2.4×¹⁰²⁷ photons/s, with a possible link to variations in Loki's infrared brightness. The band-shapes imply rotational temperatures of 550-1000 K for the emitting gas, lending further evidence to a volcanic origin for sulfur monoxide. An attempt to detect the B¹Σ→X³Σ⁻ transition of SO at 0.97 μm was unsuccessful; simultaneous detection with the 1.707 μm band would permit determination of the SO column abundance.     

Multi-wavelength simulations of atmospheric radiation from Io with a 3-D spherical-shell backward Monte Carlo radiative transfer model

Conflicting observations regarding the dominance of either sublimation or volcanism as the source of the atmosphere on Io and disparate reports on the extent of its spatial distribution and the absolute column abundance invite the development of detailed computational models capable of improving our understanding of Io’s unique atmospheric structure and origin. Improving upon previous models, Walker et al. (Walker, A.C., Gratiy, S.L., Levin, D.A., Goldstein, D.B., Varghese, P.L., Trafton, L.M., Moore, C.H., Stewart, B. [2009]. Icarus) developed a fully 3-D global rarefied gas dynamics model of Io’s atmosphere including both sublimation and volcanic sources of SO₂ gas. The fidelity of the model is tested by simulating remote observations at selected wavelength bands and comparing them to the corresponding astronomical observations of Io’s atmosphere. The simulations are performed with a new 3-D spherical-shell radiative transfer code utilizing a backward Monte Carlo method. We present: (1) simulations of the mid-infrared disk-integrated spectra of Io’s sunlit hemisphere at 19 μm, obtained with TEXES during 2001-2004; (2) simulations of disk-resolved images at Lyman-α obtained with the Hubble Space Telescope (HST), Space Telescope Imaging Spectrograph (STIS) during 1997-2001; and (3) disk-integrated simulations of emission line profiles in the millimeter wavelength range obtained with the IRAM-30 m telescope in October-November 1999. We found that the atmospheric model generally reproduces the longitudinal variation in band depth from the mid-infrared data; however, the best match is obtained when our simulation results are shifted ∼30° toward lower orbital longitudes. The simulations of Lyman-α images do not reproduce the mid-to-high latitude bright patches seen in the observations, suggesting that the model atmosphere sustains columns that are too high at those latitudes. The simulations of emission line profiles in the millimeter spectral region support the hypothesis that the atmospheric dynamics favorably explains the observed line widths, which are too wide to be formed by thermal Doppler broadening alone.     

Ion cyclotron waves at Io: implications for the temporal variation of Io's atmosphere

When the flowing torus plasma encounters the upper atmosphere of Jupiter's moon, Io, newly created ions are rapidly accelerated by the motional electric field. Many of these ions are reneutralized and form a spray of fast neutrals that travel far away from Io before being reionized by photoionization and impact. These ions, now far from Io, are unstable to the generation of ion cyclotron waves. These waves in turn act as a mass spectrometer allowing Galileo magnetic measurements to be used to probe the composition of the atmosphere of Io and how it varies in time and in space. We now have six Galileo passes by Io on which we have measurements with sufficient cadence to examine the ion cyclotron waves. One of these passes, on Galileo's 32nd orbit has not been discussed previously. These passes provide sufficient observations to begin to distinguish the sources of variability. We find that while the atmosphere of Io varies temporally throughout the mission, it also has a spatial variation in composition at any instant of time.     

Relative contributions of sublimation and volcanoes to Io's atmosphere inferred from its plasma interaction during solar eclipse

We present a model that describes Io's delayed electrodynamic response to a temporal change in Io's atmosphere. Our model incorporates the relevant physical processes involved in Io's atmosphere-ionosphere-magnetosphere electrodynamic interaction to predict the far-ultraviolet (FUV) radiation as Io enters Jupiter's shadow and re-emerges into sunlight. The predicted FUV brightnesses are highly nonlinear as the strength of the electrodynamic interaction depends on the ratios of ionospheric conductances to the torus Alfvén conductance, but the former are functions of electrodynamics and the atmospheric density, which decays rapidly upon entering eclipse. Key factors governing the time evolution are the column density due to sublimation and the column density due to volcanoes, which maintain the background atmosphere during eclipse. The plasma interaction does not react instantaneously, but lags to a temporarily changing atmosphere. We find three qualitatively different scenarios with two of them including a post-eclipse brightening. The brightness ratio of in-sunlight/in-eclipse coupled with the existence of a sub-jovian equatorial spot constrains the volcanic column density to several times 10¹⁸ m⁻², based on the currently available observations. Thus in sunlight, the sublimation driven part of Io's atmosphere dominates the volcanically driven contribution by roughly a factor of 10 or more.     

Lava lakes on Io: observations of Io's volcanic activity from Galileo NIMS during the 2001 fly-bys

Galileo's Near-Infrared Mapping Spectrometer (NIMS) obtained its final observations of Io during the spacecraft's fly-bys in August (I31) and October 2001 (I32). We present a summary of the observations and results from these last two fly-bys, focusing on the distribution of thermal emission from Io's many volcanic regions that give insights into the eruption styles of individual hot spots. We include a compilation of hot spot data obtained from Galileo, Voyager, and ground-based observations. At least 152 active volcanic centers are now known on Io, 104 of which were discovered or confirmed by Galileo observations, including 23 from the I31 and I32 Io fly-by observations presented here. We modify the classification scheme of Keszthelyi et al. (2001, J. Geophys. Res. 106 (E12) 33 025-33 052) of Io eruption styles to include three primary types: promethean (lava flow fields emplaced as compound pahoehoe flows with small plumes <200 km high originating from flow fronts), pillanian (violent eruptions generally accompanied by large outbursts, >200 km high plumes and rapidly-emplaced flow fields), and a new style we call “lokian” that includes all eruptions confined within paterae with or without associated plume eruptions). Thermal maps of active paterae from NIMS data reveal hot edges that are characteristic of lava lakes. Comparisons with terrestrial analogs show that Io's lava lakes have thermal properties consistent with relatively inactive lava lakes. The majority of activity on Io, based on locations and longevity of hot spots, appears to be of this third type. This finding has implications for how Io is being resurfaced as our results imply that eruptions of lava are predominantly confined within paterae, thus making it unlikely that resurfacing is done primarily by extensive lava flows. Our conclusion is consistent with the findings of Geissler et al. (2004, Icarus, this issue) that plume eruptions and deposits, rather than the eruption of copious amounts of effusive lavas, are responsible for Io's high resurfacing rates. The origin and longevity of islands within ionian lava lakes remains enigmatic.     

Observations and temperatures of Io's Pele Patera from Cassini and Galileo spacecraft images

Pele has been the most intense high-temperature hotspot on Io to be continuously active during the Galileo monitoring from 1996-2001. A suite of characteristics suggests that Pele is an active lava lake inside a volcanic depression. In 2000-2001, Pele was observed by two spacecraft, Cassini and Galileo. The Cassini observations revealed that Pele is variable in activity over timescales of minutes, typical of active lava lakes in Hawaii and Ethiopia. These observations also revealed that the short-wavelength thermal emission from Pele decreases with rotation of Io by a factor significantly greater than the cosine of the emission angle, and that the color temperature becomes more variable and hotter at high emission angles. This behavior suggests that a significant portion of the visible thermal emission from Pele comes from lava fountains within a topographically confined lava body. High spatial resolution, nightside images from a Galileo flyby in October 2001 revealed a large, relatively cool (<800 K) region, ringed by bright hotspots, and a central region of high thermal emission, which is hypothesized to be due to fountaining and convection in the lava lake. Images taken through different filters revealed color temperatures of 1500±80 K from Cassini ISS data and 1605±220 and 1420±100 K from small portions of Galileo SSI data. Such temperatures are near the upper limit for basaltic compositions. Given the limitations of deriving lava eruption temperature in the absence of in situ measurement, it is possible that Pele has lavas with ultramafic compositions. The long-lived, vigorous activity of what is most likely an actively overturning lava lake in Pele Patera indicates that there is a strong connection to a large, stable magma source region.     

A primordial origin for the atmospheric methane of Saturn’s moon Titan

The origin of Titan’s atmospheric methane is a key issue for understanding the origin of the saturnian satellite system. It has been proposed that serpentinization reactions in Titan’s interior could lead to the formation of the observed methane. Meanwhile, alternative scenarios suggest that methane was incorporated in Titan’s planetesimals before its formation. Here, we point out that serpentinization reactions in Titan’s interior are not able to reproduce the deuterium over hydrogen (D/H) ratio observed at present in methane in its atmosphere, and would require a maximum D/H ratio in Titan’s water ice 30% lower than the value likely acquired by the satellite during its formation, based on Cassini observations at Enceladus. Alternatively, production of methane in Titan’s interior via radiolytic reactions with water can be envisaged but the associated production rates remain uncertain. On the other hand, a mechanism that easily explains the presence of large amounts of methane trapped in Titan in a way consistent with its measured atmospheric D/H ratio is its direct capture in the satellite’s planetesimals at the time of their formation in the solar nebula. In this case, the mass of methane trapped in Titan’s interior can be up to ∼1300 times the current mass of atmospheric methane.     

Enceladus' plume: Compositional evidence for a hot interior

The combined observations of Saturn's moon Enceladus by the Cassini CAPS, INMS and UVIS instruments detected water vapor geysers in which were present molecular nitrogen (N₂), carbon dioxide (CO₂), methane (CH₄), propane (C₃H₈), acetylene (C₂H₂), and several other species, together with all of the decomposition products of water. We propose that the presence of N₂ in the plume indicates thermal decomposition of ammonia, and hence high temperatures in the interior of the moon (e.g., 500 to 800 K). Such an environment also appears to be suitable for the production of methane (CH₄) from carbon monoxide (CO), or carbon dioxide (CO₂). The presence of C₂H₂ and C₃H₈ strongly suggest that catalytic reactions took place within a very hot environment. The internal environment of Enceladus is inferred to be or have been favorable for aqueous, catalytic chemistry. This permits the synthesis of many complex organic compounds that could be detected in future Cassini observations.     

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

The near-IR spectrum of Titan modeled with an improved methane line list

We have obtained spatially resolved spectra of Titan in the near-infrared J, H and K bands at a resolving power of ∼5000 using the near-infrared integral field spectrometer (NIFS) on the Gemini North 8 m telescope. Using recent data from the Cassini/Huygens mission on the atmospheric composition and surface and aerosol properties, we develop a multiple-scattering radiative transfer model for the Titan atmosphere. The Titan spectrum at these wavelengths is dominated by absorption due to methane with a series of strong absorption band systems separated by window regions where the surface of Titan can be seen. We use a line-by-line approach to derive the methane absorption coefficients. The methane spectrum is only accurately represented in standard line lists down to ∼2.1 μm. However, by making use of recent laboratory data and modeling of the methane spectrum we are able to construct a new line list that can be used down to 1.3 μm. The new line list allows us to generate spectra that are a good match to the observations at all wavelengths longer than 1.3 μm and allow us to model regions, such as the 1.55 μm window that could not be studied usefully with previous line lists such as HITRAN 2008. We point out the importance of the far-wing line shape of strong methane lines in determining the shape of the methane windows. Line shapes with Lorentzian, and sub-Lorentzian regions are needed to match the shape of the windows, but different shape parameters are needed for the 1.55 μm and 2 μm windows. After the methane lines are modeled our observations are sensitive to additional absorptions, and we use the data in the 1.55 μm region to determine a D/H ratio of 1.77 ± 0.20 × 10⁻⁴, and a CO mixing ratio of 50 ± 11 ppmv. In the 2 μm window we detect absorption features that can be identified with the ν₅ + 3ν₆ and 2ν₃ + 2ν₆ bands of CH₃D.     

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.     

Photometry of Triton 1992-2004: Surface volatile transport and discovery of a remarkable opposition surge

Triton, the large satellite of Neptune, was imaged by the Voyager 2 spacecraft in 1989 with dark plumes originating in its volatile-rich south polar region. Southern summer solstice, a time when seasonal volatile transport should be at a maximum, occurred in 2001. Ground-based observations of Triton’s rotational light curve obtained from Table Mountain Observatory in 2000-2004 reveal volatile transport on its surface. When compared with a static frost model constructed from Voyager images, the light curve shows an increase in total amplitude. An earlier light curve obtained in 1992 from Mauna Kea Observatory is consistent with the static frost model. This movement of volatiles on the surface agrees with recent imaging results from the Hubble Space Telescope (Bauer, J.M., Buratti, B.J., Li, J.-Y., Mosher, J.A., Hicks, M.D., Schmidt, B.E., Goguen, J.D. [2010]. Astrophys. J. 723, L49-L52). The changes in the light curve can be explained by the transport of nitrogen frost on the surface or by the uncovering of bedrock of less volatile methane. We also find that Triton exhibits a large opposition surge at solar phase angles less than 0.1°. This surge cannot be entirely explained by the effects of coherent backscatter.