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

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

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.     

Formation of Jupiter using opacities based on detailed grain physics

Numerical simulations, based on the core-nucleated accretion model, are presented for the formation of Jupiter at 5.2 AU in three primordial disks with three different assumed values of the surface density of solid particles. The grain opacities in the envelope of the protoplanet are computed using a detailed model that includes settling and coagulation of grains and that incorporates a recalculation of the grain size distribution at each point in time and space. We generally find lower opacities than the 2% of interstellar values used in previous calculations (Hubickyj, O., Bodenheimer, P., Lissauer, J.J. [2005]. Icarus 179, 415-431; Lissauer, J.J., Hubickyj, O., D’Angelo, G., Bodenheimer, P. [2009]. Icarus 199, 338-350). These lower opacities result in more rapid heat loss from and more rapid contraction of the protoplanetary envelope. For a given surface density of solids, the new calculations result in a substantial speedup in formation time as compared with those previous calculations. Formation times are calculated to be 1.0, 1.9, and 4.0 Myr, and solid core masses are found to be 16.8, 8.9, and 4.7 M_⊕, for solid surface densities, σ, of 10, 6, and 4 g cm⁻², respectively. For σ = 10 and σ = 6 g cm⁻², respectively, these formation times are reduced by more than 50% and more than 80% compared with those in a previously published calculation with the old approximation to the opacity.     

A radar survey of M- and X-class asteroids II. Summary and synthesis

Using the S-band radar at Arecibo Observatory, we observed six new M-class main-belt asteroids (MBAs), and re-observed one, bringing the total number of Tholen M-class asteroids observed with radar to 19. The mean radar albedo for all our targets is , significantly higher than the mean radar albedo of every other class (Magri, C., Nolan, M.C., Ostro, S.J., Giorgini, J.D. [2007]. Icarus 186, 126-151). Seven of these objects (Asteroids 16 Psyche, 129 Antigone, 216 Kleopatra, 347 Pariana, 758 Mancunia, 779 Nina, 785 Zwetana) have radar albedos indicative of a very high metal content , and consistent with a remnant iron/nickel core interpretation (irons) or exotic high metal meteorite types such as CB. We propose designating these high radar albedo objects as Mm. Two asteroids, 110 Lydia and 678 Fredegundis, have more moderate radar albedos , but exhibit high values at some rotation phases suggesting a significant metal content. The remaining 10 objects have moderate radar albedos at all rotation phases. Most of our targets have visible/near-infrared spectra (Hardersen, P.S., Gaffey, M.J., Abell, P.A. [2005]. Icarus 175, 141-158; Fornasier, S., Clark, B.E., Dotto, E., Migliorini, A., Ockert-Bell, M., Barucci, M.A. [2009]. Icarus, submitted for publication) that indicate the presence of at least some silicate phases. All of the non-Mm asteroids show a positive correlation between visual and radar albedo but the reasons for this are not clear. All of the higher radar albedo targets (the 7 Mm asteroids, Lydia, and Fredegundis) show moderate to large variations in radar albedo with rotation phase. We suggest that their high radar reflectivity exaggerates irregularities in the asteroid shape to cause this behavior. One-third of our targets show evidence for asteroid-scale concavities or bifurcation. Based on all the evidence available, we suggest that most Tholen M-class asteroids are not remnant iron cores or enstatite chondrites, but rather collisional composites of silicates and irons with compositions more analogous to stony-iron meteorites and high-iron carbonaceous chondrites.     

Ice lines, planetesimal composition and solid surface density in the solar nebula

To date, there is no core accretion simulation that can successfully account for the formation of Uranus or Neptune within the observed 2–3 Myr lifetimes of protoplanetary disks. Since solid accretion rate is directly proportional to the available planetesimal surface density, one way to speed up planet formation is to take a full accounting of all the planetesimal-forming solids present in the solar nebula. By combining a viscously evolving protostellar disk with a kinetic model of ice formation, which includes not just water but methane, ammonia, CO and 54 minor ices, we calculate the solid surface density of a possible giant planet-forming solar nebula as a function of heliocentric distance and time. Our results can be used to provide the starting planetesimal surface density and evolving solar nebula conditions for core accretion simulations, or to predict the composition of planetesimals as a function of radius. We find three effects that favor giant planet formation by the core accretion mechanism: (1) a decretion flow that brings mass from the inner solar nebula to the giant planet-forming region, (2) the fact that the ammonia and water ice lines should coincide, according to recent lab results from Collings et al. [Collings, M.P., Anderson, M.A., Chen, R., Dever, J.W., Viti, S., Williams, D.A., McCoustra, M.R.S., 2004. Mon. Not. R. Astron. Soc. 354, 1133–1140], and (3) the presence of a substantial amount of methane ice in the trans-saturnian region. Our results show higher solid surface densities than assumed in the core accretion models of Pollack et al. [Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y., 1996. Icarus 124, 62–85] by a factor of 3–4 throughout the trans-saturnian region. We also discuss the location of ice lines and their movement through the solar nebula, and provide new constraints on the possible initial disk configurations from gravitational stability arguments.     

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.     

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.     

