MESSENGER observations of the plasma environment near Mercury

The MESSENGER Fast Imaging Plasma Spectrometer (FIPS) measured the bulk plasma characteristics of Mercury's magnetosphere and solar wind environment during the spacecraft's first two flybys of the planet on 14 January 2008 (M1) and 6 October 2008 (M2), producing the first measurements of thermal ions in Mercury's magnetosphere. In this work, we identify major features of the Mercury magnetosphere in the FIPS proton data and describe the data analysis process used for recovery of proton density (n_p) and temperature (T_p) with a forward modeling technique, required because of limitations in measurement geometry. We focus on three regions where the magnetospheric flow speed is likely to be low and meets our criteria for the recovery process: the M1 plasma sheet and the M1 and M2 dayside and nightside boundary-layer regions. Interplanetary magnetic field (IMF) conditions were substantially different between the two flybys, with intense reconnection signatures observed by the Magnetometer during M2 versus a relatively quiet magnetosphere during M1. The recovered ion density and temperature values for the M1 quiet-time plasma sheet yielded n_p∼1–10 cm⁻³, T_p∼2×10⁶ K, and plasma β∼2. The nightside boundary-layer proton densities during M1 and M2 were similar, at n_p∼4–5 cm⁻³, but the temperature during M1 (T_p∼4–8×10⁶ K) was 50% less than during M2 (T_p∼8×10⁶ K), presumably due to reconnection in the tail. The dayside boundary layer observed during M1 had a density of ∼16 cm⁻³ and temperature of 2×10⁶ K, whereas during M2 this region was less dense and hotter (n_p∼8 cm⁻³ and T_p∼10×10⁶ K), again, most likely due to magnetopause reconnection. Overall, the southward interplanetary magnetic field during M2 clearly produced higher T_p in the dayside and nightside magnetosphere, as well as higher plasma β in the nightside boundary, ∼20 during M2 compared with ∼2 during M1. The proton plasma pressure accounts for only a fraction (24% for M1 and 64% for M2) of the drop in magnetic pressure upon entry into the dayside boundary layer. This result suggests that heavy ions of planetary origin, not considered in this analysis, may provide the “missing” pressure. If these planetary ions were hot due to “pickup” in the magnetosheath, the required density for pressure balance would be an ion density of ∼1 cm⁻³ for an ion temperature of ∼10⁸ K.     

Vertical structure of Jupiter's troposphere from nonlinear simulations of long-lived vortices

Numerical simulations of jovian vortices at tropical and temperate latitudes, under different atmospheric conditions, have been performed using the EPIC code [Dowling, T.E., Fisher, A.S., Gierasch, P.J., Harrington, J., LeBeau, R.P., Santori, C.M., 1998. Icarus 132, 221-238] to simulate the high-resolution observations of motions and of the lifetimes presented in a previous work [Legarreta, J., Sánchez-Lavega, A., 2005. Icarus 174, 178-191] and infer the vertical structure of Jupiter's troposphere. We first find that in order to reproduce the longevity and drift rate of the vortices, the Brunt-Väisälä frequency of the atmosphere in the upper troposphere (pressures P∼1 to 7 bar) should have a lower limit value of 5×¹⁰⁻³ s⁻¹, increasing upward up to 1.25×¹⁰⁻² s⁻¹ at pressures P∼0.5 bar (latitudes between 15° and 45° in both hemispheres). Second, the vortices drift also depend on the vertical structure of the zonal wind speed in the same range of altitudes. Simulations of the slowly drifting Southern hemisphere vortices (GRS, White Ovals and anticyclones at 40° S) require a vertically-constant zonal-wind with depth, but Northern hemisphere vortices (cyclonic “barges” and anticyclones at 19, 41 and 45° N) require decreasing winds at a rate of ∼5 m s⁻¹ per scale height. However vortices drifting at a high speed, close to or in the peak of East or West jets and in both hemispheres, require the wind speed slightly increasing with depth, as is the case for the anticyclones at 20° S and at 34° N. We deduce that the maximum absolute vertical shear of the zonal wind from P∼1 bar up to P∼7 bar in these jets is ∼15 m s⁻¹ per scale height. Intense vortices with tangential velocity at their periphery ∼100 m s⁻¹ tend to decay asymptotically to velocities ∼40 to 60 m s⁻¹ with a characteristic time that depends on the vortex intensity and static stability of the atmosphere. The vortices adjust their tangential velocity to the averaged peak to peak velocity of the opposed eastward and westward jets at their boundary. We show through our simulations that large-scale and long-lived vortices whose maximum tangential velocity is ∼100 m s⁻¹ can survive by absorbing smaller intense vortices.     

