Robust atmospheric and radiative transfer modeling will be required to properly interpret reflected light and thermal emission spectra of terrestrial exoplanets. This will help break observational degeneracies between the numerous atmospheric, planetary, and stellar factors that drive planetary climate. Here we simulate the climates of Earth-like worlds around the Sun with increasingly slow rotation periods, from Earth-like to fully Sun-synchronous, using the ROCKE-3D general circulation model. We then provide these results as input to the Spectral Planet Model (SPM), which employs the SMART radiative transfer model to simulate the spectra of a planet as it would be observed from a future space-based telescope. We find that the primary observable effects of slowing planetary rotation rate are the altered cloud distributions, altitudes, and opacities which subsequently drive many changes to the spectra by altering the absorption band depths of biologically-relevant gas species (e.g., H2O, O2, and O3). We also identify a potentially diagnostic feature of synchronously rotating worlds in mid-infrared H2O absorption/emission lines.

Planetary rotation rate has a significant effect on atmospheric circulation, where the strength of the Coriolis effect in part determines the efficiency of latitudinal heat transport, altering cloud distributions, surface temperatures, and precipitation patterns. In this study we use the ROCKE-3D dynamic-ocean general circulation model to study the effects of slow rotations and increased insolations on the 'fractional habitability' and silicate weathering rate of an Earth-like world. Defining the fractional habitability f_h to be the percentage of a planet's surface which falls in the 0 <= T <= 100 C temperature regime, we find a moderate increase in f_h with a 10% and 20% increase in insolation and a possible maximum in f_h at sidereal day-lengths between 8 and 32 times that of the modern Earth. By tracking precipitation and run-off we further determine that there is a rotational regime centered on a 4-day period in which the silicate weathering rate is maximized and is particularly strongly peaked at higher overall insolations. Because of weathering's integral role in the long-term carbonate-silicate cycle, we suggest that climate stability may be strongly affected by the anticipated rotational evolution of temperate terrestrial-type worlds, and should be considered a major factor in their study. In light of our results we argue that planetary rotation period is an important factor to consider when determining the habitability of terrestrial worlds.

One popular view of Venus' climate history describes a world that has spent much of its life with surface liquid water, plate tectonics, and a stable temperate climate. Part of the basis for this optimistic scenario is the high deuterium to hydrogen ratio from the Pioneer Venus mission that was interpreted to imply Venus had a shallow ocean's worth of water throughout much of its history. Another view is that Venus had a long lived (approximately 100 million year) primordial magma ocean with a CO2 and steam atmosphere. Venus' long lived steam atmosphere would sufficient time to dissociate most of the water vapor, allow significant hydrogen escape and oxidize the magma ocean. A third scenario is that Venus had surface water and habitable conditions early in its history for a short period of time (< 1 Gyr), but that a moist/runaway greenhouse took effect because of a gradually warming sun, leaving the planet desiccated ever since. Using a general circulation model we demonstrate the viability of the first scenario using the few observational constraints available. We further speculate that Large Igneous Provinces and the global resurfacing 100s of millions of years ago played key roles in ending the clement period in its history and presenting the Venus we see today. The results have implications for what astronomers term "the habitable zone," and if Venus-like exoplanets exist with clement conditions akin to modern Earth we propose to place them in what we term the "optimistic Venus zone."

Upcoming telescopes such as the James Webb Space Telescope (JWST), or the Extremely Large Telescope (ELTs), may soon be able to characterize, through transmission, emission or reflection spectroscopy, the atmospheres of rocky exoplanets orbiting nearby M dwarfs. One of the most promising candidates is the late M dwarf system TRAPPIST-1 which has seven known transiting planets for which Transit Timing Variation (TTV) measurements suggest that they are terrestrial in nature, with a possible enrichment in volatiles. Among these seven planets, TRAPPIST-1e seems to be the most promising candidate to have habitable surface conditions, receiving ~66 % of the Earth's incident radiation, and thus needing only modest greenhouse gas inventories to raise surface temperatures to allow surface liquid water to exist. TRAPPIST-1e is therefore one of the prime targets for JWST atmospheric characterization. In this context, the modeling of its potential atmosphere is an essential step prior to observation. Global Climate Models (GCMs) offer the most detailed way to simulate planetary atmospheres. However, intrinsic differences exist between GCMs which can lead to different climate prediction and thus observability of gas and/or cloud features in transmission and thermal emission spectra. Such differences should preferably be known prior to observations. In this paper we present a protocol to inter-compare planetary GCMs. Four testing cases are considered for TRAPPIST-1e but the methodology is applicable to other rocky exoplanets in the Habitable Zone. The four test cases included two land planets composed with a modern Earth and pure CO2 atmospheres, respectively, and two aqua planets with the same atmospheric compositions. Currently, there are four participating models (LMDG, ROCKE-3D, ExoCAM, UM), however this protocol is intended to let other teams participate as well.

The habitable fraction of a planet's surface is important for the detectability of surface biosignatures. The extent and distribution of habitable areas is influenced by external parameters that control the planet's climate, atmospheric circulation, and hydrological cycle. We explore these issues using the ROCKE-3D General Circulation Model, focusing on terrestrial water fluxes and thus the potential for the existence of complex life on land. Habitability is examined as a function of insolation and planet rotation for an Earth-like world with zero obliquity and eccentricity orbiting the Sun. We assess fractional habitability using an aridity index that measures the net supply of water to the land. Earth-like planets become ``superhabitable'' (a larger habitable surface area than Earth) as insolation and day-length increase because their climates become more equable, reminiscent of past warm periods on Earth when complex life was abundant and widespread. The most slowly rotating, most highly irradiated planets, though, occupy a hydrological regime unlike any on Earth, with extremely warm, humid conditions at high latitudes but little rain and subsurface water storage. Clouds increasingly obscure the surface as insolation increases, but visibility improves for modest increases in rotation period. Thus, moderately slowly rotating rocky planets with insolation near or somewhat greater than modern Earth's appear to be promising targets for surface characterization by a future direct imaging mission.

The current goals of the astrobiology community are focused on developing a framework for the detection of biosignatures, or evidence thereof, on objects inside and outside of our solar system. A fundamental aspect of understanding the limits of habitable environments (surface liquid water) and detectable signatures thereof is the study of where the boundaries of such environments can occur. Such studies provide the basis for understanding how a once inhabitable planet might come to be uninhabitable. The archetype of such a planet is arguably Earth's sibling planet, Venus. Given the need to define the conditions that can rule out biorelated signatures of exoplanets, Venus provides a unique opportunity to explore the processes that led to a completely uninhabitable environment by our current definition of the term. Here we review the current state of knowledge regarding Venus, particularly in the context of remote sensing techniques that are being or will be employed in the search for and characterization of exoplanets. We discuss candidate Venus analogs identified by the Kepler and TESS exoplanet missions and provide an update to exoplanet demographics that can be placed in the potential runaway greenhouse regime where Venus analogs are thought to reside. We list several major outstanding questions regarding the Venus environment and the relevance of those questions to understanding the atmospheres and interior structure of exoplanets. Finally, we outline the path towards a deeper analysis of our sibling planet and the synergy to exoplanetary science.

High obliquity planets represent potentially extreme limits of terrestrial climate, as they exhibit large seasonality, a reversed annual-mean pole-to-equator gradient of stellar heating, and novel cryospheres. A suite of 3-D global climate model simulations with a dynamic ocean is performed with Earthlike atmospheres for low and high obliquity planets with various stellar fluxes, CO2 concentrations, and initial conditions to explore the propensity for high obliquity climates approaching the outer edge of the Habitable Zone to undergo global glaciation. We also simulate planets with thick CO2 or H2 atmospheres, such as those expected to develop near or beyond the outer edge of the Habitable Zone. We show that high obliquity planets are hotter than their low obliquity counterparts due to ice-albedo feedbacks for cold climates, and water vapor in warm climates. We suggest that the water vapor greenhouse trapping is greater on high obliquity bodies due to the different dynamical regimes that occur between the two states. While equatorial ice-belts are stable at high obliquity in some climate regimes, it is harder to achieve global glaciation than for a low obliquity planet. Temperate polar conditions can be present at high obliquity at forcings for which low obliquity planets would be in a hard snowball state. We suggest the conditions on high obliquity planets are likely to be more favorable for a robust biosphere to develop approaching the outer edge of the HZ. However, the influence of obliquity diminishes for dense atmospheres, in agreement with calculations from 1-D Energy Balance Models.

The current progress in the detection of terrestrial type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles, and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favorable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro) physical, chemical, and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the habitable zone to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

The lunar history of water deposition, loss, and transport post-accretion has become an important consideration in relation to the possibility of a human outpost on the Moon. Very recent work has shown that a secondary primordial atmosphere of up to 10 mb could have been emplaced ~3.5 billion years ago due to volcanic outgassing from the maria. Using a zero dimensional chemistry model we demonstrate the temperature dependence of the resulting major atmospheric components (CO or CO$_2$). We use a three dimensional general circulation model to test the viability of such an atmosphere and derive its climatological characteristics. Based on these results we then conjecture on its capability to transport volatiles outgassed from the maria to the permanently shadowed regions at the poles. Our preliminary results demonstrate that atmospheres as low as 1 mb are viable and that permanent cold trapping of volatiles is only possible at the poles.

The solar tide in an ancient Venusian ocean is simulated using a dedicated numerical tidal model. Simulations with varying ocean depth and rotational periods ranging from -243 to 64 sidereal Earth days are used to calculate the tidal dissipation rates and associated tidal torque. The results show that the tidal dissipation could have varied by more than 5 orders of magnitude, from 0.001--780 Gigawatts, depending on rotational period and ocean depth. The associated tidal torque is about 2 orders of magnitude below the present day Venusian atmospheric torque, and could change the Venusian day length by up to 72 days per million years depending on rotation rate. Consequently, an ocean tide on ancient Venus could have had significant effects on the rotational history of the planet. These calculations have implications for the rotational periods of similarly close-in exoplanetary worlds and the location of the inner edge of the liquid water habitable zone.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

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.

The potential habitability of known exoplanets is often categorized by a nominal equilibrium temperature assuming a Bond albedo of either 0.3, similar to Earth, or 0. As an indicator of habitability, this leaves much to be desired, because albedo on other planets can be very different, and because surface temperature exceeds equilibrium temperature due to the atmospheric greenhouse effect. We use an ensemble of 3-dimensional general circulation model simulations to show that for a range of habitable planets, much of the variability of Bond albedo, equilibrium temperature, and even surface temperature can be predicted with useful accuracy from incident stellar flux and stellar temperature, two known external parameters for every confirmed exoplanet.

The planetary mass and radius sensitivity of exoplanet discovery capabilities has reached into the terrestrial regime. The focus of such investigations is to search within the Habitable Zone where a modern Earth-like atmosphere may be a viable comparison. However, the detection bias of the transit and radial velocity methods lies close to the host star where the received flux at the planet may push the atmosphere into a runaway greenhouse state. One such exoplanet discovery, Kepler-1649b, receives a similar flux from its star as modern Venus does from the Sun, and so was categorized as a possible exoVenus. Here we discuss the planetary parameters of Kepler-1649b with relation to Venus to establish its potential as a Venus analog. We utilize the general circulation model ROCKE-3D to simulate the evolution of the surface temperature of Kepler-1649b under various assumptions, including relative atmospheric abundances. We show that in all our simulations the atmospheric model rapidly diverges from temperate surface conditions towards a runaway greenhouse with rapidly escalating surface temperatures. We calculate transmission spectra for the evolved atmosphere and discuss these spectra within the context of the James Webb Space Telescope (JWST) Near-Infrared Spectrograph (NIRSpec) capabilities. We thus demonstrate the detectability of the key atmospheric signatures of possible runaway greenhouse transition states and outline the future prospects of characterizing potential Venus analogs.

We present a large ensemble of simulations of an Earth-like world with increasing insolation and rotation rate. Unlike previous work utilizing idealized aquaplanet configurations we focus our simulations on modern Earth-like topography. The orbital period is the same as modern Earth, but with zero obliquity and eccentricity. The atmosphere is 1 bar N$_{2}$-dominated with CO$_{2}$=400 ppmv and CH$_{4}$=1 ppmv. The simulations include two types of oceans; one without ocean heat transport (OHT) between grid cells as has been commonly used in the exoplanet literature, while the other is a fully coupled dynamic bathtub type ocean. The dynamical regime transitions that occur as day length increases induce climate feedbacks producing cooler temperatures, first via the reduction of water vapor with increasing rotation period despite decreasing shortwave cooling by clouds, and then via decreasing water vapor and increasing shortwave cloud cooling, except at the highest insolations. Simulations without OHT are more sensitive to insolation changes for fast rotations while slower rotations are relatively insensitive to ocean choice. OHT runs with faster rotations tend to be similar with gyres transporting heat poleward making them warmer than those without OHT. For slower rotations OHT is directed equator-ward and no high latitude gyres are apparent. Uncertainties in cloud parameterization preclude a precise determination of habitability but do not affect robust aspects of exoplanet climate sensitivity. This is the first paper in a series that will investigate aspects of habitability in the simulations presented herein. The datasets from this study are opensource and publicly available.

This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux.

It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years.

We explore the atmospheric and surface history of a hypothetical paleo-Venus climate using a 3-D General Circulation Model. We constrain our model with the in-situ and remote sensing Venus data available today. Given that Venus and Earth are believed to be similar geochemically some aspects of Earth's history are also utilized. We demonstrate that it is possible for ancient Venus and Venus-like exoplanetary worlds to exist within the liquid water habitable zone with insolations up to nearly 2 times that of modern Earth.

The goals of the astrobiology community are focussed on developing a framework for the detection of biosignatures, or evidence thereof, on objects inside and outside of our solar system. A fundamental aspect of understanding the limits of habitable environments and detectable signatures is the study of where the boundaries of such environments can occur. Thus, the need to study the creation, evolution, and frequency of hostile environments for habitability is an integral part of the astrobiology story. These provide the opportunity to understand the bifurcation, between habitable and uninhabitable. The archetype of such a planet is the Earth's sister planet, Venus, and provides a unique opportunity to explore the processes that created a completely uninhabitable environment and thus define the conditions that can rule out bio-related signatures. We advocate a continued comprehensive study of our sister planet, including models of early atmospheres, compositional abundances, and Venus-analog frequency analysis from current and future exoplanet data. Moreover, new missions to Venus to provide in-situ data are necessary.

We explain how GCM boundary condition choices such as ocean-lake coverage-depth, rotation rate, atmospheric constituents and other factors influence surface conditions in ROCKE-3D paleo-Venus simulations. Studies such as these should also be considered when examining liquid water habitability in similar exoplanet experiments. We also describe how one ingests 3D topographic data from NASAs Venus Magellan Spacecraft radar observations into the ROCKE-3D Planetary General Circulation Model.

The nearby exoplanet Proxima Centauri b will be a prime future target for characterization, despite questions about its retention of water. Climate models with static oceans suggest that an Earth-like Proxima b could harbor a small dayside region of surface liquid water at fairly warm temperatures despite its weak instellation. We present the first 3-dimensional climate simulations of Proxima b with a dynamic ocean. We find that an ocean-covered Proxima b could have a much broader area of surface liquid water but at much colder temperatures than previously suggested, due to ocean heat transport and depression of the freezing point by salinity. Elevated greenhouse gas concentrations do not necessarily produce more open ocean area because of possible dynamic regime transitions. For an evolutionary path leading to a highly saline present ocean, Proxima b could conceivably be an inhabited, mostly open ocean planet dominated by halophilic life. For an ocean planet in 3:2 spin-orbit resonance, a permanent tropical waterbelt exists for moderate eccentricity. Simulations of Proxima Centauri b may also be a model for the habitability of planets receiving similar instellation from slightly cooler or warmer stars, e.g., in the TRAPPIST-1, LHS 1140, GJ 273, and GJ 3293 systems.

We demonstrate the effectiveness of a relatively straightforward analysis of the complex 3D Fourier transform of galaxy coordinates derived from redshift surveys. Numerical demonstrations of this approach are carried out on a volume-limited sample of the Sloan Digital Sky Survey redshift survey. The direct unbinned transform yields a complex 3D data cube quite similar to that from the Fast Fourier Transform (FFT) of finely binned galaxy positions. In both cases deconvolution of the sampling window function yields estimates of the true transform. Simple power spectrum estimates from these transforms are roughly consistent with those using more elaborate methods. However we concentrate on the less often studied Fourier phase spectrum, a simple and general framework for characterizing non-Gaussianity, more easily interpretable than the tangled, incomplete multi-point methods conventionally used. No significant signature of non-Gaussianity has been found in the relatively small data set analyzed, but we identify some threads of modern large scale inference methodology that will presumably yield detections in new wider and deeper surveys.

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) is a 3-Dimensional General Circulation Model (GCM) developed at the NASA Goddard Institute for Space Studies for the modeling of atmospheres of Solar System and exoplanetary terrestrial planets. Its parent model, known as ModelE2 (Schmidt et al. 2014), is used to simulate modern and 21st Century Earth and near-term paleo-Earth climates. ROCKE-3D is an ongoing effort to expand the capabilities of ModelE2 to handle a broader range of atmospheric conditions including higher and lower atmospheric pressures, more diverse chemistries and compositions, larger and smaller planet radii and gravity, different rotation rates (slowly rotating to more rapidly rotating than modern Earth, including synchronous rotation), diverse ocean and land distributions and topographies, and potential basic biosphere functions. The first aim of ROCKE-3D is to model planetary atmospheres on terrestrial worlds within the Solar System such as paleo-Earth, modern and paleo-Mars, paleo-Venus, and Saturn's moon Titan. By validating the model for a broad range of temperatures, pressures, and atmospheric constituents we can then expand its capabilities further to those exoplanetary rocky worlds that have been discovered in the past and those to be discovered in the future. We discuss the current and near-future capabilities of ROCKE-3D as a community model for studying planetary and exoplanetary atmospheres.

Rotation in planetary atmospheres plays an important role in regulating atmospheric and oceanic heat flow, cloud formation and precipitation. Using the Goddard Institute for Space Studies (GISS) three dimension General Circulation Model (3D-GCM) we investigate how the effects of varying rotation rate and increasing the incident stellar flux on a planet set bounds on a planet's habitable zone with its parent star. From ensemble climate simulations we identify which factors are the primary controllers of uncertainty in setting these bounds. This is shown in particular for fully coupled ocean (FCO) runs -- some of the first that have been utilized in this context. Results with a Slab Ocean (SO) of 100m mixed layer depth are compared with a similar study by Yang et al. 2014, which demonstrates consistency across models. However, there are clear differences for rotations rates of 1-16x present Earth sidereal day lengths between the 100m SO and FCO models, which points to the necessity of using FCOs whenever possible. The latter was recently demonstrated quite clearly by Hu & Yang 2014 in their aquaworld study with a FCO when compared with similar mixed layer ocean studies and by Cullum et al. 2014. We also show how these results have implications for Venus in the early history of our Solar System since even at this time Venus received more solar flux than Earth does today while it may still have had a slow retrograde rotation. The Venus runs utilize a 2.9Gya solar spectrum generated with the code of Claire et al. 2012, a modern Venus topography with an ocean filling the lowlands (giving an equivalent depth of 310 meters if spread across the entire surface), atmosphere of 1 bar N2, CO2=0.4mb, CH4=0.001mb and present day orbital parameters, radius, & gravity. We demonstrate that ancient Venus could have had quite moderate surface temperatures given these assumptions.

The Kepler era of exoplanetary discovery has presented the Astronomical community with a cornucopia of planetary systems very different from the one which we inhabit. It has long been known that Jupiter plays a major role in the orbital parameters of Mars and it's climate, but there is also a long-standing belief that Jupiter would play a similar role for Earth if not for its large moon. Using a three dimensional general circulation model (3-D GCM) with a fully-coupled ocean we simulate what would happen to the climate of an Earth-like world if Mars did not exist, but a Jupiter-like planet was much closer to Earth's orbit. We investigate two scenarios that involve evolution of the Earth-like planet's orbital eccentricity from 0--0.283 over 6500 years, and from 0--0.066 on a time scale of 4500 years. In both cases we discover that they would maintain relatively temperate climates over the time-scales simulated. More Earth-like planets in multi-planet systems will be discovered as we continue to survey the skies and the results herein show that the proximity of large gas giant planets may play an important role in the habitability of these worlds. These are the first such 3-D GCM simulations using a fully-coupled ocean with a planetary orbit that evolves over time due to the presence of a giant planet.

