Researcher profile

Natasha E. Batalha

Natasha E. Batalha contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
9works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2024arXiv

Modeling the day-night temperature variations of ultra-hot Jupiters: confronting non-grey general circulation models and observations

Ultra-hot Jupiters (UHJs) are natural laboratories to study extreme physics in planetary atmospheres and their rich observational data sets are yet to be confronted with models with varying complexities at a population level. In this work, we update the general circulation model of Tan & Komacek (2019) to include a non-grey radiative transfer scheme and apply it to simulate the realistic thermal structures, phase-dependent spectra, and wavelength-dependent phase curves of UHJs. We performed grids of models over a large range of equilibrium temperatures and rotation periods for varying assumptions, showing that the fractional day-night brightness temperature differences remain almost constant or slightly increase with increasing equilibrium temperature from the visible to mid-infrared wavelengths. This differs from previous work primarily due to the increasing planetary rotation rate with increasing equilibrium temperature for fixed host star type. Radiative effects of varying atmospheric compositions become more significant in dayside brightness temperature in longer wavelengths. Data-model comparisons of dayside brightness temperatures and phase curve amplitudes as a function of equilibrium temperature are in broad agreement. Observations show a large scatter compared to models even with a range of different assumptions, indicating significantly varying intrinsic properties in the hot Jupiter population. Our cloud-free models generally struggle to match all observations for individual targets with a single set of parameter choices, indicating the need for extra processes for understanding the heat transport of UHJs.

preprint2024arXiv

The impact of spectral line wing cut-off: Recommended standard method with application to MAESTRO opacity database

When computing cross-sections from a line list, the result depends not only on the line strength, but also the line shape, pressure-broadening parameters, and line wing cut-off (i.e., the maximum distance calculated from each line centre). Pressure-broadening can be described using the Lorentz lineshape, but it is known to not represent the true absorption in the far wings. Both theory and experiment have shown that far from the line centre, non-Lorentzian behaviour controls the shape of the wings and the Lorentz lineshape fails to accurately characterize the absorption, leading to an underestimation or overestimation of the opacity continuum depending on the molecular species involved. The line wing cut-off is an often overlooked parameter when calculating absorption cross sections, but can have a significant effect on the appearance of the spectrum since it dictates the extent of the line wing that contributes to the calculation either side of every line centre. Therefore, when used to analyse exoplanet and brown dwarf spectra, an inaccurate choice for the line wing cut-off can result in errors in the opacity continuum, which propagate into the modeled transit spectra, and ultimately impact/bias the interpretation of observational spectra, and the derived composition and thermal structure. Here, we examine the different methods commonly utilized to calculate the wing cut-off and propose a standard practice procedure (i.e., absolute value of 25~cm$^{-1}$ for $P\leqslant$~200~bar and 100~cm$^{-1}$ for $P >$ ~200~bar) to generate molecular opacities which will be used by the open-access {\tt MAESTRO} (Molecules and Atoms in Exoplanet Science: Tools and Resources for Opacities) database. The pressing need for new measurements and theoretical studies of the far-wings is highlighted.

preprint2022arXiv

Cloudy and Cloud-free Thermal Phase Curves with PICASO: Applications to WASP-43b

We present new functionality within PICASO, a state-of-the-art radiative transfer model for exoplanet and brown dwarf atmospheres, by developing a new pipeline that computes phase-resolved thermal emission (thermal phase curves) from three-dimensional (3D) models. Because PICASO is coupled to Virga, an open-source cloud code, we are able to produce cloudy phase curves with different sedimentation efficiencies ($f_{sed}$) and cloud condensate species. We present the first application of this new algorithm to hot Jupiter WASP-43b. Previous studies of the thermal emission of WASP-43b from Kataria et al. found good agreement between cloud-free models and dayside thermal emission, but an overestimation of the nightside flux, for which clouds have been suggested as a possible explanation. We use the temperature and vertical wind structure from the cloud-free 3D general circulation models of Kataria et al. and post-process it using PICASO, assuming that clouds form and affect the spectra. We compare our models to results from Kataria et al., including Hubble Space Telescope Wide-Field Camera 3 (WFC3) observations of WASP-43b from Stevenson et al. In addition, we compute phase curves for Spitzer at 3.6 and 4.5 $μm$ and compare them to observations from Stevenson et al. We are able to closely recover the cloud-free results, even though PICASO utilizes a coarse spatial grid. We find that cloudy phase curves provide much better agreement with the WFC3 and Spitzer nightside data, while still closely matching the dayside emission. This work provides the community with a convenient, user-friendly tool to interpret phase-resolved observations of exoplanet atmospheres using 3D models.

preprint2022arXiv

PICASO 3.0: A One-Dimensional Climate Model for Giant Planets and Brown Dwarfs

Upcoming James Webb Space Telescope (JWST) observations will allow us to study exoplanet and brown dwarf atmospheres in great detail. The physical interpretation of these upcoming high signal-to-noise observations requires precise atmospheric models of exoplanets and brown dwarfs. While several one-dimensional and three-dimensional atmospheric models have been developed in the past three decades, these models have often relied on simplified assumptions like chemical equilibrium and are also often not open-source, which limits their usage and development by the wider community. We present a python-based one-dimensional atmospheric radiative-convective equilibrium model. This model has heritage from the Fortran-based code (Marley et al.,1996} which has been widely used to model the atmospheres of Solar System objects, brown dwarfs, and exoplanets. In short, the basic capability of the original model is to compute the atmospheric state of the object under radiative-convective equilibrium given its effective or internal temperature, gravity, and host--star properties (if relevant). In the new model, which has been included within the well-utilized code-base PICASO, we have added these original features as well as the new capability of self-consistently treating disequilibrium chemistry. This code is widely applicable to Hydrogen-dominated atmospheres (e.g., brown dwarfs and giant planets).

preprint2021arXiv

A New Sedimentation Model for Greater Cloud Diversity in Giant Exoplanets and Brown Dwarfs

The observed atmospheric spectrum of exoplanets and brown dwarfs depends critically on the presence and distribution of atmospheric condensates. The Ackerman & Marley (2001) methodology for predicting the vertical distribution of condensate particles is widely used to study cloudy atmospheres and has recently been implemented in an open-source python package virga. The model relies upon input parameter $f_{\text{sed}}$, the sedimentation efficiency, which until now has been held constant. The relative simplicity of this model renders it useful for retrieval studies due to its rapidly attainable solutions. However, comparisons with more complex microphysical models such as CARMA have highlighted inconsistencies between the two approaches, namely that the cloud parameters needed for radiative transfer produced by virga are dissimilar to those produced by CARMA. To address these discrepancies, we have extended the original Ackerman and Marley methodology in virga to allow for non-constant $f_{\text{sed}}$ values, in particular those that vary with altitude. We discuss one such parameterization and compare the cloud mass mixing ratio produced by virga with constant and variable $f_{\text{sed}}$ profiles to that produced by CARMA. We find that the variable $f_{\text{sed}}$ formulation better captures the profile produced by CARMA with heterogeneous nucleation, yet performs comparatively to constant $f_{\text{sed}}$ for homogeneous nucleation. In general, virga has the capacity to handle any $f_{\text{sed}}$ with an explicit anti-derivative, permitting a plethora of alternative cloud profiles that are otherwise unattainable by constant $f_{\text{sed}}$ values. The ensuing flexibility has the potential to better agree with increasingly complex models and observed data.

preprint2020arXiv

Warming Early Mars with Climate Cycling: The Effect of CO2-H2 Collision-induced Absorption

Explaining the evidence for surface liquid water on early Mars has been a challenge for climate modelers, as the sun was ~30% less luminous during the late-Noachian. We propose that the additional greenhouse forcing of CO2-H2 collision-induced absorption is capable of bringing the surface temperature above freezing and can put early Mars into a limit-cycling regime. Limit cycles occur when insolation is low and CO2 outgassing rates are unable to balance with the rapid drawdown of CO2 during warm weathering periods. Planets in this regime will alternate between global glaciation and transient warm climate phases. This mechanism is capable of explaining the geomorphological evidence for transient warm periods in the martian record. Previous work has shown that collision-induced absorption of CO2-H2 was capable of deglaciating early Mars, but only with high H2 outgassing rates (greater than ~600 Tmol/yr) and at high surface pressures (between 3 to 4 bars). We used new theoretically derived collision-induced absorption coefficients for CO2-H2 to reevaluate the climate limit cycling hypothesis for early Mars. Using the new and stronger absorption coefficients in our 1-dimensional radiative convective model as well as our energy balance model, we find that limit cycling can occur with an H2 outgassing rate as low as ~300 Tmol/yr at surface pressures below 3 bars. Our results agree more closely with paleoparameters for early martian surface pressure and hydrogen abundance.

preprint2019arXiv

The Precision of Mass Measurements Required for Robust Atmospheric Characterization of Transiting Exoplanets

Two of TESS&#39;s major science goals are to measure masses for 50 planets smaller than 4 Earth radii and to discover high-quality targets for atmospheric characterization efforts. It is important that these two goals are linked by quantifying what precision of mass constraint is required to yield robust atmospheric properties of planets. Here, we address this by conducting retrievals on simulated JWST transmission spectra under various assumptions for the degree of uncertainty in the planet&#39;s mass for a representative population of seven planets ranging from terrestrials to warm Neptunes to hot Jupiters. Only for the cloud-free, low metallicity gas giants are we able to infer exoplanet mass from transmission spectroscopy alone, to ~10% accuracy. For low metallicity cases (<4xSolar) we are able to accurately constrain atmospheric properties without prior knowledge of the planet&#39;s mass. For all other cases (including terrestrial-like planets), atmospheric properties can only be inferred with a mass precision of better than ~50%. At this level, though, the widths of the posterior distributions of the atmospheric properties are dominated by the uncertainties in mass. With a precision of ~20%, the widths of the posterior distributions are dominated by the spectroscopic data quality. Therefore, as a rule-of-thumb, we recommend: a ~50% mass precision for initial atmospheric characterization and a ~20% mass precision for more detailed atmospheric analyses.

preprint2016arXiv

Climate Cycling on Early Mars Caused by the Carbonate-Silicate Cycle

For decades, scientists have tried to explain the evidence for fluvial activity on early Mars, but a consensus has yet to emerge regarding the mechanism for producing it. One hypothesis suggests early Mars was warmed by a thick greenhouse atmosphere. Another suggests that early Mars was generally cold but was warmed occasionally by impacts or by episodes of enhanced volcanism. These latter hypotheses struggle to produce the amounts of rainfall needed to form the martian valleys, but are consistent with inferred low rates of weathering compared to Earth. Here, we provide a geophysical mechanism that could have induced cycles of glaciation and deglaciation on early Mars. Our model produces dramatic climate cycles with extended periods of glaciation punctuated by warm periods lasting up to 10 Myr, much longer than those generated in other episodic warming models. The cycles occur because stellar insolation was low, and because CO2 outgassing is not able to keep pace with CO2 consumption by silicate weathering followed by deposition of carbonates. While CO2 by itself is not able to deglaciate early Mars in our model, we assume that the greenhouse effect is enhanced by substantial amounts of H2 outgassed from Mars&#39; reduced crust and mantle. Our hypothesis can be tested by future Mars exploration that better establishes the time scale for valley formation.

preprint2016arXiv

Limit cycles can reduce the width of the habitable zone

The liquid water habitable zone (HZ) describes the orbital distance at which a terrestrial planet can maintain above-freezing conditions through regulation by the carbonate-silicate cycle. Recent calculations have suggested that planets in the outer regions of the habitable zone cannot maintain stable, warm climates, but rather should oscillate between long, globally glaciated states and shorter periods of climatic warmth. Such conditions, similar to &#39;Snowball Earth&#39; episodes experienced on Earth, would be inimical to the development of complex land life, including intelligent life. Here, we build upon previous studies with an updated an energy balance climate model to calculate this &#39;limit cycle&#39; region of the habitable zone where such cycling would occur. We argue that an abiotic Earth would have a greater CO$_2$ partial pressure than today because plants and other biota help to enhance the storage of CO$_2$ in soil. When we tune our abiotic model accordingly, we find that limit cycles can occur but that previous calculations have overestimated their importance. For G stars like the Sun, limit cycles occur only for planets with CO$_2$ outgassing rates less than that on modern Earth. For K and M star planets, limit cycles should not occur; however, M-star planets may be inhospitable to life for other reasons. Planets orbiting late G-type and early K-type stars retain the greatest potential for maintaining warm, stable conditions. Our results suggest that host star type, planetary volcanic activity, and seafloor weathering are all important factors in determining whether planets will be prone to limit cycling.