Orbital evolution of small binary asteroids

About 15% of both near-Earth and main-belt asteroids with diameters below 10 km are now known to be binary. These small asteroid binaries are relatively uniform and typically contain a fast-spinning, flattened primary and a synchronously rotating, elongated secondary that is 20-40% as large (in diameter) as the primary. The principal formation mechanism for these binaries is now thought to be YORP (Yarkovsky-O’Keefe-Radzievskii-Paddack) effect induced spin-up of the primary followed by mass loss and accretion of the secondary from the released material. It has previously been suggested (?uk, M. [2007]. Astrophys. J. 659, L57-L60) that the present population of small binary asteroids is in a steady state between production through YORP and destruction through binary YORP (BYORP), which should increase or decrease secondary’s orbit, depending on the satellite’s shape. However, BYORP-driven evolution has not been directly modeled until now. Here we construct a simple numerical model of the binary’s orbital as well the secondary’s rotational dynamics which includes BYORP and selected terms representing main solar perturbations. We find that many secondaries should be vulnerable to chaotic rotation even for relatively low-eccentricity mutual orbits. We also find that the precession of the mutual orbit for typical small binary asteroids might be dominated by the perturbations from the prolate and librating secondary, rather than the oblate primary. When we evolve the mutual orbit by BYORP we find that the indirect effects on the binary’s eccentricity (through the coupling between the orbit and the secondary’s spin) dominate over direct ones caused by the BYORP acceleration. In particular, outward evolution causes eccentricity to increase and eventually triggers chaotic rotation of the secondary. We conclude that the most likely outcome will be reestablishing of the synchronous lock with a “flipped” secondary which would then evolve back in. For inward evolution we find an initial decrease of eccentricity and secondary’s librations, to be followed by later increase. We think that it is likely that various forms of dissipation we did not model may damp the secondary’s librations close to the primary, allowing for further inward evolution and a possible merger. We conclude that a merger or a tidal disruption of the secondary are the most likely outcomes of the BYORP evolution. Dissociation into heliocentric pairs by BYORP alone should be very difficult, and satellite loss might be restricted to the minority of systems containing more than one satellite at the time.     

Main belt binary asteroidal systems with eccentric mutual orbits

Using 8-10-m class telescopes and their Adaptive Optics (AO) systems, we conducted a long-term adaptive optics campaign initiated in 2003 focusing on four binary asteroid systems: (130) Elektra, (283) Emma, (379) Huenna, and (3749) Balam. The analysis of these data confirms the presence of their asteroidal satellite. We did not detect any additional satellite around these systems even though we have the capability of detecting a loosely-bound fragment (located at 1/4×RHill) ∼40 times smaller in diameter than the primary. The orbits derived for their satellites display significant eccentricity, ranging from 0.1 to 0.9, suggesting a different origin. Based on AO size estimate, we show that (130) Elektra and (283) Emma, G-type and P-type asteroids, respectively, have a significant porosity (30-60% considering CI-CO meteorites as analogs) and their satellite's eccentricities (e∼0.1) are possibly due to excitation by tidal effects. (379) Huenna and (3749) Balam, two loosely bound binary systems, are most likely formed by mutual capture. (3749) Balam's possible high bulk density is similar to (433) Eros, another S-type asteroid, and should be poorly fractured as well. (379) Huenna seems to display both characteristics: the moonlet orbits far away from the primary in term of stability (20%×RHill), but the primary's porosity is significant (30-60%).     

Particle size and abundance of HC3N ice in Titan’s lower stratosphere at high northern latitudes

Up to now, there has been no corroboration from Cassini CIRS of the Voyager IRIS-discovery of cyanoacetylene (HC₃N) ice in Titan’s thermal infrared spectrum. We report the first compelling spectral evidence from CIRS for the ν₆ HC₃N ice feature at 506 cm⁻¹ at latitudes 62°N and 70°N, from which we derive particle sizes and column abundances in Titan’s lower stratosphere. We find mean particle radii of 3.0 μm and 2.3 μm for condensed HC₃N at 62°N and 70°N, respectively, and corresponding ice phase molecular column abundances in the range 1-10 × 10¹⁶ mol cm⁻². Only upper limits for cloud abundances can be established at latitudes of 85°N, 55°N, 30°N, 10°N, and 15°S. Under the assumption that cloud tops coincide with the uppermost levels at which HC₃N vapor saturates, we infer geometric thicknesses for the clouds equivalent to 10-20 km or so, with tops at 165 km and 150 km at 70°N and 62°N, respectively.     

