Geologic mapping of the Zal region of Io

We produced a regional geologic map of the Zal region of Io's antijovian hemisphere using Galileo mission data. We discuss the geologic features, summarize the map units and structures that are present, discuss the nature of volcanic activity, and present an analysis of the volcanic, tectonic, and gradational processes that affect the region. The Zal region consists of five primary types of geologic materials: plains, mountains, paterae floors, flows, and diffuse deposits. The flows and patera floors are similar, but are subdivided based on uncertainties regarding emplacement environments and mechanisms. The Zal region includes two hotspots detected by Galileo: one along the western scarp of the Zal Patera volcano and one at the Rustam Patera volcano (name submitted to IAU). A third hotspot at the nearby At'am Patera volcano (name submitted to IAU) is the source of diffuse and pyroclastic materials that blanket north Zal Mons. The western bounding scarp of Zal Patera is the location of a fissure vent that is the source of multiple silicate lava flows. The floor of Zal Patera has been partially resurfaced by dark lava flows, although portions of the patera floor appear bright and unchanged during the Galileo mission. This suggests that the floor did not undergo complete resurfacing as a flooding lava lake but does contain a compound flow field. Mountain materials exhibit stages of degradation; lineated material degrades into mottled material. We have explored the possibility that north and south Zal Mons were originally one structure. We propose that strike-slip faulting and subsequent rifting separated the mountain units, opened a fissure which serves as a vent for lava flow, and created a depression which, by further extension during the rifting event, became Zal Patera. With comparison to other regional maps of Io, this work provides insight into the general geologic evolution of Io.     

Geologic mapping of the Amirani-Gish Bar region of Io: Implications for the global geologic mapping of Io

We produced the first geologic map of the Amirani-Gish Bar region of Io, the last of four regional maps generated from Galileo mission data. The Amirani-Gish Bar region has five primary types of geologic materials: plains, mountains, patera floors, flows, and diffuse deposits. The flows and patera floors are thought to be compositionally similar, but are subdivided based on interpretations regarding their emplacement environments and mechanisms. Our mapping shows that volcanic activity in the Amirani-Gish Bar region is dominated by the Amirani Eruptive Center (AEC), now recognized to be part of an extensive, combined Amirani-Maui flow field. A mappable flow connects Amirani and Maui, suggesting that Maui is fed from Amirani, such that the post-Voyager designation “Maui Eruptive Center” should be revised. Amirani contains at least four hot spots detected by Galileo, and is the source of widespread bright (sulfur?) flows and active dark (silicate?) flows being emplaced in the Promethean style (slowly emplaced, compound flow fields). The floor of Gish Bar Patera has been partially resurfaced by dark lava flows, although other parts of its floor are bright and appeared unchanged during the Galileo mission. This suggests that the floor did not undergo complete resurfacing as a lava lake as proposed for other ionian paterae. There are several other hot spots in the region that are the sources of both active dark flows (confined within paterae), and SO₂- and S₂-rich diffuse deposits. Mapped diffuse deposits around fractures on mountains and in the plains appear to serve as the source for gas venting without the release of magma, an association previously unrecognized in this region. The six mountains mapped in this region exhibit various states of degradation. In addition to gaining insight into this region of Io, all four maps are studied to assess the best methodology to use to produce a new global geologic map of Io based on the newly released, combined Galileo-Voyager global mosaics. To convey the complexity of ionian surface geology, we find that a new global geologic map of Io should include a map sheet displaying the global abundances and types of surface features as well as a complementary GIS database as a means to catalog the record of surface changes observed since the Voyager flybys and during the Galileo mission.     

Volcanism on Io: New insights from global geologic mapping

