Short-term variability of Jupiter’s extended sodium nebula

Ground-based optical observations of D₁ and D₂ line emissions from Jupiter’s sodium nebula, which extend over several hundreds of jovian radii, were carried out at Mt. Haleakala, Maui, Hawaii using a wide field filter imager from May 19 to June 21, 2007. During this observation, the east-west asymmetry of the nebula with respect to the Io’s orbital motion was clearly identified. Particularly, the D₁+D₂ brightness on the western side of Jupiter is strongly controlled by the Io phase angle. The following scenario was developed to explain this phenomenon as follows: First, more ionospheric ions like NaX⁺, which are thought to produce fast neutral sodium atoms due to a dissociative recombination process, are expected to exist in Io’s dayside hemisphere rather than in the nightside one. Second, it is expected that more NaX⁺ ionospheric ions are picked up by the jovian co-rotating magnetic field when Io’s leading hemisphere is illuminated by the Sun. Third, the sodium atom ejection rate varies with respect to Io’s orbital position as a result of the first two points. Model simulations were performed using this scenario. The model results were consistent with the observation results, suggesting that Io’s ionosphere is expected to be controlled by solar radiation just like Earth.     

Jupiter’s moment of inertia: A possible determination by Juno

The moment of inertia of a giant planet reveals important information about the planet’s internal density structure and this information is not identical to that contained in the gravitational moments. The forthcoming Juno mission to Jupiter might determine Jupiter’s normalized moment of inertia NMoI = C/MR² by measuring Jupiter’s pole precession and the Lense–Thirring acceleration of the spacecraft (C is the axial moment of inertia, and M and R are Jupiter’s mass and mean radius, respectively). We investigate the possible range of NMoI values for Jupiter based on its measured gravitational field using a simple core/envelope model of the planet assuming that J₂ and J₄ are perfectly known and are equal to their measured values. The model suggests that for fixed values of J₂ and J₄ a range of NMoI values between 0.2629 and 0.2645 can be found. The Radau–Darwin relation gives a NMoI value that is larger than the model values by less than 1%. A low NMoI of ∼0.236, inferred from a dynamical model (Ward, W.R., Canup, R.M. [2006]. Astrophys. J. 640, L91–L94) is inconsistent with this range, but the range is model dependent. Although we conclude that the NMoI is tightly constrained by the gravity coefficients, a measurement of Jupiter’s NMoI to a few tenths of percent by Juno could provide an important constraint on Jupiter’s internal structure. We carry out a simplified assessment of the error involved in Juno’s possible determination of Jupiter’s NMoI.     

The formation of Uranus and Neptune in solid-rich feeding zones: Connecting chemistry and dynamics

The core accretion theory of planet formation has at least two fundamental problems explaining the origins of Uranus and Neptune: (1) dynamical times in the trans-saturnian solar nebula are so long that core growth can take >15 Myr and (2) the onset of runaway gas accretion that begins when cores reach ∼10M_⊕ necessitates a sudden gas accretion cutoff just as Uranus and Neptune’s cores reach critical mass. Both problems may be resolved by allowing the ice giants to migrate outward after their formation in solid-rich feeding zones with planetesimal surface densities well above the minimum-mass solar nebula. We present new simulations of the formation of Uranus and Neptune in the solid-rich disk of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) using the initial semimajor axis distribution of the Nice model (Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A. [2005]. Nature 435, 466-469; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R. [2005]. Nature 435, 462-465; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F. [2005]. Nature 435, 459-461), with one ice giant forming at 12 AU and the other at 15 AU. The innermost ice giant reaches its present mass after 3.8-4.0 Myr and the outermost after 5.3-6 Myr, a considerable time decrease from previous one-dimensional simulations (e.g. Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y. [1996]. Icarus 124, 62-85). The core masses stay subcritical, eliminating the need for a sudden gas accretion cutoff.Our calculated carbon mass fractions of 22% are in excellent agreement with the ice giant interior models of Podolak et al. (Podolak, M., Weizman, A., Marley, M. [1995]. Planet. Space Sci. 43, 1517-1522) and Marley et al. (Marley, M.S., Gómez, P., Podolak, M. [1995]. J. Geophys. Res. 100, 23349-23354). Based on the requirement that the ice giant-forming planetesimals contain >10% mass fractions of methane ice, we can reject any Solar System formation model that initially places Uranus and Neptune inside of Saturn’s orbit. We also demonstrate that a large population of planetesimals must be present in both ice giant feeding zones throughout the lifetime of the gaseous nebula. This research marks a substantial step forward in connecting both the dynamical and chemical aspects of planet formation. Although we cannot say that the solid-rich solar nebula model of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) gives exactly the appropriate initial conditions for planet formation, rigorous chemical and dynamical tests have at least revealed it to be a viable model of the early Solar System.     

