Full numerical simulations of catastrophic small body collisions

The outcome of collisions between small icy bodies, such as Kuiper belt objects, is poorly understood and yet a critical component of the evolution of the trans-neptunian region. The expected physical properties of outer Solar System materials (high porosity, mixed ice-rock composition, and low material strength) pose significant computational challenges. We present results from catastrophic small body collisions using a new hybrid hydrocode to N-body code computational technique. This method allows detailed modeling of shock propagation and material modification as well as gravitational reaccumulation. Here, we consider a wide range of material strengths to span the possible range of Kuiper belt objects. We find that the shear strength of the target is important in determining the collision outcome for 2 to 50-km radius bodies, which are traditionally thought to be in a pure gravity regime. The catastrophic disruption and dispersal criteria, , can vary by up to a factor of three between strong crystalline and weak aggregate materials. The material within the largest reaccumulated remnants experiences a wide range of shock pressures. The dispersal and reaccumulation process results in the material on the surfaces of the largest remnants having experienced a wider range of shock pressures compared to material in the interior. Hence, depending on the initial structure and composition, the surface materials on large, reaccumulated bodies in the outer Solar System may exhibit complex spectral and albedo variations. Finally, we present revised catastrophic disruption criteria for a range of impact velocities and material strengths for outer Solar System bodies.     

Thermal evolution of Kuiper belt objects, with implications for cryovolcanism

We investigate the internal thermal evolution of Kuiper belt objects (KBOs), small (radii <1000 km), icy (mean densities ) bodies orbiting beyond Neptune, focusing on Pluto's moon Charon in particular. Our calculations are time-dependent and account for differentiation. We review evidence for ammonia hydrates in the ices of KBOs, and include their effects on the thermal evolution. A key finding is that the production of the first melt, at the melting point of ammonia dihydrate, ≈176 K, triggers differentiation of rock and ice. The resulting structure comprises a rocky core surrounded by liquids and ice, enclosed within a >100-km thick undifferentiated crust of rock and ice. This structure is especially conducive to the retention of subsurface liquid, and bodies the size of Charon or larger (radii >600 km) are predicted to retain some subsurface liquid to the present day. We discuss the possibility that this liquid can be brought to the surface rapidly via self-propagating cracks. We conclude that cryovolcanism is a viable process expected to affect the surfaces of large KBOs including Charon.     

Photochemistry of methane-water ices

We report a study on the broadband ultraviolet photolysis of methane-water ice mixtures, at low methane concentrations and temperatures relevant to the icy satellites of the outer Solar System. The photochemistry of these mixtures is dominated by the action of hydroxyl radicals on methane and the resulting products. This implies that, given sufficient exposure time, the methane will eventually be completely oxidized to carbon dioxide. The presence of methane inhibits the formation of hydrogen peroxide by serving as a trap for hydroxyl radicals. The distribution of photochemical products is broadly similar to that previously conducted using ion and electron sources, with some differences possibly attributable to the difference in radiation source. The results are applicable to a variety of icy bodies in the Solar System. On Enceladus, where methane mixed with water is measured in the plumes, methane in the surface ices is subject to oxidation and will eventually be converted to CO₂. The CH stretch feature detected in the VIMS spectra of the Enceladus surface ice suggests that methane is currently being supplied to the surface ice, likely from re-condensation of the plume gas.     

Colors and collision rates within the Kuiper belt: problems with the collisional resurfacing scenario

