arxiv: 2604.03031 · v1 · submitted 2026-04-03 · 🌊 nlin.PS

Recognition: 2 theorem links

· Lean Theorem

Vegetation Pattern Formation via Energy-Balance-Constrained Modeling

Authors on Pith no claims yet

Pith reviewed 2026-05-13 18:36 UTC · model grok-4.3

classification 🌊 nlin.PS

keywords vegetation pattern formationsemi-arid environmentsenergy balancewater conservationhillslopereaction-diffusionlinear stabilitypattern migration

0 comments

The pith

Energy-balance and water-conservation principles constrain vegetation models that reproduce uphill band migration and aridity-dependent wavelengths on slopes.

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

The paper starts from conservation laws to restrict the possible forms of reaction-diffusion models for vegetation in semi-arid zones before selecting any specific nonlinearities. This produces a family of semilinear equations whose representative member is a fourth-order vegetation PDE coupled to quasi-steady water transport on a hillslope. Linear analysis then isolates three instability routes whose dominance shifts with slope angle. On sloped ground the water-mediated route prevails, yielding the two observed trends of longer wavelengths at higher aridity and uphill band motion; on flat ground the energy-balance route can act alone. A reader would care because the constraints make explicit which physical ingredients are required and how patterns must change when rainfall or terrain varies.

Core claim

Energy-balance and water-conservation principles constrain the admissible class of reaction-diffusion models for vegetation before any specific closure is chosen. A representative Euler-Lagrange semilinear closure produces a fourth-order vegetation equation coupled to quasi-steady water transport. Linear stability analysis identifies three instability mechanisms—classical water-mediated feedback, energy-balance spatial coupling, and water deflection by vegetation gradients—whose balance depends on terrain geometry. On slopes, water-mediated coupling dominates and reproduces the empirical trends of increasing wavelength with aridity and uphill band migration. On flat terrain, energy-balance耦合

What carries the argument

The family of semilinear closures derived from energy-balance and water-conservation principles, represented by an Euler-Lagrange form that yields a fourth-order vegetation PDE coupled to quasi-steady water transport.

If this is right

On slopes the dominant water-mediated coupling produces vegetation bands whose wavelength increases with aridity.
Vegetation bands on slopes migrate uphill under the same mechanism.
On flat terrain energy-balance spatial coupling alone suffices to drive pattern-forming instability.
Numerical simulations of the full nonlinear system confirm the linear stability predictions.
Continuation methods reveal a narrow hysteresis region consistent with subcritical bifurcation.

Where Pith is reading between the lines

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

The same conservation-law approach could be extended to two-dimensional domains to predict spots, labyrinths, and gaps.
Climate-driven increases in aridity would shift pattern wavelengths and stability thresholds in directions given by the model.
Analogous energy-balance and conservation constraints might reduce arbitrary choices in models of other ecological patterns such as fairy circles or mussel beds.
Relaxing the quasi-steady water assumption could expose transient dynamics or additional instabilities not visible in the present reduction.

Load-bearing premise

The family of semilinear closures motivated by energy-balance and water-conservation principles includes forms that accurately represent the physical processes, together with the quasi-steady water transport approximation.

What would settle it

Field measurements showing that vegetation bands on slopes migrate downhill or that pattern wavelengths decrease rather than increase with aridity would falsify the predicted dominance of water-mediated coupling.

read the original abstract

Vegetation in semi-arid environments self-organizes into striking spatial patterns -- bands, spots, labyrinths, and gaps -- with characteristic wavelengths on the order of tens to hundreds of meters. Existing reaction-diffusion models postulate nonlinearities and transport laws from qualitative physical reasoning, making it hard to distinguish essential structural features from artifacts of the chosen forms. Here we show how energy-balance and water-conservation principles can constrain the admissible model class before a specific closure is chosen. These constraints motivate a family of semilinear closures; an Euler--Lagrange representative yields a fourth-order vegetation equation coupled to quasi-steady water transport on a one-dimensional hillslope. Linear stability analysis identifies three instability mechanisms: classical water-mediated feedback, energy-balance spatial coupling, and water deflection by vegetation gradients. Their balance depends on terrain geometry. On slopes, the water-mediated coupling dominates and the model reproduces two empirical observations: pattern wavelength increases with aridity, and vegetation bands migrate uphill. On flat terrain, the energy-balance spatial coupling can drive instability independently. Numerical simulations confirm the linear predictions, and exploratory continuation reveals a narrow hysteresis region consistent with subcritical bifurcation.

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 constrains vegetation pattern models via energy and water conservation to get a fourth-order PDE that recovers observed trends on slopes and flat ground.

