Radar observations of E-class Asteroids 44 Nysa and 434 Hungaria

We observed the E-class main-belt Asteroids (MBAs) 44 Nysa and 434 Hungaria with Arecibo Observatory's S-band (12.6 cm) radar. Both asteroids exhibit polarization ratios higher than those measured for any other MBA: Nysa, μc=0.50±0.02 and Hungaria, μc=0.8±0.1. This is consistent with the high polarization ratios measured for every E-class near-Earth asteroid (NEA) observed by Benner et al. [Benner, L.A.M., and 10 collegues, 2008. Icarus, submitted for publication] and suggests a common cause. Our estimates of radar albedo are 0.19±0.06 for Nysa and 0.22±0.06 for Hungaria. These values are higher than those of most MBAs and, when combined with their high polarization ratios, suggest that the surface bulk density of both asteroids is high. We model Nysa as an ellipsoid of dimension 113×67×65 km (±15%) giving an effective diameter Deff=79±10 km, consistent with previous estimates. The echo waveforms are not consistent with a contact binary as suggested by Kaasalainen et al. [Kaasalainen, M., Torppa, J., Piironen, J., 2002. Astron. Astrophys. 383, L19-L22]. We place a constraint on Hungaria's maximum diameter, Dmax?11 km consistent with previous size estimates.     

Physical properties of Asteroid (25143) Itokawa—Target of the Hayabusa sample return mission

We present results of a ground-based observational study of the Hayabusa mission target near-Earth Asteroid (25143) Itokawa. Our data consist of BVRI-filter CCD photometry and low resolution CCD spectroscopy, from which the asteroid's rotation period, axial ratio, broadband colors, and taxonomic classification are derived. Analysis of the R-filter lightcurve data shows a synodic rotation period of 12.12±0.02 h, consistent with results from other observers. We observed a maximum peak-to-peak amplitude of 1.05 magnitudes, which—depending on the taxonomic class assumed when correcting for phase angle effects—implies a minimum axial ratio of 2.14. The shape of the rotation lightcurves varies considerably between data sets due to the changing viewing geometry. The lightcurve data from this study has been included in the shape model analysis of Kaasalainen et al. (2003 Astron. Astrophys, 405, L29-L32) and the Hapke analysis of Lederer et al. (2005 Icarus 173,153-165). Color variations were also observed, with the interpolated color indices at lightcurve midpoint being: (B-V)=0.94±0.05, (V-R)=0.40±0.06, and (V-I)=0.74±0.07. Our low resolution Palomar spectra from March 2001 covered a wavelength range of 0.3-1.0 μm. We measured a spectral slope of 9.3±0.3%/100 nm between 0.55-0.70 μm and a deep 1 μm absorption (equivalent ECAS color: w-x=−0.111±0.003, v-x=0.031±0.003). Comparison of our spectra with published ECAS data from Zellner et al. (1985 Icarus 61, 355-416) indicates that this object is most likely of Q- or S-type, similar to ordinary chondrite meteorites. Our data are more consistent with a Q-type body when both the spectroscopic data and the available BVRI photometry are taken into account.     

Spin axis of (2953) Vysheslavia and its implications

Photometric observations made during the years 2000-2005 are used to determine the pole orientation of (2953) Vysheslavia, a ?15-km size member of the Koronis family. We find admissible solutions for ecliptic latitude and longitude of the rotation pole P₃: βp=−64°±10° and λp=11°±8° or P₄: βp=−68°±8° and λp=192°±8°. These imply obliquity values γ=154°±14° and γ=157°±11°, respectively. The sidereal rotation period is Psid=0.2622722±0.0000018 day. This result is interesting for two reasons: (i) the obliquity value between 90° and 180° is consistent with a prediction done by Vokrouhlický et al. [Vokrouhlický, D., Bro?, M., Farinella, P., Kne?evi?, Z., 2001. Icarus 150, 78-93] that Vysheslavia might have been transported to its unstable orbit by the Yarkovsky effect, and (ii) with the obliquity close to 180°, Vysheslavia seems to belong to one of the two distinct groups in the Koronis family found recently by Slivan [Slivan, S.M., 2002. Nature 419, 49-51], further supporting the case of dichotomy in the spin axis distribution in this family. We also argue against the possibility that Vysheslavia reached its current orbit by a recent collisional breakup.     

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.     

