Abstract

A dengue disease epidemic model with nonlinear incidence is formulated and analyzed. The equilibria and threshold of the model are found. The stability of the system is analyzed through a geometric approach to stability. The proposed model also exhibits backward bifurcation under suitable conditions on parameters. Our results imply that a nonlinear incidence produces rich dynamics and they should be studied carefully in order to analyze the spread of disease more accurately. Finally, numerical simulations are presented to illustrate the analytical findings.

1. Introduction

Dengue disease is transmitted in humans by mosquitoes, called vectors, who carry the disease without getting it themselves. Every year many people die from dengue disease. It has become a global health problem due to its rapid spread. In order to predict the development extent of such epidemics, it is useful to find the suitable mathematical models, which can provide insight into the dynamics of a disease and can help us to make right decisions on public health policies.

However in the region of high infective prevalence such as Africa, some tropical countries, and so forth, there are lots of people infected with dengue, and so forth. and hence there are lots of infected mosquitoes. So in these endemic regions residents are frequently bitten by infected mosquitoes, thus making individuals exposed with multiple bites of mosquito which causes individuals to move to infected class with increased rate compared to people who are exposed with single effective bite of mosquito. So it is reasonable to consider nonlinear saturation incidence in place of bilinear incidence or standard incidence in the model as it also incorporates the effects of increased exposure to mosquito bites.

In fact the incidence rate is very important in the transmission dynamics of the disease and the qualitative behaviour of the disease changes with the change in the incidence rate. In a recent paper by Xiao and Tang [1], interaction of the nonlinear incidence and the partial immunity for a simple SIV epidemic model is investigated and it is shown that due to the nonlinear incidence, vaccination may contribute to disease spread rather than to its elimination. A vector-host epidemic model with nonlinear incidence is also analyzed by Cai and Li [2]. Here they considered nonlinear saturating type incidence and they studied the effect of this type of incidence rate on the basic reproduction number and the equilibrium level of the infective population. Hu et al. [3] analyzed an SIR epidemic model with nonlinear incidence rate and treatment. Here it is shown that reducing the basic reproduction number below one is not enough to eliminate the disease as backward bifurcation occurs for this model. Also when the basic reproduction number is greater than unity, there are multiple endemic equilibria. The conditions for the stability of endemic equilibrium and the existence of limit cycle are discussed in this paper. A discrete-time epidemic model with nonlinear incidence rates is studied by Li et al. [4], where a very complex dynamical behaviors are demonstrated. It is shown that there are possibilities of the transcritical bifurcation, flip bifurcation, Hopf bifurcation and chaos. A vector-borne model with nonlinear incidence is analyzed by Ozair et al. [5]. Here they considered the conditions for the stability of endemic equilibrium.

In recent years, dengue has become a major public health problem and is endemic in several countries. It is affecting almost one-third of the world's population. Dengue is an infection caused by a virus known as flavivirus and is spread by their vector Aedes aegypti. There are four different serotypes of flavivirus but infection with one serotype gives life-long immunity to that serotype but not to other three serotypes. So one can formulate dengue disease model as SIS or SIR model depending upon serotypes of flavivirus. In this paper we have formulated an SIR epidemic model with nonlinear incidence assuming the infection of dengue is caused by the flavivirus which is giving life-long immunity. In dengue-endemic area a person can be infected with more than one serotype of flavivirus in his/her lifetime. But in the current study we have ignored this fact as it will complicate the model and its analysis.

The paper is organized as follows. In Section 2 the model and some preliminary properties are presented. In Section 3 equilibria and bifurcation are studied. Section 4 is devoted to the global stability analysis of the endemic equilibrium, when it is unique. Numerical simulations are given in Section 5, we end the paper with conclusions in Section 6.

2. Model Formulation

The vector-host epidemic model with nonlinear incidence can be described by the following system of differential equations: The total host population size at time , given by , is partitioned into subclasses of individuals who are susceptible, infectious and recovered, with sizes denoted by , and , respectively. Furthermore, the host population dies at a natural death rate . In addition, the host population is recruited at a rate . We assume that vertical transmission in the host does not occur so that all newly recruited individuals are susceptible. The per capita recovery rate of the hosts is given by . The recovered individuals are assumed to acquire permanent immunity and there is no transfer from the class back to the class. The nonlinear incidence rate is , where and are positive constants, is the number of vectors at time who carry the pathogen. It corresponds to an increased rate of infection due to two exposures over a short time period. The single contacts lead to infection at the rate , whereas the new infective individuals arise from double exposures at a rate .