We present a numerical check of the collisional resurfacing (CR) hypothesis proposed to explain the observed color diversity within the Kuiper Belt (where surface reddening due to space weathering is counteracted by regular resurfacing of neutral material after mutual collisions). Deterministic simulations are performed in order to estimate the relative spatial distribution of kinetic energy received by collisions, , for a population of target Kuiper Belt objects (KBOs) embedded in a swarm of impactors distributed within the belt. Four different impactor disks have been considered, depending on the excitation and the external limit of the belt and the density of the scattered KBOs (SKBOs) population. The obtained results are compared to the relative color index distribution within the observed Kuiper Belt, in order to derive possible similarities between the high vs low objects spatial distribution in our simulations and the bluer vs redder KBOs distribution in the “real” Kuiper Belt. Such similarities are found for several important features, in particular the general correlations between highly impacted objects and high rms excitation and low perihelion q values that are in good agreement with equivalent correlations found for the bluest objects of the observed belt. Nevertheless, simulations disagree with observations on two crucial points. (1) The plutinos are significantly more collisionally affected than the rest of our test KBO population, whereas there is no observed tendency toward bluer plutinos. (2) There is always a much stronger correlation between and eccentricities than inclinations, whereas observations show just the opposite feature. The presence of numerous SKBO impactors could significantly damp these problematic features, but cannot erase them. Whether these contradictions invalidate the whole CR scenario or not remains yet uncertain, since the physical processes at play are still far from being fully understood and the sample of available observational data is still relatively limited. But it seems nevertheless that the scenario might not hold in its simple present form.     

On the collisional disruption of porous icy targets simulating Kuiper belt objects

We present results from 27 impact experiments using porous (porosity ranging from 0.39 to 0.54) ice targets and solid ice projectiles at impact speeds ranging from 90 to 155 m/s. These targets were designed to simulate Kuiper Belt Objects (KBOs) in structure. We measured a specific energy for shattering, , of 2.1×¹⁰⁵ erg/g for those snowball targets hit by intact ice projectiles; this is of the same order as that measured for solid ice targets. The fragment mass distribution follows a power law, although the exponent is not simply related to the largest fragment size as assumed by fragmentation models. We provide the first measurement of the three-dimensional mass-velocity distribution for disrupted ice targets and find that fragment speeds range from ∼2 to ∼20 m/s. The fraction of collisional kinetic energy that is partitioned into ejecta speeds is between 1 and 15% (although it should be noted that the lower limit is more reliable than the upper).     

Ion irradiation of H2O ice on top of sulfurous solid residues and its relevance to the Galilean satellites

In this paper we present the results of new experiments of ion irradiation of water ice deposited on top of a solid sulfurous residue to study the potential formation of SO₂ at the interface ice/refractory material and discuss the possibility that this mechanism accounts for the sulfur dioxide ice detected on the surfaces of the Galilean satellites. In situ infrared spectroscopy was the used experimental technique. We have irradiated a thin film of H₂O frost on a sulfurous layer with 200 keV of He⁺ at 80 K. The used sulfurous residue was obtained by irradiation of frozen SO₂ at 16 K and it is used as a template of sulfur bearing solid materials. We have not found evidences of the efficient formation of SO₂ after irradiation of H₂O ice on top of the sulfurous residue. An upper limit to the production yield of SO₂, of interface area for each 100 eV of energy absorbed in 1 cm³ of ice-covered residue, has been estimated. These results have relevance in the context of the surfaces of the icy Galilean satellites in which SO₂ was detected. Our results show that radiolysis of mixtures of water ice and refractory sulfurous materials is not the primary formation mechanism responsible for the SO₂ present on the surfaces of the Galilean satellites.     

The orbit, mass, size, albedo, and density of (65489) Ceto/Phorcys: A tidally-evolved binary Centaur

Hubble Space Telescope observations of Uranus- and Neptune-crossing object (65489) Ceto/Phorcys (provisionally designated 2003 FX₁₂₈) reveal it to be a close binary system. The mutual orbit has a period of 9.554±0.011 days and a semimajor axis of 1840±48 km. These values enable computation of a system mass of (5.41±0.42)×¹⁰¹⁸ kg. Spitzer Space Telescope observations of thermal emission at 24 and 70 μm are combined with visible photometry to constrain the system's effective radius and geometric albedo . We estimate the average bulk density to be , consistent with ice plus rocky and/or carbonaceous materials. This density contrasts with lower densities recently measured with the same technique for three other comparably-sized outer Solar System binaries (617) Patroclus, (26308) 1998 SM₁₆₅, and (47171) 1999 TC₃₆, and is closer to the density of the saturnian irregular satellite Phoebe. The mutual orbit of Ceto and Phorcys is nearly circular, with an eccentricity ?0.015. This observation is consistent with calculations suggesting that the system should tidally evolve on a timescale shorter than the age of the Solar System.     

