A planet's climate can be strongly affected by its orbital eccentricity and obliquity. Here we use a 1-dimensional energy balance model modified to include a simple runaway greenhouse (RGH) parameterization to explore the effects of these two parameters on the climate of Earth-like aqua planets - completely ocean-covered planets - orbiting F-, G-, K-, and M-dwarf stars. We find that the range of instellations for which planets exhibit habitable surface conditions throughout an orbit decreases with increasing eccentricity. However, the appearance of temporarily habitable conditions during an orbit creates an eccentric habitable zone (EHZ) that is sensitive to orbital eccentricity and obliquity, planetary latitude, and host star spectral type. We find that the fraction of a planet's orbit over which it exhibits habitable surface conditions is larger on eccentric planets orbiting M-dwarf stars, due to the lower broadband planetary albedos of these planets. Planets with larger obliquities have smaller EHZs, but exhibit warmer climates if they do not enter a snowball state during their orbits. We also find no transient runaway greenhouse state on planets at all eccentricities. Rather, planets spend their entire orbits either in a RGH or not. For G-dwarf planets receiving 100% of the modern solar constant and with eccentricities above 0.55, an entire Earth ocean inventory can be lost in 3.6 Gyr. M-dwarf planets, due to their larger incident XUV flux, can become desiccated in only 690 Myr with eccentricities above 0.38. This work has important implications for eccentric planets that may exhibit surface habitability despite technically departing from the traditional habitable zone as they orbit their host stars.

The energy balance and climate of planets can be affected by the reflective properties of their land, ocean, and frozen surfaces. Here we investigate the effect of host star spectral energy distribution (SED) on the albedo of these surfaces using a one-dimensional (1-D) energy balance model (EBM). Incorporating spectra of M-, K-, G- and F-dwarf stars, we determined the effect of varying fractional and latitudinal distribution of land and ocean surfaces as a function of host star SED on the overall planetary albedo, climate, and ice-albedo feedback response. While noting that the spatial distribution of land masses on a given planet will have an effect on the overall planetary energy balance, we find that terrestrial planets with higher average land/ocean fractions are relatively cooler and have higher albedo regardless of star type. For Earth-like planets orbiting M-dwarf stars the increased absorption of water ice in the near-infrared (NIR), where M-dwarf stars emit much of their energy, resulted in warmer global mean surface temperatures, ice lines at higher latitudes, and increased climate stability as the ice-albedo feedback became negative at high land fractions. Conversely, planets covered largely by ocean, and especially those orbiting bright stars, had a considerably different energy balance due to the contrast between the reflective land and the absorptive ocean surface, which in turn resulted in warmer average surface temperatures than land-covered planets and a stronger potential ice-albedo feedback. While dependent on the properties of individual planetary systems, our results place so constraints on a range of climate states of terrestrial exoplanets based on albedo and incident flux.

The pathways through which incoming energy is distributed between the surface and atmosphere has been analyzed for the Earth. However, the effect of the spectral energy distribution of a host star on the energy budget of an orbiting planet may be significant given the wavelength-dependent absorption properties of atmospheric CO2 and water vapor, and surface ice and snow. We have quantified the flow of energy on aqua planets orbiting M-, G-, and F-dwarf stars, using a 3D Global Climate Model with a static ocean. The atmosphere and surface of an M-dwarf planet receiving an instellation equal to 88% of the modern solar constant at the top of the atmosphere absorb 12% more incoming stellar radiation than those of a G-dwarf planet receiving 100% of the modern solar constant, and 17% more radiation than a F-dwarf planet receiving 108% of the modern solar constant, resulting in climates similar to modern-day Earth on all three planets, assuming a 24-hr rotation period and fixed CO2. At 100% instellation, a synchronously-rotating M-dwarf planet exhibits smaller flux absorption in the atmosphere and on the surface of the dayside, and a dayside mean surface temperature that is 37 K colder than its rapidly-rotating counterpart. Energy budget diagrams are included to illustrate the variations in global energy budgets as a function of host star spectral class, and can contribute to habitability assessments of planets as they are discovered.

