3D Transport-Induced Disequilibrium Chemistry

• Three-dimensional transport-induced disequilibrium chemistry is the process where atmospheric dynamics outpace local chemical kinetics, freezing chemical abundances away from equilibrium.
• State-of-the-art 3D GCMs couple advection, diffusion, and chemical kinetics to accurately predict quench levels and their impact on observable spectral features.
• Observable signatures include enhanced opacities in CO, CH4, and NH3 bands, with transport processes modifying temperature profiles and radiative feedbacks.

Three-dimensional (3D) transport-induced disequilibrium chemistry refers to the set of processes by which large-scale atmospheric motions (winds, convection, circulation) drive the chemical composition of an atmospheric environment—such as a giant planet, brown dwarf, sub-Neptune, or exoplanet—away from its local thermochemical equilibrium state. In this regime, the timescales for dynamical mixing in multiple spatial directions are shorter than the characteristic timescales for chemical reactions to restore equilibrium, resulting in prominent departures from local equilibrium that are inherently three-dimensional in structure and have profound consequences for observed spectra and atmospheric thermal profiles.

1. Fundamental Principles and Timescale Hierarchies

The core mechanism underlying 3D transport-induced disequilibrium chemistry is the competition between dynamical transport and chemical kinetics. The evolution of the abundance of a chemical species $n_i(x,y,z,t)$ in an atmosphere is governed by the advection-diffusion-reaction equation: $$\frac{\partial n_i}{\partial t} + \nabla\cdot(n_i \mathbf{u}) = \nabla\cdot(D \nabla n_i) + P_i - L_i n_i$$ where $\mathbf{u}$ is the local 3D wind field, $D$ is an eddy-diffusion tensor (with components $K_{zz}$, $K_h$ for vertical/horizontal mixing), and $P_i$, $L_i$ are production and loss rates from chemical reactions.

A critical control parameter is the local ratio of transport timescale $\tau_{\rm dyn}$ and chemical timescale $\tau_{\rm chem}$.

• The characteristic horizontal and vertical advection timescales for planetary atmospheres are $$\frac{\partial n_i}{\partial t} + \nabla\cdot(n_i \mathbf{u}) = \nabla\cdot(D \nabla n_i) + P_i - L_i n_i$$0 and $$\frac{\partial n_i}{\partial t} + \nabla\cdot(n_i \mathbf{u}) = \nabla\cdot(D \nabla n_i) + P_i - L_i n_i$$1, where $$\frac{\partial n_i}{\partial t} + \nabla\cdot(n_i \mathbf{u}) = \nabla\cdot(D \nabla n_i) + P_i - L_i n_i$$2 is a planetary scale, $$\frac{\partial n_i}{\partial t} + \nabla\cdot(n_i \mathbf{u}) = \nabla\cdot(D \nabla n_i) + P_i - L_i n_i$$3 is horizontal wind speed, $$\frac{\partial n_i}{\partial t} + \nabla\cdot(n_i \mathbf{u}) = \nabla\cdot(D \nabla n_i) + P_i - L_i n_i$$4 is atmospheric scale height, and $$\frac{\partial n_i}{\partial t} + \nabla\cdot(n_i \mathbf{u}) = \nabla\cdot(D \nabla n_i) + P_i - L_i n_i$$5 is vertical velocity (e.g., for HD 189733b, $$\frac{\partial n_i}{\partial t} + \nabla\cdot(n_i \mathbf{u}) = \nabla\cdot(D \nabla n_i) + P_i - L_i n_i$$6 s, $$\frac{\partial n_i}{\partial t} + \nabla\cdot(n_i \mathbf{u}) = \nabla\cdot(D \nabla n_i) + P_i - L_i n_i$$7--$$\frac{\partial n_i}{\partial t} + \nabla\cdot(n_i \mathbf{u}) = \nabla\cdot(D \nabla n_i) + P_i - L_i n_i$$8 s for the photosphere (Steinrueck et al., 2018)).
• The chemical timescale for interconversion between, e.g., CO/CH$$\frac{\partial n_i}{\partial t} + \nabla\cdot(n_i \mathbf{u}) = \nabla\cdot(D \nabla n_i) + P_i - L_i n_i$$9 is $\mathbf{u}$0, with $\mathbf{u}$1 and $\mathbf{u}$2 as Arrhenius-law rate coefficients, potentially $\mathbf{u}$3 s in photospheric conditions.
• The "quench level" (altitude or pressure) is defined where $\mathbf{u}$4; above this point, abundances are "frozen" at the quenched value, and local temperature or pressure variations cannot restore equilibrium (Mendonça et al., 2018, Bordwell et al., 2018).

The transport process acts in all dimensions: longitudinal (zonal) mixing by equatorial superrotating jets, meridional (latitudinal) mixing by overturning circulation and eddies, and vertical mixing by convection and turbulence. The Damköhler number, $\mathbf{u}$5, quantifies the dynamical vs. chemical control: Da $\mathbf{u}$6 implies transport dominance and quenching (Zamyatina et al., 2022).

