Finding all multiple stable fixpoints of \(n\)-species Lotka-Volterra competition models.

*(English)*Zbl 1381.92084Summary: One way to explore assembly of extant and novel communities from species pools, and by that biodiversity and species ranges, is to study the equilibrium behavior of dynamic competition models such as the Lotka-Volterra competition (LVC) model. We present a novel method (COMMUSTIX) to determine all stable fixpoints of the general LVC model with abundances \(x\) from a given pool of \(n\) species. To that purpose, we split the species in potentially surviving species (\(x_{\mathrm{i}} > 0\)) and in others going extinct (\(x_{\mathrm{i}} = 0\)). We derived criteria for the stability of \(x_{\mathrm{i}} = 0\) and for the equilibrium of \(x_{\mathrm{i}} > 0\) to determine possible combinations of extinct and surviving species by iteratively applying a mixed binary linear optimization algorithm. We tested this new method against (a) the numerical solution at equilibrium of the LVC ordinary differential equations (ODEs) and (b) the fixpoints of all combinations of surviving and extinct species (possible only for small \(n\)), tested for stability and non-negativity. The tests revealed that COMMUSTIX is reliable, it detects all multiple stable fixpoints (SFPs), which is not guaranteed by solving the ODEs, and more efficient than the combinations method. With COMMUSTIX, we studied the dependence of the fixpoint behavior on the competition strengths relative to the intra-specific competition. If inter-specific competition was considerably lower than intra-specific competition, only globally SFPs occurred. In contrast, if all inter-specific was higher than intra-specific competition, multiple SFPs consisting of only one species occurred. If competition strengths in the species pool ranged from below to above the intra-specific competition, either global or multiple SFPs strongly differing in species composition occurred. The species richness over all SFPs was high for pools of species with similar, either weak or strong competition, and lower for species with dissimilar or close to intra-specific competition strengths. The new approach is a reliable and efficient tool for further extensive examinations of the dependence of community compositions on parameter settings of the LVC model.

##### MSC:

92D25 | Population dynamics (general) |

34D20 | Stability of solutions to ordinary differential equations |

90C05 | Linear programming |

##### Keywords:

Lotka-Volterra competition model; multiple stable fixpoints; stability; non-negativity; mixed binary linear programming; relative competition strength
PDF
BibTeX
Cite

\textit{H. Lischke} and \textit{T. J. Löffler}, Theor. Popul. Biol. 115, 24--34 (2017; Zbl 1381.92084)

Full Text:
DOI

##### References:

[1] | Allesina, S.; Tang, S., Stability criteria for complex ecosystems, Nature, 483, 205-208, (2012) |

[2] | Brauer, A., Limits for the characteristic roots of a matrix.7, Duke Math. J., 25, 583-600, (1958) |

[3] | Carrara, F., Experimental evidence for strong stabilizing forces at high functional diversity of aquatic microbial communities, Ecology, 96, 1340-1350, (2015) |

[4] | Carrara, F., Inferring species interactions in ecological communities: a comparison of methods at different levels of complexity, Methods Ecol. Evol., (2015) |

[5] | Chase, J. M., Stochastic community assembly causes higher biodiversity in more productive environments, Science, 328, 1388-1391, (2010) |

[6] | Chesson, P., Mechanisms of maintenance of species diversity, Ann. Rev. Ecol. Syst., 31, 343-366, (2000) |

[7] | Corrada Bravo, H., et al. 2015. Package ‘Rcplex’: R interface to CPLEX. |

[8] | DeOca, F. M.; Zeeman, M. L., Extinction in nonautonomous competitive Lotka-Volterra systems, Proc. Amer. Math. Soc., 124, 3677-3687, (1996) · Zbl 0866.34029 |

[9] | Fukami, T., Community assembly along a species pool gradient: implications for multiple-scale patterns of species diversity, Popul. Ecol., 46, 137-147, (2004) |

[10] | Fukami, T., 2015. Historical contingency in community assembly: integrating niches, species pools, and priority effects. In: D. J. Futuyma, (Ed.), Annual Review of Ecology, Evolution, and Systematics, Vol 46. Annual Reviews, Palo Alto, pp. 1-23. |

[11] | Gao, Q., Scaling simulation models for spatially heterogeneous ecosystems with diffusive transportation, Landsc. Ecol., 16, 289, (2001) |

[12] | Garren, K. R., Bounds for the eigenvalues of a matrix. NASA TECHNICALNOTE, vol. NASA TN D-4373, (1968), National Aeromautics and Space Administration, Langley Research Center |

[13] | Gilpin, M. E.; Case, T. J., Multiple domains of attraction in competition communities, Nature, 261, 40-42, (1976) |

[14] | Giorgi, G., Zuccotti, C., 2015. An overview on D-stable matrices. DEM Working Paper Series, Vol. # 97 (02-15) Universita di Pavia Department of Economics and Management. |

[15] | Godsoe, W., The effect of competition on species’ distributions depends on coexistence, rather than scale alone, Ecography, 38, 1071-1079, (2015) |

[16] | Goh, B. S., Sector stability of a complex ecosystem model, Math. Biosci., 40, 157-166, (1978) · Zbl 0386.92013 |

[17] | Goh, B. S.; Jennings, L. S., Feasibility and stability in randomly assembled Lotka-Volterra models, Ecol. Modell., 3, 63-71, (1977) |

[18] | Gravel, D., Stability and complexity in model meta-ecosystems, Nature Commun., 7, (2016) |

[19] | Griva, I., Linear and nonlinear optimization, (2009), Society for Industrial and Applied Mathematics Philadelphia |

[20] | Hastings, A.; Wysham, D. B., Regime shifts in ecological systems can occur with no warning, Ecol. Lett., 13, 464-472, (2010) |

[21] | Hofbauer, J.; Siegmund, K., Stability of \(n\)-species communities, (Hofbauer, J.; Siegmund, K., Evolutionary Games and Population Dynamics, (1998), Cambridge University Press Cambridge), 181-202 |

[22] | Horn, R. A.; Johnson, C. R., Matrix analysis, (2012), Cambridge University Press |

[23] | Kallrath, J., Gemischt-ganzzahlige optimierung: modellierung in der praxis. mit fallstudien aus chemie, energiewirtschaft, papierindustrie, metallgewerbe, produktion und logistik, (2013), Springer Spektrum Wiesbaden · Zbl 1261.90001 |

[24] | Koch, T., MIPLIB 2010: mixed integer programming library version 5, Math. Program. Comput., 3, 103-163, (2011) |

[25] | Kunstler, G., Plant functional traits have globally consistent effects on competition, Nature, 529, 204-207, (2016) |

[26] | Law, R.; Morton, R. D., Permanence and the assembly of ecological communities, Ecology, 77, 762-775, (1996) |

[27] | Li, X. Z., The criteria for globally stable equilibrium in \(n\)-dimensional Lotka-Volterra systems, J. Math. Anal. Appl., 240, 600-606, (1999) · Zbl 0947.34044 |

[28] | Logofet, D. O., Stronger-than-Lyapunov notions of matrix stability, or how “flowers” help solve problems in mathematical ecology, Linear Algebra Appl., 398, 75-100, (2005) · Zbl 1062.92072 |

[29] | May, R. M., Will a large complex system be stable?, Nature, 238, 413-414, (1972) |

[30] | May, R. M., Thresholds and breakpoints in ecosystems with a multiplicity of stable states, Nature, 269, 471-477, (1977) |

[31] | Moore, G. E., Cramming more components onto integrated circuits, Proc. IEEE, 86, 82-85, (1998), Reprinted from Electronics, pg 114-117, April 19, 1965 |

[32] | Neilson, R. P., Forecasting regional to global plant migration in response to climate change, Bioscience, 55, 749-759, (2005) |

[33] | Petraitis, P., Multiple stable states in natural ecosystems, (2013), Oxford University Press USA |

[34] | Poisot, T., Beyond species: why ecological interaction networks vary through space and time, Oikos, 124, 243-251, (2015) |

[35] | Pu, Z.; Jiang, L., Dispersal among local communities does not reduce historical contingencies during metacommunity assembly, Oikos, 124, 1327-1336, (2015) |

[36] | Ricklefs, R. E., Intrinsic dynamics of the regional community, Ecol. Lett., 18, 497-503, (2015) |

[37] | Roberts, A., Stability of a feasible random ecosystem, Nature, 251, 607-608, (1974) |

[38] | Rodriguez, J., Structural continuity and multiple alternative stable states in middle pleistocene european Mammalian communities, Palaeogeogr. Palaeoclimatol. Palaeoecol., 239, 355-373, (2006) |

[39] | Rohr, R. P., On the structural stability of mutualistic systems, Science, 345, 1253497-1-1253497-9, (2014) |

[40] | Scheffer, M., Catastrophic shifts in ecosystems, Nature, 413, 591-596, (2001) |

[41] | Scheffer, M., Early-warning signals for critical transitions, Nature, 461, 53-59, (2009) |

[42] | Scheffer, M., Thresholds for boreal biome transitions, Proc. Natl. Acad. Sci., 109, 21384-21389, (2012) |

[43] | Scheffer, M.; Carpenter, S. R., Catastrophic regime shifts in ecosystems: linking theory to observation, Trends Ecol. Evol., 18, 648-656, (2003) |

[44] | Shigesada, N., The effects of interference competition on stability, structure and invasion of a multi-species system, J. Math. Biol., 21, 97-113, (1984) · Zbl 0568.92017 |

[45] | Shurin, J. B., Alternative stable states and regional community structure, J. Theoret. Biol., 227, 359-368, (2004) |

[46] | Soetaert, K., Solving differential equations in R: package desolve, J. Stat. Softw., 33, 1-25, (2010) |

[47] | Svenning, J.-C., The influence of interspecific interactions on species range expansion rates, Ecography, 37, 01-012, (2014) |

[48] | Svirezhev, Y., Lotka-Volterra models and the global vegetation pattern, Ecol. Modell., 135, 135-146, (2000) |

[49] | Svirezhev, Y. M., Nonlinearities in mathematical ecology: phenomena and models — would we live in volterra’s world?, Ecol. Modell., 216, 89-101, (2008) |

[50] | Takeuchi, Y.; Adachi, N., The existence of globally stable equilibria of ecosystems of the generalized Volterra type, J. Math. Biol., 10, 401-415, (1980) · Zbl 0458.92019 |

[51] | Van den Driessche, P.; Zeeman, M. L., Three-dimensional competitive Lotka-Volterra systems with no periodic orbits, SIAM J. Appl. Math., 58, 227-234, (1998) · Zbl 0910.34050 |

[52] | van Nes, E. H.; Scheffer, M., Large species shifts triggered by small forces, Am. Natur., 164, 255-266, (2004) |

[53] | Wang, Y.; Wu, H., Dynamics of competitive Lotka-Volterra systems that can be projected to a line, Comput. Math. Appl., 47, 1263-1271, (2004) · Zbl 1154.34354 |

[54] | Wisz, M. S., The role of biotic interactions in shaping distributions and realised assemblages of species: implications for species distribution modelling, Biol. Rev., 88, 15-30, (2013) |

[55] | Wolkowicz, H.; Styan, G. P.H., Bounds for eigenvalues using traces, Linear Algebra Appl., 29, 471-506, (1980) · Zbl 0435.15015 |

[56] | Wolkowicz, H.; Styan, G. P.H., More bounds for eigenvalues using traces, Linear Algebra Appl., 31, 1-17, (1980) · Zbl 0434.15003 |

[57] | Zeeman, M. L., Hopf bifurcations in competitive threedimensional Lotka-Volterra systems, Dyn. Stab. Syst., 8, 189-216, (1993) · Zbl 0797.92025 |

[58] | Zeeman, M. L., Extinction in competitive Lotka-Volterra systems, Proc. Amer. Math. Soc., 123, 87-96, (1995) · Zbl 0815.34039 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.