pith. sign in

arxiv: 2506.03981 · v1 · submitted 2025-06-04 · 🧮 math.DS · nlin.PS

Beyond water limitation in vegetation-autotoxicity patterning: a cross-diffusion model

Pith reviewed 2026-05-19 10:51 UTC · model grok-4.3

classification 🧮 math.DS nlin.PS MSC 35K5792D40
keywords vegetation patternsautotoxicitycross-diffusionTuring instabilityreaction-diffusion modelsbiomass-toxicity feedbackspatial patternsfast-reaction limit
0
0 comments X

The pith

A biomass-toxicity cross-diffusion model produces stable Turing vegetation patterns without water dynamics.

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

The paper derives a two-equation cross-diffusion model for biomass and toxicity as the fast-reaction limit of a three-species system. This reduced model incorporates autotoxicity effects of growth inhibition, extra mortality, and propagation reduction. Linear stability analysis combined with simulations and continuation shows that the cross-diffusion term produces stable spatial patterns across a wide parameter range. A sympathetic reader would care because the result demonstrates that biomass-toxicity feedback alone can drive vegetation patterning, removing the need for explicit water limitation.

Core claim

The central claim is that a cross-diffusion model derived as the fast-reaction limit of a three-species system with biomass, toxicity, and an auxiliary variable supports the formation of stable Turing patterns in vegetation for a wide range of parameter values, thanks to the cross-diffusion term, even without explicit water dynamics.

What carries the argument

The cross-diffusion term in the reduced biomass-toxicity equations, obtained from the fast-reaction limit of the three-species system, which allows toxicity to influence biomass movement and thereby generate the instability required for Turing patterns.

If this is right

  • Vegetation patterns can form in the absence of water limitation driven only by autotoxicity feedback.
  • The model adds propagation reduction to the previously considered effects of growth inhibition and increased mortality.
  • Pattern formation occurs robustly for broad parameter ranges as verified by linear analysis, simulations, and numerical continuation.
  • The framework applies to any ecological interaction reducible to a fast-reaction limit with cross-diffusion.

Where Pith is reading between the lines

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

  • The same reduction technique could generate cross-diffusion models for autotoxicity in microbial or algal systems.
  • Controlled experiments with uniform moisture but varying toxin levels could test whether patterns appear as predicted.
  • Reducing autotoxin production might offer a practical way to suppress unwanted spatial patterning in crops.
  • The approach connects to other fast-limit derivations in chemical and physical pattern-forming systems.

Load-bearing premise

The three-species system with separated time scales must admit a well-defined fast-reaction limit that produces exactly the cross-diffusion terms used in the two-equation model.

What would settle it

Numerical simulations of the two-equation model with the cross-diffusion coefficient set to zero should produce no spatial patterns while nonzero values should produce them; equivalently, field data showing stable vegetation patterns in autotoxic systems with uniform water supply would support the claim.

Figures

Figures reproduced from arXiv: 2506.03981 by Annalisa Iuorio, Cinzia Soresina, Francesco Giannino.

