Hostname: page-component-586b7cd67f-t7fkt Total loading time: 0 Render date: 2024-11-22T15:08:44.966Z Has data issue: false hasContentIssue false

Simple dynamics in non-monotone Kolmogorov systems

Published online by Cambridge University Press:  02 December 2021

Lei Niu
Affiliation:
Department of Applied Mathematics, Donghua University, Shanghai 201620, China ([email protected]) Institute for Nonlinear Science, Donghua University, Shanghai 201620, China
Alfonso Ruiz-Herrera
Affiliation:
Department of Mathematics, Faculty of Science, University of Oviedo, Oviedo, Spain ([email protected])
Rights & Permissions [Opens in a new window]

Abstract

In this paper we analyse the global dynamical behaviour of some classical models in the plane. Informally speaking we prove that the folkloric criteria based on the relative positions of the nullclines for Lotka–Volterra systems are also valid in a wide class of discrete systems. The method of proof consists of dividing the plane into suitable positively invariant regions and applying the theory of translation arcs in a subtle manner. Our approach allows us to extend several results of the theory of monotone systems to nonmonotone systems. Applications in models with weak Allee effect, population models for pioneer-climax species, and predator–prey systems are given.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s), 2021. Published by Cambridge University Press on behalf of The Royal Society of Edinburgh

1. Introduction

The discrete equation

(1.1)\begin{equation} x_{n+1}=x_ng(x_n),\quad n\in \mathbb{Z}_+{=}\{0,1,2,\ldots\}, \end{equation}

is a popular modeling framework for analysing the dynamical behaviour of a single species. In (1.1), $x_n\geq 0$ is the population density in the $n$-th generation and $g(x_n)\geq 0$ represents the density-dependent growth rate (or fitness function) from generation to generation. A common assumption in population dynamics is that $g$ is decreasing. This means that the growth rate is mainly determined by negative density-dependent mechanisms such as intra-specific competition [Reference Fishman5] or cannibalism. However, cooperative predation, resource defense, increased availability of mates and conspecific enhancement of reproduction are other biological mechanisms producing non-monotone growth rates, see [Reference Kang and Yakubu20]. For example, alders, big leaf maples, poplars and some pine trees thrive at low densities and decrease at high ones due to overcrowding effects and other ecosystem constraints [Reference Franke and Yakubu6, Reference Franke and Yakubu8, Reference Franke and Yakubu10, Reference Selgrade and Namkoong33].

In describing the interactions of $k$ species, a natural extension of equation (1.1) is

(1.2)\begin{equation} x_i(n+1)=x_i(n)H_i(x_1(n),\ldots,x_k(n)),\quad i=1,\ldots,k. \end{equation}

The growth rate $H_i$ is typically of the form

(1.3)\begin{equation} H_i(x_1(n),\ldots,x_k(n))=g_i\left(\sum_{j=1}^{k}a_{ij}x_j(n)\right) \end{equation}

with $g_i:[0,\,+\infty )\longrightarrow [0,\,+\infty )$ a continuous function. In particular, if $a_{ij}>0$ and $g_{i}$ is decreasing for all $i,\,j=1,\,\ldots,\,n$, system (1.2) describes the evolution of $k$ competing species. See [Reference Franke and Yakubu7, Reference Franke and Yakubu9, Reference Franke and Yakubu10, Reference Gyllenberg, Jiang, Niu and Yan13, Reference May23] and the references therein for a detailed discussion on these models.

Understading the dynamical behaviour of (1.2) is of critical importance from an applied point of view. There are several approaches for analysing this issue for competitive systems. For example, the convexity arguments given by Kon [Reference Kon22], the split Liapunov function method given by Baigent and Hou [Reference Baigent and Hou2, Reference Hou and Baigent18] or the theory of carrying simplex [Reference Ruiz-Herrera31]. Recently, Hou [Reference Hou16, Reference Hou17] has provided a criterion of global attraction based on the relative position of the nullclines reminiscent to the classical results for the Lotka–Volterra system

(1.4)\begin{equation} \frac{dx_i}{dt}=x_i\left(r_i-\sum_{j=1}^{k}a_{ij}x_j\right),\quad i=1,\ldots,k, \end{equation}

see [Reference Ahmad and Lazer1] and the references therein. In contrast with the monotone case, the literature on non-monotone systems is relatively scarce. It is worth noting that system (1.2) can exhibit chaotic dynamics.

In this paper we describe the global dynamical picture of model (1.2) when the functions $g_{i}$ are not necessarily decreasing. In particular, our criteria could be perceived as an extension of Hou's results to non-monotone systems. The method of proof is completely different from those papers mentioned above. First, we divide the phase space in suitable positively invariant regions and then we apply the theory of translation arcs [Reference Campos, Ortega and Tineo3, Reference Graff and Ruiz-Herrera12, Reference Ortega27] in a subtle manner. As emphasized in § 3, our conclusions are rather sharp.

The organization of the paper is as follows. In § 2, we introduce some notation and definitions. In § 3, we give the main results of the paper. Applications to population models for pioneer-climax species [Reference Franke and Yakubu7, Reference Franke and Yakubu9, Reference Franke and Yakubu10] and species with weak Allee effects [Reference Kang and Yakubu20] are discussed. To show the versatility of our results with different interactions, we discuss the dynamical behaviour of several classical predator–prey systems. We conclude the paper with a discussion on our findings.

2. Mathematical framework

The Euclidean disk with centre at $p=(p_{1},\,p_{2})\in \mathbb {R}^{2}$ and radius $R>0$ is denoted by

\[ D(p,R)=\{(x,y)\in\mathbb{R}^{2}:\|(x,y)-(p_{1},p_{2})\|\leq R\} \]

with $\|(x,\,y)\|=\sqrt {x^{2}+y^{2}}$. A subset of the plane homeomorphic to a (closed) Euclidean disk is called a topological disk. A map $h:V\subset \mathbb {R}^{2}\longrightarrow V$ which is continuous and injective is called an embedding. Notice that if $h(V)=V$, then $h$ is a homeomorphism. In this section, we work, without further mention, with embeddings defined on topological disks. We recall that $h:D\longrightarrow D$ is an orientation preserving embedding if it has degree one, that is,

\[ \text{deg}(h-q_{0},U)=1 \]

where $h(p_{0})=q_{0}$ with $p_{0}\in Int D$ and $U\subset D$ is any open neighbourhood of $p_{0}$. The reader can consult the appendix in [Reference Ortega27] for a detailed discussion on topological degree and index. When $h: D\longrightarrow D$ is an embedding of class $\mathcal {C}^{1}$, the sign of the determinant of the Jacobian matrix of $h$ can determine whether $h$ preserves the orientation. Specifically, if

(2.1)\begin{equation} \det D h(x)> 0 \end{equation}

for all $x\in D$, then $h$ is an orientation preserving embedding. Now we give a practical criterion to deduce when a map is an orientation preserving embedding. We have employed the notation $\partial D$ to refer to the boundary of $D$ in $\mathbb {R}^{2}$.

Proposition 2.1 Let $h:D\longrightarrow D$ be a map of class $\mathcal {C}^{1}$ with $D\subset \mathbb {R}^{2}$ a topological disk. Assume the following conditions:

  1. (i) $h(\partial D)\subset D$.

  2. (ii) $\det Dh(x)>0$ for all $x\in D$.

  3. (iii) There is a point $q\in D$ with a unique pre-image, that is, $h^{-1}(\{q\})=\{p\}$.

Then, $h(D)\subset D$ and

\[ h:D\longrightarrow D \]

is an orientation preserving embedding.

Proof. Using that $h^{-1}(\{q\})=\{p\}$ and that $h$ is locally injective (by (ii)), we deduce that

\[ h:D\longrightarrow h(D) \]

is a homeomorphism (see lemma 2.3.4 in [Reference Chow and Hale4]). We know that $h(D)$ is a topological disk with boundary contained in the topological disk $D$. This implies that $h(D)\subset D$ (see lemma 6 in [Reference Ortega27, p. 39]). Hence,

\[ h: D\longrightarrow D \]

is an embedding. Notice that $h$ is an orientation preserving embedding because

\[ \det Dh(x)>0 \]

for all $x\in D$.

Given $h:D\longrightarrow D$ an embedding, the $\omega$-limit set of a point $p\in D$ is defined as

\[ \omega(h,p)=\{q\in D: h^{\sigma(n)}(p)\longrightarrow q\ for \ some\ sequence\ \sigma(n)\longrightarrow +\infty\}. \]

We say that $h$ has trivial dynamics if for all $x\in D$,

\[ \omega(h,x)\subset Fix(h) \]

where $Fix(h)$ denotes the fixed point set of $h$. It is worth noting that if $h$ is dissipative and $\omega (h,\,x)\subset Fix(h)$ for all $x\in D$, then $\omega (h,\,x)$ is a connected set (see proposition 2 in [Reference Ortega27, p. 42]). In particular, if $h$ has a finite number of fixed points, then for all $x\in D$, there is $q\in Fix(h)$ so that $\omega (h,\,x)=\{q\}$. We stress that an embedding $h:D\longrightarrow D$ defined on a topological disk is always dissipative.

Next we recall two results based on the theory of translation arcs.

Theorem 2.2 (Corollary 2.1 in [Reference Ruiz-Herrera30], [Reference Campos, Ortega and Tineo3])

Let $h: D\longrightarrow D$ be an orientation preserving embedding defined on the topological disk $D$. If $Fix(h)\subset \partial D,$ then $h$ has trivial dynamics.

Theorem 2.3 (Theorem 2.1 in [Reference Niu and Ruiz-Herrera26], [Reference Ortega27])

Let $h: D\longrightarrow D$ be an orientation preserving embedding defined on the topological disk $D$. If $Fix(h)\cap Int D=\{q\}$ with $index(h,\,q)=-1$ then $h$ has trivial dynamics.

In the previous theorem, $index(h,\,q)$ denotes the usual topological index of $h$ at $q$. We mention that theorem 2.1 in [Reference Niu and Ruiz-Herrera26] deals with homeomorphisms. However, the same proof works for embeddings.

3. Geometric criteria of trivial dynamics

In this section we consider the system