Sedna and the Oort Cloud around a migrating Sun

Recent numerical simulations have demonstrated that the Sun’s dynamical history within the Milky Way may be much more complex than that suggested by its current low peculiar velocity (Sellwood, J.A., Binney, J.J. [2002]. Mon. Not. R. Astron. Soc. 336, 785-796; Roškar, R., Debattista, V.P., Quinn, T.R., Stinson, G.S., Wadsley, J. [2008]. Astrophys. J. 684, L79-L82). In particular, the Sun may have radially migrated through the galactic disk by up to 5-6 kpc (Roškar, R., Debattista, V.P., Quinn, T.R., Stinson, G.S., Wadsley, J. [2008]. Astrophys. J. 684, L79-L82). This has important ramifications for the structure of the Oort Cloud, as it means that the Solar System may have experienced tidal and stellar perturbations that were significantly different from its current local galactic environment. To characterize the effects of solar migration within the Milky Way, we use direct numerical simulations to model the formation of an Oort Cloud around stars that end up on solar-type orbits in a galactic-scale simulation of a Milky Way-like disk formation. Surprisingly, our simulations indicate that Sedna’s orbit may belong to the classical Oort Cloud. Contrary to previous understanding, we show that field star encounters play a pivotal role in setting the Oort Cloud’s extreme inner edge, and due to their stochastic nature this inner edge sometimes extends to Sedna’s orbit. The Sun’s galactic migration heightens the chance of powerful stellar passages, and Sedna production occurs around ∼20-30% of the solar-like stars we study. Considering the entire Oort Cloud, we find its median distance depends on the minimum galactocentric distance attained during the Sun’s orbital history. The inner edge also shows a similar dependence but with increased scatter due to the effects of powerful stellar encounters. Both of these Oort Cloud parameters can vary by an order of magnitude and are usually overestimated by an Oort Cloud formation model that assumes a fixed galactic environment. In addition, the amount of material trapped in outer Oort Cloud orbits (a > 20,000 AU) can be extremely low and may present difficulties for traditional models of Oort Cloud formation and long-period comet production.     

Dynamic features of Io's extended sodium distributions

We report on two-dimensional imaging observations of D-line emissions from the extended distribution of iogenic sodium atoms with two fields of view (±20 RJ (narrow FOV) and ±400 RJ (wide FOV)) simultaneously by using a portable small telescope or camera lens. We derived dynamic feature of the band-shaped and spray-shaped distributions near Io's orbit by means of continuous observation. The observations confirm the phenomenological behavior of the sodium cloud on two spatial scales, as previously observed by Pilcher et al. [Pilcher, C.B., Smyth, W.H., Combi, M.R., Fertel, J.H., 1984. Astrophys. J. 287, 427-444], Schneider et al. [Schneider, N.M., Trauger, J.T., Wilson, J.K., Brown, D.I., Evans, R.W., Shemansky, D.E., 1991. Science 253, 1394-1397], and Mendillo et al. [Mendillo, M., Baumgartner, J., Flynn, B., Hughes, W.S., 1990. Nature 348, 312-314]. We also confirm an elongated oval emission distribution of the sodium nebula and derivation of its detailed east-west asymmetry depending on Io's phase angle, which was first noted by Flynn et al. [Flynn, B., Mendillo, M., Baumgartner, J., 1994. J. Geophys. Res. 99, 8403-8409]. We then did model analyses to investigate the source process for sodium atoms and the dynamics behind their distribution. We conclude that the essential of molecular ion mechanisms to the band-shaped distribution is in agreement with Wilson and Schneider [Wilson, J.K., Schneider, N.M., 1999. J. Geophys. Res. 104, 16567-16583]. We differ from Wilson et al. [Wilson, J.K., Mendillo, M., Baumgartner, J., Schneider, N.M., Trauger, J.T., Flynn, B., 2002. Icarus 157, 476-489] in finding that charge exchange process contributes more to the spray-shaped distribution and sodium nebula than sputtering does. These results derived the double-peaked velocity distribution of released sodium atoms, and re-confirmed the source rates in agreement with past studies.     

