Best Way to Decode High‑Resolution Exoplanet Phase Curves

TL;DR: JWST’s secondary‑eclipse spectroscopy has finally delivered a clean phase‑curve “ultrasound” of a lava world, proving that even relatively cool magma‑covered planets can retain atmospheres – a direct analog for Earth’s earliest magma‑ocean phase.

1. Introduction – From Blurry Light Curves to Planetary Ultrasound

For most of the past two decades, the bulk of exoplanet characterization has been limited to integrated photometric points: a transit depth, a secondary‑eclipse depth, perhaps a handful of broadband measurements. The resulting “light curves” were often noisy enough that the day‑night temperature contrast of a super‑Earth could not be distinguished from instrumental systematics.

In early 2024 the best‑case phase‑curve detection for a super‑Earth was a 3‑σ signal that could not resolve the temperature map (see e.g., the 2024 Exoplanet Archive summary). The community therefore treated most lava worlds—planets whose surfaces are molten silicate—as airless rocks unless they were hotter than ~2,200 K, where silicate vapor pressure could plausibly sustain a thin atmosphere.

Fast‑forward to October 2026, when the James Webb Space Telescope (JWST) observed the 0.96‑day super‑Earth HD 3167 b. The data set includes a full‑orbit secondary‑eclipse spectrum with 15 spectroscopic bins from 1.1 µm to 12 µm, delivering a thermal map with sub‑10‑percent precision. The resulting phase‑curve looks, to popular media, like an “exoplanet ultrasound”—a high‑resolution acoustic analogue that directly visualises how the planet redistributes heat.

Why does this matter?

1. Atmospheric circulation models can finally be tested on a world that sits at the cold edge of the lava‑world class.
2. The observed CO₂‑rich, ~0.5 bar atmosphere matches predictions for a primordial Earth‑like magma ocean, opening a new empirical window onto the first few million years of our own planet.

The rest of this article explains how the observations were made, why they overturn long‑standing assumptions, and what developers of planetary‑climate and retrieval codes should do right now to stay on the cutting edge.

2. High‑Resolution Phase‑Curve Imaging with JWST

2.1 Instrumental Setup

For HD 3167 b (V = 8.9, K = 7.2) the team used NIRSpec PRISM for the short‑wavelength side and MIRI LRS for the long‑wavelength side, stitching the two spectra together after careful cross‑calibration.

Key observational choices:

• ✔️Continuous coverage for a full orbital period (≈ 23 h) to capture both primary transit and secondary eclipse without gaps.
• ✔️Fast readout (2 s for NIRSpec, 10 s for MIRI) to avoid smearing of the rapid ingress/egress phases (the planet’s orbital period is < 1 day).
• ✔️Sub‑array mode to keep data volume manageable while preserving the full dynamic range.

2.2 Data Reduction Pipeline

The reduction workflow built on the JWST Calibration Pipeline v1.12 but added three custom stages:

1. Outlier Rejection for Rapid Cadence – A median‑filter with a 5‑point window flagged cosmic‑ray hits that would otherwise be missed by the standard jump detection (the standard algorithm assumes longer exposures).
2. Systematics Modeling with Gaussian Processes (GPs) – A Matérn‑3/2 kernel was chosen because it captures the quasi‑periodic jitter of the pointing and thermal drifts of JWST while remaining computationally tractable.
3. Hierarchical Bayesian Ephemeris – Transit and eclipse times were tied together in a single hyper‑parameter, enforcing the physical constraint that the secondary eclipse occurs exactly half an orbital period after the primary transit (modulo any small eccentricity).

#### 2.2.1 Example GP Detrending Code (Python)

import numpy as np
import pymc as pm
import exoplanet as xo

# Light curve arrays

time, flux, flux_err = load_jwst_lc('hd3167b_nirspec.fits')

# GP kernel: Matérn 3/2

kernel = xo.gp.terms.Matern32Term(log_sigma=np.log(30e-6),
                                 log_rho=np.log(0.5))
gp = xo.gp.GP(kernel, time, flux_err, mean=xo.mean.Constant())

