Chaos Induced by Snap-Back Repeller in a Two Species Competitive Model

In this paper, we investigate the complex dynamics of two-species Ricker-type discrete-time competitive model. We perform a local stability analysis for the fixed points and we will discuss about its persistence for boundary fixed points. This system inherits the dynamics of one-dimensional Ricker model such as cascade of period-doubling bifurcation, periodic windows and chaos. We explore the existence of chaos for the equilibrium points for a specific case of this system using Marotto theorem and proving the existence of snap-back repeller. We use several dynamical systems tools to demonstrate the qualitative behaviors of the system.

Keywords

Share and Cite:

Azizi, T. and Alali, B. (2020) Chaos Induced by Snap-Back Repeller in a Two Species Competitive Model. American Journal of Computational Mathematics, 10, 311-328. doi: 10.4236/ajcm.2020.102017.

1. Introduction

When we study the evolution of population dynamics, two major types of mathematical modelings can be used: the continuous-time dynamical systems and the discrete-time dynamical systems. For the purpose of modeling small size population and non-overlapping generations, the discrete time systems are the appropriate model . There are so many studies that have been worked on discovering complex behaviors of discrete competitive model during the last decades    . There are not many of these studies which are concerning about the existence of chaos in higher dimensional discrete dynamical systems. Chaos and chaos synchronizations have attracted many researchers for many years  . In 1975, Li and York provided a simple criterion for chaos in one dimensional discrete dynamical systems, “period three implies chaos” . This definition is the first description of chaos. Although, a precise definition of chaos was presented by their work, however, F.R. Marotto mentioned that the essential properties of chaos are the following: 1) there exist an infinite number of periodic solutions of various periods; 2) there exists an uncountably infinite set of points which exhibit random behavior; and 3) there is a highly sensitivity to initial conditions   . Marotto extended Li-York’s chaos in one-dimension to multi-dimension through introducing the notion of snapback repeller by his famous theorem in 1978 a few years after Li and York definition for chaos. To explain more, we have mentioned the Marotto’s definition for “Snap-back rappeler” and then his theorem :

Definition 1.1 (Marotto-1978) Let f be differentiable in ${B}_{{r}^{\prime }}\left(z\right)$. The point $z\in {ℝ}^{n}$ is an expanding fixed point of f in ${B}_{{r}^{\prime }}\left(z\right)$, if $f\left(z\right)=z$ and all eigenvalues of $Df\left(x\right)$ exceed 1 in norm for all $x\in {B}_{{r}^{\prime }}\left(z\right)$.

Definition 1.2 (Marotto-1978) Assume that z is an expanding fixed point of f in ${B}_{{r}^{\prime }}\left(z\right)$ for some ${r}^{\prime }>0$. Then z is said to be an snap-back repeller of f if there exists a point ${z}_{0}\in {B}_{{r}^{\prime }}\left(z\right)$ with ${z}_{0}\ne z$ and ${f}^{M}\left({z}_{0}\right)=z$ and $|D{f}^{M}\left({z}_{0}\right)|\ne 0$ for some positive integer M .

Figure 1 demonstrates the schematic diagram of snap-back repeller point.

Under the assumptions for definitions (1.1) and (1.2), the following theorem by Marotto holds.

Theorem 1.3 (Marotto-1978) If f possesses a snap back repeller, then f is chaotic in the following sense: There exist 1) a positive integer N, such that f has a point of period p, for each integer $p\ge N$, 2) a scrambled set of f, i.e., an uncountable set S containing no periodic points of f, such that

a) $f\left(S\right)\subset S$,

b) $\mathrm{lim}{\mathrm{sup}}_{n\to \infty }‖{f}^{n}\left(x\right)-{f}^{n}\left(y\right)‖>0$, for all $x,y\in S$, with $x\ne y$,

c) $\mathrm{lim}{\mathrm{sup}}_{n\to \infty }‖{f}^{n}\left(x\right)-{f}^{n}\left(y\right)‖>0$, for all $x\in S$ and periodic point y of f,

3) an uncountable subset ${S}_{0}$ of S, such that $\mathrm{lim}{\mathrm{inf}}_{n\to \infty }‖{f}^{n}\left(x\right)-{f}^{n}\left(y\right)‖=0$, for every $x,y\in {S}_{0}$ .

However, there was a minor technical flaw in his work   . Although he wanted to apply his theorem to any repelling fixed point, some of the conditions that he considered in the proof of his theorem were associated with only expanding fixed points. He incorrectly mentioned that if the absolute value for all eigenvalues of $Df\left(z\right)$ is larger than 1, then the fixed point z is an expanding fixed point of f. As we know all expanding fixed points are repelling and its converse is not true. Therefore, Marotto definition for snap-back repeller and then his proof about existence of snap-back repeller implies chaos had a minor error. Chen et al., 1998; Lin et al., 2002; Li and Chen, 2003a; discussed about the flaws of Marotto’s theorem and some of them provided several counterexamples to say that if all eigenvalues of the Jacobian $Df\left(z\right)$ at the fixed point z are greater than one in norm, we cannot say always there exists some $s>1$ and ${r}^{\prime }>0$ such that for all $x,y\in {B}_{{r}^{\prime }}\left(z\right)$, $‖f\left(x\right)-f\left(y\right)‖>s‖x-y‖$. Then they redefined the Marotto’s Theorem as the following form :

Figure 1. Snap-Back repeller schematic diagram.

Theorem 1.4 (Marotto-Li-Chen Theorem (2003)) Consider the following n-dimensional discrete dynamical system:

${x}_{n+1}=f\left({x}_{n}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}_{n}\in {ℝ}^{n},\text{\hspace{0.17em}}\text{ }n=0,1,2,\cdots$

where $f:{ℝ}^{n}\to {ℝ}^{n}$ and z is a fixed point. Also assume that

1) $f\left(x\right)$ is continuously differentiable in ${B}_{{r}^{\prime }}\left(z\right)$ for some ${r}^{\prime }>0$,

2) All eigenvalues of ${\left(Df\left(z\right)\right)}^{\text{T}}Df\left(z\right)$ are greater than 1,

3) There exists a point ${z}_{0}=\left\{x|‖x-z‖\le {r}^{\prime }$ and all eigenvalues of ${\left(Df\left(x\right)\right)}^{\text{T}}Df\left(x\right)$ are larger than 1}, with ${z}_{0}\ne z$, such that ${f}^{M}\left({z}_{0}\right)=z$ where ${f}^{i}\left({z}_{0}\right)\in {B}_{{r}^{\prime }}\left(z\right),i=0,1,2,\cdots ,M$, and the determinant $|D{f}^{M}\left({z}_{0}\right)|\ne 0$, for some positive integer M.

Then, the system is chaotic in the sense of Li-York .

Marotto refined his theorem in 2005 and he explained that a fixed point z is called a repelling fixed point under differentiable function $f:{ℝ}^{n}\to {ℝ}^{n}$ if all eigenvalues of $Df\left(z\right)$ exceed 1 in magnitude, but z is expanding only if

$‖f\left(x\right)-f\left(y\right)‖>s‖x-y‖$

where $s>1$, for all $x,y$ sufficiently close to z with $x\ne y$ (for $x,y\in {B}_{{r}^{\prime }}\left(z\right)$ ). This implies that f is a 1-1 function in ${B}_{{r}^{\prime }}\left(z\right)$ .

