ORIGINAL RESEARCH article

Front. Ecol. Evol., 10 November 2022

Sec. Models in Ecology and Evolution

Volume 10 - 2022 | https://doi.org/10.3389/fevo.2022.1068830

Contrasting effects of dispersal network heterogeneity on ecosystem stability in rock-paper-scissors games

  • 1. Research Center for Theoretical Ecology, Ministry of Education’s Key Laboratory of Poyang Lake Wetland and Watershed Research, School of Geography and Environment, Jiangxi Normal University, Nanchang, China

  • 2. School of Mathematics, Statistics and Actuarial Sciences, University of Kent, Canterbury, United Kingdom

Abstract

Intransitive competition, typically represented by the classic rock-paper-scissors game, provides an endogenous mechanism promoting species coexistence. As well known, species dispersal and interaction in nature might occur on complex patch networks, with species interacting in diverse ways. However, the effects of different interaction modes, combined with spatial heterogeneity in patch connectivities, have not been well integrated into our general understanding of how stable coexistence emerges in cyclic competition. We thus incorporate network heterogeneity into the classic rock-paper-scissors game, in order to compare ecosystem stability under two typical modes of interaction: species compete to fill empty sites, and species seize each other’s colony sites. On lattice-structured regular networks, the two interaction modes produce similar stability patterns through forming conspecific clusters to reduce interspecific competition. However, for heterogeneous networks, the interaction modes have contrasting effects on ecosystem stability. Specifically, if species compete for colony sites, increasing network heterogeneity stabilizes competitive dynamics. When species compete to fill empty sites, an increase in network heterogeneity leads to larger population fluctuations and therefore a higher risk of stochastic extinctions, in stark contrast to current knowledge. Our findings strongly suggest that particular attention should be devoted to testing which mode of interaction is more appropriate for modeling a given system.

Introduction

Understanding the mechanisms of ecosystem stability is a fundamental issue in ecology (May, 1972; ; Levine and HilleRisLambers, 2009; ). Many proposed mechanisms rely heavily on exogenous factors mitigating the effects of competition, such as niche differentiation (; Levine and HilleRisLambers, 2009; ) and intermediate disturbance hypothesis (; Roxburgh et al., 2004; Liao et al., 2022). In contrast to these exogenous mechanisms, intransitive competition provides an endogenous mechanism promoting coexistence (Laird and Schamp, 2006; ; Soliveres et al., 2015; Levine et al., 2017), and, in particular, the simple rock-paper-scissors game has become a typical case to explain ecosystem stability (; ; Reichenbach et al., 2007). The classic, cyclic game of rock-paper-scissors usually leads to species’ abundances neutrally cycling without converging to a stable equilibrium point, which is rarely observed in nature. To overcome this problem, many approaches have been proposed to explain the robust persistence of the cyclically competing system, such as higher-order interactions () and spatially structured interactions (; ; Rojas-Echenique and Allesina, 2011; ; Zhang et al., 2022). In particular, the latter mechanism of local structured interactions stabilizing coexistence has been documented experimentally (). Thus, it is widely believed that the inclusion of spatial structure, where the interactions and dispersal of individuals are local, can result in stable coexistence in rock-paper-scissors games. Yet, these conclusions are mostly drawn from lattice-based models, where each individual is assumed to only interact with its surrounding neighbors. This assumption is relatively restrictive, as species dispersal and interactions in nature might occur on complex networks with variation in patch connectivities (i.e., spatial heterogeneity in dispersal networks; Urban and Keitt, 2001; ; ; ; ; ; ; Li et al., 2021; Zhang et al., 2021).

There is abundant evidence that landscape structure, and other factors (e.g., patch quality; Liao et al., 2013), can result in anisotropic (i.e., directionally biased) dispersal behavior and therefore spatial heterogeneity in patch connectivities (Urban and Keitt, 2001; ; ; ; ; ; ). Since different patches in the landscape might be perceived differently by species (; ), and the resulting dispersal network would display diverse patterns of patch connectivity (Yeaton and Bond, 1991; ; Nicholson and Possingham, 2006; ; ; ; ). For instance, species dispersal between sub-reefs within the Great Barrier Reef has been characterized with scale-free networks (), while seed dispersal by birds, as opposed by winds, is better described by an irregular network than a spatially uniform network. In addition, identified a large spatial dispersal network of temporary ponds, which are used as breeding sites for amphibian species, following a power-law degree distribution. As such, there has been an increasing interest in exploring the effects of network heterogeneity on ecosystem stability using graph theory (Szabó et al., 2004; Szolnoki and Szabó, 2004; Masuda and Konno, 2006; ; Schütt and Claussen, 2010; ; ; Nagatani et al., 2018; ; ; Zhang et al., 2021). Many studies have found that increasing network heterogeneity (i.e., increasing variation in patch connectivities) can promote stable coexistence in cyclic competition (Masuda and Konno, 2006; Schütt and Claussen, 2010; Nagatani et al., 2018).

