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.     

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.     

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.     

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.     

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.     

Near-surface flow of volcanic gases on Io

Significant near-surface flow of gas several hundred kilometers from Pele (Plume 1) on Io is indicated by a series of bright, elongate albedo markings. Particles produced at small, local vents are apparently carried as much as 70 km farther “downwind” from Pele. The gas densities and velocities necessary to suspend 0.1 to 10 μm particles at such a distance imply mass flow rates of 10⁷ to 10⁹ g/sec. Such flow rates are consistent with other estimates of mass transport by the plume. The large flow rates so far from the source allow an estimate of the rate of resurfacing of Io by lava flows and pyroclastics that is independent of estimates based on meteorite flux or on the amount of solids carried within the plumes themselves.     

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.     

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.     

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.     

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.     

Photochemistry of a Volcanically Driven Atmosphere on Io: Sulfur and Oxygen Species from a Pele-Type Eruption

To determine how active volcanism might affect the standard picture of sulfur dioxide photochemistry on Io, we have developed a one-dimensional atmospheric model in which a variety of sulfur-, oxygen-, sodium-, potassium-, and chlorine-bearing volatiles are volcanically outgassed at Io's surface and then evolve due to photolysis, chemical kinetics, and diffusion. Thermochemical equilibrium calculations in combination with recent observations of gases in the Pele plume are used to help constrain the composition and physical properties of the exsolved volcanic vapors. Both thermochemical equilibrium calculations (Zolotov and Fegley 1999, Icarus141, 40-52) and the Pele plume observations of Spencer et al. (2000; Science288, 1208-1210) suggest that S₂ may be a common gas emitted in volcanic eruptions on Io. If so, our photochemical models indicate that the composition of Io's atmosphere could differ significantly from the case of an atmosphere in equilibrium with SO₂ frost. The major differences as they relate to oxygen and sulfur species are an increased abundance of S, S₂, S₃, S₄, SO, and S₂O and a decreased abundance of O and O₂ in the Pele-type volcanic models as compared with frost sublimation models. The high observed SO/SO₂ ratio on Io might reflect the importance of a contribution from volcanic SO rather than indicate low eddy diffusion coefficients in Io's atmosphere or low SO “sticking” probabilities at Io's surface; in that case, the SO/SO₂ ratio could be temporally and/or spatially variable as volcanic activity fluctuates. Many of the interesting volcanic species (e.g., S₂, S₃, S₄, and S₂O) are short lived and will be rapidly destroyed once the volcanic plumes shut off; condensation of these species near the source vent is also likely. The diffuse red deposits associated with active volcanic centers on Io may be caused by S₄ radicals that are created and temporarily preserved when sulfur vapor (predominantly S₂) condenses around the volcanic vent. Condensation of SO across the surface and, in particular, in the polar regions might also affect the surface spectral properties. We predict that the S/O ratio in the torus and neutral clouds might be correlated with volcanic activity—during periods when volcanic outgassing of S₂ (or other molecular sulfur vapors) is prevalent, we would expect the escape of sulfur to be enhanced relative to that of oxygen, and the S/O ratio in the torus and neutral clouds could be correspondingly increased.     

The final Galileo SSI observations of Io: orbits G28-I33

We present the observations of Io acquired by the Solid State Imaging (SSI) experiment during the Galileo Millennium Mission (GMM) and the strategy we used to plan the exploration of Io. Despite Galileo's tight restrictions on data volume and downlink capability and several spacecraft and camera anomalies due to the intense radiation close to Jupiter, there were many successful SSI observations during GMM. Four giant, high-latitude plumes, including the largest plume ever observed on Io, were documented over a period of eight months; only faint evidence of such plumes had been seen since the Voyager 2 encounter, despite monitoring by Galileo during the previous five years. Moreover, the source of one of the plumes was Tvashtar Catena, demonstrating that a single site can exhibit remarkably diverse eruption styles—from a curtain of lava fountains, to extensive surface flows, and finally a ∼400 km high plume—over a relatively short period of time (∼13 months between orbits I25 and G29). Despite this substantial activity, no evidence of any truly new volcanic center was seen during the six years of Galileo observations. The recent observations also revealed details of mass wasting processes acting on Io. Slumping and landsliding dominate and occur in close proximity to each other, demonstrating spatial variation in material properties over distances of several kilometers. However, despite the ubiquitous evidence for mass wasting, the rate of volcanic resurfacing seems to dominate; the floors of paterae in proximity to mountains are generally free of debris. Finally, the highest resolution observations obtained during Galileo's final encounters with Io provided further evidence for a wide diversity of surface processes at work on Io.     

