Showing 1–42 of 42 results for author: Zahnle, K

Noble Gas Planetology and the Xenon Clouds of Uranus

Authors: Kevin J Zahnle

Abstract: Noble gases provide tracers of cosmic provenance that are accessible to a future Uranus Atmospheric Probe. Argon and krypton are expected to be well-mixed on Uranus with respect to H$_2$ and He, although condensation at the winter pole may be possible. The Ar/H$_2$ and Ar/Kr ratios address whether the materials accreted by Uranus resembled the extremely cold materials accreted by Jupiter's atmosph… ▽ More Noble gases provide tracers of cosmic provenance that are accessible to a future Uranus Atmospheric Probe. Argon and krypton are expected to be well-mixed on Uranus with respect to H$_2$ and He, although condensation at the winter pole may be possible. The Ar/H$_2$ and Ar/Kr ratios address whether the materials accreted by Uranus resembled the extremely cold materials accreted by Jupiter's atmosphere, or whether they were warmer like comet 67P/Churyumov-Gerasimenko, or whether Uranus is like neither. Xenon condenses as an ice, probably on methane ice, in Uranus's upper troposphere. Condensation may complicate interpretation of Xe/H$_2$, but it also presents an opportunity to collect concentrated xenon samples suitable for measuring isotopes. Solar System Xe tracks three distinct nucleosynthetic xenon reservoirs, one evident in the Sun and in chondritic meteorites, a second in refractory presolar grains, and a third evident in comet 67P/C-G and in Earth's air. The first and third reservoirs appear to have been captured from different clouds of gas. The two gases do not appear to have been well-mixed; moreover, the high $^{129}$Xe/$^{132}$Xe ratio in 67P/C-G implies that the gas was captured before the initial nucleosynthetic complement of $^{129}$I (15.7 Myr half-life) had decayed. Xenon's isotopic peculiarities, if seen in Uranus, could usefully upset our understanding of planetary origins. Krypton's isotopic anomalies are more subtle and may prove hard to measure. There is a slight chance that neon and helium fractionations can be used to constrain how Uranus acquired its nebular envelope. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 26 pages, 8 figures

JWST observations of K2-18b can be explained by a gas-rich mini-Neptune with no habitable surface

Authors: Nicholas F. Wogan, Natasha E. Batalha, Kevin Zahnle, Joshua Krissansen-Totton, Shang-Min Tsai, Renyu Hu

Abstract: JWST recently measured the transmission spectrum of K2-18b, a habitable-zone sub-Neptune exoplanet, detecting CH$_4$ and CO$_2$ in its atmosphere. The discovery paper argued the data are best explained by a habitable "Hycean" world, consisting of a relatively thin H$_2$-dominated atmosphere overlying a liquid water ocean. Here, we use photochemical and climate models to simulate K2-18b as both a H… ▽ More JWST recently measured the transmission spectrum of K2-18b, a habitable-zone sub-Neptune exoplanet, detecting CH$_4$ and CO$_2$ in its atmosphere. The discovery paper argued the data are best explained by a habitable "Hycean" world, consisting of a relatively thin H$_2$-dominated atmosphere overlying a liquid water ocean. Here, we use photochemical and climate models to simulate K2-18b as both a Hycean planet and a gas-rich mini-Neptune with no defined surface. We find that a lifeless Hycean world is hard to reconcile with the JWST observations because photochemistry only supports $< 1$ part-per-million CH$_4$ in such an atmosphere while the data suggest about $\sim 1\%$ of the gas is present. Sustaining %-level CH$_4$ on a Hycean K2-18b may require the presence of a methane-producing biosphere, similar to microbial life on Earth $\sim 3$ billion years ago. On the other hand, we predict that a gas-rich mini-Neptune with $100 \times$ solar metallicity should have 4% CH$_4$ and nearly 0.1% CO$_2$, which are compatible with the JWST data. The CH$_4$ and CO$_2$ are produced thermochemically in the deep atmosphere and mixed upward to the low pressures sensitive to transmission spectroscopy. The model predicts H$_2$O, NH$_3$ and CO abundances broadly consistent with the non-detections. Given the additional obstacles to maintaining a stable temperate climate on Hycean worlds due to H$_2$ escape and potential supercriticality at depth, we favor the mini-Neptune interpretation because of its relative simplicity and because it does not need a biosphere or other unknown source of methane to explain the data. △ Less

Submitted 3 February, 2024; v1 submitted 19 January, 2024; originally announced January 2024.

Comments: Accepted for publication at ApJL

Impact Rates in the Outer Solar System

Authors: David Nesvorny, Luke Dones, Mario De Pra, Maria Womack, Kevin J. Zahnle

Abstract: Previous studies of cometary impacts in the outer Solar System used the spatial distribution of ecliptic comets (ECs) from dynamical models that assumed ECs began on low-inclination orbits (<5 deg) in the Kuiper belt. In reality, the source population of ECs - the trans-Neptunian scattered disk - has orbital inclinations reaching up to ~30 deg. In Nesvorny et al. (2017), we developed a new dynamic… ▽ More Previous studies of cometary impacts in the outer Solar System used the spatial distribution of ecliptic comets (ECs) from dynamical models that assumed ECs began on low-inclination orbits (<5 deg) in the Kuiper belt. In reality, the source population of ECs - the trans-Neptunian scattered disk - has orbital inclinations reaching up to ~30 deg. In Nesvorny et al. (2017), we developed a new dynamical model of ECs by following comets as they evolved from the scattered disk to the inner Solar System. The model was absolutely calibrated from the population of Centaurs and active ECs. Here we use our EC model to determine the steady-state impact flux of cometary/Centaur impactors on Jupiter, Saturn, Uranus, and their moons. Relative to previous work (Zahnle et al. 2003), we find slightly higher impact probabilities on the outer moons and lower impact probabilities on the inner moons. The impact probabilities are smaller when comet disruption is accounted for. The results provide a modern framework for the interpretation of the cratering record in the outer Solar System. △ Less

Submitted 28 July, 2023; originally announced July 2023.

Comments: PSJ, in press

Origin of Life Molecules in the Atmosphere After Big Impacts on the Early Earth

Authors: Nicholas F. Wogan, David C. Catling, Kevin J. Zahnle, Roxana Lupu

Abstract: The origin of life on Earth would benefit from a prebiotic atmosphere that produced nitriles, like HCN, which enable ribonucleotide synthesis. However, geochemical evidence suggests that Hadean air was relatively oxidizing with negligible photochemical production of prebiotic molecules. These paradoxes are resolved by iron-rich asteroid impacts that transiently reduced the entire atmosphere, allow… ▽ More The origin of life on Earth would benefit from a prebiotic atmosphere that produced nitriles, like HCN, which enable ribonucleotide synthesis. However, geochemical evidence suggests that Hadean air was relatively oxidizing with negligible photochemical production of prebiotic molecules. These paradoxes are resolved by iron-rich asteroid impacts that transiently reduced the entire atmosphere, allowing nitriles to form in subsequent photochemistry. Here, we investigate impact-generated reducing atmospheres using new time-dependent, coupled atmospheric chemistry and climate models, which account for gas-phase reactions and surface-catalysis. The resulting H$_2$-, CH$_4$- and NH$_3$-rich atmospheres persist for millions of years, until hydrogen escapes to space. HCN and HCCCN production and rainout to the surface can reach $10^9$ molecules cm$^{-2}$ s$^{-1}$ in hazy atmospheres with a mole ratio of $\mathrm{CH_4} / \mathrm{CO_2} > 0.1$. Smaller $\mathrm{CH_4} / \mathrm{CO_2}$ ratios produce HCN rainout rates $< 10^5$ molecules cm$^{-2}$ s$^{-1}$, and negligible HCCCN. The minimum impactor mass that creates atmospheric $\mathrm{CH_4} / \mathrm{CO_2} > 0.1$ is $4 \times 10^{20}$ to $5 \times 10^{21}$ kg (570 to 1330 km diameter), depending on how efficiently iron reacts with a steam atmosphere, the extent of atmospheric equilibration with an impact-induced melt pond, and the surface area of nickel that catalyzes CH$_4$ production. Alternatively, if steam permeates and deeply oxidizes crust, impactors $\sim 10^{20}$ kg could be effective. Atmospheres with copious nitriles have $> 360$ K surface temperatures, perhaps posing a challenge for RNA longevity, although cloud albedo can produce cooler climates. Regardless, post-impact cyanide can be stockpiled and used in prebiotic schemes after hydrogen has escaped to space. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: Accepted by Planetary Science Journal on 7/4/23

Revealing the Mysteries of Venus: The DAVINCI Mission

Authors: James B. Garvin, Stephanie A. Getty, Giada N. Arney, Natasha M. Johnson, Erika Kohler, Kenneth O. Schwer, Michael Sekerak, Arlin Bartels, Richard S. Saylor, Vincent E. Elliott, Colby S. Goodloe, Matthew B. Garrison, Valeria Cottini, Noam Izenberg, Ralph Lorenz, Charles A. Malespin, Michael Ravine, Christopher R. Webster, David H. Atkinson, Shahid Aslam, Sushil Atreya, Brent J. Bos, William B. Brinckerhoff, Bruce Campbell, David Crisp , et al. (14 additional authors not shown)

Abstract: The Deep Atmosphere Venus Investigation of Noble gases, Chemistry, and Imaging (DAVINCI) mission described herein has been selected for flight to Venus as part of the NASA Discovery Program. DAVINCI will be the first mission to Venus to incorporate science-driven flybys and an instrumented descent sphere into a unified architecture. The anticipated scientific outcome will be a new understanding of… ▽ More The Deep Atmosphere Venus Investigation of Noble gases, Chemistry, and Imaging (DAVINCI) mission described herein has been selected for flight to Venus as part of the NASA Discovery Program. DAVINCI will be the first mission to Venus to incorporate science-driven flybys and an instrumented descent sphere into a unified architecture. The anticipated scientific outcome will be a new understanding of the atmosphere, surface, and evolutionary path of Venus as a possibly once-habitable planet and analog to hot terrestrial exoplanets. The primary mission design for DAVINCI as selected features a preferred launch in summer/fall 2029, two flybys in 2030, and descent sphere atmospheric entry by the end of 2031. The in situ atmospheric descent phase subsequently delivers definitive chemical and isotopic composition of the Venus atmosphere during a cloud-top to surface transect above Alpha Regio. These in situ investigations of the atmosphere and near infrared descent imaging of the surface will complement remote flyby observations of the dynamic atmosphere, cloud deck, and surface near infrared emissivity. The overall mission yield will be at least 60 Gbits (compressed) new data about the atmosphere and near surface, as well as first unique characterization of the deep atmosphere environment and chemistry, including trace gases, key stable isotopes, oxygen fugacity, constraints on local rock compositions, and topography of a tessera. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 41 pages, 14 figures, accepted for publication in the Planetary Science Journal

Journal ref: 2022, PSJ, 3, 117

Anoxic Atmospheres on Mars Driven by Volcanism: Implications for Past Environments and Life

Authors: Steven F. Sholes, Megan L. Smith, Mark W. Claire, Kevin J. Zahnle, David C. Catling

