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.     

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

We present the results of additional calculations involving the collisions of km-scale rubble piles. In new work, we used the Open Dynamics Engine (ODE), an open-source library for the simulation of rigid-body dynamics that incorporates a sophisticated collision-detection and resolution routine. We found that using ODE resulted in a speed-up of approximately a factor of 30 compared with previous code. In this paper we report on the results of almost 1200 separate runs, the bulk of which were carried out with 1000-2000 elements. We carried out calculations with three different combinations of the coefficients of friction η and (normal) restitution ?: low (η=0,?=0.8), medium (η=0,?=0.5), and high (η=0.5,?=0.5) dissipation.For target objects of ∼1 km in radius, we found reduced critical disruption energy values in head-on collisions from 2 to 100 J kg⁻¹ depending on dissipation and impactor/target mass ratio. Monodisperse objects disrupted somewhat more easily than power-law objects in general. For oblique collisions of equal-mass objects, mildly off-center collisions (b/b₀=0.5) seemed to be as efficient or possibly more efficient at collisional disruption as head-on collisions. More oblique collisions were less efficient and the most oblique collisions we tried (b/b₀=0.866) required up to ∼200 J kg⁻¹ for high-dissipation power-law objects. For calculations with smaller numbers of elements (total impactor or 200 elements) we found that collisions were more efficient for smaller numbers of more massive elements, with values as low as for low-dissipation cases. We also analyzed our results in terms of the relations proposed by Stewart and Leinhardt [Stewart, S.T., Leinhardt, Z.M., 2009. Astrophys. J. 691, L133-L137] where where Q_R is the impact kinetic energy per unit total mass m_i+m_T. Although there is a significant amount of scatter, our results generally bear out the suggested relation.     

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.     

Dynamical erosion of the asteroid belt and implications for large impacts in the inner Solar System

The cumulative effects of weak resonant and secular perturbations by the major planets produce chaotic behavior of asteroids on long timescales. Dynamical chaos is the dominant loss mechanism for asteroids with diameters in the current asteroid belt. In a numerical analysis of the long-term evolution of test particles in the main asteroid belt region, we find that the dynamical loss history of test particles from this region is well described with a logarithmic decay law. In our simulations the loss rate function that is established at persists with little deviation to at least . Our study indicates that the asteroid belt region has experienced a significant amount of depletion due to this dynamical erosion—having lost as much as ∼50% of the large asteroids—since 1 Myr after the establishment of the current dynamical structure of the asteroid belt. Because the dynamical depletion of asteroids from the main belt is approximately logarithmic, an equal amount of depletion occurred in the time interval 10-200 Myr as in 0.2-4 Gyr, roughly ∼30% of the current number of large asteroids in the main belt over each interval. We find that asteroids escaping from the main belt due to dynamical chaos have an Earth-impact probability of ∼0.3%. Our model suggests that the rate of impacts from large asteroids has declined by a factor of 3 over the last 3 Gyr, and that the present-day impact flux of objects on the terrestrial planets is roughly an order of magnitude less than estimates currently in use in crater chronologies and impact hazard risk assessments.     

Eclipsing binary Trojan asteroid Patroclus: Thermal inertia from Spitzer observations

We present mid-infrared observations of the binary L5-Trojan system (617) Patroclus-Menoetius before, during, and after two shadowing events, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. For the first time, we effectively observe changes in asteroid surface temperature in real time, allowing the thermal inertia to be determined very directly. A new detailed binary thermophysical model is presented which accounts for the system’s known mutual orbit, arbitrary component shapes, and thermal conduction in the presence of eclipses.We obtain two local thermal-inertia values, representative of the respective shadowed areas: and . The average thermal inertia is estimated to be , potentially with significant surface heterogeneity. This first thermal-inertia measurement for a Trojan asteroid indicates a surface covered in fine regolith. Independently, we establish the presence of fine-grained (<a few μm) silicates on the surface, based on emissivity features near 10 and similar to those previously found on other Trojans.We also report V-band observations and report a lightcurve with complete rotational coverage. The lightcurve has a low amplitude of peak-to-peak, implying a roughly spherical shape for both components, and is single-periodic with a period equal to the period of the mutual orbit, indicating that the system is fully synchronized.The diameters of Patroclus and Menoetius are 106±11 and , respectively, in agreement with previous findings. Taken together with the system’s known total mass, this implies a bulk mass density of , significantly below the mass density of L4-Trojan asteroid (624) Hektor and suggesting a bulk composition dominated by water ice.All known physical properties of Patroclus, arguably the best studied Trojan asteroid, are consistent with those expected in icy objects with devolatilized surface (extinct comets), consistent with what might be implied by recent dynamical modeling in the framework of the Nice Model.     