CCD photometry of 23 minor planets

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

Low-speed impacts between rubble piles modeled as collections of polyhedra

We present results of modeling rubble piles as collections of polyhedra. The use of polyhedra allows more realistic (irregular) shapes and interactions (e.g. collisions), particularly for objects of different sizes. Rotational degrees of freedom are included in the modeling, which may be important components of the motion. We solved the equations of rigid-body dynamics, including frictional/inelastic collisions, for collections of up to several hundred elements. As a demonstration of the methods and to compare with previous work by other researchers, we simulated low-speed collisions between km-scale bodies with the same general parameters as those simulated by Leinhardt et al. [Leinhardt, Z.M., Richardson, D.C., Quinn, T., 2000. Icarus 146, 133-151]. High-speed collisions appropriate to present-day asteroid encounters require additional treatment of shock effects and fragmentation and are the subject of future work; here we study regimes appropriate to planetesimal accretion and re-accretion in the aftermath of catastrophic events. Collisions between equal-mass objects at low speeds () were simulated for both head-on and off-center collisions between rubble piles made of a power-law mass spectrum of sub-elements. Very low-speed head-on collisions produce single objects from the coalescence of the impactors. For slightly higher speeds, extensive disruption occurs, but re-accretion produces a single object with most of the total mass. For increasingly higher speeds, the re-accreted object has smaller mass, finally resulting in complete catastrophic disruption with all sub-elements on escape trajectories and only small amounts of mass in re-accreted bodies. Off-center collisions at moderately low speeds produce two re-accreted objects of approximately equal mass, separating at greater than escape speed. At high speed, complete disruption occurs as with the high-speed head-on collisions. Head-on collisions at low to moderate speeds result in objects of mostly oblate shape, while higher speed collisions produce mostly prolate objects, as do off-center collisions at moderate and high speeds. Collisions carried out with the same dissipative coefficients (coefficient of restitution ?n=0.8, zero friction) as used by Leinhardt et al. [Leinhardt, Z.M., Richardson, D.C., Quinn, T., 2000. Icarus 146, 133-151] result in a value for specific energy for disruption , somewhat lower than the value of 2 J/kg found by them, while collisions with a lower coefficient of restitution and friction [?n=0.5, ?t=0, μ=0.5, similar to those used by Michel, et al. [Michel, P., Benz, W., Richardson, D.C., 2004. Planet. Space Sci. 52, 1109-1117] for SPH + N-body calculations] yield .     

The size and thermal properties of the nucleus of Comet 22P/Kopff

