How collectively integrated are ecological communities?

Beyond abiotic conditions, do population dynamics mostly depend on a species' direct predators, preys and conspecifics? Or can indirect feedback that ripples across the whole community be equally important? Determining where ecological communities sit on the spectrum between these two characterizations requires a metric able to capture the difference between them. Here we show that the spectral radius of a community's interaction matrix provides such a metric, thus a measure of ecological collectivity, which is accessible from imperfect knowledge of biotic interactions and related to observable signatures. This measure of collectivity integrates existing approaches to complexity, interaction structure and indirect interactions. Our work thus provides an original perspective on the question of to what degree communities are more than loose collections of species or simple interaction motifs and explains when pragmatic reductionist approaches ought to suffice or fail when applied to ecological communities.

The coordination of green-brown food webs and their disruption by anthropogenic nutrient inputs.

Aim: Our goal was to quantify nitrogen flows and stocks in green-brown food webs in different ecosystems, how they differ across ecosystems and how they respond to nutrient enrichment. Location: Global. Time period: Contemporary. Major taxa studied: Plants, phytoplankton, macroalgae, invertebrates, vertebrates and zooplankton. Methods: Data from >500 studies were combined to estimate nitrogen stocks and fluxes in green-brown food webs in forests, grasslands, brackish environments, seagrass meadows, lakes and oceans. We compared the stocks, fluxes and metabolic rates of different functional groups within each food web. We also used these estimates to build a dynamical model to test the response of the ecosystems to nutrient enrichment. Results: We found surprising symmetries between the green and brown channels across ecosystems, in their stocks, fluxes and consumption coefficients and mortality rates. We also found that nitrogen enrichment, either organic or inorganic, can disrupt this balance between the green and brown channels. Main conclusions: Linking green and brown food webs reveals a previously hidden symmetry between herbivory and detritivory, which appears to be a widespread property of natural ecosystems but can be disrupted by anthropogenic nitrogen additions.

High-integrity human intervention in ecosystems: Tracking self-organization modes.

Humans play major roles in shaping and transforming the ecology of Earth. Unlike natural drivers of ecosystem change, which are erratic and unpredictable, human intervention in ecosystems generally involves planning and management, but often results in detrimental outcomes. Using model studies and aerial-image analysis, we argue that the design of a successful human intervention form calls for the identification of the self-organization modes that drive ecosystem change, and for studying their dynamics. We demonstrate this approach with two examples: grazing management in drought-prone ecosystems, and rehabilitation of degraded vegetation by water harvesting. We show that grazing can increase the resilience to droughts, rather than imposing an additional stress, if managed in a spatially non-uniform manner, and that fragmental restoration along contour bunds is more resilient than the common practice of continuous restoration in vegetation stripes. We conclude by discussing the need for additional studies of self-organization modes and their dynamics.

General statistical scaling laws for stability in ecological systems.

Ecological stability refers to a family of concepts used to describe how systems of interacting species vary through time and respond to disturbances. Because observed ecological stability depends on sampling scales and environmental context, it is notoriously difficult to compare measurements across sites and systems. Here, we apply stochastic dynamical systems theory to derive general statistical scaling relationships across time, space, and ecological level of organisation for three fundamental stability aspects: resilience, resistance, and invariance. These relationships can be calibrated using random or representative samples measured at individual scales, and projected to predict average stability at other scales across a wide range of contexts. Moreover deviations between observed vs. extrapolated scaling relationships can reveal information about unobserved heterogeneity across time, space, or species. We anticipate that these methods will be useful for cross-study synthesis of stability data, extrapolating measurements to unobserved scales, and identifying underlying causes and consequences of heterogeneity.

Scaling up biodiversity-ecosystem functioning relationships: the role of environmental heterogeneity in space and time.

The biodiversity and ecosystem functioning (BEF) relationship is expected to be scale-dependent. The autocorrelation of environmental heterogeneity is hypothesized to explain this scale dependence because it influences how quickly biodiversity accumulates over space or time. However, this link has yet to be demonstrated in a formal model. Here, we use a Lotka-Volterra competition model to simulate community dynamics when environmental conditions vary across either space or time. Species differ in their optimal environmental conditions, which results in turnover in community composition. We vary biodiversity by modelling communities with different sized regional species pools and ask how the amount of biomass per unit area depends on the number of species present, and the spatial or temporal scale at which it is measured. We find that more biodiversity is required to maintain functioning at larger temporal and spatial scales. The number of species required increases quickly when environmental autocorrelation is low, and slowly when autocorrelation is high. Both spatial and temporal environmental heterogeneity lead to scale dependence in BEF, but autocorrelation has larger impacts when environmental change is temporal. These findings show how the biodiversity required to maintain functioning is expected to increase over space and time.

