Publications

The Martian South Polar Layered Deposits (SPLD) are composed mostly of ice and dust with a thin perennial CO2 cover and some internal CO2 ice layers. In the North, the seasonal CO2 cap is lost during summer, allowing H2O ice to sublimate into the atmosphere. In the South, the perennial CO2 cover prevents H2O ice sublimation. This work uses the Mars Planetary Climate Model to investigate how the H2O and CO2 cycles are affected if the thin perennial CO2 SPLD cover is lost. We find that during southern summer, the atmospheric water content will more than double in the south polar region. However, on a global scale, the NPLD is still the dominant source of humidity because of its larger surface area. When exposing some of the South Polar Cap buried water ice, the south polar cap becomes the dominant source of atmospheric humidity due to Mars's spin-orbital alignment.

(Formula presented.) is the primary component of the martian atmosphere and its seasonal surface-atmosphere exchange is responsible for many of the climate phenomena on the planet. Near-surface ground water ice (GI) has been found to inhibit seasonal (Formula presented.) ice accumulations. Previous studies concerning the response of the (Formula presented.) cycle to orbital variations did not take into account the redistribution of GI arising from the same orbital variations. This work aims to analyze the effect of GI redistribution on the (Formula presented.) cycle in past climates. We use the LMD Planetary Climate Model to simulate the full (Formula presented.) cycle at different orbital configurations and compare simulations with reference modern GI as observed by the Mars Odyssey Neutron Spectrometer (\u201cMONS GI\u201d scenario) to simulations with equilibrium GI produced by the Mars Subsurface Ice Model (\u201cEq. GI\u201d scenario). In the Eq. GI scenario, equilibrium GI underlies 0.80.9 of the seasonal caps area at high obliquity periods, whereas in the reference MONS GI scenario, the overlap between GI and the seasonal cap is reduced, reaching less than 0.3 by obliquity (Formula presented.). The mass and duration of seasonal (Formula presented.) ice are significantly reduced relative to the reference scenario, especially in the mid-latitudes, and the expected increase of seasonal pressure amplitude with obliquity is attenuated (by a factor of (Formula presented.) 2 at obliquity (Formula presented.)). Interpolating the seasonal variations and mean annual pressure over the past 20,000 kyr, we highlight the influence of GI migration in attenuating both seasonal pressure variations and long-term oscillations of the martian atmospheric pressure.

We fit a dynamical model to Kepler systems that contain four or more transiting planets using the analytic method AnalyticLC and obtain physical and orbital parameters for 101 planets in 23 systems, of which 95 are of mass significance better than 3σ, and 46 are without previously reported mass constraints or upper limits. In addition, we compile a list of 71 Kepler objects of interest that display significant transit impact parameter variations (TbVs), complementing our previously published work on two- and three-transiting-planet systems. Together, these works include the detection of significant TbV signals of 130 planets, which is, to our knowledge, the largest catalog of this type to date. The results indicate that the typical detectable TbV rate in the Kepler population is of order 10−2 yr−1 and that rapid TbV rates (≳0.05 yr−1) are observed only in systems that contain a transiting planet with an orbital period less than ∼20 days. The observed TbV rates are only weakly correlated with orbital period within Keplers ≲100-day-period planets. If this extends to longer periods, it implies a limit on the utility of the transit technique for long-period planets. The TbVs we find may not be detectable in direct impact parameter measurements, but rather are inferred from the full dynamics of the system, encoded in all types of transit variations. Finally, we find evidence that the mutual inclination distribution is qualitatively consistent with the previously suggested angular momentum deficit model using an independent approach.

Estimation of planetary orbital and physical parameters from light-curve data relies heavily on the accurate interpretation of Transit Timing Variation (TTV) measurements. In this letter, we review the process of TTV measurement and compare two fitting paradigms-one that relies on making transit-by-transit timing estimates and then fitting a TTV model to the observed timings, and one that relies on fitting a global flux model to the entire light-curve data set simultaneously. The latter method is achieved either by solving for the underlying planetary motion (often referred to as "photodynamics"), or by using an approximate or empirical shape of the TTV signal. We show that, across a large range of the transit S/N regime, the probability distribution function of the mid-transit time significantly deviates from a Gaussian, even if the flux errors do distribute normally. Treating the timing uncertainties as if they are distributed normally leads, in such a case, to a wrong interpretation of the TTV measurements. We illustrate these points using numerical experiments and conclude that a fitting process that relies on a global flux fitting, rather than the derived TTVs, should be preferred.

Mars harbors ice deposits in several forms, on the surface and in the subsurface, which exchange with each other on various timescales. We seek to study the pore ice evolution over millennial time scales and how it contributes to and affects the Polar cap's evolution. We calculate the evolution of SubSurface Ice (SSI) pore filling by coupling two models, the Mars LMD Global Climate Model, which calculates the atmospheric and surface evolution on an annual timescale, and the dynamical version of the Mars Subsurface Ice Model, which calculates the evolution of the SSI on a millennial timescale. The SSI latitudinal boundary fluctuates over more than 25° in one obliquity cycle, overall extending equatorward of latitude ±35° at high obliquity, and receding to about ±60° at low obliquity. In locations where the SSI is stable continuously over orbital cycles, the simulations predict layering caused by a sublimation front at the SSI top boundary. Between 5 and 2.5 Myr ago, the subsurface lost at least ∼95 m of polar equivalent layer ice. The SSI flux routinely reaches ∼1 mm/Mars year. In addition to the direct contribution to the growth of the North Polar Layered Deposits (NPLD), the SSI causes variations in the NPLD accumulation rate due to the changes in the SSI distribution that affect the seasonal energy budget. These variations are comparable to the change in rate due to variations in orbital elements. When running paleo-climate simulations, particularly to reconstruct the NPLD profile, changes in the SSI distribution should be considered.

Exposed scarps images and ice-penetrating radar measurements in the North Polar Layered Deposits (NPLD) of Mars show alternating layers that provide an archive of past climate oscillations, that are thought to be linked to orbital variations, akin to Milankovitch cycles on Earth. We use the Laboratoire de Météorologie Dynamique Martian Global Climate Model to study paleoclimate states to enable a better interpretation of the NPLD physical and chemical stratigraphy. When a tropical ice reservoir is present, water vapor transport from the tropics to the poles at low obliquity is modulated by the intensity of summer. At times of low and relatively constant obliquity, the flux still varies due to other orbital elements, promoting polar layer formation. Ice migrates from the tropics toward the poles in two stages. First, when surface ice is present in the tropics, and second, when the equatorial deposit is exhausted, from ice that was previously deposited in mid-high latitudes. The polar accumulation rate is significantly higher when tropical ice is available, forming thicker layers per orbital cycle. However, the majority of the NPLD is sourced from ice that temporary resided in the mid-high latitudes and the layers become thinner as the source location moves poleward. The migration stages imprint different D/H ratios in different sections in the PLDs. The NPLD is isotopically depleted compared to the South Polar Layered Deposits in all simulations. Thus we predict the D/H ratio of the atmosphere in contact with NPLD upper layers is biased relative to the average global ice reservoirs.

We apply AnalyticLC, an analytic model described in an accompanying paper, to interpret Kepler data of systems that contain two or three transiting planets. We perform tests to verify that the obtained solutions agree with full N-body integrations, and that the number of model parameters is statistically justified. We probe non-co-planar interactions via impact parameter variations (TbVs), enabled by our analytic model. The subset of systems with a valid solution includes 54 systems composed of 140 planets, more than half of which are without previously reported mass constraints. Overall we provide: (i) estimates of physical and orbital properties for all systems analyzed; (ii) 102 planets with mass detections significant to better than 3σ, 43 of which have masses below 5 m ⊕; and (iii) 35 TbVs significant to better than 3σ. We focus on select systems showing strong TbVs, which can result from either interaction among the known transiting planets, or with a nontransiting object, and provide: (iv) a method to constrain the parameters of such unseen companions. These results are enabled by an accurate 3D photodynamical model, of a kind expected to become increasingly important for modeling multidecade photometric and composite (radial velocity, astrometry) data sets.

As human exploration missions to Mars are on the horizon, microbial cross-contamination remains a key issue to address. These issues can be approached today using advances in molecular metagenomics methods, which include rapid and sensitive sequencing platforms for characterizing microbial populations. Combined with analog missions, these methods provide powerful tools for assessing the challenges associated with planetary exploration. Here, we designed a protocol to monitor forward and backward contamination events and progression in an 11-days Mars analog mission in the Ramon crater in Israel. Forward contamination soil samples were collected daily from three sitestwo sites in close proximity to the habitat and one isolated site. Backward contamination was determined in samples from nitrile gloves of six analog astronauts before and after extravehicular activities Temperature, relative humidity and soil composition data were also collected for all sites. Environmental DNA samples were extracted in the main habitat and 16S (bacterial) and 18S (eukaryotic, fungal) rRNA gene amplicons were sequenced and analyzed to study microbial population diversity and composition. Shannon Diversity index analysis and Principal Coordinates analysis (PCoA) of rRNA genes indicated that differences in the diversity and population composition were significant in sites closer to the habitat when compared to a reference site. These samples also demonstrated the introduction of human-associated taxa to the environment. Backward contamination consisted of bacterial taxa found on gloves upon return from EVA and also detected in soil, altogether 44 genera, indicating backward contamination events. To our knowledge, this is the first protocol to utilize advanced molecular technologies to investigate forward and backward contamination in a Mars analog mission.

Mars Polar Science is an integrated, compelling system that serves as a nearby analogue to numerous other planets, supports human exploration, and habitability. Mars possesses the closest and most easily accessible layered ice deposits outside of Earth, and accessing those layers to read the climate record would be a triumph for planetary science.

Water ice is thought to be trapped in large permanently shadowed regions in the Moons polar regions, due to their extremely low temperatures. Here, we show that many unmapped cold traps exist on small spatial scales, substantially augmenting the areas where ice may accumulate. Using theoretical models and data from the Lunar Reconnaissance Orbiter, we estimate the contribution of shadows on scales from 1 km to 1 cm, the smallest distance over which we find cold-trapping to be effective for water ice. Approximately 1020% of the permanent cold-trap area for water is found to be contained in these micro cold traps, which are the most numerous cold traps on the Moon. Consideration of all spatial scales therefore substantially increases the number of cold traps over previous estimates, for a total area of ~40,000 km², about 60% of which is in the south. A majority of cold traps for water ice is found at latitudes > 80° because permanent shadows equatorward of 80° are typically too warm to support ice accumulation. Our results suggest that water trapped at the lunar poles may be more widely distributed and accessible as a resource for future missions than previously thought.

Non-Keplerian dynamics of planetary orbits manifest in the transit light curve as variations of different types. In addition to transit timing variations, the shape of the transits contains additional information on variations in the geometry of the orbit. This study presents an analytic approach to light-curve fitting: dynamical variations in the orbital elements are transformed to a light curve using an analytic function with a restricted set of fitting parameters. Our method requires no N-body integration, resulting in a smaller number of degrees of freedom and a faster calculation. The approach described here is for the case of secular perturbations. By assuming that the orbital motion is dominated by nodal and apsidal precessions, analytic expressions for the light-curve transit parameters are derived as a function of the orbital variations. Detecting and characterizing such dynamical scenarios provides information regarding the possible existence of nontransiting companions, or the nonspherical mass distribution of the host star. The variations may imply forces out of the orbital plane, and thus probe mutual inclinations among components of the system. The derived models successfully reproduce the vanishing transit signals of KOI 120.01, and suggest a possible interesting scenario of a planet orbiting one member of a close-in binary system undergoing unusually rapid nodal regression. The model parameters are degenerate, so we provide relevant information for follow-up observations, which are suggested in order to place further constraints on this unique Kepler object.

Asteroid pairs are those found to share similar heliocentric orbits but not positions. The leading theory suggests that each had a single progenitor that split due to rotational-fission of a weak, rubble-pile structured body. By constructing shape models of asteroid pairs from multiple-apparition observations and using a lightcurve inversion technique, we mapped the gravitational and rotational accelerations on the surfaces of these asteroids. This allows us to construct a map of local slopes on the asteroids' surfaces. In order to test for frictional failure, we determine the maximum rotation rate at which an area larger than half the surface area of the secondary member (assumed to be the ejected component) has a slope value >40 degrees, the angle of friction of lunar regolith, where loose material will begin sliding. We use this criterion to constrain the failure stress operating on the body, just before disruption at the commonly observed spin barrier of 2.2 h. Our current sample includes shape models of eleven primary members of asteroid pairs, observed from the Wise Observatory in the last decade. In the studied parameter space we find that the shape models only reach the spin barrier when their bulk density is larger than the ~2 g cm
^−
3 measured for the rubble pile structured 25143 Itokawa, and better matches 433 Eros' value of 2.7 g cm
^−
3, suggesting that km-sized asteroid pairs are dense compared to sub-km bodies. Assuming ejection of secondary components that are larger than those observed (up to the maximal size allowing separation), can also increase the spin barrier of the asteroids, thus supporting the previously suggested scenario of continuous disruption of the secondary. In addition, cohesion levels of hundreds of Pascals are also required to prevent these shape models from disrupting at spin rates slower than the usual spin barrier.

The layered polar caps of Mars have long been thought to be related to variations in orbit and axial tilt. We dynamically link Mars's past climate variations with the stratigraphy and isotopic composition of its ice by modeling the exchange of H2O and HDO among three reservoirs. The model shows that the interplay among equatorial, mid-latitude, and north-polar layered deposits (NPLD) induces significant isotopic changes in the cap. The diffusive properties of the sublimation lags and dust content in our model result in a cap size consistent with current Mars. The layer thicknesses are mostly controlled by obliquity variations, but the precession period of 50 kyr dominates the variations in the isotopic composition during epochs of relatively low and nearly constant obliquity such as at present. Isotopic sampling of the top 100 m may reveal climate oscillations unseen in the layer thicknesses and would thus probe recent precession-driven climate cycles.

The Earth-Moon system is suggested to have formed through a single giant collision, in which the Moon accreted from the impact-generated debris disk. However, such giant impacts are rare, and during its evolution, the Earth experienced many more smaller impacts, producing smaller satellites that potentially coevolved. In the multiple-impact hypothesis of lunar formation, the current Moon was produced from the mergers of several smaller satellites (moonlets), each formed from debris disks produced by successive large impacts. In the Myr between impacts, a pre-existing moonlet tidally evolves outward until a subsequent impact forms a new moonlet, at which point both moonlets will tidally evolve until a merger or system disruption. In this work, we examine the likelihood that preexisting moonlets survive subsequent impact events, and explore the dynamics of Earth-moonlet systems that contain two moonlets generated Myr apart. We demonstrate that pre-existing moonlets can tidally migrate outward, remain stable during subsequent impacts, and later merge with newly created moonlets (or re-collide with the Earth). Formation of the Moon from the mergers of several moonlets could therefore be a natural byproduct of the Earth's growth through multiple impacts. More generally, we examine the likelihood and consequences of Earth having prior moons, and find that the stability of moonlets against disruption by subsequent impacts implies that several large impacts could post-date Moon formation.

Turbet and Tran (2017) have identified that we mistakenly used air-broadened CO2 absorption spectra to generate coefficients for a rapid radiative transfer code, instead of self-broadened spectra. The resulting underestimation of absorption by CO2 led us to suggest that when the effects of CO2 line mixing are taken into account, surface temperatures on early Mars are up to 15K colder than those calculated in previous studies, in which empirical corrections of the CO2 spectrum were used ( factors). Using the correct CO2 absorption spectra, including the effects of line mixing, Turbet and Tran (2017) found that the surface temperature is colder by only 2K than temperatures calculated with the factor approach. While we acknowledge the mistake, and thank Turbet and Tran (2017) for finding and correcting it, we note that our main conclusions hold that surface temperatures are overestimated when empirical factors are used instead of a full account of CO2 line mixing in radiative transfer calculations in CO2-rich planetary atmospheres. Additionally, we emphasize that the rapid radiative transfer model we developed and tested is useful for studies of planetary climate, provided the correct absorption spectra are used to generate new k coefficients.