Present-day Venus is an inhospitable place with surface temperatures approaching 750K and an atmosphere over 90 times as thick as present day Earth's. Billions of years ago the picture may have been very different. We have created a suite of 3D climate simulations using topographic data from the Magellan mission, solar spectral irradiance estimates for 2.9 and 0.715 billion years ago, present day Venus orbital parameters, an ocean volume consistent with current theory and measurements, and an atmospheric composition estimated for early Venus. Using these parameters we find that such a world could have had moderate temperatures if Venus had a rotation period slower than about 16 Earth days, despite an incident solar flux 46-70% higher than modern Earth receives. At its current rotation period of 243 days, Venus's climate could have remained habitable until at least 715 million years ago if it hosted a shallow primordial ocean. These results demonstrate the vital role that rotation and topography play in understanding the climatic history of exoplanetary Venus-like worlds being discovered in the present epoch.

We study the long term orbital evolution of a terrestrial planet under the gravitational perturbations of a giant planet. In particular, we are interested in situations where the two planets are in the same plane and are relatively close. We examine both possible configurations: the giant planet orbit being either outside or inside the orbit of the smaller planet. The perturbing potential is expanded to high orders and an analytical solution of the terrestrial planetary orbit is derived. The analytical estimates are then compared against results from the numerical integration of the full equations of motion and we find that the analytical solution works reasonably well. An interesting finding is that the new analytical estimates improve greatly the predictions for the timescales of the orbital evolution of the terrestrial planet compared to an octupole order expansion. Finally, we briefly discuss possible applications of the analytical estimates in astrophysical problems.

Despite strong interest in the starburst (hereafter SB) phenomenon, the concept remains ill-defined. We use a strict definition of SB to examine the statistical properties of local SB and post-starburst (hereafter PB) galaxies. We also seek relationships to active galaxies. Potential SB galaxies are selected from the SDSS DR7 and their stellar content is analysed. We apply an age dependent dust attenuation correction and derive star formation rates (SFR), ages and masses of the young and old populations. The photometric masses nicely agree with dynamical masses derived from the H-alpha emission line width. To select SB galaxies, we use the birthrate parameter b=SFR/<SFR>, requiring b>=3. The PB sample is selected from the citerion EW(Hdelta_abs)>=6 A. Only 1% of star-forming galaxies are found to be SB galaxies. They contribute 3-6% to the stellar production and are therefore unimportant for the local star formation activity. The median SB age is 70 Myr, roughly independent of mass. The b-parameter strongly depends on burst age. Values close to b=60 are found at ages ~10 Myr, while almost no SBs are found at ages >1 Gyr. The median baryonic burst mass fraction of sub-L* galaxies is 5%, decreasing slowly with mass. The median mass fraction of the recent burst in the PB sample is 5-10%. The age-mass distribution of the progenitors of the PBs is bimodal with a break at log(M)~10.6 above which the ages are doubled. The SB and PB luminosity functions (hereafter LFs) follow each other closely until M_r~-21, when AGNs begin to dominate. The PB LF continues to follow the AGN LF while SB loose significance. This suggests that the number of luminous SBs is underestimated by about one dex at high luminosities, due to large amounts of dust and/or AGN blending. It also indicates that the SB phase preceded the AGN phase. We also discuss the conditions for global gas outflow caused by stellar feedback.

The major uncertainties in studies of the multi-scale structure of the Universe arise not from observational errors but from the variety of legitimate definitions and detection methods for individual structures. To facilitate the study of these methodological dependencies we have carried out 12 different analyses defining structures in various ways. This has been done in a purely geometrical way by utilizing the HOP algorithm as a unique parameter-free method of assigning groups of galaxies to local density maxima or minima. From three density estimation techniques (smoothing kernels, Bayesian Blocks and self organizing maps) applied to three data sets (the Sloan Digital Sky Survey Data Release 7, the Millennium Simulation and randomly distributed points) we tabulate information that can be used to construct catalogs of structures connected to local density maxima and minima. The resulting sizes follow continuous multi-scale distributions with no indication of the presence of a discrete hierarchy. We also introduce a novel void finder that utilizes a method to assemble Delaunay tetrahedra into connected structures and characterizes regions very nearly empty of galaxies in the source catalog.

Edwin Hubble is famous for a number of discoveries that are well known to amateur and professional astronomers, students and the general public. The origins of these discoveries are examined and it is demonstrated that, in each case, a great deal of supporting evidence was already in place. In some cases the discoveries had either already been made, or competing versions were not adopted for complex scientific and sociological reasons.

We present an unsupervised machine learning approach that can be employed for estimating photometric redshifts. The proposed method is based on a vector quantization approach called Self--Organizing Mapping (SOM). A variety of photometrically derived input values were utilized from the Sloan Digital Sky Survey's Main Galaxy Sample, Luminous Red Galaxy, and Quasar samples along with the PHAT0 data set from the PHoto-z Accuracy Testing project. Regression results obtained with this new approach were evaluated in terms of root mean square error (RMSE) to estimate the accuracy of the photometric redshift estimates. The results demonstrate competitive RMSE and outlier percentages when compared with several other popular approaches such as Artificial Neural Networks and Gaussian Process Regression. SOM RMSE--results (using $\Delta$z=z$_{phot}$--z$_{spec}$) for the Main Galaxy Sample are 0.023, for the Luminous Red Galaxy sample 0.027, Quasars are 0.418, and PHAT0 synthetic data are 0.022. The results demonstrate that there are non--unique solutions for estimating SOM RMSEs. Further research is needed in order to find more robust estimation techniques using SOMs, but the results herein are a positive indication of their capabilities when compared with other well-known methods.

Edwin Hubble is often credited with discovering the expanding Universe based on spectra taken by him. This statement is incorrect and we feel that it is the responsibility of those who are aware of the historical facts to set the record straight.

The Large Synoptic Survey Telescope (LSST) simulator being built by Andy Connolly and collaborators is an impressive undertaking and should make working with LSST in the beginning stages far more easy than it was initially with the Sloan Digital Sky Survey. However, I would like to focus on an equally important problem that has not yet been discussed here, but in the coming years the community will need to address -- can we deal with the flood of data from LSST and will we need to rethink the way we work?

The purpose of the New York Workshop on Computer, Earth and Space Sciences is to bring together the New York area's finest Astronomers, Statisticians, Computer Scientists, Space and Earth Scientists to explore potential synergies between their respective fields. The 2011 edition (CESS2011) was a great success, and we would like to thank all of the presenters and participants for attending. This year was also special as it included authors from the upcoming book titled "Advances in Machine Learning and Data Mining for Astronomy". Over two days, the latest advanced techniques used to analyze the vast amounts of information now available for the understanding of our universe and our planet were presented. These proceedings attempt to provide a small window into what the current state of research is in this vast interdisciplinary field and we'd like to thank the speakers who spent the time to contribute to this volume.

It has recently been demonstrated that one can accurately derive galaxy morphology from particular primary and secondary isophotal shape estimates in the Sloan Digital Sky Survey imaging catalog. This was accomplished by applying Machine Learning techniques to the Galaxy Zoo morphology catalog. Using the broad bandpass photometry of the Sloan Digital Sky Survey in combination with with precise knowledge of galaxy morphology should help in estimating more accurate photometric redshifts for galaxies. Using the Galaxy Zoo separation for spirals and ellipticals in combination with Sloan Digital Sky Survey photometry we attempt to calculate photometric redshifts. In the best case we find that the root mean square error for Luminous Red Galaxies classified as ellipticals is as low as 0.0118. Given these promising results we believe better photometric redshift estimates for all galaxies in the Sloan Digital Sky Survey ($\sim$350 million) will be feasible if researchers can also leverage their derived morphologies via Machine Learning. These initial results look to be promising for those interested in estimating Weak-Lensing, Baryonic Acoustic Oscillation, and other fields dependent upon accurate photometric redshifts.

Three methods for detecting and characterizing structure in point data, such as that generated by redshift surveys, are described: classification using self-organizing maps, segmentation using Bayesian blocks, and density estimation using adaptive kernels. The first two methods are new, and allow detection and characterization of structures of arbitrary shape and at a wide range of spatial scales. These methods should elucidate not only clusters, but also the more distributed, wide-ranging filaments and sheets, and further allow the possibility of detecting and characterizing an even broader class of shapes. The methods are demonstrated and compared in application to three data sets: a carefully selected volume-limited sample from the Sloan Digital Sky Survey redshift data, a similarly selected sample from the Millennium Simulation, and a set of points independently drawn from a uniform probability distribution -- a so-called Poisson distribution. We demonstrate a few of the many ways in which these methods elucidate large scale structure in the distribution of galaxies in the nearby Universe.

Scientific data sets continue to increase in both size and complexity. In the past, dedicated graphics systems at supercomputing centers were required to visualize large data sets, but as the price of commodity graphics hardware has dropped and its capability has increased, it is now possible, in principle, to view large complex data sets on a single workstation. To do this in practice, an investigator will need software that is written to take advantage of the relevant graphics hardware. The Viewpoints visualization package described herein is an example of such software. Viewpoints is an interactive tool for exploratory visual analysis of large, high-dimensional (multivariate) data. It leverages the capabilities of modern graphics boards (GPUs) to run on a single workstation or laptop. Viewpoints is minimalist: it attempts to do a small set of useful things very well (or at least very quickly) in comparison with similar packages today. Its basic feature set includes linked scatter plots with brushing, dynamic histograms, normalization and outlier detection/removal. Viewpoints was originally designed for astrophysicists, but it has since been used in a variety of fields that range from astronomy, quantum chemistry, fluid dynamics, machine learning, bioinformatics, and finance to information technology server log mining. In this article, we describe the Viewpoints package and show examples of its usage.

Expanding upon the work of Way and Srivastava 2006 we demonstrate how the use of training sets of comparable size continue to make Gaussian process regression (GPR) a competitive approach to that of neural networks and other least-squares fitting methods. This is possible via new large size matrix inversion techniques developed for Gaussian processes (GPs) that do not require that the kernel matrix be sparse. This development, combined with a neural-network kernel function appears to give superior results for this problem. Our best fit results for the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample using u,g,r,i,z filters gives an rms error of 0.0201 while our results for the same filters in the luminous red galaxy sample yield 0.0220. We also demonstrate that there appears to be a minimum number of training-set galaxies needed to obtain the optimal fit when using our GPR rank-reduction methods. We find that morphological information included with many photometric surveys appears, for the most part, to make the photometric redshift evaluation slightly worse rather than better. This would indicate that most morphological information simply adds noise from the GP point of view in the data used herein. In addition, we show that cross-match catalog results involving combinations of the Two Micron All Sky Survey, SDSS, and Galaxy Evolution Explorer have to be evaluated in the context of the resulting cross-match magnitude and redshift distribution. Otherwise one may be misled into overly optimistic conclusions.

Terrestrial exoplanets are on the verge of joining the ranks of astronomically accessible objects. Interpreting their observable characteristics, and informing decisions on instrument design and use, will hinge on the ability to model these planets successfully across a vast range of configurations and climate forcings. A hierarchical approach that addresses fundamental behaviors as well as more complex, specific, situations is crucial to this endeavor and is presented here. Incorporating Earth-centric knowledge, and continued cross-disciplinary work will be critical, but ultimately the astrophysical study of terrestrial exoplanets must be encouraged to develop as its own field.

We calculate photometric redshifts from the Sloan Digital Sky Survey Main Galaxy Sample, The Galaxy Evolution Explorer All Sky Survey, and The Two Micron All Sky Survey using two new training-set methods. We utilize the broad-band photometry from the three surveys alongside Sloan Digital Sky Survey measures of photometric quality and galaxy morphology. Our first training-set method draws from the theory of ensemble learning while the second employs Gaussian process regression both of which allow for the estimation of redshift along with a measure of uncertainty in the estimation. The Gaussian process models the data very effectively with small training samples of approximately 1000 points or less. These two methods are compared to a well known Artificial Neural Network training-set method and to simple linear and quadratic regression. Our results show that robust photometric redshift errors as low as 0.02 RMS can regularly be obtained. We also demonstrate the need to provide confidence bands on the error estimation made by both classes of models. Our results indicate that variations due to the optimization procedure used for almost all neural networks, combined with the variations due to the data sample, can produce models with variations in accuracy that span an order of magnitude. A key contribution of this paper is to quantify the variability in the quality of results as a function of model and training sample. We show how simply choosing the "best" model given a data set and model class can produce misleading results.

The Southern Abell Redshift Survey contains 39 clusters of galaxies with redshifts in the range 0.0 < z < 0.31 and a median redshift depth of z = 0.0845. SARS covers the region 0 < decl < -65 deg, R.A. < 5h, R.A. > 21h (while avoiding the LMC and SMC) with b > 40. Cluster locations were chosen from the Abell and Abell-Corwin-Olowin catalogs while galaxy positions were selected from the Automatic Plate Measuring Facility galaxy catalog with extinction-corrected magnitudes in the range 15 <= b_j < 19. SARS utilized the Las Campanas 2.5 m duPont telescope, observing either 65 or 128 objects concurrently over a 1.5 sq deg field. New redshifts for 3440 galaxies are reported in the fields of these 39 clusters of galaxies.

The Large Scale Structure in the galaxy distribution is investigated using The First Data Release of the Sloan Digital Sky Survey. Using the Minimal Spanning Tree technique we have extracted sets of filaments, of wall--like structures, of galaxy groups, and of rich clusters from this unique sample. The physical properties of these structures were then measured and compared with the statistical expectations based on the Zel'dovich' theory. The measured characteristics of galaxy walls were found to be consistent with those for a spatially flat $\Lambda$CDM cosmological model with $\Omega_m\approx$ 0.3 and $\Omega_\Lambda \approx$ 0.7, and for Gaussian initial perturbations with a Harrison -- Zel'dovich power spectrum. Furthermore, we found that the mass functions of groups and of unrelaxed structure elements generally fit well with the expectations from Zel'dovich' theory. We also note that both groups and rich clusters tend to prefer the environments of walls, which tend to be of higher density, rather than the environments of filaments, which tend to be of lower density.

We present non-Gaussianity testing on recently-released derived maps from the first-year WMAP data by Tegmark, de Oliveria-Costa and Hamilton. Our test is based on a phase mapping technique which has the advantage of testing non-Gaussianity at separate multipole bands. We show that their foreground-cleaned map is against the random-phase hypothesis at all 4 multipole bands centered around l = 150, 290, 400 and 500. Their Wiener-filtered map, on the other hand, is Gaussian for l < 250, and marginally Gaussian for 224 < l < 350. However, we see the evidence of non-Gaussianity for l > 350 as we detect certain degrees of phase coupling, hence against the random-phase hypothesis. Our phase mapping technique is particularly useful for testing the accuracy of component separation methods.

We report an analysis of the dynamical structure of clusters of galaxies from a survey of photometric and spectroscopic observations in the fields of southern Abell Clusters. We analyze the galaxy velocity field in extended regions up to 7 h^-1 Mpc from cluster centers and we estimate mean velocity dispersions and their radial dependence. Only one from a total number of 41 Abell clusters does not correspond to a dynamically bound system. However, four of these bound objects are double clusters. We estimate that 20 % (7 clusters) of the 35 remaining are subject to serious projection effects. Normalizing the clustercentric distances by means of the overdensity radius r_200, and the velocity dispersion profiles (VDPs) by the corresponding mean cluster velocity dispersion, we computed the average VDP. Our results indicate a flat behavior of the mean VDP at large distances from the cluster center. Nevertheless, we found that for the inner part of the clusters (r/r_200 leq 1) the VDP is up to a 10% smaller than at larger radii.