Orbits and masses of Saturn's coorbital and F-ring shepherding satellites

We have obtained numerically integrated orbits for Saturn's coorbital satellites, Janus and Epimetheus, together with Saturn's F-ring shepherding satellites, Prometheus and Pandora. The orbits are fit to astrometric observations acquired with the Hubble Space Telescope and from Earth-based observatories and to imaging data acquired from the Voyager spacecraft. The observations cover the 38 year period from the 1966 Saturn ring plane crossing to the spring of 2004. In the process of determining the orbits we have found masses for all four satellites. The densities derived from the masses for Janus, Epimetheus, Prometheus, and Pandora in units of g cm⁻³ are , , , and , respectively.     

Ultraviolet photolysis of amino acids in a 100 K water ice matrix: Application to the outer Solar System bodies

We report the rates of decomposition by ultraviolet (UV) photolysis of four amino acids in millimeter-thick crystalline water ice matrices at 100 K to constrain the survivability of these important organic molecules within ice lying near the surfaces of outer Solar System bodies. We UV-irradiated crystalline ice samples containing known concentrations of the amino acids glycine, aspartic acid, glutamic acid, and phenylalanine, then we measured the surviving concentrations using high performance liquid chromatography (HPLC) with fluorescence detection. From these experiments, we determine photolytic decomposition rates and half-lives. The half-life varies linearly with the ice thickness for all acids studied here. For example, glycine is the most resistant to photolytic destruction with a half-life of 50, 12, and 3.7 h in 1.6, 0.28, and 0.14 mm thick ices, respectively. We explain this linear variation of half-life with thickness as a consequence of extinction, mostly due to scattering, within these macroscopically thick ice samples. Applied to low latitude surface ice on Jupiter's satellite Europa, this analysis indicates that the concentration of any of these amino acids within the top meter of similar ice will be halved within a ∼10 year timescale.     

Orbital resonances in the inner neptunian system: II. Resonant history of Proteus, Larissa, Galatea, and Despina

We investigate the orbital history of the small neptunian satellites discovered by Voyager 2. Over the age of the Solar System, tidal forces have caused the satellites to migrate radially, bringing them through mean-motion resonances with one another. In this paper, we extend our study of the largest satellites Proteus and Larissa [Zhang, K., Hamilton, D.P., 2007. Icarus 188, 386-399] by adding in mid-sized Galatea and Despina. We test the hypothesis that these moons all formed with zero inclinations, and that orbital resonances excited their tilts during tidal migration. We find that the current orbital inclinations of Proteus, Galatea, and Despina are consistent with resonant excitation if they have a common density . Larissa's inclination, however, is too large to have been caused by resonant kicks between these four satellites; we suggest that a prior resonant capture event involving either Naiad or Thalassa is responsible. Our solution requires at least three past resonances with Proteus, which helps constrain the tidal migration timescale and thus Neptune's tidal quality factor: 9000<QN<36,000. We also improve our determination of Qs for Proteus and Larissa, finding 36<QP<700 and 18<QL<200. Finally, we derive a more general resonant capture condition, and work out a resonant overlap criterion relevant to satellite orbital evolution around an oblate primary.     

Radiation-induced amorphization of crystalline ice

We study radiation-induced amorphization of crystalline ice, analyzing the results of three decades of experiments with a variety of projectiles, irradiation energy, and ice temperature, finding a similar trend of increasing resistance of amorphization with temperature and inconsistencies in results from different laboratories. We discuss the temperature dependence of amorphization in terms of the ‘thermal spike’ model. We then discuss the common use of the 1.65 μm infrared absorption band of water as a measure of degree of crystallinity, an increasingly common procedure to analyze remote sensing data of astronomical icy bodies. The discussion is based on new, high quality near-infrared reflectance absorption spectra measured between 1.4 and 2.2 μm for amorphous and crystalline ices irradiated with 225 keV protons at 80 K. We found that, after irradiation with 10¹⁵ protons cm⁻², crystalline ice films thinner than the ion range become fully amorphous, and that the infrared absorption spectra show no significant changes upon further irradiation. The complete amorphization suggests that crystalline ice observed in the outer Solar System, including trans-neptunian objects, may results from heat from internal sources or from the impact of icy meteorites or comets.     

