key: cord-0009018-1e7qm4x2
authors: Mishra, Arti; Gakkhar, Sunita
title: Non-linear Dynamics of Two-Patch Model Incorporating Secondary Dengue Infection
date: 2017-11-30
journal: Int J Appl Comput Math
DOI: 10.1007/s40819-017-0460-z
sha: 1e94fc4a1744356edf7ac876a5295aa0294adce5
doc_id: 9018
cord_uid: 1e7qm4x2

In this paper, the impact of human migration on the dynamics of dengue epidemic has been discussed. The vector-host model considers two patches with different dengue serotype in each patch. The model considers the constant rate of migration in susceptible and recovered class from one patch to other. Recovered migrants from prior infection are exposed to secondary infection in the patch where different serotype is present. The basic reproduction number is computed and analyzed in terms of migration parameters. The model is analyzed for the existence and local stability of various equilibrium states in terms of migration parameters. The numerical simulations for the choice of relevant data from literature have been performed to verify analytical results and to further explore the dynamics of the system. The sensitivity analysis of basic reproduction number with respect to migration parameters is carried out. It is found that immigration in a patch increases the basic reproduction in respective patch and vice-versa. The basic reproduction number has been estimated for the two states of Brazil which verifies the occurrence of severe epidemic in one of the states of Brazil.

Dengue was confined to only nine countries before 1970. However, it has now been spread to more than 100 countries in the tropical, subtropical and temperate areas of North America, South America, Africa and Southeast Asia. According to one estimate, it affects 50-100 million people every year [1]. In fact, World Health Organization has declared dengue as the fastest communicable mosquito-borne disease in the world in terms of human morbidity and mortality [2, 3] . The causes of emergence and resurgence of dengue disease are attributable to multiple factors including urbanization, fast transportation, economic development and changes in human behavior. The fast increase in modes of transportation, their reach, speed and efficiency led to the global movement of human population [2, 4] . This may be responsible for the spread of infectious diseases from isolated locations to the new regions. It has been reported that another 100 million people travel to dengue affected areas where they not only increase the risk of contracting the disease but also spreading it further afield.

Dengue fever exhibits up to four closely related and distinct dengue viruses or serotypes. The origin of all four dengue serotypes is found to have in Asian forests and subsequently they have spread worldwide because of the migration of human and commerce [5] . The mosquitoes Aedes aegypti and Aedes albopictus are the main carriers of these viruses. It is transmitted through infected mosquito bites. Once a person becomes infected by a serotype, he/she will never be reinfected by that serotype due to induced lifelong immunity against it. However, only temporal cross immunity is developed against the other serotypes [6] . Therefore, the recovered individuals from a serotype may get secondary infection when exposed to a different serotype. Primary infections are generally asymptomatic and nonfatal, while secondary infections sometimes lead to the fatal forms of dengue infection, namely, dengue hemorrhagic fever (DHF) or dengue shock syndrome (DSS) [1, 7] . Mathematical modeling has been applied to understand and control the many vector-borne diseases like dengue, malaria, leishmaniasis, chagas diseases etc [8] [9] [10] [11] [12] [13] . Some mathematical models on dengue dynamics incorporating primary and secondary infection have been proposed and analyzed by various investigators [11] [12] [13] [14] [15] [16] .

Due to short flying range of mosquitoes, the role of human migration becomes crucial in dissemination of dengue serotypes from one part to other parts of the region [17] . Some studies on dengue have also confirmed that human migration is the main cause of its spatial spread [18, 19] . In particular, human movement was identified as an important factor for the spread of dengue serotype-1 from the southern region of Vietnam to its northern and central regions [20] . Also, human migration from one patch (area) to another increases the chances of secondary infection [6, 21] . During 1977 During -1978 , the dengue epidemic in Cuba occurred due to DEN-1 serotype with not a single case of DHF [22] . However, human migration from endemic countries (with DEN-2 serotype) to municipalities in Cuba has been reported to be responsible for DHF/DSS cases of subsequent epidemics in 1981 and 1997 [22, 23] . These facts exemplify the two largest outbreaks of secondary infection cases by DEN-2 in Cuba in 1981 and 1997 [22] [23] [24] . Further, human migration may be one of the reasons for the spread of secondary infection.

