Computational and Mathematical Methods in Medicine

Computational and Mathematical Methods in Medicine / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 6862516 | https://doi.org/10.1155/2020/6862516

Chinwendu E. Madubueze, Sambo Dachollom, Isaac Obiajulu Onwubuya, "Controlling the Spread of COVID-19: Optimal Control Analysis", Computational and Mathematical Methods in Medicine, vol. 2020, Article ID 6862516, 14 pages, 2020. https://doi.org/10.1155/2020/6862516

Controlling the Spread of COVID-19: Optimal Control Analysis

Academic Editor: Zoran Bursac
Received06 May 2020
Revised25 Aug 2020
Accepted05 Sep 2020
Published17 Sep 2020

Abstract

Coronavirus disease 2019 (COVID-19) is a disease caused by severe acute respiratory syndrome coronavirus 2 (SARS CoV-2). It was declared on March 11, 2020, by the World Health Organization as pandemic disease. The disease has neither approved medicine nor vaccine and has made governments and scholars search for drastic measures in combating the pandemic. Regrettably, the spread of the virus and mortality due to COVID-19 has continued to increase daily. Hence, it is imperative to control the spread of the disease particularly using nonpharmacological strategies such as quarantine, isolation, and public health education. This work studied the effect of these different control strategies as time-dependent interventions using mathematical modeling and optimal control approach to ascertain their contributions in the dynamic transmission of COVID-19. The model was proven to have an invariant region and was well-posed. The basic reproduction number and effective reproduction numbers were computed with and without interventions, respectively, and were used to carry out the sensitivity analysis that identified the critical parameters contributing to the spread of COVID-19. The optimal control analysis was carried out using the Pontryagin’s maximum principle to figure out the optimal strategy necessary to curtail the disease. The findings of the optimal control analysis and numerical simulations revealed that time-dependent interventions reduced the number of exposed and infected individuals compared to time-independent interventions. These interventions were time-bound and best implemented within the first 100 days of the outbreak. Again, the combined implementation of only two of these interventions produced a good result in reducing infection in the population. While, the combined implementation of all three interventions performed better, even though zero infection was not achieved in the population. This implied that multiple interventions need to be deployed early in order to reduce the virus to the barest minimum.

1. Introduction

The Novel Coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a new strain of coronaviruses that cause the coronavirus disease 2019 (COVID-19) and was declared a pandemic by the World Health Organization (WHO) on March 11, 2020 [1]. The virus was discovered in December 2019 in Wuhan City of Hubei Province, China [2, 3]. SARS-CoV-2 belongs to the order of Nidovirales, a family of Coronaviridae, and subfamily of Orthocoronavirinae [4]. Coronaviruses are a group of enveloped viruses with a positive-sense, single-stranded RNA and viral particles resembling a crown from which the name was derived [3].

The COVID-19 is a highly infectious disease that can be spread directly or indirectly from an infectious person to a healthy person through the eye, nose, and mouth via droplets produced when coughing or sneezing [2, 5, 6]. The exact source of the disease is uncertain. However, rodents and bats have been suspected by many researchers [5, 7]. The SARS-CoV-2 can survive up to 8-10 hours over porous surfaces (like paper, wood, sponge, and fabric) and a little more than 8-10 hours on nonpermeable surfaces (glass, plastics, metals, etc.) [2]. It has an incubation period of usually 2-14 days [2, 8]. Its symptoms are similar to that of common cold or flu. Also, others include fever, dry cough, shortness of breath, and pneumonia [5]. The severity of the illness can vary in different people from mild to severe symptoms based on their age and health status [2, 5, 6]. Almost 80% of COVID-19 patients are either asymptomatic or have mild symptoms and usually recover from the disease within 2 weeks. However, high mortality is recorded among the aged people and people with underlying chronic diseases, 2% of COVID-19 sufferers are under 18 years of age, out of which, fewer than 3% developed severe conditions [2]. COVID-19 has a low mortality rate that ranges from 2%-3%, which is significantly less than 10% of the severe acute respiratory syndrome (SARS) in 2003 and 35% of Middle-East respiratory syndrome (MERS) in 2012 [2, 3, 5, 9]. Due to its high infectivity, COVID-19 spread exponentially to virtually every part of the world within three months [1, 10].

As of April 2020, almost every country of the world has recorded at least one positive case despite speculation that the virus does not thrive in regions with hot weather. Entire Europe (especially, Italy and Spain) has become the epicentre of the outbreak, while the United States of America, Asia, and Australia record hundreds of new infections daily, with thousands of disease mortalities recorded [10]. In Africa, Nigeria recorded her first case on February 27, 2020, but as of April 26, 2020, just within 60 days, the figure has risen to 1364 infections with 45 deaths so far [11]. It is sad to note that most of the mortalities of COVID-19 especially in developing countries are attributed to poor medical facilities and medical personnel.

