Interaction between eddies and mean flow in Jupiter's atmosphere: Analysis of Cassini imaging data

Beebe et al. [Beebe, R.F., et al., 1980. Geophys. Res. Lett. 17, 1-4] and Ingersoll et al. [Ingersoll, A.P., et al., 1981. J. Geophys. Res. 86, 8733-8743] used images from Voyagers 1 and 2 to analyze the interaction between zonal winds and eddies in Jupiter's atmosphere. They reported a high positive correlation between Jupiter's eddy momentum flux, , and the variation of zonal velocity with latitude, . This correlation implied a surprisingly high rate of conversion of energy from eddies to zonal flow: , a value more than 10% of Jupiter's thermal flux emission. However, Sromovsky et al. [Sromovsky, L.A., et al., 1982. J. Atmos. Sci. 39, 1413-1432] argued that possible biases in the analysis could have caused an artificially high correlation. In addition, significant differences in the derived eddy flux between datasets put into question the robustness of any one result. We return to this long-standing puzzle using images of Jupiter from the Cassini flyby of December 2000. Our method is similar to previous analyses, but utilizes an automatic feature tracker instead of the human eye. The number of velocity vectors used in this analysis is over 200,000, compared to the 14,000 vectors used by Ingersoll et al. We also find a positive correlation between and and derive a global average power per unit mass, , ranging from . Utilizing Ingersoll et al.'s estimate of the mass per unit area involved in the transport, this would imply a rate of energy conversion of . We discuss the implications of this result and employ several tests to demonstrate its robustness.     

Orbits and masses of Saturn's coorbital and F-ring shepherding satellites

We have obtained numerically integrated orbits for Saturn's coorbital satellites, Janus and Epimetheus, together with Saturn's F-ring shepherding satellites, Prometheus and Pandora. The orbits are fit to astrometric observations acquired with the Hubble Space Telescope and from Earth-based observatories and to imaging data acquired from the Voyager spacecraft. The observations cover the 38 year period from the 1966 Saturn ring plane crossing to the spring of 2004. In the process of determining the orbits we have found masses for all four satellites. The densities derived from the masses for Janus, Epimetheus, Prometheus, and Pandora in units of g cm⁻³ are , , , and , respectively.     

Gravitational tidal waves in Titan's upper atmosphere

Tidal waves driven by Titan's orbital eccentricity through the time-dependent component of Saturn's gravitational potential attain nonlinear, saturation amplitudes (|T^′|>10 K, , and ) in the upper atmosphere (?500 km) due to the approximate exponential growth as the inverse square root of pressure. The gravitational tides, with vertical wavelengths of ∼100-150 km above 500 km altitude, carry energy fluxes sufficient in magnitude to affect the energy balance of the upper atmosphere with heating rates in the altitude range of 500-900 km.     

CCD photometry of 23 minor planets

We present CCD photometric observations of 23 main-belt asteroids, of which 8 have never been observed before; thus, the data of these objects are the first in the literature. The majority showed well-detectable light variations, exceeding 0^m1. We have determined synodic periods for 756 Lilliana (936), 1270 Datura (34), 1400 Tirela (1336), 1503 Kuopio (998), 3682 Welther (359), 7505 Furushu (414) and 11436 1969 QR (123), while uncertain period estimates were possible for 469 Argentina (123), 546 Herodias (104) and 1026 Ingrid (53). The shape of the lightcurves of 3682 Welther changed on a short time-scale and showed dimmings that might be attributed to eclipses in a binary system. For the remaining objects, only lower limits of the periods and amplitudes were concluded.     

CASSINI/VIMS-V at Jupiter: Radiometric calibration test and data results

During the Cassini-Huygens flyby of Jupiter in December 2000, VIMS-V acquired multispectral data cubes of Jupiter's atmosphere. The visual and infrared imaging spectrometer-visual channel (VIMS-V) is one of the principal contributions of Italian Space Agency (ASI) to the Cassini-Huygens mission to Saturn. VIMS-V is an imaging spectrometer operating in the wavelength range 300-, with a (nominal) spectral resolution of , and a (nominal) spatial resolution of . VIMS-V is boresighted with the VIMS-IR channel operating in the wavelength range 0.8-. During the early phases of the Cassini mission, the spacecraft encountered Venus (June 23, 1999), followed shortly thereafter by a flyby of the Earth. During the Earth flyby the Moon (August 17, 1999) was observed. Following the Earth-Moon flyby, the spacecraft encountered Jupiter (closest approach on December 31, 2000), and during the roughly 6 months prior to Jupiter closest approach a series of observations were made of most of the objects in the Jovian system. We have determined the instrumental transfer function of VIMS-V using the Moon and Venus day side data. This transfer function was then used to remove instrumental effects from the Jupiter data and to convert raw instrumental response numbers to spectral radiance from the target. It was thus possible to study the spectral variability of Jupiter's atmosphere across its disk using data from both the visual (V) and infrared (IR) channels of VIMS. In this paper we discuss the main results obtained by the V channel. We have analyzed the principal spectral features of Jupiter atmosphere, and in particular, the spatial variation of methane and ammonia absorption bands over the Jovian disk. Using the instrument's spatial mapping capabilities we have investigated the nature of the absorption band in the spectrum of Jupiter's atmosphere at that is consistent with the presence of ammonia or water vapor. After comet Shoemaker-Levy 9 impacted Jupiter, water vapor was considered the most likely cause of the absorption feature, but our data indicate that ammonia is the source of this band. Other analyses were performed using standard techniques such as forming band ratios and removal of the continuum. Our analyses confirm previous ground or satellite based observations. We were also able to verify the instrument radiometric calibration, using observations conducted during the close encounters with Venus and the Moon.     

