Population dynamics of fisheries

Predator bluefin trevally sizing up schooling anchovies, in the Maldives

A fishery is an area with an associated fish or aquatic population which is harvested for its commercial or recreational value. Fisheries can be wild or farmed. Population dynamics describes the ways in which a given population grows and shrinks over time, as controlled by birth, death, and migration. It is the basis for understanding changing fishery patterns and issues such as habitat destruction, predation and optimal harvesting rates. The population dynamics of fisheries is used by fisheries scientists to determine sustainable yields.[1][2][3]

The basic accounting relation for population dynamics is the BIDE (Birth, Immigration, Death, Emigration) model, shown as:[4]

N1 = N0 + B D + I E

where N1 is the number of individuals at time 1, N0 is the number of individuals at time 0, B is the number of individuals born, D the number that died, I the number that immigrated, and E the number that emigrated between time 0 and time 1. While immigration and emigration can be present in wild fisheries, they are usually not measured.

A fishery population is affected by three dynamic rate functions:

If these rates are measured over different time intervals, the harvestable surplus of a fishery can be determined. The harvestable surplus is the number of individuals that can be harvested from the population without affecting long term stability (average population size). The harvest within the harvestable surplus is called compensatory mortality, where the harvest deaths are substituting for the deaths that would otherwise occur naturally. Harvest beyond that is additive mortality, harvest in addition to all the animals that would have died naturally.

Care is needed when applying population dynamics to real world fisheries. Over-simplistic modelling of fisheries has resulted in the collapse of key stocks.[5][6]

History

The first principle of population dynamics is widely regarded as the exponential law of Malthus, as modelled by the Malthusian growth model. The early period was dominated by demographic studies such as the work of Benjamin Gompertz and Pierre François Verhulst in the early 19th century, who refined and adjusted the Malthusian demographic model. A more general model formulation was proposed by F.J. Richards in 1959, by which the models of Gompertz, Verhulst and also Ludwig von Bertalanffy are covered as special cases of the general formulation.[7]

Population size

The population size (usually denoted by N) is the number of individual organisms in a population.

The effective population size (Ne) was defined by Sewall Wright, who wrote two landmark papers on it (Wright 1931, 1938). He defined it as "the number of breeding individuals in an idealized population that would show the same amount of dispersion of allele frequencies under random genetic drift or the same amount of inbreeding as the population under consideration". It is a basic parameter in many models in population genetics. Ne is usually less than N (the absolute population size).

Small population size results in increased genetic drift. Population bottlenecks are when population size reduces for a short period of time.

Overpopulation may indicate any case in which the population of any species of animal may exceed the carrying capacity of its ecological niche.

Virtual population analysis

Virtual population analysis (VPA) is a cohort modeling technique commonly used in fisheries science for reconstructing historical fish numbers at age using information on death of individuals each year. This death is usually partitioned into catch by fisheries and natural mortality. VPA is virtual in the sense that the population size is not observed or measured directly but is inferred or back-calculated to have been a certain size in the past in order to support the observed fish catches and an assumed death rate owing to non-fishery related causes.

Minimum viable population

The minimum viable population (MVP) is a lower bound on the population of a species, such that it can survive in the wild. More specifically MVP is the smallest possible size at which a biological population can exist without facing extinction from natural disasters or demographic, environmental, or genetic stochasticity.[8] The term "population" refers to the population of a species in the wild.

As a reference standard, MVP is usually given with a population survival probability of somewhere between ninety and ninety-five percent and calculated for between one hundred and one thousand years into the future.

The MVP can be calculated using computer simulations known as population viability analyses (PVA), where populations are modelled and future population dynamics are projected.

Maximum sustainable yield

In population ecology and economics, the maximum sustainable yield or MSY is, theoretically, the largest catch that can be taken from a fishery stock over an indefinite period.[9][10] Under the assumption of logistic growth, the MSY will be exactly at half the carrying capacity of a species, as this is the stage at when population growth is highest. The maximum sustainable yield is usually higher than the optimum sustainable yield.

This logistic model of growth is produced by a population introduced to a new habitat or with very poor numbers going through a lag phase of slow growth at first. Once it reaches a foothold population it will go through a rapid growth rate that will start to level off once the species approaches carrying capacity. The idea of maximum sustained yield is to decrease population density to the point of highest growth rate possible. This changes the number of the population, but the new number can be maintained indefinitely, ideally.