Despite these advances, several lattice-based models (Laird and Schamp, 2008; Rojas-Echenique and Allesina, 2011; Zhang et al., 2022) have observed that local intransitive competition can reduce species coexistence compared to long-range competition, in stark contrast to current knowledge of local interactions stabilizing coexistence (; ; ; ; ). Rojas-Echenique and Allesina (2011) and Zhang et al. (2022) attributed the opposite outcomes to different interaction modes, which can induce distinct stabilizing mechanisms in lattice-structured models. Indeed, species in diverse natural ecosystems might interact in different ways. For example, there are two typical interaction modes often observed in nature: seedlings of tree species or propagules of grass species compete to fill gaps, while animal species fight directly for colony sites (Rojas-Echenique and Allesina, 2011; ; Zhang et al., 2022). However, the effects of these different interaction modes, in combination with spatial heterogeneity in dispersal networks mentioned above, have not been well integrated into our general understanding of how stable coexistence emerges in cyclic competition. In this study, we thus incorporate dispersal network heterogeneity into the classic rock-paper-scissors games, in order to make a comparative analysis of ecosystem stability between the two typical interaction modes.

Materials and methods

Dispersal network heterogeneity

We consider a landscape consisting of a finite number (N) of patches (so-called colony sites), with each accommodating only one individual (or a subpopulation) of a species. In the landscape, individuals can move and interact between patches only along a predefined set of dispersal pathways. This shapes a dispersal network, with patches and dispersal pathways being represented by network nodes and links, respectively. To model the effect of spatial heterogeneity in dispersal networks (with fixed average patch degree 4), we generate four typical network structures with contrasting heterogeneities (i.e., the extent of variation in patch connectivities):

  • (1)

    A lattice-structured regular network with all patches having the same degree (illustrated in Figure 1A with each patch linked to other four patches).

  • (2)

    A randomly structured network with randomly connected patches (Watts and Strogatz, 1998), yielding small variation in patch degrees (Figure 1B).

  • (3)

    An exponential network constructed based on the algorithm of random attachment (). This produces a greater variation in patch connectivities than the random network (Figure 1C).

  • (4)

    A scale-free network structured according to the algorithm of preferential attachment (), producing the highest variation in patch connectivities (Figure 1D).

FIGURE 1

In these networks, species are assumed to use dispersal links in either direction without preference (i.e., undirected dispersal).

Two typical interaction modes

We consider a system of three cyclically competing species on a landscape of size N = 10,000 patches connected by dispersal networks. Each patch can only accommodate one individual of a species. The competitive relationships between these three species (i, j, and k) follow the rock-paper-scissors game, an example of intransitive competition which usually yields species coexistence with oscillations (; Li et al., 2020).

To focus solely on the effect of species competition on ecosystem stability, we perform simulations as follows: (i) according to the typical assumption in previous work (Masuda and Konno, 2006; Rojas-Echenique and Allesina, 2011; ; Nagatani et al., 2018; ; Zhang et al., 2022), initially all patches are populated with individuals randomly drawn from the three species; (ii) in each time step, we perform a competition event using the two interaction modes specified below; (iii) we repeat step (ii) for a long time, finding that 1,000 generations (1 generation = 10,000 time steps at N = 10,000 patches) are sufficient for the system to achieve steady state; (iv) at steady state, we record the number of individuals for each species and the spatial patterns at every generation.

In this study, we consider two interaction modes separately in the simulations.

Mode 1: Similar to and , in each time step, we randomly select a focal individual (with probability 1/N in the whole network) to die, and immediately choose two individuals randomly from its directly linked neighbors for pairwise competition, with the offspring of the winner occupying this empty patch. If the two individuals belong to the same species, then this species directly occupies the empty patch. If there is only one neighbor for the empty patch, then it is occupied directly by this neighbor.

Mode 2: Similar to Rojas-Echenique and Allesina (2011), in each time step, we randomly select two directly linked individuals for competition, with the offspring of the superior competitor substituting the inferior one (otherwise keeping the original state).

Under these two different interaction modes, we firstly explore how increasing network heterogeneity affects ecosystem stability in classic rock-paper-scissors games, and then test whether these outcomes are robust to varying network size.

Results

