pith. sign in

arxiv: 2509.04766 · v2 · submitted 2025-09-05 · 🧮 math.AP

Stability and Self-Organized Patterns in Coupled Ecohydrological--Fire Dynamics: A Model of Vegetation--Rainfall--Bushfire Interactions

Serena Dipierro , Enrico Valdinoci This is my paper

Pith reviewed 2026-05-18 19:31 UTC · model grok-4.3

classification 🧮 math.AP MSC 35K5735B3535C0792D40
keywords reaction-diffusion systemsstability analysistraveling wavesecohydrologybushfire dynamicsTuring instabilityvegetation patternseigenvalue crossing
0
0 comments X

The pith

Diffusion stabilizes an unstable uniform state in a vegetation-rainfall-bushfire model and produces stable traveling wave trains via a nonsingular eigenvalue crossing.

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

The paper studies a three-component reaction-diffusion system that couples water reservoirs, vegetation biomass, and bushfire activity. It establishes that diffusion can turn an unstable homogeneous equilibrium stable against spatially nonuniform perturbations while also creating periodic traveling trains. These trains arise as orbits in the linearized system from a nonsingular crossing of the imaginary axis, with a third eigenvalue staying negative to guarantee linear stability of the waves. The mechanism differs from classical Turing or Hopf bifurcations because the model lacks distinct activator-inhibitor roles and does not require unequal diffusion speeds. The work further identifies a Turing instability that produces slow-frequency oscillations under plant competition in low-rainfall regimes.

Core claim

In the three-component reaction-diffusion model for vegetation, rainfall, and bushfire, diffusion stabilizes the unstable homogeneous equilibrium with respect to spatially nonuniform perturbations and generates traveling trains as periodic orbits of the linearized system arising from a nonsingular eigenvalue crossing of the imaginary axis while a third eigenvalue remains real and negative, thereby ensuring linear stability for monochromatic waves.

What carries the argument

Three-component reaction-diffusion system whose linearized operator exhibits a nonsingular pair of eigenvalues crossing the imaginary axis while a third eigenvalue stays negative.

If this is right

  • Homogeneous vegetation-fire equilibria can remain stable under diffusion for ranges of rainfall and ignition parameters.
  • Periodic traveling wave trains emerge as linearly stable solutions without requiring classical activator-inhibitor separation.
  • In the plant-competition variant, slow-frequency spatial oscillations appear via Turing instability when rainfall is small.
  • Stability regions in parameter space can be mapped explicitly by tracking eigenvalue crossings.

Where Pith is reading between the lines

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

  • Observations of banded vegetation or periodic fire fronts in semi-arid landscapes could be re-examined as possible outcomes of this diffusion-induced stabilization rather than traditional Turing mechanisms.
  • The eigenvalue-crossing construction might extend to other three-species ecological models where diffusion is expected to regularize rather than destabilize uniform states.
  • Long-term simulations of the full nonlinear system could test whether the linearly stable traveling trains persist or interact to form more complex patterns.

Load-bearing premise

The specific functional forms chosen for vegetation growth, fire ignition, and water uptake accurately represent the dominant processes in real ecosystems.

What would settle it

A numerical simulation or field measurement in the predicted stable parameter regime in which small spatial perturbations grow exponentially instead of decaying.

Figures

Figures reproduced from arXiv: 2509.04766 by Enrico Valdinoci, Serena Dipierro.

Figure 1
Figure 1. Figure 1: Plot of 27 trajectories of (1.2) corresponding to the parameters α := β := γ := δ := ε := η := ζ := 1, with final time T ∈ {5, 10, 50}. These assumptions lead to the equation ∂tf = f(αv − βw) + c∆f. Notice that the product of fv (respectively, fw) on the right-hand side of the equation above corresponds to a “random encounter” between fire and vegetation (respectively, fire and water). Also, to model the e… view at source ↗
Figure 2
Figure 2. Figure 2: Plot of 27 trajectories of (1.2) corresponding to the parameters α := 2, β := γ := δ := η := ζ := 1, and ε := 1 10 , with final time T ∈ {5, 10, 50}. We stress that, to keep the analysis as simple as possible, we have assumed the terrain to be flat and homogeneous (variations of fuel and slope can be taken into account, for example, by replacing the Laplacian with an inhomogeneous diffusion operator) [PIT… view at source ↗
Figure 3
Figure 3. Figure 3: A “tiger bush” plateau in Niger, with vegetation dominated by Combretum micranthum and Guiera senegalensis. United States Geological Survey (Public Domain). While similar in some respects, this system of equations is structurally very different from the one for three competing species engaged in a rock-paper-scissors game that has been used in mathematical biology, for example, to describe dynamics in liza… view at source ↗
Figure 4
Figure 4. Figure 4: Central Manitoba and Quebec, Canada: images captured by the Operational Land Imager-2 (OLI-2 ) on NASA’s satellite Landsat 9 (May 23, 2025, and July 14, 2025) (Public Domain). We now return to the analysis of the system of partial differential equations (1.1), with the aim of detecting both the formation of patterns created by diffusion and the stabilization of the [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: False-color images acquired by the Moderate Resolution Imaging Spectroradiome￾ter (MODIS) on board NASA’s Terra satellite: burn scars in Australia (February 3, 2020), Ore￾gon (September 27, 2020), Kansas (April 12, 2023), and Bolivia (December 9, 2023) (Public Domain). We observe that the result showcased by Theorem 1.2 exhibits several structural differences compared to “classical” cases of Turing pattern… view at source ↗
Figure 6
Figure 6. Figure 6: Displacement of the eigenvalues of the matrix A in (3.2) when α := 2, β := γ := δ := η := ζ := 1, ε := 1 10 , and c|k| 2 = d|k| 2 ∈ [0, 2]. The periodic solutions for the linearized system in Theorem 1.3 arise when a pair of eigenvalues crosses the imaginary axis from right to left as µ increases. Note that the third eigenvalue remains real and negative. The characteristic polynomial of A is P(λ) = λ 3 + a… view at source ↗
Figure 7
Figure 7. Figure 7: Displacement of the eigenvalues of the matrix L in (5.1) when α := β := δ := ε := η := ζ := 1, µ := γ, ς := ε/2, ℓ := ς µv⋆ . The complex eigenvalues are displayed as the parameter γ varies over [0, 1]. Colors follow a rainbow scheme: red for high γ and violet for low γ (in the situation in which we are reducing γ, the image must be interpreted “from red to violet”). The instability established in Theorem … view at source ↗
read the original abstract

This paper investigates the conditions for the stability and emergence of patterns in a new three-component reaction-diffusion system. The system describes the coexistence and interaction of water reservoirs, vegetation, and bushfire activity in a given ecosystem. We perform a detailed stability analysis to determine the parameter space where an unstable homogeneous equilibrium becomes stable with respect to spatially nonuniform perturbations. We also use diffusion to generate traveling trains in the form of periodic orbits of the linearized system. These orbits are remnants of an unstable equilibrium in the absence of diffusion and arise from a nonsingular eigenvalue crossing of the imaginary axis, while a third eigenvalue remains real and negative, thereby ensuring linear stability for monocromatic waves. These phenomena differ from ``classical'' Turing and Hopf bifurcations, as the model does not involve distinct ``activators'' and ``inhibitors'', and the effects observed are not the byproduct of diffusion with necessarily differing speeds. Also, differently from the classical Turing pattern, the role of diffusion in this context is to stabilize, rather than destabilize, homogeneous equilibria. We also consider the case of plant competition, showing a suitable form of Turing instability for slow-frequency oscillations in a small rainfall regime.

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

Summary. The manuscript introduces a three-component reaction-diffusion system for vegetation, water reservoirs, and bushfire activity. It performs a linear stability analysis to identify parameter regimes in which an unstable homogeneous equilibrium is stabilized against spatially nonuniform perturbations by diffusion. The work further shows that diffusion generates traveling wave trains as periodic orbits of the linearized system via a nonsingular eigenvalue crossing of the imaginary axis (with a third eigenvalue remaining real and negative), ensuring linear stability of monochromatic waves. These mechanisms are contrasted with classical Turing and Hopf bifurcations, and a separate Turing instability is reported for plant competition under low-rainfall conditions.

Significance. If the central claims are rigorously established, the results would provide a non-classical route to pattern formation in which diffusion stabilizes rather than destabilizes selected modes of an ODE-unstable equilibrium, without requiring activator-inhibitor structure or disparate diffusion rates. This could inform modeling of self-organized vegetation-fire patterns in ecohydrological systems and would constitute a mathematically interesting extension of reaction-diffusion theory.

major comments (2)
  1. [Abstract] Abstract: The claim that diffusion causes an 'unstable homogeneous equilibrium [to] become stable with respect to spatially nonuniform perturbations' while the equilibrium remains unstable in the absence of diffusion is load-bearing for the stability regions and traveling-wave results. In any reaction-diffusion system the dispersion relation satisfies λ(k=0) = eigenvalues of the reaction Jacobian J; diffusion cannot change Re(λ(0)). If Re(λ(0)) > 0, then sup_k Re(λ(k)) ≥ Re(λ(0)) > 0 and the homogeneous state is unstable in the PDE. The manuscript must clarify whether 'stable with respect to spatially nonuniform perturbations' is intended only as a statement about the sign of Re(λ(k)) for k > 0, and must explain how this affects the interpretation of the claimed stability regions and the linear stability of the monochromatic waves.
  2. [Abstract] Abstract: The assertion that traveling trains arise from a 'nonsingular eigenvalue crossing of the imaginary axis, while a third eigenvalue remains real and negative' requires an explicit dispersion relation or characteristic equation (including the precise functional forms of the reaction terms and the diffusion matrix) to verify that the crossing is indeed nonsingular and that the remaining eigenvalue stays negative for the relevant wave numbers. Without these, the linear stability conclusion for the monochromatic waves cannot be confirmed.
minor comments (1)
  1. [Abstract] The abstract states that the phenomena 'differ from classical Turing and Hopf bifurcations' because the model lacks distinct activators and inhibitors and does not rely on differing diffusion speeds. A brief comparison with the standard dispersion-relation conditions for those bifurcations would help readers locate the novelty.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and insightful comments, which help clarify the presentation of our results on diffusion-induced effects in the three-component ecohydrological-fire model. We respond to each major comment below.

read point-by-point responses

  1. Referee: The claim that diffusion causes an 'unstable homogeneous equilibrium [to] become stable with respect to spatially nonuniform perturbations' while the equilibrium remains unstable in the absence of diffusion is load-bearing for the stability regions and traveling-wave results. In any reaction-diffusion system the dispersion relation satisfies λ(k=0) = eigenvalues of the reaction Jacobian J; diffusion cannot change Re(λ(0)). If Re(λ(0)) > 0, then sup_k Re(λ(k)) ≥ Re(λ(0)) > 0 and the homogeneous state is unstable in the PDE. The manuscript must clarify whether 'stable with respect to spatially nonuniform perturbations' is intended only as a statement about the sign of Re(λ(k)) for k > 0, and must explain how this affects the interpretation of the claimed stability regions and the linear stability of the monochromatic waves.

    Authors: We agree that the dispersion relation is continuous at k=0 and that diffusion cannot alter the eigenvalues at k=0, so an ODE-unstable equilibrium remains unstable in the PDE. Our phrasing 'stable with respect to spatially nonuniform perturbations' was meant to refer specifically to the possibility that Re(λ(k)) < 0 for k > 0 due to the diffusive terms, even while Re(λ(0)) > 0. This contrasts with classical Turing instability (where diffusion destabilizes k > 0 while k = 0 is stable) and is the non-classical feature we highlight. We will revise the abstract and the linear stability section to state this distinction explicitly and to note that the claimed stability regions concern the damping of nonuniform modes, which can still support the emergence of traveling waves or other structures despite uniform instability. The linear stability analysis of the monochromatic waves is performed separately by linearizing about the wave profile, not the homogeneous equilibrium, so that part of the interpretation is unaffected. We have added a clarifying paragraph in Section 2. revision: yes

  2. Referee: The assertion that traveling trains arise from a 'nonsingular eigenvalue crossing of the imaginary axis, while a third eigenvalue remains real and negative' requires an explicit dispersion relation or characteristic equation (including the precise functional forms of the reaction terms and the diffusion matrix) to verify that the crossing is indeed nonsingular and that the remaining eigenvalue stays negative for the relevant wave numbers. Without these, the linear stability conclusion for the monochromatic waves cannot be confirmed.

    Authors: We accept that the manuscript would be strengthened by including the explicit dispersion relation. In the revised version we will add the full characteristic equation obtained from the linearized reaction-diffusion system, together with the precise reaction Jacobian and diagonal diffusion matrix. The cubic polynomial in λ is derived in the usual way by substituting Fourier modes (or the appropriate traveling-wave ansatz) and we will verify the nonsingular crossing by showing that the derivative of the real part of the complex-conjugate pair with respect to the bifurcation parameter is nonzero at the critical value. The negativity of the remaining real eigenvalue will be confirmed by direct evaluation or the Routh-Hurwitz criterion over the relevant interval of wave numbers. These details will appear in a new subsection of the stability analysis and in an appendix containing the explicit reaction terms. revision: yes

Circularity Check

0 steps flagged

No circularity: direct linearization of stated reaction-diffusion system yields eigenvalue conditions without reduction to inputs or self-citations

full rationale

The paper conducts a standard linear stability analysis of a three-component reaction-diffusion PDE by computing the spectrum of the operator obtained from the reaction Jacobian minus k² times the diffusion matrix. The reported conditions for stabilization with respect to nonuniform perturbations and the existence of traveling trains follow directly from solving the resulting cubic characteristic equation for each wave number k and identifying parameter regimes where Re(λ(k)) changes sign for k > 0 while the k = 0 mode remains as given by the ODE Jacobian. No parameter is fitted to data and then relabeled as a prediction, no ansatz is smuggled via prior self-citation, and the derivation does not define any quantity in terms of itself. The analysis is therefore self-contained against the model's explicit equations and does not reduce to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The abstract does not enumerate free parameters, background axioms, or new postulated entities. The model is presented as a new system whose reaction terms and diffusion coefficients are taken as given inputs for the stability analysis.

pith-pipeline@v0.9.0 · 5749 in / 1230 out tokens · 33285 ms · 2026-05-18T19:31:38.135658+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

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

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

26 extracted references · 26 canonical work pages

  1. [1]

    Banasiak, Y

    J. Banasiak, Y. Dumont, I. V. Yatat Djeumen, Spreading Speeds and Traveling Waves for Monotone Systems of Impulsive Reaction–Diffusion Equations: Application to Tree–Grass Interactions in Fire-prone Savannas . Different. Eq. Dynam. Syst. 31 , 547--580, (2023)

    work page 2023

  2. [2]

    Borgogno, P

    F. Borgogno, P. D'Odorico, F. Laio, L. Ridolfi, Mathematical models of vegetation pattern formation in ecohydrology . Rev. Geophys 47 , RG1005 (2009)

    work page 2009

  3. [3]

    Carter, A

    P. Carter, A. Doelman, Traveling stripes in the Klausmeier model of vegetation pattern formation . SIAM J. Appl. Math. 78 , 6, 3213--3237 (2018)

    work page 2018

  4. [4]

    Carter, A

    P. Carter, A. Doelman, A. Iuorio, F. Veerman, Travelling pulses on three spatial scales in a Klausmeier-type vegetation-autotoxicity model . Nonlinearity 37 , 9, 095008 (2024)

    work page 2024

  5. [5]

    S. Clar, B. Drossel, F. Schwabl, Forest fires and other examples of self-organized criticality . J. Phys.: Condens. Matter 8 , 6803 (1996)

    work page 1996

  6. [6]

    Dipierro, E

    S. Dipierro, E. Valdinoci, Elliptic partial differential equations from an elementary viewpoint---a fresh glance at the classical theory . World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, xvii+651 pp., ISBN 978-981-12-9079-4 (2024)

    work page 2024

  7. [7]

    E. J. B. van Etten, R. A. Davis, T. S. Doherty, Fire in Semi-Arid Shrublands and Woodlands: Spatial and Temporal Patterns in an Australian Landscape , Front. Ecol. Evol. 9 (2021)

    work page 2021

  8. [8]

    M. O. Gani, Md. A. I. Arif, A. S. Howladar, A. Bashar, Stability analysis of periodic traveling waves in a model of vegetation patterns in semi-arid ecosystems . Model. Earth Systems Environm. 7 , 1511--1522 (2021)

    work page 2021

  9. [9]

    Ge, The hidden order of Turing patterns in arid and semi-arid vegetation ecosystems , Proce

    Z. Ge, The hidden order of Turing patterns in arid and semi-arid vegetation ecosystems , Proce. Nat. Acad. Sciences United States of America 120 , 42, e2306514120 (2023)

    work page 2023

  10. [10]

    von Hardenberg, E

    J. von Hardenberg, E. Meron, M. Shachak, Y. Zarmi, Diversity of Vegetation Patterns and Desertification . Phys. Rev. Lett. 87 , 198101 (2001)

    work page 2001

  11. [11]

    C. A. Hiemstra, G. E. Liston, W. A. Reiners, Observing, modelling, and validating snow redistribution by wind in a Wyoming upper treeline landscape . Ecol. Modelling 197 , 35--51 (2006)

    work page 2006

  12. [12]

    R. E. Keane, R. A. Loehman, L. M. Holsinger, The FireBGCv2 Landscape Fire Succession Model: A Research Simulation Platform for Exploring Fire and Vegetation Dynamics . Rocky Mountain Research Station General Technical Report, RMRS-GTR-255 (2011)

    work page 2011

  13. [13]

    B. C. Kirkup, M. A. Riley, Antibiotic-mediated antagonism leads to a bacterial game of rock--paper--scissors in vivo . Nature 428 , 412 (2004)

    work page 2004

  14. [14]

    C. A. Klausmeier, Regular and Irregular Patterns in Semiarid Vegetation . Science 284 , 5421, 1826--1828 (1999)

    work page 1999

  15. [15]

    Lefever, N

    R. Lefever, N. Barbier, P. Couteron, O. Lejeune, Deeply gapped vegetation patterns: on crown/root allometry, criticality and desertification . J. Theoret. Biol. 261 , 2, 194—209 (2009)

    work page 2009

  16. [16]

    E. L. Loudermilk, J. J. O'Brien, S. L. Goodrick, R. R. Linn, N. S. Skowronski, J. K. Hiers, Vegetation's influence on fire behavior goes beyond just being fuel . Fire Ecol. 18 , 9 (2022)

    work page 2022

  17. [17]

    Meron, Pattern-formation approach to modelling spatially extended ecosystems

    E. Meron, Pattern-formation approach to modelling spatially extended ecosystems . Ecol. Modelling 234 , 70--82 (2012)

    work page 2012

  18. [18]

    Moustakas, Fire acting as an increasing spatial autocorrelation force: Implications for pattern formation and ecological facilitation

    A. Moustakas, Fire acting as an increasing spatial autocorrelation force: Implications for pattern formation and ecological facilitation . Ecolog. Complexity 21 , 142--149 (2015)

    work page 2015

  19. [19]

    S. M. Owen, C. H. Sieg, A. J. S\'anchez Meador, P. Z. Ful\'e, J. M. Iniguez, L. S. Baggett, P. J. Fornwalt, M. A. Battaglia, Spatial patterns of ponderosa pine regeneration in high-severity burn patches . Forest Ecol. Manag. 405 , 134--149 (2017)

    work page 2017

  20. [20]

    Rietkerk, M

    M. Rietkerk, M. C. Boerlijst, F. van Langevelde, R. HilleRisLambers, J. van de Koppel, L. Kumar, H. H. T. Prins, A. M. de Roos, Self‐Organization of Vegetation in Arid Ecosystems . Amer. Naturalist 160 , 4, 524--530 (2002)

    work page 2002

  21. [21]

    Rietkerk, S

    M. Rietkerk, S. C. Dekker, P C. de Ruiter, J. van de Koppel, Self-Organized Patchiness and Catastrophic Shifts in Ecosystems . Science 305 , 5692, 1926--1929 (2004)

    work page 1926

  22. [22]

    T. M. Scanlon, K. K. Caylor, S. A. Levin, I. Rodriguez-Iturbe, Positive feedbacks promote power-law clustering of Kalahari vegetation . Nature 449 , 209--212 (2007)

    work page 2007

  23. [23]

    J. A. Sherratt, Turing Patterns in Deserts . In: S. B. Cooper, A. Dawar, B. L\"owe (eds), How the World Computes . Lecture Notes in Computer Science, vol 7318. Springer, Berlin (2012)

    work page 2012

  24. [24]

    Sinervo, C

    B. Sinervo, C. M. Lively, The rock--paper--scissors game and the evolution of alternative male strategies . Nature 380 , 240 (1996)

    work page 1996

  25. [25]

    D. G. Sprugel, Dynamic structure of wave-regenerated Abies balsamea forests in the northeastern United States . J. Ecology 64 , 3, 889--911 (1976)

    work page 1976

  26. [26]

    X. B. Wu, S.R. Archer, Scale-Dependent Influence of Topography-Based Hydrologic Features on Patterns of Woody Plant Encroachment in Savanna Landscapes . Landscape Ecol. 20 , 733--742 (2005)

    work page 2005