HST BVI photometry of Triton and Proteus

BVI photometry of Triton and Proteus was derived from HST images taken in 1997. The VEGAMAG photometric technique was used. Triton was found to be brighter by a few percent than observations of the 1970's and 1980's, as expected due to the increasingly greater exposure of the bright south polar region. The leading side was also found to be brighter than the trailing side by 0.09 mag in all filters—50% larger than reported by Franz [Franz, O.G., 1981. Icarus 45, 602-606]. Contrary to our previous results [Pascu, D., et al., 1998. Bull. Am. Astron. Soc. 30, 1101], we found no episodic reddening. Our previous conclusions were based on an inaccurate early version of the Charge Transfer Efficiency (CTE) correction. The present result limits the start of the reddening event reported by Hicks and Buratti [Hicks, M.D., Buratti, B.J., 2004. Icarus 171, 210-218]. Our (B-V) result of 0.70±0.01 supports the global blueing described by Buratti et al. [Buratti, B.J., Goguen, J.D., Gibson, J., Mosher, J., 1994. Icarus 110, 303-314]. Our observations of July 1997 agree with the Voyager results and are among the bluest colors seen. We found Proteus somewhat brighter than earlier studies, but in good agreement with the recent value given by Karkoschka [Karkoschka, E., 2003. Icarus 162, 400-407]. A leading/trailing brightness asymmetry was detected for Proteus, with the leading side 0.1 mag brighter. The unique differences in action of the endogenic and exogenic processes on Triton and Proteus provides an opportunity to separate the endogenic and exogenic effects on Triton.     

The roughness of the dark side of Iapetus from the 2004 to 2005 flyby

The roughness of a planetary surface offers clues to its past geologic history. We apply a surface roughness model developed by Buratti and Veverka (Buratti, B.J., Veverka, J. [1985]. Icarus 64, 320-328) to Cassini ISS data from the January 1st, 2005 flyby of Iapetus. This model uses the observed scattering behavior to provide a depth to radius factor q quantifying the size of idealized craters on the surface. Our findings indicate that the surface on the dark side is significantly smoother than the surfaces of other icy low-albedo saturnian satellites. We have found that the average depth to radius on the leading (dark) side is 0.084, corresponding to a Hapke mean slope angle of 6°. As compared to the 13-33° Hapke mean slope angle of other icy satellites (Buratti, B.J., and 10 colleagues [2008]. Icarus 193, 309-322), our results present a clearly different picture for the leading surface of Iapetus, suggesting that the dark deposit contributes to the decrease in macroscopic surface roughness of the leading side. Attempts were made to obtain an average depth to radius value for the trailing (bright) side; however the scans of the bright side from this flyby exhibited large variations in albedo, resulting in results that were physically unrealistic.     

Dynamic features of Io's extended sodium distributions

We report on two-dimensional imaging observations of D-line emissions from the extended distribution of iogenic sodium atoms with two fields of view (±20 RJ (narrow FOV) and ±400 RJ (wide FOV)) simultaneously by using a portable small telescope or camera lens. We derived dynamic feature of the band-shaped and spray-shaped distributions near Io's orbit by means of continuous observation. The observations confirm the phenomenological behavior of the sodium cloud on two spatial scales, as previously observed by Pilcher et al. [Pilcher, C.B., Smyth, W.H., Combi, M.R., Fertel, J.H., 1984. Astrophys. J. 287, 427-444], Schneider et al. [Schneider, N.M., Trauger, J.T., Wilson, J.K., Brown, D.I., Evans, R.W., Shemansky, D.E., 1991. Science 253, 1394-1397], and Mendillo et al. [Mendillo, M., Baumgartner, J., Flynn, B., Hughes, W.S., 1990. Nature 348, 312-314]. We also confirm an elongated oval emission distribution of the sodium nebula and derivation of its detailed east-west asymmetry depending on Io's phase angle, which was first noted by Flynn et al. [Flynn, B., Mendillo, M., Baumgartner, J., 1994. J. Geophys. Res. 99, 8403-8409]. We then did model analyses to investigate the source process for sodium atoms and the dynamics behind their distribution. We conclude that the essential of molecular ion mechanisms to the band-shaped distribution is in agreement with Wilson and Schneider [Wilson, J.K., Schneider, N.M., 1999. J. Geophys. Res. 104, 16567-16583]. We differ from Wilson et al. [Wilson, J.K., Mendillo, M., Baumgartner, J., Schneider, N.M., Trauger, J.T., Flynn, B., 2002. Icarus 157, 476-489] in finding that charge exchange process contributes more to the spray-shaped distribution and sodium nebula than sputtering does. These results derived the double-peaked velocity distribution of released sodium atoms, and re-confirmed the source rates in agreement with past studies.     

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.     