Definition 1.5 (Marotto-2005) Suppose z is a fixed point of f with all eigenvalues of $Df\left(z\right)$ exceeding 1 in magnitude and suppose that there exists a point ${z}_{0}\ne z$ in a repelling neighborhood of z and an integer $M>1$, such that ${x}_{M}=z$ and $\mathrm{det}\left(Df\left({x}_{k}\right)\right)\ne 0$ for $1\le k\le M$ where ${x}_{k}={f}^{k}\left({z}_{0}\right)$. Then z is called a snapback repeller of f .

He claimed that since $\mathrm{det}\left(Df\left({x}_{k}\right)\right)\ne 0$ for all $1\le k\le M$, then the homoclinic orbit is transversal in the sense that f for all $k\le M$ is 1-1 map in a neighborhood of ${x}_{k}$.

As Marotto explained in 1978, the condition $\mathrm{det}\left(Df\left({x}_{k}\right)\right)\ne 0$ guarantees the existence of the inverse of ${f}^{M}$ in ${B}_{{r}^{\prime }}\left(z\right)$. He mentioned that functions exhibit chaos and complex behavior when they possess snap-back repeller.

But what will happen that existence of a transverse homoclinic map convince us that we have chaos? As it is mentioned by many authors, a point which is in intersection of stable manifold and unstable manifold of a hyperbolic fixed point is called homoclinic point      . If stable manifolds and unstable manifold of the hyperbolic fixed point, intersect transversally, then we have transverse homoclinic point in the intersection of both manifolds. In a neighborhood of a transverse homoclinic point, our map possesses an invariant cantor set on which it is topologically conjugate to a shift map. Shift map acting on the space of bi-infinite sequences of 0’s and 1’s and it has the following properties:

A countable infinity of periodic orbits consists of orbits of all periods.

1) An uncountable infinity of non-periodic orbits.

2) A dense orbit.

Although, Wiggins in  mentioned that understanding the orbit structure of a map in that invariant Cantor set is impossible, he could show that the map in that invariant set behaves the same as shift map.

There are some researches which have more details about small neighborhood of a point on the homoclinic orbit     . The homoclinic orbits and homoclinic bifurcations which occur in continuous time dynamical systems has been studied widely by   are using in discrete time systems by defining the Poincare map . In 2011, L. Gardini et al. showed that critical homoclinic orbits lead to snap-back repellers and chaos too .

As Gardini et al. discussed, in non-invertible maps homoclinic orbits may be associate with expanding fixed points and or expanding cycles. Also, they mentioned that in the neighborhood of such homoclinic orbits, there exists an invariant set on which the map is chaotic. They even for the case that They proved that even if $\mathrm{det}\left(Df\left({x}_{k}\right)\right)=0$, there are some situations in which the map is chaotic although Marotto theorem does not work. Laura et al., provide a definition for non-critical expanding fixed points and then they defined when a homoclinic orbit is critical. They used those definitions to prove a generalization of Marotto theorem in the case that we do not need the homoclinic orbit to be non-degenerate :

Theorem 1.6 (L. Gardini. et al., (2011)) Let f be a piecewise smooth non-invertible map, $f:X\to X,\text{ }X\in {ℝ}^{n}$. Let p be an expanding fixed point of f and $O\left(p\right)$ a noncritical homoclinic orbit of p. Then in any neighborhood of $O\left(p\right)$, there exists an invariant cantor like set $\Lambda$ on which the f is chaotic .

In , Gardini studied the homoclinic bifurcations in n dimensional endomorphisms (maps with a nonunique inverse) which are associated to expanding periodic orbits. The study of chaos for these kinds of map in one dimension was studied by Mira in 1987 . Since, this topic is out of the discussion for this paper, so we avoid going through that. In this paper, we study the local dynamics of a two-species Ricker competitive model with four biological parameters. We will conduct a local stability analysis to study the local dynamics of the steady states of the system. We will use the persistence theory to study the global dynamics of the system. To study the chaotic dynamics of the system, we focus on a specific case with only three biological parameters. We provide the condition under which Marotto theorem works for positive fixed points of this new system. Furthermore, this model does not have a Neimark-Sacker bifurcation and inherits the same dynamics as one dimensional Ricker model. We will numerically demonstrate the local and qualitative dynamics of the system using several dynamical system tools.

2. The Two-Species Ricker Competitive Model and Its Local Dynamics

The Ricker model is a well known population model which demonstrates stable, periodic and non-periodic and complex nonlinear dynamics  . Here, we consider a two-species Ricker model which is a special case of model (2) in  and has the following form:

${f}_{1}={X}_{1}\left(n+1\right)={X}_{1}\left(n\right){\text{e}}^{{r}_{1}\left(1-\frac{{X}_{1}\left(n\right)}{k}-{X}_{2}\left(n\right)\right)}$ (1)

${f}_{2}={X}_{2}\left(n+1\right)={X}_{2}\left(n\right){\text{e}}^{{r}_{2}\left(1-\frac{{X}_{2}\left(n\right)}{l}-{X}_{1}\left(n\right)\right)}$ (2)

here, ${X}_{1}$ demonstrates the population size of the first species, ${X}_{2}$ represents the population size of the second species, ${r}_{1}$ and ${r}_{2}$ are the intrinsic growth rate, k and l the carrying capacity of the environment.

The Jacobian matrix for (1)-(2) has the form

$J:=\left[\begin{array}{cc}\frac{\partial {f}_{1}}{\partial {X}_{1}}& \frac{\partial {f}_{1}}{\partial {X}_{2}}\\ \frac{\partial {f}_{2}}{\partial {X}_{1}}& \frac{\partial {f}_{2}}{\partial {X}_{2}}\end{array}\right]$ (3)

where

$\frac{\partial {f}_{1}}{\partial {X}_{1}}=\left(1-\frac{{r}_{1}{X}_{1}}{k}\right)\mathrm{exp}\left({r}_{1}\left(1-\frac{{X}_{1}}{k}-{X}_{2}\right)\right)$

$\frac{\partial {f}_{1}}{\partial {X}_{2}}=-{r}_{1}{X}_{1}\mathrm{exp}\left({r}_{1}\left(1-\frac{{X}_{1}}{k}-{X}_{2}\right)\right)$

$\frac{\partial {f}_{2}}{\partial {X}_{1}}=-{r}_{2}{X}_{2}\mathrm{exp}\left({r}_{2}\left(1-\frac{{X}_{2}}{l}-{X}_{1}\right)\right)$

$\frac{\partial {f}_{2}}{\partial {X}_{2}}=\left(1-\frac{{r}_{2}{X}_{2}}{l}\right)\mathrm{exp}\left({r}_{2}\left(1-\frac{{X}_{2}}{l}-{X}_{1}\right)\right)$

Then, at the origin we have

${J|}_{\left(0,0\right)}=\left(\begin{array}{cc}{\text{e}}^{{r}_{1}}& 0\\ 0& {\text{e}}^{{r}_{2}}\end{array}\right)$

and for the fixed point $\left(k,0\right)$ we have

${J|}_{\left(k,0\right)}=\left(\begin{array}{cc}1-{r}_{1}& -k{r}_{1}\\ 0& {\text{e}}^{{r}_{2}\left(1-l\right)}\end{array}\right)$

and for the fixed point $\left(0,l\right)$ we have

${J|}_{\left(0,l\right)}=\left(\begin{array}{cc}{\text{e}}^{{r}_{1}\left(1-k\right)}& 0\\ -l{r}_{2}& 1-{r}_{2}\end{array}\right)$

and for the positive fixed point $\left({X}_{1}^{*},{X}_{2}^{*}\right)=\left(\frac{k\left(1-l\right)}{1-kl},\frac{l\left(1-k\right)}{1-kl}\right)$, we have

