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.     

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.     

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.     

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.     

The atmospheric signature of Io's Prometheus plume and anti-jovian hemisphere: evidence for a sublimation atmosphere

Using the Hubble Space Telescope's Space Telescope Imaging Spectrograph we have obtained for the first time spatially resolved 2000-3000 Å spectra of Io's Prometheus plume and adjoining regions on Io's anti-jovian hemisphere in the latitude range 60° N-60° S, using a 0.1″ slit centered on Prometheus and tilted roughly 45° to the spin axis. The SO₂ column density peaked at 1.25×10¹⁷ cm⁻² near the equator, with an additional 5×10¹⁶ cm⁻² enhancement over Prometheus corresponding to a model volcanic SO₂ output of 10⁵ kg s⁻¹. Apart from the Prometheus peak, the SO₂ column density dropped fairly smoothly away from the subsolar point, even over regions that included potential volcanic sources. At latitudes less than ±30°, the dropoff rate was consistent with control by vapor pressure equilibrium with surface frost with subsolar temperature 117.3±0.6 K, though SO₂ abundance was higher than predicted by vapor pressure control at mid-latitudes, especially in the northern hemisphere. We conclude that, at least at low latitudes on the anti-jovian hemisphere where there are extensive deposits of optically-thick SO₂ frost, the atmosphere is probably primarily supported by sublimation of surface frost. Although the 45° tilt of our slit prevents us from separating the dependence of atmospheric density on solar zenith angle from its dependence on latitude, the pattern is consistent with a sublimation atmosphere regardless of which parameter is the dominant control. The observed drop in gas abundance towards higher latitudes is consistent with the interpretation of previous Lyman alpha images of Io as indicating an atmosphere concentrated at low latitudes. Comparison with previous disk-resolved UV spectroscopy, Lyman-alpha images, and mid-infrared spectroscopy suggests that Io's atmosphere is denser and more widespread on the anti-jovian hemisphere than at other longitudes. SO₂ gas temperatures were in the range of 150-250 K over the majority of the anti-jovian hemisphere, consistent with previous observations. SO was not definitively detected in our spectra, with upper limits to the SO/SO₂ ratio in the range 1-10%, roughly consistent with previous observations. S₂ gas was not seen anywhere, with an upper limit of 7.5×10¹⁴ cm⁻² for the Prometheus plume, confirming that this plume is significantly poorer in S₂ than the Pele plume (S₂ /SO₂<0.005, compared to 0.08-0.3 at Pele). In addition to the gas absorption signatures, we have observed continuum emission in the near ultraviolet (near 2800 Å) for the first time. The brightness of the observed emission was directly correlated with the SO₂ abundance, strongly peaking in the equatorial region over Prometheus. Emission brightness was modestly anti-correlated with the jovian magnetic latitude, decreasing when Io intersected the torus centrifugal equator.     

The Dual Sources of Io's Sodium Clouds

Io's sodium clouds result mostly from a combination of two atmospheric escape processes at Io. Neutralization of Na⁺ and/or NaX⁺ pickup ions produces the “stream” and the “jet” and results in a rectangular-shaped sodium nebula around Jupiter. Atmospheric sputtering of Na by plasma torus ions produces the “banana cloud” near Io and a diamond-shaped sodium nebula. Charge exchange of thermal Na⁺ with Na in Io's atmosphere does not appear to be a major atmospheric ejection process. The total ejection rate of sodium from Io varied from 3×10²⁶ to 25×10²⁶ atoms/s over seven years of observations. Our results provide further evidence that Io's atmospheric escape is driven from collisionally thick regions of the atmosphere rather than from the exosphere.     

The post-eclipse brightening of Io

We review the photometric work on eclipse reappearances of Io. New observations of eclipse reappearances of Io confirm the post-eclipse brightness anomaly reported by Binder and Cruikshank (1964) but testify to its intermittent nature. A post-eclipse anomaly of approximately 0.07 mag was observed on two occasions in 1972, while observations of Europa and Ganymede showed no brightness anomaly greater than 0.01 mag. The atmospheric condensation model for the anomaly on Io is reviewed in terms of the quantity of frost required to produce the effect and the corresponding amount of gas liberated to the atmosphere upon sublimation. The observational data and the results from a stellar occultation are in general accord with the theoretical predictions of the stability of heavy gases on Io, while both observational and theoretical criteria are satisfied by a tenuous atmosphere of a heavy gas such as methane or ammonia having a surface pressure ～10^?7 bar.     

The ion mass loading rate at Io

The Io plasma torus, composed of mostly heavy ions of oxygen and sulfur, is sustained by an Iogenic mass loading rate of ∼10³⁰ amu s⁻¹ = 1.6 × 10²⁸ SO₂ s⁻¹ or approximately 10³ kg s⁻¹(A.L. Broadfoot et al., 1979, Science 204, 979-982). We argue on the basis of available power sources, reanalysis of F. Bagenal (1997, Geophys. Res. Lett. 24, 2111-2114), HST UV remote sensing, and detailed model calculations that at most 20% of this mass leaves Io in the form of ions, i.e., ≤3 × 10²⁷ × (n_e,0/3600 cm⁻³) ions s⁻¹, where n_e,0 is the average torus electron density. For the Galileo spacecraft Io pass in December 1995, the ion mass loading rate was ≤3 × 10²⁷ ions s⁻¹, whereas for the Voyager epoch with lower n_e,0 (=2000 cm⁻³), this rate would be ≤1.7 × 10²⁷ ions s⁻¹, consistent with the D.E. Shemansky (1980, Astrophys. J. 242, 1266-1277) mass loading limit of ≤1 × 10²⁷ ions s⁻¹. We investigate the processes that control Io’s large scale electrodynamic interaction and find that the elastic collision rate exceeds the ionization/pickup rate by at least a factor of 5 for all atmospheric column densities considered (10¹⁶-10²¹ m⁻²) and by a factor of ∼100 for the most realistic column density. Consequently, elastic collisions are mostly responsible for Io’s high conductances and thus generate Io’s large scale electrodynamic interaction such as the generation of Io’s electric current system and the slowing of the plasma flow. The electrodynamic part of Io’s interaction is thus best described as an ionosphere-like interaction rather than a comet-like interaction. An analytic expression for total electron impact rates is derived for Io’s atmosphere, which is independent of any particular model for the 3D interaction of torus electrons with its atmosphere.     

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.     

The polar contribution to the heat flow of Io

Polar brightness temperatures on Io are higher than expected for any passive surface heated by absorbed sunlight. This discrepancy implies large scale volcanic activity from which we derive a new component of Io's heat flow. We present a ‘Three Component’ thermal background, infrared emission model for Io that includes active polar regions. The widespread polar activity contributes an additional ∼0.6 W m⁻² to Io's heat flow budget above the ∼2.5 W m⁻² from thermal anomalies. Our estimate for Io's global average heat flow increases to ∼3±1 W m⁻² and ∼1.3±0.4×10¹⁴ watts total.     

A cometary ionosphere model for Io

The ionosphere of Jupiter's satellite Io, discovered by the Pioneer 10 radio-occultation experiment, cannot easily be understood in terms of a model of a gravitationally bound, Earth-like ionosphere. Io's gravitational field is so weak that a gravitationally bound ionosphere would probably be blown away by the ram force of the Jovian magnetospheric wind — i.e., the plasma corotating in the Jovian magnetosphere. We propose here a model in which the material for Io's atmosphere and ionosphere is drawn from the ionosphere of Jupiter through a Birkeland current system that is driven by the potential induced across Io by the Jovian corotation electric field. We argue that the ionization near Io is caused by a comet-like interaction between the corotating plasma and Io's atmosphere. The initial interaction employs the critical velocity phenomenon proposed many years ago by Alfvén. Further ionization is produced by the impact of Jovian trapped energetic electrons, and the ionization thus created is swept out ahead of Io in its orbit. Thus, we suggest that what has been reported as a day-night ionospheric asymmetry is in fact an upstream-downstream asymmetry caused by the Jovian magnetospheric wind.Paper dedicated to Professor Hannes Alfvén on the occasion of his 70th birthday, 30th May, 1978.     