Jupiter's Cloud Structure as Constrained by Galileo Probe and HST Observations

The Galileo Probe sampled Jupiter's atmosphere at the edge of a 5-μm hot spot, where it found very little cloud opacity above the 700 mb level. Only τ=1-2 at λ=0.5 μm was inferred from Net Flux Radiometer observations (Sromovsky et al. 1998, J. Geophys. Res.103, 22,929-22,977), in seeming conflict with Chanover et al. (1997, Icarus128, 294-305) who inferred τ=6-8 above the 700 mb level (at λ∼0.9 μm) from 893-nm and 953-nm WFPC2 observations of a group of hot spots. Postulating a heterogeneous cloud structure is one way to resolve the conflict. We obtained a more satisfying resolution by reinterpretation of the HST observations with Probe-compatible assumptions about the vertical distribution of cloud particles. Assuming a physically thin upper (putative NH₃) cloud with adjustable optical depth and effective pressure (p_eff<440 mb) and a physically thin midlevel (putative NH₄SH) cloud with adjustable optical depth but a fixed pressure of 1.2 bars, we are able to fit WPFC2 observations with probe-consistent opacities in hot spot regions. With the same cloud pressures, but higher middle cloud opacities, we are even able to fit the visibly bright regions. Little variability is seen in the upper cloud. Best fits to October 1995 WFPC2 observations in dark regions (5-μm hot spots) yielded τ_upper=1.3-1.9 at 0.9 μm and p_eff=240 mb−270 mb, while in visibly bright regions between hot spots we obtained τ_upper=1.6-2.2 and p_eff=250 mb−290 mb. May 1996 observations yielded slightly higher values of τ_upper (1.8-2.3 and 2.0-2.8) and p_eff (250 mb−310 mb and 265 mb−320 mb). We found that the most important variable parameter is the opacity of the middle cloud, which ra nged from τ=1, 2 in dark regions, to τ=8-30 in bright regions. From limb darkening characteristics, we inferred a wavelength-dependent haze opacity ranging from 0.2±0.05 at 660 nm to 0.35±0.05 at 953 nm, and an effective haze pressure near 120 mb. We did not find it necessary to use low single scattering albedos that require effective imaginary indices, that are several orders of magnitude larger than the values of the main putative cloud components.     

Numerical models of Saturn's long-lived anticyclones

New measurements of the dynamical properties of the long-lived Saturn's anticyclonic vortex known as “Brown Spot” (BS), discovered during the Voyager 1 and 2 flybys in 1980-1981 at latitude 43.1° N, and model simulations using the EPIC code, have allowed us to constrain the vertical wind shear and static stability in Saturn's atmosphere (vertically from pressure levels from 10 mbar to 10 bars) at this latitude. BS dynamical parameters from Voyager images include its size as derived from cloud albedo gradient (6100 km East-West times 4300 km North-South), mean tangential velocity ( at 2400 km from center) and mean vorticity (4.0±1.5×¹⁰⁻⁵ s⁻¹), lifetime >1 year, drift velocity relative to Voyager's System III rotation rate, mean meridional atmospheric wind profile at cloud level at its latitude and interactions with nearby vortices (pair orbiting and merging). An extensive set of numerical experiments have been performed to try to reproduce this single vortex properties and its observed mergers with smaller anticyclones by varying the vertical structure of the zonal wind and adjusting the static stability of the lower stratosphere and upper troposphere. Within the context of the EPIC model atmosphere, our simulations indicate that BS's drift velocity, longevity and merging behavior are very sensitive to these two atmospheric properties. The best results at the BS latitude occur for static stability conditions that use a Brunt-Väisäla frequency constant in the upper troposphere (from 0.5 to 10 bar) above 3.2×¹⁰⁻³ s⁻¹ and suggest that the wind speed slightly decays below the visible cloud deck from ∼0.5 to 10 bar at a rate per scale height. Changing the vortex latitude within the band domain introduces latitude oscillations in the vortex but not a significant meridional migration. Simulated mergers always showed orbiting movements with a typical merging time of about three days, very close to the time-span observed in the interaction of real vortices. Although these results are not unique in view of the unknowns of Saturn's deep atmosphere, they serve to constrain realistically its structure for ongoing Cassini observations.     

