The Third Option: Color Phase Curves to Characterize the Atmospheres of Temperate Rocky Exoplanets

Detecting and characterizing the atmospheres of rocky exoplanets has proven to be challenging for JWST. Transit spectroscopy of the TRAPPIST-1 planets has been impacted by the effects of spots and faculae on the host star. Secondary eclipses have detected hot rocks, but evidence for atmospheres has been difficult to obtain. However, there is a third option that we call color phase curves. This method will apply to synchronously rotating non-transiting planets as well as transiting planets. A color phase curve uses photometry at a long-IR wavelength near the peak of the planetary thermal emission (e.g., 21 microns) divided by photometry at a shorter wavelength where the star dominates more strongly (e.g., 12 microns). We avoid wavelengths having potentially strong molecular absorption (e.g., 15 microns) to minimize degeneracies in the color phase curve, and we aim to detect and characterize the planetary atmosphere via its longitudinal heat transfer. The ratio of two wavelengths observed nearly simultaneously is designed to isolate thermal emission from the planet, discriminate against the star, and largely cancel instrumental systematic effects. Moreover, we show that invoking mass-radius relations, and using self-consistent physical models, will permit the longitudinal heat transfer to be measured independent of the orbital inclination. Radial velocity surveys are detecting many new exoplanets, including temperate rocky worlds with Earth-like masses. Most of those planets will not transit, but color phase curves have the potential to detect and characterize their atmospheres.

💡 Research Summary

The paper introduces “color phase curves” (CPC) as a novel technique for detecting and characterizing the atmospheres of temperate rocky exoplanets, especially those that do not transit. The method relies on simultaneous photometry in two mid‑infrared bands: a long‑wavelength band near the planet’s thermal emission peak (≈21 µm) and a shorter band where stellar emission dominates (≈12.8 µm). By taking the ratio of the fluxes in these two bands, the authors show that planetary thermal emission can be isolated while stellar variability (spots, faculae) and many instrumental systematics largely cancel out.

Key physical arguments: In the Rayleigh‑Jeans limit the Planck function is linear in temperature, so the effect of a temperature perturbation (e.g., a starspot) is wavelength‑independent. Consequently, the flux ratio of a spotted versus an unspotted star is essentially the same at 12.8 µm and 21 µm, differing by only a few hundred parts per million for realistic M‑dwarf spot parameters. This dramatically reduces the “transit light‑source effect” that hampers transit spectroscopy. Stellar flares, which are non‑thermal in the mid‑IR, can be identified by their rapid rise and short duration and either modeled out or excluded from the analysis.

Instrumental systematics are mitigated because both bands are recorded with the same detector array on JWST’s MIRI instrument, allowing many detector‑related drifts to cancel in the ratio. The authors propose a “densified switching” observing strategy: consecutive short exposures in the two filters within a single visit, possibly repeated at strategic orbital phases (e.g., near expected maxima and minima) to maximize signal‑to‑noise.

The scientific payoff hinges on measuring the day‑night temperature contrast (ΔT) of a synchronously rotating planet. In the absence of an atmosphere, the contrast is large, producing a high‑amplitude phase curve; an atmosphere redistributes heat, reducing the amplitude. By combining CPC amplitude with mass–radius relations and self‑consistent radiative‑convective models, the authors demonstrate that ΔT can be inferred even without knowledge of the orbital inclination, a crucial advantage for non‑transiting planets.

A detailed simulation of Proxima Centauri b illustrates feasibility: a 30‑hour campaign (15 h per band) yields a detectable ≈200 ppm flux ratio modulation, sufficient to distinguish an atmosphere from a bare rock. The analysis assumes the planet is tidally locked, the host star’s effective temperature is 3000 K, and spot coverage is 8 % with a 200 K temperature deficit. The resulting spot‑induced ratio variation is only 634 ppm, well below the planetary signal.

The authors discuss broader applicability. Many nearby M‑dwarfs host RV‑detected rocky planets that do not transit; CPC could provide the first atmospheric constraints for these worlds. Limitations include the requirement of tidal locking (or a known resonance), accurate knowledge of planetary radius and albedo to convert flux ratios into temperatures, and the relatively low sensitivity of MIRI at 21 µm, which restricts the method to the brightest nearby systems (≤10 pc). Long‑term stellar variability and residual detector drifts remain sources of uncertainty that will need careful calibration.

In summary, color phase curves offer a promising, under‑exploited avenue for atmospheric characterization of temperate rocky exoplanets, complementing traditional transit and eclipse spectroscopy. By leveraging JWST/MIRI’s stability and the physics of the Rayleigh‑Jeans regime, CPC can isolate planetary thermal emission, mitigate stellar contamination, and measure heat redistribution without requiring a transit geometry. If incorporated into future JWST observing programs, CPC could dramatically expand the sample of rocky exoplanets with atmospheric constraints, especially among the growing catalog of RV‑discovered, non‑transiting worlds.