Saturn’s F ring grains: Aggregates made of crystalline water ice

We present models of the near-infrared (1-5 μm) spectra of Saturn’s F ring obtained by Cassini’s Visual and Infrared Mapping Spectrometer (VIMS) at ultra-high phase angles (177.4-178.5°). Modeling this spectrum constrains the size distribution, composition, and structure of F ring particles in the 0.1-100 μm size range. These spectra are very different from those obtained at lower phase angles; they lack the familiar 1.5 and 2 μm absorption bands, and the expected 3 μm water ice primary absorption appears as an unusually narrow dip at 2.87 μm. We have modeled these data using multiple approaches. First, we use a simple Mie scattering model to constrain the size distribution and composition of the particles. The Mie model allows us to understand the overall shapes of the spectra in terms of dominance by diffraction at these ultra-high phase angles, and also to demonstrate that the 2.87 μm dip is associated with the Christiansen frequency of water ice (where the real refractive index passes unity). Second, we use a combination of Mie scattering with Effective Medium Theory to probe the effect of porous (but structureless) particles on the overall shape of the spectrum and depth of the 2.87 μm band. Such simple models are not able to capture the shape of this absorption feature well. Finally, we model each particle as an aggregate of discrete monomers, using the Discrete Dipole Approximation (DDA) model, and find a better fit for the depth of the 2.87 μm feature. The DDA models imply a slightly different overall size distribution. We present a simple heuristic model which explains the differences between the Mie and DDA model results. We conclude that the F ring contains aggregate particles with a size distribution that is distinctly narrower than a typical power law, and that the particles are predominantly crystalline water ice.     

Radar and optical observations and physical modeling of triple near-Earth Asteroid (136617) 1994 CC

We report radar, photometric, and spectroscopic observations of near-Earth Asteroid (136617) 1994 CC. The radar measurements were obtained at Goldstone (8560 MHz, 3.5 cm) and Arecibo (2380 MHz, 12.6 cm) on 9 days following the asteroid’s approach within 0.0168 AU on June 10, 2009. 1994 CC was also observed with the Panchromatic Robotic Optical Monitoring and Polarimetry Telescopes (PROMPT) on May 21 and June 1-3. Visible-wavelength spectroscopy was obtained with the 5-m Hale telescope at Palomar on August 25. Delay-Doppler radar images reveal that 1994 CC is a triple system; along with (153591) 2001 SN263, this is only the second confirmed triple in the near-Earth population. Photometry obtained with PROMPT yields a rotation period for the primary P = 2.38860 ± 0.00009 h and a lightcurve amplitude of ∼0.1 mag suggesting a shape with low elongation. Hale telescope spectroscopy indicates that 1994 CC is an Sq-class object. Delay-Doppler radar images and shape modeling reveal that the primary has an effective diameter of 0.62 ± 0.06 km, low pole-on elongation, few obvious surface features, and a prominent equatorial ridge and sloped hemispheres that closely resemble those seen on the primary of binary near-Earth Asteroid (66391) 1999 KW4. Detailed orbit fitting reported separately by Fang et al. (Fang, J., Margot, J.-L., Brozovic, M., Nolan, M.C., Benner, L.A.M., Taylor, P.A. [2011]. Astron. J. 141, 154-168) gives a mass of the primary of 2.6 × 10¹¹ kg that, coupled with the effective diameter, yields a bulk density of 2.1 ± 0.6 g cm⁻³. The images constrain the diameters of the inner and outer satellites to be 113 ± 30 m and 80 ± 30 m, respectively. The inner satellite has a semimajor axis of ∼1.7 km (∼5.5 primary radii), an orbital period of ∼30 h, and its Doppler dispersion suggests relatively slow rotation, 26 ± 12 h, consistent with spin-orbit lock. The outer satellite has an orbital period of ∼9 days and a rotation period of 14 ± 7 h, establishing that the rotation is not spin-orbit locked. Among all binary and triple systems observed by radar, at least 25% (7/28) have a satellite that rotates more rapidly than its orbital period. This suggests that asynchronous configurations with P_rotation < P_orbital are relatively common among multiple systems in the near-Earth population. 1994 CC’s outer satellite has an observed maximum separation from the primary of ∼5.7 km (∼18.4 primary radii) that is the largest separation relative to primary radius seen to date among all 36 known binary and triple NEA systems. 1994 CC, (153591) 2001 SN263, and 1998 ST27 are the only triple and binary systems known with satellite separations >10 primary radii, suggesting either a detection bias, or that such widely-separated satellites are relatively uncommon in NEA multiple systems.     