We begin our analysis by inspecting the temporal evolution of species abundances (Figure 2). When species compete to fill an empty patch (Mode 1), increasing network heterogeneity (from lattice-structured regular to scale-free networks) increases the magnitude of population fluctuations and speeds up stochastic species extinctions. In contrast, if the superior competitor directly substitutes the inferior one (Mode 2), all three species can coexist, and the magnitude of population fluctuations around the equilibrium point (≈1/3) decreases as network heterogeneity increases. In lattice-structured regular networks, both interaction modes yield similar fluctuation size (Figure 2A vs. Figure 2E with different ranges of y-axis). This is more clear when looking at their dynamic trajectories (Figure 2I vs. Figure 2M). In these dynamic trajectories (Figures 2I–P), the abundances of the species at each generation represent a point in the 3-simplex, whose vertices correspond to a monospecific population. Over time, the point follows a trajectory which will reflect the state of the system. If the system is close to the equilibrium abundances, the trajectory eventually occupies a small area near it, whilst larger fluctuations cover larger areas. Similar to Figures 2A–H, the interaction mode 1 results in greater fluctuations and ultimately monoculture in more heterogeneous networks (Figures 2I–L), as opposed to the interaction mode 2 where ecosystems are more stable in networks with higher heterogeneity.

FIGURE 2

To test whether these outcomes are robust, we simulate 20 replicates for each case, by regenerating the dispersal network in each replicate (Figure 3). In Mode 1, three species can coexist with oscillations around the equilibrium point in the lattice-structured regular network (see Figure 2A), unlike other three heterogeneous networks where only one species ultimately dominates the system (Figure 3A). As such, the outcome in the regular network is absent in Figure 3A. In other three heterogeneous networks, we record the extinction time of the first species in each replicate, with more time required for species extinction implying more stable coexistence. As shown in Figure 3A, increasing network heterogeneity generally leads to species going extinct sooner. In Mode 2, all species can coexist with fluctuations around the equilibrium point in all dispersal networks. Thus, we use variation in fluctuations size (C.V. around the mean species abundances at steady state) to characterize ecosystem stability, with lower values of C.V. indicating more stable coexistence. As shown in Figure 3B, an increase in network heterogeneity generally results in smaller fluctuations in species abundances, yielding a more stable system. These outcomes further confirm our previous conclusions in Figure 2.

FIGURE 3

To better understand the mechanisms behind these dynamic behaviors, we display several snapshots of the spatial patterns of the three competitors in complex networks (Figure 4). In the lattice-structured regular network, both interaction modes produce similar conspecific clumping patterns (Figures 4A,B). To clearly illustrate the spatial organization in dispersal networks with contrasting heterogeneities, we take a system of N = 100 patches for example by displaying their snapshots at the 5th generation in Mode 1 and the 20th generation in Mode 2 (Figures 4C–J). In Mode 1 (Figures 4C–F), when three species have almost equal abundance in the regular network, increasing network heterogeneity which tends to increase network modularity, enlarges differences in species abundances, most likely accelerating stochastic extinctions. In Mode 2 (Figures 4G–J), if three species have unequal abundances in the regular network, dispersal network heterogeneity (i.e., network modularity) acts as a driving force to equalize their abundances, thereby stabilizing the ecosystem. We also observe that species can form self-organized conspecific clusters with the most connected patches at the core in these heterogeneous networks.

FIGURE 4

Under both interaction modes, we finally focus on how system size (i.e., increasing network size) affects fluctuations in species abundance in networks with contrasting heterogeneities (Figure 5). Intuitively, increasing network size promotes ecosystem stability, regardless of interaction mode and network heterogeneity. Specifically, stochastic extinctions in Mode 1 require more time to occur as network size increases, especially in networks with less heterogeneity (Figure 5A). This also demonstrates that increasing network heterogeneity greatly destabilizes the system. In Mode 2, the size of population fluctuation declines linearly with increasing network size (with log-log scale in Figure 5B), thereby promoting ecosystem stability. Furthermore, increasing network heterogeneity leads to smaller fluctuations, further confirming that network heterogeneity can stabilize the cyclically competing system.

FIGURE 5

Discussion

Incorporating complex networks into the classic rock-paper-scissors game, we find that whether network heterogeneity can stabilize competitive dynamics depends on the interaction mode. Specifically, if species compete directly for colony sites, increasing spatial heterogeneity in dispersal networks stabilizes the cyclically competing ecosystems, further confirming previous theoretical arguments (Masuda and Konno, 2006; Schütt and Claussen, 2010; Nagatani et al., 2018). In contrast, when species compete to fill an empty site, an increase in network heterogeneity leads to stronger population fluctuations and thus increases the risk of stochastic extinctions. Interestingly, in the lattice-structured regular network, both interaction modes display similar coexistence patterns based on the same mechanism: local interactions allow species to survive by forming conspecific clusters (via self-organization) where interspecific competition only takes place at the borders between heterospecific clusters, thereby decreasing the effective interspecific competition (by reducing interspecific encounter rate) and therefore population fluctuations (Figures 4A,B).

