# Theoretical studies on the effects of dispersal corridors on the permanence of discrete predator-prey models in patchy environment.

1. IntroductionThe conception of permanence in ecological communities is first discussed in [1]. This notion of permanence is also known under different aliases, such as permanent coexistence or uniform persistence. Aptly termed, permanence of a population model reflects that all the species in the ecosystem will coexist for a relatively long time; that is, no extinction of any species will occur in the near future--this is very important in biological conservation. A permanent model may attain globally asymptotically stable positive equilibrium [2] or may have periodic solution [3, 4], and even a chaotic model can be permanent [5]. The sufficient and/or necessary conditions established to guarantee the permanence of various population models may offer helpful suggestions for environmental managers to take applicable measures for the sustainability of the ecosystems [6-8].

It is common that in the real world today the habitats of the species are separated into isolated patches due to spatial heterogeneity or the development of mankind, such as the construction of highways or railways [9-16]. This patchy environment may restrict the activities of the species and lead to adverse effects on the permanence of the ecological communities; therefore, some dispersal corridors maybe built to allow movement between the patches. For example, the construction of Qinghai-Tibet railway in China separates the habitats of the Tibetan antelopes to pieces of small patches. In order to reduce the impact of patchy environment on the Tibetan antelopes, several underpasses (dispersal corridors) along the railway are constructed for the Tibetan antelopes to move from one side of the railway to the other side [17]. Indeed some researches [18] have verified that the Tibetan antelopes have really made use of these underpasses to cross the railway.

In the literature there are quite many papers that discuss the permanence of continuous time predator-prey or competitive systems in patchy environment: for example, one can refer to [19-23] and the many references cited therein. Though some continuous time models have their discrete time analogues, it is common knowledge that these discrete time analogues may present different dynamics, which are possibly more complicated [24-26]. To our best knowledge, there are few references on discrete time predator-prey or competitive models in patchy environment; see, for instance, [22, 27-30]. In [28, 30], the effects of dispersion on discrete models for single species in patchy environment have been discussed; two patches have been considered in [30] while n patches have been tackled in [28]. In [27, 29], the effects of dispersion on discrete competitive models in patchy environment have been analyzed; in particular discrete Leslie-Gower type competitive model has been adopted in [27]. Motivated by the above papers as well as the importance of permanence in ecosystems, in this paper we will investigate the permanence of discrete time predator-prey models of Ricker type in patchy environment. To be specific, we will tackle two discrete predator-prey models in patchy environment; the first model has no dispersal corridors while the second model has dispersal corridors. Sufficient conditions to guarantee permanence will be established for both models. The effects of dispersal corridors on permanence are then analyzed by comparing the two sets of sufficient conditions. We remark that in the literature the numerous research papers on continuous/discrete predator-prey models in patchy environment have focused mainly on the permanence criteria and/or the existence of periodic solutions of the models; however, the effects of dispersal corridors on the permanence of the models are rarely studied until the present paper.

The paper is organized as follows. In Section 2, we construct the two discrete predator-prey models in patchy environment with/without dispersal corridors; also in this section we give some preliminaries such as the definition of permanence. In Section 3, sufficient conditions for the permanence of the two discrete predator-prey models in patchy environment with/without dispersal corridors are offered. The effects of dispersal corridors on the permanence of the models are analyzed in Section 4. In Section 5, we present five examples where numerical simulations are performed to confirm the theoretical results and to illustrate the effects of dispersal corridors. Finally, we conclude in Section 6. The mathematical proofs of the theorems in Section 3 are presented in Appendix.

2. Models and Preliminaries

There are many types of discrete time mathematical models, which are governed by difference equations, established to reflect the fluctuation of the single-species population [31]. Among them, the Ricker type model, given as

x(n + 1) = x (n) exp [r - ax (n)], (1)

is first established for some kind of fishery population [32] and subsequently used to model other kinds of single-species population [31, 33, 34]. In model (1), x(n) is the population density of the species x at the nth time step, r is the intrinsic growth rate of x, a is the strength of intraspecific competition of x, and r, a are positive constants. Model (1) is logistic; that is, the fluctuation of the population of x is undergoing the effect of density dependency. One can refer to [31] for more details about (1).

In the predator-prey model that we are considering in this paper, the prey is x and the predator is y, and the predator y only feeds on the prey x. Adopting model (1) to reflect the population growth of the prey x, we have the following discrete predator-prey model of Lotka-Volterra type:

x(n + 1) = x(n) exp [[r.sub.1] - [a.sub.11]x(n) - [a.sub.12]y(n)],

y(n + 1) = y(n) exp [-[r.sub.2] + [a.sub.21]x(n) - [a.sub.22]y(n)]. (2)

In model (2), [r.sub.1] is the intrinsic growth rate of the prey x, [a.sub.11] is the intraspecific competition strength of the prey x, [a.sub.12] is the predation strength, [r.sub.2] is the death rate of the predator y if there is no prey x for y to feed on, [a.sub.21] reflects the predation efficiency or the rate of energy transferred from the prey x to the predator y through predation, and [a.sub.22] is the intraspecific competition strength of the predator y. The parameters [r.sub.i], [a.sub.ij], i, and j [member of] {1, 2} are all positive constants. The permanence and the global asymptotic stability of the positive equilibrium of model (2) have been discussed in [35, 36].

Now, suppose that the habitats of the predator-prey system (2) are separated into two patches, patch 1 and patch 2. We denote [x.sub.i](n) and [y.sub.i](n) as the respective population density of the prey x and the predator y in patch i at the nth time step, i [member of] {1,2}. Let [r.sub.1i] be the intrinsic growth rate of the prey x in patch i and let [r.sub.2i] be the death rate of the predator y in patch i, i [member of] {1, 2}. Also, [a.sup.(k).sub.ij] has the same biological meaning as [a.sub.ij] in model (2) except that now it is the case in patch k, i, j, k [member of] {1,2}. At this moment, there are no dispersal corridors between the two patches. Thus, we have the following discrete predator-prey model in patchy environment without dispersal corridors:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (3)

Next, we introduce dispersal corridors to the patchy model (3). It is assumed that the prey x can migrate from one patch to another via these dispersal corridors, but the predator y cannot use the dispersal corridors to migrate. We make some remarks to this assumption. Firstly, this assumption of only prey (but not predator) dispersal is usual in continuous/discrete time predator-prey models in patchy environment; see, for example, [22, 23, 37, 38]. Moreover, this assumption fits the common definition of prey refuge if the predator is fixed in one patch, and intensive studies have been done; refer to [39] and the references cited therein. Our assumption here is that the prey can move among patches while the predator is fixed in each patch--an example of such a system can be found in [40], where the prey (mayflies) is mobile among different patches, but the predator (charr) is confined to its patch since the river is fenced into patches. Indeed, systems of this type may exist naturally in swamp fields, where the prey (such as flying insects) can migrate to different ponds while the predator (such as fish) is fixed in each pond. Denote [d.sub.ij] as the dispersal rate of the prey x migrating from patch i to patch j, i, j [member of] {1,2}, i [not equal to] j. We have the following discrete predator-prey model in patchy environment with dispersal corridors:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (4)

Here, [alpha] and [beta] represent the predation rates of the predator on the prey that migrates from another patch. Note that model (4) is reduced to model (3) if the dispersal rates [d.sub.12] and [d.sub.21] are both zero. Model (4) may be viewed as a discrete analogue of the continuous time predator-prey reaction-diffusion model for patchy environment proposed by Allen [41]. In [41] an extinction result, whereby some sufficient conditions are established to ensure the extinction of the species, has been developed. We remark that, from a practical standpoint, permanence results that guarantee species survival are more useful than extinction results; as such, we will focus on the permanence of the discrete model.