Post-equinox observations of Uranus: Berg’s evolution, vertical structure, and track towards the equator

We present observations of Uranus taken with the near-infrared camera NIRC2 on the 10-m W.M. Keck II telescope, the Wide Field Planetary Camera 2 (WFPC2) and the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) from July 2007 through November 2009. In this paper we focus on a bright southern feature, referred to as the “Berg.” In Sromovsky et al. (Sromovsky, L.A., Fry, P.M., Hammel, H.B., Ahue, A.W., de Pater, I., Rages, K.A., Showalter, M.R., van Dam, M. [2009]. Icarus 203, 265-286), we reported that this feature, which oscillated between latitudes of −32° and −36° for several decades, suddenly started on a northward track in 2005. In this paper we show the complete record of observations of this feature’s track towards the equator, including its demise. After an initially slow linear drift, the feature’s drift rate accelerated at latitudes ∣θ∣ < 25°. By late 2009 the feature, very faint by then, was spotted at a latitude of −5° before disappearing from view. During its northward track, the feature’s morphology changed dramatically, and several small bright unresolved features were occasionally visible poleward of the main “streak.” These small features were sometimes visible at a wavelength of 2.2 μm, indicative that the clouds reached altitudes of ∼0.6 bar. The main part of the Berg, which is generally a long sometimes multipart streak, is estimated to be much deeper in the atmosphere, near 3.5 bars in 2004, but rising to 1.8-2.5 bars in 2007 after it began its northward drift. Through comparisons with Neptune’s Great Dark Spot and simulations of the latter, we discuss why the Berg may be tied to a vortex, an anticyclone deeper in the atmosphere that is visible only through orographic companion clouds.     

On Titan’s Xanadu region

A large, circular marking ∼1800 km across is seen in near-infrared images of Titan. The feature is centered at 10°S, 120°W on Titan, encompasses much of Titan’s western Xanadu region, and has an off-center, quasi-circular, inner margin about 700 km across, with lobate outer margins extending 200-500 km from the inner margin. On the feature’s southern flank is Tui Regio, an area that has very high reflectivity at 5 μm, and is hypothesized to exhibit geologically recent cryovolcanic flows (Barnes, J.W. et al. [2006]. Geophys. Res. Lett. 33), similar to flows seen in Hotei Regio, a cryovolcanic area whose morphology may be controlled by pre-existing, crustal fractures resulting from an ancient impact (Soderblom, L.A. et al. [2009]. Icarus, 204). The spectral reflectivity of the large, circular feature is quite different than that of its surroundings, making it compositionally distinct, and radar measurements of its topography, brightness temperature and volume scattering also suggest that the feature is quite distinct from its surroundings. These and several other lines of evidence, in addition to the feature’s morphology, suggest that it may occupy the site of an ancient impact.     

Initial sizes of planetesimals and accretion of the asteroids

