Journal of Applied Mathematics

Volume 2014 (2014), Article ID 285350, 14 pages

http://dx.doi.org/10.1155/2014/285350

Research Article

Stability Analysis of Competing Insect Species for a Single Resource

Sizah Mwalusepo,1,2 Henri E. Z. Tonnang,¹ Estomih S. Massawe,² Tino Johansson,³ and Bruno Pierre Le Ru1,4,5

¹African Insect Science for Food and Health (icipe), P.O. Box 30772-00100, Nairobi, Kenya

²Department of Mathematics, University of Dar es Salaam, P.O. Box 35062, Dar es Salaam, Tanzania

³CHIESA Project, African Insect Science for Food and Health (icipe), P.O. Box 30772-00100, Nairobi, Kenya

⁴Unité de Recherche IRD 072, Laboratoire Evolution, Génomes et Spéciation, UPR 9034, 22 CNRS, 91198 Gif/Yvette, France

⁵Université Paris-Sud 11, 91405 Orsay, France

Received 22 February 2014; Accepted 12 June 2014; Published 6 July 2014

Academic Editor: Carlos Conca

Copyright © 2014 Sizah Mwalusepo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The models explore the effects of resource and temperature on competition between insect species. A system of differential equations is proposed and analysed qualitatively using stability theory. A local study of the models is performed around axial, planar, and interior equilibrium points to successively estimate the effect of (i) one species interacting with a resource, (ii) two competing species for a single resource, and (iii) three competing species for a single resource. The local stability analysis of the equilibrium is discussed using Routh-Hurwitz criteria. Numerical simulation of the models is performed to investigate the sensitivity of certain key parameters. The models are used to predict population dynamics in the selected cases studied. The results show that when a single species interacts with a resource, the species will be able to establish and sustain a stable population. However, in competing situation, it is observed that the combinations of three parameters (half-saturation, growth rate, and mortality rate) determine which species wins for any given resource. Moreover, our results indicate that each species is the superior competitor for the resource for the range of temperature for which it has the lowest equilibrium resource.

1. Introduction

During the past decades the theory of competition has played an increasingly important role in both development of general ecological theory and the interpretation of field data [1]. Competition is one of the primary biotic factors that shape patterns of distribution, abundance, and diversity in ecological communities [2]. The prominent status of competition in modern ecology is undoubtedly linked to views of pioneers in ecology. They strongly advocate the central role of competitive interactions between species [3–5]. The term “species” in the present study represents insect. Plant-feeding insects have become focal organisms for studies on population dynamics and community diversity because competition is deeply ingrained in their activities. It became accepted as a factor organizing insect assemblages, despite the fact that empirical evidence supporting competition was mostly derived from plant, vertebrate, and marine systems [6, 7]. In this circumstance the importance of interspecific competition is a highly controversial and unresolved issue for community ecology in general and for phytophagous insect in particular. Many theoretical studies of interspecific competition have dealt with the dependence of competition on the resource utilization abilities of each species [8–11]. The models used in these studies usually belong to two groups; the first group uses the classical theory of ecological competition between two or more species [12–14] that is attributed to Lotka [15] and Volterra [16]. These models are often “phenomenological” because the competition parameters are not independently derived values that can allow direct prediction of coexistence [17]. The second group employs resource-based theory in ecological competition; it considers the dynamic of the resource explicitly as well as the population dynamics of the competing species [8, 9, 17, 18]. The simplest form of resource-based competition occurs in laboratory apparatus, called a chemostat or continuous culture [3, 11]. In comparison with the classical models, the resource-based models may be less general and more difficult to analyse [8, 19, 20]. However, the resource-based models are often more predictive because the parameters can be measured on species alone [21]. Although competition among species is often cited as a major determinant of natural distributional patterns, there have been few studies in which knowledge of the mechanism of competition is used to predict these patterns [17]. The many biotic and abiotic factors that influence competition and possible interactions among them make such prediction difficult.

Climate, especially temperature, has a strong and direct effect on insect development, reproduction, and survival [22, 23]. This abiotic factor can be used to define ecological suitability of insect species and can thus dictate composition of pest communities in different agroecological regions. It affects life histories, insect distributions, interactions with host-plants, and their antagonists [22, 23]. However, this important parameter (temperature) has not much been incorporated in the theoretical studies relating resource and insect competition. Therefore, the objective of the present work is to determine and analyse the stability of the equilibrium points of insect species competing for a single resource while accounting for temperature through qualitatively analysis. We analyse the system as a combination of three submodels: (i) a single species with a single resource, (ii) two species competing for one resource, and (iii) three species competing for one resource.