${J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}=\left(\begin{array}{cc}\frac{-1+kl+{r}_{1}-{r}_{1}l}{-1+kl}& \frac{-k\left(-1+l\right){r}_{1}}{-1+kl}\\ \frac{-k\left(-1+k\right){r}_{2}}{-1+kl}& \frac{-1+kl+{r}_{2}-{r}_{2}k}{-1+kl}\end{array}\right)$ (4)

Proposition 2.1 The local stability analysis results for the fixed points $\left(0,0\right)$, $\left(k,0\right)$, $\left(0,l\right)$ of (1)-(2) are summarized as below:

1) The equilibrium point $\left(0,0\right)$ is always an unstable fixed point.

2) The equilibrium point $\left(k,0\right)$ for $l<1$ and $0<{r}_{1}<2$, has a stable manifold in ${X}_{1}$ direction and an unstable manifold in ${X}_{2}$ direction and is a saddle point. Also, $\left(k,0\right)$ for $l>1$ and $0<{r}_{1}<2$, has a stable manifold in ${X}_{1}$ direction and a stable manifold in ${X}_{2}$ direction and is a stable node. Moreover, $\left(k,0\right)$ for $l<1$ and ${r}_{1}>2$, has an unstable manifold in ${X}_{1}$ direction and an unstable manifold in ${X}_{2}$ direction and is an unstable node. Finally, $\left(k,0\right)$ for $l>1$ and ${r}_{1}>2$, has an unstable manifold in ${X}_{1}$ direction and a stable manifold in ${X}_{2}$ direction and is a saddle point.

3) The equilibrium point $\left(0,l\right)$ for $k<1$ and $0<{r}_{1}<2$, has a stable manifold in ${X}_{2}$ direction and an unstable manifold in ${X}_{1}$ direction and is a saddle point. Also, $\left(0,l\right)$ for $k>1$ and $0<{r}_{2}<2$, has a stable manifold in ${X}_{1}$ direction and a stable manifold in ${X}_{2}$ direction and is a stable node. Moreover, $\left(0,l\right)$ for $k<1$ and ${r}_{2}>2$, has an unstable manifold in ${X}_{1}$ direction and an unstable manifold in ${X}_{2}$ direction and is an unstable node. Finally, $\left(0,l\right)$ for $k>1$ and ${r}_{1}>2$, has an unstable manifold in ${X}_{2}$ direction and a stable manifold in ${X}_{1}$ direction and is a saddle point.

3. Global Stability Analysis Using Persistence Theory

3.1. Boundedness of the System Solutions

To study the global stability of the equilibrium points of system, at first we prove that all solutions in the first quadrant ${ℝ}_{+}^{2}$ are eventually bounded.

Theorem 3.1 For ${r}_{1},{r}_{2}>0$, $k,l>0$ and initial conditions in the first quadrant ${ℝ}_{+}^{2}$, i.e. ${X}_{1}\left(0\right)>0$ and ${X}_{2}\left(0\right)>0$, for the system of (1)-(2) we have: ${X}_{1}>0$ and ${X}_{2}>0$ for all $n\in {ℤ}^{+}$. In addition, we can find some positive number M, such that $ma{x}_{n\in {ℤ}^{+}}\left\{{X}_{1}\left(n\right),{X}_{2}\left(n\right)\right\}\le M$.

Proof. By induction.

Since ${X}_{1}\left(0\right)>0$ we have $\mathrm{exp}\left({r}_{1}\left(1-\frac{{X}_{1}\left(0\right)}{k}\right)\right)>0$, hence

${X}_{1}\left(1\right)={X}_{1}\left(0\right){\text{e}}^{{r}_{1}\left(1-\frac{{X}_{1}\left(0\right)}{k}-{X}_{2}\left(0\right)\right)}<{X}_{1}\left(0\right){\text{e}}^{{r}_{1}\left(1-\frac{{X}_{1}\left(0\right)}{k}\right)}>0$

Assume that for $n\le j$, we have ${X}_{1}\left(j\right)>0$. Then for $n=j+1$ we have

${X}_{1}\left(j+1\right)={X}_{1}\left(j\right){\text{e}}^{{r}_{1}\left(1-\frac{{X}_{1}\left(j\right)}{k}-{X}_{2}\left(j\right)\right)}>0$

Therefore ${X}_{1}\left(n\right)>0$ for any $n\in {ℤ}^{+}$. Similarly, since ${X}_{1}\left(0\right)>0$ and ${X}_{2}\left(0\right)>0$, we automatically have $\mathrm{exp}\left({r}_{2}\left(1-\frac{{X}_{2}\left(0\right)}{j}\right)\right)>0$ is positive. Hence,

${X}_{2}\left(1\right)={X}_{2}\left(0\right){\text{e}}^{{r}_{2}\left(1-\frac{{X}_{2}\left(0\right)}{j}-{X}_{1}\left(0\right)\right)}<{X}_{2}\left(0\right){\text{e}}^{{r}_{2}\left(1-\frac{{X}_{2}\left(0\right)}{j}\right)}>0$

Assume that for $n\le j$, we have ${X}_{2}\left(j\right)>0$. Then for $n=j+1$ we have

${X}_{2}\left(j+1\right)={X}_{2}\left(j\right){\text{e}}^{{r}_{2}\left(1-\frac{{X}_{2}\left(l\right)}{j}-{X}_{1}\left(j\right)\right)}>0$

Therefore ${X}_{2}\left(n\right)>0$ for any $n\in {ℤ}^{+}$.

To find an upper bound, we know,

${X}_{1}\left(n+1\right)={X}_{1}\left(n\right){\text{e}}^{{r}_{1}\left(1-\frac{{X}_{1}\left(n\right)}{k}\right)}\le \underset{x\in {ℝ}^{+}}{max}\left\{f\left(x\right)\right\}$

If we define ${f}_{1}\left(x\right)=x{\text{e}}^{{r}_{1}\left(1-\frac{x}{k}\right)}$, then ${{f}^{\prime }}_{1}\left(x\right)=\left(1-\frac{{r}_{1}x}{k}\right){\text{e}}^{{r}_{1}\left(1-\frac{x}{k}\right)}$ and ${f}_{1}\left(x\right)$ has critical points at $x=\frac{k}{{r}_{1}}$. Since ${{f}^{\prime }}_{1}\left(x\right)>0$ if $x<\frac{k}{{r}_{1}}$ and ${{f}^{\prime }}_{1}\left(x\right)<0$ if $x>\frac{k}{{r}_{1}}$, then $x=\frac{k}{{r}_{1}}$ is the maximal point of ${f}_{1}\left(x\right)$, i.e. ${\mathrm{max}}_{x\in {ℝ}^{+}}\left\{{f}_{1}\left(x\right)\right\}={f}_{1}\left(\frac{k}{{r}_{1}}\right)$. Hence,

${x}_{1}\left(n+1\right)={X}_{1}\left(n\right){\text{e}}^{{r}_{1}\left(1-\frac{{X}_{1}\left(n\right)}{k}-{X}_{2}\left(n\right)\right)}\le {f}_{1}\left(\frac{k}{{r}_{1}}\right)=\frac{k{\text{e}}^{{r}_{1}-1}}{{r}_{1}}={M}_{1}$