In view of the biological background, we let the initial values of models (3) and (4) be

[x.sub.1](0) = [x.sub.1] > 0, [x.sub.2](0) = [x.sub.2] > 0;

[y.sub.1](0) = [y.sub.1] > 0, [y.sub.2](0) = [y.sub.2] > 0. (5)

Also, throughout the paper we assume the following:

([H.sub.1]) [r.sub.ij], [a.sup.(k).sub.ij], i, j, k [member of] {1, 2} and [alpha], [beta] are all positive constants, and [d.sub.12], [d.sub.21] are nonnegative constants.

It is obvious that all the solutions of either model (3) or model (4) with initial values (5) are positive, which corresponds well to the biological meanings of [x.sub.i](n), [y.sub.i](n), i [member of] {1,2}, noting that population density of any species is always nonnegative.

The permanence results of models (3) and (4) will be presented in Section 3. To be precise, we state the definition of permanence as follows.

Definition 1. Model (3) (or 4)) with initial values (5) is said to be permanent if there exist positive constants [m.sub.i], [m.sup.*.sub.i], [M.sub.i], and [M.sup.*.sub.i], i [member of] {1, 2}, such that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (6)

hold for each solution of model (3) (or (4)) with initial values (5).

3. Permanence of the Models

In this section, we will establish the permanence of models (3) and (4) with initial values (5). The mathematical proofs of Theorems 2, 4, and 5 will be presented in Appendix.

3.1. Permanence without Dispersal Corridors. We consider model (3) together with initial values (5). Notice that in (3) the first two equations are independent of the last two equations; therefore, we only need to consider the permanence of the model comprising just the first two equations of (3); namely,

[x.sub.1](n + 1) = [x.sub.1](n)exp[[r.sub.11] - [a.sup.(1).sub.11][x.sub.1](n) - [a.sup.(1).sub.12][y.sub.1](n)],

[y.sub.1](n + 1) = [y.sub.1](n)exp[-[r.sub.21] + [a.sup.(1).sub.21][x.sub.1](n) - [a.sup.(1).sub.22][y.sub.1](n)] (7)

with initial values

[x.sub.1](0) = [x.sub.1] > 0, [y.sub.1](0) = [y.sub.1] > 0. (8)

Theorem 2. Suppose that

[r.sub.11] - [a.sup.(1).sub.12][M.sub.1] > 0, (9)

-[r.sub.21] + [a.sup.(1).sub.21][m.sub.1] > 0 (10)

are satisfied, where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (11)

Then, model (7) with initial values (8) is permanent.

Remark 3. (a) The permanence or global asymptotical stability for the positive equilibrium of model (7) has been studied in [35, 36, 42]. In [35] it is shown that, for sufficiently small [r.sub.11] and [r.sub.21], model (7) has a positive equilibrium which is globally attractive, and as a consequence model (7) is permanent. So the conditions for the permanence of (7) are not directly derived in [35]. On the other hand, based on the study of a nonautonomous version of (7), it is proved in [36] (Corollary 4) that the positive equilibrium of (7) is globally asymptotically stable if [r.sub.11] [less than or equal to] 1 together with other conditions. Further, in [42] it is proved that the positive equilibrium of (7) can also be globally asymptotically stable if [r.sub.11] > 1. It is known that a system is permanent if its positive equilibrium is globally asymptotically stable [8]. Hence, in [36, 42] once again permanence is a consequence but the conditions to ensure it are indirectly obtained.

(b) Our method to obtain the permanence of (7) (or (3)) is different from that in [35, 36, 42]. In particular, the inequalities [m.sub.1] [less than or equal to] [M.sub.1] and [m.sup.*.sub.1] [less than or equal to] [M.sup.*.sub.1] in Definition 1 can be verified here (refer to (A.17), (A.18), (A.26), and (A.27) in Appendix). In most cases if it is proved that each solution of a model is bounded both below and above by some positive constants, then the model is considered to be permanent. However, it is not easy to show the inequalities such as (A.17) and (A.26).

Using Theorem 2, we give the following theorem on the permanence of model (3).

Theorem 4. Suppose that conditions (9) and (10) are satisfied, and

[r.sub.12] - [a.sup.(2).sub.12][M.sub.2] > 0, (12)

-[r.sub.22] + [a.sup.(2).sub.21][m.sub.2] > 0 (13)

hold, where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)

Then, model (3) with initial values (5) is permanent.

3.2. Permanence with Dispersal Corridors. The following theorem gives the sufficient conditions to guarantee the permanence of predator-prey model (4) in patchy environment with dispersal corridors.

Theorem 5. Suppose that

[r.sub.11] - [a.sup.(1).sub.12][S.sub.1] > 0, (15)

-[r.sub.21] + [a.sup.(1).sub.21][s.sub.1] > 0, (16)

[r.sub.12] - [a.sup.(2).sub.12][S.sub.2] > 0, (17)

-[r.sub.22] + [a.sup.(2).sub.21][s.sub.2] > 0 (18)

are satisfied, where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (19)

Then, model (4) with initial values (5) is permanent.

Remark 6. Unlike the proof of Theorem 4, in the proof of Theorem 5 (refer to Appendix) the inequalities [m.sub.1] [less than or equal to] [M.sub.1], [m.sup.*.sub.1] [less than or equal to] [M.sup.*.sub.1], [m.sub.2] [less than or equal to] [M.sub.2], and [m.sup.*.sub.2] [less than or equal to] [M.sup.*.sub.2] in Definition 1 cannot be verified. This also reinforces Remark 3(b) where we indicate that such inequalities are not easy to prove in most cases.

4. Effects of Dispersal Corridors

In Section 3, sufficient conditions for the permanence of models (3) and (4) with initial values (5) are offered in Theorems 4 and 5, respectively. In this section, we will analyze the effects of dispersal corridors on the permanence of discrete predator-prey system by comparing the two sets of sufficient conditions obtained. Recall that the sufficient conditions to guarantee the permanence of predator-prey system without dispersal corridors in patchy environment are (9)-(13), while (15)-(18) are the sufficient conditions to guarantee the permanence of the similar system with dispersal corridors. From the proofs of Theorems 4 and 5 (see Appendix), these conditions can be classified in Table 1.

4.1. Effects of Dispersal Corridors on the Permanence of the Prey. We begin with conditions (9) and (15), which are mainly concerned with the permanence of the prey x in patch 1 under the situation of without or with dispersal corridors. Rewrite (9) as

[r.sub.11] > [a.sup.(1).sub.12][M.sub.1]

= [[a.sup.(1).sub.12]/[a.sup.(1).sub.22]]exp[-[r.sub.21] + [[a.sup.(1).sub.21]/[a.sup.(1).sub.11]]exp([r.sub.11] - 1) - 1] [??] [B.sub.1]. (20)

Also, rewrite (15) as

[r.sub.11] > [a.sup.(1).sub.12][S.sub.1]

= [[a.sup.(1).sub.12]/[a.sup.(1).sub.22]]exp[-[r.sub.21] + [[a.sup.(1).sub.21]/[[a.sup.(1).sub.11] + [d.sub.12]]]exp([r.sub.11] - 1) - 1] [??] [B.sub.2]. (21)

