# Modeling and Analysis of Epidemic Diffusion with Population Migration

Authors: Ming Liu; Yihong Xiao
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

An improved Susceptible-Infected-Susceptible (SIS) epidemic diffusion model with population migration between two cities is modeled. Global stability conditions for both the disease-free equilibrium and the endemic equilibrium are analyzed and proved. The main contribution of this paper is  reflected in epidemic modeling and analysis which considers unequal migration rates, and only susceptible individuals can migrate between the two cities. Numerical simulation shows when the epidemic diffusion system is stable, number of infected individuals in one city can reach zero, while the number of infected individuals in the other city is still positive. On the other hand, decreasing population migration in only one city seems not as effective as improving the recovery rate for controlling the epidemic diffusion.


---

## Body

## 1. Introduction

As mentioned in Rachaniotis et al. [1], a serious epidemic is a problem that tests the ability of a nation to effectively protect its population, to reduce human loss, and to rapidly recover. Sometime such a problem may acquire worldwide dimensions. For example, during the period from November 2002 to August 2003, 8422 people in 29 countries were infected with SARS; 916 of them were dead at last for the effective medical resources appeared late. Other diseases, such as HIV and H1N1, can also cause significant numbers of direct infectious disease deaths. Epidemic diffusion is a typical complex dynamic system problem in Gao et al. [2], for we do not know what kind of epidemic outbreaks, when it outbreaks, and how it diffuses. Generally, after an epidemic outbreaks, public officials are faced with many critical and complex issues, the most important of which is to make certain how the epidemic diffuses so that the rescue operation efficiency is maximized.Traditionally, analytical works on epidemic diffusion are concentrated on the compartmental epidemic models of ordinary differential equations (Mishra and Saini [3]; Sun and Hsieh [4]; Li et al. [5]; Zhang et al. [6]; Zhang and Ma [7]). In these models, the total population is divided into several independence classes, and each class of individuals is closed into a compartment. The sizes of the compartments are large enough, and the mixing of members is homogeneous. In other words, the models based on the differential equations are always under the assumption of both homogeneous infectivity and homogeneous connectivity of each individual. Recently, Kim et al. [8] described the transmission of avian influenza among birds and humans. The behavior of positive solutions to a reaction diffusion system with homogeneous Neumann boundary conditions was investigated in their work. Liu and Zhang [9] presented an SEIRS epidemic model on the scale-free networks, where the active contact number of each vertex was assumed to be either constant or proportional to its degree in their model. Samsuzzoha et al. [10] used a diffusive epidemic model to describe the transmission of influenza. The equations were solved numerically by using the splitting method under different initial distribution of population density. Further, Samsuzzoha et al. [11] presented a vaccinated diffusive compartmental epidemic model to explore the impact of vaccination as well as diffusion on the transmission dynamics of influenza. A very recent research by Shi and Dong [12] formulates and discusses models for the spread of infectious diseases with variable population sizes and vaccinations on the susceptible individuals.The other stream of related research to our work is on the epidemic diffusion with population migration. For instance, Hethcote [13] proposed that deterministic communicable disease models were initial value problems for a system of ordinary differential equations, and thus he considered the asymptotic stability for the equilibrium points for models involving temporary immunity, disease-related fatalities, carriers, migration, dissimilar interacting groups, and transmission by vectors. In his work, both susceptible individuals and infected individuals in each population could migrate (only equal rates were considered), which led to different equilibriums. Another model that considers two interacting populations undergoing SIS dynamics was presented in Kribs-Zaleta and Velasco-Hernández [14]. The authors considered that the two groups may have different values for model parameters especially those dealing with vaccination. Liebovitch and Schwartz [15] proposed that classical disease models always use a mass action term as the interaction between infected and susceptible people in separate patches, and they derived the equations when this interaction is a migration of people between patches. Sani and Kroese [16] formulated various mathematical control problems for HIV spread in mobile heterosexual populations. They applied the cross-entropy method to solve these highly multimodal and nonlinear optimization problems, and demonstrated the effectiveness of the method via a range of experiments and illustrated how the form of the optimal control function depends on the mathematical model used for the HIV spread. Yang et al. [17] considered SIR and SIS epidemic models with bilinear incidence and migration between two patches, where infected individuals cannot migrate from one patch to another due to medical screening. They found the thresholds classifying the global dynamics of the models in terms of the model parameters, and they obtained the global asymptotical stability of the disease free and the disease-endemic equilibrium. Wolkewitz and Schumacher [18] pointed out that the main limitation of the compartmental models is that several parameters are based on uncertain expert guesses (default values) and are not estimated from the study data. Lee et al. [19] extended the SEIR model to incorporate population migration between cities and investigated the effectiveness of travel restrictions as a control against the spread of influenza.It is worth mentioning that majority of the existing studies relies on different kinds of differential equations. For instance, first-order partial differential equations are used to integrate the age structures; second-order partial differential equations are suitable when a diffusion term exists; integral differential equations or differential equations are often used when time delay or delay factors are considered. These are the typical methods to model and analyze the complex dynamic system of epidemic diffusion. As a continued work of our previous paper in [20], this paper presents an SIS epidemic model with population migration between two cities. We consider unequal migration rates for these two populations, and only susceptible individuals can migrate, which is different from the whole existing works. The remainder of the paper is organized as follows. Section 2 proposes the SIS epidemic diffusion model, which considers susceptible individuals migration between two cities. In Section 3, the dynamic characteristics of the proposed SIS model are analyzed. Numerical simulation and key parameters analysis are presented in Section 4. Finally, Section 5 discusses the conclusions and suggests for the future research.


