Space weathering of small Koronis family asteroids in the SDSS Moving Object Catalog

Following the work of Rivkin et al. (Rivkin, A.S., Thomas, C.A., Trilling, D.E., Enga, M., Grier, J.A. [2011]. Icarus 211, 1294–1297) and Thomas et al. (Thomas, C.A., Rivkin, A.S, Trilling, D.E., Enga, M., Grier, J.A. [2011a]. Icarus 212, 158–166), we investigate space weathering trends in the Koronis family using the larger sample size of the Sloan Digital Sky Survey Moving Object Catalog. We confirm the trend in spectral slope seen in our earlier work and extend our results by investigating the trend in band depth (i ? z color index) to show that Koronis family asteroids smaller than 4 km show the transition from fresh Q-type to weathered S-type surfaces.     

Rotation periods of binary asteroids with large separations – Confronting the Escaping Ejecta Binaries model with observations

Durda et al. (Durda, D.D., Bottke, W.F., Enke, B.L., Merline, W.J., Asphaug, E., Richardson, D.C., Leinhardt, Z.M. [2004]. Icarus 170, 243–257), using numerical models, suggested that binary asteroids with large separation, called Escaping Ejecta Binaries (EEBs), can be created by fragments ejected from a disruptive impact event. It is thought that six binary asteroids recently discovered might be EEBs because of the high separation between their components (～100 > a/R_p > ～20).However, the rotation periods of four out of the six objects measured by our group and others and presented here show that these suspected EEBs have fast rotation rates of 2.5–4 h. Because of the small size of the components of these binary asteroids, linked with this fast spinning, we conclude that the rotational-fission mechanism, which is a result of the thermal YORP effect, is the most likely formation scenario. Moreover, scaling the YORP effect for these objects shows that its timescale is shorter than the estimated ages of the three relevant Hirayama families hosting these binary asteroids. Therefore, only the largest (D ～ 19 km) suspected asteroid, (317) Roxane, could be, in fact, the only known EEB.In addition, our results confirm the triple nature of (3749) Balam by measuring mutual events on its lightcurve that match the orbital period of a nearby satellite in addition to its distant companion. Measurements of (1509) Esclangona at different apparitions show a unique shape of the lightcurve that might be explained by color variations.     

Absolute magnitudes of asteroids and a revision of asteroid albedo estimates from WISE thermal observations

We obtained estimates of the Johnson V absolute magnitudes (H) and slope parameters (G) for 583 main-belt and near-Earth asteroids observed at Ond?ejov and Table Mountain Observatory from 1978 to 2011. Uncertainties of the absolute magnitudes in our sample are <0.21 mag, with a median value of 0.10 mag. We compared the H data with absolute magnitude values given in the MPCORB, Pisa AstDyS and JPL Horizons orbit catalogs. We found that while the catalog absolute magnitudes for large asteroids are relatively good on average, showing only little biases smaller than 0.1 mag, there is a systematic offset of the catalog values for smaller asteroids that becomes prominent in a range of H greater than ～10 and is particularly big above H ～ 12. The mean (H_catalog ? H) value is negative, i.e., the catalog H values are systematically too bright. This systematic negative offset of the catalog values reaches a maximum around H = 14 where the mean (H_catalog ? H) is ?0.4 to ?0.5. We found also smaller correlations of the offset of the catalog H values with taxonomic types and with lightcurve amplitude, up to ～0.1 mag or less. We discuss a few possible observational causes for the observed correlations, but the reason for the large bias of the catalog absolute magnitudes peaking around H = 14 is unknown; we suspect that the problem lies in the magnitude estimates reported by asteroid surveys. With our photometric H and G data, we revised the preliminary WISE albedo estimates made by Masiero et al. (Masired, J.R. et al. [2011]. Astrophys. J. 741, 68–89) and Mainzer et al. (Mainzer, A. et al. [2011b]. Astrophys. J. 743, 156–172) for asteroids in our sample. We found that the mean geometric albedo of Tholen/Bus/DeMeo C/G/B/F/P/D types with sizes of 25–300 km is p_V = 0.057 with the standard deviation (dispersion) of the sample of 0.013 and the mean albedo of S/A/L types with sizes 0.6–200 km is 0.197 with the standard deviation of the sample of 0.051. The standard errors of the mean albedos are 0.002 and 0.006, respectively; systematic observational or modeling errors can predominate over the quoted formal errors. There is apparent only a small, marginally significant difference of 0.031 ± 0.011 between the mean albedos of sub-samples of large and small (divided at diameter 25 km) S/A/L asteroids, with the smaller ones having a higher albedo. The difference will have to be confirmed and explained; we speculate that it may be either a real size dependence of surface properties of S type asteroids or a small size-dependent bias in the data (e.g., a bias towards higher albedos in the optically-selected sample of asteroids). A trend of the mean of the preliminary WISE albedo estimates increasing with asteroid size decreasing from D ～ 30 down to ～5 km (for S types) showed in Mainzer et al. (Mainzer, A. et al. [2011a]. Astrophys. J. 741, 90–114) appears to be mainly due to the systematic bias in the MPCORB absolute magnitudes that progressively increases with H in the corresponding range H = 10–14.     