Abstract: Mars today has no active volcanism and its atmosphere is oxidizing, dominated by the photochemistry of CO2 and H2O. Using a one-dimensional photochemical model, we consider whether plausible volcanic gas fluxes could have switched the redox-state of the past martian atmosphere to reducing conditions. In our model, the total quantity and proportions of volcanic gases depend on the water content, ou… ▽ More Mars today has no active volcanism and its atmosphere is oxidizing, dominated by the photochemistry of CO2 and H2O. Using a one-dimensional photochemical model, we consider whether plausible volcanic gas fluxes could have switched the redox-state of the past martian atmosphere to reducing conditions. In our model, the total quantity and proportions of volcanic gases depend on the water content, outgassing pressure, and oxygen fugacity of the source melt. We find that with reasonable melt parameters the past martian atmosphere (~3.5 Gyr to present) could have easily reached reducing and anoxic conditions with modest levels of volcanism, >0.14 km^3/yr, well within the range of prior estimates. Counter-intuitively we also find that more reducing melts with lower oxygen fugacity require greater amounts of volcanism to switch a paleo-atmosphere from oxidizing to reducing. The reason is that sulfur is more stable in such melts and lower absolute fluxes of sulfur-bearing gases more than compensate for increases in the proportions of H2 and CO. These results imply that ancient Mars should have experienced periods with anoxic and reducing atmospheres even through the mid-Amazonian whenever volcanic outgassing was sustained at sufficient levels. Reducing anoxic conditions are potentially conducive to the synthesis of prebiotic organic compounds, such as amino acids, and are therefore relevant to the possibility of life on Mars. Also, anoxic reducing conditions should have influenced the type of minerals that were formed on the surface or deposited from the atmosphere such as elemental polysulfur (S8) as a signature of past reducing atmospheres. Finally, our models allow us to estimate the amount of volcanically sourced atmospheric sulfate deposited over Mars' history, approximately 10^6 to 10^9 Tmol, with a spread depending on assumed outgassing rate history and magmatic source conditions. △ Less

Submitted 24 March, 2021; originally announced March 2021.

Comments: 6 Figures, 1 Table, 3 Appendices

Journal ref: 2017, Icarus 290 pgs. 46-62

Atmospheric Erosion by Giant Impacts onto Terrestrial Planets: A Scaling Law for any Speed, Angle, Mass, and Density

Authors: Jacob A. Kegerreis, Vincent R. Eke, David C. Catling, Richard J. Massey, Luis F. A. Teodoro, Kevin J. Zahnle

Abstract: We present a new scaling law to predict the loss of atmosphere from planetary collisions for any speed, angle, impactor mass, target mass, and body compositions, in the regime of giant impacts onto broadly terrestrial planets with relatively thin atmospheres. To this end, we examine the erosion caused by a wide range of impacts, using 3D smoothed particle hydrodynamics simulations with sufficientl… ▽ More We present a new scaling law to predict the loss of atmosphere from planetary collisions for any speed, angle, impactor mass, target mass, and body compositions, in the regime of giant impacts onto broadly terrestrial planets with relatively thin atmospheres. To this end, we examine the erosion caused by a wide range of impacts, using 3D smoothed particle hydrodynamics simulations with sufficiently high resolution to directly model the fate of low-mass atmospheres around 1% of the target's mass. Different collision scenarios lead to extremely different behaviours and consequences for the planets. In spite of this complexity, the fraction of lost atmosphere is fitted well by a power law. Scaling is independent of the system mass for a constant impactor mass ratio. Slow atmosphere-hosting impactors can also deliver a significant mass of atmosphere, but always accompanied by larger proportions of their mantle and core. Different Moon-forming impact hypotheses suggest that around 10 to 60% of a primordial atmosphere could have been removed directly, depending on the scenario. We find no evident departure from the scaling trends at the extremes of the parameters explored. The scaling law can be incorporated readily into models of planet formation. △ Less

Submitted 30 September, 2020; v1 submitted 8 July, 2020; originally announced July 2020.

Comments: Published in ApJL. 12 pages, 6 figures

Creation and Evolution of Impact-generated Reduced Atmospheres of Early Earth

Authors: Kevin Zahnle, Roxana Lupu, David Catling, Nick Wogan

Abstract: The origin of life on Earth seems to demand a highly reduced early atmosphere, rich in CH4, H2, and NH3, but geological evidence suggests that Earth's mantle has always been relatively oxidized and its emissions dominated by CO2 H2O, and N2. The paradox can be resolved by exploiting the reducing power inherent in the "late veneer," i.e., material accreted by Earth after the Moon-forming impact. Is… ▽ More The origin of life on Earth seems to demand a highly reduced early atmosphere, rich in CH4, H2, and NH3, but geological evidence suggests that Earth's mantle has always been relatively oxidized and its emissions dominated by CO2 H2O, and N2. The paradox can be resolved by exploiting the reducing power inherent in the "late veneer," i.e., material accreted by Earth after the Moon-forming impact. Isotopic evidence indicates that the late veneer consisted of extremely dry, highly reduced inner solar system materials, suggesting that Earth's oceans were already present when the late veneer came. The major primary product of reaction between the late veneer's iron and Earth's water was H2. Ocean vaporizing impacts generate high pressures and long cooling times that favor CH4 and NH3. Impacts too small to vaporize the oceans are much less productive of CH4 and NH3, unless (i) catalysts were available to speed their formation, or (ii) additional reducing power was extracted from pre-existing crustal or mantle materials. The transient H2-CH4 atmospheres evolve photochemically to generate nitrogenated hydrocarbons at rates determined by solar radiation and hydrogen escape, on timescales ranging up to tens of millions of years and with cumulative organic production ranging up to half a kilometer. Roughly one ocean of hydrogen escapes. The atmosphere after the methane's gone is typically H2 and CO rich, with eventual oxidation to CO2 rate-limited by water photolysis and hydrogen escape. △ Less

Submitted 25 February, 2020; v1 submitted 31 December, 2019; originally announced January 2020.

Comments: Submitted 31 Dec 2019 to "The Planetary Science Journal." Accepted 10 Feb 2020

Venus as a Laboratory for Exoplanetary Science

Authors: Stephen R. Kane, Giada Arney, David Crisp, Shawn Domagal-Goldman2, Lori S. Glaze, Colin Goldblatt, David Grinspoon, James W. Head, Adrian Lenardic, Cayman Unterborn, Michael J. Way, Kevin J. Zahnle

Abstract: 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… ▽ More 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. △ Less

Submitted 26 August, 2019; v1 submitted 7 August, 2019; originally announced August 2019.

Comments: 21 pages, 4 figures, accepted for publication in JGR Planets

O($^3P$)+CO$_2$ scattering cross sections at superthermal collision energies for planetary aeronomy

Authors: Marko Gacesa, Robert J. Lillis, Kevin J. Zahnle

Abstract: We report new elastic and inelastic cross sections for O($^3P$)+CO$_2$ scattering at collision energies from 0.03 to 5 eV, of major importance to O escape from Mars, Venus, and CO$_2$-rich atmospheres. The cross sections were calculated from first principles using three newly constructed ab-initio potential energy surfaces correlating to the lowest energy asymptote of the complex. The surfaces wer… ▽ More We report new elastic and inelastic cross sections for O($^3P$)+CO$_2$ scattering at collision energies from 0.03 to 5 eV, of major importance to O escape from Mars, Venus, and CO$_2$-rich atmospheres. The cross sections were calculated from first principles using three newly constructed ab-initio potential energy surfaces correlating to the lowest energy asymptote of the complex. The surfaces were restricted to a planar geometry with the CO$_2$ molecule assumed to be in linear configuration fixed at equilibrium. Quantum-mechanical coupled-channel formalism with a large basis set was used to compute state-to-state integral and differential cross sections for elastic and inelastic O($^3P$)+CO$_2$ scattering between all pairs of rotational states of CO$_2$ molecule. The elastic cross sections are 35\% lower at 0.5 eV and more than 50\% lower at 4+ eV than values commonly used in studies of processes in upper and middle planetary atmospheres of Mars, Earth, Venus, and CO$_2$-rich planets. Momentum transfer cross sections, of interest for energy transport, were found to be lower than predicted by mass-scaling. △ Less

Submitted 16 January, 2021; v1 submitted 26 June, 2019; originally announced June 2019.

Comments: 10 pages, 8 figures; published in MNRAS

MSC Class: 81U35; 85A20

Journal ref: Monthly Notices of the Royal Astronomical Society 491 (4), 5650-5659 (2020)

An optical transmission spectrum for the ultra-hot Jupiter WASP-121b measured with the Hubble Space Telescope

Authors: Thomas M. Evans, David K. Sing, Jayesh Goyal, Nikolay Nikolov, Mark S. Marley, Kevin Zahnle, Gregory W. Henry, Joanna K. Barstow, Munazza K. Alam, Jorge Sanz-Forcada, Tiffany Kataria, Nikole K. Lewis, Panayotis Lavvas, Gilda E. Ballester, Lotfi Ben-Jaffel, Sarah D. Blumenthal, Vincent Bourrier, Benjamin Drummond, Antonio Garcia Munoz, Mercedes Lopez-Morales, Pascal Tremblin, David Ehrenreich, Hannah R. Wakeford, Lars A. Buchhave, Alain Lecavelier des Etangs , et al. (2 additional authors not shown)

Abstract: We present an atmospheric transmission spectrum for the ultra-hot Jupiter WASP-121b, measured using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). Across the 0.47-1 micron wavelength range, the data imply an atmospheric opacity comparable to - and in some spectroscopic channels exceeding - that previously measured at near-infrared wavelengths (1.15-1.65 m… ▽ More We present an atmospheric transmission spectrum for the ultra-hot Jupiter WASP-121b, measured using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). Across the 0.47-1 micron wavelength range, the data imply an atmospheric opacity comparable to - and in some spectroscopic channels exceeding - that previously measured at near-infrared wavelengths (1.15-1.65 micron). Wavelength-dependent variations in the opacity rule out a gray cloud deck at a confidence level of 3.8-sigma and may instead be explained by VO spectral bands. We find a cloud-free model assuming chemical equilibrium for a temperature of 1500K and metal enrichment of 10-30x solar matches these data well. Using a free-chemistry retrieval analysis, we estimate a VO abundance of -6.6(-0.3,+0.2) dex. We find no evidence for TiO and place a 3-sigma upper limit of -7.9 dex on its abundance, suggesting TiO may have condensed from the gas phase at the day-night limb. The opacity rises steeply at the shortest wavelengths, increasing by approximately five pressure scale heights from 0.47 to 0.3 micron in wavelength. If this feature is caused by Rayleigh scattering due to uniformly-distributed aerosols, it would imply an unphysically high temperature of 6810+/-1530K. One alternative explanation for the short-wavelength rise is absorption due to SH (mercapto radical), which has been predicted as an important product of non-equilibrium chemistry in hot Jupiter atmospheres. Irrespective of the identity of the NUV absorber, it likely captures a significant amount of incident stellar radiation at low pressures, thus playing a significant role in the overall energy budget, thermal structure, and circulation of the atmosphere. △ Less