On the left sides of (20) and (21), [r.sub.11] is the intrinsic growth rate of the prey x in patch 1, which makes positive contribution to the population growth of the Prey since it is greater than 0. The right sides of (20) and (21) can be understood as the negative contribution, caused by the intracompetition of species x in patch 1, the predation by the predator y in patch 1, and also the dispersion to patch 2 if dispersal corridors exist, to the population growth of the prey. Thus, ([r.sub.11] - [B.sub.1]) and ([r.sub.11] - [B.sub.2]) can be interpreted as the net reproduction rates of the prey x in patch 1 without or with dispersal corridors, respectively. Consequently, the biological interpretation of (20) or (21) is that the prey x in patch 1 is permanent if its net reproduction rate is positive.

The only difference between (20) and (21) is whether the dispersal rate [d.sub.12] is involved or not. Recall that [d.sub.12] (>0) is the dispersion rate of the prey x from patch 1 to patch 2. It is obvious that

[B.sub.1] > [B.sub.2], (22)

or equivalently

[r.sub.11] - [B.sub.1] < [r.sub.11] - [B.sub.2]. (23)

This inequality shows that the net reproduction rate of the prey in patch 1 is larger when dispersal corridors are introduced. Hence, the existence of dispersal corridors for the prey x to migrate from patch 1 to patch 2 is helpful to the permanence of the prey in patch 1. Moreover, if [r.sub.11] - [B.sub.1] > 0, then [r.sub.11] - [B.sub.2] > 0 always holds; but if [r.sub.11] - [B.sub.2] > 0, then we are not certain whether [r.sub.11] - [B.sub.1] > 0. In other words, with the intrinsic growth rate [r.sub.11] of the prey x in patch 1 keeping the same value regardless of the existence or nonexistence of dispersal corridors, if the prey in patch 1 can be permanent without dispersal corridors, then it is also permanent in patch 1 with dispersal corridors; but if the prey in patch 1 is permanent with dispersal corridors, it may or may not be permanent without the dispersal corridors. In fact, the population growth of the prey in patch 1 is density dependent; hence, an appropriate dispersion from patch 1 to patch 2 may counterbalance the effect of density dependency. Certainly only the dispersion rate [d.sub.12] in an appropriate range is helpful to the permanence of the prey in patch 1. If the dispersion rate is too large, then the residence of the prey in patch 1 will decrease to a very small number, which will lead to a decrease in the intrinsic growth rate [r.sub.11] and will consequently increase the value of [a.sup.(1).sub.22] the intraspecific competition strength of the predator in patch 1, due to limited food sources. Hence, in the case of large dispersion rate [d.sub.12], the value of [B.sub.2] will decrease. If the decrease in [r.sub.11] is larger than the decrease in [B.sub.2], then the net reproduction rate of the prey in patch 1 may decrease and (21) may not hold; that is, the permanence of the prey in patch 1 may not be guaranteed.

In summary, an appropriate degree of dispersion of the prey from patch 1 to patch 2, measured by the dispersion rate [d.sub.12], is helpful to the permanence of the prey in patch 1. Similar conclusion can be drawn in patch 2 by comparing the condition (12) with (17). An appropriate dispersion rate [d.sub.21] of the prey from patch 2 to patch 1 is helpful to the permanence of the prey in patch 2.

4.2. Effects of Dispersal Corridors on the Permanence of the Predator. We begin with conditions (10) and 16), which are mainly concerned with the permanence of the predator y in patch 1 under the situation of without or with dispersal corridors. Rewrite (10) as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (24)

and rewrite (16) as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (25)

For convenience, we also list the expressions of [M.sub.1] and [S.sub.1] (refer to Theorems 2 and 5) as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (26)

On the left sides of (24) and (25), [r.sub.21] is the death rate of the predator y in patch 1. The right sides of (24) and (25) can be understood as the positive contribution to the population growth of the predator in patch 1 due to the predation. Hence, ([B.sub.3] - [r.sub.21]) and ([B.sub.4] - [r.sub.21]) can be interpreted as the net reproduction rates of the predator y in patch 1 without or with dispersal corridors, respectively. Consequently, the biological interpretation of (24) or (25) is that the predator y in patch 1 is permanent if its net reproduction rate is positive.

There are two "places" in (25) where the dispersal rate [d.sub.12] appears; one is in the denominator of [B.sub.4], while another is in the term [S.sub.1]. We will discuss the biological implications of them. First, it is clear that the [d.sub.12] that is explicit in the denominator of (25) has a diminishing effect on [B.sub.4]. Therefore, this dispersal rate in the denominator is disadvantageous to the net reproduction rate of the population growth, further the permanence, of the predator in patch 1. In fact, the migration of the prey from patch 1 to patch 2, while the predator remains fixed to patch 1, causes the decrease in the food sources of the predator, which in turn is disadvantageous to the permanence of the predator in patch 1.

Next, we consider the effect of the dispersal rate [d.sub.12] in the term [S.sub.1]. It is obvious from 26) that [S.sub.1] < [M.sub.1]; thus

[r.sub.11] - [a.sup.(1).sub.12][S.sub.1] > [r.sub.11] - [a.sup.(1).sub.12][M.sub.1], (27)

or equivalently (see (20) and (21))

[r.sub.11] - [B.sub.2] > [r.sub.11] - [B.sub.1], (28)

which is exactly (23). We have observed earlier that the net reproduction rate of the prey in patch 1 with dispersal corridors is larger than that without dispersal corridors. In the following, we continue to analyze how this larger net reproduction rate of the prey in patch 1 impacts the net reproduction rate of the predator in patch 1.

Consider the function g(x) = x exp(x - a) where a > 0 in the interval (0,+ra). The derivative of g(x) is g'(x) = (1 + x) exp(x - a) > 0, so g(x) is monotonically increasing in (0, +to). Thus, noting (27) we have

([r.sub.11] - [a.sup.(1).sub.12][S.sub.1])exp[[r.sub.11] - [a.sup.(1).sub.12][S.sub.1] - exp([r.sub.11] - 1)]

> ([r.sub.11] - [a.sup.(1).sub.12][M.sub.1])exp[[r.sub.11] - [a.sup.(1).sub.12][M.sub.1] - exp([r.sub.11] - 1)]. (29)

Noting the definitions of [B.sub.3] and [B.sub.4], inequality (29) shows that the dispersal rate [d.sub.12] which appeared in [S.sub.1] is advantageous to the increase of the net reproduction rate of the predator in patch 1, which in turn is helpful to the permanence of the predator in patch 1. In fact, from (27) (or (23)) this dispersion rate in [S.sub.1] increases the net reproduction rate of the prey in patch 1, which consequently increases the food supply for the predator in patch 1. As a result, this contributes positively to the net reproduction rate of the predator in patch 1, which further helps to sustain the permanence of the predator in patch 1.

In summary, the effects of the dispersion of the prey from patch 1 to patch 2, measured by the dispersal rate [d.sub.12], on the permanence of the predator in patch 1 are twofold: one is disadvantageous--contributed by the presence of [d.sub.12] in the denominator of [B.sub.4], while the other is advantageous--contributed by the presence of [d.sub.12] in [S.sub.1]. Similar conclusion can be drawn in patch 2 by comparing the condition (13) with (18). The dispersion of the prey from patch 2 to patch 1, measured by the dispersal rate [d.sub.21], has multiple effects, both advantageous and disadvantageous, on the permanence of the predator in patch 2.

