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.     

Stratigraphical evidence of late Amazonian periglaciation and glaciation in the Astapus Colles region of Mars

Recent modeling of the meteorological conditions during and following times of high obliquity suggests that an icy mantle could have been emplaced in western Utopia Planitia by atmospheric deposition during the late Amazonian period [Costard, F.M., Forget, F., Madeleine, J.B., Soare, R.J., Kargel, J.S., 2008. Lunar Planet. Sci. 39. Abstract 1274; Madeleine, B., Forget, F., Head, J.W., Levrard, B., Montmessin, F., 2007. Lunar Planet. Sci. 38. Abstract 1778]. Astapus Colles (ABa) is a late Amazonian geological unit — located in this hypothesized area of accumulation — that comprises an icy mantle tens of meters thick [Tanaka, K.L., Skinner, J.A., Hare, T.M., 2005. US Geol. Surv. Sci. Invest., Map 2888]. For the most part, this unit drapes the early Amazonian Vastitas Borealis interior unit (ABv_i); to a lesser degree it overlies the early Amazonian Vastitas Borealis marginal unit (ABv_m) and the early to late Hesperian UP plains unit HBu₂ [Tanaka, K.L., Skinner, J.A., Hare, T.M., 2005. US Geol. Surv. Sci. Invest., Map 2888]. Landscapes possibly modified by late-Amazonian periglacial processes [Costard, F.M., Kargel, J.S., 1995. Icarus 114, 93-112; McBride, S.A., Allen, C.C., Bell, M.S., 2005. Lunar Planet. Sci. 36. Abstract 1090; Morgenstern, A., Hauber, E., Reiss, D., van Gasselt, S., Grosse, G., Schirrmeister, L., 2007. J. Geophys. Res. 112, doi:10.1029/2006JE002869. E06010; Seibert, N.M., Kargel, J.S., 2001. Geophys. Res. Lett. 28, 899-902; Soare, R.J., Kargel, J.S., Osinski, G.R., Costard, F., 2007. Icarus 191, 95-112; Soare, R.J., Osinski, G.R., Roehm, C.L., 2008. Earth Planet. Sci. Lett. 272, 382-393] and glacial processes [Milliken, R.E., Mustard, J.F., Goldsby, D.L., 2003. J. Geophys. Res. 108 (E6), doi:10.1029/2002JE002005. 5057; Mustard, J.F., Cooper, C.D., Rifkin, M.K., 2001. Nature 412, 411-414; Tanaka, K.L., Skinner, J.A., Hare, T.M., 2005. US Geol. Surv. Sci. Invest., Map 2888] have been reported within the region. Researchers have assumed that the periglacial and glacial landscapes occur within the same geological unit, the ABa [i.e., Morgenstern, A., Hauber, E., Reiss, D., van Gasselt, S., Grosse, G., Schirrmeister, L., 2007. J. Geophys. Res. 112; doi:10.1029/2006JE002869. E06010; Tanaka, K.L., Skinner, J.A., Hare, T.M., 2005. US Geol. Surv. Sci. Invest., Map 2888]. In this study we use HiRISE (High Resolution Image Science Experiment, Mars Reconnaissance Orbiter) imagery to identify the stratigraphical separation of the two landscapes and show that periglacial landscape modification has occurred in the geological units that underlie the ABa, not in the ABa itself. Moreover, we suggest that the periglacial landscape extends well beyond the perimeter of the ABa and could be the product of “wet” cold-climate processes. These processes involve freeze-thaw cycles and intermittently stable liquid-water at or near the surface. By contrast, we propose that the ABa is a very recent late-Amazonian geological unit formed principally by “dry” cold-climate processes. These processes comprise accumulation (by atmospheric deposition) and ablation (by sublimation).     

Dielectric constant estimation of the uppermost Basal Unit layer in the martian Boreales Scopuli region

