Water in HD 209458bs atmosphere from 3.6 - 8 microns IRAC photometric observations in primary transit

The hot Jupiter HD 209458b was observed during primary transit at 3.6, 4.5, 5.8 and 8.0 microns using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We detail here the procedures we adopted to correct for the systematic trends present in the IRAC data. The light curves were fitted including limb darkening effects and fitted using Markov Chain Monte Carlo and prayer-bead Monte Carlo techniques, finding almost identical results. The final depth measurements obtained by a combined Markov Chain Monte Carlo fit are at 3.6 microns, 1.469 +- 0.013 % and 1.448 +- 0.013 %; at 4.5 microns, 1.478 +- 0.017 % ; at 5.8 microns, 1.549 +- 0.015 % and at 8.0 microns 1.535 +- 0.011 %. Our results clearly indicate the presence of water in the planetary atmosphere. Our broad band photometric measurements with IRAC prevent us from determining the additional presence of other other molecules such as CO, CO2 and methane for which spectroscopy is needed. While water vapour with a mixing ratio of 10^-4-10^-3 combined with thermal profiles retrieved from the day-side may provide a very good fit to our observations, this data set alone is unable to resolve completely the degeneracy between water abundance and atmospheric thermal profile.

💡 Research Summary

The paper presents a detailed analysis of primary‑transit observations of the hot‑Jupiter HD 209458b obtained with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope at four wavelengths: 3.6, 4.5, 5.8 and 8.0 µm. The authors first describe the raw data reduction and the extensive systematic corrections required for each IRAC channel. These corrections address intra‑pixel sensitivity variations, pixel‑phase effects, temporal drifts, and background fluctuations. A position‑dependent correction function is applied to the 3.6 µm channel, while time‑dependent baseline models are used for the longer‑wavelength channels.

After systematic removal, the light curves are fitted with a transit model that incorporates stellar limb darkening. To extract the transit depth and other orbital parameters, two independent statistical techniques are employed: a Markov Chain Monte Carlo (MCMC) analysis with broad priors, and a residual‑permutation (prayer‑bead) Monte Carlo approach that reshuffles the observed residuals to test the robustness of the fit against correlated noise. Both methods converge on virtually identical depth values, confirming the reliability of the results.

The measured transit depths are 1.469 ± 0.013 % at 3.6 µm, 1.478 ± 0.017 % at 4.5 µm, 1.549 ± 0.015 % at 5.8 µm, and 1.535 ± 0.011 % at 8.0 µm. These values are consistent with the presence of a strong water‑vapor absorption band across the IRAC bandpasses, especially the slight increase in depth at 5.8 µm and 8.0 µm where water opacity is highest. The authors model the atmosphere using a water mixing ratio between 10⁻⁴ and 10⁻³ combined with temperature–pressure profiles derived from day‑side emission studies. This model reproduces the observed broadband depths, but the authors emphasize a degeneracy: similar transit depths can be obtained by adjusting either the water abundance or the thermal structure. Consequently, the data alone cannot uniquely constrain the water mixing ratio.

Because IRAC provides only broadband photometry, the presence of other molecules such as CO, CO₂, or CH₄ cannot be definitively assessed. These species have spectral features that overlap with water or lie in narrower wavelength intervals that require spectroscopic resolution. The authors therefore advocate for follow‑up observations with higher spectral resolution instruments (e.g., JWST/NIRSpec, HST/WFC3) to break the water‑abundance/temperature degeneracy and to search for additional absorbers.

In summary, the study demonstrates that precise IRAC photometry, when coupled with rigorous systematic correction and robust statistical fitting, can detect water vapor in the atmosphere of a transiting exoplanet. It also highlights the limitations of broadband photometry for detailed compositional analysis and underscores the need for future spectroscopic measurements to fully characterize the atmospheric chemistry and thermal structure of HD 209458b.