3D Atmospheric Structure Research Articles

Informed Systematic Method to Identify Variable Mid- and Late-T Dwarfs

The majority of brown dwarfs show some level of photometric or spectrophotometric variability in different wavelength ranges. This variability allow us to trace the 3D atmospheric structures of variable brown dwarfs and directly-imaged exoplanets with radiative-transfer models and mapping codes. Nevertheless, to date, we do not have an informed method to preselect the brown dwarfs that might show a higher variability amplitude for a thorough variability study. In this work, we designed and tested near-infrared spectral indices to preselect the most likely variable mid- and late-T dwarfs, which overlap in effective temperatures with directly-imaged exoplanets. We used time-resolved near-infrared Hubble Space Telescope Wide Field Camera 3 spectra of a T6.5 dwarf, 2MASS J22282889–431026, to design our novel spectral indices. We tested these spectral indices on 26 T5.5–T7.5 near-infrared SpeX/IRTF spectra, and we provided eight new mid- and late-T variable candidates. We estimated the variability fraction of our sample as %, which agrees with the variability fractions provided by Metchev for mid- to late-T dwarfs. In addition, two of the three previously known variables in our sample of SpeX spectra are flagged as variable candidates by our indices. Similarly, all seven known nonvariables in our sample are flagged as nonvariable objects by our indices. These results suggest that our spectral indices might be used to find variable mid- and late-T brown dwarf variables. These indices may be crucial in the future to select cool directly-imaged exoplanets for variability studies.

ThERESA: Three-dimensional Eclipse Mapping with Application to Synthetic JWST Data

Abstract Spectroscopic eclipse observations, like those possible with the James Webb Space Telescope, should enable 3D mapping of exoplanet day sides. However, fully flexible 3D planet models are overly complex for the data and computationally infeasible for data-fitting purposes. Here, we present ThERESA, a method to retrieve the 3D thermal structure of an exoplanet from eclipse observations by first retrieving 2D thermal maps at each wavelength and then placing them vertically in the atmosphere. This approach allows the 3D model to include complex thermal structures with a manageable number of parameters, hastening fit convergence and limiting overfitting. An analysis runs in a matter of days. We enforce consistency of the 3D model by comparing the vertical placement of the 2D maps with their corresponding contribution functions. To test this approach, we generated a synthetic JWST NIRISS-like observation of a single hot-Jupiter eclipse using a global circulation model of WASP-76b and retrieved its 3D thermal structure. We find that a model that places the 2D maps at different depths depending on latitude and longitude is preferred over a model with a single pressure for each 2D map, indicating that ThERESA is able to retrieve 3D atmospheric structure from JWST observations. We successfully recover the temperatures of the planet’s day side, the eastward shift of its hot spot, and the thermal inversion. ThERESA is open source and publicly available as a tool for the community.

Toward a multidimensional analysis of transmission spectroscopy

Considering the relatively high precision that will be reached by future observatories, it has recently become clear that one dimensional (1D) atmospheric models, in which the atmospheric temperature and composition of a planet are considered to vary only in the vertical, will be unable to represent exoplanetary transmission spectra with a sufficient accuracy. This is particularly true for warm to (ultra-) hot exoplanets because the atmosphere is unable to redistribute all the energy deposited on the dayside, creating a strong thermal and often compositional dichotomy on the planet. This situation is exacerbated by transmission spectroscopy, which probes the terminator region. This is the most heterogeneous region of the atmosphere. However, if being able to compute transmission spectra from 3D atmospheric structures (from a global climate model, e.g.) is necessary to predict realistic observables, it is too computationally expensive to be used in a data inversion framework. For this reason, there is a need for a medium-complexity 2D approach that captures the most salient features of the 3D model in a sufficiently fast implementation. With this in mind, we present a new open-source documented version of Pytmosph3R that handles the computation of transmission spectra for atmospheres with up to three spatial dimensions and can account for time variability. Taking the example of an ultrahot Jupiter, we illustrate how the changing orientation of the planet during the transit can allow us to probe the horizontal variations in the atmosphere. We further implement our algorithm in TauREx to allow the community to perform 2D retrievals. We describe our extensive cross-validation benchmarks and discuss the accuracy and numerical performance of each model.

A Significant Increase in Detection of High-resolution Emission Spectra Using a Three-dimensional Atmospheric Model of a Hot Jupiter