UV photolysis, organic molecules in young disks, and the origin of meteoritic amino acids

The origin of complex organic molecules such as amino acids and their precursors found in meteorites and comets is unknown. Previous studies have accounted for the complex organic inventory of the Solar System by aqueous chemistry on warm meteoritic parent bodies, or by accretion of organics formed in the interstellar medium. This paper proposes a third possibility: that complex organics were created in situ by ultraviolet light from nearby O/B stars irradiating ices already in the Sun’s protoplanetary disk. If the Sun was born in a dense cluster near UV-bright stars, the flux hitting the disk from external stars could be many orders of magnitude higher than that from the Sun alone. Such photolysis of ices in the laboratory can rapidly produce amino acid precursors and other complex organic molecules. I present a simple model coupling grain growth and UV exposure in a young circumstellar disk. It is shown that the production may be sufficient to create the Solar System’s entire complex organic inventory within 10⁶ yr. Subsequent aqueous alteration on meteoritic parent bodies is not ruled out.     

Embedded star clusters and the formation of the Oort cloud: II. The effect of the primordial solar nebula

This paper deals with Oort cloud formation while the Sun was in an embedded cluster and surrounded by its primordial nebula. This work is a continuation of Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82], building on the model presented therein, and adding the aerodynamic drag and gravitational potential of the primordial solar nebula. Results are presented of numerical simulations of comets subject to the gravitational influence of the Sun, Jupiter, Saturn, star cluster and primordial solar nebula; some of the simulations included the gravitational influence of Uranus and Neptune as well. The primordial solar nebula was approximated by the minimum-mass Hayashi model [Hayashi, C., Nakozawa, K., Nakagawa, Y., 1985. In: Black, D.C., Matthews, M.S. (Eds.). Protostars and Planets II. Univ. of Arizona Press, Tucson, AZ] whose inner and outer radii have been truncated at various distances from the Sun. A comet size of 1.7 km was used for most of our simulations. In all of our simulations, the density of the primordial solar nebula decayed exponentially with an e-folding time of 2 Myr. It turns out that when the primordial solar nebula extends much beyond Saturn or Neptune, virtually no material will end up in the Oort cloud (OC) during this phase. Instead, the majority of the material will be on circular orbits inside of Jupiter if the inner edge of the disk is well inside Jupiter's orbit. If the disk's inner edge is beyond Jupiter's orbit, most comets end up on orbits in exterior mean-motion resonances with Saturn when Uranus and Neptune are not present. In those cases where the outer edge of the disk is close to Saturn or Neptune, the fraction of material that ends up in the subsequently formed OC is much less than that found in Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82] for the same cluster densities. This implies that for comets of roughly 2 km in size, the presence of the primordial solar nebula hinders OC formation. A byproduct of some of our simulations are endresults with a substantial fraction of the comets in the Uranus-Neptune scattered disk. A subsequent followup of this material is planned for the near future. In order to determine the effect of the size of the comets on OC formation efficiency, a set of runs with the same initial conditions but different cometary radii have been performed as well, from which it is determined that the threshold comet size to begin producing significant Oort clouds is roughly 20 km. This implies that the presence of the primordial solar nebula acts as a size-sorting mechanism, with large bodies unaffected by the gas drag and ending up in the OC while small bodies remain trapped in the planetary region, in the models studied.     

The source of widespread 3-μm absorption in Jupiter’s clouds: Constraints from 2000 Cassini VIMS observations

