Showing 1–50 of 64 results for author: Abbot, D

Can AI weather models predict out-of-distribution gray swan tropical cyclones?

Authors: Y. Qiang Sun, Pedram Hassanzadeh, Mohsen Zand, Ashesh Chattopadhyay, Jonathan Weare, Dorian S. Abbot

Abstract: Predicting gray swan weather extremes, which are possible but so rare that they are absent from the training dataset, is a major concern for AI weather/climate models. An important open question is whether AI models can extrapolate from weaker weather events present in the training set to stronger, unseen weather extremes. To test this, we train independent versions of the AI model FourCastNet on… ▽ More Predicting gray swan weather extremes, which are possible but so rare that they are absent from the training dataset, is a major concern for AI weather/climate models. An important open question is whether AI models can extrapolate from weaker weather events present in the training set to stronger, unseen weather extremes. To test this, we train independent versions of the AI model FourCastNet on the 1979-2015 ERA5 dataset with all data, or with Category 3-5 tropical cyclones (TCs) removed, either globally or only over the North Atlantic or Western Pacific basin. We then test these versions of FourCastNet on 2018-2023 Category 5 TCs (gray swans). All versions yield similar accuracy for global weather, but the one trained without Category 3-5 TCs cannot accurately forecast Category 5 TCs, indicating that these models cannot extrapolate from weaker storms. The versions trained without Category 3-5 TCs in one basin show some skill forecasting Category 5 TCs in that basin, suggesting that FourCastNet can generalize across tropical basins. This is encouraging and surprising because regional information is implicitly encoded in inputs. No version satisfies gradient-wind balance, implying that enforcing such physical constraints may not improve generalizability to gray swans. Given that current state-of-the-art AI weather/climate models have similar learning strategies, we expect our findings to apply to other models and extreme events. Our work demonstrates that novel learning strategies are needed for AI weather/climate models to provide early warning or estimated statistics for the rarest, most impactful weather extremes. △ Less

Submitted 22 October, 2024; v1 submitted 18 October, 2024; originally announced October 2024.

Substantial extension of the lifetime of the terrestrial biosphere

Authors: R. J. Graham, Itay Halevy, Dorian Abbot

Abstract: Approximately one billion years (Gyr) in the future, as the Sun brightens, Earth's carbonate-silicate cycle is expected to drive CO$_2$ below the minimum level required by vascular land plants, eliminating most macroscopic land life. Here, we couple global-mean models of temperature- and CO$_2$-dependent plant productivity for C$_3$ and C$_4$ plants, silicate weathering, and climate to re-examine… ▽ More Approximately one billion years (Gyr) in the future, as the Sun brightens, Earth's carbonate-silicate cycle is expected to drive CO$_2$ below the minimum level required by vascular land plants, eliminating most macroscopic land life. Here, we couple global-mean models of temperature- and CO$_2$-dependent plant productivity for C$_3$ and C$_4$ plants, silicate weathering, and climate to re-examine the time remaining for terrestrial plants. If weathering is weakly temperature-dependent (as recent data suggest) and/or strongly CO$_2$-dependent, we find that the interplay between climate, productivity, and weathering causes the future luminosity-driven CO$_2$ decrease to slow and temporarily reverse, averting plant CO$_2$ starvation. This dramatically lengthens plant survival from 1 Gyr up to $\sim$1.6-1.86 Gyr, until extreme temperatures halt photosynthesis, suggesting a revised kill mechanism for land plants and potential doubling of the future lifespan of Earth's land macrobiota. An increased future lifespan for the complex biosphere may imply that Earth life had to achieve a smaller number of ``hard steps'' (unlikely evolutionary transitions) to produce intelligent life than previously estimated. These results also suggest that complex photosynthetic land life on Earth and exoplanets may be able to persist until the onset of the moist greenhouse transition. △ Less

Submitted 16 September, 2024; originally announced September 2024.

Comments: 19 pages, 6 figures. Accepted for publication in Planetary Science Journal

Climate Bistability at the Inner Edge of the Habitable Zone due to Runaway Greenhouse and Cloud Feedbacks

Authors: Bowen Fan, Da Yang, Dorian S. Abbot

Abstract: Understanding the climate dynamics at the inner edge of the habitable zone (HZ) is crucial for predicting the habitability of rocky exoplanets. Previous studies using Global Climate Models (GCMs) have indicated that planets receiving high stellar flux can exhibit climate bifurcations, leading to bistability between a cold (temperate) and a hot (runaway) climate. However, the mechanism causing this… ▽ More Understanding the climate dynamics at the inner edge of the habitable zone (HZ) is crucial for predicting the habitability of rocky exoplanets. Previous studies using Global Climate Models (GCMs) have indicated that planets receiving high stellar flux can exhibit climate bifurcations, leading to bistability between a cold (temperate) and a hot (runaway) climate. However, the mechanism causing this bistability has not been fully explained, in part due to the difficulty associated with inferring mechanisms from small numbers of expensive numerical simulations in GCMs. In this study, we employ a two-column (dayside and nightside), two-layer climate model to investigate the physical mechanisms driving this bistability. Through mechanism-denial experiments, we demonstrate that the runaway greenhouse effect, coupled with a cloud feedback on either the dayside or nightside, leads to climate bistability. We also map out the parameters that control the location of the bifurcations and size of the bistability. This work identifies which mechanisms and GCM parameters control the stellar flux at which rocky planets are likely to retain a hot, thick atmosphere if they experience a hot start. This is critical for the prioritization of targets and interpretation of observations by the James Webb Space Telescope (JWST). Furthermore, our modeling framework can be extended to planets with different condensable species and cloud types. △ Less

Submitted 22 August, 2024; originally announced August 2024.

Comments: 4 figures, submitted to ApJL

Report number: Accepted for publication in ApJL on 09.05.2024

Predator and Prey: A Minimum Recipe for the Transition from Steady to Oscillating Precipitation in Hothouse Climates

Authors: Da Yang, Dorian S. Abbot, Seth Seidel

Abstract: In the present tropical atmosphere, precipitation typically exhibits noisy, small-amplitude fluctuations about an average. However, recent cloud-resolving simulations show that in a hothouse climate, precipitation can shift to a regime characterized by nonlinear oscillations. In this regime, intense precipitation events are separated by several dry days. This raises questions about what triggers t… ▽ More In the present tropical atmosphere, precipitation typically exhibits noisy, small-amplitude fluctuations about an average. However, recent cloud-resolving simulations show that in a hothouse climate, precipitation can shift to a regime characterized by nonlinear oscillations. In this regime, intense precipitation events are separated by several dry days. This raises questions about what triggers the shift from a quasi-equilibrium state of precipitation to nonlinear precipitation oscillations and what factors determine the characteristics of these oscillations. To address these questions, we present a low-order model that includes two nonlinear ordinary differential equations, one for precipitation and the other for convective inhibition (CIN). Three processes govern the development of precipitation: a convective trigger that enhances precipitation, a self-limiting mechanism that reduces intense precipitation, and the effect of CIN in suppressing precipitation. On the other hand, CIN increases due to compensating subsidence from convection and decays exponentially over time due to radiation. A critical parameter in our model is the time-mean CIN (CIN*). As CIN* gradually increases, precipitation transitions from a quasi-equilibrium state to a nonlinear oscillation via a supercritical Hopf bifurcation. In the high CIN* limit, our model reduces to predator-prey dynamics, with CIN as the predator and precipitation as the prey. Here, the nonlinear oscillation's amplitude (maximum precipitation) grows with CIN*, and its period increases with the oscillation amplitude. A suite of cloud-resolving simulations corroborates these predictions from our low-order model, highlighting the effect of convective triggering and inhibition in regulating precipitation variability. △ Less

Submitted 21 August, 2024; originally announced August 2024.

Impact of Planetary Parameters on Water Clouds Microphysics

Authors: Huanzhou Yang, Thaddeus D. Komacek, Owen B. Toon, Eric T. Wolf, Tyler D. Robinson, Caroline Chael, Dorian S. Abbot

Abstract: Potentially habitable exoplanets are targets of great interest for the James Webb Space Telescope and upcoming mission concepts such as the Habitable Worlds Observatory. Clouds strongly affect climate and habitability, but predicting their properties is difficult. In Global Climate Models (GCMs), especially those aiming at simulating Earth, cloud microphysics is often crudely approximated by assum… ▽ More Potentially habitable exoplanets are targets of great interest for the James Webb Space Telescope and upcoming mission concepts such as the Habitable Worlds Observatory. Clouds strongly affect climate and habitability, but predicting their properties is difficult. In Global Climate Models (GCMs), especially those aiming at simulating Earth, cloud microphysics is often crudely approximated by assuming that all cloud particles have a single, constant size or a prescribed size distribution and that all clouds in a grid cell are identical. For exoplanets that range over a large phase space of planetary properties, this method could result in large errors. In this work, our goal is to determine how cloud microphysics on terrestrial exoplanets, whose condensable is mainly water vapor, depend on aerosol properties and planetary parameters such as surface pressure, surface gravity, and incident stellar radiation. We use the Community Aerosol and Radiation Model for Atmospheres as a 1D microphysical model to simulate the formation and evolution of clouds including the processes of nucleation, condensation, evaporation, coagulation, and vertical transfer. In these 1D idealized experiments, we find that the parameters that determine the macrophysical thermal structure of the atmospheres, including surface pressure and stellar flux, impact cloud radiative effect (CRE) most significantly. Parameters such as gravity and number density of aerosols working as cloud condensation nuclei affect the microphysical processes of cloud formation, including activation and vertical transfer. They also have a significant, though weaker effect on CRE. This work motivates the development of more accurate GCM cloud schemes and should aid in the interpretation of future observations. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 17 pages, 10 figures

Journal ref: Huanzhou Yang et al 2024 ApJ 966 152

Using Explainable AI and Transfer Learning to understand and predict the maintenance of Atlantic blocking with limited observational data

Authors: Huan Zhang, Justin Finkel, Dorian S. Abbot, Edwin P. Gerber, Jonathan Weare

Abstract: Blocking events are an important cause of extreme weather, especially long-lasting blocking events that trap weather systems in place. The duration of blocking events is, however, underestimated in climate models. Explainable Artificial Intelligence are a class of data analysis methods that can help identify physical causes of prolonged blocking events and diagnose model deficiencies. We demonstra… ▽ More Blocking events are an important cause of extreme weather, especially long-lasting blocking events that trap weather systems in place. The duration of blocking events is, however, underestimated in climate models. Explainable Artificial Intelligence are a class of data analysis methods that can help identify physical causes of prolonged blocking events and diagnose model deficiencies. We demonstrate this approach on an idealized quasigeostrophic model developed by Marshall and Molteni (1993). We train a convolutional neural network (CNN), and subsequently, build a sparse predictive model for the persistence of Atlantic blocking, conditioned on an initial high-pressure anomaly. Shapley Additive ExPlanation (SHAP) analysis reveals that high-pressure anomalies in the American Southeast and North Atlantic, separated by a trough over Atlantic Canada, contribute significantly to prediction of sustained blocking events in the Atlantic region. This agrees with previous work that identified precursors in the same regions via wave train analysis. When we apply the same CNN to blockings in the ERA5 atmospheric reanalysis, there is insufficient data to accurately predict persistent blocks. We partially overcome this limitation by pre-training the CNN on the plentiful data of the Marshall-Molteni model, and then using Transfer Learning to achieve better predictions than direct training. SHAP analysis before and after transfer learning allows a comparison between the predictive features in the reanalysis and the quasigeostrophic model, quantifying dynamical biases in the idealized model. This work demonstrates the potential for machine learning methods to extract meaningful precursors of extreme weather events and achieve better prediction using limited observational data. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: 29 pages, 10 figures

AI can identify Solar System instability billions of years in advance

Authors: Dorian S. Abbot, J. D. Laurence-Chasen, Robert J. Webber, David M. Hernandez, Jonathan Weare

Abstract: Rare event schemes require an approximation of the probability of the rare event as a function of system state. Finding an appropriate reaction coordinate is typically the most challenging aspect of applying a rare event scheme. Here we develop an artificial intelligence (AI) based reaction coordinate that effectively predicts which of a limited number of simulations of the Solar System will go un… ▽ More Rare event schemes require an approximation of the probability of the rare event as a function of system state. Finding an appropriate reaction coordinate is typically the most challenging aspect of applying a rare event scheme. Here we develop an artificial intelligence (AI) based reaction coordinate that effectively predicts which of a limited number of simulations of the Solar System will go unstable using a convolutional neural network classifier. The performance of the algorithm does not degrade significantly even 3.5 billion years before the instability. We overcome the class imbalance intrinsic to rare event problems using a combination of minority class oversampling, increased minority class weighting, and pulling multiple non-overlapping training sequences from simulations. Our success suggests that AI may provide a promising avenue for developing reaction coordinates without detailed theoretical knowledge of the system. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: 3 pages, 1 figure

Journal ref: Res. Notes AAS 8 3 (2024)

Critical role of vertical radiative cooling contrast in triggering episodic deluges in small-domain hothouse climates

Authors: Xinyi Song, Dorian S. Abbot, Jun Yang

Abstract: Seeley and Wordsworth (2021) showed that in small-domain cloud-resolving simulations the pattern of precipitation transforms in extremely hot climates ($\ge$ 320 K) from quasi-steady to organized episodic deluges, with outbursts of heavy rain alternating with several dry days. They proposed a mechanism for this transition involving increased water vapor absorption of solar radiation leading to net… ▽ More Seeley and Wordsworth (2021) showed that in small-domain cloud-resolving simulations the pattern of precipitation transforms in extremely hot climates ($\ge$ 320 K) from quasi-steady to organized episodic deluges, with outbursts of heavy rain alternating with several dry days. They proposed a mechanism for this transition involving increased water vapor absorption of solar radiation leading to net lower-tropospheric radiative heating. This heating inhibits lower-tropospheric convection and decouples the boundary layer from the upper troposphere during the dry phase, allowing lower-tropospheric moist static energy to build until it discharges, resulting in a deluge. We perform cloud-resolving simulations in polar night and show that the same transition occurs, implying that some revision of their mechanism is necessary. We show that episodic deluges can occur even if the lower-tropospheric radiative heating rate is negative, as long as the magnitude of the upper-tropospheric radiative cooling is about twice as large. We find that in the episodic deluge regime the mean precipitation can be inferred from the atmospheric column energy budget and the period can be predicted from the time for radiation and reevaporation to cool the lower atmosphere. △ Less

Submitted 2 July, 2023; originally announced July 2023.

Mercury's chaotic secular evolution as a subdiffusive process

Authors: Dorian S. Abbot, Robert J. Webber, David M. Hernandez, Sam Hadden, Jonathan Weare

Abstract: Mercury's orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can cause g1 to decrease to the approximately constant value of g5 and create a resonance. Previous work has approximated the variation in g1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular and N-body mod… ▽ More Mercury's orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can cause g1 to decrease to the approximately constant value of g5 and create a resonance. Previous work has approximated the variation in g1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular and N-body models on timescales longer than 10 Gyr. Here we show that the diffusive model underpredicts the Mercury instability probability by a factor of 3-10,000 on timescales less than 5 Gyr, the remaining lifespan of the Solar System. This is because g1 exhibits larger variations on short timescales than the diffusive model would suggest. To better model the variations on short timescales, we build a new subdiffusive phenomological model for g1. Subdiffusion is similar to diffusion but exhibits larger displacements on short timescales and smaller displacements on long timescales. We choose model parameters based on the behavior of the g1 trajectories in the N-body simulations, leading to a tuned model that can reproduce Mercury instability statistics from 1-40 Gyr. This work motivates fundamental questions in Solar System dynamics: Why does subdiffusion better approximate the variation in g1 than standard diffusion? Why is there an upper bound on g1, but not a lower bound that would prevent it from reaching g5? △ Less

Submitted 12 April, 2024; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: accepted at ApJ

Simple physics and integrators accurately reproduce Mercury instability statistics

Authors: Dorian S. Abbot, David M. Hernandez, Sam Hadden, Robert J. Webber, Georgios P. Afentakis, Jonathan Weare

Abstract: The long-term stability of the Solar System is an issue of significant scientific and philosophical interest. The mechanism leading to instability is Mercury's eccentricity being pumped up so high that Mercury either collides with Venus or is scattered into the Sun. Previously, only three five-billion-year $N$-body ensembles of the Solar System with thousands of simulations have been run to assess… ▽ More The long-term stability of the Solar System is an issue of significant scientific and philosophical interest. The mechanism leading to instability is Mercury's eccentricity being pumped up so high that Mercury either collides with Venus or is scattered into the Sun. Previously, only three five-billion-year $N$-body ensembles of the Solar System with thousands of simulations have been run to assess long-term stability. We generate two additional ensembles, each with 2750 members, and make them publicly available at \texttt{https://archive.org/details/@dorianabbot}. We find that accurate Mercury instability statistics can be obtained by (1) including only the Sun and the 8 planets, (2) using a simple Wisdom-Holman scheme without correctors, (3) using a basic representation of general relativity, and (4) using a time step of 3.16 days. By combining our Solar System ensembles with previous ensembles we form a 9,601-member ensemble of ensembles. In this ensemble of ensembles, the logarithm of the frequency of a Mercury instability event increases linearly with time between 1.3 and 5 Gyr, suggesting that a single mechanism is responsible for Mercury instabilities in this time range and that this mechanism becomes more active as time progresses. Our work provides a robust estimate of Mercury instability statistics over the next five billion years, outlines methodologies that may be useful for exoplanet system investigations, and provides two large ensembles of publicly available Solar System integrations that can serve as testbeds for theoretical ideas as well as training sets for artificial intelligence schemes. △ Less

Submitted 21 February, 2023; v1 submitted 30 December, 2022; originally announced December 2022.

Comments: accepted at ApJ

Journal ref: Ap.J. 944 (2023) 190

Inner Habitable Zone Boundary for Eccentric Exoplanets

Authors: Xuan Ji, Nora Bailey, Daniel Fabrycky, Edwin S. Kite, Jonathan H. Jiang, Dorian S. Abbot

Abstract: The climate of a planet can be strongly affected by its eccentricity due to variations in the stellar flux. There are two limits for the dependence of the inner habitable zone boundary (IHZ) on eccentricity: (1) the mean-stellar flux approximation ($S_{\mbox{IHZ}} \propto \sqrt{1-e^2}$), in which the temperature is approximately constant throughout the orbit, and (2) the maximum-stellar flux appro… ▽ More The climate of a planet can be strongly affected by its eccentricity due to variations in the stellar flux. There are two limits for the dependence of the inner habitable zone boundary (IHZ) on eccentricity: (1) the mean-stellar flux approximation ($S_{\mbox{IHZ}} \propto \sqrt{1-e^2}$), in which the temperature is approximately constant throughout the orbit, and (2) the maximum-stellar flux approximation ($S_{\mbox{IHZ}} \propto (1-e)^2$), in which the temperature adjusts instantaneously to the stellar flux. Which limit is appropriate is determined by the dimensionless parameter $Π= \frac{C}{BP}$, where $C$ is the heat capacity of the planet, $P$ is the orbital period, and $B=\frac{\partial Ω}{\partial T_s}$, where $Ω$ is the outgoing longwave radiation and $T_s$ is the surface temperature. We use the Buckingham $Π$ theorem to derive an analytical function for the IHZ in terms of eccentricity and $Π$. We then build a time-dependent energy balance model to resolve the surface temperature evolution and constrain our analytical result. We find that $Π$ must be greater than about $\sim 1$ for the mean-stellar flux approximation to be nearly exact and less than about $\sim 0.01$ for the maximum-stellar flux approximation to be nearly exact. In addition to assuming a constant heat capacity, we also consider the effective heat capacity including latent heat (evaporation and precipitation). We find that for planets with an Earth-like ocean, the IHZ should follow the mean-stellar flux limit for all eccentricities. This work will aid in the prioritization of potentially habitable exoplanets with non-zero eccentricity for follow-up characterization. △ Less

Submitted 11 January, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted to ApJL

Revealing the statistics of extreme events hidden in short weather forecast data

Authors: Justin Finkel, Edwin P. Gerber, Dorian S. Abbot, Jonathan Weare

Abstract: Extreme weather events have significant consequences, dominating the impact of climate on society. While high-resolution weather models can forecast many types of extreme events on synoptic timescales, long-term climatological risk assessment is an altogether different problem. A once-in-a-century event takes, on average, 100 years of simulation time to appear just once, far beyond the typical int… ▽ More Extreme weather events have significant consequences, dominating the impact of climate on society. While high-resolution weather models can forecast many types of extreme events on synoptic timescales, long-term climatological risk assessment is an altogether different problem. A once-in-a-century event takes, on average, 100 years of simulation time to appear just once, far beyond the typical integration length of a weather forecast model. Therefore, this task is left to cheaper, but less accurate, low-resolution or statistical models. But there is untapped potential in weather model output: despite being short in duration, weather forecast ensembles are produced multiple times a week. Integrations are launched with independent perturbations, causing them to spread apart over time and broadly sample phase space. Collectively, these integrations add up to thousands of years of data. We establish methods to extract climatological information from these short weather simulations. Using ensemble hindcasts by the European Center for Medium-range Weather Forecasting (ECMWF) archived in the subseasonal-to-seasonal (S2S) database, we characterize sudden stratospheric warming (SSW) events with multi-centennial return times. Consistent results are found between alternative methods, including basic counting strategies and Markov state modeling. By carefully combining trajectories together, we obtain estimates of SSW frequencies and their seasonal distributions that are consistent with reanalysis-derived estimates for moderately rare events, but with much tighter uncertainty bounds, and which can be extended to events of unprecedented severity that have not yet been observed historically. These methods hold potential for assessing extreme events throughout the climate system, beyond this example of stratospheric extremes. △ Less

Submitted 23 January, 2023; v1 submitted 10 June, 2022; originally announced June 2022.

Comments: 28 pages, 6 figures. Resubmitted for publication

The Effect of Ocean Salinity on Climate and Its Implications for Earth's Habitability

Authors: Stephanie L. Olson, Malte F. Jansen, Dorian S. Abbot, Itay Halevy, Colin Goldblatt

Abstract: The influence of atmospheric composition on the climates of present-day and early Earth has been studied extensively, but the role of ocean composition has received less attention. We use the ROCKE-3D ocean-atmosphere general circulation model to investigate the response of Earth's present-day and Archean climate system to low vs. high ocean salinity. We find that saltier oceans yield warmer clima… ▽ More The influence of atmospheric composition on the climates of present-day and early Earth has been studied extensively, but the role of ocean composition has received less attention. We use the ROCKE-3D ocean-atmosphere general circulation model to investigate the response of Earth's present-day and Archean climate system to low vs. high ocean salinity. We find that saltier oceans yield warmer climates in large part due to changes in ocean dynamics. Increasing ocean salinity from 20 g/kg to 50 g/kg results in a 71% reduction in sea ice cover in our present-day Earth scenario. This same salinity change also halves the pCO$_2$ threshold at which Snowball glaciation occurs in our Archean scenarios. In combination with higher levels of greenhouse gases such as CO$_2$ and CH$_4$, a saltier ocean may allow for a warm Archean Earth with only seasonal ice at the poles despite receiving 20% less energy from the Sun. △ Less

Submitted 13 May, 2022; originally announced May 2022.

Comments: Accepted for publication in Geophysical Research Letters 04/19/22

Critical percolation threshold restricts late-summer Arctic sea ice melt pond coverage

Authors: Predrag Popović, Mary C. Silber, Dorian S. Abbot

Abstract: During the summer, vast regions of Arctic sea ice are covered by meltwater ponds that significantly lower the ice reflectivity and accelerate melting. Ponds develop over the melt season through an initial rapid growth stage followed by drainage through macroscopic holes. Recent analysis of melt pond photographs indicates that late-summer ponds exist near the critical percolation threshold, a speci… ▽ More During the summer, vast regions of Arctic sea ice are covered by meltwater ponds that significantly lower the ice reflectivity and accelerate melting. Ponds develop over the melt season through an initial rapid growth stage followed by drainage through macroscopic holes. Recent analysis of melt pond photographs indicates that late-summer ponds exist near the critical percolation threshold, a special pond coverage fraction below which the ponds are largely disconnected and above which they are highly connected. Here, we show that the percolation threshold, a statistical property of ice topography, constrains pond evolution due to pond drainage through macroscopic holes. We show that it sets the approximate upper limit and scales pond coverage throughout its evolution after the beginning of drainage. Furthermore, we show that the rescaled pond coverage during drainage is a universal function of a single non-dimensional parameter, $η$, roughly interpreted as the number of drainage holes per characteristic area of the surface. This universal curve allows us to formulate an equation for pond coverage time-evolution during and after pond drainage that captures the dependence on environmental parameters and is supported by observations on undeformed first-year ice. This equation reveals that pond coverage is highly sensitive to environmental parameters, suggesting that modeling uncertainties could be reduced by more directly parameterizing the ponds' natural parameter, $η$. Our work uncovers previously unrecognized constraints on melt pond physics and places ponds within a broader context of phase transitions and critical phenomena. Therefore, it holds promise for improving ice-albedo parameterizations in large-scale models. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Journal ref: Popović, P., Silber, M. C., and Abbot, D. S. (2020). Critical percolation threshold restricts late-summer Arctic sea ice melt pond coverage. Journal of Geophysical Research: Oceans, 125(8), e2019JC016029

Snow topography on undeformed Arctic sea ice captured by an idealized "snow dune" model

Authors: Predrag Popović, Justin Finkel, Mary C. Silber, Dorian S. Abbot

Abstract: Our ability to predict the future of Arctic sea ice is limited by ice's sensitivity to detailed surface conditions such as the distribution of snow and melt ponds. Snow on top of the ice decreases ice's thermal conductivity, increases its reflectivity (albedo), and provides a source of meltwater for melt ponds during summer that decrease the ice's albedo. In this paper, we develop a simple model o… ▽ More Our ability to predict the future of Arctic sea ice is limited by ice's sensitivity to detailed surface conditions such as the distribution of snow and melt ponds. Snow on top of the ice decreases ice's thermal conductivity, increases its reflectivity (albedo), and provides a source of meltwater for melt ponds during summer that decrease the ice's albedo. In this paper, we develop a simple model of pre-melt snow topography that accurately describes snow cover of flat, undeformed Arctic sea ice on several study sites for which data was available. The model considers a surface that is a sum of randomly sized and placed "snow dunes" represented as Gaussian mounds. This model generalizes the "void model" of Popović et al. (2018) and, as such, accurately describes the statistics of melt pond geometry. We test this model against detailed LiDAR measurements of the pre-melt snow topography. We show that the model snow-depth distribution is statistically indistinguishable from the measurements on flat ice, while small disagreement exists if the ice is deformed. We then use this model to determine analytic expressions for the conductive heat flux through the ice and for melt pond coverage evolution during an early stage of pond formation. We also formulate a criterion for ice to remain pond-free throughout the summer. Results from our model could be directly included in large-scale models, thereby improving our understanding of energy balance on sea ice and allowing for more reliable predictions of Arctic sea ice in a future climate. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Journal ref: Popović, P., Finkel, J., Silber, M. C., & Abbot, D. S. (2020). Snow topography on undeformed Arctic sea ice captured by an idealized "snow dune" model. Journal of Geophysical Research: Oceans, 125(9), e2019JC016034

Data-driven transition path analysis yields a statistical understanding of sudden stratospheric warming events in an idealized model

Authors: Justin Finkel, Robert J. Webber, Edwin P. Gerber, Dorian S. Abbot, Jonathan Weare

Abstract: Atmospheric regime transitions are highly impactful as drivers of extreme weather events, but pose two formidable modeling challenges: predicting the next event (weather forecasting), and characterizing the statistics of events of a given severity (the risk climatology). Each event has a different duration and spatial structure, making it hard to define an objective "average event." We argue here… ▽ More Atmospheric regime transitions are highly impactful as drivers of extreme weather events, but pose two formidable modeling challenges: predicting the next event (weather forecasting), and characterizing the statistics of events of a given severity (the risk climatology). Each event has a different duration and spatial structure, making it hard to define an objective "average event." We argue here that transition path theory (TPT), a stochastic process framework, is an appropriate tool for the task. We demonstrate TPT's capacities on a wave-mean flow model of sudden stratospheric warmings (SSWs) developed by Holton and Mass (1976), which is idealized enough for transparent TPT analysis but complex enough to demonstrate computational scalability. Whereas a recent article (Finkel et al. 2021) studied near-term SSW predictability, the present article uses TPT to link predictability to long-term SSW frequency. This requires not only forecasting forward in time from an initial condition, but also \emph{backward in time} to assess the probability of the initial conditions themselves. TPT enables one to condition the dynamics on the regime transition occurring, and thus visualize its physical drivers with a vector field called the \emph{reactive current}. The reactive current shows that before an SSW, dissipation and stochastic forcing drive a slow decay of vortex strength at lower altitudes. The response of upper-level winds is late and sudden, occurring only after the transition is almost complete from a probabilistic point of view. This case study demonstrates that TPT quantities, visualized in a space of physically meaningful variables, can help one understand the dynamics of regime transitions. △ Less

Submitted 19 October, 2022; v1 submitted 28 August, 2021; originally announced August 2021.

Comments: 18 pages, 7 figures (main text), 19 pages, 1 figure (supplement). Accepted for publication in the Journal of the Atmospheric Sciences

Rare Event Sampling Improves Mercury Instability Statistics

Authors: Dorian S. Abbot, Robert J. Webber, Sam Hadden, Darryl Seligman, Jonathan Weare

Abstract: Due to the chaotic nature of planetary dynamics, there is a non-zero probability that Mercury's orbit will become unstable in the future. Previous efforts have estimated the probability of this happening between 3 and 5 billion years in the future using a large number of direct numerical simulations with an N-body code, but were not able to obtain accurate estimates before 3 billion years in the f… ▽ More Due to the chaotic nature of planetary dynamics, there is a non-zero probability that Mercury's orbit will become unstable in the future. Previous efforts have estimated the probability of this happening between 3 and 5 billion years in the future using a large number of direct numerical simulations with an N-body code, but were not able to obtain accurate estimates before 3 billion years in the future because Mercury instability events are too rare. In this paper we use a new rare event sampling technique, Quantile Diffusion Monte Carlo (QDMC), to estimate that the probability of a Mercury instability event in the next 2 billion years is approximately $10^{-4}$ in the REBOUND N-body code. We show that QDMC provides unbiased probability estimates at a computational cost of up to 100 times less than direct numerical simulation. QDMC is easy to implement and could be applied to many problems in planetary dynamics in which it is necessary to estimate the probability of a rare event. △ Less

Submitted 27 December, 2021; v1 submitted 16 June, 2021; originally announced June 2021.

Comments: Abbot et al 2021 ApJ 923 236

Learning forecasts of rare stratospheric transitions from short simulations

Authors: Justin Finkel, Robert J. Webber, Dorian S. Abbot, Edwin P. Gerber, Jonathan Weare

Abstract: Rare events arising in nonlinear atmospheric dynamics remain hard to predict and attribute. We address the problem of forecasting rare events in a prototypical example, Sudden Stratospheric Warmings (SSWs). Approximately once every other winter, the boreal stratospheric polar vortex rapidly breaks down, shifting midlatitude surface weather patterns for months. We focus on two key quantities of int… ▽ More Rare events arising in nonlinear atmospheric dynamics remain hard to predict and attribute. We address the problem of forecasting rare events in a prototypical example, Sudden Stratospheric Warmings (SSWs). Approximately once every other winter, the boreal stratospheric polar vortex rapidly breaks down, shifting midlatitude surface weather patterns for months. We focus on two key quantities of interest: the probability of an SSW occurring, and the expected lead time if it does occur, as functions of initial condition. These \emph{optimal forecasts} concretely measure the event's progress. Direct numerical simulation can estimate them in principle, but is prohibitively expensive in practice: each rare event requires a long integration to observe, and the cost of each integration grows with model complexity. We describe an alternative approach using integrations that are \emph{short} compared to the timescale of the warming event. We compute the probability and lead time efficiently by solving equations involving the transition operator, which encodes all information about the dynamics. We relate these optimal forecasts to a small number of interpretable physical variables, suggesting optimal measurements for forecasting. We illustrate the methodology on a prototype SSW model developed by Holton and Mass (1976) and modified by stochastic forcing. While highly idealized, this model captures the essential nonlinear dynamics of SSWs and exhibits the key forecasting challenge: the dramatic separation in timescales between a single event and the return time between successive events. Our methodology is designed to fully exploit high-dimensional data from models and observations, and has the potential to identify detailed predictors of many complex rare events in meteorology. △ Less

Submitted 28 August, 2021; v1 submitted 15 February, 2021; originally announced February 2021.

Comments: 26 pages, 7 figures, major revision after original. Accepted to Monthly Weather Review, American Meteorological Society

Probing the capability of future direct imaging missions to spectrally constrain the frequency of Earth-like planets

Authors: Jade H. Checlair, Geronimo L. Villanueva, Benjamin P. C. Hayworth, Stephanie L. Olson, Thaddeus D. Komacek, Tyler D. Robinson, Predrag Popovic, Huanzhou Yang, Dorian S. Abbot

Abstract: A critical question in astrobiology is whether exoEarth candidates (EECs) are Earth-like, in that they originate life that progressively oxygenates their atmospheres similarly to Earth. We propose answering this question statistically by searching for O2 and O3 on EECs with missions such as HabEx or LUVOIR. We explore the ability of these missions to constrain the fraction, fE, of EECs that are Ea… ▽ More A critical question in astrobiology is whether exoEarth candidates (EECs) are Earth-like, in that they originate life that progressively oxygenates their atmospheres similarly to Earth. We propose answering this question statistically by searching for O2 and O3 on EECs with missions such as HabEx or LUVOIR. We explore the ability of these missions to constrain the fraction, fE, of EECs that are Earth-like in the event of a null detection of O2 or O3 on all observed EECs. We use the Planetary Spectrum Generator to simulate observations of EECs with O2 and O3 levels based on Earth's history. We consider four instrument designs: LUVOIR-A (15m), LUVOIR-B (8m), HabEx with a starshade (4m, "HabEx/SS"), HabEx without a starshade (4m, "HabEx/no-SS"); as well as three estimates of the occurrence rate of EECs (eta_earth): 24%, 5%, and 0.5%. In the case of a null-detection, we find that for eta_earth = 24%, LUVOIR-A, LUVOIR-B, and HabEx/SS would constrain fE to <= 0.094, <= 0.18, and <= 0.56, respectively. This also indicates that if fE is greater than these upper limits, we are likely to detect O3 on at least 1 EEC. Conversely, we find that HabEx/no-SS cannot constrain fE, due to the lack of an coronagraph ultraviolet channel. For eta_earth = 5%, only LUVOIR-A and LUVOIR-B would be able to constrain fE, to <= 0.45 and <= 0.85, respectively. For eta_earth = 0.5%, none of the missions would allow us to constrain fE, due to the low number of detectable EECs. We conclude that the ability to constrain fE is more robust to uncertainties in eta_earth for missions with larger aperture mirrors. However all missions are susceptible to an inconclusive null detection if eta_earth is sufficiently low. △ Less

Submitted 18 January, 2021; originally announced January 2021.

Comments: Accepted at The Astronomical Journal, January 7 2021. arXiv admin note: substantial text overlap with arXiv:2008.03952

Non-detection of O$_2$/O$_3$ informs frequency of Earth-like planets with LUVOIR but not HabEx

Authors: Jade H. Checlair, Benjamin P. C. Hayworth, Stephanie L. Olson, Thaddeus D. Komacek, Geronimo L. Villanueva, Predrag Popović, Huanzhou Yang, Dorian S. Abbot

Abstract: A critical question in the search for extraterrestrial life is whether exoEarths are Earth-like, in that they host life that progressively oxygenates their atmospheres roughly following Earth's oxygenation history. This question could be answered statistically by searching for O$_2$ and O$_3$ on exoEarths detected by HabEx or LUVOIR. The point of this paper is to compare the ability of HabEx and L… ▽ More A critical question in the search for extraterrestrial life is whether exoEarths are Earth-like, in that they host life that progressively oxygenates their atmospheres roughly following Earth's oxygenation history. This question could be answered statistically by searching for O$_2$ and O$_3$ on exoEarths detected by HabEx or LUVOIR. The point of this paper is to compare the ability of HabEx and LUVOIR to prevent a false negative answer to this question, in which we do not detect O$_2$ or O$_3$ on any planet even if all exoEarths are Earth-like. Our approach is to assign O$_2$ and O$_3$ values drawn from Earth's history to a distribution of detectable exoEarths and determine whether O$_2$ and O$_3$ would be detectable using the Planetary Spectrum Generator. We find that if exoEarths tend to be Earth-like, we expect to detect O$_3$ with a LUVOIR-sized instrument. We also find that LUVOIR is unlikely to have a false negative scenario in the context of searching for Earth-like life on its targeted exoEarths. Because of that, if LUVOIR does not detect O$_2$ or O$_3$ on any exoEarths, we will be able to constrain the maximum number of exoEarths that could be Earth-like. In contrast, we find that even if all exoEarths are Earth-like, HabEx has up to a 22% chance of not detecting O$_2$ or O$_3$ on any of them. This is because HabEx will detect less planets and cannot reliably detect O$_2$ and O$_3$ at all potential Proterozoic levels. This is a strong argument for building a larger telescope such as LUVOIR if we want to determine whether exoEarths tend to be Earth-like. △ Less

Submitted 16 August, 2020; v1 submitted 10 August, 2020; originally announced August 2020.

Comments: Must revise noise parameters used for HabEx

Hurricane genesis is favorable on terrestrial exoplanets orbiting late-type M dwarf stars

Authors: Thaddeus D. Komacek, Daniel R. Chavas, Dorian S. Abbot

Abstract: Hurricanes are one of the most extreme storm systems that occur on Earth, characterized by strong rainfall and fast winds. The terrestrial exoplanets that will be characterized with future infrared space telescopes orbit M dwarf stars. As a result, the best observable terrestrial exoplanets have vastly different climates than Earth, with a large dayside-to-nightside irradiation contrast and relati… ▽ More Hurricanes are one of the most extreme storm systems that occur on Earth, characterized by strong rainfall and fast winds. The terrestrial exoplanets that will be characterized with future infrared space telescopes orbit M dwarf stars. As a result, the best observable terrestrial exoplanets have vastly different climates than Earth, with a large dayside-to-nightside irradiation contrast and relatively slow rotation. Hurricanes may affect future observations of terrestrial exoplanets because they enhance the vertical transport of water vapor and could influence ocean heat transport. In this work, we explore how the environment of terrestrial exoplanets orbiting M dwarf stars affects the favorability of hurricane genesis (formation). To do so, we apply metrics developed to understand hurricane genesis on Earth to three-dimensional climate models of ocean-covered exoplanets orbiting M dwarf stars. We find that hurricane genesis is most favorable on intermediate-rotating tidally locked terrestrial exoplanets with rotation periods of $\sim 8-10~\mathrm{days}$. As a result, hurricane genesis is most favorable for terrestrial exoplanets in the habitable zones of late-type M dwarf stars. The peak in the favorability of hurricane genesis at intermediate rotation occurs because sufficient spin is required for hurricane genesis, but the vertical wind shear on fast-rotating terrestrial exoplanets disrupts hurricane genesis. We find that hurricane genesis is less favorable on slowly rotating terrestrial exoplanets, which agrees with previous work. Future work using simulations that resolve hurricane genesis and evolution can test our expectations for how the environment affects the favorability of hurricane genesis on tidally locked terrestrial exoplanets. △ Less

Submitted 26 June, 2020; originally announced June 2020.

Comments: Accepted at ApJ, 13 pages, 5 figures, 2 tables

Journal ref: The Astrophysical Journal, 898:115, 2020 August 1

The Effect of Substellar Continent Size on Ocean Dynamics of Proxima Centauri b

Authors: Andrea M. Salazar, Stephanie L. Olson, Thaddeus D. Komacek, Haynes Stephens, Dorian S. Abbot

Abstract: The potential habitability of tidally locked planets orbiting M-dwarf stars has been widely investigated in recent work, typically with a non-dynamic ocean and without continents. On Earth, ocean dynamics are a primary means of heat and nutrient distribution. Continents are a critical source of nutrients, strongly influence ocean dynamics, and participate in climate regulation. In this work, we in… ▽ More The potential habitability of tidally locked planets orbiting M-dwarf stars has been widely investigated in recent work, typically with a non-dynamic ocean and without continents. On Earth, ocean dynamics are a primary means of heat and nutrient distribution. Continents are a critical source of nutrients, strongly influence ocean dynamics, and participate in climate regulation. In this work, we investigate how the size of a substellar land mass affects the oceans ability to transport heat and upwell nutrients on the tidally locked planet Proxima Centauri b using the ROCKE-3D coupled ocean-atmosphere General Circulation Model (GCM). We find that dayside ice-free ocean and nutrient delivery to the mixed layer via upwelling are maintained across all continent sizes. We also find that Proxima Centauri bs climate is more sensitive to differences among atmospheric GCMs than to the inclusion of ocean dynamics in ROCKE-3D. Finally, we find that Proxima Centauri b transitions from a lobster state where ocean heat transport distributes heat away from the substellar point to an eyeball state where heat transport is restricted and surface temperature decreases symmetrically from the substellar point when the continent size exceeds about 20 percent of the surface area. Our work suggests that both a dynamic ocean and continents are unlikely to decrease the habitability prospects of nearby tidally locked targets like Proxima Centauri b that could be investigated with future observations by the James Webb Space Telescope (JWST). △ Less

Submitted 28 May, 2020; originally announced May 2020.

Comments: Accepted to ApJ Letters May 19th, 2020

Path properties of atmospheric transitions: illustration with a low-order sudden stratospheric warming model

Authors: Justin Finkel, Dorian Abbot, Jonathan Weare

Abstract: Many rare weather events, including hurricanes, droughts, and floods, dramatically impact human life. To accurately forecast these events and characterize their climatology requires specialized mathematical techniques to fully leverage the limited data that are available. Here we describe \emph{transition path theory} (TPT), a framework originally developed for molecular simulation, and argue that… ▽ More Many rare weather events, including hurricanes, droughts, and floods, dramatically impact human life. To accurately forecast these events and characterize their climatology requires specialized mathematical techniques to fully leverage the limited data that are available. Here we describe \emph{transition path theory} (TPT), a framework originally developed for molecular simulation, and argue that it is a useful paradigm for developing mechanistic understanding of rare climate events. TPT provides a method to calculate statistical properties of the paths into the event. As an initial demonstration of the utility of TPT, we analyze a low-order model of sudden stratospheric warming (SSW), a dramatic disturbance to the polar vortex which can induce extreme cold spells at the surface in the midlatitudes. SSW events pose a major challenge for seasonal weather prediction because of their rapid, complex onset and development. Climate models struggle to capture the long-term statistics of SSW, owing to their diversity and intermittent nature. We use a stochastically forced Holton-Mass-type model with two stable states, corresponding to radiative equilibrium and a vacillating SSW-like regime. In this stochastic bistable setting, from certain probabilistic forecasts TPT facilitates estimation of dominant transition pathways and return times of transitions. These "dynamical statistics" are obtained by solving partial differential equations in the model's phase space. With future application to more complex models, TPT and its constituent quantities promise to improve the predictability of extreme weather events, through both generation and principled evaluation of forecasts. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: 17 pages, 12 figures. To be published in the Journal of the Atmospheric Sciences

Simulations of Water Vapor and Clouds on Rapidly Rotating and Tidally Locked Planets: a 3D Model Intercomparison

Authors: Jun Yang, Jeremy Leconte, Eric T. Wolf, Timonthy Merlis, Daniel D. B. Koll, Francois Forget, Dorian S. Abbot

Abstract: Robustly modeling the inner edge of the habitable zone is essential for determining the most promising potentially habitable exoplanets for atmospheric characterization. Global climate models (GCMs) have become the standard tool for calculating this boundary, but divergent results have emerged among the various GCMs. In this study, we perform an inter-comparison of standard GCMs used in the field… ▽ More Robustly modeling the inner edge of the habitable zone is essential for determining the most promising potentially habitable exoplanets for atmospheric characterization. Global climate models (GCMs) have become the standard tool for calculating this boundary, but divergent results have emerged among the various GCMs. In this study, we perform an inter-comparison of standard GCMs used in the field on a rapidly rotating planet receiving a G-star spectral energy distribution and on a tidally locked planet receiving an M-star spectral energy distribution. Experiments both with and without clouds are examined. We find relatively small difference (within 8 K) in global-mean surface temperature simulation among the models in the G-star case with clouds. In contrast, the global-mean surface temperature simulation in the M-star case is highly divergent (20-30 K). Moreover, even differences in the simulated surface temperature when clouds are turned off are significant. These differences are caused by differences in cloud simulation and/or radiative transfer, as well as complex interactions between atmospheric dynamics and these two processes. For example we find that an increase in atmospheric absorption of shortwave radiation can lead to higher relative humidity at high altitudes globally and, therefore, a significant decrease in planetary radiation emitted to space. This study emphasizes the importance of basing conclusions about planetary climate on simulations from a variety of GCMs and motivates the eventual comparison of GCM results with terrestrial exoplanet observations to improve their performance. △ Less

Submitted 24 December, 2019; originally announced December 2019.

Comments: 34 pages, 17 figures

Journal ref: The Astrophysical Journal 2019

Clouds will likely prevent the detection of water vapor in JWST transmission spectra of terrestrial exoplanets

Authors: Thaddeus D. Komacek, Thomas J. Fauchez, Eric T. Wolf, Dorian S. Abbot

Abstract: We are on the verge of characterizing the atmospheres of terrestrial exoplanets in the habitable zones of M dwarf stars. Due to their large planet-to-star radius ratios and higher frequency of transits, terrestrial exoplanets orbiting M dwarf stars are favorable for transmission spectroscopy. In this work, we quantify the effect that water clouds have on the amplitude of water vapor transmission s… ▽ More We are on the verge of characterizing the atmospheres of terrestrial exoplanets in the habitable zones of M dwarf stars. Due to their large planet-to-star radius ratios and higher frequency of transits, terrestrial exoplanets orbiting M dwarf stars are favorable for transmission spectroscopy. In this work, we quantify the effect that water clouds have on the amplitude of water vapor transmission spectral features of terrestrial exoplanets orbiting M dwarf stars. To do so, we make synthetic transmission spectra from general circulation model (GCM) experiments of tidally locked planets. We improve upon previous work by considering how varying a broad range of planetary parameters affects transmission spectra. We find that clouds lead to a 10-100 times increase in the number of transits required to detect water features with the James Webb Space Telescope (JWST) with varying rotation period, incident stellar flux, surface pressure, planetary radius, and surface gravity. We also find that there is a strong increase in the dayside cloud coverage in our GCM simulations with rotation periods $\gtrsim 12 \ \mathrm{days}$ for planets with Earth's radius. This increase in cloud coverage leads to even stronger muting of spectral features for slowly rotating exoplanets orbiting M dwarf stars. We predict that it will be extremely challenging to detect water transmission features in the atmospheres of terrestrial exoplanets in the habitable zone of M dwarf stars with JWST. However, species that are well-mixed above the cloud deck (e.g., CO$_2$ and CH$_4$) may still be detectable on these planets with JWST. △ Less

Submitted 18 December, 2019; originally announced December 2019.

Comments: 7 pages, 3 figures, accepted at ApJ Letters

Journal ref: The Astrophysical Journal Letters, 888:L20, 2020 January 10

No Snowball on Habitable Tidally Locked Planets with a Dynamic Ocean

Authors: Jade H. Checlair, Stephanie L. Olson, Malte F. Jansen, Dorian S. Abbot

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

Submitted 14 October, 2019; originally announced October 2019.

Comments: Accepted at ApJL Sept 25 2019

The Snowball Stratosphere

Authors: R. J. Graham, Tiffany Shaw, Dorian Abbot

Abstract: According to the Snowball Earth hypothesis, Earth has experienced periods of low-latitude glaciation in its deep past. Prior studies have used general circulation models (GCMs) to examine the effects such an extreme climate state might have on the structure and dynamics of Earth's troposphere, but the behavior of the stratosphere has not been studied in detail. Understanding the snowball stratosph… ▽ More According to the Snowball Earth hypothesis, Earth has experienced periods of low-latitude glaciation in its deep past. Prior studies have used general circulation models (GCMs) to examine the effects such an extreme climate state might have on the structure and dynamics of Earth's troposphere, but the behavior of the stratosphere has not been studied in detail. Understanding the snowball stratosphere is important for developing an accurate account of the Earth's radiative and chemical properties during these episodes. Here we conduct the first analysis of the stratospheric circulation of the Snowball Earth using ECHAM6 general circulation model simulations. In order to understand the factors contributing to the stratospheric circulation, we extend the Statistical Transformed Eulerian Mean framework. We find that the stratosphere during a snowball with prescribed modern ozone levels exhibits a weaker meridional overturning circulation, reduced wave activity, stronger zonal jets, and is extremely cold relative to modern conditions. Notably, the snowball stratosphere displays no sudden stratospheric warmings. Without ozone, the stratosphere displays slightly weaker circulation, a complete lack of polar vortex, and even colder temperatures. We also explicitly quantify for the first time the cross-tropopause mass exchange rate and stratospheric mixing efficiency during the snowball and show that our values do not change the constraints on CO$_2$ inferred from geochemical proxies during the Marinoan glaciation ($\sim$635 Ma), unless the O$_2$ concentration during the snowball was orders of magnitude less than the CO$_2$ concentration. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Comments: 28 pages, 11 figures, submitted to JGR-Atmospheres

Oceanographic Considerations for Exoplanet Life Detection

Authors: Stephanie L. Olson, Malte Jansen, Dorian S. Abbot

Abstract: Liquid water oceans are at the center of our search for life on exoplanets because water is a strict requirement for life as we know it. However, oceans are dynamic habitats---and some oceans may be better hosts for life than others. In Earth's ocean, circulation transports essential nutrients such as phosphate and is a first-order control on the distribution and productivity of life. Of particula… ▽ More Liquid water oceans are at the center of our search for life on exoplanets because water is a strict requirement for life as we know it. However, oceans are dynamic habitats---and some oceans may be better hosts for life than others. In Earth's ocean, circulation transports essential nutrients such as phosphate and is a first-order control on the distribution and productivity of life. Of particular importance is upward flow from the dark depths of the ocean in response to wind-driven divergence in surface layers. This `upwelling' returns essential nutrients that tend to accumulate at depth via sinking of organic particulates back to the sunlit regions where photosynthetic life thrives. Ocean dynamics are likely to impose constraints on the activity and atmospheric expression of photosynthetic life in exo-oceans as well, but we lack an understanding of how ocean dynamics may differ on other planets. We address this issue by exploring the sensitivity of ocean dynamics to a suite of planetary parameters using ROCKE-3D, a fully coupled ocean-atmosphere GCM. Our results suggest that planets that rotate slower and have higher surface pressure than Earth may be the most attractive targets for remote life detection because upwelling is enhanced under these conditions, resulting in greater nutrient supply to the surface biosphere. Seasonal deepening of the mixed layer on high obliquity planets may also enhance nutrient replenishment from depth into the surface mixed layer. Efficient nutrient recycling favors greater biological activity, more biosignature production, and thus more detectable life. More generally, our results demonstrate the importance of considering oceanographic phenomena for exoplanet life detection and motivate future interdisciplinary contributions to the emerging field of exo-oceanography. △ Less

Submitted 20 May, 2020; v1 submitted 6 September, 2019; originally announced September 2019.

Comments: 21 pages, 13 figures. Published May 20, 2020 in ApJ

Scaling Relations for Terrestrial Exoplanet Atmospheres from Baroclinic Criticality

Authors: Thaddeus D. Komacek, Malte F. Jansen, Eric T. Wolf, Dorian S. Abbot

Abstract: The macroturbulent atmospheric circulation of Earth-like planets mediates their equator-to-pole heat transport. For fast-rotating terrestrial planets, baroclinic instabilities in the mid-latitudes lead to turbulent eddies that act to transport heat poleward. In this work, we derive a scaling theory for the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanet atmosphere… ▽ More The macroturbulent atmospheric circulation of Earth-like planets mediates their equator-to-pole heat transport. For fast-rotating terrestrial planets, baroclinic instabilities in the mid-latitudes lead to turbulent eddies that act to transport heat poleward. In this work, we derive a scaling theory for the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanet atmospheres. This theory is built on the work of Jansen & Ferrari (2013), and determines how unstable the atmosphere is to baroclinic instability (the baroclinic "criticality") through a balance between the baroclinic eddy heat flux and radiative heating/cooling. We compare our scaling theory to General Circulation Model (GCM) simulations and find that the theoretical predictions for equator-to-pole temperature contrast and bulk lapse rate broadly agree with GCM experiments with varying rotation rate and surface pressure throughout the baroclincally unstable regime. Our theoretical results show that baroclinic instabilities are a strong control of heat transport in the atmospheres of Earth-like exoplanets, and our scalings can be used to estimate the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanets. These scalings can be tested by spectroscopic retrievals and full-phase light curves of terrestrial exoplanets with future space telescopes. △ Less

Submitted 7 August, 2019; originally announced August 2019.

Comments: Accepted at ApJ, 8 pages, 4 figures

Journal ref: The Astrophysical Journal, Volume 883, Issue 1, article id. 46, 8 pp. (2019)

Identifying Candidate Atmospheres on Rocky M dwarf Planets via Eclipse Photometry

Authors: Daniel D. B. Koll, Matej Malik, Megan Mansfield, Eliza M. -R. Kempton, Edwin Kite, Dorian Abbot, Jacob L. Bean

Abstract: Most rocky planets in the galaxy orbit a cool host star, and there is large uncertainty among theoretical models whether these planets can retain an atmosphere. The James Webb Space Telescope (JWST) might be able to settle this question empirically, but most proposals for doing so require large observational effort because they are based on spectroscopy. Here we show that infrared photometry of se… ▽ More Most rocky planets in the galaxy orbit a cool host star, and there is large uncertainty among theoretical models whether these planets can retain an atmosphere. The James Webb Space Telescope (JWST) might be able to settle this question empirically, but most proposals for doing so require large observational effort because they are based on spectroscopy. Here we show that infrared photometry of secondary eclipses could quickly identify "candidate" atmospheres, by searching for rocky planets with atmospheres thick enough that atmospheric heat transport noticeably reduces their dayside thermal emission compared to that of a bare rock. For a planet amenable to atmospheric follow-up, we find that JWST should be able to confidently detect the heat redistribution signal of an O(1) bar atmosphere with one to two eclipses. One to two eclipses is generally much less than the effort needed to infer an atmosphere via transmission or emission spectroscopy. Candidate atmospheres can be further validated via follow-up spectroscopy or phase curves. In addition, because this technique is fast it could enable a first atmospheric survey of rocky exoplanets with JWST. We estimate that the TESS mission will find ~100 planets that are too hot to be habitable but that can be quickly probed via eclipse photometry. Knowing whether hot, rocky planets around M dwarfs have atmospheres is important not only for understanding the evolution of uninhabitable worlds: if atmospheres are common on hot planets, then cooler, potentially habitable planets around M dwarfs are also likely to have atmospheres. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: Submitted to ApJ. Also see these three companion papers: 1. Mansfield et al (submitted), "Identifying Atmospheres on Rocky Exoplanets Through Inferred High Albedo", 2. Malik et al (submitted), "Analyzing Atmospheric Temperature Profiles and Spectra of M dwarf Rocky Planets", and 3. Koll (submitted) "A Scaling Theory for Atmospheric Heat Redistribution on Rocky Exoplanets"

Robustness of Gaian Feedbacks to Climate Perturbations

Authors: Olivia D. N. Alcabes, Stephanie Olson, Dorian S. Abbot

Abstract: The Gaia hypothesis postulates that life regulates its environment to be favorable for its own survival. Most planets experience numerous perturbations throughout their lifetimes such as asteroid impacts, volcanism, and the evolution of their host star's luminosity. For the Gaia hypothesis to be viable, life must be able to keep the conditions of its host planet habitable, even in the face of thes… ▽ More The Gaia hypothesis postulates that life regulates its environment to be favorable for its own survival. Most planets experience numerous perturbations throughout their lifetimes such as asteroid impacts, volcanism, and the evolution of their host star's luminosity. For the Gaia hypothesis to be viable, life must be able to keep the conditions of its host planet habitable, even in the face of these challenges. ExoGaia, a model created to investigate the Gaia hypothesis, has been previously used to demonstrate that a randomly mutating biosphere is in some cases capable of maintaining planetary habitability. However, those model scenarios assumed that all non-biological planetary parameters were static, neglecting the inevitable perturbations that real planets would experience. To see how life responds to climate perturbations to its host planet, we created three climate perturbations in ExoGaia: one rapid cooling of a planet and two heating events, one rapid and one gradual. The planets on which Gaian feedbacks emerge without climate perturbations are the same planets on which life is most likely to survive each of our perturbation scenarios. Biospheres experiencing gradual changes to the environment are able to survive changes of larger magnitude than those experiencing rapid perturbations, and the magnitude of change matters more than the sign. These findings suggest that if the Gaia hypothesis is correct, then typical perturbations that a planet would experience may be unlikely to disrupt Gaian systems. △ Less

Submitted 14 January, 2020; v1 submitted 3 June, 2019; originally announced June 2019.

Comments: 6 pages, 6 figures; accepted for publication in MNRAS

Maximizing simulated tropical cyclone intensity with action minimization

Authors: David A. Plotkin, Robert J. Webber, Morgan E O'Neill, Jonathan Weare, Dorian S. Abbot

Abstract: Direct computer simulation of intense tropical cyclones (TCs) in weather models is limited by computational expense. Intense TCs are rare and have small-scale structures, making it difficult to produce large ensembles of storms at high resolution. Further, models often fail to capture the process of rapid intensification, which is a distinguishing feature of many intense TCs. Understanding rapid i… ▽ More Direct computer simulation of intense tropical cyclones (TCs) in weather models is limited by computational expense. Intense TCs are rare and have small-scale structures, making it difficult to produce large ensembles of storms at high resolution. Further, models often fail to capture the process of rapid intensification, which is a distinguishing feature of many intense TCs. Understanding rapid intensification is especially important in the context of global warming, which may increase the frequency of intense TCs. To better leverage computational resources for the study of rapid intensification, we introduce an action minimization algorithm applied to the WRF and WRFPLUS models. Action minimization nudges the model into forming more intense TCs than it otherwise would; it does so via the maximum likelihood path in a stochastic formulation of the model, thereby allowing targeted study of intensification mechanisms. We apply action minimization to simulations of Hurricanes Danny (2015) and Fred (2009) at 6 km resolution to demonstrate that the algorithm consistently intensifies TCs via physically plausible pathways. We show an approximately ten-fold computational savings using action minimization to study the tail of the TC intensification distribution. Further, for Hurricanes Danny and Fred, action minimization produces perturbations that preferentially reduce low-level shear as compared to upper-level shear, at least above a threshold of approximately $4 \mathrm{\ m \ s^{-1}}$. We also demonstrate that asymmetric, time-dependent patterns of heating can cause significant TC intensification beyond symmetric, azimuthally-averaged heating and find a regime of non-linear response to asymmetric heating that has not been extensively studied in previous work. △ Less

Submitted 1 May, 2019; originally announced May 2019.

Effects of Radius and Gravity on the Inner Edge of the Habitable Zone

Authors: Huanzhou Yang, Thaddeus D. Komacek, Dorian S. Abbot

Abstract: A rigorous definition of the habitable zone and its dependence on planetary properties is part of the search for habitable exoplanets. In this work, we use the general circulation model ExoCAM to determine how the inner edge of the habitable zone of tidally locked planets orbiting M dwarf stars depends on planetary radius, surface gravity, and surface pressure. We find that the inner edge of the h… ▽ More A rigorous definition of the habitable zone and its dependence on planetary properties is part of the search for habitable exoplanets. In this work, we use the general circulation model ExoCAM to determine how the inner edge of the habitable zone of tidally locked planets orbiting M dwarf stars depends on planetary radius, surface gravity, and surface pressure. We find that the inner edge of the habitable zone for more massive planets occurs at higher stellar irradiation, as found in previous one-dimensional simulations. We also determine the relative effects of varying planetary radius and surface gravity. Increasing the planetary radius leads to a lower planetary albedo and warmer climate, pushing the inner edge of the habitable zone to lower stellar irradiation. This results from a change in circulation regime that leads to the disruption of the thick, reflective cloud deck around the substellar point. Increasing gravity increases the outgoing longwave radiation, which moves the inner edge of the habitable zone to higher stellar irradiation. This is because the column mass of water vapor decreases with increasing gravity, leading to a reduction in the greenhouse effect. The effect of gravity on the outgoing longwave radiation is stronger than the effect of radius on the planetary albedo, so that increasing gravity and radius together causes the inner edge of the habitable zone to move to higher stellar irradiation. Our results show that the inner edge of the habitable zone for more massive terrestrial planets occurs at a larger stellar irradiation. △ Less

Submitted 28 April, 2019; originally announced April 2019.

Comments: 7 pages, 4 figures, 1 table, Accepted at ApJL

Practical rare event sampling for extreme mesoscale weather

Authors: Robert J. Webber, David A. Plotkin, Morgan E O'Neill, Dorian S. Abbot, Jonathan Weare

Abstract: Extreme mesoscale weather, including tropical cyclones, squall lines, and floods, can be enormously damaging and yet challenging to simulate; hence, there is a pressing need for more efficient simulation strategies. Here we present a new rare event sampling algorithm called Quantile Diffusion Monte Carlo (Quantile DMC). Quantile DMC is a simple-to-use algorithm that can sample extreme tail behavio… ▽ More Extreme mesoscale weather, including tropical cyclones, squall lines, and floods, can be enormously damaging and yet challenging to simulate; hence, there is a pressing need for more efficient simulation strategies. Here we present a new rare event sampling algorithm called Quantile Diffusion Monte Carlo (Quantile DMC). Quantile DMC is a simple-to-use algorithm that can sample extreme tail behavior for a wide class of processes. We demonstrate the advantages of Quantile DMC compared to other sampling methods and discuss practical aspects of implementing Quantile DMC. To test the feasibility of Quantile DMC for extreme mesoscale weather, we sample extremely intense realizations of two historical tropical cyclones, 2010 Hurricane Earl and 2015 Hurricane Joaquin. Our results demonstrate Quantile DMC's potential to provide low-variance extreme weather statistics while highlighting the work that is necessary for Quantile DMC to attain greater efficiency in future applications. △ Less

Submitted 6 April, 2019; originally announced April 2019.

Comments: 18 pages, 9 figures

A Statistical Comparative Planetology Approach to Maximize the Scientific Return of Future Exoplanet Characterization Efforts

Authors: Jade H. Checlair, Dorian S. Abbot, Robert J. Webber, Y. Katherina Feng, Jacob L. Bean, Edward W. Schwieterman, Christopher C. Stark, Tyler D. Robinson, Eliza Kempton, Olivia D. N. Alcabes, Daniel Apai, Giada Arney, Nicolas Cowan, Shawn Domagal-Goldman, Chuanfei Dong, David P. Fleming, Yuka Fujii, R. J. Graham, Scott D. Guzewich, Yasuhiro Hasegawa, Benjamin P. C. Hayworth, Stephen R. Kane, Edwin S. Kite, Thaddeus D. Komacek, Ravi K. Kopparapu , et al. (15 additional authors not shown)

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

Submitted 12 March, 2019; originally announced March 2019.

Comments: White paper submitted in response to the solicitation of feedback for the "2020 Decadal Survey" by the National Academy of Sciences

The Importance of 3D General Circulation Models for Characterizing the Climate and Habitability of Terrestrial Extrasolar Planets

Authors: Eric T. Wolf, Ravi Kopparapu, Vladimir Airapetian, Thomas Fauchez, Scott D. Guzewich, Stephen R. Kane, Daria Pidhorodetska, Michael J. Way, Dorian S. Abbot, Jade H. Checlair, Christopher E. Davis, Anthony Del Genio, Chaunfei Dong, Siegfried Eggl, David P. Fleming, Yuka Fujii, Nader Haghighipour, Nicholas Heavens, Wade G. Henning, Nancy Y. Kiang, Mercedes Lopez-Morales, Jacob Lustig-Yaeger, Vikki Meadows, Christopher T. Reinhard, Sarah Rugheimer , et al. (5 additional authors not shown)

Abstract: While recently discovered exotic new planet-types have both challenged our imaginations and broadened our knowledge of planetary system workings, perhaps the most compelling objective of exoplanet science is to detect and characterize habitable and possibly inhabited worlds orbiting in other star systems. For the foreseeable future, characterizations of extrasolar planets will be made via remote s… ▽ More While recently discovered exotic new planet-types have both challenged our imaginations and broadened our knowledge of planetary system workings, perhaps the most compelling objective of exoplanet science is to detect and characterize habitable and possibly inhabited worlds orbiting in other star systems. For the foreseeable future, characterizations of extrasolar planets will be made via remote sensing of planetary spectroscopic and temporal signals, along with careful fitting of this data to advanced models of planets and their atmospheres. Terrestrial planets are small and significantly more challenging to observe compared to their larger gaseous brethren; however observatories coming on-line in the coming decade will begin to allow their characterization. Still, it is not enough to invest only in observational endeavors. Comprehensive modeling of planetary atmospheres is required in order to fully understand what it is that our grand telescopes see in the night-sky. In our quest to characterize habitable, and possibly inhabited worlds, 3D general circulation models (GCMs) should be used to evaluate potential climate states and their associated temporal and spatial dependent observable signals. 3D models allow for coupled, self-consistent, multi-dimensional simulations, which can realistically simulate the climates of terrestrial extrasolar planets. A complete theoretical understanding of terrestrial exoplanetary atmospheres, gained through comprehensive 3D modeling, is critical for interpreting spectra of exoplanets taken from current and planned instruments, and is critical for designing future missions that aim to measure spectra of potentially habitable exoplanets as one of their key science goals. We recommend continued institutional support for 3D GCM modeling teams that focus on planetary and exoplanetary applications. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: NAS ASTRO 2020, White Paper! 2 figures, and 10 pages

Habitable zone predictions and how to test them

Authors: Ramses M. Ramirez, Dorian S. Abbot, Yuka Fujii, Keiko Hamano, Edwin Kite, Amit Levi, Manasvi Lingam, Theresa Lueftinger, Tyler D. Robinson, Andrew Rushby, Laura Schaefer, Elizabeth Tasker, Giovanni Vladilo, Robin D. Wordsworth

Abstract: The habitable zone (HZ) is the region around a star(s) where standing bodies of water could exist on the surface of a rocky planet. The classical HZ definition makes a number of assumptions common to the Earth, including assuming that the most important greenhouse gases for habitable planets are CO2 and H2O, habitable planets orbit main-sequence stars, and that the carbonate-silicate cycle is a un… ▽ More The habitable zone (HZ) is the region around a star(s) where standing bodies of water could exist on the surface of a rocky planet. The classical HZ definition makes a number of assumptions common to the Earth, including assuming that the most important greenhouse gases for habitable planets are CO2 and H2O, habitable planets orbit main-sequence stars, and that the carbonate-silicate cycle is a universal process on potentially habitable planets. Here, we discuss these and other predictions for the habitable zone and the observations that are needed to test them. We also, for the first time, argue why A-stars may be interesting HZ prospects. Instead of relying on unverified extrapolations from our Earth, we argue that future habitability studies require first principles approaches where temporal, spatial, physical, chemical, and biological systems are dynamically coupled. We also suggest that next-generation missions are only the beginning of a much more data-filled era in the not-too-distant future, when possibly hundreds to thousands of HZ planets will yield the statistical data we need to go beyond just finding habitable zone planets to actually determining which ones are most likely to exhibit life. △ Less

Submitted 11 March, 2019; v1 submitted 8 March, 2019; originally announced March 2019.

Comments: White paper submitted to the NAS Astro 2020 decadal survey(5 pages, 1 figure + cover page, 13 co-authors + 10 co-signers). First author name added to description. In this arxiv version, one more co-signer was added after the NAS version was submitted. Link: http://surveygizmoresponseuploads.s3.amazonaws.com/fileuploads/623127/4458621/64-2dc3fcc201a3075e3da6233f1a91c668_RamirezRamsesM.pdf

Ocean Dynamics and the Inner Edge of the Habitable Zone for Tidally Locked Terrestrial Planets

Authors: Jun Yang, Dorian S. Abbot, Daniel D. B. Koll, Yongyun Hu, Adam P. Showman

Abstract: Recent studies have shown that ocean dynamics can have a significant warming effect on the permanent night sides of 1 to 1 tidally locked terrestrial exoplanets with Earth-like atmospheres and oceans in the middle of the habitable zone. However, the impact of ocean dynamics on the habitable zone's boundaries (inner edge and outer edge) is still unknown and represents a major gap in our understandi… ▽ More Recent studies have shown that ocean dynamics can have a significant warming effect on the permanent night sides of 1 to 1 tidally locked terrestrial exoplanets with Earth-like atmospheres and oceans in the middle of the habitable zone. However, the impact of ocean dynamics on the habitable zone's boundaries (inner edge and outer edge) is still unknown and represents a major gap in our understanding of this type of planets. Here we use a coupled atmosphere-ocean global climate model to show that planetary heat transport from the day to night side is dominated by the ocean at lower stellar fluxes and by the atmosphere near the inner edge of the habitable zone. This decrease in oceanic heat transport (OHT) at high stellar fluxes is mainly due to weakening of surface wind stress and a decrease in surface shortwave energy deposition. We further show that ocean dynamics have almost no effect on the observational thermal phase curves of planets near the inner edge of the habitable zone. For planets in the habitable zone's middle range, ocean dynamics moves the hottest spot on the surface eastward from the substellar point. These results suggest that future studies of the inner edge may devote computational resources to atmosphere-only processes such as clouds and radiation. For studies of the middle range and outer edge of the habitable zone, however, fully coupled ocean-atmosphere modeling will be necessary. Note that due to computational resource limitations, only one rotation period (60 Earth days) has been systematically examined in this study; future work varying rotation period as well as other parameters such as atmospheric mass and composition is required. △ Less

Submitted 6 February, 2019; originally announced February 2019.

Comments: 34 pages, 13 figures, and 1 table

Journal ref: The Astrophysical Journal, 2019

The atmospheric circulation and climate of terrestrial planets orbiting Sun-like and M-dwarf stars over a broad range of planetary parameters

Authors: Thaddeus D. Komacek, Dorian S. Abbot

Abstract: The recent detections of temperate terrestrial planets orbiting nearby stars and the promise of characterizing their atmospheres motivates a need to understand how the diversity of possible planetary parameters affects the climate of terrestrial planets. In this work, we investigate the atmospheric circulation and climate of terrestrial exoplanets orbiting both Sun-like and M-dwarf stars over a wi… ▽ More The recent detections of temperate terrestrial planets orbiting nearby stars and the promise of characterizing their atmospheres motivates a need to understand how the diversity of possible planetary parameters affects the climate of terrestrial planets. In this work, we investigate the atmospheric circulation and climate of terrestrial exoplanets orbiting both Sun-like and M-dwarf stars over a wide swath of possible planetary parameters, including the planetary rotation period, surface pressure, incident stellar flux, surface gravity, planetary radius, and cloud particle size. We do so using a general circulation model (GCM) that includes non-grey radiative transfer and the effects of clouds. The results from this suite of simulations generally show qualitatively similar dependencies of circulation and climate on planetary parameters as idealized GCMs, with quantitative differences due to the inclusion of additional model physics. Notably, we find that the effective cloud particle size is a key unknown parameter that can greatly affect the climate of terrestrial exoplanets. We confirm a transition between low and high dayside cloud coverage of synchronously rotating terrestrial planets with increasing rotation period. We determine that this cloud transition is due to eddy-driven convergence near the substellar point and should not be parameterization-dependent. Finally, we compute full-phase light curves from our simulations of planets orbiting M-dwarf stars, finding that changing incident stellar flux and rotation period affect observable properties of terrestrial exoplanets. Our GCM results can guide expectations for planetary climate over the broad range of possible terrestrial exoplanets that will be observed with future space telescopes. △ Less

Submitted 8 February, 2019; v1 submitted 2 January, 2019; originally announced January 2019.

Comments: 21 pages, 19 figures, 5 tables. Updated to reflect published version

Journal ref: The Astrophysical Journal, 871:245 (20pp), 2019 February 1

Differences in water vapor radiative transfer among 1D models can significantly affect the inner edge of the habitable zone

Authors: Jun Yang, Jeremy Leconte, Eric T. Wolf, Colin Goldblatt, Nicole Feldl, Timothy Merlis, Yuwei Wang, Daniel D. B. Koll, Feng Ding, Francois Forget, Dorian S. Abbot

Abstract: An accurate estimate of the inner edge of the habitable zone is critical for determining which exoplanets are potentially habitable and for designing future telescopes to observe them. Here, we explore differences in estimating the inner edge among seven one-dimensional (1D) radiative transfer models: two line-by-line codes (SMART and LBLRTM) as well as five band codes (CAM3, CAM4_Wolf, LMDG, SBDA… ▽ More An accurate estimate of the inner edge of the habitable zone is critical for determining which exoplanets are potentially habitable and for designing future telescopes to observe them. Here, we explore differences in estimating the inner edge among seven one-dimensional (1D) radiative transfer models: two line-by-line codes (SMART and LBLRTM) as well as five band codes (CAM3, CAM4_Wolf, LMDG, SBDART, and AM2) that are currently being used in global climate models. We compare radiative fluxes and spectra in clear-sky conditions around G- and M-stars, with fixed moist adiabatic profiles for surface temperatures from 250 to 360 K. We find that divergences among the models arise mainly from large uncertainties in water vapor absorption in the window region (10 um) and in the region between 0.2 and 1.5 um. Differences in outgoing longwave radiation increase with surface temperature and reach 10-20 Wm^-2; differences in shortwave reach up to 60 Wm^-2, especially at the surface and in the troposphere, and are larger for an M-dwarf spectrum than a solar spectrum. Differences between the two line-by-line models are significant, although smaller than among the band models. Our results imply that the uncertainty in estimating the insolation threshold of the inner edge (the runaway greenhouse limit) due only to clear-sky radiative transfer is ~10% of modern Earth's solar constant (i.e., ~34 Wm^-2 in global mean) among band models and ~3% between the two line-by-line models. These comparisons show that future work is needed focusing on improving water vapor absorption coefficients in both shortwave and longwave, as well as on increasing the resolution of stellar spectra in broadband models. △ Less

Submitted 5 September, 2018; originally announced September 2018.

Comments: 25 pages, 1 table, and 10 figures

Journal ref: The Astrophysical Journal, 826:222, 2016;

The Atmosphere

Authors: Dorian S. Abbot

Abstract: These notes contain everything necessary to run a flipped course on "The Atmosphere" at an introductory undergraduate level. There are notes for the students to read before each course meeting and problems for them to work on in small groups during course meetings. Topics include (1) atmospheric composition, structure, and thermodynamics; (2) solar and terrestrial radiation in the atmospheric ener… ▽ More These notes contain everything necessary to run a flipped course on "The Atmosphere" at an introductory undergraduate level. There are notes for the students to read before each course meeting and problems for them to work on in small groups during course meetings. Topics include (1) atmospheric composition, structure, and thermodynamics; (2) solar and terrestrial radiation in the atmospheric energy balance; (3) atmospheric dynamics and circulation. I include 10 problem sets, six practice midterms, and three practice finals. Problems are drawn from the atmospheres of modern and past Earth, solar system planets, and extrasolar planets. I can provide solutions to the in-class problems and problem sets to teachers upon request. You are free to use these notes in your classes, and to expand them as you please. If you catch any typos or errors, please send them to me. Enjoy! △ Less

Submitted 25 June, 2019; v1 submitted 8 June, 2018; originally announced June 2018.

Comments: 215 pages

Simple rules govern the patterns of Arctic sea ice melt ponds

Authors: Predrag Popović, B. B. Cael, Mary Silber, Dorian S. Abbot

Abstract: Climate change, amplified in the far north, has led to rapid sea ice decline in recent years. In the summer, melt ponds form on the surface of Arctic sea ice, significantly lowering the ice reflectivity (albedo) and thereby accelerating ice melt. Pond geometry controls the details of this crucial feedback; however, a reliable model of pond geometry does not currently exist. Here we show that a sim… ▽ More Climate change, amplified in the far north, has led to rapid sea ice decline in recent years. In the summer, melt ponds form on the surface of Arctic sea ice, significantly lowering the ice reflectivity (albedo) and thereby accelerating ice melt. Pond geometry controls the details of this crucial feedback; however, a reliable model of pond geometry does not currently exist. Here we show that a simple model of voids surrounding randomly sized and placed overlapping circles reproduces the essential features of pond patterns. The only two model parameters, characteristic circle radius and coverage fraction, are chosen by comparing, between the model and the aerial photographs of the ponds, two correlation functions which determine the typical pond size and their connectedness. Using these parameters, the void model robustly reproduces the ponds' area-perimeter and area-abundance relationships over more than 6 orders of magnitude. By analyzing the correlation functions of ponds on several dates, we also find that the pond scale and the connectedness are surprisingly constant across different years and ice types. Moreover, we find that ponds resemble percolation clusters near the percolation threshold. These results demonstrate that the geometry and abundance of Arctic melt ponds can be simply described, which can be exploited in future models of Arctic melt ponds that would improve predictions of the response of sea ice to Arctic warming. △ Less

Submitted 12 April, 2018; originally announced April 2018.

Comments: Main article: 5 pages, 3 figures; Supplementary Material: 19 pages, 8 figures

Journal ref: Phys. Rev. Lett., 120, no. 14 (2018): 148701

The continued importance of habitability studies

Authors: Ramses M. Ramirez, Dorian S. Abbot, Vladimir Airapetian, Yuka Fujii, Keiko Hamano, Amit Levi, Tyler D. Robinson, Laura Schaefer, Eric T. Wolf, Robin D. Wordsworth

Abstract: This is a white paper in response to the National Academy of Sciences "Exoplanet Science Strategy" call. We summarize recent advances in theoretical habitability studies and argue that such studies will remain important for guiding and interpreting observations. Interactions between 1-D and 3-D climate modelers will be necessary to resolve recent discrepancies in model results and improve habita… ▽ More This is a white paper in response to the National Academy of Sciences "Exoplanet Science Strategy" call. We summarize recent advances in theoretical habitability studies and argue that such studies will remain important for guiding and interpreting observations. Interactions between 1-D and 3-D climate modelers will be necessary to resolve recent discrepancies in model results and improve habitability studies. Observational capabilities will also need improvement. Although basic observations can be performed with present capabilities, technological advances will be necessary to improve climate models to the level needed for planetary habitability studies. △ Less

Submitted 3 March, 2019; v1 submitted 1 March, 2018; originally announced March 2018.

Comments: This is a white paper submitted to the National Academies 2018 "Exoplanet Science Strategy" call(6 pages total, including cover page). Corrected references section in this version. http://sites.nationalacademies.org/SSB/CurrentProjects/SSB_180659

Decrease in hysteresis of planetary climate for planets with long solar days

Authors: Dorian S. Abbot, Jonah Bloch-Johnson, Jade Checlair, Navah X. Farahat, R. J. Graham, David Plotkin, Predrag Popovic, and Francisco Spaulding-Astudillo

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

Submitted 31 January, 2018; originally announced January 2018.

Comments: accepted at ApJ

No snowball on habitable tidally locked planets

Authors: Jade Checlair, Kristen Menou, Dorian S. Abbot

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

Submitted 14 July, 2017; v1 submitted 24 May, 2017; originally announced May 2017.

Comments: Accepted at the Astrophysical Journal on 07/11/17

A Statistical Comparative Planetology Approach to the Hunt for Habitable Exoplanets and Life Beyond the Solar System

Authors: Jacob L. Bean, Dorian S. Abbot, Eliza M. -R. Kempton

Abstract: The search for habitable exoplanets and life beyond the Solar System is one of the most compelling scientific opportunities of our time. Nevertheless, the high cost of building facilities that can address this topic and the keen public interest in the results of such research requires the rigorous development of experiments that can deliver a definitive advance in our understanding. Most work to d… ▽ More The search for habitable exoplanets and life beyond the Solar System is one of the most compelling scientific opportunities of our time. Nevertheless, the high cost of building facilities that can address this topic and the keen public interest in the results of such research requires the rigorous development of experiments that can deliver a definitive advance in our understanding. Most work to date in this area has focused on a "systems science" approach of obtaining and interpreting comprehensive data for individual planets to make statements about their habitability and the possibility that they harbor life. This strategy is challenging because of the diversity of exoplanets, both observed and expected, and the limited information that can be obtained with astronomical instruments. Here we propose a complementary approach that is based on performing surveys of key planetary characteristics and using statistical marginalization to answer broader questions than can be addressed with a small sample of objects. The fundamental principle of this comparative planetology approach is maximizing what can be learned from each type of measurement by applying it widely rather than requiring that multiple kinds of observations be brought to bear on a single object. As a proof of concept, we outline a survey of terrestrial exoplanet atmospheric water and carbon dioxide abundances that would test the habitable zone hypothesis and lead to a deeper understanding of the frequency of habitable planets. We also discuss ideas for additional surveys that could be developed to test other foundational hypotheses is this area. △ Less

Submitted 17 May, 2017; originally announced May 2017.

Comments: ApJL in press

Effect of surface-mantle water exchange parameterizations on exoplanet ocean depths

Authors: Thaddeus D. Komacek, Dorian S. Abbot

Abstract: Terrestrial exoplanets in the canonical habitable zone may have a variety of initial water fractions due to random volatile delivery by planetesimals. If the total planetary water complement is high, the entire surface may be covered in water, forming a "waterworld." On a planet with active tectonics, competing mechanisms act to regulate the abundance of water on the surface by determining the par… ▽ More Terrestrial exoplanets in the canonical habitable zone may have a variety of initial water fractions due to random volatile delivery by planetesimals. If the total planetary water complement is high, the entire surface may be covered in water, forming a "waterworld." On a planet with active tectonics, competing mechanisms act to regulate the abundance of water on the surface by determining the partitioning of water between interior and surface. Here we explore how the incorporation of different mechanisms for the degassing and regassing of water changes the volatile evolution of a planet. For all of the models considered, volatile cycling reaches an approximate steady-state after $\sim 2 \ \mathrm{Gyr}$. Using these steady-states, we find that if volatile cycling is either solely dependent on temperature or seafloor pressure, exoplanets require a high abundance ($\gtrsim 0.3\%$ of total mass) of water to have fully inundated surfaces. However, if degassing is more dependent on seafloor pressure and regassing mainly dependent on mantle temperature, the degassing rate is relatively large at late times and a steady-state between degassing and regassing is reached with a substantial surface water fraction. If this hybrid model is physical, super-Earths with a total water fraction similar to that of the Earth can become waterworlds. As a result, further understanding of the processes that drive volatile cycling on terrestrial planets is needed to determine the water fraction at which they are likely to become waterworlds. △ Less

Submitted 15 September, 2016; originally announced September 2016.

Comments: Accepted at ApJ

Journal ref: ApJ, 832, 54 (2016)

Analytical investigation of the decrease in the size of the habitable zone due to limited CO$_2$ outgassing rate

Authors: Dorian S. Abbot

Abstract: The habitable zone concept is important because it focuses the scientific search for extraterrestrial life and aids the planning of future telescopes. Recent work has shown that planets near the outer edge of the habitable zone might not actually be able to stay warm and habitable if CO$_2$ outgassing rates are not large enough to maintain high CO$_2$ partial pressures against removal by silicate… ▽ More The habitable zone concept is important because it focuses the scientific search for extraterrestrial life and aids the planning of future telescopes. Recent work has shown that planets near the outer edge of the habitable zone might not actually be able to stay warm and habitable if CO$_2$ outgassing rates are not large enough to maintain high CO$_2$ partial pressures against removal by silicate weathering. In this paper I use simple equations for the climate and CO$_2$ budget of a planet in the habitable zone that can capture the qualitative behavior of the system. With these equations I derive an analytical formula for an effective outer edge of the habitable zone, including limitations imposed by the CO$_2$ outgassing rate. I then show that climate cycles between a Snowball state and a warm climate are only possible beyond this limit if the weathering rate in the Snowball climate is smaller than the CO$_2$ outgassing rate (otherwise stable Snowball states result). I derive an analytical solution for the climate cycles including a formula for their period in this limit. This work allows us to explore the qualitative effects of weathering processes on the effective outer edge of the habitable zone, which is important because weathering parameterizations are uncertain. △ Less

Submitted 9 June, 2016; originally announced June 2016.

Comments: accepted at ApJ

Temperature Structure and Atmospheric Circulation of Dry, Tidally Locked Rocky Exoplanets

Authors: Daniel D. B. Koll, Dorian S. Abbot

Abstract: Next-generation space telescopes will observe the atmospheres of rocky planets orbiting nearby M-dwarfs. Understanding these observations will require well-developed theory in addition to numerical simulations. Here we present theoretical models for the temperature structure and atmospheric circulation of dry, tidally locked rocky exoplanets with grey radiative transfer and test them using a gener… ▽ More Next-generation space telescopes will observe the atmospheres of rocky planets orbiting nearby M-dwarfs. Understanding these observations will require well-developed theory in addition to numerical simulations. Here we present theoretical models for the temperature structure and atmospheric circulation of dry, tidally locked rocky exoplanets with grey radiative transfer and test them using a general circulation model (GCM). First, we develop a radiative-convective model that captures surface temperatures of slowly rotating and cool atmospheres. Second, we show that the atmospheric circulation acts as a global heat engine, which places strong constraints on large-scale wind speeds. Third, we develop a radiative-convective-subsiding model which extends our radiative-convective model to hot and thin atmospheres. We find that rocky planets develop large day-night temperature gradients at a ratio of wave-to-radiative timescales up to two orders of magnitude smaller than the value suggested by work on hot Jupiters. The small ratio is due to the heat engine inefficiency and asymmetry between updrafts and subsidence in convecting atmospheres. Fourth, we show using GCM simulations that rotation only has a strong effect on temperature structure if the atmosphere is hot or thin. Our models let us map out atmospheric scenarios for planets such as GJ 1132b and show how thermal phase curves could constrain them. Measuring phase curves of short-period planets will require similar amounts of time on the James Webb Space Telescope as detecting molecules via transit spectroscopy, so future observations should pursue both techniques. △ Less

Submitted 3 May, 2016; originally announced May 2016.

Comments: Accepted in ApJ

Journal ref: Astrophysical Journal 825 (2016)

A proposal for climate stability on H2-greenhouse planets

Authors: Dorian S. Abbot

Abstract: A terrestrial planet in an orbit far outside of the standard habitable zone could maintain surface liquid water as a result of H2-H2 collision-induced absorption by a thick H2 atmosphere. Without a stabilizing climate feedback, however, habitability would be accidental and likely brief. In this letter I propose stabilizing climate feedbacks for such a planet that require only that biological funct… ▽ More A terrestrial planet in an orbit far outside of the standard habitable zone could maintain surface liquid water as a result of H2-H2 collision-induced absorption by a thick H2 atmosphere. Without a stabilizing climate feedback, however, habitability would be accidental and likely brief. In this letter I propose stabilizing climate feedbacks for such a planet that require only that biological functions have an optimal temperature and operate less efficiently at other temperatures. For example, on a planet with a net source of H2 from its interior, H2-consuming life (such as methanogens) could establish a stable climate. If a positive perturbation is added to the equilibrium temperature, H2 consumption by life will increase (cooling the planet) until the equilibrium climate is reestablished. The potential existence of such feedbacks makes H2-warmed planets more attractive astrobiological targets. △ Less

Submitted 24 November, 2015; originally announced November 2015.

Comments: 5 pages, 3 figures, accepted at ApJL