The next step on the path toward another Earth is to find atmospheres similar to those of Earth and Venus - high-molecular-weight (secondary) atmospheres - on rocky exoplanets. Many rocky exoplanets are born with thick (> 10 kbar) H$_2$-dominated atmospheres but subsequently lose their H$_2$; this process has no known Solar System analog. We study the consequences of early loss of a thick H$_2$ atmosphere for subsequent occurrence of a high-molecular-weight atmosphere using a simple model of atmosphere evolution (including atmosphere loss to space, magma ocean crystallization, and volcanic outgassing). We also calculate atmosphere survival for rocky worlds that start with no H$_2$. Our results imply that most rocky exoplanets orbiting closer to their star than the Habitable Zone that were formed with thick H$_2$-dominated atmospheres lack high-molecular-weight atmospheres today. During early magma ocean crystallization, high-molecular-weight species usually do not form long-lived high-molecular-weight atmospheres; instead they are lost to space alongside H$_2$. This early volatile depletion also makes it more difficult for later volcanic outgassing to revive the atmosphere. However, atmospheres should persist on worlds that start with abundant volatiles (for example, waterworlds). Our results imply that in order to find high-molecular-weight atmospheres on warm exoplanets orbiting M-stars, we should target worlds that formed H$_2$-poor, that have anomalously large radii, or which orbit less active stars.

Planets with 2 $R_{\oplus}$ < $R$ < 3 $R_{\oplus}$ and orbital period $<$100 d are abundant; these sub-Neptune exoplanets are not well understood. For example, $Kepler$ sub-Neptunes are likely to have deep magma oceans in contact with their atmospheres, but little is known about the effect of the magma on the atmosphere. Here we study this effect using a basic model, assuming that volatiles equilibrate with magma at $T$ $\sim$ 3000 K. For our Fe-Mg-Si-O-H model system, we find that chemical reactions between the magma and the atmosphere and dissolution of volatiles into the magma are both important. Thus, magma matters. For H, most moles go into the magma, so the mass target for both H$_2$ accretion and H$_2$ loss models is weightier than is usually assumed. The known span of magma oxidation states can produce sub-Neptunes that have identical radius but with total volatile masses varying by 20-fold. Thus, planet radius is a proxy for atmospheric composition but not for total volatile content. This redox diversity degeneracy can be broken by measurements of atmosphere mean molecular weight. We emphasise H$_2$ supply by nebula gas, but also consider solid-derived H$_2$O. We find that adding H$_2$O to Fe probably cannot make enough H$_2$ to explain sub-Neptune radii because $>$10$^3$-km thick outgassed atmospheres have high mean molecular weight. The hypothesis of magma-atmosphere equilibration links observables such as atmosphere H$_2$O/H$_2$ ratio to magma FeO content and planet formation processes. Our model's accuracy is limited by the lack of experiments (lab and/or numerical) that are specific to sub-Neptunes; we advocate for such experiments.

Mars has a thin (6 mbar) CO2 atmosphere currently. There is strong evidence for paleolakes and rivers formed by warm climates on Mars, including after 3.5 billion years (Ga) ago, which indicates that a CO2 atmosphere thick enough to permit a warm climate was present at these times. Since Mars no longer has a thick CO2 atmosphere, it must have been lost. One possibility is that Martian CO2 was lost to space. Oxygen escape rates from Mars are high enough to account for loss of a thick CO2 atmosphere, if CO2 was the main source of escaping O. But here, using H isotope ratios, O escape calculations, and quantification of the surface O sinks on Mars, we show for the first time that O escape from Mars after 3.5 Ga must have been predominantly associated with the loss of H2O, not CO2, and therefore it is unlikely that >250 mbar Martian CO2 has been lost to space in the last 3.5 Ga, because such results require highly unfavored O loss scenarios. It is possible that the presence of young rivers and lakes on Mars could be reconciled with limited CO2 loss to space if crater chronologies on Mars are sufficiently incorrect that all apparently young rivers and lakes are actually older than 3.5 Ga, or if climate solutions exist for sustained runoff on Mars with atmospheric CO2 pressure <250 mbar. However, our preferred solution to reconcile the presence of <3.5 Gya rivers and lakes on Mars with the limited potential for CO2 loss to space is a large, as yet undiscovered, geological C sink on Mars.

A candidate explanation for Early Mars rivers is atmospheric warming due to surface release of H$_2$ or CH$_4$ gas. However, it remains unknown how much gas could be released in a single event. We model the CH$_4$ release by one mechanism for rapid release of CH$_4$ from clathrate. By modeling how CH$_4$-clathrate release is affected by changes in Mars' obliquity and atmospheric composition, we find that a large fraction of total outgassing from CH$_4$ clathrate occurs following Mars' first prolonged atmospheric collapse. This atmosphere-collapse-initiated CH$_4$-release mechanism has three stages. (1) Rapid collapse of Early Mars' carbon dioxide atmosphere initiates a slower shift of water ice from high ground to the poles. (2) Upon subsequent CO$_2$-atmosphere re-inflation and CO$_2$-greenhouse warming, low-latitude clathrate decomposes and releases methane gas. (3) Methane can then perturb atmospheric chemistry and surface temperature, until photochemical processes destroy the methane. Within our model, we find that under some circumstances a Titan-like haze layer would be expected to form, consistent with transient deposition of abundant complex abiotic organic matter on the Early Mars surface. We also find that this CH$_4$-release mechanism can warm Early Mars, but special circumstances are required in order to uncork 10$^{17}$ kg of CH$_4$, the minimum needed for strong warming. Specifically, strong warming only occurs when the fraction of the hydrate stability zone that is initially occupied by clathrate exceeds 10%, and when Mars' first prolonged atmospheric collapse occurs for atmospheric pressure > 1 bar.

The James Webb Space Telescope (JWST) will open up the possibility of comprehensively measuring the emission spectra of rocky exoplanets orbiting M dwarfs to detect and characterize their atmospheres. In preparation for this opportunity, we present model atmospheres for three M-dwarf planets particularly amenable to secondary eclipse spectroscopy -- TRAPPIST-1b, GJ 1132b, and LHS 3844b. Using three limiting cases of candidate atmospheric compositions (pure H2O, pure CO2 and solar abundances) we calculate temperature-pressure profiles and emission and reflection spectra in radiative-convective equilibrium, including the effects of a solid surface. We find that the atmospheric radiative transfer is significantly influenced by the cool M-star irradiation; H2O and CO2 absorption bands in the near-infrared are strong enough to absorb a sizeable fraction of the incoming stellar light at low pressures, which leads to temperature inversions in the upper atmosphere. The non-gray band structure of gaseous opacities in the infrared is hereby an important factor. Opacity windows are muted at higher atmospheric temperatures, so we expect temperature inversions to be common only for sufficiently cool planets. We also find that pure CO2 atmospheres exhibit lower overall temperatures and stronger reflection spectra compared to models of the other two compositions. We estimate that for GJ 1132b and LHS 3844b we should be able to distinguish between different atmospheric compositions with JWST. The emission lines from the predicted temperature inversions are currently hard to measure, but high resolution spectroscopy with future Extremely Large Telescopes may be able to detect them.

Transiting planets with radii 2-3 $R_\bigoplus$ are much more numerous than larger planets. We propose that this drop-off is so abrupt because at $R$ $\sim$ 3 $R_\bigoplus$, base-of-atmosphere pressure is high enough for the atmosphere to readily dissolve into magma, and this sequestration acts as a strong brake on further growth. The viability of this idea is demonstrated using a simple model. Our results support extensive magma-atmosphere equilibration on sub-Neptunes, with numerous implications for sub-Neptune formation and atmospheric chemistry.

The upcoming launch of the James Webb Space Telescope (JWST) means that we will soon have the capability to characterize the atmospheres of rocky exoplanets. However, it is still unknown whether such planets orbiting close to M dwarf stars can retain their atmospheres, or whether high-energy irradiation from the star will strip the gaseous envelopes from these objects. We present a new method to detect an atmosphere on a synchronously rotating rocky exoplanet around a K/M dwarf, by using thermal emission during secondary eclipse to infer a high dayside albedo that could only be explained by bright clouds. Based on calculations for plausible surface conditions, we conclude that a high albedo could be unambiguously interpreted as a signal of an atmosphere for planets with substellar temperatures of $T_{sub}=$ 410-1250 K. This range corresponds to equilibrium temperatures of $T_{eq}=$ 300-880 K. We compare the inferred albedos of eight possible planet surface compositions to cloud albedo calculations. We determine that a layer of clouds with optical depths greater than $\tau=0.5$-$7$, would have high enough albedos to be distinguishable from a bare rock surface. This method of detecting an atmosphere on a rocky planet is complementary to existing methods for detecting atmospheres, because it provides a way to detect atmospheres with pressures below 1 bar (e.g. Mars), which are too tenuous to transport significant heat but thick enough to host high-albedo clouds.

Most rocky planets in the galaxy orbit a cool host star, and there is large uncertainty among theoretical models whether these planets can retain an atmosphere. The James Webb Space Telescope (JWST) might be able to settle this question empirically, but most proposals for doing so require large observational effort because they are based on spectroscopy. Here we show that infrared photometry of secondary eclipses could quickly identify "candidate" atmospheres, by searching for rocky planets with atmospheres thick enough that atmospheric heat transport noticeably reduces their dayside thermal emission compared to that of a bare rock. For a planet amenable to atmospheric follow-up, we find that JWST should be able to confidently detect the heat redistribution signal of an O(1) bar atmosphere with one to two eclipses. One to two eclipses is generally much less than the effort needed to infer an atmosphere via transmission or emission spectroscopy. Candidate atmospheres can be further validated via follow-up spectroscopy or phase curves. In addition, because this technique is fast it could enable a first atmospheric survey of rocky exoplanets with JWST. We estimate that the TESS mission will find ~100 planets that are too hot to be habitable but that can be quickly probed via eclipse photometry. Knowing whether hot, rocky planets around M dwarfs have atmospheres is important not only for understanding the evolution of uninhabitable worlds: if atmospheres are common on hot planets, then cooler, potentially habitable planets around M dwarfs are also likely to have atmospheres.

Ancient hydrology is recorded by sedimentary rocks on Mars. The most voluminous sedimentary rocks that formed during Mars' Hesperian period are sulfate-rich rocks, explored by the $Opportunity$ rover from 2004-2012 and soon to be investigated by the $Curiosity$ rover at Gale crater. A leading hypothesis for the origin of these sulfates is that the cations were derived from evaporation of deep-sourced groundwater, as part of a global circulation of groundwater. Global groundwater circulation would imply sustained warm Earthlike conditions on Early Mars. Global circulation of groundwater including infiltration of water initially in equilibrium with Mars' CO$_2$ atmosphere implies subsurface formation of carbonate. We find that the CO$_2$ sequestration implied by the global groundwater hypothesis for the origin of sulfate-rich rocks on Mars is 30-5000 bars if the $Opportunity$ data are representative of Hesperian sulfate-rich rocks, which is so large that (even accounting for volcanic outgassing) it would bury the atmosphere. This disfavors the hypothesis that the cations for Mars' Hesperian sulfates were derived from upwelling of deep sourced groundwater. If, instead, Hesperian sulfate-rich rocks are approximated as pure Mg-sulfate (no Fe), then the CO$_2$ sequestration is 0.3-400 bars. The low end of this range is consistent with the hypothesis that the cations for Mars' Hesperian sulfates were derived from upwelling of deep sourced groundwater. In both cases, carbon sequestration by global groundwater circulation actively works to terminate surface habitability, rather than being a passive marker of warm Earthlike conditions. $Curiosity$ will soon be in a position to discriminate between these two hypotheses. Our work links Mars sulfate cation composition, carbon isotopes, and climate change.

The dynamics of Mars' obliquity are believed to be chaotic, and the historical ~3.5 Gyr (late-Hesperian onward) obliquity probability density function (PDF) is high uncertain and cannot be inferred from direct simulation alone. Obliquity is also a strong control on post-Noachian Martian climate, enhancing the potential for equatorial ice/snow melting and runoff at high obliquities (> 40{\deg}) and enhancing the potential for desiccation of deep aquifers at low obliquities (< 25{\deg}). We developed a new technique using the orientations of elliptical craters to constrain the true late-Hesperian-onward obliquity PDF. To do so, we developed a forward model of the effect of obliquity on elliptic crater orientations using ensembles of simulated Mars impactors and ~3.5 Gyr-long Mars obliquity simulations. In our model, the inclinations and speeds of Mars crossing objects bias the preferred orientation of elliptic craters which are formed by low-angle impacts. Comparison of our simulation predictions with a validated database of elliptic crater orientations allowed us to invert for best-fitting obliquity history. We found that since the onset of the late-Hesperian, Mars' mean obliquity was likely low, between ~10{\deg} and ~30{\deg}, and the fraction of time spent at high obliquities > 40{\deg} was likely < 20%.

The recent discovery of a staggering diversity of planets beyond the Solar System has brought with it a greatly expanded search space for habitable worlds. The Kepler exoplanet survey has revealed that most planets in our interstellar neighborhood are larger than Earth and smaller than Neptune. Collectively termed super-Earths and mini-Neptunes, some of these planets may have the conditions to support liquid water oceans, and thus Earth-like biology, despite differing in many ways from our own planet. In addition to their quantitative abundance, super-Earths are relatively large and are thus more easily detected than true Earth twins. As a result, super-Earths represent a uniquely powerful opportunity to discover and explore a panoply of fascinating and potentially habitable planets in 2020 - 2030 and beyond.

Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets.

The habitable zone (HZ) is the region around a star(s) where standing bodies of water could exist on the surface of a rocky planet. The classical HZ definition makes a number of assumptions common to the Earth, including assuming that the most important greenhouse gases for habitable planets are CO2 and H2O, habitable planets orbit main-sequence stars, and that the carbonate-silicate cycle is a universal process on potentially habitable planets. Here, we discuss these and other predictions for the habitable zone and the observations that are needed to test them. We also, for the first time, argue why A-stars may be interesting HZ prospects. Instead of relying on unverified extrapolations from our Earth, we argue that future habitability studies require first principles approaches where temporal, spatial, physical, chemical, and biological systems are dynamically coupled. We also suggest that next-generation missions are only the beginning of a much more data-filled era in the not-too-distant future, when possibly hundreds to thousands of HZ planets will yield the statistical data we need to go beyond just finding habitable zone planets to actually determining which ones are most likely to exhibit life.

Early Mars climate research has well-defined goals (Mars Exploration Program Analysis Group 2018). Achieving these goals requires geologists and climate modelers to coordinate. Coordination is easier if results are expressed in terms of well-defined parameters. Key parameters include the following quantitative geologic constraints. (1) Cumulative post-3.4 Ga precipitation-sourced water runoff in some places exceeded 1 km column. (2) There is no single Early Mars climate problem: the traces of $\ge$2 river-forming periods are seen. Relative to rivers that formed earlier in Mars history, rivers that formed later in Mars history are found preferentially at lower elevations, and show a stronger dependence on latitude. (3) The duration of the longest individual river-forming climate was >(10$^2$-10$^3$) yr, based on paleolake hydrology. (4) Peak runoff production was >0.1 mm/hr. However, (5) peak runoff production was intermittent, sustained (in a given catchment) for only <10% of the duration of river-forming climates. (6) The cumulative number of wet years during the valley-network-forming period was >10$^5$ yr. (7) Post-Noachian light-toned, layered sedimentary rocks took >10$^7$ yr to accumulate. However, (8) an "average" place on Mars saw water for <10$^7$ yr after the Noachian, suggesting that the river-forming climates were interspersed with long globally-dry intervals. (9) Geologic proxies for Early Mars atmospheric pressure indicate pressure was not less than 0.012 bar but not much more than 1 bar. A truth table of these geologic constraints versus currently published climate models shows that the late persistence of river-forming climates, combined with the long duration of individual lake forming climates, is a challenge for most models.

Lakes existed on Mars later than 3.6 billion years ago, according to sedimentary evidence for deltaic deposition. The observed fluvio-lacustrine deposits suggest that individual lake-forming climates persisted for at least several thousand years (assuming dilute flow). But the lake watersheds' little weathered soils indicate a largely dry climate history, with intermittent runoff events. Here we show that these observational constraints, while inconsistent with many previously-proposed triggers for lake-forming climates, are consistent with a methane burst scenario. In this scenario, chaotic transitions in mean obliquity drive latitudinal shifts in temperature and ice loading that destabilize methane clathrate. Using numerical simulations, we find that outgassed methane can build up to atmospheric levels sufficient for lake forming climates, for past clathrate hydrate stability zone occupancy fractions >0.04. Such occupancy fractions are consistent with methane production by water-rock reactions due to hydrothermal circulation on early Mars. We further estimate that photochemical destruction of atmospheric methane curtails the duration of individual lake-forming climates to less than a million years, consistent with observations. We conclude that methane bursts represent a potential pathway for intermittent excursions to a warm, wet climate state on early Mars.

Large alluvial fan deposits on Mars record relatively recent habitable surface conditions ($\lesssim$3.5 Ga, Late Hesperian - Amazonian). We find net sedimentation rate <(4-8) {\mu}m/yr in the alluvial-fan deposits, using the frequency of craters that are interbedded with alluvial-fan deposits as a fluvial-process chronometer. Considering only the observed interbedded craters sets a lower bound of >20 Myr on the total time interval spanned by alluvial-fan aggradation, >10^3-fold longer than previous lower limits. A more realistic approach that corrects for craters fully entombed in the fan deposits raises the lower bound to >(100-300) Myr. Several factors not included in our calculations would further increase the lower bound. The lower bound rules out fan-formation by a brief climate anomaly. Therefore, during the Late Hesperian - Amazonian on Mars, persistent or repeated processes permitted habitable surface conditions.

Many habitable zone exoplanets are expected to form with water mass fractions higher than that of the Earth. For rocky exoplanets with 10-1000x Earth's H2O but without H2, we model the multi-Gyr evolution of ocean temperature and chemistry, taking into account C partitioning, high-pressure ice phases, and atmosphere-lithosphere exchange. Within our model, for Sun-like stars, we find that: (1)~the duration of habitable surface water is strongly affected by ocean chemistry; (2)~possible ocean pH spans a wide range; (3)~surprisingly, many waterworlds retain habitable surface water for >1 Gyr, and (contrary to previous claims) this longevity does not necessarily involve geochemical cycling. The key to this cycle-independent planetary habitability is that C exchange between the convecting mantle and the water ocean is curtailed by seafloor pressure on waterworlds, so the planet is stuck with the ocean mass and ocean cations that it acquires during the first 1% of its history. In our model, the sum of positive charges leached from the planetary crust by early water-rock interactions is - coincidentally - often within an order of magnitude of the early-acquired atmosphere+ocean inorganic C inventory overlaps. As a result, pCO2 is frequently in the "sweet spot" (0.2-20 bar) for which the range of semimajor axis that permits surface liquid water is about as wide as it can be. Because the width of the HZ in semimajor axis defines (for Sun-like stars) the maximum possible time span of surface habitability, this effect allows for Gyr of habitability as the star brightens. We illustrate our findings by using the output of an ensemble of N-body simulations as input to our waterworld evolution code. Thus (for the first time in an end-to-end calculation) we show that chance variation of initial conditions, with no need for geochemical cycling, can yield multi-Gyr surface habitability on waterworlds.

Post-Noachian Martian paleochannels indicate the existence of liquid water on the surface of Mars after about 3.5 Gya (Irwin et al., 2015; Palucis et al., 2016). In order to explore the effects of variations in CO$_{2}$ partial pressure and obliquity on the possibility of surface water, we created a zero-dimensional surface energy balance model. We combine this model with physically consistent orbital histories to track conditions over the last 3.5 Gyr of Martian history. We find that melting is allowed for atmospheric pressures corresponding to exponential loss rates of $dP/dt \propto t^{-3.73}$ or faster, but this rate is within $0.5 \sigma$ of the rate calculated from initial measurements made by the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission, if we assume all the escaping oxygen measured by MAVEN comes from atmospheric CO$_{2}$ (Lillis et al., 2017; Tu et al., 2015). Melting at this loss rate matches selected key geologic constraints on the formation of Hesperian river networks, assuming optimal melt conditions during the warmest part of each Mars year (Irwin et al., 2015; Stopar et al., 2006; Kite et al., 2017a,b). The atmospheric pressure has a larger effect on the surface energy than changes in Mars's mean obliquity. These results show that initial measurements of atmosphere loss by MAVEN are consistent with atmospheric loss being the dominant process that switched Mars from a melt-permitting to a melt-absent climate (Jakosky et al., 2017), but non-CO$_{2}$ warming will be required if $<2$ Gya paleochannels are confirmed, or if most of the escaping oxygen measured by MAVEN comes from H$_{2}$O.

Recent discoveries imply that Early Mars was habitable for life-as-we-know-it; that Enceladus might be habitable; and that many stars have Earth-sized exoplanets whose insolation favors surface liquid water. These exciting discoveries make it more likely that spacecraft now under construction - Mars 2020, ExoMars rover, JWST, Europa Clipper - will find habitable, or formerly habitable, environments. Did these environments see life? Given finite resources (\$10bn/decade for the US ), how could we best test the hypothesis of a second origin of life? Here, we first state the case for and against flying life detection missions soon. Next, we assume that life detection missions will happen soon, and propose a framework for comparing the value of different life detection missions: Scientific value = (Reach x grasp x certainty x payoff) / \$ After discussing each term in this framework, we conclude that scientific value is maximized if life detection missions are flown as hypothesis tests. With hypothesis testing, even a nondetection is scientifically valuable.

Most of Mars' ancient sedimentary rocks by volume are in wind-eroded sedimentary mounds, but the connections between mound form and wind erosion are unclear. We perform mesoscale simulations of different crater and mound morphologies to understand the formation of sedimentary mounds. As crater depth increases, slope winds produce increased erosion near the base of the crater wall, forming mounds. Peak erosion rates occur when the crater depth is ~2 km. Mound evolution depends on the size of the host crater. In smaller craters mounds preferentially erode at the top, becoming more squat, while in larger craters mounds become steeper-sided. This agrees with observations where smaller craters tend to have proportionally shorter mounds, and larger craters have mounds encircled by moats. If a large-scale sedimentary layer blankets a crater, then as the layer recedes across the crater it will erode more towards the edges of the crater, resulting in a crescent-shaped moat. When a 160 km diameter mound-hosting crater is subject to a prevailing wind, the surface wind stress is stronger on the leeward side than on the windward side. This results in the center of the mound appearing to `march upwind' over time, and forming a `bat-wing' shape, as is observed for Mt. Sharp in Gale crater.

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 around the time of fluvial activity. Here we use high-resolution orthophotos and digital terrain models from the Mars Reconnaissance Orbiter to identify ancient craters 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$\pm$0.1 bar, rising to 1.9$\pm$0.2 bar if rimmed circular mesas - interpreted to be erosionally-resistant fills of floors of impact craters - are excluded. We assume target properties appropriate for desert alluvium: if sediment had rock-like rock-mass strength similar to bedrock at the time of impact, the 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 the warm and wet CO2/H2O greenhouse is ruled out, and long-term average temperatures were most likely below freezing.

Unraveling the stratigraphic record is the key to understanding ancient climate and past climate changes on Mars. River deposits when placed in stratigraphic order could constrain the number, magnitudes, and durations of the wettest climates in Mars history. We establish the stratigraphic context of river deposits in Aeolis Dorsa sedimentary basin, 10E of Gale crater. Here, wind has exhumed a stratigraphic section of >=4 unconformity-bounded sedimentary rock packages, recording >=3 distinct episodes of surface runoff. Early deposits (>700m thick) are embayed by river deposits (>400m), which are in turn unconformably draped by fan-shaped deposits (<100m) which we interpret as alluvial fans. Yardang-forming deposits (>900 m) unconformably drape all previous deposits. River deposits embay a dissected sedimentary-rock landscape, and comprise >=2 distinguishable units. The total interval spanned by river deposits is >(1x10^6-2x10^7) yr; more if we include alluvial-fan deposits. Alluvial-fan deposits unconformably postdate thrust faults which crosscut river deposits. We infer a relatively dry interval of >4x10^7 yr after river deposits formed and before fan-shaped deposits formed. The time gap between the end of river deposition and the onset of yardang-forming deposits is constrained to >10^8 yr by the density of impact craters embedded at the unconformity. We correlate yardang-forming deposits to the upper layers of Gale crater's mound (Mt. Sharp/Aeolis Mons), and fan-shaped deposits to Peace Vallis fan. Alternations between periods of low vs. high mean obliquity may have modulated erosion-deposition cycling in Aeolis. This is consistent with results from an ensemble of simulations of Solar System orbital evolution and the resulting history of Mars obliquity. Almost all simulations yield intervals of continuously low mean Mars obliquity that are long enough to match our unconformity data.

Chloride-bearing deposits on Mars record high-elevation lakes during the waning stages of Mars' wet era (mid-Noachian to late Hesperian). The water source pathways, seasonality, salinity, depth, lifetime, and paleoclimatic drivers of these widespread lakes are all unknown. Here we combine reaction-transport modeling, orbital spectroscopy, and new volume estimates from high-resolution digital terrain models, in order to constrain the hydrologic boundary conditions for forming the chlorides. Considering a T = 0 degrees C system, we find: (1) individual lakes were >100 m deep and lasted decades or longer; (2) if volcanic degassing was the source of chlorine, then the water-to-rock ratio or the total water volume were probably low, consistent with brief excursions above the melting point and/or arid climate; (3) if the chlorine source was igneous chlorapatite, then Cl-leaching events would require a (cumulative) time of >10 yr at the melting point; (4) Cl masses, divided by catchment area, give column densities 0.1 - 50 kg Cl/m^2, and these column densities bracket the expected chlorapatite-Cl content for a seasonally-warm active layer. Deep groundwater was not required. Taken together, our results are consistent with Mars having a usually cold, horizontally segregated hydrosphere by the time chlorides formed.

River deposits are one of the main lines of evidence that tell us that Mars once had a climate different from today, and so changes in river deposits with time tell us something about how Mars climate changed with time. In this study, we focus in on one sedimentary basin - Aeolis Dorsa - which contains an exceptionally high number of exceptionally well-preserved river deposits that appear to have formed over an interval of $>$0.5 Myr. We use changes in the river deposits' scale with stratigraphic elevation as a proxy for changes in river paleodischarge. Meander wavelengths tighten upwards and channel widths narrow upwards, and there is some evidence for a return to wide large-wavelength channels higher in the stratigraphy. Meander wavelength and channel width covary with stratigraphic elevation. The factor of 1.5-2 variations in paleochannel dimensions with stratigraphic elevation correspond to $\sim$2.6-fold variability in bank-forming discharge (using standard wavelength-discharge scalings and width-discharge scalings). Taken together with evidence from a marker bed for discharge variability at $\sim$10m stratigraphic distances, the variation in the scale of river deposits indicates that bank-forming discharge varied at both 10m stratigraphic (10$^2$ - 10$^6$ yr) and $\sim$100 m stratigraphic (10$^3$ - 10$^9$ yr) scales. Because these variations are correlated across the basin, they record a change in basin-scale forcing, rather than smaller-scale internal feedbacks. Changing sediment input leading to a change in characteristic slopes and/or drainage area could be responsible, and another possibility is changing climate ($\pm$50 W/m$^2$ in peak energy available for snow/ice melt).

We present a new database of $>$300 layer-orientations from sedimentary mounds on Mars. These layer orientations, together with draped landslides, and draping of rocks over differentially-eroded paleo-domes, indicate that for the stratigraphically-uppermost $\sim$1 km, the mounds formed by the accretion of draping strata in a mound-shape. The layer-orientation data further suggest that layers lower down in the stratigraphy also formed by the accretion of draping strata in a mound-shape. The data are consistent with terrain-influenced wind erosion, but inconsistent with tilting by flexure, differential compaction over basement, or viscoelastic rebound. We use a simple landscape evolution model to show how the erosion and deposition of mound strata can be modulated by shifts in obliquity. The model is driven by multi-Gyr calculations of Mars' chaotic obliquity and a parameterization of terrain-influenced wind erosion that is derived from mesoscale modeling. Our results suggest that mound-spanning unconformities with kilometers of relief emerge as the result of chaotic obliquity shifts. Our results support the interpretation that Mars' rocks record intermittent liquid-water runoff during a $>$10$^8$-yr interval of sedimentary rock emplacement.

Small-crater counts on Mars light-toned sedimentary rock are often inconsistent with any isochron; these data are usually plotted then ignored. We show (using an 18-HiRISE-image, >10^4 crater dataset) that these non-isochron crater counts are often well-fit by a model where crater production is balanced by crater obliteration via steady exhumation. For these regions, we fit erosion rates. We infer that Mars light-toned sedimentary rocks typically erode at ~10^2 nm/yr, when averaged over 10 km^2 scales and 10^7-10^8 yr timescales. Crater-based erosion-rate determination is consistent with independent techniques, but can be applied to nearly all light-toned sedimentary rocks on Mars. Erosion is swift enough that radiolysis cannot destroy complex organic matter at some locations (e.g. paleolake deposits at SW Melas), but radiolysis is a severe problem at other locations (e.g. Oxia Planum). The data suggest that the relief of the Valles Marineris mounds is currently being reduced by wind erosion, and that dust production on Mars <3 Gya greatly exceeds the modern reservoir of mobile dust.

We provide estimates of atmospheric pressure and surface composition on short-period rocky exoplanets with dayside magma pools and silicate vapor atmospheres. Atmospheric pressure tends toward vapor-pressure equilibrium with surface magma, and magma-surface composition is set by the competing effects of fractional vaporization and surface-interior exchange. We use basic models to show how surface-interior exchange is controlled by the planet's temperature, mass, and initial composition. We assume that mantle rock undergoes bulk melting to form the magma pool, and that winds flow radially away from the substellar point. With these assumptions, we find that: (1) atmosphere-interior exchange is fast when the planet's bulk-silicate FeO concentration is low, and slow when FeO concentration is high; (2) magma pools are compositionally well-mixed for substellar temperatures $\lesssim$ 2400 K, but compositionally variegated and rapidly variable for substellar temperatures $\gtrsim$ 2400 K; (3) currents within the magma pool tend to cool the top of the solid mantle ("tectonic refrigeration"); (4) contrary to earlier work, many magma planets have time-variable surface compositions.

Spacecraft observations suggest that the plumes of Saturn's moon Enceladus draw water from a subsurface ocean, but the sustainability of conduits linking ocean and surface is not understood. Observations show sustained (though tidally modulated) fissure eruptions throughout each orbit, and since the 2005 discovery of the plumes. Peak plume flux lags peak tidal extension by $\sim$1 radian, suggestive of resonance. Here we show that a model of the tiger stripes as tidally-flexed slots that puncture the ice shell can simultaneously explain the persistence of the eruptions through the tidal cycle, the phase lag, and the total power output of the tiger stripe terrain, while suggesting that the eruptions are maintained over geological timescales. The delay associated with flushing and refilling of \emph{O}(1) m-wide slots with ocean water causes erupted flux to lag tidal forcing and helps to buttress slots against closure, while tidally pumped in-slot flow leads to heating and mechanical disruption that staves off slot freeze-out. Much narrower and much wider slots cannot be sustained. In the presence of long-lived slots, the 10$^6$-yr average power output of the tiger stripes is buffered by a feedback between ice melt-back and subsidence to \emph{O}(10$^{10}$) W, which is similar to the observed power output, suggesting long-term stability. Turbulent dissipation makes testable predictions for the final flybys of Enceladus by the \emph{Cassini} spacecraft. Our model shows how open connections to an ocean can be reconciled with, and sustain, long-lived eruptions. Turbulent dissipation in long-lived slots helps maintain the ocean against freezing, maintains access by future Enceladus missions to ocean materials, and is plausibly the major energy source for tiger stripe activity.

The impactor flux early in Mars history was much higher than today, so sedimentary sequences include many buried craters. In combination with models for the impactor flux, observations of the number of buried craters can constrain sedimentation rates. Using the frequency of crater-river interactions, we find net sedimentation rate \lesssim 20-300 {\mu}m/yr at Aeolis Dorsa. This sets a lower bound of 1-15 Myr on the total interval spanned by fluvial activity around the Noachian-Hesperian transition. We predict that Gale Crater's mound (Aeolis Mons) took at least 10-100 Myr to accumulate, which is testable by the Mars Science Laboratory.

Ancient sediments provide archives of climate and habitability on Mars. Gale Crater, the landing site for the Mars Science Laboratory (MSL), hosts a 5 km high sedimentary mound. Hypotheses for mound formation include evaporitic, lacustrine, fluviodeltaic, and aeolian processes, but the origin and original extent of Gale's mound is unknown. Here we show new measurements of sedimentary strata within the mound that indicate ~3 degree outward dips oriented radially away from the mound center, inconsistent with the first three hypotheses. Moreover, although mounds are widely considered to be erosional remnants of a once crater-filling unit, we find that the Gale mound's current form is close to its maximal extent. Instead we propose that the mound's structure, stratigraphy, and current shape can be explained by growth in place near the center of the crater mediated by wind-topography feedbacks. Our model shows how sediment can initially accrete near the crater center far from crater-wall katabatic winds, until the increasing relief of the resulting mound generates mound-flank slope-winds strong enough to erode the mound. Our results indicate mound formation by airfall-dominated deposition with a limited role for lacustrine and fluvial activity, and potentially limited organic carbon preservation. Morphodynamic feedbacks between wind and topography are widely applicable to a range of sedimentary mounds and ice mounds across the Martian surface, and possibly other planets.

A model for the formation and distribution of sedimentary rocks on Mars is proposed. The rate-limiting step is supply of liquid water from seasonal melting of snow or ice. The model is run for a O(10^2) mbar pure CO2 atmosphere, dusty snow, and solar luminosity reduced by 23%. For these conditions snow only melts near the equator, and only when obliquity >40 degrees, eccentricity >0.12, and perihelion occurs near equinox. These requirements for melting are satisfied by 0.01-20% of the probability distribution of Mars' past spin-orbit parameters. Total melt production is sufficient to account for aqueous alteration of the sedimentary rocks. The pattern of seasonal snowmelt is integrated over all spin-orbit parameters and compared to the observed distribution of sedimentary rocks. The global distribution of snowmelt has maxima in Valles Marineris, Meridiani Planum and Gale Crater. These correspond to maxima in the sedimentary-rock distribution. Higher pressures and especially higher temperatures lead to melting over a broader range of spin-orbit parameters. The pattern of sedimentary rocks on Mars is most consistent with a Mars paleoclimate that only rarely produced enough meltwater to precipitate aqueous cements and indurate sediment. The results suggest intermittency of snowmelt and long globally-dry intervals, unfavorable for past life on Mars. This model makes testable predictions for the Mars Science Laboratory rover at Gale Crater. Gale Crater is predicted to be a hemispheric maximum for snowmelt on Mars.

We report here on the discovery of stellar occultations, observed with Kepler, that recur periodically at 15.685 hour intervals, but which vary in depth from a maximum of 1.3% to a minimum that can be less than 0.2%. The star that is apparently being occulted is KIC 12557548, a K dwarf with T_eff = 4400 K and V = 16. Because the eclipse depths are highly variable, they cannot be due solely to transits of a single planet with a fixed size. We discuss but dismiss a scenario involving a binary giant planet whose mutual orbit plane precesses, bringing one of the planets into and out of a grazing transit. We also briefly consider an eclipsing binary, that either orbits KIC 12557548 in a hierarchical triple configuration or is nearby on the sky, but we find such a scenario inadequate to reproduce the observations. We come down in favor of an explanation that involves macroscopic particles escaping the atmosphere of a slowly disintegrating planet not much larger than Mercury. The particles could take the form of micron-sized pyroxene or aluminum oxide dust grains. The planetary surface is hot enough to sublimate and create a high-Z atmosphere; this atmosphere may be loaded with dust via cloud condensation or explosive volcanism. Atmospheric gas escapes the planet via a Parker-type thermal wind, dragging dust grains with it. We infer a mass loss rate from the observations of order 1 M_E/Gyr, with a dust-to-gas ratio possibly of order unity. For our fiducial 0.1 M_E planet, the evaporation timescale may be ~0.2 Gyr. Smaller mass planets are disfavored because they evaporate still more quickly, as are larger mass planets because they have surface gravities too strong to sustain outflows with the requisite mass-loss rates. The occultation profile evinces an ingress-egress asymmetry that could reflect a comet-like dust tail trailing the planet; we present simulations of such a tail.

Feedbacks that can destabilize the climates of synchronously-rotating rocky planets may arise on planets with strong day-night surface temperature contrasts. Earth-like habitable-zone (HZ) planets maintain stable surface liquid water over geological time. This requires equilibrium between the temperature-dependent rate of greenhouse-gas consumption by weathering,and greenhouse-gas resupply by other processes. Detected small-radius exoplanets, and anticipated M-dwarf HZ rocky planets, are expected to be tidally locked. We investigate two feedbacks that can destabilize climate on tidally-locked planets. (1) If small changes in pressure alter the temperature distribution across a planet's surface such that the weathering rate increases when the pressure decreases, a positive feedback occurs involving increasing weathering rate near the substellar point, decreasing pressure, and increasing substellar surface temperature. (2) When decreases in pressure increase the surface area above the melting point (through reduced advective cooling of the substellar point), and the corresponding increase in volume of liquid causes net dissolution of the atmosphere, a further decrease in pressure occurs. We use an idealized energy balance model to map out the conditions under which these instabilities may occur. The weathering runaway can shrink the habitable zone, and cause geologically rapid 10^3-fold pressure shifts within the HZ. Mars may have undergone a weathering runaway in the past. Substellar dissolution is usually a negative feedback or weak positive feedback on changes in pressure. Both instabilities are suppressed if the atmosphere has a high radiative efficiency. Our results are most relevant for atmospheres that are thin and have low greenhouse-gas radiative efficiency. These results identify a new pathway by which HZ planets can undergo rapid climate shifts and become uninhabitable.

Channel networks on the plateau adjacent to Juventae Chasma have the highest drainage densities reported on Mars.We model frozen precipitation on the Juventae plateau,finding that the trigger for forming these channel networks could have been ephemeral lakeshore precipitation,and that they do not require past temperatures higher than today.If short-lived and localized events explain some dendritic channel networks on Mars, this would weaken the link between dendritic valley networks and surface climate conditions that could sustain life. Our analysis uses MRAMS simulations and HiRISE DTMs.We model localized weather systems driven by water vapor release from ephemeral lakes during outflow channel formation.At Juventae Chasma,mean snowfall reaches a maximum of 0.9mm/hr water equivalent on the SW rim of the chasm.Radiative effects of the thick cloud cover raise maximum (minimum, mean) plateau surface temperatures by up to 24K(9K, 17K)locally.The key result is that the area of maximum modeled precipitation shows a striking correspondence to the mapped Juventae plateau channel networks.Three independent methods show this fit is unlikely to be due to chance.We use a snowpack energy balance model to show that if the snow has the albedo of dust(0.28), and for a solar luminosity of 0.8($\equiv$3.0Gya), then if the atmospheric greenhouse effect is unchanged from(6K warmer than)today only 0.4%(21%)of lake-induced precipitation events produce snowpack that undergoes melting.However, warming from associated dense cloud cover would allow melting over a wider range of conditions.In these localized precipitation scenarios, global temperatures need not be higher than today, and the rest of the planet remains dry.

We use the Mars Regional Atmospheric Modeling System (MRAMS) to simulate lake storms on Mars, finding that intense localized precipitation will occur for lake size >=10^3 km^2. Mars has a low-density atmosphere, so deep convection can be triggered by small amounts of latent heat release. In our reference simulation, the buoyant plume lifts vapor above condensation level, forming a 20km-high optically-thick cloud. Ice grains grow to 200 microns radius and fall near (or in) the lake at mean rates up to 1.5 mm/hr water equivalent (maximum rates up to 6 mm/hr water equivalent). Because atmospheric temperatures outside the surface layer are always well below 273K, supersaturation and condensation begin at low altitudes above lakes on Mars. In contrast to Earth lake-effect storms, lake storms on Mars involve continuous precipitation, and their vertical velocities and plume heights exceed those of tropical thunderstorms on Earth. Convection does not reach above the planetary boundary layer for lakes <<10^3 km^2 or for atmospheric pressure >O(10^2) mbar. Instead, vapor is advected downwind with little cloud formation. Precipitation occurs as snow, and the daytime radiative forcing at the land surface due to plume vapor and storm clouds is too small to melt snow directly (<+10 W/m^2). However, if orbital conditions are favorable, then the snow may be seasonally unstable to melting and produce runoff to form channels. We calculate the probability of melting by running thermal models over all possible orbital conditions and weighting their outcomes by probabilities given by Laskar et al., 2004. We determine that for an equatorial vapor source, sunlight 15% fainter than at present, and snowpack with albedo 0.28 (0.35), melting may occur with 4%(0.1%) probability. This rises to 56%(12%) if the ancient greenhouse effect was modestly (6K) greater than today.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

We provide estimates of volcanism versus time for planets with Earth-like composition and masses from 0.25 to 25 times Earth, as a step toward predicting atmospheric mass on extrasolar rocky planets. Volcanism requires melting of the silicate mantle. We use a thermal evolution model, calibrated against Earth, in combination with standard melting models, to explore the dependence of convection-driven decompression mantle melting on planet mass. Here we show that (1) volcanism is likely to proceed on massive planets with plate tectonics over the main-sequence lifetime of the parent star; (2) crustal thickness (and melting rate normalized to planet mass) is weakly dependent on planet mass; (3) stagnant lid planets live fast (they have higher rates of melting than their plate tectonic counterparts early in their thermal evolution) but die young (melting shuts down after a few Gyr); (4) plate tectonics may not operate on high mass planets because of the production of buoyant crust which is difficult to subduct; and (5) melting is necessary but insufficient for efficient volcanic degassing - volatiles partition into the earliest, deepest melts, which may be denser than the residue and sink to the base of the mantle on young, massive planets. Magma must also crystallize at or near the surface, and the pressure of overlying volatiles must be fairly low, if volatiles are to reach the surface. If volcanism is detected in the Tau Ceti system, and tidal forcing can be shown to be weak, this would be evidence for plate tectonics.

Following the optical imaging of the exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhaut's debris disk is gravitationally shaped by an interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass < 3 Jupiter masses. If the belt and planet orbits are apsidally aligned, our model predicts a planet mass of 0.5 Jupiter masses. The inner edge of the debris disk at 133 AU lies at the periphery of Fom b's chaotic zone, and the mean disk eccentricity of 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planet's chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of 100 Myr, and model them separately from their dust grain progeny; the latter's orbits are strongly affected by radiation pressure and their lifetimes are limited to 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fom b's nominal space velocity does not bear this out, but the astrometric uncertainties may be large. If the apsidal misalignment is real, our upper mass limit of 3 Jupiter masses still holds. The belt contains at least 3 Earth masses of solids that are grinding down to dust. Such a large mass in solids is consistent with Fom b having formed in situ.

The areal centroids of the youngest polar deposits on Mars are offset from those of adjacent paleopolar deposits by 5-10 degrees. We test the hypothesis that the offset is the result of true polar wander (TPW), the motion of the solid surface with respect to the spin axis, caused by a mass redistribution within or on the surface of Mars. In particular, we consider TPW driven by late-stage volcanism during the late Hesperian to Amazonian. There is observational and qualitative support for this hypothesis: in both North and South, observed offsets lie close to a great circle 90 degrees from Tharsis, as expected for polar wander after Tharsis formed. We calculate the magnitude and direction of TPW produced by mapped late-stage lavas for a range of lithospheric thicknesses, lava thicknesses, eruption histories, and prior polar wander events. If Tharsis formed close to the equator, the stabilizing effect of a fossil rotational bulge located close to the equator leads to predicted TPW of <2 degrees, too small to account for observed offsets. If, however, Tharsis formed far from the equator, late-stage TPW driven by low-latitude, late-stage volcanism would be 6-33 degrees, similar to that inferred from the location of paleopolar deposits. 4.4+/-1.3x10^19 kg of young erupted lava can account for the offset of the Dorsa Argentea Formation from the present-day south rotation pole. This mass is consistent with prior mapping-based estimates and would imply a mass release of CO2 by volcanic degassing similar to that in the atmosphere at the present time. The South Polar Layered Deposits are offset from the spin axis in the opposite sense to the other paleopolar deposits. This can be explained by an additional contribution from a plume beneath Elysium. We conclude with a list of observational tests of the TPW hypothesis.

Fomalhaut is a bright star 7.7 parsecs (25 light years) from Earth that harbors a belt of cold dust with a structure consistent with gravitational sculpting by an orbiting planet. Here, we present optical observations of an exoplanet candidate, Fomalhaut b. In the plane of the belt, Fomalhaut b lies approximately 119 astronomical units (AU) from the star and 18 AU from the belt, matching predictions. We detect counterclockwise orbital motion using Hubble Space Telescope observations separated by 1.73 years. Dynamical models of the interaction between the planet and the belt indicate that the planet's mass is at most three times that of Jupiter for the belt to avoid gravitational disruption. The flux detected at 800 nm is also consistent with that of a planet with mass no greater than a few times that of Jupiter. The brightness at 600 nm and the lack of detection at longer wavelengths suggest that the detected flux may include starlight reflected off a circumplanetary disk, with dimension comparable to the orbits of the Galilean satellites. We also observed variability of unknown origin at 600 nm.