There is no specific treatment or vaccine available for COVID-19 probably because it is a new disease, and vaccine development usually takes up to 18 months [9, 12]. There is no approved medicine that eradicates the virus; however, treatment is mainly supportive [2, 9, 13]. It is because of these realities that governments across the world have resorted to nonpharmacologic measures. For instance, the Nigerian government has sensitized its citizenry on the need to adopt safety measures such as wearing of disposable surgical face masks, regular hand-washing with plenty of soap under running water, and the use of alcohol-based hand sanitizer in the absence of soap and water among others as recommended by the WHO [1, 10, 11]. Also, many governments worldwide are spending billions of the United States’ dollars as well as soliciting aids from well-spirited individuals and organizations towards combating the COVID-19 pandemic. Furthermore, many countries have imposed compulsory self-quarantine and restricted movements of their citizenries (lockdown/sit at home), closure of businesses, and borders as preventive measures [1, 10]. These interventions have succeeded greatly in curtailing the transborder spread of the SARS-CoV-2 from country to country. Nevertheless, the emerging major problem in the spread of COVID-19 is human-to-human transmission in a heterogeneous community. Sadly, the implementation of these interventional policies of governments (e.g., total lockdown of movement, businesses, and fear of quarantine/isolation) has thrown up another new challenge in the fight of the disease because of hunger and poverty especially in developing countries in sub-Saharan Africa where governments lack social securities. Therefore, there is the need to find cost-effective ways of halting the COVID-19 pandemic with minimal economic and social disruptions to avert impending catastrophic economic rupture.

Scholars are approaching this pursuit from two broad but complementary aspects of sciences: the medical and natural sciences. The medical scientists are busy trying to identify the source(s) of the disease, quicker ways of detecting the disease, treatment, and vaccine production [1, 12, 13]. The natural scientists are busy trying to proffer interventional measures through the development of mathematical models that will control the disease transmission especially now that there is no vaccine or known treatment.

Mathematical models have over the years proven to be reliable and efficient tools employed in formulating control strategies towards suppressing and mitigating the effects of infectious diseases, epidemics, and pandemics such as Ebola, SARS, and MERS [1416]. For COVID-19, some mathematical models have been produced which aim at halting the spread of the disease and forecasting its transmission through simulations. Some scholars focused on calculating the basic reproduction number, [8, 1719], and failed to consider the effect of public health education, quarantine, and isolation on the transmission of COVID-19. In the work of Imai et al. [9], they assumed COVID-19 is highly inconsistent in terms of the number of new infections just like SARS. They also investigated the consistency of their model with realities of the outbreak size using the set of simulated epidemic paths. Their findings affirm that without the implementation of holistic control measures, human-to-human transmissibility of COVID-19 is enough to sustain the pandemic and postulated that COVID-19 will have a diminutive generation time if the majority of COVID-19 cases have mild to moderate symptoms. Shen et al. [20] credited the high case detection rate and quick response by China and the world at large, to experiences from fighting the previous coronaviruses. Their findings postulate that COVID-19 may be a weak species in the coronavirus family, using the national epidemic of Wuhan in China with a fatality rate of 11.02% (9.26-12.78%) which is less than to those of SARS (14-15%) and MERS (34.4%) with a total of 8042 (95% CI: 4199-11884) infections and 898 (368-1429) death, respectively. Chen et al. [5] in their research simulated the potency of transmission of COVID-19 from bat (probable) source to humans using their Bats-Hosts-Reservoir-People transmission model [6]. They calculated the basic reproduction number using the next-generation matrix approach, and their results revealed that COVID-19 has higher transmissibility than MERs in the Middle East countries. Rabajante [21] reveals that more havoc and transmission/spread of COVID-19 is being perpetrated within the period an infected person is exposed. Rabajante [21] stated that such an infected person can transmit the virus, especially in a social/public gathering in a remote community within 14 days infectious period. Using the early models of COVID-19, Rabajante [21] recommends the maximum observance of control measures in any public/social events. Tang et al. [3] updated their previous model to a time-dependent model. They took into cognizance new interventional advances made in the COVID-19 fight. Their updated findings reported that the best control measure is persistent and constant strict self-isolation. They predicted that the pandemic will peak if the public health measures are adhered to. He et al. [6] studied the transmission of COVID-19 with binomial distributions in their discrete-time stochastic epidemic model. Their model parameters were derived from fitted reported data of China from January 11 to February 13, 2020. Their basic reproduction number affirms the positive contributions of various control measures recommended by WHO. While the result of numerical simulation suggests that the disease will peak around February 2020 with contact rate as a paramount factor in the control of COVID-19.

Bordered on how best to control the disease with minimal risk (optimization theory), assessing the consequences of some interventional measures and the risk involved especially now that antiviral treatment and vaccines are not yet available, it is crucial to investigate the optimal control of some control measures. Optimal control is the generalization of the classical calculus of variation in optimization theory. It involves minimizing the cost function and converting a given optimal control model into a Hamiltonian function and apply the Pontryagin’s maximum principle. Optimal control has been successfully applied to infectious diseases like HIV, Ebola, Tuberculosis, and SARS. For a disease like COVID-19 that spreads fast, the timing of implementing control measures is important. Unfortunately, very few researchers like [22, 23] considered an optimal control analysis of the COVID-19 transmission and suggest that more researches should be directed in this regard. Djidjou-Demasse et al. [22] in their work employed the concept of optimal control theory to explore the best control strategy to implement while awaiting the vaccine. They deduce that the only end to COVID-19 is when humans develop immunity. They weigh the options of humans developing natural immunity after infection or after they have been vaccinated. Their findings reveal that vaccination will best minimize the cost of loss of human lives, while maximal implementation of control strategies will peak the pandemic in four months after onset. Their results forecasted the possibility of having an efficient vaccine to be in 18 months. Moore and Okyere [23] attributed the rapid spread of COVID-19 to poor medical amenities. Their optimal control analysis focused on the controls: personal protection, treatment, and environmental spraying (environmental hygiene) as time-dependent control functions. Their numerical simulation reveals that optimal implementation of all the control measures greatly reduces exposed and infectious individuals in the population.

