Models of Jupiter and Saturn after Galileo mission

New models of Jupiter are based on observational data provided by the Galileo spaceprobe, which considerably improved previously existing estimates of the helium abundance in the atmosphere of Jupiter. These data yield for Jupiter’s atmosphere 20% of the solar oxygen abundance and do not agree with the results of the analysis of the collision of comet Shoemaker-Levy 9 with Jupiter (10 times the solar value). Therefore, both the models of Jupiter with water-depleted and water-enriched atmosphere are considered. By analogy with Jupiter, trial models of Saturn with a water-depleted external envelope are also developed. The molecular-metallic phase transition pressure of hydrogen P_m was taken to be 1.5, 2 and 3 Mbar. Since Saturn’s internal molecular envelope is noticeably enriched in the IR-component (its weight concentration, 0.25–0.30, being by a factor of 3–4 higher than in Jupiter), the phase transition pressure in Saturn can be lower than in Jupiter. In the constructed models, the IR-core masses are 3–3.5 M_⊕ for Jupiter and 3–5.5 M_⊕ for Saturn. Jupiter’s and Saturn’s IR-cores can be considered embryos onto which the accretion of the gas occurred during the formation of the planets. The mass of the hydrogen–helium component dispersed in the zone of planetary formation constitutes ≈2–5 planetary masses for Jupiter and ≈11–14 planetary masses for Saturn.     

Galileo Photopolarimetry of Jupiter at 678.5 nm

Brightness and linear polarization measurements at 678.5 nm for four south-north strips of Jupiter are studied. These measurements were obtained in 1997 by the Galileo photopolarimeter/radiometer. The observed brightness exhibits latitudinal variations consistent with the belt/zone structure of Jupiter. The observed degree of linear polarization is small at low latitudes and increases steeply toward higher latitudes. No clear correlations were observed between the degree of linear polarization and the brightness. The observed direction of polarization changes from approximately parallel to the local scattering plane at low latitudes to perpendicular at higher latitudes. For our studies, we used atmospheric models that include a haze layer above a cloud layer. Parameterized scattering matrices were employed for the haze and cloud particles. On a pixel-wise basis, the haze optical thickness and the single-scattering albedo of the cloud particles were derived from the observed brightness and degree of linear polarization; results were accepted only if they were compatible with the observed direction of polarization. Using atmospheric parameter values obtained from Pioneer 10 and 11 photopolarimetry for the South Tropical Zone and the north component of the South Equatorial Belt, this analysis yielded acceptable results for very few pixels, particularly at small phase angles. However, for almost all pixels, acceptable results were found when the parameterized scattering matrix of the cloud particles was adjusted to produce more negative polarization for single scattering of unpolarized light, especially at large scattering angles, similar to some laboratory measurements of ammonia ice crystals. Using this adjusted model, it was found that the derived latitudinal variation of the single-scattering albedo of the cloud particles is consistent with the belt/zone structure, and that the haze optical thickness steeply increases toward higher latitudes.     

Implications from Galileo Observations on the Interior Structure and Chemistry of the Galilean Satellites

Data from the recent gravity measurements by the Galileo mission are used to construct wide ranges of interior structure and composition models for the Galilean satellites of Jupiter. These models show that mantle densities of Io and Europa are consistent with an olivine-dominated mineralogy with the ratios of Mg to Fe components depending on mantle temperature for Io and on ice shell thickness for Europa. The mantle density and composition depend relatively little on core composition. The size of the core is largely determined by the core's composition with core radius increasing with the concentration of a light component such as sulfur. For Io, the range of possible core sizes is between 38 and 53% of the satellite's radius. For Europa, there is also a substantial effect of the thickness of the ice layer which is varied between 120 and 170 km on the core size. Core sizes are between 10 and 45% of Europa's radius. The core size of Ganymede ranges between one-quarter and one-third of the surface radius depending on its sulfur content and the thickness of the ice shell. A subset of the Ganymede models is consistent with an olivine-dominated mantle mineralogy. The thickness of the silicate mantle above the core varies between 900 and 1100 km. The outermost ice shell is about 900 km in thickness and is further subdivided by pressure-induced phase transitions into ice I, ice III, ice V, and ice VI layers. Callisto should be differentiated, albeit incompletely. It is proposed that this satellite was never molten at a large scale but differentiated through the convective gradual unmixing of the ice and the metal/rock component. Bulk iron-to-silicon ratios Fe/Si calculated for the inner pair of satellites, Io and Europa, are less than the CI carbonaceous chondrite value of 1.7±0.1, whereas ratios for the outer pair, Ganymede and Callisto, cover a broad range above the chondritic value. Although the ratios are uncertain, in particular for Ganymede and Callisto, the values are sufficiently distinct to suggest a difference in composition between these two pairs of satellites. This may indicate a difference in iron-silicon fractionation during the formation of both classes of satellites in the protojovian nebula.     

