# Analysis of a Single Species Model with Dissymmetric Bidirectional Impulsive Diffusion and Dispersal Delay

Authors: Haiyun Wan; Long Zhang; Zhidong Teng
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

In most models of population dynamics, diffusion between two patches is assumed to be either continuous or discrete, but in the real natural ecosystem, impulsive diffusion provides a more suitable manner to model the actual dispersal (or migration) behavior for many ecological species. In addition, the species not only requires some time to disperse or migrate among the patches but also has some possibility of loss during dispersal. In view of these facts, a single species model with dissymmetric bidirectional impulsive diffusion and dispersal delay is formulated. Criteria on the permanence and extinction of species are established. Furthermore, the realistic conditions for the existence, uniqueness, and the global stability of the positive periodic solution are obtained. Finally, numerical simulations and
discussion are presented to illustrate our theoretical results.


---

## Body

## 1. Introduction

In the last few years, mathematicians and ecologists have been actively investigating the dispersal of populations, a ubiquitous phenomenon in population dynamics. Levin [1] showed that both spatial dispersal of populations and population dynamics are much affected by spatial heterogeneity. In real life, dispersal often occurs among patches in ecological environments; because of the ecological effects of human activities and industries, such as the location of manufacturing industries and the pollution of the atmosphere, soil, and rivers, reproduction- and population-based territories and other habitats have been broken into patches. Thus, realistic models should include dispersal processes that take into consideration the effects of spatial heterogeneity.In recent years, increasing attention has been paid to the dynamics of a large number of mathematical models with diffusion, and many nice results have been obtained. The persistence and extinction for ordinary differential equation and delayed differential equation models were investigated in [2–6]. Global stability of equilibrium and periodic solution for diffusing model were studied in [7–12]. However, in all of above population dispersing systems, it is always assumed that the dispersal occurs at every time. For example, in [7], Beretta and Takeuchi proposed the following single-species diffusion Volterra models with continuous time delays:
(1)
x
˙
i
=
x
i
(
e
i
-
a
i
x
i
+
γ
i
∫
-
∞
t
‍
F
i
(
t
-
τ
)
x
i
(
τ
)
d
τ
)
+
∑
μ
=
1
n
‍
D
i
μ
(
t
)
(
x
μ
-
x
i
)
,
i
∈
N
,
where N
=
{
1
,
…
,
n
}
n is the number of patches, and x
i is the population density in the ith patch. The form of the dispersal established in this model is continuous; that is, the dispersal is always happening at any time.Actually, real dispersal behavior is very complicated and is always influenced by environmental change and human activities. In many practical situations, it is often the case that maybe one of the species suffers a significant loss or increase in density for some reason at some transitory time slots. These short-term perturbations are often assumed to be in the form of impulses in the modeling process. For example, when winter comes, birds will migrate between patches in search for a better environment, whereas they do not diffuse in other seasons, and the excursion of foliage seeds occurs during a fixed period of time every year. Therefore, impulsive differential equations [13] provide a natural description of such system. With the developments and applications of impulsive differential equations, theories of impulsive differential equations have been introduced into population dynamics, and many important studies have been performed [14–20].In [14], the authors studied the following autonomous single-species model with impulsively bidirectional diffusion:
(2)
x
˙
1
(
t
)
=
x
1
(
t
)
[
a
1
-
b
1
x
1
(
t
)
]
,
x
˙
2
(
t
)
=
x
2
(
t
)
[
a
2
-
b
2
x
2
(
t
)
]
,
t
≠
n
τ
,
Δ
x
1
(
t
)
=
d
1
(
x
2
(
t
)
-
x
1
(
t
)
)
,
Δ
x
2
(
t
)
=
d
2
(
x
1
(
t
)
-
x
2
(
t
)
)
,
t
=
n
τ
,
where a
i, b
i
(
i
=
1,2
) are the intrinsic growth rate and density-dependent parameters of the population x
i and d
i is the dispersal rate in the ith patch. Consider Δ
x
i
=
x
i
(
n
τ
+
)
-
x
i
(
n
τ
-
), where x
i
(
n
τ
+
)
=
lim
⁡
t
→
n
τ
+
x
i
(
t
) represents the density of population in the ith patch immediately after the nth diffusion pulse at time t
=
n
τ, while x
i
(
n
τ
-
)
=
lim
⁡
t
→
n
τ
-
x
i
(
t
)
=
x
i
(
τ
) represents the density of population in the ith patch before the nth diffusion pulse at time t
=
n
τ (τ the period of dispersal between any two pulse events is a positive constant, n
=
1,2
,
…). It is assumed here that the net exchange from the jth patch to ith patch is proportional to the difference x
j
-
x
i of population densities. The dispersal behavior of populations between two patches occurs only at the impulsive instants n
τ. Obviously, in this model, species x inhabits, respectively, two patches before the pulse appears, and when the time at the pulse comes, species x in two patches disperses from one patch to the other. The boundedness and global stability of positive periodic solution were obtained.Time delay often appears in many control systems (such as aircraft, chemical, or process control systems) either in the state, the control input, or the measurements. In order to reflect the dynamical behaviors of models that depend on the past history of system, it is often necessary to incorporate time delays into systems [21]. There have been extensive theoretical works on delay differential equations in the past three decades. The research topics include global asymptotic stability of equilibria, existence of periodic solutions, complicated behavior, and chaos (e.g., [8, 22–24]).Takeuchi et al. in [25] studied the following population model with time delays that introduced the dispersal time for individuals to move from one patch to other patches:
(3)
x
˙
i
(
t
)
=
x
i
(
t
)
[
a
i
-
b
i
x
i
(
t
)
]
+
∑
j
=
1
n
‍
[
ε
i
j
d
i
j
x
j
(
t
-
τ
i
j
(
t
)
)
-
d
j
i
(
t
)
x
i
(
t
)
]
,
hhhhhhhhhhhhhhhhhhhhh
i
=
1,2
,
…
,
n
,
where b
i, i
=
1
,
…
,
n, are positive constants, d
i
j, i, j
=
1
,
…
,
n, are nonnegative constants, ε
i
j
=
e
γ
i
j
τ
i
j satisfy 0
<
ε
i
j
≤
1, a
i, i
=
1
,
…
,
n, are constants, and some of them may be negative. In this paper, the authors took account of dispersal delay; however, they assumed that the dispersal is continuous.It is well known that the application of impulsive delay differential equations to population dynamics has been an interesting topic since it is reasonable and correct in modelling the evolution of population, such as pest management [26].However, in all of the impulsive dispersal models studied up till now, there are few papers considering the dispersal delay, which is really a pity. Actually, in the real world, the migration between patches is usually not immediate; that is, dispersal processes often involve time delay. For example, elks move from higher to lower elevations to escape cold in winter, and ungulates migrate annually among grazing areas to follow spatiotemporal changes in rainfall. Obviously, this kind of dispersal delay between patches extensively exists in the real world. Therefore, it is a very basilic problem to research this kind of population dynamic systems.Moreover, in the above impulsive dispersal models, it is assumed that the dispersal occurs between homogeneous habitat patches; that is, the dispersal rate between any two patches is equal or symmetrical [1, 11], which is really too idealized for a real ecosystem. Actually, in the real world, due to the heterogeneity of the spatiotemporal distributions in nature, movement between fragments of patches is usually not the same rate in both directions. In addition, once the individuals leave their present habitat, they may not successfully reach a new one, due to predation, harvesting, or other reasons, so that there are traveling losses. Thus, the dispersal rates among these patches are not always the same. Rather, in real ecological situations they are different (or dissymmetrical) [27, 28].Therefore, it is our basilic goal to investigate a single species model with dissymmetric impulse dispersal and dispersal delay. Motivated by the calculation hereinbefore, in this paper, we extend system (2) with dispersal delay and dispersal loss and consider it
(4)
d
x
1
(
t
)
d
t
=
x
1
(
t
)
[
a
1
-
b
1
x
1
(
t
)
]
,
d
x
2
(
t
)
d
t
=
x
2
(
t
)
[
a
2
-
b
2
x
2
(
t
)
]
,
t
≠
n
τ
,
Δ
x
1
(
t
)
=
d
2
x
2
(
t
-
τ
0
)
-
D
1
x
1
(
t
)
,
Δ
x
2
(
t
)
=
d
1
x
1
(
t
-
τ
0
)
-
D
2
x
2
(
t
)
,
t
=
n
τ
,
where D
i
(
i
=
1,2
) is the rate of population x
i emigrating from the ith patch and d
i is the rate of population x
i immigrating from the ith patch. Here we assume 0
≤
d
i
≤
D
i
≤
1, which means that there possibly exists mortality during migration between two patches. τ
0
≤
τ stands for the time delay; that is, a period of time of species x dispersing between patches.The organization of this paper is as follows. In Section2, as preliminaries, the definition of permanence and some useful lemmas are introduced. From discrete dynamic system theory, we establish the stroboscopic map of system (4), by which we can obtain the dynamical behaviors of it. In Section 3, the results of permanence and extinction for the system are presented. The existence and the uniqueness of the positive periodic solution for system (4) are established in Section 4. In Section 5, using the discrete dynamic system theory in [29], we can get the global stability of the positive periodic solution for the system. Finally, we give a brief discussion and our theoretical results are conformed by numerical simulations.


## 2. Preliminaries

In this section we introduce a definition and some notations and state some results which will be useful in subsequent sections.For any fixedt
∈
R
+
=
[
0
,
∞
), let ϕ
(
s
)
:
[
t
-
τ
0
,
t
]
→
R
2 be a piecewise continuous function such that ϕ
(
s
) is continuous in s
≠
n
τ, ϕ
(
n
τ
-
)
=
lim
⁡
s
→
n
τ
-
ϕ
(
s
)
=
ϕ
(
n
τ
) and ϕ
(
n
τ
+
)
=
lim
⁡
s
→
n
τ
+
ϕ
(
s
) exist. For any fixed t
∈
R
+, let P
C
[
t
-
τ
0
,
t
] denote the Banach space of all such piecewise continuous functions ϕ
(
θ
)
:
[
t
-
τ
0
,
t
]
→
R
2 with the norm ∥
ϕ
∥
=
sup
⁡
t
-
τ
0
≤
θ
≤
t
|
ϕ
(
θ
)
|. Further, let P
C
+
[
t
-
τ
0
,
t
] = {
ϕ
=
(
ϕ
1
,
ϕ
2
)
∈
P
C
[
t
-
τ
0
,
t
]
:
ϕ
i
(
θ
)
≥
0 for all θ
∈
[
t
-
τ
0
,
t
] and ϕ
i
(
t
)
>
0 for i
=
1,2
}.Motivated by the biological background of system (4), in this paper we always assume that all solutions of system (4) satisfy the following initial conditions:
(5)
x
i
(
θ
)
=
ϕ
i
(
θ
)
∀
θ
∈
[
-
τ
0
,
0
]
,
i
=
1
,
2
,
where ϕ
=
(
ϕ
1
,
ϕ
2
)
∈
P
C
+
[
-
τ
0
,
0
]. For any ϕ
∈
P
C
+
[
-
τ
0
,
0
], by the fundamental theory of impulsive functional differential equations [30, 31], system (4) has a unique solution x
(
t
,
ϕ
)
=
(
x
1
(
t
,
ϕ
)
,
x
2
(
t
,
ϕ
)
) satisfying the initial conditions (5).Definition 1.
System (4) is said to be permanent, if there are positive constants m
i and M
i such that
(6)
m
i
≤
liminf
⁡
t
→
∞
x
i
(
t
)
≤
limsup
⁡
t
→
∞
x
i
(
t
)
≤
M
i
,
i
=
1,2
,
for any positive solution x
(
t
)
=
(
x
1
(
t
)
,
x
2
(
t
)
) of system (4).We consider the following scalar impulsive differential equation:(7)
x
˙
(
t
)
=
f
(
t
,
x
(
t
)
)
,
t
≠
n
τ
,
x
(
n
τ
+
)
=
I
(
x
(
n
τ
)
)
,
t
=
n
τ
,
where t
∈
R
+, n
τ
>
0, is an impulsive time sequence, x
∈
R, f
(
t
,
x
)
:
R
+
×
R
→
R is continuous, and I
(
x
)
:
R
→
R is a nondecreasing function. We have the following comparison theorem [30, 31] for (7).Lemma 2.
Letx
(
t
) be a solution of system (7) defined on [
t
0
,
T
] and satisfy
(8)
u
˙
(
t
)
≤
(
≥
)
f
(
t
,
u
(
t
)
)
,
t
≠
n
τ
,
u
(
n
τ
+
)
≤
(
≥
)
I
(
u
(
n
τ
)
)
,
t
=
n
τ
.
Ifu
(
t
0
)
≤
(
≥
)
x
(
t
0
), then u
(
t
)
≤
(
≥
)
x
(
t
) for all t
∈
[
t
0
,
T
].Lemma 3 (see [32]).
Consider the following nonlinear impulsive system:(9)
x
˙
(
t
)
=
x
(
t
)
(
a
(
t
)
-
b
(
t
)
x
(
t
)
)
,
t
≠
t
k
,
x
(
t
k
+
)
=
(
1
-
θ
k
)
x
(
t
k
)
,
t
=
t
k
,
k
∈
N
,
x
(
0
+
)
=
x
0
,
where a
(
t
) and b
(
t
) are bounded and continuous ω-periodic functions defined on R
+, b
(
t
)
≥
0 for all t
∈
R
+ and impulsive coefficients θ
k
∈
(
0,1
)
,
θ
k
=
θ
k
+
q, q is a fixed positive integer, and t
k
+
q
=
t
k
+
ω. If ∫
0
ω
‍
a
(
s
)
d
s
+
ln
⁡
∏
i
=
1
q
‍
(
1
-
θ
i
)
>
0, then there exists a unique positive periodic solution of system (9), which is globally asymptotically stable.Remark 4.
If system (9) degenerates into the following autonomous impulsive differential equation:
(10)
x
˙
(
t
)
=
x
(
t
)
(
a
-
b
x
(
t
)
)
,
t
≠
t
k
,
x
(
t
k
+
)
=
(
1
-
θ
)
x
(
t
k
)
,
t
=
t
k
,
k
∈
N
,
x
(
0
+
)
=
x
0
,
where a, b, θ are positive constants, t
k
+
1
=
t
k
+
ω. As a consequence of Lemma 3, we have the following result: if 1
-
θ
-
e
-
a
ω
>
0, then system (10) has a unique positive periodic solution x
*
(
t
), which is globally asymptotically stable. In fact, here ∫
0
ω
‍
a
d
s
+
ln
⁡
(
1
-
θ
)
>
0
≡
1
-
θ
-
e
-
a
ω
>
0.Lemma 5 (see [29]).
LetF
:
R
+
n
→
R
+
n be continuous, C
1 in int (
R
+
n
), and suppose D
F
(
0
) exists with lim
⁡
x
→
0
+
D
F
(
z
)
=
D
F
(
0
). In addition, assume(a)
D
F
(
x
)
>
0,if
x
>
0;
(b)
D
F
(
y
)
<
D
F
(
x
),if
0
<
x
<
y.
IfF
(
0
)
=
0, let λ
=
ρ
(
D
F
(
0
)
). If λ
≤
1, then, for every x
≥
0
,
F
n
(
x
)
→
0 as n
→
∞; if λ
>
1, then either F
n
(
x
)
→
∞ as n
→
∞ for every x
>
0 or there exists a unique nonzero fixed point q of F. In the latter case, q
>
0 and, for every x
>
0, F
n
(
x
)
→
q as n
→
∞.
IfF
(
0
)
≠
0, then either F
n
(
x
)
→
∞ as n
→
∞ for every x
≥
0 or there exists a unique fixed point q of F. In the latter case, q
>
0 and, for every x
>
0, F
n
(
x
)
→
q as n
→
∞.Next, we analyze system (4). Integrating and solving the first two equations of system (4) between pulses, we have
(11)
x
i
(
t
)
=
[
b
i
a
i
+
(
1
x
i
(
n
τ
+
)
-
b
i
a
i
)
e
-
a
i
(
t
-
n
τ
)
]
-
1
,
h
h
h
h
h
h
h
h
h
n
τ
<
t
≤
(
n
+
1
)
τ
,
i
=
1,2
.Similarly, considering the last two equations of system (4), we obtain the following stroboscopic map:
(12)
x
1
n
+
1
=
(
1
-
D
1
)
x
1
n
h
1
+
c
1
x
1
n
+
d
2
x
2
n
h
2
e
a
2
τ
0
+
l
2
x
2
n
,
x
2
n
+
1
=
(
1
-
D
2
)
x
2
n
h
2
+
c
2
x
2
n
+
d
1
x
1
n
h
1
e
a
1
τ
0
+
l
1
x
1
n
;
here x
i
n
+
1
=
x
i
[
(
n
+
1
)
τ
+
], c
i
=
(
b
i
/
a
i
)
(
1
-
e
-
a
i
τ
)
>
0, l
i
=
(
b
i
/
a
i
)
(
1
-
e
-
a
i
(
τ
-
τ
0
)
)
>
0, 0
<
h
i
=
e
-
a
i
τ
<
1.Remark 6.
System (12) is a difference system, which means that densities of population in two patches have values at the previous pulse. We are, in other words, stroboscopically sampling at its pulsing period. The dynamical behavior of system (12), coupled with (11), determines the dynamical behavior of system (4). In the following sections, we will focus our attention on system (12) and investigate various aspects of its dynamical behavior.To write system (12) as a map, we define the map :
R
+
2
→
R
+
2:
(13)
F
1
(
x
1
,
x
2
)
=
(
1
-
D
1
)
x
1
h
1
+
c
1
x
1
+
d
2
x
2
h
2
e
a
2
τ
0
+
l
2
x
2
,
F
2
(
x
1
,
x
2
)
=
(
1
-
D
2
)
x
2
h
2
+
c
2
x
2
+
d
1
x
1
h
1
e
a
1
τ
0
+
l
1
x
1
.The set of all iterations of the mapF is equivalent to the set of all density sequences generated by system (12); F
(
x
) is the map evaluated at the point x
=
(
x
1
,
x
2
)
∈
R
+
2. Consequently, in system (12), F
n describes the population densities in the time n
τ.On the positivity of solutions of system (4) we have the following result.Lemma 7.
The solutionx
(
t
,
t
0
,
ϕ
) of system (4) with initial condition (5) is positive, that is, x
(
t
,
t
0
,
ϕ
)
>
0 on the interval of the existence.The proof of Lemma7 is simple; we hence omit it here.


## 3. Permanence and Extinction

In this section, we present conditions to ensure that system (12) is permanent and extinct which will imply the permanence and extinction of system (4). The permanence plays an important role in mathematical ecology since the criterion of permanence for ecological systems is a condition ensuring the long-term survival of all species. So, we firstly prove system (12) is permanent.Theorem 8.
Suppose(
H
1
)
1
-
D
i
-
h
i
>
0
,
i
=
1,2
,
hold; then system (12) is permanent.Proof.
Letx
i
(
t
)
∈
P
C
′ be the solution of system (4) satisfying the initial conditions (5). From the first equation of system (12), we have
(14)
x
1
n
+
1
=
1
-
D
1
h
1
(
x
1
n
)
-
1
+
c
1
+
d
2
h
2
e
a
2
τ
0
(
x
2
n
)
-
1
+
l
2
<
1
-
D
1
c
1
+
d
2
l
2
.
Similarly, we have(15)
x
2
n
+
1
=
1
-
D
2
h
2
(
x
2
n
)
-
1
+
c
2
+
d
1
h
1
e
a
1
τ
0
(
x
1
n
)
-
1
+
l
1
<
1
-
D
2
c
2
+
d
1
l
1
.
Hence, by (14) and (15) we know that system (12) has an ultimately upper bound.
Next, we prove that all the solutions of system (4) are ultimately below bounded. Since d
i
≥
0
(
i
=
1,2
), from the third equation of system (4), we have
(16)
x
1
(
n
τ
+
)
=
x
1
(
n
τ
-
)
+
Δ
x
1
(
n
τ
)
=
(
1
-
D
1
)
x
1
(
n
τ
)
+
d
2
x
2
(
n
τ
-
τ
0
)
≥
(
1
-
D
1
)
x
1
(
n
τ
)
.
Similarly,(17)
x
2
(
n
τ
+
)
≥
(
1
-
D
2
)
x
2
(
n
τ
)
.
Thus, system (4) becomes
(18)
x
˙
1
(
t
)
=
x
1
(
t
)
[
a
1
-
b
1
x
1
(
t
)
]
,
x
˙
2
(
t
)
=
x
2
(
t
)
[
a
2
-
b
2
x
2
(
t
)
]
,
t
≠
n
τ
,
x
1
(
n
τ
+
)
≥
(
1
-
D
1
)
x
1
(
n
τ
)
,
x
2
(
n
τ
+
)
≥
(
1
-
D
2
)
x
2
(
n
τ
)
,
t
=
n
τ
.
From (18), we find that there is no relation between x
1
(
t
) and x
2
(
t
). Therefore, we will discuss them, respectively;
(19)
x
˙
1
(
t
)
=
x
1
(
t
)
[
a
1
-
b
1
x
1
(
t
)
]
,
t
≠
n
τ
,
x
1
(
n
τ
+
)
≥
(
1
-
D
1
)
x
1
(
n
τ
)
,
t
=
n
τ
,
x
˙
2
(
t
)
=
x
2
(
t
)
[
a
2
-
b
2
x
2
(
t
)
]
,
t
≠
n
τ
,
x
2
(
n
τ
+
)
≥
(
1
-
D
2
)
x
2
(
n
τ
)
,
t
=
n
τ
.
If(
H
1
) holds, from Remark 4, we can obtain that the auxiliary system
(20)
u
˙
1
(
t
)
=
u
1
(
t
)
[
a
1
-
b
1
u
1
(
t
)
]
,
t
≠
n
τ
,
u
1
(
n
τ
+
)
=
(
1
-
D
1
)
u
1
(
n
τ
)
,
t
=
n
τ
,
has a unique positive periodic solution u
1
*
(
t
)
=
x
1
*
(
t
) which is globally asymptotically stable.
Letu
1
(
t
) be the solution of system (20) with initial value u
1
(
0
+
)
=
x
1
(
0
+
). By Lemma 2, we have
(21)
x
1
(
t
)
≥
u
1
(
t
)
,
∀
t
≥
0
.
Hence, for anyε
>
0 sufficiently small, there exists a T
1
>
0 such that
(22)
x
1
(
t
)
≥
x
1
*
(
t
)
-
ε
≜
m
1
,
for
t
≥
T
1
.
Similarly, if(
H
1
) holds, for above ε
>
0, there exists a T
2
>
0 such that
(23)
x
2
(
t
)
≥
x
2
*
(
t
)
-
ε
≜
m
2
,
for
t
≥
T
2
.
Denotem
=
min
⁡
{
m
1
,
m
2
} and T
=
max
⁡
{
T
1
,
T
2
}; then we have x
1
(
t
)
>
m and x
2
(
t
)
>
m, t
≥
T. Finally, we can determine that there exist constants α
i, β
i
(
0
<
α
i
<
β
i
)
(
i
=
1,2
), such that α
1
≤
liminf
⁡
n
→
∞
x
1
n
≤
limsup
⁡
n
→
∞
x
1
n
≤
β
1, and α
2
≤
liminf
⁡
n
→
∞
x
2
n
≤
limsup
⁡
n
→
∞
x
2
n
≤
β
2. The proof of Theorem 8 is completed.Next, we present condition to ensure that system (12) is extinct.Theorem 9.
System (12) is extinct if
(24)
[
d
1
d
2
e
(
a
1
+
a
2
)
τ
0
-
(
1
-
D
1
)
(
1
-
D
2
)
]
e
(
a
1
+
a
2
)
τ
+
e
a
1
τ
(
1
-
D
1
)
+
e
a
2
τ
(
1
-
D
2
)
≤
1
.Proof.
Let us consider the system (13). Obviously, F
(
x
1
,
x
2
) is continuous C
1 in int
⁡
(
R
+
2
), and F
(
0,0
)
=
0. We obtain
(25)
D
F
(
x
1
,
x
2
)
=
(
(
1
-
D
1
)
h
1
(
h
1
+
c
1
x
1
)
2
d
2
h
2
e
a
2
τ
0
(
h
2
e
a
2
τ
0
+
l
2
x
2
)
2
d
1
h
1
e
a
1
τ
0
(
h
1
e
a
1
τ
0
+
l
1
x
1
)
2
(
1
-
D
2
)
h
2
(
h
2
+
c
2
x
2
)
2
)
,
D
F
(
0,0
)
=
(
1
-
D
1
h
1
d
2
h
2
e
a
2
τ
0
d
1
h
1
e
a
1
τ
0
1
-
D
2
h
2
)
.
Obviously,lim
⁡
(
x
1
,
x
2
)
→
(
0,0
)
D
F
(
x
1
,
x
2
)
=
D
F
(
0,0
); if x
>
0
,
D
F
(
x
)
>
0; if 0
<
x
<
y, D
F
(
y
)
<
D
F
(
x
). We have the characteristic equation of D
F
(
0,0
):
(26)
λ
2
-
(
1
-
D
1
h
1
+
1
-
D
2
h
2
)
λ
+
(
1
-
D
1
)
(
1
-
D
2
)
h
1
h
2
-
d
1
d
2
h
1
h
2
e
(
a
1
+
a
2
)
τ
0
=
0
.
Letλ
=
ρ
(
D
F
(
0,0
)
); then we have
(27)
λ
=
(
+
4
d
1
d
2
h
1
h
2
e
(
a
1
+
a
2
)
τ
0
)
1
/
2
1
-
D
1
h
1
+
1
-
D
2
h
2
+
(
(
1
-
D
1
h
1
+
1
-
D
2
h
2
)
2
-
4
(
1
-
D
1
)
(
1
-
D
2
)
h
1
h
2
+
4
d
1
d
2
h
1
h
2
e
(
a
1
+
a
2
)
τ
0
(
1
-
D
1
h
1
+
1
-
D
2
h
2
)
2
)
1
/
2
)
×
(
2
)
-
1
=
(
+
(
1
-
D
1
h
1
-
1
-
D
2
h
2
)
2
+
4
d
1
d
2
h
1
h
2
e
(
a
1
+
a
2
)
τ
0
1
-
D
1
h
1
+
1
-
D
2
h
2
+
(
1
-
D
1
h
1
-
1
-
D
2
h
2
)
2
+
4
d
1
d
2
h
1
h
2
e
(
a
1
+
a
2
)
τ
0
)
×
(
2
)
-
1
.
Assumeλ
>
1; then by (27) we can obtain
(28)
1
-
D
1
h
1
+
1
-
D
2
h
2
+
d
1
d
2
h
1
h
2
e
(
a
1
+
a
2
)
τ
0
-
(
1
-
D
1
)
(
1
-
D
2
)
h
1
h
2
>
1
;
that is,
(29)
[
d
1
d
2
e
(
a
1
+
a
2
)
τ
0
-
(
1
-
D
1
)
(
1
-
D
2
)
]
e
(
a
1
+
a
2
)
τ
+
e
a
1
τ
(
1
-
D
1
)
+
e
a
2
τ
(
1
-
D
2
)
>
1
,
which contradicts with (24). Therefore we have λ
≤
1. By Lemma 5, we can get F
n
(
x
1
,
x
2
)
→
(
0,0
) as n
→
∞, which means that system (12) is extinct. This completes the proof.