Multi-wavelength observations of Asteroid 2100 Ra-Shalom

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

The orbit and size distribution of small Solar System objects orbiting the Sun interior to the Earth's orbit

We present the first observational measurement of the orbit and size distribution of small Solar System objects whose orbits are wholly interior to the Earth's (Inner Earth Objects, IEOs, with aphelion <0.983 AU). We show that we are able to model the detections of near-Earth objects (NEO) by the Catalina Sky Survey (CSS) using a detailed parameterization of the CSS survey cadence and detection efficiencies as implemented within the Jedicke et al. [Jedicke, R., Morbidelli, A., Spahr, T., Petit, J.M., Bottke, W.F., 2003. Icarus 161, 17-33] survey simulator and utilizing the Bottke et al. [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] model of the NEO population's size and orbit distribution. We then show that the CSS detections of 4 IEOs are consistent with the Bottke et al. [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] IEO model. Observational selection effects for the IEOs discovered by the CSS were then determined using the survey simulator in order to calculate the corrected number and H distribution of the IEOs. The actual number of IEOs with H<18 (21) is 36±26 (530±240) and the slope of the H magnitude distribution (∝¹⁰αH) for the IEOs is . The slope is consistent with previous measurements for the NEO population of αNEO=0.35±0.02 [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] and αNEO=0.39±0.013 [Stuart, J.S., Binzel, R.P., 2004. Icarus 170, 295-311]. Based on the agreement between the predicted and observed IEO orbit and absolute magnitude distributions there is no indication of any non-gravitational effects (e.g. Yarkovsky, tidal disruption) affecting the known IEO population.     

Rotation and color properties of the nucleus of Comet 2P/Encke

We present results from CCD observations of Comet 2P/Encke acquired at Steward Observatory's 2.3 m Bok Telescope on Kitt Peak. The observations were carried out in October 2002 when the comet was near aphelion. Rotational lightcurves in B-, V-, and R-filters were acquired over two nights of observations, and analysed to study the physical and color properties of the nucleus. The average apparent R-filter magnitude across both nights corresponds to a mean effective radius of 3.95±0.06 km, and this value is similar to that found for the V- and B-filters. Taking the observed brightness range, we obtain a/b?1.44±0.06 for the semi-axial ratio of Encke's nucleus. Applying the axial ratio to the R-filter photometry gives nucleus semi-axes of [3.60±0.09]×[5.20±0.13] km, using the empirically-derived albedo and phase coefficient. No coma or tail was seen despite deep imaging of the comet, and flux limits from potential unresolved coma do not exceed a few percent of the total measured flux, for standard coma models. This is consistent with many other published data sets taken when the comet was near aphelion. Our data includes the first detailed time series multi-color measurements of a cometary nucleus, and significant color variations were seen on October 3, though not repeated on October 4. The average color indices across both nights are: (V−R)=0.39±0.06 and (B−V)=0.73±0.06 (). We analysed the R-filter time-series photometry using the method of Harris et al. [Harris, A.W., Young, J.W., Bowell, E., Martin, L.J., Millis, R.L., Poutanen, M., Scaltriti, F., Zappala, V., Schober, H.J., Debehogne, H., Zeigler, K.W., 1989. Icarus 77, 171-186] to constrain the rotation period of the comet's nucleus, and find that a period of ∼11.45 h will satisfy the data, however the errors bars are large. We have successfully linked our data with the September 2002 data from Fernández et al. [Fernández, Y.R., Lowry, S.C., Weissman, P.R., Mueller, B.E.A., Samarasinha, N.H., Belton, M.J.S., Meech, K.J., 2005. Icarus 175, 194-214]—taken just 2-3 weeks before the current data set—and we show that a rotation period of just over 11 h works extremely well for the combined data set. The resulting best-fit period is 11.083±0.003 h, consistent with the Fernández et al. value.     

Determination of physical properties of the Asteroid (41) Daphne from interferometric observations in the thermal infrared