Constraints on deep-seated zonal winds inside Jupiter and Saturn

The atmospheres of Jupiter and Saturn exhibit strong and stable zonal winds. How deep the winds penetrate unabated into each planet is unknown. Our investigation favors shallow winds. It consists of two parts. The first part makes use of an Ohmic constraint; Ohmic dissipation associated with the planet's magnetic field cannot exceed the planet's net luminosity. Application to Jupiter (J) and Saturn (S) shows that the observed zonal winds cannot penetrate below a depth at which the electrical conductivity is about six orders of magnitude smaller than its value at the molecular-metallic transition. Measured values of the electrical conductivity of molecular hydrogen yield radii of maximum penetration of 0.96RJ and 0.86RS, with uncertainties of a few percent of R. At these radii, the magnetic Reynolds number based on the zonal wind velocity and the scale height of the magnetic diffusivity is of order unity. These limits are insensitive to difficulties in modeling turbulent convection. They permit complete penetration along cylinders of the equatorial jets observed in the atmospheres of Jupiter and Saturn. The second part investigates how deep the observed zonal winds actually do penetrate. As it applies heuristic models of turbulent convection, its conclusions must be regarded as tentative. Truncation of the winds in the planet's convective envelope would involve breaking the Taylor-Proudman constraint on cylindrical flow. This would require a suitable nonpotential acceleration which none of the obvious candidates appears able to provide. Accelerations arising from entropy gradients, magnetic stresses, and Reynolds stresses appear to be much too weak. These considerations suggest that strong zonal winds are confined to shallow, stably stratified layers, with equatorial jets being the possible exception.     

Measurements of hot plasmas in the magnetosphere of Jupiter

This paper presents an overview of a number of the principal findings regarding the hot plasmas (E 50 keV) in Jupiter's magnetosphere by the HISCALE instrument during the encounter of the Ulysses spacecraft with the planet in February 1992. The hot plasma ion fluxes measured by HI-SCALE in the dayside magnetosphere are similar to those measured in the same energy range in this region by the Voyager spacecraft in 1979. Within the dayside plasma sheet, the hot-ion energy densities are comparable with, or larger than, the magnetic field energy densities; these hot ions are found to corotate at about one-half the planetary corotational speed. For ions of energies 500 keV/nucleon, the protons contributed from 50–60% to as much as 80% of the energy content of these plasmas. Strong, magnetic-field-aligned streaming was found for both the ions and electrons in the high-latitude duskside magnetosphere. The ion and electron pitch-angle distributions could be characterized by cos²⁵ α throughout many of the high anisotropy intervals of the outbound pass. There is some evidence in the ion pitch-angle distributions for a corotational component in the hot plasmas at high Jovian latitudes. While there are limitations owing to the finite geometries of the detector telescope systems on the determination of the angular spreads of the ion and electron beams, the measurements show that there are intervals when the particle distributions are not bidirectional. At such times, locally the hot plasmas could be carrying currents of 10⁻⁴μAm⁻². The temporal variations in the streaming electron fluxes are substantially larger than the variations measured for the fluxes that are more locally mirroring. The temporal variations contain periodicities that may correspond to hydromagnetic wave frequencies in the magnetosphere as well as to larger scale motions of magnetospheric plasmas. On nearly half of the days for about a 130 day interval around the time of the Ulysses encounter with the planet, particles of Jovian origin were measured in the interplanetary medium. An event discussed herein shows evidence of an energy dependence of the particle release process from the planetary magnetosphere into the interplanetary medium.     

Simulations of high-latitude spots on Jupiter: Constraints on vortex strength and the deep wind

