# Theoretical ecology

### Theoretical ecology

Mathematical models developed in theoretical ecology predict complex food webs are less stable than simple webs.[1]:75–77[2]:64

Theoretical ecology is the scientific discipline devoted to the study of ecological systems using theoretical methods such as simple conceptual models, mathematical models, computational simulations, and advanced data analysis. Effective models improve understanding of the natural world by revealing how the dynamics of species populations are often based on fundamental biological conditions and processes. Further, the field aims to unify a diverse range of empirical observations by assuming that common, mechanistic processes generate observable phenomena across species and ecological environments. Based on biologically realistic assumptions, theoretical ecologists are able to uncover novel, non-intuitive insights about natural processes. Theoretical results are often verified by empirical observation, revealing the power of theoretical methods in both predicting and understanding the noisy, diverse biological world.

The field is broad and includes foundations in applied mathematics, computer science, biology, statistical physics, genetics, chemistry, evolution, and conservation biology. Theoretical ecology aims to explain a diverse range of phenomena in the life sciences, such as population growth and dynamics, fisheries, competition, evolutionary theory, epidemiology, animal behavior and group dynamics, food webs, ecosystems, spatial ecology, and the effects of climate change.

Theoretical ecology has further benefited from the advent of fast computing power, allowing the analysis and visualization of large-scale computational simulations of ecological phenomena. Importantly, these modern tools provide quantitative predictions about the effects of human induced environmental change on a diverse variety of ecological phenomena, such as: species invasions, climate change, the effect of fishing and hunting on food network stability, and the global carbon cycle.

## Contents

• Modelling approaches 1
• Population ecology 2
• Exponential growth 2.1
• Logistic growth 2.2
• Structured population growth 2.3
• Community ecology 3
• Predator-prey 3.1
• Host-pathogen 3.2
• Host-parasitoid 3.3
• Competition and mutualism 3.4
• Spatial ecology 4
• Biogeography 4.1
• Niche theory 4.2
• Neutral theory 4.3
• Metapopulations 4.4
• Ecosystem ecology 5
• Food webs 5.1
• Systems ecology 5.2
• Ecophysiology 6
• Behavioral ecology 7
• Swarm behaviour 7.1
• Evolutionary ecology 8
• Other theories 9
• History 10
• Theoretical and mathematical ecologists 11
• Journals 12
• References 14

## Modelling approaches

As in most other sciences, mathematical models form the foundation of modern ecological theory.

• Phenomenological models: distill the functional and distributional shapes from observed patterns in the data, or researchers decide on functions and distribution that are flexible enough to match the patterns they or others (field or experimental ecologists) have found in the field or through experimentation.[3]
• Mechanistic models: model the underlying processes directly, with functions and distributions that are based on theoretical reasoning about ecological processes of interest.[3]

Ecological models can be deterministic or stochastic.[3]

• Deterministic models always evolve in the same way from a given starting point.[4] They represent the average, expected behavior of a system, but lack random variation. Many system dynamics models are deterministic.
• Stochastic models allow for the direct modeling of the random perturbations that underlie real world ecological systems. Markov chain models are stochastic.

Species can be modelled in continuous or discrete time.[5]

• Continuous time is modelled using differential equations.
• Discrete time is modelled using difference equations. These model ecological processes that can be described as occurring over discrete time steps. Matrix algebra is often used to investigate the evolution of age-structured or stage-structured populations. The Leslie matrix, for example, mathematically represents the discrete time change of an age structured population.[6][7][8]

Models are often used to describe real ecological reproduction processes of single or multiple species. These can be modelled using stochastic branching processes. Examples are the dynamics of interacting populations (predation competition and mutualism), which, depending on the species of interest, may best be modeled over either continuous or discrete time. Other examples of such models may be found in the field of mathematical epidemiology where the dynamic relationships that are to be modeled are host-pathogen interactions.[5]

Bifurcation diagram of the logistic map

Bifurcation theory is used to illustrate how small changes in parameter values can give rise to dramatically different long run outcomes, a mathematical fact that may be used to explain drastic ecological differences that come about in qualitatively very similar systems.[9] Logistic maps are polynomial mappings, and are often cited as providing archetypal examples of how chaotic behaviour can arise from very simple non-linear dynamical equations. The maps were popularized in a seminal 1976 paper by the theoretical ecologist Robert May.[10] The difference equation is intended to capture the two effects of reproduction and starvation.

In 1930,

A distinction can be made between mathematical ecologists, ecologists who apply mathematics to ecological problems, and mathematicians who develop the mathematics itself that arises out of ecological problems.

## Theoretical and mathematical ecologists

Stephen Hubbell and Michael Rosenzweig combined theoretical and practical elements into works that extended MacArthur and Wilson's Island Biogeography Theory - Hubbell with his Unified Neutral Theory of Biodiversity and Biogeography and Rosenzweig with his Species Diversity in Space and Time.

Simberloff added statistical rigour to experimental ecology and was a key figure in the SLOSS debate, about whether it is preferable to protect a single large or several small reserves.[63] This resulted in the supporters of Jared Diamond's community assembly rules defending their ideas through Neutral Model Analysis.[63] Simberloff also played a key role in the (still ongoing) debate on the utility of corridors for connecting isolated reserves.

Theoretical ecology draws on pioneering work done by G. Evelyn Hutchinson and his students. Brothers H.T. Odum and E.P. Odum are generally recognised as the founders of modern theoretical ecology. Robert MacArthur brought theory to community ecology. Daniel Simberloff was the student of E.O. Wilson, with whom MacArthur collaborated on The Theory of Island Biogeography, a seminal work in the development of theoretical ecology.[62]

## History

In contrast to previous ecological theories which considered floods to be catastrophic events, the river flood pulse concept argues that the annual flood pulse is the most important aspect and the most biologically productive feature of a river's ecosystem.[60][61]

## Other theories

The cybernetician and anthropologist Gregory Bateson observed in the 1970s that, though writing it only as an example, Wallace had "probably said the most powerful thing that’d been said in the 19th Century".[59] Subsequently, the connection between natural selection and systems theory has become an area of active research.[57]

The action of this principle is exactly like that of the centrifugal governor of the steam engine, which checks and corrects any irregularities almost before they become evident; and in like manner no unbalanced deficiency in the animal kingdom can ever reach any conspicuous magnitude, because it would make itself felt at the very first step, by rendering existence difficult and extinction almost sure soon to follow.[58]

The British biologist Alfred Russel Wallace is best known for independently proposing a theory of evolution due to natural selection that prompted Charles Darwin to publish his own theory. In his famous 1858 paper, Wallace proposed natural selection as a kind of feedback mechanism which keeps species and varieties adapted to their environment.[57]

## Evolutionary ecology

Recently, a number of mathematical models have been discovered which explain many aspects of the emergent behaviour. Swarm algorithms follow a Lagrangian approach or an Eulerian approach.[52] The Eulerian approach views the swarm as a field, working with the density of the swarm and deriving mean field properties. It is a hydrodynamic approach, and can be useful for modelling the overall dynamics of large swarms.[53][54][55] However, most models work with the Lagrangian approach, which is an agent-based model following the individual agents (points or particles) that make up the swarm. Individual particle models can follow information on heading and spacing that is lost in the Eulerian approach.[52][56] Examples include ant colony optimization, self-propelled particles and particle swarm optimization

emergent behaviour that occurs when individual agents follow simple behavioral rules.

Flocks of birds can abruptly change their direction in unison, and then, just as suddenly, make a unanimous group decision to land.[51]

## Ecophysiology

