Research Article | Open Access

Ganga Ram Phaijoo, Dil Bahadur Gurung, "Modeling Impact of Temperature and Human Movement on the Persistence of Dengue Disease", *Computational and Mathematical Methods in Medicine*, vol. 2017, Article ID 1747134, 9 pages, 2017. https://doi.org/10.1155/2017/1747134

# Modeling Impact of Temperature and Human Movement on the Persistence of Dengue Disease

**Academic Editor:**Gul Zaman

#### Abstract

Dengue is a vector-borne infectious disease endemic in many parts of the world. The disease is spreading in new places due to human movement into the dengue disease supporting areas. Temperature is the major climatic factor which affects the biological processes of the mosquitoes and their interaction with the viruses. In the present work, we propose a multipatch model to assess the impact of temperature and human movement in the transmission dynamics of dengue disease. The work consists of system of ordinary differential equations that describe the transmission dynamics of dengue disease between humans and mosquitoes. Human population is divided into four classes: susceptible, exposed, infectious, and recovered. Mosquito population is divided into three classes: susceptible, exposed, and infectious. Basic reproduction number of the model is obtained using Next-Generation Matrix method. The qualitative analysis of the model is made in terms of the basic reproduction number. Parameters used in the model are considered temperature dependent. Dynamics of vector and host populations are investigated with different human movement rates and different temperature levels. Numerical results show that proper management of human movement between patches helps reducing the burden of dengue disease. It is also seen that the temperature affects the transmission dynamics of the disease significantly.

#### 1. Introduction

Dengue disease is a vector-borne viral infection that usually occurs in tropical and subtropical countries. Nowadays, the disease has been recognized in over 100 countries and an estimated 50–100 million dengue cases occur annually. The disease is threatening about 40% of the world’s population [1]. The disease is transmitted by the bites of infected mosquitoes named* Aedes aegypti* and* Aedes albopictus*. Four serologically different viruses* DEN* 1–*DEN* 4 cause the disease. Infection from one serotype grants life-long immunity to that strain and also shows temporary cross-immunity to the others. However, ultimately the recovered patient will become more susceptible to the other three forms [2, 3]. We assume that the infectivity of the mosquitoes ends with their death since they have a short lifespan.

Mathematical modeling has become an interesting tool for the understanding of epidemic diseases and to propose strategies to control the transmission of the disease. In 1927, Kermack and McKendrick developed an SIR model to describe epidemic diseases [4]. The model is being followed by many researchers to investigate the transmission dynamics of infectious diseases with some modifications. Esteva and Vargas proposed the SIR model to address dengue disease transmission considering constant and variable human populations [5, 6]. Since then many mathematical models have been proposed to study different aspects of dengue disease transmission. Authors [7, 8] discussed the role of awareness in controlling dengue disease transmission. Pinho et al. used a mathematical model for dengue disease transmission analysis comparing two dengue epidemics [9]. Authors [10–14] focused on incubation period to study dengue disease transmission. Sardar et al. discussed a mathematical model of dengue disease transmission with memory. They incorporated memory in the model by using a fractional differential operator [15].

Dengue infections are sensitive to the climate. Changing climate factors affect the potential for the geographic spread and future dengue disease. One of the principal determinants of* Aedes* mosquitoes’ survival is temperature which has been associated with seasonal changes. The temperature plays an important role in the life cycle and behavior of the mosquitoes. So, mathematical studies have been made to understand the role of temperature in transmission dynamics of dengue disease. Brady et al. modeled* Aedes aegypti* and* Aedes albopictus* survival at different temperature levels in laboratory and field settings [16]. Liu-Helmersson et al. studied the vectorial capacity of* Aedes aegypti* and made investigations on the effects of temperature and implications for global dengue epidemic potential [17]. Polwiang discussed the seasonal basic reproduction number of dengue and impacts of climate on transmission of the disease [18].

Travel and transport contribute to the spread of infectious diseases like dengue in new places. So, one of the major factors contributing to the reemergence of infectious diseases is human movement from one place to the other. They help the disease in expanding their geographic range. Many mathematical models are proposed to address the impact of movement of humans and dispersal of vectors in the transmission dynamics of infectious diseases. Wang and Zhao discussed an epidemic model in patchy environment to describe the dynamics of disease spread among patches due to population dispersal [19]. An epidemic model was proposed by Wang and Mulone to describe the dynamics of disease spread between two patches due to population dispersal. They proved that reproduction number is a threshold of the uniform persistence and disappearance of the disease [20]. Arino and van den Driessche gave some analytical results for a model that describes the propagation of a disease in a population of individuals who travel between patches [21]. Hsieh et al. proposed a multipatch epidemic model to study the impact of travel of humans on the spread of disease between patches with different level of disease prevalence [22]. Cosner et al. investigated the effects of human movement on the persistence of vector-borne diseases [23]. Dynamics of malaria disease was studied in patchy environment by Auger et al. They generalized Ross-Macdonald model to -patches to describe the transmission dynamics of the disease [24]. Lee and Castillo-Chavez [25] and Phaijoo and Gurung [26] discussed dengue disease transmission dynamics in patchy environment.

Temperature influences dengue disease dynamics by affecting dynamics of mosquitoes and vector host interactions. Dengue disease has been spreading rapidly to new areas via human movement. So, in the present work, we propose a multipatch SEIR-SEI model of dengue disease considering the temperature dependent model parameters to study the impact of temperature and movement of humans on the persistence of dengue disease. We have considered different temperature levels and different movement rates in different patches. Basic reproduction number of the individual patches and a combined basic reproduction number are computed. Local stability of disease-free equilibrium point is proved by basic reproduction number.

#### 2. Model Description and Formulation

The total human (host) population in each patch is subdivided into the classes: susceptible , exposed , infectious , and recovered . Mosquito (vector) population is subdivided into the classes: susceptible , exposed , and infectious , . Recovered class in the mosquito population is not considered due to their short lifespan.

The recruitment rate of host population is . Susceptible hosts get infected by infectious vectors at the rate , where is the biting rate and is the transmission probability from vector to host. The exposed host becomes infectious at the rate after developing the symptoms. Infectious host recovers at the rate . Host dies naturally with the rate . In case of vector population, susceptible vector gets infected by interaction with infectious hosts at the rate . The exposed vector becomes infectious at the rate developing the symptoms of the disease. is the natural death rate of vectors.

Here, the model parameters , , , , and are temperature dependent. The temperature dependency relations are discussed below [17, 18]:We consider human movement between the patches. Human of patch moves to patch at the rate and the human of patch moves to patch at the rate . Here and represents , , , and , respectively for susceptible, exposed, infectious, and recovered human movement rates.

The system of ordinary differential equations describing the present multipatch model [22] is given by where

#### 3. Equilibrium Point and Stability Analysis

In this section, we find disease-free equilibrium (DFE) of the system of (2) and discuss its stability. An equilibrium is said to be disease-free if there is no infective population in both host and vector populations.

Theorem 1. *Model (2) has a unique disease-free equilibrium.*

*Proof. *In disease-free situation, , and other variables , , , , and for .

System of (2) for host population in disease-free situation can be written aswhereSystem of (2) for vector population in disease-free situation can be written aswhereSince the matrix has all off-diagonal entries negative and each column sum is positive, is nonsingular -matrix. Matrix is irreducible as the matrix has nonzero diagonal elements. So, must have positive inverse [27]. Hence, the system of (4) has a unique solution .

Again, matrix is a diagonal matrix with positive diagonal elements. So, exists with positive diagonal elements. Hence, the system of (6) has a unique solution and system (2) has a unique disease-free equilibrium.

*Basic Reproduction Number*. When a typical infective is introduced into a completely susceptible population, the expected number of new infections produced by this single infective during its infectious period is called basic reproduction number.

To find the mathematical expression for the basic reproduction number, we order the variables related to the infections by , , , . We use Next-Generation Matrix method [28, 29] to find transmission matrix, , and transition matrix, , and we find basic reproduction number asFor the system of (2),Here,Matrices and are irreducible nonnegative -matrices. So, and exist. Also, , , , and are diagonal matrices. So, their inverses exist. Hence, exists and basic reproduction number, , is given by

Theorem 2 (local stability). *The disease-free equilibrium point of the system of (2) is locally asymptotically stable if and unstable if .*

*Proof. *Jacobian matrix for the system of (2) at disease-free equilibrium is given by Matrix is triangular matrix. So, the stability of the system of (2) depends on matrices and . Matrix can be written as Matrices and (defined in Theorem 1) are nonsingular -matrices. So, the matrix has eigenvalues with negative real parts [27]. Hence, the stability of model (2) depends on the matrix only. Here, matrix is nonnegative matrix and is a nonsingular -matrix. So, the matrix will have eigenvalues with negative real parts if [29]; that is, . Thus, the disease-free equilibrium is locally asymptotically stable if . If , then . Which shows that at least one eigenvalue lies in right half plane. So, the disease-free equilibrium is unstable if .

When only the two patches are taken into the consideration, the basic reproduction is given bywhereHere, is the basic reproduction number of patch 1 and is the basic reproduction number of patch 2.

#### 4. Simulations and Discussion

Temperature plays a significant role in the transmission dynamics of dengue disease. Small change in temperature can affect whole dynamics of the disease. Human movement from one place to the other helps spreading disease into new areas and influences the prevalence of the disease. Thus, both temperature and human movement have a significant influence on the transmission dynamics of dengue disease. For the simulation purpose, the following data are used: , , . The parameters , , , , are considered temperature dependent following [17].

Figures 1–4 are drawn with different temperature levels to investigate the dynamics of infectious hosts of patch 1 and patch 2. Figures 1 and 2 are drawn when there is no human movement between the patches. Here, patch 1 is high disease prevalent compared to patch 2. With the human movement, it is seen that infectious host population is decreased in patch 1 and the population is increased in patch 2. Thus, the human movement can cause the low endemic patch to be high endemic and high endemic patch to be low endemic patch (Figures 1–4). Also, the figures show that the burden of disease is increased with temperature. Again, the number of infectious hosts is seen increasing initially due to interaction of hosts with infectious vectors. Afterwards the number is seen decreased due to recovery from the disease and natural death of hosts (Figures 1–4).

Figures 5 and 6 show the impact of movement of infectious and susceptible hosts, respectively, on basic reproduction number . Infectious host can infect the mosquitoes of the patch where the hosts are travelling and the susceptible hosts can get infected of the disease from the mosquitoes of the patch where the hosts have travelled. It is observed that movement of both infectious and susceptible hosts from low prevalent patch to the high prevalent patch increases the endemic level of the disease. But their movement from high prevalent patch to the low prevalent patch decreases the endemic level of the disease.

Temperature has a significant influence on basic reproduction number (Figures 7 and 8). In patch 1, the prevalence of disease is seen increasing with temperature and the maximum disease prevalence has occurred at C temperature as in [17]. In case of patch 2 where average temperature range is C to C, disease prevalence increases with the increase in temperature and the maximum disease prevalence has occurred at C.

##### 4.1. Dynamics with Unidirectional Movement

In this section, we investigate the impact of host movement in one direction only with different temperature levels. Figures 9–12 show the dynamics of infectious host population of patch 1 and patch 2 when there is host movement from patch 1 to patch 2 only or patch 2 to patch 1 only. When the hosts from patch 2 are not allowed to move to patch 1 (Figures 9 and 10), burden of disease is decreased in patch 1 and increased in patch 2. When the hosts from patch 1 are restricted to travel to patch 2, the burden of disease is increased in patch 1 and decreased in patch 2 (Figures 11 and 12). In each case, the dynamics of infectious hosts are seen temperature dependent. Disease prevalence is observed increasing with temperature. Thus, movement of hosts can cause the patch to be less disease prevalent (Figures 9 and 12) and more disease prevalent (Figures 10 and 11).

When only the hosts from patch 2 are allowed to move to patch 1, basic reproduction number of patch 1 increases and that of patch 2 decreases with the increase in movement rate (Figure 13). Also, basic reproduction number of patch 1 decreases and that of patch 2 increases when only the hosts from patch 1 are allowed to move to patch 2 (Figure 14). So, the host population can be moved from one patch to the other to bring the disease under control.

#### 5. Conclusion

Temperature plays a significant role in dynamics of dengue disease transmission. It affects the lifecycle and biting behavior of mosquitoes. Human movements contribute in spreading the disease in new places. We have proposed multipatch model of dengue disease with the human movement between patches considering temperature dependent model parameters. In the present work, we explored the impact of temperature and host movement between patches on the transmission dynamics of dengue disease. We have investigated the stability of disease-free equilibrium point. It is observed that the point is locally asymptotically stable when basic reproduction number and unstable when . Simulated results show that basic reproduction number depends on temperature and host movement. The prevalence of disease can increase or decrease with temperature and mobility of hosts from one patch to the other. Present work shows that the burden of the disease can be reduced by managing the host movement and the temperature can enhance the strength of the disease. These pieces of information can be helpful to the concerned authorities to bring dengue disease under control.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### References

*Global strategy for dengue prevention and control 2012-2020*, World Health Organization, 2012.- N. M. Ferguson, C. A. Donnelly, and R. M. Anderson, “Transmission dynamics and epidemiology of dengue: Insights from age-stratified sero-prevalence surveys,”
*Philosophical Transactions of the Royal Society B: Biological Sciences*, vol. 354, no. 1384, pp. 757–768, 1999. View at: Publisher Site | Google Scholar - D. J. Gubler, “Dengue and dengue hemorrhagic fever,”
*Clinical Microbiology Reviews*, vol. 11, no. 3, pp. 480–496, 1998. View at: Google Scholar - W. O. Kermack and A. G. McKendrick, “A contribution to the mathematical theory of epidemics,”
*Proceedings of the Royal Society A*, vol. 115, no. 772, pp. 700–721, 1927. View at: Publisher Site | Google Scholar - L. Esteva and C. Vargas, “Analysis of a dengue disease transmission model,”
*Mathematical Biosciences*, vol. 150, no. 2, pp. 131–151, 1998. View at: Publisher Site | Google Scholar | Zentralblatt MATH - L. Esteva and C. Vargas, “A model for dengue disease with variable human population,”
*Journal of Mathematical Biology*, vol. 38, no. 3, pp. 220–240, 1999. View at: Publisher Site | Google Scholar | MathSciNet - S. Gakkhar and N. C. Chavda, “Impact of Awareness on the Spread of Dengue Infection in Human Population,”
*Applied Mathematics*, vol. 04, no. 08, pp. 142–147, 2013. View at: Publisher Site | Google Scholar - G. R. Phaijoo and D. B. Gurung, “Mathematical model of dengue disease transmission dynamics with control measures,”
*Journal of Advances in Mathematics and Computer Science*, vol. 23, no. 3, pp. 1–12, 2017. View at: Google Scholar - S. T. Pinho, C. P. Ferreira, L. Esteva, F. R. Barreto, V. C. Morato e Silva, and M. G. Teixeira, “Modelling the dynamics of dengue real epidemics,”
*Philosophical Transactions of the Royal Society of London. Series A. Mathematical, Physical and Engineering Sciences*, vol. 368, no. 1933, pp. 5679–5693, 2010. View at: Publisher Site | Google Scholar | MathSciNet - M. Chan and M. A. Johansson, “The Incubation Periods of Dengue Viruses,”
*PLoS ONE*, vol. 7, no. 11, Article ID e50972, 2012. View at: Publisher Site | Google Scholar - P. Pongsumpun, “Transmission model for dengue disease with and without the effect of extinsic incubation period,”
*KMITL Science and Technology Journal*, vol. 6, pp. 74–82, 2006. View at: Google Scholar - P. Pongsumpun, “Mathematical model of dengue disease with the incubation period of virus,”
*World Academy of Science, Engineering and Technology*, vol. 44, pp. 328–332, 2008. View at: Google Scholar - S. Side and M. S. M. Noorani, “SEIR model for transmission of dengue fever in selangor Malaysia,”
*International Journal of Modern Physics*, vol. 9, pp. 380–389, 2012. View at: Publisher Site | Google Scholar - S. Side and S. M. Noorani, “A SIR model for spread of dengue fever disease (simulation for South Sulawesi, Indonesia and Selangor, Malaysia),”
*World Journal of Modelling and Simulation*, vol. 9, no. 2, pp. 96–105, 2013. View at: Google Scholar - T. Sardar, S. Rana, and J. Chattopadhyay, “A mathematical model of dengue transmission with memory,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 22, no. 1-3, pp. 511–525, 2015. View at: Publisher Site | Google Scholar | MathSciNet - O. J. Brady, M. A. Johansson, C. A. Guerra et al., “Modelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settings,”
*Parasites and Vectors*, vol. 6, no. 1, article no. 351, 2013. View at: Publisher Site | Google Scholar - J. Liu-Helmersson, H. Stenlund, A. Wilder-Smith, and J. Rocklöv, “Vectorial capacity of
*Aedes aegypti*: effects of temperature and implications for global dengue epidemic potential,”*PLoS ONE*, vol. 9, no. 3, Article ID e89783, 2014. View at: Publisher Site | Google Scholar - S. Polwiang, “The seasonal reproduction number of dengue fever: Impacts of climate on transmission,”
*PeerJ*, vol. 2015, no. 7, Article ID e1069, 2015. View at: Publisher Site | Google Scholar - W. Wang and X.-Q. Zhao, “An epidemic model in a patchy environment,”
*Mathematical Biosciences*, vol. 190, no. 1, pp. 97–112, 2004. View at: Publisher Site | Google Scholar | MathSciNet - W. Wang and G. Mulone, “Threshold of disease transmission on a patch environment,”
*Journal of Mathematical Analysis and Applications*, vol. 285, no. 1, pp. 321–335, 2003. View at: Publisher Site | Google Scholar | MathSciNet - J. Arino and P. van den Driessche, “A multi-city epidemic model,”
*Mathematical Population Studies. An International Journal of Mathematical Demography*, vol. 10, no. 3, pp. 175–193, 2003. View at: Publisher Site | Google Scholar | MathSciNet - Y.-H. Hsieh, P. van den Driessche, and L. Wang, “Impact of travel between patches for spatial spread of disease,”
*Bulletin of Mathematical Biology*, vol. 69, no. 4, pp. 1355–1375, 2007. View at: Publisher Site | Google Scholar | MathSciNet - C. Cosner, J. C. Beier, R. S. Cantrell et al., “The effects of human movement on the persistence of vector-borne diseases,”
*Journal of Theoretical Biology*, vol. 258, no. 4, pp. 550–560, 2009. View at: Publisher Site | Google Scholar | MathSciNet - P. Auger, E. Kouokam, G. Sallet, M. Tchuente, and B. Tsanou, “The Ross-Macdonald model in a patchy environment,”
*Mathematical Biosciences*, vol. 216, no. 2, pp. 123–131, 2008. View at: Publisher Site | Google Scholar - S. Lee and C. Castillo-Chavez, “The role of residence times in two-patch dengue transmission dynamics and optimal strategies,”
*Journal of Theoretical Biology*, vol. 374, pp. 152–164, 2015. View at: Google Scholar - G. R. Phaijoo and D. B. Gurung, “Mathematical Study of Dengue Disease Transmission in Multi-Patch Environment,”
*Applied Mathematics*, vol. 07, no. 14, pp. 1521–1533, 2016. View at: Publisher Site | Google Scholar - A. Berman and R. J. Plemmons,
*Nonnegative Matrices in Mathematical Sciences*, Academic Press, New York, NY, USA, 1979. View at: MathSciNet - O. Diekmann, J. A. P. Heesterbeek, and J. A. J. Metz, “On the definition and 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 Site | Google Scholar | MathSciNet - P. van den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,”
*Mathematical Biosciences*, vol. 180, pp. 29–48, 2002. View at: Publisher Site | Google Scholar | MathSciNet

#### Copyright

Copyright © 2017 Ganga Ram Phaijoo and Dil Bahadur Gurung. 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.