Pith. sign in

REVIEW 2 major objections 4 minor 68 references

A flow-kick model of dryland vegetation patterns: the impact of rainfall variability on resilience

T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Rainfall randomness narrows the window for banded vegetation patterns.

desk verdict Solid stochastic extension of the flow-kick dryland model; direction of the result is convincing, but the headline MAP=34.8 threshold rests on an unvalidated Lyapunov-exponent estimate. read the letter →

arxiv 2501.01569 v2 pith:NA5O5RJI submitted 2025-01-02 nlin.PS

Punit Gandhi , Matthew Oline , Mary Silber This is my paper · ORCID
classification nlin.PS MSC 34A3735K5737N2592D40
keywords drylandvegetationpatternsflow-kickmodelrainfallvariabilityrandomdynamicalsystemsLyapunovexponentsreaction-diffusionresiliencebanded
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Dryland vegetation self-organizes into bands that harvest runoff from bare upslope areas, but whether this survival mechanism works may depend on rainfall variability, not just total rainfall. This paper compares a model with perfectly periodic storms to a model with randomly timed and randomly sized storms at the same mean annual precipitation. Under random rainfall the uniform vegetation state loses stability to banded patterns only at a lower mean precipitation--$\mathrm{MAP}=34.8$ cm/year in both protocols studied, versus $52.4$ or $75.4$ cm/year for periodic rainfall--and the collapse boundary for uniform vegetation moves upward. Numerical simulations show the whole interval of mean rainfall over which patterns exist shrinks, and time-averaged biomass drops. The conclusion is that increasing storm variability, even without changing mean rainfall, can reduce the resilience that banded vegetation patterns provide.

What carries the argument

The load-bearing object is the kick-flow stability map $M_k[w_0,b_0,h_0,\tau_d]=\Psi_k\,\Omega_k$, whose product over successive random storms gives the maximal Lyapunov exponent $\lambda_k$ from Eq. (20). The kick factor $\Omega_k$ contains the positive feedbacks that concentrate storm water: infiltration-enhancement and flow-speed-reduction terms, with the overland travel distance $\ell_0=\nu(b)h_0/\iota(b)$ setting the resonance wavelengths $k_n^*=n\pi/\ell_0$. The flow factor $\Psi_k$ integrates the slow reaction-diffusion system linearized about the post-kick uniform state. The paper uses the sign of $\lambda_k$ as the pattern-forming instability criterion, and in the periodic case it reduces to Floquet theory.

What would settle it

Run both ramped simulation protocols with a different long random sequence (for example $10^6$ cycles instead of $10^5$) and check whether the zero crossing of the maximal Lyapunov exponent still lands at $\mathrm{MAP}=34.8$ cm/year in both the $H_0$-fixed and $T_d$-fixed paths; if the onset shifts by more than a few cm/year, the fixed-onset claim is falsified.

Watch

Extended reading notes

Core claim

The paper's central claim is that random variability in storm timing and storm depth, at the same mean annual precipitation as a periodic reference, shrinks the region in rainfall-parameter space where dryland vegetation patterns can establish and persist. Under periodic rainfall the uniform vegetation state loses stability to banded patterns at $\mathrm{MAP}=52.4$ cm/year when $H_0=1$ cm and at $75.4$ cm/year when $T_d=15$ days; under exponentially distributed storms and interstorm times the same instability occurs at $\mathrm{MAP}=34.8$ cm/year in both protocols, and the lower transcritical threshold below which the uniform state collapses moves to higher precipitation. Simulations confirm the narrowed interval and show that random rainfall lowers time-averaged biomass, destroys the resonance-tongue structure of the linear instability, and moves pattern collapse to higher mean rainfall, so variability acts as an extra stressor even when mean rainfall is fixed.

Load-bearing premise

The load-bearing premise is that the numerically estimated maximal Lyapunov exponent of the random storm sequence--computed from one sequence of $10^5$ cycles with rescaling every $100$ cycles--marks the true onset of pattern formation; the paper itself notes that negative Lyapunov exponents do not always imply stability, and the check covers only two paths in the $(H_0,T_d)$ plane.

Editorial extensions

If this is right

  • Under random rainfall the uniform vegetation state is stable over a wider range of low and high mean precipitation than under periodic rainfall, so patterns appear only in a narrower band of $\mathrm{MAP}$.
  • The linear instability's resonance tongues, tied to the distance surface water travels, disappear with randomness; simulations then show simple traveling bands migrating uphill.
  • The lower collapse boundary (transcritical point) rises, and the bistable region where patterns and bare soil coexist shrinks; patterns persist below it only for a finite stochastic lifetime.
  • The random pattern-onset boundary is insensitive to whether $\mathrm{MAP}$ is set by changing storm depth or interstorm time, sitting near $\mathrm{MAP}=34.8$ cm/year.
  • Random rainfall lowers the time-averaged biomass of both uniform and patterned states compared with periodic rainfall at the same mean annual precipitation.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the model draws storm depths and waiting times from exponential distributions, the near-universal onset at $\mathrm{MAP}=34.8$ cm/year may be specific to those distributions; testing gamma-distributed storms or correlated storm clusters would show whether the threshold depends on the shape of the variability.
  • The finding that linear predictions match nonlinear simulations only when randomness is added suggests a time-averaged or effective-parameter model might reproduce the random behavior with much cheaper computation--an avenue the paper does not pursue.
  • If this result carries to real hillslopes, monitoring changes in rainfall variability, not just mean annual precipitation, could help forecast where banded vegetation will disappear; remote sensing of band widths and migration speeds might serve as an early indicator.
  • The paper's periodic-boundary setup omits water loss through runoff; including open boundaries could further shrink or shift the random-rainfall window, likely strengthening the reported effect.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The paper studies a flow-kick model of banded vegetation on gentle hillslopes, in which rainstorms are instantaneous kicks to soil water and the dry intervals are governed by a two-component reaction-diffusion system. It compares idealized periodic rainfall with random rainfall (exponentially distributed storm depths and inter-storm times) at the same mean annual precipitation. The main claims are that randomness shifts the transcritical collapse threshold upward, lowers the MAP at which the uniform vegetation state becomes unstable to spatial perturbations, and narrows the MAP interval over which patterns exist, thereby reducing the resilience of pattern-forming dryland ecosystems. Quantitative anchors are the random-rainfall pattern onset at MAP = 34.8 cm/year versus periodic onsets of 52.4 cm/year (H0 = 1 cm) and 75.4 cm/year (Td = 15 days) along the two parameter paths studied.

Significance. The question addressed is timely and the flow-kick framework is well matched to it. If the quantitative thresholds hold, the paper delivers a clear, falsifiable prediction: increased rainfall variability with unchanged mean precipitation reduces the resilience of banded vegetation. Strengths of the manuscript include an exact analytic transcritical threshold for random kick strength with a proved q-Pochhammer identity (Appendix B), a transparent Lyapunov-exponent formulation for the random kick-flow map (Section 4.2), and an extensive battery of ramped simulations with an explicit state-classification protocol (Section 5.2). The authors also honestly flag the limitation that negative Lyapunov exponents do not in general imply stability (Section 6, citing Leonov and Kuznetsov). The main quantitative claim, however, currently rests on a single-sequence Lyapunov estimator and on only two validated simulation paths, so the headline number 34.8 cm/year and the associated invariance statement need additional support before the strongest version of the resilience conclusion can be accepted.

major comments (2)
  1. [Section 4.2, Eq. (21)] The maximal Lyapunov exponent for the random kick-flow map is approximated from a single sequence of n = 10^5 cycles with rescaling every m = 100 iterations, and the paper reports no convergence study in n or m, no seed-to-seed variability, and no confidence interval. Because the zero crossing of this estimator sets the random-rainfall pattern onset at MAP = 34.8 cm/year (Figures 6 and 8), a finite-time bias or slow convergence in the estimator would directly move the paper's headline number. I recommend adding a convergence study (e.g., the inferred onset as a function of n for n = 10^4, 10^5, 10^6 and of the rescaling interval m), repeating the calculation over several independent random sequences, and reporting a bootstrap or inter-quantile range for the inferred MAP. The uncertainty should be stated next to the 34.8 cm/year value.
  2. [Section 5.1 and Figure 8] The statement that with random rainfall 'the onset of the pattern forming instability always occurs at this MAP, independent of the particular choice of H0 and Td associated with it' is an extrapolation from two one-parameter paths (H0 = 1 cm with Td varied, and Td = 15 days with H0 varied) together with the threshold curve in Figure 7(a), which is computed with the same Lyapunov estimator. Two validated paths do not establish invariance across the (H0, Td) plane, and Section 6 itself notes that negative Lyapunov exponents do not imply stability, so the linear zero crossing is not by itself a rigorous onset. I recommend either testing additional (H0, Td) paths or explicitly softening the claim to the tested paths and adjusting the abstract and conclusion to state the uncertainty in the universal value.