We describe interferometric observations of the Asteroid (41) Daphne in the thermal infrared obtained with the Mid-Infrared Interferometric Instrument (MIDI) and the Auxiliary Telescopes (ATs) of the European Southern Observatory (ESO) Very Large Telescope Interferometer (VLTI). We derived the size and the surface thermal properties of (41) Daphne by means of a thermophysical model (TPM), which is used for the interpretation of interferometric data for the first time. From our TPM analysis, we derived a volume equivalent diameter for (41) Daphne of 189 km, using a non-convex 3-D shape model derived from optical lightcurves and adaptive optics images (B. Carry, private communication). On the other hand, when using the convex shape of Kaasalainen et al. (Kaasalainen, M., Mottola, S., Fulchignoni, M. [2002]. Icarus 159, 369-395) in our TPM analysis, the resulting volume equivalent diameter of (41) Daphne is between 194 and 209 km, depending on the surface roughness. The shape of the asteroid is used as an a priori information in our TPM analysis. No attempt is made to adjust the shape to the data. Only the size of the asteroid and its thermal parameters such as, albedo, thermal inertia and roughness are adjusted to the data. We estimated our model systematic uncertainty to be of 4% and of 7% on the determination of the asteroid volume equivalent diameter depending on whether the non-convex or the convex shape is used, respectively. In terms of thermal properties, we derived a value of the surface thermal inertia smaller than 50 J m⁻² s^−0.5 K⁻¹ and preferably in the range between 0 and ∼30 J m⁻² s^−0.5 K⁻¹. Our TPM analysis also shows that Daphne has a moderate macroscopic surface roughness.     

Looking a gift horse in the mouth: Evaluation of wide-field asteroid photometric surveys

It has recently become possible to do a photometric survey of many asteroids at once, rather than observing single asteroids one (or occasionally a couple) at a time. We evaluate two such surveys. Dermawan et al. (Dermawan et al. [2011]. Publ. Astron. Soc. Jpn. 63, S555–S576) observed one night on the Subaru 8.2 m telescope, and Masiero et al. (Masiero, J., Jedicke, R., Durech, J., Gwen, S., Denneau, L., Larsen, J. [2009]. Icarus 204, 145–171) observed six nights over 2 weeks with the 3.6 m CFHT. Dermawan claimed 83 rotation periods from 127 detected asteroids; Masiero et al. claimed 218 rotation periods from 828 detections. Both teams claim a number of super-fast rotators (P < 2.2 h) among main belt asteroids larger than 250 m diameter, some up to several km in diameter. This would imply that the spin rate distribution of main belt asteroids differs from like-sized NEAs, that there are larger super-fast rotators (monolithic asteroids) in the main belt than among NEAs. Here we evaluate these survey results, applying the same criteria for reliability of results that we apply to all results listed in our Lightcurve Database (Warner, B.D., Harris, A.W., Pravec, P. [2009a]. Icarus 202, 134–146). In doing so, we assigned reliability estimates judged sufficient for inclusion in statistical studies for only 27 out of 83 (33%) periods claimed by Dermawan, and only 87 out of 218 (40%) periods reported by Masiero et al.; none of the super-fast rotators larger than about 250 m diameter claimed by either survey received a reliability rating judged sufficient for analysis. We find no reliable basis for the claim of different rotation properties between main belt and near-Earth asteroids. Our analysis presents a cautionary message for future surveys.     

Cassini ISS astrometric observations of the inner jovian satellites, Amalthea and Thebe

We present a total of 289 new astrometric observations of the inner jovian satellites, Amalthea and Thebe, obtained using the Cassini ISS narrow angle camera. Observations were made using image sequences from 2000 December 11-12 (inbound) and 2001 January 15-16 (outbound), at phase angles of approximately 2° and 122°, respectively. Target distances were of order 284 R_J, giving a maximum resolution of approximately 100 km/pixel. Centroided line and sample values for 239 observations of Amalthea and 50 of Thebe are provided, together with estimated camera pointing information for each image. Orbit fitting using a uniformly precessing Keplerian ellipse model, taking into account the oblateness of Jupiter up to terms in J₆, gave RMS fit residuals of 0.364 and 0.443 pixel for Amalthea and Thebe, respectively (equivalent to 0.450 and 0.547 arcsec). RMS residuals relative to the JPL JUP230 ephemeris were 0.306 and 0.604 pixel (equivalent to 0.378 and 0.746 arcsec), for Amalthea and Thebe. The fitted orbital parameters confirm the relatively high inclinations of these satellites (0.374°±0.002° and 1.076°±0.003°, respectively), equivalent to maximum vertical displacements above Jupiter's equatorial plane of 1188±6 and 4240±12 km, respectively, consistent with current estimates of the half-thicknesses of the Amalthea and Thebe gossamer rings [Ockert-Bell, M.E., Burns, J.A., Dauber, I.J., Thomas, P.C., Veverka, J., Belton, M.J.S., Klaasen, K.P., 1999. Icarus 138, 188-213].     

