Abstract and Applied Analysis

Volume 2014 (2014), Article ID 689573, 15 pages

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

## Global Stability Analysis of a Nonautonomous Stage-Structured Competitive System with Toxic Effect and Double Maturation Delays

^{1}Institute of Systems Science, Northeastern University, Shenyang 110004, China^{2}State Key Laboratory of Integrated Automation of Process Industry, Northeastern University, Shenyang 110004, China^{3}Institute of Information and Computational Science, Northeastern University at Qinhuangdao, Qinhuangdao 066004, China

Received 25 February 2014; Revised 19 August 2014; Accepted 22 August 2014; Published 19 October 2014

Academic Editor: Zhiming Guo

Copyright © 2014 Chao Liu and Yuanke Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We investigate a nonautonomous two-species competitive system with stage structure and double time delays due to maturation for two species, where toxic effect of toxin liberating species on nontoxic species is considered and the inhibiting effect is zero in absence of either species. Positivity and boundedness of solutions are analytically studied. By utilizing some comparison arguments, an iterative technique is proposed to discuss permanence of the species within competitive system. Furthermore, existence of positive periodic solutions is investigated based on continuation theorem of coincidence degree theory. By constructing an appropriate Lyapunov functional, sufficient conditions for global stability of the unique positive periodic solution are analyzed. Numerical simulations are carried out to show consistency with theoretical analysis.

#### 1. Introduction

In recent years, many research efforts have been made on competitive Lotka-Volterra system with stage structure and time delay. By incorporating a constant time delay into single species model, a stage-structured model is proposed in the pioneering work [1], where time delay reflects a delayed birth of immature population and a reduced survival of immature population to their maturity. The model system takes the following form: where and represent the immature population and mature population density at time , respectively. denotes the birth rate of immature population; stands for the death rate of immature population. is the death and overcrowding rate of mature population. denotes time of immature population to maturity. The term represents the immature species which are born at time and survive at time with immature death rate and therefore represents transformation of immature species to mature species. It is found that all ecologically relevant solutions tend to the positive equilibrium solution as time , and various aspects of the above proposed system including positivity and boundedness of solutions are discussed in [1].

Zeng et al. propose a nonautonomous competitive two-species model with stage structure in one species in [2], where conditions of permanence are obtained. Furthermore, existence and asymptotic stability of periodic solution are proved under some assumptions if the proposed model turns out to be a periodic system. A two-species Lotka-Volterra type competition model with stage structure for both species is proposed and investigated in [3], where the individuals of each species are classified as immature and mature. By constructing a suitable Lyapunov function, sufficient conditions are derived for the global stability of nonnegative equilibria of the proposed model in the case of constant coefficients. Furthermore, a set of easily verifiable sufficient conditions are obtained for the existence of positive periodic solution when coefficients are assumed to be positively continuous periodic functions. In [4], there is a time delayed periodic system which describes the competition among mature populations. The evolutionary behavior of model system is analyzed and some sufficient conditions for competitive coexistence and exclusion are obtained.

A nonautonomous competitive Lotka-Volterra system is studied in [5]; it reveals a computable necessary and sufficient condition for the system to be totally permanent when the growth rates have averages and the interaction coefficients are nonnegative constants. Along with this research, permanence for a class of competitive Lotka-Volterra systems is discussed in [6] which extends the work done in [5], and a computable necessary and sufficient condition is found for the permanence of all subsystems of the system and its small perturbation on the interaction matrix. In [7], a two-species competitive model with stage structure is discussed, and the dynamics of coupled system of semilinear parabolic equations with time delays are investigated, which show that the introduction of diffusion does not affect the permanence and extinction of the species, though the introduction of stage structure brings negative effect on it. In [8], sufficient conditions are obtained for the existence of a unique, globally attractive, strictly positive (componentwise), almost periodic solution of a nonautonomous, almost periodic competitive two-species model with a stage structure in one species. An example together with its numeric simulations shows the feasibility of our main results, which generalize the main results of Zeng et al. [2]. According to two types of well-known periodic single species population growth models with time delay, two corresponding periodic competitive systems with multiple delays are proposed in [9], and the same criteria for the existence and globally asymptotic stability of positive periodic solutions of the above two competitive systems are derived. In [10], a discrete periodic competitive model with stage structure is established, and some sufficient and realistic conditions are obtained for existence of a positive periodic solution of the proposed system. In [11], a periodic nonautonomous competitive stage-structured system with infinite delay is considered, where the adult members of -species are in competition. For each of the -species the model incorporates a time delay which represents the time from birth to maturity of that species. Infinite delay is introduced which denotes the influential effect of the entire past history of the system on the current competition interactions. By using the comparison principle, if the growth rates are sufficiently large, then the solutions are uniformly permanent. Then, by using Horn’s fixed point theorem, the existence of positive periodic solution of the system with finite delay is discussed. Finally, it is proved that even the system with infinite delay admits a positive periodic solution.

