The heartbeat of the volcano: The discovery of episodic activity at Prometheus on Io

The temporal signature of thermal emission from a volcano is a valuable clue to the processes taking place both at and beneath the surface. The Galileo Near Infrared Mapping Spectrometer (NIMS) observed the volcano Prometheus, on the jovian moon Io, on multiple occasions between 1996 and 2002. The 5 micron (μm) brightness of this volcano shows considerable variation from orbit to orbit. Prometheus exhibits increases in thermal emission that indicate episodic (though non-periodic) effusive activity in a manner akin to the current Pu'u 'O'o-Kupaianaha (afterwards referred to as the Pu'u 'O'o) eruption of Kilauea, Hawai'i. The volume of material erupted during one Prometheus eruption episode (defined as the interval from minimum thermal emission to peak and back to minimum) from 6 November 1996 to 7 May 1997 is estimated to be ∼0.8 km³, with a peak instantaneous volumetric flux (effusion rate) of ∼140 m³ s⁻¹, and an averaged volumetric flux (eruption rate) of ∼49 m³ s⁻¹. These quantities are used to model subsurface structure, magma storage and magma supply mechanisms, and likely magma chamber depth. Prometheus appears to be supplied by magma from a relatively shallow magma chamber, with a roof at a minimum depth of ∼2-3 km and a maximum depth of ∼14 km. This is a much shallower depth range than sources of supply proposed for explosive, possibly ultramafic, eruptions at Pillan and Tvashtar. As Prometheus-type effusive activity is widespread on Io, shallow magma chambers containing magma of basaltic or near-basaltic composition and density may be common. This analysis strengthens the analogy between Prometheus and Pu'u 'O'o, at least in terms of eruption style. Even though the style of eruption appears to be similar (effusive emplacement of thin, insulated, compound pahoehoe flows) the scale of activity at Prometheus greatly exceeds current activity at Pu'u 'O'o in terms of volume erupted, area covered, and magma flux. Whereas the estimated magma chamber at Prometheus dwarfs the Pu'u 'O'o magma chamber, it fits within expectations if the Pu'u 'O'o chamber were scaled for the greater volumetric flux and lower gravity of Io. Recent volumetric eruption rates derived from Galileo data for Prometheus were considerably smaller than the rate that produced the extensive flows formed in the ∼17 years between the Voyager and Galileo missions. These smaller eruption rates, coupled with the fact that flows are not expanding laterally, may mean that the immediate heat source that generates the Prometheus plume is simultaneously running out of available volatiles and the thermal energy that drives mobilization of volatiles. This raises the question of whether the current Prometheus eruption is in its last throes.     

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.     

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 evolving flow of Jupiter’s White Ovals and adjacent cyclones