## 2. Epidemic Diffusion Model (SIS)

Since the compartment model of epidemic diffusion is a mature theory, herein we omit the verbose introduction of the framework process. In this paper, we divide people in epidemic areas into two groups: susceptible individuals(S) and infected individuals (I). The transfer diagram of individuals in the epidemic areas can be illustrated as Figure 1.Figure 1
The transfer diagram of SIS model with population migration.To smooth the formulation progress of the SIS epidemic diffusion model in the following subsections, some assumptions and parameters are specified as follows.(1)
The susceptible individuals and the infected individuals in cityi at time t are denoted as Si(t) and Ii(t), respectively. Thus, the total individuals in city i are Ni(t)=Si(t)+Ii(t),   i=1,2.(2)
b and d are the natural birth rate and the natural death rate, respectively. γ is the recovery rate. β is the propagation coefficient. To facilitate the process in the following sections, we assume that b=d. Moreover, disease-related death rate is not considered in this work.(3)
Only the susceptible individuals can migrate in this paper.ai represents the migrating-out rate of susceptible individuals in city i(ai>0 for i=1,2 and a1≠a2).(4)
Using the notationN to represent the total number of the population in these two cities, N=N1+N2. Note that N is a constant.Hence, the ordinary differential equations for the SIS epidemic diffusion model can be formulated as(1)dS1dt=bN1-a1S1+a2S2-dS1-βS1I1+γI1,dI1dt=βS1I1-γI1-dI1,dS2dt=bN2+a1S1-a2S2-dS2-βS2I2+γI2,dI2dt=βS2I2-γI2-dI2.ODE (1) describes the following dynamics of epidemic diffusion among the population groups. (1) The change rate of the susceptible population in both city 1 and city 2 is determined by the entry population, the exiting population, and the losing population who actually gets exposed to the disease and thus is counted towards the class of infected population. The last one is in proportion to the propagation coefficient β and both of the current mass of the susceptible individuals and the current mass of the infected individuals. (2) The change rate of the infected population is determined by the difference between the entering population, those of the susceptible population who get sick, the exiting population, and the losing population. All parameters β,b,γ,a1, anda2 are positive, and initial conditions for the model are demonstrated as follows:
(2)I1(0)=i10≪N,I2(0)=i20≪N,S1(0)=s10,S2(0)=N-s10-i10-i20.


## 3. Model Analysis

### 3.1. Condition of the Epidemic Diffusion