Submitted 25 October, 2018; originally announced October 2018.

Comments: Accepted for publication in The Astronomical Journal

Strange messenger: A new history of hydrogen on Earth, as told by Xenon

Authors: Kevin J. Zahnle, Marko Gacesa, David C. Catling

Abstract: Atmospheric xenon is strongly mass fractionated, the result of a process that apparently continued through the Archean and perhaps beyond. Previous models that explain Xe fractionation by hydrodynamic hydrogen escape cannot gracefully explain how Xe escaped when Ar and Kr did not, nor allow Xe to escape in the Archean. Here we show that Xe is the only noble gas that can escape as an ion in a photo… ▽ More Atmospheric xenon is strongly mass fractionated, the result of a process that apparently continued through the Archean and perhaps beyond. Previous models that explain Xe fractionation by hydrodynamic hydrogen escape cannot gracefully explain how Xe escaped when Ar and Kr did not, nor allow Xe to escape in the Archean. Here we show that Xe is the only noble gas that can escape as an ion in a photo-ionized hydrogen wind, possible in the absence of a geomagnetic field or along polar magnetic field lines that open into interplanetary space. To quantify the hypothesis we construct new 1-D models of hydrodynamic diffusion-limited hydrogen escape from highly-irradiated CO$_2$-H$_2$-H atmospheres. The models reveal three minimum requirements for Xe escape: solar EUV irradiation needs to exceed $10\times$ that of the modern Sun; the total hydrogen mixing ratio in the atmosphere needs to exceed 1% (equiv. to 0.5% CH$_4$); and transport amongst the ions in the lower ionosphere needs to lift the Xe ions to the base of the outflowing hydrogen corona. The long duration of Xe escape implies that, if a constant process, Earth lost the hydrogen from at least one ocean of water, roughly evenly split between the Hadean and the Archean. However, to account for both Xe's fractionation and also its depletion with respect to Kr and primordial $^{244}$Pu, Xe escape must have been limited to small apertures or short episodes, which suggests that Xe escape was restricted to polar windows by a geomagnetic field, or dominated by outbursts of high solar activity, or limited to transient episodes of abundant hydrogen, or a combination of these. Xenon escape stopped when the hydrogen (or methane) mixing ratio became too small, or EUV radiation from the aging Sun became too weak, or charge exchange between Xe$^+$ and O$_2$ rendered Xe neutral. △ Less

Submitted 18 September, 2018; originally announced September 2018.

Comments: 58 pages 18 figures 5 tables

Consequences of Giant Impacts on Early Uranus for Rotation, Internal Structure, Debris, and Atmospheric Erosion

Authors: J. A. Kegerreis, L. F. A. Teodoro, V. R. Eke, R. J. Massey, D. C. Catling, C. L. Fryer, D. G. Korycansky, M. S. Warren, K. J. Zahnle

Abstract: We perform a suite of smoothed particle hydrodynamics simulations to investigate in detail the results of a giant impact on the young Uranus. We study the internal structure, rotation rate, and atmospheric retention of the post-impact planet, as well as the composition of material ejected into orbit. Most of the material from the impactor's rocky core falls in to the core of the target. However, f… ▽ More We perform a suite of smoothed particle hydrodynamics simulations to investigate in detail the results of a giant impact on the young Uranus. We study the internal structure, rotation rate, and atmospheric retention of the post-impact planet, as well as the composition of material ejected into orbit. Most of the material from the impactor's rocky core falls in to the core of the target. However, for higher angular momentum impacts, significant amounts become embedded anisotropically as lumps in the ice layer. Furthermore, most of the impactor's ice and energy is deposited in a hot, high-entropy shell at a radius of ~3 Earth radii. This could explain Uranus' observed lack of heat flow from the interior and be relevant for understanding its asymmetric magnetic field. We verify the results from the single previous study of lower resolution simulations that an impactor with a mass of at least 2 Earth masses can produce sufficiently rapid rotation in the post-impact Uranus for a range of angular momenta. At least 90% of the atmosphere remains bound to the final planet after the collision, but over half can be ejected beyond the Roche radius by a 2 or 3 Earth mass impactor. This atmospheric erosion peaks for intermediate impactor angular momenta (~3*10^36 kg m^2 s^-1). Rock is more efficiently placed into orbit and made available for satellite formation by 2 Earth mass impactors than 3 Earth mass ones, because it requires tidal disruption that is suppressed by the more massive impactors. △ Less

Submitted 3 July, 2018; v1 submitted 19 March, 2018; originally announced March 2018.

Comments: Published in ApJ. 14 pages, 11 figures

Report number: LA-UR-18-22206

Journal ref: ApJ 861 (2018) 52

The longevity of water ice on Ganymedes and Europas around migrated giant planets

Authors: Owen R. Lehmer, David C. Catling, Kevin J. Zahnle

Abstract: The gas giant planets in the Solar System have a retinue of icy moons, and we expect giant exoplanets to have similar satellite systems. If a Jupiter-like planet were to migrate toward its parent star the icy moons orbiting it would evaporate, creating atmospheres and possible habitable surface oceans. Here, we examine how long the surface ice and possible oceans would last before being hydrodynam… ▽ More The gas giant planets in the Solar System have a retinue of icy moons, and we expect giant exoplanets to have similar satellite systems. If a Jupiter-like planet were to migrate toward its parent star the icy moons orbiting it would evaporate, creating atmospheres and possible habitable surface oceans. Here, we examine how long the surface ice and possible oceans would last before being hydrodynamically lost to space. The hydrodynamic loss rate from the moons is determined, in large part, by the stellar flux available for absorption, which increases as the giant planet and icy moons migrate closer to the star. At some planet-star distance the stellar flux incident on the icy moons becomes so great that they enter a runaway greenhouse state. This runaway greenhouse state rapidly transfers all available surface water to the atmosphere as vapor, where it is easily lost from the small moons. However, for icy moons of Ganymede's size around a Sun-like star we found that surface water (either ice or liquid) can persist indefinitely outside the runaway greenhouse orbital distance. In contrast, the surface water on smaller moons of Europa's size will only persist on timescales greater than 1 Gyr at distances ranging 1.49 to 0.74 AU around a Sun-like star for Bond albedos of 0.2 and 0.8, where the lower albedo becomes relevant if ice melts. Consequently, small moons can lose their icy shells, which would create a torus of H atoms around their host planet that might be detectable in future observations. △ Less

Submitted 5 April, 2017; v1 submitted 20 March, 2017; originally announced March 2017.

Comments: Accepted by The Astrophysical Journal

The cosmic shoreline: the evidence that escape determines which planets have atmospheres, and what this may mean for Proxima Centauri b

Authors: Kevin J. Zahnle, David C. Catling

Abstract: The planets of the Solar System divide neatly between those with atmospheres and those without when arranged by insolation ($I$) and escape velocity ($v_{\mathrm{esc}}$). The dividing line goes as $I \propto v_{\mathrm{esc}}^4$. Exoplanets with reported masses and radii are shown to crowd against the extrapolation of the Solar System trend, making a metaphorical cosmic shoreline that unites all th… ▽ More The planets of the Solar System divide neatly between those with atmospheres and those without when arranged by insolation ($I$) and escape velocity ($v_{\mathrm{esc}}$). The dividing line goes as $I \propto v_{\mathrm{esc}}^4$. Exoplanets with reported masses and radii are shown to crowd against the extrapolation of the Solar System trend, making a metaphorical cosmic shoreline that unites all the planets. The $I \propto v_{\mathrm{esc}}^4$ relation may implicate thermal escape. We therefore address the general behavior of hydrodynamic thermal escape models ranging from Pluto to highly-irradiated Extrasolar Giant Planets (EGPs). Energy-limited escape is harder to test because copious XUV radiation is mostly a feature of young stars, and hence requires extrapolating to historic XUV fluences ($I_{\mathrm{xuv}}$) using proxies and power laws. An energy-limited shoreline should scale as $I_{\mathrm{xuv}} \propto v_{\mathrm{esc}}^3\sqrtρ$, which differs distinctly from the apparent $I_{\mathrm{xuv}} \propto v_{\mathrm{esc}}^4$ relation. Energy-limited escape does provide good quantitative agreement to the highly irradiated EGPs. Diffusion-limited escape implies that no planet can lose more than 1% of its mass as H$_2$. Impact erosion, to the extent that impact velocities $v_{\mathrm{imp}}$ can be estimated for exoplanets, fits to a $v_{\mathrm{imp}} \approx 4\,-\,5\, v_{\mathrm{esc}}$ shoreline. The proportionality constant is consistent with what the collision of comet Shoemaker-Levy 9 showed us we should expect of modest impacts in deep atmospheres. With respect to the shoreline, Proxima Centauri b is on the metaphorical beach. Known hazards include its rapid energetic accretion, high impact velocities, its early life on the wrong side of the runaway greenhouse, and Proxima Centauri's XUV radiation. In its favor is a vast phase space of unknown unknowns. △ Less

Submitted 13 July, 2017; v1 submitted 11 February, 2017; originally announced February 2017.

Comments: 19 pages, 4 figures, 1 table

Journal ref: The Astrophysical Journal, 843:122 (23pp), 2017 July 10

Sulfur Hazes in Giant Exoplanet Atmospheres: Impacts on Reflected Light Spectra

Authors: Peter Gao, Mark S Marley, Kevin Zahnle, Tyler D Robinson, Nikole K Lewis

Abstract: Recent work has shown that sulfur hazes may arise in the atmospheres of some giant exoplanets due to the photolysis of H$_{2}$S. We investigate the impact such a haze would have on an exoplanet's geometric albedo spectrum and how it may affect the direct imaging results of WFIRST, a planned NASA space telescope. For temperate (250 K $<$ T$_{\rm eq}$ $<$ 700 K) Jupiter--mass planets, photochemical… ▽ More Recent work has shown that sulfur hazes may arise in the atmospheres of some giant exoplanets due to the photolysis of H$_{2}$S. We investigate the impact such a haze would have on an exoplanet's geometric albedo spectrum and how it may affect the direct imaging results of WFIRST, a planned NASA space telescope. For temperate (250 K $<$ T$_{\rm eq}$ $<$ 700 K) Jupiter--mass planets, photochemical destruction of H$_{2}$S results in the production of $\sim$1 ppmv of \seight between 100 and 0.1 mbar, which, if cool enough, will condense to form a haze. Nominal haze masses are found to drastically alter a planet's geometric albedo spectrum: whereas a clear atmosphere is dark at wavelengths between 0.5 and 1 $μ$m due to molecular absorption, the addition of a sulfur haze boosts the albedo there to $\sim$0.7 due to scattering. Strong absorption by the haze shortward of 0.4 $μ$m results in albedos $<$0.1, in contrast to the high albedos produced by Rayleigh scattering in a clear atmosphere. As a result, the color of the planet shifts from blue to orange. The existence of a sulfur haze masks the molecular signatures of methane and water, thereby complicating the characterization of atmospheric composition. Detection of such a haze by WFIRST is possible, though discriminating between a sulfur haze and any other highly reflective, high altitude scatterer will require observations shortward of 0.4 $μ$m, which is currently beyond WFIRST's design. △ Less