minor comments (4)
  1. [Throughout] There are several typographical errors: 'Lyaponuv' should be 'Lyapunov' (Sections 4.2 and 6), 'trails' should be 'trials' (Section 2.5 and Figure 12), 'obeserved' should be 'observed' (Section 5.2), 'purpler dotted' should be 'purple dotted' (Section 4.3), and 'initialed' should be 'initialized' (Section 2.5).
  2. [Figure 8(c)] The text describing Figure 8(c) refers to the second distribution as 'magenta', while the caption says 'purple'; please make the color names consistent.
  3. [Section 2.4 and Section 3] The notation h0 and τd is used both for dimensionless variables and for dimensioned quantities, and the formula MAP = 3.65 h0/τd mixes these conventions. Clarify the units in the sentence preceding Eq. (3) and in the dimensionless map (13).
  4. [Section 3, Eq. (17)] For the fully random case and for random-timing-only case, the expectation in Eq. (17) is evaluated numerically over 10^6 cycles; reporting the sensitivity of the inferred threshold to the number of cycles would make the lower-boundary results easier to assess.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the random-versus-periodic rainfall comparison is a new computation that does not reduce to the model inputs or to prior results.

full rationale

Walked the derivation chain. Section 3 derives the transcritical threshold from the linearized uniform flow (Eqs. 14-17) and validates it against direct flow-kick cycle averages; Section 4 constructs the kick-flow stability map and computes maximal Lyapunov exponents for both periodic and random rainfall; Section 5 checks the predicted onset and wavenumbers with ramped nonlinear simulations. No predicted quantity is defined in terms of the output it claims to explain. The random-rainfall instability onset at MAP~34.8 cm/year is a computed zero crossing of the Lyapunov-exponent estimator, not a fitted parameter, and the claimed H0/Td independence is an observed alignment of the computed threshold with a constant-MAP line rather than an imposed constraint. The only notable self-citation is use of the kick sensitivity formula J_k and the resonance/wavelength picture from Gandhi et al. [26]; that prior work derives the formula from the stated hydrological PDE and does not contain the random-rainfall result, so it is independent support rather than load-bearing circularity. The Section 6 caveat about negative Lyapunov exponents and Perron effects, together with the two-path simulation validation, bears on numerical reliability and extrapolation strength, not on whether the derivation is circular.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The model relies on the flow-kick framework of Gandhi et al. [26] with exponential rainfall statistics. No new entities are introduced. The free parameters listed are numerical protocol choices, not fitted physical constants; the central claim is robust to their exact values, but the quoted MAP thresholds may shift with different choices.

free parameters (3)
  • State classification threshold B_epsilon = 0.02 kg/m^2
    Hand-chosen threshold in the ramped simulations (Section 5.2, Fig. 9) for distinguishing spatially uniform from patterned states; the shaded MAP intervals in Fig. 8 could shift if this value were different.
  • Lyapunov exponent estimation length n and rescaling interval m = n = 10^5, m = 10^2
    Numerical hyperparameters chosen for the estimate in Eq. (21); the onset MAP = 34.8 cm/year is quoted without confidence intervals.
  • Ramped simulation step size and dwell time = 0.1 cm/year every 25 to 100 years
    Choice of ramp protocol in Section 5.2; the paper shows collapse ranges are roughly independent of ramp rate (Fig. 14), but pattern onset intervals in Fig. 12 may exhibit hysteresis.
assumptions (4)
  • domain assumption Storm arrival times follow a Poisson process and storm depths are exponentially distributed (Section 2.3)
    This is a simplified rainfall model from the hydrology literature; the quantitative size of the variability effect may depend on this distributional choice.
  • domain assumption The kick is computed as the quasi-steady state of the surface water PDE (2) with infiltration rate and flow speed depending on biomass via (8), and with flow speed independent of water depth
    This fast hydrology idealization inherited from Gandhi et al. [26] determines the spatial profile of the soil water kick, which is the only water transport in the model.
  • domain assumption The linearized kick contribution J_k in Eq. (19) is taken from Gandhi et al. [26] without re-derivation
    The paper's Lyapunov exponent analysis uses this formula as a building block; it is a previously published result by overlapping authors.
  • domain assumption Periodic boundary conditions are used on a 1D domain representing a mid-slope section (Section 2.5 and limitations in Section 6)
    This neglects water loss through runoff at the domain edges, which the authors acknowledge may be important for high-intensity storms.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A flow-kick model of dryland vegetation patterns: the impact of rainfall variability on resilience." pith.science (2026). https://pith.science/paper/NA5O5RJI