As shown above,I1(0)=i10≪N, I2(0)=i20≪N, S1(0)=s10, and S2(0)=N-s10-i10-i20 are initial conditions for the proposed model, which symbolize the initial number of susceptible and infected individuals. Then, it is easy to obtain the initial condition for epidemic diffusion, which should satisfy the following premise:
(3)dI1dt|t=0>0ordI2dt|t=0>0.Taking it into (1), we can obtain the initial condition of the susceptible individuals in city 1 and city 2:
(4)s10>b+γβors20<N-i10-i20-b+γβ.Equation (4) shows that the spread of epidemic happens only when s10 and s20 meet the previous initial conditions.


### 3.2. Existence of the System Equilibrium Solution

Generally, it is difficult to obtain the analytic solution of (1). To analyze the epidemic diffusion, we consider the stable state of (1). Considering that b=d and expunging S2, (1) can be rewritten as
(5)dS1dt=-a1S1+a2(N-S1-I1-I2)-βS1I1+(b+γ)I1,dI1dt=βS1I1-(b+γ)I1,dI2dt=β(N-S1-I1-I2)I2-(b+γ)I2.LettingdI1/dt=0, we can get I1=0 or S1=(b+γ)/β. Similarly, letting dI2/dt=0, we can obtain I2=0 or S1+I1+I2=N-((b+γ)/β). With the partial derivatives dS1/dt=0, dI1/dt=0, and dI2/dt=0, we can obtain one equilibrium point of the SIS epidemic diffusion system intuitively when I1=0 and I2=0:
(6)P1=(S1,I1,I2)=(a2a1+a2N,0,0).From (6) we can see that both numbers of infected individuals in city 1 and city 2 are zero, which indicate that epidemic diffusion in these two cities does not happen, and all the individuals in these two cities are susceptible individuals at last. Herein, we call it disease-free equilibrium point.WhenI1=0 and S1+I1+I2=N-((b+γ)/β), we can obtain the second equilibrium point of the SIS epidemic diffusion system:
(7)P2=(S1,I1,I2)=(a2a1·b+γβ,0,N-a1+a2a1·b+γβ).From (7), when the SIS epidemic diffusion system is stable, the number of infected individuals in city 1 is zero, and some infected individuals in city 2 still exist. In this condition, we call it the endemic equilibrium point.Likewise, whenS1=(b+γ)/β and I2=0, we can obtain the third equilibrium point of the SIS epidemic diffusion system:
(8)P3=(S1,I1,I2)=(b+γβ,N-a1+a2a2·b+γβ,0).In line with the above work, when the SIS epidemic diffusion system is stable, the number of infected individuals in city 2 is zero, and some infected individuals in city 1 still exist. So it is called endemic equilibrium point as well.It is worth mentioning that, whenS1=(b+γ)/β and S1+I1+I2=N-((b+γ)/β), there is dS1/dt=(a2-a1)·(b+γ)/β≠0 for that a1≠a2. That means, under the conditions of S1=(b+γ)/β and S1+I1+I2=N-((b+γ)/β), there is no solution for the simultaneous equations dS1/dt=0, dI1/dt=0, and dI2/dt=0.


### 3.3. Stability of the System Equilibrium Solution