## 4. Existence and Uniqueness of Positive Periodic Solution

In this part, we will prove the existence and uniqueness of the fixed points of system (12), which means that system (4) has a uniquely positive periodic solution.Theorem 10.
If(
H
1
) holds, then there exists a unique positive fixed point x
*
=
(
x
1
*
,
x
2
*
) of system (12).Proof.
Corresponding to (12), let us consider the following system:
(30)
x
1
=
(
1
-
D
1
)
x
1
h
1
+
c
1
x
1
+
d
2
x
2
h
2
e
a
2
τ
0
+
l
2
x
2
,
x
2
=
(
1
-
D
2
)
x
2
h
2
+
c
2
x
2
+
d
1
x
1
h
1
e
a
1
τ
0
+
l
1
x
1
.
From (30), we have
(31)
x
1
-
(
1
-
D
1
)
x
1
h
1
+
c
1
x
1
>
0
,
x
2
-
(
1
-
D
2
)
x
2
h
2
+
c
2
x
2
>
0
;
hence
(32)
x
1
>
1
c
1
(
1
-
D
1
-
h
1
)
=
ξ
,
x
2
>
1
c
2
(
1
-
D
2
-
h
2
)
=
η
.
From (30), we also obtain
(33)
x
2
=
h
2
e
a
2
τ
0
[
x
1
-
(
(
1
-
D
1
)
x
1
)
/
(
h
1
+
c
1
x
1
)
]
d
2
-
l
2
[
x
1
-
(
(
1
-
D
1
)
x
1
)
/
(
h
1
+
c
1
x
1
)
]
,
G
(
x
1
,
x
2
)
=
d
1
x
1
h
1
e
a
1
τ
0
+
l
1
x
1
+
(
1
-
D
2
)
x
2
h
2
+
c
2
x
2
-
x
2
.
Thus(34)
x
2
(
ξ
)
=
0
,
G
(
ξ
)
=
d
1
ξ
h
1
e
a
1
τ
0
+
l
1
ξ
>
0
.
Letx
2
→
+
∞; then f
(
x
1
)
=
x
1
-
(
(
1
-
D
1
)
x
1
)
/
(
h
1
+
c
1
x
1
)
→
d
2
/
l
2. And f
(
x
1
) is an increasing function on the interval [
ξ
,
+
∞
). Since f
(
ξ
)
=
0, so there exists x
-
1
>
ξ such that f
(
x
-
1
)
=
d
2
/
l
2. We can easily find that G
(
x
-
1
)
<
0. By the zero theory of continuous function, there exists (
x
1
*
,
x
2
*
) such that
(35)
ξ
<
x
1
*
<
x
-
1
,
G
(
x
1
*
,
x
2
*
)
=
0
,
(36)
x
2
*
=
h
2
e
a
2
τ
0
[
x
1
*
-
(
(
1
-
D
1
)
x
1
*
)
/
(
h
1
+
c
1
x
1
*
)
]
d
2
-
l
2
[
x
1
*
-
(
(
1
-
D
1
)
x
1
*
)
/
(
h
1
+
c
1
x
1
*
)
]
.
Next, we will prove the uniqueness of the fixed point.
It follows from (33) that we obtain
(37)
d
x
2
d
x
1
=
d
2
h
2
e
a
2
τ
0
[
1
-
h
1
(
1
-
D
1
)
/
(
h
1
+
c
1
x
1
)
2
]
{
d
2
-
l
2
[
x
1
-
(
(
1
-
D
1
)
x
1
)
/
(
h
1
+
c
1
x
1
)
]
}
2
,
d
G
d
x
1
=
d
1
h
1
e
a
1
τ
0
(
h
1
e
a
1
τ
0
+
l
1
x
1
)
2
+
[
h
2
(
1
-
D
2
)
(
h
2
+
c
2
x
2
)
2
-
1
]
d
x
2
d
x
1
,
d
G
d
x
1
=
1
×
(
{
d
2
-
l
2
[
x
1
-
(
1
-
D
1
)
x
1
h
1
+
c
1
x
1
]
}
2
×
(
h
1
e
a
1
τ
0
+
l
1
x
1
)
2
{
d
2
-
l
2
[
x
1
-
(
1
-
D
1
)
x
1
h
1
+
c
1
x
1
]
}
2
)
-
1
*
{
d
1
h
1
e
a
1
τ
0
{
d
2
-
l
2
[
x
1
-
(
1
-
D
1
)
x
1
h
1
+
c
1
x
1
]
}
2
+
d
2
h
2
e
a
2
τ
0
(
h
1
e
a
1
τ
0
+
l
1
x
1
)
2
×
[
1
-
h
1
(
1
-
D
1
)
(
h
1
+
c
1
x
1
)
2
]
[
h
2
(
1
-
D
2
)
(
h
2
+
c
2
x
2
)
2
-
1
]
}
.
Let(38)
φ
(
x
)
=
d
1
h
1
e
a
1
τ
0
{
d
2
-
l
2
[
x
1
-
(
1
-
D
1
)
x
1
h
1
+
c
1
x
1
]
}
2
+
d
2
h
2
e
a
2
τ
0
(
h
1
e
a
1
τ
0
+
l
1
x
1
)
2
×
[
1
-
h
1
(
1
-
D
1
)
(
h
1
+
c
1
x
1
)
2
]
[
h
2
(
1
-
D
2
)
(
h
2
+
c
2
x
2
)
2
-
1
]
;
then
(39)
d
φ
(
x
)
d
x
1
=
2
d
1
h
1
e
a
1
τ
0
{
d
2
-
l
2
[
x
1
-
(
1
-
D
1
)
x
1
h
1
+
c
1
x
1
]
}
×
{
-
l
2
[
1
-
h
1
(
1
-
D
1
)
(
h
1
+
c
1
x
1
)
2
]
}
+
2
l
1
d
2
h
2
e
a
2
τ
0
(
h
1
e
a
1
τ
0
+
l
1
x
1
)
×
[
1
-
h
1
(
1
-
D
1
)
(
h
1
+
c
1
x
1
)
2
]
[
h
2
(
1
-
D
2
)
(
h
2
+
c
2
x
2
)
2
-
1
]
+
d
2
h
2
e
a
2
τ
0
(
h
1
e
a
1
τ
0
+
l
1
x
1
)
2
×
2
c
1
h
1
(
1
-
D
1
)
(
h
1
+
c
1
x
1
)
3
[
h
2
(
1
-
D
2
)
(
h
2
+
c
2
x
2
)
2
-
1
]
+
d
2
h
2
e
a
2
τ
0
(
h
1
e
a
1
τ
0
+
l
1
x
1
)
2
×
[
1
-
h
1
(
1
-
D
1
)
(
h
1
+
c
1
x
1
)
2
]
[
-
2
c
2
h
2
(
1
-
D
2
)
(
h
2
+
c
2
x
2
)
3
]
.
By (31), we have (
1
-
D
1
)
/
(
h
1
+
c
1
x
1
)
<
1; since h
1
/
(
h
1
+
c
1
x
1
)
<
1, so h
1
(
1
-
D
1
)
/
(
h
1
+
c
1
x
1
)
2
<
1. Similarly, we have h
2
(
1
-
D
2
)
/
(
h
2
+
c
2
x
2
)
2
<
1. Therefore, we obtain d
φ
(
x
)
/
d
x
<
0, which implies that φ
(
x
) is a decreasing function on the interval [
ξ
,
+
∞
).
Since(40)
φ
(
ξ
)
=
d
1
d
2
2
h
1
e
a
1
τ
0
+
d
2
h
2
e
a
2
τ
0
×
[
h
1
e
a
1
τ
0
+
l
1
c
1
(
1
-
D
1
-
h
1
)
]
2
×
1
-
D
1
-
h
1
1
-
D
1
1
-
D
2
-
h
2
h
2
>
0
;
(41)
φ
(
x
-
1
)
=
d
1
h
1
e
a
1
τ
0
{
d
2
-
l
2
[
x
-
1
-
(
1
-
D
1
)
x
-
1
h
1
+
c
1
x
-
1
]
}
2
+
d
2
h
2
e
a
2
τ
0
(
h
1
e
a
1
τ
0
+
l
1
x
-
1
)
2
×
[
1
-
h
1
(
1
-
D
1
)
(
h
1
+
c
1
x
-
1
)
2
]
[
h
2
(
1
-
D
2
)
(
h
2
+
c
2
x
2
(
x
-
1
)
)
2
-
1
]
=
0
+
d
2
h
2
e
a
2
τ
0
(
h
1
e
a
1
τ
0
+
l
1
x
-
1
)
2
×
[
1
-
h
1
(
1
-
D
1
)
(
h
1
+
c
1
x
-
1
)
2
]
(
-
1
)
<
0
,
using the zero theory of continuous function, there exists a unique point ξ
1
∈
(
ξ
,
x
-
1
) such that φ
(
ξ
1
)
=
0. Besides,
(42)
φ
(
x
1
)
>
0
,
∀
x
1
∈
(
ξ
,
ξ
1
)
,
φ
(
x
1
)
<
0
,
∀
x
1
∈
(
ξ
1
,
+
∞
)
;
thus
(43)
d
G
(
x
1
)
d
x
1
>
0
,
∀
x
1
∈
(
ξ
,
ξ
1
)
,
d
G
(
x
1
)
d
x
1
<
0
,
∀
x
1
∈
(
ξ
1
,
+
∞
)
,
which, together with G
(
ξ
)
>
0, leads to G
(
x
1
)
>
0
,
∀
x
1
∈
(
ξ
,
ξ
1
). By G
(
ξ
1
)
>
0
,
G
(
x
-
1
)
<
0, we have that there exists a unique point x
1
*
∈
(
ξ
1
,
x
-
1
) such that G
(
x
1
*
,
x
2
*
)
=
0. The proof is completed.


## 5. Global Stability

Now, we prove that the positive fixed points(
x
1
*
,
x
2
*
) of system (30) are globally stable by using Lemma 5, which means that the positive periodic solution of system (4) is globally stable.Theorem 11.
If(
H
1
) holds, then there exists a unique positive fixed point x
*
=
(
x
1
*
,
x
2
*
) of the map F, and, for every x
=
(
x
1
,
x
2
)
>
0
,
F
n
(
x
)
→
x
* as n
→
∞.Proof.
For any smallε
1
>
0, ε
2
>
0, we make the change of variable
(44)
x
1
=
u
+
ε
1
,
x
2
=
v
+
ε
2
.
By (13), we get the map F
(
u
,
v
)
=
(
F
1
(
u
,
v
)
,
F
2
(
u
,
v
)
); that is,
(45)
u
=
(
1
-
D
1
)
(
u
+
ε
1
)
h
1
+
c
1
(
u
+
ε
1
)
+
d
2
(
v
+
ε
2
)
h
2
e
a
2
τ
0
+
l
2
(
v
+
ε
2
)
-
ε
1
=
F
1
(
u
,
v
)
,
v
=
(
1
-
D
2
)
(
v
+
ε
2
)
h
2
+
c
2
(
v
+
ε
2
)
+
d
1
(
u
+
ε
1
)
h
1
e
a
1
τ
0
+
l
1
(
u
+
ε
1
)
-
ε
2
=
F
2
(
u
,
v
)
.
Now, we show thatF
(
u
,
v
) satisfies the hypotheses of Lemma 5. It is easy to see that F
(
u
,
v
) is continuous, C
1 in int
⁡
(
R
+
2
), and F
(
0,0
)
≠
0.
Since(46)
D
F
(
u
,
v
)
=
(
(
1
-
D
1
)
h
1
[
h
1
+
c
1
(
u
+
ε
1
)
]
2
d
2
h
2
e
a
2
τ
0
[
h
2
e
a
2
τ
0
+
l
2
(
v
+
ε
2
)
]
2
d
1
h
1
e
a
1
τ
0
[
h
1
e
a
1
τ
0
+
l
1
(
u
+
ε
1
)
]
2
(
1
-
D
2
)
h
2
[
h
2
+
c
2
(
v
+
ε
2
)
]
2
)
,