UV spectrum of Enceladus

We present a far ultraviolet (FUV) spectrum of Saturn’s moon Enceladus from the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST). We have put upper limits on emission from C, N, and O lines in Enceladus’ atmosphere and column densities for the C lines assuming solar resonance scattering. We find these upper limits to be relatively low—on the order of tens to thousands of Rayleighs and with C column densities on the order of 10⁸–10¹⁵ cm^?2, depending on the assumed source size. We also present a segment of a reflectance spectrum in the FUV from ～1900–2130 Å. This region was sensitive to the different ice mixtures in the model spectra reported by Hendrix et al. (Hendrix, A.R., Hansen, C.J., Holsclaw, G.M. [2010]. Icarus, 206, 608). We find the spectrum brightens quickly longward of ～1900 Å, constraining the absorption band observed by Hendrix et al. from ～170 to 190 nm. We find our data is consistent with the suggestion of Hendrix et al. of the presence of ammonia ice (or ammonia hydrate) to darken that region, and also possibly tholins to darken the mid-UV, as reported by Verbiscer et al. (Verbiscer, A.J., French, R.G., McGhee, C.A. [2005]. Icarus, 173, 66).     

A survey of small fast rotating asteroids among the near-Earth asteroid population

A survey of 62 small near-Earth asteroids was conducted to determine the rotation state of these objects and to search for rapid rotation. Since results for 9 of the asteroids were previously published (Pravec, P., Hergenrother, C.W., Whiteley, R.J., Šarounová, L., Kušnirák, P., Wolf, M. [2000]. Icarus 147, 477-486; Pravec, P. et al. [2005] Icarus 173, 108-131; Whiteley, R.J., Tholen, D.J., Hergenrother, C.W. [2002a]. Icarus 157, 139-154; Hergenrother, C.W., Whiteley, R.J., Christensen, E.J. [2009]. Minor Planet Bull. 36, 16-18.), this paper will present results for the remaining 53 objects. Rotation periods significantly less than 2 h are indicative of intrinsic strength in the asteroids, while periods longer than 2 h are typically associated with gravitationally bound aggregates. Asteroids with absolute magnitude (H) values ranging from 20.4 to 27.4 were characterized. The slowest rotator with a definite period is 2004 BW₁₈ with a period of 8.3 h, while 2000 DO₈ and 2000 WH₁₀ are the fastest with periods of 1.3 min. A minimum of two-thirds of asteroids with H > 20 are fast rotating and have periods significantly faster than 2.0 h. The percentage of rapid rotators increases with decreasing size and a minimum of 79% of H ? 24 objects are rapid rotators. Slowly-rotating objects, some with periods as long as 10-20 h, make up a small though significant fraction of the small asteroid population. There are three fast rotators with relatively large possible diameters (D): 2001 OE₈₄ with 470 ? D ? 820 m (Pravec, P., Kušnirák, P., Šarounová, L., Harris, A.W., Binzel, R.P., Rivkin, A.S. [2002b]. Large coherent Asteroid 2001 OE84. In: Warmbein, B. (Eds.), Proceedings of Asteroids, Comets, Meteors - ACM 2002. Springer, Berlin, pp. 743-745), 2001 FE₉₀ with 265 ? D ? 594 m (Hicks, M., Lawrence, K., Rhoades, H., Somers, J., McAuley, A., Barajas, T. [2009]. The Astronomer’s Telegrams, # 2116), and 2001 VF₂ with a possible D of 145 ? D ? 665 m. Using the diameters derived from nominal absolute magnitudes and albedos, the remainder of the fast rotating population is completely consistent with D ? 200 m. Even when taking into account the largest possible uncertainties in the determination of diameters, the remainder must all have D ? 400 m. With the exceptions of 2001 OE₈₄, this result agrees with previous upper diameter limits for fast rotators in Pravec and Harris (Pravec, P., Harris, A.W. [2000]. Icarus 148, 589-593) and Whiteley et al. (Whiteley, R.J, Tholen, D.J., Hergenrother, C.W. [2002a]. Icarus 157, 139-154.     

Constraining albedo,diameter and composition of near-Earth asteroids via near-infrared spectroscopy

We present a method to constrain the albedo and diameters of near-Earth asteroids (NEAs) based on thermal flux in their near-infrared spectra (0.7–2.5 μm) using the Standard Thermal Model. Near-infrared spectra obtained with the SpeX instrument on NASA Infrared Telescope Facility are used to estimate the albedo and diameters of 12 NEAs (1992 JE, 1992 UY4, 1999 JD6, 2004 XP14, 2005 YY93, 2007 DS84, 2005 AD13, 2005 WJ56, 1999 JM8, 2005 RC34, 2003 YE45, and 2008 QS11). Albedo estimates were compared with average albedo for various taxonomic classes outlined by Thomas et al. (Thomas, C.A. et al. [2011]. Astron. J. 142(3)) and are consistent with their results. Spectral band parameters, like band centers, are derived and compared to spectra of laboratory mineral mixtures and meteorites to constrain their composition and possible meteorite analogs. Based on our study we estimate the albedos and diameters of these NEAs and compare them with those obtained by other techniques such as ground-based mid-infrared, Spitzer thermal infrared and Arecibo radar. Our results are broadly consistent with the results from other direct methods like radar. Determining the compositions of low albedo asteroids is a challenge due to the lack of deep silicate absorption features. However, based on weak absorption features and albedo, we suggest possible meteorite analogs for these NEAs, which include black chondrites, CM2 carbonaceous chondrites and enstatite achondrites. We did not find any specific trends in albedo and composition among the asteroids we observed.     

Space weathering of small Koronis family members

The space weathering process and its implications for the relationships between S- and Q-type asteroids and ordinary chondrite meteorites is an often debated topic in asteroid science. Q-type asteroids have been shown to display the best spectral match to ordinary chondrites (McFadden, L.A., Gaffey, M.J., McCord, T.B. [1985]. Science 229, 160–163). While the Q-types and ordinary chondrites share some spectral features with S-type asteroids, the S-types have significantly redder spectral slopes than the Q-types in visible and near-infrared wavelengths. This reddening of spectral slope is attributed to the effects of space weathering on the observed surface composition. The analysis by Binzel et al. (Binzel, R.P., Rivkin, A.S., Stuart, J.S., Harris, A.W., Bus, S.J., Burbine, T.H. [2004]. Icarus 170, 259–294) provided a missing link between the Q- and S-type bodies in near-Earth space by showing a reddening of spectral slope in objects from 0.1 to 5 km that corresponded to a transition from Q-type to S-type asteroid spectra, implying that size, and therefore surface age, is related to the relationship between S- and Q-types. The existence of Q-type asteroids in the main-belt was not confirmed until Mothé-Diniz and Nesvorny (Mothé-Diniz, T., Nesvorny, D. [2008]. Astron. Astrophys. 486, L9–L12) found them in young S-type clusters. The young age of these families suggest that the unweathered surface could date to the formation of the family. This leads to the question of whether older S-type main-belt families can contain Q-type objects and display evidence of a transition from Q- to S-type. To answer this question we have carried out a photometric survey of the Koronis family using the Kitt Peak 2.1 m telescope. This provides a unique opportunity to compare the effects of the space weathering process on potentially ordinary chondrite-like bodies within a population of identical initial conditions. We find a trend in spectral slope for objects 1–5 km that shows the transition from Q- to S-type in the main-belt. This data set will prove crucial to our understanding of the space weathering process and its relevant timescales.     

Ground-based near-infrared observations of water vapour in the Venus troposphere

We present a study of water vapour in the Venus troposphere obtained by modelling specific water vapour absorption bands within the 1.18 μm window. We compare the results with the normal technique of obtaining the abundance by matching the peak of the 1.18 μm window. Ground-based infrared imaging spectroscopy of the night side of Venus was obtained with the Anglo-Australian Telescope and IRIS2 instrument with a spectral resolving power of R ～ 2400. The spectra have been fitted with modelled spectra simulated using the radiative transfer model VSTAR. We find a best fit abundance of 31 ppmv (?6 +9 ppmv), which is in agreement with recent results by Bézard et al. (Bézard, B., Fedorova, A., Bertaux, J.-L., Rodin, A., Korablev, O. [2011]. Icarus, 216, 173–183) using VEX/SPICAV (R ～ 1700) and contrary to prior results by Bézard et al. (Bézard, B., de Bergh, C., Crisp, D., Maillard, J.P. [1990]. Nature, 345, 508–511) of 44 ppmv (±9 ppmv) using VEX/VIRTIS-M (R ～ 200) data analyses. Comparison studies are made between water vapour abundances determined from the peak of the 1.18 μm window and abundances determined from different water vapour absorption features within the near infrared window. We find that water vapour abundances determined over the peak of the 1. 18 μm window results in plots with less scatter than those of the individual water vapour features and that analyses conducted over some individual water vapour features are more sensitive to variation in water vapour than those over the peak of the 1. 18 μm window. No evidence for horizontal spatial variations across the night side of the disk are found within the limits of our data with the exception of a possible small decrease in water vapour from the equator to the north pole. We present spectral ratios that show water vapour absorption from within the lowest 4 km of the Venus atmosphere only, and discuss the possible existence of a decreasing water vapour concentration towards the surface.     

Thermal structure and CO distribution for the Venus mesosphere/lower thermosphere: 2001–2009 inferior conjunction sub-millimeter CO absorption line observations

Sub-millimeter ¹²CO (346 GHz) and ¹³CO (330 GHz) line absorptions, formed in the mesosphere and lower thermosphere of Venus (70–120 km), have been mapped across the nightside Venus disk during 2001–2009 inferior conjunctions, employing the James Clerk Maxwell Telescope (JCMT). Radiative transfer analysis of these thermal line absorptions supports temperature and CO mixing profile retrievals, as well as Doppler wind fields (described in the companion paper, Clancy et al., 2012). Temporal sampling over the hourly, daily, weekly and interannual timescales was obtained over 2001–2009. On timescales inferred as several weeks, we observe changes between very distinctive CO and temperature nightside distributions. Retrieved nightside CO, temperature distributions for January 2006 and August 2007 observations display strong local time, latitudinal gradients consistent with early morning (2–3 am), low-to-mid latitude (0–40NS) peaks of 100–200% in CO and 20–30 K in temperature. The temperature increases are most pronounced above 100 km altitudes, whereas CO variations extend from 105 km (top altitude of retrieval) down to below 80 km in the mesosphere. In contrast, the 2004 and 2009 periods of observation display modest temperature (5–10 K) and CO (30–60%) increases, that are centered on antisolar (midnight) local times and equatorial latitudes. Doppler wind derived global (zonal and should be SSAS) circulations from the same data do not exhibit variations correlated with these CO, temperature short-term variations. However, large-scale residual wind fields not fit by the zonal, SSAS circulations are observed in concert with the strong temperature, CO gradients observed in 2006 and 2007 (Clancy et al., 2010). These short term variations in nightside CO, temperature distributions may also be related to observed nightside variations in O₂ airglow (Hueso, H., Sánchez-Lavega, A., Piccioni, G., Drossart, P., Gérard, J.C., Khatuntsev, I., Zasova, L., Migliorini, A. [2008]. J. Geophys. Res. 113, E00B02. doi:10.1029/2008JE003081) and upper mesospheric SO and SO₂ layers (Sandor, B.J., Clancy, R.T., Moriarty-Schieven, G.H., Mills, F.P. [2010]. Icarus 208, 49–60).The retrieved temperature profiles also exhibit 20 K long-term (2001–2009) variations in nightside (whole disk) average mesospheric (80–95 km) temperatures, similar to 1982–1991 variations identified in previous millimeter CO line observations (Clancy et al., 1991). Global average diurnal variations in lower thermospheric temperatures and mesospheric CO abundances decreased by a factor-of-two between 2000–2002 versus 2007–2009 periods of combined dayside and nightside observations. The infrequency and still limited temporal extent of the observations make it difficult to assign specific timescales to such longer term variations, which may be associated with longer term variations observed for cloud top SO₂ (Esposito, L.W., Bertaux, J.-L., Krasnopolsky, V., Moroz, V.I., Zasova, L.V. [1997]. Chemistry of lower atmosphere and clouds. In: Bougher, S.W., Hunten, D.M., Phillips, R.J. (Eds.), VENUS II, 1362pp) and mesospheric water vapor (Sandor, B.J., Clancy, R.T. [2005]. Icarus 177, 129–143) abundances.     

Main Belt Comet P/2008 R1 Garradd: Duration of activity

In this work we investigated changes of the water emission from a model comet of the size and orbital elements of Comet P/2008 R1 (Garradd). We performed simulations for model cometary nuclei of different compositions and two different orientations in space. Our simulations indicate, that the emission of water decreases from one orbital period to another one, but in some cases slowly. When the rotation axis of the nucleus lies in the orbital plane the seasonal maximum of water production during the first two orbital periods can be as high as about 10²⁶ mol s^?1, but decreases by two orders of magnitude during only 50 orbital periods. The highest rate of water production after many orbital periods is expected when the rotation axis is perpendicular to the orbital plane – the seasonal maximum of water production can be about 5 × 10²⁵ mol s^?1 during the first two orbital periods after activation of the comet and no more than 0.8 × 10²⁵ mol s^?1 500 orbital periods later. The upper estimate for the production of water derived from observations of P/2008 R1 (Garradd) by Jewitt et al. (Jewitt, D., Yang, B., Haghighipour, N. [2009]. Astron. J. 137, 4313–4321) is 5 × 10²⁵.     

Near-infrared spectroscopic survey of B-type asteroids: Compositional analysis

We present near-infrared spectra of 23 B-type asteroids obtained with the NICS camera-spectrograph at the 3.56 m Telescopio Nazionale Galileo. We also compile additional visible and near-infrared spectra of another 22 B-type asteroids from the literature. A total of 45 B-types are analyzed. No significant trends in orbital properties of our sample were detected when compared with all known B-types and all known asteroids. The reflectance spectra of the asteroids in the 0.8–2.5 μm range show a continuous shape variation, from a monotonic negative (blue) slope to a positive (red) slope. This continuous spectral trend is filling the gap between the two main groups of B-types published by Clark et al. ([2010]. J. Geophys. Res., 115, 6005–6027). We found no clear correlations between the spectral slope and the asteroids’ sizes or heliocentric distances. We apply a clustering technique to reduce the volume of data to six optimized “average spectra” or “centroids”, representative of the whole sample. These centroids are then compared against meteorite spectra from the RELAB database. We found carbonaceous chondrites as the best meteorite analogs for the six centroids. There is a progressive change in analogs that correlates with the spectral slope: from CM2 chondrites (water-rich, aqueously altered) for the reddest centroid, to CK4 chondrites (dry, heated/thermally altered) for the bluest one.     

Spin rate distribution of small asteroids

The spin rate distribution of main belt/Mars crossing (MB/MC) asteroids with diameters 3-15 km is uniform in the range from f=1 to 9.5 d⁻¹, and there is an excess of slow rotators with f<1 d⁻¹. The observed distribution appears to be controlled by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. The magnitude of the excess of slow rotators is related to the residence time of slowed down asteroids in the excess and the rate of spin rate change outside the excess. We estimated a median YORP spin rate change of ≈0.022 d⁻¹/Myr for asteroids in our sample (i.e., a median time in which the spin rate changes by 1 d⁻¹ is ≈45 Myr), thus the residence time of slowed down asteroids in the excess is ≈110 Myr. The spin rate distribution of near-Earth asteroids (NEAs) with sizes in the range 0.2-3 km (∼5 times smaller in median diameter than the MB/MC asteroids sample) shows a similar excess of slow rotators, but there is also a concentration of NEAs at fast spin rates with f=9-10 d⁻¹. The concentration at fast spin rates is correlated with a narrower distribution of spin rates of primaries of binary systems among NEAs; the difference may be due to the apparently more evolved population of binaries among MB/MC asteroids.     

Investigating the Vesta-vestoid-HED connection

We present spectral observations of Minor Planet 4 Vesta, of five V-type asteroids which are physically near Vesta, and of two V-type NEAs. We use these spectra to determine the presence or absence of a weak feature at 506.5 nm which is indicative of the presence of spin-forbidden Fe²⁺ in sixfold coordination. As with our earlier observations [Cochran and Vilas, Icarus 134 (1998) 207-212], we find this feature at all observed rotational phases of Vesta and again see the trend that spectra at longitudes between 240° and 360° have a smaller 506.5 nm feature equivalent width than spectra obtained at other longitudes. Additionally, we searched for this feature in V-class main-belt and NEA asteroids and positively detected the feature in main-belt Asteroid 2579 Spartacus and possibly in 3376 Armandhammer. The other objects lacked the feature. Our results are compared with previous observations of this feature by Vilas et al. [Icarus 147 (2000) 119-128]. The spatial distribution of the bodies as a function of the presence of this feature was investigated. We discuss the implication of the presence of this feature and the depth of the 0.9 μm pyroxene band for the scenario that pieces of Vesta were transported, via the 3:1 and ν₆ resonances, to the NEAs, and thence to inclusion in our meteorite collections as HED meteorites.     

Radiative transfer analyses of Titan’s tropical atmosphere

Titan’s optical and near-IR spectra result primarily from the scattering of sunlight by haze and its absorption by methane. With a column abundance of 92 km amagat (11 times that of Earth), Titan’s atmosphere is optically thick and only ～10% of the incident solar radiation reaches the surface, compared to 57% on Earth. Such a formidable atmosphere obstructs investigations of the moon’s lower troposphere and surface, which are highly sensitive to the radiative transfer treatment of methane absorption and haze scattering. The absorption and scattering characteristics of Titan’s atmosphere have been constrained by the Huygens Probe Descent Imager/Spectral Radiometer (DISR) experiment for conditions at the probe landing site (Tomasko, M.G., Bézard, B., Doose, L., Engel, S., Karkoschka, E. [2008a]. Planet. Space Sci. 56, 624–247; Tomasko, M.G. et al. [2008b]. Planet. Space Sci. 56, 669–707). Cassini’s Visual and Infrared Mapping Spectrometer (VIMS) data indicate that the rest of the atmosphere (except for the polar regions) can be understood with small perturbations in the high haze structure determined at the landing site (Penteado, P.F., Griffith, C.A., Tomasko, M.G., Engel, S., See, C., Doose, L., Baines, K.H., Brown, R.H., Buratti, B.J., Clark, R., Nicholson, P., Sotin, C. [2010]. Icarus 206, 352–365). However the in situ measurements were analyzed with a doubling and adding radiative transfer calculation that differs considerably from the discrete ordinates codes used to interpret remote data from Cassini and ground-based measurements. In addition, the calibration of the VIMS data with respect to the DISR data has not yet been tested. Here, VIMS data of the probe landing site are analyzed with the DISR radiative transfer method and the faster discrete ordinates radiative transfer calculation; both models are consistent (to within 0.3%) and reproduce the scattering and absorption characteristics derived from in situ measurements. Constraints on the atmospheric opacity at wavelengths outside those measured by DISR, that is from 1.6 to 5.0 μm, are derived using clouds as diffuse reflectors in order to derive Titan’s surface albedo to within a few percent error and cloud altitudes to within 5 km error. VIMS spectra of Titan at 2.6–3.2 μm indicate not only spectral features due to CH₄ and CH₃D (Rannou, P., Cours, T., Le Mouélic, S., Rodriguez, S., Sotin, C., Drossart, P., Brown, R. [2010]. Icarus 208, 850–867), but also a fairly uniform absorption of unknown source, equivalent to the effects of a darkening of the haze to a single scattering albedo of 0.63 ± 0.05. Titan’s 4.8 μm spectrum point to a haze optical depth of 0.2 at that wavelength. Cloud spectra at 2 μm indicate that the far wings of the Voigt profile extend 460 cm^?1 from methane line centers. This paper releases the doubling and adding radiative transfer code developed by the DISR team, so that future studies of Titan’s atmosphere and surface are consistent with the findings by the Huygens Probe. We derive the surface albedo at eight spectral regions of the 8 × 12 km² area surrounding the Huygens landing site. Within the 0.4–1.6 μm spectral region our surface albedos match DISR measurements, indicating that DISR and VIMS measurements are consistently calibrated. These values together with albedos at longer 1.9–5.0 μm wavelengths, not sampled by DISR, resemble a dark version of the spectrum of Ganymede’s icy leading hemisphere. The eight surface albedos of the landing site are consistent with, but not deterministic of, exposed water ice with dark impurities.     

Obliquity variations of a moonless Earth

We numerically explore the obliquity (axial tilt) variations of a hypothetical moonless Earth. Previous work has shown that the Earth’s Moon stabilizes Earth’s obliquity such that it remains within a narrow range, between 22.1° and 24.5°. Without lunar influence, a frequency map analysis by Laskar et al. (Laskar, J., Joutel, F., Robutel, P. [1993]. Nature 361, 615–617) showed that the obliquity could vary between 0° and 85°. This has left an impression in the astrobiology community that a big moon is necessary to maintain a habitable climate on an Earth-like planet. Using a modified version of the orbital integrator mercury, we calculate the obliquity evolution for moonless Earths with various initial conditions for up to 4 Gyr. We find that while obliquity varies significantly more than that of the actual Earth over 100,000 year timescales, the obliquity remains within a constrained range, typically 20–25° in extent, for timescales of hundreds of millions of years. None of our Solar System integrations in which planetary orbits behave in a typical manner show obliquity accessing more than 65% of the full range allowed by frequency-map analysis. The obliquities of moonless Earths that rotate in the retrograde direction are more stable than those of prograde rotators. The total obliquity range explored for moonless Earths with rotation periods less than 12 h is much less than that for slower-rotating moonless Earths. A large moon thus does not seem to be needed to stabilize the obliquity of an Earth-like planet on timescales relevant to the development of advanced life.     

The remarkable surface homogeneity of the Dawn mission target (1) Ceres

Dwarf-planet (1) Ceres is one of the two targets, along with (4) Vesta, that will be studied by the NASA Dawn spacecraft via imaging, visible and near-infrared spectroscopy, and gamma-ray and neutron spectroscopy. While Ceres’ visible and near-infrared disk-integrated spectra have been well characterized, little has been done about quantifying spectral variations over the surface. Any spectral variation would give us insights on the geographical variation of the composition and/or the surface age. The only work so far was that of Rivkin and Volquardsen ([2010], Icarus 206, 327) who reported rotationally-resolved spectroscopic (disk-integrated) observations in the 2.2–4.0 μm range; their observations showed evidence for a relatively uniform surface.Here, we report disk-resolved observations of Ceres with SINFONI (ESO VLT) in the 1.17–1.32 μm and 1.45–2.35 μm wavelength ranges. The observations were made under excellent seeing conditions (0.6″), allowing us to reach a spatial resolution of ～75 km on Ceres’ surface. We do not find any spectral variation above a 3% level, suggesting a homogeneous surface at our spatial resolution. Slight variations (about 2%) of the spectral slope are detected, geographically correlated with the albedo markings reported from the analysis of the HST and Keck disk-resolved images of Ceres (Li et al. [2006], Icarus 182, 143; Carry et al. [2008], Astron. Astrophys. 478, 235). Given the lack of constraints on the surface composition of Ceres, however, we cannot assert the causes of these variations.     

Rebuttal to the comment by Malhotra and Strom on “Constraints on the source of lunar cataclysm impactors”

?uk et al. (?uk, M. Gladman, B.J., Stewart, S.T. [2010]. Icarus 207 590-594) concluded that the the lunar cataclysm (late heavy bombardment) was recorded in lunar Imbrian era craters, and that their size distribution is different from that of main belt asteroids (which may have been the dominant pre-Imbrian impactors). This result would likely preclude the asteroid belt as the direct source of lunar cataclysm impactors. Malhotra and Strom (Malhotra, R., Strom, R.G. [2011]. Icarus) maintain that the lunar impactor population in the Imbrian era was the same as in Nectarian and pre-Nectarian periods, and this population had a size distribution identical to that of main belt asteroids. In support of this claim, they present an Imbrian size distribution made from two data sets published by Wilhelms et al. (Wilhelms, D.E., Oberbeck, V.R., Aggarwal, H.R. [1978]. Proc. Lunar Sci. Conf. 9, 3735-3762). However, these two data sets cannot be simply combined as they represent areas of different ages and therefore crater densities. Malhotra and Strom (Malhotra, R., Strom, R.G. [2011]. Icarus) differ with the main conclusion of Wilhelms et al. (Wilhelms, D.E., Oberbeck, V.R., Aggarwal, H.R. [1978]. Proc. Lunar Sci. Conf. 9, 3735-3762) that the Nectarian and Imbrian crater size distributions were different. We conclude that the available data indicate that the lunar Imbrian-era impactors had a different size distribution from the older ones, with the Imbrian impactor distribution being significantly richer in small impactors than that of older lunar impactors or current main-belt asteroids.     

Comet 17P/Holmes: Possibility of a CO driven explosion

This work is a continuation of our previous paper about brightening of Comet 17P/Holmes (Kossacki, K.J., Szutowicz, S. [2010]. Icarus 207, 320–340). In that paper we presented results of simulations indicating that the nonuniform crystallization of amorphous water ice itself is probably not sufficient for an explosion. In the present work we investigate the possibility that the explosion is caused by a rapid sublimation of the CO ice leading to the rise of gas pressure above the tensile strength of the nucleus. We simulated evolution of a model nucleus in the orbit of Comet 17P/Holmes. The nucleus is composed of water ice, carbon monoxide ice and dust and has the shape of an elongated ellipsoid. The simulations include crystallization of amorphous ice in the nucleus, changes of the dust mantle thickness, and sublimation of the CO ice. In our model CO is mantling grains composed of dust and amorphous water ice. Orientation of the nuclear spin axis in space is the same as derived in Moreno et al. (Moreno, F., Ortiz, J.L., Santos-Sanz, P., Morales, N., Vidal-Nunez, M.J., Lara, L.M., Gutierrez, P.J. [2008]. Astrophys. J. 677, L63–L66) for Comet Holmes during recent brightening event. Hence, the angle between the orbital and the equatorial planes of the comet is I = 95°, and the cometocentric solar longitude at perihelion is Φ = 210°. The calculations are performed for the south pole being the sub-solar point close to time of the outburst. Our computations indicate, that the CO pressure within the comet nucleus can rise to high values. When the layer between the dust mantle and the crystallization front of the amorphous water ice is very fine grained, few microns in radius, the CO pressure within the nucleus can exceed 10 kPa. This value is the lowest estimate for the tensile strength of the nucleus of Comet Holmes (Reach, W.T., Vaubaillon, J., Lisse, C.M., Holloway, M., Rho, J. [2010]. Icarus 208, 276–292). Hence, when the gas pressure reaches this value the nucleus may explode.     

Photometric,spectral phase and temperature effects on 4 Vesta and HED meteorites: Implications for the Dawn mission

Phase angle and temperature are two important parameters that affect the photometric and spectral behavior of planetary surfaces in telescopic and spacecraft data. We have derived photometric and spectral phase functions for the Asteroid 4 Vesta, the first target of the Dawn mission, using ground-based telescopes operating at visible and near-infrared wavelengths (0.4–2.5 μm). Photometric lightcurve observations of Vesta were conducted on 15 nights at a phase angle range of 3.8–25.7° using duplicates of the seven narrowband Dawn Framing Camera filters (0.4–1.0 μm). Rotationally resolved visible (0.4–0.7 μm) and near-IR spectral observations (0.7–2.5 μm) were obtained on four nights over a similar phase angle range. Our Vesta photometric observations suggest the phase slope is between 0.019 and 0.029 mag/deg. The G parameter ranges from 0.22 to 0.37 consistent with previous results (e.g., Lagerkvist, C.-I., Magnusson, P., Williams, I.P., Buontempo, M.E., Argyle, R.W., Morrison, L.V. [1992]. Astron. Astrophys. Suppl. Ser. 94, 43–71; Piironen, J., Magnusson, P., Lagerkvist, C.-I., Williams, I.P., Buontempo, M.E., Morrison, L.V. [1997]. Astron. Astrophys. Suppl. Ser. 121, 489–497; Hasegawa, S. et al. [2009]. Lunar Planet. Sci. 40. ID 1503) within the uncertainty. We found that in the phase angle range of 0° < α ? 25° for every 10° increase in phase angle Vesta’s visible slope (0.5–0.7 μm) increases 20%, Band I and Band II depths increase 2.35% and 1.5% respectively, and the BAR value increase 0.30. Phase angle spectral measurements of the eucrite Moama in the lab show a decrease in Band I and Band II depths and BAR from the lowest phase angle 13° to 30°, followed by possible small increases up to 90°, and then a dramatic drop between 90° and 120° phase angle. Temperature-induced spectral effects shift the Band I and II centers of the pyroxene bands to longer wavelengths with increasing temperature. We have derived new correction equations using a temperature series (80–400 K) of HED meteorite spectra that will enable interpretation of telescopic and spacecraft spectral data using laboratory calibrations at room temperature (300 K).     

Upper limits for H2SO4 in the mesosphere of Venus

Rapid temporal variability of SO₂ and SO in the Venus 85–100 km mesosphere (Sandor, B.J., Clancy, R.T., Moriarty-Schieven G.H. [2007]. Bull. Am. Astron. Soc. 39, 503; Sandor, B.J., Clancy, R.T., Moriarty-Schieven, G.H., Mills, F.P. [2010]. Icarus 208, 49–60) requires in situ sources and sinks for these molecules. While many loss mechanisms are recognized, no process for in situ production is known. Observational investigations to find, or constrain other potential sulfur reservoirs offer one method toward understanding the applicable photochemistry. Here, we report upper limits for gas-phase H₂SO₄ (sulfuric acid) abundances in Venus’ 85–100 km upper mesosphere, derived from 16 ground-based sub-mm spectroscopic observations in the period 2004–2008. Unlike the ubiquitous sulfuric acid solid/liquid aerosol, the gas phase would be photochemically active, potentially both source and sink for SO and SO₂. H₂SO₄ is retrieved from sub-mm lines located in the same bandpass as the SO₂ and SO lines described by Sandor et al. (Sandor, B.J., Clancy, R.T., Moriarty-Schieven, G.H., Mills, F.P. [2010]. Icarus 208, 49–60). H₂SO₄ upper limits reported here are thus simultaneous and spatially coincident with measurements of SO₂ and SO, providing for analysis of the three sulfur species collectively. The average H₂SO₄ abundance over 16 observations is 1 ± 2 ppb (i.e. <3 ppb). Upper limits for individual observations range from 3 to 44 ppb, where quality of the observing weather is the dominant constraint on measurement precision. The sum of H₂SO₄, SO₂ and SO varies widely. In one comparison, the sum [H₂SO₄ + SO₂ + SO] measured on one date differs by 10-σ from the sum measured 2 months later. We conclude that upper mesospheric sulfur atoms are not conserved among the three molecules, that H₂SO₄ is not a significant sulfur reservoir for balancing the observed variations of [SO₂ + SO], and is not relevant to the (still unknown) photochemistry responsible for observed behavior of SO₂ and SO. Having ruled out H₂SO₄, we infer that elemental sulfur is the most probable candidate for the needed third reservoir.