The present size frequency distribution (SFD) of bodies in the asteroid belt appears to have preserved some record of the primordial population, with an excess of bodies of diameter D ∼ 100 km relative to a simple power law. The survival of Vesta’s basaltic crust also implies that the early SFD had a shallow slope in the range ∼10-100 km. (Morbidelli, A., Bottke, W.F., Nesvorny, D., Levison, H.F. [2009]. Icarus 204, 558-573) were unable to produce these features by accretion from an initial population of km-sized planetesimals. They concluded that bodies with sizes in the range ∼100-1000 km and a SFD similar to the current population were produced directly from solid particles of sub-meter scale, without experiencing accretion through intermediate sizes. We present results of new accretion simulations in the primordial asteroid region. The requisite SFD can be produced from an initial population of planetesimals of sizes ?0.1 km, smaller than the usual assumption of km-sized bodies. The bump at D ∼ 100 km is produced by a transition from dispersion-dominated runaway growth to a regime dominated by Keplerian shear, before the formation of large protoplanetary embryos. Thus, accretion of the asteroids from an initial population of small (sub-km) planetesimals cannot be ruled out.     

Titan impacts and escape

We report on hydrodynamic calculations of impacts of large (multi-kilometer) objects on Saturn’s moon Titan. We assess escape from Titan, and evaluate the hypothesis that escaping ejecta blackened the leading hemisphere of Iapetus and peppered the surface of Hyperion.We carried out two- and three-dimensional simulations of impactors ranging in size from 4 to 100 km diameter, impact velocities between 7 and 15 km s⁻¹, and impact angles from 0° to 75° from the vertical. We used the ZEUSMP2 hydrocode for the calculations. Simulations were made using three different geometries: three-dimensional Cartesian, two-dimensional axisymmetric spherical polar, and two-dimensional plane polar. Three-dimensional Cartesian geometry calculations were carried out over a limited domain (e.g. 240 km on a side for an impactor of size d_i = 10 km), and the results compared to ones with the same parameters done by Artemieva and Lunine (2005); in general the comparison was good. Being computationally less demanding, two-dimensional calculations were possible for much larger domains, covering global regions of the satellite (from 800 km below Titan’s surface to the exobase altitude 1700 km above the surface). Axisymmetric spherical polar calculations were carried out for vertical impacts. Two-dimensional plane-polar geometry calculations were made for both vertical and oblique impacts. In general, calculations among all three geometries gave consistent results.Our basic result is that the amount of escaping material is less than or approximately equal to the impactor mass even for the most favorable cases. Amounts of escaping material scaled most strongly as a function of velocity, with high-velocity impacts generating the largest amount, as expected. Dependence of the relative amount of escaping mass f_esc = m_esc/M_i on impactor diameter d_i was weak. Oblique impacts (impact angle θ_i > 45°) were more effective than vertical or near-vertical impacts; ratios of m_esc/M_i ∼ 1-2 were found in the simulations.     

Analysis of Titan CH4 3.3 μm upper atmospheric emission as measured by Cassini/VIMS

After molecular nitrogen, methane is the most abundant species in Titan’s atmosphere and plays a major role in its energy budget and its chemistry. Methane has strong bands at 3.3 μm emitting mainly at daytime after absorption of solar radiation. This emission is strongly affected by non-local thermodynamic equilibrium (non-LTE) in Titan’s upper atmosphere and, hence, an accurate modeling of the non-LTE populations of the emitting vibrational levels is necessary for its analysis. We present a sophisticated and extensive non-LTE model which considers 22 CH₄ levels and takes into account all known excitation mechanisms in which they take part. Solar absorption is the major excitation process controlling the population of the v₃-quanta levels above 1000 km whereas the distribution of the vibrational energy within levels of similar energy through collisions with N₂ is also of importance below that altitude. CH₄-CH₄ vibrational exchange of v₄-quanta affects their population below 500 km. We found that the ν₃ → ground band dominates Titan’s 3.3 μm daytime limb radiance above 750 km whereas the ν₃ + ν₄ → ν₄ band does below that altitude and down to 300 km. The ν₃ + ν₂ → ν₂, the 2ν₃ → ν₃, and the ¹³CH₄ν₃ → ground bands each contribute from 5% to 8% at regions below 800 km. The ν₃ + 2ν₄ → 2ν₄and ν₂ + ν₃ + ν₄ → ν₂ + ν₄ bands each contribute from 2% to 5% below 650 km. Contributions from other CH₄ bands are negligible. We have used the non-LTE model to retrieve the CH₄ abundance from 500 to 1100 km in the southern hemisphere from Cassini-VIMS daytime measurements near 3.3 μm. Our retrievals show good agreement with previous measurements and model results, supporting a weak deviation from well mixed values from the lower atmosphere up to 1000 km.     

Eclipse reappearances of Io: Time-resolved spectroscopy (1.9-4.2 μm)

We obtained time-resolved, near-infrared spectra of Io during the 60-90 min following its reappearance from eclipse by Jupiter on five occasions in 2004. The purpose was to search for spectral changes, particularly in the well-known SO₂ frost absorption bands, that would indicate surface-atmosphere exchange of gaseous SO₂ induced by temperature changes during eclipse. These observations were a follow-on to eclipse spectroscopy observations in which Bellucci et al. [Bellucci et al., 2004. Icarus 172, 141-148] reported significant changes in the strengths of two strong SO₂ bands in data acquired with the VIMS instrument aboard the Cassini spacecraft. One of the bands (4.07 μm [ν₁ + ν₃]) observed by Bellucci et al. is visible from ground-based observatories and is included in our data. We detected no changes in Io’s spectrum at any of the five observed events during the approximately 60-90 min during which spectra were obtained following Io’s emergence from Jupiter’s shadow. The areas of the three strongest SO₂ bands in the region 3.5-4.15 μm were measured for each spectrum; the variation of the band areas with time does not exceed that which can be explained by the Io’s few degrees of axial rotation during the intervals of observation, and in no case does the change in band strength approach that seen in the Cassini VIMS data. Our data are of sufficient quality and resolution to show the weak 2.198 μm (4549.6 cm⁻¹) 4ν₁ band of SO₂ frost on Io for what we believe is the first time. At one of the events (June 22, 2004), we began the acquisition of spectra ∼6 min before Io reappeared from Jupiter’s shadow, during which time it was detected through its own thermal emission. No SO₂ bands were superimposed on the purely thermal spectrum on this occasion, suggesting that the upper limit to condensed SO₂ in the vertical column above Io’s surface was ∼4 × 10⁻⁵ g cm⁻².     

Radar observations and a physical model of contact binary Asteroid 4486 Mithra

Arecibo (2380 MHz, 12.6 cm) and Goldstone (8560 MHz, 3.5 cm) delay-Doppler radar images obtained in July and August of 2000 reveal that 4486 Mithra is an irregular, significantly bifurcated object, with a central valley ∼380 m deep and a long axis potentially exceeding 2 km. With its bimodal appearance, Mithra is a strong candidate for a contact binary asteroid. Sequences of Goldstone images spanning up to 3 h per day show very little rotation and establish that Mithra is an unusually slow rotator. We used Goldstone and Arecibo data to estimate Mithra’s 3D shape and spin state. We obtain prograde (λ = 337°, β = 19°) and retrograde (λ = 154°, β = −19°) models that give comparable fits, have very similar shapes roughly resembling an hourglass, and have a rotation period of 67.5 ± 6.0 h. The dimensions of these two models are very similar; for the prograde solution the maximum dimensions are X = 2.35 ± 0.15 km, Y = 1.65 ± 0.10 km, Z = 1.44 ± 0.10 km. Dynamical analysis of our models suggests that in the past, Mithra most likely went through a period of even slower rotation with its obliquity close to 90°. The spin rate is predicted to be increasing due to thermal torque (YORP), while the obliquity, which is currently +68° and +106° for the prograde and retrograde models, respectively, is predicted to move away from 90°.     

The cool surfaces of binary near-Earth asteroids

