The savanna biome is characterised by a continuous vegetation cover, comprised of herbaceous and woody plants. The coexistence of species in arid savannas, where water availability is the main limiting resource for plant growth, provides an apparent contradiction to the classical principle of competitive exclusion. Previous theoretical work using nonspatial models has focussed on the development of an understanding of coexistence mechanisms through the consideration of resource niche separation and ecosystem disturbances. In this paper, we propose that a spatial self-organisation principle, caused by a positive feedback between local vegetation growth and water redistribution, is sufficient for species coexistence in savanna ecosystems. We propose a spatiotemporal ecohydrological model of partial differential equations, based on the Klausmeier reaction-advection-diffusion model for vegetation patterns, to investigate the effects of spatial interactions on species coexistence on sloped terrain. Our results suggest that species coexistence is a possible model outcome, if a balance is kept between the species' average fitness (a measure of a species' competitive abilities in a spatially uniform setting) and their colonisation abilities. Spatial heterogeneities in resource availability are utilised by the superior coloniser (grasses), before it is outcompeted by the species of higher average fitness (trees). A stability analysis of the spatially nonuniform coexistence solutions further suggests that grasses act as ecosystem engineers and facilitate the formation of a continuous tree cover for precipitation levels unable to support a uniform tree density in the absence of a grass species.

Vegetation patterns are a ubiquitous feature of water-deprived ecosystems. Despite the competition for the same limiting resource, coexistence of several plant species is commonly observed. We propose a two-species reaction-diffusion model based on the single-species Klausmeier model, to analytically investigate the existence of states in which both species coexist. Ecologically, the study finds that coexistence is supported if there is a small difference in the plant species' average fitness, measured by the ratio of a species' capabilities to convert water into new biomass to its mortality rate. Mathematically, coexistence is not a stable solution of the system, but both spatially uniform and patterned coexistence states occur as metastable states. In this context, a metastable solution in which both species coexist corresponds to a long transient (exceeding 10 3 years in dimensional parameters) to a stable one-species state. This behaviour is characterised by the small size of a positive eigenvalue which has the same order of magnitude as the average fitness difference between the two species. Two mechanisms causing the occurrence of metastable solutions are established: a spatially uniform unstable equilibrium and a stable one-species pattern which is unstable to the introduction of a competitor. We further discuss effects of asymmetric interspecific competition (e.g. shading) on the metastability property.

Vegetation patterns are a characteristic feature of semi-arid regions. On hillsides these patterns occur as stripes running parallel to the contours. The Klausmeier model, a coupled reaction-advection-diffusion system, is a deliberately simple model describing the phenomenon. In this paper, we replace the diffusion term describing plant dispersal by a more realistic nonlocal convolution integral to account for the possibility of long-range dispersal of seeds. Our analysis focuses on the rainfall level at which there is a transition between uniform vegetation and pattern formation. We obtain results, valid to leading order in the large parameter comparing the rate of water flow downhill to the rate of plant dispersal, for a negative exponential dispersal kernel. Our results indicate that both a wider dispersal of seeds and an increase in dispersal rate inhibit the formation of patterns. Assuming an evolutionary trade-off between these two quantities, mathematically motivated by the limiting behaviour of the convolution term, allows us to make comparisons to existing results for the original reaction-advection-diffusion system. These comparisons show that the nonlocal model always predicts a larger parameter region supporting pattern formation. We then numerically extend the results to other dispersal kernels, showing that the tendency to form patterns depends on the type of decay of the kernel.