Limb Asymmetries on WASP-39b: A Multi-GCM Comparison of Chemistry, Clouds, and Hazes

arXiv:2509.21588 · astro-ph.EP · Submitted 2025-09-25 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Today's paper: "Limb Asymmetries on WASP-39b".

Jocelyn: Given the potential for millions in cost associated with error, this synthesis must be precise, structured, and exhaustive.

Vera: First, who's behind it and why it matters.

Paper summary: Vera: So this paper, "Limb Asymmetries on WASP-39b: A Multi-GCM Comparison of Chemistry, Clouds, and Hazes," it looks at the observations from JWST on WASP-39b where they finally observed separate spectra for the morning and evening limbs <ref:2509.21588#pg1,Limb Asymmetries on WASP-39b: A Multi-GCM Comparison of Chemistry, Clouds>.

Jocelyn: The main thing this paper is doing is taking those observations and comparing them against five different General Circulation Models, or GCMs, to see how well those models can explain what we’re seeing.

Vera: The authors are looking at specific asymmetries they found: a transit radius difference of about four hundred ppm and a stronger four point three mu m CO two feature in the evening terminator compared to the morning one.

Jocelyn: It’s important because they want to figure out what causes those asymmetries, whether it’s temperature differences, cloud effects, or something else entirely.

Subrahmanyan: From a theoretical standpoint, this is about testing how atmospheric physics—the chemistry and the dynamics—actually works on these hot Jupiters. They're using these models as tools to see which physical processes match the observational data we get from JWST (<ref:2509.21588#pg1>).

Vera: They are using different model setups, including one that simulates disequilibrium thermochemistry and another that includes cloud radiative feedback, to check the results.

Jocelyn: And they actually post-process the temperature structures from all those simulations using a 2D photochemical model and one simulation with a cloud microphysics model <ref:2509.21588#pg1>.

Subrahmanyan: What’s striking is that even though the predicted temperatures in these different models vary quite a bit, they seem remarkably consistent when it comes to predicting those morning–evening temperature differences (<ref:2509.21588#pg2>).

Vera: That consistency in the temperature structure is an interesting finding, especially since we know that’s one factor causing these asymmetries.

Jocelyn: The paper also looks at how chemistry plays a role, specifically seeing what equilibrium-chemistry simulations predict versus what we actually observe.

Conclusion: Vera: So looking at the whole picture from "Limb Asymmetries on WASP-39b: A Multi-GCM Comparison of Chemistry, Clouds, and Hazes," it seems like the authors are really stressing how complex this asymmetry is <ref:2509.21588#pg1,Limb Asymmetries on WASP-39b: A Multi-GCM Comparison of Chemistry, Clouds>.

Jocelyn: The core puzzle they’re highlighting is that while their temperature models match up pretty well, they still can’t get the big difference in the four point three mu m CO two feature down to about two thousand ppm <ref:2509.21588#pg1>.

Subrahmanyan: That large CO two asymmetry is where things get tricky; it suggests that just having a temperature gradient isn't enough to explain the observation <ref:2509.21588#pg1>.

Vera: And they point out that including photochemistry and clouds didn't help reproduce that two thousand ppm difference in the CO two feature <ref:2509.21588#pg1>.

Jocelyn: It seems like they suggest a combination of things, like varying metallicity or unexplored parameters in cloud models, might be the reason for this mismatch.

Subrahmanyan: So what this means for us is that we probably need future models that integrate cloud microphysics and feedback across a much wider range of possibilities to really nail this observation.

Vera: Exactly, it tells us that while we have good temperature structure predictions, the cloudy parts of the atmosphere are where the real complexity lies.

Jocelyn: So next time you look at WASP-39b spectra, remember that those small asymmetries might be pointing toward a need for more detailed cloud physics in our simulations <ref:2509.21588#pg1>.

Department of Astronomy & Astrophysics, University of Chicago, Chicago, IL 60637, USA · Max Planck Institute for Astronomy, Heidelberg, Germany · Department of Astronomy, University of Maryland, College Park, MD 20742, USA · Space Research Institute, Austrian Academy of Sciences · Department of Earth and Planetary Sciences, University of California Riverside CA USA · Institute of Astronomy & Astrophysics Academia Sinica Taipei Taiwan · Center for Space and Habitability University of Bern Switzerland · Universit¨ats-Sternwarte M¨unchen Fakultät f¨ur Physik der Ludwig-Maximilians-Universit¨at München Deutschland · Department of Astronomy and Astrophysics, University of Michigan Ann Arbor MI USA · Institute of Astronomy, KU Leuven Belgium · Institute for Theoretical Physics and Computational Physics Graz University of Technology Austria · University of Texas at Austin Department of Astronomy Austin TX USA · Centre for Science at Extreme Conditions and Scottish Universities Physics Alliance School of Physics and Astronomy The University of Edinburgh UK · Department of Astronomy, University Wisconsin–Madison WI USA · Anton Pannekoek Institute for Astronomy University of Amsterdam The Netherlands · Facultad de Ingenier´ıa y Ciencias Universidad Adolfo Ib´a˜nez Santiago Chile · School of Physics and Astronomy University of Leicester United Kingdom · Space Telescope Science Institute Baltimore MD USA · William H. Miller III Department of Physics and Astronomy Johns Hopkins University Baltimore MD USA · University College London Department of Physics & Astronomy Gower St London WC1E 6BT United Kingdom · Astronomy & Astrophysics Group Department of Physics University of Warwick Coventry CV4 7AL United Kingdom · Institute of Planetary Research German Aerospace Center DLR Berlin Germany · Department of Physics Imperial College London Prince Consort Road SW7 2AZ London UK · Steward Observatory Tucson AZ USA · Department of Earth, Atmospheric and Planetary Sciences Massachusetts Institute of Technology Cambridge MA USA · Kavli Institute for Astrophysics and Space Research Massachusetts Institute of Technology Cambridge MA USA

astro-ph.EP

Submitted: 2025-09-25

Updated: 2026-10-08

Comments: revised version after peer review, accepted for publication in ApJ. Supplementary Data available on Zenodo

Code: https://github.com/exoclime/VULCAN

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 89/100

The gist: Given the potential for millions in cost associated with error, this synthesis must be precise, structured, and exhaustive.

Key concepts

Limb Asymmetries: These are differences in brightness between the morning and evening sides of a planet's atmosphere when viewed from space. For WASP-39b, this manifested as a larger transit radius and a stronger $\text{CO}_2$ feature on the evening terminator compared to the morning one.
Disequilibrium Thermochemistry: This process involves chemical reactions occurring faster than they can reach equilibrium, often due to rapid physical changes in the atmosphere. In this study, including disequilibrium processes helped explain why methane features disappeared in the morning spectrum when observed.
Cloud Microphysics: This refers to the detailed physical processes governing how clouds form and behave, such as particle size and composition. The paper suggests that current models may be underestimating terminator differences because they lack a broader exploration of these microphysics assumptions.
$4.3\ \mu\text{m}$ $\text{CO}_2$ Feature Discrepancy: This is the central puzzle: the observed difference in the strength of the infrared $\text{CO}_2$ absorption line between morning and evening terminator spectra. No model, even with photochemistry and clouds, could accurately predict this large asymmetry, indicating a missing physical driver.

Terminology

Metallicity — 1× in this paper · explained in 3 episodes

Summary

Given the potential for millions in cost associated with error, this synthesis must be precise, structured, and exhaustive.

Here is the detailed summary:

Comprehensive Research Synthesis: Limb Asymmetries on WASP-39b – A Multi-GCM Comparison of Chemistry, Clouds, and Hazes

This research investigates the observed limb asymmetries in the spectra of hot Jupiters using James Webb Space Telescope (JWST) data. The study focuses specifically on WASP-39b, where observations have revealed significant differences between the morning and evening terminators—notably a larger transit radius by approximately 400 ppm and a substantially stronger 4.3 mu m CO 2 feature in the evening terminator compared to the morning one. The core objective of the work is to interpret these asymmetries by comparing high-resolution observations with forward models derived from five different General Circulation Models (GCMs), incorporating various physical processes such as disequilibrium thermochemistry, cloud radiative feedback, and photochemical effects.

Model Consistency and Temperature Structure

A key finding is the surprising consistency among the GCM simulations regarding temperature differences between the morning and evening terminators. Despite variations in predicted temperatures across different models—which can span several hundred Kelvin within the observable atmosphere due to model differences—the clear-atmosphere models exhibit remarkable agreement in predicting these terminator temperature contrasts. Specifically, temperature differences typically peak near 0.01 bar with values ranging from 200 to 300 K for models employing full wavelength-dependent radiative transfer, and up to 400 K in models utilizing picket-fence radiative transfer. These temperature differences tend to slightly decline at lower pressures.

Chemical Features: Equilibrium vs. Disequilibrium Chemistry

The study critically examined the role of atmospheric chemistry:

Equilibrium Chemistry Predictions: Simulations assuming equilibrium chemistry often predict prominent methane (CH 4) absorption features in the morning spectrum, which are not present in the actual observations.
Disequilibrium Effects: Including disequilibrium processes significantly alters these predictions. Horizontal transport is shown to homogenize methane across the atmosphere, leading to its disappearance from the morning spectrum and thus smaller morning-evening differences at wavelengths shorter than 4 mu m.

The CO 2 Feature Discrepancy: The Central Puzzle

The most significant challenge identified by the researchers is the failure of current models to reproduce a crucial observational detail: the large amplitude difference in the 4.3 mu m CO 2 feature (about 2000 ppm) observed between the morning and evening terminators.

Model Failure: None of the tested models—even those incorporating photochemistry and clouds—were able to replicate this large CO 2 asymmetry.
Photochemistry and Transport Impact: Photochemical destruction of CH 4 at low pressures (< 1 mbar) and its subsequent homogenization via atmospheric transport at intermediate pressures effectively eliminate the predicted methane features, aligning with observations in the NIR. Furthermore, photochemical and transport processes have a moderate impact on CO 2 abundance; the amplitude of the CO 2 feature barely changes when these processes are included.
Cloud Feedback Limitations: While radiative feedback from condensate clouds can cool both terminators by several hundred Kelvin, this effect was insufficient to account for the observed CO 2 amplitude difference. Moreover, cloud microphysics models predict much larger terminator differences based on particle size, mass fraction, and composition than what is captured by the simplified equilibrium cloud models used in the study.

The Role of Clouds: A Second-Order Effect?

The research strongly suggests that clouds are not the primary driver of the observed CO 2 asymmetry. Instead, they appear to play a second-order role compared to the dominant horizontal temperature gradient. The analysis indicates that the observed evening spectrum, particularly the CO 2 feature amplitude, is more consistent with a cloud-free atmosphere.

However, a potential mechanism for explaining the CO 2 asymmetry has been proposed: Cloud formation restricted specifically to the morning terminator. This scenario could be achieved either by having an exceptionally hot evening terminator or through a broader exploration of assumptions within microphysics models.

Conclusion and Future Directions

The overarching conclusion of this work is that cloudy models, in particular, may be underestimating the true differences between the morning and evening terminators. The authors emphasize that a larger exploration of assumptions in microphysics models, coupled with studies that combine radiative feedback and cloud microphysics across a broader parameter space, is necessary to resolve the persistent discrepancy.

Furthermore, they suggest leveraging analyses from other targets (e.g., WASP107b, WASP-94Ab) and employing a variety of techniques and codes beyond the current framework (such as Catwoman) to increase confidence in low-resolution limb asymmetry measurements. The paper remains cognizant that many avenues are unexplored, necessitating significant future work in this area.

Source References Cited:

Steinrueck, M. E., Christie, D., Savel, A., et al. (2025): Model Outputs and Spectra for ”Limb Asymmetries on WASP-39b: A Multi-GCM Comparison of Chemistry, Clouds, and Hazes” (Zenodo).
Mbarek, R., & Kempton, E. M. R. (2016): ApJ, 827, 121
Savel, A. B., Kempton, E. M. R., Malik, M., et al. (2022): ApJ, 926, 85
Savel, A. B., Kempton, E. M. R., Rauscher, E., et al. (2023): ApJ, 944, 99
Steinrueck, M. E., Koskinen, T., Lavvas, P., et al. (2023): ApJ, 951, 117
Steinrueck, M. E., Showman, A. P., Lavvas, P., et al. (2021): MNRAS, 504, 2783
Wardenier, J. P., Parmentier, V., Lee, E. K. H., Line, M. R., & Gharib-Nezhad, E. (2021): MNRAS, 506, 1258
Wardenier, J. P., Parmentier, V., Line, M. R., & Lee, E. K. H. (2023): MNRAS, 525, 4942
Tsai, S.-M., Moses, J. I., Powell, D., & Lee, E. K. H. (2023): ApJL, 959, L30
Tsai, S.-M., Lee, E. K. H., Powell, D., et al. (2023a): Nature, 617, 483
Tsai, S.-M., Parmentier, V., Mendon¸ca, J. M., et al. (2024): ApJ, 963, 41
Stock, J. W., Kitzmann, D., & Patzer, A. B. C. (2022): Monthly Notices of the Royal Astronomical Society, 517, 4070
Stock, J. W., Kitzmann, D., Patzer, A. B. C., & Sedlmayr, E. (2018): Monthly Notices of the Royal Astronomical Society, 479, 865
Malsky, I., Rauscher, E., Roman, M. T., et al. (2024): ApJ, 961, 66
Roman, M. T., Kempton, E. M. R., Rauscher, E., et al.

Improvements for AI systems

Improved AI system can perform multi-GCM comparison for limb asymmetry, allowing it to quantify the extent to which GCMs of hot Jupiters in their typical configurations vary in their predictions regarding limb asymmetries as stated in Section 1.3.
Improved AI system can integrate forward modeling with photochemical and cloud microphysics models, enabling it to compare simulated outcomes against observations by assessing factors such as temperature differences, photo/thermochemistry, clouds and hazes, which is the core challenge addressed in the abstract.
Improved AI system can predict spectral features by simulating various physical regimes (e.g., clear-atmosphere models, cloud-free, equilibrium-chemistry GCM synthetic spectra) and determine which processes drive observed asymmetries, such as distinguishing between temperature differences and effects of photochemistry or cloud radiative feedback.
Improved AI system can perform sensitivity analysis by systematically exploring model parameter variations, specifically suggesting that increasing metallicity is expected to increase horizontal temperature contrasts and that changes to the assumed metallicity (and/or C/O ratio) could possibly bring the strength of the predicted CO2 feature into better agreement with the observations.
Improved AI system can assess observational limitations by determining if observed discrepancies are due to measurement uncertainty, noting that the CO2 feature is only resolved with a few data points, all of which have relatively large error bars, as discussed in Section 7.3.1.

Abstract

With JWST, observing separate spectra of the morning and evening limbs of hot Jupiters has finally become a reality. The first such observation was reported for WASP-39b, where the evening terminator was observed to have a larger transit radius by about 400 ppm and a stronger 4.3 μ m CO 2 feature than the morning terminator. Multiple factors, including temperature differences, photo/thermochemistry, clouds and hazes, could cause such limb asymmetries. To interpret these new limb asymmetry observations, a detailed understanding of how the relevant processes affect morning and evening spectra grounded in forward models is needed. Focusing on WASP-39b, we compare simulations from five different general circulation models (GCMs), including one simulating disequilibrium thermochemistry and one with cloud radiative feedback, to the recent WASP-39b limb asymmetry observations. We also post-process the temperature structures of all simulations with a 2D photochemical model and one simulation with a cloud microphysics model. Although the temperatures predicted by the different models vary considerably, the models are remarkably consistent in their predicted morning--evening temperature differences. Several equilibrium-chemistry simulations predict strong methane features in the morning spectrum, not seen in the observations. When including disequilibrium processes, horizontal transport homogenizes methane, and these methane features disappear. However, even after including photochemistry and clouds, our models still cannot reproduce the observed about 2000 ppm asymmetry in the CO 2 feature. A combination of factors, such as varying metallicity and unexplored parameters in cloud models, may explain the discrepancy, emphasizing the need for future models integrating cloud microphysics and feedback across a broader parameter space.