Remark 7. From the two sets of permanence conditions in Theorems 4 and 5, no conclusion can be drawn on the effects of [d.sub.21] on the permanence of the prey and the predator in patch 1; also the effects of [d.sub.12] on the permanence of the prey and the predator in patch 2 are unknown.

5. Numerical Simulation

In this section, we present five numerical examples to confirm the theoretical results obtained earlier as well as to illustrate the effects of the dispersal corridors. All the examples are processed by MATLAB. The first example below confirms the correctness of Theorem 4.

Example 8. Consider the following predator-prey model (3) in patchy environment without dispersal corridors:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (30)

Direct computation gives

[r.sub.11] - [a.sup.(1).sub.12][M.sub.1] = 0.1324 > 0, -[r.sub.21] + [a.sup.(1).sub.21][m.sub.1] = 0.0090 > 0,

[r.sub.12] - [a.sup.(2).sub.12][M.sub.2] = 0.1524 > 0, -[r.sub.22] + [a.sup.(2).sub.21][m.sub.2] = 0.0023 > 0. (31)

Hence, the conditions of Theorem 4 are all satisfied. By Theorem 4, model 30) is permanent and this is illustrated in Figure 1.

The next example illustrates Theorem 5.

Example 9. Consider the following predator-prey model (4) in patchy environment with dispersal corridors:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (32)

In model (32), the dispersal rates are [d.sub.12] = 0.1 and [d.sub.21] = 0.15. Also, we set [alpha] = [a.sup.(1).sub.21] = 0.15 and [beta] = [a.sup.(2).sub.21] = 0.3; that is, we assume that the predation rate of the predator on the prey which migrated from another patch is the same as the predation rate on the prey originally in the patch. By direct computation, we have all the conditions of Theorem 5 satisfied as follows:

[r.sub.11] - [a.sup.(1).sub.12][S.sub.1] = 0.1121 > 0, -[r.sub.21] + [a.sup.(1).sub.21][s.sub.1] = 0.0063 > 0,

[r.sub.12] - [a.sup.(2).sub.12][S.sub.2] = 0.2246 > 0, -[r.sub.22] + [a.sup.(2).sub.21][s.sub.2] = 0.0109 > 0. (33)

Hence, by Theorem 5 model (32) is permanent. Figure 2 illustrates the permanence of model (32).

Next, we present an example to show that dispersal corridors may be disadvantageous to the permanence of the ecosystem.

Example 10. Consider the following two models:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (34)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (35)

The difference between (34) and (35) is only whether dispersal corridors exist or not; all other coefficients are the same. Model (34) is without dispersal corridors, and a direct computation shows that all the conditions of Theorem 4 are satisfied as follows:

[r.sub.11] - [a.sup.(1).sub.12] = 0.1230 > 0, -[r.sub.21] + [a.sup.(1).sub.21][m.sub.1] = 0.0038 > 0,

[r.sub.12] - [a.sup.(2).sub.12][M.sub.2] = 0.1542 > 0, -[r.sub.22] + [a.sup.(2).sub.21][m.sub.2] = 0.0023 > 0. (36)

Hence, by Theorem 4 model (34) is permanent. Figure 3(a) illustrates the permanence of model (34).

Model (35) has dispersal corridors with [d.sub.12] = 0.61, [d.sub.21] = 0, and [beta] = [a.sup.(2).sub.21] = 0.14. Checking the conditions of Theorem 5, we find

[r.sub.11] - [a.sup.(1).sub.12][S.sub.1] = 0.1635 > 0, -[r.sub.21] + [a.sup.(1).sub.21][s.sub.1] = -0.0382 < 0,

[r.sub.12] - [a.sup.(2).sub.12][S.sub.2] = 0.1542 > 0, -[r.sub.22] + [a.sup.(2).sub.21][s.sub.2] = 0.0023 > 0. (37)

Hence, not all the conditions of Theorem 5 are satisfied and we cannot ensure the permanence of model (35). In fact, by numerical simulation we find that the predator in patch 1 is driven to extinction and so model (35) is nonpermanent. This is depicted in Figure 3(b).

This example shows that the dispersion of the prey from patch 1 to patch 2 by the dispersal corridors may be disadvantageous to the predator in patch 1, which has been discussed in Section 4.2. Further, it is observed from the above calculations that the introduction of [d.sub.12] increases the net reproduction rate ([r.sub.11] - [a.sup.(1).sub.12][S.sub.1]) of the prey in patch 1 from 0.1230 to 0.1635, which has been noted in Section 4.1.

Next, we give an example to show that dispersal corridors may be helpful to the permanence of the ecosystem.

Example 11. Consider the following two models:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (38)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (39)

The difference between (38) and (39) is only whether dispersal corridors exist or not; all other coefficients are the same. Model (38) is without dispersal corridors, and a direct computation yields

[r.sub.11] - [a.sup.(1).sub.12][M.sub.1] = -0.0847 < 0, -[r.sub.21] + [a.sup.(1).sub.21][m.sub.1] = -0.0411 < 0,

[r.sub.12] - [a.sup.(2).sub.12][M.sub.2] = 0.1737 > 0, -[r.sub.22] + [a.sup.(2).sub.21][m.sub.2] = 0.0154 > 0. (40)

Hence, not all the conditions of Theorem 4 are satisfied and we cannot guarantee the permanence of (38). Indeed, by numerical simulation it is found that the predator in patch 1 will be driven to extinction. So model (38) is nonpermanent and this is illustrated in Figure 4(a).

Model (39) has dispersal corridors with [d.sub.12] > 0, [d.sub.21] > 0, [alpha] = [a.sup.(1).sub.21] = 0.2, and [beta] = [a.sup.(2).sub.21] = 0.3. Suppose

[d.sub.12] = 0.55, [d.sub.21] = 0.21. (41)

By numerical simulation we find that the predator in patch 1 is now permanent, and model (39), as a whole, is permanent. This is depicted in Figure 4(b). Note, however, that model (39) does not satisfy some of the conditions of Theorem 5, as

[r.sub.11] - [a.sup.(1).sub.12] = -0.0742 < 0, -[r.sub.21] + [a.sup.(1).sub.21][s.sub.1] = -0.0289 < 0,

[r.sub.12] - [a.sup.(2).sub.12][S.sub.2] = 0.2362 > 0, -[r.sub.22] + [a.sup.(2).sub.21][s.sub.2] = 0.0063 > 0. (42)

This simply reinforces the fact that the conditions in Theorem 5 are sufficient conditions.

We observe that the introduction of dispersal corridors leads to (i) an increase in the net reproduction rates of both the prey and the predator in patch 1 as well as the prey in patch 2 and (ii) a decrease in the net reproduction rate of the predator in patch 2. The overall effect on the permanence of model (39) is positive. This example shows that the introduction of dispersal corridors to a patchy ecosystem may be advantageous to the permanence of the ecosystem.

We will further investigate the effects of dispersal corridors on model (39). For this, in Table 2, for different values of [d.sub.12] and [d.sub.21] we list the net reproduction rate [R.sup.x.sub.i] of the prey x in patch i, the net reproduction rate [R.sup.y.sub.i] of the predator y in patch i, i [member of] {1,2}, and the permanence status of model (39). Recall that, for a model with dispersal corridors,

[R.sup.x.sub.1] = [r.sub.11] - [a.sup.(1).sub.12][S.sub.1], [R.sup.y.sub.1] = - [r.sub.21] + [a.sup.(1).sub.21][s.sub.1],

[R.sup.x.sub.2] = [r.sub.12] - [a.sup.(2).sub.12][S.sub.2], [R.sup.y.sub.2] = - [r.sub.22] + [a.sup.(2).sub.21][s.sub.2]. (43)

From Table 2, we observe the following.

(a) The effect of dispersal corridors on the net reproduction rate of prey is positive; specifically, [R.sup.x.sub.i] increases as [d.sub.12] increases and [R.sup.x.sub.2] increases as [d.sub.21] increases. This coincides with the discussion in Section 4.1.

(b) The effects of dispersal corridors on the net reproduction rate of predator are twofold: advantageous and disadvantageous. For a fixed [d.sub.12] = 0.55, [R.sup.y.sub.2] first increases and then decreases as the value of [d.sub.21] increases from 0 to 0.28.

(c) The entire ecosystem cannot be permanent by the sole effect of [d.sub.12] (i.e., [d.sub.21] = 0). If the value of [d.sub.12] increases to sufficiently large, even the prey in patch 1 will undergo extinction due to unpredicted environmental changes, as the population density of the prey in patch 1 will become very small. If the prey in patch 1 becomes extinct, so will the predator in patch 1, since there are no food sources.

(d) Noting (c) above, we further observe that the ecosystem is also nonpermanent if [d.sub.21] is too small; an appropriate value of [d.sub.21] is needed to induce a permanent (39).

So far equilibrium dynamics have been observed in Examples 8-11. Since it is well known that Ricker type model has complex dynamics, we give an example with more complex dynamics to illustrate our results.

Example 12. Consider the following two models:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (44)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (45)

The difference between (44) and (45) is only whether dispersal corridors exist or not; all other coefficients are the same. Model (44) is without dispersal corridors, and a direct computation yields

[r.sub.11] - [a.sup.(1).sub.12][M.sub.1] = 0.0772 > 0, -[r.sub.21] + [a.sup.(1).sub.21][m.sub.1] = -0.0989 < 0,

[r.sub.12] - [a.sup.(2).sub.12][M.sub.2] = 0.3283 > 0, -[r.sub.22] + [a.sup.(2).sub.21][m.sub.2] = 01056 > 0. (46)

Though not all the conditions of Theorem 4 are satisfied, model (44) is permanent which is illustrated in Figure 5(a). It is observed from Figure 5(a) that the subsystem of patch 2 has a stable positive equilibrium while the subsystem of patch 1 does not have such stable equilibrium though it is permanent. Further illustration of this subsystem of patch 1 is given in Figure 5(b). Figure 5(b) is the phase diagram of ([x.sub.1], [y.sub.1]) and it is observed that the subsystem of patch 1 of model (44) has a periodic solution of period 3. In fact, the periodic solution is ([x.sub.1], [y.sub.1]) = (18.7232, 1.1061), (4.3689, 11.9115), (6.1079, 0.5824).

Model (45) has dispersal corridors with [d.sub.21] = 0.6, [d.sub.12] = 0, and [alpha] = [a.sup.(1).sub.21] = 0.15. The net reproduction rates are computed as follows:

[r.sub.11] - [a.sup.(1).sub.12][S.sub.1] = 0.0772 > 0, -[r.sub.21] + [a.sup.(2).sub.21][s.sub.1] = -0.0989 < 0,

[r.sub.12] - [a.sup.(2).sub.12][S.sub.2] = 0.3453 > 0, -[r.sub.22] + [a.sup.(2).sub.21][s.sub.2] = -0.0213 < 0. (47)

The dynamics of (45) are presented in Figure 6(a). It is observed that the predator in patch 2 will become extinct and consequently model (45) is not permanent. Figure 6(b) illustrates that the periodic solution observed in the subsystem of patch 1 of model (44) now vanishes due to the dispersion.

Comparing (46) and (47), we observe that when [d.sub.21] > 0, (i) the net reproduction rate of the prey in patch 2 increases, which has been noted in Section 4.1 and (ii) the net reproduction rate of the predator in patch 2 decreases its value from positive to negative; this shows that the dispersion of the prey from patch 2 to patch 1 has an adverse effect on the permanence of the predator in patch 2, which has been discussed in Section 4.2. This example shows that large dispersal rates of the prey may destroy the permanence of a predator-prey system.

6. Conclusion

The paper is mainly involved with the effects of dispersal corridors on the permanence of discrete predator-prey models in patchy environment. We have introduced two models of Ricker type undergoing the effect of density dependency--model (3) is in patchy environment without dispersal corridors and model (4) is with dispersal corridors. Two sets of sufficient conditions to guarantee the permanence of the two models are then obtained through theoretical studies, respectively. In Theorem 4, conditions (9)-(13) guarantee the permanence of model (3); this result is illustrated by Example 8. In Theorem 5, conditions (15)-(18) ensure the permanence of model (4); this result is also numerically illustrated by Example 9.

Through mathematical analysis of the two sets of sufficient conditions, combined with the actual biological background of models (3) and (4), we have discussed the effects of dispersal corridors on the permanence of discrete predator-prey model. It is found that the existence of dispersal corridors for the prey to migrate, such as from patch i to patch j, while the predator remains fixed to the original patch, is helpful to the permanence of the prey in patch i under density dependent circumstances due to an increased net reproduction rate of the prey in patch i. Moreover, it is also found that the migration of prey, such as from patch i to patch j, is both advantageous and disadvantageous to the predator in patch i because the dispersal rate [d.sub.ij] has both positive and negative effects on the net reproduction rate of the predator in patch i. Examples 10-12 are presented to illustrate the positive/negative effects of the dispersal corridors.

Appendix

In this section, we present the proofs of Theorems 2, 4, and 5. First we state a lemma which is useful in proving the theorems; this lemma is the autonomous case of Lemma 2 in [43].

Lemma 13. Suppose {z(k)} satisfies

z(k+ 1) [greater than or equal to] z (k) exp [r - az (k)], (A.1)

for k [greater than or equal to] K, where a and r are positive constants, K is a positive integer, and z(K) > 0. Further, assume that lim [sup.sub.k [right arrow] +[infinity]]z(k) [less than or equal to] A and aA/r > 1, where A is a positive constant. Then,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.2)

Proof of Theorem 2. We will first establish the first inequality in Definition 1. To begin, from the first equation of (7) we have

[x.sub.1](n + 1) [less than or equal to] [x.sub.1](n) exp ([r.sub.11] - [a.sup.(1).sub.11][x.sub.1](n)). (A.3)

It is straightforward that the maximum of the function z(x) = x exp(r - ax) where a > 0 in the interval [0, +[infinity]) is [z.sub.max] = (1/a) exp(r - 1). Therefore, it follows from (A.3) that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.4)

Let [epsilon] > 0 be arbitrarily given. From (A.4), there exists a positive integer [N.sub.1] such that

[x.sub.1](n) [less than or equal to] [A.sub.1] + [epsilon], n > [N.sub.1]. (A.5)

Using (A.5) in the second equation of (7) gives

[y.sub.1](n + 1) [less than or equal to] [y.sub.1](n)exp[-[r.sub.21] + [a.sup.(1).sub.21]([A.sub.1] + [epsilon]) - [a.sup.(1).sub.22][y.sub.1](n)],

n > [N.sub.1]. (A.6)

Hence, using the maximum of z(x) = x exp(r - ax) as earlier, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.7)

