JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 427621 10.1155/2013/427621 427621 Research Article An Epidemic Model for Tick-Borne Disease with Two Delays Li Dan Ma Wanbiao Jiang Zhichao Shen Hui-Shen Department of Applied Mathematics School of Mathematics and Physics University of Science and Technology Beijing Beijing 100083 China ustb.edu.cn 2013 7 12 2013 2013 23 09 2013 18 11 2013 2013 Copyright © 2013 Dan Li et al. 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.

We have considered an epidemic model of a tick-borne infection which has nonviraemic transmission in addition to the viremic transmission. The basic reproduction number 0, which is a threshold quantity for stability of equilibria, is calculated. If 01, then the infection-free equilibrium is globally asymptotically stable, and this is the only equilibrium. On the contrary, if 0>1, then an infection equilibrium appears which is globally asymptotically stable, when one time delay is absent. By applying a permanence theorem for infinite dimensional systems, we obtain that the disease is always present when 0>1.

1. Introduction

Vector-borne diseases are infectious diseases caused by viruses, bacteria, protozoa, or rickettsia which are primarily transmitted by disease-transmitting biological agents, called vectors, who carry the disease without getting it themselves. Vector-borne disease, in particular, tick-borne disease, such as Lyme disease, transmitted between humans by blood-feeding ticks, have been big concerns for the public health in the world.

The dynamics of the tick vector are based on the biology of the tick Ixodes ricinus that transmits the virus that causes Louping ill. Louping ill is a tick-borne, zoonotic, viral disease that is most important in sheep and red grouse. This tick has a life cycle that develops from the egg through two immature stages (larva and nymph) to the adult stage. Each immature stage requires a blood meal from a suitable vertebrate host. The adult female requires a meal before producing eggs once and then dying. Adult females can only obtain a feed from large mammals, that is, hares and not grouse, whilst the immature stages will also feed on smaller warm-blooded vertebrates such as grouse chicks. Once a questing tick locates a host it generally feeds for several days. The tick life cycle usually takes 3 to 4 years but can be longer depending on host availability and climatic conditions; ticks need high relative humidity to survive .

Mathematical models have made considerable contributions to our understanding of tick-borne infections  or tick population dynamics . Transmission from host to tick is generally believed to occur when ticks bite and feed on the blood of a viraemic host. However, in recent years it has become clear that pathogen transmission can occur through a number of other routes. For example, transmission between infected and uninfected ticks that cofeed on a host can occur in the absence of systemic infection . Another route of infection for some tick-borne pathogens is through nonviraemic transmission. This occurs via the skin aided by tick saliva when a susceptible tick feeds close to an infected tick on a host that is not displaying high levels of viraemia in the blood . Grouse chicks are known to eat ticks as part of their diet of invertebrates during the first three weeks of life. Gilbert et al.  suggest that the ingestion of infected ticks may be the main route of infection in wild grouse chicks; the authors consider the ingestion model; in this case, the virus persists over a greater range of grouse densities when ingestion is included. This is because it is now easier for the grouse to become infected with the inclusion of an extra route of infection via the oral route. In this paper, we consider the nonseasonal model; that is, the behavior of grouse and ticks was assumed constant over the course of a year. The non seasonal model can be analyzed analytically unlike the seasonal models. To investigate whether the incorporation of the seasonal biology of grouse alters the model predictions dramatically from the non seasonal model. In , they present a more realistic model which describes the chicks and adult grouse separately with chicks hatching only at a given time of the year. In this paper we build on the model by Norman et al. , which have proposed an epidemic dynamics of Louping ill virus in red grouse (the viremic host) and hares (the nonviremic host) in the following form: (1)G˙s(t)=(ag-sgGs(t))Gs(t)-β1Ti(t)Gs(t)-bgGs(t),G˙i(t)=β1Ti(t)Gs(t)-ΓGi(t),G˙z(t)=γGi(t)-bgGz(t),T˙s(t)=(aT-sTTs(t))Ts(t)-β2Ts(t)Gi(t)-θHTs(t)Ti(t)-(β3H+bT)Ts(t),T˙i(t)=β2Ts(t)Gi(t)+θHTs(t)Ti(t)-(β3H+bT)Ti(t). The model considers two hosts, a viremic host that does not support adult ticks and a non-viremic host that supports all stages. The total viremic host population size at time t, given by G(t), is partitioned into subclasses of individuals who are susceptible, infectious, and recovered, with sizes denoted by Gs(t), Gi(t), and Gz(t), respectively. Hares are assumed to be at a constant density H. Ts(t) and Ti(t) represent population density of susceptible ticks and infected ticks at any time t, with the total tick population, T(t)=Ts(t)+Ti(t).

The rate of non-viremic transmission is assumed to be proportional to the number of hares present, the number of infected ticks present and the number of susceptible ticks present, and therefore takes the form θHTiTs, where θ is a measure of the probability of non-viremic transmission occurring. Furthermore, the viremic host population dies at a natural death rate bg. The vertical transmission in the viremic host does not occur so that all newly recruited individuals are susceptible. The rate at which a grouse is bitten and infected by a tick β1 can only involve the nymphal stage as adults do not generally bite grouse and the larvae are not infected. The rate at which nymphs bite grouse was estimated in  from tick counts of nymphs on deer and grouse as well as from model predictions; γ is the rate at which infectious grouse become immune, and Γ=α+bg+γ is the rate at which infectious grouse are lost from the system and is the sum of the death rate, due to the disease, the natural death rate, and the rate at which individuals become immune. For the ticks, bT is the per capita natural death rate. β2 is the probability of an average tick biting a grouse and becoming infected per unit time, β3 is the probability of an adult tick being female, biting a hare, and then going on to reproduce per unit time before dying and being lost from the system.

