Search

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 km2, 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. A reconstruction of the distribution of cold traps on the Moon at spatial scales varying from 1 km to 1 cm shows that the smallest ones are also the most numerous, 10–20% of the total. The total surface area of the Moon that can efficiently trap water is revised substantially upward, to 40,000 km2.

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-coplanar 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 without previously reported mass constraints. Overall we provide: (1) Estimates on physical and orbital properties for all systems analyzed. (2) 102 planets with mass detections significant to better than 3 standard deviations, 43 of which are lighter than five Earth masses. (3) 35 TbVs significant to better than 3 standard deviations. We focus on select systems showing strong TbVs, which can result from either interaction among the known transiting planets, or with a non-transiting object, and provide: (4) 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 multi-decade photometric and composite (RV, astrometry) data sets., 49 pages, 10 figures, accepted for publication. Machine-readable tables are included in the source files

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 meters may reveal climate oscillations unseen in the layer thicknesses and would thus probe recent precession-driven climate cycles., Comment: Paper accepted for publication in Icarus

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 a 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 the 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., accepted for publication. 17 pages, 9 Figures

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.

Introduction: A common outcome of a giant impact event is the formation of a circumplanetary debris disk, and in some cases, the capture of the surviving impactor, which results in a system composed of a primary, secondary, and a debris disk. The material in the debris disk may be accreted by the primary body (and secondary, if it exists), escape the system, or coalesce into larger clumps. A classic example of a system thought to originate by such an impact event is the Pluto-Charon binary system (Canup, 2005, 2011), in which all six bodies (including the other smaller satellites) lie approximately on the same plane and have nearly circular orbits. Such impacts are simulated using smoothed particle hydrodynamics (SPH) codes including self gravity. The equation of state (EOS) governs the relationship among the thermodynamic variables of the simulated material. The specific EOS used in the simulation may influence the final post-impact structure. In order to quantify the typical differences, we compare two approaches, one simple analytic and one tabulated EOS. Tillotson (Tillotson, 1962) is a widely used analytic EOS and is computationally fast, but it lacks important details such as the treatment of phase changes. Sesame (Bennett et al., 1978) is a commonly used tabulated EOS and is more accurate, but is computationally slower and may be poorly sampled in the required thermodynamic phase-space. Here we show a set of SPH impact simulations that assume similar geometric and dynamic initial conditions but different EOS. Methods: We performed ~100 SPH simulations using SWIFT code (Schaller et al., 2018), simulating Pluto-like impacts with 105-106 particles. The initial bodies are assumed to be differentiated, with target to impactor mass ratio of 1 or 7/3, impact angle, ξ, of 0, 30, 45, 60°. Impact velocity was chosen to be relatively small, 1-1.1 times the escape velocity. The impactor was either spin-less or rotating with a period of 5 or 10 hours. This parameter space was motivated by previous simulations (Canup, 2011) for the formation of the Pluto-Charon system. We compare results using the Tillotson and Sesame EOS. The simulations were stopped after 4 days, the typical time for the central body to relax to a stable spherical shape. We developed an algorithm to detect post-impact clumps. Two particles were considered in contact if their mutual distance was smaller than their combined smoothing lengths. The orbital elements of each clump (“satellitesimal”, defined as 100 particles in pairwise contact, equivalent to ~10-3Mpluto) were computed and studied. Results: Disk systems were formed for impact angles >30°, in using both EOSs. Satellitesimals, when formed, showed different properties. Figure 1 panels (a,b) show the final snapshots of an impact that produced a debris disk and several satellitesimals highlighted in color corresponding to their masses. In this example, a more massive debris disk with a larger number of satellitesimals is obtained when using Tillotson EOS than using Sesame. Moreover, in addition to the target body, at least one large satellitesimal was formed in each of the simulations, but their composition, mass, and orbital elements differ between the EOSs. In comparing the Sesame run to Tillotson, the largest satellitesimal has a mass of 0.015Mpluto (with water fraction of 0.16) and 0.008Mpluto (with water fraction of 0.30) respectively. Its orbital elements are e=0.67, a=10.07Rpluto for Sesame runs, and e=0.15, a=3.2Rpluto for Tillotson. At smaller masses, using Tillotson EOS produces a greater number of clumps, as seen in Figure 1c. Three more satellitesimals were formed (some beyond the plot limits of 1b), with e=0.86, 1.35, 0.10. Note that the satellitesimal mass is an order of magnitude smaller than Charon, so alone they do not predict Charon’s formation (Canup, 2011). In the final snapshot, a greater fraction of debris disk particles lie within the Roche limit (computed using present-day Charon’s density (McKinnon et al., 2017)) in the Sesame simulation, whereas the Tillotson disk extends to a greater distance. In terms of mass, angular momentum, and composition, the debris disks are similar (0.025 versus 0.033 of the total mass; 0.19 versus 0.25 of the total angular momentum; water fraction of 0.44 versus 0.37 for Sesame and Tillotson respectively). We note the execution time was 2-3 times longer for the simulations using Sesame than Tillotson, a factor which may be considered in choosing EOS. Exploring the parameter space, we note that head-on impacts (ξ~0°) produce a merged single body, with the vast majority of ejected particles accreted by the planet, and the rest ejected to space. Oblique, faster than escape velocity impacts (vimp/vesc=1.1 and ξ=60°) resulted in two unbound bodies, with little mass transfer between the two, consistent with previous studies of planetary impacts (Leinhardt & Stewart, 2012). References: Bennett, B. I., Johnson, J. D., Kerley, G. I., & Rood, G. T. (1978). Recent developments in the Sesame equation-of-state library. https://doi.org/10.2172/5150206 Canup, R. M. (2005). A Giant Impact Origin of Pluto-Charon. Science, 307(5709), 546–550. Canup, R. M. (2011). On a Giant Impact Origin of Charon, Nix, and Hydra. The Astronomical Journal, 141(2), 35. Leinhardt, Z. M., & Stewart, S. T. (2012). Collisions Between Gravity-Dominated Bodies. I. Outcome Regimes and Scaling Laws. The Astrophysical Journal, 745(1), 79. McKinnon, W. B., Stern, S. A., Weaver, H. A., Nimmo, F., Bierson, C. J., Grundy, W. M., et al. (2017). Origin of the Pluto–Charon system: Constraints from the New Horizons flyby. Icarus. https://doi.org/10.1016/j.icarus.2016.11.019 Schaller, M., Gonnet, P., Chalk, A. B. G., & Draper, P. W. (2018, May 1). SWIFT: SPH With Inter-dependent Fine-grained Tasking. Astrophysics Source Code Library. Retrieved from https://ui.adsabs.harvard.edu/abs/2018ascl.soft05020S Tillotson, J. H. (1962). Metallic equations of state for hypervelocity impact (No. Rep. GA-3216 ). General Dynamics San Diego CA.

Impacts between planetary-sized bodies can explain the origin of satellites orbiting large ( R > 500 km) trans-Neptunian objects. Their water rich composition, along with the complex phase diagram of water, make it important to accurately model the wide range of thermodynamic conditions material experiences during an impact event and in the debris disk. Since differences in the thermodynamics may influence the system dynamics, we seek to evaluate how the choice of an equation of state (EOS) alters the system’s evolution. Specifically, we compare two EOSs that are constructed by different approaches: either by a simplified analytic description (Tillotson), or by interpolation of tabulated data (Sesame). Approximately 50 pairs of Smoothed Particle Hydrodynamics impact simulations were performed, with similar initial conditions but different EOSs, in the parameter space in which the Pluto–Charon binary is thought to form (slow impacts between Pluto-size, water rich bodies). Generally, we show that impact outcomes (e.g., circumplanetary debris disk) are consistent between EOSs. Some differences arise, importantly in the production of satellitesimals (large intact clumps) that form in the post-impact debris disk. When utilizing an analytic EOS, the emergence of satellitesimals is highly certain, while when using the tabulated EOS it is less common. This is because for the typical densities and energies experienced in these impacts, the analytic EOS predicts very low pressure values, leading to particles artificially aggregating by a tensile instability.

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.4×10-4 W m-1 K-1 at the surface, to 3.4×10-3 W m-1 K-1 at depths of ~1 m, given density values of 1100 kg m-3 at the surface, to 1800 kg m-3 at 1-m depth. On average, the scale height of these profiles is ~7 cm, corresponding to a thermal inertia of 55 ±2 J m-2 K-1 s-1/2 at 273 K, relevant to the diurnally active near-surface layer, ~4-7 cm. The temperature-dependence of thermal conductivity and heat capacity leads to a ~2x 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 100 J m-2 K-1 s-1/2) in the interiors and ejecta of Copernican-aged impact craters, and lower thermal inertia (< 50 J m-2 K-1 s-1/2) within the lunar cold spots identified by Bandfield et al. (2014). Observed trends in ejecta thermal inertia provide a potential tool for age-dating craters of previously unknown age, complementary to the approach suggested by Ghent et al. (2014). Several anomalous regions are identified in the global 128 pixels-per-degree maps presented here, including a high-thermal inertia deposit near the antipode of Tycho crater.

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.