Thermal mismatches in biological rates determine trophic control and biomass distribution under warming.

Temperature has numerous effects on the structure and dynamics of ecological communities. Yet, there is no general trend or consensus on the magnitude and directions of these effects. To fill this gap, we propose a mechanistic framework based on key biological rates that predicts how temperature influences biomass distribution and trophic control in food webs. We show that these predictions arise from thermal mismatches between biological rates and across trophic levels. We couple our theory with experimental data for a wide range of species and find that warming should lead to top-heavier terrestrial food chains and stronger top-down control in aquatic environments. We then derive predictions for the effects of temperature on herbivory and validate them with data on stream grazers. Our study provides a mechanistic explanation of thermal effects on consumer-resource systems which is crucial to better understand the biogeography and the consequences of global warming on trophic dynamics.

Scaling-up biodiversity-ecosystem functioning research.

A rich body of knowledge links biodiversity to ecosystem functioning (BEF), but it is primarily focused on small scales. We review the current theory and identify six expectations for scale dependence in the BEF relationship: (1) a nonlinear change in the slope of the BEF relationship with spatial scale; (2) a scale-dependent relationship between ecosystem stability and spatial extent; (3) coexistence within and among sites will result in a positive BEF relationship at larger scales; (4) temporal autocorrelation in environmental variability affects species turnover and thus the change in BEF slope with scale; (5) connectivity in metacommunities generates nonlinear BEF and stability relationships by affecting population  synchrony at local and regional scales; (6) spatial scaling in food web structure and diversity will generate scale dependence in ecosystem functioning. We suggest directions for synthesis that combine approaches in metaecosystem and metacommunity ecology and integrate cross-scale feedbacks. Tests of this theory may combine remote sensing with a generation of networked experiments that assess effects at multiple scales. We also show how anthropogenic land cover change may alter the scaling of the BEF relationship. New research on the role of scale in BEF will guide policy linking the goals of managing biodiversity and ecosystems.

The three regimes of spatial recovery.

An enduring challenge for ecology is identifying the drivers of ecosystem and population stability. In a spatially explicit context, key features to consider are landscape spatial structure, local interactions, and dispersal. Substantial work has been done on each of these features as a driver of stability, but little is known on the interplay between them. Missing has been a more integrative approach, able to map and identify different dynamical regimes, predicting a system's response to perturbations. Here we first consider a simple scenario, i.e., the recovery of a homogeneous metapopulation from a single localized pulse disturbance. The analysis of this scenario reveals three fundamental recovery regimes: Isolated Regime when dispersal is not significant, Rescue Regime when dispersal mediates recovery, and Mixing Regime when perturbations spread throughout the system. Despite its simplicity, our approach leads to remarkably general predictions. These include the qualitatively different outcomes of various scenarios of habitat fragmentation, the surprising benefits of local extinctions on population persistence at the transition between regimes, and the productivity shifts of metacommunities in a changing environment. This study thus provides context to known results and insight into future directions of research.

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)'.

Implications of tristability in pattern-forming ecosystems.

Many ecosystems show both self-organized spatial patterns and multistability of possible states. The combination of these two phenomena in different forms has a significant impact on the behavior of ecosystems in changing environments. One notable case is connected to tristability of two distinct uniform states together with patterned states, which has recently been found in model studies of dryland ecosystems. Using a simple model, we determine the extent of tristability in parameter space, explore its effects on the system dynamics, and consider its implications for state transitions or regime shifts. We analyze the bifurcation structure of model solutions that describe uniform states, periodic patterns, and hybrid states between the former two. We map out the parameter space where these states exist, and note how the different states interact with each other. We further focus on two special implications with ecological significance, breakdown of the snaking range and complex fronts. We find that the organization of the hybrid states within a homoclinic snaking structure breaks down as it meets a Maxwell point where simple fronts are stationary. We also discover a new series of complex fronts between the uniform states, each with its own velocity. We conclude with a brief discussion of the significance of these findings for the dynamics of regime shifts and their potential control.

The Impact of Spatial and Temporal Dimensions of Disturbances on Ecosystem Stability.