Norman et al.  consider the local stability of six biological relevant equilibria; this is done by determining the eigenvalues of the Jacobian evaluated at each equilibrium. The author chose the logistic growth function to describe the growth of grouse and ticks, obviously, it is more reasonable, but it increases the complexity of the algebraic analysis. To obtain global properties for model (1), we suppose that the host and tick population are recruited at rate that  λ  and  s. Note that the variables Gz do not appear in the other equations of the above system. This allows us to consider the following fours dimensional system: (2)G˙s(t)=λ-β1Ti(t)Gs(t)-bgGs(t),G˙i(t)=β1Ti(t)Gs(t)-ΓGi(t),T˙s(t)=s-β2Ts(t)Gi(t)-θHTs(t)Ti(t)-(β3H+bT)Ts(t),T˙i(t)=β2Ts(t)Gi(t)+θHTs(t)Ti(t)-(β3H+bT)Ti(t). The incubation period for Louping ill is 6 to 18 days in sheep. Parenterally inoculated red grouse develop clinical signs in 2 to 8 days; the dynamic of tick feeding the incubation is effectively 12 months. We introduce two types of the delays that correspond to (i) the latent period in an infected host and (ii) the latency in a vector. We know that the number of infected hosts at time t may depend on the number of susceptible hosts which is infected by the infected ticks at time t-τ, for a time delay τ>0, and e-bgτ is the survival rate. The number of infected ticks at time t may depend on the number of susceptible ticks which is feeding on the infected hosts at time t-ω, for a time delay ω>0, and e-bTω is the survival rate. Thus, the following model is proposed: (3)G˙s(t)=λ-e-bgτβ1Ti(t-τ)Gs(t)-bgGs(t),G˙i(t)=e-bgτβ1Ti(t-τ)Gs(t)-ΓGi(t),T˙s(t)=s-e-bTωβ2Ts(t)Gi(t-ω)-θHTs(t)Ti(t)-(β3H+bT)Ts(t),T˙i(t)=e-bTωβ2Ts(t)Gi(t-ω)+θHTs(t)Ti(t)-(β3H+bT)Ti(t). Although the models incorporation of time delay have more real biological meaning, they are difficult to deal with it mathematically. However, the Lyapunov functional and the LaSalle-type theorem in  can provide a direct and effective method to establish global dynamical properties for the system of those nonlinear functional differential equations. Recently, Korobeinikov , and Huang et al.  have investigated the global dynamical properties of epidemiological models with and without delay by construc suitable Lyapunov functional methods. Inspired by the work in McCluskey , we construct a global Lyapunov functional and show that the disease-free equilibrium of system (3) is globally asymptotically stable if 01, and if 0>1 the infection equilibrium is globally asymptotically stable when one delay is absent. Furthermore, by applying a permanence theorem by Smith and Zhao  for infinite dimensional systems (also see ), we obtain that the disease is always present when 0>1. For its applications, we refer the reader to [10, 11, 2123].

The paper is organized as follows. In the next section, we analyze some basic mathematical properties (positivity and boundedness of solutions). In Section 3, we give a complete global stability analysis of the disease-free equilibrium. In Section 4, we are concerned with the global stability of the infection equilibrium with one delay. In Section 5, we show that the disease is endemic in the sense of permanence whenever 0>1. A discussion of the mathematical results and of their biological implications is presented in Section 6.

2. Nonnegativity and Boundedness of Solutions

Let C=C([-σ,0];R4) be the Banach space of continuous functions from [-σ,0] to R4 equipped with the sup-norm, where σ=max{τ,ω}. The initial condition of (3) is given as (4)Gs(θ)=φ1(θ),Gi(θ)=φ2(θ),Ts(θ)=φ3(θ),Ti(θ)=φ4(θ)θ[-σ,0], where φ=(φ1,φ2,φ3,φ4)C such that φj(θ)0  (σθ0,j=1,2,3,4).

The following result establishes the feasible region of the model and shows that the model is wellposed.

Theorem 1.

Under initial conditions in (4), all solutions of system (3) are nonnegative on [0,+) and ultimately bounded.

Proof.

First, we prove that Gs(t) is positive for t0. Assuming the contrary and letting t1>0 be the first time such that Gs(t1)=0, by the first equation of system (3), we have Gs(t1)=λ>0, and hence Gs(t)<0 for t(t1-ϵ,t1), where ϵ is an arbitrarily small positive constant. This leads to a contradiction. It follows that Gs(t) is always positive. Using the same method, we have that Ts(t) is always positive. Further, from the second and the forth equations in (3), we have, respectively, (5)Gi(t)=Gi(0)e-Γt+e-bgτβ10tTi(ξ-τ)Gs(ξ)e-Γ(t-ξ)dξ,Ti(t)=Ti(0)e-0t(θHTs(γ)-(β3H+bT))dγ+e-bTωβ2×0tTs(ξ)Gi(ξ-ω)e-ξt(θHTs(γ)-(β3H+bT))dγdξ. Using (4) we know that the Gi(t) and Ti(t) of (3) are positive for all t0.

Next we show that positive solutions of (3) are ultimately bounded for t0.

Define F(t)=Gs(t)+Gi(t)+Ts(t)+Ti(t) and d=min{bg,β3H+bT}. By non-negativity of the solution, it follows that (6)F(t)=λ+s-(bgGs+ΓGi+(β3H+bT)×(Ts+Ti)F)λ+s-dF(t). This implies that F(t) is ultimately bounded, and so are Gs(t), Gi(t), Ts(t), and Ti(t).

This completes the proof of Theorem 1.

The basic reproduction number  for system (3) is (7)0=e-bgτ-bTωλsβ1β2Γbg(β3H+bT)2+sθH(β3H+bT)2.0 represents the average number of secondary infections that single infections host can generate in a totally susceptible population of hosts and ticks. The basic reproduction number consists of two terms representing two routes of infection. The tick would live for 1/(β3H+bT) units of time and produces (e-bgτλβ1)/bg(β3H+bT) infected grouse through feeding and θsH/(β3H+bT)2 infected ticks through cofeeding. Each infected grouse lives for 1/Γ units of time and produces e-bTωβ2s/Γ(β3H+bT) infected ticks.

System (3) has always an infection-free equilibrium E0=(Gs0,0,Ts0,0), that is, grouse and ticks both at their carrying capacity with no disease present, where Gs0=λ/bg, Ts0=s/(β3H+bT), and an infection equilibrium E*=(Gs*,Gi*,Ts*,Ti*), when 0>1, that is, ticks, grouse, and disease coexisting together, where Ti* is given by the following cubic: (8)f(Ti*)=ATi*2+BTi*+C, where (9)A=e-bgτΓθHβ1(β3H+bT),B=-sθHΓβ1e-bgτ+e-bgτ-bTωλβ1β2(β3H+bT)+Γ(β3H+bT)(θHbg+e-bgτβ1(β3H+bT)),C=-e-bgτ-bTωλsβ1β2-sΓθHbg+Γbg(β3H+bT)2. Hence, we can write all of the densities in terms of the infectious tick density, Ti*, (10)Gs*=λe-bgτβ1Ti*+bg,Gi*=e-bgτλβ1Ti*Γ(e-bgτβ1Ti*+bg),Ts*=sβ3H+bT-Ti*.

3. Global Stability of the Infection-Free Equilibrium of Model <bold>(<xref ref-type="disp-formula" rid="EEq1.3">3</xref>)</bold>

In this section, we would consider the global stability of the infection-free equilibrium. Next, we show that 0 determines the stability of the infection-free equilibrium.

Theorem 2.

If 01, then the infection-free equilibrium E0 is globally asymptotically stable for any time delay τ0 and ω0.

Proof.

Consider a Lyapunov functional: (11)V(t)=V1(t)+V2(t)+V3(t), with (12)V1(t)=e-bTωβ2Ts0Γ×(Gs(t)-Gs0-Gs0lnGs(t)Gs0)+e-bTωβ2Ts0ΓGi(t)+(Ts(t)-Ts0-Ts0lnTs(t)ts0)+Ti(t),V2(t)=e-bTωβ2Ts0t-ωtGi(ξ)dξ,V3(t)=e-(bgτ+bTω)β1β2Gs0Ts0Γt-τtTi(ξ)dξ. Calculate the time derivative of V: (13)dVdt=e-bTωβ2Ts0Γ(1-Gs0Gs(t))Gs(t)+e-bTωβ2Ts0ΓGi(t)+(1-Ts0Ts(t))Ts(t)+Ti(t)+e-bTωβ2Ts0(Gi(t)-Gi(t-ω))+e-(bgτ+bTω)β1β2Gs0Ts0Γ(Ti(t)-Ti(t-τ)). Since λ=bgGs0, s=(β3H+bT)Ts0, it follows that (14)dVdt=β2bgGs0Ts0Γ(2-Gs0Gs(t)-Gs(t)Gs0)+(β3H+bT)Ts0(2-Ts0Ts(t)-Ts(t)Ts0)+(β3H+bT)(0-1)Ti(t). Obviously 01 ensures that dV/dt0. It is clear that V0 and V=0 if and only if Gs(t)=Gs0, Gi(t)=0, Ts(t)=Ts0, and Ti(t)=0. Hence, it follows from the Lyapunov-LaSalle invariance principle that the infection-free equilibrium E0 is globally asymptotically stable when 01.

This completes the proof of Theorem 2.

4. Dynamical Analysis for the Infection Equilibrium

System (3) has an infection equilibrium E* when 0>1. We consider two situations: (a) τ=0, ω0 and (b) τ0, ω=0.

4.1. Global Stability of the Infection Equilibrium When <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M118"><mml:mi>τ</mml:mi><mml:mo mathvariant="bold">=</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula> and ss<inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M119"><mml:mi>ω</mml:mi><mml:mo mathvariant="bold">≥</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula>

In this subsection, the situations τ=0 and ω0 are considered.

Theorem 3.

If 0>1, then the infection equilibrium E* is globally asymptotically stable for τ=0 and ω0.

Proof.