The heat flux incident upon the surface of an airless planetary body is dominated by solar radiation during the day, and by thermal emission from topography at night. Motivated by the close relationship between this heat flux, the surface temperatures, and the stability of volatiles, we consider the effect of the slope distribution on the temperature distribution and hence prevalence of cold-traps, where volatiles may accumulate over geologic time. We develop a thermophysical model accounting for insolation, reflected and emitted radiation, and subsurface conduction, and use it to examine several idealized representations of rough topography. We show how subsurface conduction alters the temperature distribution of bowl-shaped craters compared to predictions given by past analytic models. We model the dependence of cold-traps on crater geometry and quantify the effect that while deeper depressions cast more persistent shadows, they are often too warm to trap water ice due to the smaller sky fraction and increased reflected and reemitted radiation from the walls. In order to calculate the temperature distribution outside craters, we consider rough random surfaces with a Gaussian slope distribution. Using their derived temperatures and additional volatile stability models, we estimate the potential area fraction of stable water ice on Earth's Moon. For example, surfaces with slope RMS ∼15° (corresponding to length-scales ∼10 m on the lunar surface) located near the poles are found to have a ∼10% exposed cold-trap area fraction. In the subsurface, the diffusion barrier created by the overlaying regolith increases this area fraction to ∼40%. Additionally, some buried water ice is shown to remain stable even beneath temporarily illuminated slopes, making it more readily accessible to future lunar excavation missions. Finally, due to the exponential dependence of stability of ice on temperature, we are able to constrain the maximum thickness of the unstable layer to a few decimeters.

We find that the reflectance of the lunar surface within 5 degrees of latitude of the South Pole increases rapidly with decreasing temperature, near similar to 110 K, behavior consistent with the presence of surface water ice. The North polar region does not show this behavior, nor do South polar surfaces at latitudes more than 5 degrees from the pole. This South pole reflectance anomaly persists when analysis is limited to surfaces with slopes less than 10 degrees to eliminate false detection due to the brightening effect of mass wasting, and also when the very bright south polar crater Shackleton is excluded from the analysis. We also find that south polar regions of permanent shadow that have been reported to be generally brighter at 1064 nm do not show anomalous reflectance when their annual maximum surface temperatures are too high to preserve water ice. This distinction is not observed at the North Pole. The reflectance excursion on surfaces with maximum temperatures below 110 K is superimposed on a general trend of increasing reflectance with decreasing maximum temperature that is present throughout the polar regions in the north and south; we attribute this trend to a temperature or illumination-dependent space weathering effect (e.g. Hemingway et al., 2015). We also find a sudden increase in reflectance with decreasing temperature superimposed on the general trend at 200 K and possibly at 300 K. This may indicate the presence of other volatiles such as sulfur or organics. We identified and mapped surfaces with reflectances so high as to be unlikely to be part of an ice-free population. In this south we find a similar distribution found by Hayne et al. (2015) based on UV properties. In the north a cluster of pixels near that pole may represent a limited frost exposure. (C) 2017 Elsevier Inc. All rights reserved.

In June 2009 the Lunar Reconnaissance Orbiter (LRO) spacecraft was launched to the Moon. The payload consists of 7 science instruments selected to characterize sites for future robotic and human missions. Among them, the Lunar Orbiter Laser Altimeter (LOLA) was designed to obtain altimetry, surface roughness, and reflectance measurements. The primary phase of lunar exploration lasted one year, following a 3-month commissioning phase. On completion of its exploration objectives, the LRO mission transitioned to a science mission. After 7 years in lunar orbit, the LOLA instrument continues to map the lunar surface. The LOLA dataset is one of the foundational datasets acquired by the various LRO instruments. LOLA provided a high-accuracy global geodetic reference frame to which past, present and future lunar observations can be referenced. It also obtained high-resolution and accurate global topography that were used to determine regions in permanent shadow at the lunar poles. LOLA further contributed to the study of polar volatiles through its unique measurement of surface brightness at zero phase, which revealed anomalies in several polar craters that may indicate the presence of water ice. In this paper, we describe the many LOLA accomplishments to date and its contribution to lunar and planetary science.

The rubble pile spin barrier is an upper limit on the rotation rate of asteroids larger than ~200-300. m. Among thousands of asteroids with diameters larger than ~300. m, only a handful of asteroids are known to rotate faster than 2.0. h, all are in the sub-km range (≤0.6. km). Here we present photometric measurements suggesting that (60716) 2000 GD65, an S-complex, inner-main belt asteroid with a relatively large diameter of 2.3-0.7+0.6km, completes one rotation in 1.9529. ±. 0.0002. h. Its unique diameter and rotation period allow us to examine scenarios about asteroid internal structure and evolution: a rubble pile bound only by gravity; a rubble-pile with strong cohesion; a monolithic structure; an asteroid experiencing mass shedding; an asteroid experiencing YORP spin-up/down; and an asteroid with a unique octahedron shape results with a four-peak lightcurve and a 3.9. h period. We find that the most likely scenario includes a lunar-like cohesion that can prevent (60716) 2000 GD65 from disrupting without requiring a monolithic structure or a unique shape. Due to the uniqueness of (60716) 2000 GD65, we suggest that most asteroids typically have smaller cohesion than that of lunar regolith.

Impact cratering produces characteristic variations in the topographic power spectral density (PSD) of cratered terrains, which are controlled by the size-frequency distribution of craters and the spectral content (shape) of individual features. These variations are investigated here in two parallel approaches. First, a cratered terrain model, based on Monte Carlo emplacement of craters and benchmarked by an analytical formulation of the one-dimensional PSD, is employed to generate topographic surfaces at a range of size-frequency power law exponents and shape dependencies. For self-similar craters, the slope of the PSD, , varies inversely with that of the production function, , leveling off to 0 at high (surface topography dominated by the smallest craters) and maintaining a roughly constant value ( approximate to 2) at low (surface topography dominated by the largest craters). The effects of size-dependent shape parameters and various crater emplacement rules are also considered. Second, we compare the model-derived predictions for the behavior of the PSD with values of calculated along transects from the Lunar Orbiter Laser Altimeter (LOLA). At small scales (approximate to 115 m to 1 km) model predictions agree well with the PSD slope over the observed range of lunar size-frequency distributions. Differences between global PSD slopes at subkilometer and kilometer scales reflect a scale separation in roughness consistent with prior observations using a variety of surface roughness parameters. Understanding the statistical markers left by the impact cratering process on the lunar surface is useful for distinguishing between competing geological processes on planetary surfaces throughout the solar system.

The heat flux experienced by the surface of an airless planetary body is dominated by solar insolation during the day, and by topography at night. Motivated by the close relationship between this heat flux, surface temperature and volatile stability, we consider the effect of the slope distribution on the temperature distribution and hence, prevalence of cold-traps where volatiles may accumulate over geologic time. We present a thermophysical model, accounting for four effects: insolation, reflected and emitted radiation from neighbouring slopes, and subsurface conduction. We investigate specific representative geometries, as well as generic topographies defined by their slope distributions. We find that depressions, such as craters, cast shadows that are more persistent, but warmer than those transiently cast by positive relief features convex in map view, such as hills. We find that while scarce, some briefly illuminated areas may still serve as cold-traps, potentially allowing them to be imaged from orbit. Finally, we estimate the amount of surface and subsurface ice beneath topographies with different roughness and discuss its dependency on latitude.

The Cassini Synthetic Aperture Radar has been acquiring images of Titan's surface since October 2004. To date, 59% of Titan's surface has been imaged by radar, with significant regions imaged more than once. Radar data suffer from speckle noise hindering interpretation of small-scale features and comparison of reimaged regions for change detection. We present here a new image analysis technique that combines a denoising algorithm with mapping and quantitative measurements that greatly enhance the utility of the data and offers previously unattainable insights. After validating the technique, we demonstrate the potential improvement in understanding of surface processes on Titan and defining global mapping units, focusing on specific landforms including lakes, dunes, mountains, and fluvial features. Lake shorelines are delineated with greater accuracy. Previously unrecognized dissection by fluvial channels emerges beneath shallow methane cover. Dune wavelengths and interdune extents are more precisely measured. A significant refinement in producing digital elevation models is shown. Interactions of fluvial and aeolian processes with topographic relief is more precisely observed and understood than previously. Benches in bathymetry are observed in northern sea Ligeia Mare. Submerged valleys show similar depth suggesting that they are equilibrated with marine benches. These new observations suggest a liquid level increase in the northern sea, which may be due to changes on seasonal or longer timescales.