Submitted 24 February, 2017; v1 submitted 2 January, 2017; originally announced January 2017.

Comments: 16 pages, 8 figures, 1 table, accepted for publication in AJ

On the Composition of Young, Directly Imaged Giant Planets

Authors: J. I. Moses, M. S. Marley, K. Zahnle, M. R. Line, J. J. Fortney, T. S. Barman, C. Visscher, N. K. Lewis, M. J. Wolff

Abstract: The past decade has seen significant progress on the direct detection and characterization of young, self-luminous giant planets at wide orbital separations from their host stars. Some of these planets show evidence for disequilibrium processes like transport-induced quenching in their atmospheres; photochemistry may also be important, despite the large orbital distances. These disequilibrium chem… ▽ More The past decade has seen significant progress on the direct detection and characterization of young, self-luminous giant planets at wide orbital separations from their host stars. Some of these planets show evidence for disequilibrium processes like transport-induced quenching in their atmospheres; photochemistry may also be important, despite the large orbital distances. These disequilibrium chemical processes can alter the expected composition, spectral behavior, thermal structure, and cooling history of the planets, and can potentially confuse determinations of bulk elemental ratios, which provide important insights into planet-formation mechanisms. Using a thermo/photochemical kinetics and transport model, we investigate the extent to which disequilibrium chemistry affects the composition and spectra of directly imaged giant exoplanets. Results for specific "young Jupiters" such as HR 8799 b and 51 Eri b are presented, as are general trends as a function of planetary effective temperature, surface gravity, incident ultraviolet flux, and strength of deep atmospheric convection. We find that quenching is very important on young Jupiters, leading to CO/CH4 and N2/NH3 ratios much greater than, and H2O mixing ratios a factor of a few less than, chemical-equilibrium predictions. Photochemistry can also be important on such planets, with CO2 and HCN being key photochemical products. Carbon dioxide becomes a major constituent when stratospheric temperatures are low and recycling of water via the H2 + OH reaction becomes kinetically stifled. Young Jupiters with effective temperatures <~ 700 K are in a particularly interesting photochemical regime that differs from both transiting hot Jupiters and our own solar-system giant planets. △ Less

Submitted 30 August, 2016; originally announced August 2016.

Comments: Accepted in The Astrophysical Journal

Photolytic Hazes in the Atmosphere of 51 Eri b

Authors: Kevin J. Zahnle, Mark S. Marley, Caroline V. Morley, Julianne I. Moses

Abstract: We use a 1D model to address photochemistry and possible haze formation in the irradiated warm Jupiter, 51 Eridani b. The intended focus was to be carbon, but sulfur photochemistry turns out to be important. The case for organic photochemical hazes is intriguing but falls short of being compelling. If organic hazes form, they are likeliest to do so if vertical mixing in 51 Eri b is weaker than in… ▽ More We use a 1D model to address photochemistry and possible haze formation in the irradiated warm Jupiter, 51 Eridani b. The intended focus was to be carbon, but sulfur photochemistry turns out to be important. The case for organic photochemical hazes is intriguing but falls short of being compelling. If organic hazes form, they are likeliest to do so if vertical mixing in 51 Eri b is weaker than in Jupiter, and they would be found below the altitudes where methane and water are photolyzed. The more novel result is that photochemistry turns H$_2$S into elemental sulfur, here treated as S$_8$. In the cooler models, S$_8$ is predicted to condense in optically thick clouds of solid sulfur particles, whilst in the warmer models S$_8$ remains a vapor along with several other sulfur allotropes that are both visually striking and potentially observable. For 51 Eri b, the division between models with and without condensed sulfur is at an effective temperature of 700 K, which is within error its actual effective temperature; the local temperature where sulfur condenses is between 280 and 320 K. The sulfur photochemistry we have discussed is quite general and ought to be found in a wide variety of worlds over a broad temperature range, both colder and hotter than the 650-750 K range studied here, and we show that products of sulfur photochemistry will be nearly as abundant on planets where the UV irradiation is orders of magnitude weaker than it is on 51 Eri b. △ Less

Submitted 25 April, 2016; originally announced April 2016.

Comments: 24 pages including 11 figures and a table

Developing Atmospheric Retrieval Methods for Direct Imaging Spectroscopy of Gas Giants in Reflected Light I: Methane Abundances and Basic Cloud Properties

Authors: Roxana E. Lupu, Mark S. Marley, Nikole Lewis, Michael Line, Wesley A. Traub, Kevin Zahnle

Abstract: Upcoming space-based coronagraphic instruments in the next decade will perform reflected light spectroscopy and photometry of cool, directly imaged extrasolar giant planets. We are developing a new atmospheric retrieval methodology to help assess the science return and inform the instrument design for such future missions, and ultimately interpret the resulting observations. Our retrieval techniqu… ▽ More Upcoming space-based coronagraphic instruments in the next decade will perform reflected light spectroscopy and photometry of cool, directly imaged extrasolar giant planets. We are developing a new atmospheric retrieval methodology to help assess the science return and inform the instrument design for such future missions, and ultimately interpret the resulting observations. Our retrieval technique employs a geometric albedo model coupled with both a Markov chain Monte Carlo Ensemble Sampler (emcee) and a multimodal nested sampling algorithm (MultiNest) to map the posterior distribution. This combination makes the global evidence calculation more robust for any given model, and highlights possible discrepancies in the likelihood maps. As a proof-of-concept, our current atmospheric model contains 1 or 2 cloud layers, methane as a major absorber, and a H$_2$-He background gas. This 6-to-9 parameter model is appropriate for Jupiter-like planets and can be easily expanded in the future. In addition to deriving the marginal likelihood distribution and confidence intervals for the model parameters, we perform model selection to determine the significance of methane and cloud detection as a function of expected signal-to-noise in the presence of spectral noise correlations. After internal validation, the method is applied to realistic spectra of Jupiter, Saturn, and HD 99492 c, a model observing target. We find that the presence or absence of clouds and methane can be determined with high confidence, while parameter uncertainties are model-dependent and correlated. Such general methods will also be applicable to the interpretation of direct imaging spectra of cloudy terrestrial planets. △ Less

Submitted 16 September, 2016; v1 submitted 18 April, 2016; originally announced April 2016.

Comments: 37 pages, 30 figures, accepted for publication in AJ

Thermal Emission and Albedo Spectra of Super Earths with Flat Transmission Spectra

Authors: Caroline V. Morley, Jonathan J. Fortney, Mark S. Marley, Kevin Zahnle, Michael Line, Eliza Kempton, Nikole Lewis, Kerri Cahoy

Abstract: Planets larger than Earth and smaller than Neptune are some of the most numerous in the galaxy, but observational efforts to understand this population have proved challenging because optically thick clouds or hazes at high altitudes obscure molecular features (Kreidberg et al. 2014b). We present models of super Earths that include thick clouds and hazes and predict their transmission, thermal emi… ▽ More Planets larger than Earth and smaller than Neptune are some of the most numerous in the galaxy, but observational efforts to understand this population have proved challenging because optically thick clouds or hazes at high altitudes obscure molecular features (Kreidberg et al. 2014b). We present models of super Earths that include thick clouds and hazes and predict their transmission, thermal emission, and reflected light spectra. Very thick, lofted clouds of salts or sulfides in high metallicity (1000x solar) atmospheres create featureless transmission spectra in the near-infrared. Photochemical hazes with a range of particle sizes also create featureless transmission spectra at lower metallicities. Cloudy thermal emission spectra have muted features more like blackbodies, and hazy thermal emission spectra have emission features caused by an inversion layer at altitudes where the haze forms. Close analysis of reflected light from warm (~400-800 K) planets can distinguish cloudy spectra, which have moderate albedos (0.05-0.20), from hazy models, which are very dark (0.0-0.03). Reflected light spectra of cold planets (~200 K) accessible to a space-based visible light coronagraph will have high albedos and large molecular features that will allow them to be more easily characterized than the warmer transiting planets. We suggest a number of complementary observations to characterize this population of planets, including transmission spectra of hot (>1000 K) targets, thermal emission spectra of warm targets using the James Webb Space Telescope (JWST), high spectral resolution (R~10^5) observations of cloudy targets, and reflected light spectral observations of directly-imaged cold targets. Despite the dearth of features observed in super Earth transmission spectra to date, different observations will provide rich diagnostics of their atmospheres. △ Less

Submitted 4 November, 2015; originally announced November 2015.

Comments: 23 pages, 23 figures. Revised for publication in The Astrophysical Journal

Rapid Water Loss can Extend the Lifetime of the Planetary Habitability

Authors: T. Kodama, H. Genda, Y. Abe, K. J. Zahnle

Abstract: Two habitable planetary states are proposed: an aqua planet like the Earth and a land planet that has a small amount of water. Land planets keep liquid water under larger solar radiation compared to aqua planets. Water loss may change an aqua planet into a land planet, and the planet can remain habitable for a longer time than had it stayed an aqua planet. We calculate planetary evolution with hyd… ▽ More Two habitable planetary states are proposed: an aqua planet like the Earth and a land planet that has a small amount of water. Land planets keep liquid water under larger solar radiation compared to aqua planets. Water loss may change an aqua planet into a land planet, and the planet can remain habitable for a longer time than had it stayed an aqua planet. We calculate planetary evolution with hydrogen escape for different initial water inventories and different distances from the central star. We find that there are two conditions to evolve an aqua planet into a land planet: the critical amount of water on the surface M_ml consistent with a planet being a land planet, and the critical amount of water vapor in the atmosphere M_cv that defines the onset of the runaway greenhouse state. We find that Earth-size aqua planets with initial oceans < 10 % of the Earth's can evolve into land planets if M_cv = 3 m in precipitable water and M_ml = 5 % of the Earth's ocean mass. Such planets can keep liquid water on their surface for another 2 Gyrs. The initial amount of water and M_cv are shown to be important dividing parameters of the planetary evolution path. Our results indicate that massive hydrogen escape could give a fresh start as another kind of habitable planet rather than the end of its habitability. △ Less

Submitted 12 September, 2015; originally announced September 2015.

Comments: Accepted for publication in The Astrophysical Journal, 31 pages, 11 figures

The Tethered Moon

Authors: Kevin J. Zahnle, Roxana Lupu, Anthony Dobrovolskis, Norman H. Sleep

Abstract: We address the thermal history of the Earth after the Moon-forming impact, taking tidal heating and thermal blanketing by the atmosphere into account. The atmosphere sets an upper bound of ~100 W/m^2 on how quickly the Earth can cool. The liquid magma ocean cools over 2-10 Myrs, with longer times corresponding to high angular-momentum events. Tidal heating is focused mostly in mantle materials tha… ▽ More We address the thermal history of the Earth after the Moon-forming impact, taking tidal heating and thermal blanketing by the atmosphere into account. The atmosphere sets an upper bound of ~100 W/m^2 on how quickly the Earth can cool. The liquid magma ocean cools over 2-10 Myrs, with longer times corresponding to high angular-momentum events. Tidal heating is focused mostly in mantle materials that are just beginning to freeze. The atmosphere's control over cooling sets up a negative feedback between viscosity-dependent tidal heating and temperature-dependent viscosity of the magma ocean. While the feedback holds, evolution of the Moon's orbit is limited by the modest radiative cooling rate of Earth's atmosphere. Orbital evolution is orders of magnitude slower than in conventional constant Q models, which promotes capture by resonances. The evection resonance is encountered early, when the Earth is molten. Capture by the evection resonance appears certain but unlikely to generate much eccentricity because it is encountered early when the Earth is molten and Q_Earth >> Q_Moon. Tidal dissipation in the Earth becomes more efficient (Q_Earth << Q_Moon) later when the Moon is between ~20 R_Earth and ~40 R_Earth. If lunar eccentricity grew great, this was when it did so, perhaps setting the table for some other process to leave its mark on the inclination of the Moon. △ Less