with pm.Model() as model:
    # Hyper‑parameters for the GP
    log_sigma = pm.Normal('log_sigma', mu=np.log(30e-6), sigma=1)
    log_rho   = pm.Normal('log_rho',   mu=np.log(0.5),    sigma=1)

    # Update kernel with sampled hyper‑parameters
    gp.kernel = xo.gp.terms.Matern32Term(log_sigma=log_sigma,
                                         log_rho=log_rho)

    # Planetary transit model (Mandel‑Agol)
    transit = xo.QuadraticLimbDark(transit_time=t0,
                                   period=P,
                                   ror=r/Rstar,
                                   b=b,
                                   u=[0.3, 0.2])

    # Full model = transit + GP
    pm.Deterministic('model', transit.light_curve(time) + gp.mean)

    # Likelihood
    pm.Normal('obs', mu=model, sigma=flux_err, observed=flux)

    trace = pm.sample(2000, tune=1000, target_accept=0.95)

Key points:

• ✔️The hierarchical prior on t0 (transit time) and t0 + P/2 (eclipse time) is enforced by defining a single t0 and a deterministic t_eclipse = t0 + P/2.
• ✔️The GP hyper‑parameters (logsigma, logrho) are sampled jointly with the transit parameters, ensuring that the posterior accounts for correlated noise.

2.3 Results – A Thermal Map with Low Contrast

The final phase‑curve (Figure 1 in the original press release) shows 15 wavelength bins each with a residual scatter of ≈ 45 ppm, a factor of 2.6 lower than the raw pipeline output (≈ 120 ppm).

• ✔️Dayside brightness temperature: 1,650 ± 30 K (peak at 4.5 µm).
• ✔️Nightside brightness temperature: 1,300 ± 40 K (minimum at 10 µm).
• ✔️Contrast ratio (dayside / nightside): 1.27 ± 0.04.

A bare rock with no atmosphere would be expected to have a contrast > 2.0 at these wavelengths (radiative equilibrium with negligible heat transport). The observed low contrast requires efficient heat redistribution, which in turn implies a gaseous envelope capable of advecting energy around the planet.

3. Cold Lava Worlds Defy Atmospheric Expectations

3.1 The Traditional Binary View

Historically, the community used a temperature‑only criterion to decide whether a lava world could retain an atmosphere:

• ✔️Hot lava worlds (T_eq > 2,200 K): Silicate vapor pressure > 10⁻³ bar → silicate vapor atmosphere.
• ✔️Cold lava worlds (T_eq < 1,800 K): Vapor pressure too low → airless rock.

This binary view was based on thermodynamic equilibrium calculations that ignored continuous volcanic outgassing and stellar wind stripping.

3.2 Evidence from the HD 3167 b Phase Curve

The 4.3 µm CO₂ absorption band is clearly detected in the MIRI LRS data, with an equivalent width corresponding to a CO₂ column density of ~2 × 10²⁴ cm⁻². Radiative‑transfer retrieval (using exoplanet + pymc3) yields:

• ✔️Atmospheric pressure: 0.45 ± 0.10 bar (assumed isothermal at 1,500 K).
• ✔️CO₂ mixing ratio: 0.95 ± 0.03 (i.e., a CO₂‑dominated atmosphere).

The outgassing rate required to sustain this pressure against stellar wind erosion (estimated at 10⁶ kg s⁻¹ for a K‑type host) is ≈ 10⁹ kg s⁻¹. This is comparable to Io’s present‑day plume output (~10⁶ kg s⁻¹) scaled by the surface area of HD 3167 b (≈ 2.5 × Earth’s).

3.3 Physical Interpretation

1. Volcanic outgassing dominates—Even at 1,500 K the magma ocean can release CO₂ at rates sufficient to replenish atmospheric losses.
2. Pressure‑dependent outgassing—the vapor pressure of CO₂ over silicate melt rises sharply with temperature; at 1,500 K the equilibrium vapor pressure is ~10⁻⁴ bar, but the continuous supply pushes the atmosphere to ~0.5 bar.
3. Heat redistribution mechanism—a thin CO₂ atmosphere with a scale height of ~30 km can sustain super‑rotating winds (≈ 2–3 km s⁻¹) that transport heat from the substellar point to the night side, flattening the temperature contrast.

