Discrete Dynamics in Nature and Society

Volume 2015 (2015), Article ID 847360, 16 pages

http://dx.doi.org/10.1155/2015/847360

## Basins of Attraction for Two-Species Competitive Model with Quadratic Terms and the Singular Allee Effect

Department of Mathematics, University of Rhode Island, Kingston, RI 02881-0816, USA

Received 26 August 2014; Accepted 29 October 2014

Academic Editor: Garyfalos Papashinopoulos

Copyright © 2015 A. Brett and M. R. S. Kulenović. 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

We consider the following system of difference equations: where , , , , are positive constants and are initial conditions. This system has interesting dynamics and it can have up to seven equilibrium points as well as a singular point at , which always possesses a basin of attraction. We characterize the basins of attractions of all equilibrium points as well as the singular point at and thus describe the global dynamics of this system. Since the singular point at always possesses a basin of attraction this system exhibits Allee’s effect.

#### 1. Introduction

The following difference equation is known as the Beverton-Holt model: where is the rate of change (growth or decay) and is the size of the population at the th generation.

This model was introduced by Beverton and Holt in 1957. It depicts density dependent recruitment of a population with limited resources which are not shared equally. The model assumes that the* per capita* number of offspring is inversely proportional to a linearly increasing function of the number of adults.

The Beverton-Holt model is well studied and understood and exhibits the following properties.(a)Equation (1) has two equilibrium points and when .(b)All solutions of (1) are monotonic (increasing or decreasing) sequences.(c)If , then the zero equilibrium is a global attractor; that is, , for all .(d)If , then the equilibrium point is a global attractor; that is, , for all .(e)Both equilibrium points are globally asymptotically stable in the corresponding regions of parameters and ; that is, they are global attractors with the property that small changes of initial condition result in small changes of the corresponding solution .

All these properties can be derived from the explicit form of the solution of (1): See [1–3].

The following difference equation, was introduced by Thomson [4] as a depensatory generalization of the Beverton-Holt stock-recruitment relationship used to develop a set of constraints designed to safeguard against overfishing; see [5] for further references. In view of the sigmoid shape of the function (3) is called the Sigmoid Beverton-Holt model. A very important feature of the Sigmoid Beverton-Holt model is that it exhibits the Allee effect; that is, zero equilibrium has a substantial basin of attraction, as we can see from the following results.(a)Equation (3) has a unique zero equilibrium when .(b)Equation (3) has a zero equilibrium and the positive equilibrium , when .(c) There exist a zero equilibrium and two positive equilibria, and , when .(d) All solutions of (3) are monotonic (increasing or decreasing) sequences.(e) If , then the equilibrium point is a global attractor; that is, .(f)If , then the equilibrium point is a global attractor, with the basin of attraction and is a nonhyperbolic equilibrium point with the basin of attraction .(g)If , then zero equilibrium and are locally asymptotically stable, while is repeller and the basins of attraction of the equilibrium points are given as In other words, the smaller positive equilibrium serves as the boundary between two basins of attraction. The zero equilibrium has the basin of attraction and the model exhibits the Allee effect.(h) The equilibrium points and are globally asymptotically stable in the corresponding basins of attractions and .

The two dimensional analogue of (1) is the uncoupled system where are positive parameters. The dynamics of system (5) can be derived from dynamics of each equation. Therefore, this system has an explicit solution given by (2).

Two species can interact in several different ways through competition, cooperation, or host-parasitoid interactions. For each of these interactions, we obtain variations of system (5) all of which may require different mathematical analysis.

One such variation that exhibits competitive interaction is the following model, known as the Leslie-Gower model, which was considered in Cushing et al. [6]: where all parameters are positive and the initial conditions are nonnegative. The global dynamics of system (6) was completed in [7]. Several variations of system (6) where the competition of two species was modeled by linear fractional difference equations were considered in [8–14]. An interesting fact is that none of these models exhibited the Allee effect.

The two dimensional analogue of system (3) is the following uncoupled system: where are positive parameters. The dynamics of system (7) can be derived from the dynamics of each equation in the system. Since each equation in system (7) has three possible dynamic scenarios, then system (7) possesses nine dynamic scenarios.

A variation of system (7) that exhibits competitive interactions is the system where . This system will be considered in the remainder of this paper. We will show that system (8) has similar but more complex dynamics than system (7). We will see that like system (7) the coupled system (8) may possess 1, 3, 5, or 7 equilibrium points in the hyperbolic case and 2, 4, or 6 equilibrium points in the nonhyperbolic case. In each of these cases we will show that the Allee effect is present, although is outside of the domain of definition of system (8). We will precisely describe the basins of attraction of all equilibrium points and the singular point . We will show that the boundaries of the basins of attraction of the equilibrium points are the global stable manifolds of the saddle or the nonhyperbolic equilibrium points. See [10, 11, 13–18] for related results and [19] for dynamics of competitive system with a singular point at the origin. The biological interpretation of a related system is given in [20, 21] and similar system is treated in [22]. The specific feature of our results is that no equilibrium point in the interior of the first quadrant is computable and so our analysis is based on geometric analysis of the equilibrium curves.

#### 2. Preliminaries

Our proofs use some recent general results for competitive systems of difference equations of the form: where and are continuous functions and is nondecreasing in and nonincreasing in and is nonincreasing in and nondecreasing in in some domain .

Competitive systems of the form (9) were studied by many authors in [6, 7, 9, 13, 14, 23–37] and others.

Here we give some basic notions about monotonic maps in the plane.

We define a* partial order * on (so-called South-East ordering) so that the positive cone is the fourth quadrant; that is, this partial order is defined by

Similarly, we define North-East ordering as