read the letter

The main thing to know is that this work starts from energy-balance and water-conservation principles to limit the admissible forms of a vegetation model before choosing a specific closure. That produces a family of semilinear terms, and an Euler-Lagrange pick yields a fourth-order vegetation equation coupled to quasi-steady water transport on a one-dimensional hillslope. Linear stability then isolates three instability routes whose balance shifts with terrain: classical water-mediated feedback, energy-balance spatial coupling, and deflection by vegetation gradients. On slopes the water route dominates and the model gets the two reported empirical trends—wavelength growing with aridity and bands migrating uphill—while on flat ground the energy term can drive instability by itself. Simulations match the linear predictions and show a narrow hysteresis band. This constrained route is new relative to the usual reaction-diffusion papers that begin with ad-hoc nonlinearities. It does a clean job separating structural features from closure choices and tying them to measurable terrain effects. The soft spots are modest but real. The semilinear family is motivated by the conservation laws, yet the paper still has to show that results are not overly sensitive to which member is chosen and that the forms actually represent the small-scale processes better than existing alternatives. The quasi-steady water approximation is standard but could limit applicability under strongly transient rainfall. Quantitative comparisons to field data beyond the two trends would strengthen the case. This is for researchers in mathematical ecology and pattern formation who want a more principled way to build models rather than guess terms. A reader working on conservation-law approaches to biological PDEs would find the framework useful to adapt. It deserves a serious referee because the method is distinct enough from the literature to repay detailed checking of the derivations and numerics.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a constrained modeling framework for vegetation pattern formation in semi-arid environments by imposing energy-balance and water-conservation principles to restrict admissible reaction-diffusion models to a family of semilinear closures. An Euler-Lagrange representative produces a fourth-order vegetation equation coupled to quasi-steady water transport on a one-dimensional hillslope. Linear stability analysis identifies three distinct instability mechanisms (water-mediated feedback, energy-balance spatial coupling, and water deflection by vegetation gradients) whose relative importance depends on terrain slope. On slopes the water-mediated mechanism dominates and the model recovers the empirical trends of increasing pattern wavelength with aridity and uphill migration of bands; on flat terrain the energy-balance coupling alone can drive instability. Numerical simulations and continuation methods confirm the linear predictions and reveal a narrow hysteresis region.

Significance. If the central derivation and closure assumptions hold, the work supplies a principled route to reduce arbitrariness in vegetation pattern models, replacing qualitative choices of nonlinearities with conservation-law constraints. The explicit separation of three instability mechanisms and their terrain dependence offers mechanistic insight that is difficult to obtain from conventional ad-hoc models. Reproduction of two well-documented empirical trends without heavy parameter tuning, together with the identification of an independent flat-terrain instability, constitutes a falsifiable advance that could improve predictions of ecosystem response to changing aridity.

major comments (2)

[Linear stability analysis] The linear stability analysis section must derive the dispersion relation for the flat-terrain limit explicitly (setting slope to zero) to demonstrate that the energy-balance spatial coupling produces a positive growth rate at finite wavenumber independently of the water-mediated term; without this explicit reduction the claim that the mechanism operates independently remains unverified.
[Model formulation] The quasi-steady water-transport approximation is load-bearing for both the slope and flat cases; the manuscript should quantify the timescale separation (vegetation vs. water) and show that relaxing it does not alter the sign of the dominant eigenvalues or the reported migration direction.

minor comments (2)

[Numerical methods] The abstract states that numerical simulations confirm the linear predictions, yet no information is given on the spatial discretization, time-stepping method, or domain size used; these details belong in the methods section.
[Model derivation] The family of semilinear closures is introduced but only one Euler-Lagrange representative is analyzed; a brief sensitivity test to alternative closures (e.g., different polynomial degrees) would strengthen the claim that the reported instabilities are structural rather than closure-specific.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and positive report. We address each major comment below and will revise the manuscript to incorporate the requested clarifications.

read point-by-point responses

Referee: [Linear stability analysis] The linear stability analysis section must derive the dispersion relation for the flat-terrain limit explicitly (setting slope to zero) to demonstrate that the energy-balance spatial coupling produces a positive growth rate at finite wavenumber independently of the water-mediated term; without this explicit reduction the claim that the mechanism operates independently remains unverified.

Authors: We agree that an explicit derivation strengthens the claim. In the revised manuscript we will set the slope parameter to zero, isolate the energy-balance spatial-coupling terms, and derive the resulting dispersion relation, confirming a band of positive growth rates at finite wavenumbers that persists when the water-mediated coefficients are removed. revision: yes
Referee: [Model formulation] The quasi-steady water-transport approximation is load-bearing for both the slope and flat cases; the manuscript should quantify the timescale separation (vegetation vs. water) and show that relaxing it does not alter the sign of the dominant eigenvalues or the reported migration direction.

Authors: We will add a paragraph quantifying the timescale separation using representative parameter values from the semi-arid literature, showing vegetation growth timescales are two to three orders of magnitude slower than water redistribution. A complete relaxation of the quasi-steady assumption would require solving the fully time-dependent coupled system, which lies outside the present scope; we will note this limitation while confirming that the sign of the dominant eigenvalues and uphill migration direction are preserved under the standard separation. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation from external conservation principles

full rationale

The paper starts from energy-balance and water-conservation principles as external physical constraints to restrict the model class to semilinear closures. An Euler-Lagrange representative then produces the fourth-order vegetation equation, whose linear stability analysis recovers the reported trends on slopes and flat terrain. No load-bearing step reduces by construction to a fitted parameter, self-citation chain, or renaming of known results; the central claims follow directly from the constrained equations and standard stability methods without internal reduction to inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

Review based on abstract only; specific free parameters and invented entities are not enumerated in the provided text. The central approach rests on conservation principles as domain assumptions.

free parameters (1)

parameters within chosen semilinear closure
The specific representative from the admissible family likely requires parameter choices, though none are detailed in the abstract.

axioms (2)

domain assumption Energy-balance and water-conservation principles constrain the admissible class of vegetation models before a specific closure is chosen
Stated directly in the abstract as the motivating step for the family of semilinear closures.
domain assumption Quasi-steady approximation for water transport on the hillslope
Invoked in the abstract to couple the fourth-order vegetation equation to water transport.

pith-pipeline@v0.9.0 · 5491 in / 1407 out tokens · 75486 ms · 2026-05-13T18:36:58.507335+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

energy-balance and water-conservation principles constrain the admissible model class... semilinear closures; an Euler–Lagrange representative yields a fourth-order vegetation equation
IndisputableMonolith/Foundation/BranchSelection.lean branch_selection unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

score functional E[u]=∫(u−κ/2 G²)... Euler–Lagrange operator... even local symbol; nonpositive quartic

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages

[1]

To obtain a concrete PDE, we need an evolution law for the vegetation field that is built from these ingredients

Semilinear variational closure and hillslope model.The preceding section developed a constrained class of admissible model ingredients, including sign-constrained energy mismatches, gradient-expansion spatial structure, and quasi- steady water, without specifying the vegetation dynamics. To obtain a concrete PDE, we need an evolution law for the vegetatio...

work page
[2]

no kernel

Linear stability analysis.We now ask whether the uniform vegetated state derived in sections 2 and 3 is stable to spatially periodic perturbations. We compute the uniform steady state, linearize the vegetation–water system in Fourier space, and, for the Euler–Lagrange representative, obtain an exact decomposition of the growth rate into a local even polyn...