Modelling meteor ablation in the venusian atmosphere

A comparative study of meteor ablation in the atmospheres of the Earth and Venus is presented. The classical single body meteor ablation model is extended to incorporate a heat penetration depth estimate allowing the simulation of larger meteoroids, than would an isothermal model. The ablation of icy and rocky meteoroids, with densities of 1.0 and 3.4 g cm⁻³, respectively, and initial radii of up to for rock and for ice (equivalent to an initial mass of in both cases), was simulated in both atmospheres. In general venusian meteors are brighter than terrestrial equivalents. Large, slow, rocky objects may be up to 0.7 mag brighter on Venus, while small, icy particles with entry speeds in the range 30-60 km s⁻¹, are found to be upwards of 2.7 mag brighter than at the Earth. Venusian meteors reach maximum brightness at greater altitudes than would similar particles at the Earth. Rocky meteoroids have their points of maximum brightness some 15-35 km higher up at Venus, between 90 and 120 km, whereas, for icy particles this altitude difference is about 5-25 km higher up than at the Earth, in the range 100-125 km. These findings agree, for the most part, with recent analytical studies. Venusian meteors, which last from 100 ms to , tend to be shorter-lived than terrestrial meteors, with correspondingly shorter visible trails. Large (), slow () icy particles reach a maximum magnitude of ∼−2 at Venus and remain visible for about one second, with a large section of the smaller faster meteoroids simulated here remaining visible for several hundred milliseconds. In light of recent space-based meteor observations at the Earth [Jenniskens, P., Tedesco, E., Muthry, J., Laux, C.O., Price, S., 2002. Meteorit. Planet. Sci. 37, 1071-1078], such brightness, height and duration estimates as suggested in this work, may be used in developing future observational campaigns to be carried out from Venus orbit.     

Calculation of the enrichment of the giant planet envelopes during the “late heavy bombardment”

The giant planets of our solar system possess envelopes consisting mainly of hydrogen and helium but are also significantly enriched in heavier elements relatively to our Sun. In order to better constrain how these heavy elements have been delivered, we quantify the amount accreted during the so-called “late heavy bombardment”, at a time when planets were fully formed and planetesimals could not sink deep into the planets. On the basis of the “Nice model”, we obtain accreted masses (in terrestrial units) equal to for Jupiter, and for Saturn. For the two other giant planets, the results are found to depend mostly on whether they switched position during the instability phase. For Uranus, the accreted mass is with an inversion and without an inversion. Neptune accretes in models in which it is initially closer to the Sun than Uranus, and otherwise. With well-mixed envelopes, this corresponds to an increase in the enrichment over the solar value of 0.033±0.001 and 0.074±0.007 for Jupiter and Saturn, respectively. For the two other planets, we find the enrichments to be 2.1±1.4 (w/ inversion) or 1.2±0.7 (w/o inversion) for Uranus, and 2.0±1.2 (w/ inversion) or 2.7±1.6 (w/o inversion) for Neptune. This is clearly insufficient to explain the inferred enrichments of ∼4 for Jupiter, ∼7 for Saturn and ∼45 for Uranus and Neptune.     

A ballistics analysis of the Deep Impact ejecta plume: Determining Comet Tempel 1's gravity, mass, and density