In [12], a nonautonomous predator-prey system with discrete time delay is studied, where there is epidemic disease in the predator. By using some techniques of the differential inequalities and delay differential inequalities, the permanence of system is discussed under some appropriate conditions. When all the coefficients of the system are periodic, the existence and global attractivity of the positive periodic solution are studied by Mawhin’s continuation theorem and constructing a suitable Lyapunov functional. Furthermore, when the coefficients of the system are not absolutely periodic but almost periodic, sufficient conditions are also derived for the existence and asymptotic stability of the almost periodic solution. In [13], general -species nonautonomous Lotka-Volterra competitive systems with pure-delays and feedback controls are discussed. New sufficient conditions, for which a part of the -species remains permanent, are established by applying the method of multiple Lyapunov functionals and introducing a new analysis technique.

By utilizing Brouwer fixed point theorem and constructing a suitable Lyapunov function, the periodic solution and global stability for a nonautonomous competitive Lotka-Volterra diffusion system are investigated in [14]; it can be found that the system has a unique periodic solution which is globally stable under some appropriate conditions. In [15], a delay differential equation model for the interaction between two species is investigated. The maturation delay for each species is modelled as a distribution, to allow for the possibility that individuals may take different amount of time to maturity. Positivity and boundedness of the solutions are studied, and global stability is analyzed for each equilibrium. A Lotka-Volterra competitive system with infinite delay and feedback controls is proposed in [16]. By using Mawhin’s continuation theorem of coincidence degree theory, an impulsive nonautonomous Lotka-Volterra predator-prey system with harvesting terms is investigated in [17]. Some sufficient conditions for the existence of multiple positive almost periodic solutions for the system under consideration are discussed. Furthermore, existence of multiple positive almost periodic solutions to other types of population systems can be studied by using the same method obtained in this paper. By using the method of multiple Lyapunov functionals and by developing a new analysis technique, some sufficient conditions are obtained that guarantee that some of the species are driven to extinction. A three-dimensional nonautonomous competitive Lotka-Volterra system is considered in [18]; it is shown that if the growth rates are positive, bounded, and continuous functions, and the averages of the growth rates satisfy certain inequalities, then any positive solution has the property that one of its components vanishes. In [19], an almost periodic multispecies Lotka-Volterra mutualism system with time delays and impulsive effects is investigated. By using the theory of comparison theorem and constructing a suitable Lyapunov functional, sufficient conditions which guarantee the existence and uniqueness and global asymptotical stability of almost periodic solution of this system are obtained.

It is well known that the effect of toxins on ecological systems is an important issue from mathematical and experimental points of view [20, 21]. The first mathematical model to represent the toxic liberating interaction between two competing species is introduced by Maynard Smith [22]. The model is based upon a two-species Lotka-Volterra competition model with an additional term to take into account the effect of toxic substances released by one species to another, which takes the following form: where and represent the density of two competing species at time , respectively. and denote the birth rate of species and species, respectively. and are the rate of intraspecific competition term for the first and second species, respectively. and stand for the rate of interspecific competition, respectively. and represent the toxic inhibition rate for the first species by the second species and vice versa. By considering that denotes the rate of toxic inhibition for the nontoxic species released by the toxin liberating species and all other parameters share the same biological interpretations mentioned in model system (2), work done in [22] is extended in [23] and the generalized model system is as follows: where the toxic substance producing action follows the mathematical term [23].