Cratering of icy targets by different impactors: Laboratory experiments and implications for cratering in the Solar System

Studies of impacts (impactor velocity about 5 km s⁻¹) on icy targets were performed. The prime goal was to study the response of solid CO₂ targets to impacts and to find the differences between the results of impacts on CO₂ targets with those on H₂O ice targets. The crater dimensions in CO₂ ice were found to scale with impact energy, with little dependence on projectile density (which ranged from nylon to copper, i.e., 1150-8930 kg m⁻³). At equal temperatures, craters in CO₂ ice were the same diameter as those in water ice, but were shallower and smaller in volume. In addition, the shape of the radial profiles of the craters was found to depend strongly on the type of ice and to change with impact energy. The impact speed of the data is comparable to that for impacts on many types of icy bodies in the outer Solar System (e.g., the satellites of the giant planets, the cometary nuclei and the Kuiper Belt objects), but the size and thus energy of the impactors is lower. Scaling with impact energy is demonstrated for the impacts on CO₂ ice. The issue of impact disruption (rather than cratering) is discussed by analogy with that on water ice. Expressions for the critical energy density for the onset of disruption rather than cratering are established for water ice as a function of porosity and silicate content. Although the critical energy density for disruption of CO₂ ice is not established, it is argued that the critical energy to disrupt a CO₂ ice body will be greater than that for a (non-porous) water ice body of the similar mass.     

Formation of radical species in photolyzed CH4:N2 ices

We report photochemical studies of thin cryogenic ice films composed of N₂, CH₄ and CO in ratios analogous to those on the surfaces of Neptune’s largest satellite, Triton, and on Pluto. Experiments were performed using a hydrogen discharge lamp, which provides an intense source of ultraviolet light to simulate the sunlight-induced photochemistry on these icy bodies. Characterization via infrared spectroscopy showed that C₂H₆ and C₂H₂, and HCO are formed by the dissociation of CH₄ into H, CH₂ and CH₃ and the subsequent reaction of these radicals within the ice. Other radical species, such as C₂, , CN, and CNN, are observed in the visible and ultraviolet regions of the spectrum. These species imply a rich chemistry based on formation of radicals from methane and their subsequent reaction with the N₂ matrix. We discuss the implications of the formation of these radicals for the chemical evolution of Triton and Pluto. Ultimately, this work suggests that , CN, HCO, and CNN may be found in significant quantities on the surfaces of Triton and Pluto and that new observations of these objects in the appropriate wavelength regions are warranted.     

Accurate absolute magnitudes for Kuiper belt objects and Centaurs

Accurate absolute optical magnitude values (HV and HR) for Kuiper belt objects (KBOs) and Centaurs are becoming increasingly important as observations in other wavelengths, particularly SIRTF thermal infrared measurements, become available for large samples of objects. We present accurate HV and HR values for 90 KBOs and Centaurs, based on our published optical photometry. We find that our HV values are in good agreement with those available from the European photometric survey of minor bodies in the outer Solar System. Comparison with HV values from the JPL Horizons database and the Minor Planet Center database shows that these sources are systematically brighter than ours by about 0.3 mag.     

Spectrophotometry of the small satellites of Saturn and their relationship to Iapetus, Phoebe, and Hyperion

We have obtained broadband spectrophotometric observations of four of the recently discovered small satellites of Saturn (Gladman et al., 2001, Nature 412, 163-166). The new data enable an understanding of the provenance, composition, and interrelationships among these satellites and the other satellites of Saturn, particularly Iapetus, Phoebe, and Hyperion. Temporal coverage of one satellite (S21 Tarvos) was sufficient to determine a partial rotational lightcurve. Our major findings include: (1) the satellites are red and similar in color, comparable to D-type asteroids, some KBOs, Iapetus, and Hyperion; (2) none of the satellites, including those from the “Phoebe Group” has any spectrophotometric relationship to Phoebe; and (3) S21 Tarvos exhibits a rotational lightcurve, although the data are not well-constrained and more observations are required to fit a definitive period. Dust created by meteoritic impacts and ejected from these satellites and additional undiscovered ones may be the source of the exogenous material deposited on the low-albedo side of Iapetus. Recent work which states that the small irregular satellites of Saturn have impacted Phoebe at least 6-7 times in the age of the Solar System (Nesvorny et al., 2003, Astron. J. 126, 398-429), suggests that such collisions may have propelled additional material from both Phoebe and the small irregular satellites toward Iapetus. The accretion of material from outer retrograde satellites may be a process that also occurs on Callisto and the uranian satellites.     