The second component of the vector population is the number of pathogen-free (susceptible) vectors at time , given by . The total size of the vector population at time , given by and is subdivided into these two vector-population classes, the susceptible vectors and infectious vectors. The vector population dies at a natural death rate . In addition, the vector population is recruited at a birth rate . Although there are evidences that the pathogen of several vector-borne diseases (e.g., West Nile fever, yellow fever and Lyme disease) can be transmitted from (female) parent to offspring in the vector population, we will assume that all newborn vectors are susceptible and vertical transmission is neglected. Susceptible vectors start carrying the pathogen after getting into contact (biting) an infective host at a rate so that the incidence of newly infected vectors is given by a mass-action term . In contrast to the host population, once the vectors become carriers of the microparasite, they carry it for life.

We make some reasonable technical assumptions on the parameters of the model, namely that ,   and for . The above systems for the host population and the vector are also equipped with initial conditions as follows: , , , and .

The total population size can be determined by or from the differential equation which is derived by adding the first three equations in (1). The total number of vectors can be determined by or from the differential equation which is derived by adding the last two equations in (1).

It is easily seen that both for the host population and for the vector population the corresponding total population sizes are asymptotically constant: So we can assume without loss of generality that ,   for all .

Then the dynamics of system (1) is qualitatively equivalent to the dynamics of system given by The values of and can be determined correspondingly from , and respectively.

For biological reasons we need the solutions nonnegative. Mathematical properties of the solutions lead us to study the system (5) in the closed set where denotes the non-negative cone of including its lower dimensional faces. It can be verified that is positively invariant with respect to (5). We denote by and the boundary and the interior of in , respectively.

3. Equilibria and Bifurcation Analysis

The disease-free equilibrium of system (5) is . The endemic equilibria of system (5) can be deduced by the system, which gives, and is the positive solution of the following equation, where Thus, we observe that By the Descartes’ rule of signs, it follows that there is a unique endemic equilibrium whenever and there are two endemic equilibria whenever , and .

Furthermore, we observe that there is a bifurcation point when , and . In fact the quantity can be expressed in terms of new threshold , where By a little algebraic calculation, we have whenever . Taking into account the above consideration, we have the following theorem.

Theorem 1. System (7) admits a unique endemic equilibrium when ; there are no endemic equilibria for ; there are two endemic equilibria, and , for and .

Next we analyze the stability of the disease-free equilibrium . The Jacobian matrix corresponding to (7) is First, we observe that so that the eigenvalues are given by the roots of the following qubic equation Hence if , then the eigenvalues are all negatives and is locally asymptotically stable. If , then two eigenvalues are negative and one is positive, so that is unstable.

Since there are no explicit expressions for equilibria and and the existence conditions for the two equilibria are very complex, it is very difficult to prove their stability by using the existing analytic methods.

Theorem 1 establishes that is a bifurcation value. In fact, across the disease free equilibrium changes its stability properties. Now, we investigate what kind of bifurcation occurs at . In order to do that, we will make use of the result summarized below, which has been obtained in [6] and is based on the use of the center manifold theory [7].

Let us consider a general system of ODEs with a parameter : Without loss of generality, we assume that is an equilibrium for (16).

Theorem 2 (see (Castillo-Chavez and Song [6])). Assume:(A1) is the linearization matrix of system (16) around the equilibrium with evaluated at . Zero is a simple eigenvalue of and all other eigenvalues of have negative real parts.(A2)Matrix has a right eigenvector and a left eigenvector corresponding to the zero eigenvalue.
Let denotes the th component of and, Then the local dynamics of system (16) around are totally determined by and (i), . When , with , is locally asymptotically stable and there exists a positive unstable equilibrium; when , is unstable and there exists a negative and locally asymptotically stable equilibrium;(ii), . When , with , is unstable; when , is locally asymptotically stable and there exists a positive unstable equilibrium;(iii), . When , with , is unstable and there exists a locally asymptotically stable negative equilibrium; when , is stable and a positive unstable equilibrium appears;(iv), . When changes from negative to positive, changes its stability from stable to unstable. Correspondingly, a negative unstable equilibrium becomes positive and locally asymptotically stable.

It clearly appears that, at a transcritical bifurcation takes place: more precisely, when and , such a bifurcation is forward; when and the bifurcation at is backward. We will apply Theorem 2 to show that system (7) may exhibit a backward bifurcation when . We consider the disease-free equilibrium and observe that the condition can be seen, in terms of the parameter , as . The eigenvalues of the matrix are given by Thus is a simple zero eigenvalue of the matrix and the other eigenvalues are real and negative. Hence, when (or equivalently when ), the disease-free equilibrium is a nonhyperbolic equilibrium: the assumption (A1) of Theorem 2 is then verified.