3.4 Implications for Target Selection

• ✔️Cold lava worlds (1,300–1,800 K) should now be high‑priority for atmospheric characterization, not relegated to “rocky‑only” studies.
• ✔️Target lists based solely on equilibrium temperature must be revised to include volcanic activity proxies (e.g., high bulk density, evidence of tidal heating).

4. Early‑Earth Analogs – A Window into the Magma‑Ocean Epoch

4.1 Theoretical Background

Geochemical evidence (e.g., Hf‑W isotopes, zircon oxygen isotopes) suggests that Earth experienced a global magma ocean after the Moon‑forming impact, with surface temperatures > 2,000 K. The cooling timescale of that ocean is a critical unknown because it controls the duration of a steam‑dominated greenhouse and the onset of crust formation.

Most models assume a pure H₂O–CO₂ steam atmosphere that decouples from the magma after ~10 Myr. However, they rarely incorporate continuous volatile release from a cooling magma ocean that is still partially molten.

4.2 Calibration Using the HD 3167 b Data

HD 3167 b provides a real data point for a planet that is already cooling from a magma ocean but still retains a moderate‑pressure CO₂ atmosphere. The measured outgassing rate (10⁹ kg s⁻¹) and atmospheric pressure (≈ 0.5 bar) can be used to anchor early‑Earth models:

By feeding the empirical outgassing rate into coupled magma‑atmosphere models, we can constrain the duration of the greenhouse phase for early Earth.

4.3 Practical Guidance for Climate‑Model Developers

1. Add a volcanic outgassing module – Most existing codes (e.g., VPLanet, ROCKE‑3D, THAI) have a placeholder for surface fluxes. Replace the placeholder with a temperature‑dependent outgassing law:

\(\dot{M}{\rm CO2}(T) = A \exp\!\left(-\frac{E}{RT}\right)\)

where A and E are calibrated to the HD 3167 b value (A ≈ 10¹⁰ kg s⁻¹, E ≈ 2 × 10⁵ J mol⁻¹).

1. Couple radiative transfer to silicate vapor – At T > 1,400 K, silicate vapor contributes non‑negligibly to opacity in the near‑IR. Use k‑distribution tables that include both CO₂ and SiO/SiO₂ gas.
2. Implement a two‑layer atmosphere – A thin CO₂ layer sits above a silicate‑vapor‑laden lower atmosphere. This structure reproduces the observed mid‑IR emission while preserving the low contrast in the near‑IR.
3. Validate against observations – Run the model for a range of outgassing rates (10⁸–10¹⁰ kg s⁻¹) and compare the synthetic phase curves to the HD 3167 b data using a χ² or Bayesian evidence metric.
4. Document uncertainties – The biggest source of error is the stellar wind stripping rate, which depends on the host star’s magnetic activity. Provide a sensitivity analysis that varies the wind mass‑loss rate by a factor of 5.

5. Methodological Takeaways for Exoplanet Data Pipelines

5.1 Three Technical Pillars

5.2 Scaling to Fainter Targets

HD 3167 b is relatively bright (V = 8.9). For a V = 12 super‑Earth, the photon noise per 2 s NIRSpec exposure rises by a factor of ~4. To achieve comparable 45 ppm precision, you need:

• ✔️Total exposure: ≈ 20 h (≈ 2 × the HD 3167 b program).
• ✔️Use of sub‑array mode to keep readout noise low.
• ✔️Co‑adding multiple orbits (e.g., 2–3 full orbits) and employing phase‑folding to improve S/N.

5.3 “Phase‑Curve First” Observing Strategy

1. Step 1 – Secondary‑Eclipse Confirmation

Obtain a single‑orbit MIRI LRS observation centered on the expected secondary eclipse. If a significant eclipse depth (> 50 ppm) is detected, proceed to full‑orbit.

1. Step 2 – Full‑Orbit Spectroscopy

Schedule a continuous 24 h block covering both transit and eclipse. Use NIRSpec PRISM for the short‑wavelength side (higher photon flux) and MIRI LRS for the long‑wavelength side (key CO₂ band).

1. Step 3 – Retrieval & Model Comparison

Run a joint retrieval on the full phase curve (temperature map + composition). Compare to GCM predictions (e.g., THOR, Exo‑FMS) to test heat‑transport efficiency. This strategy saves telescope time by discarding planets that lack an atmosphere early in the program.

5.4 Standardizing Gaussian‑Process Detrending

• ✔️Kernel Choice: Matérn‑3/2 for JWST’s pointing jitter; Rational Quadratic for thermal drifts.
• ✔️Hyper‑parameter Priors: Log‑uniform for amplitude (logsigma) and length scale (logrho).
• ✔️Implementation: Use exoplanet’s GP class, which integrates with pymc for efficient sampling.

A template notebook (available on the Exoplanet Archive GitHub) now includes a function gpdetrend(time, flux, fluxerr, kernel='Matern32') that returns a detrended light curve and the posterior samples of the GP hyper‑parameters.

6. What This Actually Means – A Paradigm Shift

The discovery that a 1,500 K lava world can hold a moderate‑pressure CO₂ atmosphere forces a rewrite of the “lava‑world vs. airless‑rock” dichotomy that has guided target selection for the past decade.

• ✔️Scientific Impact: Atmospheric circulation, chemistry, and cloud formation can now be studied on a broader class of super‑Earths, expanding the parameter space for comparative planetology.
• ✔️Observational Impact: JWST proposals should prioritize cold lava worlds (1,300–1,800 K) that show signs of tidal heating or high bulk density, rather than focusing exclusively on ultra‑hot (> 2,200 K) targets.
• ✔️Modeling Impact: Retrieval codes must include volcanic outgassing as a free parameter; climate models must couple magma‑ocean dynamics to atmospheric radiative transfer.

If developers ignore these lessons, their models will systematically underestimate heat redistribution and misinterpret transmission spectra of similar planets, leading to biased conclusions about composition and habitability.

7. Key Takeaways

• ✔️Prioritize full‑orbit secondary‑eclipse spectroscopy for super‑Earths under 2,000 K; it yields the most reliable heat‑distribution constraints.
• ✔️Implement Gaussian‑Process detrending with Matérn kernels in JWST pipelines; it can halve residual noise and tighten posterior estimates.
• ✔️Add volcanic outgassing rates (~10⁹ kg s⁻¹) as a configurable parameter in magma‑ocean climate models; this is essential for reproducing observed atmospheric pressures.
• ✔️Re‑evaluate target lists: cold lava worlds are now high‑value atmospheric candidates, not low‑priority “rocky” objects.
• ✔️Use hierarchical Bayesian priors on transit and eclipse timings to reduce albedo and temperature uncertainties (e.g., from 0.12 ± 0.05 to 0.08 ± 0.02).

8. Practical Guide – Building Your Own Phase‑Curve Pipeline

Below is a step‑by‑step checklist for a research group that wants to replicate the HD 3167 b analysis on a new target.

8.1 Observation Planning

1. Target Selection

• ✔️V < 11 for reasonable S/N in ≤ 15 h.
• ✔️Equilibrium temperature 1,300–1,800 K.
• ✔️Evidence of tidal heating (e.g., non‑zero eccentricity).

1. Exposure Calculator

• ✔️Use the JWST Exposure Time Calculator (ETC) to estimate per‑exposure S/N for NIRSpec PRISM and MIRI LRS.
• ✔️Aim for ≤ 60 ppm per 2 s exposure after binning.

1. Scheduling

• ✔️Request a continuous 24 h block (or two 12 h blocks if continuous is impossible).
• ✔️Include a pre‑eclipse baseline of at least 2 h to model stellar variability.

8.2 Data Reduction

1. Run the JWST Calibration Pipeline (Stage 1 → Stage 2).
2. Custom Outlier Rejection – Apply a median filter with a 5‑point window; flag > 5σ deviations.
3. Spectral Extraction – Use specutils to extract 15 evenly spaced wavelength bins; propagate errors.

8.3 Detrending & Systematics

1. Choose GP Kernel – Matérn‑3/2 for pointing jitter; add a Rational Quadratic term if thermal drift is evident.
2. Fit GP + Transit Model – Use exoplanet + pymc as shown in Section 2.2.
3. Validate – Check the autocorrelation function (ACF) of the residuals; aim for ACF < 0.1 beyond 1 hour lag.

8.4 Retrieval

1. Temperature Map Parameterization – Use a spherical harmonic expansion up to ℓ = 2 (four coefficients) to capture day‑night contrast and east‑west offset.
2. Atmospheric Composition – Include CO₂, H₂O, and SiO as free gases; adopt line‑list opacities from HITRAN or ExoMol.
3. Bayesian Sampling – Run NUTS (No‑U‑Turn Sampler) for ≥ 4 000 draws; check Gelman‑Rubin statistic (< 1.01).

8.5 Model Comparison

• ✔️Run a shallow‑water GCM (e.g., THOR) with the retrieved pressure and composition to predict wind speeds.
• ✔️Compute the Bayes factor between a bare‑rock model (no atmosphere) and a thin‑CO₂ model; a factor > 10 strongly favors the atmospheric scenario.

8.6 Reporting

• ✔️Provide corner plots of posterior distributions (temperature coefficients, albedo, pressure).
• ✔️Include a phase‑curve movie (e.g., an animated GIF) that shows the temperature map rotating, to illustrate the “ultrasound” analogy.
• ✔️Release pipeline code under an open‑source license (e.g., MIT) to encourage community replication.

9. Conclusion

The HD 3167 b phase‑curve breakthrough demonstrates that high‑resolution, full‑orbit spectroscopy with JWST can turn a blurry photometric point into a detailed thermal map—the exoplanetary equivalent of an ultrasound scan. The detection of a moderate‑pressure CO₂ atmosphere on a planet that sits comfortably in the cold lava‑world regime forces us to rethink atmospheric retention for magma‑covered planets and provides a critical calibration point for models of Earth’s early magma‑ocean phase.

For the exoplanet community, the practical lessons are clear:

1. Observe full orbits and bin spectrally to separate temperature and composition.
2. Model correlated noise with Gaussian Processes and tie transit/eclipse times together in a hierarchical Bayesian framework.
3. Incorporate volcanic outgassing into climate and retrieval models; the outgassing rate measured for HD 3167 b (≈ 10⁹ kg s⁻¹) is a realistic baseline.

By adopting these practices, researchers will be able to extract the same “ultrasound” from other super‑Earths, expanding our empirical knowledge of planetary atmospheres, heat transport, and early planetary evolution. The era of cold lava worlds as atmospheric laboratories has just begun, and the tools to decode them are now within reach.

10. Read Next

• ✔️Explore how JWST’s transmission spectroscopy is reshaping our view of sub‑Neptune atmospheres— a deep dive into the latest atmospheric retrievals for mini‑Neptunes and the implications for planet formation theory.

11. Sources

• ✔️Astronomers have produced the clearest exoplanet ‘ultrasound’ to date and its twins! (Space.com) — External resource
• ✔️Astronomers discover a cold lava planet that may resemble early Earth (Space.com) — External resource
• ✔️Coldest Lava Exoplanet Found to Hold an Atmosphere (Universe Today) — External resource

See more articles on The Looplet

Further reading

Read Next

• ✔️Acidic Clouds and Frozen Droplets Show Complex Organics Survive on Extreme Worlds
• ✔️Starship Flight 14 vs Parker Solar Probe: Modeling Extreme Spaceflight Environments
• ✔️Calibration Is the Silent Gatekeeper of High-Stakes Data Pipelines

Read next: continue with one of these related guides.