From the foregoing literature, interventions have been invested in, advocated for, and implemented by various stakeholders and still ongoing in the fight of COVID-19. These have cost a huge sum of money and time, casualties in businesses, economies, lives, etc. Unfortunately, the world is still recording high mortality and morbidity due to the disease. The few mathematical models that abound on COVID-19 suggesting diverse interventional control measures are yet to explore critically the optimal control analysis of those control parameters. This is necessary to ascertain their contributions in the dynamic transmission of COVID-19 for guidance in formulating better policies on the fight against COVID-19.

The model by Gumel et al. [16] forms the motivation for this study. Gumel et al. studied the impact of quarantine and isolation on the transmission dynamics of SARS. They assumed that everyone quarantined progress to isolation. The control measures in their work were assumed to be time-independent control measures. However, there is a possibility that some people will not develop symptoms after the quarantine. So they return to susceptible class to avoid being infected in the isolation centre. Also, it is well known that behavioral change played a very important role in the spread of diseases. Public health education contributed to people’s behavioral changes towards infectious diseases such as Cholera [24] and Ebola virus disease [25]. It will help the health personnel to reach out to people and influence them to adopt new behavioral changes and practice personal hygiene. Thus, this study seeks to ascertain the effectiveness of public health education, quarantine and isolation in reducing the infection of COVID-19, and the time taken to achieve that. It will establish the optimal control strategies required and the proportion of exposed individuals that will be quarantine to curtail the disease. It will seek the effect of time-dependent control variables and control constants on the transmission dynamics of COVID-19. The effect of constant controls will be explored using sensitivity analysis which will be used to identify the most sensitive model parameter that will be targeted. Pontryagin’s maximum principle will be applied to the optimal control model.

The rest of the paper is organized as follows: Section 2 is the model formulation for the COVID-19 with control measures. The model analysis for the COVID-19 is discussed in Section 3 with a sensitivity analysis of the model parameters. The formulation of an optimal control of the COVID-19 model and its analysis is done in Section 4, while Section 5 is the numerical simulations and its discussion. Section 6 is the conclusion.

2. Model Formulation

In this section, the formulation of a deterministic model for COVID-19 is presented. The model by Gumel et al. [16] used for the control of the SARS outbreak is extended for the control of COVID-19 in this study. The total population, , at time, , is divided into subpopulations: Susceptible, ; Exposed, ; Quarantined,; Infectious not hospitalized, ; Hospitalized/Isolated Infectious, ; and Recovered, . We further extend their model by incorporating public health education and the possibility of persons in the Quarantined,, who test negative for COVID-19 to return to the Susceptible, . The quarantined compartment comprises persons who come from high-risk regions and contacts of those who tested positive for COVID-19. These persons are kept for the incubation period of the virus. They are tested for the virus within this period. Those who test negative returned to Susceptible, , while those who test positive are taken to the compartment of Hospitalized/Isolated Infectious, . Those who miss quarantine but test positive are in the Infectious not hospitalized, , compartment from where they either recover because of their strong immunity or enter the compartment of Hospitalized/Isolated Infectious.

The human population at any given time, , is given by

In the Susceptible compartment, , a proportion of humans are recruited into the population at a constant rate, , through immigration/birth of no risk population and through a proportion, , of quarantine individuals that return to susceptible compartment after 14 days without symptoms at the rate, . People exit the susceptible compartment either through infection induced by the disease with the force of infection, . Infection is acquired via direct contact with infectious human contaminants or droplets. The public health education/awareness campaign, reduces the force of infection, , and it is time-dependent. The force of infection is given aswhere and are the modification factors for the exposed, quarantined, and hospitalized/isolated individuals. The parameters, and , are associated with the hygiene consciousness of the quarantine and the hospitalized/isolated individuals, respectively. The exposed compartment, , gains population through infection induced by the disease at the rate of and from a proportion, , of the recruitment of people immigrating from a high-risk population of COVID-19 at the rate of . A proportion, , of the exposed individuals exits through quarantine at the rate , and the remaining proportion, , of the exposed individuals escapes the quarantine and progresses to the infected compartment at the rate, due to ignorance or fear of being quarantined. We assumed that some of the immigrants in were either infected and at presymptomatic stage or had no infection [26]. This implies that not all individuals in will develop symptoms. Thus, individuals in who did not develop symptoms after the incubation period will progress to , while those who developed symptoms will progress to . This implies that a proportion, , of the quarantined individuals exits back to susceptible class after 14 days of no symptoms and reexamination at the rate, , while a proportion, , of the quarantined individuals that test positive progresses to Hospitalized/Isolated compartment at the rate. Also, the infectious not hospitalized individuals are isolated at the rate, , and can recover due to a boost in immunity at the rate, , and progress to the recovered compartment or die of the virus at the rate, . For the Isolated/Hospitalized compartment, , they gain population from a proportion, , of quarantined humans that become infectious during the 14 days quarantine period and the infectious not hospitalized individuals that are isolated. People exit the hospitalized/isolated compartment through recovery at the rate, , or death-induced rate, . Furthermore, the compartment of recovered, , gains population from the infectious not hospitalized individuals that miss isolation but recover due to boost in immunity, and from the hospitalized/isolated individuals at the rates of and , respectively. The recovered individuals are assumed to develop permanent immunity to COVID-19, and all the compartments exit their compartments through natural death rate, . The description of the parameters used in the COVID-19 model is given in Table 1. The system diagram for the transmission of COVID-19 is shown in Figure 1.