Lemma 1.
Disease-free equilibrium pointP1 in the SIS epidemic diffusion system is stable only when β<min{(a1+a2)(b+γ)/a1N,(a1+a2)(b+γ)/a2N}.Proof.
LettingP=dS1/dt, Q=dI1/dt, and R=dI2/dt, the Jacobi matrix of (5) can be obtained as follows: (9)J=(∂P∂S1∂P∂I1∂P∂I2∂Q∂S1∂Q∂I1∂Q∂I2∂R∂S1∂R∂I1∂R∂I2)=(-a1-a2-βI1b+γ-a2-βS1-a2βI1βS1-(b+γ)0-βI2-βI2β(N-S1-I1-2I2)-(b+γ)).
ForP1=(S1,I1,I2)=(a2/(a1+a2)N,0,0), the Jacobi matrix J can be rewritten as follows:(10)
J
P
1
=
(
-
a
1
-
a
2
b
+
γ
-
a
2
-
a
2
N
a
1
+
a
2
β
-
a
2
0
a
2
N
a
1
+
a
2
β
-
(
b
+
γ
)
0
0
0
a
1
N
a
1
+
a
2
β
-
(
b
+
γ
)
)According to the Jacobi matrixJP1, it is easy to obtain the secular equation of (5):
(11)(λ+a1+a2)(λ+b+γ-a2Na1+a2β)×(λ+b+γ-a1Na1+a2β)=0.The three latent roots of this secular equation are-a1-a2, a2N/(a1+a2)β-b-γ, and a1N/(a1+a2)β-b-γ. Based on Routh-Hurwitz stability criterion, only when β<(a1+a2)(b+γ)/a1N and β<(a1+a2)(b+γ)/a2N, three latent roots of the secular equation would have negative real parts, simultaneously, and then P1=(S1,I1,I2)=(a2/(a1+a2)N,0,0) is the stable solution of the differential equations.Lemma 2.
Endemic equilibrium pointP2 in the SIS epidemic diffusion system is stable only when a2<a1 and β>(a1+a2)(b+γ)/a1N.Proof.
As far as we are concerned, if the endemic equilibrium pointP2=(S1,I1,I2)=(a2/a1·(b+γ)/β,0,N-(a1+a2)/a1·(b+γ)/β) exists, it should satisfy condition I2>0 firstly. Namely, the propagation coefficient β should satisfy β>(a1+a2)(b+γ)/a1N. Then, similar to Lemma 1, we can obtain the Jacobi matrix for P2 as follows:(12)
J
P
2
=
(
-
a
1
-
a
2
b
+
γ
-
a
2
-
a
2
a
1
(
b
+
γ
)
-
a
2
0
a
2
a
1
(
b
+
γ
)
-
(
b
+
γ
)
0
a
1
+
a
2
a
1
(
b
+
γ
)
-
β
N
a
1
+
a
2
a
1
(
b
+
γ
)
-
β
N
a
1
+
a
2
a
1
(
b
+
γ
)
-
β
N
)
.According to the Jacobi matrixJP2, we can obtain the secular equation of (5) again:
(13)[λ+b+γ-a2a1(b+γ)](λ2+A1λ+A0)=0,
where A0=a1[βN-((a1+a2)/a1)(b+γ)] and A1=a1+a2+βN-((a1+a2)/a1)(b+γ). Obviously, one of the latent roots of (13) is λ1*=((a2-a1)/a1)(b+γ). Only when a2<a1, the latent root λ1* is negative. On the other hand, when β>(a1+a2)(b+γ)/a1N, there is A0>0 and A1>0. Based on Routh-Hurwitz stability criterion, the other two latent roots of (13) will be with negative real parts. Therefore, P2 is the stable solution of the simultaneous differential equations only when a2<a1 and β>(a1+a2)(b+γ)/a1N.Lemma 3.
Endemic equilibrium pointP3 in the SIS epidemic diffusion system is stable only when a1<a2 and β>((a1+a2)(b+γ))/a2N.Proof.
Similar to Lemma2, if the endemic equilibrium point P3=(S1,I1,I2)=((b+γ)/β,N-(a1+a2)/a2·(b+γ)/β,0) exists, it should satisfy condition I1>0. That is, the propagation coefficient β should satisfy β>(a1+a2)(b+γ)/a2N. Then, we can obtain the Jacobi matrix for P3 as follows:(14)
J
P
3
=
(
-
a
1
-
a
2
-
β
N
+
a
1
+
a
2
a
2
(
b
+
γ
)
-
a
2
-
a
2
β
N
-
a
1
+
a
2
a
2
(
b
+
γ
)
0
0
0
0
a
1
-
a
2
a
2
(
b
+
γ
)
)
.Again, according to the Jacobi matrixJP3, we can obtain the secular equation of (5):
(15)[λ-a1-a2a2(b+γ)](λ2+B1λ+B0)=0,
where B0=a2[βN-((a1+a2)/a2)(b+γ)] and B1=a1+a2+βN-((a1+a2)/a2)(b+γ). Obviously, one of the latent roots of (15) is λ1*=((a1-a2)/a2)(b+γ). Only when a1<a2, the latent root λ1* is negative. On the other hand, when β>(a1+a2)(b+γ)/a2N, there is B0>0 and B1>0. Based on Routh-Hurwitz stability criterion, the other two latent roots of (15) will be with negative real parts. Therefore, P3 is the stable solution of the simultaneous differential equations only when a1<a2 and β>(a1+a2)(b+γ)/a2N.Remark 4.
From Lemmas1, 2, and 3, we can draw a conclusion that the diffusion threshold of the SIS epidemic diffusion model relies on the migrating-out coefficients of susceptible individuals of the two cities ai(i=1,2); it also depends on the three key parameters: the total individuals of the two cities N, the birth rate b, and the recovery rate γ.


