Klausmeier-Gray-Scott Model

• The Klausmeier-Gray-Scott model is a vegetation-water interaction framework that uses an autocatalytic uptake term and water diffusion to simulate dryland pattern dynamics.
• It interpolates between slope-driven and diffusion-driven regimes by incorporating fractional and nonlocal transport to capture terrain effects.
• Bifurcation analysis reveals critical Turing instabilities and multistability phenomena, offering insights into pattern control and the onset of vegetation band formation.

The Klausmeier-Gray-Scott model is a vegetation-water interaction model for dryland pattern formation in which the vegetation equation retains the autocatalytic uptake term $u^{2}w$, while the water equation uses diffusion rather than the downhill advection term of the classical Klausmeier model. In the flat-terrain formulation summarized in "Fractional Vegetation-Water Model in Arid and Semi-Arid Environments: Pattern Formation and Numerical Simulations" (Speciale et al., 6 Oct 2025), it is the diffusion-dominated endpoint of a continuum connecting slope-driven Klausmeier dynamics to diffusion-driven dynamics. Closely related literature also studies a quasilinear "Gray-Scott-Klausmeier vegetation-water interaction model" with porous-medium water transport and advection, indicating a broader family of Klausmeier/Gray-Scott-type systems used to analyze Turing instability, amplitude equations, and multistability in patterned vegetation (Belin et al., 22 Jan 2026).

1. Canonical formulations and model lineage

In the formulation explicitly identified as the Klausmeier-Gray-Scott model, the plant biomass $u$ and water $w$ satisfy

$$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial^2 w}{\partial x^2} - w - u^2 w + a, \end{cases}$$

where slope is absent and water spreading is via diffusion. The corresponding Klausmeier model replaces the water-diffusion term by advection,

$$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial w}{\partial x} - w - u^2 w + a, \end{cases}$$

with $\nu$ encoding slope as water velocity down the slope (Speciale et al., 6 Oct 2025).

Model	Water transport term	Physical regime
Klausmeier	$\nu \frac{\partial w}{\partial x}$	High slope
Klausmeier-Gray-Scott	$\nu \frac{\partial^2 w}{\partial x^2}$	Flat
Fractional bridge	$\nu D_z^{\alpha+1}W$	Intermediate/anomalous

A related quasilinear formulation studies vegetation density $v$ and water density $u$0 through

$u$1

and is described as a Gray-Scott-Klausmeier model intended to capture vegetation bands on landscapes, especially on slopes (Belin et al., 22 Jan 2026). This suggests that the term "Klausmeier-Gray-Scott model" is often used in practice for a model class rather than a single canonical PDE, with the common structural feature being the Gray-Scott-type autocatalytic interaction between vegetation and water.

2. Spatially homogeneous kinetics and the Gray-Scott bifurcation template

A rigorous bifurcation benchmark for this model class is provided by the spatially homogeneous Gray-Scott reduction

$u$2

obtained by dropping diffusion from the Gray-Scott PDE. This two-dimensional dynamical system has a trivial steady state at $u$3, always stable, and may have two non-trivial critical points whose existence is determined by

$u$4

When $u$5, the two nontrivial points coalesce in a saddle-node bifurcation (Delgado et al., 2017).

The global bifurcation map contains several codimension-one and codimension-two structures. A Takens-Bogdanov point occurs at

$u$6

where the Jacobian has a zero eigenvalue of multiplicity two; locally the dynamics are equivalent to

$u$7

corresponding to the classical subcritical BT scenario. A Bautin point occurs at

$u$8

where the first Lyapunov coefficient vanishes and the Hopf bifurcation changes from supercritical to subcritical. From these organizing centers emanate Hopf, homoclinic, saddle-node, and fold-of-cycles curves, partitioning parameter space into regions with distinct phase portraits, including coexistence of two limit cycles of different stability (Delgado et al., 2017).

For Klausmeier-Gray-Scott research, the importance of this Gray-Scott result is not that the same codimension-two points have been established for the vegetation-water PDE, but that the paper explicitly states that the methods and the bifurcation map also apply to related autocatalytic systems, such as the Klausmeier-Gray-Scott class. The homogeneous Gray-Scott analysis therefore serves as a reference bifurcation template for the local kinetics underlying Klausmeier/Gray-Scott vegetation models.

3. Fractional interpolation between slope-driven and diffusion-driven regimes

A direct bridge between the classical Klausmeier model and the Klausmeier-Gray-Scott model is provided by the fractional traveling-wave system

$u$9

with $w$0, migration speed $w$1, and a space-fractional Caputo derivative of order $w$2: $w$3 The parameter $w$4 is explicitly linked to the slope of the domain. At $w$5, $w$6 and the model recovers classical advection, i.e. the Klausmeier regime. At $w$7, $w$8 and the model recovers classical diffusion, i.e. the Klausmeier-Gray-Scott regime. Intermediate values $w$9 model anomalous advection and nonlocal water redistribution (Speciale et al., 6 Oct 2025).

The linearized Hopf analysis yields a critical migration speed $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial^2 w}{\partial x^2} - w - u^2 w + a, \end{cases}$$0 as a function of $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial^2 w}{\partial x^2} - w - u^2 w + a, \end{cases}$$1. The full formula is given in Eq. (3.8) of the paper, while the limiting cases are

$$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial^2 w}{\partial x^2} - w - u^2 w + a, \end{cases}$$2

Accordingly, the Hopf bifurcation migrates to lower speeds as $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial^2 w}{\partial x^2} - w - u^2 w + a, \end{cases}$$3 increases, mirroring the suppression of pattern migration as terrain flattens. Numerical simulations using an explicit method adapted to the Caputo derivative show oscillatory vegetation patterns for $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial^2 w}{\partial x^2} - w - u^2 w + a, \end{cases}$$4 across a range of $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial^2 w}{\partial x^2} - w - u^2 w + a, \end{cases}$$5: for $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial^2 w}{\partial x^2} - w - u^2 w + a, \end{cases}$$6 the bands migrate uphill, for intermediate $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial^2 w}{\partial x^2} - w - u^2 w + a, \end{cases}$$7 migration persists but slows, and for $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial^2 w}{\partial x^2} - w - u^2 w + a, \end{cases}$$8 the migration speed vanishes and stationary periodic patterns emerge (Speciale et al., 6 Oct 2025).

A common simplification is to regard the Klausmeier and Klausmeier-Gray-Scott models as qualitatively separate systems. The fractional formulation shows instead that, in the summarized literature, they can be treated as endpoints of a continuous interpolation in which the terrain-coupled parameter $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial^2 w}{\partial x^2} - w - u^2 w + a, \end{cases}$$9 controls the transition from advection-dominated to diffusion-dominated pattern dynamics.

4. Pattern onset and Ginzburg-Landau reduction

Near the first instability, the Gray-Scott-Klausmeier vegetation-water model admits a rigorous weakly nonlinear reduction to the Ginzburg-Landau equation. The full quasilinear system

$$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial w}{\partial x} - w - u^2 w + a, \end{cases}$$0

has homogeneous steady states, and linear stability analysis around non-trivial equilibria reveals conditions for Turing and Turing-Hopf instabilities. For suitable parameters, lowering rainfall $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial w}{\partial x} - w - u^2 w + a, \end{cases}$$1 can induce a Turing instability at critical wavenumber $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial w}{\partial x} - w - u^2 w + a, \end{cases}$$2, marking the onset of pattern formation (Belin et al., 22 Jan 2026).

The slow modulation of the critical mode is governed by the amplitude equation

$$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial w}{\partial x} - w - u^2 w + a, \end{cases}$$3

where $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial w}{\partial x} - w - u^2 w + a, \end{cases}$$4, $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial w}{\partial x} - w - u^2 w + a, \end{cases}$$5, and the coefficients $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial w}{\partial x} - w - u^2 w + a, \end{cases}$$6 are computable from the original model parameters and the linearized operator. The central theorem states that, for sufficiently smooth initial data close to the steady state and for $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial w}{\partial x} - w - u^2 w + a, \end{cases}$$7 small enough, there exists a solution of the full Gray-Scott-Klausmeier system satisfying

$$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial w}{\partial x} - w - u^2 w + a, \end{cases}$$8

so the Ginzburg-Landau approximation is valid on timescales $$\begin{cases} \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} - m u + u^2 w,\ \frac{\partial w}{\partial t} = \nu \frac{\partial w}{\partial x} - w - u^2 w + a, \end{cases}$$9 with error $\nu$0 (Belin et al., 22 Jan 2026).

This result fixes the scope of amplitude-equation reasoning in the Klausmeier-Gray-Scott setting. The Ginzburg-Landau equation is justified close to the Turing instability and for initial conditions sufficiently close to the homogeneous state; it is not presented as a global reduction valid arbitrarily far from onset. Within that regime, however, it provides a mathematically rigorous reduced description of band selection, slow envelope dynamics, and weakly nonlinear saturation even in the presence of the quasilinear porous-medium term $\nu$1.

5. Nonlocal diffusion, weak-solution theory, and analytical obstacles

Nonlocal variants of Klausmeier-type vegetation models replace classical plant diffusion by an integro-differential operator,

$\nu$2

with $\nu$3 symmetric, positive, radial, $\nu$4, and finite second moment. On a bounded domain $\nu$5, the resulting system is

$\nu$6

with nonlocal Dirichlet volume constraints $\nu$7 for $\nu$8 and homogeneous Dirichlet boundary conditions for the water equation (Jaramillo et al., 2024).

The mathematical distinction from classical diffusion is decisive. The Laplacian generates a bilinear form on $\nu$9, whereas the nonlocal operator defines a bilinear form only on a subspace of $\nu \frac{\partial w}{\partial x}$0,

$\nu \frac{\partial w}{\partial x}$1

Because the Galerkin approximations for $\nu \frac{\partial w}{\partial x}$2 are only inherently bounded in $\nu \frac{\partial w}{\partial x}$3 rather than $\nu \frac{\partial w}{\partial x}$4, Aubin’s compactness theorem does not apply directly. The existence proof therefore augments the system with equations for the spatial derivatives $\nu \frac{\partial w}{\partial x}$5 and $\nu \frac{\partial w}{\partial x}$6, derives energy estimates for $\nu \frac{\partial w}{\partial x}$7, and uses the additional regularity to recover compactness and pass to the limit in the nonlinear terms (Jaramillo et al., 2024).

For the Klausmeier-Gray-Scott literature, the significance is methodological. The paper places these techniques in the broader context of other nonlocal reaction-diffusion systems, such as the Gray-Scott model and the Klausmeier-Gray-Scott model, where long-range dispersal or nonlocal transport destroys the standard $\nu \frac{\partial w}{\partial x}$8-based parabolic compactness framework. The need for volume constraints, augmented derivative equations, and carefully chosen bilinear forms is therefore part of the analytical infrastructure of generalized Klausmeier-Gray-Scott theory.

6. Multistability, Busse balloons, and non-local grazing

The Busse balloon is the two-dimensional region in parameter-wavenumber space within which spatially periodic patterns are stable. In dryland vegetation models, it encodes the coexistence of multiple stable wavelengths for the same precipitation level. For a Klausmeier model with non-local grazing, recent work extends singular perturbation methods previously used for Klausmeier/Gray-Scott models to include non-local interaction, obtains analytical control near the homoclinic limit, and combines this with numerical continuation to study the deformation of the Busse balloon away from the typical banana shape (Siero et al., 19 Jun 2026).

The model considered has the form

$\nu \frac{\partial w}{\partial x}$9

with grazing pressure depending on spatial averages $\nu \frac{\partial^2 w}{\partial x^2}$0. Four regimes are distinguished: proportional/disproportionate and sustained/natural grazing. The Busse-balloon response depends strongly on regime. Under proportional sustained grazing, the balloon lifts off from the $\nu \frac{\partial^2 w}{\partial x^2}$1 axis at low precipitation and low wavenumber; under proportional natural grazing, it becomes pinched at intermediate wavenumbers and precipitation; under disproportionate grazing, it can split or disconnect into two regions (Siero et al., 19 Jun 2026).

These results modify standard intuition about resilience in patterned drylands. The paper emphasizes that large-wavelength patterns do not persist the longest under decreasing precipitation in the sustained grazing scenarios, and it underscores the importance of Busse-balloon shape when inferring ecosystem response to climate change. For the Klausmeier-Gray-Scott framework, the main lesson is that multistability is not adequately summarized by the mere existence of a Turing threshold; the geometry of the stable-wavelength set is itself a dynamical observable, and non-local interaction can reorganize it qualitatively.

7. Control and rigorous pattern construction in related Gray-Scott systems

Two additional Gray-Scott results delimit techniques that are plausibly transferable to Klausmeier-Gray-Scott systems. First, time-delayed feedback control in the Gray-Scott reaction-diffusion model uses

$\nu \frac{\partial^2 w}{\partial x^2}$2

with single-species, diagonal, and mixed control schemes. In the summarized results, activator control with $\nu \frac{\partial^2 w}{\partial x^2}$3 or inhibitor control with $\nu \frac{\partial^2 w}{\partial x^2}$4 stabilizes the non-trivial homogeneous steady state or converts traveling waves to stationary Turing patterns, whereas diagonal control fails to control spatiotemporal chaos and is ineffective for traveling waves for any $\nu \frac{\partial^2 w}{\partial x^2}$5 and $\nu \frac{\partial^2 w}{\partial x^2}$6. The paper explicitly states that these results are generally applicable to other spatially extended systems with similar kinetic architecture, such as the related Klausmeier-Gray-Scott models (Kyrychko et al., 2012).

Second, a constructive Newton-Kantorovich framework has been developed for the stationary 2D Gray-Scott equations in a Hilbert space $\nu \frac{\partial^2 w}{\partial x^2}$7 with symmetry-reduced subspace $\nu \frac{\partial^2 w}{\partial x^2}$8. With an approximate inverse $\nu \frac{\partial^2 w}{\partial x^2}$9 for the linearization, existence and uniqueness of a localized stationary solution follow from radii-polynomial conditions involving bounds $\nu D_z^{\alpha+1}W$0, $\nu D_z^{\alpha+1}W$1, and $\nu D_z^{\alpha+1}W$2. This framework yields computer-assisted proofs of four different localized $\nu D_z^{\alpha+1}W$3-symmetric stationary patterns and an additional theorem showing that such localized structures arise as limits of periodic branches as the period tends to infinity (Cadiot et al., 2024).

Neither result is a direct theorem about the Klausmeier-Gray-Scott PDE itself, but together they delineate an active methodological perimeter around the model class: delayed-feedback control for steering spatiotemporal regimes, and computer-assisted functional-analytic verification of localized patterns. Since the constructive framework is stated to extend readily to a large class of reaction-diffusion systems, a plausible implication is that analogous proof architectures may become available for stationary Klausmeier-Gray-Scott patterns once the relevant invertibility and product estimates are established.