Phase curves of materials of Io: Interpretation in terms of Hapke's function

The photometric function developed by B. Hapke (1981,J. Geophys. Res.86, 3039–3054; 1984, Icarus59, 41–59) has been applied to near-opposition (α = 2–8°) disk-resolved phase curves for three color classes on Io, and the disk-integrated phase curve (α = 2–159°) of the satellite as a whole. Derived values of the Hapke compaction parameter h suggest that (1) a large percentage of the material on Io's surface has a porosity significantly greater significantly greater than 60%, supporting the estimate of high porosity made by D.L. Matson and D.B. Nash (1983,J. Geophys. Res.88, 4771–4783) and Nelson et al. (1984, Bull. Amer. Astron. Soc.16, 683–685; 1984,EOS65, 982–983); and (2) Average (“orange”) and Polar (“brown”) materials are significantly more porous than Bright (“white”) materials, a cottrast consistent with the Matson and Nash (1983) SO₂ cold trap model. The best-fit single particle phase function becomes more backscattering on moving from Polar to Average to Bright materials, with the surface of Io on average exhibiting significant backscattering comparable in magnitude to that of the lunar surface. For the color classes, and for Io as a whole, the degree of backscattering tends to increase toward longer wavelengths. The average macroscopic roughness of the Ionian surface, characterized by a mean slope angle of Ø ≃ 25°, is similar to that of other solid surface in the solar system. Consistency between observed limb darkening and that predicted by the Hapke model requires the presence of significant macroscopic roughness (Ø ≥ 20°) for the Average regions, but not necessarily for the Bright and Polar materials.     

The global distribution of sulfur dioxide ice on Io, observed with OSIRIS on the W.M. Keck telescope

We present ground based observations of Io taken with a high spatial resolution imaging spectrometer on 1 and 2 June 2006. We mapped the 1.98 and 2.12 μm absorptions of SO₂ frost, across Io's surface. We analyze these data with surface reflectance modeling using the Hapke method to determine the general frost distribution. This analysis also determined a lower limit of 700 μm on the grain size for the areas of strongest absorption. We incorporate our findings of a predominantly equatorial distribution of SO₂ frost, with the maps of Carlson et al. [Carlson, R.W., Smythe, W.D., Lopes-Gautier, R.M.C., Davies, A.G., Kamp, L.W., Mosher, J.A., Soderblom, L.A., Leader, F.E., Mehlman, R., Clark, R.N., Fanale, F.P., 1997. Geophys. Res. Lett. 24, 2479-2482], McEwen [McEwen, A.S., 1988. Icarus 73, 385-426] and Douté et al. [Douté, S., Schmitt, B., Lopes-Gautier, R., Carlson, R., Soderblom, L., Shirley, J., and The Galileo NIMS Team, 2001. Icarus 149, 107-132] to produce a self consistent explanation of the global distribution of SO₂. We propose that the differences between the above maps is attributable, in part, to the different bands that were studied by the investigators.     

Io's surface composition based on reflectance spectra of sulfur/salt mixtures and proton-irradiation experiments

The available full-disk reflectance spectra of Io in the range 0.3 to 2.5 μm have been interpreted by comparison with new laboratory spectra of a wide variety of natural and synthetic mineral phases in order to determine a surface compositional model for Io that is consistent with Io's other known chemical and physical properties. Our results indicate that the dominant mineral phases are sulfates and free sulfur derived from them, which points toward a low temperature and initially water-rich surface assemblage. Our current preferred mineral phase mixture that best matches the Io data and is simultaneously most consistent with other constraints, consists of a fine-grained particulate mixture of free sulfur (55 vol%), dehydrated bloedite [Na₂Mg(SO₄)₂·xH₂O] (30 vol%) ferric sulfate [Fe₂(SO₄)₃·xH₂O] (15 vol%), and trace amounts of hematite [Fe₂O₃]. Other salts may be present, such as halite and sodium nitrate, as well as clay minerals. Such a model is consistent with a probable pre- and post-accretion thermal history of Io-forming material and Io's observed Na emission and other properties. These results further support the evaporite surface hypothesis of Fanale et al'; while not precluding the presence of certain silicate phases such as montmorillonite.The average surface of Io's leading hemisphere appears to contain less free sulfur and more salts and to be finer grained than that of the trailing hemisphere. Since Io is immersed in Jupiter's magnetosphere, irradiation damage effects from low-energy proton bombardment were studied. Irradiation damage of lattices is estimated to be a relatively minor but operative process on the surface of Io; irradiation darkening by sulfate reduction to free sulfur and by F-center production in salts may be partly responsible for the differences in albedo of leading and trailing hemispheres and equatorial and polar regions of Io, but slight regional differences in relative intrinsic phase concentration on the surface may likewise account for these global variations in albedo.Possible unusual surface properties predicted by this model include: posteclipse darkening in certain wavelenghts, limb brightening in certain wavelengths, and unusual surface electrical properties. Further refinement of Io's surface composition model and better understanding of surface irradiation effects will be possible when observational data in the range 0.20 to 0.30 μm are obtained and when improved spectra in the range 0.30 to 5.0 μm are obtained having increased spectral, spatial, and temporal resolution.     

