Characterization of Titan’s Ontario Lacus region from Cassini/VIMS observations

Liquid hydrocarbons were long predicted on Titan’s surface before the RADAR instrument onboard Cassini detected lakes poleward of 70°N in July 2006. Before that the Cassini Imaging Science Subsystem (ISS) observed a lake-like feature in the South Pole, named Ontario Lacus, in July 2004. Here we analyze one observation of Ontario Lacus taken by the Visual and Infrared Mapping Spectrometer (VIMS) on 2007 December 5, during the T 38 flyby. This is the best spatially resolved image of a Titan lake to date by an imaging spectrometer, and has been previously reported in Brown et al. (Brown, R.H., Soderblom, L.A., Soderblom, J.M., Clark, R.N., Jaumann, R., Barnes, J.W., Sotin, C., Buratti, B., Baines, K.H., Nicholson, P.D. [2008]. Nature 454, 607–610) and in Barnes et al. (Barnes, J.W. et al. [2009]. Icarus 201, 217–225). The observing geometry and our data processing will be explained, followed by a discussion of the main characteristics of the image. The analyzed image covers a small portion of Ontario Lacus and shows what appears from RADAR data to be a region of modest slope (“ramp”) adjacent to the dark lake itself. Our analysis of 5.0 μm spectral data suggests that the previously reported absorption feature of ethane seen at shorter wavelengths may be produced by damp sediments adjacent to the main liquid basin. The latter appears to be absorbing all of the photons incident upon it in the 5 μm spectral region and shows no discernible absorption features. A characterization of the basin composition and morphology is developed with the help of ISS and RADAR observations. The simplest model consistent with the data is an optically deep lake surrounded by a region in which ethane, propane, possibly methane, and other, less volatile hydrocarbons and nitriles are present mixed into spectroscopically neutral sediments. The dominance of relatively low vapor pressure organics outside the lake itself suggests a retreat of Ontario Lacus associated with evaporation on seasonal or longer timescales, consistent with analysis of RADAR and ISS images.     

Solar phase curves and phase integrals for the leading and trailing hemispheres of Iapetus from the Cassini Visual Infrared Mapping Spectrometer

We performed photometry of Cassini Visual Infrared Mapping Spectrometer observations of Iapetus to produce the first phase integrals calculated directly from solar phase curves of Iapetus for the leading hemisphere and to estimate the phase integrals for the trailing hemisphere. We also explored the phase integral dependence on wavelength and geometric albedo. The extreme dichotomy of the brightness of the leading and trailing sides of Iapetus is reflected in their phase integrals. Our phase integrals, which are lower than the results of Morrison et al. (Morrison, D., Jones, T.J., Cruikshank, D.P., Murphy, R.E. [1975]. Icarus 24, 157-171) and Squyres et al. (Squyres, S.W., Buratti, B.J., Veverka, J., Sagan, C. [1984]. Icarus 59, 426-435), have profound implications on the energy balance and volatile transport on this icy satellite.     

A newly discovered impact crater in Titan's Senkyo: Cassini VIMS observations and comparison with other impact features

Senkyo is an equatorial plain on Titan filled with dunes and surrounded by hummocky plateaus. During the Titan targeted flyby T61 on August 25, 2009, the Cassini Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft observed a circular feature, centered at 5.4° N and 341°W, that superimposes the dune fields and a bright plateau. This circular feature, which has been named Paxsi by the International Astronomical Union, is 120±10 km in diameter (measured from the outer edge of the crater rim) and exhibits a central bright area that can be interpreted as the central peak or pit of an impact crater. Although there are only a handful of certain impact craters on Titan, there are two other craters that are of similar size to this newly discovered feature and that have been studied by VIMS: Sinlap (Le Mouélic et al., 2008) and Selk (Soderblom et al., 2010). Sinlap is associated with a large downwind, fan-like feature that may have been formed from an impact plume that rapidly expanded and deposited icy particles onto the surface. Although much of the surrounding region is covered with dunes, the plume region is devoid of dunes. The formation process of Selk also appears to have removed (or covered up) dunes from parts of the adjacent dune-filled terrain. The circular feature on Senkyo is quite different: there is no evidence of an ejecta blanket and the crater itself appears to be infilled with dune material. The rim of the crater appears to be eroded by fluvial processes; at one point the rim is breached. The rim is unusually narrow, which may be due to mass wasting on its inside and subsequent infill by dunes. Based on these observations, we interpret this newly discovered feature to be a more eroded crater than both Sinlap and Selk. Paxsi may have formed during a period when Titan was warmer and more ductile than it is currently.     

Radarclinometry of the sand seas of Africa’s Namibia and Saturn’s moon Titan

