Journal of Applied Mathematics

Volume 2013 (2013), Article ID 427621, 11 pages

http://dx.doi.org/10.1155/2013/427621

## An Epidemic Model for Tick-Borne Disease with Two Delays

Department of Applied Mathematics, School of Mathematics and Physics, University of Science and Technology Beijing, Beijing 100083, China

Received 23 September 2013; Accepted 18 November 2013

Academic Editor: Hui-Shen Shen

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.

#### Abstract

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 , which is a threshold quantity for stability of equilibria, is calculated. If , then the infection-free equilibrium is globally asymptotically stable, and this is the only equilibrium. On the contrary, if , 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 .

#### 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 [1].

Mathematical models have made considerable contributions to our understanding of tick-borne infections [1–8] or tick population dynamics [9–11]. 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 [1]. 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 [12]. Grouse chicks are known to eat ticks as part of their diet of invertebrates during the first three weeks of life. Gilbert et al. [13] 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 [13], 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. [1], 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: 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 , given by , is partitioned into subclasses of individuals who are susceptible, infectious, and recovered, with sizes denoted by , , and , respectively. Hares are assumed to be at a constant density . and represent population density of susceptible ticks and infected ticks at any time , with the total tick population, .

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 , where is a measure of the probability of non-viremic transmission occurring. Furthermore, the viremic host population dies at a natural death rate . 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 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 [14] 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 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, is the per capita natural death rate. is the probability of an average tick biting a grouse and becoming infected per unit time, 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. [1] 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 thatand. Note that the variables do not appear in the other equations of the above system. This allows us to consider the following fours dimensional system: 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 may depend on the number of susceptible hosts which is infected by the infected ticks at time , for a time delay , and is the survival rate. The number of infected ticks at time may depend on the number of susceptible ticks which is feeding on the infected hosts at time , for a time delay , and is the survival rate. Thus, the following model is proposed: 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 [15] can provide a direct and effective method to establish global dynamical properties for the system of those nonlinear functional differential equations. Recently, Korobeinikov [16–18], and Huang et al. [19] 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 [20], we construct a global Lyapunov functional and show that the disease-free equilibrium of system (3) is globally asymptotically stable if , and if the infection equilibrium is globally asymptotically stable when one delay is absent. Furthermore, by applying a permanence theorem by Smith and Zhao [21] for infinite dimensional systems (also see [22]), we obtain that the disease is always present when . For its applications, we refer the reader to [10, 11, 21–23].

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 . A discussion of the mathematical results and of their biological implications is presented in Section 6.

#### 2. Nonnegativity and Boundedness of Solutions

Let be the Banach space of continuous functions from to equipped with the sup-norm, where . The initial condition of (3) is given as where such that .

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 and ultimately bounded.*

* Proof. *First, we prove that is positive for . Assuming the contrary and letting be the first time such that , by the first equation of system (3), we have , and hence for , where is an arbitrarily small positive constant. This leads to a contradiction. It follows that is always positive. Using the same method, we have that is always positive. Further, from the second and the forth equations in (3), we have, respectively,
Using (4) we know that the and of (3) are positive for all .

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

Define and . By non-negativity of the solution, it follows that
This implies that is ultimately bounded, and so are , , , and .

This completes the proof of Theorem 1.

The basic reproduction number [24] for system (3) is 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 units of time and produces infected grouse through feeding and infected ticks through cofeeding. Each infected grouse lives for units of time and produces infected ticks.

System (3) has always an infection-free equilibrium , that is, grouse and ticks both at their carrying capacity with no disease present, where , , and an infection equilibrium , when , that is, ticks, grouse, and disease coexisting together, where is given by the following cubic: where Hence, we can write all of the densities in terms of the infectious tick density, ,

#### 3. Global Stability of the Infection-Free Equilibrium of Model **(3)**

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

Theorem 2. *If , then the infection-free equilibrium is globally asymptotically stable for any time delay and .*

* Proof. *Consider a Lyapunov functional:
with
Calculate the time derivative of :
Since , , it follows that
Obviously ensures that . It is clear that and if and only if , , , and . Hence, it follows from the Lyapunov-LaSalle invariance principle that the infection-free equilibrium is globally asymptotically stable when .

This completes the proof of Theorem 2.

#### 4. Dynamical Analysis for the Infection Equilibrium

System (3) has an infection equilibrium when . We consider two situations: (a) , and (b) , .

##### 4.1. Global Stability of the Infection Equilibrium When and ss

In this subsection, the situations and are considered.

Theorem 3. *If , then the infection equilibrium is globally asymptotically stable for and .*

* Proof. *Define a Lyapunov functional:
where
The derivatives of are given by
Hence, we obtain
Since , , and
it follows that
Using equality
we obtain
where
Since the function
is always nonpositive for any function , and if and only if . Therefore, we have .