Abstract High-resolution spectroscopy has opened the way for new, detailed study of exoplanet atmospheres. There is evidence that this technique can be sensitive to the complex, three-dimensional (3D) atmospheric structure of these planets. In this work, we perform cross-correlation analysis of high-resolution (R ∼ 100,000) CRIRES/VLT emission spectra of the hot Jupiter HD 209458b. We generate template emission spectra from a 3D atmospheric circulation model of the planet, accounting for temperature structure and atmospheric motions—winds and planetary rotation—missed by spectra calculated from one-dimensional models. In this first-of-its-kind analysis, we find that using template spectra generated from a 3D model produces a more significant detection (6.9σ) of the planet’s signal than any of the hundreds of one-dimensional models we tested (maximum of 5.1σ). We recover the planet’s thermal emission, its orbital motion, and the presence of CO in its atmosphere at high significance. Additionally, we analyzed the relative influences of 3D temperature and chemical structures in this improved detection, including the contributions from CO and H2O, as well as the role of atmospheric Doppler signatures from winds and rotation. This work shows that the hot Jupiter’s 3D atmospheric structure has a first-order influence on its emission spectra at high resolution and motivates the use of multidimensional atmospheric models in high-resolution spectral analysis.

Exoplanetary Monte Carlo radiative transfer with correlated-k – I. Benchmarking transit and emission observables

Abstract Current observational data of exoplanets are providing increasing detail of their 3D atmospheric structures. As characterization efforts expand in scope, the need to develop consistent 3D radiative-transfer methods becomes more pertinent as the complex atmospheric properties of exoplanets are required to be modelled together consistently. We aim to compare the transmission and emission spectra results of a 3D Monte Carlo radiative transfer (MCRT) model to contemporary radiative-transfer suites. We perform several benchmarking tests of an MCRT code, Cloudy Monte Carlo Radiative Transfer (cmcrt), to transmission and emission spectra model output. We add flexibility to the model through the use of k-distribution tables as input opacities. We present a hybrid MCRT and ray tracing methodology for the calculation of transmission spectra with a multiple scattering component. cmcrt compares well to the transmission spectra benchmarks at the 10s of ppm level. Emission spectra benchmarks are consistent to within 10 per cent of the 1D models. We suggest that differences in the benchmark results are likely caused by geometric effects between plane-parallel and spherical models. In a practical application, we post-process a cloudy 3D HD 189733b GCM model and compare to available observational data. Our results suggest the core methodology and algorithms of cmcrt produce consistent results to contemporary radiative transfer suites. 3D MCRT methods are highly suitable for detailed post-processing of cloudy and non-cloudy 1D and 3D exoplanet atmosphere simulations in instances where atmospheric inhomogeneities, significant limb effects/geometry or multiple scattering components are important considerations.

Effects of a fully 3D atmospheric structure on exoplanet transmission spectra: retrieval biases due to day–night temperature gradients

Transmission spectroscopy provides us with information on the atmospheric properties at the limb, which is often intuitively assumed to be a narrow annulus around the planet. Consequently, studies have focused on the effect of atmospheric horizontal heterogeneities along the limb. Here we demonstrate that the region probed in transmission – the limb – actually extends significantly towards the day and night sides of the planet. We show that the strong day–night thermal and compositional gradients expected on synchronous exoplanets create sufficient heterogeneities across the limb that result in important systematic effects on the spectrum and bias its interpretation. To quantify these effects, we developed a 3D radiative-transfer model able to generate transmission spectra of atmospheres based on 3D atmospheric structures. We first apply this tool to a simulation of the atmosphere of GJ 1214 b to produce synthetic JWST observations and show that producing a spectrum using only atmospheric columns at the terminator results in errors greater than expected noise. This demonstrates the necessity for a real 3D approach to model data for such precise observatories. Secondly, we investigate how day–night temperature gradients cause a systematic bias in retrieval analysis performed with 1D forward models. For that purpose we synthesise a large set of forward spectra for prototypical HD 209458 b- and GJ 1214 b-type planets varying the temperatures of the day and night sides as well as the width of the transition region. We then perform typical retrieval analyses and compare the retrieved parameters to the ground truth of the input model. This study reveals systematic biases on the retrieved temperature (found to be higher than the terminator temperature) and abundances. This is due to the fact that the hotter dayside is more extended vertically and screens the nightside – a result of the non-linear properties of atmospheric transmission. These biases will be difficult to detect as the 1D profiles used in the retrieval procedure are found to provide an excellent match to the observed spectra based on standard fitting criteria. This must be kept in mind when interpreting current and future data.

Simulating the cloudy atmospheres of HD 209458 b and HD 189733 b with the 3D Met Office Unified Model