MSY is extensively used for fisheries management.[11][12] Unlike the logistic (Schaefer) model, MSY in most modern fisheries models occurs at around 30% of the unexploited population size. This fraction differs among populations depending on the life history of the species and the age-specific selectivity of the fishing method.

However, the approach has been widely criticized as ignoring several key factors involved in fisheries management and has led to the devastating collapse of many fisheries. As a simple calculation, it ignores the size and age of the animal being taken, its reproductive status, and it focuses solely on the species in question, ignoring the damage to the ecosystem caused by the designated level of exploitation and the issue of bycatch. Among conservation biologists it is widely regarded as dangerous and misused.[5][6]

Recruitment

Recruitment is the number of new young fish that enter a population in a given year. The size of fish populations can fluctuate by orders of magnitude over time, and five to 10-fold variations in abundance are usual. This variability applies across time spans ranging from a year to hundreds of years. Year to year fluctuations in the abundance of short lived forage fish can be nearly as great as the fluctuations that occur over decades or centuries. This suggests that fluctuations in reproductive and recruitment success are prime factors behind fluctuations in abundance. Annual fluctuations often seem random, and recruitment success often has a poor relationship to adult stock levels and fishing effort. This makes prediction difficult.[13]

The recruitment problem is the problem of predicting the number of fish larvae in one season that will survive and become juvenile fish in the next season. It has been called "the central problem of fish population dynamics"[14] and “the major problem in fisheries science".[15] Fish produce huge volumes of larvae, but the volumes are very variable and mortality is high. This makes good predictions difficult.[16]

According to Daniel Pauly,[15][17] the definitive study was made in 1999 by Ransom Myers.[18] Myers solved the problem "by assembling a large base of stock data and developing a complex mathematical model to sort it out. Out of that came the conclusion that a female in general produced three to five recruits per year for most fish.”[15]

Overfishing

The Traffic Light colour convention, showing the concept of Harvest Control Rule (HCR), specifying when a rebuilding plan is mandatory in terms of precautionary and limit reference points for spawning biomass and fishing mortality rate.
Main article: Overfishing

The notion of overfishing hinges on what is meant by an acceptable level of fishing.

A current operational model used by some fisheries for predicting acceptable levels is the Harvest Control Rule (HCR). This formalizes and summarizes a management strategy which can actively adapt to subsequent feedback. The HCR is a variable over which the management has some direct control and describes how the harvest is intended to be controlled by management in relation to the state of some indicator of stock status. For example, a harvest control rule can describe the various values of fishing mortality which will be aimed at for various values of the stock abundance. Constant catch and constant fishing mortality are two types of simple harvest control rules.[19]

Metapopulation

Main article: Metapopulation

A metapopulation is a group of spatially separated populations of the same species which interact at some level. The term was coined by Richard Levins in 1969. 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".[20]

A metapopulation generally consists of several distinct populations together with areas of suitable habitat which are currently unoccupied. 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 prone it is to extinction.

Although individual populations have finite life-spans, the population 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).

Age class structure

Main article: Age class structure

Age can be determined by counting growth rings in fish scales, otoliths, cross-sections of fin spines for species with thick spines such as triggerfish, or teeth for a few species. Each method has its merits and drawbacks. Fish scales are easiest to obtain, but may be unreliable if scales have fallen off of the fish and new ones grown in their places. Fin spines may be unreliable for the same reason, and most fish do not have spines of sufficient thickness for clear rings to be visible. Otoliths will have stayed with the fish throughout its life history, but obtaining them requires killing the fish. Also, otoliths often require more preparation before ageing can occur.

An age class structure with gaps in it, for instance a regular bell curve for the population of 1-5 year-old fish, excepting a very low population for the 3-year-olds, implies a bad spawning year 3 years ago in that species.

Often fish in younger age class structures have very low numbers because they were small enough to slip through the sampling nets, and may in fact have a very healthy population.

Population cycle

Main article: Population cycle

A population cycle occurs where populations rise and fall over a predictable period of time. There are some species where population numbers have reasonably predictable patterns of change although the full reasons for population cycles is one of the major unsolved ecological problems. There are a number of factors which influence population change such as availability of food, predators, diseases and climate.

Trophic cascades

Main article: Trophic cascade

Trophic cascades occur when predators in a food chain suppress the abundance of their prey, thereby releasing the next lower trophic level from predation (or herbivory if the intermediate trophic level is an herbivore). For example, if the abundance of large piscivorous fish is increased in a lake, the abundance of their prey, zooplanktivorous fish, should decrease, large zooplankton abundance should increase, and phytoplankton biomass should decrease. This theory has stimulated new research in many areas of ecology. Trophic cascades may also be important for understanding the effects of removing top predators from food webs, as humans have done in many places through hunting and fishing activities.

Classic examples
  1. In lakes, piscivorous fish can dramatically reduce populations of zooplanktivorous fish, zooplanktivorous fish can dramatically alter freshwater zooplankton communities, and zooplankton grazing can in turn have large impacts on phytoplankton communities. Removal of piscivorous fish can change lake water from clear to green by allowing phytoplankton to flourish.[21]
  2. In the Eel River, in Northern California, fish (steelhead and roach) consume fish larvae and predatory insects. These smaller predators prey on midge larvae, which feed on algae. Removal of the larger fish increases the abundance of algae.[22]
  3. In Pacific kelp forests, sea otters feed on sea urchins. In areas where sea otters have been hunted to extinction, sea urchins increase in abundance and decimate kelp[23]

A recent theory, the mesopredator release hypothesis, states that the decline of top predators in an ecosystem results in increased populations of medium-sized predators (mesopredators).

Basic models

 \frac{dN}{dt} = r N \left(1 - \frac {N}{K} \right)
where N(t) represents number of individuals at time t, r the intrinsic growth rate and K is the carrying capacity, or the maximum number of individuals that the environment can support.
L'(t) = r_B \left( L_\infty - L(t) \right)
where rB is the von Bertalanffy growth rate and L the ultimate length of the individual.
H(E,X)=q E X\!
where the variables are; H, referring to catch (harvest) over a given period of time (e.g. a year); E, the fishing effort over the given period; X, the fish stock biomass at the beginning of the period (or the average biomass), and the parameter q represents the catchability of the stock.
Assuming the catch to equal the net natural growth in the population over the same period (\dot{X}=0), the equilibrium catch is a function of the long term fishing effort E:
H(E)=q K E (1-\frac{qE}{r})
r and K being biological parameters representing intrinsic growth rate and natural equilibrium biomass respectively.
C = \frac{F}{F+M} (1-e^{-(F+M)T}) N_0
where T is the time period and is usually left out (i.e. T=1 is assumed). The equation assumes that fishing and natural mortality occur simultaneously and thus "compete" with each other. The first term expresses the proportion of deaths that are caused by fishing, and the second and third term the total number of deaths.[25]
N_{t+1} = N_t e^{r(1-\frac{N_t}{k})}
Here r is interpreted as an intrinsic growth rate and k as the carrying capacity of the environment. The Ricker model was introduced in the context of the fisheries by Ricker (1954).
n_{t+1} = \frac{R_0 n_t}{1+ n_t/M}.
Here R0 is interpreted as the proliferation rate per generation and K = (R0  1) M is the carrying capacity of the environment.
{dN \over dt} = aN^2 - bN
where N is the size of a population, a is a half of the average probability of a birth of a male (the same for females) of a potential arbitrary parents pair within a year, and b is an average probability of a death of a fish within a year.

Predator-prey equations

The classic predator-prey equations are a pair of first order, non-linear, differential equations used to describe the dynamics of biological systems in which two species interact, one a predator and one its prey. They were proposed independently by Alfred J. Lotka in 1925 and Vito Volterra in 1926.

An extension to these are the competitive Lotka-Volterra equations, which provide a simple model of the population dynamics of species competing for some common resource.

In the 1930s Alexander Nicholson and Victor Bailey developed a model to describe the population dynamics of a coupled predator-prey system. The model assumes that predators search for prey at random, and that both predators and prey are assumed to be distributed in a non-contiguous ("clumped") fashion in the environment.[26]

