Mathematical Model of the Spread of the Coronavirus Disease 2019 (COVID-19) in Burkina Faso ()
1. Introduction
Since December 2019, the Coronavirus disease 2019 (COVID-19) has been raving all the continents. Outbreak in Wuhan, China, the disease passed through Europe to finally reach Africa south of the Sahara. The Coronavirus disease 2019 (COVID-19) is an infectious disease that erupted in China in December 2019 and quickly spread into other countries in Europe and America. In West Africa, the first cases were reported in Senegal, Nigeria and Burkina Faso. In Burkina Faso, the first two cases were reported on March 9, 2020 ( [1] [2] [3] ); they were a Burkinabè couple who returned from travel from Mulhouse (France) after a few days stay there.
Currently, many mathematical models of the COVID-19 have been developed, mainly for the Coronavirus epidemic in China (see [1] [3] - [8] ). We have developed a model describing the Coronavirus epidemic in Burkina Faso, focussing on the effects of public policies imposed by the government to contain this epidemic and the number of reported occurred cases.
2. Mathematical Model
Considering the known characteristics of the Coronavirus disease 2019 (COVID-19) pandemic, we assume that each person is in one of the following compartments:
• S (Susceptible) means the number of persons who are not infected by the disease pathogen at time t, so who are susceptible to infection.
• E (Exposed) represents the number of persons who are in the incubation period after being infected by the disease pathogen. These persons have no visible clinical sign of the disease. They can infect other people but with lower probability than people in the infectious comportments.
• I (Infectious) means the number of persons who start developing clinical signs, these persons are asymptomatic infectious cases.
• Ir (Reported symptomatic infectious cases) represents the number of persons who are reported and isolated at time t.
• Iu (Unreported symptomatic infectious cases) stands for the number of unreported symptomatic persons at time t.
• R (Recovered) represents a person who has survived the disease, is no longer infectious and has developed a natural immunity to the disease pathogen.
This leads to the following transfer diagram (Figure 1).
The mathematical model consists of the following system of ordinary differential equations:
(1)
![]()
Figure 1. Transfer diagram for the mathematical model of COVID-19.
where
•
is the contact rate of a person in state S (day−1) at time t.
•
is the transition rate of a person in state E (day−1).
•
is such that
is the average time during which asymptomatic infectious cases are symptomatic and
.
•
is the transition rate of a person in state
to the state R (day−1).
•
is the transition rate of a person in state
to the state R (day−1).
• N is the number of people in the territory before the start of the pandemic.
The initial data of the system is supplemented by
The time t is in day, the asymptomatic infectious individuals
are infectious for an average time period of
days. We also suppose that the population is constant i.e.
.
Proposition 2.1. The basic reproduction number for the model system (1) is defined by
(2)
Proof. We use the method of next generation matrix in [9] to compute the reproduction number
.
We get
(3)
We have
(4)
and
(5)
On
, we get
(6)
and
(7)
Thus, we obtain
(8)
The basic reproduction number is defined as the dominant eigeinvalue of the matrix
.
Therefore,
The basic reproduction number
is defined as the number of cases that one infected person generates on average during his infectious period, in an uninfected population and without any special control measures. This number does not change during the spread of the disease.
The effective reproduction number
is defined as the number of cases that one infected person generates during his infectious period. This effective reproduction number depends on time, so, on public policies (change during the spread of the disease). Furthermore,
and the spread of the disease slows when
.
Since
therefore
(9)
3. Data
In this part, we use real data generated by the Burkina Faso National Health Commission against Coronavirus Disease 2019 (COVID-19) called CORUS. Recall that the first cases were reported in Burkina Faso on March 9, 2020, From that date, we witnessed an exponential growth in the number of cumulated cases. Table 1 represents the daily cumulative number of reported infected cases and Table 2 the daily reported infected cases.
With these data, we can see the daily dispersion of the infected case (Figure 2), and the cumulative reported infected cases (Figure 3).
![]()
Figure 2. Daily infected reported number.
![]()
Figure 3. Cumulative infected reported number.
![]()
Table 1. Cumulative daily reported cases data from March 9, 2020 to March 27, 2020 for Burkina Faso by Burkina Faso National Health Commission against Coronavirus Disease 2019 (COVID-19).
![]()
Table 2. Daily reported cases data from March 9, 2020 to March 27, 2020 for Burkina Faso by Burkina Faso National Health Commission against Coronavirus Disease 2019 (COVID-19).
4. Model Parameters for COVID-19 in Burkina Faso
Some of the parameters used in the simulations presented in Section 5 are from the literature. However, due to specificity of our real data, we fit some other parameters mainly to adjust the cumulative curve
. This curve represents the cumulative number of reported symptomatic infectious cases at time t. So using the fact that
, we have:
For the beginning of the epidemic, we assume that the growth of the cumulative curve is exponential as in [10], i.e.
We use the cumulative number of reported symptomatic infectious cases from Table 1 to evaluate
and
.
So from [4] and [10], we have
With the data in Section 3, we fit the cumulative infectious cases reported number by :
so we have
and
.
5. Simulation and Comments
The data we use here in Table 3 were calibrated for the situation of China, the first big source of COVID-19 [1] [3] [4] [7]. Due to the spread of the disease,
![]()
Table 3. Summary of some parameters calibrated for COVID-19 from the literature and data fit in Burkina Faso.
new studies should be performed to analyze its behavior in other sanitary population and climatic conditions. Only a limited number of studies accepted by the scientific community are available and these are what we use in Table 3.
For the parameters used in this model, we present some important parameters and threshold values related to the Coronavirus Epidemic in Burkina Faso. In particular, we observe that the basic reproduction number
is bigger than other reproduction number values reported in the literature [3] [10] [12]. This could be explained by the fact that we have taken into account unreported infected persons.
We assume that the exponential increase phase of the epidemic is intrinsic to the population of each region. Also, the Susceptible population
is not significantly reduced over the time. We suppose that the entire population of Burkina Faso at the date
are susceptible so,
, the exposed population at the same date
. Applying the model to Ouagadougou, the capital city, we set
and
.
5.1. Without Any Public Policies until March 26, 2020
We use the cumulative curve to fit the parameter
. From the beginning of the epidemic until March 9, 2020, we assume that there was no public policy so we fit
. This could be corroborated by the fact that the cumulative infected data curve fits well with the component of the reported infectious. Figure 4 shows the evolution of the cumulative infected reported cases and the forecast
![]()
Figure 4. Fitting data with the model (1).
relating to the model. Figure 5 shows the evolutions of different types of infected people without any public policy in the country. For the entire Burkina Faso, we can reach the bar of six million infected people in the long term with five hundred thousand people in Ouagadougou as is shown in Figure 6.
5.2. With Pubic Policies Started on March 27, 2020
From this date, we decrease the rate of contact
due to different types of measures taken by the public authorities. We can see the effect of these measures on the spikes of different types of infectious cases which decrease according to the degrees of the measures taken. We can see the curve of the contact rate in Figure 7.
(10)
Depending on the public measures taken,
increases, so the contact rate
decreases and it is possible to limit the number of infected persons.
is chosen in such a way that the simulation for the time interval aligns with the cumulative reported case data. So we are able to predict the future values of the epidemic from the early cumulative reported data. The earlier the decisions public policies are made, the better the management of the epidemic. So for instance, Figures 8-13 show the evolution of different types of infectious cases in Burkina for
,
,
,
,
,
. Figure 14 and Figure 15 represent the evolution of the cumulative reported case when
and
.
![]()
Figure 5. Forcast of the evolution of different types of infectious in Burkina.
![]()
Figure 6. Forcast of the evolution of different types of infectious in Ouagadougou.
![]()
Figure 7. Contact rate
for
.
![]()
Figure 8. Forcast of the evolution of different types of infectious cases in Burkina for
.
![]()
Figure 9. Forcast of the evolution of different types of infectious cases in Burkina for
.
![]()
Figure 10. Forcast of the evolution of different types of infectious cases in Burkina for
.
![]()
Figure 11. Forcast of the evolution of different types of infectious cases in Burkina for
.
![]()
Figure 12. Forcast of the evolution of different types of infectious cases in Burkina for
.
![]()
Figure 13. Forcast of the evolution of different types of infectious cases in Burkina for
.
![]()
Figure 14. Evolution of the cumulative infectious cases in Burkina for
.
![]()
Figure 15. Evolution of the cumulative infectious cases in Burkina for
.
6. Conclusion
In this paper, we have developed a mathematical model of COVID-19 for Burkina Faso, inspired by models in [3] and [10]. We have been able to estimate some parameters which have made it possible to fit the model to real data from the start of the Epidemic up to March 27, 2020 (when public policies were introduced). It emerges from this model that the most important parameter here is the contact rate which is a time dependent function (with respect to the public policies taken). A drastic reduction of the contact rate can lead to a considerable reduction in the number of infectious and of the duration of the epidemic.