ParametersDescription

Public health education rate at time,
Quarantine rate at time,
The isolation rate for infectious not hospitalized individuals at time,
Transmission rate
The recovery rate for infectious not hospitalized individuals
Recovery rate of hospitalized/isolated individuals
Proportion of exposed individuals quarantined
The death-induced rate for infectious not hospitalized individuals
The death-induced rate for hospitalized/isolated individuals
Progression rate from quarantine to susceptible compartment after the incubation period
Progression rate for exposed individuals that missed quarantine to infectious not hospitalized compartment
The isolation rate of those that developed symptoms during the quarantine period
The proportion of persons coming from a high-risk area of COVID-19
Immigration rate

From the schematic diagram in Figure 1, the model equations are derived as followswith as the initial conditions.

3. Model Analysis

For the sake of model analysis, the controls, , are considered as constants, that is, . Investigation of some properties of model analysis will be carried out in order to understand the impact of the constant control parameters on the transmission dynamics of the COVID-19.

3.1. Positivity and Boundedness of Solutions
3.1.1. Invariant Region

The solutions of the model are uniformly bounded in a positive invariant region,

The total population at any time,, is given by (1) and

Solving equation (5) using Groonwall’s inequality gives . This means as , . Hence, the nonnegative solution set of the model equations (3) enters the feasible region, , which is a positively invariant set.

3.1.2. Positivity of the Solutions

The following theorem proves that the solution of the model is nonnegative for .

Theorem 1. Let the initial solutions satisfy . The model has nonnegative solutions which are contained in the feasible region, .

Proof. From the first equation of (3),Solving (6) gives .
In the same way,This shows that the solution set is nonnegative for all , since exponential functions and initial solutions are nonnegative.

3.2. Existence of Disease-Free Equilibrium State

From equation (3), and are represented as follows:

The model equations (3) become

The disease-free equilibrium state, , is established when there are no infective immigrants into the population (i.e., ) and when there is no disease in the community (i.e., ).

The equilibrium state for the model equations (9) is at the state when , and these are solved simultaneously to give the disease-free equilibrium state. , as

3.3. Basic Reproduction Number

The basic reproduction number, , is a threshold quantity that predicts the spread of disease in the population. It is an average number an infective will infect people in a wholly susceptible population. If , the infection will die out. If , the infection will persist in the population. The approach of Next-generation method by Driessche and Watmough [27] is used to compute .

The rates of new infection and the transfer from in and out the infected compartments are given by

The partial derivatives of and at the DFE, , yield

The basic reproduction number, , which is the spectral radius of the matrix, , is given aswhere

If , we have a perfect quarantine with no infectious not hospitalized individuals, i.e., . The reproduction number with perfect quarantine, iswhere .

If , there is no exposed individual in quarantine, the basic reproduction number, , is given bywhere .

3.4. Existence of Endemic Equilibrium State

The equilibrium state for the COVID-19 model (9) is obtained by solving

Solving these simultaneously, these equations give

Substituting as and simplifying yields the following quadratic equationwhere

If , we have and , in equation (19) where .

So, equation (19) will become

This implies from equation (21) when ,

corresponds to disease-free equilibrium (DFE) state, , in equation (10), while represents the endemic equilibrium state, when .

Substituting into gives the endemic equilibrium state, where

When , we have from equation (19) that

Substituting as equation (25) into , we have

Equation (24) gives the endemic equilibrium state when while equations (25) and (28) give the endemic equilibrium state when provided that equations (25) and (28) satisfy the inequality,

3.5. Sensitivity Analysis of the Model Parameters

It is important to know the relative contribution of different model parameters responsible for the transmission and prevalence of any disease. This will help to identify where to focus interventions that will reduce the mortality and morbidity due to the disease. In this study, the sensitivity analysis is examined to identify crucial model parameters that will reduce the burden of the disease and also quantify the impact of each input parameter on the value of an outcome. The initial disease transmission is directly related to the basic reproduction number, . Therefore, we perform a sensitivity analysis on to identify the most critical parameters that will curtail the spread of COVID-19. We use forward normalized sensitivity index of to measure the relative change in , to the relative change in the model parameter . This is also defined using partial derivatives if is a differentiable function of the model parameter, , as is defined in Chitnis et al. [28]; Sanchez and Blowe [29] bywhere is the sensitivity index of with respect to parameter, .