Since [epsilon] is arbitrary, letting [epsilon] [right arrow] 0 in (A.7) yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.8)

Now, in view of inequality (9), we can choose [[epsilon].sub.1] > 0 such that

[r.sub.11] - [a.sup.(1).sub.12]([M.sub.1] + [[epsilon].sub.1]) > 0. (A.9)

Further, noting (A.8) we have from the first equation of (7)

[x.sub.1](n + 1) [greater than or equal to] [x.sub.1](n) exp [[r.sub.11] - [a.sup.(1).sub.11][x.sub.1](n) - [a.sup.(1).sub.12]([M.sub.1] + [[epsilon].sub.1])], (A.10)

for sufficiently large n. Using the inequality exp(x - 1) [greater than or equal to] x which holds for x > 0, together with the definition of [A.sub.1] (see (A.4)), we find

[a.sup.(1).sub.11][A.sub.1] = exp([r.sub.11] - 1) [greater than or equal to] [r.sub.11] > [r.sub.11] - [a.sup.(1).sub.12]([M.sub.1] + [[epsilon].sub.1]) > 0, (A.11)

or

[[a.sup.(1).sub.11][A.sub.1]/[[r.sub.11] - [a.sup.(1).sub.12]([M.sub.1] + [[epsilon].sub.1])]] > 1. (A.12)

Inequalities (A.4) and (A.12) imply that Lemma 13 can be applied to (A.10). Thus, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.13)

Again let [[epsilon].sub.1] [right arrow] 0 in (A.13) and it follows that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.14)

It is clear from (9) that [m.sub.1] > 0.

We will further show that [m.sub.1] < [A.sub.1]. From the derivation of (A.12), we have

[a.sup.(1).sub.11][A.sub.1] > [r.sub.11] - [a.sup.(1).sub.12][M.sub.1]; (A.15)

that is, [r.sub.11] - [a.sup.(1).sub.12][M.sub.1] - [a.sup.(1).sub.11][A.sub.1] < 0, and so

exp[[r.sub.11] - [a.sup.(1).sub.12][M.sub.1] - [a.sup.(1).sub.11][A.sub.1]] < 1. (A.16)

Therefore, from the definition of [m.sub.1] (see (A.14)) we find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.17)

In summary, combining (A.4), (A.14), and (A.17) we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (a 18)

In the following, we will establish the second inequality in Definition 1. In fact, we will prove that there exists [m.sub.3] > 0 such that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.19)

Recalling that [m.sub.1] > 0 and noting inequality (10), we can choose [[epsilon].sub.2] > 0 sufficiently small such that [m.sub.1] - [[epsilon].sub.2] > 0 and -[r.sub.21] + [a.sup.(1).sub.21]([m.sub.1] - [[epsilon].sub.2]) > 0. From the second equation of (7) and (A.14), we have

[y.sub.1](n + 1) [greater than or equal to] [y.sub.1](n)exp[-[r.sub.21] + [a.sup.(1).sub.21]([m.sub.1] - [[epsilon].sub.2]) - [a.sup.(1).sub.22][y.sub.1](n), (A.20)

for sufficiently large n. Using the inequality exp(x - 1) [greater than or equal to] x which holds, for x > 0, and the fact that [m.sub.1] < [A.sub.1], we find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.21)

Hence,

[[a.sup.(1).sub.22][M.sub.1]/[-[r.sub.21] + [a.sup.(1).sub.21]([m.sub.1] - [[epsilon].sub.2]]] > 1, (A.22)

which, together with (A.8), indicates that we can apply Lemma 13 to (A.20) to give

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.23)

Letting [[epsilon].sub.2] [right arrow] 0, we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (A.24)

It is obvious from (10) that [m.sub.3] > 0.

We will further show that [m.sub.3] < [M.sub.1]. From (A.21), we have

-[r.sub.21] + [a.sup.(1).sub.21][m.sub.1] - [a.sup.(1).sub.22][M.sub.1] < 0, (A.25)

and so

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.26)

Combining (A.8), (A.24), and (A.26), we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.27)

In view of (A.18) and (A.27), we have obtained both inequalities in Definition 1 and so the permanence of model (7) with initial values (8) is established. This completes the proof of Theorem 2. []

Proof of Theorem 4. Using a similar argument as in the proof of Theorem 2, we see that, under the assumptions (12) and (13), the following model

[x.sub.2](n + 1) = [x.sub.2](n)exp[[r.sub.12] - [a.sup.(2).sub.11][x.sub.2](n) - [a.sup.(2).sub.12][y.sub.2](n)],

[y.sub.2](n + 1) = [y.sub.2](n)exp[-[r.sub.22] + [a.sup.(2).sub.21][x.sub.2](n) - [a.sup.(2).sub.22][y.sub.2](n)] (A.28)

with initial values [x.sub.2](0) = [x.sub.2] > 0, [y.sub.2](0) = [y.sub.2] > 0 is permanent. In fact, analogous to (A.18) and (A.27), we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (A.29)

where [m.sub.2], [A.sub.2], and [M.sub.2] are given in Theorem 4 and

[m.sub.4] = [[-[r.sub.22] + [a.sup.(2).sub.21][m.sub.2]]/[a.sup.(2).sub.22]]exp(-[r.sub.22] + [a.sup.(2).sub.21][m.sub.2] - [a.sup.(2).sub.22][M.sub.2]). (A.30)

Coupling with Theorem 2, the permanence of model (3) with initial values (5) follows immediately. This completes the proof of Theorem 4. []

Proof of Theorem 5. Consider the following auxiliary system:

u(n + 1) = w(n) exp [[r.sub.11] - [a.sup.(1).sub.11] u(n) - [a.sup.(1).sub.12]v(n) - [d.sub.12]u(n)],

v(n + 1) = v(n) exp [-[r.sub.21] + [a.sup.(1).sub.21]u n) - [a.sup.(1).sub.22]v(n)] (A.31)

with initial values u(0) = [x.sub.1] > 0, v(0) = [y.sub.1] >0.

Comparing (7) and (A.31), it follows from Theorem 2 that model (A.31) is permanent under conditions (15) and (16). Therefore, there exist positive constants [m.sup.(1).sub.1] and [m.sup.(2).sub.1] such that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.32)

From the first two equations of (4), it is clear that, for any N, where N represents the set of all positive integers, we have [x.sub.1](n) [greater than or equal to] u(n) and [y.sub.1](n)[greater than or equal to]> v(n). Hence,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.33)

Next, we consider the following auxiliary system:

P(n + 1) = p(n)exp[[r.sub.12] - [a.sup.(2).sub.11]p(n) - [a.sup.(2).sub.12]q(n) - [d.sub.21]p(n)],

q(n + 1) = q(n)exp[-[r.sub.22] + a(2).sub.21p(n) - [a.sup.(2).sub.22]q(n)] (A.34)

with initial values p(0) = [x.sub.2] > 0, q(0) = [y.sub.2] > 0. By a similar argument as above, we see that (A.34) is permanent under conditions (17) and (18). Thus, there exist positive constants [m.sup.(1).sub.2] and [m.sup.(2).sub.2] such that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.35)

In view of the last two equations of (4), we have [x.sub.2](n) [greater than or equal to] p(n) and [y.sub.2](n) [greater than or equal to] q(n), for any n [member of] N. It follows that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.36)

To complete the proof, in the following we will prove that there exist positive constants [M.sup.(j).sub.i], i, j [member of] {1,2}, such that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.37)