Similarly, we define ${f}_{2}\left(x\right)=x{\text{e}}^{{r}_{2}\left(1-\frac{x}{l}\right)}$, then ${{f}^{\prime }}_{2}\left(x\right)=\left(1-\frac{{r}_{2}x}{l}\right){\text{e}}^{{r}_{2}\left(1-\frac{x}{l}\right)}$ and ${f}_{2}\left(x\right)$ has critical points at $x=\frac{l}{{r}_{2}}$. Since ${{f}^{\prime }}_{2}\left(x\right)>0$ if $x<\frac{l}{{r}_{2}}$ and ${{f}^{\prime }}_{2}\left(x\right)<0$ if $x>\frac{l}{{r}_{2}}$, then $x=\frac{l}{{r}_{2}}$ is the maximal point of ${f}_{2}\left(x\right)$, i.e. ${\mathrm{max}}_{x\in {ℝ}^{+}}\left\{{f}_{2}\left(x\right)\right\}={f}_{2}\left(\frac{l}{{r}_{2}}\right)$.

${X}_{2}\left(n+1\right)={X}_{2}\left(n\right){\text{e}}^{{r}_{2}\left(1-\frac{{X}_{2}\left(n\right)}{l}-{X}_{1}\left(n\right)\right)}\le {f}_{2}\left(\frac{l}{{r}_{2}}\right)=\frac{l{\text{e}}^{{r}_{2}-1}}{{r}_{2}}={M}_{2}$

Therefore, we can find some positive number $M=\mathrm{max}\left\{{M}_{1},{M}_{2}\right\}$, such that $ma{x}_{n\in {ℤ}^{+}}\left\{{X}_{1}\left(n\right),{X}_{2}\left(n\right)\right\}\le M$.

3.2. Persistence of the Species

To work on global stability, we need to study the persistence theory  . Here, we consider two cases:

1) Persistence of system corresponding to $\left(k,0\right)$.

2) Persistence of system corresponding to $\left(0,l\right)$.

3.2.1. Case 1: Persistence of System Corresponding to (k, 0)

For the first case, we have:

$P=\left\{\left({X}_{1},{X}_{2}\right):{X}_{1}\ge 0,{X}_{2}\ge 0\right\}$

${P}_{k,0}=\left\{\left({X}_{1},{X}_{2}\right)\in P:{X}_{1}>0\right\}$

$\partial {P}_{k,0}=P\{P}_{k,0}$

Proposition 3.2 The system is uniformly persistent with respect to $\left({P}_{k,0},\partial {P}_{k,0}\right)$.

Proof. Here, $\partial {P}_{k,0}$ is closed in P. For any positive solution of $\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)$ of the system, as we proved in theorem (3.1), we have

${X}_{1}\left(n+1\right)\le {X}_{1}\left(n\right){\text{e}}^{{r}_{1}\left(1-\frac{{X}_{1}\left(n\right)}{k}\right)}\le \underset{{X}_{1}\in {ℝ}^{+}}{\mathrm{max}}\left\{{f}_{1}\left({X}_{1}\right)\right\}=\frac{k{\text{e}}^{{r}_{1}-1}}{{r}_{1}}={M}_{1}$

And for large enough n

${X}_{2}\left(n+1\right)\le {X}_{2}\left(n\right){\text{e}}^{{r}_{2}\left(1-\frac{{X}_{2}\left(n\right)}{l}\right)}\le \underset{{X}_{2}\in {ℝ}^{+}}{\mathrm{max}}\left\{{f}_{2}\left({X}_{2}\right)\right\}=\frac{l{\text{e}}^{{r}_{2}-1}}{{r}_{2}}={M}_{2}$

Therefore, system (1)-(2) is point dissipative. Assume for all $n\ge 0$

$\begin{array}{l}{Y}_{\partial }=\left\{\left({X}_{1}\left(0\right),{X}_{2}\left(0\right)\right):\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)\text{\hspace{0.17em}}\text{satisfies}\text{\hspace{0.17em}}\text{the}\text{\hspace{0.17em}}\text{system}\text{\hspace{0.17em}}\text{equations}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)\in \partial {P}_{k,0}\right\}\end{array}$

We see that

${Y}_{\partial }=\left\{\left(0,{X}_{2}\right):{X}_{2}\ge 0\right\}=\partial {P}_{k,0}$

Moreover, $\left(0,0\right)$ is the unique equilibrium in ${Y}_{\partial }$. Define ${W}^{s}\left(0,0\right)$ to be the stable manifold for $\left(0,0\right)$. We show that

${W}^{s}\left(0,0\right)\cap {P}_{k,0}=\varnothing$

Assume that in contradiction, there exist a solution $\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)$ of system with ${X}_{1}\left(n\right)>0$ such that

$\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)\to \left(0,0\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{as}\text{\hspace{0.17em}}\text{ }n\to \infty$

Then, for large n we have

${X}_{1}\left(n+1\right)>{X}_{1}\left(n\right){\text{e}}^{{r}_{1}/2}$

Since ${r}_{1}>0$, it follows that ${X}_{1}\left(n\right)\to \infty$ as $n\in \infty$ and contradiction. Also, every orbit in ${Y}_{\partial }$ tends to $\left(0,0\right)$ as $n\to \infty$. It means that $\left(0,0\right)$ is an isolated invariant set in P and acyclic in ${Y}_{\partial }$. Note that ${Y}_{\partial }$ repels uniformly the solution of systems with positive ${X}_{1}\left(n\right)$  . It follows that there is ${s}_{1}>0$ such that ${X}_{1}\left(n\right)>{s}_{1}$ for large enough n.

Theorem 3.3. There exists ${s}_{1}>0$ such that for any ${X}_{1}\left(0\right)>0$ we have

${s}_{1}<{X}_{1}\left(n\right)<\frac{k{\text{e}}^{{r}_{1}-1}}{{r}_{1}}$

Proof. By proposition (3.2).

Theorem 3.4 All solutions $\left\{\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)\right\}$ of system with ${X}_{1}\left(0\right)>0$ and ${X}_{2}\left(0\right)\ge 0$, for $l>1$ and $0<{r}_{1}<2$, are decreasing to the fixed point $\left(k,0\right)$, i.e.

$\underset{n\to +\infty }{\mathrm{lim}}{X}_{1}\left(n\right)=k,\text{ }\underset{n\to +\infty }{\mathrm{lim}}{X}_{2}\left(n\right)=0$

Proof. By proposition (3.2) and theorem (3.3).

3.2.2. Case 2: Persistence of System Corresponding to (0, l)

For this case, we have:

$Q=\left\{\left({X}_{1},{X}_{2}\right):{X}_{1}\ge 0,{X}_{2}\ge 0\right\}$

${Q}_{0,l}=\left\{\left({X}_{1},{X}_{2}\right)\in Q:{X}_{2}>0\right\}$

$\partial {Q}_{0,l}=Q\{Q}_{0,l}$

Proposition 3.5 The system is uniformly persistent with respect to $\left({Q}_{0,l},\partial {Q}_{0,l}\right)$.

Proof. Here, $\partial {Q}_{0,l}$ is closed in Q. Similarly, for any positive solution of $\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)$ of the system (1)-(2), similar to theorem (3.1), we can write

${X}_{1}\left(n+1\right)\le {X}_{1}\left(n\right){\text{e}}^{{r}_{1}\left(1-\frac{{X}_{1}\left(n\right)}{k}\right)}\le \underset{{X}_{1}\in {ℝ}^{+}}{\mathrm{max}}\left\{f\left({X}_{1}\right)\right\}=\frac{k{\text{e}}^{{r}_{1}-1}}{{r}_{1}}={M}_{1}$

For large enough n