Spin vectors in the Koronis family: II. Additional clustered spins, and one stray

We recorded 101 new rotation lightcurves of five Koronis family members, and then combined the new observations with previous data to determine the objects' sidereal rotation periods, spin vector orientations, and model shape solutions. The observing program was tailored specifically for spin vector analyses by determining single-apparition Lumme–Bowell solar phase coefficients, and by measuring synodic rotation periods precisely enough to unambiguously count the rotations between two consecutive oppositions, which is a prerequisite for identifying the correct sidereal period. The new data make possible first pole and shape determinations for (263) Dresda, (462) Eriphyla, and (1289) Kutaïssi, and they improve the models for (277) Elvira and (534) Nassovia, two objects previously studied by Slivan et al. [Slivan, S.M., Binzel, R.P., Crespo da Silva, L.D., Kaasalainen, M., Lyndaker, M.M., Kr?o, M., 2003. Icarus 162, 285–307]. Our results increase the number of Koronis family spin vectors reported in the literature to fourteen, a sample which now includes the seven largest family members. The spin properties of Eriphyla (rotation period , spin vector obliquity ε=51°) and Kutaïssi (P=3.62 h, ε=165°) are consistent with the markedly nonrandom distribution reported by Slivan [Slivan, S.M., 2002. Nature 419, 49–51], and explained by Vokrouhlický et al. [Vokrouhlický, D., Nesvorný, D., Bottke, W.F., 2003. Nature 425, 147–151] as the result of the effects of thermal “YORP” torques combined with solar and planetary gravitational torques. Dresda (P=16.81 h, ε=16°) is the first prograde Koronis member whose spin obliquity and spin rate significantly differ from the clustered spin properties previously found for other prograde Koronis members; nevertheless, its spin vector is consistent with several of the spin evolution possibilities that were identified in the YORP modeling.     

Bias-corrected population, size distribution, and impact hazard for the near-Earth objects

Utilizing the largest available data sets for the observed taxonomic (Binzel et al., 2004, Icarus 170, 259-294) and albedo (Delbo et al., 2003, Icarus 166, 116-130) distributions of the near-Earth object population, we model the bias-corrected population. Diameter-limited fractional abundances of the taxonomic complexes are A-0.2%; C-10%, D-17%, O-0.5%, Q-14%, R-0.1%, S-22%, U-0.4%, V-1%, X-34%. In a diameter-limited sample, ∼30% of the NEO population has jovian Tisserand parameter less than 3, where the D-types and X-types dominate. The large contribution from the X-types is surprising and highlights the need to better understand this group with more albedo measurements. Combining the C, D, and X complexes into a “dark” group and the others into a “bright” group yields a debiased dark-to-bright ratio of ∼1.6. Overall, the bias-corrected mean albedo for the NEO population is 0.14±0.02, for which an H magnitude of 17.8±0.1 translates to a diameter of 1 km, in close agreement with Morbidelli et al. (2002, Icarus 158 (2), 329-342). Coupling this bias corrected taxonomic and albedo model with the H magnitude dependent size distribution of (Stuart, 2001, Science 294, 1691-1693) yields a diameter distribution with 1090±180 NEOs with diameters larger than 1 km. As of 2004 June, the Spaceguard Survey has discovered 56% of the NEOs larger than 1 km. Using our size distribution model, and orbital distribution of (Stuart, 2001, Science 294, 1691-1693) we calculate the frequency of impacts into the Earth and the Moon. Globally destructive collisions (∼10²¹ J) of asteroids 1 km or larger strike the Earth once every 0.60±0.1 Myr on average. Regionally destructive collisions with impact energy greater than 4×10¹⁸ J (∼200 m diameter) strike the Earth every 56,000±6000 yr. Collisions in the range of the Tunguska event (4-8×10¹⁶ J) occur every 2000-3000 yr. These values represent the average time between randomly spaced impacts; actual impacts could occur more or less closely spaced solely by chance. As a verification of these impact rates, the crater production function of Shoemaker et al. (1990, Geological Society of American Special Paper 247) has been updated by combining this new population model with a crater formation model to find that the observed crater production function on both the Earth and Moon agrees with the rate of crater production expected from the current population of NEOs.     