We present results regarding the dynamical meteorology of Jupiter’s White Ovals at different points in their evolution. Starting from the era with three White Ovals FA, BC, and DE (Galileo), continuing to the post-merger epoch with only one Oval BA (Cassini), and finally to Oval BA’s current reddened state (New Horizons), we demonstrate that the dynamics of their flow have similarly evolved along with their appearance. In the Galileo epoch, Oval DE had an elliptical shape with peak zonal wind speeds of ∼90 m s⁻¹ in both its northern and southern peripheries. During the post-merger epoch, Oval BA’s shape was more triangular and less elliptical than Oval DE; in addition to widening in the north-south direction, its northern periphery was 20 m s⁻¹ slower, and its southern periphery was 20 m s⁻¹ faster than Oval DE’s flow during the Galileo era. Finally, in the New Horizons era, the reddened Oval BA had evolved back to a classical elliptical form. The northern periphery of Oval BA increased in speed by 20 m s⁻¹ from Cassini to New Horizons, ending up at a speed nearly identical to that of the northern periphery of Oval DE during Galileo. However, the peak speeds along the southern rim of the newly formed Oval BA were consistently faster than the corresponding speeds in Oval DE, and they increased still further between Cassini and New Horizons, ending up at ∼140-150 m s⁻¹. Relative vorticity maps of Oval BA reveal a cyclonic ring surrounding its outer periphery, similar to the ring present around the Great Red Spot. The cyclonic ring around Oval BA in 2007 appears to be moderately stronger than observed in 1997 and 2001, suggesting that this may be associated with the coloration of the vortex. The modest strengthening of the winds in Oval BA, the appearance of red aerosols, and the appearance of a turbulent, cyclonic feature to Oval BA’s northwest create a strong resemblance with the Great Red Spot from both a dynamical and morphological perspective.In addition to the White Ovals, we also measure the winds within two compact cyclonic regions, one in the Galileo data set and one in the Cassini data set. In the images, these cyclonic features appear turbulent and filamentary, but our wind field reveals that the flow manifests as a coherent high-speed collar surrounding relatively quiescent interiors. Our relative vorticity maps show that the vorticity likewise concentrates in a collar near the outermost periphery, unlike the White Ovals which have peak relative vorticity magnitudes near the center of the vortex. The cyclones contain several localized bright regions consistent with the characteristics of thunderstorms identified in other studies. Although less studied than their anticyclonic cousins, these cyclones may offer crucial insights into the planet’s cloud-level energetics and dynamical meteorology.     

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.     

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.     

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.     

Simulation of Io’s auroral emission: Constraints on the atmosphere in eclipse

We study the morphology of Io’s aurora by comparing simulation results of a three-dimensional (3D) two-fluid plasma model to observations by the high-resolution Long-Range Reconnaissance Imager (LORRI) on-board the New Horizons spacecraft and by the Hubble Space Telescope Advanced Camera for Surveys (HST/ACS). In 2007, Io’s auroral emission in eclipse has been observed simultaneously by LORRI and ACS and the observations revealed detailed features of the aurora, such as a huge glowing plume at the Tvashtar paterae close to the North pole. The auroral radiation is generated in Io’s atmosphere by collisions between impinging magnetospheric electrons and various neutral gas components. We calculate the interaction of the magnetospheric plasma with Io’s atmosphere-ionosphere and simulate the auroral emission. Our aurora model takes into account not only the direct influence of the atmospheric distribution on the morphology and intensity of the emission, but also the indirect influence of the atmosphere on the plasma environment and thus on the exciting electrons. We find that the observed morphology in eclipse can be explained by a smooth (non-patchy) equatorial atmosphere with a vertical column density that corresponds to ∼10% of the column density of the sunlit atmosphere. The atmosphere is asymmetric with two times higher density and extension on the downstream hemisphere. The auroral emission from the Tvashtar volcano enables us to constrain the plume gas content for the first time. According to our model, the observed intensity of the Tvashtar plume implies a mean column density of ∼5 × 10¹⁵ cm⁻² for the plume region.     

Geology and activity around volcanoes on Io from the analysis of NIMS spectral images

The last two successful flybys of Io by Galileo in 2001 (orbits I31, I32) allowed the Near Infrared Mapping Spectrometer to enrich its collection of IR spectral image cubes of the satellite. These data cover hemispheric portions of Io, several volcanic centers as well as their surroundings with a spatial resolution ranging from 2 to 93 km pixel⁻¹. They map thermal emission from the hot-spots and the distribution of solid SO₂ in the 1.0-4.7 μm spectral range. We obtain maps of SO₂ abundance and granularity from the NIMS data using the method of Douté et al. (2002, Icarus 158, 460-482). The maps are correlated to distinguish four different physical units that indicate zones of SO₂ condensation, metamorphism and sublimation. We relate these information with visible images from Galileo's Solid State Imaging System and with detailed mapping of the thermal emission produced by Io's surface. Our principal goal is to understand the mechanisms controlling how lava, pyroclastics and gas are emitted by different types of volcanoes and how these products evolve. The 800 km diameter white ring of fallout created by a violent “Pillanian” eruption during summer of 2001 is at least partly composed of solid SO₂ and has enriched preexisting regional deposits. Orange materials have been recently or are currently emplaced 240 km south from the main eruption site, possibly as sulfur flows. A similar event may have taken place in the past at Ababinili Patera (12.5° N, 142° W). Carefull study of SO₂ maps covering the Emakong region also suggests that sulfur forms the bright channel-fed flow emerging from the south eastern side of the caldera. Within the main caldera of Tvashtar Catena completely cooled patches of crust exist. Elsewhere, the caldera is still cooling from previous episodes of flooding. We confirm that Amirani emits constantly large amount of SO₂ gas by interaction of fresh lava with the volatiles of the underlying plains. Nevertheless SO₂ frost is not the major component of the bright white ring seen in the SSI images. Over the whole Gish Bar region, SO₂ frost seems barely stable and is constantly regenerated. The stability increases along gray filamentary structures which could be faults filled with materials having peculiar thermal properties. Northwest of Gish Bar Patera, a localized bright deposit shows an unusual spectral signature potentially indicative of H₂O molecules forming ice crystals or being trapped in a nonidentified matrix. The Chaac region may present a thickened old crust reducing the geothermal flux to levels lower than 0.5 W m⁻² and thus creating a cold trap for SO₂. Looking at the abundance and degree of metamorphose of SO₂, we establish the relative age of different flows and ejecta for the Sobo Fluctus. Finally the assumption that the white patches in visible images indicate SO₂ rich deposits is once again challenged. In the Camaxtli region we identify a topographically controlled compact white deposit showing only moderate SO₂ abundance. In contrast, we detect two spots of quite pure SO₂ ice on the gray flanks of Emakong. Furthermore, the close association of fumarolic SO₂ and red S₂ already noted for several volcanic centers is observed at Tupan.     

The dust trail of Comet 67P/Churyumov-Gerasimenko

We report the detection of Comet 67P/Churyumov-Gerasimenko's dust trail and nucleus in 24 μm Spitzer Space Telescope images taken February 2004. The dust trail is not found in optical Palomar images taken June 2003. Both the optical and infrared images show a distinct neck-line tail structure, offset from the projected orbit of the comet. We compare our observations to simulated images using a Monte Carlo approach and a dynamical model for comet dust. We estimate the trail to be at least one orbit old (6.6 years) and consist of particles of size ?100 μm. The neck-line is composed of similar sized particles, but younger in age. Together, our observations and simulations suggest grains 100 μm and larger in size dominate the total mass ejected from the comet. The radiometric effective radius of the nucleus is 1.87±0.08 km, derived from the Spitzer observation. The Rosetta spacecraft is expected to arrive at and orbit this comet in 2014. Assuming the trail is comprised solely of 1 mm radius grains, we compute a low probability (∼¹⁰⁻³) of a trail grain impacting with Rosetta during approach and orbit insertion.     

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.     

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.     

Keck AO survey of Io global volcanic activity between 2 and 5 μm

We present in this Keck AO paper the first global high angular resolution observations of Io in three broadband near-infrared filters: Kc (2.3 μm), Lp (3.8 μm), and Ms (4.7 μm). The Keck AO observations are composed of 13 data sets taken during short time intervals spanning 10 nights in December, 2001. The MISTRAL deconvolution process, which is specifically aimed for planetary images, was applied to each image. The spatial resolution achieved with those ground-based observations is 150, 240, and 300 km in the Kc, Lp, and Ms band, respectively, making them similar in quality to most of the distant observations of the Galileo/NIMS instrument. Eleven images per filter were selected and stitched together after being deprojected to build a cylindrical map of the entire surface of the satellite. In Kc-band, surface albedo features, such as paterae (R>60 km) are easily identifiable. The Babbar region is characterized by extremely low albedo at 2.2 μm, and shows an absorption band at 0.9 μm in Galileo/SSI data. These suggest that this region is covered by dark silicate deposits, possibly made of orthopyroxene. In the Lp-Ms (3-5 μm) bands, the thermal emission from active centers is easily identified. We detected 26 hot spots in both broadband filters over the entire surface of the minor planet; two have never been seen active before, nine more are seen in the Ms band. We focused our study on the hot spots detected in both broadband filters. Using the measurements of their brightness, we derived the temperature and area covered by 100 brightness measurements. Loki displayed a relatively quiescent activity. Dazhbog, a new eruption detected by Galileo/NIMS in August 2001, is a major feature in our survey. We also point out the fading of Tvashtar volcanic activity after more than two years of energetic activity, and the presence of a hot, but small, active center at the location of Surt, possibly a remnant of its exceptional eruption detected in February 2001. Two new active centers, labeled F and V on our data, are detected. Using the best temperature and the surface area derived from the L and M band intensities, we calculated the thermal output of each active center. The most energetic hot spots are Loki and Dazhbog, representing respectively 36 and 19% of the total output calculated from a temperature fit of all hot spots (20.6×¹⁰¹² W). Based on the temperature derived from hot spots (∼400 K), our measurement can unambiguously identify the contribution to the heat flux from the silicate portion of the surface. Because the entire surface was observed, no extrapolation was required to calculate that flux. It is also important to note that we measured the brightness of the individual hot spots when they were located close to the Central Meridian. This minimizes the line-of-sight effect which does not follow strictly a classical cosine law. Finally, we argue that despite the widespread volcanic activity detected, Io was relatively quiescent in December 2001, with a minimum mean total output of 0.4-1.2 W m⁻². This output is at least a factor of two lower than those inferred from observations made at longer wavelengths and at different epochs.     

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.     

Astrometry and dynamics of Anthe (S/2007 S 4), a new satellite of Saturn

We describe the astrometry and dynamics of Anthe (S/2007 S 4), a new satellite of Saturn discovered in images obtained using the Imaging Science Subsystem (ISS) of the Cassini spacecraft. Included are details of 63 observations, of which 28 were obtained with Cassini's narrow-angle camera (NAC) and 35 using its wide-angle camera (WAC), covering an observation time-span of approximately 3 years. We estimate the diameter of Anthe to be ∼1.8 km. Orbit modeling based on a numerical integration of the full equations of motion fitted to the observations show that Anthe is in a first-order 11:10 mean motion resonance with Mimas. Two resonant arguments are librating: ?₁=11λ^′−10λ−? and ?₂=11λ^′−10λ−?^′−Ω^′+Ω, where λ, ? and Ω refer to the mean longitude, longitude of pericenter and longitude of ascending node of Mimas and Anthe, with the primed quantities corresponding to Anthe. These resonances cause periodic variations in the orbital elements. The semi-major axis varies by ±26 km over a 913-day period. Anthe is also close to a second-order eccentricity-type mean motion resonant relationship of the form 77:75 with Methone. Since Methone is also in a first-order resonance with Mimas [Spitale, J.N., Jacobson, R.A., Porco, C.C., Owen, W.M., 2006. Astron. J. 132, 692-710], an additional indirect perturbation exists between Methone and Anthe via Mimas. Neither effect is detectable in the orbit fitting and the short-term dynamical evolution of Anthe is dominated by the Mimas-Anthe resonances alone. The expected modulation effect from the Mimas-Tethys 4:2 inclination resonance is also insignificant over this time period. By including Cassini ISS observations of Mimas in the numerical integration fit, we estimate the GM of Mimas to be , consistent with Jacobson et al. [Jacobson, R.A., Spitale, J., Porco, C.C., Owen, W.M., 2006. Astron. J. 132, 711-713].     

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

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

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.     

Search for and limits on plume activity on Mimas, Tethys, and Dione with the Cassini Visual Infrared Mapping Spectrometer (VIMS)

