Striped pattern selection by advective reaction-diffusion systems: resilience of banded vegetation on slopes

Phenotypic plasticity: A missing element in the theory of vegetation pattern formation

Regular spatial patterns of vegetation are a common sight in drylands. Their formation is a population-level response to water stress that increases water availability for the few via partial plant mortality. At the individual level, plants can also adapt to water stress by changing their phenotype. Phenotypic plasticity of individual plants and spatial patterning of plant populations have extensively been studied independently, but the likely interplay between the two robust mechanisms has remained unexplored. In this paper, we incorporate phenotypic plasticity into a multi-level theory of vegetation pattern formation and use a fascinating ecological phenomenon, the Namibian "fairy circles," to demonstrate the need for such a theory. We show that phenotypic changes in the root structure of plants, coupled with pattern-forming feedback within soil layers, can resolve two puzzles that the current theory fails to explain: observations of multi-scale patterns and the absence of theoretically predicted large-scale stripe and spot patterns along the rainfall gradient. Importantly, we find that multi-level responses to stress unveil a wide variety of more effective stress-relaxation pathways, compared to single-level responses, implying a previously underestimated resilience of dryland ecosystems.

Pattern Formation in Mesic Savannas

We analyze a spatially extended version of a well-known model of forest-savanna dynamics, which presents as a system of nonlinear partial integro-differential equations, and study necessary conditions for pattern-forming bifurcations. Homogeneous solutions dominate the dynamics of the standard forest-savanna model, regardless of the length scales of the various spatial processes considered. However, several different pattern-forming scenarios are possible upon including spatial resource limitation, such as competition for water, soil nutrients, or herbivory effects. Using numerical simulations and continuation, we study the nature of the resulting patterns as a function of system parameters and length scales, uncovering subcritical pattern-forming bifurcations and observing significant regions of multistability for realistic parameter regimes. Finally, we discuss our results in the context of extant savanna-forest modeling efforts and highlight ongoing challenges in building a unifying mathematical model for savannas across different rainfall levels.

Effect of nonlocal grazing on dry-land vegetation dynamics

Dry-land ecosystems have become a matter of grave concern, due to the growing threat of land degradation and bioproductivity loss. Self-organized vegetation patterns are a remarkable characteristic of these ecosystems; apart from being visually captivating, patterns modulate the system response to increasing environmental stress. Empirical studies hinted that herbivory is one the key regulatory mechanisms behind pattern formation and overall ecosystem functioning. However, most of the mathematical models have taken a mean-field strategy to grazing; foraging has been considered to be independent of spatial distribution of vegetation. To this end, an extended version of the celebrated plant-water model due to Klausmeier has been taken as the base here. To encompass the effect of heterogeneous vegetation distribution on foraging intensity and subsequent impact on entire ecosystem, grazing is considered here to depend on spatially weighted average vegetation density instead of density at a particular point. Moreover, varying influence of vegetation at any location over gazing elsewhere is incorporated by choosing a suitable averaging function. A comprehensive analysis demonstrates that inclusion of spatial nonlocality alters the understanding of system dynamics significantly. The grazing ecosystem is found to be more resilient to increasing aridity than it was anticipated to be in earlier studies on nonlocal grazing. The system response to rising environmental pressure is also observed to vary depending on the grazer. Obtained results also suggest the possibility of multistability due to the history dependence of the system response. Overall, this work indicates that the spatial heterogeneity in grazing intensity has a decisive role to play in the functioning of water-limited ecosystems.

Beyond Turing: far-from-equilibrium patterns and mechano-chemical feedback

Turing patterns are commonly understood as specific instabilities of a spatially homogeneous steady state, resulting from activator-inhibitor interaction destabilized by diffusion. We argue that this view is restrictive and its agreement with biological observations is problematic. We present two alternatives to the classical Turing analysis of patterns. First, we employ the abstract framework of evolution equations to enable the study of far-from-equilibrium patterns. Second, we introduce a mechano-chemical model, with the surface on which the pattern forms being dynamic and playing an active role in the pattern formation, effectively replacing the inhibitor. We highlight the advantages of these two alternatives vis-à-vis the classical Turing analysis, and give an overview of recent results and future challenges for both approaches. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Coupling Turing stripes to active flows