Ecosystems constantly face disturbances which vary in their spatial and temporal features, yet little is known on how these features affect ecosystem recovery and persistence, i.e., ecosystem stability. We address this issue by considering three ecosystem models with different local dynamics, and ask how their stability properties depend on the spatial and temporal properties of disturbances. We measure the spatial dimension of disturbances by their spatial extent while controlling for their overall strength, and their temporal dimension by the average frequency of random disturbance events. Our models show that the return to equilibrium following a disturbance depends strongly on the disturbance's extent, due to rescue effects mediated by dispersal. We then reveal a direct relation between the temporal variability caused by repeated disturbances and the recovery from an isolated disturbance event. Although this could suggest a trivial dependency of ecosystem response on disturbance frequency, we find that this is true only up to a frequency threshold, which depends on both the disturbance spatial features and the ecosystem dynamics. Beyond this threshold the response changes qualitatively, displaying spatial clusters of disturbed regions, causing an increase in variability, and even a system-wide collapse for ecosystems with alternative stable states. Thus, spanning the spatial dimension of disturbances is a way to probe the underlying dynamics of an ecosystem. Furthermore, considering spatial and temporal dimensions of disturbances in conjunction is necessary to predict ecosystem responses with dramatic ecological consequences, such as regime shifts or population extinction.

Asymmetric disease dynamics in multihost interconnected networks.

Epidemic spread in single-host systems strongly depends on the population's transmission network. However, little is known regarding the spread of epidemics across networks representing populations of multiple hosts. We explored cross-species transmission in a multilayer network where layers represent populations of two distinct hosts, and disease can spread across intralayer (within-host) and interlayer (between-host) edges. We developed an analytic framework for the SIR epidemic model to examine the effect of (i) source of infection and (ii) between-host asymmetry in infection probabilities, on disease risk. We measured risk as outbreak probability and outbreak size in a focal host, represented by one network layer. Numeric simulations were used to validate the analytic formulations. We found that outbreak probability is determined by a complex interaction between source of infection and between-host infection probabilities, whereas outbreak size is mainly affected by the non-focal host to focal host infection probability. Hence, inter-specific asymmetry in infection probabilities shapes disease dynamics in multihost networks. These results highlight the importance of considering multiple measures of disease risk and advance our understanding of disease spread in multihost systems. The study provides a flexible way to model disease dynamics in multiple hosts while considering contact heterogeneity within and between species. We strongly encourage empirical studies that include information on both cross-species infection rates and network structure of multiple hosts. Such studies are necessary to corroborate our theoretical results and to improve our understanding of multihost epidemiology.

Desertification by front propagation?

Understanding how desertification takes place in different ecosystems is an important step in attempting to forecast and prevent such transitions. Dryland ecosystems often exhibit patchy vegetation, which has been shown to be an important factor on the possible regime shifts that occur in arid regions in several model studies. In particular, both gradual shifts that occur by front propagation, and abrupt shifts where patches of vegetation vanish at once, are a possibility in dryland ecosystems due to their emergent spatial heterogeneity. However, recent theoretical work has suggested that the final step of desertification - the transition from spotted vegetation to bare soil - occurs only as an abrupt shift, but the generality of this result, and its underlying origin, remain unclear. We investigate two models that detail the dynamics of dryland vegetation using a markedly different functional structure, and find that in both models the final step of desertification can only be abrupt. Using a careful numerical analysis, we show that this behavior is associated with the disappearance of confined spot-pattern domains as stationary states, and identify the mathematical origin of this behavior. Our findings show that a gradual desertification to bare soil due to a front propagation process can not occur in these and similar models, and opens the question of whether these dynamics can take place in nature.

Discovery of fairy circles in Australia supports self-organization theory.

Vegetation gap patterns in arid grasslands, such as the "fairy circles" of Namibia, are one of nature's greatest mysteries and subject to a lively debate on their origin. They are characterized by small-scale hexagonal ordering of circular bare-soil gaps that persists uniformly in the landscape scale to form a homogeneous distribution. Pattern-formation theory predicts that such highly ordered gap patterns should be found also in other water-limited systems across the globe, even if the mechanisms of their formation are different. Here we report that so far unknown fairy circles with the same spatial structure exist 10,000 km away from Namibia in the remote outback of Australia. Combining fieldwork, remote sensing, spatial pattern analysis, and process-based mathematical modeling, we demonstrate that these patterns emerge by self-organization, with no correlation with termite activity; the driving mechanism is a positive biomass-water feedback associated with water runoff and biomass-dependent infiltration rates. The remarkable match between the patterns of Australian and Namibian fairy circles and model results indicate that both patterns emerge from a nonuniform stationary instability, supporting a central universality principle of pattern-formation theory. Applied to the context of dryland vegetation, this principle predicts that different systems that go through the same instability type will show similar vegetation patterns even if the feedback mechanisms and resulting soil-water distributions are different, as we indeed found by comparing the Australian and the Namibian fairy-circle ecosystems. These results suggest that biomass-water feedbacks and resultant vegetation gap patterns are likely more common in remote drylands than is currently known.

