Morphology, temperature, and eruption dynamics at Pele

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

Sulfur volcanism on Io

In February 2003, March 2003 and January 2004 Pele plume transmission spectra were obtained during Jupiter transit with Hubble's Space Telescope Imaging Spectrograph (STIS), using the 0.1″ wide slit and the G230LB grating. The STIS spectra covered the 2100-3100 Å wavelength regions and extended spatially along Io's limb encompassing the region directly above and northward of the vent of the Pele volcano. The S₂ and SO₂ absorption signatures evident in these data indicate that the gas signature at Pele was temporally variable, and that an S₂ absorption signature was present ∼12° from the Pele vent near 6±5 S and 264±15 W, suggesting the presence of another S₂ bearing plume on Io. Contemporaneous with the spectral data, UV and visible-wavelength images of the plume were obtained in reflected sunlight with the Advanced Camera for Surveys (ACS) prior to Jupiter transit. The dust scattering recorded in these data provide an additional qualitative measure of plume activity on Io, indicating that the degree of dust scattering over Pele varied as a function of the date of observation, and that there were several other dust bearing plumes active during the observations. We present constraints on the composition and variability of the gas abundances of the Pele plume as well as the plumes detected by ACS and recorded within the STIS data, as a function of time.     

Numerical modeling of ionian volcanic plumes with entrained particulates

Volcanic plumes on Jupiter's moon Io are modeled using the direct simulation Monte Carlo (DSMC) method. The modeled volcanic vent is interpreted as a “virtual” vent. A parametric study of the “virtual” vent gas temperature and velocity is performed to constrain the gas properties at the vent by observables, particularly the plume height and the surrounding condensate deposition ring radius. Also, the flow of refractory nano-size particulates entrained in the gas is modeled with “overlay” techniques which assume that the background gas flow is not altered by the particulates. The column density along the tangential line-of-sight and the shadow cast by the plume are calculated and compared with Voyager and Galileo images. The parametric study indicates that it is possible to obtain a unique solution for the vent temperature and velocity for a large plume like Pele. However, for a small Prometheus-type plume, several different possible combinations of vent temperature and velocity result in both the same shock height and peak deposition ring radius. Pele and Prometheus plume particulates are examined in detail. Encouraging matches with observations are obtained for each plume by varying both the gas and particle parameters. The calculated tangential gas column density of Pele agrees with that obtained from HST observations. An upper limit on the size of particles that track the gas flow well is found to be ∼10 nm, consistent with Voyager observations of Loki. While it is certainly possible for the plumes to contain refractory dust or pyroclastic particles, especially in the vent vicinity, we find that the conditions are favorable for SO₂ condensation into particles away from the vent vicinity for Prometheus. The shadow cast by Prometheus as seen in Galileo images is also reproduced by our simulation. A time averaged frost deposition profile is calculated for Prometheus in an effort to explain the multiple ring structure observed around the source region. However, this multiple ring structure may be better explained by the calculated deposition of entrained particles. The possibility of forming a dust cloud on Io is examined and, based on a lack of any such observed clouds, a subsolar frost temperature of less than 118 K is suggested.     

Characterizing Io’s Pele,Tvashtar and Pillan plumes: Lessons learned from Hubble

