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.     

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.     

A comprehensive numerical simulation of Io’s sublimation-driven atmosphere

Io’s sublimation-driven atmosphere is modeled using the direct simulation Monte Carlo (DSMC) method. These rarefied gas dynamics simulations improve upon earlier models by using a three-dimensional domain encompassing the entire planet computed in parallel. The effects of plasma heating, planetary rotation, inhomogeneous surface frost, molecular residence time of SO₂ on the exposed (non-volatile) rocky surface, and surface temperature distribution are investigated. Circumplanetary flow is predicted to develop from the warm dayside toward the cooler nightside. Io’s rotation leads to a highly asymmetric frost surface temperature distribution (due to the frost’s high thermal inertia) which results in circumplanetary flow that is not axi-symmetric about the subsolar point. The non-equilibrium thermal structure of the atmosphere, specifically vibrational and rotational temperatures, is also examined. Plasma heating is found to significantly inflate the atmosphere on both the dayside and nightside. The plasma energy flux causes high temperatures at high altitudes but plasma energy depletion through the dense gas column above the warmest frost permits gas temperatures cooler than the surface at low altitudes. A frost map (Douté, S., Schmitt, B., Lopes-Gautier, R., Carlson, R., Soderblom, L., Shirley, J., and the Galileo NIMS Team [2001]. Icarus 149, 107-132) is used to control the sublimated flux of SO₂ which can result in inhomogeneous column densities that vary by nearly a factor of four for the same surface temperature. A short residence time for SO₂ molecules on the “rock” component is found to smooth lateral atmospheric inhomogeneities caused by variations in the surface frost distribution, creating an atmosphere that looks nearly identical to one with uniform frost coverage. A longer residence time is found to agree better with mid-infrared observations (Spencer, J.R., Lellouch, E., Richter, M.J., López-Valverde, M.A., Jessup, K.L, Greathouse, T.K., Flaud, J. [2005]. Icarus 176, 283-304) and reproduce the observed anti-jovian/sub-jovian column density asymmetry. The computed peak dayside column density for Io assuming a surface frost temperature of 115 K agrees with those suggested by Lyman-α observations (Feaga, L.M., McGrath, M., Feldman, P.D. [2009]. Icarus 201, 570-584). On the other hand, the peak dayside column density at 120 K is a factor of five larger and is higher than the upper range of observations (Jessup, K.L., Spencer, J.R., Ballester, G.E., Howell, R.R., Roesler, F., Vigel, M., Yelle, R. [2004]. Icarus 169, 197-215; Spencer et al., 2005).     

Observations and temperatures of Io's Pele Patera from Cassini and Galileo spacecraft images

Pele has been the most intense high-temperature hotspot on Io to be continuously active during the Galileo monitoring from 1996-2001. A suite of characteristics suggests that Pele is an active lava lake inside a volcanic depression. In 2000-2001, Pele was observed by two spacecraft, Cassini and Galileo. The Cassini observations revealed that Pele is variable in activity over timescales of minutes, typical of active lava lakes in Hawaii and Ethiopia. These observations also revealed that the short-wavelength thermal emission from Pele decreases with rotation of Io by a factor significantly greater than the cosine of the emission angle, and that the color temperature becomes more variable and hotter at high emission angles. This behavior suggests that a significant portion of the visible thermal emission from Pele comes from lava fountains within a topographically confined lava body. High spatial resolution, nightside images from a Galileo flyby in October 2001 revealed a large, relatively cool (<800 K) region, ringed by bright hotspots, and a central region of high thermal emission, which is hypothesized to be due to fountaining and convection in the lava lake. Images taken through different filters revealed color temperatures of 1500±80 K from Cassini ISS data and 1605±220 and 1420±100 K from small portions of Galileo SSI data. Such temperatures are near the upper limit for basaltic compositions. Given the limitations of deriving lava eruption temperature in the absence of in situ measurement, it is possible that Pele has lavas with ultramafic compositions. The long-lived, vigorous activity of what is most likely an actively overturning lava lake in Pele Patera indicates that there is a strong connection to a large, stable magma source region.     

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

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

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.     

High-Resolution Keck Adaptive Optics Imaging of Violent Volcanic Activity on Io

Io, the innermost Galilean satellite of Jupiter, is a fascinating world. Data taken by Voyager and Galileo instruments have established that it is by far the most volcanic body in the Solar System and suggest that the nature of this volcanism could radically differ from volcanism on Earth. We report on near-IR observations taken in February 2001 from the Earth-based 10-m W. M. Keck II telescope using its adaptive optics system. After application of an appropriate deconvolution technique (MISTRAL), the resolution, ∼100 km on Io's disk, compares well with the best Galileo/NIMS resolution for global imaging and allows us for the first time to investigate the very nature of individual eruptions. On 19 February, we detected two volcanoes, Amirani and Tvashtar, with temperatures differing from the Galileo observations. On 20 February, we noticed a slight brightening near the Surt volcano. Two days later it had turned into an extremely bright volcanic outburst. The hot spot temperatures (>1400 K) are consistent with a basaltic eruption and, being lower limits, do not exclude an ultramafic eruption. These outburst data have been fitted with a silicate-cooling model, which indicates that this is a highly vigorous eruption with a highly dynamic emplacement mechanism, akin to fire-fountaining. Its integrated thermal output was close to the total estimated output of Io, making this the largest ionian thermal outburst yet witnessed.     

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

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

Eclipse reappearances of Io: Time-resolved spectroscopy (1.9-4.2 μm)

We obtained time-resolved, near-infrared spectra of Io during the 60-90 min following its reappearance from eclipse by Jupiter on five occasions in 2004. The purpose was to search for spectral changes, particularly in the well-known SO₂ frost absorption bands, that would indicate surface-atmosphere exchange of gaseous SO₂ induced by temperature changes during eclipse. These observations were a follow-on to eclipse spectroscopy observations in which Bellucci et al. [Bellucci et al., 2004. Icarus 172, 141-148] reported significant changes in the strengths of two strong SO₂ bands in data acquired with the VIMS instrument aboard the Cassini spacecraft. One of the bands (4.07 μm [ν₁ + ν₃]) observed by Bellucci et al. is visible from ground-based observatories and is included in our data. We detected no changes in Io’s spectrum at any of the five observed events during the approximately 60-90 min during which spectra were obtained following Io’s emergence from Jupiter’s shadow. The areas of the three strongest SO₂ bands in the region 3.5-4.15 μm were measured for each spectrum; the variation of the band areas with time does not exceed that which can be explained by the Io’s few degrees of axial rotation during the intervals of observation, and in no case does the change in band strength approach that seen in the Cassini VIMS data. Our data are of sufficient quality and resolution to show the weak 2.198 μm (4549.6 cm⁻¹) 4ν₁ band of SO₂ frost on Io for what we believe is the first time. At one of the events (June 22, 2004), we began the acquisition of spectra ∼6 min before Io reappeared from Jupiter’s shadow, during which time it was detected through its own thermal emission. No SO₂ bands were superimposed on the purely thermal spectrum on this occasion, suggesting that the upper limit to condensed SO₂ in the vertical column above Io’s surface was ∼4 × 10⁻⁵ g cm⁻².     

The nature of the volcanic activity at Loki: Insights from Galileo NIMS and PPR data

Loki is the largest patera and the most energetic hotspot on Jupiter's moon Io, in turn the most volcanically active body in the Solar System, but the nature of the activity remains enigmatic. We present detailed analysis of Galileo Near-Infrared Mapping Spectrometer (NIMS) and PhotoPolarimeter/Radiometer (PPR) observations covering the 1.5-100 μm wavelength range during the I24, I27, and I32 flybys. The general pattern of activity during these flybys is consistent with previously proposed models of a resurfacing wave periodically crossing a silicate lava lake. In particular our analysis of the I32 NIMS observations shows, over much of the observed patera, surface temperatures and implied ages closely matching those expected for a wave advancing counterclockwise at 0.94-1.38 km/day. The age pattern is different than other published analyses which do not show as clearly this azimuthal pattern. Our analysis also shows two additional distinctly different patera surfaces. The first is located along the inner and outer margins where components with a 3.00-4.70-μm color temperature of 425 K exist. The second is located at the southwestern margin where components with a 550-K color temperature exist. Although the high temperatures could be caused by disruption of a lava lake crust, some additional mechanism is required to explain why the southwest margin is different from the inner or outer ones. Finally, analysis of the temperature profiles across the patera reveal a smoothness that is difficult to explain by simple lava cooling models. Paradoxically, at a subpixel level, wide temperature distributions exist which may be difficult to explain by just the presence of hot cracks in the lava crust. The resurfacing wave and lava cooling models explain well the overall characteristics of the observations. However, additional physical processes, perhaps involving heat transport by volatiles, are needed to explain the more subtle features.     

Hydrated minerals on Europa’s surface: An improved look from the Galileo NIMS investigation

The surface composition of Europa is of great importance for understanding both the internal evolution of Europa and its putative ocean. The Near Infrared Mapping Spectrometer (NIMS) investigation on Galileo observed Europa and the other Galilean satellites from 0.7 to 5.2 μm with spatial resolution down to a few kilometers during flybys by the spacecraft as it orbited Jupiter. These data have been analyzed and results published over the life of the Galileo mission and afterward. One result was the discovery of hydrated minerals at some locations on Europa and Ganymede. The data are noisy, especially for Europa, due to radiation affecting the NIMS electronics and detectors, and other artifacts are also present. The NIMS data are now being reprocessed using the accumulated knowledge gained over the entire missions to remove noise spikes and compensate for some other defects in the data. We are analyzing these reprocessed data in an attempt to defined better the nature of the hydrate spectral features and improve their interpretation. We report here on analyses of two NIMS reprocessed observations for the 0.7-3-μm region. A revised hydrate spectrum is calculated and mapped in detail across two lineaments. The spectrum shows the expected distorted water features but little or no spectral structure in these features. A narrow, weak spectral feature appears at 1.344 μm, which is weakly correlated with lower albedo. Several other weak features may be present but are difficult to confirm in these limited data sets. The hydrate signature shows the greatest strength within and toward the center of the lineaments, confirming and strengthening the association of the hydrate with these endogenic features. This trend may indicate that the material in the lineaments is youngest toward the center and has more water frost coverage toward the edge. A small, visually dark, circular feature has a spectrum that shows both hydrate and crystalline water ice features and perhaps contains a hydrate different in spectral characteristics and perhaps composition than found in the lineament.     

Ground-based observations of time variability in multiple active volcanoes on Io

Since before the beginning of the Galileo spacecraft’s Jupiter orbital tour, we have observed Io from the ground using NASA’s Infrared Telescope Facility (IRTF). We obtained images of Io in reflected sunlight and in-eclipse at 2.3, 3.5, and 4.8 μm. In addition, we have measured the 3.5 μm brightness of an eclipsed Io as it is occulted by Jupiter. These lightcurves enable us to measure the brightness and one-dimensional location of active volcanoes on the surface. During the Galileo era, two volcanoes were observed to be regularly active: Loki and either Kanehekili and/or Janus. At least 12 other active volcanoes were observed for shorter periods of time, including one distinguishable in images that include reflected sunlight. These data can be used to compare volcano types and test volcano eruption models, such as the lava lake model for Loki.     

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.     

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.     

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.     

New insights on the binary Asteroid 121 Hermione

We report on the results of a 6-month photometric study of the main-belt binary C-type Asteroid 121 Hermione, performed during its 2007 opposition. We took advantage of the rare observational opportunity afforded by one of the annual equinoxes of Hermione occurring close to its opposition in June 2007. The equinox provides an edge-on aspect for an Earth-based observer, which is well suited to a thorough study of Hermione’s physical characteristics. The catalog of observations carried out with small telescopes is presented in this work, together with new adaptive optics (AO) imaging obtained between 2005 and 2008 with the Yepun 8-m VLT telescope and the 10-m Keck telescope. The most striking result is confirmation that Hermione is a bifurcated and elongated body, as suggested by Marchis, et al. [Marchis, F., Hestroffer, D., Descamps, P., Berthier, J., Laver, C., de Pater, I., 2005. Icarus 178, 450-464]. A new effective diameter of 187 ± 6 km was calculated from the combination of AO, photometric and thermal observations. The new diameter is some 10% smaller than the hitherto accepted radiometric diameter based on IRAS data. The reason for the discrepancy is that IRAS viewed the system almost pole-on. New thermal observations with the Spitzer Space Telescope agree with the diameter derived from AO and lightcurve observations. On the basis of the new AO astrometric observations of the small 32-km diameter satellite we have refined the orbit solution and derived a new value of the bulk density of Hermione of 1.4 + 0.5/−0.2 g cm⁻³. We infer a macroscopic porosity of ∼33 + 5/−20%.     

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.     

Radar observations of Asteroids 64 Angelina and 69 Hesperia

We report new radar observations of E-class Asteroid 64 Angelina and M-class Asteroid 69 Hesperia obtained with the Arecibo Observatory S-band radar (2480 MHz, 12.6 cm). Our measurements of Angelina’s radar bandwidth are consistent with reported diameters and poles. We find Angelina’s circular polarization ratio to be 0.8 ± 0.1, tied with 434 Hungaria for the highest value observed for main-belt asteroids and consistent with the high values observed for all E-class asteroids (Benner, L.A.M., Ostro, S.J., Magri, C., Nolan, M.C., Howell, E.S., Giorgini, J.D., Jurgens, R.F., Margot, J.L., Taylor, P.A., Busch, M.W., Shepard, M.K. [2008]. Icarus 198, 294-304; Shepard, M.K., Kressler, K.M., Clark, B.E., Ockert-Bell, M.E., Nolan, M.C., Howell, E.S., Magri, C., Giorgini, J.D., Benner, L.A.M., Ostro, S.J. [2008b]. Icarus 195, 220-225). Our radar observations of 69 Hesperia, combined with lightcurve-based shape models, lead to a diameter estimate, D_eff = 110 ± 15 km, approximately 20% smaller than the reported IRAS value. We estimate Hesperia to have a radar albedo of , consistent with a high-metal content. We therefore add 69 Hesperia to the Mm-class (high metal M) (Shepard, M.K., Clark, B.E., Ockert-Bell, M., Nolan, M.C., Howell, E.S., Magri, C., Giorgini, J.D., Benner, L.A.M., Ostro, S.J., Harris, A.W., Warner, B.D., Stephens, R.D., Mueller, M. [2010]. Icarus 208, 221-237), bringing the total number of Mm-class objects to eight; this is 40% of all M-class asteroids observed by radar to date.