Gradual regime shifts in fairy circles.

Large responses of ecosystems to small changes in the conditions--regime shifts--are of great interest and importance. In spatially extended ecosystems, these shifts may be local or global. Using empirical data and mathematical modeling, we investigated the dynamics of the Namibian fairy circle ecosystem as a case study of regime shifts in a pattern-forming ecosystem. Our results provide new support, based on the dynamics of the ecosystem, for the view of fairy circles as a self-organization phenomenon driven by water-vegetation interactions. The study further suggests that fairy circle birth and death processes correspond to spatially confined transitions between alternative stable states. Cascades of such transitions, possible in various pattern-forming systems, result in gradual rather than abrupt regime shifts.

Species survival emerge from rare events of individual migration.

Ecosystems greatly vary in their species composition and interactions, yet they all show remarkable resilience to external influences. Recent experiments have highlighted the significant effects of spatial structure and connectivity on the extinction and survival of species. It has also been emphasized lately that in order to study extinction dynamics reliably, it is essential to incorporate stochasticity, and in particular the discrete nature of populations, into the model. Accordingly, we applied a bottom-up modeling approach that includes both spatial features and stochastic interactions to study survival mechanisms of species. Using the simplest spatial extension of the Lotka-Volterra predator-prey model with competition, subject to demographic and environmental noise, we were able to systematically study emergent properties of this rich system. By scanning the relevant parameter space, we show that both survival and extinction processes often result from a combination of habitat fragmentation and individual rare events of recolonization.

Interplay between Turing mechanisms can increase pattern diversity.

We use the context of dryland vegetation to study a general problem of complex pattern-forming systems: multiple pattern-forming instabilities that are driven by distinct mechanisms but share the same spectral properties. We find that the co-occurrence of two Turing instabilities when the driving mechanisms counteract each other in some region of the parameter space results in the growth of a single mode rather than two interacting modes. The interplay between the two mechanisms compensates for the simpler dynamics of a single mode by inducing a wider variety of patterns, which implies higher biodiversity in dryland ecosystems.

Regime shifts in models of dryland vegetation.

Drylands are pattern-forming systems showing self-organized vegetation patchiness, multiplicity of stable states and fronts separating domains of alternative stable states. Pattern dynamics, induced by droughts or disturbances, can result in desertification shifts from patterned vegetation to bare soil. Pattern formation theory suggests various scenarios for such dynamics: an abrupt global shift involving a fast collapse to bare soil, a gradual global shift involving the expansion and coalescence of bare-soil domains and an incipient shift to a hybrid state consisting of stationary bare-soil domains in an otherwise periodic pattern. Using models of dryland vegetation, we address the question of which of these scenarios can be realized. We found that the models can be split into two groups: models that exhibit multiplicity of periodic-pattern and bare-soil states, and models that exhibit, in addition, multiplicity of hybrid states. Furthermore, in all models, we could not identify parameter regimes in which bare-soil domains expand into vegetated domains. The significance of these findings is that, while models belonging to the first group can only exhibit abrupt shifts, models belonging to the second group can also exhibit gradual and incipient shifts. A discussion of open problems concludes the paper.

Regime shifts in models of dryland vegetation.

Drylands are pattern-forming systems showing self-organized vegetation patchiness, multiplicity of stable states and fronts separating domains of alternative stable states. Pattern dynamics, induced by droughts or disturbances, can result in desertification shifts from patterned vegetation to bare soil. Pattern formation theory suggests various scenarios for such dynamics: an abrupt global shift involving a fast collapse to bare soil, a gradual global shift involving the expansion and coalescence of bare-soil domains and an incipient shift to a hybrid state consisting of stationary bare-soil domains in an otherwise periodic pattern. Using models of dryland vegetation, we address the question of which of these scenarios can be realized. We found that the models can be split into two groups: models that exhibit multiplicity of periodic-pattern and bare-soil states, and models that exhibit, in addition, multiplicity of hybrid states. Furthermore, in all models, we could not identify parameter regimes in which bare-soil domains expand into vegetated domains. The significance of these findings is that, while models belonging to the first group can only exhibit abrupt shifts, models belonging to the second group can also exhibit gradual and incipient shifts. A discussion of open problems concludes the paper.