We compute the sensitivity indices for each parameter in . For instance, the sensitivity index of for is given as

Table 2 shows the sensitivity indices of for other parameters in and their parameter values. The parameter values are taken from the literature on COVID-19, SARS, and MERS.


ParameterValueSourceIndex signSensitivity index value

[30]
0.3[31]
0.1Assumed
0.1[31]
[29]+
[29]
Assumed
0.07143 [29]
0.0079 [31]
0.03521 [21]
[21]
[31]
[29]
[29]

From Table 2, the sensitivity indices with negative signs indicate that the value of decreases when they are increasing, while the sensitivity indices with positive signs show that the value of increases when they are increasing. The sensitivity analysis shows that the most sensitive parameters are in the descending order of and so on. These parameters will halt the spread of COVID-19 by reducing and increasing . It implies that the control parameters, , will reduce the spread of COVID-19 if they are increased. This is also shown in Figure 2 for the impact of and on . This implies that increasing the rate of implementation of interventions such as awareness, quarantine, and isolation in the exposed and infected not hospitalized population will halt the spread of COVID-19. In reducing we may consider the behavioral change in the transmission rate for further research.

4. Optimal Control Analysis

The control effort, , represents the public health education effort in educating people about the importance of social distancing, stay at home, and hand-washing in halting the spread of COVID-19. The control effort, , represents the effort used to quarantine the exposed individuals, and the control effort, , represents the effort used to isolate the infected individuals. The public health education effort involves educating the public through social media, television, radio, and traditional rulers in the community on how to observe social distancing and hand washing. The efforts used to quarantine the exposed individuals and isolate infected individuals include recruiting and training of the health workers on how to wear personal protective equipment (PPE), tracing the contacts of those exposed to the COVID-19 through home visits and phone calls, counseling, provision of an ambulance to convey the infected individuals to the isolation centre, general/COVID-19 tests, provision of isolation centres for treatment, and other related logistics.

Our goal is to minimize the cost function given assubject to the system of differential equations (3). All control efforts, , are assumed to be bounded and Lebesgue measurable time-dependent functions on the interval , where is the final time. The control effort set is defined as

The parameters, , , and , are the balancing cost factors for the public health education effort, , the quarantine efforts, , and the isolation effort, , respectively. The terms, , , and , represent the costs associated with public health education, quarantine, and isolation, respectively. Based on the literature for the optimal control of epidemics, the cost of the controls is assumed to be nonlinear and quadratic [32, 33]. If , then 100% effort is applied in public health education, quarantine, and isolation, respectively, at time, . Conversely, if, then no public health education for the people, no quarantine is carried out for the exposed (latent) individuals, and no isolation for the infected not hospitalized individuals.

The control time-dependent parameters will be considered in this section. Our goal is to find an optimal control for public health education effort, , quarantine effort for exposed individuals, , and isolation effort for infected individuals, , such that

The necessary conditions that an optimal solution must satisfy are obtained by applying the Pontryagin’s Maximum Principle to the COVID-19 model of equation (3). This principle converts system (3) and equation (32) into a problem of minimizing pointwise Hamiltonian, , given as:where , denote the associated costate variables for the state variables .

Using equation (35), we state the following theorem.

Theorem 2. Given an optimal control and solutions of the corresponding state system (3) that minimizes over , there exist costate variables ,, that satisfy the following systems of equationswhere is defined in equation (2) and final time conditionsAlso, the optimality conditions, , and are given by

Proof. Differentiating the Hamiltonian function, , at the respective solutions of equations (3) and the optimal control with final time conditions, the differential equations governing the costate variables are obtained as follows:This gives the costate system in equations (36).
The optimality conditions are given in the interior of the control setasSolving for as , as , and as , yield

Thus, using the bounds of the controls, , , and , the optimal control efforts in the compact form are given by equation (38).

The equations (3), (36) with optimality conditions (38) and the initial conditions, , and final time conditions (37) gives the optimality system.

Owing to the priori boundedness of the state variables, the costate functions, and the resulting Lipschitz structure of the ODEs, the uniqueness of the solutions of the optimality system is obtained for the small-time interval, . Hence, the bounded solutions to the optimality system are unique for .

5. Numerical Simulations and Discussion

5.1. Numerical Simulations

We carried out numerical simulations to investigate the impact of public education, quarantine, isolation, and the proportion of exposed individuals that will be quarantined. This is implemented using the parameter values and initial conditions from the literature on COVID-19, SARS, and MERS [31, 34]. The initial conditions for the state variables are as follows: , , , , , and , while , . The balance costs in the objective function are , , and . The other parameter values are given in Table 2. The optimality system is solved using the forward-backward sweep scheme. The details of the scheme are presented by Lenhart and Workman [35]. Many researchers have computed different values of the basic reproduction number for the person to person transmission, reservoir to person transmission, environmental transmission, and some of their results have been compared with other types of coronaviruses, SARS, and MERS, and their results show almost similar results [5, 36]. Hence, we focus our numerical simulation on the impact of a different combination of control interventions with their different control profiles on the transmission dynamics of COVID-19.

6. Discussion