We numerically solve the active nematohydrodynamic equations of motion, coupled to a Turing reaction-diffusion model, to study the effect of active nematic flow on the stripe patterns resulting from a Turing instability. If the activity is uniform across the system, the Turing patterns dissociate when the flux from active advection balances that from the reaction-diffusion process. If the activity is coupled to the concentration of Turing morphogens, and neighbouring stripes have equal and opposite activity, the system self organises into a pattern of shearing flows, with stripes tending to fracture and slip sideways to join their neighbours. We discuss the role of active instabilities in controlling the crossover between these limits. Our results are of relevance to mechanochemical coupling in biological systems.

A Comparison of the "Reduced Losses" and "Increased Production" Models for Mussel Bed Dynamics

Self-organised regular pattern formation is one of the foremost examples of the development of complexity in ecosystems. Despite the wide array of mechanistic models that have been proposed to understand pattern formation, there is limited general understanding of the feedback processes causing pattern formation in ecosystems, and how these affect ecosystem patterning and functioning. Here we propose a generalised model for pattern formation that integrates two types of within-patch feedback: amplification of growth and reduction of losses. Both of these mechanisms have been proposed as causing pattern formation in mussel beds in intertidal regions, where dense clusters of mussels form, separated by regions of bare sediment. We investigate how a relative change from one feedback to the other affects the stability of uniform steady states and the existence of spatial patterns. We conclude that there are important differences between the patterns generated by the two mechanisms, concerning both biomass distribution in the patterns and the resilience of the ecosystems to disturbances.

An integrodifference model for vegetation patterns in semi-arid environments with seasonality

Vegetation patterns are a characteristic feature of semi-deserts occurring on all continents except Antarctica. In some semi-arid regions, the climate is characterised by seasonality, which yields a synchronisation of seed dispersal with the dry season or the beginning of the wet season. We reformulate the Klausmeier model, a reaction–advection–diffusion system that describes the plant–water dynamics in semi-arid environments, as an integrodifference model to account for the temporal separation of plant growth processes during the wet season and seed dispersal processes during the dry season. The model further accounts for nonlocal processes involved in the dispersal of seeds. Our analysis focusses on the onset of spatial patterns. The Klausmeier partial differential equations (PDE) model is linked to the integrodifference model in an appropriate limit, which yields a control parameter for the temporal separation of seed dispersal events. We find that the conditions for pattern onset in the integrodifference model are equivalent to those for the continuous PDE model and hence independent of the time between seed dispersal events. We thus conclude that in the context of seed dispersal, a PDE model provides a sufficiently accurate description, even if the environment is seasonal. This emphasises the validity of results that have previously been obtained for the PDE model. Further, we numerically investigate the effects of changes to seed dispersal behaviour on the onset of patterns. We find that long-range seed dispersal inhibits the formation of spatial patterns and that the seed dispersal kernel’s decay at infinity is a significant regulator of patterning.

Flow Induced Symmetry Breaking in a Conceptual Polarity Model

Important cellular processes, such as cell motility and cell division, are coordinated by cell polarity, which is determined by the non-uniform distribution of certain proteins. Such protein patterns form via an interplay of protein reactions and protein transport. Since Turing’s seminal work, the formation of protein patterns resulting from the interplay between reactions and diffusive transport has been widely studied. Over the last few years, increasing evidence shows that also advective transport, resulting from cytosolic and cortical flows, is present in many cells. However, it remains unclear how and whether these flows contribute to protein-pattern formation. To address this question, we use a minimal model that conserves the total protein mass to characterize the effects of cytosolic flow on pattern formation. Combining a linear stability analysis with numerical simulations, we find that membrane-bound protein patterns propagate against the direction of cytoplasmic flow with a speed that is maximal for intermediate flow speed. We show that the mechanism underlying this pattern propagation relies on a higher protein influx on the upstream side of the pattern compared to the downstream side. Furthermore, we find that cytosolic flow can change the membrane pattern qualitatively from a peak pattern to a mesa pattern. Finally, our study shows that a non-uniform flow profile can induce pattern formation by triggering a regional lateral instability.

The effect of climate change on the resilience of ecosystems with adaptive spatial pattern formation