Stratigraphical and morphological evidence for pingo genesis in the Cerberus plains

‘Rootless’ debris cones (or pseudocraters) occur in platy, patterned ground throughout the Cerberus plains of Mars and are thought to represent the products of explosive magma-ice interaction [Lanagan et al., 2001. Geophys. Res. Lett. 28, 2365-2368; Fagents et al., 2002. In: Smellie, J.L., Chapman, M.G. (Eds.), Volcano-Ice Interaction on Earth and Mars. In: Geol. Soc. Spec. Publ., vol. 202, pp. 295-317]. Requiring lava and water interspersed, they are central to theories of multiple magmatic and aqueous flood events [Burr et al., 2002. Icarus 159, 53-73; Berman, D.C., Hartmann, W.K., 2002. Icarus 159, 1-17] and widespread sheet volcanism [Keszthelyi et al., 2000. J. Geophys. Res. 105, 15027-15049] in the region during the late Amazonian (a region reported to have been occupied by water bodies ranging from lakes to oceans [Scott et al., 1995. Map of Mars showing channels and possible paleolake basins. USGS Miscellaneous Investigations Series, Map I-2461 (1:30,000,000)]). The nature of the platy substrate is the subject of debate, with evidence given for lava [Keszthelyi et al., 2000. J. Geophys. Res. 105, 15027-15049; Plescia, J.B., 2003. Icarus 164, 79-95] and ice [Brakenridge, G.R., 1993. Lunar Planet. Sci. XXIV (Part 1), 175-176; Rice et al., 2002. Lunar Planet. Sci. XXXIII. Abstract #2026; Murray et al., 2005. Nature 434, 352-355]. The superposition relationships of cones and platy deposits in the channels of the Athabasca Valles precludes a magmatic origin, indicating later formation as permafrost mounds (or ‘pingos’), with implications for geologically recent flood volcanism, age constraints on young surfaces and recent climate change on Mars.     

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.     

Simulations of planet migration driven by planetesimal scattering

Evidence has mounted for some time that planet migration is an important part of the formation of planetary systems, both in the Solar System [Malhotra, R., 1993. Nature 365, 819-821] and in extrasolar systems [Mayor, M., Queloz, D., 1995. Nature 378, 355-359; Lin, D.N.C., Bodenheimer, P., Richardson, D.C., 1996. Nature 380, 606-607]. One mechanism that produces migration (the change in a planet's semi-major axis a over time) is the scattering of comet- and asteroid-size bodies called planetesimals [Fernandez, J.A., Ip, W.-H., 1984. Icarus 58, 109-120]. Significant angular momentum exchange can occur between the planets and the planetesimals during local scattering, enough to cause a rapid, self-sustained migration of the planet [Ida, S., Bryden, G., Lin, D.N.C., Tanaka, H., 2000. Astrophys. J. 534, 428-445]. This migration has been studied for the particular case of the four outer planets of the Solar System (as in Gomes et al. [Gomes, R.S., Morbidelli, A., Levison, H.F., 2004. Icarus 170, 492-507]), but is not well understood in general. We have used the Miranda [McNeil, D., Duncan, M., Levison, H.F., 2005. Astron. J. 130, 2884-2899] computer simulation code to perform a broad parameter-space survey of the physical variables that determine the migration of a single planet in a planetesimal disk. Migration is found to be predominantly inwards, and the migration rate is found to be independent of planet mass for low-mass planets in relatively high-mass disks. Indeed, a simple scaling relation from Ida et al. [Ida, S., Bryden, G., Lin, D.N.C., Tanaka, H., 2000. Astrophys. J. 534, 428-445] matches well with the dependencies of the migration rate:

(1)     

Compositional heterogeneity of Asteroid 4 Vesta’s southern hemisphere: Implications for the Dawn mission

High signal-to-noise, rotationally-resolved spectra of Asteroid 4 Vesta’s southern hemisphere from the 2007 opposition were used to constrain its compositional and mineralogical variations. The spectra were rotationally-phased using closely timed HST observations of Vesta by Li et al. (Li, J.-Y., McFadden, L.A., Thomas, P.C., Mutchler, M.J., Parker, J.Wm., Young, E.F., Russell, C.T., Sykes, M.V., Schmidt, B.E. [2010]. Icarus 208, 238–251). The average surface of Vesta’s southern hemisphere is analogous to a howardite or polymict eucrite assemblage similar to the northern hemisphere, although the band parameters are distinctly shifted towards the diogenite zone on the Band–Band plot. A few distinct compositional units were detected and they might be related to albedo features detected by Hubble Space Telescope (Li et al., 2010). We have identified two compositionally distinct regions overlaying the background surface. The first unit is a polymict eucrite and/or low-Ca eucrite compositional unit at 143° longitude that border the eucrite zone on the Band–Band plot and the second is a diogenite unit at 159°. While we did not detect any distinct olivine units as suggested by Gaffey (Gaffey, M.J. [1997]. Icarus 127, 130–157), we cannot rule out the possibility of smaller olivine-rich units that are below the detection limit of the instrumentation we used. Based on the analysis and the limitations of the data, we do not suggest that Vesta’s surface is olivine-free. Mean pyroxene chemistry estimates for both hemispheres broadly agree with one another (to within one-sigma) with the northern hemisphere ferrosilite (Fs) and wollastonite (Wo) values being slightly higher than southern hemisphere.     

The formation of Uranus and Neptune in solid-rich feeding zones: Connecting chemistry and dynamics

The core accretion theory of planet formation has at least two fundamental problems explaining the origins of Uranus and Neptune: (1) dynamical times in the trans-saturnian solar nebula are so long that core growth can take >15 Myr and (2) the onset of runaway gas accretion that begins when cores reach ∼10M_⊕ necessitates a sudden gas accretion cutoff just as Uranus and Neptune’s cores reach critical mass. Both problems may be resolved by allowing the ice giants to migrate outward after their formation in solid-rich feeding zones with planetesimal surface densities well above the minimum-mass solar nebula. We present new simulations of the formation of Uranus and Neptune in the solid-rich disk of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) using the initial semimajor axis distribution of the Nice model (Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A. [2005]. Nature 435, 466-469; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R. [2005]. Nature 435, 462-465; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F. [2005]. Nature 435, 459-461), with one ice giant forming at 12 AU and the other at 15 AU. The innermost ice giant reaches its present mass after 3.8-4.0 Myr and the outermost after 5.3-6 Myr, a considerable time decrease from previous one-dimensional simulations (e.g. Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y. [1996]. Icarus 124, 62-85). The core masses stay subcritical, eliminating the need for a sudden gas accretion cutoff.Our calculated carbon mass fractions of 22% are in excellent agreement with the ice giant interior models of Podolak et al. (Podolak, M., Weizman, A., Marley, M. [1995]. Planet. Space Sci. 43, 1517-1522) and Marley et al. (Marley, M.S., Gómez, P., Podolak, M. [1995]. J. Geophys. Res. 100, 23349-23354). Based on the requirement that the ice giant-forming planetesimals contain >10% mass fractions of methane ice, we can reject any Solar System formation model that initially places Uranus and Neptune inside of Saturn’s orbit. We also demonstrate that a large population of planetesimals must be present in both ice giant feeding zones throughout the lifetime of the gaseous nebula. This research marks a substantial step forward in connecting both the dynamical and chemical aspects of planet formation. Although we cannot say that the solid-rich solar nebula model of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) gives exactly the appropriate initial conditions for planet formation, rigorous chemical and dynamical tests have at least revealed it to be a viable model of the early Solar System.     

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.     

