pith. sign in

arxiv: 2604.14342 · v1 · submitted 2026-04-15 · 🌌 astro-ph.EP · astro-ph.IM· physics.ao-ph

A Fast and Physically Grounded Ocean Model for GCMs: The Dynamical Slab Ocean Model of the Generic-PCM (rev. 3423)

Siddharth Bhatnagar , Francis Codron , Ehouarn Millour , Emeline Bolmont , Maura Brunetti , J\'er\^ome Kasparian , Martin Turbet , Guillaume Chaverot This is my paper

Pith reviewed 2026-05-10 11:47 UTC · model grok-4.3

classification 🌌 astro-ph.EP astro-ph.IMphysics.ao-ph
keywords slab ocean modelocean heat transportGCMexoplanet climateSverdrup balanceGent-McWilliams parameterizationsea ice albedo
0
0 comments X

The pith

A dynamical slab ocean model adds realistic heat transport to fast GCMs at almost no extra cost.

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

The paper develops a slab ocean model that includes ocean heat transport through a Sverdrup balance for wind-driven currents and the Gent-McWilliams parameterization for mesoscale eddies. This addresses the gap in exoplanet GCMs where full ocean dynamics are too slow but static slabs omit key heat redistribution effects. Aquaplanet tests show cooler tropics, less sea ice, and a new double precipitation band, while Earth runs match observed global temperature within 1°C and albedo within 0.01. The approach stays computationally light due to parallelization, enabling long runs needed for broad parameter sweeps. If the approximations hold, they let modelers explore how ocean transport shapes climate and observables on terrestrial planets without full 3D ocean costs.

Core claim

The authors build a single-layer dynamical slab ocean that calculates horizontal heat transport from wind stress via Sverdrup balance and from eddy mixing via the first use of Gent-McWilliams in a slab setting, plus a spectrally dependent sea-ice albedo. In aquaplanet simulations this produces ocean heat transport profiles that agree at first order with fully coupled atmosphere-ocean GCMs, cooler tropical sea-surface temperatures, reduced sea ice, and an Ekman-driven double-banded equatorial precipitation pattern. Applied to present-day Earth the model yields a global mean surface temperature of 13°C, planetary albedo of 0.32, and sea-ice extent with smaller seasonal biases than static-slab

What carries the argument

Sverdrup balance plus Gent-McWilliams parameterization inside a single-layer slab, which computes large-scale ocean heat transport directly from surface wind stress and eddy effects without vertical structure.

If this is right

  • Enabling ocean heat transport cools tropical sea-surface temperatures and reduces overall sea ice compared with static slabs.
  • Ekman upwelling produces a double-banded equatorial precipitation pattern absent in no-OHT runs.
  • Earth simulations match observed global temperature to within 1°C and albedo to within 0.01.
  • Sea-ice extent seasonal cycle improves markedly relative to simulations that omit ocean heat transport.
  • The added physics costs almost nothing extra in wall-clock time, allowing the same number of model years as static-slab runs.

Where Pith is reading between the lines

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

  • The model could be used to test how changes in rotation rate or atmospheric mass alter ocean-driven climate patterns on exoplanets.
  • Its sea-ice albedo treatment might reveal sensitivities in paleoclimate runs where ice extent strongly feedbacks on temperature.
  • If the first-order OHT agreement persists across a wider range of stellar fluxes, the framework becomes a practical default for terrestrial exoplanet ensembles.

Load-bearing premise

A single-layer slab with Sverdrup balance and Gent-McWilliams parameterization is enough to capture the dominant large-scale ocean heat transport.

What would settle it

Run the model on an aquaplanet and compare its zonally averaged ocean heat transport profile directly against a high-resolution fully coupled GCM using identical atmosphere and forcing; first-order mismatch would show the single-layer approximations fail.

Figures