It has long been suggested that water ice can exist in extremely cold regions near the lunar poles, where sublimation loss is negligible. The geographic distribution of H-bearing regolith shows only a partial or ambiguous correlation with permanently shadowed areas, thus suggesting that another mechanism may contribute to locally enhancing water concentrations. We show that under suitable conditions, water molecules can be pumped down into the regolith by day-night temperature cycles, leading to an enrichment of H₂O in excess of the surface concentration. Ideal conditions for pumping are estimated and found to occur where the mean surface temperature is below 105 K and the peak surface temperature is above 120 K. These conditions complement those of the classical cold traps that are roughly defined by peak temperatures lower than 120 K. On the present-day Moon, an estimated 0.8% of the global surface area experiences such temperature variations. Typically, pumping occurs on pole-facing slopes in small areas, but within a few degrees of each pole the equator-facing slopes are preferred. Although pumping of water molecules is expected over cumulatively large areas, the absolute yield of this pump is low; at best, a few percent of the H₂O delivered to the surface could have accumulated in the near-surface layer in this way. The amount of ice increases with vapor diffusivity and is thus higher in the regolith with large pore spaces.

We model the primary crater production of small (D20. km higher than Zunil), a factor ~2 younger than estimated using the Hartmann production function which assumes 6. mbar surface pressure. We estimate ages of 52.3. Ma and 23.9. Ma for North Ray and Cone crater respectively, consistent with cosmic ray exposure ages from Apollo samples. Our results indicate that the average cratering rate has been constant on these bodies over these time periods. Since our Monte Carlo simulations demonstrate that the existing crater chronology systems can be applied to date young surfaces using small craters on the Moon and Mars, we conclude that the signal from secondary craters in the isochrons must be relatively small at these locations, as our Monte Carlo model only generates primary craters.

The decay of the martian atmosphere - which is dominated by carbon dioxide - is a component of the long-term environmental change on Mars from a climate that once allowed rivers to flow to the cold and dry conditions of today. The minimum size of craters serves as a proxy for palaeopressure of planetary atmospheres, because thinner atmospheres permit smaller objects to reach the surface at high velocities and form craters. The Aeolis Dorsa region near Gale crater on Mars contains a high density of preserved ancient craters interbedded with river deposits and thus can provide constraints on atmospheric density at the time of fluvial activity. Here we use high-resolution images and digital terrain models from the Mars Reconnaissance Orbiter to identify ancient craters in deposits in Aeolis Dorsa that date to about 3.6 Gyr ago and compare their size distribution with models of atmospheric filtering of impactors. We obtain an upper limit of 0.9 ± 0.1 bar for the martian atmospheric palaeopressure, rising to 1.9 ± 0.2 bar if rimmed circular mesas - interpreted to be erosionally-resistant fills or floors of impact craters - are excluded. We assume target properties appropriate for desert alluvium: if sediment had rock-mass strength similar to bedrock at the time of impact, the paleopressure upper limit increases by a factor of up to two. If Mars did not have a stable multibar atmosphere at the time that the rivers were flowing - as suggested by our results - then a warm and wet CO₂/H₂O greenhouse is ruled out, and long-term average temperatures were most likely below freezing.

The presence of an atmosphere, initially suggested based on limb darkening by Sola (1904) and later by the presence of methane spectral lines by Kuiper (1944), has long given Titan a special place in the minds of planetary geologists. The first close-up images were obtained by Pioneer 11 in 1979 (Gehrels et al., 1980), confirming a substantial atmosphere. These early observations led to the diversion of the trajectory of the Voyager I spacecraft to a closer encounter with Titan in 1980. Although the visible cameras on Voyager also had difficulty seeing Titan's surface (Richardson et al., 2004), radio occultation experiments suggested a surface pressure of 1.5 bars and temperature near 95 K (Lindal et al., 1983). These results were exciting because, for a methane mixing ratio of a few percent at the surface (Hunten, 1978), they placed methane's partial pressure near its triple point. Thus, like water on Earth, solid, liquid, and gaseous methane could potentially exist in Titan's environment. Ethane, which is the main product of methane photolysis, can also be liquid under these conditions. The presence of condensable volatiles in Titan's thick atmosphere opens the door for active fluvial, lacustrine, and pluvial processes that can shape its landscape with similar morphologies to those we find on Earth.

Cometary activity in main-belt asteroids probes the ice content of these objects and provides clues to the history of volatiles in the inner Solar system. We search the Palomar Transient Factory survey to derive upper limits on the population size of activemain-belt comets (MBCs). From data collected from 2009 March through 2012 July, we extracted~2million observations of ~220 thousand known main-belt objects (40 per cent of the known population, down to ~1-km diameter) and discovered 626 new objects in multinight linked detections. We formally quantify the 'extendedness' of a small-body observation, account for systematic variation in this metric (e.g. due to on-sky motion) and evaluate this method's robustness in identifying cometary activity using observations of 115 comets, including two known candidate MBCs and six newly discovered non-MBCs (two of which were originally designated as asteroids by other surveys). We demonstrate a 66 per cent detection efficiency with respect to the extendedness distribution of the 115 sampled comets, and a 100 per cent detection efficiency with respect to extendedness levels greater than or equal to those we observed in the known candidateMBCs P/2010 R2 (La Sagra) and P/2006VW₁₃₉.Using a log-constant prior, we infer 95 per cent confidence upper limits of 33 and 22 activeMBCs (per million main-belt asteroids down to ~1-km diameter), for detection efficiencies of 66 and 100 per cent, respectively. In a follow-up to this morphological search, we will perform a photometric (disc-integrated brightening) search for MBCs.

Saturn's moon Titan has lakes and seas of liquid hydrocarbon and a dense atmosphere, an environment conducive to generating wind waves. Cassini observations thus far, however, show no indication of waves. We apply models for wind wave generation and detection to the Titan environment. Results suggest wind speed thresholds at a reference altitude of 10. m of 0.4-0.7. m/s for liquid compositions varying between pure methane and equilibrium mixtures with the atmosphere (ethane has a threshold of 0.6. m/s), varying primarily with liquid viscosity. This reduced threshold, as compared to Earth, results from Titan's increased atmosphere-to-liquid density ratio, reduced gravity and lower surface tension. General Circulation Models (GCMs) predict wind speeds below derived thresholds near equinox, when available observations of lake surfaces have been acquired. Predicted increases in winds as Titan approaches summer solstice, however, will exceed expected thresholds and may provide constraints on lake composition and/or GCM accuracy through the presence or absence of waves during the Cassini Solstice Mission. A two-scale microwave backscatter model suggests that returns from wave-modified liquid hydrocarbon surfaces may be below the pixel-scale noise floor of Cassini radar images, but can be detectable using real-aperture scatterometry, pixel binning and/or observations obtained in a specular geometry.

Shackleton crater is nearly coincident with the Moons south pole. Its interior receives almost no direct sunlight and is a perennial cold trap, making Shackleton a promising candidate location in which to seek sequestered volatiles. However, previous orbital and Earth-based radar mapping and orbital optical imaging have yielded conflicting interpretations about the existence of volatiles. Here we present observations from the Lunar Orbiter Laser Altimeter on board the Lunar Reconnaissance Orbiter, revealing Shackleton to be an ancient, unusually well-preserved simple crater whose interior walls are fresher than its floor and rim. Shackleton floor deposits are nearly the same age as the rim, suggesting that little floor deposition has occurred since the crater formed more than three billion years ago. At a wavelength of 1,064 nanometres, the floor of Shackleton is brighter than the surrounding terrain and the interiors of nearby craters, but not as bright as the interior walls. The combined observations are explicable primarily by downslope movement of regolith on the walls exposing fresher underlying material. The relatively brighter crater floor is most simply explained by decreased space weathering due to shadowing, but a one-micrometre-thick layer containing about 20 per cent surficial ice is an alternative possibility.

The Palomar Transient Factory (PTF) is a synoptic survey designed to explore the transient and variable sky in a wide variety of cadences. We use PTF observations of fields that were observed multiple times (≳10) per night, for several nights, to find asteroids, construct their light curves and measure their rotation periods. Here we describe the pipeline we use to achieve these goals and present the results from the first four (overlapping) PTF fields analysed as part of this programme. These fields, which cover an area of 21deg ², were observed on four nights with a cadence of ∼20min. Our pipeline was able to detect 624 asteroids, of which 145 (≈20 per cent) were previously unknown. We present high-quality rotation periods for 88 main-belt asteroids and possible period or lower limit on the period for an additional 85 asteroids. For the remaining 451 asteroids, we present lower limits on their photometric amplitudes. Three of the asteroids have light curves that are characteristic of binary asteroids. We estimate that implementing our search for all existing high-cadence PTF data will provide rotation periods for about 10000 asteroids mainly in the magnitude range ≈14 to ≈20.

The Mars polar layered deposits (PLD) likely hold an extensive record of recent climate during a period of high-amplitude orbit and obliquity cycles. Previous work has detected limited evidence for orbital signatures within PLD stratigraphy, but data from the High Resolution Imaging Science Experiment (HiRISE) permit renewed analysis of PLD stratigraphy at sub-meter scale. Topography derived from HiRISE images using stereogrammetry resolves beds previously detectable only as alternating light and dark bands in visible images. We utilize these data to measure the thickness of individual beds within the PLD, corrected for non-horizontal bed orientation. Stratigraphic columns and bed thickness profiles are presented for two sites within the NPLD, and show several sets of finely bedded units 1-2 m thick; isolated marker beds 3-4 m thick; and undifferentiated sections. Bed thickness measurements for three sites within the SPLD exhibit only one bed type based on albedo and morphology, and bed thicknesses have a larger mean and variance compared to measurements for the NPLD. Power spectra of brightness and slope derived along the measured stratigraphic sections confirm the regularity of NPLD fine bed thickness, and the lack of a dominant SPLD bed thickness. The regularity of fine bed thickness of the NPLD is consistent with quasiperiodic bed formation, albeit with unknown temporal period; the SPLD thickness measurements show no such regularity.

Mercury's coupled 3:2 spin-orbit resonance in conjunction with its relatively high eccentricity of ∼0.2 and near-zero obliquity results in both a latitudinal and longitudinal variation in annual average solar insolation and thus equatorial hot and cold regions. This results in an asymmetric temperature distribution in the lithosphere and a long wavelength lateral variation in lithosphere structure and strength that mirrors the insolation pattern. We employ a thermal evolution model for Mercury generating strength envelopes of the lithosphere to demonstrate and quantify the possible effects the insolation pattern has on Mercury's lithosphere. We find the heterogeneity in lithosphere strength is substantial and increases with time. We also find that a crust thicker than that of the Moon or Mars and dry rheologies for the crust and mantle are favorable when compared with estimates of brittle-ductile transition depths derived from lobate scarps. Regions of stronger and weaker compressive strength imply that the accommodation of radial contraction of Mercury as its interior cooled, manifest as lobate scarps, may not be isotropic, imparting a preferential orientation and distribution to the lobate scarps.

Cassini RADAR topography data are used to evaluate Titan's hypsometric profile, and to make comparisons with other planetary bodies. Titan's hypsogram is unimodal and strikingly narrow compared with the terrestrial planets. To investigate topographic extremes, a novel variant on the classic hypsogram is introduced, with a logarithmic abscissa to highlight mountainous terrain. In such a plot, the top of the terrestrial hypsogram is quite distinct from those of Mars and Venus due to the 'glacial buzz-saw' that clips terrestrial topography above the snowline. In contrast to the positive skew seen in other hypsograms, with a long tail of positive relief due to mountains, there is an indication (weak, given the limited data for Titan so far) that the Titan hypsogram appears slightly negatively skewed, suggesting a significant population of unfilled depressions. Limited data permit only a simplistic comparison of Titan topography with other icy satellites but we find that the standard deviation of terrain height (albeit at different scales) is similar to those of Ganymede and Europa.

The accuracy of geologic chronologies can, in principle, be improved through orbital tuning, the systematic adjustment of a chronology to bring the associated record into greater alignment with an orbitally derived signal. It would be useful to have a general test for the success of orbital tuning, and one proposal has been that eccentricity ought to covary with the amplitude envelope associated with precession variability recorded in tuned geologic records. A common procedure is to filter a tuned geologic record so as to pass precession period variability and compare the amplitude modulation of the resulting signal against eccentricity. There is a reasonable expectation for such a relationship to be found in paleoclimate records because the amplitude of precession forcing depends upon eccentricity. However, there also exists a relationship between eccentricity and the frequency of precession such that orbital tuning generates eccentricity-like amplitude modulation in filtered signals, regardless of the accuracy of the chronology or the actual presence of precession. This relationship results from the celestial mechanics governing eccentricity and precession and from the interaction between frequency modulation and amplitude modulation caused by filtering. When the eccentricity of Earth's orbit is small, the frequency of climatic precession undergoes large variations and less precession energy is passed through a narrow-band filter. Furthermore, eccentricity-like amplitude modulation is routinely obtained from pure noise records that are orbitally tuned to precession and then filtered. We conclude that the presence of eccentricity-like amplitude modulation in precession-filtered records does not support the accuracy of orbitally tuned time scales.

Of more than 400 filled lakes now identified on Titan, the first and largest reported in the southern latitudes is Ontario Lacus, which is dark in both infrared and microwave. Here we describe recent observations including synthetic aperture radar (SAR) images by Cassini's radar instrument (λ= 2 cm) and show morphological evidence for active material transport and erosion. Ontario Lacus lies in a shallow depression, with greater relief on the southwestern shore and a gently sloping, possibly wave-generated beach to the northeast. The lake has a closed internal drainage system fed by Earth-like rivers, deltas and alluvial fans. Evidence for active shoreline processes, including the wave-modified lakefront and deltaic deposition, indicates that Ontario is a dynamic feature undergoing typical terrestrial forms of littoral modification.

Ontario Lacus is the largest and best characterized lake in Titan's south polar region. In June and July 2009, the Cassini RADAR acquired its first Synthetic Aperture Radar (SAR) images of the area. Together with closest approach altimetry acquired in December 2008, these observations provide a unique opportunity to study the lake's nearshore bathymetry and complex refractive properties. Average radar backscatter is observed to decrease exponentially with distance from the local shoreline. This behavior is consistent with attenuation through a deepening layer of liquid and, if local topography is known, can be used to derive absorptive dielectric properties. Accordingly, we estimate nearshore topography from a radar altimetry profile that intersects the shoreline on the East and West sides of the lake. We then analyze SAR backscatter in these regions to determine the imaginary component of the liquid's complex index of refraction (κ). The derived value, κ = (6.1_-1.3^+1.7) × 10^-4, corresponds to a loss tangent of tan Δ = (9.2_-2.0^+2.5) × 10^-4 and is consistent with a composition dominated by liquid hydrocarbons. This value can be used to test compositional models once the microwave optical properties of candidate materials have been measured. In areas that do not intersect altimetry profiles, relative slopes can be calculated assuming the index of refraction is constant throughout the liquid. Accordingly, we construct a coarse bathymetry map for the nearshore region by measuring bathymetric slopes for eleven additional areas around the lake. These slopes vary by a factor of ∼5 and correlate well with observed shoreline morphologies.

The East Candor Interior Layered Deposit (ILD) has signatures of mono- and polyhydrated sulfate in alternating layers that give insight into the processes which formed these layered deposits and on the environmental conditions acting on them since then. We use orbital data to explore multiple hypotheses for how these deposits formed: (1) sulfate-bearing ILDs experience hydration changes on seasonal to a few years timescales under current Mars environmental conditions; (2) the deposits experience hydration under recent Mars conditions but require the wetter climate of high obliquity; and (3) the kieserite could be an original or diagenetic part of a complex evaporite mineral assemblage. Modeled climatology shows recent Mars environmental conditions might pass between multiple sulfate fields. However, comparison of Observatoire pour la Minéralogie, l'Eau, les Glaces et l'Activité (OMEGA) and Compact Reconnaissance Imaging Spectrometer (CRISM) observations of the same ILD do not show changes in hydration over 2 Mars years. Low temperatures might slow the kinetics of that transition; it is likely that more clement conditions during periods of high obliquity are needed to overcome mineral metastability and hydrate kieserite-bearing deposits. We find the alternate model, that the deposit is a cyclic evaporite sequence of mono- and polyhydrated sulfates, also plausible but with an unexplained dearth of Fe sulfates.

Experiments demonstrate for the first time the deposition of subsurface ice directly from atmospheric water vapor under Mars surface conditions. Deposition occurs at pressures below the triple point of water and therefore in the absence of a bulk liquid phase. Significant quantities of ice are observed to deposit in porous medium interstices; the maximum filling fraction observed in our experiments is ∼90%, but the evidence is consistent with ice density in pore spaces asymptotically approaching 100% filling. The micromorphology of the deposited ice reveals several noteworthy characteristics including preferential early deposition at grain contact points, complete pore filling between some grains, and captured atmospheric gas bubbles. The boundary between ice-bearing and ice-free soil, the "ice table," is a sharp interface, consistent with theoretical investigations of subsurface ice dynamics. Changes of surface radiative properties are shown to affect ice table morphology through their modulation of the local temperature profile. Accumulation of ice is shown to reduce the diffusive flux and thus reduce the rate of further ice deposition. Numerical models of the experiments based on diffusion physics are able to reproduce experimental ice contents if the parameterization of growth rate reduction has expected contributions in addition to reduced porosity. Several phenomena related to the evolution of subsurface ice are discussed in light of these results, and interpretations are given for a range of potential observations being made by the Phoenix Scout Lander.

Home Plate is a layered plateau observed by the Mars Exploration Rover Spirit in the Columbia Hills of Gusev Crater. The structure is roughly 80 m in diameter, and the raised margin exposes a stratigraphic section roughly 1.5 m in thickness. Previous work has proposed a pyroclastic surge, possibly followed by aeolian reworking of the ash, for the depositional origin for these beds. We have performed a quantitative analysis of the structure, stratigraphy, and sedimentology at this location. Our results are consistent with an explosive volcaniclastic origin for the layered sediments. Analysis of bedding orientations over half of the circumference of Home Plate reveals a radially inward dipping structure, consistent with deposition in the volcanic vent, or topographic draping of a preexisting depression. Detailed observations of the sedimentology show that grain sorting varies significantly between outcrops on the east and west sides. Observations on the western side show a well-sorted population of sand sized grains which comprise the bedrock, while the eastern margin shows a wider range of grain sizes, including some coarse granules. These observations are consistent with primary deposition by a pyroclastic surge. However, aeolian reworking of the upper stratigraphic unit is not ruled out. Identification of explosive volcanic products on Mars may implicate magma interaction with subsurface hydrologic reservoirs in the past.

The diffusion coefficient of water vapor through porous media at Mars-like surface conditions is measured for a variety of complex particle size distributions and soil compositions. Micron-sized dust simulants, mixtures of sand- and dust-sized particles, and salt-encrusted sand are examined. We find that while the value of the diffusion coefficient, D, can be reduced by up to a factor of 10 for heavily salt-encrusted soils (minimum observed D = 0.4 ± 0.04 cm² s¹), moderate amounts of salt only produce minor reductions in D. Mechanical packing of pure dust can lower D by a similar amount, while mixtures of dust with sand-sized particles produce at most a factor of ~4 reduction. We conclude that present-day processes of aeolian redistribution, moderate levels of salt encrustation, and volatile loss from dirty ice would be inefficient at producing soil deposits and lags on Mars that pose significant barriers to diffusion. Therefore, subsurface ice deposits that are thermally unstable would not be protected against sublimative loss by such materials.

The Dynamic Albedo of Neutrons (DAN) instrument, a neutron scattering instrument currently being constructed by IKI, is a component of the science payload of the 2009 Mars Science Laboratory (MSL) mission. Based on simulations using the Monte-Carlo N-Particle Extended (MCNPX) particle physics code, DAN is able to measure bulk water content and to detect variations in water concentration up to ∼50 cm beneath the ground surface, assuming perfect detector performance. Data from DAN, combined with measurements from other instruments on MSL, allow derivation of profiles of water content to a depth of ∼15 cm.