Hubble Space Telescope/Wide Field and Planetary Camera 2 (HST/WFPC2) images of Io obtained between 1995 and 2007 between 0.24 and 0.42 μm led to the detection of the Pele plume in reflected sunlight in 1995 and 1999; imaging of the Pele plume via absorption of jovian light in 1996 and 1999; detection of the Prometheus-type Pillan plume in reflected sunlight in 1997; and detection of the 2007 Pele-type Tvashtar plume eruption in reflected sunlight and via absorption of jovian light. Based on a detailed analysis of these observations we characterize and compare the gas and dust properties of each of the detected plumes. In each case, the brightness of the plumes in reflected sunlight is less at 0.26 μm than at 0.33 μm. Mie scattering analysis of the wavelength dependence of each plume’s reflectance signature suggests that range of particle sizes within the plumes is quite narrow. Assuming a normal distribution of particle sizes, the range of mean particle sizes is ～0.035–0.12 μm for the 1997 Pillan eruption, ～0.05–0.08 μm for the 1999 Pele and 2007 Tvasthar plumes, and ～0.05–0.11 μm for the 1995 Pele plume, and in each case the standard deviation in the particle size distribution is <15%. The Mie analysis also suggests that the 2007 Tvashtar eruption released ～10⁹ g of sulfur dust, the 1999 Pele eruption released ～10⁹ g of SO₂ dust, the 1997 Pillan eruption released ～10¹⁰ g of SO₂ dust, and the 1995 Pele plume may have released ～10¹⁰ g of SO₂ dust. Analysis of the plume absorption signatures recorded in the F255W filter bandpass (0.24–0.28 μm) indicates that the opacity of the 2007 Tvashtar plume was 2× that of the 1996 and 1999 Pele plume eruptions. While the sulfur dust density estimated for the Tvashtar from the reflected sunlight data could have produced 61% of the observed plume opacity, <10% of the 1999 Pele F255W plume opacity could have resulted from the SO₂ dust detected in the eruption. Accounting for the remaining F255W opacity level of the Pele and Tvasthar plumes based on SO₂ and S₂ gas absorption, the SO₂ and S₂ gas density inferred for each plume is almost equivalent corresponding to ～2–6 × 10¹⁶ cm^?2 and 3–5 × 10¹⁵ cm^?2, respectively, producing SO₂ and S₂ gas resurfacing rates ～0.04–0.2 cm yr^?1 and 0.007–0.01 cm yr^?1; and SO₂ and S₂ gas masses ～1–4 × 10¹⁰ g and ～2–3 × 10⁹ g; for a total dust to gas ratio in the plumes ～10^?1–10^?2. The 2007 Tvashtar plume was detected by HST at ～380 ± 40 km in both reflected sunlight and absorbed jovian light; in 1999, the detected Pele plume altitude was 500 km in absorbed jovian light, but in reflected sunlight the detected height was ～2× lower. Thus, for the 1999 Pele plume, similar to the 1979 Voyager Pele plume observations, the most efficient dust reflections occurred in the region closest to the plume vent. The 0.33–0.42 μm brightness of the 1997 Pillan plume was 10–20× greater than the Pele or Tvashtar plumes, exceeding by a factor of 3 the average brightness levels observed within 200 km of 1979 Loki eruption vent. But, the 0.26 μm brightness of the 1997 Pillan plume in reflected sunlight was significantly lower than would be predicted by the dust scattering model. Presuming that the 0.26 μm brightness of the 1997 Pillan plume was attenuated by the eruption plume’s gas component, then an SO₂ gas density ～3–6 × 10¹⁸ cm^?2 is inferred from the data (for S₂/SO₂ ratios ?4%), comparable to the 0.3–2 × 10¹⁸ cm^?2 SO₂ density detected at Loki in 1979 (Pearl, J.C. et al. [1979]. Nature 280, 755; Lellouch et al., 1992), and producing an SO₂ gas mass ～3–8 × 10¹¹ g and an SO₂ resurfacing rate ～8–23 cm yr^?1. These results confirm the connection between high (?10¹⁷ cm^?2) SO₂ gas content and plumes that scatter strongly at nearly blue wavelengths, and it validates the occurrence of high density SO₂ gas eruptions on Io. Noting that the SO₂ gas content inferred from a spectrum of the 2003 Pillan plume was significantly lower ～2 × 10¹⁶ cm^?2 (Jessup, K.L., Spencer, J., Yelle, R. [2007]. Icarus 192, 24–40); and that the Pillan caldera was flooded with fresh SO₂ frost/slush just prior to the 1997 Pillan plume eruption (Geissler, P., McEwen, A., Phillips, C., Keszthelyi, L., Spencer, J. [2004a]. Icarus 169, 29–64; Phillips, C.B. [2000]. Voyager and Galileo SSI Views of Volcanic Resurfacing on Io and the Search for Geologic Activity at Europa. Ph.D. Thesis, Univ. of Ariz., Tucson); we propose that the density of SO₂ gas released by this volcano is directly linked to the local SO₂ frost abundance at the time of eruption.     

Charting thermal emission variability at Pele,Janus Patera and Kanehekili Fluctus with the Galileo NIMS Io Thermal Emission Database (NITED)

Using the NIMS Io Thermal Emission Database (NITED), a collection of over 1000 measurements of radiant flux from Io’s volcanoes (Davies, A.G. et al. [2012]. Geophys. Res. Lett. 39, L01201. doi:10.1029/2011GL049999), we have examined the variability of thermal emission from three of Io’s volcanoes: Pele, Janus Patera and Kanehekili Fluctus. At Pele, the 5-μm thermal emission as derived from 28 night time observations is remarkably steady at 37 ± 10 GW μm^?1, re-affirming previous analyses that suggested that Pele an active, rapidly overturning silicate lava lake. Janus Patera also exhibits relatively steady 5-μm thermal emission (≈20 ± 3 GW μm^?1) in the four observations where Janus is resolved from nearby Kanehekili Fluctus. Janus Patera might contain a Pele-like lava lake with an effusion rate (Q_F) of ≈40–70 m³ s^?1. It should be a prime target for a future mission to Io in order to obtain data to determine lava eruption temperature. Kanehekili Fluctus has a thermal emission spectrum that is indicative of the emplacement of lava flows with insulated crusts. Effusion rate at Kanehekili Fluctus dropped by an order of magnitude from ≈95 m³ s^?1 in mid-1997 to ≈4 m³ s^?1 in late 2001.     

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.     

Sulfur flows of Ra Patera,Io

Voyager 1 imaging data have been used to investigate the color and morphology of several radial flow-like features at Ra Patera, a broad volcanic structure at approximately 8° latitude and 325° longitude on the Galilean satellite Io (J1). It was found that downstream progressions of flow color and morphology are consistent with lava of a predominately sulfur composition cooling radiatively and erupting in the range of 470 to 520°K at effusion rates at 10¹⁰ to 10¹¹ cm³/sec. This implies global resurfacing rates by volcanic flows on Io of the order of 1 cm/year. Calculated energy content and effusion rates for flows at Ra Patera, using the physical parameters of sulfur, are of the order of the largest known terrestial basaltic eruptions and are consistent with calculations of globally available energy.     

Supersonic meteorology of Io: Sublimation-driven flow of SO2

The horizontal flow of SO₂ gas from day side to night side of Io is calculated. The surface is assumed to be covered by a frost whose vapor pressure at the subsolar point is orders of magnitude larger than that on the night side. Temperature of the frost is controlled by radiation. The flow is hydrostatic and turbulent, with velocity and entropy per particle independent of height. The vertically integrated conservation equations for mass, momentum, and energy are solved for atmospheric pressure, temperature, and horizontal velocity as functions of solar zenith angle. Formulas from boundary layer theory govern the interaction between atmosphere and surface. The flow becomes supersonic as it expands away from the subsolar point, as in the theory of rocket nozzles and the solar wind. Within 35° of the subsolar point atmospheric pressureis less than the frost vapor pressure, and the frost sublimes. Elsewhere, atmospheric pressure is greater than the frost vapor pressure, and the frost condenses. The two pressures seldom differ by more than a factor of 2. The sublimation rate at the subsolar point is proportional to the frost vapor pressure, which is a sensitive function of temperature. For a subsolar temperature of 130°K, the sublimation rate is 10¹⁵ molecules/cm²/sec. Diurnally averaged sublimation rates at the equator are comparable to the 0.1 cm/year resurfacing rate required for burial of impact craters. At the poles where both the vapor pressures and atmospheric pressures are low, the condensation rates are 100 times smaller. Surface pressures near the terminator are generally too low to account for the ionosphere discovered by Pioneer 10. The possibility of a noncondensable gas in addition to SO₂ must be seriously considered.     

Simulation of gas dynamics and radiation in volcanic plumes on Io

Modeling results of volcanic plumes on Jupiter’s moon Io are presented. Two types of low density axisymmetric SO₂ plume flows are modeled using the direct simulation Monte Carlo (DSMC) method. Thermal radiation from all three vibrational bands and overall rotational lines of SO₂ molecules is modeled. A high resolution computation of the flow in the vicinity of the vent was obtained by multidomain sequential calculation to improve the modeling of the radiation signature. The radiation features are examined both by calculating infrared emission spectra along different lines-of-sight through the plume and with the DSMC modeled emission images of the whole flow field. It is found that most of the radiation originates in the vicinity of the vent, and non-LTE (non-local-thermodynamic equilibrium) cooling by SO₂ rotation lines exceeds cooling in the v₂ vibrational band at high altitude.In addition to the general shape of the plumes, the calculated average SO₂ column density (∼10¹⁶ cm⁻²) over a Pele-type plume and the related frost-deposition ring structure (at R ∼ 500 km from the vent) are in agreement with observations. These comparisons partially validate the modeling. It is suggested that an observation with spatial resolution of less than 30 km is needed to measure the large spatial variation of SO₂ near a Pele-type plume center. It is also found that an influx of 1.1 × 10²⁹ SO₂ s⁻¹ (or 1.1 × 10⁴ kg s⁻¹) is sufficient to reproduce the observed SO₂ column density at Pele. The simulation results also show some interesting features such as a multiple bounce shock structure around Prometheus-type plumes and the frost depletion by plume-induced erosion on the sunlit side of Io. The model predicts the existence of a canopy shock, a ballistic region inside the Pele-type plume, and the negligible effect of surface heating by plume emission.     

Speckle imaging of volcanic hotspots on Io with the Keck telescope

Using speckle imaging techniques on the 10-m W.M. Keck I telescope, we observed near-infrared emission at 2.2 μm from volcanic hotspots on Io in July-August 1998. Using several hundreds of short-exposure images we reconstructed diffraction-limited images of Io on each of three nights. We measured the positions of individual hotspots to ±0.004″ or better, corresponding to a relative positional error of ∼20 km on Io's surface. The sensitivity of normal ground-based images of Io is limited by confusion between overlapping sources; by resolving these multiple points we detected up to 17 distinct hotspots, the largest number ever seen in a single image.During the month-long span of our 1998 observations, several events occurred. Loki was at the end of a long brightening, and we observed it to fade in flux by a factor of 2.8 over the course of one month. At the 3-sigma level we see evidence that Loki's position shifts by ∼100 km. This suggests that the brightening may not have been located at the “primary” Loki emission center but at a different source within the Loki caldera. We also see a bright transient source near Loki. Among many other sources we detect a dim source on the limb of Io at the latitude of Pele; this source is consistent with 2.7% of the thermal emission from the Pele volcano complex being scattered by the Pele plume, which would be the first detection of a plume through scattered infrared hotspot emission.     

Shield volcano topography and the rheology of lava flows on Io

Stereo and photoclinometry derived topography of shield-like volcanoes on Io indicate little relief (<3 km) and very low slopes (0.2° to 0.6°). Several shield volcanoes appear to be associated with broad rises of 1 to 3 km, but only 5 shield volcanoes have been identified with steep flank slopes (between 4° and 10°). These steep slopes are restricted to within 20-30 km of the summit, but where discernable, most of the lava flows observed on these edifices occur on the outer flanks where slopes are less than a degree. Despite their abundance, ionian shield volcanoes are among the flattest in the Solar System. The steepest volcanoes on Io are most comparable to large venusian shield volcanoes. Using simplistic Bingham rheologies we estimate the viscosity and yield strengths of ionian lavas. Yield strengths are estimated at 10¹-10² Pa, lower than most basaltic lavas. Viscosity estimates range from 10³ to 10⁵ Pa s, although these are probably upper limits. Actual values may have been as low as 10⁰ Pa s. Viscosity is sensitive to flow velocity, which is poorly known on Io. The best constraint on flow velocity comes from observations of the 1997 Pillan eruption, which bracket the eruptive phase to 132 day maximum, and more probably less than 50 days. Low slopes, long run-out distances and our estimated rheologic properties are consistent with (but not proof of) a low silica, low viscosity, high temperature composition for ionian lavas, supporting arguments for low-silica lava compositions such as basalt or komatiite. We cannot eliminate sulfur on rheologic grounds, however.     

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.     