A six compartmental model for the transmission dynamics of COVID-19 with quarantine, isolation, and public health education as time-dependent control measures is examined using the work of Gumel et al. [14] as a guide. The model is for human-human transmission that involves imported cases and community spread. The model is proven to have an invariant region. This region is where the model is well-posed and makes biological sense to be carried out for the human population. The basic reproduction number, , is when none of the exposed individuals are quarantined compared to when all of the exposed individuals are quarantined. This means that a single infected person can transmit the infection to approximately two other persons when there is no quarantine, while there is a possibility of stopping further transmission of infection when quarantine is implemented. However, there is a chance that some of the exposed individuals evade quarantine due to fear of stigmatization and death. Therefore, public health education/awareness will help to correct their misconceptions and encourage them to accept the control measures. Furthermore, people that have recovered need to share their experiences in quarantine and isolation centres with members of their community to enlist the cooperation of the entire community. When there is no isolation of the infected not hospitalized individuals in the population, the basic reproduction number, . This means that one infected not hospitalized person will infect approximately three persons in the population. The presence of isolation will help to reduce the number of infected not hospitalized individuals in the population. The simultaneous implementation of the three interventions reduces the number of infected individuals compared to the implementation of two interventions in the infected population (Figure 3(a)). This implementation takes about 64% input of awareness for 80 days, 58% input of quarantine for 95 days, and 100% input of isolation for 98 days before they drop slowly to their lower bound (see Figure 4(a)). This does not achieve zero infection in the population which implies that more interventions are needed to eradicate the virus. On the other hand, the combined implementation of public health education and quarantine measures produces a better result for the exposed population (Figure 3(b)). It takes about 100% input of public health education for 90 days and 100% input of quarantine for 95 days to trace 2000 contacts in the exposed population (Figure 4(b)). It means that public health education/awareness should reach all the hooks and corners of the population. People need to be aware of the virus, and also, the creation of adequate awareness of COVID-19 among the population will facilitate contact tracing and quarantine of high-risk individuals. It will also help to identify those who do not qualify for quarantine but tested positive for COVID-19 to be isolated. This will reduce further transmission of COVID-19 in the population.

Furthermore, time is of importance in implementing these interventions (see Figure 5). The number of exposed and infected individuals in the time-dependent interventions is 3,602 and 695 compared to 10,690 and 2,531 in the time-independent interventions, respectively. This implies that 7,088 and 1,836 individuals will not be exposed and infected, respectively, if interventions are implemented timely. These interventions are good to implement early which is the first 2-10 days of the outbreak. This will keep the burden of COVID-19 low. The virus will remain in the population for a prolonged time if there are no adequate interventions in place, but it will eventually drop over time (see Figure 6). With the interventions such as public health education/awareness, quarantine, and isolation, the number of exposed and infected individuals will reduce drastically within a short time but not to zero, leaving a residue of infected individuals with the potential to cause a further outbreak. This implies that COVID-19 will not be eradicated even with timely implementation of interventions unless a vaccine is developed.

When the proportion of the exposed individuals that are quarantine is increasing, it reduces the number of exposed individuals and infected individuals in the population (see Figure 7). When no exposed person is traced and quarantined in the population, the virus will remain in the population even when public health education/awareness and isolation interventions are present. When 30% of the exposed individuals are traced and quarantined immediately, the number of exposed individuals and infected individuals reduces to 2500 and 290 persons, respectively. The number of exposed individuals and infected individuals is about 1200 and 90 persons when 50% of the exposed individuals are quarantined. Again, quarantining 80% of exposed individuals will result in about 700 exposed persons and 20 infected persons in the population. This does not eradicate the infection in the population. To achieve zero infection in the population, we postulate that, additional interventions such as mass testing, and vaccination need to be incorporated. These were not done in this work and would be the focus of further research. This is in line with the new directive by WHO for research.

7. Conclusion

In this paper, a new deterministic mathematical model of COVID-19 was formulated with quarantine, isolation, and public health education as interventions. The model was also used as a prototyped to extensively investigate the contributions of these control measures to ascertain their individual and combined contributions in curbing the transmission and spread of COVID-19. The model analysis includes the establishment of the Invariant region and positivity of the model, the existence of disease-free equilibrium, and computation of the basic reproduction number . It was found that the basic reproduction number, , is when none of the exposed individuals are quarantined compared to when all of the exposed individuals are quarantined, . This means that a single infected person can transmit the infection to approximately two other persons when there is no quarantine, while there is a possibility of stopping further transmission of infection when there is quarantine. It was also shown that when there is no isolation of the infected not hospitalized individuals in the population, the basic reproduction number, . This means that one infected not hospitalized person will infect approximately three persons in the population. The presence of isolation will help to reduce the number of infected not hospitalized individuals in the population. The simultaneous implementation of the three interventions reduces the number of infected individuals compared to the implementation of two interventions in the infected population. Furthermore, it is observed that the time-dependent interventions reduce the number of exposed and infected individuals by 7,088 and 1,836, respectively. With the interventions such as quarantine, isolation, and public health education, the number of exposed and infected individuals will reduce drastically within a short time but not to zero, leaving a residue of infected individuals with the potential to cause a further outbreak. This implies that COVID-19 may not be eradicated even with the timely implementation of these interventions. Therefore, further interventions are needed to stop the spread of COVID-19.