The discovery of planets orbiting stars other than the Sun has accelerated over the past decade, and this trend will continue as new space- and ground-based observatories employ next-generation instrumentation to search the skies for habitable worlds. However, many factors and processes can affect planetary habitability and must be understood to accurately determine a planet's habitability potential. While climate models have long been used to understand and predict climate and weather patterns on the Earth, a growing community of researchers has begun to apply these models to extrasolar planets. This work has provided a better understanding of how orbital, surface, and atmospheric properties affect planetary climate and habitability; how these climatic effects might change for different stellar and planetary environments; and how the habitability and observational signatures of newly discovered planets might be influenced by these climatic factors. This review summarizes the origins and evolution of the burgeoning field of exoplanet climatology, discusses recent work using a hierarchy of computer models to identify those planets most capable of supporting life, and offers a glimpse into future directions of this quickly evolving subfield of exoplanet science.

Over the past three decades instruments on the ground and in space have discovered thousands of planets outside the solar system. These observations have given rise to an astonishingly detailed picture of the demographics of short-period planets, but are incomplete at longer periods where both the sensitivity of transit surveys and radial velocity signals plummet. Even more glaring is that the spectra of planets discovered with these indirect methods are either inaccessible (radial velocity detections) or only available for a small subclass of transiting planets with thick, clear atmospheres. Direct detection can be used to discover and characterize the atmospheres of planets at intermediate and wide separations, including non-transiting exoplanets. Today, a small number of exoplanets have been directly imaged, but they represent only a rare class of young, self-luminous super-Jovian-mass objects orbiting tens to hundreds of AU from their host stars. Atmospheric characterization of planets in the <5 AU regime, where radial velocity (RV) surveys have revealed an abundance of other worlds, is technically feasible with 30-m class apertures in combination with an advanced AO system, coronagraph, and suite of spectrometers and imagers. There is a vast range of unexplored science accessible through astrometry, photometry, and spectroscopy of rocky planets, ice giants, and gas giants. In this whitepaper we will focus on one of the most ambitious science goals --- detecting for the first time habitable-zone rocky (<1.6 R_Earth) exoplanets in reflected light around nearby M-dwarfs

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.

A possible surface type that may form in the environments of M-dwarf planets is sodium chloride dihydrate, or "hydrohalite" (NaCl $\cdot$ 2H$_2$O), which can precipitate in bare sea ice at low temperatures. Unlike salt-free water ice, hydrohalite is highly reflective in the near-infrared, where M-dwarf stars emit strongly, making the effect of the interaction between hydrohalite and the M-dwarf SED necessary to quantify. We carried out the first exploration of the climatic effect of hydrohalite-induced salt-albedo feedback on extrasolar planets, using a three-dimensional global climate model. Under fixed CO$_2$ conditions, rapidly-rotating habitable-zone M-dwarf planets receiving 65% or less of the modern solar constant from their host stars exhibit cooler temperatures when an albedo parameterization for hydrohalite is included in climate simulations, compared to simulations without such a parameterization. Differences in global mean surface temperature with and without this parameterization increase as the instellation is lowered, which may increase CO$_2$ build-up requirements for habitable conditions on planets with active carbon cycles. Synchronously-rotating habitable-zone M-dwarf planets appear susceptible to salt-albedo feedback at higher levels of instellation (90% or less of the modern solar constant) than planets with Earth-like rotation periods, due to their cooler minimum day-side temperatures. These instellation levels where hydrohalite seems most relevant correspond to several recently-discovered potentially habitable M-dwarf planets, including Proxima Centauri b, TRAPPIST-1e, and LHS 1140b, making an albedo parameterization for hydrohalite of immediate importance in future climate simulations.