Figures reproduced from arXiv: 2604.14342 by Ehouarn Millour, Emeline Bolmont, Francis Codron, Guillaume Chaverot, J\'er\^ome Kasparian, Martin Turbet, Maura Brunetti, Siddharth Bhatnagar.

Figure 1
Figure 1. Figure 1: (a) Bare sea ice albedo as a function of ice thickness in the Visible (VIS; green) and Near￾Infrared (NIR; red) bands in our study, with curves fit to Antarctic observations from Brandt et al. (2005, filled symbols: measurements; open: interpolated). The spectrally-independent albedo used in Charnay et al. (2013) is shown for reference (dashed green). (b) Broadband spectral distribution of snow and ice alb… view at source ↗
Figure 2
Figure 2. Figure 2: Cross-sectional representation of the ocean along latitude and depth. The top section [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (a) Zonal mean sea surface temperature (SST, 10-year average) of the aquaplanet for Case 1 (no oceanic heat transport, OHT) and Case 7 (all OHT on) simulations. (c) Difference in SST (ON – OFF), showing the effect of OHT on the same. (b) Zonal mean precipitation (rain + snow, 10-year average) for the same simulations. (d) Difference in precipitation (ON – OFF) resulting in the pronounced peak in equatorial… view at source ↗
Figure 4
Figure 4. Figure 4: Zonal-mean sea surface temperature (SST) differences of the aquaplanet between selected [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Contribution of each ocean heat transport (OHT) term towards the northward meridional [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Example of the causal chain of the coupled ocean-atmosphere system: Impact of ocean [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Total northward meridional ocean heat transport (in Petawatts) for modern Earth, as [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: (a) Seasonal evolution of extrapolar sea surface temperatures (SST; 60°S–60°N), and, (b) global sea ice coverage for NOAA/NSIDC observations (black), for the model with ocean heat transport enabled (red), and disabled (blue). Shaded regions represent the 2σ inter-annual variabil￾ity (30–40 years). Vertical dashed lines indicate the timing of the March and September equinoxes. The seasonal evolution of extr… view at source ↗
Figure 9
Figure 9. Figure 9: Gent-McWilliams (GM) restratification relative to convective adjustment: Zonally aver [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Total and decomposed ocean heat transport (OHT) for [PITH_FULL_IMAGE:figures/full_fig_p023_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Ocean heat transport (OHT) suppresses hemispheric climate asymmetries caused by [PITH_FULL_IMAGE:figures/full_fig_p024_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: The stabilising influence of the Sverdrup balance on the climate: [PITH_FULL_IMAGE:figures/full_fig_p025_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Zonally averaged surface air temperature for modern Earth simulations under two [PITH_FULL_IMAGE:figures/full_fig_p027_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Annually averaged sea surface temperature (SST) maps for the zero-obliquity aquaplanet [PITH_FULL_IMAGE:figures/full_fig_p034_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Annually averaged atmospheric profiles for the zero-obliquity aquaplanet simulations. [PITH_FULL_IMAGE:figures/full_fig_p035_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Seasonal evolution of global sea ice coverage for NSIDC observations (black), the model [PITH_FULL_IMAGE:figures/full_fig_p036_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Decadal-mean net surface wind vectors for simulations [PITH_FULL_IMAGE:figures/full_fig_p037_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Decadal-mean zonally averaged plots for (a) sea surface temperature (SST) for OHT￾off (blue) and OHT-on (red) configurations in the present study (solid lines) and Codron (2012) (dashed lines), and (b) precipitation. Differences between the two implementations reflect changes in the atmospheric model and resolution between the studies. The large difference in equatorial temperatures in the OHT-off cases (… view at source ↗
Figure 19
Figure 19. Figure 19: Influence of Gent–McWilliams (GM) transport on deep-ocean thermal structure: [PITH_FULL_IMAGE:figures/full_fig_p039_19.png] view at source ↗
read the original abstract

Ocean dynamics are often sidelined in exoplanet climate studies due to the high computational cost of fully coupled atmosphere-ocean general circulation models (GCMs). However, ocean heat transport (OHT) can play a critical role in shaping the climate and observables of terrestrial planets. As a compromise, most exoplanet GCMs rely on slab ocean models without OHT. Here, we present an improved compromise - a fast and physically grounded dynamical slab ocean model, implemented in the Generic Planetary Climate Model (Generic-PCM). The model extends previous frameworks by incorporating a Sverdrup balance formulation for wind-driven Ekman transport, the first application of the Gent-McWilliams parameterisation of mesoscale eddies in a slab ocean model, and a spectrally and thickness-dependent treatment of sea ice and snow albedo. In aquaplanet simulations, enabling OHT produces substantial changes in both surface climate and atmospheric circulation, including cooler tropical sea surface temperatures, reduced sea ice, and the emergence of a double-banded equatorial precipitation pattern driven by Ekman-induced upwelling. The resulting OHT profiles show first-order agreement with fully coupled atmosphere-ocean GCMs. Applied to modern Earth, the model reproduces key large-scale climate properties, including a global mean surface temperature of 13{\deg}C (within 1{\deg}C of observations), planetary albedo of 0.32 (within 0.01), and sea ice extent with significantly reduced seasonal biases relative to simulations without OHT. Due to model parallelisation, these improvements are achieved at almost no additional computation cost compared to OHT-disabled simulations run over the same number of model years. This enables long integrations, making the model particularly well suited for exoplanet and paleoclimate studies where broad parameter exploration is essential.

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 / 3 minor

Summary. The paper introduces a dynamical slab ocean model for the Generic-PCM that incorporates Sverdrup balance for wind-driven Ekman transport, the Gent-McWilliams parameterization for mesoscale eddies, and spectrally/thickness-dependent sea ice and snow albedo. In aquaplanet simulations, enabling OHT alters tropical SSTs, reduces sea ice, and produces a double-banded equatorial precipitation pattern, with resulting OHT profiles showing first-order agreement with full atmosphere-ocean GCMs. Applied to modern Earth, the model reproduces a global mean surface temperature of 13°C (within 1°C of observations), planetary albedo of 0.32 (within 0.01), and sea ice extent with reduced seasonal biases relative to no-OHT runs, all at nearly zero additional computational cost due to parallelization.

Significance. If the central claims hold, the model offers a computationally efficient and physically motivated way to include ocean heat transport in exoplanet GCMs, enabling long integrations and broad parameter sweeps that are impractical with full 3D ocean models. The parallelization and lack of added cost are clear strengths for paleoclimate and exoplanet applications.

major comments (2)
  1. [Abstract and Earth validation section] Abstract and results section on Earth validation: The reported agreement (GMST of 13°C within 1°C, albedo 0.32 within 0.01, reduced sea-ice bias) is presented without quantitative error bars, RMS differences, or a description of the parameter selection/tuning procedure. This makes it difficult to evaluate whether the match constitutes an independent test or a post-hoc fit to observed global temperature and albedo.
  2. [Model formulation section] Model formulation section: The single-layer Sverdrup + Gent-McWilliams formulation omits vertical density structure and thermohaline/meridional overturning contributions to OHT. Since these components can dominate (~0.5–1 PW) in Earth-like regimes, the first-order agreement demonstrated in aquaplanet cases may not generalize, weakening the claim that the slab is 'physically grounded' for arbitrary exoplanet climates.
minor comments (3)
  1. [Abstract] The abstract states 'significantly reduced seasonal biases' for sea ice but provides no quantitative measure (e.g., seasonal amplitude or RMS error) to support this.
  2. [Figures and results] Figure captions and results text should explicitly label which curves correspond to the dynamical slab OHT versus reference GCMs and observations for direct comparison.
  3. [Methods] The manuscript would benefit from a brief statement on the horizontal resolution and integration length used in the aquaplanet and Earth simulations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their insightful and constructive comments. We address each major comment below and have revised the manuscript accordingly to improve clarity and transparency.

read point-by-point responses

  1. Referee: [Abstract and Earth validation section] Abstract and results section on Earth validation: The reported agreement (GMST of 13°C within 1°C, albedo 0.32 within 0.01, reduced sea-ice bias) is presented without quantitative error bars, RMS differences, or a description of the parameter selection/tuning procedure. This makes it difficult to evaluate whether the match constitutes an independent test or a post-hoc fit to observed global temperature and albedo.

    Authors: We agree that additional quantitative metrics would strengthen the Earth validation section. In the revised manuscript, we now report RMS differences for GMST, planetary albedo, and sea ice extent relative to observations. We also clarify the parameter selection: the Gent-McWilliams coefficient and other constants were taken from standard literature values without iterative tuning or optimization to match global temperature or albedo. The agreement is therefore presented as an a posteriori validation. revision: yes

  2. Referee: [Model formulation section] Model formulation section: The single-layer Sverdrup + Gent-McWilliams formulation omits vertical density structure and thermohaline/meridional overturning contributions to OHT. Since these components can dominate (~0.5–1 PW) in Earth-like regimes, the first-order agreement demonstrated in aquaplanet cases may not generalize, weakening the claim that the slab is 'physically grounded' for arbitrary exoplanet climates.

    Authors: We acknowledge that the single-layer formulation necessarily omits thermohaline circulation and full vertical density structure. This is a deliberate simplification for computational efficiency in exoplanet applications. We have added an expanded discussion section that explicitly states these limitations, notes their importance in Earth-like regimes, and specifies the conditions (e.g., wind-driven dominated transport in aquaplanets) under which the model provides a physically grounded first-order approximation. We retain the 'physically grounded' phrasing only in reference to the included parameterizations while being transparent about the omissions. revision: partial

Circularity Check

0 steps flagged

No significant circularity in dynamical slab ocean derivation

full rationale

The paper derives OHT via explicit dynamical equations (Sverdrup balance for Ekman transport plus Gent-McWilliams eddy parameterization) applied to a single-layer slab, then compares the resulting profiles to full 3D GCMs and Earth observations. This chain is self-contained: the transport is an output of the stated physics rather than a fitted input, self-defined quantity, or result imported solely via author self-citation. Earth reproduction of GMST, albedo, and sea-ice metrics is presented as validation of the independent computation, not as a construction that forces the outcome. No load-bearing self-citation, ansatz smuggling, or renaming of known results is exhibited in the provided text. Score 0 is the appropriate default for a model whose central claim rests on standard dynamical closures without reduction to its own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of the Sverdrup balance and Gent-McWilliams closure when applied to a single-layer ocean; these are standard but approximate closures whose accuracy in a slab setting is not independently verified here.

axioms (2)
  • domain assumption Sverdrup balance holds for wind-driven transport in the slab layer
    Invoked to compute Ekman transport from wind stress without solving full momentum equations.
  • domain assumption Gent-McWilliams parameterization can be applied directly to a slab ocean without vertical resolution
    Used to represent mesoscale eddy heat transport; standard in 3D models but its slab adaptation is new.

pith-pipeline@v0.9.0 · 5685 in / 1417 out tokens · 25929 ms · 2026-05-10T11:47:51.072214+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

13 extracted references · 13 canonical work pages

  1. [1]

    Bhatnagar: Slab ocean model documentation,https://lmdz-forge.lmd.jussieu.fr/mediawiki/ Planets/index.php/Slab_ocean_model, last access: 8 August 2025,

    work page 2025

  2. [2]

    librosa/librosa: 0.6.3,

    Bhatnagar, S., Codron, F., Millour, E., Bolmont, E., Brunetti, M., Kasparian, J., Turbet, M., and Chaverot, G.: Dataset for "A Fast and Physically Grounded Ocean Model for GCMs: The Dynamical Slab Ocean Model of the Generic-PCM (rev. 3423)", https://doi.org/10.5281/zenodo. 18771594,

    work page doi:10.5281/zenodo

  3. [3]

    Influence of the internal structure, Astronomy & Astrophysics, https://doi.org/10.1051/0004-6361/202554529,

    Bolmont, E., Sastre, M., Revol, A., Kervazo, M., and Tobie, G.: Estimation of the tidal heating in the TRAPPIST-1 planets. Influence of the internal structure, Astronomy & Astrophysics, https://doi.org/10.1051/0004-6361/202554529,

    work page doi:10.1051/0004-6361/202554529

  4. [4]

    N., Savoie, M., and Windnagel, A

    Fetterer, F., Knowles, K., Meier, W. N., Savoie, M., and Windnagel, A. K.: Sea Ice Index, Version 3, [Data set], https://doi.org/10.7265/N5K072F8,

    work page doi:10.7265/n5k072f8

  5. [5]

    S., Seager, S., Mennesson, B., et al

    Gaudi, B. S., Seager, S., Mennesson, B., Kiessling, A., Warfield, K., Cahoy, K., Clarke, J. T., Domagal-Goldman, S., Feinberg, L., Guyon, O., et al.: The Habitable Exoplanet Observatory (HabEx) mission concept study final report, arXiv preprint arXiv:2001.06683,

    work page arXiv 2001

  6. [6]

    Generic-PCM-Team: Overview of the Generic PCM,https://lmdz-forge.lmd.jussieu.fr/ mediawiki/Planets/index.php/Overview_of_the_Generic_PCM, last access: 8 August 2025,

    work page 2025

  7. [7]

    Generic-PCM-Team: The Generic-PCM model,http://svn.lmd.jussieu.fr/Planeto/trunk, last access: 8 August 2025,

    work page 2025

  8. [8]

    and Peale, S.: Spin-orbit coupling in the solar system, Astronomical Journal, Vol

    Goldreich, P. and Peale, S.: Spin-orbit coupling in the solar system, Astronomical Journal, Vol. 71, p. 425 (1966), 71, 425,

    work page 1966

  9. [9]

    Kalnay, E., Kanamitsu, M., Kistler, R., Collins, W., Deaven, D., Gandin, L., Iredell, M., Saha, S., White, G., Woollen, J., Zhu, Y., Leetmaa, A., Reynolds, B., Chelliah, M., Ebisuzaki, W., Higgins, W., Janowiak, J., Mo, K. C., Ropelewski, C., Wang, J., Jenne, R., and Joseph, D.: The NCEP/NCAR 40-Year Reanalysis Project., Bulletin of the American Meteorolo...

    work page doi:10.1175/1520-0477(1996)077 1996

  10. [10]

    Marshall, J., Donohoe, A., Ferreira, D., and McGee, D.: The ocean’s role in setting the mean position of the Inter-Tropical Convergence Zone, Climate Dynamics, 42, 1967–1979,

    work page 1967

  11. [11]

    Rose, B. E. J.: Stable ‘Waterbelt’ climates controlled by tropical ocean heat transport: A nonlinear coupled climate mechanism of relevance to Snowball Earth, Journal of Geophysical Research: Atmospheres, 120, 1404–1423, https://doi.org/https://doi.org/10.1002/2014JD022659,

    work page doi:10.1002/2014jd022659

  12. [12]

    46 Turbet, M. and Forget, F.: 3-D Global modelling of the early martian climate under a dense CO2+ H2 atmosphere and for a wide range of surface water inventories, arXiv preprint arXiv:2103.10301,

    work page arXiv

  13. [13]

    and Abbot, D.: Sea-ice dynamics strongly promote Snowball Earth initiation and desta- bilize tropical sea-ice margins, Climate of the Past, 8, 2079–2092,

    Voigt, A. and Abbot, D.: Sea-ice dynamics strongly promote Snowball Earth initiation and desta- bilize tropical sea-ice margins, Climate of the Past, 8, 2079–2092,

    work page 2079