Publications

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 sites–two 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 Moon’s 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 10–20% 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.

Current lunar origin scenarios suggest that Earth's Moon may have resulted from the merger of two (or more) smaller moonlets. Dynamical studies of multiple moons find that these satellite systems are not stable, resulting in moonlet collision or loss of one or more of the moonlets. We perform Smoothed Particle Hydrodynamic (SPH) impact simulations of two orbiting moonlets inside the planetary gravitational potential and find that the classical outcome of two bodies impacting in free space is altered as erosive mass loss is more significant with decreasing distance to the planet. Depending on the conditions of accretion, each moonlet could have a distinct isotopic signature; therefore, we assess the initial mixing during their merger, in order to estimate whether future measurements of surface variations could distinguish between lunar origin scenarios (single vs. multiple moonlets). We find that for comparable-size impacting bodies in the accretionary regime, surface mixing is efficient, but in the hit-and-run regime, only small amount of material is transferred between the bodies. However, sequences of hit-and-run impacts are expected, which will enhance the surface mixing. Overall, our results show that large-scale heterogeneities can arise only from the merger of drastically different component masses. Surfaces of moons resulting from merger of comparable-sized components have little material heterogeneities, and such impacts are preferred, as the relatively massive impactor generates more melt, extending the lunar magma ocean phase.

The high planetary multiplicity revealed by Kepler implies that transit timing variations (TTVs) are intrinsically common. The usual procedure for detecting these TTVs is biased to long-period, deep transit planets, whereas most transiting planets have short periods and shallow transits. Here we introduce the Spectral Approach technique to TTVs that allows expanding the TTV catalog toward lower TTV amplitude, shorter orbital period, and shallower transit depth. In the spectral approach, we assume that a sinusoidal TTV exists in the data and then calculate the improvement to χ
² that this model allows over that of the linear-ephemeris model. This enables detection of TTVs even in cases where the transits are too shallow, so that individual transits cannot be timed. The spectral approach is more sensitive because it has fewer free parameters in its model. Using the spectral approach, we (a) detect 129 new periodic TTVs in Kepler data (an increase of ∼2/3 over a previous TTV catalog); (b) constrain the TTV periods of 34 long-period TTVs and reduce amplitude errors of known TTVs; and (c) identify cases of multi-periodic TTVs, for which absolute planetary mass determination may be possible. We further extend our analysis by using perturbation theory assuming a small TTV amplitude at the detection stage, which greatly speeds up our detection (to a level of few seconds per star). Our extended TTV sample shows no deficit of short-period or low-amplitude transits, in contrast to previous surveys, in which the detection schemes were significantly biased against such systems.

We used infrared data from the Lunar Reconnaissance Orbiter (LRO) Diviner Lunar Radiometer Experiment to globally map thermophysical properties of the Moon's regolith fines layer. Thermal conductivity varies from 7.4x10(-4)Wm(-1)K(-1) at the surface to 3.4x10(-3)Wm(-1)K(-1) at depths of similar to 1m, given density values of 1,100kgm(-3) at the surface to 1,800kgm(-3) at 1m depth. On average, the scale height of these profiles is similar to 7cm, corresponding to a thermal inertia of 552Jm(-2)K(-1)s(-1/2) at 273K, relevant to the diurnally active near-surface layer, similar to 4-7cm. The temperature dependence of thermal conductivity and heat capacity leads to an similar to 2 times diurnal variation in thermal inertia at the equator. On global scales, the regolith fines are remarkably uniform, implying rapid homogenization by impact gardening of this layer on timescales 100Jm(-2)K(-1)s(-1/2)) in the interiors and ejecta of Copernican-aged impact craters and lower thermal inertia (

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.