Here we show results from thermal-infrared observations of km-sized binary near-Earth asteroids (NEAs). We combine previously published thermal properties for NEAs with newly derived values for three binary NEAs. The η value derived from the near-Earth asteroid thermal model (NEATM) for each object is then used to estimate an average thermal inertia for the population of binary NEAs and compared against similar estimates for the population of non-binaries. We find that these objects have, in general, surface temperatures cooler than the average values for non-binary NEAs as suggested by elevated η values. We discuss how this may be evidence of higher-than-average surface thermal inertia. This latter physical parameter is a sensitive indicator of the presence or absence of regolith: bodies covered with fine regolith, such as the Earth’s moon, have low thermal inertia, whereas a surface with little or no regolith displays high thermal inertia. Our results are suggestive of a binary formation mechanism capable of altering surface properties, possibly removing regolith: an obvious candidate is the YORP effect.We present also newly determined sizes and geometric visible albedos derived from thermal-infrared observations of three binary NEAs: (5381) Sekhmet, (153591) 2001 SN₂₆₃, and (164121) 2003 YT₁. The diameters of these asteroids are 1.41 ± 0.21 km, 1.56 ± 0.31 km, and 2.63 ± 0.40 km, respectively. Their albedos are 0.23 ± 0.13, 0.24 ± 0.16, and 0.048 ± 0.015, respectively.     

Coupled convection and tidal dissipation in Europa’s ice shell using non-Newtonian grain-size-sensitive (GSS) creep rheology

We present self-consistent, fully coupled two-dimensional (2D) numerical models of thermal evolution and tidal heating to investigate how convection interacts with tidal dissipation under the influence of non-Newtonian grain-size-sensitive creep rheology (plausibly resulting from grain boundary sliding) in Europa’s ice shell. To determine the thermal evolution, we solved the convection equations (using finite-element code ConMan) with the tidal dissipation as a heat source. For a given heterogeneous temperature field at a given time, we determined the tidal dissipation rate throughout the ice shell by solving for the tidal stresses and strains subject to Maxwell viscoelastic rheology (using finite-element code Tekton). In this way, the convection and tidal heating are fully coupled and evolve together. Our simulations show that the tidal dissipation rate can have a strong impact on the onset of thermal convection in Europa’s ice shell under non-Newtonian GSS rheology. By varying the ice grain size (1-10 mm), ice-shell thickness (20-120 km), and tidal-strain amplitude (0-4 × 10⁻⁵), we study the interrelationship of convection and conduction regimes in Europa’s ice shell. Under non-Newtonian grain-size-sensitive creep rheology and ice grain size larger than 1 mm, no thermal convection can initiate in Europa’s ice shell (for thicknesses <100 km) without tidal dissipation. However, thermal convection can start in thinner ice shells under the influence of tidal dissipation. The required tidal-strain amplitude for convection to occur decreases as the ice-shell thickness increases. For grain sizes of 1-10 mm, convection can occur in ice shells as thin as 20-40 km with the estimated tidal-strain amplitude of 2 × 10⁻⁵ on Europa.     

Waves in Cassini UVIS stellar occultations: 2. The C ring

We performed a complete wavelet analysis of Saturn’s C ring on 62 stellar occultation profiles. These profiles were obtained by Cassini’s Ultraviolet Imaging Spectrograph High Speed Photometer. We used a WWZ wavelet power transform to analyze them. With a co-adding process, we found evidence of 40 wavelike structures, 18 of which are reported here for the first time. Seventeen of these appear to be propagating waves (wavelength changing systematically with distance from Saturn). The longest new wavetrain in the C ring is a 52-km-long wave in a plateau at 86,397 km. We produced a complete map of resonances with external satellites and possible structures rotating with Saturn’s rotation period up to the eighth order, allowing us to associate a previously observed wave with the Atlas 2:1 inner Lindblad resonance (ILR) and newly detected waves with the Mimas 6:2 ILR and the Pandora 4:2 ILR. We derived surface mass densities and mass extinction coefficients, finding σ = 0.22(±0.03) g cm⁻² for the Atlas 2:1 ILR, σ = 1.31(±0.20) g cm⁻² for the Mimas 6:2 ILR, and σ = 1.42(±0.21) g cm⁻² for the Pandora 4:2 ILR. We determined a range of mass extinction coefficients (κ = τ/σ) for the waves associated with resonances with κ = 0.13 (±0.03) to 0.28(±0.06) cm² g⁻¹, where τ is the optical depth. These values are higher than the reported values for the A ring (0.01-0.02 cm² g⁻¹) and the Cassini Division (0.07-0.12 cm² g⁻¹ from Colwell et al. (Colwell, J.E., Cooney, J.H., Esposito, L.W., Srem?evi?, M. [2009]. Icarus 200, 574-580)). We also note that the mass extinction coefficient is probably not constant across the C ring (in contrast to the A ring and the Cassini Division): it is systematically higher in the plateaus than elsewhere, suggesting smaller particles in the plateaus. We present the results of our analysis of these waves in the C ring and estimate the mass of the C ring to be between3.7(±0.9) × 10¹⁶ kg and 7.9(±2.0) × 10¹⁶ kg (equivalent to an icy satellite of radius between 28.0(±2.3) km and 36.2(±3.0) km with a density of 400 kg m⁻³, close to that of Pan or Atlas). Using the ring viscosity derived from the wave damping length, we also estimate the vertical thickness of the C ring between 1.9(±0.4) m and 5.6(±1.4) m, comparable to the vertical thickness of the Cassini Division.     

First disk-resolved spectroscopy of (4) Vesta

Vesta, the second largest Main-Belt Asteroid, will be the first to be explored in 2011 by NASA’s Dawn mission. It is a dry, likely differentiated body with spectrum suggesting that is has been resurfaced by basaltic lava flows, not too different from the lunar maria.Here we present the first disk-resolved spectroscopic observations of an asteroid from the ground. We observed (4) Vesta with the ESO-VLT adaptive optics equipped integral-field near-infrared spectrograph SINFONI, as part of its science verification campaign. The highest spatial resolution of ∼90 km on Vesta’s surface was obtained during excellent seeing conditions (0.5^″) in October 2004.We observe spectral variations across Vesta’ surface that can be interpreted as variations of either the pyroxene composition, or the effect of surface aging. We compare Vesta’s 2 μm absorption band to that of howardite-eucrite-diogenite (HED) meteorites that are thought to originate from Vesta, and establish particular links between specific regions and HED subclasses. The overall composition is found to be mostly compatible with howardite meteorites, although a small area around 180°E longitude could be attributed to a diogenite-rich spot. We finally focus our spectral analysis on the characteristics of Vesta’s bright and dark regions as seen from Hubble Space Telescope’s visible and Keck-II’s near-infrared images.     

Latitudinal variations in Titan’s methane and haze from Cassini VIMS observations

We analyze observations taken with Cassini’s Visual and Infrared Mapping Spectrometer (VIMS), to determine the current methane and haze latitudinal distribution between 60°S and 40°N. The methane variation was measured primarily from its absorption band at 0.61 μm, which is optically thin enough to be sensitive to the methane abundance at 20-50 km altitude. Haze characteristics were determined from Titan’s 0.4-1.6 μm spectra, which sample Titan’s atmosphere from the surface to 200 km altitude. Radiative transfer models based on the haze properties and methane absorption profiles at the Huygens site reproduced the observed VIMS spectra and allowed us to retrieve latitude variations in the methane abundance and haze. We find the haze variations can be reproduced by varying only the density and single scattering albedo above 80 km altitude. There is an ambiguity between methane abundance and haze optical depth, because higher haze optical depth causes shallower methane bands; thus a family of solutions is allowed by the data. We find that haze variations alone, with a constant methane abundance, can reproduce the spatial variation in the methane bands if the haze density increases by 60% between 20°S and 10°S (roughly the sub-solar latitude) and single scattering absorption increases by 20% between 60°S and 40°N. On the other hand, a higher abundance of methane between 20 and 50 km in the summer hemisphere, as much as two times that of the winter hemisphere, is also possible, if the haze variations are minimized. The range of possible methane variations between 27°S and 19°N is consistent with condensation as a result of temperature variations of 0-1.5 K at 20-30 km. Our analysis indicates that the latitudinal variations in Titan’s visible to near-IR albedo, the north/south asymmetry (NSA), result primarily from variations in the thickness of the darker haze layer, detected by Huygens DISR, above 80 km altitude. If we assume little to no latitudinal methane variations we can reproduce the NSA wavelength signatures with the derived haze characteristics. We calculate the solar heating rate as a function of latitude and derive variations of ∼10-15% near the sub-solar latitude resulting from the NSA. Most of the latitudinal variations in the heating rate stem from changes in solar zenith angle rather than compositional variations.     

Titan’s atomic hydrogen corona

Based on measurements performed by the Hydrogen Deuterium Absorption Cell (HDAC) aboard the Cassini orbiter, Titan’s atomic hydrogen exosphere is investigated. Data obtained during the T9 encounter are used to infer the distribution of atomic hydrogen throughout Titan’s exosphere, as well as the exospheric temperature.The measurements performed during the flyby are modeled by performing Monte Carlo radiative transfer calculations of solar Lyman-α radiation, which is resonantly scattered on atomic hydrogen in Titan’s exosphere. Two different atomic hydrogen distribution models are applied to determine the best fitting density profile. One model is a static model that uses the Chamberlain formalism to calculate the distribution of atomic hydrogen throughout the exosphere, whereas the second model is a Particle model, which can also be applied to non-Maxwellian velocity distributions.The density distributions provided by both models are able to fit the measurements although both models differ at the exobase: best fitting exobase atomic hydrogen densities of n_H = (1.5 ± 0.5) × 10⁴ cm⁻³ and n_H = (7 ± 1) × 10⁴ cm⁻³ were found using the density distribution provided by both models, respectively. This is based on the fact that during the encounter, HDAC was sensitive to altitudes above about 3000 km, hence well above the exobase at about 1500 km. Above 3000 km, both models produce densities which are comparable, when taking into account the measurement uncertainty.The inferred exobase density using the Chamberlain profile is a factor of about 2.6 lower than the density obtained from Voyager 1 measurements and much lower than the values inferred from current photochemical models. However, when taking into account the higher solar activity during the Voyager flyby, this is consistent with the Voyager measurements. When using the density profile provided by the particle model, the best fitting exobase density is in perfect agreement with the densities inferred by current photochemical models.Furthermore, a best fitting exospheric temperature of atomic hydrogen in the range of T_H = (150-175) ± 25 K was obtained when assuming an isothermal exosphere for the calculations. The required exospheric temperature depends on the density distribution chosen. This result is within the temperature range determined by different instruments aboard Cassini. The inferred temperature is close to the critical temperature for atomic hydrogen, above which it can escape hydrodynamically after it diffused through the heavier background gas.     

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

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

A geomorphic analysis of Hale crater, Mars: The effects of impact into ice-rich crust

Hale crater, a 125 × 150 km impact crater located near the intersection of Uzboi Vallis and the northern rim of Argyre basin at 35.7°S, 323.6°E, is surrounded by channels that radiate from, incise, and transport material within Hale’s ejecta. The spatial and temporal relationship between the channels and Hale’s ejecta strongly suggests the impact event created or modified the channels and emplaced fluidized debris flow lobes over an extensive area (>200,000 km²). We estimate ∼10¹⁰ m³ of liquid water was required to form some of Hale’s smaller channels, a volume we propose was supplied by subsurface ice melted and mobilized by the Hale-forming impact. If 10% of the subsurface volume was ice, based on a conservative porosity estimate for the upper martian crust, 10¹² m³ of liquid water could have been present in the ejecta. We determine a crater-retention age of 1 Ga inside the primary cavity, providing a minimum age for Hale and a time at which we propose the subsurface was volatile-rich. Hale crater demonstrates the important role impacts may play in supplying liquid water to the martian surface: they are capable of producing fluvially-modified terrains that may be analogous to some landforms of Noachian Mars.