(3.1)\begin{equation} \left\{\begin{array}{@{}l} x_{n+1}=x_{n}g_{1}(x_{n}+\alpha y_{n})\\ y_{n+1}=y_{n}g_{2}(y_{n}+\beta x_{n}) \end{array}\right. \end{equation}

with $\alpha,\,\beta >0$ and $g_{i}:[0,\,+\infty )\longrightarrow (0,\,+\infty )$ a function of class $\mathcal {C}^{1}$ for $i=1,\,2$. Denote by

\[ G(x,y)=(G_{1}(x,y),G_{2}(x,y))=(xg_{1}(x+\alpha y),y g_{2}(y+\beta x)). \]

Notice that $G([0,\,+\infty )^{2})\subset [0,\,+\infty )^{2}$. A common assumption for the single-species model (1.1), already stated by Ricker [Reference Ricker28] and by Moran [Reference Moran25], is that there is a positive equilibrium $p$ (the carrying capacity) so that

  1. (C) $g(x)>1$ for all $x\in (0,\,p)$ and $g(x)<1$ for all $x>p$.

The previous condition can be adapted in model (3.1) as follows:

  1. (P) For $i=1,\,2$, there exists $r_{i}>0$ so that $g_{i}(x)>1$ if $x\in (0,\,r_{i})$ and $g_{i}(x)<1$ if $x>r_{i}$.

Biologically, condition (P) means that the density of population of species $i$ increases (resp. decreases) in the next generation when the weighted total biomass of both species is below (resp. above) a threshold.

In this section, we will make the assumption

  1. (S) $\det DG(x,\,y)>0$ for all $(x,\,y)\in \mathcal {R}$, where $\mathcal {R}=[0,\, A_{1}]\times [0,\,A_{2}]$ with

    \[ A_{1}=\max\left\{r_{1},\frac{r_{2}}{\beta}\right\},\quad A_{2}=\max\left\{r_{2},\frac{r_{1}}{\alpha}\right\}. \]

Condition (S) implies that $G$ is one-to-one on $\mathcal {R}$. The biological meaning of this is simple. If we take two different initial data, the densities of population are different each other in any future generation. To avoid technical problems, we also suppose:

  1. (H) $g_{i}'(r_{i})\not =0$ for $i=1,\,2$.

Define

\begin{align*} L_1& =\{(x,y)\in [0,+\infty)^{2}:x+\alpha y= r_1\},\\ L_2& =\{(x,y)\in [0,+\infty)^{2}:\beta x+ y= r_2\},\\ \mathcal{L}_1^{-}& =\{(x,y)\in [0,+\infty)^{2}:x+\alpha y\leq r_1\} \end{align*}

and

\[ \mathcal{L}_2^{-}=\{(x,y)\in [0,+\infty)^{2}:\beta x+ y\leq r_2\}. \]

The system has, at most, four equilibria namely $(0,\,0)$, $(r_{1},\,0)$, $(0,\,r_{2})$ and

(3.2)\begin{equation} (x^{*},y^{*})=\left(\frac{r_{1}-\alpha r_{2}}{1-\alpha\beta},\frac{r_{2}-\beta r_{1}}{1-\alpha \beta}\right). \end{equation}

Obviously, the last equilibrium is located at the intersection of $L_{1}$ and $L_{2}$. In our analysis, we exclude the case $L_{1}=L_{2}$ ($\alpha =\beta =1$ and $r_{1}=r_{2}$).

We prove that the relative position of $L_{1}$ and $L_{2}$ completely determines the dynamical behaviour of (3.1). This is a well-known result when the system is monotone. Our contribution will be to show that it is also true for non-monotone systems. We stress that if $(x,\,y)\in (0,\,+\infty )^{2}\backslash \mathcal {R}$, then

(3.3)\begin{equation} \left\{\begin{array}{@{}l} x>G_{1}(x,y)\\ y>G_{2}(x,y). \end{array}\right. \end{equation}

Proposition 3.1 Assume that (P) and (S) are satisfied. Then, $G(\mathcal {R})\subset \mathcal {R}$ and

\[ G|_{\mathcal{R}}: \mathcal{R}\longrightarrow \mathcal{R} \]

is an orientation preserving embedding.

Proof. First we notice that $G^{-1}(\{0\})=\{0\}$. Next we prove that $G(\partial \mathcal {R})\subset \mathcal {R}$. Take a point of the form $(A_{1},\,y)$ with $y\in [0,\, A_{2}]$. We observe that $A_{1}+\alpha y\geq r_{1}$ and $y+\beta A_{1}\geq r_{2}$. Then, by condition (P),

\begin{align*} & G_{1}(A_{1},y)=A_{1}g_{1}(A_{1}+\alpha y)\leq A_{1}\\ & G_{2}(A_{1},y)=yg_{2}(\beta A_{1}+ y)\leq y\leq A_{2}. \end{align*}

Now, it is clear that

\[ \{G(A_{1},y):y\in [0,A_{2}]\}\subset \mathcal{R}. \]

Analogously, we can prove that

\[ \{G(x,A_{2}):x\in [0,A_{1}]\}\subset \mathcal{R}. \]

We also observe that condition (S) ensures that $G_{1}(x,\,0)=xg_{1}(x)$ and $G_{2}(0,\,y)=yg_{2}(y)$ are strictly increasing (they are locally injective and have two fixed points). Thus, $G([0,\,A_{1}]\times \{0\})\subset [0,\,A_{1}]\times \{0\}$ and $G(\{0\}\times [0,\,A_{2}])\subset \{0\}\times [0,\,A_{2}]$. Collecting the above information, we deduce that

\[ G(\partial \mathcal{R})\subset\mathcal{R}. \]

The conclusion now follows from proposition 2.1 immediately.

Lemma 3.2 Let $\{(x_{n},\,y_{n})\}$ be a sequence of (3.1) with initial condition $(x_{0},\,y_{0})\not \in \mathcal {R}$. Assume that (P) and (S) hold. Then, one of the following properties is satisfied:

  1. (i) $\{(x_{n},\,y_{n})\}\longrightarrow q$ with $q$ an equilibrium of (3.1).

  2. (ii) There exists $n_{0}\in \mathbb {N}$ so that $(x_{n},\,y_{n})\in \mathcal {R}$ for all $n\geq n_{0}$.

Proof. Assume that $(x_{0},\,y_{0})\in (0,\,+\infty )^{2}$. If $(x_{n},\,y_{n})\not \in \mathcal {R}$ for all $n\in \mathbb {N}$, condition (P) ensures that $\{x_{n}\}$ and $\{y_{n}\}$ are strictly decreasing, see (3.3). Then, (i) holds. If there is $n_{0}\in \mathbb {N}$ so that $(x_{n_{0}},\,y_{n_{0}})\in \mathcal {R}$, we deduce that $(x_{n},\,y_{n})\in \mathcal {R}$ for all $n\geq n_{0}$. Notice that $G(\mathcal {R})\subset \mathcal {R}$ by the previous proposition. The case $x_{0}=0$ (resp. $y_{0}=0$) is treated in an analogous manner. Observe that in such a case $x_{n}=0$ (resp. $y_{n}=0$) for all $n\in \mathbb {N}$ and $y_{n}$ (resp. $x_{n}$) is decreasing provided $y_{n}>A_{2}$ (resp. $x_{n}>A_{1}$).

Remark 3.3 Lemma 3.2 says that it is enough for analysing the dynamical behaviour of (3.1) in $\mathcal {R}$. We repeatedly use this fact in the subsequent results.

The following result describes the behaviour of (3.1) on the axes.

Lemma 3.4 Assume that (P) and (S) hold. Then, $r_{1}$ and $r_{2}$ are globally asymptotically stable for

\[ x_{n+1}=x_{n}g_{1}(x_{n}) \]

and

\[ y_{n+1}=y_{n}g_{2}(y_{n}) \]

respectively in $(0,\,+\infty )$.

Proof. We study the equilibrium $r_{1}$. The other case is analogous and we omit the details. By remark 3.3, we can restrict our attention on $[0,\,A_{1}]$. As mentioned in the proof of proposition 3.1, $\varphi (x)=xg_{1}(x)$ is strictly increasing in $[0,\,A_{1}]$. We remark that the unique fixed points of $\varphi$ are $0$ and $r_{1}$, and the fact $\varphi (A_{1})\leq A_{1}$. By condition (P), $g_1(x)>1$ for all $x\in (0,\,r_{1})$. Hence, $\varphi (x)>x$ for all $x\in (0,\,r_{1})$. Now it is clear that $0$ is unstable. Thus, $r_{1}$ is a global attractor for $\varphi$ in $(0,\,A_{1}]$.

Now we are in a position to give the main result of this section.

Theorem 3.5 Assume (P), (H) and (S) hold. Then the dynamical behaviour of (3.1) is as follows:

  1. (i) $(r_{1},\,0)$ is a global attractor in $(0,\,+\infty )^{2}$ provided $\frac {r_{2}}{\beta }\leq r_{1}$ and $\frac {r_{1}}{\alpha }>r_{2}$.

  2. (ii) $(0,\,r_{2})$ is a global attractor in $(0,\,+\infty )^{2}$ provided $\frac {r_{2}}{\beta }> r_{1}$ and $\frac {r_{1}}{\alpha }\leq r_{2}$.

  3. (iii) $(x^{*},\,y^{*})$ is a global attractor in $(0,\,+\infty )^{2}$ provided $\frac {r_{2}}{\beta }>r_{1}$ and $\frac {r_{1}}{\alpha }>r_{2}$.

  4. (iv) There is trivial dynamics in (3.1) with $(r_{1},\,0)$ and $(0,\,r_{2})$ as local attractors in $(0,\,+\infty )^{2}$ provided $\frac {r_{2}}{\beta }< r_{1}$ and $\frac {r_{1}}{\alpha }< r_{2}$.

Proof. First we realize that the origin is always a local repeller in $(0,\,+\infty )^{2}$. Indeed, by condition (P), there is a neighbourhood $U$ of the origin in $\mathcal {R}$ so that $g_{1}(x+\alpha y)>1$ and $g_{2}(\beta x+y)>1$ for all $(x,\,y)\in U$ with $x\not =0$ and $y\not =0$. This implies that $G_{1}(x,\,y)>x$ and $G_{2}(x,\,y)>y$ for all $(x,\,y)\in U$ with $x\not =0$ and $y\not =0$. Since $G((0,\,+\infty )^{2})\subset (0,\,+\infty )^{2}$, it is clear that the origin is a local repeller in $(0,\,+\infty )^{2}$. We also stress that by remark 3.3, it is enough to study the dynamical behaviour in $\mathcal {R}$. Now we are ready to prove the theorem.

(i) In this case, $L_{1}$ and $L_{2}$ do not intersect in $(0,\,+\infty )^{2}$. Then, $G|_{\mathcal {R}}$ has three fixed points, namely, $(0,\,0)$, $(r_{1},\,0)$ and $(0,\,r_{2})$. Since $Fix(G|_{\mathcal {R}})\subset \partial \mathcal {R}$, theorem 2.2 and proposition 3.1 imply that for each $(x_{0},\,y_{0})\in Int\mathcal {R}$, there exists $q\in Fix(G|_{\mathcal {R}})$ so that

\[ \omega(G|_{\mathcal{R}}, (x_{0},y_{0}))=\{q\}. \]

Next we prove that $q\not =(0,\,0)$ and $q\not =(0,\,r_{2})$. To see this, we check that both fixed points are local repellers in $(0,\,+\infty )^{2}$. At the beginning, we have already mentioned this fact for the origin. On the other hand, the eigenvalues of the linearized system at $(0,\,r_{2})$ are

(3.4)\begin{equation} \{1+r_{2}g_{2}'(r_{2}),g_{1}(\alpha r_{2})\} \end{equation}

where the associated eigenvectors are $(0,\,1)$ and $(w_{1},\,w_{2})$ with $w_{1}\not =0$ respectively. Since $r_{1}>\alpha r_{2}$ and (P), we conclude that $g_{1}(\alpha r_{2})>1$. On the other hand, we have already seen in the proof of proposition 3.1 that $\phi (y)=yg_{2}(y)$ is an increasing function. Moreover, $r_{2}>0$ is a global attractor of $\phi$ by lemma 3.4. Using these two facts together with (H), we conclude that $\phi '(r_{2})=1+r_{2}g_{2}'(r_{2})\in [0,\,1)$. Observe that by (S), $\phi '(r_{2})\not =0$. Now, it is clear that $(0,\,r_{2})$ is a hyperbolic saddle point. In particular, it is a local repeller in $(0,\,+\infty )^{2}$.

The proof of (ii) is analogous and we omit the details.

(iii) In this case, $G|_{\mathcal {R}}$ has four fixed points, namely, $(0,\,0)$, $(r_{1},\,0)$, $(0,\,r_{2})$ and $(x^{*},\,y^{*})$. Define $\mathcal {A}=\mathcal {L}_{1}^{-}\cup \mathcal {L}_{2}^{-}$. The boundary of $\mathcal {A}$ is made of four segments: $\{0\}\times [0,\,\frac {r_{1}}{\alpha }]$, $[0,\,\frac {r_{2}}{\beta }]\times \{0\}$, $l_{1}$ and $l_{2}$ where $l_{1}$ is the segment on $L_{1}$ with ends at $(0,\,\frac {r_{1}}{\alpha })$ and $(x^{*},\,y^{*})$ and $l_{2}$ is the segment on $L_{2}$ with ends at $(\frac {r_{2}}{\beta },\,0)$ and $(x^{*},\,y^{*})$. Since $\frac {r_{2}}{\beta }>r_{1}$ and $\frac {r_{1}}{\alpha }>r_{2}$, we have that

(3.5)\begin{equation} l_{1}\backslash\{(x^{*},y^{*})\}\subset [0,+\infty)^{2}\backslash \mathcal{L}_{2}^{-} \end{equation}

and

(3.6)\begin{equation} l_{2}\backslash\{(x^{*},y^{*})\}\subset [0,+\infty)^{2}\backslash \mathcal{L}_{1}^{-}. \end{equation}

Then, given $(x,\,y)\in l_{1}$ different from the ends $(x^{*},\,y^{*})$ and $(0,\,\frac {r_{1}}{\alpha })$, we have that $G_{1}(x,\,y)=x$ because $(x,\,y)\in l_{1}\subset L_{1}$ and $G_{2}(x,\,y)< y$ by (3.5). Analogously,

\begin{align*} G_{1}(x,y)< x \\ G_{2}(x,y)=y \end{align*}

for all $(x,\,y)\in l_{2}\backslash \{(x^{*},\,y^{*}),\,(\frac {r_{2}}{\beta },\,0)\}$. These expressions guarantee that

\[ G|_{\mathcal{R}}(\partial \mathcal{A})\subset \mathcal{A}. \]

Recall that the intervals $[0,\,A_{1}]\times \{0\}$ and $\{0\}\times [0,\,A_{2}]$ are positively invariant under $G$. Now, arguing as in the proof of proposition 2.1, we have that $G(\mathcal {A})\subset \mathcal {A}$ and

\[ G|_{\mathcal{A}}: \mathcal{A}\longrightarrow \mathcal{A} \]

is an orientation preserving embedding. A critical fact is that

\[ Fix(G|_{\mathcal{A}})\subset \partial \mathcal{A}. \]

Thus, if we apply theorem 2.2 to $G|_{\mathcal {A}}$, we conclude that for each $p\in \mathcal {A}$, there exists $q\in Fix(G|_{\mathcal {A}})$ so that

\[ \omega(G|_{\mathcal{A}},p)=\{q\}. \]

On the other hand, we notice that $\frac {r_{2}}{\beta }>r_{1},\, \frac {r_{1}}{\alpha }>r_{2}$ and (P) imply that

\[ g_{2}(\beta r_{1})>1\text{ and } g_{1}(\alpha r_{2})>1. \]

Repeating the argument of the proof of (i), we can prove $(r_{1},\,0)$ and $(0,\,r_{2})$ are local saddle points. Recall that the origin is always a local repeller in $(0,\,+\infty )^{2}$. Consequently, for all $p\in \mathcal {A}\cap (0,\,+\infty )^{2}$,

\[ \omega(G|_{\mathcal{A}},p)=\{(x^{*},y^{*})\}. \]

Finally, we consider a sequence $\{(x_{n},\,y_{n})\}$ obtained from (3.1) so that $(x_{0},\,y_{0})\in \mathcal {R}\backslash \mathcal {A}$ with $x_{0}>0$ and $y_{0}>0$. By the same argument as that in lemma 3.2, one of the following facts holds:

  1. h1 $(x_{n},\,y_{n})$ tends to a fixed point of $G$.

  2. h2 $(x_{n},\,y_{n})\in \mathcal {A}$ for all $n\geq n_{0}$ with $n_{0}$ large enough.

If h1 holds, then the fixed point has to be $(x^{*},\,y^{*})$ because the other fixed points are local repellers in $(0,\,+\infty )^{2}$. If h2 holds, then we apply the above argument. The proof of (iii) is now completed.

(iv) In this case, $G|_{\mathcal {R}}$ has four fixed points, namely, $(0,\,0)$, $(r_{1},\,0)$, $(0,\,r_{2})$ and $(x^{*},\,y^{*})$. By theorem 5.1 in [Reference Ruiz-Herrera32], we deduce that

\[ index (G,(x^{*},y^{*}))={-}1. \]

Obviously, $index (G|_{\mathcal {R}},\,(x^{*},\,y^{*}))=index (G,\,(x^{*},\,y^{*}))$. Then, applying theorem 2.3 to $G|_{\mathcal {R}}$, we conclude that $G|_{\mathcal {R}}$ has trivial dynamics. Finally, linearizing the system and arguing as in (i), we conclude that $(r_{1},\,0)$ and $(0,\,r_{2})$ are local attractors.

Example 3.6 Models with weak Allee effect.

The predation by a generalist predator with a saturating functional response is a common mechanism associated with the presence of weak Allee effects [Reference Kang and Yakubu20]. A natural extension in the plane of the single species models with this Allee effect is

(3.7)\begin{equation} \left\{\begin{array}{@{}l} x_{n+1}= x_{n} g(x_{n}+\alpha y_{n})\\ y_{n+1}=y_{n} g(y_{n}+\beta x_{n}), \end{array}\right. \end{equation}

where $\alpha,\,\beta >0,$

\[ g(x)= \exp\left(r(1-x)-\frac{a}{1+b x}\right) \]

with $r>a>0$ and $b>0$, see lemma 1.1 in [Reference Kang and Yakubu20].

Although $g$ is not always monotone (see Fig. 1), we can apply theorem 3.5. Some particular choices of parameters are $r=0.5$, $a=0.45$ and $b=3$ together with $\alpha =0.5$ and $\beta =0.6$ (i); $\alpha =1.1$ and $\beta =0.6$ (iii); $\alpha =1.1$ and $\beta =1.4$ (iv). In general, condition (S) is difficult to verify in (3.7) because $\det D G(x,\,y)$ has a very complex expression. Notice that $det DG(x,\,y)$ can be negative for some points $(x,\,y)\in \mathcal {R}$ and for some choice of the parameters. Actually, system (3.7) can exhibit chaotic dynamics.

FIGURE 1. Representation of the function $g$ in system (3.7) with $r=0.5$, $a=0.45$ and $b=3$.

Example 3.7 A nonmonotone version of May–Oster model [Reference May and Oster24].

Consider the system

(3.8)\begin{equation} \left\{\begin{array}{@{}l} x_{n+1}= x_{n} g(x_{n}+\alpha y_{n})\\ y_{n+1}=y_{n} g(y_{n}+\beta x_{n}), \end{array}\right. \end{equation}

with $g(x)=\exp (x(r-x))$. After tedious computation, we can check that

\[ \det DG(x,y)>0 \]

for all $(x,\,y)\in \mathcal {R}$ in the following range of parameters:

  1. Case 1: $r=0.1$, $\alpha \in (0.5,\,1)$ and $\beta \in (3,\,4.2)$.

  2. Case 2: $r=0.1$, $\alpha \in (0.5,\,1)$ and $\beta \in (0.2,\,0.5)$.

  3. Case 3: $r=0.1$, $\alpha \in (1,\,6)$ and $\beta \in (1,\,6)$.

As a direct application of theorem 3.5, we obtain that $(0.1,\,0)$ is a global attractor in case 1, there is an interior fixed point $(x^{*},\,y^{*})$ that is a global attractor in case 2 and there is trivial dynamics with $(0.1,\,0)$ and $(0,\,0.1)$ as local attractors in case 3. As mentioned in the title of the example, (3.8) is a variant of the classical model discussed in [Reference May and Oster24]. We mention that this type of growth rates also appear in the evolution of climax species, see [Reference Franke and Yakubu7, Reference Franke and Yakubu9, Reference Franke and Yakubu10].

Example 3.8 The conclusions in theorem 3.5 might not be true when the condition (S) does not hold.

Consider

(3.9)\begin{equation} \left\{\begin{array}{@{}l} x_{n+1}= x_{n} g_1(x_{n}+\alpha y_{n})\\ y_{n+1}=y_{n} g_2(y_{n}+\beta x_{n}) \end{array}\right. \end{equation}

with $g_1(x)=\frac {r}{x+a}$ and $g_2(x)=\exp (p-qx)$, where $\alpha =\beta =a=1$, $r=21$, $p=2.5$, $q=0.1$. In this case, we have $\mathcal {R}=[0,\,25]\times [0,\,25]$. One can check that

\[ \det DG(0,20)={-}\sqrt{e}<0, \]

so condition (S) does not hold for model (3.9) because $(0,\,20)\in \mathcal {R}$. As shown in [Reference Franke and Yakubu7], there is an asymptotically stable 2-cycle for this model so that the two species can coexist. This indicates that the conclusions in theorem 3.5 is not true for the model (3.9) by noticing that $\frac {r_{2}}{\beta }> r_{1}$ and $\frac {r_{1}}{\alpha }\leq r_{2}$.

4. Extinction in planar systems beyond (3.1)

Predator–prey models are prototypes of Kolmogorov systems in which the results of § 3 are not directly applicable. However, some ideas developed in theorem 3.5 also work in these models. This shows the versatility of the mathematical framework given in § 2.

Consider the system

(4.1)\begin{equation} \left\{\begin{array}{@{}l} x_{n+1}=x_{n}f_{1}(x_{n},y_{n})\\ y_{n+1}=y_{n}f_{2}(x_{n},y_{n}) \end{array}\right. \end{equation}

with $f_{1},\, f_{2}:[0,\,+\infty )^{2}\longrightarrow (0,\,+\infty )$ functions of class $\mathcal {C}^{1}$. We denote by

\[ F(x,y)=(F_{1}(x,y),F_{2}(x,y))=(xf_{1}(x,y),y f_{2}(x,y)) \]

the map associated with (4.1). To model the predator–prey iteraction, we assume the following conditions ($x$ predator and $y$ prey):

  1. (PP1) $\frac {\partial F_{1}}{\partial x}(x,\,y)>0$ and $\frac {\partial F_{2}}{\partial y}(x,\,y)>0$ for all $(x,\,y)\in [0,\,+\infty )^{2}$.

  2. (PP2) $\frac {\partial f_{1}(x,\,y)}{\partial y}>0$ and $\frac {\partial f_{1}(x,\,y)}{\partial x}< 0$ for all $(x,\,y)\in [0,\,+\infty )^{2}$.

  3. (PP3) $\frac {\partial f_{2}(x,\,y)}{\partial x}<0$ and $\frac {\partial f_{2}(x,\,y)}{\partial y}< 0$ for all $(x,\,y)\in [0,\,+\infty )^{2}$.

  4. (LG) There are two fixed points $(p^{*},\,0)$ and $(0,\,q^{*})$ with $p^{*},\,q^{*}>0$ which are globally asymptotically stable on the $x$-axis and $y$-axis respectively.

  5. (PP4) There is a constant $\widetilde {K}>0$ so that $x f_{1}(x,\,q^{*}+1)< x$ for all $x>\widetilde {K}$.

Condition (PP1) means that the intra-specific competition is contest (see [Reference Hassell15]). (PP2) indicates that the growth rate of the predator is the result of the conjunction of two biological facts: the intra-specific competition and the consumption of the prey. (PP3) has an analogous meaning for the prey. (LG) says that each species in isolation has logistic growth rate with carrying capacities $p^{*}$ and $q^{*}$ respectively. (PP4) says that the predator density decreases in the next generation when it is above a suitable threshold.

We say that system (4.1) is permanent if there are two constants $\varepsilon,\, M>0$ so that

\begin{align*} \varepsilon & \leq \liminf x_{n}\leq \limsup x_{n}\leq M\\ \varepsilon & \leq \liminf y_{n}\leq \limsup y_{n}\leq M \end{align*}

for all sequence $\{(x_{n},\,y_{n})\}$ of (4.1) with initial condition $(x_{0},\,y_{0})\in (0,\,+\infty )^{2}$. Informally speaking, the notion of permanence excludes the extinction of some species in the system. To apply the classical theory of permanent systems we need to recall the notion of absorbing set. We say that a positively invariant set $\mathcal {R}\subset [0,\,+\infty )^{2}$ is an absorbing set for (4.1) if for all $(x_{0},\,y_{0})\in [0,\,+\infty )^{2}$, there is $n_{0}\in \mathbb {N}$ so that $F^{n}(x_{0},\,y_{0})\in \mathcal {R}$ for all $n\geq n_{0}$.

The following lemma is an immediate consequence of lemma 2.1 in [Reference Hutson and Moran19] by considering the average Liapunov function $V(x,\,y)=x^{\mu _1}y^{\mu _2}$. See also theorem 3.1 in [Reference Ruiz-Herrera29].

Lemma 4.1 Assume that the map $F$ associated with (4.1) has an absorbing set $\mathcal {R}:=[0,\,K_1]\times [0,\,K_2]$ with $K_1,\,K_2>0$. If there are real numbers $\mu _1,\,\mu _2>0$ such that

(4.2)\begin{equation} \mu_1\ln f_1(q)+\mu_2\ln f_2(q)>0 \end{equation}

for each $q$ in $\Omega (\partial \mathcal {R}):=\displaystyle \bigcup _{p\in \partial \mathcal {R}}\omega (F,\,p),$ then (4.1) is permanent.

Next we give the main result of this section.

Theorem 4.2 Suppose that (PP1)–(PP4) and (LG) hold. Then every sequence of (4.1) is bounded. Moreover, we have the following conclusion:

  1. (i) If $f_{2}(p^{*},\,0)>1,$ (4.1) is permanent.

  2. (ii) If $f_{2}(p^{*},\,0)\leq 1,$ $(p^{*},\,0)$ is a global attractor in $(0,\,+\infty )^{2}$.

Proof. First, let us prove that $\mathcal {R}=[0,\,\widetilde {K}]\times [0,\,q^{*}+1]$ is an absorbing set for (4.1). Take $\{(x_{n},\,y_{n})\}$ a sequence of (4.1). Define $\phi (y)=F_{2}(0,\,y)$. It follows from (PP1) and (PP3) that $\phi (y)$ is strictly increasing and $\phi (y)\geq F_{2}(x,\,y)$ for all $(x,\,y)\in [0,\,+\infty )^{2}$. From these properties we can deduce that

(4.3)\begin{equation} y_{n}\leq \phi^{n}(y_{0}) \end{equation}

for all $n\in \mathbb {N}$. Indeed,

\[ y_{1}=F_{2}(x_{0},y_{0})\leq \phi(y_{0}). \]

Repeating this argument, we obtain that

\[ y_{2}\leq \phi(y_{1}). \]

Using that $\phi$ is increasing and the previous step,

\[ y_{2}\leq \phi^{2}(y_{0}). \]

We complete the argument after a simple induction. Since $q^{*}>0$ is a global attractor of

\[ z_{n+1}=\phi(z_{n}) \]

in $(0,\,+\infty )$ by (LG) it is clear that there exists $n_{1}\in \mathbb {N}$ so that $y_{n}\in [0,\,q^{*}+1]$ for all $n\geq n_{1}$ by (4.3). Next we define $\varphi (x)=F_{1}(x,\,q^{*}+1)$. By (PP1), this function is strictly increasing. Using (PP4), we have that for all $z_{0}\in [0,\,+\infty )$, the sequence $\{z_{n}\}$ obtained from

\[ z_{n+1}=F_{1}(z_{n},q^{*}+1) \]

satisfies that $z_{n}\longrightarrow q$ with $q<\widetilde {K}$. We also know by (PP2) that $F_{1}(x_{n},\,y_{n})\leq F_{1}(x_{n},\,q^{*}+1)=\varphi (x_{n})$ for all $n\geq n_{1}$. Arguing as above, we can find $n_{2}\geq n_{1}$ so that

\[ x_{n}\in [0,\widetilde{K}] \]

for all $n\geq n_{2}$. Now, it is clear that $[0,\,\widetilde {K}]\times [0,\,q^{*}+1]$ is an absorbing set for (4.1). Next we properly prove the theorem.

(i) It follows from (LG) that

\[ \Omega(\partial \mathcal{R})=\{(0,0),(p^{*},0),(0,q^{*})\}. \]

Note that $f_1(p^{*},\,0)=1$ and $f_2(0,\,q^{*})=1$ because $(p^{*},\,0)$ and $(0,\,q^{*})$ are fixed points. Then, (PP2) and (PP3) imply that $f_i(0,\,0)>1$, $i=1,\,2$, and $f_1(0,\,q^{*})>1$. We know by assumptions that $f_{2}(p^{*},\,0)>1$. At this moment, it is clear there are $\mu _1,\,\mu _2>0$ such that (4.2) holds for each $q\in \{(0,\,0),\,(p^{*},\,0),\,(0,\,q^{*})\}$. Hence (4.1) is permanent by lemma 4.1. (ii) It suffices to prove that every sequence of (4.1) with initial condition in $\mathcal {R}$ converges to $(p^{*},\,0)$. First we prove that (4.1) does not admit fixed points in $(0,\,+\infty )^{2}$. Assume, by contradiction, that $(p,\,q)\in (0,\,+\infty )^{2}$ is a fixed point of $F$. Then,

(4.4)\begin{equation} \left\{\begin{array}{@{}l} 1=f_{1}(p,q),\\ 1=f_{2}(p,q). \end{array}\right. \end{equation}

By (PP2), we obtain that

\[ f_{1}(p^{*},0)=f_{1}(p,q)=1>f_{1}(p,0). \]

Consequently, $p>p^{*}$. On the other hand, by (PP3), we deduce that

\[ 1=f_{2}(p,q)< f_{2}(p,0)< f_{2}(p^{*},0)\leq 1, \]

a contradiction. After simple computations, we can prove that

\[ \det DF(x,y)>0 \]

for all $(x,\,y)\in [0,\,+\infty )^{2}$. Moreover, it is clear that $F^{-1}(\{0\})=\{0\}$. Then, by proposition 2.1, we have that

\[ F|_{\mathcal{R}}:\mathcal{R}\longrightarrow \mathcal{R} \]

is an orientation preserving embedding. We notice that

\[ Fix(F|_{\mathcal{R}})\subset \partial \mathcal{R}. \]

Moreover, $f_{1}(0,\,0)>1$, $f_{2}(0,\,0)>1$ and $f_{1}(0,\,q^{*})>1$, (see the proof of (i)). Linearizing the system, we easily obtain that $(0,\,0)$ is a local repeller. Regarding $(0,\,q^{*})$, we have that $\{f_{1}(0,\,q^{*}),\,1+q^{*} \frac {\partial f_{2}}{\partial y}(0,\,q^{*})\}$ are the eigenvalues of the linearized system. By (PP1) and (PP3), $0<\frac {\partial F_{2}}{\partial y}(0,\,q^{*})=1+q^{*} \frac {\partial f_{2}}{\partial y}(0,\,q^{*})<1$. In other words, $(0,\,q^{*})$ is a saddle point so that the repelling direction is $(w_{1},\,w_{2})$ with $w_{1}>0$. The conclusion follows from theorem 2.2.

There are many classical models that satisfy conditions (PP1)(PP4) and (LG). For example,

(4.5)\begin{equation} \left\{\begin{array}{@{}l} x_{n+1}= \dfrac{ r x_{n}}{1+x_{n}+h_{1}(y_{n})}\\ y_{n+1}=\dfrac{s y_{n}}{1+ y_{n}+h_{2}(x_{n})} \end{array}\right. \end{equation}

where $r>1+h_{1}(0)$, $s>1+h_{2}(0)$ and the functions $h_{1},\,h_{2}:[0,\,+\infty )\longrightarrow (0,\,+\infty )$ are of class $\mathcal {C}^{1}$ with $h_{1}'(x)<0$ and $h_{2}'(x)>0$ for all $x\in (0,\,+\infty )$.

5. Discussion

It is well known that the relative position of the nullclines determines the dynamical behaviour of the classical Lotka–Volterra model

(5.1)\begin{equation} \left\{\begin{array}{@{}l} x'=x (r_{1}-x-\alpha y)\\ y'=y(r_{2}-y-\beta x). \end{array}\right. \end{equation}

In this paper we have proved that the same results remain true in a broad family of discrete systems, namely

(5.2)\begin{equation} \left\{\begin{array}{@{}l} x_{n+1}=x_{n} g_{1}(x_{n}+\alpha y_{n})\\ y_{n+1}=y_{n} g_{2}(y_{n}+\beta x_{n}). \end{array}\right. \end{equation}

In our analysis, we have imposed mainly the following:

  1. (P) For $i=1,\,2$, there exists $r_{i}>0$ so that $g_{i}(x)>1$ if $x\in (0,\,r_{i})$ and $g_{i}(x)<1$ if $x>r_{i}$.

  2. (S) $\det DG(x,\,y)>0$ for all $(x,\,y)\in \mathcal {R}$ with $\mathcal {R}=[0,\, A_{1}]\times [0,\,A_{2}]$ with

    \[ A_{1}=\max\left\{r_{1},\frac{r_{2}}{\beta}\right\},\quad A_{2}=\max\left\{r_{2},\frac{r_{1}}{\alpha}\right\} \]
    and $G(x,\,y)=(x g_{1}(x+\alpha y),\,y g_{2}(y+\beta x))$.

The role of these conditions is critical. Notice that (P) is a necessary condition to maintain the dynamical behaviour of (5.1) in (5.2). On the other hand, if we drop (S), as discussed in § 3, new phenomena emerge in (5.2) in comparison with (5.1). For example, the presence of 2-cycles or chaotic dynamics. It is worth stressing that condition (P) encompasses functions that are not monotone. Particularly, we can describe the dynamical behaviour of species subject to Allee effects. Other marked examples are the population models for pioneer-climax species that appears when $g_{1}$ is decreasing and $g_{2}$ is one-humped, see [Reference Harper14, Reference Kim and Marlin21, Reference Selgrade and Namkoong33, Reference Sumner34]. For a long time, the dominant topic in these models was mainly the exclusion of the pioneer species, ignoring other dynamical patterns. In contrast with this point of view, theorem 3.5 describes all the possible dynamical patterns. In particular, we study the exclusion for the climax species, which has also been analysed recently by Gilbertson and Kot [Reference Gilbertson and Kot11].

There are many models in population dynamics that are not of the form (5.2), i.e., the growth rates are not a scalar function composed by a linear combination of the densities of the species. Nevertheless, many results are valid when the map associated with the model satisfies condition (S). This is the case of most predator–prey systems with a generalist predator. For these models, we have proved that the absence of coexistence states leads to the exclusion of the prey. A direct consequence of this is that the presence of any oscillatory behaviour in $(0,\,+\infty )^{2}$ implies the existence of a coexistence state.

Acknowledgments

The authors would like to thank the editor and the anonymous referees for their feedback and comments, which improved the quality of the manuscript. The work of L. Niu was supported by the National Natural Science Foundation of China (NSFC) under Grant No. 12001096 and the Fundamental Research Funds for the Central Universities (No. 2232021G-13).

References

Ahmad, S. and Lazer, A. C., One species extinction in an autonomous competition model, In Proc. First World Congress Nonlinear Analysis (Walter de Gruyter, Berlin, 1995).CrossRefGoogle Scholar
Baigent, S. and Hou, Z.. Global stability of discrete-time competitive population models. J. Differ. Equ. Appl. 23 (2017), 13781396.CrossRefGoogle Scholar
Campos, J., Ortega, R. and Tineo, A.. Homeomorphisms of the disk with trivial dynamics and extinction of competitive systems. J. Differ. Eq. 138 (1997), 157170.CrossRefGoogle Scholar
Chow, S. N. and Hale, J. K.. Methods of Bifurcation Theory (Springer-Verlag, New York, 1982).CrossRefGoogle Scholar
Fishman, M. A.. Density effects in population growth: an exploration. Biosystems 40 (1997), 219236.CrossRefGoogle ScholarPubMed
Franke, J. E. and Yakubu, A.. Global attractors in competitive systems. Nonlinear Anal.: Theory, Methods Appl. 16 (1991a), 111129.CrossRefGoogle Scholar
Franke, J. E. and Yakubu, A.. Mutual exclusion versus coexistence for discrete competitive systems. J. Math. Biol. 30 (1991b), 161168.CrossRefGoogle Scholar
Franke, J. E. and Yakubu, A.. Geometry of exclusion principles in discrete systems. J. Math. Anal. Appl. 168 (1992), 385400.CrossRefGoogle Scholar
Franke, J. E. and Yakubu, A.. Exclusion principles for density-dependent discrete pioneer-climax models. J. Math. Anal. Appl. 187 (1994b), 10191046.CrossRefGoogle Scholar
Franke, J. E. and Yakubu, A.. Pioneer exclusion in a one-hump discrete pioneer-climax competitive system. J. Math. Biol. 32 (1994a), 771787.CrossRefGoogle Scholar
Gilbertson, N. and Kot, M., Dynamics of a discrete-time pioneer–climax model, Theoretical Ecology .Google Scholar
Graff, G. and Ruiz-Herrera, A.. A strategy to locate fixed points and global perturbations of ODE's: mixing topology with metric conditions. J. Dynam. Differ. Eq. 26 (2014), 93107.CrossRefGoogle Scholar
Gyllenberg, M., Jiang, J., Niu, L. and Yan, P.. Permanence and universal classification for discrete-time competitive systems via the carrying simplex. Discrete Contin. Dyn. Syst. 40 (2020), 16211663.CrossRefGoogle Scholar
Harper, J. L.. Population Biology of Plants (Academic Press, London, 1977).Google Scholar
Hassell, M. P.. Density-dependence in single-species populations. J. Anim. Ecol. 44 (1975), 283295.CrossRefGoogle Scholar
Hou, Z.. Geometric method for global stability and repulsion in Kolmogorov systems. Dynam. Syst. 34 (2019), 456483.CrossRefGoogle Scholar
Hou, Z.. Geometric method for global stability of discrete population models. Discrete Contin. Dyn. Syst. Ser. B 25 (2020), 33053334.Google Scholar
Hou, Z. and Baigent, S.. Fixed point global attractors and repellors in competitive Lotka–Volterra systems. Dynam. Syst. 26 (2011), 367390.CrossRefGoogle Scholar
Hutson, V. and Moran, W.. Persistence of species obeying difference equations. J. Math. Biol. 15 (1982), 203213.CrossRefGoogle Scholar
Kang, Y. and Yakubu, A.-A.. Weak Allee effects and species coexistence. Nonlinear Anal.: Real World Appl. 12 (2011), 33293345.Google Scholar
Kim, H. and Marlin, J. A.. Solutions of a pioneer-climax model. Can. Appl. Math. Q. 7 (1999), 143158.Google Scholar
Kon, R.. Convex dominates concave: an exclusion principle in discrete-time Kolmogorov systems. Proc. Am. Math. Soc. 134 (2006), 30253034.CrossRefGoogle Scholar
May, R. M., Ecosystem patterns in randomly fluctuating environments, In: Progress in Theoretical Biology (ed., Rosen and Snell) (Academic Press, New York, 1974).Google Scholar
May, R. M. and Oster, G. F.. Bifurcations and dynamic complexity in simple ecological models. Am. Nat. 110 (1976), 573599.CrossRefGoogle Scholar
Moran, P.. Some remarks on animal population dynamics. Biometrics 6 (1950), 250258.CrossRefGoogle ScholarPubMed
Niu, L. and Ruiz-Herrera, A.. Trivial dynamics in discrete-time systems: carrying simplex and translation arcs. Nonlinearity 31 (2018), 26332650.CrossRefGoogle Scholar
Ortega, R.. Periodic Differential Equations in the Plane: A Topological Perspective (Walter de Gruyter GmbH & Co KG, 2019).CrossRefGoogle Scholar
Ricker, W. E.. Stock and recruitment. J. Fisheries Board Canada 11 (1954), 559623.CrossRefGoogle Scholar
Ruiz-Herrera, A.. Permanence of two species and fixed point index. Nonlinear Anal.: Theory, Methods Appl. 74 (2011), 146153.CrossRefGoogle Scholar
Ruiz-Herrera, A.. Topological criteria of global attraction with applications in population dynamics. Nonlinearity 25 (2012), 28232841.CrossRefGoogle Scholar
Ruiz-Herrera, A.. Exclusion and dominance in discrete population models via the carrying simplex. J. Differ. Eq. Appl. 19 (2013), 96113.CrossRefGoogle Scholar
Ruiz-Herrera, A.. Attraction in nonmonotone planar systems and real-life models. J. Dyn. Diff. Equat. (2020).Google Scholar
Selgrade, J. F. and Namkoong, G.. Stable periodic behavior in a pioneer-climax model. Nat. Resour. Model. 4 (1990), 215227.CrossRefGoogle Scholar
Sumner, S.. Hopf bifurcation in pioneer-climax competing species models. Math. Biosci. 137 (1996), 124.CrossRefGoogle ScholarPubMed
Figure 0

FIGURE 1. Representation of the function $g$ in system (3.7) with $r=0.5$, $a=0.45$ and $b=3$.