In a rapidly changing world, quantifying ecosystem resilience is an important challenge. Historically, resilience has been defined via models that do not take spatial effects into account. These systems can only adapt via uniform adjustments. In reality, however, the response is not necessarily uniform, and can lead to the formation of (self-organised) spatial patterns - typically localised vegetation patches. Classical measures of resilience cannot capture the emerging dynamics in spatially self-organised systems, including transitions between patterned states that have limited impact on ecosystem structure and productivity. We present a framework of interlinked phase portraits that appropriately quantifies the resilience of patterned states, which depends on the number of patches, the distances between them and environmental conditions. We show how classical resilience concepts fail to distinguish between small and large pattern transitions, and find that the variance in interpatch distances provides a suitable indicator for the type of imminent transition. Subsequently, we describe the dependency of ecosystem degradation based on the rate of climatic change: slow change leads to sporadic, large transitions, whereas fast change causes a rapid sequence of smaller transitions. Finally, we discuss how pre-emptive removal of patches can minimise productivity losses during pattern transitions, constituting a viable conservation strategy.

Model of pattern formation in marsh ecosystems with nonlocal interactions

Smooth cordgrass Spartina alterniflora is a grass species commonly found in tidal marshes. It is an ecosystem engineer, capable of modifying the structure of its surrounding environment through various feedbacks. The scale-dependent feedback between marsh grass and sediment volume is particularly of interest. Locally, the marsh vegetation attenuates hydrodynamic energy, enhancing sediment accretion and promoting further vegetation growth. In turn, the diverted water flow promotes the formation of erosion troughs over longer distances. This scale-dependent feedback may explain the characteristic spatially varying marsh shoreline, commonly observed in nature. We propose a mathematical framework to model grass-sediment dynamics as a system of reaction-diffusion equations with an additional nonlocal term quantifying the short-range positive and long-range negative grass-sediment interactions. We use a Mexican-hat kernel function to model this scale-dependent feedback. We perform a steady state biharmonic approximation of our system and derive conditions for the emergence of spatial patterns, corresponding to a spatially varying marsh shoreline. We find that the emergence of such patterns depends on the spatial scale and strength of the scale-dependent feedback, specified by the width and amplitude of the Mexican-hat kernel function.

Labyrinth ice pattern formation induced by near-infrared irradiation

Patterns are broad phenomena that relate to biology, chemistry, and physics. The dendritic growth of crystals is the most well-known ice pattern formation process. Tyndall figures are water-melting patterns that occur when ice absorbs light and becomes superheated. Here, we report a previously undescribed ice and water pattern formation process induced by near-infrared irradiation that heats one phase more than the other in a two-phase system. The pattern formed during the irradiation of ice crystals tens of micrometers thick in solution near equilibrium. Dynamic holes and a microchannel labyrinth then formed in specific regions and were characterized by a typical distance between melted points. We concluded that the differential absorption of water and ice was the driving force for the pattern formation. Heating ice by laser absorption might be useful in applications such as the cryopreservation of biological samples.

Large scale patterns in mussel beds: stripes or spots?

An aerial view of an intertidal mussel bed often reveals large scale striped patterns aligned perpendicular to the direction of the tide; dense bands of mussels alternate periodically with near bare sediment. Experimental work led to the formulation of a set of coupled partial differential equations modelling a mussel-algae interaction, which proved pivotal in explaining the phenomenon. The key class of model solutions to consider are one-dimensional periodic travelling waves (wavetrains) that encapsulate the abundance of peak and trough mussel densities observed in practice. These solutions may, or may not, be stable to small perturbations, and previous work has focused on determining the ecologically relevant (stable) wavetrain solutions in terms of model parameters. The aim of this paper is to extend this analysis to two space dimensions by considering the full stripe pattern solution in order to study the effect of transverse two-dimensional perturbations-a more true to life problem. Using numerical continuation techniques, we find that some striped patterns that were previously deemed stable via the consideration of the associated wavetrain solution, are in fact unstable to transverse two-dimensional perturbations; and numerical simulation of the model shows that they break up to form regular spotted patterns. In particular, we show that break up of stripes into spots is a consequence of low tidal flow rates. Our consideration of random algal movement via a dispersal term allows us to show that a higher algal dispersal rate facilitates the formation of stripes at lower flow rates, but also encourages their break up into spots. We identify a novel hysteresis effect in mussel beds that is a consequence of transverse perturbations.

Grazing Away the Resilience of Patterned Ecosystems