Systems ecology can be seen as an application of general systems theory to ecology. It takes a holistic and interdisciplinary approach to the study of ecological systems, and particularly ecosystems. Systems ecology is especially concerned with the way the functioning of ecosystems can be influenced by human interventions. Like other fields in theoretical ecology, it uses and extends concepts from thermodynamics and develops other macroscopic descriptions of complex systems. It also takes account of the energy flows through the different trophic levels in the ecological networks. In systems ecology the principles of ecosystem energy flows are considered formally analogous to the principles of energetics. Systems ecology also considers the external influence of ecological economics, which usually is not otherwise considered in ecosystem ecology.[50] For the most part, systems ecology is a subfield of ecosystem ecology.

### Systems ecology

In 1927, Charles Elton published an influential synthesis on the use of food webs, which resulted in them becoming a central concept in ecology.[48] In 1966, interest in food webs increased after Robert Paine's experimental and descriptive study of intertidal shores, suggesting that food web complexity was key to maintaining species diversity and ecological stability.[49] Many theoretical ecologists, including Sir Robert May and Stuart Pimm, were prompted by this discovery and others to examine the mathematical properties of food webs. According to their analyses, complex food webs should be less stable than simple food webs.[1]:75–77[2]:64 The apparent paradox between the complexity of food webs observed in nature and the mathematical fragility of food web models is currently an area of intensive study and debate. The paradox may be due partially to conceptual differences between persistence of a food web and equilibrial stability of a food web.[1][2]

[47]

### Food webs

Ecopath is a free ecosystem modelling software suite, initially developed by NOAA, and widely used in fisheries management as a tool for modelling and visualising the complex relationships that exist in real world marine ecosystems.