Suppose the contrary; there exists no such [M.sup.(1).sub.1] > 0 to satisfy lim [sup.sub.n [right arrow] +[infinity]][x.sub.1](n) [less than or equal to] [M.sup.(1).sub.1], and then we can select a subsequence [{[x.sub.1]([n.sub.k])}.sup.[infinity].sub.k=1] such that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.38)

Multiplying the first equation with the third equation of (4), we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.39)

Note that the two-variable function f(x, y) = xy exp(a - bx - cy) is a bounded function in the first quadrant {(x, y) | x > 0, y > 0} provided that a, b, c > 0. Hence, it follows from (A.39) that lim [sup.sub.n [right arrow] +[infinity]][x.sub.1](n)[x.sub.2](n) is finite. Noting (A.38), this implies [lim.sub.k [right arrow] +[infinity]][x.sub.2]([n.sub.k]) = 0, which is a contradiction to (A.36). This shows that lim [sup.sub.n [right arrow] +[infinity]][x.sub.1](n) [less than or equal to] [M.sup.(1).sub.1]. Similarly, we can prove that lim [sup.sub.n [right arrow] +[infinity]] [x.sub.2](n) [less than or equal to] [M.sup.(1).sub.2]. Summarizing, we have established the following inequalities:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.40)

Now, from the second equation of (4) and (A.40), we have, for sufficiently small [[epsilon].sub.3] > 0,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (A.41)

for sufficiently large n. Noting that the maximum of the function z(x) = x exp(r - ax) where a > 0 in the interval [0, +[infinity]) is [z.sub.max] = (1/a) exp(r - 1), it follows from (A.41) that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.42)

Letting [[epsilon].sub.3] [right arrow] 0 in the above inequality, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (A.43)

Likewise, a similar deduction as above gives

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.44)

Combining (A.40), (A.43), and (A.44), we have shown that (A.37) holds. The permanence of model (4) with initial values (5) now follows from (A.33), (A.36), and (A.37). This completes the proof of Theorem 5. []

http://dx.doi.org/10.1155/2014/140902

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank the reviewers for their comments which help to improve the paper. The research is sponsored by Jiangsu Overseas Research & Training Program.

References

[1] P. Schuster, K. Sigmund, and R. Wolff, "Dynamical systems under constant organization. III. Cooperative and competitive behavior of hypercycles," Journal of Differential Equations, vol. 32, no. 3, pp. 357-368, 1979.

[2] C. Wu and J. Cui, "Global dynamics of discrete competitive models with large intrinsic growth rates," Discrete Dynamics in Nature and Society, vol. 2009, Article ID 710353, 15 pages, 2009.

[3] C. Wu, "Permanence and stable periodic solution for a discrete competitive system with multidelays," Advances in Difference Equations, vol. 2009, Article ID 375486, 12 pages, 2009.

[4] Z. Zhou and X. Zou, "Stable periodic solutions in a discrete periodic logistic equation," Applied Mathematics Letters, vol. 16, no. 2, pp. 165-171, 2003.

[5] C. Celik and O. Duman, "Allee effect in a discrete-time predator-prey system," Chaos, Solitons and Fractals, vol. 40, no. 4, pp. 1956-1962, 2009.

[6] J. Bastinec, L. Berezansky, J. Diblik, and Z. Smarda, "On a delay population model with a quadratic nonlinearity without positive steady state," Applied Mathematics and Computation, vol. 227, pp. 622-629, 2014.

[7] L. Berezansky, J. Bastinec, J. Diblik, and Z. Smarda, "On a delay population model with quadratic nonlinearity," Advances in Difference Equations, vol. 2012, article 230, 2012.

[8] V. A. A. Jansen and K. Sigmund, "Shaken not stirred: on permanence in ecological communities," Theoretical Population Biology, vol. 54, no. 3, pp. 195-201, 1998.

[9] A. Hastings, "Complex interactions between dispersal and dynamics: lessons from coupled logistic equations," Ecology, vol. 74, no. 5, pp. 1362-1372, 1993.

[10] M. D. Holland and A. Hastings, "Strong effect of dispersal network structure on ecological dynamics," Nature, vol. 456, no. 7223, pp. 792-795, 2008.

[11] Y. Kim, O. Kwon, and F. Li, "Evolution of dispersal toward fitness," Bulletin of Mathematical Biology, vol. 75, no. 12, pp. 2474-2498, 2013.

[12] S. A. Levin, "Dispersal and population interactions," American Naturalist, vol. 108, pp. 207-228, 1974.

[13] Y. Lou and F. Lutscher, "Evolution of dispersal in open advective environments," Journal of Mathematical Biology, 2013.

[14] N. Moquet, M. F. Hoopes, and P. Amarasekare, "The world is patchy and heterogeneous! Trade-off and source-sink dynamics in competitive metacommunities," in Metacommunities: Spatial Dynamics and Ecological Communities, chapter 10, pp. 237-262, University of Chicago Press, Chicago, 111, USA, 2005.

[15] S. J. Schreiber and T P. Killingback, "Spatial heterogeneity promotes coexistence of rock-paper-scissors metacommunities," Theoretical Population Biology, vol. 86, pp. 1-11, 2013.

[16] X. Zhang and W. D. Wang, "Importance of dispersal adaptations of two competitive populations between patches," Ecological Modelling, vol. 222, no. 1, pp. 11-20, 2011.

[17] http://www.china.com.cn/city/txt/2007-07/02/content8466486.htm.

[18] F. Kong, Studies on the adaptability and behavior time assignments of Tibtan antelope using wild life passageways of Qinghai-Tibet railway, [Master dissertation], Northwest University, Shaanxi, China, 2009, (Chinese).

[19] J. Cui and L. Chen, "Permanence and extinction in logistic and Lotka-Volterra systems with diffusion," Journal of Mathematical Analysis and Applications, vol. 258, no. 2, pp. 512-535, 2001.

[20] J. Cui, Y. Takeuchi, and Z. Lin, "Permanence and extinction for dispersal population systems," Journal of Mathematical Analysis and Applications, vol. 298, no. 1, pp. 73-93, 2004.

[21] Z. Lu and Y. Zhou, Advances in Mathematical Biology, Science Press, Beijing, China, 2006, (Chinese).

[22] C. Wu and J. Cui, "Permanence for a delayed discrete predator-prey model with prey dispersal," International Journal of Biomathematics, vol. 2, no. 3, pp. 311-320, 2009.

[23] L. Zhang and Z. Teng, "Permanence for a delayed periodic predator-prey model with prey dispersal in multi-patches and predator density-independent," Journal of Mathematical Analysis and Applications, vol. 338, no. 1, pp. 175-193, 2008.

[24] E. Braverman and D. Kinzebulatov, "On linear perturbations of the Ricker model," Mathematical Biosciences, vol. 202, no. 2, pp. 323-339, 2006.

[25] Y. Kang, D. Armbruster, and Y. Kuang, "Dynamics of a plant-herbivore model," Journal of Biological Dynamics, vol. 2, no. 2, pp. 89-101, 2008.

[26] G. Sun, G. Zhang, and Z. Jin, "Dynamic behavior of a discrete modified Ricker and Beverton-Holt model," Computers and Mathematics with Applications, vol. 57, no. 8, pp. 1400-1412, 2009.

[27] Y. Chow and S. R. Jang, "Coexistence in a discrete competition model with dispersal," Journal of Difference Equations and Applications, vol. 19, no. 4, pp. 615-632, 2013.

