Sulfur dioxide in the mid-infrared transmission spectrum of WASP-39b.

The recent inference of sulfur dioxide (SO2) in the atmosphere of the hot (approximately 1,100 K), Saturn-mass exoplanet WASP-39b from near-infrared JWST observations1-3 suggests that photochemistry is a key process in high-temperature exoplanet atmospheres4. This is because of the low (<1 ppb) abundance of SO2 under thermochemical equilibrium compared with that produced from the photochemistry of H2O and H2S (1-10 ppm)4-9. However, the SO2 inference was made from a single, small molecular feature in the transmission spectrum of WASP-39b at 4.05 µm and, therefore, the detection of other SO2 absorption bands at different wavelengths is needed to better constrain the SO2 abundance. Here we report the detection of SO2 spectral features at 7.7 and 8.5 µm in the 5-12-µm transmission spectrum of WASP-39b measured by the JWST Mid-Infrared Instrument (MIRI) Low Resolution Spectrometer (LRS)10. Our observations suggest an abundance of SO2 of 0.5-25 ppm (1σ range), consistent with previous findings4. As well as SO2, we find broad water-vapour absorption features, as well as an unexplained decrease in the transit depth at wavelengths longer than 10 µm. Fitting the spectrum with a grid of atmospheric forward models, we derive an atmospheric heavy-element content (metallicity) for WASP-39b of approximately 7.1-8.0 times solar and demonstrate that photochemistry shapes the spectra of WASP-39b across a broad wavelength range.

Photochemically produced SO2 in the atmosphere of WASP-39b.

Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability1. However, no unambiguous photochemical products have been detected in exoplanet atmospheres so far. Recent observations from the JWST Transiting Exoplanet Community Early Release Science Program2,3 found a spectral absorption feature at 4.05 µm arising from sulfur dioxide (SO2) in the atmosphere of WASP-39b. WASP-39b is a 1.27-Jupiter-radii, Saturn-mass (0.28 MJ) gas giant exoplanet orbiting a Sun-like star with an equilibrium temperature of around 1,100 K (ref. 4). The most plausible way of generating SO2 in such an atmosphere is through photochemical processes5,6. Here we show that the SO2 distribution computed by a suite of photochemical models robustly explains the 4.05-µm spectral feature identified by JWST transmission observations7 with NIRSpec PRISM (2.7σ)8 and G395H (4.5σ)9. SO2 is produced by successive oxidation of sulfur radicals freed when hydrogen sulfide (H2S) is destroyed. The sensitivity of the SO2 feature to the enrichment of the atmosphere by heavy elements (metallicity) suggests that it can be used as a tracer of atmospheric properties, with WASP-39b exhibiting an inferred metallicity of about 10× solar. We further point out that SO2 also shows observable features at ultraviolet and thermal infrared wavelengths not available from the existing observations.

High atmospheric metal enrichment for a Saturn-mass planet.

Atmospheric metal enrichment (that is, elements heavier than helium, also called 'metallicity') is a key diagnostic of the formation of giant planets1-3. The giant planets of the Solar System show an inverse relationship between mass and both their bulk metallicities and atmospheric metallicities. Extrasolar giant planets also display an inverse relationship between mass and bulk metallicity4. However, there is significant scatter in the relationship and it is not known how atmospheric metallicity correlates with either planet mass or bulk metallicity. Here we show that the Saturn-mass exoplanet HD 149026b (refs. 5-9) has an atmospheric metallicity 59-276 times solar (at 1σ), which is greater than Saturn's atmospheric metallicity of roughly 7.5 times solar10 at more than 4σ confidence. This result is based on modelling CO2 and H2O absorption features in the thermal emission spectrum of the planet measured by the James Webb Space Telescope. HD 149026b is the most metal-rich giant planet known, with an estimated bulk heavy element abundance of 66 ± 2% by mass11,12. We find that the atmospheric metallicities of both HD 149026b and the Solar System giant planets are more correlated with bulk metallicity than planet mass.

Early Release Science of the exoplanet WASP-39b with JWST NIRCam.

Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy (for example, refs. 1,2) provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength coverage, moderate spectral resolution and high precision, which, together, are not achievable with previous observatories. Now that JWST has commenced science operations, we are able to observe exoplanets at previously uncharted wavelengths and spectral resolutions. Here we report time-series observations of the transiting exoplanet WASP-39b using JWST's Near InfraRed Camera (NIRCam). The long-wavelength spectroscopic and short-wavelength photometric light curves span 2.0-4.0 micrometres, exhibit minimal systematics and reveal well defined molecular absorption features in the planet's spectrum. Specifically, we detect gaseous water in the atmosphere and place an upper limit on the abundance of methane. The otherwise prominent carbon dioxide feature at 2.8 micrometres is largely masked by water. The best-fit chemical equilibrium models favour an atmospheric metallicity of 1-100-times solar (that is, an enrichment of elements heavier than helium relative to the Sun) and a substellar C/O ratio. The inferred high metallicity and low C/O ratio may indicate significant accretion of solid materials during planet formation (for example, refs. 3,4,) or disequilibrium processes in the upper atmosphere (for example, refs. 5,6).

Vertically Resolved Magma Ocean-Protoatmosphere Evolution: H2, H2O, CO2, CH4, CO, O2, and N2 as Primary Absorbers.

The earliest atmospheres of rocky planets originate from extensive volatile release during magma ocean epochs that occur during assembly of the planet. These establish the initial distribution of the major volatile elements between different chemical reservoirs that subsequently evolve via geological cycles. Current theoretical techniques are limited in exploring the anticipated range of compositional and thermal scenarios of early planetary evolution, even though these are of prime importance to aid astronomical inferences on the environmental context and geological history of extrasolar planets. Here, we present a coupled numerical framework that links an evolutionary, vertically resolved model of the planetary silicate mantle with a radiative-convective model of the atmosphere. Using this method, we investigate the early evolution of idealized Earth-sized rocky planets with end-member, clear-sky atmospheres dominated by either H2, H2O, CO2, CH4, CO, O2, or N2. We find central metrics of early planetary evolution, such as energy gradient, sequence of mantle solidification, surface pressure, or vertical stratification of the atmosphere, to be intimately controlled by the dominant volatile and outgassing history of the planet. Thermal sequences fall into three general classes with increasing cooling timescale: CO, N2, and O2 with minimal effect, H2O, CO2, and CH4 with intermediate influence, and H2 with several orders of magnitude increase in solidification time and atmosphere vertical stratification. Our numerical experiments exemplify the capabilities of the presented modeling framework and link the interior and atmospheric evolution of rocky exoplanets with multiwavelength astronomical observations.