2. Model Formulation

To formulate the model, the following assumptions are taken into consideration.

(i)The growth of the insect species is limited by a single resource.(ii)All insect species require the same resource to survive and reproduce.(iii)The resource availability depends on the rate of resource supply.(iv)No biological control measure is applied; that is, no parasitoid is applied.(v)The dynamic is assumed to follow “Monod’s” type function.(vi)All values of the competition coefficient are nonnegative.(vii)The growth rate is a function of both resource and temperature.

Taking into account the above considerations, then we have the schematic flow charts (Figure 1).

(a)

(b)

(c)

Figure 1: Flow charts summarizing the conceptual framework of the formulated models; the boxes represent state variables, which describe the state or condition of the system components. The arrows illustrate relationships among state variables as the movement of ecological interactions; (a) shows one species and one resource interactions, (b) represents two species competing for one resource, and (c) illustrates three species competing for one resource.

From the flow charts, three systems of ordinary differential equations governing the growth and the dependence of the species to the resource were deduced.

2.1. Case 1: One Species and One Resource Model

Consider

2.2. Case 2: Two Species and One Resource Model

Consider

2.3. Case 3: Three Species and One Resource Model

Consider

Let , in these systems.

The Michaelis-Menten equation is employed to describe the relationship between cellular metabolism and resource. For every large value of the resource, this equation gives approximately a value of 1. For a resource equal to , it gives approximately a value of 1/2, which justifies why is referred to as the half-saturation constant. To be biologically relevant, the following conditions must be satisfied: and . The variables and parameters used in the models are described in Abbreviation Section.

3. Models Analysis

The mathematical models equations (1)–(3) are analysed qualitatively. An elaborative analysis of the equilibrium is performed and conditions for existence of the proposed systems are provided. The stability analysis of the equilibrium is studied by computing the variational matrices for each equilibrium point. However, the local stability is established through Jacobian matrix of the systems using the Routh-Hurwitz criteria. For linearized systems, the Jacobian matrix is given by

This is associated with the characteristic equation:

In general, the characteristic equation is a polynomial expression with degree equal to the number of species competing, is the eigenvalues of the linearized system, and is the identity matrix.

Theorem 1 (Routh-Hurwitz criteria). Given the polynomial , where the coefficients are real constants, , define the Hurwitz matrices using the coefficients of the characteristic polynomial. Then, all of the roots of the polynomial are negative or have negative real parts if and only if the determinants of all Hurwitz matrices are positive. For polynomial of degree and 4, the Hurwitz criteria are summarized as follows: if : :  , : , and .

Theorem 2. For an matrix the determinant is given by the explicit formula , where is the permutation group, is the signum function of permutation in the permutation group, which returns and for even and odd permutation, respectively, and represents the value in th position after the reordering.

3.1. Models Implementation and Population Size Predictions

The calculations are carried out with computer program written in . In the program, model equations are solved using the 4th order Runge-Kutta algorithm with step size of 0.01. Once the models equations are solved, graphs are generated displaying predictions of future species population size at fixed initial number of the resource. For a random selection of species, initial population size, the models are used to estimate the values of species at a given time interval; that may be day, week, or month depending on the species life cycle.

3.2. Qualitative Analysis of the Developed Models

3.2.1. One Species and One Resource

The model for this case has at most two equilibrium points, namely, axial equilibrium point and interior equilibrium point . The asterisk on and indicates that this is a nonzero equilibrium value of and , respectively. These equilibrium points are obtained when , and . The local stability of each equilibrium point of model system (1) is ascertained by the following theorem.

Theorem 3. (a) By using the Routh-Hurwitz criterion, is locally asymptotically stable if

(b) The interior point is locally asymptotically stable if

3.2.2. Two Species One Resource

Model (2) has three equilibrium points, which are

Theorem 4. Through Routh-Hurwitz criteria, the equilibrium points , , and will be locally asymptotically stable if the following conditions hold: for the polynomial degree , the equilibrium points can be determined without explicitly solving the characteristic polynomial   +    +    +   . Then if ; otherwise, the equilibrium points will be unstable.

3.2.3. Three Species One Resource

For this case, the system of ordinary differential equations has four equilibrium points:

Theorem 5. The equilibrium point is locally asymptotically stable if (i) , (ii) , (iii) , and (iv) ; otherwise, it will be unstable.