## 3.1. Condition of the Epidemic Diffusion

As shown above,I1(0)=i10≪N, I2(0)=i20≪N, S1(0)=s10, and S2(0)=N-s10-i10-i20 are initial conditions for the proposed model, which symbolize the initial number of susceptible and infected individuals. Then, it is easy to obtain the initial condition for epidemic diffusion, which should satisfy the following premise:
(3)dI1dt|t=0>0ordI2dt|t=0>0.Taking it into (1), we can obtain the initial condition of the susceptible individuals in city 1 and city 2:
(4)s10>b+γβors20<N-i10-i20-b+γβ.Equation (4) shows that the spread of epidemic happens only when s10 and s20 meet the previous initial conditions.


## 3.2. Existence of the System Equilibrium Solution

Generally, it is difficult to obtain the analytic solution of (1). To analyze the epidemic diffusion, we consider the stable state of (1). Considering that b=d and expunging S2, (1) can be rewritten as
(5)dS1dt=-a1S1+a2(N-S1-I1-I2)-βS1I1+(b+γ)I1,dI1dt=βS1I1-(b+γ)I1,dI2dt=β(N-S1-I1-I2)I2-(b+γ)I2.LettingdI1/dt=0, we can get I1=0 or S1=(b+γ)/β. Similarly, letting dI2/dt=0, we can obtain I2=0 or S1+I1+I2=N-((b+γ)/β). With the partial derivatives dS1/dt=0, dI1/dt=0, and dI2/dt=0, we can obtain one equilibrium point of the SIS epidemic diffusion system intuitively when I1=0 and I2=0:
(6)P1=(S1,I1,I2)=(a2a1+a2N,0,0).From (6) we can see that both numbers of infected individuals in city 1 and city 2 are zero, which indicate that epidemic diffusion in these two cities does not happen, and all the individuals in these two cities are susceptible individuals at last. Herein, we call it disease-free equilibrium point.WhenI1=0 and S1+I1+I2=N-((b+γ)/β), we can obtain the second equilibrium point of the SIS epidemic diffusion system:
(7)P2=(S1,I1,I2)=(a2a1·b+γβ,0,N-a1+a2a1·b+γβ).From (7), when the SIS epidemic diffusion system is stable, the number of infected individuals in city 1 is zero, and some infected individuals in city 2 still exist. In this condition, we call it the endemic equilibrium point.Likewise, whenS1=(b+γ)/β and I2=0, we can obtain the third equilibrium point of the SIS epidemic diffusion system:
(8)P3=(S1,I1,I2)=(b+γβ,N-a1+a2a2·b+γβ,0).In line with the above work, when the SIS epidemic diffusion system is stable, the number of infected individuals in city 2 is zero, and some infected individuals in city 1 still exist. So it is called endemic equilibrium point as well.It is worth mentioning that, whenS1=(b+γ)/β and S1+I1+I2=N-((b+γ)/β), there is dS1/dt=(a2-a1)·(b+γ)/β≠0 for that a1≠a2. That means, under the conditions of S1=(b+γ)/β and S1+I1+I2=N-((b+γ)/β), there is no solution for the simultaneous equations dS1/dt=0, dI1/dt=0, and dI2/dt=0.


## 3.3. Stability of the System Equilibrium Solution