Cassini/VIMS observation of an Io post-eclipse brightening event

During the Cassini-Jupiter flyby, VIMS observed Io at different phase angles, both in full sunlight and in eclipse. By using the sunlight measurements, we were able to produce phase curves in the visual through all the near infrared wavelengths covered by the VIMS instrument (0.85-5.1 μm). The phase angle spanned from ∼2° to ∼120°. The measurements, done just after Io emerged from Jupiter's shadow, show an increase of about 15% in Io's reflectance with respect to what would be predicted by the phase curve. This behavior is observed at wavelengths >1.2 μm. Moreover, just after emergence from eclipse an increase of about 25% is observed in the depth of SO₂ frost bands at 4.07 and 4.35 μm. At 0.879<λ<1.04 μm the brightening is 10-24%. Below λ=0.879 μm the brightening, if present, should be less than the precision of our measurements (∼5%). Apparently, these observations are not explained neither by a diverse spatial distribution of SO₂ on the Io' surface nor by atmospheric SO₂ condensation on the surface during the eclipse.     

UV spectrum of Enceladus

We present a far ultraviolet (FUV) spectrum of Saturn’s moon Enceladus from the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST). We have put upper limits on emission from C, N, and O lines in Enceladus’ atmosphere and column densities for the C lines assuming solar resonance scattering. We find these upper limits to be relatively low—on the order of tens to thousands of Rayleighs and with C column densities on the order of 10⁸–10¹⁵ cm^?2, depending on the assumed source size. We also present a segment of a reflectance spectrum in the FUV from ～1900–2130 Å. This region was sensitive to the different ice mixtures in the model spectra reported by Hendrix et al. (Hendrix, A.R., Hansen, C.J., Holsclaw, G.M. [2010]. Icarus, 206, 608). We find the spectrum brightens quickly longward of ～1900 Å, constraining the absorption band observed by Hendrix et al. from ～170 to 190 nm. We find our data is consistent with the suggestion of Hendrix et al. of the presence of ammonia ice (or ammonia hydrate) to darken that region, and also possibly tholins to darken the mid-UV, as reported by Verbiscer et al. (Verbiscer, A.J., French, R.G., McGhee, C.A. [2005]. Icarus, 173, 66).     

Io's dayside SO2 atmosphere

An extensive set of HI Lyman-α images obtained with the Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) from 1997-2001 has been analyzed to provide information about the spatial and temporal character of Io's SO₂ atmosphere. An atmospheric distribution map derived from the observations reveals that the sunlit SO₂ atmosphere is temporally stable on a global scale, with only small local changes. An anti-/sub-jovian asymmetry in the SO₂ distribution is present in all 5 years of the observations. The average daytime atmosphere is densest on the anti-jovian hemisphere in the equatorial regions, with a maximum equatorial column density of 5.0×¹⁰¹⁶ cm−2 at 140° longitude. The SO₂ atmosphere also has greater latitudinal extent on the anti-jovian hemisphere as compared to the sub-jovian. The atmospheric distribution appears to be best correlated with the location of hot spots and known volcanic plumes, although small number statistics for the plumes limits the correlation.     

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

Io’s atmosphere: Constraints on sublimation support from density variations on seasonal timescales using NASA IRTF/TEXES observations from 2001 to 2010

