Methane absorption in the atmosphere of Jupiter from 1800 to 9500 cm and implications for vertical cloud structure

New measurements of the low-temperature near-infrared absorption of methane (Sihra, 1998, Laboratory measurements of near-infrared methane bands for remote sensing of the jovian atmosphere, Ph.D. thesis, University of Oxford) have been combined with existing, longer path-length, higher-temperature data of Strong et al. (1993, Spectral parameters of self- and hydrogen-broadened methane from 2000 to 9500 cm⁻¹ for remote sounding of the atmosphere of Jupiter, J. Quant. Spectrosc. Radiat. Trans. 50, 309-325) and fitted with band models. The combined data set is found to be more consistent with previous low-temperature methane absorption measurements than that of Strong et al. (1993, J. Quant. Spectrosc. Radiat. Trans. 50, 309-325) but covers the same wider wavelength range and accounts for both self- and hydrogen-broadening conditions. These data have been fitted with k-coefficients in the manner described by Irwin et al. (1996, Calculated k-distribution coefficients for hydrogen- and self-broadened methane in the range 2000-9500 cm⁻¹ from exponential sum fitting to band modelled spectra, J. Geophys. Res. 101, 26,137-26,154) and have been used in multiple-scattering radiative transfer models to assess their impact on our previous estimates of the jovian cloud structure obtained from Galileo Near-Infrared Mapping Spectrometer (NIMS) observations (Irwin et al., 1998, Cloud structure and atmospheric composition of Jupiter retrieved from Galileo NIMS real-time spectra, J. Geophys. Res. 103, 23,001-23,021; Irwin et al., 2001, The origin of belt/zone contrasts in the atmosphere of Jupiter and their correlation with 5-μm opacity, Icarus 149, 397-415; Irwin and Dyudina, 2002, The retrieval of cloud structure maps in the equatorial region of Jupiter using a principal component analysis of Galileo/NIMS data, Icarus 156, 52-63). Although significant differences in methane opacity are found at cooler temperatures, the difference in the optical depth of the atmosphere due to methane is found to diminish rapidly with increasing pressure and temperature and thus has negligible effect on the cloud structure inferred at deeper levels. Hence the main cloud opacity variation is still found to peak at around 1-2 bar using our previous analytical approach, and is thus still in disagreement with Galileo Solid State Imager (SSI) determinations (Banfield et al., 1998, Jupiter's cloud structure from Galileo imaging data, Icarus 135, 230-250; Simon-Miller et al., 2001, Color and the vertical structure in Jupiter's belts, zones and weather systems, Icarus 154, 459-474) which place the main cloud deck near 0.9 bar. Further analysis of our retrievals reveals that this discrepancy is probably due to the different assumptions of the two analyses. Our retrievals use a smooth vertically extended cloud profile while the SSI determinations assume a thin NH₃ cloud below an extended haze. When the main opacity in our model is similarly assumed to be due to a thin cloud below an extended haze, we find the main level of cloud opacity variation to be near the 1 bar level—close to that determined by SSI and moderately close to the expected condensation level of ammonia ice of 0.85 bar, assuming that the abundance of ammonia on Jupiter is (7±1)×¹⁰⁻⁴ (Folkner et al., 1998, Ammonia abundance in Jupiter's atmosphere derived from the attenuation of the Galileo probe's radio signal, J. Geophys. Res. 103, 22,847-22,855; Atreya et al., 1999, A comparison of the atmospheres of Jupiter and Saturn: deep atmospheric composition, cloud structure, vertical mixing, and origin, Planet. Space Sci. 47, 1243-1262). However our data in the 1-2.5 μm range have good height discrimination and our lowest estimate of the cloud base pressure of 1 bar is still too great to be consistent with the most recent estimates of the ammonia abundance of 3.5 × solar. Furthermore the observed limited spatial distribution of ammonia ice absorption features on Jupiter suggests that pure ammonia ice is only present in regions of localised vigorous uplift (Baines et al., 2002, Fresh ammonia ice clouds in Jupiter: spectroscopic identification, spatial distribution, and dynamical implications, Icarus 159, 74-94) and is subsequently rapidly modified in some way which masks its pure absorption features. Hence we conclude that the main cloud deck on Jupiter is unlikely to be composed of pure ammonia ice and instead find that it must be composed of either NH₄SH or some other unknown combination of ammonia, water, and hydrogen sulphide and exists at pressures of between 1 and 2 bar.     