2. Modeling Approaches and Chemical Kinetics Integration

State-of-the-art 3D general circulation models (GCMs) now couple non-hydrostatic dynamical cores to radiative transfer modules and explicit chemical kinetics solvers. Chemical networks, e.g., the 30-species Venot reduced network, are integrated via operator splitting: advection and mixing steps are followed by local kinetic updates for each grid cell, with rate coefficients $\mathbf{u}$7 sourced from laboratory or theoretical studies (Drummond et al., 2020, Liu et al., 9 Apr 2026).

A common approach, justified by timescale arguments, is the "chemical relaxation" or "quench" scheme, in which selected tracers (e.g., CO, CH$\mathbf{u}$8) are advected by the resolved flow, subject to a chemical source/sink term that exponentially relaxes toward the local equilibrium value on a prescribed $\mathbf{u}$9 (Mendonça et al., 2018, Bell et al., 2024).

For parameterization and comparison to 1D column models, the effective vertical eddy diffusivity $D$0 can be diagnosed a posteriori from GCM passive tracer fields via flux–gradient relationships: $D$1 or directly from mixing-length theory $D$2 (Liu et al., 9 Apr 2026, Liu et al., 30 Jun 2025).

A key advance is the identification of the "chemical scale height" $D$3, derived from the vertical profiles of reaction rate, density, and equilibrium abundance, leading to improved quench-point estimates over traditional density-scale-height-based prescriptions (Bordwell et al., 2018).

3. Atmospheric Regimes and Spatial Structure of Disequilibrium

The manifestation of 3D transport-induced disequilibrium chemistry depends strongly on planetary regime:

• Hot Jupiters: Strong horizontal winds ($D$4–3 km/s) and convection homogenize CO/CH$D$5 abundances above the quench level, often rendering the dayside CO-dominated and the nightside, contrary to equilibrium expectations, also CO-dominated due to advection, with latitudinal/polar isolation possible in cold vortices (Steinrueck et al., 2018, Mendonça et al., 2018).
• Temperate sub-Neptunes and mini-Neptunes: Longer chemical timescales and deep detached convective zones (e.g., 1–5 bar) set the quench level deep in the atmosphere; vertical transport sustains super-equilibrium CO/CO$D$6 at upper levels, while strong zonal flows homogenize compositions longitudinally. Rotation rate modulates the meridional gradients, especially at high latitudes via transient eddies (Liu et al., 9 Apr 2026, Liu et al., 30 Jun 2025).
• Ice Giants (Uranus, Neptune): Deep mixing and slow kinetics produce substantial equator–pole gradients in quenched species (e.g., CH$D$7), controlled by large-scale overturning circulation; $D$8 values and internal heat flux govern the vertical gradients and depth of quench (Moses et al., 2020).
• Porous media and solid–liquid interfaces: At Darcy and pore scales, incomplete mixing and "lamella folding" produce strong local concentration gradients, yielding velocity- and mixing-rate-dependent departures from complete-mixing predictions, requiring upscaling models that resolve the chemical interface structure (Izumoto et al., 2023).

Spatially, the post-quench distribution exhibits nearly uniform abundance above $D$9 except where latitudinal mixing is weak (cold poles, vortices). In GCMs, longitude–latitude maps show "rings" of uniformity broken only by residual 3D circulation barriers (Zamyatina et al., 2022, Drummond et al., 2020).

4. Spectroscopic and Observational Implications

Transport-induced disequilibrium chemistry produces observational signatures in both transmission and emission spectra, as well as in phase curves:

• Transit spectra: Enhanced CO/CH$K_{zz}$0/NH$K_{zz}$1 abundances from vertical quenching deepen spectral features at their band centers, e.g., CH$K_{zz}$2 at 3.3 $K_{zz}$3m, CO$K_{zz}$4 at 4.3 $K_{zz}$5m, NH$K_{zz}$6 at 10 $K_{zz}$7m. Transmission depth changes of 60–300 ppm are predicted for several targets, with the largest effects occurring near the "sweet spot" where quench pressures coincide with photospheric levels (Zamyatina et al., 2022, Liu et al., 9 Apr 2026).
• Emission spectra and phase curves: Quenched (enhanced) opacity redistributes the IR photosphere to lower pressures (cooler temperatures), reducing dayside-to-nightside flux contrast in bands dominated by key species (e.g., CH$K_{zz}$8, CO). For example, uniform CH$K_{zz}$9 reduces the day–night contrast at 3.6 and 8 $K_h$0m by up to 30%. In some cases (e.g., at 4.5 $K_h$1m) compensation between CO and H$K_h$2O opacities neutralizes the expected disequilibrium signature (Steinrueck et al., 2018).
• JWST/ARIEL detectability: Simulation studies predict that future instrumentation can distinguish disequilibrium models at $K_h$3 across multiple modes and spectral bands, particularly for those species whose abundance is most amplified by transport (e.g., CH$K_h$4 and NH$K_h$5 in specific bands) (Steinrueck et al., 2018, Liu et al., 9 Apr 2026, Zamyatina et al., 2022).
• Terminator and meridional asymmetries: Day–night and east–west limb differences in tracer and temperature profiles produce distinct absorption signatures, e.g., evening limb absorption exceeding morning limb by $K_h$620% in mass mixing ratio and $K_h$712 ppm spectral depth (Liu et al., 30 Jun 2025, Liu et al., 9 Apr 2026).

5. Comparison to 1D and 2D Models: Necessity of Full 3D Treatment

One-dimensional and two-dimensional chemistry–transport models fail to capture crucial aspects of transport-induced disequilibrium seen in full 3D:

• 1D models parameterize vertical mixing via $K_h$8 and can reproduce the global-mean profile with a fitted $K_h$9, but they miss longitudinal and latitudinal transport, which is typically more rapid and typically dominates quenching in strongly zonally banded atmospheres (Steinrueck et al., 2018, Drummond et al., 2020).
• The new chemical-scale-height method ($P_i$0) enables improved 1D estimates for the quench point where full 3D simulation is infeasible (Bordwell et al., 2018).
• 2D and 1D approaches do not capture effects such as equatorial jets, latitudinal/terminator gradients, upwelling or subsidence zones, or eddy-driven meridional mixing, leading to failure in predicting asymmetric spectral or spatial signatures (Liu et al., 30 Jun 2025, Zamyatina et al., 2022).
• Only by self-consistently integrating chemistry, dynamics, and radiative transfer in 3D can feedbacks such as local heating/cooling, zonal wind modifications, and cloud–chemistry interactions be captured (Steinrueck et al., 2018, Bell et al., 2024, Drummond et al., 2020).

6. Feedbacks: Radiative, Dynamical, and Cloud-Chemistry Interactions

Disequilibrium chemistry has quantitative feedbacks onto atmospheric structure:

• Thermal structure is directly affected: in HD 189733b, quenched CO reduces dayside temperatures by 50–100 K and warms large portions of the nightside by $P_i$150 K; CH$P_i$2-dominated regimes can increase sub-photospheric temperatures by 200–400 K (Steinrueck et al., 2018).
• Altered opacity profiles impact the radiative heating rates, shifting jet strengths, day–night temperature contrasts, and hotspot offsets.
• In WASP-43b, nightside cloud condensation, becoming optically thick above 100 mbar, further cools/isothermalizes the upper atmosphere, suppressing chemical kinetics, and amplifying disequilibrium (and matching observed CH$P_i$3 non-detections) (Bell et al., 2024).
• Cloud microphysics and photochemistry (e.g., UV-driven haze formation) are often coupled to transport; both can enhance or obscure disequilibrium spectral features (Zamyatina et al., 2022, Bell et al., 2024).

7. Broader Implications, Limitations, and Prospects

Accurate retrieval of atmospheric compositions in planetary and exoplanetary atmospheres now requires accounting for 3D transport-induced disequilibrium chemistry. This includes:

• Use of $P_i$4 from 3D GCMs as input to 1D retrievals for temperate sub-Neptunes and mini-Neptunes (Liu et al., 9 Apr 2026).
• Consideration of latitudinal circulation and its quenching effects for Ice Giants, with observed meridional gradients constraining global upwelling/downwelling (Moses et al., 2020).
• In porous or complex flows, explicit modeling of pore-scale lamella folding is necessary for predicting reaction rates and disequilibrium hot spots, requiring detailed upscaling for natural geological or industrial systems (Izumoto et al., 2023).
• Future models must merge photochemistry, cloud microphysics, and 3D dynamics, allowing non-equilibrium abundances and feedbacks to be predicted and compared directly to high-precision observations.

Limitations of current models include incomplete coupling of cloud radiative effects, photochemistry at low pressures, and insufficient spatial resolution for turbulent/eddy mixing. Recommendations from recent work include the extension to fully composition-coupled radiative transfer, explicit eddy diagnostics, and adoption of chemically informed scaling laws for all mixing schemes (Bordwell et al., 2018, Steinrueck et al., 2018, Zamyatina et al., 2022).

In summary, three-dimensional transport-induced disequilibrium chemistry is a dynamically and chemically controlled regime where atmospheric transport timescales outpace local reaction kinetics, producing observable and dynamically significant deviations from chemical equilibrium—a fundamental process shaping the atmospheric structures and spectra of a diverse range of planetary environments (Steinrueck et al., 2018, Liu et al., 9 Apr 2026, Zamyatina et al., 2022, Bell et al., 2024, Liu et al., 30 Jun 2025, Drummond et al., 2020, Mendonça et al., 2018, Bordwell et al., 2018, Moses et al., 2020, Izumoto et al., 2023).