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. 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. We find no evident departure from the trend at the extremes of the parameters explored. The scaling law can readily be incorporated into models of planet formation.

Chemical disequilibrium in exoplanetary atmospheres (detectable with remote spectroscopy) can indicate life. The modern Earth's atmosphere-ocean system has a much larger chemical disequilibrium than other solar system planets with atmospheres because of oxygenic photosynthesis. However, no analysis exists comparing disequilibrium on lifeless, prebiotic planets to disequilibrium on worlds with primitive chemotrophic biospheres that live off chemicals and not light. Here, we use a photochemical-microbial ecosystem model to calculate the atmosphere-ocean disequilibria of Earth with no life and with a chemotrophic biosphere. We show that the prebiotic Earth likely had a relatively large atmosphere-ocean disequilibrium due to the coexistence of water and volcanic H2, CO2, and CO. Subsequent chemotrophic life probably destroyed nearly all of the prebiotic disequilibrium through its metabolism, leaving a likely smaller disequilibrium between N2, CO2, CH4, and liquid water. So, disequilibrium fell with the rise of chemotrophic life then later rose with atmospheric oxygenation due to oxygenic photosynthesis. We conclude that big prebiotic disequilibrium between H2 and CO2 or CO and water is an anti-biosignature because these easily metabolized species can be eaten due to redox reactions with low activation energy barriers. However, large chemical disequilibrium can also be a biosignature when the disequilibrium arises from a chemical mixture with biologically insurmountable activation energy barriers, and clearly identifiable biogenic gases. Earth's modern disequilibrium between O2, N2, and liquid water along with minor CH4 is such a case. Thus, the interpretation of disequilibrium requires context. With context, disequilibrium can be used to infer dead or living worlds.

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.

Using Gaussian Process regression to analyze the Martian surface methane Tunable Laser Spectrometer (TLS) data reported by Webster (2018), we find that the TLS data, taken as a whole, are not statistically consistent with seasonal variability. The subset of data derived from an enrichment protocol of TLS, if considered in isolation, are equally consistent with either stochastic processes or periodic variability, but the latter does not favour seasonal variation.

Whether extant life exists in the martian subsurface is an open question. High concentrations of photochemically produced CO and H2 in the otherwise oxidizing martian atmosphere represent untapped sources of biologically useful free energy. These out-of-equilibrium species diffuse into the regolith, so subsurface microbes could use them as a source of energy and carbon. Indeed, CO oxidation and methanogenesis are relatively simple and evolutionarily ancient metabolisms on Earth. Consequently, assuming CO- or H2- consuming metabolisms would evolve on Mars, the persistence of CO and H2 in the martian atmosphere set limits on subsurface metabolic activity. Here, we constrain such maximum subsurface metabolic activity on Mars using a 1-D photochemical model with a hypothetical global biological sink on atmospheric CO and H2. We increase the biological sink until the model atmospheric composition diverges from observed abundances. We find maximum biological downward subsurface sinks of 1.5e8 molecules/cm^2/s for CO and 1.9e8 molecules/cm^2/s for H2. These covert to a maximum metabolizing biomass of <10^27 cells or <2e11 kg, equivalent to <1e-4 to 1e-5 of Earth's biomass, depending on the terrestrial estimate. Diffusion calculations suggest this upper biomass limit applies to the top few kilometers of the martian crust in communication with the atmosphere at low to mid latitudes. This biomass limit is more robust than previous estimates because we test multiple possible chemoautotrophic ecosystems over a broad parameter space of tunable model variables using an updated photochemical model with precise atmospheric concentrations and uncertainties from Curiosity. Our results of sparse or absent life in the martian subsurface also demonstrate how the atmospheric redox pairs of CO-O2 and H2-O2 may constitute antibiosignatures, which may be relevant to excluding life on exoplanets.

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.

The James Webb Space Telescope (JWST) may be capable of finding biogenic gases in the atmospheres of habitable exoplanets around low mass stars. Considerable attention has been given to the detectability of biogenic oxygen, which could be found using an ozone proxy, but ozone detection with JWST will be extremely challenging, even for the most favorable targets. Here, we investigate the detectability of biosignatures in anoxic atmospheres analogous to those that likely existed on the early Earth. Arguably, such anoxic biosignatures could be more prevalent than oxygen biosignatures if life exists elsewhere. Specifically, we simulate JWST retrievals of TRAPPIST-1e to determine whether the methane plus carbon dioxide disequilibrium biosignature pair is detectable in transit transmission. We find that ~10 transits using the Near InfraRed Spectrograph (NIRSpec) prism instrument may be sufficient to detect carbon dioxide and constrain methane abundances sufficiently well to rule out known, non-biological CH$_{4}$ production scenarios to ~90% confidence. Furthermore, it might be possible to put an upper limit on carbon monoxide abundances that would help rule out non-biological methane-production scenarios, assuming the surface biosphere would efficiently drawdown atmospheric CO. Our results are relatively insensitive to high altitude clouds and instrument noise floor assumptions, although stellar heterogeneity and variability may present challenges.

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.