We produced the first complete, 1:15 M-scale global geologic map of Jupiter’s moon Io, based on a set of monochrome and color Galileo-Voyager image mosaics produced at a spatial resolution of 1 km/pixel. The surface of Io was mapped into 19 units based on albedo, color and surface morphology, and is subdivided as follows: plains (65.8% of surface), lava flow fields (28.5%), mountains (3.2%), and patera floors (2.5%). Diffuse deposits (DD) that mantle the other units cover ∼18% of Io’s surface, and are distributed as follows: red (8.6% of surface), white (6.9%), yellow (2.1%), black (0.6%), and green (∼0.01%). Analyses of the geographical and areal distribution of these units yield a number of results, summarized below. (1) The distribution of plains units of different colors is generally geographically constrained: Red-brown plains occur >±30° latitude, and are thought to result from enhanced alteration of other units induced by radiation coming in from the poles. White plains (possibly dominated by SO₂ + contaminants) occur mostly in the equatorial antijovian region (±30°, 90-230°W), possibly indicative of a regional cold trap. Outliers of white, yellow, and red-brown plains in other regions may result from long-term accumulation of white, yellow, and red diffuse deposits, respectively. (2) Bright (possibly sulfur-rich) flow fields make up 30% more lava flow fields than dark (presumably silicate) flows (56.5% vs. 43.5%), and only 18% of bright flow fields occur within 10 km of dark flow fields. These results suggest that secondary sulfurous volcanism (where a bright-dark association is expected) could be responsible for only a fraction of Io’s recent bright flows, and that primary sulfur-rich effusions could be an important component of Io’s recent volcanism. An unusual concentration of bright flows at ∼45-75°N, ∼60-120°W could be indicative of more extensive primary sulfurous volcanism in the recent past. However, it remains unclear whether most bright flows are bright because they are sulfur flows, or because they are cold silicate flows covered in sulfur-rich particles from plume fallout. (3) We mapped 425 paterae (volcano-tectonic depressions), up from 417 previously identified by Radebaugh et al. (Radebaugh, J., Keszthelyi, L.P., McEwen, A.S., Turtle, E.P., Jaeger, W., Milazzo, M. [2001]. J. Geophys. Res. 106, 33005-33020). Although these features cover only 2.5% of Io’s surface, they correspond to 64% of all detected hot spots; 45% of all hot spots are associated with the freshest dark patera floors, reflecting the importance of active silicate volcanism to Io’s heat flow. (4) Mountains cover only ∼3% of the surface, although the transition from mountains to plains is gradational with the available imagery. 49% of all mountains are lineated and presumably layered, showing evidence of linear structures supportive of a tectonic origin. In contrast, only 6% of visible mountains are mottled (showing hummocks indicative of mass wasting) and 4% are tholi (domes or shields), consistent with a volcanic origin. (5) Initial analyses of the geographic distributions of map units show no significant longitudinal variation in the quantity of Io’s mountains or paterae, in contrast to earlier studies. This is because we use the area of mountain and patera materials as opposed to the number of structures, and our result suggests that the previously proposed anti-correlation of mountains and paterae (Schenk, P., Hargitai, H., Wilson, R., McEwen, A., Thomas, P. [2001]. J. Geophys. Res. 106, 33201-33222; Kirchoff, M.R., McKinnon, W.B., Schenk, P.M. [2011]. Earth Planet. Sci. Lett. 301, 22-30) is more complex than previously thought. There is also a slight decrease in surface area of lava flows toward the poles of Io, perhaps indicative of variations in volcanic activity. (6) The freshest bright and dark flows make up about 29% of all of Io’s flow fields, suggesting active emplacement is occurring in less than a third of Io’s visible lava fields. (7) About 47% of Io’s diffuse deposits (by area) are red, presumably deriving their color from condensed sulfur gas, and ∼38% are white, presumably dominated by condensed SO₂. The much greater areal extent of gas-derived diffuse deposits (red + white, 85%) compared to presumably pyroclast-bearing diffuse deposits (dark (silicate tephra) + yellow (sulfur-rich tephra), 15%) indicates that there is effective separation between the transport of tephra and gas in many Ionian explosive eruptions. Future improvements in the geologic mapping of Io can be obtained via (a) investigating the relationships between different color/material units that are geographically and temporally associated, (b) better analysis of the temporal variations in the map units, and (c) additional high-resolution images (spatial resolutions ∼200 m/pixel or better). These improvements would be greatly facilitated by new data, which could be obtained by future missions.     