Surface changes on Io during the Galileo mission

A careful survey of Galileo SSI global monitoring images revealed more than 80 apparent surface changes that took place on Io during the 5 year period of observation, ranging from giant plume deposits to subtle changes in the color or albedo of patera surfaces. Explosive volcanic activity was discovered at four previously unrecognized centers: an unnamed patera to the south of Karei that produced a Pele-sized red ring, a patera to the west of Zal that produced a small circular bright deposit, a large orange ring detected near the north pole of Io, and a small bright ring near Io's south pole. Only a handful of Io's many active volcanoes produced large scale explosive eruptions, and several of these erupted repeatedly, leaving at least 83% of Io's surface unaltered throughout the Galileo mission. Most of the hot spots detected from SSI, NIMS and ground-based thermal observations caused no noticeable surface changes greater than 10 km in extent over the five year period. Surface changes were found at every location where active plumes were identified, including Acala which was never seen in sunlight and was only detected through auroral emissions during eclipse. Two types of plumes are distinguished on the basis of the size and color of their deposits, confirming post-Voyager suggestions by McEwen and Soderblom [Icarus 55 (1983) 191]. Smaller plumes produce near-circular rings typically 150-200 km in radius that are white or yellow in color unless contaminated with silicates, and frequently coat their surroundings with frosts of fine-grained SO₂. The larger plumes are much less numerous, limited to a half dozen examples, and produce oval, orange or red, sulfur-rich rings with maximum radii in the north-south direction that are typically in the range from 500 to 550 km. Both types of plumes can be either episodic or quasi-continuous over a five year period. Repeated eruptions of the smaller SO₂-rich plumes likely contribute significantly to Io's resurfacing rate, whereas dust ejection is likely dominated by the tenuous giant plumes. Both types of plume deposits fade on time-scales of months to years through burial and alteration. Episodic seepages of SO₂ at Haemus Montes, Zal Montes, Dorian Montes, and the plateau to the north of Pillan Patera may have been triggered by activity at nearby volcanic centers.     

Two classes of volcanic plumes on Io