An electromagnetic inversion model has been applied to echoes from the subsurface sounding Shallow Radar (SHARAD) to retrieve the dielectric properties of the uppermost Basal Unit (BU) beneath the North Polar Layered Deposits of Mars. SHARAD data have been carefully selected to satisfy the assumption of the inversion model which requires a stratigraphy consisting of mostly plane parallel layers. The resulting values of the dielectric constant have been interpreted in terms of a variable percentage of dust in an ice–dust mixture through the use of a mixing model for dielectric properties. The resulting dust content exceeds 65%, reaching perhaps 95%, depending on the permittivity values assumed for the dust. Such a concentration is higher than that obtained by Selvans et al. (Selvans, M.M., Plaut, J.J., Aharonson, O. [2010]. J. Geophys. Res, 115, E09003). This discrepancy could be justified considering that our observations refer to the uppermost BU layer, whereas Selvans et al. (Selvans, M.M., Plaut, J.J., Aharonson, O. [2010]. J. Geophys. Res, 115, E09003) probed the BU full thickness. Moreover, if the BU is considered spatially inhomogeneous, with very different dust content and thickness (Tanaka, K.L., Skinner, J.A., Fortezzo, C.M., Herkenhoff, K.E., Rodriguez, J.A.P., Bourke, M.C., Kolb, E.J., Okubo, C.H. [2008]. Icarus, 196, 318–358), the discrepancy could be furtherly reconciled.     

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).     

How the Enceladus dust plume feeds Saturn’s E ring

Pre-Cassini models of Saturn’s E ring [Horányi, M., Burns, J., Hamilton, D., 1992. Icarus 97, 248-259; Juhász, A., Horányi, M., 2002. J. Geophys. Res. 107, 1-10] failed to reproduce its peculiar vertical structure inferred from Earth-bound observations [de Pater, I., Martin, S.C., Showalter, M.R., 2004. Icarus 172, 446-454]. After the discovery of an active ice-volcanism of Saturn’s icy moon Enceladus the relevance of the directed injection of particles for the vertical ring structure of the E ring was swiftly recognised [Juhász, A., Horányi, M., Morfill, G.E., 2007. Geophys. Res. Lett. 34, L09104; Kempf, S., Beckmann, U., Moragas-Klostermeyer, G., Postberg, F., Srama, R., Economou, T., Schmidt, J., Spahn, F., Grün, E., 2008. Icarus 193, 420-437]. However, simple models for the delivery of particles from the plume to the ring predict a too small vertical ring thickness and overestimate the amount of the injected dust.Here we report on numerical simulations of grains leaving the plume and populating the dust torus of Enceladus. We run a large number of dynamical simulations including gravity and Lorentz force to investigate the earliest phase of the ring particle life span. The evolution of the electrostatic charge carried by the initially uncharged grains is treated selfconsistently. Freshly ejected plume particles are moving in almost circular orbits because the Enceladus orbital speed exceeds the particles’ ejection speeds by far. Only a small fraction of grains that leave the Hill sphere of Enceladus survive the next encounter with the moon. Thus, the flux and size distribution of the surviving grains, replenishing the ring particle reservoir, differs significantly from the flux and size distribution of the particles freshly ejected from the plume. Our numerical simulations reproduce the vertical ring profile measured by the Cassini Cosmic Dust Analyzer (CDA) [Kempf, S., Beckmann, U., Moragas-Klostermeyer, G., Postberg, F., Srama, R., EconoDmou, T., Smchmidt, J., Spahn, F., Grün, E., 2008. Icarus 193, 420-437]. From our simulations we calculate the deposition rates of plume particles hitting Enceladus’ surface. We find that at a distance of 100 m from a jet a 10 m sized ice boulder should be covered by plume particles in 10⁵-10⁶ years.     

Ancient melting of mid-latitude snowpacks on Mars as a water source for gullies

We hypothesize that during past epochs of high obliquity seasonal snowfields at mid-latitudes melted to produce springtime sediment-rich surface flows resulting in gully formation. Significant seasonal mid-latitude snowfall does not occur on Mars today. General Circulation Model (GCM) results, however, suggest that under past climate conditions there may have been centimeters of seasonal mid-latitude snowfall [Mischna, M.A., Richardson, M.I., Wilson, R.J., McCleese, D.J., 2003. J. Geophys. Res. Planets 108, doi:10.1029/2003JE002051. 5062]. Gully locations have been tabulated by several researchers (e.g. [Heldmann, J.L., Mellon, M.T., 2004. Icarus 168, 285–304; Heldmann, J.L., Carlsson, E., Johansson, H., Mellon, M.T., Toon, O.B., 2007. Icarus 188, 324–344; Malin, M.C., Edgett, K.S., 2000. Science 288, 2330–2335]) and found to correspond to mid-latitude bands. A natural question is whether the latitudinal bands where the gullies are located correspond to areas where the ancient snowfalls may have melted, producing runoff which may have incised gullies. In this study we model thin snowpacks with thicknesses similar to those predicted by [Mischna, M.A., Richardson, M.I., Wilson, R.J., McCleese, D.J., 2003. J. Geophys. Res. Planets 108, doi:10.1029/2003JE002051. 5062]. We model these snowpacks under past climate regimes in order to determine whether snowmelt runoff could have occurred, and whether significant amounts of warm soil (T>273 K) existed on both poleward and equatorward slopes in the regions where gullies exist. Both warm soil and water amounts are modeled because soil and water may have mixed to form a sediment-rich flow. We begin by applying the snowpack model of Williams et al. [Williams, K.E., Toon, O.B., Heldmann, J.E., Mellon, M., 2008. Icarus 196, 565–577] to past climate regimes characterized by obliquities of 35° (600 ka before present) and 45° (5.5 ma before present), and to all latitudes between 70° N and 70° S. We find that the regions containing significant snowmelt runoff correspond to the regions identified by Heldmann and Mellon [Heldmann, J.L., Mellon, M.T., 2004. Icarus 168, 285–304], Heldmann et al. [Heldmann, J.L., Carlsson, E., Johansson, H., Mellon, M.T., Toon, O.B., 2007. Icarus 188, 324–344] and Malin and Edgett [Malin, M.C., Edgett, K.S., 2000. Science 288, 2330–2335] as containing large numbers of gullies. We find that the snowmelt runoff (>1 mm, with equivalent rainfall rates of 0.25 mm/h) and warm soil (>1 cm depth) would have occurred on slopes within the gullied latitudinal bands. The snowfall amounts modeled are predicted to be seasonal [Mischna, M.A., Richardson, M.I., Wilson, R.J., McCleese, D.J., 2003. J. Geophys. Res. Planets 108, doi:10.1029/2003JE002051. 5062], and our modeling finds that under the previous climate regimes there would have been meltwater present on the slopes in question for brief periods of time, on the order of days, each year. Our model provides a simple explanation for the latitudinal distribution of the gullies, and also suggests that the gullies date to times when water migrated away from the present poles to the mid-latitudes.     

Investigations of the variability of dust particle sizes in the martian atmosphere using the NASA Ames General Circulation Model

We present a Mars General Circulation Model (GCM) numerical investigation of the physical processes (i.e., wind stress and dust devil dust lifting and atmospheric transport) responsible for temporal and spatial variability of suspended dust particle sizes. Measurements of spatial and temporal variations in airborne dust particles sizes in the martian atmosphere have been derived from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) spectral and emission phase function data [Wolff, M.J., Clancy, R.T., 2003. J. Geophys. Res. (Planets) 108 (E9), doi:10.1029/2003JE002057. 1-1; Clancy, R.T., Wolff, M.J., Christensen, P.R., 2003. J. Geophys. Res. (Planets) 108 (E9), doi:10.1029/2003JE002058. 2-1]. The range of dust particle sizes simulated by the NASA Ames GCM is qualitatively consistent with TES-derived observations of effective dust particle size variability. Model results suggest that the wind stress dust lifting scheme (which produces regionally confined dust lifting) is the process responsible for the majority of the dust particle size variability in the martian atmosphere. Additionally, model results suggest that atmospheric transport processes play an important role in the evolution of atmospheric dust particles sizes during substantial dust storms on Mars. Finally, we show that including the radiative effects of a spatially variable particle size distribution significantly influences thermal and dynamical fields during the dissipation phase of the simulated global dust storm.     

The effects of the martian regolith on GCM water cycle simulations

This paper describes General Circulation Model (GCM) simulations of the martian water cycle focusing on the effects of an adsorbing regolith. We describe the 10-layer regolith model used in this study which has been adapted from the 1-D model developed by Zent, A.P., Haberle, R.M., Houben, H.C., Jakosky, B.M. [1993. A coupled subsurface-boundary layer model of water on Mars. J. Geophys. Res. 98 (E2), 3319-3337, February]. Even with a 30-min timestep and taking into account the effect of surface water ice, our fully implicit scheme compares well with the results obtained by Zent, A.P., Haberle, R.M., Houben, H.C., Jakosky, B.M. [1993. A coupled subsurface-boundary layer model of water on Mars. J. Geophys. Res. 98 (E2), 3319-3337, February]. This means, however, that the regolith is not able to reproduce the diurnal variations in column water vapour abundance of up to a factor of 2-3 as seen in some observations, with only about 10% of the atmospheric water vapour column exchanging with the subsurface on a daily basis. In 3-D simulations we find that the regolith adsorbs water preferentially in high latitudes. This is especially true in the northern hemisphere, where perennial subsurface water ice builds up poleward of 60° N at depths which are comparable to the Odyssey observations. Much less ice forms in the southern high latitudes, which suggests that the water ice currently present in the martian subsurface is not stable under present conditions and is slowly subliming and being deposited in the northern hemisphere. When initialising the model with an Odyssey-like subsurface water ice distribution the model is capable of forcing the simulated water cycle from an arbitrary state close to the Mars Global Surveyor Thermal Emission Spectrometer observations. Without the actions of the adsorbing regolith the equilibrated water cycle is found to be a factor of 2-4 too wet. The process by which this occurs is by adsorption of water during northern hemisphere summer in northern mid and high latitudes where it remains locked in until northern spring when the seasonal CO₂ ice cap retreats. At this time the water diffuses out of the regolith in response to increased temperature and is returned to the residual water ice cap by eddie transport.     

Subsurface structure of Planum Boreum from Mars Reconnaissance Orbiter Shallow Radar soundings

We map the subsurface structure of Planum Boreum using sounding data from the Shallow Radar (SHARAD) instrument onboard the Mars Reconnaissance Orbiter. Radar coverage throughout the 1,000,000-km² area reveals widespread reflections from basal and internal interfaces of the north polar layered deposits (NPLD). A dome-shaped zone of diffuse reflectivity up to 12 μs (∼1-km thick) underlies two-thirds of the NPLD, predominantly in the main lobe but also extending into the Gemina Lingula lobe across Chasma Boreale. We equate this zone with a basal unit identified in image data as Amazonian sand-rich layered deposits [Byrne, S., Murray, B.C., 2002. J. Geophys. Res. 107, 5044, 12 pp. doi:10.1029/2001JE001615; Fishbaugh, K.E., Head, J.W., 2005. Icarus 174, 444-474; Tanaka, K.L., Rodriguez, J.A.P., Skinner, J.A., Bourke, M.C., Fortezzo, C.M., Herkenhoff, K.E., Kolb, E.J., Okubo, C.H., 2008. Icarus 196, 318-358]. Elsewhere, the NPLD base is remarkably flat-lying and co-planar with the exposed surface of the surrounding Vastitas Borealis materials. Within the NPLD, we delineate and map four units based on the radar-layer packets of Phillips et al. [Phillips, R.J., and 26 colleagues, 2008. Science 320, 1182-1185] that extend throughout the deposits and a fifth unit confined to eastern Gemina Lingula. We estimate the volume of each internal unit and of the entire NPLD stack (821,000 km³), exclusive of the basal unit. Correlation of these units to models of insolation cycles and polar deposition [Laskar, J., Levrard, B., Mustard, J.F., 2002. Nature 419, 375-377; Levrard, B., Forget, F., Montmessin, F., Laskar, J., 2007. J. Geophys. Res. 112, E06012, 18 pp. doi:10.1029/2006JE002772] is consistent with the 4.2-Ma age of the oldest preserved NPLD obtained by Levrard et al. [Levrard, B., Forget, F., Montmessin, F., Laskar, J., 2007. J. Geophys. Res. 112, E06012, 18 pp. doi:10.1029/2006JE002772]. We suggest a dominant layering mechanism of dust-content variation during accumulation rather than one of lag production during periods of sublimation.     

Interfacial liquid water on Mars and its potential role in formation of hill and dune gullies

Gullies are among the most intriguing structures identified on the surface of Mars. Most common are gullies located on the slopes of craters which are probably formed by liquid water transported by shallow aquifers (Heldmann, J.L., Carlsson, E., Johansson, H., Mellon, M.T., Toon, O.B. [2007]. Icarus 188, 324-344). Two particular types of gullies are found on slopes of isolated hills and dunes. The hill-slope gullies are located mostly at 50°S, which is at the high end of latitudes of bulk of the gullies found so far. The dune gullies are found in several locations up to 65°S (Reiss, D., Jaumann, R., Kereszturi, A., Sik, A., Neukum, G. [2007]. Lunar Planet. Sci. XXXVIII. Abstract 1993), but the best known are those in Russel crater at 54°S. The hill and dune gullies are longer than others making the aquifers explanation for their formation unlikely (Balme, M., Mangold, N., Baratoux, D., Costard, F., Gosselin, M., Masson, P., Pnet, P., Neukum, G. [2006]. J. Geophys. Res. 111. doi:10.1029/2005JE002607). Recently it has been noted that thin liquid films of interfacial water can play a role in rheological processes on the surface of Mars (Moehlmann, D. [2008]. Icarus 195, 131-139. Kereszturi, A., Moehlmann, D., Berczi, Sz., Ganti, T., Kuti, A., Sik, A., Horvath, A. [2009]. Icarus 201, 492-503.). Here we try to answer the question whether interfacial liquid water may occur on Mars in quantities large enough to play a role in formation of gullies. To verify this hypothesis we have calculated thermal models for hills and dunes of various steepness, orientation and physical properties. We find that within a range of average expected values of parameters it is not possible to have more than a few monolayers of liquid water at depths greater than a centimeter. To create subsurface interfacial water film significantly thicker and hence to produce conditions for the slope instability, parameters have to be chosen to have their extreme realistic values or an additional source of surface heating is needed. One possibility for additional heating is a change of atmospheric conditions due to a local dust storm. We conclude that if interfacial water is responsible for the formation of the hill-slope gullies, our results may explain why the hill gullies are rare.     

Digging deep for ice in Isidis Planitia—New constraints on subsurface volatile concentrations from thermal modeling

The Isidis Planitia region on Mars usually is regarded as a comparably attractive site for landing missions based on engineering constraints such as elevation and smooth regional topography. The Mars Express landed element Beagle 2 was deployed to this area, and the southern margin of the basin was selected as one of the backup landing sites for the NASA Mars Exploration Rovers.Especially in the context of the Beagle 2 mission, Isidis Planitia has been discussed as a place which might have experienced a volatile-rich history with associated potential for biological activity [e.g. Bridges et al., 2003. Selection of the landing site in Isidis Planitia of Mars Probe Beagle 2. J. Geophys. Res. 108(E1), 5001, doi: 10.1029/2001JE001820]. However the measurements of by the GRS instrument on Mars Odyssey indicate a maximum inferred water abundance of only 3 wt% in the upper few meters of the surface [Feldman et al., 2004. Global distribution of near-surface hydrogen on Mars. J. Geophys. Res. 109, E09006, doi: 10.1029/2003JE002160]. Based on these measurements this area seems to be one of the driest spots in the equatorial region of Mars.To support future landing site selections we took a more detailed look at the minimum burial depth of stable ice deposits in this area, focusing as an example on the planned Beagle 2 landing site. We are especially interested in the likelihood of ground ice deposits within the range of proposed subsurface sampling tools as drills or ‘mole’-like devices [Richter et al., 2002. Development and testing of subsurface sampling devices for the Beagle 2 Lander. Planet. Space Sci. 50, 903-913] given reasonable physical constraints for the surface and near surface material.For a mission like ExoMars [Kminek, G., Vago, J.L., 2005. The Aurora Exploration Program—The ExoMars Mission. In: Proceedings of the 35th Lunar and Planetary Science Conference, abstract no. 1111, 15-19 March 2004, League City, TX] with a focus on finding traces of fossil life the area might be of potential interest, because these traces would be better conserved in the dry soil. Modeling and measurement indicate that Isidis Planitia is indeed a dry place and any hypothetical ground ice deposits in this region are out of range of currently proposed sampling devices.     

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.     

A laboratory model of Saturn’s North Polar Hexagon

A hexagonal structure has been observed at ∼76°N on Saturn since the 1980s (Godfrey, D.A. [1988]. Icarus 76, 335-356). Recent images by Cassini (Baines, K., Momary, T., Roos-Serote, M., Atreya, S., Brown, R., Buratti, B., Clark, R., Nicholson, P. [2007]. Geophys. Res. Abstr. 9, 02109; Baines, K., Momary, T., Fletcher, L., Kim, J., Showman, A., Atreya, S., Brown, R., Buratti, B., Clark, R., Nicholson, P. [2009]. Geophys. Res. Abstr. 11, 3375) have shown that the feature is still visible and largely unchanged. Its long lifespan and geometry has puzzled the planetary physics community for many years and its origin remains unclear. The measured rotation rate of the hexagon may be very close to that of the interior of the planet (Godfrey, D.A. [1990]. Science 247, 1206-1208; Caldwell, J., Hua, X., Turgeon, B., Westphal, J.A., Barnet, C.D. [1993]. Science 206, 326-329; Sánchez-Lavega, A., Lecacheux, J., Colas, F., Laques, P. [1993]. Science 260, 329-332), leading to earlier interpretations of the pattern as a stationary planetary wave, continuously forced by a nearby vortex (Allison, M., Godfrey, D.A., Beebe, R.F. [1990]. Science 247, 1061-1063). Here we present an alternative explanation, based on an analysis of both spacecraft observations of Saturn and observations from laboratory experiments where the instability of quasi-geostrophic barotropic (vertically uniform) jets and shear layers is studied. We also present results from a barotropic linear instability analysis of the saturnian zonal wind profile, which are consistent with the presence of the hexagon in the North Pole and absence of its counter-part in the South Pole. We propose that Saturn’s long-lived polygonal structures correspond to wavemodes caused by the nonlinear equilibration of barotropically unstable zonal jets.     

Principal components analysis of Mars in the near-infrared

Principal components analysis and target transformation are applied to near-infrared image cubes of Mars in a study to disentangle the spectra into a small number of spectral endmembers and characterize the spectral information. The image cubes are ground-based telescopic data from the NASA Infrared Telescope Facility during the 1995 and 1999 near-aphelion oppositions when ice clouds were plentiful [ [Clancy et al., 1996] and [56]], and the 2003 near-perihelion opposition when ice clouds are generally limited to topographically high regions (volcano cap clouds) but airborne dust is more common [Martin, L.J., Zurek, R.W., 1993. J. Geophys. Res. 98 (E2), 3221-3246]. The heart of the technique is to transform the data into a vector space along the dimensions of greatest spectral variance and then choose endmembers based on these new “trait” dimensions. This is done through a target transformation technique, comparing linear combinations of the principal components to a mineral spectral library. In general Mars can be modeled, on the whole, with only three spectral endmembers which account for almost 99% of the data variance. This is similar to results in the thermal infrared with Mars Global Surveyor Thermal Emission Spectrometer data [Bandfield, J.L., Hamilton, V.E., Christensen, P.R., 2000. Science 287, 1626-1630]. The globally recovered surface endmembers can be used as inputs to radiative transfer modeling in order to measure ice abundance in martian clouds [Klassen, D.R., Bell III, J.F., 2002. Bull. Am. Astron. Soc. 34, 865] and a preliminary test of this technique is also presented.     

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.     

The effect of surface topography on the lunar photoelectron sheath and electrostatic dust transport