The Hanle effect of the coronal Lα line of hydrogen: Theoretical investigation

This paper is devoted to a computation of the effect of a magnetic field on the linear polarization of the coronal L line of hydrogen. Recent works (Gabriel et al., 1971) have shown that the linear polarization of this line is due to resonant scattering of the incident chromospheric L line. The Hanle effect is the modification of this linear polarization, due to the magnetic field. After having briefly recalled the main features of this effect and the conditions of the coronal L line formation, we present the theoretical formalism to be used for Hanle effect computations. The effect of the hyperfine structure of the line is included. Then the results of our computations are given in terms of linear polarization as a function of the magnetic field. We get that the effect of the hyperfine structure on these results is negligible, although this is not evident a priori. When the hyperfine structure is neglected, the line structure is simplified and the Hanle effect can be expressed with analytical formulae, which we give in the last part of this paper. After integration along the line of sight, these formulae could be used for magnetic field determination in the solar corona from measurements of the linear polarization of the L line.     

Polarization studies of the Venus uv contrasts: Cloud height and haze variability

Polarimetry is able to show direct evidence for compositional differences in the Venus clouds. We present observations (collected during $\text{2}\text{1}\text{2}$ Venus years by the Pioneer Venus Orbiter) of the polarization in four colors of the bright and dark ultraviolet features. We find that the polarization is significantly different between the bright and dark areas. The data show that the “null” model of L. W. Esposito (1980, J. Geophys. Res.85, 8151–8157) and the “overlying haze” model of J. B. Pollack et al. (1980, J. Geophys. Res.85, 8223–8231) are insufficient. Exact calculations of the polarization, including multiple scattering and vertical inhomogeneity near the Venus cloud tops, are able to match the observations. Our results give a straightforward interpretation of the polarization differences in terms of known constituents of the Venus atmosphere. The submicron haze and uv absorbers are anticorrelated: for haze properties as given by K. Kawabata et al. (1980, J. Geophys. Res.85, 8129–8140) the excess haze depth at 9350 Å over the bright regions is Δτ_h = 0.03 ± 0.02. The cloud top is slightly lower in the dark features: the extra optical depth at 2700 Å in Rayleigh scattering above the darker areas is Δτ_R = 0.010 ± 0.005. This corresponds to a height difference of 1.2 ± 0.6 km at the cloud tops. The calculated polarization which matches our data also explains the relative polarization of bright and dark features observed by Mariner 10. The observed differential polarization cannot be explained by differential distribution of haze, if the haze aerosols have an effective size of 0.49 μm, as determined by K. Kawabata et al. (1982, submitted) for the aerosols overlying the Venus equator. We propose two models for the uv contrasts consistent with our results. In a physical model, the dark uv regions are locations of vertical convergence and horizontal divergence. In a chemical model, we propose that the photochemistry is limited by local variations in water vapor and molecular oxygen. The portions of the atmosphere where these constituents are depleted at the cloud tops are the dark uv features. Strong support for this chemical explanation is the observation that the number of sulfur atoms above the cloud tops is equal over both the bright and dark areas. The mass budget of sulfur at these altitudes is balanced between excess sulfuric acid haze over the bright regions and excess SO₂ in the dark regions.     

Vector magnetic fields in prominences