THE RETURN OF MARTIAN SAMPLES TO EARTH has long been recognized as an essential component of a cycle of exploration that begins with orbital reconnaissance and in situ surface investigations. Major questions about life, climate, and geology require answers from state-of-the-art laboratories on Earth. Spacecraft instrumentation cannot perform critical measurements such as precise radiometric age dating, sophisticated stable isotopic analyses, and definitive life-detection assays. Returned sample studies could respond radically to unexpected findings, and returned materials could be archived for study by future investigators with even more capable laboratories. Unlike martian meteorites, returned samples could be acquired with known context from selected sites on Mars according to the prioritized exploration goals and objectives. The ND-MSR-SAG formulated the following 11 high-level scientific objectives that indicate how a balanced program of ongoing MSR missions could help to achieve the objectives and investigations described by MEPAG (2006). (1) Determine the chemical, mineralogical, and isotopic composition of the crustal reservoirs of carbon, nitrogen, sulfur, and other elements with which they have interacted and characterize carbon-, nitrogen-, and sulfur-bearing phases down to submicron spatial scales in order to document processes that could sustain habitable environments on Mars both today and in the past. (2) Assess the evidence for prebiotic processes, past life, and extant life on Mars by characterizing the signatures of these phenomena in the form of structure/morphology, biominerals, organic molecular and isotopic compositions, and other evidence within their geologic contexts. (3) Interpret the conditions of martian water-rock interactions through the study of their mineral products. (4) Constrain the absolute ages of major martian crustal geologic processes, including sedimentation, diagenesis, volcanism/plutonism, regolith formation, hydrothermal alteration, weathering, and cratering. (5) Understand paleoenvironments and the history of nearsurface water on Mars by characterizing the clastic and chemical components, depositional processes, and postdepositional histories of sedimentary sequences. (6) Constrain the mechanism and timing of planetary accretion, differentiation, and the subsequent evolution of the martian crust, mantle, and core. (7) Determine how the martian regolith was formed and modified and how and why it differs from place to place. (8) Characterize the risks to future human explorers in the areas of biohazards, material toxicity, and dust/granular materials and contribute to the assessment of potential in situ resources to aid in establishing a human presence on Mars. (9) For the present-day martian surface and accessible shallow subsurface environments, determine the preservation potential for the chemical signatures of extant life and prebiotic chemistry by evaluating the state of oxidation as a function of depth, permeability, and other factors. (10) Interpret the initial composition of the martian atmosphere, the rates and processes of atmospheric loss/gain over geologic time, and the rates and processes of atmospheric exchange with surface condensed species. (11) For martian climate-modulated polar deposits, determine their age, geochemistry, conditions of formation, and evolution through the detailed examination of the composition of water, CO₂, and dust constituents, as well as isotopic ratios and detailed stratigraphy of the upper layers of the surface. MSR would attain its greatest value if samples are collected as sample suites that represent the diversity of the products of various planetary processes. Sedimentary materials likely contain complex mixtures of chemical precipitates, volcaniclastics, impact glass, igneous rock fragments, and phyllosilicates. Aqueous sedimentary deposits are important for performing measurements of life detection, observations of critical mineralogy and geochemical patterns, and trapped gases. On Earth, hydrothermally altered rocks can preserve a record of hydrothermal systems that provided water, nutrients, and chemical energy necessary to sustain microorganisms. They also might have preserved fossils in their mineral deposits. Hydrothermal processes alter the mineralogy of crustal rocks and inject CO₂ and reduced gases into the atmosphere. Chemical alteration that occurs at nearsurface ambient conditions (typically 2O and CO₂ between near-surface solid materials and the atmosphere, and processes that involve fluids and sublimation. Regolith studies would help facilitate future human exploration by assessing toxicity and potential resources. Polar ices would constrain present and past climatic conditions and help elucidate water cycling. Surface ice samples from the Polar Layered Deposits (PLD) or seasonal frost deposits would help to quantify surface/atmosphere interactions. Short cores could help to resolve recent climate variability. Atmospheric gas samples would constrain the composition of the atmosphere and processes that influenced its origin and evolution. Trace organic gases (e.g., methane and ethane) could be analyzed for abundances, distribution, and relationships to a potential martian biosphere. Returned atmospheric samples that contain Ne, Kr, CO₂, CH₄, and C₂H₆ would confer major scientific benefits. Chemical and mineralogical analyses of martian dust would help to elucidate the weathering and alteration history of Mars. Given the global homogeneity of martian dust, a single sample is likely to be representative of the planet. A depth-resolved suite of samples should be obtained from depths that range from cm to several m within regolith or from rock outcrop to investigate trends in the abundance of oxidants (e.g., OH, HO₂, H₂O₂, and peroxy radicals), the effects of radiation, and the preservation of organic matter. Other sample suites include impact breccias that might sample rock types that are otherwise not available locally, tephra consisting of fine-grained regolith material or layers and beds possibly delivered from beyond the landing site, and meteorites whose alteration history could provide insights into martian climatic history. The following factors would affect our ability to achieve MSR's scientific objectives: (1) Sample size. A full program of scientific investigations would likely require samples of ≥8 g for bedrock, loose rocks, and finer-grained regolith. To support required biohazard testing, each sample requires an additional 2 g, leading to an optimal size of 10 g. Textural studies of some rock types might require one or more larger samples of ∼20 g. Material should remain to be archived for future investigations. (2) Number of samples. Studies of differences between samples could provide more information than detailed studies of a single sample. The number of samples needed to address MSR scientific objectives effectively is 35 (28 rock, 4 regolith, 1 dust, 2 gas). If the MSR mission recovers the MSL cache, it should also collect 26 additional samples (20 rock, 3 regolith, 1 dust and 2 atmospheric gas). The total mass of these samples is expected to be about 345 g (or 380 g with the MSL cache). The total returned mass with sample packaging would be about 700 g. (3) Sample encapsulation. To retain scientific value, returned samples must not commingle, each sample must be linked uniquely to its documented field context, and rocks should be protected against fragmentation during transport. A smaller number or mass of carefully managed samples is far more valuable than a larger number or mass of poorly managed samples. The encapsulation of at least some samples must retain any released volatile components. (4) Diversity of the returned collection. The diversity of returned samples must be commensurate with the diversity of rocks and regolith encountered. This guideline substantially influences landing site selection and rover operation protocols. It is scientifically acceptable for MSR to visit only a single site, but visiting 2 independent landing sites would be much more valuable. (5) In situ measurements for sample selection and documentation of field context. Relatively few samples could be returned from the vast array of materials the MSR rover would encounter; thus we must be able to choose wisely. At least 3 kinds of in situ observations are needed (color imaging, microscopic imaging, and mineralogy measurement) and possibly as many as 5 (also elemental analysis and reduced-carbon analysis). No significant difference exists in the observations needed for sample selection vs. sample documentation. Revisiting a previously occupied site might result in a reduction in the number of instruments.

Synthetic Aperture Radar (SAR) images of Titan's north polar region reveal quasi-circular to complex features which are interpreted to be liquid hydrocarbon lakes. We investigate methane transport in Titan's hydrologic cycle using the global distribution of lake features. As of May 2007, the SAR data set covers ∼22% of the surface and indicates multiple lake morphologies which are correlated across the polar region. Lakes are limited to latitudes above 55°N and vary from

We have derived the real and imaginary indices of refraction for 10 phyllosilicate minerals-montmorillonite, beidellite, nontronite, hectorite, saponite, illite, illite-smectite (60/40 interlayered) kaolinite, halloysite, and serpentine-from 100-2000 cm^-1 (5-100 μm) at 2 cm^-1 spectral sampling using classical Lorentz-Lorenz dispersion theory. We present the real and imaginary indices and the oscillator parameters with which they were modeled. Use of these optical constants will aid in the modeling of thermal infrared spectra of planets, asteroids, interplanetary and interstellar dust, and protoplanetary disks around nearby stars.

Slope streaks are surficial mass movements that are abundant in the dust-covered regions of Mars. Targeting of slope streaks seen in Viking images with the Mars Orbiter Camera provides observations of slope streak dust activity over two to three decades. In all study areas, new and persisting dark slope streaks are observed. Slope streaks disappeared in one area, with persisting streaks nearby. New slope streaks are found to be systematically darker than persisting streaks, which indicates gradual fading. Far more slope streaks formed at the study sites than have faded from visibility. The rate of formation at the study sites was 0.03 new slope streaks per existing streak per Mars year. Bright slope streaks do not presently form in sudden events as dark slope streaks do. Instead, bright streaks might form from old dark slope streaks, perhaps transitioning through a partially faded stage.

The diffusion coefficient of water vapor in unconsolidated porous media is measured for various soil simulants at Mars-like pressures and subzero temperatures. An experimental chamber which simultaneously reproduces a low-pressure, low-temperature, and low-humidity environment is used to monitor water flux from an ice source through a porous diffusion barrier. Experiments are performed on four types of simulants: 40-70 μm glass beads, sintered glass filter disks, 1-3 μm dust (both loose and packed), and JSC Mars-1. A theoretical framework is presented that applies to environments that are not necessarily isothermal or isobaric. For most of our samples, we find diffusion coefficients in the range of 2.8 to 5.4 cm² s^-1 at 600 Pascal and 260 K. This range becomes 1.9-4.7 cm² s^-1 when extrapolated to a Mars-like temperature of 200 K. Our preferred value for JSC Mars-1 at 600 Pa and 200 K is 3.7 ± 0.5 cm² s^-1. The tortuosities of the glass beads is about 1.8. Packed dust displays a lower mean diffusion coefficient of 0.38 ± 0.26 cm² s^-1, which can be attributed to transition to the Knudsen regime where molecular collisions with the pore walls dominate. Values for the diffusion coefficient and the variation of the diffusion coefficient with pressure are well matched by existing models. The survival of shallow subsurface ice on Mars and the providence of diffusion barriers are considered in light of these measurements.

High-latitude ground ice on Mars discovered by the Gamma Ray Spectrometer suite is thought to be thermally stable owing to the presence of vapor in the Martian atmosphere. However, local slopes can alter surface and subsurface temperatures substantially, and hence allow ground ice to persist at locations where it would otherwise be unstable. Global statistics of the topography of Mars are computed, processed, and extrapolated to derive a description of surface roughness on spatial scales to which ground ice should be sensitive. This slope distribution is convolved with a new thermal model for the dependence of subsurface ice on slope, to produce a prediction of the global ice distribution that includes the effect of topographic roughness. In the highest latitudes, slopes reduce the amount of buried ice, while in lower latitudes the ice fraction increases, widening the geographic boundary of the ice table. At the high latitudes, where ice is stable beneath horizontal ground, the estimated reduction of ice is small compared to the existing ice volume. Areas in the midlatitudes with high surface roughness that have previously been predicted to be ice free are predicted to contain quantities of ice that may be detectable at present and accessible in the future. Slopes cause ground ice to be stable to latitudes of about 25 degrees in both hemispheres, including, for example, areas within the northern Olympus Mons aureole deposits, Hecates Tholus, and Hellas basin. Ice is unstable at equatorial latitudes, even when accounting for surface slopes.

The eroded remains of a fluvial distributary network in Eberswalde crater are uniquely well preserved among similar structures on Mars. A quantitative analysis of the exposed stratigraphy has been performed to investigate the internal structure of the deposit. Using topographic information derived from stereo pairs of high-resolution Mars Orbiter Camera images, we have for the first time quantified the orientation of individual layers exposed along the distal end of the distributary network. In combination with topographic data from the Mars Orbiter Laser Altimeter, we have examined plausible scenarios for the formation of this structure. We find that the evidence is inconsistent with formation both as an alluvial fan and as a progradational delta. Instead, we find that an aggradational delta best fits the observed characteristics of the channel network and the Eberswalde basin as a whole. We conclude that the delta likely formed not in a stable long-lived lake but over the course of a small number of shorter lacustrine episodes, which were not sustained at equilibrium conditions.

Mars was most active during its first billion years. The core, mantle, and crust formed within ∼50 million years of solar system formation. A magnetic dynamo in a convecting fluid core magnetized the crust, and the global field shielded a more massive early atmosphere against solar wind stripping. The Tharsis province became a focus for volcanism, deformation, and outgassing of water and carbon dioxide in quantities possibly sufficient to induce episodes of climate warming. Surficial and near-surface water contributed to regionally extensive erosion, sediment transport, and chemical alteration. Deep hydrothermal circulation accelerated crustal cooling, preserved variations in crustal thickness, and modified patterns of crustal magnetization.

Observations by the Mars Orbiter Laser Altimeter have been used to detect subtle changes of the polar surface height during the course of seasonal cycles that correlate with the expected pattern of CO₂ deposition and sublimation. Using altimetric crossover residuals from the Mars Orbiter Laser Altimeter, we show that while zonally averaged data capture the global behavior of CO₂ exchange, there is a dependence of the pattern on longitude. At the highest latitudes the surface height change is as high as 1.5-2 m peak to peak, and it decreases equatorward. Decomposition of the signal into harmonics in time allows inspection of the spatial pattern and shows that the annual component is strongly correlated with the residual south polar cap deposits and, to a lesser extent, with the north polar cap. In the north, the second harmonic (semiannual) component correlates with the location of the ice deposits. The phases of the annual cycles are in agreement with observations by the Thermal Emission Spectrometer of the timing of the annual disappearance of CO₂ frost from the surface at the high latitudes. At lower latitudes, frost sublimation ("Crocus date") predates the mean depositional minima, as expected. These global-scale, volumetric measurements of the distribution of condensed CO₂ can be combined with measurements of the deposited column mass density derived from the Neutron Spectrometer on board Mars Odyssey to yield an estimate of the density of the seasonally exchanging material of 0.5 ± 0.1 g/cm³ . These constraints should be considered in models of the Martian climate system and volatile cycles.

The Mars Orbiter Camera (MOC) has imaged, sometimes repeatedly, mass movements known as slope streaks, which are abundant in the dust-covered regions on Mars. They are among the few known examples of contemporary surface changes. A survey of 173 collocated image pairs indicates that these features are currently forming at a high rate of ∼7% per existing streak, per Martian year. Either there is a complete turnover within a few decades or the streak population is currently increasing rapidly. Large spatial, as well as possible temporal, variations in the formation rate are obtained from these data. Streaks do not appear to fade over time periods comparable to their inverse formation rate of ∼28 years, as seen by analysis of Viking Orbiter images containing streaks that are still visible in MOC images. Gradual or stochastic variations in dust deposition may be needed to explain observations of changes in the formation rate, and its current imbalance with the fading rate.

From gravity and topography data collected by the Mars Global Surveyor spacecraft we calculate gravity/topography admittances and correlations in the spectral domain and compare them to those predicted from models of lithospheric flexure. On the basis of these comparisons we estimate the thickness of the Martian elastic lithosphere (T_e) required to support the observed topographic load since the time of loading. We convert T_e to estimates of heat flux and thermal gradient in the lithosphere through a consideration of the response of an elastic/plastic shell. In regions of high topography on Mars (e.g., the Tharsis rise and associated shield volcanoes), the mass-sheet (small-amplitude) approximation for the calculation of gravity from topography is inadequate. A correction that accounts for finite-amplitude topography tends to increase the amplitude of the predicted gravity signal at spacecraft altitudes. Proper implementation of this correction requires the use of radii from the center of mass (collectively known as the planetary "shape") in lieu of "topography" referenced to a gravitational equipotential. Anomalously dense surface layers or buried excess masses are not required to explain the observed admittances for the Tharsis Montes or Olympus Mons volcanoes when this correction is applied. Derived T_e values generally decrease with increasing age of the lithospheric load, in a manner consistent with a rapid decline of mantle heat flux during the Noachian and more modest rates of decline during subsequent epochs.

Measurements acquired by the Mars Orbiter Laser Altimeter on board the Mars Global Surveyor indicate that large drainage systems on Mars have geomorphic characteristics inconsistent with prolonged erosion by surface runoff. We find the topography has not evolved to an expected equilibrium terrain form, even in areas where runoff incision has been previously interpreted. By analogy with terrestrial examples, groundwater sapping may have played an important role in the incision. Longitudinally flat floor segments may provide a direct indication of lithologic layers in the bedrock, altering subsurface hydrology. However, it is unlikely that floor levels are entirely due to inherited structures due to their planar cross-cutting relations. These conclusions are based on previously unavailable observations, including extensive piece-wise linear longitudinal profiles, frequent knickpoints, hanging valleys, and small basin concavity exponents.

The Mars Orbiter Laser Altimeter (MOLA), an instrument on the Mars Global Surveyor spacecraft, has measured the topography, surface roughness, and 1.064-μm reflectivity of Mars and the heights of volatile and dust clouds. This paper discusses the function of the MOLA instrument and the acquisition, processing, and correction of observations to produce global data sets. The altimeter measurements have been converted to both gridded and spherical harmonic models for the topography and shape of Mars that have vertical and radial accuracies of ∼1 m with respect to the planet's center of mass. The current global topographic grid has a resolution of 1/64° in latitude x 1/32° in longitude (1 x 2 km² at the equator). Reconstruction of the locations of incident laser pulses on the Martian surface appears to be at the 100-m spatial accuracy level and results in 2 orders of magnitude improvement in the global geodetic grid of Mars. Global maps of optical pulse width indicative of 100-m-scale surface roughness and 1.064-μm reflectivity with an accuracy of 5% have also been obtained.

Measurements from the Near Earth Asteroid Rendezvous (NEAR)-Shoemaker Laser Rangefinder (NLR) indicate that asteroid 433 Eros is a consolidated body with a complex shape dominated by collisions. The offset between the asteroid's center of mass and center of figure indicates a small deviation from a homogeneous internal structure that is most simply explained by variations in mechanical structure. Regional-scale relief and slope distributions show evidence for control of topography by a competent substrate. Impact crater morphology is influenced by both gravity and structural control. Small-scale topography reveals ridges and grooves that may be generated by impact-related fracturing.

Elevations measured by the Mars Orbiter Laser Altimeter have yielded a high-accuracy global map of the topography of Mars. Dominant features include the low northern hemisphere, the Tharsis province, and the Hellas impact basin. The northern hemisphere depression is primarily a long-wavelength effect that has been shaped by an internal mechanism. The topography of Tharsis consists of two broad rises. Material excavated from Hellas contributes to the high elevation of the southern hemisphere and to the scarp along the hemispheric boundary. The present topography has three major drainage centers, with the northern lowlands being the largest. The two polar cap volumes yield an upper limit of the present surface water inventory of 3.2 to 4.7 million cubic kilometers.

Eighteen profiles of approximately N-S-trending topography from the Mars Orbiter Laser Altimeter (MOLA) are used to analyze the shape of Mars' northern hemisphere. MOLA observations show smaller northern hemisphere flattening than previously thought. The hypsometric distribution is narrowly peaked with >20% of the surface lying within 200 m of the mean elevation. Low elevation correlates with low surface roughness, but the elevation and roughness may reflect different mechanisms. Bouguer gravity indicates less variability in crustal thickness and/or lateral density structure than previously expected. The 3.1-km offset between centers of mass and figure along the polar axis results in a pole-to-equator slope at all longitudes. The N-S slope distribution also shows a subtle longitude-dependent variation that may represent the antipodal effect of the formation of Tharsis.