Submitted 6 August, 2015; originally announced August 2015.

Comments: 34 pages, 10 figures, preprint format

Journal ref: Earth and Planetary Science Letters 427, 74-82 (2015)

HST hot-Jupiter transmission spectral survey: Haze in the atmosphere of WASP-6b

Authors: N. Nikolov, D. K. Sing, A. S. Burrows, J. J. Fortney, G. W. Henry, F. Pont, G. E. Ballester, S. Aigrain, P. A. Wilson, C. M. Huitson, N. P. Gibson, J. -M. Desert, A. Lecavelier des Etangs, A. P. Showman, A. Vidal-Madjar, H. R. Wakeford, K. Zahnle

Abstract: We report Hubble Space Telescope (HST) optical to near-infrared transmission spectroscopy of the hot Jupiter WASP-6b, measured with the Space Telescope Imaging Spectrograph (STIS) and Spitzer's InfraRed Array Camera (IRAC). The resulting spectrum covers the range $0.29-4.5\,μ$m. We find evidence for modest stellar activity of WASP-6b and take it into account in the transmission spectrum. The overa… ▽ More We report Hubble Space Telescope (HST) optical to near-infrared transmission spectroscopy of the hot Jupiter WASP-6b, measured with the Space Telescope Imaging Spectrograph (STIS) and Spitzer's InfraRed Array Camera (IRAC). The resulting spectrum covers the range $0.29-4.5\,μ$m. We find evidence for modest stellar activity of WASP-6b and take it into account in the transmission spectrum. The overall main characteristic of the spectrum is an increasing radius as a function of decreasing wavelength corresponding to a change of $Δ(R_p/R_{\ast})=0.0071$ from 0.33 to $4.5\,μ$m. The spectrum suggests an effective extinction cross-section with a power law of index consistent with Rayleigh scattering, with temperatures of $973\pm144$ K at the planetary terminator. We compare the transmission spectrum with hot-Jupiter atmospheric models including condensate-free and aerosol-dominated models incorporating Mie theory. While none of the clear-atmosphere models is found to be in good agreement with the data, we find that the complete spectrum can be described by models that include significant opacity from aerosols including Fe-poor Mg$_2$SiO$_4$, MgSiO$_3$, KCl and Na$_2$S dust condensates. WASP-6b is the second planet after HD189733b which has equilibrium temperatures near $\sim1200$ K and shows prominent atmospheric scattering in the optical. △ Less

Submitted 17 November, 2014; originally announced November 2014.

Comments: 18 pages, 15 figures, 7 tables

Methane, Carbon Monoxide, and Ammonia in Brown Dwarfs and Self-Luminous Giant Planets

Authors: Kevin J. Zahnle, Mark S. Marley

Abstract: We address disequilibrum abundances of some simple molecules in the atmospheres of solar composition brown dwarfs and self-luminous extrasolar giant planets using a kinetics-based 1D atmospheric chemistry model. Our approach is to use the full kinetics model to survey the parameter space with effective temperatures between 500 K and 1100 K. In all of these worlds equilibrium chemistry favors CH4 o… ▽ More We address disequilibrum abundances of some simple molecules in the atmospheres of solar composition brown dwarfs and self-luminous extrasolar giant planets using a kinetics-based 1D atmospheric chemistry model. Our approach is to use the full kinetics model to survey the parameter space with effective temperatures between 500 K and 1100 K. In all of these worlds equilibrium chemistry favors CH4 over CO in the parts of the atmosphere that can be seen from Earth, but in most disequilibrium favors CO. The small surface gravity of a planet strongly discriminates against CH4 when compared to an otherwise comparable brown dwarf. If vertical mixing is like Jupiter's, the transition from methane to CO occurs at 500 K in a planet. Sluggish vertical mixing can raise this to 600 K; but clouds or more vigorous vertical mixing could lower this to 400 K. The comparable thresholds in brown dwarfs are $1100\pm100$ K. Ammonia is also sensitive to gravity, but unlike CH4/CO, the NH3/N2 ratio is insensitive to mixing, which makes NH3 a potential proxy for gravity. HCN may become interesting in high gravity brown dwarfs with very strong vertical mixing. Detailed analysis of the CO-CH4 reaction network reveals that the bottleneck to CO hydrogenation goes through methanol, in partial agreement with previous work. Simple, easy to use quenching relations are derived by fitting to the complete chemistry of the full ensemble of models. These relations are valid for determining CO, CH4, NH3, HCN, and CO2 abundances in the range of self-luminous worlds we have studied but may not apply if atmospheres are strongly heated at high altitudes by processes not considered here (e.g., wave breaking). △ Less

Submitted 10 October, 2014; v1 submitted 26 August, 2014; originally announced August 2014.

Comments: Astrophysical Journal, in press. Clarity improvements throughout and one new figure. 17 figures, 20 pages

The Atmospheres of Earth-like Planets after Giant Impact Events

Authors: R. E. Lupu, Kevin Zahnle, Mark S. Marley, Laura Schaefer, Bruce Fegley, Caroline Morley, Kerri Cahoy, Richard Freedman, Jonathan J. Fortney

Abstract: It is now understood that the accretion of terrestrial planets naturally involves giant collisions, the moon-forming impact being a well known example. In the aftermath of such collisions the surface of the surviving planet is very hot and potentially detectable. Here we explore the atmospheric chemistry, photochemistry, and spectral signatures of post-giant-impact terrestrial planets enveloped by… ▽ More It is now understood that the accretion of terrestrial planets naturally involves giant collisions, the moon-forming impact being a well known example. In the aftermath of such collisions the surface of the surviving planet is very hot and potentially detectable. Here we explore the atmospheric chemistry, photochemistry, and spectral signatures of post-giant-impact terrestrial planets enveloped by thick atmospheres consisting predominantly of CO2, and H2O. The atmospheric chemistry and structure are computed self-consistently for atmospheres in equilibrium with hot surfaces with composition reflecting either the bulk silicate Earth (which includes the crust, mantle, atmosphere and oceans) or Earth's continental crust. We account for all major molecular and atomic opacity sources including collision-induced absorption. We find that these atmospheres are dominated by H2O and CO2, while the formation of CH4, and NH3 is quenched due to short dynamical timescales. Other important constituents are HF, HCl, NaCl, and SO2. These are apparent in the emerging spectra, and can be indicative that an impact has occurred. The use of comprehensive opacities results in spectra that are a factor of 2 lower in surface brightness in the spectral windows than predicted by previous models. The estimated luminosities show that the hottest post-giant-impact planets will be detectable with near-infrared coronagraphs on the planned 30m-class telescopes. The 1-4um region will be most favorable for such detections, offering bright features and better contrast between the planet and a potential debris disk. We derive cooling timescales on the order of 10^5-10^6 Myrs, based on the modeled effective temperatures. This leads to the possibility of discovering tens of such planets in future surveys. △ Less

Submitted 6 January, 2014; originally announced January 2014.

Comments: 51 pages, 16 figures, preprint format; ApJ submitted

HST hot Jupiter Transmission Spectral Survey: A detection of Na and strong optical absorption in HAT-P-1b

Authors: N. Nikolov, D. K. Sing, F. Pont, A. S. Burrows, J. J. Fortney, G. E. Ballester, T. M. Evans, C. M. Huitson, H. R. Wakeford, P. A. Wilson, S. Aigrain, D. Deming, N. P. Gibson, G. W. Henry, H. Knutson, A. Lecavelier des Etangs, A. P. Showman, A. Vidal-Madjar, K. Zahnle

Abstract: We present an optical to near-infrared transmission spectrum of the hot Jupiter HAT-P-1b, based on HST observations, covering the spectral regime from 0.29 to 1.027μm with STIS, which is coupled with a recent WFC3 transit (1.087 to 1.687μm). We derive refined physical parameters of the HAT-P-1 system, including an improved orbital ephemeris. The transmission spectrum shows a strong absorption sign… ▽ More We present an optical to near-infrared transmission spectrum of the hot Jupiter HAT-P-1b, based on HST observations, covering the spectral regime from 0.29 to 1.027μm with STIS, which is coupled with a recent WFC3 transit (1.087 to 1.687μm). We derive refined physical parameters of the HAT-P-1 system, including an improved orbital ephemeris. The transmission spectrum shows a strong absorption signature shortward of 0.55μm, with a strong blueward slope into the near-ultraviolet. We detect atmospheric sodium absorption at a 3.3σ significance level, but find no evidence for the potassium feature. The red data implies a marginally flat spectrum with a tentative absorption enhancement at wavelength longer than ~0.85μm. The STIS and WFC3 spectra differ significantly in absolute radius level (4.3 +/- 1.6 pressure scale heights), implying strong optical absorption in the atmosphere of HAT-P-1b. The optical to near-infrared difference cannot be explained by stellar activity, as simulta- neous stellar activity monitoring of the G0V HAT-P-1b host star and its identical companion show no significant activity that could explain the result. We compare the complete STIS and WFC3 transmission spectrum with theoretical atmospheric mod- els which include haze, sodium and an extra optical absorber. We find that both an optical absorber and a super-solar sodium to water abundance ratio might be a sce- nario explaining the HAT-P-1b observations. Our results suggest that strong optical absorbers may be a dominant atmospheric feature in some hot Jupiter exoplanets. △ Less

Submitted 30 September, 2013; originally announced October 2013.

HST hot Jupiter transmission spectral survey: evidence for aerosols and lack of TiO in the atmosphere of WASP-12b

Authors: D. K. Sing, A. Lecavelier des Etangs, J. J. Fortney, A. S. Burrows, F. Pont, H. R. Wakeford, G. E. Ballester, N. Nikolov, G. W. Henry, S. Aigrain, D. Deming, T. M. Evans, N. P. Gibson, C. M. Huitson, H. Knutson, A. P. Showman, A. Vidal-Madjar, P. A. Wilson, M. H. Williamson, K. Zahnle