Lemma 1.
Disease-free equilibrium pointP1 in the SIS epidemic diffusion system is stable only when β<min{(a1+a2)(b+γ)/a1N,(a1+a2)(b+γ)/a2N}.Proof.
LettingP=dS1/dt, Q=dI1/dt, and R=dI2/dt, the Jacobi matrix of (5) can be obtained as follows: (9)J=(∂P∂S1∂P∂I1∂P∂I2∂Q∂S1∂Q∂I1∂Q∂I2∂R∂S1∂R∂I1∂R∂I2)=(-a1-a2-βI1b+γ-a2-βS1-a2βI1βS1-(b+γ)0-βI2-βI2β(N-S1-I1-2I2)-(b+γ)).
ForP1=(S1,I1,I2)=(a2/(a1+a2)N,0,0), the Jacobi matrix J can be rewritten as follows:(10)
J
P
1
=
(
-
a
1
-
a
2
b
+
γ
-
a
2
-
a
2
N
a
1
+
a
2
β
-
a
2
0
a
2
N
a
1
+
a
2
β
-
(
b
+
γ
)
0
0
0
a
1
N
a
1
+
a
2
β
-
(
b
+
γ
)
)According to the Jacobi matrixJP1, it is easy to obtain the secular equation of (5):
(11)(λ+a1+a2)(λ+b+γ-a2Na1+a2β)×(λ+b+γ-a1Na1+a2β)=0.The three latent roots of this secular equation are-a1-a2, a2N/(a1+a2)β-b-γ, and a1N/(a1+a2)β-b-γ. Based on Routh-Hurwitz stability criterion, only when β<(a1+a2)(b+γ)/a1N and β<(a1+a2)(b+γ)/a2N, three latent roots of the secular equation would have negative real parts, simultaneously, and then P1=(S1,I1,I2)=(a2/(a1+a2)N,0,0) is the stable solution of the differential equations.Lemma 2.
Endemic equilibrium pointP2 in the SIS epidemic diffusion system is stable only when a2<a1 and β>(a1+a2)(b+γ)/a1N.Proof.
As far as we are concerned, if the endemic equilibrium pointP2=(S1,I1,I2)=(a2/a1·(b+γ)/β,0,N-(a1+a2)/a1·(b+γ)/β) exists, it should satisfy condition I2>0 firstly. Namely, the propagation coefficient β should satisfy β>(a1+a2)(b+γ)/a1N. Then, similar to Lemma 1, we can obtain the Jacobi matrix for P2 as follows:(12)
J
P
2
=
(
-
a
1
-
a
2
b
+
γ
-
a
2
-
a
2
a
1
(
b
+
γ
)
-
a
2
0
a
2
a
1
(
b
+
γ
)
-
(
b
+
γ
)
0
a
1
+
a
2
a
1
(
b
+
γ
)
-
β
N
a
1
+
a
2
a
1
(
b
+
γ
)
-
β
N
a
1
+
a
2
a
1
(
b
+
γ
)
-
β
N
)
.According to the Jacobi matrixJP2, we can obtain the secular equation of (5) again:
(13)[λ+b+γ-a2a1(b+γ)](λ2+A1λ+A0)=0,
where A0=a1[βN-((a1+a2)/a1)(b+γ)] and A1=a1+a2+βN-((a1+a2)/a1)(b+γ). Obviously, one of the latent roots of (13) is λ1*=((a2-a1)/a1)(b+γ). Only when a2<a1, the latent root λ1* is negative. On the other hand, when β>(a1+a2)(b+γ)/a1N, there is A0>0 and A1>0. Based on Routh-Hurwitz stability criterion, the other two latent roots of (13) will be with negative real parts. Therefore, P2 is the stable solution of the simultaneous differential equations only when a2<a1 and β>(a1+a2)(b+γ)/a1N.Lemma 3.
Endemic equilibrium pointP3 in the SIS epidemic diffusion system is stable only when a1<a2 and β>((a1+a2)(b+γ))/a2N.Proof.
Similar to Lemma2, if the endemic equilibrium point P3=(S1,I1,I2)=((b+γ)/β,N-(a1+a2)/a2·(b+γ)/β,0) exists, it should satisfy condition I1>0. That is, the propagation coefficient β should satisfy β>(a1+a2)(b+γ)/a2N. Then, we can obtain the Jacobi matrix for P3 as follows:(14)
J
P
3
=
(
-
a
1
-
a
2
-
β
N
+
a
1
+
a
2
a
2
(
b
+
γ
)
-
a
2
-
a
2
β
N
-
a
1
+
a
2
a
2
(
b
+
γ
)
0
0
0
0
a
1
-
a
2
a
2
(
b
+
γ
)
)
.Again, according to the Jacobi matrixJP3, we can obtain the secular equation of (5):
(15)[λ-a1-a2a2(b+γ)](λ2+B1λ+B0)=0,
where B0=a2[βN-((a1+a2)/a2)(b+γ)] and B1=a1+a2+βN-((a1+a2)/a2)(b+γ). Obviously, one of the latent roots of (15) is λ1*=((a1-a2)/a2)(b+γ). Only when a1<a2, the latent root λ1* is negative. On the other hand, when β>(a1+a2)(b+γ)/a2N, there is B0>0 and B1>0. Based on Routh-Hurwitz stability criterion, the other two latent roots of (15) will be with negative real parts. Therefore, P3 is the stable solution of the simultaneous differential equations only when a1<a2 and β>(a1+a2)(b+γ)/a2N.Remark 4.
From Lemmas1, 2, and 3, we can draw a conclusion that the diffusion threshold of the SIS epidemic diffusion model relies on the migrating-out coefficients of susceptible individuals of the two cities ai(i=1,2); it also depends on the three key parameters: the total individuals of the two cities N, the birth rate b, and the recovery rate γ.


