![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
This paper studies a 2-dimensional discrete-time competition model of Ricker type with reproductive delay. The model is examined under the assumption that species 1 and 2 have the same properties except that a fraction η of species 1 individuals delays the initiation of reproduction. This assumption ensures that species 1 is dominated by species 2 in the sense that species 2 is increasing whenever species 1 is increasing. It is shown that, even under this assumption, delayed reproduction can be adaptive, i.e. species 1 can invade the monoculture system of species 2 while species 2 cannot invade the monoculture system of species 1, if the population is fluctuating. The result is obtained by analytically examining the species invasibility at boundary 2-cycles, whose coordinates can be estimated by assuming .
1. Introduction
This paper studies the following discrete-time competition model:
where
and
denote the number of juveniles of species 1 and 2, respectively. It is assumed that both species are semelparous, i.e. individuals reproduce only once in their lifetime and die immediately after reproduction. The parameter s denotes the survival probability of juveniles. The most important assumption is that a fraction η of survived juveniles of species 1 postpones maturity and reproduction while all survived juveniles of species 2 mature and reproduce in unit time. The density dependence is assumed to act on reproduction, so the factor
represents the number of offspring produced by an adult individual. Here b>0 and
. It is intuitive that delayed reproduction is maladaptive since delayed reproduction increases the possibility of death before reproduction. However, as shown in several papers (e.g. see [Citation18]), delayed reproduction can be adaptive in variable environments. Although these studies assume external factors modelling variable environments, we show that delayed reproduction can be adaptive even in a simple nonlinear deterministic model that exhibits population fluctuation.
By writing ,
, and
, the above model equation becomes
(1)
(1) where the tilde is dropped for convenience. This system shall be studied in this paper under the condition that
,
. It is clear that this system is a special case of the following system
(2)
(2) where
,
,
,
, and β are positive. This competition system was studied by Yakubu [Citation19], who considered that the first terms of the right-hand sides of both equations represent the density dependent self-reproduction while the term
indicates that species 1 is planted (or stocked) at a fixed rate and found several important properties of system (Equation2
(2)
(2) ). To understand these properties, we need to recall the dynamics of the following competition model of Ricker type:
(3)
(3) where
,
,
,
,
, and
are positive. System (Equation2
(2)
(2) ) with
is a special case of system (Equation3
(3)
(3) ). This system has been studied in many papers (e.g. see [Citation5,Citation10,Citation11,Citation15]) since Hassell and Comins [Citation8] to understand competition between two species whose reproductions occur discretely in time. Franke and Yakubu [Citation6, Theorem 3.1] studied system (Equation3
(3)
(3) ) as a special case of
(4)
(4) and showed that extinction of species i occurs in system (Equation3
(3)
(3) ) irrespective of initial conditions if species i is dominated by the other species in the sense that the other species is increasing whenever species i is increasing. More precisely, species i is said to be dominated by species j if
is a proper subset of
, where
and
. In the case of system (Equation2
(2)
(2) ) with
, except for the special case
, either species 1 or 2 is dominated by the other. Thus, under the assumption that species 1 is dominated by species 2 in system (Equation2
(2)
(2) ) with
, Yakubu [Citation19] examined the influence of planting of species 1 on its persistence and revealed that two species of system (Equation2
(2)
(2) ) can coexist at a stable positive 2-cycle. The general result of Franke and Yakubu [Citation5, Theorem 6.1] on system (Equation4
(4)
(4) ) shows that if each single-species dynamics within the coordinate axes has a globally attractive positive fixed point, then the dominated species goes extinct irrespective of initial conditions. Therefore, for species coexistence in system (Equation2
(2)
(2) ) oscillation of single-species population is crucial. This result is fascinating since coexistence can occur even though system (Equation2
(2)
(2) ) has no positive fixed points.
After the paper [Citation19], Yakubu [Citation20] posed open problems concerning the global attractivity of a positive 2-cycle of (Equation2(2)
(2) ) and the dynamics of system (Equation2
(2)
(2) ) and its variant have been studied in several papers. By applying a general result, Elaydi and Yakubu [Citation4] showed that the positive 2-cycle of (Equation2
(2)
(2) ) cannot be globally attracting in the interior of
. Kon [Citation14] showed that species 2 goes extinct in system (Equation2
(2)
(2) ) irrespective of initial conditions if species 2 is dominated by species 1. Kang and Smith [Citation13] studied a variant of system (Equation2
(2)
(2) ) and showed that the system has a positive 2-cycle that attracts all positive points except Lebesgue measure zero set (see also [Citation12]). Although these studies advanced our understanding of system (Equation2
(2)
(2) ), its dynamics when species 1 is dominated by species 2 is still not well understood. In this paper, by focussing on the stability of boundary 2-cycles of system (Equation1
(1)
(1) ), which is a special case of system (Equation2
(2)
(2) ), we reveal its dynamics and give mathematical evidence that delayed reproduction can be adaptive when the population is fluctuating.
The rest of this paper is organized as follows. Section 2 studies the dynamics of system (Equation1(1)
(1) ) on the boundary of
. We are concerned with boundary 2-cycles. Since the dynamics on the
-axis obeys the Ricker map, we review the analytical expression of its 2-cycles. By using the fact that the dynamics on the
-axis can be seen as a perturbation of the Ricker map, we obtain an approximate expression of 2-cycles on the
-axis by assuming
. Section 3 reviews the fixed-point stability of system (Equation1
(1)
(1) ) and confirms that species 1 is dominated by species 2. Section 4 examines the stability of the boundary 2-cycles constructed in Section 2 and analytically obtain the stability criterion in terms of parameters. By the result of Section 4, we find that delayed reproduction can be adaptive. Furthermore, we find that it is unlikely that two species can either mutually invade or mutually prevent invasion of each other at boundary 2-cycles if η is small. However, we numerically show that the parameter region for mutual invasion at boundary 2-cycles gets larger as η increases. In such a parameter region, we find that two species can coexist. Section 5 examines the global dynamics of (Equation1
(1)
(1) ) both analytically and numerically. We analytically show that it is unlikely that system (Equation1
(1)
(1) ) has a positive 2-cycle if η is small. Numerical simulations suggest that if system (Equation1
(1)
(1) ) has a stable boundary 2-cycle, then it attracts almost all positive points as long as
. Section 6 includes some concluding remarks.
2. Boundary dynamics and boundary 2-cycles
On the - and
-axes, system (Equation1
(1)
(1) ) is reduced to the one-dimensional maps
respectively. Note that
for all x if
.
The map is the well-known Ricker map (e.g. see [Citation16]). This map has two fixed points 0 and
. The inequality
always holds and
holds if and only if
. Thus, the origin is always unstable and the positive fixed point
is asymptotically stable if
and is unstable if
. It is known that, at the critical value of
, a period-doubling bifurcation occurs resulting in the birth of an asymptotically stable 2-cycle. Let
and
(
) be 2-periodic points of the 2-cycle. Then they have the following analytical expression [Citation17]:
(5)
(5) where
and h is the strictly increasing function defined by
(6)
(6) for
. Note that
as
. Since
is an increasing function of λ, the amplitude of the 2-cycle
increases as λ increases. The 2-cycle
is asymptotically stable if
, which is equivalent to
, where
(e.g. see [Citation3]).
The map has two fixed points 0 and
. The inequality
always holds and
holds if and only if
. Concerning a 2-cycle of
, we obtain the following lemma.
Lemma 2.1
The map has a positive 2-cycle
for all
sufficiently small if
. The 2-cycle
is asymptotically stable (resp. unstable) for all
sufficiently small if
(resp.
). The 2-periodic points have the following expression:
(7)
(7)
Proof.
Let and
(
) be 2-periodic points of the map
. Then they are given by solving the following equation for x:
where
. Since
is reduced to the Ricker map
when
, both
and
hold. Thus, by the implicit function theorem, the equation
has the solutions
(8)
(8) in neighbourhoods of
and
, respectively. It is straightforward to show that Equation (Equation8
(8)
(8) ) is equivalent to (Equation7
(7)
(7) ). Since the 2-cycle
of the map
exists if
, the 2-cycle
of the map
exists for all
sufficiently small if
. Furthermore, the 2-cycle
of the map
is asymptotically stable (resp. unstable) for all
sufficiently small if
(resp.
) since
,
, and
continuously depend on η in a neighbourhood of
and
holds if
.
3. Stability of fixed points
The fixed-point stability of system (Equation2(2)
(2) ) is studied by Yakubu [Citation19] (see also [Citation4]). In this section, we briefly review the fixed-point stability of a special case of system (Equation2
(2)
(2) ), i.e. (Equation1
(1)
(1) ). The fixed points of system (Equation1
(1)
(1) ) are given by solving the system of equations
(9)
(9) It is clear that system (Equation1
(1)
(1) ) has the three boundary fixed points
,
, and
, where
and
as defined in the previous section. Since
, the inequality
always holds. If both
and
are not zero, the first and second equations of (Equation9
(9)
(9) ) are reduced to the incompatible equations
and
, respectively. Thus, there are no positive fixed points. Furthermore,
implies that
is a proper subset of
, i.e. species 1 is dominated by species 2 (see Figure ).
Figure 1. The phase plane of system (Equation1
(1)
(1) ). The set
is a proper subset of
. Species 1 is dominated by species 2. The unstable manifold of
and the stable manifold of
have nonempty intersections with the interior of
.
![Figure 1. The phase plane (x1,x2) of system (Equation1(1) {x1(n+1)=x1(n){(1−η)exp(λ−x1(n)−x2(n))+ηs}x2(n+1)=x2(n)exp(λ−x1(n)−x2(n)),(1) ). The set D1 is a proper subset of D2. Species 1 is dominated by species 2. The unstable manifold of (X1,0) and the stable manifold of (0,X2) have nonempty intersections with the interior of R+2.](/cms/asset/1647794a-7aed-4bd7-952d-a529bb5dbe3a/tjbd_a_2248171_f0001_ob.jpg)
Let be the Jacobi matrix of system (Equation1
(1)
(1) ) evaluated at
. Then
has the eigenvalues
and
, both of which are larger than one for all
sufficiently small since
. Thus, the origin is a source for all
sufficiently small.
has the eigenvalues
and
. The former eigenvalue determines the stability of the fixed point
of the map
defined in the previous section. Note that
is destabilized at
as λ is increased. The latter eigenvalue is transversal in the sense that it has an eigenvector transversal to the
-axis. Thus the unstable manifold of
always has a nonempty intersection with the interior of
.
has the eigenvalues
and
. The former eigenvalue determines the stability of the fixed point
of the map
defined in the previous section. Note that
is destabilized at
as λ is increased. The latter eigenvalue is transversal. Thus the stable manifold of
always has a nonempty intersection with the interior of
. Figure summarizes some of the above results.
4. Stability of boundary 2-cycles
Concerning the stability of boundary 2-cycles of system (Equation1(1)
(1) ), we obtain the following theorem.
Theorem 4.1
If , then system (Equation1
(1)
(1) ) with
sufficiently small has the 2-cycles
and
.
(a) | If | ||||
(b) | If | ||||
(c) | If | ||||
(d) | If |
Proof.
Since system (Equation1(1)
(1) ) is reduced to the Ricker map on the
-axis, system (Equation1
(1)
(1) ) has the 2-cycle
for all
if
. Recall that
denotes the Jacobi matrix of system (Equation1
(1)
(1) ) evaluated at
. Then the Jacobi matrix of the second iterate of system (Equation1
(1)
(1) ) evaluated at
is given by
. This matrix has the eigenvalues
(10)
(10) As mentioned in Section 2, the absolute value of the former eigenvalue is less than one if
. The latter eigenvalue takes the form
Using (Equation5
(5)
(5) ), we find that the absolute value of this eigenvalue is less (resp. larger) than one for all
sufficiently small if
(resp.
). Therefore, the case where the 2-cycle
is hyperbolic can be divided into four cases: (a)
and
, (b)
and
, (c)
and
, and (d)
and
. The 2-cycle
is a sink in case (a), a saddle in cases (b) and (c), and a source in case (d). Furthermore, we find that
has an unstable manifold that intersects the interior of
in case (b) since the latter eigenvalue is transversal. Similarly, we find that
has a stable manifold that intersects the interior of
in case (c).
By Lemma 2.1, system (Equation1(1)
(1) ) has the 2-cycle
on the
-axis for all
sufficiently small if
. The Jacobi matrix of the second iterate of system (Equation1
(1)
(1) ) evaluated at
is given by
. This matrix has the eigenvalues
(11)
(11) As shown in Lemma 2.1, the absolute value of the former eigenvalue is less than one for all
sufficiently small if
. Equation (Equation5
(5)
(5) ) reduces the latter eigenvalue into
whose absolute value is less (resp. larger) than one for all
sufficiently small if
(reps.
) since
. Therefore, the case where the 2-cycle
is hyperbolic can be divided into four cases (a), (b), (c), and (d) as above. The 2-cycle
is a saddle in case (a), a sink in cases (b), a source in case (c), and a saddle in case (d). Furthermore, we find that
has an unstable manifold that intersects the interior of
in case (a) since the latter eigenvalue is transversal. Similarly, we find that
has a stable manifold that intersects the interior of
in case (d).
The inequalities and
can be rewritten as follows:
(12)
(12) In fact, since the function
defined by (Equation6
(6)
(6) ) is strictly increasing for x>0,
if and only if
and
if and only if
. Thus we obtain (Equation12
(12)
(12) ).
Note that is asymptotically stable for all
sufficiently small if
and
and is unstable for all
sufficiently small if either
or
and
is asymptotically stable for all
sufficiently small if
and
and is unstable for all
sufficiently small if either
or
.
The above result is summarized in Table . The table shows that each boundary fixed point or 2-cycle has the indicated stability if they satisfy the indicated parameter conditions. Note, however, that the fixed point and the 2-cycles
and
have the indicated stability if
is sufficiently small. Since
, we have the following expression:
Figure shows the parameter plane demarcated with the stability of the boundary fixed points and the boundary 2-cycles of system (Equation1
(1)
(1) ) with
.
Table 1. Stability of the boundary fixed points and the boundary 2-cycles of system (Equation1(1)
(1) ).
Figure 2. The parameter plane of system (Equation1
(1)
(1) ) with
. In regions A, B, and C,
,
, and
are asymptotically stable, respectively. There are no stable fixed points nor stable 2-cycles on the boundary of
if
.
![Figure 2. The parameter plane (λ,s) of system (Equation1(1) {x1(n+1)=x1(n){(1−η)exp(λ−x1(n)−x2(n))+ηs}x2(n+1)=x2(n)exp(λ−x1(n)−x2(n)),(1) ) with η≈0. In regions A, B, and C, (0,X2), P2, and P1 are asymptotically stable, respectively. There are no stable fixed points nor stable 2-cycles on the boundary of R+2 if λ>λ2.](/cms/asset/40126809-698d-4761-945b-81b021fb84ab/tjbd_a_2248171_f0002_ob.jpg)
Figure 3. The parameter plane of system (Equation1
(1)
(1) ). The green and red regions correspond to the region B and C in Figure , respectively. In the yellow region, both
and
have unstable transversal eigenvalues, i.e. both species can mutually invade at the boundary 2-cycles. On the left and right panels,
and
, respectively.
![Figure 3. The parameter plane (λ,s) of system (Equation1(1) {x1(n+1)=x1(n){(1−η)exp(λ−x1(n)−x2(n))+ηs}x2(n+1)=x2(n)exp(λ−x1(n)−x2(n)),(1) ). The green and red regions correspond to the region B and C in Figure 2, respectively. In the yellow region, both P1 and P2 have unstable transversal eigenvalues, i.e. both species can mutually invade at the boundary 2-cycles. On the left and right panels, η=0.01 and η=0.1, respectively.](/cms/asset/b8b611c3-9a46-40b8-9e4b-fd93fc932718/tjbd_a_2248171_f0003_oc.jpg)
From Table , we find that if belongs to the region C of Figure and
, then the monoculture system of species 2 is settled in
and cannot prevent invasion of species 1 while the monoculture system of species 1 is settled in
and can prevent invasion of species 2. This implies that delayed reproduction can be adaptive if a fraction η of juveniles delays the initiation of reproduction and both the juvenile survival probability and the amplitude of population fluctuation are sufficiently large. Although we do not have a mathematical estimate how small η should be, the numerical result shown in Figure suggests that the region C persists even if
. In Figure , the eigenvalues shown in (Equation10
(10)
(10) ) and (Equation11
(11)
(11) ) are numerically calculated. In the yellow region of Figure , two species can mutually invade each other at
and
and coexist. We see that the yellow region gets larger as η increases.
5. Global dynamics
In this section, we examine the global dynamics of system (Equation1(1)
(1) ) analytically and numerically. First, we show that it is unlikely that system (Equation1
(1)
(1) ) has a positive 2-cycle as long as
.
Theorem 5.1
Suppose that either or
and
. Then system (Equation1
(1)
(1) ) has no positive 2-cycles for all
sufficiently small.
Proof.
Let be a positive 2-cycle of system (Equation1
(1)
(1) ). Then it satisfies the system of equations
(13)
(13) where
(14)
(14) By using the second equation of (Equation13
(13)
(13) ), we can remove
from the first equation of (Equation13
(13)
(13) ) as follows:
(15)
(15) where
(16)
(16) Since
, the quadratic Equation (Equation15
(15)
(15) ) has the two positive distinct real roots
By using (Equation14
(14)
(14) ) and (Equation16
(16)
(16) ), we can remove
from the second equation of (Equation13
(13)
(13) ) as follows:
Thus,
and
must satisfy the following system of linear equations:
(17)
(17) Since the two equations of (Equation17
(17)
(17) ) represent parallel lines on the plane
if
and have negative slopes, (Equation17
(17)
(17) ) does not have positive solutions
for all
sufficiently small if
Since
for
, the above condition is equivalent to
where
is used to obtain the first equality.
Since is strictly decreasing for
, we have
, which implies that
holds for all
if
.
Since system (Equation1(1)
(1) ) has no positive fixed points, any orbit of system (Equation1
(1)
(1) ) cannot converge to a positive point. Furthermore, by the above theorem, we see that it is unlikely that an orbit staring at a positive point converges to a positive 2-cycle if
. Although system (Equation1
(1)
(1) ) might have more complicated positive attractors, numerical simulations suggest that any orbit of system (Equation1
(1)
(1) ) starting at a positive point cannot remain in a compact set bounded away from the coordinate axes as long as
and
belongs to the region B or C of Figure (see Figure ). Figure (a) shows the basin of attraction of
of system (Equation1
(1)
(1) ) whose
belongs to the region C of Figure and
. The points in the green and yellow regions are attracted by
. The two colours indicate whether a point is attracted by
or
when looking at the orbit every other time. It is unlikely that the points on the boundaries of regions with different colours are attracted by
. At least it is clear that the fixed point
has a stable manifold that intersects with the interior of
. Therefore, the points on the boundaries of two different colours seems to be attracted by
. Figure shows a single forward orbit of system (Equation1
(1)
(1) ) with the same parameters as in Figure (a). The orbit starts in a neighbourhood of
. Since the fixed point
is unstable, the orbit leaves its neighbourhood and approaches
after passing by the fixed point
. This behaviour suggests that there is a heteroclinic orbit from
to
. Figure (b) shows the basin of attraction of
of system (Equation1
(1)
(1) ) whose
belongs to the region B of Figure and
.
Figure 4. The basins of attraction of the boundary 2-cycles of system (Equation1(1)
(1) ). The vertical and horizontal axes are
and
, respectively. (a) The points in the green and yellow regions are attracted by
and
, respectively, under the second iterate of system (Equation1
(1)
(1) ). The parameters are
, s = 0.6, and
. (b) The points in the black and red regions are attracted by
and
, respectively, under the second iterate of system (Equation1
(1)
(1) ). The parameters are
, s = 0.1, and
.
![Figure 4. The basins of attraction of the boundary 2-cycles of system (Equation1(1) {x1(n+1)=x1(n){(1−η)exp(λ−x1(n)−x2(n))+ηs}x2(n+1)=x2(n)exp(λ−x1(n)−x2(n)),(1) ). The vertical and horizontal axes are x1 and x2, respectively. (a) The points in the green and yellow regions are attracted by (P1,0) and (P1∗,0), respectively, under the second iterate of system (Equation1(1) {x1(n+1)=x1(n){(1−η)exp(λ−x1(n)−x2(n))+ηs}x2(n+1)=x2(n)exp(λ−x1(n)−x2(n)),(1) ). The parameters are λ=2.5, s = 0.6, and η=0.1. (b) The points in the black and red regions are attracted by (0,P2) and (0,P2∗), respectively, under the second iterate of system (Equation1(1) {x1(n+1)=x1(n){(1−η)exp(λ−x1(n)−x2(n))+ηs}x2(n+1)=x2(n)exp(λ−x1(n)−x2(n)),(1) ). The parameters are λ=2.5, s = 0.1, and η=0.1.](/cms/asset/2e2cf4fd-9f1d-4a90-8cc1-7e19f014e353/tjbd_a_2248171_f0004_oc.jpg)
Figure 5. A single forward orbit of system (Equation1(1)
(1) ). The points of the orbit are connected with black and red lines if they are points of
th and 2n + 1th iteration, respectively. The parameters are the same as in Figure (a).
![Figure 5. A single forward orbit of system (Equation1(1) {x1(n+1)=x1(n){(1−η)exp(λ−x1(n)−x2(n))+ηs}x2(n+1)=x2(n)exp(λ−x1(n)−x2(n)),(1) ). The points of the orbit are connected with black and red lines if they are points of 2nth and 2n + 1th iteration, respectively. The parameters are the same as in Figure 4 (a).](/cms/asset/91158bb6-0e13-4f01-b196-3b60131bfeba/tjbd_a_2248171_f0005_oc.jpg)
6. Concluding remarks
We studied the dynamics of system (Equation1(1)
(1) ), which describes competition between two semelparous species. One species is assumed to have a fraction η of individuals who delay the initiation of reproduction while the other species is assumed to reproduce every unit time. All other properties of these two species are assumed to be identical. Under the assumption
, we analytically obtained a parameter plane
demarcated with the advantage of these two species (see Figure ). The parameter plane shows that delayed reproduction can be adaptive if both juvenile survival probability and amplitude of population fluctuation are sufficiently large as long as a fraction of juveniles delays reproduction.
System (Equation1(1)
(1) ) is a special case of system (Equation2
(2)
(2) ) studied by Yakubu [Citation19], who considers competition between two semelparous species and assumes that one species is endangered and planted (or stocked) with a constant rate. In this interpretation, system (Equation1
(1)
(1) ) considers the case where the endangered species whose fertility is reduced by the proportion η in compared with the competitor is planted (or stocked) at a rate
(0<s<1). Our result shows that if the population is fluctuating, then the endangered species can be saved by poor planting not covering the reduction of fertility. Note that a similar interpretation is possible if species 1 is considered iteroparous and species 2 semelparous (see [Citation1,Citation7]).
By applying the general result of Franke and Yakubu [Citation5, Theorem 6.1] to system (Equation1(1)
(1) ), we can show that if the fixed points
and
are global attractors within each positive coordinate axis, then
attracts all points in the interior of
since species 1 is dominated by species 2. However, as we showed, if the coordinate axes have the 2-cycles
and
, then both
and
can be asymptotically stable in system (Equation1
(1)
(1) ) with
. Note, however, that a stable boundary 2-cycle cannot attract all points in the interior of
since
always has a stable manifold that intersects with the interior of
. Our numerical simulation suggests that heteroclinic orbits from
to
belong to the stable manifold. Furthermore, as shown by Kang and Smith [Citation13], who examines a variant of system (Equation2
(2)
(2) ), it is likely that all the pre-images of the closure of the union of all the heteroclinic orbits also cannot be attracted by
(see also [Citation12]). It is a future problem to reveal the basin of attraction of
when it is asymptotically stable. It is worth noting that riddled basins of attraction are observed in a slightly more general system that includes system (Equation2
(2)
(2) ) as a special case [Citation2] and intermingled basins of attraction are constructed in a certain discrete-time competition model [Citation9].
Disclosure statement
No potential conflict of interest was reported by the author(s).
Additional information
Funding
References
- M. Bulmer, Theoretical Evolutionary Ecology, Sinauer Associates, Sunderland, Massachusetts, 1994.
- O. de Feo and R. Ferriere, Bifurcation analysis of population invasion: On-off intermittency and basin riddling, Internat. J. Bifur. Chaos Appl. Sci. Engrg. 10(2) (2000), pp. 443–452.
- S. Elaydi and R.J. Sacker, Basin of attraction of periodic orbits of maps on the real line, J. Differ. Equ. Appl. 10(10) (2004), pp. 881–888.
- S. Elaydi and A.-A. Yakubu, Global stability of cycles: Lotka-Volterra competition model with stocking, J. Differ. Equ. Appl. 8(6) (2002), pp. 537–549.
- J.E. Franke and A.-A. Yakubu, Global attractors in competitive systems, Nonlinear Anal. 16(2) (1991), pp. 111–129.
- J.E. Franke and A.-A. Yakubu, Mutual exclusion versus coexistence for discrete competitive systems, J. Math. Bio. 30(2) (1991), pp. 161–168.
- M. Gatto, The evolutionary optimality of oscillatory and chaotic dynamics in simple population models, Theor. Popul. Biol. 43(3) (1993), pp. 310–336.
- M. Hassell and H.N. Comins, Discrete time models for two-species competition, Theor. Popul. Biol. 9(2) (1976), pp. 202–221.
- F. Hofbauer, J. Hofbauer, P. Raith, and T. Steinberger, Intermingled basins in a two species system, J. Math. Bio. 49(3) (2004), pp. 293–309.
- J. Hofbauer, V. Hutson, and W. Jansen, Coexistence for systems governed by difference equations of lotka-volterra type, J. Math. Bio. 25(5) (1987), pp. 553–570.
- H. Jiang and T.D. Rogers, The discrete dynamics of symmetric competition in the plane, J. Math. Bio. 25(6) (1987), pp. 573–596.
- Y. Kang, Pre-images of invariant sets of a discrete-time two-species competition model, J. Differ. Equ. Appl. 18(10) (2011), pp. 1709–1733.
- Y. Kang and H. Smith, Global dynamics of a discrete two-species Lottery-Ricker competition model, J. Biol. Dyn. 6(2) (2012), pp. 358–376.
- R. Kon, Convex dominates concave: An exclusion principle in discrete-time Kolmogorov systems, Proc. Amer. Math. Soc. 134(10) (2006), pp. 3025–3034.
- R. Luís, S. Elaydi, and H. Oliveira, Stability of a Ricker-type competition model and the competitive exclusion principle, J. Biol. Dyn. 5(6) (2011), pp. 636–660.
- R.M. May and G.F. Oster, Bifurcations and dynamic complexity in simple ecological models, Am. Nat. 110(974) (1976), pp. 573–599.
- S.D. Mylius and O. Diekmann, The resident strikes back: Invader-induced switching of resident attractor, J. Theoret. Biol. 211(4) (2001), pp. 297–311.
- S. Tuljapurkar, Delayed reproduction and fitness in variable environments, Proc. Natl. Acad. Sci. USA 87(3) (1990), pp. 1139–1143.
- A.-A. Yakubu, The effects of planting and harvesting on endangered species in discrete competitive systems, Math. Biosc. 126(1) (1995), pp. 1–20.
- A.-A. Yakubu, A discrete competition system with planting, J. Differ. Equ. Appl. 4(2) (1998), pp. 213–214.