A comparison of detection sensitivity between ALTAIR and Arecibo meteor observations: Can high power and large aperture radars detect low velocity meteor head-echoes

Meteor head-echo observations using High Power and Large Aperture (HPLA) radars have been routinely used for micrometeor studies for over a decade. The head-echo is a signal from the radar-reflective plasma region traveling with the meteoroid and its detection allows for very precise determination of instantaneous meteor altitude, velocity and deceleration. Unlike specular meteor radars (SMR), HPLA radars are diverse instruments when compared one to another. The operating frequencies range from 46 MHz to 1.29 GHz while the antenna configurations changes from 18,000 dipoles in a 300 m×300 m square array, phase arrays of dipoles to single spherical or parabolic dishes of various dimensions. Hunt et al. [Hunt, S.M., Oppenheim, M., Close, S., Brown, P.G., McKeen, F., Minardi, M., 2004. Icarus 168, 34-42] and Close et al. [Close, S., Brown, P., Campbell-Brown, M., Oppenheim, M., Colestock, P., 2007. Icarus, doi:10.1016/j.icarus.2006.09.07] recently showed, by utilizing a head-echo plasma-based model, the presence of instrumental biases in the ALTAIR VHF radar system against detecting meteors produced by very small particles (<1 μg) moving at slow (∼20 km/s) velocities due to the low head echo radar cross-section (RCS) associated with these particles. In this paper we apply the same methodology to the Arecibo 430 MHz radar and compare the results with those presented by Close et al. [Close, S., Brown, P., Campbell-Brown, M., Oppenheim, M., Colestock, P., 2007. Icarus, doi:10.1016/j.icarus.2006.09.07]. We show that, if the methodology applied by Hunt et al. [Hunt, S.M., Oppenheim, M., Close, S., Brown, P.G., McKeen, F., Minardi, M., 2004. Icarus 168, 34-42] and Close et al. [Close, S., Brown, P., Campbell-Brown, M., Oppenheim, M., Colestock, P., 2007. Icarus, doi:10.1016/j.icarus.2006.09.07] is accurate, for particles at least 1 μg or heavier, while the bias may exist for the ALTAIR measurements, it does not exist in the Arecibo data due to its greater sensitivity.     

Radar observations of Asteroids 64 Angelina and 69 Hesperia

We report new radar observations of E-class Asteroid 64 Angelina and M-class Asteroid 69 Hesperia obtained with the Arecibo Observatory S-band radar (2480 MHz, 12.6 cm). Our measurements of Angelina’s radar bandwidth are consistent with reported diameters and poles. We find Angelina’s circular polarization ratio to be 0.8 ± 0.1, tied with 434 Hungaria for the highest value observed for main-belt asteroids and consistent with the high values observed for all E-class asteroids (Benner, L.A.M., Ostro, S.J., Magri, C., Nolan, M.C., Howell, E.S., Giorgini, J.D., Jurgens, R.F., Margot, J.L., Taylor, P.A., Busch, M.W., Shepard, M.K. [2008]. Icarus 198, 294-304; Shepard, M.K., Kressler, K.M., Clark, B.E., Ockert-Bell, M.E., Nolan, M.C., Howell, E.S., Magri, C., Giorgini, J.D., Benner, L.A.M., Ostro, S.J. [2008b]. Icarus 195, 220-225). Our radar observations of 69 Hesperia, combined with lightcurve-based shape models, lead to a diameter estimate, D_eff = 110 ± 15 km, approximately 20% smaller than the reported IRAS value. We estimate Hesperia to have a radar albedo of , consistent with a high-metal content. We therefore add 69 Hesperia to the Mm-class (high metal M) (Shepard, M.K., Clark, B.E., Ockert-Bell, M., Nolan, M.C., Howell, E.S., Magri, C., Giorgini, J.D., Benner, L.A.M., Ostro, S.J., Harris, A.W., Warner, B.D., Stephens, R.D., Mueller, M. [2010]. Icarus 208, 221-237), bringing the total number of Mm-class objects to eight; this is 40% of all M-class asteroids observed by radar to date.     

Study of the vertical structure of Saturn's atmosphere using HST/WFPC2 images

