# A Malaria Model with Two Delays

Authors: Hui Wan; Jing-an Cui
Journal: Discrete Dynamics in Nature and Society
 (2013)

---

## Abstract

A transmission model of malaria with two delays is formulated. We calculate the basic reproduction numberR0 for the model. It is shown that the basic reproduction number is a decreasing function of two time delays. The existence of the equilibria is studied. Our results suggest that the model undergoes a backward bifurcation, which implies that bringing the basic reproduction number below 1 is not enough to eradicate malaria.


---

## Body

## 1. Introduction

Malaria is a mosquito-borne disease which is due to the four species of the genus Plasmodium, namely,Plasmodium falciparum, Plasmodium vivax, Plasmodium malariae, and Plasmodium ovale. These parasites are transmitted to the human host through a bite by an infected female anopheles mosquito.Sporozoites are injected into a human host, which are carried through the blood to the liver within 30 min [1]. They invade hepatocytes and undergo a process of asexual replication (exoerythrocytic schizogony) to give rise to 10–40 thousand merozoites per sporozoite. Up to this point, the infection is nonpathogenic and clinically silent. After about 7–9 days, the liver schizonts rupture to release the merozoites into the blood and clinical symptoms, such as fever, pain, chills, and sweats, may develop. Each merozoite invades an erythrocyte and divides to form an erythrocytic schizont containing about 16 daughter merozoites [2]. These merozoites either reinfect fresh erythrocytes, giving rise to cyclical blood-stage infection with a periodicity of 48–72 h, depending on the Plasmodium species, or differentiate into sexual transmission stages called gametocytes. When a second mosquito bites the infected human, the gametocytes are ingested, giving rise to extracellular gametes. In the mosquito midgut, the gametes fuse to form a motile zygote (ookinete), which penetrates the mid-gut wall and forms an oocyst, within which meiosis takes place and haploid sporozoites develop. These sporozoites migrate to the salivary glands. The incubation period within the mosquito may last 8–22 days [3]. The variation in the length of time is due to the environmental temperature. For P. falciparum, the average time is 12 days [4, 5]. Malaria can also be transmitted through blood transfusion, organ transplantation and transplacental malaria (i.e., congenital malaria) can also be significant in populations which are partially immune to malaria.The application of mathematics to study of infectious disease appear to have been initiated by Berniulli, 1760. Presently, a lot of mathematical models have been proposed to evaluate and compare control procedures and preventive strategies, and to investigate the relative effects of various sociological, biological, and environmental factors on the spread of diseases ([6–14], etc.). These models have played a very important role in the history and development of epidemiology.The transmission process involves considerable time delay both in human and in mosquitoes due to the incubation periods of the several forms of the parasites [9]. Several works have been done to consider the effect of delay. Taking account of the incubation period in vector (mosquito), Takeuchi et al. [10] considered a differential-delay model for vector-borne diseases and get the global stability of the endemic equilibrium under appropriate conditions. While Wei et al. [15] proposed a similar model of vector-borne disease which has direct mode of transmission in addition to the vector-mediated transmission and proved that the introduction of a time delay in the host-to-vector transmission term can destabilize the system and periodic solutions can arise through Hopf bifurcation. In [9], Ruan et al. proposed a delayed Ross-Macdonald model with delays both in human and in mosquitoes. The effect of time delays on the basic reproduction number and the dynamics of the transmission was studied.Both the models in [10, 15] adopted mass action formulation, where the contact rate depends on the size of the total host population while the model in [9] adopted standard incidence formulation, where the contact rate is assumed to be constant (see [16] for detailed derivation of these incidence functions). The choice of one formulation over the other really depends on the disease being modeled and, in some cases, the need for mathematical (analytical) tractability [17]. A crucial question is which of the two incidence functions is more suitable for modeling infectious disease in general, and malaria in particular? As noted in [16], data ([4, 18]) strongly suggests that standard incidence is better for modeling human diseases than mass action incidence. Consequently, standard incidence rate is adopted in this paper.