Proof. The determinant of the Jacobian matrix at is given by

Let

Then the characteristic equation of the matrix is calculated as

where

Applying the Routh-Hurwitz criterion applied on a 4-by-4 matrix, will be locally asymptotically stable if ,  ,  , and and unstable otherwise.

Theorem 6. The planar equilibrium points, , , and are locally asymptotically stable if (i) , (ii) , (iii) , and (iv) .

Proof. The Jacobian matrix at is given by

where

The determinant of is given by

The characteristic equation is

where

Applying the Routh-Hurwitz criterion, the planar equilibrium point will be locally asymptotically stable if (i) , (ii) , (iii) , and (iv) .

Similarly, for and , again through the Routh-Hurwitz criteria, we have the following: from the polynomial , with if , and the system will be conditionally stable.

3.3. Resource Response to Temperature

It is useful to determine the quantitative and qualitative behaviour of planar equilibrium in response to temperature, since the quantity of resource can highly vary within any particular region of qualitatively similar parameters space. In this effect, let us consider the resource as a function of temperature with the following expression:

where are constants.

We determine the response of to temperature using the derivate .

The derivative of the equilibrium of the temperature dependent resource function with respect to temperature is a dampened exponential function given by

This generates a rather simple relationship governing the directional response of equilibrium resource to temperature variation. For numerical simulation, the parameters are chosen to give the behaviour of the isoclines for the species and the resource. This is done using three different sets of parameters (unequal growth rate, equal growth rate, and equal consumption rates). Due to the lack of parameters values in the literature and from the field study we assumed the following values in Table 1 to simulate the model.

Table 1: Parameters used in simulations.

Parameters
Species one
Species two
Species three

0.1
0.17
0.2

0.02
0.02
0.02

0.01
0.01
0.01

70
70
70

0.02
0.1
0.7

3
0.5
0.4

7
7
7

22
20
23

4. Results

This section presents the results of the numerical simulations that were carried out to investigate the dynamical behaviour of the formulated models. In order to verify the theoretical predications of the models, the numerical simulations are carried out using fixed values of the model parameters of the models. These parameters values are chosen for simulation purposes only. However, these parameters may be biologically feasible. Numerical simulations for three species are done singly at five different temperatures (10°C, 15°C, 20°C, 25°C, and 30°C) and the results are represented in Figures 2, 3, 4, 5, 6, 7, 8, and 9.

(a)

(b)

Figure 2: Behavior of the isoclines for one species and the resource at different temperatures; (a) shows the behavior of population size and the resource at temperatures 10°C, 15°C, 20°C, 25°C, and 30°C, respectively; (b) shows the behavior of population size and resource at different initial population size at temperature 25°C.

(a)

(b)

(c)

Figure 3: Behavior of the isoclines for two species competing for a single resource when both species initial population size are the same; curves (a)–(c) show the behavior of population size for each species and resource at temperatures: (a) 10°C and 15°C, (b) 20°C and 25°C, and (c) and 30°C, respectively.

(a)

(b)

(c)

Figure 4: Behavior of the isoclines when the population size of species one is greater than species two; curves (a)–(c) show the behavior of population size for each species and the resource at temperatures: (a) 10°C and 15°C, (b) 20°C and 25°C, and (c) at 30°C, respectively.

(a)

(b)

(c)

Figure 5: Behavior of the isoclines when the population size of species two is greater than species one; curves (a)–(c) show the behavior of population size for each species and the resource at temperatures: (a) 10°C and 15°C, (b) 20°C and 25°C, and (c) at 30°C, respectively.

(a)

(b)

(c)

(d)

(e)

Figure 6: Behavior of the isoclines for three species competing for a single resource when both species initial population size are the same; curves (a)–(e) show the behavior of population size for each resource at temperatures 10°C, 15°C, 20°C, 25°C, and 30°C, respectively. The solid lines represent the resource; species one are represented with round dot lines, species two are represented with dash lines, and species three are represented with long dash dot lines and all species have different growth rate.

(a)

(b)

(c)

(d)

(e)

Figure 7: Behavior of the isoclines when initial, population size of species one is greater than species two and three; curves (a)–(e) show the behavior of population size for each species and resource at temperatures 10°C, 15°C, 20°C, 25°C, and 30°C, respectively.

(a)

(b)

(c)

(d)

(e)