@misc{pith2026250101569,
  author       = {Pith},
  title        = {Pith review of: A flow-kick model of dryland vegetation patterns: the impact of rainfall variability on resilience},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NA5O5RJI}},
  note         = {Machine review of arXiv:2501.01569}
}
read the original abstract

In many drylands around the globe, vegetation self-organizes into regular spatial patterns in response to aridity stress. We consider the regularly-spaced vegetation bands, on gentle hill-slopes, that survive low rainfall conditions by harvesting additional stormwater from upslope low-infiltration bare zones. We are interested in the robustness of this pattern formation survival mechanism to changes in rainfall variability. For this, we use a flow-kick modeling framework that treats storms as instantaneous kicks to the soil water. The positive feedbacks in the storm-level hydrology, that act to concentrate water within the vegetation bands, are captured through the spatial profiles of the soil water kicks. Between storms, the soil water and vegetation, modeled by a two-component reaction-diffusion system, evolve together. We use a combination of linear stability analysis and numerical simulation to compare predictions of idealized periodic rainfall, with no variability, to predictions when there is randomness in the timing and magnitude of water input from storms. We show that including these random elements leads to a decrease in the parameter range over which patterns appear. This suggests that an increase in storm variability, even with the same mean annual rainfall, may negatively impact the resilience of these pattern-forming dryland ecosystems.

Figures

Figures reproduced from arXiv: 2501.01569 by the authors.

Figure 1
Figure 1. (a) Satellite images of banded vegetation patterns from sites in Mexico (27 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. (a) Distribution of storms as a function of mean depth [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Two examples of a sequence of surface water kicks [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figures from the paper (14 more)
Figure 4
Figure 4. Figure 4: (a) Plot of biomass level B0 at the start of a flow-kick cycle vs. mean annual precipitation rate (MAP). The blue solid curve is for the case of periodic kicks with H0 = 5 cm, and incrementing the time between kicks to decrease MAP. There is a transcritical bifurcation…
Figure 5
Figure 5. Figure 5: (a) Zero-biomass stability boundaries for periodic kick, random kick strength only, random [PITH_FULL_IMAGE:figures/full_fig_p013_5.png]
Figure 1
Figure 1. Figure 1: Random Rainfall With random rainfall the largest Lyaponuv exponent λk, given by (20), can be numerically approximated by computing perturbation amplitudes after n repeated iterations 16 [PITH_FULL_IMAGE:figures/full_fig_p016_1.png]
Figure 6
Figure 6. Figure 6: Maximal Lyaponuv exponent λk given by equation (20) as a function of MAP and (dimensioned) K. We control MAP by (a) fixing mean storm depth at H0 = 1 cm and varying the mean time between storms Td, and by (b) fixing the mean time between storms to be Td = 15 days and v…
Figure 7
Figure 7. Figure 7: (a) Thresholds for pattern forming instability in the flow-kick parameter space under both [PITH_FULL_IMAGE:figures/full_fig_p018_7.png]
Figure 8
Figure 8. Figure 8: Mean biomass level prior to storm for uniform states with random rainfall is shown in [PITH_FULL_IMAGE:figures/full_fig_p020_8.png]
Figure 9
Figure 9. Figure 9: (a) Schematic of simulation protocol for ramped numerical experiments. (b) Decision [PITH_FULL_IMAGE:figures/full_fig_p021_9.png]
Figure 10
Figure 10. Figure 10: Ramped Simulations with Periodic Rainfall. Biomass levels (upper panels) and mean wavenumbers (lower panels) for patterns from simulations in which MAP is decreased (blue) and increased (red) by (a) changing the dry period Td between storms with fixed H0 = 1 cm and by…
Figure 11
Figure 11. Figure 11: Spacetime plots of annually-averaged biomass distribution from 200 year simulations on [PITH_FULL_IMAGE:figures/full_fig_p024_11.png]
Figure 12
Figure 12. Figure 12: Ramped simulations near upper instability with random rainfall. Simulations results from 50 trials in which MAP is decreased (blue) and increased (red) near the pattern forming instability (vertical black dotted line). MAP is varied by (a) changing the dry period Td b…
Figure 13
Figure 13. Figure 13: Ramped simulations near lower instability with random rainfall. Simulation results from 50 trials in which MAP is decreased (blue) and increased (red) near the transcritical point (vertical black dotted line). MAP is varied by (a) changing the dry period Td between st…
Figure 14
Figure 14. Figure 14: (a) Fraction of simulations with patterns in simulations where [PITH_FULL_IMAGE:figures/full_fig_p028_14.png]
Figure 15
Figure 15. Figure 15: Computing the kick ∆w(x) for storm depth h0 = 1 and ˆb(x) = 2(1 + cos(4x)), indicated respectively by the blue line and green shading in the top panel. (top) Characteristics hˆ(x; y) from (29). For x ∈ [x ∗ 1 , x∗ 2 ] we see there is only one piece where hˆ ≤ h0, so t…
Figure 16
Figure 16. Figure 16: (a) Comparison of (36) with the distribution of [PITH_FULL_IMAGE:figures/full_fig_p036_16.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

68 extracted references · 68 canonical work pages

  1. [1]

    F. T. Maestre, B. M. Benito, M. Berdugo, L. Concostrina-Zubiri, M. Delgado-Baquerizo, D. J. Eldridge, E. Guirado, N. Gross, S. K´ efi, Y. Le Bagousse-Pinguet, R. Ochoa-Hueso, and Soliveres. S. Biogeography of global drylands. New Phytologist, 231(2):540–558, 2021

    work page 2021

  2. [2]

    F. T. Maestre, D. J. Eldridge, S. Soliveres, S. K´ efi, M. Delgado-Baquerizo, M. A. Bowker, P. Garc ´ ıa-Palacios, J. Gait´ an, A. Gallardo, R. L´ azaro, and M. Berdugo. Structure and func- tioning of dryland ecosystems in a changing world.Annu. Rev. Ecol. Evol. Syst., 47(1):215–237, 2016

    work page 2016

  3. [3]

    Menaut and B

    J.-C. Menaut and B. Walker. Banded vegetation patterning in arid and semiarid environ- ments: ecological processes and consequences for management , volume 149. Springer Science & Business Media, 2001

    work page 2001

  4. [4]

    Deblauwe, N

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

    work page 2008

  5. [5]

    Gandhi, S

    P. Gandhi, S. Iams, S. Bonetti, and M. Silber. Vegetation pattern formation in drylands. In Dryland Ecohydrology, pages 469–509. Springer, 2019

    work page 2019

  6. [6]

    W. A. Macfadyen. Vegetation patterns in the semi-desert plains of British Somaliland. Geogr. J., 116(4/6):199–211, 1950

    work page 1950

  7. [7]

    J. E. G. W. Greenwood. The development of vegetation patterns in somaliland protectorate. Geogr. J., pages 465–473, 1957

    work page 1957

  8. [8]

    C. F. Hemming. Vegetation arcs in somaliland. J. Ecol., pages 57–67, 1965

    work page 1965

Show all 68 references
  1. [9]

    Deblauwe, P

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

    2012

  2. [10]

    G. G. Penny, K. E. Daniels, and S. E. Thompson. Local properties of patterned vegetation: quantifying endogenous and exogenous effects. Phil. Trans. R. Soc. A , 371(2004):20120359, 2013

    2004

  3. [11]

    Gowda, S

    K. Gowda, S. Iams, and M. Silber. Signatures of human impact on self-organized vegetation in the Horn of Africa. Sci. Rep., 8:3622, 2018

    2018

  4. [12]

    Bastiaansen, O

    R. Bastiaansen, O. Ja ¨ ıbi, V. Deblauwe, M. B. Eppinga, K. Siteur, E. Siero, S. Mermoz, A. Bou- vet, A. Doelman, and M. Rietkerk. Multistability of model and real dryland ecosystems through spatial self-organization. Proc. Natl. Acad. Sci. , 115(44):11256–11261, 2018

    2018

  5. [13]

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

    1999

  6. [14]

    Lejeune, P

    O. Lejeune, P. Couteron, and R. Lefever. Short range co-operativity competing with long range inhibition explains vegetation patterns. Acta Oecologica, 20(3):171–183, 1999

    1999

  7. [15]

    Rietkerk, M

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

    2002

  8. [16]

    Gilad, J

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

    2004

  9. [17]

    Bastiaansen, A

    R. Bastiaansen, A. Doelman, M. B. Eppinga, and M. Rietkerk. The effect of climate change on the resilience of ecosystems with adaptive spatial pattern formation.Ecol. Lett., 23(3):414–429, 2020

    2020

  10. [18]

    Rietkerk, S

    M. Rietkerk, S. C. Dekker, P. C. de Ruiter, and J. van de Koppel. Self-organized patchiness and catastrophic shifts in ecosystems. Science, 305(5692):1926–1929, 2004

    1926

  11. [19]

    E. Meron. Nonlinear physics of ecosystems . CRC Press Boca Raton, FL, 2015

    2015

  12. [20]

    Martinez-Garcia, C

    R. Martinez-Garcia, C. Cabal, J. M. Calabrese, E. Hern´ andez-Garc ´ ıa, C. E. Tarnita, C. L´ opez, and J. A. Bonachela. Integrating theory and experiments to link local mechanisms and ecosystem-level consequences of vegetation patterns in drylands. Chaos, Solitons & Fractals ...

    2023

  13. [21]

    Noy-Meir

    I. Noy-Meir. Desert ecosystems: environment and producers. Annu. Rev. Ecol. Evol. Syst. , 4(1):25–51, 1973

    1973

  14. [22]

    Airbus, CNES/Airbus, Landsat/Copernicus, Maxar Technologies, USDA/FPAC/GEO

    Google Map. Airbus, CNES/Airbus, Landsat/Copernicus, Maxar Technologies, USDA/FPAC/GEO. Accessed Dec. 30, 2024

    2024

  15. [23]

    Cooperative observer network (COOP) data

    National Oceanic and Atmospheric Administration (NOAA). Cooperative observer network (COOP) data. Accessed via NOAA National Centers for Environmental Information website, 2024

    2024

  16. [24]

    Masson-Delmotte, P

    V. Masson-Delmotte, P. Zhai, A. Pirani, S. L. Connors, C. C. B. Pell, A. G. Hope, A. C. Allen, M. Tignor, and E. Poloczanska. Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Clima...

    2021

  17. [25]

    F. T. Maestre, R. Salguero-G´ omez, and J. L. Quero. It is getting hotter in here: determining and projecting the impacts of global environmental change on drylands, 2012

    2012

  18. [26]

    Gandhi, L

    P. Gandhi, L. Liu, and M. Silber. A pulsed-precipitation model of dryland vegetation pattern formation. SIAM J. Appl. Dyn. Syst. , 22:657, 2023

    2023

  19. [27]

    von Hardenberg, E

    J. von Hardenberg, E. Meron, M. Shachak, and Y. Zarmi. Diversity of vegetation patterns and desertification. Phys. Rev. Lett. , 87(19):198101, 2001

    2001

  20. [28]

    Siteur, E

    K. Siteur, E. Siero, M. B Eppinga, J. D. M. Rademacher, A. Doelman, and M. Rietkerk. Beyond Turing: The response of patterned ecosystems to environmental change. Ecol. Complex. , 20:81–96, 2014

    2014

  21. [29]

    Gowda, Y

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

    2016

  22. [30]

    Meyer, A

    K. Meyer, A. Hoyer-Leitzel, S. Iams, I. Klasky, V. Lee, S. Ligtenberg, E. Bussmann, and M. L. Zeeman. Quantifying resilience to recurrent ecosystem disturbances using flow–kick dynamics. Nature Sustain., 1(11):671–678, 2018

    2018

  23. [31]

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

    2005

  24. [32]

    van der Stelt, A

    S. van der Stelt, A. Doelman, G. Hek, and J. D. M. Rademacher. Rise and fall of periodic patterns for a generalized Klausmeier–Gray–Scott model. J. Nonlin. Sci. , 23(1):39–95, 2013

    2013

  25. [33]

    Striped pattern selection by advective reaction-diffusion systems: Resilience of banded vege- tation on slopes

    E Siero, Arjen Doelman, MB Eppinga, Jens DM Rademacher, M Rietkerk, and K Siteur. Striped pattern selection by advective reaction-diffusion systems: Resilience of banded vege- tation on slopes. Chaos, 25(3):036411, 2015

    2015

  26. [34]

    Carter and A

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

    2018

  27. [35]

    D’Odorico, F

    P. D’Odorico, F. Laio, and L. Ridolfi. Vegetation patterns induced by random climate fluctu- ations. Geophys. Res. Lett., 33(19), 2006

    2006

  28. [36]

    Ursino and S

    N. Ursino and S. Contarini. Stability of banded vegetation patterns under seasonal rainfall and limited soil moisture storage capacity. Adv. Water Resour. , 29(10):1556–1564, 2006

    2006

  29. [37]

    Guttal and C

    V. Guttal and C. Jayaprakash. Self-organization and productivity in semi-arid ecosystems: Implications of seasonality in rainfall. J. Theor. Biol. , 248(3):490–500, 2007

    2007

  30. [38]

    A. G. Konings, S. C. Dekker, M. Rietkerk, and G. G. Katul. Drought sensitivity of patterned vegetation determined by rainfall-land surface feedbacks. J. Geophys. Res. G: Biogeosciences, 116(G4), 2011

    2011

  31. [39]

    Siteur, M

    K. Siteur, M. B. Eppinga, D. Karssenberg, M. Baudena, M. F. P. Bierkens, and M. Rietkerk. How will increases in rainfall intensity affect semiarid ecosystems? Water Resour. Res. , 50(7):5980–6001, 2014

    2014

  32. [40]

    Gandhi, S

    P. Gandhi, S. Bonetti, S. Iams, A. Porporato, and M. Silber. A fast–slow model of banded vegetation pattern formation in drylands. Physica D, 410:132534, 2020

    2020

  33. [41]

    L. F. Gordillo and P. E. Greenwood. Intermittent precipitation-dependent interactions, en- compassing allee effect, may yield vegetation patterns in a transitional parameter range. Bull. Math. Biol. , 85(10):86, 2023

    2023

  34. [42]

    O. V. Crompton and S. E. Thompson. Sensitivity of dryland vegetation patterns to storm characteristics. Ecohydrology, 14(2):e2269, 2021

    2021

  35. [43]

    Crompton, A

    O. Crompton, A. Sytsma, and S. Thompson. Emulation of the saint venant equations enables rapid and accurate predictions of infiltration and overland flow velocity on spatially heteroge- neous surfaces. Water Resour. Res. , 55(8):7108–7129, 2019

    2019

  36. [44]

    A. Y. Kletter, J. Von Hardenberg, E. Meron, and A. Provenzale. Patterned vegetation and rainfall intermittency. J. Theor. Biol. , 256(4):574–583, 2009

    2009

  37. [45]

    Yizhaq, S

    H. Yizhaq, S. Sela, T. Svoray, S. Assouline, and G. Bel. Effects of heterogeneous soil-water diffusivity on vegetation pattern formation. Water Resour. Res. , 50(7):5743–5758, 2014

    2014

  38. [46]

    Yizhaq, I

    H. Yizhaq, I. Stavi, Mo. Shachak, and G. Bel. Geodiversity increases ecosystem durability to prolonged droughts. Ecol. Complex., 31:96–103, 2017

    2017

  39. [47]

    Hamster, P

    C. Hamster, P. van Heijster, and E. Siero. Blurring the busse balloon: Patterns in a stochastic klausmeier model. arXiv preprint arXiv:2411.13238 , 2024. 39

    2024 arXiv

  40. [48]

    Rodriguez-Iturbe, D

    I. Rodriguez-Iturbe, D. R. Cox, and V. Isham. Some models for rainfall based on stochastic point processes. Proc. R. Soc. Lond. A: Math. Phys. Sci. , 410(1839):269–288, 1987

    1987

  41. [49]

    L. F. Shampine and M. W. Reichelt. The Matlab ODE suite. SIAM J. Sci. Comput. , 18(1):1– 22, 1997

    1997

  42. [50]

    J. D. Meiss. Differential Dynamical Systems . SIAM, 2007

    2007

  43. [51]

    K. J. Meyer, H. Fusco, C. Smith, and A. Hoyer-Leitzel. Continuation of fixed points and bifurcations from ode to flow-kick disturbance models. SIAM J. Appl. Dyn. Syst. , 23(4):2983– 3012, 2024

    2024

  44. [52]

    Rietkerk, R

    M. Rietkerk, R. Bastiaansen, S. Banerjee, J. van de Koppel, M. Baudena, and A. Doel- man. Evasion of tipping in complex systems through spatial pattern formation. Science, 374(6564):eabj0359, 2021

    2021

  45. [53]

    N. J. Middleton and D. S. G. Thomas. World atlas of desertification. Publication Office of the European Union, 1992

    1992

  46. [54]

    D’Odorico and A

    P. D’Odorico and A. Bhattachan. Hydrologic variability in dryland regions: impacts on ecosys- tem dynamics and food security. Phi. Trans. R. Soc. B , 367(1606):3145–3157, 2012

    2012

  47. [55]

    L. A. Gherardi and O. E. Sala. Effect of interannual precipitation variability on dryland productivity: A global synthesis. Glob. Change Biol. , 25(1):269–276, 2019

    2019

  48. [56]

    G. A. Leonov and N. V. Kuznetsov. Time-varying linearization and the perron effects. Int. J. Bifurcat. Chaos, 17(04):1079–1107, 2007

    2007

  49. [57]

    Butler and N

    T. Butler and N. Goldenfeld. Robust ecological pattern formation induced by demographic noise. Phys. Rev. E , 80(3):030902, 2009

    2009

  50. [58]

    Karig, K

    D. Karig, K. M. Martini, T. Lu, N. A. DeLateur, N. Goldenfeld, and R. Weiss. Stochastic turing patterns in a synthetic bacterial population. Proc. Natl. Acad. Sci., 115(26):6572–6577, 2018

    2018

  51. [59]

    Crompton, G

    O. Crompton, G. G. Katul, and S. Thompson. Resistance formulations in shallow overland flow along a hillslope covered with patchy vegetation. Water Resour. Res., 56(5):e2020WR027194, 2020

    2020

  52. [60]

    B. K. Bera, O. Tzuk, J. J. R. Bennett, and E. Meron. Linking spatial self-organization to community assembly and biodiversity. Elife, 10:e73819, 2021

    2021

  53. [61]

    J. J. R. Bennett, B. K. Bera, M. Ferr´ e, H. Yizhaq, S. Getzin, and E. Meron. Phenotypic plasticity: A missing element in the theory of vegetation pattern formation. Proc. Natl. Acad. Sci., 120(50):e2311528120, 2023

    2023

  54. [62]

    Hoyer-Leitzel and S

    A. Hoyer-Leitzel and S. Iams. Impulsive fire disturbance in a savanna model: Tree–grass coexistence states, multiple stable system states, and resilience. Bull. Math. Biol. , 83(11):113, 2021

    2021

  55. [63]

    M. L. Zeeman, K. Meyer, E. Bussmann, A. Hoyer-Leitzel, S. Iams, I. J. Klasky, V. Lee, and S. Ligtenberg. Resilience of socially valued properties of natural systems to repeated disturbance: A framework to support value-laden management decisions. Nat. Resour. Model., 31(3):e12...

    2018

  56. [64]

    Patel and P

    S. Patel and P. De Leenheer. On the spectral radius properties of a key matrix in periodic impulse control. Syst. Control Lett. , 187:105781, 2024

    2024

  57. [65]

    E. H. Colombo, C. L´ opez, and E. Hern´ andez-Garc ´ ıa. Pulsed interaction signals as a route to biological pattern formation. Phys. Rev. Lett. , 130(5):058401, 2023

    2023

  58. [66]

    E. H. Colombo, R. Martinez-Garcia, J. M. Calabrese, C. L´ opez, and E. Hern´ andez-Garc ´ ıa. Pulsed interactions unify reaction–diffusion and spatial nonlocal models for biological pattern formation. J. Stat. Mech.: Theory Exp. , 2024(3):034001, 2024

    2024

  59. [67]

    Hoyer-Leitzel, S

    A. Hoyer-Leitzel, S. M. Iams, A. J. Haslam-Hyde, M. L. Zeeman, and N. H. Fefferman. An immuno-epidemiological model for transient immune protection: A case study for viral respi- ratory infections. Infect. Dis. Model. , 8(3):855–864, 2023

    2023

  60. [68]

    Gasper and M

    G. Gasper and M. Rahman. Basic Hypergeometric Series , volume 96. Cambridge university press, 2011. 41

    2011

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.