The Cassini flyby of Jupiter in 2000 provided spatially resolved spectra of Jupiter’s atmosphere using the Visual and Infrared Mapping Spectrometer (VIMS). A prominent characteristic of these spectra is the presence of a strong absorption at wavelengths from about 2.9 μm to 3.1 μm, previously noticed in a 3-μm spectrum obtained by the Infrared Space Observatory (ISO) in 1996. While Brooke et al. (Brooke, T.Y., Knacke, R.F., Encrenaz, T., Drossart, P., Crisp, D., Feuchtgruber, H. [1998]. Icarus 136, 1-13) were able to fit the ISO spectrum very well using ammonia ice as the sole source of particulate absorption, Sromovsky and Fry (Sromovsky, L.A., Fry, P.M. [2010]. Icarus 210, 211-229), using significantly revised NH₃ gas absorption models, showed that ammonium hydrosulfide (NH₄SH) provided a better fit to the ISO spectrum than NH₃, but that the best fit was obtained when both NH₃ and NH₄SH were present in the clouds. Although the large FOV of the ISO instrument precluded identification of the spatial distribution of these two components, the VIMS spectra at low and intermediate phase angles show that 3-μm absorption is present in zones and belts, in every region investigated, and both low- and high-opacity samples are best fit with a combination of NH₄SH and NH₃ particles at all locations. The best fits are obtained with a layer of small ammonia-coated particles (r ∼ 0.3 μm) overlying but often close to an optically thicker but still modest layer of much larger NH₄SH particles (r ∼ 10 μm), with a deeper optically thicker layer, which might also be composed of NH₄SH. Although these fits put NH₃ ice at pressures less than 500 mb, this is not inconsistent with the lack of prominent NH₃ features in Jupiter’s longwave spectrum because the reflectivity of the core particles strongly suppresses the NH₃ absorption features, at both near-IR and thermal wavelengths. Unlike Jupiter, Saturn lacks the broad 3-μm absorption feature, but does exhibit a small absorption near 2.965 μm, which resembles a similar jovian feature and suggests that both planets contain upper tropospheric clouds of sub-micron particles containing ammonia as a minor fraction.     

A comprehensive numerical simulation of Io’s sublimation-driven atmosphere

Io’s sublimation-driven atmosphere is modeled using the direct simulation Monte Carlo (DSMC) method. These rarefied gas dynamics simulations improve upon earlier models by using a three-dimensional domain encompassing the entire planet computed in parallel. The effects of plasma heating, planetary rotation, inhomogeneous surface frost, molecular residence time of SO₂ on the exposed (non-volatile) rocky surface, and surface temperature distribution are investigated. Circumplanetary flow is predicted to develop from the warm dayside toward the cooler nightside. Io’s rotation leads to a highly asymmetric frost surface temperature distribution (due to the frost’s high thermal inertia) which results in circumplanetary flow that is not axi-symmetric about the subsolar point. The non-equilibrium thermal structure of the atmosphere, specifically vibrational and rotational temperatures, is also examined. Plasma heating is found to significantly inflate the atmosphere on both the dayside and nightside. The plasma energy flux causes high temperatures at high altitudes but plasma energy depletion through the dense gas column above the warmest frost permits gas temperatures cooler than the surface at low altitudes. A frost map (Douté, S., Schmitt, B., Lopes-Gautier, R., Carlson, R., Soderblom, L., Shirley, J., and the Galileo NIMS Team [2001]. Icarus 149, 107-132) is used to control the sublimated flux of SO₂ which can result in inhomogeneous column densities that vary by nearly a factor of four for the same surface temperature. A short residence time for SO₂ molecules on the “rock” component is found to smooth lateral atmospheric inhomogeneities caused by variations in the surface frost distribution, creating an atmosphere that looks nearly identical to one with uniform frost coverage. A longer residence time is found to agree better with mid-infrared observations (Spencer, J.R., Lellouch, E., Richter, M.J., López-Valverde, M.A., Jessup, K.L, Greathouse, T.K., Flaud, J. [2005]. Icarus 176, 283-304) and reproduce the observed anti-jovian/sub-jovian column density asymmetry. The computed peak dayside column density for Io assuming a surface frost temperature of 115 K agrees with those suggested by Lyman-α observations (Feaga, L.M., McGrath, M., Feldman, P.D. [2009]. Icarus 201, 570-584). On the other hand, the peak dayside column density at 120 K is a factor of five larger and is higher than the upper range of observations (Jessup, K.L., Spencer, J.R., Ballester, G.E., Howell, R.R., Roesler, F., Vigel, M., Yelle, R. [2004]. Icarus 169, 197-215; Spencer et al., 2005).     