Transfer of mass from Io to Europa and beyond due to cometary impacts

We simulate the production and orbital evolution of escaping ejecta due to cometary impacts on Io. The model includes the four Galilean satellites, Amalthea, Thebe, Jupiter's gravitational moments, Saturn and the Sun. Five scenarios are examined: an impact at the apex, the sub-jovian point, the anti-jovian point, the antapex, and at the south pole of Io. We estimate that on average a cometary impact injects thrice its mass (in the form of Io surface material) into jovicentric orbit. The majority of the escaping debris comes back to Io, but a sizeable fraction (between 5.0 and 8.7%) manages to reach Europa, and a smaller fraction Ganymede (between 1.5 and 4.6%). Smaller fractions reached Amalthea Thebe, Callisto, and Jupiter itself. For million year time scales, the mass transfer to Europa is estimated as 1.8-3.1×¹⁰¹⁴ g/Myr. The median time for transfer of ejecta from Io to Europa is ∼56 years.     

Sputtering of sulfur by kiloelectronvolt ions: Application to the magnetospheric plasma interaction with Io

Measurements of total yields, temperature dependences, mass spectra, and energy spectra of molecules sputtered from condensed sulfur (S₈) at low temperatures by keV ions are reported and results are given for Jovian plasma ion bombardment of Io. A change in the reflectance of the sulfur, which can be removed by annealing, is produced by the most penetrating ions and may be connected with the darker, colder polar regions on Io. The measured sputtering yields are much lower than those estimated earlier for room temperature sulfur films but are comparable to previous measurements of keV ion sputtering of SO₂ at low temperatures. The corrected mass spectrum indicates that ≈66% of the total yield corresponds to S₂ ejection while only 5 and 16% correspond to S and S₃, respectively. Therefore, if ions reach the surface of Io its atmosphere will have a non-negligible sulfur component of primarily S₂. The ejection of S and S₂ is temperature independent for temperatures characteristic of most of the surface of Io. The energy spectrum for S has an approximate 1/E² dependence at high ejection energies, whereas S₂ and S₃ fall off more rapidly. Assuming 50% coverage of both sulfur and SO₂ and a thin atmosphere (e.g., nightside and polar region) the direct sputter injection of sulfur atoms and molecules into the Jovian plasma torus and the indirect injection due to coronal processes are estimated. These injection rates for sulfur are compared to those for SO₂ showing that injection from sulfur deposits contributes 13% to the total mass injection rate of ∼2–3 × 10²⁹ amu/sec.     

Alkali and Chlorine Photochemistry in a Volcanically Driven Atmosphere on Io

Observations of the Io plasma torus and neutral clouds indicate that the extended ionian atmosphere must contain sodium, potassium, and chlorine in atomic and/or molecular form. Models that consider sublimation of pure sulfur dioxide frost as the sole mechanism for generating an atmosphere on Io cannot explain the presence of alkali and halogen species in the atmosphere—active volcanoes or surface sputtering must also be considered, or the alkali and halide species must be discharged along with the SO₂ as the frost sublimates. To determine how volcanic outgassing can affect the chemistry of Io's atmosphere, we have developed a one-dimensional photochemical model in which active volcanoes release a rich suite of S-, O-, Na-, K-, and Cl-bearing vapor and in which photolysis, chemical reactions, condensation, and vertical eddy and molecular diffusion affect the subsequent evolution of the volcanic gases. Observations of Pele plume constituents, along with thermochemical equilibrium calculations of the composition of volcanic gases exsolved from high-temperature silicate magmas on Io, are used to constrain the composition of the volcanic vapor. We find that NaCl, Na, Cl, KCl, and K will be the dominant alkali and chlorine gases in atmospheres generated from Pele-like plume eruptions on Io. Although the relative abundances of these species will depend on uncertain model parameters and initial conditions, these five species remain dominant for a wide variety of realistic conditions. Other sodium and chlorine molecules such as NaS, NaO, Na₂, NaS₂, NaO₂, NaOS, NaSO₂, SCl, ClO, Cl₂, S₂Cl, and SO₂Cl₂ will be only minor constituents in the ionian atmosphere because of their low volcanic emission rates and their efficient photochemical destruction mechanisms. Our modeling has implications for the general appearance, properties, and variability of the neutral sodium clouds and jets observed near Io. The neutral NaCl molecules present at high altitudes in atmosph eres generated by active volcanoes might provide the NaX⁺ ion needed to help explain the morphology of the high-velocity sodium “stream” feature observed near Io.     

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.     

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.     