Comparison of Voyager 1 and Voyager 2 images of the south polar region of Io has revealed that a major volcanic eruption occured there during the period between the two spacecraft encounters. An annular deposit ～1400 km in diameter formed around the Aten Patera caldera (311°W, 48°S), the floor of which changed from orange to red-black. The characteristics of this eruption are remarkably similar to those described earlier for an eruption centered on Surt caldera (338°W, 45°N) that occured during the same period, also at high latitude, but in the north. Both volcanic centers were evidently inactive during the Voyager 1 and 2 encounters but were active sometime between the two. The geometric and colorimetric characteristics, as well as scale of the two annular deposits, are virtually identical; both resemble the surface features formed by the eruption of Pele (255°W, 18°S). These three very large plume eruptions suggest a class of eruption distinct from that of six smaller plumes observed to be continously active by both Voyagers 1 and 2. The smaller plumes, of which Prometheus is the type example, are longer-lived, deposit bright, whitish material, erupt at velocities of ～0.5 km sec^?1, and are concentrated at low latitudes in an equatorial belt around the satellite. The very large Pele-type plumes, on the other hand, are relatively short-lived, deposit darker red materials, erupt at ～1.0 km sec^?1, and (rather than restricted to a latitudinal band) are restricted in longitude from 240° to 360°W. Both direct thermal infrared temperature measurements and the implied color temperatures for quenched liquid sulfur suggest that hot spot temperatures of ～650°K are associated with the large plumes and temperatures <400°K with the small plumes. The typical eruption duration of the small plumes is at least several years; that of the large plumes appears to be of the order of days to weeks. The two classes therefore differ by more than two orders of magnitude in duration of eruption. Based on uv, visible, and infrared spectra, the small plumes seem to contain and deposit SO₂ in their annuli whereas the large plumes apparently do not. Two other plumes that occur at either end of the linear feature Loki may be intermediate or hybrid between the two classes, exhibiting attributes of both. Additionally, Loki occurs in the area of overlap in the regional distributions of the two plume classes. Two distinct volcanic systems involving different volatiles may be responsible for the two classes. We propose that the discrete temperatures associated with the two classes are a direct reflection of sulfur's peculiar variation in viscosity with temperature. Over two temperature ranges (～400 to 430°K and >650°K), sulfur is a low-viscosity fluid (orange and black, respectively); at other temperatures it is either solid or has a high viscosity. As a result, there will be two zones in Io's crust in which liquid sulfur will flow freely: a shallow zone of orange sulfur and a deeper zone of black sulfur. A low-temperature system driven by SO₂ heated to 400 to 400°K by the orange sulfur zone seems the best model for the small plumes; a system driven by sulfur heated to >650°K by hot or even molten silicates in the black sulfur zone seems the best explanation for the large plume class. The large Pele-type plumes are apparently concentrated in a region of the satellite in which a thinner sulfur-rich crust overlies the tidally heated silicate lithosphere, so the black sulfur zone may be fairly shallow in this region. The Prometheus-type plumes are possibly confined to the equatorial belt by some process that concentrates SO₂ fluid in the equatorial crust.     

Evolution of Io's volatile inventory

Using Voyager results, we have made crude estimates of the rate at which Io loses volatiles by a variety of processes to the surrounding magnetosphere for both the current SO₂-dominated atmosphere as well as hypothetical paleoatmospheres in which other gases, such as N₂, may have been the dominant constituent. Loss rates are strongly influenced by the surface pressure on the night side, the relationship between the exobase and the Jovian magnetospheric boundary, the exospheric temperature, and the peak altitudes reached by volcanic plumes. Several mechanisms make significant contributions to the prodigious rate at which Io is currently losing volatiles. These include: interaction of the magnetospheric plasma with volcanic plume particles and the background atmosphere; sputtering of ices on the surface, if the nightside atmospheric pressure is low enough; and Jeans' escape of O, a dissociation product of SO₂ gas. For paleoatmospheres, only the first two of these mechanisms would have been effective. However, they are capable of eliminating large amounts of N₂ and other volatiles from Io over the satellite's lifetime. Io could have also lost large amounts of water over its lifetime due to the extensive recycling of water between its upper and lower crust, with the partial dissociation of water vapor in silicate magma chambers initiating this loss process. Significant amounts of water may also have been lost as a result of the interaction of the magnetospheric plasma with water ice particles in volcanic plumes. Once an SO₂-dominated atmosphere becomes established, much water may have also been lost through the sputtering of surface water ice.     

The size distribution of Jupiter's main ring from Galileo imaging and spectroscopy