The dayside near-surface lunar plasma environment is electrostatically complex, due to the interaction between solar UV-induced photoemission, the collection of ambient ions and electrons, and the presence of micron and sub-micron sized dust grains. Further complicating this environment, although less well understood in effect, is the presence of surface relief, typically in the form of craters and/or boulders. It has been suggested that such non-trivial surface topography can lead to complex electrostatic potentials and fields, including “mini-wakes” behind small obstacles to the solar wind flow and “supercharging” near sunlit-shadowed boundaries (Criswell, D.R., De, B.R. [1977]. J. Geophys. Res. 82 (7); De, B.R., Criswell, D.R. [1977]. J. Geophys. Res. 82 (7); Farrell, W.M., Stubbs, T.J., Vondrak, R.R., Delory, G.T., Halekas, J.S. [2007]. Geophys. Res. Lett. 34; Wang, X., Horányi, M., Sternovsky, Z., Robertson, S., Morfill, G.E. [2007]. Geophys. Res. Lett. 34, L16104). In this paper, we present results from a three-dimensional, self-consistent, electrostatic particle-in-cell code used to model the dayside near-surface lunar plasma environment over a variety of local times with the presence of a crater. Additionally, we use the particle-in-cell model output to study the effect of surface topography on the dynamics of electrostatic dust transport, with the goal of understanding previous observations of dust dynamics on the Moon and dust ponding on various asteroids.     

An assessment and test of Enceladus as an important source of Saturn’s ring atmosphere and ionosphere

The existence of an oxygen exosphere and ionosphere in Saturn’s main ring region has been confirmed by the Saturn Orbital Insertion (SOI) observations of the Cassini spacecraft. Through the ion-molecule collisions, the ring atmosphere could serve as a source of ions throughout Saturn’s magnetosphere. If photolysis of ice in the main rings is the dominant source of O₂, then the complex structure of the ring atmosphere/ionosphere and the injection rate of neutral O₂ will be subject to modulation by the seasonal variation of Saturn along its orbit (Tseng, Wei-Ling, Ip, W.-H., Johnson, R.E., Cassidy, T.A., Erlod, M.K. [2010]. Icarus 206, 382-389). In addition, the radio and plasma wave science (RPWS) instrument onboard Cassini found that a large amount of the Enceladus-originated water-group plasma would be deposited on the outer edge of the A ring (Farrell, W.M., Kaiser, M.L., Gurnett, D.A., Kurth, W.S., Persoon, A.M., Wahlund, J.E., Canu, P. [2008]. Geophys. Res. Lett. 35, L02203). A large amount of Enceladus’ plume neutrals (water-group neutrals) would collide with the main rings through collisional interaction with the ambient neutrals and plasma ions (Jurac, S., Richardson, J.D. [2007]. Geophys. Res. Lett. 34, L08102; Cassidy, T.A., Johnson, R.E. [2010]. Icarus, in press). These absorbed ions and neutrals could be recycled to neutral oxygen molecules via grain-surface chemistry to contribute the ring oxygen atmosphere. In this work, we have examined the mass budget of the ring oxygen atmosphere of Saturn taking into account such an “exogenic” source. The maximum O₂ source rate from recycling of Enceladus-originated plasma and neutrals is probably comparable or higher to the one from photolytic decomposition of ices. In the above case, the neutral O₂ source rate would be independent of the solar insolation angle. Therefore, even at Saturn’s Equinox, the extended oxygen atmosphere still could be an important supplier of oxygen ions in the saturnian magnetosphere. We have performed several studies for different recycling source rates from Enceladus. These predictions need further the Cassini Plasma Spectrometer (CAPS) and the Magnetospheric Imaging Instrument (MIMI) observations to be verified in future.     

Surface composition and physical properties of several trans-neptunian objects from the Hapke scattering theory and Shkuratov model