Geomorphologic mapping of the Menrva region of Titan using Cassini RADAR data

We made a detailed geomorphologic map of the Menrva region of Titan, using Cassini RADAR data as our map base. Using similar techniques and approaches that were applied to mapping Magellan radar images of Venus, and earlier, more generalized Titan maps, we were able to define and characterize 10 radar morphologic units, along with inferred dunes and fluvial channels, from the RADAR data. Structural features, such as scarps, ridges, and lineaments were also identified. Using principles of superposition, cross-cutting, and embayment relations we created a sequence of map units for this region. We interpret Menrva to be a 440 km wide degraded impact basin, in agreement with earlier studies by Elachi et al. (Elachi, C. et al. [2006]. Nature 441, 709-713) and Wood et al. (Wood, C.A., Lorenz, R., Kirk, R., Lopes, R., Mitchell, K., Stofan, E., and the Cassini RADAR Team [2010]. Icarus 206, 334-344), and identify it as the oldest feature in the map region. Exogenic processes including hydrocarbon fluid channelization forming the Elivagar Flumina channel network and dune fields resulting from aeolian activity are the current geologic processes dominating our map area, and these processes have contributed to the erosion of the crater’s ejecta field. There is evidence of multiple episodes of channel formation, erosion and burial by aeolian deposits, as observed elsewhere on Titan by e.g., Barnes et al. (Barnes, J.W. et al. [2005]. Icarus 195, 400-414). Channel outflow regions have morphologies suggestive of streams formed by flash floods, and dune fields are small and restricted rather than forming large dune seas, consistent with a desert-like environment for this region with low supply of hydrocarbon particles, also consistent with other studies by e.g., Lorenz et al. (Lorenz, R.D. et al. [2008a]. Planet. Space Sci. 56, 1132-1144). There is no evidence of cryovolcanism or non-impact-related tectonic activity in the Menrva region, although this region is too small to infer anything about the roles of these processes elsewhere on Titan. This work suggests detailed geomorphologic mapping can confidently be applied to Cassini RADAR data, and we suggest that more extensive mapping should be done using RADAR, ISS, and VIMS data geographically distributed across Titan to assess its usefulness for a future combined RADAR-ISS-VIMS-based global geologic map.     

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.     