A preliminary discussion is presented of measurements of the polarization of the He i D₃ multiplet in a quiescent prominence, observed with a wavelength-scanning Stokes polarimeter. For a series of 43 observations in the same prominence, the linear polarization of the major component of D₃ lies primarily in the range 1 to 2% and of the wing component, the range 2 to 5%; the polarization vector angle lies primarily in the range 10–25° for the major component, and 25–35° for the other component. From a more limited data set, the polarization of both components is found to first increase as a function of height in the prominence, and then to decrease; the polarization angles of the major component vary in a random-like way with height, while the wing component shows a systematic change. The amount of polarization and the angle of polarization are governed by the Hanle effect. The collective effect of the group of lines at the peak of D₃ evidently has a different sensitivity to the Hanle effect than does the wing component, thus yielding at least four independent measurements - two polarizations and two angles. With some redundancy, the vector magnetic field can then be established using the detailed theory of the Hanle effect. Since the wing component of D₃ is a simple triplet, an initial estimate of the magnetic field strength and its horizontal orientation, 0, relative to the line of sight, is simply obtained. Examples of such calculations are presented.The National Center for Atmospheric Research is sponsored by the National Science Foundation.Operated by the Association of Universities for Research in Astronomy, Inc. under contract AST 78-17292 with the National Science Foundation.     

Observation of the hydrogen corona with SPICAM on Mars Express

A series of seven dedicated Lyman-α observations of exospheric atomic hydrogen in the martian corona were performed in March 2005 with the ultraviolet spectrometer SPICAM on board Mars Express. Two types of observations are analyzed, observations at high illumination (for a solar zenith angle SZA equal to 30°) and observations at low illumination (for a solar zenith angle equal to 90° (terminator), and near the south pole). The measured Lyman-α emission is interpreted as purely resonant scattering of solar photons. Because the Lyman-α emission is optically thick, we use a forward model to analyze this data. Below the exobase, the hydrogen density is described by a diffusive model and above the exobase, it follows Chamberlain's approach without satellite particles. For different hydrogen density profiles between 80 and 50,000 km, the volume emission rates are computed by solving the radiative transfer equation. Such an approach has been used to analyze the Mariner 6, 7 exospheric Lyman-α data during the late 1960s. A reasonable fit of the set of observations is obtained assuming an exobase temperature between 200 and 250 K and an exobase density of ∼1-4 × 10⁵ cm⁻³ in good agreement with photochemical models. However, when considering the average exospheric temperature of 200 K measured by other methods [Leblanc, F., Chaufray, J.Y., Witasse, O., Lilensten, J., Bertaux, J.-L., 2006a. J. Geophys. Res. 111 (E9), doi:10.1029/2005JE002664. E09S11; Leblanc, F., Chaufray, J.-Y., Bertaux, J.-L., 2007. Geophys. Res. Lett. 34, doi:10.1029/2006GL028437. L02206; Bougher, S.W., Engel, S., Roble, R.G., Foster, B., 2000. J. Geophys. Res. 105, 17669-17692; Mazarico, E., Zuber, M.T., Lemoine, F.G., Smith, D.E., 2007. J. Geophys. Res. 112, doi:10.1029/2006JE002734. E05014] a supplementary hot population is needed above the exobase to reconcile Lyman-α measurements with these previous measurements, particularly for the observations at low SZA. In this case, the densities and temperatures at the exobase for the two populations are 1.0±0.2×¹⁰⁵ cm⁻³ and T=200 K and 1.9±0.5×¹⁰⁴ cm⁻³ and T>500 K for the cold and hot populations respectively at low SZA. At high SZA, the densities and temperatures are equal to 2±0.4×¹⁰⁵ cm⁻³ and T=200 K and n=1.2±0.5×¹⁰⁴ cm⁻³ and T>500 K. These high values of the hot hydrogen component are not presently supported by the theory. Moreover, it is important to underline that the two population model remains relatively poorly constrained by the limited range of altitude covered by the present set of SPICAM measurements and cannot be unambiguously identified because of the global uncertainty of our model and of SPICAM measurements. For a one population solution, an average water escape rate of 7.5 × 10⁻⁴ precipitable μm/yr is estimated, yielding a lifetime of 13,000 years for the average present water vapor content of 10 precipitable microns.     

Circular polarization of the 1.4 GHz emission from Jupiter's radiation belts

