pith. sign in

arxiv: 2604.07987 · v2 · submitted 2026-04-09 · 🌌 astro-ph.EP

Three-dimensional transport-induced chemistry on temperate sub-Neptune K2-18b, Part II: the combined effects of atmospheric dynamics and chemical reactions

Jiachen Liu , Duncan Christie , Jun Yang , Krisztian Kohary This is my paper

Pith reviewed 2026-05-10 17:21 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords K2-18bsub-Neptunetransport-induced chemistry3D GCMdisequilibrium chemistryvertical mixingJWST spectra
0
0 comments X

The pith

Vertical transport on K2-18b lifts CO and CO2 abundances in the upper atmosphere to ~10^{-3} while equilibrium chemistry predicts levels below 10^{-15}.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper models the atmosphere of the temperate sub-Neptune K2-18b with a three-dimensional general circulation model that couples dynamics to a chemical network at 180 times solar metallicity. Vertical motions driven by the circulation carry molecules upward faster than reactions can restore equilibrium, enriching carbon monoxide and carbon dioxide by many orders of magnitude in the upper layers. Zonal winds then mix these enhanced abundances evenly around the planet at those heights. The photospheric abundances stay similar for rotation periods from synchronous to 10:1 resonance, and the model supplies an effective vertical eddy diffusion coefficient for simpler one-dimensional calculations. Synthetic spectra generated from the three-dimensional fields match existing JWST transmission data at a comparable level.

Core claim

Vertical transport affects the chemical structure significantly, making CO2 and CO more abundant (~10^{-3}) in the upper atmosphere compared to the chemical equilibrium abundance (<10^{-15}), and horizontal winds further homogenize the chemical composition zonally in this region. Molecular abundances in the photosphere generally agree across different rotation periods. A passive tracer yields a one-dimensional equivalent eddy-diffusion coefficient K_zz, and the resulting transmission spectra provide a comparable fit to JWST observations.

What carries the argument

The three-dimensional general circulation model coupled to a chemical reaction network, with a passive tracer used to derive the equivalent vertical eddy diffusion coefficient K_zz.

If this is right

  • Molecular abundances in the photosphere remain consistent across synchronous and asynchronous rotations with 2:1, 6:1, and 10:1 spin-orbit resonances.
  • The derived K_zz value supplies a ready parameter for one-dimensional atmospheric models of similar temperate sub-Neptunes.
  • Synthetic transmission spectra from the three-dimensional fields achieve a fit to JWST observations that is comparable to existing data.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same vertical enrichment mechanism is expected to operate on other cool, hydrogen-rich exoplanets whose chemical timescales exceed dynamical timescales.
  • One-dimensional retrieval analyses of spectra may systematically underestimate carbon dioxide and monoxide if they do not incorporate the transport-induced enhancements reported here.
  • Applying the same modeling approach to planets with different metallicities would test whether the enrichment factor remains roughly constant.
  • The homogenization by zonal winds implies that dayside and nightside spectra should look similar at the altitudes probed by transmission.

Load-bearing premise

The chosen 180 times solar metallicity, the specific chemical network, and the model resolution produce realistic chemical timescales and transport without missing key processes such as clouds or photochemistry.

What would settle it

A direct measurement showing CO or CO2 mixing ratios in the upper atmosphere of K2-18b remain below 10^{-4}, or transmission spectra that cannot be reproduced without assuming chemical equilibrium, would contradict the enrichment result.

Figures

Figures reproduced from arXiv: 2604.07987 by Duncan Christie, Jiachen Liu, Jun Yang, Krisztian Kohary.

Figure 1
Figure 1. Figure 1: Global-mean vertical profiles of CO (a), CO2 (b), NH3 (c), CH4 (d), and H2O (e) mole fractions and passive tracer mass mixing ratio (f). Coloured lines show results from simulations assuming K2-18b is a synchronous rotator (blue), or in 2:1 SOR (orange), 6:1 SOR (purple), and 10:1 SOR (red). Black dashed lines depict the chemical equilibrium abundances calculated via Gibbs free energy minimisation using th… view at source ↗
Figure 2
Figure 2. Figure 2: Horizontal distribution of passive tracer (a), CO (b), CO2 (c), NH3 (d), CH4 (e), and H2O (f) at 0.001 bar. The results are transformed into the heliocentric frame, keeping the substellar point at 0◦ longitude and 0◦ latitude. The black star-shaped markers indicate the location of the substellar point. The streamlines indicate the direction of the flow. To highlight the horizontal gradient, the contours in… view at source ↗
Figure 3
Figure 3. Figure 3: Column-integrated CO and CO2 mole fraction contrast between high- (> 75◦ ) and low- (< 15◦ ) latitude regions around quench level. The vertical axis shows the difference between the high- and low-latitude column abundances normalised by the global mean. Orange squares and blue circles denote CO and CO2, respectively. Positive values indicate that CO and CO2 are more abundant at high-latitude regions. mosph… view at source ↗
Figure 4
Figure 4. Figure 4: 𝐾𝑧𝑧 derived with the flux-gradient relationship I (a) and II (b) and the mixing-length theory (c). Panels (a)–(c) show results from the synchronous simulation, while panels (d)–(f) correspond to the 2:1 SOR simulation. Coloured lines represent the 𝐾𝑧𝑧 profiles calculated every thirty model days from the instantaneous model output , where bluer lines represent earlier stages of the simulation and redder lin… view at source ↗
Figure 5
Figure 5. Figure 5: The spatial distribution of 𝐾𝑧𝑧 estimated using the flux-gradient relationship II (equation 4). (a) The latitude versus pressure distribution of zonal mean 𝐾𝑧𝑧 . (b) The horizontal distribution of 𝐾𝑧𝑧 at 0.001 bar. The black star-shaped marker indicates the location of the substellar point. 𝐾𝑧𝑧 shown in this figure is derived from the synchronous simulation, averaged over the last 3600 days after it spins … view at source ↗
Figure 6
Figure 6. Figure 6: The comparison between the chemical species abundances from 1D ATMO simulations (solid) employing different 𝐾𝑧𝑧 profiles and the global-mean results in the synchronous 3D simulation (dashed). The three panels are results with the average 𝐾𝑧𝑧 profile derived from the flux-gradient relationship I (a), flux-gradient relationship II (b), and the mixing length theory (c) as shown in [PITH_FULL_IMAGE:figures/fu… view at source ↗
Figure 7
Figure 7. Figure 7: K2-18b synthetic transmission spectra predicted by the GCM simulations (coloured lines) compared to JWST NIRISS and NIRSpec data (black and grey error bars) from fig. 3 in Madhusudhan et al. (2023). Panel (a) shows the full (morning plus evening terminators) transmission spectra from chemical kinetics–transport simulations, assuming K2-18b as a synchronous rotator (blue), or has a 2:1 SOR (orange), 6:1 SOR… view at source ↗
Figure 8
Figure 8. Figure 8: K2-18b synthetic transmission spectrum predicted by the GCM simulation compared to JWST NIRISS and NIRSpec data (grey and black error bars) from fig. 3 in Hu et al. (2025). Panel (a) shows the full (morning plus evening terminators) transmission spectrum from the synchronous chemical kinetics–transport simulation as an example. Panel (b) shows the contributions from the dominant chemical species, CO, CO2, … view at source ↗
read the original abstract

The upper atmospheres of temperate sub-Neptunes are strongly influenced by atmospheric dynamics due to their cool equilibrium temperature and thereby longer chemical timescales than the atmospheric dynamical timescales. In this study, we used a three-dimensional (3D) general circulation model to investigate the transport-induced disequilibrium chemistry and vertical mixing on temperate gas-rich mini-Neptunes, using K2-18b as an example. We model K2-18b assuming 180 times solar metallicity and consider it as either a synchronous or an asynchronous rotator, exploring spin-orbit resonances of 2:1, 6:1, and 10:1. We find that the vertical transport affects the chemical structure significantly, making CO$_2$ and CO more abundant ($\sim$10$^{-3}$) in the upper atmosphere compared to the chemical equilibrium abundance (<10$^{-15}$), and horizontal winds further homogenize the chemical composition zonally in this region. Molecular abundances in the photosphere generally agree across different rotation periods. We employ a passive tracer in the model to estimate the one-dimensional (1D) equivalent eddy-diffusion coefficient ($K_{zz}$) of K2-18b, providing a parameter useful for future 1D atmospheric models. Additionally, synthetic transmission spectra generated from our model are compared with the JWST observations, and we find that our model can provide a comparable fit to the observations. This work offers a 3D perspective on transport-induced chemistry on a temperate sub-Neptune and derives vertical mixing parameters to support 1D modelling.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper uses 3D GCM simulations of temperate sub-Neptune K2-18b at 180x solar metallicity for synchronous and asynchronous rotators (2:1, 6:1, 10:1 resonances) to study transport-induced disequilibrium chemistry. Vertical transport is shown to quench CO2 and CO to ~10^{-3} in the upper atmosphere (vs. chemical equilibrium <10^{-15}), with horizontal winds providing zonal homogenization; photospheric abundances agree across periods, a 1D-equivalent K_zz is derived from a passive tracer, and synthetic transmission spectra provide a comparable fit to JWST data.

Significance. If the central quenching and homogenization results hold, the work supplies a useful 3D perspective on dynamics-chemistry coupling for temperate sub-Neptunes and a practical K_zz value for 1D models, directly aiding JWST interpretation. The strength is the explicit comparison of 3D transport timescales to chemical equilibrium; however, the significance is reduced by the absence of photochemistry and limited validation of the reported abundance contrasts.

major comments (2)
  1. [Abstract/Methods] Abstract and methods description: the central claim that vertical transport quenches CO2 and CO to ~10^{-3} (versus equilibrium <10^{-15}) rests on a chemical network that excludes photolysis reactions. At the low pressures and cool temperatures of the upper atmosphere, UV-driven dissociation timescales for CO2 and CO can be comparable to or shorter than the GCM dynamical mixing times, so the reported disequilibrium abundances are computed under an incomplete kinetic model that may artificially preserve the quenched values.
  2. [Results] Results section (abundance and spectral comparisons): no error bars, sensitivity tests, or quantitative metrics are provided for the ~10^{-3} vs. <10^{-15} contrast, nor for the statement that 'molecular abundances in the photosphere generally agree across different rotation periods.' Details on post-processing, chemical network choices, and how they affect the contrast are also absent, undermining the robustness of the cross-rotation and JWST-fit claims.
minor comments (2)
  1. [Abstract] The abstract states that the model 'can provide a comparable fit' to JWST observations but does not specify which spectral features, wavelength range, or goodness-of-fit metric is used.
  2. [Methods/Results] Notation for the derived K_zz should be clarified (e.g., whether it is vertically averaged or pressure-dependent) to aid use in future 1D models.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed review of our manuscript on transport-induced disequilibrium chemistry for K2-18b. We address each major comment below and have prepared revisions to improve clarity, robustness, and completeness where possible.

read point-by-point responses

  1. Referee: [Abstract/Methods] Abstract and methods description: the central claim that vertical transport quenches CO2 and CO to ~10^{-3} (versus equilibrium <10^{-15}) rests on a chemical network that excludes photolysis reactions. At the low pressures and cool temperatures of the upper atmosphere, UV-driven dissociation timescales for CO2 and CO can be comparable to or shorter than the GCM dynamical mixing times, so the reported disequilibrium abundances are computed under an incomplete kinetic model that may artificially preserve the quenched values.

    Authors: We acknowledge that the chemical network employed is restricted to thermal reactions and omits photolysis, as the study centers on the coupling between atmospheric dynamics and transport-induced quenching. At the low temperatures of K2-18b, thermal chemical timescales for CO and CO2 are indeed long compared to dynamical timescales, which underpins the reported quenching. However, we agree that photolysis could modify upper-atmosphere abundances. In the revised manuscript we will (i) explicitly state in the abstract and methods that results are based on a thermal kinetic network, (ii) add a discussion paragraph comparing photolysis and dynamical timescales with order-of-magnitude estimates, and (iii) note photochemistry as a limitation for future modeling. These changes clarify the scope without altering the core transport findings. revision: partial

  2. Referee: [Results] Results section (abundance and spectral comparisons): no error bars, sensitivity tests, or quantitative metrics are provided for the ~10^{-3} vs. <10^{-15} contrast, nor for the statement that 'molecular abundances in the photosphere generally agree across different rotation periods.' Details on post-processing, chemical network choices, and how they affect the contrast are also absent, undermining the robustness of the cross-rotation and JWST-fit claims.

    Authors: We agree that additional quantitative support and methodological detail will strengthen the presentation. The revised manuscript will include: error bars on abundance profiles derived from GCM temporal and spatial variability; quantitative metrics (e.g., standard deviation or range) demonstrating photospheric abundance agreement across rotation periods; expanded methods text describing post-processing steps for spectra and the specific chemical network (including key reactions and references); and a short sensitivity discussion on how network assumptions influence the reported contrasts. These additions will better substantiate the robustness of the cross-rotation results and JWST spectral comparisons. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper's central results on CO2/CO abundances (~10^{-3} vs. equilibrium <10^{-15}) and zonal homogenization are direct outputs of a 3D GCM coupled to a chemical network, with vertical transport and horizontal winds applied as dynamical forcings. The K_zz value is obtained as a post-hoc diagnostic by injecting a passive tracer into the same 3D simulation and extracting an equivalent 1D eddy coefficient; this diagnostic is not fed back into the chemistry integration or used to redefine the reported abundances. No self-citations are invoked to justify uniqueness or to close a definitional loop, and the model is compared against external JWST data without parameter tuning that would force the disequilibrium result. The derivation therefore remains self-contained and independent of its own fitted quantities.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claims rest on the assumption that the chosen metallicity, chemical network, and GCM dynamics produce realistic disequilibrium without additional processes; no new entities are postulated.

free parameters (2)
  • metallicity = 180x solar
    Set to 180 times solar; directly controls equilibrium abundances and therefore the contrast with transport-enhanced values.
  • spin-orbit resonance ratios
    Tested at 2:1, 6:1, 10:1; chosen to explore asynchronous rotation effects.
axioms (1)
  • domain assumption Chemical timescales exceed dynamical timescales in the upper atmosphere of temperate sub-Neptunes
    Invoked to justify why transport-induced disequilibrium dominates over equilibrium chemistry.

pith-pipeline@v0.9.0 · 5596 in / 1469 out tokens · 55556 ms · 2026-05-10T17:21:35.433536+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

  1. [1]

    S., et al., 2016, A&A, 595, A36 BaeyensR.,KoningsT.,VenotO.,CaroneL.,DecinL.,2022,MNRAS,512, 4877 Barrier E

    Amundsen D. S., et al., 2016, A&A, 595, A36 BaeyensR.,KoningsT.,VenotO.,CaroneL.,DecinL.,2022,MNRAS,512, 4877 Barrier E. F. L., Madhusudhan N., 2025, MNRAS, 538, 2463 Batalha N. M., 2014, PNAS, 111, 12647 Batygin K., 2018, AJ, 155, 178 Benneke B., et al., 2019, ApJ, 887, L14 Benneke B., et al., 2024, arXiv preprint arXiv:2403.03325 Blain D., Charnay B., B...

    work page doi:10.5281/zenodo.8305232 2016

  2. [2]

    Cycle 1, ID

    ed RS Stepleman et al Hu R., Damiano M., 2021, Deep Characterization of the Atmosphere of a Temperate Sub-Neptune, JWST Proposal. Cycle 1, ID. #2372 Hu R., Damiano M., Scheucher M., Kite E., Seager S., Rauer H., 2021, ApJ, 921, L8 Hu R., et al., 2025, arXiv e-prints, p. arXiv:2507.12622 Hunten D. M., 1975, in McCormac B. M., ed., , Atmospheres of Earth an...

    work page doi:10.1007/978-94-010-1799-2_4 2021

  3. [3]

    Results are derived from simulated data over the first 100 d, in which the horizontal temperature contrast begins to form, with an output interval of 10 d

    (d) 0.990 0.995 1.000 1.005 1.010 5.0 2.5 0.0 2.5 5.0 1.0 0.5 0.0 0.5 1.0 5.0 2.5 0.0 2.5 5.0 Eddy Heat Transport Divergence Components Figure B2.Pressure versus latitude distribution of air temperature contrast (a), transient eddy heat transport (b), stationary eddy heat transport (c), and total eddy heat transport (d) for the 6:1 SOR simulation. Results...

    work page 2026