Now we denote by , a right eigenvector associated with the zero eigenvalue . It follows: so that Furthermore, the left eigenvector satisfying is given by: Then, the left eigenvector turns out to be: We can thus compute the coefficient and defined in Theorem 2, that is, Taking into account of system (5) and considering in and only the nonzero derivatives for the terms and , it follows that: In view of (21) and (23), we get Observe that the coefficient is always positive so that, according to Theorem 2, it is the sign of the coefficient which decides the local dynamics around the disease-free equilibrium for .

Let us introduce The coefficient is positive if and only if . In this case, the direction of the bifurcation of system (5) at is backward. So we summarize the above results in the following theorem.

Theorem 3. If , system (5) exhibits a backward bifurcation when . If , system (5) exhibits a forward bifurcation when .

The local dynamics in the neighborhood of the bifurcation value is described in the former case by the (i) of Theorem 2 and in the later case by the (iv) of the same theorem.

4. Global Stability Analysis

In Section 3 it has been shown that implies the existence and uniqueness of the endemic equilibrium . The stability analysis of will be here performed through the geometric approach to global stability due to Li and Muldowney [8]. The method has been summarized in the Appendix.

System (5), under the assumption , satisfies conditions (H1) and (H2) in Appendix. In fact, when , then is unstable. The instability of , together with , implies the uniform persistence, that is, there exists a constant such that: The uniform persistence is equivalent to the existence of a compact set in the interior of which is absorbing for (5). Thus, (H1) is verified. Moreover, as previously shown, is the only equilibrium in the interior of , so that (H2) is verified, too. Now it remains to find conditions for which the Bendixson criterion given by Theorem A.4 is verified. Taking into account of the Jacobian matrix (13), we obtain the second additive compound matrix , where We consider the following function So that Therefore where

Consider now the norm in as where denotes the vector in and denote by the Lozinski measure with respect to this norm. It follows where ,   are matrix norms with respect to the vector norm, and denotes the Lozinski measure with respect to the norm. More specifically, Therefore, Observe that system (5) provides the following equalities Substituting (40) into (38) and (40) into (39), respectively, we have Thus, (36) implies where In view of ,   and , we can deduce that if then where Hence and the Bendixson criterion given by Theorem A.4 is thus verified. The discussions above may be summarized as follows.

Theorem 4. If , then system (5) admits a unique endemic equilibrium . It is globally asymptotically stable with respect to solutions of (5) initiating in the interior of , provided that inequality (44) holds true.

5. Simulation

The system (5) is simulated for various sets of parameters using the package XPP [9]. In Figures 14, phase planes are drawn which confer the existence and the stability of different equilibria of the system (5). Here Figure 1 corresponds to the situation stated in Theorem 1, where and we get two endemic equilibria and and one disease free equilibria . And for , it is found that disease free equilibrium and the endemic equilibrium point with largest infective population are stable and the equilibrium point is unstable.

Figure 2 corresponds to the situation when and also and in this case we get unique stable endemic equilibrium point and unstable disease free equilibrium point . Figure 3 shows the stability of unique feasible disease free equilibrium for . For , and we get two feasible equilibria: the unstable disease free equilibria and a unique endemic equilibrium which is shown in Figure 4. This fact is more clear from the bifurcation Figure 5, which is obtained by considering as the critical parameter. The horizontal axis is labelled with the appropriate value of the reproduction number corresponding to this bifurcation parameter . For we get stable infection free equilibrium point and at we get forward bifurcation, that is, for we get unstable disease free equilibrium and unique stable endemic equilibrium point. Figure 6 is showing the backward bifurcation for the parameter values satisfying . This diagram too is obtained by considering as the critical parameter and horizontal axis is labelled with the appropriate value of the reproduction number corresponding to this bifurcation parameter . It is observed that when the reproduction number is between to , the infection free equilibrium alone is stable, for we have bistability where either the infection free equilibrium is stable or the equilibrium is stable. Again by increasing we get and we have unique stable endemic equilibrium point and unstable disease free equilibrium . The equilibrium when it exists is always saddle. As backward bifurcation is quite significant in this figure, so it seems that it can arise for biologically reasonable parameter values too.

6. Conclusions