Finding life on exoplanets from telescopic observations is an ultimate goal of exoplanet science. Life produces gases and other substances, such as pigments, which can have distinct spectral or photometric signatures. Whether or not life is found with future data must be expressed with probabilities, requiring a framework of biosignature assessment. We present a framework in which we advocate using biogeochemical "Exo-Earth System" models to simulate potential biosignatures in spectra or photometry. Given actual observations, simulations are used to find the Bayesian likelihoods of those data occurring for scenarios with and without life. The latter includes "false positives" where abiotic sources mimic biosignatures. Prior knowledge of factors influencing planetary inhabitation, including previous observations, is combined with the likelihoods to give the Bayesian posterior probability of life existing on a given exoplanet. Four components of observation and analysis are necessary. 1) Characterization of stellar (e.g., age and spectrum) and exoplanetary system properties, including "external" exoplanet parameters (e.g., mass and radius) to determine an exoplanet's suitability for life. 2) Characterization of "internal" exoplanet parameters (e.g., climate) to evaluate habitability. 3) Assessment of potential biosignatures within the environmental context (components 1-2) and any corroborating evidence. 4) Exclusion of false positives. The resulting posterior Bayesian probabilities of life's existence map to five confidence levels, ranging from "very likely" (90-100%) to "very unlikely" ($\le$10%) inhabited.

The early Earth's environment is controversial. Climatic estimates range from hot to glacial, and inferred marine pH spans strongly alkaline to acidic. Better understanding of early climate and ocean chemistry would improve our knowledge of the origin of life and its coevolution with the environment. Here, we use a geological carbon cycle model with ocean chemistry to calculate self-consistent histories of climate and ocean pH. Our carbon cycle model includes an empirically justified temperature and pH dependence of seafloor weathering, allowing the relative importance of continental and seafloor weathering to be evaluated. We find that the Archean climate was likely temperate (0-50 {\deg}C) due to the combined negative feedbacks of continental and seafloor weathering. Ocean pH evolves monotonically from 6.6 (+0.6,-0.4) (2{\sigma}) at 4.0 Ga to 7.0 (+0.7,-0.5) (2{\sigma}) at the Archean-Proterozoic boundary, and to 7.9 (+0.1,-0.2) (2{\sigma}) at the Proterozoic-Phanerozoic boundary. This evolution is driven by the secular decline of pCO2, which in turn is a consequence of increasing solar luminosity, but is moderated by carbonate alkalinity delivered from continental and seafloor weathering. Archean seafloor weathering may have been a comparable carbon sink to continental weathering, but is less dominant than previously assumed, and would not have induced global glaciation. We show how these conclusions are robust to a wide range of scenarios for continental growth, internal heat flow evolution and outgassing history, greenhouse gas abundances, and changes in the biotic enhancement of weathering.

How can scientists conclude with high confidence that an exoplanet hosts life? As telescopes come on line over the next 20 years that can directly observe photons from terrestrial exoplanets, this question will dictate the activities of many scientists across many fields. The expected data will be sparse and with low signal-to-noise, which will make disentangling biosignatures from abiotic features challenging. Our Earth is not just unique in that it hosts life, it is also the only terrestrial planet with direct observations of its interior through seismic waves, and compositional evolution through field and laboratory measurements. This extensive research reveals a planet born from collisions between worlds, followed by a complicated biogeochemical evolution. Exoplanet interiors, on the other hand, can only be constrained by the following observations: 1) photometric and spectroscopic analysis of the planet's atmosphere, 2) spectroscopic and photometric analysis of the host star, and 3) companion planet properties. From these (future) data, astrobiologists must generate plausible compositional and evolutionary models that constrain a potentially habitable exoplanet's internal properties and history, provide environmental context, and rule out geochemical explanations for any putative biosignatures. The goal of this white paper is to frame the role of geophysical and geochemical processes relevant to the search for life beyond the Solar System and to identify critical, but understudied, areas of future research.