Figure 8: Behavior of the isoclines when the population size of species two is greater than that of species one and three; curves (a)–(e) show the behavior of population size for each species and resource at temperatures 10°C, 15°C, 20°C, 25°C, and 30°C, respectively.

(a)

(b)

(c)

(d)

(e)

Figure 9: Behavior of the isoclines when the population size of species three is greater than species two and one; curves (a)–(e) show the behavior of population size for each species and the resource at temperatures 10°C, 15°C, 20°C, 25°C, and 30°C, respectively.

At 10°C (Figure 2(a)), for a single species, the results show a slow increase in the curve of the species population size and a slow decrease in the resource quantity. At 15°C, 20°C, 25°C, and 30°C, respectively, similar curves behaviour shape is obtained. It starts with a rapid increase of the population size of the species, followed by slow decreases and a rapid decrease of the resource. Figure 2(b) shows the relationship between the different initial population size (10, 20, and 30) with the resource at 25°C; it was observed that the more the initial population size of the species the resource quantity decrease faster.

For two species competing for a single resource (Figures 3–5), the results show that species one is a superior competitor at all temperatures tested when the initial population size of species one is greater than species two and within 20–30°C when both species started at the same initial population size. Species two is a superior competitor at all temperatures tested when the initial population size of species two is greater than species one and at 10°C when both species started at the same initial population size. No difference is observed between the two species at 15°C when the initial population size of the two species is identical.

In the case of three species competing for a single resource, the species with higher initial population size is the superior competitor at all temperatures. However, when all species have equal initial population size, species two grew faster than its two counterparts at lower temperatures (from 10°C to 15°C). Between 20°C and 25°C species one is superior competitor. With the same growth rates, the species with the lower mortality is the superior competitor when all species have equal initial population size. It is observed that one species is the superior competitor at 10°C to 15°C and with increase in temperature (°C) another species became superior when all species have equal consumption rates.

Figure 10 represents a combination of temperature and resource for which each species has no net population growth for a given mortality rate. Herein, the zero net growth isoclines represent equilibrium boundaries of individual species based on their requirements in resource availability. The curves show the amount of the resource that each species needs at equilibrium to have a reproductive rate that can balance its mortality rate. In general, within the unshaded region, the population size of the species will increase, while it would decrease in the shaded region.

(a)

(b)

(c)

(d)

Figure 10: Temperature and equilibrium resource isoclines; curves (a)–(c) represent equilibrium boundaries of each species based on the requirements of equilibrium resource, within the unshaded region, the population size of the species will increase when it would decrease in shaded region, and (d) illustrated when three species are superimposed.

5. Discussion

We have conducted a qualitative analysis of the effects of resource and temperature on competition between insect species. For this purpose, temperature dependent system of differential equations was formulated. Both qualitative and numerical analyses of the models were considered. In respect to each model, the axial, planar, and interior equilibrium points were obtained and their stabilities were investigated. Conditions on the parameters that insure the stability and instability of the boundary equilibrium points were presented. For a single species interacting with the resource, the local stability analyses revealed that the species could maintain a stable population when growing alone under all conditions used in the analysis. It was illustrated that as the population size of the species grows, the resource depletes to the level of equilibrium of the resource availability where the population size of the species stops growing. Numerical simulation also revealed that the dynamical outcomes of the interactions are very sensitive to the model parameter values and to the initial population size of the species. In axial equilibrium point, it was observed that if the value of the equilibrium of the resource availability is too large than the supply point of the resource, the population of the species cannot persist. In one hand, this is due to the fact that the resource availability directly depends on the rate of resource supply. On the other hand, if the supply of the resource is sufficiently large, the species can establish a stable population and equally control the resource to the equilibrium of the resource availability.

Generally, the obtained results highlight an essential assumption of the resource competition theory, which suggested that when several species compete for the same resource, the species with the lowest requirement for the resource is the superior competitor at equilibrium [24]. This species competitively displaces all its counterparts, independent of the initial size of the population for other competing species. This means that at the equilibrium points, in a habitat with a single limiting resource, the best competitor is the species that has the lowest equilibrium resource value . It implies that a given species is able to continue growing as long as the resource is greater than the lowest equilibrium resource value . Reciprocally, species can go to local extinction if the resource is less than the lowest equilibrium resource . For these values of , the resource-dependent growth rates of the species is less than their mortality rates, and all species will experience negative exponential decline towards the equilibrium point. These results are in agreement with Tilman findings reported by Tilman [8].