If ecosystems are governed primarily by stochastic processes, through which its subsequent state would be determined by both predictable and random actions, they may be more resilient to sudden change than each species individually. In the absence of a balance of nature, the species composition of ecosystems would undergo shifts that would depend on the nature of the change, but entire ecological collapse would probably be infrequent events. In 1997, Robert Ulanowicz used information theory tools to describe the structure of ecosystems, emphasizing mutual information (correlations) in studied systems. Drawing on this methodology and prior observations of complex ecosystems, Ulanowicz depicts approaches to determining the stress levels on ecosystems and predicting system reactions to defined types of alteration in their settings (such as increased or reduced energy flow, and eutrophication.[45]

Introducing new elements, whether biotic or abiotic, into ecosystems can be disruptive. In some cases, it leads to ecological collapse, trophic cascades and the death of many species within the ecosystem. The abstract notion of ecological health attempts to measure the robustness and recovery capacity for an ecosystem; i.e. how far the ecosystem is away from its steady state. Often, however, ecosystems rebound from a disruptive agent. The difference between collapse or rebound depends on the toxicity of the introduced element and the resiliency of the original ecosystem.

## Ecosystem ecology

The model can also be extended to combinations of the four possible linear or non-linear dependencies of colonization and extinction on p are described in more detail in.[44]

1. Colonization may be dependent on p linearly (m*(1-p)) as opposed to the non-linear m*p*(1-p) regime described above. This mode of replication of a species is called the “rain of propagules”, where there is an abundance of new individuals entering the population at every generation. In such a scenario, the steady state where the population is zero is usually unstable.[43]
2. Extinction may depend non-linearly on p (e*p*(1-p)) as opposed to the linear (e*p) regime described above. This is referred to as the “rescue effect” and it is again harder to drive a population extinct under this regime.[43]

The form of the differential equations used in this simplistic modelling approach can be modified. For example:

(p*1 may be inferred by symmetry). It is interesting to note that if e is zero, the dynamics of the system favor the species that is better at colonizing (i.e. has the higher m value). This leads to a very important result in theoretical ecology known as the Intermediate Disturbance Hypothesis, where the biodiversity (the number of species that coexist in the population) is maximized when the disturbance (of which e is a proxy here) is not too high or too low, but at intermediate levels.[42]

p^*_2 = \frac{e}{m_1} - \frac{m_1}{m_2}

In this case, if e is too high, p1 and p2 will be zero at steady state. However, when the rate of extinction is moderate, p1 and p2 can stably coexist. The steady state value of p2 is given by

\frac{dp_1}{dt} = m_1 p_1 (1 - p_1) - e p_1
\frac{dp_2}{dt} = m_2 p_2 (1 - p_1 - p_2) - e p_2 - m p_1 p_2

where m is the rate of colonization, and e is the rate of extinction.[41] In this model, if e < m, the steady state value of p is 1 – (e/m) while in the other case, all the patches will eventually be left empty. This model may be made more complex by addition of another species in several different ways, including but not limited to game theoretic approaches, predator-prey interactions, etc. We will consider here an extension of the previous one-species system for simplicity. Let us denote the proportion of patches occupied by the first population as p1, and that by the second as p2. Then,

\frac{dp}{dt} = m p (1-p) - e p

Spatial analysis of ecological systems often reveals that assumptions that are valid for spatially homogenous populations – and indeed, intuitive – may no longer be valid when migratory subpopulations moving from one patch to another are considered.[40] In a simple one-species formulation, a subpopulation may occupy a patch, move from one patch to another empty patch, or die out leaving an empty patch behind. In such a case, the proportion of occupied patches may be represented as

### Metapopulations

Hubbell built on earlier neutral concepts, including MacArthur & Wilson's theory of island biogeography[20] and Gould's concepts of symmetry and null models.[38]

Under unified neutral theory, complex ecological interactions are permitted among individuals of an ecological community (such as competition and cooperation), providing all individuals obey the same rules. Asymmetric phenomena such as parasitism and predation are ruled out by the terms of reference; but cooperative strategies such as swarming, and negative interaction such as competing for limited food or light are allowed, so long as all individuals behave the same way. The theory makes predictions that have implications for the management of biodiversity, especially the management of rare species. It predicts the existence of a fundamental biodiversity constant, conventionally written θ, that appears to govern species richness on a wide variety of spatial and temporal scales.

Unified neutral theory is a hypothesis proposed by Stephen Hubbell in 2001.[20] The hypothesis aims to explain the diversity and relative abundance of species in ecological communities, although like other neutral theories in ecology, Hubbell's hypothesis assumes that the differences between members of an ecological community of trophically similar species are "neutral," or irrelevant to their success. Neutrality means that at a given trophic level in a food web, species are equivalent in birth rates, death rates, dispersal rates and speciation rates, when measured on a per-capita basis.[38] This implies that biodiversity arises at random, as each species follows a random walk.[39] This can be considered a null hypothesis to niche theory. The hypothesis has sparked controversy, and some authors consider it a more complex version of other null models that fit the data better.

### Niche theory

The diversity and containment of coral reef systems make them good sites for testing niche and neutral theories.[37]

Biogeography is most keenly observed on islands, which has led to the development of the subdiscipline of island biogeography. These habitats are often a more manageable areas of study because they are more condensed than larger ecosystems on the mainland. In 1967, Robert MacArthur and E.O. Wilson published The Theory of Island Biogeography. This showed that the species richness in an area could be predicted in terms of factors such as habitat area, immigration rate and extinction rate.[34] The theory is considered one of the fundamentals of ecological theory.[35] The application of island biogeography theory to habitat fragments spurred the development of the fields of conservation biology and landscape ecology.[36]

Biogeography is the study of the distribution of species in space and time. It aims to reveal where organisms live, at what abundance, and why they are (or are not) found in a certain geographical area.

## Spatial ecology

See Competitive Lotka-Volterra equations for further extensions of this model.

The r coefficients give a “base” growth rate to each species, while K coefficients correspond to the carrying capacity. What can really change the dynamics of a system, however are the α terms. These describe the nature of the relationship between the two species. When α12 is negative, it means that N2 has a negative effect on N1, by competing with it, preying on it, or any number of other possibilities. When α12 is positive, however, it means that N2 has a positive effect on N1, through some kind of mutualistic interaction between the two. When both α12 and α21 are negative, the relationship is described as competitive. In this case, each species detracts from the other, potentially over competition for scarce resources. When both α12 and α21 are positive, the relationship becomes one of mutualism. In this case, each species provides a benefit to the other, such that the presence of one aids the population growth of the other.

\frac{dN_1}{dt} = \frac{r_1N_1}{K_1}\left( K_1 - N_1 + \alpha_{12}N_2 \right)
\frac{dN_2}{dt} = \frac{r_2N_2}{K_2}\left( K_2 - N_2 + \alpha_{21}N_1 \right)

In studies of the populations of two species, the Lotka-Volterra system of equations has been extensively used to describe dynamics of behavior between two species, N1 and N2. Examples include relations between D. discoiderum and E. coli,[32] as well as theoretical analysis of the behavior of the system.[33]

### Competition and mutualism

where f(Nt, Pt) describes the probability of infection (typically, Poisson distribution), λ is the per-capita growth rate of hosts in the absence of parasitoids, and c is the conversion efficiency, as in the Lotka-Volterra model.[21]

N_{t+1} = \lambda \ N_t \ [ 1 - f(N_t, P_t) ]
P_{t+1} = c \ N_t \ f(N_t, p_t)

The third interaction, that of host and parasitoid, can be analyzed by the Nicholson-Bailey model, which differs from Lotka-Volterra and SIR models in that it is discrete in time. This model, like that of Lotka-Volterra, tracks both populations explicitly. Typically, in its general form, it states:

### Host-parasitoid

The second interaction, that of host and pathogen, differs from predator-prey interactions in that pathogens are much smaller, have much faster generation times, and require a host to reproduce. Therefore, only the host population is tracked in host-pathogen models. Compartmental models that categorize host population into groups such as susceptible, infected, and recovered (SIR) are commonly used.[31]

### Host-pathogen

A credible, simple alternative to the Lotka-Volterra predator-prey model and their common prey dependent generalizations is the ratio dependent or Arditi-Ginzburg model.[29] The two are the extremes of the spectrum of predator interference models. According to the authors of the alternative view, the data show that true interactions in nature are so far from the Lotka-Volterra extreme on the interference spectrum that the model can simply be discounted as wrong. They are much closer to the ratio dependent extreme, so if a simple model is needed one can use the Arditi-Ginzburg model as the first approximation.[30]

Volterra originally used the model to explain fluctuations in fish and shark populations after fishing was curtailed during the First World War. However, the equations have subsequently been applied more generally.[25] Other examples of these models include the Lotka-Volterra model of the snowshoe hare and Canadian lynx in North America,[26] any infectious disease modeling such as the recent outbreak of SARS [27] and biological control of California red scale by the introduction of its parasitoid, Aphytis melinus .[28]

where N is the prey and P is the predator population sizes, r is the rate for prey growth, taken to be exponential in the absence of any predators, α is the prey mortality rate for per-capita predation (also called ‘attack rate’), c is the efficiency of conversion from prey to predator, and d is the exponential death rate for predators in the absence of any prey.

\frac{dN(t)}{dt} = N(t)(r-\alpha P(t))
\frac{dP(t)}{dt} = P(t)(c \alpha N(t) -d)

Predator-prey interactions exhibit natural oscillations in the populations of both predator and the prey.[21] In 1925, the American mathematician Alfred J. Lotka developed simple equations for predator-prey interactions in his book on biomathematics.[22] The following year, the Italian mathematician Vito Volterra, made a statistical analysis of fish catches in the Adriatic[23] and independently developed the same equations.[24] It is one of the earliest and most recognised ecological models, known as the Lotka-Volterra model:

### Predator-prey

Lotka-Volterra model of cheetah-baboon interactions. Starting with 80 baboons (green) and 40 cheetahs, this graph shows how the model predicts the two species numbers will progress over time.

An ecological community is a group of trophically similar, sympatric species that actually or potentially compete in a local area for the same or similar resources.[20] Interactions between these species form the first steps in analyzing more complex dynamics of ecosystems. These interactions shape the distribution and dynamics of species. Of these interactions, predation is one of the most widespread population activities.[21] Taken in its most general sense, predation comprises predator-prey, host-pathogen, and host parasitoid interactions.

## Community ecology

If parameter values in L are estimated from demographic data on a specific population, a structured model can then be used to predict whether this population is expected to grow or decline in the long-term, and what the expected age distribution within the population will be. This has been done for a number of species including loggerhead sea turtles and right whales.[18][19]

where Nt is a vector of the number of individuals in each class at time t and L is a matrix that contains the survival probability and fecundity for each class. The matrix L is referred to as the Leslie matrix for age-structured models, and as the Lefkovitch matrix for stage-structured models.[17]

\mathbf{N}_{t+1} = \mathbf{L}\mathbf{N}_t

Another assumption of the exponential growth model is that all individuals within a population are identical and have the same probabilities of surviving and of reproducing. This is not a valid assumption for species with complex life histories. The exponential growth model can be modified to account for this, by tracking the number of individuals in different age classes (e.g. one-, two-, and three-year-olds) or different stage classes (juveniles, sub-adults, and adults) separately, and allowing individuals in each group to have their own survival and reproduction rates. The general form of this model is

### Structured population growth

The biological significance of K becomes apparent when stabilities of the equilibria of the system are considered. It is the carrying capacity of the population. The equilibria of the system are N = 0 and N = K. If the system is linearized, it can be seen that N = 0 is an unstable equilibrium while K is a stable equilibrium.[16]

where r = b-d and K = (b-d)/(a+c).

\frac{dN(t)}{dt}=rN(t) \left(1-\frac{N}{K}\right)

This can be rewritten as:[16]

\frac{dN(t)}{dt}=((b-aN(t))-(d-cN(t)))N(t)

The exponential growth model makes a number of assumptions, many of which often do not hold. For example, many factors affect the intrinsic growth rate and is often not time-invariant. A simple modification of the exponential growth is to assume that the intrinsic growth rate varies with population size. This is reasonable: the larger the population size, the fewer resources available, which can result in a lower birth rate and higher death rate. Hence, we can replace the time-invariant r with r’(t) = (b –a*N(t)) – (d + c*N(t)), where a and c are constants that modulate birth and death rates in a population dependent manner (e.g. intraspecific competition). Both a and c will depend on other environmental factors which, we can for now, assume to be constant in this approximated model. The differential equation is now:[16]

### Logistic growth

a trajectory known as Malthusian growth, after Thomas Malthus, who first described its dynamics in 1798. A population experiencing Malthusian growth follows an exponential curve, where N(0) is the initial population size. The population grows when r > 0, and declines when r < 0. The model is most applicable in cases where a few organisms have begun a colony and are rapidly growing without any limitations or restrictions impeding their growth (e.g. bacteria inoculated in rich media).

N(t) = N(0) \ e^{rt} ,

where r is the per capita growth rate, or the intrinsic growth rate of the organism. It can also be described as r = b-d, where b and d are the per capita time-invariant birth and death rates, respectively. This first order linear differential equation can be solved to yield the solution

\frac{dN(t)}{dt}=rN(t)

The most basic way of modeling population dynamics is to assume that the rate of growth of a population depends only upon the population size at that time and the per capita growth rate of the organism. In other words, if the number of individuals in a population at a time t, is N(t), then the rate of population growth is given by:

### Exponential growth

Population ecology is a sub-field of ecology that deals with the dynamics of species populations and how these populations interact with the environment.[15] It is the study of how the population sizes of species living together in groups change over time and space, and was one of the first aspects of ecology to be studied and modelled mathematically.

## Population ecology

Because ecological systems are typically fisheries.[14]

. evolutionarily stable strategy defined the concept of the [13], who in his seminal 1972 paper, “Game Theory and the Evolution of Fighting",John Maynard Smith followed closely by [12],Evolution and the Theory of Games applied game theory to evolutionary biology in his Richard Lewontin In 1961, [11]