# A Fractional SEIR Epidemic Model for Spatial and Temporal Spread of Measles in Metapopulations

Authors: Emile Franc Doungmo Goufo; Suares Clovis Oukouomi Noutchie; Stella Mugisha
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

Measles is a higher contagious disease that can spread in a community population depending
on the number of people (children) susceptible or infected and also depending on their movement in the
community. In this paper we present a fractional SEIR metapopulation system modeling the
spread of measles. We restrict ourselves to the dynamics between four distinct cities (patches).
We prove that the fractional metapopulation model is well posed (nonnegative solutions) and
we provide the condition for the stability of the disease-free equilibrium. Numerical simulations
show that infection will be proportional to the size of population in each city, but the disease
will die out. This is an expected result since it is well known for measles (Bartlett (1957)) that, in communities
which generate insufficient new hosts, the disease will die out.


---

## Body

## 1. Introduction and Important Facts

Mathematical modeling in epidemiology is concerned with describing the spread of disease and its effect on people. This itself encompasses a range of disciplines, from biology, mathematics, and engineering to sociology and philosophy, all of which are utilized to a better understanding and containing of the spread of infection. The models developed in epidemiology provide a wider understanding of the mechanisms that influence the spread of a disease and they suggest control strategies [1]. One of the early triumphs of mathematical epidemiology was a formulation of a model in [2] to predict the behavior of a disease. The total population in this model was assumed to be constant and divided into three compartments: suspended, infectious, and recovered. Over the years, more complex models have been derived. Among them, we count models based on dynamics with derivative of fractional order and metapopulation models.Fractional differential equations and their applications have been comprehensively studied in a large number of works [3–6]. Although numerous analyses have been done in modeling the dynamics of epidemiological diseases, most of them have been restricted to integer order (delay) differential equations. In recent years, it has turned out that many phenomena in different fields can be described very successfully by the models using fractional order differential equations. Adomian decomposition method has been employed to obtain solutions of a system of nonlinear fractional differential equations in number of articles and the method is extended in [5], where the authors used a complementary method, Frobenius decomposition method, to generalize the classical Darcy law. Epidemiology combined with fractional calculus was applied in [6] for an outbreak of dengue fever.Metapopulation models are defined as system of differential equations generated by discrete spatial models with continuous time. Indeed, they describe movement of individuals between discrete spatial patches that can be cities, towns, places, and so forth. Metapopulation models have been thoroughly analyzed in numeral articles [7, 8].Measles is a higher contagious viral disease caused by infection of Paramyxovirus, generally of the genusMorbillivirus. Its incubation period is located somewhere between 9 and 12 days and its infectivity period between 4 and 9 days. Measles is highly present in early childhood and its epidemics are commonly related to aggregation of children at schools or childcare centres. It is recommended to get vaccinated against it at around 18 months of age and have a booster at 4 to 5 years of age. The disease is particularly characterized by its low mortality and high morbidity. Measles will continue to circulate in a community with a higher number of susceptible hosts by birth of children. However, in communities which generate insufficient new hosts, measles will die out. This theory was introduced in 1957 by Bartlett [9], who brought out the critical population size for a community and referred it to be the minimum number supporting measles.Our work is motivated by an article by the online magazineOtago Daily Time [10], following the emergence of measles after a student brought the disease from Sydney where he participated in a hip hop competition. On 1 March 2014, the online magazineOtago Daily Time released a warning over measles outbreak in New Zealand after the disease was discovered on the student. After that, many cases of the disease in different cities all over the country were noticed. The metapopulation model, presented in Section 3, is then suitable to describe and analyze the spread of the disease.


## 2. Formulation of Fractional Temporal SEIR Measles Model

We formulate here the system modeling the fractional temporal spread of measles in a human population. In this model, a population supposed constant is divided into different classes, disjoint and based on their disease status. At timet, S
=
S
(
t
) is the fraction of population representing individuals susceptible to measles, E
=
E
(
t
) is the fraction of population representing individuals exposed to measles, I
=
I
(
t
) is the fraction of population representing individuals infectious with measles, and R
=
R
(
t
) is the fraction of population representing individuals that recovered from measles. We assume that all recruitment is done by birth into the class of susceptible and occurs at constant birth rate b. The rate constant for nondisease related death is μ; thus 1
/
μ is the average lifetime. We use the standard mass balance incidence expressions β
(
t
)
S
I to indicate successful transmission of measles due to effective contacts dynamics in the population by infectious individuals. Once infected, a fraction of exposed people becomes infectious with a constant rate σ, so that 1
/
σ is the average incubation period. Some infectious individuals will recover after a treatment or a certain period of time at a rate constant ζ, making 1
/
ζ the average infectious period. The transfer diagram for this model is described by Figure 1 and formulated by the following differential equations:
(1)
D
t
γ
S
=
b
-
(
β
(
t
)
I
+
μ
)
S
,
D
t
γ
E
=
β
(
t
)
S
I
-
(
σ
+
μ
)
E
,
D
t
γ
I
=
σ
E
-
(
ζ
+
μ
)
I
,
D
t
γ
R
=
ζ
I
-
μ
R
,
where
(2)
D
t
γ
u
(
t
)
=
∂
γ
∂
t
γ
u
(
t
)
=
1
Γ
(
1
-
γ
)
∫
0
t
(
t
-
r
)
-
γ
∂
∂
r
u
(
r
)
d
r
,
with 0
≤
γ
<
1 being the fractional derivative of the function u
(
t
) in the sense of Caputo [3], with Γ being the Gama function.Figure 1
Transfer diagram for temporal dynamics transmission of measles.


## 3. Metapopulation Model

In this section we formulate the metapopulation dynamics [7] for four patches representing four cities in New Zealand: Auckland, Bay of Plenty, Wellington, and Dunedin. Those cities are chosen since, on 1 March 2014, there were 43 cases of individuals found infected with measles in Auckland (among which 20 students were of the same school), 15 cases in the city of Bay of Plenty, and 3 cases in Wellington. There was no case found in Dunedin, a city in the southern region of Otago in South Island of New Zealand. The aim is to deal with prevention against measles to reach the region of Otago. For reasons of commodity, we denote the cities Auckland, Bay of Plenty, Wellington, and Dunedin by the letters A, B, W, and D, respectively. Let us consider the set P
=
{
A
;
B
;
W
;
D
} of the four patches. We call m
x
y
c the rate of travel from city x to city y in compartment c with c
=
S
,
E
,
I, or R; that is, it represents the transfer rate of individuals in the compartment c of city x moving to the same compartment c in city y. It is clear that m
x
x
c
=
0, for all x
∈
P and c
∈
C
=
{
S
,
E
,
I
,
R
}. Then the travel model is given in Figure 2 and, for each city x
∈
P, we obtain the system of four differential equations given by
(3)
D
t
γ
S
x
=
b
x
-
(
β
x
(
t
)
I
x
+
μ
x
)
S
x
+
∑
y
∈
P
‍
S
y
m
y
x
S
-
S
x
∑
y
∈
P
‍
m
x
y
S
,
D
t
γ
E
x
=
β
x
(
t
)
S
x
I
x
-
(
σ
x
+
μ
x
)
E
x
+
∑
y
∈
P
‍
E
y
m
y
x
E
-
E
x
∑
y
∈
P
‍
m
x
y
E
,
D
t
γ
I
x
=
σ
x
E
x
-
(
ζ
x
+
μ
x
)
I
x
+
∑
y
∈
P
‍
I
y
m
y
x
I
-
I
x
∑
y
∈
P
‍
m
x
y
I
,
D
t
γ
R
x
=
ζ
x
I
x
-
μ
x
R
x
+
∑
y
∈
P
‍
R
y
m
y
x
R
-
R
x
∑
y
∈
P
‍
m
x
y
R
,
with initial conditions
(4)
S
x
(
0
)
>
0
,
E
x
(
0
)
≥
0
,
I
x
(
0
)
≥
0
,
R
x
(
0
)
≥
0
,
with
∑
x
∈
P
‍
E
x
(
0
)
+
I
x
(
0
)
≥
0
.
These initial conditions ensure that, at the beginning, there are infected individuals in the system. They also imply the well posedness of the model. Indeed, we have the generalized version of the Mean Value Theorem that is stated as follows.Figure 2
Travel model dynamics between cities. The arrows symbolize the travel ratesm
x
y between cities x and y.Theorem 1.
Let the functionJ
∈
C
[
t
1
,
t
2
] and its fractional derivative D
t
α
J
∈
C
(
t
1
,
t
2
] for 0
≤
α
<
1, and t
1
,
t
2
∈
R; then one has
(5)
J
(
t
)
=
J
(
t
1
)
+
1
Γ
(
α
)
D
t
α
J
(
τ
)
(
t
-
t
1
)
α
∀
t
∈
(
t
1
,
t
2
]
,
where 0
≤
τ
<
t.Proof.
See the “Generalized Mean Value Theorem” proved in [11].Thus, consider the interval[
0
,
t
2
] for any t
2
>
o; this theorem implies that the function J
:
[
0
,
t
2
]
→
R
+ is nonincreasing on [
0
,
t
2
] if D
t
α
J
(
t
)
≤
0 for all t
∈
(
0
,
t
2
) and nondecreasing on [
0
,
t
2
] if D
t
α
J
(
t
)
≥
0 for all t
∈
(
0
,
t
2
). Furthermore, the nonnegativity of solutions to the problem (3)-(4) is obvious since the total population in city x, say N
x, evolves according to the sum of the four equations in (3) such that a positive N
x yields nonnegative solutions of (3)-(4).


## 4. Equilibrium Points, Basic Reproduction Ratio, and Stability

We evaluate the equilibrium points of (3) by letting, for each city x
∈
P,
(6)
D
t
γ
S
x
=
0
,
D
t
γ
E
x
=
0
,
D
t
γ
I
x
=
0
,
D
t
γ
R
x
=
0
.
The disease-free equilibrium (DFE) of the system is obtained by setting E
x
(
0
)
=
0
,
I
x
(
0
)
=
0 for each city x
∈
P and its stability is investigated using the next generation operator [12]. Thus, at the DFE for all the cities x
∈
P, we have S
x
=
N
x and adding the equation given in (3) yields at the DFE
(7)
b
x
-
μ
x
N
x
+
∑
y
∈
P
‍
N
y
m
y
x
S
-
∑
y
∈
P
‍
N
x
m
x
y
S
=
0
.
Following the same approach as [7], we compute the basic reproduction number as the spectral radius
(8)
R
0
=
ρ
(
F
V
-
1
)
,
where F and V are, respectively, the matrix of new infections and the matrix of transfer between compartments. The disease-free equilibrium of system (3) is linearly stable if R
0
<
1 and unstable if R
0
>
1.


## 5. Numerical Simulations

To provide numerical approximations and simulations of solutions to the fractional model (3)-(4) with different values of γ, we use the implementation code of the predictor-corrector PECE method of Adams-Bashforth-Moulton type described in [13]. We restrict our analysis to the case of symmetry in the intercities movements, meaning m
x
y
=
m
y
x. We assume that m
x
y is inversely proportional to the distance D
x
y between the cities x and y. That is,
(9)
m
x
y
D
x
y
=
k
,
for some k
∈
R. Table 1 shows the different distances between cities. Figure 3 shows the dynamics of the disease in the city of Auckland (A), for the four compartments with different values of the fractional derivative order γ
=
1
;
0.95 and 0.90. Auckland has been chosen just as an example since it contains the highest number of people among the four cities as shown in Table 1. In Figure 4, we can see the global evolution the infection in the four cities, with different values of γ. It is clear that the disease dies out irrespective of the three values of γ. In the simulation for the city x
=
A (Auckland), for example, we have taken E
A
(
0
)
=
1, I
A
(
0
)
=
0.01, and R
A
(
0
)
=
0.02 and have considered the parametersζ
A
=
0.14,
μ
A
=
0,
β
A
=
0.1,
k
=
1,
σ
A
=
0.09,
b
A
=
0.Table 1
Auckland
Bay of Plenty
Wellington
Dunedin
Initial population (×105)
S
A
(
0
) = 1417,9
S
B
(
0
) = 124,1
S
W
(
0
) = 397,5
S
D
(
0
) = 118,25
Auckland
0
Bay of Plenty
200
0
Wellington
650
470
0
Dunedin
1060
1020
610
0
Approximate distances between cities (km).Measles dynamics in Auckland (A) for the 4 compartments for γ
=
1; 0.95; 0.90.
(a)
(b)
(c)
(d)Global dynamics of the infection in the 4 cities forγ
=
1; 0.95; 0.90.
(a)
(b)
(c)
(d)


## 6. Concluding Remarks

In this paper we have presented the well posedness of fractional SEIR metapopulation dynamics of measles between four cities in New Zealand. We used real data according to the article by the online magazineOtago Daily Time, following the emergence of the disease after a student brought it from Sydney where he participated in a hip hop competition. We have started by showing nonnegativity of solutions to the fractional metapopulation model, thereby addressing the problem of its well posedness. We have also shown that the disease-free equilibrium of the model is linearly stable if R
0
<
1 and unstable if R
0
>
1. Numerical simulations have shown that, even in fractional dynamics of measles in metapopulation, the epidemic will not occur in communities which generate insufficient new hosts, which is in accordance with the theory of Bartlett [9]. This work generalizes the preceding ones with the inclusion of the fractional dynamics to a combined SEIR and metapopulation model, giving at the same time one of the multiple applications of fractional differential equations.



---

*Source: 781028-2014-06-22.xml*