Io: Heat flow from dark paterae

Dark paterae on the jovian satellite Io are evidence of recent volcanic activity. Some paterae appear to be entirely filled with dark volcanic material, while others have only partially darkened floors. Dark paterae have area and heat flow longitudinal distributions that are bimodal as well as anti-correlated with the longitudinal distribution of mountains on Io at a global scale. As part of our study of Io’s total heat flow, we have examined the darkest paterae and quantified their thermal emission in order to assess their contribution. This is the first time that the areas of the dark material in these paterae have been measured with such precision and correlated with their thermal emission. Dark paterae yield a significantly larger contribution to Io’s heat flow than dark volcanic fields. Dark paterae (including Loki Patera) yield at least ∼4 × 10¹³ W or ∼40% of Io’s total heat flow. In comparison, dark flow fields yield ∼10¹³ W or ∼10% of Io’s total heat flow. Of the total heat loss from dark paterae, Loki Patera alone yields ∼10¹³ W or ∼10% of Io’s total thermal 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.     

The Zamama-Thor region of Io: Insights from a synthesis of mapping, topography, and Galileo spacecraft data

We have studied data from the Galileo spacecraft's three remote sensing instruments (Solid-State Imager (SSI), Near-Infrared Mapping Spectrometer (NIMS), and Photopolarimeter-Radiometer (PPR)) covering the Zamama-Thor region of Io's antijovian hemisphere, and produced a geomorphological map of this region. This is the third of three regional maps we are producing from the Galileo spacecraft data. Our goal is to assess the variety of volcanic and tectonic materials and their interrelationships on Io using planetary mapping techniques, supplemented with all available Galileo remote sensing data. Based on the Galileo data analysis and our mapping, we have determined that the most recent geologic activity in the Zamama-Thor region has been dominated by two sites of large-scale volcanic surface changes. The Zamama Eruptive Center is a site of both explosive and effusive eruptions, which emanate from two relatively steep edifices (Zamama Tholi A and B) that appear to be built by both silicate and sulfur volcanism. A ∼100-km long flow field formed sometime after the 1979 Voyager flybys, which appears to be a site of promethean-style compound flows, flow-front SO₂ plumes, and adjacent sulfur flows. Larger, possibly stealthy, plumes have on at least one occasion during the Galileo mission tapped a source that probably includes S and/or Cl to produce a red pyroclastic deposit from the same vent from which silicate lavas were erupted. The Thor Eruptive Center, which may have been active prior to Voyager, became active again during the Galileo mission between May and August 2001. A pillanian-style eruption at Thor included the tallest plume observed to date on Io (at least 500 km high) and new dark lava flows. The plume produced a central dark pyroclastic deposit (probably silicate-rich) and an outlying white diffuse ring that is SO₂-rich. Mapping shows that several of the new dark lava flows around the plume vent have reoccupied sites of earlier flows. Unlike most of the other pillanian eruptions observed during the Galileo mission, the 2001 Thor eruption did not produce a large red ring deposit, indicating a relative lack of S and/or Cl gases interacting with the magma during that eruption. Between these two eruptive centers are two paterae, Thomagata and Reshef. Thomagata Patera is located on a large shield-like mesa and shows no signs of activity. In contrast, Reshef Patera is located on a large, irregular mesa that is apparently undergoing degradation through erosion (perhaps from SO₂-sapping or chemical decomposition of sulfur-rich material) from multiple secondary volcanic centers.     

Post-solidification cooling and the age of Io's lava flows

The modeling of thermal emission from active lava flows must account for the cooling of the lava after solidification. Models of lava cooling applied to data collected by the Galileo spacecraft have, until now, not taken this into consideration. This is a flaw as lava flows on Io are thought to be relatively thin with a range in thickness from ∼1 to 13 m. Once a flow is completely solidified (a rapid process on a geological time scale), the surface cools faster than the surface of a partially molten flow. Cooling via the base of the lava flow is also important and accelerates the solidification of the flow compared to the rate for the ‘semi-infinite’ approximation (which is only valid for very deep lava bodies). We introduce a new model which incorporates the solidification and basal cooling features. This model gives a superior reproduction of the cooling of the 1997 Pillan lava flows on Io observed by the Galileo spacecraft. We also use the new model to determine what observations are necessary to constrain lava emplacement style at Loki Patera. Flows exhibit different cooling profiles from that expected from a lava lake. We model cooling with a finite-element code and make quantitative predictions for the behavior of lava flows and other lava bodies that can be tested against observations both on Io and Earth. For example, a 10-m-thick ultramafic flow, like those emplaced at Pillan Patera in 1997, solidifies in ∼450 days (at which point the surface temperature has cooled to ∼280 K) and takes another 390 days to cool to 249 K. Observations over a sufficient period of time reveal divergent cooling trends for different lava bodies [examples: lava flows and lava lakes have different cooling trends after the flow has solidified (flows cool faster)]. Thin flows solidify and cool faster than flows of greater thickness. The model can therefore be used as a diagnostic tool for constraining possible emplacement mechanisms and compositions of bodies of lava in remote-sensing data.     

Io: Heat flow from dark volcanic fields

Dark flow fields on the jovian satellite Io are evidence of current or recent volcanic activity. We have examined the darkest volcanic fields and quantified their thermal emission in order to assess their contribution to Io’s total heat flow. Loki Patera, the largest single source of heat flow on Io, is a convenient point of reference. We find that dark volcanic fields are more common in the hemisphere opposite Loki Patera and this large scale concentration is manifested as a maximum in the longitudinal distribution (near ∼200 °W), consistent with USGS global geologic mapping results. In spite of their relatively cool temperatures, dark volcanic fields contribute almost as much to Io’s heat flow as Loki Patera itself because of their larger areal extent. As a group, dark volcanic fields provide an asymmetric component of ∼5% of Io’s global heat flow or ∼5 × 10¹² W.     

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.     

The variability of volcanic activity at Zamama, Culann, and Tupan Patera on Io as seen by the Galileo Near Infrared Mapping Spectrometer

Zamama, Culann, and Tupan Patera are three large, persistent volcanic centers on the jovian moon Io. As part of an ongoing project to quantify contributions from individual volcanic centers to Io’s thermal budget, we have quantified the radiant flux from all suitable observations made by the Galileo Near Infrared Mapping Spectrometer (NIMS) of these volcanoes, in some cases filling omissions in previous analyses. At Zamama, after a long period of cooling, we see a peak in thermal emission that corresponds with new plume activity. Subsequently, toward the end of the Galileo epoch, thermal emission from Zamama drops off in a manner consistent with a greatly reduced eruption rate and the cooling of emplaced flows. Culann exhibits possible episodic activity. We present the full Tupan Patera NIMS dataset and derive new estimates of thermal output and temporal behavior. Eruption rates at these three volcanoes are on the order of 30 m³ s⁻¹, consistent with a previous analysis of NIMS observations of Prometheus, and nearly an order of magnitude greater than Kilauea volcano, Hawai’i, Earth’s most active volcano. We propose that future missions to the jovian system could better constrain activity at these volcanoes and others where similar styles of activity are taking place by obtaining data on a time scale of, ideally, at least one observation per day. Observations at similar or even shorter timescales are desirable during initial waxing phases of eruption episodes. These eruptions are identifiable from their characteristic spectral signatures and temporal behavior.     

Volcanic activity at Tvashtar Catena, Io

Galileo's Solid State Imager (SSI) observed Tvashtar Catena four times between November 1999 and October 2001, providing a unique look at a distinctive high latitude volcanic complex on Io. The first observation (orbit I25, November 1999) resolved, for the first time, an active extraterrestrial fissure eruption; the brightness temperature was at least 1300 K. The second observation (orbit I27, February 2000) showed a large (∼500 km²) region with many, small, hot, regions of active lava. The third observation was taken in conjunction with Cassini imaging in December 2000 and showed a Pele-like, annular plume deposit. The Cassini images revealed an ∼400 km high Pele-type plume above Tvashtar Catena. The final Galileo SSI observation of Tvashtar (orbit I32, October 2001), revealed that obvious (to SSI) activity had ceased, although data from Galileo's Near Infrared Mapping Spectrometer (NIMS) indicated that there was still significant thermal emission from the Tvashtar region. In this paper, we primarily analyze the style of eruption during orbit I27 (February 2000). Comparison with a lava flow cooling model indicates that the behavior of the Tvashtar eruption during I27 does not match that of simple advancing lava flows. Instead, it may be an active lava lake or a complex set of lava flows with episodic, overlapping eruptions. The highest reliable color temperature is ∼1300 K. Although higher temperatures cannot be ruled out, they do not need to be invoked to fit the observed data. The total power output from the active lavas in February 2000 was at least ¹⁰¹¹ W.     

The surface of Io: Geologic units,morphology, and tectonics

A prelimanary geologic map, representing 26.5% of the surface of Io, has been compiled using best-resolution (0.5 to 5 km/line pair) Voyager 1 images and (as a base) a preliminary pictorial map of Io. Nine volcanic units are identified, including materials of mountains (1.9% of total area), plains (49.6%), flows (31.1%), cones (0.1%), and crater vents (4.0%), in addition to seven types of structural features. Photogeologic evidence indicates a dominantly silicate composition for the mountain material, which supports heights of at least 9 ± 1 km. Sulfur flows of diverse viscosity and sulfur-silicate mixtures are thought to compose the pervasive plains. Pit crater and shield crater vent wall scarps reach heights of 2 km and layered plains boundary scarps have estimated heights of 150 to 1700 m; such scarps indicate a material with considerable strenght. A cumulative, volcanic crater size-frequency distribution plot has been prepared using 170 mapped Ionian vents with diameters > 14 km; the shape and slope of the curve are like those for impact craters on other bodies in the solar system, attesting to a similar nonrandom distribution to crater diameters and a surplus of small craters. Io's equatorial zone has six times the number of vents per unit area as the south polar zone. No craters of unequivocal impact origin have been identified on Io to date. A total of 151 lineaments and grabens are recognized with four dominant azimuthal trends forming two nearly orthogonal sets spaces 110° apart (N 85° E, N 25° W and N 45° E, N 55°W). The mapped area lies within the longitudinal zone (250 to 323°) of least-abundant SO₂ frost, indicating that other sulfurous components dominate the upper surface layers in this area.     

Io: Volcanic thermal sources and global heat flow

We have examined thermal emission from 240 active or recently-active volcanic features on Io and quantified the magnitude and distribution of their volcanic heat flow during the Galileo epoch. We use spacecraft data and a geological map of Io to derive an estimate of the maximum possible contribution from small dark areas not detected as thermally active but which nevertheless appear to be sites of recent volcanic activity. We utilize a trend analysis to extrapolate from the smallest detectable volcanic heat sources to these smallest mapped dark areas. Including the additional heat from estimates for “outburst” eruptions and for a multitude of very small (“myriad”) hot spots, we account for ～62 × 10¹² W (～59 ± 7% of Io’s total thermal emission). Loki Patera contributes, on average, 9.6 × 10¹² W (～9.1 ± 1%). All dark paterae contribute 45.3 × 10¹² W (～43 ± 5%). Although dark flow fields cover a much larger area than dark paterae, they contribute only 5.6 × 10¹² W (～5.3 ± 0.6%). Bright paterae contribute ～2.6 × 10¹² W (～2.5 ± 0.3%). Outburst eruption phases and very small hot spots contribute no more than ～4% of Io’s total thermal emission: this is probably a maximum value. About 50% of Io’s volcanic heat flow emanates from only 1.2% of Io’s surface. Of Io’s heat flow, 41 ± 7.0% remains unaccounted for in terms of identified sources. Globally, volcanic heat flow is not uniformly distributed. Power output per unit surface area is slightly biased towards mid-latitudes, although there is a stronger bias toward the northern hemisphere when Loki Patera is included. There is a slight favoring of the northern hemisphere for outbursts where locations were well constrained. Globally, we find peaks in thermal emission at ～315°W and ～105°W (using 30° bins). There is a minimum in thermal emission at around 200°W (almost at the anti-jovian longitude) which is a significant regional difference. These peaks and troughs suggest a shift to the east from predicted global heat flow patterns resulting from tidal heating in an asthenosphere. Global volcanic heat flow is dominated by thermal emission from paterae, especially from Loki Patera (312°W, 12°N). Thermal emission from dark flows maximises between 165°W and 225°W. Finally, it is possible that a multitude of very small hot spots, smaller than the present angular resolution detection limits, and/or cooler, secondary volcanic processes involving sulphurous compounds, may be responsible for at least part of the heat flow that is not associated with known sources. Such activity should be sought out during the next mission to Io.     

Short-term variability of Jupiter’s extended sodium nebula

Ground-based optical observations of D₁ and D₂ line emissions from Jupiter’s sodium nebula, which extend over several hundreds of jovian radii, were carried out at Mt. Haleakala, Maui, Hawaii using a wide field filter imager from May 19 to June 21, 2007. During this observation, the east-west asymmetry of the nebula with respect to the Io’s orbital motion was clearly identified. Particularly, the D₁+D₂ brightness on the western side of Jupiter is strongly controlled by the Io phase angle. The following scenario was developed to explain this phenomenon as follows: First, more ionospheric ions like NaX⁺, which are thought to produce fast neutral sodium atoms due to a dissociative recombination process, are expected to exist in Io’s dayside hemisphere rather than in the nightside one. Second, it is expected that more NaX⁺ ionospheric ions are picked up by the jovian co-rotating magnetic field when Io’s leading hemisphere is illuminated by the Sun. Third, the sodium atom ejection rate varies with respect to Io’s orbital position as a result of the first two points. Model simulations were performed using this scenario. The model results were consistent with the observation results, suggesting that Io’s ionosphere is expected to be controlled by solar radiation just like Earth.     

The distribution of sodium in Io's cloud: Implications

Models for the distribution of sodium in Io's vicinity and in a disk in Io's orbital plane, compared with observational data, support arguments (1) that Io is the source of the sodium, (2) that sodium is ejected from the inside hemisphere and most of the high velocity sodium which is observed is ejected from the leading inside quadrant, (3) that most of the sodium leads Io in Io's vicinity but follows Io at distances of more than 7R_j from Jupiter, (4) that a significant fraction of the sodium flux is ejected at large angles with respect to Io's orbital plane, (5) that the source velocity distribution has a pronounced high-velocity tail, and (6) that impact ionization by electrons is significant at large distances from Io.     

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

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

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.     

The geothermal gradient of Io: Consequences for lithosphere structure and volcanic eruptive activity