Ecosystems' responses to changing environmental conditions can be modulated by spatial self-organization. A prominent example of this can be found in drylands, where formation of vegetation patterns attenuates the magnitude of degradation events in response to decreasing rainfall. In model studies, the pattern wavelength responds to changing conditions, which is reflected by a rather gradual decline in biomass in response to decreasing rainfall. Although these models are spatially explicit, they have adopted a mean-field approach to grazing. By taking into account spatial variability when modeling grazing, we find that (over)grazing can lead to a dramatic shift in biomass, so that degradation occurs at rainfall rates that would otherwise still maintain a relatively productive ecosystem. Moreover, grazing increases the resilience of degraded ecosystem states. Consequently, restoration of degraded ecosystems could benefit from the introduction of temporary small-scale exclosures to escape from the basin of attraction of degraded states.

Stabilizing a homoclinic stripe

For a large class of reaction-diffusion systems with large diffusivity ratio, it is well known that a two-dimensional stripe (whose cross-section is a one-dimensional homoclinic spike) is unstable and breaks up into spots. Here, we study two effects that can stabilize such a homoclinic stripe. First, we consider the addition of anisotropy to the model. For the Schnakenberg model, we show that (an infinite) stripe can be stabilized if the fast-diffusing variable (substrate) is sufficiently anisotropic. Two types of instability thresholds are derived: zigzag (or bending) and break-up instabilities. The instability boundaries subdivide parameter space into three distinct zones: stable stripe, unstable stripe due to bending and unstable due to break-up instability. Numerical experiments indicate that the break-up instability is supercritical leading to a 'spotted-stripe' solution. Finally, we perform a similar analysis for the Klausmeier model of vegetation patterns on a steep hill, and examine transition from spots to stripes.This article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 2)'.

Spatially localized structures in the Gray-Scott model

Spatially localized structures in the one-dimensional Gray-Scott reaction-diffusion model are studied using a combination of numerical continuation techniques and weakly nonlinear theory, focusing on the regime in which the activator and substrate diffusivities are different but comparable. Localized states arise in three different ways: in a subcritical Turing instability present in this regime, and from folds in the branch of spatially periodic Turing states. They also arise from the fold of spatially uniform states. These three solution branches interconnect in complex ways. We use numerical continuation techniques to explore their global behaviour within a formulation of the model that has been used to describe dryland vegetation patterns on a flat terrain.This article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 2)'.

Multistability of model and real dryland ecosystems through spatial self-organization

Spatial self-organization of dryland vegetation constitutes one of the most promising indicators for an ecosystem's proximity to desertification. This insight is based on studies of reaction-diffusion models that reproduce visual characteristics of vegetation patterns observed on aerial photographs. However, until now, the development of reliable early warning systems has been hampered by the lack of more in-depth comparisons between model predictions and real ecosystem patterns. In this paper, we combined topographical data, (remotely sensed) optical data, and in situ biomass measurements from two sites in Somalia to generate a multilevel description of dryland vegetation patterns. We performed an in-depth comparison between these observed vegetation pattern characteristics and predictions made by the extended-Klausmeier model for dryland vegetation patterning. Consistent with model predictions, we found that for a given topography, there is multistability of ecosystem states with different pattern wavenumbers. Furthermore, observations corroborated model predictions regarding the relationships between pattern wavenumber, total biomass, and maximum biomass. In contrast, model predictions regarding the role of slope angles were not corroborated by the empirical data, suggesting that inclusion of small-scale topographical heterogeneity is a promising avenue for future model development. Our findings suggest that patterned dryland ecosystems may be more resilient to environmental change than previously anticipated, but this enhanced resilience crucially depends on the adaptive capacity of vegetation patterns.

A topographic mechanism for arcing of dryland vegetation bands

Banded patterns consisting of alternating bare soil and dense vegetation have been observed in water-limited ecosystems across the globe, often appearing along gently sloped terrain with the stripes aligned transverse to the elevation gradient. In many cases, these vegetation bands are arced, with field observations suggesting a link between the orientation of arcing relative to the grade and the curvature of the underlying terrain. We modify the water transport in the Klausmeier model of water-biomass interactions, originally posed on a uniform hillslope, to qualitatively capture the influence of terrain curvature on the vegetation patterns. Numerical simulations of this modified model indicate that the vegetation bands arc convex-downslope when growing on top of a ridge, and convex-upslope when growing in a valley. This behaviour is consistent with observations from remote sensing data that we present here. Model simulations show further that whether bands grow on ridges, valleys or both depends on the precipitation level. A survey of three banded vegetation sites, each with a different aridity level, indicates qualitatively similar behaviour.