It should be noted that models of the persistence and extinction of a population or community in a polluted environment have been investigated in [23]. But all of those papers have a basic assumption that the capacity of the environment is so large that the change of toxicant in the environment that comes from uptake and egestation by the organisms can be neglected. This assumption is not made in [24, 25], some sufficient conditions on persistence or extinction of a population have been obtained, and the threshold between the two has also been obtained for most situations. In [26, 27], there are modified delay differential equation models of the growth of two species of plankton having competitive and allelopathic effects on each other. By using the continuation theorem of coincidence degree theory, a set of easily verifiable sufficient conditions are obtained for the existence of positive periodic solutions for this model. Recently, some discussions and investigations of the nonautonomous competitive model with toxic effects are made. A periodic competitive stage-structured Lotka-Volterra model with the effects of toxic substances is investigated in [28]. It is shown that toxic substances play an important role in the extinction of species. A set of sufficient conditions guarantee that one of the components is driven to extinction while the other is globally attractive. The dynamical behavior of a two-species competitive system affected by toxic substances is investigated in [21], where each species produces a substance toxic to the other species. Boundedness and local and global stabilities are also addressed. It should be noted that toxic interaction follows the mathematical term suggested in model system (2) and each mature individual produces a substance toxic to the other mature individuals only when the other mature individual is present, and the immature individual is not affected by the toxicant [21]. However, to the author’s best knowledge, dynamical behavior and stability analysis of nonautonomous stage-structured competitive system with toxin liberating species and nontoxic species have not been investigated. Generally speaking, it takes some time for a species to reach maturity to produce the toxicant; then toxin liberating mature individual produces a substance toxic to the nontoxic mature individuals only. The inhibiting effect is zero in absence of either species, and the immature individual of each species is not affected by the toxicant. Furthermore, the species compete each other for the limited life resource within closed environment, but this competition only happens among the mature individuals and does not involve the immature individuals. Consequently, it is necessary to investigate the dynamic effect of stage structure and toxic effect on the population dynamics of two-species competitive system with toxin liberating species and nontoxic species.

The rest section of this paper is organized as follows: a nonautonomous two-species competitive model is established in the second section. Stage structure and maturation delay for each species are introduced, and toxic effect of toxin liberating species on nontoxic species is considered. In the third section, qualitative analyses are performed to investigate the effect of stage structure and toxic substances on the dynamical behavior of two-species competitive model system. The positivity and boundedness of solutions are analytically studied. By utilizing some comparison arguments, an iterative technique is proposed to discuss permanence of the species within competitive system. Furthermore, existence of positive periodic solutions is considered based on continuation theorem of coincidence degree theory. By constructing an appropriate Lyapunov functional, sufficient conditions for global stability of the unique positive periodic solution are analyzed. Numerical simulations are provided to support the theoretical findings obtained in this paper. Finally, this paper ends with a conclusion.

#### 2. Model Formulation

In this paper, the effect of stage structure and toxic substances on the dynamical behavior of two-species competitive model system is investigated under the following five hypotheses, which are given as follows.(H1)Two competing species, that is, nontoxic species and toxin liberating species, are considered in this paper. It is assumed that each species is divided into two-stage groups, and the immature and mature individuals are divided by a fixed period. and represent immature population density of nontoxic species and toxin liberating species at time , respectively; and denote mature population density of nontoxic species and toxin liberating species at time , respectively.(H2)-periodic continuous functions and denote the birth rate of immature population of nontoxic species and toxin liberating species at time , respectively. -periodic continuous functions and stand for the death rate of immature population of nontoxic species and toxin liberating species at time , respectively. -periodic continuous functions and are the death and overcrowding rate of mature population of nontoxic species and toxin liberating species at time , respectively.(H3) denotes time of immature nontoxic species to maturity. The term represents the immature nontoxic species which are born at time and survive at time with immature death rate. denotes time of immature toxin liberating species to maturity. The term represents the immature toxin liberating species which are born at time and survive at time with immature death rate.(H4)For toxin liberating species, it takes some time to attain its level of maturity to produce the toxic substances, and toxin liberating mature individual produces a substance toxic to the nontoxic mature individuals only. The inhibiting effect is zero in absence of either species, and the immature individual of each species is not affected by the toxicant. Based on model system (3), the toxic effect released by toxin liberating species on nontoxic species is described by the mathematical term , where toxic inhibition rate is represented by an -periodic continuous function .(H5)Nontoxic species and toxin liberating species compete each other for the common resource within closed environment, but this competition only happens among the mature individuals and does not involve the immature individuals. -periodic continuous function represents interspecific competition rate for the mature nontoxic species by the mature toxin liberating species, and -periodic continuous function represents interspecific competition rate for the mature toxin liberating species by the mature nontoxic species.

Based on hypotheses (H1)–(H5), a nonautonomous stage-structured competitive model with toxic effect and double maturation delays is established as follows:

In this paper, model system (4) is investigated with the following initial conditions:

For the continuity of the initial conditions, it is required that