A very limited but relevant metapopulation epidemic models involving movement of individuals between discrete spatial patches have been proposed by investigators [25] [26] [27] [28] [29] [30] [31] .

In particular, such compartment models for multi-patches have been discussed for measles, influenza and SARS [25, 28, 29, 31] . Prosper et al. [32] have discussed Ross-Macdonald type two patch SI metapopulation model for malaria by considering different degrees of transmission in patches. Gao and Ruan [27] have derived a multi-patch model to study the impact of population dispersal on spatial spread of malaria between patches. A limited work on metapopulation model incorporating host-vector dynamics has been carried out. A two patch vector-host model incorporating secondary infection of dengue disease has been discussed by Gakkhar and Mishra [33] . They have assumed no vector migration and human migration is allowed only in one direction.

In this paper, a two-patch vector-host model for the dengue epidemic incorporating secondary infection has been proposed and analyzed. The model considers human migration between two discrete patches having different dengue serotypes. Due to imperfect cross immunity towards heterologous serotype, the recovered migrants of patch-1 may get sec-ondary infection in patch-2. The impact of human migration on persistence or extinction of disease in patches has been discussed by considering all possible combinations of migration parameters. Stability analysis and numerical simulations have been carried out for various equilibrium states. Further, sensitivity indices of basic reproduction numbers for migration parameters have been computed. In the case study, the estimation of basic reproduction has also been performed for the two states of Brazil. The conclusion is given in the last section.

Keeping the above facts into consideration, let there be two distinct serotypes of dengue namely serotype-1 (DEN-1/DEN-3/DEN-4) and serotype-2 (DEN-2) prevalent in two discrete patches: patch-1 and patch-2 respectively. Human (host) population is compartmentalized into S i , I i and R i for susceptible, infected and recovered individuals respectively in the ith patch, i = 1, 2. It is assumed that the infected individuals are unable to migrate due to disease-induced weakness. The susceptible individuals of patch-1 (patch-2) may become primary infected with serotype-2 (serotype-1) when they migrate to patch-2 (patch-1). The recovered migrants from patch-1 to patch-2 may become secondary infected when exposed to serotype-2 and form a new compartment of secondary infection I 12 in patch-2. As, the serotype DEN-2 is more associated with secondary infection [7, 23, 24, 34, 35] , the almost perfect cross immunity has been assumed for the recovered migrants from patch-2 to patch-1. The probability of getting secondary infection from serotype-1 is assumed to be negligibly small in patch-1. Accordingly, the recovered migrants from patch-2 to patch-1 will not get secondary infection. The secondary infected population I 21 is not being considered (I 21 ∼ 0) in patch-1. Further, this simplifying assumption reduces the dimension of the model.

Let ω be the constant rate of recruitment in susceptible class and μ be the natural death rate of host in both the patches. The disease-induced death rate e is considered only in the secondary infection compartment of patch-2. Let β i and γ i ; i = 1, 2 be the rates of transmission of infection and recovery of the host population respectively in the ith patch. The ρ proportion of recovered migrants from serotype-1 become secondary infected by serotype-2 at the rate β 12 in patch-2. Let m 1 and m 2 be the rates of migration from patch-1 to patch-2 and vice-versa respectively.

LetŪ i (t) andV i (t) be the susceptible and infected vector population in the ith patch. Due to short life span of vector, once infected they never recover. Accordingly, no recovery class has been considered for them. Further, let ω 1 be the recruitment rate of mosquitoes in the absence of vertical transmission. Let μ 1 and μ 2 be the constant natural death rates of mosquitoes in patch-1 and patch-2 respectively. The transmission rate of infection to mosquitoes by infected host is assumed to be σ . The mosquitoes are not migrating and their dynamics remains the same irrespective of the patch. Considering all parameters to be positive, the following model has been formulated:

The model is associated with following non-negative initial conditions:

The schematic diagram has been drawn in Fig. 1 showing the transmission dynamics of disease in two patches. Note that for simplicity, the mortality rates of mosquitoes are considered to be equal in both the patches i.e. (μ 1 = μ 2 ). Consider the following non-dimensional parameters: 

The dimensionless model is given as:

The model is associated with following non-negative initial conditions:

Consider a system of differential equations in R n + as

is positive invariant for the non-linear system (12)- (22) . By adding the host dynamics from the model (12)- (22) gives,

Using standard comparison theorem [36] ,

Again, by adding vector dynamics in from the model (12)- (22) gives,

Using standard comparison theorem [36] ,

Also,Ṅ < 0 andṀ < 0 for N > 2 and M > 2 respectively. This shows that the solutions of the system (12)- (22) converge towards the set Ω.

The Lemma 1 shows that the all solutions of the non-linear system (12)-(22) are non-negative and bounded. Therefore, the model is mathematically as well as biologically well behaved. In next subsection, the existence of equilibria for the model (12)- (22) has been discussed.

The non-linear system of Eqs. (12)- (22) have four equilibrium states (S 1 ,

1. The disease-free state (E 0 ) = (Ŝ 1 , 0, 0, 1, 0,Ŝ 2 , 0, 0, 0, 1, 0),

The expressions for non-zero state variables at equilibrium point E 1 are given as follows:

The expressions for non-zero state variables at equilibrium point E 2 are obtained as follows:S 1 =

4. The endemic state (E )=(S 1 , I 1 , R 1 , U 1 , V 1 , S 2 , I 2 , R 2 , I 12 , U 2 , V 2 ) exists under the following conditions:

where,

The expressions for the state variables at equilibrium level are omitted as they are lengthy and complex.

The basic reproduction number is computed by next generation approach [37] . The details of the new infections and transfer matrices are given in "Appendix". The basic reproduction number, R 0 (say) is given as

,

or 10 and R 01 are the basic reproduction numbers of patch-1 and patch-2 respectively.

Observe that there will be no effect of migration on the basic reproduction numbers R 10 and R 01 when p 2 = p 3 .

This section analyzes the stability of the non-linear system (12)- (22) for various equilibrium states by computing eigenvalues of the Jacobian matrix J [E]. The general matrix J [E] of the system is given in "Appendix".

Theorem 1 The disease-free state E 0 is locally asymptotically stable for R 2 10 < 1 and R 2 01 < 1.

Proof About disease-free state, seven of the eigenvalues of the Jacobian matrix (J [E 0 ]) are negative and computed as −1(multiplicity 4), −1 − p 2 − p 3 (multiplicity 2) and -1 -

The remaining four eigenvalues are given as:

Further simplifications give all eigenvalues with negative real part under condition (27) . Therefore, the disease-free state is locally asymptotically stable under condition (27) .

Further, for global stability of disease-free state, the following theorem is concluded:

The locally asymptotically stable state E 0 is also globally stable for

Proof For arbitrarily chosen positive constants A, B, C and D, consider the positive definite function L(I 1 , I 2 , V 1 , V 2 ) as:

Taking derivative of L(I 1 , I 2 , V 1 , V 2 ) with respect to t and its simplifications yield, 

(< 1) for positive p 2 and p 3 .

then the inequalities can be combined to give R 2 10 < C < α < 1; R 2 01 < D < β < 1;

The positive arbitrary constants C and D can now be chosen to satisfy the above inequality ensuring thatL(S, I 1 , I 2 , V 1 , V 2 ) is negative. Accordingly, the function L(S, I 1 ,

is a Lyapunov function for the condition (28) . At V 1 = 0, V 2 = 0, I 1 = 0 and I 2 = 0,L becomes zero. If V 1 = 0, V 2 = 0, I 1 = 0 and I 2 = 0, then {E 0 } is the only largest invariant set that contains a subset in which all these variables are zero. By applying LaSalle's invariance principle [38] , all trajectories in the closed set Ω approach the equilibrium point E 0 . Hence, the locally asymptotically stable disease-free state E 0 is also globally asymptotically stable. Hence, the result is proved.

It is observed that E 0 may still be globally stable even when one or both of the conditions in (28) are not satisfied. Further, when the state E 0 is unstable, existence of some other states may become possible. These possibilities are explored next.

It may be noted that the existence of E 1 requires R 2 10 > 1 [see condition (23) ] which implies the instability of disease-free state E 0 . 2 10 > 1 then ξ = (1 + p 2 + p 3 )( p 2 p 5 + p 1 p 9 + p 2 + p 9 + 2 p 2 p 9 )

Proof Assuming ξ > 1 leads to R 2 10 < 1 which is the contradiction. Therefore, ξ < 1 whenever R 2 10 > 1. (29), the state E 1 will be locally asymptotically stable when

Proof For the local stability of E 1 = (S 1 , I 1 , R 1 , U 1 , V 1 , S 2 , 0, R 2 , 0, 1, 0), the seven of the eigenvalues of the Jacobian matrix J [E 1 ] are computed as:

. Note that all of these seven eigenvalues are negative provided condition (30) is satisfied. The other four eigenvalues are the roots of the polynomial

where,

Now, substituting the values of state variables at E 1 and simplifying, it has been found that all the four roots of the polynomial are having negative real part for R 2 10 > 1, if the Routh-Hurwitz conditions satisfy. Hence, the state E 1 will be locally asymptotically stable if the Routh-Hurwitz conditions hold.

Remark 1 Accordingly, state E 1 is stable when R 2 01 < 1. However, when R 2 01 > 1, the state may be stable/unstable subject to the condition (30) . These are further explored in terms of migration:

• Observe that ξ = 1 in absence of migration ( p 2 = p 3 = 0). Therefore, in absence of migration the state E 1 will be locally asymptotically stable when R 2 01 < 1. • When migration is allowed only in patch-1 i.e. p 2 = 0 but p 3 > 0, it is observed that ξ = 1 and the state E 1 is locally asymptotically stable for R 2 01 < 1. • When p 2 > 0 and p 3 > 0 then by simplifying the expression for ξ , it may be noted that ξ still remains smaller than 1. Therefore, the state E 1 is again locally stable for R 2 01 < 1.

Remark 2 Stability of the state E 1 is possible even though R 2 01 > 1 since the stability condition ξ R 2 01 < 1 may still be satisfied for sufficiently small ξ < 1 in presence of migration.

Local stability of E 2 (disease in patch-2 state) state is discussed below:

Lemma 3 When R 2 01 > 1 then η = (1 + p 2 + p 3 )( p 3 + p 9 + 2 p 3 p 9 + p 4 p 9 + p 3 p 6 ) (1 + 2 p 3 )( p 9 + p 2 p 9 + p 3 p 9 + p 4 p 9 + p 2 p 4 p 9 ) < 1

Proof Assuming η > 1 leads to R 2 01 < 1 which is the contradiction. Therefore, η < 1 whenever R 2 01 > 1. Theorem 4 For η given in (31) , the state E 2 will be locally asymptotically stable when

Proof For the local stability of E 2 =(S 1 , 0, R 1 , 1, 0, S 2 , I 2 , R 2 , I 12 , U 2 , V 2 ), the seven eigenvalues of the 11 × 11 Jacobian matrix are given as follows:

and

. Note that all of these seven eigenvalues are negative for condition (32) . The other eigenvalues are the roots of following polynomial of degree four

Now putting the values of variables at state E 2 , all the four roots of the above polynomial are having negative real part for R 2 01 > 1 provided the Routh-Hurwitz conditions satisfy. Hence, the state E 2 will be locally asymptotically stable if the Routh-Hurwitz conditions hold.

Remark 3 Accordingly, state E 2 is stable when R 2 10 < 1. However, when R 2 10 > 1, the state may be stable/unstable. These are further explored in terms of migration:

• When both the patches are isolated i.e. p 2 = 0 and p 3 = 0 Assuming R 2 01 = p 4 p 9 1 + p 6 > 1, the state E 2 will be locally asymptotically stable when R 2 10 < 1. • When migration is allowed only in patch-2 i.e. p 3 =0 but p 2 > 0, then it is observed that η = 1 gives the local stability of E 2 for R 2 10 < 1. • When p 3 > 0 and p 2 > 0 then by simplification yields η < 1. Therefore, the state E 2 is again locally stable for R 2 10 < 1. Remark 4 Stability of the state E 2 is possible even though R 2 10 > 1 since the stability condition η R 2 10 < 1 may still be satisfied for sufficiently small η < 1 in presence of migration. Remark 5 When R 2 10 > 1 and R 2 01 > 1, the states E 1 and E 2 both exist. when R 2 10 η < 1 or R 2 01 ξ < 1, the states E 1 or E 2 respectively are locally stable while the state E * does not exist. The state E * exists for R 2 10 η > 1 and R 2 01 ξ > 1.

For the local stability of the endemic state E , the three of the eigenvalues of Jacobian matrix are −1(multiplicity two) and −1 − p 6 − p 8 . Two of the eigenvalues are the negative roots of second degree polynomial

Further, the remaining eigenvalues are the roots of the polynomial of degree 6:

The expressions for the coefficients of above six degree polynomial omitted from the text as they are lengthy and complex. The numerical simulations have been performed for the stability of the endemic state in the next section.

Numerical simulation of the system (1)- (11) are carried out for the data given in Table 1 . For different combination of migration parameters m 1 and m 2 , the thresholds responsible for the existence and stability of states E 1 , E 2 and E * have been computed in Table 2 .

For case 1, the existence condition (25) of the state E * is violated while both the states E 1 and E 2 exist and are given below: According to stability conditions (30) and (32) of the states E 1 and E 2 respectively, E 1 is found to be locally stable while E 2 gets unstable. A projection of phase plot in S 1 − I 1 − V 1 hyperplane has been drawn in Fig. 2 for the set of initial conditions Y 1 , Y 2 , Y 3 and Y 4 chosen in the neighborhood of E 2 . The solution trajectories are found to be converging to the state E 1 in hyperplane S 1 − I 1 − V 1 .

Again, for case 2 in Table 2 , the existence condition (25) of the state E * is violated while both the states E 1 and E 2 exist. Fig. 3 shows a projection of phase plot in S 2 − I 2 − I 12 hyperplane for the initial conditions Z 1 , Z 2 , Z 3 and Z 4 chosen in the neighborhood of E 1 . It is observed that the all the trajectories starting from the neighborhood of E 1 converge to the state E 2 .

For case 3, the states E 1 , E 2 and E * exist: The initial conditions C 1 and C 2 have been considered in the neighborhood of E 1 and the C 3 and C 4 have been considered in the neighborhood of E 2 . It is observed from a 3D phase plot I 1 − I 2 − I 12 drawn in Fig. 4 that all the solution trajectories converge to E * in hyperplane I 1 − I 2 − I 12 . This verifies the instability of the states E 1 (R 2 01 ξ > 1) and E 2 (R 2 10 η > 1). Further, this confirms the stability of E * . 

In disease modeling, sensitivity analysis is important to know the change in basic reproduction number (R 0 ) in response to the parameters it involves. The control measures can then be targeted on the basis of the sign and magnitude of the sensitivity indices with respect to various parameters. To compute the sensitivity of parameters of different scale, the normalized Table 3 Sensitivity indices of R 10 forward sensitivity index has been used [42] . The normalized forward sensitivity index of R 0 for the parameter p (say) is defined as

Consider basic reproduction numbers R 10 and R 01 in terms of dimensional parameters as

and R 2 01 = β 2 ωω 1 σ (μ + 2m 1 ) μμ 2 1 (μ + γ 2 )(μ + m 1 + m 2 ) The sensitivity indices with respect to the migration parameters (m 1 = m 2 ) for the two patches are obtained as follows:

;

Due to emigration from patch-1 Υ R 10 m 1 is found to be negative. The same is true for Υ R 01 m 2 . However, Υ R 10 m 2 and Υ R 01 m 1 are observed to be positive. This is due to immigration in respective patches.

From the above expressions, it can be observed that sensitivity indices for the migration parameters depend on host death rate μ also. For μ = 0.00004, the sensitivity indices with respect to migration parameters are given in Table 3 . It can be concluded that when rate of immigration is higher than emigration in a patch, the basic reproduction number of respective patch increases. Particularly, from given Table 3 , if m 1 increases by 10%, the basic reproduction number for patch-1 decreases by 1.238%. However, there will be increase of 3.664% in patch-2.

In this paper, a two patch model incorporating human migration is applied for 2003 dengue outbreak in the two states of Brazil namely, Ceara and Rio de Janeiro. The first case of DEN-2 serotype in Rio de Janeiro was reported in 1990 then spread over to other states of the country namely Ceara, Bahia, Rio Grande do Norte, Alagoas, So Paulo, Mato Grosso do Sul, Mato Grosso, with the highest incidence rate in northeast region [43] . During summer 2002 a large outbreak of dengue due to DEN-3 serotype occurred in Rio de Janeiro [44] . In 2003, the human migration rate in Ceara was quite high as compared to Rio de Janeiro [45] . The sequential infection by DEN-2 serotype led to severe secondary infection cases in the form of DHF in the state Ceara [46, 47] . Keeping these facts in mind, let us consider the state Rio de Janeiro as patch-1 with DEN-3 serotype prevalent at time t and the state Ceara as patch-2 where DEN-2 serotype is present. The estimation of R 0 has been performed for the two patches from the initial growth phase of the epidemics [48] . Let us consider that at the beginning of the epidemic, the cumulative number of cases, K (t), varies exponentially as K (t) ∝ K 1 exp(λt) (33) where, λ is the force of infection and K 1 is constant. Accordingly, the infected host and vector population of both the patches are assumed as

Again, the number of non-susceptible hosts and vectors can be assumed negligible,

Now using Eq. (34) for infected host and vector population in the model (1)- (11) gives,

Multiplying Eqs. (36) and (37) together for patch-1 and Eqs. (38) and (39) for patch-2 respectively give the following:

The R E 10 and R E 01 are estimated reproduction numbers for the patch-1 and patch-2 respectively. Now, from the Eq. (33), it can be easily seen that, at the beginning of the epidemic number of new cases in a month in the patch-1 and patch-2 (say J 1 and J 2 respectively) would be proportional to the cumulative number of cases, i.e. J 1 ∼ λ 1 I 10 K 1 (t) and J 2 ∼ λ 2 I 20 K 2 (t) . Again by plotting the number of new cases per month against the cumulative number of cases K (t) for the two patches separately, the force of infection which would be the slope of respective curves, can be obtained [49] .

Using the monthly data of the two states of Brazil namely, Rio de Janeiro (say patch-1) and Ceara (say patch-2) of 2003 from "WHO Dengue NET", the commutative number of cases have been plotted against new number of cases in Figs. 5 and 6 respectively. The force of infection of patch-1 and patch-2 have been given in Table 4 by computing the slope of respective curves. The time series for the data of respective states has been given in Tables 5  and 6 in "Appendix". The host parameters (host death rate and recovery rate) have been taken from literature [40] . The estimation of mosquito mortality rates are described in next subsection. 

The mortality rate of mosquitoes is temperature dependent. By using the average temperature records of nine months of the states Rio de Janeiro and Ceara, the estimation of mortality rate of the mosquito Aedes aegypti has been explored by the enzyme experiment [50, 51] .

The following formula is obtained to calculate the mortality rate of Aedes aegypti mosquito: μ m (T ) = 0.8692 − 0.159T + 0.01116T 2 − 3.408 × 10 −4 T 3 + 3.809 × 10 −6 T 4 (42) where, T is the temperature and μ m (T ) is the mortality rate of mosquito. By using the monthly temperature records of the states Ceara [52] and Rio de Janeiro [53] and substituting in Eq. (42) , the range of mortality rate for the two respective states has been calculated. The mortality rate of Aedes aegypti mosquito is found to be in the range of 0.02987767-0.0363767 for the state Rio de Janeiro. However, for the state Ceara, it is coming in the range of 0.02714614-0.03403341. Using the average value of the mortality rates, the estimated value of basic reproduction number of the respective states are given in Table 4 .

It is concluded from the Table 4 that the basic reproduction number for the Ceara state of Brazil is greater than that of the Rio de Janeiro state. This verifies the fact that severe epidemic was occurred in Ceara in 2003. It supports that human migration is one of the reasons for the severe outbreak of DHF.

In this paper, a non-linear two patch dynamic model has been proposed to study the dynamics of dengue transmission. It is assumed that the two distinct serotypes are predominant in respective patches. The main emphasis is given on inter-patch migration of human pop-ulation and its impact on spatial spread of primary as well as secondary infection. The existence conditions for equilibrium states are obtained in terms of R 0 . The disease-free state is found to be globally stable provided the basic reproduction numbers (R 10 , R 01 ) of both the patches are less than one. It is found from analysis that for the equilibrium state E 1 (E 2 ), where disease persists in patch-1 (patch-2), it will persist there irrespective of migration if the basic reproduction number of patch-2 (patch-1) is less than 1. Numerical implementations have been carried out for the relevant data from the literature to explore the stability of endemic state. For the set of data where E * exists, it is found that starting from the initial conditions in the neighborhood of the states E 1 and E 2 , the solution trajectories converge to the state E * . This shows the stability of the state E * .

The magnitude of migration parameters and their mutual relationships are very crucial to predict the stability of various states. When equal rates of migration are considered in both the patches, the status of disease will not be affected due to migration. In other words, the disease-free state will not become endemic or vice-versa in any patch due to migration. However, the level of primary/secondary infection changes without affecting the disease-free/endemic status. On the other hand, when rate of migration is different in the patches, the persistence/extinction of infection will depend on the magnitude of the basic reproduction number of patches as well as the migration parameters.

Sensitivity indices of the basic reproduction number for the migration parameters are computed to analyze the effect of emigration and immigration in respective patches. Accordingly, it is found that emigration from a patch decreases the basic reproduction number of respective patch which consequently increases the chances to establish the infectionfree state. However, higher rates of immigration in a patch may increase the infection level.

The case study for the two states of Brazil namely, Rio de Janeiro and Ceara is carried out to validate the model for the sequential infection by DEN-2 and DEN-3 serotypes occurred during 2003. It has been found from literature that Ceara had more immigration as compared to Rio de Janeiro in 2003 because of this Ceara had more severe dengue epidemic in 2003. By estimation of basic reproduction numbers for these two states, the basic reproduction number for Ceara state is found to be higher than the Rio de Janeiro state. This is in line with the fact that severe epidemic was occurred in Ceara in 2003.

The details of the Jacobian matrices for computing the basic reproduction number and for the system (12)- (22) are as follows:

The Jacobian matrices of the system (12)- (22) for the new infections (F) and transfer from one compartment to another (Y ) are given below: 

Resurgent vector-borne, diseases as a global health problem

Emerging flaviviruses: the spread and resurgence of Japanese encephalitis, West Nile and dengue viruses

Dengue viral infections

Dengue and Dengue Hemorrhagic Fever

Emerging and resurging vector-borne diseases

Clinical and laboratory presentations of dengue patients with different serotypes

Role of macrophage in the disease dynamics of cutaneous Leishmaniasis: a delay induced mathematical study

A mathematical approach to control cutaneous leishmaniasis through insecticide spraying

Impulsive application to design effective therapies against cutaneous Leishmaniasis under mathematical perceptive

A mathematical model of dengue transmission with memory

The effects of awareness and vector control on two strains dengue dynamics

Modelling the use of wolbachia to control dengue fever transmission

The role of seasonality and import in a minimalistic multi-strain dengue model capturing differences between primary and secondary infections: complex dynamics and its implications for data analysis

Dengue transmission during the FIFA World Cup in Brazil

How much complexity is needed to describe the fluctuations observed in dengue hemorrhagic fever incidence data?

Dengue diversity across spatial and temporal scales: local structure and the effect of host population size

Socially structured human movement shapes dengue transmission despite the diffusive effect of mosquito dispersal

The role of human movement in the transmission of vector-borne pathogens

Dengue virus in sub-tropical Northern and Central Viet Nam: population immunity and climate shape patterns of viral invasion and maintenance

Epidemiological risk factors associated with high global frequency of inapparent dengue virus infections

DHF epidemics in Cuba, 1981 and 1997: some interesting observations

Enhanced severity of secondary dengue-2 infections: death rates in 1981 and 1997 Cuban outbreaks

Fatal dengue hemorrhagic fever in Cuba, 1997

A multi-city epidemic model

Disease spread in metapopulations, nonlinear dynamics and evolution equations

A multipatch malaria model with logistic growth populations

Impact of travel between patches for spatial spread of disease

Modeling the spread of influenza among cities

An epidemic model in a patchy environment

The effect of global travel on the spread of SARS

Assessing the role of spatial heterogeneity and human movement in malaria dynamics and control

A dengue model incorporating saturation incidence and human migration

Serotype-specific differences in the risk of dengue hemorrhagic fever: an analysis of data collected in Bangkok

Serotype-specific dengue virus circulation and dengue disease in Bangkok, Thailand from 1973 to 1999

The Theory of the Chemostat

Notes on R 0

The stability of dynamical systems

Competitive exclusion in a vector-host model for the dengue fever

Modelling the dynamics of dengue real epidemics

Analysis of a dengue disease transmission model

Sensitivity Analysis in a Dengue Epidemiological Model

Epidemiological Scenario of Dengue in Brazil

Dengue Virus Type 3, Brazil

Climate Change Impacts on Agriculture and Internal Migration in Brazil

Change in age pattern of persons with dengue, northeastern Brazil

Clinical and epidemiological characterization of dengue hemorrhagic fever cases in northeastern Brazil

Estimation of R 0 from the initial phase of an outbreak of a vector-borne infection

Early determination of the reproductive number for vector-borne diseases: the case of dengue in Brazil

Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue

The seasonal reproduction number of dengue fever: impacts of climate on transmission

The next generation matrix FY −1 at the disease-free equilibrium point E 0 is evaluated as:The dominant eigenvalue of the above matrix is the basic reproduction number, R 0 (say). It is given as,The general matrix J [E] of the system (12)-(22) about any equilibrium state (E) is computed as:j 4,4 = −p 3 − 1 − p 4 V 2 ; j 5,5 = −1 − p 6 ; j 6,6 = −p 3 − 1; j 6,11 = −ρp 2 p 7 R 1 ; j 7,7 = −1 − p 6 − p 8 ; j 7,11 = ρp 2 p 7 R 1 ; j 8,8 = −1 − p 9 I 1 ; j 9,9 = −1 − p 9 I 2The details of number of dengue cases of Rio de Janeiro and Ceara are as follows: