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

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

Impact of electron chemistry on the structure and composition of Io's atmosphere

Two-dimensional model calculations (altitude and solar zenith angle) are performed to investigate the impact of electron chemistry on the composition and structure of Io's atmosphere. The calculations are based upon the model of Wong and Smyth (2000, Icarus 146, 60-74) for Io's SO₂ sublimation atmosphere with the addition of new electron chemistry, where the interactions of the electrons and neutrals are treated in a simple fashion. The model calculations are presented for Io's atmosphere at western elongation (dusk ansa) for both a low-density case (subsolar temperature of 113 K) and a high-density case (subsolar temperature of 120 K). The impact of electron-neutral chemistry on the composition and structure of Io's atmosphere is confined primarily to an interaction layer. The penetration depth of the interaction layer is limited to high altitudes in the thicker dayside atmosphere but reaches the surface in the thinner dayside and/or nightside atmosphere at larger solar zenith angles. Within most of the thicker dayside atmosphere, the column density of SO₂ is not significantly altered by electrons, but in the interaction layer all number densities are significantly altered: SO₂ is reduced, O, SO, S, and O₂ are greatly enhanced, and O, SO, and S become comparable to SO₂ at high altitudes. For the thinner nightside atmosphere, the species number densities are dramatically altered: SO₂ is drastically reduced to the least abundant species of the SO₂ family, SO and O₂ are significantly reduced at all altitudes, and O and S are dramatically enhanced and become the dominant species at all altitudes except near the surface. The interaction layer also defines the location of the emission layer for neutrals excited by electron impact and hence determines the fraction of the total neutral column density that is visible in remote observation. Electron chemistry may also impact the ratio of the equatorial to polar SO₂ column density deduced from Lyman-α images and the north-south alternating and System III longitude-dependent asymmetry observed in polar O and S emissions.     

Temporal behavior of the SO 1.707 μm ro-vibronic emission band in Io's atmosphere

We report observations of the ro-vibronic a¹Δ→X³Σ⁻ transition of SO at 1.707 μm on Io. These data were taken while Io was eclipsed by Jupiter, on four nights between July 2000 and March 2003. We analyze these results in conjunction with a previously published night to investigate the temporal behavior of these emissions. The observations were all conducted using the near-infrared spectrometer NIRSPEC on the W.M. Keck II telescope. The integrated emitted intensity for this band varies from 0.8×¹⁰²⁷ to 2.4×¹⁰²⁷ photons/s, with a possible link to variations in Loki's infrared brightness. The band-shapes imply rotational temperatures of 550-1000 K for the emitting gas, lending further evidence to a volcanic origin for sulfur monoxide. An attempt to detect the B¹Σ→X³Σ⁻ transition of SO at 0.97 μm was unsuccessful; simultaneous detection with the 1.707 μm band would permit determination of the SO column abundance.     

Short-term variability of Jupiter’s extended sodium nebula

Ground-based optical observations of D₁ and D₂ line emissions from Jupiter’s sodium nebula, which extend over several hundreds of jovian radii, were carried out at Mt. Haleakala, Maui, Hawaii using a wide field filter imager from May 19 to June 21, 2007. During this observation, the east-west asymmetry of the nebula with respect to the Io’s orbital motion was clearly identified. Particularly, the D₁+D₂ brightness on the western side of Jupiter is strongly controlled by the Io phase angle. The following scenario was developed to explain this phenomenon as follows: First, more ionospheric ions like NaX⁺, which are thought to produce fast neutral sodium atoms due to a dissociative recombination process, are expected to exist in Io’s dayside hemisphere rather than in the nightside one. Second, it is expected that more NaX⁺ ionospheric ions are picked up by the jovian co-rotating magnetic field when Io’s leading hemisphere is illuminated by the Sun. Third, the sodium atom ejection rate varies with respect to Io’s orbital position as a result of the first two points. Model simulations were performed using this scenario. The model results were consistent with the observation results, suggesting that Io’s ionosphere is expected to be controlled by solar radiation just like Earth.     

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.     

Jupiter’s South South Temperate Zone vortices: Observations and simulations

The region in Jupiter’s atmosphere with the highest density of anticyclonic spot-like vortices is the region known as the South South Temperate Zone (SSTZ), which is located between the eastward jet at ≈−42.9° latitude and the westward jet at ≈−39.2° latitude. We present a characterization of the spots found in this region based on ground-based and Hubble Space Telescope observations from the years 1993 to 2007. Mergers have been reported between spots in this region, similar to those observed for the White Ovals in the latitudinal domain immediately equatorward (northward). We use a multilayer model to perform numerical simulations that capture the details of a well-observed merger event involving multiple interacting vortices. We find that the vertical stratification has an important effect in the outcome of the interaction between spots. In particular it can play a determining role on whether or not a cyclone embedded between two approaching anticyclones can inhibit their merging. From our simulations we conclude that the background static stability of the atmosphere in the SSTZ is better characterized by an average value of .     

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.     

Numerical modeling of ionian volcanic plumes with entrained particulates

Volcanic plumes on Jupiter's moon Io are modeled using the direct simulation Monte Carlo (DSMC) method. The modeled volcanic vent is interpreted as a “virtual” vent. A parametric study of the “virtual” vent gas temperature and velocity is performed to constrain the gas properties at the vent by observables, particularly the plume height and the surrounding condensate deposition ring radius. Also, the flow of refractory nano-size particulates entrained in the gas is modeled with “overlay” techniques which assume that the background gas flow is not altered by the particulates. The column density along the tangential line-of-sight and the shadow cast by the plume are calculated and compared with Voyager and Galileo images. The parametric study indicates that it is possible to obtain a unique solution for the vent temperature and velocity for a large plume like Pele. However, for a small Prometheus-type plume, several different possible combinations of vent temperature and velocity result in both the same shock height and peak deposition ring radius. Pele and Prometheus plume particulates are examined in detail. Encouraging matches with observations are obtained for each plume by varying both the gas and particle parameters. The calculated tangential gas column density of Pele agrees with that obtained from HST observations. An upper limit on the size of particles that track the gas flow well is found to be ∼10 nm, consistent with Voyager observations of Loki. While it is certainly possible for the plumes to contain refractory dust or pyroclastic particles, especially in the vent vicinity, we find that the conditions are favorable for SO₂ condensation into particles away from the vent vicinity for Prometheus. The shadow cast by Prometheus as seen in Galileo images is also reproduced by our simulation. A time averaged frost deposition profile is calculated for Prometheus in an effort to explain the multiple ring structure observed around the source region. However, this multiple ring structure may be better explained by the calculated deposition of entrained particles. The possibility of forming a dust cloud on Io is examined and, based on a lack of any such observed clouds, a subsolar frost temperature of less than 118 K is suggested.     

Thermal structure and composition of Jupiter’s Great Red Spot from high-resolution thermal imaging

Thermal-IR imaging from space-borne and ground-based observatories was used to investigate the temperature, composition and aerosol structure of Jupiter’s Great Red Spot (GRS) and its temporal variability between 1995 and 2008. An elliptical warm core, extending over 8° of longitude and 3° of latitude, was observed within the cold anticyclonic vortex at 21°S. The warm airmass is co-located with the deepest red coloration of the GRS interior. The maximum contrast between the core and the coldest regions of the GRS was 3.0-3.5 K in the north-south direction at 400 mbar atmospheric pressure, although the warmer temperatures are present throughout the 150-500 mbar range. The resulting thermal gradients cause counter-rotating flow in the GRS center to decay with altitude into the lower stratosphere. The elliptical warm airmass was too small to be observed in IRTF imaging prior to 2006, but was present throughout the 2006-2008 period in VLT, Subaru and Gemini imaging.Spatially-resolved maps of mid-IR tropospheric aerosol opacity revealed a well-defined lane of depleted aerosols around the GRS periphery, and a correlation with visibly-dark jovian clouds and bright 4.8-μm emission. Ammonia showed a similar but broader ring of depletion encircling the GRS. This narrow lane of subsidence keeps red aerosols physically separate from white aerosols external to the GRS. The visibility of the 4.8-μm bright periphery varies with the mid-IR aerosol opacity of the upper troposphere. Compositional maps of ammonia, phosphine and para-H₂ within the GRS interior all exhibit north-south asymmetries, with evidence for higher concentrations north of the warm central core and the strongest depletions in a symmetric arc near the southern periphery. Small-scale enhancements in temperature, NH₃ and aerosol opacity associated with localized convection are observed within the generally-warm and aerosol-free South Equatorial Belt (SEB) northwest of the GRS. The extent of 4.8-μm emission from the SEB varied as a part of the 2007 ‘global upheaval,’ though changes during this period were restricted to pressures greater than 500 mbar. Finally, a region of enhanced temperatures extended southwest of the GRS during the survey, restricted to the 100-400 mbar range and with no counterpart in visible imaging or compositional mapping. The warm airmass was perturbed by frequent encounters with the cold airmass of Oval BA, but no internal thermal or compositional effects were noted in either vortex during the close encounters.     

The evolving flow of Jupiter’s White Ovals and adjacent cyclones

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

Titan’s atomic hydrogen corona

Based on measurements performed by the Hydrogen Deuterium Absorption Cell (HDAC) aboard the Cassini orbiter, Titan’s atomic hydrogen exosphere is investigated. Data obtained during the T9 encounter are used to infer the distribution of atomic hydrogen throughout Titan’s exosphere, as well as the exospheric temperature.The measurements performed during the flyby are modeled by performing Monte Carlo radiative transfer calculations of solar Lyman-α radiation, which is resonantly scattered on atomic hydrogen in Titan’s exosphere. Two different atomic hydrogen distribution models are applied to determine the best fitting density profile. One model is a static model that uses the Chamberlain formalism to calculate the distribution of atomic hydrogen throughout the exosphere, whereas the second model is a Particle model, which can also be applied to non-Maxwellian velocity distributions.The density distributions provided by both models are able to fit the measurements although both models differ at the exobase: best fitting exobase atomic hydrogen densities of n_H = (1.5 ± 0.5) × 10⁴ cm⁻³ and n_H = (7 ± 1) × 10⁴ cm⁻³ were found using the density distribution provided by both models, respectively. This is based on the fact that during the encounter, HDAC was sensitive to altitudes above about 3000 km, hence well above the exobase at about 1500 km. Above 3000 km, both models produce densities which are comparable, when taking into account the measurement uncertainty.The inferred exobase density using the Chamberlain profile is a factor of about 2.6 lower than the density obtained from Voyager 1 measurements and much lower than the values inferred from current photochemical models. However, when taking into account the higher solar activity during the Voyager flyby, this is consistent with the Voyager measurements. When using the density profile provided by the particle model, the best fitting exobase density is in perfect agreement with the densities inferred by current photochemical models.Furthermore, a best fitting exospheric temperature of atomic hydrogen in the range of T_H = (150-175) ± 25 K was obtained when assuming an isothermal exosphere for the calculations. The required exospheric temperature depends on the density distribution chosen. This result is within the temperature range determined by different instruments aboard Cassini. The inferred temperature is close to the critical temperature for atomic hydrogen, above which it can escape hydrodynamically after it diffused through the heavier background gas.     

The atmospheric signature of Io's Prometheus plume and anti-jovian hemisphere: evidence for a sublimation atmosphere

Using the Hubble Space Telescope's Space Telescope Imaging Spectrograph we have obtained for the first time spatially resolved 2000-3000 Å spectra of Io's Prometheus plume and adjoining regions on Io's anti-jovian hemisphere in the latitude range 60° N-60° S, using a 0.1″ slit centered on Prometheus and tilted roughly 45° to the spin axis. The SO₂ column density peaked at 1.25×10¹⁷ cm⁻² near the equator, with an additional 5×10¹⁶ cm⁻² enhancement over Prometheus corresponding to a model volcanic SO₂ output of 10⁵ kg s⁻¹. Apart from the Prometheus peak, the SO₂ column density dropped fairly smoothly away from the subsolar point, even over regions that included potential volcanic sources. At latitudes less than ±30°, the dropoff rate was consistent with control by vapor pressure equilibrium with surface frost with subsolar temperature 117.3±0.6 K, though SO₂ abundance was higher than predicted by vapor pressure control at mid-latitudes, especially in the northern hemisphere. We conclude that, at least at low latitudes on the anti-jovian hemisphere where there are extensive deposits of optically-thick SO₂ frost, the atmosphere is probably primarily supported by sublimation of surface frost. Although the 45° tilt of our slit prevents us from separating the dependence of atmospheric density on solar zenith angle from its dependence on latitude, the pattern is consistent with a sublimation atmosphere regardless of which parameter is the dominant control. The observed drop in gas abundance towards higher latitudes is consistent with the interpretation of previous Lyman alpha images of Io as indicating an atmosphere concentrated at low latitudes. Comparison with previous disk-resolved UV spectroscopy, Lyman-alpha images, and mid-infrared spectroscopy suggests that Io's atmosphere is denser and more widespread on the anti-jovian hemisphere than at other longitudes. SO₂ gas temperatures were in the range of 150-250 K over the majority of the anti-jovian hemisphere, consistent with previous observations. SO was not definitively detected in our spectra, with upper limits to the SO/SO₂ ratio in the range 1-10%, roughly consistent with previous observations. S₂ gas was not seen anywhere, with an upper limit of 7.5×10¹⁴ cm⁻² for the Prometheus plume, confirming that this plume is significantly poorer in S₂ than the Pele plume (S₂ /SO₂<0.005, compared to 0.08-0.3 at Pele). In addition to the gas absorption signatures, we have observed continuum emission in the near ultraviolet (near 2800 Å) for the first time. The brightness of the observed emission was directly correlated with the SO₂ abundance, strongly peaking in the equatorial region over Prometheus. Emission brightness was modestly anti-correlated with the jovian magnetic latitude, decreasing when Io intersected the torus centrifugal equator.     

Power spectral analysis of Jupiter’s clouds and kinetic energy from Cassini

We present suggestive evidence for an inverse energy cascade within Jupiter’s atmosphere through a calculation of the power spectrum of its kinetic energy and its cloud patterns. Using Cassini observations, we composed full-longitudinal mosaics of Jupiter’s atmosphere at several wavelengths. We also utilized image pairs derived from these observations to generate full-longitudinal maps of wind vectors and atmospheric kinetic energy within Jupiter’s troposphere. We computed power spectra of the image mosaics and kinetic energy maps using spherical harmonic analysis. Power spectra of Jupiter’s cloud patterns imaged at certain wavelengths resemble theoretical spectra of two-dimensional turbulence, with power-law slopes near −5/3 and −3 at low and high wavenumbers, respectively. The slopes of the kinetic energy power spectrum are also near −5/3 at low wavenumbers. At high wavenumbers, however, the spectral slopes are relatively flatter than the theoretical prediction of −3. In addition, the image mosaic and kinetic energy power spectra differ with respect to the location of the transition in slopes. The transition in slope is near planetary wavenumber 70 for the kinetic energy spectra, but is typically above 200 for the image mosaic spectra. Our results also show the importance of calculating spectral slopes from full 2D velocity maps rather than 1D zonal mean velocity profiles, since at large wavenumbers the spectra differ significantly, though at low wavenumbers, the 1D zonal and full 2D kinetic energy spectra are practically indistinguishable. Furthermore, the difference between the image and kinetic energy spectra suggests some caution in the interpretation of power spectrum results solely from image mosaics and its significance for the underlying dynamics. Finally, we also report prominent variations in kinetic energy within the equatorial jet stream that appear to be associated with the 5 μm hotspots. Other eddies are present within the flow collar of the Great Red Spot, suggesting caution when interpreting snapshots of the flow inside these features as representative of a time-averaged state.     

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.     

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.     

Nature of the iogenic plasma source in Jupiter's magnetosphere: I. Circumplanetary distribution

Three-dimensional calculations are presented for the circumplanetary nature of the iogenic plasma source (pickup ions produced by electron and charge exchange processes in the plasma torus) created by O and S gases located above Io's exobase in its corona and escaping extended neutral clouds (designated as the “Outer Region”). These calculations are undertaken using neutral cloud models for O and S with realistic incomplete collisional cascade source velocity distributions and rates at Io's exobase and realistic spacetime loss processes in the plasma torus. The resulting spatial distributions for O and S about Jupiter are highly peaked at Io but extend at much lower density levels all about the planet, particularly within Io's orbit where they may play a role in the pitch angle scattering and energy loss of radially inward diffusing energetic electrons for the synchrotron radiation belts of Jupiter, in producing bite-outs in the energy distribution of energetic heavy ions near Io's orbit, and in providing a charge exchange source for energetic neutral atoms (ENAs) detected both near and far from Jupiter. For the iogenic plasma source created by these neutrals, two-dimensional distributions produced by integrating the three-dimensional information along the magnetic field lines are presented for the instantaneous values of the pickup ion rates, the total- and net-mass loading rates, the mass-per-unit-magnetic-flux source rate, the pickup conductivity, the pickup radial current, and the pickup ion power (or energy rate). On the circumplanetary spatial scale, the instantaneous iogenic plasma source is highly peaked about Io's position on its orbit around Jupiter. The degree of orbital asymmetry and its physical origin are discussed, and overall spatially integrated rates are presented. The spatially integrated net-mass loading rate is 154 kg s⁻¹ and the total (electron impact and charge exchange) mass loading rate is 275 kg s⁻¹. Rough minimum estimates are made for the spatially integrated total-mass loading rate created by the “Inner Region” (spatial region below Io's exobase) and are at least ∼1 to 2.5 times larger than that for the Outer Region. Implications of the iogenic plasma source created by the Outer Region and the Inner Region are discussed.     

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.     

Nitrile gas chemistry in Titan’s atmosphere

This work presents the first study of the gaseous products resulting from the partial dissociation of methane and nitrogen in the PAMPRE experimental setup simulating Titan’s atmospheric chemistry.Using cryogenic trapping, the gaseous products generated from the chemical reactions occurring in the reactor have been trapped. Analyses of these products by gas chromatography coupled to mass spectrometry have allowed the detection and identification of more than 30 reaction products. Most of them are identified as nitrile species, accompanied by aliphatic hydrocarbons and a few aromatics compounds. The observed species are in agreement with the data from the recent Cassini-Huygens mission as well as from other laboratory setups capable of dissociating nitrogen and methane. This work emphasizes the probable importance of nitrogen-bearing compounds in the chemistry taking place in Titan’s atmosphere.Furthermore, a quantification of mono-nitriles with saturated alkyl chains has been performed relatively to hydrogen cyanide and shows a power law dependence in their concentration. This dependence is consistent with the Cassini-INMS data and Titan’s photochemical models.An empirical relationship has been extracted from our experimental data: [C_xH_2x−1N] = 100x⁻⁵, where x is the number of carbon atoms in the nitrile molecule. This relationship can be directly used in order to foretell the concentration of heavier nitriles induced by chemistry in Titan’s atmosphere.