Asteroid 21 Lutetia at 3 μm: Observations with IRTF SpeX

We present observations of Asteroid 21 Lutetia collected 2003–2008 using the SpeX instrument on the NASA Infrared Telescope Facility (IRTF) covering 2–4 μm. We also reevaluate NSFCam observations obtained in 1996 (Rivkin, A.S., Lebofsky, L.A., Clark, B.E., Howell, E.S., Britt, D.T. [2000]. Icarus 145, 351–368). Taken together, these show deeper 3-μm band depths (of order 3–5%) in the southern hemisphere of Lutetia, and shallower band depths (of order 2% or less) in the north. Such variation is consistent with observations at shorter wavelength by previous workers (Nedelcu, D.A. et al. [2007]. Astron. Astrophys. 470, 1157–1164; Lazzarin, M. et al. [2010]. Mon. Not. R. Astron. Soc. 408, 1433–1437), who observed hemispheric-level variations from C-like spectra to X-like spectra.While the shallowness of absorption bands on Lutetia hinders identification of its surface composition, goethite appears plausible as a constituent in its southern hemisphere (Beck, P., Quirico, E., Sevestre, D., Montes-Hernandez, G., Pommerol, A., Schmitt, B. [2011]. Astron. Astrophys. 526, A85–A89). Mathematical models of space weathered goethite are most consistent with Lutetia’s southern hemisphere spectrum, but more work and further observations are necessary to confirm this suggestion.     

Methane on Uranus: The case for a compact CH4 cloud layer at low latitudes and a severe CH4 depletion at high-latitudes based on re-analysis of Voyager occultation measurements and STIS spectroscopy

Lindal et al. (Lindal, G.F., Lyons, J.R., Sweetnam, D.N., Eshleman, V.R., Hinson, D.P. [1987]. J. Geophys. Res. 92 (11), 14987-15001) presented a range of temperature and methane profiles for Uranus that were consistent with 1986 Voyager radio occultation measurements of refractivity versus altitude. A localized refractivity slope variation near 1.2 bars was interpreted to be the result of a condensed methane cloud layer. However, models fit to near-IR spectra found particle concentrations much deeper in the atmosphere, in the 1.5-3 bar range (Sromovsky, L.A., Irwin, P.G.J., Fry, P.M. [2006]. Icarus 182, 577-593; Sromovsky, L.A., Fry, P.M. [2010]. Icarus 210, 211-229; Irwin, P.G.J., Teanby, N.A., Davis, G.R. [2010]. Icarus 208, 913-926), and a recent analysis of STIS spectra argued for a model in which aerosol particles formed diffusely distributed hazes, with no compact condensation layer (Karkoschka, E., Tomasko, M. [2009]. Icarus 202, 287-309). To try to reconcile these results, we reanalyzed the occultation observations with the He volume mixing ratio reduced from 0.15 to 0.116, which is near the edge of the 0.033 uncertainty range given by Conrath et al. (Conrath, B., Hanel, R., Gautier, D., Marten, A., Lindal, G. [1987]. J. Geophys. Res. 92 (11), 15003-15010). This allowed us to obtain saturated mixing ratios within the putative cloud layer and to reach above-cloud and deep methane mixing ratios compatible with STIS spectral constraints. Using a 5-layer vertical aerosol model with two compact cloud layers in the 1-3 bar region, we find that the best fit pressure for the upper compact layer is virtually identical to the pressure range inferred from the occultation analysis for a methane mixing ratio near 4% at 5°S. This strongly argues that Uranus does indeed have a compact methane cloud layer. In addition, our cloud model can fit the latitudinal variations in spectra between 30°S and 20°N, using the same profiles of temperature and methane mixing ratio. But closer to the pole, the model fails to provide accurate fits without introducing an increasingly strong upper tropospheric depletion of methane at increased latitudes, in rough agreement with the trend identified by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M. [2009]. Icarus 202, 287-309).     

Comment on “Constraints on the source of lunar cataclysm impactors” (Cuk et al., 2010, Icarus 207, 590-594)

Cuk et al. (Cuk, M., Gladman, B.J., Stewart, S.T. [2010]. Icarus 207, 590-594) argue that the projectiles bombarding the Moon at the time of the so-called lunar cataclysm could not have been mainbelt asteroids ejected by purely gravitational means, in contradiction with a conclusion that was reached by Strom et al. (Strom, R.G., Malhotra, R., Ito, T., Yoshida, F., Kring, D.A. [2005]. Science 309, 1847-1850). We demonstrate that Cuk et al.’s argument is erroneous because, contrary to their arguments, the lunar highlands do register the cataclysm impacts, lunar class 1 craters do not represent the size distribution of the cataclysm craters, and the crater size distributions on the late-forming basins are quite similar to those of the highlands craters, albeit at a lower number density due to the rapid decline of the impact flux during the cataclysm.     

Do young martian ray craters have ages consistent with the crater count system?

McEwen et al. (McEwen, A.S., Preblich, B.S., Turtle, E.P., Artemieva, N.A., Golombek, M.P., Hurst, M., Kirk, R.L., Burr, D.M., Christensen, P. [2005]. Icarus 176, 351-381) developed a useful test for the internal consistency of crater-count chronometry systems. They argued that certain multi-kilometer, fresh-looking martian craters with prominent rays should be the youngest or near-youngest craters in their size range. The “McEwen et al. test” is that the ages determined from crater densities of the smallest superimposed craters (typically diameter D ∼ 5-20 m) should thus be comparable to the expected formation intervals of the host primary. McEwen et al. concluded from MOC data that crater chronometry failed this test by factors of 700-2000. We apply HiRISE and other imagery to eight different young craters in order to re-evaluate their arguments. We use existing crater chronology systems as well as the reported observed production rate of 16 m craters (Malin, M.C., Edgett, K., Posiolova, L., McColley, S., Noe Dobrea, E. [2006]. Science 314, 1573-1557; Hartmann, W.K., Quantin, C., Mangold, N. [2007]. Icarus 186, 11-23; Kreslavsky [2007]. Seventh International Conference on Mars, 3325). Every case passes the McEwen et al. test. We conclude that the huge inconsistencies suggested by McEwen et al. are spurious. Many of these craters show evidence of impact into ice-rich material, and appear to have ice-flow features and sublimation pits on their floors. As production rate data improve, decameter-scale craters will provide a valuable way of dating these young martian geological formations and the processes that modify them.     