In the late 1980s, a credible, simple alternative to the Lotka-Volterra predator-prey model (and its common prey dependent generalizations) emerged, the ratio dependent or Arditi-Ginzburg model.[27] 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.[28]

See also

Notes

  1. Wilderbuera, Thomas K and Zhang, Chang Ik (1999) Evaluation of the population dynamics and yield characteristics of Alaska plaice, Pleuronectes quadrituberculatus, in the eastern Bering Sea. Fisheries Research. Volume 41, Issue 2.
  2. Richard W Zabel, Chris J Harvey, Steven L Katz, Thomas P Good, Phillip S Levin (2003) Ecologically Sustainable Yield. American Scientist, March–April.
  3. Kapur V, Troy D and Oris J (1997) A Sustainable Fishing Simulation Using Mathematical Modeling Crossroads.
  4. Caswell, H. 2001. Matrix population models: Construction, analysis and interpretation, 2nd Edition. Sinauer Associates, Sunderland, Massachusetts. ISBN 0-87893-096-5.
  5. 5.0 5.1 Larkin PA (1977) "An epitaph for the concept of maximum sustained yield" Transactions of the American Fisheries Society, 106: 1–11.
  6. 6.0 6.1 Walters C and Maguire J (1996) "Lessons for stock assessment from the northern cod collapse", Reviews in Fish Biology and Fisheries, 6:125–137.
  7. Richards F. J. (1959) "A Flexible Growth Function for Empirical Use", Journal of Experimental Botany, 10: 290-301.
  8. Holsinger (2007)
  9. Europa: European Union (2006) Management based on maximum sustainable yield
  10. Europa: European Union (2006) Questions and Answers on Maximum Sustainable Yield (MSY)
  11. WWF Publications (2007) The Maximum Sustainable Yield objective in Fisheries
  12. New Zealand Ministry of Fisheries. MSY Harvest Strategies.
  13. Houde DE (2009) Recruitment variability in Jakobsen T, Fogarty MJ, Megrey BA and Moksness E. Fish Reproductive Biology, Chapter 3, Wiley-Blackwell. ISBN 978-1-4051-2126-2.
  14. Beyer, J.E. (1981) Aquatic ecosystems-an operational research approach. University of Washington Press. ISBN 0-295-95719-0.
  15. 15.0 15.1 15.2 Ransom A. Myers, 54, Dies; Expert on Loss of Fish Stocks The New York Times, 29 March 2007.
  16. Bakun A (1985) "Comparative studies and the recruitment problem: Searching for generalizations" CalCOFI Report Vol 26.
  17. Ransom Aldrich Myers (1952-2007) by Daniel Pauly, Nature, 10 May 2007.
  18. Myers R.A. (1995) "Recruitment of marine fish: the relative roles of density-dependent and density-independent mortality in the egg, larval, and juvenile stages" Marine Ecology Progress Series, 128: 305-310
  19. Coad, Brian W.; McAllister, Don E. (2008). "Dictionary of Ichthyology". Retrieved 2 January 2015.
  20. Levins, R. 1969. "Some demographic and genetic consequences of environmental heterogeneity for biological control." Bulletin of the Entomological Society of America, 15, 237-240
  21. Carpenter SR, Kitchell JF, Hodgson JR (1985) Cascading trophic interactions and lake productivity. Bioscience 35:634-639
  22. Power ME (1990) Effects of fish in river food webs. Science 250: 811-814
  23. Estes JA, Palmisano JF (1974) Sea otters: their role in structuring nearshore communities. Science 185: 1058-1060
  24. Quinn, Terrance J. II (2003). "Ruminations on the development and future of population dynamics models in fisheries". Natural Resource Modeling 16 (4): 341–392. doi:10.1111/j.1939-7445.2003.tb00119.x.
  25. Baranov, FI (1918) "On the question of the biological basis of fisheries" Izvestiya, 1: 81–128. (Translated from Russian by W.E. Ricker, 1945)
  26. Hopper
  27. Arditi, R. and Ginzburg, L.R. 1989. Coupling in predator-prey dynamics: ratio dependence. Journal of Theoretical Biology 139: 311-326.
  28. Arditi, R. and Ginzburg, L.R. 2012. How Species Interact: Altering the Standard View on Trophic Ecology. Oxford University Press, New York, NY.

References

Further reading

External links