## 2. Derivation of the Model

In order to study the impact of incubation periods in both human and mosquitoes on the basic reproduction number and the transmission dynamics of malaria over long periods, we propose a model based on SIRS in human population and SI for the mosquito vector population. Since the mosquito dynamics operates on a much faster time scale than the human dynamics and the turnover of the mosquito population is very high and the total size of vector population is largely exceeding the human total size, the mosquito population can be considered to be at a equilibrium with respect to changes in the human population. Hence, the total number of mosquito population is assumed to be constant.We formulate an SIRS-SI model for the spread of malaria in the human and mosquito population with the total population size at timet is given by Nh(t) and Nv(t), correspondingly. For the human population, the three compartments represent individuals who are susceptible, infectious, and partially immune, with sizes at time t denoted by Sh(t), Ih(t), and Rh(t), respectively. The incubation period in a human has duration τ1≥0. A susceptible individual infected by a mosquito will not become infectious until τ1 time units later and we assume that no one recovers during this incubation period. At the infectious stage a host may die from the disease, recover into the susceptible class, or may recover with acquired partial immunity. Infectious human hosts are those with infectious gametocytes in the blood stream. Partially immune hosts still have protective antibodies and other immune effectors at low levels. Λh>0 is the human input (birth) rate. μh>0 and αh≥0 are the natural and disease-induced death rates, respectively. γh>0 is the rate at which human hosts acquire immunity. θ≥0 is the per capita rate of recovery into the susceptible class from being infectious. ρ≥0 is the per capita rate of loss of immunity in human hosts. βh>0 is the proportion of bites on man that produce an infection. For the mosquito population, the two compartments represent susceptible and infectious mosquitoes, with sizes at time t denoted by Sv(t) and Iv(t), respectively. The vector component of the model does not include immune class as mosquitoes never recover from infection that is, their infective period ends with their death due to their relatively short life cycle. Thus the immune class in the mosquito population is negligible and death occurs equally in all groups. Our model also excludes the immature mosquitoes since they do not participate in the infection cycle and are, thus, in the waiting period, which limits the vector population growth. We assume that the mosquito population is constant and equal to Nv, with birth and death rate constants equal to μv. βv is the probability that a mosquito becomes infectious when it bites an infectious human. The biting rate b of mosquitoes is the average number of bites per mosquito per day. This rate depends on a number of factors, in particular, climatic ones, but for simplicity in this paper, we assume b constant. The incubation interval in the mosquito has duration τ2≥0.Considering the assumption made above, the interaction between human hosts and the mosquito vector population with standard incidence rate is described as shown below:(1)dShdt=Λh-bβhNh(t-τ1)Iv(t-τ1)Sh(t-τ1)e-μhτ1-μhSh+θIh+ρRh,dIhdt=bβhNh(t-τ1)Iv(t-τ1)Sh(t-τ1)e-μhτ1-(γh+μh+αh+θ)Ih,dRhdt=γhIh-μhRh-ρRh,dSvdt=μvSv+μvIv-bβvNh(t-τ2)Ih(t-τ2)Sv(t-τ2)e-μvτ2-μvSv,dIvdt=bβvNh(t-τ2)Ih(t-τ2)Sv(t-τ2)e-μvτ2-μvIv,dNhdt=Λh-μhNh-αhIh,
where Sh(t), Ih(t), and Rh(t) denote the number of susceptible, infective, and recovered in the human population, and Sv(t), Iv(t) denote the number of susceptible and infective in the mosquito vector population. In addition, Sh+Ih+Rh=Nh and Sv+Iv=Nv.System (1) satisfies the initial conditions: Sh(θ)=Sh0(θ), Ih(θ)=Ih0(θ), Rh(θ)=Rh0(θ), Sv(θ)=Sv0(θ), Iv(θ)=Iv0(θ), θ∈[-τ,0], where τ=max{τ1,τ2}. The nonnegative cone of R6 space ShIhRhSvIvNh is positively invariant for system (1) since the vector field on the boundary does not point to the exterior. What is more, since dNh/dt<0 for Nh>Λh/μh and Nv is constant, all trajectories in the first quadrant enter or stay inside the region
(2)Z+={Sh+Ih+Rh=Nh⩽Λhμh,Sv+Iv=Nv}.
The continuity of the right-hand side of (1) implies that unique solutions exist on a maximal interval. Since solutions approach, enter or stay in D+, they are eventually bounded and hence exist for t>0. Therefore, the initial value problem for system (1) is mathematically well posed and biologically reasonable since all variables remain nonnegative.In order to reduce the number of parameters and simplify system (1), we normalize the human and mosquito vector population sh=Sh/(Λh/μh), ih=Ih/(Λh/μh), rh=Rh/(Λh/μh), nh=Nh/(Λh/μh), sv= Sv/Nv, iv=Iv/Nv, and m=Nv/(Λh/μh).Sincerh=nh-sh-ih and sv=1-iv, the dynamics of system (1) is qualitatively equivalent to the dynamics of system given by
(3)dshdt=μh-bβhmnh(t-τ1)iv(t-τ1)sh(t-τ1)e-μhτ1-μhsh+θih+ρ(nh-sh-ih),dihdt=bβhmnh(t-τ1)iv(t-τ1)sh(t-τ1)e-μhτ1-(γh+μh+αh+θ)ih,divdt=bβvnh(t-τ2)ih(t-τ2)(1-iv(t-τ2))e-μvτ2-μviv,dnhdt=μh-μhnh-αhih,
and all trajectories in the nonnegative cone R+4 enter or stay inside the region D={0⩽sh, 0⩽ih, sh+ih⩽nh⩽1, 0⩽iv⩽1}.Define the basic reproduction number by(4)R0=mb2βhβve-μhτ1e-μvτ2μv(γh+μh+αh+θ).
It represents the expected number of secondary cases caused by a single infected individual introduced into an otherwise susceptible population of hosts and vectors. Take a primary case of host with a natural death rate of μh, the average time spent in an infectious state is 1/(γh+μh+αh+θ). Since the incubation period in mosquitoes has duration τ2, during which period some mosquitoes may die, the average number of mosquito bites received from m susceptible mosquitoes each with a biting rate b gives a total of bβvme-μvτ2/(γh+μh+αh+θ) mosquitoes successfully infected by the primary human case. Each of these mosquitoes survives for an average time 1/μv. Because of another incubation period τ1 in human, during which period some individuals may die, totally bβhe-μhτ1/μv hosts will be infected by a single mosquito. Therefore, the total number of secondary cases is thus mb2βhβve-μhτ1e-μvτ2/ μv(γh+μh+αh+θ), which is the basic reproduction number R0.


## 3. Existence of Equilibria

Equating the derivatives on the left-hand side to zero and solving the resulting algebraic equations. The points of equilibriumP^(s^h,i^h,i^v,n^h) satisfy the following relations:
(5)s^h=μh+ρ-Hi^hμh+ρ,n^h=μh-αhi^hμh,i^v=bβvμhi^he-μvτ2(bβvμhe-μvτ2-αhμv)i^h+μvμh,
where H=γh+μh+αh+ρ+ρ(αh/μh). Substituting (5) in the corresponding second equilibrium equation of (3), we obtain that the solutions are i^h=0, and the zero points of the quadratic polynomial
(6)r(ih)=a1ih2+a2ih+a3,
where
(7)a1=(bβvμhe-μvτ2-αhμv)αhμh,a2=2αhμv-bβvμhe-μvτ2-mb2βhβvμhH(μh+ρ)(γh+μh+αh+θ)e-μhτ1e-μvτ2,a3=μhμv(R0-1),R0=mb2βhβve-μhτ1e-μvτ2μv(γh+μh+αh+θ).The solutioni^h=0 gives the disease-free equilibrium point P0(1,0,0,1). We are looking for nontrivial equilibrium solutions in the interior of Ω. From (5), in order to keep s^h and n^h positive, i^h∈(0,(μh+ρ)/H) must be satisfied. Evaluating r(ih) at the end points of the interval, we obtain
(8)r(0)=μhμv(R0-1),r(μh+ρH)<0.WhenR0>1, then r(0)>0; therefore, there exists a unique root in the interval (0,(μh+ρ)/H), which implies the existence of a unique equilibrium point P*(s^h*,i^h*,i^v*,n^h*), where i^h* satisfies: if a1=0,
(9)i^h*=-a3a2.
If a1≠0,
(10)i^h*=-a2-Δ2a1,
where Δ=a22-4a1a3.WhenR0=1, the roots of (6) are 0 and -a2/a1. It is easy to see that there exists a unique root in the interval (0,(μh+ρ)/H) if and only if a2>0, which implies the existence of a unique equilibrium point P*.WhenR0<1, then r(0)<0. The conditions to have at least one root in the mentioned interval are a1<0, 0<-a2/2a1<(μh+ρ)/H and Δ=a22-4a1a3⩾0. If Δ>0, there exist two roots in the interval mentioned above, which implies the existence of two equilibria P1(s^h1,i^h1,i^v1,n^h1) and P2(s^h2,i^h2,i^v2,n^h2), where
(11)i^h1=-a2+Δ2a1,i^h2=-a2-Δ2a1.IfΔ=0, there is a unique root in the interval mentioned above, which implies the existence of unique equilibrium.Then, we can conclude the above results in the following theorem.Theorem 1.
(1) The disease-free equilibriumP0(1,0,0,1) always exists.
(2) IfR0>1, there exists a unique positive equilibrium P*(s^h*,i^h*,i^v*,n^h*).
(3) IfR0=1, then there is a positive equilibrium P*(s^h*,i^h*,i^v*,n^h*) when a2>0; otherwise, there is no positive equilibrium.
(4) IfR0<1, then (a) if a1⩾0, there is no positive equilibrium; (b) if a1<0, the system (3) has two positive equilibria P1(s^h1,i^h1,i^v1,n^h1) and P2(s^h2,i^h2,i^v2,n^h2) if and only if Δ>0 and 0<-a2/2a1<(μh+ρ)/H, and these two equilibria merge with each other if and only if 0<-a2/2a1<(μh+ρ)/H and Δ=0; otherwise, there is no positive equilibrium.


