Lava lakes on Io: New perspectives from modeling

Loki Patera (310° W, 12° N) is Io's largest patera at ∼180 km in diameter. Its morphology and distinct thermal behavior have led researchers to hypothesize that Loki Patera may either be an active lava lake that experiences periodic overturn, or a shallow depression whose floor is episodically resurfaced with thin flows. Using results from mathematical models, we suggest that a better model for Loki's behavior is the terrestrial superfast spreading East Pacific Rise (EPR), near 17°^30′ south. We propose that, like at the southern EPR, Loki Patera is underlain by a thin, persistent magma “lens” that feeds thin, temporary lava lakes within the patera. Also like the southern EPR, overspilling of the volcanic depression is rare, with most of the lava volume being emplaced via a subsurface network of lava tubes.     

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.     

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.     

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.     

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.     

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.     

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.     

Thermal relaxation in young neutron stars

The internal properties of the neutron star crust can be probed by observing the epoch of thermal relaxation. After the supernova explosion, powerful neutrino emission quickly cools the stellar core, while the crust stays hot. The cooling wave then propagates through the crust, as a result of its finite thermal conductivity. When the cooling wave reaches the surface (age 10–100 yr) , the effective temperature drops sharply from 250 eV to 30 or 100 eV, depending on the cooling model. The crust relaxation time is sensitive to the (poorly known) microscopic properties of matter of subnuclear density, such as the heat capacity, thermal conductivity, and superfluidity of free neutrons. We calculate the cooling models with the new values of the electron thermal conductivity in the inner crust, based on a realistic treatment of the shapes of atomic nuclei. Superfluid effects may shorten the relaxation time by a factor of 4. The comparison of theoretical cooling curves with observations provides a potentially powerful method of studying the properties of the neutron superfluid and highly unusual atomic nuclei in the inner crust.     

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.     

Shoreline features of Titan's Ontario Lacus from Cassini/VIMS observations

We analyze observations of Titan's south polar lake Ontario Lacus obtained by Cassini's Visual and Infrared Mapping Spectrometer during the 38th flyby of Titan (T38; 2007 December 5). These near-closest-approach observations have the highest signal-to-noise, the finest spatial resolution, and the least atmospheric influence of any near-infrared lake observation to date. We use the large, spatially flat, and low-albedo interior of Ontario Lacus as a calibration target allowing us to derive an analytical atmospheric correction for emission angle. The dark lake interior is surrounded by two separate annuli that follow the lake interior's contours. The inner annulus is uniformly dark, but not so much as the interior lake, and is generally 5-10 kilometers wide at the lake's southeastern margin. We propose that it represents wet lakebed sediments exposed by either tidal sloshing of the lake or seasonal methane loss leading to lower lake-volume. The exterior annulus is bright and shows a spectrum consistent with a relatively low water-ice content relative to the rest of Titan. It may represent fine-grained condensate deposits from a past era of higher lake level. Together, the annuli seem to indicate that the lake level for Ontario Lacus has changed over time. This hypothesis can be tested with observations scheduled for future Titan flybys.     

The equatorial shape and gravity field of Mercury from MESSENGER flybys 1 and 2

On 14 January and 6 October 2008 the MESSENGER spacecraft passed within 200 km of the surface of Mercury. These flybys by MESSENGER provided the first observations of Mercury from a spacecraft since the Mariner 10 flybys in 1974 and 1975. Data from the Mercury Laser Altimeter (MLA) provided new information on the equatorial shape of Mercury, and Doppler tracking of the spacecraft through the flybys provided new data on the planet’s gravity field. The MLA passes were on opposite hemispheres of the planet and span collectively ∼40% of the equatorial circumference. The mean elevation of topography observed during flyby 1, in the longitude range 0-90°E, is greater than that seen during flyby 2 in the longitude range 180-270°E, indicating an offset between centers of mass and figure having a magnitude and phase in general agreement with topography determined by Earth-based radar. Both MLA profiles are characterized by slopes of ∼0.015° downward to the east, which is consistent with a long-wavelength equatorial shape defined by a best-fitting ellipse. The Doppler tracking data show sensitivity to the gravitational structure of Mercury. The equatorial ellipticity of the gravitational field, C_2,2, is well determined and correlates with the equatorial shape. The S_2,2 coefficient is ∼0, as would be expected if Mercury’s coordinate system, defined by its rotational state, is aligned along its principal axes of inertia. The recovered value of the polar flattening of the gravitational potential, J₂, is considerably lower in magnitude than the value obtained from Mariner 10 tracking, a result that is problematic for internal structure models. This parameter is not as well constrained as the equatorial ellipticity because the flyby trajectories were nearly in the planet’s equatorial plane. The residuals from the Doppler tracking data suggest the possibility of mascons on Mercury, but flyby observations are of insufficient resolution for confident recovery. For a range of assumptions on degree of compensation and crustal and mantle densities, the allowable crustal thickness is consistent with the upper limit of about 100 km estimated from the inferred depth of faulting beneath a prominent lobate scarp, an assumed ductile flow law for crustal material, and the condition that temperature at the base of the crust does not exceed the solidus temperature. The MESSENGER value of C_2,2 has allowed an improved estimate of the ratio of the polar moment of inertia of the mantle and crust to the full polar moment (C_m/C), a refinement that strengthens the conclusion that Mercury has at present a fluid outer core.     

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.     

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.     