Abstract: We present HST optical transmission spectra of the transiting hot Jupiter WASP-12b, taken with the STIS instrument. From the transmission spectra, we are able to decisively rule out prominent absorption by TiO in the exoplanet's atmosphere. Strong pressure-broadened Na and K absorption signatures are also excluded, as are significant metal-hydride features. We compare our combined broadband spectr… ▽ More We present HST optical transmission spectra of the transiting hot Jupiter WASP-12b, taken with the STIS instrument. From the transmission spectra, we are able to decisively rule out prominent absorption by TiO in the exoplanet's atmosphere. Strong pressure-broadened Na and K absorption signatures are also excluded, as are significant metal-hydride features. We compare our combined broadband spectrum to a wide variety of existing aerosol-free atmospheric models, though none are satisfactory fits. However, we do find that the full transmission spectrum can be described by models which include significant opacity from aerosols: including Rayleigh scattering, Mie scattering, tholin haze, and settling dust profiles. The transmission spectrum follows an effective extinction cross section with a power-law of index alpha, with the slope of the transmission spectrum constraining the quantity alphaT = -3528+/-660 K, where T is the atmospheric temperature. Rayleigh scattering (alpha=-4) is among the best fitting models, though requires low terminator temperatures near 900 K. Sub-micron size aerosol particles can provide equally good fits to the entire transmission spectrum for a wide range of temperatures, and we explore corundum as a plausible dust aerosol. The presence of atmospheric aerosols also helps to explain the modestly bright albedo implied by Spitzer observations, as well as the near black body nature of the emission spectrum. Ti-bearing condensates on the cooler night-side is the most natural explanation for the overall lack of TiO signatures in WASP-12b, indicating the day/night cold-trap is an important effect for very hot Jupiters. These finding indicate that aerosols can play a significant atmospheric role for the entire wide range of hot-Jupiter atmospheres, potentially affecting their overall spectrum and energy balance.(abridged) △ Less

Submitted 20 September, 2013; originally announced September 2013.

Comments: 19 pages, 14 figures, 5 tables. Accepted for publication in MNRAS

Hubble Space Telescope hot Jupiter Transmission Spectral Survey: detection of water in HAT-P-1b from Wide Field Camera 3 near-infrared spatial scan observations

Authors: H. R. Wakeford, D. K. Sing, D. Deming, N. P. Gibson, J. J. Fortney, A. S. Burrows, G. Ballester, N. Nikolov, S. Aigrain, G. Henry, H. Knutson, A. Lecavelier des Etangs, F. Pont, A. P. Showman, A. Vidal-Madjar, K. Zahnle

Abstract: We present Hubble Space Telescope near-infrared transmission spectroscopy of the transiting hot-Jupiter HAT-P-1b. We observed one transit with Wide Field Camera 3 using the G141 low-resolution grism to cover the wavelength range 1.087- 1.678 μm. These time series observations were taken with the newly available spatial scan mode that increases the duty cycle by nearly a factor of two, thus improvi… ▽ More We present Hubble Space Telescope near-infrared transmission spectroscopy of the transiting hot-Jupiter HAT-P-1b. We observed one transit with Wide Field Camera 3 using the G141 low-resolution grism to cover the wavelength range 1.087- 1.678 μm. These time series observations were taken with the newly available spatial scan mode that increases the duty cycle by nearly a factor of two, thus improving the resulting photometric precision of the data. We measure a planet-to-star radius ratio of Rp/R*=0.11709+/-0.00038 in the white light curve with the centre of transit occurring at 2456114.345+/-0.000133 (JD). We achieve S/N levels per exposure of 1840 (0.061%) at a resolution of {Δλ}=19.2nm (R~70) in the 1.1173 - 1.6549μm spectral region, providing the precision necessary to probe the transmission spectrum of the planet at close to the resolution limit of the instrument. We compute the transmission spectrum using both single target and differential photometry with similar results. The resultant transmission spectrum shows a significant absorption above the 5-σ level matching the 1.4μm water absorption band. In solar composition models, the water absorption is sensitive to the ~1 mbar pressure levels at the terminator. The detected absorption agrees with that predicted by an 1000 K isothermal model, as well as with that predicted by a planetary-averaged temperature model. △ Less

Submitted 9 August, 2013; originally announced August 2013.

Comments: 14 pages, 17 figures, 3 tables, submitted April 4th to MNRAS, minor revisions

An HST Optical to Near-IR Transmission Spectrum of the Hot Jupiter WASP-19b: Detection of Atmospheric Water and Likely Absence of TiO

Authors: C. M. Huitson, D. K. Sing, F. Pont, J. J. Fortney, A. S. Burrows, P. A. Wilson, G. E. Ballester, N. Nikolov, N. P. Gibson, D. Deming, S. Aigrain, T. M. Evans, G. W. Henry, A. Lecavelier des Etangs, A. P. Showman, A. Vidal-Madjar, K. Zahnle

Abstract: We measure the transmission spectrum of WASP-19b from 3 transits using low-resolution optical spectroscopy from the HST Space Telescope Imaging Spectrograph (STIS). The STIS spectra cover a wavelength range of 0.29-1.03 microns with resolving power R=500. The optical data are combined with archival near-IR data from the HST Wide Field Camera 3 (WFC3) G141 grism, covering the wavelength range 1.087… ▽ More We measure the transmission spectrum of WASP-19b from 3 transits using low-resolution optical spectroscopy from the HST Space Telescope Imaging Spectrograph (STIS). The STIS spectra cover a wavelength range of 0.29-1.03 microns with resolving power R=500. The optical data are combined with archival near-IR data from the HST Wide Field Camera 3 (WFC3) G141 grism, covering the wavelength range 1.087-1.687 micron, with resolving power R=130. We obtain the transmission spectrum from 0.53-1.687 microns with S/N levels between 3000 and 11,000 in 0.1 micron bins. WASP-19 is a very active star, with optical stellar flux varying by a few per cent over time. We correct the transit light curves for the effects of stellar activity using ground-based activity monitoring with the Cerro Tololo Inter-American Observatory (CTIO). While we were not able to construct a transmission spectrum using the blue optical data due to the presence of large occulted star spots, we were able to use the spot crossings to help constrain the mean stellar spot temperature. To search for predicted features in the hot-Jupiter atmosphere, we also define spectral indices for differential radius measurements to specifically search for the presence of TiO and alkali line features. Our measurements rule out TiO features predicted for a planet of WASP-19b's equilibrium temperature (2050 K) in the transmission spectrum at the 2.7-2.9 sigma confidence level, depending on atmospheric model formalism. The WFC3 transmission spectrum shows strong absorption features due to the presence of H2O, which is detected at the 4 sigma confidence level between 1.1 and 1.4 microns. The results indicate that WASP-19b is a planet with no or low levels of TiO and without a high C/O ratio. The lack of observable TiO features are possibly due to rainout, breakdown from stellar activity or the presence of other absorbers in the optical. △ Less

Submitted 11 July, 2013; v1 submitted 8 July, 2013; originally announced July 2013.

Comments: Submitted 2013 February 28, accepted 2013 July 5

A Framework for Characterizing the Atmospheres of Low-Mass Low-Density Transiting Planets

Authors: Jonathan J. Fortney, Christoph Mordasini, Nadine Nettelmann, Eliza M. -R. Kempton, Thomas P. Greene, Kevin Zahnle

Abstract: We perform modeling investigations to aid in understanding the atmospheres and composition of small planets of ~2-4 Earth radii, which are now known to be common in our galaxy. GJ 1214b is a well studied example whose atmospheric transmission spectrum has been observed by many investigators. Here we take a step back from GJ 1214b to investigate the role that planetary mass, composition, and temper… ▽ More We perform modeling investigations to aid in understanding the atmospheres and composition of small planets of ~2-4 Earth radii, which are now known to be common in our galaxy. GJ 1214b is a well studied example whose atmospheric transmission spectrum has been observed by many investigators. Here we take a step back from GJ 1214b to investigate the role that planetary mass, composition, and temperature play in impacting the transmission spectra of these low-mass low-density (LMLD) planets. Under the assumption that these planets accrete modest hydrogen-dominated atmospheres and planetesimals, we use population synthesis models to show that predicted metal enrichments of the H/He envelope are high, with metal mass fraction Z_env values commonly 0.6 to 0.9, or ~100 to 400+ times solar. The high mean molecular weight of such atmospheres (μ~ 5-12) would naturally help to flatten the transmission spectrum of most LMLD planets. The high metal abundance would also provide significant condensible material for cloud formation. It is known that the H/He abundance in Uranus and Neptune decreases with depth, and we show that atmospheric evaporation of LMLD planets could expose atmospheric layers with gradually higher Z_env. However, values of Z_env close to solar composition can also arise, so diversity should be expected. Photochemically produced hazes, potentially due to methane photolysis, are another possibility for obscuring transmission spectra. Such hazes may not form above T_eq of ~800-1100 K, which is testable if such warm, otherwise low mean molecular weight atmospheres are stable against atmospheric evaporation. We find that available transmission data are consistent with relatively high mean molecular weight atmospheres for GJ 1214b and "warm Neptune" GJ 436b. We examine future prospects for characterizing GJ 1214b with Hubble and the James Webb Space Telescope. △ Less

Submitted 3 September, 2013; v1 submitted 18 June, 2013; originally announced June 2013.

Comments: Accepted to The Astrophysical Journal on August 12, 2013. Minor changes, including new-subsection on mass-radius of LMLD planets

Quantitatively Assessing the Role of Clouds in the Transmission Spectrum of GJ 1214b

Authors: Caroline V. Morley, Jonathan J. Fortney, Eliza M. -R. Kempton, Mark S. Marley, Channon Visscher, Kevin Zahnle

Abstract: Recent observations of the super-Earth GJ 1214b show that it has a relatively featureless transmission spectrum. One suggestion is that these observations indicate that the planet's atmosphere is vertically compact, perhaps due to a water-rich composition that yields a large mean molecular weight. Another suggestion is that the atmosphere is hydrogen/helium-rich with clouds that obscure predicted… ▽ More Recent observations of the super-Earth GJ 1214b show that it has a relatively featureless transmission spectrum. One suggestion is that these observations indicate that the planet's atmosphere is vertically compact, perhaps due to a water-rich composition that yields a large mean molecular weight. Another suggestion is that the atmosphere is hydrogen/helium-rich with clouds that obscure predicted absorption features. Previous models that incorporate clouds have included their effect without a strong physical motivation for their existence. Here, we present model atmospheres of GJ 1214b that include physically-motivated clouds of two types. We model the clouds that form as a result of condensation in chemical equilibrium, as they likely do on brown dwarfs, which include KCl and ZnS for this planet. We also include clouds that form as a result of photochemistry, forming a hydrocarbon haze layer. We use a photochemical kinetics model to understand the vertical distribution and available mass of haze-forming molecules. We model both solar and enhanced-metallicity cloudy models and determine the cloud properties necessary to match observations. In enhanced-metallicity atmospheres, we find that the equilibrium clouds can match the observations of GJ 1214b if they are lofted high into the atmosphere and have a low sedimentation efficiency (fsed=0.1). We find that models with a variety of hydrocarbon haze properties can match the observations. Particle sizes from 0.01 to 0.25 micron can match the transmission spectrum with haze-forming efficiencies as low as 1-5%. △ Less

Submitted 14 June, 2013; v1 submitted 17 May, 2013; originally announced May 2013.

Comments: Revised for ApJ, Figure 8 fixed in v2

A photochemical model for the carbon-rich planet WASP-12b

Authors: Ravi kumar Kopparapu, James F. Kasting, Kevin J. Zahnle