We combine high-resolution observations of the dynamical behavior of small vortices (diameters ?5000 km) located at latitude 60°N on Jupiter with forward modeling, using the EPIC atmospheric model, to address two open questions: the dependence of the zonal winds with depth, and the strength of vortices that are too small to apply cloud tracking to their internal structure. The observed drift rates of the vortices can only be reproduced in the model when the zonal winds increase slightly with depth below the cloud tops, with a vertical shear that is less than was measured at 7°N at the southern rim of a 5-μm hotspot by the Galileo Probe Doppler Wind Experiment (DWE). This supports the idea that Jupiter's vertical shear may vary significantly with latitude. Our simulations suggest that the morphology of the mergers between vortices mainly depends on their maximum tangential velocities, the best results occurring when the tangential velocity is close to the velocity difference of the alternating jets constraining the zone in which the vortices are embedded. We use this correlation, together with the high-resolution data available for the White Ovals, to derive an empirical relationship between the maximum tangential velocity of a jovian vortex and its size, normalized by the strength and size of the encompassing shear zone. The Great Red Spot stands out as a significant anomaly to this relationship, but interestingly it is becoming less so with time.     

Potential Jupiter astmospheric constituents: candidates for the mass spectrometer in the Galileo atmospheric probe

An array of carbon-hydrogen-nitrogen (CHN) compounds are formed when mixtures of ammonia and acetylene are photolyzed with ultraviolet (UV) light. Some of these compounds, formed in our investigation of pathways for HCN synthesis on Jupiter, have not been encountered previously in observational, theoretical or laboratory photochemical studies. Knowledge of these structures may facilitate the interpretation of the mass spectral data obtained by the Galileo probe as it passed through the atmosphere of Jupiter.     

Monte Carlo Radiative Transfer Modeling of Lightning Observed in Galileo Images of Jupiter

We study lightning on Jupiter and the clouds illuminated by the lightning using images taken by the Galileo orbiter. The Galileo images have a resolution of ∼25 km/pixel and are able to resolve the shape of single lightning spots, which have half widths (radii) at half the maximum intensity in the range 45-80 km. We compare the shape and width of lightning flashes in the images with simulated flashes produced by our 3D Monte Carlo light-scattering model.The model calculates Monte Carlo scattering of photons in a 3D opacity distribution. During each scattering event, light is partially absorbed. The new direction of the photon after scattering is chosen according to a Henyey-Greenstein phase function. An image from each direction is produced by accumulating photons emerging from the cloud in a small range (bins) of emission angles. The light source is modeled either as a point or a vertical line.A plane-parallel cloud layer does not always fit the data. In some cases the cloud over the light source appears to resemble cumulus clouds on Earth. Lightning is estimated to occur at least as deep as the bottom of the expected water cloud. For the six flashes studied, we find that the clouds above the lightning are optically thick (τ>5). Jovian flashes are more regular and circular than the largest terrestrial flashes observed from space. On Jupiter there is nothing equivalent to the 30-40-km horizontal flashes that are seen on Earth.     

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.     

Composition and formation of Mercury: Constraints from future electrical conductivity measurements

Among the terrestrial planets, Mercury's composition is characterized by two specific features: a high density and a low surface FeO content. Based on these two constraints, different geochemical models have been proposed, according to different formation scenarios. Here thermodynamical modeling is used to derive the mantle and crust mineralogy associated with these geochemical models. For each mineralogy, the electrical conductivity profile and associated electromagnetic data are computed. Due to the very different oxide/silicate ratios, most geochemical models proposed for Mercury's formation show very different electromagnetic signatures. As a result, future measurements with MESSENGER and BepiColombo missions will help differentiating between different interior models and different formation scenarios.     

The deep water abundance on Jupiter: New constraints from thermochemical kinetics and diffusion modeling

We have developed a one-dimensional thermochemical kinetics and diffusion model for Jupiter’s atmosphere that accurately describes the transition from the thermochemical regime in the deep troposphere (where chemical equilibrium is established) to the quenched regime in the upper troposphere (where chemical equilibrium is disrupted). The model is used to calculate chemical abundances of tropospheric constituents and to identify important chemical pathways for CO-CH₄ interconversion in hydrogen-dominated atmospheres. In particular, the observed mole fraction and chemical behavior of CO is used to indirectly constrain the jovian water inventory. Our model can reproduce the observed tropospheric CO abundance provided that the water mole fraction lies in the range (0.25-6.0) × 10⁻³ in Jupiter’s deep troposphere, corresponding to an enrichment of 0.3-7.3 times the protosolar abundance (assumed to be H₂O/H₂ = 9.61 × 10⁻⁴). Our results suggest that Jupiter’s oxygen enrichment is roughly similar to that for carbon, nitrogen, and other heavy elements, and we conclude that formation scenarios that require very large (>8× solar) enrichments in water can be ruled out. We also evaluate and refine the simple time-constant arguments currently used to predict the quenched CO abundance on Jupiter, other giant planets, and brown dwarfs.     