## 4. Numerical Simulation and Key Parameters Discussion

In this section, we take a numerical simulation to test how well the proposed model may be applied in practice. The initial values of parameters in the proposed epidemic diffusion model are acquired from Liu and Zhao [21] and by interviews with public health care administrative personnel, which are given as follows: β=8×10-6, b=2×10-4, γ=0.4, a1=0.02, a2=0.01,N=105, S1(0)=0.7×105, I1(0)=600, and I2(0)=400. We use MATLAB 7.0 mathematical solver together with Runge-Kutta method to simulate the epidemic model. The test is performed on an Intel (R) Core (TM) i3 CPU 2.4 GHz with 2 GB RAM under Microsoft Windows XP. Figure 2 is the evolution trajectories of the epidemic model. The curves, respectively, represent the different groups of people over time in these two cities.Figure 2
Evolution trajectories of the SIS epidemic model.From Figure2, one can see that the evolution trajectories of the SIS epidemic model with population migration between two cities are complicated. The number of susceptible individuals in city 1 decreases gradually with time increasing, while the number of susceptible individuals in city 2 increases at first and then decreases. On the other hand, the number of infected individuals in city 1 increases at first and then decreases, while the number of infected individuals in city 2 increases gradually with time increasing. However, all the susceptible individuals and infected individuals in city 1 and city 2 tend to the fixed values when time is long enough (t>300). Meanwhile, the limit value of the SIS epidemic diffusion model with population migration between two cities is Q1=(S1,S2,I1,I2)=(2.5013×104,5.0025×104,0,2.4961×104).In line with the initial values we defined previously, we havea2<a1 and β>(a1+a2)(b+γ)/a1N. According to Lemma 2 in Section 3, one can get that the number of susceptible and infected individuals will be converged at Q2=(S1*,S2*,I1*,I2*)=(2.50125×104,5.00245×104,0,2.4963×104). One can see that Q1 is very close to Q2, which is not a surprise, as it is consistent with the analytical conclusion in the last section. Once an epidemic outbreak occurs, we are more concerned with the change regularity of the infected individuals in practice. Therefore, in the following subsections, we will discuss the relationship between the key parameters and the number of the infected individuals.Figure3 demonstrates the change of the number of the infected individuals in both cities with different propagation coefficients. It is easy to know that the evolution trajectories of infected individuals in city 1 are different from those in city 2 for any certain propagation coefficient β. From Figure 3(a), we can see that, in the first some days, the larger β is, the faster spread of the epidemic in city 1 is. However, from Figure 3(b), we cannot get the similar conclusion. Number of infected individuals in city 2 is not in direct proportion to the propagation coefficient β. It is worth mentioning that, when the number of infected individuals in city 1 reaches zero, the number of infected individuals in the other city is still positive when the epidemic diffusion system is stable. This is consistent with Lemma 2 in Section 3. Similarly, if the initial conditions changed, we can also test and verify the other lemmas.Number of the infected individuals versus differentβ.
(a)
City 1
(b)
City 2Figure4 demonstrates the change of the number of the infected individualss in both cities with different recovery rate γ. As shown in Figure 4(a), when γ=0.2, the maximum number of infected individuals in city 1 is about 3.3×104. When γ=0.3, the number is about 1.5×104. When γ=0.4, the maximum number of infected individual in city 1 is less than 0.5×104. It informs us that the larger the recovery rate constant is, the smaller of the maximum number of infected individuals in city 1 is. Similar phenomenon can also be observed from Figure 4(b) in city 2. Such figure informs us that it is important to improve the recovery rate as much as possible when controlling an epidemic spread.Number of the infected individuals versus differentγ.
(a)
City 1
(b)
City 2Figure5 shows the change of the number of infected individuals in both cities with different migrating-out coefficients a1 and a2. According to Figure 5, one can observe that no matter in city 1 or in city 2, the evolution trajectories of the infected individuals may generate a serious change when the migrating-out coefficient of susceptible individuals changed. For example, in city 1, when a2<a1(a1=0.02, a2=0.01 and a1=0.02, a2=0.015), the number of infected individuals tends to be zero. However, when a1<a2(a1=0.02, a2=0.025 and a1=0.02, a2=0.03), the number of infected individuals tends to be a positive constant above 1×104. In other words, with the increment of migrating-out coefficient in city 2, the limit number of infected individuals in city 1 may become positive from zero. The larger the migrating-out coefficient in city 2 is, the larger the limit number of infected individuals in city 1 is. Opposite to city 1, when a2<a1(a1=0.02, a2=0.01 and a1=0.02, a2=0.015), the number of infected individuals in city 2 tends to be a positive constant above 1×104. When a1<a2(a1=0.02, a2=0.025 and a1=0.02, a2=0.03), the number of infected individuals is very small and tends to be zero at last. It informs us that with the increment of migrating-out coefficient in city 2, the limit number of infected individuals in city 2 may become zero from a positive value. The larger the migrating-out coefficient in city 2 is, the smaller the limit number of infected individuals in city 2 is. To summarize, decreasing the migration population in only one city is not as effective as improving the recovery rate for controlling the epidemic diffusion. However, we can find a trade-off between the migrating-out coefficients in these two cities and hence can control the infected individuals in both cities at last.Number of the infected individuals versus differenta1 and a2.
(a)
City 1
(b)
City 2


## 5. Conclusions

In this paper, we considered an SIS epidemic diffusion model with population migration between two cities, and proved that the model has one disease-free equilibrium point and two endemic equilibrium points. To summarize, our model differs from past researches in at least two aspects. Firstly, unequal migration rates between the two cities are considered; secondly, only the susceptible individuals can migrate between the two cities. Numerical simulation shows that number of infected individuals in one city can reach zero, while the number of infected individuals in the other city is still positive when the epidemic diffusion system is stable. One the other hand, decreasing population migration in only one city is not as effective as improving the recovery rate for controlling the epidemic diffusion.As a limitation of the model, disease-related death rate of infected individuals in both cities was not considered. Actually, some epidemics would cause death, and thus consideration of a disease-related death rate would make it more accurate. On the other hand, epidemic diffusion model with population migration among more than two cities should also be taken into account. A discussion of possible extension to more than two cities would be more realistic and the research would be more practical.



---

*Source: 583648-2013-08-29.xml*