Gravitational signature of rotationally distorted Jupiter caused by deep zonal winds

Both deep zonal winds, if they exist, and the basic rotational distortion of Jupiter contribute to its zonal gravity coefficients J_n for n ? 2. In order to capture the gravitational signature of Jupiter that is caused solely by its deep zonal winds, one must take into account the full effect of rotational distortion by computing the coefficients J_n in non-spherical geometry. This represents a difficult and challenging problem because the widely-used spherical-harmonic-expansion method becomes no longer suitable. Based on the model of a polytropic Jupiter with index unity, we compute Jupiter’s gravity coefficients J₂, J₄, J₆, … , J₁₂ taking into account the full effect of rotational distortion of the gaseous planet using a finite element method. For the model of deep zonal winds on cylinders parallel to the rotation axis, we also compute the variation of the gravity coefficients ΔJ₂, ΔJ₄, ΔJ₆, … , ΔJ₁₂ caused solely by the effect of the winds in non-spherical geometry. It is found that the effect of the zonal winds on lower-order coefficients is weak, ∣ΔJ_n/J_n∣ < 1%, for n = 2, 4, 6, but it is substantial for the high-degree coefficients with n ? 8.     

Three-body capture of irregular satellites: Application to Jupiter

We investigate a new theory of the origin of the irregular satellites of the giant planets: capture of one member of a ∼100-km binary asteroid after tidal disruption. The energy loss from disruption is sufficient for capture, but it cannot deliver the bodies directly to the observed orbits of the irregular satellites. Instead, the long-lived capture orbits subsequently evolve inward due to interactions with a tenuous circumplanetary gas disk.We focus on the capture by Jupiter, which, due to its large mass, provides a stringent test of our model. We investigate the possible fates of disrupted bodies, the differences between prograde and retrograde captures, and the effects of Callisto on captured objects. We make an impulse approximation and discuss how it allows us to generalize capture results from equal-mass binaries to binaries with arbitrary mass ratios.We find that at Jupiter, binaries offer an increase of a factor of ∼10 in the capture rate of 100-km objects as compared to single bodies, for objects separated by tens of radii that approach the planet on relatively low-energy trajectories. These bodies are at risk of collision with Callisto, but may be preserved by gas drag if their pericenters are raised quickly enough. We conclude that our mechanism is as capable of producing large irregular satellites as previous suggestions, and it avoids several problems faced by alternative models.     

Subsurface Zonal Flows

We study the zonal flow in solar subsurface layers, analyzing about six years of GONG++ high-resolution Doppler data with ring-diagram analysis. We focus on the variation of zonal flow with magnetic activity over a range of depths from the surface to about 16 Mm. There is a positive correlation between unsigned magnetic flux and zonal flow at most depths. We calculate the average zonal flow for a quiet- and an active-region subset defined as dense-pack locations with an unsigned magnetic flux less than 3.4 G and locations with greater than 65.0 G, respectively. The average zonal flow of active regions is about 4 m s⁻¹ larger than the average flow of quiet regions. This difference increases slightly with increasing depth, which might be explained by a nonradial inclination of the flux tubes or a different extent in depth of different magnetic features. The difference shows no apparent pattern in time and latitude, which makes it unlikely that it is simply a manifestation of the torsional-oscillation pattern. As a byproduct, we find that the size of the North – South asymmetry of the rotation rate decreases during the same epoch.     

An analytical model of the Comet's collision with Jupiter

An important feature observed in the wake of the Jupiter-comet clash was the appearance of the ring structure axisymmetrically positioned around the center of the impact. The persistent expansion of the dark rings and its speed indicated an outward propagating gravity wave (Benka, 1995). We employ an analytical model of constant density, uniform finite depth and inviscid fluid layer to investigate the wave motion produced by the impact of Comet Shoemaker-Levy 9 on the Jovian atmosphere. The relevant thermal effects are neglected and an explosion resulting from the collision is then described by an initial impulsive pressure at the surface of the Jovian atmosphere. Under the assumption that all the kinetic energy of a comet fragment is completely converted into the energy of wave motions in the Jovian atmosphere, an analytical formula describing the relationship between the resulting wave motion in the atmosphere and the parameters of a comet fragment (the radius, density and speed) is derived. Results from the present simple analytical model give a qualitative agreement with observations regarding the distance and speed of the waves.     

The Cassini-Huygens flyby of Jupiter

The Cassini-Huygens spacecraft flew by Jupiter on December 30, 2000. The instruments aboard the spacecraft started making scientific observations three months earlier. Joint, collaborative observations were carried out with the teams of other spacecraft, notably Galileo, and with Earth-based observers. An operational overview of the flyby is presented and attention drawn to contributions of the eleven papers of this series which follow. Prime achievements of this campaign have been to better define the present state of fundamental elements of the jovian system, confirming many previously tentative conclusions. Particularly noteworthy is that the interactions between the solar wind and the jovian magnetosphere have been explored far deeper than before, along with the link to the morphology and dynamics of the jovian aurora.     