${X}_{2}\left(n+1\right)\le {X}_{2}\left(n\right){\text{e}}^{{r}_{2}\left(1-\frac{{X}_{2}\left(n\right)}{l}\right)}\le \underset{{X}_{2}\in {ℝ}^{+}}{\mathrm{max}}\left\{f\left({X}_{2}\right)\right\}=\frac{l{\text{e}}^{{r}_{2}-1}}{{r}_{2}}={M}_{2}$

Thus, system (1)-(2) is point dissipative. Now, for all $n\ge 0$, we set

$\begin{array}{l}{L}_{\partial }=\left\{\left({X}_{1}\left(0\right),{X}_{2}\left(0\right)\right):\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)\text{\hspace{0.17em}}\text{satisfies}\text{\hspace{0.17em}}\text{the}\text{\hspace{0.17em}}\text{systeme}\text{\hspace{0.17em}}\text{quations}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)\in \partial {Q}_{0,l}\right\}\end{array}$

for which

${L}_{\partial }=\left\{\left({X}_{1},0\right):{X}_{2}\ge 0\right\}=\partial {Q}_{0,l}$

Moreover, $\left(0,0\right)$ is the unique equilibrium in ${L}_{\partial }$. Set ${W}^{s}\left(0,0\right)$ to be the stable manifold for $\left(0,0\right)$. We prove that

${W}^{s}\left(0,0\right)\cap {Q}_{0,l}=\varnothing$

By contradiction, there exist a solution $\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)$ of system with ${X}_{2}\left(n\right)>0$ such that

$\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)\to \left(0,0\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{as}\text{\hspace{0.17em}}\text{ }n\to \infty$

For large n we have

${X}_{2}\left(n+1\right)>{X}_{2}\left(n\right){\text{e}}^{{r}_{2}/2}$

Since ${r}_{2}>0$, it leads to ${X}_{2}\left(n\right)\to \infty$ as $n\in \infty$ which is a contradiction. Also, every orbit in ${L}_{\partial }$ tends to $\left(0,0\right)$ as $n\to \infty$. It implies that $\left(0,0\right)$ is an isolated invariant set in Q and acyclic in ${L}_{\partial }$. Here, ${l}_{\partial }$ repels uniformly the solutions of system with positive ${X}_{2}\left(n\right)$  . It follows that there is ${s}_{2}>0$ such that ${X}_{2}\left(n\right)>{s}_{2}$ for large enough n.

Theorem 3.6 There exists ${s}_{2}>0$ such that for any ${X}_{2}\left(0\right)>0$ we have

${s}_{2}<{X}_{2}\left(n\right)<\frac{l{\text{e}}^{{r}_{2}-1}}{{r}_{2}}$

Proof. By proposition (3.5).

Theorem 3.7 All solutions $\left\{\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)\right\}$ of system with ${X}_{1}\left(0\right)\ge 0$ and ${X}_{2}\left(0\right)>0$, for $k>1$ and $0<{r}_{2}<2$, are decreasing to the fixed point $\left(0,l\right)$, i.e.

$\underset{n\to +\infty }{\mathrm{lim}}{X}_{1}\left(n\right)=0,\text{ }\underset{n\to +\infty }{\mathrm{lim}}{X}_{2}\left(n\right)=l$

Proof. By proposition (3.5) and theorem (3.6).

Finally, we have the following result

Theorem 3.8 If there are positive constants ${s}_{1},{s}_{2}>0$ and ${M}_{1},{M}_{2}>0$ such that the solution $\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)$ of system satisfies

$0<{s}_{1}\le \underset{n\to +\infty }{\mathrm{lim}}\mathrm{inf}{X}_{1}\left(n\right)\le \underset{n\to +\infty }{\mathrm{lim}}\mathrm{sup}{X}_{1}\left(n\right)\le {M}_{1}=\frac{k{\text{e}}^{{r}_{1}-1}}{{r}_{1}}$

$0<{s}_{2}\le \underset{n\to +\infty }{\mathrm{lim}}\mathrm{inf}{X}_{2}\left(n\right)\le \underset{n\to +\infty }{\mathrm{lim}}\mathrm{sup}{X}_{2}\left(n\right)\le {M}_{2}=\frac{l{\text{e}}^{{r}_{2}-1}}{{r}_{2}}$

Then, system (1)-(2) is persistent. If system is not persistent, it is called non-persistent smith2011dynamical.

4. Application of Snap-Back Repeller and Marroto Chaos in Study of Chaotic Dynamics of System

In this section, we explore analytically chaos in the sense of Marotto for a specific case of model (1)-(2). Without loss of generality, we consider $k=l$, then we have

$F:=\left\{\begin{array}{l}{g}_{1}\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)={X}_{1}\left(n\right)\mathrm{exp}\left({r}_{1}\left(1-\frac{{X}_{1}\left(n\right)}{k}-{X}_{2}\left(n\right)\right)\right)\\ {g}_{2}\left({X}_{1}\left(n\right),{X}_{2}\left(n\right)\right)={X}_{2}\left(n\right)\mathrm{exp}\left({r}_{2}\left(1-\frac{{X}_{2}\left(n\right)}{k}-{X}_{1}\left(n\right)\right)\right)\end{array}$ (5)

The Jacobian matrix for (5) has the form

$J:=\left[\begin{array}{cc}\frac{\partial {g}_{1}}{\partial {X}_{1}}& \frac{\partial {g}_{1}}{\partial {X}_{2}}\\ \frac{\partial {g}_{2}}{\partial {X}_{1}}& \frac{\partial {g}_{2}}{\partial {X}_{2}}\end{array}\right]$ (6)

where

$\frac{\partial {g}_{1}}{\partial {X}_{1}}=\left(1-\frac{{r}_{1}{X}_{1}}{k}\right)\mathrm{exp}\left({r}_{1}\left(1-\frac{{X}_{1}}{k}-{X}_{2}\right)\right)$ (7)

$\frac{\partial {g}_{1}}{\partial {X}_{2}}=-{r}_{1}{X}_{1}\mathrm{exp}\left({r}_{1}\left(1-\frac{{X}_{1}}{k}-{X}_{2}\right)\right)$ (8)

$\frac{\partial {g}_{2}}{\partial {X}_{1}}=-{r}_{2}{X}_{2}\mathrm{exp}\left({r}_{2}\left(1-\frac{{X}_{2}}{k}-{X}_{1}\right)\right)$ (9)

$\frac{\partial {g}_{2}}{\partial {X}_{2}}=\left(1-\frac{{r}_{2}{X}_{2}}{k}\right)\mathrm{exp}\left({r}_{2}\left(1-\frac{{X}_{2}}{k}-{X}_{1}\right)\right)$ (10)

For this specific case, we have four fixed points $\left(0,0\right)$, $\left(k,0\right)$, $\left(0,k\right)$ and $\left({X}_{1}^{*},{X}_{2}^{*}\right)=\left(\frac{k}{k+1},\frac{k}{k+1}\right)$. At $\left(0,0\right)$ we have

${J|}_{\left(0,0\right)}=\left(\begin{array}{cc}{\text{e}}^{{r}_{1}}& 0\\ 0& {\text{e}}^{{r}_{2}}\end{array}\right)$

and at $\left(k,0\right)$ we have

${J|}_{\left(k,0\right)}=\left(\begin{array}{cc}1-{r}_{1}& -k{r}_{1}\\ 0& {\text{e}}^{{r}_{2}\left(1-k\right)}\end{array}\right)$

and also for the fixed point $\left(0,k\right)$ we have