D
F
(
0,0
)
=
(
(
1
-
D
1
)
h
1
(
h
1
+
c
1
ε
1
)
2
d
2
h
2
e
a
2
τ
0
(
h
2
e
a
2
τ
0
+
l
2
ε
2
)
2
d
1
h
1
e
a
1
τ
0
(
h
1
e
a
1
τ
0
+
l
1
ε
1
)
2
(
1
-
D
2
)
h
2
(
h
2
+
c
2
ε
2
)
2
)
,
so, lim
⁡
(
u
,
v
)
→
(
0
+
,
0
+
)
D
F
(
u
,
v
)
=
D
F
(
0,0
). Obviously, if (
u
,
v
)
>
0, then D
F
(
u
,
v
)
>
0; if (
0,0
)
<
(
u
1
,
v
1
)
<
(
u
2
,
v
2
), then D
F
(
u
1
,
v
1
)
>
D
F
(
u
2
,
v
2
). It satisfies all the conditions of Lemma 5; then, for every u
>
0
,
v
>
0, we have F
n
(
u
,
v
)
→
(
x
1
*
-
ε
1
,
x
2
*
-
ε
2
) as n
→
∞. Corresponding to x
-
y coordinate, this means, for x
1
>
ε
1, x
2
>
ε
2, the system (30) tends to the unique fixed point.
It follows from the permanence of system (12) that we have x
1
n
>
ε
1, x
2
n
>
ε
2 for any initial value (
x
1
(
0
+
)
,
x
2
(
0
+
)
)
>
(
0,0
).
From the above analysis, we can know that, for everyx
1
(
0
+
)
>
0, x
2
(
0
+
)
>
0, the trajectory of system (12) will tend to (
x
1
*
,
x
2
*
). The proof is completed.


## 6. Numerical Simulation and Discussion

In order to test the validity of our results, first, for (4) we use the parameters values (Val. 1) in Table 1. We can easily test that the assumptions in Theorems 8 and 10 hold, which means the populations x
(
t
)
=
(
x
1
(
t
)
,
x
2
(
t
)
) in the two patches are permanent and have a unique periodic solution x
*
(
t
)
=
(
x
1
*
(
t
)
,
x
2
*
(
t
)
) which is globally stable (see Figure 1(a)). Moreover, if, in Table 1, we consider the influence of time delay, then we can see that the permanence and stability for species x unchanged. The details are given in Table 2. However, the longer the duration of the time delay, the lower the limit inferior and the limit superior of x (see Figures 1(b), 1(c), and 1(d)). This implies that the case with dispersal delay is harmful to live for species x.Table 1
Parameter values used in the simulations of Model (4).
Parameter
Interpretation
Val. 1
Val. 2
a
1
Intrinsic growth rate of populations in patch 1
0.9
0.3
a
2
Intrinsic growth rate of populations in patch 2
0.8
0.4
b
1
Density dependence rate of populations in patch 1
0.75
0.1
b
2
Density dependence rate of populations in patch 2
0.5
0.2
D
1
Successfully emigrant rate of populations from patch 1
0.8
0.9
D
2
Successfully emigrant rate of populations from patch 2
0.7
0.8
d
1
Successfully immigrant rate from patch 1 to patch 2
0.6
0.7
d
2
Successfully immigrant rate from patch 2 to patch 1
0.1
0.5
τ
The period of dispersal between two pulse events
2
2
τ
0
Dispersal delay between two pulse events
0
1
Val: value.Table 2
Simulations of Model (4).
Case
τ
0
x
1
x
2
Figure
1
0
Permanent
Permanent
Figure1(a)
2
1
Permanent
Permanent
Figure1(b)
3
1.5
Permanent
Permanent
Figure1(c)
4
2
Permanent
Permanent
Figure1(d)(a), (b), (c), and (d) Dynamical behavior of system (4). Here, we take three sets of initial values (1, 1.2), (0.6, 0.8), and (0.2, 1).
(a)
(b)
(c)
(d)Next, we take the parameters values (Val. 2) in Table1. We can easily test that assumption in Theorem 9
[
(
d
1
d
2
/
e
(
a
1
+
a
2
)
τ
0
)
-
(
1
-
D
1
)
(
1
-
D
2
)
]
e
(
a
1
+
a
2
)
τ
+
e
a
1
τ
(
1
-
D
1
)
+
e
a
2
τ
(
1
-
D
2
)
=
0.5865
<
1 holds, which implies the populations x
(
t
)
=
(
x
1
(
t
)
,
x
2
(
t
)
) in the two patches are extinct (see Figure 2). Comparing Figure 2(a) and Figure 2(b), we realize that system (4) with time delay accelerates the extinction comparing with no delay (see Figure 2). This is reasonable from a biological point of view. Without delay means less loss during dispersion, which implies that more members can arrive to other patches. Otherwise, populations go extinct due to much loss. The details are given in Table 3.Table 3
Simulations of Model (4).
Case
τ
0
x
1
x
2
Figure
1
0
Extinct
Extinct
Figure2(a)
2
1
Extinct
Extinct
Figure2(b)(a) and (b) Dynamical behavior of system (4). Here, we take three sets of initial values (1.7, 2.1), (1.9, 1.5), and (1.4, 2.2).
(a)
(b)Furthermore, we takea
1
=
1.1
+
0.04
cos
⁡
(
π
t
), a
2
=
0.8
+
0.01
sin
(
π
t
), b
1
=
0.5
+
0.04
sin
(
π
t
), b
2
=
0.4
+
0.12
cos
⁡
(
π
t
), d
2
=
0.55 and keep other parameters unchanged from in Val. 1. Here we can see that the period of individual intrinsic growth rate and density dependence rate T
=
2 is equal to the period of migration τ. In this case, we have let the period of the environment match the period of migration. If we take τ
0
=
0, from simulation (see Figure 3(a)), we can see the populations x
(
t
)
=
(
x
1
(
t
)
,
x
2
(
t
)
) are permanent and have a uniquely periodic solution which is globally stable. However, if we take the migration period τ
=
2.009
,
2.022, respectively, with other parameters unchanged, from numerical simulations (see Figures 3(b) and 3(c)), we can see that population dynamics change from almost periodic to chaotic. For Case 4 of Table 4, we consider the influence of time delay (see Figure 3(d)). Comparing Figures 3(c) and 3(d), we realize that system (4) with time delay is more complicated than without. The details are given in Table 4.Table 4
Simulations of Model (4).
Case
τ
τ
0
x
1
x
2
Figure
1
2
0
PP
PP
Figure3(a)
2
2.009
0
PAP
PAP
Figure3(b)
3
2.022
0
PC
PC
Figure3(c)
4
2.05
1.1
PC
PC
Figure3(d)
PP: permanent and periodic, PAP: permanent and almost periodic, and PC: permanent and chaotic.(a), (b), (c), and (d) Dynamical behavior of system (4). Here, we take the initial value x
0
=
(
x
10
,
x
20
) = (1.735, 1.873).
(a)
(b)
(c)
(d)Lastly, if we takea
1
=
0.9
+
0.04
cos
⁡
(
π
t
) and keep other parameters unchanged with Figure 3, by numerical simulations (see Figure 4), we find that all of the solutions of system (4) which through the initial points will converge to the positive periodic solution (
x
1
*
,
x
2
*
). Therefore, we can guess that under the assumptions of Theorem 10 system (4) has a unique positive periodic solution which is globally asymptotically stable. In addition, the periodic solution with time delay is larger than without delay which indicates that the duration of the time delay is beneficial to live for species x
2 compared with species x
1.(a) and (b) Dynamical behavior of system (4). We take a series of initial points, such as (1.935,2.3), (1.939,2.34), and (1.943,2.38).
(a)
(b)



---

*Source: 701545-2014-04-09.xml*