It is easy to see that . It is clear that if and only if , , , . Hence, it is also known from stability theorem [25] that the infection equilibrium is globally asymptotically stable for any time delay under the condition .

This completes the proof of Theorem 3.

##### 4.2. Global Stability of the Infection Equilibrium When and

In this subsection, the situations and are considered.

Theorem 4. *If , then the infection equilibrium is globally asymptotically stable for and .*

* Proof. *Define a Lyapunov functional:
where
The derivatives of are given by
Hence, we obtain
Using equality
we obtain
where
Since the function
is always non-positive for any function , and if and only if . Therefore, we have .

It is easy to see that . It is clear that if and only if , , , . Hence, it is also known from stability theorem [25] that the infection equilibrium is globally asymptotically stable for any time delay under the condition .

This completes the proof of Theorem 4.

#### 5. Permanence

System (3) is said to be uniformly persistent if there exists a constant such that any solution of (3) satisfies

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

Theorem 5. *System (3) is uniformly persistent if it satisfies .*

* Proof . *Let
which is relatively closed in .

Now we prove that is positively invariant for . By the second and forth equations of (3), we have
Since and , it follows from (35) that
Thus, is positively invariant for .

We set
We claim that
Assuming , for all , it suffices to show that , for all . If it is not true, then there exists such that either (a) , or (b) , .

For case (a), from the forth equation of (3), we have
Since , we have
By the invariant of , we have , for all . Then from the second equation of (3), we have , for all.

From the first equation of (3), we obtain , thus . By the invariant of , we have .

Then from the second equation of (3), we obtain ; by the invariant of , we have , for all . From the first and third equations of (3), we obtain ; thus . By the invariant of , we have , , which contradicts the assumption that , for all . Similarly, we can obtain a contradiction for case (b). This proves the claim (38).

Let , where is the global attractor of restricted to . We show that . In fact, from and the first and third equations of (3), we have and . Thus, is the isolated invariant set in .

Next, we show that . If this is not true, then there exists a solution such that
For any sufficiently small constant , there exists a positive constant such that
For the constant given above, it follows from the second and forth equations of (3) that
If , as , then by a standard comparison argument and the nonnegativity, the solution of the following monotone system
with the initial condition , , for all converges to as well. Thus, , where is defined by
Differentiating with respect to time gives
Because , we have for a sufficiently small . Therefore, goes to either infinity or a positive number , which leads to a contradiction with . Thus we have .

Define by
It is clear that and . Thus by [20, Theorem 3], we have
for some constant . Let , where is the constant such that and . We showed that , , , and .

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 3–5. 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 .

For system (3), we choose , , , and , and other values are shown in Table 1. By computation, we obtain , and the infection equilibrium does not exist. By the analysis of Section 2, we know that the infection-free equilibrium is globally asymptotically stable (see Figure 1). Furthermore, we choose , , , and , and other values are shown in Table 1. In Figure 2, we know that under the condition indicated in the Theorem 3 with the infection equilibrium will be globally asymptotically stable. In Figure 3, we know that under the condition indicated in the Theorem 4 with , the infection equilibrium will be globally asymptotically stable with and . In Figure 4, we know that under the condition indicated in Theorem 5 with , the infection equilibrium will be globally asymptotically stable with and .

#### 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 of system (3) is globally asymptotically stable when . If then the infection equilibrium 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 .

The reproductive number plays a crucial role for epidemic dynamics. Considering the basic reproductive number 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 with two time delays; from the numerical simulations, it is shown that the infection equilibrium is globally asymptotically stable when ; 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).

#### References

- R. Norman, D. Ross, M. K. Laurenson, and P. J. Hudson, “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*, vol. 48, no. 2, pp. 119–134, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Norman, R. G. Bowers, M. Begon, and P. J. Hudson, “Persistence of tick-borne virus in the presence of multiple host species: tick reservoirs and parasite mediated competition,”
*Journal of Theoretical Biology*, vol. 200, no. 1, pp. 111–118, 1999. View at Publisher · View at Google Scholar · View at Scopus - T. Caraco, G. Gardner, W. Maniatty, E. Deelman, and B. K. Szymanski, “Lyme disease: self-regulation and pathogen invasion,”
*Journal of Theoretical Biology*, vol. 193, no. 4, pp. 561–575, 1998. View at Publisher · View at Google Scholar · View at Scopus - R. Rosà, A. Pugliese, R. Norman, and P. J. Hudson, “Thresholds for disease persistence in models for tick-borne infections including non-viraemic transmission, extended feeding and tick aggregation,”
*Journal of Theoretical Biology*, vol. 224, no. 3, pp. 359–376, 2003. View at Publisher · View at Google Scholar · View at MathSciNet - S. Sandberg, T. E. Awerbuch, and A. Spielman, “A comprehensive multiple matrix model representing the life cycle of the tick that transmits the agent of Lyme disease,”
*Journal of Theoretical Biology*, vol. 157, no. 2, pp. 203–220, 1992. View at Google Scholar · View at Scopus - T. Awerbuch-Friedlander, R. Levins, and M. Predescu, “The role of seasonality in the dynamics of deer tick populations,”
*Bulletin of Mathematical Biology*, vol. 67, no. 3, pp. 467–486, 2005. View at Publisher · View at Google Scholar · View at MathSciNet - P. J. Hudson, R. Norman, M. K. Laurenson et al., “Persistence and transmission of tick-borne viruses: ixodes ricinus and louping-ill virus in red grouse populations,”
*Parasitology*, vol. 111, pp. S49–S58, 1995. View at Google Scholar · View at Scopus - P. J. Hudson, “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.
- H. Zhang and L. Chen, “A model for two species with stage structure and feedback control,”
*International Journal of Biomathematics*, vol. 1, no. 3, pp. 267–286, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. Song and H. Guo, “Global stability of a stage-structured predator-prey system,”
*International Journal of Biomathematics*, vol. 1, no. 3, pp. 313–326, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. G. Mwambi, J. Baumgartner, and K. P. Hadeler, “Ticks and tick-borne diseases: a vector-host interaction model for the brown ear tick (Rhipicephalus appendiculatus),”
*Statistical Methods in Medical Research*, vol. 9, no. 3, pp. 279–301, 2000. View at Google Scholar · View at Scopus - M. Labuda, L. D. Jones, T. Williams, V. Danielova, and P. A. Nuttall, “Efficient transmission of tick-borne encephalitis virus between cofeeding ticks,”
*Journal of Medical Entomology*, vol. 30, no. 1, pp. 295–299, 1993. View at Google Scholar · View at Scopus - L. Gilbert, L. D. Jones, M. K. Laurenson, E. A. Gould, H. W. Reid, and P. J. Hudson, “Ticks need not bite their red grouse hosts to infect them with louping ill virus,”
*Proceedings of the Royal Society B*, vol. 271, no. 4, pp. S202–S205, 2004. View at Google Scholar · View at Scopus - L. Gilbert, R. Norman, K. M. Laurenson, H. W. Reid, and P. J. Hudson, “Disease persistence and apparent competition in a three-host community: an empirical and analytical study of large-scale, wild populations,”
*Journal of Animal Ecology*, vol. 70, no. 6, pp. 1053–1061, 2001. View at Publisher · View at Google Scholar · View at Scopus - Y. Kuang,
*Delay Differential Equations with Applications in Population Dynamics*, vol. 191 of*Mathematics in Science and Engineering*, Academic Press, Boston, Mass, USA, 1993. View at MathSciNet - A. Korobeinikov, “Lyapunov functions and global properties for SEIR and SEIS epidemic models,”
*Mathematical Medicine and Biology*, vol. 21, no. 2, pp. 75–83, 2004. View at Google Scholar · View at Scopus - A. Korobeinikov and P. K. Maini, “A Lyapunov function and global properties for SIR and SEIR epidemiological models with nonlinear incidence,”
*Mathematical Biosciences and Engineering*, vol. 1, no. 1, pp. 57–60, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Korobeinikov, “Global properties of infectious disease models with nonlinear incidence,”
*Bulletin of Mathematical Biology*, vol. 69, no. 6, pp. 1871–1886, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. Huang, Y. Takeuchi, W. Ma, and D. Wei, “Global stability for delay SIR and SEIR epidemic models with nonlinear incidence rate,”
*Bulletin of Mathematical Biology*, vol. 72, no. 5, pp. 1192–1207, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. C. McCluskey, “Complete global stability for an SIR epidemic model with delay—distributed or discrete,”
*Nonlinear Analysis: Real World Applications*, vol. 11, no. 1, pp. 55–59, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. L. Smith and X.-Q. Zhao, “Robust persistence for semidynamical systems,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 47, no. 9, pp. 6169–6179, 2001. View at Publisher · View at Google Scholar · View at MathSciNet - X.-Q. Zhao,
*Dynamical Systems in Population Biology*, CMS Books in Mathematics/Ouvrages de Mathématiques de la SMC, 16, Springer, New York, NY, USA, 2003. View at MathSciNet - H. R. Thieme, “Persistence under relaxed point-dissipativity (with application to an endemic model),”
*SIAM Journal on Mathematical Analysis*, vol. 24, no. 2, pp. 407–435, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - O. Diekmann, J. A. P. Heesterbeek, and J. A. J. Metz, “On the definition and the computation of the basic reproduction ratio ${R}_{0}$ in models for infectious diseases in heterogeneous populations,”
*Journal of Mathematical Biology*, vol. 28, no. 4, pp. 365–382, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. K. Hale and S. M. Verduyn Lunel,
*Introduction to Functional-Differential Equations*, vol. 99 of*Applied Mathematical Sciences*, Springer, New York, NY, USA, 1993. View at MathSciNet