We detected the nucleus of Comet 22P/Kopff at 4.87 AU from the Sun with the two IRS peak-up cameras of the Spitzer Space Telescope on April 19, 2007. Using the thermal model of [Groussin, O., and 15 colleagues, 2007. Icarus 187, 16-25], we derive a nucleus size of 1.89±0.16 km, in agreement with [Lamy, P., Toth, I., Jorda, L., Groussin, O., A'Hearn, M.F., Weaver, H.A., 2002. Icarus 156, 442-455], and a thermal inertia .     

Thermal infrared and optical observations of four near-Earth asteroids

We present thermal infrared photometry and spectrophotometry of four near-Earth asteroids (NEAs), namely (433) Eros, (66063) 1998 RO₁, (137032) 1998 UO₁, and (138258) 2000 GD₂, using the United Kingdom Infrared Telescope (UKIRT) in 2002. For two objects, i.e. (433) Eros and (137032) 1998 UO₁, quasi-simultaneous optical observations were also obtained, using the Jacobus Kapteyn Telescope (JKT). For (127032) 1998 UO₁, we obtain a rotation period P=3.0±0.1 h and an absolute visual magnitude HV=16.7±0.4. The Standard Thermal Model (STM), Fast Rotating Model (FRM) and near-Earth asteroid Thermal Model (NEATM) have been fitted to the IR fluxes to determine effective diameters Deff, geometric albedos pv, and beaming parameters η. The derived values are (433) Eros: Deff=23.3±3.5 km (at lightcurve maximum), pv=0.24±0.07, η=0.95±0.19; (66063) 1998 RO₁: , ; (137032) 1998 UO₁: Deff<1.13 km, pv>0.29; (138258) 2000 GD₂: Deff=0.27±0.04 km, , η=0.74±0.15. (66063) 1998 RO₁ is a binary asteroid from lightcurve characteristics [Pravec, P., and 56 colleagues, 2006. Icarus 181, 63-93] and we estimate the effective diameter of the primary (Dp) and secondary (Ds) components: and . The diameter and albedo of (138258) 2000 GD₂ are consistent with the trend of decreasing diameter for S- and Q-type asteroids found by Delbó et al. [Delbó, M., Harris, A.W., Binzel, R.P., Pravec, P., Davies, J.K., 2003. Icarus 166, 116-130]. A possible trend of increasing beaming parameter with diameter for small (less than about 3 km) S- and Q-type asteroids is found.     

Estimating the masses of Saturn’s A and B rings from high-optical depth N-body simulations and stellar occultations

We have completed a series of local N-body simulations of Saturn’s B and A rings in order to identify systematic differences in the degree of particle clumping into self-gravity wakes as a function of orbital distance from Saturn and dynamical optical depth (a function of surface density). These simulations revealed that the normal optical depth of the final configuration can be substantially lower than one would infer from a uniform distribution of particles. Adding more particles to the simulation simply piles more particles onto the self-gravity wakes while leaving relatively clear gaps between the wakes. Estimating the mass from the observed optical depth is therefore a non-linear problem. These simulations may explain why the Cassini UVIS instrument has detected starlight at low incidence angles through regions of the B ring that have average normal optical depths substantially greater than unity at some observation geometries [Colwell, J.E., Esposito, L.W., Srem?evi?, M., Stewart, G.R., McClintock, W.E., 2007. Icarus 190, 127-144]. We provide a plausible internal density of the particles in the A and B rings based upon fitting the results of our simulations with Cassini UVIS stellar occultation data. We simulated Cassini-like occultations through our simulation cells, calculated optical depths, and attempted to extrapolate to the values that Cassini observes. We needed to extrapolate because even initial optical depths of >4 (σ > 240 g cm⁻²) only yielded final optical depths no greater than 2.8, smaller than the largest measured B ring optical depths. This extrapolation introduces a significant amount of uncertainty, and we chose to be conservative in our overall mass estimates. From our simulations, we infer the surface density of the A ring to be , which corresponds to a mass of . We infer a minimum surface density of for Saturn’s B ring, which corresponds to a minimum mass estimate of . The A ring mass estimate agrees well with previous analyses, while the B ring is at least 40% larger. In sum, our lower limit estimate is that the total mass of Saturn’s ring system is 120-200% the mass of the moon Mimas, but significantly larger values would be plausible given the limitations of our simulations. A significantly larger mass for Saturn’s rings favors a primordial origin for the rings because the disruption of a former satellite of the required mass would be unlikely after the decay of the late heavy bombardment of planetary surfaces.     

Gemini-N mid-IR observations of the dust properties of the ejecta excavated from Comet 9P/Tempel 1 during Deep Impact

We present mid-infrared spectra and images from the Gemini-N (+MICHELLE) observational campaign of Comet 9P/Tempel 1 before, during, and after its encounter with Deep Impact. We use our thermal grain model to probe the 10 μm properties of the dust grains in the coma of the comet. Before impact (3 July 2005 UT), and more than 24 h after impact (5, 16, and 28 July 2005 UT), the comet dust grains were composed mostly of amorphous olivine, and were relatively large (peak of the grain size distribution ). For the night of impact, we extract spectra by centering on the nucleus, and offset 1″ from the nucleus in the direction of the impact ejecta plume. We find small dust grains (∼0.2 μm) of a diverse mineralogy (amorphous olivine, amorphous pyroxene, amorphous carbon, and crystalline olivine) populating the ejecta. The submicron sized dust grains move faster than the other, larger grains (?0.7 μm), with amorphous olivine and amorphous carbon traveling together, and amorphous pyroxene and crystalline olivine dispersing at a similar rate. Deriving a velocity law from a time-of-flight analysis, we find that the material traveled with a velocity law scaled by and with a power of p=0.5. This velocity power-law requires a sustained release of grains for the duration of 45-60 min after impact. Since the mineral species are traveling at different speeds, and there was a sustained release of grains due to a possible “gas-plume,” we conclude that the different minerals did not originate from grain aggregates destroyed by the impact, but instead arise from an inhomogeneous nucleus.     

Titan's hydrodynamically escaping atmosphere: Escape rates and the structure of the exobase region

In Strobel [Strobel, D.F., 2008. Icarus, 193, 588-594] a mass loss rate from Titan's upper atmosphere, , was calculated for a single constituent, N₂ atmosphere by hydrodynamic escape as a high density, slow outward expansion driven principally by solar UV heating due to CH₄ absorption. It was estimated, but not proven, that the hydrodynamic mass loss is essentially CH₄ and H₂ escape. Here the individual conservation of momentum equations for the three major components of the upper atmosphere (N₂, CH₄, H₂) are solved in the low Mach number limit and compared with Cassini Ion Neutral Mass Spectrometer (INMS) measurements to demonstrate that light gases (CH₄, H₂) preferentially escape over the heavy gas (N₂). The lightest gas (H₂) escapes with a flux 99% of its limiting flux, whereas CH₄ is restricted to ?75% of its limiting flux because there is insufficient solar power to support escape at the limiting rate. The respective calculated H₂ and CH₄ escape rates are 9.2×¹⁰²⁷ and 1.7×¹⁰²⁷ s−1, for a total of . From the calculated densities, mean free paths of N₂, CH₄, H₂, and macroscopic length scales, an extended region above the classic exobase is inferred where frequent collisions are still occurring and thermal heat conduction can deliver power to lift the escaping gas out of the gravitational potential well. In this region rapid acceleration of CH₄ outflow occurs. With the thermal structure of Titan's thermosphere inferred from INMS data by Müller-Wodarg et al. [Müller-Wodarg, I.C.F., Yelle, R.V., Cui, J., Waite Jr., J.H., 2008. J. Geophys. Res. 113, doi:10.1029/2007JE003033. E10005], in combination with calculated temperature profiles that include sputter induced plasma heating at the exobase, it is concluded that on average that the integrated, globally average, orbit-averaged, plasma heating rate during the Cassini epoch does not exceed ().     

Mass loss due to sputtering and thermal processes in meteoroid ablation

Conventional meteoroid theory assumes that the dominant mode of ablation (which we will refer to as thermal ablation) is by evaporation following intense heating during atmospheric flight. Light production results from excitation of ablated meteoroid atoms following collisions with atmospheric constituents. In this paper, we consider the question of whether sputtering may provide an alternative disintegration process of some importance. For meteoroids in the mass range from 10^-3 to and covering a meteor velocity range from 11 to , we numerically modeled both thermal ablation and sputtering ablation during atmospheric flight. We considered three meteoroid models believed to be representative of asteroidal ( mass density), cometary () and porous cometary () meteoroid structures. Atmospheric profiles which considered the molecular compositions at different heights were use in the sputtering calculations. We find that while in many cases (particularly at low velocities and for relatively large meteoroid masses) sputtering contributes only a small amount of mass loss during atmospheric flight, in some cases sputtering is very important. For example, a porous meteoroid at will lose nearly 51% of its mass by sputtering, while a asteroidal meteoroid at will lose nearly 83% of its mass by sputtering. We argue that sputtering may explain the light production observed at very great heights in some Leonid meteors. We discuss methods to observationally test the predictions of these computations. A search for early gradual tails on meteor light curves prior to the commencement of intense thermal ablation possibly represents the most promising approach. The impact of this work will be most dramatic for very small meteoroids such as those observed with large aperture radars. The heights of ablation and decelerations observed using these systems may provide evidence for the importance of sputtering.     

Multi-wavelength observations of Asteroid 2100 Ra-Shalom

We observed near-Earth asteroid (NEA) 2100 Ra-Shalom over a six-year period, obtaining rotationally resolved spectra in the visible, near-infrared, thermal-infrared, and radar wavelengths. We find that Ra-Shalom has an effective diameter of Deff=2.3±0.2 km, rotation period P=19.793±0.001 h, visual albedo pv=0.13±0.03, radar albedo , and polarization ratio μc=0.25±0.04. We used our radar observations to generate a three-dimensional shape model which shows several structural features of interest. Based on our thermal observations, Ra-Shalom has a high thermal inertia of ∼10³ J m⁻² s^−0.5 K⁻¹, consistent with a coarse or rocky surface and the inferences of others [Harris, A.W., Davies, J.K., Green, S.F., 1998. Icarus 135, 441-450; Delbo, M., Harris, A.W., Binzel, R.P., Pravec, P., Davies, J.K., 2003. Icarus 166, 116-130]. Our spectral data indicate that Ra-Shalom is a K-class asteroid and we find excellent agreement between our spectra and laboratory spectra of the CV3 meteorite Grosnaja. Our spectra show rotation-dependent variations consistent with global variations in grain size. Our radar observations show rotation-dependent variations in radar albedo consistent with global variations in the thickness of a relatively thin regolith.     

Diffuse reflectance of altered olivine grains: Remote sensing detection and implications for Mars studies

The study of Mars linked to its climatic history, the presence of water and possible forms of life on the planet, are becoming more and more important in recent years. In this paper we present laboratory reflection measurements on olivine samples which are believed to be the product of alteration induced by their coexistence with water. The reflection spectra have been obtained in the wavelength range from 0.2 to and for different grain sizes; they exhibit, among various bands, some of which due to H₂O and OH^-, a characteristic absorption feature at about , whose emission wings at 0.56 and near have relative strengths which are size dependent and well correlated. The identification of such feature in the observed spectra of Mars would provide useful information about the grain size of the Martian regolith and, in association with other bands, also about the possible presence of olivine altered by water in the past.