#### 3. Qualitative Analysis of Model System

In this section, qualitative analysis of the nonautonomous model system (4) is performed, which is utilized to discuss dynamic effect of toxic effect and maturation delay on population dynamics. The positivity and boundedness of solutions are analytically studied. By utilizing some comparison arguments, an iterative technique is proposed to discuss permanence of the species within competitive system. Furthermore, existence of positive periodic solutions is investigated based on continuation theorem of coincidence degree theory. By constructing an appropriate Lyapunov functional, sufficient conditions for global stability of the unique positive periodic solution are analyzed.

Some mathematical notations are adopted for convenience of the following statement: where is a -periodic continuous function.

##### 3.1. Permanence of Solutions

Theorem 1. *Solutions of model system (4) with initial conditions (5) and (6) are positive for all .*

*Proof. *Firstly, we show that for all . Otherwise, if it is false, since for all , then it can be derived that there exists a such that .

Define . According to the definition of , it can be obtained that

It follows from the second equation of model system (4) that
and it is easy to show that , which is a contradiction to (8). Hence, for all .

By a direct computation, it follows from the first equation of model system (4) that
since for all ; it is easy to show that for all based on (10).

By utilizing the similar proof, it can be obtained that and for all . Consequently, solutions of model system (4) with initial conditions (5) and (6) are positive for all .

Theorem 2. *Solutions of model system (4) with initial conditions (5) and (6) are ultimately bounded.*

*Proof. *Let , where is an arbitrary positive solution of model system (4) with the initial conditions (5) and (6).

Calculating the derivative of along the solution of model system (4) gives that
where .

By using the standard comparison principle [20], it follows from (11) that
which implies that any solution of model system (4) with initial conditions (5) and (6) is ultimately bounded.

Lemma 3 (see [29]). *Consider the following differential equation:
**
where and for ; we have that*(i)*if , then ;*(ii)*if , then .*

*Lemma 4 (see [29]). Consider the following differential equation:
where and for ; we have
*

*Definition 5 (see [30]). *Model system
where , , . Model system (16) is said to be permanent if, for any solution , there exists a constant and such that for all .

*Definition 6 (see [30]). *The domain is said to be an ultimately bounded domain, if is a closed, bounded subset of , and there exists a constant such that for all .

*Theorem 7. If and , then model system (4) is permanent with initial conditions (5) and (6).*

*Proof. *According to the second equation of model system (4) and Theorem 1, we get that

By virtue of Lemma 4 and (17), there exists a positive time such that, for sufficiently small and , it yields

By rearranging (10), it can be obtained that

For any , it follows from (18) and (19) that

Based on the fourth equation of model system (4) and Theorem 1, it can be obtained that
holds for .

By virtue of Lemma 4 and (21), there exists such that, for sufficiently small and , it yields

By direct computing, it follows from the third equation of model system (4) that

For any , it follows from (22) and (23) that

Furthermore, it follows from the second equation of model system (4) that

Based on Lemma 3 and (25), there exists and for any and sufficiently small ,
holds provided that

According to (19), for any , we get that

For any , it follows from the fourth equation of model system (4) that

Based on Lemma 3 and (29), there exists and for any and sufficiently small ,
holds provided that

According to (23) and (30), it can be obtained that

According to the second equation of model system (4), we get that

By virtue of Lemma 3 and (33), there exists such that, for sufficiently small and , it yields
which holds provided that

For any , it follows from (19) and (34) that

Based on the fourth equation of model system (4), it can be obtained that
holds for .

By virtue of Lemma 3 and (37), there exists such that, for sufficiently small and , it yields
which holds provided that

For any , it follows from (23) and (38) that

Furthermore, it follows from the second equation of model system (4) that

Based on Lemma 3 and (41), there exists and, for any and sufficiently small ,
holds provided that

According to (19) and (42), for any , we get that

For any , it follows from the fourth equation of model system (4) that

Based on Lemma 3 and (45), there exists and, for any and sufficiently small ,
holds provided that

According to (23) and (46), it can be obtained that

By using simple computation, it is easy to show that six inequalities (27), (31), (35), (39), (43), and (47) hold if the following two inequalities and hold.

Furthermore, eight sequences will be obtained by repeating the discussion in this manner, which are given as follows:

It is easy to show that and the sequences () are decreasing as increases, which implies that exists; furthermore, it is easy to show that and the sequences () are increasing as increases, which implies that exists. Consequently, it follows from (49) that