Define a Lyapunov functional: (15)U1(t)=e-bTωβ2Ts*Gi*β1Ti*Gs*(Gs(t)-Gs*-Gs*lnGs(t)Gs*)+e-bTωβ2Ts*Gi*β1Ti*Gs*(Gi(t)-Gi*-Gi*lnGi(t)Gi*)+(Ts(t)-Ts*-Ts*lnTs(t)Ts*)+(Ti(t)-Ti*-Ti*lnTi(t)Ti*)+e-bTωβ2Ts*Gi*U+(t), where (16)U+(t)=0ω(Gi(t-θ)Gi*-1-lnGi(t-θ)Gi*)dθ. The derivatives of U+(t) are given by (17)dU+(t)dt=ddt0ω(Gi(t-θ)Gi*-1-lnGi(t-θ)Gi*)dθ=0ωddt(Gi(t-θ)Gi*-1-lnGi(t-θ)Gi*)dθ=-0ωddθ(Gi(t-θ)Gi*-1-lnGi(t-θ)Gi*)dθ=Gi(t)Gi*-Gi(t-ω)Gi*+lnGi(t-ω)Gi(t). Hence, we obtain (18)dU1dt=e-bTωβ2Ts*Gi*β1Ti*Gs*(1-Gs*Gs(t))Gs(t)+e-bTωβ2Ts*Gi*β1Ti*S*(1-Gi*Gi(t))Gi(t)+(1-Ts*Ts(t))Ts(t)+(1-Ti*Ti(t))Ti(t)+e-bTωβ2Ts*Gi*dU+(t)dt=e-bTωβ2Ts*Gi*β1Ti*Gs*(Gs*Gs(t)(λ-β1Ti(t)Gs(t)-bgGs(t))λ-β1Ti(t)Gs(t)-bgGs(t)+Gs*Gs(t)(λ-β1Ti(t)Gs(t)-bgGs(t)))+e-bTωβ2Ts*Gi*β1Ti*Gs*(+Gi*Gi(t)(β1Ti(t)Gs(t)-ΓGi(t))β1Ti(t)Gs(t)-ΓGi(t)+Gi*Gi(t)(β1Ti(t)Gs(t)-ΓGi(t)))+(-θHTs(t)Ti(t)-(β3H+bT)Ts(t)s-e-bTωβ2Ts(t)Gi(t-ω)-θHTs(t)Ti(t)  -(β3H+bT)Ts(t)+Ts*Ts(t)(s-e-bTωβ2Ts(t)Gi(t-ω)-θHTs(t)Ti(t)-(β3H+bT)Ts(t)e-bTω))+(+  θHTs(t)Ti(t)-(β3H+bT)Ti(t)e-bTω)e-bTωβ2Ts(t)Gi(t-ω)+θHTs(t)Ti(t)-(β3H+bT)Ti(t)+Ti*Ti(t)(e-bTωβ2Ts(t)Gi(t-ω)+  θHTs(t)Ti(t)-(β3H+bT)Ti(t)e-bTω))+e-bTωβ2Ts*Gi*(Gi(t)Gi*-Gi(t-ω)Gi*+lnGi(t-ω)Gi(t)Gi(t)Gi*). Since λ=β1Ti*Gs*+bgGs*, s=e-bTωβ2Ts*Gi*-θHTs*Ti*-(β3H+bT)Ts*, and (19)ΓGi*=β1Ti*Gs*,(β3H+bT)Ti*=e-bTωβ2Ts*Gi*+θHTs*Ti*, it follows that (20)dU1dt=e-bTωβ2bgGs*Ts*β1Ti*×(2-Gs*Gs(t)-Gs(t)Gs*)+θHTi*Ts*(2-Ts*Ts(t)-Ts(t)Ts*)+(β3H+bT)Ts*(2-Ts*Ts(t)-Ts(t)Ts*)+e-bTωβ2Ts*Gi*×(4-Gs*Gs(t)-Gi*Gs(t)Ti(t)Gs*Ti*Gi(t)-Ts*Ts(t)-Ti*Ts(t)Gi(t-ω)Ts*Gi*Ti(t)). Using equality (21)lnGi(t-ω)Gi(t)=lnGs*Gs(t)+lnTs*Ts(t)+lnTi*Ts(t)Gi(t-ω)Ts*Gi*Ti(t)+lnGi*Gs(t)Ti(t)Gs*Ti*Gi(t), we obtain (22)dU1dt=-bgGs*1Gs(t)(Gs(t)-Gs*)2-θHTi*Ts*1Ts(t)(Ts(t)-Ts*)2-(β3H+bT)Ts*1Ts(t)(Ts(t)-Ts*)2+e-bTωβ2Ts*Gi*(Ω1+Ω2+Ω3+Ω4), where (23)Ω1=1-Gs*Gs(t)+lnGs*Gs(t),Ω2=1-Gi*Gs(t)Ti(t)Gs*Ti*Gi(t)+lnGi*Gs(t)Ti(t)Gs*Ti*Gi(t),Ω3=1-Ts*Ts(t)+lnTs*Ts(t),Ω4=1-Ti*Ts(t)Gi(t-ω)Ts*Gi*Ti(t)+lnTi*Ts(t)Gi(t-ω)Ts*Gi*Ti(t). Since the function (24)g(t)=1-f(t)+lnf(t) is always nonpositive for any function f(t)>0, and g(t)=0 if and only if f(t)=1. Therefore, we have Ω1,Ω2,Ω3,Ω40.

It is easy to see that dU1/dt0. It is clear that dU1/dt=0 if and only if Gs(t)=Gs*, Gi(t)=Gi*, Ts(t)=Ts*, Ti(t)=Ti*. Hence, it is also known from stability theorem  that the infection equilibrium E* is globally asymptotically stable for any time delay ω0 under the condition 0>1.

This completes the proof of Theorem 3.

4.2. Global Stability of the Infection Equilibrium When <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M152"><mml:mi>τ</mml:mi><mml:mo mathvariant="bold"> </mml:mo><mml:mo mathvariant="bold">≥</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula> and <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M153"><mml:mi>ω</mml:mi><mml:mo> </mml:mo><mml:mo mathvariant="bold">=</mml:mo><mml:mo> </mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula>

In this subsection, the situations τ0 and ω=0 are considered.

Theorem 4.

If 0>1, then the infection equilibrium E* is globally asymptotically stable for τ0 and ω=0.

Proof.

Define a Lyapunov functional: (25)U2(t)=ebgτβ2Ts*Gi*β1Ti*Gs*(Gs(t)-Gs*-Gs*lnGs(t)Gs*)+ebgτβ2Ts*Gi*β1Ti*Gs*(Gi(t)-Gi*-Gi*lnGi(t)Gi*)+(Ts(t)-Ts*-Ts*lnTs(t)Ts*)+(Ti(t)-Ti*-Ti*lnTi(t)Ti*)+β2Ts*Gi*U-(t), where (26)U-(t)=0τ(Ti(t-θ)Ti*-1-lnTi(t-θ)Ti*)dθ. The derivatives of U-(t) are given by (27)dU-(t)dt=ddt0τ(Ti(t-θ)Ti*-1-lnTi(t-θ)Ti*)dθ=0τddt(Ti(t-θ)Ti*-1-lnTi(t-θ)Ti*)dθ=-0τddθ(Ti(t-θ)Ti*-1-lnTi(t-θ)Ti*)dθ=Ti(t)Ti*-Ti(t-τ)Ti*+lnTi(t-τ)Ti(t). Hence, we obtain (28)dU2dt=ebgτβ2Ts*Gi*β1Ti*Gs*(1-Gs*Gs(t))Gs(t)+ebgτβ2Ts*Gi*β1Ti*S*(1-Gi*Gi(t))Gi(t)+(1-Ts*Ts(t))Ts(t)+(1-Ti*Ti(t))Ti(t)+β2Ts*Gi*dU-(t)dt=bgβ2Gs*Ts*Gi*β1Ti*(2-Gs*Gs(t)-Gs(t)Gs*)+θHTi*Ts*(2-Ts*Ts(t)-Ts(t)Ts*)+(β3H+bT)Ts*(2-Ts*Ts(t)-Ts(t)Ts*)+β2Ts*Gi*lnTi(t-τ)Ti(t)+β2Ts*Gi*(4-Gs*Gs(t)-Gi*Gs(t)Ti(t-τ)Gs*Ti*Gi(t)-Ts*Ts(t)-Ti*Ts(t)Gi(t)Ts*Gi*Ti(t)). Using equality (29)lnTi(t-τ)Ti(t)=lnGs*Gs(t)+lnTs*Ts(t)+lnTi*Ts(t)Gi(t)Ts*Gi*Ti(t)+lnGi*Gs(t)Ti(t-τ)Gs*Ti*Gi(t), we obtain (30)dU2dt=-bgβ2Gs*Ts*Gi*β1Ti*1Gs(t)(Gs(t)-Gs*)2-θHTi*Ts*1Ts(t)(Ts(t)-Ts*)2-(β3H+bT)Ts*1Ts(t)(Ts(t)-Ts*)2+β2Ts*Gi*(Ω1+Ω5+Ω3+Ω6), where (31)Ω5=1-Gi*Gs(t)Ti(t-τ)Gs*Ti*Gi(t)+lnGi*Gs(t)Ti(t-τ)Gs*Ti*Gi(t),Ω6=1-Ti*Ts(t)Gi(t)Ts*Gi*Ti(t)+lnTi*Ts(t)Gi(t)Ts*Gi*Ti(t). Since the function (32)g(t)=1-f(t)+lnf(t) is always non-positive for any function f(t)>0, and g(t)=0 if and only if f(t)=1. Therefore, we have Ω1,Ω5,Ω3,Ω60.

It is easy to see that dU2/dt0. It is clear that dU2/dt=0 if and only if Gs(t)=Gs*, Gi(t)=Gi*, Ts(t)=Ts*, Ti(t)=Ti*. Hence, it is also known from stability theorem  that the infection equilibrium E* is globally asymptotically stable for any time delay τ0 under the condition 0>1.

This completes the proof of Theorem 4.

5. Permanence

System (3) is said to be uniformly persistent if there exists a constant η>0 such that any solution (Gs(t),Gi(t),Ts(t),Ti(t)) of (3) satisfies (33)liminftGs(t)η,liminftGi(t)η,liminftTs(t)η,liminftTi(t)η.

Now we give a result on the uniform persistence of system (3). To proceed, we introduce the following notation and terminology. Denote by P(t), t0 the family of solution operators corresponding to (3). The ω-limit set ω(x) of x consists of yX such that there exists a sequence tn as n with P(tn)xy as n.