Radarclinometry is a powerful technique for estimating heights of landforms in synthetic aperture radar (SAR) images of planetary surfaces. In particular, it has been used to estimate heights of dunes in the sand seas of Saturn’s moon Titan (Lorenz, R.D., and 39 colleagues [2006]. Science 312, 724-727). In this work, we verify the technique by comparing dune heights derived from radarclinometry to known topography of dune fields in the Namib sand sea of western Africa. We compared results from three different image grid spacings, and found that 350 m/pixel (the same spacing at which the Cassini RADAR data was processed) is sufficient to determine dune height for dunes of similar morphometry to those of the Namib sand sea. At this grid spacing, height estimates derived from radarclinometry are largely representative of, though may underestimate by as much as 30%, or overestimate by as much as 40%, true dune height. Applying the technique to three regions on Titan, we estimate dune heights of 45-180 m, and dune spacings of 2.3-3.3 km. Obtaining accurate heights of Titan’s dunes will help to constrain the total organic inventory on Titan.     

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.     

Anomalous radar backscatter from Titan’s surface?

Since Cassini arrived at Saturn in 2004, its moon Titan has been thoroughly mapped by the RADAR instrument at 2-cm wavelength, in both active and passive modes. Some regions on Titan, including Xanadu and various bright hummocky bright terrains, contain surfaces that are among the most radar-bright encountered in the Solar System. This high brightness has been generally attributed to volume scattering processes in the inhomogeneous, low-loss medium expected for a cold, icy satellite surface. We can test this assumption now that the emissivity has been obtained from the concurrent radiometric measurements for nearly all the surface, with unprecedented accuracy (Janssen et al., and the Cassini RADAR Team [2009]. Icarus 200, 222-239). Kirchhoff’s law of thermal radiation relates the radar and radiometric properties in a way that has never been fully exploited. In this paper we examine here how this law may be applied in this case to better understand the nature of Titan’s radar-bright regions. We develop a quantitative model that, when compared to the observational data, allows us to conclude that either the reflective characteristics of the putative volume scattering subsurface must be highly constrained, or, more likely, organized structure on or in the surface is present that enhances the backscatter.     

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.     

Velocity and vorticity measurements of Jupiter's Great Red Spot using automated cloud feature tracking

We have produced mosaics of the Great Red Spot (GRS) using images taken by the Galileo spacecraft in May 2000, and have measured the winds of the GRS using an automated algorithm that does not require manual cloud tracking. Our technique yields a high-density, regular grid of wind velocity vectors that is advantageous over a limited number of scattered wind vectors that result from manual cloud tracking. The high-velocity collar of the GRS is clearly seen from our velocity vector map, and highest wind velocities are measured to be around 170 m s⁻¹. The high resolution of the mosaics has also enabled us to map turbulent eddies inside the chaotic central region of the GRS, similar to those mapped by Sada et al. [Sada, P.V., Beebe, R.F., Conrath, B.J., 1996. Icarus 119, 311-335]. Using the wind velocity measurements, we computed particle trajectories around the GRS as well as maps of relative and absolute vorticities. We have discovered a narrow ring of cyclonic vorticity that surrounds the main anti-cyclonic high-velocity collar. This narrow ring appears to correspond to a ring surrounding the GRS that is bright in 5 μm [Terrile, R.J., Beebe, R.F., 1979. Science 204, 948-951]. It appears that this cyclonic ring is not a transient feature of the GRS, as we have discovered it in a re-analysis of Galileo data taken in 1996 first analyzed by Vasavada et al. [Vasavada, A.R., and 13 colleagues, 1998. Icarus 135, 265-275]. We also calculate how absolute vorticity changes as a function of latitude along a trajectory around the GRS and compare these measurements to similar ones performed by Dowling and Ingersoll [Dowling, T.E., Ingersoll, A.P., 1988. J. Atmos. Sci. 45, 1380-1396] using Voyager data. We show no dramatic evolution in the structure of the GRS since the Voyager era except for additional evidence for a counter-rotating GRS core, an increase in velocity in the main velocity collar, and an overall decrease in the length of the GRS.     

On Titan’s Xanadu region

A large, circular marking ∼1800 km across is seen in near-infrared images of Titan. The feature is centered at 10°S, 120°W on Titan, encompasses much of Titan’s western Xanadu region, and has an off-center, quasi-circular, inner margin about 700 km across, with lobate outer margins extending 200-500 km from the inner margin. On the feature’s southern flank is Tui Regio, an area that has very high reflectivity at 5 μm, and is hypothesized to exhibit geologically recent cryovolcanic flows (Barnes, J.W. et al. [2006]. Geophys. Res. Lett. 33), similar to flows seen in Hotei Regio, a cryovolcanic area whose morphology may be controlled by pre-existing, crustal fractures resulting from an ancient impact (Soderblom, L.A. et al. [2009]. Icarus, 204). The spectral reflectivity of the large, circular feature is quite different than that of its surroundings, making it compositionally distinct, and radar measurements of its topography, brightness temperature and volume scattering also suggest that the feature is quite distinct from its surroundings. These and several other lines of evidence, in addition to the feature’s morphology, suggest that it may occupy the site of an ancient impact.     

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.     

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.     

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.     

Comment on “Evidence of electrical activity on Titan drawn from the Schumann resonances sent by Huygens probe” by J.A. Morente, J.A. Portí, A. Salinas, E.A. Navarro [2008, Icarus, 195, 802-811]