A map is called* competitive* if it is nondecreasing with respect to , that is, if the following holds:

For each , define for to be the usual four quadrants based on and numbered in a counterclockwise direction; for example, .

For let denote the* interior* of .

The following definition is from [35].

*Definition 1. *Let be a nonempty subset of . A competitive map is said to satisfy condition () if for every , in , implies , and is said to satisfy condition () if for every , in , implies .

The following theorem was proved by de Mottoni and Schiaffino [38] for the Poincaré map of a periodic competitive Lotka-Volterra system of differential equations. Smith generalized the proof to competitive and cooperative maps [34].

Theorem 2. *Let be a nonempty subset of . If is a competitive map for which () holds, then for all , is eventually componentwise monotone. If the orbit of has compact closure, then it converges to a fixed point of . If instead () holds, then for all , is eventually componentwise monotone. If the orbit of has compact closure in , then its omega limit set is either a period-two orbit or a fixed point.*

It is well known that a stable period-two orbit and a stable fixed point may coexist; see Hess [39].

The following result is from [35], with the domain of the map specialized to be the cartesian product of intervals of real numbers. It gives a sufficient condition for conditions () and ().

Theorem 3. *Let be the cartesian product of two intervals in . Let be a competitive map. If is injective and for all then satisfies (). If is injective and for all then satisfies ().*

Theorems 2 and 3 are quite applicable as we have shown in [40], in the case of competitive systems in the plane consisting of rational equations.

The following result is from [18], which generalizes the corresponding result for hyperbolic case from [7]. Related results have been obtained by Smith in [34].

Theorem 4. *Let be a rectangular subset of and let be a competitive map on . Let be a fixed point of such that has nonempty interior (i.e., is not the NW or SE vertex of ).**Suppose that the following statements are true. *(a)*The map is strongly competitive on .*(b)* is on a relative neighborhood of .*(c) *The Jacobian matrix of at has real eigenvalues , such that , where is stable and the eigenspace associated with is not a coordinate axis.*(d)*Either and
* *or and
**Then there exists a curve in such that *(i)* is invariant and a subset of ;*(ii) *the endpoints of lie on ;*(iii) *;*(iv) * is the graph of a strictly increasing continuous function of the first variable;*(v) * is differentiable at if or one sided differentiable if , and in all cases is tangential to at ;*(vi)* separates into two connected components, namely,
*(vii)* is invariant, and as for every ;*(viii)* is invariant, and as for every .*

*The following result is a direct consequence of the Trichotomy Theorem of Dancer and Hess (see [7, 39]) and is helpful for determining the basins of attraction of the equilibrium points.*

*Corollary 5. If the nonnegative cone of is a generalized quadrant in , and if has no fixed points in the ordered interval other than and , then the interior of is either a subset of the basin of attraction of or a subset of the basin of attraction of .*

*The next results give the existence and uniqueness of invariant curves emanating from a nonhyperbolic point of unstable type, that is, a nonhyperbolic point where second eigenvalue is outside interval . Similar result for a nonhyperbolic point of stable type, that is, a nonhyperbolic point where second eigenvalue is in the interval , follows from Theorem 4. See Kulenović and Merino, Invariant Curves of Planar Competitive and Cooperative Maps.*

*Theorem 6. Let and let be a strongly competitive map with a unique fixed point , such that is continuously differentiable in a neighborhood of . Assume further that at the point the map has associated characteristic values and satisfying and .Then there exist curves , in and there exist with such that (i)for , is invariant, north-east strongly linearly ordered, such that and ; the endpoints , of , where , belong to the boundary of . For with , is a subset of the closure of one of the components of . Both and are tangential at to the eigenspace associated with ;(ii)for , let be the component of whose closure contains . Then is invariant. Also, for , accumulates on , and for , accumulates on .(iii)Let and .Then is invariant.*

*Corollary 7. Let a map with fixed point be as in Theorem 6. Let , be the sets as in Theorem 6. If satisfies (), then for , is invariant, and for every , the iterates converge to or to a point of . If satisfies (), then and . For every , the iterates either converge to or converge to a period-two point or to a point of .*

*3. Local Stability of Equilibrium Points*

*3. Local Stability of Equilibrium Points*

*First we present the local stability analysis of the equilibrium points. It is interesting that the local stability analysis is the more difficult part of our analysis.*

*The equilibrium points of system (8) satisfy the following system of equations:
*

*All solutions of system (16) with at least one zero component are given as where , where , and where . The equilibrium point exists when , and exists when .*

*The equilibrium points with strictly positive coordinates satisfy the following system of equations:
*

*From (17) we have that all real solutions of the system (17) belong to the positive quadrant, since and . By eliminating from (17) we obtain
*

*The next result gives the necessary and sufficient conditions for (18) and so system (16) to have between zero and 4 solutions. As we show in Section 4.2 the global dynamics depends on the number of the equilibrium points with positive coordinates.*

*Lemma 8. Let
*

Assume that . Then the following holds. (a)If , , and , then (18) has four simple real roots.(b)If and , then (18) has no real roots.(c)If , then (18) has two simple real roots.(d)If and , then (18) has one real double root.(e)If and , then (18) has two real simple roots and one real double root.(f)If , , and , then (18) has two real double roots.(g)If , , and , then (18) has no real roots.(h)If , , and , then (18) has one real root of multiplicity four.

*Proof. *The discrimination matrix [41] of and is given by
Let denote the determinant of the submatrix of , formed by the first rows and the first columns, for where
So, by straightforward calculation one can see that
The rest of the proof follows in view of Theorem 1 in [41].

*Geometrically solutions of system (17) are intersections of two ellipses that satisfy the equations
with respective vertices and . See Figure 1.*