There exist different mechanisms between the two interaction modes that induce contrasting effects of network heterogeneity on ecosystem stability. In our model, the stabilizing role of network heterogeneity is observed when competition occurs through the superior competitor displacing the inferior one in a randomly chosen pair of neighbors. This competition mode has been commonly used and is symmetric in a cyclical way. However, increasing variation in patch connectivities implies that a few patches are highly connected while most have only few connections. A patch with more connections is more likely to be selected as a neighbor, while another neighbor can be sampled from multiple patches that might accommodate any of the three species. As such, which species can occupy this highly connected patch will change more frequently in networks with higher heterogeneity. To some extent, this prohibits the growth of self-organized conspecific clusters with the most connected patches at the core. In fact, increasing network heterogeneity breaks the symmetry in species interactions, that is, those individuals in highly connected patches have a higher chance to participate in competition than other individuals in poorly connected patches. This results in asymmetric interactions that can produce negative frequency dependence (NFD) which is absent in symmetric interactions (Rojas-Echenique and Allesina, 2011). Specifically, the NFD can decrease the average fitness of individuals when they become more common (cf. Zhang et al., 2022). Thus, the NFD, which is more significant in more heterogeneous networks, can suppress the population growth of dominant species but promote the growth of rare species, thereby increasing the frequency of oscillations and stabilizing competitive dynamics.

Network heterogeneity is destabilizing for the interaction mode where species compete for an empty site. This is because, in this mode, each species has the same probability to be selected to die, regardless of patch connectivities, so NFD cannot emerge. As such, conspecific clusters, with highly connected patches at the core, can grow (via self-organization) without restriction especially in networks with higher heterogeneity. Within these conspecific clusters, a competition event does not contribute to variation in species abundance, as local competition occurs between conspecific individuals. Thus, changes in species abundances can only take place along the borders between heterospecific clusters. However, the length of these borders increases more rapidly in more heterogeneous networks (i.e., a rapid increase in the number of “active” individuals that can change the state at the borders), promoting interspecific encounter rate and therefore interspecific competition. In addition, these heterospecific clusters can be treated as compartments dominated by different species, but there are lots of links connecting these compartments. If a superior competitor invades a compartment dominated by an inferior one, then it can colonize this compartment via the highly connected core. Ultimately, this results in larger population fluctuations and makes stochastic extinction events more likely.

This study incorporates complex networks into the rock-paper-scissors game under two commonly observed modes of interaction. In the lattice-structured regular network, both interaction modes yield similar system stability through forming conspecific clusters to reduce interspecific competition (Figure 4A vs. Figure 4B). Interestingly, we find that the inclusion of network heterogeneity can induce contrasting coexistence patterns between the two interaction modes, due to different mechanisms as explained above. This strongly suggests that particular attention should be devoted to testing, theoretically and experimentally, which mode of interaction is more appropriate for modeling a given competing system. For example, whereas the death of a tree or grass creates a gap in plant communities (i.e., Mode 1), animal species more typically fight for territory (i.e., Mode 2). Thus, these communities should be modeled with different interaction modes and so can be expected to be affected differently by the structure of the landscape they inhabit. Overall, even if these outcomes are obtained using an extremely simplified model, our findings can help identify different mechanisms to explain the role of different interaction modes in stabilizing competitive dynamics in complex networks with contrasting heterogeneities.

Statements

Data availability statement

The original contributions presented in this study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

GG, ZZ, and HZ built the model, performed simulations, and analyzed the results. JL conceived the study and wrote the manuscript. DB contributed substantially to revisions. All authors contributed to the article and approved the submitted version.

Funding

JL was supported by the National Natural Science Foundation of China (Nos. 32271548 and 31901175).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Summary

Keywords

competitive intransitivity, cyclically competing ecosystems, dispersal network heterogeneity, ecosystem stability, rock-paper-scissors games

Citation

Guo G, Zhang Z, Zhang H, Bearup D and Liao J (2022) Contrasting effects of dispersal network heterogeneity on ecosystem stability in rock-paper-scissors games. Front. Ecol. Evol. 10:1068830. doi: 10.3389/fevo.2022.1068830

Received

13 October 2022

Accepted

28 October 2022

Published

10 November 2022

Volume

10 - 2022

Edited by

Rui-Wu Wang, Northwestern Polytechnical University, China

Reviewed by

Min Su, Hefei University of Technology, China; Xiaojie Chen, University of Electronic Science and Technology of China, China

Updates

Copyright

*Correspondence: Jinbao Liao,

This article was submitted to Models in Ecology and Evolution, a section of the journal Frontiers in Ecology and Evolution

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics