Lotka-Volterra equations: Decomposition,stability, and structure

Summary The major objective of this paper is to propose a new decomposition-aggregation framework for stability analysis of Lotka-Volterra equations employing the concept of vector Liapunov functions. Both the disjoint and the overlapping decompositions are introduced to increase flexibility in constructing Liapunov functions for the overall system. Our second objective is to consider the Lotka-Volterra equations under structural perturbations, and derive conditions under which a positive equilibrium is connectively stable. Both objectives of this paper are directed towards a better understanding of the intricate interplay between stability and complexity in the context of robustness of model ecosystems represented by Lotka-Volterra equations. Only stability of equilibria in models with constant parameters is considered here. Nonequilibrium analysis of models with nonlinear time-varying parameters is the subject of a companion paper.Research supported by U.S. Department of Energy under the Contract EC-77-S-03-1493.On leave from Kobe University, Kobe, Japan.     

On the generality of stability-complexity relationships in Lotka-Volterra ecosystems

Understanding how complexity persists in nature is a long-standing goal of ecologists. In theoretical ecology, local stability is a widely used measure of ecosystem persistence and has made a major contribution to the ecosystem stability-complexity debate over the last few decades. However, permanence is coming to be regarded as a more satisfactory definition of ecosystem persistence and has relatively recently become available as a tool for assessing the global stability of Lotka-Volterra communities. Here we document positive relationships between permanence and Lotka-Volterra food web complexity and report a positive correlation between the probability of local stability and permanence. We investigate further the frequency of discrepancy (attributed to fragile systems that are locally stable but not permanent or locally unstable systems that are permanent and have cyclic or chaotic dynamics), associate non-permanence with the local stability or instability of equilibria on the boundary of the state-space, and investigate how these vary with aspects of ecosystem complexity. We find that locally stable interior equilibria tend to have all locally unstable boundary equilibria. Since a locally stable boundary is inconsistent with permanent dynamics, this can explain the observed positive correlation between local interior stability and permanence. Our key finding is that, at least in Lotka-Volterra model ecosystems, local stability may be a better measure of persistence than previously thought.     

非自治Lotka-Volterra型捕食扩散系统中有害时滞对该系统持久性的影响

研究了时滞对一类非自治Lotka-Volterra型捕食扩散系统的影响,该系统由n个斑块组成,食饵种群可以在斑块间迁移,而摘食者限制在某一个斑块不能扩散．我们假设密度制约项系数并不总是严格正的．通过运用比较定理及时滞泛函微分方程的基本原理,分两种情况表明了在一定条件下系统是一致持久的．两种情况的结果表明时滞的引入和变化即可能是“有害”,也可能是”无害”．进一步还说明了系统在一致持久性的条件下至少存在一个正周期解．这些结果是对已知的非自治Lotka-Volterra系统的一些结果的推广与改进．     

General equilibrium of an ecosystem

Ecosystems and economies are inextricably linked: ecosystem models and economic models are not linked. Consequently, using either type of model to design policies for preserving ecosystems or improving economic performance omits important information. Improved policies would follow from a model that links the systems and accounts for the mutual feedbacks by recognizing how key ecosystem variables influence key economic variables, and vice versa. Because general equilibrium economic models already are widely used for policy making, the approach used here is to develop a general equilibrium ecosystem model which captures salient biological functions and which can be integrated with extant economic models. In the ecosystem model, each organism is assumed to be a net energy maximizer that must exert energy to capture biomass from other organisms. The exerted energies are the "prices" that are paid to biomass, and each organism takes the prices as signals over which it has no control. The maximization problem yields the organism's demand for and supply of biomass to other organisms as functions of the prices. The demands and supplies for each biomass are aggregated over all organisms in each species which establishes biomass markets wherein biomass prices are determined. A short-run equilibrium is established when all organisms are maximizing and demand equals supply in every biomass market. If a species exhibits positive (negative) net energy in equilibrium, its population increases (decreases) and a new equilibrium follows. The demand and supply forces in the biomass markets drive each species toward zero stored energy and a long-run equilibrium. Population adjustments are not based on typical Lotka-Volterra differential equations in which one entire population adjusts to another entire population thereby masking organism behavior; instead, individual organism behavior is central to population adjustments. Numerical simulations use a marine food web in Alaska to illustrate the model and to show several simultaneous predator/prey relationships, prey switching by the top predator, and energy flows through the web.     

Emergence of a complex and stable network in a model ecosystem with extinction and mutation

We propose a minimal model of the dynamics of diversity-replicator equations with extinction, invasion and mutation. We numerically study the behavior of this simple model and show that it displays completely different behavior from the conventional replicator equation and the generalized Lotka-Volterra equation. We reach several significant conclusions as follows: (1) a complex ecosystem can emerge when mutants with respect to species-specific interaction are introduced; (2) such an ecosystem possesses strong resistance to invasion; (3) a typical fixation process of mutants is realized through the rapid growth of a group of mutualistic mutants with higher fitness than majority species; (4) a hierarchical taxonomic structure (like family-genus-species) emerges; and (5) the relative abundance of species exhibits a typical pattern widely observed in nature. Several implications of these results are discussed in connection with the relationship of the present model to the generalized Lotka-Volterra equation.     

The effects of interference competition on stability,structure and invasion of a multi-species system

Summary An interference competition model for a many species system is presented, based on Lotka-Volterra equations in which some restrictions are imposed on the parameters. The competition coefficients of the Lotka-Volterra equations are assumed to be expressed as products of two factors: the intrinsic interference to other individuals and the defensive ability against such interference. All the equilibrium points of the model are obtained explicitly in terms of its parameters, and these equilibria are classified according to the concept of sector stability. Thus survival or extinction of species at a stable equilibrium point can be determined analytically.The result of the analysis is extended to the successional processes of a community. A criterion for invasion of a new species is obtained and it is also shown that there are some characteristic quantities which show directional changes as succession proceeds.     

Multiple hysteretic patterns from elementary population models

Critical transitions whereby small changes in conditions can cause large and irreversible changes in ecosystem states are a cause of increasing concern in ecology. Here, we focus on the irreversibility of these transitions, formally known as hysteresis. We explore how simple correlations between parameters in Lotka-Volterra predator-prey equations result in a variety of complicated hysteretic patterns. These patterns include “unattainable” stable states that once lost may never be recovered. We suspect these patterns to be common in natural systems, where interactions between diverse assemblages are unavoidable. Thus, understanding underlying hysteretic structures may be necessary for rescuing lost ecosystem states and avoiding future losses.     

Qualitative permanence of Lotka–Volterra equations

In this paper, we consider permanence of Lotka-Volterra equations. We investigate the sign structure of the interaction matrix that guarantees the permanence of a Lotka-Volterra equation whenever it has a positive equilibrium point. An interaction matrix with this property is said to be qualitatively permanent. Our results provide both necessary and sufficient conditions for qualitative permanence.     

When the feasibility of an ecosystem is sufficient for global stability?

We show via a Liapunov function that in every model ecosystem governed by generalized Lotka-Volterra equations, a feasible steady state is globally asymptotically stable if the number of interaction branches equals n-1, where n is the number of species. This means that the representative graph for which the theorem holds is a 'tree' and not only an alimentary chain. Our result is valid also in the case of non-homogeneous systems, which model situations in which input fluxes are present.     

Invariant distributions for multi-population models in random environments

We obtain the existence of a solution and invariant distribution for systems of stochastic differential equations which represent populations in random environments. The method used is a stochastic Lyapunov function, based on a theorem of Kushner. The method is applied to a system of two populations exchainging individuals through migration, and to a generalized n-dimensional Lotka-Volterra system.     

Successful invasion of a food web in a chemostat

A food web in a chemostat is considered in which an arbitrary number of competitor populations compete for a single, essential, nonreproducing, growth-limiting substrate, and an arbitrary number of predator populations prey on some or all of the competitor populations. Although any number of predator populations may prey on the same competitor population, each predator population preys on only one competitor population. The dynamics of substrate uptake is modeled by Lotka-Volterra or Michaelis-Menten (Holling type I or II), but the dynamics of competitor uptake is restricted to Lotka-Volterra. Based on certain parameters, the model predicts the asymptotic survival or extinction of each of the different populations and suggests how competitor and/or predator populations could successfully invade the chemostat with or without causing a diverse ecosystem to crash. Similarly, it suggests how the elimination of certain populations could result in a more diverse or less diverse system.     

Effective competition determines the global stability of model ecosystems

We investigate the stability of Lotka-Volterra (LV) models constituted by two groups of species such as plants and animals in terms of the intragroup effective competition matrix, which allows separating the equilibrium equations of the two groups. In matrix analysis, the effective competition matrix represents the Schur complement of the species interaction matrix. It has been previously shown that the main eigenvalue of this effective competition matrix strongly influences the structural stability of the model ecosystem. Here, we show that the spectral properties of the effective competition matrix also strongly influence the dynamical stability of the model ecosystem. In particular, a necessary condition for diagonal stability of the full system, which guarantees global stability, is that the effective competition matrix is diagonally stable, which means that intergroup interactions must be weaker than intra-group competition in appropriate units. For mutualistic or competitive interactions, diagonal stability of the effective competition is a sufficient condition for global stability if the inter-group interactions are suitably correlated, in the sense that the biomass that each species provides to (removes from) the other group must be proportional to the biomass that it receives from (is removed by) it. For a non-LV mutualistic system with saturating interactions, we show that the diagonal stability of the corresponding LV system close to the fixed point is a sufficient condition for global stability.     

Almost periodic solution of non-autonomous Lotka-Volterra predator-prey dispersal system with delays