Geologic mapping of the Zal region of Io

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

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.     

Flybys in the planar, circular, restricted, three-body problem

An analysis is presented of gravity assisted flybys in the planar, circular, restricted three-body problem (pcr3bp) that is inspired by the Keplerian map and by the Tisserand- Poincaré graph. The new Flyby map is defined and used to give insight on the flyby dynamics and on the accuracy of the linked-conics model. The first main result of this work is using the Flyby map to extend the functionality of the Tisserand graph to low energies beyond the validity of linked conics. Two families of flybys are identified: Type I (direct) flybys and Type II (retrograde) flybys. The second main result of this work shows that Type I flybys exist at all energies and are more efficient than Type II flybys, when both exist. The third main result of this work is the introduction of a new model, called ??Conics, When I Can??, which mixes numerical integration and patched conics formulas, and has applications beyond the scope of this work. The last main result is an example trajectory with multiple flybys at Ganymede, all outside the linked-conics domain of applicability. The trajectory is computed with the pcr3bp, and connects an initial orbit around Jupiter intersecting the Callisto orbit, to an approach transfer to Europa. Although the trajectory presented has similar time of flight and radiation dose of other solutions found in literature, the orbit insertion ??v is 150 m/s lower. For this reason, the transfer is included in the lander option of the Europa Habitability Mission Study.     

The Retrieval of Cloud Structure Maps in the Equatorial Region of Jupiter Using a Principal Component Analysis of Galileo/NIMS Data

The cloud structure of the jovian atmosphere at pressures less than 2 bars has previously been estimated using near-infrared observations such as those by both the Solid State Imager (SSI) and Near-Infrared Mapping Spectrometer (NIMS) instruments on board the Galileo spacecraft. Unfortunately, complete near-infrared spectra, such as those measured by NIMS, take a long time to be analyzed with multiple-scattering radiative transfer models and thus it has until now been rather difficult to use these data to produce wide-area cloud maps.In this paper we show how principal component analysis may be employed to isolate a small number of empirical orthogonal functions (EOFs) from spectra of Jupiter made by Galileo/NIMS. These EOFs may be used to represent the variance of real NIMS spectra to a high degree of accuracy and with good noise and “drop-out” discrimination. Because of this, a small set of representative spectra may then be calculated using these EOFs and input into a retrieval model that generates a table of fitted cloud profiles for each case. This approach avoids the long times required for analyzing a large number of spectra with full multiple-scattering radiative transfer models and allows us to represent the variability of the 3000 spectra contained in the observations of the North Equatorial Belt (NEB) used in this study with only 75 representative spectra. The cloud structures fitted to these representative spectra were interpolated for the spectra found at individual locations in the measured NIMS data set to produce maps of cloud opacity and mean particle size. We find that the dominant opacity variation, anticorrelated with 5-μm brightness, exists in the 1-2 bar pressure range. The distribution of the cloud at 0.72 bars is mapped and found to be more zonally diffuse than the lower clouds. We find at least one 2000-km-sized deep convective cloud in the NEB vertically extending to all the pressure levels that can be sensed by NIMS.     

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.     

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.     

Time-variable emission from transiently accreting neutron stars in quiescence due to deep crustal heating

Transiently accreting neutron stars in quiescence ( L _X ≲10³⁴ erg s⁻¹) have been observed to vary in intensity by factors of few, over time-scales of days to years. If the quiescent luminosity is powered by a hot neutron star core, the core cooling time-scale is much longer than the recurrence time, and cannot explain the observed, more rapid variability. However, the non-equilibrium reactions which occur in the crust during outbursts deposit energy in isodensity shells, from which the thermal diffusion time-scale to the photosphere is days to years. The predicted magnitude of variability is too low to explain the observed variability unless – as is widely believed – the neutrons beyond the neutron-drip density are superfluid. Even then, the variability due to this mechanism in models with standard core neutrino cooling processes is less than 50 per cent – still too low to explain the reported variability. However, models with rapid core neutrino cooling can produce a variability by a factor as great as 20, on time-scales of days to years following an outburst. Thus, the factors of ∼ few intensity variability observed from transiently accreting neutron stars can be accounted for by this mechanism only if rapid core cooling processes are active.