3D Thermal-Driven Ocean Circulation

• 3D thermal-driven ocean circulation is a system of global-scale, stratified flows driven by thermal gradients, creating conveyor-belt structures and multi-scale heat/mass transport.
• It employs innovative diagnostic techniques like Lagrangian Coherent Structures and multi-layer isopycnal frameworks to reveal complex vertical and meridional exchange pathways.
• The study integrates reduced-order models and full 3D simulations to quantify the influences of mixing, rotation, and boundary forcing on deep and surface current dynamics.

Three-dimensional (3D) thermal-driven ocean circulation encompasses the global-scale, stratified, and rotating flows sustained primarily by gradients in heat input, forming a central component of planetary interiors, atmospheric–oceanic climate, and even subsurface seas of icy worlds. Thermal gradients—such as those imposed by insolation, geothermal sources, or phase-change boundaries—drive vertical and meridional overturning, baroclinic gyres, and deep/interior boundary currents. Recent research leverages Lagrangian Coherent Structures (LCS), multi-layer isopycnal frameworks, and reduced-order geostrophic theory to elucidate how 3D geometry, mixing, and boundary forcing interact to produce conveyor-belt–like structures, diverse upwelling/downwelling regimes, and multi-scale heat/mass transport in both terrestrial and extraterrestrial oceans (Bruera et al., 2024, Vanderborght et al., 22 Oct 2025, Bhagtani et al., 2023, Yang et al., 1 Feb 2026, Lobo et al., 2020, Lemasquerier et al., 2023).

1. Governing Frameworks for Thermal-Driven Circulation

The fundamental equations for 3D thermal-driven flows are the rotating Boussinesq or anelastic Navier–Stokes equations with buoyancy (thermal) forcing:

$$\frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla)\mathbf{u} + 2\boldsymbol{\Omega} \times \mathbf{u} = -\frac{1}{\rho_0}\nabla p' + b \mathbf{\hat{z}} + \mathcal{V}(\mathbf{u})$$

$$\frac{\partial b}{\partial t} + \mathbf{u} \cdot \nabla b = \kappa \nabla^2 b + S_b$$

$\nabla \cdot \mathbf{u} = 0$

where $\mathbf{u}$ is the velocity, $b = g(1-\rho/\rho_0)$ is the buoyancy, $p'$ is pressure perturbation, $\Omega$ is planetary rotation, $\mathcal{V}$ denotes viscous/stress terms, $\kappa$ is the diffusivity, and $S_b$ incorporates sources such as surface heat, melting, or bottom heating (Bhagtani et al., 2023, Bruera et al., 2024, Yang et al., 1 Feb 2026). Stratification arises from vertical and horizontal density gradients forced by temperature (and, on Earth or icy moons, also salinity or pressure) variations.

Boundary conditions vary: free surface with wind and buoyancy fluxes (Earth), fixed-temperature or fixed-flux (ocean world interiors), and moving domain boundaries where phase transitions (e.g., lava, ice) are encountered (Lobo et al., 2020, Yang et al., 1 Feb 2026, Lemasquerier et al., 2023). Dynamics are further shaped by domain geometry, mixing rates, rotation (Coriolis effects), and planetary parameters such as stratification and Rayleigh numbers.

2. Conveyor-Belt Geometry, Mixing, and Lagrangian Transport

Thermal gradients and 3D geometry generate complex conveyor-like pathways. Lagrangian Coherent Structures (LCS) provide a rigorous diagnostic for the architecture of 3D transport: ridges of the finite-time Lyapunov exponent (FTLE) or sharp gradients in arc-length–based Lagrangian descriptors $M$ reveal intersecting surfaces that act as material barriers for vertical and horizontal exchange (Bruera et al., 2024). Repelling LCS (backward-time) separate regions whose parcels diverge, while attracting LCS (forward-time) delineate convergence surfaces.

Key findings in the North Atlantic Meridional Overturning Circulation (AMOC):

• LCS “curtains” or “ridges” vertically link deep and surface layers, forming 3D conveyor belts.
• Rapid upwelling (∼80 days from 1.3 km) emerges at shelf regions such as Flemish Cap due to sharp LCS ridges.
• Abyssal-to-surface upwelling pathways, taking ∼840 days from 3.5 km, are organized along mesoscale eddy boundaries in regions like the Irminger Sea.
• Downwelling times for surface waters into the AMOC branches at several locations are O(1000) days (Bruera et al., 2024).

Clustered parcel trajectories within these LCS-defined corridors maintain coherence over years, confirming the power of this geometric approach in tracing and quantifying 3D vertical/meridional interchange.

3. Scaling Laws, Reduced-Order Models, and Global Structure

Reduced-dimensional and analytic models clarify how thermal gradients, mixing, and rotation set the spatial/temporal scales of overturning:

• Below the Ekman layer, geostrophic balance ($f\mathbf{k} \times \mathbf{u}_g = -\nabla p/\rho_0$) and thermal wind ($$f \partial \mathbf{u}/\partial z = -g/\rho_0 \nabla_h \rho$$) govern vertical shears and boundary currents (Vanderborght et al., 22 Oct 2025).
• Overturning strength $\Psi$ scales as $\Psi \propto \kappa_b^{2/3} \Delta T^{1/3}$ in diffusive regimes; upwelling in re-entrant channels is set by adiabatic wind-induced residuals or Gent–McWilliams-like eddy parameterizations.
• Cross-equatorial flows and boundary-interior exchanges are diagnosed by contrasts in boundary temperature/density and boundary-intensified mixing (Vanderborght et al., 22 Oct 2025, Jonathan, 2022).
• Both analytic and full 3D GCM (MITgcm) simulations confirm these scalings within 10–20% for the streamfunction, up/downwelling, and stratification (Vanderborght et al., 22 Oct 2025).

In global Earth oceans, baroclinic vorticity input from surface buoyancy gradients (∂Q/∂y) adds quasi-linearly to the classical wind-driven Sverdrup transport on decadal timescales (anomaly in gyre strength ΔΨ ~ 0.15 Sv per W m⁻² gradient in surface heat flux), but interacts nonlinearly via interior stratification and eddy feedbacks on longer timescales (Bhagtani et al., 2023).

4. Planetary Contexts: Icy and Magma Ocean Worlds

Thermal-driven ocean circulation is pervasive across diverse planetary bodies:

Enceladus and Icy Worlds

• Stratified, meridionally overturning cells are forced by meridional variations in ice-ocean interface thermal/mass flux, salinity, and seafloor heating patterns.
• Isopycnal-layer models and MITgcm simulations demonstrate that spatially separated regions of melting (pole) and freezing (equator) enforce a pole-to-equator overturning, with strengths 0.1–0.5 Sv and shallow stratification near the surface (Lobo et al., 2020, Zeng et al., 2023).
• Positive feedbacks between cross-equatorial ocean heat transport and basal melt can induce bistability and hemisphere-scale symmetry breaking in the ice shell (Kang et al., 2022).
• Convective translation of seafloor tidal heating through the ocean can amplify or suppress pole-to-equator ice thickness contrasts, with the efficiency controlled by the regime (thermal winds versus zonal jets) (Lemasquerier et al., 2023).

Magma Oceans on Lava Worlds

• 3D thermal-only–forced magma oceans on tidally locked exoplanets exhibit rotation-dominated, shallow gyres and overturning, with heat transport ($10^3$–$10^4$ W m⁻²) far smaller than incident stellar flux ($10^6$ W m⁻²).
• Coriolis effects induce western intensification, hemispheric gyres, and displacement of the deepest thermocline away from substellar points; vertical diffusivity strongly regulates mean depth and overturning strength (Yang et al., 1 Feb 2026).

5. Observational and Computational Diagnostics

Analysis and prediction of 3D thermal-driven circulation rely on a spectrum of methodologies:

• Lagrangian Coherent Structures (LCS): FTLE and $M$-based diagnostics quantify transport barriers and upwelling/downwelling corridors (Bruera et al., 2024).
• Boundary-tracing/Decomposition: AMOC and ACC streamfunctions can be robustly reconstructed from densities along sloping western/eastern (or northern/southern) boundaries plus Ekman and depth-independent terms, reducing 3D overturning characterization to a handful of boundary profiles (Jonathan, 2022).
• Data-driven 3D Models: Transformer-based models such as ORCA-DL, trained on OGCM outputs, replicate stratification, overturning, and heat content anomalies, matching or exceeding predictive skill for ENSO and decadal variability compared to classical OGCMs (Guo et al., 2024).

6. Implications, Feedbacks, and Theoretical Synthesis

Three-dimensional, thermally forced ocean circulation regulates the storage and exchange of heat, carbon, and biogeochemical tracers, and underpins the climatic roles of the AMOC, gyres, and planetary energy budgets. Crucial aspects include:

• The geometry of surface–deep exchange pathways, their coherence, and their mixing/lateral coupling control the timescales and regions by which thermal properties (and by extension, dissolved gases and nutrients) are redistributed.
• Boundary densities and stratification set not only overturning rates but also the sensitivity of circulation to external forcing, facilitating regime shifts, feedback loops, and multi-stability (especially prominent in icy satellite oceans).
• On Earth, combining wind-driven (Sverdrup) and buoyancy-driven (baroclinic) vorticity inputs yields a more complete predictive theory of gyre and overturning variability, especially as interior stratification and eddy transport are perturbed under climate change (Bhagtani et al., 2023).

Quantifying, diagnosing, and predicting 3D thermal-driven circulation thus demands integrated perspectives: geometric Lagrangian diagnostics, analytic and reduced-order models, high-resolution numerics, and data-driven emulation—each validated by boundary measurements, stratification structure, and transport observation (Bruera et al., 2024, Vanderborght et al., 22 Oct 2025, Bhagtani et al., 2023, Yang et al., 1 Feb 2026, Guo et al., 2024, Lobo et al., 2020, Zeng et al., 2023, Kang et al., 2022, Lemasquerier et al., 2023, Jonathan, 2022).