This paper studies a non-autonomous Lotka-Volterra almost periodic predator-prey dispersal system with discrete and continuous time delays which consists of n-patches, the prey species can disperse among n-patches, but the predator species is confined to one patch and cannot disperse. By using comparison theorem and delay differential equation basic theory, we prove the system is uniformly persistent under some appropriate conditions. Further, by constructing suitable Lyapunov functional, we show that the system is globally asymptotically stable under some appropriate conditions. By using almost periodic functional hull theory, we show that the almost periodic system has a unique globally asymptotical stable strictly positive almost periodic solution. The conditions for the permanence, global stability of system and the existence, uniqueness of positive almost periodic solution depend on delays, so, time delays are "profitless". Finally, conclusions and two particular cases are given. These results are basically an extension of the known results for non-autonomous Lotka-Volterra systems.     

种间竞争的一个新的数学模型—对经典的Lotka—Vol—terra竞争方程的扩充

本文根据营养动力学理论,建立了一类种间竞争的新的数学模型:它是单种群增长的Cui-Lawson模型,在种间竞争上的推广。新的种间竞争模型克服了经典的种间竞争的Lotka-Volteira方程的局限与不足,具有更广泛和复杂的行为,并在特殊条件下以Lotka-Volterra竞争方程为其特例。因此,新的种间竞争的数学模型是更一般的解释性模型,是对经典的Lotka-Voterra竞争方程的扩充。     

Context dependency of relationships between biodiversity and ecosystem functioning is different for multiple ecosystem functions

Increasing concern over the loss of biodiversity has led to attempts to quantify relationships between biodiversity and ecosystem functioning. While manipulative investigations have accumulated substantial evidence to support the notion that decreasing biodiversity can be detrimental to the functioning of ecosystems, recent investigations have identified the potential importance of physical processes in moderating biodiversity – ecosystem function relationships at larger geographical scales. In this study, the relationship between the genus richness of benthic macro‐invertebrates and five measures of ecosystem functioning (macrofaunal biomass, depth of the apparent redox discontinuity, fluxes of ammonium and NO_x and the abundance of nematodes) was determined over a large scale wave‐induced bed shear stress gradient on the seabed of the northern Irish Sea. Ecosystem functioning was significantly correlated to genus richness for four out of five ecosystem functions. However, wave stress moderated the genus richness – ecosystem functioning relationship for only one of the ecosystem functions; genus richness had a positive effect on the depth of the apparent redox discontinuity in the sediment at high wave stress but not at low wave stress. These results indicate that the effects of biodiversity on some ecosystem functions may be sufficiently strong to generate patterns in ecosystems where other factors are also affecting ecosystem processes, but that the biodiversity–ecosystem function relationship for can be dependent on environmental conditions for specific ecosystem functions.     

Modeling arbuscular mycorrhizal infection: is % infection an appropriate variable?

Several available models of arbuscular mycorrhizal infection are based on fitting % infection to a logistic curve and then relating the various parameters to biological functions. I suggest here that this direction is misleading. Percent infection is a value derived from the growth of two interdependent but distinct organisms, each of which is seeking to maximize its own growth and survival. I suggest that two-organism models, such as those derived from Lotka-Volterra equations, are more useful for understanding the biology and functioning of mycorrhizae. Accepted: 22 October 2000     

一类具分段常数变量的捕食-食饵系统的Neimark-Sacker分支

讨论了具时滞与分段常数变量的捕食-食饵生态模型的稳定性及Neimark-Sacker分支;通过计算得到连续模型对应的差分模型,基于特征值理论和Schur-Cohn判据得到正平衡态局部渐进稳定的充分条件;以食饵的内禀增长率为分支参数,运用分支理论和中心流形定理分析了Neimark-Sacker分支的存在性与稳定性条件;通过举例和数值模拟验证了理论的正确性。     

Individual based modeling and parameter estimation for a Lotka-Volterra system

Stochastic component, inevitable in biological systems, makes problematic the estimation of the model parameters from a single sequence of measurements, despite the complete knowledge of the system. We studied the problem of parameter estimation using individual-based computer simulations of a 'Lotka-Volterra world'. Two kinds (species) of particles--X (preys) and Y (predators)--moved on a sphere according to deterministic rules and at the collision (interaction) of X and Y the particle X was changed to a new particle Y. Birth of preys and death of predators were simulated by addition of X and removal of Y, respectively, according to exponential probability distributions. With this arrangement of the system, the numbers of particles of each kind might be described by the Lotka-Volterra equations. The simulations of the system with low (200-400 particles on average) number of individuals showed unstable oscillations of the population size. In some simulation runs one of the species became extinct. Nevertheless, the oscillations had some generic properties (e.g. mean, in one simulation run, oscillation period, mean ratio of the amplitudes of the consecutive maxima of X and Y numbers, etc.) characteristic for the solutions of the Lotka-Volterra equations. This observation made it possible to estimate the four parameters of the Lotka-Volterra model with high accuracy and good precision. The estimation was performed using the integral form of the Lotka-Volterra equations and two parameter linear regression for each oscillation cycle separately. We conclude that in spite of the irregular time course of the number of individuals in each population due to stochastic intraspecies component, the generic features of the simulated system evolution can provide enough information for quantitative estimation of the system parameters.