Measurements of the 1.4 GHz emission from Jupiter made when D_E was 3°·1 show the circular polarization to vary from +0.8 to ?1.1% as the planet rotates. The rms scatter of the points about the mean curve is only 0.09%. Expressed as a function of Jovian magnetic latitude the polarization at first increases linearly but beyond latitudes ～7° the curve flattens. This shape requires that the radiating electrons have a pitch angle distribution similar to that inferred earlier from the beaming and linear polarization. The magnitude of the circular polarization requires an equatorial magnetic flux density in the belt of about 0.3G, consistent with the Pioneer results.Compared with measurements made one orbital period earlier, the total flux density has decreased by 15%, but the beaming has not changed appreciably.     

Influence of the Atmospheric Model on Hanle Diagnostics

We clarify the uncertainty in the inferred magnetic field vector via the Hanle diagnostics of the hydrogen Lyman-\(\upalpha\) line when the stratification of the underlying atmosphere is unknown. We calculate the anisotropy of the radiation field with plane-parallel semi-empirical models under the nonlocal thermal equilibrium condition and derive linear polarization signals for all possible parameters of magnetic field vectors based on an analytical solution of the atomic polarization and Hanle effect. We find that the semi-empirical models of the inter-network region (FAL-A) and network region (FAL-F) show similar degrees of anisotropy in the radiation field, and this similarity results in an acceptable inversion error (e.g., \({\sim}\, 40~\mbox{G}\) instead of 50 G in field strength and \({\sim}\,100^{\circ}\) instead of \(90^{\circ}\) in inclination) when FAL-A and FAL-F are swapped. However, the semi-empirical models of FAL-C (averaged quiet-Sun model including both inter-network and network regions) and FAL-P (plage regions) yield an atomic polarization that deviates from all other models, which makes it difficult to precisely determine the magnetic field vector if the correct atmospheric model is not known (e.g., the inversion error is much larger than 40% of the field strength; \({>}\,70~\mbox{G}\) instead of \(50~\mbox{G}\)). These results clearly demonstrate that the choice of model atmosphere is important for Hanle diagnostics. As is well known, one way to constrain the average atmospheric stratification is to measure the center-to-limb variation of the linear polarization signals. The dependence of the center-to-limb variations on the atmospheric model is also presented in this paper.     

Local Lyman α emitters and their relevance to high-redshift star-forming galaxies

The Lyα line is an important diagnostic of star formation and activity in galaxies. The analysis of Lyα is complicated due to the resonant nature of the line and radiative transfer effects. High spectral resolution studies of local starburst galaxies with the unprecedented UV capabilities of the HST have shown that this line is either seen in absorption or in emission and in the latter case with a P?Cygni profile indicative of a large scale outflow of neutral gas. Moreover, HST imaging obtained with HST-ACS of a sample of 6 star-forming galaxies has revealed that a substantial fraction of the Lyα photons are diffused far away from the emissive knots. Since the importance of Lyα for tracing large scale structure, correlation functions, and galaxy formation is recognized, Lyα will remain a very important probe of the distant universe for the foreseeable future, and it is therefore imperative to acquire a better understanding of what mechanisms regulate our ability to detect this line.     

Recent observations of the decimetric radio emission from Jupiter

The position angle of linear polarization and the degree of circular polarization of the decimetric emission show no marked variations with time, although there may be an indication in the position angle data of a secular decrease in the inclination of the magnetic field structure. High resolution observations of the radiation belts show that the asymmetry in the emission found in synthesis maps originates from an extended region. The presence of this features accounts well for the observed variation in the position of the centroid of the emission with rotation. The radio observations suggest that any equatorial displacement of the main magnetic field structure from the center of the planet is less than 0.05 R_J.     

The determination of vector magnetic fields in prominences from the observations of the Stokes profiles in the D3 line of helium