[28] J. A. L. Silva and F. T. Giordani, "Density-dependent migration and synchronism in metapopulations," Bulletin of Mathematical Biology, vol. 68, no. 2, pp. 451-465, 2006.

[29] J. A. L. Silva and F. T. Giordani, "Density-dependent dispersal in multiple species metapopulations," Mathematical Biosciences and Engineering, vol. 5, no. 4, pp. 843-857, 2008.

[30] A. Yakubu, "Asynchronous and synchronous dispersals in spatially discrete population models," SIAM Journal on Applied Dynamical Systems, vol. 7, no. 2, pp. 284-310, 2008.

[31] R. M. May, "Nonlinear problems in ecology," in Chaotic Behavior of Deterministic Systems, G. Iooss, R. Helleman, and R. Stora, Eds., North-Holland, Amsterdam, The Netherlands, 1983.

[32] E. W. Ricker, "Stock and recruitment," Journal of the Fisheries Research Board of Canada, vol. 11, no. 5, pp. 559-623, 1954.

[33] A. Hastings, Population Biology: Concepts and Models, Springer, New York, NY, USA, 1996.

[34] H. L. Smith, Monotone Dynamical Systems: An Introduction to the Theory of Competitive and Cooperative Systems, vol. 41 of Mathematical Surveys and Monographs, American Mathematical Society, Providence, RI, USA, 1995.

[35] Z. Lu and W. Wang, "Permanence and global attractivity for Lotka-Volterra difference systems," Journal of Mathematical Biology, vol. 39, no. 3, pp. 269-282, 1999.

[36] W. Wang and Z. Lu, "Global stability of discrete models of Lotka-Voltera type," Nonlinear Analysis: Theory, Methods & Applications, vol. 35, no. 8, pp. 1019-1030, 1999.

[37] M. Y. Li and Z. Shuai, "Global-stability problem for coupled systems of differential equations on networks," Journal of Differential Equations, vol. 248, no. 1, pp. 1-20, 2010.

[38] L. Zu, D. Q. Jiang, and F. Q. Jiang, "Existence, stationary distribution, and extinction of predator-prey system of prey dispersal with stochastic perturbation," Abstract and Applied Analysis, vol. 2012, Article ID 547152, 24 pages, 2012.

[39] F. S. Berezovskaya, B. J. Song, and C. Castillo-Chavez, "Role of prey dispersal and refuges on predator-prey dynamics," SIAM Journal on Applied Mathematics, vol. 70, no. 6, pp. 1821-1839, 2010.

[40] G. E. Forrester, "Influences of predatory fish on the drift dispersal and local density of stream insects," Ecology, vol. 75, no. 5, pp. 1208-1218, 1994.

[41] L. J. S. Allen, "Persistence and extinction in single-species reaction-diffusion models," Bulletin of Mathematical Biology, vol. 45, no. 2, pp. 209-227, 1983.

[42] C. Wu, "On the global asymptotic behavior of a discrete predator-prey model," Advances and Applications in Mathematical Sciences, vol. 10, no. 2, pp. 143-157, 2011.

[43] X. Yang, "Uniform persistence and periodic solutions for a discrete predator-prey system with delays," Journal of Mathematical Analysis and Applications, vol. 316, no. 1, pp. 161-177, 2006.

Chunqing Wu, (1) Shengming Fan, (1) and Patricia J. Y. Wong (2)

(1) School of Mathematics and Physics, Changzhou University, Changzhou 213164, China

(2) School of Electrical and Electronic Engineering, Nanyang Technological University, Singapore 639798

Correspondence should be addressed to Patricia J. Y. Wong; ejywong@ntu.edu.sg

Received 25 March 2014; Revised 11 June 2014; Accepted 12 June 2014; Published 17 July 2014

Academic Editor: Josef Diblik

TABLE 1: Sufficient conditions for the permanence of models with-out/with dispersal corridors. Concerned with the Without dispersal With dispersal permanence corridors corridors The prey a in patch 1 (9) (15) The predator y in patch 1 (10) (16) The prey a in patch 2 (12) (17) The predator y in patch 2 (13) (18) TABLE 2: Net reproduction rates and the permanence of model (39) under different dispersal rates. Dispersal rates [R.sup.x.sub.1] [R.sup.y.sub.1] [d.sub.12] = 0, -0.0847 -0.0411 [d.sub.21] = 0 [d.sub.12] = 0.18, -0.0793 -0.0346 [d.sub.21] = 0 [d.sub.12] = 0.36, -0.0762 -0.0312 [d.sub.21] = 0 [d.sub.12] = 0.55, -0.0742 -0.0289 [d.sub.21] = 0 [d.sub.12] = 1, -0.0715 -0.0260 [d.sub.21] = 0 [d.sub.12] = 0.55, -0.0742 -0.0289 [d.sub.21] = 0.02 [d.sub.12] = 0.55, -0.0742 -0.0289 [d.sub.21] = 0.07 [d.sub.12] = 0.55, -0.0742 -0.0289 [d.sub.21] = 0.14 [d.sub.12] = 0.55, -0.0742 -0.0289 [d.sub.21] = 0.21 [d.sub.12] = 0.55, -0.0742 -0.0289 [d.sub.21] = 0.28 Dispersal rates [R.sup.x.sub.2] [R.sup.y.sub.2] [d.sub.12] = 0, 0.1737 0.0154 [d.sub.21] = 0 [d.sub.12] = 0.18, 0.1737 0.0154 [d.sub.21] = 0 [d.sub.12] = 0.36, 0.1737 0.0154 [d.sub.21] = 0 [d.sub.12] = 0.55, 0.1737 0.0154 [d.sub.21] = 0 [d.sub.12] = 1, 0.1737 0.0154 [d.sub.21] = 0 [d.sub.12] = 0.55, 0.1835 0.0163 [d.sub.21] = 0.02 [d.sub.12] = 0.55, 0.2030 0.0157 [d.sub.21] = 0.07 [d.sub.12] = 0.55, 0.2224 0.0116 [d.sub.21] = 0.14 [d.sub.12] = 0.55, 0.2362 0.0063 [d.sub.21] = 0.21 [d.sub.12] = 0.55, 0.2464 0.00086 [d.sub.21] = 0.28 Dispersal rates Permanence of model (39) [d.sub.12] = 0, No, predator in [d.sub.21] = 0 patch 1 is extinct [d.sub.12] = 0.18, No, predator in [d.sub.21] = 0 patch 1 is extinct [d.sub.12] = 0.36, No, predator in [d.sub.21] = 0 patch 1 is extinct [d.sub.12] = 0.55, No, predator in [d.sub.21] = 0 patch 1 is extinct [d.sub.12] = 1, No, predator in [d.sub.21] = 0 patch 1 is extinct [d.sub.12] = 0.55, No, predator in [d.sub.21] = 0.02 patch 1 is extinct [d.sub.12] = 0.55, Yes [d.sub.21] = 0.07 [d.sub.12] = 0.55, Yes [d.sub.21] = 0.14 [d.sub.12] = 0.55, Yes [d.sub.21] = 0.21 [d.sub.12] = 0.55, Yes [d.sub.21] = 0.28

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Wu, Chunqing; Fan, Shengming; Wong, Patricia J.Y. |

Publication: | Abstract and Applied Analysis |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 10445 |

Previous Article: | Adaptive backstepping control based on floating offshore high temperature superconductor generator for wind turbines. |

Next Article: | Study of a biparametric family of iterative methods. |

Topics: |