Abstract
Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (that is, elements heavier than helium, also called âmetallicityâ)1,2,3, and thus the formation processes of the primary atmospheres of hot gas giants4,5,6. It is also one of the most promising species to detect in the secondary atmospheres of terrestrial exoplanets7,8,9. Previous photometric measurements of transiting planets with the Spitzer Space Telescope have given hints of the presence of CO2, but have not yielded definitive detections owing to the lack of unambiguous spectroscopic identification10,11,12. Here we present the detection of CO2 in the atmosphere of the gas giant exoplanet WASP-39b from transmission spectroscopy observations obtained with JWST as part of the Early Release Science programme13,14. The data used in this study span 3.0â5.5âmicrometres in wavelength and show a prominent CO2 absorption feature at 4.3âmicrometres (26-sigma significance). The overall spectrum is well matched by one-dimensional, ten-times solar metallicity models that assume radiativeâconvectiveâthermochemical equilibrium and have moderate cloud opacity. These models predict that the atmosphere should have water, carbon monoxide and hydrogen sulfide in addition to CO2, but little methane. Furthermore, we also tentatively detect a small absorption feature near 4.0âmicrometres that is not reproduced by these models.
Similar content being viewed by others
Main
WASP-39b is a hot (planetary equilibrium temperature of 1,170âK assuming zero albedo and full heat redistribution), transiting exoplanet that orbits a G7-type star with a period of 4.055âdays15. The planet has approximately the same mass as Saturn (Mâ=â0.28âMJ, where MJ is the mass of Jupiter) but is about 50% larger (radius Râ=â1.28âRJ, where RJ is the radius of Jupiter), probably owing to the high level of irradiation that it receives from its host star16,17,18. We chose this planet for the JWST Early Release Science (ERS) transmission spectroscopy observations because analyses of existing space- and ground-based data detected large spectral features and showed that there was minimal contamination of the planetary signal from stellar activity10,19,20,21. The main spectral features previously detected were confidently attributed to sodium, potassium and water vapour absorption10,19,20, whereas carbon dioxide (CO2) was suggested to explain the deep transit at 4.5âµm seen with Spitzer10.
Atmospheric metallicity has long been thought to be a diagnostic of the relative accretion of solids and gas during the formation of gas giant planets, both of which bring heavy elements to the hydrogen-dominated envelope and visible atmosphere4,5,6. The metallicity of WASP-39bâs host star, which is a proxy for the metal enrichment of the protoplanetary disk that the planet formed in, is approximately solar15,22,23,24. Therefore, the planet massâatmospheric metallicity trend observed in the Solar System giants25,26 predicts that it has an enhancement of about ten-times solar (like that of Saturn; ref. 27). In addition, interior structure models that match WASP-39bâs low density predict a 95th percentile upper limit for the atmospheric metallicity of 55-times solar, under the limiting assumption that the planet has no heavy-element core and that all the metals are evenly distributed throughout the envelope28.
Despite having some of the highest signal-to-noise detections of spectral features in its transmission spectrum, modelling of the existing data for WASP-39b has resulted in metallicity estimates ranging across five orders of magnitude, from 0.003-times solar to 300-times solar10,29,30,31,32,33. The wide range of values stems from the data being of insufficient quality to break the degeneracy between clouds and metallicity in transmission spectra models34, as well as uncertainty over the interpretation of the photometric measurements by the Spitzer Space Telescope at 3.6âµm and 4.5âµm. Thus, spectroscopic data with greater precision, finer spectral channels and wider wavelength coverage were needed to better constrain the metallicity of this (and other) giant exoplanet atmospheres.
The first JWST ERS observation of WASP-39b was obtained using the Near Infrared Spectrograph (NIRSpec)35,36 on 10 July 2022, between 15:24 and 23:37âutc. We used the Bright Object Time Series (BOTS) mode with the 1.6â³âÃâ1.6â³ fixed-slit aperture and the PRISM disperser to capture spectra between 0.5âµm and 5.5âµm. The data were recorded using the SUB512 subarray with five groups per integration and the NRSRAPID readout pattern, which gave integration times of 1.38âs. NIRSpec obtained a total of 21,500 integrations over 8.23âh of observations centred on the 2.8-h transit duration of WASP-39b.
The count rate in the PRISM mode varies significantly over the bandpass owing to the spectral energy distribution of the star and the wavelength dependency of the spectrograph dispersion. Therefore, the observations were designed to saturate at shorter wavelengths in order to obtain sufficient signal-to-noise ratio at the longer wavelengths in the bandpass that have not previously been studied spectroscopically. Wavelengths between 0.71âµm and 2.09âµm have at least one group saturated in the pixel at the centre of the spectral trace. We concentrate here on the analysis of the data longwards of 3.0âµm that are not impacted by saturation to investigate the spectrum overlapping with the previous 3.6âµm and 4.5âµm Spitzer photometric measurements. The subset of the PRISM data described herein has a native spectral resolving power (Râ=âλ/Îλ, where λ is wavelength) of 100â350. For this study, we binned the data to lower resolving powers (values range from 60 to 200 depending on wavelength and reduction). The binning is done at the light-curve level before the fitting of the transit depths that constitute the transmission spectrum. Analyses of JWST/NIRSpec transit observations obtained during commissioning have shown that similar levels of binning as we use here results in minimal systematics37. An analysis of the complete PRISM dataset at full resolution, including recovery of the saturated part of the spectrum, is ongoing.
We reduced the NIRSpec PRISM data for WASP-39b using the JWST Science Calibration Pipeline along with customized routines to minimize noise in the time-series spectra (Methods). We performed four different reductions of the transmission spectrum starting from the uncalibrated data21,38,39,40. Figure 1 shows derived spectroscopic transit light curves from one of the reductions. We confirm with our analysis of the WASP-39b data that NIRSpec transit observations at a resolving power of 60â200 are nearly free of systematics. We achieved close to photon-noise-limited measurements in the spectroscopic light curves after trimming the first 10âmin of data and removing a linear trend in time with an average rate of about 190âppmâhâ1 across the bandpass. We also obtained similar results by fitting the full time series with a downwards trending exponential ramp (timescale about 100âmin) combined with a quadratic function of time. The lack of large systematics in these data stands in contrast to previous transit spectroscopy observations with space- or ground-based telescopes41.
a, Spectroscopic light curves for WASP-39bâs transit with a spectral resolving power of 20 and a time cadence of 1âmin (data are binned and offset vertically for display purposes only). An exoplanet light-curve model was fitted to the data using a quadratic limb-darkening law with an exponential ramp and a quadratic function of time removed. b, Residuals of the binned light curve after subtracting the transit model scaled up by a factor of five to show the structure. The r.m.s. of the residuals are given in units of ppm. The numbers in brackets are the ratio of the r.m.s. to the predicted photon-limited noise.
The transmission spectra derived from the different reductions, shown in Fig. 2, have excellent agreement. They all show a large feature at 4.3âµm, as well as a smaller feature near 4.0âµm (discussed below). Detailed modelling of the Fast InfraRed Exoplanet Fitting Lyghtcurve (FIREFLy)-reduced data yields a statistical significance of 26Ï for the large feature (Methods). We attribute this feature to CO2 absorption based on a comparison of the resolved band shape with theoretical models and the spectra of brown dwarfs42. Figure 2 also includes Spitzerâs two broadband photometric measurements10, which are consistent with the JWST data to better than 2Ï after integrating the transmission spectrum over the Spitzer bandpasses. We also see good agreement (better than 2Ï for all reductions) in the relative transit depths between the 3.6-µm and 4.5-µm channels. The comparison shown in Fig. 2 demonstrates both the consistency in the derived spectra from multiple, independent analyses and the reliability of the previous Spitzer measurements.
The JWST data (small coloured points) are compared with Spitzerâs two Infrared Array Camera (IRAC) broadband photometric measurements (grey circles and corresponding sensitivity curves labelled IRAC1 and IRAC2). The axis on the right shows equivalent scale heights (750â1,000âkm) in WASP-39bâs atmosphere; for plotting purposes, we assume that one scale height corresponds to 800âkm. The JWST data are consistent with the Spitzer points (within 2Ï) when integrated over the broad bandpasses (indicated by the horizontal lines). The relative transit depths between the 3.6-µm and 4.5-µm channels are also consistent within 2Ï between independent reductions of the JWST data, with most of the deviation coming from the 3.6-µm bandpass. Vertical error bars indicate 1Ï uncertainties.
We compared the data with a suite of one-dimensional atmospheric structure and transmission spectrum models to constrain the composition of WASP-39bâs atmosphere. These models assume radiativeâconvectiveâthermochemical equilibrium, and they adopt a scaled solar abundance pattern. We calculated planet-specific grids of these models over a range of atmospheric metallicities, carbon-to-oxygen ratios and cloud properties using four different codes. These grids of self-consistent model transmission spectra were then fitted to the FIREFLy-reduced data (the fit results are independent of which dataset we use) while also adjusting for a reference radius at 1âbar. The results are illustrated in Fig. 3; see Methods for further details.
Top: a comparison of the FIREFLy reduction and its 1Ï uncertainties (labelled âDataâ) to the best-fit ScCHIMERA theoretical model binned to the resolution of the data (blue curve; Methods). The key parameters of the model are 10-times solar metallicity, a carbon-to-oxygen ratio of 0.35 and cloud opacity of 7âÃâ10â3âcm2âgâ1. The impact of the opacity sources expected from thermochemical equilibrium over the full bandpass are indicated by removing the opacity contribution from individual gases one at a time. As in Fig. 2, the axis on the right shows equivalent scale heights in WASP-39bâs atmosphere. Bottom: the molecular absorption cross-sections for each gas in the best-fit model. The model is well matched to the data (Χ2/Ndataâ=â1.3), suggesting that our assumptions broadly capture the important physics and chemistry in WASP-39bâs atmosphere. However, there is a feature near 4.0âµm that cannot be reproduced by the models used here. The strong CO2 absorption (4.1â4.6âµm) and the apparent lack of methane (3.0â3.5âµm) is what drives the solution to an elevated atmospheric metal enrichment, ruling out previous low-metallicity estimates29,30,31. The other reductions and models give similar results.
Under similar assumptions, all four model grids are able to match the dominant spectral morphologiesânamely the strong CO2 feature between 4.1âµm and 4.6âµm and the rise in transit depth bluewards of 3.6âµm owing to water (H2O) vapour (a species that had been detected previously at shorter wavelengths10). More subtle modulations over the whole bandpass are potentially owing to contributions from clouds, carbon monoxide (CO) and hydrogen sulfide (H2S), although the degree to which the two gas species contribute is unknown pending further study.
Several models for warm gas giant atmospheres predict that the CO2 abundance scales quadratically with atmospheric metallicity, becoming detectable at 4.3âµm for metallicities above that of the Sun1,2,3. The representative best-fit model shown in Fig. 3 is consistent with this scenario. It has a 10-times-solar metal enrichment and a slightly subsolar carbon-to-oxygen ratio (0.35, compared with the solar value of 0.55; ref. 43). The moderate contribution of cloud opacity predicted by the best-fit model is consistent with interpretations of previous population-level studies of planets that have similar temperatures and gravities to WASP-39b44,45. It is also consistent with the predictions of aerosol microphysics and global circulation models of hot giant planets46,47,48.
In addition to the large CO2 feature, we also identify a smaller spectral feature near 4.0âµm that is not matched by our thermochemical equilibrium models (Fig. 3). This feature is present in all four independent reductions and has a significance of 2Ï (Methods). Further data analysis and modelling including non-equilibrium chemistry are needed to fully assess the robustness of this feature and to identify the chemical species that gives rise to it. Additional JWST ERS observations of WASP-39b that will use the G395H grating on NIRSpec also have the potential to confirm the 4.0-µm feature and resolve it in greater detail.
The grid fits explored here favour lower metallicities than refs. 10,21, and higher metallicities than ref. 31, even though the Spitzer data that their studies included are consistent with our JWST data. The higher precision and more resolved measurement of the CO2 feature enabled by JWST pulls the models of refs. 10,21 to lower metallicity and increased cloudiness. Nevertheless, it is not possible to obtain a robust confidence interval on this inference without more rigorous Bayesian analyses, which is left to future work (Methods). Continued modelling of WASP-39b will also be aided by the future measurements of the planetâs transmission spectrum from 0.5âµm to 5.5âµm that are also being obtained by this ERS programme. The final transmission spectrum will ultimately have higher spectral resolution than the data presented here (more than four times over most of the bandpass), and will be validated using multiple JWST instruments.
Methods
Data reduction
We reduced the JWST NIRSpec PRISM data for WASP-39b using four separate pipelines to confirm that the results did not depend on the specifics of the analyses, as was sometimes the case for results from the Spitzer Space Telescope (for example, ref. 49). The descriptions below refer to calibration pipelines and other software whose code and citations appear in âCode availabilityâ.
tshirt pipeline
We used the Time Series Helper and Integration Reduction Tool40 (tshirt) to extract light curves of the spectrum. This pipeline modifies the JWST Calibration pipeline steps to improve the precision of the reduction. tshirt has been used to successfully analyse the JWST transit observations of HAT-P-14b that were obtained during commissioning with the Near Infrared Camera (NIRCam)37. First, we used an updated bias frame from commissioning programme 1130 observation 29 and ran the JWST Calibration pipeline until the reference pixels step. We then applied a correction for 1/f noise (so named since the noise power is inversely proportional to the signal frequency, f), which varies for odd and even rows and for each column. We use background pixels for the calibration as reference pixels are not available in this subarray. We skipped the jump and dark subtraction steps because they were seen to add noise to the light curves. tshirt fits the profile of the spectrum with splines and rejects outlier pixels that are more than 50Ï from the spline fits. We used covariance-weighted extraction50 with an assumed pixel correlation of 0.08. For spectral extraction, we used a background region no closer than 7âpixels on either side of the source and an extraction region width of 16âpixels. The scatter in the light curve was consistent with the theoretical limit of photon and read noise over short timescales.
We fit the light curves with a second-order (quadratic) polynomial baseline, uninformative quadratic limb-darkening priors and an exponential start-up ramp with 10Ï clipping of outliers. To begin, we fit the white-light curve with priors on the transit centre, inclination and period from ref. 22. We also used the ratio a/R* (where a is the semi-major axis and R* is the stellar radius) from ref. 22 but widened the uncertainty on this parameter because the enforced prior resulted in significant residuals. Next, we fit each spectroscopic light curve individually with the orbital parameters fixed at the value from the white-light posterior medians. We modelled the light curves using the âexoplanetâ code51 and the pymc352 sampler. We evaluated the wavelengths using the JWST Calibration pipeline at pixel row 16 (Yâ=â16) from the world coordinate solution. This uses an instrument model and could not be verified owing to a lack of strong stellar absorption features at the NIRSpec resolution. All the other reductions adopted this wavelength calibration. As shown in Fig. 1, the standard deviation in the out-of-transit light curve approaches the theoretical limit of photon and read noise at short wavelengths, but is 20% to 40% higher at longer wavelengths, which may be related to uncorrected 1/f noise.
Eureka! pipeline
Eureka!39 is a data reduction and analysis pipeline for time-series observations with the JWST or the Hubble Space Telescope Its modular, multi-stage design provides flexibility and ease of comparison at any step, starting from uncalibrated FITS files and resulting in precise transmission or emission spectra. Eureka! has been used to successfully analyse the JWST transit observations of HAT-P-14b that were obtained during commissioning with NIRCam37.
We began the data reduction process using the uncalibrated raw data files with the âuncalâ suffix available from the Mikulski Archive for Space Telescopes (MAST) archive. The first stage of the Eureka! pipeline is primarily a wrapper for Stage 1 of the JWST Calibration pipeline, which converts groups to slopes. For this dataset, we skipped the jump detection step as it led to a large fraction of detector pixels being incorrectly flagged as outliers. We did, however, search for and flag outliers at multiple points in subsequent stages. We also manually updated the bad-pixel map to include identified hot pixels on the detector that were not provided in the current (July 2022) full-detector STScI data-quality map. As part of Eureka!, we performed a custom background subtraction at the group level before Stage 1 ramp fitting to account for 1/f noise introduced during detector readout. We set the top and bottom six rows of the detector as our background region and flagged pixels deemed outliers at >3Ï. We then subtracted the mean flux per pixel column and repeated this for each group and integration in the observation. Similarly to Stage 1, the second stage of the Eureka! pipeline is a wrapper for Stage 2 of the JWST Calibration pipeline, which calibrates the two-dimensional time series of fitted slopes. Here, we skipped the flux calibration step, thus leaving the data in units of digital number (DN) per second (DNâsâ1).
For Stage 3, we performed background subtraction and optimal extraction of the stellar spectrum for each integration with Eureka!. We used only pixels 14 to 495 in the dispersion direction of the 512âÃâ32-pixel subarray, as NIRSpecâs throughput is negligible beyond this range. We also masked pixels that have a non-zero data-quality flag to avoid any impact of outlier pixels on the extracted spectra or background subtraction. The position of the source on the detector along the cross-dispersion dimension is located by fitting a Gaussian to the pixel values summed over all detector columns. For each pixel, we examined its flux variation in time and performed a double-iteration, 10Ï-outlier rejection test. We then executed a second column-by-column background subtraction, this time at the integration level, using pixels located at least 8âpixels away from the source position to compute the mean background per column. Performing this additional background subtraction reduced the number of outliers in the measured light curves and accounted for the residual background and/or noise introduced during the ramp fitting procedure. As with Stage 1, we exclude 3Ï outliers from our background region. We adopted an aperture half-width of 7âpixels for our optimal spectral extraction step, constructing the profile from the median frame. At the end of this stage, we obtained a time series of one-dimensional spectra.
For the remaining stages, we used multiple pipelines (Eureka!39 and ExoTEP53,54,55) to generate and fit the light curves. We first generated median-normalized light curves at the instrumentâs native resolution (that is, from each detector column) using our Stage 3 outputs. We then clipped additional outliers in time for the white and spectroscopic light curves. For this step, we first rejected integrations that were more than 3Ï outliers for the source position in the cross-dispersion direction, the width of the fitted Gaussian to the spatial profile or the drift in the dispersion direction. Next, we produced a median-filtered version of the light curve and clipped out 3Ï outliers in flux. We jointly fit astrophysical and systematics model parameters to the white and individual spectroscopic light curves. Our astrophysical transit model used the batman package56 with uniform priors, fitting for the following astrophysical parameters: the two coefficients of a stellar quadratic limb-darkening law, impact parameter, semi-major axis, transit time and the planet-to-star radius ratio (Rp/R*) in each of the wavelength channels. Although the limb-darkening coefficients and planet-to-star radius ratio were fit independently in each spectroscopic channel, we used the best-fitting value of the planetâs impact parameter, semi-major axis and transit time from a white-light curve fit as a fixed value in the wavelength-dependent fits. For the systematics model, we assumed a linear trend in time for each wavelength channel, fitting for both the slope and y intercept. Last, we fit a single-point scatter to each light curve, which illustrates the level of additional noise required for our joint model to reach a reduced chi-squared (Χ2) of unity. The white-light curve residuals have a root mean square (r.m.s.) of 3,013âppm, and the spectroscopic light curves above 3âµm have a median r.m.s. of 5,779âppm. Similar to the reduction shown in Fig. 1, both pipelines reach near photon noise. The Eureka! and ExoTEP transmission spectra appear nearly identical; therefore, only one (Eureka!) is shown in Fig. 2.
Tiberius pipeline
We built on the pipeline developed for the analysis of the Low Resolution Ground-Based Exoplanet Survey using Transmission Spectroscopy (LRG-BEASTS) data21,57,58 to provide an independent reduction of the data. We began with the outputs of the JWST Calibration Stage 1 pipeline with the jump step correction turned off. We created bad-pixel and cosmic-ray masks by identifying 5Ï outliers in running medians operating along pixel rows and along individual pixels in time. Before tracing the spectra, we interpolated each column of the detector onto a finer grid, ten times the initial spatial resolution, to improve the extraction of flux at the subpixel level. We used a fourth-order polynomial to trace the spectra and a four-pixel-wide aperture. To remove the 1/f noise, we fit a linear polynomial to 21 background pixels along each column in the cross-dispersion direction. Next, to correct for shifts in the dispersion direction, we cross-correlated each stellar spectrum with the first spectrum of the observation to account for very small (0.003â0.005) subpixel shifts. Our white-light curve spans a wavelength range of 0.518â5.348âμm after masking saturated pixels, and our 147 spectroscopic light curves used 3-pixel-wide bins across this same wavelength range. We masked frames 20751â20765 owing to a high-gain-antenna move that led to increased noise in the light curves.
We fit our light curves with a combination of a quadratically limb-darkened transit model (through batman56) with a linear-in-time polynomial. We began by fitting the white-light curve to derive the system parameters: inclination, i, time of mid-transit, TC, the semi-major axis scaled to the stellar radius, a/R*, and the linear limb-darkening coefficient, u1. We placed wide boundaries on the parameter values only to prevent unphysical values. In practice, the parameter values did not get close to the boundaries. We fixed the planetâs orbital period to 4.0552941âd and the eccentricity to 0 from ref. 22. We fixed the quadratic coefficient, u2, to theoretical values determined by ExoTiC-LD59,60 with three-dimensional stellar models61, and fit for u1. We used a LevenbergâMarquardt algorithm to fit our light curves, rescaled our photometric uncertainties to give a reduced Χ2â=â1 for our best-fit model and then re-ran the fits. For the spectroscopic light curves, the system parameters (i, TC and a/R*) were held fixed to the best-fit values found from the white-light curve. The white-light curve residuals had an r.m.s. of 2,761âppm and the spectroscopic light curve residuals had a median r.m.s. of 6,731âppm. In both cases, the variance of the residuals scales upon binning as expected for Poisson noise.
FIREFLy pipeline
We also reduced the data using the FIREFLy reduction routines38. These routines utilize the JWST Calibration pipeline with custom modifications. This pipeline has been used to successfully analyse the JWST transit observations of HAT-P-14b that were obtained during commissioning with NIRSpec G39537. We removed 1/f noise (see ref. 36) at the group level, as the 1/f noise changes from group to group. We also skipped the jump step and instead flagged and removed cosmic rays, bad pixels, hot pixels and other outliers using median filtering of the data both spatially and in time, flagging pixels using a 5Ï-outlier threshold algorithm. The time series of two-dimensional spectra were aligned using cross-correlation and interpolation, with the time-series spectra exhibiting an r.m.s. jitter of 0.005âpixels in the x-axis direction and 0.0026âpixels in the y-axis direction. We found a small inverse ramp in the light curves, which settled down after the first 2,000 exposures, which we discarded. We fit the light curves with the batman56 transit model along with a linear baseline and a second-order jitter detrending polynomial of x and y detector position as described by ref. 38, which are present in the spectrophotometry at the 53â±â2-ppm level in the x direction and at the 140â±â3-ppm level in the y direction. We applied a fixed quadratic limb-darkening law using the three-dimensional models61 computed using the methods of ref. 62 from ExoTiC-LD59,60. In fitting the 3âμm to 5.5âμm white-light curve, we allowed the semi-major axis in units of stellar radii a/R*, inclination i and central transit time T0 to freely vary along with the transit depth and systematics model. We used the Markov chain Monte Carlo sampling routine emcee63 to find the best-fit parameters and measure the posterior distribution. We find the 3â5.5âμm white-light curve has a transit depth of 2.1368â±â0.0014% and achieves 808-ppm scatter in the residuals. This is within 6% of the expected noise limit of 758âppm as calculated by the JWST Calibration pipeline, with the scatter of the residuals decreasing to below 40âppm upon binning with no detectable red noise. We fit each spectroscopic light curve shown in Fig. 1 with the same astrophysical and systematic models as the white-light curve, except fixing the system parameters (a/R*, i and T0). The transmission spectral light-curve residuals for each bin are typically within 5% of pipeline error or better, also with no detectable red noise.
Dataâmodel comparison
We compared the extracted transmission spectral data to a suite of one-dimensional self-consistent radiativeâconvectiveâthermochemical equilibrium model atmospheres (see, for example, refs. 64,65 for a general description of such models) described below. In short, all models are able to fit the 3â5.5âμm spectra consistently (with Χ2/Ndataâ<â1.4, where Ndata is the number of spectral data points) with a 10-times-solar metal enrichment and varying grey cloud opacity for their single best estimate. Comparisons of the model fits from each grid are shown in Extended Data Fig. 1. For additional parameters within the grid (for example, carbon-to-oxygen ratio (C/O) and heat redistribution), there is some discrepancy between each model gridâs single best estimate values. Additional Bayesian analyses are needed to rigorously quantify confidence intervals on atmospheric properties of interest, which is beyond the scope of this work. Future works will focus on modelling that includes the effects of disequilibrium chemistry, aerosol microphysics and three-dimensional circulation effects. We assumed the following parameters in the modelling: stellar effective temperature, Teffâ=â5,512âK, stellar radiusâ=â0.932âRâ, planet massâ=â0.281âMJ, planet radiusâ=â1.279âRJ and planet orbital semi-major axisâ=â0.04828âau.
ScCHIMERA
This framework was first described in refs. 66,67, with the most recent updates, methods and opacity sources described in ref. 68. We compute the converged atmospheric structure (temperatureâpressure and thermochemical equilibrium gas mixing ratio profiles) over a grid of atmospheric metallicity ([M/H], where the square brackets indicate log10 enrichment relative to solar43) spaced at 0.25-dex intervals between 0 and 2.25 (1- to 175-times solar) and C/O at values of 0.20, 0.35, 0.55, 0.70, 0.75 and 0.80. We assume full day-to-night temperature redistribution69 as planets in this temperature regime are unlikely to possess strong day-to-night temperature contrast70,71. We then compute transmission spectra72,73 from these converged atmospheric structures. To match the models to the data, the DYNESTY74 fitting routine is used to search for the optimal [M/H] and C/O (via nearest neighbour) while simultaneously adjusting the 1-bar planetary radius (which controls the absolute transit depth) and an opaque, grey, uniformly vertically distributed, cloud opacity (κcld). The optimal model resulting from this process is [M/H]â=â+1.0, C/Oâ=â0.35 and log10κcldâ=ââ2.15âcm2âgâ1. The metallicity and cloud opacity are primarily driven by the strength of the 4.3-µm CO2 feature and lack of methane (CH4) absorption near 3.3âµm. This result is what is shown in the main text (Fig. 3), which also illustrates the relative contribution of the key opacity sources (H2O (refs. 75,76), CO (refs. 77,78), CO2 (refs. 79,80), H2S (refs. 78,81) and CH4 (refs. 78,82)) to the overall spectral shape. Extended Data Fig. 2 shows the atmospheric structure (temperature profile and gas mixing ratio profiles) for this best-fit model.
PICASO
The core one-dimensional radiativeâconvective model is based on the legacy âExtrasolar Giant Planetâ code described in refs. 69,80,83 and since updated and modernized within the PICASO84 framework described in ref. 85 (PICASO 3.0). The PICASO 3.0 model uses gaseous opacities created from the references listed in ref. 80. The grid of PICASO models contains metallicity points at 0.1-, 0.3-, 1-, 3-, 10-, 30-, 50- and 100-times solar; C/O at 0.23, 0.46, 0.69 and 0.92; and also assumes full dayânight heat redistribution. The clouds are modelled using the Virga86 implementation of the Eddysed87 framework, which requires a vertical mixing coefficient, Kzz (constant with altitude; log10Kzzâ=â5, 7, 9 and 11 (cgs units)) and a vertically constant sedimentation parameter (fsedâ=â0.6, 1, 3, 6 and 10), with optical/material properties for clouds thought to exist at WASP-39bâs pressures and temperatures (Na2S, MnS and MgSiO3). The fsed parameter controls the vertical extent of the cloud, and Kzz and fsed together control the mean droplet sizes with altitude in the atmosphere. A Χ2 grid search along the described dimensions is performed to identify the best fit. Within this grid, the nominal best fit (Χ2/Ndataâ=â1.34) is 10-times-solar metallicity, a subsolar C/O (0.23), with an extended large droplet cloud (fsedâ=â0.6, log10Kzzâ=â9) that produces a grey continuum over these wavelengths, consistent with the ScCHIMERA results above.
ATMO
The ATMO radiativeâconvectiveâthermochemical equilibrium solver is described in refs. 88,89,90,91. This grid consists of model transmission spectra for four different dayânight energy redistribution factors (0.25, 0.5, 0.75 and 1.0, where 0.5 is âfullâ and 1.0 is âdayside onlyâ), six metallicities (0.1-, 1-, 10-, 50-, 100- and 200-times solar), six C/O ratios (0.35, 0.55, 0.70, 0.75, 1.0 and 1.5), two haze factors (no haze and 10-times multi-gas Rayleigh scattering) and four grey cloud factors (no cloud, 0.5-, 1- and 5-times the strength of H2 Rayleigh scattering at 350ânm between 1-mbar and 50-mbar pressure levels). Each model transmission spectrum from the grid is binned to the same resolution as that of the observations to compute Χ2, with a (wavelength independent) transit depth offset as the free parameter. Within this grid, we find a best-fit model (Χ2/Ndataâ=â1.39) spectrum arising from a redistribution factor of 0.75 (slightly hotter than a full dayânight redistribution would produce), a metallicity of 10-times solar, a super-solar C/O ratio of 0.7, a haze factor of 10 and a cloud factor of 5.
PHOENIX
This model originates from the PHOENIX stellar atmosphere code92 adapted for exoplanets93 with additional modelling and opacity updates described in refs. 94,95. The model grid is computed for an array of irradiation temperatures (920âK, 1,020âK, 1,120âK and 1,220âK), metallicities (0.1-, 1-, 10- and 100-times solar) and C/O (0.3, 0.54, 0.7 and 1.0), and includes a sampling of opaque, grey clouds at specified cloud-top pressures. The nominal best-fit model (Χ2/Ndataâ=â1.32) from this grid set-up results in a 10-times-solar metallicity and subsolar C/O (0.3) atmosphere with a cloud-top pressure of 0.3âmbar.
Quantifying feature detection significance
We quantified the detection significance96 of CO2 with the following steps. The best-fit grid model without CO2 (that is, the âno CO2â black curve shown in Fig. 3) is first subtracted from the data, leaving behind a strong residual feature due to CO2 (Extended Data Fig. 3). The peak per-spectral-bin mean signal-to-noise ratio of this residual feature is about 10Ï. To utilize the full line/band shape we then fit the residual peak with (1) a four-parameter Gaussian model (centroid, amplitude, width and vertical offset), shown as red curves in Extended Data Fig. 3, and (2) a âno featureâ constant using a nested sampling routine74. The Bayesian evidence between the Gaussian model and constant model were then used to compute a Bayes factor, B, and corresponding detection significance97. For the CO2 residual feature, ln(B) is 340.5, which equates to a 26.2Ï detection. From this analysis, we conclude that the CO2 feature is robustly detected.
On inspecting Figs. 2 and 3, there appears to be a feature near 4.0âµm (just short of the major CO2 feature). We repeated the same analysis as above, but instead compared the Bayesian evidence from a two-component Gaussian model fit (to accommodate for both the CO2 feature and the unknown absorber) to that of the single component Gaussian model fit above. On doing so, we find ln(B)â=â0.98, which equates to a 2Ï significance. Restricting the prior range for the second Gaussian to be localized near the 4-µm feature boosts the significance to 2.3Ï. Future analyses will focus on the nature of this feature and more rigorous quantification via nested Bayesian model comparison within atmospheric retrieval frameworks (for example, ref. 34).
Data availability
The data used in this paper are associated with JWST programme ERS 1366 (observation #4) and are available from the Mikulski Archive for Space Telescopes (https://mast.stsci.edu). Science data processing version (SDP_VER) 2022_2a generated the uncalibrated data that we downloaded from MAST. We used JWST calibration software version (CAL_VER) 1.5.3 with modifications described in the text. We used calibration reference data from context (CRDS_CTX) 0916, except as noted in the text. All the data and models presented in this publication can be found at https://doi.org/10.5281/zenodo.6959427. Source data are provided with this paper.
Code availability
The codes used in this publication to extract, reduce and analyse the data are as follows: STScI JWST Calibration pipeline37 (https://github.com/spacetelescope/jwst), tshirt40, Eureka!39 (https://eurekadocs.readthedocs.io/en/latest/), Tiberius21,56,57 and FIREFLy38. In addition, these made use of Exoplanet51 (https://docs.exoplanet.codes/en/latest/), Pymc352 (https://docs.pymc.io/en/v3/index.html), ExoTEP53,54,55, Batman56 (http://lkreidberg.github.io/batman/docs/html/index.html), ExoTiC-ISM59 (https://github.com/Exo-TiC/ExoTiC-ISM), ExoTiC-LD60 (https://exotic-ld.readthedocs.io/en/latest/), Emcee63 (https://emcee.readthedocs.io/en/stable/), DYNESTY74 (https://dynesty.readthedocs.io/en/stable/index.html) and chromatic (https://zkbt.github.io/chromatic/), each of which use the standard Python libraries scipy98, numpy99, astropy100,101 and matplotlib102. The atmospheric models used to fit the data are as follows: PICASO84 (https://natashabatalha.github.io/picaso/), Virga86 (https://natashabatalha.github.io/virga/), ScCHIMERA68 (https://github.com/mrline/CHIMERA), ATMO88,89,90,91 and PHOENIX94.
References
Lodders, K. & Fegley, B. Atmospheric chemistry in giant planets, brown dwarfs, and low-mass dwarf stars. I. Carbon, nitrogen, and oxygen. Icarus 155, 393â424 (2002).
Zahnle, K., Marley, M. S., Freedman, R. S., Lodders, K. & Fortney, J. J. Atmospheric sulfur photochemistry on hot Jupiters. Astrophys. J. Lett. 701, L20âL24 (2009).
Moses, J. I. et al. Compositional diversity in the atmospheres of hot Neptunes, with application to GJ 436b. Astrophys. J. 777, 34 (2013).
Pollack, J. B. et al. Formation of the giant planets by concurrent accretion of solids and gas. Icarus 124, 62â85 (1996).
Fortney, J. J. et al. A framework for characterizing the atmospheres of low-mass low-density transiting planets. Astrophys. J. 775, 80 (2013).
Venturini, J., Alibert, Y. & Benz, W. Planet formation with envelope enrichment: new insights on planetary diversity. Astron. Astrophys. 596, A90 (2016).
Kaltenegger, L. & Traub, W. A. Transits of Earth-like planets. Astrophys. J. 698, 519â527 (2009).
Deming, D. et al. Discovery and characterization of transiting super Earths using an all-sky transit survey and follow-up by the James Webb Space Telescope. Publ. Astron. Soc. Pac. 121, 952â967 (2009).
Lustig-Yaeger, J., Meadows, V. S. & Lincowski, A. P. The detectability and characterization of the TRAPPIST-1 exoplanet atmospheres with JWST. Astron. J. 158, 27 (2019).
Wakeford, H. R. et al. The complete transmission spectrum of WASP-39b with a precise water constraint. Astron. J. 155, 29 (2018).
Spake, J. J. et al. Abundance measurements of H2O and carbon-bearing species in the atmosphere of WASP-127b confirm its supersolar metallicity. Mon. Not. R. Astron. Soc. 500, 4042â4064 (2021).
Alderson, L. et al. A comprehensive analysis of WASP-17bâs transmission spectrum from space-based observations. Mon. Not. R. Astron. Soc. 512, 4185â4209 (2022).
Stevenson, K. B. et al. Transiting exoplanet studies and community targets for JWSTâs Early Release Science Program. Publ. Astron. Soc. Pac. 128, 094401 (2016).
Bean, J. L. et al. The Transiting Exoplanet Community Early Release Science Program for JWST. Publ. Astron. Soc. Pac. 130, 114402 (2018).
Faedi, F. et al. WASP-39b: a highly inflated Saturn-mass planet orbiting a late G-type star. Astron. Astrophys. 531, A40 (2011).
Laughlin, G., Crismani, M. & Adams, F. C. On the anomalous radii of the transiting extrasolar planets. Astrophys. J. Lett. 729, L7 (2011).
Thorngren, D. P. & Fortney, J. J. Bayesian analysis of hot-Jupiter radius anomalies: evidence for ohmic dissipation? Astron. J. 155, 214 (2018).
Sarkis, P., Mordasini, C., Henning, T., Marleau, G. D. & Mollière, P. Evidence of three mechanisms explaining the radius anomaly of hot Jupiters. Astron. Astrophys. 645, A79 (2021).
Fischer, P. D. et al. HST hot-Jupiter transmission spectral survey: clear skies for cool Saturn WASP-39b. Astrophys. J. 827, 19 (2016).
Nikolov, N. et al. VLT FORS2 comparative transmission spectroscopy: detection of Na in the atmosphere of WASP-39b from the ground. Astrophys. J. 832, 191 (2016).
Kirk, J. et al. LRG-BEASTS: transmission spectroscopy and retrieval analysis of the highly inflated Saturn-mass planet WASP-39b. Astron. J. 158, 144 (2019).
Mancini, L. et al. The GAPS programme with HARPS-N at TNG. XVI. Measurement of the RossiterâMcLaughlin effect of transiting planetary systems HAT-P-3, HAT-P-12, HAT-P-22, WASP-39, and WASP-60. Astron. Astrophys. 613, A41 (2018).
Biazzo, K. et al. The GAPS programme with HARPS-N at TNG. XXXV. Fundamental properties of transiting exoplanet host stars. Astron. Astrophys. 664, A61 (2022).
Polanski, A. S., Crossfield, I. J. M., Howard, A. W., Isaacson, H. & Rice, M. Chemical abundances for 25 JWST exoplanet host stars with KeckSpec. Res. Note. AAS 6, 155 (2022).
Lodders K. & Fegley B. The Planetary Scientistâs Companion (Oxford Univ. Press, 1998).
Kreidberg, L. et al. A precise water abundance measurement for the hot Jupiter WASP-43b. Astrophys. J. Lett. 793, L27 (2014).
Atreya, S. K. et al. The origin and evolution of Saturn: a post-Cassini perspective. Preprint at https://arxiv.org/abs/2205.06914 (2022).
Thorngren, D. & Fortney, J. J. Connecting giant planet atmosphere and interior modeling: constraints on atmospheric metal enrichment. Astrophys. J. Lett. 874, L31 (2019).
Tsiaras, A. et al. A population study of gaseous exoplanets. Astron. J. 155, 156 (2018).
Fisher, C. & Heng, K. Retrieval analysis of 38 WFC3 transmission spectra and resolution of the normalization degeneracy. Mon. Not. R. Astron. Soc. 481, 4698â4727 (2018).
Pinhas, A., Madhusudhan, N., Gandhi, S. & MacDonald, R. H2O abundances and cloud properties in ten hot giant exoplanets. Mon. Not. R. Astron. Soc. 482, 1485â1498 (2019).
Welbanks, L. et al. Massâmetallicity trends in transiting exoplanets from atmospheric abundances of H2O, Na, and K. Astrophys. J. Lett. 887, L20 (2019).
Min, M., Ormel, C. W., Chubb, K., Helling, C. & Kawashima, Y. The ARCiS framework for exoplanet atmospheres. Modelling philosophy and retrieval. Astron. Astrophys. 642, A28 (2020).
Benneke, B. & Seager, S. How to distinguish between cloudy mini-Neptunes and water/volatile-dominated super-Earths. Astrophys. J. 778, 153 (2013).
Jakobsen, P. et al. The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope. I. Overview of the instrument and its capabilities. Astron. Astrophys. 661, A80 (2022).
Birkmann, S. M. et al. The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope. IV. Capabilities and predicted performance for exoplanet characterization. Astron. Astrophys. 661, A83 (2022).
Rigby, J. et al. Characterization of JWST science performance from commissioning. Preprint at https://arxiv.org/abs/2207.05632 (2022).
Rustamkulov, Z., Sing, D. K., Liu, R. & Wang, A. Analysis of a JWST NIRSpec lab time series: characterizing systematics, recovering exoplanet transit spectroscopy, and constraining a noise floor. Astrophys. J. Lett. 928, L7 (2022).
Bell, T. J. et al. Eureka!: an end-to-end pipeline for JWST time-series observations. J. Open Source Sci. 7, 4503â4508 (2022).
Glidic, K. et al. Atmospheric characterization of hot Jupiter CoRoT-1 b using the Wide Field Camera 3 on the Hubble Space Telescope. Astron. J. 164, 19 (2022).
Beichman, C. et al. Observations of transiting exoplanets with the James Webb Space Telescope (JWST). Publ. Astron. Soc. Pac. 126, 1134â1173 (2014).
Yamamura, I., Tsuji, T. & Tanabé, T. AKARI observations of brown dwarfs. I. CO and CO2 bands in the near-infrared spectra. Astrophys. J. 722, 682â698 (2010).
Lodders, K., Palme, H. & Gail, H.-P. in Solar System: Astronomy, Astrophysics, and Cosmology Vol. 4B (ed. Trümper, J. E.) 712 (Landolt-Börnstein - Group VI âAstronomy and Astrophysicsâ, Springer, 2009).
Stevenson, K. B. Quantifying and predicting the presence of clouds in exoplanet atmospheres. Astrophys. J. Lett. 817, L16 (2016).
Fu, G. et al. Statistical analysis of Hubble/WFC3 transit spectroscopy of extrasolar planets. Astrophys. J. Lett. 847, L22 (2017).
Parmentier, V., Fortney, J. J., Showman, A. P., Morley, C. & Marley, M. S. Transitions in the cloud composition of hot Jupiters. Astrophys. J. 828, 22 (2016).
Kataria, T. et al. The atmospheric circulation of a nine-hot-Jupiter sample: probing circulation and chemistry over a wide phase space. Astrophys. J. 821, 9 (2016).
Gao, P. et al. Aerosol composition of hot giant exoplanets dominated by silicates and hydrocarbon hazes. Nat. Astron. 4, 951â956 (2020).
Himes, M. D. & Harrington, J. On the dayside atmosphere of WASP-12b. Astrophys. J. 931, 86 (2022).
Schlawin, E. et al. JWST noise floor. II. Systematic error sources in JWST NIRCam time series. Astron. J. 161, 115 (2021).
Foreman-Mackey, D. et al. exoplanet: gradient-based probabilistic inference for exoplanet data & other astronomical time series. J.Open Source Softw. 6, 3285 (2021).
Salvatier, J., Wiecki, T. V. & Fonnesbeck, C. PyMC3: python probabilistic programming framework. Astrophysics Source Code Library ascl:1610.016 (2016).
Benneke, B. et al. Spitzer observations confirm and rescue the habitable-zone super-Earth K2-18b for future characterization. Astrophys. J. 834, 187 (2017).
Benneke, B. et al. Water vapor and clouds on the habitable-zone sub-Neptune exoplanet K2-18b. Astrophys. J. Lett. 887, L14 (2019).
Benneke, B. et al. A sub-Neptune exoplanet with a low-metallicity methane-depleted atmosphere and Mie-scattering clouds. Nat. Astron. 3, 813â821 (2019).
Kreidberg, L. batman: basic transit model calculation in Python. Publ. Astron. Soc. Pac. 127, 1161â1165 (2015).
Kirk, J. et al. LRG-BEASTS III: ground-based transmission spectrum of the gas giant orbiting the cool dwarf WASP-80. Mon. Not. R. Astron. Soc. 474, 876â885 (2018).
Kirk, J. et al. ACCESS and LRG-BEASTS: a precise new optical transmission spectrum of the ultrahot Jupiter WASP-103b. Astron. J. 162, 34 (2021).
Laginja, I. & Wakeford, H. ExoTiC-ISM: a Python package for marginalised exoplanet transit parameters across a grid of systematic instrument models. J Open Source Softw 5, 2281 (2020).
Wakeford, H. R. & Grant, D. Exo-TiC/ExoTiC-LD: ExoTiC-LD v2.1 Zenodo Release. Zenodo https://doi.org/10.5281/zenodo.6809899 (2022).
Magic, Z., Chiavassa, A., Collet, R. & Asplund, M. The Stagger-grid: a grid of 3D stellar atmosphere models. IV Limb darkening coefficients. Astron. Astrophys. 573, A90 (2015).
Sing, D. K. Stellar limb-darkening coefficients for CoRot and Kepler. Astron. Astrophys. 510, A21 (2010).
Foreman-Mackey, D. et al. emcee v3: a Python ensemble sampling toolkit for affine-invariant MCMC. J. Open Source Softw. 4, 1864 (2019).
Marley, M. S. & Robinson, T. D. On the cool side: modeling the atmospheres of brown dwarfs and giant planets. Ann. Rev. Astron. Astrophys. 53, 279â323 (2015).
Hubeny, I. Model atmospheres of sub-stellar mass objects. Mon. Not. R. Astron. Soc. 469, 841â869 (2017).
Piskorz, D. et al. Ground- and space-based detection of the thermal emission spectrum of the transiting hot Jupiter KELT-2Ab. Astron. J. 156, 133 (2018).
Arcangeli, J. et al. Hâ opacity and water dissociation in the dayside atmosphere of the very hot gas giant WASP-18b. Astrophys. J. Lett. 855, L30 (2018).
Mansfield, M. et al. A unique hot Jupiter spectral sequence with evidence for compositional diversity. Nat. Astron. 5, 1224â1232 (2021).
Fortney, J. J., Marley, M. S., Lodders, K., Saumon, D. & Freedman, R. Comparative planetary atmospheres: models of TrES-1 and HD 209458b. Astrophys. J. Lett. 627, L69âL72 (2005).
Perez-Becker, D. & Showman, A. P. Atmospheric heat redistribution on hot Jupiters. Astrophys. J. 776, 134 (2013).
Zhang, X. Atmospheric regimes and trends on exoplanets and brown dwarfs. Res. Astron. Astrophys. 20, 099 (2020).
Mai, C. & Line, M. R. Exploring exoplanet cloud assumptions in JWST transmission spectra. Astrophys. J. 883, 144 (2019).
Iyer, A. R. & Line, M. R. The influence of stellar contamination on the interpretation of near-infrared transmission spectra of sub-Neptune worlds around M-dwarfs. Astrophys. J. 889, 78 (2020).
Speagle, J. S. DYNESTY: a dynamic nested sampling package for estimating Bayesian posteriors and evidences. Mon. Not. R. Astron. Soc. 493, 3132â3158 (2020).
Polyansky, O. L. et al. ExoMol molecular line lists XXX: a complete high-accuracy line list for water. Mon. Not. R. Astron. Soc. 480, 2597â2608 (2018).
Gharib-Nezhad, E. et al. EXOPLINES: molecular absorption cross-section database for brown dwarf and giant exoplanet atmospheres. Astrophys. J. Suppl. Ser. 254, 34 (2021).
Li, G. et al. Rovibrational line lists for nine isotopologues of the CO molecule in the X1Σ+ ground electronic state. Astrophys. J. Suppl. Ser. 216, 15 (2015).
Grimm, S. L. et al. HELIOS-K 2.0 opacity calculator and open-source opacity database for exoplanetary atmospheres. Astrophys. J. Suppl. Ser. 253, 30 (2021).
Huang, X., Gamache, R. R., Freedman, R. S., Schwenke, D. W. & Lee, T. J. Reliable infrared line lists for 13 CO2 isotopologues up to Eâ²=18,000 cmâ1 and 1500 K, with line shape parameters. J. Quant. Spectrosc. Radiat. Transf. 147, 134â144 (2014).
Marley, M. S. et al. The Sonora brown dwarf atmosphere and evolution models. I. Model description and application to cloudless atmospheres in rainout chemical equilibrium. Astrophys. J. 920, 85 (2021).
Azzam, A. A. A., Tennyson, J., Yurchenko, S. N. & Naumenko, O. V. ExoMol molecular line listsâXVI. The rotationâvibration spectrum of hot H2S. Mon. Not. R. Astron. Soc. 460, 4063â4074 (2016).
Hargreaves, R. J. et al. An accurate, extensive, and practical line list of methane for the HITEMP database. Astrophys. J. Suppl. Ser. 247, 55 (2020).
Marley, M. S. & McKay, C. P. Thermal structure of Uranusâ atmosphere. Icarus 138, 268â286 (1999).
Batalha, N. E., Marley, M. S., Lewis, N. K. & Fortney, J. J. Exoplanet reflected-light spectroscopy with PICASO. Astrophys. J. 878, 70 (2019).
Mukherjee, S., Batalha, N. E., Fortney, J. J. & Marley, M. S. PICASO 3.0: a one-dimensional climate model for giant planets and brown dwarfs. Preprint at https://arXiv.org/abs/2208.07836 (2022).
Rooney, C. M., Batalha, N. E., Gao, P. & Marley, M. S. A new sedimentation model for greater cloud diversity in giant exoplanets and brown dwarfs. Astrophys. J. 925, 33 (2022).
Ackerman, A. S. & Marley, M. S. Precipitating condensation clouds in substellar atmospheres. Astrophys. J. 556, 872 (2001).
Tremblin, P. et al. Fingering convection and cloudless models for cool brown dwarf atmospheres. Astrophys. J. 804, L17 (2015).
Drummond, B. et al. The effects of consistent chemical kinetics calculations on the pressureâtemperature profiles and emission spectra of hot Jupiters. Astron. Astrophys. 594, A69 (2016).
Goyal, J. M. et al. A library of ATMO forward model transmission spectra for hot Jupiter exoplanets. Mon. Not. R. Astron. Soc. 474, 5158â5185 (2018).
Goyal, J. M. et al. A library of self-consistent simulated exoplanet atmospheres. Mon. Not. R. Astron. Soc. 498, 4680â4704 (2020).
Hauschildt, P. H., Allard, F. & Baron, E. The NextGen Model atmosphere grid for 3000 ⤠Teff ⤠10,000 K. Astrophys. J. 512, 377 (1999).
Barman, T. S., Hauschildt, P. H. & Allard, F. Irradiated planets. Astrophys. J. 556, 885 (2001).
Lothringer, J. D. & Barman, T. The PHOENIX exoplanet retrieval algorithm and using Hâ opacity as a probe in ultrahot Jupiters. Astron. J. 159, 289 (2020).
Lothringer, J. D. et al. New window into planet formation and migration: refractory-to-volatile elemental ratios in ultra-hot Jupiters. Astrophys. J. 914, 12 (2021).
Gregory, P. C. (ed.) Bayesian Logical Data Analysis for the Physical Sciences: A Comparative Approach with âMathematicaâ Support (Cambridge Univ. Pres, 2005).
Trotta, R. Bayes in the sky: Bayesian inference and model selection in cosmology. Contemp. Phys. 49, 71â104 (2008).
Virtanen, P. et al. SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat. Methods 17, 261â272 (2020).
Harris, C. R. et al. Array programming with NumPy. Nature 585, 357â362 (2020).
The Astropy Collaboration Astropy: a community Python package for astronomy. Astron. Astrophys. 558, A33 (2013).
The Astropy Collaboration The Astropy Project: building an open-science project and status of the v2.0 core package. Astron. J. 156, 123 (2018).
Hunter, J. D. Matplotlib: a 2D graphics environment. Comput. Sci. Eng. 9, 90â95 (2007).
Acknowledgements
This work is based on observations made with the NASA/ESA/CSA James Webb Space Telescope. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with programme JWST-ERS-01366. Support for programme JWST-ERS-01366 was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127. The results reported herein benefited during the design phase from collaborations and/or information exchange within NASAâs Nexus for Exoplanet System Science (NExSS) research coordination network sponsored by NASAâs Science Mission Directorate. N.M.B. acknowledges funding from NASAâs Interdisciplinary Consortia for Astrobiology Research (NNH19ZDA001N-ICAR) under award number 19-ICAR19_2-0041.
Author information
Authors and Affiliations
Consortia
Contributions
N.M.B., J.L.B. and K.B.S. provided overall programme leadership and management. D.K.S., E.M.-R.K., H.R.W., I.J.M.C., J.L.B., K.B.S., L.K., M.L.-M., M.R.L., N.M.B., V.P. and Z.K.B.-T. made significant contributions to the design of the programme. K.B.S. generated the observing plan with input from the team. E.S., N.E. and T.G.B. provided instrument expertise. B.B., E.M.-R.K., H.R.W., I.J.M.C., J.L.B., L.K., M.L.-M., M.R.L., N.M.B. and Z.K.B.-T. led or co-led working groups and/or contributed to significant strategic planning efforts like the design and implementation of the pre-launch Data Challenges. A.L.C., D.K.S., E.S., N.E., N.P.G., T.G.B. and V.P. generated simulated data for pre-launch testing of methods. A.L.C., A.D.F., C.P., E.-M.A., E.S., J.K., L.A., T.J.B. and Z.R. contributed to the development of data analysis pipelines and/or provided the data analysis products used in this analysis, reduced the data, modelled the light curves, and produced the planetary spectrum. J.M.G., J.D.L., K.O., M.R.L., N.E.B., S.M. and S.E.M. generated theoretical model grids for comparison with data. A.L.C., D.K.S., E.S., H.R.W., J.K., J.J.F., J.M.G., J.L.B., J.D.L., K.B.S., K.O., M.R.L., N.E.B., N.M.B., S.M., S.E.M., Z.K.B.-T. and Z.R. contributed significantly to the writing of this manuscript. E.-M.A., E.S., M.R.L. and Z.K.B.-T. generated figures for this manuscript. Over 100 of the co-authors participated in the design of this programme and contributed to the ERS observing proposal. Nearly all co-authors participated in pre-launch working groups and programme events. Nearly all authors provided scientific and technical input on the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature thanks Masahiro Ikoma and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Comparison of transmission spectrum modelling results from different codes for WASP-39b.
Despite different radiativeâconvective equilibrium and chemical solvers, treatments of clouds, grid spacing and grid-fitting approaches, all four grids arrive at the same 10-times solar metallicity point solution. Additionally, all four provide an acceptable fit to the data, with best-fitting Χ2/Ndata < 1.4.
Extended Data Fig. 2 Atmospheric structure arising from the best-fit model.
The thick red curve (and corresponding top x axis) shows the resulting 1D radiativeâconvective equilibrium temperature profile. The dashed lines (and bottom x axis) show the vertical gas mixing ratio profiles under the assumption of thermochemical equilibrium. These abundances, along with the absorption cross-sections shown in the bottom panel of Fig. 3, are what control the relative contributions of each gaseous opacity to the total transmission spectrum.
Extended Data Fig. 3 Assessment of the strength of spectral features for WASP-39b.
Residual features (blue data points) after subtracting the continuum best model (black âno CO2â model curve in Fig. 3). A best-fitting ensemble of a two-component Gaussian model to both the CO2 feature and the unknown absorber feature (~4âµm) is shown in red.
Supplementary information
Supplementary Data
This file contains files to recreate the main text figures. It includes: TRANSMISSION_SPECTRA_DATA, containing the data from each reduction; MODEL_FITS, containing the best-fitting model spectra from each of the four model grids described in the manuscript as well as the âremove one gas at a timeâ spectra to recreate the top panel of Fig. 3; and the data necessary to reproduce the light curves shown in Fig. 1 are contained in the CSV file.
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
JWST Transiting Exoplanet Community Early Release Science Team. Identification of carbon dioxide in an exoplanet atmosphere. Nature 614, 649â652 (2023). https://doi.org/10.1038/s41586-022-05269-w
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-022-05269-w
This article is cited by
-
Inferring chemical disequilibrium biosignatures for Proterozoic Earth-like exoplanets
Nature Astronomy (2024)
-
A warm Neptuneâs methane reveals core mass and vigorous atmospheric mixing
Nature (2024)
-
Sulfur dioxide in the mid-infrared transmission spectrum of WASP-39b
Nature (2024)
-
A benchmark JWST near-infrared spectrum for the exoplanet WASP-39 b
Nature Astronomy (2024)
-
Atmospheric carbon depletion as a tracer of water oceans and biomass on temperate terrestrial exoplanets
Nature Astronomy (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.