A metapopulation consists of a group of spatially separated populations of the same species which interact at some level. The term metapopulation was coined by Richard Levins in 1969 to describe a model of population dynamics of insect pests in agricultural fields, but the idea has been most broadly applied to species in naturally or artificially fragmented habitats. In Levins' own words, it consists of "a population of populations".
A metapopulation is generally considered to consist of several distinct populations together with areas of suitable habitat which are currently unoccupied. In classical metapopulation theory, each population cycles in relative independence of the other populations and eventually goes extinct as a consequence of demographic stochasticity (fluctuations in population size due to random demographic events); the smaller the population, the more chances of inbreeding depression and prone to extinction.
Although individual populations have finite life-spans, the metapopulation as a whole is often stable because immigrants from one population (which may, for example, be experiencing a population boom) are likely to re-colonize habitat which has been left open by the extinction of another population. They may also emigrate to a small population and rescue that population from extinction (called the rescue effect). Such a rescue effect may occur because declining populations leave niche opportunities open to the "rescuers".
The development of metapopulation theory, in conjunction with the development of source-sink dynamics, emphasised the importance of connectivity between seemingly isolated populations. Although no single population may be able to guarantee the long-term survival of a given species, the combined effect of many populations may be able to do this.
Metapopulation theory was first developed for terrestrial ecosystems, and subsequently applied to the marine realm. In fisheries science, the term "sub-population" is equivalent to the metapopulation science term "local population". Most marine examples are provided by relatively sedentary species occupying discrete patches of habitat, with both local recruitment and recruitment from other local populations in the larger metapopulation. Kritzer & Sale have argued against strict application of the metapopulation definitional criteria that extinction risks to local populations must be non-negligible.:32
The first experiments with predation and spatial heterogeneity were conducted by G.F. Gause in the 1930s, based on the Lotka-Volterra equation, which was formulated in the mid-1920s, but no further application had been conducted. The Lotka-Volterra equation suggested that the relationship between predators and their prey would result in population oscillations over time based on the initial densities of predator and prey. Gause's early experiments to prove the predicted oscillations of this theory failed because the predator–prey interactions were not influenced by immigration. However, once immigration was introduced, the population cycles accurately depicted the oscillations predicted by the Lotka-Volterra equation, with the peaks in prey abundance shifted slightly to the left of the peaks of the predator densities. Huffaker's experiments expanded on those of Gause by examining how both the factors of migration and spatial heterogeneity lead to predator–prey oscillations.
In order to study predation and population oscillations, Huffaker used mite species, one being the predator and the other being the prey. He set up a controlled experiment using oranges, which the prey fed on, as the spatially structured habitat in which the predator and prey would interact. At first, Huffaker experienced difficulties similar to those of Gause in creating a stable predator–prey interaction. By using oranges only, the prey species quickly became extinct followed consequently with predator extinction. However, he discovered that by modifying the spatial structure of the habitat, he could manipulate the population dynamics and allow the overall survival rate for both species to increase. He did this by altering the distance between the prey and oranges (their food), establishing barriers to predator movement, and creating corridors for the prey to disperse. These changes resulted in increased habitat patches and in turn provided more areas for the prey to seek temporary protection. When the prey would become extinct locally at one habitat patch, they were able to reestablish by migrating to new patches before being attacked by predators. This habitat spatial structure of patches allowed for coexistence between the predator and prey species and promoted a stable population oscillation model. Although the term metapopulation had not yet been coined, the environmental factors of spatial heterogeneity and habitat patchiness would later describe the conditions of a metapopulation relating to how groups of spatially separated populations of species interact with one another. Huffaker's experiment is significant because it showed how metapopulations can directly affect the predator–prey interactions and in turn influence population dynamics.
Levins' original model applied to a metapopulation distributed over many patches of suitable habitat with significantly less interaction between patches than within a patch. Population dynamics within a patch were simplified to the point where only presence and absence were considered. Each patch in his model is either populated or not.
Let N be the fraction of patches occupied at a given time. During a time dt, each occupied patch can become unoccupied with an extinction probability edt. Additionally, 1 − N of the patches are unoccupied. Assuming a constant rate c of propagule generation from each of the N occupied patches, during a time dt, each unoccupied patch can become occupied with a colonization probability cNdt . Accordingly, the time rate of change of occupied patches, dN/dt, is
This equation is mathematically equivalent to the logistic model, with a carrying capacity K given by
and growth rate r
At equilibrium, therefore, some fraction of the species's habitat will always be unoccupied.
Huffaker's studies of spatial structure and species interactions are an example of early experimentation in metapopulation dynamics. Since the experiments of Huffaker and Levins, models have been created which integrate stochastic factors. These models have shown that the combination of environmental variability (stochasticity) and relatively small migration rates cause indefinite or unpredictable persistence. However, Huffaker's experiment almost guaranteed infinite persistence because of the controlled immigration variable.
One major drawback of the Levins model is that it is deterministic, whereas the fundamental metapopulation processes are stochastic. Metapopulations are particularly useful when discussing species in disturbed habitats, and the viability of their populations, i.e., how likely they are to become extinct in a given time interval. The Levins model cannot address this issue. A simple way to extend the Levins' model to incorporate space and stochastic considerations is by using the contact process. Simple modifications to this model can also incorporate for patch dynamics. At a given percolation threshold, habitat fragmentation effects take place in these configurations predicting more drastic extinction thresholds.
For conservation biology purposes, metapopulation models must include (a) the finite nature of metapopulations (how many patches are suitable for habitat), and (b) the probabilistic nature of extinction and colonisation. Also, note that in order to apply these models, the extinctions and colonisations of the patches must be asynchronous.
By combining nanotechnology with landscape ecology, a habitat landscape can be nanofabricated on-chip by building a collection of nanofabricated bacterial habitats, and connecting them by corridors in different topological arrangements and with nano-scale channels providing them with the local ecosystem service of habitat renewal. These landscapes of MHPs can be used as physical implementations of an adaptive landscape: by generating a spatial mosaic of patches of opportunity distributed in space and time. The patchy nature of these fluidic landscapes allows for the study of adapting bacterial cells in a metapopulation system operating on-chip within a synthetic ecosystem. The metapopulation biology and evolutionary ecology of these bacterial systems, in these synthetic ecosystems, can be addressed using experimental biophysics.
Metapopulation models have been used to explain life-history evolution, such as the ecological stability of amphibian metamorphosis in small vernal ponds. Alternative ecological strategies have evolved. For example, some salamanders forgo metamorphosis and sexually mature as aquatic neotenes. The seasonal duration of wetlands and the migratory range of the species determines which ponds are connected and if they form a metapopulation. The duration of the life history stages of amphibians relative to the duration of the vernal pool before it dries up regulates the ecological development of metapopulations connecting aquatic patches to terrestrial patches.