Terrestrial planets orbiting within the habitable zones of M-stars are likely to become tidally locked in a 1:1 spin:orbit configuration and are prime targets for future characterization efforts. An issue of importance for the potential habitability of terrestrial planets is whether they could experience snowball events (periods of global glaciation). Previous work using an intermediate complexity atmospheric Global Climate Model (GCM) with no ocean heat transport suggested that tidally locked planets would smoothly transition to a snowball, in contrast with Earth, which has bifurcations and hysteresis in climate state associated with global glaciation. In this paper, we use a coupled ocean-atmosphere GCM (ROCKE-3D) to model tidally locked planets with no continents. We chose this configuration in order to consider a case that we expect to have high ocean heat transport. We show that including ocean heat transport does not reintroduce the snowball bifurcation. An implication of this result is that a tidally locked planet in the habitable zone is unlikely to be found in a snowball state for a geologically significant period of time.

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.

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

The ice-albedo feedback on rapidly-rotating terrestrial planets in the habitable zone can lead to abrupt transitions (bifurcations) between a warm and a snowball (ice-covered) state, bistability between these states, and hysteresis in planetary climate. This is important for planetary habitability because snowball events may trigger rises in the complexity of life, but could also endanger complex life that already exists. Recent work has shown that planets tidally locked in synchronous rotation states will transition smoothly into the snowball state rather than experiencing bifurcations. Here we investigate the structure of snowball bifurcations on planets that are tidally influenced, but not synchronously rotating, so that they experience long solar days. We use PlaSIM, an intermediate-complexity global climate model, with a thermodynamic mixed layer ocean and the Sun's spectrum. We find that the amount of hysteresis (range in stellar flux for which there is bistability in climate) is significantly reduced for solar days with lengths of tens of Earth days, and disappears for solar days of hundreds of Earth days. These results suggest that tidally influenced planets orbiting M and K-stars that are not synchronously rotating could have much less hysteresis associated with the snowball bifurcations than they would if they were rapidly rotating. This implies that the amount of time it takes them to escape a snowball state via CO$_2$ outgassing would be greatly reduced, as would the period of cycling between the warm and snowball state if they have a low CO$_2$ outgassing rate.

The TRAPPIST-1, Proxima Centauri, and LHS 1140 systems are the most exciting prospects for future follow-up observations of potentially inhabited planets. All orbit nearby M-stars and are likely tidally locked in 1:1 spin-orbit states, which motivates the consideration of the effects that tidal locking might have on planetary habitability. On Earth, periods of global glaciation (snowballs) may have been essential for habitability and remote signs of life (biosignatures) because they are correlated with increases in the complexity of life and in the atmospheric oxygen concentration. In this paper we investigate the snowball bifurcation (sudden onset of global glaciation) on tidally locked planets using both an energy balance model and an intermediate-complexity global climate model. We show that tidally locked planets are unlikely to exhibit a snowball bifurcation as a direct result of the spatial pattern of insolation they receive. Instead they will smoothly transition from partial to complete ice coverage and back. A major implication of this work is that tidally locked planets with an active carbon cycle should not be found in a snowball state. Moreover, this work implies that tidally locked planets near the outer edge of the habitable zone with low CO2 outgassing fluxes will equilibrate with a small unglaciated substellar region rather than cycling between warm and snowball states. More work is needed to determine how the lack of a snowball bifurcation might affect the development of life on a tidally locked planet.

We use Hubble Space Telescope Wide-Field Camera 3 (HST/WFC3) rest-frame optical imaging to select a pilot sample of star-forming galaxies in the redshift range z = 2.00-2.65 whose multi-component morphologies are consistent with expectations for major mergers. We follow up this sample of major merger candidates with Keck/NIRSPEC longslit spectroscopy obtained in excellent seeing conditions (FWHM ~ 0.5 arcsec) to obtain Halpha-based redshifts of each of the morphological components in order to distinguish spectroscopic pairs from false pairs created by projection along the line of sight. Of six pair candidates observed, companions (estimated mass ratios 5:1 and 7:1) are detected for two galaxies down to a 3sigma limiting emission-line flux of ~ 10^{-17} erg/s/cm2. This detection rate is consistent with a ~ 50% false pair fraction at such angular separations (1-2 arcsec), and with recent claims that the star-formation rate (SFR) can differ by an order of magnitude between the components in such mergers. The two spectroscopic pairs identified have total SFR, SFR surface densities, and stellar masses consistent on average with the overall z ~ 2 star forming galaxy population.