Aims.To understand and compare the 3D atmospheric structure of HD 209458 b and HD 189733 b, focusing on the formation and distribution of cloud particles, as well as their feedback on the dynamics and thermal profile.Methods.We coupled the 3D Met Office Unified Model (UM), including detailed treatments of atmospheric radiative transfer and dynamics, to a kinetic cloud formation scheme. The resulting model self-consistently solves for the formation of condensation seeds, surface growth and evaporation, gravitational settling and advection, cloud radiative feedback via absorption, and crucially, scattering. We used fluxes directly obtained from the UM to produce synthetic spectral energy distributions and phase curves.Results.Our simulations show extensive cloud formation in both HD 209458 b and HD 189733 b. However, cooler temperatures in the latter result in higher cloud particle number densities. Large particles, reaching 1μm in diameter, can form due to high particle growth velocities, and sub-μm particles are suspended by vertical flows leading to extensive upper-atmosphere cloud cover. A combination of meridional advection and efficient cloud formation in cooler high latitude regions, results in enhanced cloud coverage for latitudes above 30° and leads to a zonally banded structure for all our simulations. The cloud bands extend around the entire planet, for HD 209458 b and HD 189733 b, as the temperatures, even on the day side, remain below the condensation temperature of silicates and oxides. Therefore, the simulated optical phase curve for HD 209458 b shows no “offset”, in contrast to observations. Efficient scattering of stellar irradiation by cloud particles results in a local maximum cooling of up to 250 K in the upper atmosphere, and an advection-driven fluctuating cloud opacity causes temporal variability in the thermal emission. The inclusion of this fundamental cloud-atmosphere radiative feedback leads to significant differences with approaches neglecting these physical elements, which have been employed to interpret observations and determine thermal profiles for these planets. This suggests that readers should be cautious of interpretations neglecting such cloud feedback and scattering, and that the subject merits further study.

Reflected light from 3D exoplanetary atmospheres and simulation of HD 209458b

We present radiation transfer models that demonstrate that reflected light levels from 3D exoplanetary atmospheres can be more than 50 per cent lower than those predicted by models of homogeneous or smooth atmospheres. Compared to smooth models, 3D atmospheres enable starlight to penetrate to larger depths resulting in a decreased probability for the photons to scatter back out of the atmosphere before being absorbed. The increased depth of penetration of starlight in a 3D medium is a well-known result from theoretical studies of molecular clouds and planetary atmospheres. For the first time we study the reflectivity of 3D atmospheres as a possible explanation for the apparent low geometric albedos inferred for extrasolar planetary atmospheres. Our models indicate that 3D atmospheric structure may be an important contributing factor to the non-detections of scattered light from exoplanetary atmospheres. We investigate the self-shadowing radiation transfer effects of patchy cloud cover in 3D scattered light simulations of the atmosphere of HD 209458b. We find that, for a generic planet, geometric albedos can be as high as 0.45 in some limited situations, but that in general the geometric albedo is much lower. We conclude with some explanations on why extrasolar planets are likely dark at optical wavelengths.

3D water-vapor tomography with Shanghai GPS network to improve forecasted moisture field

The vertical structure of water vapor in atmosphere is one of the initial information of numerical weather forecast model. Because of the strong variation of water vapor in atmosphere and limited spatio-temporal solutions of traditional observation technique, the initial water vapor field of numerical weather forecast model can not accurately be described. At present, using GIPS slant observations to study water vapor profile is very popular in the world. Using slant water vapor(SWV) observations from Shanghai GIPS network,we diagnose the three-dimensional(3D) water vapor structure over Shanghai area firstly in China. In water vapor tomography, Gauss weighted function is used as horizontal constraint, the output of numerical forecast is used as apriori information, and boundary condition is also considered. For the problem without exact apriori weights for observations, estimation of variance components is introduced firstly in water vapor tomography to determine posteriori weights. Robust estimation is chosen for reducing the effect of blunders on solutions. For the descending characteristic of water vapor with height increasing, non-equal weights are used along vertical direction. Comparisons between tomography results and the profile provided by numerical model (MM5) show that the forecasted moisture fields of MM5 can be improved obviously by GPS slant water vapor. Using GIPS slant observations to study 3D structure of atmosphere in near real-time is very important for improving initial water vapor field of short-term weather forecast and enhancing the accuracy of numerical weather forecast.

A model for atmospheric background radiance structures

Infrared radiance fluctuations in the atmosphere result from fluctuations in the density of atmospheric species, individual molecular state populations, and kinetic temperatures along the sensor line of sight (LOS). As part of the SHARC-4 development program, a model for the prediction of atmospheric background radiance fluctuations was constructed. It predicts a two dimensional radiance spatial covariance function and power spectral density (PSD) from the underlying 3D atmospheric structures. Inputs to the model include the parameters of an atmospheric temperature fluctuation model and vertical atmospheric profiles. The model can be applied for arbitrary sensor viewing geometries, including limb viewing. In the upper atmosphere non-equilibrium effects are important. Fluctuations in kinetic temperature can result in correlated or anti-correlated fluctuations in vibrational state temperatures. The model accounts for these effects and predicts spatial covariance functions and PSD's for molecular state number densities and vibrational temperatures. SHARC is used to predict the non-equilibrium dependence of molecular state number density fluctuations on kinetic temperature and density fluctuations, and to calculate mean LOS radiances and radiance derivatives. The goal is to understand the origin of LOS radiance fluctuations and the prediction of fluctuation statistics based on local fluctuations in gas kinetic temperatures, and densities. The modeling capabilities are illustrated with a sample image prediction for the Midcourse Space Experiment (MSX) having an MSX sensor bandpass and field-of-view (FOV).