Abstract: The hot Jupiter WASP-12b is a heavily irradiated exoplanet in a short period orbit around a G0-star with twice the metallicity of the Sun. A recent thermochemical equilibrium analysis based on Spitzer and ground-based infrared observations suggests that the presence of $\ch4$ in its atmosphere and the lack of $\h2o$ features can only be explained if the carbon-to-oxygen ratio in the planet's atmos… ▽ More The hot Jupiter WASP-12b is a heavily irradiated exoplanet in a short period orbit around a G0-star with twice the metallicity of the Sun. A recent thermochemical equilibrium analysis based on Spitzer and ground-based infrared observations suggests that the presence of $\ch4$ in its atmosphere and the lack of $\h2o$ features can only be explained if the carbon-to-oxygen ratio in the planet's atmosphere is much greater than the solar ratio ($\ctoo = 0.54$). Here, we use a 1-D photochemical model to study the effect of disequilibrium chemistry on the observed abundances of $\h2o, \com, \co2$ and $\ch4$ in the WASP-12b atmosphere. We consider two cases: one with solar $\ctoo$ and another with $\ctoo = 1.08$. The solar case predicts that $\h2o$ and $\com$ are more abundant than $\co2$ and $\ch4$, as expected, whereas the high $\ctoo$ model shows that $\com$, C$_{2}$H$_{2}$ and HCN are more abundant. This indicates that the extra carbon from the high $\ctoo$ model is in hydrocarbon species. $\h2o$ photolysis is the dominant disequilibrium mechanism that alters the chemistry at higher altitudes in the solar $\ctoo$ case, whereas photodissociation of C$_{2}$H$_{2}$ and HCN is significant in the super-solar case. Furthermore, our analysis indicates that $\c2h2$ is the major absorber in the atmosphere of WASP-12b and the absorption features detected near 1.6 and 8 micron may be arising from C$_{2}$H$_{2}$ rather than $\ch4$. The Hubble Space Telescope's WFC3 can resolve this discrepancy, as $\c2h2$ has absorption between $1.51 - 1.54$ microns, while $\ch4$ does not. △ Less

Submitted 12 October, 2011; originally announced October 2011.

Comments: Accepted to Astrophysical Journal. Reaction list available with the journal version of the paper

Faint young Sun paradox remains

Authors: Colin Goldblatt, Kevin J. Zahnle

Abstract: The Sun was fainter when the Earth was young, but the climate was generally at least as warm as today; this is known as the `faint young Sun paradox'. Rosing et al. [1] claim that the paradox can be resolved by making the early Earth's clouds and surface less reflective. We show that, even with the strongest plausible assumptions, reducing cloud and surface albedos falls short by a factor of two o… ▽ More The Sun was fainter when the Earth was young, but the climate was generally at least as warm as today; this is known as the `faint young Sun paradox'. Rosing et al. [1] claim that the paradox can be resolved by making the early Earth's clouds and surface less reflective. We show that, even with the strongest plausible assumptions, reducing cloud and surface albedos falls short by a factor of two of resolving the paradox. A temperate Archean climate cannot be reconciled with the low level of CO2 suggested by Rosing et al. [1]; a stronger greenhouse effect is needed. △ Less

Submitted 30 May, 2011; v1 submitted 26 May, 2011; originally announced May 2011.

Comments: 3 pages, no figures. In press in Nature. v2 corrects typo in author list in original submission

The Role of Ejecta in the Small Crater Populations on the Mid-Sized Saturnian Satellites

Authors: Edward B. Bierhaus, Luke Dones, José Luis Alvarellos, Kevin Zahnle

Abstract: We find evidence that crater ejecta play an important role in the small crater populations on the Saturnian satellites, and more broadly, on cratered surfaces throughout the Solar System. We measure crater populations in Cassini images of Enceladus, Rhea, and Mimas, focusing on image data with scales less than 500 m/pixel. We use recent updates to crater scaling laws and their constants to estimat… ▽ More We find evidence that crater ejecta play an important role in the small crater populations on the Saturnian satellites, and more broadly, on cratered surfaces throughout the Solar System. We measure crater populations in Cassini images of Enceladus, Rhea, and Mimas, focusing on image data with scales less than 500 m/pixel. We use recent updates to crater scaling laws and their constants to estimate the amount of mass ejected in three different velocity ranges: (i) greater than escape velocity, (ii) less than escape velocity and faster than the minimum velocity required to make a secondary crater (v_min), and (iii) velocities less than v_min. Although the vast majority of mass on each satellite is ejected at speeds less than v_min, our calculations demonstrate that the differences in mass available in the other two categories should lead to observable differences in the small crater populations; the predictions are borne out by the measurements we have made to date. Rhea, Tethys, and Dione have sufficient surface gravities to retain ejecta moving fast enough to make secondary crater populations. The smaller satellites, such as Enceladus but especially Mimas, are expected to have little or no traditional secondary populations because their escape velocities are near the threshold velocity necessary to make a secondary crater. Our work clarifies why the Galilean satellites have extensive secondary crater populations relative to the Saturnian satellites. The presence, extent, and sizes of sesquinary craters (craters formed by ejecta that escape into temporary orbits around Saturn before re-impacting the surface) is not yet well understood. Finally, our work provides further evidence for a "shallow" size-frequency distribution (slope index of ~2 for a differential power-law) for comets a few km diameter and smaller. [slightly abbreviated] △ Less

Submitted 12 May, 2011; originally announced May 2011.

Comments: Submitted to Icarus. 77 double-spaced pages, including 25 figures and 5 tables

The Atmospheric Chemistry of GJ 1214b: Photochemistry and Clouds

Authors: Eliza Miller-Ricci Kempton, Kevin Zahnle, Jonathan J. Fortney

Abstract: Recent observations of the transiting super-Earth GJ 1214b reveal that its atmosphere may be hydrogen-rich or water-rich in nature, with clouds or hazes potentially affecting its transmission spectrum in the optical and very-near-IR. Here we further examine the possibility that GJ 1214b does indeed possess a hydrogen-dominated atmosphere, which is the hypothesis that is favored by models of the bu… ▽ More Recent observations of the transiting super-Earth GJ 1214b reveal that its atmosphere may be hydrogen-rich or water-rich in nature, with clouds or hazes potentially affecting its transmission spectrum in the optical and very-near-IR. Here we further examine the possibility that GJ 1214b does indeed possess a hydrogen-dominated atmosphere, which is the hypothesis that is favored by models of the bulk composition of the planet. We study the effects of non-equilibrium chemistry (photochemistry, thermal chemistry, and mixing) on the planet's transmission spectrum. We furthermore examine the possibility that clouds could play a significant role in attenuating GJ 1214b's transmission spectrum at short wavelengths. We find that non-equilibrium chemistry can have a large effect on the overall chemical composition of GJ 1214b's atmosphere, however these changes mostly take place above the height in the atmosphere that is probed by transmission spectroscopy. The effects of non-equilibrium chemistry on GJ 1214b's transmission spectrum are therefore minimal, with the largest effects taking place if the planet's atmosphere has super-solar metallicity and a low rate of vertical mixing. Interestingly, we find that the best fit to the observations of GJ 1214b's atmosphere in transmission occur if the planet's atmosphere is deficient in CH4, and possesses a cloud layer at a pressure of ~200 mbar. This is consistent with a picture of efficient methane photolysis, accompanied by formation of organic haze that obscures the lower atmosphere of GJ 1214b at optical wavelengths. However, for methane to be absent from GJ 1214b's transmission spectrum, UV photolysis of this molecule must be efficient at pressures of greater than ~1 mbar, whereas we find that methane only photolyzes to pressures less than 0.1 mbar, even under the most optimistic assumptions. (Abridged) △ Less

Submitted 22 December, 2011; v1 submitted 28 April, 2011; originally announced April 2011.

Comments: Accepted to ApJ; 32 pages, 8 figures, 1 table

Clouds and the Faint Young Sun Paradox

Authors: Colin Goldblatt, Kevin J. Zahnle

Abstract: We investigate the role which clouds could play in resolving the Faint Young Sun Paradox (FYSP). Lower solar luminosity in the past means that less energy was absorbed on Earth (a forcing of -50 Wm-2 during the late Archean), but geological evidence points to the Earth being at least as warm as it is today, with only very occasional glaciations. We perform radiative calculations on a single global… ▽ More We investigate the role which clouds could play in resolving the Faint Young Sun Paradox (FYSP). Lower solar luminosity in the past means that less energy was absorbed on Earth (a forcing of -50 Wm-2 during the late Archean), but geological evidence points to the Earth being at least as warm as it is today, with only very occasional glaciations. We perform radiative calculations on a single global mean atmospheric column. We select a nominal set of three layered, randomly overlapping clouds, which are both consistent with observed cloud climatologies and reproduce the observed global mean energy budget of Earth. By varying the fraction, thickness, height and particle size of these clouds we conduct a wide exploration of how changed clouds could affect climate, thus constraining how clouds could contribute to resolving the FYSP. Low clouds reflect sunlight but have little greenhouse effect. Removing them entirely gives a~forcing of +25 Wm-2 whilst more modest reduction in their efficacy gives a forcing of +10 to +15 Wm-2. For high clouds, the greenhouse effect dominates. It is possible to generate +50 Wm-2 forcing from enhancing these, but this requires making them 3.5 times thicker and 14 K colder than the standard high cloud in our nominal set and expanding their coverage to 100% of the sky. Such changes are not credible. More plausible changes would generate no more that +15 Wm-2 forcing. Thus neither fewer low clouds nor more high clouds can provide enough forcing to resolve the FYSP. Decreased surface albedo can contribute no more than +5 Wm-2 forcing. Some models which have been applied to the FYSP do not include clouds at all. These overestimate the forcing due to increased CO2 by 20 to 25% when CO2 is 0.01 to 0.1 bar. △ Less

Submitted 15 February, 2011; originally announced February 2011.

Comments: 18 pages, 13 figures. Accepted for publication in "Climate of the Past". CP is a two stage journal: for the preceding "Climate of the Past Discussions" article, see http://www.clim-past-discuss.net/6/1163/2010/cpd-6-1163-2010.html

On the volatile enrichments and heavy element content in HD 189733b

Authors: O. Mousis, J. I. Lunine, J. -M. Petit, K. Zahnle, L. Biennier, S. Picaud, T. V. Johnson, J. B. A. Mitchell, V. Boudon, D. Cordier, M. Devel, R. Georges, C. Griffith, N. Iro, M. S. Marley, U. Marboeuf

Abstract: Favored theories of giant planet formation center around two main paradigms, namely the core accretion model and the gravitational instability model. These two formation scenarios support the hypothesis that the giant planet metallicities should be higher or equal to that of the parent star. Meanwhile, spectra of the transiting hot Jupiter HD189733b suggest that carbon and oxygen abundances range… ▽ More Favored theories of giant planet formation center around two main paradigms, namely the core accretion model and the gravitational instability model. These two formation scenarios support the hypothesis that the giant planet metallicities should be higher or equal to that of the parent star. Meanwhile, spectra of the transiting hot Jupiter HD189733b suggest that carbon and oxygen abundances range from depleted to enriched with respect to the star. Here, using a model describing the formation sequence and composition of planetesimals in the protoplanetary disk, we determine the range of volatile abundances in the envelope of HD189733b that is consistent with the 20--80 Earth-masses of heavy elements estimated to be present in the planet's envelope. We then compare the inferred carbon and oxygen abundances to those retrieved from spectroscopy and we find a range of supersolar values that directly fit both spectra and internal structure models. In some cases, we find that the apparent contradiction between the subsolar elemental abundances and the {mass of heavy elements predicted in HD189733b by internal structure models} can be explained by the presence of large amounts of carbon molecules in the form of polycyclic aromatic hydrocarbons and soots in the upper layers of the envelope, as suggested by recent photochemical models. A diagnostic test that would confirm the presence of these compounds in the envelope is the detection of acetylene. △ Less

Submitted 18 November, 2010; originally announced November 2010.

Comments: Accepted for publication in The Astrophysical Journal

Exchange of ejecta between Telesto and Calypso: Tadpoles, horseshoes, and passing orbits

Authors: Anthony R. Dobrovolskis, Jose' Luis Alvarellos, Kevin J. Zahnle, Jack J. Lissauer

Abstract: We have numerically integrated the orbits of ejecta from Telesto and Calypso, the two small Trojan companions of Saturn's major satellite Tethys. Ejecta were launched with speeds comparable to or exceeding their parent's escape velocity, consistent with impacts into regolith surfaces. We find that the fates of ejecta fall into several distinct categories, depending on both the speed and direction… ▽ More We have numerically integrated the orbits of ejecta from Telesto and Calypso, the two small Trojan companions of Saturn's major satellite Tethys. Ejecta were launched with speeds comparable to or exceeding their parent's escape velocity, consistent with impacts into regolith surfaces. We find that the fates of ejecta fall into several distinct categories, depending on both the speed and direction of launch. The slowest ejecta follow sub-orbital trajectories and re-impact their source moon in less than one day. Slightly faster debris barely escape their parent's Hill sphere and are confined to tadpole orbits, librating about Tethys' triangular Lagrange points L4 (leading, near Telesto) or L5 (trailing, near Calypso) with nearly the same orbital semi-major axis as Tethys, Telesto, and Calypso. These ejecta too eventually re-impact their source moon, but with a median lifetime of a few dozen years. Those which re-impact within the first ten years or so have lifetimes near integer multiples of 348.6 days (half the tadpole period). Still faster debris with azimuthal velocity components >~ 10 m/s enter horseshoe orbits which enclose both L4 and L5 as well as L3, but which avoid Tethys and its Hill sphere. These ejecta impact either Telesto or Calypso at comparable rates, with median lifetimes of several thousand years. However, they cannot reach Tethys itself; only the fastest ejecta, with azimuthal velocities >~ 40 m/s, achieve "passing orbits" which are able to encounter Tethys. Tethys accretes most of these ejecta within several years, but some 1 % of them are scattered either inward to hit Enceladus or outward to strike Dione, over timescales on the order of a few hundred years. △ Less

Submitted 24 May, 2010; originally announced May 2010.

Planetary Formation and Evolution Revealed with a Saturn Entry Probe: The Importance of Noble Gases

Authors: Jonathan J. Fortney, Kevin Zahnle, Isabelle Baraffe, Adam Burrows, Sarah E. Dodson-Robinson, Gilles Chabrier, Tristan Guillot, Ravit Helled, Franck Hersant, William B. Hubbard, Jack J. Lissauer, Mark S. Marley

Abstract: The determination of Saturn's atmospheric noble gas abundances are critical to understanding the formation and evolution of Saturn, and giant planets in general. These measurements can only be performed with an entry probe. A Saturn probe will address whether enhancement in heavy noble gases, as was found in Jupiter, are a general feature of giant planets, and their ratios will be a powerful con… ▽ More The determination of Saturn's atmospheric noble gas abundances are critical to understanding the formation and evolution of Saturn, and giant planets in general. These measurements can only be performed with an entry probe. A Saturn probe will address whether enhancement in heavy noble gases, as was found in Jupiter, are a general feature of giant planets, and their ratios will be a powerful constraint on how they form. The helium abundance will show the extent to which helium has phase separated from hydrogen in the planet's deep interior. Jupiter's striking neon depletion may also be tied to its helium depletion, and must be confirmed or refuted in Saturn. Together with Jupiter's measured atmospheric helium abundance, a consistent evolutionary theory for both planets, including "helium rain" will be possible. We will then be able to calibrate the theory of the evolution of all giant planets, including exoplanets. In addition, high pressure H/He mixtures under giant planet conditions are an important area of condensed matter physics that are beyond the realm of experiment. △ Less

Submitted 18 November, 2009; originally announced November 2009.

Comments: 7 pages. Submitted to the Giant Planets panel of the 2013-2022 Planetary Science Decadal Survey

Thermometric Soots on Warm Jupiters?

Authors: K. Zahnle, M. S. Marley, J. J. Fortney

Abstract: We use a 1D thermochemical and photochemical kinetics model to predict the disequilibrium stratospheric chemistries of warm and hot Jupiters (800 < T < 1200 K). Thermal chemistry and vertical mixing are generally more important than photochemistry. At 1200 K, methane is oxidized to CO and CO2 by OH radicals from thermal decomposition of water. At T < 1000 K, methane is reactive but stable enough t… ▽ More We use a 1D thermochemical and photochemical kinetics model to predict the disequilibrium stratospheric chemistries of warm and hot Jupiters (800 < T < 1200 K). Thermal chemistry and vertical mixing are generally more important than photochemistry. At 1200 K, methane is oxidized to CO and CO2 by OH radicals from thermal decomposition of water. At T < 1000 K, methane is reactive but stable enough to reach the stratosphere, while water is stable enough that OH levels are suppressed by reaction with H2. These trends raise the effective C/O ratio in the reacting gases above unity. Reduced products such as ethylene, acetylene, and hydrogen cyanide become abundant; further polymerization should lead to formation of PAHs (Poly-Aromatic Hydrocarbons) and soots. Parallel shifts are seen in the sulfur chemistry, with CS and CS2 displacing S2 and HS as the interesting disequilibrium products. Although lower temperature is a leading factor favoring hydrocarbons, higher rates of vertical mixing, lower metallicities, and lower incident UV radiation also favor organic synthesis. Acetylene (the first step toward PAHs) formation is especially favored by high eddy diffusion coefficients Kzz > 10^10 cm2/s. In most cases planetary compositions inferred from transit observations will differ markedly from those inferred from reflected or emitted light from the same planet. The peculiar properties of HD 189733b compared to other hot Jupiters - a broadband blue haze, little sign of Na or K, and hints of low metallicity - can be explained by an organic haze if the planet is cool enough. Whether this interpretation applies to HD 189733b itself, many organic-rich warm Jupiters are sure to be discovered in the near future. △ Less

Submitted 15 July, 2010; v1 submitted 3 November, 2009; originally announced November 2009.

Comments: Revision of paper originally submitted to Ap J Letters, now redirected to Ap J. Revised version now features much more background and discussion, a slightly changed title, and a helpful 29 page table of reaction rates not included here. Conclusions are essentially unchanged

Atmospheric Sulfur Photochemistry on Hot Jupiters

Authors: K. Zahnle, M. S. Marley, R. S. Freedman, K. Lodders, J. J. Fortney

Abstract: We develop a new 1D photochemical kinetics code to address stratospheric chemistry and stratospheric heating in hot Jupiters. Here we address optically active S-containing species and CO2 at 1200 < T < 2000 K. HS (mercapto) and S2 are highly reactive species that are generated photochemically and thermochemically from H2S with peak abundances between 1-10 mbar. S2 absorbs UV between 240 and 340… ▽ More We develop a new 1D photochemical kinetics code to address stratospheric chemistry and stratospheric heating in hot Jupiters. Here we address optically active S-containing species and CO2 at 1200 < T < 2000 K. HS (mercapto) and S2 are highly reactive species that are generated photochemically and thermochemically from H2S with peak abundances between 1-10 mbar. S2 absorbs UV between 240 and 340 nm and is optically thick for metallicities [SH] > 0 at T > 1200 K. HS is probably more important than S2, as it is generally more abundant than S2 under hot Jupiter conditions and it absorbs at somewhat redder wavelengths. We use molecular theory to compute an HS absorption spectrum from sparse available data and find that HS should absorb strongly between 300 and 460 nm, with absorption at the longer wavelengths being temperature sensitive. When the two absorbers are combined, radiative heating (per kg of gas) peaks at 100 microbars, with a total stratospheric heating of about 8 x 10^4 W/m^2 for a jovian planet orbiting a solar-twin at 0.032 AU. Total heating is insensitive to metallicity. The CO2 mixing ratio is a well-behaved quadratic function of metallicity, ranging from 1.6 x 10^-8 to 1.6 x 10^-4 for -0.3 < [M/H] < 1.7. CO2 is insensitive to insolation, vertical mixing, temperature (1200 < T <2000 K), and gravity. The photochemical calculations confirm that CO2 should prove a useful probe of planetary metallicity. △ Less

Submitted 26 June, 2009; v1 submitted 9 March, 2009; originally announced March 2009.

Comments: Astrophysical Journal Lett. in press; important revision includes effect of updated thermodynamic data and a new opacity source

Journal ref: Astrophys.J.701:L20-L24,2009

Material enhancement in protoplanetary nebulae by particle drift through evaporation fronts

Authors: J. N. Cuzzi, K. J. Zahnle

Abstract: Solid material in a protoplanetary nebula is subject to vigorous redistribution processes relative to the nebula gas. Meter-sized particles drift rapidly inwards near the nebula midplane, and material evaporates when the particles cross a condensation/evaporation boundary. The material cannot be removed as fast in its vapor form as it is being supplied in solid form, so its concentration increas… ▽ More Solid material in a protoplanetary nebula is subject to vigorous redistribution processes relative to the nebula gas. Meter-sized particles drift rapidly inwards near the nebula midplane, and material evaporates when the particles cross a condensation/evaporation boundary. The material cannot be removed as fast in its vapor form as it is being supplied in solid form, so its concentration increases locally by a large factor (more than an order of magnitude under nominal conditions). As time goes on, the vapor phase enhancement propagates for long distances inside the evaporation boundary (potentially all the way in to the star). Meanwhile, material is enhanced in its solid form over a characteristic lengthscale outside the evaporation boundary. This effect is applicable to any condensible (water, silicates, {\it etc.}). Three distinct radial enhancement/depletion regimes can be discerned by use of a simple model. Meteoritics applications include oxygen fugacity and isotopic variations, as well as isotopic homogenization in silicates. Planetary system applications include more robust enhancement of solids in Jupiter's core formation region than previously suggested. Astrophysical applications include differential, time-dependent enhancement of vapor phase CO and H$_2$O in the terrestrial planet regions of actively accreting protoplanetary disks. △ Less

Submitted 12 September, 2004; originally announced September 2004.

Comments: To appear in Astrophys. J., vol 614, Oct 10 2004 issue

Journal ref: Astrophys.J. 614 (2004) 490-496