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

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.     

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⁻².     

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.     

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

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

The ion mass loading rate at Io

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

1-D DSMC simulation of Io's atmospheric collapse and reformation during and after eclipse

A one-dimensional Direct Simulation Monte Carlo (DSMC) model is used to examine the effects of a non-condensable species on Io's sulfur dioxide sublimation atmosphere during eclipse and just after egress. Since the vapor pressure of SO₂ is extremely sensitive to temperature, the frost-supported dayside sublimation atmosphere had generally been expected to collapse during eclipse as the surface temperature dropped. For a pure SO₂ atmosphere, however, it was found that during the first 10 min of eclipse, essentially no change in the atmospheric properties occurs at altitudes above ∼100 km due to the finite ballistic/acoustic time. Hence immediately after ingress the auroral emission morphology above 100 km should resemble that of the immediate pre-eclipse state. Furthermore, the collapse dynamics are found to be greatly altered by the presence of even a small amount of a non-condensable species which forms a diffusion layer near the surface that prevents rapid collapse. It is found that after 10 min essentially no collapse has occurred at altitudes above ∼20 km when a nominal mole fraction of non-condensable gas is present. Collapse near the surface occurs relatively quickly until a static diffusion layer many mean free paths thick of the non-condensable gas builds up which then retards further collapse of the SO₂ atmosphere. For example, for an initial surface temperature of 110 K and 35% non-condensable mole-fraction, the ratio of the SO₂ column density to the initial column density was found to be 0.73 after 10 min, 0.50 after 30 min, and 0.18 at the end of eclipse. However, real gas species (SO, O₂) may not be perfectly non-condensable at Io's surface temperatures. If the gas species was even weakly condensable (non-zero sticking/reaction coefficient) then the effect of the diffusion layer on the dynamics was dramatically reduced. In fact, if the sticking coefficient of the non-condensable exceeds ∼0.25, the collapse dynamics are effectively the same as if there were no non-condensable present. This sensitivity results because the loss of non-condensable to the surface reduces the effective diffusion layer size, and the formation of an effective diffusion layer requires that the layer be stationary; this does not occur if the surface is a sink. Upon egress, vertical stratification of the condensable and non-condensable species occurs, with the non-condensable species being lifted (or pushed) to higher altitudes by the sublimating SO₂ after the sublimating atmosphere becomes collisional. Stratification should affect the morphology and intensity of auroral glows shortly after egress.     

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

Io's sodium emission cloud

Strong evidence that Io's sodium emission is due to resonant scattering is given by our observations which show a monotonic increase of emission intensity with residual solar intensity. In addition we detected no emission during three eclipse observations of Io. We propose a resonant scattering model with two spacial components comprising an optically thick atmosphere extending 10³ km above Io's surface surrounded by an optically thin cloud which forms a partial torus around Jupiter. In this model a flux of 10⁷ cm^?2 sec^?1 sodium atoms are sputtered from Io's surface by heavy energetic ions which are accelerated in a plasma sheath around Io. The atoms sputtered from the surface collide with atoms in Io's atmosphere so the equipartition of kinetic energy is established. The total sodium abundance is about 3 × 10¹³ cm^?2. During Io's day, sodium and other atmospheric constituents are ionized, giving rise to the ionosphere observed by Pioneer 10. Atoms escape by means of Jeans escape from the critical level, which is at the top of the atmosphere and the base of the cloud. We have observed sodium emission 6arcsec (6 Io diameters) above and below Io's orbital plane and 23arcsec toward Jupiter in Io's orbital plane. No emission was detected at maximum elongation 180° from Io. We interpret these results to mean that atoms escaping from Io form a partial torus whose thickness is about 12 arcsec and whose length is at least one-fifth of Io's orbital circumference.     

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.     

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.     

Hubble observations of Jupiter’s north–south conjugate ultraviolet aurora

Comparisons of the northern and southern far ultraviolet (UV) auroral emissions of Jupiter from the Hubble Space Telescope (HST) or any other ultraviolet imager have mostly been made so far on a statistical basis or were not obtained with high sensitivity and resolution. Such observations are important to discriminate between different mechanisms responsible for the electron acceleration of the different components of the aurora such as the satellite footprints, the «main oval» or the polar emissions. The field of view of the ACS and STIS cameras on board HST is not wide enough to provide images of the full jovian disk. We thus compare the morphology of the north and south aurora observed 55 min apart and we point out similarities and differences. On one occasion HST pointed successively the two polar regions and auroral images were seen separated by only 3 min. This makes it possible to compare the emission structure and the emitted FUV power of corresponding regions. We find that most morphological features identified in one hemisphere have a conjugate counterpart in the other hemisphere. However, the power associated with conjugate regions of the main oval, diffuse or discrete equatoward emission observed quasi-simultaneously may be different in the two hemispheres. It is not directly nor inversely proportional to the strength of the B-field as one might expect for diffuse precipitation or field-aligned acceleration with equal ionospheric electron density in both hemispheres. Finally, the lack of symmetry of some polar emissions suggests that some of them could be located on open magnetic field lines.     

X-ray probes of magnetospheric interactions with Jupiter's auroral zones, the Galilean satellites, and the Io plasma torus

Remote observations with the Chandra X-ray Observatory and the XMM-Newton Observatory have shown that the jovian system is a source of X-rays with a rich and complicated structure. The planet's polar auroral zones and its disk are both powerful sources of X-ray emission. Chandra observations revealed X-ray emission from the Io plasma torus and from the Galilean moons Io, Europa, and possibly Ganymede. The emission from the moons is due to bombardment of their surfaces by highly energetic magnetospheric protons, and oxygen and sulfur ions. These ions excite atoms in their surfaces leading to fluorescent X-ray emission lines. These lines are produced against an intense background continuum, including bremsstrahlung radiation from surface interactions of primary magnetospheric and secondary electrons. Although the X-ray emission from the Galilean moons is faint when observed from Earth orbit, an imaging X-ray spectrometer in orbit around one or more of these moons, operating from 200 eV to 8 keV with 150 eV energy resolution, would provide a detailed mapping of the elemental composition in their surfaces. Surface resolution of 40 m for small features could be achieved in a 100-km orbit around one moon while also remotely imaging surfaces of other moons and Jupiter's upper atmosphere at maximum regional resolutions of hundreds of kilometers. Due to its relatively more benign magnetospheric radiation environment, its intrinsic interest as the largest moon in the Solar System, and its mini-magnetosphere, Ganymede would be the ideal orbital location for long-term observational studies of the jovian system. Here we describe the physical processes leading to X-ray emission from the surfaces of Jupiter's moons and the properties required for the technique of imaging X-ray spectroscopy to map the elemental composition of their surfaces, as well as studies of the X-ray emission from the planet's aurora and disk and from the Io plasma torus.     

Diffuse auroral precipitation in the jovian upper atmosphere and magnetospheric electron flux variability

The deposition of energetic electrons in Jupiter's upper atmosphere provides a means, via auroral observations, of monitoring electron and plasma wave activity within the magnetosphere. Not only does particle precipitation indicate a potential change in atmospheric chemistry, it allows for the study of episodic, pronounced flux enhancements in the energetic electron population. A study has been made of the effects of such electron injections into the jovian magnetosphere and of their ability to provide the source population for variations in diffuse auroral emissions. To identify the source region of precipitating auroral electrons, we have investigated the pitch-angle distributions of high-resolution Galileo Energetic Particle Detector (EPD) data that indicate strong flux levels near the loss cone. The equatorial source region of precipitating electrons has been determined from the locations of Galileo's in situ measurements by tracing magnetic field lines using the KK97 model. The primary source region for Jupiter's diffuse aurora appears to lie in the magnetic equator at 15-40 RJ, with the predominant contribution to precipitation flux (tens of ergs cm⁻² s⁻¹ sr⁻¹) stemming from <30 RJ. Variability of flux for energetic electrons in this region is also important to the irradiation of surfaces and atmospheres for the Galilean moons: Europa, Ganymede, and Callisto. The average diffuse auroral precipitation flux has been shown to vary by as much as a factor of six at a given radial location. This variability appears to be associated with electron injection events that have been identified in high-resolution Galileo EPD data. These electron flux enhancements are also associated with increased whistler-mode wave activity and magnetic field perturbations, as detected by the Galileo Plasma Wave Subsystem (PWS) and Magnetometer (MAG), respectively. Resonant interactions with the whistler-mode waves cause electron pitch-angle scattering and lead to pitch-angle isotropization and precipitation.     

Simultaneous mapping of SO2, SO, NaCl in Io’s atmosphere with the Submillimeter Array

Many of the key properties of Io’s atmosphere, such as its spatial distribution, temperature, column density and composition, are still not fully assessed despite decades of extensive observations. The contribution of the possible gas sources to the atmospheric replenishment are then still unclear.This paper presents disk-resolved observations performed with the Submillimeter Array (SMA) at 345 GHz of atmospheric rotational lines of the main atmospheric species SO₂, and, for the first time, of the minor species SO and NaCl. All these species appear concentrated on the anti-jovian hemisphere, but do not share the same spatial distribution. The obtained maps and line-averaged fluxes are compared to realistic models simulating gas sources including volcanic plume outgassing, SO₂ frost sublimation and photolysis. Arguments in favor of each sources are examined and compared to observations, putting constraints on their relative roles for each species.While sublimation clearly appears as the favored major source for SO₂, SO₂ photolysis may account for most of the production of SO. Using constraints on the volcanic plumes distribution from Galileo results, we find that direct volcanic input can only contribute for a minor fraction of atmospheric SO₂, but represent a more significant source for SO atmosphere, and is likely to be the only source for NaCl. Temperature and column densities findings are also presented for SO₂, and compare well to previously published observations and atmospheric models.     

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.     

NN detectability in Pluto’s atmosphere

Based on the vapor pressure behavior of Pluto’s surface ices, Pluto’s atmosphere is expected to be predominantly composed of N₂ gas. Measurement of the N₂ isotopologue ¹⁵N/¹⁴N ratio within Pluto’s atmosphere would provide important clues to the evolution of Pluto’s atmosphere from the time of formation to its present state. The most straightforward way of determining the N₂ isotopologue ¹⁵N/¹⁴N ratio in Pluto’s atmosphere is via spectroscopic observation of the ¹⁴N¹⁵N gas species. Recent calculations of the 80–100 nm absorption behavior of the ¹⁴N₂ and ¹⁴N¹⁵N isotopologues by Heays et al. (Heays, A.N. et al. [2011]. J. Chem. Phys. 135, 244301), Lewis et al. (Lewis, B.R., Heays, A.N., Gibson, S.T., Lefebvre-Brion, H., Lefebvre, R. [2008]. J. Chem. Phys. 129, 164306); Lewis et al. (Lewis, B.R., Gibson, S.T., Zhang, W., Lefebvre-Brion, H., Robbe, J.-M. [2005]. J. Chem. Phys. 122, 144302), and Haverd et al. (Haverd, V.E., Lewis, B.R., Gibson, S.T., Stark, G. [2005]. J. Chem. Phys. 123, 214304) show that the peak magnitudes of the ¹⁴N₂ and ¹⁴N¹⁵N absorption bandhead cross-sections are similar, but the locations of the bandhead peaks are offset in wavelength by ∼0.05–0.1 nm. These offsets make the segregation of the ¹⁴N₂ and ¹⁴N¹⁵N absorption signatures possible. We use the most recent N₂ isotopologue absorption cross-section calculations and the atmospheric density profiles resulting from photochemical models developed by Krasnopolsky and Cruickshank (Krasnopolsky, V.A., Cruickshank, D.P. [1999]. J. Geophys. Res. 104, 21979–21996) to predict the level of solar light that will be transmitted through Pluto’s atmosphere as a function of altitude during a Pluto solar occultation. We characterize the detectability of the isotopic absorption signature per altitude assuming ¹⁴N¹⁵N concentrations ranging from 0.1% to 2% of the ¹⁴N₂ density and instrumental spectral resolutions ranging from 0.01 to 0.3 nm. Our simulations indicate that optical depth of unity is attained in the key ¹⁴N¹⁵N absorption bands located between 85 and 90 nm at altitudes ∼1100–1600 km above Pluto’s surface. Additionally, an ¹⁴N¹⁵N isotope absorption depth ∼4–15% is predicted for observations obtained at these altitudes at a spectral resolution of ∼0.2–0.3 nm, if the N₂ isotopologue ¹⁵N/¹⁴N percent ratio is comparable to the 0.37–0.6% ratio observed at Earth, Titan and Mars. If we presume that the predicted absorption depth must be at least 25% greater than the expected observational uncertainty, then it follows that a statistically significant detection of these signatures and constraint of the N₂ isotopologue ¹⁴N/¹⁵N ratio within Pluto’s atmosphere will be possible if the attainable observational signal-to noise (S/N) ratio is ?9. The New Horizons (NH) Mission will be able to obtain high S/N, 0.27–0.35 nm full-width half-max 80–100 nm spectral observations of Pluto using the Alice spectrograph. Based on the NH/Alice specifications we have simulated 0.3 nm spectral resolution solar occultation spectra for the 1100–1600 km altitude range, assuming 30 s integration times. These simulations indicate that NH/Alice will obtain spectral observations within this altitude range with a S/N ratio ∼25–50, and should be able to reliably detect the ¹⁴N¹⁵N gas absorption signature between 85 and 90 nm if the ¹⁴N¹⁵N concentration is ∼0.3% or greater. This, additionally, implies that the non-detection of the ¹⁴N¹⁵N species in the 1100–1600 km range by NH/Alice may be used to reliably establish an upper limit to the N₂ isotopologue ¹⁵N/¹⁴N ratio within Pluto’s atmosphere. Similar results may be derived from 0.2 to 0.3 nm spectral resolution observations of any other N₂-rich Solar System or exoplanet atmosphere, provided the observations are attained with similar S/N levels.     

Electron precipitation equatorward of the auroral oval and the mantle aurora in the midday sector

In the midday sector, the hard electron precipitation and the associated patchy aurora at geomagnetic latitude ～65° are the only auroral features (? 20 keV) located equatorward of the dayside auroral oval during intense and moderately disturbed geomagnetic conditions. We identify the patchy luminosity in the midday and late morning sectors as the active mantle aurora. The mantle aurora was found by Sanford (1964) using the IGY-IGC auroral patrol spectrographs and which was thought to be non-visual. The precipitating electrons reside mostly at energies greater than several keV with an energy flux of ? 0.1 erg cm^?2 s^?1 sr^?1 during geomagnetic active periods. This hard precipitation occurs in a region which is asymmetric in L.T. with respect to the noon meridian. The region extends from the morning sector to only early afternoon (13–14 M.L.T.) along the geomagnetic latitude circle of about 65–70°. The model calculation indicates that the mantle aurora is produced by the precipitation of the energetic electrons which drift azimuthally from the plasma sheet at the midnight sector to the dayside magnetopause during magnetospheric substorms.     

Time-dependent studies of the aurora: Effects of particle precipitation on the dynamic morphology of ionospheric and atmospheric properties

A self-consistent, time-dependent numerical model of the aurora and high-latitude ionos-phere has been developed. It is used to study the response of ionospheric and atmospheric properties in regions subjected to electron bombardment. The time history of precipitation events is arbitrarily specified and computations are made for a variety of electron spectral energy distributions and flux magnitudes. These include soft electron precipitation, such as might occur on the poleward edge of the auroral oval and within the magnetospheric cleft, and harder spectra representative of particle precipitation commonly observed within and on the equatorward edge of the auroral oval. Both daytime and night-time aurorae are considered. The results of the calculations show that the response of various ionospheric and atmospheric parameters depends upon the spectral energy distribution and flux magnitudes of the precipitating electrons during the auroral event. Various properties respond with different time constants that are influenced by coupling processes described by the interactive model. The soft spectrum aurora affects mainly the ionospheric F region, where it causes increases in the electron density, electron temperature and the 6300 Å red line intensity from normal quiet background levels during both daytime and night-time aurora. The fractional variation is greater for the night-time aurora. The hard spectrum aurorae, in general, do not greatly affect the F-2 region of the ionosphere; however, in the F-1 and E regions, large increases from background conditions are shown to occur in the electron and ion temperatures, electron and ion densities, airglow emission rates and minor neutral constituent densities during the build-up phase of the auroral event. During the decay phase of the aurora, most of these properties decrease at nearly the same rate as the specified particle precipitation flux. However, some ionospheric and atmospheric species have a long memory of the auroral event. The odd nitrogen species N(⁴S) and NO probably do not ever reach steady-state densities between auroral storms.     

UV and IR auroral emission model for the outer planets: Jupiter and Saturn comparison

Planetary aurora display the dynamic behavior of the plasma gas surrounding a planet. The outer planetary aurora are most often observed in the ultraviolet (UV) and the infrared (IR) wavelengths. How the emissions in these different wavelengths are connected with the background physical conditions are not yet well understood. Here we investigate the sensitivity of UV and IR emissions to the incident precipitating auroral electrons and the background atmospheric temperature, and compare the results obtained for Jupiter and Saturn. We develop a model which estimates UV and IR emission rates accounting for UV absorption by hydrocarbons, ion chemistry, and non-LTE effects. Parameterization equations are applied to estimate the ionization and excitation profiles in the H₂ atmosphere caused by auroral electron precipitation. The dependences of UV and IR emissions on electron flux are found to be similar at Jupiter and Saturn. However, the dependences of the emissions on electron energy are different at the two planets, especially for low energy (<10 keV) electrons; the UV and IR emissions both decrease with decreasing electron energy, but this effect in the IR is less at Saturn than at Jupiter. The temperature sensitivity of the IR emission is also greater at Saturn than at Jupiter. These dependences are interpreted as results of non-LTE effects on the atmospheric temperature and density profiles. The different dependences of the UV and IR emissions on temperature and electron energy at Saturn may explain the different appearance of polar emissions observed at UV and IR wavelengths, and the differences from those observed at Jupiter. These results lead to the prediction that the differences between the IR and UV aurora at Saturn may be more significant than those at Jupiter. We consider in particular the occurrence of bright polar infrared emissions at Saturn and quantitatively estimate the conditions for such IR-only emissions to appear.