Computational and Mathematical Methods in Medicine

Volume 2017, Article ID 1747134, 9 pages

https://doi.org/10.1155/2017/1747134

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

Department of Natural Sciences (Mathematics), School of Science, Kathmandu University, Dhulikhel, Kavre, Nepal

Correspondence should be addressed to Ganga Ram Phaijoo; pn.ude.uk@maragnag

Received 15 April 2017; Revised 5 June 2017; Accepted 8 August 2017; Published 19 September 2017

Academic Editor: Gul Zaman

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.

#### 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).