We present results of simple N-body simulations that strengthen the suggestion that A3266 is composed of two subunits of comparable mass that have recently merged. Both the real cluster and the N-body dark-matter cluster show mixed signals of substructure under statistical tests. However, in a decidedly non-statistical approach allowed by the wide-area coverage and large number of redshifts they measured in A3266, Quintana, Ramirez, and Way (1996; QRW) sliced the real cluster in redshift space to uncover a peculiar spatial distribution of galaxies that they suggested was the result of a recent merger. In our simulations, a similar distribution is the result of an ongoing merger between two comparable-mass units that started about 2x10^9 years ago in the N-body simulations. We also find that the distribution of emission line galaxies in A3266 traces the same structure. We discuss further tests of our merger hypothesis, and speculate on the possibility that a similar process might be occurring in other, apparently-relaxed clusters at the present epoch.

We detect four isolated, X-ray over-luminous (Lx>2e43 [h/0.5]**-2 erg/s) elliptical galaxies (OLEGs) in our 160 square degree ROSAT PSPC survey. The extent of their X-ray emission, total X-ray luminosity, total mass, and mass of the hot gas in these systems correspond to poor clusters, and the optical luminosity of the central galaxies (M_R<-22.5 + 5 lg h) is comparable to that of cluster cDs. However, there are no detectable fainter galaxy concentrations around the central elliptical. The mass-to-light ratio within the radius of detectable X-ray emission is in the range 250-450 Msun/Lsun, which is 2-3 times higher than typically found in clusters or groups. These objects can be the result of galaxy merging within a group. However, their high M/L values are difficult to explain in this scenario. OLEGs must have been undisturbed for a very long time, which makes them the ultimate examples of systmes in hydrostatic equilibrium. The number density of OLEGs is n=2.4(+3.1-1.2}x10**-7 (h/0.5)**-3 Mpc**-3 at the 90% confidence. They comprise 20% of all clusters and groups of comparable X-ray luminosity, and nearly all galaxies brighter than M_R=-22.5. The estimated contirubution of OLEGs to the total mass density in the Universe is close to that of T>7 keV clusters.

A dynamical analysis of the structure of the clusters of galaxies Abell 119 and Abell 133 is presented, using new redshift data combined with existing data from the literature. We compare our results with those from the X-ray data for these clusters, and with those from radio data for the central cD galaxy in each cluster. A comparison of the mass estimate based on X-ray data and that obtained here after subgroups are eliminated shows them to be comparable. After the elimination of subgroups, 125 galaxy members in Abell 119 and 120 in Abell 133 give dispersions of 472 km/s and 735 km/s respectively. However, our dynamical analysis of the optical data shows little substructure in the velocity field of Abell 133, conflicting with what is seen in the Rosat X-ray map. Abell 119 seems to have multiple structures along the line of sight. We derive virial mass estimates of 3.05 X 10E14 Solar Masses for Abell 119 and 7.79 X 10E14 Solar Masses for Abell 133 within 1.5 Mpc, which agree well with the X-ray-derived masses within errors.

Using spectrophotometry of a large sample of galaxies in 19 Abell Clusters, we have selected 42 candidate AGN by the criteria used by Dressler et.al. 1985 (DTS) in their analysis of the statistics of 22 AGN in 14 rich cluster fields, which are based on the equivalent width of [OII]3727, Hbeta, and [OIII]5007 emission. We have then discriminated AGN from HII region-like galaxies (hereafter HII galaxies) in the manner developed by Veilleux & Osterbrok 1987 using the additional information provided by Halpha and [NII]6583 or Halpha and [SII]6716+6731 emission, in order to test the reliability of the selection criteria used by DTS. We find that our sample is very similar to that of DTS before we discriminate AGN from HII galaxies, and would lead to similar conclusions. However, we find that their method inevitably mixes HII galaxies with AGN, even for the most luminous objects in our sample. Other authors have attempted to quantify the relative fraction of cluster to field AGN since the study of DTS (Hill & Oegerle 1993, Biviano et al. 1997) reaching similar conclusions, but using similar criteria to theirs to select AGN (or using the [OIII]5007/Hbeta flux ratio test that also mixes HII galaxies with AGN). Our sample of true AGN is left too small to reach statistically meaningful conclusions, therefore a new study with the more time-consuming method that includes the other lines will be required to quantify the true relative fraction of cluster to field AGN.

We present magnitudes and redshifts of galaxies in the cluster of galaxies CL0017-20 and in the double cluster CL0500-24. Luminosity functions and velocity distributions have been derived. CL0017-20 possesses a virial mass of $2.9 \times 10^{14}$ M$_\odot$ and a \vmag band luminosity of $2.3 \pm 0.5 \times 10^{12}$ L$_\odot$ ($M/L_{\rm V} = 127$), much of it in a very compact core of size 113 kpc. Our new velocities and photometry of CL0500-24 confirm its double nature. Our group algorithm clearly detects two subconcentrations with masses 1.9 and $2.2 \times 10^{14}$ M$_{\odot}$. Its total \vmag band luminosity amounts to $2.7 \pm 0.8 \times 10^{12}$ L$_{\odot}$ ($M/L_{\rm V}$ of 157). Since the crossing and collapse times in both clusters are smaller than the age of the Universe, they appear bound and virialized. It is surprising that although $M/L$ is small in these clusters, giant arcs have been detected in both.