A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris

With the collection of six years of MGS tracking data and three years of Mars Odyssey tracking data, there has been a continual improvement in the JPL Mars gravity field determination. This includes the measurement of the seasonal changes in the gravity coefficients (e.g., , , , , , ) caused by the mass exchange between the polar ice caps and atmosphere. This paper describes the latest gravity field MGS95J to degree and order 95. The improvement comes from additional tracking data and the adoption of a more complete Mars orientation model with nutation, instead of the IAU 2000 model. Free wobble of the Mars' spin axis, i.e. polar motion, has been constrained to be less than 10 mas by looking at the temporal history of and . A strong annual signature is observed in , and this is a mixture of polar motion and ice mass redistribution. The Love number solution with a subset of Odyssey tracking data is consistent with the previous liquid outer core determination from MGS tracking data [Yoder et al., 2003. Science 300, 299-303], giving a combined solution of k₂=0.152±0.009 using MGS and Odyssey tracking data. The solutions for the masses of the Mars' moons show consistency between MGS, Odyssey, and Viking data sets; Phobos GM=(7.16±0.005)×¹⁰⁻⁴ km³/s² and Deimos GM=(0.98±0.07)×¹⁰⁻⁴ km³/s². Average MGS orbit errors, determined from differences in the overlaps of orbit solutions, have been reduced to 10-cm in the radial direction and 1.5 m along the spacecraft velocity and normal to the orbit plane. Hence, the ranging to the MGS and Odyssey spacecraft has resulted in position measurements of the Mars system center-of-mass relative to the Earth to an accuracy of one meter, greatly reducing the Mars ephemeris errors by several orders of magnitude, and providing mass estimates for Asteroids 1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, and 324 Bamberga.     

The dust trail of Comet 67P/Churyumov-Gerasimenko between 2004 and 2006

We report on observations of the dust trail of Comet 67P/Churyumov-Gerasimenko (CG) in visible light with the Wide Field Imager at the ESO/MPG 2.2 m telescope at 4.7 AU before aphelion, and at with the MIPS instrument on board the Spitzer Space Telescope at 5.7 AU both before and after aphelion. The comet did not appear to be active during our observations. Our images probe large dust grains emitted from the comet that have a radiation pressure parameter β<0.01. We compare our observations with simulated images generated with a dynamical model of the cometary dust environment and constrain the emission speeds, size distribution, production rate and geometric albedo of the dust. We achieve the best fit to our data with a differential size distribution exponent of −4.1, and emission speeds for a β=0.01 particle of 25 m/s at perihelion and 2 m/s at 3 AU. The dust production rate in our model is on the order of 1000 kg/s at perihelion and 1 kg/s at 3 AU, and we require a dust geometric albedo between 0.022 and 0.044. The production rates of large (>) particles required to reproduce the brightness of the trail are sufficient to also account for the coma brightness observed while the comet was inside 3 AU, and we infer that the cross-section in the coma of CG may be dominated by grains of the order of .     

Fragment properties at the catastrophic disruption threshold: The effect of the parent body’s internal structure

Numerical simulations of asteroid breakups, including both the fragmentation of the parent body and the gravitational interactions between the fragments, have allowed us to reproduce successfully the main properties of asteroid families formed in different regimes of impact energy, starting from a non-porous parent body. In this paper, using the same approach, we concentrate on a single regime of impact energy, the so-called catastrophic threshold usually designated by , which results in the escape of half of the target’s mass. Thanks to our recent implementation of a model of fragmentation of porous materials, we can characterize for both porous and non-porous targets with a wide range of diameters. We can then analyze the potential influence of porosity on the value of , and by computing the gravitational phase of the collision in the gravity regime, we can characterize the collisional outcome in terms of the fragment size and ejection speed distributions, which are the main outcome properties used by collisional models to study the evolutions of the different populations of small bodies. We also check the dependency of on the impact speed of the projectile.In the strength regime, which corresponds to target sizes below a few hundreds of meters, we find that porous targets are more difficult to disrupt than non-porous ones. In the gravity regime, the outcome is controlled purely by gravity and porosity in the case of porous targets. In the case of non-porous targets, the outcome also depends on strength. Indeed, decreasing the strength of non-porous targets make them easier to disrupt in this regime, while increasing the strength of porous targets has much less influence on the value of . Therefore, one cannot say that non-porous targets are systematically easier or more difficult to disrupt than porous ones, as the outcome highly depends on the assumed strength values. In the gravity regime, we also confirm that the process of gravitational reaccumulation is at the origin of the largest remnant’s mass in both cases. We then propose some power-law relationships between and both target’s size and impact speed that can be used in collisional evolution models. The resulting fragment size distributions can also be reasonably fitted by a power-law whose exponent ranges between −2.2 and −2.7 for all target diameters in both cases and independently on the impact velocity (at least in the small range investigated between 3 and 5 km/s). Then, although ejection velocities in the gravity regime tend to be higher from porous targets, they remain on the same order as the ones from non-porous targets.     

On the population, physical decay and orbital distribution of Jupiter family comets: Numerical simulations

We study the Jupiter family comet (JFC) population assumed to come from the Scattered Disk and transferred to the Jupiter’s zone through gravitational interactions with the Jovian planets. We shall define as JFCs those with orbital periods and Tisserand parameters in the range 2<T?3.1, while those comets coming from the same source, but that do not fulfill the previous criteria (mainly because they have periods ) will be called ‘non-JFCs’. We performed a series of numerical simulations of fictitious comets with a purely dynamical model and also with a more complete dynamical-physical model that includes besides nongravitational forces, sublimation and splitting mechanisms. With the dynamical model, we obtain a poor match between the computed distributions of orbital elements and the observed ones. However with the inclusion of physical effects in the complete model we are able to obtain good fits to observations. The best fits are attained with four splitting models with a relative weak dependence on q, and a mass loss in every splitting event that is less when the frequency is high and vice versa. The mean lifetime of JFCs with radii and is found to be of about 150-200 revolutions (∼. The total population of JFCs with radii within Jupiter’s zone is found to be of 450±50. Yet, the population of non-JFCs with radii in Jupiter-crossing orbits may be ∼4 times greater, thus leading to a whole population of JFCs + non-JFCs of ∼2250±250. Most of these comets have perihelia close to Jupiter’s orbit. On the other hand, very few non-JFCs reach the Earth’s vicinity (perihelion distances ) which gives additional support to the idea that JFCs and Halley-type comets have different dynamical origins. Our model allows us to define the zones of the orbital element space in which we would expect to find a large number of JFCs. This is the first time, to our knowledge, that a physico-dynamical model is presented that includes sublimation and different splitting laws. Our work helps to understand the role played by these erosion effects in the distribution of the orbital elements and lifetimes of JFCs.     

Titan haze distribution and optical properties retrieved from recent observations

The VIMS instrument onboard Cassini observed the north polar region of Titan at 113° phase angle, 28 December 2006. On this spectral image, a vast polar cloud can be seen northward to 62°N, and elsewhere, the haze appears as the dominant source of scattering. Because the surface does not appear in the wavelength range between 0.3 and , this spectro-image is ideal to study airborn scatterers both in methane bands and windows. In this work, we study this image, along with another image taken at 13° phase angle. This image probe both the atmosphere and the surface from pole to pole. First, we characterise the spatial distribution of the haze layer above 100 km between 80°S and 70°N. We find a north south asymmetry with a haze opacity increasing by a factor 3 from the south pole to the equator, then a constant value up to about 30°N and a decrease of a factor 2 between 30°N and about 60°N. Beyond 60°N, we can see the influence of the north polar cloud, even in the band, but no polar haze accumulation. The fact that the north polar region is still in the polar night is a possible explanation. No haze accumulation is observed in the southern polar region. Secondly, we partly identify the origin of spectral features in the 2.8-μm methane window, which are found to be due to deuterated methane (CH₃D). This allows the analyse of this window and to retrieve the opacity of scatterers layer below 80 km (presumably made of aerosols and condensate droplets) between 35°N and 60°N. Finally, we constrained the values and the spectral behaviour of the imaginary part of the aerosol refractive index in the range between 0.3 and . To do so, we studied the 2.8-μm window with the image taken at 113° phase angle. To complete the analysis, we studied the transmission through the haze layer in the 3.4-μm band observed in solar occultation mode with VIMS, and we analysed the single scattering albedo retrieved with DISR instrument between 0.4 and . The imaginary part of the refractive index that we find for Titan aerosols follows Khare et al. (Khare, B.N. et al. [1984]. Icarus 60, 127-137) optical constant up to and becomes constant beyond this wavelength at least up to . It also has a prominent peak at and a secondary peak at , which indicates material rich in C-H bonds, with much less N-H bonds than in Khare et al. (1984) tholins.     

3782 Celle: Discovery of a binary system within the Vesta family of asteroids

Photometric observations of the minor planet (3782) Celle, which has been associated both dynamically and spectroscopically with the Vesta asteroid family, were obtained using the 1.8-m Vatican Advanced Technology Telescope during September 2001 and December 2002-January 2003. Analysis of these data reveals a normal rotational lightcurve (, amplitude =0.10-0.15 mag). During the 2002-2003 run, anomalous attenuation events were observed lasting for about 2.6-3.5 h that varied in amplitude from 0.15-0.3 mag. The attenuations were of two distinct types that can clearly be identified as primary and secondary occultation/eclipses similar to those that have been previously observed in known minor planet binary systems (Pravec et al., 2000). We therefore interpret our data as clear evidence that (3782) Celle is actually an asynchronous binary system with an orbital period of (Ryan et al., 2003). A preliminary model, based on spherical components, yields a primary-to-secondary diameter ratio of 0.43±0.01 and a combined bulk density of for the two components. Because these objects are likely to be composed of basaltic fragments, this density is indicative of a moderate to a highly fractured internal structure for at least one, if not both, of the binary components. Since the Vesta family is believed to have been created via a cratering event, this finding has important implications for understanding possible ejecta re-accumulation and satellite formation in subcatastrophic collisions.     

Thermal disturbances caused by lander shadowing and the measurability of the martian planetary heat flow

A measurement of the martian planetary heat flow requires the determination of the subsurface temperature gradient, which is affected by surface insolation. I investigate the propagation of thermal disturbances caused by lander shadowing and derive measurement requirements for in situ heat flow experiments. I find that for short term measurements spanning 180 sol, a measurement depth of at least 2 m is needed to guarantee a stable thermal environment directly underneath the lander for Moon-like thermal conductivities of . For extremely large conductivities of , this depth needs to be increased to 4 m, but if the probe can be deployed outside the lander structure, the respective depths can be decreased by 1 m. For long term measurements spanning at least a full martian year heat flow perturbations are smaller than 5% below a depth of 3 m directly underneath the lander. Outside the lander structure, essentially unperturbed measurements may be conducted at depths of 0.5 and 1.5 m for thermal conductivities of 0.02 and , respectively.     

Chemical and mineralogical size segregation in the impact disruption of inhomogeneous, anhydrous meteorites

We performed impact disruption experiments on pieces from eight different anhydrous chondritic meteorites—four weathered ordinary chondrite finds from North Africa (NWA791, NWA620, NWA869 and MOR001), three almost unweathered ordinary chondrite falls (Mbale, Gao, and Saratov), and an almost unweathered carbonaceous chondrite fall (Allende). In each case the impactor was a small (1/8 or 1/4 in) aluminum sphere fired at the meteorite target at , comparable to the mean collision speed in the main-belt. Some of the ∼5 to debris from each disruption was collected in aerogel capture cells, and the captured particles were analyzed by in situ synchrotron-based X-ray fluorescence. For each meteorite, many of the smallest particles ( up to in size, depending on the meteorite) exhibit very high Ni/Fe ratios compared to the Ni/Fe ratios measured in the larger particles , a composition consistent with the smallest debris being dominated by matrix material while the larger debris is dominated by fragments from olivine chondrules. These results may explain why the interplanetary dust particles (IDPs) collected from the Earth's stratosphere are C-rich and volatile-rich compared to the presumed solar nebula composition. The IDPs may simply sample the matrix of an inhomogeneous parent body, structurally and mineralogically similar to the chondritic meteorites, which are inhomogeneous assemblages of compact, strong, C- and volatile-poor chondrules that are distributed in a more porous, C- and volatile-rich matrix. In addition, these results may explain why the micrometeorites, which are to millimeters in size, recovered from the polar ices are Ni- and S-poor compared to chondritic meteorites, since these polar micrometeorites may preferentially sample fragments from the Ni- and S-poor olivine chondrules. These results indicate that the average composition of the IDPs may be biased towards the composition of the matrix of the parent body while the average composition of the polar micrometeorites may be more heavily weighted towards the composition of the chondrules and clasts. Thus, neither the IDPs nor the polar micrometeorites may sample the bulk composition of their respective parent bodies.We determined the threshold collisional specific energy for these chondritic meteorites to be 1419 J/kg, about twice the value for terrestrial basalt. Comparison of the mass of the largest fragment produced in the disruption of an sample of the porous ordinary chondrite Saratov with the largest fragment produced in the disruption of an sample of the compact ordinary chondrite MOR001 when each was struck by an impactor having approximately the same kinetic energy confirms that it requires significantly more energy to disrupt a porous target than a non-porous target.These results may also have important implications for the design of spacecraft missions intended to sample the composition and mineralogy of the chondritic asteroids and other inhomogeneous bodies. A Stardust-like spacecraft intended to sample asteroids by collecting only the small debris from a man-made impact onto the asteroid may collect particles that over-sample the matrix of the target and do not provide a representative sample of the bulk composition. The impact collection technique to be employed by the Japanese HAYABUSA (formerly MUSES-C) spacecraft to sample the asteroid Itokawa may result in similar mineral segregation.     

Calculation of the enrichment of the giant planet envelopes during the “late heavy bombardment”

The giant planets of our solar system possess envelopes consisting mainly of hydrogen and helium but are also significantly enriched in heavier elements relatively to our Sun. In order to better constrain how these heavy elements have been delivered, we quantify the amount accreted during the so-called “late heavy bombardment”, at a time when planets were fully formed and planetesimals could not sink deep into the planets. On the basis of the “Nice model”, we obtain accreted masses (in terrestrial units) equal to for Jupiter, and for Saturn. For the two other giant planets, the results are found to depend mostly on whether they switched position during the instability phase. For Uranus, the accreted mass is with an inversion and without an inversion. Neptune accretes in models in which it is initially closer to the Sun than Uranus, and otherwise. With well-mixed envelopes, this corresponds to an increase in the enrichment over the solar value of 0.033±0.001 and 0.074±0.007 for Jupiter and Saturn, respectively. For the two other planets, we find the enrichments to be 2.1±1.4 (w/ inversion) or 1.2±0.7 (w/o inversion) for Uranus, and 2.0±1.2 (w/ inversion) or 2.7±1.6 (w/o inversion) for Neptune. This is clearly insufficient to explain the inferred enrichments of ∼4 for Jupiter, ∼7 for Saturn and ∼45 for Uranus and Neptune.     

Thermal emission spectroscopy (5.2-38 μm) of three Trojan asteroids with the Spitzer Space Telescope: Detection of fine-grained silicates

We present thermal emission spectra (5.2-38 μm) of the Trojan asteroids 624 Hektor, 911 Agamemnon, and 1172 Aneas. The observations used the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. Emissivity spectra are created by dividing the measured Spectral Energy Distribution (SED) by a model of the thermal continuum. We employ the Standard Thermal Model (STM), allowing physical parameters (e.g., radius and albedo) to vary in order to find the best thermal continuum fit to the SED. The best-fit effective radius (R) and visible geometric albedo (pv) for Hektor (R=110.0±7.3, ) and Aneas (R=69.1±5.1, ) agree very well with previous estimates, and for Agamemnon (R=71.5±5.2, ) we find slightly a smaller size and higher albedo than previously derived. Other thermal models (e.g., thermophysical) result in estimates of R and pv that vary a few percent from the STM, but the resulting emissivity spectra are identical. The emissivity spectra of all three asteroids display an emissivity plateau near 10-μm and another broader rise from ∼18 to 28 μm. We interpret these as indications of fine-grained silicates on the surfaces of these asteroids. The emissivity spectra more closely resemble emission spectra from cometary comae than those from solid surfaces and measured in the laboratory for powdered meteorites and regolith analogs. We hypothesize that the coma-like emission from the solid surfaces of trojans may be due to small silicate grains being imbedded in a relatively transparent matrix, or to a very under-dense (fairy-castle) surface structure. These hypotheses need to be tested by further laboratory and theoretical scattering work as well as continued thermal emission observations of asteroids.     

Photochemical escape of oxygen from Mars: A comparison of the exobase approximation to a Monte Carlo method

The non-thermal escape of neutral O atoms from Mars at the current epoch is largely due to dissociative recombination of :

     

Venus wind map at cloud top level with the MTR/THEMIS visible spectrometer, I: Instrumental performance and first results

Solar light gets scattered at cloud top level in Venus’ atmosphere, in the visible range, which corresponds to the altitude of 67 km. We present Doppler velocity measurements performed with the high resolution spectrometer MTR of the Solar telescope THEMIS (Teide Observatory, Canary Island) on the sodium D2 solar line . Observations lasted only 49 min because of cloudy weather. However, we could assess the instrumental velocity sensitivity, per pixel of 1 arcsec, and give a value of the amplitude of zonal wind at equator at .     

Impact of atmospheric gravity waves on the jovian ionosphere

We investigate the effects of atmospheric gravity waves on the vertical and horizontal structure of the ionosphere of Jupiter. The presented non-linear, two-dimensional model of the jovian ionosphere allows for spatially and temporally varying neutral wind and temperature fields and tracks the time evolution of six ionospheric species, , and . An analytical approach is used to validate the model results for linear, small-amplitude waves and to elucidate the mechanisms that leads to perturbations in the density of the main ion species, H⁺ and . We demonstrate that the long-lived H⁺ ions are perturbed directly by wave dynamics whereas short-lived ions such as are perturbed by chemical interactions with other perturbed ion species. The model is then applied using larger gravity wave amplitudes consistent with observations. Atmospheric gravity waves propagating at high altitudes create layers of enhanced electron density similar to the system of layers observed during the J0-ingress radio occultation of the Galileo spacecraft. Our best fit to the J0-ingress observation is achieved using an 82 min period forcing wave with horizontal and vertical wavelengths of 500 km and 60 km respectively, and peaks at 510 km above the 1 bar pressure level. We further investigate the effects of the wave-induced ion flux on the background ionospheric structure and demonstrate that in the presence of a gravity wave the background density profiles of the H⁺ and ions are significantly modified. We also find that the column density of has variations that can exceed 10% as the wave propagates.     

Isotope fractionation in the photochemical escape of O from Mars

The production of energetic and escaping neutral O atoms at the current epoch in the martian thermosphere is thought to be dominated by the dissociative recombination process: