Search

The two innermost planets of the Proxima Centauri system are separated by just 0.02 au, inducing strong gravitational interactions between them. We assess this interaction by leveraging fast orbital stability indicators and find that orbital stability is very likely if the initial eccentricities of planets b and d are less than $\sim 0.2$, but cannot confirm stability at larger values. We find that stability is not strongly affected by the true masses of the planets or by the distant planet c. However, mutual inclinations between 95$^\circ$ and 142$^\circ$ often result in unstable motion. We further explore the long-term evolution of the orbits in these stable regions of parameter space and find that circularization can take over 5 Gyr. This tidal evolution could support surface energy fluxes in excess of 1 W m$^{-2}$ for over 1 Gyr, possibly affecting planet b's habitability., Comment: 15 pages, 7 figures, 4 tables. Accepted for publication in ApJ

The atmospheric dynamics of tidally locked hot Jupiters is characterized by strong equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photochemical (horizontal and vertical) transport model, VULCAN 2D, which improves on the pseudo-2D approaches by allowing for nonuniform zonal winds. We extensively validate our VULCAN 2D with analytical solutions and benchmark comparisons. Applications to HD 189733 b and HD 209458 b reveal a transition in mixing regimes: horizontal transport predominates below ∼0.1 mbar, while vertical mixing is more important at higher altitudes above 0.1 mbar. Motivated by the previously inferred carbon-rich atmosphere, we find that HD 209458 b with supersolar carbon-to-oxygen ratio (C/O) exhibits pronounced C2H4 absorption on the morning limb but not on the evening limb, due to horizontal transport from the nightside. We discuss when a pseudo-2D approach is a valid assumption and its inherent limitations. Finally, we demonstrate the effect of horizontal transport in transmission observations and its impact on the morning−evening limb asymmetry with synthetic spectra, highlighting the need to consider global transport when interpreting exoplanet atmospheres.

Rotational dynamics of the Earth, over geological timescales, have profoundly affected local and global climatic evolution, probably contributing to the evolution of life. To better retrieve the Earth's rotational history, and motivated by the published hypothesis of a stabilized length of day during the Precambrian, we examine the effect of thermal tides on the evolution of planetary rotational motion. The hypothesized scenario is contingent upon encountering a resonance in atmospheric Lamb waves, whereby an amplified thermotidal torque cancels the opposing torque of the oceans and solid interior, driving the Earth into a rotational equilibrium. With this scenario in mind, we construct an ab initio model of thermal tides on rocky planets describing a neutrally stratified atmosphere. The model takes into account dissipative processes with Newtonian cooling and diffusive processes in the planetary boundary layer. We retrieve from this model a closed-form solution for the frequency-dependent tidal torque which captures the main spectral features previously computed using 3D general circulation models. In particular, under longwave heating, diffusive processes near the surface and the delayed thermal response of the ground prove to be responsible for attenuating, and possibly annihilating, the accelerating effect of the thermotidal torque at the resonance. When applied to the Earth, our model prediction suggests the occurrence of the Lamb resonance in the Phanerozoic, but with an amplitude that is insufficient for the rotational equilibrium. Interestingly, though our study was motivated by the Earth's history, the generic tidal solution can be straightforwardly and efficiently applied in exoplanetary settings., Comment: 20 pages (+14 for appendices), 6 figures

Global circulation models (GCMs) play an important role in contemporary investigations of exoplanet atmospheres. Different GCMs evolve various sets of dynamical equations which can result in obtaining different atmospheric properties between models. In this study, we investigate the effect of different dynamical equation sets on the atmospheres of hot Jupiter exoplanets. We compare GCM simulations using the quasi-primitive dynamical equations (QHD) and the deep Navier-Stokes equations (NHD) in the GCM THOR. We utilise a two-stream non-grey "picket-fence" scheme to increase the realism of the radiative transfer calculations. We perform GCM simulations covering a wide parameter range grid of system parameters in the population of exoplanets. Our results show significant differences between simulations with the NHD and QHD equation sets at lower gravity, higher rotation rates or at higher irradiation temperatures. The chosen parameter range shows the relevance of choosing dynamical equation sets dependent on system and planetary properties. Our results show the climate states of hot Jupiters seem to be very diverse, where exceptions to prograde superrotation can often occur. Overall, our study shows the evolution of different climate states which arise just due to different selections of Navier-Stokes equations and approximations. We show the divergent behaviour of approximations used in GCMs for Earth, but applied for non Earth-like planets., Comment: accepted MNRAS (21 June 2023)

Molecular oxygen in our atmosphere has increased from less than a part per million in the Archean Eon, to a fraction of a percent in the Proterozoic, and finally to modern levels during the Phanerozoic. The ozone layer formed with the early Proterozoic oxygenation. While oxygen itself has only minor radiative and climatic effects, the accompanying ozone has important consequences for Earth climate. Using the Community Earth System Model (CESM), a 3-D general circulation model, we test the effects of various levels of ozone on Earth's climate. When CO2 is held constant, the global mean surface temperature decreases with decreasing ozone, with a maximum drop of ~3.5 K at near total ozone removal. By supplementing our GCM results with 1-D radiative flux calculations, we are able to test which changes to the atmosphere are responsible for this temperature change. We find that the surface temperature change is caused mostly by the stratosphere being much colder when ozone is absent; this makes it drier, substantially weakening the greenhouse effect. We also examine the effect of the structure of the upper troposphere and lower stratosphere on the formation of clouds, and on the global circulation. At low ozone, both high and low clouds become more abundant, due to changes in the tropospheric stability. These generate opposing short-wave and long-wave radiative forcings that are nearly equal. The Hadley circulation and tropospheric jet streams are strengthened, while the stratospheric polar jets are weakened, the latter being a direct consequence of the change in stratospheric temperatures. This work identifies the major climatic impacts of ozone, an important piece of the evolution of Earth's atmosphere., Comment: 20 pages, 10 figures, accepted for publication at Climate of the Past

Energy balance models (EBMs) are one- or two-dimensional climate models that can provide insight into planetary atmospheres, particularly with regard to habitability. Because EBMs are far less computationally intensive than three-dimensional general circulation models (GCMs), they can be run over large, uncertain parameter spaces and can be used to explore long-period phenomena like carbon and Milankovitch cycles. Because horizontal dimensions are incorporated in EBMs, they can explore processes that are beyond the reach of one-dimensional radiative-convective models (RCMs). EBMs are, however, dependent on parameterizations and tunings to account for physical processes that are neglected. Thus, EBMs rely on observations and results from GCMs and RCMs. Different EBMs have included a wide range of parameterizations (for albedo, radiation, and heat diffusion) and additional physics, such as carbon cycling and ice sheets. This CUISINES exoplanet model intercomparison project (exoMIP) will compare various EBMs across a set of numerical experiments. The set of experiments will include Earth-like planets at different obliquities, parameter sweeps across obliquity, and variations in instellation and CO$_2$ abundance to produce hysteresis diagrams. We expect a range of different results due to the choices made in the various codes, highlighting which results are robust across models and which are dependent on parameterizations or other modeling choices. Additionally, it will allow developers to identify model defects and determine which parameterizations are most useful or relevant to the problem of interest. Ultimately, this exoMIP will allow us to improve the consistency between EBMs and accelerate the process of discovering habitable exoplanets., Comment: 19 pages, 5 figures, accepted for publication at Planetary Science Journal

The atmospheric dynamics of tidally-locked hot Jupiters is characterized by strong equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photochemical (horizontal and vertical) transport model, VULCAN 2D, which improves on the pseudo-2D approaches by allowing for non-uniform zonal winds. We extensively validate our VULCAN 2D with analytical solutions and benchmark comparisons. Applications to HD 189733 b and HD 209458 b reveal a transition in mixing regimes: horizontal transport predominates below $\sim$0.1 mbar while vertical mixing is more important at higher altitudes above 0.1 mbar. Motivated by the previously inferred carbon-rich atmosphere, we find that HD 209458 b with super-solar carbon-to-oxygen ratio (C/O) exhibits pronounced C$_2$H$_4$ absorption on the morning limb but not on the evening limb, owing to horizontal transport from the nightside. We discuss when a pseudo-2D approach is a valid assumption and its inherent limitations. Finally, we demonstrate the effect of horizontal transport in transmission observations and its impact on the morning-evening limb asymmetry with synthetic spectra, highlighting the need to consider global transport when interpreting exoplanet atmospheres., Comment: 19 pages, 20 figures, published in ApJ

General circulation models (GCMs) provide context for interpreting multi-wavelength, multi-phase data of the atmospheres of tidally locked exoplanets. In the current study, the non-hydrostatic THOR GCM is coupled with the HELIOS radiative transfer solver for the first time, supported by an equilibrium chemistry solver (FastChem), opacity calculator (HELIOS-K) and Mie scattering code (LX-MIE). To accurately treat the scattering of radiation by medium-sized to large aerosols/condensates, improved two-stream radiative transfer is implemented within a GCM for the first time. Multiple scattering is implemented using a Thomas algorithm formulation of the two-stream flux solutions, which decreases the computational time by about 2 orders of magnitude compared to the iterative method used in past versions of HELIOS. As a case study, we present four GCMs of the hot Jupiter WASP-43b, where we compare the temperature, velocity, entropy, and streamfunction, as well as the synthetic spectra and phase curves, of runs using regular versus improved two-stream radiative transfer and isothermal versus non-isothermal layers. While the global climate is qualitatively robust, the synthetic spectra and phase curves are sensitive to these details. A THOR+HELIOS WASP-43b GCM (horizontal resolution of about 4 degrees on the sphere and with 40 radial points) with multi-wavelength radiative transfer (30 k-table bins) running for 3000 Earth days (864,000 time steps) takes about 19-26 days to complete depending on the type of GPU., Comment: 31 pages, 24 figures, accepted for publication at MNRAS

Rocky planets hosted by close-in extrasolar systems are likely to be tidally locked in 1:1 spin-orbit resonance, a configuration where they exhibit permanent dayside and nightside. Because of the resulting day-night temperature gradient, the climate and large-scale circulation of these planets are strongly determined by their atmospheric stability against collapse, which designates the runaway condensation of greenhouse gases. To better constrain the surface conditions of rocky planets located in the habitable zone of their host star, it is therefore crucial to elucidate the mechanisms that govern the day-night heat redistribution. As a first attempt to bridge the gap between multiple modelling approaches ranging from idealised models to 3-D General Circulation Models (GCM), we developed a General Circulation Meta-Model (GCMM) able to reproduce both the closed-form solutions provided by analytical models and the numerical solutions obtained from GCM simulations. We used this approach to characterise the atmospheric stability of Earth-sized rocky planets with dry atmospheres containing CO2, and we benchmarked it against 3-D GCM simulations using THOR GCM. We observe that the collapse pressure below which collapse occurs can vary by ~40% around the value predicted by analytical scaling laws depending on the mechanisms taken into account among radiative transfer, atmospheric dynamics, and turbulent diffusion. Particularly, we find (i) that the turbulent diffusion taking place in the dayside planetary boundary layer (PBL) globally tends to warm up the nightside surface hemisphere except in the transition zone between optically thin and optically thick regimes, (ii) that the PBL also significantly affects the day-night advection timescale, and (iii) that the slow rotator approximation holds from the moment that the normalised equatorial Rossby deformation radius is greater than 2., Comment: 32 pages, 16 figures, accepted for publication in Astronomy & Astrophysics

With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the exoplanetary science community to simulate the atmospheres of mini-Neptunes. We focus on two targets well studied both observationally and theoretically with planned JWST Cycle 1 observations: the warm GJ~1214b and the cooler K2-18b. For each target, we consider a temperature-forced case, a clear sky dual-grey radiative transfer case, and a clear sky multi band radiative transfer case, covering a range of complexities and configurations where we know differences exist between GCMs in the literature. This paper presents all the details necessary to participate in the intercomparison, with the intention of presenting the results in future papers. Currently, there are eight GCMs participating (ExoCAM, Exo-FMS, FMS PCM, Generic PCM, MITgcm, RM-GCM, THOR, and the UM), and membership in the project remains open. Those interested in participating are invited to contact the authors., Comment: Accepted to PSJ

Rocky planets hosted by close-in extrasolar systems are likely to be tidally locked in 1:1 spin-orbit resonance, a configuration where they exhibit permanent dayside and nightside. Because of the resulting day-night temperature gradient, the climate and large-scale circulation of these planets are strongly determined by their atmospheric stability against collapse, which designates the runaway condensation of greenhouse gases. To better constrain the surface conditions of rocky planets located in the habitable zone of their host star, it is therefore crucial to elucidate the mechanisms that govern the day-night heat redistribution. As a first attempt to bridge the gap between multiple modelling approaches ranging from idealised models to 3-D General Circulation Models (GCM), we developed a General Circulation Meta-Model (GCMM) able to reproduce both the closed-form solutions provided by analytical models and the numerical solutions obtained from GCM simulations. We used this approach to characterise the atmospheric stability of Earth-sized rocky planets with dry atmospheres containing CO2, and we benchmarked it against 3-D GCM simulations using THOR GCM. We observe that the collapse pressure below which collapse occurs can vary by ~40% around the value predicted by analytical scaling laws depending on the mechanisms taken into account among radiative transfer, atmospheric dynamics, and turbulent diffusion. Particularly, we find (i) that the turbulent diffusion taking place in the dayside planetary boundary layer (PBL) globally tends to warm up the nightside surface hemisphere except in the transition zone between optically thin and optically thick regimes, (ii) that the PBL also significantly affects the day-night advection timescale, and (iii) that the slow rotator approximation holds from the moment that the normalised equatorial Rossby deformation radius is greater than 2., Comment: 32 pages, 16 figures, accepted for publication in Astronomy & Astrophysics