We present an analysis of 19 μm spectra of Io’s SO₂ atmosphere from the TEXES mid-infrared high spectral resolution spectrograph on NASA’s Infrared Telescope Facility, incorporating new data taken between January 2005 and June 2010 and a re-analysis of earlier data taken from November 2001 to January 2004. This is the longest set of contiguous observations of Io’s atmosphere using the same instrument and technique thus far. We have fitted all 16 detected blended absorption lines of the ν₂ SO₂ vibrational band to retrieve the subsolar values of SO₂ column abundance and the gas kinetic temperature. By incorporating an existing model of Io’s surface temperatures and atmosphere, we retrieve sub-solar column densities from the disk-integrated data. Spectra from all years are best fit by atmospheric temperatures <150 K. Best-fit gas kinetic temperatures on the anti-Jupiter hemisphere, where SO₂ gas abundance is highest, are low and stable, with a mean of 108 (±18) K. The sub-solar SO₂ column density between longitudes of 90–220° varies from a low of 0.61 (±0.145) × 10^?17 cm^?2, near aphelion in 2004, to a high of 1.51 (±0.215) × 10¹⁷ cm^?2 in 2010 when Jupiter was approaching its early 2011 perihelion. No correlation in the gas temperature was seen with the increasing SO₂ column densities outside the errors.Assuming that any volcanic component of the atmosphere is constant with time, the correlation of increasing SO₂ abundance with decreasing heliocentric distance provides good evidence that the atmosphere is at least partially supported by frost sublimation. The SO₂ frost thermal inertias and albedos that fit the variation in atmospheric density best are between 150–1250 W m^?2 s^?1/2 K^?1 and 0.613–0.425 respectively. Photometric evidence favors albedos near the upper end of this range, corresponding to thermal inertias near the lower end. This relatively low frost thermal inertia produces larger amplitude seasonal variations than are observed, which in turn implies a substantial additional volcanic atmospheric component to moderate the amplitude of the seasonal variations of the total atmosphere on the anti-Jupiter hemisphere. The seasonal thermal inertia we measure is unique both because it refers exclusively to the SO₂ frost surface component, and also because it refers to relatively deep subsurface layers (few meters) due to the timescales of many years, while previous studies have determined thermal inertias at shallower levels (few centimeters), relevant for timescales of ～2 h (eclipse) or ～2 days (diurnal curves).     

Galileo observations of volcanic plumes on Io

Io's volcanic plumes erupt in a dazzling variety of sizes, shapes, colors and opacities. In general, the plumes fall into two classes, representing distinct source gas temperatures. Most of the Galileo imaging observations were of the smaller, more numerous Prometheus-type plumes that are produced when hot flows of silicate lava impinge on volatile surface ices of SO₂. Few detections were made of the giant, Pele-type plumes that vent high temperature, sulfur-rich gases from the interior of Io; this was partly because of the insensitivity of Galileo's camera to ultraviolet wavelengths. Both gas and dust spout from plumes of each class. Favorably located gas plumes were detected during eclipse, when Io was in Jupiter's shadow. Dense dust columns were imaged in daylight above several Prometheus-type eruptions, reaching heights typically less than 100 km. Comparisons between eclipse observations, sunlit images, and the record of surface changes show that these optically thick dust columns are much smaller in stature than the corresponding gas plumes but are adequate to produce the observed surface deposits. Mie scattering calculations suggest that these conspicuous dust plumes are made up of coarse grained “ash” particles with radii on the order of 100 nm, and total masses on the order of 10⁶ kg per plume. Long exposure images of Thor in sunlight show a faint outer envelope apparently populated by particles small enough to be carried along with the gas flow, perhaps formed by condensation of sulfurous “snowflakes” as suggested by the plasma instrumentation aboard Galileo as it flew through Thor's plume [Frank, L.A., Paterson, W.R., 2002. J. Geophys. Res. (Space Phys.) 107, doi:10.1029/2002JA009240. 31-1]. If so, the total mass of these fine, nearly invisible particles may be comparable to the mass of the gas, and could account for much of Io's rapid resurfacing.     

Photochemistry and Vertical Transport in Io's Atmosphere and Ionosphere

We present an updated model for the photochemistry of Io's atmosphere and ionosphere and use this model to investigate the sensitivity of the chemical structure to vertical transport rates. SO₂is assumed to be the dominant atmospheric gas, with minor molecular sodium species such as Na₂S or Na₂O released by sputtering or venting from the surface. Photochemical products include SO, O₂, S, O, Na, NaO, NaS, and Na₂. We consider both “thick” and “thin” SO₂atmospheres that encompass the range allowed by recent HST and millimeter-wave observations, and evaluate the possibility that O₂and/or SO may be significant minor dayside constituents and therefore likely dominant nightside gases. The fast reaction between S and O₂limits the column abundance of O₂to ∼10⁴less than that calculated by Kumar (J. Geophys. Res.87, 1677–1684, 1982; 89(A9), 7399–7406, 1984) for a pure sulfur/oxygen atmosphere. If a significant source of NaO₂or Na₂O were supplied by the surface and mixed rapidly upward, then oxygen liberated in the chemical reactions which also liberate free Na would provide an additional source of O₂. Fast eddy mixing will enhance the transport of molecular sodium species to the exobase, in addition to increasing the vertical transport rate of ions. Ions produced in the atmosphere will be accelerated by the reduced corotation electric field penetrating the atmosphere. These ions experience collisions with the neutral gas, leading to enhanced vertical ion diffusion. The dominant ion, Na⁺, is lost primarily by charge exchange with Na₂O and/or Na₂S in the lower atmosphere and by diffusion through the ionopause in the upper atmosphere. The atmospheric column abundance of SO, O₂, and the upper atmosphere escape rates of Na, S, O, and molecular sodium species are all strong functions of the eddy mixing rate. Most atmospheric escape, including that of molecular sodium species, probably occurs from the low density “background” SO₂atmosphere, while a localized high density “volcanic” SO₂atmosphere can yield an ionosphere consistent with that detected by the Pioneer 10 spacecraft.     

The Atmosphere of Io: Abundances and Sources of Sulfur Dioxide and Atomic Hydrogen

An analysis and interpretation of reflected solar Lyman α intensity data acquired with the Hubble Space Telescope (HST) implies an equatorially confined atmosphere with SO₂ column densities ∼ 1–2 × 10¹⁶ cm^-2. Poleward of 30° the SO₂ density must decrease sharply reaching an asymptotic polar value of < 10¹⁵ cm^-2 at 45° to achieve the observed 2 kR intensity peaks. The corresponding surface reflectivities must be either a constant 0.047 for higher equatorial SO₂ or a variable reflectivity of 0.027 with lower SO₂ densities at the equator increasing to a polar value of ∼ 0.05. The average residence time for an atmospheric SO₂ molecule is ∼ 2–3 days for the canonical mass loading rate of the Io plasma torus = 10³⁰ amu s^-1. With atomic hydrogen in the atmosphere and corona constrained by the HST observations, it is estimated that a pickup proton density ratio of 0.25–0.4% can be sustained by a supply of Io plasma torus protons neutralized in Io's atmosphere/exosphere, if protons constitute 7% of the total torus ion density, which is close to the Chust et al. (1999) pickup proton density ratio and under the widely quoted 10% proton content of the torus. This revised version was published online in July 2006 with corrections to the Cover Date.     

Energy deposition and primary chemical products in Titan’s upper atmosphere

Cassini results indicate that solar photons dominate energy deposition in Titan’s upper atmosphere. These dissociate and ionize nitrogen and methane and drive the subsequent complex organic chemistry. The improved constraints on the atmospheric composition from Cassini measurements demand greater precision in the photochemical modeling. Therefore, in order to quantify the role of solar radiation in the primary chemical production, we have performed detailed calculations for the energy deposition of photons and photoelectrons in the atmosphere of Titan and we validate our results with the Cassini measurements for the electron fluxes and the EUV/FUV emissions. We use high-resolution cross sections for the neutral photodissociation of N₂, which we present here, and show that they provide a different picture of energy deposition compared to results based on low-resolution cross sections. Furthermore, we introduce a simple model for the energy degradation of photoelectrons based on the local deposition approximation and show that our results are in agreement with detailed calculations including transport, in the altitude region below 1200 km, where the effects of transport are negligible. Our calculated, daytime, electron fluxes are in good agreement with the measured fluxes by the Cassini Plasma Spectrometer (CAPS), and the same holds for the measured FUV emissions by the Ultraviolet Imaging Spectrometer (UVIS). Finally, we present the vertical production profiles of radicals and ions originating from the interaction of photons and electrons with the main components of Titan’s atmosphere, along with the column integrated production rates at different solar zenith angles. These can be used as basis for any further photochemical calculations.