## 4. Disease-Free Equilibrium

It is easy to see the disease-free equilibriumP0(1,0,0,1) always exists. In this section, we study the stability of the disease-free equilibrium P0.Theorem 2.
P
0 is locally stable if R0<1. P0 is a degenerate equilibrium of codimension one and is stable except in one direction if R0=1. P0 is unstable if R0>1.Proof.
Linearizing the system (3) at P0. The eigenvalues of the Jacobian matrix J(P0) are -μh-ρ, -μh and the solutions of the following transcendental equation:
(12)(λ+μv)(γh+μh+αh+θ)-mb2βhβve-μhτ1-μvτ2e-(τ1+τ2)λ=0.
Let
(13)F(λ,τ1,τ2)=(λ+μv)(γh+μh+αh+θ)-mb2βhβve-μhτ1-μvτ2e-(τ1+τ2)λ.
Clearly, F(λ,τ1,τ2) is an analytic function. F(0,τ1,τ2)=μv(γh+μh+αh+θ)(1-R0), and F(λ,0,0)= (λ+μv)(γh+μh+αh+θ)-mb2βhβv. In the following, we discuss the distribution of the solutions of (12) in three cases.(i)
R0<1, then F(0,τ1,τ2)>0. Since Fλ′(λ,τ1,τ2)>0 for λ≥0, τ1>0 and τ2>0, (12) has no zero root and positive real roots for all positive τ1 and τ2. Now we claim that (12) does not have any purely imaginary roots. Suppose that (12) has a pair of purely imaginary roots ωi, ω>0 for some τ1 and τ2. Then, by calculation, ω must be a positive real root of
(14)ω4+[(γh+μh+αh+θ)2+μv2]ω2+(γh+μh+αh+θ)2μv2-(mb2βhβve-μhτ1-μvτ2)2=0.
However, it is easy to see that (14) does not have nonnegative real roots when R0<1. Hence, (12) does not have any purely imaginary roots. On the other hand, one can easily get that the roots of F(λ,0,0)=0 all have negative real parts when R0<1. By the implicit function theorem and the continuity of F(λ,τ1,τ2), we know that all roots of (12) have negative real parts for positive τ1 and τ2, which implies that P0 is stable.(ii)
R0=1, then F(0,τ1,τ2)=0. Since Fλ′(λ,τ1,τ2)>0 for λ≥0, τ1>0 and τ2>0, (12) has a simple zero root and no positive root for all positive τ1 and τ2. Using a similar argument as in (i), we can obtain that except a zero root, all roots of (12) have negative real parts for positive τ1 and τ2. Thus, P0 is a degenerate equilibrium of codimension one and is stable except in one direction.(iii)
R0>1, then F(0,τ1,τ2)<0. Since limλ→∞F(λ,τ1,τ2)=∞ and Fλ′(λ,τ1,τ2)>0 for λ≥0, τ1>0 and τ2>0, (12) has a unique positive real root for all positive τ1 and τ2 and P0 is unstable.Theorem 3.
IfR0≤1, the disease-free equilibrium P0 is globally asymptotically stable when α=0.Proof.
We denote byxt the translation of the solution of the system (3), that is, xt=(sh(t+θ), ih(t+θ), iv(t+θ)), nh(t+θ) where θ∈[-τ,0], τ=max{τ1,τ2}. In order to prove the globally stability of P0, we define a Lyapunov function:
(15)L(xt)=eμhτ1ih+mbβhμveμvτ2iv(t)+mbβh∫t-τ1tiv(t)sh(t)nh(t)dt+mb2βhβvμv∫t-τ2tih(t)1-iv(t)nh(t)dt.
The derivative of L along the solutions of (3) is given by
(16)dL(xt)dt=-(μh+γh+θ)eμhτ1ih(t)(1-R0e-μhτ11-iv(t)nh(t))-mbβheμvτ2iv(t)(1-sh(t)nh(t)e-μvτ2).
Since limt→∞nh=1 and R0≤1, we can easily get that dF/dt≤0 and the set {P0} is the largest invariant set within the set where dL(xt)/dt=0. By LaSalle-Lyapunov Theorem, all trajectories approach P0 as t→∞.Usually, the disease can be controlled if the basic reproduction number is smaller than one. Nevertheless, Theorem1 indicates the possibility of backward bifurcation (where the locally asymptotically stable DFE coexists with endemic equilibrium when R0<1) when a1<0 for the model (3) and it becomes impossible to control the disease by just reducing the basic reproduction number below one. Therefore a further threshold condition beyond the basic reproduction is essential for the control of the spread of malaria.To check for this, the discriminantΔ is set to zero and solved for the critical value of R0, denoted by R^0, given by
(17)R^0=1+a224a1μhμv.
Thus, backward bifurcation would occur for values of R0 such that R^0<R0<1. This is illustrated by simulating the model with the following set of parameter values (it should be stated that these parameters are chosen for illustrative purpose only, and may not necessarily be realistic epidemiologically).Letμh=1/22000, μv=1/14, βh=0.3, βv=0.95, m=0.5, αh=0.05, γh=0.26, b=0.5, ρ=0.01, θ=0.08, τ1=15. With this set of parameters, R^0=0.093970<1. When τ2=10, R0=0.62555<1, so that R^0<R0<1. There are two positive roots for (6), which are i^h1=.3501527770e-3 and i^h2=.8864197396e-3.Thus, the following result is established.Lemma 4.
The model (3) undergoes backward bifurcation when a1<0 and 0<-a2/2a1<(μh+ρ)/H holds and R^0<R0<1.Epidemiologically, because of the existence of backward bifurcation, whether malaria will prevail or not depends on the initial states.R0<1 is no longer sufficient for disease elimination. In order to eradicate the disease, we have a subthreshold number R^0 and it is now necessary to reduce R0 to a value less than R^0.The existence of the subthreshold conditionR^0 also has some implications. In some areas where malaria dies out (R0<R^0), it is possible for the disease to reestablish itself in the population because of a small change in environmental or control variables which increase the basic reproduction number above. On the other hand, in some areas with endemic malaria, it may be possible to eradicate the disease with small increasing in control programs such that the basic reproduction number less than R^0.From the expression ofR0, we can see that the basic reproduction number is a decreasing function of both time delays. The time delays have important effect on the basic reproduction number. Especially, since global warming decreases the duration of incubation period and R0 increases with the decreasing of τ1 or τ2, global warming will affect the transmission of malaria a lot.From Theorem3, if there is no disease-induced death, backward bifurcation will not occur and there is a unique endemic equilibrium when R0>1. Therefore, the disease-induced death is the cause of the backward bifurcation phenomenon in the model (3).


## 5. Endemic Equilibria

To obtain precise results, we first assume that malaria does not produce significant mortality (α=0). The assumption is not justifiable in all regions where malaria is endemic but it is a useful first approximation. In this part, we will determine the stability of the unique endemic equilibria P* when R0>1.Linearizing the system (3) at P*. The eigenvalues of the Jacobian matrix J(P*) are -μh and the solutions of the following transcendental equation:
(18)λ3+A1λ2+A2λ+A3=0,
where
(19)A1=2μh+γh+ρ+θ+μv+mbβhi^ve-μhτ1e-λτ1+bβvi^he-μvτ2e-λτ2,A2=(μh+ρ)(γh+μh+θ)+(μh+ρ)μv+(γh+μh+θ)μv+(μh+ρ+γh+μh+θ)bβvi^he-μvτ2e-λτ2-μv(γh+μh+θ)e-(τ1+τ2)λ+(γh+μh+ρ+μv+bβvi^he-μvτ2e-λτ2)×mbβhi^ve-μhτ1e-λτ1,A3=(μh+ρ)(γh+μh+θ)μv+(μh+ρ)(γh+μh+θ)bβvi^he-μvτ2e-λτ2-(μh+ρ)(γh+μh+θ)μve-(τ1+τ2)λ+(γh+μh+ρ)mbβhi^ve-μhτ1e-λτ1×(μv+bβvi^he-μvτ2e-λτ2).
For any nonnegative τ1 and τ2, we have the following proposition.Proposition 5.
For the endemic equilibriumP*(s^h*,i^h*,i^v*,n^h*) of the system with characteristic equation (18), one always has
(20)A1>0,A2>0,A3>0,A1A2-A3>0.Proof.
It is clear thatA1>0,  A2>0 and A3>0. For A1A2-A3, we have
(21)A1A2-A3=(γh+μh+ρ)×{e-(τ1+τ2)λ(μh+ρ)(γh+μh+θ)+(μh+ρ)μv+(γh+μh+θ)μv+(μh+ρ+γh+μh+θ)bβvi^he-μvτ2e-λτ2-μv(γh+μh+θ)e-(τ1+τ2)λ+(γh+μh+ρ+μv)mbβhi^ve-μhτ1e-λτ1}+(mbβhi^ve-μhτ1e-λτ1+μh+θ)A2+(μv+bβvi^he-μvτ2e-λτ2)×{e-(τ1+τ2)λ(μh+ρ)μv+(γh+μh+θ)μv+(μh+ρ+γh+μh+θ)bβvi^he-μvτ2e-λτ2+(γh+μh+ρ+μv+bβvi^he-μvτ2e-λτ2)×mbβhi^ve-μhτ1e-λτ1-μv(γh+μh+θ)e-(τ1+τ2)λ}+(μh+ρ)(γh+μh+θ)μve-(τ1+τ2)λ>0.Case 1.
Whenτ1=τ2=0, as a result of Proposition 5 and Hurwitz criterion, all roots of the characteristic equation (18) have negative real parts and the endemic equilibrium P* of (3) is stable when τ1=τ2=0.Case 2.
Whenτ1>0,  τ2=0, the characteristic equation (18) becomes
(22)λ3+A11λ2+A21λ+A31=e-λτ1(T11λ2+T21λ+T31),
where
(23)A11=2μh+γh+ρ+θ+μv+bβvi^h,A21=(μh+ρ)(γh+μh+θ)+(μh+ρ)μv+(γh+μh+θ)μv+(μh+ρ+γh+μh+θ)bβvi^h,A31=(μh+ρ)(γh+μh+θ)μv+(μh+ρ)(γh+μh+θ)bβvi^h,T11=-mbβhi^ve-μhτ1,T21=(γh+μh+θ)μv-(γh+μh+ρ+μv+bβvi^h)mbβhi^ve-μhτ1,T31=(μh+ρ)(γh+μh+θ)μv-(γh+μh+ρ)(μv+bβvi^h)mbβhi^ve-μhτ1.
By the implicit function theorem and the continuity of the left-hand side function of (18), all roots of (22) have negative real parts for small τ1. Notice that the condition R0>1 is equivalent to
(24)τ1<τ1*=1μhlnmb2βhβvμv(γh+μh+θ).
Furthermore, we claim that (22) does not have any nonnegative real roots for any τ1>0. Rewrite (22) by moving the positive terms from the right-hand side to the left-hand side. The rewritten (22) takes the form
(25)λ3+A11λ2+A~21λ+A~31=e-λτ1(T11λ2+T~21λ+T~31),
where A~21=A21-(γh+μh+θ)μve-λτ1, A~31=A31-(μh+ρ)(γh+μh+θ)μve-λτ1. It is easy to see A~21>0 and A~31>0 for all λ≥0 and τ1∈(0,τ1*). Consequently, the left-hand side in (25) is positive for all λ≥0 while the right-hand side is negative for all λ≥0 and the two cannot be equal for any λ≥0. Therefore, (22) does not have any nonnegative real roots for any τ1∈(0,τ1*). Now we want to show that all roots of (22) have negative real parts for all τ1∈(0,τ1*). To do so, we show that (22) does not have any purely imaginary roots for all τ1∈(0,τ1*).
We assume thatλ=iω with ω>0 being a root of (22). Then ω must satisfy the following system:
(26)A21ω-ω3=T21ωcos(ωτ1)-(T31-T11ω2)sin(ωτ1),A31-A11ω2=(T31-T11ω2)cos(ωτ1)+T21ωsin(ωτ1).
Thus, ω must be a positive root of
(27)ω6+B1ω4+B2ω2+B3=0,
where
(28)B1=A112-2A21-T112,B2=A21-2A11A31+2T11T31-T212,B3=A312-T312.
Let z=ω2, then (27) becomes
(29)z3+B1z2+B2z+B3=0.
Clearly, IfB1≥0, B2≥0 and B3≥0, then (29) has no positive real roots. Therefore, (22) does not have any purely imaginary roots for all τ1∈(0,τ1*) so that all roots of the characteristic equation (22) have negative real parts and the endemic equilibrium P* of (3) is stable.Case 3.
Whenτ2>0,  τ1=0, the characteristic equation (18) becomes
(30)λ3+A12λ2+A22λ+A32=e-λτ2(T12λ2+T22λ+T32),
where
(31)A12=2μh+γh+ρ+θ+μv+mbβhi^v,A22=(μh+ρ)(γh+μh+θ)+(μh+ρ)μv+(γh+μh+θ)μv+(μh+ρ+γh+μv)mbβhi^v,A32=(μh+ρ)(γh+μh+θ)μv+μv(γh+μh+ρ)mbβhi^v,T12=-bβvi^he-μvτ2,T22=(γh+μh+θ)μv-(2μh+γh+ρ+θ+mbβhi^v)bβvi^he-μvτ2,T32=(μh+ρ)(γh+μh+θ)μv-(γh+μh+ρ)mb2βhβvi^hi^ve-μvτ2-(μh+ρ)(γh+μh+θ)bβvi^he-μvτ2.
The condition R0>1 is equivalent to
(32)τ2<τ2*=1μvlnmb2βhβvμv(γh+μh+θ).
Using a similar argument as in Case 2, we know that all roots of (30) have negative real parts for τ2∈(0,τ2*) when C1≥0,  C2≥0, and C3≥0, where
(33)C1=A122-2A22-T122,C2=A22-2A12A32+2T12T32-T222,C3=A322-T322.Case 4.
Whenτ2>0,  τ1>0, the condition R0>1 is equivalent to
(34)τ2<τ2*(τ1)=1μvlnmb2βhβve-μhτ1μv(γh+μh+θ).
From Cases 1, and 2 the roots of (18) only have negative real parts for τ1∈[0,τ1*) and τ2=0 under the condition of Bi≥0, i=1,2,3. By the implicit function theorem and the continuity of the left-hand side function of (18), there is a τ¯2(τ1) satisfying 0<τ¯2(τ1)≤τ2*(τ1), such that all roots of (18) have negative real parts for 0<τ2<τ¯2(τ1). We show that τ¯2(τ1)=τ2*(τ1) when Bi≥0 and Ci≥0, i=1,2,3. Suppose τ¯2(τ1)<τ2*(τ1) for τ1∈[0,τ1*), then there must be a τ~2(τ1), τ¯2(τ1)<τ~2(τ1)<τ2*(τ1), such that one root of (18) has nonnegative real part for τ2=τ~2(τ1). As a result of the continuity of τ2 in τ1, we have τ~2(0)<τ2*(0)=τ2*. However, from the argument in Case 3, we know that all the roots of (18) have negative real parts for τ1=0 and τ2∈(0,τ2*). Contradict. Thus, τ¯2(τ1)=τ2*(τ1) which implies that the endemic equilibrium P* is stable when τ1∈[0,τ1*), τ2∈[0,τ2*(τ1)), Bi≥0, and Ci≥0, i=1,2,3. The above analysis can be summarized into the following theorem.Theorem 6.
IfR0>1, Bi≥0 and Ci≥0, i=1,2,3, the unique endemic equilibrium P* of system (3) is stable.