New laboratory measurements of CH4 in Titan's conditions and a reanalysis of the DISR near-surface spectra at the Huygens landing site

Laboratory spectra of methane-nitrogen mixtures have been recorded in the near-infrared range (1.0-1.65 μm) in conditions similar to Titan's near surface, to facilitate the interpretation of the DISR/DLIS (DISR—Descent Imager/Spectral Radiometer) spectra taken during the last phase of the descent of the Huygens Probe, when the surface was illuminated by a surface-science lamp. We used a 0.03 cm⁻¹ spectral resolution, adequate to resolve the lines at high pressure (p_N2∼1.5 bar). By comparing the laboratory spectra with synthetic calculations in the well-studied ν₂+2ν₃ band (7515-7620 cm⁻¹), we determine a methane absorption column density of 178±20 cm atm and a temperature of 118±10 K in our experiment. From this, we derive the methane absorption coefficients over 1.0-1.65 μm with a 0.03 cm⁻¹ sampling, allowing for the extrapolation of the results to any other methane column density under the relevant pressure and temperature conditions. We then revisit the calibration and analysis of the Titan “lamp-on” DLIS spectra. We infer a 5.1±0.8% methane-mixing ratio in the first 25 m of Titan's atmosphere. The CH₄ mixing ratio measured 90 s after landing from a distance of 45 cm is found to be 0.92±0.25 times this value, thus showing no post-landing outgassing of methane in excess of ∼20%. Finally, we determine the surface reflectivity as seen between 25 m and 45 cm and find that the 1500 nm absorption band is deeper in the post-landing spectrum as compared to pre-landing.     

Effects of a large convective storm on Saturn's equatorial jet

Hubble Space Telescope observations revealed that Saturn's equatorial jet at the cloud level blows at ∼275 m s⁻¹ today, approximately half the ∼470 m s⁻¹ wind during the Voyager flybys in 1980-1981. Radiative transfer calculations estimate the clouds to be significantly higher today than in 1980. The higher clouds make it difficult to observationally isolate any true slowdown from the vertical wind shear because Voyager and Cassini observations show that the winds become slower with altitude. Here, we test the hypothesis that the large equatorial storm in 1990 called the Great White Spot (GWS) decelerated the equatorial jet. We first use order of magnitude estimates to show: (1) if the GWS triggers vertical momentum redistribution, a minor speed change in the troposphere can lead to a substantial stratospheric wind speed change; (2) storm-triggered turbulent mixing slows a prograde equatorial jet; and (3) a prograde equatorial jet inhibits turbulent mixing in latitude. To test whether a GWS-like large storm decelerates the equatorial jet, we perform numerical experiments using the Explicit Planetary Isentropic Coordinate (EPIC) atmosphere model. Our simulation results are consistent with our order of magnitude predictions. We show that the storm excites waves, and the waves transport westward momentum from the troposphere to the stratosphere and decelerate the equatorial jet by as much as ∼40 m s⁻¹ at the 10-mbar level. However, our results show that the storm's effect is too weak at the cloud levels to halve the jet's speed from ∼470 m s⁻¹. Our results suggest that a combination of higher clouds and a true slowdown is necessary to explain the apparent equatorial jet slowdown. We also analyze the effect of waves on the apparent cloud motions, and show that waves can influence cloud-tracking wind speed measurements.     

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.     

Laboratory kinetics of C2H radical reactions with ethane, propane, and n-butane at T = 96-296 K: implications for Titan