Orbital resonances in the inner neptunian system: I. The 2:1 Proteus-Larissa mean-motion resonance

We investigate the orbital resonant history of Proteus and Larissa, the two largest inner neptunian satellites discovered by Voyager 2. Due to tidal migration, these two satellites probably passed through their 2:1 mean-motion resonance a few hundred million years ago. We explore this resonance passage as a method to excite orbital eccentricities and inclinations, and find interesting constraints on the satellites' mean density () and their tidal dissipation parameters (Qs>10). Through numerical study of this mean-motion resonance passage, we identify a new type of three-body resonance between the satellite pair and Triton. These new resonances occur near the traditional two-body resonances between the small satellites and, surprisingly, are much stronger than their two-body counterparts due to Triton's large mass and orbital inclination. We determine the relevant resonant arguments and derive a mathematical framework for analyzing resonances in this special system.     

Embedded star clusters and the formation of the Oort Cloud

Observations suggest most stars originate in clusters embedded in giant molecular clouds [Lada, C.J., Lada, E.A., 2003. Annu. Rev. Astron. Astrophys. 41, 57-115]. Our Solar System likely spent 1-5 Myrs in such regions just after it formed. Thus the Oort Cloud (OC) possibly retains evidence of the Sun's early dynamical history and of the stellar and tidal influence of the cluster. Indeed, the newly found objects (90377) Sedna and 2000 CR₁₀₅ may have been put on their present orbits by such processes [Morbidelli, A., Levison, H.F., 2004. Astron. J. 128, 2564-2576]. Results are presented here of numerical simulations of the orbital evolution of comets subject to the influence of the Sun, Jupiter and Saturn (with their current masses on orbits appropriate to the period before the Late Heavy Bombardment (LHB) [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461]), passing stars and tidal force associated with the gas and stars of an embedded star cluster. The cluster was taken to be a Plummer model with 200-400 stars, with a range of initial central densities. The Sun's orbit was integrated in the cluster potential together with Jupiter and Saturn and the test particles. Stellar encounters were incorporated by directly integrating the effects of stars passing within a sphere centred on the Sun of radius equal to the Plummer radius for low-density clusters and half a Plummer radius for high-density clusters. The gravitational influence of the gas was modeled using the tidal force of the cluster potential. For a given solar orbit, the mean density, 〈ρ〉, was computed by orbit-averaging the density of material encountered. This parameter proved to be a good measure for predicting the properties of the OC. On average 2-18% of our initial sample of comets end up in the OC after 1-3 Myr. A comet is defined to be part of the OC if it is bound and has q>35 AU. Our models show that the median distance of an object in the OC scales approximately as 〈ρ〉^−1/2 when . Our models easily produce objects on orbits like that of (90377) Sedna [Brown, M.E., Trujillo, C., Rabinowitz, D., 2004. Astrophys. J. 617, 645-649] within ∼1 Myr in cases where the mean density is or higher; one needs mean densities of order to create objects like 2000 CR₁₀₅ by this mechanism, which are reasonable (see, e.g., Guthermuth, R.A., Megeath, S.T., Pipher, J.L., Williams, J.P., Allen, L.E., Myers, P.C., Raines, S.N., 2005. Astrophys. J. 632, 397-420). Thus the latter object may also be part of the OC. Close stellar passages can stir the primordial Kuiper Belt to sufficiently high eccentricities (e?0.05; Kenyon, S.J., Bromley, B.C., 2002. Astron. J. 123, 1757-1775) that collisions become destructive. From the simulations performed it is determined that there is a 50% or better chance to stir the primordial Kuiper Belt to eccentricities e?0.05 at 50 AU when . The orbit of the new object 2003 UB₃₁₃ [Brown, M.E., Trujillo, C.A., Rabinowitz, D.L., 2005. Astrophys. J. 635, L97-L100] is only reproduced for mean cluster densities of the order of , but in the simulations it could not come to be on its current orbit by this mechanism without disrupting the formation of bodies in the primordial Kuiper Belt down to 20 AU. It is therefore improbable that the latter object is created by this mechanism.     

Constraining the surface properties of Saturn's icy moons, using Cassini/CIRS emissivity spectra

We have conducted a search for emissivity features in the thermal infrared spectrum of the icy satellites of Saturn, Phoebe, Iapetus, Enceladus, Tethys, and Hyperion, observed by the Composite Infrared Spectrometer (CIRS) on board the Cassini spacecraft. Despite the heterogeneity of the composition of these bodies depicted by Earth-based and Cassini/VIMS observations, the CIRS spectra of all satellites are undistinguishable from black-body spectra, with no detectable emissivity feature. However, several materials, which have been detected on the surface of the same bodies, present emissivity features in the analyzed spectral range. In particular, water ice presents features with sufficient contrast to be detected by CIRS. Here we study the physical causes of the absence of features by simulating the effects of intimate mixtures using models of directional emissivity for optically thick surfaces for different particle sizes and abundances, and porosities. The simulations include a set of materials detected on the Phoebe's surface, like water ice, hydrated silicates, and organics. We find that featureless spectra can be produced in three scenarios: (1) ice particles with large sizes, (2) mixtures of ices dominated by dark contaminants, and (3) small particles with large porosity. Constraints imposed by the NIR spectra of the satellites favors the latter scenario as the more likely explanation to the absence of emissivity features on the icy satellites of Saturn.     

Thermal inertia and bolometric Bond albedo values for Mimas, Enceladus, Tethys, Dione, Rhea and Iapetus as derived from Cassini/CIRS measurements

Spectra taken by Cassini’s Composite Infrared Spectrometer (CIRS) between 10 and 600 cm⁻¹ (17-1000 μm) of surface thermal emission of Mimas, Enceladus, Tethys, Dione, Rhea and Iapetus have been used to derive the thermal inertia and bolometric Bond albedo values. Only an upper limit for the bolometric Bond albedo of Iapetus’ dark leading side could be determined due to the insensitivity of the thermal model to albedo when albedos are very low. The thermal inertia in this region however is better constrained. The CIRS coverage of Enceladus is extensive enough that the latitudinal variation in these values from 60°S to 70°N has been determined in 10° wide bins. The bolometric Bond albedos determined here are consistent with literature values which show the surface of the saturnian icy moons to be covered in ice contaminated to varying degrees. The thermal inertia of the moons is shown to be in the range 9-, approximately 2-6 times lower than that of the Galilean satellites, implying a less well consolidated and more porous surface. The thermal inertias of Iapetus and Phoebe are somewhat higher, suggesting that the very low thermal inertias of satellites from Rhea inwards may be related to their probable coating of E-ring material. Latitudinal variations on the surface of Enceladus show that the bolometric Bond albedo and thermal inertia increase towards the active plume source at the south pole.     

The global shape of Europa: Constraints on lateral shell thickness variations

The global shape of Europa is controlled by tidal and rotational potentials and possibly by lateral variations in ice shell thickness. We use limb profiles from four Galileo images to determine the best-fit hydrostatic shape, yielding a mean radius of 1560.8±0.3 km and a radius difference a−c of 3.0±0.9 km, consistent with previous determinations and inferences from gravity observations. Adding long-wavelength topography due to proposed lateral variations in shell thickness results in poorer fits to the limb profiles. We conclude that lateral shell thickness variations and long-wavelength isostatically supported topographic variations do not exceed 7 and 0.7 km, respectively. For the range of rheologies investigated (basal viscosities from 10¹⁴ to ) the maximum permissible (conductive) shell thickness is 35 km. The relative uniformity of Europa's shell thickness is due to either a heat flux from the silicate interior, lateral ice flow at the base of the shell, or convection within the shell.