Corona-like atmospheric escape from KBOs: I. Gas dynamics

We show that for low temperatures (T∼30 K) and small, but non-negligible, gravitational fields the hydrodynamic escape of gas can be treated by Parker's theory of coronal expansion [Parker, E.N., 1963. Interplanetary Dynamical Processes. Interscience Publishers, New York]. We apply this theory to gas escape from Kuiper belt objects. We derive limits on the density and radius of the bodies for which this theory is applicable, and show how the flow depends on the mean molecular weight and internal degrees of freedom of the gas molecules. We use these results to explain the CH₄ dichotomy seen on KBOs [Schaller, E.L., Brown, M.E., 2007. Astrophys. J., 659, L61-L64].     

Effects of mass loss for highly-irradiated giant planets

We present calculations for the evolution and surviving mass of highly-irradiated extrasolar giant planets (EGPs) at orbital semimajor axes ranging from 0.023 to 0.057 AU using a generalized scaled theory for mass loss, together with new surface-condition grids for hot EGPs and a consistent treatment of tidal truncation. Theoretical estimates for the rate of energy-limited hydrogen escape from giant-planet atmospheres differ by two orders of magnitude, when one holds planetary mass, composition, and irradiation constant. Baraffe et al. [Baraffe, I., Selsis, F., Chabrier, G., Barman, T.S., Allard, F., Hauschildt, P.H., Lammer, H., 2004. Astron. Astrophys. 419, L13-L16] predict the highest rate, based on the theory of Lammer et al. [Lammer, H., Selsis, F., Ribas, I., Guinan, E.F., Bauer, S.J., Weiss, W.W., 2003. Astrophys. J. 598, L121-L124]. Scaling the theory of Watson et al. [Watson, A.J., Donahue, T.M., Walker, J.C.G., 1981. Icarus 48, 150-166] to parameters for a highly-irradiated exoplanet, we find an escape rate ∼¹⁰² lower than Baraffe's. With the scaled Watson theory we find modest mass loss, occurring early in the history of a hot EGP. In this theory, mass loss including the effect of Roche-lobe overflow becomes significant primarily for masses below a Saturn mass, for semimajor axes ?0.023 AU. This contrasts with the Baraffe model, where hot EGPs are claimed to be remnants of much more massive bodies, originally several times Jupiter and still losing substantial mass fractions at present.     

Trapping of N2, CO and Ar in amorphous ice—Application to comets