work page
[3]

Theβ 1uwterm in our model plays the same structural role as the Klausmeier uptakewn 2; the dif- ferent power ofnreflects different density-dependence assumptions

Comparison with existing models.A structural comparison with the Klausmeier model (1.1), detailed in E, highlights several differences. Theβ 1uwterm in our model plays the same structural role as the Klausmeier uptakewn 2; the dif- ferent power ofnreflects different density-dependence assumptions. Both forms are consistent with the sign constraints from T...

work page
[4]

We use pseudospectral spatial discretization [5] on a periodic domain with ETDRK4 time integration [6, 15]; full details of the numerical scheme are given in G

Numerical simulations.We solve the nonlinear system numerically to test whether the linear predictions survive into the finite-amplitude regime and to explore the nonlinear branch structure. We use pseudospectral spatial discretization [5] on a periodic domain with ETDRK4 time integration [6, 15]; full details of the numerical scheme are given in G. Table...

work page
[5]

Discussion.The framework operates at three layers of constraint: energy- balance sign constraints, gradient-expansion form for spatial terms, and a semilinear variational closure. The first two layers restrict the admissible model class indepen- dently of the closure; the semilinear family built from them shares a common water transfer function Φ(k) and a...

work page
[6]

For slow, overdamped systems, variational closure is a natural modeling choice

models are constructed by writing a free energy from thermodynamic arguments and postulating gradient dynamics. For slow, overdamped systems, variational closure is a natural modeling choice. Vegetation dynamics, governed by growth and mortality on seasonal to decadal timescales, are slow, though slowness alone does not guarantee gradient-flow structure; ...

work page
[7]

Conclusion.Energy balance, water conservation, and gradient expansion sharply constrain the admissible class of dryland vegetation models before any clo- sure is chosen. The Euler–Lagrange representative of this constrained class yields an exactly decomposable dispersion relation with three instability mechanisms, one of which—energy-balance spatial coupl...

work page
[8]

Allen, S. M. and Cahn, J. W. (1979). A microscopic theory for antiphase bound- ary motion and its application to antiphase domain coarsening.Acta Metallurgica, 27(6):1085–1095

work page 1979
[9]

Borgogno, F., D’Odorico, P., Laio, F., and Ridolfi, L. (2009). Mathematical models of vegetation pattern formation in ecohydrology.Reviews of Geophysics, 47(1):RG1005

work page 2009
[10]

Cahn, J. W. and Hilliard, J. E. (1958). Free energy of a nonuniform system. I. Interfacial free energy.The Journal of Chemical Physics, 28(2):258–267. VEGETATION PATTERNS VIA ENERGY-BALANCE-CONSTRAINED MODELING25

work page 1958
[11]

Campbell, G. S. and Norman, J. M. (1998).An Introduction to Environmental Biophysics. Springer, New York, 2nd edition

work page 1998
[12]

Y., Quarteroni, A., and Zang, T

Canuto, C., Hussaini, M. Y., Quarteroni, A., and Zang, T. A. (2006).Spectral Methods: Fundamentals in Single Domains. Springer, Berlin

work page 2006
[13]

Cox, S. M. and Matthews, P. C. (2002). Exponential time differencing for stiff systems.Journal of Computational Physics, 176(2):430–455

work page 2002
[14]

Deblauwe, V., Barbier, N., Couteron, P., Lejeune, O., and Bogaert, J. (2008). The global biogeography of semi-arid periodic vegetation patterns.Global Ecology and Biogeography, 17(6):715–723

work page 2008
[15]

Deblauwe, V., Couteron, P., Bogaert, J., and Barbier, N. (2012). Determinants and dynamics of banded vegetation pattern migration in arid climates.Ecological Monographs, 82(1):3–21

work page 2012
[16]

Gandhi, P., Iams, S., Bonetti, S., and Silber, M. (2019). Vegetation pattern formation in drylands. In D’Odorico, P., Porporato, A., and Wilkinson Runyan, C., editors,Dryland Ecohydrology, pages 469–509. Springer, Cham