${J|}_{\left(0,k\right)}=\left(\begin{array}{cc}{\text{e}}^{{r}_{1}\left(1-k\right)}& 0\\ -k{r}_{2}& 1-{r}_{2}\end{array}\right)$

and finally for the positive fixed point $\left({X}_{1}^{*},{X}_{2}^{*}\right)=\left(\frac{k}{k+1},\frac{k}{k+1}\right)$, we have

${J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}=\left(\begin{array}{cc}\frac{k+1-{r}_{1}}{k+1}& \frac{-k{r}_{1}}{k+1}\\ \frac{-k{r}_{2}}{k+1}& \frac{k+1-{r}_{2}}{k+1}\end{array}\right)$ (11)

where

$det\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)=-\frac{k{r}_{1}{r}_{2}-{r}_{1}{r}_{2}-k+{r}_{1}+{r}_{2}-1}{k+1}$ (12)

$\text{tr}\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)=\frac{2k+2-{r}_{2}-{r}_{1}}{k+1}$ (13)

and also, characteristic polynomial has the form

$P\left(X\right):={X}^{2}-\frac{2k+2-{r}_{2}-{r}_{1}}{k+1}X-\frac{k{r}_{1}{r}_{2}-{r}_{1}{r}_{2}-k+{r}_{1}+{r}_{2}-1}{k+1}$ (14)

Proposition 4.1 The local stability analysis results for the fixed points $\left(0,0\right)$, $\left(k,0\right)$, $\left(0,k\right)$ of (5) are summarized as below:

1) The equilibrium point $\left(0,0\right)$ is always an unstable fixed point.

2) The equilibrium point $\left(k,0\right)$ for $k<1$ and $0, has a stable manifold in ${X}_{1}$ direction and an unstable manifold in ${X}_{2}$ direction and is a saddle point. Also, $\left(k,0\right)$ for $k>1$ and $0, has a stable manifold in ${X}_{1}$ direction and a stable manifold in ${X}_{2}$ direction and is a stable node. Moreover, $\left(k,0\right)$ for $k<1$ and $r>2$, has an unstable manifold in ${X}_{1}$ direction and an unstable manifold in ${X}_{2}$ direction and is an unstable node. Finally, $\left(k,0\right)$ for $k>1$ and $r>2$, has an unstable manifold in ${X}_{1}$ direction and a stable manifold in ${X}_{2}$ direction and is a saddle point.

3) The equilibrium point $\left(0,k\right)$ for $k<1$ and $0, has a stable manifold in ${X}_{2}$ direction and an unstable manifold in ${X}_{1}$ direction and is a saddle point. Also, $\left(0,k\right)$ for $k>1$ and $0, has a stable manifold in ${X}_{1}$ direction and a stable manifold in ${X}_{2}$ direction and is a stable node. Moreover, $\left(0,k\right)$ for $k<1$ and $r>2$, has an unstable manifold in ${X}_{1}$ direction and an unstable manifold in ${X}_{2}$ direction and is an unstable node. Finally, $\left(0,k\right)$ for $k>1$ and $r>2$, has an unstable manifold in ${X}_{2}$ direction and a stable manifold in ${X}_{1}$ direction and is a saddle point.

Proposition 4.2 The local stability analysis results for the fixed points $\left({X}_{1}^{*},{X}_{2}^{*}\right)=\left(\frac{k}{k+1},\frac{k}{k+1}\right)$ of (5) are summarized as below:

1) The equilibrium point $\left({X}_{1}^{*},{X}_{2}^{*}\right)$ is an unstable fixed point if and only if

${r}_{1}{r}_{2}\left(1-k\right)+2\left(k+1\right)<\left({r}_{1}+{r}_{2}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}4\left(k+1\right)-2\left({r}_{1}+{r}_{2}\right)+{r}_{1}{r}_{2}\left(1-k\right)>0\text{ },\text{\hspace{0.17em}}\text{\hspace{0.17em}}k<1$

or

$k<\frac{{r}_{1}{r}_{2}-{r}_{1}-{r}_{2}}{{r}_{1}{r}_{2}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}4\left(k+1\right)-2\left({r}_{1}+{r}_{2}\right)+{r}_{1}{r}_{2}\left(1-k\right)>0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}k<1$

2) The equilibrium point $\left({X}_{1}^{*},{X}_{2}^{*}\right)$ is a stable fixed point if and only if

$k>\frac{{r}_{1}{r}_{2}-{r}_{1}-{r}_{2}}{{r}_{1}{r}_{2}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}4\left(k+1\right)-2\left({r}_{1}+{r}_{2}\right)+{r}_{1}{r}_{2}\left(1-k\right)>0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}k<1$

3) The equilibrium point $\left({X}_{1}^{*},{X}_{2}^{*}\right)$ is a saddle point if and only if

$4\left(k+1\right)-2\left({r}_{1}+{r}_{2}\right)+{r}_{1}{r}_{2}\left(1-k\right)<0\text{ },\text{\hspace{0.17em}}\text{\hspace{0.17em}}k>1$

Proof. Using Theorem 1.1.1 (Linearized Stability) in .

The equilibrium point $\left({X}_{1}^{*},{X}_{2}^{*}\right)$ is an unstable fixed point if and only if $|det\left(J\right)|>1$ and $|\text{tr}\left(J\right)|<|1+det\left(J\right)|$. $\text{tr}\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)-det\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)-1<0$ gives us:

$\frac{{r}_{1}{r}_{2}\left(k-1\right)}{k+1}<0\text{ }\to \text{ }k<1$ (15)

Also, $\text{tr}\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)+det\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)+1<0$ gives us:

$\frac{4\left(k+1\right)-2\left({r}_{1}+{r}_{2}\right)+{r}_{1}{r}_{2}\left(1-k\right)}{k+1}>0$ (16)

and $det\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)>1$ gives us

$\frac{{r}_{1}{r}_{2}\left(1-k\right)-\left({r}_{1}+{r}_{2}\right)}{k+1}<0$

that is to say

$k>\frac{{r}_{1}{r}_{2}-{r}_{1}-{r}_{2}}{{r}_{1}{r}_{2}}$

Moreover, $det\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)<-1$ gives us

$\frac{{r}_{1}{r}_{2}\left(1-k\right)+2\left(k+1\right)-\left({r}_{1}+{r}_{2}\right)}{k+1}<0$

The positive fixed point of system (5) is asymptotically stable if and only if

$|\text{tr}\left(J\right)|<1+det\left(J\right)<2$ (17)

We check (17) using (12) and (13). $\text{tr}\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)-det\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)-1<0$ and $\text{tr}\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)+det\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)+1<0$ give us (15) and (16). and $det\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)<1$ gives us

$\frac{{r}_{1}{r}_{2}\left(1-k\right)-\left({r}_{1}+{r}_{2}\right)}{k+1}<0$

that is to say

$k>\frac{{r}_{1}{r}_{2}-{r}_{1}-{r}_{2}}{{r}_{1}{r}_{2}}$

Finally, The equilibrium point $\left({X}_{1}^{*},{X}_{2}^{*}\right)$ is a saddle point if and only if ${\text{tr}}^{2}\left(J\right)-4det\left(J\right)>0$ and $|\text{tr}\left(J\right)|>|1+det\left(J\right)|$. The first condition gives us

$\frac{{\left({r}_{1}-{r}_{2}\right)}^{2}+4{k}^{2}}{{\left(k+1\right)}^{2}}>0$

which is always true. Another conditions to check are: $\text{tr}\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)-det\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)-1>0$ gives us:

$\frac{{r}_{1}{r}_{2}\left(k-1\right)}{k+1}>0$ (18)

and, $\text{tr}\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)+det\left({J|}_{\left({X}_{1}^{*},{X}_{2}^{*}\right)}\right)+1<0$ which gives us:

$\frac{4\left(k+1\right)-2\left({r}_{1}+{r}_{2}\right)+{r}_{1}{r}_{2}\left(1-k\right)}{k+1}<0$ (19)

Numerical simulations, including bifurcation diagrams and time series display that this model demonstrates chaotic oscillations after a cascade of period-doubling bifurcations. As we can see in Figure 2, there are chaotic regions which are embedded in periodic windows regions. The periodic behaviors which appear alternately in the chaotic area, contain a copy of bifurcation diagram and it is repeating when we are changing the bifurcation parameter r. The bifurcation diagram for system (5) with respect to r displays the same qualitative dynamics for different values of k. Moreover, we have run bifurcation analysis with respect to k with different r values in Figure 3.

Also, if we look at Figure 4, at first, the equilibrium point is stable, when we increase r, it loses stability, from one cycle to two cycles, and produces a flip bifurcation. As r continues to increase, periodic oscillations are observed with periods 4, …, which eventually leads to chaos.

To prove the existence of chaos for the map (5) in the sense of Marotto, we need to find the conditions under which the fixed point ${Z}^{*}=\left({X}_{1}^{*},{X}_{2}^{*}\right)$ of the system is a snap-back repeller. According to definition (1.5) and Figure 1, we need to find a neighborhood ${B}_{{r}^{\prime }}\left({Z}^{*}\right)$ of ${Z}^{*}$ in which all eigenvalues have absolute value more than one. Now, we give the following lemma which we need that to prove chaos in the sense of Marotto for positive fixed point ${Z}^{*}=\left({X}_{1}^{*},{X}_{2}^{*}\right)$ of map (5).

Figure 2. Bifurcation diagram of system (4.1) when $k=10$ and ${r}_{1}={r}_{2}=r$.

Figure 3. Bifurcation diagram of system (4.1) when ${r}_{1}=2.75$ and ${r}_{2}=2.5$.

Figure 4. Evolution of host population ${X}_{1}$ and its coupled ${X}_{2}$ in time for system (4.1) when $k=10$.

Lemma 4.3 Assume that the conditions of the first part of the proposition (4.2) are satisfied. The fixed point ${Z}^{*}=\left({X}_{1}^{*},{X}_{2}^{*}\right)$ of map F is called snap-back repeller if there exists a point ${Z}_{0}=\left({X}_{1},{X}_{2}\right)$ in the neighborhood of ${Z}^{*}$ such that ${Z}_{0}\ne {Z}^{*}$, $F\left({Z}_{0}\right)={Z}^{*}$, $|det\left({J|}_{\left({X}_{1},{X}_{2}\right)}\right)|\ne 0$, that is to say, at first, the following system of equations has a unique solution

$\left\{\begin{array}{l}{X}_{1}^{*}={X}_{1}\mathrm{exp}\left({r}_{1}\left(1-\frac{{X}_{1}}{k}-{X}_{2}\right)\right)\\ {X}_{2}^{*}={X}_{2}\mathrm{exp}\left({r}_{2}\left(1-\frac{{X}_{2}}{k}-{X}_{1}\right)\right)\end{array}$ (20)

and

${k}^{2}-k\text{ }\left({r}_{2}{X}_{2}+{r}_{1}{X}_{1}+{r}_{1}{r}_{2}{X}_{1}{X}_{2}\right)+{r}_{1}{r}_{2}{X}_{1}{X}_{2}\ne 0$ (21)

Then ${Z}^{*}$ for some parameter values $\left({r}_{1},{r}_{2}\right)$ and k, is a snap-back repeller for map (5).

Proof. From $|det\left({J|}_{\left({X}_{1},{X}_{2}\right)}\right)|\ne 0$, we have:

$\frac{\partial {g}_{1}}{\partial {X}_{1}}\text{ }\frac{\partial {g}_{2}}{\partial {X}_{2}}-\frac{\partial {g}_{1}}{\partial {X}_{2}}\frac{\partial {g}_{2}}{\partial {X}_{1}}\ne 0$

which

$\left(\left(1-\frac{{r}_{1}{X}_{1}}{k}\right)\left(1-\frac{{r}_{2}{X}_{2}}{k}\right)-{r}_{1}{r}_{2}{X}_{1}{X}_{2}\right){\text{e}}^{\left({r}_{1}\left(1-\frac{{X}_{1}}{k}-{X}_{2}\right)+{r}_{2}\left(1-\frac{{X}_{2}}{k}-{X}_{1}\right)\right)}\ne 0$

and it gives us

${k}^{2}-k\text{ }\left({r}_{2}{X}_{2}+{r}_{1}{X}_{1}+{r}_{1}{r}_{2}{X}_{1}{X}_{2}\right)+{r}_{1}{r}_{2}{X}_{1}{X}_{2}\ne 0$

Therefore, any solution ${Z}^{*}=\left({X}_{1}^{*},{X}_{2}^{*}\right)\ne \left({X}_{1},{X}_{2}\right)={Z}_{0}$ of system (20) which satisfies the first part of the proposition (4.2) and (21), is snap-back repeller for system (5).

Theorem 4.4 Under the assumptions of the first part of proposition (4.2) and lemma (4.3), the map (4.1) is chaotic in the sense of Li-York, which means that: There exist 1) a positive integer N, such that map (4.1) has a point of period p, for each integer $p\ge N$, 2) a scrambled set of F, i.e., an uncountable set S containing no periodic points of F, such that

a) $F\left(S\right)\subset S$,

b) $\mathrm{lim}{\mathrm{sup}}_{n\to \infty }‖{F}^{n}\left(x\right)-{F}^{n}\left(y\right)‖>0$, for all $x,y\in S$, with $x\ne y$,

c) $\mathrm{lim}{\mathrm{sup}}_{n\to \infty }‖{F}^{n}\left(x\right)-{F}^{n}\left(y\right)‖>0$, for all $x\in S$ and periodic point y of f,

3) an uncountable subset ${S}_{0}$ of S, such that $\mathrm{lim}{\mathrm{inf}}_{n\to \infty }‖{F}^{n}\left(x\right)-{F}^{n}\left(y\right)‖=0$, for every $x,y\in {S}_{0}$.

Proof. By theorem (1.3).

5. Conclusion