New Horizons Alice ultraviolet observations of a stellar occultation by Jupiter’s atmosphere

The Alice ultraviolet spectrograph onboard the New Horizons spacecraft observed two occultations of the bright star χ Ophiucus by Jupiter’s atmosphere on February 22 and 23, 2007 during the approach phase of the Jupiter flyby. The ingress occultation probed the atmosphere at 32°N latitude near the dawn terminator, while egress probed 18°N latitude near the dusk terminator. A detailed analysis of both the ingress and egress occultations, including the effects of molecular hydrogen, methane, acetylene, ethylene, and ethane absorptions in the far ultraviolet (FUV), constrains the eddy diffusion coefficient at the homopause level to be  cm² s⁻¹, consistent with Voyager measurements and other analyses (Festou, M.C., Atreya, S.K., Donahue, T.M., Sandel, B.R., Shemansky, D.E., Broadfoot, A.L. [1981]. J. Geophys. Res. 86, 5717-5725; Vervack Jr., R.J., Sandel, B.R., Gladstone, G.R., McConnell, J.C., Parkinson, C.D. [1995]. Icarus 114, 163-173; Yelle, R.V., Young, L.A., Vervack Jr., R.J., Young, R., Pfister, L., Sandel, B.R. [1996]. J. Geophys. Res. 101 (E1), 2149-2162). However, the actual derived pressure level of the methane homopause for both occultations differs from that derived by [Festou et al., 1981] and [Yelle et al., 1996] from the Voyager ultraviolet occultations, suggesting possible changes in the strength of atmospheric mixing with time. We find that at 32°N latitude, the methane concentration is  cm⁻³ at 70,397 km, the methane concentration is  cm⁻³ at 70,383 km, the acetylene concentration is  cm⁻³ at 70,364 km, and the ethane concentration is  cm⁻³ at 70,360 km. At 18°N latitude, the methane concentration is  cm⁻³ at 71,345 km, the methane concentration is  cm⁻³ at 71,332 km, the acetylene concentration is cm⁻³ at 71,318 km, and the ethane concentration is  cm⁻³ at 71,315 km. We also find that the H₂ occultation light curve is best reproduced if the atmosphere remains cold in the microbar region such that the base of the thermosphere is located at a lower pressure level than that determined by in situ instruments aboard the Galileo probe (Seiff, A., Kirk, D.B., Knight, T.C.D., Young, R.E., Mihalov, J.D., Young, L.A., Milos, F.S., Schubert, G., Blanchard, R.C., Atkinson, D. [1998]. J. Geophys. Res. 103 (E10), 22857-22889) - the Sieff et al. temperature profile leads to too much absorption from H₂ at high altitudes. However, this result is highly model dependent and non-unique. The observations and analysis help constrain photochemical models of Jupiter’s atmosphere.     

Galileo in-situ dust measurements in Jupiter’s gossamer rings

Galileo was the first artificial satellite to orbit Jupiter. During its late orbital mission the spacecraft made two passages through the giant planet’s gossamer ring system. The impact-ionization dust detector on board successfully recorded dust impacts during both ring passages and provided the first in-situ measurements from a dusty planetary ring. During the first passage—on 5 November 2002 while Galileo was approaching Jupiter—dust measurements were collected until a spacecraft anomaly at 2.33R_J (Jupiter radii) just 16 min after a close flyby of Amalthea put the spacecraft into a safing mode. The second ring passage on 21 September 2003 provided ring dust measurements down to about 2.5R_J and the Galileo spacecraft was destroyed shortly thereafter in a planned impact with Jupiter. In all, a few thousand dust impacts were counted with the instrument accumulators during both ring passages, but only a total of 110 complete data sets of dust impacts were transmitted to Earth. Detected particle sizes range from about 0.2 to 5 μm, extending the known size distribution by an order of magnitude towards smaller particles than previously derived from optical imaging [Showalter, M.R., de Pater, I., Verbanac, G., Hamilton, D.P., Burns, J.A., 2008. Icarus 195, 361-377; de Pater, I., Showalter, M.R., Macintosh, B., 2008. Icarus 195, 348-360]. The grain size distribution increases towards smaller particles and shows an excess of these tiny motes in the Amalthea gossamer ring compared to the Thebe ring. The size distribution for the Amalthea ring derived from our in-situ measurements for the small grains agrees very well with the one obtained from images for large grains. Our analysis shows that particles contributing most to the optical cross-section are about 5 μm in radius, in agreement with imaging results. The measurements indicate a large drop in particle flux immediately interior to Thebe’s orbit and some detected particles seem to be on highly-tilted orbits with inclinations up to 20°. Finally, the faint Thebe ring extension was detected out to at least 5R_J, indicating that grains attain higher eccentricities than previously thought. The drop interior to Thebe, the excess of submicron grains at Amalthea, and the faint ring extension indicate that grain dynamics is strongly influenced by electromagnetic forces. These findings can all be explained by a shadow resonance as detailed by Hamilton and Krüger [Hamilton, D.P., Krüger, H., 2008. Nature 453, 72-75].     

Keck adaptive optics images of Jupiter’s north polar cap and Northern Red Oval

We present observations at near-infrared wavelengths (1-5 μm) of Jupiter’s north polar region and Northern Red Oval (NN-LRS-1). The observations were taken with the near-infrared camera NIRC2 coupled to the adaptive optics system on the 10-m W.M. Keck Telescope on UT 21 August 2010. At 5-μm Jupiter’s disk reveals considerable structure, including small bright rings which appear to surround all small vortices. It is striking, though, that no such ring is seen around the Northern Red Oval. In de Pater et al. [2010a. Icarus 210, 742-762], we showed that such rings also exist around all small vortices in Jupiter’s southern hemisphere, and are absent around the Great Red Spot and Red Oval BA. We show here that the vertical structure and extent of the Northern Red Oval is very similar to that of Jupiter’s Red Oval BA. These new observations of the Northern Red Oval, therefore, support the idea of a dichotomy between small and large anticyclones, in which ovals larger than about two Rossby deformation radii do not have 5-μm bright rings. In de Pater et al. [2010a. Icarus 210, 742-762], we explained this difference in terms of the secondary circulations within the vortices. We further compare the brightness distribution of our new 5-μm images with previously published radio observations of Jupiter, highlighting the depletion of NH₃ gas over areas that are bright at 5 μm.     

Power spectral analysis of Jupiter’s clouds and kinetic energy from Cassini

We present suggestive evidence for an inverse energy cascade within Jupiter’s atmosphere through a calculation of the power spectrum of its kinetic energy and its cloud patterns. Using Cassini observations, we composed full-longitudinal mosaics of Jupiter’s atmosphere at several wavelengths. We also utilized image pairs derived from these observations to generate full-longitudinal maps of wind vectors and atmospheric kinetic energy within Jupiter’s troposphere. We computed power spectra of the image mosaics and kinetic energy maps using spherical harmonic analysis. Power spectra of Jupiter’s cloud patterns imaged at certain wavelengths resemble theoretical spectra of two-dimensional turbulence, with power-law slopes near −5/3 and −3 at low and high wavenumbers, respectively. The slopes of the kinetic energy power spectrum are also near −5/3 at low wavenumbers. At high wavenumbers, however, the spectral slopes are relatively flatter than the theoretical prediction of −3. In addition, the image mosaic and kinetic energy power spectra differ with respect to the location of the transition in slopes. The transition in slope is near planetary wavenumber 70 for the kinetic energy spectra, but is typically above 200 for the image mosaic spectra. Our results also show the importance of calculating spectral slopes from full 2D velocity maps rather than 1D zonal mean velocity profiles, since at large wavenumbers the spectra differ significantly, though at low wavenumbers, the 1D zonal and full 2D kinetic energy spectra are practically indistinguishable. Furthermore, the difference between the image and kinetic energy spectra suggests some caution in the interpretation of power spectrum results solely from image mosaics and its significance for the underlying dynamics. Finally, we also report prominent variations in kinetic energy within the equatorial jet stream that appear to be associated with the 5 μm hotspots. Other eddies are present within the flow collar of the Great Red Spot, suggesting caution when interpreting snapshots of the flow inside these features as representative of a time-averaged state.     

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.     

Accretion of Saturn’s mid-sized moons during the viscous spreading of young massive rings: Solving the paradox of silicate-poor rings versus silicate-rich moons

The origin of Saturn’s inner mid-sized moons (Mimas, Enceladus, Tethys, Dione and Rhea) and Saturn’s rings is debated. Charnoz et al. [Charnoz, S., Salmon J., Crida A., 2010. Nature 465, 752–754] introduced the idea that the smallest inner moons could form from the spreading of the rings’ edge while Salmon et al. [Salmon, J., Charnoz, S., Crida, A., Brahic, A., 2010. Icarus 209, 771–785] showed that the rings could have been initially massive, and so was the ring’s progenitor itself. One may wonder if the mid-sized moons may have formed also from the debris of a massive ring progenitor, as also suggested by Canup [Canup, R., 2010. Nature 468, 943–946]. However, the process driving mid-sized moon accretion from the icy debris disks has not been investigated in details. In particular, Canup’s (2010) model does not seem able to explain the varying silicate contents of the mid-sized moons (from 6% to 57% in mass). Here, we explore the formation of large objects from a massive ice-rich ring (a few times Rhea’s mass) and describe the fundamental properties and implications of this new process. Using a hybrid computer model, we show that accretion within massive icy rings can form all mid-sized moons from Mimas to Rhea. However in order to explain their current locations, intense dissipation within Saturn (with Q_p < 2000) is required. Our results are consistent with a satellite origin tied to the rings formation at least 2.5 Gy ago, both compatible with either a formation concurrent to Saturn or during the Late Heavy Bombardment. Tidal heating related to high-eccentricity post-accretional episodes may induce early geological activity. If some massive irregular chunks of silicates were initially present within the rings, they would be present today inside the satellites’ cores which would have accreted icy shells while being tidally expelled from the rings (via a heterogeneous accretion process). These moons may be either mostly icy, or, if they contain a significant amount of rock, already differentiated from the ice without the need for radiogenic heating. The resulting inner mid-sized moons may be significantly younger than the Solar System and a ∼1 Gyr formation delay is possible between Mimas and Rhea. The rings resulting from this process would evolve to a state compatible with current mass estimates of Saturn’s rings, and nearly devoid of silicates, apart from isolated silicate chunks coated with ice, interpreted as today Saturn’s rings’ propellers and ring-moons (like Pan or Daphnis).     

The primordial excitation and clearing of the asteroid belt—Revisited

We have performed new simulations of two different scenarios for the excitation and depletion of the primordial asteroid belt, assuming Jupiter and Saturn on initially circular orbits as predicted by the Nice Model of the evolution of the outer Solar System [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465]. First, we study the effects of sweeping secular resonances driven by the depletion of the solar nebula. We find that these sweeping secular resonances are incapable of giving sufficient dynamical excitation to the asteroids for nebula depletion timescales consistent with estimates for solar-type stars, and in addition cannot cause significant mass depletion in the asteroid belt or produce the observed radial mixing of different asteroid taxonomic types. Second, we study the effects of planetary embryos embedded in the primordial asteroid belt. These embedded planetary embryos, combined with the action of jovian and saturnian resonances, can lead to dynamical excitation and radial mixing comparable to the current asteroid belt. The mass depletion driven by embedded planetary embryos alone, even in the case of an eccentric Jupiter and Saturn, is roughly 10-20× less than necessary to explain the current mass of the main belt, and thus a secondary depletion event, such as that which occurs naturally in the Nice Model, is required. We discuss the implications of our new simulations for the dynamical and collisional evolution of the main belt.     

Atmospheric mass deposition by captured planetesimals

We examine the deposition of heavy elements in the envelope of a protoplanet growing according to the core accretion scenario of Pollack et al. [Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y., 1996. Icarus 124, 62-85]. We use their atmospheric models and the deposition rates that they computed, and we calculate the amount of heavy elements that can be dissolved in the envelope. For planetesimals composed of a mixture of water, CHON, and rock, we find that almost all of the water is dissolved in the atmosphere. A substantial amount of CHON is also dissolved but it remains sequestered at high temperatures near the core. Some fraction of the rock is also dissolved in the very high temperature region near the core envelope boundary. If this dissolved material can be mixed upward later in the planet's evolution, the resulting structure would be much closer to that determined by matching the moments of Jupiter's gravitational field.     

Determining a tilt in Titan’s north-south albedo asymmetry from Cassini images

Analysis of Titan’s hemispheric brightness asymmetry from mapped Cassini images reveals an axis of symmetry that is tilted with respect to the rotational axis of the solid body. Twenty images taken from 2004 through 2007 show a mean axial offset of 3.8 ± 0.9° relative to the solid body’s pole, directed 79 ± 24° to the west of the sub-solar longitude. These values are consistent with recent measurements of an implied atmospheric spin axis determined from isothermal mapping by [Achterberg, R.K., Conrath, B.J., Gierasch, P.J., Flasar, F.M., Nixon, C.A., 2008. Icarus 197, 549-555].     

The distributions and ages of refractory objects in the solar nebula

Refractory objects such as Calcium, Aluminum-rich Inclusions, Amoeboid Olivine Aggregates, and crystalline silicates, are found in primitive bodies throughout our Solar System. It is believed that these objects formed in the hot, inner solar nebula and were redistributed during the mass and angular momentum transport that took place during its early evolution. The ages of these objects thus offer possible clues about the timing and duration of this transport. Here we study how the dynamics of these refractory objects in the evolving solar nebula affected the age distribution of the grains that were available to be incorporated into planetesimals throughout the Solar System. It is found that while the high temperatures and conditions needed to form these refractory objects may have persisted for millions of years, it is those objects that formed in the first 10⁵ years that dominate (make up over 90%) those that survive throughout most of the nebula. This is due to two effects: (1) the largest numbers of refractory grains are formed at this time period, as the disk is rapidly drained of mass during subsequent evolution and (2) the initially rapid spreading of the disk due to angular momentum transport helps preserve this early generation of grains as opposed to later generations. This implies that most refractory objects found in meteorites and comets formed in the first 10⁵ years after the nebula formed. As these objects contained live ²⁶Al, this constrains the time when short-lived radionuclides were introduced to the Solar System to no later than 10⁵ years after the nebula formed. Further, this implies that the t=0 as defined by meteoritic materials represents at most, the instant when the solar nebula finished accreting significant amounts of materials from its parent molecular cloud.     

Gas-drag-assisted capture of Himalia's family

To elucidate the capture of Jupiter's outer moons, we reverse-evolve satellites from their present orbits to their original heliocentric paths in the presence of Jupiter's primordial circumplanetary disk (Lubow et al., 1999, Astrophys. J. 526, 1001-1012; Canup and Ward, 2003, Astron. J. 124, 3404-3423). Our orbital histories use a symplectic integrator that allows dissipation. We assume that the present satellites Himalia, Elara, Lysithea, Leda, and S/2000 J11 are collisional fragments of a single parent. Our simulations show that this “prograde-cluster progenitor” (PCP) could be derived from objects with heliocentric orbits like those of the Hilda asteroid group. We show analytically that this capture is energetically possible. We also compare the spectroscopic characteristics of the prograde cluster members (Grav et al., 2003, Icarus, submitted for publication) with those of the Hildas, and conclude that the surface color of the prograde-cluster progenitor is consistent with an origin within the Hilda group. Accordingly, gas drag in the primordial jovian nebula is found to offer a plausible explanation for the origin of the prograde cluster. A similar capture mechanism is proposed for Saturn's Phoebe.