Mass-injection rate from Io into the Io plasma torus

Theoretical arguments have been presented to the effect that both plasma and energy are supplied to the Jovian magnetosphere primarily from internal sources. If we assume that Io is the source of plasma for the Jovian magnetosphere and that outward flow of plasma from the torus is the means of drawing from the kinetic energy of rotation of Jupiter to drive magnetospheric phenomena, we can obtain a new, independent estimate of the rate of mass injection from Io into the Io plasma torus. We explicitly assume the solar wind supplies neither plasma nor energy to the Jovian magnetosphere in significant amounts. The power expended by the Jovian magnetosphere is supplied by torus plasma falling outward through the corotational-centrifugal-potential field. A lower limit to the rate of mass injection into the torus, which on the average must equal the rate of mass loss from the torus, is therefore derivable if we adopt a value for the power expended to drive the various magnetospheric phenomena. This method yields an injection rate of at least 10³ kg/sec, a value in agreement with the results obtained by two other independent methods of estimating mass injection rate. If this injection rate from Io and extraction of energy from Jupiter's kinetic energy of rotation has been maintained over geologic time, then approximately 0.1% of Io's mass (principally in the form of sulfur and oxygen) has been lost to the Jovian magnetosphere, and Jupiter's spin rate has been reduced by less than 0.1%.     

The sodium and hydrogen gas clouds of Io

Models are developed to describe the spatial distribution of gases emitted by Io and are applied to recent observations which indicate extensive gas clouds of hydrogen and sodium in orbit around Jupiter. Hydrogen and sodium atoms are emitted from Io with velocities in the range 2 to 3 km sec^?1, with fluxes of about 10¹⁰ and 10⁸cm^?2sec^?1 for hydrogen and sodium respectively. Hydrogen atoms may be formed by photodecomposition of gases such as NH₃ or H₂S released from the satellite surface and may escape thermally from an exosphere whose temperature is about 500 K. Sodium may be ejected from the surface by energetic particles or by ultraviolet radiation and it appears that a non-thermal mechanism drawing energy from Jupiter's magnetic field is required in order to account for its release to space.     

Mapping of Io's thermal radiation by the Galileo photopolarimeter-radiometer (PPR) instrument

Between 1999 and 2002, the Galileo spacecraft made 6 close flybys of Io during which many observations of Io's thermal radiation were made with the photopolarimeter-radiometer (PPR). While the NIMS instrument could measure thermal emission from hot spots with T>200 K, PPR was the only Galileo instrument capable of mapping the lower temperatures of older, cooling lava flows, and the passive background. We tabulate all data taken by PPR of Io during these flybys and describe some scientific highlights revealed by the data. The data include almost complete coverage of Io at better than 250 km resolution, with extensive regional coverage at higher resolutions. We found a modest poleward drop in nighttime background temperatures and evidence of thermal inertia variations across the surface. Comparison of high spatial resolution temperature measurements with observed daytime SO₂ gas pressures on Io provides evidence for local cold trapping of SO₂ frost on scales smaller than the 60 km resolution of the PPR data. We also calculated the power output from several hot spots and estimated total global heat flow to be about 2.0-2.6 W m⁻². The low-latitude diurnal temperature variations for the regions between obvious hot spots are well matched by a laterally-inhomogeneous thermal model with less than 1 W m⁻² endogenic heat flow.