Studying the evolution of population models and complex dynamics of competitive models has attracted many researchers during several past decades. In this paper, we studied the complex dynamics of a two-species Ricker model which consists of four different biological parameters. We explored the stability of the origin and two other boundary fixed points using local stability theorem. Also, we provided the condition under which the solutions are bounded. We have seen that this model undergoes period doubling bifurcation but it does not show Neimark-Sacker bifurcation. We used the persistence theory to reveal the global behavior of system and we discovered the persistence of the system for two boundary fixed points. Afterward, we changed the model to a specific case with only three biological parameters and we discussed about the local stability of extinction and boundary fixed points of the system. Moreover, we discovered the chaotic dynamics of the new model using Marotto theorem. As we discussed, Marotto theorem is a rigorous theorem to study chaotic dynamics for systems with higher dimensions and can be used to study the chaotic dynamics of competitive models. We presented the conditions under which the new system undergoes snap-back repeller and as a result, it is chaotic in the sense of Li-York. Finally, we used bifurcation diagram to demonstrate the interesting dynamics of new system and the role of biological parameters r and k in appearance of different types of complicated dynamics. The new system has the same number of fixed points as the first system and the bifurcation analysis displayed the same qualitative dynamics for both species as we expected.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

  Liu, H., Zhang, K., Ye, Y., Wei, Y.M. and Ma, M. (2019) Dynamic Complexity and Bifurcation Analysis of a Host-Parasitoid Model with Allee Effect and Holling Type III Functional Response. Advances in Difference Equations, 2019, Article No.: 507. https://doi.org/10.1186/s13662-019-2430-8  Smith, H.L. (1998) Planar Competitive and Cooperative Difference Equations. Journal of Difference Equations and Applications, 3, 335-357. https://doi.org/10.1080/10236199708808108  Lus, R., Elaydi, S. and Oliveira, H. (2011) Stability of a Ricker-Type Competitive Model and the Competitive Exclusion Principle. Journal of Difference Equations and Applications, 5, 636-660. https://doi.org/10.1080/17513758.2011.581764  Gyllenberg, M., Jiang, J. and Niu, L. (2019) A Note on Global Stability of Three-Dimensional Ricker Models. Journal of Difference Equations and Applications, 25, 142-150. https://doi.org/10.1080/10236198.2019.1566459  Wu, D.Y. and Zhang, H. (2014) Bifurcation Analysis of a Two-Species Competitive Discrete Model of Plankton Allelopathy. Advances in Difference Equations, 2014, Article No.: 70. https://doi.org/10.1186/1687-1847-2014-70  Azizi, T. and Kerr, G. (2020) Chaos Synchronization in Discrete-Time Dynamical Systems with Application in Population Dynamics. Journal of Applied Mathematics and Physics, 8, 406-423. https://doi.org/10.4236/jamp.2020.83031  Azizi, T. and Kerr, G. (2020) Synchronized Cycles of Generalized Nicholson-Bailey Model. American Journal of Computational Mathematics, 10, 147-166. https://doi.org/10.4236/ajcm.2020.101009  Li, T.-Y. and Yorke, J.A. (1975) Period Three Implies Chaos. The American Mathematical Monthly, 82, 985-992. https://doi.org/10.1080/00029890.1975.11994008  Marotto, F.R. (1978) Snap-Back Repellers Imply Chaos in Rn. Journal of Mathematical Analysis and Applications, 63, 199-223. https://doi.org/10.1016/0022-247X(78)90115-4  Marotto, F.R. (1979) Chaotic Behavior in the Hénon Mapping. Journal of Mathematical Analysis and Applications, 68, 187-194. https://doi.org/10.1007/BF01418128  Marotto, F.R. (2005) On Redefining a Snap-Back Repeller. Chaos, Solitons & Fractals, 25, 25-28. https://doi.org/10.1016/j.chaos.2004.10.003  Lin, W., Ruan, J. and Zhao, W.R. (2002) On the Mathematical Clarification of the Snap-Back-Repeller in High-Dimensional Systems and Chaos in a Discrete Neural Network Model. International Journal of Bifurcation and Chaos, 12, 1129-1139. https://doi.org/10.1142/S0218127402004966  Li, C.P. and Chen, G.R. (2003) An Improved Version of the Marotto Theorem. Chaos, Solitons & Fractals, 18, 69-77. https://doi.org/10.1016/S0960-0779(02)00605-7  Gardini, L. (1994) Homoclinic Bifurcations in n-Dimensional Endomorphisms, Due to Expanding Periodic Points. Nonlinear Analysis: Theory, Methods & Applications, 23, 1039-1089. https://doi.org/10.1016/0362-546X(94)90198-8  Gardini, L., Sushko, I., Avrutin, V. and Schanz, M. (2011) Critical Homoclinic Orbits Lead to Snap-Back Repellers. Chaos, Solitons & Fractals, 44, 433-449. https://doi.org/10.1016/j.chaos.2011.03.004  Wiggins, S. (2003) Introduction to Applied Nonlinear Dynamical Systems and Chaos. Princeton University Press, Springer Science & Business Media, Princeton, 2.  Balibrea, F. and Cascales, A. (2012) Li-Yorke Chaos in Perturbed Rational Difference Equations. International Conference on Difference Equations and Applications, 44, 49-61. https://doi.org/10.1007/978-3-662-52927-0_4  Panchuk, A. (2016) Some Aspects on Global Analysis of Discrete Time Dynamical Systems. In: Qualitative Theory of Dynamical Systems, Tools and Applications for Economic Modelling, Springer, Cham, 161-186. https://doi.org/10.1007/978-3-319-33276-5_2  Smale, S. (2000) Diffeomorphisms with Many Periodic Points. The Collected Papers of Stephen Smale: Volume 2. World Scientific, 44, 636-653.  Moser, J. (2001) Stable and Random Motions in Dynamical Systems: With Special Emphasis on Celestial Mechanics. Princeton University Press, Princeton, 1. https://doi.org/10.1515/9781400882694  Robinson, C. (1998) Dynamical Systems: Stability, Symbolic Dynamics, and Chaos. CRC Press, Boca Raton.  Guckenheimer, J. and Holmes, P. (1984) Nonlinear Oscillations, Dynamical Systems and Bifurcations of Vector Fields. Journal of Applied Mechanics, 51, 947. https://doi.org/10.1115/1.3167759  Shilnikov, L.P. (1970) A Contribution to the Problem of the Structure of an Extended Neighborhood of a Rough Equilibrium State of Saddle-Focus Type. Mathematics of the USSR-Sbornik, 10, 91. https://doi.org/10.1070/SM1970v010n01ABEH001588  Shilnikov, L.P., Shilnikov, A.L., Turaev, D.V. and Chua, L.O. (2001) Methods Of Qualitative Theory In Nonlinear Dynamics (Part II). World Scientific, Singapore, New Jersey, London, Hong Kong. https://doi.org/10.1142/4221  Mira, C. (1987) Chaotic Dynamics: From the One-Dimensional Endomorphism to the Two-Dimensional Diffeomorphism. World Scientific, Singapore. https://doi.org/10.1142/0413  Ricker, W.E. (1954) Stock and Recruitment. Journal of the Fisheries Board of Canada, 11, 559-623. https://doi.org/10.1139/f54-039  Azizi, T., et al. (2015) Dynamics of a Discrete-Time Plant-Herbivore Model. Caspian Journal of Mathematical Sciences (CJMS), 4, 241-256.  Smith, H.L. and Thieme, H.R. (2011) Dynamical Systems and Population Persistence. American Mathematical Society, Providence, 118.  Kong, X.Z., Chen, L.P. and Yang, W.S. (2010) Note on the Persistent Property of a Discrete Lotka-Volterra Competitive System with Delays and Feedback Controls. Advances in Difference Equations, 2010, Article No.: 249364. https://doi.org/10.1155/2010/249364  Thieme, H.R. (1993) Persistence under Relaxed Point-Dissipativity (with Application to an Endemic Model). SIAM Journal on Mathematical Analysis, 24, 407-435. https://doi.org/10.1137/0524026  Hirsch, M.W., Smith, H.L. and Zhao, X.Q. (2001) Chain Transitivity, Attractivity, and Strong Repellors for Semidynamical Systems. Journal of Dynamics and Differential Equations, 13, 107-131. https://doi.org/10.1023/A:1009044515567  Kulenovic, M.R.S. and Ladas, G. (2001) Dynamics of Second Order Rational Difference Equations: With Open Problems and Conjectures. Chapman and Hall/CRC, New York. https://doi.org/10.1201/9781420035384 