Search

The increased brightness temperature of young rocky protoplanets during their magma ocean epoch makes them potentially amenable to atmospheric characterization to distances from the solar system far greater than thermally equilibrated terrestrial exoplanets, offering observational opportunities for unique insights into the origin of secondary atmospheres and the near surface conditions of prebiotic environments. The Large Interferometer For Exoplanets (LIFE) mission will employ a space-based mid-infrared nulling interferometer to directly measure the thermal emission of terrestrial exoplanets. Here, we seek to assess the capabilities of various instrumental design choices of the LIFE mission concept for the detection of cooling protoplanets with transient high-temperature magma ocean atmospheres, in young stellar associations in particular. Using the LIFE mission instrument simulator (LIFEsim) we assess how specific instrumental parameters and design choices, such as wavelength coverage, aperture diameter, and photon throughput, facilitate or disadvantage the detection of protoplanets. We focus on the observational sensitivities of distance to the observed planetary system, protoplanet brightness temperature using a blackbody assumption, and orbital distance of the potential protoplanets around both G- and M-dwarf stars. Our simulations suggest that LIFE will be able to detect (S/N $\geq$ 7) hot protoplanets in young stellar associations up to distances of $\approx$100 pc from the solar system for reasonable integration times (up to $\sim$hours). Detection of an Earth-sized protoplanet orbiting a solar-sized host star at 1 AU requires less than 30 minutes of integration time. M-dwarfs generally need shorter integration times. The contribution from wavelength regions $<$6 $\mu$m is important for decreasing the detection threshold and discriminating emission temperatures., Comment: 18 pages, 19 figures; accepted for publication in A&A

The equilibrium rotation rate of a planet is determined by the sum of torques acting on its solid body. For planets with atmospheres, the dominant torques are usually the gravitational tide, which acts to slow the planet's rotation rate, and the atmospheric thermal tide, which acts to spin up the planet. Previous work demonstrated that rocky planets with thick atmospheres may produce strong enough thermal tides to avoid tidal locking, but a study of how the strength of the thermal tide depends on atmospheric properties has not been done. In this work, we use a combination of simulations from a global climate model and analytic theory to explore how the thermal tide depends on the shortwave and longwave optical depth of the atmosphere, the surface pressure, and the absorbed stellar radiation. We find that for planets in the habitable zones of M stars only high-pressure but low-opacity atmospheres permit asynchronous rotation owing to the weakening of the thermal tide at high longwave and shortwave optical depths. We conclude that asynchronous rotation may be very unlikely around low-mass stars, which may limit the potential habitability of planets around M stars.

We investigate the impact of nightside cloud formation on the observable night-day contrast of tidally-locked terrestrial planet atmospheres. We demonstrate that, in the case where the planetary dayside is only 10s of Kelvin hotter than the planetary nightside, the presence of optically thick nightside clouds can lead to observations that mimic a planet without an atmosphere, despite the planet actually hosting a significant (10 bar) atmosphere. The scenario presented in this work requires a level of intrinsic atmospheric day/night temperature contrast such that the nightside can form clouds while the dayside is too hot for cloud formation to occur. This scenario is most likely for hotter terrestrials and terrestrials with low volatile inventories. We thus note that a substantial dayside/nightside temperature difference alone does not robustly indicate that a planet does not host an atmosphere and additional observations and modeling are essential for characterization. We further discuss several avenues for future study to improve our understanding of the terrestrial planets and how best to characterize them with JWST., Comment: Accepted to ApJL, 12 pages, 9 figures

Atmospheric retrievals are widely used to constrain exoplanet properties from observed spectra. We investigate how the common nonphysical retrieval assumptions of vertically constant molecule abundances and cloud-free atmospheres affect our characterization of an exo-Earth (an Earth-twin orbiting a Sun-like star). Specifically, we use a state-of-the-art retrieval framework to explore how assumptions for the $\mathrm{H_2O}$ profile and clouds affect retrievals. In a first step, we validate different retrieval models on a low-noise simulated 1D mid-infrared (MIR) spectrum of Earth. Thereafter, we study how these assumptions affect the characterization of Earth with the Large Interferometer For Exoplanets (LIFE). We run retrievals on LIFE mock observations based on real disk-integrated MIR Earth spectra. The performance of different retrieval models is benchmarked against ground truths derived from remote sensing data. We show that assumptions for the $\mathrm{H_2O}$ abundance and clouds directly affect our characterization. Overall, retrievals that use physically motivated models for the $\mathrm{H_2O}$ profile and clouds perform better on the empirical Earth data. For observations of Earth with LIFE, they yield accurate estimates for the radius, pressure-temperature structure, and the abundances of $\mathrm{CO_2}$, $\mathrm{H_2O}$, and $\mathrm{O_3}$. Further, at $R=100$, a reliable and bias-free detection of the biosignature $\mathrm{CH_4}$ becomes feasible. We conclude that the community must use a diverse range of models for temperate exoplanet atmospheres to build an understanding of how different retrieval assumptions can affect the interpretation of exoplanet spectra. This will enable the characterization of distant habitable worlds and the search for life with future space-based instruments., Comment: Published in ApJ; 17 pages (main text incl. 9 figures and 5 tables) + appendix; comments are welcome

We simulate atmospheric fractionation in escaping planetary atmospheres using IsoFATE, a new open-source numerical model. We expand the parameter space studied previously to planets with tenuous atmospheres that exhibit the greatest helium and deuterium enhancement. We simulate the effects of EUV-driven photoevaporation and core-powered mass loss on deuterium-hydrogen and helium-hydrogen fractionation of sub-Neptune atmospheres around G, K, and M stars. Our simulations predict prominent populations of deuterium- and helium-enhanced planets along the upper edge of the radius valley with mean equilibrium temperatures of 370 K and as low as 150 K across stellar types. We find that fractionation is mechanism-dependent, so constraining He/H and D/H abundances in sub-Neptune atmospheres offers a unique strategy to investigate the origin of the radius valley around low-mass stars. Fractionation is also strongly dependent on retained atmospheric mass, offering a proxy for planetary surface pressure as well as a way to distinguish between desiccated enveloped terrestrials and water worlds. Deuterium-enhanced planets tend to be helium-dominated and CH4-depleted, providing a promising strategy to observe HDO in the 3.7 um window. We present a list of promising targets for observational follow-up., Comment: Accepted by ApJ

Although the scientific principles of anthropogenic climate change are well-established, existing calculations of the warming effect of carbon dioxide rely on spectral absorption databases, which obscures the physical foundations of the climate problem. Here we show how CO2 radiative forcing can be expressed via a first-principles description of the molecule's key vibrational-rotational transitions. Our analysis elucidates the dependence of carbon dioxide's effectiveness as a greenhouse gas on the Fermi resonance between the symmetric stretch mode $\nu_1$ and bending mode $\nu_2$. It is remarkable that an apparently accidental quantum resonance in an otherwise ordinary three-atom molecule has had such a large impact on our planet's climate over geologic time, and will also help determine its future warming due to human activity. In addition to providing a simple explanation of CO2 radiative forcing on Earth, our results may have implications for understanding radiation and climate on other planets., Comment: 15 pages, 7 figures, accepted to The Planetary Science Journal

In Earth's current climate, moist convective updraft speeds increase with surface warming. This trend suggests that very vigorous convection might be the norm in extremely hot and humid atmospheres, such as those undergoing a runaway greenhouse transition. However, theoretical and numerical evidence suggests that convection is actually gentle in water vapor-dominated atmospheres, implying that convective vigor may peak at some intermediate humidity level. Here, we perform small-domain convection-resolving simulations of an Earth-like atmosphere over a wide range of surface temperatures and confirm that there is indeed a peak in convective vigor, which we show occurs near T_s ~ 330 K. We show that a similar peak in convective vigor exists when the relative abundance of water vapor is changed by varying the amount of background (non-condensing) gas at fixed T_s, which may have implications for Earth's climate and atmospheric chemistry during the Hadean and Archean. We also show that Titan-like thermodynamics (i.e., a thick nitrogen atmosphere with condensing methane and low gravity) produce a peak in convective vigor at T_s ~ 95 K, which is curiously close to the current surface temperature of Titan. Plotted as functions of the saturation specific humidity at cloud base, metrics of convective vigor from both Earth-like and Titan-like experiments all peak when cloud-base air contains roughly 10% of the condensible gas by mass. Our results point to a potentially common phenomenon in terrestrial atmospheres: that moist convection is most vigorous when the condensible component is between dilute and non-dilute abundance.

Detection of water vapor in the atmosphere of temperate rocky exoplanets would be a major milestone on the path towards characterization of exoplanet habitability. Past modeling work has shown that cloud formation may prevent the detection of water vapor on Earth-like planets with surface oceans using the James Webb Space Telescope (JWST). Here we analyze the potential for atmospheric detection of H2O on a different class of targets: arid planets. Using transit spectrum simulations, we show that atmospheric H2O may be easier to be detected on arid planets with cold-trapped ice deposits on the surface, because such planets will not possess thick H2O cloud decks that limit the transit depth of spectral features. However, additional factors such as band overlap with CO2 and other gases, extinction by mineral dust, overlap of stellar and planetary H2O lines, and the ultimate noise floor obtainable by JWST still pose important challenges. For this reason, combination of space- and ground-based spectroscopic observations will be essential for reliable detection of H2O on rocky exoplanets in the future., Comment: 14 pages, 4 figures, accepted for publication in ApJL

Earth's distant past and potentially its future include extremely warm "hothouse" climate states, but little is known about how the atmosphere behaves in such states. One distinguishing characteristic of hothouse climates is that they feature lower-tropospheric radiative heating, rather than cooling, due to the closing of the water vapor infrared window regions. Previous work has suggested that this could lead to temperature inversions and significant changes in cloud cover, but no previous modeling of the hothouse regime has resolved convective-scale turbulent air motions and cloud cover directly, thus leaving many questions about hothouse radiative heating unanswered. Here, we conduct simulations that explicitly resolve convection and find that lower-tropospheric radiative heating in hothouse climates causes the hydrologic cycle to shift from a quasi-steady regime to a "relaxation oscillator" regime, in which precipitation occurs in short and intense outbursts separated by multi-day dry spells. The transition to the oscillatory regime is accompanied by strongly enhanced local precipitation fluxes, a significant increase in cloud cover, and a transiently positive (unstable) climate feedback parameter. Our results indicate that hothouse climates may feature a novel form of "temporal" convective self-organization, with implications for both cloud coverage and erosion processes.

Terrestrial-type exoplanets orbiting nearby red dwarf stars (M dwarfs) are the first potentially habitable exoplanets suitable for atmospheric characterization in the near future. Understanding the stability of water in cold-trap regions on such planets is critical because it directly impacts transmission spectroscopy observations, the global energy budget, and long-term surface water evolution. Here we diagnose the humidity distribution in idealized general circulation model (GCM) simulations of terrestrial-type exoplanets. We use the `tracer of last saturation' technique to study the saturation statistics of air parcels. We find that on synchronously rotating planets, the water vapor abundance in the nightside upper troposphere depends weakly on planetary rotation, while more water vapor builds up in the nightside lower troposphere on fast-rotating planets. We then discuss how last-saturation statistics can elucidate the multiple moist climate equilibrium states on synchronously and asynchronously rotating arid planets. We show that the multiple moist climate states arise from the cold-trapping competition between the substellar upper atmosphere and cold surface regions. We find that fast synchronously rotating planets tend to trap surface water on the nightside as a result of their weak atmospheric and strong surface cold traps compared to the slow rotating case. These results elucidate the nature of the water cycle on arid rocky exoplanets and will aid interpretation of atmospheric observations in the future., Comment: 28 pages, 10 figures, 2 tables, 1 movie, accepted for publication in The Planetary Science Journal

Reconciling the geology of Mars with models of atmospheric evolution remains a major challenge. Martian geology is characterized by past evidence for episodic surface liquid water, and geochemistry indicating a slow and intermittent transition from wetter to drier and more oxidizing surface conditions. Here we present a new model incorporating randomized injection of reducing greenhouse gases and oxidation due to hydrogen escape, to investigate the conditions responsible for these diverse observations. We find that Mars could have transitioned repeatedly from reducing (H2-rich) to oxidizing (O2-rich) atmospheric conditions in its early history. Our model predicts a generally cold early Mars, with mean annual temperatures below 240 K. If peak reducing gas release rates and background CO2 levels are high enough, it nonetheless exhibits episodic warm intervals sufficient to degrade crater walls, form valley networks and create other fluvial/lacustrine features. Our model also predicts transient buildup of atmospheric O2, which can help explain the occurrence of oxidized mineral species such as manganese oxides at Gale Crater. We suggest that the apparent Noachian--Hesperian transition from phyllosilicate deposition to sulfate deposition around 3.5 billion years ago can be explained as a combined outcome of increasing planetary oxidation, decreasing groundwater availability and a waning bolide impactor flux, which dramatically slowed the remobilization and thermochemical destruction of surface sulfates. Ultimately, rapid and repeated variations in Mars' early climate and surface chemistry would have presented both challenges and opportunities for any emergent microbial life., Comment: Published March 8 2021 in Nature Geoscience

The evolution of a single raindrop falling below a cloud is governed by fluid dynamics and thermodynamics fundamentally transferable to planetary atmospheres beyond modern Earth's. Here, we show how three properties that characterize falling raindrops -- raindrop shape, terminal velocity, and evaporation rate -- can be calculated as a function of raindrop size in any planetary atmosphere. We demonstrate that these simple, interrelated characteristics tightly bound the possible size range of raindrops in a given atmosphere, independently of poorly understood growth mechanisms. Starting from the equations governing raindrop falling and evaporation, we demonstrate that raindrop ability to vertically transport latent heat and condensible mass can be well captured by a new dimensionless number. Our results have implications for precipitation efficiency, convective storm dynamics, and rainfall rates, which are properties of interest for understanding planetary radiative balance and (in the case of terrestrial planets) rainfall-driven surface erosion., Comment: submitted to JGR: Planets; 40 pages, 8 figures, 3 tables, 4 appendices; supporting information with 9 pages, 8 figures, 1 table; associated code at https://github.com/kaitlyn-loftus/rainprops

Ultra-hot disintegrating exoplanets have been detected with tails trailing behind and/or shooting ahead of them. These tails are believed to be made of dusts that are carried upward by the supersonic flow escaping the planet's gravity field from the fiercely heated permanent day-side. Conserving angular momentum, this day-side escape flux would lead the planet in orbit. In order to explain the trailing tails in observation, radiation pressure, a repulsive force pushing the escape flow away from the host star is considered to be necessary. We here investigate whether escape could occur on the night-side as the escape flow is deflected by the pressure gradient force. We demonstrate in an idealized framework that escape flux from the night-side could dominate that from the day-side; and the former may naturally explain the commonly-observed trailing tails based on angular momentum conservation, without the need to invoke radiation pressure, which has previously been thought to be the key. We also find analytical approximations for both dayside and nightside escape fluxes, which may be applied to study planetary evolution of disintegrating planets and to infer planetary sizes from observations of the properties of their dusty tails.

We present gray gas general circulation model (GCM) simulations of the tidally locked mini-Neptune GJ 1214b. On timescales of 1,000-10,000 Earth days, our results are comparable to previous studies of the same planet, in the sense that they all exhibit two off-equatorial eastward jets. Over much longer integration times (50,000-250,000 Earth days) we find a significantly different circulation and observational features. The zonal-mean flow transitions from two off-equatorial jets to a single wide equatorial jet that has higher velocity and extends deeper. The hot spot location also shifts eastward over the integration time. Our results imply a convergence time far longer than the typical integration time used in previous studies. We demonstrate that this long convergence time is related to the long radiative timescale of the deep atmosphere and can be understood through a series of simple arguments. Our results indicate that particular attention must be paid to model convergence time in exoplanet GCM simulations, and that other results on the circulation of tidally locked exoplanets with thick atmospheres may need to be revisited., Comment: 36 pages, 17 figures, 2 tables

A key factor in determining the potential habitability of synchronously rotating planets is the strength of the atmospheric boundary layer inversion between the dark side surface and the free atmosphere. Here we analyse data obtained from polar night measurements at the South Pole and Alert Canada, which are the closest analogues on Earth to conditions on the dark sides of synchronously rotating exoplanets without and with a maritime influence, respectively. On Earth, such inversions rarely exceed 30 K in strength, because of the effect of turbulent mixing induced by phenomena such as so-called mesoscale slope winds, which have horizontal scales of 10s to 100s of km, suggesting a similar constraint to near-surface dark side inversions. We discuss the sensitivity of inversion strength to factors such as orography and the global-scale circulation, and compare them to a simulation of the planet Proxima Centauri b. Our results demonstrate the importance of comparisons with Earth data in exoplanet research, and highlight the need for further studies of the exoplanet atmospheric collapse problem using mesoscale and eddy-resolving models., Comment: Accepted by Astrophysical Journal Letters; 14 pages, 3 figures

Terrestrial-type exoplanets orbiting nearby red dwarf stars (M-dwarfs) are among the best targets for atmospheric characterization and biosignature searches in the near future. Recent evolutionary studies have suggested that terrestrial planets in the habitable zone of M-dwarfs are probably tidally locked and have limited surface water inventories as a result of their host stars' high early luminosities. Several previous climate simulations of such planets have indicated that their remaining water would be transported to the planet's permanent nightside and become trapped as surface ice, leaving the dayside devoid of water. Here we use a three-dimensional general circulation model with a water cycle and accurate radiative transfer scheme to investigate the surface water evolution on slowly rotating tidally locked terrestrial planets with limited surface water inventories. We show that there is a competition for water trapping between the nightside surface and the substellar tropopause in this type of climate system. Although under some conditions the surface water remains trapped on the nightside as an ice sheet, in other cases liquid water stabilizes in a circular area in the substellar region as a wetland. Planets with 1 bar N$_2$ and atmospheric CO$_2$ levels greater than 0.1 bar retain stable dayside liquid water, even with very small surface water inventories. Our results reveal the diversity of possible climate states on terrestrial-type exoplanets and highlight the importance of surface liquid water detection techniques for future characterization efforts., Comment: accepted for publication in ApJL

The response of the general circulation in a dry atmosphere to various atmospheric shortwave absorptivities is investigated in a three-dimensional general circulation model with grey radiation. Shortwave absorption in the atmosphere reduces the incoming radiation reaching the surface but warms the upper atmosphere, significantly shifting the habitable zone toward the star. The strong stratification under high shortwave absorptivity suppresses the Hadley cell in a manner that matches previous Hadley cell scalings. General circulation changes may be observable through cloud coverage and superrotation. The equatorial superrotation in the upper atmosphere strengthens with the shortwave opacity, as predicted based on the gradient wind of the radiative-convective equilibrium profile. There is a sudden drop of equatorial superrotation at very low shortwave opacity. This is because the Hadley cell in those cases are strong enough to fill the entire troposphere with zero momentum air from the surface. A diurnal cycle (westward motion of substellar point relative to the planet) leads to acceleration of the equatorial westerlies in general, through the enhancement of the equatorward eddy momentum transport, but the response is not completely monotonic, perhaps due to the resonance of tropical waves and the diurnal forcing.

Jupiter's atmosphere is enriched in C, N, S, P, Ar, Kr and Xe with respect to solar abundances by a factor of ~3. Gas Giant envelopes are mainly enriched through the dissolution of solids in the atmosphere, and this constant enrichment factor is puzzling since several of the above elements are not expected to have been in the solid phase in Jupiter's feeding zone; most seriously, Ar and the main carrier of N, N2, only condense at the very low temperatures, 21-26 K, associated with the outer solar nebula. We propose that a plausible solution to the enigma of Jupiter's uniform enrichment pattern is that Jupiter's core formed exterior to the N2 and Ar snowlines, beyond 30 au, resulting in a Solar composition core in all volatiles heavier than Ne. During envelope accretion and planetesimal bombardment, some of the core mixed in with the envelope causing the observed enrichment pattern. We show that this scenario naturally produces the observed atmosphere composition, even with substantial pollution from N-poor pebble and planetesimal accretion in Jupiter's final feeding zone. We note that giant core formation at large nebular radii is consistent with recent models of gas giant core formation through pebble accretion, which requires the core to form exterior to Jupiter's current location to counter rapid inward migration during the core and envelope formation process. If this scenario is common, gas giant core formation may account for many of the gaps observed in protoplanetary disks between 10s and 100 au., Comment: Accepted for publication by ApJ

Despite surface liquid water's importance to habitability, observationally diagnosing its presence or absence on exoplanets is still an open problem. Inspired within the Solar System by the differing sulfur cycles on Venus and Earth, we investigate thick sulfate (H2SO4-H2O) aerosol haze and high trace mixing ratios of SO2 gas as observable atmospheric features whose sustained existence is linked to the near absence of surface liquid water. We examine the fundamentals of the sulfur cycle on a rocky planet with an ocean and an atmosphere in which the dominant forms of sulfur are SO2 gas and H2SO4-H2O aerosols (as on Earth and Venus). We build a simple but robust model of the wet, oxidized sulfur cycle to determine the critical amounts of sulfur in the atmosphere-ocean system required for detectable levels of SO2 and a detectable haze layer. We demonstrate that for physically realistic ocean pH values (pH > 6) and conservative assumptions on volcanic outgassing, chemistry, and aerosol microphysics, surface liquid water reservoirs with greater than 1e-3 Earth oceans are incompatible with a sustained observable H2SO4-H2O haze layer and sustained observable levels of SO2. Thus, we propose the observational detection of a H2SO4-H2O haze layer and of SO2 gas as two new remote indicators that a planet does not host significant surface liquid water., Comment: published in ApJ; associated code at https://github.com/kaitlyn-loftus/alien-sulfur-cycles

Low-gravity waterworlds ($M\lesssim 0.1 M_{\oplus}$) are of interest for their potential habitability. The weakly bound atmospheres of such worlds have proportionally larger radiative surfaces and are more susceptible to escape. We conduct a unified investigation into these phenomena, combining analytical energy balance and hydrodynamic escape with line-by-line radiative transfer calculations. Because outgoing radiation is forced to increase with surface temperature by the expansion of the radiative surface, we find that these worlds do not experience a runaway greenhouse. Furthermore, we show that a long-lived liquid water habitable zone is possible for low-gravity waterworlds of sufficient mass. Its inner edge is set by the rate of atmospheric escape, because a short-lived atmosphere limits the time available for life to evolve. In describing the physics of the parameter space transition from "planet-like" to "comet-like", our model produces a lower bound for habitability in terms of gravity. These results provide valuable insights in the context of the ongoing hunt for habitable exoplanets and exomoons., Comment: published in ApJ; text updated to match published version

The radii and orbital periods of 4000+ confirmed/candidate exoplanets have been precisely measured by the Kepler mission. The radii show a bimodal distribution, with two peaks corresponding to smaller planets (likely rocky) and larger intermediate-size planets, respectively. While only the masses of the planets orbiting the brightest stars can be determined by ground-based spectroscopic observations, these observations allow calculation of their average densities placing constraints on the bulk compositions and internal structures. Yet an important question about the composition of planets ranging from 2 to 4 Earth radii still remains. They may either have a rocky core enveloped in a H2-He gaseous envelope (gas dwarfs) or contain a significant amount of multi-component, H2O-dominated ices/fluids (water worlds). Planets in the mass range of 10-15 Earth masses, if half-ice and half-rock by mass, have radii of 2.5 Earth radii, which exactly match the second peak of the exoplanet radius bimodal distribution. Any planet in the 2-4 Earth radii range requires a gas envelope of at most a few mass percentage points, regardless of the core composition. To resolve the ambiguity of internal compositions, we use a growth model and conduct Monte Carlo simulations to demonstrate that many intermediate-size planets are water worlds., Comment: PNAS link: https://www.pnas.org/content/116/20/9723 Complete data and mass-radius tables are available at: https://www.cfa.harvard.edu/~lzeng/planetmodels.html

We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets., Comment: Submitted as an Astro2020 Science White Paper

Exploring diverse planetary atmospheres requires modeling tools that are both accurate and flexible. Here, we develop a three-dimensional general circulation model (3D GCM) that for the first time uses a line-by-line approach to describe the radiative transfer. We validate our GCM by comparing with published results done by different 1D and 3D models. To demonstrate the versatility of the model, we apply the GCM to the hot Earth-sized exoplanet GJ 1132b and study its climate and circulation assuming an atmosphere dominated by abiotic oxygen (O$_2$). Our simulations show that a minor CO$_2$ composition can change the circulation pattern substantially, intensifying the equatorial superrotation in particular. Computation of the phase-resolved spectroscopy indicates that the vertical profile of the superrotating jet could be inferred in future spectro-photometric observations by the phase shift of the hotspot in the CO$_2$ principle absorption band centered at 667 cm$^{-1}$. We also show that atmospheric mass could potentially be constrained by the phase amplitude in the O$_2$ vibrational fundamental band for planets with O$_2$-rich atmospheres, although further experimental and/or theoretical O$_2$-O$_2$ collision-induced absorption data at high temperatures is needed to confirm this. More physical schemes such as moist dynamics will be implemented in the GCM in the future so that it can be used to tackle a wide variety of planetary climate problems., Comment: accepted for publication in ApJ

While recently discovered exotic new planet-types have both challenged our imaginations and broadened our knowledge of planetary system workings, perhaps the most compelling objective of exoplanet science is to detect and characterize habitable and possibly inhabited worlds orbiting in other star systems. For the foreseeable future, characterizations of extrasolar planets will be made via remote sensing of planetary spectroscopic and temporal signals, along with careful fitting of this data to advanced models of planets and their atmospheres. Terrestrial planets are small and significantly more challenging to observe compared to their larger gaseous brethren; however observatories coming on-line in the coming decade will begin to allow their characterization. Still, it is not enough to invest only in observational endeavors. Comprehensive modeling of planetary atmospheres is required in order to fully understand what it is that our grand telescopes see in the night-sky. In our quest to characterize habitable, and possibly inhabited worlds, 3D general circulation models (GCMs) should be used to evaluate potential climate states and their associated temporal and spatial dependent observable signals. 3D models allow for coupled, self-consistent, multi-dimensional simulations, which can realistically simulate the climates of terrestrial extrasolar planets. A complete theoretical understanding of terrestrial exoplanetary atmospheres, gained through comprehensive 3D modeling, is critical for interpreting spectra of exoplanets taken from current and planned instruments, and is critical for designing future missions that aim to measure spectra of potentially habitable exoplanets as one of their key science goals. We recommend continued institutional support for 3D GCM modeling teams that focus on planetary and exoplanetary applications., Comment: NAS ASTRO 2020, White Paper! 2 figures, and 10 pages

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., Comment: White paper submitted to the NAS Astro 2020 decadal survey(5 pages, 1 figure + cover page, 13 co-authors + 10 co-signers). First author name added to description. In this arxiv version, one more co-signer was added after the NAS version was submitted. Link: http://surveygizmoresponseuploads.s3.amazonaws.com/fileuploads/623127/4458621/64-2dc3fcc201a3075e3da6233f1a91c668_RamirezRamsesM.pdf

The Origins Space Telescope (OST) is one of four mission concepts currently being studied by NASA in preparation for the Astrophysics 2020 Decadal Survey. With active cooling (~4 K), OST will be sensitive in mid- to far-IR wavelengths, using imaging and spectroscopy to probe the furthest reaches of our galaxies, trace the path of water through star and planet formation, and place thermochemical constraints on the atmospheres of exoplanets ranging in size from Jupiter to Earth. This contribution to the Exoplanet Science Strategy committee discusses the significant advancements that the OST Mid-Infrared Imager, Spectrometer, and Coronagraph (MISC) instrument can make in studying cool planetary atmospheres. We particularly focus on the atmospheres of transiting rocky planets in the habitable zones of mid-to-late M stars. We discuss how OST thermal infrared observations can significantly enhance our understanding of the temperature structure and molecular abundances of biologically interesting gases on these worlds, including O3, CH4, H2O, and CO2., Comment: White paper submitted to The National Academies of Science, Engineering, and Medicine Exoplanet Science Strategy Committee

This is a white paper in response to the National Academy of Sciences "Exoplanet Science Strategy" call. We summarize recent advances in theoretical habitability studies and argue that such studies will remain important for guiding and interpreting observations. Interactions between 1-D and 3-D climate modelers will be necessary to resolve recent discrepancies in model results and improve habitability studies. Observational capabilities will also need improvement. Although basic observations can be performed with present capabilities, technological advances will be necessary to improve climate models to the level needed for planetary habitability studies., Comment: This is a white paper submitted to the National Academies 2018 "Exoplanet Science Strategy" call(6 pages total, including cover page). Corrected references section in this version. http://sites.nationalacademies.org/SSB/CurrentProjects/SSB_180659

Potentially-habitable planets orbiting M-dwarfs are of intense astrobiological interest because they are the only rocky worlds accessible to biosignature search over the next 10+ years due to a confluence of observational effects. Simultaneously, recent experimental and theoretical work suggests that UV light may have played a key role in the origin of life on Earth, and especially the origin of RNA. Characterizing the UV environment on M-dwarfs planets is important to understanding whether life as we know it could emerge on such worlds. In this work, we couple radiative transfer models to observed M-dwarf spectra to determine the UV environment on prebiotic Earth-analog planets orbiting M-dwarfs. We calculate dose rates to quantify the impact of different host stars on prebiotically-important photoprocesses. We find that M-dwarf planets have access to 100-1000 times less bioactive UV fluence than the young Earth. It is unclear whether UV-sensitive prebiotic chemistry that may have been important to abiogenesis, such as the only known prebiotically plausible pathways for pyrimidine ribonucleotide synthesis, could function on M-dwarf planets. This uncertainty affects objects like the recently-discovered habitable-zone planets orbiting Proxima Centauri, TRAPPIST-1, and LHS 1140. Laboratory studies of the sensitivity of putative prebiotic pathways to irradiation level are required to resolve this uncertainty. If steady-state M-dwarf UV output is insufficient to power these pathways, transient elevated UV irradiation due to flares may suffice; laboratory studies can constrain this possibility as well., Comment: 18 pages; 6 figures; published in ApJ

Recent findings suggest Mars may have been a clement environment for the emergence of life, and may even have compared favorably to Earth in this regard. These findings have revived interest in the hypothesis that prebiotically important molecules or even nascent life may have formed on Mars and been transferred to Earth. UV light plays a key role in prebiotic chemistry. Characterizing the early Martian surface UV environment is key to understanding how Mars compares to Earth as a venue for prebiotic chemistry. Here, we present two-stream multi-layer calculations of the UV surface radiance on Mars at 3.9 Ga, to constrain the surface UV environment as a function of atmospheric state. We explore a wide range of atmospheric pressures, temperatures and compositions, corresponding to the diversity of Martian atmospheric states consistent with available constraints. We include the effects of clouds and dust. We calculate dose rates to quantify the effect of different atmospheric states on UV-sensitive prebiotic chemistry. We find that for normative clear-sky CO2-H2O atmospheres, the UV environment on young Mars is comparable to young Earth. This similarity is robust to moderate cloud cover: thick clouds ($\tau$>100) are required to significantly affect the Martian UV environment, because cloud absorption is degenerate with atmospheric CO2. On the other hand, absorption from SO2, H2S, and dust is nondegenerate with CO2, meaning if they can build up to high levels, surface UV fluence will be suppressed. These absorbers have spectrally variable absorption, meaning that their presence affects prebiotic pathways in different ways. In particular, high SO2 environments may admit UV fluence that favors pathways conducive to abiogenesis over pathways unfavorable to it. However, better measurements of the spectral quantum yields of these pathways are required to evaluate this hypothesis definitively., Comment: In press at Astrobiology; minor errors fixed relative to v1,2

The evidence for abundant liquid water on early Mars despite the faint young Sun is a long-standing problem in planetary research. Here we present new ab initio spectroscopic and line-by-line climate calculations of the warming potential of reduced atmospheres on early Mars. We show that the strength of both CO2-H2 and CO2-CH4 collision-induced absorption (CIA) has previously been significantly underestimated. Contrary to previous expectations, methane could have acted as a powerful greenhouse gas on early Mars due to CO2-CH4 CIA in the critical 250-500 cm^-1 spectral window region. In atmospheres of 0.5 bar CO2 or more, percent levels of H2 or CH4 raise annual mean surface temperatures by tens of degrees, with temperatures reaching 273 K for pressures of 1.25-2~bar and 2-10% of H2 and CH4. Methane and hydrogen produced following aqueous alteration of Mars' crust could have combined with volcanically outgassed CO2 to form transient atmospheres of this composition 4.5-3.5 Ga. Our results also suggest that inhabited exoplanets could retain surface liquid water at significant distances from their host stars., Comment: published in GRL Jan 2017, this is final preprint version

GJ 1132 b is a nearby Earth-sized exoplanet transiting an M dwarf, and is amongst the most highly characterizable small exoplanets currently known. In this paper we study the interaction of a magma ocean with a water-rich atmosphere on GJ 1132b and determine that it must have begun with more than 5 wt% initial water in order to still retain a water-based atmosphere. We also determine the amount of O2 that can build up in the atmosphere as a result of hydrogen dissociation and loss. We find that the magma ocean absorbs at most ~10% of the O2 produced, whereas more than 90% is lost to space through hydrodynamic drag. The most common outcome for GJ 1132 b from our simulations is a tenuous atmosphere dominated by O2, although for very large initial water abundances atmospheres with several thousands of bars of O2 are possible. A substantial steam envelope would indicate either the existence of an earlier H2 envelope or low XUV flux over the system's lifetime. A steam atmosphere would also imply the continued existence of a magma ocean on GJ 1132 b. Further modeling is needed to study the evolution of CO2 or N2-rich atmospheres on GJ 1132 b., Comment: 14 pages, 11 figures, accepted at ApJ

The nature of the early Martian climate is one of the major unanswered questions of planetary science. Key challenges remain, but a new wave of orbital and in situ observations and improvements in climate modeling have led to significant advances over the last decade. Multiple lines of geologic evidence now point to an episodically warm surface during the late Noachian and early Hesperian periods 3-4 Ga. The low solar flux received by Mars in its first billion years and inefficiency of plausible greenhouse gases such as CO2 means that the steady-state early Martian climate was likely cold. A denser CO2 atmosphere would have caused adiabatic cooling of the surface and hence migration of water ice to the higher altitude equatorial and southern regions of the planet. Transient warming caused melting of snow and ice deposits and a temporarily active hydrological cycle, leading to erosion of the valley networks and other fluvial features. Precise details of the warming mechanisms remain unclear, but impacts, volcanism and orbital forcing all likely played an important role. The lack of evidence for glaciation across much of Mars' ancient terrain suggests the late Noachian surface water inventory was not sufficient to sustain a northern ocean. While mainly inhospitable on the surface, early Mars may nonetheless have presented significant opportunities for the development of microbial life., Comment: In press. See published version for enhanced figures and sidenotes (http://www.annualreviews.org/doi/abs/10.1146/annurev-earth-060115-012355)

We use a 3D general circulation model to compare the primitive Martian hydrological cycle in "warm and wet" and "cold and icy" scenarios. In the warm and wet scenario, an anomalously high solar flux or intense greenhouse warming artificially added to the climate model are required to maintain warm conditions and an ice-free northern ocean. Precipitation shows strong surface variations, with high rates around Hellas basin and west of Tharsis but low rates around Margaritifer Sinus (where the observed valley network drainage density is nonetheless high). In the cold and icy scenario, snow migration is a function of both obliquity and surface pressure, and limited episodic melting is possible through combinations of seasonal, volcanic and impact forcing. At surface pressures above those required to avoid atmospheric collapse (~0.5 bar) and moderate to high obliquity, snow is transported to the equatorial highland regions where the concentration of valley networks is highest. Snow accumulation in the Aeolis quadrangle is high, indicating an ice-free northern ocean is not required to supply water to Gale crater. At lower surface pressures and obliquities, both H2O and CO2 are trapped as ice at the poles and the equatorial regions become extremely dry. The valley network distribution is positively correlated with snow accumulation produced by the cold and icy simulation at 41.8 degrees obliquity but uncorrelated with precipitation produced by the warm and wet simulation. Because our simulations make specific predictions for precipitation patterns under different climate scenarios, they motivate future targeted geological studies., Comment: 20 pages, 14 figures, 2 tables. In press at JGR (Planets)

Atmospheric collapse is likely to be of fundamental importance to tidally locked rocky exoplanets but remains understudied. Here, general results on the heat transport and stability of tidally locked terrestrial-type atmospheres are reported. First, the problem is modeled with an idealized 3D general circulation model (GCM) with gray gas radiative transfer. It is shown that over a wide range of parameters the atmospheric boundary layer, rather than the large-scale circulation, is the key to understanding the planetary energy balance. Through a scaling analysis of the interhemispheric energy transfer, theoretical expressions for the day-night temperature difference and surface wind speed are created that reproduce the GCM results without tuning. Next, the GCM is used with correlated-k radiative transfer to study heat transport for two real gases (CO2 and CO). For CO2, empirical formulae for the collapse pressure as a function of planetary mass and stellar flux are produced, and critical pressures for atmospheric collapse at Earth's stellar flux are obtained that are around five times higher (0.14 bar) than previous gray gas estimates. These results provide constraints on atmospheric stability that will aid in future interpretation of observations and exoplanet habitability modeling., Comment: 48 pages double-spaced, 11 figures, 2 tables

Several clues indicate that Titan's atmosphere has been depleted in methane during some period of its history, possibly as recently as 0.5-1 billion years ago. It could also happen in the future. Under these conditions, the atmosphere becomes only composed of nitrogen with a range of temperature and pressure allowing liquid or solid nitrogen to condense. Here, we explore these exotic climates throughout Titan's history with a 3D Global Climate Model (GCM) including the nitrogen cycle and the radiative effect of nitrogen clouds. We show that for the last billion years, only small polar nitrogen lakes should have formed. Yet, before 1 Ga, a significant part of the atmosphere could have condensed, forming deep nitrogen polar seas, which could have flowed and flooded the equatorial regions. Alternatively, nitrogen could be frozen on the surface like on Triton, but this would require an initial surface albedo higher than 0.65 at 4 Ga. Such a state could be stable even today if nitrogen ice albedo is higher than this value. According to our model, nitrogen flows and rain may have been efficient to erode the surface. Thus, we can speculate that a paleo-nitrogen cycle may explain the erosion and the age of Titan's surface, and may have produced some of the present valley networks and shorelines. Moreover, by diffusion of liquid nitrogen in the crust, a paleo-nitrogen cycle could be responsible of the flattening of the polar regions and be at the origin of the methane outgassing on Titan., Comment: Accepted for publication in Icarus on July 7, 2014

Detection of life on other planets requires identification of biosignatures, i.e., observable planetary properties that robustly indicate the presence of a biosphere. One of the most widely accepted biosignatures for an Earth-like planet is an atmosphere where oxygen is a major constituent. Here we show that lifeless habitable zone terrestrial planets around any star type may develop oxygen-dominated atmospheres as a result of water photolysis, because the cold trap mechanism that protects H2O on Earth is ineffective when the atmospheric inventory of non-condensing gases (e.g., N2, Ar) is low. Hence the spectral features of O2 and O3 alone cannot be regarded as robust signs of extraterrestrial life., Comment: 5 pages, 3 figures, accepted for publication in ApJ Letters. Minor changes from previous version

Because the solar luminosity increases over geological timescales, Earth climate is expected to warm, increasing water evaporation which, in turn, enhances the atmospheric greenhouse effect. Above a certain critical insolation, this destabilizing greenhouse feedback can "runaway" until all the oceans are evaporated. Through increases in stratospheric humidity, warming may also cause oceans to escape to space before the runaway greenhouse occurs. The critical insolation thresholds for these processes, however, remain uncertain because they have so far been evaluated with unidimensional models that cannot account for the dynamical and cloud feedback effects that are key stabilizing features of Earth's climate. Here we use a 3D global climate model to show that the threshold for the runaway greenhouse is about 375 W/m$^2$, significantly higher than previously thought. Our model is specifically developed to quantify the climate response of Earth-like planets to increased insolation in hot and extremely moist atmospheres. In contrast with previous studies, we find that clouds have a destabilizing feedback on the long term warming. However, subsident, unsaturated regions created by the Hadley circulation have a stabilizing effect that is strong enough to defer the runaway greenhouse limit to higher insolation than inferred from 1D models. Furthermore, because of wavelength-dependent radiative effects, the stratosphere remains cold and dry enough to hamper atmospheric water escape, even at large fluxes. This has strong implications for Venus early water history and extends the size of the habitable zone around other stars., Comment: Published in Nature. Online publication date: December 12, 2013. Accepted version before journal editing and with Supplementary Information

Different solutions have been proposed to solve the "faint young Sun problem", defined by the fact that the Earth was not fully frozen during the Archean despite the fainter Sun. Most previous studies were performed with simple 1-D radiative convective models and did not account well for the clouds and ice-albedo feedback or the atmospheric and oceanic transport of energy. We apply a global climate model (GCM) to test the different solutions to the faint young Sun problem. We explore the effect of greenhouse gases (CO2 and CH4), atmospheric pressure, cloud droplet size, land distribution, and Earth's rotation rate. We show that neglecting organic haze, 100 mbar of CO2 with 2 mbar of CH4 at 3.8 Ga and 10 mbar of CO2 with 2 mbar of CH4 at 2.5 Ga allow a temperate climate (mean surface temperature between 10{\deg}C and 20{\deg}C). Such amounts of greenhouse gases remain consistent with the geological data. Removing continents produces a warming lower than +4{\deg}C. The effect of rotation rate is even more limited. Larger droplets (radii of 17 microns versus 12 microns) and a doubling of the atmospheric pressure produce a similar warming of around +7{\deg}C. In our model, ice-free water belts can be maintained up to 25{\deg}N/S with less than 1 mbar of CO2 and no methane. An interesting cloud feedback appears above cold oceans, stopping the glaciation. Such a resistance against full glaciation tends to strongly mitigate the faint young Sun problem., Comment: Journal of Geophysical Research Atmosphere. Article published online 19 SEP 2013

Water photolysis and hydrogen loss from the upper atmospheres of terrestrial planets is of fundamental importance to climate evolution but remains poorly understood in general. Here we present a range of calculations we performed to study the dependence of water loss rates from terrestrial planets on a range of atmospheric and external parameters. We show that CO2 can only cause significant water loss by increasing surface temperatures over a narrow range of conditions, with cooling of the middle and upper atmosphere acting as a bottleneck on escape in other circumstances. Around G-stars, efficient loss only occurs on planets with intermediate CO2 atmospheric partial pressures (0.1 to 1 bar) that receive a net flux close to the critical runaway greenhouse limit. Because G-star total luminosity increases with time but XUV/UV luminosity decreases, this places strong limits on water loss for planets like Earth. In contrast, for a CO2-rich early Venus, diffusion limits on water loss are only important if clouds caused strong cooling, implying that scenarios where the planet never had surface liquid water are indeed plausible. Around M-stars, water loss is primarily a function of orbital distance, with planets that absorb less flux than ~270 W/m^2 (global mean) unlikely to lose more than one Earth ocean of H2O over their lifetimes unless they lose all their atmospheric N2/CO2 early on. Because of the variability of H2O delivery during accretion, our results suggest that many 'Earth-like' exoplanets in the habitable zone may have ocean-covered surfaces, stable CO2/H2O-rich atmospheres, and high mean surface temperatures., Comment: Accepted for publication in ApJ with minor revisions

The investigation of planets around other stars began with the study of gas giants, but is now extending to the discovery and characterization of super-Earths and terrestrial planets. Motivated by this observational tide, we survey the basic dynamical principles governing the atmospheric circulation of terrestrial exoplanets, and discuss the interaction of their circulation with the hydrological cycle and global-scale climate feedbacks. Terrestrial exoplanets occupy a wide range of physical and dynamical conditions, only a small fraction of which have yet been explored in detail. Our approach is to lay out the fundamental dynamical principles governing the atmospheric circulation on terrestrial planets--broadly defined--and show how they can provide a foundation for understanding the atmospheric behavior of these worlds. We first survey basic atmospheric dynamics, including the role of geostrophy, baroclinic instabilities, and jets in the strongly rotating regime (the "extratropics") and the role of the Hadley circulation, wave adjustment of the thermal structure, and the tendency toward equatorial superrotation in the slowly rotating regime (the "tropics"). We then survey key elements of the hydrological cycle, including the factors that control precipitation, humidity, and cloudiness. Next, we summarize key mechanisms by which the circulation affects the global-mean climate, and hence planetary habitability. In particular, we discuss the runaway greenhouse, transitions to snowball states, atmospheric collapse, and the links between atmospheric circulation and CO2 weathering rates. We finish by summarizing the key questions and challenges for this emerging field in the future., Comment: Invited review, in press for the Arizona Space Science Series book "Comparative Climatology of Terrestrial Planets" (S. Mackwell, M. Bullock, and J. Harder, editors). 56 pages, 26 figures

The inner edge of the classical habitable zone is often defined by the critical flux needed to trigger the runaway greenhouse instability. This 1D notion of a critical flux, however, may not be so relevant for inhomogeneously irradiated planets, or when the water content is limited (land planets). Here, based on results from our 3D global climate model, we find that the circulation pattern can shift from super-rotation to stellar/anti stellar circulation when the equatorial Rossby deformation radius significantly exceeds the planetary radius. Using analytical and numerical arguments, we also demonstrate the presence of systematic biases between mean surface temperatures or temperature profiles predicted from either 1D or 3D simulations. Including a complete modeling of the water cycle, we further demonstrate that for land planets closer than the inner edge of the classical habitable zone, two stable climate regimes can exist. One is the classical runaway state, and the other is a collapsed state where water is captured in permanent cold traps. We identify this "moist" bistability as the result of a competition between the greenhouse effect of water vapor and its condensation. We also present synthetic spectra showing the observable signature of these two states. Taking the example of two prototype planets in this regime, namely Gl581c and HD85512b, we argue that they could accumulate a significant amount of water ice at their surface. If such a thick ice cap is present, gravity driven ice flows and geothermal flux should come into play to produce long-lived liquid water at the edge and/or bottom of the ice cap. Consequently, the habitability of planets at smaller orbital distance than the inner edge of the classical habitable zone cannot be ruled out. Transiting planets in this regime represent promising targets for upcoming observatories like EChO and JWST., Comment: Accepted for publication in Astronomy and Astrophysics, complete abstract in the pdf, 18 pages, 18 figures

On the basis of geological evidence, it is often stated that the early martian climate was warm enough for liquid water to flow on the surface thanks to the greenhouse effect of a thick atmosphere. We present 3D global climate simulations of the early martian climate performed assuming a faint young sun and a CO2 atmosphere with pressure between 0.1 and 7 bars. The model includes a detailed radiative transfer model using revised CO2 gas collision induced absorption properties, and a parameterisation of the CO2 ice cloud microphysical and radiative properties. A wide range of possible climates is explored by using various values of obliquities, orbital parameters, cloud microphysic parameters, atmospheric dust loading, and surface properties. Unlike on present day Mars, for pressures higher than a fraction of a bar, surface temperatures vary with altitude because of the adiabatic cooling and warming of the atmosphere when it moves vertically. In most simulations, CO2 ice clouds cover a major part of the planet but greenhouse effect does not exceed +15 K. We find that a CO2 atmosphere could not have raised the annual mean temperature above 0{\deg}C anywhere on the planet. The collapse of the atmosphere into permanent CO2 ice caps is predicted for pressures higher than 3 bar, or conversely at pressure lower than one bar if the obliquity is low enough. Summertime diurnal mean surface temperatures above 0{\deg}C (a condition which could have allowed rivers to form) are predicted for obliquity larger than 40{\deg} at high latitudes but not in locations where most valley networks are observed. In the absence of other warming mechanisms, our climate model results are thus consistent with a cold early Mars scenario in which non climatic mechanisms must occur to explain the evidence for liquid water. In a companion paper by Wordsworth et al., we simulate the hydrological cycle on such a planet., Comment: Icarus, In press. 23 pages, 18 figures

It has been suggested that the recently discovered exoplanet GJ581d might be able to support liquid water due to its relatively low mass and orbital distance. However, GJ581d receives 35% less stellar energy than Mars and is probably locked in tidal resonance, with extremely low insolation at the poles and possibly a permanent night side. Under such conditions, it is unknown whether any habitable climate on the planet would be able to withstand global glaciation and / or atmospheric collapse. Here we present three-dimensional climate simulations that demonstrate GJ581d will have a stable atmosphere and surface liquid water for a wide range of plausible cases, making it the first confirmed super-Earth (exoplanet of 2-10 Earth masses) in the habitable zone. We find that atmospheres with over 10 bar CO2 and varying amounts of background gas (e.g., N2) yield global mean temperatures above 0 degrees Celsius for both land and ocean-covered surfaces. Based on the emitted IR radiation calculated by the model, we propose observational tests that will allow these cases to be distinguished from other possible scenarios in the future., Comment: Accepted to the Astrophysical Journal Letters; 9 pages, 1 table, 4 figures