In July of 2005, the Deep Impact mission collided a 366 kg impactor with the nucleus of Comet 9P/Tempel 1, at a closing speed of 10.2 km s⁻¹. In this work, we develop a first-order, three-dimensional, forward model of the ejecta plume behavior resulting from this cratering event, and then adjust the model parameters to match the flyby-spacecraft observations of the actual ejecta plume, image by image. This modeling exercise indicates Deep Impact to have been a reasonably “well-behaved” oblique impact, in which the impactor-spacecraft apparently struck a small, westward-facing slope of roughly 1/3-1/2 the size of the final crater produced (determined from initial ejecta plume geometry), and possessing an effective strength of not more than . The resulting ejecta plume followed well-established scaling relationships for cratering in a medium-to-high porosity target, consistent with a transient crater of not more than 85-140 m diameter, formed in not more than 250-550 s, for the case of (gravity-dominated cratering); and not less than 22-26 m diameter, formed in not less than 1-3 s, for the case of (strength-dominated cratering). At , an upper limit to the total ejected mass of 1.8×¹⁰⁷ kg (1.5-2.2×¹⁰⁷ kg) is consistent with measurements made via long-range remote sensing, after taking into account that 90% of this mass would have stayed close to the surface and then landed within 45 min of the impact. However, at , a lower limit to the total ejected mass of 2.3×¹⁰⁵ kg (1.5-2.9×¹⁰⁵ kg) is also consistent with these measurements. The expansion rate of the ejecta plume imaged during the look-back phase of observations leads to an estimate of the comet's mean surface gravity of (0.17-0.90 mm s⁻²), which corresponds to a comet mass of mt=4.5×¹⁰¹³ kg (2.3-12.0×¹⁰¹³ kg) and a bulk density of (200-1000 kg m⁻³), where the large high-end error is due to uncertainties in the magnitude of coma gas pressure effects on the ejecta particles in flight.     

Tvashtar awakening detected in April 2006 with OSIRIS at the W.M. Keck Observatory

We detected a volcanic outburst in Io's northern hemisphere on 17 April 2006 with the OSIRIS imaging spectrometer at Keck, and confirmed it was still erupting on 2 June 2006. The eruption, which we name 060417A, was located in Tvashtar Paterae, ∼100 km southeast of the February 2000 eruption. The observed temperature was , over a surface area of , providing a total thermal output of .     

Accretion of the gaseous envelope of Jupiter around a 5-10 Earth-mass core

New numerical simulations of the formation and evolution of Jupiter are presented. The formation model assumes that first a solid core of several M_⊕ accretes from the planetesimals in the protoplanetary disk, and then the core captures a massive gaseous envelope from the protoplanetary disk. Earlier studies of the core accretion-gas capture model [Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y., 1996. Icarus 124, 62-85] demonstrated that it was possible for Jupiter to accrete with a solid core of 10-30 M_⊕ in a total formation time comparable to the observed lifetime of protoplanetary disks. Recent interior models of Jupiter and Saturn that agree with all observational constraints suggest that Jupiter's core mass is 0-11 M_⊕ and Saturn's is 9-22 M_⊕ [Saumon, G., Guillot, T., 2004. Astrophys. J. 609, 1170-1180]. We have computed simulations of the growth of Jupiter using various values for the opacity produced by grains in the protoplanet's atmosphere and for the initial planetesimal surface density, σinit,Z, in the protoplanetary disk. We also explore the implications of halting the solid accretion at selected core mass values during the protoplanet's growth. Halting planetesimal accretion at low core mass simulates the presence of a competing embryo, and decreasing the atmospheric opacity due to grains emulates the settling and coagulation of grains within the protoplanet's atmosphere. We examine the effects of adjusting these parameters to determine whether or not gas runaway can occur for small mass cores on a reasonable timescale. We compute four series of simulations with the latest version of our code, which contains updated equation of state and opacity tables as well as other improvements. Each series consists of a run without a cutoff in planetesimal accretion, plus up to three runs with a cutoff at a particular core mass. The first series of runs is computed with an atmospheric opacity due to grains (hereafter referred to as ‘grain opacity’) that is 2% of the interstellar value and . Cutoff runs are computed for core masses of 10, 5, and 3 M_⊕. The second series of Jupiter models is computed with the grain opacity at the full interstellar value and . Cutoff runs are computed for core masses of 10 and 5 M_⊕. The third series of runs is computed with the grain opacity at 2% of the interstellar value and . One cutoff run is computed with a core mass of 5 M_⊕. The final series consists of one run, without a cutoff, which is computed with a temperature dependent grain opacity (i.e., 2% of the interstellar value for ramping up to the full interstellar value for ) and . Our results demonstrate that reducing grain opacities results in formation times less than half of those for models computed with full interstellar grain opacity values. The reduction of opacity due to grains in the upper portion of the envelope with has the largest effect on the lowering of the formation time. If the accretion of planetesimals is not cut off prior to the accretion of gas, then decreasing the surface density of planetesimals lowers the final core mass of the protoplanet, but increases the formation timescale considerably. Finally, a core mass cutoff results in a reduction of the time needed for a protoplanet to evolve to the stage of runaway gas accretion, provided the cutoff mass is sufficiently large. The overall results indicate that, with reasonable parameters, it is possible that Jupiter formed at 5 AU via the core accretion process in 1 Myr with a core of 10 M_⊕ or in 5 Myr with a core of 5 M_⊕.     

Detection of CH3D on Titan

We report the detection of ¹³CH₃D in Titan's stratosphere from Cassini/CIRS infrared spectra near 8.7 μm. Fitting simultaneously the ν₆ bands of both ¹³CH₃D and ¹²CH₃D and the ν₄ band of CH₄, we derive a D/H ratio equal to and a ¹²C/¹³C ratio in deuterated methane of , consistent with that measured in normal methane.