Galileo's Solid State Imaging experiment (SSI) obtained 36 visible wavelength images of Jupiter's ring system during the nominal mission (Ockert-Bell et al., 1999, Icarus 138, 188-213) and another 21 during the extended mission. The Near Infrared Mapping Spectrometer (NIMS) recorded an observation of Jupiter's main ring during orbit C3 at wavelengths from 0.7 to 5.2 μm; a second observation was attempted during orbit E4. We analyze the high phase angle NIMS and SSI observations to constrain the size distribution of the main ring's micron-sized dust population. This portion of the population is best constrained at high phase angles, as the light scattering behavior of small dust grains dominates at these geometries and contributions from larger ring particles are negligible. High phase angle images of the main ring obtained by the Voyager spacecraft covered phase angles between 173.8° and 176.9° (Showalter et al., 1987, Icarus 69, 458-498). Galileo images extend this range up to 178.6°. We model the Galileo phase curve and the ring spectra from the C3 NIMS ring observation as the combination of two power law distributions. Our analysis of the main ring phase curve and the NIMS spectra suggests the size distribution of the smallest ring particles is a power law with an index of 2.0±0.3 below a size of ∼15 μm that transitions to a power law with an index of 5.0±1.5 at larger sizes. This combined power law distribution, or “broken power law” distribution, yields a better fit to the NIMS data than do the power law distributions that have previously been fit to the Voyager imaging data (Showalter et al., 1987, Icarus 69, 458-498). The broken power law distribution reconciles the results of Showalter et al. (1987, Icarus 69, 458-498) and McMuldroch et al. (2000, Icarus 146, 1-11), who also analyzed the NIMS data, and can be considered as an obvious extension of a simple power law. This more complex size distribution could indicate that ring particle production rates and/or lifetimes vary with size and may relate to the physical processes that control their evolution. The significant near arm/far arm asymmetry reported elsewhere (see Showalter et al., 1987, Icarus 69, 458-498; Ockert-Bell et al., 1999, Icarus 138, 188-213) persists in the data even after the main ring is isolated in the SSI images. However, the sense of the asymmetry seen in Galileo images differs from that seen in Voyager images. We interpret this asymmetry as a broad-scale, azimuthal brightness variation. No consistent association with the magnetic field of Jupiter has been observed. It is possible that these longitudinal variations may be similar to the random brightness fluctuations observed in Saturn's F ring by Voyager (Smith et al., 1982, Science 215, 504-537) and during the 1995 ring plane crossings (Nicholson et al., 1996, Science 272, 509-515; Bosh and Rivkin, 1996, Science 272, 518-521; Poulet et al., 2000, Icarus 144, 135-148). Stochastic events may thus play a significant role in the evolution of the jovian main ring.     

Geysers of Enceladus: Quantitative analysis of qualitative models

Aspects of two qualitative models of Enceladus’ dust plume—the so-called “Cold Faithful” [Porco, C.C., et al., 2006. Cassini observes the active south pole of Enceladus. Science 311, 1393-1401; Ingersoll, A.P., et al., 2006. Models of the Enceladus plumes. In: Bulletin of the American Astronomical Society, vol. 38, p. 508] and “Frigid Faithful” [Kieffer, S.W., et al., 2006. A clathrate reservoir hypothesis for Enceladus’ south polar plume. Science 314, 1764; Gioia, G., et al., 2007. Unified model of tectonics and heat transport in a Frigid Enceladus. Proc. Natl. Acad. Sci. 104, 13578-13591] models—are analyzed quantitatively. The former model assumes an explosive boiling of subsurface liquid water, when pressure exerted by the ice crust is suddenly released due to an opening crack. In the latter model the existence of a deep shell of clathrates below Enceladus’ south pole is conjectured; clathrates can decompose explosively when exposed to vacuum through a fracture in the outer icy shell. For the Cold Faithful model we estimate the maximal velocity of ice grains, originating from water splashing in explosive boiling. We find that for water near the triple point this velocity is far too small to explain the observed plume properties. For the Frigid Faithful model we consider the problem of momentum transfer from gas to ice particles. It arises since any change in the direction of the gas flow in the cracks of the shell requires re-acceleration of the entrained grains. While this effect may explain the observed speed difference of gas and grains if the gas evaporates from triple point temperature (273.15 K) [Schmidt, J., et al., 2008. Formation of Enceladus dust plume. Nature 451, 685], the low temperatures of the Frigid Faithful model imply a too dilute vapor to support the observed high particle fluxes in Enceladus’ plume.     