From the formulated models, it was observed that the relationships between the available resource and resulting population size are nonlinear. This makes the analysis of the developed mathematical equations in varying temperature difficult. We therefore opted to conduct the analysis at fixed temperatures and then generalize by intervals [25–28]. The study of insects’ physiology has demonstrated that different species have different temperature optimum that affects their competitive ability for resource exploitation [10]. To predict the outcome of competition between insect species within a range of temperatures, we superimpose the zero net growth isoclines of the species and results show that one species has an equilibrium resource significantly lower than its counterparts.

Furthermore, in a competing situation it is observed that three parameters, half-saturation constant, growth rate, and the mortality rate determine which species wins for any given resource. A species with a high affinity for the resource can still lose if it has low growth rate and a high death rate. This relationship is simple as it demonstrates that resource-based models bring more understanding than the classical theory in ecology for species resource competition [8, 11]. Moreover, the use of resource-based models emphasizes the role of resource in determining competitive interactions and the abilities of each species to acquire and use the limiting resource. Classical models of competition [14, 17] are useful tools for understanding insect population growth, whereas resource-based models create direct link between the species and the resource. These types of models are more predictive.

Further experimental or field investigation in this direction may be helpful to verify whether or not the process prescribed by the formulated models actually functions in a similar manner under practical consideration. For example, we should verify the accuracy of our models on the maize-lepidopteran stem borer pests system. In fact, maize is the most important staple food in most sub-Saharan Africa [29]. In spite of its importance, maize production is constrained by lepidopteran stem borer pests, the most damaging being the crambid Chilo partellus (Swinhoe) and the noctuid Busseola fusca (Fuller) and Sesamia calamistis Hampson [30]; more and in-depth research is needed on its predominant insect pests with different optimal temperatures to guide the formulation of appropriate adaptation strategies to the current climate change [31]. Therefore, our models can be a useful tool to test how temperature can affect the competition between these insect pests, their future distribution, and incidence on maize production.

In conclusion, the present study qualitatively analysed a series of mathematical equations to demonstrate competition of insect species feeding on one resource. Therefore, it would be interesting to conduct further quantitative and qualitative studies for all kind of insect communities with different feeding behaviour (boring insects, chewing insects, and phloem-feeding insects) competing for one resource for possible improvement of integrated pest management strategies. Furthermore, many questions remain open including applying or improving the model and obtaining the results when the resource and insect species are in different stages (neonates and old instars larvae). Thus, purely observational methods may give an indication of whether competition is occurring, particularly if time-series data are used.

Abbreviations

: | Amount of resource available at time
: | Population size of species at time
: | The physical factor (temperature)
: | The optimal temperature of species
: | Maximal gross productivity rate of species at temperature
: | Population growth rate of species
: | Consumption rate of the resource by species
: | Mortality rate of species
: | The half-saturation constant of species
: | The flow or supply rate of the resource
: | The supply point
: | The width environmental tolerance.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are grateful to the entire members of the Noctuid Stem Borer Biodiversity (NSBB) team of the African Insect Science for Food and Health (icipe), the Institut de Recherche pour le Développment (IRD), and University of Dar es salaam, through the Department of Mathematics for their support. This research was conducted under the Climate Change Impacts on Ecosystem Services and Food Security in Eastern Africa (CHIESA) project. CHIESA is funded by the Ministry for Foreign Affairs of Finland.

References

D. H. Robert, Predation, Apparent Competition, and the Structure of Prey Communities, Academic Press, New York, NY, USA, 1977.

I. Kaplan and R. F. Denno, “Interspecific interactions in phytophagous insects revisited: a quantitative assessment of competition theory,” Ecology Letters, vol. 10, no. 10, pp. 977–994, 2007.

G. S. K. Wolkowicz and Z. Lu, “Global dynamics of a mathematical model of competition in the chemostat: general response functions and differential death rates,” SIAM Journal on Applied Mathematics, vol. 52, no. 1, pp. 222–233, 1992.

A. J. Lotka, “The growth of mixed populations: two species competing for a common food supply,” Journal of Washington Academic of Science, vol. 22, no. 16, pp. 461–469, 1932.

G. F. Gause, The Struggle for Existence, Williams and Wilkins, Baltimore, Md, USA, 1934.

T. Jermy, “Is there competition between phytophagous insects?” Journal of Zoological Systematic and Evolutionary Research, vol. 23, no. 4, pp. 275–285, 1985.

R. F. Denno, M. S. McClure, and J. R. Ott, “Interspecific interactions in phytophagous insects: competition reexamined and resurrected,” Annual Review of Entomology, vol. 40, pp. 297–331, 1995.

D. Tilman, “Resource competition between planktonic algae: an experimental and theoretical approach,” Ecology, vol. 58, no. 2, pp. 338–348, 1977.

D. Tilman, “The importance of the mechanisms of interspecific competition,” The American Naturalist, vol. 129, no. 5, pp. 769–774, 1978.

J. P. Grover, Resource Competition, Chapman & Hall, London, UK, 1997.

I. P. Martines, H. V. Kojouharov, and J. P. Grover, “A chemostat model of resource competition and allelopathy,” Applied Mathematics and Computation, vol. 215, no. 2, pp. 573–582, 2009.

B. Liu and L. Zhang, “Dynamics of a two-species Lotka-Volterra competition system in a polluted environment with pulse toxicant input,” Applied Mathematics and Computation, vol. 214, no. 1, pp. 155–162, 2009.

O. A. Ebraheen, A. D. Fordyce, and D. Niall, “A 3-species competition model for bio-control,” Applied Mathematics and Computation, vol. 218, no. 18, pp. 9690–9698, 2012.

A. A. S. Zaghrout and F. M. Kandil, “Competition between three microbial populations for a single limiting resource in continuous culture,” Applied Mathematics and Computation, vol. 92, no. 2-3, pp. 271–281, 1998.

A. J. Lotka, Elementary of Physical Biology, Williams and Wilkins, Baltimore, Md, USA, 1925.

V. Volterra, “Fluctuations in the abundance of a species considered mathematically,” Nature, vol. 118, no. 2972, pp. 558–560, 1926.

D. Tilman, “Competition and biodiversity in spatially structured habitats,” Ecology, vol. 75, no. 1, pp. 2–16, 1994.

J. A. Leon and D. B. Tumpson, “Competition between two species for two complementary or substitutable resources,” Journal of Theoretical Biology, vol. 50, no. 1, pp. 185–201, 1975.

S. B. Hsu, “Limiting behavior for competing species,” SIAM Journal on Applied Mathematics, vol. 34, no. 4, pp. 760–763, 1978.

S. B. Hsu, S. Hubbell, and P. Waltman, “A mathematical theory for single-nutrient competition in continuous cultures of micro-organisms,” SIAM Journal on Applied Mathematics, vol. 32, no. 2, pp. 366–383, 1977.

S. R. Hansen and S. P. Hubbell, “Single-nutrient microbial competition: qualitative agreement between experimental and theoretically forecast outcomes,” Science, vol. 207, no. 4438, pp. 1491–1493, 1980.

J. S. Bale, G. J. Masters, I. D. Hodkinson et al., “Herbivory in global climate change research: direct effects of rising temperature on insect herbivores,” Global Change Biology, vol. 8, no. 1, pp. 1–16, 2002.

M. Ladányi and L. Horváth, “A review of the potential climate change impact on insect populations,” Applied Ecology and Environmental Research, vol. 8, no. 1, pp. 143–152, 2010.

W. J. O'Brien, “The dynamic of nutrient limitation of phytoplankton algae: a model reconsidered,” Ecology, vol. 55, no. 1, pp. 135–141, 1974.

R. Levins, “Coexistence in a variable environment,” The American Naturalist, vol. 114, no. 6, pp. 765–783, 1979.

P. Chesson, “Multispecies competition in variable environments,” Theoretical Population Biology, vol. 45, no. 3, pp. 227–276, 1994.

S. B. Hsu and P. Waltman, “Competition in the chemostat when one competitor produces a toxin,” Japan Journal of Industrial and Applied Mathematics, vol. 15, no. 3, pp. 471–490, 1998.

S. B. Hsu, Y.-S. Li, and P. Waltman, “Competition in the presence of a lethal external inhibitor,” Mathematical Biosciences, vol. 167, no. 2, pp. 177–199, 2000.

M. Smale, D. Byerlee, and T. Jayne, “Maize Revolution in Sub-Saharan Africa,” The World Bank, 2001

R. Kfir, W. A. Overholt, Z. R. Khan, and A. Polaszek, “Biology and management of economically important lepidopteran cereal stem borers in Africa,” Annual Review of Entomology, vol. 47, pp. 701–731, 2002.

R. James and R. Washington, “Changes in African temperature and precipitation associated with degrees of global warming,” Climatic Change, vol. 117, no. 4, pp. 859–872, 2013.