General circulation models (GCMs) provide context for interpreting multi-wavelength, multi-phase data of the atmospheres of tidally locked exoplanets. In the current study, the non-hydrostatic THOR GCM is coupled with the HELIOS radiative transfer solver for the first time, supported by an equilibrium chemistry solver (FastChem), opacity calculator (HELIOS-K) and Mie scattering code (LX-MIE). To accurately treat the scattering of radiation by medium-sized to large aerosols/condensates, improved two-stream radiative transfer is implemented within a GCM for the first time. Multiple scattering is implemented using a Thomas algorithm formulation of the two-stream flux solutions, which decreases the computational time by about 2 orders of magnitude compared to the iterative method used in past versions of HELIOS. As a case study, we present four GCMs of the hot Jupiter WASP-43b, where we compare the temperature, velocity, entropy, and streamfunction, as well as the synthetic spectra and phase curves, of runs using regular versus improved two-stream radiative transfer and isothermal versus non-isothermal layers. While the global climate is qualitatively robust, the synthetic spectra and phase curves are sensitive to these details. A THOR+HELIOS WASP-43b GCM (horizontal resolution of about 4 degrees on the sphere and with 40 radial points) with multi-wavelength radiative transfer (30 k-table bins) running for 3000 Earth days (864,000 time steps) takes about 19-26 days to complete depending on the type of GPU., Comment: 31 pages, 24 figures, accepted for publication at MNRAS

General circulation models (GCMs) provide context for interpreting multiwavelength, multiphase data of the atmospheres of tidally locked exoplanets. In the current study, the non-hydrostatic THOR GCM is coupled with the HELIOS radiative transfer solver for the first time, supported by an equilibrium chemistry solver (FastChem), opacity calculator (HELIOS-K), and Mie scattering code (LX-MIE). To accurately treat the scattering of radiation by medium-sized to large aerosols/condensates, improved two-stream radiative transfer is implemented within a GCM for the first time. Multiple scattering is implemented using a Thomas algorithm formulation of the two-stream flux solutions, which decreases the computational time by about 2 orders of magnitude compared to the iterative method used in past versions of HELIOS. As a case study, we present four GCMs of the hot Jupiter WASP-43b, where we compare the temperature, velocity, entropy, and streamfunction, as well as the synthetic spectra and phase curves, of runs using regular versus improved two-stream radiative transfer and isothermal versus non-isothermal layers. While the global climate is qualitatively robust, the synthetic spectra and phase curves are sensitive to these details. A THOR + HELIOS WASP-43b GCM (horizontal resolution of about 4 deg on the sphere and with 40 radial points) with multiwavelength radiative transfer (30 k-table bins) running for 3000 Earth days (864 000 time-steps) takes about 19–26 d to complete depending on the type of GPU

The photometric and spectroscopic signatures of habitable planets orbiting FGK stars may be modulated by surface ice coverage. To estimate its frequency and locations, we simulated the climates of hypothetical planets with a 1D energy balance model and assumed that the planets possess properties similar to modern Earth (mass, geography, atmosphere). We first simulated planets with fixed rotational axes and circular orbits, finding that the vast majority (>90%) of planets with habitable surfaces are free of ice. For planets with partial ice coverage, the parameter space for ice caps (interannual ice located at the poles) is about as large as that for "ice belts" (interannual ice located at the equator), but belts only persist on land. We then performed simulations that mimicked perturbations from other planets by forcing sinusoidal orbital and rotational oscillations over a range of frequencies and amplitudes. We assume initially ice-free surfaces and set the initial eccentricity distribution to mirror known exoplanets, while the initial obliquity distribution matches planet formation predictions, ie favoring 90 degrees. For these dynamic cases, we find again that ~90% of habitable planets are free of surface ice for a range of assumptions for ice's albedo. Planets orbiting F dwarfs are three times as likely to have ice caps than belts, but for planets orbiting K and G dwarfs ice belts are twice as likely as caps. In some cases, a planet's surface ice can cycle between the equatorial and polar regions. Future direct imaging surveys of habitable planets may be able to test these predictions., Comment: 23 pages, 13 figures, 4 tables, accepted to Planetary Science Journal. Source code available at https://github.com/VirtualPlanetaryLaboratory/vplanet, and scripts to generate data and figures available at https://github.com/caitlyn-wilhelm/IceCoverage

Radiative-transfer (RT) is a key component for investigating atmospheres of planetary bodies. With the 3D nature of exoplanet atmospheres being important in giving rise to their observable properties, accurate and fast 3D methods are required to be developed to meet future multi-dimensional and temporal data sets. We develop an open source GPU RT code, gCMCRT, a Monte Carlo RT forward model for general use in planetary atmosphere RT problems. We aim to automate the post-processing pipeline, starting from direct global circulation model (GCM) output to synthetic spectra. We develop albedo, emission and transmission spectra modes for 3D and 1D input structures. We include capability to use correlated-k and high-resolution opacity tables, the latter of which can be Doppler shifted inside the model. We post-process results from several GCM groups including ExoRad, SPARC/MITgcm THOR, UK Met Office UM, Exo-FMS and the Rauscher model. Users can therefore take advantage of desktop and HPC GPU computing solutions. gCMCRT is well suited for post-processing large GCM model grids produced by members of the community and for high-resolution 3D investigations., Comment: 20 pages, 9 figures, submitted to ApJ 29/10/21, revised post-reviewer comments 18/3/22

General circulation models are often used to explore exoclimate parameter spaces and classify atmospheric circulation regimes. Models are tuned to give reasonable climate states for standard test cases, such as the Held-Suarez test, and then used to simulate diverse exoclimates by varying input parameters such as rotation rates, instellation, atmospheric optical properties, frictional timescales and so on. In such studies, there is an implicit assumption that the model which works reasonably well for the standard test case will be credible at all points in an arbitrarily wide parameter space. Here, we test this assumption using the open-source general circulation model THOR to simulate atmospheric circulation on tidally locked Earth-like planets with rotation periods of 0.1 to 100 days. We find that the model error, as quantified by the ratio between physical and spurious numerical contributions to the angular momentum balance, is extremely variable across this range of rotation periods with some cases where numerical errors are the dominant component. Increasing model grid resolution does improve errors but using a higher-order numerical diffusion scheme can sometimes magnify errors for finite-volume dynamical solvers. We further show that to minimize error and make the angular momentum balance more physical within our model, the surface friction timescale must be smaller than the rotational timescale., Comment: 8 pages, 5 figures, accepted by ApJ

The era of atmospheric characterization of terrestrial exoplanets is just around the corner. Modeling prior to observations is crucial in order to predict the observational challenges and to prepare for the data interpretation. This paper presents the report of the TRAPPIST Habitable Atmosphere Intercomparison workshop (2020 September 14-16). A review of the climate models and parameterizations of the atmospheric processes on terrestrial exoplanets, model advancements, and limitations, as well as direction for future model development, was discussed. We hope that this report will be used as a roadmap for future numerical simulations of exoplanet atmospheres and maintaining strong connections to the astronomical community.

The era of atmospheric characterization of terrestrial exoplanets is just around the corner. Modeling prior to observations is crucial in order to predict the observational challenges and to prepare for the data interpretation. This paper presents the report of the TRAPPIST Habitable Atmosphere Intercomparison workshop (2020 September 14-16). A review of the climate models and parameterizations of the atmospheric processes on terrestrial exoplanets, model advancements, and limitations, as well as direction for future model development, was discussed. We hope that this report will be used as a roadmap for future numerical simulations of exoplanet atmospheres and maintaining strong connections to the astronomical community.

The era of atmospheric characterization of terrestrial exoplanets is just around the corner. Modeling prior to observations is crucial in order to predict the observational challenges and to prepare for the data interpretation. This paper presents the report of the TRAPPIST Habitable Atmosphere Intercomparison workshop (2020 September 14-16). A review of the climate models and parameterizations of the atmospheric processes on terrestrial exoplanets, model advancements, and limitations, as well as direction for future model development, was discussed. We hope that this report will be used as a roadmap for future numerical simulations of exoplanet atmospheres and maintaining strong connections to the astronomical community.

The era of atmospheric characterization of terrestrial exoplanets is just around the corner. Modeling prior to observations is crucial in order to predict the observational challenges and to prepare for the data interpretation. This paper presents the report of the TRAPPIST Habitable Atmosphere Intercomparison workshop (2020 September 14-16). A review of the climate models and parameterizations of the atmospheric processes on terrestrial exoplanets, model advancements, and limitations, as well as direction for future model development, was discussed. We hope that this report will be used as a roadmap for future numerical simulations of exoplanet atmospheres and maintaining strong connections to the astronomical community.

General circulation models (GCMs) are often used to explore exoclimate parameter spaces and classify atmospheric circulation regimes. Models are tuned to give reasonable climate states for standard test cases, such as the Held–Suarez test, and then used to simulate diverse exoclimates by varying input parameters such as rotation rates, instellation, atmospheric optical properties, frictional timescales, and so on. In such studies, there is an implicit assumption that the model works reasonably well for the standard test case will be credible at all points in an arbitrarily wide parameter space. Here, we test this assumption using the open-source GCM THOR to simulate atmospheric circulation on tidally locked Earth-like planets with rotation periods of 0.1–100 days. We find that the model error, as quantified by the ratio between physical and spurious numerical contributions to the angular momentum balance, is extremely variable across this range of rotation periods with some cases where numerical errors are the dominant component. Increasing model grid resolution does improve errors, but using a higher-order numerical diffusion scheme can sometimes magnify errors for finite-volume dynamical solvers. We further show that to minimize error and make the angular momentum balance more physical within our model, the surface friction timescale must be smaller than the rotational timescale.

The era of atmospheric characterization of terrestrial exoplanets is just around the corner. Modeling prior to observations is crucial in order to predict the observational challenges and to prepare for the data interpretation. This paper presents the report of the TRAPPIST Habitable Atmosphere Intercomparison workshop (2020 September 14-16). A review of the climate models and parameterizations of the atmospheric processes on terrestrial exoplanets, model advancements, and limitations, as well as direction for future model development, was discussed. We hope that this report will be used as a roadmap for future numerical simulations of exoplanet atmospheres and maintaining strong connections to the astronomical community.

The photometric and spectroscopic signatures of habitable planets orbiting FGK stars may be modulated by surface ice coverage. To estimate its frequency and locations, we simulated the climates of hypothetical planets with a 1D energy balance model and assumed that the planets possess properties similar to modern Earth (mass, geography, atmosphere). We first simulated planets with fixed rotational axes and circular orbits, finding that the vast majority (>90%) of planets with habitable surfaces are free of ice. For planets with partial ice coverage, the parameter space for ice caps (interannual ice located at the poles) is about as large as that for "ice belts" (interannual ice located at the equator), but belts only persist on land. We then performed simulations that mimicked perturbations from other planets by forcing sinusoidal orbital and rotational oscillations over a range of frequencies and amplitudes. We assume initially ice-free surfaces and set the initial eccentricity distribution to mirror known exoplanets, while the initial obliquity distribution matches planet formation predictions, ie favoring 90 degrees. For these dynamic cases, we find again that ~90% of habitable planets are free of surface ice for a range of assumptions for ice's albedo. Planets orbiting F dwarfs are three times as likely to have ice caps than belts, but for planets orbiting K and G dwarfs ice belts are twice as likely as caps. In some cases, a planet's surface ice can cycle between the equatorial and polar regions. Future direct imaging surveys of habitable planets may be able to test these predictions., Comment: 23 pages, 13 figures, 4 tables, accepted to Planetary Science Journal. Source code available at https://github.com/VirtualPlanetaryLaboratory/vplanet, and scripts to generate data and figures available at https://github.com/caitlyn-wilhelm/IceCoverage

General circulation models are often used to explore exoclimate parameter spaces and classify atmospheric circulation regimes. Models are tuned to give reasonable climate states for standard test cases, such as the Held-Suarez test, and then used to simulate diverse exoclimates by varying input parameters such as rotation rates, instellation, atmospheric optical properties, frictional timescales and so on. In such studies, there is an implicit assumption that the model which works reasonably well for the standard test case will be credible at all points in an arbitrarily wide parameter space. Here, we test this assumption using the open-source general circulation model THOR to simulate atmospheric circulation on tidally locked Earth-like planets with rotation periods of 0.1 to 100 days. We find that the model error, as quantified by the ratio between physical and spurious numerical contributions to the angular momentum balance, is extremely variable across this range of rotation periods with some cases where numerical errors are the dominant component. Increasing model grid resolution does improve errors but using a higher-order numerical diffusion scheme can sometimes magnify errors for finite-volume dynamical solvers. We further show that to minimize error and make the angular momentum balance more physical within our model, the surface friction timescale must be smaller than the rotational timescale., Comment: 8 pages, 5 figures, accepted by ApJ

Radiative-transfer (RT) is a key component for investigating atmospheres of planetary bodies. With the 3D nature of exoplanet atmospheres being important in giving rise to their observable properties, accurate and fast 3D methods are required to be developed to meet future multi-dimensional and temporal data sets. We develop an open source GPU RT code, gCMCRT, a Monte Carlo RT forward model for general use in planetary atmosphere RT problems. We aim to automate the post-processing pipeline, starting from direct global circulation model (GCM) output to synthetic spectra. We develop albedo, emission and transmission spectra modes for 3D and 1D input structures. We include capability to use correlated-k and high-resolution opacity tables, the latter of which can be Doppler shifted inside the model. We post-process results from several GCM groups including ExoRad, SPARC/MITgcm THOR, UK Met Office UM, Exo-FMS and the Rauscher model. Users can therefore take advantage of desktop and HPC GPU computing solutions. gCMCRT is well suited for post-processing large GCM model grids produced by members of the community and for high-resolution 3D investigations., Comment: 20 pages, 9 figures, submitted to ApJ 29/10/21, revised post-reviewer comments 18/3/22

The era of atmospheric characterization of terrestrial exoplanets is just around the corner. Modeling prior to observations is crucial in order to predict the observational challenges and to prepare for the data interpretation. This paper presents the report of the TRAPPIST Habitable Atmosphere Intercomparison (THAI) workshop (14-16 September 2020). A review of the climate models and parameterizations of the atmospheric processes on terrestrial exoplanets, model advancements and limitations, as well as direction for future model development was discussed. We hope that this report will be used as a roadmap for future numerical simulations of exoplanet atmospheres and maintaining strong connections to the astronomical community., Comment: Accepted in the AAS Planetary Science Journal (PSJ)

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., Comment: 22 pages, 9 figures, 2 tables, accepted for publication in the Astrophysical Journal

Weathering of silicate rocks at a planetary surface can draw down CO$_2$ from the atmosphere for eventual burial and long-term storage in the planetary interior. This process is thought to provide an essential negative feedback to the carbonate-silicate cycle (carbon cycle) to maintain clement climates on Earth and potentially similar temperate exoplanets. We implement thermodynamics to determine weathering rates as a function of surface lithology (rock type). These rates provide upper limits that allow estimating the maximum rate of weathering in regulating climate. This modeling shows that the weathering of mineral assemblages in a given rock, rather than individual minerals, is crucial to determine weathering rates at planetary surfaces. By implementing a fluid-transport controlled approach, we further mimic chemical kinetics and thermodynamics to determine weathering rates for three types of rocks inspired by the lithologies of Earth's continental and oceanic crust, and its upper mantle. We find that thermodynamic weathering rates of a continental crust-like lithology are about one to two orders of magnitude lower than those of a lithology characteristic of the oceanic crust. We show that when the CO$_2$ partial pressure decreases or surface temperature increases, thermodynamics rather than kinetics exerts a strong control on weathering. The kinetically- and thermodynamically-limited regimes of weathering depend on lithology, whereas, the supply-limited weathering is independent of lithology. Our results imply that the temperature-sensitivity of thermodynamically-limited silicate weathering may instigate a positive feedback to the carbon cycle, in which the weathering rate decreases as the surface temperature increases., Comment: Accepted for publication in The Planetary Science Journal

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., Comment: 22 pages, 9 figures, 2 tables, accepted for publication in the Astrophysical Journal

THOR is the first open-source general circulation model (GCM) developed from scratch to study the atmospheres and climates of exoplanets, free from Earth- or Solar System-centric tunings. It solves the general non-hydrostatic Euler equations (instead of the primitive equations) on a sphere using the icosahedral grid. In the current study, we report major upgrades to THOR, building upon the work of Mendonça et al. (2016). First, while the Horizontally Explicit Vertically Implicit (HEVI) integration scheme is the same as that described in Mendonça et al. (2016), we provide a clearer description of the scheme and improved its implementation in the code. The differences in implementation between the hydrostatic shallow (HSS), quasi-hydrostatic deep (QHD) and non-hydrostatic deep (NHD) treatments are fully detailed. Second, standard physics modules are added: two-stream, double-gray radiative transfer and dry convective adjustment. Third, THOR is tested on additional benchmarks: tidally-locked Earth, deep hot Jupiter, acoustic wave, and gravity wave. Fourth, we report that differences between the hydrostatic and non-hydrostatic simulations are negligible in the Earth case, but pronounced in the hot Jupiter case. Finally, the effects of the so-called "sponge layer", a form of drag implemented in most GCMs to provide numerical stability, are examined. Overall, these upgrades have improved the flexibility, user-friendliness, and stability of THOR.

The atmospheres of temperate planets may be regulated by geochemical cycles. Silicate weathering provides essential negative feedback to the carbonate-silicate cycle (carbon cycle) to maintain temperate climates on Earth and possibly on Earth-sized temperate exoplanets. The intensity of weathering is normally attributed to the kinetics of weathering reactions of individual minerals. The implementation of a transport-controlled weathering model shows that when the CO2 volume mixing ratio decreases or surface temperature increases, equilibrium thermodynamics rather than kinetics exerts a strong control on weathering. Modeling the weathering of all minerals in a given rock instead of individual minerals is crucial. The transition between kinetically- and thermodynamically-limited regimes of weathering is strongly sensitive to rock lithology. Application of this model to Earth suggests that global mean continental granite and seafloor basalt weathering rates are likely limited by the supply of fresh rocks, yet regional weathering rates can be influenced by both kinetics and thermodynamics. Consideration of total dissolved inorganic carbon as a proxy for weathering results in another CO2 drawdown regime: CO2 dissolution, where aqueous bicarbonate and carbonate ions produced by rock weathering are lower in concentration than aqueous CO2. Upper limits to weathering as a function of lithology are provided to calculate the maximum impact of weathering on the carbon cycle. The temperature-sensitivity of the thermodynamically-limited silicate weathering provides a potential positive feedback to the carbon cycle which may shift the inner edge of the habitable zone.

THOR is the first open-source general circulation model (GCM) developed from scratch to study the atmospheres and climates of exoplanets, free from Earth- or Solar System-centric tunings. It solves the general non-hydrostatic Euler equations (instead of the primitive equations) on a sphere using the icosahedral grid. In the current study, we report major upgrades to THOR, building upon the work of Mendon\c{c}a et al. (2016). First, while the Horizontally Explicit Vertically Implicit (HEVI) integration scheme is the same as that described in Mendon\c{c}a et al. (2016), we provide a clearer description of the scheme and improved its implementation in the code. The differences in implementation between the hydrostatic shallow (HSS), quasi-hydrostatic deep (QHD) and non-hydrostatic deep (NHD) treatments are fully detailed. Second, standard physics modules are added: two-stream, double-gray radiative transfer and dry convective adjustment. Third, THOR is tested on additional benchmarks: tidally-locked Earth, deep hot Jupiter, acoustic wave, and gravity wave. Fourth, we report that differences between the hydrostatic and non-hydrostatic simulations are negligible in the Earth case, but pronounced in the hot Jupiter case. Finally, the effects of the so-called "sponge layer", a form of drag implemented in most GCMs to provide numerical stability, are examined. Overall, these upgrades have improved the flexibility, user-friendliness, and stability of THOR., Comment: 57 pages, 31 figures, revised, accepted for publication in ApJS

We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled simultaneously to simulate the evolution of terrestrial planets, gaseous planets, and stars. The code is validated by reproducing a selection of observations and past results. VPLanet is written in C and designed so that the user can choose the physics modules to apply to an individual object at runtime without recompiling, i.e., a single executable can simulate the diverse phenomena that are relevant to a wide range of planetary and stellar systems. This feature is enabled by matrices and vectors of function pointers that are dynamically allocated and populated based on user input. The speed and modularity of VPLanet enables large parameter sweeps and the versatility to add/remove physical phenomena to assess their importance. VPLanet is publicly available from a repository that contains extensive documentation, numerous examples, Python scripts for plotting and data management, and infrastructure for community input and future development., Comment: 75 pages, 34 figures, 10 tables, accepted to the Proceedings of the Astronomical Society of the Pacific. Source code, documentation, and examples available at https://github.com/VirtualPlanetaryLaboratory/vplanet

We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled simultaneously to simulate the evolution of terrestrial planets, gaseous planets, and stars. The code is validated by reproducing a selection of observations and past results. VPLanet is written in C and designed so that the user can choose the physics modules to apply to an individual object at runtime without recompiling, i.e., a single executable can simulate the diverse phenomena that are relevant to a wide range of planetary and stellar systems. This feature is enabled by matrices and vectors of function pointers that are dynamically allocated and populated based on user input. The speed and modularity of VPLanet enables large parameter sweeps and the versatility to add/remove physical phenomena to assess their importance. VPLanet is publicly available from a repository that contains extensive documentation, numerous examples, Python scripts for plotting and data management, and infrastructure for community input and future development., Comment: 75 pages, 34 figures, 10 tables, accepted to the Proceedings of the Astronomical Society of the Pacific. Source code, documentation, and examples available at https://github.com/VirtualPlanetaryLaboratory/vplanet

THOR is the first open-source general circulation model (GCM) developed from scratch to study the atmospheres and climates of exoplanets, free from Earth- or Solar System-centric tunings. It solves the general non-hydrostatic Euler equations (instead of the primitive equations) on a sphere using the icosahedral grid. In the current study, we report major upgrades to THOR, building upon the work of Mendon\c{c}a et al. (2016). First, while the Horizontally Explicit Vertically Implicit (HEVI) integration scheme is the same as that described in Mendon\c{c}a et al. (2016), we provide a clearer description of the scheme and improved its implementation in the code. The differences in implementation between the hydrostatic shallow (HSS), quasi-hydrostatic deep (QHD) and non-hydrostatic deep (NHD) treatments are fully detailed. Second, standard physics modules are added: two-stream, double-gray radiative transfer and dry convective adjustment. Third, THOR is tested on additional benchmarks: tidally-locked Earth, deep hot Jupiter, acoustic wave, and gravity wave. Fourth, we report that differences between the hydrostatic and non-hydrostatic simulations are negligible in the Earth case, but pronounced in the hot Jupiter case. Finally, the effects of the so-called "sponge layer", a form of drag implemented in most GCMs to provide numerical stability, are examined. Overall, these upgrades have improved the flexibility, user-friendliness, and stability of THOR., Comment: 57 pages, 31 figures, revised, accepted for publication in ApJS

Exoplanet hunting efforts have revealed the prevalence of exotic worlds with diverse properties, including Earth-sized bodies, which has fueled our endeavor to search for life beyond the Solar System. Accumulating experiences in astrophysical, chemical, and climatological characterization of uninhabitable planets are paving the way to characterization of potentially habitable planets. In this paper, we review our possibilities and limitations in characterizing temperate terrestrial planets with future observational capabilities through the 2030s and beyond, as a basis of a broad range of discussions on how to advance "astrobiology" with exoplanets. We discuss the observability of not only the proposed biosignature candidates themselves but also of more general planetary properties that provide circumstantial evidence, since the evaluation of any biosignature candidate relies on its context. Characterization of temperate Earth-sized planets in the coming years will focus on those around nearby late-type stars. The James Webb Space Telescope (JWST) and later 30-meter-class ground-based telescopes will empower their chemical investigations. Spectroscopic studies of potentially habitable planets around solar-type stars will likely require a designated spacecraft mission for direct imaging, leveraging technologies that are already being developed and tested as part of the Wide Field InfraRed Survey Telescope (WFIRST) mission. Successful initial characterization of a few nearby targets will be an important touchstone toward a more detailed scrutiny and a larger survey that are envisioned beyond 2030. The broad outlook this paper presents may help develop new observational techniques to detect relevant features as well as frameworks to diagnose planets based on the observables. Key Words: Exoplanets-Biosignatures-Characterization-Planetary atmospheres-Planetary surfaces. Astrobiology 18, 739-778.

The chemical composition of an exoplanet is a key ingredient in constraining its formation history. Iron is the most abundant transition metal, but has never been directly detected in an exoplanet due to its highly refractory nature. KELT-9b (HD 195689b) is the archetype of the class of ultra-hot Jupiters that straddle the transition between stars and gas-giant exoplanets and serve as distinctive laboratories for studying atmospheric chemistry, because of its high equilibrium temperature of 4050 +/- 180 K. These properties imply that its atmosphere is a tightly constrained chemical system that is expected to be nearly in chemical equilibrium and cloud-free. It was previously predicted that the spectral lines of iron will be detectable in the visible range of wavelengths. At these high temperatures, iron and several other transition metals are not sequestered in molecules or cloud particles and exist solely in their atomic forms. Here, we report the direct detection of atomic neutral and singly-ionized iron (Fe and Fe+), and singly-ionized titanium (Ti+) in KELT-9b via the cross-correlation technique applied to high-resolution spectra obtained during the primary transit of the exoplanet., Comment: 7 pages, 4 figures

We derive analytic, closed form, numerically stable solutions for the total flux received from a spherical planet, moon or star during an occultation if the specific intensity map of the body is expressed as a sum of spherical harmonics. Our expressions are valid to arbitrary degree and may be computed recursively for speed. The formalism we develop here applies to the computation of stellar transit light curves, planetary secondary eclipse light curves, and planet-planet/planet-moon occultation light curves, as well as thermal (rotational) phase curves. In this paper we also introduce STARRY, an open-source package written in C++ and wrapped in Python that computes these light curves. The algorithm in STARRY is six orders of magnitude faster than direct numerical integration and several orders of magnitude more precise. STARRY also computes analytic derivatives of the light curves with respect to all input parameters for use in gradient-based optimization and inference, such as Hamiltonian Monte Carlo (HMC), allowing users to quickly and efficiently fit observed light curves to infer properties of a celestial body's surface map., Comment: 55 pages, 20 figures. Accepted to the Astronomical Journal. Check out the code at https://github.com/rodluger/starry

The atmospheric temperatures of the ultra-hot Jupiter KELT-9b straddle the transition between gas giants and stars, and therefore between two traditionally distinct regimes of atmospheric chemistry. Previous theoretical studies assume the atmosphere of KELT-9b to be in chemical equilibrium. Despite the high ultraviolet flux from KELT-9, we show using photochemical kinetics calculations that the observable atmosphere of KELT-9b is predicted to be close to chemical equilibrium, which greatly simplifies any theoretical interpretation of its spectra. It also makes the atmosphere of KELT-9b, which is expected to be cloudfree, a tightly constrained chemical system that lends itself to a clean set of theoretical predictions. Due to the lower pressures probed in transmission (compared to emission) spectroscopy, we predict the abundance of water to vary by several orders of magnitude across the atmospheric limb depending on temperature, which makes water a sensitive thermometer. Carbon monoxide is predicted to be the dominant molecule under a wide range of scenarios, rendering it a robust diagnostic of the metallicity when analyzed in tandem with water. All of the other usual suspects (acetylene, ammonia, carbon dioxide, hydrogen cyanide, methane) are predicted to be subdominant at solar metallicity, while atomic oxygen, iron and magnesium are predicted to have relative abundances as high as 1 part in 10,000. Neutral atomic iron is predicted to be seen through a forest of optical and near-infrared lines, which makes KELT-9b suitable for high-resolution ground-based spectroscopy with HARPS-N or CARMENES. We summarize future observational prospects of characterizing the atmosphere of KELT-9b., Comment: Accepted by ApJ. 9 pages, 6 figures. Corrected minor errors in Figures 1a and 1b (some line styles were switched by accident), text and conclusions unchanged, these minor changes will be updated in final ApJ proof

Using an energy balance model with ice sheets, we examine the climate response of an Earth-like planet orbiting a G dwarf star and experiencing large orbital and obliquity variations. We find that ice caps couple strongly to the orbital forcing, leading to extreme ice ages. In contrast with previous studies, we find that such exo-Milankovitch cycles tend to impair habitability by inducing snowball states within the habitable zone. The large amplitude changes in obliquity and eccentricity cause the ice edge, the lowest latitude extent of the ice caps, to become unstable and grow to the equator. We apply an analytical theory of the ice edge latitude to show that obliquity is the primary driver of the instability. The thermal inertia of the ice sheets and the spectral energy distribution of the G dwarf star increase the sensitivity of the model to triggering runaway glaciation. Finally, we apply a machine learning algorithm to demonstrate how this technique can be used to extend the power of climate models. This work illustrates the importance of orbital evolution for habitability in dynamically rich planetary systems. We emphasize that as potentially habitable planets are discovered around G dwarfs, we need to consider orbital dynamics., Comment: 37 pages, 26 figures, accepted at the Astronomical Journal

For the first time in human history, we will soon be able to apply the scientific method to the question "Are We Alone?" The rapid advance of exoplanet discovery, planetary systems science, and telescope technology will soon allow scientists to search for life beyond our Solar System through direct observation of extrasolar planets. This endeavor will occur alongside searches for habitable environments and signs of life within our Solar System. While the searches are thematically related and will inform each other, they will require separate observational techniques. The search for life on exoplanets holds potential through the great diversity of worlds to be explored beyond our Solar System. However, there are also unique challenges related to the relatively limited data this search will obtain on any individual world. This white paper reviews the scientific community's ability to use data from future telescopes to search for life on exoplanets. This material summarizes products from the Exoplanet Biosignatures Workshop Without Walls (EBWWW). The EBWWW was constituted by a series of online and in person activities, with participation from the international exoplanet and astrobiology communities, to assess state of the science and future research needs for the remote detection of life on planets outside our Solar System., Comment: This is a white paper that was submitted to the National Academies of Sciences Study: Astrobiology Science Strategy for the Search for Life in the Universe

We describe how environmental context can help determine whether oxygen (O2) detected in extrasolar planetary observations is more likely to have a biological source. Here we provide an in-depth, interdisciplinary example of O2 biosignature identification and observation, which serves as the prototype for the development of a general framework for biosignature assessment. Photosynthetically generated O2 is a potentially strong biosignature, and at high abundance, it was originally thought to be an unambiguous indicator for life. However, as a biosignature, O2 faces two major challenges: (1) it was only present at high abundance for a relatively short period of Earth's history and (2) we now know of several potential planetary mechanisms that can generate abundant O2 without life being present. Consequently, our ability to interpret both the presence and absence of O2 in an exoplanetary spectrum relies on understanding the environmental context. Here we examine the coevolution of life with the early Earth's environment to identify how the interplay of sources and sinks may have suppressed O2 release into the atmosphere for several billion years, producing a false negative for biologically generated O2. These studies suggest that planetary characteristics that may enhance false negatives should be considered when selecting targets for biosignature searches. We review the most recent knowledge of false positives for O2, planetary processes that may generate abundant atmospheric O2 without a biosphere. We provide examples of how future photometric, spectroscopic, and time-dependent observations of O2 and other aspects of the planetary environment can be used to rule out false positives and thereby increase our confidence that any observed O2 is indeed a biosignature. These insights will guide and inform the development of future exoplanet characterization missions. Key Words: Biosignatures-Oxygenic photosynthesis-Exoplanets-Planetary atmospheres. Astrobiology 18, 630-66

For the first time in human history, we will soon be able to apply the scientific method to the question "Are We Alone?" The rapid advance of exoplanet discovery, planetary systems science, and telescope technology will soon allow scientists to search for life beyond our Solar System through direct observation of extrasolar planets. This endeavor will occur alongside searches for habitable environments and signs of life within our Solar System. While the searches are thematically related and will inform each other, they will require separate observational techniques. The search for life on exoplanets holds potential through the great diversity of worlds to be explored beyond our Solar System. However, there are also unique challenges related to the relatively limited data this search will obtain on any individual world. This white paper reviews the scientific community's ability to use data from future telescopes to search for life on exoplanets. This material summarizes products from the Exoplanet Biosignatures Workshop Without Walls (EBWWW). The EBWWW was constituted by a series of online and in person activities, with participation from the international exoplanet and astrobiology communities, to assess state of the science and future research needs for the remote detection of life on planets outside our Solar System.

We describe how environmental context can help determine whether oxygen (O2) detected in extrasolar planetary observations is more likely to have a biological source. Here we provide an in-depth, interdisciplinary example of O2 biosignature identification and observation, which serves as the prototype for the development of a general framework for biosignature assessment. Photosynthetically generated O2 is a potentially strong biosignature, and at high abundance, it was originally thought to be an unambiguous indicator for life. However, as a biosignature, O2 faces two major challenges: (1) it was only present at high abundance for a relatively short period of Earths history and (2) we now know of several potential planetary mechanisms that can generate abundant O2 without life being present. Consequently, our ability to interpret both the presence and absence of O2 in an exoplanetary spectrum relies on understanding the environmental context. Here we examine the coevolution of life with the early Earths environment to identify how the interplay of sources and sinks may have suppressed O2 release into the atmosphere for several billion years, producing a false negative for biologically generated O2. These studies suggest that planetary characteristics that may enhance false negatives should be considered when selecting targets for biosignature searches. We review the most recent knowledge of false positives for O2, planetary processes that may generate abundant atmospheric O2 without a biosphere. We provide examples of how future photometric, spectroscopic, and time-dependent observations of O2 and other aspects of the planetary environment can be used to rule out false positives and thereby increase our confidence that any observed O2 is indeed a biosignature. These insights will guide and inform the development of future exoplanet characterization missions. Key Words: Biosignatures-Oxygenic photosynthesis-Exoplanets-Planetary atmospheres. Astrobiology 18, 630-662.

We describe how environmental context can help determine whether oxygen (O2) detected in extrasolar planetary observations is more likely to have a biological source. Here we provide an in-depth, interdisciplinary example of O2 biosignature identification and observation, which serves as the prototype for the development of a general framework for biosignature assessment. Photosynthetically generated O2 is a potentially strong biosignature, and at high abundance, it was originally thought to be an unambiguous indicator for life. However, as a biosignature, O2 faces two major challenges: (1) it was only present at high abundance for a relatively short period of Earth's history and (2) we now know of several potential planetary mechanisms that can generate abundant O2 without life being present. Consequently, our ability to interpret both the presence and absence of O2 in an exoplanetary spectrum relies on understanding the environmental context. Here we examine the coevolution of life with the early Earth's environment to identify how the interplay of sources and sinks may have suppressed O2 release into the atmosphere for several billion years, producing a false negative for biologically generated O2. These studies suggest that planetary characteristics that may enhance false negatives should be considered when selecting targets for biosignature searches. We review the most recent knowledge of false positives for O2, planetary processes that may generate abundant atmospheric O2 without a biosphere. We provide examples of how future photometric, spectroscopic, and time-dependent observations of O2 and other aspects of the planetary environment can be used to rule out false positives and thereby increase our confidence that any observed O2 is indeed a biosignature. These insights will guide and inform the development of future exoplanet characterization missions. Key Words: Biosignatures-Oxygenic photosynthesis-Exoplanets-Planetary atmospheres. Astrobiology 18, 630-66

For the first time in human history, we will soon be able to apply the scientific method to the question "Are We Alone?" The rapid advance of exoplanet discovery, planetary systems science, and telescope technology will soon allow scientists to search for life beyond our Solar System through direct observation of extrasolar planets. This endeavor will occur alongside searches for habitable environments and signs of life within our Solar System. While the searches are thematically related and will inform each other, they will require separate observational techniques. The search for life on exoplanets holds potential through the great diversity of worlds to be explored beyond our Solar System. However, there are also unique challenges related to the relatively limited data this search will obtain on any individual world. This white paper reviews the scientific community's ability to use data from future telescopes to search for life on exoplanets. This material summarizes products from the Exoplanet Biosignatures Workshop Without Walls (EBWWW). The EBWWW was constituted by a series of online and in person activities, with participation from the international exoplanet and astrobiology communities, to assess state of the science and future research needs for the remote detection of life on planets outside our Solar System.

We describe how environmental context can help determine whether oxygen (O2) detected in extrasolar planetary observations is more likely to have a biological source. Here we provide an in-depth, interdisciplinary example of O2 biosignature identification and observation, which serves as the prototype for the development of a general framework for biosignature assessment. Photosynthetically generated O2 is a potentially strong biosignature, and at high abundance, it was originally thought to be an unambiguous indicator for life. However, as a biosignature, O2 faces two major challenges: (1) it was only present at high abundance for a relatively short period of Earth's history and (2) we now know of several potential planetary mechanisms that can generate abundant O2 without life being present. Consequently, our ability to interpret both the presence and absence of O2 in an exoplanetary spectrum relies on understanding the environmental context. Here we examine the coevolution of life with the early Earth's environment to identify how the interplay of sources and sinks may have suppressed O2 release into the atmosphere for several billion years, producing a false negative for biologically generated O2. These studies suggest that planetary characteristics that may enhance false negatives should be considered when selecting targets for biosignature searches. We review the most recent knowledge of false positives for O2, planetary processes that may generate abundant atmospheric O2 without a biosphere. We provide examples of how future photometric, spectroscopic, and time-dependent observations of O2 and other aspects of the planetary environment can be used to rule out false positives and thereby increase our confidence that any observed O2 is indeed a biosignature. These insights will guide and inform the development of future exoplanet characterization missions. Key Words: Biosignatures-Oxygenic photosynthesis-Exoplanets-Planetary atmospheres. Astrobiology 18, 630-66

Exoplanet hunting efforts have revealed the prevalence of exotic worlds with diverse properties, including Earth-sized bodies, which has fueled our endeavor to search for life beyond the Solar System. Accumulating experiences in astrophysical, chemical, and climatological characterization of uninhabitable planets are paving the way to characterization of potentially habitable planets. In this paper, we review our possibilities and limitations in characterizing temperate terrestrial planets with future observational capabilities through the 2030s and beyond, as a basis of a broad range of discussions on how to advance "astrobiology" with exoplanets. We discuss the observability of not only the proposed biosignature candidates themselves but also of more general planetary properties that provide circumstantial evidence, since the evaluation of any biosignature candidate relies on its context. Characterization of temperate Earth-sized planets in the coming years will focus on those around nearby late-type stars. The James Webb Space Telescope (JWST) and later 30-meter-class ground-based telescopes will empower their chemical investigations. Spectroscopic studies of potentially habitable planets around solar-type stars will likely require a designated spacecraft mission for direct imaging, leveraging technologies that are already being developed and tested as part of the Wide Field InfraRed Survey Telescope (WFIRST) mission. Successful initial characterization of a few nearby targets will be an important touchstone toward a more detailed scrutiny and a larger survey that are envisioned beyond 2030. The broad outlook this paper presents may help develop new observational techniques to detect relevant features as well as frameworks to diagnose planets based on the observables. Key Words: Exoplanets-Biosignatures-Characterization-Planetary atmospheres-Planetary surfaces. Astrobiology 18, 739-778.

We derive analytic, closed form, numerically stable solutions for the total flux received from a spherical planet, moon or star during an occultation if the specific intensity map of the body is expressed as a sum of spherical harmonics. Our expressions are valid to arbitrary degree and may be computed recursively for speed. The formalism we develop here applies to the computation of stellar transit light curves, planetary secondary eclipse light curves, and planet-planet/planet-moon occultation light curves, as well as thermal (rotational) phase curves. In this paper we also introduce STARRY, an open-source package written in C++ and wrapped in Python that computes these light curves. The algorithm in STARRY is six orders of magnitude faster than direct numerical integration and several orders of magnitude more precise. STARRY also computes analytic derivatives of the light curves with respect to all input parameters for use in gradient-based optimization and inference, such as Hamiltonian Monte Carlo (HMC), allowing users to quickly and efficiently fit observed light curves to infer properties of a celestial body's surface map., Comment: 55 pages, 20 figures. Accepted to the Astronomical Journal. Check out the code at https://github.com/rodluger/starry

The chemical composition of an exoplanet is a key ingredient in constraining its formation history. Iron is the most abundant transition metal, but has never been directly detected in an exoplanet due to its highly refractory nature. KELT-9b (HD 195689b) is the archetype of the class of ultra-hot Jupiters that straddle the transition between stars and gas-giant exoplanets and serve as distinctive laboratories for studying atmospheric chemistry, because of its high equilibrium temperature of 4050 +/- 180 K. These properties imply that its atmosphere is a tightly constrained chemical system that is expected to be nearly in chemical equilibrium and cloud-free. It was previously predicted that the spectral lines of iron will be detectable in the visible range of wavelengths. At these high temperatures, iron and several other transition metals are not sequestered in molecules or cloud particles and exist solely in their atomic forms. Here, we report the direct detection of atomic neutral and singly-ionized iron (Fe and Fe+), and singly-ionized titanium (Ti+) in KELT-9b via the cross-correlation technique applied to high-resolution spectra obtained during the primary transit of the exoplanet., Comment: 7 pages, 4 figures