Chemical disequilibrium in planetary atmospheres has been proposed as a generalized method for detecting life on exoplanets through remote spectroscopy. Among solar system planets with substantial atmospheres, the modern Earth has the largest thermodynamic chemical disequilibrium due to the presence of life. However, how this disequilibrium changed over time and, in particular, the biogenic disequilibria maintained in the anoxic Archean or less oxic Proterozoic eons are unknown. We calculate the atmosphere-ocean disequilibrium in the Precambrian using conservative proxy- and model-based estimates of early atmospheric and oceanic compositions. We omit crustal solids because subsurface composition is not detectable on exoplanets, unlike above-surface volatiles. We find that (i) disequilibrium increased through time in step with the rise of oxygen; (ii) both the Proterozoic and Phanerozoic may have had remotely detectable biogenic disequilibria due to the coexistence of O$_{2}$, N$_{2}$, and liquid water; and (iii) the Archean had a biogenic disequilibrium caused by the coexistence of N$_2$, CH$_4$, CO$_2$, and liquid water, which, for an exoplanet twin, may be remotely detectable. On the basis of this disequilibrium, we argue that the simultaneous detection of abundant CH$_{4}$ and CO$_{2}$ in a habitable exoplanet's atmosphere is a potential biosignature. Specifically, we show that methane mixing ratios greater than 0.001 are potentially biogenic, whereas those exceeding 0.01 are likely biogenic due to the difficulty in maintaining large abiotic methane fluxes to support high methane levels in anoxic atmospheres. Biogenicity would be strengthened by the absence of abundant CO, which should not coexist in a biological scenario.

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

Recent observations and analysis of low mass (<10$M_{\oplus}$), exoplanets have found that rocky planets only have radii up to 1.5-2$R_{\oplus}$. Two general hypotheses exist for the cause of the dichotomy between rocky and gas-enveloped planets (or possible water worlds): either low mass planets do not necessarily form thick atmospheres of a few wt. %, or the thick atmospheres on these planets easily escape driven by x-ray and extreme ultraviolet (XUV) emissions from young parent stars. Here we show that a cutoff between rocky and gas-enveloped planets due to hydrodynamic escape is most likely to occur at a mean radius of 1.76$\pm$0.38 (2$\sigma$) $R_{\oplus}$ around Sun-like stars. We examine the limit in rocky planet radii predicted by hydrodynamic escape across a wide range of possible model inputs using 10,000 parameter combinations drawn randomly from plausible parameter ranges. We find a cutoff between rocky and gas-enveloped planets that agrees with the observed cutoff. The large cross-section available for XUV absorption in the extremely distended primitive atmospheres of low mass planets results in complete loss of atmospheres during the ~100 Myr phase of stellar XUV saturation. In contrast, more massive planets have less distended atmospheres and less escape, and so retain thick atmospheres through XUV saturation and then indefinitely as the XUV and escape fluxes drop over time. The agreement between our model and exoplanet data leads us to conclude that hydrodynamic escape plausibly explains the observed upper limit on rocky planet size and few planets (a "valley" or "radius gap") in the 1.5-2$R_{\oplus}$ range.

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{\rho}$, 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.

Oxygen isotopes in marine cherts have been used to infer hot oceans during the Archean with temperatures between 60{\deg}C (333 K) and 80{\deg}C (353 K). Such climates are challenging for the early Earth warmed by the faint young Sun. The interpretation of the data has therefore been controversial. 1D climate modeling inferred that such hot climates would require very high levels of CO2 (2-6 bars). Previous carbon cycle modeling concluded that such stable hot climates were impossible and that the carbon cycle should lead to cold climates during the Hadean and the Archean. Here, we revisit the climate and carbon cycle of the early Earth at 3.8 Ga using a 3D climate-carbon model. We find that CO2 partial pressures of around 1 bar could have produced hot climates given a low land fraction and cloud feedback effects. However, such high CO2 partial pressures should not have been stable because of the weathering of terrestrial and oceanic basalts, producing an efficient stabilizing feedback. Moreover, the weathering of impact ejecta during the Late Heavy Bombardment (LHB) would have strongly reduced the CO2 partial pressure leading to cold climates and potentially snowball Earth events after large impacts. Our results therefore favor cold or temperate climates with global mean temperatures between around 8{\deg}C (281 K) and 30{\deg}C (303 K) and with 0.1-0.36 bar of CO2 for the late Hadean and early Archean. Finally, our model suggests that the carbon cycle was efficient for preserving clement conditions on the early Earth without necessarily requiring any other greenhouse gas or warming process.

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.

The next generation of ground and space-based telescopes will image habitable planets around nearby stars. A growing literature describes how to characterize such planets with spectroscopy, but less consideration has been given to the usefulness of planet colors. Here, we investigate whether potentially Earth-like exoplanets could be identified using UV-visible-to-NIR wavelength broadband photometry (350-1000 nm). Specifically, we calculate optimal photometric bins for identifying an exo-Earth and distinguishing it from uninhabitable planets including both Solar System objects and model exoplanets. The color of some hypothetical exoplanets - particularly icy terrestrial worlds with thick atmospheres - is similar to Earth's because of Rayleigh scattering in the blue region of the spectrum. Nevertheless, subtle features in Earth's reflectance spectrum appear to be unique. In particular, Earth's reflectance spectrum has a 'U-shape' unlike all our hypothetical, uninhabitable planets. This shape is partly biogenic because O2-rich, oxidizing air is transparent to sunlight, allowing prominent Rayleigh scattering, while ozone absorbs visible light, creating the bottom of the 'U'. Whether such uniqueness has practical utility depends on observational noise. If observations are photon limited or dominated by astrophysical sources (zodiacal light or imperfect starlight suppression), then the use of broadband visible wavelength photometry to identify Earth twins has little practical advantage over obtaining detailed spectra. However, if observations are dominated by dark current then optimized photometry could greatly assist preliminary characterization. We also calculate the optimal photometric bins for identifying extrasolar Archean Earths, and find that the Archean Earth is more difficult to unambiguously identify than a modern Earth twin.

Atmospheric chemical disequilibrium has been proposed as a method for detecting extraterrestrial biospheres from exoplanet observations. Chemical disequilibrium is potentially a generalized biosignature since it makes no assumptions about particular biogenic gases or metabolisms. Here, we present the first rigorous calculations of the thermodynamic chemical disequilibrium in Solar System atmospheres, in which we quantify the available Gibbs energy: the Gibbs free energy of an observed atmosphere minus that of atmospheric gases reacted to equilibrium. The purely gas phase disequilibrium in Earth's atmosphere is mostly attributable to O2 and CH4. The available Gibbs energy is not unusual compared to other Solar System atmospheres and smaller than that of Mars. However, Earth's fluid envelope contains an ocean, allowing gases to react with water and requiring a multiphase calculation with aqueous species. The disequilibrium in Earth's atmosphere-ocean system (in joules per mole of atmosphere) ranges from 20 to 2E6 times larger than the disequilibria of other atmospheres in the Solar System. Only on Earth is the chemical disequilibrium energy comparable to the thermal energy per mole of atmosphere. Earth's disequilibrium is biogenic, mainly caused by the coexistence of N2, O2 and liquid water instead of more stable nitrate. In comparison, the O2-CH4 disequilibrium is minor. We identify abiotic processes that cause disequilibrium in the other atmospheres. Our metric requires minimal assumptions and could potentially be calculated using observations of exoplanet atmospheres. However, further work is needed to establish whether thermodynamic disequilibrium is a practical exoplanet biosignature, requiring an assessment of false positives, noisy observations, and other detection challenges. Our Matlab code and databases for these calculations are available, open source.

A minimum atmospheric temperature, or tropopause, occurs at a pressure of around 0.1 bar in the atmospheres of Earth, Titan, Jupiter, Saturn, Uranus and Neptune, despite great differences in atmospheric composition, gravity, internal heat and sunlight. In all these bodies, the tropopause separates a stratosphere with a temperature profile that is controlled by the absorption of shortwave solar radiation, from a region below characterised by convection, weather, and clouds. However, it is not obvious why the tropopause occurs at the specific pressure near 0.1 bar. Here we use a physically-based model to demonstrate that, at atmospheric pressures lower than 0.1 bar, transparency to thermal radiation allows shortwave heating to dominate, creating a stratosphere. At higher pressures, atmospheres become opaque to thermal radiation, causing temperatures to increase with depth and convection to ensue. A common dependence of infrared opacity on pressure, arising from the shared physics of molecular absorption, sets the 0.1 bar tropopause. We hypothesize that a tropopause at a pressure of approximately 0.1 bar is characteristic of many thick atmospheres, including exoplanets and exomoons in our galaxy and beyond. Judicious use of this rule could help constrain the atmospheric structure, and thus the surface environments and habitability, of exoplanets.

We present an analytic 1-D radiative-convective model of the thermal structure of planetary atmospheres. Our model assumes that thermal radiative transfer is gray and can be represented by the two-stream approximation. Model atmospheres are assumed to be in hydrostatic equilibrium, with a power law scaling between the atmospheric pressure and the gray thermal optical depth. The convective portions of our models are taken to follow adiabats that account for condensation of volatiles through a scaling parameter to the dry adiabat. By combining these assumptions, we produce simple, analytic expressions that allow calculations of the atmospheric pressure-temperature profile, as well as expressions for the profiles of thermal radiative flux and convective flux. We explore the general behaviors of our model. These investigations encompass (1) worlds where atmospheric attenuation of sunlight is weak, which we show tend to have relatively high radiative-convective boundaries, (2) worlds with some attenuation of sunlight throughout the atmosphere, which we show can produce either shallow or deep radiative-convective boundaries, depending on the strength of sunlight attenuation, and (3) strongly irradiated giant planets (including Hot Jupiters), where we explore the conditions under which these worlds acquire detached convective regions in their mid-tropospheres. Finally, we validate our model and demonstrate its utility through comparisons to the average observed thermal structure of Venus, Jupiter, and Titan, and by comparing computed flux profiles to more complex models.