The kinetics of the reactions of C₂H radical with ethane (k₁), propane (k₂), and n-butane (k₃) are studied over the temperature range of T = 96-296 K with a pulsed Laval nozzle apparatus that utilizes a pulsed laser photolysis-chemiluminescence technique. The C₂H decay profiles in the presence of both the alkane reactant and O₂ are monitored by the CH(A²Δ) chemiluminescence tracer method. The results, together with available literature data, yield the following Arrhenius expressions: k₁(T) = (0.51 ± 0.06) × 10⁻¹⁰ exp[(−76 ± 30)K/T] cm³ molecule⁻¹ s⁻¹ (T = 96-800 K), k₂(T) = (0.98 ± 0.32) × 10⁻¹⁰exp[(−71 ± 60)K/T] cm³ molecule⁻¹ s⁻¹ (T = 96-361 K), and k₃(T) = (1.23 ± 0.26) × 10⁻¹⁰ cm³ molecule⁻¹ s⁻¹ (T = 96-297 K). At T = 296 K, k₁ is measured as a function of total pressure and has little or no pressure dependence. The results from this work support a direct hydrogen abstraction mechanism for the title reactions. Implications to the atmospheric chemistry of Titan are discussed.     

Keck Adaptive Optics Images of Uranus and Its Rings

We present adaptive optic images of Uranus obtained with the 10-m W. M. Keck II telescope in June 2000, at wavelengths between 1 and 2.4 μm. The angular resolution of the images is ∼0.06-0.09″. We identified eight small cloud features on Uranus's disk, four of which were in the northern hemisphere. The latter features are ∼1000-2000 km in extent and located in the upper troposphere, above the methane cloud, at pressures between 0.5 and 1 bar. Our data have been combined with HST data by Hammel et al. (2001, Icarus153, 229-235); the combination of Keck and HST data allowed derivation of an accurate wind velocity profile. Our images further show Uranus's entire ring system: the asymmetric ? ring, as well as the three groups of inner rings (outward from Uranus): the rings 6+5+4, α+β, and the η+γ+δ rings. We derived the equivalent I/F width and ring particle reflectivity for each group of rings. Typical particle albedos are ∼0.04-0.05, in good agreement with HST data at 0.9 μm.     

Titan’s 5-μm window: observations with the Very Large Telescope

We report on mid-resolution (R∼2000) spectroscopic observations of Titan, acquired in November 2000 with the Very Large Telescope and covering the range 4.75-5.07 μm. These observations provide a detailed characterization of the CO (1-0) vibrational band, clearly separating for the first time individual CO lines (P10 to P19 lines of ¹³CO). They indicate that the CO/N₂ mixing ratio in Titan’s troposphere is 32±10 ppm. Comparison with photochemical models indicates that CO is not in a steady state in Titan’s atmosphere. The observations confirm that Titan’s 5-μm continuum geometric albedo is ∼0.06, and further indicates a ∼20% albedo decrease over 4.98-5.07 μm. Nonzero flux is detected at the 0.01 geometric albedo level in the saturated core of the ¹²CO (1-0) band, at 4.75-4.85 μm, providing evidence for backscattering on the stratospheric haze. Finally, emission lines are detected at 4.75-4.835 μm, coinciding in position with lines from the CO(1-0) and/or CO(2-1) bands. Matching them by thermal emission would require Titan’s stratosphere to be much warmer (by ∼ 25 K at 0.1 mbar) than indicated by the methane 7.7-μm emission and the Voyager radio-occultation. We show instead that a nonthermal mechanism, namely solar-excited fluorescence, is a more plausible source for these emissions. Improved observations and laboratory measurements on the vibrational-translational relaxation of CO are needed for further interpretation of these emissions in terms of a CO stratospheric mixing ratio.     

A Study on the Periodicity of Hardness Ratio of the X-ray Binary Cyg X-1 as a Black Hole Candidate

Using the All-Sky Monitor (ASM, 1.5∼12 keV) data of Rossi X-ray Timing Explorer (RXTE) from January 1996 to October 2011, we have analyzed in detail the power spectrum of the hardness ratio (HR) (5∼12 keV/3∼5 keV) of the X-ray binary Cyg X-1 as a black hole candidate. The results show that the HR exhibits the following periodical variations: (1) During MJD = 50087∼55841, the HR presented the the periods T ≈ 5.6 d, T ≈ 40.0 d, T ≈ 78.4 d, T ≈ 173.8 d, and T ≈ 400/800 d; (2) When Cyg X-1 was in the hard state, the HR exhibited the periods T ≈ 5.6 d, T ≈ 33.7/67.6 d, T ≈ 45.3 d, and T ≈ 165.3 d; (3) When Cyg X-1 was in the soft state, the HR exhibited the periods T ≈ 38.5 d, T ≈ d, and T ≈ 128.3 d. Moreover, using the viscosity theory and Zdziarski accretion disk model, we have made a discussion on the physical mechanism of this kind of periodicity.     

Dust devil sediment flux on Earth and Mars: Laboratory simulations

Laboratory simulations using the Arizona State University Vortex Generator (ASUVG) were run to simulate sediment flux in dust devils in terrestrial ambient and Mars-analog conditions. The objective of this study was to measure vortex sediment flux in the laboratory to yield estimations of natural dust devils on Earth and Mars, where all parameters may not be measured. These tests used particles ranging from 2 to 2000 μm in diameter and 1300 to 4800 kg m⁻³ in density, and the results were compared with data from natural dust devils on Earth and Mars. Typically, the cores of dust devils (regardless of planetary environment) have a pressure decrease of ∼0.1-1.5% of ambient atmospheric pressure, which enhances the lifting of particles from the surface. Core pressure decreases in our experiments ranged from ∼0.01% to 5.00% of ambient pressure (10 mbar Mars cases and 1000 mbar for Earth cases) corresponding to a few tenths of a millibar for Mars cases and a few millibars for Earth cases. Sediment flux experiments were run at vortex tangential wind velocities of 1-45 m s⁻¹, which typically correspond to ∼30-70% above vortex threshold values for the test particle sizes and densities. Sediment flux was determined by time-averaged measurements of mass loss for a given vortex size. Sediment fluxes of ∼10⁻⁶-10⁰ kg m⁻² s⁻¹ were obtained, similar to estimates and measurements for fluxes in dust devils on Earth and Mars. Sediment flux is closely related to the vortex intensity, which depends on the strength of the pressure decrease in the core (ΔP). This study found vortex size is less important for lifting materials because many different diameters can have the same ΔP. This finding is critical in scaling the laboratory results to natural dust devils that can be several orders of magnitude larger than the laboratory counterparts.     

Global contraction of planetary bodies due to despinning: Application to Mercury and Iapetus

We extend previous work on the global tectonic patterns generated by despinning with a self-consistent treatment of the isotropic despinning contraction that has been ignored. We provide simple analytic approximations that quantify the effect of the isotropic despinning contraction on the global shape and tectonic pattern. The isotropic despinning contraction of Mercury is ∼93 m (T/1 day)⁻², where T is the initial rotation period. If we take into account both the isotropic contraction and the degree-2 deformations associated with despinning, the preponderance of compressional tectonic features on Mercury’s surface requires an additional isotropic contraction ?1 km (T/1 day)⁻², presumably due to cooling of the interior and growth of the solid inner core. The isotropic despinning contraction of Iapetus is ∼9 m (T/16 h)⁻², and it is not sensitive to the presence of a core or the thickness of the elastic lithosphere. The tectonic pattern expected for despinning, including the isotropic contraction, does not explain Iapetus’ ridge. Furthermore, the ridge remains unexplained with the addition of any isotropic compressional stresses, including those generating by cooling.     

Titan's aerosols: Comparison between our model and DISR findings