Morente et al. [Morente, J.A., Portí, J.A., Salinas, A., Navarro, E.A., 2008. Icarus. doi:10.1016/j.icarus.2008.02.004] have recently presented a new analysis of the Permittivity, Wave and Altimetry (PWA) measurements made during the descent of the Huygens Probe through the atmosphere of Titan. They claimed the identification of several Schumann resonance harmonics and concluded in favor of a lightning activity on Titan. We report here several reasons for not endorsing this paper.     

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.     

Titan's 3-micron spectral region from ISO high-resolution spectroscopy

The near-infrared spectrum of Titan, Saturn's largest moon and one of the Cassini/Huygens' space mission primary targets, covers the 0.8 to 5 micron region in which it shows several weak CH₄ absorption regions, and in particular one centered near 2.75 micron. Due to the interference of telluric absorption, only part of this window region (2.9-3.1 μm) has previously been observed from the ground [Noll, K.S., Geballe, T.R., Knacke, R., Pendleton, F., Yvonne, J., 1996. Icarus 124, 625-631; Griffith, C.A., Owen, T., Miller, G.A., Geballe, T., 1998. Nature 395, 575-578; Griffith, C.A., Owen, T., Geballe, T.R., Rayner, J., Rannou, P., 2003. Science 300, 628-630; Geballe, T.R., Kim, S.J., Noll, K.S., Griffith, C.A., 2003. Astrophys. J. 583, L39-L42]. We report here on the first spectroscopic observations of Titan covering the whole 2.4-4.9 μm region by two instruments on board the Infrared Space Observatory (ISO) in 1997. These observations show the 2.75-μm window in its complete extent for the first time. In this study we have also used a high-resolution Titan spectrum in the 2.9-3.6 μm region taken with the Keck [Geballe, T.R., Kim, S.J., Noll, K.S., Griffith, C.A., 2003. Astrophys. J. 583, L39-L42; Kim, S.J., Geballe, T.R., Noll, K.S., Courtin, R., 2005. Icarus 173, 522-532] to infer information on the atmospheric parameters (haze extinction, single scattering albedo, methane abundance, etc.) by fitting the methane bands with a detailed microphysical model of Titan's atmosphere (updated from Rannou, P., McKay, C.P., Lorenz, R.D., 2003. Planet. Space Sci. 51, 963-976). We have included in this study an updated version of a database for the CH₄ absorption coefficients [STDS, Wenger, Ch., Champion, J.-P., 1998. J. Quant. Spectrosc. Radiat. Transfer 59, 471-480. See also http://www.u-bourgogne.fr/LPUB/TSM/sTDS.html for latest updates; Boudon, V., Champion, J.-P., Gabard, T., Loëte, M., Michelot, F., Pierre, G., Rotger, M., Wenger, Ch., Rey, M., 2004. J. Mol. Spectrosc. 228, 620-634]. For the atmosphere we find that (a) the haze extinction profile that best matches the data is one with higher (by 40%) extinction in the atmosphere with respect to Rannou et al. (2003) down to about 30 km where a complete cut-off occurs; (b) the methane mixing ratio at Titan's surface cannot exceed 3% on a disk-average basis, yielding a maximum CH₄ column abundance of 2.27 km-am in Titan's atmosphere. From the derived surface albedo spectrum in the 2.7-3.08 micron region, we bring some constraints on Titan's surface composition. The albedo in the center of the methane window varies from 0.01 to 0.08. These values, compared to others reported in the other methane windows, show a strong compatibility with the water ice spectrum in the near-infrared. Without confirming its existence from this work alone, our data then appear to be compatible with water ice. A variety of other ices, such as CO₂, NH₃, tholin material or hydrocarbon liquid cannot be excluded from our data, but an additional unidentified component with a signature around 2.74 micron is required to satisfy the data.     

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.     

UV spectrum of Enceladus

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

Optical properties of Titan and early Earth haze laboratory analogs in the mid-visible

Scattering and absorption of sunlight by aerosols are integral to understanding the radiative balance of any planetary atmosphere covered in a haze, such as Titan and possibly the early Earth. One key optical parameter of an aerosol is its refractive index. We have simulated both Titan and early Earth organic haze aerosols in the laboratory and measured the real and imaginary portion of their refractive index at λ = 532 nm using cavity ringdown aerosol extinction spectroscopy. This novel technique allows analysis on freely-floating particles minutes after formation. For our Titan analog particles, we find a real refractive index of n = 1.35 ± 0.01 and an imaginary refractive index k = 0.023 ± 0.007, and for the early Earth analog particles we find n = 1.81 ± 0.02 and k = 0.055 ± 0.020. The Titan analog refractive index has a smaller real and similar imaginary refractive index compared to most previous laboratory measurements of Titan analog films, including values from Khare et al. (Khare, B.N., Sagan, C., Arakawa, E.T., Suits, F., Callcott, T.A., Williams, M.W. [1984]. Icarus 60, 127-137). These newly measured Titan analog values have implications for spacecraft retrievals of aerosol properties on Titan. The early Earth analog has a significantly higher real and imaginary refractive index than Titan analogs reported in the literature. These differences suggest that, for a given amount of aerosol, the early Earth analog would act as a stronger anti-greenhouse agent than the Titan analog.     

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.