Seven of the nine known Mars Trojan asteroids belong to an orbital cluster 2 named after its largest member, (5261) Eureka. Eureka is probably the progenitor of the whole cluster, which formed at least 1Gyr ago(3). It has been suggested 3 that the thermal YORP (Yarkovsky-O'Keefe-Radzievskii-Paddack) effect spun up Eureka, resulting in fragments being ejected by the rotational-fission mechanism. Eureka's spectrum exhibits a broad and deep absorption band around 1 mu m, indicating an olivine-rich composition(4). Here we show evidence that the Trojan Eureka cluster progenitor could have originated as impact debris excavated from the Martian mantle. We present new near-infrared observations of two Trojans ((311999) 2007 NS2 and (385250) 2001 DH47) and find that both exhibit an olivine-rich reflectance spectrum similar to Eureka's. These measurements confirm that the progenitor of the cluster has an achondritic composition(4). Olivine-rich reflectance spectra are rare amongst asteroids' but are seen around the largest basins on Mars(6). They are also consistent with some Martian meteorites (for example, Chassigny(7)) and with the material comprising much of the Martian mantle(8,9). Using numerical simulations, we show that the Mars Trojans are more likely to be impact ejecta from Mars than captured olivine-rich asteroids transported from the main belt. This result directly links specific asteroids to debris from the forming planets.

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.

Fast and accurate radiative transfer methods are essential for modeling CO₂-rich atmospheres, relevant to the climate of early Earth and Mars, present-day Venus, and some exoplanets. Although such models already exist, their accuracy may be improved as better theoretical and experimental constraints become available. Here we develop a unidimensional radiative transfer code for CO₂-rich atmospheres, using the correlated k approach and with a focus on modeling early Mars. Our model differs from existing models in that it includes the effects of CO₂ collisional line mixing in the calculation of the line-by-line absorption coefficients. Inclusion of these effects results in model atmospheres that are more transparent to infrared radiation and, therefore, in colder surface temperatures at radiative-convective equilibrium, compared with results of previous studies. Inclusion of water vapor in the model atmosphere results in negligible warming due to the low atmospheric temperatures under a weaker early Sun, which translate into climatically unimportant concentrations of water vapor. Overall, the results imply that sustained warmth on early Mars would not have been possible with an atmosphere containing only CO₂ and water vapor, suggesting that other components of the early Martian climate system are missing from current models or that warm conditions were not long lived.

The dwarf planet Ceres may have an ice‐rich crust, and subsurface ice exposed by impacts or endogenic activity would be subject to sublimation. We model surface and subsurface temperatures on Ceres to assess lifetimes of water ice and other volatiles. Topographic shadowing allows a small but nonnegligible fraction (∼0.4%) of Ceres' surface to be perennially below the ∼110 K criterion for 1 Gyr of stability. These areas are found above 60° latitude. Other molecules (CH3OH, NH3, SO2, and CO2) may be cold trapped in smaller abundances. A model for the transport, gravitational escape, and photoionization of H2O molecules suggests net accumulation in the cold traps. Buried ice is stable within a meter for > 1 Gyr at latitudes higher than ∼50°. An illuminated polar cap of water ice would be stable within a few degrees of the poles only if it maintained a high albedo (>0.5) at present obliquity. If the obliquity exceeded 5° in the geologically recent past, then a putative polar cap would have been erased. At latitudes 0°–30°, ice is stable under solar illumination only briefly (∼10–100 years), unless it has high albedo and thermal inertia, in which case lifetimes of > 104 years are possible. Finally, a small hemispheric asymmetry exists due to the timing of Ceres' perihelion passage, which would lead to a detectable enhancement of ice in the northern hemisphere if the orbital elements vary slowly relative to the ice accumulation rate. Our model results are potentially testable during the Dawn science mission.

Unraveling the stratigraphic record is the key to understanding ancient climate and past climate changes on Mars (Grotzinger, J. et al. [2011]. Astrobiology 11, 77-87). Stratigraphic records of river deposits hold particular promise because rain or snowmelt must exceed infiltration plus evaporation to allow sediment transport by rivers. Therefore, river deposits when placed in stratigraphic order could constrain the number, magnitudes, and durations of the wettest (and presumably most habitable) climates in Mars history. We use crosscutting relationships to establish the stratigraphic context of river and alluvial-fan deposits in the Aeolis Dorsa sedimentary basin, 10°E of Gale crater. At Aeolis Dorsa, wind erosion has exhumed a stratigraphic section of sedimentary rocks consisting of at least four unconformity-bounded rock packages, recording three or more distinct episodes of surface runoff. Early deposits (>700. m thick) are embayed by river deposits (>400. m thick), which are in turn unconformably draped by fan-shaped deposits (900. m thick) unconformably drape all previous deposits.River deposits embay a dissected landscape formed of sedimentary rock. The river deposits are eroding out of at least two distinguishable units. There is evidence for pulses of erosion during the interval of river deposition. The total interval spanned by river deposits is >(1×10⁶-2×10⁷) yr, and this is extended if we include alluvial-fan deposits. Alluvial-fan deposits unconformably postdate thrust faults which crosscut the river deposits. This relationship suggests a relatively dry interval of >4×10⁷yr after the river deposits formed and before the fan-shaped deposits formed, based on probability arguments. Yardang-forming layered deposits unconformably postdate all of the earlier deposits. They contain rhythmite and their induration suggests a damp or wet (near-) surface environment. The time gap between the end of river deposition and the onset of yardang-forming layered deposits is constrained to >1×10⁸yr by the high density of impact craters embedded at the unconformity. The time gap between the end of alluvial-fan deposition and the onset of yardang-forming layered deposits was at least long enough for wind-induced saltation abrasion to erode 20-30m into the alluvial-fan deposits. We correlate the yardang-forming layered deposits to the upper layers of Gale crater's mound (Mt. Sharp/Aeolis Mons), and the fan-shaped deposits to Peace Vallis fan in Gale crater. Alternations between periods of low mean obliquity and periods of high mean obliquity may have modulated erosion-deposition cycling in Aeolis. This is consistent with the results from an ensemble of simulations of Solar System orbital evolution and the resulting history of the obliquity of Mars. 57 of our 61 simulations produce one or more intervals of continuously low mean Mars obliquity that are long enough to match our Aeolis Dorsa unconformity data.

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.

Cassini RADAR images of Titan's south polar region acquired during southern summer contain lake features which disappear between observations. These features show a tenfold increases in backscatter cross-section between images acquired one year apart, which is inconsistent with common scattering models without invoking temporal variability. The morphologic boundaries are transient, further supporting changes in lake level. These observations are consistent with the exposure of diffusely scattering lakebeds that were previously hidden by an attenuating liquid medium. We use a two-layer model to explain backscatter variations and estimate a drop in liquid depth of approximately 1-m-per-year. On larger scales, we observe shoreline recession between ISS and RADAR images of Ontario Lacus, the largest lake in Titan's south polar region. The recession, occurring between June 2005 and July 2009, is inversely proportional to slopes estimated from altimetric profiles and the exponential decay of near-shore backscatter, consistent with a uniform reduction of 4 ± 1.3. m in lake depth.Of the potential explanations for observed surface changes, we favor evaporation and infiltration. The disappearance of dark features and the recession of Ontario's shoreline represents volatile transport in an active methane-based hydrologic cycle. Observed loss rates are compared and shown to be consistent with available global circulation models. To date, no unambiguous changes in lake level have been observed between repeat images in the north polar region, although further investigation is warranted. These observations constrain volatile flux rates in Titan's hydrologic system and demonstrate that the surface plays an active role in its evolution. Constraining these seasonal changes represents the first step toward our understanding of longer climate cycles that may determine liquid distribution on Titan over orbital time periods.

The acquisition of new global elevation data from the Lunar Orbiter Laser Altimeter, carried on the Lunar Reconnaissance Orbiter, permits quantification of the surface roughness properties of the Moon at unprecedented scales and resolution. We map lunar surface roughness using a range of parameters: median absolute slope, both directional (along-track) and bidirectional (in two dimensions); median differential slope; and Hurst exponent, over baselines ranging from ∼17 m to ∼2.7 km. We find that the lunar highlands and the mare plains show vastly different roughness properties, with subtler variations within mare and highlands. Most of the surface exhibits fractal-like behavior, with a single or two different Hurst exponents over the given baseline range; when a transition exists, it typically occurs near the 1 km baseline, indicating a significant characteristic spatial scale for competing surface processes. The Hurst exponent is high within the lunar highlands, with a median value of 0.95, and lower in the maria (with a median value of 0.76). The median differential slope is a powerful tool for discriminating between roughness units and is useful in characterizing, among other things, the ejecta surrounding large basins, particularly Orientale, as well as the ray systems surrounding young, Copernican-age craters. In addition, it allows a quantitative exploration on mare surfaces of the evolution of surface roughness with age.

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^–1), 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.

We explore the capability of a method of mapping the depth distribution of a hydrogen-rich layer in the top meter of Mars from the neutron currents measured by the Mars Odyssey Neutron Spectrometer. Assuming the soil can be modeled by two layers of known composition having different hydrogen contents, simulations allow an inversion of the neutron data into knowledge of depth and hydrogen content of the lower layer. The determination of these variables is sensitive to the hypothesis of chemical composition of the soil. We quantify this contribution to the uncertainty in the method first in terms of individual chemical elements and then in terms of macroscopic absorption cross sections. To minimize this source of error, an average composition was inferred from Mars Exploration Rover data. Possible compositions having a wide range of macroscopic absorption cross sections were used to evaluate the uncertainty associated with our calculations. We finally compare our results to ice table depth estimates predicted by two published theoretical models at locations where the composition is relatively well known. The fit is excellent in the southern high latitudes but questionable in the northern high latitudes. Possible explanations of these differences include the high geographical variations of the neutron currents relative to the spatial width of the response function of the instrument and the overly simple model we, of necessity, used for surface layering.

Constraints on the mineralogical composition of low-albedo, low-sulfur sands at Meridiani Planum are determined from Mars Exploration Rover (MER) Opportunity Miniature Thermal Emission Spectrometer (Mini-TES), Mössbauer, and Alpha Proton X-Ray Spectrometer measurements. Results of this work show that the sand is olivine basaltic in composition, with minor amounts of sulfate and a high-silica phase (glass or secondary amorphous silica). Measurements from all three instruments indicate that pyroxene is twice as abundant as olivine, and that the pyroxene composition is dominated by the low-calcium variety. The volume abundance of olivine is constrained to be 10-15%. Results from detailed analyses of MER data are used to ground truth the spectral emissivity and mineral abundances derived from orbit with Mars Global Surveyor TES data. TES-derived mineral abundances are within 5% of those derived from MER data, which is generally within the statistical errors associated with TES-derived phase abundances. The agreement lends support to global- and regional-scale variations in mineralogical composition determined from TES data in previous studies. An alternative method of least squares minimization is used for modeling the TES and Mini-TES data; the benefits of this method are demonstrated by comparison with conventional least squares techniques previously used by TES data users.

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.

The Mars Orbiter Camera on board the Mars Global Surveyor spacecraft has returned images of numerous dark streaks that are the result of down-slope mass movement occurring under present-day martian climatic conditions. We systematically analyze over 23,000 high-resolution images and demonstrate that slope streaks form exclusively in regions of low thermal inertia (confirming earlier results), steep slopes, and, remarkably, only where peak temperatures exceed 275 K. The northernmost streaks, which form in the coldest environment, form preferentially on warmer south-facing slopes. Repeat images of sites with slope streaks show changes only if the time interval between the two images includes the warm season. Surprisingly (in light of the theoretically short residence time of H₂O close to the surface), the data support the possibility that small amounts of water are transiently present in low-latitude near-surface regions of Mars and undergo phase transitions at times of high insolation, triggering the observed mass movements.

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.