We have studied the vertical structure of hazes at six different latitudes (−60°, −50°, −30°, −10°, +30°, and +50°) on Saturn's atmosphere. For that purpose we have compared the results of our forward radiative transfer model to limb-to-limb reflectivity scans at four different wavelengths (230, 275, 673.2, and 893 nm). The images were obtained with the Hubble Space Telescope Wide Field Planetary Camera 2 in September 1997, during fall on Saturn's northern hemisphere. The spatial distribution of particles appears to be very variable with latitude both in the stratosphere and troposphere. For the latitude range +50° to −50°, an atmospheric structure consisting of a stratospheric haze and a tropospheric haze interspersed by clear gas regions has been found adequate to explain the center to limb reflectivities at the different wavelengths. This atmospheric structure has been previously used by Ortiz et al. (1996, Icarus 119, 53-66) and Stam et al. (2001, Icarus 152, 407-422). In this work the top of the tropospheric haze is found to be higher at the southern latitudes than at northern latitudes. This hemispherical asymmetry seems to be related to seasonal effects. Different latitudes experience different amount of solar insolation that can affect the atmospheric structure as the season varies with time. The haze optical thickness is largest (about 30 at 673.2 nm) at latitudes ±50 and −10 degrees, and smallest (about 18) at ±30 degrees. The stratospheric haze is found to be optically thin at all studied latitudes from −50 to +50 degrees being maximum at −10° (τ=0.033). At −60° latitude, where the UV images show a strong darkening compared to other regions on the planet, the cloud structure is remarkably different when compared to the other latitudes. Here, aerosol and gas are found to be uniformly mixed down to the 400 mbar level.     

Simulation of Titan haze formation using a photochemical flow reactor: The optical constants of the polymer

Solar UV is the principal energy source impinging the atmosphere of Titan while the energy from the electrons in Saturn's magnetosphere is less than 0.5% of the UV light. Titan haze analogs were prepared by the photolysis of a mixture of gases that simulate the composition of its atmosphere (nitrogen, methane, hydrogen, acetylene, ethylene, and cyanoacetylene). The real (n) and imaginary (k) parts of the complex refractive index of haze analogs formed from four different gas mixtures were calculated from the spectral properties of the solid polymer in UV-visible, near infrared and infrared wavelength spectral regions. The value of n was constant at 1.6±0.1 throughout the 0.2-2.5 μm region. The variation of k with wavelength for the values derived for Titan has a lower error than the absolute values of k so the more significant comparisons are with the slopes of the k(λ) plots in the UV-VIS region. Three of the photochemical Titan haze analogs had slopes comparable to those derived for Titan from the Voyager data (Rages and Pollack, 1980, Icarus 41, 119-130; McKay and Toon, 1992, in: Proceedings of the Symposium on Titan, in: ESA SP, Vol. 338, pp. 185-190). The slopes of the k(λ) plots for haze analogs prepared by spark discharge (Khare et al., 1984, Icarus 60, 127-137) and plasma discharge (Ramirez et al., 2002, Icarus 156, 515-529) were also comparable to Titan's. These finding show that the k(λ) plots do not differentiate between different laboratory simulations of atmospheric chemistry on Titan in the UV-VIS near IR region (0.2-2.5 microns). There is a large difference between the k(λ) in the infrared between the haze analogs prepared photochemically and analogs prepared using a plasma discharges (Khare et al., 1984, Icarus 60, 127-137; Coll et al., 1999, Planet. Space Sci. 47, 1331-1340; Khare et al., 2002, Icarus 160, 172-182). The C/N ratio in the haze analog prepared by discharges is in the 2-11 range while that of the photochemical analogs is in the 18-24 range. The use of discharges and UV light for initiating the chemistry in Titan's atmosphere is discussed.     

Scaling of melt production in hypervelocity impacts from high-resolution numerical simulations