We solve numerically the equations describing the transfer of heat through the lithosphere of Io by a mixture of conduction and volcanic advection as proposed by O’Reilly and Davies (O’Reilly, T.C., Davies, G.F. [1981]. Geophys. Res. Lett. 8, 313-316), removing the requirement that average material properties must be used. As expected, the dominance of advective heat transfer by volcanic eruptions means that Io’s geothermal gradient well away from volcanic centres is very small, of order 1 K km⁻¹. This result is independent of any reasonable assumptions about the radiogenic heating rate in the lithosphere. The lithosphere temperature does not increase greatly above the surface temperature until the base of the lithosphere is approached, except in limited areas around shallow magma bodies. As a consequence, solid volatile sulphur compounds mobilized by volcanic processes and re-deposited on the surface of Io commonly remain solid until they reach great depths as they are progressively buried by ongoing activity. For current estimates of the volcanic heat transfer rate, melting of SO₂ does not begin until a depth of ∼20 km and sulphur remains solid to a depth of ∼26 km in a 30 km thick lithosphere. Rising magmas can incorporate fluids from these deep sulphur compound aquifers, and we quantify the major influence that this can have on the bulk density of the magma and hence the resulting possible intrusion and eruption styles.     

Monte Carlo modeling of Io’s [OI] 6300 Å and [SII] 6716 Å auroral emission in eclipse

We present a Monte Carlo (MC) model of [OI] 6300 Å and [SII] 6716 Å emission from Io entering eclipse. The simulation accounts for the 3-D distribution of SO₂, O, SO, S, and O₂ in Io’s atmosphere, several volcanic plumes, and the magnetic field around Io. Thermal electrons from the jovian plasma torus are input along the simulation domain boundaries and move along the magnetic field lines distorted by Io, occasionally participating in collisions with neutrals. We find that the atmospheric asymmetry resulting from varying degrees of atmospheric collapse across Io (due to eclipse ingress) and the presence of volcanoes contributes significantly to the unique morphology of the [OI] 6300 Å emission. The [OI] radiation lifetime of ∼134 s limits the emission to regions that have a sufficiently low neutral density so that intermolecular collisions are rare. We find that at low altitudes (typically <40 km) and in volcanic plumes (Pele, Prometheus, etc.) the number density is large enough (>4 × 10⁹ cm⁻³) to collisionally quench nearly all (>95%) of the excited oxygen for reasonable quenching efficiencies. Upstream (relative to the plasma flow), Io’s perturbation of the jovian magnetic field mirrors electrons with high pitch angles, while downstream collisions can trap the electrons. This magnetic field perturbation is one of the main physical mechanisms that results in the upstream/downstream brightness asymmetry in [OI] emission seen in the observation by Trauger et al. (Trauger, J.T., Stapelfeldt, K.R., Ballester, G.E., Clarke, J.I., 1997. HST observations of [OI] emissions from Io in eclipse. AAS-DPS Abstract (1997DPS29.1802T)). There are two other main causes for the observed brightness asymmetry. First, the observation’s viewing geometry of the wake spot crosses the dayside atmosphere and therefore the wake’s observational field of view includes higher oxygen column density than the upstream side. Second, the phased entry into eclipse results in less atmospheric collapse and thus higher collisional quenching on the upstream side relative to the wake. We compute a location (both in altitude and latitude) for the intense wake emission feature that agrees reasonably well with this observation. Furthermore, the peak intensity of the simulated wake feature is less than that observed by a factor of ∼3, most likely because our model does not include direct dissociation-excitation of SO₂ and SO. We find that the latitudinal location of the emission feature depends not so much on the tilt of the magnetic field as on the relative north/south flux tube depletion that occurs due to Io’s changing magnetic latitude in the plasma torus. From 1-D simulations, we also find that the intensity of [SII] 6716 and 6731 Å emission is much weaker than that of [OI] even if the [SII] excitation cross section is 10³ times larger than excitation to [OI]. This is because the density of S⁺ is much less than that of O and because the Einstein-A coefficient of the [SII] emission is a factor of ∼10 smaller than that of [OI].     

The ion mass loading rate at Io

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