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 Sunlike 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.

Planetary climate can be affected by the interaction of the host star spectral energy distribution with the wavelength-dependent reflectivity of ice and snow. In this study, we explored this effect with 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 (3-D) general circulation model was 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 simulated 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 exhibited a stronger ice-albedo feedback. We found that ice extent was much greater on a planet orbiting an F-dwarf star than on a planet orbiting a G-dwarf star at an equivalent flux distance, and that ice-covered conditions occurred on an F-dwarf planet with only a 2% reduction in instellation relative to the present instellation on Earth, assuming fixed CO 2 (present atmospheric level on Earth). A similar planet orbiting the Sun at an equivalent flux distance required an 8% reduction in instellation, while a planet orbiting an Mdwarf star required an additional 19% reduction in instellation to become ice-covered, equivalent to 73% of the modern solar constant. The reduction in instellation must be larger for planets orbiting cooler stars due in large part to the stronger absorption of longer-wavelength radiation by icy surfaces on these planets in addition to stronger absorption by water vapor and CO 2 in their atmospheres, which provides increased downwelling longwave radiation. Lowering the IR and visible-band surface ice and snow albedos for an M-dwarf planet increased the planet's climate stability against changes in instellation and slowed the descent into global ice coverage. The surface ice-albedo feedback effect becomes less important at the outer edge of the habitable zone, where atmospheric CO 2 could be expected to be high such that it maintains clement conditions for surface liquid water. We showed that *3-10 bar of CO 2 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 CO 2 is needed to maintain open water for a planet orbiting an M-dwarf star than would be the case for hotter main-sequence stars.

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 Mdwarf 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 CO 2 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 CO 2 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 CO 2 . 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.

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 Without long timescale regulation of non-condensable greenhouse species at Earth-like temperatures and pressures, such as CO 2 , 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.

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.

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 (Borucki et al., 2013), 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 ≤ e ≤ 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 the 3-D Laboratoire de Météorologie Dynamique (LMD) Generic global climate model (GCM) indicate that the surface habitability of this planet is sensitive to orbital configuration. With 3 bar of CO2 in its atmosphere, we find that 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° and 90°). A climate similar to that of modern-day Earth is possible for the entire range of stable eccentricities if atmospheric CO2 is increased to 5 bar levels. In a low-CO2 case (Earth-like levels), simulations with version 4 of the Community Climate System Model (CCSM4) GCM 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, CO2 levels above 3 bar would be required to distribute enough heat to the nightside 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. Key Words: Extrasolar planets—Habitability—Planetary environments. Astrobiology 16, 443–464.

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. 2 Rushby et al.planet, which describes the total proportion or fraction of incident stellar flux, across all wavelengths, that is reflected back to space. This is in contrast to the geometric, spherical, or V-band albedo, which is wavelength and phase angle dependent. The relationship between albedo, temperature, and global ice cover represents a positive feedback within the planetary system, known as the ice-albedo feedback, that operates as ice cover advances equator-ward from the poles (due to a reduction in temperature) and planetary broadband albedo increases, thereby compounding the reduction in temperature and the further growth of icy surfaces of higher albedo (Lian & Cess, 1977;von Paris et al., 2013). Frozen surface features comprised of the same overlying material on planets around stars of different spectral types have considerably different Bond albedos, due to their spectral energy distributions (SED), which differ significantly. For example, smaller, redder stars, have peak output in the ∼0.8 to 1.2 µm range, where water ice and snow are particularly absorptive (Shields et al., 2013;Joshi & Haberle, 2012;von Paris et al., 2013). The geometric albedo of ice and snow begins to decrease at wavelengths greater than ∼0.5 µm (see figure 1), and therefore the albedo of snow and ice covered surfaces on planets orb...