The volume of melt produced in hypervelocity planetary impacts and the size and shape of the melted region are key to understanding the impact histories of solid planetary bodies and the geological effects of impacts on their surfaces and interiors. Prior work of Pierazzo et al. (Pierazzo, E., Vickery, A.M., Melosh, H.J. [1997]. Icarus 127, 408-423) gave the first estimates of impact melt production in geological materials using a modern hydrocode and equation of state. However, computational limits at the time forced use of low resolution, which may have resulted in low melt volumes. Our simulations with 50 times higher resolution provide independent confirmation of the Pierazzo et al. (Pierazzo, E., Vickery, A.M., Melosh, H.J. [1997]. Icarus 127, 408-423) melt volumes in aluminum, iron, dunite, and granite impacts at velocities between 20 and 80 km/s. In ice/ice impacts, we find that melt volumes depend on target temperature and are lower than predicted by Pierazzo et al. (Pierazzo, E., Vickery, A.M., Melosh, H.J. [1997]. Icarus 127, 408-423). Our melt volumes are directly proportional to impact energy for all materials, over a wide range of impact velocity. We also report new data for melt volume scalings for ice/dunite and iron/dunite impacts and the size and shape of melted region, valuable for interpretation of cratering records and studies of impact-induced differentiation.     

Near-infrared spectral observations and interpretations for S-asteroids 138 Tolosa, 306 Unitas, 346 Hermentaria, and 480 Hansa