Dense planetary rings and the viscous overstability

This paper examines the onset of the viscous overstability in dense particulate rings. First, we formulate a dense gas kinetic theory that is applicable to the saturnian system. Our model is essentially that of Araki and Tremaine [Araki, S., Tremaine, S., 1986. Icarus 65, 83-109], which we show can be both simplified and generalised. Second, we put this model to work computing the equilibrium properties of dense planetary rings, which we subsequently compare with the results of N-body simulations, namely those of Salo [Salo, H., 1991. Icarus 90, 254-270]. Finally, we present the linear stability analyses of these equilibrium states, and derive criteria for the onset of viscous overstability in the self-gravitating and non-self-gravitating cases. These are framed in terms of particle size, orbital frequency, optical depth, and the parameters of the collision law. Our results compare favourably with the simulations of Salo et al. [Salo, H., Schmidt, J., Spahn, F., 2001. Icarus 153, 295-315]. The accuracy and practicality of the continuum model we develop encourages its general use in future investigations of nonlinear phenomena.     

On Hapke photometric model predictions on reflectance of closely packed particulate surfaces

In a recent paper Hapke et al. (Hapke, B., Shepard, M., Nelson, R., Smythe, W., Piatek, J. [2009]. Icarus 199, 210-218) performed bi-directional reflectance measurements on closely-packed particulate surfaces of micrometer-sized particles and compared these with both the Hapke IMSA photometric model, and a numerical radiative transfer algorithm, the MDYZ (Mishchenko, M., Dlugach, J., Yanovitskij, E., Zakharova, N. [1999a]. J. Quant. Spectrosc. Radiat. Trans. 63, 409-432). To account for the effects of close packing, Hapke et al. applied a diffraction truncation scheme to remove the diffraction spike and supplied the renormalized single scattering phase function to the IMSA. They found that the IMSA prediction is a better match with measurement than that of MDYZ. In this work we demonstrate that the diffraction truncation procedure outlined by Hapke et al. contains an error. By following Hapke et al.’s intended truncation scheme, we have found that the IMSA model is not sufficiently anisotropic to describe the reflectance pattern of measurements on surface reflectance of closely packed large spherical particles.     

Photometric modeling of Asteroid 5535 Annefrank from Stardust observations

In this paper the Stardust disk-integrated phase curve at phase 47.2-134.6° of the Asteroid 5535 Annefrank, combined with groundbased observations (at phase 2.3-18.3°), are fit with Hapke’s photometric model. We confirm Newburn et al.’s (Newburn, R.L. et al. [2003]. J. Geophys. Res. 108 (E11), 5117. doi:10.1029/2003JE002106) observation that Annefrank exhibits a steep phase curve. This manifests itself in an unusually high fit surface roughness parameter of 49°. The single particle scattering albedo is 0.62, also high for an S-asteroid, while the fit phase function is more forward scattering than the typical S-asteroid being nearly isotropic with an asymmetry parameter of −0.09. The fit opposition surge width (h = 0.015) is typical of S-asteroids. However these fits assume a spherical shape to the asteroid. Li et al. (Li, J., A’Hearn, M.F., McFadden, L.A. [2004]. Icarus, 415-431) have shown that this assumption may lead to significant errors particularly at high phase angles leading to higher modeled single particle scattering albedos, macroscopic roughnesses and more forward scattering phase functions than actually exhibited. Our results confirm this finding—fitting only the data below 90° phase yields lower particle albedos (0.41) and roughnesses (20°) and more backscattering particles (−0.19) than the fit including the high phase angle data. Overall Annefrank appears to be on the bright side but otherwise is typical for an S-type asteroid suggesting that it may be a recent collisional fragment with a relatively immature surface which has had relatively little time to be weathered.     

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.