In this paper, a dengue disease epidemic model with nonlinear incidence has been analyzed. If , the disease-free equilibrium is locally asymptotically stable, that is, the disease dies out. We find the condition under which the system exhibits backward bifurcation at . The stability of two endemic equilibria in the case of backward bifurcation for is not discussed analytically but we have shown these by numerical simulations. It is found that the endemic equilibrium with the largest infective population is stable, while the endemic equilibrium with less infective population is unstable. In the case of a unique endemic equilibrium, the global stability analysis has been performed. A generalization of the Poincaré-Bendixson criterion has been used. The sufficient conditions that we found are not completely satisfactory, as they involve the constant of uniform persistence . How to implement the geometric approach to give better stability conditions for epidemic models with convex incidence rate appears to be an open and stimulating challenge.

Appendix

Let be a function for in an open set . Consider the differential equation Denote by the solution of (A.1) such that . We take the following two assumptions.(H1) There exits a compact absorbing set .(H2)Equation (A.1) has a unique equilibrium in .

The equilibrium is said to be globally stable in if it is locally stable and all trajectories in converge to . The following global-stability problem is formulated in [8].

Theorem A.1. Under the assumptions (H1) and (H2), the global stability of with respect to is implied by its local stability.

Assumptions (H1) and (H2) are satisfied if is global stability in . For , a Bendixson criterion is a condition satisfied by which precludes the existence of nonconstant periodic solution of (A.1). A Bendixson criterion is said to be robust under local perturbations of at if, for sufficiently small and neighborhood of , it is also satisfied by such that the support and , where Such will be called local -perturbations of    at . It is easy to see that the classic Bendixson’s condition for is robust under local perturbations of at each . Bendixson’s condition for higher dimensional system that are robust are discussed in [8, 10, 11].

A point is wandering for (A.1) if there exists a neighborhood of and such that is empty for all . Thus, for example, all equilibria and limit points are nonwandering. The following is a version of the local Closing lemma of Pugh (see [12, 13]).

Lemma A.2. Let . Assume that all positive semitrajectories of (A.1) are bounded. Suppose that is a nonwandering point of (A.1) and that . Then, for each neighborhood    of and , there exists a local- perturbation    of    at such that(a) and(b)the perturbed system has a nonconstant periodic solution whose trajectory passes through .

The following general global-stability principle is established in [8].

Theorem A.3. Suppose that assumptions (H1) and (H2) hold. Assume that (A.1) satisfies a Bendixson criterion that is robust under local perturbations of at all nonequilibrium nonwandering points for (A.1). Then is globally stable in provided it is stable.

The main idea of the proof in [8] for Theorem A.3 is as follows: suppose that system (A.1) satisfies a Bendixson criterion. Then it does not have any nonconstant periodic solutions. Moreover, the robustness of the Bendixson criterion implies that all nearby differential equations have nonconstant periodic solutions. Thus by Lemma A.2, all of nonwandering points of (A.1) in must be equilibria. In particular, each omega limit point in must be an equilibrium. Therefore for all since is the only equilibrium in .

A method of deriving a Bendixson criterion in is developed in [14]. The idea is to show that second compound equation with respect to a solution to (A.1), is uniformly asymptotically stable, and the exponential decay rate of all solutions to (A.3) is uniform for in each compact subset of . Here is second additive compound matrix of the Jacobian matrix . Generally speaking, for an matrix , is a matrix and in the special case , one has

is a matrix, and thus (A.3) is a linear system of dimension . If is simply connected, the above mentioned stability property of (A.3) implies the exponential decay of the surface area of any compact surface in , which in turn precludes the existence of any invariant simply closed rectifiable curve in , including periodic orbits. The required uniform asymptotic stability of system (A.3) can be proved by constructing a suitable Lyapunov function.

Let matrix-valued function that is for . Assume that exists and is continuous for , the compact absorbing set. A quantity is defined as where and the matrix is obtained by replacing each entry of by its derivative in the direction of , . The quantity is the Lozinski measure of with respect to a vector norm in , , defined by see [14, p. 41]. It is shown in [8] that, if is simply connected, the condition rules out the presence of any orbit that gives rise to a simple closed rectifiable curve that is invariant for (A.1), such as periodic orbits, homoclinic orbits, and heteroclinic cycles. Moreover, it is robust under local perturbations of near any nonequilibrium point that is nonwandering. In particular, the following global stability result is proved in Theorem 3.5 of [8].

Theorem A.4. Assume that is simply connected and that assumptions (H1), (H2) hold. Then the unique equilibrium of (A.1) is globally stable in if .

Acknowledgments

The authors are very grateful to two referees for their valuable comments and suggestions which led to an improvement of their original paper. They also thank the Handling Editor and all the editorial staff for the help and support provided to them. This work is supported by the NSF of China (no. 11271314) and NSF of Henan Province (no. 112300410058).