Non-Keplerian dynamics of planetary orbits manifest in the transit light-curve as variations of different types. In addition to Transit Timing Variations (TTV's), 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 non-transiting companions, or the non-spherical 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 followup observations, which are suggested in order to place further constraints on this unique Kepler object., Comment: 16 pages, 5 figures

Inferring planetary parameters from transit timing variations (TTVs) is challenging for small exoplanets because their transits may be so weak that determination of individual transit timing is difficult or impossible. We implement a useful combination of tools that together provide a numerically fast global photodynamical model. This is used to fit the TTV-bearing light curve, in order to constrain the masses of transiting exoplanets in low-eccentricity, multiplanet systems—and small planets in particular. We present inferred dynamical masses and orbital eccentricities in four multi-planet systems from Kepler's complete long-cadence data set. We test our model against Kepler-36/KOI-277, a system with some of the most precisely determined planetary masses through TTV inversion methods, and find masses of 5.56+0.41 −0.45 and 9.76+0.79 −0.89 for Kepler-36 b and c, respectively—consistent with literature in both value and error. We then improve the mass determination of the four planets in Kepler-79/KOI-152, where literature values were physically problematic to 12.5+4.5 −3.6, 9.5+2.3 −2.1, 11.3+2.2 −2.2 and 6.3+1.0 −1.0 for Kepler-79 b, c, d, and e, respectively. We provide new mass constraints where none existed before for two systems. These are 12.5+3.2 −2.6 for Kepler-450 c, and 3.3+1.7 −1.0 and 17.4+7.1 −3.8 for Kepler-595 c (previously KOI-547.03) and b, respectively. The photodynamical code used here, called PyDynamicaLC, is made publicly available.

The hypothesis of lunar origin by a single giant impact can explain some aspects of the Earth-Moon system. However, it is difficult to reconcile giant impact models with the compositional similarity of the Earth and Moon without violating angular momentum constraints. Furthermore, successful giant impact scenarios require very specific conditions such that they have a low probability of occurring. Here we present numerical simulations suggesting that the Moon could instead be the product of a succession of a variety of smaller collisions. In this scenario, each collision forms a debris disk around the proto-Earth that then accretes to form a moonlet. The moonlets tidally advance outward, and may coalesce to form the Moon. We find that sub-lunar moonlets are a common result of impacts expected onto the proto-Earth in the early solar system and find that the planetary rotation is limited by impact angular momentum drain. We conclude that, assuming efficient merger of moonlets, a multiple impact scenario can account for the formation of the Earth-Moon system with its present properties.

Fast and accurate radiative transfer methods are essential for modeling CO2-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 CO2-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 CO2 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 CO2 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 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 Myrs 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 pre-existing moonlets survive subsequent impact events, and explore the dynamics of Earth-moonlet systems that contain two moonlets generated Myrs 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., 12 pages, 8 figures, accepted by ApJ

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 non-negligible 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, 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 yr), unless it has high albedo and thermal inertia, in which case lifetimes of > 104 yr 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.

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 (β ∼ 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 (∼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.

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^6–2 × 10^7) 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^7 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^8 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–30 m 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 high planetary multiplicity revealed by Kepler implies that Transit Time 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 to TTVs technique that allows expanding the TTVs catalog towards 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 $\chi^2$ this model allows over that of linear ephemeris model. This enables detection of TTVs even in cases where the transits are too shallow so individual transits cannot be timed. The Spectral Approach is more sensitive due to the reduced number of free parameters in its model. Using the Spectral Approach, we: (a) detect 131 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; (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 small TTV amplitude at the detection stage, which greatly speeds up our detection (to a level of few seconds per star). Our extended TTVs 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., Comment: ApJS accepted. 23 pages, 40 figures

Seven of the nine known Mars Trojan asteroids belong to an orbital cluster named after its largest member 5261 Eureka. Eureka is likely the progenitor of the whole cluster, which formed at least 1 Gyr ago. It was suggested that the thermal YORP effect spun-up Eureka resulting with 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. 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. Olivine-rich reflectance spectra are rare amongst asteroids but are seen around the largest basins on Mars. They are also consistent with some Martian meteorites (e.g. Chassigny), and with the material comprising much of the Martian mantle. 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 links directly specific asteroids to debris from the forming planets., Comment: 10 text pages, 3 figures, 3 SI pages

We find that the reflectance of the lunar surface within 5° of latitude of the South Pole increases rapidly with decreasing temperature, near ∼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° from the pole. This South pole reflectance anomaly persists when analysis is limited to surfaces with slopes less than 10° 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.

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.

Asteroid pairs 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 greater than 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 gr cm-3 measured for the rubble pile structured 25143 Itokawa, 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., Comment: Manuscript under review in Icarus

International audience; Dune fields on Titan cover more than 17% of the moon's surface, constituting the largest known surface reservoir of organics. Their confinement to the equatorial belt, shape, and eastward direction of propagation offer crucial information regarding both the wind regime and sediment supply. Herein, we present a comprehensive analysis of Titan's dune orientations using automated detection techniques on nonlocal denoised radar images. By coupling a new dune growth mechanism with wind fields generated by climate modeling, we find that Titan's dunes grow by sediment transport on a nonmobile substratum. To be fully consistent with both the local crestline orientations and the eastward propagation of Titan's dunes, the sediment should be predominantly transported by strong eastward winds, most likely generated by equinoctial storms or occasional fast westerly gusts. Additionally, convergence of the meridional transport predicted in models can explain why Titan's dunes are confined within ±30 • latitudes, where sediment fluxes converge.

Sedimentary rocks preserved on the surface of Mars represent a natural archive of past climate conditions. Although the details of their formation often remain poorly constrained, the recent detection of rhythmic bedding patterns in the Arabia Terra region suggests the influence of orbital variations on sedimentary deposition. Here we detail a number of new sites which exhibit quasiperiodic stratigraphic variations, demonstrating their occurrence throughout the equatorial region of the planet. We characterize these recorded signals as well as the local geomorphic context and structural attributes. Two cyclic units are identified within Gale crater, the landing site of the Mars Science Laboratory mission, enabling estimation of possible formation timescales for the geologic units that may be studied in situ by the rover. We find a general lack of fluvial features in connection with rhythmic geologic units, contrasting these sites with the aperiodic deltaic stratigraphy found at Eberswalde crater. Possible formation scenarios and their climatic implications are discussed for the diverse set of quasiperiodic sedimentary units. We propose multiple depositional pathways for recording cyclic climate changes, including repeated evaporitic precipitation from groundwater discharge in topographic lows as well as largely anhydrous accumulation of atmospheric dust for deposits outside of confined basins. The preservation of orbital signals in sediments distributed across a wide range of geographic settings suggests a pervasive influence on Martian climate conditions through time.

We present maps of the topographic roughness of the Moon at hectometer and kilometer scales. The maps are derived from range profiles obtained by the Lunar Orbiter Laser Altimeter (LOLA) instrument onboard the Lunar Reconnaissance Orbiter (LRO) spacecraft. As roughness measures, we used the interquartile range of profile curvature at several baselines, from 115 m to 1.8 km, and plotted these in a global map format. The maps provide a synoptic overview of variations of typical topographic textures and utilize the exceptional ranging precision of the LOLA instrument. We found that hectometer-scale roughness poorly correlates with kilometer-scale roughness, because they reflect different sets of processes and time scales. Hectometer-scale roughness is controlled by regolith accumulation and modification processes and affected by the most recent events, primarily, geologically recent (1–2 Ga) meteoritic impacts. Kilometer-scale roughness reflects major geological (impact, volcanic and tectonic) events in earlier geological history. Young large impact craters are rough, and their roughness decreases with age. The global roughness maps revealed a few unusually dense clusters of hectometer- and decameter-size impact craters that differ in their morphology and settings from typical secondary crater clusters and chains; the origin of these features is enigmatic. The maps can assist in the geological mapping of the lunar maria by revealing contacts between volcanic plain units. The global roughness maps also clearly reveal cryptomaria, old volcanic plains superposed by younger materials, primarily crater and basin ejecta.

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.

Cassini RADAR SARtopo and altimetry data are used to construct a global gridded 1 × 1° elevation map, for use in Global Circulation Models, hydrological models and correlative studies. The data are sparse, and so most of the map domain (∼90%) is populated with interpolated values using a spline algorithm. The highest (∼+520 m) gridded point observed is at 48°S, 12°W. The lowest point observed (∼1700 m below a 2575 km sphere) is at 59°S, 317°W: this may be a basin where liquids presently in the north could have resided in the past. If the deepest point were once a sea with the areal extent of present-day Ligeia Mare, it would be ∼1000 m deep. We find four prominent topographic rises, each ∼200 km wide, radar-bright and heavily dissected, distributed over a ∼3000 km arc in the southeastern quadrant of Titan (∼40–60°S, 15–150°W).

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 (PTF) survey to derive upper limits on the population size of active main-belt comets (MBCs). From data collected March 2009 through July 2012, we extracted 2 million observations of 220 thousand known main-belt objects (40% of the known population, down to 1-km diameter) and discovered 626 new objects in multi-night 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-main-belt comets (two of which were originally designated as asteroids by other surveys). We demonstrate a 66% detection efficiency with respect to the extendedness distribution of the 115 sampled comets, and a 100% detection efficiency with respect to extendedness levels greater than or equal to those we observed in the known candidate MBCs P/2010 R2 (La Sagra) and P/2006 VW139. Using a log-constant prior, we infer 95% confidence upper limits of 33 and 22 active MBCs (per million main-belt asteroids down to 1-km diameter), for detection efficiencies of 66% and 100%, respectively. In a follow-up to this morphological search, we will perform a photometric (disk-integrated brightening) search for MBCs., 19 pages, 15 figures. Accepted 30-May-2013 to Monthly Notices of the Royal Astronomical Society (MNRAS)

At a mean diameter of ~650 m, the near-Earth asteroid (455213) 2001 OE84 (OE84 for short) has a rapid rotation period of 0.486542+-0.000002 hours, which is uncommon for asteroids larger than ~200 m. We revisited OE84 14 years after it was first, and last, observed by Pravec et al. (2002) in order to measure again its spin rate and to search for changes. We have confirmed the rapid rotation and, by fitting the photometric data from 2001 and 2016 using the lightcurve inversion technique, we determined a retrograde sense of rotation, with the spin axis close to the ecliptic south pole; an oblate shape model of a/b=1.32+-0.04 and b/c=1.8+-0.2; and no change in spin rate between 2001 and 2016. Using these parameters we constrained the body's internal strength, and found that current estimations of asteroid cohesion (up to ~80 Pa) are insufficient to maintain an intact rubble pile at the measured spin rate of OE84. Therefore, we argue that a monolithic asteroid, that can rotate at the rate of OE84 without shedding mass and without slowing down its spin rate, is the most plausible for OE84, and we give constraints on its age, since the time it was liberated from its parent body, between 2-10 million years., Comment: 11 pages, 5 figures, 2 tables. Accepted for publication in Icarus

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 lightcurves 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 analyzed as part of this program. These fields, which cover an area of 21 deg^2, were observed on four nights with a cadence of ~20 min. Our pipeline was able to detect 624 asteroids, of which 145 (~20%) 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 lightcurves 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 10,000 asteroids mainly in the magnitude range ~14 to ~20.

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 15 K 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 2 K 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.

All planetary bodies with old surfaces exhibit planetary-scale impact craters: vast scars caused by the large impacts at the end of Solar System accretion or the late heavy bombardment. Here we investigate the geophysical consequences of planetary-scale impacts into a Mars-like planet, by simulating the events using a smoothed particle hydrodynamics (SPH) model. Our simulations probe impact energies over two orders of magnitude (2 × 10^(27)–6 × 10^(29) J), impact velocities from the planet’s escape velocity to twice Mars’ orbital velocity (6–50 km/s), and impact angles from head-on to highly oblique (0–75°). The simulation results confirm that for planetary-scale impacts, surface curvature, radial gravity, the large relative size of the impactor to the planet, and the greater penetration of the impactor, contribute to significant differences in the geophysical expression compared to small craters, which can effectively be treated as acting in a half-space. The results show that the excavated crustal cavity size and the total melt production scale similarly for both small and planetary-scale impacts as a function of impact energy. However, in planetary-scale impacts a significant fraction of the melt is sequestered at depth and thus does not contribute to resetting the planetary surface; complete surface resetting is likely only in the most energetic (6 × 10^(29) J), slow, and head-on impacts simulated. A crater rim is not present for planetary-scale impacts with energies >10^(29) J and angles ≤45°, but rather the ejecta is more uniformly distributed over the planetary surface. Antipodal crustal removal and melting is present for energetic (>10^(29) J), fast (>6 km/s), and low angle (≤45°) impacts. The most massive impactors (with both high impact energy and low velocity) contribute sufficient angular momentum to increase the rotation period of the Mars-sized target to about a day. Impact velocities of >20 km/s result in net mass erosion from the target, for all simulated energies and angles. The hypothesized impact origin of planetary structures may be tested by the presence and distribution of the geochemically-distinct impactor material.

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.

A set of lakes filled or partially filled with liquid hydrocarbon and empty lake basins have been discovered in the high latitudes of Saturn's moon Titan. These features were mapped by the radar instrument on the Cassini orbiter. Here we quantify the distribution of the lakes and basins, and show a pronounced hemispheric asymmetry in their occurrence. Whereas significant fractions of the northern high latitudes are covered by filled and empty lakes, the same latitudes in the southern hemisphere are largely devoid of such features. We propose that in addition to known seasonal changes, the observed difference in lake distribution may be caused by an asymmetry in the seasons on Titan that results from the eccentricity of Saturn's orbit around the Sun. We suggest that the consequent hemispheric difference in the balance between evaporation and precipitation could lead to an accumulation of lakes in one of Titan's hemispheres. This effect would be modulated by, and reverse with, dynamical variations in the orbit. We propose that much like in the Earth's glacial cycles, the resulting vigorous hydrologic cycle has a period of tens of thousands of years and leads to active geologic surface modification in the polar latitudes.

Global acquisition of infrared spectra and high-resolution visible and infrared imagery has enabled the placement of compositional information within stratigraphic and geologic context. Mare Serpentis, a low albedo region located northwest of Hellas Basin, is rich in spectral and thermophysical diversity and host to numerous isolated exposures of in situ rocky material. Most martian surfaces are dominated by fine-grained particulate materials that bear an uncertain compositional and spatial relationship to their source. Thus location and characterization of in situ rock exposures is important for understanding the origin of highland materials and the processes which have modified those materials. Using spectral, thermophysical and morphologic information, we assess the local and regional stratigraphy of the Mare Serpentis surface in an effort to reconstruct the geologic history of the region. The martian highlands in Mare Serpentis are dominated by two interspersed surface units, which have distinct compositional and thermophysical properties: (1) rock-dominated surfaces relatively enriched in olivine and pyroxene, and depleted in high-silica phases, and (2) sediment or indurated material depleted in olivine and pyroxene, with relatively higher abundance of high-silica phases. This is a major, previously unrecognized trend which appears to be pervasive in the Mare Serpentis region and possibly in other highland areas. The detailed observations have led us to form two hypotheses for the relationship between these two units: either (1) they are related through a widespread mechanical and/or chemical alteration process, where less-mafic plains materials are derived from the mafic bedrock, but have been compositionally altered in the process of regolith formation, or (2) they are stratigraphically distinct units representing separate episodes of upper crust formation. Existing observations suggest that the second scenario is more likely. In this scenario, plains materials represent older, degraded, and possibly altered, “basement” rock, whereas the rocky exposures represent later additions to the crust and are probably volcanic in origin. These hypotheses should be further testable with decimeter-resolution imagery and meter-resolution short wavelength infrared spectra.

Current understanding of weather, climate and global atmospheric circulation on Mars is incomplete, in particular at altitudes above about 30 km. General circulation models for Mars are similar to those developed for weather and climate forecasting on Earth and require more martian observations to allow testing and model improvements. However, the available measurements of martian atmospheric temperatures, winds, water vapour and airborne dust are generally restricted to the region close to the surface and lack the vertical resolution and global coverage that is necessary to shed light on the dynamics of Mars middle atmosphere at altitudes between 30 and 80 km (ref.7). Here we report high-resolution observations from the Mars Climate Sounder instrument on the Mars Reconnaissance Orbiter. These observations show an intense warming of the middle atmosphere over the south polar region in winter that is at least 10-20 K warmer than predicted by current model simulations. To explain this finding, we suggest that the atmospheric downwelling circulation over the pole, which is part of the equator-to-pole Hadley circulation, may be as much as 50 more vigorous than expected, with consequences for the cycles of water, dust and CO"2 that regulate the present-day climate on Mars. © 2008 Macmillan Publishers Limited.

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.

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.

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.