A method is presented to measure the magnetic field vector in prominences by means of the polarimetric observations in the D₃ line of He obtained with the High Altitude Observatory Stokes polarimeter. The characteristics of the observed Stokes profiles are discussed. The theory of the Hanle effect is reformulated in the representation of the irreducible tensors of the density matrix, and is generalized to derive the circular polarization profiles across the spectral line in terms of the intensity and direction of the prominence magnetic field. The circular polarization profile so deduced can be employed to obtain useful information which adds to that carried by the linear polarization observations. A non-linear least-squares algorithm is proposed to derive the measurement of the magnetic field from the observations, and a consistency check is suggested to test the adequacy of the theoretical model to describe the physics of the He I atomic excitation in prominences.On leave from: Astrophysical Observatory of Arcetri, Largo E. Fermi, 5, 50125 Firenze, Italy.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Vector magnetic fields in prominences

The Stokes components of He i D₃ emission in two quiescent prominences, using full spectral profile measurements, are analyzed to derive vector magnetic fields. Two independently developed schemes, based on the Hanle effect, are used for interpretation. They involve solutions of the statistical equilibrium equations for the He i D₃ multiplet, including the effect of coherency and full level crossing, which predict the magnetic field dependence of the observed polarization. Derived magnetic field vector solutions for each pair of linear polarization Stokes profiles corresponding to an observational point in the prominence are, intrinsically, not uniquely determined, and a set of possible solutions is usually obtained. However, mutual consistency of these solutions with those independently predicted by the form of the circular polarized component, allow, in almost all cases, rejection of all solutions of a set except one symmetrical pair. Of such a pair, a unique solution can be determined with a high confidence level by reference to independent potential field information. Field vectors are found usually to be close to horizontal and normal to the prominence surface, but extreme exceptions are found. Field values range from 6 G to 60 G. The derived vectorfield configurations and their magnitudes are briefly discussed relative to these prominences and to different quiescent prominence models.The National Center for Atmospheric Research is sponsored by the National Science Foundation.Operated by the Association of Universities for Research in Astronomy, Inc. under contract with the National Science Foundation.     

The Hanle effect and the diagnostics of turbulent magnetic fields in the solar atmosphere

The theory of the Hanle effect is used to interpret the linear polarization measured in a number of spectral lines on the solar disk near the heliographic north and south poles, in search for a turbulent magnetic field in the solar atmosphere. The Hanle depolarization is separated from a number of other effects, including collisional depolarization and scattering geometry. Although the main aim of the paper is to elucidate the physics of the Hanle effect as applied to the Sun, our results indicate the existence of hidden or turbulent magnetic flux near the temperature minimum of the solar atmosphere, with a field strength between 10 and 100 G. This field is hidden in the sense that it is not seen in measurements of the longitudinal Zeeman effect (solar magnetograms). It carries more total magnetic flux than the kG network fields.     

Planetary dynamos: effects of electrically conducting flows overlying turbulent regions of magnetic field generation

A fully three-dimensional, nonlinear, time-dependent, multi-layered spherical kinematic dynamo model is used to study the effect on the observable external magnetic field of flow in an electrically conducting layer above a spherical turbulent dynamo region in which the α effect generates the magnetic field. It is shown that the amplitude and structure of an observable planetary magnetic field are largely determined by the magnitude and structure of the flow in the overlying layer. It is also shown that a strong-field planetary dynamo can be readily produced by the effect of an electrically conducting flow layer at the top of a convective core. The overlying layer and the underlying convective region constitute a magnetically strongly coupled system. Such overlying layers might exist at the top of the Earth's core due to chemical or thermal causes, in the cores of other terrestrial planets for similar reasons, and in Saturn due to the differentiation of helium from hydrogen. An electrically conducting and differentially rotating layer could exist above the metallic hydrogen region in Jupiter and affect the jovian magnetic field similar to the overlying layers in other planets. Lateral temperature gradients resulting in thermal winds drive the flow in the overlying layers. All planetary magnetic fields could be maintained by similar turbulent convective dynamos in the field-generation regions of planets with the differences among observable magnetic fields due to different circulations in the overlying electrically conducting layers.     

Deuterium chemistry and airglow in the jovian thermosphere

We present a detailed study of the distribution of key deuterated species (viz., atomic D and HD) and the associated deuterium Lyman-α airglow in the jovian thermosphere. The reactions that appear to govern the abundances of these deuterated species are used in conjunction with C₂-chemistry in a 1-D photochemical-diffusion model. While the D abundance is mainly sensitive to H densities and the vibrational temperature profile, the D vertical distribution also depends on other parameters such as eddy mixing and the uncertain values of some of the reaction rate constants. We consider different scenarios by varying several parameters controlling the D distribution in the thermosphere. A radiative transfer model with coupling of the H and D Lyman-α lines is employed to obtain line profiles and total intensities at disk center for these scenarios. This allows a comparison of the impact of various parameters on the jovian D Lyman-α emission. A consequence of these chemical processes in the jovian thermosphere is the formation of CH₂D, CH₃D, and C₂H₅D, and other deuterated species. We also discuss the source of these deuterated hydrocarbons and their abundance. We find that HD vibrational chemistry impacts D in the thermosphere, CH₃D and C₂H₅D are vibrationally enhanced in the thermosphere, and variations in abundance of CH₃D and C₂H₅D in the thermosphere may reflect dynamical activity (i.e., Kh) in the jovian upper atmosphere. An observing program dedicated to providing such measurements of these testable phenomena would provide further insight into the synergistic coupling between chemistry, energetics and airglow in the jovian upper atmosphere.     

Updated Galileo probe mass spectrometer measurements of carbon, oxygen, nitrogen, and sulfur on Jupiter

The in situ measurements of the Galileo Probe Mass Spectrometer (GPMS) were expected to constrain the abundances of the cloud-forming condensible volatile gases: H₂O, H₂S, and NH₃. However, since the probe entry site (PES) was an unusually dry meteorological system—a 5-μm hotspot—the measured condensible volatile abundances did not follow the canonical condensation-limited vertical profiles of equilibrium cloud condensation models (ECCMs) such as Weidenschilling and Lewis (1973, Icarus 20, 465-476). Instead, the mixing ratios of H₂S and NH₃ increased with depth, finally reaching well-mixed equilibration levels at pressures far greater than the lifting condensation levels, whereas the mixing ratio of H₂O in the deep well-mixed atmosphere could not be measured. The deep NH₃ mixing ratio (with respect to H₂) of (6.64±2.54)×10⁻⁴ from 8.9-11.7 bar GPMS data is consistent with the NH₃ profile from probe-to-orbiter signal attenuation (Folkner et al., 1998, J. Geophys. Res. 103, 22847-22856), which had an equilibration level of about 8 bar. The GPMS deep atmosphere H₂S mixing ratio of (8.9±2.1)×10⁻⁵ is the only measurement of Jupiter's sulfur abundance, with a PES equilibration level somewhere between 12 and 15.5 bar. The deepest water mixing ratio measurement is (4.9±1.6)×10⁻⁴ (corresponding to only about 30% of the solar abundance) at 17.6-20.9 bar, a value that is probably much smaller than Jupiter's bulk water abundance. The ¹⁵N/¹⁴N ratio in jovian NH₃ was measured at (2.3±0.3)×10⁻³ and may provide the best estimate of the protosolar nitrogen isotopic ratio. The GPMS methane mixing ratio is (2.37±0.57)×10⁻³; although methane does not condense on Jupiter, we include its updated analysis in this report because like the condensible volatiles, it was presumably brought to Jupiter in icy planetesimals. Our detailed discussion of calibration and error analysis supplements previously reported GPMS measurements of condensible volatile mixing ratios (Niemann et al., 1998, J. Geophys. Res. 103, 22831-22846; Atreya et al., 1999, Planet. Space Sci. 47, 1243-1262; Atreya et al., 2003, Planet. Space Sci. 51, 105-112) and the nitrogen isotopic ratio (Owen et al., 2001b, Astrophys. J. Lett. 553, L77-L79). The approximately three times solar abundance of NH₃ (along with CH₄ and H₂S) is consistent with enrichment of Jupiter's atmosphere by icy planetesimals formed at temperatures <40 K (Owen et al., 1999, Nature 402 (6759), 269-270), but would imply that H₂O should be at least 3×solar as well. An alternate model, using clathrate hydrates to deliver the nitrogen component to Jupiter, predicts O/H?9×solar (Gautier et al., 2001, Astrophys. J. 550 (2), L227-L230). Finally we show that the measured condensible volatile vertical profiles in the PES are consistent with column-stretching or entraining downdraft scenarios only if the basic state (the pre-stretched column or the entrainment source region) is described by condensible volatile vertical profiles that are drier than those in the equilibrium cloud condensation models. This dryness is supported by numerous remote sensing results but seems to disagree with observations of widespread clouds on Jupiter at pressure levels predicted by equilibrium cloud condensation models for ammonia and H₂S.     

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.     

Different eigenproblem models for field line resonances in cold plasma: Effect on magnetospheric density estimates

Magnetospheric plasma density can be remotely sensed through ground-based magnetometer data using a suitable model for field line resonances (FLRs) formed by standing shear Alfvén wave on closed geomagnetic field lines. The simplest type of FLR model, which is also the most relevant for magnetometer data inversion purposes, is based on solving a certain eigenvalue problem. Over the years a number of such models have been developed [Singer, H.J., Southwood, D.J., Walker, R.J., Kivelson, M.G., 1981. Alfvén wave resonances in a realistic magnetospheric magnetic field geometry. J. Geophys. Res. 86, 4589–4596; Rankin, R., Fenrich, F., Tikhonchuk, V.T., 2000. Shear Alfvén waves on stretched magnetic field lines near midnight in Earth's magnetosphere. Geophys. Res. Lett. 27, 3265–3268; Rankin, R., Kabin, K., Marchand, R., 2006. Alfvénic field line resonances in arbitrary magnetic field topology. Adv. Space Res. 38, 1720–1729]. In this paper we summarize the properties of these models and investigate the effect of using these different models on the magnetospheric density inferred from the ground-based magnetometer measurements. We also formulate a simple criterion which can be used to determine which one of these models should be used for a particular field line.     

Partial frequency redistribution of polarized radiation

Resonance polarization, which is created by the scattering of an anisotropic radiation field in regions of zero or weak magnetic fields, is strongly dependent on the frequency redistribution taking place during the scatterings. Here we discuss the frequency redistribution matrix relevant to resonance lines, concentrating on linear polarization. First we analyze in detail the redistribution matrix in a zero magnetic field given by the theory of Omont, Smith and Cooper (1972), revisited by Domke and Hubeny (1988). We explain that the linear polarization maxima which may appear in the wings of the Stokes Q profiles of strong resonances lines such as the Ca I 4227 Å line are coherent frequency redistribution effects. Various approximate forms of the frequency redistribution matrix are also examined. For resonance polarization in a weak magnetic field, we suggest a new expression for the redistribution matrix which can be used at all line frequencies, and is consistent with the condition that the Hanle effect acts only in the line core.     

Joule heating of Io's ionosphere by unipolar induction currents

The unexpectedly large scale height of Io's ionosphere (Kliore, A., et al., 1975, Icarus24, 407–410) together with the relatively large molecular weight of the likely principal constituent, SO₂ (Pearl, J., et al., 1979, Nature280, 755–758), suggest a high ionospheric temperature. Electrical induction in Io's ionosphere due to the corotating plasma bound to the Jovian magnetosphere is one possible source for attainment of such high temperatures. Accordingly, unipolar induction models were constructed to calculate ionospheric joule heating numerically. Heating rates produced by highly simplified models lie in the range 10^?9 to 10^?8 W/m³. These heating rates are lower than those determined from uv photodissociative heating models (Kumar, S., 1980, Geophys. Res. Lett.7, 9–12) at low levels in the ionosphere but are comparable in the upper ionosphere. The low electrical heating rate throughout most of the ionosphere is due to the power limitation imposed by the Alfvén wings which complete the electrical circuit (Neubauer, F.M., 1980, J. Geophys. Res.85, 1171–1178). Contrary to the pre-Voyager calculations of Cloutier, P. A., et al. (1978, Astrophys. Space Sci.55, 93–112), our numerical results show that the J × B force density due to unipolar induction currents in the ionosphere is much less than the gravitational force density when the combined mass of the neutral species is included. The binding and coupling of the ionosphere is principally due to the relatively dense (possibly localized) neutral SO₂ atmosphere. In regions where the ions and neutrals are collisionally coupled the ionosphere will not be stripped off by the J × B forces. However at a level above that (to which the ions move by diffusion only) the charged species would be removed. Thus there appears to be no need to postulate the existence of an intrinsic Ionian magnetic field as suggested by Kivelson, M. G., et al. (79, Science 205, 491–493) and Southwood, S. J., et al. (1980, J. Geophys. Res., in press) in order to retain the observed ionosphere.     

The variation of Io's auroral footprint brightness with the location of Io in the plasma torus

Ultraviolet and near-infrared observations of auroral emissions from the footprint of Io's magnetic Flux Tube (IFT) mapping to Jupiter's ionosphere have been interpreted via a combination of the unipolar inductor model [Goldreich, P., Lynden-Bell, D., 1969. Astrophys. J. 156, 59-78] and the multiply-reflected Alfvén wave model [Belcher, J.W., 1987. Science 238, 170-176]. While both models successfully explain the general nature of the auroral footprint and corotational wake, and both predict the presence of multiple footprints, the details of the interaction near Io are complicated [Saur, J., Neubauer, F.M., Connerney, J.E.P., Zarka, P., Kivelson, M.G., 2004. In: Bagenal, F., Dowling, T.E., McKinnon, W.B. (Eds.), Jupiter: The Planet, Satellites and Magnetosphere. Cambridge University Press, Cambridge, UK, pp. 537-560; Kivelson, M.G., Bagenal, F., Kurth, W.S., Neubauer, F.M., Paranicas, C., Saur, J., 2004. In: Bagenal, F., Dowling, T.E., McKinnon, W.B. (Eds.), Jupiter: The Planet, Satellites and Magnetosphere. Cambridge University Press, Cambridge, UK, pp. 513-536]. The auroral footprint brightness is believed to be a good remote indicator of the strength of the interaction near Io, indicating the energy and current strength linking Io with Jupiter's ionosphere. The brightness may also depend in part on local auroral acceleration processes near Jupiter. The relative importance of different physical processes in this interaction can be tested as Jupiter's rotation and Io's orbital motion shift Jupiter's magnetic centrifugal equator past Io, leading to longitudinal variations in the plasma density near Io and functionally different variations in the local field strength near Jupiter where the auroral emissions are produced. Initial HST WFPC2 observations found a high degree of variability in the footprint brightness with time, and some evidence for systematic variations with longitude [Clarke, J.T., Ben Jaffel, L., Gérard, J.-C., 1998. J. Geophys. Res. 103, 20217-20236], however the data were not of sufficient quality to determine functional relationships. In this paper we report the results from a second, more thorough study, using a series of higher resolution and sensitivity HST STIS observations and a model for the center to limb dependence of the optically thin auroral emission brightness based on measurements of the auroral curtain emission distribution with altitude. A search for correlations between numerous parameters has revealed a strong dependence between Io's position in the plasma torus and the resulting footprint brightness that persists over several years of observations. The local magnetic field strength near Jupiter (i.e. the size of the loss cone) and the expected north/south asymmetry in auroral brightness related to the path of currents generated near Io through the plasma torus en route to Jupiter appear to be less important than the total plasma density near Io. This is consistent with the near-Io interaction being dominated by collisions of corotating plasma and mass pickup, a long-standing view which has been subject to considerable debate. The brightness of the auroral footprint emissions, however, does not appear to be proportional to the incident plasma density or energy, and the interpretation of this result will require detailed modeling of the interaction near Io.