Near-infrared (∼0.7 to ∼2.5 μm) spectra of S-asteroids 138 Tolosa, 306 Unitas, 346 Hermentaria, and 480 Hansa suggest the presence of variable amounts of orthopyroxene ± clinopyroxene ± olivine ± plagioclase feldspar on the surfaces of these asteroids. The spectra of these asteroids were compared to laboratory mineral mixtures of orthopyroxene, clinopyroxene, and olivine [Singer, R.B., 1981. J. Geophys. Res. 86 (B9), 7967-7982; Cloutis, E.A., 1985. Master's thesis]. The band parameters (band centers, band areas) were quantified and temperature-corrected [Moroz et al., 2000. Icarus 147, 79-93; Gaffey et al., 2002. In: Bottke Jr., W.F., Cellino, A., Paolicchi, P., Binzel, R.P. (Eds.), Asteroids III. The University of Arizona Press, Tucson, pp. 183-204]. Each S-asteroid in this paper exhibits an overall spectral shape with band parameters that are inconsistent with ordinary chondrite near-infrared spectra and their inferred mineral abundances and/or pyroxene chemistries. 138 Tolosa displays a complex spectrum with a broad ∼1 μm absorption feature that displays a double Band I minimum, a well-defined absorption at ∼1.3 μm, and a broad, but weak absorption in the ∼2 μm region. Although different interpretations exist, the Tolosa spectrum is most consistent with a ∼60/40 mixture of Type B clinopyroxene and orthopyroxene. Spectra of 306 Unitas suggest a surface with variable amounts of low-Ca pyroxene and olivine. Unitas is located in the S-(IV) and S-(VI) subtype regions in Gaffey et al. [1993. Icarus 106, 573-602]. 346 Hermentaria exhibits a complex, broad Band I absorption feature and a weak Band II feature, which suggests a ∼50/50 mixture of clinopyroxene and orthopyroxene. Hermentaria is classified as an S-(III). Spectra of 480 Hansa suggest a dominant low-Ca pyroxene component with lesser amounts of olivine. Based on these characterizations, these four S-asteroids should not be considered as potential ordinary chondrite parent bodies. Furthermore, these results suggest that these S-asteroids experienced at least partial melting temperatures [T?∼950 °C: Gaffey et al., 1993. Icarus 106, 573-602; Keil, K., 2000. Planet. Space Sci. 48, 887-903] during the formation epoch in the early Solar System. Continuing spectroscopic investigations will discern the relative abundance of chondritic and thermally-evolved objects among the S-type asteroids that have survived since the formation epoch ∼4.56 billion years ago.     

Linking the collisional history of the main asteroid belt to its dynamical excitation and depletion

The main belt is believed to have originally contained an Earth mass or more of material, enough to allow the asteroids to accrete on relatively short timescales. The present-day main belt, however, only contains ∼5×¹⁰⁻⁴ Earth masses. Numerical simulations suggest that this mass loss can be explained by the dynamical depletion of main belt material via gravitational perturbations from planetary embryos and a newly-formed Jupiter. To explore this scenario, we combined dynamical results from Petit et al. [Petit, J. Morbidelli, A., Chambers, J., 2001. The primordial excitation and clearing of the asteroid belt. Icarus 153, 338-347] with a collisional evolution code capable of tracking how the main belt undergoes comminution and dynamical depletion over 4.6 Gyr [Bottke, W.F., Durda, D., Nesvorny, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H., 2005. The fossilized size distribution of the main asteroid belt. Icarus 175, 111-140]. Our results were constrained by the main belt's size-frequency distribution, the number of asteroid families produced by disruption events from diameter D>100 km parent bodies over the last 3-4 Gyr, the presence of a single large impact crater on Vesta's intact basaltic crust, and the relatively constant lunar and terrestrial impactor flux over the last 3 Gyr. We used our model to set limits on the initial size of the main belt as well as Jupiter's formation time. We find the most likely formation time for Jupiter was 3.3±2.6 Myr after the onset of fragmentation in the main belt. These results are consistent with the estimated mean disk lifetime of 3 Myr predicted by Haisch et al. [Haisch, K.E., Lada, E.A., Lada, C.J., 2001. Disk frequencies and lifetimes in young clusters. Astrophys. J. 553, L153-L156]. The post-accretion main belt population, in the form of diameter D?1000 km planetesimals, was likely to have been 160±40 times the current main belt's mass. This corresponds to 0.06-0.1 Earth masses, only a small fraction of the total mass thought to have existed in the main belt zone during planet formation. The remaining mass was most likely taken up by planetary embryos formed in the same region. Our results suggest that numerous D>200 km planetesimals disrupted early in Solar System history, but only a small fraction of their fragments survived the dynamical depletion event described above. We believe this may explain the limited presence of iron-rich M-type, olivine-rich A-type, and non-Vesta V-type asteroids in the main belt today. The collisional lifetimes determined for main belt asteroids agree with the cosmic ray exposure ages of stony meteorites and are consistent with the limited collisional evolution detected among large Koronis family members. Using the same model, we investigated the near-Earth object (NEO) population. We show the shape of the NEO size distribution is a reflection of the main belt population, with main belt asteroids driven to resonances by Yarkovsky thermal forces. We used our model of the NEO population over the last 3 Gyr, which is consistent with the current population determined by telescopic and satellite data, to explore whether the majority of small craters (D<0.1-1 km) formed on Mercury, the Moon, and Mars were produced by primary impacts or by secondary impacts generated by ejecta from large craters. Our results suggest that most small craters formed on these worlds were a by-product of secondary rather than primary impacts.     

Near-nucleus photometry of comets using archived NEAT data

Though optimized to discover and track fast moving Near-Earth Objects (NEOs), the Near-Earth Asteroid Tracking (NEAT) survey dataset can be mined to obtain information on the comet population observed serendipitously during the asteroid survey. We have completed analysis of over 400 CCD images of comets obtained during the autonomous operations of two 1.2-m telescopes: the first on the summit of Haleakala on the Hawaiian island of Maui and the second on Palomar Mountain in southern California. Photometric calibrations of each frame were derived using background catalog stars and the near-nucleus comet photometry measured. We measured dust production and normalized magnitudes for the coma and nucleus in order to explore cometary activity and comet size-frequency distributions. Our data over an approximately two-year time frame (2001 August-2003 February) include 52 comets: 12 periodic, 19 numbered, and 21 non-periodic, obtained over a wide range of viewing geometries and helio/geocentric distances. Nuclear magnitudes were estimated for a subset of comets observed. We found that for low-activity comets (Afρ<100 cm) our model gave reasonable estimates for nuclear size and magnitude. The slope of the cumulative luminosity function of our sample of low-activity comets was 0.33 ± 0.04, consistent with the slope we measured for the Jupiter-family cometary nuclei collected by Fernández et al. [Fernández, J.A., Tancredi, G., Rickman, H., Licandro, J., 1999. Astron. Astrophys. 392, 327-340] of 0.38 ± 0.02. Our slopes of the cumulative size distribution α=1.50±0.08 agree well with the slopes measured by Whitman et al. [Whitman, K., Morbidelli, A., Jedicke, R., 2006. Icarus 183, 101-114], Meech et al. [Meech, K.J., Hainaut, O.R., Marsden, B.G., 2004. Icarus 170, 463-491], Lowry et al. [Lowry, S.C., Fitzsimmons, A., Collander-Brown, S., 2003. Astron. Astrophys. 397, 329-343], and Weissman and Lowry [Weissman, P.R., Lowry, S.C., 2003. Lunar Planet. Sci. 34. Abstract 34].