Our model [Dimitrov, V., Bar-Nun, A., 1999. A model of energy dependent agglomeration of hydrocarbon aerosol particles and implication to Titan's aerosol. J. Aerosol. Sci. 30(1), 35-49] describes the experimentally found polymerization of C₂H₂ and HCN to form aerosol embryos, their growth and adherence to form various aerosols objects [Bar-Nun, A., Kleinfeld, I., Ganor, E., 1988. Shape and optical properties of aerosols formed by photolysis of C₂H₂, C₂H₄ and HCN. J. Geophys. Res. 93, 8383-8387]. These loose fractal objects describe well the findings of DISR on the Huygens probe [Tomasko, M.G., Bézard, B., Doose, L., Engel, S., Karkoschka, E., 2008. Measurements of methane absorption by the descent imager/spectral radiometer (DISR) during its descent through Titan's atmosphere. Planet. Space Sci., this issue, doi:10.1016/j.pss.2007]. These include (1) various regular objects of R=(0.035-0.064)×10⁻⁶ m, as compared with DISR's 0.05×10⁻⁶ m; (2) diverse low and high fractal structures composed of random combinations of various regular and irregular objects; (3) the number density of fractal particles is 6.9×10⁶ m⁻³ at Z=100 km, as compared with DISR's finding of 5.0×10⁶ m⁻³ at Z=80 km; (4) the number of structural units per higher fractals in the atmosphere at Z∼100 km is (2400-2700), as compared with DISR's 3000, and their size being of R=(5.4-6.4)×10⁻⁶ m will satisfy this value and (5) condensation of CH₄ on the highly fractal structures could begin at the altitude where thin methane clouds were observed, filling somewhat the new open fractal structures.     

The two basic components of the interstellar neutral hydrogen and its relation with the galactic structure

The two basic components of the neutral hydrogen, cool dense clouds merged in a hotter tenuous medium, are studied using 21 cm absorption data of the Parkes Survey. The mean parameters obtained for the typical clouds next to the galactic plane are τ_p = 1.7, velocity half-width=3.3 km s^?1. Their temperatures areT _sc ≥40 K with a meanT _sc =63±12 K and the obtained hot gas density isn _HH=(0.15±0.05) atom cm^?3. Theoretical analysis following Giovanelli and Brown (1973) reveals that the pressure equilibrium condition (n _HH+2n _e )·T _SH≈n _HC·T _sc is compatible with the quoted values if it is assumed that the cosmic abundances in the interstellar medium are below the adopted normal solar abundance. This lack of heavy elements suggests accretion to grains which is consistent with the observed narrow concentration of the dark matter on the galactic layer (≤100 pc halfwidth). The same pressure condition leads to a mean cool cloud density ofn _HC～30 atom cm^?3 and a hot gas temperature ofT _SH～10 500 K. Comparison with data from Hii regions suggests that the cool clouds are somewhat denser and less extensive than such regions. An explanation for it is the expansion that the Hii regions went through in their origin. Comparison with 21 cm emission data shows that the cloud galactic layer is only about a quarter as thick as the hot gas layer. All the present results suggest that only such clouds can be spatially related with the typical I population associated with the spiral structure.     

Constraints on gas giant satellite formation from the interior states of partially differentiated satellites

An icy satellite whose interior is composed of a homogeneous ice/rock mixture must avoid melting during its entire history, including during its formation when it was heated by deposition of accretional energy and short-lived radioisotopes. Estimates of the temperature rise associated with radiogenic and accretional heating, coupled with limits on satellite melting can be used to constrain the timing of formation of a partially differentiated satellite relative to the origin of the calcium-aluminum-rich inclusions (CAI's) as a function of its accretion timescale and the protosatellite disk temperature (Td). Geological characterization and spacecraft radio tracking data suggest that Callisto, the outermost regular satellite of Jupiter, and Saturn's mid-sized satellite Rhea, are partially differentiated if their interiors are in hydrostatic equilibrium. Because the specific conditions during the satellites' formation are uncertain, we determine accretional temperature profiles for a range of values for Td and accretion time scales with the limiting assumption that impactors deposit all their energy close to the surface, leading to maximally effective radiative cooling. We find that Callisto can remain unmelted during formation if it accreted on a time scale longer than 0.6 Myr. Considering both radiogenic and accretional heating, Callisto must have finished accreting no earlier than ∼4 Myr after formation of CAI's for Td=100 K. Warmer disks or larger impactors that deposit their energy at depth in the satellite would require longer and/or later formation times. If Rhea accreted slowly (in 10⁵ to 10⁶ years), its growth must have finished no earlier than ∼2 Myr after CAI's for 70 K?Td?250 K to avoid early melting. If Rhea formed quickly (?10³ yr), its formation must have been delayed until at least 2 to 7 Myr after CAI's and in a disk with Td<190 K in the small impactor limit. If the satellites form in slow-inflow-supplied disks as proposed by Canup and Ward [Canup, R.M., Ward, W.R., 2002. Astron. J. 124, 3404-3423], the implied satellite ages suggest that gas inflow to the giant planets ceased no earlier than ∼4 Myr after CAI's, comparable to average nebular lifetimes inferred from observations of circumstellar disks.     

Saturn's magnetospheric interaction with Titan as defined by Cassini encounters T9 and T18: New results

We present new results of Cassini's T9 flyby with complementary observations from T18. Based on Cassini plasma spectrometer (CAPS) and Cassini magnetometer (MAG), compositional evidence shows the upstream flow for both T9 and T18 appears composed of light ions (H⁺ and H₂⁺), with external pressures ∼30 times lower than that for the earlier TA flyby where heavy ions dominated the magnetospheric plasma. When describing the plasma heating and sputtering of Titan's atmosphere, T9 and T18 can be considered interactions of low magnetospheric energy input. On the other hand, T5, when heavy ion fluxes are observed to be higher than typical (i.e., TA), represents the limiting case of high magnetospheric energy input to Titan's upper atmosphere. Anisotropy estimates of the upstream flow are 1<T_⊥/T_‖<3 and the flow is perpendicular to B, indicative of local picked up ions from Titan's H and H₂ coronae extending to Titan's Hill sphere radius. Beyond this distance the corona forms a neutral torus that surrounds Saturn. The T9 flyby unexpectedly resulted in observation of two “wake” crossings referred to as Events 1 and 2. Event 2 was evidently caused by draped magnetosphere field lines, which are scavenging pickup ions from Titan's induced magnetopause boundary with outward flux ∼2×10⁶ ions/cm²/s. The composition of this out flow is dominated by H₂⁺ and H⁺ ions. Ionospheric flow away from Titan with ion flux ∼7×10⁶ ion/cm²/s is observed for Event 1. In between Events 1 and 2 are high energy field aligned flows of magnetosphere protons that may have been accelerated by the convective electric field across Titan's topside ionosphere. T18 observations are much closer to Titan than T9, allowing one to probe this type of interaction down to altitudes ∼950 km. Comparisons with previously reported hybrid simulations are made.     

On the Vertical Thermal Structure of Pluto's Atmosphere

A radiative–conductive model for the vertical thermal structure of Pluto's atmosphere is developed with a non-LTE treatment of solar heating in the CH₄3.3 μm and 2.3 μm bands, non-LTE radiative exchange and cooling in the CH₄7.6 μm band, and LTE cooling by CO rotational line emission. The model includes the effects of opacity and vibrational energy transfer in the CH₄molecule. Partial thermalization of absorbed solar radiation in the CH₄3.3 and 2.3 μm bands by rapid vibrational energy transfer from the stretch modes to the bending modes generates high altitude heating at sub-microbar pressures. Heating in the 2.3 μm bands exceeds heating in 3.3 μm bands by approximately a factor of 6 and occurs predominantly at microbar pressures to generate steep temperature gradients ∼10–20 K km⁻¹forp> 2 μbar when the surface or tropopause pressure is ∼3 μbar and the CH₄mixing ratio is a constant 3%. This calculated structure may account for the “knee” in the stellar occultation lightcurve. The vertical temperature structure in the first 100 km above the surface is similar for atmospheres with Ar, CO, and N₂individually as the major constituent. If a steep temperature gradient ∼20 K km⁻¹is required near the surface or above the tropopause, then the preferred major constituent is Ar with 3% CH₄mixing ratio to attain a calculated ratio ofT/(= 3.5 K amu⁻¹) in agreement with inferred values from stellar occultation data. However, pure Ar and N₂ices at the same temperature yield an Ar vapor pressure of only ∼0.04 times the N₂vapor pressure. Alternative scenarios are discussed that may yield acceptable fits with N₂as the dominant constituent. One possibility is a 3 μbar N₂atmosphere with 0.3% CH₄that has 106 K isothermal region (T/= 3.8 K amu⁻¹) and ∼8 K km⁻¹surface/tropopause temperature gradient. Another possibility would be a higher surface pressure ∼10 μbar with a scattering haze forp> 2 μbar. Our model with appropriate adjustments in the CH₄density profile to Triton's inferred profile yields a temperature profile consistent with the UVS solar occultation data (Krasnopolsky, V. A., B. R. Sandel, and F. Herbert 1992.J. Geophys. Res.98, 3065–3078.) and ground-based stellar occultation data (Elliot, J. L., E. W. Dunham, and C. B. Olkin 1993.Bull. Am. Astron. Soc.25, 1106.).     

Baroclinic instabilities in Jupiter's zonal flow

The baroclinic stability of Jupiter's zonal flow is investigated using a model consisting of two continuously stratified fluid layers. The upper layer, containing a zonal shear flow and representing the Jovian cloudy regions above p ～ 5 bars, is the same as Eady's (1949) model for the Earth. The lower layer has a relatively large but finite depth with a quiescent basic state, representing the deep Jovian fluid bulk below p ～ 5 bars. Due to the presence of the lower layer, the linearized non-dimensional growth rates are drastically reduced from the O(1) growth rates of the original Early model. Only very long wavelengths relative to the upper fluid's radius of deformation L₁ are unstable. Eddy transports of heat are also reduced relative to estimates based on scaling arguments alone. Since the hydrostatic approximation for the lower-layer perturbation breaks down at great depths, a second model is presented in which energy propagates downward in an infinitely deep lower fluid obeying the full linearized fluid equations. In this model, the growth rates are again very small, but now all wavelengths are unstable with maximum growth rates occurring for wavelengths O(1) relative to L₁. These results illustrate the importance for the upper-layer meteorology of the interface boundary condition with the lower fluid, which is radically different from the rigid lower boundary of the Earth's troposphere.     

Revised vertical cloud structure of Uranus from UKIRT/UIST observations and changes seen during Uranus’ Northern Spring Equinox from 2006 to 2008: Application of new methane absorption data and comparison with Neptune

Long-slit spectroscopy observations of Uranus by the United Kingdom InfraRed Telescope UIST instrument in 2006, 2007 and 2008 have been used to monitor the change in Uranus’ vertical and latitudinal cloud structure through the planet’s Northern Spring Equinox in December 2007.These spectra were analysed and presented by Irwin et al. (Irwin, P.G.J., Teanby, N.A., Davis, G.R. [2009]. Icarus 203, 287-302), but since publication, a new set of methane absorption data has become available (Karkoschka, E., Tomasko, M. [2010]. Methane absorption coefficients for the jovian planets from laboratory, Huygens, and HST data. Icarus 205, 674-694.), which appears to be more reliable at the cold temperatures and high pressures of Uranus’ deep atmosphere. We have fitted k-coefficients to these new methane absorption data and we find that although the latitudinal variation and inter-annual changes reported by Irwin et al. (2009) stand, the new k-data place the main cloud deck at lower pressures (2-3 bars) than derived previously in the H-band of ∼3-4 bars and ∼3 bars compared with ∼6 bars in the J-band. Indeed, we find that using the new k-data it is possible to reproduce satisfactorily the entire observed centre-of-disc Uranus spectrum from 1 to 1.75 μm with a single cloud at 2-3 bars provided that we make the particles more back-scattering at wavelengths less than 1.2 μm by, for example, increasing the assumed single-scattering albedo from 0.75 (assumed in the J and H-bands) to near 1.0. In addition, we find that using a deep methane mole fraction of 4% in combination with the associated warm ‘F’ temperature profile of Lindal et al. (Lindal, G.F., Lyons, J.R., Sweetnam, D.N., Eshleman, V.R., Hinson, D.P. [1987]. J. Geophys. Res. 92, 14987-15001), the retrieved cloud deck using the new (Karkoschka and Tomasko, 2010) methane absorption data moves to between 1 and 2 bars.The same methane absorption data and retrieval algorithm were applied to observations of Neptune made during the same programme and we find that we can again fit the entire 1-1.75 μm centre-of-disc spectrum with a single cloud model, providing that we make the stratospheric haze particles (of much greater opacity than for Uranus) conservatively scattering (i.e. ω = 1) and we also make the deeper cloud particles, again at around the 2 bar level more reflective for wavelengths less than 1.2 μm. Hence, apart from the increased opacity of stratospheric hazes in Neptune’s atmosphere, the deeper cloud structure and cloud composition of Uranus and Neptune would appear to be very similar.