Several different trans-neptunian objects have been studied in order to investigate their physical and chemical properties. New observations in the 1.1-1.4 μm range, obtained with the ISAAC instrument, are presented in order to complete previous observations carried out with FORS1 in the visible and SINFONI in the near infrared. All of the observations have been performed at the ESO/Very Large Telescope. We analyze the spectra of six different objects (2003 AZ₈₃, Echeclus, Ixion, 2002 AW₁₉₇, 1999 DE₉ and 2003 FY₁₂₈) in the 0.45-2.3 μm range with the model of Hapke (Hapke, B. [1981]. J. Geophys. Res. 86, 4571-4586) and the method of Shkuratov et al. (Shkuratov, Y., Starukhina, L., Hoffmann, H., Arnold, G. [1999]. Icarus 137, 235-246). Water ice is found on two objects, and in particular it is confirmed in its amorphous and crystalline states on 2003 AZ₈₄ surface. Upper limits on the water ice content are given for the other four TNOs investigated, confirming previous results (Barkume, K.M., Brown, M.E., Schaller, E.L. [2008]. Astron. J. 135, 55-67; Guilbert, A., Alvarez-Candal, A., Merlin, F., Barucci, M.A., Dumas, C., de Bergh, C., Delsanti, A. [2009]. Icarus 201, 272-283). Whatever the absorption features in the near infrared, all objects but one exhibit a moderate red slope in the visible, as most TNOs and Centaurs. We discuss the implications of the presence of water ice and the probable sources of the red slope.     

Experimental test of a radiative transfer model of the optical effects of space weathering

We compare laboratory measurements of the optical effects of nanophase iron on near-IR reflectance spectra of transparent silica gel infused with small iron particles [Noble, S.K., Pieters, C.M., Keller, L.P., 2007. Icarus 192, 629-642] with a radiative transfer model of the process [Hapke, B., 2001. J. Geophys. Res. 106 (E5), 10039-10074]. We find that the measurements exhibit reddening and darkening effects of nanophase (<50 nm) iron particles, a darkening effect of somewhat larger particles (>50 nm) and mixing effects of silica gel particles of varying total iron abundance. The radiative transfer model reproduces the effects of nanophase iron within the experimental uncertainties.     

Exponential Gaussian approach for spectral modeling: The EGO algorithm I. Band saturation

Curve fitting techniques are a widespread approach to spectral modeling in the VNIR range [Burns, R.G., 1970. Am. Mineral. 55, 1608-1632; Singer, R.B., 1981. J. Geophys. Res. 86, 7967-7982; Roush, T.L., Singer, R.B., 1986. J. Geophys. Res. 91, 10301-10308; Sunshine, J.M., Pieters, C.M., Pratt, S.F., 1990. J. Geophys. Res. 95, 6955-6966]. They have been successfully used to model reflectance spectra of powdered minerals and mixtures, natural rock samples and meteorites, and unknown remote spectra of the Moon, Mars and asteroids. Here, we test a new decomposition algorithm to model VNIR reflectance spectra and call it Exponential Gaussian Optimization (EGO). The EGO algorithm is derived from and complementary to the MGM of Sunshine et al. [Sunshine, J.M., Pieters, C.M., Pratt, S.F., 1990. J. Geophys. Res. 95, 6955-6966]. The general EGO equation has been especially designed to account for absorption bands affected by saturation and asymmetry. Here we present a special case of EGO and address it to model saturated electronic transition bands. Our main goals are: (1) to recognize and model band saturation in reflectance spectra; (2) to develop a basic approach for decomposition of rock spectra, where effects due to saturation are most prevalent; (3) to reduce the uncertainty related to quantitative estimation when band saturation is occurring. In order to accomplish these objectives, we simulate flat bands starting from pure Gaussians and test the EGO algorithm on those simulated spectra first. Then we test the EGO algorithm on a number of measurements acquired on powdered pyroxenes having different compositions and average grain size and binary mixtures of orthopyroxenes with barium sulfate. The main results arising from this study are: (1) EGO model is able to numerically account for the occurrence of saturation effects on reflectance spectra of powdered minerals and mixtures; (2) the systematic dilution of a strong absorber using a bright neutral material is not responsible for band deformation. Further work is still required in order to analyze the behavior of the EGO algorithm with respect to the saturation phenomena using more complex band shapes than pyroxene bands.