Figure 1
Figure 1. Figure 1: Schematic representation of the fast-reaction model in Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: shows the steady states of the roots biomass (left) and toxicity (right) obtained in the simulations of Eq. (5)–(6) over the one-dimensional domain Ω = [0, L] for the four scenarios (i)–(iv) outlined above. Here, scenario (i) is represented by a dotted-dashed line, scenario (ii) by a solid line, scenario (iii) by a dashed line, and scenario (iv) by a dotted line. We note that for this specific set of param… view at source ↗
Figure 3
Figure 3. Figure 3: Result of simulation of Eq. (5)–(6) on a two-dimensional domain squared domain Ω (L = 8 metres) at t = tfin for the roots biomass R (upper panel) and toxicity concentration T (lower panel) for the four different scenar￾ios (i), (ii), (iii), (iv) corresponding to the first, second, third, and fourth column, respectively. Other parameter values are fixed as in Eq. (21). 4.2. Bifurcations and numerical contin… view at source ↗
Figure 4
Figure 4. Figure 4: Region of admissible σ-values as a function of γ and s (describing the growth-inhibition and extra-mortality effect, respectively) in (γ, s, σ)-space (all other parameter values are fixed as in Eq. (21)). This grey region is delimited by the two-dimensional surfaces σ = σL (as defined in the Turing condition (18)) and σ = dR, in blue and cyan respectively. software pde2path to compute the bifurcation diagr… view at source ↗
Figure 5
Figure 5. Figure 5: Bifurcation diagram with respect to (a) σ in scenario (ii) and (b) parameter s in scenario (iii). The black line corresponds to the homogenous branch, while branches of nonhomogenous solutions are shown in grey and orange. In (a), the orange branch is the second branch bifurcating subcritically from the homogeneous branch. In (b), the orange branch is the first branch bifurcating subcritically from the hom… view at source ↗
read the original abstract

Many mathematical models describing vegetation patterns are based on biomass--water interactions, due to the impact of this limited resource in arid and semi-arid environments. However, in recent years, a novel biological factor called autotoxicity has proved to play a key role in vegetation spatiotemporal dynamics, particularly by inhibiting biomass growth and increasing its natural mortality rate. In a standard reaction-diffusion framework, biomass-toxicity dynamics alone are unable to support the emergence of stable spatial patterns. In this paper, we derive a cross-diffusion model for biomass and toxicity dynamics as the fast-reaction limit of a three-species system involving dichotomy and different time scales. Within this general framework, in addition to growth inhibition and extra-mortality already considered in previous studies, the additional effect of ''propagation reduction'' induced by autotoxicity on vegetation dynamics is obtained. By combining linearised analysis, simulations, and continuation, we investigate the formation of spatial patterns. Thanks to the cross-diffusion term, for the first time, a spatial model based solely on biomass-toxicity feedback without explicit water dynamics supports the formation of stable (Turing) vegetation patterns for a wide range of parameter values.

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

Summary. The paper derives a two-equation cross-diffusion system for biomass and toxicity as the fast-reaction limit of a three-species model with dichotomy and separated timescales. It performs linear stability analysis, numerical simulations, and continuation to show that the cross-diffusion terms enable stable Turing patterns over a wide parameter range, providing the first such spatial model based solely on biomass-toxicity feedback without explicit water dynamics.

Significance. If the reduced model is a faithful limit, the work supplies a new mechanism for vegetation patterning driven by autotoxicity alone. The combination of linearised analysis, simulations, and continuation provides concrete evidence for pattern formation across parameter values, which is a methodological strength. The explicit fast-reaction derivation from an independent three-species system (rather than ad-hoc fitting) adds transparency to the modelling assumptions.

major comments (2)
  1. [Model derivation] Model derivation section: the reduction to the specific cross-diffusion terms in the two-equation system is presented as the fast-reaction limit of the three-species model, but no convergence theorem, uniform estimates, or error bounds are supplied to guarantee that the limit exists and reproduces exactly those terms (including the propagation-reduction effect). This is load-bearing for the central claim, because the Turing analysis and all subsequent simulations are performed on the reduced system; if additional reaction terms survive or the cross terms differ, the reported patterns do not correspond to the stated biological mechanism.
  2. [Abstract and §4] Abstract and §4 (numerical results): the claim of patterns 'for a wide range of parameter values' is not accompanied by explicit intervals, robustness metrics, or error measures on the continuation or simulation results. Without these, the quantitative support for the 'wide range' assertion remains only moderate.
minor comments (2)
  1. [Model derivation] Notation for the cross-diffusion coefficients should be introduced with a clear table or explicit definitions immediately after the reduced equations to improve readability.
  2. [Numerical results] Figure captions for the simulation panels should state the precise parameter values and initial conditions used, rather than referring only to 'typical' values.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation of the significance of our work and for the detailed, constructive comments. We respond to each major comment below.

read point-by-point responses
  1. Referee: [Model derivation] Model derivation section: the reduction to the specific cross-diffusion terms in the two-equation system is presented as the fast-reaction limit of the three-species model, but no convergence theorem, uniform estimates, or error bounds are supplied to guarantee that the limit exists and reproduces exactly those terms (including the propagation-reduction effect). This is load-bearing for the central claim, because the Turing analysis and all subsequent simulations are performed on the reduced system; if additional reaction terms survive or the cross terms differ, the reported patterns do not correspond to the stated biological mechanism.

    Authors: We acknowledge that the derivation in the manuscript is formal and does not include a rigorous convergence theorem, uniform estimates, or error bounds. The reduction is obtained by the standard singular-perturbation procedure for fast-reaction limits with separated timescales, which produces the stated cross-diffusion terms and the propagation-reduction effect. In the revision we will expand the model-derivation section with an explicit step-by-step asymptotic calculation and add a short discussion of the modelling assumptions under which the formal limit is expected to hold, together with references to analogous reductions in the literature. A complete rigorous convergence analysis lies outside the scope of the present biologically focused study. revision: partial

  2. Referee: [Abstract and §4] Abstract and §4 (numerical results): the claim of patterns 'for a wide range of parameter values' is not accompanied by explicit intervals, robustness metrics, or error measures on the continuation or simulation results. Without these, the quantitative support for the 'wide range' assertion remains only moderate.

    Authors: We agree that the quantitative support for the 'wide range' claim can be strengthened. In the revised manuscript we will augment §4 with explicit intervals for the key parameters (cross-diffusion coefficients, reaction rates, and diffusion ratios) over which stable Turing patterns are found by numerical continuation. We will also report robustness metrics (e.g., the fraction of sampled parameter space yielding patterns) and include numerical error indicators from the simulations to make the evidence more precise. revision: yes

Circularity Check

0 steps flagged

Derivation via fast-reaction limit of independent three-species system is self-contained and non-circular

full rationale

The paper derives the two-equation cross-diffusion model explicitly as the fast-reaction limit of a separate three-species system with dichotomy and separated timescales. This limit process supplies the cross terms as an output rather than defining them in terms of the target Turing patterns or fitting them to pattern data. The subsequent linear stability analysis, simulations, and continuation are performed on the derived reduced system; they do not feed back into the derivation step. No self-definitional loop, fitted-input-as-prediction, or load-bearing self-citation chain appears in the model-construction chain. The cited prior work on autotoxicity effects is used only for the reaction terms already studied elsewhere and is not required to justify the cross-diffusion structure itself. The derivation therefore remains independent of the final pattern-formation claims.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the fast-reaction limit and the specific functional forms chosen for autotoxicity effects.

axioms (1)
  • domain assumption The three-species system with dichotomy and different time scales possesses a well-defined fast-reaction limit yielding the cross-diffusion model.
    Invoked to obtain the two-equation system from the underlying three-species dynamics.

pith-pipeline@v0.9.0 · 5738 in / 1212 out tokens · 28239 ms · 2026-05-19T10:51:26.226858+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

52 extracted references · 52 canonical work pages

  1. [1]

    Abbas, F

    M. Abbas, F. Giannino, A. Iuorio, Z. Ahmad, and F. Calabrò. PDE models for vegetation biomass and autotoxicity.Mathematics and Computers in Simulation, 228:386–401, 2025

  2. [2]

    Abbas, F

    M. Abbas, F. Giannino, A. Iuorio, and F. Calabrò. Enhanced forecasting of biomass-toxicity- water models using numerical simulations.Advances in Computational Science and Engineering, 2(4):380–404, 2024

  3. [3]

    Allegrezza, G

    M. Allegrezza, G. Bonanomi, M. Zotti, F. Cartenì, M. Moreno, L. Olivieri, M. Garbarino, G. Tesei, F. Giannino, and S. Mazzoleni. Biogeography and shape of fungal fairy rings in the Apennine mountains, Italy.Journal of Biogeography, 49(2):353–363, 2022

  4. [4]

    Bastiaansen, M

    R. Bastiaansen, M. Chirilus-Bruckner, and A. Doelman. Pulse solutions for an extended Klaus- meier model with spatially varying coefficients.SIAM Journal on Applied Dynamical Systems, 19(1):1–57, 2020

  5. [5]

    Bastiaansen and A

    R. Bastiaansen and A. Doelman. The dynamics of disappearing pulses in a singularly perturbed reaction–diffusion system with parameters that vary in time and space.Physica D: Nonlinear Phenomena, 388:45–72, 2019

  6. [6]

    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.Ecology Letters, 23(3):414–429, 2020

  7. [7]

    B.K. Bera, O. Tzuk, J.J.R. Bennett, U. Dieckmann, and E. Meron. Can spatial self-organization inhibit evolutionary adaptation?Journal of The Royal Society Interface, 22(222):20240454, 2025

  8. [8]

    Bonanomi, G

    G. Bonanomi, G. Incerti, A. Stinca, F. Cartenì, F. Giannino, and S. Mazzoleni. Ring formation in clonal plants.Community Ecology, 15(1):77–86, 2014

  9. [9]

    Bonanomi, M

    G. Bonanomi, M. Zotti, M. Idbella, P. Termolino, V. De Micco, and S. Mazzoleni. Field ev- idence for litter and self-DNA inhibitory effects on Alnus glutinosa roots. New Phytologist, 236(2):399–412, 2022. 14

  10. [10]

    Cartenì, A

    F. Cartenì, A. Marasco, G. Bonanomi, S. Mazzoleni, M. Rietkerk, and F. Giannino. Negative plant soil feedback explaining ring formation in clonal plants. Journal of Theoretical Biology, 313:153–161, 2012

  11. [11]

    Carter and A

    P. Carter and A. Doelman. Traveling stripes in the Klausmeier model of vegetation pattern formation. SIAM Journal on Applied Mathematics, 78(6):3213–3237, 2018

  12. [12]

    Carter, A

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

  13. [13]

    Consolo, G

    G. Consolo, G. Grifò, and G. Valenti. Modelling vegetation patterning on sloped terrains: The role of toxic compounds.Physica D: Nonlinear Phenomena, 459:134020, 2023

  14. [14]

    De Baets, J

    S. De Baets, J. Poesen, A. Knapen, G.G. Barberá, and J.A. Navarro. Root characteristics of rep- resentative Mediterranean plant species and their erosion-reducing potential during concentrated runoff. Plant and Soil, 294(1–2):169–183, 2007

  15. [15]

    Non-triangularcross-diffusionsystemswithpredator–preyreaction terms

    L.DesvillettesandC.Soresina. Non-triangularcross-diffusionsystemswithpredator–preyreaction terms. Ricerche di Matematica, 68:295–314, 2019

  16. [16]

    Eigentler

    L. Eigentler. Intraspecific competition in models for vegetation patterns: Decrease in resilience to aridity and facilitation of species coexistence.Ecological Complexity, 42:100835, 2020

  17. [17]

    Eliaš, D

    J. Eliaš, D. Hilhorst, M. Mimura, and Y. Morita. Singular limit for a reaction-diffusion-ODE system in a neolithic transition model.Journal of Differential Equations, 295:39–69, 2021

  18. [18]

    Eliaš, H

    J. Eliaš, H. Izuhara, M. Mimura, and B.Q. Tang. An aggregation model of cockroaches with fast-or-slow motion dichotomy.Journal of Mathematical Biology, 85(3):28, 2022

  19. [19]

    Fanelli, C

    D. Fanelli, C. Cianci, and F. Di Patti. Turing instabilities in reaction-diffusion systems with cross diffusion. The European Physical Journal B, 86(4):1–8, 2013

  20. [20]

    Gandhi, L

    P. Gandhi, L. Liu, and M. Silber. Stochastic pulsed precipitation model of dryland vegetation pattern formation. arXiv preprint arXiv:2112.14851, 2021

  21. [21]

    Gandhi, L

    P. Gandhi, L. Werner, S. Iams, K. Gowda, and M. Silber. A topographic mechanism for arcing of dryland vegetation bands.Journal of The Royal Society Interface, 15(147):20180508, 2018

  22. [22]

    Gilad, J

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

  23. [23]

    Grifò, A

    G. Grifò, A. Iuorio, and F. Veerman. Far-from-equilibrium traveling pulses in sloped semiarid environments driven by autotoxicity effects.SIAM Journal on Applied Mathematics, 85(1):188– 209, 2025

  24. [24]

    M. Iida, M. Mimura, and H. Ninomiya. Diffusion, cross-diffusion and competitive interaction. Journal of Mathematical Biology, 53(4):617–641, 2006

  25. [25]

    Iuorio and F

    A. Iuorio and F. Veerman. The influence of autotoxicity on the dynamics of vegetation spots. Physica D: Nonlinear Phenomena, 427:133015, 2021

  26. [26]

    Jaïbi, A

    O. Jaïbi, A. Doelman, M. Chirilus-Bruckner, and E. Meron. The existence of localized vegeta- tion patterns in a systematically reduced model for dryland vegetation.Physica D: Nonlinear Phenomena, 412:132637, 2020. 15

  27. [27]

    S. Kéfi, M. Rietkerk, M. van Baalen, and M. Loreau. Local facilitation, bistability and transitions in arid ecosystems.Theoretical Population Biology, 71(3):367–379, 2007

  28. [28]

    S. Kéfi, V. Guttal, W.A. Brock, S.R. Carpenter, V. Dakos, et al. Early warning signals of ecological transitions: Methods for spatial patterns.PLoS ONE, 9(3):e92097, 2014

  29. [29]

    Klausmeier

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

  30. [30]

    Kuehn and C

    C. Kuehn and C. Soresina. Numerical continuation for a fast-reaction system and its cross-diffusion limit. SN Partial Differential Equations and Applications, 1(2):1–26, 2020

  31. [31]

    Marasco, F

    A. Marasco, F. Giannino, and A. Iuorio. Modelling competitive interactions and plant–soil feed- back in vegetation dynamics.Ricerche di Matematica, 69(2):553–577, 2020

  32. [32]

    Marasco, A

    A. Marasco, A. Iuorio, F. Cartenì, G. Bonanomi, D.M. Tartakovsky, S. Mazzoleni, and F. Gian- nino. Vegetation pattern formation due to interactions between water availability and toxicity in plant-soil feedback.Bulletin of Mathematical Biology, 76(11):2866–2883, 2014

  33. [33]

    Martínez, O

    F. Martínez, O. Merino, A. Martín, D. García Martín, and J. Merino. Belowground structure and production in a mediterranean sand dune shrub community.Plant and Soil, 201:209–216, 1998

  34. [34]

    Martínez-García, C

    R. Martínez-García, C. Cabal, J.M. Calabrese, E. Hernández-García, C.E. Tarnita, C. López, and J.A. Bonachela. Integrating theory and experiments to link local mechanisms and ecosystem-level consequences of vegetation patterns in drylands.Chaos, Solitons & Fractals, 166:112881, 2023

  35. [35]

    Mazzoleni

    S. Mazzoleni. Composition comprising nucleic acids of parasitic, pathogenic or weed biological systems for inhibiting and/or controlling the growth of said systems, 2015. Patent App. WO 02/0624 A2

  36. [36]

    Mazzoleni, G

    S. Mazzoleni, G. Bonanomi, F. Giannino, G. Incerti, S.C. Dekker, and M. Rietkerk. Modelling the effects of litter decomposition on tree diversity patterns.Ecological Modelling, 221(23):2784–2792, 2010

  37. [37]

    Mazzoleni, G

    S. Mazzoleni, G. Bonanomi, G. Incerti, M.L. Chiusano, V. Lanzotti, et al. Inhibitory and toxic effects of extracellular self-DNA in litter: a mechanism for negative plant-soil feedbacks?New Phytologist, 205(3):1195–1210, 2014

  38. [38]

    E. Meron. Pattern-formation approach to modelling spatially extended ecosystems. Ecological Modelling, 234:70–82, 2012

  39. [39]

    Morgan, C

    J. Morgan, C. Soresina, B.Q. Tang, and B.-N. Tran. Singular limit and convergence rate via projection method in a model for plant-growth dynamics with autotoxicity. arXiv preprint arXiv:2408.06177, 2024

  40. [40]

    Evasion of tipping in complex systems through spatial pattern formation.Science, 374(6564):eabj0359, 2021

    M.Rietkerk, R.Bastiaansen, S.Banerjee, J.vandeKoppel, M.Baudena, andA.Doelman. Evasion of tipping in complex systems through spatial pattern formation.Science, 374(6564):eabj0359, 2021

  41. [41]

    Ritchie, A.L

    J.S. Ritchie, A.L. Krause, and R.A. Van Gorder. Turing and wave instabilities in hyperbolic reaction–diffusion systems: The role of second-order time derivatives and cross-diffusion terms on pattern formation. Annals of Physics, 444:169033, 2022

  42. [42]

    Rottschäfer and A

    V. Rottschäfer and A. Doelman. On the transition from the Ginzburg-Landau equation to the extended Fisher-Kolmogorov equation.Physica D: Nonlinear Phenomena, 118(3-4):261–292, 1998. 16

  43. [43]

    Salvatori, M

    N. Salvatori, M. Moreno, M. Zotti, A. Iuorio, F. Cartenì, G. Bonanomi, S. Mazzoleni, and F. Gian- nino. Process based modelling of plants-fungus interactions explains fairy ring types and dynamics. Scientific Reports, 13(1):19918, 2023

  44. [44]

    Sherratt

    J.A. Sherratt. Pattern solutions of the Klausmeier model for banded vegetation in semi-arid environments I.Nonlinearity, 23(10):2657–2675, 2010

  45. [45]

    Shigesada, K

    N. Shigesada, K. Kawasaki, and E. Teramoto. Spatial segregation of interacting species.Journal of Theoretical Biology, 79(1):83–99, 1979

  46. [46]

    Soresina

    C. Soresina. Hopf bifurations in the full SKT model and where to find them. Discrete and Continuous Dynamical Systems Series S, 15(9):2673–2693, 2022

  47. [47]

    H. Uecker. Hopf bifurcation and time periodic orbits withpde2path– Algorithms and applications. Communications in Computational Physics, 25(3):812–852, 2019

  48. [48]

    H. Uecker. Numerical Continuation and Bifurcation in Nonlinear PDEs. Society for Industrial and Applied Mathematics, 2021

  49. [49]

    Villar-Sepúlveda, A.R

    E. Villar-Sepúlveda, A.R. Champneys, and A.L. Krause. Designing reaction-cross-diffusion sys- tems with Turing and wave instabilities.arXiv preprint arXiv:2409.06860, 2024

  50. [50]

    von Hardenberg, E

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

  51. [51]

    C. Yang, J. Li, F. Zhang, N. Liu, and Y. Zhang. The optimal Redfield N: P ratio caused by fairy ring fungi stimulates plant productivity in the temperate steppe of China.Fungal Ecology, 34:91–98, 2018

  52. [52]

    Zotti, G

    M. Zotti, G. Bonanomi, and S. Mazzoleni. Fungal fairy rings: history, ecology, dynamics and engineering functions. IMA Fungus, 16:e138320, 2025. 17