The prospects for the habitability of M-dwarf planets have long been debated, due to key differences between the unique stellar and planetary environments around these low-mass stars, as compared to hotter, more luminous Sun-like stars. Over the past decade, significant progress has been made by both space- and ground-based observatories to measure the likelihood of small planets to orbit in the habitable zones of M-dwarf stars. We now know that most M dwarfs are hosts to closely-packed planetary systems characterized by a paucity of Jupiter-mass planets and the presence of multiple rocky planets, with roughly a third of these rocky M-dwarf planets orbiting within the habitable zone, where they have the potential to support liquid water on their surfaces. Theoretical studies have also quantified the effect on climate and habitability of the interaction between the spectral energy distribution of M-dwarf stars and the atmospheres and surfaces of their planets. These and other recent results fill in knowledge gaps that existed at the time of the previous overview papers published nearly a decade ago by Tarter et al. (2007) and Scalo et al. (2007). In this review we provide a comprehensive picture of the current knowledge of M-dwarf planet occurrence and habitability based on work done in this area over the past decade, and summarize future directions planned in this quickly evolving field.

Conventional definitions of habitability require abundant liquid surface water to exist continuously over geologic timescales. Water in each of its thermodynamic phases interacts with solar and thermal radiation and is the cause for strong climatic feedbacks. Thus, assessments of the habitable zone require models to include a complete treatment of the hydrological cycle over geologic time. Here, we use the Community Atmosphere Model from the National Center for Atmospheric Research to study the evolution of climate for an Earth-like planet at constant CO2, under a wide range of stellar fluxes from F-, G-, and K-dwarf main sequence stars. Around each star we find four stable climate states defined by mutually exclusive global mean surface temperatures (Ts); snowball (Ts < 235 K), waterbelt (235 K < Ts < 250 K), temperate (275 K < Ts < 315 K), and moist greenhouse (Ts > 330 K). Each is separated by abrupt climatic transitions. Waterbelt, temperate, and cooler moist greenhouse climates can maintain open-ocean against both sea-ice albedo and hydrogen escape processes respectively, and thus constitute habitable worlds. We consider the warmest possible habitable planet as having Ts ~ 355 K, at which point diffusion limited water-loss could remove an Earth ocean in ~1 Gyr. Without long timescale regulation of non-condensable greenhouse species at Earth-like temperatures and pressures, such as CO2, habitability can be maintained for an upper limit of ~2.2, ~2.4 and ~4.7 Gyr around F-, G- and K-dwarf stars respectively due to main sequence brightening.

As lower-mass stars often host multiple rocky planets, gravitational interactions among planets can have significant effects on climate and habitability over long timescales. Here we explore a specific case, Kepler-62f, a potentially habitable planet in a five-planet system with a K2V host star. N-body integrations reveal the stable range of initial eccentricities for Kepler-62f is $0.00 \leqslant e \leqslant 0.32$, absent the effect of additional, undetected planets. We simulate the tidal evolution of Kepler-62f in this range and find that, for certain assumptions, the planet can be locked in a synchronous rotation state. Simulations using LMD Generic GCM indicate that with 3 bars of CO$_2$ in its atmosphere, Kepler-62f would only be warm enough for surface liquid water at the upper limit of this eccentricity range, providing it has a high planetary obliquity (between 60$^\circ$ and 90$^\circ$). A climate similar to modern-day Earth is possible for the entire range of stable eccentricities if atmospheric CO$_2$ is increased to 5-bar levels. In a low-CO$_2$ case, simulations with CCSM4 and LMD Generic GCM indicate that increases in planetary obliquity and orbital eccentricity coupled with an orbital configuration that places the summer solstice at or near pericenter permit regions of the planet with above-freezing surface temperatures. This may melt ice sheets formed during colder seasons. If Kepler-62f is synchronously rotating and has an ocean, CO$_2$ levels above 3 bars would be required to distribute enough heat to the night side of the planet to avoid atmospheric freeze-out and permit a large enough region of open water at the planet's substellar point to remain stable. Overall, we find multiple plausible combinations of orbital and atmospheric properties that permit surface liquid water on Kepler-62f.