## 6. Discussions

In this paper, in order to study the impact of the incubation periods in human and mosquitoes, we formulate a model with two delays for the transmission dynamics of malaria.Existence of equilibria is obtained under different conditions and their stabilities are analyzed too. We have also identified the basic reproduction numberR0, which gives the expected number of new infections (in mosquitoes or humans) from one infectious individual (human or mosquito) over the duration of the infectious period, given that all other members of the population are susceptible, in terms of the model parameters. Especially, we have proved mathematically backward bifurcation may occur for R0<1 which implies that bringing the basic reproduction number below 1 is not enough to eradicate malaria.The parametersτ1 and τ2 are the incubation periods in human and mosquitoes, which are temperature dependent. With the increasing of temperature in a range, they become shorter. Since the basic reproduction number is a monotone decreasing function of both time delays, global warming will then exacerbate the transmission of malaria.The reason inducing backward bifurcation in [6] is the mosquito biting rate. While, the reasons inducing the backward bifurcation in our paper are the standard incidence rate and the disease-induced death rate which is big enough. By our results, the disease-free equilibrium is globally stable when R0<1 if the disease-induced death is omitted. However, the result about backward bifurcation in this paper has important implications for malaria control, which implies that bringing the basic reproduction number below 1 is not enough to eradicate malaria.Malaria transmission can be affected by a lot of aspects. In this paper, we are trying to model the impact of temperature increasing on malaria transmission. In some extent, our results can provide a theoretical principle for allocating and using medical health resource reasonably, which can be applied in the practice of malaria prevention and control.



---

*Source: 601265-2013-03-10.xml*