Data Availability

The data used in this article are included within.

Additional Points

Recommendation. Governments should ensure prompt implementation of quarantine, isolation and public health education in the COVID-19 fight in order to suppress its infectivity and mitigate the disease burden. Further Research. Incorporation of mass testing and/or vaccination in the current model to ascertain its potential to eradicate COVID-19 in the population.

Conflicts of Interest

The authors declared there are no conflicts of interest.

Acknowledgments

We want to appreciate Dr. Sunday Madubueze for proofreading the manuscript and the respective departments for their support. We also thank the reviewers for their comments that improve this paper.

References

  1. World Health Organization, “Coronavirus Disease 2019 (COVID-19), Situation Report -51, Data as reported by 11 March 2020,” Tech. Rep., https://www.who.int/emergencies/disease/novel-coronavirus-2019/situation-reports. View at: Google Scholar
  2. 1mg, “Coronavirus No-Panic Help guide,” https://smef.org.uk/wp-content/uploads/2020/03/Corona-Ebook.pdf.pdf.pdf. View at: Google Scholar
  3. B. Tang, N. L. Bragazzi, Q. Li, S. Tang, Y. Xiao, and J. Wu, “An updated estimation of the risk of transmission of the novel coronavirus (2019-nCov),” Infectious Disease Modelling, vol. 5, pp. 248–255, 2020. View at: Publisher Site | Google Scholar
  4. W. G. Carlos, C. S. dela Cruz, B. Cao, S. Pasnick, and S. Jamil, “Novel Wuhan (2019-nCoV) Coronavirus,” American Journal of Respiratory and Critical Care Medicine, vol. 201, no. 4, pp. P7–P8, 2020. View at: Publisher Site | Google Scholar
  5. T. Chen, J. Rui, Q. Wang, Z. Zhao, J. Cui, and L. Yin, “A Mathematical Model for Simulating the Phase-Based Transmissibility of a Novel Coronavirus,” Journal of Infectious Disease of Poverty, vol. 9, no. 1, 2020. View at: Publisher Site | Google Scholar
  6. S. He, S. Tang, and L. Rong, “A discrete stochastic model of the COVID-19 outbreak: forecast and control,” AIMS Journal of Mathematical Biosciences and Engineering, vol. 17, no. 4, pp. 2792–2804, 2020. View at: Publisher Site | Google Scholar
  7. P. Zhou, X. L. Yang, X. G. Wang et al., “A pneumonia outbreak associated with a new coronavirus of probable bat origin,” Nature, vol. 579, no. 7798, pp. 270–273, 2020. View at: Publisher Site | Google Scholar
  8. Q. Li, X. Guan, P. Wu et al., “Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia,” The New England Journal of Medicine, vol. 382, no. 13, pp. 1199–1207, 2020. View at: Publisher Site | Google Scholar
  9. N. Imai, A. Cori, I. Dorigatti et al., Report 3: Transmissibility of 2019-nCoV, Imperial College London, 2020. View at: Publisher Site
  10. World Health Organization, “Coronavirus Disease 2019 (COVID-19), Situation Report -25,” Tech. Rep., https://www.who.int/emergencies/disease/novel-coronavirus-2019/situation-reports. View at: Google Scholar
  11. Nigerian Centre for Disease Control (NCDC), “COVID-19 SITUATION REPORT: Situation Report 1 and Report 58,” Tech. Rep., https://ncdc.gov.ng/disease/sitreps/?cat=14&name=An. View at: Google Scholar
  12. P. Surico and A. Galeotti, “The Economics of a Pandemic: The Case of COVID-19: Lectures,” in Wheeler Institute for Business and J-IDEA, Imperial College London Development, European Research Council, London Business School, 2020. View at: Google Scholar
  13. H. Lu, “Drug Treatment Options for the 2019-New Coronavirus (2019-nCoV),” Bioscience Trends, vol. 14, no. 1, pp. 69–71, 2020. View at: Publisher Site | Google Scholar
  14. C. E. Madubueze, A. R. Kimbir, and T. Aboiyar, “Global Stability of Ebola Virus Disease Model with Contact Tracing and Quarantine,” Applications & Applied Mathematics, vol. 13, no. 1, pp. 382–403, 2018. View at: Google Scholar
  15. Y. Kim, S. Lee, C. Chu, S. Choe, S. Hong, and Y. Shin, “The characteristics of middle eastern respiratory syndrome coronavirus transmission dynamics in South Korea,” Osong Public Health and Research Perspectives, vol. 7, no. 1, pp. 49–55, 2016. View at: Publisher Site | Google Scholar
  16. A. B. Gumel, S. Ruan, T. Day et al., “Modelling Strategies for Controlling SARS Outbreaks,” Proceedings of the Royal Society of London. Series B: Biological Sciences, vol. 271, no. 1554, pp. 2223–2232, 2004. View at: Publisher Site | Google Scholar
  17. J. T. Wu, K. Leung, and G. M. Leung, “Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study,” The Lancet, vol. 395, no. 10225, pp. 689–697, 2020. View at: Publisher Site | Google Scholar
  18. S. Zhao, Q. Lin, J. Ran et al., “Preliminary Estimation of the Basic Reproduction Number of Novel Coronavirus (2019-nCoV) in China, from 2019 to 2020: A Data-Driven Analysis in the Early Phase of the Outbreak,” International Journal of Infectious Diseases, vol. 92, pp. 214–217, 2020. View at: Publisher Site | Google Scholar
  19. S. Zhao, S. S. Musa, Q. Lin et al., “Estimating the unreported number of novel coronavirus (2019-nCoV) cases in China in the first half of January 2020: a data-driven modelling analysis of the early outbreak,” Journal of Clinical Medicine, vol. 9, no. 2, 2020. View at: Publisher Site | Google Scholar
  20. M. Shen, Z. Peng, Y. Xiao, and L. Zhang, Modelling the Epidemic Trend of the 2019 Novel Coronavirus Outbreak in China, bioRxiv preprint first posted online, 2020.
  21. J. F. Rabajante, “Insights from Early Mathematical Models of 2019-nCoV Acute Respiratory Disease (COVID-19) Dynamics,” 2020, https://arxiv.org/abs/2002.05296. View at: Google Scholar
  22. R. Djidjou-Demasse, Y. Michalakis, M. Choisy, M. T. Sofonea, and S. Alizona, Optimal COVID-19 Epidemic Control Until Vaccine Deployment, 2020, Preprint.
  23. S. E. Moore and E. Okyere, “Controlling the Transmission Dynamics of COVID-19,” https://arxiv.org/abs/2004.00443. View at: Google Scholar
  24. S. Edward and N. Nyerere, “A Mathematical Model for the Dynamics of Cholera with Control Measures,” Applied and Computational Mathematics, vol. 4, no. 2, pp. 53–63, 2015. View at: Publisher Site | Google Scholar
  25. S. Edward, E. M. Lusekelo, D. M. Ndidi, and E. Simanjilo, “Mathematical Modelling of the Transmission Dynamics of Ebola Virus Disease with Control Strategies,” International Journal of Sciences: Basic and Applied Research, vol. 33, no. 1, pp. 112–130, 2017. View at: Google Scholar
  26. World Health Organization, “Coronavirus Disease 2019 (COVID-19), Situation Report -73, Data as reported by 11 March 2020,” Tech. Rep., https://www.who.int/emergencies/disease/novel-coronavirus-2019/situation-reports. View at: Google Scholar
  27. P. Van Den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, no. 1-2, pp. 29–48, 2002. View at: Publisher Site | Google Scholar
  28. N. Chitnis, J. M. Hyman, and J. M. Cushing, “Determining important Parameters in the spread of malaria through the sensitivity analysis of a mathematical model,” Bulletin of mathematical biology, vol. 70, no. 5, pp. 1272–1296, 2008. View at: Publisher Site | Google Scholar
  29. M. A. Sanchez and S. M. Blower, “Uncertainty and sensitivity analysis of the basic reproductive rate: tuberculosis as an example,” American Journal of Epidemiology, vol. 145, no. 12, pp. 1127–1137, 1997. View at: Publisher Site | Google Scholar
  30. J. Gjorgjieva, K. Smith, G. Chowell, F. Sánchez, J. Snyder, and C. Castillo-Chavez, “The role of vaccination in the control of SARS,” Mathematical Biosciences and Engineering, vol. 2, no. 4, pp. 753–769, 2005. View at: Publisher Site | Google Scholar
  31. X. Yan and Y. Zou, “Optimal and sub-optimal quarantine and isolation control in SARS epidemics,” Mathematical and Computer Modelling, vol. 47, no. 1-2, pp. 235–245, 2008. View at: Publisher Site | Google Scholar
  32. A. V. Kamyad, R. Akbari, A. A. Heydari, and A. Heydari, “Mathematical modeling of transmission dynamics and optimal control of vaccination and treatment for hepatitis B virus,” Computational and Mathematical Methods in Medicine, vol. 2014, Article ID 475451, 15 pages, 2014. View at: Publisher Site | Google Scholar
  33. S. F. Sadiq, M. A. Khan, S. Islam, G. Zaman, I. H. Jung, and S. A. Khan, “Optimal control of an epidemic model of leptospirosis with nonlinear saturated incidences,” Annual Research and Review in Biology, vol. 4, no. 3, pp. 560–576, 2014. View at: Publisher Site | Google Scholar
  34. N. al-Asuoad, L. Rong, S. Alaswad, and M. Shillor, “Mathematical Model and Simulations of MERS Outbreak: Predictions and Implications for Control Measures,” Biomath, vol. 5, no. 2, article 1612141, 2016. View at: Publisher Site | Google Scholar
  35. S. Lenhart and J. T. Workman, Optimal Control Applied to Biological Models, Chapman and Hall/CRC, 2007. View at: Publisher Site
  36. C. Yang and J. Wang, “A Mathematical Model for the Novel Coronavirus Epidemic in Wuhan, China,” Mathematical Biosciences and Engineering, vol. 17, no. 3, pp. 2708–2724, 2020. View at: Publisher Site | Google Scholar

Copyright © 2020 Chinwendu E. Madubueze et al. 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.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views156
Downloads50
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.