Impact of Parameter Variability and Environmental Noise on the Klausmeier Model of Vegetation Pattern Formation

Semi-arid ecosystems made up of patterned vegetation, for instance, are thought to be highly sensitive. This highlights the importance of understanding the dynamics of the formation of vegetation patterns. The most renowned mathematical model describing such pattern formation consists of two partial differential equations and is often referred to as the Klausmeier model. This paper provides analytical and numerical investigations regarding the influence of different parameters, including the so-far not contemplated evaporation, on the long-term model results. Another focus is set on the influence of different initial conditions and on environmental noise, which has been added to the model. It is shown that patterning is beneficial for semi-arid ecosystems, that is, vegetation is present for a broader parameter range. Both parameter variability and environmental noise have only minor impacts on the model results. Increasing mortality has a high, nonlinear impact underlining the importance of further studies in order to gain a sufficient understanding allowing for suitable management strategies of this natural phenomenon.

Nonlocal grazing in patterned ecosystems

Many ecosystems exhibit gapped, labyrinthine, striped or spotted patterns. Important examples are vegetation patterns in drylands: these patterns are viewed as precursors of a catastrophic transition to a degraded state. A possible source of degradation is overgrazing, but many current spatially extended models include grazing in a local linear way. In this article nonlocal grazing responses are derived, taking into account (1) how many consumers there are (demographic response) (2) where they are (aggregative response) and (3) how much they forage (functional response). Different assumptions lead to different grazing responses, the type of grazing has a large influence on how ecosystems adapt to changing environmental conditions. In dryland simulations the different types of grazing are shown to alter the desertification process driven by decreasing rainfall. A sufficiently strong aggregative response leads to the suppression of vegetation patterns, nuancing their role as generic early warning signals.

Assessing the robustness of spatial pattern sequences in a dryland vegetation model

A particular sequence of patterns, 'gaps→labyrinth→spots', occurs with decreasing precipitation in previously reported numerical simulations of partial differential equation dryland vegetation models. These observations have led to the suggestion that this sequence of patterns can serve as an early indicator of desertification in some ecosystems. Because parameter values in the vegetation models can take on a range of plausible values, it is important to investigate whether the pattern sequence prediction is robust to variation. For a particular model, we find that a quantity calculated via bifurcation-theoretic analysis appears to serve as a proxy for the pattern sequences that occur in numerical simulations across a range of parameter values. We find in further analysis that the quantity takes on values consistent with the standard sequence in an ecologically relevant limit of the model parameter values. This suggests that the standard sequence is a robust prediction of the model, and we conclude by proposing a methodology for assessing the robustness of the standard sequence in other models and formulations.

When does colonisation of a semi-arid hillslope generate vegetation patterns?

Patterned vegetation occurs in many semi-arid regions of the world. Most previous studies have assumed that patterns form from a starting point of uniform vegetation, for example as a response to a decrease in mean annual rainfall. However an alternative possibility is that patterns are generated when bare ground is colonised. This paper investigates the conditions under which colonisation leads to patterning on sloping ground. The slope gradient plays an important role because of the downhill flow of rainwater. One long-established consequence of this is that patterns are organised into stripes running parallel to the contours; such patterns are known as banded vegetation or tiger bush. This paper shows that the slope also has an important effect on colonisation, since the uphill and downhill edges of an isolated vegetation patch have different dynamics. For the much-used Klausmeier model for semi-arid vegetation, the author shows that without a term representing water diffusion, colonisation always generates uniform vegetation rather than a pattern. However the combination of a sufficiently large water diffusion term and a sufficiently low slope gradient does lead to colonisation-induced patterning. The author goes on to consider colonisation in the Rietkerk model, which is also in widespread use: the same conclusions apply for this model provided that a small threshold is imposed on vegetation biomass, below which plant growth is set to zero. Since the two models are quite different mathematically, this suggests that the predictions are a consequence of the basic underlying assumption of water redistribution as the pattern generation mechanism.