木卫二的内部结构模型和海洋层--伽利略探测的最新研究进展

伽利略探测获得了木卫二的最新探测资料，其更精确的重力探测数据使深入研究木卫二的内部结构成为可能，研究表明，木卫二有着比较复杂的内部结构，最外层是水冰层，中层是岩石幔，中心是金属核（Fe或Fe-FeS)。但重力探测数据并不能确定水冰层中水的物态。多方面的证据（如理论，地学，磁场等证据）都表明木卫二的表面冰层下存在液体的导电海洋层，其中磁场方面的证据可能最为有力，概述了对木卫地的内部结构模型和海洋层两方面的研究，并推断潮汐演化是木卫二海洋层存在重要原因。     

Phosphine on Jupiter and Saturn from Cassini/CIRS

The global distribution of phosphine (PH₃) on Jupiter and Saturn is derived using 2.5 cm⁻¹ spectral resolution Cassini/CIRS observations. We extend the preliminary PH₃ analyses on the gas giants [Irwin, P.G.J., and 6 colleagues, 2004. Icarus 172, 37-49; Fletcher, L.N., and 9 colleagues, 2007a. Icarus 188, 72-88] by (a) incorporating a wider range of Cassini/CIRS datasets and by considering a broader spectral range; (b) direct incorporation of thermal infrared opacities due to tropospheric aerosols and (c) using a common retrieval algorithm and spectroscopic line database to allow direct comparison between these two gas giants.The results suggest striking similarities between the tropospheric dynamics in the 100-1000 mbar regions of the giant planets: both demonstrate enhanced PH₃ at the equator, depletion over neighbouring equatorial belts and mid-latitude belt/zone structures. Saturn's polar PH₃ shows depletion within the hot cyclonic polar vortices. Jovian aerosol distributions are consistent with previous independent studies, and on Saturn we demonstrate that CIRS spectra are most consistent with a haze in the 100-400 mbar range with a mean optical depth of 0.1 at 10 μm. Unlike Jupiter, Saturn's tropospheric haze shows a hemispherical asymmetry, being more opaque in the southern summer hemisphere than in the north. Thermal-IR haze opacity is not enhanced at Saturn's equator as it is on Jupiter.Small-scale perturbations to the mean PH₃ abundance are discussed both in terms of a model of meridional overturning and parameterisation as eddy mixing. The large-scale structure of the PH₃ distributions is likely to be related to changes in the photochemical lifetimes and the shielding due to aerosol opacities. On Saturn, the enhanced summer opacity results in shielding and extended photochemical lifetimes for PH₃, permitting elevated PH₃ levels over Saturn's summer hemisphere.     

Gas Dissipation from the Protosatellite Disk of Jupiter

We consider the dissipation of the gaseous component from the gas–dust accretion disk of Jupiter in which the Galilean satellites were formed. The thermal dissipation of hydrogen and helium is shown to be ineffective. It could ensure the loss of gas only for a low-mass disk and only if the rarefied outer layers of the disk are heated to 10⁴ K. Such a high disk temperature is not reached through Jupiter's radiation in existing models of its formation, but it could be provided by UV radiation of the early Sun after the dissipation of the protoplanetary disk. The viscous dissipation (with a viscosity parameter 10^–3 in the -disk model) related to disk accretion onto Jupiter could disperse a low-mass disk in 10⁷ years. A magnetocentrifugal mechanism, which produced a disk wind during accretion capable of carrying away 0.1 of the accreted gas mass, was probably also involved in the dispersal of the Jovian disk. Differential dispersion, with the loss of only hydrogen and helium and the retention of water vapor and heavier gases in the disk, is possible only in a low-mass disk model. We conclude that the water contained in the Galilean satellites was brought in mainly by solid planetesimals captured into the disk during mutual inelastic collisions in Jupiter's sphere of influence.     

Life cycles of spots on Jupiter from Cassini images

Using the sequence of 70-day continuum-band (751 nm) images from the Cassini Imaging Science System (ISS), we record over 500 compact oval spots and study their relation to the large-scale motions. The ∼100 spots whose vorticity could be measured—the large spots in most cases—were all anticyclonic. We exclude cyclonic features (chaotic regions) because they do not have a compact oval shape, but we do record their interactions with spots. We distinguish probable convective storms from other spots because they appear suddenly, grow rapidly, and are much brighter than their surroundings. The distribution of lifetimes for spots that appeared and disappeared during the 70-day period follows a decaying exponential with time constant (mean lifetime) of 3.5 days for probable convective storms and 16.8 days for all other spots. Extrapolating the exponential beyond 70 days seriously underestimates the number of spots that existed for the entire 70-day period. This and other evidences (size, shape, distribution in latitude) suggest that these long-lived spots with lifetime larger than 70 days are from a separate population. The zonal wind profile obtained manually by tracking individual features (this study) agrees with that obtained automatically by correlating brightness variations in narrow latitude bands (Porco et al., 2003). Some westward jets have developed more curvature and some have developed less curvature since Voyager times, but the number of westward jets that violate the barotropic stability criterion is about the same. In the northern hemisphere the number of spots is greatest at the latitudes of the westward jets, which are the most unstable regions according to the barotropic stability criterion. During the 70-day observation period the Great Red Spot (GRS) absorbed nine westward-moving spots that originated in the South Equatorial Belt (SEB), where most of the probable convective storms originate. Although the probable convective storms do not directly transform themselves into westward-moving spots, their common origin in the SEB suggests that moist convection and the westward jet compose a system that has maintained the GRS over its long lifetime.     

Composition of Comets: Observations and Models

We analyze the chemical composition and abundances of comets based on in situ measurements of Comet 1P/Halley and remote sensing observations of several recent bright comets including Hale-Bopp (C/1995 O1) and Hyakutake (C/1996 B2), in light of the elemental abundances of the solar system. Nitrogen is underabundant in comets relative to the solar system because nitrogen tends to be in N₂, which is chemically relatively inert. While many details remain uncertain, some gross features are emerging. The abundance of water : silicates: carbonaceous molecules (CO, CO₂, and hydrocarbons) by mass is approximately 1 : 1 : 1. Furthermore, the mass abundance of ice : dust (silicates and hydrocarbon polycondensates) is about1 : 1. We compare a list of identified comet molecules with molecules detected in the interstellar medium, although a comparison with their relative abundances, particularly in the ice phase, would be more meaningful. However, ice-phase abundances are not yet available. One can expect a variation of the abundances of carbon-bearing molecules in comets to be associated with their place of origin in the solar nebula. However, we also note that comets are heterogeneous. Thus, observed differences may be related to the place of origin, heterogeneity of the nucleus, or acquired through evolution. The molecular and elemental compositions of the coma are most likely not the same as those in the nucleus. This is particularly true for volatile ices and their gases and for the dust-to-ice and dust-to-gas ratios. Analyses must carefully consider the three sources of gas: Water from the surface of the nucleus, gases more volatile than water from the interior of the nucleus, and gases from the sublimation of the dust distributed in the coma. Topography on the surface of the nucleus may cause important evolutionary differences in the dust-to-gas mass ratio. Relatively inactive areas on the surface of the nucleus are probably associated with convex topography. Gas sublimated from convex areas (hills and mountains) diverges more strongly relative to gas sublimated from concave areas, which can entrain dust more efficiently. Thus, the entrainment of dust from convex areas is poor and dust may fall back to the surface of the nucleus creating a dust mantle, which further inhibits outgassing.     

Jupiter: New retrieved thermal profiles and ammonia distributions

New thermal profiles of Jupiter are retrieved from recent far infrared spectral measurements and for H₂ mixing ratios varying from 0.8 to 0.94. The effective temperature corresponding to the inferred thermal profile is 123.15 ± 0.35°K. Far-infrared brightness temperature spectra computed from these profiles are compared to experimental data including measurements made at high spectral resolution in the NH₃ν₂ band at 10 μm and in NH₃ pure rotational bands between 40 and 110 μm. It is found that a strong depletion of NH₃ does occur in the Jovian stratosphere and that ammonia seems to be undersaturated in the upper troposphere.