Nonlocally coupled moisture model for convective self-aggregation

Hirofumi Tomita; Seiya Nishizawa; Shin-ichiro Shima; Tomoro Yanase

arxiv: 2404.04146 · v2 · pith:U5ETY6BAnew · submitted 2024-04-05 · ⚛️ physics.ao-ph · nlin.AO· nlin.PS· physics.flu-dyn

Nonlocally coupled moisture model for convective self-aggregation

Tomoro Yanase , Shin-ichiro Shima , Seiya Nishizawa , Hirofumi Tomita This is my paper

Pith reviewed 2026-05-24 02:32 UTC · model grok-4.3

classification ⚛️ physics.ao-ph nlin.AOnlin.PSphysics.flu-dyn

keywords convective self-aggregationmoisture modelnonlocal couplingbistabilityweak temperature gradientatmospheric columnsdiffusion stabilizationmoisture mode

0 comments

The pith

Nonlocal coupling of moisture columns induces bistability between dry and moist equilibria that triggers convective self-aggregation.

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

The paper introduces a reduced mathematical model in which the atmosphere is represented solely by the vertically integrated water vapor in each column under the weak temperature gradient approximation. Analysis of the resulting nonlinear system demonstrates that nonlocal coupling between columns generates two stable fixed points, one dry and one moist. Self-aggregation occurs once finite-amplitude disturbances allow the destabilizing effect of this coupling to exceed the stabilizing effect of diffusion. The dominant wavelength of the resulting pattern is set by the spatial range of the coupling: global coupling favors the longest possible scale while finite-range coupling selects scales comparable to its filter length. The work therefore identifies the coupling-diffusion balance as the controlling mechanism for both the onset and the spatial structure of aggregated convection.

Core claim

In the nonlocally coupled moisture model, nonlocal coupling between atmospheric columns induces bistability, leading to dry and moist equilibria. This reflects the circulation effects driven by horizontal differential heating due to convection and radiation. The bistable self-aggregated state realizes when destabilization by nonlocal coupling, triggered by finite-amplitude disturbances in the uniform state, overcomes stabilization by diffusion. For globally coupled systems where all columns are equally coupled, perturbations with the maximum wavelength exhibit the highest growth rate, resulting in a solution with an infinitely long wavelength understood as the dynamical system's heteroclinic

What carries the argument

Nonlocal coupling operator applied to column water vapor that generates bistable dry and moist equilibria.

If this is right

Finite-amplitude perturbations, not infinitesimal noise, are sufficient to initiate aggregation once coupling strength exceeds diffusion.
Globally coupled columns evolve toward domain-scale aggregation corresponding to heteroclinic connections between dry and moist states.
Finite-range coupling selects aggregation wavelengths near the characteristic coupling length rather than the longest available scale.
The uniform moist state remains stable against small disturbances but can be destabilized by larger ones when nonlocal effects dominate.
The same coupling-diffusion competition that sets the aggregation scale also governs the transition between aggregated and uniform regimes.

Where Pith is reading between the lines

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

The reduction to a single moisture field may allow similar bistable behavior to be identified in other tropical convective regimes where temperature gradients are weak.
Parameterizations that incorporate nonlocal moisture coupling could reproduce self-aggregation without resolving full three-dimensional dynamics.
The preference for finite wavelengths under nonlocal coupling offers a possible explanation for the observed mesoscale size of aggregated cloud clusters.
Extending the model to include time-dependent coupling strength could predict how aggregation responds to changing large-scale conditions.

Load-bearing premise

The weak temperature gradient approximation holds so that the state of each atmospheric column is fully captured by its vertically integrated water vapor alone.

What would settle it

A simulation or observation in which self-aggregation develops while column water vapor remains spatially uniform would contradict the claim that bistability in moisture is required.

Figures

Figures reproduced from arXiv: 2404.04146 by Hirofumi Tomita, Seiya Nishizawa, Shin-ichiro Shima, Tomoro Yanase.

**Figure 1.** Figure 1: Schematic of the moisture model. (a) The budget equation for moisture leads to the prognostic equation for vertically integrated water vapor content that includes surface evaporation (Evp.), surface precipitation (Prec.), vertical advection (Vadv.), horizontal advection (Hadv.), and horizontal eddy diffusion (Hdiff.). (b) The budget equation for heat results in the diagnostic equation for vertically integr… view at source ↗

**Figure 2.** Figure 2: Dynamical features of the reaction term. (a) Potential function 𝑉(𝑞") for different 𝑞%". Filled and open circles represent stable and unstable fixed points, respectively. (b) Bifurcation diagram of 𝑞" against 𝑞%". Solid and dashed lines represent stable and unstable fixed points, respectively. The dotted line represents 𝑞" = 𝑞%". The water balance 𝑃( = 𝐸 is assumed. Green and brown arrows indicate the dire… view at source ↗

**Figure 6.** Figure 6: An example of the trajectory and horizontal distribution of periodic solutions based on the dynamical system described by a set of Eq. (16). (a) Trajectory in a 3D phase space. (b)–(d) Horizontal variation of the state variables 𝑢, 𝑞", and 𝑤. (e)–(g) Trajectories of pairs of state variables projected onto a 2D phase plane. (b) (c) (d) (e) (f) (g) x u qv w u u qv qv w w (a) qv u 0.00 w x , u x , qv x , w u … view at source ↗

read the original abstract

Clouds play a central role in climate physics by interacting with precipitation, radiation, and circulation. Despite being a fundamental issue in convective organization, the self-aggregation of clouds lacks a theoretical explanation due to its complexity. In this study, we introduce an idealized mathematical model where the system's state is represented solely by the vertically integrated water vapor content of atmospheric columns under the weak temperature gradient approximation. By analyzing the nonlinear dynamics of this simplified system, we mathematically elucidate the mechanisms that determine the onset of self-aggregation and the spatial scale of the self-aggregated state. Nonlocal coupling between atmospheric columns induces bistability, leading to dry and moist equilibria. This reflects the circulation effects driven by horizontal differential heating due to convection and radiation. The bistable self-aggregated state realizes when destabilization by nonlocal coupling, triggered by finite-amplitude disturbances in the uniform state, overcomes stabilization by diffusion. For globally coupled systems where all columns are equally coupled, perturbations with the maximum wavelength exhibit the highest growth rate. This results in a solution with an infinitely long wavelength, understood as the dynamical system's heteroclinic trajectories describing the steady state's spatial evolution. Conversely, for nonlocally coupled systems with finite filter lengths, perturbations with wavelengths close to the characteristic length of the coupling are preferred. The results reveal that the balance between nonlocal coupling and diffusion is essential for understanding convective self-aggregation. Moreover, this study suggests a physical similarity between convective self-aggregation and moisture mode.

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.

Desk Editor's Note private letter to a colleague

The paper gives a clean moisture-only model under WTG that produces bistability and wavelength selection from nonlocal coupling, but the WTG closure itself looks untested once aggregation creates large contrasts.

read the letter

The main contribution is a reduced model that keeps only column-integrated moisture, adds explicit nonlocal coupling, and works out the conditions for bistability between dry and moist equilibria plus the preferred scale set by the filter length. That combination of bistability and explicit wavelength selection from the coupling-diffusion balance is not standard in the earlier literature on this setup. The analysis also ties the finite-amplitude trigger to the point where nonlocal effects overcome diffusion, which is a useful way to frame the onset question. The math is laid out directly from the stated approximations, and the global-coupling limit recovering an infinite-wavelength heteroclinic solution is a nice limiting case. The stress-test concern lands. Once the model produces strong horizontal moisture gradients, the WTG diagnostic relation used to close the moisture equation requires temperature perturbations to stay small, yet the paper does not report any a-posteriori check on the size of those perturbations in the aggregated solutions. If the gradients become order-one, the nonlocal operator loses its justification. That is the central modeling choice, so the lack of verification is a real gap rather than a minor detail. The free parameters (filter length and diffusion coefficient) are stated openly, and nothing in the construction looks circular. This is aimed at people who build and analyze idealized tropical-convection models. A reader already working on moisture-mode or self-aggregation theory will find the derivations worth following even if they ultimately disagree with the WTG step. It is worth sending to referees because the construction is explicit and the questions it raises are sharp enough to repay detailed review.

Referee Report

2 major / 2 minor

Summary. The paper presents an idealized single-variable model for convective self-aggregation in which the state is reduced to vertically integrated column moisture under the weak-temperature-gradient (WTG) approximation. Nonlocal coupling between columns is shown to produce bistability between dry and moist equilibria; the aggregated state emerges when nonlocal destabilization overcomes diffusive stabilization. For globally coupled systems the longest-wavelength mode grows fastest, yielding an infinite-wavelength heteroclinic solution, while finite-filter nonlocal coupling selects wavelengths near the filter scale. The work claims a mathematical explanation for the onset and scale selection of self-aggregation and a similarity to moisture-mode dynamics.

Significance. If the derivations and closure remain valid, the model supplies a transparent, low-dimensional dynamical system that isolates the competing roles of nonlocal radiative-convective coupling and diffusion in producing bistability and scale selection. This offers a useful theoretical benchmark for more comprehensive CRM or GCM studies of aggregation and may clarify the link between aggregation and moisture-mode instability.

major comments (2)

[Model formulation and WTG closure] The central claim that nonlocal coupling induces bistability rests on the WTG closure that reduces the system to a single moisture field. No a-posteriori diagnostic is presented showing that diagnosed temperature perturbations remain small once large horizontal moisture contrasts develop in the aggregated equilibria; if those perturbations become O(1), the diagnostic relation used to close the moisture equation and to derive the effective nonlocal operator loses its justification.
[Linear stability analysis] The linear stability analysis that concludes 'perturbations with the maximum wavelength exhibit the highest growth rate' for the globally coupled case (leading to the infinite-wavelength heteroclinic solution) is not accompanied by an explicit dispersion relation or growth-rate formula. Without the explicit eigenvalue problem it is impossible to verify that the claimed wavelength selection follows rigorously rather than from post-hoc inspection of numerical solutions.

minor comments (2)

[Abstract and §2] The abstract states that the model is 'represented solely by the vertically integrated water vapor content' but does not list the explicit governing PDE or the precise form of the nonlocal filter; these should appear in the main text with numbered equations.
[Throughout] Notation for the nonlocal filter length and diffusion coefficient should be introduced once and used consistently; the current text refers to them only descriptively.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their detailed review and constructive comments on our manuscript. We address each major comment below and have prepared revisions to the manuscript accordingly.

read point-by-point responses

Referee: The central claim that nonlocal coupling induces bistability rests on the WTG closure that reduces the system to a single moisture field. No a-posteriori diagnostic is presented showing that diagnosed temperature perturbations remain small once large horizontal moisture contrasts develop in the aggregated equilibria; if those perturbations become O(1), the diagnostic relation used to close the moisture equation and to derive the effective nonlocal operator loses its justification.

Authors: We agree that validating the WTG approximation in the presence of large moisture contrasts is important for the model's credibility. Although the WTG is an approximation commonly employed in tropical dynamics studies, we will add in the revised manuscript an a-posteriori diagnostic of the temperature field in the steady aggregated states to show that temperature perturbations remain small, thereby justifying the closure. revision: yes
Referee: The linear stability analysis that concludes 'perturbations with the maximum wavelength exhibit the highest growth rate' for the globally coupled case (leading to the infinite-wavelength heteroclinic solution) is not accompanied by an explicit dispersion relation or growth-rate formula. Without the explicit eigenvalue problem it is impossible to verify that the claimed wavelength selection follows rigorously rather than from post-hoc inspection of numerical solutions.

Authors: We appreciate this observation. The linear analysis in the manuscript is based on the eigenvalue problem derived from the linearized equation, but the explicit dispersion relation was not presented. We will include the full derivation of the growth rate formula as a function of wavenumber in the revised manuscript to rigorously demonstrate that the maximum wavelength mode has the highest growth rate. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation proceeds from stated WTG closure to bistability analysis

full rationale

The paper defines a reduced moisture-only model explicitly via the weak-temperature-gradient approximation, then analyzes its nonlinear PDE dynamics to obtain bistability and wavelength selection. No parameters are fitted to the model's own output, no self-citation supplies a load-bearing uniqueness theorem, and the equilibria are not defined in terms of themselves. The derivation chain is self-contained against the model's own equations.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the weak temperature gradient approximation and an unspecified form of the nonlocal coupling operator; no free parameters are explicitly named in the abstract, but the filter length and diffusion coefficient are implicit choices required for the scale selection result.

free parameters (2)

nonlocal filter length
Determines the preferred wavelength in finite-range coupling; its value is chosen to match the characteristic scale of interest.
diffusion coefficient
Controls the stabilizing term that must be overcome by nonlocal destabilization.

axioms (2)

domain assumption Weak temperature gradient approximation holds throughout the domain
Reduces the system to a single moisture variable and is invoked to justify the column-integrated formulation.
domain assumption Nonlocal coupling represents circulation effects from horizontal differential heating
Provides the physical interpretation of the coupling term.

pith-pipeline@v0.9.0 · 5813 in / 1335 out tokens · 20222 ms · 2026-05-24T02:32:29.218120+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

[1]

global coupling

Introduction Clouds play a central role in determining Earth's climate (Siebesma et al. 2020). Among the various cloud phenomena, the roles and mechanisms of organized clouds are not well understood (Neelin et al. 2008; Sherwood et al. 2010; Bony et al. 2015; Mauritsen and Stevens 2015; Reed and Medeiros 2016; Sherwood et al. 2020). Recently, convective s...

work page 2020
[2]

=𝑞"$, 𝑤=0, and 𝑢 can be any value). Green open circles represent the fixed points for the heteroclinic trajectory at 𝑞

Steady state solution of the globally coupled system Finally, we investigate why the distance between clusters approaches infinity in the globally coupled system to highlight the necessity of nonlocal coupling for maintaining a finite inter-cluster distance. To understand the characteristics of the horizontal scale in a globally coupled system, we need to...

work page doi:10.1175/jas-d-15-0170.1 2003

[1] [1]

global coupling

Introduction Clouds play a central role in determining Earth's climate (Siebesma et al. 2020). Among the various cloud phenomena, the roles and mechanisms of organized clouds are not well understood (Neelin et al. 2008; Sherwood et al. 2010; Bony et al. 2015; Mauritsen and Stevens 2015; Reed and Medeiros 2016; Sherwood et al. 2020). Recently, convective s...

work page 2020

[2] [2]

=𝑞"$, 𝑤=0, and 𝑢 can be any value). Green open circles represent the fixed points for the heteroclinic trajectory at 𝑞

Steady state solution of the globally coupled system Finally, we investigate why the distance between clusters approaches infinity in the globally coupled system to highlight the necessity of nonlocal coupling for maintaining a finite inter-cluster distance. To understand the characteristics of the horizontal scale in a globally coupled system, we need to...

work page doi:10.1175/jas-d-15-0170.1 2003