Intermittent precipitation and spatial Allee effects drive irregular vegetation patterns in semiarid ecosystems

Vegetation in semi-arid ecosystems frequently organizes into spatially heterogeneous mosaics that regulate ecosystem functioning, productivity, and resilience. These patterns arise from local biological interactions, including facilitation among neighboring plants and competition for limiting resources. Classical theoretical approaches have attributed such organization to scale-dependent feedbacks, predicting regular spatial patterns and abrupt transitions to collapse. However, growing empirical and theoretical evidence reveal that environmental variability and demographic stochasticity can fundamentally reshape spatial organization, driving irregular clusters, dynamic mosaics, and gradual rather than catastrophic vegetation declines. In drylands, rainfall variability is a dominant source of environmental forcing: precipitation typically occurs in short, irregular pulses that transiently enhance survival and recruitment before competitive interactions again dominate. Near persistence thresholds, ecosystem dynamics are therefore governed not only by average climatic conditions but also by the timing and spatial coincidence of favorable events. Under these conditions, positive density dependence and local facilitation can critically determine whether vegetation patches persist, expand, or collapse. Here, we develop an individual-based model that integrates intermittent precipitation with local Allee effects to examine how stochastic rainfall shapes spatial organization and persistence. We show that the interaction between pulsed resource availability and density-dependent survival generates irregular cluster structures and strongly modulates extinction risk, with resilience emerging from local spatial covariance and neighborhood density rather than from total biomass alone. These results highlight the importance of individual-level, stochastic processes in determining ecosystem resilience.

💡 Research Summary

The paper addresses a long‑standing discrepancy between classic deterministic models of dry‑land vegetation, which predict regular, crystalline patterns and abrupt catastrophic shifts, and the irregular, dynamic mosaics observed in nature. The authors develop an individual‑based model (IBM) that explicitly incorporates two key stochastic processes: intermittent precipitation events and a local Allee effect. Plants are represented as discrete agents occupying a two‑dimensional periodic domain. Competition and facilitation are modeled with Gaussian interaction kernels (ψ_c and ψ_f) whose amplitudes (γ_c, γ_f) and spatial scales (σ_c, σ_f) determine how the presence of neighbours influences an individual’s death and birth rates. Death follows a quadratic dependence on cumulative competition (D_n = X_n² + c), while birth is linearly enhanced by cumulative facilitation (B_n = ρ + Y_n). By setting ρ = 0 the model forces population growth to rely entirely on positive density dependence, i.e., the Allee effect.

Precipitation is introduced as a Poisson process with rate λ; each rain pulse lasts for a fixed duration d. During a pulse the interaction kernels are temporarily altered: either facilitation strength γ_f is increased to γ_f⁺ or competition strength γ_c is reduced to γ_c⁻. Thus, a rain event temporarily tilts the balance toward positive feedback, allowing small clusters to survive and expand. Between pulses the system reverts to its baseline kernels, re‑establishing competition‑dominated dynamics.

The authors quantify clustering using a pair‑correlation function C(r, t). Values above one indicate aggregation at distance r, while values below one signal regular spacing. Simulations explore a wide range of λ and d, revealing several robust patterns:

1. Pulse frequency and duration control cluster size and persistence. Low λ and short d keep most individuals below the Allee threshold, leading to rapid extinction of small patches and low overall biomass. High λ or long d sustain elevated facilitation, allowing clusters to grow, merge, and form large, irregular patches.

2. Local density, not total biomass, predicts resilience. Even when total plant cover is identical across scenarios, those with higher local neighbourhood density (i.e., where the Allee threshold is exceeded) exhibit dramatically lower extinction risk. This highlights spatial covariance as the true driver of ecosystem stability.

3. Gradual rather than abrupt transitions. The IBM replaces the saddle‑node bifurcations of PDE models with smooth declines in mean biomass as precipitation becomes scarcer. Cluster‑size distributions shift from scale‑free power laws under frequent rain to truncated distributions under drought, matching field observations of irregular patterns.

4. Early‑warning signals are obscured. Because stochasticity and local Allee dynamics dominate, traditional indicators such as increasing spatial variance or autocorrelation become unreliable predictors of collapse.

The discussion situates these findings within the broader literature. It emphasizes that deterministic reaction‑diffusion or advection‑diffusion frameworks, while analytically tractable, miss essential demographic noise and the temporal structure of water input. By contrast, the IBM captures finite‑size effects, stochastic recruitment, and the crucial role of intermittent resource pulses. The authors argue that management actions aiming to enhance local facilitation—through artificial watering, soil amendments, or planting of nurse species—could effectively raise the local Allee threshold and improve resilience, even if mean annual precipitation does not change.

In conclusion, the study demonstrates that the interplay of intermittent precipitation and spatial Allee effects can generate the irregular vegetation mosaics characteristic of semi‑arid ecosystems and that resilience emerges from local spatial structure rather than from bulk productivity. The IBM provides a mechanistic, bottom‑up platform for exploring climate‑change scenarios, informing conservation strategies, and refining theoretical predictions about pattern formation and regime shifts in drylands.