Optical properties of solid and liquid sulfur at visible and infrared wavelengths

Knowledge of the optical constants of elemental sulfur has potential applications to Venus, Jupiter, Io, Amalthea, and the Earth. The real part, n, of the index of refraction of liquid sulfur (at 133°C) and of solid orthorhombic sulfur (at 25°C) for the wavelength range 0.4–2.0 μm were measured ellipsometrically. The imaginary part, k, of the refractive index of liquid sulfur was obtained by transmittance measurements at the same temperature and wavelength range. The reflectance of semi-infinite slabs of solid and liquid sulfur is calculated using the measured n and k values. We confirm that sulfur melts on Io would be classified as “black” by the Voyager imaging system.     

Near-opposition limb darkening of solids of planetary interest

This paper presents a laboratory study of the limb darkening near opposition, of particulate materials of planetary interest and concentrates on the wavelength dependence of this limb darkening. We find that near zero phase the scattering properties of most particulate materials can be described adequately by Minnaert's law. However, there are materials for which such a representation is totally inadequate. Examples are bronzite and graphite, materials that tend to fracture into flakes having mirrorlike surfaces. In addition, there are materials, such as olivine, whose scattering properties within deep absorption bands show definite departures from Minnaert's law at large angles of incidence or emission. Our Minnaert parameters, k and B₀, were measured at a phase angle of α = 4°. For samples of comparable surface texture and roughness, k and B₀ are approximately linearly related, k increasing as B₀ increases. Very dark materials tend to have k ～ 0.5 to 0.6, while very bright materials tend to have k ～ 1. The linear relationship between k and B₀ can be explained in terms of the varying importance of multiple scattering in the surface layer. Thus for materials for which B₀ is strongly wavelength dependent, so is k. For example, for olivine, k varies from 0.73 to 0.87 between 0.4 and 1.2 μm. These variations are closely correlated with those in B₀: the value of k is relatively high outside of absorption bands and relatively low within them. For bright materials, k is very sensitive to surface roughness. For example, for BaSO₄ powder, k can be changed from ～1.0 to ～0.8 by this effect alone, a fact which has relevance to the photometry of frost-covered satellites. For dark materials, the effects of surface roughness on k are smaller and more subtle.     

Surface changes on Io during the Galileo mission

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

Limb darkening of meteorites and asteroids

This paper presents the results of a laboratory study of the limb darkening, near opposition, of the carbonaceous chondrites Orgueil (C1), Murchison (C2), and Allende (C3), the ordinary chondrite Bruderheim (L6), and a stainless-steel powder. These materials represent possible analogs for the surface materials of C, S, and M asteroids respectively. At low phase angles, the limb-darkening behavior of all materials studied is well represented by Minnaert's law. For carbonaceous chondrites, the Minnaert limb-darkening parameter k is nearly independent of wavelength for wavelengths between 0.4 and 0.9 μm, with a typical value of k = 0.55. The reflectance parameter, B₀, varies from 0.045 to 0.065 over the same range of wavelengths. Both k and B₀ are larger for the stainless-steel powder and the ordinary chondrite, due to the increased importance of multiple scattering in the surface layer. If no limb darkening were present, k would equal $\text{1}\text{2}$ and the geometric albedo (p) of an asteroid would equal the normal reflectance $\text{(r}{}_{\mathrm{<mtext>n</mtext>}}\text{? B}{}_0\text{)}$ of its surface material. For bodies whose surface material is appreciably limb darkened, the geometric albedo measured at the telescope will be lower than the true normal reflectance of surface material; we estimate that for S and M objects r_n ? 1.05 p. In the case of nonspherical asteroids, because the distribution of incidence and emission angles varies as the asteroid rotates, the geometric albedo must change with aspect. If limb darkening is not considered when interpreting asteroid light curves, the values of b/a derived will be too extreme. This effect is probably too small to be observed for C asteroids, because of their intrinsically low reflectances, but could be appreciable for S and M objects.     

The reflection spectrum of liquid sulfur: Implications for Io

The spectral reflectance from 0.38 to 0.75 μm of a column of liquid sulfur has been measured at several temperatures between the melting point (～118°C) and 173°C. Below 160°C the spectral reflectance was observed to vary reversibly as a function of temperature, independent of the previous thermal history of the column. Once the temperature exceeded 160°C, the spectrum would not change given a subsequent decrease in temperature. The spectral reflectance of the liquid-sulfur column at all temperatures was very low (10–19%). Combining this information with Voyager spectrophotometry of Jupiter's satellite Io, it is concluded that liquid sulfur at any temperature on Io's surface would be classified as a “black area” according to the standards used by the Voyager imaging team in their spectrophotometric analysis (L. Soderblom, T. V. Johnson, D. Morrison, E. Danielson, B. L. Smith, J. Veverka, A. Cook, C. Sagan, P. Kupferman, D. Pieri, J. Mosher, C. Avis, J. Gradie, and T. Clancy (1980). Geophys. Res. Lett.7, 963–966).     

Voyager disk-integrated photometry of Io

Voyager full-disk images of Io, available at solar phase angle of α = 2?29° and 101?159°, allow comparisons of the satellite's near-opposition photometric behavior with Earth-based results and the determination of the phase curve out to very high phase angles. The near-opposition data were reduced iteratively for self-consistent phase and rotation curves in each Voyager filter; the resulting phase coefficients, geometric albedos, and rotational lightcurves are consistent with Earth-based findings, except for a previously noted tendency for Voyager to yield somewhat redder spectral information. The derived near-opposition phase coefficients, ranging between 0.016 and 0.024 mag/ deg, decrease with increasing wavelength, a trend weakly noted in some Earth-based observations. The full, α = 2?159° phase curves allow the first direct determination of the phase integral of Io at several wavelengths: q rises from ≈0.7 in the ultraviolet to ≈0.8 in the orange. Combination of the Voyager phase integrals with Earth-based albedo information leads to a best estimate of the bolometric Bond albedo of 0.50 ± 0.10, a value consistent with, but slightly below, previous estimates.     

Sodium remote from Io

We find that faint sodium emission originating in the middle Jupiter magnetosphere has two distinct kinematical components. The “normal” signature of atoms on bound orbits with large apojoves seems always to be present, and we suggest these atoms are an extension of the bright, near-Io sodium cloud. The “fast” signature, with speeds up to at least 100 km sec^?1, is seen only occasionally, and we suggest it is due to an interaction of the near-Io sodium cloud with the corotating, heavy-ion plasma. Both elastic and charge-exchange collisions seem consistent with the observed kinematical and temporal signatures. Elastic collisions seem marginally more capable of producing the high observed sodium atom speeds. We predict observable occurences of the fast component in the hours following passage of the Io sodium cloud through the plasma centrifugal symmetry surface if Io is at a favorable orbital longitude. Between 10 and 20 R_J we find an atomic sodium density ～10^?2 cm^?3. If the photoionization lifetime applies, an Io source of at least 10²⁶ sodium atoms sec^? is required to maintain this remote sodium population.     

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.     

The Sun and the transcendental world of binaries

The theory of gravity says that a binary with orbital frequency ν should be a source of gravitational waves at the double frequency and higher harmonics. This implies that long-term exposure of an ensemble of binaries to gravity waves with frequency ν _G can result in (a) a lack of binaries with frequencies near frequency ν _G /2 and its higher harmonics (the effect of unstable orbits) and/or (b) an excess of binaries whose orbital frequencies are “absolutely” incommensurable with ν _G /2 and its higher harmonics (the effect of stable orbits). It is assumed that the stable-orbit frequencies are almost equal to multiples of πν _G /2 and ν _G /2π, where π plays the role of a “perfect” factor ensuring the best antiresonance of binaries. The statistical analysis of frequencies of 5774 Galactic close binary systems (CBSs) with periods P less than 10 days is based on calculating the resonance spectrum that indicates the best common multiple for a given set of frequencies with allowance for the factor π. The CBS distribution turns out to be modulated by the frequency ν _G = 104.4(5) μHz, and this effect is the most pronounced for superfast and compact rotators, such as cataclysmic variables (CVs) and related objects. The frequency ν _G is, within the error, equal to the “enigmatic” frequency ν₀ = 104.160(1) μHz com discovered earlier in the power spectra of the Sun and brightness variations of some extragalactic sources. This confirms the existence of a “coherent cosmic oscillation” of the Universe with frequency ν₀(ν _G ). The new astrophysical phenomenon naturally explains an “CV period gap” at frequencies ≈πν _G /3 (P ≈ 153 min) and maxima at the neighboring frequencies ≈πν _G /2 and ≈πν _G /4 (P ≈ 102 and ≈204 min, respectively). The remarkable and “mysterious” role of the transcendental number π for the world of binaries is emphasized, and the mystery of physical nature of the “universal” oscillation ν₀(ν _G ) is highlighted.     

Moments of inertia and the Chandlerian period for two-and three-layer models of the Galilean satellite Io

We consider two-layer (Fe-FeS core+silicate mantle) and three-layer (Fe-FeS core+silicate mantle+crust) models of the Galilean satellite Io. Two parameters are known from observations for the equilibrium figure of the satellite, the mean density ρ₀ and the Love number k₂. Previously, the Radau-Darwin formula was used to determine the mean moment of inertia. Using formulas of the Figure Theory, we calculated the principal moments of inertia A, B, and C and the mean moment of inertia I for the two-and three-layer models of Io using ρ₀ and k₂ as the boundary conditions. We concluded that when modeling the internal structure of Io, it is better to use the observed value of k₂ than the moment of inertia I derived from k₂ using the Radau-Darwin formula. For the models under consideration, we calculated the Chandlerian wobble periods of Io. For the three-layer model, this period is approximately 460 days.     

Voyager disk resolved photometry of the Saturnian satellites

Photometric analysis of Voyager images of the medium-sized icy satellites of Saturn shows that their surfaces exhibit a wide range of scattering properties. At low phase angles, Rhea and Dione closely follow lunar behavior with almost no limb darkening. Mimas, Tethys, and especially Enceladus shiw significant limb darkening at low phase angles, which suggests multiple scattering is important for their surfaces. A simple photometric function of the form I/F = f(α)Aμ₀/(μ + μ₀) + (1 ? A)μ₀ has been fit to the observations. For normal reflectances <0.6, we find lunar-like scattering properties (A = 1). No satellite's surface can be described by Lambert's Law (A = 0). Dione exhibits the widest albedo variations (about 50%). A longitudinal dark stripe which represents a 15% decrease in albedo is situated near the center of the trailing side of Tethys. A correlation is found between the albedo and color of the satellites: the darker objects are redder. Similarly, darker areas of each satellite are redder. Spectral reflectances of Mimas and Enceladus can be derived for the first time. After the proper calibrations to the Voyager color images are made, it is found that both satellites have remarkably flat spectra into the ultraviolet.     

Io's 4-μm band and the role of adsorbed SO2

The role of adsorbed SO₂ on Io's surface particles in producing the observed spectral absorption band near 4 μm in Io's reflectance spectrum is explored. Calculations show that a modest 50% monolayer coating of adsorbed SO₂ molecules on submicron grains of sulfur of alkali sulfide, assumed to make up Io's uppermost optical surface (“radialith”), will result in a ν₁ + ν₃ absorption band near 4 μm with depth ～30% below the adjacent continuum, consistent with the observed strength of the Io band. The precise wavelength position of the ν₁ + ν₃ band of SO₂ in different phase states such as frost, ice, adsorbate, and gas are summarized from the experimental literature and compared with the available telescopic measurements of the Io band position. The results suggest that the 4-μm band in Io's full disk spectrum can best be explained by the presence on Io's surface of widespread SO₂ in the form of adsorbate rather than ice or frost.     

Observational tests for sulfur allotropes on Io

One of the intrinsic properties of particulate sulfur allotropes is a change in UV-visible reflectivity with temperature change of the material. The surface of Io experiences temperature changes during eclipse which are sufficient to cause a detectable change in the spectral reflectivity of sulfur; thus, if the surface of Io is composed primarily of sulfur allotropes, a change in reflectivity at certain wavelengths should be observable shortly after eclipse reappearance. We observed four eclipse reappearances during July and August of 1983 and saw no posteclipse brightening effects in filter bands selected for sensitivity to color changes in sulfur. Our model of the brightness change for S₈ (“yellow” sulfur) implies that this material covers less than 50% of Io's surface. Negative posteclipse brightening observations were also obtained with a filter chosen for the high contrast between SO₂ frost and the average albedo of the surface of Io at that wavelength. We conclude that no significant condensation of optically thick SO₂ occurred on the surface of Io during these eclipses.     

Sputtering of sulfur by kiloelectronvolt ions: Application to the magnetospheric plasma interaction with Io

Measurements of total yields, temperature dependences, mass spectra, and energy spectra of molecules sputtered from condensed sulfur (S₈) at low temperatures by keV ions are reported and results are given for Jovian plasma ion bombardment of Io. A change in the reflectance of the sulfur, which can be removed by annealing, is produced by the most penetrating ions and may be connected with the darker, colder polar regions on Io. The measured sputtering yields are much lower than those estimated earlier for room temperature sulfur films but are comparable to previous measurements of keV ion sputtering of SO₂ at low temperatures. The corrected mass spectrum indicates that ≈66% of the total yield corresponds to S₂ ejection while only 5 and 16% correspond to S and S₃, respectively. Therefore, if ions reach the surface of Io its atmosphere will have a non-negligible sulfur component of primarily S₂. The ejection of S and S₂ is temperature independent for temperatures characteristic of most of the surface of Io. The energy spectrum for S has an approximate 1/E² dependence at high ejection energies, whereas S₂ and S₃ fall off more rapidly. Assuming 50% coverage of both sulfur and SO₂ and a thin atmosphere (e.g., nightside and polar region) the direct sputter injection of sulfur atoms and molecules into the Jovian plasma torus and the indirect injection due to coronal processes are estimated. These injection rates for sulfur are compared to those for SO₂ showing that injection from sulfur deposits contributes 13% to the total mass injection rate of ∼2–3 × 10²⁹ amu/sec.     

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.     

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.     

Thermal history of the Jovian satellite Io

The discovery of large volcanic eruptions on Io suggests that Io is one of the most geologically active planetary bodies. The energy source of this geologic activity is believed to be tidal heating induced by Jupiter. A number of thermal history calculations were done to investigate the effect of tidal heating on the thermal history of Io taking into account solid state convection and advective heat transfer. These simulations show that the total tidal heating energy in Io is almost equal to the advectively transferred heat, indicating that the observed heat flow from Io is nearly equal to the total tidal heating energy. Since total tidal heating energy is dependent on the radius of the liquid mantle and the internal dissipation factor (Q), the radius of the liquid mantle can be estimated for a given value of Q. Some reasonable thermal history models of Io were obtained using a model with Q ≈ 25–50 in which the magma source of Ionian volcanism is at a depth of 100–300 km. The models satisfy the heat flow data and the existence of a thick lithosphere. Using a model with Q = 25 and L = 300 km (thickness of the advective region) as the standard model (model II), we then studied the effect of convective heat transfer and the initial temperature distribution on the Ionian thermal history. In these calculations, the other parameters are the same as in the standard model (model II). These calculations show that although the temperature distribution in the central region reflects the difference in the efficiency of convective heat transfer and initial temperature distribution, the temperature distribution in the outer region does not changes appreciably.