HiRISE observations of gas sublimation-driven activity in Mars’ southern polar regions: IV. Fluid dynamics models of CO2 jets

The High Resolution Imaging Science Experiment (HiRISE) onboard Mars Reconnaissance Orbiter (MRO) has been used to monitor the seasonal evolution of several regions at high southern latitudes on Mars and, in particular, the jet-like activity which may result from the process described by Kieffer (Kieffer, H.H. [2007]. J. Geophys. Res. (Planets) 112, E08005. doi:10.1029/2006JE002816) involving translucent CO₂ ice. In this work, we concentrate on attempting to model the dusty CO₂ gas jets using a computational fluid dynamics code. Models that included surface slopes of up to 20° (as an analogy to the jet activity seen in “Inca City”, 81°S, 296°E), wind (from 0 to 6 m s⁻¹), variable vent cross-section and length, particles (including a particle size distribution) and mass loading (with dust to gas ratios exceeding 1) were investigated. The structure of the resulting gas jets, the particle distribution within the jets, the deposition patterns (including their dependence on particle size), and the appearance of jets when viewed from different orientations (including from a nadir-pointing camera) have been investigated for a range of input parameters. The results provide predictions for the size-dependency of altitudes of particles within a plume and the distribution of particle sizes in the deposition fans. Where slopes are a dominant influence, larger particles are expected to be seen furthest from the vent. Where wind is dominant, smaller particles should travel to larger distances. Models producing deposition patterns consistent in length (∼80 m) and form with fans observed by HiRISE on MRO have been demonstrated. The models also suggest that downward flow of gas produced by drag effects from particles falling from the jet under gravity could provide a mechanism for the production of bright haloes which are observed to surround dark fan deposits in MOC, HiRISE and CRISM.     

Filamentation of volcanic plumes on the Jovian satellite Io

Volcanic plumes on the Jovian satellite Io may be a visible manifestation of a plasma-arc discharge phenomenon. The amount of power in the plasma arc (10¹¹ W) is not enough to account for all the energy dissipated by the volcanoes. However, once a volcano is initiated by tidal and geologic processes, the dynamics of the volcanic plumes can be influenced by the plasma arcs. As initially pointed out by Gold (1979), plasma arcs are expected because of 10⁶ A currents and 400 kV potentials generated by the flow past Io of a torus of relatively dense magnetospheric plasma. We utilize our experience with laboratory plasma arcs to investigate the plume dynamics. The filamentation in the plume of the volcano Prometheus and its cross-sectional shape is quantitatively consistent with theories developed from laboratory observation.Paper dedicated to Professor Hannes Alfvén on the occasion of his 80th birthday, 30 May 1988.     

Properties and dynamics of Jupiter's gossamer rings from Galileo, Voyager, Hubble and Keck images

We present a comprehensive examination of Jupiter's “gossamer” rings based on images from Voyager, Galileo, the Hubble Space Telescope and the W.M. Keck Telescope. We compare our results to the simple dynamical model of Burns et al. [Burns, J.A., Showalter, M.R., Hamilton, D.P., Nicholson, P.D., de Pater, I., Ockert-Bell, M., Thomas, P., 1999. Science 284, 1146-1150] in which dust is ejected from Amalthea and Thebe and then evolves inward under Poynting-Robertson drag. The ring follows many predictions of the model rather well, including a linear reduction in thickness with decreasing radius. However, some deviations from the model are noted. For example, additional material appears to be concentrated just interior to the orbits of the two moons. At least in the case of Amalthea's ring, that material is in the same orbital plane as Amalthea's inclined orbit and may be trapped at the Lagrange points. Thebe's ring shows much larger vertical excursions from the model, which may be related to perturbations by several strong Lorentz resonances. Photometry is consistent with the dust obeying a relatively flat power-law size distribution, very similar to dust in the main ring. However, the very low backscatter reflectivity of the ring, and the flat phase curve of the ring at low phase angles, require that the ring be composed of distinctly non-spherical particles.