work page 2019
[17]

Gandhi, P., Oline, M., and Silber, M. (2025). A flow-kick model of dryland vegetation patterns: the impact of rainfall variability on resilience.SIAM Journal on Applied Dynamical Systems, 24(4):3092–3126

work page 2025
[18]

Gilad, E., von Hardenberg, J., Provenzale, A., Shachak, M., and Meron, E. (2004). Ecosystem engineers: from pattern formation to habitat creation.Phys- ical Review Letters, 93(9):098105

work page 2004
[19]

Gilad, E., von Hardenberg, J., Provenzale, A., Shachak, M., and Meron, E. (2007). A mathematical model for plants as ecosystem engineers.Journal of The- oretical Biology, 244(4):680–691

work page 2007
[20]

Gowda, K., Chen, Y., Iams, S., and Silber, M. (2016). Assessing the robustness of spatial pattern sequences in a dryland vegetation model.Proceedings of the Royal Society A, 472(2187):20150893

work page 2016
[21]

Jones, H. G. (2013).Plants and Microclimate: A Quantitative Approach to En- vironmental Plant Physiology. Cambridge University Press, 3rd edition

work page 2013
[22]

and Trefethen, L

Kassam, A.-K. and Trefethen, L. N. (2005). Fourth-order time-stepping for stiff PDEs.SIAM Journal on Scientific Computing, 26(4):1214–1233

work page 2005
[23]

Klausmeier, C. A. (1999). Regular and irregular patterns in semiarid vegetation. Science, 284(5421):1826–1828

work page 1999
[24]

(2012).Nonlinear Physics of Ecosystems

Meron, E. (2012).Nonlinear Physics of Ecosystems. CRC Press

work page 2012
[25]

Meron, E. (2018). From patterns to function in living systems: Dryland ecosys- tems as a case study.Annual Review of Condensed Matter Physics, 9:79–103

work page 2018
[26]

C., van Langevelde, F., HilleRisLambers, R., van de Koppel, J., Kumar, L., Prins, H

Rietkerk, M., Boerlijst, M. C., van Langevelde, F., HilleRisLambers, R., van de Koppel, J., Kumar, L., Prins, H. H. T., and de Roos, A. M. (2002). Self-organization of vegetation in arid ecosystems.The American Naturalist, 160(4):524–530

work page 2002
[27]

and Doelman, A

Sewalt, L. and Doelman, A. (2017). Spatially periodic multipulse patterns in a generalized Klausmeier–Gray–Scott model.SIAM Journal on Applied Dynamical Systems, 16(2):1113–1163

work page 2017
[28]

Sherratt, J. A. (2005). An analysis of vegetation stripe formation in semi-arid landscapes.Journal of Mathematical Biology, 51(2):183–197

work page 2005
[29]

and Hohenberg, P

Swift, J. and Hohenberg, P. C. (1977). Hydrodynamic fluctuations at the con- vective instability.Physical Review A, 15(1):319–328

work page 1977
[30]

Valentin, C., d’Herb` es, J.-M., and Poesen, J. (1999). Soil and water components of banded vegetation patterns.Catena, 37(1–2):1–24

work page 1999
[31]

von Hardenberg, J., Meron, E., Shachak, M., and Zarmi, Y. (2001). Diversity of 26C. M. TOPAZ vegetation patterns and desertification.Physical Review Letters, 87(19):198101. Appendix A. Linearization of the water equation. To obtain the water transfer function Φ(k), we substituteU=U 0 +ϵ ˆU eikX and W=W 0 +ϵ ˆW eikX into the quasi-steady water equation (3....

work page 2001
[32]

nonlinear remainder

It modifies the balance between destabilizing and stabilizing contributions to the local part ofσ(k), not the propagation speed. Appendix C. Cancellations in the linearized Euler–Lagrange operator. This section verifies two structural properties of the linearized Euler–Lagrange operator. First, all odd-in-kcontributions from the local part cancel. Second,...

work page