Theorem 5.

System (3) is uniformly persistent if it satisfies 0>1.

Proof .

Let (34)X0={ϕX:ϕ2(0)>0,ϕ4(0)>0},X=XX0={ϕX:ϕ2(0)=0  or  ϕ4(0)=0}, which is relatively closed in X.

Now we prove that X0 is positively invariant for P(t). By the second and forth equations of (3), we have (35)dGi(t)dt-ΓGi(t),dTi(t)dt-(β3H+bT)Ti(t),kkkkkkkkkkkkkkkkit0. Since Gi(0,ϕ)=ϕ2(0)>0 and Ti(0,ϕ)=ϕ4(0)>0, it follows from (35) that (36)Gi(t,ϕ)ϕ2(0)e-Γt,Ti(t,ϕ)ϕ4(0)e-(β3H+bT)t,kkkkkkkkkkkkkkkkkt0. Thus, X0 is positively invariant for P(t).

We set (37)M={ϕX:Y(t)ϕ  satisfies (3),Y(t)ϕX,t0}. We claim that (38)M={(Gs,0,Ts,0)}. Assuming Y(t)M, for all t0, it suffices to show that Gi(t)=Ti(t)=0, for all t0. If it is not true, then there exists t0>0 such that either (a) Gi(t0)=0, Ti(t0)>0 or (b) Gi(t0)>0, Ti(t0)=0.

For case (a), from the forth equation of (3), we have (39)dTi(t)dt-(β3H+bT)Ti(t). Since Ti(t0)>0, we have (40)Ti(t)Ti(t0)e-(β3H+bT)(t-t0)>0,tt0. By the invariant of M, we have Gi(t)0, for all t0. Then from the second equation of (3), we have Ti(t-τ)Gs(t)0, for all  t0.

From the first equation of (3), we obtain G˙s(t)=λ-bgGs(t), thus limt+Gs(t)=λ/bg. By the invariant of M, we have Gs(t)=λ/bg.

Then from the second equation of (3), we obtain Ti(t-τ)=0; by the invariant of M, we have Ti(t)0, for all t0. From the first and third equations of (3), we obtain T˙s(t)=s-(β3H+bT)Ts(t); thus limt+Ts(t)=s/(β3H+bT). By the invariant of M, we have Ts(t)=s/(β3H+bT  ), t0, which contradicts the assumption that (Gs(t),Gi(t),Ts(t),Ti(t))M, for all t0. Similarly, we can obtain a contradiction for case (b). This proves the claim (38).

Let A=xAbω(x), where Ab is the global attractor of P(t) restricted to X. We show that A=E0. In fact, from AM and the first and third equations of (3), we have limtGs(t)=Gs0 and limtTs(t)=Ts0. Thus, E0 is the isolated invariant set in X.

Next, we show that Ws(E0)X0=. If this is not true, then there exists a solution (Gst,Git,Tst,Tit)X0 such that (41)limtGs(t)=Gs0,limtGi(t)=0,limtTs(t)=Ts0,limtTi(t)=0. For any sufficiently small constant ɛ>0, there exists a positive constant T0=T0(ɛ) such that (42)Gs(t)>Gs0-ɛ>0,tT0,Ts(t)>Ts0-ɛ>0,tT0. For the constant ɛ given above, it follows from the second and forth equations of (3) that (43)G˙i(t)e-bgτβ1Ti(t-τ)×(Gs0-ɛ)-ΓGi(t),T˙i(t)e-bTωβ2(Ts0-ɛ)×Gi(t-ω)+θH(Ts0-ɛ)Ti(t)-(β3H+bT)Ti(t). If Gi(t),Ti(t)0, as t, then by a standard comparison argument and the nonnegativity, the solution (x(t),y(t)) of the following monotone system (44)x˙(t)=e-bgτβ1y(t-τ)(Gs0-ɛ)-Γx(t),y˙(t)=e-bTωβ2(Ts0-ɛ)x(t-ω)+θH(Ts0-ɛ)y(t)-(β3H+bT)y(t) with the initial condition x(t)=Gi(t), y(t)=Ti(t), for all t[t0,t0+T0] converges to (0,0) as well. Thus, limtW(t)=0, where W(t)>0 is defined by (45)W(t)=e-bTωβ2(Ts0-ɛ)Γx(t)+y(t)+e-bTωβ2(Ts0-ɛ)t-ωtx(ξ)dξ+e-bgτ-bTωβ1β2(Gs0-ɛ)(Ts0-ɛ)Γt-τty(ξ)dξ. Differentiating W(t) with respect to time t gives (46)dW(t)dt|(44)=[θH(Ts0-ɛ)+e-bgτ-bTωβ1β2(Gs0-ɛ)(Ts0-ɛ)Γ-(β3H+bT)θH(Ts0-ɛ)]y(t). Because 0>1, we have θH(Ts0-ɛ)+e-bgτ-bTωβ1β2(Gs0-ɛ)(Ts0-ɛ)/Γ-(β3H+bT)>0 for a sufficiently small ɛ. Therefore, W(t) goes to either infinity or a positive number t, which leads to a contradiction with limtW(t)=0. Thus we have Ws(E0)X0=.

Define p:XR+ by (47)p(ϕ)=min{ϕ2(0),ϕ4(0)},ϕX. It is clear that X0=p-1(0,) and X=p-1(0). Thus by [20, Theorem 3], we have (48)liminft(Gi(t),Ti(t))(η0,η0) for some constant η0>0. Let η=min{η0,ɛ}, where ɛ is the constant such that  liminftGs(t)ɛ>0 and liminftTs(t)ɛ>0. We showed that liminftGs(t)η, liminftGi(t)η, liminftTs(t)η, and liminftTi(t)η.

This completes the proof of Theorem 5.

6. Numerical Simulations

In this section, by using suitable numerical methods, we simulate the systems to illustrate the analytical results obtained in Sections 35. Moreover, we would verify the performance as shown in Section 5 and find that the system (3) is globally asymptotically stable when two delays are present simultaneously.

Table 1 shows the literature reported values for parameters which we used for the simulations. Besides the parameter value shown in Table 1, we can find η by using the equation η=α+γ+bg.

Parameter values used in the simulations shown in figures.

Parameter Value per month Source
λ varies Assumed
b g 0.087 20
β 1 0.01 2
β 2 0.02525 2
α 2.31 20
γ 0.5775 2
s varies Assumed
θ 0.00001 2
H 10 2
b T 0.0277 Unpublished date
β 3 0.0113 2

For system (3), we choose λ=3, s=2, τ=10, and ω=15, and other values are shown in Table 1. By computation, we obtain R0=0.0818<1, and the infection equilibrium does not exist. By the analysis of Section 2, we know that the infection-free equilibrium E0=(34.4828,0,14.2146,0) is globally asymptotically stable (see Figure 1). Furthermore, we choose λ=8, s=6, τ=0, and ω=15, and other values are shown in Table 1. In Figure 2, we know that under the condition indicated in the Theorem 3 with R0=1.5645>1 the infection equilibrium E*=(63.7684,1.324,38.421,4.15) will be globally asymptotically stable. In Figure 3, we know that under the condition indicated in the Theorem 4 with R0=1.8425>1, the infection equilibrium E*=(75.346,1.237,39.96,3.856) will be globally asymptotically stable with τ=3 and ω=0. In Figure 4, we know that under the condition indicated in Theorem 5 with R0=1.1638>1, the infection equilibrium E*=(68.72,1.0472,39.456,4.385) will be globally asymptotically stable with τ=5 and ω=10.

The infection-free equilibrium E0 of system (3) is globally asymptotically stable when τ=10, ω=15, and 01. The initial functions are (60,10,40, and 10).

The infection equilibrium E* of system (3) is asymptotically stable when τ=0, ω=15, and 0>1. The initial functions are (100,10,100, and 10).

The infection equilibrium E* of system (3) is asymptotically stable when  τ=3, ω=0, and 0>1. The initial functions are (100,10,100, and 10).

The infection equilibrium E* of system (3) is asymptotically stable when τ=5, ω=10, and 0>1. The initial functions are (100,10,100, and 10).

7. Conclusion

In the present paper, we investigate an epidemic model for tick-borne disease with non-viraemic transmission and two time delays. By constructing Lyapunov functionals, we prove that the infection-free equilibrium E0 of system (3) is globally asymptotically stable when 01. If 0>1 then the infection equilibrium E* is globally asymptotically stable when one delay is absent. By applying the persistence theory for infinite dimensional systems, we obtain that the populations of system (3) can coexist permanently if 0>1.

The reproductive number plays a crucial role for epidemic dynamics. Considering the basic reproductive number 0=e-bgτ-bTωλsβ1β2/Γbg(β3H+bT)2+sθH/(β3H+bT)2 as a function of τ and ω, we can find that it is decreasing in τ and ω and it tends to 0 if the time delay τ or ω tends to 1. Therefore, the basic reproductive number can be reduced by increasing the intracellular delay. This may provide an insight for developing vaccination and regular treatment with acaricides that try to increase the intracellular delay.

It is interesting to consider the global stability of the infection equilibrium E* with two time delays; from the numerical simulations, it is shown that the infection equilibrium E* is globally asymptotically stable when 0>1; we leave this as a future work.

The nonseasonal model can be analysed analytically unlike the seasonal models. The non seasonal model is also simpler and quicker to implement in Mathematica. However, the effect of the seasonal behavior of ticks will be explored before deciding whether any form of seasonality will be incorporated in future work.

Acknowledgment

The research was partially supported by the Funds of USTB and NNSF of China (11071013).

Norman R. Ross D. Laurenson M. K. Hudson P. J. The role of non-viraemic transmission on the persistence and dynamics of a tick borne virus—Louping ill in red grouse (Lagopus lagopus scoticus) and mountain hares (Lepus timidus) Journal of Mathematical Biology 2004 48 2 119 134 10.1007/s00285-002-0183-5 MR2041497 ZBL1058.92040 Norman R. Bowers R. G. Begon M. Hudson P. J. Persistence of tick-borne virus in the presence of multiple host species: tick reservoirs and parasite mediated competition Journal of Theoretical Biology 1999 200 1 111 118 2-s2.0-0033533378 10.1006/jtbi.1999.0982 Caraco T. Gardner G. Maniatty W. Deelman E. Szymanski B. K. Lyme disease: self-regulation and pathogen invasion Journal of Theoretical Biology 1998 193 4 561 575 2-s2.0-0032555474 10.1006/jtbi.1998.0722 Rosà R. Pugliese A. Norman R. Hudson P. J. Thresholds for disease persistence in models for tick-borne infections including non-viraemic transmission, extended feeding and tick aggregation Journal of Theoretical Biology 2003 224 3 359 376 10.1016/S0022-5193(03)00173-5 MR2067243 Sandberg S. Awerbuch T. E. Spielman A. A comprehensive multiple matrix model representing the life cycle of the tick that transmits the agent of Lyme disease Journal of Theoretical Biology 1992 157 2 203 220 2-s2.0-0026615693 Awerbuch-Friedlander T. Levins R. Predescu M. The role of seasonality in the dynamics of deer tick populations Bulletin of Mathematical Biology 2005 67 3 467 486 10.1016/j.bulm.2004.08.003 MR2215551 Hudson P. J. Norman R. Laurenson M. K. Newborn D. Gaunt M. Jones L. Reid H. Gould E. Bowers R. Dobson A. Persistence and transmission of tick-borne viruses: ixodes ricinus and louping-ill virus in red grouse populations Parasitology 1995 111 S49 S58 2-s2.0-0029550954 Hudson P. J. Grouse in space and time: the population biology of a managed gamebird The report of the Game conservancys Scottish Grouse Research project and North of England Grouse Research project. Game Conservancy Trust Fordingbridge, 1992 Zhang H. Chen L. A model for two species with stage structure and feedback control International Journal of Biomathematics 2008 1 3 267 286 10.1142/S1793524508000230 MR2449829 ZBL1155.92359 Song X. Guo H. Global stability of a stage-structured predator-prey system International Journal of Biomathematics 2008 1 3 313 326 10.1142/S1793524508000266 MR2449832 ZBL1173.34043 Mwambi H. G. Baumgartner J. Hadeler K. P. Ticks and tick-borne diseases: a vector-host interaction model for the brown ear tick (Rhipicephalus appendiculatus) Statistical Methods in Medical Research 2000 9 3 279 301 2-s2.0-0033781112 Labuda M. Jones L. D. Williams T. Danielova V. Nuttall P. A. Efficient transmission of tick-borne encephalitis virus between cofeeding ticks Journal of Medical Entomology 1993 30 1 295 299 2-s2.0-0027344629 Gilbert L. Jones L. D. Laurenson M. K. Gould E. A. Reid H. W. Hudson P. J. Ticks need not bite their red grouse hosts to infect them with louping ill virus Proceedings of the Royal Society B 2004 271 4 S202 S205 2-s2.0-2342473275 Gilbert L. Norman R. Laurenson K. M. Reid H. W. Hudson P. J. Disease persistence and apparent competition in a three-host community: an empirical and analytical study of large-scale, wild populations Journal of Animal Ecology 2001 70 6 1053 1061 2-s2.0-0035214189 10.1046/j.0021-8790.2001.00558.x Kuang Y. Delay Differential Equations with Applications in Population Dynamics 1993 191 Boston, Mass, USA Academic Press xii+398 Mathematics in Science and Engineering MR1218880 Korobeinikov A. Lyapunov functions and global properties for SEIR and SEIS epidemic models Mathematical Medicine and Biology 2004 21 2 75 83 2-s2.0-3042703290 Korobeinikov A. Maini P. K. A Lyapunov function and global properties for SIR and SEIR epidemiological models with nonlinear incidence Mathematical Biosciences and Engineering 2004 1 1 57 60 10.3934/mbe.2004.1.57 MR2130156 ZBL1062.92061 Korobeinikov A. Global properties of infectious disease models with nonlinear incidence Bulletin of Mathematical Biology 2007 69 6 1871 1886 10.1007/s11538-007-9196-y MR2329184 ZBL05265715 Huang G. Takeuchi Y. Ma W. Wei D. Global stability for delay SIR and SEIR epidemic models with nonlinear incidence rate Bulletin of Mathematical Biology 2010 72 5 1192 1207 10.1007/s11538-009-9487-6 MR2669941 ZBL1197.92040 McCluskey C. C. Complete global stability for an SIR epidemic model with delay—distributed or discrete Nonlinear Analysis: Real World Applications 2010 11 1 55 59 10.1016/j.nonrwa.2008.10.014 MR2570523 ZBL1185.37209 Smith H. L. Zhao X.-Q. Robust persistence for semidynamical systems Nonlinear Analysis: Theory, Methods & Applications 2001 47 9 6169 6179 10.1016/S0362-546X(01)00678-2 MR1971507 Zhao X.-Q. Dynamical Systems in Population Biology 2003 New York, NY, USA Springer xiv+276 CMS Books in Mathematics/Ouvrages de Mathématiques de la SMC, 16 MR1980821 Thieme H. R. Persistence under relaxed point-dissipativity (with application to an endemic model) SIAM Journal on Mathematical Analysis 1993 24 2 407 435 10.1137/0524026 MR1205534 ZBL0774.34030 Diekmann O. Heesterbeek J. A. P. Metz J. A. J. On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations Journal of Mathematical Biology 1990 28 4 365 382 10.1007/BF00178324 MR1057044 ZBL0726.92018 Hale J. K. Verduyn Lunel S. M. Introduction to Functional-Differential Equations 1993 99 New York, NY, USA Springer x+447 Applied Mathematical Sciences MR1243878