Cassini Visual Infrared Mapping Spectrometer (VIMS) observations of Mimas, Tethys, and Dione obtained during the nominal and extended missions at large solar phase angles were analyzed to search for plume activity. No forward scattered peaks in the solar phase curves of these satellites were detected. The upper limit on water vapor production for Mimas and Tethys is one order of magnitude less than the production for Enceladus. For Dione, the upper limit is two orders of magnitude less, suggesting this world is as inert as Rhea (Pitman, K.M., Buratti, B.J., Mosher, J.A., Bauer, J.M., Momary, T., Brown, R.H., Nicholson, P.D., Hedman, M.M. [2008]. Astrophys. J. Lett. 680, L65-L68). Although the plumes are best seen at ∼2.0 μm, Imaging Science Subsystem (ISS) Narrow Angle Camera images obtained at the same time as the VIMS data were also inspected for these features. None of the Cassini ISS images shows evidence for plumes. The absence of evidence for any Enceladus-like plumes on the medium-sized saturnian satellites cannot absolutely rule out current geologic activity. The activity may below our threshold of detection, or it may be occurring but not captured on the handful of observations at large solar phase angles obtained for each moon. Many VIMS and ISS images of Enceladus at large solar phase angles, for example, do not contain plumes, as the active “tiger stripes” in the south pole region are pointed away from the spacecraft at these times. The 7-year Cassini Solstice Mission is scheduled to gather additional measurements at large solar phase angles that are capable of revealing activity on the saturnian moons.     

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.     

Mapping of the Culann-Tohil region of Io from Galileo imaging data

We have used Galileo spacecraft data to produce a geomorphologic map of the Culann-Tohil region of Io's antijovian hemisphere. This region includes a newly discovered shield volcano, Ts?i Goab Tholus and a neighboring bright flow field, Ts?i Goab Fluctus, the active Culann Patera and the enigmatic Tohil Mons-Radegast Patera-Tohil Patera complex. Analysis of Voyager global color and Galileo Solid-State Imaging (SSI) high-resolution, regional (50-330 m/pixel), and global color (1.4 km/pixel) images, along with available Galileo Near-Infrared Mapping Spectrometer (NIMS) data, suggests that 16 distinct geologic units can be defined and characterized in this region, including 5 types of diffuse deposits. Ts?i Goab Fluctus is the center of a low-temperature hotspot detected by NIMS late during the Galileo mission, and could represent the best case for active effusive sulfur volcanism detected by Galileo. The Culann volcanic center has produced a range of explosive and effusive deposits, including an outer yellowish ring of enhanced sulfur dioxide (SO₂), an inner red ring of SO₂ with short-chain sulfur (S₃-S₄) contaminants, and two irregular green diffuse deposits (one in Tohil Patera) apparently produced by the interaction of dark, silicate lava flows with sulfurous contaminants ballistically-emplaced from Culann's eruption plume(s). Fresh and red-mantled dark lava flows west of the Culann vent can be contrasted with unusual red-brown flows east of the vent. These red-brown flows have a distinct color that is suggestive of a compositional difference, although whether this is due to surface alteration or distinct lava compositions cannot be determined. The main massif of Tohil Mons is covered with ridges and grooves, defining a unit of tectonically disrupted crustal materials. Tohil Mons also contains a younger unit of mottled crustal materials that were displaced by mass wasting processes. Neighboring Radegast Patera contains a NIMS hotspot and a young lava lake of dark silicate flows, whereas the southwest portion of Tohil Patera contains white flow-like units, perhaps consisting of ‘ponds’ of effusively emplaced SO₂. From 0°-15° S the hummocky bright plains unit away from volcanic centers contains scarps, grooves, pits, graben, and channel-like features, some of which have been modified by erosion. Although the most active volcanic centers appear to be found in structural lows (as indicated by mapping of scarps), DEMs derived from stereo images show that, with the exception of Tohil Mons, there is less than 1 km of relief in the Culann-Tohil region. There is no discernable correlation between centers of active volcanism and topography.