Distant planets in globally ice-covered, "snowball", states may depend on increases in their host stars' luminosity to become hospitable for surface life. Using a General Circulation Model (GCM), we simulated the equilibrium climate response of a planet to a range of instellations from an F-, G-, or M- dwarf star. The range of instellation that permits both complete ice cover and at least partially ice-free climate states is a measure of the climate hysteresis that a planet can exhibit. An ice-covered planet with high climate hysteresis would show a higher resistance to the initial loss of surface ice coverage with increases in instellation, and abrupt, extreme ice loss once deglaciation begins. Our simulations indicate that the climate hysteresis depends sensitively on the host star spectral energy distribution. Under fixed CO2 conditions, a planet orbiting an M-dwarf star exhibits a smaller climate hysteresis, requiring a smaller instellation to initiate deglaciation than planets orbiting hotter, brighter stars. This is due to the higher absorption of near-IR radiation by ice on the surfaces and greenhouse gases and clouds in the atmosphere of an M-dwarf planet. Increases in atmospheric CO2 further lower the climate hysteresis, as M-dwarf snowball planets exhibit a larger radiative response than G-dwarf snowball planets for the same increase in CO2. For a smaller hysteresis, planets near the outer edge of the habitable zone will thaw earlier in their evolutionary history, and will experience a less abrupt transition out of global ice cover.

Planetary climate can be affected by the interaction of the host star spectral energy distribution with the wavelength-dependent reflectivity of ice and snow. Here we explore this effect using a one dimensional (1-D), line-by-line, radiative-transfer model to calculate broadband planetary albedos as input to a seasonally varying, 1-D energy-balance climate model. A three-dimensional general circulation model is also used to explore the atmosphere's response to changes in incoming stellar radiation, or instellation, and surface albedo. Using this hierarchy of models we simulate planets covered by ocean, land, and water ice of varying grain size, with incident radiation from stars of different spectral types. Terrestrial planets orbiting stars with higher near-UV radiation exhibit a stronger ice-albedo feedback. We find that ice-covered conditions occur on an F-dwarf planet with only a 2% reduction in instellation relative to the present instellation on Earth, assuming fixed CO2 (present atmospheric level on Earth). A similar planet orbiting the Sun at an equivalent flux distance requires an 8% reduction in instellation, while a planet orbiting an M-dwarf star requires an additional 19% reduction in instellation to become ice-covered, equivalent to 73% of the modern solar constant. The surface ice-albedo feedback effect becomes less important at the outer edge of the habitable zone, where atmospheric CO2 can be expected to be high in order to maintain clement conditions for surface liquid water. We show that 3-10 bars of CO2 will entirely mask the climatic effect of ice and snow, leaving the outer limits of the habitable zone unaffected by the spectral dependence of water ice and snow albedo. However, less CO2 is needed to maintain open water for a planet orbiting an M-dwarf star, than would be the case for hotter main-sequence stars.

We have obtained the first time-resolved, disc-integrated observations of Earth's poles with the Deep Impact spacecraft as part of the EPOXI Mission of Opportunity. These data mimic what we will see when we point next-generation space telescopes at nearby exoplanets. We use principal component analysis (PCA) and rotational lightcurve inversion to characterize color inhomogeneities and map their spatial distribution from these unusual vantage points, as a complement to the equatorial views presented in Cowan et al. (2009). We also perform the same PCA on a suite of simulated rotational multi-band lightcurves from NASA's Virtual Planetary Laboratory 3D spectral Earth model. This numerical experiment allows us to understand what sorts of surface features PCA can robustly identify. We find that the EPOXI polar observations have similar broadband colors as the equatorial Earth, but with 20-30% greater apparent albedo. This is because the polar observations are most sensitive to mid-latitudes, which tend to be more cloudy than the equatorial latitudes emphasized by the original EPOXI Earth observations. The cloudiness of the mid-latitudes also manifests itself in the form of increased variability at short wavelengths in the polar observations, and as a dominant gray eigencolor in the south polar observation. We construct a simple reflectance model for a snowball Earth. By construction, our model has a higher Bond albedo than the modern Earth; its surface albedo is so high that Rayleigh scattering does not noticeably affect its spectrum. The rotational color variations occur at short wavelengths due to the large contrast between glacier ice and bare land in those wavebands. Thus we find that both the broadband colors and diurnal color variations of such a planet would be easily distinguishable from the modern-day Earth, regardless of viewing angle.