Based on Definition 5 and (50), it can be concluded that model system (4) is persistent if and hold.

*3.2. Existence of Positive Periodic Solutions*

*3.2. Existence of Positive Periodic Solutions*

*Definition 8 (see [31]). *Let be a linear mapping and let be a continuous mapping, where and are real Banach spaces. If and is closed in , then is called a Fredholm mapping of index zero.

*If is Fredholm mapping of index zero and there exist continuous projectors and such that
then restriction of to is invertible.*

*Definition 9 (see [31]). *Denote the inverse of by . Supposing that is an open bounded subset of , if is bounded and is compact, then the mapping is called -compact on . Since is isomorphic to , there exists an isomorphism .

*Lemma 10 (see [31]). Let be an open bounded set, let be a Fredholm mapping of index zero, and let be -compact on . If the following three conditions hold: (i) for any and ,(ii) for any ,(iii),then has at least one solution in .*

*Theorem 11. If , , then model system (4) with initial conditions (5) and (6) has at least one positive -periodic solution.*

*Proof. *Consider the subsystem of model system (4):

Let , .

By substituting and into (52), it can be obtained that

It should be noted that if model system (53) has one -periodic solution , then is a positive -periodic solution of model system (52).

In order to utilize Lemma 10 in a straightforward manner, we define
where denotes the Euclidean norm; it is easy to show that both and are Banach spaces with the norm ; then define
where , , and

Furthermore, we define

According to the above definitions, it is not difficult to verify that , are closed in , , and both and are continuous projectors such that and .

Based on the above analysis, it can be obtained that is a Fredholm mapping of index zero.

Furthermore, the inverse of exists and takes the following form:

Hence, and can be defined as follows, respectively,

It is easy to show that and are continuous. In order to facilitate the proof based on Lemma 10, we also need to find an appropriate open and bounded subset , which can be found by the following two steps.*Step 1*. According to the operator equation for , the upper and lower bound of and will be estimated as follows:
where and have been defined in (56).

Suppose that is a solution of model system (60) for some . By integrating (60) over the interval , it can be obtained that

Based on definition , there exist such that

Multiplying the first equation of (60) by and integrating it over give that

It follows from (64) that

On the other hand, by using the inequality,

Based on (65) and (66), it can be obtained that
which derives that

It follows from (60) and (68) that

According to (68) and (69), it can be obtained that

Multiplying the second equation of (60) by and integrating it over give that

It follows from (71) that

Based on (66) and (72), it can be obtained that
which derives that

It follows from (60) and (74) that

According to (74) and (75), it can be obtained that

It should be noted that

Based on (64), it can be obtained that
which derives that
hold provided that .

According to (69) and (79), it can be obtained that

By virtue of (70) and (80), if , then

Similarly, it is easy to show that
which derives that
holds provided that .

According to (75) and (83), it is derived that

By virtue of (76) and (84), if , then

It is obvious that and in (81) and (85) are independent of .*Step 2*. In order to construct an appropriate open and bounded subset , denote , where is sufficiently large such that the unique solution of the following algebraic equations:
satisfies .

Select , which implies that condition (i) of Lemma 10 holds.

When , is a constant vector in with . Consequently, it can be concluded that
which implies that condition (ii) of Lemma 10 is satisfied.

Take , . It follows from straightforward computation that
where is the unique solution of (86). Hence, the condition (iii) of Lemma 10 holds.

Furthermore, it is easy to see that the set is equicontinuous and uniformly bounded. By using the Arzela-Ascoli theorem [31], it can be shown that is compact and is -compact.

Consequently, all conditions (i)–(iii) of Lemma 10 hold for . It follows from Lemma 10 that model system (53) has at least one -periodic solution , and model system (52) has at least one -periodic solution .

Let be a positive -periodic solution of model system (52); it follows from (19) and (23) that
are -periodic continuous function.

Based on the above analysis, if the following two inequalities hold:
then model system (4) with initial conditions (5) and (6) has at least one positive -periodic solution .

*3.3. Global Stability Analysis*

*3.3. Global Stability Analysis*

*Theorem 12. If , , then model system (4) with initial conditions (5) and (6) has a unique positive -periodic globally stable solution , where
and is a positive constant and , , , have been defined in (50).*

*Proof. *Suppose that is a positive -periodic solution of model system (4) with initial conditions (5) and (6).

Construct a Lyapunov functional as follows:

By calculating the upper right derivative of along the positive -periodic solution of model system (4), it can be obtained that