Seasonal Variations in the North-South Asymmetry of Polarized Light of Jupiter

Over an 18-year period, seasonal changes in the north-south asymmetry of polarization at high latitudes of Jupiter have been revealed by polarimetric observations in blue light. The average seasonal difference in the polarization degree between north and south is positive and equal to about 0.5%. There is some relationship between seasonal variations in the observed difference and the seasonal north-south asymmetry in solar radiation incident on Jupiter's atmosphere. There are two maxima on the observed seasonal curve, falling on the jovian spring and autumn and coinciding correspondingly with positive and negative maxima of the heliographic latitude of Jupiter. Two possible explanations are discussed: seasonal changes in insolation and/or time-dependent magnetospheric influence on the polar events.     

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.     

Carbon Monoxide on Jupiter: Evidence for Both Internal and External Sources

Thirteen lines of the CO band near 4.7 μm have been observed on a jovian hot spot at a resolution of 0.045 cm⁻¹. The measured line profiles indicate that the CO mole fraction is 1.0±0.2 ppb around the 6-bar level and is larger in the upper troposphere and/or stratosphere. An external source of CO providing an abundance of 4⁺³₋₂×10¹⁶ molecules cm⁻² is implied by the observations in addition to the amount deposited at high altitude by the Shoemaker-Levy 9 collision. From a simple diffusion model, we estimate that the CO production rate is (1.5-10)×10⁶ molecules cm⁻² s⁻¹ assuming an eddy diffusion coefficient around the tropopause between 300 and 1500 cm² s⁻¹. Precipitation of oxygen atoms from the jovian magnetosphere or photochemistry of water vapor from meteoroidal material can only provide a negligible contribution to this amount. A significant fraction of the CO in Jupiter's upper atmosphere may be formed by shock chemistry due to the infall of kilometer- to subkilometer-size Jupiter family comets. Using the impact rate from Levison et al. (2000, Icarus143, 415-420) rescaled by Bottke et al. (2002, Icarus156, 399-433), this source can provide the observed stratospheric CO only if the eddy diffusion coefficient around the tropopause is 100-300 cm² s⁻¹. Higher values, ∼700 cm² s⁻¹, would require an impact rate larger by a factor of 5-10, which cannot be excluded considering uncertainties in the distribution of Jupiter family comets. Such a large rate is indeed consistent with the observed cratering record of the Galilean satellites (Zahnle et al. 1998, Icarus136, 202-222). On the other hand, the ∼1 ppb concentration in the lower troposphere requires an internal source. Revisiting the disequilibrium chemistry of CO in Jupiter, we conclude that rapid vertical mixing can provide the required amount of CO at ∼6 bar for a global oxygen abundance of 0.2-9 times the solar value considering the uncertainties in the convective mixing rate and in the chemical constants.     

Cassini-VIMS at Jupiter: solar occultation measurements using Io

We report unusual and somewhat unexpected observations of the jovian satellite Io, showing strong methane absorption bands. These observations were made by the Cassini VIMS experiment during the Jupiter flyby of December/January 2000/2001. The explanation is straightforward: Entering or exiting from Jupiter's shadow during an eclipse, Io is illuminated by solar light which has transited the atmosphere of Jupiter. This light, therefore becomes imprinted with the spectral signature of Jupiter's upper atmosphere, which includes strong atmospheric methane absorption bands. Intercepting solar light refracted by the jovian atmosphere, Io essentially becomes a “mirror” for solar occultation events of Jupiter. The thickness of the layer where refracted solar light is observed is so large (more than 3000 km at Io's orbit), that we can foresee a nearly continuous multi-year period of similar events at Saturn, utilizing the large and bright ring system. During Cassini's 4-year nominal mission, this probing technique should reveal information of Saturn's atmosphere over a large range of southern latitudes and times.     

Quasi-satellites of Jupiter

Two asteroids 2001 QQ199 and 2004 AE9 and two comets P/LINEAR-Catalina and P/LINEAR are found to be quasi-satellites of Jupiter at present time. The Tisserand parameters of these four bodies lie between 2.3 and 2.5, which means two asteroids might be Jupiter family comets and will show a cometary’s activity in the near future.     

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.     

An incompressible stratified fluid model ofComet Shoemaker-Levy 9''s collision with Jupiter

A linear model of Comet Shoemaker-Levy 9's collision with Jupiter ispresented, in which the Jovian atmosphere is taken to be a stratified, incompressible and inviscid fluid layer and the gravity wave associated with the impacts is caused by an initial impulse. Instead of adopting an initial surface deformation, we take an initial impulsive pressure P(r, z, 0) as the initial condition. It is found that this approach is more convenient for the mathematical formulation of the problem, especially for the computation of the vertical displacement of the fluid particles and the perturbation potential energy which can then be expressed by finite Fourier sine and cosine transform of P(r, z, 0).

A relationship among parameters of the impact and the Jovian atmosphereand the assumed initial impulse is deduced. We also derive the wave height of the surface wave analytically for two assumed forms of P(r, z, 0).     

Structure and evolution of the Jupiter family

Both physical and dynamical issues are important in order to judge the origin and evolution of the Jupiter family of short-period comets. The steady-state condition for maintaining this structure at its present size by captures from the classical Oort cloud is reviewed on the basis of recent results concerning the absolute number of Jupiter family comets as a function of perihelion distance as well as the coupled physical and dynamical evolutions that evidently occur. Like in most earlier investigations, a clear shortage is found in the classical Oort cloud source. The shortage seems, however, less extreme than sometimes assumed. Monte Carlo simulations are envisaged as a way to shed further light on the fate of Jupiter family comets.     

A high-resolution study of Jupiter at a frequency of 30 GHz

We present the results of our study of Jupiter and its radiation belts with a resolution of 6 arcsec at a frequency of 30 GHz using the RATAN-600 radio telescope and a MARS matrix radiometer with a sensitivity of about 6 mK ^?1/2. We monitored the integrated emission from the Jovian disk with a signal-to-noise ratio of more than 1000 for 30 days and showed its radio emission to be highly stable (≈1%). Based on daily data for the one-dimensional radio brightness distribution over the disk, we mapped the longitudinal radio brightness distribution over 100 rotation periods of Jupiter around its axis. Neither hot nor cold spots with a temperature contrast of more than 1 K were detected; their contribution to the total radio flux from the Jovian disk was no more than 0.2%. The one-dimensional latitudinal (longitude-averaged) distribution obtained on VLA with a similar resolution is shown to be an order of magnitude less uniform than the one-dimensional longitudinal (latitude-averaged) distribution obtained on RATAN-600. We have studied the radiation belts at such high frequencies for the first time and estimated their intensities and variability levels under the effect of external factors. The variable component of the radiation belts was shown to have not exceeded 0.5% of the integrated spectrum of Jupiter over the entire period of its observations. We estimated the contribution of the Galilean satellites (“Galilean noise”) in low-resolution observations; the accuracy of allowing for this noise is determined by the accuracy of estimating the temperatures of the satellites at the observing frequency. The uncertainty in the total flux does not exceed 0.1%.     

The deep wind structure of the giant planets: Results from an anelastic general circulation model

The giant gas planets have hot convective interiors, and therefore a common assumption is that these deep atmospheres are close to a barotropic state. Here we show using a new anelastic general circulation model that baroclinic vorticity contributions are not negligible, and drive the system away from an isentropic and therefore barotropic state. The motion is still aligned with the direction of the axis of rotation as in a barotropic rotating fluid, but the wind structure has a vertical shear with stronger winds in the atmosphere than in the interior. This shear is associated with baroclinic compressibility effects. Most previous convection models of giant planets have used the Boussinesq approximation, which assumes the density is constant in depth; however, Jupiter's actual density varies by four orders of magnitude through its deep molecular envelope. We therefore developed a new general circulation model (based on the MITgcm) that is anelastic and thereby incorporates this density variation. The model's geometry is a full 3D sphere down to a small inner core. It is nonhydrostatic, uses an equation of state suitable for hydrogen-helium mixtures (SCVH), and is driven by an internal heating profile. We demonstrate the effect of compressibility by comparing anelastic and Boussinesq cases. The simulations develop a mean state that is geostrophic and hydrostatic including the often neglected, but significant, vertical Coriolis contribution. This leads to modification of the standard thermal wind relation for a deep compressible atmosphere. The interior flow organizes in large cyclonically rotating columnar eddies parallel to the rotation axis, which drive upgradient angular momentum eddy fluxes, generating the observed equatorial superrotation. Heat fluxes align with the axis of rotation, and provide a mechanism for the transport of heat poleward, which can cause the observed flat meridional emission. We address the issue of over-forcing which is common in such convection models and analyze the dependence of our results on this; showing that the vertical wind structure is not very sensitive to the Rayleigh number. We also study the effect of rotation, showing how the transition from a rapidly to a slowly rotating system affects the dynamics.     

Phase separation in giant planets: inhomogeneous evolution of Saturn

We present the first models of Jupiter and Saturn to couple their evolution to both a radiative-atmosphere grid and to high-pressure phase diagrams of hydrogen with helium and other admixtures. We find that prior calculated phase diagrams in which Saturn's interior reaches a region of predicted helium immiscibility do not allow enough energy release to prolong Saturn's cooling to its known age and effective temperature. We explore modifications to published phase diagrams that would lead to greater energy release, and propose a modified H-He phase diagram that is physically reasonable, leads to the correct extension of Saturn's cooling, and predicts an atmospheric helium mass fraction Y_atmos=0.185, in agreement with recent estimates. We also explore the possibility of internal separation of elements heavier than helium, and find that, alternatively, such separation could prolong Saturn's cooling to its known age and effective temperature under a realistic phase diagram and heavy element abundance (in which case Saturn's Y_atmos would be solar but heavier elements would be depleted). In none of these scenarios does Jupiter's interior evolve to any region of helium or heavy-element immiscibility: Jupiter evolves homogeneously to the present day. We discuss the implications of our calculations for Saturn's primordial core mass.     

Estimates for the effective electrical conductivity of the core in the interior of Jupiter and Saturn

In the deep interior of the giant planets Jupiter and Saturn, ordinary hydrogen and helium are transformed into a conducting metallic liquid at extremely high pressure. It is likely that the giant planets' observed magnetic field is constantly generated in the metallic fluid core by magnetohydrodynamic processes, converting mechanic energy in the form of convection into magnetic energy. The maximum strength of their magnetic fields is likely to be limited by magnetic field instabilities which convert the magnetic energy back into convection. The parameter which governs the occurrence of magnetic instabilities is the Elsasser number, = B ²/2, where B is the field strength, is the electrical conductivity, is the rotation rate and is the density. Since magnetic instability will be very active when exceeds a critical value _c 10 (the precise value depending on the magnetic field distribution), this imposes an upper bound on the effective electrical conductivity of the metallic fluid which comprises the bulk of Jupiter's interior and much of Saturn's.Stability calculations including both toroidal (model) and poloidal (observed) components of the magnetic field in a rapidly rotating spherical shell, have been performed. The most stable configuration of the field is when the poloidal component of field is strong and the toroidal field is weak; in this case we obtain an upper bound for electrical conductivity of 3 × 10⁶ S/m; while the most unstable configuration of the field is when the toroidal and poloidal fields are comparable, giving rise to _m 3 × 10⁵ S/m. The implications of the results for general dynamo theory are also discussed.     

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

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

Long slit spectropolarimetry of Jupiter and Saturn

We present ground-based limb polarization measurements of Jupiter and Saturn consisting of full disk imaging polarimetry for the wavelength 7300 Å and spatially resolved (long-slit) spectropolarimetry covering the wavelength range 5200-9350 Å.For the polar region of Jupiter we find for λ = 6000 Å a very strong radial (perpendicular to the limb) fractional polarization with a seeing corrected maximum of about +11.5% in the South and +10.0% in the North. This indicates that the polarizing haze layer is thicker at the South pole. The polar haze layers extend down to 58° in latitude. The derived polarization values are much higher than reported in previous studies because of the better spatial resolution of our data and an appropriate consideration of the atmospheric seeing. Model calculations demonstrate that the high limb polarization can be explained by strongly polarizing (p ≈ 1.0), high albedo (ω ≈ 0.98) haze particles with a scattering asymmetry parameter of g ≈ 0.6 as expected for aggregate particles of the type described by West and Smith (West, R.A., Smith, P.H. [1991]. Icarus 90, 330-333). The deduced particle parameters are distinctively different when compared to lower latitude regions.The spectropolarimetry of Jupiter shows a decrease in the polar limb polarization towards longer wavelengths and a significantly enhanced polarization in strong methane bands when compared to the adjacent continuum. This is a natural outcome for a highly polarizing haze layer above an atmosphere where multiple scatterings are suppressed in absorption bands. For lower latitudes the fractional polarization is small, negative, and it depends only little on wavelength except for the strong CH₄-band at 8870 Å.The South pole of Saturn shows a lower polarization (p ≈ 1.0-1.5%) than the poles of Jupiter. The spectropolarimetric signal for Saturn decrease rapidly with wavelength and shows no significant enhancements in the fractional polarization in the absorption bands. These properties can be explained by a vertically extended stratospheric haze region composed of small particles <100 nm as suggested previously by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M. [2005]. Icarus 179, 195-221).In addition we find in the V- and R-band a previously not observed strong polarization feature (p = 1.5-2.0%) near the equator of Saturn. The origin of this polarization signal is unclear but it could be related to a seasonal effect.Finally we discuss the potential of ground-based limb polarization measurements for the investigation of the scattering particles in the atmospheres of Jupiter and Saturn.     

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.