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.     

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.     

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.     

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.     

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.     

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.     

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

Cassini-VIMS at Jupiter: solar occultation measurements using Io

We report unusual and somewhat unexpected observations of the jovian satellite Io, showing strong methane absorption bands. These observations were made by the Cassini VIMS experiment during the Jupiter flyby of December/January 2000/2001. The explanation is straightforward: Entering or exiting from Jupiter's shadow during an eclipse, Io is illuminated by solar light which has transited the atmosphere of Jupiter. This light, therefore becomes imprinted with the spectral signature of Jupiter's upper atmosphere, which includes strong atmospheric methane absorption bands. Intercepting solar light refracted by the jovian atmosphere, Io essentially becomes a “mirror” for solar occultation events of Jupiter. The thickness of the layer where refracted solar light is observed is so large (more than 3000 km at Io's orbit), that we can foresee a nearly continuous multi-year period of similar events at Saturn, utilizing the large and bright ring system. During Cassini's 4-year nominal mission, this probing technique should reveal information of Saturn's atmosphere over a large range of southern latitudes and times.     

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.     

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.     

Io's volcanic control of Jupiter's extended neutral clouds

Dramatic changes in the brightness and shape of Jupiter's extended sodium nebula are found to be correlated with the infrared emission brightness of Io. Previous imaging and modeling studies have shown that varying appearances of the nebula correspond to changes in the rate and the type of loss mechanism for atmospheric escape from Io. Similarly, previous IR observational studies have assumed that enhancements in infrared emissions from Io correspond to increased levels of volcanic (lava flow) activity. In linking these processes observationally and statistically, we conclude that silicate volcanism on Io controls both the rate and the means by which sodium escapes from Io's atmosphere. During active periods, molecules containing sodium become an important transient in Io's upper atmosphere, and subsequent photochemistry and molecular-ion driven dynamics enhance the high speed sodium population, leading to the brightest nebulas observed. This is not the case during volcanically quiet times when omni-present atmospheric sputtering ejects sodium to form a modest, base-level nebula. Sodium's role as a “trace gas” of the more abundant species of sulfur (S) and oxygen (O) is less certain during volcanic episodes. While we suggest that volcanism must also affect the escape rates of S and O, and consequently their extended neutral clouds, the different roles played by lava and plume sources for non-sodium species are far too uncertain to make definitive comparisons at this time.     

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.     

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.     

Observations of Electromagnetically Coupled Dust in the Jovian Magnetosphere

We report on dust measurements obtained during the seventh orbit of the Galileo spacecraft about Jupiter. The most prominent features observed are highly time variable dust streams recorded throughout the Jovian system. The impact rate varied by more than an order of magnitude with a 5 and 10 hour periodicity, which shows a correlation with Galileo's position relative to the Jovian magnetic field. This behavior can be qualitatively explained by strong coupling of nanometer-sized dust to the Jovian magnetic field. In addition to the 5 and 10 h periodicities, a longer period which is compatible with Io's orbital period is evident in the dust impact rate. This feature indicates that Io most likely is the source of the dust streams. During a close (3,095 km altitude) flyby at Ganymede on 5 April 1997 an enhanced rate of dust impacts has been observed, which suggests that Ganymede is a source of ejecta particles. Within a distance of about 25 R_J(Jupiter radius, R_J= 71,492 km) from Jupiter impacts of micrometer-sized particles have been recorded which could be particles on bound orbits about Jupiter. This revised version was published online in July 2006 with corrections to the Cover Date.     

Dynamic of Polar Plumes Observed During 2006, 2008, 2009 and 2010 Total Solar Eclipses

Based on observations of polar plume kinematics in the white-light corona during the total solar eclipse in 2006, the images obtained during multi-station observation of the eclipses of 2006, 2008, 2009 and 2010 were analysed. Several polar plumes showing similar kinematics were identified. The speeds of these dynamic features were found by comparing images obtained at different times along the path of totality. A possible connection with erupting spicules and macrospicules is discussed.     

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.     

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.     

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

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