Recent attempts using high resolution spectra to detect N⁺₂ in several comets were unsuccessful [Cochran, A.L., Cochran, W.D., Baker, E.S., 2000. Icarus 146, 583-593; Cochran, A.L., 2002. Astrophys. J. 576, L165-L168]. The upper limits on N⁺₂ in comparison with the positively detected CO⁺ for Comets C/1995 O1 Hale-Bopp, 122P/1995 S1 de Vico and 153P/2002 C1 Ikeya-Zhang range between . Ar was not detected in three recent comets [Weaver, H.A., Feldman, P.D., Combi, M.R., Krasnopolsky, V., Lisse, C.M., Shemansky, D.E., 2002. Astrophys. J. 576, L95-L98], with upper limits of Ar/CO<(3.4-7.8)×¹⁰⁻² for Comets C/1999 T1 McNaught-Hartley, C/2001 A2 LINEAR and C/2000 WM1 LINEAR. The Ar detected by Stern et al. [Stern, S.A., Slater, D.C., Festou, M.C., Parker, J.Wm., Gladstone, G.R., A'Hearn, M.F., Wilkinson, E., 2000. Astrophys. J. 544, L169-L172] for Comet C/1995 O1 Hale-Bopp, gives a ratio Ar/CO=7.25×¹⁰⁻², which was not confirmed by Cosmovici et al. [Cosmovici, C.B., Bratina, V., Schwarz, G., Tozzi, G., Mumma, M.J., Stalio, R., 2006. Astrophys. Space Sci. 301, 135-143]. Trying to solve the two problems, we studied experimentally the trapping of N₂+CO+Ar in amorphous water ice, at 24-30 K. CO was found to be trapped in the ice 20-70 times more efficiently than N₂ and with the same efficiency as Ar. The resulting Ar/CO ratio of 1.2×¹⁰⁻² is consistent with Weaver et al.'s [Weaver, H.A., Feldman, P.D., Combi, M.R., Krasnopolsky, V., Lisse, C.M., Shemansky, D.E., 2002. Astrophys. J. 576, L95-L98] non-detection of Ar. However, with an extreme starting value for N₂/CO = 0.22 in the region where the ice grains which agglomerated to produce comet nuclei were formed, the expected N₂/CO ratio in the cometary ice should be 6.6×¹⁰⁻³, much higher than its non-detection limit.     

Titan's 3-micron spectral region from ISO high-resolution spectroscopy

The near-infrared spectrum of Titan, Saturn's largest moon and one of the Cassini/Huygens' space mission primary targets, covers the 0.8 to 5 micron region in which it shows several weak CH₄ absorption regions, and in particular one centered near 2.75 micron. Due to the interference of telluric absorption, only part of this window region (2.9-3.1 μm) has previously been observed from the ground [Noll, K.S., Geballe, T.R., Knacke, R., Pendleton, F., Yvonne, J., 1996. Icarus 124, 625-631; Griffith, C.A., Owen, T., Miller, G.A., Geballe, T., 1998. Nature 395, 575-578; Griffith, C.A., Owen, T., Geballe, T.R., Rayner, J., Rannou, P., 2003. Science 300, 628-630; Geballe, T.R., Kim, S.J., Noll, K.S., Griffith, C.A., 2003. Astrophys. J. 583, L39-L42]. We report here on the first spectroscopic observations of Titan covering the whole 2.4-4.9 μm region by two instruments on board the Infrared Space Observatory (ISO) in 1997. These observations show the 2.75-μm window in its complete extent for the first time. In this study we have also used a high-resolution Titan spectrum in the 2.9-3.6 μm region taken with the Keck [Geballe, T.R., Kim, S.J., Noll, K.S., Griffith, C.A., 2003. Astrophys. J. 583, L39-L42; Kim, S.J., Geballe, T.R., Noll, K.S., Courtin, R., 2005. Icarus 173, 522-532] to infer information on the atmospheric parameters (haze extinction, single scattering albedo, methane abundance, etc.) by fitting the methane bands with a detailed microphysical model of Titan's atmosphere (updated from Rannou, P., McKay, C.P., Lorenz, R.D., 2003. Planet. Space Sci. 51, 963-976). We have included in this study an updated version of a database for the CH₄ absorption coefficients [STDS, Wenger, Ch., Champion, J.-P., 1998. J. Quant. Spectrosc. Radiat. Transfer 59, 471-480. See also http://www.u-bourgogne.fr/LPUB/TSM/sTDS.html for latest updates; Boudon, V., Champion, J.-P., Gabard, T., Loëte, M., Michelot, F., Pierre, G., Rotger, M., Wenger, Ch., Rey, M., 2004. J. Mol. Spectrosc. 228, 620-634]. For the atmosphere we find that (a) the haze extinction profile that best matches the data is one with higher (by 40%) extinction in the atmosphere with respect to Rannou et al. (2003) down to about 30 km where a complete cut-off occurs; (b) the methane mixing ratio at Titan's surface cannot exceed 3% on a disk-average basis, yielding a maximum CH₄ column abundance of 2.27 km-am in Titan's atmosphere. From the derived surface albedo spectrum in the 2.7-3.08 micron region, we bring some constraints on Titan's surface composition. The albedo in the center of the methane window varies from 0.01 to 0.08. These values, compared to others reported in the other methane windows, show a strong compatibility with the water ice spectrum in the near-infrared. Without confirming its existence from this work alone, our data then appear to be compatible with water ice. A variety of other ices, such as CO₂, NH₃, tholin material or hydrocarbon liquid cannot be excluded from our data, but an additional unidentified component with a signature around 2.74 micron is required to satisfy the data.     

Moonlet induced wakes in planetary rings: Analytical model including eccentric orbits of moon and ring particles

Saturn’s rings host two known moons, Pan and Daphnis, which are massive enough to clear circumferential gaps in the ring around their orbits. Both moons create wake patterns at the gap edges by gravitational deflection of the ring material (Cuzzi, J.N., Scargle, J.D. [1985]. Astrophys. J. 292, 276-290; Showalter, M.R., Cuzzi, J.N., Marouf, E.A., Esposito, L.W. [1986]. Icarus 66, 297-323). New Cassini observations revealed that these wavy edges deviate from the sinusoidal waveform, which one would expect from a theory that assumes a circular orbit of the perturbing moon and neglects particle interactions. Resonant perturbations of the edges by moons outside the ring system, as well as an eccentric orbit of the embedded moon, may partly explain this behavior (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S., Burns, J.A., Dones, L. [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S. [2009]. Astron. J. 138, 272-286). Here we present an extended non-collisional streamline model which accounts for both effects. We describe the resulting variations of the density structure and the modification of the nonlinearity parameter q. Furthermore, an estimate is given for the applicability of the model. We use the streamwire model introduced by Stewart (Stewart, G.R. [1991]. Icarus 94, 436-450) to plot the perturbed ring density at the gap edges.We apply our model to the Keeler gap edges undulated by Daphnis and to a faint ringlet in the Encke gap close to the orbit of Pan. The modulations of the latter ringlet, induced by the perturbations of Pan (Burns, J.A., Hedman, M.M., Tiscareno, M.S., Nicholson, P.D., Streetman, B.J., Colwell, J.E., Showalter, M.R., Murray, C.D., Cuzzi, J.N., Porco, C.C., and the Cassini ISS team [2005]. Bull. Am. Astron. Soc. 37, 766), can be well described by our analytical model. Our analysis yields a Hill radius of Pan of 17.5 km, which is 9% smaller than the value presented by Porco (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236), but fits well to the radial semi-axis of Pan of 17.4 km. This supports the idea that Pan has filled its Hill sphere with accreted material (Porco, C.C., Thomas, P.C., Weiss, J.W., Richardson, D.C. [2007]. Science 318, 1602-1607). A numerical solution of a streamline is used to estimate the parameters of the Daphnis-Keeler gap system, since the close proximity of the gap edge to the moon induces strong perturbations, not allowing an application of the analytic streamline model. We obtain a Hill radius of 5.1 km for Daphnis, an inner edge variation of 8 km, and an eccentricity for Daphnis of 1.5 × 10⁻⁵. The latter two quantities deviate by a factor of two from values gained by direct observations (Jacobson, R.A., Spitale, J., Porco, C.C., Beurle, K., Cooper, N.J., Evans, M.W., Murray, C.D. [2008]. Astron. J. 135, 261-263; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767), which might be attributed to the neglect of particle interactions and vertical motion in our model.     

Titan's hydrodynamically escaping atmosphere

The upper atmosphere of Titan is currently losing mass at a rate , by hydrodynamic escape as a high density, slow outward expansion driven principally by solar UV heating by CH₄ absorption. The hydrodynamic mass loss is essentially CH₄ and H₂ escape. Their combined escape rates are restricted by power limitations from attaining their limiting rates (and limiting fluxes). Hence they must exhibit gravitational diffusive separation in the upper atmosphere with increasing mixing ratios to eventually become major constituents in the exosphere. A theoretical model with solar EUV heating by N₂ absorption balanced by HCN rotational line cooling in the upper thermosphere yields densities and temperatures consistent with the Huygens Atmospheric Science Investigation (HASI) data [Fulchignoni, M., and 42 colleagues, 2005. Nature 438, 785-791], with a peak temperature of ∼185-190 K between 3500-3550 km. This model implies hydrodynamic escape rates of and , or some other combination with a higher H₂ escape flux, much closer to its limiting value, at the expense of a slightly lower CH₄ escape rate. Nonthermal escape processes are not required to account for the loss rates of CH₄ and H₂, inferred by the Cassini Ion Neutral Mass Spectrometer (INMS) measurements [Yelle, R.V., Borggren, N., de la Haye, V., Kasprzak, W.T., Niemann, H.B., Müller-Wodarg, I., Waite Jr., J.H., 2006. Icarus 182, 567-576].     

Simulations of planet migration driven by planetesimal scattering

Evidence has mounted for some time that planet migration is an important part of the formation of planetary systems, both in the Solar System [Malhotra, R., 1993. Nature 365, 819-821] and in extrasolar systems [Mayor, M., Queloz, D., 1995. Nature 378, 355-359; Lin, D.N.C., Bodenheimer, P., Richardson, D.C., 1996. Nature 380, 606-607]. One mechanism that produces migration (the change in a planet's semi-major axis a over time) is the scattering of comet- and asteroid-size bodies called planetesimals [Fernandez, J.A., Ip, W.-H., 1984. Icarus 58, 109-120]. Significant angular momentum exchange can occur between the planets and the planetesimals during local scattering, enough to cause a rapid, self-sustained migration of the planet [Ida, S., Bryden, G., Lin, D.N.C., Tanaka, H., 2000. Astrophys. J. 534, 428-445]. This migration has been studied for the particular case of the four outer planets of the Solar System (as in Gomes et al. [Gomes, R.S., Morbidelli, A., Levison, H.F., 2004. Icarus 170, 492-507]), but is not well understood in general. We have used the Miranda [McNeil, D., Duncan, M., Levison, H.F., 2005. Astron. J. 130, 2884-2899] computer simulation code to perform a broad parameter-space survey of the physical variables that determine the migration of a single planet in a planetesimal disk. Migration is found to be predominantly inwards, and the migration rate is found to be independent of planet mass for low-mass planets in relatively high-mass disks. Indeed, a simple scaling relation from Ida et al. [Ida, S., Bryden, G., Lin, D.N.C., Tanaka, H., 2000. Astrophys. J. 534, 428-445] matches well with the dependencies of the migration rate:

(1)     

Ion irradiation of crystalline H2O-ice: Effect on the 1.65-μm band

We have found that 0.8 MeV proton irradiation of crystalline H₂O-ice results in temperature dependent amorphization. The H₂O-ice's phase was determined using the near infrared spectrum from 1.0 μm (10,000 cm⁻¹) to 2.5 μm (4000 cm⁻¹). In crystalline H₂O-ice, the 1.65-μm (6061 cm⁻¹) band is strong while it is nearly absent in the amorphous spectrum [Schmitt, B., Quirico, E., Trotta, F., Grundy, W.M., 1998. In: Schmitt, B., de Bergh, C., Festou, M. (Eds.), Solar System Ices. Kluwer Academic, Norwell, MA, 1998, pp. 199-240]. In this experiment, at low temperatures (9, 25, and 40 K), irradiation of crystalline H₂O-ice produced the amorphous H₂O-ice's spectrum. However, at 50 K, some crystalline absorptions persisted after irradiation and at 70 and 100 K the crystalline spectrum showed only slight changes after irradiation. Our results agree with previous H₂O-ice irradiation studies examining the crystalline peaks near 44 and 62 μm by Moore and Hudson [Moore, M.H., Hudson, R.L., 1992. Astrophys. J. 401, 353-360] and near 3.07 μm by Strazzulla et al. [Strazzulla, G., Baratta, G.A., Leto, G., Foti, G., 1992. Europhys. Lett. 18, 517-522] and by Leto and Baratta [Leto, G., Baratta, G.A., 2003. Astron. Astrophys. 397, 7-13]. We present a method of measuring band areas to quantify the phase and radiation dose of icy Solar System surfaces.     

Accretion of Jupiter’s atmosphere from a supernova-contaminated molecular cloud

If Jupiter and the Sun both formed directly from the same well-mixed proto-solar nebula, then their atmospheric compositions should be similar. However, direct sampling of Jupiter’s troposphere indicates that it is enriched in elements such as C, N, S, Ar, Kr, and Xe by 2-6× relative to the Sun (Wong, M.H., Lunine, J.I., Atreya, S.K., Johnson, T., Mahaffy, P.R., Owen, T.C., Encrenaz, T. [2008]. 219-246). Most existing models to explain this enrichment require an extremely cold proto-solar nebula which allows these heavy elements to condense, and cannot easily explain the observed variations between these species. We find that Jupiter’s atmospheric composition may be explained if the Solar System’s disk heterogeneously accretes small amounts of enriched material such as supernova ejecta from the interstellar medium during Jupiter’s formation. Our results are similar to, but substantially larger than, isotopic anomalies in terrestrial material that indicate the Solar System formed from multiple distinct reservoirs of material simultaneously with one or more nearby supernovas (Trinquier, A., Birck, J.-L., Allegre, C.J. [2007]. Astrophys. J. 655, 1179-1185). Such temporal and spatial heterogeneities could have been common at the time of the Solar System’s formation, rather than the cloud having a purely well-mixed ‘solar nebula’ composition.     

Gemini observations of disks and jets in young stellar objects and in active galaxies

We present first results from the Near-infrared Integral Field Spectrograph (NIFS) located at Gemini North. For the active galaxies Cygnus A and Perseus A we observe rotationally-supported accretion disks and adduce the existence of massive central black holes and estimate their masses. In Cygnus A we also see remarkable high-excitation ionization cones dominated by photoionization from the central engine. In the T-Tauri stars HV Tau C and DG Tau we see highly-collimated bipolar outflows in the [Fe II] λ 1.644 micron line, surrounded by a slower molecular bipolar outflow seen in the H₂ lines, in accordance with the model advocated by Pyo, T.-S., et al., Astrophys. J. 570, 724 (2002).     

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.     

The role of ejecta in the small crater populations on the mid-sized saturnian satellites

We find evidence, by both observation and analysis, that primary crater ejecta play an important role in the small crater (less than a few km) populations on the saturnian satellites, and more broadly, on cratered surfaces throughout the Solar System. We measure crater populations in Cassini images of Enceladus, Rhea, and Mimas, focusing on image data with scales less than 500 m/pixel. We use recent updates to crater scaling laws and their constants (Housen, K.R., Holsapple, K.A. [2011]. Icarus 211, 856–875) to estimate the amount of mass ejected in three different velocity ranges: (i) greater than escape velocity, (ii) less than escape velocity and faster than the minimum velocity required to make a secondary crater (v_min), and (iii) velocities less than v_min. Although the vast majority of mass on each satellite is ejected at speeds less than v_min, our calculations demonstrate that the differences in mass available in the other two categories should lead to observable differences in the small crater populations; the predictions are borne out by the measurements we have made to date. In particular, Rhea, Tethys, and Dione have sufficient surface gravities to retain ejecta moving fast enough to make secondary crater populations. The smaller satellites, such as Enceladus but especially Mimas, are expected to have little or no traditional secondary populations because their escape velocities are near the threshold velocity necessary to make a secondary crater. Our work clarifies why the Galilean satellites have extensive secondary crater populations relative to the saturnian satellites. The presence, extent, and sizes of sesquinary craters (craters formed by ejecta that escape into temporary orbits around Saturn before re-impacting the surface, see Dobrovolskis, A.R., Lissauer, J.J. [2004]. Icarus 169, 462–473; Alvarellos, J.L., Zahnle, K.J., Dobrovolskis, A.R., Hamill, P. [2005]. Icarus 178, 104–123; Zahnle, K., Alvarellos, J.L., Dobrovolskis, A.R., Hamill, P. [2008]. Icarus 194, 660–674) is not yet well understood. Finally, our work provides further evidence for a “shallow” size–frequency distribution (slope index of ～2 for a differential power-law) for comets a few kilometers diameter and smaller.     

The ultraviolet reflectance of Enceladus: Implications for surface composition

The reflectance of Saturn’s moon Enceladus has been measured at far ultraviolet (FUV) wavelengths (115-190 nm) by Cassini’s Ultraviolet Imaging Spectrograph (UVIS). At visible and near infrared (VNIR) wavelengths Enceladus’ reflectance spectrum is very bright, consistent with a surface composed primarily of H₂O ice. At FUV wavelengths, however, Enceladus is surprisingly dark - darker than would be expected for pure water ice. Previous analyses have focused on the VNIR spectrum, comparing it to pure water ice (Cruikshank, D.P., Owen, T.C., Dalle Ore, C., Geballe, T.R., Roush, T.L., de Bergh, C., Sandford, S.A., Poulet, F., Benedix, G.K., Emery, J.P. [2005] Icarus, 175, 268-283) or pure water ice plus a small amount of NH₃ (Emery, J.P., Burr, D.M., Cruikshank, D.P., Brown, R.H., Dalton, J.B. [2005] Astron. Astrophys., 435, 353-362) or NH₃ hydrate (Verbiscer, A.J., Peterson, D.E., Skrutskie, M.F., Cushing, M., Helfenstein, P., Nelson, M.J., Smith, J.D., Wilson, J.C. [2006] Icarus, 182, 211-223). We compare Enceladus’ FUV spectrum to existing laboratory measurements of the reflectance spectra of candidate species, and to spectral models. We find that the low FUV reflectance of Enceladus can be explained by the presence of a small amount of NH₃ and a small amount of a tholin in addition to H₂O ice on the surface. The presence of these three species (H₂O, NH₃, and a tholin) appears to satisfy not only the low FUV reflectance and spectral shape, but also the middle-ultraviolet to visible wavelength brightness and spectral shape. We expect that ammonia in the Enceladus plume is transported across the surface to provide a global coating.     

Embedded star clusters and the formation of the Oort Cloud

Observations suggest most stars originate in clusters embedded in giant molecular clouds [Lada, C.J., Lada, E.A., 2003. Annu. Rev. Astron. Astrophys. 41, 57-115]. Our Solar System likely spent 1-5 Myrs in such regions just after it formed. Thus the Oort Cloud (OC) possibly retains evidence of the Sun's early dynamical history and of the stellar and tidal influence of the cluster. Indeed, the newly found objects (90377) Sedna and 2000 CR₁₀₅ may have been put on their present orbits by such processes [Morbidelli, A., Levison, H.F., 2004. Astron. J. 128, 2564-2576]. Results are presented here of numerical simulations of the orbital evolution of comets subject to the influence of the Sun, Jupiter and Saturn (with their current masses on orbits appropriate to the period before the Late Heavy Bombardment (LHB) [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461]), passing stars and tidal force associated with the gas and stars of an embedded star cluster. The cluster was taken to be a Plummer model with 200-400 stars, with a range of initial central densities. The Sun's orbit was integrated in the cluster potential together with Jupiter and Saturn and the test particles. Stellar encounters were incorporated by directly integrating the effects of stars passing within a sphere centred on the Sun of radius equal to the Plummer radius for low-density clusters and half a Plummer radius for high-density clusters. The gravitational influence of the gas was modeled using the tidal force of the cluster potential. For a given solar orbit, the mean density, 〈ρ〉, was computed by orbit-averaging the density of material encountered. This parameter proved to be a good measure for predicting the properties of the OC. On average 2-18% of our initial sample of comets end up in the OC after 1-3 Myr. A comet is defined to be part of the OC if it is bound and has q>35 AU. Our models show that the median distance of an object in the OC scales approximately as 〈ρ〉^−1/2 when . Our models easily produce objects on orbits like that of (90377) Sedna [Brown, M.E., Trujillo, C., Rabinowitz, D., 2004. Astrophys. J. 617, 645-649] within ∼1 Myr in cases where the mean density is or higher; one needs mean densities of order to create objects like 2000 CR₁₀₅ by this mechanism, which are reasonable (see, e.g., Guthermuth, R.A., Megeath, S.T., Pipher, J.L., Williams, J.P., Allen, L.E., Myers, P.C., Raines, S.N., 2005. Astrophys. J. 632, 397-420). Thus the latter object may also be part of the OC. Close stellar passages can stir the primordial Kuiper Belt to sufficiently high eccentricities (e?0.05; Kenyon, S.J., Bromley, B.C., 2002. Astron. J. 123, 1757-1775) that collisions become destructive. From the simulations performed it is determined that there is a 50% or better chance to stir the primordial Kuiper Belt to eccentricities e?0.05 at 50 AU when . The orbit of the new object 2003 UB₃₁₃ [Brown, M.E., Trujillo, C.A., Rabinowitz, D.L., 2005. Astrophys. J. 635, L97-L100] is only reproduced for mean cluster densities of the order of , but in the simulations it could not come to be on its current orbit by this mechanism without disrupting the formation of bodies in the primordial Kuiper Belt down to 20 AU. It is therefore improbable that the latter object is created by this mechanism.     

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.     

Redetermination of the space weathering rate using spectra of Iannini asteroid family members

We have obtained moderate S/N (∼85) spectra at a realized resolution of R∼100 for 11 members of the Iannini family, until recently the youngest known family at under 5 million years of age [Nesvorný, D., Bottke, W.F., Levison, H.F., Dones, L., 2003. Astrophys. J. 591, 486-497, 720-771]. The spectra were acquired using the Echellette Spectrograph and Imager in its low-resolution prism mode on the Keck II telescope. The family members belong to the S-complex of asteroids with perhaps some K class members. The Iannini family members' average spectral slope, defined as the slope of the best-fit line constrained to pivot about 1 at 550 nm, is (0.30±0.04)/μm, matching the (0.26±0.03)/μm reported by Jedicke et al. [Jedicke, R., Nesvorný, D., Whiteley, R.J., Ivezi?, ?., Juri?, M., 2004. Nature 429, 275-277] using SDSS [Ivezi?, ?., Juri?, M., Lupton, R.H., Tabachnik, S., Quinn, T., 2002. In: Tyson, J.A., Wolff, S. (Eds.), Survey and Other Telescope Technologies and Discoveries. In: Proc. SPIE, vol. 4836. SPIE, Bellingham, pp. 98-103] color photometry. Using our spectra for this family as well as new observations of Karin family members [Vernazza, P., Birlan, M., Rossi, A., Dotto, E., Nesvorný, D., Brunetto, R., Fornasier, S., Fulchignoni, M., Renner, S., 2006. Astron. Astrophys. 460, 945-951] and new classifications of some older families we have revised the space weathering rate of S-complex asteroids originally determined by Jedicke et al. [Jedicke, R., Nesvorný, D., Whiteley, R.J., Ivezi?, ?., Juri?, M., 2004. Nature 429, 275-277]. Following Jedicke et al. [Jedicke, R., Nesvorný, D., Whiteley, R.J., Ivezi?, ?., Juri?, M., 2004. Nature 429, 275-277] we parameterize the space weathering rate of the principal component color of the spectrum (PC₁), which is correlated with the spectral slope, as PC₁(t)=PC₁(0)+ΔPC₁[1−exp−α(t/τ)]. Our revised rate suggests that the characteristic time scale for space weathering is τ=570±220 Myr and that new S-complex clusters will have an initial color of PC₁(0)=0.31±0.04. The revised time scale is in better agreement with lab measurements and our measurements support the use of space weathering as a dating method. Under the assumption that all the spectra should be identical, since the members all derive from the same parent body and are presumably covered with similar regolith, we combined them to obtain a high-S/N composite spectrum for the family. The combined spectrum is within the S-complex.     

A laboratory model of Saturn’s North Polar Hexagon

A hexagonal structure has been observed at ∼76°N on Saturn since the 1980s (Godfrey, D.A. [1988]. Icarus 76, 335-356). Recent images by Cassini (Baines, K., Momary, T., Roos-Serote, M., Atreya, S., Brown, R., Buratti, B., Clark, R., Nicholson, P. [2007]. Geophys. Res. Abstr. 9, 02109; Baines, K., Momary, T., Fletcher, L., Kim, J., Showman, A., Atreya, S., Brown, R., Buratti, B., Clark, R., Nicholson, P. [2009]. Geophys. Res. Abstr. 11, 3375) have shown that the feature is still visible and largely unchanged. Its long lifespan and geometry has puzzled the planetary physics community for many years and its origin remains unclear. The measured rotation rate of the hexagon may be very close to that of the interior of the planet (Godfrey, D.A. [1990]. Science 247, 1206-1208; Caldwell, J., Hua, X., Turgeon, B., Westphal, J.A., Barnet, C.D. [1993]. Science 206, 326-329; Sánchez-Lavega, A., Lecacheux, J., Colas, F., Laques, P. [1993]. Science 260, 329-332), leading to earlier interpretations of the pattern as a stationary planetary wave, continuously forced by a nearby vortex (Allison, M., Godfrey, D.A., Beebe, R.F. [1990]. Science 247, 1061-1063). Here we present an alternative explanation, based on an analysis of both spacecraft observations of Saturn and observations from laboratory experiments where the instability of quasi-geostrophic barotropic (vertically uniform) jets and shear layers is studied. We also present results from a barotropic linear instability analysis of the saturnian zonal wind profile, which are consistent with the presence of the hexagon in the North Pole and absence of its counter-part in the South Pole. We propose that Saturn’s long-lived polygonal structures correspond to wavemodes caused by the nonlinear equilibration of barotropically unstable zonal jets.     

Anelastic convection-driven dynamo benchmarks: A finite element model

We present a fully non-linear, anelastic convection and dynamo model in a rotating spherical shell that is based on a finite element method. It is shown that both the hydrodynamic steadily drifting convection and convection-driven steadily drifting dynamo are in satisfactory agreement with the anelastic convection and dynamo benchmark based on the pseudo-spectral method (Jones, C.A., Boronski, P., Brun, A.S., Glatzmaier, G.A., Gastine, T., Miesch, M.S., Wicht, J. [2011]. Icarus, 216, 120–135).     

Models of the collisional damping scenario for ice-giant planets and Kuiper belt formation

Chiang et al. [Chiang, E., Lithwick, Y., Murray-Clay, R., Buie, M., Grundy, W., Holman, M., 2007. In: Protostars and Planets V, pp. 895-911] have recently proposed that the observed structure of the Kuiper belt could be the result of a dynamical instability of a system of ∼5 primordial ice-giant planets in the outer Solar System. According to this scenario, before the instability occurred, these giants were growing in a highly collisionally damped environment according to the arguments in Goldreich et al. [Goldreich, P., Lithwick, Y., Sari, R., 2004. Astrophys. J. 614, 497-507; Annu. Rev. Astron. Astrophys. 42, 549-601]. Here we test this hypothesis with a series of numerical simulations using a new code designed to incorporate the dynamical effects of collisions. We find that we cannot reproduce the observed Solar System. In particular, Goldreich et al. [Goldreich, P., Lithwick, Y., Sari, R., 2004. Astrophys. J. 614, 497-507; Annu. Rev. Astron. Astrophys. 42, 549-601] and Chiang et al. [Chiang, E., Lithwick, Y., Murray-Clay, R., Buie, M., Grundy, W., Holman, M., 2007. In: Protostars and Planets V, pp. 895-911] argue that during the instability, all but two of the ice giants would be ejected from the Solar System by Jupiter and Saturn, leaving Uranus and Neptune behind. We find that ejections are actually rare and that instead the systems spread outward. This always leads to a configuration with too many planets that are too far from the Sun. Thus, we conclude that both Goldreich et al.'s scheme for the formation of Uranus and Neptune and Chiang et al.'s Kuiper belt formation scenario are not viable in their current forms.     

The variation of Io's auroral footprint brightness with the location of Io in the plasma torus

Ultraviolet and near-infrared observations of auroral emissions from the footprint of Io's magnetic Flux Tube (IFT) mapping to Jupiter's ionosphere have been interpreted via a combination of the unipolar inductor model [Goldreich, P., Lynden-Bell, D., 1969. Astrophys. J. 156, 59-78] and the multiply-reflected Alfvén wave model [Belcher, J.W., 1987. Science 238, 170-176]. While both models successfully explain the general nature of the auroral footprint and corotational wake, and both predict the presence of multiple footprints, the details of the interaction near Io are complicated [Saur, J., Neubauer, F.M., Connerney, J.E.P., Zarka, P., Kivelson, M.G., 2004. In: Bagenal, F., Dowling, T.E., McKinnon, W.B. (Eds.), Jupiter: The Planet, Satellites and Magnetosphere. Cambridge University Press, Cambridge, UK, pp. 537-560; Kivelson, M.G., Bagenal, F., Kurth, W.S., Neubauer, F.M., Paranicas, C., Saur, J., 2004. In: Bagenal, F., Dowling, T.E., McKinnon, W.B. (Eds.), Jupiter: The Planet, Satellites and Magnetosphere. Cambridge University Press, Cambridge, UK, pp. 513-536]. The auroral footprint brightness is believed to be a good remote indicator of the strength of the interaction near Io, indicating the energy and current strength linking Io with Jupiter's ionosphere. The brightness may also depend in part on local auroral acceleration processes near Jupiter. The relative importance of different physical processes in this interaction can be tested as Jupiter's rotation and Io's orbital motion shift Jupiter's magnetic centrifugal equator past Io, leading to longitudinal variations in the plasma density near Io and functionally different variations in the local field strength near Jupiter where the auroral emissions are produced. Initial HST WFPC2 observations found a high degree of variability in the footprint brightness with time, and some evidence for systematic variations with longitude [Clarke, J.T., Ben Jaffel, L., Gérard, J.-C., 1998. J. Geophys. Res. 103, 20217-20236], however the data were not of sufficient quality to determine functional relationships. In this paper we report the results from a second, more thorough study, using a series of higher resolution and sensitivity HST STIS observations and a model for the center to limb dependence of the optically thin auroral emission brightness based on measurements of the auroral curtain emission distribution with altitude. A search for correlations between numerous parameters has revealed a strong dependence between Io's position in the plasma torus and the resulting footprint brightness that persists over several years of observations. The local magnetic field strength near Jupiter (i.e. the size of the loss cone) and the expected north/south asymmetry in auroral brightness related to the path of currents generated near Io through the plasma torus en route to Jupiter appear to be less important than the total plasma density near Io. This is consistent with the near-Io interaction being dominated by collisions of corotating plasma and mass pickup, a long-standing view which has been subject to considerable debate. The brightness of the auroral footprint emissions, however, does not appear to be proportional to the incident plasma density or energy, and the interpretation of this result will require detailed modeling of the interaction near Io.